Precise coseismic displacements from the GPS variometric approach using different precise products: Application to the 2008 MW 7.9 Wenchuan earthquake

The Global Positioning System (GPS) variometric approach has emerged as an attractive alternative to traditional well-developed positioning techniques including relative positioning and precise point positioning. Previous studies have demonstrated the capability of the variometric approach to retrieve coseismic displacements at centimeter-level precision, in a real-time manner using only readily available broadcast ephemeris. This study presents the first results comparing the performance of the variometric approach by using a variety of precise satellite orbit and clock products. Totally six kinds of products are included in our evaluation, namely the broadcast, IGS (International GNSS Service) ultra-rapid (predicted), ultra-rapid (observed), rapid, final (30-s clock) and CODE (Center for Orbit Determination in Europe) final (5-s clock) products. Static and dynamic experiments are conducted using 1-Hz GPS data covering a relatively large area in China during the 2008 Wenchuan M_W 7.9 earthquake. After removing the linear trend, the displacements using broadcast, ultra-rapid (predicted), ultra-rapid (observed) and rapid products reach nearly equivalent precisions at centimeter level. By using final and CODE final products, the precision of displacements can be significantly improved from 1.9–2.0 cm to 0.4–0.7 cm horizontally, and from 6.0–6.2 cm to 1.0–1.7 cm vertically for the dynamic experiments. The displacements using the CODE final products achieve the best precision, improved by more than 40% compared to those using the IGS final products. With the availability of IGS high-rate real-time precise products, this approach is promising to capture coseismic displacements more precisely in real time, which is crucial for earthquake and tsunami early warning.     

Performance assessment of RTPPP positioning with SSR corrections and PPP-AR positioning with FCB for multi-GNSS from MADOCA products

The state-space representation (SSR) product of satellite orbit and clock is one of the most essential corrections for real-time precise point positioning (RTPPP). When it comes to PPP ambiguity resolution (PPP-AR), the fractional cycle bias (FCB) matters. The Japan Aerospace Exploration Agency (JAXA) has developed a multi-GNSS (i.e., global navigation satellite system) advanced demonstration tool for orbit and clock analysis (MADOCA), providing free and precise orbit and clock products. Because of the shortage of relevant studies on performance evaluation, this paper focuses on the performance assessment of RTPPP and PPP-AR by real-time and offline MADOCA products. To begin with, the real-time MADOCA products are evaluated by comparing orbit and clock with JAXA final products, which gives an objective impression of the correction. Second, PPP tests in static and simulated kinematic mode are conducted to further verify the quality of real-time MADOCA products. Finally, the offline MADOCA products are assessed by PPP and PPP-AR comparisons. The results are as follows: (1) Orbit comparisons produced an average error of about 0.04–0.13 m for the global positioning system (GPS), 0.14–0.16 m for the global navigation satellite system (GLONASS), and 0.07–0.08 m for the quasi-zenith satellite system (QZSS). The G15 satellite had the most accurate orbit, with a difference of 0.04 m between the JAXA orbit products and MADOCA’s counterpart, while the R07 satellite had the least accurate orbit with a difference of 0.16 m. Clock products had an accuracy of 0.4–1.3 ns for GPS, 1.4–1.6 ns for GLONASS, and 0.7–0.8 ns for QZSS in general. The G15 satellite had the most accurate clock with a difference of only 0.40 ns between the JAXA clock products and MADOCA products, and the R07 satellite had the least accurate clock with a difference of 1.55 ns. The orbit and clock products for GLONASS performed worse than those of GPS and QZSS. (2) After convergence, the positioning accuracy was 3.0–8.1 cm for static PPP and 8.1–13.7 cm for kinematic PPP when using multi-GNSS observations and precise orbit and clock products. The PFRR station performed the good performance both in static and kinematic mode with an accuracy of 2.99 cm and 8.08 cm, respectively, whereas the CPNM station produced the worst static performance with an error of 8.09 cm, and the ANMG station produced the worst kinematic performance with a counterpart of 13.69 cm. (3) The PPP-AR solution was superior to the PPP solution, given that, with respect to PPP, post-processing PPP-AR improved the positioning accuracy and convergence time by 13–32 % (3–89 %) in GPS-only mode by 2–15 % (5–60 %) in GPS/QZSS mode. Thus, we conclude that the current MADOCA products can provide SSR corrections and FCB products with positioning accuracy at the decimeter or even centimeter level, which could meet the demands of the RTPPP and PPP-AR solutions.     

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).     

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

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

