EGM2008模型与地表观测重力数据的融合及其在川西地区的应用

依据四川盆地和其西部高原地区1549个地面重力实测点数据,利用基于‘消去-恢复’思想的最小二乘配置方法,将EGM2008地球重力场模型给出的自由空气重力异常与地面实测自由空气重力异常两种数据进行融合,并对其适用性加以验证分析.计算结果表明,在四川高原地形起伏剧烈的山地,利用采样间隔4 km的地面测点,通过该数据融合方法能将模型值与实测值的标准差从41.9 mGal(10~(-5)ms~(-2))提高到11.6 mGal,平均差异从-105.4 mGal提高到-0.5 mGal;在龙泉山以东较为平坦的盆地区域,间隔8 km的地面测点密度就可以将模型标准差从5.2 mGal提高到1.9 m Gal,平均差异从-16.8 mGal提高到-0.1 mGal.上述不同地面测点分布密度对数据融合结果影响的分析表明,在地形起伏的山地,要将模型标准差控制在10 mGal左右,必须将地面测点密度控制在2 km以内;在平坦的盆地或平原地区,地面测点60 km间隔采样就可以将EGM2008模型值标准差校正到5 mGal左右.     

用GRACE卫星跟踪数据反演地球重力场

利用141天GRACE卫星观测资料，包括K波段、星载加速度和卫星轨道数据，反演了80阶地球重力场模型IGGGRACE01S，该模型在半波长为500km的空间分辨率上，确定大地水准面的精度约为0012m，中长波(<80阶)精度优于重力卫星发射以前研制的重力场模型. 与EIGEN_GRACE02S、EIGEN_CHAMP03S和EGM96模型的位系数相比，该模型系数最接近于EIGEN_GRACE02S，与另两个模型差异较大. 比较几种模型确定的全球重力异常和大地水准面起伏，结果发现IGGGRACE01S与EIGEN_GRACE02S模型的计算结果比较接近，与EGM96模型结果差异较大，差别较大地区主要在南极等地区. 对于中国大陆，比较IGGGRACE01S模型（前72阶）计算的重力异常和NIMA重力异常数据（25°×25°网格），两者之间的标准偏差为48mGal.     

川西地区实测自由空气重力异常与EGM2008模型结果的差异

在四川盆地西部建立了一个由302个观测站组成的区域观测网,并进行了高精度流动重力与GPS观测,其目的是获取区域自由空气重力异常(简称重力异常)的分布特征,并对EGM2008模型在该区域的结果进行验证分析.研究区域实测的重力异常总体为负值,由西到东逐步从-160×10^-5ms^-2平缓变化到-60×10^-5ms^-2左右.EGM2008地球重力场模型揭示的模型重力异常较好地勾画出研究区域的总体地形分布形态,龙泉山脉以及四川盆地与青藏高原的边界皆存在明显的模型重力梯度带.研究区东南部的模型重力异常大约为-50×10^-5ms^-2左右,但在龙泉山西部成都平原地区,模型重力异常则达到-120×10^-5ms^-2左右.在区域观测网内绝大部分观测点上,模型与实测重力异常的差值几乎为一个常数(-10×10^-5ms^-2左右),说明EGM2008地球重力场模型可较好地反映实际重力场的空间分布形态.如果配以一定数量的地面观测数据进行整体调节,EGM2008地球重力场模型就可以较真实地反映实际地球重力场.     

引入最优自适应因子的状态模型法航空重力测量

引入最优自适应比例因子以改善状态模型法航空重力测量的精度,并尝试将其应用到我国困难地区的重力测量.把重力扰动当作状态量引入Kalman滤波进行最优估计,并引入最优自适应因子调节状态信息的权阵,提高重力扰动的最终解算精度.利用新疆地区不同航次和航高的实测数据,计算了垂直向下方向上的重力扰动.与全球重力场模型EGM2008的对比分析表明,差值中误差在10mGal左右,接近国家在困难地区重力测量精度的限差要求.     

应用多种来源重力异常编制中国海陆及邻区空间重力异常图及重力场解读

