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

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

GOCE gradiometer: estimation of biases and scale factors of all six individual accelerometers by precise orbit determination

A method has been implemented and tested for estimating bias and scale factor parameters for all six individual accelerometers that will fly on-board of GOCE and together form the so-called gradiometer. The method is based on inclusion of the individual accelerometer observations in precise orbit determinations, opposed to the baseline method where so-called common-mode accelerometer observations are used. The method was tested using simulated data from a detailed GOCE system simulator. It was found that the observations taken by individual accelerometers need to be corrected for (1) local satellite gravity gradient (SGG), and (2) rotational terms caused by centrifugal and angular accelerations, due to the fact that they are not located in the satellite’s center of mass. For these corrections, use is made of a reference gravity field model. In addition, the rotational terms are derived from on-board star tracker observations. With a perfect a priori gravity field model and with the estimation of not only accelerometer biases but also accelerometer drifts, scale factors can be determined with an accuracy and stability better than 0.01 for two of the three axes of each accelerometer, the exception being the axis pointing along the long axis of the satellite (more or less coinciding with the flight direction) for which the scale factor estimates are unreliable. This axis coincides with the axis of drag-free control, which results in a small variance of the signal to be calibrated and thus an inaccurate determination of its scale factor in the presence of relatively large (colored) accelerometer observation errors. In the presence of gravity field model errors, it was found that still an accuracy and stability of about 0.015 can be obtained for the accelerometer scale factors by simultaneously estimating empirical accelerations.     

Tikhonov正则化方法在GOCE重力场求解中的模拟研究

本文在阐述Tikhonov正则化方法基本原理的基础上,给出了四类可用于重力场解算的正则化矩阵(零次、一次、二次和Kaula),以及用于确定正则化参数的L曲线法和GCV方法的数学模型。基于SA方法利用模拟数据分析讨论了零次、一次以及Kaula正则化矩阵应用于GOCE全球重力场模型确定的有效性,并由Kaula正则化矩阵分析了L曲线法和GCV方法确定正则化参数的可行性。数值结果表明三类正则化矩阵获得的最优解(以大地水准面MSE最小为准则确定)的精度水平相近,关键在于相应正则化参数的确定,数值结果同时说明了GCV方法和L曲线法可用于确定正则化参数,且前者较后者具有更好的稳定性。     

融合GOCE和GRACE卫星数据的无约束重力场模型Tongji-GOGR2019S

本文在法方程层面融合GOCE卫星的V_xx、V_yy、V_zz和V_xz重力梯度分量观测数据和GRACE卫星观测数据,采用直接法解算了220阶次的重力场模型Tongji-GOGR2019S.首先利用ⅡR带通滤波器在5~41 mHz的重力梯度带宽范围内对约24个月的GOCE重力梯度观测方程进行无相移滤波处理,并组成解算220阶次重力场模型的法方程,各梯度分量根据相对于参考模型统计精度进行定权;然后与13.5 a GRACE数据建立的180阶次Tongji-Grace02s重力场模型的法方程进行叠加,解算了220阶次的无约束纯卫星重力场模型Tongji-GOGR2019S.利用EIGEN-6C4重力场模型、GNSS/水准数据、DTU15重力异常数据以及欧洲区域似大地水准面模型EGG2015等数据对Tongji-GOGR2019S模型精度进行全面的检核评定,结果表明：引入GOCE卫星梯度数据后,高于72阶的位系数精度优于Tongji-Grace02s模型,Tongji-GOGR2019S模型的整体精度接近同阶次的DIR-R6等GOCE卫星第6代模型.

     

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

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

重力及重力梯度联合反演青藏高原及邻区岩石圈三维密度结构

