A demonstration of high precision GPS orbit determination for geodetic applications

High precision orbit determination of Global Positioning System (GPS) satellites is a key requirement for GPS-based precise geodetic measurements and precise low-earth orbiter tracking, currently under study at JPL. Different strategies for orbit determination have been explored at JPL with data from a 1985 GPS field experiment. The most successful strategy uses multi-day arcs for orbit determination and includes fine tuning of spacecraft solar pressure coefficients and station zenith tropospheric delays using the GPS data. Average rms orbit repeatability values for 5 of the GPS satellites are 1.0, 1.2, and 1.7 m in altitude, cross-track, and down-track componenets when two independent 5-day fits are compared. Orbit predictions up to 24 hours outside the multi-day arcs agree within 4 m of independent solutions obtained with well tracked satellites in the prediction interval. Baseline repeatability improves with multi-day as compared to single-day arc orbit solutions. When tropospheric delay fluctuations are modeled with process noise, significant additional improvement in baseline repeatability is achieved. For a 246-km baseline, with 6-day arc solutions for GPS orbits, baseline repeatability is 2 parts in 100 million (0.4-0.6 cm) for east, north, and length components and 8 parts in 100 million for the vertical component. For 1314 and 1509 km baselines with the same orbits, baseline repeatability is 2 parts in 100 million for the north components (2-3 cm) and 4 parts in 100 million or better for east, length, and vertical components

Estimation of tropospheric wet delay from GNSS measurements

The determination of the zenith wet delay (ZWD) component can be a difficult task due to the dynamic nature of atmospheric water vapour. However, precise estimation of the ZWD is essential for high-precision Global Navigation Satellite System (GNSS) applications such as real-time positioning and Numerical Weather Prediction (NWP) modelling.The functional and stochastic models that can be used for the estimation of the tropospheric parameters from GNSS measurements are presented and discussed in this study. The focus is to determine the ZWD in an efficient manner in static mode. In GNSS, the estimation of the ZWD is directly impacted by the choice of stochastic model used in the estimation process. In this thesis, the rigorous Minimum Norm Quadratic Unbiased Estimation (MINQUE) method was investigated and compared with traditional models such as the equal-weighting model (EWM) and the elevationangle dependent model (EADM). A variation of the MINQUE method was also introduced. A simulation study of these models resulted in MINQUE outperforming the other stochastic models by at least 36% in resolving the height component. However, this superiority did not lead to better ZWD estimates. In fact, the EADM provided the most accurate set of ZWD estimates among all the models tested. The EADM also yielded the best ZWD estimates in the real data analyses for two independent baselines in Australia and in Europe, respectively.The study also assessed the validity of a baseline approach, with a reduced processing window size, to provide good ZWD estimates at Continuously Operating Reference Stations (CORS) in an efficient manner. Results show that if the a-priori station coordinates are accurately known, the baseline approach, along with a 2-hour processing window, can produce ZWD estimates that are statistically in good agreement with the estimates from external sources such as the radiosonde (RS), water vapour radiometer (WVR) and International GNSS Service (IGS) solutions. Resolving the ZWD from GNSS measurements in such a timely manner can aid NWP model in providing near real-time weather forecasts in the data assimilation process.In the real-time kinematic modelling of GNSS measurements, the first-order Gauss- Markov (GM) autocorrelation model is commonly used for the dynamic model in Kalman filtering. However, for the purpose of ZWD estimation, it was found that the GM model consistently underestimates the temporal correlations that exist among the ZWD measurements. Therefore, a new autocorrelation dynamic model is proposed in a form similar to that of a hyperbolic function. The proposed model initially requires a small number of autocorrelation estimates using the standard autocorrelation formulations. With these autocorrelation estimates, the least-squares method is then implemented to solve for the model’s parameter coefficients. Once solved, the model is then fully defined. The proposed model was shown to be able to follow the autocorrelation trend better than the GM model. Additionally, analysis of real data at an Australian IGS station has showed the proposed model performed better than the random-walk model, and just as well as the GM model. The proposed model was able to provide near real-time (i.e. 30 seconds interval) ZTD estimates to within 2 cm accuracy on average.The thesis also included an investigation into the several interpolation models for estimating missing ZWD observations that may take place during temporary breakdowns of GNSS stations, or malfunctions of RS and WVR equipments. Results indicated marginal differences between the polynomial regression models, linear interpolation, fast-Fourier transform and simple Kriging methods. However, the linear interpolation method, which is dependent on the two most recent data points, is preferable due to its simplicity. This result corresponded well with the autocorrelation analysis of the ZWD estimates where significant temporal correlations were observed for at most two hours.The study concluded with an evaluation of several trend and smoothing models to determine the best models for predicting ZWD estimates, which can help improve real-time kinematic (RTK) positioning by mitigating the tropospheric effect. The moving average (MA) and the single-exponential smoothing (SES) models were shown to be the best-performing prediction models overall. These two models were able to provide ZWD estimates with forecast errors of less 10% for up to 4 hours of prediction

