-
摘要: 对流层延迟是全球导航卫星系统(GNSS)计算的主要误差之一,其模型精度对测站坐标解算有较大影响,在高程方向尤为明显。因此,有必要对不同的对流层延迟改正模型的适用性进行评估。采用SHA解算了中国陆态网GNSS跟踪站的对流层天顶延迟数据,对常用的对流层改正模型EGNOS/UNB3m/GPT/GPT2的天顶延迟量在中国不同区域、不同季节的适用性进行了分析。结果显示,4种模型的RMS均为4~5 cm,各模型RMS之差小于1 cm,其中GPT2模型的RMS最小;4种模型的平均偏差(BIAS)为1 cm左右,GPT2模型的BIAS最大,为1.5 cm;时间上,各个模型在夏季精度普遍较低,这是因为夏季水汽丰富,对流层湿延迟变化较大;空间上,各模型在东南沿海精度较低,因为东南沿海气候湿润,湿延迟变化较大;各模型精度对测站高程不敏感,精度在比较高的测站并无明显降低。通过对不同模型在中国区域的精度分析,验证该改正模型可以为中国区域用户的对流层模型的选择提供一定的参考。Abstract: As one of the main error sources in the Global Navigation Satellite System (GNSS), the accuracy of tropospheric delay correction model influences the estimated coordinates, especially in vertical direction; therefore, it is necessary to make assessment of different tropospheric delay correction models. An assessment of four commonly used models, including EGNOS/UNB3m/GPT/GPT2, is made in this paper. The result shows that the RMS of the four models remains 4-5 cm, and the difference between each model is less than 1 cm, with GPT2 being the most accurate model; and the RMS of each model in summer is relatively much bigger than in winter, due to the abundant water vapor in summer, which makes it harder to model tropospheric delay precisely. All models show poor accuracy in southeast China, thanks to the rather changeable weather and the abundant water vapor there; and no model is sensitive to the height, since the accuracy varies little with respect to the variation of altitude. These findings provide a reference for the GNSS users in China area when choosing a suitable tropospheric delay correction model.
-
近年来,随着各类GPS监测网络的建立,尤其是中国地壳运动观测网络(简称“网络工程”)和中国大陆构造环境监测网络(简称“陆态网络”)的投入运行以来,产生了大量的GPS连续和流动观测数据。依据这些数据得到的GPS速度场及其衍生产品在地壳运动监测、地球动力学研究等领域发挥了重要作用[1-7]。在上述研究中,GPS流动水平速度场应用较多,而GPS流动垂向速度场受非构造运动影响较大且缺乏足够的观测数据来消除这些影响,使得GPS流动观测解算的垂向速度场精度不高,应用较少。
目前,我国部分区域(如京津唐、川滇地区)GPS连续观测网络已初具规模,如何利用这些较高密度的连续观测网络提高同一区域的流动观测垂向速率精度,是挖掘GPS垂向速度场价值的当务之急,具有重要的意义。鉴于此,本文提出了利用GPS连续观测修正流动观测垂向速率的方法,并选用云南地区陆态网络GPS连续和流动观测数据进行了验证。
1 GPS连续观测年周期信号提取与应用
GPS垂向时间序列存在较为显著的季节性运动,且年周期运动最为突出[8-9],因此,本文只考虑年周期运动,则有:
(1) 式中,Y为垂向时间序列;a为常数项;V0为速率;A和B为年周期运动正弦波和余弦波振幅;C为年周期运动总振幅;φ为年周期运动初相;ε为残差;t为时间。
由于网络工程和陆态网络采用的流动观测模式为每隔2~3 a观测一期,每期连续观测4 d,这种观测模式难以获取准确的年周期运动信息,因此,求解时只能顾及线性运动:
(2) 在用最小二乘理论求解参数时,如果忽略部分参数,则估计参数的平差值会包含未估计参数的影响[8],即按式(2)求得的速率V1包含了年周期运动的影响,与V0存在一定偏差。如果能够获取每个流动站的年周期运动振幅并在时间序列Y中扣除年周期运动得到Y′,再根据Y′只顾及线性运动求解速率V2,则V2应该能够扣除年周期运动的影响,等于(或者接近于)V0:
(3) 由于一些区域内GPS时间序列中的非构造噪声存在区域相关性[10-11],因此,可以利用同一个区域的连续站年周期运动信息通过空间加权内插得到流动站的年周期运动信息。对于每个流动站,年周期运动振幅A0和B0分别取其周边200 km内所有连续站的年周期运动振幅A和B的加权平均(权取距离d的平方倒数):
(4) 2 连续站数据检核
2.1 数据概述
为了检核修正方法,选取陆态网络GPS连续站密度较高的云南地区进行试验。陆态网络在该地区布设了27个连续站,其中KMIN和XIAG站从1999年开始观测,其他25个测站是从2010年陆续开始观测。本文选取这27个测站2010-01-01以后的观测数据,至2013年底这些测站观测时间均超过960 d,能够较为准确地获取该地区垂向速率和年周期运动信息。利用GAMIT/GLOBK[12]和QOCA[8]软件解算各连续站观测数据,得到各连续站2010~2013年的时间序列,然后按式(1)分别求解各连续站垂向运动速率和年周期运动振幅(表 1和图 1)。
表 1 云南地区GPS连续站垂向运动信息(2010~2013年)Table 1. Vertical Motion for GPS Fiducial Station in Yunnan Area from 2010 to 2013点名 速率 /(mm\5a-1) 年周期 观测天数 振幅/mm 初相/(°) KMIN 2.93 9.45 17.61 1 305 XIAG 0.17 9.5 357.3 1 365 YNCX* 1.52 8.55 4.11 1 108 YNDC* 2.21 9.44 9.89 1 050 YNGM -0.24 6.35 350.08 988 YNHZ 1.21 8.29 8.56 973 YNJD 1.04 10.97 8.31 1126 YNJP -1.9 8.88 17.38 1 134 YNLA -0.68 14.31 9.05 1 147 YNLC* 0.6 11.49 5.49 1 008 YNLJ* 1.39 10.56 18.84 1 115 YNMH -1.02 11.6 7.86 1 150 YNMJ -0.3 9.89 0.15 1 129 YNML 1.07 7 12.25 1 149 YNMZ* -3.17 10.65 24.8 1 163 YNRL -0.64 11.8 10.57 1 087 YNSD* 2.18 12.27 11.65 975 YNSM* -1.54 10.94 5.44 1 162 YNTC 0.84 12.73 6.87 1 149 YNTH* 1.16 8.27 17.31 1 164 YNWS 3.17 8.47 21.88 968 YNXP 2.15 7.95 6.63 1 215 YNYA 2.38 9.68 14.57 1 006 YNYL 1.95 12.85 16.35 1 163 YNYM* 1.45 10.22 17.14 1 121 YNYS 2.11 11.11 12.93 1 160 YNZD 0.86 8.97 22.98 980 注: 标*的点为检核站,未标注的点为基准站。 表 1是云南地区GPS连续站2010~2013年垂向运动信息,从中可以看出该地区各GPS连续站年周期运动较为显著,年周期运动振幅较大(7~14 mm),垂向运动速率相对较小(4 mm/a以内),因此,测站年周期运动会对垂向运动速率的求解产生较大的影响,在该地区有必要对流动观测垂向运动进行修正。
图 1是云南地区GPS连续站2010~2013年年周期运动,时间参考点为2010.0,箭头长度表示振幅,从正东逆时针旋转的方位角表示初相位方向,箭头指向东、南、西、北方向依次表示极值发生在0.25、0.50、0.75、0.0 年。结合表 1和图 1可以看出云南地区GPS连续站的年周期运动空间一致性较好,各连续站年周期运动振幅和初相均比较一致,满足本方法的使用条件。
2.2 检核结果及分析
在云南地区27个连续站中选取18个测站(图 1中蓝色点)记为基准站,9个测站(表 1标*点和图 1中黑色点)记为检核站,用以检核修正方法效果。首先根据每个连续站垂向时间序列Y按式(1)用最小二乘方法求解速率V0和年周期运动振幅A和B;然后根据每个检核站垂向时间序列Y按式(2)用最小二乘方法求解速率V1;再次,每个检核站按式(4)根据基准站的年周期运动振幅插值求得年周期运动振幅A0和B0,并在时间序列Y中扣除依据A0和B0计算得到的年周期运动,得到时间序列Y′,然后按式(3)用最小二乘方法求解速率V2。比较V1和V2与V0的差异(见表 2),如果V2与V0的差异比V1与V0的差异小,则说明修正方法有效。
表 2 检核站解算结果Table 2. Results of Check Sites检核站 年周期运动振幅/mm 速率结果/(mm\5a-1) A B A0 B0 V0 V1 V2 V1-V0 V2-V0 YNCX 8.53 0.61 9.55 1.66 1.52±0.12 2.04±0.18 1.55±0.12 0.51 0.03 YNDC 9.30 1.62 8.25 1.39 2.21±0.19 1.21±0.25 2.09±0.19 -1.01 -0.12 YNLC 11.44 1.10 8.95 0.23 0.6±0.15 1.21±0.26 0.71±0.16 0.61 0.11 YNLJ 9.99 3.41 10.43 2.38 1.39±0.11 1.29±0.2 1.24±0.1 -0.10 -0.15 YNMZ 9.66 4.47 8.17 2.36 -3.17±0.16 -3.01±0.23 -3.22±0.16 0.16 -0.06 YNSD 12.02 2.48 11.00 1.28 2.18±0.17 -0.22±0.29 1.86±0.17 -2.40 -0.32 YNSM 10.89 1.04 10.81 1.01 -1.54±0.12 -1.14±0.21 -1.54±0.12 0.40 0.00 YNTH 7.90 2.46 8.09 1.63 1.16±0.11 1.34±0.17 1.1±0.12 0.18 -0.05 YNYM 9.76 3.01 9.36 2.10 1.45±0.13 1.31±0.21 1.34±0.13 -0.14 -0.11 从表 2可以看出,不顾及年周期运动求得的速率(V1)与顾及年周期运动求得的速率(V0)差别较大,最大值达到-2.4 mm/a,说明在年周期运动较为显著的地区,在求解线性速率时必须考虑年周期运动。由于该区域年周期运动空间一致性较好,9个检核站内插得到的年周期运动振幅比较接近于实际值(差异在3 mm以内)。由于内插得到的检核站年周期运动振幅精度较高,使得本方法对垂向速率的修正效果较好,各检核站经过年周期运动修正后得到的速率(V2)比未经过修正的速率(V1)更为接近V0,9个检核站V2与V0的差异量值在0.32 mm/a以内,且V2的误差比V1小,接近V0的误差。以上结果说明,在垂向时间序列年周期运动显著且空间一致性较好的区域,通过本方法能够较好地消除年周期运动对垂向速率的影响,提高垂向速率的解算精度。
3 流动站修正结果及分析
陆态网络在云南地区布设了157个GPS流动站,分别于2011年3~8月和2013年4~9月进行了一期观测,每期观测时段长度为4 d。利用GAMIT/GLOBK和QOCA软件解算了这两期流动观测数据(解算模式与§2中GPS连续站时间序列的解算模式保持一致,即GAMIT解算参数、QOCA解算中的框架点的选取等设置一致),得到了各流动站2011~2013年的时间序列,按式(2)求得了各流动站未扣除年周期运动的2011~2013年垂向速率(图 2(a))。利用§2计算得到的云南地区27个连续站年周期运动振幅,按式(4)内插得到各流动站年周期运动振幅,然后对各流动站的时间序列扣除年周期运动,再按式(3)求得扣除年周期运动的2011~2013年垂向速率(图 2(b))。
由图 2可知,未进行年周期改正的垂向速度场(图 2(a))和改正后的垂向速度场(图 2(b))整体趋势类似,即大体以北纬24°为界,北纬24°以北的区域以隆升为主,北纬24°以南以沉降为主,大体与由GPS连续站测定的云南地区垂向速度场(图 3)特征一致。改正后的垂向速度场比未改正的速度场在空间一致性上有所改善,且更加接近GPS连续站测定的垂向速度场和水准观测结果得到的该区地壳垂直运动特征[13-14],如图 2中两个红色椭圆框内区域。未改正的速度场显得较为混乱,方向上有升也有降,量值有大有小,改正后的速度场基本上都是上升,且量值空间分布更为均匀。图 2中绿色椭圆框区域的速度场经过改正后更接近郝明等的文章中图 7结果[14]。总体而言,本方法对GPS流动垂向速度场的改进取得了较好的效果。
4 结 语
本文提出了利用连续观测站年周期运动信息通过空间插值获取流动观测站年周期运动信息、并在流动观测站时间序列中加以扣除、进而求解垂向速率的方法,并利用云南地区陆态网络2010~2013年连续观测数据和2011~2013年流动观测数据进行了验证,结果表明在年周期运动空间一致性较好的区域(云南地区)能够取得较好的效果。
本方法的前提是区域内年周期运动空间一致性较好,只有满足这个条件,才能利用连续站年周期运动信息内插得到较为准确的流动站年周期运动振幅,才能得到较为准确的垂向速率。由于本文的主要目的是验证利用连续站信息修正流动观测垂向速率的方法有效性,因此,只考虑了季节性运动中最主要的年周期运动,但本方法也适用于其他季节性运动(如半年周期运动),只需在式(1)和式(3)增加相关参数求取并加以修正。
我国国土辽阔,各区域情况不尽相同,可先处理区域内连续观测数据,分析区域垂向季节性运动机制。如果与趋势(线性)运动相比,季节性运动振幅较小,对垂向速率的影响可忽略不计,则无需进行季节性改正;如果区域内季节性运动影响较大,且有多种较为显著的周期运动(年周期、半年周期),则可以根据分析结果选取对垂向速率影响较大的季节项加以修正。为了进一步提高方法的准确性,还将进行空间插值方法的择优选取、在GPS连续站垂向季节性运动信息求取前扣除一些能够定量计算的非构造形变(如大气、积雪等负荷效应造成的非构造形变)[15-18])等分析研究。
-
表 1 SHA中国区域IGS测站对流层产品精度/mm
Table 1 Accuracy of SHA ZTD Comparing to IGS Solution/mm
BJFS SHAO CHAN LHAZ URUM KUNM TWTF MEAN BIAS 1.2 -0.7 1.5 1.5 1.0 1.0 1.0 0.9 RMS 2.0 2.0 2.3 2.6 1.6 3.8 1.6 2.3 表 2 不同模型在中国区域的精度/cm
Table 2 Accuracy of Different Models/cm
EGNOS UNB3m GPT GPT2 RMS 5.59 [1.62,8.20] 5.16 [1.52,7.79] 5.00 [1.71,11.49] 4.65 [1.42,9.56] BIAS -0.59 [-5.85,4.23] -0.33 [-5.04,3.59] 0.74 [-7.96,9.45] 1.53[-3.88,7.71] 注:“[,]”表示最小值与最大值。 -
[1] 王潜心, 许国昌, 陈正阳. 利用区域GPS网进行高海拔流动站的对流层延迟量内插[J]. 武汉大学学报·信息科学版, 2010, 35(12):1405-1408 Wang Qianxin, Xu Guochang, Chen Zhengyang. Interpolation Method of Tropospheric Delay of High Altitude Rover Based on Regional GPS Network[J].Geomatics and Information Science of Wuhan University, 2010, 35(12):1405-1408
[2] 戴吾蛟, 陈招华, 匡翠林, 等. 区域精密对流层延迟建模[J]. 武汉大学学报·信息科学版, 2011, 36(4):392-396 http://ch.whu.edu.cn/CN/abstract/abstract513.shtml Dai Wujiao, Chen Zhaohua, Kuang Cuilin, et al. Modeling Regional Precise Tropospheric Delay[J].Geomatics and Information Science of Wuhan University, 2011, 36(4):392-396 http://ch.whu.edu.cn/CN/abstract/abstract513.shtml
[3] 姚宜斌, 张瑞, 易文婷. 一种新的区域对流层拟合模型及其在PPP中的应用[J]. 武汉大学学报·信息科学版, 2012, 37(9):1024-1027 http://ch.whu.edu.cn/CN/abstract/abstract306.shtml Yao Yibin, Zhang Rui, Yi Wenting, et al. A New Regioinal Troposphere Fitting Model and Its Application to PPP[J]. Geomatics and Information Science of Wuhan University, 2012, 37(9):1024-1027 http://ch.whu.edu.cn/CN/abstract/abstract306.shtml
[4] 张双成, 叶世榕, 刘经南, 等. 动态映射函数最新进展及其在GNSS遥感水汽中的应用研究[J]. 武汉大学学报·信息科学版, 2009, 34(3):280-283 http://ch.whu.edu.cn/CN/abstract/abstract1190.shtml Zhang Shuangcheng, Ye Shirong, Liu Jingnan, et al. Latest Progress of Dynamic Mapping Functions and Its Application to GNSS Retrieved Water-Vapor[J].Geomatics and Information Science of Wuhan University, 2009, 34(3):280-283 http://ch.whu.edu.cn/CN/abstract/abstract1190.shtml
[5] Niell A E. Global Mapping Functions for the Mapping of the Atmospheric Delay at Radio Wavelengths[J]. Journal of Geophysical Research, 1996, 101(B2):3227-3246 doi: 10.1029/95JB03048
[6] Johannes B, Arthur N, Paul T, et al. The Global Mapping Function (GMF):A New Empirical Mapping Function Based on Numerical Weather Model Data[J]. Geophysical Research Letters, 2006, 33:L07304 http://cn.bing.com/academic/profile?id=2134691634&encoded=0&v=paper_preview&mkt=zh-cn
[7] Kouba J. Implementation and Testing of the Gridded Vienna Mapping Function 1(VMF1)[J]. J Geod, 2008, 82:193-205 doi: 10.1007/s00190-007-0170-0
[8] 殷海涛, 黄丁发, 熊永良,等. GPS信号对流层延迟改正新模型研究[J]. 武汉大学学报·信息科学版, 2007,32(5):454-457 http://ch.whu.edu.cn/CN/abstract/abstract1884.shtml Yin Haitao, Huang Dingfa, Xiong Yongliang, et al. New Model for Tropospheric Delay Estimation of GPS Signal[J].Geomatics and Information Science of Wuhan University, 2007,32(5):454-457 http://ch.whu.edu.cn/CN/abstract/abstract1884.shtml
[9] Vedel H, Mogensen K S, Huang X Y. Calculation of Zenith Delays from Meteorological Data Comparison of NWP Model, Radiosonde and GPS Delays[J]. Phys Chem Earth (A), 2001, 26(6-8):497-502 doi: 10.1016/S1464-1895(01)00091-6
[10] Andrei C, Chen R. Assessment of Time-series of Troposphere Zenith Delays Derived from the Global Data Assimilation System Numerical Weather Model[J]. GPS Solution, 2008, 13(2):109-117 http://cn.bing.com/academic/profile?id=1988114896&encoded=0&v=paper_preview&mkt=zh-cn
[11] Chen Qinming, Song Shuli, Heise S, et al. Assessment of ZTD Derived from ECMWF/NCEP Data with GPS ZTD over China[J]. GPS Solution, 2011, 15:415-425 doi: 10.1007/s10291-010-0200-x
[12] 陈钦明, 宋淑丽, 朱文耀. 亚洲地区ECMWF/NCEP资料计算ZTD的精度分析[J]. 地球物理学报, 2012, 55(5):1541-1548 Chen Qinming, Song Shuli, Zhu Wenyao. An Analysis of Accuracy of Zenith Tropospheric Delay Calculated from ECMWF/NCEP Data over Asian Area[J].Chinese J. Geophys, 2012, 55(5):1541-1548
[13] Saastamoinen J. Contribution to the Theory of Atmospheric Refraction Part Ⅱ, Refraction Corrections in Satellite Geodesy[J].Bulletin Geodesique, 1972, 107:13-34 doi: 10.1007%2FBF02521844
[14] Penna N, Dodson A, Chen W. Assessment of EGNOS Tropospheric Correction Model[J]. Journal of Navigation, 2000, 54(1):37-55 http://cn.bing.com/academic/profile?id=2097918869&encoded=0&v=paper_preview&mkt=zh-cn
[15] Collins J P, Langley R B. The Residual Tropospheric Propagation Delay:How Bad Can It Get?[C]. ION GPS 1998, Nashville, Tennessee, USA, 1998
[16] Leandro R F,Santos M C, Langley R B. UNB Neutral Atmosphere Models:Development and Performance[C]. ION NTM 2006, California, USA, 2006 http://www.oalib.com/references/18988338
[17] B hm J R. Heinkelmann H S. Short Note:A Global Model of Pressure and Temperature for Geodetic Applications[J]. J. Geod, 2007, 81(10):679-683 doi: 10.1007/s00190-007-0135-3
[18] Lagler K, Schindelegger M, B hm J, et al. GPT2:Empirical Slant Delay Model for Radio Space Geodetic Techniques[J]. Geophysical Research Letters, 2013,40(6):1069-1073 doi: 10.1002/grl.50288
[19] 毛健, 朱长青, 郭继发. 一种新的全球对流层天顶延迟模型[J]. 武汉大学学报·信息科学版, 2013,38(6):684-688 http://ch.whu.edu.cn/CN/abstract/abstract2667.shtml Mao Jian, Zhu Changqing, Guo Jifa. A New Global Zenith Tropospheric Delay Model[J].Geomatics and Information Science of Wuhan University, 2013,38(6):684-688 http://ch.whu.edu.cn/CN/abstract/abstract2667.shtml
[20] 李薇, 袁运斌, 欧吉坤, 等. 全球对流层天顶延迟模型IGGtrop的建立于分析[J]. 科学通报, 2012, 57(15):1317-1325 http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201215006.htm Li Wei, Yuan Yunbin, Ou Jikun, et al. A New Global Zenith Tropospheric Delay Model IGGtrop for GNSS Applications[J]. Chin Sci Bull, 2012, 57:1317-1325 http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201215006.htm
[21] 姚宜斌, 何畅勇, 张豹. 一种新的全球对流层天顶延迟模型GZTD[J]. 地球物理学报, 2013, 56(7):2218-2227 http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201307008.htm Yao Yibin, He Changyong, Zhang Bao. A New Global Tropospheric Delay Model GZTD[J]. Chinese J. Geophys, 2013, 56(7):2218-2227 http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201307008.htm
[22] Song S L, Zhu W Y, Chen Q M, et al. Establishment of a New Tropospheric Delay Correction Model over China Area[J]. Science China(Physics, Mechanics and Astronomy), 2011, 54(12):2271-2283 doi: 10.1007/s11433-011-4530-7
[23] 曲伟菁, 朱文耀, 宋淑丽. 三种对流层延迟改正模型精度评估[J]. 天文学报,2008,49(1):113-122 http://www.cnki.com.cn/Article/CJFDTOTAL-TWXB200801013.htm Qu Weijing, Zhu Wenyao, Song Shuli, et al. The Evaluation of Precision About Hopfield, Saastamoinen and EGNOS Tropospheric Delay Correction Model[J].Acta Astronmica Sinica, 2008,49(1):113-122 http://www.cnki.com.cn/Article/CJFDTOTAL-TWXB200801013.htm
[24] 杨玲, 李博峰, 楼立志. 不同对流层模型对GPS定位结果的影响[J]. 测绘通报, 2009, 4:9-11 http://www.cnki.com.cn/Article/CJFDTOTAL-CHTB200904003.htm Yang Ling, Li Bofeng, Lou Lizhi. Effects of Different Troposphere Correction Models on GPS Position[J].Bulletin of Surveying and Mapping, 2009, 4:9-11 http://www.cnki.com.cn/Article/CJFDTOTAL-CHTB200904003.htm
[25] 王君刚, 陈俊平, 王解先. GNSS对流层延迟映射模型分析[J]. 天文学进展, 2014, 32(3):383-394 http://www.cnki.com.cn/Article/CJFDTOTAL-TWJZ201403007.htm Wang Jungang, Chen Junping, Wang Jiexian.Analysis of Tropospheric Propagation Delay Mapping Function Models in GNSS[J]. Progress in Astronomy, 2014, 32(3):383-394 http://www.cnki.com.cn/Article/CJFDTOTAL-TWJZ201403007.htm
[26] Chen Junping, Wu Bin, Hu Xiaogong, et al. SHA:the GNSS Analysis Center at SHAO[J].Lecture Notes in Electrical Engineering,2012,160:213-221 doi: 10.1007/978-3-642-29175-3
[27] Chen Junping, Zhang Yize, Zhou Xuhua, et al. GNSS Clock Corrections Densification at SHAO:from 5 Minutes to 30 Seconds[J]. Science China Physics, Mechanics & Astronomy, 2013,6:1-10 http://cn.bing.com/academic/profile?id=2267711533&encoded=0&v=paper_preview&mkt=zh-cn
[28] Jin S G, Park J U, Cho J H, et al. Seasonal Variability of GPS-derived Zenith Tropospheric Delay (1994-2006) and Climate Implications[J]. J Geophys Res, 2007, 112:D09110 http://cn.bing.com/academic/profile?id=1994554146&encoded=0&v=paper_preview&mkt=zh-cn