Estimation of fractional cycle bias for GPS/BDS-2/Galileo based on international GNSS monitoring and assessment system observations using the uncombined PPP model

The Fractional Cycle Bias (FCB) product is crucial for the Ambiguity Resolution (AR) in Precise Point Positioning (PPP). Different from the traditional method using the ionospheric-free ambiguity which is formed by the Wide Lane (WL) and Narrow Lane (NL) combinations, the uncombined PPP model is flexible and effective to generate the FCB products. This study presents the FCB estimation method based on the multi-Global Navigation Satellite System (GNSS) precise satellite orbit and clock corrections from the international GNSS Monitoring and Assessment System (iGMAS) observations using the uncombined PPP model. The dual-frequency raw ambiguities are combined by the integer coefficients (4,− 3) and (1,− 1) to directly estimate the FCBs. The details of FCB estimation are described with the Global Positioning System (GPS), BeiDou-2 Navigation Satellite System (BDS-2) and Galileo Navigation Satellite System (Galileo). For the estimated FCBs, the Root Mean Squares (RMSs) of the posterior residuals are smaller than 0.1 cycles, which indicates a high consistency for the float ambiguities. The stability of the WL FCBs series is better than 0.02 cycles for the three GNSS systems, while the STandard Deviation (STD) of the NL FCBs for BDS-2 is larger than 0.139 cycles. The combined FCBs have better stability than the raw series. With the multi-GNSS FCB products, the PPP AR for GPS/BDS-2/Galileo is demonstrated using the raw observations. For hourly static positioning results, the performance of the PPP AR with the three-system observations is improved by 42.6%, but only 13.1% for kinematic positioning results. The results indicate that precise and reliable positioning can be achieved with the PPP AR of GPS/BDS-2/Galileo, supported by multi-GNSS satellite orbit, clock, and FCB products based on iGMAS.


