GOCE实测数据反演高阶重力场模型的Torus方法

不同于当前广泛使用的空域法、时域法、直接解法,本文尝试采用Torus方法处理GOCE实测数据,利用71 d的GOCE卫星引力梯度数据反演了200阶次GOCE地球重力场模型,实现了对参考模型的精化。首先,采用Butterworth零相移滤波方法加移去—恢复技术,处理引力梯度观测值中的有色噪声,并利用泰勒级数展开和Kriging方法对GOCE卫星引力梯度数据进行归算和格网化,计算得到了名义轨道上格网点处的引力梯度数据。然后,利用2D-FFT技术和块对角最小二乘方法处理名义轨道上数据,获得了200阶次的GOCE地球重力场模型GOCE_Torus。利用中国和美国的GPS/水准数据进行外部检核结果说明,GOCE_Torus与ESA发布的同期模型的精度相当;GOCE_Torus模型与200阶次的EGM2008模型相比,在美国区域精度相当,但在中国区域精度提高了4.6 cm,这充分体现了GOCE卫星观测数据对地面重力稀疏区的贡献。Torus方法拥有快速高精度反演卫星重力场模型的优势,可以在重力梯度卫星的设计、误差分析及在轨快速评估等方面得到充分应用。     

国际卫星重力梯度测量计划研究进展

本文首先阐述了重力梯度测量原理、从20世纪初到21世纪初重力梯度仪的研究历程、卫星重力梯度仪(静电悬浮重力梯度仪、超导重力梯度仪和量子重力梯度仪)的技术特征以及卫星重力梯度测量的特点;其次,介绍了基于卫星重力梯度技术恢复250阶GOCE地球重力场以及论证首先开展一维径向重力梯度仪的研制进而恢复高精度和高空间解析度中高频地球重力场可行性方面的研究进展;最后,建议我国尽早开展基于时空域混合法解算中高频地球重力场和卫星重力梯度测量系统误差分析的预先研究。     

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.     

GOCE卫星引力梯度数据滤波方法研究

针对GOCE卫星引力梯度观测值中低精度分量和低频有色噪声的处理策略问题,该文采用模型模拟值代替低精度分量Vxy和Vyz,以减弱低精度分量在坐标系转换中对高精度分量的影响。深入分析比较了多种滤波方法处理GOCE卫星引力梯度观测值中有色噪声的效果,提出采用Butterworth零相移滤波方法加移去-恢复技术的思路,实测数据的处理效果验证了该方法的有效性。     

Spherical Tikhonov regularization wavelets in satellite gravity gradiometry with random noise

 A special class of regularization methods for satellite gravity gradiometry based on Tikhonov spherical regularization wavelets is considered, with particular emphasis on the case of data blurred by random noise. A convergence rate is proved for the regularized solution, and a method is discussed for choosing the regularization level a posteriori from the gradiometer data. Received: 23 March 2000 / Accepted: 20 September 2000     

Spherical harmonic synthesis and least squares computations in satellite gravity gradiometry

The computational requirements in the simulations of geopotential estimation from satellite gravity gradiometry are discussed. Fast algorithms for spherical harmonic synthesis and least squares accumulation on a vectorizing supercomputers are presented. Using these methods, in a test case estimation of 2595 coefficients of a degree and order 50 gravity field, sustained program execution speeds of 275 Mflops (87 % peak machine speed) on a single processor of a CRAY Y-MP were achieved, with spherical harmonics computation accounting for less than 1 % of total cost. From the results, it appears that brute-force estimation of a degree and order 180 field would require 537 Million Words of memory and 85 hours of CPU time, assuming mission duration of 1 month, and execution speed of 1 Gflops. Both memory size and execution speed requirements are within the capabilities of modern multi-processor supercomputers.     

