DISSERTATION Tropospheric delay parameters derived from GNSS-tracking data
of a fast-moving fleet of trains
Ausgeführt zum Zwecke der Erlangung des akademischen Grades eines Doktors der
technischen Wissenschaften unter der Leitung von
Ao.Univ.Prof. Dipl.-Ing. Dr.techn. Robert Weber
am E 120-4
Department für Geodäsie und Geoinformation
eingereicht an der Technischen Universität Wien
Fakultät für Mathematik und Geoinformation
von
Matthias Aichinger-Rosenberger Matr.Nr: 00926550 Ried 14/7 6306 Söll
Wien, Februar 2021 (Unterschrift Verfasser/in) (Unterschrift Betreuer/in)
i
“Wonach du sehnlich ausgeschaut, Es wurde dir beschieden. Du triumphierst und jubelst laut Jetzt hab ich endlich Frieden! Ach, Freundchen, rede nicht so wild, Bez¨ahme deine Zunge! Ein jeder Wunsch, wenn er erf¨ullt, Kriegtaugenblicklich Junge.” —Wilhelm Busch (1832-1908) ii Acknowledgements
Theverse on the previous page, composed by thefamousGermanwriterWilhelm Busch, Ihad stumbled upon awhile agoinsomeold bookathome. From there on,I knew thatIneeded to add it to this thesis as an opening line since there’spractically no better way of describing myacademiccareer over the years. And as I’mwriting this and therebyfinish this chapter,there’s already somethingnew and exciting coming up forthe future. Nevertheless before that is nowthe time to thank everybodywho made this thesis and my PhD studies possible. First and foremost, Iwouldliketothank Prof. Robert Weber, who served as my supervisorand mentor not only forthis thesis but at allaspectsofuniversitylife. Thanks forgiving me the opportunity to jointhe research divisionand pursue such an interesting PhDstudy. Furthermore thanks to: Johannes B¨ohm,again forthe chance to start an academic career at TU Wien and forcompetent advice with his expertknowledge in atmosphericmodelling and space geodesy. MarcusFranzGlaner,for thegreatcamaraderieand collaborationatworkand numerous precious advices and discussions on GNSS,PPP andeven notvery ”work-specific passions” likedistance running. GregorM¨o ller, foralot of expert advice on all thingsconcerning GNSS (and GNSS Meteorology in particular) as well as the greatcollaboration not limitedbyour relatively short time at theinstitutetogether. Natalia Hanna, foralot of expertadvice on dataassimilationand carrying out the assimilation testsfor thisthesis. Zoreh Adavi,for providing ray-traced troposphericdelays forthe validation of the results. In generaltoall other colleagues at the research divisionand the institutewho Ihad thepleasuretoknowand workwith. Although my time thereended already ayearago, Ireally enjoyedevery dayofitand will alwaysbegrateful forthisexperience.
iii iv
Finally Iwould liketothankmyfamily: Eva, foryourlove andcontinuous support, withoutwhich this thesis would not have been possible/finished in such ashort timeframe. Ilove you. Paula and Fiona, the greatest giftsofmylife. Iloveyou too. My parents, brothers and all other relatives who supportedmeand kept me motivated over all those yearsI’ve nowspent studying. Abstract
Watervapourdenotes one of the most importantparametersutilizedfor describing the stateand evolution of the atmosphere. Therefore, detailed knowledgeofits distribution is of immense importance forweather forecasting. Furthermore, it is also themost effective greenhouse gas and highlyvariable in both space andtime. Thus, it is obvious that high resolutionobservations arecrucialfor accurate precipitation forecasts, especially forshort-termpredictionofsevereweather. Althoughnot intentionally built forthis purpose, GlobalNavigation Satellite Systems(GNSS)have proventomeet those requirements. The derivation of meteorologicalparameters fromGNSSobservations is based on the fact that electromagnetic signals aredelayedwhen travelling through the atmosphere. Parametrisations of these delays, most prominentlythe Zenith Total Delay (ZTD) parameter, havebeen studied extensively and proventoprovide substantial benefits foratmospheric research and especially NumericalWeather Prediction(NWP) model performance.Typically,regional/globalnetworks of static reference stationsare utilizedtoderiveZTD alongwith other parametersofinterest in GNSS analysis (e.g. station coordinates).Results areused as acontributing data source fordetermining the initialconditions of NWP models, aprocess referred to as Data Assimilation (DA).
Thisthesisgoesbeyond the approach outlined above, showing howreasonable troposphericparameters canbederived from highly-kinematic, single-frequency (SF) GNSSdata. This data wasgathered on trains ofthe Austrian FederalRailways(OBB¨ ) and processed using the Precise Point Positioning (PPP) technique. The special natureofthe observations yieldsanumber of additionalchallenges, from appropriate pre-processing and extended outlier detection, to advanced strategies in the PPP setup and forusage in DA procedures. Furthermore, since only SF dataisprovided, the treatmentofthe ionosphere represents one of the major challenges of thisstudy. Therefore different strategies havebeeninvestigated and shown to provide satisfactory results, suitable forZTD estimation. Despite thesechallenging circumstances, reasonable results forZTD estimatescould be
v vi obtained forthe analysed test cases investigating different PPP processing strategies. For validation of the results, comparisons with ZTD calculated using data from ERA5, the latest reanalysis of the European Centre forMedium-RangeWeather Forecasts (ECMWF),werecarried out. They yield very highcorrelation and an overall agreement fromthe lowmillimetre-rangeupto5centimetre, depending on solution and analysed travelling track. First tests of assimilation into aNWP modelagain showthe reasonable qualityofthe resultsasbetween 30-100 %ofthe observations areaccepted by the model. Furthermore guidelines to an operationalprocessingand possibleextensions to advancedtroposphericparameters were outlined, in order to exploitthe huge benefits (horizontal/temporal resolution)ofthis specific dataset foroperational weather forecasting. Kurzfassung
Atmosph¨a rischem Wasserdampf kommtdurch seineRolle als wichtigstesTreibhausgas und seine gleichzeitig r¨aumlichund zeitlich hoch variable Verteilung eine besondere Bedeutung inder Meteorologie zu.Dadurch ist auch verst¨a ndlich dassein m¨oglichst dichtes Beobachtungsnetz f¨ur die Niederschlagsvorhersage von entscheidender Bedeu- tung ist.Wenn auch urspr¨unglich nicht zu diesem Zweck entwickelt,k¨o nnen Global Navigation SatelliteSystems (GNSS)hier Abhilfe schaffen. Die Bestimmung von atmosph¨a rischen Wasserdampfmengen mit GNSS basiert dabei auf dem Sachverhalt dass elektromagnetische Signaleauf ihremWeg durch die Atmosph¨a re durch die Pr¨a senz vonWasserdampf verz¨o gert werden. Parametrisierungen dieser Verz¨ogerung, wie der Zenith Total Delay(ZTD),k¨o nnen aus Netzwerken fixer Stationen in der GNSSAuswertung, zusammen mit z.B. Stationskoordinaten, bestimmt werden. Die gewonnenen Resultatewerden typischerweise alskomplement¨a rer Datensatzzur Bestimmung der Anfangsbedingungen in numerischen Wettermodellen verwendet, ein Prozessder Datenassimilation bezeichnet wird.
Die vorliegende Arbeit geht ¨u berdiese Ans¨a tze hinausindem sie die Nutzbarkeit von Einfrequenz-GNSS-Daten, welche in einem hoch kinematischen Umfeld beobachtet wurden, zur Bestimmung vonZTD zeigt. Dabei handelt es sich um Lokalisierungs- daten vonZ¨u gen der Oster¨ reichischen Bandesbahnen (OBB¨ ), welche mittels Precise PointPositioning (PPP) ausgewertet wurden. Die große Anzahl der bereitsmit GNSSEmpf¨a ngernausgestatteten Z¨uge verspricht dabei eine noch nieerreichte r¨aumliche (und auch zeitliche) Aufl¨osung vonZTD Daten ¨u berganz Oster¨ reich. Die speziellen Eigenschaften der Daten bringen allerdings eine Reihe vonHerausforderungen mitsich, voneiner deutlich aufwendigeren Ausreißerkorrektur bishin zu speziellen Prozessierungs- und Assimilierungsans¨a tzenwelche in dieser Arbeit getestet wurden. Durch die Verwendung vonEinfrequenz-Daten kommt der Korrektur von ionosph¨a rischen Laufzeitverz¨ogerungen eine besondershohe Bedeutung zu.Aus diesem Grund wurden hierf¨u rmehrere verschiedene Ans¨a tze getestet, welche alle zufriedenstellende Ergebnisse
vii viii liefern konnten. Trotz dieser herausfordernden Voraussetzungenkonnten durchwegs realistische Ergeb- nisse f¨ur dieanalysierten Testf¨alleerzielt werden. Vergleiche mit ausatmosph¨arsichen Reanalyse-Daten (ERA5)des European Centre forMedium-Range Weather Forecasts (ECMWF) gerechneten ZTDs zeigen eine sehr hohe Korrelation und mittlere Abwe- ichungen vomunteren mm-Bereichbis zu 5cm, abh¨a ngigvon verwendeterSoftwareund Teststrecke. Erste Tests von Assimilation in ein numerisches Wettermodell zeigendass, abh¨a ngig vonL¨o sung und Strecke, 30 -100%der Beobachtungen vomModell akzeptiert werden.Dar¨u berhinaus wurde auch noch eine Strategief¨u rdie Implementierung dieser Idee im operationellen Betrieb entwickelt bzw. vorgeschlagen, um die großen Vorteile dieser Datens¨a tze (sehr hohe r¨aumliche/zeitliche Aufl¨osung)f¨u rdie operationelle Wettervorhersage nutzbar zu machen. Contents
Acknowledgements iii
Abstract v
Kurzfassung vii
Contents viii
1Introduction 1 1.1Thesis contribution...... 3 1.1.1 Goals ...... 4 1.2 Thesis outline ...... 5
2GlobalNavigation Satellite Systems 7 2.1Principles of GNSS positioning ...... 8 2.2 GNSS types...... 8 2.2.1 GPS ...... 9 2.2.2 GLONASS ...... 12 2.2.3 Galileo ...... 13 2.3 Precise Point Positioning ...... 14 2.3.1 Observation equations ...... 15 2.3.2 Errorsources ...... 15 2.3.3 LinearCombinationsofGNSS Observations ...... 19 2.3.4 Integerambiguityresolution ...... 22 2.3.5 Kinematic PPP...... 22
3Physicsofthe atmosphere 23 3.1General overview ...... 23 3.1.1 Optical propertiesand energy budget ...... 23 3.1.2Atmosphericcomposition ...... 24
ix x CONTENTS
3.2Vertical structure andlayers of theatmosphere ...... 25 3.2.1 Hydrostatic balance ...... 25 3.2.2 Vertical structureofpressure/density...... 26 3.2.3 Vertical distributionoftemperature ...... 27 3.2.4 Atmospheric layers ...... 27 3.3 Water vapourinthe atmosphere ...... 29 3.3.1 Watervapourmeasures ...... 30 3.4 Precipitation ...... 32 3.4.1 Numerical modellingofprecipitation ...... 33
4Signalpropagation in the atmosphere 35 4.1 Maxwellequations ...... 35 4.2 Path delays ...... 36 4.2.1 Ionospheric delay...... 37 4.2.2 Troposphericdelay...... 39
5GNSS Meteorology 45 5.1History ...... 45 5.2 GNSS-based troposphereproducts forNWP ...... 46 5.2.1 Ground-based products ...... 47 5.2.2 Space-based products...... 49 5.3 Estimation oftroposphericparameters using PPP ...... 50 5.3.1 Kinematic platforms ...... 52
6GNSS dataprocessing 53 6.1Data sources ...... 53 6.1.1 Greenlight ...... 53 6.1.2 Railway track database...... 55 6.1.3 EPOSA ...... 55 6.2 Pre-Processing...... 56 6.2.1 DilutionofPrecision ...... 56 6.2.2 Multipathanalysis...... 57 6.2.3 Cycle slip detection ...... 58 6.2.4 Data cleaning ...... 59 6.3 Parameter estimation in PPPprocessing ...... 59 6.3.1 Least-squaresadjustment ...... 60 6.3.2 Kalman Filtering ...... 62 CONTENTS xi
6.4 SEID...... 63 6.4.1 Accuracy of SEID model ...... 65 6.4.2Kinematicimplementation ...... 66 6.5 Processingscheme: CSRS ...... 66 6.6 Processingscheme: PPP-Wizard ...... 67 6.7 Processingscheme: GAMP using SEID...... 68 6.7.1 GAMP ...... 69 6.8 Processingscheme: Height-constrained PPP ...... 69 6.9 Filteringmethods and performance metrics...... 71 6.9.1 Moving averagefiltering ...... 71 6.9.2 Error metrics forperformance evaluation...... 72
7ZTD validation and assimilation 75 7.1Validation using ERA5 reanalysis products ...... 76 7.1.1 ERA5...... 76 7.1.2 ZTD calculationfromERA fields...... 77 7.1.3Quality of ERA5-derived ZTD...... 78 7.2 Basicsofdata assimilation ...... 80 7.2.1 DA techniques ...... 81 7.3 Assimilation of GNSS data in NWP...... 83 7.4 Assimilation andverification methodology...... 84 7.4.1 Model descriptionand setup ...... 84
8Results 87 8.1Overviewofconducted casestudies ...... 87 8.1.1 Test cases ...... 88 8.1.2PPP processing settings...... 89 8.2 Case study 1: 11.05.2017 ...... 92 8.2.1 Rawdataquality...... 92 8.2.2 PPP noise analysis ...... 93 8.2.3ZTD validation...... 95 8.3 Case study 2: 28.09.2017 ...... 100 8.3.1 Rawdataquality...... 100 8.3.2 PPP noise analysis...... 101 8.3.3 SEIDprocessing...... 104 8.3.4 ZTD validation ...... 106 8.4 Case study 3: 29.11.2019...... 110 xii CONTENTS
8.4.1Raw dataquality...... 110 8.4.2 PPP noiseanalysis ...... 112 8.4.3 ZTD validation ...... 113 8.5 Case study 4: 01.12.2019 ...... 116 8.5.1 Rawdataquality...... 116 8.5.2 PPP noise analysis...... 118 8.5.3 ZTD validation ...... 122 8.6 ZTD sensitivityanalysis ...... 124 8.6.1 Impact of GLONASS observations...... 124 8.6.2 Elevation cut-off ...... 127 8.6.3 Estimationofhorizontal gradients ...... 129 8.6.4 Mis-modelling of ionosphere using SEID ...... 131 8.6.5Altitude-dependency of estimation errors...... 132 8.7 Assimilation tests in WRF...... 135 8.7.1 WRFsetup ...... 135 8.7.2 Influenceofhorizontal resolution...... 135 8.7.3 Assimilation test:CS3 ...... 138 8.7.4 Assimilation test:CS2 ...... 138
9Towards real-timeprocessing 143 9.1Requiredadaptations forNRT GNSS processing ...... 144 9.1.1 Data transmission ...... 144 9.1.2 Orbit/clock products ...... 144 9.1.3 Processingsoftware ...... 144 9.1.4 Outputinterval ...... 145 9.1.5 Data format of NRTproducts...... 146 9.2 Purposedworkflowfor NRT processing ...... 146 9.2.1 SF-PPP using regionalionosphere maps ...... 146 9.2.2 DF-PPP usingSEID model ...... 146 9.3 Requirements forNWP assimilation and applications ...... 148 9.3.1 DA configuration ...... 148 9.4 Extension to advanced troposphericparameters ...... 150
10 Conclusions and outlook 151
Acronyms 157
Bibliography 171 CONTENTS xiii
Curriculum Vitae 187 xiv Chapter 1
Introduction
In times of extreme weather eventsand climate change, monitoring of the Earth’sat- mosphere is becoming amorecrucialtaskthan ever before. One of the most important, yetalso at least understood ingredientsofthis complex system is water vapour. It is not only the keyparameter in formationofcloudsand thus precipitation, but also the most importantgreenhouse gasaffecting global warming. In relationtoits importance, observationsofatmosphericwatervapourare relatively sparse. Radiosondes and other measurement techniques aretypically used to observe itsconcentration and distribution, but these techniques often suffer fromlow temporal and horizontal resolution, e.g.pro- viding measurements only once aday. This is wherethe field of Global NavigationSatelliteSystems(GNSS)can provide an alternative. Theretrieval of watervapourfrom GNSS signals relies on the fact that electromagneticwavesare delayedontheir way through theatmosphere. In terms of the neutral atmosphere, this delayoriginatesdue to the presence of waterand thus can be interpreted as ameasureofwatervapourinthe columnabove the receiver. Tropospheric parameters, derived fromGNSSsensors,are awell-established datasource forNumericalWeather Prediction (NWP) nowadays. The widespreadusage of GNSS systems formultiple fields of applicationprovides ahuge amount of observations avail- able forthe meteorological community.Basedonthe concept of GNSS Meteorology, outlined by Bevisetal. (1992), tropospheric delays areestimated operationally both in national and international networks, likethe European EGVAP network (Guerovaetal. 2016). Assimilation of these observations in NWP models is nowadays awidely-used technique to improveprecipitation andcloud cover forecasts. Several new developments havebeen achieved in the last twodecades, e.g. the establishmentofGNSS station networksused solely foratmosphericmonitoring or the usabilityofthe Galileosystem. Tropospheric estimatesare typicallyderived within adouble-difference (DD) network
1 2 Introduction solution,which can providethe high accuracy needed in NWP (∼ 1cmfor ZTD).This approach provides the most accurate results but has the downside of being cost-intensive since only high-end infrastructurecan be used and alargernumber of referencestations is needed forastable solution. Aside of this ‘classical’ approach, severalauthors havealso studied the potential of low-cost GNSS receiversinterms of tropospheric monitoring.Therefore, Precise Point Positioning (PPP)denotes the most promisingapproach forprocessingnowadays, espe- cially when resultshave to be availablein(near)real-time.The processing does not need any data fromreference stations and thus computational speed is increased remarkably. Another advantage is the fact that troposphericestimates fromPPP arenot biasedby correlations as typicallyintroduced in small-scale local station networks when applying DD processing (M¨oller 2017). These can prevent the detection of local (site-specific) watervapourtrends/anomalies and therefore can be problematicfor investigation very small-scale weather phenomena.However,when using PPP,high qualityinformation of satellite orbits and clocks (provided by institutionslikethe International GNSS Service (IGS)) is mandatory to obtain resultswith high accuracy. Low-cost receivers often only provide single-frequency (SF)observations. Amajorprob- lemwhen working with SF data is the treatment of the non-neutral atmosphere, the so-called ionosphere. For dual-frequency (DF) observations its effect can be eliminated by building the appropriatelinearcombination usingobservations on two frequencies. However, this is notpossiblewith SF observations and thus theeffect must be modelled. If ionospheric delaysand other majorerror sourcescan be modelled with sufficient accu- racy,moststudies showpromisingresults in terms of positioning as well as tropospheric monitoring. Until today, most studies concentratedonmeasurements from staticreceivers due more trivial data processingand more controllable measurement conditions. Gao et al. (2006)developed areal-time SF-PPPsystemand tested it in both road and marineenvironments. Both test scenarios achieved decimetreaccuracy in kinematic pro- cessingmode. de Bakker and Tiberius (2017)demonstrated lane identification and mappingonamo- torway using SF-PPP.Performance wasfound to be strongly dependent on the number of available satellites. Deng et al. (2009)densified an already existing station networkwith SF receivers and retrieved tropospheric delays using staticPPP in post-processing mode. In order to deal withthe ionospheric influence, they developedthe Satellite-specific Epoch-differenced Ionospheric Delay(SEID) model, which uses information fromsurrounding reference sta- 1.1Thesis contribution 3 tions to generate syntheticobservationsonasecond frequency forthe SF receiver. This approachprovided promisingresults andwill be tested in the course of this study as well. Krietemeyer et al. (2018)studiedthe feasibility of using low-cost receivers to increasethe density of GNSS networks forretrievalofPrecipitable WaterVapour(PWV). They pro- cessed oneyearofGNSS data fromanIGS station and two co-located single-frequency stations also by using the SEID algorithm to correct forthe ionospheric influence. Mendez Astudilloetal. (2018)carriedout aperformance assessment of three PPP on- lineservices and three softwarepackagesconcerningZTD estimation.Inthis context, they compared ZTDestimates from PPP processing to those obtained from the IGS tropospheric products. The resultsindicate that the online PPP services performance is better than thoseofselected PPPsoftwarepackages in all cases. Zhaoetal. (2019)estimated both PWV from low-cost multi-GNSS receivers in static mode. They retrieved ZTDs with an accuracy betterthan 10 mm compared with tro- posphere productspublished by IGS and about 3mmfor PWV compared to radiosonde results. Only very fewstudies have investigated the possibilityofmonitoring tropospheric pa- rametersfrom kinematic platforms. Most of these kinematic studies used maritime environments(shipborne (Rocken et al. (2005), Bonifaceetal. (2012), Fujita et al. (2008), Kealyetal. (2012)) and buoys(Chadwell and Bock (2001), Dodson et al. (2001)) or airborne data (Skone et al. (2006)).Uptill now, theonly investigation on train-borne GNSS data wascarriedout by Webb et al. (2016), who analysedZenith Wet Delay(ZWD) forover adecade of continuous observationsoverarangeofaltitudes on theSnowdon Mountain Railway,locatedinSnowdonia NationalPark, NorthWales. Theyfound an improvement of amulti-system GNSS solution in kinematic PPP,with an agreement of 11.6 mm forthe multi-GNSSsolution compared to 16.2 mm forthe GPS-only solutioninZWD.
1.1 Thesiscontribution
Thisstudy goes one step beyond the studies presented before by making useofthe kinematic PPP approach to derive tropospheric parametersfrom SF observationstaken on fast-movingtrains. These trains areequippedwith dual-system SF receivers, which currently areable to track GPS and GLONASS signals. In afuture state the receivers will be switched to GPS/GALILEO capability. Data is gathered on different railway tracks over Austria in the course of theproject Greenlight,led by theAustrianFederalRailways (OBB¨ ). 4 Introduction
Ahigher number of challenges is faced due to the fast movement of the sensorand a large number of obstructed areas on the track (e.g. tunnels). The potentialbenefits areavery highspatial and temporal resolution forobservations gathered attrain tracks covering huge partsofAustria and neighbouring countries. Moreover, ahighnumber of trains (∼ 1400) is already equipped withGNSS sensors and up to 200vehicles more are planned in the course of the next year. In particularthis thesis aims to answer the following questions:
• Can reliable troposphericestimates be derivedfrom GNSStrackingdataoffast- moving trains?
• Which pre-processing steps arerequired to prepare the obtained tracking data for PPPprocessing (cycle slip detection, multipath, ionospheric corrections)?
• Which dataprocessing scheme(s) is/areadequate to obtain results comparableto independent ZTD sources (NWP,ray-tracing)?
• Whatisthe quality of the derived dataand is it sufficient to be incorporatedinto an NWPdata assimilation system?
• Can the estimation and assimilationbecarried out in a(near) real-time (NRT) operational mode?
First investigations, as presented in Aichinger-Rosenberger and Weber (2020), already indicatedthe possibility to deriveZTD estimates of reasonable quality, even forrailway tracks in highlychallenging environments. This thesis extends these investigations fora number of test tracks and processing schemes, analyses specific sensitivities in the GNSS processingand presents first testsofZTD assimilation into aNWP model.
1.1.1Goals
Directly linked to thosequestions raised aboveare the following goals of the thesis:
1. Goals related to GNSS processing:
• Data cleaning:Implement astrategyfor datacleaning, mostly concerning cycle slips and multipath
• GNSS data processing: Find the best possible processing scheme(s) forthe given GNSS observations,especially taking into account adequate modelling of ionospheric corrections 1.2Thesis outline 5
• Tropospheric estimation: Estimate tropospheric parameters(ZTD) using the developedprocessingscheme(s)
• Qualityassessment: Assessthe limits of reachable qualityfor tropospheric parameter estimates using these processing scheme(s) by comparison to NWP-derived and ray-tracedtropospheric delays
• Processing of test cases (train tracks): Process data foranumber of casestudies(availabletrain tracks) in post-processing mode withfocus on tropospheric estimation
• NRT GNSS data processing: Provide strategy and theoretical basis of possible extension frompost-processing to NRT operation
2. Goals related to NWP data assimilation:
• Observation qualitycheck: Assess if qualityofestimates is sufficient to enter the model
• Assimilation casestudies: Assimilation of GNSS troposphericestimates from moving trains intoaNWP model
1.2 Thesisoutline
Chapter 2 gives an introductiononthe most importantaspects of existing GNSS,the physical/geometrical principles behind them and their usage in different scientific and real-world applications. Chapter 3 summarizes information on theEarth’satmosphere and it’s governingpro- cesses whichare of importance forthisstudy. Chapter 4 introduces the theoretical background and physical laws behind the propaga- tionofelectromagnetic signals travelling through the atmosphere. Chapter 5 takesacloser look on thehistory and thelatest scientific achievements in the specific field of GNSS Meteorology. Chapter 6 describes the different processing schemes forGNSS data used in this study. Adetailed description of processing setups is given forall deployedapproaches. Further- more the utilized dataand necessary pre-processing steps areintroduced. Chapter 7 introduces the methods and datasets used forassessing the reached quality of thetropospheric estimates derivedinChapter 6.Resultsofthis validation procedure can serveasanindication forthe usabilityofthe results forassimilation. Furthermore it dealswith assimilationprocedure of the ZTD estimates, derived fromthe processing 6 Introduction schemesintroduced in Chapter 6,intoaNWP model. Therefore ashort introduction to the mathematicalbackground ofdataassimilation andthe utilized NWP model is given. Chapter 8 providesthe resultsfrom all test cases conducted. Resultsofdata processing (in termsofZTD andheightestimatesmainly) areshown as well as validation and as- similation results. Chapter 9 aims to serveasanoutlook to apossibleNRT operational mode forGNSS tracking data processing. Therefore the theoretical frameworkand ideas forimplemen- tation as well as challenges of this approach will be examined. Chapter 10 concludes the thesis with adiscussionand summary of the main findings. Furthermoreanoutlook on futureimplementations of this research idea is provided, whichincorporatesrecommendationsfrom Chapter 9. Chapter 2
Global Navigation Satellite Systems
Global Navigation Satellite Systems (GNSS) have revolutionised anumber of scientific re- search fields as well as practical applications over thelastdecades (Aichinger-Rosenberger 2018). Originally startingasacompletely military system, it wasopened forpublic usage in 1983. Since then,ithas been extensively used and studied in the scientific andthe public societyand awidespread field of applications has emerged over thelastyears. Thefollowing points denote just asmall selection of topics which rely on or areclosely related to GNSS:
• Navigation: reaching frompedestriantoautomotive/railway/aviation usage
• High-precisionpositioning
• High-precisiontiming applications
• Assisted/Autonomous driving
• Geomonitoring: detection of plate tectonics, earthquakes, landslides
• Atmospheric monitoring:weather forecasting, climateresearch
Thefollowing section introduces the principle of GNSS positioning, which providesthe foundation forall applications mentioned above.Basedonthis principle the observation equations areformulated and thetypical error sourcesare described. Furthermore the differenttypes of GNSS areoutlined. At last,the mainprocessingscheme used in this study is outlined, Precise Point Posi- tioning (PPP). There, adetailed descriptionofthe observables, governing equations and estimated parameters of the approach is given. Benefits and weaknesses of the technique as well as typical error sourcesinfluencing its performance areoutlined.
7 8 GlobalNavigation SatelliteSystems
2.1 PrinciplesofGNSSpositioning
The basic principle of GNSS positioning relies on very precise propagation time measure- ments foranelectromagnetic signal travelling fromasatellite to areceiver.Through exact knowledge of thispropagationtime Δt thegeometric distance r from satellite to receiver results from
r = c(Tu − Ts)=cΔt, (2.1) where c denotesthe speed of light, Tu the system time at which the signalleft the satellite and Ts the system time when it reaches the receiver. However, sincethe receiver and satellite clocks arenot synchronizedthe measurementsofpropagation time will be biased. Thus thecalculateddistance will not be thereal geometric distance r but the so-called pseudorange ρ.Itisgiven by
s ρ = r + c(δtu − δt ), (2.2)
s where δtu is the offset of the receiverclock fromsystem time and δt is theoffset of the satellite clock from system time. Since on boardofthe satellites atomic clocks areused, the error in the satelliteclock is marginal compared to the receiver clock. Moreover it typicallycancels through observation differencing or can be modelled using accurate correctionmodels (see Section 2.3 formoredetails).Therefore only δtu is left in Equation 2.2,
ρ = r + cδtu, (2.3) whichnow consists of four unknowns (X,Y and Zcoordinate of receiverposition and re- ceiver clock error). Thus, pseudorangeobservations from at least four satellites (ρ1,.., ρ4) enable thecalculationofthe receiver coordinates. Additionalobservations result in an over-determined problemwhich is typically solved by least-squares adjustment. Geometrically this concept corresponds to the intersectionofthe surfaces of spheres, which havethe measured pseudoranges as radius. This is indicated in Figure 2.1,sim- plified to the 2D case (using circles instead of spheres). Typically only three satellites would be sufficient to calculate the 3D receiver position, as threespheres would intersect in one point on Earth. Nevertheless, introducingthe receiver clock error as an unknown likementioned before gives aminimumnumber of four satellitesneeded.
2.2 GNSS types
Since the start of their public availability in 1983, multipleGNSStypes havebeen built up. These include the AmericanGPS, RussiansystemGLONASS, European system 2.2GNSS types 9
Figure 2.1: Geometric conceptofGNSSpositioning, shown here inthe 2D case. Takenfrom B¨ohringer (2008)
Galileo andthe Chinese system Beidou. The first three mentioned systems aredescribed in more detail in the nextsection.
2.2.1 GPS
Since the early 1960s, several U.S governmentorganizations,including theDepartment of Defence (DOD),were interested in developing asatellite systemfor three-dimensional position determination(Kaplan and Hegarty2006)asthe successorofthe already oper- ating TransitSatelliteNavigationSystem. Theaim wastosecure the following attributes forthe system:
• global coverage
• continuous/all-weather operation
• abilitytoservehigh-dynamic platforms
• high accuracy 10 GlobalNavigation SatelliteSystems
Component Frequency
Fundamental frequency f0 =10.23MHz
Carrier L1 154·f0 =1575.42 MHz
Carrier L2 120·f0 =1227.60 MHz
Carrier L5 115·f0 =1176.45 MHz
Table 2.1: Components andfrequencies forcivil use of GPSsignal
In 1969, the Defense Navigation SatelliteSystem (DNSS) program wasestablished, an efforttoconcentratethe contributions of different organisations. As aresult, the Navstar Global Positioning System (GPS) (DoD 2020), aspace-borne positioning system for military purposes, wasfinally introduced in 1973. From thenonittookanother 10 years until the system became open to civil users in 1983.Since then GPS has been running as adual-use system, forboth military andcivil purposes.
GPS signals
GPS satellites use the so-called code divisionmultiple access (CDMA) technique to broadcastrangingcodes and navigation messages. Every satelliteemits an individual Pseudo Random Noise (PRN) code using this technique. These codes were initially solely modulatedontwo frequencies L1 (1575.42 MHz)and L2 (1227.6 MHz), which aremultipliersofthe fundamental frequency 10.23 MHz.Onto L1, Course/Acquisition (C/A) code and Precise(P1)codeismodulated as well as Precise (P2) code onto L2. In the course of the GPS modernization program anew third civil frequency L5 wasintroduced. This opened multiplenew possibilities of carrier phase combinations, especially useful forthe resolution of ambiguities and elimination of the ionosphericdelay. Table 2.1 shows an overviewofcarrierfrequencies used for(civil) GPS. In contrast to L5,the L3 =1381.05 MHz andL4=1379.913 carrier frequencies are only used formilitary purposes(Hofmann-Wellenhof et al. 2012). Since first tests of GPS revealed much betteraccuracy than originally expected, techniques likeSelective Availability(SA)were used to ensure superiorityofthe military over the civil user.This technique incorporatede.g. manipulation of orbital and satellite clock information for civil users to worsen positioning resultsuntil its use wasfinally terminated in May2000. However, encrytion of varioussignals formilitarypurposes(Anti-Spoofing)isstillinuse.
System components
In literature,GPS is usually described as asystems comprised of three different segments: 2.2GNSS types 11
• Space segment
• Control segment
• User segment
The space segment represents the satellite constellation, which consists nominallyof 24 operationalsatellites placedinsix nearly circularorbitsinclined by 55◦ with an altitude of about20200kmabove the Earth’ssurface (Karabatic 2011). This number of 24 satelliteswas reachedinJuly 1995and therebythe space segment became fully operational. Different typesofsatellites havebeen developedinthe course of the last 20-30years,which aresplit in differentblocks. By March 2019, 31 operationalsatellites areavailable anddelivering datadowntoEarth.
The control segment is build upon amaster control station, located in Colorado Springs, USA. This master station is responsible forcontrolling the satellites and keeping the system operating. This includes maintaining the satellites in their proper orbital positions (called stationkeeping)and monitoring satellite subsystem healthand status (Kaplan and Hegarty2006). In addition, five monitoring stations and anumber of ground antennas areemployed. Themonitoring stations makecontinuous measurements to all visiblesatellites and send this data to the master station. There, orbit and clock parametersare calculated and forwarded to theground antennas, which upload thepa- rameters three timesaday.Using this data, the navigation message forGPS is generated.
The user segment incorporateseveryuser of the system and is split into two dif- ferent categories, civil and military users. As already mentioned in the beginning of this chapter,the number of applicationsfor civil use haveincreased heavily over the last yearsand continuetodo so.They reach fromclassical tasks likepositioning, navigation or monitoring to meteorological and climatological applications. In principle, one distinguishes between single- and multi-frequency receivers depending on howmany frequency bandsobservations canbetaken.
GPSreference system
GPS coordinates aredefined in thethe DOD’s World Geodetic System1984 (WGS-84) (National Imagery and Mapping Agency 2000), whichrepresents the standard physical model of the Earth used forGPS applications(Kaplan and Hegarty2006). Alist of the defining parametersofthe underlyingreferenceellipsoid is givenisTable 2.2. 12 GlobalNavigation SatelliteSystems
Parameter andValue Description a =6378137.0 m Semi-majoraxisofthe ellipsoid f =1/298.257223563 Flattening of the ellipsoid −5 −1 ωE =7.292115·10 rad s Angular velocityofthe Earth µ =3986004.418·108 m3 s−2 Earth’s gravitational constant
Table 2.2: Defining parameters of WGS-84, taken from Hofmann-Wellenhof et al. (2012)
2.2.2GLONASS
Globalnaya navigatsionnaya sputnikovaya sistema (GLONASS)isthe Russian counterpart to theAmerican GPS system. LikeGPS originally designed as afully military service, its development started already in 1976.In1982 thefirst satellites were launched and in 1995 the full constellation wascompleted. Thelatterwas alsothe year in which the system became publicly available.Atthe end of the nineties, the capacityofthe system declined more and more due to lack of financial supply.This continued until 2001, when the restoration of thesystembecame one the top priorities of the Russian government. GLONASSismanaged by theRussianSpace Forces and operatedbytheCoordination Scientific Information Center of the Ministry of Defense of the RussianFederation (Xu 2007).
GLONASS signals
Just likeGPS, GLONASSsatellites transmittwo typesofsignals,astandard-precision (C/A code) and ahigh-precision signal (P code). In contrasttothe CDMA tech- niqueused by other GNSS systems, GLONASS uses frequency divisionmultiple access (FDMA), which means thateach satellitecan be identified by its ownuniquecarrier frequencies. These carrier frequencies (in MHz) aredefined by
G1=1602.0+0.5625k G2=1246.0+0.4375k where k denotesthe frequency channel. In the course of modernization, the CDMA technique is being investigated since 2008. Thisisexpected to provide multiple benefits forinteroperability with othersystems.Inaddition,anew signal on the third carrier as well as asecond civil signalonG2should be established. Since November 2018 24 GLONASS satellites areoperational andthusfulloperationalcapability(FOC) is reached. 2.2GNSS types 13
GLONASSreference system
Theterrestrialreference frame used by GLONASS is the ”PZ-90” (Earth Parameters 1990 –Parametry Zemli1990), whichuses the meanposition from1990to1995for the precise locationofthe North Pole. In contrast to this, WGS-84 as GPS reference system uses the location of the North Pole in 1984. The new version PZ-90.11 is broadcasted since 31 December 2013. It is aligned to theInternational TerrestrialReference System at epoch2011.0 at the centimeter level.
2.2.3 Galileo
Galileo is the European alternative to GPS and GLONASS, established by the European Union (EU) through acooperation of theEuropean Commission (EC), the European Space Agency(ESA) and at alaterstage the European GNSS Agency (GSA). Up to date (December 2020), thesystemconsists of 24 operationalsatellites.Infull operational mode observations from30satellitesshould become available to users. Contrary to itscounterparts GPS and GLONASS,Galileoisintended primarily forcivilian use, as basic (lower-precision) services will be free and open to everyone. From 2022 on,also high precision positioning will become available forusers with multiple signal receivers,the Galileo High-Accuracy Service(HAS).
Galileo signals
The Galileo system transmitsthree different signals: • E1
• E5 consistingofE5a andE5b
• E6
SimilartoGPS, these carriers arebased on thefundamentalfrequency f0 =10.23 MHz. In order to enable services like open,Search-And-Rescue (SAR), commercial or public regulated services,eachsatellitewill broadcastdifferent navigation messages (Karabatic 2011). An overview of all available carrier frequencies can be foundinTable 2.3. Whentakingacloser look on the chosen frequencies it can be noticed thatE1and E5a coincide withGPS signalsL1and L5, as well as E5b with the GLONASSG3carrier. This is not pure coincidence but chosen forthe sakeofbetterinteroperabilitywith the competing systems. Despite this benefit, it is also accompanied by thepossibleproblem of interference between the systems, which is avoided by using different modulation schemes forGalileo. 14 GlobalNavigation SatelliteSystems
Component Frequency
Fundamentalfrequency f0 =10.23MHz
Carrier E1 154·f0 =1575.42 MHz
Carrier E6 125·f0 =1278.75 MHz
Carrier E5 116.5·f0 =1191.795 MHz
Carrier E5a 115·f0 = 1176.45 MHz
Carrier E5b 118·f0 = 1207.14 MHz
Table 2.3: Components and frequencies of Galileo signal
2.3 Precise PointPositioning
Precise Point Positioning (PPP) is aGNSS-based positioning technique which uses un- differenced single- or dual-frequency carrier phase and/orpseudorangemeasurements from asingle receiver (Hinterberger 2016). Itsmainadvantage is independence of addi- tionaldata from (nearby)reference stations, unlikeotherprominent techniques such as Real Time Kinematic (RTK). Instead, precise orbit and clock information as provided by e.g. theInternationalGNSSservice (IGS) is utilized to ensure high-quality results. Using dual-frequency observationsincombination with finalorbit/clocksolutions, centimetre- level accuracy can be reached.This potential wasfirstshown by Zumberge et al. (1997) and since then PPPhas become one of the most promising topics in GNSSresearch. Huge scientific effortshave been put into certain aspects like ambiguityresolution(AR) or thecorrectmodelling of satellite and/orreceiver hardwarebiases. Especially AR has been the topic of anumber of studies, since in classicalPPP usuallyareal-valued con- stant is estimated as phase bias instead of the integer ambiguity. Overcoming this issue denotes one of the most promisingmeans to improvethe accuracy of results even further. Nevertheless, the technique stilldeals with anumber of restrictions which currentlylimit its potential andusage. The most prominent oneisthe often longer convergence time of the solutioncompared to RTK. In the following acomprehensivelist of advantages and disadvantages of PPPisgiven:
• Only asparse globalstation networkneeded (for computation of orbit/clock in- formation),which also implies cost-effectiveness (+)
• No simultaneousobservations (from nearbystation)required (+)
• No errorsintroduced from observationsatreference stations (networkeffects) (+)
• Error sources arerepresented as state-space correction models and notasrange 2.3Precise PointPositioning 15
corrections in the observation space (+)
• Globally available/valid corrections (+)
• Reduced computationaleffortsdue to smaller data amount which also implies faster availabilityofthe results (+)
• Longconvergence time (-)
• Complicated observationmodeling (-)
• Challengeofinteger ambiguityresolutionfor highest accuracy (-)
• High-rate and high-qualitysatellite clock information required (-)
2.3.1 Observation equations
The observation equations forthe PPP algorithmare,writtenseparately forpseudorange s observations Pr,ν
s s s Pr,ν = ρr − cδt + cδtr + δion,ν + δtro + P,ν (2.4)
s and phase observations Lr,ν
s s s s Lr,ν = ρr − cδt + cδtr − δion,ν + δtro + λνBr,ν + L,ν. (2.5)
Here thesubscript r denotes aGNSS receiver, s asatelliteand ν the corresponding s frequency.Furthermore ρr denotes the geometric distance between satelliteand receiver, s c the speed of light, δt and δtr theclock errorsofsatellite andreceiver, δion,ν the
(slant) ionosphericdelay on thecorresponding frequency, δtro the(slant) tropospheric s delay, λν Br,ν thefloat ambiguity forphase observations and P,ν/ L,ν areadditional error terms (e.g. multipath)for thecorresponding observation type.
2.3.2 Errorsources
Thissectionaims to give an overviewonsome of thedifferent error terms containedas in Equation 2.4 and 2.5.Avisual overview is given in Figure 2.2,which shows the most prominent delays/influences encountered forGNSS measurements. 16 GlobalNavigation SatelliteSystems
Figure 2.2: Component anderrorsources in pseudorange measurements, taken from Sanz Subirana et al. (2013)
Relativistic effects
Highly accurate timing capabilities arethe cornerstone of the GNSS positioningtech- nique.Thereforesatellites carryatomic clocks on board, which provide highlyaccurate timestamps forsignals leaving the satellite. Furthermore also on the receiver side, espe- cially forhigh-grade (geodetic) receivers, accuratetiming using precise clocks need to be ensured best as possible. The synchronisation of satellite and receiver clock denotes afundamentalrequirement forGNSSpositioning. Thistime transfer between satellite and receiver is influenced by relativistic effects, which can be split up intotwo parts (as originally defined by Einstein):
• General Relativityeffect: The gravitationalfieldofthe Earthinfluences theclock frequency differently forclocks in space and on Earth.
• SpecialRelativityeffect: The travelspeed of the satellites relativeto users/receivers manifests in aclock shift (time dilatation effect).
The combined effect results in ashift of the center frequency f0 given by: Δf f − f v2 ΔU v2 GM 1 1 = 0 = + = + E ( − ) (2.6) 2 2 2 2 f0 f0 2c c 2c c rs rE v denotes the satellite clock velocityrelativetothe ground receiver, c thespeed of light and ΔU the difference in gravitational potential between satellite and receiver. 2.3Precise PointPositioning 17
Integrating the aboveequationresultsinaclock correction of 38.6microseconds per day(Karabatic 2011). It is worth mentioninghere that the two relativisticeffects counteract each other.Due to GeneralRelativity, clocks on boardofthe satellites runfaster with increasingaltitude (approximately 45 microsecondsper day). SpecificRelativityonthe other handstates that the (relatively) moving clocks on boardofthe satellites run slowerthan their companions on ground. The combined effect is known as the ’factory offset’ which is appliedtothe satellite clocks before the satellitelaunch.
Another related effecttobeconsidered in the frameworkofhigh-precision appli- cations is the Sagnac effect.Itisaresult of the Earth’srotation takingplace while the GNSSsignalistraveling fromsatellitetoreceiver. The travel time correctioncan be calculated as v · ( r − r ) δt = R S (2.7) c 2 where v is thegeocentric velocityvectorofthe receiver and r R and r S arethe geocen- tricvectors ofreceiver and satellite. The correction due to the Sagnac effect can be considered software-internally in PPP processing.
Multipath
The Multipath effect occurswhen aGNSS signalisreflected off an object,such as the wall of abuilding before it is received at the antenna.These reflected signals can distort the code correlation peak andresult in arangingerror. Multipathdenotesone of the dominant error sources in GNSS positioning,especially in heavily obstructed (e.g.urban) areas. This is worsened by thefactthat modelling its influence on measurementsisa very challenging taskand completeelimination of theeffect is still practicallyimpossible forhighly dynamic environments(autonomousdriving,railway). However advancements in antenna designand specific filteralgorithms havemadeiteasier to recognizeand correct especially Non-Line-of-sight multipath. Multipath is ahighly localized error whichdoesnot cancel through differencing and therefore must be considered as themajorsource of errorinmany high precision ap- plications.Iteffects both code andcarrier phase measurements, wherethe introduced errorscan reach up to 150 meters forcodeobservations and up to 5cmfor carrierphases. Moreover, multipath can prevent coordinate convergence in PPPprocessing. Therefore it denotes one the most challenging error sources in GNSS processing, especiallyfor dynamic applications. 18 GlobalNavigation SatelliteSystems
Antenna phase center offsets
GNSS measurements refer to the phase centersofthe antenna of satelliteand receiver respectively.Offsets and variations in the locationofthe phase centerstherefore directly influence the accuracy of processing results. The effects to be considered are:
• Satellite antenna phase center (APC) offset and variation: This holdsfor orbital information in the broadcastnavigationmessages. However, precise orbitsinstead arereferred to the center of mass of the satellite, which makes acorrection to the phase center necessary forprecisePPP applications.
• Receiver APC offset and variation: As forthe satelliteantenna,the measurements arealso referred to the phase center of the receiverantenna.Since this location is not asuitable reference, apointtied to the baseofthe antenna,denoted as antenna reference point (ARP),isused instead.Information on the antenna phase center offset between theAPC and ARPnormallyare provided by the manufacturers or can be obtained from the IGS. Usually the receivercoordinates arereferred to a monument marker or to an external bench mark. The computed coordinates in the IGS SINEX files arereferred to the monumentmarker. Thesite eccentricity vectorbetween the ARP and the monument marker is also given inthe SINEX file, in up,north and east coordinates.
Phasewind-up correction
The phase wind-upeffect affects only the phase measurements and depends on the mutual orientationofthe satelliteand the receiver antennas. Phase wind-up arises due to the nature of the GPS signals as right-hand circularly polarized radio waves.For a static receiver, likethe receivers of areference network, thewind-up effect is due to the satellite orbitalmotion.Usually satellites undergo onlyslowrotationsinorderto keep their solarpanels pointing to the sun direction. This causes aphase variation, whichismisunderstood as arange variation.However,the satellitesalso undergo rapid rotationsinthe case of eclipse and noon seasons. During such an event the rotation of the antenna can reach up to one revolution, which corresponds to an increasein the carrier phase measurement of one cycle, within less than half an hour.For double- difference positioning on distances less thanfew hundreds of kilometres the carrier phase wind-up effect can be neglected. However, forlongerbaselines and in undifferenced point-positioning, thiseffect has to be taken intoaccount.Neglecting the effect and fixing orbits and clocks willleadtopositioning errors in thedm-level (Kouba andH´e roux 2001). 2.3Precise PointPositioning 19
Following Karabatic (2011), the computation of the phase wind-up can be carried out using the formulaspresentedbyWu et al. (1992):
D · D δφ = sign( )arccos( ), (2.8) |D ||D | where = k·(D ×D ), k denotes the satellite-receiver unit vectorand D , D the effective dipole vectors of satelliteand receiver. Thelattercan be derivedfromthe satellite’s (x , y , z )and receiver’s unit vectors ( x, y, z)intheir respectivebodycoordinate system:
D = x − k( k · x ) − k × y , (2.9) D = x − k ( k · x ) − k × y . (2.10)
2.3.3 LinearCombinations of GNSS Observations
In GNSS data processing it is common practicetocombine simultaneous observations us- ing linearcombinations(LC). Different typesofcombinationsexistfor different purposes, the most important ones areintroducedinthe following.Inthe subsequent formulas observations ontwo frequencies areassumed. The equations can be generalized to any pair of signals at different frequencies offered by theGNSS systems.
Ionosphere-free linearcombination
Undoubtedly the most importantLCinGNSSprocessingisthe Ionosphere-free linear combination (IF-LC), whichbuildsthe basis formost processing strategies. It utilizes the dispersive nature of the ionosphereinorder to mitigate its first-order effects (up to 99.9%)byacombination of two-frequency carrier phase (orcode) observations. For carrier phase it reads 1 L = (f2L − f2L ), (2.11) IF 2 2 1 1 2 2 f1 − f2 and forpseudorange observations 1 P = (f2P − f2P ). (2.12) IF 2 2 1 1 2 2 f1 − f2
Especially LIF denotes the commonobservable in (DF) data processingtechniques like PPP. However, there arealso problems encountered in GNSS processing when using this LC. The two most common ones are:
• Difficult ambiguity fixing 20 GlobalNavigation SatelliteSystems
• Highernoiselevel: The IF-LC shows an increased noiselevel compared to the
originalcarrier frequencies (∼ factor3comparedtoL1). Further detailsonthis aregiven in Section 2.3.3.
Geometry-free linearcombination
As acounterpart of theIF-LC,the geometry-free LC (GF-LC or L4)exposes theiono- spheric delaybyremoving all non-dispersiveparts. ThereforethisLCisoften used for the estimation of ionosphere models, as also in this study forthe SEIDalgorithm (see Chapter 6). The GF-LC is formed by thedifference between twomeasurements on dif- ferent frequencies, with allnon-frequency dependent errors, such as satellite, receiver clock and geometry(orbit and station coordinates) cancelling. For phase observations, the GF-LC can be formulated as:
LGF = L4 = L1 − L2, (2.13) and forpseudorangeobservations it reads
PGF = P4 = P2 − P1. (2.14)
As mentioned above, the GF-LC extracts the ionospheric delayand therefore the basis of many ionospheric modelling algorithms.
Wide-lane linearcombination
The Wide-lane LC (LWL)isprimarily designed to construct an observationwith alarger wavelength, whichcan be beneficial forcycleslip detectionand ambiguityfixing. It is derived by f1L1 −f2L2 LWL = (2.15) f1 −f2
Narrow-lane linearcombination
The Narrow-lane LC (LNL)provides anobservation with shorter wavelength andalower noise level compared to the originalcarriers. It allows to increase precision in ambiguity fixed processing (M¨oller 2017).
f1L1 + f2L2 LNL = (2.16) f1 + f2
Melbourne-W¨ubenna linearcombination
The Melbourne-W¨u bbena LC (LMW), proposedbyMelbourne (1985), can be seen as an extension to the WL-LC. In addition to the phase observations forming the first 2.3Precise PointPositioning 21
Observable Wavelength Noise level Multipath level
P1,2 19.0/24.2 cm σcode = ±15cm mpcode = ±150m
L1 19.0 cm σL1 = ±1 − 2mm mpL1 = ±5cm
L2 24.2 cm 1.3 · σL1 1.3 · mpL1
LIF 0.6 cm 3.2 · σL1 4.5 · mpL1
LGF ∞ 1.6 · σL1 2.3 · mpL1
LWL 86.2 cm 6.4 · σL1 9.1 · mpL1
LNL 10.7 cm 0.8 · σL1 1.1 · mpL1
LMW 86.2 cm ∼ σcode 9.1 · mpL1
Table 2.4: Overview on noise andmultipathcharacteristicsofGNSS observables, most infor- mation takenfrom M¨oller (2017). part of the equation, code observationsare used in the second part to mitigatethe ionospheric effects. Thus, the LMW is free from ionospheric delays,geometry, clock, and tropospheric effects. Thedisadvantage of theinclusion of code observations in the LC is thesignificant increase of the observation noise and the multipatheffects. However, in thefuture, advanced code modulations withsignificantly improved noise characteristics will considerably reduce the noiseofthis combination.The Melbourne-W¨u bbena LC is formed as follows:
1 1 L = (f L − f L ) − (f P + f P ) (2.17) MW 2 2 1 1 2 2 1 1 2 2 f1 − f2 f1 + f2
Expectfor multipath and measurement noise, LMW covers solely the wide-lane ambiguities. Therefore it can be used forcycleslip detection, which is exploited in the pre-processing of reference station datafor the SEIDalgorithm as described in Chapter 6.Inthe caseoflow multipath, lowcodenoiseand real-value widelane ambiguities, the LMW is directly used forthe resolutionofwide-lane ambiguities (Sch¨onemann 2014).
Anumber of other useful LCscan be found in GNSS literature, e.ginGlaner (2017).
Noise characteristics
As already seen before forthe IF-LC, some LCsshowundesirable effectsbyintroducing a higher noise level in the data processing compared to observations on theoriginal carrier frequencies (L1/P1 and L2/P2). Thereforeanoverview on thenoise characteristics of different combinations ofGNSS observationsisgiven in Table 2.4. 22 GlobalNavigation SatelliteSystems
2.3.4Integer ambiguityresolution
Although integer ambiguityestimation is not carried out within this study it should still be mentioned hereasakeyingredient forhigh-accuracy applications. If centimetre- accuracy level is required forthe position solution, fixing of the phase ambiguities to integer values is essential.For PPP processing, integer resolution is achieved by ap- plyingimproved satellite products, where phase delays areseparated from theinteger ambiguities (Hinterberger2016). In general two approaches can be distinguished:
• Ge et al. (2008): Decompositionofundifferenced ionosphere-free ambiguities into wide-lane and narrow-lane part.
• Laurichesse et al. (2009): Undifferenced ambiguities were directlyfixed to integers.
For adetailed discussion of thoseapproaches and furtherconsiderations on integer am- biguityresolution see e.g. Hinterberger (2016)orGlaner et al. (2020).
2.3.5Kinematic PPP
Kinematic PPPisapromisingprocessingtechnique foranumber of applications which do not allow forthe assumption of astatic roverposition,asdeployedin(standard) static PPP. Kinematic GPS PPP wasfirst investigated by Kouba and H´eroux (2001)and found well suited forpositioning of moving vehicles as needed e.g. in autonomous driving or agricultural applications. Although reachable accuracy is still lowerthanfor techniques likeRTK positioning, the simplicityofthe approach in terms of data requirements still is apromisingfeaturefor alot of users. Overthe last decade, different applications and processing strategies (e.g.concerning GNSS systems or orbit/clock products used) havebeen investigated. Alkan et al. (2017)used kinematic GPS/GLONASS PPP formarine applicationsand found an accuracy level of centimetre- to decimetre-level even in challenging conditions as well as improvementsinsolution availabilitythrough theusage of GLONASSobser- vations. Katsigiannietal. (2019)examined the performance of the Galileo-only PPP and PPP-AR solutionsinpost-processing kinematic mode. Resultsindicate millimetre- to centimetre-level accuracy of the PPPsolution, especially forPPP-AR processing. Yu andGao (2017)studied the benefit of using multi-GNSSobservationsand found an significant enhancement in availabilityand accuracy of the positioningsolution. More- over they also found an increase in convergence performancefor all three coordinate components. Chapter 3
Physicsofthe atmosphere
This sectionaims to give an overviewonthe principles and parameters,which areex- ploited to drawconclusions on thecurrent stateofthe atmosphere. In order to get a first glance of the medium that is investigated in the course of thisstudy,somegeneral information is provided here.
3.1 General overview
The term Atmosphere referstothe gaseous envelopeofthe Earth or anyother celestial body(M¨oller 2017). Theatmosphere of theEarth is theresult of volcanicfumigation approximately four billion yearsago.Without this formation process, life would not have been able to establish itselfonthe planet due to anumber of reasons. This layer(or at leastits bottom end) is where weather andclimatic phenomenatakeplace, which areof larger interest thanever in achanging climate as experienced today. In the following,somebasicinformation on optical properties, chemicalcomposition andvertical structure of theatmosphere is summarized. This short overview is strongly based on Wallace and Hobbs (2006), whoprovide adetailed description not only on these introduction topics but on all majorfields of atmospheric research in general.
3.1.1 Optical properties and energy budget
The optical properties ofthe Earth’s atmosphere can briefly be summarized by two main properties:
• Transparency forincoming solar(short-wave) radiation
• Opaqueness foroutgoing (long-wave) radiation emitted by theEarth’s surface
23 24 Physicsofthe atmosphere
Process Radiation fraction Solarirradiation (short-wave) absorbed at surface +161 W/m2 Thermal emissionbythe Earth’ssurface -396 W/m2 Long-wave downwardradiation (e.g. fromclouds) +333 W/m2 Net radiation +98 W/m2
Table 3.1: Net radiation budgetofthe Earth, M¨oller (2017)
These propertiesmanifest themselves in awarming of the Earth’ssurface compared to temperatures to be expected from classical(blackbody) radiative transfer theory.This behaviour is what is commonly referred to as the Greenhouse Effect. Considerations on the energy budget of thewhole Earth system aretherefore closely linked to theseproperties.The derivation of thesystems energy budget relies on the conceptofblackbodyradiation.The termblackbodydescribes asurface whereall incoming energy is absorbed.Its flux density(i.e. irradiance) F is given by the Stefan- Boltzmannlaw
F = σT 4, (3.1) with σ =5.67 · 10−8 Wm−2K−4 being the Stefan-Boltzmann constant and T the temperatureofthe surface. Inserting the relevant parametersfor theEarth givesthe netradiationbudget outlined in Table 3.1.The remaining excess of 98 W/m2 in the budget is compensated by non- radiativeprocesses likesensible (e.g. convection) or latent heatfluxes (e.g. melting processes).
3.1.2Atmospheric composition
The Earth’s atmosphere is composed of amixtureofdifferent gases. The dominant ones include hydrogen (H), carbon dioxide (CO2)and methane (CH4). Togetherwith ozone (O3), they denote the most important greenhouse gases.This means that their molecularstructure allowsthem to effectivelytrapoutgoing long-wave radiation, and therefore cause the greenhouse effect. The dominant gasesinterms of concentration in the atmosphere arediatomic nitrogen
(N2)and oxygen(O2). Togetherthey account forover 99 %offractional concentration by volume. Adetailed overviewofthe major constituents of the atmospherecan be found in Table 3.2. 3.2Vertical structure andlayersofthe atmosphere 25
Constituent Molecularweight Fractional concentration by volume
Nitrogen(N2)28.013 78.08 %
Oxygen(O2)32.000 20.95 % Argon (Ar) 39.950.93 %
Water vapour (H2O) 18.02 0-5%
Carbon dioxide (CO2) 44.01 380 ppm Neon (Ne) 20.18 18 ppm Helium (He) 4.00 5ppm
Methane (CH4) 16.04 1.75ppm
Nitrous oxide (N2O) 56.03 0.3 ppm
Ozone (O3) 56.03 0.3ppm
Table 3.2: Fractional concentrations of gaseousatmospheric constituents with respect to dry air,withimportant greenhouse gases marked bold. Takenfrom Wallaceand Hobbs (2006)
3.2 Verticalstructure and layers of the atmosphere
This section discussesthe vertical distributionofthree essentialatmosphericvariables: pressure, temperatureand watervapour. Thereforethe concept of hydrostaticbalance is introduced which allowsfinding the relationship withair density andatmosphericpres- sure. Basedonthe described distributions (mainlyonthat of temperature), distinctive layers areintroduced in which the atmosphereiscommonly divided. Due to its importance forthis study,watervapourand itsdistribution and characteristics arediscussed separately in Section 3.3.
3.2.1 Hydrostaticbalance
Assuming astaticatmosphere, pressureonevery vertical level only depends on the weight of the (fluid) layers aboveit. This is what is commonlyreferred to as hydrostatic balance.Imagineavertical columnofair as shown in Figure 3.1.Itisobvious that both pressure p and density ρ are(at least) dependent on height z.For asmall changein height δz the corresponding pressurechange δp is
∂p δp = δz. (3.2) ∂z
The assumptionofastatic(i.e not accelerating)vertical column implies that the net force on it must be zero. This netforce consistsofthe following vertical forces (negative signindicatesdownwardforce) 26 Physicsofthe atmosphere
Figure 3.1: Avertical columnofair of density ρ,horizontal cross-sectional area δA, height δzand mass M= ρδAδz. The pressure on thelower surface is p,the pressure on theupper surface is p+δp.
• gravitational force Fg = −gM = −gρδAδz,
• pressureforce acting at the topface, FT = −(p + δp)δA,and
• pressureforce acting at the bottom face, FB = pδA.
Setting the net force Fg +FT +FB =0oneobtains δp+gρδz =0andbyusingEquation 3.2 finds ∂p + gρ =0, (3.3) ∂z theequation of hydrostatic balance.Itdescribes the decrease in atmosphericpres- sure withheight, i.e with decreasing weight of theoverlyinglayers. The assumptionof negligiblevertical acceleration of the air columnisanappropriate oneinmost circum- stances. Exceptions aresmall-scale phenomena likeconvection, tornados, or turbulence in atmosphericboundary layer(ABL).
3.2.2Vertical structure of pressure/density
An importantinformation Equation 3.3 withholds is the actualvalueofp(z).Therefore the equation of state of airisused to rewriteEquation 3.3 to ∂p gp = − , (3.4) ∂z RT where T is theatmospheric temperatureand R=8314.459 J/(kmolK) the idealgas constant. 3.2Vertical structure andlayersofthe atmosphere 27
For an isothermal atmosphere (T = T0 = const.)Equation 3.4 reads ∂p gp p = − = − , (3.5) ∂z RT0 H where theconstant H = RT0/g is calledthe scale height.The resulting pressure p at height z then is (−z) p(z)=pse H , (3.6) where ps = p(z =0)is thepressure at thesurface z =0.Equation 3.6 describesthe exponentialdecrease of pressurewith height. The isothermal assumption delivers quite reasonable results forthe structurewhich is actuallyobserved. In the more realistic case of an non-isothermal atmosphere,where temperature T (z) andscale height H(z) areafunction of height too, pressure p(z) is foundas