**Journal of Geosciences and Geomatics**

## Determination of the best-fit Tropospheric Delay Model on the Nigerian Permanent GNSS Network

**Dodo Joseph Danasabe**^{1,}, **Ojigi Lazarus Mustapha**^{2}, **Tsebeje Samuel Yabayanze**^{1}

^{1}Space Geodesy and Systems Division, Centre for Geodesy and Geodynamics, Toro, Nigeria

^{2}Mission Planning, IT and Data Management, National Space Research and Development Agency, Abuja, Nigeria

### Abstract

The Federal Government of Nigeria through the Office of the Surveyor General of the Federation (OSGoF) set up surveying infrastructure throughout the country known as the NIGerian Reference GNSS NETwork (NIGNET). The NIGNET is a network of Global Navigation Satellite System (GNSS) Continuously Operating Reference Stations (CORS) set up at different locations in Nigeria for surveying and mapping. They are satellite tracking stations operating 24 hours a day providing positional solutions. As signals from the satellite pass through the different layers of the atmosphere (ionosphere and troposphere), they are refracted thus, causing delay on the arrival of the signal, which in-turns affect positioning in the horizontal and height component. The most dominant spatially correlated bias is the tropospheric effect on the GNSS satellite signals. Several global tropospheric delay models are in use by different countries to mitigate the biases cause by the troposphere. This study therefore aim to determine the best-fit tropospheric delay model for the NIGerian GNSS Reference NET work (NIGNET) using data collected from the NIGNET stations across Nigeria. Three different global tropospheric models, namely; the Saastamoinen model, Hopfield model and Niell models were used, and results compared. Four processing strategies were adopted. The first strategy was without the application of any of the models, while in the second, third and fourth strategies, the GNSS data were processed with the application of each of the models. The results indicate that, the Niell model has the lowest mean zenith tropospheric delay (ZTD) of 2.330m with root mean square error (rmse) of 0.45m, while the Hopfield and Saastamoinen models have ZTD of 2.386m and 2.398m with rmse value of 0.60m and 0.71m respectively. On the overall, the Niell model has better performance in the network. This suggests that, the application of Niell model in the processing of all GNSS data will give a more reliable result in the position domain as well as the height component. The results are very useful to surveyors and geodesist engaged in surveying and mapping, and spatial positioning of infrastructures. It will enhance the effectiveness and reliability of the tropospheric delay resolution process for regional Global Positioning System (GPS) network users.

**Keywords:** ** **tropospheric delay, Hopfield model, Saastamoinen model, Niell model, NIGNET

Received July 14, 2015; Revised August 05, 2015; Accepted August 09, 2015

**Copyright**© 2015 Science and Education Publishing. All Rights Reserved.

### Cite this article:

- Dodo Joseph Danasabe, Ojigi Lazarus Mustapha, Tsebeje Samuel Yabayanze. Determination of the best-fit Tropospheric Delay Model on the Nigerian Permanent GNSS Network.
*Journal of Geosciences and Geomatics*. Vol. 3, No. 4, 2015, pp 88-95. http://pubs.sciepub.com/jgg/3/4/1

- Danasabe, Dodo Joseph, Ojigi Lazarus Mustapha, and Tsebeje Samuel Yabayanze. "Determination of the best-fit Tropospheric Delay Model on the Nigerian Permanent GNSS Network."
*Journal of Geosciences and Geomatics*3.4 (2015): 88-95.

- Danasabe, D. J. , Mustapha, O. L. , & Yabayanze, T. S. (2015). Determination of the best-fit Tropospheric Delay Model on the Nigerian Permanent GNSS Network.
*Journal of Geosciences and Geomatics*,*3*(4), 88-95.

- Danasabe, Dodo Joseph, Ojigi Lazarus Mustapha, and Tsebeje Samuel Yabayanze. "Determination of the best-fit Tropospheric Delay Model on the Nigerian Permanent GNSS Network."
*Journal of Geosciences and Geomatics*3, no. 4 (2015): 88-95.

Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |

### At a glance: Figures

### 1. Introduction

