Spatial interpolation of regional PM2.5 concentrations in China during COVID-19 incorporating multivariate data

During specific periods when the PM_2.5 variation pattern is unusual, such as during the coronavirus disease 2019 (COVID-19) outbreak, epidemic PM_2.5 regional interpolation models have been relatively little investigated, and little consideration has been given to the residuals of optimized models and changes in model interpolation accuracy for the PM_2.5 concentration under the influence of epidemic phenomena. Therefore, this paper mainly introduces four interpolation methods (kriging, empirical Bayesian kriging, tensor spline function and complete regular spline function), constructs geographically weighted regression (GWR) models of the PM_2.5 concentration in Chinese regions for the periods from January-June 2019 and January-June 2020 by considering multiple factors, and optimizes the GWR regression residuals using these four interpolation methods, thus achieving the purpose of enhancing the model accuracy. The PM_2.5 concentrations in many regions of China showed a downward trend during the same period before and after the COVID-19 outbreak. Atmospheric pollutants, meteorological factors, elevation, zenith wet delay (ZWD), normalized difference vegetation index (NDVI) and population maintained a certain relationship with the PM_2.5 concentration in terms of linear spatial relationships, which could explain why the PM_2.5 concentration changed to a certain extent. By evaluating the model accuracy from two perspectives, i.e., the overall interpolation effect and the validation set interpolation effect, the results showed that all four interpolation methods could improve the numerical accuracy of GWR to different degrees, among which the tensor spline function and the fully regular spline function achieved the most stable effect on the correction of GWR residuals, followed by kriging and empirical Bayesian kriging.

Significant increase in water vapour over India and Indian Ocean: Implications for tropospheric warming and regional climate forcing

The increase in greenhouse gases (GHGs) due to anthropogenic activities enhances regional and global temperatures. The most abundant GHG, i.e., water vapour, has a vital positive feedback on the global warming and Earth's climate system. This study focuses on the spatial and temporal changes in water vapour in the troposphere over India and Indian Ocean as derived from the ground-based, satellite and reanalyses data, and assesses the impact on water vapour changes on the regional climate by analysing radiative effects. The analyses show that the annual mean column water vapour (CWV) is high over the northern Indian Ocean, Bay of Bengal and Peninsular India, ranging from 30 to 60 kg/m². Most regions show significant positive trends in the annual mean CWV, about 0.1-0.2 kg/m²/yr. There is a significant positive trend in water vapour in the troposphere (except 200 hPa) over the India land regions, with the highest values at 1000 hPa (0.034 g/kg/yr). The corresponding water vapour radiative effect (WVRE) is about 20-80 W/m², depending on seasons and regions. This study, therefore, indicates that the increase in tropospheric water vapour over India and Indian Ocean could affect the regional temperature and climate.

Evaluation of Water Vapor Product from TROPOMI and GOME-2 Satellites against Ground-Based GNSS Data over Europe

A novel integrated water vapor (IWV) product from TROPOspheric Monitoring Instrument (TROPOMI) is validated together with a Global Ozone Monitoring Instrument-2 (GOME-2) standard product. As reference, ground-based Global Navigation Satellite Systems (GNSS) IWV data in 235 European stations from May 2018 to May 2019 are used. Under cloud free situations, a general comparison is carried out. It suggests that TROPOMI IWV exhibits less bias than GOME-2 and better results in the dispersion and regression parameters. Moreover, TROPOMI presents more homogeneous results along the different stations. However, TROPOMI is found to be overestimating the IWV uncertainties and being, therefore, too conservative in the confidence interval considered. The dependence of satellite product performance on several variables is also discussed. TROPOMI IWV shows wet bias of 5.7% or less for IWV < 10 mm (TROPOMI) and dry bias of up to −3% (TROPOMI). In contrast, GOME-2 shows wet bias of 30% or less for IWV < 25 mm (GOME-2) and dry bias of −12.3% for IWV > 25 mm. In addition, relative standard deviation (rSD) increases as IWV increases. In addition, the dependence on solar zenith angle (SZA) was also analyzed, as solar radiation bands are used in the retrieval algorithm of both instruments. Relative mean bias error (rMBE) shows positive values for GOME-2, slightly increasing with SZA, while TROPOMI shows more stable values. However, under high SZA, GOME-2 IWV exhibits a steep increase in rMBE (overestimation), while TROPOMI IWV exhibits a moderate decrease (underestimation). rSD is slightly increasing with SZA. The influence of cloudiness on satellite IWV observations is such that TROPOMI tends to overestimate IWV more as cloudiness increases, especially for high IWV. In the case of GOME-2, the rSD slightly increases with cloudiness, but TROPOMI rSD has a marked increase with increasing cloudiness. TROPOMI IWV is an important source of information about moisture, but its algorithm could still benefit from further improvement to respond better to cloudy situations.

