快速检索        
  武汉大学学报·信息科学版  2017, Vol. 42 Issue (9): 1203-1208, 1222

文章信息

赵庆志, 姚宜斌, 罗亦泳
ZHAO Qingzhi, YAO Yibin, LUO Yiyong
附加辅助层析区域提高射线利用率的水汽反演方法
A Method to Improve the Utilization of Observation for Water Vapor Tomography by Adding Assisted Tomographic Area
武汉大学学报·信息科学版, 2017, 42(9): 1203-1208, 1222
Geomatics and Information Science of Wuhan University, 2017, 42(9): 1203-1208, 1222
http://dx.doi.org/10.13203/j.whugis20150592

文章历史

收稿日期: 2015-12-21
附加辅助层析区域提高射线利用率的水汽反演方法
赵庆志1,2, 姚宜斌1, 罗亦泳1     
1. 武汉大学测绘学院, 湖北 武汉, 430079;
2. 西安科技大学测绘科学与技术学院, 陕西 西安, 710054
摘要:在已有的全球导航卫星系统水汽反演技术中,能够利用的观测值是完整穿过整个层析区域的卫星射线。由于卫星、接收机几何位置分布以及层析区域选择的特定性,许多卫星射线是从层析区域侧面穿出,这些射线通常当作无效信息被剔除,降低了层析结果的精度。针对该缺点,提出并实现了一种附加辅助层析区域提高射线利用率的水汽反演方法,使从研究区域侧面穿出的信号信息也可以被利用。基于美国得克萨斯州(Texas)地区的连续运行参考站(Continuously Operation Reference Stations,CORS)网的实测数据,结合层析区域内的无线电探空仪数据进行实验,结果表明,该方法能够提高水汽反演结果的精度,其均方根误差的改善率为14.6%。
关键词水汽反演     CORS     射线利用率     无线电探空仪    
A Method to Improve the Utilization of Observation for Water Vapor Tomography by Adding Assisted Tomographic Area
ZHAO Qingzhi1,2, YAO Yibin1, LUO Yiyong1     
1. School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China;
2. College of Geomatics, Xi'an University of Scieace and Technology, Xi'an 710054, China
First author: ZHAO Qingzhi, PhD, specializes in GNSS data processing and GNSS Meteorology. E-mail:zhaoqingzhia@163.com
Corresponding author: YAO Yibin, PhD, professor. E-mail:regnier@ibpc.fr
Foundation support: The National Natural Science Foundation of China, Nos. 41174012, 41274022; the Jiangxi Province Natural Science Foundation of China, No. 20151BAB213031
Abstract: In existing Global Navigation Satellite System (GNSS) water vapor tomographic techniques, only signals which penetrate the entire research area are used. Due to the specific geometric distribution of satellite constellations and receivers as well as the specificity of tomography area, many signals penetrating from the side of a research area are excluded as ineffective information, which undermines the tomographic result. To solve this issue, a method is proposed for improving the utilization of signal rays by adding an assisted tomographic area, which allows the signal information crossed from the side of a research area to be utilized as well. An experiment using data from the Continuously Operation Reference Stations (CORS) network of Texas in UAS, is compared with radiosonde data of 72 249. The results show that the proposed method can enhance the accuracy of tomographic results, the improvement rate of Root Mean Square Error (RMS) was 14.6%.
Key words: water vapor tomography     CORS     signal utilization     radiosonde    

利用全球导航卫星系统(global navigation satellite system, GNSS)层析技术反演水汽已经成为获取水汽空间分布信息的重要手段之一,其高精度、全天候、高时空分辨率的优点在对暴雨等天气的监测和预报过程中得到了广泛应用[1-3]。自Bevis等[2]证明可以利用地基全球定位系统(global positioning system, GPS)获取水汽后,GNSS气象学的发展在国内外取得了许多成果[4-17]