Assessment of three numerical solution strategies for gravity field recovery from GOCE satellite gravity gradiometry implemented on a parallel platform

 The recovery of a full set of gravity field parameters from satellite gravity gradiometry (SGG) is a huge numerical and computational task. In practice, parallel computing has to be applied to estimate the more than 90 000 harmonic coefficients parameterizing the Earth's gravity field up to a maximum spherical harmonic degree of 300. Three independent solution strategies (preconditioned conjugate gradient method, semi-analytic approach, and distributed non-approximative adjustment), which are based on different concepts, are assessed and compared both theoretically and on the basis of a realistic-as-possible numerical simulation regarding the accuracy of the results, as well as the computational effort. Special concern is given to the correct treatment of the coloured noise characteristics of the gradiometer. The numerical simulations show that the three methods deliver nearly identical results—even in the case of large data gaps in the observation time series. The newly proposed distributed non-approximative adjustment approach, which is the only one of the three methods that solves the inverse problem in a strict sense, also turns out to be a feasible method for practical applications. Received: 17 December 2001 / Accepted: 17 July 2002 Acknowledgments. We would like to thank Prof. W.-D. Schuh, Institute of Theoretical Geodesy, University of Bonn, for providing us with the serial version of the PCGMA algorithm, which forms the basis for the parallel PCGMA package developed at our institute. This study was partially performed in the course of the GOCE project `From E?tv?s to mGal+', funded by the European Space Agency (ESA) under contract No. 14287/00/NL/DC. Correspondence to: R. Pail     

Improvement of the geoid in local areas by satellite gradiometry

Summary Satellite gradiometry is studied as a means to improve the geoid in local areas from a limited data coverage. Least-squares collocation is used for this purpose because it allows to combine heterogeneous data in a consistent way and to estimate the integrated effect of the attenuated spectrum. In this way accuracy studies can be performed in a general and reliable manner. It is shown that only three second-order gradients contribute significantly to the estimation of the geoidal undulations and that it is sufficient to have gradiometer data in a 5°×5° area around the estimation point. The accuracy of the geoid determination is strongly dependent on the degree and order of the reference field used. An accuracy of about ±1 m can be achieved with a reference field of (12, 12). There is an optimal satellite altitude for each reference field and this altitude may be higher than 300 km for a field of low degree and order. The influence of measuring errors is discussed and it is shown that only gradiometer data with accuracies better than ±0.05 E will give a significant improvement of the geoid. Finally, some results on the combination of satellite gradiometry and terrestrial gravity measurements are given. The proposed method seems to be well suited for local geoid determinations down to the meter range. It is especially interesting for unsurveyed and difficult areas because no terrestrial measurements are necessary. Furthermore, it has the practical advantage that only a local data coverage is needed.     

重力梯度之地形改正方法与研究

简要介绍地形起伏对地面点重力梯度的影响,推导棱柱法计算重力梯度改正的严格积分公式并和积分核累加法对比验证了其正确性;将积分区域划分为中央区和远区,中央区采用严格的棱柱法进行计算,远区采用线性积分的方式逼近真实值,简化了积分过程;根据重力梯度应用的要求针对不同类型地形设定了中央区划分准则。实验表明:根据文中提出的方法和分区准则,相比于传统的积分法和棱柱法,计算速度得到提高,且精度要求相比于积分核累加法优于0.01E。     

Computation of spherical harmonic coefficients from gravity gradiometry data to be acquired by the GOCE satellite: regularization issues

The issue of optimal regularization is investigated in the context of the processing of satellite gravity gradiometry (SGG) data that will be acquired by the GOCE (Gravity Field and Steady-State Ocean Circulation Explorer) satellite. These data are considered as the input for determination of the Earths gravity field in the form of a series of spherical harmonics. Exploitation of a recently developed fast processing algorithm allowed a very realistic setup of the numerical experiments to be specified, in particular: a non-repeat orbit; 1-s sampling rate; half-year duration of data series; and maximum degree and order set to 300. The first goal of the study is to compare different regularization techniques (regularization matrices). The conclusion is that the first-order Tikhonov regularization matrix (the elements are practically proportional to the degree squared) and the Kaula regularization matrix (the elements are proportional to the fourth power of the degree) are somewhat superior to other regularization techniques. The second goal is to assess the generalized cross-validation method for the selection of the regularization parameter. The inference is that the regularization parameter found this way is very reasonable. The time expenditure required by the generalized cross-validation method remains modest even when a half-year set of SGG data is considered. The numerical study also allows conclusions to be drawn regarding the quality of the Earths gravity field model that can be obtained from the GOCE SGG data. In particular, it is shown that the cumulative geoid height error between degrees 31 and 200 will not exceed 1 cm. AcknowledgmentsThe authors thank Dr. E. Schrama for valuable discussions and for computing the orbit used to generate the long data set. They are also grateful to Prof. Tscherning and two anonymous reviewers for numerous valuable remarks and suggestions. The orbit to generate the short data set was kindly provided by J. van den IJssel. Computing resources were provided by Stichting Nationale Computerfaciliteiten (NCF), grant SG-027.     

