Next Article in Journal
Variational Generative Adversarial Network with Crossed Spatial and Spectral Interactions for Hyperspectral Image Classification
Next Article in Special Issue
LEO-Constellation-Augmented BDS Precise Orbit Determination Considering Spaceborne Observational Errors
Previous Article in Journal
Measurement of Forest Inventory Parameters with Apple iPad Pro and Integrated LiDAR Technology
Previous Article in Special Issue
Estimation and Analysis of the Observable-Specific Code Biases Estimated Using Multi-GNSS Observations and Global Ionospheric Maps
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Comparison of Multi-GNSS Time and Frequency Transfer Performance Using Overlap-Frequency Observations

1
National Time Service Center, Chinese Academy of Sciences, Xi’an 710600, China
2
Key Laboratory of Time and Frequency Primary Standards, Chinese Academy of Sciences, Xi’an 710600, China
3
State Key Laboratory of Geo-Information Engineering, Xi’an Research Institute of Surveying and Mapping, Xi’an 710054, China
4
Key Laboratory of Precision Navigation and Timing Technology, Chinese Academy of Sciences, Xi’an 710600, China
5
University of Chinese Academy of Sciences, Yu Quan Road, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(16), 3130; https://doi.org/10.3390/rs13163130
Submission received: 5 July 2021 / Revised: 31 July 2021 / Accepted: 3 August 2021 / Published: 7 August 2021
(This article belongs to the Special Issue Beidou/GNSS Precise Positioning and Atmospheric Modeling)

Abstract

:
The modernized GPS, Galileo, and BeiDou global navigation satellite system (BDS3) offers new potential for time transfer using overlap-frequency (L1/E1/B1, L5/E5a/B2a) observations. To assess the performance of time and frequency transfer with overlap-frequency observations for GPS, Galileo, and BDS3, the mathematical models of single- and dual-frequency using the carrier-phase (CP) technique are discussed and presented. For the single-frequency CP model, the three-day average RMS values of the L5/E5a/B2a clock difference series were 0.218 ns for Galileo and 0.263 ns for BDS3, of which the improvements were 36.2% for Galileo and 43.9% for BDS3 when compared with the L1/E1/B1 solution at BRUX–PTBB. For the hydrogen–cesium time link BRUX–KIRU, the RMS values of the L5/E5a/B2a solution were 0.490 ns for Galileo and 0.608 ns for BDS3, improving Galileo by 6.4% and BDS3 by 12.5% when compared with the L1/E1/B1 solution. For the dual-frequency CP model, the average stability values of the L5/E5a/B2a solution at the BRUX–PTBB time link were 3.54 × 10−12 for GPS, 2.20 × 10−12 for Galileo, and 2.69 × 10−12 for BDS3, of which the improvements were 21.0%, 45.1%, and 52.3%, respectively, when compared with the L1/E1/B1 solution. For the BRUX–KIRU time link, the improvements were 4.2%, 30.5%, and 36.1%, respectively.

1. Introduction

Over the past few decades, the global navigation satellite system (GNSS) has been a useful and efficient spatial technique in the area of time transfer. Since the early 1980s, the global positioning system (GPS) was employed to compare remote clocks using the common-view (CV) approach [1,2], which provided precision on the order of a few nanoseconds [3]. Thereafter, the all-in-view (AV) approach was proposed to overcome the limitations of the common-view satellite condition and pseudorange measurements in the CV approach [4]. Considering that GPS carrier phase observations are two orders of magnitude more precise than that of the pseudorange measurement, a GPS carrier-phase (CP) approach was proposed for precise time and frequency transfer [5,6,7], which exhibited better performance when compared to the CV and AV approaches. Therefore, the GPS CP approach became a primary means in time laboratories worldwide.
In recent years, with the rapid development of other emerging GNSSs, such as the European Galileo and Chinese BeiDou System (BDS), research on the use of these GNSSs for precise time and frequency transfer has aroused considerable interest. Martínez-Belda et al. investigated new analysis procedures based on the E5 code plus-carrier (CPC) combination for time transfer [8,9]. Huang and Defraigne extended the standard Common GNSS Generic Time Transfer Standard (CGGTTS) software tool to the BDS regional system (BDS2) [10]. Liang et al. participated in the experiment of BDS2 time transfer in coordinated universal time (UTC), of which the stability and accuracy were evaluated by using common clock difference and multiple inter-continental baselines [11]. Tu et al. proposed a model of precise time transfer based on the BDS triple-frequency uncombined observations [12,13]. Defraigne and Verhasselt rebuilt the CGGTTS-V2E software to support the GPS, GLONASS, BeiDou, and Galileo system [14]. Zhang et al. demonstrated that the combination of GPS, BDS, and Galileo for time and frequency transfer outperforms the GPS-only, BDS-only, and Galileo-only solutions [15]. Meanwhile, the effect of the biases between the BDS global system (BDS3) and BDS2, GPS, Galileo, and GLONASS observations, on precise time and frequency transfer, were also studied [16,17,18]. Moreover, the performance of time transfer based on new-generation BDS3 triple-frequency observation was also investigated [19].
In fact, the new-generation BDS3 has a full constellation of 35 satellites, consisting of five GEO satellites, 27 MEO satellites, and three IGSO satellites. It retains the previous B1I and B3I signals of the BDS-2 satellites, but also includes the new signals B1C (1575.42 MHz) and B2a (1176.45 MHz). Considering that the E1 (1575.42 MHz) and E5a (1176.45 MHz) are equipped for the Galileo system, and the signals of L1 (1575.42 MHz) and L5 (1176.45 MHz) for the Block IIF satellites of GPS, which have formed two overlapping frequencies to enhance compatibility and interoperability between BDS3, GPS, and Galileo system. Although previous studies have focused on multi-GNSS time and frequency transfer, the performance of these overlapping frequency measurements for multi-GNSS time transfer is still unclear.
In this study, single-frequency (SF) and dual-frequency (DF) CP mathematical models were developed to assess the performance of GNSS overlap-frequency measurement time and frequency transfer. We begin with a brief description of the mathematical models of the SF and DF CP time and frequency transfer with overlapping frequency observations. Then, the experimental design and data processing strategy for assessing the performance of the models are discussed. The results and discussion are presented in this section. Finally, a summary and conclusions are drawn.

