Evaluation of Wet Mapping Functions Used in Modeling Tropospheric Propagation Delay Effect on GPS Measurements

Methods and Materials: The author compared the following mapping techniques against ray-traced delays computed for radiosonde profiles under the assumption of spherical symmetry: Saastamoinen, Hopfield, Black, Chao, Ifadis, Herring, Niell, Moffett, Black and Eisner and UNBabc mapping functions. Radiosonde data were computed from radiosonde stations at the Egyptian stations; in the south of Egypt, near the Mediterranean Sea, and near the Red Sea over a period of 5 years (2000-2005), most of the stations launched radiosonde twice daily, every day of the year. Moreover, data is received from the Egyptian Meteorology Authority.


INTRODUCTION
The tropospheric refraction, among other error sources like the ionospheric refraction, multipath and signal obstructions, is one of the main factors limiting accuracy of the Global Navigation Satellite Systems (GNSS) measurements.Nowadays, one can observe the increasing demand of reliable, accurate and fast results from precise GNSS positioning.GNSS technique is commonly used in geodesy, land surveying, civil engineering, deformation monitoring and structural monitoring where precise, centimeter-level accuracy height of sites is often compulsory.Unfortunately, the tropospheric refraction is strongly correlated with the vertical coordinate component of a site and affects this component rather than the horizontal ones [1].
Satellite signals from GNSS satellites suffer from the influence of the Earth's atmosphere, which deteriorates the accuracy of the obtained positions.The major part of this phenomenon is caused by the ionized part of the upper atmosphere and is termed the ionospheric refraction.Fortunately, an application of dual (or multi) frequency signals allows GNSS users to mitigate this impact [2].Additionally, electromagnetic waves are refracted in the neutral, lower part of the Earth's atmosphere.Influence of the neutral atmosphere is commonly referred to as the tropospheric refraction [3].In detail, the GNSS signal slows down and its path is getting curved.By contrast to the ionospheric delay, this influence cannot be eliminated using the dual frequency relationship.It is because of non-dispersive nature of the neutral atmosphere for microwave frequencies.Changes in speed and path cause delay in travel of signal in comparison to travel in vacuum.The delay due to the tropospheric refraction may be expressed as [4]: where, D trop is the tropospheric delay and n (s) is the refractive index as a function of path length.It is common to divide the troposphere into hydrostatic and wet part due to the water vapor content.The tropospheric zenith total delay may be expressed as a sum of hydrostatic and wet zenith delays: where, ZTD is the zenith total delay, ZHD the zenith hydrostatic delay, and ZWD the zenith wet delay.
Significant errors may be introduced when ZTDs are mapped into slant tropospheric delays with mapping functions.Converting (mapping) zenith delays into slant delays with the use of proper and accurate mapping function is of crucial interest because for low elevation angles (<15 degrees) the slant delay is ten times larger than for zenith direction [5].
The slant total delays toward the satellite are computed using the following formula [4]: where, e is the elevation angle towards the satellite, STD (e) the slant total delay, ZHD the zenith hydrostatic delay, md (e) the mapping function for the hydrostatic delay, ZWD the zenith wet delay, and mw (e) is the mapping function for the non-hydrostatic delay.
The hydrostatic delay is caused by atmospheric gases that are in hydrostatic equilibrium.This is usually the case for dry gases and part of the water vapor.This part of the delay can very well be modeled based on the surface air pressure.The wet delay, caused by water vapor that is not in hydrostatic equilibrium, is the source of our problems.We can derive models for this delay based on the partial pressure of water vapor or relative humidity at the surface, but these models have low accuracy and need empirical constants that may vary widely with location and time of year.
Most works concerning the assessment and development of tropospheric delay models and mapping functions has been reported by Forgues [6], Mendes [7], Langely and Guo [8] and others scientists.Forgues [6] simulated the impact of 15 mapping functions on GPS positioning as a function of a large number of factors such as the elevation angle, the site location, the duration of the observation session, and the estimation of tropospheric parameters.Using the Mapping Temperature Test (Herring [9],) function as reference, she concluded that the functions by Davis et al [10], Lanyi [11], Ifadis [12], and Niell [13] performed the best for short and long baselines.Mendes [7] presents comprehensive assessment of the most significant models and mapping functions developed in the last three decades against ray tracing data from 50 stations.He optimized the performance of mapping function based on some developed models.He recommended the Saastamoinen model to be used in predicting the zenith dry and wet delay.If meteorological measurements aren't available, He recommended Niell mapping function [13] to be used in modeling the elevation dependence of the zenith delay.If that information is available, and for elevation angle above 6°, Ifadis [12], Lanyi [11], and Herring [9] mapping functions will likely lead to identical results.For elevation angles below 6° the Lanyi mapping function [11] is no longer recommended.
Langely and Guo [8] reported on an investigation to determine the error of several mapping functions at elevation angles as low as 2° and presented a new model, UNBabc, which has better low-elevation angles performance and requires only slightly more computation time.They concluded that UNBabc is the best candidate for a Global Navigation Satellite Systems (GNSS) receiver built-in mapping function as it takes much less time to compute than Niell mapping function [13] does and has much better accuracy than Black and Eisner (B&E) [14] mapping function.
But most standard tropospheric models and mapping functions were experimentally derived using available radiosonde data, which were mostly observed on the European and North American continents.There is therefore a strong motivation to evaluate the accuracy of the current strategies proposed for wet tropospheric delay modeling especially for atmospheric conditions of Egypt.
The main research contribution from this work is the identification, classification, and comparison of alternative wet mapping function for the ray-path and the atmospheric structure employed in ray-tracing.It is a two-part contribution, parts that we now discuss.Firstly, a numerical integration based model is devoted to determine wet tropospheric delay at different zenith angles, which represent the benchmark values for the assessment of wet tropospheric delay models and mapping functions used for tropospheric modeling in GPS data in Egypt.Secondly, an assessment of ten wet mapping functions and functional formulations for describing both the elevation angle and azimuth-dependence of the wet tropospheric delay is performed.These mapping functions are Saastamoinen, Hopfield, Black, Chao, Ifadis, Herring, Niell, Moffett, Black and Eisner and UNBabc mapping functions.The basic motivation for this assessment of ten mapping functions through comparison with ray tracing through radiosonde data is determining the better wet mapping functions in the analysis of space geodetic data.

