Typhoons can be serious natural disasters for the sustainability and development of society. The development of a typhoon usually involves a pre-existing weather disturbance, warm tropical oceans, and a large amount of moisture. This implies that a large variation in the atmospheric water vapor over the path of a typhoon can be used to study the characteristics of the typhoon. This is the reason that the variation in precipitable water vapor (PWV) is often used to capture the signature of a typhoon in meteorology. This study investigates the usability of real-time PWV retrieved from global navigation satellite systems (GNSS) for typhoons’ characterizations, and especially, the following aspects were investigated: (1) The correlation between PWV and atmospheric parameters including pressure, temperature, precipitation, and wind speed; (2) water vapor transportation during a typhoon period; and (3) the correlation between the movement of a typhoon and the transportation of water vapor. The case study selected for this research was Super Typhoon Mangkhut that occurred in mid-September 2018 in Hong Kong. The PWV time series were obtained from a conversion of GNSS-derived zenith total delays (ZTDs) using observations at 10 stations selected from the Hong Kong GNSS continuously operating reference stations (CORS) network, which are also located along the path of the typhoon. The Bernese GNSS Software (ver. 5.2) was used to obtain the ZTDs; and the root mean square (RMS) of the differences between the GNSS-ZTDs and International GNSS Service post-processed ZTDs time series was less than 8 mm. The RMS of the differences between the GNSS-PWVs (i.e., the ZTDs converted PWVs) and radiosonde-derived PWVs (RS-PWVs) time series was less than 2 mm. The changes in PWV reflect the variation in wind speed during the typhoon period to a certain degree, and their correlation coefficient was 0.76, meaning a significant positive correlation. In addition, a new approach was proposed to estimate the direction and speed of a typhoon’s movement using the time difference of PWV arrival at different sites. The direction and speed estimated agreed well with the ones published by the China Meteorological Administration. These results suggest that GNSS-derived PWV has a great potential for the monitoring and even prediction of typhoon events, especially for near real-time warnings.

Real-time precise point positioning (RT-PPP) has become a powerful technique for the determination of the zenith tropospheric delay (ZTD) over a GPS (global positioning system) or GNSS (global navigation satellite systems) station of interest, and the follow-on high-precision retrieval of precipitable water vapor (PWV). The a priori zenith hydrostatic delay (ZHD) and the mapping function used in the PPP approach are the two factors that could affect the accuracy of the PPP-based ZTD significantly. If the in situ atmospheric pressure is available, the Saastamoinen model can be used to determine ZHD values, and the model-predicted ZHD results are of high accuracy. However, not all GPS/GNSS are equipped with an in situ meteorological sensor. In this research, the daily forecasting ZHD and mapping function values from VMF1 forecasting (VMF1_FC) and VMF3 forecasting (VMF3_FC) products were used for the determination of the GPS-derived PWV. The a priori ZHDs derived from VMF1_FC and VMF3_FC were first evaluated by comparing against the reference ZHDs from globally distributed radiosonde stations. GPS observations from 41 IGS stations that have co-located radiosonde stations during the period of the first half of 2020 were used to test the quality of GPS-ZTD and GPS-PWV. Three sets of ZTDs estimated from RT-PPP solutions using the a priori ZHD and mapping function from the following three VMF products were evaluated: (1) VMF1_FC; (2) VMF3_FC (resolution 5° × 5°); (3) VMF3_FC (resolution 1° × 1°). The results showed that, when the ZHDs from 443 globally distributed radiosonde stations from 1 July 2018 to 30 June 2021 were used as the reference, the mean RMSEs of the ZHDs from the three VMF products were 5.9, 5.4, and 4.3 mm, respectively. The ZTDs estimated from RT-PPP at 41 selected IGS stations were compared with those from IGS, and the results showed that the mean RMSEs of the ZTDs of the 41 stations from the three PPP solutions were 8.6, 9.0, and 8.6 mm, respectively, and the mean RMSEs of the PWV converted from their corresponding ZWDs were 1.9, 2.4, and 1.7 mm, respectively, in comparison with the reference PWV from co-located radiosonde stations. The results suggest that the a priori ZHD and mapping function from VMF1_FC and VMF3_FC can be used for the precise determination of real-time GPS/GNSS-PWV in most regions, especially the VMF3_FC (resolution 1° × 1°) product.

