基于DTU重力数据检查海洋重力测量成果的质量

在海上复杂的环境下进行重力测量,重力仪及技术人员的操作难免会出现各种各样的问题,导致测量的成果精度降低,甚至不符合规范要求。目前检查海洋重力测量成果质量的主要方法是内符合法,主要采用主测线、检测线交叉点的差值,或重复测线相邻点的差值来评估测量成果的质量。介绍了内符合法存在的弊端,比较了DTU重力数据与重力仪实测数据的精度,叙述了基于DTU重力数据检查海洋重力测量成果较大系统差的可行性,以模型数据为外符合数据源来检查实测成果的质量,作为内符合法的补充,是十分必要的。     

几种卫星测高海洋重力场模型精度比较分析

以中国南海部分海域作为试验区,依据测区船载重力测量结果分别对S&S V28.1、S&S V29.1、DTU15和DTU17卫星测高海洋重力场模型开展了精度评估与分析工作。首先利用EIGEN6C4超高阶位系数模型计算得到的参考重力异常,使用多项式拟合技术对船载重力异常进行了长波误差改正;然后借鉴3σ准则进一步控制船载重力测量结果质量;最后以RMS作为精度主要评价指标,完成了对4种海洋重力场模型精度评估。结果表明,4种海洋重力场模型检核精度整体相当,其中S&S V29.1检核精度略高于S&S V28.1,4种重力异常模型在试验海域精度大致在1.5~6mGal之间。试验结论对于后续重力场模型应用具有一定参考价值。     

卫星测高反演重力场在西太平洋海域的应用

利用高精度卫星测高数据反演重力是获取海洋重力异常的重要方法。本文介绍了自2008年以来新一代高精度测高卫星的相关信息。在此基础上,联合多颗新型卫星,进行交叉点平差,通过中心差分代替常规Sandwell方法计算垂线偏差,利用改进的逆Vening-Meinesz法和Stokes数值反解法对西太平洋某区域进行了重力场反演。反演结果与NGDC提供的船测重力值比较表明,逆Vening-Meinesz法计算误差小于Stokes数值反解法。通过添加随机噪声数值实验,进一步证明逆Vening-Meinesz法对于测高数据误差的抗差能力更强,在研究区的重力异常计算结果可信度较高。     

利用多代卫星测高数据计算中国近海及邻域重力异常

为提高海洋重力场数据的精度和空间分辨率,联合Jason-1/2、T/P、Envisat、ERS-1/2、Geosat等多代卫星测高数据计算中国近海及邻域(0°~42°N,100°~140°E)2′×2′重力异常。对卫星测高数据分别进行共线处理和自交叉点平差,并以T/P卫星测高数据为基准进行多星数据联合平差,有效削弱了卫星测高数据的时变影响和不协调性;利用逆Vening-Meinesz公式计算重力异常,与船测重力相比,均方根误差为5.4 mgal。结果表明,通过引入高精度的卫星测高数据,结合多项平差处理手段,提高了海洋重力异常的计算精度。     

海洋重力测量动态环境效应分析与补偿

针对使用小型测量船搭载摆杆型海洋重力仪获取数据质量不高的问题,在深入分析海洋环境动态效应误差特性基础上,提出了一种基于互相关分析的交叉耦合效应修正法,对高动态海洋重力测量数据实施综合误差补偿和精细处理。使用典型恶劣海况条件下的观测数据对该方法的有效性进行了验证,结果显示,重力测线成果内符合精度从原先的±9.35×10-5m/s2大幅提升到±1.43×10-5m/s2,同时使用卫星测高反演重力对精细处理结果的可靠性进行了外部检核,外符合精度也从原先的±7.73×10-5m/s2提高到±5.63×10-5m/s2。     

卫星测高中的垂线偏差法

卫星测高中的垂线偏差法是当前利用卫星测高技术研究海洋重力场的最优方法,包括利用卫星测高数据计算垂线偏差和利用该垂线偏差确定海洋重力场两部分。研究了Sandwell、Olgiafi、Hwang测高垂线偏差的计算方法和Molodenskii、Hwang利用测高垂线偏差确定海洋重力场的基本原理,分析比较了上述方法的异同,为科学地利用卫星测高资料反演海洋重力场提供理论依据。     

海洋卫星重力异常及其应用

随着空间技术的发展 ,卫星测高精度不断提高 ,目前根据卫星测高数据得到的卫星重力异常在海区具有很高的异常分辨率 ,与船测重力资料相比其精度可达到5mGal,这将十分有助于海洋沉积盆地、板块构造及地球动力学的研究。     

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

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

海洋重力测量数据精度分析

海洋重力测量仪器工作环境的特殊性,决定了其测量数据精度分析与检测难以方便实现。在分析一般数据处理方法的基础上,研究了海洋重力测量精度分析的特点,采用一种以重复测线不符值、测线网交叉点不符值的标准差对海洋重力测量的内符合精度进行分析,并根据重复测线、交叉点内符合精度评估方法,通过对某型海洋重力测量仪器的多次航行试验,获取了该型仪器测量数据的部分精度数据,分析结果表明该方法可以全面衡量仪器性能。     

