A New Method of GPS Water Vapor Tomography for Maximizing the Use of Signal Rays

The spatio-temporal distribution of atmospheric water vapor information can be obtained by global positioning system (GPS) water vapor tomography. GPS signal rays pass through the tomographic area from different boundaries because the scope of the research region (latitude, longitude, and altitude) is designated in the process of tomographic modeling, the influence of the geographic distribution of receivers, and the geometric location of satellite constellations. Traditionally, only signal rays penetrating the entire tomographic area are considered in the computation of water vapor information, whereas those passing through the sides are neglected. Therefore, the accuracy of the tomographic result, especially at the bottom of the area, does not reach its full potential. To solve this problem, this paper proposes a new method that simultaneously considers the discretized tomographic voxels and the troposphere outside the research area as unknown parameters. This method can effectively improve the utilization of existing GPS observations and increase the number of voxels crossed by satellite signals, especially by increasing the proportion of voxels penetrated. A tomographic experiment is implemented using GPS data from the Hong Kong Satellite Positioning Reference Station Network. Compared to the traditional method, the proposed method increases the number of voxels crossed by signal rays and the utilization of the observed data by 15.14% and 19.68% on average, respectively. Numerical results, including comparisons of slant water vapor (SWV), precipitable water vapor (PWV), and water vapor density profile, show that the proposed method is better than traditional methods. In comparison to the water vapor density profile, the root-mean-square error (RMS), mean absolute error (MAE), standard deviation (SD), and bias of the proposed method are 1.39, 1.07, 1.30, and −0.21 gm−3, respectively. For the SWV and PWV comparison, the RMS/MAE of the proposed method are 10.46/8.17 mm and 4.00/3.39 mm, respectively.