The satellite clocks used in the BeiDou-2 satellite navigation System (BDS) are Chinese self-developed Rb atomic clocks, and their performances and stabilities are worse than GPS and Galileo satellite clocks. Due to special periodic noises and nonlinear system errors existing in the BDS clock offset series, the GPS ultra-rapid clock model, which uses a simple quadratic polynomial plus one periodic is not suitable for BDS. Therefore, an improved prediction model for BDS satellite clocks is proposed in order to enhance the precision of ultra-rapid predicted clock offsets. First, a basic quadratic polynomial model which is fit for the rubidium (Rb) clock is constructed for BDS. Second, the main cyclic terms are detected and identified by the Fast Fourier Transform (FFT) method according to every satellite clock offset series. The detected results show that most BDS clocks have special cyclic terms which are different from the orbit periods. Therefore, two main cyclic terms are added to absorb the periodic effects. Third, after the quadratic polynomial plus two periodic fitting, some evident nonlinear system errors also exist in the model residual, and the Back Propagation (BP) neural network model is chosen to compensate for these nonlinear system errors. The simulation results show that the performance and precision using the improved model are better than that of China iGMAS ultra-rapid prediction (ISU-P) products and the Deutsches GeoForschungsZentrum GFZ BDS ultra-rapid prediction (GBU-P) products. Comparing to ISU-P products, the average improvements using the proposed model in 3 h, 6 h, 12 h and 24 h are 23.1%, 21.3%, 20.2%, and 19.8%, respectively. Meanwhile the accuracy improvements of the proposed model are 9.9%, 13.9%, 17.3%, and 21.2% compared to GBU-P products. In addition, the kinematic Precise Point Positioning (PPP) example using 8 Multi-GNSS Experiment MGEX stations shows that the precision based on the proposed clock model has improved about 16%, 14%, and 38% in the North (N), East (E) and Height (H) components.