Global Spatiotemporal Variability of Integrated Water Vapor Derived from GPS, GOME/SCIAMACHY and ERA-Interim: Annual Cycle, Frequency Distribution and Linear Trends

Atmospheric water vapor plays a prominent role in climate change and atmospheric, meteorological, and hydrological processes. Because of its high spatiotemporal variability, precise quantification of water vapor is challenging. This study investigates Integrated Water Vapor (IWV) variability for the period 1995–2010 at 118 globally distributed Global Positioning System (GPS) sites, using additional UV/VIS satellite retrievals by GOME, SCIAMACHY, and GOME-2 (denoted as GOMESCIA below), plus ERA-Interim reanalysis output. Apart from spatial representativeness differences, particularly at coastal and island sites, all three IWV datasets correlate well with the lowest mean correlation coefficient of 0.878 (averaged over all the sites) between GPS and GOMESCIA. We confirm the dominance of standard lognormal distribution of the IWV time series, which can be explained by the combination of a lower mode (dry season characterized by a standard lognormal distribution with a low median value) and an upper mode (wet season characterized by a reverse lognormal distribution with high median value) in European, Western American, and subtropical sites. Despite the relatively short length of the time series, we found a good consistency in the sign of the continental IWV trends, not only between the different datasets, but also compared to temperature and precipitation trends.

FENGYUN-4A Advanced Geosynchronous Radiation Imager Layered Precipitable Water Vapor Products’ Comprehensive Evaluation Based on Quality Control System

A physical retrieval algorithm has been developed for deriving the layered precipitable water vapor (LPWs) product from infrared radiances of the Advanced Geosynchronous Radiation Imager (AGRI) onboard FengYun-4A (FY-4A), the first of the new generation of Chinese geostationary weather satellites (FengYun-4, or FY-4 Series). The FY-4A AGRI LPWs are evaluated with different types of reference datasets based on Quality Control System (QCS), including those from Himawari-8 AHI (Advanced Himawari Imager), MODIS (Moderate Resolution Imaging Spectroradiometer), Radiosonde, ERA5 (European Centre for Medium-Range Weather Forecasts Reanalysis v5), NCEP (National Centers for Environmental Prediction) reanalysis and CMA (China Meteorological Administration) forecast product from global medium range numerical weather prediction (NWP) system. QCS is one of the important components of FY-4A ground segment, which mainly focuses on the satellite products’ evaluation and validation. It is found that the AGRI LPW product has a good agreement with different evaluating sources and the quality is favorable and stable. With the capability of frequent (5-min interval) observations over the East Asia and Western Pacific regions, the AGRI LPW products can be used to investigate the atmospheric temporal and spatial variations in the pre-landfall environment for typhoons. The QCS is a useful tool to monitor, evaluate, and validate the AGRI LPW products.

Understanding the Present-Day Spatiotemporal Variability of Precipitable Water Vapor over Ethiopia: A Comparative Study between ERA5 and GPS

Atmospheric water vapor plays a crucial role in atmospheric, climate change, meteorological, and hydrological processes. In a country like Ethiopia, with its complex topography and synoptic-scale spatiotemporal circulation patterns, the analysis of the spatiotemporal variability of precipitable water vapor (PWV) is very challenging, and is hampered by the lack of long observational datasets. In this study, we process the PWV over eight Ethiopian global positioning system (GPS) sites and one close to the Ethiopian eastern border, for the available common period 2013–2020, and compare with the PWV retrieved from the state-of-the-art ERA5 reanalysis. Both PWV datasets agree very well at our sample, with correlation coefficients between 0.96 and 0.99, GPS-PWV show a moderate wet bias compared to ERA5-PWV for the majority of the sites, and an overall root mean square error of 3.4 mm. Seasonal and diurnal cycles are also well captured by these datasets. The seasonal variations of PWV and precipitation at the sites agree very well. Maximum diurnal PWV amplitudes are observed for stations near water bodies or dense vegetation, such as Arbaminch (ARMI) and Bahir Dar (BDMT). At those stations, the PWV behavior at heavy rainfall events has been investigated and an average 25% increase (resp. decrease) from 12 h before (resp. 12 h after) the start of the rainfall event, when the PWV peaks, has been observed.