本文通过分析陆地实测空间重力异常数据、海洋船载测量空间重力异常数据、卫星测高重力异常,布格重力异常数据、EGM2008地球重力模型数据等多种来源数据的性质和精度,并对相关数据进行对比,研究了编制1:500万中国海陆空间重力异常图的数据使用方案和技术方法.在地形较为平坦、实测数据分布均匀的陆区,使用实测数据,在地形复杂,实测数据稀少以及没有实测数据的陆区或岛屿,利用布格重力异常反推空间异常的方法合成平均空间重力数据,西藏地区的数据对比实验证明合成平均空间重力异常数据是一种有效的数据补充.利用三观测列方差分解法在南海地区对船载测量空间重力数据和美国SS系列及丹麦DNSC08GRA卫星重力数据进行了方差分解计算,结果表明不同来源的卫星测高重力数据具有很大的一致性,数据精度较以往有了很大的提高.海区空间重力数据使用原则是在船载重力测量数据校准下,全面使用卫星测高重力数据进行编图.海陆过渡区的异常处理应以EGM2008地球重力模型重力场为基准参考场,实现海陆异常平缓过渡,无缝连接.对中国海陆空间重力异常场进行了小波变换处理,对空间重力异常场进行了解读,勾画出三横四竖的一级重力梯级带及其所围限的8个一级重力异常区,并划分了二级重力异常区和梯级带,为块体构造学体系中大地构造格架的建立提供了地球物理证据.     

基于矩谐分析的区域重力场建模

建立高精度高分辨率区域地球重力场模型是大地测量学科的重要科学目标之一,可为固体地球物理、地球动力学、地震学以及资源勘探等领域提供基础信息,也是当前大地测量现代化发展实现GNSS测定海拔高程的重大需求.提出了一种基于矩谐分析的区域重力场建模方法,在局部直角坐标系中求解地球引力位的Laplace方程,推导了扰动位、重力异常、重力扰动、大地水准面差距和垂线偏差的矩谐展开式以及矩谐系数阶方差和阶误差的计算公式,给出了利用重力观测值进行矩谐分析求解矩谐系数的数学模型和算法,为抑制矩谐分析中存在的周期延拓边界效应,提出了扩展计算区域范围的策略.利用EGM2008重力位模型生成模拟地面重力和航空重力观测值,加入标准差为2 mGal的高斯白噪声,分别设计基于地面重力数据和基于航空重力数据的模拟数值试验,对矩谐分析建模方法进行了验证与分析.数值结果表明基于矩谐分析构建的区域重力场模型可靠且可达很高的精度,由地面和航空重力数据计算的2.5′×2.5′大地水准面精度分别达到1 cm和1.4 cm,4 km飞行高度处航空重力观测值的向下延拓误差仅为3.1mGal,矩谐分析可为区域重力场精细结构的逼近提供一种新的选择.     

基于多种方法构建似大地水准面模型推估特性分析

针对构建似大地水准面中采用的重力水准法、EGM2008水准法、GPS水准法,结合地形复杂程度、面积大小、重力似大地水准面分辨率与精度等特点,首次评价了基于多种方法构建的似大地水准面模型的外推性与真实性.采用地形复杂区域S、平原区域N以及沿海区域Q构建区域似大地水准面模型.通过时各方法得到的似大地水准面模型推估特性分析比较,结果表明:就反应似大地水准面起伏的真实情况而言,重力水准法精度最高且外推能力最强;EGM2008水准法精度次之,外推能力较强外推范围在10 km左右;GPS水准法精度最弱,外推能力最差.构建高精度高分率1 cm级似大地水准面时须采用基于grave/level法.     

地形/均衡重力场模型构制及其拟稳分析

重力空白区数据填补的一个主要方法是基于地壳均衡理论进行的,该方法亦用于EGM系列模型的构建中.本文研究了地形数据在构制地形/均衡重力场模型中的应用,分析了补偿深度对Airy位模型和面凝聚位模型的影响,给出二者的最佳补偿深度分别为50 km和40 km.以纯卫星重力模型为参考,后者在前120阶的精度要高于前者,但在121~250阶的精度较低,组合模型精度高于单一模型精度.对地形/均衡地球重力场模型进行了EGM2008拟稳分析,研究了不同分辨率基准的拟稳效果,分析表明：30'分辨率的拟稳基准所得拟稳模型对应的阶方差与参考阶方差曲线直到360阶都有较好的一致性,以EGM2008为基准,其相对累计大地水准面高误差在140阶时为6.83cm,相对累计重力异常误差在220阶时为1.10 mGal.     

基于地面重力的卫星重力梯度检校方法

卫星重力梯度仪在轨检校是提高梯度模式重力卫星观测质量的关键.本文面向中国未来梯度模式重力卫星规划任务,研究提出一种基于地面重力的卫星在轨检校方法,该方法顾及卫星设计指标,从地面先验重力数据的精度、空间尺度以及卫星重力梯度仪的观测噪声等卫星检校要素开展分析研究,成功实现地面数据与卫星观测数据的比对检校.研究结果表明,在12°×12°的地面检校区域内,利用精度优于1 mGal的地面重力观测数据即可实现噪声低于10 mE的重力梯度仪高精度在轨检校.     

基于变差函数的插值方法计算华北地区重力场变化