本文利用GOCE L2观测重力梯度的五个独立分量（Txx,Tzz,Txy,Txz,Tyz）,联合EGM2008地球重力场模型计算垂直重力,反演计算了青藏高原及邻区0~120 km深度岩石圈三维密度结构.将经过低阶项改正、地形效应改正、沉积层界面起伏效应改正得到的剩余重力及重力梯度异常值作为观测值,以改正剩余量归一化权重作为观测权重,基于Tikhonov正则化理论建立反演目标函数.反演过程中,利用地震层析S波速度转换密度作为初始约束,通过非等权最小二乘迭代法计算得到最终反演密度.反演结果表明：（1）40 km深度,青藏高原内部为中地壳,表现为低密度,邻区为中下地壳,表现为高密度.青藏高原内部中地壳强低密度层主要分布在高原边界.其成因是印度板块俯冲和周围坚硬块体阻挡作用导致在高原边界形成的高应变积累闭锁区,为壳内低密度软弱物质的形成提供了条件.（2）80 km深度,青藏高原上地幔顶部显示出低密度的特征.高原内部东、中、西密度特征差异明显,低密度以95°E为中心线呈东西对称分布.以班公—怒江缝合带为中心,在拉萨块体和羌塘块体内从北向南出现了"低-高-低"的密度分布起伏特征.该特征与GRACE得到的莫霍面起伏特征一致,结合大地构造结果,这种起伏特征验证了印度、羌塘块体从南北两侧分别向喜马拉雅、拉萨地块挤入的双向俯冲模式.（3）四川盆地和鄂尔多斯盆地内,地壳高密度异常较地震波速异常明显偏低,表明古老的四川盆地和鄂尔多斯盆地比想象中更冷、更坚硬.塔里木盆地和柴达木盆地内壳、幔高密度的结构特征,对应地幔物质上涌.     

Non-Singular Expressions for the Gravity Gradients in the Local North-Oriented and Orbital Reference Frames

The conventional expansions of the gravity gradients in the local north-oriented reference frame have a complicated form, depending on the first- and second-order derivatives of the associated Legendre functions of the colatitude and containing factors which tend to infinity when approaching the poles. In the present paper, the general term of each of these series is transformed to a product of a geopotential coefficient and a sum of several adjacent Legendre functions of the colatitude multiplied by a function of the longitude. These transformations are performed on the basis of relations between the Legendre functions and their derivatives published by Ilk (1983). The second-order geopotential derivatives corresponding to the local orbital reference frame are presented as linear functions of the north-oriented gravity gradients. The new expansions for the latter are substituted into these functions. As a result, the orbital derivatives are also presented as series depending on the geopotential coefficients multiplied by sums of the Legendre functions whose coefficients depend on the longitude and the satellite track azimuth at an observation point. The derived expansions of the observables can be applied for constructing a geopotential model from the GOCE mission data by the time-wise and space-wise approaches. The numerical experiments demonstrate the correctness of the analytical formulas.An erratum to this article can be found at     

重力卫星和测高卫星五年来的进展

对近年升空的重力卫星CHAMP和GRACE及将于近期升空的GOCE卫星,以及测高卫星ENVISAT,Jason1,GFO,ICESat和CryoSat的技术特点进行了介绍和评估。     

Numerical study on the mixed model in the GOCE polar gap problem

Gravity gradients acquired by the Gravity field and steady-state Ocean Circulation Explorer(GOCE) do not cover the entire earth because of its sun-synchronous orbit leaving data gaps with a radius of about 6.5° in the polar regions.Previous studies showed that the loss of data in the polar regions deteriorates the accuracy of the low order(or near zonal) coefficients of the earth gravity model,which is the so-called polar gap problem in geodesy.In order to find a stable solution for the earth gravity model from the GOCE gravity gradients,three models,i.e.the Gauss-Markov model,light constraint model and the mixed model,are compared and evaluated numerically with the gravity gradient simulated with the EGM2008.The comparison shows that the Best Linear Uniformly Unbiased Estimation(BLUUE) estimator of the mixed model can solve the polar gap problem as effectively as the light constraint model;furthermore,the mixed model is more rigorous in dealing with the supplementary information and leads to a better accuracy in determining the global geoid.     

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