Rainfall Forecast Model Based on the TabNet Model

To further reduce the error rate of rainfall prediction, we used a new machine learning model for rainfall prediction and new feature engineering methods, and combined the satellite system’s method of observing rainfall with the machine learning prediction. Based on multivariate correlations among meteorological information, this study proposes a rainfall forecast model based on the Attentive Interpretable Tabular Learning neural network (TabNet). This study used self-supervised learning to help the TabNet model speed up convergence and maintain stability. We also used feature engineering methods to alleviate the uncertainty caused by seasonal changes in rainfall forecasts. The experiment used 5 years of meteorological data from 26 stations in the Beijing–Tianjin–Hebei region of China to verify the proposed rainfall forecast model. The comparative experiment proved that our proposed method improves the performance of the model, and that the basic model used is also superior to other traditional models. This research provides a high-performance method for rainfall prediction and provides a reference for similar data-mining tasks.

High-Precision GNSS PWV and Its Variation Characteristics in China Based on Individual Station Meteorological Data

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.

Retrieval of Daytime Total Column Water Vapour from OLCI Measurements over Land Surfaces

A new retrieval of total column water vapour (TCWV) from daytime measurements over land of the Ocean and Land Colour Instrument (OLCI) on-board the Copernicus Sentinel-3 missions is presented. The Copernicus Sentinel-3 OLCI Water Vapour product (COWa) retrieval algorithm is based on the differential absorption technique, relating TCWV to the radiance ratio of non-absorbing band and nearby water vapour absorbing band and was previously also successfully applied to other passive imagers Medium Resolution Imaging Spectrometer (MERIS) and Moderate Resolution Imaging Spectroradiometer (MODIS). One of the main advantages of the OLCI instrument regarding improved TCWV retrievals lies in the use of more than one absorbing band. Furthermore, the COWa retrieval algorithm is based on the full Optimal Estimation (OE) method, providing pixel-based uncertainty estimates, and transferable to other Near-Infrared (NIR) based TCWV observations. Three independent global TCWV data sets, i.e., Aerosol Robotic Network (AERONET), Atmospheric Radiation Measurement (ARM) and U.S. SuomiNet, and a German Global Navigation Satellite System (GNSS) TCWV data set, all obtained from ground-based observations, serve as reference data sets for the validation. Comparisons show an overall good agreement, with absolute biases between 0.07 and 1.31 kg/m2 and root mean square errors (RMSE) between 1.35 and 3.26 kg/m2. This is a clear improvement in comparison to the operational OLCI TCWV Level 2 product, for which the bias and RMSEs range between 1.10 and 2.55 kg/m2 and 2.08 and 3.70 kg/m2, respectively. A first evaluation of pixel-based uncertainties indicates good estimated uncertainties for lower retrieval errors, while the uncertainties seem to be overestimated for higher retrieval errors.

A Regional NWP Tropospheric Delay Inversion Method Based on a General Regression Neural Network Model

