A New Approach for the Development of Grid Models Calculating Tropospheric Key Parameters over China

: Pressure, water vapor pressure, temperature, and weighted mean temperature (T m ) are tropospheric parameters that play an important role in high-precision global navigation satellite system navigation (GNSS). As accurate tropospheric parameters are obligatory in GNSS navigation and GNSS water vapor detection, high-precision modeling of tropospheric parameters has gained widespread attention in recent years. A new approach is introduced to develop an empirical tropospheric delay model named the China Tropospheric (CTrop) model, providing meteorological parameters based on the sliding window algorithm. The radiosonde data in 2017 are treated as reference values to validate the performance of the CTrop model, which is compared to the canonical Global Pressure and Temperature 3 (GPT3) model. The accuracy of the CTrop model in regards to pressure, water vapor pressure, temperature, and weighted mean temperature are 5.51 hPa, 2.60 hPa, 3.09 K, and 3.35 K, respectively, achieving an improvement of 6%, 9%, 10%, and 13%, respectively, when compared to the GPT3 model. Moreover, three different resolutions of the CTrop model based on the sliding window algorithm are also developed to reduce the amount of gridded data provided to the users, as well as to speed up the troposphere delay computation process, for which users can access model parameters of different resolutions for their requirements. With better accuracy of estimating the tropospheric parameters than that of the GPT3 model, the CTrop model is recommended to improve the performance of GNSS positioning and navigation.


Introduction
Tropospheric delay is one of the major factors affecting global navigation satellite system (GNSS) positioning. The vertical distribution of tropospheric parameters is strongly affected by Earth's gravity. It is difficult to accurately monitor and invert the description and modeling of tropospheric parameters under the influence of many factors. In the geodetic analysis of GNSS and Very Long Baseline Interferometry (VLBI) observations, the slant tropospheric delays are generally mapped into the zenith direction (zenith tropospheric delay, ZTD) through mapping functions. The ZTD is separated into two parts, namely, the zenith hydrostatic delay (ZHD) and the zenith wet delay (ZWD). The ZHD is usually calculated by the Saastamoinen model [1] from surface pressure measurements. The ZWD is estimated in the data analysis or approximated by models such as the Saastamoinen model, which uses parameters such as the temperature and water vapor pressure as input. Furthermore, the weighted mean temperature (T m ) is also an essential parameter that is indispensable in obtaining precipitable water vapor (PWV) from the zenith wet delay. Therefore, pressure, water vapor pressure, temperature, and T m are key parameters that are used as input data to obtain tropospheric delay information. It is of great significance to Remote Sens. 2021, 13, 3546 2 of 18 monitor temperature, pressure, water vapor pressure, T m , and other tropospheric parameters for high-precision GNSS positioning and navigation, as well as GNSS meteorologic applications [2,3].
In addition, floods and droughts have been caused by extreme weather frequently. Severe weather, such as strong regional convection, and large-scale climate anomalies have brought on heavy economic losses and casualties in recent years [4][5][6]. With increasing requirements for high-precision GNSS positioning, and in-depth research in fields such as climate change and extreme weather generation mechanisms, the demand for real-time and high-spatiotemporal-resolution tropospheric parameters is increasing.
There are two main types of tropospheric parameter models according to the methods used for modeling. The Hopfield model [7], Saastamoinen model, and Black model [8] are the classic tropospheric delay models which are based on in situ meteorological parameters. Since conventional GNSS receivers are unavailable for measurements, the use of the aforementioned classic tropospheric delay models is restricted to a certain extent. As space technologies such as GNSS and VLBI are widely used in the navigation and guidance of various space vehicles as well as for the needs of climate change and weather forecasting, the real-time and high spatial resolution of meteorological parameters are particularly important [9]. Therefore, developing regional or global empirical tropospheric delay models based on atmosphere analysis data has attracted widespread attention [10][11][12][13][14][15][16]. Empirical meteorological models, such as the University of New Brunswick (UNB) models [17,18], the European Geo-stationary Navigation Overlay System (EGNOS) model [19], the TropGrid model [20,21], and the GPT models [22][23][24][25] have been developed, aimed at directly obtaining high-precision tropospheric parameters through the model without measured meteorological parameters. The UNB3 model was stored by tabulated meteorological data, which divides the latitude into five intervals. The EGNOS model was established by simplifying the UNB3 model, and the formula is different from that of the UNB3 model. Only the day of the year (DOY) and the location of the station are needed when applying this model, which is used in GNSS satellite navigation enhancement systems. The authors of [20] developed the TropGrid model considering the diurnal variations of parameters, but neglect the semiannual variations. The GPT model was established based on the ERA-40 reanalysis data. However, the GPT model only considers the annual variations of the parameters. The authors of [23] developed the GPT2 model, which considers the semiannual variations of parameters, by analyzing 10-year ECMWF ERA-Interim data with resolutions of 5 • × 5 • and 1 • × 1 • . Considered to be the tropospheric model with the highest precision for quite a long time, two meteorological parameters were added to the GPT2w model based on GPT2. GPT3 is the latest generation model for an upgraded version of GPT2w. Two parameters of gradients in the east and north direction were added to the GPT3 model [26,27].
The aforementioned tropospheric delay models have performed well in GNSS meteorology. In this work, a new approach is introduced to develop an empirical tropospheric delay model providing meteorological parameters based on the second modern-era retrospective analysis for research and applications (MERRA-2) data. The new empirical model could be used to serve real-time GNSS positioning and navigation.

