Analysis of Wide-Lane Ambiguities Derived from Geometry-Free and Geometry-Based Precise Point Positioning Models and Their Implication for Orbit and Clock Quality

Orbit and clock products are used in real-time global navigation satellite systems (GNSS) precise point positioning (PPP) without knowing their quality. This study develops a new approach to detect orbit and clock errors through comparing geometry-free and geometry-based wide-lane ambiguities in a PPP model. The reparameterization and estimation procedures of the geometry-free and geometry-based ambiguities are described in detail. The effects of orbit and clock errors on ambiguities are given in analytical expressions. The numerical similarity and differences of geometry-free and geometry-based wide-lane ambiguities are analyzed using different orbit and clock products. Furthermore, two types of typical errors in orbit and clock are simulated and their effects on wide-lane ambiguities are numerically produced and analyzed. The contribution discloses that the geometry-free and geometry-based wide-lane ambiguities are equivalent in terms of their formal errors. Although they are very close in terms of their estimates when the used orbit and clock for geometry-based ambiguities are precise enough, they are not the same, in particular, in the case that the used orbit and clock, as a combination, contain significant errors. It is discovered that the discrepancies of geometry-free and geometry-based wide-lane ambiguities coincide with the actual time-variant errors in the used orbit and clock at the line-of-sight direction. This provides a quality index for real-time users to detect the errors in real-time orbit and clock products, which potentially improves the accuracy of positioning.


Introduction
Precise point positioning (PPP) using global navigation satellite systems (GNSS) measurements was initially proposed in [1] in order to solve the high computation burden of the network approach. The PPP technique only employs a single receiver to determine its position without explicitly using reference stations. To do so, precise orbit and clock are needed to enable absolute positioning of high accuracy [2]. This technique has been intensively developed from static to kinematic positioning [3], from post-processing to real-time service [4,5], from using global positioning system (GPS) only to using multi-constellation systems [6], from ambiguity-float to integer ambiguity-fixed solution [7][8][9][10], from decimeter to centimeter accuracy [11]. The PPP technique started with the ionosphere-free combination approach [1] and then developed with the raw observation approach [11][12][13][14][15]. One of the major differences of these two approaches is that one eliminates the ionosphere delay while the other the prediction length, the variation of errors could range from a couple of centimeters to decimeters. There are also exceptional occasions where some satellites, like GPS IIF, could have large errors up to meters in their predicted orbits [22]. These errors could significantly worsen the clock bias estimation and satellite phase bias estimation as well as PPP integer ambiguity resolution in user side, in particular for real-time users.
It is common to evaluate the quality of real-time orbit and clock by comparing them with the post-processing counterpart separately [21]. However, the post-processing products are available only after one day or a couple of hours later, which is not very acceptable for real-time users. In addition, it is well known that the clock products absorb partly the errors of predicted orbits, particularly errors of the line-of-sight (LOS) direction [23][24][25]. Therefore, it is more reasonable to assess orbit and clock residual error as a whole directly from a user's perspective. The differences between GF and GB wide-lane ambiguities is a starting point to investigate.
This contribution aims at analysis of single-site ambiguities derived from both GF and GB model. We focus on numerical computation and analysis using raw observations without differential and combination operator in the measurement domain. The GB WL ambiguities are derived using different orbit and clock products. To simplify the presentation and discussion, we are focusing on GPS in this contribution, but this does not affect the generality of our analysis for other constellations. We also investigate the discrepancies between GF and GB WL ambiguities to disclose the orbit and clock quality from a user's perspective. The contribution is arranged as follows: First, the GF and GB models are respectively presented in detail. Second, the effects of orbit and clock errors on the ambiguities are given in an analytical way. Third, the results of GB and GF ambiguities are numerically analyzed for both post-processing and real-time orbit and clock products. In addition, we simulate two types of errors in orbit and clock and reprocess the data and investigate the impact on WL ambiguities for further understanding. Finally, the study is summarized with some conclusions.