Introduction
Precise Point Positioning (PPP) Ambiguity Resolution (AR) is important to obtain the positioning accuracy at centimeter-level in a short time at a station (Bisnath and Gao 2008;Wang et al. 2019). Due to the Fractional Cycle Biases (FCBs) in phase measurements which are assimilated into the undifferenced ambiguities in PPP, the integer properties of the estimated ambiguities are lost (Gabor and Nerem 1999). Currently, it is assumed that FCBs contain the biases from satellite and receiver hardware delays. Unlike the relative positioning model, where FCBs are eliminated in double-differenced Global Navigation Satellite System (GNSS) observations (Ge et al. 2008;Geng et al. 2012), the PPP processing needs the FCB products to correct the float ambiguities for recovering their integer properties.
The method of ambiguity resolution based on the satellite-satellite single differences was first proposed in PPP (Gabor and Nerem 1999), many methods have then been developed with different products to obtain the ambiguity-fixed positioning solutions (Ge et al. 2008;Collins et al. 2008;Laurichesse et al. 2009). Based on an empirical assumption that the Uncalibrated Phase Delays (UPD) are relatively stable in time, the Single-Difference (SD) FCBs of the Wide Lane (WL) and Narrow Lane (NL) float ambiguities between satellites are estimated from a reference network (Ge et al. 2008). The WL FCBs are determined by averaging the fractional parts of all WL ambiguities from the same satellites. The Melbourne-Wübbena (MW) measurements are used to derive the wide-lane ambiguities (Melbourne 1985;Wübbena 1985). Generally, the WL FCBs are stable over a few days, or even several months (Gabor and Nerem 1999). After determining the WL FCBs, the integer part of the WL ambiguities is used to estimate the float narrow-lane ambiguities. Similarly, the NL FCBs are also determined by averaging the fractional parts of the NL ambiguities. Due to the short wavelength of the narrow-lane ambiguity, for instance, about 10 cm for Global Positioning System (GPS), the NL FCBs are not as stable as the WL, which are proposed to estimate 15 min mean values.
Instead of estimating the SD FCB from the same satellite pairs, Li and Zhang (2012) proposed that the Zero-Difference (ZD) FCBs can be estimated together in the adjustment system using the least-square method with all FCB measurements (Li and Zhang 2012). This approach is adopted for the WL and NL FCBs estimation, which significantly improved the accuracy of estimated FCBs. For the NL FCBs, the NL float ambiguities derivated from the ionospheric-free float ambiguities in PPP and the WL integer values, are directly used in the FCB estimator. The low accuracy of the NL float ambiguities degrades the performance of FCB estimation. The ambiguity-fixed solutions from a GNSS network, where their double-difference ambiguities are got, are used for the NL FCB estimation (Geng et al. 2012). Additionally, the Kalman filter is also adopted for the FCB estimation epoch by epoch, which significantly speeds up the computation and is suitable for real-time applications (Xiao et al. 2018).
Distinguished from the FCB model proposed by Ge et al. (2008) and Laurichesse et al. (2009) presented the integer phase clock model, which consists of the ionospheric-free combination and the MW combination (Laurichesse et al. 2009). The NL FCBs are not estimated but assimilated into the estimates of satellite clock offsets. Hence, the NL ambiguities in the network solutions are directly fixed to the nearest integers, and the FCBcontained clock products based on the ambiguity-fixed solutions are determined (Laurichesse 2011). Similarly, Collins et al. (2008) estimated the integer recovery clocks for pseudorange and phase measurements, which is named as the decoupled clock model (Collins et al 2008). In this model, the effects from time-varying parts of code biases on the fractional parts of ambiguities are rigorously considered. The integer recovery clocks proposed by Collins et al. (2008) and Laurichesse et al. (2009) are used to recover the integer characteristics for the NL float ambiguities in PPP for ambiguity resolution. These two methods assimilated the NL FCBs into clocks estimates, which increased the computation of clock products. In the FCB method, the original International GNSS Service (IGS) precise clock products are used to separate the FCBs from ambiguities in an independent processing. Hence, the FCB method is widely adopted in many studies.
Most of the FCB estimation results are derived with the ionospheric-free PPP model. Recently, the uncombined PPP model with raw pseudorange and phase observations has attracted a great attention. The uncombined PPP model can directly estimate the raw ambiguities on each frequency and is flexible for processing multi-frequency observations. Its performance was demonstrated in the single-or dual-frequency PPP and PPP-RTK (Real-Time Kinematic) (Wübbena et al. 2005;Li et al. 2011;Zhang et al. 2011;Chen et al. 2015;Teunissen and Khodabandeh 2015;Lou et al. 2016). Li et al. (2013) presented the FCB estimates with the uncombined PPP model using the GPS L1 and L2 raw observations, which improved the performance of generating FCBs, and consequently improved the positioning accuracy of the PPP AR (Li et al. 2013). Gu et al. (2015a) determined the FCBs in the WL and NL combinations using the uncombined PPP model and demonstrated the ionosphere characteristics in the PPP AR (Gu et al. 2015a). This approach is also used for the triple-frequency FCBs estimation of BeiDou-2 Navigation Satellite System (BDS-2) (Gu et al. 2015b;Li et al. 2018). Xiao et al. (2019) further presented a triple-frequency FCB model for the PPP ambiguity resolution with raw observations, which is verified with Galileo Navigation Satellite System (Galileo) and BDS-2. The initial positioning performance of the PPP AR with the uncombined model is also demonstrated (Wang et al. 2020).
For implementing the PPP AR, the Centre National d'Etudes Spatiales (CNES) in French and the School of Geodesy and Geomatics at Wuhan University (SGG-WHU) have been routinely generating the integer recovery clocks and FCBs for public PPP applications (Loyer et al. 2012;Li et al. 2016). Recently, SGG presented the PPP AR with GPS, BeiDou Navigation Satellite System (BDS), Galileo, and Quasi-Zenith Satellite System (QZSS) with multi-GNSS FCBs using the precise satellite orbits and clocks from different IGS centers (Hu et al. 2020). Geng et al. (2019b) proposed a modified phase clock/bias model to improve the PPP AR and provided the phase clocks and daily phase biases based on the Center for Orbit Determination in Europe (CODE) precise satellite orbit products to users who process the PPP with software "PRIDE PPP-AR" (Geng et al. 2019a, b). Hence, the generation of FCBs for the PPP AR must be an essential work for supporting PPP services.
The international GNSS Monitoring and Assessment System (iGMAS), which is a scientific project initiated and led by China under the United Nations framework, provides the daily satellite precise products for PPP users to get ambiguity-float solutions (Jiao et al. 2012). To meet the requirement of stable and reliable positioning, the ambiguity-fixed solutions should be implemented with the corresponding FCB products. In our work, the uncombined PPP model is adopted to generate the FCBs with the final satellite orbit and clock corrections with GPS, Galileo, and BDS-2 observations at the Bei-Dou analysis and service center of Chang'an University, which is a member of the iGMAS. The uncombined PPP model is flexible for PPP processing with multi-frequency and multi-GNSS observations. The raw ambiguities are directly estimated in the PPP model and then combined to estimate the phase bias in the FCB estimator. The initial results of the FCB generation and the PPP AR with GPS/BDS-2/Galileo data are analyzed to present the ambiguity-fixed positioning performance based on the iGMAS. This contribution presented the great advantages of the multi-GNSS PPP in terms of convergence speed and positioning accuracy. The PPP AR based on iGMAS products showed the potential for global users in single station positioning services. In this paper, the method for FCB estimation and the PPP AR are introduced in detail. The integer coefficients (4, − 3) and (1, − 1) are adopted to decorrelate the involving parameters, reducing the effects of ionospheric delays. Then, the processing strategies and experiment data are introduced. To evaluate the PPP AR performances, the precision of the estimated FCBs for the three GNSS systems is demonstrated. The daily and hourly positioning accuracy is analyzed in static and kinematic modes. Finally, conclusions and outworks are discussed.