Radiosonde Data
The radiosonde data provide measured meteorological parameters from the ground to a height of about 30 km at more than 1500 stations around the world. There are 89 radiosonde stations over China as shown in Figure 1. The radiosonde data are obtained based on the actual measurements of meteorological sensors on the sounding balloon, which has high accuracy and credibility. However, radiosonde balloons are extremely susceptible to other disturbance factors, such as weather factors and clouds, during the ascent. Different types of sounding sensors show their own system deviations, as well as equipment failures, which lead to phenomena such as missing data. Although radiosonde data shows some disadvantages, they are still widely used to validate measurement results [27,28]. susceptible to other disturbance factors, such as weather factors and clouds, during the ascent. Different types of sounding sensors show their own system deviations, as well as equipment failures, which lead to phenomena such as missing data. Although radiosonde data shows some disadvantages, they are still widely used to validate measurement results [27,28].

MERRA-2 Reanalysis Product Data
MERRA-2 is the latest atmospheric reanalysis product comprising of data beginning in 1980. It is provided by the National Aeronautics and Space Administration (NASA) [29,30] with a spatial resolution of 0.625° × 0.5° (lon. × lat.) and a temporal resolution of 6 h. MERRA-2 replaces the original MERRA reanalysis dataset using an upgraded version of the Goddard Earth Observing System Model, Version 5 (GEOS-5) data assimilation system [31], and the Gridpoint Statistical Interpolation (GSI) analysis scheme [32][33][34]. It is also the first long-term global reanalysis product that assimilates space-based observations of aerosols in the climate system and represents their interactions with physical processes [34]. In addition, MERRA-2 tends to minimize the spurious variations related to inhomogeneity in the observational records and achieves a global balance between evaporation and precipitation through the mass conservation constraint [35]. Temperature, pressure, and specific humidity can be obtained by MERRA-2 data. The water vapor pressure and Tm can be calculated by the following equations: where ℎ denotes the specific humidity, is the pressure, ℎ is the height, is the water vapor pressure, and is the temperature.

