留言板

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

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

Cryosat-2数据的大地水准面分辨能力研究

张胜军 金涛勇 褚永海 孔祥雪

张胜军, 金涛勇, 褚永海, 孔祥雪. Cryosat-2数据的大地水准面分辨能力研究[J]. 武汉大学学报 ● 信息科学版, 2016, 41(6): 759-764. doi: 10.13203/j.whugis20140829
引用本文: 张胜军, 金涛勇, 褚永海, 孔祥雪. Cryosat-2数据的大地水准面分辨能力研究[J]. 武汉大学学报 ● 信息科学版, 2016, 41(6): 759-764. doi: 10.13203/j.whugis20140829
ZHANG Shengjun, JIN Taoyong, CHU Yonghai, KONG Xiangxue. Estimation of the Resolution Capability of the Cryosat-2 Altimeter[J]. Geomatics and Information Science of Wuhan University, 2016, 41(6): 759-764. doi: 10.13203/j.whugis20140829
Citation: ZHANG Shengjun, JIN Taoyong, CHU Yonghai, KONG Xiangxue. Estimation of the Resolution Capability of the Cryosat-2 Altimeter[J]. Geomatics and Information Science of Wuhan University, 2016, 41(6): 759-764. doi: 10.13203/j.whugis20140829

Cryosat-2数据的大地水准面分辨能力研究

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

国家重点基础研究发展计划 No.2013CB733301

国家自然科学基金 Nos.41210006,41304003

国家高技术研究发展计划 No.2013AA122502

武汉大学地球空间环境与大地测量教育部重点实验室开放基金 No.15-02-09

详细信息
    作者简介:

    张胜军,博士生,主要从事卫星测高数据处理及应用研究。zhangshengjun@whu.edu.cn

  • 中图分类号: P228

Estimation of the Resolution Capability of the Cryosat-2 Altimeter

Funds: 

The National Basic Research Program (973 Program) of China No.2013CB733301

the National Natural Science Foundation of China Nos.41210006,41304003

the National High Technology Research & Development Program (Program 863) of China No.2013AA122502

the Open Funds of the Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University No.15-02-09

