GOCE gravitational gradiometry

GOCE is the first gravitational gradiometry satellite mission. Gravitational gradiometry is the measurement of the second derivatives of the gravitational potential. The nine derivatives form a 3 × 3 matrix, which in geodesy is referred to as Marussi tensor. From the basic properties of the gravitational field, it follows that the matrix is symmetric and trace free. The latter property corresponds to Laplace equation, which gives the theoretical foundation of its representation in terms of spherical harmonic or Fourier series. At the same time, it provides the most powerful quality check of the actual measured gradients. GOCE gradiometry is based on the principle of differential accelerometry. As the satellite carries out a rotational motion in space, the accelerometer differences contain angular effects that must be removed. The GOCE gradiometer provides the components V _xx , V _yy , V _zz and V _xz with high precision, while the components V _xy and V _yz are of low precision, all expressed in the gradiometer reference frame. The best performance is achieved inside the measurement band from 5 × 10^–3 to 0.1 Hz. At lower frequencies, the noise increases with 1/f and is superimposed by cyclic distortions, which are modulated from the orbit and attitude motion into the gradient measurements. Global maps with the individual components show typical patterns related to topographic and tectonic features. The maps are separated into those for ascending and those for descending tracks as the components are expressed in the instrument frame. All results are derived from the measurements of the period from November to December 2009. While the components V _xx and V _yy reach a noise level of about \({10\;\rm{\frac{mE}{\sqrt{Hz}}}}\), that of V _zz and V _xz is about \({20\; \rm{\frac{mE}{\sqrt{Hz}}}}\). The cause of the latter’s higher noise is not yet understood. This is also the reason why the deviation from the Laplace condition is at the \({20 \;\rm{\frac{mE}{\sqrt{Hz}}}}\) level instead of the originally planned \({11\;\rm{\frac{mE}{\sqrt{Hz}}}}\). Each additional measurement cycle will improve the accuracy and to a smaller extent also the resolution of the spherical harmonic coefficients derived from the measured gradients.     

Improving GOCE cross-track gravity gradients

The GOCE gravity gradiometer measured highly accurate gravity gradients along the orbit during GOCE’s mission lifetime from March 17, 2009, to November 11, 2013. These measurements contain unique information on the gravity field at a spatial resolution of 80 km half wavelength, which is not provided to the same accuracy level by any other satellite mission now and in the foreseeable future. Unfortunately, the gravity gradient in cross-track direction is heavily perturbed in the regions around the geomagnetic poles. We show in this paper that the perturbing effect can be modeled accurately as a quadratic function of the non-gravitational acceleration of the satellite in cross-track direction. Most importantly, we can remove the perturbation from the cross-track gravity gradient to a great extent, which significantly improves the accuracy of the latter and offers opportunities for better scientific exploitation of the GOCE gravity gradient data set.     

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 orbit predictions for SLR tracking

After a descent phase of about half a year, the Gravity field and steady-state Ocean Circulation Explorer (GOCE) reached the final orbital altitude of the first measurement and operational phase (MOP-1) in September 2009. Due to this very low orbital altitude and the inactive drag compensation during descent, the generation of reliable predictions of the GOCE trajectory turned out to be a major challenge even for short prediction intervals. As predictions of good quality are a prerequisite for frequent ranging from the tracking network of the International Laser Ranging Service (ILRS), Satellite Laser Ranging (SLR) data of GOCE was very sparse at mission start and made it difficult to independently calibrate and optimize the orbit determination based on data of the Global Positioning System (GPS). In addition to the GOCE orbit predictions provided by the European Space Agency (ESA), the Astronomical Institute of the University of Bern (AIUB) started providing predictions on July 22, 2009, as part of the Level 1b to Level 2 data processing performed at AIUB. The predictions based on the 12-h ultra-rapid products of the International GNSS Service (IGS) were originally intended to primarily serve the daylight passes in the early evening hours over Europe. The corresponding along-track prediction errors were often kept below 50 m during the descent phase and allowed for the first successful SLR tracking of GOCE over Europe on July 29, 2009, by the Zimmerwald observatory. Additional predictions based on the IGS 18-h ultra-rapid products are provided by AIUB since September 20, 2009, to further optimize the GOCE SLR tracking. In this article, the development of the GOCE prediction service at AIUB is presented, and the quality of the orbit predictions is assessed for periods with and without active drag compensation. The prediction quality is discussed as a function of the prediction interval, the quality of the input products for the GPS satellite orbits and clocks, and the availability of the GOCE GPS data. From the methodological point of view, different approaches for the treatment of the non-gravitational accelerations acting on the GOCE satellite are discussed and their impact on the prediction quality is assessed, in particular during the descent phase. Eventually, an outlook is given on the significance of GOCE SLR tracking to identify systematic errors in the GPS-based orbit determination, e.g., cross-track errors induced by mismodeled GOCE GPS phase center variations (PCVs).     

GOCE引力梯度内部校准方法

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

