留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

附加辅助层析区域提高射线利用率的水汽反演方法

赵庆志 姚宜斌 罗亦泳

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

附加辅助层析区域提高射线利用率的水汽反演方法

doi: 10.13203/j.whugis20150592
基金项目: 

国家自然科学基金 41174012

国家自然科学基金 41274022

江西省自然科学基金 20151BAB213031

详细信息
    作者简介:

    赵庆志, 博士, 主要从事GNSS数据处理及GNSS气象学研究。zhaoqingzhia@163.com

    通讯作者: 姚宜斌, 博士, 教授。ybyao@whu.edu.cn
  • 中图分类号: P228

A Method to Improve the Utilization of Observation for Water Vapor Tomography by Adding Assisted Tomographic Area

Funds: 

The National Natural Science Foundation of China 41174012

The National Natural Science Foundation of China 41274022

the Jiangxi Province Natural Science Foundation of China 20151BAB213031

More Information
    Author Bio:

    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
  • 摘要: 在已有的全球导航卫星系统水汽反演技术中,能够利用的观测值是完整穿过整个层析区域的卫星射线。由于卫星、接收机几何位置分布以及层析区域选择的特定性,许多卫星射线是从层析区域侧面穿出,这些射线通常当作无效信息被剔除,降低了层析结果的精度。针对该缺点,提出并实现了一种附加辅助层析区域提高射线利用率的水汽反演方法,使从研究区域侧面穿出的信号信息也可以被利用。基于美国得克萨斯州(Texas)地区的连续运行参考站(Continuously Operation Reference Stations,CORS)网的实测数据,结合层析区域内的无线电探空仪数据进行实验,结果表明,该方法能够提高水汽反演结果的精度,其均方根误差的改善率为14.6%。
  • 图  1  层析区域三维分布图

    Figure  1.  Three-Dimensional Distribution of Tomographic Area

    图  2  辅助层析区域边长计算示意图

    Figure  2.  Schematic Representation of Distance Calculated for Assisted Tomographic Area

    图  3  层析区域内GPS测站和探空站地理位置分布

    Figure  3.  Geographic Distribution of GPS Stations and Radiosonde Station in Tomographic Area

    图  4  两种方法射线利用情况及有射线穿过网格数统计

    Figure  4.  Number of Signals Used and Number of Voxels Crossed by Rays for Two Methods

    图  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

    图  8  水汽密度廓线和RMS随高度变化情况

    Figure  8.  Water Vapor Density Profile and RMS Change with Height

    表  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
    下载: 导出CSV

    表  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
    下载: 导出CSV

    表  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
    下载: 导出CSV
  • [1] 王久珂, 韩素芹, 边海, 等.一次暴雨过程中GPS三维层析水汽场的变化特征[J].北京大学学报(自然科学版), 2014, 50(6):1053-1064 http://www.cnki.com.cn/Article/CJFDTOTAL-BJDZ201406009.htm

    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 http://www.cnki.com.cn/Article/CJFDTOTAL-BJDZ201406009.htm
    [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] 何林, 柳林涛, 苏晓庆, 等.水汽层析代数重构算法[J].测绘学报, 2015, 44(1):32-38 doi:  10.11947/j.AGCS.2015.20130308

    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 doi:  10.11947/j.AGCS.2015.20130308
    [4] 张双成, 叶世榕, 万蓉, 等.基于Kalman滤波的断层扫描初步层析水汽湿折射率分布[J].武汉大学学报·信息科学版, 2008, 33(8):796-799 http://ch.whu.edu.cn/CN/abstract/abstract1683.shtml

    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 http://ch.whu.edu.cn/CN/abstract/abstract1683.shtml
    [5] 于胜杰, 柳林涛.利用选权拟合法进行GPS水汽层析解算[J]武汉大学学报·信息科学版, 2009, 34(1):183-186 http://ch.whu.edu.cn/CN/abstract/abstract125.shtml

    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 http://ch.whu.edu.cn/CN/abstract/abstract125.shtml
    [6] 于胜杰, 柳林涛, 梁星辉.约束条件对GPS水汽层析解算的影响分析[J].测绘学报, 2010, 39(5):491-496 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201005011.htm

    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 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201005011.htm
    [7] 叶世榕, 江鹏, 刘炎炎.地基GPS网层析水汽三维分布数值积分方法[J].测绘学报, 2013, 42(5):654-660 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201305006.htm

    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 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201305006.htm
    [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] 夏鹏飞, 蔡昌盛, 戴吾蛟, 等.地基GPS联合COSMIC掩星数据的水汽三维层析研究[J]武汉大学学报·信息科学版, 2013, 38(8):892-896 http://ch.whu.edu.cn/CN/abstract/abstract2724.shtml

    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 http://ch.whu.edu.cn/CN/abstract/abstract2724.shtml
    [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 http://adsabs.harvard.edu/abs/2002GeoRL..29.2149T
    [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 http://www.sciencedirect.com/science/article/pii/S0273117710006125
    [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 Earth, 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
  • [1] 徐晓华, 张纪满, 罗佳, 高攀.  FY-3C无线电掩星折射率廓线的反演及验证 . 武汉大学学报 ● 信息科学版, 2022, 47(1): 36-44. doi: 10.13203/j.whugis20190361
    [2] 赵庆志, 姚宜斌, 辛林洋.  融合ECMWF格网数据的水汽层析精化方法 . 武汉大学学报 ● 信息科学版, 2021, 46(8): 1131-1138. doi: 10.13203/j.whugis20190323
    [3] 罗玲, 毛德华, 张柏, 王宗明, 杨桄.  基于Landsat OLI影像的典型湿地芦苇植被光能利用率遥感反演方法初探 . 武汉大学学报 ● 信息科学版, 2020, 45(4): 524-533. doi: 10.13203/j.whugis20180294
    [4] 姚宜斌, 赵庆志, 罗亦泳.  附加虚拟信号精化水汽层析模型的方法 . 武汉大学学报 ● 信息科学版, 2017, 42(11): 1658-1664. doi: 10.13203/j.whugis20150444
    [5] 王勇, 刘严萍, 李江波, 柳林涛.  GPS和无线电探空的水汽变化与PM2.5/PM10变化的相关性研究 . 武汉大学学报 ● 信息科学版, 2016, 41(12): 1626-1631. doi: 10.13203/j.whugis20140628
    [6] 袁鹏, 孙宏飞, 秦昌威, 张理想.  安徽CORS参考站三维速度场分析 . 武汉大学学报 ● 信息科学版, 2016, 41(4): 535-540. doi: 10.13203/j.whugis20140552
    [7] 王永前, 施建成, 刘志红, 冯文兰.  利用微波辐射计AMSR-E的京津冀地区大气水汽反演 . 武汉大学学报 ● 信息科学版, 2015, 40(4): 479-486. doi: 10.13203/j.whugis20130530
    [8] 钱 闯, 刘 晖, 丁志刚, 郑洪艳.  顾及非构造形变的参考站长期稳定性分析 . 武汉大学学报 ● 信息科学版, 2015, 40(9): 1259-1265. doi: 10.13203/j .whu g is20130715
    [9] 张熙, 黄丁发, 廖华, 周乐韬, 张德强.  CORS网型结构对网络PTK服务性能的影响研究 . 武汉大学学报 ● 信息科学版, 2015, 40(7): 887-893. doi: 10.13203/j.whugis20130627
    [10] 姜卫平, 袁鹏, 田挚, 肖玉钢.  区域CORS组网中的坐标基准统一方法 . 武汉大学学报 ● 信息科学版, 2014, 39(5): 566-570. doi: 10.13203/j.whugis20130383
    [11] 张小红, 朱锋, 李盼, 翟广.  区域CORS网络增强PPP天顶对流层延迟内插建模 . 武汉大学学报 ● 信息科学版, 2013, 38(6): 679-683.
    [12] 李玮, 程鹏飞, 秘金钟.  利用非组合精密单点定位提取区域电离层延迟及其精度评定 . 武汉大学学报 ● 信息科学版, 2011, 36(10): 1200-1203.
    [13] 刘经南, 赵莹, 张小红.  GNSS无线电掩星电离层反演技术现状与展望 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 631-635.
    [14] 郑作亚, 卢秀山, 张晨晰, 于志刚.  利用CORS的3DSurs系统高精度定位技术研究 . 武汉大学学报 ● 信息科学版, 2010, 35(9): 1039-1043.
    [15] 李健, 李建伟.  一种面向大规模基准站网的网格计算理论模型 . 武汉大学学报 ● 信息科学版, 2009, 34(8): 992-995.
    [16] 王晶, 姚宜斌, 刘强, 张瑞.  基于GSM蜂窝移动通信系统的CORS新工作模式 . 武汉大学学报 ● 信息科学版, 2009, 34(8): 1000-1003.
    [17] 姜卫平, 马强, 刘鸿飞.  CORS系统中坐标移动转换方法及应用 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 775-778.
    [18] 王勇, 柳林涛, 许厚泽, 黄应生.  利用GPS技术反演中国大陆水汽变化 . 武汉大学学报 ● 信息科学版, 2007, 32(2): 152-155.
    [19] 朱文泉, 陈云浩, 潘耀忠, 李京.  基于GIS和RS的中国植被光利用率估算 . 武汉大学学报 ● 信息科学版, 2004, 29(8): 694-698,714.
    [20] 刘经南, 刘晖.  连续运行卫星定位服务系统——城市空间数据的基础设施 . 武汉大学学报 ● 信息科学版, 2003, 28(3): 259-264.
  • 加载中
图(8) / 表(3)
计量
  • 文章访问数:  826
  • HTML全文浏览量:  92
  • PDF下载量:  347
  • 被引次数: 0
出版历程
  • 收稿日期:  2015-12-21
  • 刊出日期:  2017-09-05

附加辅助层析区域提高射线利用率的水汽反演方法

doi: 10.13203/j.whugis20150592
    基金项目:

    国家自然科学基金 41174012

    国家自然科学基金 41274022

    江西省自然科学基金 20151BAB213031

    作者简介:

    赵庆志, 博士, 主要从事GNSS数据处理及GNSS气象学研究。zhaoqingzhia@163.com

    通讯作者: 姚宜斌, 博士, 教授。ybyao@whu.edu.cn
  • 中图分类号: P228

摘要: 在已有的全球导航卫星系统水汽反演技术中,能够利用的观测值是完整穿过整个层析区域的卫星射线。由于卫星、接收机几何位置分布以及层析区域选择的特定性,许多卫星射线是从层析区域侧面穿出,这些射线通常当作无效信息被剔除,降低了层析结果的精度。针对该缺点,提出并实现了一种附加辅助层析区域提高射线利用率的水汽反演方法,使从研究区域侧面穿出的信号信息也可以被利用。基于美国得克萨斯州(Texas)地区的连续运行参考站(Continuously Operation Reference Stations,CORS)网的实测数据,结合层析区域内的无线电探空仪数据进行实验,结果表明,该方法能够提高水汽反演结果的精度,其均方根误差的改善率为14.6%。

English Abstract

赵庆志, 姚宜斌, 罗亦泳. 附加辅助层析区域提高射线利用率的水汽反演方法[J]. 武汉大学学报 ● 信息科学版, 2017, 42(9): 1203-1208, 1222. doi: 10.13203/j.whugis20150592
引用本文: 赵庆志, 姚宜斌, 罗亦泳. 附加辅助层析区域提高射线利用率的水汽反演方法[J]. 武汉大学学报 ● 信息科学版, 2017, 42(9): 1203-1208, 1222. doi: 10.13203/j.whugis20150592
ZHAO Qingzhi, YAO Yibin, LUO Yiyong. A Method to Improve the Utilization of Observation for Water Vapor Tomography by Adding Assisted Tomographic Area[J]. Geomatics and Information Science of Wuhan University, 2017, 42(9): 1203-1208, 1222. doi: 10.13203/j.whugis20150592
Citation: ZHAO Qingzhi, YAO Yibin, LUO Yiyong. A Method to Improve the Utilization of Observation for Water Vapor Tomography by Adding Assisted Tomographic Area[J]. Geomatics and Information Science of Wuhan University, 2017, 42(9): 1203-1208, 1222. doi: 10.13203/j.whugis20150592
  • 利用全球导航卫星系统(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]讨论了代数重构算法在水汽层析中的各种问题,并对松弛因子和迭代终止条件的确定做了详细研究。

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

    • 在水汽反演过程中,倾斜路径水汽含量(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表示待求网格水汽密度组成的列向量。

    • 由于层析区域内接收机数目有限,导致在给定的层析时段内有射线穿过的网格数目有限,因此,式(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算法求解。

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

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

    • 为了便于表达,将需要研究的层析区域(如图 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, 14],通常将其作为检验层析技术反演水汽结果的标准。本文将均方根误差(root mean square, RMS)、偏差(Bias)和平均绝对偏差(MAE)作为检验本文提出方法反演水汽精度的指标。

    • 选取美国德克萨斯地区运行参考站(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.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
    • 为了验证本文提出方法反演水汽结果的精度,对层析时段内每天两个观测历元(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。这说明,本文提出的利用辅助层析区域提高射线利用率反演水汽的方法其精度要优于传统方法。

    • 为了进一步分析本文提出的方法反演水汽密度在垂直方向上的分布情况,对层析时段内的层析结果进行统计,得到探空站所在位置上不同高度的平均水汽密度值,并与探空站计算的结果进行对比(图 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反演水汽的精度其垂直分布也优于传统方法。这进一步说明,通过本文方法,实现了在研究区域侧面穿出的射线信息对层析区域水汽反演的贡献,提高了水汽反演结果的质量。

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

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

参考文献 (23)

目录

    /

    返回文章
    返回