remote sensing A Global Conversion Factor Model for Mapping Zenith Total Delay onto Precipitable Water

: The conversion factor is a key parameter for converting zenith wet delays (ZWD) into precipitable water vapour (PWV) with a mean value of 0.15, and the traditional method of calculating it is to model the weighted average temperature in the process of conversion factor calculation. Here, we overcome the dependence on high-precision atmospheric weighted average temperature for mapping ZWD onto PWV and build a global non-meteorological parametric model for conversion factor G Π model by using the gridded data of global conversion factor time series from 2006 to 2013 provided by the Global Geodetic Observing System (GGOS) Atmosphere. Internal and external accuracy tests were performed using data from four times (UTC 00:00, 06:00, 12:00, and 18:00) per day throughout 2012 and 2014 as provided by the GGOS Atmosphere, and the statistical average root-mean-square (RMS) and mean absolute errors (MAE) on a global scale are 0.0031/0.0026 and 0.0030/0.0026, respectively, which only account for 1.5–2% of the conversion factor value. In addition, the observed GPS data are also used to validate the established G Π model, and the RMS of the PWV differences between the established model and the observed meteorological data was less than 3.2 mm. The results show that the established G Π model has a high accuracy, which can be used to calculate the PWV value where no observed meteorological parameters are available.


Introduction
It is important to understand the spatiotemporal variability of water vapour so as to study global climate change and for rainfall forecasting [1]. The traditional method of monitoring atmospheric water vapour mainly uses sounding balloons, water vapour radiometers, etc., but the data from the traditional method (such as the sounding balloons) have low temporal resolution, e.g., two or four times daily, and larger spatial resolution, e.g., 200 to 300 km [2,3]. In addition, water vapour radiometers are mainly affected by the weather and cannot work under rainfall or foggy weather conditions, and the instruments are unwieldy and expensive [4][5][6]. Therefore, traditional monitoring methods cannot meet the requirements of small-and medium-scale water vapour and weather monitoring tasks.
With the development of the global navigation satellite system (GNSS) technology, the monitoring of atmospheric water vapour by ground-based GNSS technology has become a new research hotspot. Satellite signals are delayed through the troposphere, including atmospheric dry delay and atmospheric wet delay. Askne and Nordius [4] first presented the approximate relationship between atmospheric moisture delay and precipitable water vapour (PWV) [1] and proposed a method of calculating atmospheric water vapour using ground-based GPS observations. Bevis et al. [7,8] first proposed the concept of GPS meteorology and presented basic methods for water vapour calculation based on satellite observations. The signal delay of the satellite to the ground station can be used to calculate atmospheric water vapour [9][10][11][12]. Compared to traditional methods, the GNSS technique has many advantages for atmospheric precipitation water vapour calculation: high spatial and temporal resolution, low cost, global coverage, all-weather operability, real-time data analysis, etc. Therefore, GNSS meteorological research has great practical value and development potential [13][14][15].
To obtain the PWV value, the zenith total delay (ZTD) and the zenith hydrostatic delay (ZHD) should be calculated, and then the zenith wet delay (ZWD), which is typically estimated as an unknown parameter in GNSS processing, can be extracted from the ZTD by subtracting the ZHD. Therefore, the PWV can be calculated by multiplying the conversion factor by the ZWD [16,17]. The conversion factor is mainly affected by the weighted average temperature, which is why many studies have been conducted on the weighted average temperature in the past. Ross and Rosenfeld [18] used radiosonde data from 13 stations over two years to simulate the equation of linear regression between the weighted average temperature T m and the land surface temperature T s in America [18]. The linear relationship offers good precision in mid-latitude areas; however, the equation established above has its limitations because not all GNSS stations are equipped with meteorological sensors [19,20]. Ross et al. [20] found that T m and T s are influenced by many factors such as location and season by analysing the 23-year data from 53 radiosonde stations around the world. Yao et al. [21] used global data from 135 radiosonde stations from 2005 to 2009 to establish the first-generation global weighted average temperature model GT m -I. In 2013, the secondgeneration global weighted average temperature model GT m -II was established which overcomes the data loss problem over sea areas when using the GPT model and the T m~Ts relationship proposed by Bevis [3,21].
Here, we overcome the dependence on the high-precision atmospheric weighted average temperature in the process of the ZWD conversion to the PWV and establish a global non-meteorological parametric model for the conversion factor. The values of T m used to build the model were provided by the GGOS for the period from 2008 to 2013, with a time resolution of 6 hours (UTC 00:00, 06:00, 12:00, and 18:00) and a spatial resolution of 2.5 • × 2 • .

