GRACE detection of the medium-to far-field coseismic gravity changes caused by the 2004 Mw9.3 Sumatra-Andaman earthquake

Large earthquakes cause observable changes in the Earth’s gravity field, which have been detected by the Gravity Recovery and Climate Experiment (GRACE). Since most previous studies focus on the detection of near-field gravity effects, this study provides the results from the medium- to far-field gravity changes caused by the 2004 Sumatra-Andaman earthquake that are recorded within GRACE monthly solutions. Utilizing a spherical-earth dislocation model we documented that large-scale signals predominate in the global field of the coseismic gravity changes caused by the earthquake. After removing the near-field effects, the coseismic gravity changes show a negative anomaly feature with an average magnitude of -0.18×10-8 m·s-2 in the region ranging ～40° from the epicenter, which is considered as the 'medium ffield' in this study. From the GRACE data released by Center for Space Research from August 2002 to December 2008, we retrieved the large-scale gravity changes smoothed with 3 000 km Gaussian ffilter. The results show that the coseismic gravity changes detected by GRACE in the medium field have an average of (-0.20±0.06)×10-8 m·s-2, which agrees with the model prediction. The detection confirms that GRACE is sensitive to large-scale medium-field coseismic gravitational effects of mega earthquakes, and also validates the spherical-earth dislocation model in the medium field from the perspective of satellite gravimetry.

Construction of Earth's gravitational field model from CHAMP,GRACE and GOCE data

The basic principle of spectral combination method is discussed,and the general expressions of the spectral weight and spectral combination of the united-processing of various types of gravimetric data are shown.What＇s more,based on degree error RMS of potential coefficients,the detailed expressions of spectral combination formulae and the corresponding spectral weights in the Earth＇s gravitational field model（EGM） determination using GOCE ＋ GRACE and CHAMP ＋ GRACE ＋ GOCE are derived.The fundamental situation that ulux-champ2013 s,tongji-GRACE01,go-cons-gcf-2-tim-r5 constructed respectively by CHAMP,GRACE,GOCE data and go-cons-gcf-2-dir-r5 constructed by syncretic processing of GRACE,GOCE and LAGEOS data are explained briefly,the degree error RMS,cumulative geoid height error and cumulative gravity anomaly error of these models are calculated.A syncretic model constructed from CHAMP,GRACE and GOCE data,which is expressed by champ ＋ grace ＋ goce,is obtained based on spectral combination method.Experimentation results show that the precision of CHAMP data model is the lowest in satellite-only models,so it is not needed in the determination of syncretic models.The GRACE data model can improve the GOCE data model in medium-long wavelength,so the overall precision of syncretic model can be improved.Consequently,as many types of gravimetric data as possible should be combined together in the data processing in order to strengthen the quality and reliability with widening scope and improve the precision and spatial resolution of the computational results.

A GOCE only gravity model GOSG01S and the validation of GOCE related satellite gravity models

We compile the GOCE-only satellite model GOSG01S complete to spherical harmonic degree of 220 using Satellite Gravity Gradiometry （SGG） data and the Satellite-to-Satellite Tracking （SST） observations along the GOCE orbit based on applying a least-squares analysis. The diagonal components （Vxx, Vyy, Vzz） of the gravitational gradient tensor are used to form the system of observation equations with the band-pass ARMA filter. The point-wise acceleration observations （ax, ay, az） along the orbit are used to form the system of observation equations up to the maximum spherical harmonic degree/order 130. The analysis of spectral accuracy characteristics of the newly derived gravitational model GOSG01S and the existing models GOTIM04S, GODIR04S, GOSPW04S and JYY_GOCE02S based on their comparison with the ultrahigh degree model EIGEN-6C2 reveals a significant consistency at the spectral window approximately between 80 and 190 due to the same period SGG data used to compile these models. The GOCE related satellite gravity models GOSG01S, GOTIM05S, GODIR05S, GOTIM04S, GODIR04S, GOSPW04S, JYY_- GOCE02S, EIGEN-6C2 and EGM2008 are also validated by using GPS-leveling data in China and USA. According to the truncation at degree 200, the statistic results show that all GGMs have very similar differences at GPS-leveling points in USA, and all GOCE related gravity models have better performance than EGM2008 in China. This suggests that all these models provide much more information on the gravity field than EGM2008 in areas with low terrestrial gravity coverage. And STDs of height anomaly differences in China for the selected truncation degrees show that GOCE has improved the accuracy of the global models beyond degree 90 and the accuracies of the models improve from 24 cm to 16 cm. STDs of geoid height differences in USA show that GOSG01S model has best consistency comparing with GPSleveling data for the frequency band of the degree between 20 and 160.

