Ambiguity resolution in precise point positioning with hourly data for global single receiver

Integer ambiguity resolution (IAR) can improve precise point positioning (PPP) performance significantly. IAR for PPP became a highlight topic in global positioning system (GPS) community in recent years. More and more researchers focus on this issue. Progress has been made in the latest years. In this paper, we aim at investigating and demonstrating the performance of a global zero-differenced (ZD) PPP IAR service for GPS users by providing routine ZD uncalibrated fractional offsets (UFOs) for wide-lane and narrow-lane. Data sets from all IGS stations collected on DOY 1, 100, 200 and 300 of 2010 are used to validate and demonstrate this global service. Static experiment results show that an accuracy better than 1 cm in horizontal and 1–2 cm in vertical could be achieved in ambiguity-fixed PPP solution with only hourly data. Compared with PPP float solution, an average improvement reaches 58.2% in east, 28.3% in north and 23.8% in vertical for all tested stations. Results of kinematic experiments show that the RMS of kinematic PPP solutions can be improved from 21.6, 16.6 and 37.7 mm to 12.2, 13.3 and 34.3 mm for the fixed solutions in the east, north and vertical components, respectively. Both static and kinematic experiments show that wide-lane and narrow-lane UFO products of all satellites can be generated and provided in a routine way accompanying satellite orbit and clock products for the PPP user anywhere around the world, to obtain accurate and reliable ambiguity-fixed PPP solutions.     

应用iGMAS超快速星历的实时精密单点定位研究

中国主导建设的国际GNSS监测评估系统（iGMAS）相比国际上比较成熟的IGS系统在产品精度等方面存在差别，目前实时精密单点定位应用多采用IGS实时、近实时产品。为改变这一现状，针对iGMAS产品特性以及实时精密单点定位对超快速精密星历的需求，对iGMAS超快速星历的精度和稳定性方面进行评估，设计了iGMAS产品实时/事后下载应用程序，开展了基于iGMAS超快速星历的实时精密单点定位研究，并结合NovAtel OEM617双频接收机进行了GPS实时精密单点定位试验。实验结果表明，在连续观测23min后定位误差即可收敛到分米级，较接收机原始定位精度高一个量级，且稳定性好，最终在E/N/U方向定位误差均方根分别为7.2cm、6.4cm、15.2cm，与应用IGS超快速星历实时PPP试验取得相近的结果。研究实现了iGMAS数据获取、评估和实时PPP应用的一整套方案，验证了iGMAS超快速产品的性能，对推进iGMAS产品的应用提供了借鉴。     

Evaluation of COMPASS ionospheric model in GNSS positioning

As important products of GNSS navigation message, ionospheric delay model parameters are broadcasted for single-frequency users to improve their positioning accuracy. GPS provides daily Klobuchar ionospheric model parameters based on geomagnetic reference frame, while the regional satellite navigation system of China’s COMPASS broadcasts an eight-parameter ionospheric model, COMPASS Ionospheric Model(CIM), which was generated by processing data from continuous monitoring stations, with updating the parameters every 2 h. To evaluate its performance, CIM predictions are compared to ionospheric delay measurements, along with GPS positioning accuracy comparisons. Real observed data analysis indicates that CIM provides higher correction precision in middle-latitude regions, but relatively lower correction precision for low-latitude regions where the ionosphere has much higher variability. CIM errors for some users show a common bias for in-coming COMPASS signals from different satellites, and hence ionospheric model errors are somehow translated into the receivers’ clock error estimation. In addition, the CIM from the China regional monitoring network are further evaluated for global ionospheric corrections. Results show that in the Northern Hemisphere areas including Asia, Europe and North America, the three-dimensional positioning accuracy using the CIM for ionospheric delay corrections is improved by 7.8%–35.3% when compared to GPS single-frequency positioning ionospheric delay corrections using the Klobuchar model. However, the positioning accuracy in the Southern Hemisphere is degraded due apparently to the lack of monitoring stations there.     

Mitigation of receiver biases in ionospheric observables from PPP with ambiguity resolution

PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) positioning method that requires SSR (State Space Representation) corrections in order to provide solutions with an accuracy of centimetric level. The so-called RT-PPP (Real-time PPP) is possible thanks to real-time precise SSR products, for orbits and clocks, provided by IGS (International GNSS Service) and its associate analysis centers such as CNES (Centre National d'Etudes Spatiales). CNES SSR products also enable RT-PPP with integer ambiguity resolution. In GNSS related literature, PPP with ambiguity resolution (PPP-AR) in real-time is often referred as PPP-RTK (PPP – Real Time Kinematic). PPP-WIZARD (PPP - With Integer and Zero-difference Ambiguity Resolution Demonstrator) is a software that is made available by CNES. This software is capable of performing PPP-RTK. It estimates slant ionospheric delays and other GNSS positioning parameters. Since ionospheric effects are spatially correlated by GNSS data from active networks, it is possible to model and provide ionospheric delays for any position in the network coverage area. The prior knowledge ionospheric delays can reduce positioning convergence for PPP-RTK users. Real-time ionospheric models could benefit from highly precise ionospheric delays estimated in PPP-AR. In this study, we demonstrate that ionospheric delays obtained throughout PPP-AR estimation are actu ally ionospheric observables. Ionospheric observables are biased by an order of few meters caused by the receiver hardware biases. These biases prohibit the use of PPP-WIZARD ionospheric delays to produce ionospheric models. Receiver biases correction is essential to provide ionospheric delays while using PPP-AR based ionospheric observables. In this contribution, a method was implemented to estimate and mitigate receiver hardware biases influence on slant ionospheric observables from PPP-AR. In order to assess the proposed approach, PPP-AR data from 12 GNSS stations were processed over a two-month period (March and April 2018). A comparison between IGS ionospheric products and PPP-AR based ionospheric observables corrected for receiver biases, resulted in a mean of differences of −39 cm and 51 cm standard deviation. The results are consistent with the accuracy of the IGS ionospheric products, 2–8 TECU, considering that 1 TECU is ~16 cm in L1. In another analysis, a comparison of ionospheric delays from 5 pairs of short baselines GNSS stations found an agreement of 0.001 m in mean differences with 22 cm standard deviation after receiver biases were corrected. Therefore, the proposed solution is promising and could produce high quality (1–2 TECU) slant ionospheric delays. This product can be used in a large variety of modeling approaches, since ionospheric delays after correction are unbiased. These results indicate that the proposed strategy is promising, and could benefit applications that require accuracy of 1–2 TECU (~16–32 cm in L1).     

Towards PPP-RTK: Ambiguity resolution in real-time precise point positioning

Integer ambiguity resolution at a single station can be achieved by introducing predetermined uncalibrated phase delays (UPDs) into the float ambiguity estimates of precise point positioning (PPP). This integer resolution technique has the potential of leading to a PPP-RTK (real-time kinematic) model where PPP provides rapid convergence to a reliable centimeter-level positioning accuracy based on an RTK reference network. Nonetheless, implementing this model is technically subject to how rapidly we can fix wide-lane ambiguities, stabilize narrow-lane UPD estimates, and achieve the first ambiguity-fixed solution. To investigate these issues, we used 7 days of 1-Hz sampling GPS data at 91 stations across Europe. We find that at least 10 min of observations are required for most receiver types to reliably fix about 90% of wide-lane ambiguities corresponding to high elevations, and over 20 min to fix about 90% of those corresponding to low elevations. Moreover, several tens of minutes are usually required for a regional network before a narrow-lane UPD estimate stabilizes to an accuracy of far better than 0.1 cycles. Finally, for hourly data, ambiguity resolution can significantly improve the accuracy of epoch-wise position estimates from 13.7, 7.1 and 11.4 cm to 0.8, 0.9 and 2.5 cm for the East, North and Up components, respectively, but a few tens of minutes is required to achieve the first ambiguity-fixed solution. Therefore, from the timeliness aspect, our PPP-RTK model currently cannot satisfy the critical requirement of instantaneous precise positioning where ambiguity-fixed solutions have to be achieved within at most a few seconds. However, this model can still be potentially applied to some near-real-time remote sensing applications, such as the GPS meteorology.     

基于相位偏差的实时PPP模糊度固定

采用CNES发布的实时相位偏差数据，实现包含模糊度固定的实时精密单点定位.对全球10个IGS测站10天观测数据进行RTPPP解算，分别统计模糊度首次固定时间和定位精度，结果显示利用实时相位偏差数据能在平均30min内实现模糊度首次固定，模糊度固定时水平位置误差由6cm迅速降低至2cm左右，三维位置误差由10cm迅速降低至5cm左右，同时RTPPP模糊度固定在3h观测内可保持水平3cm、三维5cm左右的定位精度.通过分析得出，基于相位偏差的RTPPP模糊度固定技术具有较高的定位精度和定位稳定性，能够快速实现cm级定位.      

基于固定模糊度PPP的低轨卫星实时精密定轨

传统动力学定轨法受制于动力学模型精度，传统几何定轨法精度受限，只能达到亚米级，而基于精密单点定位（PPP）模式的几何定轨法一般采用浮点解，定轨精度及可靠性较基于双差模式的相对定位较差。为提高PPP模式低轨定轨的定位性能，利用中国区域内外的IGS测站计算出当前所有卫星的宽巷和窄巷相位小数偏差产品，对经过中国大陆区域上空的国产低轨卫星海洋二号（HY-2）和资源三号 (ZY-3) 卫星进行固定模糊度PPP的定轨解算，与事后精密轨道结果进行比较，分析其外符合精度。结果表明：仅利用约10min弧段的HY-2和ZY-3卫星数据，切向与径向的定轨精度可达2cm左右，法向为5cm左右，较浮点解定轨精度大幅提升。基于固定模糊度PPP的定轨方法能够满足厘米级的实时精密定轨。     

Performance assessment of real-time multi-GNSS integrated PPP with uncombined and ionospheric-free combined observables

For precise position services, the real-time precise point positioning (PPP) is a promising technology. The real-time PPP performance is expected to be improved by multi-system combination. The performance of real-time multi-system PPP needs to be periodically investigated, with the increasing number of available satellites and the continuously improved quality of real-time precise products of satellite clocks and orbits. In this study, a comprehensive performance assessment is conducted for the four-system integrated real-time PPP (FSIRT-PPP) with GPS, BDS, Galileo and GLONASS in both static and kinematic modes. The datasets from 118 stations spanning approximately a month are used for analysis, and the real-time stream CLK93 is employed. The superior performance of FSIRT-PPP is validated by comparing with the results of GPS/BDS, GPS/Galileo, GPS/GLONASS, GPS-only, BDS-only, Galileo-only and GLONASS-only cases. The FSIRT-PPP using ionospheric-free (IF) combined observables can achieve a convergence time of 10.9, 4.8 and 11.8 min and a positioning accuracy of 0.4, 0.5 and 0.7 cm in the static mode in the east, north and up directions, respectively, while the derived statistic is 15.4, 7.0 and 16.4 min, and 1.6, 1.2 and 3.4 cm in the kinematic mode in the three directions, respectively. Moreover, we also compare the position solutions of real-time PPP adopting IF combined and uncombined (UC) observables, and prove the mathematical equivalence between the two PPP models in the converged stage, provided that there are no external ionospheric corrections or constraints given to the estimated ionospheric delays in the UC model. The difference between the fully converged positioning accuracy of IF-based and UC-based real-time PPP is marginal, but the UC-based real-time PPP has longer convergence time due to the influence of the significant unmodeled time-varying errors in the real-time precise products as well as the different parameterization between them. For completeness, the real-time kinematic PPP results in harsh environments and the post-processed PPP results are also presented.     

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.     

Predicting atmospheric delays for rapid ambiguity resolution in precise point positioning

Integer ambiguity resolution in precise point positioning (PPP) can shorten the initialization and re-initialization time, and ambiguity-fixed PPP solutions are also more reliable and accurate than ambiguity-float PPP solutions. However, signal interruptions are unavoidable in practical applications, particularly while operating in urban areas. Such signal interruptions can cause discontinuity of carrier phase arc, which introduces new integer ambiguities. Usually it will take approximately 15 min of continuous tracking to a reasonable number of satellites to fix new integer ambiguities. In many applications, it is impractical for a PPP user to wait for such a long time for the re-initialization. In this paper, a method for rapid ambiguity fixing in PPP is developed to avoid such a long re-initialization time. Firstly, the atmospheric delays were estimated epoch by epoch from ambiguity-fixed PPP solutions before the data gap or cycle slip occurs. A random walk procedure is then applied to predict the atmospheric delays accurately over a short time span. The predicted atmospheric delays then can be used to correct the observations which suffer from signal interruptions. Finally, the new ambiguities can be fixed with a distinct WL-LX-L3 (here LX denotes either of L1, L2) cascade ambiguity resolution strategy. Comprehensive experiments have demonstrated that the proposed method and strategy can fix zero-difference integer ambiguities successfully with only a single-epoch observation immediately after a short data gap. This technique works even when all satellites are interrupted at the same time. The duration of data gap bridged by this technique could be possibly extended if a more precise atmospheric delay prediction is found or on-the-fly (OTF) technology is applied. Based on the proposed method, real-time PPP with integer ambiguity fixing becomes more feasible in practice.