One of the main challenges of Global Navigation Satellite System (GNSS) tomography is in solving ill-conditioned system equations. Vertical constraint models are typically used in the solution procedure and play an important role in the quality of the GNSS tomography, in addition to helping resolve ill-posed problems in system equations. In this study, based on a water vapor (WV) parameter, namely IRPWV, a new vertical constraint model with six sets of coefficients for six different WV states was developed and tested throughout 2019 in the Hong Kong region with four tomographic schemes, which were carried out with the model and the traditional vertical constraint model using three different types of water vapor scale height parameters. Experimental results were numerically compared against their corresponding radiosonde-derived WV values. Compared with the tests that used the traditional model, our results showed that, first, for the daily relative error of WV density (WVD) less than 30%, the new model can lead to at least 10% and 49% improvement on average at the lower layers (below 3 km, except for the ground surface) and the upper layers (about 5–10 km), respectively. Second, the skill score of the monthly root-mean-square error (RMSE) of layered WVD above 10 accounted for about 83%, 87%, and 64%. Third, for the annual biases of layered WVD, the new model significantly decreased by 1.1–1.5 g/m³ at layers 2–3 (about 1 km), where all schemes showed the maximal bias value. Finally, for the annual RMSE of layered WVD, the new model at the lower (about 0.6–3 km) and upper layers improved by 13–42% and 5–47%, respectively. Overall, the new model performed better on GNSS tomography and significantly improved the accuracy of GNSS tomographic results, compared to the traditional model.

The tropospheric tomography is an ill-posed inversion problem due to the sparsity of global navigation satellite systems (GNSS) stations and the limitation on the projection angles of GNSS signals, which in turn affects the stability and robustness of the tomographic solution. To address this, a new tomographic algorithm, named two-step projected iterative algorithm (TSPIA), is proposed. The wet refractivity (WR) field was constructed in two steps: first, an iterative preprocessing for the initial input values was performed, and then its resultant solution was input into the projected iterative method, in which a hypothesis convex set was constructed to constrain the reconstruction based on the classical algebraic iterative reconstruction (AIR) methods. In addition, a two-dimensional normalized cumulative periodogram (2-D-NCP) termination criterion was investigated since the traditional criteria for judging the convergence of iterations use prefixed empirical thresholds, which may lead to excessive iterations and need complicated work. The TSPIA was tested using GNSS data in Hong Kong over a wet period and a dry period. Statistical results showed that, compared to the classical AIR methods, the accuracy of the reconstructed WR field of the TSPIA were improved by about 10% and 15% when radiosonde and ECMWF data were used as the reference, respectively. Moreover, experiments for the proposed 2-D-NCP criterion demonstrated noticeable computational efficiency. These results suggest that the new approaches proposed in this article can improve the performance of the iterative methods for GNSS tropospheric tomography.

For better modeling the vertical distribution and variations of water vapor, a 10-year time-series of a newly derived water vapor parameter (termed IRPWV), defined as the ratio of water vapor density (WVD) to the total precipitable water vapor (TPWV), was statistically analyzed. This research showed that the vertical distribution of IRPWV presents a periodic pattern and is highly correlated with the relative magnitude of its corresponding TPWV when compared with the other TPWVs in the same time range. Six TPWV ranges were first chosen to determine the relative magnitude and then used to classify the IRPWV vertical distributions of TPWV. For the periodic variations in each of the six classified IRPWV vertical distribution time-series, a temporal IRPWV model was developed accordingly with six sets of coefficients. The new models were validated by comparing their predictions against the reference values from sounding data at 12 radiosonde stations in China, and their performance was also evaluated against that of the commonly used exponential model. Results showed that, first, the proportions of the height range that had reduced annual root mean square error (RMSE) of WVD in all height ranges within all TPWV ranges were over 75% at the 12 stations. Then, the annual RMSEs of the WVD for all the stations were reduced by at least 11%, 20%, 43%, 48%, 40%, 38%, 32%, 35%, 32%, and 28% in each of the 10 selected height ranges, respectively.