2. Principle of Multi-GNSS Time Transfer Using Overlap-Frequency Observations

2.1. Model of the SF Time and Frequency Transfer

For the SF GNSS measurement, the observation equations of phase φ and code P may be written as follows [20]:
{ P j , k s = ρ k s + T k s + μ j l k s + d t r , k d t k s + D j d j s + ε P j , k s φ j , k s = ρ k s + T k s μ j l k s + d t r , k d t k s + B j b j s + N j s + ε φ j , k s
where the subscripts j and k denote the frequency and epoch, respectively, and the superscript s denotes the satellite. ρ k s denotes the geometric distance between the satellite and receiver, T k s is the tropospheric delay, and l k s denotes the ionospheric delay at the first frequency with μ j = f 1 / f j , d t r , k   and d t k s are the clock offsets for the receiver and satellite, respectively. B j and D j are the hardware delays of phase and code for the receiver, while b j s and d j s are the hardware delays of phase and code for the satellite. N j s denotes the ambiguity of frequency j . ε φ j , k s and ε P j , k s denote the observation noises of phase and code, respectively. It should be noted that the GNSS satellite and receiver phase center offset and variation, phase wind-up, ocean tides, and relativistic delay must be corrected in the Equation (1).
Generally, the satellite clock offset is corrected using clock products, which contain the satellite uncalibrated delay of ionosphere-free pseudorange observations.
d t s I F = d t s + f 1 2 f 1 2 f 2 2 d 1 s f 2 2 f 1 2 f 2 2 d 2 s = d t s + d 1 s β · D C B 12 s
where β = f 2 2 f 1 2 f 2 2 , D C B 12 s = d 2 s d 1 s is the satellite differential code bias (DCB) between different frequencies, which can be obtained from the International GNSS Service(IGS) or other GNSS analysis centers. Combining the two equations yields Equation (3):
{ P j , k s = ρ k s + T k s + μ j · l k s + d t r , k d t s I F + D j D C B 1 j s β · D C B 12 s + ε P j , k s φ j , k s = ρ k s + T k s μ j · l k s + d t r , k d t s I F + B j b j s + N j s + ε φ j , k s
The major challenge for GNSS time transfer with SF observation is the mitigation of the ionosphere delay, which cannot be precisely removed owing to the single-frequency data. Considering that ionospheric delay has opposite characteristics on code and CP observations, the GRoup And PHase Ionospheric Correction (GRAPHIC) observation combines the code and CP observations in Equation (1) to remove the ionospheric time delay, the ionosphere-free time and frequency model with SF observation can be written as:
G j , k s = φ j , k s + P j , k s 2 = ρ k s + T k s d t s I F + d t r , k + 1 2 D j e j s + N ¯ j s + 1 2 ( ε φ j , k s + ε P j , k s )
where G j , k s is the GRAPHIC observation; for convenience, we defined e j s = D C B 1 j s + β · D C B 12 s and N ¯ j s = 1 2 ( N j s + B j b j s ) . There are many unknown parameters in Equation (4), resulting in a rank deficiency when only the SF GRAPHIC observation is used. This means that one cannot estimate all these parameters directly using Equation (4) [21,22].
Generally, the observation of phase φ or code P in Equation (3) can also be combined to determine the unknown parameters. Then, the single-epoch, SF observation equation can be written as:
{ G j , k = A k x k + t k τ d t s I F + d t r , k + N ¯ j + 1 2 D j e j + ε G j , k P j , k = A k x k + t k τ d t s I F + d t r , k + I j , k + D j e j + ε P j , k
where G j , k = [ G j , k 1 , , G j , k m ] T and P j , k = [ P j , k 1 , , P j , k m ] T are the vectors of GRAPHIC and code observations at epoch k , respectively. x is the station coordinate vector, and the A k is corresponding design matrix,   τ is the modeling residual of tropospheric zenith delay, t k = [ t j , k 1 , , t j , k m ] T is corresponding mapping coefficients, e j = [ e j 1 , , e j m ] T is the vector of the DCB linear combination. ε G j , k and ε P j , k are the vectors of GRAPHIC and code observation noise, respectively. It should be noted that the hardware delay of code D j for the receiver is easily absorbed in the receiver clock offset d t r , k . It should be noted that I j , k = [ I j , k 1 , , I j , k m ] T is the vector of satellite ionospheric delay, which can be corrected using a global ionospheric map (GIM). The stochastic model is formulated as Q = Q r Q k , where Q r = d i a g ( σ G 2 , σ P 2 ) with σ G , and σ P is the precision of GRAPHIC and code observations. Q k is the cofactor matrix determined using the satellite-elevation-dependent model.
{ d t ¯ r , k = d t r , k + 1 2 D j N ¯ j s = 1 2 ( N j s + B j b j s )

2.2. Model of the DF Time and Frequency Transfer

For the pseudorange and carrier phase observations, the traditional CP technique using the ionosphere-free (IF) combination to eliminate the first-order ionospheric delay, which can be formulated as follows [23,24,25]:
{ P I F , k s = a m · P m , k s + a n · P n , k s = ρ I F , k s + T k s + d t r , k d t s I F + ε P I F , k s φ I F , k s = a m · φ m , k s + a n · φ n , k s = ρ I F , k s + T k s + d t r , k d t s I F + N I F s + ε φ I F , k s
where the indices m and n represent the frequency indices. a m and a n are the coefficients of the IF DF combination. For the GPS and BDS, considering that the satellite clock products are determined by using L1/L2 and B1I/B3I observations, respectively, the overlapping frequency observation E1/E5a is not used, such as Galileo, because it will absorb the satellite DCB into clock products, which will require satellite DCBs among the code to extra correct and remain compatible with the current precise clock product:
DCB cp s ( GPS ) = β L 1 , L 2 · DCB L 1 , L 2 s β L 1 , L 5 · DCB L 1 , L 5 s
DCB cp s ( BDS ) = β B 1 , B 3 · DCB B 1 , B 3 s β B 1 , B 2 a · DCB B 1 , B 2 a s
Therefore, the single-epoch, DF observation equation can be written as
{ P I F , k = A k x k + t k τ d t s I F + d t r , k + ε P I F , k φ I F , k = A k x k + t k τ d t s I F + d t r , k + N I F , k + ε P I F , k
where P I F , k = [ P I F , k 1 , , P I F , k m ] T and φ I F , k = [ φ I F , k 1 , , φ I F , k m ] T are the vectors of single-epoch IF code and carrier phase observations, respectively.
For the operation of precise GNSS time and frequency transfer, this is usually carried out between two remote time and frequency references (namely, A and B). Both the SF model and DF model for time transfer, the unknown parameter of receiver clock offset d t r , k is interested, which denotes the clock difference between the individual GNSS timescale (GNSST) and the time and frequency reference. When the hardware delays caused by the receiver, antenna, and cables are calibrated [26,27], then the formula is expressed as follows:
T ( A ) T ( B ) = ( d t r , k ( A ) GNSST ) ( d t r , k ( B ) GNSST ) = d t r , k ( A ) d t r , k ( B )
It can be noted that the GNSS time scales plays the intermediary role in time and frequency transfer, which effectively overcome the limitation of geometric distance between two time and frequency references.

3. Time and Frequency Transfer Experiment

To assess the performance of time and frequency transfer using GPS, BDS3, and Galileo overlapping-frequency observations, three multi-GNSS stations equipped with an external atomic clock were selected, of which the attribute information are summarizes in Table 1. It should be noted that the two different UTC(k) were connected to the BRUX and PTBB stations, which were maintained using a hydrogen maser. We selected BRUX as the reference station for this experiment. Hence, two time links (BRUX–PTBB and BRUX–KIRU) were formed to analyze the performance of time and frequency transfer.
During the data processing, the GNSS time and frequency transfer software named precise time transfer solution (PTTSol) was developed and employed to assess the overlap-frequency time transfer performance. Considering that most of the error models for BDS3 can be directly utilized for the other systems (GPS and Galileo), owing to the similar operation principles for these GNSSs. All observations had a sampling rate of 30 s, and the satellite elevation mask angle was set to 7°. The multi-GNSS mixed satellite orbit and clock products provided by the GNSS Research Center of Wuhan University were used at sampling intervals of 5 min and 30 s, respectively [28]. The ionospheric effect was eliminated and weakened by using the carrier phase, pseudorange combination, and CODE’s GIM for the SF model. For the DF model, the ionospheric effect was eliminated using the IF combination, which can eliminate the first-order delay [29,30]. The parameter of the satellite antenna phase center correction was temporarily unavailable for BDS3; therefore, it was not corrected. The daily multi-GNSS DCBs were derived from the Chinese Academy of Sciences [31]. Phase ambiguities were estimated as float values. The initial standard deviation values for the carrier phase and pseudorange observations were set to 0.003 and 0.3 m, respectively [32]. The tropospheric delay is typically corrected for its hydrostatic components using the Saastamoinen model, while the residual non-hydrostatic wet components are modeled as a random-walk process with noise (10−9 m2/s) [33]. The parameter of interest of the receiver clock offset was modeled as a white noise stochastic process, implying that the parameter was estimated epoch-by-epoch. The details of the observation models and data processing strategies are provided in Table 2.
Notably, the two external atomic clocks individually equipped for one time link are still free-running during the operation of time and frequency transfer, it is impossible to evaluate the performance by repeated measurement. Generally, some smooth approaches with a certain noise are used to assess the performance of time transfer [34]. Therefore, the RMS values of the smooth residuals for the time transfer results were employed to analyze the noise levels of different solutions. Here, the Kalman filter with a process noise of 0.001 ns in the forward mode is employed to determine the RMS values. Meanwhile, the Allan variance is a measure of the fractional frequency fluctuation and offers the advantage of convergence, which is the most commonly used in the area of time and frequency to measure the frequency stability. The Allan variance can be calculated as:
σ y 2 ( τ ) = 1 2 ( N 2 ) τ 2 i = 1 N 2 [ x i + 2 2 x i + 1 + x i ] 2
where x i is the i th data, and τ is the measurement interval. In fact, the value is usually expressed as the form of square root σ y ( τ ) , called the Allan deviation (ADEV). The Allan variance is the same as the ordinary variance for white frequency modulation noise, but has the advantage, for more divergent noise types, such as flicker noise, of converging to a value that is independent on the number of samples.

4. Results and Discussion

In this section, based on the results obtained from different schemes, the performance of overlapping frequency observations is analyzed based on the proposed SF-CP and DF-CP models with GPS, Galileo, and BDS3 constellations. Then, some characteristics of the L1/E1/B1 and L5/E5a/B2a observations are discussed. Finally, the performance in noise level and frequency stability of two time links, including the SF and DF for overlapping frequency observations, is discussed.

4.1. SF Time and Frequency Transfer

Figure 1 and Figure 2 present the SF time transfer results for GPS, Galileo, and BDS3 overlapping frequency observations at time link BRUX–PTBB and BRUX–KIRU, respectively. One can note that the variations in the time transfer values made a good match among GPS, Galileo, and BDS3 at each overlapping frequency solution. Moreover, the clear biases still existed among the GPS, Galileo, and BDS3 for each overlapping frequency solution, which is mainly caused by receiver clock offset being absorbed by different receiver DCBs as varying GNSS undergoing different time delays. Additionally, one can also note that the result of time transfer at L5/E5a/B2a exhibits better variation characteristics than L1/E1/B1 at both the time links of BRUX–PTBB and BRUX–KIRU.
To further compare the performance of time and frequency transfer with overlapping frequency observations for GPS, Galileo, and BDS3 systems, the indicators of noise level and frequency stability of time links are employed. The RMS of smoothed residuals for the clock difference series derived by the SF models at time link BRUX–PTBB are shown in the left panel of Figure 3. For comparison, the BRUX–KIRU solutions are also presented in the right panel. One can clearly note that the time link BRUX–PTBB shows better performance than that of the BRUX–KIRU time link, the main reason is that the former one equipped with the hydrogen–hydrogen clocks, whereas the latter one is equipped with the hydrogen–cesium clocks. As shown, the L5/E5a/B2a solutions show better performance than those of the L1/E1/B1 solutions at two time links, particularly in Galileo and BDS3 solutions. The three-day averaged RMS values of the L5/E5a/B2a solution were 0.218 ns for Galileo and 0.263 ns for BDS3, improving Galileo by 36.2%, and BDS3 by 43.9% when compared with the L1/E1/B1 solution at BRUX–PTBB. For the hydrogen–cesium time link BRUX–KIRU, the corresponding RMS values of the L5/E5a/B2a solution are 0.490 ns for Galileo and 0.608 ns for BDS3, improving Galileo by 6.4% and BDS3 by 12.5% when compared with the L1/E1/B1 solution at BRUX–KIRU. From Table 3, it can be noted that because the limited GPS satellite is able to broadcast the L5 signals, the advantage of L5/E5a/B2a solutions is not obvious, particularly in the BRUX–KIRU time link.
Figure 4 and Figure 5 present the ADEVs for SF time transfer results with GPS, Galileo, and BDS3 overlapping frequency observations at time link BRUX–PTBB and BRUX–KIRU. One can clearly see that the frequency stabilities of L1/E1/B1 solutions were also generally similar for the GPS, Galileo, and BDS3 both BRUX–PTBB and BRUX–KIRU time links, regardless of the averaging time interval. For the BRUX–PTBB time link, the average stability values of the L5/E5a/B2a solution at all time intervals are 3.54 × 10−12 for GPS, 2.20 × 10−12 for Galileo, and 2.69 × 10−12 for BDS3, of which the improvements are 21.0%, 45.1%, and 52.3%, respectively, when compared with the L1/E1/B1 solution. With respect to the BRUX–KIRU time link, the average stability values of the L5/E5a/B2a solution at all time intervals are 4.94 × 10−12 for GPS, 3.06 × 10−12 for Galileo, and 3.82 × 10−12 for BDS3, of which the improvements are 4.2%, 30.5%, and 36.1%, respectively, when compared with the L1/E1/B1 solution.

4.2. DF Time and Frequency Transfer

Figure 6 shows the DF time transfer results for GPS, Galileo, and BDS3 overlapping frequency (L1/L5, E1/E5a, B1/B2a) observations at time link BRUX–PTBB (left) and BRUX–KIRU (right), respectively. One can note that the variation of the clock difference series for GPS, Galileo, and BDS3 agreed well both at the BRUX–PTBB (Left) and BRUX–KIRU (Right), although obvious biases exist among the different series, which mainly cause by the different hardware delays for different satellite systems and can be calibrated in the time transfer campaign [35]. Meanwhile, the different estimated coordinate values of station for different satellite systems partly lead to these biases [36,37]. Therefore, we focused on the variation characteristics of the clock difference series. Table 4 presents the DF RMS of the smoothed residuals for the clock difference series. It can be noted that the RMS values are in good agreement with each GNSS at the time link of BRUX–PTBB and BRUX–KIRU, respectively. The averaged RMS improved by nearly an order of magnitude for the BRUX–PTBB time link when comparing the SF L5/E5a/B2a RMS mean values. The average improvements of RMS are 93.7%, 93.4%, and 94.0% for GPS, Galileo, and BDS3, at the BRUX–PTBB time link, and 39.7%, 25.5%, and 43.0% for GPS, Galileo, and BDS3, respectively, at the BRUX–KIRU time link.t.
Figure 7 presents the ADEVs for DF time transfer results with GPS, Galileo, and BDS3 overlapping frequency observations at time link BRUX–PTBB and BRUX–KIRU. It can be clearly seen that the variation trends are very similar for different GNSSs at the two time links. For the BRUX–PTBB time link, the average stability values of the DF solution at all time intervals were 7.80 × 10−14 for GPS, 6.95 × 10−14 for Galileo, and 7.43 × 10−14 for BDS3, of which the improvements were 97.8%, 96.8%, and 97.2%, respectively, when compared with the SF L5/E5a/B2a solution. With respect to the BRUX–KIRU time link, the average stability values of the DF solution at all time intervals were 3.82 × 10−13 for GPS, 3.21 × 10−13 for Galileo, and 3.25 × 10−13 for BDS3, with improvements of 92.3%, 89.5%, and 91.5%, respectively, when compared with the SF L5/E5a/B2a solution.

5. Conclusions

The multi-GNSS overlap-frequency observation will contribute to the compatible interoperability among different GNSSs. To further exploit the potential of those GNSS overlap-frequency observation for time and frequency transfer, we studied a SF model and DF model of GNSS precise time and frequency transfer by using the CP technique. The corresponding mathematical models are analyses and presented. Moreover, numerical analyses were conducted to assess the time and frequency transfer performance of two time transfer links.
For the SF CP time transfer model, the results show that the variations in the clock difference series from the three GNSS solutions agree well, at both the L1/E1/B1 and L5/E5a/B2a overlapping frequency. With respect to the indicators of noise level and frequency stability of the time link, the L5/E5a/B2a solutions show better performance than the L1/E1/B1 solutions at two time links. In terms of noise level, the improvements in L5/E5a/B2a were 36.2% for Galileo, and 43.9% for BDS3 when compared with the L1/E1/B1 solution at hydrogen–hydrogen BRUX–PTBB, as well as the 6.4% for Galileo, and 12.5% for BDS3 at the hydrogen–cesium time link BRUX–KIRU. For the frequency stability, the improvement in the average stability value of L5/E5a/B2a solution was 21.0% for GPS, 45.1% for Galileo, and 52.3% for BDS3 at the BRUX–PTBB time link when compared with L1/E1/B1 solutions, while the improvements were 4.2% for GPS, 30.5% for Galileo, and 36.1% for BDS3 at the BRUX–KIRU time link.
For the DF CP time transfer model, the variation of the clock difference series for GPS, Galileo, and BDS3 overlapping frequency observations agreed well at two time links. The DF result shows a significant improvement over the SF time and frequency transfer results in both the noise level and frequency stability. When compared with the L5/E5a/B2a SF solutions, the average improvements in noise level were 93.7%, 93.4%, and 94.0% for GPS, Galileo, and BDS3 at the BRUX–PTBB time link, as well as the 39.7%, 25.5%, and 43.0% for GPS, Galileo, and BDS3 at the BRUX–KIRU time link, respectively. For the frequency stability, the improvements were 97.8%, 96.8%, and 97.2% for GPS, Galileo, and BDS3, respectively at the BRUX–PTBB time link. For BRUX–KIRU, the improvements were 92.3%, 89.5%, and 91.5% for GPS, Galileo, and BDS3, respectively.
In conclusion, when using the overlapping frequency observation for multi-GNSS precise time and frequency transfer, the results exhibit equivalent performance for different GNSSs based on both SF and DF models. Meanwhile, the L5/E5a/B2a frequency observations have certain advantages in terms of time and frequency transfer. Then, the DF result shows significant improvement over the SF, of which the range can reach nearly an order of magnitude. With continuous development of the GPS, Galileo, and BeiDou satellite systems, more GNSS satellites will broadcast overlapping frequency signals, which will be helpful for the compatibility of multi-system time and frequency transfer.

Author Contributions

P.Z. and R.T. conceived and designed the experiments; P.Z. performed the experiments, analyzed the data, and wrote the paper. X.L., Y.G., J.H. (Ju Hong) and J.H. (Junqiang Han) contributed to the discussions and revisions. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the National Natural Science Foundation of China (grant no. 11903040, 41974032) and Chinese Academy of Sciences (CAS) programs of “High Level Talents”, “The Frontier Science Research Project” (grant no. QYZDB-SSW-DQC028), “Western young scholars” (grant no. XAB2019B21), China Postdoctoral Science Foundation (grant no. 2020M683763).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets analyzed in this study are managed by the MGEX and Wuhan University and are available upon request from the corresponding author.

Acknowledgments

We greatly appreciate MGEX and Wuhan University for providing multi-GNSS observation data and products.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Allan, D.W.; Weiss, M. Accurate time and frequency transfer during common-view of a GPS satellite. In Proceedings of the 1980 IEEE Frequency Control Symposium, Philadelphia, PA, USA, 28–30 May 1980; pp. 334–356. [Google Scholar]
  2. Allan, D.W.; Thomas, C. Technical directives for standardization of GPS time receiver software. Metrologia 1994, 31, 69–79. [Google Scholar] [CrossRef]
  3. Petit, G.; Jiang, Z. GPS All in View time transfer for TAI computation. Metrologia 2008, 45, 35–45. [Google Scholar] [CrossRef]
  4. Jiang, Z.; Lewandowski, W. Use of GLONASS for UTC time transfer. Metrologia 2012, 49, 57–61. [Google Scholar] [CrossRef]
  5. Orgiazzi, D.; Tavella, P.; Lahaye, F. Experimental assessment of the time transfer capability of precise point positioning (PPP). In Proceedings of the IEEE International Frequency Control Symposium and Exposition (FCS ’05), Vancouver, BC, Canada, 29–31 August 2005; pp. 337–345. [Google Scholar]
  6. Defraigne, P.; Bruyninx, C.; Guyennon, N. PPP and phase only GPS frequency transfer. In Proceedings of the IEEE International Frequency Control Symposium Jointly with the 21st European Frequency and Time Forum (EFTF ’07), Geneva, Switzerland, 29 May–1 June 2007; pp. 904–908. [Google Scholar]
  7. Zhang, P.; Tu, R.; Gao, Y.; Liu, N.; Zhang, R. Improving Galileo’s carrier-phase time transfer based on prior constraint information. J. Navig. 2019, 72, 121–139. [Google Scholar] [CrossRef]
  8. Martínez-Belda, M.C.; Defraigne, P.; Baire, Q.; Aerts, W. Single-frequency time and frequency transfer with Galileo E5. In Proceedings of the Joint Conference of the IEEE International Frequency Control and the European Frequency and Time Forum (FCS) Proceedings, San Francisco, CA, USA, 2–5 May 2011; pp. 1–6. [Google Scholar]
  9. Martinez-Belda, M.C.; Defraigne, P.; Bruyninx, C. On the potential of Galileo E5 for time transfer. IEEE Trans. Ultrason Ferroelectr. Freq. Control 2013, 60, 121–131. [Google Scholar] [CrossRef]
  10. Huang, W.; Defraigne, P. BeiDou Time Transfer with the Standard CGGTTS. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2016, 63, 1005–1012. [Google Scholar] [CrossRef]
  11. Liang, K.; Arias, F.; Petit, G.; Jiang, Z.; Tisserand, L.; Wang, Y.; Yang, Z.; Zhang, A. Evaluation of BeiDou time transfer over multiple inter-continental baselines towards UTC contribution. Metrologia 2018, 55, 513–525. [Google Scholar] [CrossRef]
  12. Tu, R.; Zhang, P.; Zhang, R.; Liu, J.; Lu, X. Modeling and performance analysis of precise time transfer based on BDS triple-frequency un-combined observations. J. Geod. 2018, 93, 837–847. [Google Scholar] [CrossRef]
  13. Tu, R.; Zhang, P.; Zhang, R.; Liu, J.; Lu, X. Modeling and Assessment of Precise Time Transfer by Using BeiDou Navigation Satellite System Triple-Frequency Signals. Sensors 2018, 18, 1017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Defraigne, P.; Verhasselt, K. Multi-GNSS time transfer with CGGTTS-V2E. In Proceedings of the European Frequency and Time Forum (EFTF), Turin, Italy, 10–12 April 2018; pp. 270–275. [Google Scholar]
  15. Zhang, P.; Tu, R.; Zhang, R.; Gao, Y.; Cai, H. Combining GPS, BeiDou, and Galileo Satellite Systems for Time and Frequency Transfer Based on Carrier Phase Observations. Remote Sens. 2018, 10, 324. [Google Scholar] [CrossRef] [Green Version]
  16. Zhang, P.; Tu, R.; Gao, Y.; Zhang, R.; Han, J. Performance of Galileo precise time and frequency transfer models using quad-frequency carrier phase observations. GPS Solut. 2020, 24, 1–18. [Google Scholar] [CrossRef]
  17. Nicolini, L.; Caporali, A. Investigation on Reference Frames and Time Systems in Multi-GNSS. Remote Sens. 2018, 10, 80. [Google Scholar] [CrossRef] [Green Version]
  18. Montenbruck, O.; Steigenberger, P.; Hauschild, A. Broadcast versus precise ephemerides: A multi-GNSS perspective. GPS Solut. 2015, 19, 321–333. [Google Scholar] [CrossRef]
  19. Su, K.; Jin, S. Triple-frequency carrier phase precise time and frequency transfer models for BDS-3. GPS Solut. 2019, 23. [Google Scholar] [CrossRef]
  20. Afifi, A.; El-Rabbany, A. Performance Analysis of Several GPS/Galileo Precise Point Positioning Models. Sensors 2015, 15, 14701–14726. [Google Scholar] [CrossRef] [Green Version]
  21. Li, B.; Zang, N.; Ge, H.; Shen, Y. Single-frequency PPP models: Analytical and numerical comparison. J. Geod. 2019, 93, 2499–2514. [Google Scholar] [CrossRef]
  22. Zhang, B.; Teunissen, P.; Yuan, Y. Joint estimation of vertical total electron content (VTEC) and satellite differential code biases (SDCBs) using low-cost receivers. J. Geod. 2018, 92, 401–413. [Google Scholar] [CrossRef] [Green Version]
  23. Zhao, Q.; Wang, G.; Liu, Z.; Hu, Z.; Dai, Z.; Liu, J. Analysis of BeiDou satellite measurements with code multipath and geometry-free ionospheric-free combinations. Sensors 2016, 16, 123. [Google Scholar] [CrossRef] [Green Version]
  24. Najibi, N.; Jin, S. Physical reflectivity and polarization characteristics for snow and ice-covered surfaces interacting with GPS signals. Remote Sens. 2013, 5, 4006–4030. [Google Scholar] [CrossRef] [Green Version]
  25. Zhang, B.; Chen, Y.; Yuan, Y. PPP-RTK based on undifferenced and uncombined observations: Theoretical and practical aspects. J. Geod. 2019, 93, 1011–1024. [Google Scholar] [CrossRef]
  26. Rovera, G.D.; Torre, J.M.; Sherwood, R.; Abgrall, M.; Courde, C.; Laas-Bourez, M.; Uhrich, P. Link calibration against receiver calibration: An assessment of GPS time transfer uncertainties. Metrologia 2014, 51, 476–490. [Google Scholar] [CrossRef]
  27. Zhang, P.; Tu, R.; Gao, Y.; Guang, W.; Zhang, R.; Cai, H. Study of time link calibration based on GPS carrier phase observation. IET Radar Sonar Navig. 2018, 12, 1330–1335. [Google Scholar] [CrossRef]
  28. Xie, X.; Geng, T.; Zhao, Q.; Cai, H.; Zhang, F.; Wang, X.; Meng, Y. Precise orbit determination for BDS-3 satellites using satellite-ground and inter-satellite link observations. GPS Solut. 2019, 23, 1–12. [Google Scholar] [CrossRef]
  29. Wang, J.; Huang, G.; Zhang, Q.; Gao, Y.; Gao, Y.; Luo, Y. GPS/BDS-2/Galileo Precise Point Positioning Ambiguity Resolution Based on the Uncombined Model. Remote Sens. 2020, 12, 1853. [Google Scholar] [CrossRef]
  30. Li, M.; Yuan, Y. Estimation and Analysis of BDS2 and BDS3 Differential Code Biases and Global Ionospheric Maps Using BDS Observations. Remote Sens. 2021, 13, 370. [Google Scholar] [CrossRef]
  31. Wang, N.; Yuan, Y.; Li, Z.; Montenbruck, O.; Tan, B. Determination of differential code biases with multi-GNSS observations. J. Geod. 2016, 90, 209–228. [Google Scholar] [CrossRef]
  32. Lin, P.; Li, X.; Zhang, X.; Li, X.; Lu, C.; Zhao, Q.; Liu, J. Considering Inter-Frequency Clock Bias for BDS Triple-Frequency Precise Point Positioning. Remote Sens. 2017, 9, 734. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, Q.; Jin, S.; Yuan, L.; Hu, Y.; Chen, J.; Guo, J. Estimation and Analysis of BDS-3 Differential Code Biases from MGEX Observations. Remote Sens. 2020, 12, 68. [Google Scholar] [CrossRef] [Green Version]
  34. 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]
  35. Zhang, P.; Tu, R.; Gao, Y.; Zhang, R.; Liu, N. Impact of BeiDou satellite-induced code bias variations on precise time and frequency transfer. Meas. Sci. Technol. 2019, 30, 035007. [Google Scholar] [CrossRef]
  36. Yao, J.; Skakun, I.; Jiang, Z.; Levine, J. A detailed comparison of two continuous GPS carrier-phase time transfer techniques. Metrologia 2015, 52, 666–676. [Google Scholar] [CrossRef]
  37. Defraigne, P.; Baire, Q. Combining GPS and GLONASS for time and frequency transfer. Adv. Space Res. 2011, 47, 265–275. [Google Scholar] [CrossRef]
Figure 1. SF time transfer result with GPS, Galileo, and BDS3 overlapping frequency observations in L1/E1/B1 (left) and L5/E5a/B2a (right) at time link BRUX–PTBB.
Figure 1. SF time transfer result with GPS, Galileo, and BDS3 overlapping frequency observations in L1/E1/B1 (left) and L5/E5a/B2a (right) at time link BRUX–PTBB.
Remotesensing 13 03130 g001
Figure 2. SF time transfer result with GPS, Galileo, and BDS3 overlapping frequency observations in L1/E1/B1 (left) and L5/E5a/B2a (right) at time link BRUX–KIRU.
Figure 2. SF time transfer result with GPS, Galileo, and BDS3 overlapping frequency observations in L1/E1/B1 (left) and L5/E5a/B2a (right) at time link BRUX–KIRU.
Remotesensing 13 03130 g002
Figure 3. RMS of smoothed residuals for the clock difference series derived by the SF models at time link BRUX–PTBB (left) and BRUX–KIRU (right) two time links.
Figure 3. RMS of smoothed residuals for the clock difference series derived by the SF models at time link BRUX–PTBB (left) and BRUX–KIRU (right) two time links.
Remotesensing 13 03130 g003
Figure 4. Comparison of ADEVs for SF time transfer result with GPS, Galileo, and BDS3 overlapping frequency observations at time link BRUX–PTBB.
Figure 4. Comparison of ADEVs for SF time transfer result with GPS, Galileo, and BDS3 overlapping frequency observations at time link BRUX–PTBB.
Remotesensing 13 03130 g004
Figure 5. Comparison of ADEVs for SF time transfer result with GPS, Galileo, and BDS3 overlapping frequency observations at time link BRUX–KIRU.
Figure 5. Comparison of ADEVs for SF time transfer result with GPS, Galileo, and BDS3 overlapping frequency observations at time link BRUX–KIRU.
Remotesensing 13 03130 g005
Figure 6. DF time transfer result with GPS, Galileo, and BDS3 overlapping frequency observations in at time link BRUX–PTBB (Left) and BRUX–KIRU (Right).
Figure 6. DF time transfer result with GPS, Galileo, and BDS3 overlapping frequency observations in at time link BRUX–PTBB (Left) and BRUX–KIRU (Right).
Remotesensing 13 03130 g006
Figure 7. Comparison of ADEVs for DF time transfer result with GPS, Galileo and BDS3 overlapping frequency observations at time link BRUX–PTBB (Left) and BRUX–KIRU (Right).
Figure 7. Comparison of ADEVs for DF time transfer result with GPS, Galileo and BDS3 overlapping frequency observations at time link BRUX–PTBB (Left) and BRUX–KIRU (Right).
Remotesensing 13 03130 g007
Table 1. Information of stations used in the experiment.
Table 1. Information of stations used in the experiment.
SiteGNSS ReceiverAntennaFrequency StandardLocation
BRUXPOLARX5TRJAVRINGANT_DMUTC(ORB)50.8° N, 4.4° E
PTBBPOLARX5TRLEIAR25.R4UTC(PTB)52.2° N, 10.4° E
KIRUPOLARX5SEPCHOKE_B3E6Industrial cesium standard67.8° N, 20.9° E
Table 2. Observation models and data processing strategies.
Table 2. Observation models and data processing strategies.
ItemsModels and Strategies
ObservationsCarrier phase and pseudorange measurements
Signal frequencyGPS: L1/L5; Galileo: E1/E5a; BDS3: B1/B2a
Satellite orbit and clockFinal multi-GNSS products from WUHN analysis center
Elevation cut-off
Sampling rate30 s
TroposphericSaastamoinen model + random-walk process for wet component (10−9 m2/s)
IonosphericDifferent Ionospheric-free linear combination for SF and DF model
EstimatorBatch least-squares estimator
Observation weightElevation-dependent model with cos2 (el) function
Receiver clock offsetEstimated with white noise
Table 3. Average number of available satellite for different GNSS frequencies at three stations.
Table 3. Average number of available satellite for different GNSS frequencies at three stations.
BRUXKIRUPTBB
GPS_L19.210.19.2
GPS_L54.55.04.5
Galileo_E17.68.17.1
Galileo_E5a7.68.17.1
BDS3_B15.25.75.2
BDS3_B2a5.25.75.2
Table 4. DF RMS of smoothed residuals for the clock difference series (ns).
Table 4. DF RMS of smoothed residuals for the clock difference series (ns).
DayBRUX–PTBBBRUX–KIRU
GPSGalileoBDS3GPSGalileoBDS3
2840.0260.0230.0280.4070.4070.406
2850.0230.0200.0200.3450.3440.344
2860.0160.0130.0140.3830.3810.386
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, P.; Tu, R.; Gao, Y.; Hong, J.; Han, J.; Lu, X. Comparison of Multi-GNSS Time and Frequency Transfer Performance Using Overlap-Frequency Observations. Remote Sens. 2021, 13, 3130. https://doi.org/10.3390/rs13163130

AMA Style

Zhang P, Tu R, Gao Y, Hong J, Han J, Lu X. Comparison of Multi-GNSS Time and Frequency Transfer Performance Using Overlap-Frequency Observations. Remote Sensing. 2021; 13(16):3130. https://doi.org/10.3390/rs13163130

Chicago/Turabian Style

Zhang, Pengfei, Rui Tu, Yuping Gao, Ju Hong, Junqiang Han, and Xiaochun Lu. 2021. "Comparison of Multi-GNSS Time and Frequency Transfer Performance Using Overlap-Frequency Observations" Remote Sensing 13, no. 16: 3130. https://doi.org/10.3390/rs13163130

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