Cong Qiu

The precise point positioning (PPP) service via the B2b signal (PPP-B2b) on the BeiDou Navigation Satellite System (BDS) provides high-accuracy orbit and clock data for global navigation satellite systems (GNSSs), enabling real-time atmospheric data acquisition without internet access. In this study, we assessed the quality of orbit, clock, and differential code bias (DCB) products from the PPP-B2b service, comparing them to post-processed products from various analysis centres. The zenith tropospheric delay (ZTD) and precipitable water vapour (PWV) were computed at 32 stations using the PPP technique with PPP-B2b corrections. These results were compared with post-processed ZTD with final orbit/clock products and ZTD/PWV values derived from the European Centre for Medium-Range Weather Forecasts Reanalysis (ERA5) and radiosonde data. For stations between 30° N and 48° N, the mean root mean square error (RMSE) of ZTD for the PPP-B2b solution was approximately 15 mm compared to ZTD from the International GNSS Service (IGS). However, accuracy declined at stations between 30° N and 38° S, with a mean RMSE of about 25 mm, performing worse than ZTD estimates using Centre National d’Études Spatiales (CNES) products. The mean RMSEs of PWV derived from PPP-B2b were 3.7 mm and 4.4 mm when compared to PWV from 11 co-located radiosonde stations and ERA5 reanalysis, respectively, and underperformed relative to CNES solutions. Seasonal variability in GNSS-derived PWV was also noted. This reduction in accuracy limits the global applicability of PPP-B2b. Despite these shortcomings, satellite-based PPP services like PPP-B2b remain viable alternatives for real-time positioning and atmospheric applications without requiring internet connectivity.

The accurate determination of Atmospheric Boundary Layer Height (ABLH) is crucial in various atmospheric studies and practical applications. In this study, we present a comprehensive comparative analysis of five distinct methods for estimating ABLH using Global Navigation Satellite System (GNSS) Radio Occultation (RO) data. These methods encompass the use of bending angle and refractivity profiles, namely Minimum Gradient methods of the Bending Angle (MGBA) and Refractivity (MGR) profiles, breaking point, Wavelet Covariance Transform (WCT), and Double-Parameter Model Function (DPMF). GNSS-RO data from COSMIC-2 and Spire are used. To establish robust validation, radiosonde data are employed as a reference, ensuring the reliability of our findings. The results reveal notable variations in the performances of these ABLH estimation methods. Specifically, the MGBA, MGR, breaking point, and DPMF methods exhibit strong correlations with the reference data. Conversely, the WCT method displays weaker correlations, higher biases, and elevated root-mean-square-errors, suggesting limitations in capturing the true ABLH. Furthermore, we remove outlier screening to facilitate a comparison of the differences among the five methods. The WCT and DPMF methods can detect strong variations in the profiles near the Earth’s surface and consider them as ABLH. However, these variations are caused by errors. The MGBA method emerges as a reliable and stable option, while the WCT and DPMF methods should be used with caution due to the lower quality of the GNSS-RO profiles near the Earth’s surface.

In this study, we investigate the performances of a commercial Global Navigation Satellite System (GNSS) Radio Occultation (RO) mission and a new-generation RO constellation, i.e., Spire and Constellation Observing System for Meteorology, Ionosphere, and Climate 2 (COSMIC-2), respectively. In the statistical comparison between Spire and COSMIC-2, the results indicate that although the average signal-to-noise ratio (SNR) of Spire is far weaker than that of COSMIC-2, the penetration of Spire is comparable to, and occasionally even better than, that of COSMIC-2. In our analysis, we find that the penetration depth is contingent upon various factors including SNR, GNSS, RO modes, topography, and latitude. With the reanalysis of the European Centre for Medium-Range Weather Forecasts and Radiosonde as the reference data, the identical error characteristics of Spire and COSMIC-2 reveal that overall, the accuracy of Spire’s neutral-atmosphere data products was found to be comparable to that of COSMIC-2.