MERRA-2 Reanalysis Product Data
MERRA-2 is the latest atmospheric reanalysis product comprising of data beginning in 1980. It is provided by the National Aeronautics and Space Administration (NASA) [29,30] with a spatial resolution of 0.625 • × 0.5 • (lon. × lat.) and a temporal resolution of 6 h. MERRA-2 replaces the original MERRA reanalysis dataset using an upgraded version of the Goddard Earth Observing System Model, Version 5 (GEOS-5) data assimilation system [31], and the Gridpoint Statistical Interpolation (GSI) analysis scheme [32][33][34]. It is also the first long-term global reanalysis product that assimilates space-based observations of aerosols in the climate system and represents their interactions with physical processes [34]. In addition, MERRA-2 tends to minimize the spurious variations related to inhomogeneity in the observational records and achieves a global balance between evaporation and precipitation through the mass conservation constraint [35]. Temperature, pressure, and specific humidity can be obtained by MERRA-2 data. The water vapor pressure and T m can be calculated by the following equations: where sh denotes the specific humidity, P is the pressure, h is the height, e is the water vapor pressure, and T is the temperature.

. Analysis of Tropospheric Parameters
Analysis of the spatiotemporal characteristics of tropospheric parameters plays an important part in modeling. Some reports have suggested that T m or temperature are related to height [36,37]. Four representative grid points over China shown in Figure 1 are chosen to determine the daily mean temperature in 2016 for analyzing the correlations between temperature and height. The temperature in each standard pressure level of MERRA-2 is interpolated upon a number of the same heights, and the result is shown in Figure 2. When geopotential height increases, the temperature decreases. There is a correlation between the temperature and the geopotential height of the four grid points over China, which can be expressed as follows: where γ denotes the temperature lapse rate, δh is the height, and k is a constant. Analysis of the spatiotemporal characteristics of tropospheric parameters plays an important part in modeling. Some reports have suggested that Tm or temperature are related to height [36,37]. Four representative grid points over China shown in Figure 1 are chosen to determine the daily mean temperature in 2016 for analyzing the correlations between temperature and height. The temperature in each standard pressure level of MERRA-2 is interpolated upon a number of the same heights, and the result is shown in Figure 2. When geopotential height increases, the temperature decreases. There is a correlation between the temperature and the geopotential height of the four grid points over China, which can be expressed as follows: where denotes the temperature lapse rate, ℎ is the height, and is a constant. The Bevis formula [38,39] shows the relationship between the surface temperature and Tm, which can be expressed as Tm = 70.2 + 0.72Ts; thus, the spatial and temporal characteristics of Tm are similar to those of Ts. However, the Bevis formula is only an approximate relationship, and separate models for Tm and temperature are developed in this paper. Tm, temperature, pressure, and water vapor pressure data provided by the MERRA-2 reanalysis data from 2012 to 2016 are divided into three intervals (15° N-30° N, 30° N-40° N, and 40° N-55° N) according to latitude in order to calculate the time series and to analyze the correlations between the parameters and latitude; the result is shown in Figure  3. All the tropospheric parameters show obvious characteristics of annual and semiannual variations. The peak value of the Tm, temperature, and water vapor pressure in one year appears in the middle of the year, showing symmetrical distribution. In addition, as the The Bevis formula [38,39] shows the relationship between the surface temperature and T m , which can be expressed as T m = 70.2 + 0.72T s ; thus, the spatial and temporal characteristics of T m are similar to those of T s . However, the Bevis formula is only an approximate relationship, and separate models for T m and temperature are developed in this paper. the aforementioned values gradually decrease, which indicates the correlation to latitude. The T m ranges from 245 K to 280 K in the high latitude and middle latitude regions, while it ranges from 280 K to 285 K in the low latitude regions. The temperature ranges from 250 K to 295 K in the high latitude regions and middle latitude areas, while it ranges from 290 K to 300 K in the low latitude regions, which displays more stable temperatures. Pressure shows the largest value in the low latitude regions, while it appears smallest in the middle latitude regions. The water vapor pressure value ranges from 0 hPa to 20 hPa in the high latitude regions and middle latitude areas, while it ranges from 10 hPa to 30 hPa in the low latitude areas, which appears more stable.
Remote Sens. 2021, 13, x FOR PEER REVIEW 5 of 17 latitude increases, the aforementioned values gradually decrease, which indicates the correlation to latitude. The Tm ranges from 245 K to 280 K in the high latitude and middle latitude regions, while it ranges from 280 K to 285 K in the low latitude regions. The temperature ranges from 250 K to 295 K in the high latitude regions and middle latitude areas, while it ranges from 290 K to 300 K in the low latitude regions, which displays more stable temperatures. Pressure shows the largest value in the low latitude regions, while it appears smallest in the middle latitude regions. The water vapor pressure value ranges from 0 hPa to 20 hPa in the high latitude regions and middle latitude areas, while it ranges from 10 hPa to 30 hPa in the low latitude areas, which appears more stable. MERRA-2 reanalysis data are used to analyze the spatial and temporal characteristics of meteorological parameters. It can be found that these tropospheric parameters are related to latitude and height, and they all have the characteristics of annual and semiannual changes. To further analyze the influence of spatial factors on the tropospheric parameters, the annual mean of the tropospheric parameters at each MERRA-2 grid point over China are determined, and the result is shown in Figure 4. The annual mean values of the tropospheric parameters in the western region are lower than those in other regions over China due to the high altitude in this area. It also has obvious characteristics indicating that the annual average values of the tropospheric parameters in high latitude areas are lower than those in the low-and middle-latitude areas. Furthermore, all the parameters show significant correlations for latitude and longitude. MERRA-2 reanalysis data are used to analyze the spatial and temporal characteristics of meteorological parameters. It can be found that these tropospheric parameters are related to latitude and height, and they all have the characteristics of annual and semiannual changes. To further analyze the influence of spatial factors on the tropospheric parameters, the annual mean of the tropospheric parameters at each MERRA-2 grid point over China are determined, and the result is shown in Figure 4. The annual mean values of the tropospheric parameters in the western region are lower than those in other regions over China due to the high altitude in this area. It also has obvious characteristics indicating that the annual average values of the tropospheric parameters in high latitude areas are lower than those in the low-and middle-latitude areas. Furthermore, all the parameters show significant correlations for latitude and longitude. Remote Sens. 2021, 13, x FOR PEER REVIEW 6 of 17