Ambiguity-float PPP
In the uncombined PPP model, the ionospheric delays can be either estimated or corrected. The linearized observation equations for pseudorange and phase observations from satellite s to receiver r are described as: where P q,s r,f and L q,s r,f are the respective pseudorange and phase measurements on the frequency f (f = 1,2), from which the computed values are removed; u q,s r is the receiver-to-satellite unit vector; x is the vector of the receiver position corrections to its preliminary position; dt q r and dt q,s are the receiver and satellite clock errors, respectively; c is the light speed in vacuum;M q,s r is the elevation-dependent mapping function for the tropospheric wet delay from the corresponding zenith one Z r ; I q,s r,1 is the ionospheric delay along the line-of-sight from a receiver to a satellite at the first frequency and γ q,s f are the satellite hardware delays of code and phase observations, respectively; ε P,f and ε L,f are the code and phase measurement noises, respectively, which include the multipath effects (Shi and Gao 2014).
Considering that the ionospheric-free combined observations which contain the satellite code biases are adopted for satellite clock error (parameter) estimates, the satellite clock parameters can be denoted as: Hence, the receiver clock parameter is denoted as: The ionospheric delays and ambiguities are reparameterized as: where DCB r = d r,1 − d r,2 and DCB s = d s 1 − d s 2 are the Differential Code Biases (DCBs) for respective receiver and satellite. After applying the satellite clock corrections, the observations for the uncombined PPP model with GPS, BDS-2, and Galileo data can be written as: where ISB C r and ISB E r are the inter-system bias for BDS-2 and Galileo with respect to GPS, respectively. Normally, the ISBs are regarded as constant which reflects the stability of receiver hardware code delays. However, the datum biases of the satellite clocks between different GNSS systems are also assimilated into the ISBs. The stability of ISBs is strongly correlated to specific satellite clock products. Correspondingly, the variations of ISBs are regarded as a white noise process. All estimated parameters in the Multi-GNSS PPP model are expressed as vector X: where Ĩ r,1 denotes the vector of the ionospheric delay parameters for all observed satellites and Ñ r,f denotes the vector of the ambiguity parameters.
Additionally, the stations in a reference network are used to extract the float ambiguities for FCBs estimation. The stations' coordinates are fixed to their references which are obtained in IGS SINEX files or from the preliminary static PPP processing for the stations not listed in SNX files.