An Improved Energy Balance Approach and Its Application in CHAMP Gravity Field Recovery

An efficient method for gravity field determination from CHAMP orbits and accelerometer data is referred to as the energy balance approach. A new CHAMP gravity field recovery strategy based on the improved energy balance approach IS developed in this paper. The method simultaneously solves the spherical harmonic coefficients, daily Integration constant, scale and bias parameters. Two 60 degree and order gravitational potential models, XISM-CHAMPO1S from the classical energy balance approach, and XISM-CHAMPO2S from the improved energy balance, are determined using about one year＇s worth of CHAMP kinematic orbits from TUM and accelerometer data from GFZ. Comparisons among XISM-CHAMPO1S, XISM-CHAMPO2S, EIGEN-CGO3C, EIGEN-CHAMPO3S, EIGEN2, ENIGNIS and EGM96 are made. The results show that the XISM-CHAMPO2S model is more accurate than EGM96, EIGENIS, EIGEN2 and XISM-CHAMPO1S at the same degree and order, and has almost the same accuracy as EIGEN-CHAMPO3S.

Gravity Field Variations Associated with the Buried Geological Structures: San Marcos Fault (NE Mexico) Case Study

Gravity data are sensitive to local vertical offsets across high-angle faults, where rocks with different densities are juxtaposed. Yet high densities in some Mesozoic sedimentary rocks just above the basement may smear out the subtle gravity signatures of basement faults. At this study the gravity data processing tends to avoid ill-described “black-box” techniques. The study area is situated in the Palomas site, Cuatrociénegas region, Coahuila, NE Mexico. The San Marcos Fault is at least 300 km long and has WNW-ESE trend from the central part of Nuevo León State through Coahuila, and finally to the eastern part of Chihuahua State. Gravimetric data shows that the lowest values of free air and Bouguer anomalies are in the southern part of the area, and the highest values are in the western and central part of the area. Between these parts exists a zone of high horizontal gravity gradient. Configuration of linear elements of gravity field (gradient zones) delimited the San Marcos Fault in the San Marcos valley below thickness of recent sedimentary cover. Two density models were carried out, which showed that the Cretaceous rocks are in discordant contact with the Paleo- zoic rocks that can be related to the San Marcos Fault. The density was determinate using to Nettleton’s method, which results highlight the presence of the San Marcos Fault. Density models showed that the Quaternary sediments are in direct contact with the San Marcos Fault.

EGM96和EGM2008地球重力场模型计算弹道扰动引力的比较

给出了利用地球重力场模型计算扰动引力以及重力场模型截断误差的公式,计算了36阶、360阶、2 160阶地球重力场模型在径向和水平方向的截断误差,并对36阶、360阶的EGM96和EGM2008重力场模型计算的弹道扰动引力矢量进行了分析和比较。结果表明,采用重力场模型计算扰动引力时,模型阶数和弹道点的高度越低,模型截断误差越大;EGM96和EGM2008模型的一致性随着模型阶数的升高而逐渐降低。总体来看,EGM2008模型的精度要好于EGM96模型。

EGM96与EGM2008地球重力场模型精度比较