Kinematic GNSS tropospheric estimation and mitigation over a range of altitudes

PhD ThesisThis thesis investigates the potential for estimating tropospheric delay from Global Navigation Satellite Systems (GNSS) stations on moving platforms experiencing a change in altitude. The ability to accurately estimate tropospheric delay in kinematic GNSS positioning has implications for improved height accuracy due to the mitigation of a major GNSS error source, and for the collection of atmospheric water vapour data for meteorology and climate studies. The potential for extending current kinematic GNSS positioning estimates of tropospheric delay from sea level based studies to airborne experiments, and the achievable height accuracy from a range of tropospheric mitigation strategies used in airborne GNSS positioning, are explored. An experiment was established at the Snowdon Mountain Railway (SMR), utilising the railway to collect a repeatable kinematic dataset, profiling 950 m of the lower atmosphere over a 50 day period. GNSS stations on stable platforms and meteorological sensors were installed at the extremities of the trajectory, allowing reference tropospheric delays and coordinates to be established. The retrieval of zenith wet delay (ZWD) from kinematic GNSS solutions using tropospheric estimation strategies is validated against an interpolated reference ZWD between GNSS stations on stable platforms, together with profiles from 100 m resolution runs of the UK Met Office Unified Model. Agreement between reference ZWD values and a combined GPS+GLONASS precise point positioning (PPP) solution is demonstrated with an accuracy of 11.6 mm (RMS), similar to a relative positioning solution and previous shipborne studies. The impact on the height accuracy from estimating tropospheric delay in kinematic GNSS positioning is examined by comparing absolute and relative GNSS positioning solutions to a reference trajectory generated from a relative GNSS positioning solution ii processed with reference to the GNSS stations on stable platforms situated at the extremities of the SMR. A height accuracy with a standard deviation of 72 mm was demonstrated for the GPS+GLONASS PPP solution, similar to a GPS-only relative solution, and providing an improvement over the GPS-only PPP solution.UK Natural Environment Research Council (NERC) studentship, and part of the work was funded by the Royal Institution of Chartered Surveyors (RICS) Education Trust

Post-Processing Precise Point Positioning Solutions with Parameter Optimization

Precise Point Positioning (PPP) technique can offer position solutions with centimeter-level accuracy by fusing precise satellite orbits and clocks with un-differenced, dual-frequency, pseudo-range, and carrier-phase observables. PPP presents a compelling alternative to Differential Global Positioning Systems, with the benefit that it only requires a single receiver and does not require simultaneous observations from many stations, making it appealing for ongoing research on hydro-graphic survey applications. The National Oceanic and Atmospheric Administration has been working on a buoy system tracked with a Global Positioning Systems receiver and Inertial Measurement Unit sensor using the PPP technique. In the interest of obtaining accurate measurements, this data is post-processed using a software package for position navigation with tight Inertial Navigation System capabilities developed by the Jet Propulsion Laboratory, this is GNSS Inferred Positioning Systems (GIPSYx). GIPSYx Software allows finely controllable user inputs for selectable models and configurations. This flexibility allows fitting the right models for different data sources but requires a tuning process to find suitable configurations. A processing strategy for buoy data with GIPSYx positioning software is described and a method to assess solutions to automatically optimize the process of finding these manually tuned model parameters is provided. Other data sources are considered to generalize this method and prove the concept of optimizing positioning software configurations from output solution evaluation using black-box optimization

