Switching and performance variations of on-orbit BDS satellite clocks

The information of the satellite clock switching and performance variations on-orbit of Chinese BeiDou-2 Navigation System (BDS) is not available for the public. In order to detect the BDS satellite clock switching and performances variation, we analyzed the precise clock offset products with a total duration of 5?years every BDS satellite equipped four atomic clocks from four different manufactures from January 2013 to October 2017. Three important contributions are concluded as follows. (1) It is found that the average time of on-orbit operation for BDS satellite clocks is about 1–2?years. There have been 22 times of clock switching for BDS satellites, of which the C05 and C08 satellites have been switched to the fourth (last) atomic clock. (2) There are frequent phase adjustments for BDS on-orbit satellite clocks, and the frequency series is relatively stable. Furthermore, there are semi-annual sinusoid cycles in the frequency drift series of C06 and C09 satellites. (3) The performances of MEO satellite clocks perform better than the IGSO and GEO satellite clocks. The average ten-thousand frequency stability of BDS satellite clocks is about 1E-13, which is worse than that of GPS and Galileo but better than that of GLONASS.     

A new stochastic model considering satellite clock interpolation errors in precise point positioning

Precise clock products are typically interpolated based on the sampling interval of the observational data when they are used for in precise point positioning. However, due to the occurrence of white noise in atomic clocks, a residual component of such noise will inevitable reside within the observations when clock errors are interpolated, and such noise will affect the resolution of the positioning results. In this paper, which is based on a twenty-one-week analysis of the atomic clock noise characteristics of numerous satellites, a new stochastic observation model that considers satellite clock interpolation errors is proposed. First, the systematic error of each satellite in the IGR clock product was extracted using a wavelet de-noising method to obtain the empirical characteristics of atomic clock noise within each clock product. Then, based on those empirical characteristics, a stochastic observation model was structured that considered the satellite clock interpolation errors. Subsequently, the IGR and IGS clock products at different time intervals were used for experimental validation. A verification using 179 stations worldwide from the IGS showed that, compared with the conventional model, the convergence times using the stochastic model proposed in this study were respectively shortened by 4.8% and 4.0% when the IGR and IGS 300-s-interval clock products were used and by 19.1% and 19.4% when the 900-s-interval clock products were used. Furthermore, the disturbances during the initial phase of the calculation were also effectively improved.     

基于实时SSR数据的LEO卫星精密单点定位研究

对目前低轨卫星实时定位的方法进行了研究,现在通常采用GPS定位,使用广播星历和普通晶振,实时定位精度一般在10m以内,不能满足高精度实时定位的需求。IGS组织在全球范围内对GPS跟踪分析,生成精密星历、精密钟差产品、按SSR格式的广播星历和钟差修正产品并在网上发布。对这些IGS产品进行了调查,提出在现有测控支持情况下,可以通过高密度上注SSR信息流实现在轨高精度定位。以某型号低轨微小卫星在轨导航增强载荷为应用背景,用IGS03产品中的1057和1058数据对双频GPS接收机的星历和钟差进行修正,采用递推最小二乘估计和LAMDA模糊度固定过对载波相位和伪距信息进行处理,在短时间内获得亚米级定位结果。     

Enhancing the kinematic precise orbit determination of low earth orbiters using GPS receiver clock modelling

Clock error estimation has been the focus of a great deal of research because of the extensive usage of clocks in GPS positioning applications. The receiver clock error in the spacecraft orbit determination is commonly estimated on an epoch-by-epoch basis, along with the spacecraft’s position. However, due to the high correlation between the spacecraft orbit altitude and the receiver clock parameters, estimates of the radial component are degraded in the kinematic approach. Using clocks with high stability, the predictable behaviour of the receiver oscillator can be exploited to improve the positioning accuracy, especially for the radial component. This paper introduces two GPS receiver clock models to describe the deterministic and stochastic property of the receiver clock, both of which can improve the accuracy of kinematic orbit determination for spacecraft in low earth orbit. In particular, the clock parameters are estimated as time offset and frequency offset in the two-state model. The frequency drift is also estimated as an unknown parameter in the three-state model. Additionally, residual non-deterministic random errors such as frequency white noise, frequency random walk noise and frequency random run noise are modelled. Test results indicate that the positioning accuracy could be improved significantly using one day of GRACE flight data. In particular, the error of the radial component was reduced by over 40.0% in the real-time scenario.     

Analysis of GNSS clock prediction performance with different interrupt intervals and application to real-time kinematic precise point positioning