Bevis等[2]首先利用倾斜路径湿延迟对对流层进行层析,并基于欧洲中期天气预报中心(European Centre for Medium Range Weather Forecasts, ECMWF)预报产品对层析结果进行验证;于胜杰等[6]分析了不同约束条件对水汽反演结果的影响,并证明层析结果受测站间高差和先验信息精度影响较大。叶世榕等[7]结合数值积分参数化思想,提出了一种适用于层析的三维分布数值积分方法,结果证明该方法能够明显提高层析结果垂直廓线的精度;刘志赵等[8]提出了合理选择层析区域垂直高度和整体移动层析区域以增加射线利用率的方法,但该方法对数据利用率的提高有很大的局限性。何林等[3]讨论了代数重构算法在水汽层析中的各种问题,并对松弛因子和迭代终止条件的确定做了详细研究。

上述研究都是利用完整穿过层析区域的射线信息建立观测方程。由于受接收机和卫星特定几何位置的影响,很多卫星射线不可避免会从研究区域侧面穿出,在传统层析观测方程构建过程中因无法利用而被剔除。针对该情况,本文提出了一种附加辅助层析区域提高射线利用率的方法,使得从研究区域侧面穿出的信号也可以被利用。

1 GNSS层析原理 1.1 层析观测方程

在水汽反演过程中,倾斜路径水汽含量(slant water vapor, SWV)是建立层析观测方程中最重要的输入量,定义为从卫星到接收机沿卫星信号传播路径上总的水汽含量[8],其表达式为:

$ {\rm{SWV}} = {10^{-6}}\cdot{\smallint _s}\rho \left( s \right){\rm{d}}s $ (1)

式中,ρ(s)表示水汽密度(g/m3);s表示卫星到接收机的信号传播路径长度,ds表示沿s路径上每段路径元素的长度;SWV表示信号传播路径上的水汽含量。层析技术将研究区域分成若干立体网格,并假定网格内的水汽密度为一常数,可将式(1) 写成下面离散化的形式:

$ {\rm{SW}}{{\rm{V}}^s} = \sum\limits_{i = 1}^n {(a_{_i}^{^s}\cdot{\rho _i})} $ (2)

式中,ais表示卫星信号s在第i个网格内的截距;ρi表示第i个网格内的水汽密度;n表示层析区域划分的立体网格总数。

利用层析区域内接收机所观测到的卫星信号可组成观测方程:

$ {\mathit{\boldsymbol{S}}_{l \times 1}} = {\mathit{\boldsymbol{A}}_{l \times n}}\cdot{\mathit{\boldsymbol{x}}_{n \times 1}} $ (3)

式中,S表示l个SWV观测值组成的列向量;l表示层析区域内所有测站SWV观测值的总个数;A表示网格截距组成的系数矩阵;x表示待求网格水汽密度组成的列向量。

1.2 约束方程建立

由于层析区域内接收机数目有限,导致在给定的层析时段内有射线穿过的网格数目有限,因此,式(3) 中A是维数较大的稀疏离散矩阵,不能完全求解出x中所有的待求水汽密度值[8]。为了克服该法方程秩亏问题,通常需要加入先验约束信息对网格内的水汽密度进行约束[16-18]。基于水汽在空间呈连续分布及网格距离与水汽密度相关性成正比的特点,本文采用水平平滑约束建立水平约束方程[3]。基于水汽密度随高度呈指数递减的特点,采用指数函数构建垂直约束方程[18]。此外,将垂直路径上的水汽信息也加入到层析模型中;附加上述约束后可以得到传统方法的层析模型:

$ \left( {\begin{array}{*{20}{c}} {\mathit{\boldsymbol{A}}_{_{l \times n}}^{^{{\rm{slant}}}}}\\ {\mathit{\boldsymbol{A}}_{_{m \times n}}^{^{{\rm{zenith}}}}}\\ {{\mathit{\boldsymbol{H}}_{l \times n}}}\\ {{\mathit{\boldsymbol{V}}_{l \times n}}} \end{array}} \right) \cdot {\mathit{\boldsymbol{x}}_{n \times 1}} = \left( \begin{array}{*{20}{c}} {\mathit{\boldsymbol{S}}_{l \times 1}}\\ {\mathit{\boldsymbol{P}}_{m \times 1}}\\ {{\bf{0}}_{l \times 1}}\\ {{\bf{0}}_{l \times 1}} \end{array} \right) $ (4)

