Precise Point Positioning Using IGS Orbit and Clock Products

The contribution details a post-processing approach that used undifferentiated dual-frequency pseudorange and carrier phase observations along with IGS procise orbit products, for stand-alone precise geodetic point positioning (static or kinematic) with cm precision. This is possible if one takes advantage of the satellite clock estimates available with the satellite coordinates in the IGS precise orbit products and models systematic effects that cause cm variations in the satelite to user range. This paper will describe the approach, summarize the adjustment procedure, and specify the earth- and space-based models that must be implementetd to achieve cm-level positioning in static mode. Furthermore, station tropospheric zenth path delays with cm precision and GPS receiver clock estimates procise to 0.1 ns are also obtained. ? 2001 John Wiley & Sons, Inc.     

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.     

Impact of sampling rate of IGS satellite clock on precise point positioning

Both static and kinematic testings are investigated by using IGS 5min, 30s and 5s-interval precise satellite clock products in precise point positioning (PPP) solution. Test results show that the sampling rate of IGS satellite clock has very little effect on the static PPP solution. All the three types of sampling intervals of precise satellite clock can satisfy mm-cm level of positioning accuracy; higher sampling rate has no significant improvement for PPP solution. However, sampling rate of satellite clock has a significant impact on the PPP solution in kinematic PPP. The higher the interval of satellite clock, the better the accuracy achieved. The accuracy of kinematic PPP achieved by using 30s-interval precise satellite clock is improved by nearly 30–50 percent with respect to the solution by using 5min-interval precise satellite clock, but using 5s and 30s-interval satellite clock can almost produce the same accuracy of kinematic solution. Moreover, the use of precise satellite clock products from different analysis centers may also produce more or less effect on the PPP solution.     

精密卫星钟差加密方法及其对星载GPS低轨卫星定轨精度影响

系统分析、比较了几种精密卫星钟差加密方法，研究了利用全球分布的IGS永久跟踪站的GPS观测数据估计高采样率卫星钟差参数的原理与方法，并将各种卫星钟差加密方法得到的结果与IGS数据分析中心估计的卫星钟差结果相比较。最后将不同加密方法得出的精密卫星钟差结果用于基于星载GPS双频非差观测值的CHAMP低轨卫星的定轨，并将不同方法得到的定轨精度进行比较。结果表明，利用地面跟踪站的GPS观测数据，可高精度、高密度地估计GPS卫星钟差，估计精度可达0．1～0．5ns。经地面GPS跟踪站数据估计的GPS卫星钟差，应用于基于PPP方法的低轨卫星定轨，其定轨精度在10cm以内。     

Impact of sampling rate of IGS satellite clock on precise point positioning

Both static and kinematic testings are investigated by using IGS 5min, 30s and 5s-interval precise satellite clock products in precise point positioning (PPP) solution. Test results show that the sampling rate of IGS satellite clock has very little effect on the static PPP solution. All the three types of sampling intervals of precise satellite clock can satisfy mm-cm level of positioning accuracy; higher sampling rate has no significant improvement for PPP solution. However, sampling rate of satellite clock has a significant impact on the PPP solution in kinematic PPP. The higher the interval of satellite clock, the better the accuracy achieved. The accuracy of kinematic PPP achieved by using 30s-interval precise satellite clock is improved by nearly 30–50 percent with respect to the solution by using 5min-interval precise satellite clock, but using 5s and 30s-interval satellite clock can almost produce the same accuracy of kinematic solution. Moreover, the use of precise satellite clock products from different analysis centers may also produce more or less effect on the PPP solution.     

Real-time clock offset prediction with an improved model

The GPS orbit precision of the IGS ultra-rapid predicted (IGU-P) products has been remarkably improved since 2007. However, the satellite clock offsets of the IGU-P products have not shown sufficient high-quality prediction to achieve sub-decimeter precision in real-time precise point positioning (RTPPP), being at the level of 1–3 ns (30–90 cm) RMS in recent years. An improved prediction model for satellite clocks is proposed in order to enhance the precision of predicted clock offsets. First, the proposed prediction model adds a few cyclic terms to absorb the periodic effects, and a time adaptive function is used to adjust the weight of the observation in the prediction model. Second, initial deviations of the predictions are reduced by using a recomputed constant term. The simulation results have shown that the proposed prediction model can give a better performance than the IGU-P clock products and can achieve precision better than 0.55 ns (16.5 cm) in real-time predictions. In addition, the RTPPP method was chosen to test the efficiency of the new model for real-time static and kinematic positioning. The numerical examples using the data set of 140 IGS stations show that the static RTPPP precision based on the proposed clock model has been improved about 22.8 and 41.5 % in the east and height components compared to the IGU-P clock products, while the precisions in the north components are the equal. The kinematic example using three IGS stations shows that the kinematic RTPPP precision based on the proposed clock model has improved about 30, 72 and 44 % in the east, north and height components.     

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

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

