利用GOCE卫星轨道数据恢复地球重力场模型方法的分析

欧空局早期公布的时域法和空域法解算的GOCE模型均采用能量守恒法处理轨道数据, 但恢复的长波重力场信号精度较低, 而且GOCE卫星在两极存在数据空白, 利用其观测数据恢复重力场模型是一个不适定问题, 导致解算的模型带谐项精度较低, 需进行正则化处理。本文分析了基于轨道数据恢复重力场模型的方法用于处理GOCE数据的精度, 对最优正则化方法和参数的选择进行研究。利用GOCE卫星2009-11-01—2010-01-31共92 d的精密轨道数据, 采用不依赖先验信息的能量守恒法、短弧积分法和平均加速度法恢复GOCE重力场模型, 利用Tikhonov正则化技术处理病态问题。结果表明, 平均加速度法恢复模型的精度最高, 能量守恒法的精度最低, 短弧积分法的精度稍差于平均加速度法。未来联合处理轨道和梯度数据时, 建议采用平均加速度法或短弧积分法处理轨道数据, 并且轨道数据可有效恢复120阶次左右的模型。Kaula正则化和SOT处理GOCE病态问题的效果最好, 并且两者对应的最优正则化参数基本一致, 但利用正则化技术不能完全抑制极空白问题的影响, 需要联合GRACE等其他数据才能获得理想的结果。     

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

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

利用最小二乘直接法反演卫星重力场模型的MPI并行算法

针对海量卫星重力数据反演高阶次地球重力场模型的密集型计算任务与高内存耗用问题,基于MPI实现了最小二乘直接法恢复高阶次位系数的并行算法。引入并行读写、分块存储与分块计算等方式完成了设计矩阵的构建、法方程的形成与求解等密集型计算任务的并行算法,数值计算结果表明三者的并行相对效率峰值可分别达到95%、68%、63%。利用GOCE轨道跟踪和径向扰动重力梯度数据(共518 400个历元)分别反演了120、240阶次地球重力场模型,计算时间仅为40 min、7 h,内存耗用峰值仅为290 MB、1.57 GB;采用与GOCE同等噪声水平的观测数据恢复的重力场模型精度与GOCE已发布模型的解算精度相一致,联合GRACE和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.     

引力场梯度张量的非奇异公式推导

传统的引力场梯度计算公式在两极附近存在奇异性,需要换用其他的非奇异计算公式。从奇异性产生的原因入手,并结合勒让德函数的有关性质,推导了一组新的计算公式。实际计算验证了该公式的正确性和有效性。     

First GOCE gravity field models derived by three different approaches

Three gravity field models, parameterized in terms of spherical harmonic coefficients, have been computed from 71 days of GOCE (Gravity field and steady-state Ocean Circulation Explorer) orbit and gradiometer data by applying independent gravity field processing methods. These gravity models are one major output of the European Space Agency (ESA) project GOCE High-level Processing Facility (HPF). The processing philosophies and architectures of these three complementary methods are presented and discussed, emphasizing the specific features of the three approaches. The resulting GOCE gravity field models, representing the first models containing the novel measurement type of gravity gradiometry ever computed, are analysed and assessed in detail. Together with the coefficient estimates, full variance-covariance matrices provide error information about the coefficient solutions. A comparison with state-of-the-art GRACE and combined gravity field models reveals the additional contribution of GOCE based on only 71 days of data. Compared with combined gravity field models, large deviations appear in regions where the terrestrial gravity data are known to be of low accuracy. The GOCE performance, assessed against the GRACE-only model ITG-Grace2010s, becomes superior at degree 150, and beyond. GOCE provides significant additional information of the global Earth gravity field, with an accuracy of the 2-month GOCE gravity field models of 10?cm in terms of geoid heights, and 3?mGal in terms of gravity anomalies, globally at a resolution of 100?km (degree/order 200).     

联合高低卫-卫跟踪和卫星重力梯度数据恢复地球重力场的谱组合法

