利用先验重力场模型对GOCE卫星重力梯度观测值进行校准分析

GOCE卫星任务搭载了高灵敏度的重力梯度仪,其观测值用于恢复高精度高分辨率的地球重力场。本文利用EIGEN-5C、EGM2008、GOTIM3、GGM03S高精度全球重力场模型,确定了GOCE引力梯度张量的对角分量观测值(Vxx、Vyy、Vzz)的校准参数,分析了比例因子的稳定性,并讨论了相同模型不同阶次、同阶次不同模型以及是否估计漂移参数对比例因子、偏差参数及校准观测值的影响。研究表明比例因子的稳定性在10-4的量级,利用250阶的EIGEN-5C模型和EGM2008模型校准得到观测值的差异小于10-4 E,远远小于观测误差,以1d为周期估计校准参数时,是否估计漂移对校准结果的影响达到0.4E。同时,校准前后观测值差异的频谱说明校准过程主要影响Vxx、Vyy、Vzz观测值的低频部分,即来自先验重力场模型的中低(150)阶次,考虑到GOCE引力梯度的观测频带,校准后的观测值可用于恢复中高频的重力场信号。     

卫星重力梯度观测数据的时变信号影响分析

系统地讨论了时变重力中潮汐信号与非潮汐信号对GOCE卫星重力梯度观测数据的影响。结果表明:(1)时变改正的量级为0.1 mE,比GOCE卫星设计精度(3.2 mE)低,但其为有色噪声,在数据预处理中必须剔除;(2)潮汐影响(0.1 mE)比非潮汐影响(0.01 mE)要高一个量级,决定着时变重力改正的精度。将本文计算结果与GOCE官方公布结果进行对比,二者具有较好一致性,验证了本文计算方法及结果的有效性。     

GOCE卫星重力梯度数据反演重力场的滤波器设计与比较分析

地球重力场和海洋环流探测（gravity field and steady-state ocean circulation explorer,GOCE）卫星重力梯度数据有色噪声和低频系统误差的滤波处理是反演高精度地球重力场的一个关键问题。针对GOCE卫星重力梯度数据的滤波处理,基于移动平均（moving average,MA）方法和CPR(circle per revolution)经验参数方法设计了两类低频系统误差滤波器,并分别将这两类滤波器与基于自回归移动平均（auto-regressive and moving average,ARMA）模型设计的有色噪声滤波器组合起来形成级联滤波器。为了分析滤波器处理的实际效果,基于空域最小二乘法采用70 d的GOCE观测数据,并联合重力恢复与气候实验（gravity recovery and climate experiment,GRACE）数据分别反演了224阶次的重力场模型GOGR-MA(MA+ARMA级联滤波)和GOGR-CPR(CPR+ARMA级联滤波)。将反演模型与采用同期数据求解的第一代GOCE系列模型及GOCE和GRACE联合模...     

联合GOCE卫星轨道和重力梯度数据严密求解重力场的模拟研究

论述了联合卫星轨道和重力梯度数据严密求解重力场的方法及数据处理方案,研究了GOCE重力场反演中有色噪声的AR去相关滤波、病态法方程的Kaula正则化和观测值最优加权的方差分量估计等关键问题。模拟结果表明:①极空白问题会降低法方程求解的稳定性,导致低次位系数的求解精度较低,而Kaula正则化可有效用于GOCE病态法方程的求解,并得到合理稳定的解;②重力梯度有色噪声会降低GOCE重力场求解的整体精度,特别是对低阶位系数的影响最为明显,而AR去相关滤波法可有效处理有色噪声,但解算结果仍含有低频误差;③方差分量估计可有效确定SST和SGG两类观测值的最优权比,并且有色噪声造成的低频误差经过联合求解后得到了抑制;④利用30d、5s采样的GOCE模拟数据恢复200阶次的重力场模型,其大地水准面和重力异常精度在纬度±83°范围内分别为±3.81cm和±1.056mGal。     

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

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

GOCE卫星重力梯度测量数据的预处理

GOCE(gravity field and steady-state ocean circulation explorer)计划的主要科学目标是以70 km空间分辨率1、mGal重力异常和1~2 cm大地水准面的精度测定全球静态地球重力场,卫星重力梯度测量数据的预处理是实现这一预期科学目标的重要任务之一。讨论了重力梯度测量数据的预处理方案、时变重力场信号改正、粗差探测和外部校准方法,为进一步开展GOCE卫星重力梯度测量数据的预处理研究提供参考和具体建议。     

利用先验重力场模型求定GOCE卫星重力梯度仪校准参数