Method of Establishing the GΠ Model
The specific expression for calculating the PWV is as follows: where Π is the conversion factor, k 2 and k 3 represent the atmospheric refraction constant with values of 22.1 k/mb and 373,900 k 2 /mb , respectively [7], R V is a constant with a value of 461.495 J · kg −1 · K −1 , and ρ w is the water vapour density (g/cm 3 ). T m represents the weighted average atmospheric temperature (unit: K). A linear relationship between the weighted average temperature and ground temperature was proposed by Bevis et al. [7,22,23] using radiosonde data over many years: where T s represents the ground temperature (unit: K) and the conversion error from the ZWD to the PWV is approximately 2% based on the above formula [7]. Π varies with season and latitude and has an approximate value of 0.15. According to Equation (1), the accuracy of atmospheric precipitation data is determined by the accuracy of the ZWD and Π [24]. The international GNSS service (IGS) can provide the ZTD with an RMS error of 1.5~5 mm while the ZHD can be estimated with an RMS error of a few millimetres by means of an empirical model using the surface-measured pressure. Therefore, a high- precision ZWD value can be obtained from GNSS observations with an RMS error of a few millimetres [25,26]. Once the ZWD is obtained, Π is the only key parameter to calculate the PWV; therefore, determining Π is a key issue in GNSS water vapour remote sensing [27].

Methods of Obtaining the Conversion Factor
As mentioned above, Π is related to T m , while T m is related to the longitude and latitude, height of the station, and DoY (day of year); in addition, the given gas constant and atmospheric refraction constants are changed for different regions and times. Generally, there are several methods of obtaining the conversion factor. For the first method, Π can be regarded as a constant, usually taking a value of approximately 0.15 with a maximum range of 15%, which means 1 mm of the PWV corresponds to approximately 6.7 mm of the ZWD [28]. The above method is the simplest way to obtain Π, but the accuracy of the estimated PWV is low. Therefore, the formula T m = 70.2 + 0.72T s or the local linear regression formula is usually used to calculate T m , and then further used to calculate the conversion factor, which is considered as the second method.
In addition, Π can also be considered as a function of latitude, DoY, and expressed as follows [29]: where θ represents the latitude of the station, t D represents the DoY, a 0 = 5.861, a 1 = 0.011, a 2 = 0.054, and a 3 = 0.138; however, the above methods are not very accurate to express the conversion factor changes in different areas. Here, a global non-meteorological parametric model of the conversion factor is proposed based on the analysis of the factors influencing the conversion factor.

Analysis of the Conversion Factor
The time series of the global conversion factor data was analysed to determine the importance of those factors influencing the conversion factor. In our following analysis, the conversion factor was obtained from the GGOS Atmosphere with the temporal resolution of four times daily at UTC 00:00, 06:00, 12:00, and 18:00, respectively.

1.
Relationship between the conversion factor and the annual change.
The conversion factors at four specific locations (N30 • E120 • , N30 • W120 • , S30 • E120 • , and S30 • W120 • ) were calculated using the 8-year data provided by the GGOS from 2006 to 2013, and Figure 1a-d shows the relationship between the conversion factor and the time in the four locations, respectively. It can be observed from Figure 1 that the conversion factor exhibits periodic changes with time, which suggests that the annual cycle is an important factor affecting the conversion factor; therefore, it should be considered when establishing the conversion factor model. The detected dependences of the conversion factor are expected since they reflect the dependences of the temperature.

2.
Relationship between the conversion factor and the semiannual change.
It is shown in Figure 1 that the conversion factor has a strong periodic relationship with the annual cycle. Therefore, the following experiments were carried out to determine whether or not the conversion factor is relative to the semiannual variation. The relationship between the conversion factor and the semiannual change was analysed using data