One of the fundamental issues in Network GPS is the ability to mitigate all potential errors and biases in the system. The term bias here refers to a physical phenomenon whereas the term error refers to the quantity remaining after the bias has been mitigated ^{[2]}. Error sources are the satellite-related errors, (satellite coordinate errors, satellite clock offsets and satellite ephemeris errors), the atmospheric-related errors (tropospheric and ionospheric errors) and the station-related errors (receiver clock offsets, antenna phase centre variations, multipath, solid earth tides and ocean tide loading). The carrier phase measurements are compromised by these errors; as such most of the errors except for troposphere, receiver clock and ionospheric delay can be mitigated to some extent through modelling ^{[12]}. The ionospheric delay, which is a function of the total electron content along the signal path, and the frequency of the propagated signal, can be eliminated because of its frequency dependency by using double-frequency ionospheric free linear combination ^{[6, 9]}.

Although the ionospheric bias can be mitigated using dual frequency receivers, tropospheric bias is currently one of the major error sources in GPS Network, which limits the full functionality of GPS Positioning. The delay of the radio signals caused by the troposphere can range from 2m at the zenith to 20m at lower elevation angles (below 10 degrees). In order to reduce the tropospheric effects, global tropospheric models derived experimentally using radiosonde data are employed today. With the establishment of the Nigerian Permanent GNSS Network as one of the latest innovation of a real-time precision positioning in meeting up with the nation’s development, security and defence; the need to investigate the impact of the different global tropospheric models became imperative.

### 2. Tropospheric Delay

The troposphere is the lower part of the atmosphere close to the earth surface. Troposphere is where the inhabitants of our planet live, it starts on the surface of the Earth and goes up to height of 9 to 18 km. The depth of it varies with the latitude; it is greatest at the equatorialregions (approximately 18 km) and minimal near the poles (about 9 km). Approximately 75-80 % of the mass of the whole atmosphere is in the troposphere including nearly all water vapour (about 99%) and dust particles ^{[22]}. The refractive index showing the ratio of the vacuum speed of light and the real speed of the signal is always bigger than 1. It means that the signals emitted by the satellites reach the receiver later, therefore a greater satellite-to-receiver range is measured. Troposphere is a non-dispersive medium for radio frequencies below 15 GHz, hence its effect is independent of GNSS frequencies. It causes delay in both GNSS carrier and code observations.

Due to the highly variable tropospheric water vapour content, it is difficult to achieve desired accuracy in this region ^{[1]}. The tropospheric delay in equation (1) is directly proportional to the refractive index, which is expressed as a function of atmospheric temperature and pressure. It is therefore expressed as ^{[6]}:

(1) |

where is the tropospheric delay, n is the refractive index and ds is the path length.

The refractivity can be divided into dry and water vapour due to the troposphere containing dry and water vapour content, hence;

(2) |

where,: refractivity of dry air: refractivity of water vapour. Expressing in terms of refractivity N; from equation (1); we obtain;

(3) |

where; : tropospheric refractivity, hence ;

(4) |

The tropospheric delay can be separated into the hydrostatic (dr) and non-hydrostatic (wet) component as shown in Figure 1.

(5) |

where; : dry tropospheric refractivity

: wet tropospheric refractivity resulting from water vapour.

The hydrostatic and wet components are caused by the dry gases (primary nitrogen and oxygen) and water vapour respectively. About 90% of the tropospheric delay is caused by the dry component, while the remaining 10% is from the wet component. The tropospheric delay is then expressed as a linear combination of the hydrostatic and wet components ^{[5]};

(6) |

or

(7) |

**Fig**

**ure**

**1.**

**Schematic diagram of the Hydrostatic and Wet components of the troposphere [9]**

Tropospheric delay is calculated in the zenith direction over the GPS station, hence the term zenith tropospheric delay, a combination of the zenith hydrostatic delay and zenith wet delay. The tropospheric delay is a function of elevation and altitude of the receiver, which depends on factors such as atmospheric temperature, pressure and relative humidity. It is not frequency-dependent as is the case with the ionosphere and cannot be eliminated through linear combination of L1 and L2 observations ^{[15]}.