重力梯度仪校准参数的确定是GOCE重力梯度观测数据处理的关键环节。本文对GOCE卫星重力梯度观测值中的时变信号与粗差进行了分析,利用高精度全球重力场模型,确定了GOCE重力梯度观测值各分量的尺度因子与偏差,并对校准结果进行了精度评定。结果表明,在测量带宽内,海潮对重力梯度观测值影响在mE量级,与重力梯度仪的精度水平相当,陆地水等非潮汐重力场时变信号略小于海潮,量级约为10^-4E;各分量重力梯度观测值的粗差比例均大于0.2%;除EGM96模型外的其他模型对GOCE重力梯度仪进行校准后,V_xx、V_yy、V_zz、V_yz分量上尺度因子的稳定性均在10^-4量级,V_xz分量能达到10^-5量级,V_xy分量为10^-2量级,这与梯度观测值各分量的精度水平一致。     

天文潮汐对卫星重力梯度观测数据的影响分析

测定地球重力场,确定高分辨率的静态地球重力场模型,是大地测量学的主要任务之一.重力场的影响主要分为潮汐部分和非潮汐部分,天文潮汐在潮汐部分中属于直接引力效应,对重力场的影响是不可忽略的.本文以一个月的星历数据为基础,分析了天文潮汐对GOCE卫星重力梯度观测数据的影响,并统计了最大值和最小值;研究了天文潮汐对地球上单点重力梯度数据的影响特征;计算了各行星对卫星重力梯度数据影响量级.研究结果表明:天文潮汐对卫星重力梯度数据的影响量级处于0.1mE,比GOCE卫星设计精度低一个量级,但是它具有周期性,属于有色噪声,因此在卫星重力梯度数据预处理中需要扣除;天文潮汐对卫星重力梯度数据各分量的影响不同,其中对角线分量Vxx,Vyy和Vzz要比其他分量略大;月球和太阳对卫星重力梯度数据的影响最大,在所有星体中占据主导地位.     

利用卫星重力梯度数据确定地球重力场的单层位模型

本文提出了处理卫星重力梯度数据以确定高分辩力重力场模型的单层位法并对其中的独立估计法进行了误差分析，数字结果显示：当卫星高度为200km，卫星数据网格宽度为15′，卫星重力梯度数据的精度为2×10~(-3)E时，利用独立估计法可得到分辩力为1°×1°（100km）的全球重力场模型，其重力异常精度小于1（mgal）；若卫星高度降至160km，卫星重力梯度数据的精度达到3×10~(-4)E，则获得的重力场模型的分辩力可提高到0.5°×0.5°（50km），其重力异常精度仍小于1（mgal）。     

卫星姿态误差对GOCE重力梯度精度的影响评估

在地球重力场和海洋环流探测卫星GOCE(Gravity field and Ocean Circulation Explorer)的观测数据中,其主要的观测量重力梯度数据不仅与搭载的6个加速度计的测量值有关,而且还与卫星自身的自转角速度存在着二次函数关系.由于加速度计测量频段的限制,这样就导致了卫星姿态的低频误差混入到...     

Efficient gravity field recovery from GOCE gravity gradient observations

 An efficient algorithm is proposed for gravity field recovery from Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) satellite gravity gradient observations. The mathematical model is formulated in the time domain, which allows the inclusion of realistic observational noise models. The algorithm combines the iterative solution of the normal equations, using a Richardson-type iteration scheme, with the fast computation of the right-hand side of the normal equations in each iteration step by a suitable approximation of the design matrix. The convergence of the iteration is investigated, error estimates are provided, and the unbiasedness of the method is proved. It is also shown that the method does not converge to the solution of the normal equations. The performance of the approach for white noise and coloured noise is demonstrated along a simulated GOCE orbit up to spherical harmonic degree and order 180. The results also indicate that the approximation error may be neglected. Received: 30 November 1999 / Accepted: 31 May 2000     

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方法拥有快速高精度反演卫星重力场模型的优势,可以在重力梯度卫星的设计、误差分析及在轨快速评估等方面得到充分应用。     

A parametric study on the impact of satellite attitude errors on GOCE gravity field recovery

In the recent design of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite mission, the gravity gradients are defined in the gradiometer reference frame (GRF), which deviates from the actual flight direction (local orbit reference frame, LORF) by up to 3–4°. The main objective of this paper is to investigate the effect of uncertainties in the knowledge of the gradiometer orientation due to attitude reconstitution errors on the gravity field solution. In the framework of several numerical simulations, which are based on a realistic mission configuration, different scenarios are investigated, to provide the accuracy requirements of the orientation information. It turns out that orientation errors have to be seriously considered, because they may represent a significant error component of the gravity field solution. While in a realistic mission scenario (colored gradiometer noise) the gravity field solutions are quite insensitive to small orientation biases, random noise applied to the attitude information can have a considerable impact on the accuracy of the resolved gravity field models.     