式中,Al×nslant为斜方向投影函数;Am×nzenith为天顶方向投影函数;Pm×1m个测站的天顶方向水汽含量组成的列向量;HV分别表示水平约束和垂直约束的系数矩阵。基于联合代数重构算法(multiplicative algebraic reconstruction technique, MART)能够短时间内快速收敛的特点[17],对于式(4) 利用MART算法求解。

2 附加辅助层析区域反演水汽的方法 2.1 附加辅助层析区域反演水汽思想

受接收机和卫星星座特定几何位置分布的影响,层析区域内总有很多射线信号会从研究区域侧面穿出,被当作无效信息剔除。这不仅降低了观测数据的利用率,也会对反演结果造成一定的影响。

为了解决上述问题,本文提出了一种附加辅助层析区域反演水汽的方法。该方法的核心思想是在经度、纬度方向上对原有层析区域进行延伸,直至扩大后的层析区域能够确保原层析区域内所有测站一定高度截止角之内的卫星射线都能完整穿过扩大后的层析区域。对扩大后的层析区域进行水汽反演,可得到原层析区域内各网格的水汽密度初值。将此初值用于迭代并结合传统方法对层析区域的水汽进行再次反演,得到更高质量的层析结果。需要指出的是,于胜杰等[6]通过多次水汽层析实验证明,先验信息对不同GPS网中的作用差异较大。当层析区域内GPS测站间的相对高差较小时,其层析方程的系数矩阵奇异程度较为严重,因此初值信息对层析结果有很大影响;扩大层析区域范围后,扩大层析区域内的水汽参数也会增多,其水汽密度层析结果的精度和可靠性可以稍差,因为只将其结果作为初值进行下一步运算。

2.2 附加辅助层析区域具体实现过程

为了便于表达,将需要研究的层析区域(如图 1中较小的红色网格区域所示)在水平方向上经过一定的延伸形成的新区域称为辅助层析区域(如图 1中较大的绿色网格区域所示),其垂直高度保持不变。

图 1 层析区域三维分布图 Figure 1 Three-Dimensional Distribution of Tomographic Area

本文提出的水汽反演方法实现过程如下。

1) 确定要延伸的水平距离,如图 1d所示。根据层析区域的垂直高度和卫星截止高度角,以研究区域的边界为起点(假定右侧站在研究区域边界上)确定出在经度和纬度方向上要延伸的水平范围,以确保层析区域内所有测站在一定卫星高度截止角范围内的卫星射线都能够完整穿过该辅助研究区域。如图 2所示,假定进行水汽层析时利用的卫星截止高度角为α,层析区域的垂直高度为H(单位为km),在层析区域边界上有一测站P,其有一条α度角的边界射线信号S,则可以通过式(5) 计算出需要延伸的水平距离d

$ d = H/{\rm{tan}}\alpha $ (5)

2) 利用原层析区域内所有测站观测到的卫星高度角大于α的射线信息建立辅助研究区域(如图 1中较大的绿色区域)的辅助观测方程,然后附加各种约束后得到辅助层析区域的层析模型。

3) 通过SVD全科分解法[20]对第2) 步建立的辅助层析模型求解,得到辅助层析区域内每个网格的水汽密度估值。然后只选取原层析区域网格(如图 1中较小的红色区域内的网格)所在位置的水汽密度估值作为初值用于下一步解算。

4) 基于传统方法利用完整穿过层析区域(如图 1中较小的红色区域)的信号射线建立观测方程,并附加约束条件得到研究区域的层析模型,如式(4) 所示。

5) 基于MART方法对层析区域的层析模型进行解算,将第3) 步中得到的水汽密度作为MART方法解算的初值进行一步修正,得到层析结果的最终解。