Fractional cycle bias estimation
The ambiguity parameter in Eq. (4) can be rewritten as: where N q,s r,f is an integer value, B q r,f and B q,s f are the respective receiver and satellite FCBs. Due to a strong correlation between ionospheric delays and ambiguities, which can be seen from Eq. (6), the accuracy of the estimated ambiguities on each frequency will be degraded by ionospheric errors. The combinations with integers (4,−3) and (1,−1), which have the low noise and long wavelengths, are selected to improve the accuracy of the estimated ambiguities and reduce the impacts of ionospheric errors (Li et al. 2018;Xiao et al. 2019). Hence, the ambiguities are combined as: Here, the combined ambiguities, Ñ q,s r,(4,−3) and Ñ q,s r,(1,−1) , are defined as new NL and WL ambiguities. The corresponding FCBs are reformed as: Since the ambiguities have the same structure in Eq. (7), the fractional parts of ambiguities can be formulated as: where n denotes the FCB measurement which is the fractional part of the real-value ambiguity solution Ñ s r . N s r presents the integer part of the real-value ambiguity, which contains the original integer ambiguity and the integer part of the code and phase delays from satellite to receiver; B r and B s are the FCBs for receiver r and satellite s, respectively. All FCB measurements in Eq. (11) from a reference network of m stations and n satellites tracked can be expressed as: In the coefficient vector R r , the rth element is 1 and the others are zero. In the coefficient vector S s , the sth element is -1 and the others are zero. One satellite for each GNSS system and new combined WL or NL ambiguities, respectively, is selected as a datum whose FCB is fixed to zero for resolving the rank deficiency in Eq. (12).
For the multi-GNSS data, the FCBs of new WL and NL combinations can be estimated together with GPS/ BDS-2/Galileo. To reduce the high computation load, the FCBs of the WL and NL combinations are estimated system by system. The individual FCBs are recovered by the inverse operation in Eq. (10). The multi-GNSS FCBs are integrated in one file by the predefined formats. The detailed flowchart is presented in Fig. 1.
Firstly, the static Uncombined PPP (UPPP) is processed for each reference station with its coordinates fixed to estimate the precise real-value ambiguities on each frequency of observed satellites. In this study, the FCBs are estimated with an interval of 15 min. Hence, in the quality checking of the float ambiguity solutions, the ambiguities estimated with the observations of less than 10 min are deleted. After the quality checking, the ambiguities are inputted into the FCB estimator. Then, the individual ambiguities on each frequency are combined using Eq. (9). The FCBs measurements in Eq. (11) are adopted in the FCBs estimation for the WL and NL combinations. Finally, the inverse operation in Eq. (10) is conducted to recover the individual FCBs on each frequency. The raw FCB on each frequency is flexible in the State Space Representation (SSR) of Radio Technical Commission for Maritime service (RTCM) for users' PPP AR (Shi 2012).

PPP AR at the user terminal
After correcting the satellite FCBs for the estimation of float ambiguities, the SD PPP ambiguities are proposed in the ambiguity resolution for removing the receiver FCB. Additionally, the combined ambiguities, as shown in Eq. (9), are recommended for the FCB estimation. Hence, the WL and NL combinations of ambiguities are sequentially fixed. To eliminate the effects of measurement noise and multipath, the satellite elevation angle adopted in the ambiguity resolution should not be less than 15°. The ambiguity fixing success rate can be further improved by the partial ambiguity resolution method. Firstly, all available ambiguities are decorrelated. The reformed ambiguities are reordered in the ascending order according to their decorrelated variances. They can also be reordered according to the satellite elevation angles. Secondly, the Least-square AMBiguity Decorrelation Adjustment (LAMBDA) method is adopted to search for the optimal integer values. Thirdly, the bootstrapped success rate P and the ratio test value R are calculated. If P<P 0 or R < R 0 , the last ambiguity with the lowest precision in the subset is removed and the second step is repeated. If the number of available ambiguities is less than four, the ambiguity-fixed solution fails. Otherwise, the integer ambiguities with a higher success rate and ratio value are confirmed as true values. Generally, the thresholds for P 0 and R 0 are set 0.999 and 2.0 (Li and Zhang 2015).
Once the integer values of ambiguities are confirmed, the tight constraint is imposed on the estimation of float ambiguities: Using Eq. (13), the ambiguity-fixed solutions for positioning can be obtained after the PPP reprocessing. Note that the constraint of integer ambiguity can also increase the precision of other estimated parameters.

Data and processing strategies
To assess the performance of the FCB estimations with the satellite orbit and clock products at Chang'an University based on the iGMAS, 316 stations from iGMAS and Multi-GNSS EXperiment (MGEX) are selected, shown in Fig. 2.
GPS satellites are tracked at all stations, and Galileo satellites are tracked at 190 stations, among which 138 stations tracked BDS-2 satellites. The uncombined PPP is implemented with GPS, Galileo, and BDS-2 observations for extracting the float ambiguities as the measurements in the FCB estimation. After generating FCBs, the uncombined PPP AR is achieved for GPS/BDS-2/Galileo combined solutions. The elevation-dependent stochastic model σ 2 =σ 2 0 /sin 2 (el) is adopted for the GNSS observations, where σ 0 is 0.003 m and 0.3 m for carrier phase and code measurements, respectively, and el is the satellite elevation angle in radian. Among GPS, Galileo and BDS-2, the ratio of their observation noise is set as 1:1:3. The Phase Center Offset (PCO) and the Phase Center Variations (PCV) published by the European Space Agency (ESA) are used to correct the observations for BDS-2 satellites, while for the others the corrections are according to the IGS14.atx file. In PPP processing, the cutoff angle of satellite elevation is 7° in the float solutions and 15° in the AR, while it is 30° in the FCB estimation for quality assurance. For the parameters estimated in PPP, the receiver clock and Inter-System Bias (ISB) are treated as white noise processes, the zenith tropospheric delay as a random walk process, and the ionospheric delays as white noise processes.

