Next Article in Journal
A Multi-Robot Coverage Path Planning Method for Maritime Search and Rescue Using Multiple AUVs
Previous Article in Journal
AutoML-Based Neural Architecture Search for Object Recognition in Satellite Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Various Analysis Center Products to Assess the Time-Frequency Transfer Performance of GPS/Galileo/BDS PPPAR Methods

1
School of Electronics and Communication Engineering, Sun Yat-Sen University, Guangzhou 510006, China
2
Shenzhen Key Laboratory of Navigation and Communication Integration, Shenzhen 518125, China
3
School of System Science and Engineering, Sun Yat-Sen University, Guangzhou 510006, China
4
College of Electronic Engineering, National University of Defense Technology, Hefei 230037, China
*
Author to whom correspondence should be addressed.
Submission received: 23 November 2022 / Revised: 16 December 2022 / Accepted: 20 December 2022 / Published: 24 December 2022

Abstract

:
Numerous organizations and Analysis Centers (AC) currently offer various Ambiguity Resolution (AR) products using various methodologies. However, there are no associated studies on their use for time-frequency transfer. This paper examines 16 Multi-GNSS Experiment (MGEX) stations with external high-precision atomic clocks to constitute 15 international time comparison links, and uses AR products data from CNES, SGG, CODE, and PRIDE laboratories, using three ambiguity-fixed strategies, to thoroughly evaluate the effects of various strategies and AR products for high-precision time-frequency transfer. We reach the following results by using the IGS final clock product as a reference and comparing it to ambiguity-float. With various ambiguity-fixed procedures, the time stability Standard Deviation (STD) of time transfer is increased for a single GPS, and the improvement ranges from 10 to 40%. The frequency stability has barely improved; up to 40%, the most notable improvement comes from FCB with GRM products. The time stability STD of combinations has improved after the addition of the Galileo system compared to the single GPS, and the improvement ranges from 2 to 9%. Most strategies have been improved, while a few techniques have been weakened with the GEC (GPS + Galileo + BDS) combination. We feel that the stability has not significantly increased with the systems’ increase in terms of short-term stability after comparing multiple groups of linkages.

1. Introduction

Common View (CV) [1,2,3], All in View (AV) [4,5], Precise Point Positioning (PPP) [6], and Two-Way Satellite Time and Frequency Transfer (TWSTFT) [7] are some of the current GNSS time transfer techniques. Since September 2009, GPS PPP has been a crucial approach in the International Atomic Time (TAI) time comparison [6]. The GNSS PPP technique has the advantages of globalization, all-weather functionality, high precision, and low cost, and its short-term stability is superior to that of TWSTFT. Its time transfer accuracy is also comparable to that of TWSTFT. One of the current hot topics in high-precision time transfer technology is the GNSS PPP approach.
Time-frequency transfer is done via pseudo-range and carrier phase observations in GNSS PPP approaches [8]. Better short-term stability can be attained compared to pseudo-range only techniques because phase observations have lower noise levels and higher measurement precision [9]. According to pertinent studies, the normal uncertainty of PPP frequency comparisons has an average day stability of roughly 1 × 10−15 and stability of 1 × 10−16 over a 30-day period. The Type A uncertainty of Circular T for the time comparison link of the BIPM is less than 0.3 ns [6].
The PPP approach was first employed for high-precision positioning [10] before being widely applied in the monitoring of crustal deformation [11], meteorology [12], high precision dynamic positioning [13], and monitoring of local seismic activity [14].
Recent advances in PPP ambiguity-fixed technology have significantly increased convergence speed, positioning accuracy, and stability. Scholars and research organizations both domestically and abroad have given it considerable attention. The ambiguity-fixed solution PPP has distinct advantages over the ambiguity-float solution PPP, because when the ambiguity converges or has a certain accuracy, it fully utilizes the integer advantage of ambiguity to convert the phase observation value that contains ambiguity into a high-precision absolute distance observation value, which quickly increases the positioning accuracy to a stable centimeter level. Because they do not adequately handle the deviation terms contained in the observed values from the model, much earlier research has concentrated on ambiguous float solution PPP. These deviation terms contaminate the ambiguity and render it unfixable. Although numerous PPP ambiguity-fixed approaches have been devised, their overarching objective is reinstating ambiguity’s fixable (integer) characteristic.
Three representative PPP ambiguity-fixed methods have been proposed: the integer phase clock method by Laurichesse et al. of the CNES (Centre National d’Etudes Spatials) in France [15], the decimal deviation method by Ge et al. of the GFZ (Geo Forschungs Zentrum) in Germany [16], and the clock difference decoupling method by Collins et al. of NRCan (Natural Resources Canada). Domestic researchers have noted that these methods are equivalent from a variety of perspectives but have also noted that they differ in the final products and algorithm implementation [17]. The observable-specific signal bias (OSBs) in Solution INdependent EXchange (SINEX) format, which includes the observable specific satellite code and phase bias corrections for each frequency, can be directly applied to the raw observations before forming any linear combination. The three methods are theoretically equivalent in positioning [18,19].
Some institutions are currently able to offer PPP AR products for the GPS, Galileo, and BDS systems. The benefits of various institutional products and various system products are therefore worth additional comparison and analysis, and then further examination of the impact of their applications, to understand the multi-GNSS (GPS, Galileo, BDS) PPP AR (whether positioning or time transfer applications). There is currently no research on time-frequency transfer, and the performance of various PPP AR products from various institutions has only been studied in terms of locating applications [20,21].
More than 80% of independent ambiguities can be consistently addressed, according to experimental findings in the literature [15], and the PPP single-day static solution with fixed ambiguity in the east direction has a positioning accuracy of about 30% greater than the PPP’s actual solution. The U direction and the clock difference parameter are correlated, as demonstrated in related research [22], eliminating any remaining uncertainty. This strengthens the position of the U direction and, in turn, the clock difference parameter.
In less than a week of solution time, the first PPPAR experiment for time-frequency transfer demonstrated that frequency stability of 1 × 10−16 was possible with PPPAR technology [9]. The GRE PPPAR mode watches the clock when RT-IPPP technology is used, detecting changes with a frequency larger than 6.1 × 10−15 [23]. Their effort is limited to a single source AR product and a single ambiguity-fixing method.
This work performs high-precision time-frequency transfer experiments with AR products from CNES, SGG, CODE, and PRIDE laboratories using three ambiguity-fixed techniques. A thorough investigation of the time-frequency transfer theory and associated GPS/Galileo/BDS PPPAR technique results for various analysis center products was carried out using the IGS 30-s final clock difference product as a reference.
The structure of the article is as follows: The PPPAR time-transfer principle and function model are initially introduced; the experimental data and comprehensive processing schemes are then presented; followed by comparison and demonstration of the utilization of several techniques when in static mode. Performance outcomes for numerous goods using multi-GNSS PPP AR time-frequency transmission are observed; finally, a summary is given and conclusions are drawn.