给出了由地球重力场模型计算重力异常和垂线偏差的公式，利用36阶、360阶EGM96、EGM2008地球重力场模型计算了国内某地区格网点重力异常和地面垂线偏差，并将其与实测数据进行比较，从而对这两种模型的精度进行了分析和比较。结果表明，在表示国内某地区格网点重力异常时，EGM2008模型精度较高；在表示地面垂线偏差时，两种模型的精度相当。

EGM2008地球重力场模型在城市正常高计算中的应用研究

基于空间分辨率1′×1′EGM2008地球重力场模型数据,用Hans-Gerd Duenck-Kerst等开发的Alltrans EGM2008Calculator 1.2软件,探讨如何获取城市高程异常及加常数,并提出了EGM2008重力场模型与我国区域采用似大地水准面之间的加常数,得到一些有益的结论,可为其他工程正常高的获取提供重要的借鉴经验。

利用EGM2008模型快速构建扰动重力梯度基准图

重力梯度张量是重力位二阶导,相比重力异常能够更好反映局部区域的细节特征。因此重力梯度导航理论上能为惯性导航提供更好的辅助。重力梯度导航的关键技术之一是背景基准图的构建,推导了扰动重力梯度张量与扰动位在局部指北坐标系中的关系式,并基于EGM2008地球重力场模型构建了一块范围的海域扰动重力梯度张量基准图。为了快速构建基准图,选取了合适的勒让德函数,并将每一个梯度张量的计算式改变求和顺序来提高同一纬度圈上的计算点的计算速度。最后利用梯度张量对角线上三元素满足拉普拉斯约束条件的原理验证了所得基准图的正确性。

A three-dimensional Moho depth model beneath the Yemeni highlands and rifted volcanic margins of the Red Sea and Gulf of Aden, Southwest Arabia

Knowing Moho discontinuity undulation is fundamental to understanding mechanisms of lithosphereasthenosphere interaction, extensional tectonism and crustal deformation in volcanic passive margins such as the study area, which is located in the southwestern corner of the Arabian Peninsula bounded by the Red Sea and the Gulf of Aden. In this work, a 3D Moho depth model of the study area is constructed for the first time by inverting gravity data from the Earth Gravitational Model(EGM2008) using the ParkerOldenburg algorithm. This model indicates the shallow zone is situated at depths of 20 km to 24 km beneath coastal plains, whereas the deep zone is located below the plateau at depths of 30 km to 35 km and its deepest part coincides mainly with the Dhamar-Rada ’a Quaternary volcanic field. The results also indicate two channels of hot magmatic materials joining both the Sana’a-Amran Quaternary volcanic field and the Late Miocene Jabal An Nar volcanic area with the Dhamar-Rada’a volcanic field. This conclusion is supported by the widespread geothermal activity(of mantle origin) distributed along these channels,isotopic data, and the upper mantle low velocity zones indicated by earlier studies.

基于ico_HEALPix网格的超高阶地球重力场建模方法

本文针对传统地理网格数据剖分的重力场模型构建在高纬度区域出现数据冗余的问题,将分层等积等纬度像素化(HEALPix)网格结构引入地球重力场解算中,提出了利用二十面体HEALPix(ico_HEALPix)网格的超高阶地球重力场建模方法,实现了全球3600阶次球谐位系数的高效构建,同时针对ico_HEALPix网格在球谐分析过程中法矩阵不是严格块对角化结构的问题,设计了迭代算法,有效提高模型构建的精度。试验表明,ico_HEALPix网格数据在数据量小于地理网格500万的前提下,通过迭代方法构建的全球地球重力场模型精度可达到优于地理网格的效果,球谐位系数误差阶RMS提升1~2个数量级,还解决了地理网格南北极点畸变和数据冗余的问题,提高了网格的数据利用率。

DQM2022系列超高阶地球重力场模型构建及其精度评估