GOCE采用的高低卫-卫跟踪和卫星重力梯度测量技术在恢复重力场方面各有所长并互为补充,如何有效利用这两类观测数据最优确定地球重力场是GOCE重力场反演的关键问题。本文研究了联合高低卫-卫跟踪和卫星重力梯度数据恢复地球重力场的最小二乘谱组合法,基于球谐分析方法推导并建立了卫星轨道面扰动位T和径向重力梯度Tzz、以及扰动位T和重力梯度分量组合{Tzz-Txx-Tyy}的谱组合计算模型与误差估计公式。数值模拟结果表明,谱组合计算模型可以有效顾及各类数据的精度和频谱特性进行最优联合求解。采用61天GOCE实测数据反演的两个180阶次地球重力场模型WHU_GOCE_SC01S（扰动位和径向重力梯度数据求解）和WHU_GOCE_SC02S（扰动位和重力梯度分量组合数据求解）,结果显示后者精度优于前者,并且它们的整体精度优于GOCE时域解,而与GOCE空域解的精度接近,验证了谱组合法的可行性与有效性。     

Mission design,operation and exploitation of the gravity field and steady-state ocean circulation explorer mission

The European Space Agency’s Gravity field and steady-state ocean circulation explorer mission (GOCE) was launched on 17 March 2009. As the first of the Earth Explorer family of satellites within the Agency’s Living Planet Programme, it is aiming at a better understanding of the Earth system. The mission objective of GOCE is the determination of the Earth’s gravity field and geoid with high accuracy and maximum spatial resolution. The geoid, combined with the de facto mean ocean surface derived from twenty-odd years of satellite radar altimetry, yields the global dynamic ocean topography. It serves ocean circulation and ocean transport studies and sea level research. GOCE geoid heights allow the conversion of global positioning system (GPS) heights to high precision heights above sea level. Gravity anomalies and also gravity gradients from GOCE are used for gravity-to-density inversion and in particular for studies of the Earth’s lithosphere and upper mantle. GOCE is the first-ever satellite to carry a gravitational gradiometer, and in order to achieve its challenging mission objectives the satellite embarks a number of world-first technologies. In essence the spacecraft together with its sensors can be regarded as a spaceborne gravimeter. In this work, we describe the mission and the way it is operated and exploited in order to make available the best-possible measurements of the Earth gravity field. The main lessons learned from the first 19 months in orbit are also provided, in as far as they affect the quality of the science data products and therefore are of specific interest for GOCE data users.     

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.     

Monitoring GOCE gradiometer calibration parameters using accelerometer and star sensor data: methodology and first results

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite, launched on 17 March 2009, is designed to measure the Earth’s mean gravity field with unprecedented accuracy at spatial resolutions down to 100?km. The accurate calibration of the gravity gradiometer on-board GOCE is of utmost importance for achieving the mission goals. ESA’s baseline method for the calibration uses star sensor and accelerometer data of a dedicated calibration procedure, which is executed every 2?months. In this paper, we describe a method for monitoring the evolution of calibration parameter during that time. The method works with star sensor and accelerometer data and does not require gravity field models, which distinguishes it from other existing methods. We present time series of calibration parameters estimated from GOCE data from 1 November 2009 to 17 May 2010. The time series confirm drifts in the calibration parameters that are present in the results of other methods, including ESA’s baseline method. Although these drifts are very small, they degrade the gravity gradients, leading to the conclusion that the calibration parameters of the ESA’s baseline method need to be linearly interpolated. Further, we find a correction of ?36 × 10^?6 for one calibration parameter (in-line differential scale factor of the cross-track gradiometer arm), which improves the gravity gradient performance. The results are validated by investigating the trace of the calibrated gravity gradients and comparing calibrated gravity gradients with reference gradients computed along the GOCE orbit using the ITG-Grace-2010s gravity field model.     

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.     

Improvements in height datum transfer expected from the GOCE mission

 One of the aims of the Earth Explorer Gravity Field and Steady-State Ocean Circulation (GOCE) mission is to provide global and regional models of the Earth's gravity field and of the geoid with high spatial resolution and accuracy. Using the GOCE error model, simulation studies were performed in order to estimate the accuracy of datum transfer in different areas of the Earth. The results showed that with the GOCE error model, the standard deviation of the height anomaly differences is about one order of magnitude better than the corresponding value with the EGM96 error model. As an example, the accuracy of the vertical datum transfer from the tide gauge of Amsterdam to New York was estimated equal to 57 cm when the EGM96 error model was used, while in the case of GOCE error model this accuracy was increased to 6 cm. The geoid undulation difference between the two places is about 76.5 m. Scaling the GOCE errors to the local gravity variance, the estimated accuracy varied between 3 and 7 cm, depending on the scaling model. Received: 1 March 2000 / Accepted: 21 February 2001     