Ionospheric tomographic common clock model of undifferenced uncombined GNSS measurements

This is a post-peer-review, pre-copyedit version of an article published in Journal of geodesy. The final authenticated version is available online at: http://dx.doi.org/10.1007/s00190-021-01568-8In this manuscript, we introduce the Ionospheric Tomographic Common Clock (ITCC) model of undifferenced uncombined GNSS measurements. It is intended for improving the Wide Area precise positioning in a consistent and simple way in the multi-GNSS context, and without the need of external precise real-time products. This is the case, in particular, of the satellite clocks, which are estimated at the Wide Area GNSS network Central Processing Facility (CPF) referred to the reference receiver one; and the precise realtime ionospheric corrections, simultaneously computed under a voxel-based tomographic model with satellite clocks and other geodetic unknowns, from the uncombined and undifferenced pseudoranges and carrier phase measurements at the CPF from the Wide Area GNSS network area. The model, without fixing the carrier phase ambiguities for the time being (just constraining them by the simultaneous solution of both ionospheric and geometric components of the uncombined GNSS model), has been successfully applied and assessed against previous precise positioning techniques. This has been done by emulating real-time conditions for Wide Area GPS users during 2018 in Poland.Peer ReviewedPostprint (published version

센티미터 급 광역 보강항법 시스템의 반송파 위상 기반 보정정보 생성 알고리즘에 관한 연구