WL Ambiguities Derived from the GF Model
The observation equation for measurements of code p and phase ϕ tracked by a receiver at location r from GPS satellites s (s = 1, 2, · · · , m), at frequency i (i = 1, 2) at epoch k, can be expressed in meters as follows [26]: E{p s i,r (k)} = ρ s r (k) + m s r · T r + δ t r (k) − δ t s (k) + µ i I s r (k) + d p i ,r − d s p i E{ϕ s i,r (k)} = ρ s r (k) + m s r · T r + δ t r (k) − δ t s (k) − µ i I s r (k) where E{·} denotes the expectation operator; ρ s r (k) is the geometric range between the receiver and satellite; T r the zenith tropospheric delay with mapping function m s r ; δ t r (k) and δ t s (k) the receiver and satellite clock bias in meters, respectively; µ i = λ 2 i /λ 2 1 , λ i the carrier wavelength of frequency i; I s r (k) the ionospheric delay; d p i ,r and d ϕ i ,r are the receiver hardware delays (or biases) on code and phase, respectively; d s p i and d s ϕ i are the satellite hardware delays (or biases) on code and phase, respectively; where it is assumed that these satellite hardware delays are relatively stable in time; N s i,r the integer carrier ambiguity. All terms with (k) are supposed to be significantly changed in time.
The so-called GF observation model dispenses with geometric a priori information. The geometry part is obtained by combining all non-dispersive terms, including the geometric range ρ s r (k), satellite δ t s (k) and receiver clock δ t r (k) and tropospheric delay m s r · T r , even though the system is rank deficient. Not all terms are estimable. The hardware delays here are chosen as S-basis, absorbed by other terms [19]. Let us define the geometry part as one common term for all observations from one satellite: Sensors 2018, 18, 1760 4 of 16 and the ionosphere term is also combined with part of satellite and receiver hardware biases: and the ambiguity terms for L1 and L2 frequency are defined as: Therefore, the functional GF model is written as: where the model is for a single satellite. Similar expressions can be seen in [19,27] with the single-difference operator between satellites. Here, we focus on using raw measurements without any difference and combination operator in observation domain. The epoch-wise solution of the ambiguity term is simply computed as follows: where the phase and code measurements can be written in a vector form as z s T and the parameters that we are interested is ambiguity term H is the coefficient matrix derived from Equation (5) only for ambiguity parameters as follows: The ionosphere and geometry terms are eliminated, which is the origin that the term of GF actually comes from.
The stochastic model for the observation is given as follows: where e s is the elevation of the satellite. According to the error propagation law, the variance-covariance of ambiguity terms can be derived as: where Equations (6) and (8)  , the variance-covariance is first updated with the assumed process noise for epoch solution at k: where Q w is computed from the assumed process noise multiplying with time interval between epoch k and k−1, and then updated with the epoch-wise variance-covariance: with the parameter being also updated with the epoch-wise estimate and its variance-covariance: where elements ofx s k and Qxs k are used to derive the WL ambiguity, i.e.: with its formal error being derived: where σ 2

WL Ambiguities Derived from the GB Model
To derive the ambiguities from the GB model, we assume that the station coordinates (X r , Y r , Z r ), satellite orbit (X s (k), Y s (k), Z s (k)) and clock errors δ t s i f (k) are precisely known. As we know that the satellite clock errors are calculated from ionosphere-free combination observations [2]: (15) and the computed range between satellite and the known station coordinates is: Add Equation (15) to Equation (1), insert ρ s r (k) the left side, and define the observed minus computed values as: (16) resulting in the following: where the ionosphere and ambiguity terms are defined as the same as in Equations (3) and (4), respectively, and the receiver clock term is defined as: The observation equations of all observed satellites in the matrix-vector form can be expressed as: ; A k is the design matrix with its elements being formed with the coefficients in Equation (17) for all visible satellites.
The stochastic model is expressed as: where It is noticed that the observation Equation (18) has redundancy when the number of satellites is larger than three, which is easy to achieve even when only GPS satellites are employed. The state transition equation is given by: where B := diag 1 0 0 m e m e m with 0 m := 0 0 · · · 0 ,e m := 1 1 · · · 1 and w ∼ N(0, Q w ), which links the troposphere and ambiguity parameters with that of previous epochs.
Assume that the filtering solution at the epoch k−1 is denoted asˆ x k−1 with its covariance matrix Qˆ x k−1 , the solution at the epoch k is computed as: with its covariance matrix being computed as: whereˆ x k,k−1 = Bˆ x k−1 , and its covariance matrix as Qˆ x k,k−1 = BQˆ x k−1 B T + Q w , and the gain matrix is The L1 and L2 ambiguities with their covariance matrix can be used to derive the GB WL ambiguity and its formal error as Equations (13) and (14) for their GF counterparts.

