基于Laplace方程的垂线偏差法反演全球海域重力异常

doi:10.11947/j.AGCS.2020.20190108

测绘学报 ›› 2020, Vol. 49 ›› Issue (4): 452-460.doi: 10.11947/j.AGCS.2020.20190108

基于Laplace方程的垂线偏差法反演全球海域重力异常

张胜军¹, 李建成², 孔祥雪³

1. 东北大学资源与土木工程学院, 辽宁沈阳 110819;
2. 武汉大学测绘学院, 湖北武汉 430079;
3. 辽宁省地震局, 辽宁沈阳 110031

收稿日期:2019-04-05 修回日期:2019-09-28 发布日期:2020-04-17
作者简介:张胜军(1987-),男,博士,讲师,研究方向为卫星测高数据处理及应用研究。E-mail:zhangshengjun@mail.neu.edu.cn
基金资助:
国家自然科学基金（41804002）；中央高校基本科研业务专项资金资助项目（170103009）；武汉大学地球空间环境与大地测量教育部重点实验室开放基金（17-01-04）

Inversion of global marine gravity anomalies with vertical deflection method deduced from Laplace equation

ZHANG Shengjun¹, LI Jiancheng², KONG Xiangxue³

1. School of Resources and Civil Engineering, Northeastern University, Shenyang 110819, China;
2. School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China;
3. Liaoning Earthquake Agency, Shenyang 110031, China

Received:2019-04-05 Revised:2019-09-28 Published:2020-04-17
Supported by:
The National Natural Science Foundation of China (No. 41804002);The Fundamental Research Funds for the Central Universities (No. 170103009);The Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University(No. 17-01-04)

摘要/Abstract

摘要： 海洋重力场模型反演的质量主要依赖于采用测高数据的精度、空间分辨率和数据分布密集程度。本文联合Geosat GM/ERM、ERS-1 GM/ERM、TOPEX/Poseidon、Envisat、Cryosat-2、Jason-1 ERM/GM和SARAL/AltiKa等多种测高观测数据集，深入比较了多种波形重跟踪算法的效果，回波数据重跟踪处理后的沿轨海面高标准差。统计表明，Sandwell算法优于MLE-4算法、Davis阈值法、改进阈值法和β参数拟合法；基于不同测高数据波形重采样的结果给出了沿轨海面梯度计算中低通滤波的参数选择方法，并采用Sandwell提出的垂线偏差法，反演了全球海域1'×1'的重力场模型。检核表明，反演结果与DTU13和SIO V23.1模型检核的差值均方根分别为3.4、1.8 mGal，与NGDC船测数据的检核精度为4~8 mGal，且本文模型在部分典型海区内精度更优。

关键词: 卫星测高, 波形重跟踪, 垂线偏差, 重力异常

Abstract: The quality of marine gravity field derived from satellite altimetry mainly relies on denser data coverage, enhanced spatial resolution as well as improved range precision. In this study, a comparison of application effect for several waveform retracking methods is firstly executed on the basis of multi-satellite altimetry dataset sincluding Geosat GM/ERM,ERS-1 GM/ERM, TOPEX/Poseidon,Envisat,Cryosat-2,Jason-1 ERM/GM and SARAL/AltiKa. According to the statistics of along-track standard deviation, the two-step retracker proposed by Sandwell is superior to other retrackers such as MLE-4, threshold, advanced threshold, and β-5 parameter fitting method. Tuned parameters for certain low-pass filtersare given within the resampling procedure for each altimeter mission in order to calculate accurate along-track slopes. Then a 1'×1' global marine gravity field model is inversed based on the remove-restore procedure by introducing the EGM2008 as reference model. The root mean squares (rms) of difference are respectively 3.4 mGal and 1.8 mGal by comparing with DTU13 and SIO V23.1 model. Meanwhile, the comparison with ship-measured data results a RMS range of 4~8 mGal, while better coincidence shows up in typical ocean areas. In summary, both the comparisons indicate that the new model has a comparable accuracy level with DTU13 and V23.1 at designed 1'×1' grid interval.

Key words: satellite altimetry, waveform retracking, vertical deflection, gravity anomaly

中图分类号:

P228.3