本文基于椭球谐分析,提出了超高阶地球重力场模型构建的局部积分改进迭代法和全球积分改进迭代法,解决了传统的局部积分改进法存在的局限性,有效提高了改进模型在中国地区的适用精度。选择EGM2008、EIGEN-6C4地球重力场模型作为初始模型,利用中国及周边地区最新的5′×5′实测格网平均重力异常数据,构建了DQM2022系列超高阶地球重力场模型,其完全阶次均为2190。采用地面实测重力异常、GNSS/水准、天文大地垂线偏差等数据对改进模型进行精度评估。结果显示:①相对于初始模型,改进模型表示中国地区重力场的精度显著提高,重力异常精度提高了2.4~2.8 mGal,高程异常精度提高了1.0~2.4 cm,垂线偏差精度提高了0.07″~0.15″;②在表示中国地区重力场时,基于EIGEN-6C4初始模型构建的改进模型精度最高,高程异常精度约为10.6 cm,垂线偏差精度约为2.1″。

Lithospheric Structure in the North China Craton Constrained from Gravity Field Model(EGM 2008)

A detailed knowledge of the thickness of the lithosphere in the North China craton（NCC） is important for understanding the significant tectonic reactivation of the craton in Mesozoic and Ce-nozoic.We achieve this goal by applying the newly proposed continuous wavelet transform theory to the Gravity Field Model（EGM 2008） data in the region.Distinct structural variations are identified in the scalogram image of profile Alxa-Datong（大同）-Qingdao（青岛）-Yellow Sea（profile ABC）,trans-versing the main units of NCC,which we interpret as mainly representing the Moho and lithosphere-asthenosphere boundary（LAB） undulations.The imaged LAB is as shallow as 60-70 km in the south-east basin and coastal areas and deepens to no more than 140 km in the northwest mountain ranges and continental interior.A rapid change of about 30 km in the LAB depth was detected at around the boundary between the Bohai（渤海） Bay basin（BBB） and the Taihang（太行） Mountains（TM）,roughly coincident with the distinct gravity decrease of more than 100 mGal that marks the North-South Grav-ity Lineament（NSGL） in the region.At last we present the gravity modeling work based on the spectral analysis results,incorporating with the observations on high-resolution seismic images and surface to-pography.The observed structural differences between the eastern and western NCC are likely associ-ated with different lithospheric tectonics across the NSGL.Combined with seismic tomography results and geochemical and petrological data,this sug-gests that complex modification of the litho-sphere probably accompanied significant litho-spheric thinning during the tectonic reactivation of the old craton.

深圳市1km高分辨率厘米级高精度大地水准面的确定

利用65个精度优于2cm的GPS水准数据、5213个实测重力点数据、100m分辨率的数字地形模型和WDM94地球重力场模型,采用移去 恢复技术计算了深圳市1km分辨率的大地水准面模型。将该模型大地水准面高与由29个GPS水准得到的大地水准面高进行比较,其差值的标准差为±1.4cm。

地球重力场模型研究进展和现状

首先回顾自1966年以来国内外基于传统重力测量技术建立的地球重力场模型。其次,介绍国际基于CHAMP实测数据建立的全球重力场模型。在重力场长波部分,EIGEN-CHAMP03S和EIGEN-3P模型符合较好;在中长波部分,EIGEN-1S、EIGEN-2、EIGEN-3P和EIGEN-CHAMP03S模型的符合性趋于一致。最后,介绍国内外已公布的GRACE卫星重力模型。由于GGM02S和WHIGG-GEGM01S模型采用的卫星观测值数量近似相等,因此二者精度符合较好;由于采用CHAMP、GRACE、卫星测高和地面重力的联合观测数据,EIGEN-CG01C、EIGEN-CG03C、EIGEN-GL04S1和EIGEN-5C模型精度均分别高于EIGEN-GRACE01S、EIGEN-GRACE02S、GGM01S、GGM02S、WHIGG-GEGM01S等GRACE-only模型。

