Tropospheric delay is one of the key factors that may affect high-precision satellite navigation and positioning and water vapor retrieval performance. Its variation in the vertical direction is much greater than that in the horizontal direction. Therefore, the vertical profile model of zenith total delay (ZTD) is important for the spatial interpolation of high-precision ZTD products and the development of ZTD models. However, in China, low spatial and temporal resolutions remain persistent in ZTD vertical profile models and limit their application. In this study, ZTD vertical profile grid models (CZTD-H model: CZTD-HM and CZTD-HE models) were developed by considering the time-varying height scale factor for China and employing ZTD layered profile information with high temporal-spatial resolution calculated using MERRA-2 data and ERA5 data based on the integration method during 2012–2016. The CZTD-H model accuracy was verified using the global navigation satellite system (GNSS) data acquired from the Crustal Movement Observation Network of China (CMONOC) and radiosonde data as reference and was compared with the canonical GPT3 model accuracy. The applicability of CZTD-HM and CZTD-HE models was discussed. The results showed that: (1) CZTD-HM and CZTD-HE models exhibited excellent performance for ZTD layered vertical interpolation in northwestern and southeastern China, respectively. Among ZTD layered profiles from 84 radiosonde stations, the RMSE of ZTD data interpolated using CZTD-HM and CZTD-HE models improved by 12.9/16.23% and 13.8/17.16% compared with GPT3-1 and GPT3-5 models, respectively. (2) The CZTD-H model maintained high performance for the spatial interpolation of GGOS grid ZTD data. Validation with ZTD data from 249 GNSS stations showed that the RMSEs of both CZTD-HM and CZTD-HE models improved by 2.8 mm (19.7%) and 2.6 mm (18.6%) compared with those of the GPT3-1 and GPT3-5 models, respectively. The CZTD-HE model showed excellent performance in summer among all the models. Only the location and day of the year were required for the application of the CZTD-H model, which showed excellent ZTD vertical correction performance in China. With the different performances of the CZTD-HE and CZTD-HM models in China, the ERA5 model can be recommended for practical applications. Therefore, these results can provide a reference for the data source selection of ZTD vertical profile model construction on the basis of high-precision reanalysis data, GNSS real-time high-precision positioning, and GNSS meteorology in China.