2.
Relationship between the conversion factor and the semiannual change.
It is shown in Figure 1 that the conversion factor has a strong periodic relationship with the annual cycle. Therefore, the following experiments were carried out to determine whether or not the conversion factor is relative to the semiannual variation. The relationship between the conversion factor and the semiannual change was analysed using data from the year 2013.
, and the X axis refers to the DoY of 1 to 365, 2013. As can be seen from Figure 2, the conversion factor exhibits semiannual changes, which indicates that the semiannual change is another influential factor. Like the annual cycle, the semiannual cycle should also be considered when establishing the conversion factor model. from the year 2013. Figure 2a-d shows the changes of the conversion factor in four positions (N30° E120°, N30° W120°, S30° E120°, and S30° W120°), and the X axis refers to the DoY of 1 to 365, 2013. As can be seen from Figure 2, the conversion factor exhibits semiannual changes, which indicates that the semiannual change is another influential factor. Like the annual cycle, the semiannual cycle should also be considered when establishing the conversion factor model.

3.
Relationship between the conversion factor and the elevation. To explore whether or not the conversion factor is related to elevation, the calculated conversion factor for the year 2013, on a global scale, was analysed for the effect of change in elevation. Figure 3 shows the conversion factor changes with altitude, and the trend is approximately linear at 0-4 km, while the trend disappears above 4 km. Due to the fact that most regions in the world lie at an elevation of less than 4 km, it is necessary to consider the elevation factor when establishing the conversion factor model [30]. In addition, some points well outside the linear relationship may be related to temperature, which remains to be further verified. It should be noted that not very many regions on Earth exist at these altitudes, and they correspond to dry regions of very low ZWD values with also low variability.

3.
Relationship between the conversion factor and the elevation. To explore whether or not the conversion factor is related to elevation, the calculat conversion factor for the year 2013, on a global scale, was analysed for the effect of chan in elevation. Figure 3 shows the conversion factor changes with altitude, and the trend approximately linear at 0-4 km, while the trend disappears above 4 km. Due to the fa that most regions in the world lie at an elevation of less than 4 km, it is necessary to co sider the elevation factor when establishing the conversion factor model [30]. In additio some points well outside the linear relationship may be related to temperature, whi remains to be further verified. It should be noted that not very many regions on Ear exist at these altitudes, and they correspond to dry regions of very low ZWD values wi also low variability.

4.
Analysis of the modelling grid division. Data with a spatial resolution of 2.5° × 2° provided by the GGOS Atmosphere we used to model the conversion factor; however, whether the conversion factor has a re tionship with the latitude and longitude has yet to be determined. Therefore, it is nece sary to analyse the correlation between the conversion factor and the latitude/longitu so as to ensure correctness of the model with regard to the adopted grid division metho Figures 4 and 5 show the relationship between the conversion factor and the longitude/la itude, respectively, based on the conversion factor obtained from different longitudes an latitudes in 2013. It can be seen from Figures 4 and 5 that the conversion factor has a stro relationship with the latitude and longitude; therefore, it is reasonable to establish t conversion factor using the grid division method. In addition, some outliers are observe mainly due to the majority of elevations above 4000-5000 m (mainly the Himalayan r gion) being located in the latitude band.

4.
Analysis of the modelling grid division.
Data with a spatial resolution of 2.5 • × 2 • provided by the GGOS Atmosphere were used to model the conversion factor; however, whether the conversion factor has a relationship with the latitude and longitude has yet to be determined. Therefore, it is necessary to analyse the correlation between the conversion factor and the latitude/longitude so as to ensure correctness of the model with regard to the adopted grid division method. Figures 4 and 5 show the relationship between the conversion factor and the longitude/latitude, respectively, based on the conversion factor obtained from different longitudes and latitudes in 2013. It can be seen from Figures 4 and 5 that the conversion factor has a strong relationship with the latitude and longitude; therefore, it is reasonable to establish the conversion factor using the grid division method. In addition, some outliers are observed, mainly due to the majority of elevations above 4000-5000 m (mainly the Himalayan region) being located in the latitude band.

