Next Article in Journal
Position Accuracy Improvement by Implementing the DGNSS-CP Algorithm in Smartphones
Previous Article in Journal
Biosensing by WGM Microspherical Resonators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mitigating BeiDou Satellite-Induced Code Bias: Taking into Account the Stochastic Model of Corrections

1
School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China
2
Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan 430079, China
3
Collaborative Innovation Center for Geospatial Technology, Wuhan 430079, China
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(6), 909; https://doi.org/10.3390/s16060909
Submission received: 8 May 2016 / Revised: 6 June 2016 / Accepted: 15 June 2016 / Published: 18 June 2016
(This article belongs to the Section Remote Sensors)

Abstract

:
The BeiDou satellite-induced code biases have been confirmed to be orbit type-, frequency-, and elevation-dependent. Such code-phase divergences (code bias variations) severely affect absolute precise applications which use code measurements. To reduce their adverse effects, an improved correction model is proposed in this paper. Different from the model proposed by Wanninger and Beer (2015), more datasets (a time span of almost two years) were used to produce the correction values. More importantly, the stochastic information, i.e., the precision indexes, were given together with correction values in the improved model. However, only correction values were given while the precision indexes were completely missing in the traditional model. With the improved correction model, users may have a better understanding of their corrections, especially the uncertainty of corrections. Thus, it is helpful for refining the stochastic model of code observations. Validation tests in precise point positioning (PPP) reveal that a proper stochastic model is critical. The actual precision of the corrected code observations can be reflected in a more objective manner if the stochastic model of the corrections is taken into account. As a consequence, PPP solutions with the improved model outperforms the traditional one in terms of positioning accuracy, as well as convergence speed. In addition, the Melbourne-Wübbena (MW) combination which serves for ambiguity fixing were verified as well. The uncorrected MW values show strong systematic variations with an amplitude of half a wide-lane cycle, which prevents precise ambiguity determination and successful ambiguity resolution. After application of the code bias correction models, the systematic variations can be greatly removed, and the resulting wide lane ambiguities are more likely to be fixed. Moreover, the code residuals show more reasonable distributions after code bias corrections with either the traditional or the improved model.

1. Introduction

The Chinese BeiDou navigation satellite system (abbreviated as BDS, or BeiDou, or COMPASS) has already launched a regional navigation service by the end of 2012, and continues to develop a global system in the near future. As of October 2015, a total of 20 BeiDou satellites including six geostationary orbit (GEO), seven inclined geosynchronous orbit (IGSO), and seven medium-altitude Earth orbit (MEO) satellites have been launched, transmitting triple-frequency signals centered at B1 (1561.098 MHz), B2 (1207.14 MHz), and B3 (1268.52 MHz). However, it is worth mentioning that only 13 operational BeiDou satellites are presently available (5GEO + 5IGSO + 3MEO) for public use.
With the advent of BeiDou, a large amount of research has been dedicated to the aspects of precise orbit and clock determination [1,2,3], regional ionospheric and tropospheric delay modeling [4,5,6], differential code bias and timing group delay corrections [7,8,9,10], multipath effects analysis [11,12], triple-frequency ambiguity resolution [13,14], and precise positioning [1,15,16,17,18,19], etc. When it comes to the characteristics of BeiDou signals, Nadarajah et al. first discovered the BeiDou Inter-Satellite-Type Biases (ISTBs), and deeply investigated the impacts of such biases on real-time kinematic (RTK) positioning and mixed receiver attitude determination [20,21]. A recent study revealed that the code-phase divergences, which are absent for GPS, GLONASS, and Galileo, however, are commonly found in BeiDou IGSO and MEO satellites [22]. The systematic variation was first mentioned by Hauschild et al. [11], but the origin of the systematic variation was not identified. The same code bias variation was also detected by Montenbruck et al. through the analysis of multipath combination [23,24]. They found that the systematic bias is elevation-dependent, which varies by 0.4–0.6 m from horizon to zenith, particularly pronounced for the B1 signal. Since such systematic bias was consistently observed with other receivers and antennas it was, therefore, attributed to the transmitting satellites. More characteristics of the systematic bias variations were well discussed in Wanninger and Beer [22]. They identified significant differences between two groups of satellites (i.e., IGSO and MEO) and among carrier frequencies (B1, B2, and B3). A correction model was then proposed to accommodate the satellite induced code bias variations, which is definitely an important contribution to the current BeiDou constellation. Nevertheless, the data used for modeling these variations is very limited. Observations from two 10-day intervals in March 2014 (DoY 070–079) and in June 2014 (DoY 160–169) were exclusively used. The long-term stability of the systematic biases could not be identified conclusively from such limited dataset. Consequently, the effectiveness of the mitigating model with fixed correction values may be a debatable issue. In addition, only the correction values on the specific nodes are provided, whereas the stochastic information, i.e., the precision indexes, are completely missing in their correction model. Such precision indexes are, however, important for the stochastic model establishment in precise applications, such as precise point positioning (PPP). Ignoring the stochastic error of the corrections will lead to an overestimation of the actual precision of code observations. This is particularly pronounced in the case of inaccurate corrections are applied.
This paper deals with the same systematic bias variations as earlier described by Hauschild et al. [11], Montenbruck et al. [23,24], and Wanninger and Beer [22]. However, we use more datasets with a time span of almost two years (from January 2014 to October 2015) to model such orbit type-, frequency-, and elevation-dependent code biases. Moreover, the stochastic model of the corrections is considered as well. Following an analysis of the code bias variations and a brief introduction of the observation datasets used, we develop an improved correction model, which produce both correction values and their precision indexes. Afterwards, effectiveness validations for the proposed model are presented. Finally, we end with the conclusions.