Impact of Orbit and Clock Residual Errors on Ambiguity Estimation
From Equations (5) and (17), there is no difference between the GF and GB ambiguities in terms of parameterization. However, their estimates are different since they go through different estimation procedures. For the GF ambiguities, the geometry-related terms, including orbit and clock, are eliminated without effects. However it results in no redundancy for the epoch-wise estimates. The solution's accuracy is dominated by the noise of code measurements. In this case, the solution with multiple epochs is computed. For the GB ambiguities, the geometry strength is higher (as there are redundancy in each epoch), the ambiguities converge quicker than the GF ones. However, the observed minus computed values are computed from the orbit and clock, which are not always precise enough. The orbit and clock residual error at the LOS remains in the observation as follows: where oc s r indicates the slowly time-variant or time-invariant residual part of orbit and clock error combination in the LOS direction. As they are not random noise, they could be distributed to every parameters (as well as measurement residuals), but most likely to the ionosphere and ambiguity terms of the satellite as these parameters are satellite dependent. There are two cases worthy of investigation.

Orbit and Clock Residual Error Being Time-Constant Bias
It is possible that the orbit and clock contain a constant bias. For example, the code measurement type used in mobile side is not consistent with the code measurements used in the reference side when computing the clock products. It is assumed that the oc s r is a constant during the whole measurement period and there are no effects between satellites. As phase measurements are assigned much higher weight than code ones (100 times smaller in terms of standard deviation), we assume no significant orbit and clock errors are attributed to phase measurement residuals. Furthermore, the ambiguities are common terms between epochs, they likely absorb the constant bias as: The remaining error may be distributed to the ionosphere term and the code residuals. From the analytical expression of Equation (24), the effect on the WL ambiguity term is less than one tenth of oc s r in terms of cycle. In other words, the effect of constant bias in orbit and clock looks largely reduced. Bearing in mind that the reality can be more complex when our assumptions do not hold for all cases. In particular, the geometry changing may have influence as well. However, the effects of a constant bias analyzed here might be a minimum (i.e., a bottom line).

Orbit and Clock Residual Error Being Time-Variant Error
The residual orbit and clock error is rarely a constant bias in an entire tracking period, except for a mixture use of code measurements. It is very likely a time variant error that caused by the changes of LOS and the orbit error itself (as will be shown later). In this case, it is not a trivial to make a strict analytical expression, as the estimation procedure is rather complex both in current epoch and correlation between epochs. Assume that the error is directly mapped to the ambiguity terms as follows: where the effect on WL ambiguity term will be about 1.165 · oc s r cycle. In other words, one meter bias could bring 1.165 cycle error into the WL ambiguity. Even when the error is about 15 cm (often the case for real-time orbit and clock), the effects will be significantly up to 0.18 cycle, which is not acceptable for phase bias calibration for PPP IAR. In reality, the effects might be around (larger or smaller than) this level due to different geometry strength. This contribution will make a real-data numerical analysis in order to further understand the effects.