顾及码频间偏差的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。     

IGS卫星钟差产品采样间隔对PPP精度的影响

使用IGS 5 min、30 s以及CODE最新发布的5 s间隔的精密卫星钟差产品分别进行了静态和动态精密单点定位(PPP)试验。结果表明,使用三种不同采样间隔的精密卫星钟差对静态PPP定位结果的影响很小,均能满足mm至cm级的静态定位精度,采样率更高的精密卫星钟差改正对静态定位结果无显著改善;对动态PPP定位,三种采样间隔的精密卫星钟差均能满足cm至dm级的定位精度,使用30 s间隔的精密卫星钟差较使用5 min间隔的精密卫星钟差,其定位精度提高了30%~50%,而使用5 s间隔的精密卫星钟差同使用30 s间隔的精密卫星钟差获得的定位精度基本一致。     

一种基于IGS超快星历的区域性实时精密单点定位方法

实时GPS精密单点定位需要实时的卫星轨道和钟差产品,为此提出一种利用区域GPS连续运行参考站和IGS发布的IGU超快轨道进行实时精密单点定位的方法.该方法首先利用连续运行参考站观测数据与IGU超快轨道预报部分进行实时GPS卫星钟差的估计,然后利用估计得到的实时GPS卫星钟差产品和IGU超快轨道预报部分,进行用户GPS接...     

几种精密卫星钟差加密方法的比较与分析

卫星钟差是影响GPS高精度单点定位的一个重要因素,通过分析几种加密GPS精密卫星钟差的方法,将加密结果与GFZ(German Research Centre for Geosciences)提供的数据进行比较,并通过精度分析,得出分段线性插值法是加密GPS精密卫星钟差一种较为可靠的方法。     

PPP与GAMIT/TRACK在地震监测中的应用

本文基于非差精密单点定位(PPP)和双差定位模式,采用研制的UNIP软件和GAMIT/TRACK运动学分析模块,对GPS技术在地震监测中的应用进行研究。利用IGS不同精度和时延的星历和钟差产品,分析了卫星轨道和钟差对两种软件动态定位结果的影响;采用不同的数据处理方案,分别提取GPS同震地表形变信息。结果表明:GAMIT/TRACK模块精密动态定位的实时性,更适用于地震等突发事件中的GPS快速响应;如果能够提供精确的卫星钟差,UNIP软件的高程精度明显优于GAMIT/TRACK模块,技术将更显优势。     

Improved GRACE kinematic orbit determination using GPS receiver clock modeling

Kinematic positions of Low Earth Orbiters based on GPS tracking are frequently used as pseudo-observations for single satellite gravity field determination. Unfortunately, the accuracy of the satellite trajectory is partly limited because the receiver synchronization error has to be estimated along with the kinematic coordinates at every observation epoch. We review the requirements for GPS receiver clock modeling in Precise Point Positioning (PPP) and analyze its impact on kinematic orbit determination for the two satellites of the Gravity Recovery and Climate Experiment (GRACE) mission using both simulated and real data. We demonstrate that a piecewise linear parameterization can be used to model the ultra-stable oscillators that drive the GPS receivers on board of the GRACE satellites. Using such a continuous clock model allows position estimation even if the number of usable GPS satellites drops to three and improves the robustness of the solution with respect to outliers. Furthermore, simulations indicate a potential accuracy improvement of the satellite trajectory of at least 40 % in the radial direction and up to 7 % in the along-track and cross-track directions when a 60-s piecewise linear clock model is estimated instead of epoch-wise independent receiver clock offsets. For PPP with real GRACE data, the accuracy evaluation is hampered by the lack of a reference orbit of significantly higher accuracy. However, comparisons with a smooth reduced-dynamic orbit indicate a significant reduction of the high-frequency noise in the radial component of the kinematic orbit.     

远海实时GPS潮汐的实时精密单点定位观测

针对传统事后精密单点定位技术的时间延迟问题,该文基于IGS RTS实时数据流产品,开展了实时精密单点定位技术在远海实时GPS验潮中的应用研究.对RTS改正的实时精密卫星轨道和钟差进行了精度验证和分析,给出了RT-PPP的数据处理策略以及实时GPS验潮的基本流程;组织和实施了渤海湾船载GPS验潮试验,以压力式验潮仪数据为参考,对远距离实时GPS潮汐测量结果进行了精度分析.结果表明:①以IGS最终卫星轨道和钟差产品为参考,RTS实时精密卫星轨道在X、y、Z方向的精度(RMS)均优于3 cm,卫星钟差的精度优于0.15 ns;②采用傅里叶低通滤波方法,消除波浪对潮汐观测的影响,进一步提取潮位信息.在忽略船体姿态改正的情况下,实时精密单点定位验潮相对于压力式验潮仪结果的最大偏差优于20 cm,RMS达到7.5 cm.     

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

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