GOCE: precise orbit determination for the entire mission

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) was the first Earth explorer core mission of the European Space Agency. It was launched on March 17, 2009 into a Sun-synchronous dusk-dawn orbit and re-entered into the Earth’s atmosphere on November 11, 2013. The satellite altitude was between 255 and 225 km for the measurement phases. The European GOCE Gravity consortium is responsible for the Level 1b to Level 2 data processing in the frame of the GOCE High-level processing facility (HPF). The Precise Science Orbit (PSO) is one Level 2 product, which was produced under the responsibility of the Astronomical Institute of the University of Bern within the HPF. This PSO product has been continuously delivered during the entire mission. Regular checks guaranteed a high consistency and quality of the orbits. A correlation between solar activity, GPS data availability and quality of the orbits was found. The accuracy of the kinematic orbit primarily suffers from this. Improvements in modeling the range corrections at the retro-reflector array for the SLR measurements were made and implemented in the independent SLR validation for the GOCE PSO products. The satellite laser ranging (SLR) validation finally states an orbit accuracy of 2.42 cm for the kinematic and 1.84 cm for the reduced-dynamic orbits over the entire mission. The common-mode accelerations from the GOCE gradiometer were not used for the official PSO product, but in addition to the operational HPF work a study was performed to investigate to which extent common-mode accelerations improve the reduced-dynamic orbit determination results. The accelerometer data may be used to derive realistic constraints for the empirical accelerations estimated for the reduced-dynamic orbit determination, which already improves the orbit quality. On top of that the accelerometer data may further improve the orbit quality if realistic constraints and state-of-the-art background models such as gravity field and ocean tide models are used for the reduced-dynamic orbit determination.     

Preprocessing of gravity gradients at the GOCE high-level processing facility

One of the products derived from the gravity field and steady-state ocean circulation explorer (GOCE) observations are the gravity gradients. These gravity gradients are provided in the gradiometer reference frame (GRF) and are calibrated in-flight using satellite shaking and star sensor data. To use these gravity gradients for application in Earth scienes and gravity field analysis, additional preprocessing needs to be done, including corrections for temporal gravity field signals to isolate the static gravity field part, screening for outliers, calibration by comparison with existing external gravity field information and error assessment. The temporal gravity gradient corrections consist of tidal and nontidal corrections. These are all generally below the gravity gradient error level, which is predicted to show a 1/f behaviour for low frequencies. In the outlier detection, the 1/f error is compensated for by subtracting a local median from the data, while the data error is assessed using the median absolute deviation. The local median acts as a high-pass filter and it is robust as is the median absolute deviation. Three different methods have been implemented for the calibration of the gravity gradients. All three methods use a high-pass filter to compensate for the 1/f gravity gradient error. The baseline method uses state-of-the-art global gravity field models and the most accurate results are obtained if star sensor misalignments are estimated along with the calibration parameters. A second calibration method uses GOCE GPS data to estimate a low-degree gravity field model as well as gravity gradient scale factors. Both methods allow to estimate gravity gradient scale factors down to the 10⁻³ level. The third calibration method uses high accurate terrestrial gravity data in selected regions to validate the gravity gradient scale factors, focussing on the measurement band. Gravity gradient scale factors may be estimated down to the 10⁻² level with this method.     

Spherical integral formulas for upward/downward continuation of gravitational gradients onto gravitational gradients

New integral formulas for upward/downward continuation of gravitational gradients onto gravitational gradients are derived in this article. They provide more options for continuation of gravitational gradient combinations and extend available mathematical apparatus formulated for this purpose up to now. The starting point represents the analytical solution of the spherical gradiometric boundary value problem in the spatial domain. Applying corresponding differential operators on the analytical solution of the spherical gradiometric boundary value problem, a total of 18 integral formulas are provided. Spatial and spectral forms of isotropic kernels are given and their behaviour for parameters of a GOCE-like satellite is investigated. Correctness of the new integral formulas and the isotropic kernels is tested in a closed-loop simulation. The derived integral formulas and the isotropic kernels form a theoretical basis for validation purposes and geophysical applications of satellite gradiometric data as provided currently by the GOCE mission. They also extend the well-known Meissl scheme.     

GPS-based precise orbit determination of the very low Earth-orbiting gravity mission GOCE

 A prerequisite for the success of future gravity missions like the European Gravity field and steady-state Ocean Circulation Explorer (GOCE) is a precise orbit determination (POD). A detailed simulation study has been carried out to assess the achievable orbit accuracy based on satellite-to-satellite tracking (SST) by the US global positioning system (GPS) and in conjunction the implications for gravity field determination. An orbit accuracy at the few centimeter level seems possible, sufficient to support the GOCE gravity mission and in particular its gravity gradiometer. Received: 21 January 2000 / Accepted: 4 July 2000     

Exploring gravity field determination from orbit perturbations of the European Gravity Mission GOCE

 A comparison was made between two methods for gravity field recovery from orbit perturbations that can be derived from global positioning system satellite-to-satellite tracking observations of the future European gravity field mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer). The first method is based on the analytical linear orbit perturbation theory that leads under certain conditions to a block-diagonal normal matrix for the gravity unknowns, significantly reducing the required computation time. The second method makes use of numerical integration to derive the observation equations, leading to a full set of normal equations requiring powerful computer facilities. Simulations were carried out for gravity field recovery experiments up to spherical harmonic degree and order 80 from 10 days of observation. It was found that the first method leads to large approximation errors as soon as the maximum degree surpasses the first resonance orders and great care has to be taken with modeling resonance orbit perturbations, thereby loosing the block-diagonal structure. The second method proved to be successful, provided a proper division of the data period into orbital arcs that are not too long. Received: 28 April 2000 / Accepted: 6 November 2000     