2. Code Bias Variation Analysis

Providing dual- or triple-frequency GNSS observations, multipath (MP) effects are commonly assessed by a linear combination of code and carrier phase observables [12,25]. That is:
M P i ( i , j ) = P i f i 2 + f j 2 f i 2 f j 2 λ i φ i + 2 f j 2 f i 2 f j 2 λ j φ j
where M P i is the multipath combination on frequency i ( i , j = 1 , 2 , 3 , i j ) , P and φ represent code range and carrier phase observables, f and λ are, respectively, the frequency and wavelength of the specific signal. If the phase noise and phase multipath are neglected, the MP combination mainly consists of code multipath, a constant ambiguity term which is a combination of the ambiguities of the two phase measurements, a combined signal delay, and thermal noise. Subtraction of the mean value from the measurements removes the phase ambiguities, which are constant if there are no cycle slips. The resulting time series of MP is dominated by code multipath and thermal noise.
The MP combination as shown in Equation (1) is an ionosphere-free and geometry-free combination. The linear coefficients are selected in such a way that ionospheric and tropospheric delay as well as all geometric contributions (e.g., clocks, orbits, and antenna movements, etc.) cancel out. Therefore, it can be used with both static and kinematic data. However, it is worth mentioning that the absolute values of multipath are unknown. Only variations of code multipath are available from the MP series.
For each code, there are at least two choices of MP combination using different frequencies among triple-frequency observations. Take BeiDou as an example; the code multipath on B1 (MP1) can be derived from Equation (1) using either (B1, B2) or (B1, B3) frequency combinations. However, one should be aware that the latter one is noisier than the former one due to a larger multiplication factor. Therefore, to keep a lower noise level, both MP1 and MP2 are derived from (B1, B2), and MP3 is computed from (B1, B3) combination in this paper.
Shown in Figure 1 are the MP series of BeiDou satellites observed on JFNG (Wuhan, China) with a Trimble Net R9 receiver. BeiDou satellites are grouped by three different orbit types, i.e., GEO, IGSO, and MEO. For comparison reason, the MP series of GPS satellites tracked by the same receiver are plotted as well in Figure 1. As expected, signals observed at low elevation angles are more likely to be affected by multipath error. On the contrary, a high elevation satellite is less likely to cause multipath error. However, unexpected code biases (drifts), which are absent for GPS, are found in BeiDou MEO and IGSO satellites. The code bias variation reaches approximately 1 m at high elevation angles, particularly for the MEOs on the B1 signal. Similar variation and its characteristics have been well discussed in Wanninger and Beer [22]. Since the code bias variation was consistently observed on other stations equipped with different receivers and antennas, such bias thus was attributed to the transmitting satellites. In addition, anomalies can be found in the MP3 series of GPS (i.e., code multipath effect on L5 signal), which show significantly larger values at high elevation angles. Similar code-phase divergences, which were attributed to satellite-internal signal reflections, had been earlier reported for GPS SVN49 and SVN62 in Spring and Dislssner [26]; Montenbruck et al. [27]; and Hauschild et al. 2012 [11].
To investigate the degree of linear dependence between MP and elevation, the Pearson correlation coefficient [28] was applied and computed as:
ρ M P , E = M P E M P E N ( M P 2 ( M P ) 2 N ) ( E 2 ( E ) 2 N )
where E and N represent satellite elevation angle and the number of samples, respectively. The correlation coefficient ρ M P , E ranges from −1 to 1. A value of 1 or −1 implies that a linear equation describes the positive or negative relationship between M P and E perfectly, with all data points lying on a line. A value of 0 implies that there is no linear correlation between the two variables.
Table 1 lists the corresponding correlation coefficients for both GPS and BeiDou satellites. In general, the correlation coefficient under 0.3 represents low correlation, while the coefficient between 0.3 and 0.5 indicates normal. If the correlation coefficient exceeding 0.5 but less than 0.7, a medium correlation exists between two components. The correlation coefficient between 0.7 and 0.9 represents a high correlation. Additionally, if the correlation coefficients were beyond 0.9, there would be a strong correlation between two components. The coefficients in Table 1 demonstrate a medium correlation between MP and elevation for BeiDou IGSOs and MEOs, whereas the correlations of GPS and BeiDou GEO satellites are rather weak or even uncorrelated. Nevertheless, one should keep in mind that the GEOs are visible at almost constant elevation angles within restricted region (i.e., Asia-Pacific region). It is unable to infer whether such code biases affect BeiDou GEO code measurements. Hence, the GEOs are excluded in the following context.

3. Observation Datasets

