1. Introduction
The BeiDou Global Navigation Satellite System (BDS-3) of China, consisting of 3 geostationary earth orbit (GEO), 3 inclined geostationary orbit (IGSO) and 24 medium earth orbit (MEO) satellites, started its service in 2020 [
1,
2,
3]. BDS-3 transmits a total of five frequency signals, including new signals, B1C (1575.420 MHz), B2a (1176.450 MHz) and B2b (1207.140 MHz), based on compatibility with the old signals, B1I (1561.098 MHz) and B3I (1268.520 MHz) [
4,
5]. The new signals, B1C and B2a, are divided into pilot and data components [
6,
7]. Some of the receivers in the International GNSS (global navigation satellite system) Service (IGS) network exclusively use the pilot component to generate pseudorange and carrier phase observations, and others use the pilot and data components jointly. Almost no receivers use the data component directly [
8]. Undoubtedly, the new frequencies together with the old frequencies enable users to possess more choices in data processing (e.g., cycle slip detection) and linear signal combinations, with more distinctive features (e.g., lower observation noise) as well [
9].
Slight time delays in the satellite signals during the transmission caused by the satellite or the receiver hardware, which are often referred to as hardware delay bias, exist in both code and carrier phase measurements and have different effects on different types of codes and carriers [
10]. To address the pseudorange hardware delay problem, a correction is usually performed on the user side using the timing group delay (TGD) parameter from either the real-time broadcast ephemeris or a post-processed high-precision differential code bias (DCB) product [
11]. The DCB is the difference between the hardware delay deviations of two code observations, which can be divided into inter-frequency and intra-frequency biases, depending on the frequency of the two code observations [
12]. Currently, the only two IGS analysis centers (IAC) that provide multi-GNSS DCB products are the Chinese Academy of Sciences (CAS) and the German Aerospace Center (DLR) [
13]. The estimation strategy of the DCB product of CAS and DLR are different; the former generates DCB products by single-site ionospheric total electron content (TEC) modeling through the IGGDCB method proposed by Li et al. [
14,
15], while the latter estimates DCB products after calculating TEC as ionospheric compensation through global ionospheric maps (GIM) [
16]. The rapid CAS product lags by approximately two to three days while DLR products are updated once every three months [
8].
The estimation and correcting methods of the DCB are hot research topics in geodesy since the DCB is one of the critical errors in GNSS precise positioning and ionospheric modeling [
17]. Guo et al. [
18] demonstrated the equivalence between BDS-2 TGD and DCB, derived a TGD/DCB correction model and performed standard point positioning (SPP) and precise point positioning (PPP) with and without DCB correction experiments using navigation messages and observations provided by the Multi-GNSS Experiment (MGEX). Ge et al. [
19] summarized the TGD correction model for each GNSS and performed SPP and PPP experiments. Before the construction of the BDS-3 system, the experimental satellite constellation (BDS-3e), with two IGSO and three MEO satellites, was established at the end of 2016 [
20]. Li et al. [
21] estimated the BDS-3e DCB and analyzed its performance by using the data provided by the MGEX and International GNSS Monitoring and Assessment System (iGMAS). After the fundamental constellation of BDS-3 started its service in 2018, Dai et al. [
22] first summarized the TGD correction of the new BDS-3 frequencies and evaluated the performance of the single- and dual-frequency SPP of BDS-2, BDS-3 and BDS-2/BDS-3, respectively. Subsequently, the DCB correction of the new frequencies of BDS-3 was analyzed and the effect of DCB correction on BDS-3 single- and dual-frequency positioning was tested using iGMAS data [
23]. Due to the broadcast of new frequencies, triple-frequency, and even multifrequency positioning, have received wider attention [
24,
25,
26,
27]. However, the TGD/DCB correction model for the triple-frequency positioning of BDS-3 has not been summarized in detail, and the effect of TGD/DCB correction on the effectiveness and positioning performance of triple-frequency positioning is unclear. This study derived the BDS-3 triple-frequency TGD/DCB correction models and aims to analyze the effect of TGD/DCB correction on BDS-3 triple-frequency SPP and PPP. Moreover, the effects of DCB on pseudorange noise, inter-frequency bias (IFB), ionospheric delay, tropospheric delay and floating ambiguity are discussed.
The rest of this paper is organized as follows. After the introduction, the BDS-3 triple-frequency positioning model is first summarized in
Section 2, and the triple-frequency TGD/DCB correction formulas are derived for different frequency combinations, including triple-frequency ionospheric free (IF) combination (IF123), two dual-frequency IF combinations (IF1213) and triple-frequency uncombined (UC123). Then, the stability of the DCB product released by the CAS is analyzed. Finally, the processing strategy is presented.
Section 3 verifies the effectiveness of TGD/DCB correction and its effect on SPP and PPP. Discussion and conclusions are provided in
Section 4 and
Section 5, respectively.
4. Discussion
Four reasonable triple-frequency combinations are selected for the experiments, namely, B1I/B2a/B3I, B1C/B2a/B3I, B1I/B2b/B3I and B1C/B2a/B3I.
The TGD/DCB correction significantly improves the positioning accuracy of SPP, especially for the IF123 and IF1213 models, because the errors are amplified, including the UCDs by the signal combinations. When the satellite TGD/DCB is corrected, the positioning accuracies of UC123 in the E and U components reaches nearly the same as those of the IF123 and IF1213 modes but slightly worse in the N component.
The DCB correction improves the accuracy of the initial epochs of PPP, but the improvement effect gradually diminishes with time. Since the high-precision carrier phase contributes more and more in subsequent epochs, the positioning accuracies of the four combinations of IF123, IF1213 and UC123 models are nearly the same after the filtering convergence. However, the convergence time of UC123 is longer than that of IF123 and IF1213 because there are more parameters to be estimated.
DCB not only affects the PPP convergence time but also affects other parameters. The pseudorange residuals without DCB have a systematic bias of several meters, and the bias value varies quite differently for different satellites. When the DCB is corrected, the pseudorange residuals are significantly more reasonably distributed, approaching zero-mean random noise.
The IFB needs to be estimated for IF1213 and UC123 due to satellite and receiver UCDs. When the satellite DCB is corrected, the IFB fluctuates less during the day and has short-term stability. The DCB also affected the estimated STEC. If DCB is not taken into account, the estimated STEC will show anomalies or even negative numbers. The floating ambiguity of each satellite is affected differently by the DCB at the earliest epoch but converges to nearly the same value after some time. Tropospheric delay is nearly unaffected by DCB, and the differences between the estimated tropospheric parameters with or without DCB correction is less than 5 mm after convergence.
In summary, the DCB correction significantly improves the positioning accuracy of BDS-3 triple-frequency SPP and reduces the convergence time of PPP. Therefore, the correction should not be neglected in BDS-3 triple-frequency positioning.