More Information
    Author Bio:

    ZHANG Shengjun, PhD candidate, specializes in satellite altimeter data processing. E-mail: zhangshengjun@whu.edu.cn

  • 摘要: Cryosat-2卫星首次提供了369 d重复周期的密集漂移轨道测高数据,可用于改进海洋重力场反演精度和分辨率。从信号和噪声两个角度对其沿轨数据进行自回归建模,给出了一种以重复轨道观测序列差分加和确定信号谱和噪声谱,并比较分析测高数据大地水准面分辨能力的方法(称为信噪比分析法)。对全球5个试验区域采用谱关联性法和信噪比法来分析Cryosat-2沿轨数据的大地水准面,其分辨能力约为23~33 km,均值约为27 km,优于单周期Geosat/ERM、ERS-1/ERM和T/P数据的38 km、43 km和37 km。研究表明,Cryosat-2数据具有更低噪声水平和更高大地水准面信号辨识度,可以为构建海洋重力场模型提供更高分辨率的短波信息。
  • 图  1  所选Cryosat-2重复轨道数据的分布

    Figure  1.  Distribution of Selected Cryosat-2 Data Series with Repeated Ground Track

    图  2  重复轨道数据序列海面高度对比及差值

    Figure  2.  Comparison and Difference Between Sea Surface Heights and Repeat Track Data Sequences

    图  3  重复轨道数据序列大地水准面分辨能力

    Figure  3.  Analysis of the Resolution Capability of the Repeat Track Data Sequence

    表  1  Cryosat-2数据概要信息

    Table  1.   General Information of Cryosat-2 Data

    周期时间范围轨道号范围
    12010-07~2010-072622~2853
    22010-08~2011-021~2853
    32011-02~2011-121~1795,2578~5344
    42011-12~2012-061~2740
    52012-07~2012-122748~5344
    62012-12~2013-121~5344
    下载: 导出CSV

    表  2  Cryosat-2数据序列的大地水准面分辨能力

    Table  2.   Resolution Capability of the Selected Cryosat-2 Data Sequence

    重复轨道号序列长度谱关联性分析/km信噪比分析/km
    北大西洋 2754A 561 22.241 23.874
    (20°~50°,310°~340°) 5281D 560 23.881 24.634
    北太平洋 3933A 558 27.568 26.385
    (10°~40°,170°~200°) 4055D 559 31.191 29.737
    印度洋 5225A 561 32.151 30.856
    (-50°~-20°, 80°~110°) 4276D 560 26.397 28.537
    南太平洋 2859A 561 31.371 28.259
    (-40°~-10°,200°~230°) 3953D 559 30.256 28.623
    南大西洋 3795A 560 24.504 25.241
    (-50°~-20°,325°~355°) 3745D 538 25.216 26.611
    平均值 -- 27.477 27.276
    下载: 导出CSV

    表  3  不同直径沿轨滤波后Cryosat-2数据序列的大地水准面分辨能力

    Table  3.   Resolution Capability of the Cryosat-2 Data Sequence through Along-Track Filtering with Different Diameter

    重复轨道号 大地水准面分辨能力/km
    172737
    2754A 16.274 14.987 15.152
    5281D 15.645 14.607 14.722
    3933A 17.506 15.694 15.892
    4055D 16.676 15.713 16.023
    5225A 16.986 15.757 15.989
    4276D 16.563 15.062 15.218
    2859A 18.846 17.636 18.548
    3953D 25.137 21.311 22.670
    3795A 16.173 15.224 15.374
    3745D 16.455 15.386 15.568
    下载: 导出CSV
  • [1] Anzenhofer M, Gruber T. Fully Reprocessed ERS-1 Altimeter Data from 1992 to 1995:Feasibility of the Detection of Long Term Sea Level Change[J]. Journal of Geophysical Research,1998, 103(C4):8089-8112
    [2] 李大炜,李建成,金涛勇,等.利用多代卫星测高资料监测1993~2011年全球海平面变化[J].武汉大学学报·信息科学版,2012,37(12):1421-1424

    Li Dawei,Li Jiancheng, Jin Taoyong, et al. Monitoring Global Sea Level Change from 1993 to 2011 Using TOPEX and Jason Altimeter Missions[J].Geomatics and Information Science of Wuhan University, 2012,37(12):1421-1424
    [3] Marsh J G, Koblinsky J, Zwally H J, et al. A Global Mean Sea Surface based Upon GEOS3 and Seasat Altimeter Data[J]. Journal of Geophysical Research, 1992, 97(B4):4915-4921
    [4] 金涛勇,李建成,姜卫平,等.基于多源卫星测高数据的新一代全球平均海面高模型[J].测绘学报,2011,40(6):723-729

    Jin Taoyong, Li Jiancheng, Jiang Weipin, et al. The New Generation of Global Mean Sea Surface Height Model Based on Multialtimetric Data[J]. Acta Geodaetica et Cartographica Sinica, 2011,40(6):723-729
    [5] 章传银,李建成,晁定波,等.联合卫星测高和海洋物理数据计算近海稳态海面地形[J].武汉测绘科技大学学报,2000,25(6):500-504

    Zhang Chuanyin, Li Jiancheng, Chao Dingbo, et al. Calculating Stationary Sea Surface Topography of Coastal Area with Altimeter and Oceanographic Data[J]. Journal of Wuhan Technical University of Surveying and Mapping, 2000,25(6):500-504
    [6] 邓凯亮,暴景阳,章传银,等.联合多代卫星测高数据确定中国近海稳态海面地形模型[J].测绘学报,2009,38(2):114-119

    Deng Kailiang, Bao Jingyang, Zhang Chuanyin, et al. The Determination of Quasi-Stationary Sea Surface Topography over China Sea by Using Multi-altimeter Data[J]. Acta Geodaetica et CartographicaSinica, 2009,38(2):114-119
    [7] Laxon S W, Giles K A, Ridout A L, et al. CryoSat-2 Estimates of Arctic Sea Ice Thickness and Volume[J].Geophysical Research Letters,2013, 40:732-737
    [8] 王虎彪,王勇,陆洋.联合多种测高数据确定中国边缘海及全球海域的垂线偏差[J].武汉大学学报·信息科学版,2007,32(9):770-773

    Wang Hubiao, Wang Yong, Lu Yang. High Precision Vertical Deflection over China Marginal Sea and Global Sea Derived from Multi-satellite Altimeter[J]. Geomatics and Information Science of Wuhan University, 2007, 32(9):770-773
    [9] Sandwell D, Smith W H F. Marine Gravity Anomaly from Geosat and ERS-1 Satellite Altimetry[J]. Journal of Geophysical Research, 1997, 102(B5):10039-10054
    [10] Andersen O B, Knudsen P. Global Marine Gravity Field from the ERS-1 and Geosat Geodetic Mission Altimetry[J]. Journal of Geophysical Research, 1998, 103(C4):8129-8137
    [11] 李建成,宁津生,陈俊勇,等.中国海域大地水准面和重力异常的确定[J]. 测绘学报,2003,32(2):114-119

    Li Jiancheng, Ning Jinsheng, Chen Junyong, et al. Geoid Determination in China Sea Areas[J]. Acta Geodaetica et Cartographica Sinica, 2003,32(2):114-119
    [12] 杨元德,鄂栋臣,黄金维,等.Geosat/GM波形重跟踪反演中国沿海区域重力异常[J].武汉大学学报·信息科学版,2008,33(12):1288-1291

    Yang Yuande, E Dongchen, Hwang Cheinway, et al. Chinese Coastal Gravity Anomalies from Waveform Retracked Geosat/GM Altimetry[J]. Geomatics and Information Science of Wuhan University, 2008,33(12):1288-1291
    [13] Yale M M, Sandwell D, Smith W H F. Comparison of Along-Track Resolution of Stacked Geosat, ERS-1 and Topex Satellite Altimeters[J]. Journal of Geophysical Research, 1995, 100(B8):15117-15127
    [14] Wingham D J, Francis C R, Baker S, et al. CryoSat:A Mission to Determine the Fluctuations in Earth's Land and Marine Ice Fields[J]. Adcances in Space Research,2006, 37:841-871
    [15] Brammer R F, Sailor R V. Preliminary Estimates of the Resolution Capability of the Seasat Radar Altimeter[J].Geophysical Research Letters, 1998, 7(3):193-196
    [16] Sailor R V, Driscoll M L. Comparison of Noise Models and Resolution Capabilities for Satellite Radar Altimeters[C]. Oceans 92 Mastering the Oceans Through Technology, Newport, 1992
    [17] Marks K M, Sailor R V. Comparison of GEOS-3 and Seasat Altimeter Resolution Capabilities[J]. Geophysical Research Letters, 1986, 13(7):697-700
    [18] Maus S, Green C M, Fairhead J D. Improved Ocean-Geoid Resolution from Retracked ERS-1 Satellite Altimeter Waveforms[J]. Geophysical Journal International, 1998, 134(1):243-253
    [19] 张玲华,郑宝玉.随机信号处理[M].北京:清华大学出版社,2003

    Zhang Linghua, Zheng Baoyu. Stochastic Signal Processing[M].Beijing:Press of Tsinghua University, 2003
    [20] Sailor R V. Signal Processing Techniques, Geoid and Its Geophysical Interpretation[M]. Boca Raton, FL:CRC Press, 1993
    [21] Labroue S, Boy F, Picot N, et al. First Quality Assessment of the Cryosat-2 Altimeter System over Ocean[J]. Advances in Space Research, 2012, 50:1030-1045
    [22] ESRIN, ESA and Mullard Space Science Laborator. Cryosat Product Handbook[EB/OL]. University College London, https://earth.esa.int/documents/10174/125272/CryoSat_Product_Handbook, 2012
    [23] Driscoll A L, Sailor R V. GFO On-Orbit Altimeter Noise Assessment[R]. NASA:NTRS, 2001
  • [1] 柯灝, 赵建虎, 周丰年, 吴敬文, 暴景阳, 赵祥伟, 谢朋朋.  联合大地水准面、海面地形和潮波运动数值模拟的长江口陆海垂直基准转换关系 . 武汉大学学报 ● 信息科学版, 2022, 47(5): 731-737. doi: 10.13203/j.whugis20200601
    [2] 李豪, 顾勇为, 韩松辉.  基于信噪比检验的双截断奇异值估计 . 武汉大学学报 ● 信息科学版, 2019, 44(2): 228-232, 239. doi: 10.13203/j.whugis20170051
    [3] 程鹏飞, 文汉江, 刘焕玲, 董杰.  卫星大地测量学的研究现状及发展趋势 . 武汉大学学报 ● 信息科学版, 2019, 44(1): 48-54. doi: 10.13203/j.whugis20180356
    [4] 郑南山, 丰秋林, 刘晨, 周晓敏.  GPS反射信号信噪比与NDVI相关性研究 . 武汉大学学报 ● 信息科学版, 2019, 44(10): 1423-1429. doi: 10.13203/j.whugis20180046
    [5] 黄雄华, 王宏霞, 蒋伟贞, 崔更申.  利用抖动量化和信噪比的自适应鲁棒音频盲水印算法 . 武汉大学学报 ● 信息科学版, 2013, 38(1): 121-125.
    [6] 游为, 范东明, 贺全兵.  利用GOCE卫星轨道反演地球重力场模型 . 武汉大学学报 ● 信息科学版, 2012, 37(3): 294-297.
    [7] 程芦颖, 许厚泽.  顾及测站点上重力场信息的大地水准面高的拟合方法 . 武汉大学学报 ● 信息科学版, 2008, 33(7): 701-705.
    [8] 吴雨航, 陈秀万, 吴才聪.  利用信噪比削弱多路径误差的方法研究 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 842-845.
    [9] 方剑, 马宗晋.  地球及其形状非对称性的重力学研究 . 武汉大学学报 ● 信息科学版, 2005, 30(8): 748-751.
    [10] 楼立志, 许厚泽.  模拟大地水准面应用于GPS水准内插的研究 . 武汉大学学报 ● 信息科学版, 2005, 30(1): 47-50.
    [11] 邵永社, 陈鹰, 李晶.  一种基于小波系数相关性的图像降噪方法 . 武汉大学学报 ● 信息科学版, 2005, 30(9): 771-774.
    [12] 申文斌, 宁津生, 李建成, 晁定波.  论相对论重力位及相对论大地水准面 . 武汉大学学报 ● 信息科学版, 2004, 29(10): 897-900.
    [13] 王正涛, 李建成, 晁定波, 姜卫平.  利用卫星测高数据研究海面高月异常变化与厄尔尼诺现象的相关性 . 武汉大学学报 ● 信息科学版, 2004, 29(8): 699-703.
    [14] 张全德, 张燕平, 张鹏, 陈现军.  精化区域大地水准面若干问题的探讨 . 武汉大学学报 ● 信息科学版, 2003, 28(S1): 94-96.
    [15] 申文斌, 宁津生, 李建成, 晁定波.  论大地水准面 . 武汉大学学报 ● 信息科学版, 2003, 28(6): 683-687.
    [16] 罗志才, 宁津生, 杨沾吉, 陈永奇.  高分辨率厘米级局部大地水准面的典型应用 . 武汉大学学报 ● 信息科学版, 2003, 28(S1): 100-103.
    [17] 李建成, 姜卫平, 章磊.  联合多种测高数据建立高分辨率中国海平均海面高模型 . 武汉大学学报 ● 信息科学版, 2001, 26(1): 40-45.
    [18] 陈俊勇, 李建成, 宁津生, 晁定波, 张燕平, 张骥.  中国新一代高精度、高分辨率大地水准面的研究和实施 . 武汉大学学报 ● 信息科学版, 2001, 26(4): 283-289,302.
    [19] 杨沾吉, 於宗俦, 于正林, 陈永奇.  GPS水准、天文重力水准与重力大地水准面多种数据联合处理的研究 . 武汉大学学报 ● 信息科学版, 1996, 21(4): 330-337.
    [20] 李建成, 宁津生, 晁定波.  卫星测高在物理大地测量应用中的若干问题 . 武汉大学学报 ● 信息科学版, 1996, 21(1): 9-14.
  • 加载中