As of August 2015, the multi-GNSS experiment (MGEX) [29] offers a global network of approximately 70 stations of GNSS receivers capable to receive BeiDou signals. Most of these stations are equipped with Trimble NetR9 receivers providing triple-frequency BeiDou observations. However, only a few stations are equipped with other receiver types, such as Septentrio POLARX4/POLARX4TR/ASTERX3 and Leica GR25 receivers, providing dual-frequency observations on frequencies B1 and B2. To achieve a global coverage of receiving stations and to be able to compare different receiver types, six stations (CHPG at Cachoeira Paulista, Brazil, CUT0 at Perth, Australia, GMSD at Nakatane town, Japan, JFNG at Wuhan, China, NKLG at Libreville, Gabon, and REUN at Le Tampon, France) equipped with Trimble receivers and four stations (CEBR at Cebreros, Spain, KOUR at Kourou, French Guiana, MAL2 at Malindi, Kenya, and MGUE at Malargue, Argentina) equipped with Septentrio receivers were selected as shown in Figure 2. Data recorded on the first day of each month from January 2014 to October 2015 were used as the core datasets for this study. The datasets spanning almost two years facilitate analysis of a long-term variation of the BeiDou satellite induced code bias. However, it is worth mentioning that some of the data may be unavailable during this period due to station related reasons, e.g., new settled devices or upgrading issues, etc. Whatever the reason, an average of 13–14 days data is usable for each station.
Figure 3 shows the number of MP samples against elevation bins with an interval of 10° for BeiDou MEO and IGSO satellites. As already mentioned, stations equipped with Septentrio receivers provide only B1 and B2 signals and, thus, the number of MP3 for MEOs is much smaller than that of MP1 and MP2. As to IGSOs, the number of MP3 is comparable to that of MP1 and MP2 since the IGSO satellites are mostly tracked by triple-frequency Trimble receivers in Asia-pacific area. Furthermore, the number of MP samples is relative smaller at the lowest and highest elevation bins due to spatial configuration reason.

4. Improved Correction Model

As shown in Figure 1, the code bias varies linearly with elevation. To obtain the best fitting results, continuous piecewise linear functions were used to model the code-phase divergences. Note that the term continuous is used in the sense that the adjacent segments of the function share the same end point. Given a set of MP series, ( E i , M P i ) , which can be divided into several consecutive segments, the following objective functions should be satisfied to find the best fits:
S = j = 1 m 1 i = 1 n j ( f j , i M P j , i ) 2 = min s .t . f j ( a j ) f j + 1 ( a j ) = 0
where
  • i = 1 , 2 , n : point index; n : number of points (samples);
  • j = 1 , 2 , m 1 : segment index; m : number of nodes;
  • a j : elevation angle on a specific node (i.e., the end point of j-th segment);
  • b j : MP value on a specific node;
  • f j , i = f j ( E j , i , b 1 , b 2 , b m ) : fitted function used in the j-th segment.
Solving the normal linear system of Equation (3) with least squares method produces the MP values on all nodes. Then, the precision of estimates, i.e., the root mean square (RMS) is derived as:
R M S j = i = 1 n j ( f j , i M P j , i ) 2 n j 1
To ensure enough samples are used and to get reliable fitting functions, we selected the nodes to be separated by 10° of elevation, and the node values for the elevation angles range from 5° to 85°. Representative results of such continuous piecewise linear models are presented in Figure 4 and Figure 5. Figure 4 shows the results for JFNG, a MGEX station located in Wuhan, China, which is equipped with a Trimble NetR9 receiver. All available BeiDou data during the test period (i.e., January 2014–October 2015) were used to model the code biases on B1, B2, and B3 frequencies. Scatterplots on each node represent the MP values of different days, which clearly reflects the variation of code bias with time. The solid lines connect the averaged MP values on consecutive segments. Likewise, the code biases estimated from all the selected stations on the same day (DoY 001/2015) were shown in Figure 5. Differently, scatter plots on the same node represent MP values of different stations, which gives us an intuitive impression of their consistency.
Obviously, such code biases are orbit type-, elevation-, and frequency-dependent, which is consistent with Wanninger and Beer [22]. The code measurements of MEO satellites tracked on B1 frequency at high elevations are seriously affected by the code bias, which reaches over 0.5 m when the elevation rises up to 70°. The code biases on B2 and B3, however, are much smaller than those on B1. As to the IGSO group of satellites, similar effects can be found among triple frequencies. However, the values of code bias as well as their frequency-dependent differences are much smaller compared to those of MEOs. In general, the code biases show a good consistency among different stations (equipped with different receivers) and different days and, thus, make it possible to be correctly modeled with limited data. Nevertheless, it is worth noting that slight differences do exist in the models if different samples are used.
Therefore, in order to obtain the best fitting model of code biases, MP series derived from all the selected stations during a period of almost two years were used to produce the correction model for BeiDou satellite induced code bias. In addition to the correction values, the corresponding precisions are also essential for the users to have a better understanding of their corrections. Notice that there is very large RMS-values at low elevation for the reason that signals observed at low elevation angles are more likely to be affected by multipath error from stations, which decrease the precision of estimation. Figure 6 and Figure 7 show the resulting correction values and their RMS for elevations range from 5° to 85° with a node separation of 10°. The final model parameters are summarized in Table 2.
To make use of them, correction values for elevations between 5° and 85° can be obtained by linear interpolation from the nearest two nodes (Equation (5)), whereas the corrections on nodes of 5° and 85° are simply used for the elevation bins of (0,5) and (85,90). Assuming the MP correction values at node points are uncorrelated and neglecting the interpolation error, the precision of a given point can be determined by the law of propagation of variance in Equation (5):
M P ( e ) = M P 0 + ( M P 1 M P 0 ) e E 0 E 1 E 0 σ M P ( e ) 2 = ( E 1 e E 1 E 0 ) 2 σ 0 2 + ( e E 0 E 1 E 0 ) 2 σ 1 2
where E j , M P j and σ j 2 ( j = 1 , 2 ) are, respectively, the elevation angles, MP values and the corresponding RMS on the nearest two nodes; M P ( e ) and σ M P ( e ) are the interpolated MP correction and RMS at the given elevation angle e .
Figure 8 shows the correction differences with respect to Wanninger and Beer [22]. In general, the differences are mostly within 0.1 m, which implies that our model agrees well with the existing one. However, the values of the differences at low elevation bins are much bigger, which reach 0.2–0.4 m. This is caused by the significantly larger noise level at low elevation angles. The precision of MP corrections as listed in Table 2 confirms this.