Experiment and Results Analysis
In order to analyze the performance of GF and GB WL ambiguities and evaluate the orbit and clock effects on the WL ambiguities, we implemented the GF and GB models using raw observations to compute L1, L2 and WL ambiguities as well as their variance-covariance as described in Sections 2 and 3. GNSS measurements with 1 s interval from selected IGS sites are downloaded from the IGS data base. These sites are deliberately selected and spread over main continents and the middle of the Pacific in terms of groups, see Figure 1. Some groups have sites located very nearby in order to have almost the same LOS. Two types of orbit and clock products are used in the experiment. One type is the determined orbit and clock products provided by ESA to the IGS, while the other is the real-time orbit and clock estimated as described in [15]. The interval of satellite orbit is 15 min, and the sample rate of satellite clock is 30 s for both types of products. The determined products are obtained by the post-processing mode, and their accuracy is better than that of real-time ones [15]. The standard GB PPP solutions, in parallel with the GF process, are computed using raw L1 and L2 observations. The parameters to be estimated are tropospheric zenith delay, ionosphere, receiver clock bias, and L1, L2 and then WL ambiguity terms for the GB process, while the GF process just produces the L1, L2 and then WL ambiguity terms. The positions of sites are fixed to the "ground truth", which are computed from many days' data with the static mode. This contribution makes use of the precise coordinates in order to determine precisely the GB-GF ambiguities by avoiding of the effects of unknown coordinates. The orbit and clock are integrated and interpolated for each epoch (1 s of data interval). Strictly speaking, the interpolation for real-time orbit and clock is not realistic in the real-time application as the future correction stream cannot be known in advance for interpolation of current epoch. However, we employ exactly the same process in order to carry out a fair comparison for both types of products by ignoring correction age effects. The elevation-dependent data weighting is used for both processes. The entire process is implemented through a strictly forward computation. the real-time orbit and clock estimated as described in [15]. The interval of satellite orbit is 15 min, and the sample rate of satellite clock is 30 s for both types of products. The determined products are obtained by the post-processing mode, and their accuracy is better than that of real-time ones [15]. The standard GB PPP solutions, in parallel with the GF process, are computed using raw L1 and L2 observations. The parameters to be estimated are tropospheric zenith delay, ionosphere, receiver clock bias, and L1, L2 and then WL ambiguity terms for the GB process, while the GF process just produces the L1, L2 and then WL ambiguity terms. The positions of sites are fixed to the "ground truth", which are computed from many days' data with the static mode. This contribution makes use of the precise coordinates in order to determine precisely the GB-GF ambiguities by avoiding of the effects of unknown coordinates. The orbit and clock are integrated and interpolated for each epoch (1 s of data interval). Strictly speaking, the interpolation for real-time orbit and clock is not realistic in the real-time application as the future correction stream cannot be known in advance for interpolation of current epoch. However, we employ exactly the same process in order to carry out a fair comparison for both types of products by ignoring correction age effects. The elevationdependent data weighting is used for both processes. The entire process is implemented through a strictly forward computation.

Analysis of GB and GF Ambiguities
This section will analyze the GB and GF ambiguities and inspect the similarity and discrepancies of GB and GF WL ambiguities. Figure 2 shows L1, L2 and WL ambiguities and their corresponding formal standard deviations for PRN 09. These ambiguities are respectively produced with the GF and GB models for site CEBR. Arbitrary offsets are removed from these ambiguities. It is obvious that both the variation of L1 and L2 GF ambiguities and their formal errors are larger than the GB counterparts, and especially they converge much slower than that of GB. However, the formal errors of WL ambiguities are very consistent and almost the same, meanwhile their estimates reveal a consistent trend despite of some small fluctuations. This confirms the analytical statement in [19] that the GF WL ambiguity precision does not differ too much from its GB counterpart. The reason is that the GF L1 and L2 ambiguities are highly correlated and even more correlated than GB ones. The WL ambiguities are therefore much precise after the differential operator, referring to (13) and (14), compared with the L1 and L2 ones. It is noticed that the GF WL ambiguities are slightly oscillated, compared with the GB ones at the first hour, but they are even stable than the GB ones afterwards. The reason is that the GB WL ambiguities may suffer some errors in orbit and clock. This will be further analyzed in other sections.