### 3. Tropospheric Delay Models

Several global tropospheric models such as the Saastamoinen model, Hopfield model, Niell model etc. have been empirically developed and employed in GPS timing receivers to correct for the tropospheric delay. These models are derived using data from available radiosonde obtained from Europe and North America continents. The global atmosphere conditions, used as constants in these models, provide a broad approximation of the tropospheric conditions, but ignore the actual atmospheric conditions on a given location, i.e., do not take into account the latitudinal and seasonal variations in the atmosphere ^{[11]}. Besides, daily variation in temperature, pressure and relative humidity can lead to error in tropospheric delays obtained using the global tropospheric models especially in the height components ^{[3]}. The location of Nigeria in the equatorial and tropical region makes it susceptible to high tropospheric effect thereby having an adverse effect on the GPS signals, which, in turn, affects positioning. In order to determine the best-fit tropospheric model for processing of data collected from the Nigerian Permanent GNSS Network, the need to investigate the impact of the different global tropospheric models on the network becomes imperative. The research investigates the performance of three global tropospheric delay models, namely Refined Saastamoinen model ^{[14]}, Modified Hopfield model ^{[7]} and Niell model ^{[10]}.

**3.1. The Refined Saastamoinen Model**

The Saastamoinen model as in ^{[12]} is expressed as a function of height of the observation station and the zenith angle. This was later modified and functionally expressed as ^{[5]}:

(8) |

Where = zenith angle of satellite

= pressure (mbar)

= temperature (K)

= partial pressure of water vapour (mbar)

=tropospheric path delay in metres

and are the corrections that depends on height of the station and .

**3.2. The Modified Hopfield Model**

^{[6]} used data from different parts of the world to develop an empirical tropospheric delay model. The Hopfield model shows dry and wet refractivity components as a function of tracking station height h above the Earth's surface and is given in the following forms:

(9) |

(10) |

where is empirically determined power of the height ratio, is the polytropic thickness for the dry part (m),is the polytropic thickness for the wet part (m), is the dry tropospheric refractivity for the stations at the Earth's surface as a function of pressure (millibars) and temperature (Kelvin), is the wet tropospheric refractivity for the station at theEarth's surface as a function of water vapor, pressure, and temperature.

Inserting equations 9 and 10 into equation 6, and integrating each element with the respective integration ranges along the vertical direction (i.e. from h = 0 to h = H_{d} and from h = 0 to h = H_{w} for the dry and wet components), we then obtain tropospheric zenith delay in units of meters ^{[6]}:

(11) |

**3.3. The Niell Model**

The Niell Model is a combination of the Saastamoinen zenith path delay with Neil mapping functions ^{[10]}. The parameters (a, b, c) used in the dry and wet components of the models as expressed in equations (10) and (11) are calculated based on the interpolation of the average and seasonal variation (amplitude) values as functions of latitude and time. For the dry component:

(10) |

For the wet component:

(11) |

where; *m*_{d} and *m*_{w} is the mapping functions for dry and wet components respectively; ε is the satellite elevation angle and *H* = orthometric height *a*_{d}, *b*_{d}, *c*_{d} are the coefficients in the dry component; *a*_{w}*, b*_{w}*, c*_{w} are the coefficients in the wet component and *a*_{ht}*, b*_{ht}*, c*_{ht} are the coefficients in the height component.

### 4. The Nigerian Permanent GNSS Reference Network

Theoffice of the Surveyor General of the Federation (OSGoF) established the NIGerian Permanent GNSS NETwork (NIGNET).The goal is to implement a new reference frame for Nigeria in line with the recommendation of the United Nation Economic Commission of Africa (UNECA) through Committee on Development, Information Science and Technology (CODIST) ^{[8]}. It is expected that, the Nigerian Permanent GNSS Reference Network as presented in Figure 2, will directly contribute to the Africa Reference Frame (AFREF).

**Fig**

**ure**

**2.**

**NIGNET CORS [Authors]**