2. Materials and Methods

2.1. PPPAR Time Transfer

The stations are typically coupled to high-precision time-frequency reference signals (1PPS and 10MHZ) for the high-precision time-frequency transfer of the PPPAR approach, such as stations A and B in Figure 1. The station leverages precision products produced outside to collect pseudo-range and carrier-phase observation data through the GNSS antenna (precision orbit, clock difference, DCB, AR products, etc.).
Firstly, determination of the difference between the clocks CLK(a) and CLK(b) connected to the receiver and the standard reference time (such as GPST, IGST, or UTC (NTSC)), performing the difference operation to obtain that difference, and then passing the calibration method (https://webtai.bipm.org/ftp/pub/tai/publication/gnss-calibration/guidelines (accessed on 1 October 2022)) to complete calibration of the hardware time delay. This paper’s clock difference is the outcome of ambiguity-fixed technology, which typically uses PPP technology. The alternative approach is as follows:
CLK(a) − CLK(b) = (CLK(a) − Ref) − (CLK(b) − Ref)

2.2. GNSS-PPPAR Model

The Equation (1) examples illustrate the main observations of GNSS raw pseudo-range and phase in Equation (2):
{ P r , j = ρ r + d t r d t s + m w , r T w , r + γ j I r , 1 + b r , j + b j s + ε P j L r , j = ρ r + d t r d t s + m w , r T w , r + γ j I r , 1 + λ j N r , j + B r , j + B j + ε L j
In Equation (2), P represents pseudo-range observations; L represents carrier-phase observations; r and j represents receiver and signal frequency, respectively; ρ represents the geometric distance between the satellite and the receiver; T w , r represents the tropospheric Zenith Wet Delay (ZWD); I represents ionospheric delay; γ j = ( f 1 / f j ) 2 and f j represents frequency j coefficient; d t r and d t s represent receiver and satellite clock difference, respectively; m w , r represents map projection function; N r , j is the integer phase ambiguity; λ j means wavelength in frequency j; ε P j s and ε L j s are the pseudo-range and carrier-phase error term which did not quantize into the function model, respectively.
In Equation (2), b r , j and B r , j are the Uncalibrated Code Delay (UCD) and Uncalibrated Phase Delay (UPD) at the receiver; b j s and B j s are UCD and UPD at the satellite end. The general function model is the combination of Ionospheric-Free (IF). We define the coefficient as follows in Equation (3):
α i j = ( f i ) 2 ( f i ) 2 ( f j ) 2 , β i j = ( f j ) 2 ( f i ) 2 ( f j ) 2
We assume Equation (4):
{ b r , j = Δ b r , j + ε b r , j b j s = Δ b j s + ε b j s B r , j = Δ B r , j + ε B r , j B j s = Δ B j s + ε B j s
in Equation (4), Δ and ε represent constant and variable parts of hardware delay, respectively [24]. IF combined pseudo-range and carrier observations can be written as follows in Equations (5)–(7):
P r , I F = α 12 P r , 1 + β 12 P r , 2 = ρ r + d t r d t s + m w , r T w , r + b r , I F + b I F s + ε P I F
L r , I F = α 12 L r , 1 + β 12 L r , 2 = ρ r + d t r d t s + m w , r T w , r + N r , I F s + B r , I F + B I F s + ε L I F
N r , I F s = λ N L N r , 1 s λ 2 γ 2 1 N r , W L s
In Equation (7), N r , W L is the Wide-Lane (WL) ambiguity and λ N L is the wavelength of the narrow-lane (NL) ambiguity. N r , 1 s is called the NL ambiguity. Hatch–Melbourne–Wubbena (HMW) combination can be expressed as follows in Equation (8):
K M W = λ W L ( L 1 λ 1 L 2 λ 2 ) f 1 P 1 + f 2 P 2 f 1 + f 2 = λ W L N r , W L s + τ r , W L + τ W L s
In Equation (8), λ W L is the WL wavelength. To eliminate the hardware delay at the receiver, PPPAR generally uses the method of an inter-satellite single difference (SD) method. The ambiguity of SD IF combination and WL ambiguity is expressed as follows in Equation (9):
{ N r , I F i j = N r , l F i N r , I F j N r , W L i j = N r , W L i N r , W L j
The IGS satellite clock difference product absorbs the pseudo-range hardware delay. The satellite clock difference can be expressed as follows in Equation (10):
{ d t P I F s = d t s Δ b I F s ε B I F s b I F s = α 12 b 1 s + β 12 b 2 s
In Equation (10), d t P I F s represents satellite clock difference products, and d t s represents the “real” clock difference of the satellite. To make use of this IF satellite clock product, we rewrite Equation (10) as Equation (11):
{ d t s + Δ b 1 s = d t P I F s + β 12 D C B 12 s d t s + Δ b 2 s = d t P I F s α 12 D C B 12 s
In Equation (11), D C B i j s is the DCB at the satellite end, which is Δ b i s Δ b j s . We substitute clock difference to obtain Equation (12):
{ P I F = ρ r + d t r , P I F + m w , r T w , r + Ψ + ε P I F L I F = ρ r + d t r , P I F + m w , r T w , r + N r , l F s + Δ B r , I F ε b r , l F + Δ B I F s + Δ b I F s + ε L I F
And we assume Equation (13):
{ d t r , P I F s = d t r s + Δ b r , I F s + ε B r , I F s Ψ = ε b r , l F s ε B r , I F s + ε b I F s + ε B I F s N r , I F s = N r , l F s + Δ B I F s + Δ b I F s + Δ B r , l F s Δ b r , I F s
In Equation (13), Ψ will be submerged by pseudo-range noise and finally enter the pseudo-range residual. The IF combination’s ambiguity is divided into a linear combination of WL and NL in Equation (7). The WL is fixed first, followed by the NL, in the fixed ambiguity sequence. Following are the fixed stages [25]:
Firstly, determine the SD IF float ambiguity. Calculate the ambiguity parameters and the accompanying Variance-Covariance Matrix (VCM) by applying the error propagation approach to the standard extended Kalman filter (EKF) calculation in Equation (14):
{ N r , I F i j = N r , I F i N r , I F j = N r , I F i j + [ Δ B I F i j ] + [ Δ b I F i ] Q N r , I F i j = ( H S D Q N r , I F H S D T ) i j
In Equation (14), N r , I F i j is a SD IF solution. H S D is SD design matrix and Q N r , I F is IF ambiguity VCM [·] indicates that this term can be eliminated by using PPPAR products.
Secondly, fixed SD WL ambiguity. The MW combination can be used to calculate the float solution to the WL ambiguity. The (OBS or FCB) product can then be used to obtain the integer solution. In general, it can also be rounded up directly in Equation (15).
N ^ r , W L i j = r o u n d ( ( Ψ ¯ m w i j + [ τ W L i j ] ) / λ W L )
In Equation (15), N ^ r , W L i j denotes fixed SD WL ambiguity and represents the average MW combination; the WL UPD τ W L i j is dispensable and can be eliminated using PPPAR products.
Thirdly, the SD NL ambiguity was fixed. Following a successful WL ambiguity correction, the expression equation for the SD NL term can be written as follows in Equation (16):
{ N ˜ r , 1 i j = ( N ˜ r , I F i j + λ 2 γ 2 1 N r , W L i j + [ τ N L i j ] ) / λ N L Q N ˜ r , 1 i j = ( 1 ( λ N L ) 2 Q N ˜ r , I F ) i j
In Equation (16), Q N ˜ r , 1 i j is the SD float NL ambiguities which are obtained by applying the covariance propagation law. The NL UPD correction τ N L i j is required for the FCB method. Since NL ambiguities are strongly correlated in PPP, the popular Least-squares Ambiguity Decorrelation Adjustment (LAMBDA) method or a modified variant serves as the foundation of a search strategy. The widely used ratio test and success rate are used to decide what must be fixed.
Fourthly, Fixed SD IF ambiguity. After the SD WL and NL are successfully fixed, the SD IF can be recalculated to obtain in Equation (17):
N ^ r , I F i j = λ N L s ( N ^ r , 1 i j [ τ N L i j ] ) λ 2 s γ 2 s 1 N ^ r , W L i j
In Equation (17), NL UPD correction τ N L i j is required for the FCB method.
Lastly, the fixed solution is updated. Other parameters can be updated by their correlation with fixed ambiguities, such as Position, ZWD, the clock difference parameters, and the remaining unfixed ambiguities addressed in Equation (18):
b ^ = b ˜ Q b ^ N ^ Q N ˜ 1 ( N ^ N ˜ )
In Equation (18), b ^ and b ˜ are fixed and unfixed parameters, respectively, N ˜ is float solution with VCM Q N ^ 1 , N ^ is the integer ambiguity, Q b ^ N ^ is the covariance matrix of b ^ and N ^ .
The receiver clock difference parameters in this paper are updated using Equations (18) and (19).
b ^ = ( p o s , c l k , Z W D , N u n f i x )
Equation (19) contains the position p o s , zenith tropospheric delay Z W D , the clock difference parameters c l k , and the remaining unfixed ambiguity N u n f i x . c l k is the clock difference solved in this paper, made up of the parameters after the ambiguity is fixed.
In order to create 15 international time comparison links, this paper chooses 16 MGEX stations with external high-precision atomic clocks, uses AR product data from CNES, SGG, CODE, and PRIDE laboratories, and conducts high-precision time-frequency transfer experiments using three ambiguity-fixed strategies. Using the IGS final clock product as a benchmark and contrasting it with ambiguity-float.

3. Experimental Data and Processing Strategies

The BIPM (Bureau International des Poids et Measures) was founded in 1875. The information from more than 500 atomic clocks provided by BIPM to time laboratories across the world is weighted to produce the International Atomic Time (TAI). The PTB (Physikalisch-Technische Bundesanstalt) in Germany serves as the key node for the current international time comparison [26].
The chosen stations are primarily based on the following criteria: they must be participating stations in the BIPM international time comparison time laboratory, have a receiver attached to an ultra-reliable atomic clock, and have at least GPS, Galileo, and BDS system observations. In order to expand the scope of the test, we also selected some stations that did not participate in the international TAI comparison and did not have an external atomic clock but had an external high-precision time and frequency source.
Figure 2 illustrates the selection of experimental observation data from 16 globally dispersed sites. The MGEX FTP server served as the source of the observation data. Table 1 displays the specific setup information for these stations, including antennas and the different kinds of external atomic clocks.
The publicly accessible PPPAR solutions produced by several agencies with three preset strategies are listed in Table 2. For the sake of conciseness, the terms “GRM”, “WUM”, “COM”, and “WHU” in the remaining sections of this article, respectively, stand for PPP AR products from CNES, SGG-WHU, CODE, and PRIDE laboratories. BDS-2 and QZSS FCB products are provided by SGG-WHU, and PRIDE laboratories supports BDS-2 and BDS-3 satellites for PPPAR [21]. This article’s data is based on information as of January 1, 2020. These three fixed techniques are referred to as FCB, IRC, and observable-specific code biases (OSBs) in the article that follows.

4. PPP Detailed Processing Strategy

The detailed PPPAR strategy of the article is shown in Table 3 below.
It should be noted that this study employs the Kalman filtering technique. As a result, the PPP outcomes at the initialization stage were not significant. The PPP results are smoothed in this research by combining forward and reverse filtering with smoothing. Three steps make up the process: forward PPP filtering results, reverse PPP filtering results, and combining the first two results using the weighted average approach. This smoothing equation has the following in Equation (20):
Q k r e s u l t = ( ( Q k f ) 1 + ( Q k b ) 1 ) 1 X k r e s u l t = Q k r e s u l t ( ( Q k f ) 1 X k f + ( Q k b ) 1 X k b ) = X k f + Q k r e s u l t ( Q k f ) 1 ( X k b X k f )
In Equation (20), Q k f and X k f are the variance matrix and state vector at time k, respectively. Q k b and X k b correspond to the forward filtering result at time k, and are, respectively, the variance matrix and state vector corresponding to the reverse filtering result at time k; Q k r e s u l t and X k r e s u l t are the variance matrix and state vector corresponding to the smoothed result at time k, respectively. The convergence time issue in the first stage of PPP or PPPAR is eliminated in this research using the bidirectional Kalman filter approach.

5. Validation and Analysis

5.1. Test Case

The ambiguity-fixed methods are FCB, IRC, and OSB. Due to the existence of precision orbit and clock products with various analysis centers, we use three fixed methods (FCB, IRC, and OSB) to create a system of six alternative ways to study our results. The following six categories for numerous items are determined by the three approaches (FCB, IRC, and OSB), as indicated in Table 4.
There is no convergence issue with the PPP approach because we employ the bidirectional filtering method for better analysis and comparison. The positioning outcome of the PTBB station is depicted in Figure 3. It is obvious that there is no convergence issue and that the convergence effect may always be obtained. It consistently maintains a 1 cm accuracy in three directions (NEU).
Using various analysis center products, some academics have carried out positioning work [20,21]. To help with a better grasp of the overall logic, we also do the positioning results here. Three directions are divided (NEU), and the RMS display is displayed in Figure 4. We can observe that the ambiguity-fixed has a positive effect on the positioning if we use the SNX weekly solution as a point of reference. After the ambiguity is resolved, nearly all solutions produce superior positioning outcomes. The performance of GBM is the finest among them and the localization outcomes are exceptional.
We are aware that there is a strong association between the U direction and clock difference [22]. We assume that the clock difference operates more effectively because the U direction is enhanced. On the basis of this concept, we continue to solve the clock difference result after the ambiguity has been resolved. With 16 stations, it is clear that the ambiguity-fixed method is quite effective at containing those diverging points. In Figure 5, we display PTBB station clock differences of the three techniques (FCB, IRC, and OSB).
No matter which fixing technique is employed, all of them are identical to the location and may successfully revert some points of clock discrepancy, as illustrated in Figure 5. We conducted further investigations to determine the magnitude of the impact effect.
We measured the STD box display of the 16 stations by comparing it to the final clock difference of the IGS week SNX file. Figure 6 illustrates how we confirmed our hypothesis that, with the exception of OSB_ com’s fixed strategy, all other fixed strategies have lower STDs for fixed clock difference than for the float solution.
As the outcome of the clock difference improves, the time transfer also does, and we can now move on to the study of time linkages made up of two stations.
We create the 15 worldwide time comparison connections mentioned above using the PTBB station as the basis. Due to space restrictions, Figure 7 shows the sequence diagrams for the two time-links, PTBB−MGUE and PTBB−USN7. It is clear that the fixed approach, particularly when using WHU goods, has a particular smoothing effect on unfixed.
We estimated the STD of the time transfer to assess the Type A uncertainty in the remote time-transfer link using the time-transfer results and calculation of station clock difference acquired from the IGS final product as a reference. For a single GPS, the time stability STD of time transfer is enhanced using various ambiguity-fixing procedures, and the improvement ranges from 10 to 40% in the Figure 8 STD box display of 15 time-links. The one with the most improvement, up to 40%, is FCB_grm. OSB_com is the unconventional strategy. In the instance of a single station (Figure 6), one exception has already occurred, and a bigger exception follows the creation of the time connection.
The frequency stability of the time-links between the PPP and PPPAR results should then be examined. The results of PPP and PPPAR alter in the same manner for the transfer at the same moment (Figure 7). The PPPAR results’ variability is minimal in comparison. We examine time linkages’ frequency stability. To examine the frequency stability of time series, we employ the Allan variance modification technique and the following equations:
δ y 2 ( τ ) = 1 N 1 i = 1 N 1 ( x n ( k + 1 ) 1 2 x n k + 1 + x n ( k 1 ) + 1 ) 2
mod σ ^ y 2 ( n τ 0 ) = 1 2 n 4 τ 0 2 ( N 3 n + 1 ) j = 1 N 3 n + 1 [ i = j j + n 1 ( x i + 2 n 2 x i + n + x i ) ] 2
The above Equations (21) and (22), τ is the interval time, and the maximum can only be half of the number of data, x n representing the phase data.
Due to space constraints, we are only able to demonstrate the features of six common time-links (PTBB−BOR1, PTBB−GMSD, PTBB−NRC1, PTBB−OHI3, PTBB−TID1, PTBB−TIDV). Figure 9 makes it clear that the frequency stability of the time link has only slightly improved as a result of the ambiguity-fixed technique. We discovered that between 1000 and 10,000 s, OSB com fix, IRC_ grm_ float, and IRC grm fix all perform poorly.
Go one step further; in an ordinary clock time link, stations participating in time transfer share the same atomic clock. The relative clock skew from station to station will be eliminated. The standard can directly calculate the deviation to evaluate the time transfer performance. Taking TIDV-TID1 as an example, as shown in Figure 10, we can’t get the conclusion that its frequency stability of fixed ambiguity is better than that of the float in terms of short-term stability.

5.2. Muti-GNSS

Researchers have carried out pertinent experiments. The results of the time-frequency transfer are enhanced by the multi-GNSS PPP approach [30]. This section’s goal is to demonstrate this for the multi-GNSS PPPAR technique time transfer. When using multiple systems, we take GPS as the benchmark, and the system differences between other systems and GPS are estimated by parameters.
We developed six schemes and fifteen particular comparison methods, as indicated in Table 5, in accordance with various ambiguity-fixed strategies and the satellite system integrated into the product. Following a set strategy-based division into FCB, IRC, and OSB, the products are then further categorized into six schemes. Additionally, there are 15 precise fixed methods in the bundled Satellite system. GEC denotes the addition of BDS based on GPS and Galileo. G indicates for single GPS. GE stands for GPS plus Galileo.
With only a minimal addition of satellite systems, we used the same PPP ambiguity fixing approach as in the preceding section. In a similar manner, we solve 16 stations, create 15 time-links centered on PTBB, and determine the STD with respect to the clock difference of the IGS weekly solution stations. The outcomes are displayed in Figure 11.
As illustrated in Figure 11, with the exception of the COM product approach employing a single GPS, the statistical results of time link STD of PPPAR, whether single-system or multi-system, can be maintained at the sub-nanosecond level; the six GE combinations have improved compared to the single GPS after the addition of the Galileo system, and the improvement ranges from two to nine percent; the GEC combo has improved the majority of methods while weakening a few others. Inspection reveals that it is caused by some stations’ subpar BDS fixes.
Due to physical constraints, we can only provide the frequency stability data of the four distinct Muti-GNSS PPPAR temporal links in Figure 12. The integration of multiple systems (Galileo, BDS) causes a significant shift in frequency stability, which is slightly different from the single system GPS. Even so, some are weakened after being corrected, while some are completely unaffected, and some are even better. We believe that despite the system’s gain in short-term stability, the stability has not been greatly enhanced based on the comparison of multiple groups of linkages.

6. Conclusions

This paper examined 16 MGEX stations with external high-precision atomic clocks to form 15 time comparison links and uses AR products data from CNES, SGG, CODE, and PRIDE laboratories. It then used three ambiguity resolution strategies (FCB, IRC, OSB) to conduct high-precision time-frequency transfer experiments in order to thoroughly evaluate the impact of various strategies and AR products.
We arrived at the following findings by using the IGS end clock product as a reference and contrasting it with the float solution. With various ambiguity-fixing procedures, the time stability STD of time transfer is increased for a single GPS, and the improvement ranges from 10 to 40%. Of these, FCB_grm has undergone the greatest improvement, up to 40%; frequency stability has barely changed. When compared to a single GPS, combinations have improved with the addition of the Galileo system, and the improvements range from 2 to 9%. The majority of methods have improved, while a few techniques have weakened with the GEC combination. We feel that the increased short-term stability of the system has not significantly enhanced stability based on the comparison of multiple groups of links. This research provides an important reference for future PPPAR time-frequency transmission engineering applications.

Author Contributions

M.O. and X.Z. conceived and designed the experiments; M.O. performed the experiments, analyzed the data, and wrote the paper; X.Z., R.Y. and D.L. contributed to discussions and revisions. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (Grant No. 2021YFA0716500), the National Natural Science Foundation of China (Grant No. 61973328, 91938301), and the Key Basic Research Projects of Shenzhen Science and Technology Commission (Grant No. 2020N259), Shenzhen Science and Technology Program (Grant No. GXWD20201231165807008, 20200830225317001), and was also supported by Shenzhen Key Laboratory of Navigation and Communication Integration (Grant No. ZDSYS20210623091807023).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, H.; Yan, M.; Guang, W. Comparison of two-way and beidou common view time transfer between ntsc and beij station. Adv. Mater. Res. 2014, 1049–1050, 1997–2000. [Google Scholar] [CrossRef]
  2. Nawrocki, J.; Lewandowski, W.; Nogas, P.; Foks, A.; Lemanski, D. An experiment of GPS + GLONASS common-view time transfer using new multi-system receivers. In Proceedings of the Frequency and Time Forum (EFTF), 2006 20th European IEEE, Braunschweig, Germany, 27–30 March 2006. [Google Scholar]
  3. Allan, D.W.; Weiss, M.A. Accurate Time and Frequency Transfer During Common-View of a GPS Satellite. In Proceedings of the IEEE Symposium on Frequency Control, Vancouver, Canada, 29–31 August 2005. [Google Scholar]
  4. Harmegnies, A.; Defraigne, P.; Petit, G. Combining gps and glonass in all-in-view for time transfer. Metrologia 2013, 50, 277–287. [Google Scholar] [CrossRef]
  5. Petit, G.; Jiang, Z. Gps all in view time transfer for tai computation. Metrologia 2007, 45, 35. [Google Scholar] [CrossRef]
  6. Petit, G. The TAIPPP pilot experiment. In Proceedings of the IEEE International Frequency Control Symposium, 2009 Joint with the 22nd European Frequency and Time Forum, Besancon, France, 20–24 April 2009. [Google Scholar]
  7. Zhang, V.S.; Parker, T.E.; Jian, Y. Long-term uncertainty in time transfer using GPS and TWSTFT techniques. In Proceedings of the Frequency Control Symposium the European Frequency Time Forum, Denver, CO, USA, 12–16 April 2015. [Google Scholar]
  8. Ray, J.; Senior, K. Geodetic techniques for time and frequency comparisons using gps phase and code measurements. Metrologia 2005, 42, 215. [Google Scholar] [CrossRef] [Green Version]
  9. Petit, G.; Kanj, A.; Loyer, S.; Delporte, J.; Mercier, F.; Perosanz, F. 1 × 10−16 frequency transfer by gps ppp with integer ambiguity resolution. Metrologia 2015, 52, 301. [Google Scholar] [CrossRef]
  10. 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]
  11. Su, K.; Jin, S.; Ge, Y. Rapid displacement determination with a stand-alone multi-GNSS receiver: GPS, Beidou, GLONASS, and Galileo. GPS Solut. 2019, 23, 54. [Google Scholar] [CrossRef]
  12. Lu, C.; Li, X.; Nilsson, T.; Ning, T.; Heinkelmann, R.; Ge, M.; Glaser, S.; Schuh, H. Real-time retrieval of precipitable water vapor from gps and beidou observations. J. Geod. 2015, 89, 843–856. [Google Scholar] [CrossRef]
  13. Gao, Z.; Zhang, H.; Ge, M.; Niu, X.; Shen, W.; Wickert, J.; Schuh, H. Tightly coupled integration of multi-gnss ppp and mems inertial measurement unit data. Gps Solut. 2017, 21, 377–391. [Google Scholar] [CrossRef]
  14. Li, X.; Ge, M.; Guo, B.; Wickert, J.; Schuh, H. Temporal point positioning approach for real-time gnss seismology using a single receiver. Geophys. Res. Lett. 2013, 40, 5677–5682. [Google Scholar] [CrossRef]
  15. 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]
  16. 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]
  17. Collins, P.; Bisnath, S.; Lahaye, F.; Héroux, P. Undifferenced gps ambiguity resolution using the decoupled clock model and ambiguity datum fixing. Navigation 2010, 57, 123–135. [Google Scholar] [CrossRef] [Green Version]
  18. Shi, J.; Gao, Y. A comparison of three ppp integer ambiguity resolution methods. Gps Solut. 2014, 18, 519–528. [Google Scholar] [CrossRef]
  19. Geng, J.; Meng, R.; Dodson, R.H.; Teferle, R.N. Integer ambiguity resolution in precise point positioning: Method comparison. J. Geod. 2010, 84, 569–581. [Google Scholar] [CrossRef] [Green Version]
  20. Glaner, M.; Weber, R. Ppp with integer ambiguity resolution for gps and galileo using satellite products from different analysis centers. GPS Solut. 2021, 25, 102. [Google Scholar] [CrossRef]
  21. Chen, C.; Xiao, G.; Chang, G.; Xu, T.; Yang, L. Assessment of gps/galileo/bds precise point positioning with ambiguity resolution using products from different analysis centers. Remote Sens. 2021, 13, 3266. [Google Scholar] [CrossRef]
  22. Ouyang, M.; Li, J.; Li, W.; Ge, Y.; Dai, Z. Research on time and frequency transfer during ppp convergence with parameters correlation comparison. Measurement 2020, 173, 108597. [Google Scholar] [CrossRef]
  23. Lyu, D.; Zeng, F.; Ouyang, X.; Zhang, H. Real-time clock comparison and monitoring with multi-gnss precise point positioning: Gps, glonass and galileo sciencedirect. Adv. Space Res. 2020, 65, 560–571. [Google Scholar] [CrossRef]
  24. Zhao, L.; Ye, S.; Jia, S. Handling the satellite inter-frequency biases in triple-frequency observations. Adv. Space Res. 2017, 59, 2048–2057. [Google Scholar] [CrossRef]
  25. Liu, S.; Sun, F.; Zhang, L.; Li, W.; Zhu, X. Tight integration of ambiguity-fixed ppp and ins: Model description and initial results. GPS Solut. 2016, 20, 39–49. [Google Scholar] [CrossRef]
  26. Defraigne, P.; Aerts, W.; Pottiaux, E. Monitoring of utc(k)’s using ppp and igs real-time products. Gps Solut. 2015, 19, 165–172. [Google Scholar] [CrossRef]
  27. Sastamoinen, J. Atmospheric correction for troposphere and stratosphere in radio ranging of satellites, in the use of artifical satellites for geodesy. Geophys. Monogr. Ser. 1972, 15, 247–252. [Google Scholar]
  28. Kouba, J. Testing of global pressure/temperature (gpt) model and global mapping function (gmf) in gps analyses. J. Geod. 2009, 83, 199–208. [Google Scholar] [CrossRef]
  29. Petit, G.; Luzum, B. Iers Conventions IERS Technical Note; Bonifatius GmbH: Paderborn, Germany, 2010. [Google Scholar]
  30. Ge, Y.; Dai, P.; Qin, W.; Yang, X.; Zhou, F.; Wang, S.; Zhao, X. Performance of multi-gnss precise point positioning time and frequency transfer with clock modeling. Remote Sens. 2019, 11, 347. [Google Scholar] [CrossRef]
Figure 1. Sketch map for PPPAR time-frequency transfer.
Figure 1. Sketch map for PPPAR time-frequency transfer.
Remotesensing 15 00092 g001
Figure 2. The distribution map of 16 stations was selected, and each station is externally connected with a high-precision atomic clock.
Figure 2. The distribution map of 16 stations was selected, and each station is externally connected with a high-precision atomic clock.
Remotesensing 15 00092 g002
Figure 3. Position results of bidirectional filtering of PTBB station.
Figure 3. Position results of bidirectional filtering of PTBB station.
Remotesensing 15 00092 g003
Figure 4. 6 methods of 3 fixed strategies’ (FCB, IRC, OSB) positioning results, including three directions of NEU and RMS value composed of three directions.
Figure 4. 6 methods of 3 fixed strategies’ (FCB, IRC, OSB) positioning results, including three directions of NEU and RMS value composed of three directions.
Remotesensing 15 00092 g004
Figure 5. Three strategies of PTBB station results with fixed and floating ambiguity clock difference, where red represents floating point and green represents fixed solution.
Figure 5. Three strategies of PTBB station results with fixed and floating ambiguity clock difference, where red represents floating point and green represents fixed solution.
Remotesensing 15 00092 g005
Figure 6. The STD box display of the 16 stations we choose to take the final clock difference of the IGS week SNX file as a reference, red ”+” stands for outliers.
Figure 6. The STD box display of the 16 stations we choose to take the final clock difference of the IGS week SNX file as a reference, red ”+” stands for outliers.
Remotesensing 15 00092 g006
Figure 7. PTBB-MGUE/PTBB-USN7 Time links with different strategies.
Figure 7. PTBB-MGUE/PTBB-USN7 Time links with different strategies.
Remotesensing 15 00092 g007
Figure 8. The STD box display of 15 time-links takes the final clock difference of the IGS week SNX file as a reference; blue “+” stands for outliers.
Figure 8. The STD box display of 15 time-links takes the final clock difference of the IGS week SNX file as a reference; blue “+” stands for outliers.
Remotesensing 15 00092 g008
Figure 9. Frequency stability statistics of 6 time-links (PTBB-BOR1, PTBB-GMSD, PTBB-NRC1, PTBB-OHI3, PTBB-TID1, PTBB-TIDV).
Figure 9. Frequency stability statistics of 6 time-links (PTBB-BOR1, PTBB-GMSD, PTBB-NRC1, PTBB-OHI3, PTBB-TID1, PTBB-TIDV).
Remotesensing 15 00092 g009
Figure 10. Frequency stability statistics of TIDV-TID1 time link with common clock.
Figure 10. Frequency stability statistics of TIDV-TID1 time link with common clock.
Remotesensing 15 00092 g010
Figure 11. Muti-GNSS PPPAR time links STD concerning the clock difference of IGS weekly solution stations; yellow ”+” stands for outliers.
Figure 11. Muti-GNSS PPPAR time links STD concerning the clock difference of IGS weekly solution stations; yellow ”+” stands for outliers.
Remotesensing 15 00092 g011
Figure 12. Frequency Stability statistics of Muti-GNSS PPPAR time links (PTBB-VILL, PTBB-USN7, PTBB-GMSD, PTBB-KOUR).
Figure 12. Frequency Stability statistics of Muti-GNSS PPPAR time links (PTBB-VILL, PTBB-USN7, PTBB-GMSD, PTBB-KOUR).
Remotesensing 15 00092 g012
Table 1. Detailed configuration information table of 16 selected stations (from igs.org).
Table 1. Detailed configuration information table of 16 selected stations (from igs.org).
SiteExternal SourceReceiverAntenna
1BOR1EXTERNAL H-MASERTRIMBLE NETR9TRM59800.00
2GMSDEXTERNAL CESIUMTRIMBLE NETR9TRM59800.00
3JFNGINTERNALTRIMBLE ALLOYTRM59800.00
4KOUREXTERNAL H-MASERSEPT POLARX5SEPCHOKE_B3E6
5MGUEEXTERNAL H-MASERSEPT
POLARX5TR
LEIAR25.R4
6MIZUINTERNALSEPT ASTERX4SEPCHOKE_B3E47
7NRC1EXTERNAL H-MASERJAVAD TRE_G3TH DELTASEPCHOKE_B3E37
8OHI3EXTERNAL H-MASERLEICA GR50SEPCHOKE_B3E49
9REDUEXTERNAL CESIUMSEPT POLARX5SEPCHOKE_B3E102
10STR1EXTERNAL CESIUMSEPT POLARX5SEPCHOKE_B3E151
11VILLEXTERNAL CESIUMSEPT POLARX5SEPCHOKE_B3E6
12YEL2EXTERNAL VCH-1008 MASERSEPT
POLARX5TR
LEIAR25.R4
13USN7EXTERNAL H-MASERSEPT
POLARX5TR
TPSCR.G5
14PTBBEXTERNAL ACTIVE H-MASERSEPT
POLARX5TR
SEPCHOKE_B3E89
15TIDVEXTERNAL H-MASERSEPT POLARX5SEPCHOKE_B3E169
16TID1EXTERNAL H-MASERSEPT POLARX5SEPCHOKE_B3E169
Table 2. An overview of the PPP AR products.
Table 2. An overview of the PPP AR products.
AgencyStrategyConstellation
CNESIRCGRE
WUMFCBGREB2J
CODEOSBsGREB2
PRIDEOSBsGREB2B3
Table 3. PPPAR Detailed Policy.
Table 3. PPPAR Detailed Policy.
ProjectStrategy
Stations16 globally distributed MGEX stations
PeriodJanuary 2020, DOY 001
ObservationsGPS: L1, L2 and Galileo: E1, E5a BDS: B1I, B3I
Interval30 s
Weightingsin elev2, cut-off angle 7°
Phase wind-upPhase wind-up Corrected
Tropospheric delay (ZHD)Global Pressure and Temperature (GPT) [27] model using the formulas of Saastamoinen
Tropospheric delay (ZWD)Estimated as a continuous piecewise linear function (2 h parameter spacing), GMF mapping function [28]
Tidal displacementsIERS conventions 2010 [29]
Relativistic effectCorrected
Sagnac Corrected
Satellite antenna PCOs and PCVsFixed to the values from igs14.atx
Table 4. 6 ways of 3 fixed strategies measures taken in this paper.
Table 4. 6 ways of 3 fixed strategies measures taken in this paper.
StrategyExtended Strategy
FCBFCB-com
FCB-gbm
FCB-grm
IRCIRC-grm
OSBOSB-com
OSB-whu
Table 5. 15 detailed fixed strategies for multi-system.
Table 5. 15 detailed fixed strategies for multi-system.
StrategyExtended StrategyExtended Extended Strategy
FCBFCB-comFCB_com_fix_G
FCB_com_fix_GE
FCB-gbmFCB_gbm_fix_G
FCB_gbm_fix_GE
FCB-grgFCB_grm_fix_G
FCB_grm_fix_GE
IRCIRC-grmIRC_grm_fix_G
IRC_grm_fix_GE
IRC_grm_fix_GEC
OSBOSB-comOSB_com_fix_G
OSB_com_fix_GE
OSB_com_fix_GEC
OSB-whuOSB_whu_fix_G
OSB_whu_fix_GE
OSB_whu_fix_GEC
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ouyang, M.; Zhu, X.; Yi, R.; Lyu, D. Using Various Analysis Center Products to Assess the Time-Frequency Transfer Performance of GPS/Galileo/BDS PPPAR Methods. Remote Sens. 2023, 15, 92. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15010092

AMA Style

Ouyang M, Zhu X, Yi R, Lyu D. Using Various Analysis Center Products to Assess the Time-Frequency Transfer Performance of GPS/Galileo/BDS PPPAR Methods. Remote Sensing. 2023; 15(1):92. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15010092

Chicago/Turabian Style

Ouyang, Mingjun, Xiangwei Zhu, Ruite Yi, and Daqian Lyu. 2023. "Using Various Analysis Center Products to Assess the Time-Frequency Transfer Performance of GPS/Galileo/BDS PPPAR Methods" Remote Sensing 15, no. 1: 92. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15010092

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