Analysis of GB and GF Ambiguities
This section will analyze the GB and GF ambiguities and inspect the similarity and discrepancies of GB and GF WL ambiguities. Figure 2 shows L1, L2 and WL ambiguities and their corresponding formal standard deviations for PRN 09. These ambiguities are respectively produced with the GF and GB models for site CEBR. Arbitrary offsets are removed from these ambiguities. It is obvious that both the variation of L1 and L2 GF ambiguities and their formal errors are larger than the GB counterparts, and especially they converge much slower than that of GB. However, the formal errors of WL ambiguities are very consistent and almost the same, meanwhile their estimates reveal a consistent trend despite of some small fluctuations. This confirms the analytical statement in [19] that the GF WL ambiguity precision does not differ too much from its GB counterpart. The reason is that the GF L1 and L2 ambiguities are highly correlated and even more correlated than GB ones. The WL ambiguities are therefore much precise after the differential operator, referring to (13) and (14), compared with the L1 and L2 ones. It is noticed that the GF WL ambiguities are slightly oscillated, compared with the GB ones at the first hour, but they are even stable than the GB ones afterwards. The reason is that the GB WL ambiguities may suffer some errors in orbit and clock. This will be further analyzed in other sections.  As can be seen from Figure 3 that the results of PRN 09 computed by another single receiver REDU have similar behavior to that of CEBR, which indicates that the WL ambiguities and their formal errors have commonality for this same satellite obtained by different receivers which are located not very far away.   As can be seen from Figure 3 that the results of PRN 09 computed by another single receiver REDU have similar behavior to that of CEBR, which indicates that the WL ambiguities and their formal errors have commonality for this same satellite obtained by different receivers which are located not very far away. As can be seen from Figure 3 that the results of PRN 09 computed by another single receiver REDU have similar behavior to that of CEBR, which indicates that the WL ambiguities and their formal errors have commonality for this same satellite obtained by different receivers which are located not very far away.      Figure 5 displays the values of GB-GF WL ambiguities computed with ESA final orbit and clock products at site KOKV for all visible GPS satellites. It is seen that these values are stable and close to zero. Again, the typical behavior of these time series is that they are noisy at the beginning of a tracking period, but become smooth and stable afterwards. According to the statistics, the biases are smaller than 0.05 cycle and the standard deviations of these values are smaller than 0.025 cycle. It is practically and commonly believed that the repeatability is good enough if differences of ambiguities are smaller than 0.1 cycle [7]. Therefore, we could state that the GB and GF WL ambiguities are equivalent in this case, though they are not the same.   Figure 5 displays the values of GB-GF WL ambiguities computed with ESA final orbit and clock products at site KOKV for all visible GPS satellites. It is seen that these values are stable and close to zero. Again, the typical behavior of these time series is that they are noisy at the beginning of a tracking period, but become smooth and stable afterwards. According to the statistics, the biases are smaller than 0.05 cycle and the standard deviations of these values are smaller than 0.025 cycle. It is practically and commonly believed that the repeatability is good enough if differences of ambiguities are smaller than 0.1 cycle [7]. Therefore, we could state that the GB and GF WL ambiguities are equivalent in this case, though they are not the same.  Figure 5 displays the values of GB-GF WL ambiguities computed with ESA final orbit and clock products at site KOKV for all visible GPS satellites. It is seen that these values are stable and close to zero. Again, the typical behavior of these time series is that they are noisy at the beginning of a tracking period, but become smooth and stable afterwards. According to the statistics, the biases are smaller than 0.05 cycle and the standard deviations of these values are smaller than 0.025 cycle. It is practically and commonly believed that the repeatability is good enough if differences of ambiguities are smaller than 0.1 cycle [7]. Therefore, we could state that the GB and GF WL ambiguities are equivalent in this case, though they are not the same.