Quality variation of GPS satellite clocks on-orbit using IGS clock products

Over 60% clocks on board of the GPS satellites are working longer than their designed life. Therefore realizing their stabilities in a long time scales is essential to GPS navigation and positioning plus IGS time scale maintaining. IGS clock products from 2001 to 2010 are used to analyze the GPS satellite clock qualities such as frequency stabilities and clock noise level. We find out that for the clocks of Block IIA satellites the frequency stabilities and clock noise are 10 times worse than that of the Block IIR and IIR-M satellites. Moreover, the linear relationships between frequency stabilities and clock residuals have been deduced with an accuracy of better than 0.02 ns. Specially, it is noticed that the clock of the PRN27 is instable and the relationship between the frequency stability and residuals is at least a quadratic curve. Therefore, we suggested that GPS satellite clocks should be weighted by their quality levels in application, and the observations of the Block IIA should not be used for real-time positioning which required precision better than one meter.     

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.     

Performance improvement of real-time PPP ambiguity resolution using a regional integer clock

Obtaining reliable GNSS uncalibrated phase delay (UPD) or integer clock products is the key to achieving ambiguity-fixed solutions for real-time (RT) precise point positioning (PPP) users. However, due to the influence of RT orbit errors, the quality of UPD/integer clock products estimated with a globally distributed GNSS network is difficult to ensure, thereby affecting the ambiguity resolution (AR) performance of RT-PPP. In this contribution, by fully utilising the consistency of orbital errors in regional GNSS network coverage areas, a method is proposed for estimating regional integer clock products to compensate for RT orbit errors. Based on Centre National d’Études Spatiales (CNES) RT precise products, the regional GPS/BDS integer clock was estimated with a CORS network in the west of China. Results showed that the difference between the estimated regional and CNES global integer clock/bias products was generally less than 5 cm for GPS, whereas clock differences of greater than 10 cm were observed for BDS due to the large RT orbit error. Compared with PPP using global integer clock/bias products, the AR performance of PPP using the regional integer clock was obviously improved for four rover stations. For single GPS, the horizontal and vertical accuracies of ambiguity-fixed PPP solutions were improved by 56.2% and 45.3% on average, respectively, whereas improvements of 67.5% and 50.5% in the horizontal and vertical directions, respectively, were observed for the combined GPS/BDS situation. Based on a regional integer clock, the RMS error of a kinematic ambiguity-fixed PPP solution in the horizontal direction could reach 0.5 cm. In terms of initialisation time, the average time to first fix (TTFF) in combined GPS/BDS PPP was shortened from 18.2 min to 12.7 min. In view of the high AR performance realised with the use of regional integer clocks, this method can be applied to scenarios that require high positioning accuracy, such as deformation monitoring.     

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.     

BDS-2/BDS-3 uncalibrated phase delay estimation considering the intra-system bias

Intra-system biases (ISBs) between BDS-2 and BDS-3 are of critical importance when combining observations from the BDS-2 and BDS-3 systems, which is meaningful to fully take advantage of the BDS positioning capability. Meanwhile, ISBs should also be considered in the estimation of BDS uncalibrated phase delays (UPDs). In this research, we present a BDS-2/BDS-3 joint-processing scheme, as well as a method for estimating BDS UPDs. The characteristics of ISBs and the quality of BDS UPDs are analyzed based on 30-day data from 130 multi-GNSS experimental (MGEX) stations. Our results indicate that the ISBs are related to the type and version of the receiver. The ISBs can be regarded as constant across the course of a given day, and the mean standard deviation (STD) values of ISBs over one month for different types of receivers are generally within 0.2 m. Moreover, to assess the quality of UPD products, the residuals of the estimated UPDs and the utilization rates of the observation data are computed. The results show that the quality of BDS UPDs can be improved by correcting the satellite-induced pseudo-range variations, and by estimating the wide-lane (WL) UPD difference between BDS-2 and BDS-3. The average RMS values of the estimated residuals of WL UPD and narrow-lane (NL) UPD are 0.07 and 0.09 cycles, respectively; moreover, the utilization rate of the observation data of WL UPD and NL UPD can reach above 90 %. The performance of BDS precise point positioning (PPP) and PPP ambiguity resolution (PPP-AR) is analyzed in terms of positioning accuracy and convergence performance in both the static and kinematic modes. Compared with PPP ambiguity-float solutions, the positioning accuracy of PPP-AR is significantly improved, especially in the east direction. The impact of ISBs on PPP and PPP-AR is also analyzed, and the results indicate that ISBs can improve the convergence speed of float PPP, but can be disregarded in PPP-AR.     

