高精度测高重力场反演南海海底地形

随着测高技术的不断发展,测高海洋重力场的精度不断提高,由其反演计算的海底地形精度也相应提高。利用Sandwell 2014年最新发布主要由Jason-1末期大地测量任务和Cryosat-2观测资料反演的测高重力场V23.1,采用重力地质法(GGM)反演了中国南海海域海底地形模型,并对结果进行了精度评价。该过程中首先直接计算一系列密度差下的反演结果,并通过船测数据检核选定优化的密度差范围,然后利用向下延拓的方法确定了最优密度差异常数为7g/cm3。与船测水深数据相比,反演得到的GGM模型与检核点船测水深数据差值的标准差达到了±70.32m;此外,还计算了由测高重力异常V15.2反演的海底模型,比较这两个模型发现:测高重力场短波部分的改善对海底地形反演精度的提升作用有限,为得到更高精度的海底地形模型,需引入短波更为敏感的重力梯度资料。     

基于频域法的留尼汪海域海底地形反演

海底地形是全球地形的重要组成部分,对地球物理科学研究、经济活动等具有重要作用。基于Parker公式,利用卫星测高重力异常和船测水深数据,采用频域的方法反演了疑似马航MH370失事区域的留尼汪海域的10°×10°的海底地形。最后将反演的水深和船测水深、国际通用的海深模型ETOPO1作比较进行精度评估,结果表明:本文反演结果与船测水深相比误差平均值为-26.038 m,标准差为176.588 m;与ETOPO1相比,差异平均值为-33.541 m,标准差为160.769 m。这表明采用重力异常数据,结合船测数据能较高精度地反演海底地形。     

测高重力场反演海底地形方法比较

根据重力异常与海底地形的相关性,利用日本西之岛附近海域内的重力异常和船测水深,分别采用SandwellSmith方法(SS,1994)和重力地质法(GGM)反演了该海域内1'×1'海底地形模型,并将反演的模型与格网化船测水深模型和ETOPO1模型进行了对比分析。在GGM中,使用了迭代法向下延拓确定密度差常数为5.7g/cm~3。与SS模型相比,GGM模型与格网化船测水深模型的差值的标准差达到了37.10m。     

基于重力地质法反演南海海底地形

基于重力地质法(GGM)采用卫星重力异常数据反演了中国南海112°E~119°E,16°N~20°N范围1'×1'的海底地形模型,通过船测水深数据的直接格网化结果和ETOP01模型的检核进行了精度评价。结果表明,GGM模型较船测水深的直接格网化结果更为精细,在船测控制点分布均匀的海域,GGM模型能获得比ETOP01模型更高的精度。为进一步验证GGM模型受海底地形特征影响的情况,采用GGM模型与ETOP01模型分别内插出9356个检核点处的水深,并统计得到两模型差值的较差结果平均值为19.560m,标准差为130.156m,相对精度为4.79%。经分析得出:GGM模型与ETOP01模型在坡度变化平缓深海区有着较好的一致性,随着水深的依次变浅,二者的差异逐渐显现,在多海山的海域差异达到最大。最后,为充分发挥重力地质法在海底地形反演中的优势总结了其最佳适用条件。     

应用重力地质法反演马里亚纳海沟地形

根据重力地质法(Gravity Geological Method, GGM),利用6 736个已知水深点数据和卫星测高重力异常反 演了马里亚纳海沟(142^°36^''~147^°18^''E,23^°~27^°N)1^''×1^''分辨率的海底地形,该计算中的密度差常数为1.20 g/cm³。将反演的地形与实测检核水深数据相比,两者较差标准差为152.9 m,平均误差为(±3.0) m,均方差为153.0 m,优于ETOPO1模型和直接把船测水深控制点格网化后的模型。将重力地质法反演的模型、ETOPO1模型、直接格网模型进行功率谱密度分析,结果显示3种模型在中长波地形能量一致,但对于短波地形,重力地质法模型能量高于其他2个模型,说明重力地质法更能描绘细致的海底地形地貌。选取2条测线与重力地质法模型进行对比研究,结果表明重力地质法在海底地形起伏小的区域优于在海底地形起伏较大区域的反演效果。     