实时精密单点定位的远海实时GPS潮汐观测

针对IGS实时数据流产品,该文开展了实时精密单点定位技术在远海实时GPS验潮中的应用研究。对RTS改正的实时精密卫星轨道和钟差进行了精度验证和分析,给出了实时精密单点定位的数据处理策略以及实时GPS验潮的基本流程;组织和实施了渤海湾船载GPS验潮试验,以压力式验潮仪数据为参考,对远距离实时GPS潮汐测量结果进行了精度分析。结果表明:以IGS最终卫星轨道和钟差产品为参考,实时数据流产品实时精密卫星轨道在X、Y、Z方向的精度均优于3cm,卫星钟差的精度优于0.15ns;采用傅里叶低通滤波方法,消除波浪对潮汐观测的影响,进一步提取潮位信息。在忽略船体姿态改正的情况下,实时精密单点定位验潮相对于压力式验潮仪结果的最大偏差优于20cm,RMS达到7.5cm。     

BDS/GPS精密单点定位收敛时间与定位精度的比较

采用武汉大学卫星导航定位技术研究中心发布的北斗精密卫星轨道和钟差,在TriP 2.0软件的基础上实现了BDS PPP定位算法,并利用大量实测数据进行了BDS/GPS静态PPP和动态PPP浮点解试验。结果表明,BDS静态PPP的收敛时间约为80min,动态PPP的收敛时间为100min;对于3h的观测数据,静态PPP收敛后定位精度优于5cm,动态PPP收敛后水平方向优于8cm,高程方向约12cm;与GPS PPP类似,东分量上定位精度较北分量稍差。当前由于BDS的全球跟踪站有限,精密轨道和钟差精度不如GPS,因此BDS PPP的收敛时间较GPS长,但收敛后可实现厘米至分米级的绝对定位。     

精密单点定位中卫星钟差加密方法的探讨与分析

本文系统地分析了多种精密卫星钟差加密方法,以IGS提供的GPSweek 1421周第2天的2号、11号、21号和28号卫星的30s间隔的精密卫星钟差为基准,然后以从中提取的15min和5min间隔的钟差为例,将其用内插或拟合等不同方法分别加密到30s,将加密结果与30s间隔的精密卫星钟差基准进行比较分析,得出样条函数内插法精度较高、较可靠,拟合法精度较差,在个别历元处,内插或拟合卫星钟差的误差仍然很大。     

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.     

Characterization of periodic variations in the GPS satellite clocks

The clock products of the International Global Navigation Satellite Systems (GNSS) Service (IGS) are used to characterize the timing performance of the GPS satellites. Using 5-min and 30-s observational samples and focusing only on the sub-daily regime, approximate power-law stochastic processes are found. The Block IIA Rb and Cs clocks obey predominantly random walk phase (or white frequency) noise processes. The Rb clocks are up to nearly an order of magnitude more stable and show a flicker phase noise component over intervals shorter than about 100 s. Due to the onboard Time Keeping System in the newer Block IIR and IIR-M satellites, their Rb clocks behave in a more complex way: as an apparent random walk phase process up to about 100 s and then changing to flicker phase up to a few thousand seconds. Superposed on this random background, periodic signals have been detected in all clock types at four harmonic frequencies, n × (2.0029 ± 0.0005) cycles per day (24 h coordinated universal time or UTC), for n = 1, 2, 3, and 4. The equivalent fundamental period is 11.9826 ± 0.0030 h, which surprisingly differs from the reported mean GPS orbital period of 11.9659 ± 0.0007 h by 60 ± 11 s. We cannot account for this apparent discrepancy but note that a clear relationship between the periodic signals and the orbital dynamics is evidenced for some satellites by modulations of the spectral amplitudes with eclipse season. All four harmonics are much smaller for the IIR and IIR-M satellites than for the older blocks. Awareness of the periodic variations can be used to improve the clock modeling, including for interpolation of tabulated IGS products for higher-rate GPS positioning and for predictions in real-time applications. This is especially true for high-accuracy uses, but could also benefit the standard GPS operational products. The observed stochastic properties of each satellite clock type are used to estimate the growth of interpolation and prediction errors with time interval.