Expression of the Conversion Factor Model
As described above, the conversion factor is affected by annual and semiannual effects and varies with height; thus, based on the above analysis, the following model was proposed: Π= a 1 + a 2 h + a 3 cos 2π doy-C 1 365.25 + a 4 cos 4π doy-C 2 365.25 + a 5 cos 2π hod-C 3 24 where a 1 represents the average value of the conversion factor, a 2 refers to the height vertical gradient, h is the station height, doy is day of year, hod is hour of day (UTC time), a 3 , a 4 , a 5 represent annual, semiannual, and day period amplitudes, respectively, while C 1 , C 2 , C 3 represent annual, semiannual, and day period initial phases, respectively. The model coefficients can be estimated by using the least squares method.

Expression of the Conversion Factor Model
As described above, the conversion factor is affected by annual and semiannual effects and varies with height; thus, based on the above analysis, the following model was proposed:

Expression of the Conversion Factor Model
As described above, the conversion factor is affected by annual and semiannual effects and varies with height; thus, based on the above analysis, the following model was proposed:  The model coefficients were estimated in this paper using the grid data provided by the GGOS Atmosphere for 2008-2013. Figure 6a-d show the coefficients a2, a3, a4, and a5, respectively, while coefficient a1 is a constant with the value of about 0.0037. As can be seen from Figure 6a, the absolute value of a2 for the global conversion factor model is smaller on land and larger in marine areas while the absolute value of a3 is smaller at lower latitudes and larger at high latitudes except for some abnormal values near the equator (see Figure 6b); it can be seen in Figure 6c that the absolute value of the semiannual amplitude is larger in polar regions and smaller in the vicinity of the equator. In addition, the semiannual absolute amplitude appears to be larger in the northern hemisphere than that in the southern hemisphere, and the largest values occur over Siberia. As can be seen from Figure 6d, the absolute value of the diurnal amplitude is relatively large on land when compared to that over sea areas. In addition, opposite values appear in the eastern hemisphere (positive values) and western hemisphere (negative values): the reason for this may be that the diurnal amplitude is closely related to the daily temperature, so the opposite trend appears in the eastern and western hemispheres.

Result
Here, the root-mean-square (RMS) and mean absolute errors (MAE) are used as the standard for evaluating the established model as follows: <!--MathType@Translator@5@5@MathML2 (no namespace).tdl@MathML (5) Figure 6. Coefficients distribution of the established global conversion factors globally estimated using the least squares method, where (a-d) represent the distribution of coefficients a2, a3, a4, and a5, respectively.

Result
Here, the root-mean-square (RMS) and mean absolute errors (MAE) are used as the standard for evaluating the established model as follows: where x i represents the estimated value of the GΠ model while x 0 i refers to the true value calculated using the GGOS data and N represents the total number of observed values; MAE can be used to embody the degree of deviation of the estimated value from the true value, while RMS is used to measure the reliability of the established model. In the following Sections 3.1 and 3.2, the internal accuracy testing is performed for the period of the data used in the determination of the conversion factor model while the external accuracy is carried out using the GGOS data over a period not used in this analysis.

Internal Accuracy Testing
To validate the internal accuracy of the established model, the conversion factor derived from the GΠ model four times daily (UTC 00:00, 06:00, 12:00, and 18:00) in 2012 was compared to those values globally derived from the GGOS data with a spatial resolution of 2.5 • × 2 • . Figure 7 shows the internal accuracy of the (a) MAE and (b) RMS of the GΠ model while Table 1 lists the statistical results on a global scale. It can be seen from Table 1 that the estimated value of the conversion factor calculated using the GΠ model was closer to that derived from the GGOS data, and the magnitude of internal accuracy was 10 −3 . It also can be seen from Figure 7 that the RMS and MAE range from 0 to 0.0317 and from 0 to 0.026 on the global scale, which indicates a good accuracy in the GΠ model; however, some abnormal values appeared near the equator, which may have been affected by the El Niño event at that time [31].

External Accuracy Testing
To assess the external accuracy of the GΠ model, the conversion factor derived from the established GΠ model in 2004 was compared with that from the GGOS data at four times each day (UTC 00:00, 06:00, 12:00, and 18:00) on a global scale. Statistical results are shown in Table 2, while Figure 8 shows the distribution of the MAE and RMS. It also can be concluded from Table 3 and Figure 8 that the validated result is similar to that arising from internal testing.