卫星测高在海洋测绘中的应用

研究了卫星测高数据在海山探测以及无图海域水深预测方面的应用。介绍了匹配滤波法检测海山的理论与方法，并研究了反演海底地形时的重力导纳理论以及滤波窗口和滤波函数的选择。实践表明，利用卫星测高数据反演海底地形的精度可优于水深的10％。     

卫星测高资料在反演海底地形中的应用

首先论述了利用卫星测高数据反演海底地形两类常用算法即解析算法和统计算法的基本原理和数学模型，在此基础上基于最小二乘配置理论，提出了统计算法的改进模型，使用新模型在南中国海地区进行了海底地形反演计算，并将反演结果与实际船测水深进行比对，进一步验证改进模型的可靠性和有效性。     

联合重力异常和重力垂直梯度异常数据反演皇帝山海域海底地形

海底地形对开展海洋科学调查和研究十分重要。以多波束为主的回声测深技术测量成本高且效率低,几十年来仅实现了全球约20%的海床测绘。对于空白区（特别是深海区域）,可以借助重力异常和重力垂直梯度异常进行回归分析反演得到,但该方法得到的比例因子鲁棒性不强。为了解决这一问题,同时考虑到两种重力数据在表征海底地形长短波长的不同优势,本文结合滑动窗口赋权和稳健回归分析来反演海底地形。在太平洋皇帝山海域（35°～45°N,165°～175°E）的实验结果表明:在船测检核点处,本文构建模型的标准差为61.02 m,相比于单一重力数据反演模型,精度分别提高了14.92%（重力异常）和2.08%（重力垂直梯度异常）,能较好地反映皇帝海山链的地形走势。     

西太平洋海域卫星测高重力异常反演与精度评估

联合T/P、ERS1/GM、ERS1/ERM、ERS2/ERM、GEOSAT/ERM、GEOSAT/GM等多源卫星测高数据,基于逆Vening-Meinesz公式和EGM2008模型,采用移去-恢复方法和快速傅里叶变换算法构建了中国西太平洋海域(0°~40°N,105°E~145°E)1'×1'重力异常模型,选取两个不同特征区域航线,将构建的重力异常模型、EGM2008重力异常模型、美国Scripps海洋研究所重力异常与船测重力数据进行比较分析。结果表明,构建的重力异常模型与船测重力总体趋势变化较为一致、与船测重力比较均方根差为4.16m Gal,总体精度与EGM2008模型重力异常、美国Scripps海洋研究所重力异常相当,反演精度和分辨率达到国际领先水平。     

重力数据融合与重力垂直梯度异常反演

重力垂直梯度异常反应了重力异常的空间变化率,在地球物理勘探等多学科中得到越来越多的应用。利用南海局部区域实测重力异常数据和Sandwell测高重力异常数据,将搜索范围、距离和精度多种因素融合考虑并对Shepard算法进行改进,给出了南海局部区域(19°N~20.5°N,114°E~115.5°E)分辨率1'×1'的重力异常并反演了对应分辨率的重力垂直梯度异常。结果表明,基于Shepard改进算法的高精度船测重力和测高重力的有机融合,增强了单一测高重力数据反演重力垂直梯度异常的细节纹理,提高了反演重力垂直梯度异常的分辨率和精度。     

Bathymetry enhancement by altimetry-derived gravity anomalies in the East Sea (Sea of Japan)

The gravity-geologic method (GGM) was used to enhance the bathymetry of the East Sea (Sea of Japan) with satellite altimetry-derived free-air gravity anomalies and shipborne depth measurements. By comparison with the bathymetry model of Smith and Sandwell’s (SAS) approach (1994), GGM was found to have an advantage with short wavelength (≤12 km) components, while SAS better predicts longer wavelength (≥25 km) components, despite its dependency on density contrast. To mitigate this limitation, a tuning density contrast of 10.25 g/cm³ between seawater and the seafloor was primarily estimated by the downward continuation method and then validated by the check points method with GGM. Similarly, SAS is limited by the “A” value in low-pass part of the Wiener filter, which defines the effective range of the wavelength components on bathymetry. As a final result, we present an enhanced GGM bathymetry model by integrating all available data.     