학위논문(박사)--서울대학교 대학원 :공과대학 기계항공공학부,2020. 2. 기창돈.Recently, the demand for high-precision navigation systems for centimeter-level service has been growing rapidly for various Global Navigation Satellite System (GNSS) applications. The network Real-Time Kinematic (RTK) is one of the candidate solution to provide high-accuracy position to user in real-time. However, the network RTK requires a lot of reference stations for nationwide service. Furthermore, it requires high-speed data-link for broadcasting their scalar-type corrections. This dissertation proposed a new concept of satellite augmentation system called Compact Wide-Area RTK, which provides centimeter-level positioning service on national or continental scales to overcoming the limitation of the legacy network RTK methods. Using the wide-area network of multiple reference stations whose distance is 200~1,000 km, the proposed system generates three types of carrier-phase-based corrections: satellite orbit corrections, satellite code/phase clock (CPC) corrections, tropospheric corrections. Through the strategy of separating the scalar-type corrections of network RTK into vector forms of each error component, it is enable to expand network RTK coverage to continental scale using a similar number of reference stations as legacy meter-level Satellite-Based Augmentation System (SBAS). Furthermore, it is possible to broadcast their corrections over a wide-area using geosynchronous (GEO) satellite with extremely low-speed datalink of 250 bps likewise of legacy SBAS. To sum up, the proposed system can improve position accuracy by centimeter-level while maintaining the hardware infrastructure of the meter-level legacy SBAS. This study mainly discussed on the overall system architecture and core algorithms for generating satellite CPC corrections and tropospheric corrections. This study proposed a new Three-Carrier Ambiguity Resolution (TCAR) algorithm using ionosphere-free combinations to correctly solve the integer ambiguity in wide-area without any ionospheric corrections. The satellite CPC corrections are calculated based on multiple stations for superior and robust performance under communication delay and outage. The proposed algorithm dramatically reduced the latency compensation errors and message amounts with compare to conventional RTK protocols. The tropospheric corrections of the compact wide-area RTK system are computed using GPS-estimated precise tropospheric delay and weather data based model together. The proposed algorithm adopts spherical harmonics function to significantly reduce the message amounts and required number of GPS reference stations than the network RTK and Precise Point Positioning-RTK (PPP-RTK), while accurately modeling the spatial characteristic of tropospheric delay with weather data together. In order to evaluate the user domain performance of the compact wide-area RTK system, this study conducted the feasibility test on mid-west and south USA using actual GPS measurements. As a result, the 95% horizontal position error is about 1.9 cm and the 95% vertical position error is 7.0 cm after the integer ambiguity is correctly fixed using GPS-only signals. The user ambiguity resolution takes about 2 minutes, and success-fix rate is about 100 % when stable tropospheric condition. In conclusion, the compact wide-area RTK system can provide centimeter-level positioning service to wide-area coverage with extremely low-speed data link via GEO satellite. We hope that this new system will consider as candidate solution for nationwide centimeter-level service such as satellite augmentation system of the Korea Positioning System (KPS).최근 자율주행자동차, 무인 드론 배송, 충돌 회피, 무인트랙터를 이용한 스마트 무인 경작 등 위성항법시스템(GNSS, Global Navigation Satellite System)을 사용하는 다양한 응용분야에서 수 cm 수준의 정밀 위치 정보에 대한 요구가 급격히 증가하고 있다. 본 학위논문에서는 1 m 급의 정확하고 신뢰성 높은 위치 서비스를 제공하는 기존의 정지궤도위성 기반 광역 보강항법 시스템(SBAS, Satellite-Based Augmentation System)의 기준국 인프라를 유지하면서 항법 성능을 수 cm 수준으로 향상시키기 위해 반송파 위상 기반의 초정밀 보정정보 생성 알고리즘에 관한 연구를 수행하였다. 실시간 정밀 측위(RTK, Real-Time Kinematic)는 반송파 위상 측정치에 포함된 미지정수를 정확하게 결정하여 수 cm 수준의 정밀 항법 서비스를 가능하게 하는 대표적인 기법이다. 그 중에서도 약 50~70 km 간격으로 분포된 다수의 기준국 정보를 활용하는 Network RTK 기법은 동적 사용자의 빠르고 정확한 위치 결정이 가능한 인프라로서 주목받고 있다. 하지만 스칼라 형태로 구성된 Network RTK 보정정보는 각 기준국 별로 관측된 위성 수에 따라 생성이 되기 때문에 보정 데이터 량이 상당히 방대하다. 메시지 전송에 필요한 데이터 량이 많을수록 고속의 통신 환경을 필요로 하며, 메시지 시간 지연이나 통신 단절에 매우 취약한 문제를 가지고 있다. 또한 스칼라 형태의 보정정보는 사용자와 기준국 간의 거리가 멀어질수록 보정 오차가 크게 발생하기 때문에 대륙 혹은 나라 규모의 광역에서 서비스하기 위해서는 수십~수백 개 이상의 기준국 인프라 구축이 필수적이다. 예를 들어, SBAS가 한반도 지역 서비스를 위해 5~7개의 기준국이 필요한 반면 Network RTK는 90~100개의 기준국이 필요하다. 즉 Network RTK는 시스템 구축 및 유지 비용이 SBAS 대비 약 15배 정도 많이 들게 된다. 본 논문에서는 기존 Network RTK의 문제점을 해결하기 위한 방법으로 대륙 급 광범위한 영역에서 실시간으로 cm급 초정밀 위치결정 서비스 제공이 가능한 Compact Wide-Area RTK 라는 새로운 개념의 광역보강항법시스템 아키텍처를 제안하였다. Compact Wide-Area RTK는 약 200~1,000 km 간격으로 넓게 분포된 기준국 네트워크를 활용하여 반송파 위상 기반의 정밀한 위성 궤도 보정정보, 위성 Code/Phase 시계 보정정보, 대류층 보정정보를 생성하는 시스템이다. 기존 스칼라 형태의 Network RTK 보정정보 대신 오차 요소 별 벡터 형태의 정밀 보정정보를 생성함으로써 데이터 량을 획기적으로 절감하고 서비스 영역을 확장할 수 있다. 최종적으로 SBAS와 마찬가지로 250 bps의 저속 통신 링크를 가진 정지궤도위성을 통해 광역으로 보정정보 방송이 가능하다. 본 논문에서는 3가지 보정정보 중 위성 Code/Phase 시계 보정정보와 대류층 보정정보 생성을 위한 핵심 알고리즘에 대해 중점적으로 연구하였다. 반송파 위상 기반의 정밀 보정정보 생성을 위해서는 먼저 미지정수를 정확하게 결정해야 한다. 본 논문에서는 삼중 주파수 반송파 위상 측정치의 무-전리층 조합을 활용하여 전리층 보정정보 없이도 정확하게 미지정수 결정 가능한 새로운 방법을 제안하였다. 위성 Code/Phase 시계 보정정보는 통신 지연 및 고장 시 우수하고 강건한 성능을 위해 다중 기준국의 모든 측정치를 활용하여 추정된다. 이 때 각 기준국 별 서로 다른 미지정수 때문에 발생하는 문제는 앞서 정확하게 결정된 기준국 간 이중차분 된 미지정수를 활용하여 수준을 조정하는 과정을 통해 해결이 가능하다. 그 결과 생성된 위성 Code/Phase 보정정보 메시지의 크기, 변화율, 잡음 수준이 크게 개선되었고, 통신 지연 시 오차 보상 성능이 기존 RTK 프로토콜 보다 99% 향상 됨을 확인하였다. 대류층 보정정보는 적은 수의 기준국 만을 활용하여 정확하게 대류층을 모델링하기 위해 자동 기상관측시스템으로부터 수집한 기상 정보를 추가로 활용하여 생성된다. 본 논문에서는 GNSS 기준국 네트워크로부터 정밀하게 추정된 반송파 위상 기반 수직 대류층 지연과 기상정보 기반으로 모델링 된 수직 대류층 지연을 함께 활용할 수 있는 새로운 알고리즘을 제안하였다. 구면조화함수를 사용하여 Network RTK 및 PPP-RTK 보다 필요한 메시지 양과 기준국 수를 크게 감소시키면서도 RMS 2 cm 수준으로 정확한 보정정보 생성이 가능함을 확인하였다. 본 논문에서 제안한 Compact Wide-Area RTK 시스템의 항법 성능을 검증하기 위해 미국 동부 지역 6개 기준국의 실측 GPS 데이터를 활용하여 테스트를 수행하였다. 그 결과 제안한 시스템은 미지정수 결정 이후 사용자의 95% 수평 위치 오차 1.9 cm, 95% 수직 위치 오차 7.0 cm 로 위치를 정확하게 결정하였다. 사용자 미지정수 결정 성능은 대류층 안정 상태에서 약 2분 내로 100% 의 성공률을 가진다. 본 논문에서 제안한 시스템이 향후 한국형 위성항법 시스템(KPS, Korean Positioning System)의 전국 단위 센티미터 급 서비스를 위한 알고리즘으로 활용되기를 기대한다.CHAPTER 1. Introduction 1 1.1 Motivation and Purpose 1 1.2 Former Research 4 1.3 Outline of the Dissertation 7 1.4 Contributions 8 CHAPTER 2. Overview of GNSS Augmentation System 11 2.1 GNSS Measurements 11 2.2 GNSS Error Sources 14 2.2.1 Traditional GNSS Error Sources 14 2.2.2 Special GNSS Error Sources 21 2.2.3 Summary 28 2.3 GNSS Augmentation System 29 2.3.1 Satellite-Based Augmentation System (SBAS) 29 2.3.2 Real-Time Kinematic (RTK) 32 2.3.3 Precise Point Positioning (PPP) 36 2.3.4 Summary 40 CHAPTER 3. Compact Wide-Area RTK System Architecture 43 3.1 Compact Wide-Area RTK Architecture 43 3.1.1 WARTK Reference Station (WRS) 48 3.1.2 WARTK Processing Facility (WPF) 51 3.1.3 WARTK User 58 3.2 Ambiguity Resolution and Validation Algorithms of Compact Wide-Area RTK System 59 3.2.1 Basic Theory of Ambiguity Resolution and Validation 60 3.2.2 A New Ambiguity Resolution Algorithms for Multi-Frequency Signals 65 3.2.3 Extra-Wide-Lane (EWL) Ambiguity Resolution 69 3.2.4 Wide-Lane (WL) Ambiguity Resolution 71 3.2.5 Narrow-Lane (NL) Ambiguity Resolution 78 3.3 Compact Wide-Area RTK Corrections 83 3.3.1 Satellite Orbit Corrections 86 3.3.2 Satellite Code/Phase Clock (CPC) Corrections 88 3.3.3 Tropospheric Corrections 89 3.3.4 Message Design for GEO Broadcasting 90 CHAPTER 4. Code/Phase Clock (CPC) Correction Generation Algorithm 93 4.1 Former Research of RTK Correction Protocol 93 4.1.1 Observation Based RTK Data Protocol 93 4.1.2 Correction Based RTK Data Protocol 95 4.1.3 Compact RTK Protocol 96 4.2 Satellite CPC Correction Generation Algorithm 100 4.2.1 Temporal Decorrelation Error Reduced Methods 102 4.2.2 Ambiguity Level Adjustment 105 4.2.3 Receiver Clock Synchronization 107 4.2.4 Averaging Filter of Satellite CPC Correction 108 4.2.5 Ambiguity Re-Initialization and Message Generation 109 4.3 Correction Performance Analysis Results 111 4.3.1 Feasibility Test Environments 111 4.3.2 Comparison of RTK Correction Protocol 113 4.3.3 Latency Compensation Performance Analysis 116 4.3.4 Message Data Bandwidth Analysis 119 CHAPTER 5. Tropospheric Correction Generation Algorithm 123 5.1 Former Research of Tropospheric Correction 123 5.1.1 Tropospheric Corrections for SBAS 124 5.1.2 Tropospheric Corrections of Network RTK 126 5.1.3 Tropospheric Corrections of PPP-RTK 130 5.2 Tropospheric Correction Generation Algorithm 136 5.2.1 ZWD Estimation Using Carrier-Phase Observations 138 5.2.2 ZWD Measurements Using Weather Data 142 5.2.3 Correction Generation Using Spherical Harmonics 149 5.2.4 Correction Applying Method for User 157 5.3 Correction Performance Analysis Results 159 5.3.1 Feasibility Test Environments 159 5.3.2 Zenith Correction Domain Analysis 161 5.3.3 Message Data Bandwidth Analysis 168 CHAPTER 6. Compact Wide-Area RTK User Test Results 169 6.1 Compact Wide-Area RTK User Process 169 6.2 User Performance Test Results 173 6.2.1 Feasibility Test Environments 173 6.2.2 User Range Domain Analysis 176 6.2.3 User Ambiguity Domain Analysis 182 6.2.4 User Position Domain Analysis 184 CHAPTER 7. Conclusions 189 Bibliography 193 초 록 207Docto