Orbit and Clock Residual Error Effects
The results demonstrated in previous section are computed with the ESA final orbit and clock products. To compare with them, this section will show results that are computed with the real-time orbit and clock products. Figure 6 displays various differences of PRN 09 WL ambiguities and orbit clock products computed at the site CEBR and REDU, respectively. In the top panels, the GF WL differences are zeros as the orbit and clock are eliminated and have no impact on the GF ambiguity estimates. The GB WL differences are computed with final and real-time orbit and clock. In this case, everything is the same except that the orbit and clock used are different. Therefore, the GB WL differences indicate that the effects of the different orbit and clock products. To further confirm this, the differences of two types of orbit and clock are also computed for the locations of the receivers, in which the orbits are mapped into the LOS direction. The values are converted to cycles by divided with the wavelength of WL ambiguities (which is 86.2 cm). It is clear that the variations of two sets of orbit and clock differences largely coincide with the differences of GB WL ambiguities. In the bottom panels, the GB minus GF WL ambiguities used two types of orbit and clock products are shown together with the orbit and clock differences themselves. The values of GB-GF computed from the final orbit and clock (OC1: GB-GF in the figure) are more stable than that of computed from the real-time orbit and clock (OC2: GB-GF in the figure). The latter one has a clear slope in the time series, which is consistent with the slope of orbit and clock differences. In this case, it is difficult to say the GF WL ambiguities are equivalent with the GB ones. However, it is promising for the real-time users that the GB-GF WL ambiguities can be used to an indicator of the orbit and clock quality from a user's perspective.

Orbit and Clock Residual Error Effects
The results demonstrated in previous section are computed with the ESA final orbit and clock products. To compare with them, this section will show results that are computed with the real-time orbit and clock products. Figure 6 displays various differences of PRN 09 WL ambiguities and orbit clock products computed at the site CEBR and REDU, respectively. In the top panels, the GF WL differences are zeros as the orbit and clock are eliminated and have no impact on the GF ambiguity estimates. The GB WL differences are computed with final and real-time orbit and clock. In this case, everything is the same except that the orbit and clock used are different. Therefore, the GB WL differences indicate that the effects of the different orbit and clock products. To further confirm this, the differences of two types of orbit and clock are also computed for the locations of the receivers, in which the orbits are mapped into the LOS direction. The values are converted to cycles by divided with the wavelength of WL ambiguities (which is 86.2 cm). It is clear that the variations of two sets of orbit and clock differences largely coincide with the differences of GB WL ambiguities. In the bottom panels, the GB minus GF WL ambiguities used two types of orbit and clock products are shown together with the orbit and clock differences themselves. The values of GB-GF computed from the final orbit and clock (OC1: GB-GF in the figure) are more stable than that of computed from the realtime orbit and clock (OC2: GB-GF in the figure). The latter one has a clear slope in the time series, which is consistent with the slope of orbit and clock differences. In this case, it is difficult to say the GF WL ambiguities are equivalent with the GB ones. However, it is promising for the real-time users that the GB-GF WL ambiguities can be used to an indicator of the orbit and clock quality from a user's perspective.    To further investigate the cause of the variation in GB-GF WL ambiguities, we plot the variations of differences between the final and real-time orbit clock (orbit is projected into the LOS) at Site KOKV for all visible GPS satellites in Figure 8. It is clearly seen that the pattern of the orbit and clock differences is very similar with that of GB-GF ambiguities, in particular, for the satellite PRN 09. The reason for such an extreme orbit and clock behavior in real-time is likely that the orientation of the satellite (mainly IIF satellites) is estimated or predicted wrongly [22]. This issue has not been solved in real-time case. Therefore, it is valuable that the variation of GB-GF values of a receiver can be taken as a quality index of orbit clock products from users' perspective. Figure 9 displays GB-GF WL ambiguities of PRN 09 computed with real-time orbit clock products for all selected sites. It is obvious that these time series are not stable. The figure discloses longer period of oscillations than one single site shown in Figure 7. We could also see that the variation patterns are almost the same for some pairs of sites, e.g., KOKB and KOKV, KIRI and NAUR, MIZU and SANT, MRO1 and PERT. These pairs of sites have different observations, but demonstrate the same GB-GF WL patterns, as they suffer almost the same orbit clock error. It is clear in this case that the GB and GF WL ambiguities are not equivalent. To further investigate the cause of the variation in GB-GF WL ambiguities, we plot the variations of differences between the final and real-time orbit clock (orbit is projected into the LOS) at Site KOKV for all visible GPS satellites in Figure 8. It is clearly seen that the pattern of the orbit and clock differences is very similar with that of GB-GF ambiguities, in particular, for the satellite PRN 09. The reason for such an extreme orbit and clock behavior in real-time is likely that the orientation of the satellite (mainly IIF satellites) is estimated or predicted wrongly [22]. This issue has not been solved in real-time case. Therefore, it is valuable that the variation of GB-GF values of a receiver can be taken as a quality index of orbit clock products from users' perspective. Figure 9 displays GB-GF WL ambiguities of PRN 09 computed with real-time orbit clock products for all selected sites. It is obvious that these time series are not stable. The figure discloses longer period of oscillations than one single site shown in Figure 7. We could also see that the variation patterns are almost the same for some pairs of sites, e.g., KOKB and KOKV, KIRI and NAUR, MIZU and SANT, MRO1 and PERT. These pairs of sites have different observations, but demonstrate the same GB-GF WL patterns, as they suffer almost the same orbit clock error. It is clear in this case that the GB and GF WL ambiguities are not equivalent.

