The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

A data-driven approach to local gravity field modelling using spherical radial basis functions

We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.     

广州市地面塌陷分布特征与人为致灾因子分析

广州市地面塌陷频繁发生,近10年来呈波动式上升趋势。其塌陷类型分为岩溶塌陷和工程地面塌陷两种,岩溶塌陷主要分布在广花盆地的花都区、白云区,工程地面塌陷主要分布在主城区。地面塌陷孕灾环境复杂、致灾因子多样,承灾体脆弱,灾害后果严重,人为因素是地面塌陷的主要诱发因素。岩溶塌陷主要是由于过量抽取地下水或矿山疏干排水、地下采空、暴雨触发所致;工程塌陷主要是人类工程行为所致,其主要致灾因子包括排水疏干与突水（突泥）作用、人工加载、人工振动、人工开挖桩基、地表渗水、地铁等地下工程盾构掘进等。     

广东省五华县地质灾害形成特征及防治对策

五华县主要地质灾害类型有滑坡、崩塌、地面塌陷、水土流失等。其中以滑坡、崩塌为主,多分布在东南、南部花岗岩区、北部花岗岩风化土区等广大中低山及丘陵区,具有点多面广,灾害点个体规模小,稳定性差,活动频繁,地质灾害发育呈明显的地域性与季节性分布等特点。五华县地质灾害的形成与发生是多种致灾因素相互作用的结果。地层岩性是其形成的内在要素,它在一定程度上决定着地质灾害的发育程度与类型;地形地貌与植被是地质灾害形成的外在条件,它制约着崩、滑、塌等致灾地质作用的形成;大气降雨是地质灾害形成与发生的激发因素,决定着地质灾害发生的速度和时间;人类工程活动是影响地质灾害形成与发生的最主要、最直接的因素。对地质灾害的防治应采用避让、预防、监测及治理措施,做到避让与治理结合,以群测群防为基本手段,点状灾害以工程治理与生物防治为主;面状灾害以生物防治为主;采用点、面结合综合治理的方法。     

隧道围岩强度不均地段塌方成因及其处理方法

由于地质构造作用以及风化作用,工程场地岩体中多存在围岩强度不均地段.在隧道施工过程中就经常遇到软硬岩石交汇而形成的地质条件复杂地段;给隧道的施工带来了很多不利情况,并且也是导致隧道塌方的原因之一.2006年1月,承德韩郭线二级公路工程,喇嘛梁隧道软硬岩石交汇地段由于地下水作用导致围岩风化程度不均,产生软弱结构面并造成了大面积塌方.基于现场的追踪调查与考察,本文详细描述了利用拱顶架设工字钢梁并配合超前注浆小导管的支护方法来处理这类围岩.     

地下矿产开采对天然气管道工程的安全影响评价

西气东输天然气管道工程途经兖州市,穿越杨庄、杨村、兴隆庄等煤矿矿区,压覆煤层埋深变化较大。该文采用预测法对各煤矿在不同开采条件下分煤层对管道运营安全的危害程度进行了评价分析。通过计算表明,煤矿开采埋深越浅,煤层越厚,对管道危害越严重。兴隆庄煤矿3煤层、杨村煤矿煤16t、煤17开采对管道影响强烈,杨庄煤矿煤16t、煤17开采对管道影响中等,兴隆庄煤矿煤16t、煤17开采对管道运营影响小。     

邹平火山岩盆地铜矿地质成因地球物理特征探讨

邹平铜矿处于齐河-广饶深大断裂带南部的邹平火山岩盆地中,形成于破火山口火山通道充填的石英正长闪长岩岩颈中央上部,包括伟晶岩型铜矿和细脉浸染状斑铜矿床两种类型。前者矿体较小,但品位高;后者品位较低,但规模中等。含矿石英正长闪长岩等密度小、磁性弱,故在火山岩系中呈现高背景重力场上的重力低和杂乱高磁场背景中的低负异常,即“重磁同低”,且高极化。重磁同低异常区和高极化率异常带,是本区寻找铜矿的有利部位。     

联合多种测高数据确定中国海及其邻域1．5′×1．5′重力异常

联合多种测高资料和Geosat/GM波形重构数据,基于EIGEN_CG01C重力场模型,采用沿轨迹加权最小二乘方法和逆Vening-Meinesz公式,确定了中国海及其邻域1.5′×1.5′重力异常。将计算结果与最新船测资料进行了比较,标准差为3.37×10-5m.s-2。     

The Glueckstadt Graben of the North-German Basin: new insights into the structure from 3D and 2D gravity analyses

The structure of the Glueckstadt Graben has been investigated by use of 3D gravity backstripping technique and by 2D gravity and magnetic modelling. Subtracting the gravity effects of the Meso-Cenozoic sediments together with Permian salt reveals a positive residual anomaly within the Glueckstadt Graben. This anomaly includes two local maxima over the Westholstein and Eastholstein Troughs. The 2D gravity models point to the presence of a high-density body within the lower crust of the Glueckstadt Graben. In addition, the results of 2D magnetic modelling indicate that the central part of the high-density body is overlain by an area with high susceptibility. Most probable, the formation of this high-density body is a result of complex poly-phase tectonic history of the study area. Finally, the results of gravity modelling indicate that Permian salt is not homogeneous. 3D gravity analysis and, especially, 2D gravity modelling have distinguished the differences in degree of salt saturation in salt-rich bodies, and elucidate the proportion of Rotliegend salt.