Results and analysis
The FCB products on Day of Year (DOY) 244, 2019 are generated to assess the performance of FCB estimation. Firstly, the stability of the FCB series over one day is analyzed with the averaged standard deviations for the combined FCBs and raw FCBs on each frequency. The accuracy of FCB estimation is evaluated by the posterior residuals of FCB measurements. After analyzing the performance of estimated FCBs, the uncombined PPP results are presented with the ambiguity-float and fixed solutions.

FCB residuals distributions
In the proposed method of FCB estimation, the FCBs for receiver and satellites are estimated separately using the WL and NL ambiguities. The FCB estimator directly provides the WL and NL combined FCBs. Hence, the accuracy assessment of the FCBs is carried out in the WL and NL combinations. The distribution of the posterior residuals is an intuitional indicator for the consistency of the FCBs results among different measurements. Generally, the residuals are close to zero indicating a high In the ionospheric-free model, the MW combination is commonly adopted to obtain the ambiguities for the WL FCB estimation, in which the accuracy of ambiguities is decreased by the averaging filter process. Compared with the FCB estimation in the ionospheric-free model, the WL combination reformed from the raw ambiguities on each frequency is free of pseudorange measurement noise and multipath. Generally, the residuals of the WL combination with longer wavelengths are smaller than the NL combination. For GPS satellites, the wavelength is about 86 cm for the WL combination and about 10 cm for the NL combination which is more sensitive to the errors. For GPS, Galileo, and BDS-2, this is verified by the RMS of the WL and NL residuals in Figs. 3, 4 and 5. The RMS of the WL residuals is 0.069, 0.046 and 0.085 cycles for GPS, Galileo, and BDS-2, while it is 0.086, 0.087, 0.106 cycles for the NL residuals, respectively. The RMS of the residuals is around or less than 0.1 cycles, which indicates a high consistency among the estimated FCBs. Additionally, for the WL and NL residuals, the RMS for BDS-2 is larger than that for GPS and Galileo, which indicates that the accuracy of satellite orbit and clock products is crucial for the FCB estimation. For the NL combination, the RMS of GPS residuals is the smallest which is reasonable because of its precise ambiguity float solutions in PPP. For the distribution of the NL combination, 92.7%, 92.4%, and 88.4% of the residuals are within [− 0.15, 0.15] (in cycles) for GPS, Galileo and BDS-2, respectively, while that is 96.1%, 99.0%, 91.1% for the WL combination. The distributions also suggest that the high consistency of the FCB measurements ensures the good accuracy of FCB products.  Compared with individual FCBs on each frequency, the stability is significantly better for the combined FCBs, especially for the WL combination. For the WL combination, the STD is 0.015, 0.005, 0.011 cycles for respective GPS, Galileo, and BDS-2. The corresponding STD for the NL combination is 0.039, 0.049, 0.139 cycles, which is larger than that for the WL combination. For individual FCB on each frequency, its fluctuation is remarkable with larger STDs. Hence, the more precise and stable FCB measurements from PPP processing are reasonable and more suitable for the FCB estimation. For the BDS-2 shown in Fig. 8, the stability of the FCBs is much poorer than other systems except for the WL combinations. At present, the precision of BDS-2 satellite ephemeris is not as good as GPS or Galileo. Consequently, the accuracy of the float ambiguities in BDS-2 PPP is degraded and the RMS of residuals in BDS-2 FCB estimation is larger than other GNSS systems as shown in Fig. 5. Hence, the mean of STDs for the BDS-2 narrowlane FCBs is larger than 0.1 cycles.

