Search

The accurate measurement of water vapor is essential for research about and the applications of meteorology, climatology, and hydrology. Based on the BeiDou PPP-B2b service, real-time precipitable water vapor (PWV) can be retrieved with the precise point positioning (PPP) software (XTW-PPP version 0.0). The experiment was conducted in Beijing in January 2023. Three solutions were designed with PPP using the BeiDou system only, the GPS system only, and the BeiDou-GPS combined solution. Real-time PWVs for the three solutions were validated with the ERA5 reanalysis data. Between the PWV values from the single BeiDou and ERA5, there was a bias of 0.7 mm and an RMSE of 1.8 mm. For the GPS case, the bias was 0.73 mm and the RMSE was 1.97 mm. The biases were less than 1 mm and RMSEs were less than 2 mm. Both the BeiDou and the GPS processing performed very well. But little improvement was found for the BeiDou-GPS combined solution, compared with the BeiDou system-only and the GPS system-only solution. This may be due to the poor handling of two different kinds of errors for the GPS and the BeiDou systems in our PPP software. A better PWV estimation with the two systems is to estimate PWV with a single system at the first step and then obtain the optimization by Bayesian model averaging.

The Atmospheric Dynamics Mission ADM-Aeolus was successfully launched in August 2018 by the European Space Agency (ESA). The Aeolus mission carried a single instrument, the first-ever Doppler wind lidar (DWL) in space, called Atmospheric LAser Doppler INstrument (ALADIN). Aeolus circled the Earth, providing vertical profiles of horizontal line-of-sight (HLOS) winds on a global scale. The Aeolus satellite’s measurements filled critical gaps in existing wind observations, particularly in remote regions such as the Brazilian Amazon. This area, characterized by dense rainforests and rich biodiversity, is essential for global climate dynamics. The weather patterns of the Amazon are influenced by atmospheric circulation driven by Hadley cells and the Intertropical Convergence Zone (ITCZ), which are crucial for the distribution of moisture and heat from the equator to the subtropics. The data provided by Aeolus can significantly enhance our understanding of these complex atmospheric processes. In this long-term validation study, we used radiosonde data collected from three stations in the Brazilian Amazon (Cruzeiro do Sul, Porto Velho, and Rio Branco) as a reference to assess the accuracy of the Level 2B (L2B) Rayleigh-clear and Mie-cloudy wind products. Statistical validation was conducted by comparing Aeolus L2B wind products and radiosonde data covering the period from October 2018 to March 2023 for Cruzeiro do Sul and Porto Velho, and from October 2018 to December 2022 for Rio Branco. Considering all available collocated winds, including all stations, a Pearson’s coefficient (r) of 0.73 was observed in Rayleigh-clear and 0.85 in Mie-cloudy wind products, revealing a strong correlation between Aeolus and radiosonde winds, suggesting that Aeolus wind products are reliable for capturing wind profiles in the studied region. The observed biases were −0.14 m/s for Rayleigh-clear and −0.40 m/s for Mie-cloudy, fulfilling the mission requirement of having absolute biases below 0.7 m/s. However, when analyzed annually, in 2022, the bias for Rayleigh-clear was −0.95 m/s, which did not meet the mission requirements.

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.

This study aims to elucidate the influence of meteorological conditions on particle levels in Cergy-Pontoise. It explores the temporal variability of PM10 pollution days by associating them with the vertical temperature profile derived from conventional radiosondes from 2013 to 2022 (regional station). The results indicate that nearly 80% of exceedance days were associated with thermal inversions, primarily observed in winter and typically lasting 1 to 3 days. Analysis of winter thermal inversion characteristics suggests that those linked to pollution primarily occur near the ground, with higher intensity in December (12.1 °C) and lower in February (10.3 °C). Persistent inversions (extended nocturnal by diurnal inversion) account for 91.4% of the total inversions associated with high concentrations. Captive balloon soundings and temperature measurements at different altitudes were conducted during the winter of 2022/2023 to clarify thermal inversion in the Oise Valley at the center of Cergy-Pontoise. The results highlight three nocturnal wind circulation mechanisms in the valley, including downslope flow, circulation influenced by an urban heat island, and mechanical air evacuation under an inversion layer towards the less steep East side of the valley. Analysis of PM with the temperature gradient in the Oise Valley shows a significant correlation, suggesting an increase in concentrations during locally detected inversions and a decrease during atmospheric disturbance.

Significant knowledge gaps exist in our understanding of urban boundary layer processes, particularly the hygrothermal state. The earth system community has successfully used microwave radiometers for several decades. However, the applicability in complex urban environments has never been adequately tested. Here, observations from a microwave radiometer are compared to radiosonde readings in a densely urbanized site in Houston, Texas. The site was influenced by both an urban heat island and the sea breeze phenomenon. The analysis showed significant disagreement between the virtual potential temperature predicted by the microwave radiometer and the radiosonde for all periods within the boundary layer. However, the values were reasonably comparable above the boundary layer. The microwave radiometer incorrectly predicted an inversion layer instead of a mixed layer during convective periods. The microwave radiometer measurements deviated from the radiosonde measurements throughout the lower troposphere for the relative humidity.

The geographical location of Yunnan province is at the upstream area of water vapor transportation from the Bay of Bengal and the South China Sea to inland China. Understanding the spatiotemporal variations of water vapor over this region holds significant importance. We utilized the Global Navigation Satellite System (GNSS) data collected from 12 stations situated in Yunnan, which are part of the Crustal Movement Observation Network of China, to retrieve hourly precipitable water vapor (PWV) data from 2011 to 2022. The retrieved PWV data at Station KMIN were evaluated by the nearby radiosonde data, and the results show that the mean bias and RMS of the differences between the two datasets are 0.08 and 1.78 mm, respectively. Average PWV values at these stations are in the range of 11.77 to 33.53 mm, which decrease from the southwest to the north of Yunnan and are negatively correlated with the stations’ heights and latitudes. Differences between average PWV in the wet season and dry season range from 12 to 27 mm. These differences tend to increase as the average PWV increases. The yearly rates of PWV variations, averaging 0.18 mm/year, are all positive for the stations, indicating a year-by-year increase in water vapor. The amplitudes of the PWV annual cycles are 9.75–20.94 mm. The spatial variation of these amplitudes is similar to that of the average PWV over the region. Generally, monthly average PWV values increase from January to July and decrease from July to December, and the growth rate is less than the decline rate. Average diurnal PWV variations show unimodal PWV distributions over the course of the day at the stations except Station YNRL, where bimodal PWV distribution was observed.

Multi-Global Navigation Satellite Systems (multi-GNSS) (including GPS, BDS, Galileo, and GLONASS) provide a significant opportunity for high-quality zenith tropospheric delay estimation and its applications in meteorology. However, the performance of zenith total delay (ZTD) retrieval from single- or multi-GNSS observations is not clear, particularly from the new, fully operating BDS-3. In this paper, zenith tropospheric delay is estimated using the single-, dual-, triple-, or four-GNSS Precise Point Positioning (PPP) technique from 55 Multi-GNSS Experiment (MGEX) stations over one year. The performance of GNSS ZTD estimation is evaluated using the International GNSS Service (IGS) standard tropospheric products, radiosonde, and the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis (ERA5). The results show that the GPS-derived ZTD time series is more consistent and reliable than those derived from BDS-only, Galileo-only, and GLONASS-only solutions. The performance of the single-GNSS ZTD solution can be enhanced with better accuracy and stability by combining multi-GNSS observations. The accuracy of the ZTD from multi-GNSS observations is improved by 13.8%, 43.8%, 27.6%, and 22.9% with respect to IGS products for the single-system solution (GPS, BDS, Galileo, and GLONASS), respectively. The ZTD from multi-GNSS observations presents higher accuracy and a significant improvement with respect to radiosonde and ERA5 data when compared to the single-system solution.

