Development of Time-dependent mean Temperature Equations for GPS Meteorology
Development of Time-dependent mean Temperature Equations for GPS Meteorology
Journal of Positioning, Navigation, and Timing. 2014. Dec, 3(4): 143-147
Copyright © 2014, The Korean GNSS Society
  • Received : September 26, 2014
  • Accepted : October 10, 2014
  • Published : December 15, 2014
Export by style
Cited by
About the Authors
Jihyun Ha

The mean temperature is one of the key parameters in computing Precipitable Water Vapor (PWV) from Global Positioning System (GPS) measurements and is usually derived as a function of surface temperature through the use of a mean temperature equation (MTE). In this study, two new types of MTEs were developed as functions solely of the observation time so that the mean temperature can be obtained without surface temperature measurements. To validate the new models, we created one-year time series of GPS-derived PWV using the new MTEs and compared them with the radiosonde-observed PWV. The bias and root-mean-square error were on the other of ~1 mm and ~2 mm, respectively.
Global Positioning System (GPS) satellite signals are delayed due to water vapor in the troposphere of the Earth before reaching GPS receivers on the ground. Based on GPS measurements, information on the water vapor distribution in the atmosphere can be obtained by computing tropospheric delays of the GPS signal. GPS meteorology provides the amount of Precipitable Water Vapor (PWV), Slant Wet Delay (SWD) and Zenith Wet Delay (ZWD). Such GPS-based meteorological information, if assimilated as a predicator or a climatic factor in Numerical Weather Prediction (NWP), can be utilized to improve the performance of NWP models ( Kuo et al. 1998 ).
Recently, to improve the accuracy of GPS-derived PWV (this quantity will be referred to as GPS PWV hereafter), a variety of studies are being conducted to develop optimized models for local meteorological conditions. Regional atmospheric model developments for GPS meteorology are usually concerned with a Mean Temperature Equation (MTE) and refraction constants. Especially, the MTE is a key parameter in translating ZWD into PWV. Bevis et al. (1992) suggested that the mean temperature of the atmosphere can be approximated by the surface temperature. The Bevis model ( Bevis et al. 1992 ) was developed using radiosonde measurements taken at 13 sites in North America. Emardson & Derks (2000) developed a model more suited to the climatic conditions of Europe. Liou et al. (2000) developed a local MTE model in addition to monthly MTEs, based on meteorological observations in Taiwan. Woo (2003) and Ha et al. (2006) developed local MTEs based on radiosonde observations in South Korea. All the MTE models mentioned in this paragraph are given as linear functions of the surface temperature and thus temperaturemeasuring devices should be co-located at the GPS site
In this study, we propose two new types of mean temperature models, in which surface temperature measurements are unnecessary. The novel concepts behind the newly developed MTE models are introduced in the following section and validation results through the error analysis of GPS PWV values obtained with new models are described in Section 3.
Through the GPS data processing, the amount of the GPS signal delay caused by water vapor can be obtained as ZWD. Then, the estimated ZWD should be converted to PWV by Eq. (1).
PPT Slide
Lager Image
PPT Slide
Lager Image
In Eq. (2), ρ is the density of liquid water (998.00879 kg/m 3 ), and Rv is the specific gas constant of water vapor (4.615 × 10 2 J/kg‧K). k3 and
PPT Slide
Lager Image
are refraction constants and their values are (3.776 ± 0.004)×10 5 K 2 /mbar and (17 ± 10) K/mbar, respectively ( Davis et al. 1985 ). Tm is the mean temperature of the atmosphere. The only variable in Eq. (2) is Tm , thus it is a key component in deriving GPS PWV from ZWD estimates. Tm is usually expressed as a function of surface temperature Ts . The most widely used MTE is the Bevis model, which is represented as Tm = 0.72 Ts + 70.2 ( Bevis et al. 1992 ).
The mean temperature can be obtained from the radiosonde measurement using the observed temperature and the partial pressure of water vapor ( Davis et al. 1985 ). Radiosonde observations used in this study span eight years from 1998 through 2005 and they are taken at seven domestic radiosonde stations in South Korea. Fig. 1 shows the mean temperature variations with respect to the year. The relation between Tm and observation time tY ( tY stands for yearly time scale) was estimated as a sinusoidal function with a linear trend, and the result is depicted in Eq. (3).
PPT Slide
Lager Image
Mean temperature variations with respect to the observation time in a yearly scale for the years 1998 through 2005.
PPT Slide
Lager Image
The RMS between Tm computed by Eq . ( 3 ) and radiosonde observations is 4.2 K. With Eq. (3), one can get the mean temperature just by plugging in the observation time as a yearly time scale.
As another method of deriving time-dependent MTEs, the observed mean temperatures against observation time tM ( tM stands for monthly time scale) were plotted in Fig. 2. In the figure, the horizontal axis is tM ; 0.0 corresponds to January 1 st and 1.0 December 31 st . Then the second-, third-, fourth-, and sixth-order polynomials were tried to find out the best fitting function. When the four fitted solutions were compared with radiosonde measurements, the RMS differences were in the range of 4.1-4.7 K. The best agreement was achieved for the fourth-order polynomial and the exact relation is shown in Eq. (4).
PPT Slide
Lager Image
Mean temperature variations with respect to the observation time in a monthly scale for the years 1998 through 2005.
PPT Slide
Lager Image
Hereafter, the two MTEs shown in Eqs. (3) and (4) will be referred to as HPt Y and HPt M , respectively.
For the error analysis, the mean temperatures observed from radiosondes and those predicted by the Bevis, HPt Y and HPt M models are plotted together in Fig. 3. In comparison with the time-dependent MTE models, mean temperature differences of the Bevis model show clear seasonal signals. For the cases of the HPt Y and HPt M models, the scatters of mean temperature differences are larger than the Bevis model, with the largest discrepancies periodically happening in the winter months.
PPT Slide
Lager Image
The left hand sides are mean temperatures observed (black solid line) and calculated (red dots) with each of three MTE models (Bevis, HPtY, and HPtM) during the years 1998 through 2005. The right hand sides are the corresponding differences.
The error statistics of the three MTE models with respect to the radiosonde measurements are listed in Table 1 . In the table, mean differences are negative, implying that mean temperatures computed by the three MTEs are larger than the radiosonde observation in the average sense. The HPt Y model has a zero bias, but its variance shows the largest value. In the case of the Bevis model, the variance and accumulated error are the lowest. Also, the Bevis model produces ~5% smaller errors in the range of >4 K differences. Thus we conclude that MTEs as functions of surface temperatures are slightly better than time-dependent ones. However, in the following section, we will show that the performance of the new models is comparable to that of the Bevis model by computing one-year time series of GPS PWV and comparing them with radiosonde-observed PWV.
The mean, variance, and accumulative Tmerror distribution over 10K, 4K, 3K, 2K and 1K differences of the Bevis, HPtY, and HPtMMean Temperature Equation models with respect to the radiosonde observation.
PPT Slide
Lager Image
The mean, variance, and accumulative Tm error distribution over 10K, 4K, 3K, 2K and 1K differences of the Bevis, HPtY, and HPtM Mean Temperature Equation models with respect to the radiosonde observation.
For validation of the HPt Y and HPt M model, PWV estimates were obtained based on GPS and radiosonde measurements taken at the Sokcho weather station for one year from March 1, 2008 to February 28, 2009. The Sokcho weather station located in the eastern Korea has a permanent GPS site and conducts radiosonde observations twice a day at 00:00 and 12:00 UTC. Thus, it is a perfect site to check the performance of GPS PWV against the radiosonde measurement. For the GPS data processing, GIPSY-OASIS II developed by Jet Propulsion Laboratory (JPL) was used ( Webb & Zumberge 1993 ). The 24-hour GPS data were processed with the standard precise point positioning technique by utilizing the JPL-provided precise orbit and clock solutions ( Zumberge et al. 1997 ). The analysis results are shown in Fig. 4 , where PWV values at 00:00 and 12:00 UTC are shown.
PPT Slide
Lager Image
(a) Radiosonde PWV (blue solid line) and GPS PWV (red dots) using the Bevis model at the Sokcho weather station from March 1, 2008 to February 28, 2009, (b) PWV differences between radiosonde PWV and GPS PWV values obtained with three different MTE models.
Fig. 4a depicts GPS PWV computed with the Bevis model, together with radiosonde PWV. Even though radiosonde and GPS PWV show the same tendency, most of the radiosonde PWV values are slightly larger than GPS estimates as can be seen in Fig 4b . Fig. 4b shows PWV differences between GPS PWV obtained with three MTE models and radiosonde PWV. From Fig. 4b , one can find that the differences are very similar regardless of the choice of the MTE model, and most of the differences are lying in the range of ±0.5 mm.
The biases between radiosonde PWV with GPS PWV, which were computed with three different MTE models, are almost at the same level. The biases for the Bevis, HPt Y and HPt M models are -1.0 mm, -1.0 mm, and -1.1 mm, respectively. The RMS values of GPS PWV computed with the Bevis, HPt Y and HPt M models are 2.4 mm, 2.3 mm, and 2.4 mm, respectively. In addition to these results, we found that the RMS and bias errors obtained in this test are comparable to previous studies of validating GPS PWV with radiosonde PWV (e.g., Liou et al. 2000 , Niell et al. 2001 ). Thus we conclude that one can replace the linear regression type MTE, which is a function of the surface temperature, with the time-dependent one without losing the accuracy and precision of GPS PWV.
Local MTE models (HPt Y and HPt M ) as functions of the observation time were developed using the eight-year-long radiosonde measurements collected at seven radiosonde sites in South Korea. The accuracies of local models and the Bevis model were analyzed with radiosonde observations assumed as the truth. The errors of the Bevis model showed seasonal signals with high errors in the summer and winter seasons. For the HPt M and HPt Y models, the errors were higher only during the winter months. For validation of the newly developed models, we compared radiosonde PWV and GPS PWV for one year. Radiosonde PWV and GPS PWV computed with three MTE models show similar trends, and the biases are ~1 mm. Also the RMS values were almost the same at ~2 mm. This result indicates that, if a locally optimized time-dependent MTE is available, it is possible to obtain accurate and precise PWV without surface temperature measurements.
This research was supported by a grant (14NSIP-B081051-01) from National Land Space Information Research Program funded by Ministry of Land, Infrastructure and Transport of Korean government.
Jihyun Ha received the Doctor’s degree in Geoinfomatic Engineering from Inha University in 2009. She was a senior researcher of Spatial Information Research Institute, Seoul, Korea, since 2012. Her research interests include GNSS, precise positioning, and GPS meteorology.
Bevis M. , Businger S. , Herring T. A. , Rocken C. , Anthes R. A. 1992 GPS meteorology: Remote sensing of atmospheric water vapor using the global positioning system JGR 97 (D14) 15784 - 15801    DOI : 10.1029/92JD01517
Davis J. L. , Herring T. A. , Shapiro I. I. , Rogers A. E. E. , Elgered G. 1985 Geodesy by ratio interferometry: Effects of atmospheric modeling errors on estimates of baseline length Radio Science 20 1593 - 1607    DOI : 10.1029/RS020i006p01593
Emardson T. R. , Derks H. J. P. 2000 On the relation between the wet delay and the integrated precipitable water vapour in the European atmosphere Meteorological Applications 7 61 - 68    DOI : 10.1017/S1350482700001377
Ha J. , Park K.-D. , Heo B.-H. 2006 Development of a local mean temperature equation for GPS-based precipitable water vapor over the Korean peninsula JASS (in Korean) 23 373 - 384    DOI : 10.5140/JASS.2006.23.4.373
Kuo Y.-H. , Zou X. , Chen S. J. , Huang W. , Guo Y.-R. 1998 A GPS/MET sounding through and intense upperlevel front Bulletin of the American Meteorological Society<0617:AGMSTA>2.0.CO;2 79 617 - 626    DOI : 10.1175/1520-0477(1998)079<0617:AGMSTA>2.0.CO;2
Liou Y.-A. , Huang C.-Y. , Teng Y.-T. 2000 Precipitable water observed by ground-based GPS receivers and microwave radiometry Earth Planets and Space 52 445 - 450    DOI : 10.1186/BF03352256
Niell A. E. , Coster A. J. , Solheim F. S. , Mendes V. B. , Toor P. C. 2001 Comparison of measurements of atmospheric wet delay by radiosonde, water vapor radiometer, GPS, and VLBI Journal of Atmospheric and Oceanic Technology<0830:comoaw>;2 18 830 - 850    DOI : 10.1175/1520-0426(2001)018<0830:comoaw>;2
Webb F. H. , Zumberge J. F. 1993 An Introduction to the GIPSY/OASIS-II JPL Publication California
Woo S.-H. 2003 Adjustment of meteorological variables for the accurate estimation of GPS PWV, MS Thesis Seoul National University (in Korean)
Zumberge J. F. , Heflin M. B. , Jefferson D. C. , Watkins M. M. , Webb F. H. 1997 Precise point positioning for the efficient and robust analysis of GPS data from lager networks JGR 102 (B3) 5005 - 5017    DOI : 10.1029/96JB03860