MAPPING FUNCTIONS
The mapping function, m (ε), is defined as the ratio of the electrical path length through the atmosphere at geometric elevation ε, to the electrical path length in the zenith direction.A mapping function is used to map the zenith delay to estimate the slant tropospheric delay.Several mapping functions have been developed in the past 20 years.The simplest mapping function is given by 1/sin(ε) [15], the cosecant of the elevation angle.In this derivation, it is assumed that spherical constant height surfaces could be approximated as planar surfaces.This is an accurate approximation only for high elevation angles and with a small degree of bending.More complex mapping functions have been developed, and different mapping functions may be used for the hydrostatic versus wet delays.Brief descriptions of the main features of various wet mapping functions are given in the following sub-sections.

Saastamoinen Model
Saastamoinen developed a total delay model for the troposphere from which a "neutral" mapping function can be derived.It is given here because it can be applied for the wet delay as well as for the dry delay and has been popular among broad parts of the GPS user community in the past.Several stages of refinement exist for Saastamoinen's approach [16].

Hopfield Model
Hopfield used real data covering the whole Earth.He has empirically found a representation of the wet refractivity as a function of the height h above the surface by, where, the mean value h w = 11000 m is used.
The result equations for wet components in (m) are [3]; where, T 0 and e 0 are measuring at the observation location and ε is the elevation angle.

Black Model
The wet mapping functions derived by Black are based on the quartic profiles developed by Hopfield and use the equivalent heights proposed by Hopfield [3].This mapping function was recommended for elevation angles above 5å nd the slant wet delay is [17]: where: e 0 : partial water vapor pressure at site in [mb], H w : upper boundary height for the wet delay/height of the tropopause, r: radial distance from earth center to site, The H w in eq. ( 7) can be approximated with help of surface temperature t as: and the scale factor L c in eq, ( 7) is: L c = 0.167 -(0.076 + 0.00015 t 0 ) .exp( -0.3 e). (9)

Chao Wet Mapping Function
Marini gave the elevation angle dependence of the atmospheric delay in the form a continued fraction, in terms of the sine of elevation angle ε [18]: where, the coefficients a, b, c… are constants or linear functions which depend on surface pressure, temperature, lapse rates, and height.
In the Chao mapping functions, the continued fraction in eq. ( 10) is truncated to second order terms and the second order sin ε is replaced by tan ε, and the coefficients a and b are determined from empirical data.The wet mapping functions are expressed as follows [19]:

Ifadis Wet Mapping Function
The global wet mapping function of Ifadis adopts the mapping function type [12]:

Herring Wet Mapping Function
Herring has developed both hydrostatic and wet mapping functions by fitting to radiosonde data from several North American stations ranging in geographic latitude from 27˚ N to 65˚ N for elevation angles down to 3˚.The mapping function's coefficients depend linearly on surface temperature, the cosine of the station latitude, and the height of the station above the geoid.The expression for the mapping function is given below [9]: where, a, b, and c are constants or linear functions, and can be determined by:

Niell Wet Mapping Function
The Neill mapping functions have no parameterization in terms of meteorological conditions, and they provide a better fit and give better accuracy for elevation angles down to 3˚.The form adopted for the mapping functions is the continued fraction with three coefficients (a, b, c).The wet mapping function depends only on the site latitude [13].The expressions for the wet mapping functions are given below: where: ε is the elevation angle, and for the wet mapping function, only an interpolation in latitude for each parameter is needed.The average values of a wet , b wet and c wet and their amplitudes are given in (

Moffett Wet Mapping Function
Simplified approximations of the Hopfield mapping functions are presented in [20].These simplified mapping functions have been used extensively, as they depend on the elevation angle only as in eq. ( 15): (15)

Black and Eisner MF
The Black and Eisner mapping function for the total delay is a further modification of Black's, and is expressed as a simple geometrical model (16), which depends on the elevation angle only, with one fitted parameter [14]: This mapping function is claimed to be valid for elevation angles greater than 7˚.

UNBabc MF
Jiming Gue and Richard B. Langley were produced a mapping Function called (UNBabc) had a 3-term continued fraction form as in eq. ( 10) [8].From a series of analyses, they concluded that parameter (a) is sensitive to the orthometric height (H) and latitude (φ) of the station where as parameters b and c could be represented by constants.The least-squares estimated parameters for the wet mapping function: The mapping functions assessed in this study and the corresponding codes and references used in the tables of results are listed in (Table 2).2003), ε, φ, H 2ẘ here: φ is the station latitude, H is the station height, r is the effective radius of the station, T is the temperature at site in K˚, t is the temperature at site in ˚C, e is the partial water vapor pressure in (mb), z is the zenith angle and ε is the elevation angle.** Means that: not specified or valid for any elevation angle.
Saastamoinen model, Hopfield model and black model give the wet slant tropospheric delay directly but in other mapping functions used for the wet zenith tropospheric delay, the Saastamoinen model eq.( 4) [21].1. Atmospheric modeling, it presents height profiles for pressure, temperature, and humidity.

𝑚(𝜀)
2. Delay modeling, it aims to determine wet tropospheric delay at different zenith angles for three stations Aswan, Helwan, and Mersa-Matrouh in different times of year in Egypt (Numerical Integration Model).
3. Assessment of prediction models or mapping functions, this stage aims to determine best models and mapping functions which present high accuracy in wet tropospheric delay prediction for Egypt.

Ray Trace Model
In order to evaluate mapping functions, we ray-traced the refractivity profiles computed from the pressure, temperature and relative humidity profiles from balloon flights launched from 5 radiosonde stations at the Egyptian stations in the south of Egypt, near the Mediterranean Sea, and near the Red Sea over a period of 5 years (2000-2005).Most of the stations launched radiosonde twice daily, every day of the year.
When ray-trace methods are used to calculate range and elevation-angle errors, integrals for the bending, electromagnetic ray path, and the straight line path or geometric path (or the sums to which these integrals are reduced by layering) must be calculated for each of the many data points of a satellite pass.The integrals, moreover, must be calculated with great precision since the range error is calculated by taking the difference between the electromagnetic path and the geometric path, two nearly equal numbers.
In the method given by Thayer [22], the number of calculations required over a satellite pass is considerably reduced since the coefficients used to evaluate the integrals are independent of the elevation angle of the ray.The elevation angle at each layer must be calculated for each new ray by using Snell's law, however, and the need for great precision remains.

Assessment of wet Tropospheric Correction Models of GPS Measurements
The Open Atmospheric Science Journal, 2017, Volume 11 9 Analysis data received from ray-traced model showed that, wet tropospheric delay in Egypt varies from 66.84 mm to 223.34 mm at zenith and slant wet delay increase slowly with high elevation angles from 30˚ -90˚.At elevation angle 60˚, wet delay error varies from min value 77.18 mm at Aswan station in April to max value 276.34 mm at Marsa Matrouh station in July month.For elevation angle 30˚, it changes from min value 133.69 at Aswan station in April to max value 478.29 mm at Marsa Matrouh station in July.
For the elevation angles lower than 30˚ the slant wet delay errors are increased suddenly.For elevation angles 15, wet delay error is 258.31 mm at Aswan station in April (min value) and 921.26 mm in July at Helwan station (max value).Wet delay errors at elevation angle 10 change from 385.09 mm to 1366.49mm. for elevation angle 7; it arrives 1688.86 mm as min value but arrives 2323.43 mm as max value at elevation angle 5. Analysis of data in Figs.(2)(3)(4) resulted that wet tropospheric delay in Egypt increases from south of Egypt to Cairo region and increases near the Mediterranean Sea.Wet delay in Egypt in summer months is more than other seasons in middle and north of Egypt but in Upper Egypt max wet delay values are measured in autumn months.