The accuracy of temperature and relative humidity (RH) profiles retrieved by the ground-based microwave radiometer (MWR) is crucial for meteorological research. In this study, the four-year measurements of brightness temperature measured by the microwave radiometer from Huangpu meteorological station in Guangzhou, China, and the radiosonde data from the Qingyuan meteorological station (70 km northwest of Huangpu station) during the years from 2018 to 2021 are compared with the sonde data. To make a detailed comparison on the performance of machine learning models in retrieving the temperature and RH profiles, four machine learning algorithms, namely Deep Learning (DL), Gradient Boosting Machine (GBM), Extreme Gradient Boosting (XGBoost) and Random Forest (RF), are employed and verified. The results show that the DL model performs the best in temperature retrieval (with the root-mean-square error and the correlation coefficient of 2.36 and 0.98, respectively), while the RH of the four machine learning methods shows different excellence at different altitude levels. The integrated machine learning (ML) RH method is proposed here, in which a certain method with the minimum RMSE is selected from the four methods of DL, GBM, XGBoost and RF for a certain altitude level. Two cases on 29 January 2021 and on 10 February 2021 are used for illustration. The case on 29 January 2021 illustrates that the DL model is suitable for temperature retrieval and the ML model is suitable for RH retrieval in Guangzhou. The case on 10 February 2021 shows that the ML RH method reaches over 85% before precipitation, implying the application of the ML RH method in pre-precipitation warnings.

In response to the nonlinear fitting difficulty of the traditional weighted average temperature (Tm) modeling, this paper proposed four machine learning (ML)-based Tm models. Based on the seven radiosondes in the Yangtze River Delta region from 2014 to 2019, four forecasting ML-based Tm models were constructed using Light Gradient Boosting Machine (LightGBM), Support Vector Machine (SVM), Random Forest (RF), and Classification and Regression Tree (CART) algorithms. The surface temperature (Ts), water vapor pressure (Es), and atmospheric pressure (Ps) were identified as crucial influencing factors after analyzing their correlations to the Tm. The ML-based Tm models were trained using seven radiosondes from 2014 to 2018. Then, the mean bias and root mean square error (RMSE) of the 2019 dataset were used to evaluate the accuracy of the ML-based Tm models. Experimental results show that the overall accuracy of the LightGBM-based Tm model is superior to the SVM, CART, and RF-based Tm models under different temporal variations. The mean RMSE of the daily LightGBM-based Tm model is reduced by 0.07 K, 0.04 K, and 0.13 K compared to the other three ML-based models, respectively. The mean RMSE of the monthly LightGBM-based Tm model is reduced by 0.09 K, 0.04 K, and 0.11 K, respectively. The mean RMSE of the quarterly LightGBM-based Tm model is reduced by 0.09 K, 0.04 K, and 0.11 K, respectively. The mean bias of the LightGBM-based Tm model is also smaller than that of the other ML-based Tm models. Therefore, the LightGBM-based Tm model can provide more accurate Tm and is more suitable for obtaining GNSS precipitable water vapor in the Yangtze River Delta region.

Global Navigation Satellite System (GNSS) Radio Occultation (RO) data play an essential role in improving numerical weather prediction (NWP) and monitoring climate change. The NOAA Commercial RO Purchase Program (CDP) purchased RO data provided by Spire Global Inc. To ensure the data quality from Spire Global Inc. is consistent with other RO missions, we need to quantify their accuracy and retrieval uncertainty carefully. In this work, Spire Wet Profile (wet temperature profile) data from 7 September 2021 to 31 October 2022, processed by the University Corporation for Atmospheric Research (UCAR), and COSMIC-2 (Constellation Observing System for Meteorology, Ionosphere, and Climate-2/Formosa Satellite Mission 7) data are evaluated through comparison with NOAA-20 Advanced Technology Microwave Sounder (ATMS) microwave sounder measurements and collocated RS41 radiosonde measurements. Through the Community Radiative Transfer Model (CRTM) simulation, we convert the Spire and COSMIC-2 RO retrievals to ATMS brightness temperature (BT) at sounding channels CH07 to CH14 (temperature channels), with weighting function peak heights from 8 km to 35 km, and CH19 to CH22 (water vapor channels), with weighting function peak heights ranging from 3.2 km to 6.7 km, and compare the simulations with the collocated NOAA-20 ATMS measurements over ocean. Using ATMS observations as references, Spire and COSMIC-2 BTs agree well with ATMS within 0.07 K for CH07-14 and 0.20 K for CH19-22. The trends between Spire and COSMIC-2 are consistent within 0.07 K/year over the oceans for ATMS CH07-CH13 and CH19-22, indicating that Spire/COSMIC-2 wet profiles are, in general, compatible with each other over oceans. The RO retrievals and RS41 radiosonde observation (RAOB) comparison shows that above 0.2 km altitude, RS41 RAOB matches Spire/COSMIC-2 temperature profiles well with a temperature difference of

In this study, we devised a constraint method, called multi-model ensemble pattern (MEP), to estimate the refractive index structure parameter (Cn2) profiles based on observational data and multiple existing models. We verified this approach against radiosonde data from field campaigns in China’s eastern and northern coastal areas. Multi-dimensional statistical evaluations for the Cn2 profiles and integrated astronomical parameters have proved MEP’s relatively reliable performance in estimating optical turbulence in the atmosphere. The correlation coefficients of MEP and measurement overall Cn2 in two areas are up to 0.65 and 0.76. A much higher correlation can be found for a single radiosonde profile. Meanwhile, the difference evaluation of integrated astronomical parameters also shows its relatively robust performance compared to a single model. The prowess of this reliable approach allows us to carry out regional investigation on optical turbulence features with routine meteorological data soon.

Constructing accurate models that provide information about water vapor content in the troposphere improves the reliability of numerical weather forecasts and the position accuracy of low-cost Global Navigation Satellite System (GNSS) receivers. However, developing models with high spatial-temporal resolution demands compact observational datasets in the regions of interest. Empirical models, such as the Global Pressure and Temperature 3 (GPT3w), have been constructed based on the monthly averaged outputs of numerical weather models. These models are based on the assimilation of existing measurements to provide estimations of atmospheric parameters. Therefore, their accuracy may be reduced over regions with a low resolution of radiosonde or continuous GNSS stations. By emerging and increasing the Low-Earth-Orbiting (LEO) satellites that measure atmospheric parameter profiles using the Radio Occultation (RO) technique, new opportunities have appeared to acquire high-resolution atmospheric observations at different altitudes. This study aims to apply these RO observations to improve the accuracy of the GPT3w model over Iran, which is sparse in terms of long-term GNSS and radiosonde measurements. The temperature, pressure, and water vapor pressure parameters from the GPT3w model have been used as the input layers of the Extremely Learning Machine (ELM) technique. The wet refractivity indices from the RO technique are considered target parameters in the output layer to train the ELM. The RO observations of 2007–2020 are applied for training, and those of 2020–2022 for evaluating the performance of the developed ELM. Our numerical results indicate that the developed ELM decreases the Root-Mean-Square Error (RMSE) values of the wet refractivity indices by about 17 percent, compared to the original GPT3w RMSE values. Additionally, the wet refractivity indices from ELM have revealed correlation coefficients of about 0.64, which is about 1.9 times those related to the original GPT3w model. The performance of ELM has also been examined by comparison with the data of six located radiosonde stations covering the year 2020. This comparison shows an improvement of about 14 percent in the average RMSE values of the estimated wet refractivity indices.

Atmospheric precipitable water vapor (PWV) is a key variable for weather forecast and climate research. Various techniques (e.g., radiosondes, global navigation satellite system, satellite remote sensing and reanalysis products by data assimilation) can be used to measure (or retrieve) PWV. However, gathering PWV data with high spatial and temporal resolutions remains a challenge. In this study, we propose a new empirical PWV grid model (called ASV-PWV) using the zenith wet delay from the Askne model and improved by the spherical harmonic function and vertical correction. Our method is convenient and enables the user to gain PWV data with only four input parameters (e.g., the longitude and latitude, time, and atmospheric pressure of the desired position). Profiles of 20 radiosonde stations in Qinghai Tibet Plateau, China, along with the latest publicly available C-PWVC2 model are used to validate the local performance. The PWV data from ASV-PWV and C-PWVC2 is generally consistent with radiosonde (the average annual bias is −0.44 mm for ASV-PWV and −1.36 mm for C-PWVC2, the root mean square error (RMSE) is 3.44 mm for ASV-PWV and 2.51 mm for C-PWVC2, respectively). Our ASV-PWV performs better than C-PWVC2 in terms of seasonal characteristics. In general, a sound consistency exists between PWV values of ASV-PWV and the fifth generation of European Centre for Medium-Range Weather Forecasts Atmospheric Reanalysis (ERA5) (total 7381 grid points in 2020). The average annual bias and RMSE are −0.73 mm and 4.28 mm, respectively. ASV-PWV has a similar performance as ERA5 reanalysis products, indicating that ASV-PWV is a potentially alternative option for rapidly gaining PWV.

