[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
A Novel Metamaterial Inspired High-Temperature Microwave Sensor in Harsh Environments
Next Article in Special Issue
An Integrated Dead Reckoning with Cooperative Positioning Solution to Assist GPS NLOS Using Vehicular Communications
Previous Article in Journal
Practical Use of Robot Manipulators as Intelligent Manufacturing Systems
Previous Article in Special Issue
A New Approach to Unwanted-Object Detection in GNSS/LiDAR-Based Navigation
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 18
Issue 9
10.3390/s18092878
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improving the Performance of Multi-GNSS Time and Frequency Transfer Using Robust Helmert Variance Component Estimation

by
Pengfei Zhang
1,2,3,
Rui Tu
1,3,4,*,
Yuping Gao
1,2,
Rui Zhang
1,4 and
Na Liu
1,2,3
1
National Time Service Center, Chinese Academy of Sciences, Shu Yuan Road, Xi’an 710600, China
2
Key Laboratory of Time and Frequency Primary Standards, Chinese Academy of Sciences, Xi’an 710600, China
3
University of Chinese Academy of Sciences, Yu Quan Road, Beijing 100049, China
4
Key Laboratory of Precision Navigation Positioning and Timing Technology, Chinese Academy of Sciences, Xi’an 710600, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(9), 2878; https://doi.org/10.3390/s18092878
Submission received: 11 July 2018 / Revised: 25 August 2018 / Accepted: 27 August 2018 / Published: 31 August 2018
(This article belongs to the Special Issue GNSS and Fusion with Other Sensors)
Browse Figures
Figure 1
<p>Primary implementation procedure using robust Helmert variance component estimation (RVCE) for time and frequency transfer.</p> "> Figure 2
<p>Layouts of multi-GNSS stations, capable of tracking GPS, BDS and Galileo constellations.</p> "> Figure 3
<p>Receiver clock time series for station ONS1 and CEDU with and without the RVCE approach.</p> "> Figure 4
<p>Results of time transfers determined by a carrier-phase (CP) technique with and without the RVCE approach for the four time links. (<b>a</b>) is for ONS1-BRUX, (<b>b</b>) is for NTS1-KAT1, (<b>c</b>) is for NTS1-CEDU and (<b>d</b>) is for BRUX-NTS1.</p> "> Figure 5
<p>Smoothing residuals of the results with and without the RVCE approach for the four time links. (<b>a</b>) is for ONS1-BRUX, (<b>b</b>) is for NTS1-KAT1, (<b>c</b>) is for NTS1-CEDU and (<b>d</b>) is for BRUX-NTS1.</p> "> Figure 6
<p>Comparison of Allan deviation (ADEV) for the four time links determined with and without RVCE. (<b>a</b>) is for ONS1-BRUX, (<b>b</b>) is for NTS1-KAT1, (<b>c</b>) is for NTS1-CEDU and (<b>d</b>) is for BRUX-NTS1.</p> "> Figure 7
<p>Improvement to stability (%) of RVCE solutions over the raw solutions within 10,000 s.</p> ">
Versions Notes

Abstract

:
The combination of multiple Global Navigation Satellite Systems (GNSSs) may improve the performance of time and frequency transfers by increasing the number of available satellites and improving the time dilution of precision. However, the receiver clock estimation is easily affected by the inappropriate weight of multi-GNSSs due to the different characteristics of individual GNSS signals as well as the outliers from observations. Thus, we utilised a robust Helmert variance component estimation (RVCE) approach to determine the appropriate weights of different GNSS observations, and to control for the influence of outliers in these observation in multi-GNSS time and frequency transfer. In order to validate the effectiveness of this approach, four time links were employed. Compared to traditional solutions, the mean improvement of smoothed residuals is 3.43% using the RVCE approach. With respect to the frequency stability of the time links, the RVCE solution outperforms the traditional solution, particularly in the short-term, and the mean improvement is markedly high at 14.89%.

1. Introduction

Since the Global Positioning System (GPS) was first utilized for precise time transfer [1], it has increasingly been applied in time frequency dissemination. From initial common-view (CV) and all-in-view (AV) techniques [2,3], which use only pseudo-range measurements, to the development of carrier-phase (CP) techniques that combine carrier phase observation [4], the accuracy of GPS time transfer has improved significantly. While earlier studies employing CV and AV techniques have recently been re-evaluated, research on GLObal Navigation Satellite System (GLONASS) time transfer has also increasingly been undertaken in those years [5,6]. Recently, other satellite systems have also been well developed, and this has provided a rich resource of Global Navigation Satellite System (GNSS) signals for time and frequency transfer. The European Union’s Galileo satellite system, which was initiated by the European Space Agency (ESA) and the European Commission in 2011 offers one promising E5 signal component with low noise in the centimeter range, and great potential to improve time transfer performance [7,8]. The BeiDou Navigation Satellite System (BDS) of China has been in continuous service since 27 December 2012, and the application of its time transfer has also been investigated [9,10].
Given that combining multi-constellation GNSS can increase the number of available satellites for data retrieval, doing so can improve the redundancy of the receiver clock parameter. While multi-GNSS time transfer has been investigated [11,12,13], such research has primarily focused on combining GPS and GLONASS. As the GLONASS satellites emit signals on individual frequencies via a frequency division multiple access (FDMA) design, the effect of inter-frequency biases (IFBs) [12] and differential code-phase biases (DCPBs) [14] for GLONASS is limited, allowing for the full exploitation of the combination potential of these satellites for precise time and frequency transfer. With the development of precise multi-GNSS satellite orbit and clock products, all satellites within the network use the same time scale for different GNSS products [15,16,17]. This provides the opportunity to combine the three code-division multiple access (CDMA) satellite constellations (i.e., GPS, Galileo and BDS) for precise time and frequency transfer. By increasing the number of available satellites and improving the time dilution of precision, the time transfer solution was determined by combining the three CDMA GNSS constellations, which outperformed single GNSS solutions [18]. However, the quality of orbit and satellite clock products and pseudo-range observation accuracy of Galileo and BDS differ from those of GPS [19], and it is therefore inappropriate to treat them equally as has been done in previous studies [18]. Meanwhile, the pseudo-range observation is an important measurement in time and frequency transfers as it carries temporal information [20,21]. However, outliers are often found in these observations owing to the environments surrounding the stations [22,23], which substantially affect time transfer performance.
In this study, a robust Helmert variance component estimation (RVCE) was introduced when using combined multi-GNSS for precise time and frequency transfer estimates based on a CP technique. This was employed to allocate the adaptive weight for different GNSSs and to suppress the influence of outliers based on their adjustment residuals. Five GNSS stations equipped with external time and frequency references, and which were capable of tracking three GNSS constellations, formed four time links that were processed. The performance of the CP technique with three GNSS observations using the RVCE approach was assessed in terms of both time transfer noise level as well as time and frequency stability [13]. The structure of this manuscript is organized as follows: Section 2 describes the observation equation of the time and frequency transfer based on multi-GNSS. The robust approach based on the Helmert variance component estimation (VCE) is discussed in Section 3, and the implementation procedure is also summarized. Thereafter, Section 4 describes the experiment that was used to assess performance, and the comparison and validation with the traditional solutions are also presented. The summary and conclusions are given in Section 5.

2. Observation Equation of Multi-GNSS CP Time Transfer

In general, the mathematic model of time transfer based on a CP technique is the same as in the precise point positioning (PPP) model. The ionosphere-free combination for individual GNSSs has typically been utilized for dual-frequency carrier phase and pseudo-range observations to eliminate the first-order effects of the ionosphere. The observation equations can be written as [21,24]:
P = ρ + c ⋅ ( d t r − d t s ) + T t r o p + ε P L = ρ + c ⋅ ( d t r − d t s ) + T t r o p + N I F + ε L
where, P and L are the ionosphere-free combinations of carrier phase and pseudo-range. The geometric distance is denoted as ρ , c is the speed of light, T t r o p is the tropospheric delay, N I F is the phase ambiguity of the ionosphere-free (IF) combination in metres, ε P and ε L are the measurement noise and multipath error for pseudo-range and carrier phase, respectively, d t s is the satellite clock offset and d t r is the receiver clock offset, which is the difference between the external time frequency reference and the GNSS time scale. It is important to note that the phase centre offset and variation, tidal loading, ocean tides, phase wind-up and relativistic delay must be corrected according to the existing models [24], although they are not explicitly expressed in Equation (1).
In the combined GPS/BDS/Galileo time and frequency transfer based on the CP technique, there are three receiver clock offsets that represent differences between the common external time frequency reference and their corresponding system time scales. When the GPS time scale (GPST) is selected as a unique reference time scale to provide the common receiver clock offset based on three CDMA GNSSs, additional inter-system bias (ISB) parameters in BDS and Galileo are introduced, which also refer to inter-system receiver hardware bias. The formula may be expressed as:
{ I S B s y s G C = d t r C − d t r G I S B s y s G E = d t r E − d t r G
The combined observation equation using GPS/BDS/Galileo measurements can therefore be expressed as [18]:
{ P i G = ρ i G + c ⋅ ( d t r G − d t i , s G ) + T t r o p + ε i , P G L i G = ρ i G + c ⋅ ( d t r G − d t i , s G ) + T t r o p + N G + ε i , P G P i C = ρ i C + c ⋅ ( d t r G + I S B s y s G C − d t i , s C ) + T t r o p + ε i , P C L i C = ρ i C + c ⋅ ( d t r G + I S B s y s G C − d t i , s C ) + T t r o p + N C + ε i , P C P i E = ρ i E + c ⋅ ( d t r G + I S B s y s G E − d t i , s E ) + T t r o p + ε i , P E L i E = ρ i E + c ⋅ ( d t r G + I S B s y s G E − d t i , s E ) + T t r o p + N E + ε i , P E
where the superscripts G , C , and E represent the GPS, BDS, and Galileo satellite systems, respectively. The subscript i denotes a single satellite. It is notable that the d t r G becomes a parameter that is jointly determined by GPS, BDS and Galileo satellite systems measurements, which are the quantities of interest when applying precise time and frequency transfer. When the two remote stations m and n equipped with different time and frequency references t m and t n , are used for time and frequency transfer, the common GPST plays a role in the intermediate quantity to achieve time and frequency transfer progress; the corresponding formula can be expressed as follows:
Δ T m , n = d t r G ( m ) − d t r G ( n ) = ( t m − G P S T ) − ( t n − G P S T ) = t m − t n
Equation (4) indicates that the receiver clock offset is a parameter of interest in time and frequency transfer. After the combined observation Equation (3) has been transformed and linearized, one least squares (LS) estimator is employed to determine the unknown parameters in the simplified Equation (5):
X = ( A T P A ) − 1 A T P L
where A is the design matrix of the parameters to be estimated, L denotes the observed value minus the computed value, and P is the corresponding weight matrix. The parameter vector X may be estimated from the expression:
X = [ x , y , z , d t r G , I S B s y s G C , I S B s y s G E , T t r o p , N G , N C , N E ]

3. Robust Helmert Variance Component Estimation (RVCE)

In Equation (5), the weight matrix is the combined weight of the three GNSSs. As the observation qualities of individual GNSSs are not the same due to the different signal structures, frequencies, accuracies of satellite orbit and clock products, signal-to-noise ratios and signal transmission paths, neither individual GNSSs nor the individual satellites should be treated equally [15,19]. Although the existing LS-VCE method can be employed to allocate the weight of different kinds of observation [25], it has lack of robust function. Moreover, in the area of precise time and frequency transfer based on multi-GNSSs, outliers inevitably exist in observations [26,27]. These outliers will further distort the results of receiver clock estimation, and do not reasonably represent the characteristics of the time and frequency references. Therefore, a robust approach based on the Helmert VCE is introduced in this study to allocate the appropriate weights and to control for the influence of outliers, which are generally referred to as the posterior measurement residuals.
The error equation of the three combined GNSSs for time and frequency can be written as:
V = A · X ^ − L = [ V G V C V E ] = [ A G A C A E ] · X ^ − [ L G L C L E ]
where, V denotes the combination posterior measurement residuals, which can be further written as the residuals of the individual systems (i.e., GPS, BDS and Galileo) and X ^ is the estimated vector of unknown parameters. Helmert VCE (i.e., posterior variance stochastic model estimation) is based on the residuals derived from the initial inappropriate weights (i.e., prior equal weights) for the observations. In the practical application of combining the three GNSSs for time and frequency transfer, the range of the GNSS residual values may differ due to their different systematic errors. The individual residuals therefore must be unified to a common range using their corresponding variances and cofactors, the formula for this unification is given as:
V ˜ ( i ) = V ( i ) σ 0 2 Q ( i )
where V ˜ denotes the standard residual matrix. The unit weight variance is σ 0 2 and Q is the corresponding cofactor, which can be determined using Equations (9) and (10), respectively:
σ 0 2 = V T P V n − m
Q = P − 1 − A N − 1 A T
where, n is the number of available GNSS satellites, m denotes the estimated parameter number and N is the matrix of the normal equation. Therefore, according to the Helmert VCE theory, the estimated variance components can be determined using the following formulae:
[ σ ^ 0 , G 2 σ ^ 0 , C 2 σ ^ 0 , E 2 ] = [ s G , G s G , C s G , E s G , C s C , C s C , E s G , E s C , E s E , E ] − 1 [ V ˜ G T P G V ˜ G V ˜ C T P C V ˜ C V ˜ E T P E V ˜ E ]
s k , k = n k − 2 t r ( N − 1 N k ) + t r ( N − 1 N k N − 1 N k ) s k , j = n k − t r ( N − 1 N k N − 1 N j ) N = A T P A }
where, t r is the trace operation for matrix. The initial weights for three GNSSs are typically inappropriate, and the corresponding variance components σ ^ 0 , G 2 , σ ^ 0 , C 2 , and σ ^ 0 , E 2 are also unequal. The relationship between the redefined weight matrix and the estimated components is as follows:
P ^ k = c 0 σ 0 , k 2 P k
where σ 0 , k 2 is the estimated variance component for the three GNSSs, c 0 is constant and usually selected from one of the estimated variance components (i.e., σ 0 , k 2 ), and P k is the initial weight for the three GNSSs. With this new weight matrix, the residuals of the three GNSSs are recalculated until the unit weight variance components of the three systems are equal or have ratios approximately equal to one.
The Helmert VCE is able to obtain an appropriate posterior weight, however, it is not robust to the influence of outliers. In the domains of time and frequency transfer, outliers are inevitably encountered due to code observation noise, equipment failures, and external interference factors [27,28], which can all result in failure of the estimated parameters to reasonably represent the characteristics of the atomic clocks. Therefore, robust estimation [29,30] is introduced to reduce the influence of outliers when combining three GNSSs for time and frequency transfer, this has also been widely applied in geodesy [31,32].
On the basis of the robust estimation theory, the LSQ estimation can be written as follows:
X ¯ = ( A T   P ¯ A ) − 1   A T P ¯ L
where, P ¯ denotes an equivalent weight matrix, which is derived from the distribution of individual GNSS residuals and their equivalent weight function. In this study, the equivalent weight function can be defined by the Institute of Geodesy and Geophysics 3 (IGG3) scheme [33,34], which is constructed with three segments [32,35] based on the following formula:
P i ¯ = { P i , | V ˜ k i | ≤ k 0 P i k 0 | V ˜ k i | ( k 1 − | V ˜ k i | k 1 − k 0 ) 2 , k 0 < | V ˜ k i | ≤ k 1 0 , | V ˜ k i | > k 1
where k 0 and k 1 denote threshold constants depending on the actual distribution, which are typically chosen as k 0 ∈ [ 1.0 ,   1.5 ] and k 1 ∈ [ 2.0 ,   8.0 ] , respectively.

4. Implementation of Multi-GNSS Time and Frequency Transfer using RVCE

Details of the four-step data processing procedure for time and frequency transfer using a RVCE algorithm are shown in Figure 1. First, when combining the three CDMA GNSSs for time and frequency transfer, the observation model can be formed according to Equation (3), and the receiver clock offset is our parameter of interest. Meanwhile, the initial weights of the three GNSSs are temporarily set as equal. The parameters in Equation (6) are then estimated using LSQ adjustment. The residuals of all satellites can be derived according to Equation (7), and may be further standardized by Equation (8). Next, the three variance components are estimated according to Equation (11). It is worth noting that when the robust function is employed, the weighs matrices for the three GNSSs in Equation (11) must be determined by their equivalent weight functions and standardized residuals before the robust variance components may be acquired. Finally, whether the three weight variance components are equal or not must be judged. If equal, the previous weights from the second step are output and used to estimate receiver clock parameters based on Equation (5), thereafter, the progress of time and frequency transfer is completed using Equation (4). If the weight variance components are not equal, the posterior weights are determined based on Equation (11), and then used to begin the next loop in LSQ adjustment until the three weight variance components are equal.

5. Experimental Design

In order to assess the effectiveness of the RVCE approach outlined above, an experiment considering the time link length was devised, the data spanned from MJD 57981 to MJD 57986. Five GNSS stations equipped with stable external time and frequency signals were examined for this experiment. These stations are all able to track the GPS, BDS and Galileo satellite constellations (Figure 2). Two of the stations (NTS1 and BRUX) are located at two time laboratories, the National Time Service Centre (NTSC), and Observatoire Royal de Belgique (ORB)), which are connected with the UTC (NTSC) and UTC (ORB) time and frequency signals, respectively. Another three stations are equipped with external H-maser references. The attributes of these five stations are given in Table 1.
When the time and frequency transfer was employed, the time link was formed between the two stations, which was the primary unit in the area of time and frequency transfer. We therefore established four time links using the five stations with geodetic distances from 883.7 to 7537.5 km, and each link has at least one UTC signal provided by UTC (NTSC) or UTC (ORB), which are summarized in Table 2.

6. Results and Discussion

The performance of time and frequency transfer determined by a CP technique with a RVCE approach (i.e., a RVCE scheme) was compared with the performance of the traditional approach that treats the three GNSSs equally (i.e., the raw scheme). In data processing, the dual-frequency observations (L1/L2 for GPS, B1/B2 for BDS and E1/E5a for Galileo) are used to establish their corresponding ionosphere-free combinations for eliminating the first-order effects of the ionosphere. The combinations for pseudo-range and carrier phase measurement are also called P3 and L3 observation in the area of time and transfer [36,37,38]. The cycle slips of carrier phase are detected by the geometry-free (GF) and Melbourne-Wübbena combination [39], which are performed according to the signal frequency of each satellite system. The satellites with elevations less than 7° are eliminated. Multi-system precise satellite orbit and clock products are obtained from the Multi-GNSS Experiment (MGEX) Center for Orbit Determination in Europe (CODE), in which system time refers to the GPS system scale [40]. The hydrostatic component of the tropospheric delay is corrected using a Saastamoinen model, and the residual non-hydrostatic component is estimated as a random-walk. It is noteworthy that the parameter of interest with respect to the receiver clock offset is estimated by combining pseudo-range and carrier phase observations following a white noise stochastic process. In brief, the receiver clock offsets at adjacent epochs are uncorrelated, and are utilized based on the characteristics of a free-running atomic clock. The details of the observation models and data processing strategies are provided in Table 3. With respect to the RVCE, we used the ratios between pairs of the three estimated variance components as the threshold to control the number of iterations in the experiment. If the ratio value was within the range of 0.9 to 1.1, the iteration was stopped. The unique receiver clock offset parameter was then estimated, and the results of the time transfer were further derived using Equation (4).
Figure 3 preliminarily presents a typical example of receiver clock time series at two stations with and without the RVCE approach. It can be noted the receiver clock value has unobvious discrepancy for two approaches apart from a few epochs. Considering that the time link is the primary analysis unit in practical work of time and frequency transfer, we therefore pay more attention to it in following text.
The results of the time transfer as determined by the multi-GNSS CP technique using the raw and RVCE schemes for all four time links are shown in Figure 4. The RVCE solutions are shown in red and the raw solutions are shown in black. To ensure the clarity of our purpose, the solutions from the RVCE scheme have been translated (10 ns for ONS1-BRUX, 20 ns for NTS1-KAT1, 20 ns for NTS1-CEDU, and −2 ns for BRUX-NTS1).
In the time link of ONS1-BRUX, it can be found that the RVCE scheme has effectively resisted some outliers. Due to the external time and frequency references at the stations, CEDU and KTA1 have not been controlled, while the trend terms are obvious, the performances are difficult to compare. Therefore, the residuals and the root mean square (RMS) value were compared, and the smoothed result is always used because it maintains both low signal noise and the variable trend of the time transfer results. The Kalman smoothing approach was employed in our experiment. Figure 5 shows the residuals of the two different schemes for the four time links, the black line is the traditional solution and the red line is RVCE solution. It can be clearly seen that the RVCE can resist more outliers when compared to the traditional approach. Table 4 presents the corresponding RMS for the time links. The performance of the RVCE solution is slightly improved relative to the raw solution, which ranges from 0.75% to 5.48%, with a mean of 3.41%.
In the domains of time and frequency, the frequency stability is an important index for time links, particularly for short term when using the GNSS technique [42]. Here, the Allan deviation (ADEV) of the time transfer solutions for two different scenarios is obtained for the four time links (Figure 6) and the corresponding ADEV values are given in Table 5. It is clear that the stability of the RVCE solution is greater than that of the raw solution at different average time (AT) τ , particularly when AT are within 10,000 s. Because the number of available data points used for computing ADEV dramatically decreases when τ is greater than 10,000 s, the corresponding measurement errors gradually increase. The error bars at different τ in the Figure 6 also can be agreed with this principle. Therefore, the results that the τ is less than 10,000 s are pay more attention and corresponding improvement for these four time links are summarized in Figure 7. It can be noted that the improvement range relative to the raw solutions differ among the four links. The mean improvement to the ONS1-BRUX link is 13.48%, while improvements are 5.97% for the NTS1-KAT1 link, 32.34% for the NTS1-CEDU link and 7.76% for the BRUX-NTS1 link. The mean value for all four time links is 14.89%.

7. Conclusions

In this study, we analyzed the combined multi-GNSS for CP time and frequency transfer using an RVCE approach, and the corresponding observation model was considered. The RVCE approach and data processing procedure were also detailed. In order to evaluate the performance of the RVCE, four time transfer links with external time and frequency references were computed with and without RVCE. Because that receiver clock estimation is easily affect by the inappropriate weights of multi-GNSSs and outliers from observations, which is closely related to the performance of precise time and frequency transfer. The RVCE approach used in this study is able to determine the observed weight for each GNSS and resist the influence of outliers and shows certain superiority to the raw solution. The smoothed residuals of this method are superior to those of the traditional multi-GNSS CP technique, and the corresponding RMS values were slightly improved for the four time links using RVCE, the range is 5.48% for ONS1-BRUX, 0.75% for NTS1-KAT1, 2.13% for NTS1-CEDU and 5.26% for BRUX-NTS1, respectively. With respect to the frequency stability of time links, the improvements over the short-term within 10,000 s are obvious. The average range is 13.48% for ONS1-BRUX link, 5.97% for NTS1-KAT1 link, 32.34% for NTS1-CEDU link, 7.76% for BRUX-NTS1 link. Therefore, combining multi-GNSS for CP time and frequency transfers using an RVCE approach can improve the validity of the solution within a certain precision range for different time links, particularly for observations that include outliers.

Author Contributions

P.Z., R.T. conceived and designed the experiments; P.Z. performed the experiments, analyzed the data, and wrote the paper under the supervision of R.T.; R.Z., Y.G. and N.L. contributed to discussions and revisions.

Funding

The work was supported by the National Key Research and Development Plan of China (grant no.: 2016YFB0501804), the National Natural Science Foundation of China (grant no.: 41504006, 41674034) and the Chinese Academy of Sciences (CAS) programs of “Pioneer Hundred Talents” (grant no: Y620YC1701) and “Light of the West” (grant no.: Y712YR4701). Further support was granted through the National Time Service Centre (NTSC) programs “Young Creative Talents” (grant no.: Y824SC1S06) and “The Frontier Science Research Project” (grant no.: QYZDB-SSW-DQC028).

Acknowledgments

The authors would like thank the International GNSS Service (IGS) and the Centre for Orbit Determination in Europe (CODE) for providing the GNSS data. Thanks also go to the reviewers for their valuable comments which have improved this paper considerably.

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. Weiss, M.A.; Petit, G.; Jiang, Z. A comparison of GPS common-view time transfer to all-in-view. In Proceedings of the 2005 IEEE International Frequency Control Symposium and Exposition, Vancouver, BC, Canada, 29–31 August 2005. [Google Scholar]
  3. Petit, G.; Jiang, Z. GPS All in View time transfer for TAI computation. Metrologia 2008, 45, 35–45. [Google Scholar] [CrossRef]
  4. Ray, J.; Senior, K. Geodetic techniques for time and frequency comparisons using GPS phase and code measurements. Metrologia 2005, 42, 215–232. [Google Scholar] [CrossRef]
  5. Jiang, Z.; Lewandowski, W. Use of GLONASS for UTC time transfer. Metrologia 2012, 49, 57–61. [Google Scholar] [CrossRef]
  6. Jiang, Z.; Lewandowski, W. Accurate GLONASS Time Transfer for Generation of the Coordinated Universal Time. Int. J. Navig. Obs. 2012, 2012. [Google Scholar] [CrossRef]
  7. Martínez-Belda, M.C.; Defraigne, P. Single-frequency time and frequency transfer with Galileo E5. In Proceedings of the IEEE Frequency Control and the European Frequency and Time Forum, San Francisco, CA, USA, 2–5 May 2011; pp. 1–6. [Google Scholar]
  8. Defraigne, P.; Martínez-Belda, M.C.; Baire, Q. Using Galileo E5 for time and frequency transfer. In Proceedings of the 3rd Galileo Scientific Colloquium, Copenhaguen, Denmark, 30 August–2 September 2011. [Google Scholar]
  9. Huang, W.; Defraigne, P. BeiDou Time Transfer with the Standard CGGTTS. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2016, 63, 1005–1012. [Google Scholar] [CrossRef] [PubMed]
  10. Guang, W.; Dong, S.; Wu, W.; Zhang, J.; Yuan, Y.; Zhang, S. Progress of BeiDou time transfer at NTSC. Metrologia 2018, 55, 175–187. [Google Scholar] [CrossRef]
  11. Dach, R.; Schaer, S.; Hugentobler, U.; Schildknecht, T.; Gäde, A. Combined multi-system GNSS analysis for time and frequency transfer. In Proceedings of the 20th European Frequency and Time Forum, Braunschweig, Germany, 27–30 March 2006; pp. 530–537. [Google Scholar]
  12. Defraigne, P.; Baire, Q. Combining GPS and GLONASS for time and frequency transfer. Adv. Space Res. 2011, 47, 265–275. [Google Scholar] [CrossRef]
  13. 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]
  14. Geng, J.; Zhao, Q.; Shi, C.; Liu, J. A review on the inter-frequency biases of GLONASS carrier-phase data. J. Geod. 2017, 91, 329–340. [Google Scholar] [CrossRef]
  15. Li, X.; Ge, M.; Dai, X.; Ren, X.; Fritsche, M.; Wickert, J.; Schuh, H. Accuracy and reliability of multi-GNSS real-time precise positioning: GPS, GLONASS, BeiDou, and Galileo. J. Geod. 2015, 89, 607–635. [Google Scholar] [CrossRef] [Green Version]
  16. Torre, A.D.; Caporali, A. An analysis of inter-system biases for multi-GNSS positioning. GPS Solut. 2015, 19, 297–307. [Google Scholar] [CrossRef]
  17. Prange, L.; Orliac, E.; Dach, R.; Arnold, D.; Beutler, G.; Schaer, S.; Jäggi, A. CODE’s five-system orbit and clock solution—the challenges of multi-GNSS data analysis. J. Geod. 2017, 91, 345–360. [Google Scholar] [CrossRef]
  18. 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]
  19. Kazmierski, K.; Sośnica, K.; Hadas, T. Quality assessment of multi-GNSS orbits and clocks for real-time precise point positioning. GPS Solut. 2018, 22. [Google Scholar] [CrossRef]
  20. Defraigne, P.; Bruyninx, C. On the link between GPS pseudorange noise and day-boundary discontinuities in geodetic time transfer solutions. GPS Solut. 2007, 11, 239–249. [Google Scholar] [CrossRef]
  21. Petit, G.; Defraigne, P. The performance of GPS time and frequency transfer: comment on ‘A detailed comparison of two continuous GPS carrier-phase time transfer techniques’. Metrologia 2016, 53, 1003–1008. [Google Scholar] [CrossRef] [Green Version]
  22. Knight, N.L.; Wang, J. A comparison of outlier detection procedures and robust estimation methods in GPS positioning. J. Navig. 2009, 62, 699–709. [Google Scholar] [CrossRef]
  23. Zair, S.; Hégarat-Mascle, S.L.; Seignez, E. Outlier Detection in GNSS Pseudo-Range/Doppler Measurements for Robust Localization. Sensors 2016, 16, 580. [Google Scholar] [CrossRef] [PubMed]
  24. Kouba, J. A guide to using International GNSS Service (IGS) Products. Available online: https://igscb.jpl.nasa.gov/igscb/resource/pubs/UsingIGSProductsVer21.pdf (accessed on 15 August 2018).
  25. Gao, Z.; Shen, W.; Zhang, H.; Ge, M.; Niu, X. Application of Helmert Variance Component Based Adaptive Kalman Filter in Multi-GNSS PPP/INS Tightly Coupled Integration. Remote Sens. 2016, 8, 553. [Google Scholar] [CrossRef]
  26. Guo, H.R. Study on the Analysis Theories and Algorithms of the Time and Frequency Characterization for Atomic Clocks of Navigation Satellites. Ph.D. Thesis, Information Engineering University, Zhengzhou, China, October 2006. [Google Scholar]
  27. Guo, H.R.; Yang, Y.X. Analyses of main error sources on time domain frequency stability for atomic clocks of navigation satellites. Geomat. Inf. Sci. Wuhan Univ. 2009, 34, 218–221. [Google Scholar]
  28. Huang, G.; Zhang, Q. Real-time estimation of satellite clock offset using adaptively robust Kalman filter with classified adaptive factors. GPS Solut. 2012, 16, 531–539. [Google Scholar] [CrossRef]
  29. Ou, J.K. Design of a new scheme of robust estimation by three steps. Acta Geod. Cartogr. Sin. 1996, 25, 173–179. [Google Scholar]
  30. Zhou, J.W. Classical theory of errors and robust estimation. Acta Geod. Cartogr. Sin. 1989, 18, 115–120. [Google Scholar]
  31. Kubik, K. An error theory for the Danish method. In Proceedings of the International Symposium on Numerical Methods for Photogrammetric Mapping, Helsinki, Finland, 26–29 September 1982; pp. 299–309. [Google Scholar]
  32. Yang, Y.X. Robust estimation for dependent observations. Manuscr. Geod. 1994, 19, 10–17. [Google Scholar]
  33. Yang, Y.X.; Song, L.J.; Xu, T.H. Robust estimator for correlated observations based on bifactor equivalent weights. J. Geod. 2002, 76, 353–358. [Google Scholar] [CrossRef]
  34. Yang, Y.X.; Song, L.J.; Xu, T.H. Robust parameter estimation for geodetic correlated observation. Acta Geod. Cartogr. Sin. 2002, 31, 95–99. [Google Scholar]
  35. Yang, Y.X.; Cheng, M.K.; Shum, C.K.; Tapley, B.D. Robust estimation of systematic errors of satellite laser range. J. Geod. 1999, 73, 345–349. [Google Scholar] [CrossRef]
  36. Defraigne, P.; Banerjee, P.; Lewandowski, W. Time transfer through GPS. Indian J. Radio Space Phys. 2007, 36, 303–312. [Google Scholar]
  37. Defraigne, P. GNSS Time Transfer. Available online: https://www.bipm.org/ws/CCTF/TAI_TRAINING/Allowed/Fundamentals/Training-2012-GNSS-Defraigne.pdf (accessed on 19 August 2018).
  38. Guidelines for Providing Data for PPP links. Available online: ftp://ftp2.bipm.org/pub/tai/timelinks/taippp/ppp_for_tai_guidelines_v2.pdf (accessed on 19 August 2018).
  39. Hatch, R. The synergism of GPS code and carrier measurements. In Proceedings of the Third International Symposium on Satellite Doppler Positioning, Las Cruces, NM, USA, 8–12 February 1982; pp. 1213–1231. [Google Scholar]
  40. Prange, L.; Dach, R.; Lutz, S.; Schaer, S.; Jäggi, A. The CODE solution for the IGS MGEX. In Proceedings of the IGS workshop, Pasadena, CA, USA, 23–27 June 2014. [Google Scholar]
  41. Niell, A.E. Global mapping functions for the atmosphere delay at radio wavelengths. J. Geophys. Res. 1996, 101, 3227–3246. [Google Scholar] [CrossRef]
  42. Jiang, Z.; Petit, G. Combination of TWSTFT and GNSS for accurate UTC time transfer. Metrologia 2009, 46, 305–314. [Google Scholar] [CrossRef]