图 2 辅助层析区域边长计算示意图 Figure 2 Schematic Representation of Distance Calculated for Assisted Tomographic Area
3 方法有效性检验

利用无线电探空仪获取的数据能够得到垂直方向上精确的水汽密度信息[3, 14],通常将其作为检验层析技术反演水汽结果的标准。本文将均方根误差(root mean square, RMS)、偏差(Bias)和平均绝对偏差(MAE)作为检验本文提出方法反演水汽精度的指标。

3.1 实验数据及层析策略

选取美国德克萨斯地区运行参考站(Continuously Operation Reference Stations, CORS)网中13个GPS测站(如图 3中▲所示)的观测数据进行层析实验,各测站具体站名和相关信息如表 1所示,选取时间为2015-05-10~2015-05-31共22 d,需要说明的是在此后的数据分析均是基于该数据进行的。其中,在层析区域有一个72 249探空站(radiosonde station),如图 3中·所示,该探空站每天在UTC 00:00和12:00发射探空气球获取数据,本文将该探空数据计算结果作为检核层析结果的标准。

图 3 层析区域内GPS测站和探空站地理位置分布 Figure 3 Geographic Distribution of GPS Stations and Radiosonde Station in Tomographic Area
表 1 美国得克萨斯州地区CORS网中使用的测站信息 Table 1 Information of Stations Using CORS Network in Texas
测站名 纬度/(°) 经度/(°) 高/km 接收机类型
TXCO 33.150 0 96.616 7 161.893 TRIMBLE NETR5
TXDA 32.783 3 96.666 7 160.642 TRIMBLE NETR9
TXDC 33.233 3 97.600 0 255.282 TRIMBLE NETR5
TXDE 33.200 0 97.150 0 178.789 TRIMBLE NETR5
TXES 32.366 7 96.850 0 163.749 TRIMBLE NETR9
TXGR 32.233 3 97.750 0 177.482 TRIMBLE NETR5
TXJA 33.183 3 98.133 3 326.043 TRIMBLE NETR5
TXKE 32.400 0 97.316 7 227.949 TRIMBLE NETR5
TXMW 32.800 0 98.133 3 246.447 TRIMBLE NETR5
TXSG 32.850 0 97.333 3 181.667 TRIMBLE NETR5
TXST 32.216 7 98.166 7 376.574 TRIMBLE NETR5
TXWE 32.750 0 97.816 7 337.386 TRIMBLE NETR9
ZFW1 32.816 7 97.050 0 155.193 NOV WAASGII

本文实验所选取的研究范围(如图 3中红色区域范围所示):纬度方向为32.1°N~33.3°N,分辨率为0.2°;经度方向为96.5°W~98.3°W,分辨率为0.3°;垂直方向为0~10 km,分辨率为1 km;研究区域共有6×6×10个网格。在进行层析时,选取的卫星截止高度角为10°,因此根据§2.2部分可以确定出辅助层析区域的范围(如图 3整个区域范围所示),以确保原层析区域内的13个接收机在10°~90°卫星截止高度角之间的所有射线信号都能完整穿过辅助研究区域;其具体范围:纬度31.7°N~33.7°N,经度96.2°W~98.6°W;辅助层析区域内网格个数为10×8×10。

实验中采用两种方法获取原层析区域内的水汽信息。

方法1  采用传统方法建立的层析模型反演水汽,即式(4) 建立的层析模型;

方法2  利用本文提出的层析方法反演水汽,即利用§2.2给出的解算步骤。

GAMIT (v10.5) 软件对GPS数据处理的天顶对流层延迟参数估计精度优于±1 cm[21],结合气象参数得到的PWV只有±1.5 mm,因此基于投影函数得到的SWV其精度足以满足水汽层析的需求。本文在利用GAMIT(v10.5) 软件[22]对实验数据进行处理时,为了消除局域网中各测站对流层参数相关性的影响,引入了INEG、NIST、NILB和PIE1 4个网外IGS辅助站参与解算[23]