RESULTS AND DISCUSSION
The results of the assessment are summarized in Tables (3-7) and the mean biases and rms errors of these mapping functions are showed in (Figs.5-9).Table 3. Mean and rms of the wet tropospheric delay errors of the differences between the delays computed using the mapping functions and the ray trace results at 75˚ and 85˚ elevation angles.

Fig. (5).
Mean difference of the wet tropospheric delay errors between the delays computed using the mapping functions and the ray trace results for elevation angles from 30˚ to 90˚.

Fig. (6).
Mean difference of the wet tropospheric delay errors between the delays computed using the mapping functions and the ray trace results for elevation angles from 5˚ to 30˚.  7).Root mean square wet tropospheric delay errors about the mean of the differences between the delays computed using the mapping functions and the ray trace results for 30˚ and 60˚ elevation angles.

Fig. (8).
Root mean square wet tropospheric delay errors about the mean of the differences between the delays computed using the mapping functions and the ray trace results for 10˚ and 15˚ elevation angles.9).Root mean square wet tropospheric delay errors about the mean of the differences between the delays computed using the mapping functions and the ray trace results for 5˚, 7˚ and 9˚ elevation angles.
For elevation angles above 30˚, virtually all mapping functions except Hopfield MF and Black wet delay model yield errors of less than 23 mm with rms errors no more than 7 mm.for elevation angle 30˚, mean error difference between Hopfield MF and ray-trace model arrives 35.93 mm with rms 10.46 mm but for Black wet delay model the mean difference is -38.15 mm and rms is 14.21 mm, means that the Hopfield model and Black wet delay model is not citable for atmospheric conditions of Egypt.From analysis of the figures, SA, CH, IF, HE, NI, MO, B & E and UNBabc mapping functions provide sub-centimeter accuracy for angles above 15˚ with minimum difference between them no more than 4 mm in mean difference and 2 mm in rms values.For elevation angles below 10˚, only a few of the functions are found to adequately meet the requirements currently imposed by space geodetic techniques.
The SA, CH and MO mapping functions are quite accurate at elevation angles above 10˚.Since, at elevation angle 10˚, the mean difference for SA mapping functions arrives 57.64 mm with rms 16.83 mm, and 64.93 mm with rms 18.93 mm for CH mapping functions and 65.38 mm with rms 19.06 mm for MO mapping functions.For low elevation angles less than 10˚, The Saastamoinen mapping functions perform well.However, the mean difference by SA mapping function compared to the ray tracing for the data set used at elevation angle 7˚ is 76.52 mm with rms 22.97 mm and 76.52 mm with rms 22.97 mm for elevation angle 5˚, comparing with 92.51 mm with rms 26.94 mm for CH mapping functions at 7˚ elevation angle and 128.85 mm with rms 37.44 mm at elevation angle 5˚.
Both the Hopfield and Black mapping functions differences with respect to ray tracing indicate some seasonal and/or latitudinal dependence.This might be caused by the use of nominal values for the tropopause height and temperature lapse rate.However, it seems likely that such nominal values will be used in space geodetic software since values of such parameters for specific sites and atmospheric conditions are generally not exactly known.
Table 8 summarizes the best five mapping functions based on the mean difference, rms and relative error as compared to the delay computed using ray trace model for Egypt.For this data set, the Saastamoinen, Chao, Moffett, Ifadis and Herring mapping functions provided the same level of accuracy for elevation angle more than 10˚.For elevation angle less than 10˚, as can be seen in Table 8b, and 8c, the Saastamoinen, Chao and Moffett functions consecutively place first, second and third, respectively.The results states that if information is available on the vertical temperature and water vapor pressure distribution in the atmosphere, the Saastamoinen mapping function should be used.Otherwise, one of the mapping functions derived by Chao and Moffet should be used.

Fig. ( 1 )
Fig. (1) represents the flow chart followed in this Research.It contains several stages as follows:
Fig.(7).Root mean square wet tropospheric delay errors about the mean of the differences between the delays computed using the mapping functions and the ray trace results for 30˚ and 60˚ elevation angles.
Fig.(9).Root mean square wet tropospheric delay errors about the mean of the differences between the delays computed using the mapping functions and the ray trace results for 5˚, 7˚ and 9˚ elevation angles.