Estimating integrated water vapor trends from VLBI, GPS,and numerical weather models: sensitivity totropospheric parameterization

©2018. American Geophysical UnionIn this study, we estimate integrated water vapor (IWV) trends from very long baseline interferometry (VLBI) and global navigation satellite systems (GNSS) data analysis, as well as from numerical weather models (NWMs). We study the impact of modeling and parameterization of the tropospheric delay from VLBI on IWV trends. We address the impact of the meteorological data source utilized to model the hydrostatic delay and the thermal deformation of antennas, as well as the mapping functions employed to project zenith delays to arbitrary directions. To do so, we derive a new mapping function, called Potsdam mapping functions based on NWM data and a new empirical model, GFZ‐PT. GFZ‐PT differs from previous realizations as it describes diurnal and subdiurnal in addition to long‐wavelength variations, it provides harmonic functions of ray tracing‐derived gradients, and it features robustly estimated rates. We find that alternating the mapping functions in VLBI data analysis yields no statistically significant differences in the IWV rates, whereas alternating the meteorological data source distorts the trends significantly. Moreover, we explore methods to extract IWV given a NWM. The rigorously estimated IWV rates from the different VLBI setups, GNSS, and ERA‐Interim are intercompared, and a good agreement is found. We find a quite good agreement comparing ERA‐Interim to VLBI and GNSS, separately, at the level of 75%.DFG, 255986470, GGOS-SIM-2: Simulation des "Global Geodetic Observing System