Tropospheric delay is a major error source that affects the initialization and re-initialization speed of the Global Navigation Satellite System’s (GNSS) medium-/long-range baseline in Network Real-Time Kinematic (NRTK) positioning. Fusing the meteorological data from the Numerical Weather Prediction (NWP) model to estimate the zenith tropospheric delay (ZTD) is one of the current research hotspots. However, research has shown that the ZTD derived from NWP models is still not accurate enough for high-precision GNSS positioning applications without the estimation of the residual tropospheric delay. To date, General Regression Neural Network (GRNN) has been applied in many fields. It has a high learning speed and simple structure, and can approximate any function with arbitrary precision. In this study, we developed a regional NWP tropospheric delay inversion method based on a GRNN model to improve the accuracy of the tropospheric delay derived from the NWP model. The accuracy of the tropospheric delays derived from reanalysis data of the European Center for Medium-Range Weather Forecasts (ECMWF) and the US National Centers for Environmental Prediction (NCEP) was assessed through comparisons with the results of the International GPS Service (IGS). The variation characteristics of the residual of the ZTD inverted by NWP data were analyzed considering the factors of temperature, humidity, latitude, and season. To evaluate the performance of this new method, the National Center Atmospheric Research (NCAR) troposphere data of 650 stations in Japan in 2005 were collected as a reference to compare the accuracy of the ZTD before and after using the new method. The experimental results showed that the GRNN model has obvious advantages in fitting the NWP ZTD residual. The mean residual and the root mean square deviation (RMSD) of the ZTD inverted using the method of this study were 9.5 mm and 12.7 mm, respectively, showing reductions of 20.8% and 19.1%, respectively, as compared to the standard NWP model. For long-range baseline (155 km and 207 km), the corrected NWP-constrained RTK showed a reduction of over 43% in the initialization time compared with the standard RTK, and showed a reduction of over 24% in the initialization time compared with the standard NWP-constrained RTK.

Column Integrated Water Vapor and Aerosol Load Characterization with the New ZEN-R52 Radiometer

The study shows the first results of the column-integrated water vapor retrieved by the new ZEN-R52 radiometer. This new radiometer has been specifically designed to monitor aerosols and atmospheric water vapor with a high degree of autonomy and robustness in order to allow the expansion of the observations of these parameters to remote desert areas from ground-based platforms. The ZEN-R52 device shows substantial improvements compared to the previous ZEN-R41 prototype: a smaller field of view, an increased signal-to-noise ratio, better stray light rejection, and an additional channel (940 nm) for precipitable water vapor (PWV) retrieval. PWV is inferred from the ZEN-R52 Zenith Sky Radiance (ZSR) measurements using a lookup table (LUT) methodology. The improvement of the new ZEN-R52 in terms of ZSR was verified by means of a comparison with the ZEN-R41, and with the Aerosol Robotic Network (AERONET) Cimel CE318 (CE318-AERONET) at Izaña Observatory, a Global Atmosphere Watch (GAW) high mountain station (Tenerife, Canary Islands, Spain), over a 10-month period (August 2017 to June 2018). ZEN-R52 aerosol optical depth (AOD) was extracted by means of the ZEN–AOD–LUT method with an uncertainty of ±0.01 ± 0.13*AOD. ZEN-R52 PWV extracted using a new LUT technique was compared with quasi-simultaneous (±30 s) Fourier Transform Infrared (FTIR) spectrometer measurements as reference. A good agreement was found between the two instruments (PWV means a relative difference of 9.1% and an uncertainty of ±0.089 cm or ±0.036 + 0.061*PWV for PWV <1 cm). This comparison analysis was extended using two PWV datasets from the same CE318 reference instrument at Izaña Observatory: one obtained from AERONET (CE318-AERONET), and another one using a specific calibration of the 940-nm channel performed in this work at Izaña Atmospheric Research Center Observatory (CE318-IARC), which improves the PWV product.

Hourly PWV Dataset Derived from GNSS Observations in China

