Precise modeling of tropospheric delay and weighted mean temperature (T m ) is critical for Global Navigation Satellite System (GNSS) positioning and meteorology. However, the model data in previous models cover a limited time span, which limits the accuracy of these models. Besides, the vertical variations of tropospheric delay and T m are not perfectly modeled in previous studies, which affects the performance of height corrections. In this study, we used the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalysis from 1979 to 2017 to build a new empirical model. We first carefully modeled the lapse rates of tropospheric delay and T m . Then we considered the temporal variations by linear trends, annual, and semi-annual variations and the spatial variations by grids. This new model can provide zenith hydrostatic delay (ZHD), zenith wet delay (ZWD), and T m worldwide with a spatial resolution of 1 • × 1 • . We used the ECMWF ERA-Interim data and the radiosonde data in 2018 to validate this new model in comparison with the canonical GPT2w model. The results show that the new model has higher accuracies than the GPT2w model in all parameters. Particularly, this new model largely improves the accuracy in estimating ZHD and T m at high-altitude (relative to the grid point height) regions.