Evaluation and calibration of MODIS near-infrared precipitable water vapor over China using GNSS observations and ERA-5 reanalysis dataset

Water vapor is one of the most important parameters in climatic studies. MODerate-resolution Imaging Spectroradiometer (MODIS) is a key instrument and can provide spatially continuous precipitable water vapor (PWV) products. This study was focused on the performance evaluation of the MODIS near-infrared PWV product (MOD-NIR-PWV) over China. For a comprehensive assessment of the performance of MOD-NIR-PWV, PWV retrieved from the measurements at the global navigation satellite systems (GNSS) stations (i.e., GNSS-PWV) and the ERA5 reanalysis dataset (ERA-PWV) from 2013 to 2018 were used as the reference. To investigate the suitability of using ERA-PWV as the reference for the evaluation, ERA-PWV was compared to the high-accuracy GNSS-PWV at 246 GNSS stations and PWV retrieved from radiosonde observations (RS-PWV) at 78 radiosonde stations over China. The results showed that the mean bias and mean root-mean-square (RMS) of the differences between ERA-PWV and GNSS-PWV across all the stations were 0.5 and 1.7 mm, respectively, and the mean correlation coefficient of the two datasets was above 0.96. The values were 0.4 and 1.9 mm and 0.97, respectively, for the differences between ERA-PWV and RS-PWV. This suggests the suitability of ERA-PWV as the reference for the evaluation of MOD-NIR-PWV. In addition, MOD-NIR-PWV was compared with both GNSS-PWV and ERA-PWV, and their mean bias and mean RMS were 2.9 and 3.8 mm (compared to GNSS-PWV) and 2.1 and 3.0 mm (compared to ERA-PWV), respectively. The positive bias values and the non-normal distribution of the differences between MOD-NIR-PWV and both reference datasets imply that a considerable systematic overestimation of MOD-NIR-PWV over China may exist. To mitigate the systematic bias, ERA-PWV was utilized as the sample data due to its spatial continuities, and a grid-based calibration model was developed based on the annual and semiannual periodicities in the differences between MOD-NIR-PWV and ERA-PWV at each grid point. Afte