The Wanshan calibration site (WSCS) is the first in-situ field for calibration and validation (Cal/Val) of HY-2 satellite series in China. It was built in December, 2018 and began business operation in 2020. In order to define an accurate datum for Cal/Val of altimeters, the permanent GNSS station (PGS) data of the WSCS observed on Zhiwan (ZWAN) and Wailingding (WLDD) islands were processed using GAMIT/GLOBK software in a regional solution, combined with 61 GNSS stations distributed nearby, collected from the GNSS Research Center, Wuhan University (GRC). The Hector software was used to analyze the trend of North (N), East (E), and Up (U) directions using six different noise models with criteria of maximum likelihood estimation (MLE), Akaike Information Criteria (AIC), and the Bayesian Information Criteria (BIC). We found that the favorite noise models were white noise plus generalized Gauss–Markov noise (WN + GGM), followed by generalized Gauss–Markov noise (GGM). Then, we compared the PGS velocities of each direction with the Scripps Orbit and Permanent Array Center (SOPAC) output parameters and found that there was good agreement between them. The PGSs in the WSCS had velocities in the N, E, and U directions of −10.20 ± 0.39 mm/year, 31.09 ± 0.36 mm/year, and −2.24 ± 0.66 mm/year for WLDD, and −10.85 ± 0.38 mm/year, 30.67 ± 0.30 mm/year, and −3.81 ± 0.66 mm/year for ZWAN, respectively. The accurate datum was defined for Cal/Val of altimeters for WSCS as a professional in-situ site. Moreover, the zenith wet delay (ZWD) of the coastal PGSs in the regional and sub-regional solutions was calculated and used to validate the microwave radiometers (MWRs) of Jason-3, Haiyang-2B (HY-2B), and Haiyang-2C (HY-2C). A sub-regional PGS solution was processed using 19 continuous operational reference stations (CORS) of Hong Kong Geodetic Survey Services to derive the ZWD and validate the MWRs of the altimeters. The ZWD of the PGSs were compared with the radiosonde-derived data in the regional and sub-regional solutions. The difference between them was −7.72~2.79 mm with an RMS of 14.53~18.62 mm, which showed good consistency between the two. Then, the PGSs’ ZWD was used to validate the MWRs. To reduce the land contamination of the MWR, we determined validation distances of 6~30 km, 16~28 km, and 18~30 km for Jason-3, HY-2B, and HY-2C, respectively. The ZWD differences between PGSs and the Jason-3, HY-2B, and HY-2C altimeters were −2.30 ± 16.13 mm, 9.22 ± 22.73 mm, and −3.02 ± 22.07 mm, respectively.

A Doppler lidar mounted on a rotary platform has been developed for measuring wind fields in the upper troposphere and stratosphere. The rotating platform was used to support a large system for the detection of wind velocities of sight (VOS) in four directions. The principle, structure, and parameters of the lidar system are introduced. The Fabry-Perot interferometer (FPI), the core component of the wind measurement system, was designed after comprehensively considering the measurement uncertainty and the influence of Mie scattering. Its dual-edge channel bandwidth is 1.05 GHz with 3.48 GHz spacing. In operation, the FPI channels are locked to the laser frequency with a stability of 14.8 MHz. Compared with the local radiosonde, it was found that the deviation in wind speed below 28 km was generally less than 10 m/s, and the deviation in wind direction below 19 km was less than 10 degrees. The 42-day profile comparison between lidar in Hefei and radiosondes in Anqing and Fuyang was analyzed. The statistical results show that the wind speed and wind direction deviations between lidar and radiosondes below 20 km were approximately 10 m/s and 20 degrees, respectively, which are comparable to the regional differences in the wind field. However, as altitudes exceed 20 km, the deviations increased rapidly with height. The experiments indicate that the Doppler lidar could measure wind fields from 7 km to 30 km, with better detection accuracy below 20 km.

The weighted mean temperature of the troposphere, Tm, is a key parameter in GNSS meteorology. It can be routinely derived based on meteorological data from radiosonde (RS) or numerical weather models. Alternatively, it can be also derived through a least-squares model of the ratio between the precipitable water vapor from RS data and the zenith wet delay estimates from GNSS measurement in the precise point positioning mode. In this last case, we found anomalous Tm values for the remote sub-tropical humid location of the Tahiti Island in the South Pacific Ocean and traced these anomalous values to anomalous zenith total delays (ZTD) that seem to have an accuracy poorer by one order of magnitude than the claimed accuracy of ZTD delays from worldwide databases. The possible causes of these discrepancies are discussed.

Although ground-based precipitable water vapor (PWV) can be obtained with a high temporal resolution and spatial resolution of tens of kilometers in an urban area using Global Navigation Satellite System (GNSS) observation, it remains fairly sparse in the vast regions over Mainland China. Satellite-derived PWV has a high spatial resolution, thereby enabling the accurate investigation of regional climate change. However, understanding the quality of satellite-derived PWV products is a prerequisite before use, which has become the focus of this study. PWV products, namely, MOD05_L2 and MYD05_L2, over the entirety of Mainland China derived from MODerate Resolution Imaging Spectroradiometer (MODIS) onboard Terra and Aqua satellites are validated over the period 2000–2017 using multiple sources, including GNSS, radiosonde, AErosol RObotic NETwork (AERONET), and European Center for Medium-Range Weather Forecasts ERA-Interim (ECMWF). The accuracy of MODIS PWV products is less than 4.3 mm over the entirety of Mainland China; however, it varies in four areas of South China (SC), North China (NC), Northwest China (NWC), and Tibetan Plateau (TP), separately with values ranging from 3 mm to 6 mm. A linear fit model is applied to calibrate the MODIS PWV products, and the accuracies of the corrected PWV from MODIS infrared (IR) and near-infrared (NIR) products have been improved by approximately 7.5% and 50.6%, respectively. The MODIS PWV is compared and calibrated in the four areas, and the improved accuracies vary widely. The root mean square error (RMSE) of IR PWV is approximately less than 4 mm over China except for the SC area with a value of approximately 5.3 mm after calibration, whereas the values of NIR PWV are approximately 2 mm over the entirety of Mainland China, except for the TP area, with a value of approximately 2.6 mm. The MODIS NIR PWV performs better than that of IR PWV data in most areas of Mainland China regardless of whether with calibration. The validation and calibration of MODIS water vapor products over Mainland China also indicate their capability to investigate the seasonal and annual variations as well as long-term trend changes in water vapor in China.

Accurately determining the atmospheric boundary layer height (ABLH) is needed when one is addressing the air quality-related issues in highly urbanized areas, as well as when one is investigating issues that are related to the emission and transport of dust aerosols over the source region. In this study, we propose a new ABLH retrieval method, which is named ADEILP (ABLH that is determined by polarization lidar); it is based on the short-term polarized lidar observation that took place during the intensive field campaign in July 2021 in Tazhong, the hinterland of Taklimakan Desert. Furthermore, we conducted comparisons between the ABLH that was identified using a radiosonde (ABLHsonde), the ABLH that was identified by ERA5 (ABLHERA5) and the ABHL that was identified by ADELIP (ABLHADELIP), and we discussed the implications of the dust events. The ADELIP method boasts remarkable advancements in two parts: (1) the lidar volume linear depolarization ratio (VLDR) that represented the aerosol type was adopted, which is very effective in distinguishing between the different types of boundary layers (e.g., mixing layer and residual layer); (2) the idea of breaking up the entire layer into sub-layers was applied on the basis of the continues wavelet transform (CWT) method, which is favorable when one is considering the effect of fine stratification in an aerosol layer. By combining the appropriate height limitations, these factors ensured that there was good robustness of the ADELIP method, thereby enabling it to deal with complex boundary layer structures. The comparisons revealed that ABLHADELIP shows good consistency with ABLHsonde and ABLHERA5 for non-dust events. Nevertheless, the ADELIP method overestimated the stable boundary layer and underestimated the heights of the mixing layer. The dust events seem to be a possible reason for the great difference between ABLHERA5 and ABLHsonde. Thus, it is worth suggesting that the influence that is caused by the differences of the vertical profile in the ERA5 product should be carefully considered when the issues on dust events are involved. Overall, these findings support the climatological analysis of the atmosphere boundary layer and the vertical distribution characteristics of aerosols over typical climatic zones.