Continuous and timely real-time satellite orbit and clock products are mandatory for real-time precise point positioning (RT-PPP). Real-time high-precision satellite orbit and clock products should be predicted within a short time in case of communication delay or connection breakdown in practical applications. For prediction, historical data describing the characteristics of the real-time orbit and clock can be used as the basis for performing the prediction. When historical data are scarce, it is difficult for many existing models to perform precise predictions. In this paper, a linear regression model is used to predict clock products. Seven-day GeoForschungsZentrum (GFZ) final clock products sampled at 30 s are used to analyze the characteristics of GNSS clocks. It is shown that the linear regression model can be used as the prediction model for the satellite clock products. In addition, the accuracy of the clock prediction for different satellites are analyzed using historical data with different periods (such as 2 and 10 epochs). Experimental results show that the accuracy of the clock with the linear regression prediction model using historical data with 10 epochs is 1.0 ns within 900 s. This is higher accuracy than that achieved using historical data of 2 epochs. Finally, the performance analysis for real-time kinematic precise point positioning (PPP) is provided using GFZ final clock prediction results and state space representation (SSR) clock prediction results when communication delay or connection breakdown occur. Experimental results show that the positioning accuracy without prediction is better than that with prediction in general, whether using the final clock product or the SSR clock product. For the final clock product, the positioning accuracy in the north (N), east (E), and up (U) directions is better than 10.0 cm with all visible GNSS satellites with prediction. In comparison, the 3D positioning accuracy of N, E, and U directions with visible GNSS satellites whose prediction accuracy is better than 0.1 ns using historical data of 10 epochs is improved from 15.0 cm to 7.0 cm. For the SSR clock product, the positioning accuracy of N, E, and U directions is better than 12.0 cm with visible GNSS satellites with prediction. In comparison, the 3D positioning accuracy of N, E, and U directions with visible GNSS satellites whose prediction accuracy is better than 0.1 ns using historical data of 10 epochs is improved from 12.0 cm to 9.0 cm.     

GPS授时校频方法研究与试验结果

为了解决多目标综合测量系统各测站之间时间同步和频率校准问题,提出了利用GPS(Global Positioning System)单星或多星共视方法进行站间时间同步与校频,给出了这两种方法的计算公式,分析了星历误差、星钟误差、电离层折射误差、对流层折射误差、多径效应和接收机硬件延迟对时间同步精度的影响.为了验证GPS授时校频精度,进行了相关试验.通过与铯原子钟比对,表明利用GPS可实现纳秒级时间同步,校频精度也优于5.0×10-11,多星共视具有更高的同步校频精度.      

The GLONASS-M satellite yaw-attitude model

The proper modeling of the satellites’ yaw-attitude is a prerequisite for high-precision Global Navigation Satellite System (GNSS) positioning and poses a particular challenge during periods when the satellite orbital planes are partially eclipsed. Whereas a lot of effort has been put in to examine the yaw-attitude control of GPS satellites that are in eclipsing orbits, hardly anything is known about the yaw-attitude behavior of eclipsing GLONASS-M satellites. However, systematic variations of the carrier phase observation residuals in the vicinity of the orbit’s noon and midnight points of up to ±27 cm indicate significant attitude-related modeling issues. In order to explore the GLONASS-M attitude laws during eclipse seasons, we have studied the evolution of the horizontal satellite antenna offset estimates during orbit noon and orbit midnight using a technique that we refer to as “reverse kinematic precise point positioning”. In this approach, we keep all relevant global geodetic parameters fixed and estimate the satellite clock and antenna phase center positions epoch-by-epoch using 30-second observation and clock data from a global multi-GNSS ground station network. The estimated horizontal antenna phase center offsets implicitly provide the spacecraft’s yaw-attitude. The insights gained from studying the yaw angle behavior have led to the development of the very first yaw-attitude model for eclipsing GLONASS-M satellites. The derived yaw-attitude model proves to be much better than the nominal yaw-attitude model commonly being used by today’s GLONASS-capable GNSS software packages as it reduces the observation residuals of eclipsing satellites down to the normal level of non-eclipsing satellites and thereby prevents a multitude of measurements from being incorrectly identified as outliers. It facilitates continuous satellite clock estimation during eclipse and improves in particular the results of kinematic precise point positioning of ground-based receivers.     

基于GPS/北斗组合系统的高轨卫星定位技术研究

针对目前高轨GPS信号可用性差及定位精度低的特点, 对GPS/北斗组合系统的 高轨卫星定位技术进行研究, 对比分析了单GPS系统与GPS/北斗组合系统的卫 星可见性和几何精度因子. 结果表明, GPS/北斗组合系统比单GPS系统的卫星可 见性好, 且定位精度高. 同时通过提出在星载接收机上采用高精度原子钟, 可实现三星定位, 降低对接收机的技术要求.      