Improvement in the Determination of the Marine Geoid by Estimating the Bathymetry from Altimetry and Depth Soundings

Abstract

The contribution of bathymetry to the estimation of gravity field related quantities is investigated in an extended test area in the Mediterranean Sea. The region is located southwest of the island of Crete, Greece, bounded between 33^? ≤ ? ≤ 35^? and 15^? ≤ λ ≤ 25^?. Gravity anomalies from the KMS99 gravity field and shipborne depth soundings are used with a priori statistical characteristics of depths in a least-squares collocation procedure to estimate a new bathymetry model. Two different global bathymetry models, namely JGP95E and Sandwell and Smith V8, are used to derive the depth a priori statistical information, while the estimated model is compared against both the global ones and the shipborne depth soundings to assess whether there is an improvement. Various marine geoid models are estimated using ERS1 and GEOSAT Geodetic Mission altimetry and shipborne gravity data. In that process, the effect of the bathymetry is computed using both the estimated and the original depths through a residual terrain modeling reduction. The TOPEX/Poseidon Sea Surface Heights, known for their high accuracy and precision, and the GEOMED solution for the geoid in the Mediterranean are used as control for the validation of the new geoid models and to assess the improvement that the estimated depths offer to geoid modeling. The results show that the newly estimated bathymetry agrees better (by about 30 to 300 m) with the shipborne depth soundings and provides smoother residual geoid heights and gravity anomalies (by about 8–20%) than those from global models. Finally, the achieved accuracy in geoid modeling ranges between 6 and 10 cm (1σ).     

Improvement in the Determination of the Marine Geoid by Estimating the Bathymetry from Altimetry and Depth Soundings

The contribution of bathymetry to the estimation of gravity field related quantities is investigated in an extended test area in the Mediterranean Sea. The region is located southwest of the island of Crete, Greece, bounded between 33^ˆ ≤ ϕ ≤ 35^ˆ and 15^ˆ ≤ λ ≤ 25^ˆ. Gravity anomalies from the KMS99 gravity field and shipborne depth soundings are used with a priori statistical characteristics of depths in a least-squares collocation procedure to estimate a new bathymetry model. Two different global bathymetry models, namely JGP95E and Sandwell and Smith V8, are used to derive the depth a priori statistical information, while the estimated model is compared against both the global ones and the shipborne depth soundings to assess whether there is an improvement. Various marine geoid models are estimated using ERS1 and GEOSAT Geodetic Mission altimetry and shipborne gravity data. In that process, the effect of the bathymetry is computed using both the estimated and the original depths through a residual terrain modeling reduction. The TOPEX/Poseidon Sea Surface Heights, known for their high accuracy and precision, and the GEOMED solution for the geoid in the Mediterranean are used as control for the validation of the new geoid models and to assess the improvement that the estimated depths offer to geoid modeling. The results show that the newly estimated bathymetry agrees better (by about 30 to 300 m) with the shipborne depth soundings and provides smoother residual geoid heights and gravity anomalies (by about 8-20%) than those from global models. Finally, the achieved accuracy in geoid modeling ranges between 6 and 10 cm (1σ).     

重力异常的视密度反演推断扬子与华北块体在南黄海地区的地质界线

