At present, the global reliability and accuracy of Precipitable Water Vapor (PWV) from different reanalysis products have not been comprehensively evaluated. In this study, PWV values derived by 268 Global Navigation Satellite Systems (GNSS) stations around the world covering the period from 2016 to 2018 are used to evaluate the accuracies of PWV values from five reanalysis products. The temporal and spatial evolution is not taken into account in this analysis, although the temporal and spatial evolution of atmospheric flows is one of the most important information elements available in numerical weather prediction products. The evaluation results present that five reanalysis products with PWV accuracy from high to low are in the order of the fifth generation of European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA5), ERA-Interim, Japanese 55-year Reanalysis (JRA-55), National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR), and NCEP/DOE (Department of Energy) according to root mean square error (RMSE), bias and correlation coefficient. The ERA5 has the smallest RMSE value of 1.84 mm, while NCEP/NCAR and NCEP/DOE have bigger RMSE values of 3.34 mm and 3.51 mm, respectively. The findings demonstrate that ERA5 and two NCEP reanalysis products have the best and worst performance, respectively, among five reanalysis products. The differences in the accuracy of the five reanalysis products are mainly attributed to the differences in the spatial resolution of reanalysis products. There are some large absolute biases greater than 4 mm between GNSS PWV values and the PWV values of five reanalysis products in the southwest of South America and western China due to the limit of terrains and fewer observations. The accuracies of five reanalysis products are compared in different climatic zones. The results indicate that the absolute accuracies of five reanalysis products are highest in the polar regions and lowest in the tropics. Furthermore, the effects of different seasons on the accuracies of five reanalysis products are also analyzed, which indicates that RMSE values of five reanalysis products in summer and in winter are the largest and the smallest in the temperate regions. Evaluation results from five reanalysis products can help us to learn more about the advantages and disadvantages of the five released water vapor products and promote their applications.

In this study, we focus on the kinematic precise point positioning (PPP) solutions at high‐latitudes during the March 2015 great geomagnetic storm. We aim to discover the mechanism behind the positioning degradation from the perspective of the impacts of the storm‐induced ionospheric disturbance on the global navigation satellite system (GNSS) data processing. We observed that the phase scintillation dominated the amplitude scintillation at high‐latitudes and the variation pattern of the rate of total electron content index (ROTI) was consistent with that of the phase scintillation during the storm. The kinematic PPP errors at high‐latitudes were almost three times larger than those at the middle‐ and low‐latitude, which were accompanied by large ROTI variations. From the perspective of GNSS data processing, the large positioning errors were also found to be related to the large number of satellites experiencing cycle slips (CSs). Based on the lock time from the ionospheric scintillation monitoring receiver, we found that a large amount of the CSs was falsely detected under the conventional threshold of the CS detector. By increasing such threshold, the kinematic positioning accuracy at high‐latitudes can be improved to obtain similar magnitude as at middle‐ and low‐latitude. The improved positioning accuracy may suggest that the ionospheric disturbance induced by the geomagnetic storm at high‐latitudes has minor effects on triggering the CSs. Therefore, precise positioning can be achieved at high‐latitudes under geomagnetic storms, given that the CS problem is well addressed.

The high-precision prediction of Earth orientation parameters (EOPs) is essential for astro-geodynamics, high-precision space navigation and positioning, on-board autonomous orbits determination and deep space exploration. However, the prediction accuracy of existing models is much lower than the estimation accuracy of geodetic technical solutions, which affects certain high-precision real-time users. To improve the prediction accuracy of EOP in short- and long-term period, we propose a hybrid model by combining the singular spectrum analysis (SSA), least squares (LSs) and support vector machine (SVM) in the study. Through SSA algorithm, the deterministic time-varying signal of EOP time series can be more precisely and reasonably detected and modeled. Based on the optimization theory, we reconstruct the EOP sequences using SSA and establish the LS extrapolation model based on the reconstructed series. Then, the residuals from SSA reconstruction and those from the LS model, are used for SVM training and prediction. The results of two-year prediction experiments based on the EOP 14 C04 series show that the proposed hybrid model has significant improvements in polar motion (PM) and length of day (LOD) for different prediction intervals (1–360 d) compared with the LS + autoregression (AR) model. The prediction error for x-component of polar motion (PMX) is reduced by 40.2%, 31.0% and 51.4% while that for y-component of polar motion (PMY) is 22.1%, 23.3% and 55.6% for prediction period of 30, 90 and 180 d respectively. For LOD, the maximum prediction improvement can reach to 53.8% during the predicted 360 d. In addition, the proposed method has better accuracy in mid- and long-term PM( x, y ) predictions compared to the Bulletin A, with a 360 d prediction error of 27.273 and 21.741 mas.