The rapid variation of atmospheric water vapor is important for a regional hydrologic cycle and climate change. However, it is rarely investigated in China, due to the lack of a precipitable water vapor (PWV) dataset with high temporal resolution. Therefore, this study focuses on the generation of an hourly PWV dataset using Global Navigation Satellite System (GNSS) observations derived from the Crustal Movement Observation Network of China. The zenith total delay parameters estimated by GAMIT/GLOBK software are used and validated with an average root mean square (RMS) error of 4-5 mm. The pressure (P) and temperature (T) parameters used to calculate the zenith hydrostatic delay (ZHD) and weighted average temperature of atmospheric water vapor (T_m) are derived from the fifth-generation reanalysis dataset of the European Centre for Medium-Range Weather Forecasting (ECMWF ERA5) products. The values of P and T at the GNSS stations are obtained by interpolation in the horizontal and vertical directions using empirical formulas. T_m is calculated at the GNSS stations using the improved global pressure and temperature 2 wet (IGPT2w) model in China with an RMS of 3.32 K. The interpolated P and T are validated by interpolating the grid-based ERA5 data into radiosonde stations. The average RMS and bias of P and T in China are 2.71/-1.11 hPa and 1.88/-0.51 K, respectively. Therefore, the error in PWV with a theoretical RMS of 1.85 mm over the period of 2011-2017 in China can be obtained. Finally, the hourly PWV dataset in China is generated and the practical accuracy of the generated PWV dataset is validated using the corresponding AERONET and radiosonde data at specific stations. Numerical results reveal that the average RMS values of the PWV dataset in the four geographical regions of China are less than 3 mm. A case analysis of the PWV diurnal variations as a response to the EI Niño event of 2015-2016 is performed. Results indicate the capability of the hourly PWV dataset of monitoring the rapid water vapor changes in China.

Analysis and Comparison of GPS Precipitable Water Estimates between Two Nearby Stations on Tahiti Island

Since Bevis first proposed Global Positioning System (GPS) meteorology in 1992, the precipitable water (PW) estimates retrieved from Global Navigation Satellite System (GNSS) networks with high accuracy have been widely used in many meteorological applications. The proper estimation of GNSS PW can be affected by the GNSS processing strategy as well as the local geographical properties of GNSS sites. To better understand the impact of these factors, we compare PW estimates from two nearby permanent GPS stations (THTI and FAA1) in the tropical Tahiti Island, a basalt shield volcano located in the South Pacific, with a mean slope of 8% and a diameter of 30 km. The altitude difference between the two stations is 86.14 m, and their horizontal distance difference is 2.56 km. In this paper, Bernese GNSS Software Version 5.2 with precise point positioning (PPP) and Vienna mapping function 1 (VMF1) was applied to estimate the zenith tropospheric delay (ZTD), which was compared with the International GNSS Service (IGS) Final products. The meteorological parameters sourced from the European Center for Medium-Range Weather Forecasts (ECMWF) and the local weighted mean temperature ( T m ) model were used to estimate the GPS PW for three years (May 2016 to April 2019). The results show that the differences of PW between two nearby GPS stations is nearly a constant with value 1.73 mm. In our case, this difference is mainly driven by insolation differences, the difference in altitude and the wind being only second factors.

GPS-PWV based Improved Long-Term Rainfall Prediction Algorithm for Tropical Regions

Global positioning system (GPS) satellite delay is extensively used in deriving the precipitable water vapor (PWV) with high spatio–temporal resolution. One of the recent applications of GPS derived PWV values are to predict rainfall events. In the literature, there are rainfall prediction algorithms based on GPS-PWV values. Most of the algorithms are developed using data from temperate and sub-tropical regions. Mostly these algorithms use maximum PWV rate, maximum PWV variation and monthly PWV values as a criterion to predict the rain events. This paper examines these algorithms using data from the tropical stations and proposes the use of maximum PWV value for better prediction. When maximum PWV value and maximum rate of increment criteria are implemented on the data from the tropical stations, the false alarm ( F A ) rate is reduced by almost 17% as compared to the results from the literature. There is a significant reduction in F A rates while maintaining the true detection ( T D ) rates as high as that of the literature. A study done on the varying historical length of data and lead time values shows that almost 80% of the rainfall can be predicted with a false alarm of 26 . 4 % for a historical data length of 2 hours and a lead time of 45 min to 1 hour.

Comparison Analysis of Total Precipitable Water of Satellite-Borne Microwave Radiometer Retrievals and Island Radiosondes

Total precipitable water (TPW) of satellite-borne microwave radiometer retrievals is compared with the data that were collected from 49 island radiosonde stations for the period 2007–2015. Great consistency was found between TPW measurements made by radiosonde and eight satellite-borne microwave radiometers, including SSMI-F13, SSMI-F14, SSMIS-F16, SSMIS-F17, AMSR-E, AMSR-2, GMI, and WindSat. Mean values of the TPW differences for eight satellites ranged from −0.51 to 0.38mm, both root mean square errors and standard deviations were around 3mm, and all of the correlation coefficients between satellite TPW retrievals and radiosonde TPW for each satellite can reach 0.99. Subsequently, an analysis of the comparison results was conducted, which revealed three problems in the satellite TPW retrieval and two problems in radiosonde data. For TPW retrievals of satellite, when the values are above 60 mm, the precision of TPW retrieval significantly decreases with a distinct dry bias, which can reach 4 mm; additionally, abias related to wind speed and the uncertainty with the TPW retrieval in the presence of rain, which is stronger than 1mm/h, was found. The TPW measurements of radiosonde made by the type of IM-MK3 from India were quite unreliable, and almost all of the radiosonde data during the daytime were plagued by a dry bias.