A macro-vertical structure is closely related to weather evolution and the energy budget balance of the atmospheric system of the Earth. In this study, radiosonde data were used to identify a cloud vertical structure (CVS) using the adjusted relative humidity threshold method. To evaluate the reliability and stability of this method, the results obtained based on the spatiotemporal matching criteria established in this study were compared with Ka-band millimetre-wave cloud radar (MMCR) observation data. This comparison showed that both devices exhibit high consistency in low-level cloud detection. With the increase in the cloud height, the frequency of the cloud appearance detection by the radiosonde became higher than that by the MMCR. In spring, the results of the CVS detection by the two devices were in good agreement. Specifically, the determination coefficients of the modified degrees of freedom (adjusted R-square) of the cloud base height (CBH) and cloud top height (CTH) detected by the two devices were 0.934 and 0.879, respectively. The horizontal drift of the radiosonde was the smallest in summer, and the adj. R-square values of the CBH and CTH were 0.814 and 0.852, respectively. The CVS observation results by the radiosonde and the MMCR were significantly different in autumn (the adj. R-Square values of the CBH and CTH were 0.715 and 0.629, respectively). In winter, the adj. R-Square values of the CBH and CTH observed by the radiosonde and the MMCR were 0.958 and 0.710, respectively. The statistics and analysis of the results of the distribution characteristics of the CVSs using radiosonde data from 2019 to 2021 from Xi’an showed that the average CTH and CBH were at 7–10 km and 3–5 km, respectively. The frequencies of the cloud absence, rainfall, and two- and three-layer clouds were the highest in the winter (34.36%), autumn (12.99%), and summer, respectively.

Atmospheric water vapor is one of the important parameters for weather and climate studies. Generally, atmospheric water vapor can be monitored by some techniques, such as the Global Navigation Satellite System (GNSS), radiosonde (RS), remote sensing and numerical weather forecast (NWF). However, the comprehensive retrieval and application of precipitable water vapor (PWV) using multi techniques has been hardly performed before, which becomes the focus of this study. A comprehensive PWV retrieval and application platform (CPRAP) is first established by combing the ground-based (GNSS), space-based (Fengyun-3A, Sentinel-3A) and reanalysis-based (the fifth-generation reanalysis dataset of the European Centre for Medium-Range Weather Forecasting, ERA5) techniques. Additionally, its applications are then extended to drought and rainfall monitoring using the CPRAP-derived PWV. The statistical result shows that PWV derived from ground-based GNSS has high accuracy in China, with the root mean square (RMS), Bias and mean absolute error (MAE) of 2.15, 0.05 and 1.65 mm, respectively, when the RS-derived PWV is regarded as the reference. In addition, the accuracy of PWV derived from the space-based (FY-3A and Sentinel-3A) techniques technique is also validated and the RMS, Bias and MAE of a Medium Resolution Spectral Imager (MERSI) onboard Fengyun-3A (FY-3A) and an Ocean and Land Color Instrument (OLCI) onboard Sentinel-3A are 4.46/0.56/3.61 mm and 2.95/0.01/1.37 mm, respectively. Then, the performance of ERA5-derived PWV is evaluated based on GNSS-derived and RS-derived PWV. The result also shows good accuracy of ERA5-provided PWV with the averaged RMS, Bias and MAE of 1.86/0.11/1.48 mm and 0.90/−0.05/1.51 mm, respectively. Finally, the PWV data derived from the established CPRAP are further used for drought and rainfall monitoring. The applied results reveal that the calculated the standardized precipitation evapotranspiration index (SPEI) using the CPRAP-derived PWV can monitor the drought and the correlation coefficient ranges from 0.83 to 0.9 when compared with the SPEI. Furthermore, in this paper correlation analysis between PWV derived from the CPRAP and rainfall, and its potential for rainfall monitoring was also validated. Such results verify the significance of the established CPRAP for weather and climate studies.

Methods for now-casting adverse weather conditions with the potential to cause disruption to aircraft landings often make use of real-time measurements at high temporal resolution. This paper describes processing methodologies developed to derive meteorological parameters from such measurements recorded in the vicinity of Hong Kong International Airport, specifically a radiometer in King’s Park, a wind profiler and surface anemometer on Cheung Chau Island and weather buoys in the Pearl River estuary. These parameters are suitable for use as input to a now-casting application of the computationally efficient airflow model, FLOWSTAR, which has previously been shown to predict mountain waves generated by flow over Lantau Island to the southeast of the airport. Radiosonde data from King’s Park have been used to test the radiometer processing method; the novel approach of using minimum and maximum potential temperature deviations from a series of height-dependent linear profiles to derive radiometer inversion layer parameters generates data that compares well with values derived from corresponding radiosonde profiles. Mountain wave strength depends on the magnitude of wind speed in the inversion layer; wind profiler data can be used to estimate typical and maximum wind speeds and associated wind directions using estimates of inversion layer depth derived from the radiometer data. With estimates of surface sensible heat flux appropriate for the airport’s coastal location calculated using a marine boundary layer scheme, a dataset of meteorological parameters at 20-min resolution has been derived for input into the FLOWSTAR model. The combination of automated meteorological data processing methods and flow field modelling has the potential to form part of a now-casting system for determining strong wind shear conditions at the airport.

Water vapor is the most dominant greenhouse gas in the atmosphere and plays a critical role in Earth’s energy budget and hydrological cycle. This study aims to characterize the long-term seasonal variation of relative humidity (RH), convective available potential energy (CAPE), and convective inhibition (CIN) from surface and radiosonde observations from 1980–2020. The results show that during the monsoon season, very high RH values are depicted while low values are depicted during the pre-monsoon season. West Coast stations represent large RH values compared to other stations throughout the year. Irrespective of the season, the coastal regions show higher RH values during monsoon season. Regardless of season, the coastal regions have higher RH values during the monsoon season. During the pre-monsoon season, the coastal region has high RH values, whereas other regions have high RH values during the monsoon season. The rate of increase in RH in North-West India is 5.4%, followed by the West Coast, Central, and Southern parts of India. An increase in water vapor leads to raised temperature, which alters the instability conditions. In terms of seasonal variation, our findings show that CAPE follows a similar RH pattern. CAPE increases sharply in Central India and the West Coast region, while it declines in South India. Opposite features are observed in CIN with respect to CAPE variability over India. The results of the study provide additional evidence with respect to the role of RH as an influencing factor for an increase in CAPE over India.

In Taiwan, the frequency of afternoon convection increases in summer (July and August), and the peak hour of afternoon convection occurs at 1500–1600 local solar time (LST). Afternoon convection events are forecasted based on the atmospheric stability index, as computed from the 0800 LST radiosonde data. However, the temporal and spatial resolution and forecast precision are not satisfactory. This study used the observation data of Aqua satellite overpass near Taiwan around 1–3 h before the occurrence of afternoon convection. Its advantages are that it improves the prediction accuracy and increases the data coverage area, which means that more airports can use results of this research, especially those without radiosondes. In order to determine the availability of Atmospheric Infrared Sounder (AIRS) in Taiwan, 2010–2016 AIRS and radiosonde-sounding data were used to determine the accuracy of AIRS. This study also used 2017–2018 AIRS data to establish K index (KI) and total precipitable water (TPW) thresholds for the occurrence of afternoon convection of four airports in Taiwan. Finally, the KI and TPW were calculated using the independent AIRS atmospheric sounding (2019–2020) to forecast the occurrence of afternoon convection at each airport. The average predictive accuracy rate of the four airports is 84%. Case studies at Hualien Airport show the average predictive accuracy rate of this study is 81.8%, which is 9.1% higher than that of the traditional sounding forecast (72.7%) during the same period. Research results show that using AIRS data to predict afternoon convection in this study could not only increase data coverage area but also improve the accuracy of the prediction effectively.

