A Zenith Tropospheric Delay Modeling Method Based on the UNB3m Model and Kriging Spatial Interpolation

To accurately estimate Zenith Tropospheric Delay (ZTD) for high-precision positioning of the Global Navigation Satellite System (GNSS), this study proposes a modeling method of ZTD based on the UNB3m model and Kriging spatial interpolation, in which the optimal spatial interpolation parameters are determined based on the errors corresponding to different combinations of the interpolation parameters, and the spatial distribution of the GNSS modeling stations is determined by the interpolation errors of the randomly selected GNSS stations for several times. To verify the accuracy and reliability of the proposed model, the ZTD estimates of 132,685 epochs with 1 h or 2 h temporal resolution for 28 years from 1997 to 2025 of the global network of continuously operating GNSS tracking stations are used as inputs; the ZTD results at any position and the corresponding observation moment can be obtained with the proposed model. The experimental results show that the model error is less than 30 mm in more than 85% of the observation epochs, the ZTD estimation results are less affected by the horizontal position and height of the GNSS stations than traditional models, and the ZTD interpolation error is improved by 10–40 mm compared to the GPT3 and UNB3m models at the four GNSS checking stations. Therefore, this technology can provide ZTD estimation results for single- and dual-frequency hybrid deformation monitoring, as well as dense ZTD data for Precipitable Water Vapor (PWV) inversion. Since the proposed method has the advantages of simple implementation, high accuracy, high reliability, and ease of promotion, it is expected to be fully applied in other high-precision positioning applications.