Simulation Study of Orbit/Clock Error Effects
In order to further investigate the effects of orbit and clock errors on the WL ambiguities, we simulate two types of errors directly in the LOS for one satellite of an entire tracking period at Site KOKV: (1) a constant bias (e.g., 1.0 m); (2) time-variant errors with 0.15 m amplitude and 2 h period of sine or cosine changes. The errors are added to only PRN 30 of the ESA final orbit and clock. The GB-GF WL ambiguities are then recomputed and compared with the previous ones. The purpose for such a simulation is to see how effects of precisely "known errors" on the GB-GF ambiguities. Figure 10 displays differences of recomputed GB-GF WL ambiguities with respect to their counterparts without the simulated errors. The differences are purely caused by the simulated errors. The top panel displays the GB-GF WL effects of a 1.0 m constant bias is added into the PRN 30 during 5:00 h-11:30 h. The effect is about a constant bias of 0.084 cycle. There is no visible effects on other satellites. This indicates that the GB-GF WL ambiguity is not very sensible to a constant bias, which confirms the analysis in "orbit and clock residual error being constant bias".

Simulation Study of Orbit/Clock Error Effects
In order to further investigate the effects of orbit and clock errors on the WL ambiguities, we simulate two types of errors directly in the LOS for one satellite of an entire tracking period at Site KOKV: (1) a constant bias (e.g., 1.0 m); (2) time-variant errors with 0.15 m amplitude and 2 h period of sine or cosine changes. The errors are added to only PRN 30 of the ESA final orbit and clock. The GB-GF WL ambiguities are then recomputed and compared with the previous ones. The purpose for such a simulation is to see how effects of precisely "known errors" on the GB-GF ambiguities. Figure 10 displays differences of recomputed GB-GF WL ambiguities with respect to their counterparts without the simulated errors. The differences are purely caused by the simulated errors. The top panel displays the GB-GF WL effects of a 1.0 m constant bias is added into the PRN 30 during 5:00 h-11:30 h. The effect is about a constant bias of 0.084 cycle. There is no visible effects on other satellites. This indicates that the GB-GF WL ambiguity is not very sensible to a constant bias, which confirms the analysis in "orbit and clock residual error being constant bias". The bottom panel demonstrates the GB-GF WL effects after the time-variant errors are added into the PRN 30 during 5:00 h-11:30 h. Firstly, the effects of PRN 30 show also the sine or cosine behavior. Secondly, there are also similar patterns with smaller magnitude in many other satellites tracked in the same period. These effects are brought by the added errors with only 0.15 m in (absolute) size, which is much smaller than 1.0 m. This indicates that the GB-GF WL ambiguity is very sensible to time-variant errors as analyzed in Section 4.1, which provides an opportunity for real-time user to evaluate the quality of real-time orbit and clock.

