Next Article in Journal
Incorporating the Plant Phenological Trajectory into Mangrove Species Mapping with Dense Time Series Sentinel-2 Imagery and the Google Earth Engine Platform
Next Article in Special Issue
Integrated Precise Orbit Determination of Multi-GNSS and Large LEO Constellations
Previous Article in Journal
Frost Damage Assessment in Wheat Using Spectral Mixture Analysis
Previous Article in Special Issue
Mitigation of Unmodeled Error to Improve the Accuracy of Multi-GNSS PPP for Crustal Deformation Monitoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

PPP and PPP-AR Kinematic Post-Processed Performance of GPS-Only, Galileo-Only and Multi-GNSS

1
Centre National d’Etudes Spatiales (CNES), Observatoire Midi-Pyrénées (OMP), Géosciences Environnement Toulouse (GET), 14 avenue Edouard Belin, 31400 Toulouse, France
2
Collecte Localisation Satellites (CLS Group), 8-10 rue Hermès, 31520 Ramonville Saint Agne, France
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(21), 2477; https://0-doi-org.brum.beds.ac.uk/10.3390/rs11212477
Submission received: 31 August 2019 / Revised: 14 October 2019 / Accepted: 21 October 2019 / Published: 24 October 2019

Abstract

:
Precise point positioning (PPP) has been used for decades not only for general positioning needs but also for geodetic and other scientific applications. The CNES-CLS Analysis Centre (AC) of the International GNSS Service (IGS) is performing PPP with phase ambiguity resolution (PPP-AR) using the zero-difference ambiguity fixing approach also known as “Integer PPP” (IPPP). In this paper we examine the postprocessed kinematic PPP and PPP-AR using Galileo-only, GPS-only and Multi-GNSS (GPS + Galileo) constellations. The interest is to examine the accuracy for each GNSS system individually but also of their combination to measure the current benefits of using Galileo within a Multi-GNSS PPP and PPP-AR. Results show that Galileo-only positioning is nearly at the same level as GPS-only; around 2–4 mm horizontal and aound 10 mm vertical repeatability (example station of BRUX). In addition, the use of Galileo system—even uncompleted—improves the performance of the positioning when combined with GPS giving mm level repeatability (improvement of around 30% in East, North and Up components). Repeatabilities observed for Multi-GNSS (GPS + GAL) PPP-AR, taking into account the global network statistics, are a little larger, with 8 mm in horizontal and 17 mm in vertical directions. This result shows that including Galileo ameliorates the best positioning accuracy achieved until today with GPS PPP-AR.

Graphical Abstract

1. Introduction

The precise point positioning method (PPP) [1] is a well-known and widely used method for positioning using zero-difference observations. This method is used for calculating the coordinates of a station without the need of a reference station nearby as a control station. It has been used in numerous scientific applications; namely general static positioning [2] for local or global networks, kinematic positioning [3], time transfer [4], etc. The PPP accuracy reaches the highest levels once the carrier phase ambiguities are resolved—the so-called PPP-AR (PPP with ambiguity resolution) or less frequently integer PPP (IPPP). The French Space Agency, Centre National d’Etudes Spatiales (CNES) and the Collecte Localisation Satellites company (CLS Group) Analysis Centre (AC) of the International GNSS Service (IGS) use an undifferenced ambiguity resolution solution for GPS [5] and Galileo [6]. Such a method can be applied also for the determination of precise orbits and integer recovery clocks (IRC). In order to apply this method for PPP and PPP-AR for GPS and Galileo, satellite orbits and IRC products as well as satellite biases have to be consistent. In our study they are taken from the CNES-CLS AC. There are two conditions for this method: for PPP or PPP-AR, orbit and clock products have to be consistent and of the best quality; for PPP-AR IRCs and associated methods, biases need to be available for the products used. For the first condition, any precise orbit like the ones made available by the International GNSS Service (IGS) or Analysis Centers (ACs) in the frame of the Multi GNSS pilot project (MGEX) [7] can be used. For the second we can use the products delivered routinely by the CNES-CLS AC for GPS [8] and Galileo [6].
So far there are already publications giving the possibilities of Galileo positioning [9] or studying the performance the benefit of using multi-GNSS measurements in PPP [7,8]. In a previous study [6] we presented the zero-differences approach used to compute IRC products for Galileo and GPS. Recently we presented the positioning capabilities of Galileo-only solutions [10]. In this article we examine the performance and precision of PPP and PPP-AR in a comparison of Galileo-only, GPS-only and GPS + Galileo AR solutions.
This study is done using the GPS and the incomplete Galileo constellation. As of 2019, the Galileo constellation comprises 26 satellites in total (4 In-Orbit Validation (IOV) and 22 Full Operational Capability (FOC)). Two of the FOC satellites are in eccentric orbits, one IOV is not available and one FOC is not usable. The total number of usable satellites is 22: 3 IOV and 21 FOC satellites [11]. The full constellation is only 4 satellites away and it is scheduled for 2020 [12]. The following hypothesis is examined and verified: if GPS and Galileo measurements are of the same quality, then the simultaneous use of the two systems will give better accuracy than when using each system separately. If the two systems are compatible, then the number of measurements is nearly double than in a single-system solution, and the number of parameters (apart from biases for the different frequencies used) is the same.
This publication is organized in the following sections. Firstly, the undifferenced ambiguity resolution and PPP-AR processing is briefly presented. Section 3 is devoted to the experiments, the processing and the results. Finally, in Section 4, some conclusions are given together with suggestions for further work and perspectives.

