GOCE Data Processing: The Spherical Cap Regularization Approach

Due to the sun-synchronous orbit of the satellite gravity gradiometry mission GOCE, the measurements will not be globally available. As a consequence, using a set of base functions with global support such as spherical harmonics, the matrix of normal equations tends to be ill-conditioned, leading to weakly determined low-order spherical harmonic coefficients. The corresponding geopotential strongly oscillates at the poles. Considering the special configuration of the GOCE mission, in order to stabilize the normal equations matrix, the Spherical Cap Regularization Approach (SCRA) has been developed. In this approach the geopotential function at the poles is predescribed by an analytical continuous function, which is defined solely in the spatially restricted polar regions. This function could either be based on an existing gravity field model or, alternatively, a low-degree gravity field solution which is adjusted from GOCE observations. Consequently the inversion process is stabilized. The feasibility of the SCRA is evaluated based on a numerical closed-loop simulation, using a realistic GOCE mission scenario. Compared with standard methods such as Kaula and Tikhonov regularization, the SCRA shows a considerably improved performance.     

GOCE Quick-Look Gravity Solution: Application of the Semianalytic Approach in the Case of Data Gaps and Non-Repeat Orbits

The satellite mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer), the first Core Mission of the Earth Explorer Programme funded by ESA (European Space Agency), is dedicated to the precise modelling of the Earth's gravity field, with its launch planned for 2006. The mathematical models for parameterizing the Earth's gravity field are based on a series expansion into spherical harmonics, yielding a huge number of unknown coefficients. Their computation leads to the solution of very large normal equation systems. An efficient way to handle these equation systems is the so-called semianalytic or lumped coefficients approach, which theoretically requires an uninterrupted, continuous time series of observations, recorded along an exact circular repeat orbit. In this paper the consequences of violating these conditions are analyzed. The effects of an interrupted observation stream onto the estimated spherical harmonic coefficients are demonstrated, and an iterative strategy, which reduces the negative influence depending on the characteristics of the data gaps, is proposed. Additionally, the impact of an imperfectly closing orbit (non-repeat orbit) on the gravity field model is analyzed, and a strategy to minimize the corresponding errors is presented. The applicability of the semianalytic approach also to a joint inversion of satellite-to-satellite tracking data in high-low mode (hl-SST) and satellite gravity gradiometry (SGG) observations is demonstrated, where the analysis of the former component is based on the energy conservation law. Several realistic case studies prove that the semianalytic approach is a feasible tool to generate quick-look gravity solutions, i.e. fast coefficient estimates using only partial data sets. This quick-look analysis shall be able to detect potential distortions of statistical significance (e.g. systematic errors) in the input data, and to give a fast feedback to the GOCE mission control.     

辽宁地区时变重力场源变化特征分析

地震重力分析通过研究时变重力场变化获取地球内部介质物性变化信息。采用贝叶斯重力网平差方法对辽宁地区2011—2014年共计四年7期的流动重力观测资料进行处理,对研究区重力观测网总体监测能力做出分析,选用欧拉反褶积对研究区的重力变化场源深度及空间分布特征进行反演和解释。通过反演计算发现,在2013年灯塔M_S5.1地震前沈阳—辽阳地区重力变化较大,形成高梯度带,反映该地区地震活动性增强;而欧拉反褶积计算反演地下介质场源深度在10～40 km,与地震实际发生深度相符,且集中在断裂带附近。本研究可为研究辽宁地区深部孕震环境特征提供一定参考。     

昆明地区重力异常数据处理

为了有效地划分出重力异常,利用小波多尺度分析和奇异值分解重建两种方法对昆明地区布格重力异常数据进行处理,实现了位场的分离.分离结果反映了不同深度和尺度地质体密度的不均匀性,从而较好地圈定了盆地,并揭示了几条主要断裂的空间延展特征.运用这两种方法对重力数据的处理都具有可行性.     

2014年山西平遥M_S4.1地震前后重力场变化特征分析

利用2012~2014年山西省绝对重力和流动重力观测资料,平差计算得到该时间段重力场变化值,分析区域重力场演化特征与2014年4月平遥MS4.1地震的关系。研究发现:(1)地震发生前研究区域出现了区域性的重力异常变化,在震中附近产生局部重力异常变化,形成了与活动断裂走向基本一致的重力变化高梯度带,这可能预示着与地震发生相关的构造运动或应力增强作用。(2)地震前后研究区重力场累积变化表现为应力场增强—持续增强—局部运动受阻超过临界值—能量突然释放,震后重力场变化恢复到平稳状态。     

2014年康定M_S6.3级地震前后的重力变化特征

利用四川地区2010年8月至2015年3月的流动重力观测数据作出重力场等值线变化图像,研究分析2014年康定MS6.3级地震前后的重力变化特征,结果表明:(1)三岔口地区长期保持区域性重力正值变化异常及重力变化高梯度带;(2)在康定MS6.3级地震前后三岔口地区的区域性重力正值变化异常及重力变化高梯度带变化并不是特别显著;(3)在康定MS6.3级地震前震区重力变化等值线出现四象限区域或者类似于四象限区域的分布特征;(4)康定MS6.3级地震发震位置处于0值线附近且0值线在此处发生明显的转折、畸变;(5)康定MS6.3级地震发生在重力减小后回调增加的过程中,震后重力继续回调增加。     

Derivation of the CHAMP-only global gravity field model TUG-CHAMP04 applying the energy integral approach

A global gravity field model TUG-CHAMP04, derived from CHAMP (CHAllenging Minisatellite Payload) satellite-to-satellite GPS tracking observations in the high-low mode (SST-hl) in combination with CHAMP accelerometry, is presented and described in detail in this paper. For this purpose the energy integral approach was applied to precise kinematic orbits and accelerometer data. The advantage of these kinds of orbits is that they are derived from purely geometrical information, hence no external gravity field information is used for the determination of the positions. The disadvantage of precise kinematic orbit information is, that no velocities are delivered and hence a procedure has to be elaborated to deduce the velocities from kinematic positions. This work is done in preparation for ESA’s GOCE (Gravity field and steady state Ocean Circulation Explorer) satellite mission (scheduled launch November 2006), aiming at a high precision and high-resolution gravity field model on a global scale. This paper concentrates on the CHAMP data processing, where, in contrast to the usual standard method (processing in the Earth fixed frame), an approach in the inertial frame is chosen. Focus is taken on the data preprocessing of both accelerometer and orbit data, emphasising on the correct treatment of data-gaps and outlier detection. Furthermore an arc-wise weighting strategy is introduced and the advantages/disadvantages of this approach are discussed. Finally, the TUG-CHAMP04 model, calculated from one year of CHAMP data is compared with the official CHAMP gravity field model EIGEN-3p and terrestrial data (GPS levelling data).     

宁夏地区重力场变化与2012年永宁MS4.6地震

利用宁夏地区2011—2013年期间多期重力观测资料,采用中国地震局推广的LG-ADJ程序对观测资料进行平差处理,给出地震前后3年的区域重力场变化和点值变化图像。结合构造活动分析,对2012永宁4.6级地震前后宁夏地区重力异常变化进行研究。认为：2012年永宁地震前,测区重力场出现较好的中期前兆性变化图像,临震前短期异常突增,重力场在空间和量值上均出现大幅的上升变化,重力点值震前趋势下降-转折上升-加速上升-发震的变化,是非常值得关注的前兆异常。