External Accuracy Testing
To assess the external accuracy of the G<!--MathType@Translator@5@5@MathML2 (no namespace).tdl@MathML 2.0 (no namespace)@ --> <math> <semantics>  The compared result was further analysed to obtain the distribution of the RMS. Figure 9 shows the histogram of the average RMS of external accuracy testing for the GΠ model. It can be observed that most RMS values were less than 0.005, which suggests that the GΠ model provides a good model to the GGOS data.
To analyse the accuracy of the GΠ model at different heights, longitudes, and latitudes, the external accuracy of the MAE and RMS was calculated, and the statistical results are shown in Tables 3-5 and Figure 9. aabeqaamaabaabauaakeaacqqHGoauaaa@4080@ </annotation> </semantics> </math> <!--MathType@End@5@5@ --> model provides a good model to the GGOS data. To analyse the accuracy of the G<!--MathType@Translator@5@5@MathML namespace).tdl@MathML 2.0 (no namespace)@ --> <math> <semantics> <mi>&#x03A0;</mi> <annotation encoding='MathType-MTEF'>MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbWexLMBbXgBd9gzLbvyNv2CaeHbl7mZLdGeaGqipu0Je9sqqr pepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs  It can be observed from Table 3 that the average MAE and RMS of the GΠ model increased with height, which suggests that the conversion factor is influenced by altitude. Table 4 shows that the average MAE and RMS of the GΠ model presented an increasing trend with longitude and then decreased, while the values of the MAE and RMS were largest in the 90 •~1 20 • longitude. Table 5 shows that the average values of the MAE and RMS increased with latitude. In Figure 10a-f, a few larger MAE and RMS values are found at heights of up to 1000 m, within 60 • W to 60 • E longitudes and at 30 • N to 30 • S latitudes: this phenomenon may be caused by the El Niño event. In addition, the overall accuracy of the GΠ model is high at low latitudes and low longitudes and relatively poor over high latitude and high longitude ranges.

Case Study
To verify the accuracy of the established G<!--MathType@Translator@5@5@MathML2 (no namespace).tdl@MathML 2.0 (no namespace)@ --> <math> <semantics> <mi>&#x03A0;</mi> <annotation encoding='MathType-MTEF'>MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbWexLMBbXgBd9gzLbvyNv2CaeHbl7mZLdGeaGqipu0Je9sqqr pepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs 0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaai aabeqaamaabaabauaakeaacqqHGoauaaa@4080@ </annotation> </semantics> </math> <!--MathType@End@5@5@ --> model, comparisons for stations HKSC and ANJI in the CORS network of Hong Kong and Zhejiang Province were conducted. THKSC and ANJI were selected because the large PWV values existed for the whole year at those stations, which can reflect the real performance of proposed model. In addition, the collocated meteorological stations existed at those two stations, and the high-precision PWV could be obtained as the reference. The data for the whole year of 2014 from HKSC and for the period from 1 May to 20 July 2015 from the ANJI site were selected. The specific location and height of the two