Filter design for GOCE gravity gradients

The GOCE satellite observes gravity gradients with unprecedented accuracy and resolution. The GOCE observations are reliable within a well-defined measurement bandwidth. In this study, different finite and infinite impulse response filters have been designed to obtain the demanded pass. Exhaustive time and frequency domain investigations prove that the proposed infinite impulse response filter can be a real competitor of the existing solution of the filtering problem.     

GOCE引力梯度内部校准方法

GOCE卫星引力梯度仪的精确校准是反演高精度重力场的前提之一,本文利用GOCE卫星L1b数据中的引力梯度仪及恒星敏感器数据实现了卫星引力梯度的内部校准。以最小二乘联合多个恒星敏感器观测数据确定内部校准使用的角速度,有效避免了单个恒星敏感器低精度角速度分量对坐标转换过程的影响。考虑到恒星敏感器坐标系与梯度仪坐标系间旋转矩阵随时间的变化,本文在ESA官方内部校准方法的基础上,提出了顾及旋转矩阵校准参数的内部校准模型,并利用2009年11月的GOCE实测数据验证了该方法的效果。结果表明,该旋转矩阵校准参数数值约100″,且在该月存在3″~30″的漂移;与GOCE官方内部校准方法对比,从卫星引力梯度精度结果来看,在低于0.005 Hz频段内,同时解算旋转矩阵的校准参数与梯度仪内3个加速度计对的校准参数的内部校准模型优于仅考虑加速度计对校准参数的模型;除此之外,本文讨论了以该模型为基础的GOCE梯度仪数据校准的可能方法,为GOCE及后续重力卫星的数据处理工作提供参考。     

The static gravity field model DGM-1S from GRACE and GOCE data: computation, validation and an analysis of GOCE mission’s added value

We present a global static model of the Earth’s gravity field entitled DGM-1S based on GRACE and GOCE data. The collection of used data sets includes nearly 7 years of GRACE KBR data and 10 months of GOCE gravity gradient data. The KBR data are transformed with a 3-point differentiation into quantities that are approximately inter-satellite accelerations. Gravity gradients are processed in the instrumental frame. Noise is handled with a frequency-dependent data weighting. DGM-1S is complete to spherical harmonic degree 250 with a Kaula regularization being applied above degree 179. Its performance is compared with a number of other satellite-only GRACE/GOCE models by confronting them with (i) an independent model of the oceanic mean dynamic topography, and (ii) independent KBR and gravity gradient data. The tests reveal a competitive quality for DGM-1S. Importantly, we study added value of GOCE data by comparing the performance of satellite-only GRACE/GOCE models with models produced without GOCE data: either ITG-Grace2010s or EGM2008 depending on which of the two performs better in a given region. The test executed based on independent gravity gradients quantifies this added value as 25–38 % in the continental areas poorly covered with terrestrial gravimetry data (Equatorial Africa, Himalayas, and South America), 7–17 % in those with a good coverage with these data (Australia, North America, and North Eurasia), and 14 % in the oceans. This added value is shown to be almost entirely related to coefficients below degree 200. It is shown that this gain must be entirely attributed to gravity gradients acquired by the mission. The test executed based on an independent model of the mean dynamic topography suggests that problems still seem to exist in satellite-only GRACE/GOCE models over the Pacific ocean, where noticeable deviations between these models and EGM2008 are detected, too.     

Evaluation of the first GOCE static gravity field models using terrestrial gravity,vertical deflections and EGM2008 quasigeoid heights