3.2 射线利用及射线穿过网格数分析

基于§3.1给出的两种方法,统计了层析时段内每天一次不同方法射线利用情况及与有射线穿过的网格个数(图 4);表 2给出了每天的平均统计信息。由图 4表 2可以看出,方法2在射线使用条数和有射线穿过的网格数方面均大于方法1。通过计算,本文提出的方法使射线的平均利用率提高了18.94%。

图 4 两种方法射线利用情况及有射线穿过网格数统计 Figure 4 Number of Signals Used and Number of Voxels Crossed by Rays for Two Methods
表 2 两种方法射线利用情况及有射线穿过网格数统计信息 Table 2 Statistical Information of the Number of Signal Used and Voxels Crossed by Rays
方法 信号使用条数 信号穿过网格数
均值 最大值 最小值 均值 最大值 最小值
1 676 786 655 392 307 277
2 803 905 776 399 312 284
3.3 水汽反演精度检验

为了验证本文提出方法反演水汽结果的精度,对层析时段内每天两个观测历元(UTC 00:00和12:00) 不同方法反演的水汽结果进行对比。首先利用两种方法反演水汽得到探空站所在位置上相应历元的水汽密度估值,然后与探空站数据计算结果进行对比,统计得到每天的均方根误差、偏差和平均绝对误差(见图 5~图 7);表 3给出了每天的平均统计结果。

图 5 实验时段内两种方法层析结果的RMS对比 Figure 5 RMS Comparison of Tomographic Result Derived from Two Methods During Experimental Period
图 6 实验时段内两种方法层析结果的偏差对比 Figure 6 Bias Comparison of Tomographic Result Derived from Two Methods During Experimental Period
图 7 实验时段内两种方法层析结果的平均绝对偏差对比 Figure 7 MAE Comparison of Tomographic Result Derived from Two Methods During Experimental Period
表 3 与探空数据对比22天的统计信息/(g·m-3) Table 3 Statistical Information Compared with Radiosonde Data for 22 Days/(g·m-3)
方法 RMS Bias MAE
均值 最大值 最小值 均值 最大值 最小值 均值 最大值 最小值
1 1.78 2.93 0.54 0.11 1.58 -1.29 1.26 2.14 0.42
2 1.52 2.75 0.42 0.02 1.25 -1.27 1.04 1.97 0.34

图 5~图 7可以看出,方法2的RMS、Bias和MAE均优于方法1;由表 3也可以看出,方法2的平均RMS、Bias和MAE均小于方法1。这说明,本文提出的利用辅助层析区域提高射线利用率反演水汽的方法其精度要优于传统方法。

3.4 水汽密度廓线对比

为了进一步分析本文提出的方法反演水汽密度在垂直方向上的分布情况,对层析时段内的层析结果进行统计,得到探空站所在位置上不同高度的平均水汽密度值,并与探空站计算的结果进行对比(图 8)。图 8(a)给出了两种方法与无线电探空仪对比水汽密度随高度的变化情况,图 8(b)给出了不同高度上的RMS。

图 8 水汽密度廓线和RMS随高度变化情况 Figure 8 Water Vapor Density Profile and RMS Change with Height

图 8(a)可以看出,方法2反演水汽密度得到的廓线信息在不同高度上与探空数据计算结果具有更好的一致性。由图 8(b)也可以看出,方法2反演水汽的精度其垂直分布也优于传统方法。这进一步说明,通过本文方法,实现了在研究区域侧面穿出的射线信息对层析区域水汽反演的贡献,提高了水汽反演结果的质量。

4 结语

传统层析方法中,从研究区域侧面穿出的射线信息因无法被使用而导致了观测数据的浪费。针对该情况,本文提出了一种附加辅助层析区域提高射线利用率的水汽反演方法。基于美国得克萨斯地区的CORS网2015-05-10~2015-05-31共22天的观测数据进行实验,验证了本文提出方法反演水汽的有效性及精度。

