A Comprehensive Evaluation of Key Tropospheric Parameters from ERA5 and MERRA-2 Reanalysis Products Using Radiosonde Data and GNSS Measurements

Tropospheric delay is a major error source in the Global Navigation Satellite System (GNSS), and the weighted mean temperature (T-m) 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, T-m, 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), T-m, 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 T-m 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.