Recently, four global geopotential models (GGMs) were computed and released based on the first 2 months of data collected by the Gravity field and steady-state Ocean Circulation Explorer (GOCE) dedicated satellite gravity field mission. Given that GOCE is a technologically complex mission and different processing strategies were applied to real space-collected GOCE data for the first time, evaluation of the new models is an important aspect. As a first assessment strategy, we use terrestrial gravity data over Switzerland and Australia and astrogeodetic vertical deflections over Europe and Australia as ground-truth data sets for GOCE model evaluation. We apply a spectral enhancement method (SEM) to the truncated GOCE GGMs to make their spectral content more comparable with the terrestrial data. The SEM utilises the high-degree bands of EGM2008 and residual terrain model data as a data source to widely bridge the spectral gap between the satellite and terrestrial data. Analysis of root mean square (RMS) errors is carried out as a function of (i) the GOCE GGM expansion degree and (ii) the four different GOCE GGMs. The RMS curves are also compared against those from EGM2008 and GRACE-based GGMs. As a second assessment strategy, we compare global grids of GOCE GGM and EGM2008 quasigeoid heights. In connection with EGM2008 error estimates, this allows location of regions where GOCE is likely to deliver improved knowledge on the Earth’s gravity field. Our ground truth data sets, together with the EGM2008 quasigeoid comparisons, signal clear improvements in the spectral band ~160–165 to ~180–185 in terms of spherical harmonic degrees for the GOCE-based GGMs, fairly independently of the individual GOCE model used. The results from both assessments together provide strong evidence that the first 2 months of GOCE observations improve the knowledge of the Earth’s static gravity field at spatial scales between ~125 and ~110 km, particularly over parts of Asia, Africa, South America and Antarctica, in comparison with the pre-GOCE-era.     

现代低轨卫星对地球重力场探测的实践和进展

综述了现代低轨卫星对地球重力场测量的特点和近况，介绍了已经和即将发射的重力卫星CHAMP、GRACE、GOCE和新型测高卫星，讨论了作为现代重力卫星首次实践－－CHAMP卫星的进展和目前尚待解决的问题。     

Geoid and high resolution sea surface topography modelling in the mediterranean from gravimetry,altimetry and GOCE data: evaluation by simulation

The determination of local geoid models has traditionally been carried out on land and at sea using gravity anomaly and satellite altimetry data, while it will be aided by the data expected from satellite missions such as those from the Gravity field and steady-state ocean circulation explorer (GOCE). To assess the performance of heterogeneous data combination to local geoid determination, simulated data for the central Mediterranean Sea are analyzed. These data include marine and land gravity anomalies, altimetric sea surface heights, and GOCE observations processed with the space-wise approach. A spectral analysis of the aforementioned data shows their complementary character. GOCE data cover long wavelengths and account for the lack of such information from gravity anomalies. This is exploited for the estimation of local covariance function models, where it is seen that models computed with GOCE data and gravity anomaly empirical covariance functions perform better than models computed without GOCE data. The geoid is estimated by different data combinations and the results show that GOCE data improve the solutions for areas covered poorly with other data types, while also accounting for any long wavelength errors of the adopted reference model that exist even when the ground gravity data are dense. At sea, the altimetric data provide the dominant geoid information. However, the geoid accuracy is sensitive to orbit calibration errors and unmodeled sea surface topography (SST) effects. If such effects are present, the combination of GOCE and gravity anomaly data can improve the geoid accuracy. The present work also presents results from simulations for the recovery of the stationary SST, which show that the combination of geoid heights obtained from a spherical harmonic geopotential model derived from GOCE with satellite altimetry data can provide SST models with some centimeters of error. However, combining data from GOCE with gravity anomalies in a collocation approach can result in the estimation of a higher resolution geoid, more suitable for high resolution mean dynamic SST modeling. Such simulations can be performed toward the development and evaluation of SST recovery methods.     

GOCE卫星重力梯度校准与无阻尼控制效果分析

GOCE卫星是首颗搭载高精度梯度仪,通过加速度计差分测量确定地球重力场的现代重力卫星。该卫星设计为无阻尼飞行状态（沿轨方向）,加速度计并未安置在卫星质心,这些特点使得GOCE与标准的卫星跟踪卫星重力测量模式有着显著的区别。本文首先指出GOCE任务中普通模式加速度校准存在不严密性问题,并提出了分别校准6个加速度计,分离偏差参数的方案。利用GOCE任务期内的几何法精密轨道,采用动力法完成校准,并分析了无阻尼控制的效果,发现：①虽然GOCE所在轨道高度的中性大气密度较GRACE高两到三个量级,但GOCE卫星在沿轨方向的残余非保守力比GRACE卫星的对应分量小一个量级,充分显示了无阻尼控制系统的补偿效果;②通过精密轨道内插的轨道速度与动力法轨道速度的比较可以得出,卫星无阻尼控制系统对GOCE卫星速度的显著影响;③计算了GOCE卫星所受的非保守力。获得了GOCE任务期间的加速度计校准参数,并讨论了利用其辅助重力梯度仪数据预处理的可能方法。