Next Article in Journal
Special Issue “Advanced Signal Processing in Intelligent Systems for Health Monitoring”
Previous Article in Journal
A Simple and Practical Single-Camera Stereo-Digital Image Correlation Using a Color Camera and X-Cube Prism
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of BDS Fractional Cycle Biases and PPP Ambiguity Resolution

1
School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China
2
GNSS Research Center, Wuhan University, Wuhan 430079, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(21), 4725; https://doi.org/10.3390/s19214725
Submission received: 5 October 2019 / Revised: 26 October 2019 / Accepted: 29 October 2019 / Published: 31 October 2019
(This article belongs to the Section Remote Sensors)

Abstract

:
It is difficult to enable traditional precise point positioning (PPP) with ambiguity resolution (AR) due to fractional cycle biases (FCBs), which limit the accuracy and reliability of positioning results. The BeiDou Navigation Satellite System (BDS) has been providing continuous positioning, navigation, and timing (PNT) services in the global region since the end of 2018. The BDS constellation includes geostationary earth orbit (GEO), inclined geostationary orbit (IGSO), and medium earth orbit (MEO) satellites. However, its hybrid constellation structure and the satellite-side multipath effect have hindered the BDS PPP AR for two main reasons: (1) some receivers have half-cycle biases between GEO and non-GEO satellites, which result in the inconsistency of hardware delays for each satellite type; (2) the correction model for elevation-dependent satellite-side multipath effect is only applicable to IGSO and MEO, while in the case of GEO the effect cannot be effectively weakened or eliminated. To rectify these problems, a suitable strategy is proposed for estimating BDS FCBs, whereby the GEO FCBs and non-GEO FCBs are estimated independently. Results show that the FCBs estimated by the new strategy for GEO and non-GEO are more stable compared to the traditional strategy. The GEO wide-lane (WL) FCBs fluctuate less than 0.3 cycle in one month, except for C05, while the variation of non-GEO WL FCBs is about 0.1 cycle. In addition, compared to the traditional strategy, the fractions of GEO WL ambiguities after the removal of FCBs estimated by the new strategy can be improved noticeably from 53.5% to 78.5%, and from 71.8% to 92.3% for <0.15 cycle and <0.25 cycle respectively, which could be comparable with non-GEO. Simultaneously, the improvement of GEO narrow-lane (NL) ambiguities is from 28.9% to 40.2%, and from 40.4% to 53.3% for <0.10 cycle and <0.15 cycle respectively, are less noticeable. This is mainly due to the low precision IGS products for GEO. After PPP AR, the mean convergence time is shorted from 56.0 min to 43.6 min, and from 71.6 min to 62.7 min for static PPP and kinematic PPP, respectively.

1. Introduction

Precise point positioning (PPP) can obtain high precision and absolute positioning based on precise satellite orbit and clock offset products [1]. The BeiDou Navigation Satellite System (BDS) has been providing global services since 27 December 2018, and the BDS satellites for precise positioning include geostationary earth orbit (GEO), inclined geostationary orbit (IGSO), and medium earth orbit (MEO). Due to the unique constellation design of BDS, the number of available satellites in the Asia-Pacific region is higher than in other areas. As we know, ambiguity resolution (AR) technology can shorten the convergence time of PPP and improve positioning accuracy and stability. However, the BDS hybrid constellation and the multipath effect in the satellite side significantly hinder ambiguity resolution, especially for PPP, which limits its application scope substantially. Many studies have found that hardware delays exist at the satellite and receiver terminals, which are highly coupled with the ambiguity parameters. It is very difficult to separate them from each other [2,3,4], leading to the loss of the integer property of ambiguity parameters. In relative positioning, the hardware delays can be eliminated by double-difference (DD), so the DD ambiguity has the integer characteristics. For PPP, it is necessary to restore integer characteristics for ambiguities by using satellite fractional cycle biases (FCBs) products [3] or precise clock offset products with FCBs corrections [4]. In addition, the ambiguity of PPP can also be fixed by using the decoupled clock model [5] or the zero-differenced hardware delays model [6]. In theory, these above PPP AR methods are also applicable to BDS.
However, some studies have found that there are some unique errors in the BDS satellite terminals, such as the severe multipath effect [7,8], which makes it difficult to fix the ambiguities for BDS PPP directly by using the above methods. Although the multipath of BDS MEO and IGSO satellites can be better weakened by the piecewise linear correction model based on the elevation angle [9], this model is not applicable to GEO. Because the elevation angle of a station-GEO pair is nearly unchanged, the deviation is an approximate constant which is difficult to be modeled. It is worth noting that the influence of the deviation on the observations is not the same for different station-GEO pairs. Therefore, the multipath effect at the GEO satellite end remains difficult to be effectively weakened or eliminated. Furthermore, the BDS GEO satellite signal is also affected by the multipath at the station, which shows a strong correlation between different days [10]. Moreover, for some commonly used receivers, such as the Trimble NETR9 and JAVAD TRE G3T DELTA, there is Inter-satellite-type bias (ISTB) between the BDS GEO satellites and the non-GEO satellites, which may cause half-cycle jumps for the observations of these receiver types [11,12]. For these reasons, we investigate the FCB estimation strategy of the BDS in depth, analyze the BDS ISTB and the satellite-side multipath effect in detail, and propose a suitable strategy for BDS FCBs estimation to achieve PPP fix solution.
After Section 1, the article is presented as follows: The BDS PPP model and FCBs estimation method are described in Section 2. Section 3 discusses the ISTB, BDS satellite multipath and BDS FCBs estimation strategies, and the characteristics of BDS FCBs in detail. Section 4 uses experiments to analyze the positioning accuracy and stability of the BDS static and kinematic PPP. Finally, conclusions are provided in Section 5.

2. Methods

