Abstract

Global Navigation Satellite System (GNSS) troposphere tomography has become one of the most cost-effective means to obtain three-dimensional (3-d) image of the tropospheric water vapour field. Traditional methods divide the tomography area into a number of 3-d voxels and assume that the water vapour density at any voxel is a constant during the given period. However, such behaviour breaks the spatial continuity of water vapour density in a horizontal direction and the number of unknown parameters needing to be estimated is very large. This is the focus of the paper, which tries to reconstruct the water vapor field using the tomographic technique without imposing empirical horizontal and vertical constraints. The proposed approach introduces the layered functional model in each layer vertically and only an a priori constraint is imposed for the water vapor information at the location of the radiosonde station. The elevation angle mask of 30° is determined according to the distribution of intersections between the satellite rays and different layers, which avoids the impact of ray bending and the error in slant water vapor (SWV) at low elevation angles on the tomographic result. Additionally, an optimal weighting strategy is applied to the established tomographic model to obtain a reasonable result. The tomographic experiment is performed using Global Positioning System (GPS) data of 12 receivers derived from the Satellite Positioning Reference Station Network (SatRef) in Hong Kong. The quality of the established tomographic model is validated under different weather conditions and compared with the conventional tomography method using 31-day data, respectively. The numerical result shows that the proposed method is applicable and superior to the traditional one. Comparisons of integrated water vapour (IWV) of the proposed method with that derived from radiosonde and European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim data show that the root mean square (RMS)/Bias of their differences are 3.2/−0.8 mm and 3.3/−1.7 mm, respectively, while the values of traditional method are 5.1/−3.9 mm and 6.3/−5.9 mm, respectively. Furthermore, the water vapour density profiles are also compared with radiosonde and ECMWF data, and the values of RMS/Bias error for the proposed method are 0.88/0.06 g/m3 and 0.92/−0.08 g/m3, respectively, while the values of the traditional method are 1.33/0.38 g/m3 and 1.59/0.40 g/m3, respectively.

Highlights

  • Global Navigation Satellite System (GNSS) tropospheric tomography has become a new tool for studying the three-dimensional (3-d) distribution of atmospheric water vapour information. 3-d water vapour information can reflect the changes in water vapour images at different altitudes, and show the vertical motion of water vapour to provide more information about tropospheric environmental variations [1,2,3]

  • Some studies show that the integrated water vapour (IWV) map derived from interferometric synthetic aperture radar (InSAR) can be used for improving the accuracy of the troposphere tomographic results [15,16,17]. [11] proposed a new parameterised method that can yield a superior result compared to the constant parameterization method

  • The calculated slant wet delay (SWD) value is projected by the estimated parameters that are derived from the Global Positioning System (GPS) observations based on the GAMIT/GLOBK (v10.5) software [23] and the formula is expressed by [2] as: SWD = mw(Ele) · ZWD + mw(Ele) · cot(Ele) · (GNw · cos(Azi) + GWw · sin(Azi))

Read more

Summary

Introduction

Global Navigation Satellite System (GNSS) tropospheric tomography has become a new tool for studying the three-dimensional (3-d) distribution of atmospheric water vapour information. 3-d water vapour information can reflect the changes in water vapour images at different altitudes, and show the vertical motion of water vapour to provide more information about tropospheric environmental variations [1,2,3]. Subject to the constellation of GNSS satellites and the geometric distribution of ground-based receivers in a regional network, some voxels are not crossed by any satellite ray when following the method proposed by [2], which leads to an inversion problem of the tomographic observation equation. To solve this problem, some constraints (horizontal, vertical, and boundary constraints) were imposed on the tomographic model in previous studies [2,6,7,8,9,10,11,12,13].

Principle of the Proposed Parameterised Approach for Tropospheric Tomography
Horizontal Parameterised Approach for Tropospheric Tomography
Prior Constraint
The Weights Strategy for Proposed Tomography Approach
Selection of the Weights for the Same Type of Equation
Determination of the Weights for Different Types of Equations
Processing Strategy
Nu8mber
External Accuracy Validation
Comparison with the Traditional Tomographic Method
Uniformity of the Height System
Water Vapour Profile Comparison
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call