The precise point positioning service on B2b signal (PPP-B2b) is a real-time decimeter-level positioning service provided by the BeiDou-3 Global Navigation Satellite System (BDS-3). The service provides users with high-precision orbit and clock corrections through geostationary orbit (GEO) satellites, which means that the PPP-B2b service would be unusable if GEO satellites were blocked. In this study, the performance of PPP-B2b corrections and real-time positioning results during outages of the PPP-B2b service are comprehensively investigated. The results showed that PPP can achieve satisfactory accuracy during outages of the PPP-B2b service by extending the nominal validity of the received PPP-B2b corrections. After extending the PPP-B2b corrections for 10 min, for BDS-3 medium earth orbit (MEO) satellites, the mean root-mean-square error (RMSE) values of the extended orbit were 0.16 m, 0.26 m, and 0.23 m in the radial, along-, and cross-track directions, respectively. The accuracy of the BDS-3 inclined geostationary orbit (IGSO) satellites was slightly worse than that of the BDS-3 MEO satellites; for Global Positioning System (GPS) satellites, the mean RMSE values of the extended orbit were 0.11 m, 0.45 m, and 0.33 m in the radial, along-, and cross-track directions, respectively. In terms of the extended clock, the mean standard deviation (STD) reached 0.17 ns, 0.20 ns, and 0.22 ns after 10 min for the BDS-3 MEO, BDS-3 IGSO, and GPS satellites, respectively. The positioning performance maintained with the extended corrections during the PPP-B2b service outage was evaluated based on five stations in and around China. Our experiments showed that, as long as the interruption time does not exceed 10 min, the real-time positioning with extended PPP-B2b corrections can achieve a comparable accuracy with that obtained following PPP-B2b correction.

Nowadays, precipitable water vapor (PWV) retrieved from ground-based Global Navigation Satellite Systems (GNSS) tracking stations has heralded a new era of GNSS meteorological applications, especially for severe weather prediction. Among the existing models that use PWV timeseries to predict heavy precipitation, the “threshold-based” models, which are based on a set of predefined thresholds for the predictors used in the model for predictions, are effective in heavy precipitation nowcasting. In previous studies, monthly thresholds have been widely accepted due to the monthly patterns of different predictors being fully considered. However, the primary weakness of this type of thresholds lies in their poor prediction results in the transitional periods between two consecutive months. Therefore, in this study, a new method for the determination of an optimal set of diurnal thresholds by adopting a 31-day sliding window was first proposed. Both the monthly and diurnal variation characteristics of the predictors were taken into consideration in the new method. Then, on the strength of the new method, an improved PWV-based model for heavy precipitation prediction was developed using the optimal set of diurnal thresholds determined based on the hourly PWV and precipitation records for the summer over the period 2010–2017 at the co-located HKSC–KP (King’s Park) stations in Hong Kong. The new model was evaluated by comparing its prediction results against the hourly precipitation records for the summer in 2018 and 2019. It is shown that 96.9% of heavy precipitation events were correctly predicted with a lead time of 4.86 h, and the false alarms resulting from the new model were reduced to 25.3%. These results suggest that the inclusion of the diurnal thresholds can significantly improve the prediction performance of the model.

Nowadays, the Global Navigation Satellite Systems (GNSS) have become an effective atmospheric observing technique to remotely sense precipitable water vapor (PWV) mainly due to their high spatiotemporal resolutions. In this study, from an investigation for the relationship between GNSS-derived PWV (GNSS-PWV) and heavy precipitation, it was found that from several hours before heavy precipitation, PWV was probably to start with a noticeable increase followed by a steep drop. Based on this finding, a new model including five predictors for heavy precipitation prediction is proposed. Compared with the existing 3-factor model that uses three predictors derived from the ascending trend of PWV time series (i.e., PWV value, PWV increment and rate of the PWV increment), the new model also includes two new predictors derived from the descending trend: PWV decrement and rate of PWV decrement. The use of the two new predictors for reducing the number of misdiagnosis predictions is proposed for the first time. The optimal set of monthly thresholds for the new five-predictor model in each summer month were determined based on hourly GNSS-PWV time series and precipitation records at three co-located GNSS/weather stations during the 8-year period 2010–2017 in the Hong Kong region. The new model was tested using hourly GNSS-PWV and precipitation records obtained at the above three co-located stations during the summer months in 2018 and 2019. Results showed that 189 of the 198 heavy precipitation events were correctly predicted with a lead time of 5.15 h, and the probability of detection reached 95.5%. Compared with the 3-factor method, the new model reduced the FAR score by 32.9%. The improvements made by the new model have great significance for early detection and predictions of heavy precipitation in near real-time.