The atmospheric humidity in the Polar Regions is an important factor for the global budget of water vapour, which is a significant indicator of Earth’s climate state and evolution. The Global Navigation Satellite System (GNSS) can make a valuable contribution in the calculation of the amount of Precipitable Water Vapour (PW). The PW values retrieved from Global Positioning System (GPS), hereafter PWGPS, refer to 20-year observations acquired by more than 40 GNSS geodetic stations located in the polar regions. For GNSS stations co-located with radio-sounding stations (RS), which operate Vaisala radiosondes, we estimated the PW from RS observations (PWRS). The PW values from the ERA-Interim global atmospheric reanalysis were used for validation and comparison of the results for all the selected GPS and RS stations. The correlation coefficients between times series are very high: 0.96 for RS and GPS, 0.98 for RS and ERA in the Arctic; 0.89 for RS and GPS, 0.97 for RS and ERA in Antarctica. The Root-Mean-Square of the Error (RMSE) is 0.9 mm on average for both RS vs. GPS and RS vs. ERA in the Arctic, and 0.6 mm for RS vs. GPS and 0.4 mm for RS vs. ERA in Antarctica. After validation, long-term trends, both for Arctic and Antarctic regions, were estimated using Hector scientific software. Positive PWGPS trends dominate at Arctic sites near the borders of the Atlantic Ocean. Sites located at higher latitudes show no significant values (at 1σ level). Negative PWGPS trends were observed in the Arctic region of Greenland and North America. A similar behaviour was found in the Arctic for PWRS trends. The stations in the West Antarctic sector show a general positive PWGPS trend, while the sites on the coastal area of East Antarctica exhibit some significant negative PWGPS trends, but in most cases, no significant PWRS trends were found. The present work confirms that GPS is able to provide reliable estimates of water vapour content in Arctic and Antarctic regions too, where data are sparse and not easy to collect. These preliminary results can give a valid contribution to climate change studies.

Global Navigation Satellite System (GNSS) water vapor tomography provides a four-dimensional (4-D) distribution of water vapor in the atmosphere for weather monitoring. It has developed into a widely used technique in numerical weather prediction (NWP). Vertical stratification is essential in discretizing the tomographic region. Traditional discretization methods divide the tomographic area into regular voxels with an equal height interval, which ignores the dynamic exponential distribution of water vapor. In recent years, non-uniform stratification methods have been widely validated by tomographic experiments. However, such experiments have not proposed a specific calculation method for stratification thickness. Therefore, in this paper, we introduced an adaptive non-uniform stratification method that follows the exponential distribution of water vapor in the tomographic region and presented the process of iterative calculation to acquire the optimal stratification interval. The proposed approach was applied based on the exponential decreasing trend in water vapor with increasing altitude. Moreover, it could adaptively calculate the interval of stratification height according to water vapor content. The tomographic experiments were performed using Global Positioning System (GPS) data from 19 ground-based stations in the Hong Kong Satellite Positioning Reference Station Network (SatRef) from 1 to 31 August 2019. The results indicated that, compared to the traditional stratification method, the root mean square error derived from the proposed approach was reduced by 0.26 g/m3. Additionally, severe weather can negatively affect the accuracy of the tomographic results. The results also showed that the accuracy of the tomographic results was reduced with increasing altitude. Moreover, the performance of the tomographic water vapor fields below 3000 m was improved by the proposed approach.

Precipitable water vapor (PWV) plays a vital role in climate research, especially for Antarctica in which meteorological observations are insufficient due to the adverse climate and topography therein. Reanalysis data sets provide a great opportunity for Antarctic water vapor research. This study investigates the climatological PWV means, variability and trends over Antarctica from four reanalyses, including the fifth generation of European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA5), the Second Modern-Era Retrospective analysis for Research and Applications (MERRA-2), Japanese 55-year Reanalysis (JRA-55) and National Centers for Environmental Prediction/Department of Energy (NCEP/DOE), in the period of 2001–2018 based on radiosonde and GNSS observations. PWV data from the ERA5, MERRA-2, JRA-55 and NCEP/DOE have been evaluated by radiosonde and GNSS observations, showing that ERA5 and MERRA-2 perform better than JRA-55 and NCEP/DOE with mean root mean square (RMS) errors below 1.2 mm. The climatological PWV mean distribution over Antarctica roughly shows a decreasing trend from west to east, with the highest content in summer and the lowest content in winter. The PWV variability is generally small over Antarctica, showing a seasonal dependence that is larger in the cold season and smaller in the warm season. PWV trends for all reanalyses at most Antarctic regions are insignificant and most reanalyses present overall drying trends from 2001 to 2018, except for ERA5 exhibiting a moistening trend. PWV trends also show seasonal and regional dependence. All reanalyses are generally consistent with radiosonde and GNSS observations in reproducing the PWV means (mean differences within 1.1 mm), variability (mean differences within 3%) and trends (mean differences within 6.4% decade−1) over Antarctica, except for NCEP/DOE showing spurious variability and trends in East Antarctica. Results can help us further understand these four reanalysis PWV products and promote climate research in Antarctica.

Tropospheric delay is a major error source in the Global Navigation Satellite System (GNSS), and the weighted mean temperature (Tm) is a key parameter in precipitable water vapor (PWV) retrieval. Although reanalysis products like the National Centers for Environmental Prediction (NCEP) and the European Center for Medium-Range Weather Forecasts (ECMWF) Re-Analysis-Interim (ERA-Interim) data have been used to calculate and model the tropospheric delay, Tm, and PWV, the limitations of the temporal and spatial resolutions of the reanalysis data have affected their performance. The release of the fifth-generation accurate global atmospheric reanalysis (ERA5) and the second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) provide the opportunity to overcome these limitations. The performances of the zenith tropospheric delay (ZTD), zenith wet delay (ZWD), Tm, and zenith hydrostatic delay (ZHD) of ERA5 and MERRA-2 data from 2016 to 2017 were evaluated in this work using GNSS ZTD and radiosonde data over the globe. Taking GNSS ZTD as a reference, the ZTD calculated from MERRA-2 and ERA5 pressure-level data were evaluated in temporal and spatial scales, with an annual mean bias and root mean square (RMS) of 2.3 and 10.9 mm for ERA5 and 4.5 and 13.1 mm for MERRA-2, respectively. Compared to radiosonde data, the ZHD, ZWD, and Tm derived from ERA5 and MERRA-2 data were also evaluated on temporal and spatial scales, with annual mean bias values of 1.1, 1.7 mm, and 0.14 K for ERA5 and 0.5, 4.8 mm, and –0.08 K for MERRA-2, respectively. Meanwhile, the annual mean RMS was 4.5, 10.5 mm, and 1.03 K for ERA5 and 4.4, 13.6 mm, and 1.17 K for MERRA-2, respectively. Tropospheric parameters derived from MERRA-2 and ERA5, with improved temporal and spatial resolutions, can provide a reference for GNSS positioning and PWV retrieval.

Atmospheric data obtained from the radio occultation (RO) technique are a well-recognized source of information for weather and climate studies. From the Challenging Minisatellite Payload (CHAMP) mission launched in July 2000 to the most recent Constellation Observing System for Meteorology, Ionosphere, and Climate follow-on (COSMIC-2) program, a continuous RO dataset of about 20 years has been collected, and a new opportunity for long-term climate analyses using multi-mission RO observations has subsequently arisen. Therefore, assessments of the long-term homogeneities of multi-mission RO data have become a necessary research task. For this purpose, in this study, we identified systematic discrepancies between the RO temperature profiles from the CHAMP, COSMIC, and Meteorological Operational Polar Satellite (METOP) missions. The results show that the temperature profiles from all three RO missions agree well in the upper troposphere and lower stratosphere (UTLS, 9–20 km altitude) regions, while some systematic discrepancies are found in the lower troposphere (2–8 km) and the high-altitude region (21–30 km). The homogeneities of long-term RO temperature climatologies were assessed by comparing them with radiosonde temperature records. The results of this comparison show obvious temporal inhomogeneities in the lower troposphere. The reasons for these temporal inhomogeneities include the systematic discrepancies between multi-mission RO profiles, the different monthly numbers of RO profiles, and the residual sampling error. The results of this study suggest that the systematic discrepancies between different RO missions should be thoroughly considered in the development of long-term multi-mission RO-based climatologies.