About Eleven (11) NIGNET CORS stations were established at the time of conducting this research as shown in Figure 1. With the growing capabilities of GPS as a high precision positioning system for surveying and mapping, monitoring geophysical hazards, sea level change and as well as coordinating geodetic activities; there is a necessity for the NIGNET stations to be defined on the precise reference system such as International Terrestrial Reference System (ITRS) that managed the International Terrestrial Reference frame (ITRF).

### 5. Materials and Methods

The use of network of reference stations, instead of the single reference station, has become widely acceptable within the GNSS community as solution for high precision satellite positioning applications ^{[17]}. This allows modelling of the atmospheric errors such as the tropospheric propagation delays that complicate the process of ambiguity fixing, which is often considered necessary for high-precision positioning and thus, significantly reducing the errors for long baselines thereby enhancing positioning accuracy.

**5.1. Study Area**

Six (6) stations of the Nigerian Permanent Network of Continuously Operating Reference Stations (CORS) were used. The choice of these stations is due to data availability from these stations for the period of study compared to other stations of the network. Table 1 shows description of the six (6) NIGNET stations used in this research.

**5.2. Data Acquisition**

Twenty-four hours (24hrs) raw GPS data at 30-second data rate in RINEX format for the stations shown in Figure 2 and precise satellite ephemeris data for GPS week 1409 were downloaded from the International GNSS Service (IGS) for the day of the year (DoY) 01/2011 to 07/2011. The ocean tide loading data for each station was obtained from ^{[18]}. Similarly, the Earth Orientation Parameters and the Ionosphere models were downloaded from the Bernese website ^{[19]}. Summary of the parameters used are given in Table 2.

**5.3. Processing Strategy**

Four processing strategies using the Bernese GPS Scientific Software V.5.0 were employed. They include:

Strategy I: In this strategy, the processing is done without the application of the tropospheric model. The ionosphere-free double difference (IF DD) residuals and final coordinates are extracted for analysis.

Strategy II: Processing with the application of the Niell model and standard atmosphere, the IF DD residuals, final coordinates and the zenith tropospheric delay are extracted for analysis.

Strategy III: Processing with the application of the Modified Hopfield model and standard atmosphere; the IF DD residuals, final coordinates and the zenith tropospheric delay are extracted for analysis.

Strategy IV: Processing with the application of the Refined Saastaminen model and standard atmosphere; the IF DD residuals and final coordinates are extracted for analysis purpose.

The coordinates of all stations were estimated. This retains the flexibility for later changes in the realization of the reference frame. However, to check the consistency of the data used in the processing with the coordinates of the IGS core sites, a minimal constraint solution was generated for the network.

### 6. Results and Discussions

The analysis of the results was done based on the Ionospheric Free Double Difference (IF DD) residuals, the final station coordinates and the zenith tropospheric delay obtained from each of the global tropospheric delay model, in order to ascertain the best fit tropospheric delay model for the network.

**6.1. Assessment of the Tropospheric Delay Models on the basis of the Baseline Ionospheric-free Double Difference (IF DD) Residual**

One of the tools used in the assessment of tropospheric model in a GPS network is the comparison of the baseline IF DD residuals over which the tropospheric models are being assessed ^{[4]}. The Root Mean Square Error (RMSE) characterizes the performance of the models. Fifteen baselines were formed from where the RMSE were computed for all satellites.

Table 3 summarises the numerical results for all the baselines in terms of the RMS IF DD residuals. The IF DD residuals of strategy I (No model) has largeresiduals compare to strategies II, III and IV respectively. This is expected because no model is applied. The result indicates that, the three models are able to reduce the size of the residuals. However, no significant residual difference in the three models is noticed. The UNEC-RUSThaving the shortest baseline of 188.8km hasrmse value of 11.9cm when no model is applied, while the longest baseline ULAG-FUTY has RMS value of 17.0cm. This presupposes that, the tropospheric delay is distance-dependent error. This result is in agreement with ^{[21]}. The longer the baseline, the more the effects of the troposphere. The Niell model gave a better result.

Table 4 provides the percentile improvement in the RMS DD IF residuals for strategies II, III and IV respectively. From the table, the result reveals that baseline percentage improvement varies from 8% to 91%. The Niell model had the most percentage improvement in the network with an average of 46%. The modified Hopfield and Saastaminen models had 41.3% and 35.6% respectively. All stations connect to station RUST have large rmse, this could be attributed to data gabs at station RUST.