The normal PPP model is described in the following equations [13,14]:
P r , i s = ρ r s + c ( d t r d t s ) + f 1 2 f i 2 I r , 1 s + T r s + d r , i d i s + ε P r , i s
L r , i s = ρ r s + c ( d t r d t s ) f 1 2 f i 2 I r , 1 s + T r s + λ i ( b r , i b i s ) + λ i N r , i s + ε L r , i s
where r , s , i represent the receiver, satellite and frequency number, respectively; f i , λ i donate the frequency i and wavelength;   d t r , d t s are the clock offsets of receiver and satellite in seconds; I r , 1 s is the slant ionospheric delay at 1-frequency in meters; T r s is the slant troposphere delay in meters; d r , i ,   d i s are the code hardware delays of receiver and satellite at the i-frequency in meters; b r , i , b i s are the phase delays of receiver and satellite at the i-frequency in cycles, respectively; N r , i s is the integer ambiguity in cycles; ε P r , i s , ε L r , i s are the measurement noises of pseudorange and carrier phase in meters, respectively. For simplicity, we assume that phase windup, phase center variation, relativistic effect, and tidal load have been corrected in Equations (1) and (2).
The ionospheric-free (IF) linear combination can be written as follows:
P r , i f s = ρ r s + c ( d t r d t s ) + T r s + d r , i f d i f s + ε P r , i f s
L r , i f s = ρ r s + c ( d t r d t s ) + T r s + λ i f ( b r , i f b i f s ) + λ i f N r , i f s + ε L r , i f s
In PPP, the IGS precise clock products are generally used and the above equations could be re-written as follows:
P r , i f s = ρ r s + c ( d t r d t s ) + T r s + ε P r , i f s
L r , i f s = ρ r s + c ( d t r d t s ) + T r s + λ i f B r , i f s + ε L r , i f s
d t r = d t r + d r , i f / c ,     d t s = d t s + d i f s / c
B r , i f s = f 1 f 1 + f 2 B r , 1 s + f 1 f 2 f 1 2 f 2 2 N r , w l s
B r , 1 s = N r , 1 s + b r , 1 s
b r , 1 s = f 1 f 1 f 2 ( b r , 1 b 1 s ) f 2 f 1 f 2 ( b r , 2 b 2 s ) f 1 + f 2 f 1 · d r , i f d i f s λ i f
In Equations (7)–(10), the actual estimated IF ambiguity will absorb the pseudorange and phase hardware delay parameters, which is the reason for traditional PPP to only obtain float solutions.
In order to obtain the PPP fixed solution, it is necessary to separate the corresponding hardware delay parameters from the ambiguity parameters. The un-differenced ambiguity can be expressed as the following equation:
B r , i s = N r , i s + b r , i s
B r , i s and N r , i s represent the float and integer ambiguity in cycles, respectively;   b r , i s = b r , i b i s ; b r , i , b i s are the actual estimated FCBs of receiver and satellite in cycles, respectively.
Generally, wide-lane (WL) ambiguity is fixed by the Melbourne-Wübbena (M-W) combination [15,16] using the following equations:
B r , w l s = N r , w l s + b r , w l s
b r , w l s = b r , w l + b w l s
b r , w l = ( b r , 1 b r , 2 ) f 1 f 2 f 1 + f 2 ( d r , 1 λ 1 + d r , 2 λ 2 )
b w l s = ( b 1 s b 2 s ) f 1 f 2 f 1 + f 2 ( d 1 s λ 1 + d 2 s λ 2 )
It should be noted that the integer part of b r , w l cannot be separated from the ambiguity parameter, which is contained by the N r , w l s . But the effect of it will be balanced by narrow-lane (NL) FCB [3].
Assuming that n stations observe m visible satellites, for station k with m satellites, Equation (11) can be expressed as matrix in the following equation:
[ B k , i 1 B k , i 2 B k , i m ] = [ I m × m   R m × 1   S m × m ] [ N m × 1 b k , i b i   m × 1 s ] ,   Q k , i ( m × m )
where I m × m donates the unit matrix; R m × 1 donates the m dimension column vector and all elements are 1; S m × m is the m dimension diagonal matrix, all elements are −1; N m × 1 is the m dimension integer ambiguity column vector; b k , i is the FCB of receiver k; b i   m × 1 s is the m dimension vector of satellites FCBs; Q k , i ( m × m ) is the variance-covariance matrix.
The mathematical model of the FCB estimation has been discussed in detail in Section 2. However, many studies have found that there are some unique errors in BDS, such as the severe multipath effect in the satellite side and Inter-satellite-type bias, which may affect the accuracy of the BDS FCB. In Section 3, the characteristic of these biases will be analyzed.

3. BDS FCB Estimation

3.1. Inter-Satellite-Type Bias