Analysis of the Characteristics of the Lapse Rate
Analysis of the spatiotemporal characteristics of the lapse rate parameters also plays an important part in modeling. When ruling out the differences in elevation data between the different data sources, such as the ellipsoidal height and the geopotential height, the elevation of GNSS and radiosonde stations is inconsistent with the height of the grid point. It must be considered that height correction in the vertical dimension plays an important role in modeling, because these tropospheric parameters are sensitive to height, showing notable changes in the vertical direction. If their vertical changes are properly considered, the model will be able to have better performance at different heights.
The changes in temperature and Tm along the vertical direction are usually adopted to be expressed by the linear function [20,40]. The exponential function is usually adopted to express the changes in pressure [24,25] and water vapor pressure along the vertical direction [21]. In this paper, conventional methods are used to carry out height correction, and they have been adopted in the GPT3 model and other tropospheric delay models. Equations (4)-(7) are for Tm, temperature, pressure, water vapor pressure, and pressure:

Analysis of the Characteristics of the Lapse Rate
Analysis of the spatiotemporal characteristics of the lapse rate parameters also plays an important part in modeling. When ruling out the differences in elevation data between the different data sources, such as the ellipsoidal height and the geopotential height, the elevation of GNSS and radiosonde stations is inconsistent with the height of the grid point. It must be considered that height correction in the vertical dimension plays an important role in modeling, because these tropospheric parameters are sensitive to height, showing notable changes in the vertical direction. If their vertical changes are properly considered, the model will be able to have better performance at different heights.
The changes in temperature and T m along the vertical direction are usually adopted to be expressed by the linear function [20,40]. The exponential function is usually adopted to express the changes in pressure [24,25] and water vapor pressure along the vertical direction [21]. In this paper, conventional methods are used to carry out height correction, and they have been adopted in the GPT3 model and other tropospheric delay models. Equations (4)-(7) are for T m , temperature, pressure, water vapor pressure, and pressure: where T U m , T U , e U , and P U refer to the meteorological values at the station; T G m , T G , e G , and P G represent the meteorological values at the grid point; H U and H G denote the height at the station and the grid point, respectively; and γ, β, and λ are the lapse rates. M refers to the molar mass of dry air (28.965·10 −3 kg/mol), R represents the universal gas constant (8.31432 J/K·mol), and g denotes the gravitational coefficient. A new parameter τ is introduced to show the proportion of temperature and temperature lapse rate, which can be expressed as: g·M R·β (8) λ can be obtained by Equation (6) or (9). In this work, λ is obtained by fitting Equation (9): where k 2 and k 3 denote the refractive index constants. R d refers to the specific gas constant for the dry constituents [41].
To further analyze the distribution of the annual mean of the lapse rate parameters, the data over China are determined from 2012 to 2016, and the result is shown in Figure 5. The annual mean value of the pressure lapse rate, temperature lapse rate, and T m lapse rate shows obvious geographical characteristics over China. The values are smaller in the west of China than in other regions, which is due to the higher elevation in the western region. The distribution of the water vapor decrease factor is different from the others, which shows that the values in the northeast and southwest areas are larger than that in the other regions. Furthermore, all the parameters show significant correlations for latitude and longitude.  values at the station;  ,  ,  ,  and  represent the meteorological values at the grid point;  and  denote the  height at the station and the grid point, respectively; and , , and are the lapse rates. refers to the molar mass of dry air (28.965·10 −3 kg/mol), R represents the universal gas constant (8.31432 J/K·mol), and denotes the gravitational coefficient. A new parameter is introduced to show the proportion of temperature and temperature lapse rate, which can be expressed as: λ can be obtained by Equation (6) or (9). In this work, λ is obtained by fitting Equation (9): where and denote the refractive index constants. refers to the specific gas constant for the dry constituents [41].
To further analyze the distribution of the annual mean of the lapse rate parameters, the data over China are determined from 2012 to 2016, and the result is shown in Figure  5. The annual mean value of the pressure lapse rate, temperature lapse rate, and Tm lapse rate shows obvious geographical characteristics over China. The values are smaller in the west of China than in other regions, which is due to the higher elevation in the western region. The distribution of the water vapor decrease factor is different from the others, which shows that the values in the northeast and southwest areas are larger than that in the other regions. Furthermore, all the parameters show significant correlations for latitude and longitude.