GPS/Galileo/BDS-2 PPP results
The GPS/BDS-2/Galileo PPP performance is investigated by the ambiguity-float or -fixed solutions with our estimated FCB products. 24 stations in Fig. 9 from the iGMAS/MGEX network are selected to process the uncombined PPP and get the ambiguity resolutions with GPS, Galileo and BDS-2 dual-frequency observations. The daily static PPP solutions are presented to show the final positioning accuracy and the superiority of the PPP AR. Then, the hourly solutions are obtained to investigate the convergence time with different time lengths as well as the ambiguity success fixing rate. The kinematic PPP AR performance is also evaluated with daily observations. Note that the reference coordinates for each station are the means of one-week static daily positioning solutions.
To study the stability of the PPP AR solutions, 3-h positioning errors from the one-day solution series at station YARR are presented in Fig. 10. Compared with the ambiguity-float solutions, the PPP AR achieves higher accuracy in a short period and gives more stable solutions. This indicates that the PPP AR can fast provide a reliable solution to users.
The RMS of the static daily positioning results for GPS/ BDS-2/Galileo PPP is shown in Fig. 11.
The RMS of the ambiguity-float solutions is 0.15, 0.06, and 0.20 cm for respective east, north, and up directions, while that is 0.14, 0.06, and 0.17 cm for the ambiguity-fixed solutions. The results indicate that the PPP Hourly solutions are obtained for the PPP AR in a short period of time, shown in Fig. 12.
The ambiguity success fixing rate is defined as the ratio of number of fixed solutions to that of all solutions. At different epochs, significant improvements in positioning RMS are achieved with the PPP AR. With the observations of 10 min, the Three-Dimensional (3D) RMS is 12.72 cm for the ambiguity-float solutions, and 11.65 cm for the PPP AR. The ambiguity success fixing rate is only 72.3%. For the observations of 15 min, the 3D RMS is 10.37 and 8.41 cm for the respective ambiguity-float and -fixed solutions, and the ambiguity success fixing rate is 82.4%. Finally, in hourly static solutions, the 3D RMS is 3.66 cm for ambiguity-float solutions and 2.10 cm for the PPP AR, improved by 42.6%, and the ambiguity success fixing rate is 97.7%. The positioning RMSs and improvements for respective east, north, and up directions are presented in Table 1.
The performance of the kinematic PPP with GPS/ BDS-2/Galileo data is also analyzed at the 24 stations. In Fig. 13, the initial 3-h kinematic positioning errors from the one-day solutions of station YARR show that the ambiguity-fixed solutions reduce significantly the fluctuation, realizing stable positioning services. In Fig. 14   improve 3D RMS by only 13.1% not as prominent as the static solutions.

Conclusions
To achieve the fast and precise single station positioning service, the PPP AR technique is the key to speed up the positioning bias convergence and improve the accuracy. Therefore, the FCB product is as essential as satellite orbit and clock products for PPP processing. We first introduced the FCB estimation method with the uncombined PPP model and the PPP AR strategy. The precision of the FCBs products and the PPP AR solutions with GPS/BDS-2/Galileo are evaluated. The uncombined PPP method is used to estimate the raw ambiguities on each frequency which are for the FCBs estimation. Although we do not have to sequentially estimate the WL, NL FCBs as in the IF PPP model, the integer coefficients are used to reduce the correlation between the first two raw ambiguities in the uncombined PPP model. The more stable FCBs on the new WL and NL combinations which are formulated with integer vectors, (1, − 1) and (4, − 3), are used in the FCBs estimator. The raw FCBs on each frequency recovered from the new WL and NL combinations are provided to users to flexibly achieve the ambiguity-fixed solutions with a specific PPP model.
With the proposed method, we estimated the precise FCB products based on iGMAS satellite products. For the posterior residuals of the FCB estimation, the RMS of the WL combination is significantly smaller than the NL combination. The longer wavelength of the WL combination is helpful to resist errors and keep a high consistency among the FCB measurements. The accuracy of BDS FCBs is degraded by the lower accuracy of satellite precise products. This will be improved in the BeiDou-3 Navigation Satellite System (BDS-3).
After estimating the FCB products, the PPP AR with GPS, BDS-2, and Galileo raw observations are processed in static and kinematic modes. In daily static positioning, the ambiguity-fixed solutions have not significant improvements compared with the ambiguity-float solutions, while the stability and reliability of positioning are enhanced with the PPP AR which is also found in kinematic mode. With hourly static solutions, 42.6% improvement and 97.7% ambiguity success fixing rate are achieved for the ambiguity-fixed solutions compared with the float solutions. The significant improvement of faster convergence time of positioning is also found. Finally, the 3D RMS for the ambiguity-fixed solutions is 2.10 cm.
Currently, only the BDS-2 data is adopted in our estimation which will be improved with the available BDS-3 satellites in the next work. The uncombined PPP model is flexible to process the triple-, quadruple-frequency or more frequency observables. Hence, the multi-frequency FCB products will be generated in future work.