Researches present that there are half-cycle jumps between GEO and non-GEO satellites for some receiver types. Thus, in this case, the DD ambiguities of GEO and non-GEO satellites for BDS zero-baseline cannot be fixed [11,12]. Due to different decoding methods of GEO and non-GEO signals in different receivers, ISTB may be introduced into the phase observations for BDS, which is related to frequency, receiver type and firmware version at the receiver side. For two receivers of the same type, ISTB can be eliminated by double-difference. But for two different types of receivers, ISTB may cause a half-cycle jump in DD ambiguity, which seriously damages its integer characteristic. Considering the influence of ISTB, observation Equations (1) and (2) can be rewritten as:
P r , i s = ρ r s + c ( d t r d t s ) + f 1 2 f i 2 I r , 1 s + T r s + d r , i d i s + ε P r , i s
L r , i s = ρ r s + c ( d t r d t s ) f 1 2 f i 2 I r , 1 s + T r s + λ i ( b r , i β b i s ) + λ i N r , i s + ε L r , i s
where b r , i β = b r , i + I S T B r , i β ; I S T B r , i β is the ISTB at frequency i of receiver r in cycles; β is the number of satellite types for BDS. We assumed that the GEO satellites are type 1, and the non-GEO satellites are type 2. ISTB will be absorbed by b r , i β .
As we know, DD ambiguity can eliminate the influence of hardware delays perfectly. In order to show the influence of ISTB, a new analysis method is proposed for ISTB in this paper and DD-ambiguity is chosen as the research object. From the (8), (9), (10) and (12), the DD ambiguities are expressed in the following formulas:
λ i f N r 1 r 2 , i f s 1 s 2 = c f 1 f 1 2 f 2 2 ( N r 1 r 2 , 1 s 1 s 2 + I S T B r 1 r 2 , 1 β 1 β 2 ) c f 2 f 1 2 f 2 2 ( N r 1 r 2 , 2 s 1 s 2 + I S T B r 1 r 2 , 2 β 1 β 2 )
N r 1 r 2 , w l s 1 s 2 = N r 1 r 2 , w l s 1 s 2 + I S T B r 1 r 2 , w l β 1 β 2
N r 1 r 2 , n l s 1 s 2 = f 1 + f 2 f 1 N r 1 r 2 , i f s 1 s 2 f 2 f 1 f 2 N r 1 r 2 , w l s 1 s 2
N r 1 r 2 , i f s 1 s 2 = N r 1 r 2 , i f s 1 s 2 + I S T B r 1 r 2 , i f β 1 β 2
I S T B r 1 r 2 , i β 1 β 2 = ( I S T B r 1 , i β 1 I S T B r 2 , i β 1 ) ( I S T B r 1 , i β 2 I S T B r 2 , i β 2 )
Clearly, ISTB does not affect the DD ambiguity for two stations with the identical receiver type or the same observed satellite type. When it is an odd multiple of half-cycle, the fraction of   I S T B r 1 r 2 , i β 1 β 2 is not equal to zero and N r 1 r 2 , i s 1 s 2 equals N r 1 r 2 , i s 1 s 2 plus I S T B r 1 r 2 , i β 1 β 2 , which will seriously destroy the integer characteristics of DD ambiguity. From the (19)–(23), it can be deduced that the difference between the actual estimated DD ambiguity N r 1 r 2 , i s 1 s 2 and its true value N r 1 r 2 , i s 1 s 2 is as follows:
Δ I S T B , i f = N r 1 r 2 , i f s 1 s 2 N r 1 r 2 , i f s 1 s 2 = I S T B r 1 r 2 , i f β 1 β 2 = f 1 2 f 1 2 f 2 2 I S T B r 1 r 2 , 1 β 1 β 2 f 1 f 2 f 1 2 f 2 2 I S T B r 1 r 2 , 2 β 1 β 2
Δ I S T B , w l = N r 1 r 2 , w l s 1 s 2 N r 1 r 2 , w l s 1 s 2 = I S T B r 1 r 2 , w l β 1 β 2 = I S T B r 1 r 2 , 1 β 1 β 2 I S T B r 1 r 2 , 2 β 1 β 2  
Δ I S T B , n l = N r 1 r 2 , n l s 1 s 2 N r 1 r 2 , n l s 1 s 2 = f 1 + f 2 f 1 ( N r 1 r 2 , i f ' s 1 s 2 N r 1 r 2 , i f s 1 s 2 ) f 2 f 1 f 2 ( N r 1 r 2 , w l ' s 1 s 2 N r 1 r 2 , w l s 1 s 2 ) = f 1 + f 2 f 1 Δ I S T B , i f f 2 f 1 f 2 i n t ( Δ I S T B , w l )
Only fractions of ISTB could be acquired from ambiguities, so all the possible ISTB values for GEO and non-GEO combination are I S T B r 1 r 2 , 1 β 1 β 2 = ( 0.5 , 0 , 0.5 ) ,   I S T B r 1 r 2 , 2 β 1 β 2 = ( 0.5 , 0 , 0.5 ) .

Analyses