An assessment of the precise products on static Precise Point Positioning using Multi-Constellation GNSS

Precise point positioning (PPP) is highly dependent on the precise ephemerides and satellite clock products that are used. Different ephemeris and clock products are available from a variety of different organizations. The aim of this paper is to assess the achievable static positioning accuracy and precision when using different precise ephemerides from three analysis centres Natural Resources Canada (EMX), European Space Agency (ESA) and GeoForschungsZentrum (GFZ), using GPS alone, GLONASS alone, and GPS and GLONASS combined. It will be shown in this paper that the precise products are significantly affected by the time-base of the reference stations, and that this is propagated through to all the estimated satellite clocks. In order to overcome the combined biases in the estimated satellite clock, in the PPP processing, these clocks errors need to be handled with an appropriate variation in the estimated receiver clock. It will also be shown that the precise coordinates of the satellites differ between the analysis centres, and this affects the PPP position estimation at the millimetre level. However, all those products will be shown to result in the same level of precision for all coordinate components and are equivalent to the horizontal precision from a Global Double Difference (GDD) solution. For the horizontal coordinate component, the level of agreement between the PPP solutions, and with the GDD solution, is at the millimetre level. There is a notable, but small, bias in the north coordinate components of the PPP solutions, from the corresponding north component of the GDD solutions. It is shown that this difference is due to the different strategy adopted for the GDD and PPP solutions, with PPP being more affected by the changing satellite systems. The precision of the heights of the receiver sites will be shown to be almost the same across all the PPP scenarios, with all three products. Finally, it will be concluded that accuracy of the height component is system dependent and is related to the behaviour of antenna phase centre with the different constellation type