时钟偏差辅助的GPS完整性监测算法

对伪距残差最小二乘的接收机自主完整性监测(RAIM)算法进行了分析.在此基础上,通过对接收机时钟偏差的建模,将模型预测的时钟偏差引入伪距残差最小二乘的RAIM算法中,保证了在可见星仅为4颗的情况下,仍能利用χ2检验法对全球定位系统(GPS)进行完整性监测,从而达到提高完整性监测算法有效性的目的.计算机仿真的结果显示,这种辅助方式不仅计算简单,而且有效可行.      

导航卫星原子钟舱温度控制方法及其验证

原子钟是导航卫星的重要组成部分,可为卫星系统提供高准确及高稳定度的时间频率源.原子钟工作性能与环境温度变化密切相关,为保证其在轨连续、稳定运行,热控系统需为其提供良好的工作温度环境.本文以某导航卫星原子钟舱温度控制为研究内容,给出原子钟舱热控设计方案、控制算法,并进行仿真分析和试验验证.在轨遥测数据表明,卫星原子钟舱热控方案和控制算法设计合理,仿真分析及试验结果有效,各原子钟在轨工作温度满足要求,原子钟温度稳定度满足并优于设计指标近一个数量级.     

BeiDou-3 orbit and clock quality of the IGS Multi-GNSS Pilot Project

Within the Multi-GNSS Pilot Project (MGEX) of the International GNSS Service (IGS), precise orbit and clock products for the BeiDou-3 global navigation satellite system (BDS-3) are routinely generated by a total of five analysis centers. The processing standards and specific properties of the individual products are reviewed and the BDS-3 orbit and clock product performance is assessed through direct inter-comparison, satellite laser ranging (SLR) residuals, clock stability analysis, and precise point positioning solutions. The orbit consistency evaluated by the signal-in-space range error is on the level of 4–8 cm for the medium Earth orbit satellites whereas SLR residuals have RMS values between 3 and 9 cm. The clock analysis reveals sytematic effects related to the elevation of the Sun above the orbital plane for all ACs pointing to deficiencies in solar radiation pressure modeling. Nevertheless, precise point positioning with the BDS-3 MGEX orbit and clock products results in 3D RMS values between 7 and 8 mm.     

Early analysis of precise orbit and clock offset determination for the satellites of the global BeiDou-3 system

Eight new-generation BeiDou satellites (BeiDou-3) have been launched into Medium Earth Orbit (MEO), allowing for global coverage since March 2018, and they are equipped with new hydrogen atomic clocks and updated rubidium clocks. Firstly, we analyzed the signals for the carrier-to-noise-density ratio (C/N0) and pseudorange multipath (MP) by using international GNSS (Global Navigation Satellite System) Monitoring and Assessment System (iGMAS) station data, and found that B1C has a lower C/N0, and B2a has the same level of C/N0 as the B1I and B3I signals. For pseudorange multipath, compared with the BeiDou-2 satellites, the obvious systematic variation of MP scatters related to the elevation angle is greatly improved for the BeiDou-3 and BeiDou-3e satellites signals. For the signals of the BeiDou-3 satellites, the order of the Root Mean Square (RMS) values of multipath and noise is B3I?<?B1I?<?B2a?<?B1C. Then, the comparison of the precise orbit determination and clock offset determination for the BeiDou-2, BeiDou-3, and BeiDou-3 experimental (BeiDou-3e) satellites was done by using 10 stations from iGMAS. The 3D precision of the 24?h orbit overlap is 24.55, 25.61, and 23.35?cm for the BeiDou-3, BeiDou-3e, and BeiDou-2 satellites, respectively. BeiDou-3 satellite has a comparable precision to that of the BeiDou-2 satellite. For the precision of clock offset estimation, the Standard Deviation (STD) of the BeiDou-3 MEO satellite is 0.350?ns, which is an improvement of 0.042?ns over that of the BeiDou-2 MEO satellite. The stabilities of the BeiDou-3 and BeiDou-3e onboard clocks are better than those of BeiDou-2 by factors of 2.84 and 1.61 at an averaging time of 1000 and 10,000?s, respectively.     

基于三频数据的北斗卫星导航系统DCB参数精度评估方法

差分码偏差（Differential Code Biases,DCB）参数作为导航电文中重要的一项,是影响用户PNT服务的主要误差源之一。北斗卫星导航系统（后文简称“北斗系统”）发射三个频点的导航信号,在导航电文中需要发播卫星的2个TGD（Timing Group Delay）参数。文章首先介绍了北斗系统卫星DCB参数最小二乘解算与形式误差评估;其次根据北斗系统三频特点,提出了不同频点组合计算垂直方向电离层电子总含量（VTEC）互差的DCB精度定量评估方法,并与IGS(International GNSS Service)提供的GPS卫星DCB精度进行比较;最后,详细分析了DCB参数精度对用户等效距离误差（UERE）计算和定位计算的影响,分别采用卫星出场标定DCB参数和经过解算DCB参数进行评估。实测数据分析结果表明,北斗系统卫星DCB参数解算形式误差与IGS解算GPS卫星DCB参数形式误差相当,但受卫星类型和解算测站的几何分布限制,北斗系统卫星DCB参数解算不确定度相比IGS略差,估计精度优于0.5ns,不同频率组合计算VTEC互差绝对值均值优于0.6TECU。相比采用卫星出场标定值,采用系统解算DCB参数后,双频用户三维位置误差改善13.80%～47.42%。     

氢钟在守时钟组配置中的应用研究

氢钟因其良好的低噪声特性,被广泛应用于高精度时间/ 频率测量和比对。对于物理信号噪声特性要求较高的时间实验室,氢钟已经成为守时钟组不可或缺的频率源。然而由于物理结构和工作方式等多方面的原因,氢钟具有明显的频率漂移特性,此外氢钟的价格较铯钟更为昂贵,因此氢钟的数量成为守时钟组配置中急需解决的核心问题。本文首先分析了氢钟在守时钟组配置中的必要性,然后通过建立守时钟组优化配置数学模型,利用多个时间实验室的原子钟数据资料,采用统计学的方法,提出守时钟组的多种配置方式。并通过分析比较不同配置方式下的守时钟组性能,提出稳定度一定的条件下,中小型时间实验室的守时钟组最优配置方案。     

GPS P(Y)码功率增强监测分析

全球定位系统(GPS)的Block IIR-M和Block IIF卫星具备可编程功率输出能力, 可以灵活增强单个信号分量的发射功率。为了系统评估GPS P(Y)码的功率增强能力, 对弹性功率的原理进行了理论分析, 提出了GPS信号功率增强的监测分析方法, 利用国际GNSS监测评估系统(iGMAS)和国际GNSS服务(IGS)监测站数据、高增益天线监测数据、事后精密星历对GPS增强P(Y)码的覆盖性及星座性能、空间信号和用户端性能进行了分析。结果表明:在保持发射总功率和民用信号功率不变的情况下, Block IIF和Block IIR-M卫星的L1 P(Y)码和L2 P(Y)码功率相比正常水平分别增强约6 dB和5 dB;在功率增强信号覆盖区内仅利用19颗增强卫星进行双频单点定位, 位置误差不大于15 m(95%);当可见增强卫星数为6, 增强后的P(Y)码载噪比为55 dB·Hz时, P(Y)码之间的多址干扰引起的等效载噪比下降量为0.4 dB。      

Beidou wide-area augmentation system clock error correction and performance verification

On December 27, 2018, the Beidou-3 System (BDS-3) has completed the deployment of 18 Medium Earth Orbit (MEO) satellites combined as a space constellation. In addition to the augmentation information for the new system signals B1C and B2a, the BDS-3 is compatible with the three augmentation information broadcast by the BDS-2 system for B1I and B3I signals: equivalent clock error correction, User Differential Ranging Error (UDRE) and Grid Ionosphere Vertical Error (GIVE). In this paper, the observation data of Beidou monitoring network are used to analyze the pseudo-range observation quality of the smooth transition signals B1I and B3I of BDS-2 and BDS-3. At the same time, the relationship between the equivalent clock error correction and the prediction error of the satellite clock is analyzed by using the Two-Way Satellite Time and Frequency Transfer (TWSTFT) data. The results show that the correlation between the equivalent clock error correction and the monitored clock error by using the TWSTFT data is greater than 60%. We calculate the UDRE by using the equivalent clock error correction. The results show that the satellite equivalent clock error correction can improve the accuracy of User Equivalent Range Error (UERE) by about 50%. This paper also compares the positioning accuracy of the BDS-2 satellites with the BDS-2 satellites combined BDS-3 satellites. The results show that the three-dimensional positioning accuracy is improved by about 30% after the BDS-3 satellites are added.     

Network based real-time precise point positioning

Network based real-time precise point positioning system includes two stages, i.e. real-time estimation of satellite clocks based on a reference network and real-time precise point positioning thereafter. In this paper, a satellite- and epoch-differenced approach, adopted from what is introduced by Han et al. (2001), is presented for the determination of satellite clocks and for the precise point positioning. One important refinement of our approach is the implementation of the robust clock estimation. A prototype software system is developed, and data from the European Reference Frame Permanent Network on September 19, 2009 is used to evaluate the approach. Results show that our approach is 3 times and 90 times faster than the epoch-difference approach and the zero-difference approach, respectively, which demonstrates a significant improvement in the computation efficiency. The RMS of the estimated clocks is at the level of 0.1 ns (3 cm) compared to the IGS final clocks. The clocks estimates are then applied to the precise point positioning in both kinematic and static mode. In static mode, the 2-h estimated coordinates have a mean accuracy of 3.08, 5.79, 6.32 cm in the North, East and Up directions. In kinematic mode, the mean kinematic coordinates accuracy is of 4.63, 5.82, 9.20 cm.     

Precise orbit determination and precision comparison for FY3C and FY3D with receiver antenna GPS and BDS PCV using spaceborne BDS and GPS observation data

The FY3C and FY3D satellites were equipped with global navigation satellite occultation detector (GNOS) receivers that received both GPS and BDS-2 signals. For further improving precise orbit determination (POD) precisions, we estimated receiver GPS and BDS signal phase center variations (PCV) models with 2° and 5° resolutions and set the different weights for GPS and BDS-2 observations in the combined POD. The BDS-based POD precision using BDS-2 satellite antenna phase center offset (PCO) values from the China Satellite Navigation Office (CSNO) are not as accurate as those obtained from the International GNSS Service (IGS) Multi-GNSS experiments project (MGEX). The estimated receiver GPS and BDS PCV models with 2° and 5° resolutions were estimated from the GPS phase residuals of GPS-based POD and BDS phase residuals of combined POD, respectively. In most cases, the POD precisions using the estimated PCVs with 2° resolution are superior to those with 5° resolution. The precisions of the BDS-based POD and combined POD were both improved by introducing the receiver BDS PCV models. The weighting for GPS and BDS-2 observations can further improve the precision of the combined POD. The tested results of selected weights are better than those with equal weight in the combined POD. The experiment results show that orbital precisions of FY3C are worse than those of FY3D.     

Orbit determination performances using single- and double-differenced methods: SAC-C and KOMPSAT-2

In this paper, Global Positioning System-based (GPS) Orbit Determination (OD) for the KOrea-Multi-Purpose-SATellite (KOMPSAT)-2 using single- and double-differenced methods is studied. The requirement of KOMPSAT-2 orbit accuracy is to allow 1 m positioning error to generate 1-m panchromatic images. KOMPSAT-2 OD is computed using real on-board GPS data. However, the local time of the KOMPSAT-2 GPS receiver is not synchronized with the zero fractional seconds of the GPS time internally, and it continuously drifts according to the pseudorange epochs. In order to resolve this problem, an OD based on single-differenced GPS data from the KOMPSAT-2 uses the tagged time of the GPS receiver, and the accuracy of the OD result is assessed using the overlapping orbit solution between two adjacent days. The clock error of the GPS satellites in the KOMPSAT-2 single-differenced method is corrected using International GNSS Service (IGS) clock information at 5-min intervals. KOMPSAT-2 OD using both double- and single-differenced methods satisfies the requirement of 1-m accuracy in overlapping three dimensional orbit solutions. The results of the SAC-C OD compared with JPL’s POE (Precise Orbit Ephemeris) are also illustrated to demonstrate the implementation of the single- and double-differenced methods using a satellite that has independent orbit information available for validation.     

Monitoring atomic clocks on board GNSS satellites

A method for monitoring atomic clocks on board Global Navigation Satellites System (GNSS) satellites is described to address the issue of clock related signal integrity in safety–critical applications of GNSS. The carrier-phase time transfer is employed in the clock monitoring method which enables tight tracking of the satellite onboard clocks and thus improves detectability of clock anomalies. Detecting onboard clock anomalies requires the ability to monitor clocks in real time, and a Kalman filter can then be utilized to estimate the phase offsets between the satellite clocks and ground clocks. This study, using the difference between the measured and predicted phase offset as a test statistic, sets a threshold for clock anomalies based on the prediction interval approach. Finally the validity of the monitoring method is examined by processing a set of real GNSS data that includes two recent incidents of clock anomalies in GNSS satellites.