The numerical analyses in this contribution are based on BDS observations collected with two different receiver types, namely Trimble and Javad. Two zero baselines were selected from IGS ZIM2-ZIM3 on DOY 150, 2015 and CUCC-CUTC on DOY 151, 2018, in Curtin University (http://saegnss2.curtin.edu.au/ldc/). The receiver information is listed in Table 1 and Table 2.
According to Table 3, when I S T B r 1 r 2 , 1 β 1 β 2 equals ± 0.5 and I S T B r 1 r 2 , 2 β 1 β 2 equals 0, Δ I S T B , w l equals ± 0.5 . After Δ I S T B , w l is rounded, we may get two possible values because of numerical error ε , so int ( Δ I S T B , w l ) = i n t ( 0.5 ± ε ) = 0   o r   1 , Δ I S T B , n l = f 1 2 ( f 1 f 2 )   o r   f 1 2 f 2 2 ( f 1 f 2 ) . Considering f1 = 1561.098, f2 = 1207.140,   Δ I S T B , n l   = 2.205   o r 1.205 can be derived. In this case, the fractions of WL DD ambiguity between GEO and non-GEO satellites will be 0.5 cycle, while that of NL DD ambiguity will be −0.8, −0.2, 0.2 and 0.8 (Figure 1). Similarly, when I S T B r 1 r 2 , 1 β 1 β 2 = ± 0.5   and   I S T B r 1 r 2 , 2 β 1 β 2 = ± 0.5 , int ( Δ I S T B , w l ) = 1   o r   1   and   Δ I S T B , n l = ( 0.5 , 0.5 ) could be derived from Table 4. That is, when I S T B r 1 r 2 , 1 β 1 β 2 and I S T B r 1 r 2 , 2 β 1 β 2 equal ± 0.5 , the fractional parts of NL DD ambiguity between GEO and non-GEO will be 0.5 cycle, while that of WL DD ambiguity will be an integer (Figure 2).
From Equation (11), if I S T B r , i 1 is not equal to I S T B r , i 2 , GEO and non-GEO FCBs are different at receiver side, thus two FCB parameters should be considered. For Multi-GNSS Experiment (MGEX) networks, there are dozens of types and versions of receivers, among which some receivers may be I S T B r , i 1 I S T B r , i 2 . Obviously, it is essential to accurately identify ISTB inconsistent receiver types from MGEX networks. However, this is not only onerous but is also unrealistic.

3.2. BDS Satellite-Induced Code Pseudorange Variations

Hauschild et al. found that the BDS pseudorange observations are unstable and have systematic deviation, which is strongly related to the satellite elevation angle [7,8]. Wanninger and Beer found that these code variations are elevation-dependent and there are no dependencies on receiver type, observation time, or satellite azimuth [9]. An elevation-dependent model was proposed for IGSO and MEO, which can weaken the fluctuation of IGSO/MEO pseudorange observations well. However, an absolute determination of the correction model is not feasible due to the ambiguous carrier phase observations and unknown signal delays, which will affect the FCBs [9]. This model is not applicable to the GEO, because the elevation angle of a GEO satellite relative to a certain station is basically unchanged, which is difficult to be modeled. Due to the relationship between the satellite-induced code variations of BDS and elevation, formulas are shown as:
P r , i s = ρ r s + c ( d t r d t s ) + f 1 2 f i 2 I r , 1 s + T r s + d r , i d i s + c o d e _ b i a s P r , i s ( elev ) + ε P r , i s
L r , i s = ρ r s + c ( d t r d t s ) f 1 2 f i 2 I r , 1 s + T r s + λ i ( b r , i β b i s ) + λ i N r , i s + ε L r , i s
c o d e _ b i a s P r , i s ( elev ) is the satellite-induced code variation of receiver r with satellite s in meters.
B r , w l s = N r , w l s + b r , w l β b w l s c o d e _ b i a s P r , n l s ( elev )
b r , w l β = ( b r , 1 β b r , 2 β ) f 1 f 2 f 1 + f 2 { d r , 1 d 1 s λ 1 + d r , 2 d 2 s λ 2 }
c o d e b i a s P r , n l s ( elev ) = f 1 f 2 f 1 + f 2 { c o d e b i a s P r , n l s ( elev ) λ 1 + c o d e b i a s P r , n l s ( elev ) λ 2 }
Equation (29) shows that the M-W observations will be affected by the pseudorange deviation. The effect is verified by M-W observations of B1 and B2 on DOY 083 2018 of MGEX station XMIS for the BDS. Figure 3 shows the variation of the M-W observations of C01 and C02 with elevation angle, while Figure 4 and Figure 5 show the variation of the M-W observations of C06 and C11.
From Figure 3, although the elevation of the BDS GEO satellite is almost unchanged, its M-W values show periodic fluctuations.
Figure 4 demonstrates that the fluctuation of the M-W value for C06 likes the shape of M, and the change of its elevation angle is also similar to M. They are very similar and highly correlated. In Figure 5, the M-W value of C11 has a very similar shape with its elevation angle. It can be seen that the M-W values of the BDS IGSO and MEO satellites are strongly correlated with elevation angles, which is mainly caused by the satellite-induced code pseudorange variations. Satellite-induced code variation is a unique bias in BDS, so it must be considered in BDS PPP AR.
The piecewise linear correction model proposed by Wanninger and Beer [9] is used to correct the multipath effect for BDS IGSO and MEO. However, this model is not applicable to GEO satellites. Figure 6 and Figure 7 display the corrected M-W values of C06 and C11 vary with the satellite elevation angle.
In Figure 6 and Figure 7, the corrected M-W values tend to be stable within an arc. However, it can be seen that at a low elevation angle, there are still some fluctuations, indicating that some errors still exist after the correction, and the code multipath effect has not been completely eliminated. Due to the limitation of the piecewise correction model, it is necessary to explore a more accurate model. Since there is no better correction model to be used, the piecewise model is used in this study.

3.3. Strategy Comparison

From the above analysis in Section 3.1 and Section 3.2, we can see that there are ISTB and pseudorange deviation problems in BDS. The ISTB results in the inconsistency of GEO and non-GEO hardware delay of some stations, which need to be estimated separately. For the BDS satellite-induced code variations, it has to be emphasized that an absolute determination of the correction models is not feasible due to the ambiguous phase observations and unknown signal delays. Hence, the corrections will affect the FCB estimations, which are essential for fixing ambiguity in PPP [9]. Meanwhile, because of the special orbit type of the GEO satellite, its elevation angle is almost unchanged, so the variation should be an approximate constant. But for different stations and different GEO, the elevation angles are not the same, so the variations are not the same consistency for different station-GEO pairs. Therefore, it will not be completely absorbed by the receiver or satellite hardware delay. In addition, the precision of BDS GEO products released by IGS is lower than that of IGSO and MEO.
Based on the above analyses, the paper proposes a new strategy for the classified estimation of BDS FCBs: GEO FCBs, and non-GEO FCBs. This strategy can effectively avoid ISTB problems, datum inconsistency introduced by the code variations correction model, and low precise GEO products, which can improve the accuracy of the BDS FCB. Meanwhile, we will call the unclassified estimation [14] the traditional strategy to avoid confusion. In order to verify the validity of the new strategy in this study, the BDS WL FCBs and NL FCBs are estimated by the new strategy and traditional strategy respectively. Then the results of the two strategies are compared and analyzed in detail. The data are collected by 103 stations from the MGEX network, on DOY 061-090, 2018. The distribution of the 103 MGEX stations is shown in Figure 8.

3.3.1. BDS WL FCB

The BDS WL FCBs are calculated by a new strategy and traditional strategy based on the stations in Figure 8, with an interval of one day. The time series of WL FCBs are shown in Figure 9 and Figure 10. To eliminate FCBs datum, C02 and C08 are selected as reference satellites for GEO and non-GEO, respectively.
It can be seen from Figure 9 that within one month, the WL FCBs of IGSO and MEO estimated by the traditional strategy, except for C13, are relatively stable, and the fluctuation range is less than 0.2 cycle. The fluctuations are larger than 0.2 cycle for C13. On the other hand, GEO FCBs fluctuate greatly. From Figure 10, by the new strategy, the WL FCBs of IGSO and MEO are quite stable with fluctuations of less than 0.1 cycle. For GEO satellites, it is worth noting that the WL FCBs stability of C01, C03, C04 are comparable to that of IGSO and MEO in DOY 066-078, with fluctuations less than 0.1 cycle, but C05 still fluctuates greatly. For other time spans, GEO FCBs vary greatly. The reason for the instability may be the lack of IGS precision products for some GEO satellites during these days, resulting in insufficient GEO satellites that are actually used for FCBs estimation.
Meanwhile, Figure 11 and Figure 12 show the distribution of BDS WL fractional parts after the removal of FCBs by the traditional and new strategies, respectively. We find that compared to the traditional strategy, the ambiguity residuals of the new strategy are improved from 53.5% to 78.5%, and from 71.8% to 92.3% for <0.15 cycle and <0.25 cycle respectively for GEO, while the improvements for IGSO and MEO are very small. Furthermore, the precision of GEO WL FCBs could be comparable with IGSO and MEO in the new strategy. It shows that the new strategy can effectively improve the quality of BDS GEO FCBs.

3.3.2. BDS NL FCB

Because of the rapid variation of NL FCBs, it is estimated every 5 min during a day. It can be seen from Figure 13 and Figure 14 that the NL FCBs of IGSO and MEO are relatively stable except for C11, and the daily variation range of NL FCB is less than 0.2 cycles. The variation is very small for the adjacent 30–50 epochs (150–250 min), basically about 0.1 cycle or less. Compared to the IGSO and MEO, BDS GEO NL FCB is unstable, which is mainly due to the low precise IGS products for GEO. Moreover, the stability of BDS GEO NL FCB by the new strategy is still much better than that of the traditional strategy.
Figure 15 and Figure 16 show the distribution of BDS NL fractional parts after the removal of FCBs estimated by the traditional and the new strategy, respectively. Although it is not very obvious, there is still an improvement on the NL ambiguity residuals of IGSO and MEO by the new strategy. However, the NL ambiguity residuals of GEO by the new strategy are improved from 28.9% to 40.2%, and from 40.4% to 53.3% for <0.10 cycle and <0.15 cycle respectively, compared to the traditional strategy. The improvements of NL FCBs in the new strategy are not as obvious as that of WL FCBs. Besides, the precision of GEO NL FCBs is much worse than that of IGSO and MEO. It is mainly because of the low accuracy of current IGS products of GEO, which is worse than IGSO and MEO. What’s worse is that the wavelength of NL is only about 11cm, which is much shorter than WL, so GEO NL FCBs will be badly affected.

4. Experimental Results and Analysis

Using the new FCB strategy mentioned above, the BDS WL and NL FCBs were estimated to correct the fractional parts of BDS PPP ambiguities and obtain the fixed solutions. In this study, 10 stations distributed in the Asia-Pacific region (Figure 8) were selected to analyze and verify the accuracy of fixed solutions of the BDS PPP. The fastest time to first fix (TTFF) and positioning accuracy for PPP AR are analyzed. In this paper, convergence time is defined as the time required to achieve an accuracy of less than 5 cm, 5 cm, 10 cm for static PPP and less than 10 cm, 10 cm, 20 cm for Kinematic PPP in east, north and up directions, respectively. The TTFF is defined as the time taken for the first ambiguity to be successfully fixed [14]. The data duration is from DOY 061 to 090 in 2018, and the sample interval is 30 s. The estimated results are compared with the IGS weekly solutions.

4.1. Static PPP AR

Due to large number of stations and long observation time, this paper chose the KAT1 and XMIS stations to display the results at DOY 083 of 2018 which are shown in Figure 17, Figure 18 and Figure 19.
The result of the fixed solution converges rapidly in east, north, and up directions when the ambiguity is fixed. After convergence, the horizontal accuracy of the fixed solution is less than 1 cm, while that of the float solution is approximately 2 cm. The fixed solution in the up direction is approximately 3 cm, while that of the float solution is approximately 4 cm. In addition, it can be seen that there is a small systematic deviation in the up direction, which may be due to no accurate phase center offset (PCO) and phase center variation (PCV) corrections of receiver antenna for BDS currently.
Figure 20 and Table 5 show the results of BDS static PPP on DOY 061-091, 2018 for 10 selected stations in the Asia-Pacific region. As shown in Figure 20, the mean value of TTFFs for the 10 selected stations is about 40 min for static PPP AR. The convergence times are shortened for all the stations based on static PPP AR. The mean values of convergence time are 56.0 min, 43.6 min for float and fixed solution, respectively. The improvement is about 22.1%. Meanwhile, from Table 5 we can see that the RMS values of the fixed solution are less than 2 cm, 1 cm in the east, north directions. The accuracy of the north direction is the highest, followed by the east, while the up direction is the lowest for both float and fixed solutions. After the ambiguities are fixed, the accuracies are improved by about 16.7%, 23.0%, and 6.7% in the east, north, and up directions, respectively.

4.2. Kinematic PPP AR

The observations of 30 days from 10 stations in Asia-Pacific are used to evaluate the BDS kinematic PPP. The results of XMIS and KAT1 on DOY 083, 2018 are shown in Figure 21 and Figure 22.
The convergence time is clearly shortened when the ambiguities are fixed, especially in the east and north directions. The RMS of position error can be improved from 6 cm, 4 cm to 5 cm, 3.5 cm in the east and north directions, respectively. In the up direction, the accuracy of the fixed solution is slightly better than that of float solution after convergence, but the improvement is not obvious. Furthermore, the precision of the up direction is poor, and the jitter after convergence is significant, exceeding 20 cm at times. This may be due to no accurate PCO and PCV corrections of the receiver antenna for the BDS and the current configuration of the BDS constellation.
The results of BDS kinematic PPP during DOY 061-091, 2018 are shown in Figure 23 and Table 6 for 10 selected stations in the Asia-Pacific region. As shown in Figure 23, the average TTFF is about 58.7 min for kinematic PPP AR. Except for the station YAR3, the convergence time is shortened for all the other stations after kinematic PPP AR. The mean convergence times are 71.6 min and 62.7 min for float and fixed solution, respectively. The improvement is about 12.4%. From Table 6, it can be seen that the RMS values of the fixed solution are less than 5 cm, and 3.5 cm in the east, north directions. After the ambiguities are fixed, the accuracies are improved by about 18.1%, 21.2% and 7.9% in the east, north and up directions, respectively. The improvement of accuracy in the east and north directions is very obvious, while the improvement is slight in the up direction.

5. Conclusions

In this paper, several BDS biases are studied and analyzed in detail to enable ambiguity resolution. Research shows that there are half-cycle jumps between BDS GEO and non-GEO satellites for some receiver types, which results in inconsistencies in hardware delay between GEO and non-GEO for these receivers. Furthermore, there are serious satellite-induced code biases for the BDS satellites. The biases in IGSO and MEO can be corrected with existing models, while there is no correction model available to GEO.
In order to obtain high precision FCBs, this paper proposes a new strategy for estimating BDS GEO and non-GEO FCBs, respectively. This strategy could effectively avoid ISTB problems. Results show that the GEO and non-GEO FCBs estimated by the new strategy are more stable compared to the traditional strategy. The variation of GEO FCBs is less than 0.3 cycle in one month, except for C05. The FCBs of non-GEO satellites in one month have fluctuated 0.1 cycle. In addition, compared to the traditional strategy, the fractions of GEO WL ambiguities after removal of FCBs estimated by new strategy can be obviously improved from 53.5%, 71.8% to 78.5%, 92.3% for <0.15 cycle and <0.25 cycle respectively, which could be comparable with non-GEO. Simultaneously, the improvement of the GEO NL ambiguities is from 28.9%, 40.4% to 40.2%, 53.3% for <0.10 cycle and <0.15 cycle, which is not very obvious. This is because the current IGS products for GEO are in low precision. However, the NL FCBs of the adjacent 10–20 epochs (50–100 min) are less than 0.1 cycle. The daily variation of IGSO and MEO NL FCBs is much smaller, less than 0.2 cycle, except for C13. The variation of the adjacent 30–50 epochs (150–250 min) is very small, being less than 0.1 cycle.
In this study, the fixed solution of the BDS PPP is analyzed by using 10 stations distributed in the Asia-Pacific region. The results show that the convergence speed is faster when ambiguities are fixed. The accuracy and stability of static BDS PPP AR in the east, north, and up directions are improved compared with float solutions. The mean convergence times are shortened from 56.0 min to 43.6 min, and from 71.6 min to 62.7 min for static and kinematic PPP AR, respectively. The east and north directions are improved greatly, while the up direction is improved slightly. For static PPP, the position biases of ambiguity-fixed PPP are improved by about 16.7%, 23.0%, and 6.7% in the east, north, and up directions, respectively. Similarly, the improvements are approximately 18.1%, 21.2%, and 7.9% in the east, north, and up directions respectively, for BDS kinematic PPP.

Author Contributions

Conceptualization, W.J. and W.Z.; methodology, W.Z., X.L. and Q.C.; software, W.Z., H.C. and X.L.; validation, W.Z., X.L. and X.A.; formal analysis, W.Z.; investigation, W.Z.; resources, X.L.; data curation, X.L.; writing—original draft preparation, W.Z.; writing—review and editing, W.J. and H.C.; visualization, X.L. and X.A.; supervision, W.J.; project administration, Q.C.

Funding

This research is supported by the Major Technology Innovation Project of Hubei Province of China (2018AAA066), the National Science Foundation for Distinguished Young Scholars of China (Grant No.41525014), the Natural Science Innovation Group Foundation of China (No. 41721003), and the Major Project of Beijing Future Urban Design Innovation Center, Beijing University of Civil Engineering and Architecture (UDC 2018031321).

Acknowledgments

The authors gratefully acknowledge GFZ for providing GNSS precise orbit and clock products; many thanks to the IGS for providing the MGEX stations data and SNX products.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zumberge, J.F.; Heflin, M.B.; Jefferson, D.C.; Watkins, M.M.; Webb, F.H. Precise point positioning for the efficient and robust analysis of GPS data from large networks. J. Geophys. Res. Solid Earth 1997, 102, 5005–5017. [Google Scholar] [CrossRef] [Green Version]
  2. Blewitt, G. Carrier phase ambiguity resolution for the Global Positioning System applied to geodetic baselines up to 2000 km. J. Geophys. Res. Solid Earth 1989, 94, 10187–10203. [Google Scholar] [CrossRef] [Green Version]
  3. Ge, M.; Gendt, G.; Rothacher, M.; Shi, C.; Liu, J. Resolution of GPS carrier-phase ambiguities in Precise Point Positioning (PPP) with daily observations. J. Geod. 2008, 82, 389–399. [Google Scholar] [CrossRef]
  4. Laurichesse, D.; Mercier, F.; Berthias, J.P.; Broca, P.; Cerri, L. Integer Ambiguity Resolution on Undifferenced GPS Phase Measurements and Its Application to PPP and Satellite Precise Orbit Determination. Navigation 2009, 56, 135–149. [Google Scholar] [CrossRef]
  5. Collins, P.; Lahaye, F.; He´rous, P.; Bisnath, S. Precise point positioning with AR using the decoupled clock model. In Proceedings of the ION GNSS, Savannah, GA, USA, 16–19 September 2008; pp. 1315–1322. [Google Scholar]
  6. Li, X.; Zhang, X. Improving the estimation of uncalibrated fractional phase offsets for PPP ambiguity resolution. Navigation 2012, 65, 513–529. [Google Scholar] [CrossRef]
  7. Hauschild, A.; Montenbruck, O.; Sleewaegen, J.M.; Huisman, L.; Teunissen, P.J.G. Characterization of Compass M-1 signals. GPS Solut. 2012, 16, 117–126. [Google Scholar] [CrossRef]
  8. Montenbruck, O.; Hauschild, A.; Steigenberger, P.; Hugentobler, U.; Teunissen, P.J.G.; Nakamura, S. Initial assessment of the COMPASS/BeiDou-2 regional navigation satellite system. GPS Solut. 2013, 17, 211–222. [Google Scholar] [CrossRef]
  9. Wanninger, L.; Beer, S. BeiDou satellite-induced code pseudorange variations: Diagnosis and therapy. GPS Solut. 2015, 19, 639–648. [Google Scholar] [CrossRef]
  10. 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]
  11. Nadarajah, N.; Teunissen, P.J.G.; Raziq, N. BeiDou inter-satellitetype bias evaluation and calibration for mixed receiver attitude determination. Sensors 2013, 13, 9435–9463. [Google Scholar] [CrossRef] [PubMed]
  12. 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. 2015, 19, 357–368. [Google Scholar] [CrossRef]
  13. Hofmann-Wellenhof, B.; Lichtenegger, H.; Wasle, E. GNSS: Global Navigation Satellite Systems: GPS, Glonass, Galileo, and More; Springer: New York, NY, USA, 2008. [Google Scholar]
  14. Li, X.; Li, X.; Yuan, Y.; Zhang, K.; Zhang, X.; Wickert, J. Multi-GNSS phase delay estimation and PPP ambiguity resolution: GPS, BDS, GLONASS, Galileo. J. Geod. 2018, 92, 579–608. [Google Scholar] [CrossRef]
  15. 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. [Google Scholar]
  16. Wübbena, G. Software developments for geodetic positioning with GPS using TI-4100 code and carrier measurements. In Proceedings of the First International Symposium on Precise Posi-Tioning with the Global Positioning System, Rockville, MD, USA, 15–19 April 1985. [Google Scholar]