2. Materials and Methods

The PPP with ambiguity resolution processing has the advantage that the processing can be done directly at the user level without any reference station around. The requirement, however, is that precise and consistent satellite products (derived from a global network of stations) must be available. In this study products are taken from the CNES-CLS AC: satellite orbit file (.sp3), clock file (.clk) and the Wide-Lane satellite biases (.wsb) (given either in the header of the clock files or taken from the CNES/CLS portal. Files Wide_lane_GAL_satellite_biais.wsb and Wide_lane_GPS_satellite_biais.wsb are available at ftp://ftpsedr.cls.fr/pub/igsac/).
The used zero-difference ambiguity fixing method [8,13] equation model bellow is given for the pseudorange (code) and carrier phase measurements for two frequencies:
P r , i s =   ρ r s + c Δ t + T r s + I r s + b i s + b r , i + E r , i s
P r , j s =   ρ r s + c Δ t + T r s + f i 2 f j 2 I r s + b j s + b r , j + E r , j s
L r , i s = λ i φ r , i s =   ρ r s + c Δ t + T r s I r s + λ i N r , i s + λ i W r s + β i s + β r , i + ε r , i s
L r , j s = λ j φ r , j s =   ρ r s + c Δ t + T r s f i 2 f j 2 I r s + λ j N r , j s + λ j W r s + β j s + β r , j + ε r , j s
where:
  • P r , i s , P r , j s are the code measurement at receiver r from satellite s on frequency i or j (m)
  • L r , i s , L r , j s are the phase measurement at receiver r from satellite s on frequency i or j (m)
  • ρ r s is the geometric distance between receiver and satellite (m)
  • Δ t ( Δ t = δ t r δ t s ) is the clock correction related to the satellite ( δ t s ) and the receiver ( δ t r ) with respect to the synchronization to the GPS time (s)
  • T r s is the troposphere delay (m)
  • I r s is the ionosphere delay (m)
  • E r , i s , E r , j s are the code measurement errors at receiver r from satellite s on frequency i or j (m) including all sources of code errors: multipath and noise.
  • f i , f j are the carrier frequency i or j (Hz)
  • c is the speed of light in vacuum (m/s)
  • λ i , λ j are the nominal wavelength of the carrier frequency i or j (m)
  • φ r , i s , φ r , j s are the carrier phase measurement at receiver r from satellite s on frequency i or j (cycles)
  • N r , i s , N r , j s are the integer carrier phase ambiguity at receiver r from satellite s on frequency i or j
  • W r s is the carrier phase wind up effect (cycles)
  • b s , b r are the code phase biases of satellite and receiver (m)
  • β s , β r are the carrier phase biases of satellite and receiver (m)
  • ε r , i s , ε r , j s are the carrier phase measurement error at receiver r from satellite s on frequency i or j (m) including all sources of phase errors, remaining uncorrected phase center offset and phase center variation, multipath and noise.
In this study the frequencies used for the GPS system are in the band of L 1 and L 2 ( f L 1 = 154 × f 0 and   f L 2 = 120 × f 0 ). For the Galileo system are in the band of E 1 and E 5 a ( f E 1 = 154 × f 0 and   f E 5 a = 115 × f 0 ), where f 0 = 10.23 MHz.
From all these four equations it is possible to form a Melbourne–Wübbena linear combination that has the identity to reduce measurement noise and to cancel out any geometric, ionospheric and clock terms [14,15]:
M W r s =   λ w l φ r s = ( f i f i f j L r , i s f j f i f j L r , j s ) ( f i f i + f j P r , i s + f j f i + f j P r , j s ) = λ w l ( N r , i s   N r , j s μ s + μ r ( t ) ) =   λ w l ( N w l , r s μ s + μ r ( t ) )
where:
  • M W r s is the Melbourne-Wübbena linear combination at receiver r from satellite s (m)
  • λ w l ( λ w l = c / ( f i f j ) =   λ i λ j / ( λ j λ i ) ) is the wide-lane (WL) wavelength (m)
  • N w l , r s ( N w l , r s = N r , i s   N r , j s ) is the WL ambiguity at receiver r from satellite s
  • μ s is the delay coming from the satellite (also known in the bibliography as WL satellite bias (WSB))
  • μ r ( t ) is the delay coming from the receiver (also known in the bibliography as WL receiver bias (WRB))
It has been observed that for the GPS system the μ s are stable over long periods of time and can be considered as constant during at least one day [16]. For the Galileo system they are stable for longer periods; up to months [10]. The μ r ( t ) delay is considered to vary over time because it depends on the behavior of each receiver.
As it is seen from the Equation (5), the terms N w l , r s , μ s and μ r ( t ) are totally correlated. Normally, the μ s is known; i.e., it is calculated and provided from the CNES/CLS AC to the users [6,8]. The μ r (one per epoch of measurements) and the N w l , r s (one per satellite pass) can be separated and solved using all available equations (Equation (5) corrected by the μ s ) from all satellites in view, using a Least Squares Estimation (LSE) processing associated with a bootstrap method [17,18]. The result after this step are the N w l , r s and μ r ( t ) .
Once the N w l , r s is determined, the following step is to form an ionosphere-free linear combination that has the property to cancel out the first order of ionospheric effects. These combinations use the equations in the two frequencies where they apply to them a coefficient: α i to the frequency i and α j to the frequency j respectively [19].
α i =   f i 2 f i 2 f j 2
α j =   f i 2 f i 2 f j 2
Using the above coefficients and Equations (6) and (7), the ionosphere-free linear combinations for code and carrier phase become:
P r , I F s =   α i P r , i s + α j P r , j s = ρ r s + c Δ t + T r s + E r , I F s
L r , I F s =   α i L r , i s + α j L r , j s = ρ r s + c Δ t + T r s + c f i 2 f j 2 ( f i N r , i s f j N r , j s ) + c f i 2 f j 2 ( f i W r s f j W r s ) + ε r , I F s = ρ r s + c Δ t + T r s + λ n l N r , i s + λ n l λ w l λ j N w l , r s + λ n l W r s + ε r , I F s
where:
  • P r , I F s is the ionosphere-free code measurement at receiver r from satellite s (m)
  • E r , I F s is the ionosphere-free code measurement error at receiver r from satellite s (m)
  • L r , I F s is the ionosphere-free carrier measurement at receiver r from satellite s (m)
  • ε r , I F s is the ionosphere-free carrier measurement error at receiver r from satellite s (m)
  • λ n l ( λ n l = c / ( f i + f j ) =   λ i λ j / ( λ i + λ j ) ) is the narrow-lane (NL) wavelength (m)
The following Table 1 gives the wide-lane and the narrow-lane wavelengths for GPS and Galileo frequencies used in this article.
The term W r s is calculated by using the models proposed by Kouba [20]. All equations are then gathered to form another LSE processing. The system of equations for the ionosphere-free code and phase measurements are modelled according to the models given in Table 2 and used to estimate all the parameters: stations position, tropospheric parameters, receiver clocks and inter system bias in case of multi-GNSS measurements. For the PPP case the N r , i s are solved with the other parameters as real values. For the PPP-AR case they are fixed to integer numbers within a bootstrap processing [8,13]. Once the integer N r , i s is determined, the final PPP-AR solution is performed.
Figure 1 is giving a graphic overview of the PPP and PPP-AR method used. In the beginning, the satellite μ s are needed together with the station RINEX file. The pre-processing phase uses the M W r s observations formed from the individual RINEX observations and the GPS and Galileo biases μ s to solve the N w l , r s for integers. Then the N w l , r s and satellite orbit and clock products are used as inputs to form the ionosphere-free combination measurements (code and carrier phase) for the first processing. The system of equations is solved for the PPP solution. In the case of PPP-AR processing the proceedure continues further after having fixed the N r , i s to integers; in this second processing to give the PPP-AR solution mostly ionosphere-free ambiguity-free carrier phase measurements are used for GPS and/or Galileo.
For the experimentation, one week of data (11–17 February 2019, Day of Year (DOY): 42–49/2019) is chosen. During that period, 31 GPS satellites and 24 Galileo satellites (including the ecliptic E14 and E18) were processed. All models and processing parameters and given in Table 2.
The IGS compiles a consistent set of absolute Antenna Phase Centre (APC) (i.e., Phase Center Offset (PCO) and Phase Center Variation (PCV)) corrections for both ground stations and satellites antennas, which are provided in so-called the Antenna Exchange Format (ANTEX) files [20]. These are very important for the calculation of the geometrical distance between satellite and receiver. For the APC corrections, it was decided to use the ANTEX14 file from the IGS. During the processing of the present article, ANTEX14 file included the APC values of the receivers for GPS L1/L2 frequencies, but not of Galileo E1/E5a frequencies. Delivery of the respective receiver APC values for the Galileo E1/E5a frequencies was underway from the IGS Antenna WG [33].
For the ambiguity fixing step biases, we do not use any additional bias for Galileo ambiguity fixing since a previous study [6] showed that the Galileo μ s biases were compatible with all kind of receivers and modulations (such as L1C, L1W, L1X, etc.).

3. Results

The following section shows some examples of PPP and PPP-AR positioning in detail as well as global network summary graphs in East (E), North (N) and Up (U) components.

3.1. Some Station Examples

The station BRUX from the IGS network has been chosen to show in detail the temporal series for the entire week, as examples of good PPP and PPP-AR repeatability. Figure 2 gives examples for PPP and PPP-AR solutions for Galileo-only, GPS-only and their combinations. For each temporal series the mean value is calculated and subtracted to center the series to zero. In addition, Table 3 shows the percentages of ambiguity resolution for the PPP-AR mode.
From Figure 2, it is observed that certain irregularities of the repeatability patterns in the PPP solutions do almost disappear in a PPP-AR solution. There is a little jump at the end of DOY 043 for north and east directions coming from the contribution of the GPS system (also seen in GPS-only PPP of around 2 cm in East and 1 cm in North) to the Multi-GNSS solution. Nevertheless, it is seen that in Multi-GNSS PPP-AR this jump is reduced to less than 0.5 cm.
The importance of AR in precise positioning is also seen when comparing the PPP solution to the PPP-AR solution. It is seen that some jumps in the PPP solution no longer appear when performing PPP-AR processing. The PPP-AR mode is shown to be smoother and more linear than the PPP mode (in particular for the east direction). Another example is the one seen for Galileo-only around DOY 046. There is a downward jump of around 2 cm for east and north directions that is eliminated for the PPP-AR mode.
The 1-σ values of the above temporal series for BRUX and another two examples of stations (CAS1 and NYA) for PPP and PPP-AR are gathered in the following Table 4, Table 5 and Table 6.
From these examples it is clear that level of accuracy achieved varies for each individual station. This could be because of several parameters such as multipath, station ANTEX parameters (i.e., for Galileo station, antenna ANTEX files are not yet provided so therefore the GPS station antenna ANTEX files were used), etc. Notwithstanding, it is seen that the level of accuracy from Galileo-only solutions is nearly comparable to the one of GPS-only solutions. For the Galileo-only solutions there were less measurements used than for the GPS-only solution due to the fact that the Galileo constellation has less satellites than GPS. It is expected that once the Galileo constellation is complete, the accuracy of PPP and PPP-AR respectively will improve. Ambiguity resolution improves the solution about 1–2 mm in East and North directions (around 10–45% improvement) (Here and elsewhere in the text, improvements in % are computed according to the formula: 100 % · ( v 2 v 1 ) / v 1 , where v 1 and v 2 are the values of before (hence reference) and after the change.) and about 1–2 mm in up direction (around 5–20% improvement). It is observed that the up component is improving less than east and north components when comparing the PPP and the PPP-AR cases using a single GNSS system. This is explained due to the fact that highly correlated parameters: i.e., up component, the tropospheric parameters (i.e., Zenith Tropospheric Delay (ZTD)) and station clocks are better de-correlated in the Multi-GNSS solutions. This happens mainly because more satellite measurements are participating in the LSE solutions.
It is interesting to observe that even though the Galileo-only solution is not better than the GPS-only solution (but it is of the same order), when using both systems the combined solution is improved with respect to the GPS only one. This means that adding Galileo ameliorates the overall performance of positioning both in PPP as well as in PPP-AR mode.

3.2. Global Network of Stations

The previous graphs from the example stations give some indications:
  • The PPP-AR mode gives better repeatability for the timeseries than the PPP mode.
  • Galileo only solution gives similar level of 1-σ values repeatability than the GPS only solution.
  • The use of Galileo can improve the current precise positioning situation of GPS when used in a Multi-GNSS combination in a PPP-AR mode.
Nevertheless, the graphs shown in detail referred to a very small specimen of stations. It was considered important to process a network of 50 the IGS stations used in MGEX. In this way it can be investigated whether these assumptions are valid globally and whether there are any potential geographical dependencies that affect the positioning accuracy.
Figure 3 and Figure 4 show the overall network performance of the 1-σ values (as computed for Figure 2) as well and the global RMS values (the RMS of all the values plotted) for GPS PPP-AR and GPS+Galileo PPP-AR scenarios for east, north and up components.
The following Table 7 is giving the RMS of all AR percentages for the entire network and the RMS for the whole week of processing. It is observed that in general the Galileo percentages are a little higher than the GPS ones.
The global positioning RMS values are gathered in Table 8 for the three directions from the global maps. As it is seen, adding Galileo to the constellation can improve the positioning globally both for PPP and PPP-AR modes. Even the PPP mode of Multi-GNSS gives better accuracy than the GPS PPP-AR (which is considered as the best positioning that can be achieved until now). Ambiguity resolution improves the solution about 1–3 mm in East and North directions (around 10–20% improvement) and about 0–2 mm in Up direction (around 2–8% improvement). Combining Galileo with GPS can improve the timeseries by around 3–4 mm in East and North directions and by 9 mm in Up direction. Comparing the current best positioning performance (i.e., GPS-only PPP-AR) with the one when adding Galileo (i.e., GPS+Galileo PPP-AR), we see that the results improve from 9.3 mm to 6.7 mm (28%) for the East component, 8.3 mm to 5.6 mm (33%) for the North component and 24 mm to 16.8 mm (30%) for the Up component.

4. Discussion and Conclusions

In this paper we investigated the performance (in term of repeatability) of kinematic PPP and PPP-AR positioning using GPS, Galileo and combined measurements of the two systems on a global network.
The Galileo-only PPP-AR solution performance repeatability reaches 14 mm in horizontal direction and 31 mm in vertical direction. This level is of similar order of magnitude and just below the GPS-only solution one. It is expected that once the Galileo constellation is complete, the accuracy of Galileo-only PPP and PPP-AR will improve relative to today (due to the addition of four more satellites and to the availability of receiver antenna calibration for Galileo L5 frequency).
The Multi-GNSS solutions for both PPP and PPP-AR give much better results in term of repeatability than both systems used separately (we observe a gain in repeatability of around 30% in horizontal and in vertical directions); and this even though the Galileo-only solutions are by little not better than the GPS-only solutions. This is logically explained by the increased number of measurements and the highest satellite geometry diversity in the Multi-GNSS solutions relatively to the single systems ones (with only one inter-system bias parameter added). This result proves also that despite the unavailable receiver’s antenna patterns for the E5 Galileo frequencies used in this study, the two systems are already compatible at the sub-centimeter level. We assume, as a consequence, that parameters that are highly correlated such as the vertical components, the tropospheric ZTD parameter and the station clocks are more de-correlated in the Multi-GNSS solutions than in single-systems ones.
The best repeatabilities observed in this study for a 300 s dynamic Multi-GNSS GPS + GAL PPP-AR solution reach less than 3 mm in horizontal direction and 7 mm in the vertical one for the BRUX receiver. The global network statistics are a little larger with 8 mm in horizontal and 17 mm in vertical directions. For the global network, no particular geographical pattern has been seen among the 50 stations used.
These results show that Galileo can really contribute to Multi-GNSS precise positioning and improve the best solutions obtained today with GPS. Scientific, geodetic and geophysical applications can already benefit from the combined processing of GPS and Galileo observations.

Author Contributions

Conceptualization, G.K., S.L. and F.P.; methodology, G.K.; software, S.L.; G.K. writing—original draft preparation, S.L. and F.P.; writing—review and editing.

Funding

This research received no external funding.

Acknowledgments

All authors would like to acknowledge the work of Mini Gupta and Jean-Charles Marty for their work changing the EXE-PPP and the GINS software. These were extensively used for the experiments of this article. In addition, fruitful discussions with Flavien Mercier and Alvaro Santamaria were very much appreciated.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zumberge, J.F.; Heflin, M.B.; Jefferson, D.C.; Watkins, M.M.; Webb, F.H. Precise point positioning for the efficient and robust analysis of GPS data from large networks. J. Geophys. Res. 1997, 102, 5005–5017. [Google Scholar] [CrossRef] [Green Version]
  2. Lescarmontier, L.; Legrésy, B.; Coleman, R.; Perosanz, F.; Mayet, C.; Testut, L. Vibrations of Mertz Glacier ice tongue, East Antarctica. J. Glaciol. 2012, 58, 665–676. [Google Scholar] [CrossRef]
  3. Fund, F.; Perosanz, F.; Testut, L.; Loyer, S. An Integer Precise Point Positioning technique for sea surface observations using a GPS buoy. Adv. Space Res. 2013, 51, 1311–1322. [Google Scholar] [CrossRef]
  4. Petit, G.; Kanj, A.; Loyer, S.; Delporte, J.; Mercier, F.; Perosanz, F. 1 × 10−16 frequency transfer by GPS PPP with integer ambiguity resolution. Metrologia 2015, 52, 301–309. [Google Scholar] [CrossRef]
  5. Laurichesse, D.; Mercier, F.; Berthias, J.; Broca, P.; Cerri, L. Integer ambiguity resolution on undifferenced GPS phase measurements and its application to PPP and satellite precise orbit determination. Navig. J. Inst. Navig. 2009, 2, 135–149. [Google Scholar] [CrossRef]
  6. Katsigianni, G.; Loyer, S.; Perosanz, F.; Mercier, F.; Zajdel, R.; Sośnica, K. Improving Galileo orbit determination using zero-differenceambiguity fixing in a Multi-GNSS processing. Adv. Space Res. 2019, 63, 2952–2963. [Google Scholar] [CrossRef]
  7. Montenbruck, O.; Steigenberger, P.; Prange, L.; Deng, Z.; Zhao, Q.; Perosanz, F.; Romero, I.; Noll, C.; Stürze, A.; Weber, G.; et al. The Multi-GNSS Experiment (MGEX) of the International GNSS Service (IGS)—Achievements, prospects and challenges. Adv. Space Res. 2017, 59, 1671–1697. [Google Scholar] [CrossRef]
  8. Loyer, S.; Perosanz, F.; Mercier, F.; Capdeville, H.; Marty, J. Zero-difference GPS ambiguity resolution at CNES-CLS IGS analysis center. J. Geod. 2012, 991–1003. [Google Scholar] [CrossRef]
  9. Hadas, T.; Kazmierski, K.; Sośnica, K. Performance of Galileo-only dual-frequency absolute positioning using the fully serviceable Galileo constellation. GPS Solut. 2019, 23, 12. [Google Scholar] [CrossRef]
  10. Katsigianni, G.; Perosanz, F.; Loyer, S.; Gupta, M. Galileo millimeter-level kinematic precise point positioning with ambiguity resolution. Earth Planets Space 2019, 71–76. [Google Scholar] [CrossRef]
  11. GSA, Constellation Information. Available online: https://www.gsc-europa.eu/system-service-status/constellation-information (accessed on 4 October 2019).
  12. ESA, Four New Galileos Join Europe’s Largest Satellite Constellation. Available online: https://www.esa.int/Our_Activities/Navigation/Four_new_Galileos_join_Europe_s_largest_satellite_constellation (accessed on 4 October 2019).
  13. Laurichesse, D.; Mercier, F. Integer ambiguity resolution on undifferenced GPS phase measurements and its application to PPP. In Proceedings of the ION GNSS 2007 20th International Technical Meeting of the Satellite Division, Fort Worth, TX, USA, 25–28 September 2007. [Google Scholar]
  14. Melbourne, W. The case for ranging in GPS based geodetic system. In Proceedings of the 1st International Symposium on Precise Positioning With the Global Positioning System. U.S. Department of Commerce, Rockville, MD, USA, 15–19 April 1985. [Google Scholar]
  15. Wübbena, G. Software developments for geodetic positioning with GPS using TI-4100 code and carrier measurements. In Proceedings of the 1st International Symposium on Precise Positioning With the Global Positioning System. U.S. Department of Commerce, Rockville, MD, USA, 15–19 April 1985. [Google Scholar]
  16. Mercier, F.; Laurichesse, D. Receiver/Payload hardware biases stability requirements for undifferenced Widelane ambiguity blocking. In Proceedings of the Scientific and Fundamental Aspects of the Galileo Program, Toulouse, France, 1–4 October 2007. [Google Scholar]
  17. Verhagen, S. The GNSS Integer Ambiguities: Estimation and Validation; Publications on Geodesy: Delft, The Netherlands, 2004; pp. 27–32. [Google Scholar]
  18. Teunissen, P. Success probability of integer GPS ambiguity rounding and bootstrapping. J. Geod. 1998, 72, 606–612. [Google Scholar] [CrossRef] [Green Version]
  19. Hofmann-Wellenhof, B.; Lichtenegger, H.; Wasle, E. GNSS-Global Navigation Satellite Systems GPS, GLONASS, Galileo and More; Springer: Vienna, Austria, 2008; p. 59. [Google Scholar]
  20. Kouba, J. A Guide to Using International GNSS Service (IGS) Products. Available online: http://acc.igs.org/UsingIGSProductsVer21.pdf (accessed on 4 October 2019).
  21. Marty, J.C. Algorithmic Documentation of the GINS Software. Available online: https://www5.obs-mip.fr/wp-content-omp/uploads/sites/28/2017/11/GINS_Algo_2013.pdf (accessed on 4 October 2019).
  22. CNES/CLS Analysis Center for IGS. Available online: https://igsac-cnes.cls.fr/ (accessed on 4 October 2019).
  23. GSA, Galileo Satellite Metadata. Available online: https://www.gsc-europa.eu/support-to-developers/galileo-satellite-metadata (accessed on 4 October 2019).
  24. Böhm, J.; Möller, G.; Schindelegger, M.; Pain, G.; Weber, R. Development of an improved empirical model for slant delays in the troposphere (GPT2w). GPS Solut. 2015, 19, 433–441. [Google Scholar] [CrossRef]
  25. Lagler, K.; Schindelegger, M.; Böhm, J.; Krásná, H.; Nilsson, T. GPT2: Empirical slant delay model for radio space geodetic techniques. Geophys. Res. Lett. 2013, 40, 1069–1073. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Hernández-Pajares, M.; Juan, J.M.; Sanz, J.; Orús, R. Second-order ionospheric term in GPS: Implementation and impact on geodetic estimates. J. Geophys. Res. 2007, 112. [Google Scholar] [CrossRef]
  27. Altamimi, Z.; Rebischung, P.; Métivier, L.; Collilieux, X. ITRF2014: A new release of the International Terrestrial Reference Frame modeling nonlinear station motions. J. Geophys. Res. Solid Earth 2016, 12, 6109–6131. [Google Scholar] [CrossRef]
  28. Kouba, J. A simplified yaw-attitude model for eclipsing GPS satellites. GPS Solut. 2009, 13, 1–12. [Google Scholar] [CrossRef]
  29. Carrère, L.; Lyard, M.; Cancet, M.; Guillot, A.; Roblou, L. FES2012: A new global tidal model taking advantage of nearly 20 years of altimetry. In Proceedings of the 20 Years of Progress in Radar Altimetry Symposium, Venice, Italy, 24–29 September 2012. [Google Scholar]
  30. Petit, G.; Luzum, B. IERS Conventions; Verlag des Bundesamts für Kartographie und Geodäsie: Frankfurt am Main, Germany, 2010; Volume 36, p. 179. [Google Scholar]
  31. Bizouard, C.; Lambert, S.; Gattano, C.; Becker, O.; Richard, J.-Y. The IERS EOP 14C04 solution for Earth orientation parameters consistent with ITRF 2014. J. Geod. 2019, 93, 621–633. [Google Scholar] [CrossRef]
  32. Chen, G.; Herring, T.A. Effects of atmospheric azimuthal asymmetry on the analysis of space geodetic data. J. Geophys. Res. 1997, 102, 20489–20502. [Google Scholar] [CrossRef]
  33. IGS, Antenna Working Group Charter and Members. Available online: https://kb.igs.org/hc/en-us/articles/202023633-Antenna-Working-Group-Charter-and-Members (accessed on 4 October 2019).
Figure 1. Steps of the procedure to perform precise point positioning (PPP) and PPP with phase ambiguity resolution (PPP-AR) of a combined Multi-GNSS solution (GPS + Galileo).
Figure 1. Steps of the procedure to perform precise point positioning (PPP) and PPP with phase ambiguity resolution (PPP-AR) of a combined Multi-GNSS solution (GPS + Galileo).
Remotesensing 11 02477 g001
Figure 2. PPP and PPP-AR solutions for BRUX station for Galileo-only, GPS-only and Multi-GNSS for the period 11–19 February 2019, DOY: 42–49/2019. For Galileo-only and GPS + Galileo, additional biases of +0.05 m and −0.05 m have been added respectively for better representation in the graph.
Figure 2. PPP and PPP-AR solutions for BRUX station for Galileo-only, GPS-only and Multi-GNSS for the period 11–19 February 2019, DOY: 42–49/2019. For Galileo-only and GPS + Galileo, additional biases of +0.05 m and −0.05 m have been added respectively for better representation in the graph.
Remotesensing 11 02477 g002
Figure 3. GPS PPP-AR solutions for the network of IGS stations in East, North and Up components and their global RMS (mm).
Figure 3. GPS PPP-AR solutions for the network of IGS stations in East, North and Up components and their global RMS (mm).
Remotesensing 11 02477 g003
Figure 4. GPS + Galileo PPP-AR solutions for the network of IGS stations in East, North and Up components and their global RMS (mm).
Figure 4. GPS + Galileo PPP-AR solutions for the network of IGS stations in East, North and Up components and their global RMS (mm).
Remotesensing 11 02477 g004
Table 1. Values for wide-lane and narrow-lane wavelength.
Table 1. Values for wide-lane and narrow-lane wavelength.
GNSS Frequency λ w l ( m ) λ n l ( m )
GPS ( L 1 , L 2 )0.8620.107
Galileo ( E 1 , E 5 a )0.7510.109
Table 2. Parameters, models and strategy of experiments.
Table 2. Parameters, models and strategy of experiments.
Processing Strategy
SoftwareGINS, DYNAMO, EXE_PPP [21]
StrategyPPP, PPP-AR zero-difference
EstimationStatic with Kalman, 300 s sampling
Orbit, Clocks and Satellites Biaises
Orbits and clocksCNES-CLS orbits (‘grg’) [8]
Satellite biasesCNES-CLS wide-lane satellite biases [22]
GNSS relative weightingEqual weighting for GPS and Galileo
GNSS measurement sigmas
(at 0° of elevation)
Code: 60 cm; Phase: 3.5 mm
Elevation cut-off8 deg
Elevation weighting function, where ϑ is the elevation angle σ ( ϑ ) = 0.0035 0.15 + 0.85 sin ϑ
Models for Processing
Antenna phase center correctionsANTEX14 PCO/PCV [23]
Troposphere modelVMF1 [24] + GPT2 [25] (A priori local meteorological parameters (pressure, temperature, and wet mapping function coefficients) of GPT2 model are used to compute hydrostatic delays and for the wet mapping function VMF1. We then adjust 1 zenithal tropospheric delay per two hours in factor of the wet mapping function).
Ionosphere modelIonosphere-free combination and second ordercorrections [26]
Reference frameITRF 2014 [27]
Attitude modelKouba [28] for GPS and GSA [23] for Galileo
Ocean loading effectsFES2012 [29]
Earth orientation modellingIERS Conventions 2010 [30]
Earth orientation parametersEOP C04 [31]
Phase windupModels used by Kouba [28]
Estimated Parameters
Troposphere1 ZTD/2 h, 1 pair of gradients (E, N)/day (1 couple of gradients in north and east direction are also adjusted per day following [31,32]
Observation sampling300 s
Inter-system biases (phase obs.)1 per station (zero mean condition)
Station coordinates estimatesX, Y, Z transformed to East, North, Up
Table 3. Ambiguity resolution percentages for Galileo and GPS system for BRUX station.
Table 3. Ambiguity resolution percentages for Galileo and GPS system for BRUX station.
BRUX AR (%)042043044045046047048
Galileo95.2410010097.6710010095.45
GPS98.6391.6791.3098.4410098.4894.03
Table 4. 1-σ values PPP and PPP-AR for BRUX station.
Table 4. 1-σ values PPP and PPP-AR for BRUX station.
BRUXModeEast (mm)North (mm)Up (mm)
GalileoPPP4.74.611.7
PPP-AR2.62.910.3
GPSPPP4.74.19.2
PPP-AR2.43.48.5
GPS + GalileoPPP3.42.79.1
PPP-AR2.12.47.3
Table 5. 1-σ values PPP and PPP-AR for CAS1 station.
Table 5. 1-σ values PPP and PPP-AR for CAS1 station.
CAS1ModeEast (mm)North (mm)Up (mm)
GalileoPPP7.16.816.6
PPP-AR4.25.215.6
GPSPPP6.46.815.2
PPP-AR3.85.214.3
GPS + GalileoPPP4.64.511.3
PPP-AR3.13.710.3
Table 6. 1-σ values PPP and PPP-AR for NYA2 station.
Table 6. 1-σ values PPP and PPP-AR for NYA2 station.
NYA2ModeEast (mm)North (mm)Up (mm)
GalileoPPP4.95.013.9
PPP-AR2.82.915.5
GPSPPP4.24.116.0
PPP-AR2.42.211.7
GPS + GalileoPPP3.33.110.3
PPP-AR2.52.19.8
Table 7. Ambiguity resolution percentages for Galileo and GPS system for BRUX station.
Table 7. Ambiguity resolution percentages for Galileo and GPS system for BRUX station.
Network AR (%)042043044045046047048RMS
Galileo90.7195.3895.3995.5696.1695.6789.3994.07
GPS93.5592.2589.8593.9094.0994.5489.6292.56
Table 8. Global RMS of 1-σ values for PPP and PPP-AR for the entire network examined.
Table 8. Global RMS of 1-σ values for PPP and PPP-AR for the entire network examined.
GlobalModeEast (mm)North (mm)Up (mm)
GalileoPPP17.014.633.1
PPP-AR13.712.230.8
GPSPPP11.89.426.0
PPP-AR9.38.324.0
GPS + GalileoPPP7.96.117.2
PPP-AR6.75.616.8

Share and Cite

MDPI and ACS Style

Katsigianni, G.; Loyer, S.; Perosanz, F. PPP and PPP-AR Kinematic Post-Processed Performance of GPS-Only, Galileo-Only and Multi-GNSS. Remote Sens. 2019, 11, 2477. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11212477

AMA Style

Katsigianni G, Loyer S, Perosanz F. PPP and PPP-AR Kinematic Post-Processed Performance of GPS-Only, Galileo-Only and Multi-GNSS. Remote Sensing. 2019; 11(21):2477. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11212477

Chicago/Turabian Style

Katsigianni, Georgia, Sylvain Loyer, and Felix Perosanz. 2019. "PPP and PPP-AR Kinematic Post-Processed Performance of GPS-Only, Galileo-Only and Multi-GNSS" Remote Sensing 11, no. 21: 2477. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11212477

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop