A modified mix-differenced approach for estimating multi-GNSS real-time satellite clock offsets

A modified mixed-differenced approach for estimating multi-GNSS real-time clock offsets is presented. This approach, as compared to the earlier presented mixed-differenced approach which uses epoch-differenced and undifferenced observations, further adds a satellite-differenced process. The proposed approach, based on real-time orbit products and a mix of epoch-differenced and satellite-differenced observations to estimate only satellite clock offsets and tropospheric zenith wet delays, has fewer estimated parameters than other approaches, and thus its implementing procedure is efficient and can be performed and extended easily. To obtain high accuracy, the approach involves three steps. First, the high-accuracy tropospheric zenith wet delay of each station is estimated using mixed-differenced carrier phase observations. Second, satellite clock offset changes between adjacent epochs are estimated using also mixed-differenced carrier phase observations. Third, the satellite clock offsets at the initial epoch are estimated using satellite-differenced pseudorange observations. Finally, the initial epoch clock results and clock offset changes are concatenated to obtain the clock results of the current epoch. To validate the real-time satellite clock results, multi-GNSS post-processing clock products from IGS ACs were selected for comparison. From the comparison, the standard deviations of the GPS, GLONASS, BeiDou and Galileo systems clock results are approximately 0.1–0.4 ns, except for the BeiDou GEO satellites. The root mean squares are about 0.4–2.3 ns, which are similar to those of other international real-time products. When the clock estimates were assessed based on a pseudo-kinematic PPP procedure, the positioning accuracies in the East, North and Up components reach 5.6, 5.5 and 7.6 cm, respectively, which meet the centimeter level and are comparable to the application of other products.     

Satellite clock estimation at 1 Hz for realtime kinematic PPP applications

Realtime kinematic precise point positioning (PPP) requires 1 Hz GPS satellite clock corrections. An efficient clock estimation approach is presented. It applies a combined dual-thread algorithm consisting of an undifferenced (UD) and epoch-differenced (ED) engine. The UD engine produces absolute clock values every 5 s, and the ED engine produces relative clock values between neighboring epochs at 1-s interval. A final 1-Hz satellite clock can be generated by combining the UD absolute clock and ED relative clock efficiently and accurately. Forty stations from a global tracking network are used to estimate the realtime 1-Hz clock with the proposed method. Both the efficiency and accuracy of the resultant clock corrections are validated. Efficiency test shows that the UD processing thread requires an average time of 1.88 s on a 1-GHz CPU PC for one epoch of data, while ED processing requires only 0.25 s. Accuracy validation test shows that the estimated 1-Hz clock agrees with IGS final clock accurately. The RMS values of all the available GPS satellite clock bias are less than 0.2 ns (6 cm), and most of them are less than 0.1 ns (3 cm). All the RMS values of Signal in Space Range Error (SISRE) are at centimeter level. Applying the accurate and realtime clock to realtime PPP, an accuracy of 10 cm in the horizontal and 20 cm in the vertical is achieved after a short period of initialization.     

多模GNSS精密卫星钟差估计与分析

对基于历元间差分相位和非差伪距观测值的混合差分卫星钟差估计方法进行了改进,实现了多模全球导航卫星系统（Global Navigation Satellite System,GNSS）卫星钟差联合快速估计。选择了全球分布的50个跟踪站进行实验,对卫星钟差精度进行了分析和精密单点定位（Precise Point Positioning,PPP）验证。结果表明:多模卫星钟差与武汉大学提供的最终精密卫星钟差互差优于0.2 ns,精密单点定位结果与武汉大学发布的最终精密卫星轨道和钟差产品的定位精度相当。     

北斗/GPS实时精密卫星钟差融合解算模型及精度分析

实时GNSS精密单点定位(PPP)技术必须使用实时的高精度卫星精密轨道和钟差。本文研究了精密卫星钟差融合解算模型及策略,并利用滤波算法实现了北斗/GPS实时精密卫星钟差融合估计算法。仿真实时试验结果显示:获得的北斗/GPS实时钟差与GFZ事后多GNSS精密钟差(GBM)的标准差在0.15 ns左右;使用该钟差进行GPS动态PPP试验,收敛后水平精度优于5 cm,高程精度优于10 cm;使用仿真实时钟差进行的北斗动态PPP与使用GFZ事后多GNSS精密钟差开展的试验相比精度相当,可实现分米级定位。     

卫星钟差解算及其星间单差模糊度固定

整数相位模糊度解算可以显著提高GNSS精密单点定位（PPP）的精度。本文提出一种解算卫星钟差的方法,通过固定星间单差模糊度恢复出能够支持单台接收机进行整数模糊度解算的卫星钟差,即所谓的“整数”钟差。为了实现星间单差模糊度固定,分别通过卫星端宽巷FCB解算和模糊度基准的选择与固定恢复出宽巷和窄巷模糊度的整数性质。为了证明本文方法的可行性,采用IGS测站的GPS数据进行卫星钟差解算试验。结果表明,在解算钟差时,星间单差模糊度固定的平均成功率为73%。得到的卫星钟差与IGS最终钟差产品相比,平均的RMS和STD分别为0.170和0.012 ns。448个IGS测站的星间单差宽巷和窄巷模糊度小数部分的分布表明本文得到的卫星钟差和FCB产品具备支持PPP用户进行模糊度固定的能力。基于以上产品开展了模拟动态PPP定位试验,结果表明模糊度固定之后,N、E、U和3D的定位精度（RMS）分别达到0.009、0.010、0.023和0.027 m,与不固定模糊度或采用IGS钟差的结果相比,分别提高了30.8%、61.5%、23.3%和37.2%。     

GBM快速轨道产品及非差模糊度固定对其精度的改进

GNSS是实时定位导航最重要的方法,精密卫星轨道钟差产品是GNSS高精度服务的前提。国际GNSS服务中心(IGS)及其分析中心长期致力于GNSS数据处理的研究及高精度轨道和钟差产品的提供。GFZ作为分析中心之一,提供GBM多系统快速产品。本文基于2015—2021年GBM提供的精密轨道产品,阐述了数据处理策略,分析了轨道的精度,介绍了非差模糊度固定的原理和对精密定轨的影响。结果表明:GBM快速产品中的GPS轨道精度与IGS后处理精密轨道相比的精度约为11~13 mm,轨道6 h预报精度约为6 cm;GLONASS预报精度约为12 cm,Galileo在该时期的精度均值为10 cm,但是在2016年底以后精度提升到5 cm左右;北斗系统的中轨卫星(medium earth orbit,MEO)在2020年以后预报精度约为10 cm;北斗的静止轨道卫星(geostationary earth orbit,GEO)卫星和QZSS卫星的预报精度在米级;卫星激光测距检核表明,Galileo、GLONASS、BDS-3 MEO卫星轨道精度分别为23、41、47 mm;此外,采用150 d观测值的试验结果表明,采用非差模糊度固定能显著改善MEO卫星轨道精度,对GPS、GLONASS、Galileo、BDS-2和BDS-3的MEO卫星的6 h时预报精度改善率分别为9%~15%、15%~18%、11%~13%、6%~17%和14%~25%。     

Single epoch estimation of the Galileo integrity chain sensor station clock offsets

The Galileo integrity chain depends on a number of key factors, one of which is contamination of the signal-in-space errors with residual errors other than imperfect modelling of satellite orbits and clocks. A potential consequence of this is that the user protection limit is driven not by the errors associated with the imperfect orbit and clock modelling, but by the distortions induced by noise and bias in the integrity chain. These distortions increase the minimum bias the integrity chain can guarantee to detect, which is reflected in the user protection limit. A contributor to this distortion is the inaccuracy associated with the estimation of the offset between the Galileo sensor station (GSS) receiver clocks and the Galileo system time (GST). This offset is termed the receiver clock synchronization error (CSE). This paper describes the research carried out to determine both the CSE and its associated error using GPS data as captured with the Galileo System Test Bed Version 1 (GSTB-V1). In the study we simulate open access to a time datum using IGS data. Two methods are compared for determining CSE and the corresponding uncertainty (noise) across a global network of tracking stations. The single-epoch single-station method is an ‘averaging’ technique that uses a single epoch of data, and is carried out at individual sensor stations, without recourse to the data from other stations. The global network solution method is also single epoch based, but uses the inversion of a linearised model of the global system to solve for the CSE simultaneously at all GSS along with a number of other parameters that would otherwise be absorbed into the CSE estimate in the averaging technique. To test the effectiveness of various configurations in the two methods the estimated synchronisation errors across the GSS network (comprising 25 stations) are compared to the same values as estimated by the International GPS Service (IGS) using a global tracking network of around 150 stations, as well as precise orbit and satellite clock models determined by a combination of global analysis centres. The results show that the averaging technique is vulnerable to unmodelled errors in the satellite clock offsets from system time, leading to receiver CSE errors in the region of 12 ns (3.7 m), this value being largely driven by the satellite CSE errors. The global network approach is capable of delivering CSE errors at the level of 1.5 ns (46 cm) depending on the number of parameters in the linearised model. The International GNSS Service (IGS) receiver clock estimates were used as a truth model for comparative assessment.     

High-rate GPS clock corrections from CODE: support of 1 Hz applications

GPS zero-difference applications with a sampling rate up to 1 Hz require corresponding high-rate GPS clock corrections. The determination of the clock corrections in a full network solution is a time-consuming task. The Center for Orbit Determination in Europe (CODE) has developed an efficient algorithm based on epoch-differenced phase observations, which allows to generate high-rate clock corrections within reasonably short time (< 2 h) and with sufficient accuracy (on the same level as the CODE rapid or final clock corrections, respectively). The clock determination procedure at CODE and the new algorithm is described in detail. It is shown that the simplifications to speed up the processing are not causing a significant loss of accuracy for the clock corrections. The high-rate clock corrections have in essence the same quality as clock corrections determined in a full network solution. In order to support 1 Hz applications 1-s clock corrections would be needed. The computation time, even for the efficient algorithm, is not negligible, however. Therefore, we studied whether a reduced sampling is sufficient for the GPS satellite clock corrections to reach the same or only slightly inferior level of accuracy as for the full 1-s clock correction set. We show that high-rate satellite clock corrections with a spacing of 5 s may be linearly interpolated resulting in less than 2% degradation of accuracy.     

顾及卫星姿态的多系统精密钟差产品综合

国际GNSS服务(IGS)提供的GPS综合产品被广泛应用于各种高精度科学研究中. 随着各国卫星导航系统的发展,亟需研究针对多系统全球卫星导航系统(GNSS)产品的综合策略. 由于卫星姿态与钟差相互耦合,综合钟差时额外考虑姿态改正将进一步提高综合产品精度,因此研究了一种顾及卫星姿态的GNSS钟差综合策略,改正姿态后GPS综合残差最大可减小80%. 对142个IGS测站进行精密单点定位(PPP)解算发现,综合产品比单个分析中心产品更加稳定,东(E)、北(N)、高(U)方向的动态定位精度最大可提升22.7%、16.7%和18.3%. 相对于未顾及姿态改正的综合产品,顾及姿态改正的综合产品的动态定位精度最大可提升65.3%.      

Fast estimation and analysis of the inter-frequency clock bias for Block IIF satellites

The inter-frequency bias of PRN25 was noticed by the scientific community and considered to be caused by thermal variations. The inter-frequency bias leads to an apparent inter-frequency clock bias (IFCB), which could be obtained using the difference of two ionosphere-free phase combinations (L1/L2 and L1/L5). We present an efficient approach derived from the epoch-differenced strategy for fast estimation of IFCBs for Block IIF satellites. For the analysis, data from 32 stations from the IGS network spanning 10 months (DOY 213, 2011–153, 2012) are processed. The processing times show that the epoch-differenced method is more efficient than the undifferenced one. In order to study the features of IFCB, a harmonic analysis is performed by using a FFT (fast Fourier transformation), and significant periodic variations with the periods of 12, 6 and 8 h are noticed. The fourth-order period is determined by comparing the performances of the model with different periods. After determination, a harmonics-based function of order 4 is used to model the IFCB, and the single-day amplitudes and phases are estimated for the 10 months from a least squares fit. Based on the estimated results, the characterization of IFCB is discussed. The algorithm is incorporated into the MGPSS software developed at SHAO (Shanghai Astronomical Observatory, Chinese Academy of Sciences) and used to monitor the IFCB variations of GPS and COMPASS systems in near real time.     

顾及码频间偏差的GPS/GLONASS实时卫星钟差估计

在进行GPS/GLONASS联合卫星钟差估计时,GLONASS码频间偏差（inter-frequency bias,IFB）因卫星频率间的差异而无法被测站接收机钟差参数吸收,其一部分将进入GLONASS卫星钟差估值中。通过引入多个"时频偏差"参数（inter-system and inter-frequency bias,ISFB）及附加基准约束对测站GLONASS码IFB进行函数模型补偿,实现其与待估卫星钟差参数的有效分离,并对所估计实时卫星钟差和实时精度单点定位（real-time precise point positioning,RT-PPP）进行精度评估。结果表明,在卫星钟差估计观测方程中忽略码IFB,会明显降低GLONASS卫星钟差估值精度;新方法能有效避免码IFB对卫星钟差估值的影响,所获得GPS、GLONASS卫星钟差与ESA（European Space Agency）事后精密钟差产品偏差平均均方根值分别小于0.2 ns、0.3 ns。利用实时估计卫星钟差进行静态RT-PPP,当观测时段长为2 h时,GPS单系统、GPS/GLONASS组合系统的3D定位精度优于10 cm,GLONASS单系统3D定位精度约为15 cm;三种模式24 h单天解的3D定位精度均优于5 cm。     

An enhanced strategy for GNSS data processing of massive networks

Although the computational burden of global navigation satellite systems (GNSS) data processing is nowadays already a big challenge, especially for huge networks, integrated processing of denser networks with data of multi-GNSS and multi-frequency is desired in the expectation of more accurate and reliable products. Based on the concept of carrier range, in this study, the precise point positioning with integer ambiguity resolution is engaged to obtain the integer ambiguities for converting carrier phases to carrier ranges. With such carrier ranges and pseudo-ranges, rigorous integrated processing is realized computational efficiently for the orbit and clock estimation using massive networks. The strategy is validated in terms of computational efficiency and product quality using data of the IGS network with about 460 stations. The experimental validation shows that the computation time of the new strategy increases gradually with the number of stations. It takes about 14 min for precise orbit and clock determination with 460 stations, while the current strategy needs about 82 min. The overlapping orbit RMS is reduced from 27.6 mm with 100 stations to 24.8 mm using the proposed strategy, and the RMS could be further reduced to 23.2 mm by including all 460 stations. Therefore, the new strategy could be applied to massive networks of multi-GNSS and multi-frequency receivers and possibly to achieve GNSS data products of higher quality.     

Improved antenna phase center models for GLONASS

Thanks to the increasing number of active GLONASS satellites and the increasing number of multi-GNSS tracking stations in the network of the International GNSS Service (IGS), the quality of the GLONASS orbits has become significantly better over the last few years. By the end of 2008, the orbit RMS error had reached a level of 3–4 cm. Nevertheless, the strategy to process GLONASS observations still has deficiencies: one simplification, as applied within the IGS today, is the use of phase center models for receiver antennas for the GLONASS observations, which were derived from GPS measurements only, by ignoring the different frequency range. Geo++ GmbH calibrates GNSS receiver antennas using a robot in the field. This procedure yields now separate corrections for the receiver antenna phase centers for each navigation satellite system, provided its constellation is sufficiently populated. With a limited set of GLONASS calibrations, it is possible to assess the impact of GNSS-specific receiver antenna corrections that are ignored within the IGS so far. The antenna phase center model for the GLONASS satellites was derived in early 2006, when the multi-GNSS tracking network of the IGS was much sparser than it is today. Furthermore, many satellites of the constellation at that time have in the meantime been replaced by the latest generation of GLONASS-M satellites. For that reason, this paper also provides an update and extension of the presently used correction tables for the GLONASS satellite antenna phase centers for the current constellation of GLONASS satellites. The updated GLONASS antenna phase center model helps to improve the orbit quality.     

单频到五频多系统GNSS精密单点定位参数估计与应用

全球导航卫星系统(GNSS)已发展至多频多系统时代,特别以我国北斗卫星导航系统(BDS)为代表的四大全球导航卫星系统可全天时、全天候播发十几个频率的伪距、相位和多普勒等观测信息。多频多系统GNSS为用户提供更多的观测数据和组合选择,为精密定位、导航和授时(PNT)应用带来了新的机遇,如高精度位置服务、大地测量、空间天气和灾害监测等。但多频多系统GNSS观测为精密单点定位(PPP)组合模型和系统偏差及大气延迟估计等带来诸多问题和挑战。本文给出了单频到五频多系统GNSS精密单点定位(PPP)模型,估计和评估了单频到五频多系统GNSS PPP定位精度、接收机钟差、对流层延迟、卫星和接收机硬件延迟,以及频间偏差。给出了GNSS PPP最新应用进展,包括GNSS气象学、电离层模拟、时间频率传递、建筑物安全和地震监测及其应用。结果表明,多频多系统极大地提高了GNSS PPP参数估计的精度和可靠性,具有重要的应用价值。最后给出了多频多系统GNSS PPP应用前景与展望。     

GPS satellite clock determination in case of inter-frequency clock biases for triple-frequency precise point positioning

Significant time-varying inter-frequency clock biases (IFCBs) within GPS observations prevent the application of the legacy L1/L2 ionosphere-free clock products on L5 signals. Conventional approaches overcoming this problem are to estimate L1/L5 ionosphere-free clocks in addition to their L1/L2 counterparts or to compute IFCBs between the L1/L2 and L1/L5 clocks which are later modeled through a harmonic analysis. In contrast, we start from the undifferenced uncombined GNSS model and propose an alternative approach where a second satellite clock parameter dedicated to the L5 signals is estimated along with the legacy L1/L2 clock. In this manner, we do not need to rely on the correlated L1/L2 and L1/L5 ionosphere-free observables which complicates triple-frequency GPS stochastic models, or account for the unfavorable time-varying hardware biases in undifferenced GPS functional models since they can be absorbed by the L5 clocks. An extra advantage over the ionosphere-free model is that external ionosphere constraints can potentially be introduced to improve PPP. With 27 days of triple-frequency GPS data from globally distributed stations, we find that the RMS of the positioning differences between our GPS model and all conventional models is below 1 mm for all east, north and up components, demonstrating the effectiveness of our model in addressing triple-frequency observations and time-varying IFCBs. Moreover, we can combine the L1/L2 and L5 clocks derived from our model to calculate precisely the L1/L5 clocks which in practice only depart from their legacy counterparts by less than 0.006 ns in RMS. Our triple-frequency GPS model proves convenient and efficient in combating time-varying IFCBs and can be generalized to more than three frequency signals for satellite clock determination.     

Kalman-filter-based GPS clock estimation for near real-time positioning

In this article, an algorithm for clock offset estimation of the GPS satellites is presented. The algorithm is based on a Kalman-filter and processes undifferenced code and carrier-phase measurements of a global tracking network. The clock offset and drift of the satellite clocks are estimated along with tracking station clock offsets, tropospheric zenith path delay and carrier-phase ambiguities. The article provides a brief overview of already existing near-real-time and real-time clock products. The filter algorithm and data processing scheme is presented. Finally, the accuracy of the orbit and clock product is assessed with a precise orbit determination of the MetOp satellite and compared to results gained with other real-time products.