In this study, we validated the consistency of the GRUAN RS92 and RS41 datasets, versions EDT.1 and GDP.2, in the upper troposphere and lower stratosphere (200–20 hPa), through dual launch campaigns at the GRUAN site and using the radio occultation (RO) product and the ERA5 reanalysis from ECMWF as standards for double difference comparison. Separate comparisons with the references were also performed in order to trace the origin of the bias between the two instruments. Then, the performance of the GRUAN raw temperature correction algorithm was evaluated, from the aspects of day–night, the solar zenith angle, and the pressure level, for GDP.2 version products. The results show that RS92.EDT.1 has a warm bias of 0.355 K, compared to RS41.EDT.1, at 20 hPa, during daytime. This bias was found to mainly originate from RS92.EDT.1, based on the separate comparison with RO or ECMWF ERA5 data. RS92.GDP.2 is consistent with RS41.GDP.2, but a separate comparison indicated that the two original GDP.2 products have a ~1 K warm bias at 20 hPa during daytime, compared with RO or ECMWF ERA5 data. The GRUAN correction method can reduce the warm bias up to 0.5 K at 20 hPa during daytime. As a result, this GRUAN correction method is efficient, and it is dependent on the solar zenith angle and pressure level.

The Global Navigation Satellite System (GNSS) plays an important role in retrieving high temporal–spatial resolution precipitable water vapor (PWV) and its applications. The weighted mean temperature (Tm) is a key parameter for the GNSS PWV estimation, which acts as the conversion factor from the zenith wet delay (ZWD) to the PWV. The Tm is determined by the air pressure and water vapor pressure, while it is not available nearby most GNSS stations. The empirical formular is often applied for the GNSS station surface temperature (Ts) but has a lower accuracy. In this paper, the temporal and spatial distribution characteristics of the coefficients of the linear Tm-Ts model are analyzed, and then a piecewise-linear Tm-Ts relationship is established for each GPS station using radiosonde data collected from 2011 to 2019. The Tm accuracy was increased by more than 10% and 20% for 86 and 52 radiosonde stations, respectively. The PWV time series at 377 GNSS stations from the infrastructure construction of national geodetic datum modernization and Crustal Movement Observation Network of China (CMONC) were further obtained from the GPS observations and meteorological data from 2011 to 2019. The PWV accuracy was improved when compared with the Bevis model. Furthermore, the daily and monthly average values, long-term trend, and its change characteristics of the PWV were analyzed using the high-precision inversion model. The results showed that the averaged PWV was higher in Central-Eastern China and Southern China and lower in Northwest China, Northeast China, and North China. The PWV is increasing in most parts of China, while the some PWVs in North China show a downward trend.

From late 2019 to early 2020, forest fires in southeastern Australia caused huge economic losses and huge environmental pollution. Monitoring forest fires has become increasingly important. A new method of fire detection using the difference between global navigation satellite system (GNSS)-derived precipitable water vapor and radiosonde-derived precipitable water vapor (ΔPWV) is proposed. To study the feasibility of the new method, the relationship is studied between particulate matter 10 (PM10) (2.5 to 10 microns particulate matter) and ΔPWV based on Global Positioning System (GPS) data, radiosonde data, and PM10 data from 1 June 2019 to 1 June 2020 in southeastern Australia. The results show that before the forest fire, ΔPWV and PM10 were smaller and less fluctuating. When the forest fire happened, ΔPWV and PM10 were increasing. Then after the forest fire, PM10 became small with relatively smooth fluctuations, but ΔPWV was larger and more fluctuating. Correlation between the 15-day moving standard deviation (STD) time series of ΔPWV and PM10 after the fire was significantly higher than that before the fire. This study shows that ΔPWV is effective in monitoring forest fires based on GNSS technique before and during forest fires in climates with more uniform precipitation, and using ΔPWV to detect forest fires based on GNSS needs to be further investigated in climates with more precipitation and severe climate change.

Global Positioning System (GPS) radio occultation (RO) and radiosonde (RS) observations are two major types of observations assimilated in numerical weather prediction (NWP) systems. Observation error variances are required input that determines the weightings given to observations in data assimilation. This study estimates the error variances of global GPS RO refractivity and bending angle and RS temperature and humidity observations at 521 selected RS stations using the three-cornered hat method with additional ERA-Interim reanalysis and Global Forecast System forecast data available from 1 January 2016 to 31 August 2019. The global distributions, of both RO and RS observation error variances, are analyzed in terms of vertical and latitudinal variations. Error variances of RO refractivity and bending angle and RS specific humidity in the lower troposphere, such as at 850 hPa (3.5 km impact height for the bending angle), all increase with decreasing latitude. The error variances of RO refractivity and bending angle and RS specific humidity can reach about 30 N-unit2, 3 × 10−6 rad2, and 2 (g kg−1)2, respectively. There is also a good symmetry of the error variances of both RO refractivity and bending angle with respect to the equator between the Northern and Southern Hemispheres at all vertical levels. In this study, we provide the mean error variances of refractivity and bending angle in every 5°-latitude band between the equator and 60°N, as well as every interval of 10 hPa pressure or 0.2 km impact height. The RS temperature error variance distribution differs from those of refractivity, bending angle, and humidity, which, at low latitudes, are smaller (less than 1 K2) than those in the midlatitudes (more than 3 K2). In the midlatitudes, the RS temperature error variances in North America are larger than those in East Asia and Europe, which may arise from different radiosonde types among the above three regions.

In this paper, we assess, in the framework of Global Navigation Satellite System (GNSS) meteorology, the accuracy of GNSS propagation delays corresponding to the Saastamoinen zenith hydrostatic delay (ZHD) model and the Vienna Mapping function VMF1/VMF3 (hydrostatic and wet), with reference to radiosonde ray-tracing delays over a three-year period on 28 globally distributed sites. The results show that the Saastamoinen ZHD estimates have a mean root mean square (RMS) error of 1.7 mm with respect to the radiosonde. We also detected some seasonal signatures in these Saastamoinen ZHD estimates. This indicates that the Saastamoinen model, based on the hydrostatic assumption and the ground pressure, is insufficient to capture the full variability of the ZHD estimates over time with the accuracy needed for GNSS meteorology. Furthermore, we found that VMF3 slant hydrostatic delay (SHD) estimates outperform the corresponding VMF1 SHD estimates (equivalent SHD RMS error of 4.8 mm for VMF3 versus 7.1 mm for VMF1 at 5° elevation angle), with respect to the radiosonde SHD estimates. Unexpectedly, the situation is opposite for the VMF3 slant wet delay (SWD) estimates compared to VMF1 SWD estimates (equivalent SWD RMS error of 11.4 mm for VMF3 versus 7.0 mm for VMF1 at 5° elevation angle). Our general conclusion is that the joint approach using ZHD models and mapping functions must be revisited, at least in the framework of GNSS meteorology.

A new concept is proposed for estimating the zenith wet delay (ZWD) and atmospheric weighted average temperature by inputting the temperature, total pressure, and specific humidity from surface weather data. In addition, a new ZWD integral method is described for highly accurate calculation of the ZWD from radiosonde observation. To evaluate the advantages of the new discrete integral formula, we utilized the 8-year radiosonde profiles of 85 stations in China from 2010 to 2017 to validate the accuracy of the radiosonde-derived ZWD. The results showed that the mean accuracy of the ZWD derived from radiosonde data was 4.28 mm. Next, the new ZWD model was assessed using two sets of reference values derived from radiosonde data and GNSS precise point positioning in China. The results confirm that the new development improved the accuracy of the estimation of the tropospheric wet delay from the surface meteorological data. The performance of this new model can be seen as an important step toward accurately correcting the tropospheric delay in Global Navigation Satellite System (GNSS) real-time navigation and positioning. It can also be used in GNSS meteorology for weather forecasting and climate research.

