Multi-GNSS satellite clock estimation constrained with oscillator noise model in the existence of data discontinuity

During the past years, real-time precise point positioning has been proven to be an efficient tool in the applications of navigation, precise orbit determination of LEO as well as earthquake and tsunami early warning, etc. One of the most crucial issues of these applications is the high-precision real-time GNSS satellite clock. Though the performance and character of the GNSS onboard atomic frequency standard have been widely studied, the white noise model is still the most popular hypothesis that employed in the real-time GNSS satellite clock estimation. However, concerning the real-time applications, significant data discontinuity may arise either due to the fact that only regional stations involved, or the failure in the stations, satellites and network connections. These data discontinuity would result in an arbitrary clock jump between adjacent arcs when the clock offsets are modeled as white noise. In addition, it is also expected that the detection and identification of outliers would be benefited from the constrains of the satellite oscillator noise model. Thus in this contribution, based on the statistic analysis of almost 2-year multi-GNSS precise clock products, we developed the oscillator noise model for the satellites of GPS, GLONASS, BDS and Galileo according to the oscillator type as well as the block type. Then, the efficiency of this oscillator noise model in multi-GNSS satellite clock estimation is demonstrated with 2-months data for both regional and global networks in simultaneous real-time mode. For the regional network, the results suggest that compared with the traditional solution based on white noise model, the improvement is 44.4 and 12.1% on average for STD and RMS, respectively, and the improvement is mainly attributed to the efficiency of the oscillator noise model during the convergence period and the gross error resistance. Concerning the global experiment, since the stations guarantee the continuous tracking of the satellites with redundant observable, the improvement is not as evident as that of regional experiment for GPS, GLONASS and BDS. The STD of Galileo clock improves from 0.28 to 0.19 ns due to that, the satellites E14 and E18 still suffer significant data discontinuity during our experimental period.