卫星重力梯度数据重力异常的精度分析

针对GOCE卫星确定的地球重力场模型精度的不确定性,对比分析GOCE位模型与多个不同重力场模型确定的重力异常,并将其分别与船测重力数据、南极航空重力数据、北极重力数据以及美国和中国台湾地面重力数据比较研究。结果表明:GOCE位模型的内符合精度最高,与地面重力观测数据符合最优;与船测以及航空重力测量符合相对较差、精度较低。研究表明,在一定精度前提下,GOCE卫星确定的重力数据可用于无人区,从而提高重力观测数据的覆盖率。     

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

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

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     

The gravity field model IGGT_R1 based on the second invariant of the GOCE gravitational gradient tensor

Based on tensor theory, three invariants of the gravitational gradient tensor (IGGT) are independent of the gradiometer reference frame (GRF). Compared to traditional methods for calculation of gravity field models based on the gravity field and steady-state ocean circulation explorer (GOCE) data, which are affected by errors in the attitude indicator, using IGGT and least squares method avoids the problem of inaccurate rotation matrices. The IGGT approach as studied in this paper is a quadratic function of the gravity field model’s spherical harmonic coefficients. The linearized observation equations for the least squares method are obtained using a Taylor expansion, and the weighting equation is derived using the law of error propagation. We also investigate the linearization errors using existing gravity field models and find that this error can be ignored since the used a-priori model EIGEN-5C is sufficiently accurate. One problem when using this approach is that it needs all six independent gravitational gradients (GGs), but the components \(V_{xy}\) and \(V_{yz}\) of GOCE are worse due to the non-sensitive axes of the GOCE gradiometer. Therefore, we use synthetic GGs for both inaccurate gravitational gradient components derived from the a-priori gravity field model EIGEN-5C. Another problem is that the GOCE GGs are measured in a band-limited manner. Therefore, a forward and backward finite impulse response band-pass filter is applied to the data, which can also eliminate filter caused phase change. The spherical cap regularization approach (SCRA) and the Kaula rule are then applied to solve the polar gap problem caused by GOCE’s inclination of \(96.7^{\circ }\). With the techniques described above, a degree/order 240 gravity field model called IGGT_R1 is computed. Since the synthetic components of \(V_{xy}\) and \(V_{yz}\) are not band-pass filtered, the signals outside the measurement bandwidth are replaced by the a-priori model EIGEN-5C. Therefore, this model is practically a combined gravity field model which contains GOCE GGs signals and long wavelength signals from the a-priori model EIGEN-5C. Finally, IGGT_R1’s accuracy is evaluated by comparison with other gravity field models in terms of difference degree amplitudes, the geostrophic velocity in the Agulhas current area, gravity anomaly differences as well as by comparison to GNSS/leveling data.     

GOCE地球重力场模型在青藏地区的评价及分析

2009年GOCE卫星升空以后,卫星重力梯度数据参与解算的GOCE系列重力场模型已有多家研究机构相继公布。本文分别采用青藏地区的GPS/水准和重力异常实测数据对GOCE重力场模型进行了外部测试,并在重力异常验证过程中引入了一种新的滤波方法,验证结果表明在青藏地区GOCE重力场模型相比其它系列模型的优势在于中波段。同时,探讨了GOCE重力场模型与其他系列模型在青藏地区主要差异值的空间分布以及首次利用统计分析方法找出模型之间主要差异值的阶次分布,得出如下结论：模型之间的较大差异值在空间水平方向上主要分布在喜马拉雅山脉、天山等地形起伏较大的区域,在垂直方向上主要集中在岩石圈。     

GRACE和GOCE卫星监测的青藏高原重力场特征

文章阐述了对青藏高原重力场进行研究的意义,并进一步利用重力卫星GRACE和GOCE的数据对该区域的重力场特征进行了描述.通过对该区域的重力异常、径向引力梯度的计算和分析,可以得出:在青藏高原的西部,有明显的3条重力异常区,这与当地的地形有关,也与断层的位置有关;引力梯度比重力异常具有更高的空间分辨率;重力变化剧烈的区域与梯度的异常区有一定的对应关系,同时也是地球动力活动变化剧烈的区域.