Universal time (UT1-UTC) is a key component of Earth orientation parameters (EOP), which is important for the study of monitoring the changes in the Earth’s rotation rate, climatic variation, and the characteristics of the Earth. Many existing UT1-UTC prediction models are based on the combination of least squares (LS) and stochastic models such as the Autoregressive (AR) model. However, due to the complex periodic characteristics in the UT1-UTC series, LS fitting produces large residuals and edge distortion, affecting extrapolation accuracy and thus prediction accuracy. In this study, we propose a combined prediction model based on polynomial curve fitting (PCF), weighted least squares (WLS), and AR, namely, the PCF+WLS+AR model. The PCF algorithm is used to obtain accurate extrapolation values, and then the residuals of PCF are predicted by the WLS+AR model. To obtain more accurate extrapolation results, annual and interval constraints are introduced in this work to determine the optimal degree of PCF. Finally, the multiple sets prediction experiments based on the International Earth Rotation and Reference Systems Service (IERS) EOP 14C04 series are carried out. The comparison results indicate that the constrained PCF+WLS+AR model can efficiently and precisely predict the UT1-UTC in the mid and long term. Compared to Bulletin A, the proposed model can improve accuracy by up to 33.2% in mid- and long-term UT1-UTC prediction.

The precise orbit determination (POD) accuracy of the Chinese BeiDou Navigation Satellite System (BDS) is still not comparable to that of the Global Positioning System because of the unfavorable geometry of the BDS and the uneven distribution of BDS ground monitoring stations. Fortunately, low Earth orbit (LEO) satellites, serving as fast moving stations, can efficiently improve BDS geometry. Nearly all studies on Global Navigation Satellite System POD enhancement using large LEO constellations are based on simulations and their results are usually overly optimistic. The receivers mounted on a spacecraft or an LEO satellite are usually different from geodetic receivers and the observation conditions in space are more challenging than those on the ground. The noise level of spaceborne observations needs to be carefully calibrated. Moreover, spaceborne observational errors caused by space weather events, i.e., solar geomagnetic storms, are usually ignored. Accordingly, in this study, the actual spaceborne observation noises are first analyzed and then used in subsequent observation simulations. Then, the observation residuals from the actual-processed LEO POD during a solar storm on 8 September 2017 are extracted and added to the simulated spaceborne observations. The effect of the observational errors on the BDS POD augmented with different LEO constellation configurations is analyzed. The results indicate that the noise levels from the Swarm-A, GRACE-A, and Sentinel-3A satellites are different and that the carrier-phase measurement noise ranges from 2 mm to 6 mm. Such different noise levels for LEO spaceborne observations cause considerable differences in the BDS POD solutions. Experiments calculating the augmented BDS POD for different LEO constellations considering spaceborne observational errors extracted from the solar storm indicate that these errors have a significant influence on the accuracy of the BDS POD. The 3D root mean squares of the BDS GEO, IGSO, and MEO satellite orbits are 1.30 m, 1.16 m, and 1.02 m, respectively, with a Walker 2/1/0 LEO constellation, and increase to 1.57 m, 1.72 m, and 1.32 m, respectively, with a Walker 12/3/1 constellation. When the number of LEO satellites increases to 60, the precision of the BDS POD improves significantly to 0.89 m, 0.77 m, and 0.69 m for the GEO, IGSO, and MEO satellites, respectively. While 12 satellites are sufficient to enhance the BDS POD to the sub-decimeter level, up to 60 satellites can effectively reduce the influence of large spaceborne observational errors, i.e., from solar storms.

Global Navigation Satellite System (GNSS) ultra-rapid orbit is critical for geoscience and real-time engineering applications. To improve the computational efficiency and the accuracy of predicted orbit, a parallel approach for multi-GNSS ultra-rapid orbit determination is proposed based on Message Passing Interface (MPI)/Open Multi Processing (OpenMP). This approach, compared with earlier efficient methods, can improve the efficiency of multi-GNSS ultra-rapid orbit solution without changing the original observation data and retaining the continuity and consistency of the original parameters to be estimated. To obtain high efficiency, three steps are involved in the approach. First and foremost, the normal equation construction is optimized in parallel based on MPI. Second, equivalent reduction of the estimated parameters is optimized using OpenMP parallel method. Third, multithreading is used for parallel orbit extrapolation. Thus, GNSS ultra-rapid orbit determination is comprehensively optimized in parallel, and the computation efficiency is greatly improved. Based on the data from MGEX and IGS stations, experiments are carried out to analyze the performance of the proposed approach in computational efficiency, accuracy and stability. The results show that the approach greatly improves the efficiency of satellite orbit determination. It can realize 1-h update frequency for the multi-GNSS ultra-rapid orbit determination using 88 stations with four-system observations. The accuracy of the GPS, GLONASS, Galileo and BDS ultra-rapid orbit with 1-h update frequency using the parallel approach is approximately 33.4%,31.4%,40.1% and 32.8% higher than that of the original orbit, respectively. The root mean squares (RMS) of GPS, GLONASS, Galileo and BDS predicted orbit are about 3.2 cm, 5.1 cm, 5.6 cm and 11.8 cm. Moreover, the orbit provided by the proposed method has a better stability. The precision loss of all parallel optimization can be negligible and the original correlation between the parameters is fully preserved.