Case Study
To verify the accuracy of the established GΠ model, comparisons for stations HKSC and ANJI in the CORS network of Hong Kong and Zhejiang Province were conducted. THKSC and ANJI were selected because the large PWV values existed for the whole year at those stations, which can reflect the real performance of proposed model. In addition, the collocated meteorological stations existed at those two stations, and the high-precision PWV could be obtained as the reference. The data for the whole year of 2014 from HKSC and for the period from 1 May to 20 July 2015 from the ANJI site were selected. The specific location and height of the two selected stations are given in Table 6. The high-precision PWV data were calculated using GPS observations and the observed meteorological data based on the method proposed in Section 2 (met-derived PWV). In our experiment, the precise pointing positioning (PPP) technique [32,33] was used to process the GPS observations using the conventional ionosphere-free (IF) code and phase measurements. The receiver coordinates, the tropospheric zenith wet delay, the receiver clock error, and the IF phase ambiguities were regarded as unknown parameters, which were estimated using the extended Kalman filter. The global pressure and temperature (GPT) model [34] was used to obtain the meteorological parameter (pressure) while the Saastamoinen model was introduced to calculate the tropospheric hydrostatic delay [35]. The corresponding meteorological data were checked before use. In addition, the PWV values were also calculated using the values of parameter Π estimated on the basis of the established model (model-derived PWV). It is worth mentioning that the final estimated parameter Π was based on the re-established model, which excludes the outliers mentioned in Section 3.2. Two kinds of PWV values were compared as shown in Figure 11, where (a) and (b) provide comparisons of the PWV time series. It can be observed that the model-derived PWV time series was consistent with the met-derived PWV time series. In addition, it can be concluded that the PWV difference increased with increasing PWV values. If the met-derived PWV is considered as a reference, the RMS of the model-derived PWV for the ANJI site is 2.95 mm while the STD and bias are 1 mm and 2.77 mm, respectively. For the HKSC station, the RMS of the model-derived PWV was 3.18 mm while the STD and bias were 1.6 mm and 2.75 mm, respectively. It can also be seen from Figure 10a,b that the PWV differences for the two sites were small when the PWV values were relatively low (i.e., below 40 mm) while the differences increased with increasing PWV. The analysis shows that there was a small systemic bias between the met-derived PWV and the model-derived PWV, which may be due to the failure to consider all factors when building the GΠ model. Generally, the GΠ model established in this research has a high precision and can be used for the calculation of the PWV in the absence of meteorological parameters.

Discussion
It can be observed from the internal and external validations that the proposed mode has a good performance on the global scale. However, some outliers appeared around the equator both in internal and external experiments, and such result may be related to the

Discussion
It can be observed from the internal and external validations that the proposed model has a good performance on the global scale. However, some outliers appeared around the equator both in internal and external experiments, and such result may be related to the GGOS data used for modeling and affected by the El Niño event at that time [31]. It can also be found that coefficient a3 of the established global conversion factors globally estimated using the least squares method (Figure 6) also shows the similar distribution as that of the MAE and RMS in internal and external validations; it means that the outliers are mainly affected by coefficient a3 of the established global conversion factors in Equation (4). Coefficient a3 is the annual phases of the global conversion factors highly related to the global climate change; therefore, the occurrence of the El Niño event may have an evident influence on coefficient a3.
Two stations in Zhejiang province and Hong Kong were used to further validate the proposed global conversion factor model, and the experimental result also shows the good performance of the proposed model. Although those two stations cannot fully evaluate the performance of the proposed global conversion factor model, the selected two stations have relatively high PWV values with the average PWV of approximately 40 mm. Due to this paper being mainly focused on the establishment of the global conversion factor model, more attentions was paid to the modeling. In addition, internal and external accuracy validations of the proposed model on a global scale using the GGOS data were performed, which proved the good performance of the proposed model. Therefore, only the case with two stations was used to perform the application experiment.

Conclusions
A global conversion factor model of non-meteorological parameters (GΠ model) was established by using gridded data provided by the GGOS Atmosphere from 2006 to 2014. The established model considers the changes of the annual, semiannual, and day atmospheric water vapour; in addition, this model is also related to longitude, latitude, and height. The internal accuracy of the established model was tested using the data from 2012 provided by the GGOS Atmosphere, and the result shows that the average values of the MAE and RMS of the GΠ model were 0.0026 and 0.0031 respectively. The external accuracy was also tested by using data provided by the GGOS in 2014, and the result shows that the average values of the MAE and RMS of the established GΠ model were 0.0026 and 0.0030, respectively. Those validated results reveal that the established GΠ model based on the GGOS Atmosphere data offers high precision, and the internal and external accuracies are of the order of magnitude of 10 −3 , which accounts for only 1.5-2% of the conversion factor value. To validate the applicability of the established GΠ model, a case study was assessed that involved two sites in different regions with large PWV values for the whole year over different time periods. The statistical results from both sites revealed that the RMS error of the PWV calculated based on the established GΠ model was about 3 mm; in addition, the model coefficient could be easily obtained at a specific grid point and the conversion factor was directly obtained, which was of significance for the calculation of the PWV, especially when faced with the unavailability of meteorological parametric data.