图(3) / 表(3)
计量
  • 文章访问数:  1086
  • HTML全文浏览量:  30
  • PDF下载量:  317
  • 被引次数: 0
出版历程
  • 收稿日期:  2015-03-27
  • 刊出日期:  2016-06-05

Cryosat-2数据的大地水准面分辨能力研究

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

    国家重点基础研究发展计划 No.2013CB733301

    国家自然科学基金 Nos.41210006,41304003

    国家高技术研究发展计划 No.2013AA122502

    武汉大学地球空间环境与大地测量教育部重点实验室开放基金 No.15-02-09

    作者简介:

    张胜军,博士生,主要从事卫星测高数据处理及应用研究。zhangshengjun@whu.edu.cn

  • 中图分类号: P228

摘要: Cryosat-2卫星首次提供了369 d重复周期的密集漂移轨道测高数据,可用于改进海洋重力场反演精度和分辨率。从信号和噪声两个角度对其沿轨数据进行自回归建模,给出了一种以重复轨道观测序列差分加和确定信号谱和噪声谱,并比较分析测高数据大地水准面分辨能力的方法(称为信噪比分析法)。对全球5个试验区域采用谱关联性法和信噪比法来分析Cryosat-2沿轨数据的大地水准面,其分辨能力约为23~33 km,均值约为27 km,优于单周期Geosat/ERM、ERS-1/ERM和T/P数据的38 km、43 km和37 km。研究表明,Cryosat-2数据具有更低噪声水平和更高大地水准面信号辨识度,可以为构建海洋重力场模型提供更高分辨率的短波信息。