定点重复重力测量是获取区域重力场变化的主要手段之一, 重力场特征与地形起伏、 构造走向等因素相关。 以华北地区为例, 考虑区域内地形和构造的北东向分布规律, 从EGM2008重力模型中拟合各向异性变差函数参数, 利用变差函数网格化插值, 对华北地区2009—2013年期间重力场观测数据进行网格化重建, 获取华北地区重力场时空变化结果。 研究结果表明基于区域地形、 构造特征的各向异性变差函数插值方法, 获得的空间重力场变化在重力异常梯级带上更加明显, 重力变化与活动构造分布具有更好的一致性。 本文研究方法对于恢复区域时空重力场异常具有重要意义, 有助于提高应用重力资料划分潜在地震危险区空间位置的精度, 为华北地区震情研判和构建地震预报定量指标体系提供可靠的地球物理场数据。     

Error sources and data limitations for the prediction of surface gravity: a case study using benchmarks

Gravity-based heights require gravity values at levelled benchmarks (BMs), which sometimes have to be predicted from surrounding observations. We use the Earth Gravitational Model 2008 (EGM2008) and the Australian National Gravity Database (ANGD) as examples of model and terrestrial observed data respectively to predict gravity at Australian National Levelling Network (ANLN) BMs. The aim is to quantify errors that may propagate into the predicted BM gravity values and then into gravimetric height corrections (HCs). Our results indicate that an approximate ±1 arc-min horizontal position error of the BMs causes maximum errors in EGM2008 BM gravity of ~22 mGal (~55 mm in the HC at ~2200 m elevation) and ~18 mGal for ANGD BM gravity because the values are not computed at the true location of the BM. We use RTM (residual terrain modelling) techniques to show that ~50% of EGM2008 BM gravity error in a moderately mountainous region can be accounted for by signal omission. Non-representative sampling of ANGD gravity in this region may cause errors of up to 50 mGals (~120 mm for the Helmert orthometric correction at ~2200 m elevation). For modelled gravity at BMs to be viable, levelling networks need horizontal BM positions accurate to a few metres, while RTM techniques can be used to reduce signal omission error. Unrepresentative gravity sampling in mountains can be remedied by denser and more representative re-surveys, and/or gravity can be forward modelled into regions of sparser gravity.     

A study of regional gravity field recovery from GOCE vertical gravity gradient data in the Auvergne test area using collocation

Gravity field and steady-state Ocean Circulation Explorer (GOCE) is the first satellite mission that observes gravity gradients from the space, to be primarily used for the determination of high precision global gravity field models. However, the GOCE gradients, having a dense data distribution, may potentially provide better predictions of the regional gravity field than those obtained using a spherical harmonic Earth Geopotential Model (EGM). This is investigated in Auvergne test area using Least Squares Collocation (LSC) with GOCE vertical gravity gradient anomalies (T_zz), removing the long wavelength part from EGM2008 and the short wavelength part by residual terrain modelling (RTM). The results show that terrain effects on the vertical gravity gradient are significant at satellite altitude, reaching a level of 0.11 E?tv?s unit (E.U.) in the mountainous areas. Removing the RTM effects from GOCE T_zz leads to significant improvements on the LSC predictions of surface gravity anomalies and quasigeoid heights. Comparison with ground truth data shows that using LSC surface free air gravity anomalies and quasi-geoid heights are recovered from GOCE T_zz with standard deviations of 11 mGal and 18 cm, which is better than those obtained by using GOCE EGMs, demonstrating that information beyond the maximal degree of the GOCE EGMs is present. Investigation of using covariance functions created separately from GOCE T_zz and terrestrial free air gravity anomalies, suggests that both covariance functions give almost identical predictions. However, using covariance function obtained from GOCE T_zz has the effect that the predicted formal average error estimates are considerably larger than the standard deviations of predicted minus observed gravity anomalies. Therefore, GOCE T_zz should be used with caution to determine the covariance functions in areas where surface gravity anomalies are not available, if error estimates are needed.     

Combined Regional Gravity Model of the Andean Convergent Subduction Zone and Its Application to Crustal Density Modelling in Active Plate Margins