Performance analysis of PPP ambiguity resolution with UPD products estimated from different scales of reference station networks

Integer ambiguity fixing with uncalibrated phase delay (UPD) products can significantly shorten the initialization time and improve the accuracy of precise point positioning (PPP). Since the tracking arcs of satellites and the behavior of atmospheric biases can be very different for the reference networks with different scales, the qualities of corresponding UPD products may be also various. The purpose of this paper is to comparatively investigate the influence of different scales of reference station networks on UPD estimation and user ambiguity resolution. Three reference station networks with global, wide-area and local scales are used to compute the UPD products and analyze their impact on the PPP-AR. The time-to-first-fix, the unfix rate and the incorrect fix rate of PPP-AR are analyzed. Moreover, in order to further shorten the convergence time for obtaining precise positioning, a modified partial ambiguity resolution (PAR) and corresponding validation strategy are presented. In this PAR method, the ambiguity subset is determined by removing the ambiguity one by one in the order of ascending elevations. Besides, for static positioning mode, a coordinate validation strategy is employed to enhance the reliability of the fixed coordinate. The experiment results show that UPD products computed by smaller station network are more accurate and lead to a better coordinate solution; the PAR method used in this paper can shorten the convergence time and the coordinate validation strategy can improve the availability of high precision positioning.     

The update of BDS-2 TGD and its impact on positioning

Timing group delay (TGD) is an important parameter that affects the positioning performance of global navigation satellite systems (GNSS). The BeiDou navigation satellite system (BDS) broadcasts TGD corrections from B3I frequency to B1I and B2I frequencies, namely TGD₁ and TGD₂. On July 21, 2017, BDS updated TGD values with a maximum change of more than 4 ns. In this contribution, we explain the motivation for the BDS TGD update, which is due to the systematic bias between narrowly correlated and widely correlated pseudo-ranges in BDS monitoring receivers. To investigate the impact of the updated TGD, BDS signal-in-space range error (SISRE) and user positioning performance regarding single point positioning (SPP) and precise point positioning (PPP) are analyzed. Results show that after the update of TGD, the difference between the new TGD and multi-GNSS experiment (MGEX) differential code bias (DCB) decreases from 1.38 ns to 0.29 ns on TGD₁ and from 0.40 ns to 0.25 ns on TGD₂. With the contribution of more accurate TGD, the systematic bias of BDS radial SISRE no longer exists, and the overall BDS SISRE also reduces from 1.33 m to 0.87 m on B1I/B2I frequency, from 1.05 m to 0.89 m on B1I frequency, from 0.92 m to 0.91 m on B2I frequency, respectively, which proves the similar precision of BDS TGD and MGEX DCB. One week of statistical results from 28 globally distributed MGEX stations shows that the SPP performance improves on non-B3I frequencies after the TGD update, with a maximum improvement of more than 22% for the B1I/B2I or B1I/B3I combination. The new TGD mainly reduces SPP positioning bias in the East component. The updated TGD also slightly improves the PPP convergence performance for the B1I/B3I combination, but mostly contributes to a more accurate estimation of the receiver clock and ambiguities.     

The GFZ real-time GNSS precise positioning service system and its adaption for COMPASS