**6.2. Assessment of the Tropospheric Models in the Position Domain**

To study the tropospheric delay models in the position domain, the coordinate differences of the station in the North, East and Height (horizontal and height) components were computed and analysed. Figures 3.0, 4.0, and 5.0 shows the standard deviation of the coordinates in North, East and Height components.

The result reveals that, the Niell and Hopfield models show no significant deviations in the North and East components respectively.

However, the differences in the application of the tropospheric delay models reveals that, the Niell model shows considerable improvement with network standard deviation of 5.02m, 3.72m in the north and east component respectively, while the Hopfield Model followed closely with network standard deviation of 5.22m and 3.8m in the north and east components respectively.

**Fig**

**ure**

**3**

**.**Standard deviation in the North component

**Fig**

**ure**

**4.**Standard deviation in the East component

**Fig**

**ure**

**5.**Standard deviation in the Height component

**6.3. Mean Zenith Tropospheric Delay (ZTD) at each Station**

The tropospheric delay is calculated in the zenith direction over the GNSS station. The Zenith Tropospheric Delay (ZTD) gives insight into the tropospheric conditions above the GPS site. Table 5 show the statistics of the zenith tropospheric delay for each of the GNSS station based on the application of each tropospheric delay model.

The mean ZTD computed at each station reveals that, station RUST has the highest ZTD value of 2.503m, 2.562m and 2.592m for Niell, Hopfield and Saastamoinenmodels respectively. This is followed closely by station ULAG having 2.57m, 2.51 and 2.57m for Niell, Hopfield and Saastamoinen models respectively. This result presupposes that stations at low latitude are highly susceptible to tropospheric delay ^{[21]}.

The mean Network ZTD produced by the three tropospheric delay shows that the Niell model has the lowest network ZTD of 2.36m with mean RMS value of 0.45m, follow by the Hopfield model with network ZTD of 2.386m and mean RMS value of 0.61m. The Saastamoinen model has the highest mean network ZTD of 2.398m with mean RMS value of 0.71m.

### 7. Conclusion

This research has demonstrated the influence of different tropospheric delay models on the Nigerian Permanent GNSS Network. The result indicates that the residual tropospheric delay affects the position precision.

Increase in baseline length results in higher tropospheric effect, this is noticed on baseline ULAG-FUTY with the highest baseline length of 1060.5km. Tropospheric delay increases during the morning hours and decreases at sunset. The three models investigated i.e. the Saastamoinen, Hopfield and Niell models show no significance difference in their performance; better improvements in the position domain were achieved by the application of the Niell model compared to Hopfield and Saastamoinen models. The Niell model produced a better mitigation of the tropospheric delay, with an average percentage improvement of 46.0% while; Hopfield and Saastamoinen models have 41.3% and 35.6% percentage improvement respectively. The result also indicates that, the Niell has the lowest mean average zenith tropospheric delay (ZTD) of 2.33m with RMS of 0.45m. On the overall, the Niell model has better performance in the network in this research. This result is in agreement with ^{[20]}.

### Acknowledgement

The authors acknowledge the Office of the Surveyor General of Nigeria for providing the data used in this research.

### References

[1] | Ahn, Y. W., Lachapelle, G., Skone, S. and Sahm, S. “Analysis of GPS RTK Performance using External NOAA Tropospheric Corrections Integrated with a Multiple Reference Station Approach”, GPS Solution, 10, 171-186. 2006. | ||

In article | View Article | ||

[2] | Bingley, R, GNSS Principles and Observables, Nottingham, IESSG, University of Nottingham UK, M.Sc. Lecture Monograph. 2004. | ||

In article | |||

[3] | Dodo, J. D., and Kamarudin, M. N. “Investigation on the Impact of Tropospheric Models on baseline precision in a local GPS network: Case of the Malaysian RTKnet”, Journal of Geomatics 2 (1) 137-142. 2008. | ||

In article | |||