Global navigation satellite systems (GNSSs) have become an important tool to derive atmospheric products, such as the total zenith tropospheric delay (ZTD) and precipitable water vapor (PWV) for weather and climate studies. The ocean tide loading (OTL) effect is one of the primary errors that affects the accuracy of GNSS-derived ZTD/PWV, which means the study and choice of the OTL model is an important issue for high-accuracy ZTD estimation. In this study, GNSS data from 1 January 2019 to 31 January 2019 are processed using precise point positioning (PPP) at globally distributed stations. The performance of seven widely used global OTL models is assessed and their impact on the GNSS-derived ZTD is investigated by comparing them against the ZTD calculated from co-located radiosonde observations. The results indicate that the inclusion or exclusion of the OTL effect will lead to a difference in ZTD of up to 3–15 mm for island stations, and up to 1–2 mm for inland stations. The difference of the ZTD determined with different OTL models is quite small, with a root-mean-square (RMS) value below 1.5 mm at most stations. The comparison between the GNSS-derived ZTD and the radiosonde-derived ZTD indicates that the adoption of OTL models can improve the accuracy of GNSS-derived ZTD. The results also indicate that the adoption of a smaller cutoff elevation, e.g., 3° or 7°, can significantly reduce the difference between the ZTDs determined by GNSS and radiosonde, when compared against a 15° cutoff elevation. Compared to the radiosonde-derived ZTD, the RMS error of GNSS-derived ZTD is approximately 25–35 mm at a cutoff elevation of 15°, and 15–25 mm when the cutoff elevation is set to 3°.

Global navigation satellite systems (GNSSs) have become an important tool for remotely sensing water vapor in the atmosphere. In GNSS data processing, mapping functions and gradient models are needed to map the zenith tropospheric delay (ZTD) to the slant total tropospheric delay (STD) along a signal path. Therefore, it is essential to investigate the spatial–temporal performance of various mapping functions and gradient models in the determination of STD. In this study, the STDs at nine elevations were first calculated by applying the ray-tracing method to the atmospheric European Reanalysis-Interim (ERA—Interim) dataset. These STDs were then used as the reference to study the accuracy of the STDs that determined the ZTD together with mapping functions and gradient models. The performance of three mapping functions (i.e., Niell mapping function (NMF), global mapping function (GMF), and Vienna mapping function (VMF1)) and three gradient models (i.e., Chen, MacMillan, and Meindl) in six regions (the temperate zone, Qinghai–Tibet Plateau, Equator, Sahara Desert, Amazon Rainforest, and North Pole) in determining slant tropospheric delay was investigated in this study. The results indicate that the three mapping functions have relatively similar performance above a 15° elevation, but below a 15° elevation, VMF1 clearly performed better than the GMF and NMF. The results also show that, if no gradient model is included, the root-mean-square (RMS) of the STD is smaller than 2 mm above the 30° elevation and smaller than 9 mm above the 15° elevation but shows a significant increase below the 15° elevation. For example, in the temperate zone, the RMS increases from approximately 35 mm at the 10° elevation to approximately 160 mm at the 3° elevation. The inclusion of gradient models can significantly improve the accuracy of STDs by 50%. All three gradient models performed similarly at all elevations and in all regions. The bending effect was also investigated, and the results indicate that the tropospheric delay caused by the bending effect is normally below 13 mm above a 15° elevation, but this delay increases dramatically from approximately 40 mm at a 10° elevation to approximately 200 mm at a 5° elevation, and even reaches 500–700 mm at a 3° elevation in most studied regions.