月球重力场模型研究进展和我国将来月球卫星重力梯度计划实施

本文介绍了基于国际探月观测数据建立的月球重力场模型:8×4、15×8、13×13、5×5、7×7、16×16-1/2/3、Lun60d、GLGM-1/2、LP75D/G、LP100K/J、LP165P、LP150Q和SGM90d;通过对比SST-HL/LL-Doppler-VLBI和SST-HL/SGG-Doppler-VLBI跟踪观测模式的优缺点,建议我国将来首期月球卫星重力测量计划采用SST-HL/SGG-Doppler-VLBI较优;其次,通过对比静电悬浮、超导和量子卫星重力梯度仪的优缺点,建议我国将来首期月球卫星重力梯度计划采用静电悬浮重力梯度仪;并建议我国将来首颗月球重力梯度卫星的轨道高度(50～100 km)选择在已有月球探测卫星的测量盲区,轨道倾角(90°±3°)设计为有利于月球卫星观测数据全球覆盖的近极轨模式。     

联合卫星重力和卫星测高数据确定稳态海洋动力地形

新一代卫星重力探测任务GRACE大大提高了地球重力场模型中长波分量的精度,使得联合卫星测高平均海面分离更精细稳态海洋动力地形成为可能。本文利用T/P(1994年~2003年)和JASON-1(2003年~2005年)卫星测高数据确定了全球30′×30′平均海面高;基于重力场模型WHU-GM-05,计算得到对应于海面高的GRACE海洋大地水准面格网值;利用“移去-恢复法”和高斯滤波求得全球稳态海面地形。与EGM96、R io05、ECCO和GGM02模型进行比较,检验结果表明GRACE任务有效的改善了海洋大地水准面的精度,使得稳态海洋动力地形能够呈现更多细部。     

GOCE卫星重力梯度观测值的时变重力场变化改正及影响分析

GOCE卫星重力梯度观测值为高阶静态重力场反演提供了重要的数据支撑,但其在使用前需考虑扣除时变重力场变化的影响.本文研究了GOCE卫星重力梯度观测值的时变重力场变化改正方法,更新了ESA标准和背景模型,以更好地扣除时变重力场变化的影响,自主实现了由GOCE卫星Level1b重力梯度数据直接进行重力场反演.本文通过3种时...     

联合卫星重力和GNSS观测反演地震断层参数的研究EI北大核心CSCD

地震的发生会造成巨大的破坏和损失,大地测量技术能观测地震形变和反演地震断层的错动情况,对于减少地震灾害具有非常重要的现实意义。GNSS技术已广泛用于震源参数的反演工作,但对于发生在海域的地震,由于GNSS站点往往位于陆地一侧,导致对断层的约束 能力有限。卫星重力技术因其全天候、全球覆盖、连续性,以及不受地域地形限制等诸多特点,可弥补海洋一侧常规观测不足。联合GNSS和GRACE两种观测手段反演震源机制,可进一步提高反演海域地震的震源参数。为此,论文基于GRACE和GNSS观测的同震和震后形变,联合反演了断层参数和地幔黏滞性系数等,主要成果如下。     

A comparison of different mass elements for use in gravity gradiometry