5. Validation of the Models

To verify the proposed correction model, absolute positioning techniques, such as standard point positioning (SPP) and PPP, are good choices since they both are affected by the satellite-induced code biases described above. It should be noted that, however, such code biases may be obscured by the residual ionospheric delay and broadcast ephemeris (i.e., satellite orbits and clock corrections) errors, which are regarded as the dominating errors in SPP. Consequently, the code bias correction has marginal effect on the positioning results. To clearly demonstrate the effectiveness of the improved correction model, PPP tests were performed with three different processing schemes. In the first scheme, the satellite-induced code biases were kept without any treatment. In the second scheme, the traditional correction model provided in Wanninger and Beer [22] were used to remove the code biases. The improved correction model which contains correction values as well as precision indexes were applied for mitigating the same code biases in the last scheme. In the interest of brevity, “None”, “Traditional”, and “Improved” are used to denote the above-mentioned three different schemes throughout the rest of this article, if there is no additional explanation. The common processing strategies for PPP are summarized in Table 3. Precise BeiDou orbits and clock corrections provided by German Research Center for Geoscience (GFZ, Potsdam, Germany) were used. The positioning results were compared with the IGS weekly solutions [30]. In general, the reference coordinates have an accuracy of few millimeters.

5.1. Positioning Error

The validation of the correction is verified by PPP solution using five stations: CUT0, JFNG, XMIS (Christmas Island, Australia), GMSD (Nakatane town, Japan), NNOR (New Norcia, Australia), during days 45–50 in 2014 and days 130–135 in 2015 while the XMIS and SIN1 (Singapore) stations were excluded from the code bias modeling. Two representative PPP solutions on different stations and different days are shown in Figure 9 and Figure 10. Figure 9 shows the epoch-wise positioning error of CUT0 and JFNG recorded on DoY 047/2014, and Figure 10 shows the positioning error of CUT0 and XMIS recorded on DoY 130/2015. Owing to the highly-weighted carrier phase, the positioning results are not seriously affected by the un-modeled code biases. In general, an accuracy of 0.1 in horizontal and 0.2 m vertically can be achieved after convergence with the current BeiDou constellation even if we leave the code biases uncorrected. By comparing the first two schemes, the code bias corrections with the traditional model does not really yield noticeable improvements in terms of positioning. Sometimes inaccurate correction values may even degrade the positioning accuracy. However, the third scheme outperforms the former two in terms of positioning accuracy as well as convergence speed. This implies that the improved model, which provides not only correction values but also stochastic information, greatly enhances PPP performance, particularly at the initial stage.
To gain more insight into the PPP convergence, the three dimensional (3D) positioning errors of the first few hours on CUT0 and JFNG are shown in Figure 11. At the initial stage, large differences can be observed for both CUT0 and JFNG, indicating that the positioning results are quite sensitive to the different processing schemes. After a short time convergence (15–30 min), the positioning accuracy of CUT0 is noticeably improved by the code bias corrections with the traditional model, whereas the situation is totally different for JFNG. Note that, the convergence time means obtaining a 3D positioning error less than one decimeter, which has been adopted by Li and Zhang [17]. The convergence time for five stations under three schemes and the average convergence time are shown in Table 4.
In fact, the correction values may be not accurate enough. The inaccurate corrections may even degrade the precision of code observations. However, as already mentioned, the stochastic information of the code bias corrections are completely missing in the traditional model. In other words, the precision of the code observables has been overestimated in this case. This will, unfortunately, lead to an uncertainty of their corrections, and eventually the positioning accuracy cannot be always enhanced. Compared against the former two schemes, solutions with the improved model show the best positioning performance with an accuracy of ~0.1 m, suggesting that a proper stochastic model for the corrections is critical. Reviewing the RMS of the code bias corrections in Table 2, we can find that the precision of the provided corrections range 0.2–0.7 m. Such an uncertainty is comparable to the assumed noise level of code measurements. Therefore, the actual precision of the corrected code observations can be reflected in an objective manner once the stochastic model of the corrections is taken into account.

5.2. Wide-Lane Ambiguity

In addition to the positioning error, the wide-lane ambiguities, i.e., the so-called Melbourne-Wübbena (MW) linear combination [36,37], were derived from the B1/B2 dual-frequency code and carrier phase measurements. Figure 12 and Figure 13 shows the WL ambiguity and EWL ambiguity of representative IGSO and MEO satellites with different code bias correction schemes. As shown in the upper two plots of Figure 12 and Figure 13, the WL and EWL values derived from original code observables reveal strong elevation-dependent systematic variations with amplitudes of half a wide-lane cycle, which prevents precise ambiguity determination and successful ambiguity resolution. After application of the code bias correction models, the systematic variations have been greatly removed, and the resulting MW series run much more stable and closer to the nearest integers. The fixing rate of MW ambiguities increases from 80.4% to 91.8%. Additionally, the MW values of all involved IGSO (C6, C7, C8, C9, and C10) and MEO (C11, C12, and C14) satellites were calculated and the statistics were given in Figure 14. Likewise, the MW derivations from the latter two schemes agree well with each other. This is reasonable since they have consistent code bias correction values as shown in Figure 8. Moreover, the distributions of the corrected MW values are much more concentrated compared to the one calculated from uncorrected code measurements. Obviously, the corrected MW values are more likely to be fixed, which definitely contribute to the subsequent undifferenced ambiguity resolution.