Global navigation satellite system (GNSS) tomography can effectively sense the three-dimensional structure of tropospheric water vapor (WV) using the GNSS observations. Numerous studies have utilized a tomographic window to include more epochs of observations, which significantly increases the number of valid signals. However, considering the tomography grid limits, a massive number of valid signals inevitably exhibits linear dependence. This dependence makes it impossible to improve the rank score of the tomography coefficient matrix by blindly introducing a large number of valid rays. Furthermore, excessive valid signals may lead to a high condition number in the coefficient matrix (ill-condition problem), which causes unstable results using the GNSS-WV tomography. Considering these problems, we proposed an improved tomographic approach, which applies a refined tomographic window. It differs from the general tomographic window in that the window is refined to traverse the valid signals available 15 min before and after the target epoch while retaining only the linearly independent parts (characteristic signal). Compared to the conventional method, the proposed method can filter the characteristic signal, which increases the rank score of the coefficient matrix and improves the stability of the tomography model. In this paper, we used GNSS observations from the Hong Kong Satellite Positioning Reference Station Network (SatRef) to validate the performance of the proposed method over the day-of-year (DOY) periods of 130–132, 2019 and 146–148, 2019. The numerical results showed that, by using a refined tomographic window, the proposed method obtained superior WV products in comparison with that of the conventional method.

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

The vertical profiles of temperature and water vapour from the Atmospheric InfraRed Sounder (AIRS) have been validated across various regions of the globe as an effort to provide a substitute for radiosonde observations. However, there is a paucity of inter-comparisons over West Africa where local convective processes dominate and radiosonde observations (RAOBs) are limited. This study validates AIRS temperature and relative humidity profiles for selected radiosonde stations in West Africa. Radiosonde data were obtained from the AMMA and DACCIWA campaigns which spanned 2006–2008 and June–July 2016 respectively and offered a period of prolonged radiosonde observations in West Africa. AIRS performance was evaluated with the bias and root mean square difference (RMSD) at seven RAOB stations which were grouped into coastal and inland. Evaluation was performed on diurnal and seasonal timescales, cloud screening conditions and derived thunderstorm instability indices. At all timescales, the temperature RMSD was higher than the AIRS accuracy mission goal of ±1 K. Relative humidity RMSD was satisfactory with deviations

The Global Navigation Satellite System (GNSS) Radio Occultation (RO) retrieved temperature and specific humidity profiles can be widely used for weather and climate studies in troposphere. However, some aspects, such as the influences of background data on these retrieved moist profiles have not been discussed yet. This research evaluates RO retrieved temperature and specific humidity profiles from Wegener Center for Climate and Global Change (WEGC), Radio Occultation Meteorology Satellite Application Facility (ROM SAF) and University Corporation for Atmospheric Research (UCAR) Boulder RO processing centers by comparing with measurements from 10 selected Integrated Global Radiosonde Archive (IGRA) radiosonde stations in different latitudinal bands over 2007 to 2010. The background profiles used for producing their moist profiles are also compared with radiosonde. We found that RO retrieved temperature profiles from all centers agree well with radiosonde. Mean differences at polar, mid-latitudinal and tropical stations are varying within ±0.2 K, ±0.5 K and from −1 to 0.2 K, respectively, with standard deviations varying from 1 to 2 K for most pressure levels. The differences between RO retrieved and their background temperature profiles for WEGC are varying within ±0.5 K at altitudes above 300 hPa, and the differences for ROM SAF are within ±0.2 K, and that for UCAR are within 0.5 K at altitudes below 300 hPa. Both RO retrieved and background specific humidity above 600 hPa are found to have large positive differences (up to 40%) against most radiosonde measurements. Discrepancies of moist profiles among the three centers are overall minor at altitudes above 300 hPa for temperature and at altitudes above 700 hPa for specific humidity. Specific humidity standard deviations are largest at tropical stations in June July August months. It is expected that the outcome of this research can help readers to understand the characteristics of moist products among centers.

The atmospheric boundary layer height is important for constraining air pollution and meteorological models. This study attempted to validate the MODIS-estimated atmospheric boundary layer height (ABLH), and variation in the ABLH in Uganda was evaluated. The ABLH was estimated from MODIS data using the mixing ratio profile gradient method and compared to the ABLH estimated from radiosonde data using three different methods. Unlike in studies in other regions of the world, correlations between ABLH estimated using MODIS and radiosonde data were weak, implying limited usefulness of MODIS data for determining ABLH. However, the diurnal variation in MODIS-derived ABLH and particulate matter (PM10) was consistent with the expected inverse relationship between PM10 mass concentration and ABLH, and the mean MODIS-derived ABLH values were significantly lower during wet seasons than dry seasons, as expected.

We investigated the accuracy of operational medium-range ensemble forecasts for 29 Atlantic hurricanes between 2007 and 2019. Upper-level troughs with strong wind promoted northward movement of hurricanes over the mid-latitudes. For hurricanes with upper-level troughs, relatively large errors in the prediction of troughs result in large ensemble spreads, which result in failure to forecast hurricane track. In contrast, for hurricanes without upper-level troughs, mean central position errors are relatively small in all operational forecasts because of the absence of upper-level strong wind around troughs over the mid-latitudes. Hurricane Irma in September 2017 was accompanied by upper-level strong wind around a trough; errors and ensemble spreads for the predicted upper-level trough are small, contributing to smaller errors and small ensemble spreads in the predicted tracks of Irma. Our observing system experiment reveals that inclusion of additional Arctic radiosonde observation data obtained from research vessel Mirai in 2017 improves error and ensemble spread in upper-level trough with strong wind at initial time for forecast, increasing the accuracy of the forecast of the track of Irma in 2017.

Conventional radiosondes can be used to measure the relative humidity over liquid (RHL) by assuming a saturated vapor pressure over the liquid. However, this assumption results in significant errors with respect to measurements in the upper troposphere, where the effect of ice is dominant. Therefore, this study presents a novel method that considers the effects of ice to determine the relative humidity from radiosonde RHL data for the last 40 years (1979–2018) over the upper layers of the Korean peninsula. Even though the relative humidity obtained from the reanalysis data was significantly different from the radiosonde-based RHL, the difference was much reduced when relative humidity was calculated using the novel method proposed in this study. Such improvements in the estimated relative humidity could be attributed to the consideration of the ice effect at temperatures above freezing level. Additionally, the validity of the relative humidity estimated in this study was established based on a two-week case analysis of data from Boseong station. Furthermore, two peak relative humidity modes for the lower and upper layers were clearly identified in the mean climatology profiles, which indirectly suggested the absence of mid-level clouds around the 700-hPa level and 500-hPa level in winter and summer, respectively. This study is meaningful as it is the first study to determine the relative humidity distribution over the Korean peninsula using radiosonde observations. The scientific value obtained can potentially be expanded by applying the proposed method to other radiosonde observation networks, which are widely distributed worldwide.

The planetary boundary layer height (PBLH) is a vital parameter to characterize the surface convection, which determines the diffusion of air pollutants. The accurate inversion of PBLH is extremely important for the study of aerosol concentrations, in order to predict air quality and provide weather forecast. Aerosol lidar, a powerful remote sensing instrument for detecting the characteristics of atmospheric temporal and spatial evolution, can continuously retrieve the planetary boundary layer (PBL) and obtain high resolution measurements. However, multi-layer conditions, including one or more layers of aerosol, or cloud above the PBL, can seriously interfere the accuracy of PBLH determined by lidar. A new technique of maximum limited height initialization and range restriction (MLHI-RR) is proposed to eliminate the impact of multi-layer conditions on PBLH determination. Four widely used methods for deriving PBLH are utilized, in addition to the MLHI-RR constraint. Comparisons demonstrate that the proposed technique can determine the PBLH in multi-layer conditions with higher accuracy. The proposed technique requires no affiliate information besides lidar signals, which provide a convenient method for PBLH determination under complicated conditions.

The Global Navigation Satellite System (GNSS) tomographic technique can be used for remote sensing of the three-dimensional water vapor (WV) distribution in the troposphere, which has attracted considerable interest. However, a significant problem in this technique is the excessive reliance on constraints (particularly in large GNSS networks). In this paper, we propose an improved tomographic method based on optimized voxel, which only considers the voxels passed by GNSS rays. The proposed method can completely prevent the tomographic algorithm interference of constraints that originated from empirical functions. Experiments in Nanjing in the periods of day-of-year (DOY) 182–184, 2019, and 244–246, 2019, show that the mean absolute error (MAE) and root mean square error (RMSE) of the WV density profile obtained using the proposed method are 0.9 and 1.3 g/m3, while those obtained using the conventional method are 1.3 and 1.8 g/m3, respectively, with respect to the radiosonde (RS) method. The numerical results show that the proposed method is reliable and has a superior accuracy to that of the conventional method.

