A numerical method to improve the spatial interpolation of water vapor from numerical weather models: a case study in South and Central America

Abstract

Commonly, numerical weather model (NWM) users can get the vertically integrated water vapor (IWV) value at a given location from the values at nearby grid points. In this study we used a validated and freely available global navigation satellite system (GNSS) IWV data set to analyze the very well-known effect of height differences. To this end, we studied the behavior of 67 GNSS stations in Central and South America with the prerequisite that they have a minimum of 5 years of data during the period from 2007 to 2013. The values of IWV from GNSS were compared with the respective values from ERA-Interim and MERRA-2 from the same period. Firstly, the total set of stations was compared in order to detect cases in which the geopotential difference between GNSS and NWM required correction. An additive integral correction to the IWV values from ERA-Interim was then proposed. For the calculation of this correction, the multilevel values of specific humidity and temperature given at 37 pressure levels by ERA-Interim were used. The performance of the numerical integration method was tested by accurately reproducing the IWV values at every individual grid point surrounding each of the GNSS sites under study. Finally, considering the IWVGNSS values as a reference, the improvement introduced to the IWVERA-Interim values after correction was analyzed. In general, the corrections were always recommended, but they are not advisable in marine coastal areas or on islands as at least two grid points of the model are usually in the water. In such cases, the additive correction could overestimate the IWV.