通过对南黄海及周边船载测量重力数据、现代卫星测高重力数据、陆地测量重力数据以及地球重力模型等多来源重力数据进行整理和融合,对南黄海地区布格重力异常特征进行了分析。利用基于切割法和插值迭代法的大深度位场向下延拓技术的视密度反演法对南黄海布格重力进行处理和解释,结合山东和朝鲜半岛的地质信息,推断了扬子块体和中朝块体在南黄海海域的界限。中朝和扬子块体在海上的分界位置为千里岩北缘断裂,向NE延伸到125°E后,折向NNW向并在36°N附近再次转为NE向进入朝鲜半岛并经过洪城南侧继续以NE向延伸。扬子的北缘为苏鲁造山带,穿过黄海与京畿造山带相接,苏鲁超高压变质带对应洪城杂岩。     

Geoid height anomalies over the Cape Verde Rise

Geoid height anomalies, as determined by satellite altimetry, suggest that the Cape Verde Rise is in local isostatic equilibrium, supported by a low-density root of altered lithosphere. A depth anomaly map shows the Cape Verde Rise to be approximately 1600 km wide and 2km high. Removal of a quadratic surface from the observed geoid heights leaves a residual positive anomaly with the same shape as the rise and an amplitude of about 8 m. The ratio of residual geoid height anomaly to depth anomaly is consistent with an isostatic root only 40 km deep on average.     

海区沉积盆地的三维变密度基底反演及其应用

建立了适合海区沉积盆地基底反演的密度模式，推导了对应的重力异常精确公式（近区计算公式）和近似公式（远区计算公式），提出模拟沉积盆地复杂变化的密度差－深度曲线的分段拟合方法。建立了一种利用地震波速度谱、重力和水深资料反演沉积盆地基底的方法。利用地震波速度谱的资料对厦澎凹陷的密度－深度关系进行了统计分析，并反演了厦澎凹陷的重力基底，认为整个厦澎凹陷可能由不同构造格局的东北和西南两部分组成，并分别讨论了这两部分的基底特征。     

Calculation of geoid undulations and gravity anomalies in the Northwest Pacific by using the Topex/Poseidon and Geosat altimeter data

INTRODUCTIONThegeoidistheiargeopotentials~econfidingmostlywiththemeanseasurfaceandisdenotedastheheightrelativetotheidealelliPSes~eoftheearth.Thegeoidundulationsinglobalaceareupto100m.TheunevenstructureOftheearthgivesrisetotheunevenfeatureofthecitysot...     

南极罗斯海重力场特征及莫霍面深度反演

本文基于中国南极考察第30航次、第32航次所获得的实测重力资料,结合NGDC资料,开展12个航次重力场数据的平差融合工作,全部386个交点平差后标准差减小为±1.53×10−5 m/s2,与卫星重力差值平均值为1.49×10−5 m/s2,均方差为±3.81×10−5 m/s2,并在此基础上采用频率域界面反演法计算莫霍面深度。研究发现,与沉积盆地对应重力异常低值相悖,在罗斯海北部盆地、维多利亚地盆地、中央海槽、东部盆地4个主要盆地腹地却表现为重力异常高值,跨度达100 km以上。莫霍面深度分布整体呈南深北浅之势,范围为10~28 km。伴随着罗斯海西部盆地的多次拉张及岩浆活动,该区域的地壳厚度和莫霍面深度高值和低值相间分布,并表现出越来越大的差异性。综合剖面结果表明,罗斯海重力异常值的长波长变化与莫霍面的起伏呈正相关关系,但是反演的莫霍面深度与区域重力场特征并非完全对应,所以岩浆底侵和地壳侵入仍不足以导致罗斯海盆地的重力异常或盆地几何形状。     

双曲面模型深海三维定位数据模拟及精度分析

结合声波射线传播规律提出了双曲面模型水声定位的数据模拟方法。利用该方法模拟一定浮标网形下声源位于某区域约2 000m水深处10km×10km范围内的声线双曲面模型定位观测数据,进行定位解算,并将截止角的概念引入水声定位中。定位结果显示,大部分区域水下DGPS定位系统的外符合定位精度在水平方向优于5m,垂直方向优于10m;网形中间区域精度较高,三维定位精度为亚米级。引入角度限制后,定位区域内外符合定位精度水平方向优于1m,垂直方向优于4m。