English Abstract

张胜军, 金涛勇, 褚永海, 孔祥雪. Cryosat-2数据的大地水准面分辨能力研究[J]. 武汉大学学报 ● 信息科学版, 2016, 41(6): 759-764. doi: 10.13203/j.whugis20140829
引用本文: 张胜军, 金涛勇, 褚永海, 孔祥雪. Cryosat-2数据的大地水准面分辨能力研究[J]. 武汉大学学报 ● 信息科学版, 2016, 41(6): 759-764. doi: 10.13203/j.whugis20140829
ZHANG Shengjun, JIN Taoyong, CHU Yonghai, KONG Xiangxue. Estimation of the Resolution Capability of the Cryosat-2 Altimeter[J]. Geomatics and Information Science of Wuhan University, 2016, 41(6): 759-764. doi: 10.13203/j.whugis20140829
Citation: ZHANG Shengjun, JIN Taoyong, CHU Yonghai, KONG Xiangxue. Estimation of the Resolution Capability of the Cryosat-2 Altimeter[J]. Geomatics and Information Science of Wuhan University, 2016, 41(6): 759-764. doi: 10.13203/j.whugis20140829
  • 卫星测高技术广泛应用于海平面变化监测、平均海平面模型的建立[1-4],动力海面地形确定、海冰研究及垂线偏差解算等科研领域[5-8],也是反演全球和局部海域的大地水准面与重力异常的主要手段[9-12]。随着测量模式、测量精度以及数据覆盖率的不断完善,如何合理综合利用多颗测高卫星的海量观测数据获取更高空间分辨率和精度的海洋重力场模型成为国际上研究的热点。而高分辨率海洋重力场模型构建的关键在于测高数据的轨道间距和大地水准面分辨能力[13]

    卫星测高数据分为两种:一种是精密重复轨道观测数据,如T/P、Envisat等,重复周期为10~35 d,对应赤道处轨道间距约为80~315 km,不能满足高分辨率海洋重力场反演要求;二是大地测量漂移轨道观测数据,如Geosat/GM(1985~1986,176 d)、ERS-1/GM(1994~1995,168 d)和Jason-1/GM(2012~2013,406 d)等,这种任务观测周期长,难以实现重复周期观测数据,在赤道处对应的轨道间距约为7~10 km,是用于海洋重力反演的主要数据源。Cryosat-2重复周期为369 d,迄今已经发生三次重复,其赤道处的轨道间距约为7.5 km[14],是有望改进海洋重力场反演的重要数据。

    大地水准面分辨能力对应测高数据可分辨的大地水准面短波特征,决定了海洋重力场的分辨率[15],对于沿轨滤波处理时的参数选择具有重要的参考价值[16]。Brammer等人提出谱关联性分析法[15],根据重复轨道数据之间存在空间频率对应的谱关联性分析大地水准面分辨能力。此后,国外研究学者使用该方法分析了多颗测高卫星数据的大地水准面分辨能力,如Geos-3和Seasat分别为75 km和50 km[17],Geosat/ERM、ERS-1/ERM以及T/P约为38 km、43 km和37 km,波形重构之后ERS-1/ERM数据可提升至35 km左右[18],共线处理之后的Geosat/ERM、ERS-1/ERM、T/P数据可提升至24 km、30 km和22 km[13]

    本文分别采用谱关联分析法和信噪比分析法推算了Cryosat-2数据的大地水准面分辨能力,为应用Cryosat-2数据改进现有海洋重力场模型提供相关研究参考。

    • 根据卫星测高原理,卫星沿星下点轨迹测量雷达脉冲反射足印区内的平均海面高度,回波足印区的直径范围约为3~5 km。不同周期的重复轨道间距通常在1 km之内,相对于足印区范围而言,重复轨道的星下点轨迹可看作完全重复,即包含相同的大地水准面信号。因此,作为大地水准面波长分辨能力分析的前提条件,重复轨道实际是要求不同数据序列包含相同的信号部分[16]

      假定z1和z2分别表示一对重复轨道的沿轨海表面高度数据序列;u表示沿轨的大地水准面信号;n1n2分别表示沿轨测量数据的噪声部分;s表示沿轨距离;存在如下关系式:

      (1)
      (2)
    • 为了描述数据序列之间在频域的谱关联性,引入均方一致性,定义为不同序列的互相关谱的平方除以单独数据序列的自相关谱的乘积,取值范围在0~1之间,1表示数据序列完全相关,0表示完全不相关,其数学表达式为:

      (3)

      式中,Sz1Sz2Sz1z2分别表示数据序列的自相关谱和互相关谱,可以采用周期图法、自相关法或参数估计法估算,其中周期图法和自相关法属于古典谱估计方法,使用过程中无法克服谱分辨率低的问题[19]。Brammer等指出针对沿轨测高数据的采样特征,自回归(auto-regressive,AR)建模可以获取令人满意的谱估计结果[15, 20]。数据序列的AR模型表达式为:

      (4)

      式中,n表示AR模型的阶数;aiwk是AR模型的系数和噪声相关参数。Marks等定义了判定大地水准面波长分辨能力的标准,即重复轨道数据的均方一致性下降为0.5时所对应的空间波长[16]。Brammer 等曾经根据均方一致性下降为0作为分析Seasat数据的判断标准[15],但0.5是更加保守和通用的选择[13, 15, 16]

    • 假定信号和噪声部分不相关,重复轨道时间序列之间作差可以消去信号部分,若由f表示空间采样频率,可以得到差分序列的功率谱和噪声谱之间的关系;顾及信号振幅与功率谱之间存在平方的关系,联合重复轨道数据加和以及差分序列的功率谱可以估算信号的功率谱。具体的关系可以表达为:

      (5)
      (6)

      式中,Sz1-z2Sz1+z2分别表示重复轨道数据差分以及加和序列的功率谱;SnSu分别表示重复轨道数据序列的噪声谱和大地水准面信号谱。在谱分析的过程中,信号与噪声的功率谱通过一对数据序列差分加和及AR建模估算得到。本文从信噪比角度考虑,当信号谱与噪声谱的比值下降为1时,大地水准面信号的功率谱与噪声谱相交处的空间频率可判定为对应大地水准面信号的分辨率,简称为信噪比分析方法,是一种更加直观判定信号分辨能力的准则。

    • Cryosat-2卫星于2010年4月8日发射,在轨运行四年多来,首次提供了重复观测的密集轨道间距测高数据[21]。欧空局发布的基线B版本的L2级地球物理数据记录中,提供了沿轨20 Hz的时间、经纬度以及顾及各校正项的表面高信息[22]。为了控制沿轨数据质量,相关数据根据回波标识、平均海平面模型检核以及三倍中误差准则判定对沿轨20 Hz数据进行剔除,针对满足筛选条件的数据,平均得到沿轨1 Hz的表面高数据。

    • 本文按照二级数据头文件中的周期号CYCLE和轨道相对编号REL_ORBIT对数据进行命名,所使用数据的时间和轨道号信息如表 1所示,需要注意的是各个周期的轨道数目并不完全一致。Cryosat-2数据手册指出2011年10月之前的B版本数据产品中包含不准确的周期编号和轨道参数[22],重复轨道数据序列不能仅依据轨道编号判定,需要辅助以判定轨道间距小于1 km的条件筛选得到。

      表 1  Cryosat-2数据概要信息

      Table 1.  General Information of Cryosat-2 Data

      周期时间范围轨道号范围
      12010-07~2010-072622~2853
      22010-08~2011-021~2853
      32011-02~2011-121~1795,2578~5344
      42011-12~2012-061~2740
      52012-07~2012-122748~5344
      62012-12~2013-121~5344

      沿轨数据的频谱分析要求采样间隔一致,因此数据处理过程中需判定沿轨采样点是否缺失,若存在缺失点则采用相邻两个采样点数据插值补齐采样间隔。为了抑制粗差数据对频谱分析的影响,对差分数据序列进行探测,若超过三倍标准差则移去该差分值并采用线性插值填补数据空白[23]

    • 在北大西洋、北太平洋、印度洋、南太平洋、南大西洋海域内选择5个30°×30°的区域作为研究区域,每个区域内分别挑选一组升弧和降弧的重复轨道数据序列对,筛选方法如§2.1所述,轨道编号及地理位置如图 1所示。以北大西洋海域内轨道编号为2754的升弧重复轨道数据序列对(Revs2754A)为例详细给出大地水准面波长分辨能力的分析过程。

      图  1  所选Cryosat-2重复轨道数据的分布

      Figure 1.  Distribution of Selected Cryosat-2 Data Series with Repeated Ground Track

      图 2(a)给出了重复轨道数据序列对Revs2754A对应的两组海面高比较,可以看出符合较好,反映出一致的大地水准面变化特征。图 2(b)是重复轨道海面高差分数据序列,相对大地水准面信号而言统称为噪声部分,主要由海洋现象、径向轨道误差以及高度计测量噪声引起。

      图  2  重复轨道数据序列海面高度对比及差值

      Figure 2.  Comparison and Difference Between Sea Surface Heights and Repeat Track Data Sequences

      基于时间序列数据分析的R语言软件,分别将沿轨海面高及其差分求和数据序列作为目标数据,首先通过增广迪基-福勒检验判定目标为稳态可建模数据序列;然后,通过绘制目标数据序列的自相关函数图、偏自相关函数图以及广义自相关函数图判定目标过程是否应该采用AR模型建模;最后结合赤池信息量准则(Akaike information criterion,AIC),选择最优模型阶数。不同数据序列的建模阶数不同,例如示例弧段Revs2754A的数据序列进行AR建模的最优阶数为7阶,Revs5281D数据序列AR建模的最优阶数为8阶。

      图 3(a)给出了AR建模之后沿轨海面高(ssh1/2)、大地水准面信号(geoid)以及噪声(noise)的功率谱图。依据信噪比方法判定信噪比下降为1的准则,即确定相应的信号谱下降至与噪声谱相交处,确定出对应的大地水准面分辨能力约为23.874 km。根据谱关联性分析方法计算重复轨道数据序列之间的均方一致性,如图 3(b)所示,按下降为0.5的判定标准确定的大地水准面分辨能力为22.241 km。

      图  3  重复轨道数据序列大地水准面分辨能力

      Figure 3.  Analysis of the Resolution Capability of the Repeat Track Data Sequence

      针对5个研究区域内的其他重复轨道数据序列对,采用与Revs2754A相同的处理方法,分析得到相应的信号波长分辨能力,结果见表 2

      表 2  Cryosat-2数据序列的大地水准面分辨能力

      Table 2.  Resolution Capability of the Selected Cryosat-2 Data Sequence

      重复轨道号序列长度谱关联性分析/km信噪比分析/km
      北大西洋 2754A 561 22.241 23.874
      (20°~50°,310°~340°) 5281D 560 23.881 24.634
      北太平洋 3933A 558 27.568 26.385
      (10°~40°,170°~200°) 4055D 559 31.191 29.737
      印度洋 5225A 561 32.151 30.856
      (-50°~-20°, 80°~110°) 4276D 560 26.397 28.537
      南太平洋 2859A 561 31.371 28.259
      (-40°~-10°,200°~230°) 3953D 559 30.256 28.623
      南大西洋 3795A 560 24.504 25.241
      (-50°~-20°,325°~355°) 3745D 538 25.216 26.611
      平均值 -- 27.477 27.276

      表 2中可以看出相同数据序列的分辨能力也存在一定的差异,因为信噪比方法针对数据序列进行AR建模和参数谱估计,依赖于建模的优劣;而谱关联性方法直接由数据序列的互相关谱和自相关谱计算,谱估计时使用一定宽度的窗函数,本文选择300~400 km的汉宁窗[13]。不同研究区域内Cryosat-2数据的大地水准面分辨能力变化范围在23~33 km之间,平均值约为27 km,优于单个周期的Geosat/ERM、ERS-1/ERM 和T/P观测数据[13]。尽管结果略逊于共线处理后的T/P和Geosat数据,但其地面轨迹间距显著优于共线数据。因此,Cryosat-2卫星可以为构建高分辨率海洋重力场模型提供密集轨道间距数据和更优的短波信息。此外,大西洋海域范围内的研究区域分析结果显著优于其他三个区域,这是因为受大西洋中脊的影响,研究区域内的大地水准面短波信号比较显著。

    • 测高数据大地水准面分辨能力为沿轨滤波处理提供滤波半径的选取依据[15],Cryosat-2数据在不同海域内的平均值约为27 km,因此本文分别选择17 km、27 km和37 km的滤波直径对重复轨道数据序列分别进行处理,高斯滤波后分析的信号分辨能力见表 3

      表 3  不同直径沿轨滤波后Cryosat-2数据序列的大地水准面分辨能力

      Table 3.  Resolution Capability of the Cryosat-2 Data Sequence through Along-Track Filtering with Different Diameter

      重复轨道号 大地水准面分辨能力/km
      172737
      2754A 16.274 14.987 15.152
      5281D 15.645 14.607 14.722
      3933A 17.506 15.694 15.892
      4055D 16.676 15.713 16.023
      5225A 16.986 15.757 15.989
      4276D 16.563 15.062 15.218
      2859A 18.846 17.636 18.548
      3953D 25.137 21.311 22.670
      3795A 16.173 15.224 15.374
      3745D 16.455 15.386 15.568

      表 3可以看出,滤波半径过小导致数据序列中的高频噪声不能有效抑制,滤波半径过大会造成大地水准面信号因过度滤波而损失,适中的滤波半径可以最大限度地提升数据序列的大地水准面分辨能力。其中南太平洋海域内海底相对平坦,而测高数据中包含相对较高的海洋变化噪声,因此大地水准面分辨能力在该海域显著降低。

    • 大地水准面分辨能力不同于测高数据的采样率,提供了沿轨滤波处理时直径选取的依据。本文联合信噪比法和谱关联性法,选取5个研究区域,分析了Cryosat-2数据的大地水准面波长分辨能力,其变化范围在23~33 km之间。

      目前,Cryosat-2数据的重复周期较少,随着在轨运行时间的累积,多个周期重复轨道数据通过共线处理和波形重构,可以有效抑制海面时变效应及地理环境干扰所致的噪声,从而进一步提升信号分辨能力。Cryosat-2数据已在全球大部分海域提供了密集的数据覆盖,因其大地水准面短波信息的更优表达特性,将是改进海洋重力场模型构建的重要数据源。

参考文献 (23)

目录

    /

    返回文章
    返回