We estimate precipitable water vapor (PWV) from data collected by the low-cost Global Navigation Satellite System (GNSS) receiver at a vessel. The dual-frequency GNSS receiver that the vessel ISABU is equipped with that is operated by the Korea Institute of Ocean Science and Technology. The ISABU served in the Pacific Ocean for scientific research during a period from August 30 to September 21, 2018. It also performs radiosonde observations to obtain a vertical profile of troposphere on the vessel’s path. The GNSS-derived PWV is compared to radiosonde observations and the Atmospheric Infrared Sounder (AIRS) on NASA’s Aqua satellite output. A bias and root-mean-square (RMS) error between shipborne GNSS-PWV and radiosonde-PWV were −1.48 and 5.22 mm, respectively. When compared to the ground GNSS-PWV, shipborne GNSS-PWV has a relatively large RMS error in comparison with radiosonde-PWV. However, the GNSS observations on the vessel are still in good agreement with radiosonde observations. On the other hand, the GNSS-PWV is not well linearly correlated with AIRS-PWV. The RMS error between the two observations was approximately 8.97 mm. In addition, we showed that the vessel on the sea surface has significantly larger carrier phase multipath error compared to the ground-based GNSS observations. This also can result in reducing the accuracy of shipborne GNSS-PWV. However, we suggest that the shipborne GNSS has sufficient potential to derive PWV with the kinematic precise point positioning (PPP) solution on the vessel.

Global Navigation Satellite System (GNSS) tomography has developed into an efficient tool for sensing the high spatiotemporal variability of atmospheric water vapor. The integration of GNSS top signals and side rays for tropospheric tomography systems using a novel height factor model (HFM) is proposed and discussed in this paper. Within the HFM, the sectional slant wet delay (SWD) of inside signals (the part of the side signal inside the tomography area), which is considered a key factor for modeling side rays, is separated into isotropic and anisotropic components. Correspondingly, two height factors are defined to calculate the isotropic and anisotropic part of tropospheric delays in the HFM. In addition, the dynamic tomography top boundary is first analyzed and determined based on 30-year radiosonde data to reasonably divide signals into top and side rays. Four special experimental schemes based on different tomography regions of Hong Kong are performed to assess the proposed HFM method, the results of which show increases of 33.42% in the mean utilization of rays, as well as decreases of 0.46 g/m3 in the average root mean square error (RMSE), compared to the traditional approach, revealing the improvement of tomography solutions when the side signals are included in the modeling. Furthermore, compared with the existing correction model for modeling side rays, the water vapor profiles retrieved from the proposed improved model are closer to the radiosonde data, which highlights the advantages of the proposed HFM for optimizing the GNSS tomography model.

The turbulent mixing and dispersion of air pollutants is strongly dependent on the vertical structure of the wind, which constitutes one of the major challenges affecting the determination of boundary layer height (BLH). Here, an adaptive method is proposed to estimate BLH from measurements of radar wind profilers (RWPs) in Beijing (BJ), Nanjing (NJ), Chongqing (CQ), and Wulumuqi (WQ), China, during the summer of 2019. Validation against simultaneous BLH estimates from radiosondes (RSs) yielded a correlation coefficient of 0.66, indicating that the method can be used to derive BLH from RWPs. Diurnal variations of BLH and the ventilation coefficient (VC) at four sites were then examined. A distinct diurnal cycle of BLH was observed over all four cities; BLH gradually increased from sunset, reached a maximum in the afternoon, and then dropped sharply after sunset. The maximum hourly average BLH (1.426 ± 0.46 km) occurred in WQ, consistent with the maximum hourly mean VC larger than 5000 m2/s observed there. By comparison, the diurnal variation of VC was not strong, with values ranging between 2000 and 3000 m2/s, likely owing to the high-humidity environment. Furthermore, surface sensible heat flux, latent heat flux, and dry mass of particulate matter with aerodynamic diameter ≤2.5 µm (PM2.5) concentrations were found to somehow affect the vertical structure of wind and thermodynamic features, leading to a difference between RS and RWP BLH estimates. This indicates that the atmospheric environment can affect BLH estimates using RWP data. The BLH results from RWPs were better in some specific cases. These findings show great potential of RWP measurements in air quality research, and will provide key data references for policy-making toward emission reductions.

In this study, sensors mounted on a new type of round-trip airborne flat-floating sounding system (RTAFSS) were used to obtain the observation data of the three stages of "rising, flat-floating and falling". This new sounding method has a good application prospect. We performed spectrum analysis on the normalized temperature fluctuation, and the vertical wavenumber spectrum from the rising and falling stages and the horizontal wavenumber spectrum from the flat-floating stage were obtained. This is the first time the complete gravity wave spectrum characteristics were obtained from three consecutive stages: rising, flat-floating and falling. The results show that the gravity wave spectrum of the three stages can be well obtained by RTAFSS. For the horizontal wavenumber spectrum, the spectral slope is basically around −2, and the difference in the spectral structure of the horizontal wave number spectrum may be due to the intermittent turbulent activity and the variable intensity of the gravitational wave during its propagation. This study aims to make experimental exploration of the spectrum characteristics of gravity waves by this new type of observation data. It is expected to reveal the spectrum characteristics of horizontal wavenumber in the stratosphere region of China, providing a theoretical basis for spectrum analysis in a wider space–time range after further network observation of RTAFSS.

Clouds are significant in the global radiation budget, atmospheric circulation, and hydrological cycle. However, knowledge regarding the observed climatology of the cloud vertical structure (CVS) over Beijing is still poor. Based on high-resolution radiosonde observations at Beijing Nanjiao Weather Observatory (BNWO) during the period 2010–2017, the method for identifying CVS depending on height-resolved relative humidity thresholds is improved, and CVS estimation by radiosonde is compared with observations by millimeter-wave cloud radar and ceilometer at the same site. Good consistency is shown between the three instruments. Then, the CVS climatology, including the frequency distribution and seasonal variation, is investigated. Overall, the occurrence frequency (OF) of cloudy cases in Beijing is slightly higher than that of clear-sky cases, and the cloud OF is highest in summer and lowest in winter. Single-layer clouds and middle-level clouds are dominant in Beijing. In addition, the average cloud top height (CTH), cloud base height (CBH), and cloud thickness in Beijing are 6.2 km, 4.0 km, and 2.2 km, respectively, and show the trend of reaching peaks in spring and minimums in winter. In terms of frequency distribution, the CTH basically resides below an altitude of 16 km, and approximately 43% of the CBHs are located at altitudes of 0.5–1.5 km. The cloud OF has only one peak located at altitudes of 4–8 km in spring, whereas it shows a trimodal distribution in other seasons. The height at which the cloud OF reaches its peak is highest in summer and lowest in winter. To the best of our knowledge, the cloud properties analyzed here are the first to elucidate the distribution and temporal variation of the CVS in Beijing from a long-term sounding perspective, and these results will provide a scientific observation basis for improving the atmospheric circulation model, as well as comparisons and verifications for measurements by ground-based remote sensing equipment.

The planetary boundary layer is the medium of energy, moisture, momentum and pollutant exchange between the surface and the atmosphere. In this paper, a method to derive the boundary layer mixing height (MH) was introduced and applied over the Heihe river basin. Atmospheric profiles from the MODerate Resolution Imaging Sepctroradiometer (MODIS) instrument onboard the NASA-Aqua satellite were used for the high spatial resolution of this method. A gap-filling method was used to replace missing MODIS data. In situ MH data were also calculated from HIWATER (Heihe Watershed Allied Telemetry Experimental Research) and WATER (Watershed Allied Telemetry Experimental Research) observational radiosonde sounding data from 2008 and 2012 using the Richardson number method combined with a subjective method. The MH occurs where there is an abrupt decrease in the MR (water vapor mixing ratio). The minimum vertical gradient of the MR is used to determine the MH. The method has an average RMSE of 370 m under clear skies and convective conditions. The seasonal variation in the MH at the Gaoya radiosonde station is also presented. The study demonstrates that remote sensing methodologies can successfully estimate the MH without the help of field measurements.