Motivated by the IGS real-time Pilot Project, GFZ has been developing its own real-time precise positioning service for various applications. An operational system at GFZ is now broadcasting real-time orbits, clocks, global ionospheric model, uncalibrated phase delays and regional atmospheric corrections for standard PPP, PPP with ambiguity fixing, single-frequency PPP and regional augmented PPP. To avoid developing various algorithms for different applications, we proposed a uniform algorithm and implemented it into our real-time software. In the new processing scheme, we employed un-differenced raw observations with atmospheric delays as parameters, which are properly constrained by real-time derived global ionospheric model or regional atmospheric corrections and by the empirical characteristics of the atmospheric delay variation in time and space. The positioning performance in terms of convergence time and ambiguity fixing depends mainly on the quality of the received atmospheric information and the spatial and temporal constraints. The un-differenced raw observation model can not only integrate PPP and NRTK into a seamless positioning service, but also syncretize these two techniques into a unique model and algorithm. Furthermore, it is suitable for both dual-frequency and sing-frequency receivers. Based on the real-time data streams from IGS, EUREF and SAPOS reference networks, we can provide services of global precise point positioning (PPP) with 5–10 cm accuracy, PPP with ambiguity-fixing of 2–5 cm accuracy, PPP using single-frequency receiver with accuracy of better than 50 cm and PPP with regional augmentation for instantaneous ambiguity resolution of 1–3 cm accuracy. We adapted the system for current COMPASS to provide PPP service. COMPASS observations from a regional network of nine stations are used for precise orbit determination and clock estimation in simulated real-time mode, the orbit and clock products are applied for real-time precise point positioning. The simulated real-time PPP service confirms that real-time positioning services of accuracy at dm-level and even cm-level is achievable with COMPASS only.     

Assessing the performance of BDS-3 for multi-GNSS static and kinematic PPP-AR

As of 2021, a total of four different GNSS constellations – namely, GPS, GLONASS, Galileo, and BDS-3 – can be used with Full Operational Capability (FOC). In this work, the contribution of BDS-3 FOC to GPS + GLONASS + Galileo (GRE) PPP-AR is investigated, considering the three different cut-off angles (7°, 30°, and 45°) and different lengths of static observation sessions (24-, 12-, 6-, 3-, 1-, 0.5-, 0.25-hour). The data of 31 IGS-MGEX stations is processed with GRE PPP-AR and GREC3 (GPS + GLONASS (using float mode) + Galileo + BDS-3) PPP-AR modes. The results showed that BDS-3 degraded the horizontal (except for 24-h sessions) and vertical accuracy of static GRE PPP-AR solutions regardless of the elevation cutoff angle and observation time. The kinematic results showed that BDS-3 significantly contributed to the accuracy of GRE kinematic PPP-AR for 30° and 45° cutoff angles. The convergence time analysis showed that BDS-3 only contributes to GRE kinematic PPP-AR for the vertical component.     

Preliminary analysis of the quality and positioning performance of BDS-3 global interoperable signal B1C&B2a

BeiDou-3 Navigation Satellite System (BDS-3) satellites are equipped with the new generation GNSS signals B1C and B2a, which support the interoperability with GPS and Galileo systems. In this study, the pseudo-range multipath error and carrier phase observation noise of the BDS-3 B1C and B2a signals were evaluated based on zero baseline measurements from the day of year (DOY) 113 to 116, 2020. Further, the precision and performance of the single point positioning (SPP) and precise point positioning (PPP) are assessed at 9 stations. This assessment manifests that the standard deviations (STDs) of the pseudo-range multipath error are about 0.09 ~ 0.22 m, while STDs of the carrier phase observation noise are about 0.075 mm. For the single-frequency SPP, its positioning precision is about 2.03 ~ 4.85 m and 3.29 ~ 10.73 m at the 99.99% confidence level in horizontal and vertical directions, respectively, while the dual-frequency SPP precision is about 1.92 ~ 8.02 m and 4.81 ~ 12.77 m in horizontal and vertical directions, respectively. For the daily static PPP, the convergence time is about 20 ~ 30 min, while the daily positioning precision can reach 1.38 ~ 4.42 cm and -1.31 ~ 4.34 cm in horizontal and vertical directions, respectively. In general, the quality and the SPP and PPP performance of the BDS-3 B1C&B2a signals are comparable to the GPS and Galileo.     

Galileo real-time orbit determination with multi-frequency raw observations

Real-time GNSS-based applications require corresponding real-time orbit products. While traditional GNSS orbits are generated with the dual-frequency IF (Ionosphere-Free) model, the increase of multi-frequency signal satellites brings new challenges for the data processing. Therefore, real-time orbit determination with the multi-frequency UC (Uncombined) model is introduced in this study considering its flexibility. With the derived mathematical model conforming to IGS (International GNSS Service) dual-frequency clock definition and one-week triple-frequency Galileo observation data from 90 IGS network stations, the convergence and accuracy of real-time orbits is assessed and the characteristics of satellite IFCB (Inter-Frequency Clock Bias) are analyzed. Results indicate that the model differences, including dual-frequency IF model, dual-frequency UC model and triple-frequency UC model, contribute to only cm-level differences with CODE (Center for Orbit Determination in Europe) final orbits after a convergence time of around 12 h. The constellation-mean RMS (Root Mean Square) differences of the converged real-time orbits with the CODE final orbits reaches about 5.0 cm, 7.0 cm and 5.0 cm for the radial, tangential and normal directions. The convergence of satellite IFCB is much faster than that of satellite orbit, which reflects a loose correlation between these two parameters. While the Galileo satellite IFCB are temporally stable, the modeling of satellite IFCB may be unreliable when over constrained and becomes even more unstable with commonly encountered datum changes. In summary, real-time GNSS orbit determination with multi-frequency raw observations is feasible and extendable with proper treatment of IFCB.     