With the high-precision products of satellite orbit and clock, uncalibrated phase delay, and the atmosphere delay corrections, Precise Point Positioning (PPP) based on a Real-Time Kinematic (RTK) network is possible to rapidly achieve centimeter-level positioning accuracy. In the ionosphere-weighted PPP–RTK model, not only the a priori value of ionosphere but also its precision affect the convergence and accuracy of positioning. This study proposes a method to determine the precision of the interpolated slant ionospheric delay by cross-validation. The new method takes the high temporal and spatial variation into consideration. A distance-dependent function is built to represent the stochastic model of the slant ionospheric delay derived from each reference station, and an error model is built for each reference station on a five-minute piecewise basis. The user can interpolate ionospheric delay correction and the corresponding precision with an error function related to the distance and time of each reference station. With the European Reference Frame (EUREF) Permanent GNSS (Global Navigation Satellite Systems) network (EPN), and SONEL (Système d'Observation du Niveau des Eaux Littorales) GNSS stations covering most of Europe, the effectiveness of our wide-area ionosphere constraint method for PPP-RTK is validated, compared with the method with a fixed ionosphere precision threshold. It is shown that although the Root Mean Square (RMS) of the interpolated ionosphere error is within 5 cm in most of the areas, it exceeds 10 cm for some areas with sparse reference stations during some periods of time. The convergence time of the 90th percentile is 4.0 and 20.5 min for horizontal and vertical directions using Global Positioning System (GPS) kinematic solution, respectively, with the proposed method. This convergence is faster than those with the fixed ionosphere precision values of 1, 8, and 30 cm. The improvement with respect to the latter three solutions ranges from 10 to 60%. After integrating the Galileo navigation satellite system (Galileo), the convergence time of the 90th percentile for combined kinematic solutions is 2.0 and 9.0 min, with an improvement of 50.0% and 56.1% for horizontal and vertical directions, respectively, compared with the GPS-only solution. The average convergence time of GPS PPP-RTK for horizontal and vertical directions are 2.0 and 5.0 min, and those of GPS + Galileo PPP-RTK are 1.4 and 3.0 min, respectively.

The performance of precise point positioning (PPP) has been significantly improved thanks to the continuous improvements in satellite orbit, clock, and ambiguity resolution (AR) technologies, but the convergence speed remains a limiting factor in real-time PPP applications. To improve the PPP precision and convergence time, tropospheric delays from a regional network can be modeled to provide precise correction for users. We focus on the precise modeling of zenith wet delay (ZWD) over a wide area with large altitude variations for improving PPP-AR. By exploiting the water vapor exponential vertical decrease, we develop a modified optimal fitting coefficients (MOFC) model based on the traditional optimal fitting coefficients (OFC) model. The proposed MOFC model provides a precision better than 1.5 cm under sparse inter-station distances over the Europe region, with a significant improvement of 70% for high-altitude stations compared to the OFC model. The MOFC model with different densities of reference stations is further evaluated in GPS and Galileo kinematic PPP-AR solutions. Compared to the PPP-AR solutions without tropospheric delay augmentation, the positioning precision of those with 100-km inter-station spacing MOFC and OFC is improved by 25.7% and 17.8%, respectively, and the corresponding time to first fix (TTFF) is improved by 36.9% and 33.0% in the high-altitude areas. On the other hand, the OFC model only slightly improves the TTFF and positioning accuracy when using the 200 km inter-station spacing modeling and even degrades the positioning for high-altitude stations, whereas using the MOFC model, the PPP-AR solutions always improve. Moreover, the positioning precision improvement of MOFC compared with OFC is about 22.1%, 21.7%, and 25.7% for the Galileo-only, GPS-only, and GPS + Galileo PPP-AR solutions, respectively.

Wide-lane (WL) uncalibrated phase delay (UPD) is usually derived from Melbourne–Wübbena (MW) linear combination and is a prerequisite in Global Navigation Satellite Systems (GNSS) precise point positioning (PPP) ambiguity resolution (AR). MW is a linear combination of pseudorange and phase, and the accuracy is limited by the larger pseudorange noise which is about one hundred times of the carrier phase noise. However, there exist inconsistent pseudorange biases which may have detrimental effect on the WL UPD estimation, and further degrade user-side ambiguity fixing. Currently, only the large part of pseudorange biases, e.g., the differential code bias (DCB), are available and corrected in PPP-AR, while the receiver-type-dependent biases have not yet been considered. Ignoring such kind of bias, which could be up to 20 cm, will cause the ambiguity fixing failure, or even worse, the incorrect ambiguity fixing. In this study, we demonstrate the receiver-type-dependent WL UPD biases and investigate their temporal and spatial stability, and further propose the method to precisely estimate these biases and apply the corrections to improve the user-side PPP-AR. Using a large data set of 1560 GNSS stations during a 30-day period, we demonstrate that the WL UPD deviations among different types of receivers can reach ± 0.3 cycles. It is also shown that such kind of deviations can be calibrated with a precision of about 0.03 cycles for all Global Positioning System (GPS) satellites. On the user side, ignoring the receiver-dependent UPD deviation can cause significant positioning error up to 10 cm. By correcting the deviations, the positioning performance can be improved by up to 50%, and the fixing rate can also be improved by 10%. This study demonstrates that for the precise and reliable PPP-AR, the receiver-dependent UPD deviations cannot be ignored and have to be handled.

Introducing precise atmosphere information into precise point positioning enables rapid ambiguity resolution and introduces a significant accuracy improvement. However, it can only be implemented in regions with dense networks and stable communication links. For larger areas, e.g., an intercontinental level, there is a conflict between the accuracy of corrections and the amount of atmosphere information to be disseminated. We develop a hierarchical augmentation mode to combine the advantages of the fitting model and region interpolation model to relieve the communication burden. Relying on the fitting model with fewer coefficients applied over large areas as the essential information, the unmodeled errors are calculated at each reference station, and further correction information is optional compensation depending on the magnitude of the unmodeled residuals. We perform the proposed models on 103 EUREF Permanent Network stations with 200-km station spacing and 84 stations as the external validation. The ionosphere and troposphere fitting models have an average accuracy of about 4.2 and 1.3 cm, respectively, under meteorologically calm conditions. The unmodeled error transmission determined by the magnitude of residuals can be reduced by 61% and 96% for the ionospheric and tropospheric delays, respectively, with respect to the legacy interpolation mode. Further compensation implemented, i.e., unmodeled residuals, can achieve instantaneous convergence for 83.6% of all solutions, and the overall initialization time is within 1.0 min. Thus, the proposed hierarchical positioning mode satisfies real-time positioning convergence requirements and significantly reduces massive corrections in communication.