张胜军, 李建成, 孔祥雪. 基于Laplace方程的垂线偏差法反演全球海域重力异常[J]. 测绘学报, 2020, 49(4): 452-460.

ZHANG Shengjun, LI Jiancheng, KONG Xiangxue. Inversion of global marine gravity anomalies with vertical deflection method deduced from Laplace equation[J]. Acta Geodaetica et Cartographica Sinica, 2020, 49(4): 452-460.

参考文献

[1] ANDERSEN O B, KNUDSEN P. Global marine gravity field from the ERS-1 and Geosat geodetic mission altimetry[J]. Journal of Geophysical Research:Oceans, 1998, 103(C4):8129-8137.
[2] ANDERSEN O B, KNUDSEN P, KENYON S, et al. Global and arctic marine gravity field from recent satellite altimetry (DTU13)[C]//Proceedings of the 76th EAGE Conference and Exhibition. Amsterdam:European Association of Geoscientists & Engineers, 2014.
[3] ANDERSEN O B, KNUDSEN P, STENSENG L. The DTU13 MSS (Mean Sea Surface) and MDT (Mean Dynamic Topography) from 20 years of satellite altimetry[C]//2014 International Association of Geodesy Symposia. Cham:Springer, 2015.
[4] SANDWELL D T, SMITH W H F. Marine gravity anomaly from Geosat and ERS1 satellite altimetry[J]. Journal of Geophysical Research:Solid Earth, 1997, 102(B5):10039-10054.
[5] SANDWELL D T, MÜLLER R D, SMITH W H F, et al. New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure[J]. Science, 2014, 346(6205):65-67.
[6] ANDERSEN O B, KNUDSEN P, BERRY P A M. The DNSC08GRA global marine gravity field from double retracked satellite altimetry[J]. Journal of Geodesy, 2010, 84(3):191-199.
[7] SANDWELL D, GARCIA E, SOOFI K, et al. Toward 1-mGal accuracy in global marine gravity from CryoSat-2, Envisat, and Jason-1[J]. The Leading Edge, 2013, 32(8):892-899.
[8] GARCIA E S, SANDWELL D T, SMITH W H F. Retracking CryoSat-2, Envisat and Jason-1 radar altimetry waveforms for improved gravity field recovery[J]. Geophysical Journal International, 2014, 196(3):1402-1422.
[9] 许厚泽, 王海瑛, 陆洋, 等. 利用卫星测高数据推求中国近海及邻域大地水准面起伏和重力异常研究[J]. 地球物理学报, 1999, 42(4):465-471. XU Houze, WANG Haiying, LU Yang, et al. Geoid undulations and gravity anomalies from T/P and ERS-1 altimeter data in the China sea and vicinity[J]. Chinese Journal of Geophysics, 1999, 42(4):465-471.
[10] 李建成, 宁津生, 陈俊勇, 等. 联合Topex/Poseidon, ERS2和Geosat卫星测高资料确定中国近海重力异常[J]. 测绘学报, 2001, 30(3):197-202. DOI:10.3321/j.issn:1001-1595.2001.03.003. LI Jiancheng, NING Jinsheng, CHEN Junyong, et al. Determination of gravity anomalies over the south China sea by combination of TOPEX/Poseidon, ERS2 and Geosat altimeter data[J]. Acta Geodaetica et Cartographica Sinica, 2001, 30(3):197-202. DOI:10.3321/j.issn:1001-1595.2001.03.003.
[11] 黄谟涛, 翟国君, 管铮, 等. 利用卫星测高数据反演海洋重力异常研究[J]. 测绘学报, 2001, 30(2):179-184. DOI:10.3321/j.issn:1001-1595.2001.02.014. HUANG Motao, ZHAI Guojun, GUAN Zheng, et al. On the recovery of gravity anomalies from altimeter data[J]. Acta Geodaetica et Cartographica Sinica, 2001, 30(2):179-184. DOI:10.3321/j.issn:1001-1595.2001.02.014.
[12] 黄谟涛, 翟国君, 欧阳永忠, 等. 利用多代卫星测高数据反演海洋重力场[J]. 测绘科学, 2006, 31(6):37-39. HUANG Motao, ZHAI Guojun, OUYANG Yongzhong, et al. Recovery of marine gravity field using integrated data from multi-satellite missions[J]. Science of Surveying and Mapping, 2006, 31(6):37-39.
[13] SMITH W H F. Resolution of seamount geoid anomalies achieved by the SARAL/AltiKa and envisat RA2 satellite radar altimeters[J]. Marine Geodesy, 2015, 38(S1):644-671.
[14] 张胜军, 李建成, 褚永海, 等. 基于Cryosat和Jason1 GM数据的垂线偏差计算与分析[J]. 武汉大学学报(信息科学版), 2015, 40(8):1012-1017. ZHANG Shengjun, LI Jiancheng, CHU Yonghai, ei al. Calculation and analysis of the deflection of vertical derived from Cryosat and Jason1 GM data[J]. Geomatics and Information Science of Wuhan University, 2015, 40(8):1012-1017.
[15] ZHANG Shengjun, SANDWELL D T. Retracking of SARAL/AltiKa radar altimetry waveforms for optimal gravity field recovery[J]. Marine Geodesy, 2017, 40(1):40-56.
[16] OLGIATI A, BALMINO G, SARRAILH M, et al. Gravity anomalies from satellite altimetry:comparison between computation via geoid heights and via deflections of the vertical[J]. Bulletin Géodésique, 1995, 69(4):252-260.
[17] HAXBY W F, KARNER G D, LABRECQUE J L, et al. Digital images of combined oceanic and continental data sets and their use in tectonic studies[J]. EOS, Transactions American Geophysical Union, 1983, 64(52):995-1004.
[18] SANDWELL D T. Along-track deflection of the vertical from Seasat:GEBCO overlays[R]. NOAA Technical Memorandum NOS NGS-40, Rockville:U.S. Department of Commerce, 1984.
[19] HWANG C. Inverse vening meinesz formula and deflection-geoid formula:applications to the predictions of gravity and geoid over the South China Sea[J]. Journal of Geodesy, 1998, 72(5):304-312.
[20] 李建成, 陈俊勇, 宁津生, 等. 地球重力场逼近理论与中国2000似大地水准面的确定[M]. 武汉:武汉大学出版社, 2003. LI Jiancheng, CHEN Junyong, NING Jinsheng, et al. Approximation Theory of the Earth Gravity Field and Determination of the Chinese Gravity Geoid Model 2000[M]. Wuhan:Wuhan University Press, 2003.
[21] SARRAILH M, BALMINO G, DOUBLET D. The arctic and antarctic oceans gravity field from ERS1 altimetric data[C]//Proceedings of the International Association of Geodesy Symposia. Berlin:Springer, 1997, 117:437-444.
[22] HWANG C, GUO Jinyun, DENG Xiaoli, et al. Coastal gravity anomalies from retracked Geosat/GM altimetry:improvement, limitation and the role of airborne gravity data[J]. Journal of Geodesy, 2006, 80(4):204-216.
[23] GUO Jinyun, GAO Yonggang, HWANG C W, et al. A multi-subwaveform parametric retracker of the radar satellite altimetric waveform and recovery of gravity anomalies over coastal oceans[J]. Science China Earth Sciences, 2010, 53(4):610-616.
[24] 杨元德, 鄂栋臣, 黄金维, 等. Geosat/GM波形重跟踪反演中国沿海区域重力异常[J]. 武汉大学学报(信息科学版), 2008, 33(12):1288-1291. YANG Yuande, E Dongchen, HUANG Jinwei, 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.
[25] SANDWELL D T, SMITH W H F. Retracking ERS-1 altimeter waveforms for optimal gravity field recovery[J]. Geophysical Journal International, 2005. 163(1):79-89.
[26] SANDWELL D T, SMITH W H F. Global marine gravity from retracked Geosat and ERS-1 altimetry:ridge segmentation versus spreading rate[J]. Journal of Geophysical Research:Solid Earth, 2009, 114(B1):B01411.
[27] 常晓涛, 李建成, 郭金运, 等. 一种多前缘多阈值的波形重构算法[J]. 地球物理学报, 2006, 49(6):1629-1634. CHANG Xiaotao, LI Jiancheng, GUO Jinyun, et al. A multi-leading edge and multi-threshold waveform retracker[J]. Chinese Journal of Geophysics, 2006, 49(6):1629-1634.
[28] BAO Lifeng, LU Yang, WANG Yong. Improved retracking algorithm for oceanic altimeter waveforms[J]. Progress in Natural Science, 2009, 19(2):195-203.
[29] 彭福凯, 沈云中. 长江口EnviSat测高数据的波形分类重构分析[J]. 测绘学报, 2015, 44(6):616-624. DOI:10.11947/j.AGCS.2015.20140374. PENG Fukai, SHEN Yunzhong. Analysis of EnviSat altimetric data around Yangtze estuary by waveform retracking[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(6):616-624. DOI:10.11947/j.AGCS.2015.20140374.
[30] SANDWELL D T. Antarctic marine gravity field from high-density satellite altimetry[J]. Geophysical Journal International, 1992, 109(2):437-448.

基于Laplace方程的垂线偏差法反演全球海域重力异常

Inversion of global marine gravity anomalies with vertical deflection method deduced from Laplace equation

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	郭金运, 金鑫, 边少锋, 常晓涛. 垂线偏差测量的固体潮和海潮改正[J]. 测绘学报, 2022, 51(7): 1215-1224.
[2]	徐天河, 穆大鹏, 闫昊明, 郭金运, 尹鹏. 近20年海平面变化成因研究进展及挑战[J]. 测绘学报, 2022, 51(7): 1294-1305.
[3]	魏子卿. 第二大地边值问题引论[J]. 测绘学报, 2022, 51(6): 797-803.
[4]	孙中苗, 管斌, 翟振和, 欧阳明达. 海洋卫星测高及其反演全球海洋重力场和海底地形模型研究进展[J]. 测绘学报, 2022, 51(6): 923-934.
[5]	董杰, 张泽宇, 文汉江, 孙文科. 超导重力数据检测到的2011年日本东北大地震(M_w 9.0)震前重力异常及同震重力变化[J]. 测绘学报, 2022, 51(1): 63-70.
[6]	范雕, 李姗姗, 欧阳永忠, 孟书宇, 陈成, 邢志斌, 张驰. 顾及海底地形非线性项的最小二乘配置反演方法[J]. 测绘学报, 2021, 50(7): 953-971.
[7]	马小辉, 孙中苗, 张志斌, 张阿丽, 袁野, 孙正雄, 王宏. 利用射电天线轴线信息测定VLBI站点垂线偏差[J]. 测绘学报, 2021, 50(3): 315-323.
[8]	宗敬文, 李厚朴, 纪兵, 欧阳永忠. 测高重力异常中央区奇异积分数值求积公式[J]. 测绘学报, 2021, 50(10): 1308-1319.
[9]	胡敏章, 张胜军, 金涛勇, 文汉江, 褚永海, 姜卫平, 李建成. 新一代全球海底地形模型BAT_WHU2020[J]. 测绘学报, 2020, 49(8): 939-954.
[10]	陈俊霖, 周春霞, 赵秋阳. 2003—2018年Byrd冰川流域冰下湖活动及水文联系——多源卫星测高数据监测结果分析[J]. 测绘学报, 2020, 49(5): 547-556.
[11]	范雕, 李姗姗, 杨军军, 孟书宇, 邢志斌, 张驰, 冯进凯. 利用多元回归分析反演西南印度洋区域海底地形[J]. 测绘学报, 2020, 49(2): 147-161.
[12]	赵鸿彬, 谷延超, 范东明, 邱春洪, 宿春鹏, 方伟浩. 联合卫星测高、GRACE与温盐数据分析红海海平面变化季节性特征[J]. 测绘学报, 2019, 48(9): 1119-1128.
[13]	管斌, 孙中苗, 刘晓刚, 翟振和. 验潮站和深海压力仪用于卫星高度计性能评估的可行性分析[J]. 测绘学报, 2019, 48(2): 161-168.
[14]	柯宝贵, 张利明, 章传银, 党亚民. 卫星测高与船载重力测量数据融合的点质量拟合法[J]. 测绘学报, 2018, 47(7): 924-929.
[15]	邢志斌, 李姗姗. 我国陆海统一似大地水准面构建的三维重力矢量法[J]. 测绘学报, 2018, 47(5): 575-583.