EGM96,WDM94和GPM98CR高阶地球重力场模型表示深圳局部重力场的比较与评价

为计算深圳精密重力大地水准面 ,利用 62个高精度GPS水准点和 4871个实测重力点数据对EGM96,WDM94和GPM98CR全球重力场模型表示深圳局部重力场进行了比较与评价。结果表明 ,由上述 3个重力场模型计算的大地水准面高和重力异常与实测值之间存在明显的系统偏差 ,当采用GPS水准数据尽可能消除系统偏差以后 ,大地水准面高的精度得到显著提高 ,若应用移去 恢复技术确定深圳高精度大地水准面 ,则WDM94应该是首选的参考重力场模型。

珠峰地区似大地水准面精化与珠峰顶正高的确定

本文使用珠峰及其周边地区的重力数据与SRTM3、1∶50000DEM、GTOPO30地形数据,以该地区的44个GPS水准点为控制,选择国内外的EGM96、WDM94、IGG05B、DQM2000D和CG03C作为参考重力场模型,采用移去-恢复技术,首次完成了珠峰地区分辨率为2.5′×2.5′高精度似大地水准面,其精度达到±9cm,并据此推算了珠峰顶高程异常值.再利用登山线路上的新测重力点与珠峰地形数据,依据严格的重力归算理论及移去-恢复技术,完成了珠峰顶似大地水准面与大地水准面差值计算,即珠峰正常高与正高的换算.

扰动重力梯度张量单分量和组合分量最小二乘配置法模型的建立

建立了利用扰动重力梯度张量Tzz分量和Txx+Tyy、Tzz-Txx-Tyy组合分量确定地球重力场的调和分析法模型,进一步推导了扰动重力梯度张量对角线三分量的自协方差和互协方差函数的级数展开式,推导了单分量、组合分量与重力位系数之间协方差函数的实用计算公式,给出了利用单分量和组合分量解算地球重力场模型的最小二乘配置法基本原理公式.结果表明,最小二乘配置法具有一定的抗差能力,随着观测数据误差的不断增大,其恢复的重力场模型有效阶次不断降低,精度也不断下降;Tzz-Txx-Tyy组合分量解算重力场模型的精度最高,其次为Tzz分量,Txx+Tyy组合分量最差.

机动加速度辅助的航姿系统扩展卡尔曼滤波

低精度航姿系统一般采用基于重力场和地磁场矢量观测的6态扩展卡尔曼滤波,载体长时间机动时,模型存在较大误差,无法保持稳定的姿态精度。将载体机动加速度描述为当前统计模型,提出了一种机动加速度辅助的9态扩展卡尔曼滤波,其状态变量包括三个姿态角误差、三轴陀螺零偏误差和三轴载体加速度误差,观测量为三轴加速度误差和三轴地磁场误差。实验表明,将航姿系统固定在转台上产生0.6 g的向心加速度时,9态扩展卡尔曼滤波在对机动加速度进行稳定估计的同时,姿态精度保持在1度以内。

多类重力场模型的精度分析及联合确定GPS点正常高的方法

针对模型尺度参数ae对高程异常计算的影响进行数值分析;利用各模型位系数方差信息及区域高精度GPS/水准数据分析国际上新近公布的ITG-GRACE2010S、AIUB-GRACE03S、EIGEN-6C、GOCO02S、DIR_R3、TIM_R3、SPW_R2、gif48以及EGM2008多类重力场模型对应不同谱域的位系数精度;提出采用简单谱组合方法和加权谱组合方法来提高高程异常的计算精度。在试验区域,计算结果表明,各模型的中低阶(约135阶)位系数整体精度要优于EGM2008模型,gif48和EIGEN-6C模型的整体精度均优于同阶次的EGM2008模型,简单谱组合方法和加权谱组合方法均能提高区域模型高程异常精度。