Advanced
Effects of Tropospheric Mapping Functions on GPS Data Processing
Effects of Tropospheric Mapping Functions on GPS Data Processing
Journal of Astronomy and Space Sciences. 2010. Mar, 27(1): 21-30
Copyright ©2010, The Korean Space Science Society
  • Received : December 12, 2009
  • Published : March 15, 2010
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Ji-Hye Won
Dept. of Geoinformatic Engineering, Inha University, Incheon 402-751, Korea
jhwon1021@inha.edu
Kwan-Dong Park
Dept. of Geoinformatic Engineering, Inha University, Incheon 402-751, Korea
Ji-Hyun Ha
Dept. of Geoinformatic Engineering, Inha University, Incheon 402-751, Korea
Jung-Ho Cho
Korea Astronomy and Space Science Institute, Daejeon 305-348, Korea
Abstract
In processing space geodetic data, mapping functions are used to convert the tropospheric signal delay along the zenith direction to the line of sight direction. In this study, we compared three mapping functions by evaluating their effects on the tropospheric signal delay and position estimates in GPS data processing. The three mapping functions tested are Niell Mapping Function (NMF), Vienna Mapping Function 1 (VMF1), and Global Mapping Function (GMF). The tropospheric delay and height estimates from VMF1 and GMF are compared with the ones obtained with NMF. The differences among mapping functions show annual signals with the maximum occurring in February or August. To quantitatively estimate the discrepancies among mapping functions, we calculated the maximum difference and the amplitude using a curve fitting technique. Both the maximum difference and amplitude have high correlations with the latitude of the site. Also, the smallest difference was found around 30 o N and the amplitudes increase toward higher latitudes. In the height estimates, the choice of mapping function did not significantly affect the vertical velocity estimate, and the precision of height estimates was improved at most of the sites when VMF1 or GMF was used instead of NMF.
Keywords
1. Introduction
Space geodetic technologies such as Global Navigation Satellite System (GNSS) and Very Long Baseline Interferometry (VLBI) correct for tropospheric delay errors to obtain more precise data processing results by using empirical or a priori models (Ha et al. 2005). And mapping functions are used to convert tropospheric delays along the zenith direction to those of the slant direction toward GNSS satellites and a VLBI point source. The relation between zenith and slant delays can be described as follows:
Lager Image
In the Equation (1), the Slant Total Delays (STD) at elevation angel e are made up of the Zenith Hydrostatic Delay (ZHD) and the Zenith Wet Delay (ZWD) with corresponding mapping functions m h ( e ) or m w ( e ) multiplied.
The names and locations of GPS stations used in this study.
Lager Image
The names and locations of GPS stations used in this study.
Lager Image
Map of GPS sites used in this study.
The first kinds of mapping functions were suggested by Saastamoinen (1972) and Marini (1972), and continuously developed by Lanyi (1984), Davis et al. (1985) and Ifadis (1986). Then, the Marini mapping function (Marini 1972) was improved into the Herring mapping function (Herring 1992) and the Niell Mapping Function (Niell 1996), which used to be the most widely used one up to recently. NMF was built on radiosonde profiles from the northern hemisphere only, so that latitudedependent biases are largest in high southern latitudes (Boehm et al. 2006a). However, lately developed mapping functions have covered this deficiency by using Numerical Weather Models (NWM) provided by European Centre for Medium-range Weather Forecasts (ECMWF). The representative mapping functions based on NWM are the Vienna Mapping Function (VMF; Boehm & Schuh 2003), the Vienna Mapping Function1 (VMF1; Boehm et al. 2006b), and the Global Mapping Function (GMF; Boehm et al. 2006a).
Some previous researches showed that replacing mapping functions from the NMF to state-ofthe- art mapping functions such as VMF1 and GMF not only affect on tropospheric delay estimates but also influence vertical position estimates. Boehm et al. (2006a) analyzed the height changes by testing NMF, GMF, and VMF1 in GPS analysis, and found that the agreement between VMF1 and GMF is very good, but NMF had the largest difference in south of 45 o S, and in the northeastern China and Japan, with station height differences up to 10 mm. As well, Munekane et al. (2008) proved that it is possible to remove annual signals from the vertical GPS time-series by correctly modeling tropospheric delays, and that the precision of height time-series does improve when the latest mapping function was utilized.
In this paper, we analyzed the variation of tropospheric delays and vertical deformations caused by different mapping function in GPS data processing. To quantitatively estimate effects of tropospheric mapping functions, we used nine GPS stations located in the Korean peninsula and surrounding areas in the range of 20 o N - 40 o N and 110 o E - 145 o E. We processed the GPS data collected at the selected stations by using GIPSY-OASIS II for four years from January 2005 through December 2008. The nine permanent GPS stations used in this study are listed in Table 1 with their locations denoted in Figure 1 .
2. Mapping Functions
The three mapping functions tested in this study commonly use the Herring type of continued fraction as shown in Equation (2), and each mapping has a different way of determining its coefficients (Herring 1992). The mapping function m i ( e ) at a specific elevation angle e can be presented as f ( e ; a i ; b i ; c i ). Here i can be marked as h for the hydrostatic mapping function, and as w the wet mapping function.
Lager Image
For NMF, VMF1, and GMF, the coefficients a , b , and c are determined from ray-tracing techniques. Input parameters for the ray-tracing are elevation angle, site height, air pressure, temperature and water vapor pressure at distinct pressure levels in the neutral atmosphere. Based on those values, the following can be calculated: refractivity at each altitude; hydrostatic and wet delays along the zenith direction; and hydrostatic and wet delays plus bending effect along the slant direction. A mapping function is defined by substituting zenith delays and slant delays produced by ray-tracing for the left side of the Equation (2), and by deciding the coefficients a i , b i , and c i based on leastsquares adjustment (Boehm & Schuh 2003). In order to apply ray-tracing, NMF used radiosonde data, while VFM1 and GMF utilized numerical models. Characteristics of the three mapping functions are explained the following section.
- 2.1 NMF
NMF sets its coefficients by applying ray-tracing to nine elevation angles from 3 o to 90 o on radiosonde data, which had been collected for one year between 1987 and 1988 (Niell 1996). The coefficients are provided only at four latitudes; 15, 30, 45, and 60 o N. Thus, for a site not located at those four latitudes, the coefficients should be determined by linear interpolation. The coefficients of the wet mapping function are defined as a constant by latitude, and those of the hydrostatic mapping function are defined as a function of latitude and observation time as follows (Niell 1996):
Lager Image
In the Equation (3), φ is the site latitude and Day-of-Year (DOY) is the date based on UT (Universal Time). a avg and a amp stand for the mean value and an amplitude, respectively; and they are given as constants.
The hydrostatic mapping function of NMF needs correction terms to make up for the heightdifference of observation site. Therefore, the hydrostatic mapping function of NMF should be redefinedas Equation (4) with correction terms:
Lager Image
In the equation, H is the height of the site (km), and the coefficients of correction terms a ht , b ht , and c ht are given as constants (Niell 1996).
- 2.2 VMF1
To determine its coefficients, VMF1 uses ray-tracing of NWM instead of radiosonde data. The coefficients a h and a w are calculated by using ray-tracing at elevation angle 3.3 o every six hours. Then, a predefined constant and formula are used for the b h and c h coefficients (Boehm et al. 2006b). For the coefficients b h and c h , the VMF1 calculated the monthly average data at 0, 6, 12 and 18 UT out of ERA40 (ECMWF Re-Analysis 40-years) in 2001. Based on the monthly average data, raytracing was applied to ten elevation angles from 3.3 o to 90 o . Through this procedure, the coefficient b h was defined as 0.0029. The coefficient c h was defined as in the following Equation (5) as a function of time and latitude:
Lager Image
The constant coefficients c 0 , c 10 , c 11 , and ψ are defined differently for the northern and southern hemisphere. The coefficients b w and c w of the wet mapping function are same as the NMF values at 45 o latitude.
- 2.3 GMF
VMF1 has a weakness that it has to directly apply ray-tracing by using a NWM in order to calculate the coefficients a h and a w , or has to download coefficients from a web site. To make up for this limitation, GMF defined the coefficients a h and a w as mathematical formulae. GMF applied ray-tracing to the monthly average values out of ERA40 from 1999 to 2002, and set the coefficients a h and a w by calculating ray-traced values for 36 months at 312 grid points around the world.
In Equation (6), the coefficient a is defined as a function of DOY (Boehm et al. 2006a). In the equation, a 0 indicates the mean, and A is the amplitude about a 0 . In both the hydrostatic and wet mapping function, a 0 and A were expanded into spherical harmonics up to degree and order nine (Boehm et al. 2006a).
Lager Image
3. Comparison of Mapping Functions
To analyze the effects on the tropospheric delay and vertical deformations caused by different mapping functions used in GPS data processing, a high-precision GPS data processing program called GIPSY-OASIS II was used (Webb & Zumberge 1993). Except for mapping functions, the same error correction model and processing techniques were applied. In the data processing, the tropospheric delay estimates produced every 10 minutes with GIPSY-OASIS II were converted into daily average values, and daily three-dimensional position estimates were converted into ITRF2005.
Lager Image
The differences of tropospheric delays.
- 3.1 Tropospheric Delays
In the GPS data processing, a priori values of the hydrostatic and wet delay were set at default ones, and then the remaining Zenith Delay Correction (ZDC) was estimated as a random walk process (Bar-Sever et al. 1998). The Zenith Total Delays (ZTD) is thus the sum of a priori values and ZDC estimates. In this study, we used ZDC estimates to assess the effect of mapping functions on tropospheric delays.
Figure 2 shows the differences of daily ZDC averages at each GPS station. The differences between NMF and VMF1 are marked with red dots and those between NMF and GMF with black ones. In VMF1 and GMF, the effect that the change of mapping functions has on tropospheric delays has a similar bias and seasonal variability compared with NMF. However, the differences in VMF1 have a relatively wide degree of scattering, while those in GMF have a relatively narrow one. Also one can see that the differences of tropospheric delays have maximum values in either February or August.
The time series of the differences shown in Figure 2 were curve-fitted with sinusoidal functions and the resulting amplitude is shown in Table 2 along with the maximum differences. The maximum values of differences in Table 2 quantitatively present the maximum changes that VMF1 and GMF have against NMF. The differences in tropospheric delays caused by different mapping functions display the largest difference and amplitude at BJFS and the least at GMSD. However, as mentioned above, the maximum differences in Table 2 are the values smoothed with curve-fitting. Thus even if the maximum difference was 2.1 mm in NMF-VMF1 at BJFS, the largest difference might be up to 5 mm according to Figure 2 . In addition, as the differences of tropospheric delay estimates were analyzed as daily average values in this experiment, the effect that the change of mapping functions have on tropospheric delays would be even larger if tropospheric delays are estimated at an interval of 1 hour or 10 minutes.
The maximum difference and amplitude of the curve-fitted tropospheric delay differences.
Lager Image
The maximum difference and amplitude of the curve-fitted tropospheric delay differences.
Lager Image
The correlation between the latitude of the GPS station and the maximum difference of troposphericdelays.
Lager Image
The correlation between the latitude of the GPS station and the amplitude of the curve-fitted tropospheric delay differences.
The differences of tropospheric delay estimates are caused because the way that coefficients of NMF, VMF1 and GMF were determined is different from one another. Boehm et al. (2006b) showed that the hydrostatic coefficient c h in VMF1 has a high correlation with the latitude and the season. Based on the result of Table 2 , we analyzed the correlation of tropospheric delays with the site latitude. And the result is shown in Figures 3 and 4 . Figure 3 displays the correlation between maximum difference and latitudes of the GPS stations, and Figure 4 between amplitude and latitude of the GPS site.
Lager Image
The time series of the vertical differences by changing mapping functions.
The correlation coefficients of the maximum difference of tropospheric delays and site latitude were larger than 0.9. As illustrated in Figure 3 , the differences between NMF and GMF are the least in the area of 30 o latitude. The periodicity for the differences of VMF1 and GMF relative to NMF reflects the influence of the season with their magnitude represented in amplitude. As seen in Figure 4 , the correlation coefficients between amplitude and the latitude are slightly larger 0.7, which are less than those between the maximum difference and the latitude. Besides, the higher latitude goes, the larger the amplitude gets. Strangely, two GPS stations located at around 30 o N and 35 o N have less amplitude than those at ∼25 o N and they are GMSD and USUD, respectively, located in Japan. Except for those two stations, it is certain that amplitude gets larger as latitude goes higher.
- 3.2 Vertical Positions
In space geodetic techniques, the precision of the vertical positions has a high correlation with tropospheric delay estimates. MacMillan & Ma (1994) showed that the height error is approximately one-fifth of the tropospheric delay errors at the elevation angle of 5 o , and Park et al. (2007) found that the error of vertical crustal deformation and the variation of ZWD had a scale factor of 3.72, depending on the correction of ocean tide loading. To analyze the variations of height estimates by different mapping functions, we used the same GPS data and processing approach as in section 3.1. Figure 5 presents the differences of vertical positions. The differences between NMF and VMF1 are marked with red dots, and those between NMF and GMF with black dots. The height differences by changing mapping functions have the same kinds of annual signals as tropospheric delays, and their magnitudes are almost twice as those of tropospheric delays in the opposite direction.
Table 3 lists the maximum difference and amplitude after curve-fitting the time series of the height differences. As shown in the table, the maximum differences and amplitudes by changing mapping functions have the greatest value at BJFS, and the least at GMSD. This is the same result as in tropospheric delays. Two largest height differences were found to be -4.2 mm at BJFS and -4.1 mm at SUWN. However, as illustrated in Figure 5 , the height differences can reach up to ≫9 mm in the negative direction at BJFS and SUWN.
The maximum difference and amplitude of the curve-fitted height differences.
Lager Image
The maximum difference and amplitude of the curve-fitted height differences.
Lager Image
The correlation between the latitude of the GPS station and the maximum difference of verticalpositions.
Lager Image
The correlation between the latitude of the GPS station and the amplitude of the curve-fitted height differences.
Figures 6 and 7 show the correlation between the height differences and the latitude of the GPS stations. Figure 6 presents the correlation between the maximum difference and latitude, and Figure 7 between amplitude and latitude. The figures show that the maximum differences and amplitudes of the height time series have high correlations with the latitude of the GPS stations. The height differences have the least value around 30 o N just as in the tropospheric delay estimates in Section 3.1.
It was anticipated that the vertical coordinate variations due to different mapping functions would be accompanied by the changes in the estimated vertical velocity and their RMS (Root Mean Square) values. Thus, we analyzed the variation of the vertical velocity and RMS by changing mapping functions and assessed their effects. The vertical time-series produced by using NMF, VMF1 and GMF were fitted as a linear regression to estimate the linear velocity and the result is listed in Table 4 . In the table, the velocity differences of VMF1 and GMF with respect to the one determined with NMF are also shown. The RMS values of each mapping function and the improvement ratios relative to the NMF are shown in Table 5 .
The vertical velocity change due to different mapping functions.
Lager Image
The vertical velocity change due to different mapping functions.
The variation of vertical velocity and its precision due to mapping functions.
Lager Image
The variation of vertical velocity and its precision due to mapping functions.
As listed in Table 4 , vertical velocities showed maximum differences of -0.2 mm/yr at AIRA and USUD, which means that the change of mapping function had little effect on the vertical velocity. From Table 5 , we found that RMS or the precision got better at all the GPS stations excluding AIRA and TSKB. For those stations where RMS got lower, VMF1 had more improvement than GMF. One interesting thing to note is that the three sites showing more than 4 % improvement of RMS are located in South Korea and China. And the improvement is less than 2 % or negative for those stations located in Japan. Consequently, it could be hypothesized that the accuracy of the latest mapping functions has a regional or spatial correlation. However, the number of sites used in this study is not large enough to validate this hypothesis. Also further investigation is due for the AIRA and TSKB stations where the vertical precision got worse.
4. Conclusion
To analyze the effect of mapping functions, nine GPS stations located in 20 o - 40 o N and 110 o - 145 o E were selected and their data were processed with the VMF1, GMF, and NMF mapping functions. As a result, the differences of tropospheric delays and vertical coordinates caused by different mapping functions displayed annual signals and showed maximum values in either February or August. Additionally, the maximum difference was observed for the case of VMF1, compared with NMF as reference. The largest difference in the daily average of tropospheric delays was ≫5 mm, and that of vertical positions was ∼9 mm. The differences in tropospheric delay and height estimates had high correlations with the latitude of GPS station, and they had the minimum correlation around 30 o N. As well, the correlation between the latitude of GPS stations and the amplitude of curve-fitted differences had a smaller correlation coefficient than that between the latitude of GPS stations and the maximum difference. But the higher latitude goes, the larger amplitude gets excluding two sites. To analyze the effect of the vertical time series induced by different mapping functions, we estimated three different sets of linear velocities at each site. The result shows that the vertical velocity change had a maximum value of -0.2 mm/yr. Therefore, we concluded the change of mapping functions had little influence on the vertical velocity. However, the precision of the vertical time-series was improved almost at all stations.
Acknowledgements
This work was supported by KRCF Grant (1345086367) of Ministry of Education,Science and Technology in 2009.
References
Bar-Sever Y. E , Kroger P. A , Borjesson J. A 1998 JGR 103 5019 -
Boehm J , Schuh H , Schwegmann W , Thorandt V 2003 Verlag des Bundesamts fur Kartographie und Geodasie Frankfurt am Main 16th Working Meeting on European VLBI for Geodesy and Astrometry 131 -
Boehm J , Niell A. E , Tregoning P , Schuh H 2006a GRL 33 L07304 -
Boehm J , Werl B , Schuh H 2006b JGR 111 B02406 -
Davis J. L , Herring T. A , Shapiro I. I , Rogers A. E , Elgered G 1985 Radio Science 20 1593 -
Ha J , Park K. D , Park P. H , Park J. U 2005 Journal of the Korean Meteorological Society 41 337 -
Herring T. A , DeMunk J. C , Spoelstra T. A 1982 Netherlands Geod commission Series Delft Proceedings of the Symposium on Refraction of Transatmospheric Signals in Geodesy 157 -
Ifadis I 1986 Technical Report No 38L Chalmers University of Technology Gothenburg
Lanyi G 1984 TDA Progress Report 42-78 Vol April-June JPL Publication Pasadena 152 -
MacMillan D. S , Ma C 1994 JGR 99 637 -
Marini J. W 1972 Radio Science 7 223 -
Munekane H , Kuroishi Y , Hatanaka Y , Yarai H 2008 Geophys J Int 175 831 -
Niell A. E 1996 JGR 101 3227 -
Park K. D , Won J , Kim H. K , Lim K. C 2007 JA&SS 24 249 -
Saastamoinen J , Henriksen S. W , Mancini A , Chovitz B. H 1972 Geophysical monograph vol15 Atmospheric correction for troposphere and stratosphere in radio ranging if satellites American Geophysical Union Washington DC 247 -
Webb F. H , Zumberge J. F 1993 in An introduction to the GIPSY/OASIS II JPLPublication Pasadena D-11088 -