实验结果表明, 本文提出的方法能够充分利用从研究区域侧面穿出的射线信息,实现了这些射线信息对最终层析结果的贡献。通过与探空数据计算结果对比发现,利用本文提出的方法得到的水汽结果有很大的改善,其RMS的改善率为14.6%,Bias的改善率为27.3%,MAE的改善率为17.5%。这说明本文提出的方法反演的水汽密度信息在RMS、Bias和MAE等方面均优于传统方法。

致谢: IGRA (Integrated Global Radiosonde Archive)提供的无线电探空数据; NOAA (National Oceanic and Atmospheric Administration)提供的实验数据。
参考文献
[1] Wang Jiuke, Han Suqin, Bian Hai, et al. Characteristics of the Three-Dimensional GPS Tomography Water Vapor Field During the Rainstorm[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2014, 50(6): 1053–1064 ( 王久珂, 韩素芹, 边海, 等. 一次暴雨过程中GPS三维层析水汽场的变化特征[J]. 北京大学学报(自然科学版), 2014, 50(6): 1053–1064. )
[2] Bevis M, Businger S, Herring T A, et al. GPS Meteorology:Remote Sensing of Atmospheric Water Vapor Using the Global Positioning System[J]. Journal of Geophysical Reseach Atmospheres, 1992, 97(D14): 15787–15801 DOI:10.1029/92JD01517
[3] He Lin, Liu Lintao, Su Xiaoqing, et al. Algebraic Reconstruction Algorithm of Vapor Tomography[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(1): 32–38 ( 何林, 柳林涛, 苏晓庆, 等. 水汽层析代数重构算法[J]. 测绘学报, 2015, 44(1): 32–38. DOI:10.11947/j.AGCS.2015.20130308 )
[4] Zhang Shuangcheng, Ye Shirong, Wan Rong, et al. Preliminary Tomography Spatial Wet Refractivity Distribution Based on Kalman Filter[J]. Geomatics and Information Science of Wuhan University, 2008, 33(8): 796–799 ( 张双成, 叶世榕, 万蓉, 等. 基于Kalman滤波的断层扫描初步层析水汽湿折射率分布[J]. 武汉大学学报·信息科学版, 2008, 33(8): 796–799. )
[5] Yu Shengjie, Liu Lintao. Application of Fitting Method by Selection of the Parameter Weights on GPS Water Vapor Tomography[J]. Geomatics and Information Science of Wuhan University, 2009, 34(1): 183–186 ( 于胜杰, 柳林涛. 利用选权拟合法进行GPS水汽层析解算[J]. 武汉大学学报·信息科学版, 2009, 34(1): 183–186. )
[6] Yu Shengjie, Liu Lintao, Liang Xinghui. Influence Analysis of Constraint Conditions on GPS Water Vapor Tomography[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(5): 491–496 ( 于胜杰, 柳林涛, 梁星辉. 约束条件对GPS水汽层析解算的影响分析[J]. 测绘学报, 2010, 39(5): 491–496. )
[7] Ye Shirong, Jiang Peng, Liu Yanyan. A Water Vapor Tomographic Numerical Quadrature Approach with Ground Based GPS Network[J]. Acta Geodaetica et Cartographica Sinica, 2013, 42(5): 654–660 ( 叶世榕, 江鹏, 刘炎炎. 地基GPS网层析水汽三维分布数值积分方法[J]. 测绘学报, 2013, 42(5): 654–660. )
[8] Chen B Y, Liu Z Z. Voxel-Optimized Regional Water Vapor Tomography and Comparison with Radiosonde and Numerical Weather Model[J]. Journal of Geodesy, 2014, 88: 691–703 DOI:10.1007/s00190-014-0715-y
[9] Xia Pengfei, Cai Changsheng, Dai Wujiao, et al. Three-Dimension Water Vapor Tomography Using Ground-Based GPS and COSMIC Occupation Observation[J]. Geomatics and Information Science of Wuhan University, 2013, 38(8): 892–896 ( 夏鹏飞, 蔡昌盛, 戴吾蛟, 等. 地基GPS联合COSMIC掩星数据的水汽三维层析研究[J]. 武汉大学学报·信息科学版, 2013, 38(8): 892–896. )
[10] Niell A E, Coter A J, Solheim F S, et al. Comparison of Measurements of Atmospheric Wet Delay by Radiosonde, Water Vapor Rdiometer, GPS, and VLBI[J]. Journal of Atmospheric and Oceanic Technology, 2001, 18(6): 830–850 DOI:10.1175/1520-0426(2001)018<0830:COMOAW>2.0.CO;2
[11] Troller M B, Urki B, Cocard M, et al. 3-D Refractivity Field from GPS Double Difference Tomography[J]. Geophys Res Lett, 2002, 29(24): 2149–2152
[12] Rohm W, Bosy J. Local Tomography Troposphere Model over Mountains Area[J]. Geophysical Research Letteus, 2009 DOI:10.1016/j.atmosres.2009.03.013
[13] Bender M, Stosius R, Zus F, et al. GNSS Water Vapor Tomography-Expected Improvements by Combing GPS, GLONASS and Galileo Observations[J]. Advances in Space Research, 2010, 47(5): 886–897
[14] Adeyemi B, Joerg S. Analysis of Water Vapor over Nigeria Using Radiosonde and Satellite Data[J]. Journal Applied Meteorology and Climatology, 2012, 51: 1855–1866 DOI:10.1175/JAMC-D-11-0119.1
[15] Liu Z Z, Wong M S, Nichol J, et al. A Multi-Sensor Study of Water Vapour from Radiosonde, MODIS and AERONET:A Case Study of Hong Kong[J]. International Journal of Climatology, 2013, 33(1): 109–120 DOI:10.1002/joc.v33.1
[16] Flores A, Ruffini G, Rius A. 4D Tropospheric Tomography Using GPS Slant Wet Delays[J]. Annales Geophysica, 2000, 18(2): 223–234 DOI:10.1007/s00585-000-0223-7
[17] Bender M, Dick G, Ge M, et al. Development of a GNSS Water Vapour Tomography System Using Algebraic Reconstruction Techniques[J]. Adv. Space Res., 2011, 47(10): 1704–1720 DOI:10.1016/j.asr.2010.05.034
[18] Bi Y M, Mao J, Li C. Preliminary Results of 4D Water Vapor Tomography in the Troposphere using GPS[J]. Advaucesin Atmospheric Sciences, 2006, 23: 551–560 DOI:10.1007/s00376-006-0551-y
[19] Elósegui P A, Ruis J L, Davis G, et al. An Experiment for Estimation of the Spatial and Temporal Variations of Water Vapor Using GPS Data[J]. Physics and Chemistry of the Eart, 1998, 23(1): 125–130 DOI:10.1016/S0079-1946(97)00254-1
[20] Ran B R, Ge W Z. Singular Value Decomposition Method Compared with Damping Least Square Method[J]. Geophysical Computing Technology, 1997, 1: 46–49
[21] Gusfarienza H, Yuwono B D, Awaluddin M, et al. Penentuan Zenith Tropospheric Delay and Precipitation Vapor Menggunakan Pejragkat Lunak GAMIT[J]. Jurnal Geodesi Undip, 2015, 4(2): 78–86
[22] Herring T A, King R W, McClusky S C. Documentation of the GAMIT GPS Analysis Software release 10.4[OL]. http://www-gpsg.mit.edu/~simon/gtgk/GAMIT_Ref.pdf, 2010
[23] Rocken C, Hove T V, Johnson J, et al. GPS/STORM-GPS Sensing of Atmospheric Water Vapor for Meteorology[J]. Journal of Atmospheric and Oceanic Technology, 1995, 12(3): 468–478 DOI:10.1175/1520-0426(1995)012<0468:GSOAWV>2.0.CO;2