Conclusions
This contribution derives L1, L2 and WL ambiguities from the GF and GB PPP model. Two types of orbit and clock products, represented different level of quality, are used to compute the WL ambiguities. The differences of GB and GF WL ambiguities are particularly analyzed. The main results and conclusions are summarized in the following statements: The GF ambiguities are obtained on an individual satellite basis, and there is no redundancy for the epoch-wise estimation. To improve GF estimates, multiple data epochs are needed. The GB The bottom panel demonstrates the GB-GF WL effects after the time-variant errors are added into the PRN 30 during 5:00 h-11:30 h. Firstly, the effects of PRN 30 show also the sine or cosine behavior. Secondly, there are also similar patterns with smaller magnitude in many other satellites tracked in the same period. These effects are brought by the added errors with only 0.15 m in (absolute) size, which is much smaller than 1.0 m. This indicates that the GB-GF WL ambiguity is very sensible to time-variant errors as analyzed in Section 4.1, which provides an opportunity for real-time user to evaluate the quality of real-time orbit and clock.

Conclusions
This contribution derives L1, L2 and WL ambiguities from the GF and GB PPP model. Two types of orbit and clock products, represented different level of quality, are used to compute the WL ambiguities. The differences of GB and GF WL ambiguities are particularly analyzed. The main results and conclusions are summarized in the following statements: The GF ambiguities are obtained on an individual satellite basis, and there is no redundancy for the epoch-wise estimation. To improve GF estimates, multiple data epochs are needed. The GB ambiguity estimation makes use of all received satellites, therefore the geometry and the redundancy, besides using multiple epochs of data, play roles on the estimation. From our implementation and experiment, there are significant differences between GB L1/L2 ambiguities and the GF ones in terms of both estimates and their precisions, however, the GB WL ambiguities do not differ much from the GF ones in terms of precision at least.
Orbit and clock have no impact on GF WL ambiguities as they are eliminated in the process, while they have effects on GB ones. Whether the GB and GF ambiguities are equivalent depends much on the quality of the used orbit and clock. In general, if GB-GF WL values are stably in between ±0.1 cycles or within three times of their formal errors, we could state that the GB and GF WL ambiguities are equivalent. This is basically not an issue for the final orbit and clock product used in our experiment, as this type of product is obtained from post-processing.
The GB and GF WL ambiguities are not always equivalent if they are computed from the real-time orbit and clock, as the type of product is sometimes not perfect. It is found in this contribution that the variation of GB-GF WL ambiguities coincide with the orbit and clock differences between the final and real-time products. Therefore, it is promising for the real-time users that the GB-GF WL ambiguities can be directly used to indicate the orbit and clock quality from a user's perspective.
Two types of orbit and clock errors are simulated and the effects on GB-GF WL ambiguities are analyzed. The constant bias in orbit and clock seems not have significant effects on GB-GF WL values if we look into the pattern of their time-series, though there is still a smaller size of bias presenting. This indicates that the GB-GF WL ambiguity is not very sensible to a constant bias in orbit and clock. However, the GB-GF WL ambiguity varies remarkably if the time-variant errors exist in the orbit and clock. Even the errors are not large, it is sensitively mapped into GB-GF WL ambiguities. As this type of errors is more often showing up than the constant bias in the real-time orbit and clock, the GB-GF WL ambiguity provides an opportunity for real-time users to evaluate the quality of real-time orbit and clock.
Author Contributions: G.C. provided the initial concept and methodology; Q.Z. helped to derive the equations in theory with G.C.; G.C. acquired the funding for the project and carried out the data processing. Q.Z. provided the software platform and results analysis. S.L. implemented the software part necessarily for this project and carried out the experimental validation; S.L. visualized the results and wrote the draft of the paper. G.C. and Q.Z. finalized the paper. All authors reviewed the paper.