GNSS增强系统中精密实时钟差高频估计及应用研究

GNSS星基差分增强系统依赖于实时轨道及钟差增强信息。本文主要研究多GNSS实时精密钟差估计模型,在传统非差基础上优化待估参数,实现了一种高效的Multi-GNSS实时钟差简化估计模型。基于PANDA软件开展了实时轨道数据处理与分析,经过验证可获得的GPS/北斗MEO/Galileo实时轨道径向精度1~5cm,北斗GEO/IGSO卫星径向精度约10cm。分析发现本文优化的实时钟差简化估计模型单历元解算效率较高,可应用于实时钟差增强信息高频(如1Hz)更新,且解算获得的实时钟差不存在常偏为绝对钟差;基于实时轨道,通过该模型可获得实时钟差精度GPS约0.22ns,北斗GEO约0.50ns、IGSO/MEO约0.24ns,Galileo约0.32ns。在此基础上,利用目前所获取的MultiGNSS实时数据流搭建了Multi-GNSS全球实时增强原型系统,并基于互联网实时播发增强信息,可初步实现实时PPP厘米级服务、伪距米级导航定位服务。     

A new differential code bias (C1–P1) estimation method and its performance evaluation

The current satellite clock products are computed using the ionosphere-free phase (L1/L2) and code (P1/P2) observations. Thus, if users conduct undifferenced positioning using these clock products together with C1 and P2 observations, the differential code bias (DCB) (C1–P1) should be properly compensated. The influence of DCB (C1–P1) on the undifferenced ambiguity solutions is investigated. Based on the investigation, we propose a new DCB (C1–P1) estimation method. Using it, the satellite DCB (C1–P1) can be computed. A 30-day (DOY 205–234, 2012) dual-frequency GPS data set is processed to estimate the DCB (C1–P1). Comparing the estimated results with that of IGS DCB products, the accuracy is better than 0.13 m. The performances of DCB (C1–P1) in the code-based single-point positioning, precise point positioning (PPP) convergence and wide-lane uncalibrated phase delay (UPD) estimation are investigated using the estimated DCB (C1–P1). The results of the code-based single-point positioning show that the influence of DCB (C1–P1) on the up direction is more evident than on the horizontal directions. The accuracy is improved by 50 % and reaches to decimeter level with DCB (C1–P1) application. The performance of DCB (C1–P1) in PPP shows that it can accelerate PPP convergence through improving the accuracy of the code observation. The computed UPD values show that influence of DCB (C1–P1) on UPD of each satellite is different, and some values are larger than 0.3 cycles.     

GPS非差相位精密单点定位方法与实现

探讨了精密单点定位的基本原理、处理方法、所涉及的误差改正及数据处理中的一些关键技术;采用直接内插IGS卫星精密星历的方法代替利用IGS跟踪站进行轨道精化方法计算卫星轨道参数,对现有精密单点定位计算方法进行了简化,使之更具有实用性。最后利用自主研发的精密非差单点定位软件计算和分析了实测数据。计算结果表明,经过大约15 min的初始化后,非差相位单历元的定位结果精确度在X,Y,Z方向上均优于20 cm。     

顾及钟差变化率的GPS卫星钟差预报法

提出了一种基于钟差变化率拟合建模的卫星钟差预报方法。以附加周期项的线性或二次多项式作为基础模型对钟差变化率序列进行拟合,最优估计卫星钟差的趋势项系数,然后直接使用精密定轨得到的相应时刻的卫星钟差计算预报初始时刻的基准项系数,来建立卫星钟差的预报模型。以IGS发布的快速星历(IGR)的卫星钟差为试验数据,对GPS星座中各种型号的所有卫星钟差进行预报。结果表明:本文方法3、6、12与24h的预报精度分别可达0.43、0.58、0.90与1.47ns,相比于传统的基于钟差拟合的预报方法,精度分别提高69.3%、61.8%、50.5%与37.2%;与IGS发布的超快速星历(IGU)的预报钟差相比,钟差精度分别提高15.7%、23.7%、27.4%与34.4%。