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

Georgia Katsigianni, Sylvain Loyer, Felix Perosanz

► To cite this version:

Georgia Katsigianni, Sylvain Loyer, Felix Perosanz. PPP and PPP-AR Kinematic Post-Processed Performance of GPS-Only, Galileo-Only and Multi-GNSS. Remote Sensing, 2019, 11 (21), pp.2477. 10.3390/rs11212477. [hal-04903318]

HAL Id: hal-04903318

https://ut3-toulouseinp.hal.science/hal-04903318v1

Submitted on 21 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

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

Georgia Katsigianni<fnref n="1" />,<fnref n="2" />,* , Sylvain Loyer<fnref n="2" /> and Felix Perosanz<fnref n="1" />

Received: 31 August 2019; Accepted: 21 October 2019; Published: 24 October 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**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 around 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.

**Keywords:** GNSS; GPS; Galileo; precise point positioning (PPP); PPP with phase ambiguity resolution (PPP-AR); integer PPP (IPPP); ambiguity resolution; zero-difference ambiguity fixing; multi-GNSS

## 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

<footnote n="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; felix.perosanz@cnes.fr</footnote>

<footnote n="2">² Collecte Localisation Satellites (CLS Group), 8-10 rue Hermès, 31520 Ramonville Saint Agne, France; sloyer@groupcls.com</footnote>

<footnote>* Correspondence: georgia.katsigianni@cnes.fr; Tel.: +33-(0)5-6133-2892</footnote>
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 below is given for the pseudorange (code) and carrier phase measurements for two frequencies:

$$P_{r,i}^{s} = \rho_{r}^{s} + c\Delta t + T_{r}^{s} + I_{r}^{s} + b_{i}^{s} + b_{r,i} + E_{r,i}^{s} \quad (1)$$

$$P_{r,j}^{s} = \rho_{r}^{s} + c\Delta t + T_{r}^{s} + \frac{f_{i}^{2}}{f_{j}^{2}} I_{r}^{s} + b_{j}^{s} + b_{r,j} + E_{r,j}^{s} \quad (2)$$

$$L_{r,i}^{s} = \lambda_{i}\varphi_{r,i}^{s} = \rho_{r}^{s} + c\Delta t + T_{r}^{s} - I_{r}^{s} + \lambda_{i}N_{r,i}^{s} + \lambda_{i}W_{r}^{s} + \beta_{i}^{s} + \beta_{r,i} + \varepsilon_{r,i}^{s} \quad (3)$$

$$L_{r,j}^{s} = \lambda_{j}\varphi_{r,j}^{s} = \rho_{r}^{s} + c\Delta t + T_{r}^{s} - \frac{f_{i}^{2}}{f_{j}^{2}} I_{r}^{s} + \lambda_{j}N_{r,j}^{s} + \lambda_{j}W_{r}^{s} + \beta_{j}^{s} + \beta_{r,j} + \varepsilon_{r,j}^{s} \quad (4)$$
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)

* $\rho_r^s$ is the geometric distance between receiver and satellite (m)

* $\Delta t$ ($\Delta t = \delta t_r - \delta t^s$) is the clock correction related to the satellite ($\delta t^s$) and the receiver ($\delta 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)

* $\lambda_i, \lambda_j$ are the nominal wavelength of the carrier frequency $i$ or $j$ (m)