Biases and Abrupt Shifts of Monthly Precipitable Water from Terra MODIS

Monthly atmospheric precipitable water (PW) from the Moderate Resolution Imaging Spectroradiometer (MODIS) aboard the Terra satellite was assessed over land at 60°S–60°N. MODIS provides two PW products by using infrared (IR) and near-IR (NIR) algorithms, respectively. An assessment was performed for both MODIS PW data from 2000 to 2014, comparing them with the measurements at international stations of the global positioning systems and with a reanalysis to detect abrupt changes through monthly variations. It is noted that MODIS IR systematically underestimated PW in over 75% of stations, and that PW estimation declines with time. MODIS NIR significantly overestimated PW for tropical land and experienced two abrupt shifts. These data defects result in large spurious decreasing trends in MODIS IR and increasing trends in MODIS NIR. The two MODIS PW products are currently not suitable for a climatic-trend analysis, highlighting the need for data reprocessing and calibration.

Real-time tropospheric delay retrieval with GPS, GLONASS, Galileo and BDS data

The precise point positioning (PPP) is a promising technology for the real-time retrieval of atmospheric parameters with a single receiver in anywhere, all-weather and any time. The real-time atmospheric parameters can be applied to the time-critical meteorology, such as the severe weather nowcasting. The PPP is a satellite-based technology. Multi-constellation integration can enhance satellite geometry and increase measurement redundancy so that the solutions of atmospheric parameters are expected to be improved. Currently, the Global Navigation Satellite System (GNSS) family includes recovered GLONASS and modernized GPS as well as the emerging Galileo and BDS. A week of GNSS observations from 160 stations are processed to retrieve the tropospheric zenith total delay (ZTD) in real time. The four-constellation mixed real-time precise products including satellite orbit and clock corrections are adopted, and their quality is evaluated. The performance of ZTD estimates is assessed in terms of accuracy and convergence time by comparing with final tropospheric ZTD products provided by two analysis centers. The ZTDs retrieved from different constellation combinations (i.e., GPS/GLONASS/Galileo/BDS, GPS/GLONASS, and GPS-only), different processing models for ionospheric delays (i.e., ionospheric-free (IF) combined PPP, and uncombined (UC) PPP), and different modes (i.e., real-time mode, and post-processing mode) are compared.

A New Zenith Tropospheric Delay Grid Product for Real-Time PPP Applications over China

Tropospheric delay is one of the major factors affecting the accuracy of electromagnetic distance measurements. To provide wide-area real-time high precision zenith tropospheric delay (ZTD), the temporal and spatial variations of ZTD with altitude were analyzed on the bases of the latest meteorological reanalysis product (ERA-Interim) provided by the European Center for Medium-Range Weather Forecasts (ECMWF). An inverse scale height model at given locations taking latitude, longitude and day of year as inputs was then developed and used to convert real-time ZTD at GPS stations in Crustal Movement Observation Network of China (CMONOC) from station height to mean sea level (MSL). The real-time ZTD grid product (RtZTD) over China was then generated with a time interval of 5 min. Compared with ZTD estimated in post-processing mode, the bias and error RMS of ZTD at test GPS stations derived from RtZTD are 0.39 and 1.56 cm, which is significantly more accurate than commonly used empirical models. In addition, simulated real-time kinematic Precise Point Positioning (PPP) tests show that using RtZTD could accelerate the BDS-PPP convergence time by up to 32% and 65% in the horizontal and vertical components (set coordinate error thresholds to 0.4 m), respectively. For GPS-PPP, the convergence time using RtZTD can be accelerated by up to 29% in the vertical component (0.2 m).