Figure 1. The fractional parts of WL and NL DD ambiguities for GEO and non-GEO satellites at zero-baseline ZIM2-ZIM3.
Figure 1. The fractional parts of WL and NL DD ambiguities for GEO and non-GEO satellites at zero-baseline ZIM2-ZIM3.
Sensors 19 04725 g001
Figure 2. The fractional parts of WL and NL DD ambiguities for GEO and non-GEO Satellites at zero-baseline CUCC-CUTC.
Figure 2. The fractional parts of WL and NL DD ambiguities for GEO and non-GEO Satellites at zero-baseline CUCC-CUTC.
Sensors 19 04725 g002
Figure 3. The M-W variation of C01, C02 at XMIS, on DOY 083, 2018.
Figure 3. The M-W variation of C01, C02 at XMIS, on DOY 083, 2018.
Sensors 19 04725 g003
Figure 4. The M-W variation of C06 at XMIS, on DOY 083, 2018.
Figure 4. The M-W variation of C06 at XMIS, on DOY 083, 2018.
Sensors 19 04725 g004
Figure 5. The M-W variation of C11 at XMIS, on DOY 083, 2018.
Figure 5. The M-W variation of C11 at XMIS, on DOY 083, 2018.
Sensors 19 04725 g005
Figure 6. The corrected M-W variation of C06 at XMIS, on DOY 083, 2018.
Figure 6. The corrected M-W variation of C06 at XMIS, on DOY 083, 2018.
Sensors 19 04725 g006
Figure 7. The corrected M-W variation of C11 at XMIS, on DOY 083, 2018.
Figure 7. The corrected M-W variation of C11 at XMIS, on DOY 083, 2018.
Sensors 19 04725 g007
Figure 8. Distribution of MGEX network stations to calculate BDS FCB. The green stars represent the station selected to evaluate the accuracy.
Figure 8. Distribution of MGEX network stations to calculate BDS FCB. The green stars represent the station selected to evaluate the accuracy.
Sensors 19 04725 g008
Figure 9. The variation of BDS WL FCB by the traditional strategy from DOY 061 to 090, 2018.
Figure 9. The variation of BDS WL FCB by the traditional strategy from DOY 061 to 090, 2018.
Sensors 19 04725 g009
Figure 10. The variation of BDS WL FCB by the new strategy from DOY 061 to 090, 2018.
Figure 10. The variation of BDS WL FCB by the new strategy from DOY 061 to 090, 2018.
Sensors 19 04725 g010
Figure 11. The distribution of BDS WL fractional parts after removal of FCBs estimated by the traditional strategy for 103 MGEX stations during DOY 061-090, 2018.
Figure 11. The distribution of BDS WL fractional parts after removal of FCBs estimated by the traditional strategy for 103 MGEX stations during DOY 061-090, 2018.
Sensors 19 04725 g011
Figure 12. The distribution of BDS WL fractional parts after removal of FCBs estimated by the new strategy for 103 MGEX stations during DOY 061-090, 2018.
Figure 12. The distribution of BDS WL fractional parts after removal of FCBs estimated by the new strategy for 103 MGEX stations during DOY 061-090, 2018.
Sensors 19 04725 g012
Figure 13. The variation of BDS NL FCB by the traditional strategy on DOY 083, 2018.
Figure 13. The variation of BDS NL FCB by the traditional strategy on DOY 083, 2018.
Sensors 19 04725 g013
Figure 14. The variation of BDS NL FCB by the new strategy on DOY 083, 2018.
Figure 14. The variation of BDS NL FCB by the new strategy on DOY 083, 2018.
Sensors 19 04725 g014
Figure 15. The distribution of BDS NL fractional parts after removal of FCBs estimated by the traditional strategy for 103 MGEX stations during DOY 061-090, 2018.
Figure 15. The distribution of BDS NL fractional parts after removal of FCBs estimated by the traditional strategy for 103 MGEX stations during DOY 061-090, 2018.
Sensors 19 04725 g015
Figure 16. The distribution of BDS NL fractional parts after removal of FCBs estimated by the new strategy for 103 MGEX stations during DOY 061-090, 2018.
Figure 16. The distribution of BDS NL fractional parts after removal of FCBs estimated by the new strategy for 103 MGEX stations during DOY 061-090, 2018.
Sensors 19 04725 g016
Figure 17. The skyplot of station: (a) KAT1; (b) XMIS.
Figure 17. The skyplot of station: (a) KAT1; (b) XMIS.
Sensors 19 04725 g017
Figure 18. The results of BDS static PPP AR at KAT1, on DOY 083 of 2018.
Figure 18. The results of BDS static PPP AR at KAT1, on DOY 083 of 2018.
Sensors 19 04725 g018
Figure 19. The results of BDS static PPP AR at XMIS, DOY 083 of 2018.
Figure 19. The results of BDS static PPP AR at XMIS, DOY 083 of 2018.
Sensors 19 04725 g019
Figure 20. Average convergence time of float and fixed solutions, and TTFFs of 10 selected stations in static PPP processing from DOY 061 to 090, 2018.
Figure 20. Average convergence time of float and fixed solutions, and TTFFs of 10 selected stations in static PPP processing from DOY 061 to 090, 2018.
Sensors 19 04725 g020
Figure 21. The results of BDS kinematic PPP AR at KAT1, on DOY 083 of 2018.
Figure 21. The results of BDS kinematic PPP AR at KAT1, on DOY 083 of 2018.
Sensors 19 04725 g021
Figure 22. The results of BDS kinematic PPP AR at XMIS, DOY 083 of 2018.
Figure 22. The results of BDS kinematic PPP AR at XMIS, DOY 083 of 2018.
Sensors 19 04725 g022
Figure 23. Average convergence time of float and fixed solutions, and TTFFs of 10 selected stations in kinematic PPP processing from DOY 061 to 090, 2018.
Figure 23. Average convergence time of float and fixed solutions, and TTFFs of 10 selected stations in kinematic PPP processing from DOY 061 to 090, 2018.
Sensors 19 04725 g023
Table 1. Receiver information of ZIM2 and ZIM3 on DOY 150, 2015.
Table 1. Receiver information of ZIM2 and ZIM3 on DOY 150, 2015.
StationReceiver TypeVersion
ZIM2TRIMBLE NETR94.85
ZIM3TRIMBLE NETR94.93
Table 2. Receiver information of CUCC and CUTC on DOY 151, 2018.
Table 2. Receiver information of CUCC and CUTC on DOY 151, 2018.
StationReceiver TypeVersion
CUCCJAVAD TRE_G3T DELTA3.73
CUTCTRIMBLE NETR95.30
Table 3. Frequency 1 with half-cycle jump (Unit: cycle).
Table 3. Frequency 1 with half-cycle jump (Unit: cycle).
I S T B r 1 r 2 , 1 β 1 β 2 I S T B r 1 r 2 , 2 β 1 β 2 Δ I S T B , w l i n t ( Δ I S T B , w l ) Δ I S T B , i f Δ I S T B , n l
0.500.50 or 1 f 1 2 2 ( f 1 2 f 2 2 ) f 1 2 ( f 1 f 2 ) o r f 1 2 f 2 2 ( f 1 f 2 )
−0.50−0.50 or −1 f 1 2 2 ( f 1 2 f 2 2 ) f 1 2 ( f 1 f 2 ) o r f 1 + 2 f 2 2 ( f 1 f 2 )
Table 4. Frequency 1 and 2 all with half-cycle jump (Unit: cycle).
Table 4. Frequency 1 and 2 all with half-cycle jump (Unit: cycle).
I S T B r 1 r 2 , 1 β 1 β 2 I S T B r 1 r 2 , 2 β 1 β 2 Δ I S T B , w l i n t ( Δ I S T B , w l ) Δ I S T B , i f Δ I S T B , n l
0.50.500 f 1 2 ( f 1 + f 2 ) 0.5
0.5−0.511 f 1 2 ( f 1 f 2 ) 0.5
−0.50.5−1−1 f 1 2 ( f 1 f 2 ) −0.5
−0.5−0.500 f 1 2 ( f 1 + f 2 ) −0.5
Table 5. The average accuracy of BDS static PPP for 10 stations (Unit: cm).
Table 5. The average accuracy of BDS static PPP for 10 stations (Unit: cm).
DirectionFixFloatAccuracy Improvement (%)
E1.792.1516.7
N0.971.2623.0
U5.816.236.7
Table 6. The average accuracy of BDS kinematic PPP for 10 stations (Unit: cm).
Table 6. The average accuracy of BDS kinematic PPP for 10 stations (Unit: cm).
DirectionFixFloatAccuracy Improvement (%)
E4.805.8618.1
N3.133.9721.2
U12.1713.217.9

Share and Cite

MDPI and ACS Style

Jiang, W.; Zhao, W.; Chen, H.; Liu, X.; An, X.; Chen, Q. Analysis of BDS Fractional Cycle Biases and PPP Ambiguity Resolution. Sensors 2019, 19, 4725. https://doi.org/10.3390/s19214725

AMA Style

Jiang W, Zhao W, Chen H, Liu X, An X, Chen Q. Analysis of BDS Fractional Cycle Biases and PPP Ambiguity Resolution. Sensors. 2019; 19(21):4725. https://doi.org/10.3390/s19214725

Chicago/Turabian Style

Jiang, Weiping, Wen Zhao, Hua Chen, Xuexi Liu, Xiangdong An, and Qusen Chen. 2019. "Analysis of BDS Fractional Cycle Biases and PPP Ambiguity Resolution" Sensors 19, no. 21: 4725. https://doi.org/10.3390/s19214725

APA Style

Jiang, W., Zhao, W., Chen, H., Liu, X., An, X., & Chen, Q. (2019). Analysis of BDS Fractional Cycle Biases and PPP Ambiguity Resolution. Sensors, 19(21), 4725. https://doi.org/10.3390/s19214725

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