5.3. Code Residual

Since the un-modeled code biases will be partly reflected in the residuals, the code residuals (i.e., residuals of the B1/B2 ionosphere-free code combination) were analyzed as well. Taking the same station, for example, Figure 15 shows the code residuals (Res) of representative IGSO and MEO satellites on XMIS. As expected, the uncorrected elevation-dependent code bias variations are visible in the code residuals for both IGSO and MEO satellites. Such variations are stronger for MEOs since they are more seriously affected by the code-phase divergence. After application of the code bias correction models, the residuals run much more like zero mean random noises, which implies that most of the code biases have been well settled. Moreover, residuals of all the involved IGSO and MEO satellites were grouped together and the statistics were shown in Figure 16 and Table 5. The code residuals show a mean of 9 mm for the IGSO group satellites, and a zero mean for the MEO group satellites for all the three schemes. For IGSOs, residual distributions of the latter two schemes (i.e., code biases corrected by the traditional and improved correction models) are slightly better than that of the first scheme (i.e., none code bias corrections applied). As for the MEOs, code residuals show much better distributions in terms of standard deviation after mitigating the code biases with either the traditional or the improved correction models.

6. Conclusions

Different from the other GNSS systems, BeiDou (IGSO and MEO) code measurements are polluted by the so-called satellite-induced code biases, or code-phase divergences, or code bias variations. The multipath combination (MP), which is a geometry-free and ionosphere-free combination, was used to investigate the code bias variations. It was confirmed that such code bias variations are orbit type-, frequency-, and elevation-dependent. In general, the code measurements of MEO satellites tracked on the B1 frequency at high elevations are more likely to be affected by the code bias variations.
To mitigate these unexpected variations, an improved correction model was developed in this paper. Approximately two years data recorded on a globally distributed MGEX stations equipped with different receiver types were used. The continuous piecewise linear functions were employed to produce the correction values as well as their precision indexes. To obtain the best-fitting results, the elevation angles were limited to 5°–85° with a node separation of 10°. To make full use of the improved model, the stochastic information of the corrections were encouraged to take into account.
PPP tests were conducted with different code bias variation treatments to verify the effectiveness of the proposed model. Positioning results show that the traditional model does not really yield noticeable improvements. Sometimes it may even degrade the positioning accuracy due to its inaccurate stochastic model. However, the positioning accuracy, as well as convergence, were obviously enhanced by the improved model. This implies that the precision of the uncorrected code measurements and even the corrected code measurements with traditional correction model may be too optimistic. However, a more real and reliable stochastic model can be set up if the precision of corrections is taken into account. Moreover, the wide-lane ambiguities and code residuals were derived for comparison, which further confirmed the effectiveness of the proposed model.

Acknowledgments

The authors gratefully acknowledge IGS Multi-GNSS Experiment (MGEX) for providing GNSS data and products. This study would not have been possible without such data. Thanks also go to the National Natural Science Foundation of China (No: 41404006, No: 41474025) and the Open Research Fund of State Key Laboratory of Information Engineering in Survey, Mapping and Remote Sensing, China (No. 15P02).

Author Contributions