[4] | Don, K., Sunil, B., Langley, R. B., Dare, P. “Performance of Long-Baseline real-Time Kinematic Applications by Improving Tropospheric Delay Modelling”, ION GNSS International Technical Meeting of the satellite Division. Long Beach, California, USA. 2004. | ||

In article | |||

[5] | Guochang, X,GPS Theory, Algorithms and Application, Springer-Verlag Berlin Heidelberg, New York. 2003. | ||

In article | |||

[6] | Hofmann-Wellenhof, B., Lichtenegger, H. and Collins, J,GPS, Theory and Practice, Springer-Verlag Wien, New York. 2001. | ||

In article | |||

[7] | Hopfield, H. S, “Two Quartic Tropospheric Refractivity Profile for Correcting Satellite Data”. Journal of Geophysical Research 74(18), 4487-4499.1969. | ||

In article | |||

[8] | Jatau, B., Fernandes, R.M.S., Adebomehin, A., and Goncalves, N “NIGNET-The New Permanent GNSS Network of Nigeria”, FIG Congress 2010 Facing the Challenges – Building the Capacity Sydney, Australia, 11-16 April 2010. | ||

In article | |||

[9] | Leick, A, GPS Satellite Surveying, USA: John Willey & Sons, Inc, 2004. | ||

In article | |||

[10] | Niell, A.E, “Global Mapping for the Atmospheric Delay at Radio Wavelenghts”,Journal of Geophysical Research 111 (B2), 3227-3246.1996. | ||

In article | |||

[11] | Roberts, C., and Rizos, C, “Mitigating Differential Troposphere for GPS-based Valcano Monitoring”, 5th International Symposium on Satellite Navigation Technology and Applications. Canberra, Australia. 2001. | ||

In article | |||

[12] | Rizos, C, “Network RTK Research and Implementation - A Geodetic Perspective”. Journal of Global Positioning System, 1(2), 144-150.2002. | ||

In article | |||

[13] | Saastamoinen, J, “Atmospheric Correction for Troposphere and Stratosphere in Radio Ranging of Satellites”,Geophysical Monograph, American Geophysical Union, Washington D.C. 247-252. 1972. | ||

In article | View Article | ||

[14] | Saastamoinen, J, “Contribution to the theory of Atmosphere Refraction”,Bulletin of Geodesique (105,106,107): 279-298, 383-397, 13-34. 1973. | ||

In article | |||

[15] | Satirapod, C., and Chalermwattanachai, P, “Impact of Different Tropospheric Models on GPS Baseline Accuracy: Case Study in Thailand”, Journal of Global Positioning Systems, 4(1-2), 36-40.2005. | ||

In article | |||

[16] | Sickel, V., J, GPS for Land Surveyors, Taylor and Francis 3^{rd,} Edition Florida, United State of America, 2008. | ||

In article | |||

[17] | Vollath, U., Buecherl, A., Landau, H., Pagels, C. and Wagner, B, “Multi-Base RTK Positioning Using Virtual Reference Stations”,Institute of Navigation, ION GPS-2000, Salt Lake City, 22-19 September, 123-131.2000. | ||

In article | |||

[18] | “Ocean Tide Loading” Available: http://www.oso.cha. [AccessedSept. 20, 2009]. | ||

In article | |||

[19] | “The Earth Orientation Parameters and the Ionosphere models” Available: http://www.bernese.unnibe.ch/ [Accessed: Sept. 21, 2009]. | ||

In article | |||

[20] | Péter, B. “The evaluation of troposphere models applied in the Hungarian Active GNSS Network” The Council of European Geodetic Surveyors 2012. | ||

In article | |||

[21] | Tajul, A. M., Lim, S., Yan, T. and Rizos, C. “Mitigation of Distance-Dependent Errors for GPS Network Positioning”. International Global Navigation Satellite Systems Society IGNSS Symposium, Surfers Paradise, Australia. | ||

In article | PubMed | ||

[22] | Randy Russel, “The Toposphere, Window to the Universe”. [http://www.windows2universe.org/earth/Atmosphere/troposphere.html] , 2010. Last accessed 23/10/2012. | ||

In article | |||