Development of the CTrop Model
Some shortcomings still exist in tropospheric models that have been developed, such as only single gridded data used for modeling. In this work, a sliding window algorithm is introduced to develop the tropospheric delay model, as shown in Figure 6. The new approach is instituted to divide the area of China into regular windows of the same size. Model parameters are estimated based on the data in each window to be taken as results of the center point of the sliding window. The realization process of the sliding window algorithm has been described in previous works [36,42].

Development of the CTrop Model
Some shortcomings still exist in tropospheric models that have been developed, such as only single gridded data used for modeling. In this work, a sliding window algorithm is introduced to develop the tropospheric delay model, as shown in Figure 6. The new approach is instituted to divide the area of China into regular windows of the same size. Model parameters are estimated based on the data in each window to be taken as results of the center point of the sliding window. The realization process of the sliding window algorithm has been described in previous works [36,42].
where is the meteorological parameters, such as temperature, Tm, pressure, and water vapor pressure; is the latitude; is the longitude; is the annual average value of the meteorological parameters; is the latitude correction; is the longitude correction; and are the annual amplitude coefficients of the meteorological parameters; and are the semiannual amplitude coefficients of the meteorological parameters, and DOY is the day of the year.
The elevation of the grid points in the atmospheric reanalysis data is inconsistent with the elevation of GNSS stations. The height correction plays an important role in As mentioned above, temperature, T m , pressure, and water vapor pressure exhibit evident seasonal, latitudinal, and longitudinal characteristics over China, which should be taken into account to obtain a high-precision model. The equation in each window is expressed as follows:  (10) where MP is the meteorological parameters, such as temperature, T m , pressure, and water vapor pressure; ϕ is the latitude; θ is the longitude; α 1 is the annual average value of the meteorological parameters; α 2 is the latitude correction; α 3 is the longitude correction; α 4 and α 5 are the annual amplitude coefficients of the meteorological parameters; α 6 and α 7 are the semiannual amplitude coefficients of the meteorological parameters, and DOY is the day of the year. The elevation of the grid points in the atmospheric reanalysis data is inconsistent with the elevation of GNSS stations. The height correction plays an important role in modeling, for which the lapse rates of meteorological parameters should be considered for height correction in the vertical dimension. The equation is written as follows:  (11) where LR is the lapse rates of parameters in the vertical dimension, such as γ, β, λ, and τ; δ 1 is the annual mean value of the lapse rate of parameters; δ 2 is the latitude correction; δ 3 is the longitude correction; δ 4 and δ 5 are the annual amplitude coefficients of the lapse rate parameters; δ 6 and δ 7 are the semiannual amplitude coefficients of the meteorological parameters. The coefficients are calculated by least-squares adjustment in each window over China from 2012 to 2016, and a grid model calculating tropospheric key parameters over China named the CTrop model is established with a spatial resolution of 1.25 For horizontal interpolation, the inverse distance weighted and bilinear methods are commonly used. Considering the fact that changes in latitude have an impact on tropospheric parameters when developing the model, the inverse distance weighted method can reduce the impact on the interpolation results in the latitude direction, and grid points farther from the user have less of an impact on the interpolation, so the inverse distance weighted method is used for horizontal interpolation.
Only the day of the year (DOY) and the location of the station are needed when applying this model, which makes it very convenient. First, the four grid points nearest to the location are identified. The parameters of these four points at the height are then calculated. Finally, inverse distance weighted interpolation is employed to interpolate the required parameters at the given location.