GNSS precise point positioning :the enhancement with GLONASS

PhD ThesisPrecise Point Positioning (PPP) provides GNSS navigation using a stand-alone receiver with no base station. As a technique PPP suffers from long convergence times and quality degradation during periods of poo satellite visibility or geometry. Many applications require reliable realtime centimetre level positioning with worldwide coverage, and a short initialisation time. To achieve these goals, this thesis considers the use of GLONASS in conjunction with GPS in kinematic PPP. This increases the number of satellites visible to the receiver, improving the geometry of the visible satellite constellation. To assess the impact of using GLONASS with PPP, it was necessary to build a real time mode PPP program. pppncl was constructed using a combination of Fortran and Python to be capable of processing GNSS observations with precise satellite ephemeris data in the standardised RINEX and SP3 formats respectively. pppncl was validated in GPS mode using both staticsites and kinematic datasets.In GPS only mode,one sigma accuracy of 6.4mm and 13mm in the horizontal and vertical respectively for 24h static positioning was seen. Kinematic horizontal and vertical accuracies of 21mm and 33mm were demonstrated. pppncl was extended to assess the impact of using GLONASS observations in addi- tion to GPS instatic and kinematic PPP. Using ESA and Veripos Apex G2 satel- lite orbit and clock products,the average time until 10cm 1D static accuracy was achieved, over arange of globally distributed sites, was seen to reduce by up to 47%. Kinematic positioning was tested for different modes of transport using real world datasets. GPS/GLONAS SPPP reduced the convergence time to decimetre accuracy by up to a factor of three. Positioning was seen to be more robust in comparison to GPS only PPP, primarily due to cycle slips not being present on both satellite systems on the occasions when they occurred,and the reduced impact of undetected outliersEngineering and Physical Sciences Research Council, Verip os/Subsea