Sensors 18 02878 g001 550
Figure 1. Primary implementation procedure using robust Helmert variance component estimation (RVCE) for time and frequency transfer.
Figure 1. Primary implementation procedure using robust Helmert variance component estimation (RVCE) for time and frequency transfer.
Sensors 18 02878 g001
Sensors 18 02878 g002 550
Figure 2. Layouts of multi-GNSS stations, capable of tracking GPS, BDS and Galileo constellations.
Figure 2. Layouts of multi-GNSS stations, capable of tracking GPS, BDS and Galileo constellations.
Sensors 18 02878 g002
Sensors 18 02878 g003 550
Figure 3. Receiver clock time series for station ONS1 and CEDU with and without the RVCE approach.
Figure 3. Receiver clock time series for station ONS1 and CEDU with and without the RVCE approach.
Sensors 18 02878 g003
Sensors 18 02878 g004 550
Figure 4. Results of time transfers determined by a carrier-phase (CP) technique with and without the RVCE approach for the four time links. (a) is for ONS1-BRUX, (b) is for NTS1-KAT1, (c) is for NTS1-CEDU and (d) is for BRUX-NTS1.
Figure 4. Results of time transfers determined by a carrier-phase (CP) technique with and without the RVCE approach for the four time links. (a) is for ONS1-BRUX, (b) is for NTS1-KAT1, (c) is for NTS1-CEDU and (d) is for BRUX-NTS1.
Sensors 18 02878 g004
Sensors 18 02878 g005 550
Figure 5. Smoothing residuals of the results with and without the RVCE approach for the four time links. (a) is for ONS1-BRUX, (b) is for NTS1-KAT1, (c) is for NTS1-CEDU and (d) is for BRUX-NTS1.
Figure 5. Smoothing residuals of the results with and without the RVCE approach for the four time links. (a) is for ONS1-BRUX, (b) is for NTS1-KAT1, (c) is for NTS1-CEDU and (d) is for BRUX-NTS1.
Sensors 18 02878 g005
Sensors 18 02878 g006 550
Figure 6. Comparison of Allan deviation (ADEV) for the four time links determined with and without RVCE. (a) is for ONS1-BRUX, (b) is for NTS1-KAT1, (c) is for NTS1-CEDU and (d) is for BRUX-NTS1.
Figure 6. Comparison of Allan deviation (ADEV) for the four time links determined with and without RVCE. (a) is for ONS1-BRUX, (b) is for NTS1-KAT1, (c) is for NTS1-CEDU and (d) is for BRUX-NTS1.
Sensors 18 02878 g006
Sensors 18 02878 g007 550
Figure 7. Improvement to stability (%) of RVCE solutions over the raw solutions within 10,000 s.
Figure 7. Improvement to stability (%) of RVCE solutions over the raw solutions within 10,000 s.
Sensors 18 02878 g007
Table
Table 1. Attributes of stations used in this study.
Table 1. Attributes of stations used in this study.
Site NameGNSS ReceiverAntennaFrequency Standard
NTS1SEPT POLARX4TRSEPCHOKE_MCUTC (NTSC)
BRUXSEPT POLARX4TRJAVRINGANT_DMUTC (ORB)
KAT1SEPT POLARX5LEIAR25.R3EXTERNAL H-MASER
CEDUSEPT POLARX5AOAD/M_TEXTERNAL H-MASER
ONS1TRIMBLE NETR9LEIAR25.R3EXTERNAL H-MASER
Table
Table 2. Time links formed and their geodetic distances.
Table 2. Time links formed and their geodetic distances.
Time Link NameGeodetic Distance (km)
ONS1-BRUX883.7
NTS1-KAT15704.2
NTS1-CEDU7294.5
BRUX-NTS17537.5
Table
Table 3. Observation models and multi-system RINEX data processing strategies.
Table 3. Observation models and multi-system RINEX data processing strategies.
ItemsModels and Strategies
ObservationsUndifferenced carrier phase and code observation
Signal selectionGPS:L1/L2; BeiDou: B1/B2; Galileo: E1/E5a
Satellite orbit and clockUsing the precise satellite products from CODE [40]
Satellite antenna phase centerCorrected using MGEX value
Initial weight between code and phase0.0001:1
Ionosphere Eliminated by ionosphere-free combination [36,37,38]
Tropospheric delayInitial model + random-walk process
Tropospheric mapping functionNeill mapping function (NMF) [41]
Elevation cutoff7°
Sampling rate30 s
Observation weightElevation dependent weight
EstimatorLS in sequential mode
Receiver clock offsetEstimated with white noise
ISBEstimated with epoch-wise variable method
Phase wind-up effectModel corrected
Table
Table 4. Improvement (%) in RMS of the result smoothing in two scenarios for the time links.
Table 4. Improvement (%) in RMS of the result smoothing in two scenarios for the time links.
Time Link NameRaw (ns)RVCE (ns)Improvement (%)
ONS1-BRUX0.4200.3975.48
NTS1-KAT10.1330.1320.75
NTS1-CEDU0.1410.1382.13
BRUX-NTS10.0190.0185.26
Table
Table 5. ADEV values of the four time links determined with and without RVCE.
Table 5. ADEV values of the four time links determined with and without RVCE.
τ
(seconds)		Time Link Name
ONS1-BRUXNTS1-KAT1NTS1-CEDUBRUX-NTS1
RawRVCERawRVCERawRVCERawRVCE
301.32 × 10−121.06 × 10−125.10 × 10−134.97 × 10−139.38 × 10−136.69 × 10−134.08 × 10−133.82 × 10−13
601.02 × 10−128.17 × 10−133.69 × 10−133.55 × 10−136.63 × 10−134.65 × 10−132.58 × 10−132.37 × 10−13
1208.27 × 10−136.46 × 10−132.69 × 10−132.55 × 10−134.75 × 10−133.30 × 10−131.60 × 10−131.46 × 10−13
2406.98 × 10−135.79 × 10−132.06 × 10−131.95 × 10−133.54 × 10−132.44 × 10−131.03 × 10−139.68 × 10−14
4805.48 × 10−134.79 × 10−131.50 × 10−131.40 × 10−132.90 × 10−131.94 × 10−138.00 × 10−147.10 × 10−14
9604.14 × 10−133.60 × 10−131.07 × 10−139.90 × 10−141.74 × 10−131.19 × 10−135.95 × 10−145.48 × 10−14
19203.03 × 10−132.60 × 10−137.88 × 10−147.22 × 10−141.62 × 10−131.04 × 10−134.33 × 10−143.80 × 10−14
38401.75 × 10−131.74 × 10−135.63 × 10−145.25 × 10−141.36 × 10−138.91 × 10−142.83 × 10−142.67 × 10−14
76807.92 × 10−147.72 × 10−144.42 × 10−144.13 × 10−141.05 × 10−136.76 × 10−142.26 × 10−142.19 × 10−14
15,3606.74 × 10−146.77 × 10−143.80 × 10−143.61 × 10−143.39 × 10−142.53 × 10−141.90 × 10−141.57 × 10−14
30,7202.11 × 10−142.17 × 10−142.46 × 10−142.58 × 10−141.82 × 10−141.95 × 10−141.22 × 10−141.19 × 10−14
61,4401.62 × 10−141.78 × 10−141.72 × 10−141.76 × 10−141.52 × 10−141.55 × 10−141.56 × 10−141.37 × 10−14

Share and Cite

MDPI and ACS Style

Zhang, P.; Tu, R.; Gao, Y.; Zhang, R.; Liu, N. Improving the Performance of Multi-GNSS Time and Frequency Transfer Using Robust Helmert Variance Component Estimation. Sensors 2018, 18, 2878. https://doi.org/10.3390/s18092878

AMA Style

Zhang P, Tu R, Gao Y, Zhang R, Liu N. Improving the Performance of Multi-GNSS Time and Frequency Transfer Using Robust Helmert Variance Component Estimation. Sensors. 2018; 18(9):2878. https://doi.org/10.3390/s18092878

Chicago/Turabian Style

Zhang, Pengfei, Rui Tu, Yuping Gao, Rui Zhang, and Na Liu. 2018. "Improving the Performance of Multi-GNSS Time and Frequency Transfer Using Robust Helmert Variance Component Estimation" Sensors 18, no. 9: 2878. https://doi.org/10.3390/s18092878

APA Style

Zhang, P., Tu, R., Gao, Y., Zhang, R., & Liu, N. (2018). Improving the Performance of Multi-GNSS Time and Frequency Transfer Using Robust Helmert Variance Component Estimation. Sensors, 18(9), 2878. https://doi.org/10.3390/s18092878

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