Validation of GOCE gravity field models by means of orbit residuals and geoid comparisons

Three GOCE-based gravity field solutions have been computed by ESA’s high-level processing facility and were released to the user community. All models are accompanied by variance-covariance information resulting either from the least squares procedure or a Monte-Carlo approach. In order to obtain independent external quality parameters and to assess the current performance of these models, a set of independent tests based on satellite orbit determination and geoid comparisons is applied. Both test methods can be regarded as complementary because they either investigate the performance in the long wavelength spectral domain (orbit determination) or in the spatial domain (geoid comparisons). The test procedure was applied to the three GOCE gravity field solutions and to a number of selected pre-launch models for comparison. Orbit determination results suggest, that a pure GOCE gravity field model does not outperform the multi-year GRACE gravity field solutions. This was expected as GOCE is designed to improve the determination of the medium to high frequencies of the Earth gravity field (in the range of degree and order 50 to 200). Nevertheless, in case of an optimal combination of GOCE and GRACE data, orbit determination results should not deteriorate. So this validation procedure can also be used for testing the optimality of the approach adopted for producing combined GOCE and GRACE models. Results from geoid comparisons indicate that with the 2 months of GOCE data a significant improvement in the determination of the spherical harmonic spectrum of the global gravity field between degree 50 and 200 can be reached. Even though the ultimate mission goal has not yet been reached, especially due to the limited time span of used GOCE data (only 2 months), it was found that existing satellite-only gravity field models, which are based on 7 years of GRACE data, can already be enhanced in terms of spatial resolution. It is expected that with the accumulation of more GOCE data the gravity field model resolution and quality can be further enhanced, and the GOCE mission goal of 1–2 cm geoid accuracy with 100 km spatial resolution can be achieved.     

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.     

Efficient propagation of error covariance matrices of gravitational models: application to GRACE and GOCE

We have applied efficient methods for computing variances and covariances of functions of a global gravity field model expanded in spherical harmonics, using the full variance–covariance matrix of the coefficients. Examples are given with recent models derived from GRACE (up to degree and order 150), and with simulated GOCE derived solutions (up to degree and order 200).     

The ITG-Goce02 gravity field model from GOCE orbit and gradiometer data based on the short arc approach

In this contribution, we describe the global GOCE-only gravity field model ITG-Goce02 derived from 7.5 months of gradiometer and orbit data. This model represents an alternative to the official ESA products as it is computed completely independently, using a different processing strategy and a separate software package. Our model is derived using the short arc approach, which allows a very effective decorrelation of the highly correlated GOCE gradiometer and orbit data noise by introducing a full empirical covariance matrix for each arc, and gives the possibility to downweight ‘bad’ arcs. For the processing of the orbit data we rely on the integral equation approach instead of the energy integral method, which has been applied in several other GOCE models. An evaluation against high-resolution global gravity field models shows very similar differences of our model compared to the official GOCE results published by ESA (release 2), especially to the model derived by the time-wise approach. This conclusion is confirmed by comparison of the GOCE models to GPS/levelling and altimetry data.     

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.     

GPS-derived orbits for the GOCE satellite

The first ESA (European Space Agency) Earth explorer core mission GOCE (Gravity field and steady-state Ocean Circulation Explorer) was launched on 17 March 2009 into a sun-synchronous dusk–dawn orbit with an exceptionally low initial altitude of about 280 km. The onboard 12-channel dual-frequency GPS (Global Positioning System) receiver delivers 1 Hz data, which provides the basis for precise orbit determination (POD) for such a very low orbiting satellite. As part of the European GOCE Gravity Consortium the Astronomical Institute of the University of Bern and the Department of Earth Observation and Space Systems are responsible for the orbit determination of the GOCE satellite within the GOCE High-level Processing Facility. Both quick-look (rapid) and very precise orbit solutions are produced with typical latencies of 1 day and 2 weeks, respectively. This article summarizes the special characteristics of the GOCE GPS data, presents POD results for about 2 months of data, and shows that both latency and accuracy requirements are met. Satellite Laser Ranging validation shows that an accuracy of 4 and 7 cm is achieved for the reduced-dynamic and kinematic Rapid Science Orbit solutions, respectively. The validation of the reduced-dynamic and kinematic Precise Science Orbit solutions is at a level of about 2 cm.     

GOCE重力卫星在军事上的应用前景分析

对GOCE重力卫星进行了详细的介绍;分析了星上载荷的配置和功能.结合GOCE任务的科学目标,详细分析了其在军事上的应用前景,并列出了军事上重要的应用方向.希望广大测量界人士对GoCE重力卫星应用领域进行更多的利用.     

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

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