A real-time ionospheric model based on GNSS Precise Point Positioning

This paper proposes a method of real-time monitoring and modeling the ionospheric Total Electron Content (TEC) by Precise Point Positioning (PPP). Firstly, the ionospheric TEC and receiver’s Differential Code Biases (DCB) are estimated with the undifferenced raw observation in real-time, then the ionospheric TEC model is established based on the Single Layer Model (SLM) assumption and the recovered ionospheric TEC. In this study, phase observations with high precision are directly used instead of phase smoothed code observations. In addition, the DCB estimation is separated from the establishment of the ionospheric model which will limit the impacts of the SLM assumption impacts. The ionospheric model is established at every epoch for real time application. The method is validated with three different GNSS networks on a local, regional, and global basis. The results show that the method is feasible and effective, the real-time ionosphere and DCB results are very consistent with the IGS final products, with a bias of 1–2 TECU and 0.4 ns respectively.     

Contribution analysis of QZSS to single-frequency PPP of GPS/BDS/GLONASS/Galileo

The Quasi-Zenith Satellite System (QZSS) established by the Japan Aerospace Exploration Agency mainly serves the Asia-Pacific region and its surrounding areas. Currently, four in-orbit satellites provide services. Most users of GNSS in the mass market use single-frequency (SF) receivers owing to the low cost. Therefore, it is meaningful to analyze and evaluate the contribution of the QZSS to SF precise point positioning (PPP) of GPS/BDS/GLONASS/Galileo systems with the emergence of GNSS and QZSS. This study compares the performances of three SF PPP models, namely the GRoup and PHase Ionospheric Correction (GRAPHIC) model, GRAPHIC with code observation model, and an ionosphere-constrained model, and evaluated the contribution of the QZSS to the SF PPP of GPS/BDS/GLONASS/Galileo systems. Moreover, the influence of code bias on the SF PPP of the BDS system is also analyzed. A two-week dataset (DOY 013–026, 2019) from 10 stations of the MGEX network is selected for validation, and the results show that: (1) For cut-off elevation angles of 15, 20, and 25°, the convergence times for the static SF PPP of GLONASS + QZSS are reduced by 4.3, 30.8, and 12.7%, respectively, and the positioning accuracy is similar compared with that of the GLONASS system. Compared with the BDS single system, the convergence times for the static SF PPP of BDS + QZSS under 15 and 25° are reduced by 37.6 and 39.2%, the horizontal positioning accuracies are improved by 18.6 and 14.1%, and the vertical components are improved by 13.9 and 21.4%, respectively. At cut-off elevation angles of 15, 20, and 25°, the positioning accuracy and precision of GPS/BDS/GLONASS/Galileo + QZSS is similar to that of GPS/BDS/GLONASS/Galileo. And the convergence times are reduced by 7.4 and 4.3% at cut-off elevation angles of 20 and 25°, respectively. In imitating dynamic PPP, the QZSS significantly improves the positioning accuracy of BDS and GLONASS. However, QZSS has little effect on the GPS-only, Galileo-only and GPS/BDS/GLONASS/Galileo. (2) The code bias of BDS IGSO and MEO cannot be ignored in SF PPP. In static SF PPP, taking the frequency band of B1I whose multipath combination is the largest among the frequency bands as an example, the vertical component has a systematic bias of approximately 0.4–1.0 m. After correcting the code bias, the positioning error in the vertical component is lower than 0.2 m, and the positioning accuracy in the horizontal component are improved accordingly. (3) The SF PPP model with ionosphere constraints has a better convergence speed, while the positioning accuracy of the three models is nearly equal. Therefore the GRAPHIC model can be used to get good positioning accuracy in the absence of external ionosphere products, but its convergence speed is slower.     

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.     

