1. Introduction
The modernization of global navigation satellite systems (GNSS) presents both opportunities and challenges. A notable aspect is the enhanced availability of observable frequencies, which substantially increases the observations available to users [
1]. This development also facilitates the provision of better resolutions in a range of fields, including wide-area precise positioning, low-orbit satellite orbiting, ionospheric monitoring, earthquake and tsunami monitoring and warning systems, and time transfer, among others. It is noteworthy that the satellite clock offset products provided by the International GNSS Service (IGS) are currently dependent on legacy dual-frequency ionosphere-free (
IF) combinations (e.g., GPS L1/L2, BDS B1/B3 with GalileoE1/E5a). Careful evaluation of the impact of additional frequencies on positioning models is a fundamental requirement in constructing a multi-frequency observation model.
Employing multi-frequency observations in precise point positioning (PPP) can offer various advantages, but also poses several challenges. Montenbruck et al. initially observed deviations in carrier phase observations that were time-, signal- and satellite-dependent, such as the discrepancy between GPS L1/L2 and L1/L5 combined carrier phase
IF observations [
2]. Meanwhile, Pan et al. showed that the time-varying phase bias is assimilated by the satellite clock offset parameter in the precise clock estimation. Consequently, when varying observations are employed, the satellite clock estimates will differ, and the distinction between them is referred to as the inter-frequency clock bias (
IFCB) [
3]. The error mentioned above could cause significant discrepancies of several tens of centimeters in the multi-GNSS multi-frequency PPP observation model [
4]. The
IFCB compensation method is a key issue in multi-GNSS multi-frequency PPP data processing. In general, there are three approaches to address the
IFCB problem: The first is to determine the satellite clock offset through the addition of uncombined L5 observations to the conventional dual-frequency
IF model [
5], so that two sets of satellite clock offset products are available, e.g., the L5 uncombined satellite clock offset and conventional L1/L2
IF satellite clock offsets, where the L1/L5
IF satellite clock offset can be obtained through the conversion of the relationship. The second is based on the combined geometry-free and ionosphere-free (
GFIF) carrier phase observations, where the
IFCB of each station is obtained using the epoch difference strategy, and the satellite
IFCB can be obtained through averaging the
IFCB of all stations [
6,
7,
8,
9]. The last, but seldom used, option is to ignore the
IFCB. Treating the
IFCB as a parameter in observation is less. Currently, the second method is widely used for
IFCB correction at the third frequency, but the receiver
IFCB is also absorbed by the satellite
IFCB. This requires substantial computational resources, while the reliability and accuracy of
IFCB estimation are limited by the number of stations in the reference network, which also increases the probability of errors in the
IFCB calculation results [
10,
11,
12,
13]. However, post and real-time phase
IFCB products are available from few organizations, which has clearly hampered the development of the multi-GNSS multi-frequency PPP. This limits the real-time applicability of the approach, especially for users who require immediate and continuous positioning solutions.
With the arrival of the multi-GNSS multi-frequency era,
IFCB has also become a hot spots of current research. Many scholars have studied and analyzed the impact of
IFCB in GPS/BDS/Galileo triple frequency (TF) PPP. Montenbruck and Lin et al. studied the
IFCB of GPS satellites using GPS TF observations, and showed that
IFCB is a GPS bias with periodic and time-varying characteristics [
14,
15,
16,
17]. Chen and Montenbruck et al. showed that the
IFCB amplitude varies up to 10–40 cm for GPS (L1/L2 and L1/L5), within 1–2 cm for Galileo (E1/E5a and E1/E5b), and within 1–3 cm for BDS-3 (B1/B3 and B1/B2a) [
2,
18]. Zhang found that the GPS Block IIF satellites’
IFCB showed a variation period of more than two years, which is apparently consistent for two satellites placed in the same orbital plane [
19]. Some scholars have also investigated the temporal characteristics of
IFCB in depth, and have built different prediction models for
IFCB correction [
20]. For the current study, it is clear that
IFCB is present at the L5 frequency of GPS Block IIF satellites, while Galileo and BDS-3 are barely affected [
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20]. An in-depth analysis of the additional signals is necessary to take full advantage of the benefits that additional frequencies bring. However, most of the studies in the literature tackle the characteristics of
IFCB and its prediction model based on datasets of a few, or dozens, of stations [
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22], and there are no relevant unified criteria to evaluate the accuracy of the products. Additionally, the accuracy of the
IFCB models and predictions may degrade over time; the prediction accuracy of
IFCB models decreases when the time span is increased from a day to a week. This means that the accuracy of
IFCB correction values may decrease over longer periods of time, potentially impacting the positioning accuracy. Most studies are limited to analyzing the temporal characteristics and periodicity of
IFCB time series, but fewer studies consider
IFCB as a parameter to be estimated and analyze the impact of stochastic models of
IFCB parameters on the positioning performance of the multi-GNSS multi-frequency PPP. Treating
IFCB as a parameter undoubtedly increases the computational load. Nonetheless, with the progression of scientific and technological advancements, this additional load does not pose a formidable challenge for future modern processors.
Most studies encountered in the literature focus primarily on the estimation and correction of
IFCB for triple-frequency PPP. The stochastic modeling of
IFCB plays a significant role in improving the understanding and impact of its parameters. The techniques proposed in this paper can be implemented directly in other positioning methods or scenarios that involve triple-frequency observations. In this study, an analysis of the power spectral density relevant to
IFCB is conducted and product correction methods are compared with parameter estimation methods. A comprehensive analysis is carried out to assess the impact (e.g., positioning accuracy, convergence time, and phase residuals) of stochastic
IFCB modeling on PPP performance based on the International GNSS Service’s (IGS) Multi-GNSS Experiment (MGEX) one-week precise-orbit and clock-offset products (Wuhan University) and data from 130 tracking stations. Additionally, the stochastic models which are employed to model the
IFCB parameters in the multi-GNSS multi-frequency PPP are explicated using the Undifferenced and Uncombined (UDUC) PPP model [
23,
24,
25].
The remainder of the article is structured as follows; the MGEX UDUC PPP model is first elaborated in
Section 2, taking into account
IFCB. In
Section 3, experiments and data processing strategies are presented. The optimal power spectral density of the
IFCB and the impact of different
IFCB stochastic modeling on PPP performance is demonstrated and compared, followed by a discussion, in
Section 4. The concluding
Section 5 contains a summary of the findings.
4. Experimental Validation
The previous theoretical analysis indicates that
IFCBs stem from both the satellite- and receiver-dependent phase time-varying parts of the hardware delay, which is the difference between the satellite clock offsets for different
IF combinations [
25]. In this study, global data from 105 stations were collected for
IFCB product estimation through methods considering different
IF combinations. Subsequently, 25 stations were used to analyze the impact of
IFCB stochastic modeling on the multi-GNSS multi-frequency PPP performance for simulated kinematic PPP. Based on the different
IFCB processes, five schemes were implemented: (a) ignoring the effect of GPS third frequency
IFCB, (b) using
IFCB product correction, (c) using random constants to estimate the
IFCB, (d) using a white noise process to estimate the
IFCB, and (e) using a random walk process to estimate the
IFCB. Then, the third-frequency residual was analyzed using the five schemes. The results are summarized and discussed after the analysis of the findings.
In this study, one-week observation data from 25 stations were used to evaluate PPP performance. The kinematic PPP performance was evaluated using indices such as positioning accuracy, convergence time, the third frequency phase residuals’ root mean square (RMS) and estimated
IFCB. The position was considered to have converged only when the positioning error (east, north, up) of the epoch and the subsequent 20 epochs were within a 1 dm threshold [
34,
35]. For convenience, the five schemes mentioned above are marked as
IFCB_NO (Neglect),
IFCB_PC (Product Correction),
IFCB_CT (random constant),
IFCB_WN (white noise), and
IFCB_RW (random walk), respectively. For the
IFCB_RW scheme, the
IFCB parameters were modeled as a random walk with a power density of 0.6 m/sqrt(s). The reason for using this value is explained in
Section 4.1.2; for the first appearance of the satellite in the arc segment, the initial value of the
IFCB was set to 0.01 m. It should be noted that the
IFCB amplitudes of BDS-3 and Galileo are less than 2.0 cm, with negligible effects on multi-GNSS multi-frequency PPP positioning [
4,
14]; therefore, the effect of the third-frequency
IFCB was ignored in this study. To facilitate the description, the
IFCBs discussed in the following sections refer to the GPS L5
IFCB.
4.1. IFCB Estimation
4.1.1. IF Combinations Estimation (IFCB_PC)
Figure 2 shows the GPS L5
IFCB products obtained using the IF combinations estimation, where different colors represent the
IFCBs of different satellites. From
Figure 2, the single-day amplitude of GPS Block IIF satellite L5
IFCB is larger among all satellites and lies within the range of 10–20 cm, which is obviously a PPP error that cannot be ignored. Additionally, it has been demonstrated that the amplitude of GPS Block III L5
IFCB ranges from 1 to 3 cm within a single day, and the standard deviation of
IFCB is about 1.5 mm in a single epoch, which is almost independent of
IFCB [
4]. Therefore, it is necessary to focus only on the effect of the GPS Block IIF satellites’ L5
IFCB and to analyze the corresponding effect on the multi-GNSS multi-frequency PPP positioning (see
Section 4.2). Since
IFCB is considered a satellite- and receiver-related hardware delay between phase frequencies, which is caused by the difference between satellite clock offset calculated with different IF combinations, a thorough analysis of the satellite’s internal and external configuration is required to fully understand the causes of its generation. The focus of this study is on the impact of
IFCB on multi-GNSS multi-frequency PPP positioning performance, so the causes of
IFCB are not discussed in this paper.
4.1.2. Parameter Estimation (IFCB_CT, IFCB_WN, IFCB_RW)
In order to analyze the effect of different stochastic models on
IFCB results, in this section the
IFCB is considered as a parameter and estimated lumped with other parameters, such as ambiguity, atmospheric error, and receiver clock offsets; corresponding time series are plotted. Taking BRUX station as an example, the
IFCB daily series (DOY 18 January 2022) derived from the GPS/BDS-3 multi-GNSS multi-frequency PPP solutions are given in
Figure 3, with different colors representing the
IFCB of different satellites.
Figure 3 shows the
IFCB results estimated via GPS/BDS-3 multi-frequency PPP using the
IFCB_CT,
IFCB_WN, and
IFCB_RW schemes. It can be seen from
Table 2 that the
IFCB differences among the three
IFCB process schemes (
IFCB_PC,
IFCB_WN,
IFCB_RW) were less than 0.02 m for BRUX. Nonetheless, the
IFCB differences between
IFCB_CT and other schemes (
IFCB_PC,
IFCB_WN, and
IFCB_RW) were significant. Comparing
Figure 2 and
Figure 3, in the
IFCB_CT scheme, the estimation of
IFCB as a random constant does not fit the characteristics of
IFCB and cannot truly reflect the change of its temporal characteristics, thus affecting the positioning accuracy and convergence performance of the multi-GNSS multi-frequency PPP (see
Section 4.2); therefore, the
IFCB_CT scheme cannot be used to model
IFCB variation.
The difference in satellite
IFCB amplitudes between
IFCB_PC and
IFCB_WN was small, with an average difference of <1 cm, and the temporal characteristic change between them was similar, with an STD of <5 mm. Obviously, considering
IFCB as a white noise process is consistent with the characteristic change of
IFCB, which reflects the temporal characteristic change of
IFCB more accurately and thus improves the multi-frequency PPP positioning accuracy and convergence. Therefore, the white noise process is suitable for describing the temporal characteristic change of
IFCB. The same conclusion applies to the
IFCB_RW scheme. It is worth noting that there are no studies on
IFCB’s power spectral density (
σIFCB) setting, so in order to obtain the optimal power spectral density value, the tests considered the range
σIFCB [60, 6 × 10
−5] at the 25 stations (DOY 18 January 2022). The corresponding statistical results are shown in
Table 3. As shown in
Table 3, the convergence times in the
σIFCB range within [60, 6 × 10
−3] are equal, and the
IFCB_PC results show that the inter-epoch variation and amplitude of
IFCB are much smaller than 0.6 m, so
σIFCB = 0.6 m/sqrt(s) is obviously a suitable choice for the
IFCB_RW scheme. Consequently, this value was adopted chosen in this paper and the same value was used for the
IFCB_WN scheme. The
IFCB_WN scheme uses the same processing scheme as
IFCB_RW. In conclusion, comparing
Figure 2 and
Figure 3 and the results in
Table 2,
IFCB_WN and
IFCB_RW reflect the variation of the temporal characteristics of
IFCB accurately, which indicates that the processing methods of the schemes used in this study are reasonable.
4.2. Three-Frequency PPP Performance
4.2.1. GPS-Only
The GPS-only results are given in
Figure 4, which shows the convergence time box diagram for five
IFCB process schemes. More than 24% of
IFCB_CT schemes had a convergence time over 60 min; this performance was comparable with that of
IFCB_NO. The average convergence times of the
IFCB_RW and
IFCB_WN schemes were comparable with that of
IFCB_PC, with the difference being close to 1 min. The median convergence time improvement for the
IFCB_RW and
IFCB_WN schemes, compared to
IFCB_NO was 43.0 min and 26.0 min, respectively. It is noteworthy that the convergence performance of GPS-Only PPP solutions was worst if the
IFCB_NO or
IFCB_CT schemes were adopted. This result was mostly caused by
IFCB (see RMS values in
Table 4).
Likewise, a clear difference in terms of weekly RMS was observed between the
IFCB_NO scheme and the other four schemes (
Table 4). The three
IFCB process schemes (
IFCB_PC,
IFCB_WN, and
IFCB_RW) were comparable in terms of positioning accuracy. Furthermore, the
IFCB_NO and
IFCB_CT schemes were generally worse than the other schemes.
4.2.2. GPS/BDS-3 and GPS/Galileo
Figure 5 shows the convergence time box diagram of GPS/BDS-3 and GPS/Galileo PPP for the five
IFCB process schemes. Similar to GPS-only, the convergence performance of the
IFCB_PC,
IFCB_RW, and
IFCB_WN schemes (about 15 min on average for GPS/BDS-3) was much shorter than that of the
IFCB_NO and
IFCB_CT schemes (about 20 min on average for GPS/BDS-3), which is a striking improvement. In general, the convergence performance of the
IFCB_PC,
IFCB_RW, and
IFCB_WN schemes is comparable. From
Figure 6, it is evident that the GPS/Galileo convergence performance is similar to that of the GPS/BDS-3 PPP solutions.
Figure 6 shows the comparisons of kinematic positioning errors for BRUX station for the five
IFCB process schemes.
Table 5 shows the weekly positioning accuracy of GPS/BDS-3 PPP at the 25 stations for the five
IFCB schemes. Not surprisingly, the positioning accuracy of the
IFCB_NO scheme is worse than that of the other schemes. Comparing
IFCB_NO with the other three schemes (
IFCB_PC,
IFCB_WN, and
IFCB_RW), the positioning accuracy is improved by 50.72% (from 2.05 cm to 1.01 cm) in the east, by 39.99% (1.55 cm to 0.93 cm) in the north, and by 52.72% (4.03 cm to 1.90 cm) in the up components. As shown in
Table 5, the positioning accuracy of the multi-GNSS multi-frequency PPP solutions with
IFCB_PC,
IFCB_WN, and
IFCB_RW is comparable among the five
IFCB process schemes.
With regard to convergence as GPS/Galileo PPP solutions, one can see from
Table 5 that the positioning accuracy of the
IFCB_PC,
IFCB_RW, and
IFCB_WN schemes is improved by 58.53% (from 2.26 cm to 0.94 cm) in the east, 45.33% (1.64 cm to 0.89 cm) in the north, and by 45.00% (3.74 cm to 2.05 cm) in the up components compared with the
IFCB_NO scheme. As shown in
Table 5, the
IFCB_PC,
IFCB_WN, and
IFCB_RW schemes are comparable among the five
IFCB process schemes for positioning accuracy of PPP solutions.
4.2.3. GPS/BDS-3/Galileo
Figure 7 shows the convergence time box diagram of GPS/BDS-3/Galileo PPP for the five
IFCB process schemes. It is evident that the differences in convergence performance are similar to those of the GPS/BDS-3 PPP solutions.
Comparing
IFCB_PC,
IFCB_WN, and
IFCB_RW convergence times of GPS/BDS-3/Galileo PPP solutions, the average convergence time of
IFCB_NO (or
IFCB_CT) is 5 min less than that of
IFCB_WN (and
IFCB_PC or
IFCB_RW). It is also evident, from
Figure 7, that the median convergence times are comparable among the five
IFCB process schemes. This is attributed to the increase in the number of BDS-3 and Galileo observations, in which the strength of the observation equation and the redundancy are increased. The positioning accuracy of the
IFCB_PC,
IFCB_RW, and
IFCB_WN schemes is improved by 41.34% (from 1.37 cm to 0.81 cm) in the east, by 23.29% (1.10 cm to 0.84 cm) in the north, and by 37.92% (2.72 cm to 1.70 cm) in the up components compared to the
IFCB_NO scheme. As shown in
Table 6, the
IFCB_PC,
IFCB_WN, and
IFCB_RW schemes are comparable among the five
IFCB process schemes for positioning accuracy of PPP solutions.
In conclusion, the above three sets of experiments show that it is feasible to use the parameter estimation method to deal with IFCB, and its effect on positioning accuracy and convergence performance is comparable to that of using product correction. Therefore, the white noise or random walk process can be used to describe IFCB parameters to accelerate the multi-GNSS multi-frequency PPP convergence in the absence of product correction or real-time applications.
4.3. Residual Analysis
Some unmodeled errors, such as
IFCB, are reflected in the post residuals due to poor observation models or incomplete corrections.
Figure 8 shows the post-phase residuals of the five
IFCB process schemes of GPS L5 at BRUX station (DOY 18 January 2022), where different colors represent the residuals of different satellites. As previously concluded, the
IFCB_NO and
IFCB_CT schemes exhibit significant systematic bias effects in the GPS L5 phase residuals, while the
IFCB_PC,
IFCB_WN, and
IFCB_RW schemes eliminate these effects. At the same time, to further evaluate the effect of
IFCB on GPS L5 phase residuals,
Figure 9 and
Table 7 give the weekly RMS values of the GPS L5 phase residuals for the 25 stations. The histogram in
Figure 9 shows the average GPS L5 phase residuals for one week at the 25 stations, while the line graph shows the percentage improvement of each station relative to the
IFCB_NO scheme.
Table 7 shows the average GPS L5 phase residual RMS values over one week of the five process schemes at the 25 stations. Compared to
IFCB_NO, the RMS of the GPS L5 phase residuals of
IFCB_PC,
IFCB_WN, and
IFCB_RW were improved from 0.3867 cm to 0.2311 cm, 0.0971 cm, and 0.0980 cm, respectively, which represent improvements of 40.25%, 74.89%, and 74.66%. The effect is especially pronounced in the case of the
IFCB_WN and
IFCB_RW schemes. Accordingly, it is concluded that the effect of the apparent systematic bias in the GPS L5 phase residuals can be eliminated using the
IFCB_PC,
IFCB_WN, or
IFCB_RW schemes; the latter two schemes have significantly better residual values than the
IFCB_PC scheme. The above conclusions indicate that the use of the estimation method results in an improved correction of the effects of
IFCB compared to the use of product correction, and can help avoid satellite false rejects due to large residuals during data processing.
4.4. Discussion
The above results show that ignoring the effect of
IFCB can seriously degrade the positioning performance of multi-GNSS multi-frequency PPP, and therefore the effect of
IFCB must be considered. The IF combinations’
IFCB estimation method averages the
IFCB over all stations, and although this method uses epoch differences in the calculation process to eliminate most of the errors, there are still some unmodeled errors remaining in the epoch-difference
GFIF observations, such as multipath, higher-order ionospheric delay, and others. This method does not consider the receiver
IFCB, so the average receiver
IFCB is absorbed by the satellite
IFCB. If the number of stations in the network is limited, the receiver
IFCB cannot be eliminated through averaging, so this method is not theoretically rigorous [
25]. This is also the reason why the residuals of
IFCB_PC are larger than those of
IFCB_RW or
IFCB_WN. This effect can be derived from the following equation:
Taking into account the unmodeled error
in Equation (10), the following equation can be obtained:
Then, the new epoch-difference
GFIF is expressed as follows:
Therefore, the current epoch satellite
IFCB becomes:
Obviously
causes an accumulation of errors, indicating that the
IFCB estimated using the
IFCB_PC method contains an unmodeled error
for
n stations, but this unmodeled error is small (RMS < 2 mm, see RMS values in
Table 7) and can be neglected. The analysis in the previous sections shows that
IFCB_NO and
IFCB_CT have comparable convergence performance, but the positioning accuracy of the latter is improved.
IFCB_PC has comparable positioning accuracy and convergence time with
IFCB_WN and
IFCB_RW, but the residuals of the latter two are smaller than that of the former. This indicates that the use of a proper
IFCB stochastic model can achieve comparable performance with the use of the product, but the method of
IFCB parameter estimation increases the computational effort on the user side and incurs higher requirements for user equipment. However, the impact on positioning performance is minimal. Therefore, the
IFCB_RW or
IFCB_WN methods are recommended for processing
IFCB, when an
IFCB product is unavailable, in the case of users requiring real-time applications.
5. Conclusions
The incorporation of multi-frequency signals into GNSS has introduced new possibilities for precise positioning and rapid ambiguity resolution. The IFCB pertains to the variation in clock offsets among distinct frequencies of satellite signals, and its proper processing is of great importance, as it significantly impacts PPP accuracy when using multi-frequency signals. Neglecting the IFCB can result in reduced positioning accuracy and increased convergence time in both static and kinematic PPP solutions. In this study, the focus lies on the appropriate modeling for phase IFCB in multi-GNSS multi-frequency PPP based on the UDUC observation model. The optimal power spectral density applicable to IFCB is analyzed and studied, and the product correction methods are compared with parameter estimation methods. The impact of IFCB stochastic modeling on the performance of Undifferenced and Uncombined PPP is also evaluated from the aspects of convergence time and positioning accuracy through five IFCB process schemes, namely IFCB_NO, IFCB_PC, IFCB_CT, IFCB_RW, and IFCB_WN.
To validate these findings, an analysis of the effect of IFCB stochastic modeling on UDUC the multi-GNSS multi-frequency PPP performance was conducted using weekly data (DOY 18–24 January 2022) from 130 IGS MGEX tracking stations. Through this analysis, the following conclusions can be stated: (i) the optimal power spectral density for IFCB is 60, 0.006 m/sqrt(s); (ii) modeling IFCB as a random walk is feasible; (iii) the multi-GNSS multi-frequency PPP performance is comparable and best for the IFCB_PC, IFCB_RW, and IFCB_WN schemes, while the performance of the IFCB_NO and IFCB_CT schemes is poor; (iv) the third frequency phase residuals of the IFCB_RW and IFCB_WN scheme are smaller than those of the IFCB_PC scheme. However, any potential impacts on the troposphere, ionosphere, and float ambiguity are not discussed.
The methodology for estimating multi-GNSS satellite phase products is currently in a developmental stage, thus these products are not yet fully mature. The commonly utilized approach, which involves the use of IFCBs corrected via products, is applicable exclusively when the server provides such products. However, IFCB products are not readily accessible, especially for users who require immediate and continuous positioning solutions. Essentially, a random walk is an expansion and accumulation of white noise, and the IFCB_RW and IFCB_WN schemes are indistinguishable. However, compared to white noise, the random walk has the advantages of exhibiting trends, correlation, and memory, making it more suitable for modeling phenomena with directionality, time dependence, and historical dependence. Therefore, in the absence of IFCB products, modeling IFCBs as random walks can result in significantly more reliable performance.