利用多代卫星测高数据计算海洋重力场

联合Jason-1/2、T/P、Envisat、ERS-1/2、Geosat等多代卫星测高数据计算中国近海及邻域（0°~42°N,100°~140°E）2′×2′重力异常。对卫星测高数据分别进行共线处理和自交叉点平差,并以T/P卫星测高数据为基准进行多星数据联合平差,有效削弱了卫星测高数据的时变影响和不协调性;利用逆Vening-Meinesz公式计算重力异常,与船测重力相比,均方根误差为5.4mGal。结果表明,通过引入高精度的卫星测高数据,结合多项平差处理手段,提高了海洋重力异常的计算精度。     

Determining the Gravitational Gradient Tensor Using Satellite-Altimetry Observations over the Persian Gulf

With the advent of satellite altimetry in 1973, new scientific applications became available in oceanography, climatology, and marine geosciences. Moreover, satellite altimetry provides a significant source of information facilitated in the geoid determination with a high accuracy and spatial resolution. The information from this approach is a sufficient alternate for marine gravity data in the high-frequency modeling of the marine gravity field quantities. The gravity gradient tensor, consisting of the second-order partial derivatives of the gravity potential, provides more localized information than gravity measurements. Marine gravity observations always carry a high noise level due to environmental effects. Moreover, it is not possible to model the high frequencies of the Earth's gravity field in a global scale using these observations. In this article, we introduce a novel approach for a determination of the gravity gradient tensor at sea level using satellite altimetry. Two numerical techniques are applied and compared for this purpose. In particular, we facilitate the radial basis functions (RBFs) and the harmonic splines. As a case study, the gravitational gradient tensor is determined and results presented in the Persian Gulf. Validation of results reveals that the solution of the harmonic spline approach has a better agreement with a theoretical zero-value of the trace of the Marussi gravitational gradient tensor. However, the data-adaptive technique in the RBF approach allows more efficient selection of the parameters and 3-D configuration of RBFs compared to a fixed parameterization by the harmonic splines.     

利用卫星测高垂线偏差计算其他重力场元的特性分析

利用卫星测高技术确定海洋重力场,垂线偏差数据作为导出观测量在实际工作中被普遍采用。利用物理大地测量边值问题的定义以及扰动位在球面边界条件下的解,给出了由垂线偏差计算大地水准面高、重力异常和扰动重力的公式。分析了不同积分计算公式在重力场阶谱表达形式下对垂线偏差误差的抑制作用,也分析了不同积分核函数的变化特性,得出基本结论:在利用卫星测高数据求解海洋重力场时,当以格网化海面垂线偏差数据计算重力场参数时,求解的大地水准面高的有效性和稳定性优于重力异常和扰动重力。     

System Ⅱ型海洋重力仪静态观测结果与分析

根据青岛海洋地质研究所新购置System II型海洋重力仪的静态观测实验,分析了重力数据的变化情况,检测到该仪器单日记录数值较为稳定,但伴随有小于±0.1mGal左右的震荡。同时分析了该设备对固体潮的反映,静态观测重力值所展现的变化特征和周期性与固体潮相一致。鉴于其能明显反应固体潮对静态观测重力值的影响,可认为其静态观测精度优于0.2mGal。并测定了重力仪月漂值小于国家标准3mGal/月。同时比较讨论了与KSS31M型重力仪在中科院海洋研究所的静态观测数据的区别与联系。     

Crustal analysis of the Ulleung Basin in the East Sea (Japan Sea) from enhanced gravity mapping

To facilitate geological analyses of the Ulleung Basin in the East Sea (Japan Sea) between Korea and Japan, shipborne and satellite altimetry-derived gravity data are combined to derive a regionally coherent anomaly field. The 2-min gridded satellite altimetry-based gravity predicted by Sandwell and Smith [Sandwell DT, Smith WHF (1997) J Geophys Res 102(B5):10,039–10,054] are used for making cross-over adjustments that reduce the errors between track segments and at the cross-over points of shipborne gravity profiles. Relative to the regionally more homogeneous satellite gravity anomalies, the longer wavelength components of the shipborne anomalies are significantly improved with minimal distortion of their shorter wavelength components. The resulting free-air gravity anomaly map yields a more coherent integration of short and long wavelength anomalies compared to that obtained from either the shipborne or satellite data sets separately. The derived free-air anomalies range over about 140 mGals or more in amplitude and regionally correspond with bathymetric undulations in the Ulleung Basin. The gravity lows and highs along the basin’s margin indicate the transition from continental to oceanic crust. However, in the northeastern and central Ulleung Basin, the negative regional correlation between the central gravity high and bathymetric low suggests the presence of shallow denser mantle beneath thinned oceanic crust. A series of gravity highs mark seamounts or volcanic terranes from the Korean Plateau to Oki Island. Gravity modeling suggests underplating by mafic igneous rocks of the northwestern margin of the Ulleung Basin and the transition between continental and oceanic crust. The crust of the central Ulleung Basin is about a 14–15 km thick with a 4–5 km thick sediment cover. It may also include a relatively weakly developed buried fossil spreading ridge with approximately 2 km of relief.     