The Central Andean subduction system is one of the most active geological structures on Earth. Although there have been a few previous studies, the structure and dynamics of the system are still not well understood. In the present study, we determine a combined regional gravity model of the Andean convergent subduction region for constraining lithospheric models. After a thorough validation and cleaning of the terrestrial gravity and height databases, the method of Least Squares Collocation was applied to consistently combine terrestrial and satellite gravity data, putting much emphasis on the stochastic modelling of the individual data components. As a result, we computed the first high-resolution regional gravity model of the study region that includes GOCE satellite gravity information. The inclusion of GOCE is an essential distinction from the independent global gravity model EGM2008. Validation against EGM2008 reveals that our regional solution is very consistent in regions where terrestrial gravity data are available, but shows systematic differences in areas with terrestrial data gaps. Artefacts in the EGM2008 of up to 150 mGal could be identified. The new combined regional model benefits from the very homogeneous error characteristics and accuracy of GOCE gravity data in the long-to-medium wavelengths down to 80–100 km. Reliable density modelling became possible also in the region of Central Andes, which lacks terrestrial gravity data. Finally, density models were adapted to fit the new regional gravity field solution. The results clearly demonstrate the capabilities of GOCE to better constrain lithospheric models.     

GRACE和地面重力测量监测到的中国大陆长期重力变化

自2002年以来,GRACE卫星探测计划可提供高精度的时变地球重力场,用以探测地球系统的物质分布.自1998年中国大陆重力监测网建立以来,利用FG5绝对重力仪和LCR-G型相对重力仪每2年对该网进行重复测量获取重力场时变信息.基于此,本文利用GRACE和地面重力测量获得了中国大陆重力场的长期年变率,利用位错理论根据USGS发布的断层模型计算了2008年汶川M_s8.0级地震的同震重力变化并进行了300 km高斯滤波.GRACE卫星重力和地面重力结果均表明华北地区地下水流失严重,在绝对重力基准站上,GRACE卫星重力与绝对重力变化率较为一致,汶川区域的地面重力变化结果可视为大地震前兆信息.     

Local geoid determination based on airborne gravity data

A mathematical model used for determination of a local geoid model by combining airborne gravity disturbances and the Earth Gravitational Model 2008 (EGM08) is shortly reviewed. The precision of the estimated local geoid model of Taiwan is tested by its comparison with the “real” geoid at Global Satellite Navigation Systems (GNSS)/levelling points. The same comparison at GNSS/levelling points is done for the geoid evaluated only by using EGM08. Conclusions concerning a rate of improvement of the “global” geoid from EGM08 using the “local” geoid from airborne gravity data are presented.     

Improving airborne strapdown vector gravimetry using stabilized horizontal components

Integrating the deflections of the vertical along the flight line can yield geoid profiles which are valuable in the study of geodesy and geophysics, fortunately, the deflections can be measured directly by vector gravimetry. Airborne vector gravimetry using a Strapdown Inertial Navigation System and the Global Navigation Satellite System (SINS/GNSS) has shown promising results in previous studies. However, the quality of the SINS and GNSS is a major limitation; in particular, the attitude errors induced by the gyros will result in large measurement errors to the horizontal components of the gravity disturbance, and these measurement errors represent the behavior of low-frequency trend. An airborne vector gravimetry method used to remove the bias and low-frequency trends in the gravity disturbance estimated for each survey line has been developed. This method uses the horizontal components of the gravity disturbance computed from EGM2008 (Earth Gravitational Model 2008) as a reference. Firstly, the horizontal measurement results obtained from the gravimeter are divided into high- and low-frequency components according to the resolution of the EGM2008, and then, the bias and low-frequency trends of the low-frequency components are corrected using a linear fit to the EGM2008 reference data. Finally, the ultimate results can be acquired after combining the high-frequency components and the corrected low-frequency components. The data used was obtained from the SGA-WZ, which is the first strapdown airborne gravimeter developed in China. The results of this method are promising. The internal accuracy of the gravity disturbance's horizontal components for repeated survey lines exceeds 3.5 mGal, and the corresponding resolution is approximately 4.8 km based on 160-s data smoothing and an airplane averaging speed of approximately 216 km/h. After applying the WCF (Wavenumber Correlation Filter), the internal accuracy of the horizontal components exceeds 2 mGal. This can satisfy the requirement of the application in geodesy and solid earth geophysics.     

Geoid determination by spectral combination of an Earth gravitational model with airborne and terrestrial gravimetry — a theoretical study

Today air-gravimetry is a versatile technique to quickly collect gravity data over large regions, where terrestrial gravity data are sparse and/or of poor quality. The method requires the data to be downward continued to sea level for use in geoid determination, an inverse problem operation that calls for smoothing of the data and/or the kernel function involved (in either spectral or space domain). In this purely theoretical study we avoid this separate computational step by performing direct geoid estimation by so-called spectral combination/filtering of the data, which includes terrestrial gravimetry, airgravimetry, an Earth Gravitational Model (EGM) as well as their signal and error degree variances. Each derived geoid estimator is presented as the sum of one or two integral formulas and the harmonic series of the EGM together with the expected mean square error of the estimator. The article is limited to a theoretical study, leaving its practical tests for future investigation.