Fei Guo and Wanke Liu conceived the idea and designed the experiments. Fei Guo and Xin Li performed the experiments and analyzed the data. All authors read, revised and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Montenbruck, O.; Hauschild, A.; Steigenberger, P.; Hugentobler, U.; Teunissen, P.; Nakamura, S. Initial Assessment of the COMPASS/BeiDou-2 Regional Navigation Satellite System. GPS Solut. 2013, 17, 211–222. [Google Scholar] [CrossRef]
  2. Zhao, Q.; Guo, J.; Li, M.; Qu, L.; Hu, Z.; Shi, C.; Liu, J. Initial results of precise orbit and clock determination for COMPASS navigation satellite system. J. Geod. 2013, 87, 475–486. [Google Scholar] [CrossRef]
  3. Deng, Z.; Zhao, Q.; Springer, T.; Prange, L.; Uhlemann, M. Orbit and Clock Determination-BeiDou. In Proceedings of the International GNSS Service (IGS) Workshop, Pasadena, CA, USA, 23–27 June 2014.
  4. Wu, X.; Hu, X.; Wang, G.; Zhong, H.; Tang, C. Evaluation of COMPASS ionospheric model in GNSS positioning. Adv. Space Res. 2013, 51, 959–968. [Google Scholar] [CrossRef]
  5. Zhang, R.; Song, W.; Yao, Y.; Shi, C.; Lou, Y.; Yi, W. Modeling regional ionospheric delay with ground-based BeiDou and GPS observations in China. GPS Solut. 2015, 19, 1–10. [Google Scholar] [CrossRef]
  6. Xu, A.; Xu, Z.; Ge, M.; Xu, X.; Zhu, H.; Sui, X. Estimating zenith tropospheric delays from BeiDou navigation satellite system observations. Sensors 2013, 13, 4514–4526. [Google Scholar] [CrossRef] [PubMed]
  7. Montenbruck, O.; Steigenberger, P.; Hauschild, A. Differential Code Bias Estimation using Multi-GNSS Observations and Global Ionosphere Maps. In Proceedings of the Institute of Navigation (ION) International Technical Meeting, San Diego, CA, USA, 26–28 January 2014.
  8. Guo, F.; Zhang, X.; Wang, J. Timing group delay and differential code bias corrections for BeiDou positioning. J. Geod. 2015, 89, 427–445. [Google Scholar] [CrossRef]
  9. Shi, C.; Fan, L.; Li, M.; Liu, Z.; Gu, S.; Zhong, S.; Song, W. An enhanced algorithm to estimate BDS satellite’s differential code biases. J. Geod. 2015, 90, 161–177. [Google Scholar] [CrossRef]
  10. Zhang, B.; Teunissen, P.J.G. Characterization of multi-GNSS between-receiver differential code biases using zero and short baselines. Chin. Sci. Bull. 2015, 21, 1840–1849. [Google Scholar] [CrossRef]
  11. Hauschild, A.; Montenbruck, O.; Sleewaegen, J.M.; Huisman, L.; Teunissen, P. Characterization of compass M-1 signals. GPS Solut. 2012, 16, 117–126. [Google Scholar] [CrossRef]
  12. Wang, G.; de Jong, K.; Zhao, Q.; Hu, Z.; Guo, J. Multipath analysis of code measurements for BeiDou geostationary satellites. GPS Solut. 2015, 19, 129–139. [Google Scholar] [CrossRef]
  13. Tang, W.; Deng, C.; Shi, C.; Liu, J. Triple-frequency carrier ambiguity resolution for Beidou navigation satellite system. GPS Solut. 2014, 18, 335–344. [Google Scholar] [CrossRef]
  14. Zhang, X.; He, X. Performance analysis of triple-frequency ambiguity resolution with BeiDou observations. GPS Solut. 2015, 20, 269–281. [Google Scholar] [CrossRef]
  15. Yang, Y.; Li, J.; Wang, A.; Xu, J.; He, H.; Guo, H.; Dai, X. Preliminary assessment of the navigation and positioning performance of BeiDou regional navigation satellite system. Sci. China Earth Sci. 2014, 57, 144–152. [Google Scholar] [CrossRef]
  16. Li, M.; Qu, L.; Zhao, Q.; Guo, J.; Su, X.; Li, X. Precise Point Positioning with the BeiDou Navigation Satellite System. Sensors 2014, 14, 927–943. [Google Scholar] [CrossRef] [PubMed]
  17. Li, X.; Zhang, X.; Ren, X.; Fritsche, M.; Wickert, J.; Schuh, H. Precise positioning with current multi-constellation Global Navigation Satellite Systems: GPS, GLONASS, Galileo and BeiDou. Sci. Rep. 2015, 5. [Google Scholar] [CrossRef] [PubMed]
  18. Li, X.; Ge, M.; Dai, X.; Ren, X.; Fritsche, M.; Wickert, J.; Schuh, H. Accuracy and reliability of multi-GNSS real-time precise positioning: GPS, GLONASS, BeiDou, and Galileo. J. Geod. 2015, 89, 607–635. [Google Scholar] [CrossRef]
  19. Guo, F.; Zhang, X.; Wang, J.; Ren, X. Modeling and assessment of triple-frequency BDS precise point positioning. J. Geod. 2016. [Google Scholar] [CrossRef]
  20. Nadarajah, N.; Teunissen, P.J.G.; Raziq, N. BeiDou inter-satellite-type bias evaluation and calibration for mixed receiver attitude determination. Sensors 2013, 7, 9435–9463. [Google Scholar] [CrossRef] [PubMed]
  21. Nadarajah, N.; Teunissen, P.J.G.; Sleewaegen, J.M.; Montenbruck, O. The mixed-receiver BeiDou inter-satellite-type bias and its impact on RTK positioning. GPS Solut. 2014, 19, 357–368. [Google Scholar] [CrossRef]
  22. Wanninger, L.; Beer, S. BeiDou satellite-induced code pseudorange variations: Diagnosis and therapy. GPS Solut. 2015, 19, 639–648. [Google Scholar] [CrossRef]
  23. Montenbruck, O.; Hauschild, A.; Steigenberger, P.; Hugentobler, U.; Riley, S. A COMPASS for Asia: First Experience with the BeiDou-2 Regional Navigation System. Poster at IGS Workshop in Olsztyn, Poland 2012. Available online: http://www.igs.org/presents/poland2012/posters (accessed on 10 May 2015).
  24. Montenbruck, O.; Rizos, C.; Weber, R.; Weber, G.; Neilan, R.; Hugentobler, U. Getting a grip on multi-GNSS-the international GNSS service MGEX campaign. GPS World 2013, 24, 44–49. [Google Scholar]
  25. De Bakker, P.F.; Tiberius, C.C.J.M.; van der Marel, H.; van Bree, R.J.P. Short and zero baseline analysis of GPS L1 C/A, L5Q, GIOVE E1B, and E5aQ signals. GPS Solut. 2012, 16, 53–64. [Google Scholar] [CrossRef]
  26. Springer, T.; Dilssner, F. SVN49 and Other GPS Anomalies; Inside GNSS: Eugene, OR, USA, 2009; pp. 32–36. [Google Scholar]
  27. Montenbruck, O.; Hauschild, A.; Steigenberger, P.; Langley, R.B. Three’s the challenge. GPS World 2010, 21, 8–19. [Google Scholar]
  28. Pearson, K. Note on regression and inheritance in the case of two parents. Proc. R. Soc. Lond. 1895, 58, 240–242. [Google Scholar] [CrossRef]
  29. Montenbruck, O.; Steigenberger, P.; Khachikyan, R.; Weber, G.; Langley, R.B.; Mervart, L.; Hugentobler, U. IGS-MGEX: Preparing the Ground for Multi-Constellation GNSS Science; Inside GNSS: Eugene, OR, USA, 2014; Volume 9, pp. 42–49. [Google Scholar]
  30. Dow, J.M.; Neilan, R.E.; Rizos, C. The International GNSS Service in a changing landscape of Global Navigation Satellite Systems. J. Geod. 2009, 83, 191–198. [Google Scholar] [CrossRef]
  31. Zhang, X.; Andersen, O.B. Surface ice flow velocity and tide retrieval of the Amery ice shelf using precise point positioning. J. Geod. 2006, 80, 171–176. [Google Scholar] [CrossRef]
  32. Boehm, J.; Heinkelmann, R.; Schuh, H. A global model of pressure and temperature for geodetic applications. J. Geod. 2007, 81, 679–683. [Google Scholar] [CrossRef]
  33. Petit, G.; Luzum, B. IERS Conventions 2010 (IERS Technical Note No. 36); International Earth Rotation and Reference Systems Service: Frankfurt, Germany, 2010. [Google Scholar]
  34. Rizos, C.; Montenbruck, O.; Weber, R.; Neilan, R.; Hugentobler, U. The IGS MGEX Experiment as a Milestone for a Comprehensive Multi-GNSS Service. In Proceedings of Institute of Navigation Pacific Positioning, Navigation and Timing (ION Pacific PNT), Honolulu, HI, USA, 23–25 April 2013; pp. 289–295.
  35. Wu, J.; Wu, S.; Hajj, G.; Bertiger, W.; Lichten, S. Effects of antenna orientation on GPS carrier phase. Manuscr. Geod. 1993, 18, 91–98. [Google Scholar]
  36. Melbourne, W.G. The case for ranging in GPS-based geodetic systems. In Proceedings of the first International Symposium on Precise Positioning with the Global Positioning System, Rockville, MD, USA, 15–19 April 1985; pp. 373–386.
  37. Wübbena, G. Software developments for geodetic positioning with GPS using TI 4100 code and carrier measurements. In Proceedings of the 1st International Symposium on Precise Positioning with the Global Positioning Systems, Rockville, MD, USA, 15–19 April 1985; pp. 403–412.