两种测高重力异常反演海底地形方法比较

针对基于测高重力异常反演海底地形理论众多、选取标准无法确定的情况,利用中国南海海域内的测高重力异常和船测水深数据研究比较了重力地质法(GGM)和SmithSandwell (SAS)法两种精度高、计算速度相对较快的海底地形反演理论。其中,GGM方法的密度差异常数Δρ由向下延拓技术确定为2.15 g·cm-3,SAS方法采用移去-恢复技术得到反演波段内重力异常和水深数据。结果表明:测线分布条件一定时,水深多在-1 000 m左右或反演区域岛礁、海山等复杂海底地形较多时选取SAS方法,水深主要在-3 000 m以深的区域或海底地形复杂程度不高时选取GGM方法则能获取更好的效果,其效果最优处与船测水深在检核点处的差值最优平均值能达-0.61 m,标准差可达14.67 m。     

太平洋卫星测高重力场与地球动力学特征

通过多卫星测高数据的综合处理，获得西太平洋卫星测高重力场，进行不同尺度、深度构造动力信息的分离，探讨诸边缘海盆的地球动力学问题。测高大地水准面反映了研究区板块相互作用的特点，其高频成分可以刻画各海盆的构造特征。测高空间重力异常也可刻画陆架构造及盆地分布，由其推算出的海底地形含有大量的海底构造信息。各边缘海盆的莫霍面埋深具有往南变浅的趋势，与菲律宾海各海盆的莫霍面埋深大致相当，说明岛弧两侧的构造动力强度基本相似。大尺度地幔流应力场总体上反映了欧亚板块向东南蠕散和太平洋板块向北西扩张的特点；日本海北侧和南海巽他陆架的中尺度上地幔对流与地幔柱之间有着密切关系，西菲律宾海的上地幔对流强化了日本-琉球-台湾-菲律宾岛弧的活动强度；小尺度地幔流主要限于软流圈层内部，在各海盆分散，而在冲绳海槽和马里亚纳海槽则会聚，可与均衡重力异常类比。还讨论了大、中、小地幔流体系的特点及相互之间的关系，籍以阐明海盆及海槽演化的地球动力学过程。     

卫星测高技术在监测海啸现象中的研究及应用

利用卫星测高技术,通过共线平差和沿轨、垂直轨迹测高大地水准面梯度的计算,得到了海啸发生2个小时后的海面波高和海面梯度异常现象,为利用卫星测高技术监测海啸、模拟海啸现象提供了依据。     

Retracking of SARAL/AltiKa Radar Altimetry Waveforms for Optimal Gravity Field Recovery

The accuracy of the marine gravity field derived from satellite altimetry depends on dense track spacing as well as high range precision. Here, we investigate the range precision that can be achieved using a new shorter wavelength Ka-band altimeter AltiKa aboard the SARAL spacecraft. We agree with a previous study that found that the range precision given in the SARAL/AltiKa Geophysical Data Records is more precise than that of Ku-band altimeter by a factor of two. Moreover, we show that two-pass retracking can further improve the range precision by a factor of 1.7 with respect to the 40 Hz-retracked data (item of range_40 hz) provided in the Geophysical Data Records. The important conclusion is that a dedicated Ka-band altimeter-mapping mission could substantially improve the global accuracy of the marine gravity field with complete coverage and a track spacing of <6 km achievable in ～1.3 years. This would reveal thousands of uncharted seamounts on the ocean floor as well as important tectonic features such as microplates and abyssal hill fabric.     

Oceanic gravity by analytical inversion of Hotine's formula

An analytical inversion of the Hotine formula is developed using fast Fourier transform techniques. Detailed mathematical derivations are used to explain the concepts behind the inverse transformation. Three modifications of the analytical inversion of the Hotine formula are compared and tested using both synthetic data from the OSU91A geopotential model and real GEOSAT altimetry data from the Exact Repeat Mission. The stability of this inverse Hotine approach is investigated using simulated data, and numerical tests are done to quantify the stability of this approach. The approach seems to be numerically stable without employing any stabilization technique. Estimated gravity information from GEOSTAT altimetry data is compared to marine gravity data from shipboard measurements in the Orphan Knoll area. The standard deviations and mean values of the differences between satellite and marine gravity disturbances are 8.2 and 2.9 mGal for the planar approximation, 9.2 and 3.7 mGal for the spherical approximation, and 9.5 and 1.9 mGal for the Molodenskii‐like approximation, respectively, indicating that latitude‐dependent errors affect the latter two approximations. Such errors could be eliminated by performing the calculations by the rigorous one‐dimensional fast Fourier transform (FFT) technique, and any data noise could be filtered out by utilizing covariance knowledge about the input geoid undulations and their errors. Simulation studies also showed that the accuracy of the techniques (for all approximations) can reach a root‐mean‐square (RMS) level of only a few mGal when proper treatment of FFT edge effects is employed and a rather wide area of results is disregarded around the edges.