Particle filter-based real-time phase line bias estimation for GNSS-based attitude determination with common-clock receivers

Global navigation satellite system (GNSS)-based attitude determination has been widely adopted in a wide variety of terrestrial, sea, air, and space applications. Recently, the emergence of commercial multi-GNSS common-clock receivers has brought new opportunities for high-precision GNSS-based attitude determination with single-differenced (SD) model. However, the key requirement of using this approach is the accurate estimation of the troublesome line bias (LB) in real-time. In this contribution, we propose a particle filter-based real-time phase LB estimation approach that apply to SD model with single-system single-frequency observations from common-clock receiver. We first analyzed the relationship between the integer ambiguity ratio value and the phase LB. It is proved that the accuracy of a given phase LB value can be qualified by the related ambiguity resolution ratio value, and the normalized ratio value can therefore be used to represent the likelihood function of observations. Then, we presented the particle filter-based real-time phase LB estimation procedure, and assessed its performance using GPS L1/BDS B1I observations from two datasets collected with different types of common-clock receivers in terms of the accuracy and convergence time of phase LB estimation, the computation load, and the positioning and attitude determination accuracy with respect to the double-differenced (DD) model. Experimental results demonstrated that the phase LB could be accurately estimated with short convergence time (generally within 15 epochs). Moreover, compared with the classical DD approach, the particle filter-based SD approach delivers comparable positioning root-mean-square (RMS) errors in the North and East components but significantly smaller RMS errors in the Up component. Accordingly, the achievable yaw accuracy is comparable whereas the pitch accuracy is remarkably improved. The improvements of positioning accuracy in the Up component and pitch accuracy are approximately 35.7 % to 63.7 %, and 33.3 % to 63.1 %, respectively. Additionally, the single-epoch computation time with our particle filter-based SD approach is generally 0.08 s, which is obviously larger than the DD approach but could still meet the requirements of real-time applications below 10 Hz sampling.     

Long-term time-varying characteristics of UPD products generated by a global and regional network and their interoperable application in PPP

The quality and availability of Uncalibrated Phase Delay (UPD) solutions are crucial to the Precise Point Positioning (PPP) service, and the long-term temporal variability and its contributing factors should be better understood. In this paper, we comprehensively investigate the long-term time-varying characteristics of each UPD product respectively generated by a global and regional network and their interoperable application in PPP-AR (ambiguity resolution), the sampling of the WL and NL UPDs are daily and 30 s, respectively. Firstly, in terms of our 30 day Wide-Lane (WL) UPD products of 31 satellites, the Standard Deviation (STD) of each satellite WL UPDs ranges from 0.04 to 0.06 cycles, indicating that the long-term prediction accuracy of satellite WL UPD is sufficient for fixing Wide-Lane ambiguities. Secondly, when a satellite in eclipsing the discontinulity may corrupt the determination of Narrow-Lane (NL) UPD in form of offset, as a result of lacking or poor satellite attitude dynamic modeling. When the influence of discontinuity is removed, the STD of our estimated satellite NL UPDs is less than 0.05 cycles. Thirdly, the STD of our estimated receiver WL UPDs is mainly below 0.2 cycles, which implies that its stability is one order poorer that of the satellite. In addition, if they are used for stations in and around the network covered region, the stability of the UPD products from the CMONOC (Crustal Movement Observation Network of China) is better than that from a global network, benefit from the fact that all the CMONOC stations are equipped with the same receiver type. Finally, the PPP-AR results show that a rate of 82.9% for stations with a WL-ambiguity-fixed rate of over 90% while 69.5% for stations with an NL-ambiguity-fixed rate of over 80% can be achieved when using UPD from the global network, which is worse than that of using UPD from the CMONOC (85.7% for stations with a WL-ambiguity-fixed rate of over 90% while 75% for stations with an NL-ambiguity-fixed rate of over 80%). The results of the experiment on the UPD interoperable application in PPP show that the global network UPD products can provide a fast AR at any single station, and the convergence time is well below 25 min. Particularly, when the location of a station is in and around the regional network, our results show that the PPP results obtained using regional UPDs enable the consistent use of global UPDs. When the location of a station is far away from the regional network, using the regional UPDs can not achieve PPP-AR. Finally, the WL UPDs of the previous day is used for forecasting to estimate the NL UPDs, the stability analysis results of NL UPDs solution and positioning results are demonstrate the validity of forecasted UPD products.