Analysis of the Accuracy of the CTrop Model
The performance of the CTrop model is validated by radiosonde data over China in 2017 compared with the GPT3 model at a resolution of 1 • × 1 • . The results are summarized in Table 1.  Table 1 lists the accuracy of the CTrop model for meteorological parameters in comparison with the GPT3 model. In terms of water vapor pressure, the RMS of the CTrop model is 2.60 hPa and is smaller than that of the GPT3 model, which decreases by 9%. As for pressure, both models reveal a negative bias. Although the CTrop model shows a larger bias than that of the GPT3 model, it attains a smaller RMS, which decreases by 6%. In terms of temperature, it also reveals a negative bias, which indicates that the value calculated by the CTrop and GPT3 model is smaller than that of the radiosonde data. The bias of the CTrop model is -0.11 K, which decreases by 91% compared to the GPT3 model. The performance of the CTrop model is superior to that of the GPT3 model, attaining an improvement of 10%. In terms of T m , the CTrop model exhibits a smaller error than that of the GPT3 model. The bias of the CTrop model is 0.19 K and the RMS of the CTrop model is 3.35 K, with reductions of 87% and 13% compared with that of the GPT3 model, respectively.
To analyze the spatial characteristics of the performance of the CTrop and GPT3 models over China, the accuracy of the tropospheric parameters at each radiosonde site are calculated, and the results are shown in Figures 7-10.
As can be observed in Figure 7, the overall distribution of performance of the pressure of the CTrop model is consistent with that of the GPT3 model over China, and the accuracy of both models is high in most radiosonde sites. The largest error of radiosonde sites for the CTrop model is the same as that of the GPT3 model. Furthermore, the performance of the CTrop model in the Taibei radiosonde site (25.03 • N, 121.51 • E) is better than that of the GPT3 model. Figure 8 shows that the GPT3 model attains the largest positive and negative bias in terms of water vapor pressure. The CTrop model shows a smaller RMS in the western area than in the eastern area. The performance of the CTrop and GPT3 models in the coastal areas is relatively low, the reason for which may be that the rainfall in the coastal areas is high, affected by the ocean climate. Figure 9 shows that both the GPT3 and CTrop models perform well in the southeast area for T m over China while displaying large errors in the northeast region. The CTrop model attains a smaller error than that of the GPT3 model in the west area over China. Figure 10 shows that most stations show positive bias in the southeast area and negative bias in the northwest area. As the latitude increases, the RMS accuracy shows a downward trend. The CTrop model performs better in high latitude than does the GPT3 model. In short, the CTrop model shows better accuracy in estimating the tropospheric parameters than that of the GPT3 model. To analyze the spatial characteristics of the performance of the CTrop and GPT models over China, the accuracy of the tropospheric parameters at each radiosonde sit are calculated, and the results are shown in Figures 7-10.
As can be observed in Figure 7, the overall distribution of performance of th pressure of the CTrop model is consistent with that of the GPT3 model over China, an the accuracy of both models is high in most radiosonde sites. The largest error o radiosonde sites for the CTrop model is the same as that of the GPT3 model. Furthermore the performance of the CTrop model in the Taibei radiosonde site (25.03°N, 121.51°E) i better than that of the GPT3 model. Figure 8 shows that the GPT3 model attains the larges positive and negative bias in terms of water vapor pressure. The CTrop model shows smaller RMS in the western area than in the eastern area. The performance of the CTro and GPT3 models in the coastal areas is relatively low, the reason for which may be tha the rainfall in the coastal areas is high, affected by the ocean climate. Figure 9 shows tha both the GPT3 and CTrop models perform well in the southeast area for Tm over Chin while displaying large errors in the northeast region. The CTrop model attains a smalle error than that of the GPT3 model in the west area over China. Figure 10 shows that mos stations show positive bias in the southeast area and negative bias in the northwest area As the latitude increases, the RMS accuracy shows a downward trend. The CTrop mode performs better in high latitude than does the GPT3 model. In short, the CTrop mode shows better accuracy in estimating the tropospheric parameters than that of the GPT model.      The height of the radiosonde site is divided into three bands, and the performance at different heights of the CTrop and GPT3 models are calculated. The results are listed in Table 2. The precision of the two models is related to height. With the increase in height, the performance of pressure and water vapor pressure is improved, but the temperature decreases. The accuracy of the Tm of the CTrop model is improved with the increase in height, which is contrary to that of the GPT3 model.

Analysis of the Accuracy of Different Resolutions of the CTrop Model
From the performance of the CTrop model compared with the GPT3 model, it can be observed that the CTrop model based on the sliding window algorithm performs superiorly to the GPT3 model. Since the resolutions of the GPT3 model are 1° × 1° and 5° × 5°, the applicability of the model is not sufficiently abundant. Based on the sliding window algorithm, users can access model parameters of different resolutions for their requirements, which represents one of its advantages. The goal behind the development of a lower horizontal resolution version of the CTrop model is the reduction in the amount The height of the radiosonde site is divided into three bands, and the performance at different heights of the CTrop and GPT3 models are calculated. The results are listed in Table 2. The precision of the two models is related to height. With the increase in height, the performance of pressure and water vapor pressure is improved, but the temperature decreases. The accuracy of the T m of the CTrop model is improved with the increase in height, which is contrary to that of the GPT3 model.   has better accuracy in the west of China than that of the GPT3-5 model. Consequently, the lower horizontal resolution version of the CTrop model shows better accuracy in estimating the tropospheric parameters than that of the GPT3 model.

Conclusions
In this work, the distribution characteristics of meteorological parameters are analyzed, and it is observed that the meteorological parameters exhibit major annual and semiannual periodic variations that are also related to latitude and longitude. Considering the spatial distribution and time-varying characteristics of the meteorological parameters a refined regional empirical model (CTrop) based on the sliding window algorithm is developed for the estimation of tropospheric key parameters over China. Only the day of the year (DOY) and the location of the station are needed when applying this model which makes it very convenient. The performance of the CTrop and GPT3 models are validated by radiosonde data. Validation results demonstrate that the CTrop model shows higher precision than that of the GPT3 model in all meteorological parameters. The improvements are 6%, 10%, 9%, and 13% for pressure, temperature, water vapor pressure and weighted mean temperature, respectively. Three different resolutions of the CTrop model are also developed based on the sliding window algorithm, for which users can Tables 3 and 4 indicate that as the spatial resolution of the CTrop model decreases, the accuracy gradually decreases. The Ctrop model also has better accuracy than the GPT3-5 model in lower spatial resolution. Regarding the water vapor pressure, the RMS of the Ctrop-5 model is 2.71 hPa, showing an improvement of 5% compared with the GPT3-5 model. Considering the pressure, the CTrop-5 model achieves accuracy improvements of 10% when compared to the GPT3-5 model. Regarding temperature, the bias of the CTrop-5 model appears as a negative value, while the GPT3-5 model displays a positive value. The RMS of the CTrop-5 model is 3.26 K, showing an improvement of 22% when compared with the GPT3-5 model. In terms of T m , the precision of the CTrop-5 model is 3.43 K, also presenting a slight improvement when compared with the GPT3-5 model. Figure 11 shows that all models display a large RMS of water vapor pressure in the southeast of China. The CTrop-5 model shows an accuracy improvement in the northeast of China when compared with the GPT3-5 model. Both models show large errors in the southeast of China, which is subject to subtropical monsoons. Water vapor is transported by monsoons from the sea to these areas, resulting in a higher error for modeling [12]. Figure 12 shows that the CTrop-5 model and the GPT3-5 model present positive and negative biases of pressure in the majority of China. The CTrop-5 model presents a slight improvement at a few stations in the northeast of China when compared with the GPT3-5 model. Figure 13 shows that all models have better accuracy of T m in the south of China than in the north of China. The CTrop-5 model presents an improvement in the west of China when compared with the GPT3-5 model. Figure 14 shows that the CTrop-5 model has better accuracy in the west of China than that of the GPT3-5 model. Consequently, the lower horizontal resolution version of the CTrop model shows better accuracy in estimating the tropospheric parameters than that of the GPT3 model.

Conclusions
In this work, the distribution characteristics of meteorological parameters are analyzed, and it is observed that the meteorological parameters exhibit major annual and semiannual periodic variations that are also related to latitude and longitude. Considering the spatial distribution and time-varying characteristics of the meteorological parameters, a refined regional empirical model (CTrop) based on the sliding window algorithm is developed for the estimation of tropospheric key parameters over China. Only the day of the year (DOY) and the location of the station are needed when applying this model, which makes it very convenient. The performance of the CTrop and GPT3 models are validated by radiosonde data. Validation results demonstrate that the CTrop model shows higher precision than that of the GPT3 model in all meteorological parameters. The improvements are 6%, 10%, 9%, and 13% for pressure, temperature, water vapor pressure, and weighted mean temperature, respectively. Three different resolutions of the CTrop model are also developed based on the sliding window algorithm, for which users can access model parameters of different resolutions for their requirements to reduce the amount of gridded data provided to the users, as well as to speed up the troposphere delay computation process. The results demonstrate that the CTrop model realizes better performance than that of the GPT3 model over China, which can meet the needs for GNSS meteorology and GNSS positioning. Although the CTrop model is only developed in regions over China with the aim of meeting the needs for GNSS positioning in China, we intend to expand the model to a global model in future work.