Topographic and isostatic mass anomalies affect the external gravity field of the Earth. Therefore, these effects also exist in the gravity gradients observed, e.g., by the satellite gravity gradiometry mission GOCE (Gravity and Steady-State Ocean Circulation Experiment). The downward continuation of the gravitational signals is rather difficult because of the high-frequency behaviour of the combined topographic and isostatic effects. Thus, it is preferable to smooth the gravity field by some topographic-isostatic reduction. In this paper the focus is on the modelling of masses in the space domain, which can be subdivided into different mass elements and evaluated with analytical, semi-analytical and numerical methods. Five alternative mass elements are reviewed and discussed: the tesseroid, the point mass, the prism, the mass layer and the mass line. The formulae for the potential, the attraction components and the Marussi tensor of second-order potential derivatives are provided. The formulae for different mass elements and computation methods are checked by assuming a synthetic topography of constant height over a spherical cap and the position of the computation point on the polar axis. For this special situation an exact analytical solution for the tesseroid exists and a comparison between the analytical solution of a spherical cap and the modelling of different mass elements is possible. A comparison of the computation times shows that modelling by tesseroids with different methods produces the most accurate results in an acceptable computation time. As a numerical example, the Marussi tensor of the topographic effect is computed globally using tesseroids calculated by Gauss–Legendre cubature (3D) on the basis of a digital height model. The order of magnitude in the radial-radial component is about  ± 8 E.U. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

On the feasibility of using satellite gravity observations for detecting large-scale solid mass transfer events

The main focus of this paper is to assess the feasibility of utilizing dedicated satellite gravity missions in order to detect large-scale solid mass transfer events (e.g. landslides). Specifically, a sensitivity analysis of Gravity Recovery and Climate Experiment (GRACE) gravity field solutions in conjunction with simulated case studies is employed to predict gravity changes due to past subaerial and submarine mass transfer events, namely the Agulhas slump in southeastern Africa and the Heart Mountain Landslide in northwestern Wyoming. The detectability of these events is evaluated by taking into account the expected noise level in the GRACE gravity field solutions and simulating their impact on the gravity field through forward modelling of the mass transfer. The spectral content of the estimated gravity changes induced by a simulated large-scale landslide event is estimated for the known spatial resolution of the GRACE observations using wavelet multiresolution analysis. The results indicate that both the Agulhas slump and the Heart Mountain Landslide could have been detected by GRACE, resulting in \({\vert }0.4{\vert }\) and \({\vert }0.18{\vert }\) mGal change on GRACE solutions, respectively. The suggested methodology is further extended to the case studies of the submarine landslide in Tohoku, Japan, and the Grand Banks landslide in Newfoundland, Canada. The detectability of these events using GRACE solutions is assessed through their impact on the gravity field.     

Polar motion determined by Doppler satellite observations

Doppler observations of Navy Navigation Satellites have been used to comput pole positions on a daily basis since 1969. Limited computations have been performed using data on file for the period 1964–1969. Results of recent computations give a standard error in pole position based on 48 hours of Doppler observations of 7 cm. However, effects of errors in the orbit due to uncertainties in the gravity field prevent the attainment of this precision; the standard deviation of pole position for this time span is 60 cm, giving a standard error for a five day mean based on observations of two satellites of 25 cm.     

利用卫星重力探测南极冰盖质量变化

针对不同GIA模型在南极地区的改正差异较大,分析了采用不同GIA模型分析得到的南极冰盖融化误差值,研究了南极冰盖质量变化加速度。研究结果表明,扣除IJ05模型影响后,整个南极冰盖质量以(-81.5±4.2)Gt/a的速度和(-12.3±6.5)Gt/a2的加速度融化,对海平面上升贡献分别为(0.23±0.01)mm/a和(0.03±0.02)mm/a2。西南极处于加速融化状态,速度由2003—2008年间的-37.7Gt/a增加至2009—2013年的-156.0Gt/a。东南极冰盖在2009年前处于稳定,2009年后逐渐积累,且积累速度不断增加。2003-2013期间内东南极,西南极冰盖质量变化速度及加速度分别为(28±13.4)Gt/a、(9.8±2.8)Gt/a2、-(108.1±3.5)Gt/a和-(21.1±2.9)Gt/a2,南极冰盖融化主要来自西南极。南极冰盖周年变化影响较强,每年10月左右其质量变化振幅达到最大,半周年变化影响最大在西南极,S2潮波振幅在Ronne冰架附近较大,振幅最大值达到25mm。