* $\varphi_{r,i'}^s, \varphi_{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)

* $\beta^s, \beta_r$ are the carrier phase biases of satellite and receiver (m)

* $\epsilon_{r,i'}^s, \epsilon_{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 L1 and L2 ($f_{L1} = 154 \times f_0$ and $f_{L2} = 120 \times f_0$). For the Galileo system are in the band of E1 and E5a ($f_{E1} = 154 \times f_0$ and $f_{E5a} = 115 \times 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]:

$$ 
\begin{aligned}
MW_r^s &= \lambda_{wl} \varphi_r^s = \left( \frac{f_i}{f_i - f_j} L_{r,i}^s - \frac{f_j}{f_i - f_j} L_{r,j}^s \right) - \left( \frac{f_i}{f_i + f_j} P_{r,i}^s + \frac{f_j}{f_i + f_j} P_{r,j}^s \right) \\
&= \lambda_{wl} (N_{r,i}^s - N_{r,j}^s - \mu^s + \mu_r(t)) = \lambda_{wl} (N_{wl,r}^s - \mu^s + \mu_r(t))
\end{aligned}
\quad (5) $$

where:

* $MW_r^s$ is the Melbourne-Wübbena linear combination at receiver $r$ from satellite $s$ (m)

* $\lambda_{wl}$ ($\lambda_{wl} = c/(f_i - f_j) = \lambda_i \lambda_j / (\lambda_j - \lambda_i)$) is the wide-lane (WL) wavelength (m)

* $N_{wl,r}^s$ ($N_{wl,r}^s = N_{r,i}^s - N_{r,j}^s$) is the WL ambiguity at receiver $r$ from satellite $s$

* $\mu^s$ is the delay coming from the satellite (also known in the bibliography as WL satellite bias (WSB))

* $\mu_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 $\mu^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 $\mu_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_{wl,r}^s, \mu^s$ and $\mu_r(t)$ are totally correlated. Normally, the $\mu^s$ is known; i.e., it is calculated and provided from the CNES/CLS AC to the users [6,8]. The $\mu_r$ (one per epoch of measurements) and the $N_{wl,r}^s$ (one per satellite pass) can be separated and solved
using all available equations (Equation (5) corrected by the $\mu^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_{wl,r}^s$ and $\mu_r(t)$.

Once the $N_{wl,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: $\alpha_i$ to the frequency i and $\alpha_j$ to the frequency j respectively [19].

$$ \alpha_i = \frac{f_i^2}{f_i^2 - f_j^2} \qquad (6) $$

$$ \alpha_j = \frac{-f_i^2}{f_i^2 - f_j^2} \qquad (7) $$

Using the above coefficients and Equations (6) and (7), the ionosphere-free linear combinations for code and carrier phase become:

$$ P_{r,IF}^s = \alpha_i P_{r,i}^s + \alpha_j P_{r,j}^s = \rho_r^s + c\Delta t + T_r^s + E_{r,IF}^s \qquad (8) $$

$$ \begin{aligned} L_{r,IF}^s &= \alpha_i L_{r,i}^s + \alpha_j L_{r,j}^s \\ &= \rho_r^s + c\Delta t + T_r^s + \frac{c}{f_i^2 - f_j^2} (f_i N_{r,i}^s - f_j N_{r,j}^s) + \frac{c}{f_i^2 - f_j^2} (f_i W_r^s - f_j W_r^s) + \varepsilon_{r,IF}^s \qquad (9) \\ &= \rho_r^s + c\Delta t + T_r^s + \lambda_{nl} N_{r,i}^s + \frac{\lambda_{nl} \lambda_{wl}}{\lambda_j} N_{wl,r}^s + \lambda_{nl} W_r^s + \varepsilon_{r,IF}^s \end{aligned} $$

where:

* $P_{r,IF}^s$ is the ionosphere-free code measurement at receiver r from satellite s (m)

* $E_{r,IF}^s$ is the ionosphere-free code measurement error at receiver r from satellite s (m)

* $L_{r,IF}^s$ is the ionosphere-free carrier measurement at receiver r from satellite s (m)

* $\varepsilon_{r,IF}^s$ is the ionosphere-free carrier measurement error at receiver r from satellite s (m)

* $\lambda_{nl}$ ($\lambda_{nl} = c/(f_i + f_j) = \lambda_i \lambda_j / (\lambda_i + \lambda_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.

<table><caption>**Table 1.** Values for wide-lane and narrow-lane wavelength.</caption><thead><tr><td>GNSS Frequency</td><td>λ<sub>wl</sub>(m)</td><td>λ<sub>nl</sub>(m)</td></tr></thead><tbody><tr><td>GPS (L1, L2)</td><td>0.862</td><td>0.107</td></tr><tr><td>Galileo (E1, E5a)</td><td>0.751</td><td>0.109</td></tr></tbody></table>

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.
<table><caption>**Table 2.** Parameters, models and strategy of experiments.</caption><thead><tr><th colspan="2">Processing Strategy</th></tr></thead><tbody><tr><td>Software</td><td>GINS, DYNAMO, EXE_PPP [21]</td></tr><tr><td>Strategy</td><td>PPP, PPP-AR zero-difference</td></tr><tr><td>Estimation</td><td>Static with Kalman, 300 s sampling</td></tr><tr><th colspan="2">Orbit, Clocks and Satellites Biases</th></tr><tr><td>Orbits and clocks</td><td>CNES-CLS orbits ('grg') [8]</td></tr><tr><td>Satellite biases</td><td>CNES-CLS wide-lane satellite biases [22]</td></tr><tr><td>GNSS relative weighting</td><td>Equal weighting for GPS and Galileo</td></tr><tr><td>GNSS measurement sigmas (at 0° of elevation)</td><td>Code: 60 cm; Phase: 3.5 mm</td></tr><tr><td>Elevation cut-off</td><td>8 deg</td></tr><tr><td>Elevation weighting function, where θ is the elevation angle</td><td>σ(θ) = <sup>0.0035</sup>/<sub>0.15+0.85 sin θ</sub></td></tr><tr><th colspan="2">Models for Processing</th></tr><tr><td>Antenna phase center corrections</td><td>ANTEX14 PCO/PCV [23]</td></tr><tr><td>Troposphere model</td><td>VMF1 [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 zenith tropospheric delay per two hours in factor of the wet mapping function).</td></tr><tr><td>Ionosphere model</td><td>Ionosphere-free combination and second ordercorrections [26]</td></tr><tr><td>Reference frame</td><td>ITRF 2014 [27]</td></tr><tr><td>Attitude model</td><td>Kouba [28] for GPS and GSA [23] for Galileo</td></tr><tr><td>Ocean loading effects</td><td>FES2012 [29]</td></tr><tr><td>Earth orientation modelling</td><td>IERS Conventions 2010 [30]</td></tr><tr><td>Earth orientation parameters</td><td>EOP C04 [31]</td></tr><tr><td>Phase windup</td><td>Models used by Kouba [28]</td></tr><tr><th colspan="2">Estimated Parameters</th></tr><tr><td>Troposphere</td><td>1 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]</td></tr><tr><td>Observation sampling</td><td>300 s</td></tr><tr><td>Inter-system biases (phase obs.)</td><td>1 per station (zero mean condition)</td></tr><tr><td>Station coordinates estimates</td><td>X, Y, Z transformed to East, North, Up</td></tr></tbody></table>

Figure 1 is giving a graphic overview of the PPP and PPP-AR method used. In the beginning, the satellite $\mu^s$ are needed together with the station RINEX file. The pre-processing phase uses the $MW_r^s$ observations formed from the individual RINEX observations and the GPS and Galileo biases $\mu^s$ to solve the $N_{wl,r}^s$ for integers. Then the $N_{wl,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 procedure 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.
<figure><img src="image_4.png" /><figcaption>**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).</figcaption></figure>

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 $\mu^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.

<figure><img src="image_5.png" /><figcaption>**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.</figcaption></figure>

<table><caption>**Table 3.** Ambiguity resolution percentages for Galileo and GPS system for BRUX station.</caption><thead><tr><th>BRUX AR (%)</th><th>042</th><th>043</th><th>044</th><th>045</th><th>046</th><th>047</th><th>048</th></tr></thead><tbody><tr><td>Galileo</td><td>95.24</td><td>100</td><td>100</td><td>97.67</td><td>100</td><td>100</td><td>95.45</td></tr><tr><td>GPS</td><td>98.63</td><td>91.67</td><td>91.30</td><td>98.44</td><td>100</td><td>98.48</td><td>94.03</td></tr></tbody></table>

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 Tables 4–6.
<table><caption>**Table 4.** 1-σ values PPP and PPP-AR for BRUX station.</caption><thead><tr><th>BRUX</th><th>Mode</th><th>East (mm)</th><th>North (mm)</th><th>Up (mm)</th></tr></thead><tbody><tr><td rowspan="2">Galileo</td><td>PPP</td><td>4.7</td><td>4.6</td><td>11.7</td></tr><tr><td>PPP-AR</td><td>2.6</td><td>2.9</td><td>10.3</td></tr><tr><td rowspan="2">GPS</td><td>PPP</td><td>4.7</td><td>4.1</td><td>9.2</td></tr><tr><td>PPP-AR</td><td>2.4</td><td>3.4</td><td>8.5</td></tr><tr><td rowspan="2">GPS + Galileo</td><td>PPP</td><td>3.4</td><td>2.7</td><td>9.1</td></tr><tr><td>PPP-AR</td><td>2.1</td><td>2.4</td><td>7.3</td></tr></tbody></table>

<table><caption>**Table 5.** 1-σ values PPP and PPP-AR for CAS1 station.</caption><thead><tr><th>CAS1</th><th>Mode</th><th>East (mm)</th><th>North (mm)</th><th>Up (mm)</th></tr></thead><tbody><tr><td rowspan="2">Galileo</td><td>PPP</td><td>7.1</td><td>6.8</td><td>16.6</td></tr><tr><td>PPP-AR</td><td>4.2</td><td>5.2</td><td>15.6</td></tr><tr><td rowspan="2">GPS</td><td>PPP</td><td>6.4</td><td>6.8</td><td>15.2</td></tr><tr><td>PPP-AR</td><td>3.8</td><td>5.2</td><td>14.3</td></tr><tr><td rowspan="2">GPS + Galileo</td><td>PPP</td><td>4.6</td><td>4.5</td><td>11.3</td></tr><tr><td>PPP-AR</td><td>3.1</td><td>3.7</td><td>10.3</td></tr></tbody></table>

<table><caption>**Table 6.** 1-σ values PPP and PPP-AR for NYA2 station.</caption><thead><tr><th>NYA2</th><th>Mode</th><th>East (mm)</th><th>North (mm)</th><th>Up (mm)</th></tr></thead><tbody><tr><td rowspan="2">Galileo</td><td>PPP</td><td>4.9</td><td>5.0</td><td>13.9</td></tr><tr><td>PPP-AR</td><td>2.8</td><td>2.9</td><td>15.5</td></tr><tr><td rowspan="2">GPS</td><td>PPP</td><td>4.2</td><td>4.1</td><td>16.0</td></tr><tr><td>PPP-AR</td><td>2.4</td><td>2.2</td><td>11.7</td></tr><tr><td rowspan="2">GPS + Galileo</td><td>PPP</td><td>3.3</td><td>3.1</td><td>10.3</td></tr><tr><td>PPP-AR</td><td>2.5</td><td>2.1</td><td>9.8</td></tr></tbody></table>

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.

Figures 3 and 4 show the overall network performance of the 1-σ values (as computed for Figure 2) as well as 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.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** GPS PPP-AR solutions for the network of IGS stations in East, North and Up components and their global RMS (mm).</figcaption></figure>
GPS+Galileo PPP-AR

<figure><img src="image_7.png" /><figcaption>**Figure 4.** GPS + Galileo PPP-AR solutions for the network of IGS stations in East, North and Up components and their global RMS (mm).</figcaption></figure>

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.

<table><caption>**Table 7.** Ambiguity resolution percentages for Galileo and GPS system for BRUX station.</caption><thead><tr><th>Network AR (%)</th><th>042</th><th>043</th><th>044</th><th>045</th><th>046</th><th>047</th><th>048</th><th>RMS</th></tr></thead><tbody><tr><td>Galileo</td><td>90.71</td><td>95.38</td><td>95.39</td><td>95.56</td><td>96.16</td><td>95.67</td><td>89.39</td><td>94.07</td></tr><tr><td>GPS</td><td>93.55</td><td>92.25</td><td>89.85</td><td>93.90</td><td>94.09</td><td>94.54</td><td>89.62</td><td>92.56</td></tr></tbody></table>

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.

<table><caption>**Table 8.** Global RMS of 1-σ values for PPP and PPP-AR for the entire network examined.</caption><thead><tr><th>Global</th><th>Mode</th><th>East (mm)</th><th>North (mm)</th><th>Up (mm)</th></tr></thead><tbody><tr><td rowspan="2">Galileo</td><td>PPP</td><td>17.0</td><td>14.6</td><td>33.1</td></tr><tr><td>PPP-AR</td><td>13.7</td><td>12.2</td><td>30.8</td></tr><tr><td rowspan="2">GPS</td><td>PPP</td><td>11.8</td><td>9.4</td><td>26.0</td></tr><tr><td>PPP-AR</td><td>9.3</td><td>8.3</td><td>24.0</td></tr><tr><td rowspan="2">GPS + Galileo</td><td>PPP</td><td>7.9</td><td>6.1</td><td>17.2</td></tr><tr><td>PPP-AR</td><td>6.7</td><td>5.6</td><td>16.8</td></tr></tbody></table>

## 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. [CrossRef]

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. [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. [CrossRef]

4. Petit, G.; Kanj, A.; Loyer, S.; Delporte, J.; Mercier, F.; Perosanz, F. $1 \times 10^{-16}$ frequency transfer by GPS PPP with integer ambiguity resolution. *Metrologia* **2015**, *52*, 301–309. [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. [CrossRef]

6. Katsigianni, G.; Loyer, S.; Perosanz, F.; Mercier, F.; Zajdel, R.; Sośnica, K. Improving Galileo orbit determination using zero-difference ambiguity fixing in a Multi-GNSS processing. *Adv. Space Res.* **2019**, *63*, 2952–2963. [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. [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. [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. [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. [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.

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.

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.

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.

17. Verhagen, S. *The GNSS Integer Ambiguities: Estimation and Validation*; Publications on Geodesy: Delft, The Netherlands, 2004; pp. 27–32.

18. Teunissen, P. Success probability of integer GPS ambiguity rounding and bootstrapping. *J. Geod.* **1998**, *72*, 606–612. [CrossRef]
19. Hofmann-Wellenhof, B.; Lichtenegger, H.; Wasle, E. GNSS-Global Navigation Satellite Systems GPS, GLONASS, Galileo and More; Springer: Vienna, Austria, 2008; p. 59.

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. [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. [CrossRef] [PubMed]

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*. [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. [CrossRef]

28. Kouba, J. A simplified yaw-attitude model for eclipsing GPS satellites. *GPS Solut.* **2009**, *13*, 1–12. [CrossRef]

29. Carrère, L.; Lyard, M.; Cancet, M.; Guillot, A.; Roblon, 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.

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.

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. [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. [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><img src="image_8.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).