Figure 1. BeiDou and GPS code MP series against elevation angles observed from JFNG (DoY 081/2015).
Figure 1. BeiDou and GPS code MP series against elevation angles observed from JFNG (DoY 081/2015).
Sensors 16 00909 g001
Figure 2. Distribution of the selected MGEX stations used in this study.
Figure 2. Distribution of the selected MGEX stations used in this study.
Sensors 16 00909 g002
Figure 3. Number of MP samples against elevation bins with an interval of 10° for BeiDou MEO (left) and IGSO satellites (right).
Figure 3. Number of MP samples against elevation bins with an interval of 10° for BeiDou MEO (left) and IGSO satellites (right).
Sensors 16 00909 g003
Figure 4. Piecewise linear MP models for BeiDou MEO (left) and IGSO satellites (right) on JFNG station using data from different days (during January 2014–October 2015).
Figure 4. Piecewise linear MP models for BeiDou MEO (left) and IGSO satellites (right) on JFNG station using data from different days (during January 2014–October 2015).
Sensors 16 00909 g004
Figure 5. Piecewise linear MP models for BeiDou MEO (left) and IGSO satellites (right) using data from different stations (equipped with six Trimble receivers and four Septentrio receivers).
Figure 5. Piecewise linear MP models for BeiDou MEO (left) and IGSO satellites (right) using data from different stations (equipped with six Trimble receivers and four Septentrio receivers).
Sensors 16 00909 g005
Figure 6. Suggested correction values of code bias (orbit type-, frequency-, and elevation-dependent).
Figure 6. Suggested correction values of code bias (orbit type-, frequency-, and elevation-dependent).
Sensors 16 00909 g006
Figure 7. Suggested RMS of correction values (orbit type-, frequency-, and elevation-dependent).
Figure 7. Suggested RMS of correction values (orbit type-, frequency-, and elevation-dependent).
Sensors 16 00909 g007
Figure 8. Differences of the code bias correction values compared with those of Wanninger and Beer.
Figure 8. Differences of the code bias correction values compared with those of Wanninger and Beer.
Sensors 16 00909 g008
Figure 9. Epoch-wise PPP solutions of CUT0 and JFNG (DoY 047/2014).
Figure 9. Epoch-wise PPP solutions of CUT0 and JFNG (DoY 047/2014).
Sensors 16 00909 g009
Figure 10. Epoch-wise PPP solutions of CUT0 and XMIS (DoY 130/2015).
Figure 10. Epoch-wise PPP solutions of CUT0 and XMIS (DoY 130/2015).
Sensors 16 00909 g010
Figure 11. 3D positioning error of the first two hours on CUT0 and JFNG (DoY 047/2014).
Figure 11. 3D positioning error of the first two hours on CUT0 and JFNG (DoY 047/2014).
Sensors 16 00909 g011
Figure 12. Time series of the derived MW value on XMIS (DoY 130/2015).
Figure 12. Time series of the derived MW value on XMIS (DoY 130/2015).
Sensors 16 00909 g012
Figure 13. Time series of the derived EWL value on CUT0 (DoY 122/2015).
Figure 13. Time series of the derived EWL value on CUT0 (DoY 122/2015).
Sensors 16 00909 g013
Figure 14. Histograms of all the derived MW values (left: IGSOs; right: MEOs) on XMIS (DoY 130/2015).
Figure 14. Histograms of all the derived MW values (left: IGSOs; right: MEOs) on XMIS (DoY 130/2015).
Sensors 16 00909 g014
Figure 15. Time series of representative code residuals on XMIS (DoY 130/2015).
Figure 15. Time series of representative code residuals on XMIS (DoY 130/2015).
Sensors 16 00909 g015
Figure 16. Histograms of all the code residuals (left: IGSOs; right: MEOs) on XMIS (DoY 130/2015).
Figure 16. Histograms of all the code residuals (left: IGSOs; right: MEOs) on XMIS (DoY 130/2015).
Sensors 16 00909 g016
Table 1. Correlations between MP and elevation series.
Table 1. Correlations between MP and elevation series.
GPSBeiDou
GEOIGSOMEO
MP1−0.04−0.01−0.49−0.61
MP20.080.00−0.48−0.60
MP30.12−0.02v0.39−0.36
Table 2. Correction values and their precisions for BeiDou MEO/IGSO code measurements.
Table 2. Correction values and their precisions for BeiDou MEO/IGSO code measurements.
Elevation Nodes (°)Correction Values (m)RMS of Corrections (m)
MEOsIGSOsMEOsIGSOs
B1B2B3B1B2B3B1B2B3B1B2B3
5−0.109−0.140−0.060−0.101−0.148−0.0650.7210.5880.5800.7090.5640.576
15−0.169−0.148−0.087−0.203−0.250−0.1620.6050.4800.4990.6510.5320.582
25−0.150−0.121−0.070−0.222−0.224−0.1680.4760.3730.4010.5000.3710.409
35−0.105−0.062−0.053−0.123−0.110−0.0780.3880.2910.2900.4030.2970.303
450.0040.0470.022−0.066−0.043−0.0490.3330.2540.2580.3890.2780.244
550.1810.1850.0960.0360.0440.0210.2930.2200.2410.3080.2300.223
650.4110.3260.1800.1070.1060.0680.2750.1940.2110.2620.2100.208
750.6740.4770.2800.1630.1780.1300.2610.1880.2060.2510.2130.212
850.8530.6000.3730.2450.2600.2080.2330.1730.1980.2170.1950.190
Table 3. PPP processing strategy.
Table 3. PPP processing strategy.
ItemsModels
Estimator (engine)Kalman filter, TriP software [31]
ObservationsB1/B2 code and carrier phase measurements
Sampling rate30 s
Elevation cutoff10°
Weighting schemeElevation dependent; 3 mm and 0.3 m for raw phase and code, respectively
Ionospheric delayEliminated by Ionosphere-free combination (s)
Tropospheric delayDry component: corrected with GPT model [32] wet component: estimated as random-walk process, GMF mapping function
Relativistic EffectApplied
Station displacementCorrected by IERS Convention 2010, including Solid Earth tide and ocean tide loading [33]
Satellite antenna phase center offsetCorrected with conventional PCO values from MGEX [34]
Receiver antenna phase center offsetCorrected by the same PCO values as GPS
Phase-windup effectCorrected [35]
Receiver clockEstimated, epoch-wise solution
Station coordinateEstimated, epoch-wise solution
Phase ambiguitiesEstimated, constant for each arc; float value
Table 4. Convergence time under three schemes for five stations.
Table 4. Convergence time under three schemes for five stations.
Convergence Time (min)Schemes
UncorrectedTraditionalImproved
GMSD125108.529.5
CUT013013436
JFNG7613632.5
XIMS15716255
NNOR13012733.5
Average123.6133.537.3
Table 5. Mean and standard deviation (STD) values of code residuals with different processing schemes (unit: m).
Table 5. Mean and standard deviation (STD) values of code residuals with different processing schemes (unit: m).
IGSOMEO
NoneTraditionalImprovedNoneTraditionalImproved
Mean0.0090.0090.0090.0000.0000.000
STD1.1661.1551.1521.8021.7301.740

Share and Cite

MDPI and ACS Style

Guo, F.; Li, X.; Liu, W. Mitigating BeiDou Satellite-Induced Code Bias: Taking into Account the Stochastic Model of Corrections. Sensors 2016, 16, 909. https://doi.org/10.3390/s16060909

AMA Style

Guo F, Li X, Liu W. Mitigating BeiDou Satellite-Induced Code Bias: Taking into Account the Stochastic Model of Corrections. Sensors. 2016; 16(6):909. https://doi.org/10.3390/s16060909

Chicago/Turabian Style

Guo, Fei, Xin Li, and Wanke Liu. 2016. "Mitigating BeiDou Satellite-Induced Code Bias: Taking into Account the Stochastic Model of Corrections" Sensors 16, no. 6: 909. https://doi.org/10.3390/s16060909

APA Style

Guo, F., Li, X., & Liu, W. (2016). Mitigating BeiDou Satellite-Induced Code Bias: Taking into Account the Stochastic Model of Corrections. Sensors, 16(6), 909. https://doi.org/10.3390/s16060909

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop