Error analyses of CHAMP data for recovery of the Earth’s gravity field

A preliminary commission error analysis whereby orbit perturbation theory and other techniques are used to assess and predict the recovery of the Earths gravity field from the challenging microsatellite payload (CHAMP) mission is developed and implemented. With CHAMP launched in July 2000, accumulated evidence is now available to quantify the errors in the recovery procedure including the orbital precision from GPS, attitude errors, accelerometer noise and thruster mismatch/misalignment. For the latter, numerical integrations using a variable length single-step Runge–Kutta integrator and a fixed length multi-step method are compared to assess the error associated with assuming that the thruster misalignment can be spread uniformly across a step interval. Error degree variances from simulated studies are compared to results from a recently released CHAMP-based gravity field, EIGEN-1S. It is seen that the orbital positioning, as derived from the onboard GPS receiver, is critical, with accelerometer noise contributing at a lower level. Attitude error, at currently quoted accuracy, is not significant as an error source. AcknowledgementsThe authors would like to thank the UK Natural Environment Research Council (Grant No. NER/A/0000/00612) for financing this study and GFZ for supplying the data and technical support.     

A technique for modeling the Earth’s gravity field on the basis of satellite accelerations

A technique is proposed for Earths gravity field modeling on the basis of satellite accelerations that are derived from precise orbit data. The functional model rests on Newtons second law. The computational procedure is based on the pre-conditioned conjugate-gradient (PCCG) method. The data are treated as weighted average accelerations rather than as point-wise ones. As a result, a simple three-point numerical differentiation scheme can be used to derive them. Noise in the orbit-derived accelerations is strongly dependent on frequency. Therefore, the key element of the proposed technique is frequency-dependent data weighting. Fast convergence of the PCCG procedure is ensured by a block-diagonal pre-conditioner (approximation of the normal matrix), which is derived under the so-called Colombo assumptions. Both uninterrupted data sets and data with gaps can be handled. The developed technique is compared with other approaches: (1) the energy balance approach (based on the energy conservation law) and (2) the traditional approach (based on the integration of variational equations). Theoretical considerations, supported by a numerical study, show that the proposed technique is more accurate than the energy balance approach and leads to approximately the same results as the traditional one. The former finding is explained by the fact that the energy balance approach is only sensitive to the along-track force component. Information about the cross-track and the radial component of the gravitational potential gradient is lost because the corresponding force components do no work and do not contribute to the energy balance. Furthermore, it is shown that the proposed technique is much (possibly, orders of magnitude) faster than the traditional one because it does not require the computation of the normal matrix. Hints are given on how the proposed technique can be adapted to the explicit assembling of the normal matrix if the latter is needed for the computation of the model covariance matrix.Acknowledgments. Professor R. Klees is thanked for support of the project and for numerous fruitful discussions. The authors are also thankful to Dr. J. Kusche for useful remarks and to Dr. E. Schrama, his solid background in satellite geodesy proved to be very helpful. A large number of valuable comments were made by Dr. S.-C. Han, Dr. P. Schwintzer, and an anonymous reviewer; their contribution is greatly acknowledged. The satellite orbits used in the numerical study were kindly provided by Dr. P. Visser (Aerospace Department, Delft University of Technology). Access to the SGI Origin 3800 computer was provided by Stichting Nationale Computerfaciliteiten (NCF), grant SG-027.     

Efficient and accurate high-degree spherical harmonic synthesis of gravity field functionals at the Earth’s surface using the gradient approach

Spherical harmonic synthesis (SHS) of gravity field functionals at the Earth’s surface requires the use of heights. The present study investigates the gradient approach as an efficient yet accurate strategy to incorporate height information in SHS at densely spaced multiple points. Taylor series expansions of commonly used functionals quasigeoid heights, gravity disturbances and vertical deflections are formulated, and expressions of their radial derivatives are presented to arbitrary order. Numerical tests show that first-order gradients, as introduced by Rapp (J Geod 71(5):282–289, 1997) for degree 360 models, produce cm- to dm-level RMS approximation errors over rugged terrain when applied with EGM2008 to degree 2190. Instead, higher-order Taylor expansions are recommended that are capable of reducing approximation errors to insignificance for practical applications. Because the height information is separated from the actual synthesis, the gradient approach can be applied along with existing highly efficient SHS routines to compute surface functionals at arbitrarily dense grid points. This confers considerable computational savings (above or well above one order of magnitude) over conventional point-by-point SHS. As an application example, an ultra-high resolution model of surface gravity functionals (EurAlpGM2011) is constructed over the entire European Alps that incorporates height information in the SHS at 12,000,000 surface points. Based on EGM2008 and residual topography data, quasigeoid heights, gravity disturbances and vertical deflections are estimated at ~200m resolution. As a conclusion, the gradient approach is efficient and accurate for high-degree SHS at multiple points at the Earth’s surface.     

Earth’s gravity field modelling based on satellite accelerations derived from onboard GPS phase measurements

GPS data collected by satellite gravity missions can be used for extracting the long-wavelength part of the Earth’s gravity field. We propose a new data processing method which makes use of the ‘average acceleration’ approach to gravity field modelling. In this method, satellite accelerations are directly derived from GPS carrier phase measurements with an epoch-differenced scheme. As a result, no ambiguity solutions are needed and the systematic errors that do not change much from epoch to epoch are largely eliminated. The GPS data collected by the Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) satellite mission are used to demonstrate the added value of the proposed method. An analysis of the residual accelerations shows that accelerations derived in this way are more precise, with noise being reduced by about 20 and 5% at the cross-track component and the other two components, respectively, as compared to those based on kinematic orbits. The accelerations obtained in this way allow the recovery of the gravity field to a slightly higher maximum degree compared to the solution based on kinematic orbits. Furthermore, the gravity field solution has an overall better performance. Errors in spherical harmonic coefficients are smaller, especially at low degrees. The cumulative geoid height error is reduced by about 15 and 5% up to degree 50 and 150, respectively. An analysis in the spatial domain shows that large errors along the geomagnetic equator, which are caused by a high electron density coupled with large short-term variations, are substantially reduced. Finally, the new method allows for a better observation of mass transport signals. In particular, sufficiently realistic signatures of regional mass anomalies in North America and south-west Africa are obtained.     

Surface gravitational field and topography changes induced by the Earth’s fluid core motions

Using a Love number formalism, the elastic deformations of the mantle and the mass redistribution gravitational potential within the Earth induced by the fluid pressure acting at the core–mantle boundary (CMB) are computed. This pressure field changes at a decadal time scale and may be estimated from observations of the surface magnetic field and its secular variation. First, using a spherical harmonic expansion, the poloidal and toroidal part of the fluid velocity field at the CMB for the last 40 years is computed, under the hypothesis of tangential geostrophy. Then the associated geostrophic pressure, whose order of magnitude is about 1000 Pa, is computed. The surface topography induced by this pressure field is computed using Love numbers, and is a few millimetres. The mass redistribution gravitational potential induced by these deformations and, in particular, the zonal components of the related surface gravitational potential perturbation (J₂, J₃ and J₄ coefficients), are calculated. Overall perturbations for the J₂ coefficient of about 10^–10, for J₃ of about 10^–11 and for J₄ are found of about 0.3×10^–11. Finally, these theoretical results are compared with recent observations of the decadal variation of J₂ from satellite laser ranging. Results concerning J₂ can be described as follows: first, they are one order of magnitude too small to explain the observed decadal variation of J₂ and, second, they show a significant linear trend over the last 40 years, whose rate of decrease amounts to 7% of the observed value.     

Conversion of time-varying Stokes coefficients into mass anomalies at the Earth’s surface considering the Earth’s oblateness

Time-varying Stokes coefficients estimated from GRACE satellite data are routinely converted into mass anomalies at the Earth’s surface with the expression proposed for that purpose by Wahr et al. (J Geophys Res 103(B12):30,205–30,229, 1998). However, the results obtained with it represent mass transport at the spherical surface of 6378 km radius. We show that the accuracy of such conversion may be insufficient, especially if the target area is located in a polar region and the signal-to-noise ratio is high. For instance, the peak values of mean linear trends in 2003–2015 estimated over Greenland and Amundsen Sea embayment of West Antarctica may be underestimated in this way by about 15%. As a solution, we propose an updated expression for the conversion of Stokes coefficients into mass anomalies. This expression is based on the assumptions that: (i) mass transport takes place at the reference ellipsoid and (ii) at each point of interest, the ellipsoidal surface is approximated by the sphere with a radius equal to the current radial distance from the Earth’s center (“locally spherical approximation”). The updated expression is nearly as simple as the traditionally used one but reduces the inaccuracies of the conversion procedure by an order of magnitude. In addition, we remind the reader that the conversion expressions are defined in spherical (geocentric) coordinates. We demonstrate that the difference between mass anomalies computed in spherical and ellipsoidal (geodetic) coordinates may not be negligible, so that a conversion of geodetic colatitudes into geocentric ones should not be omitted.     

Nonlinear diffusion filtering of data on the Earth’s surface

The paper deals with data filtering on closed surfaces using linear and nonlinear diffusion equations. We define a surface finite-volume method to approximate numerically parabolic partial differential equations on closed surfaces, namely on a sphere, ellipsoid or the Earth’s surface. The closed surface as a computational domain is approximated by a polyhedral surface created by planar triangles and we construct a dual co-volume grid. On the co-volumes we define a weak formulation of the problem by applying Green’s theorem to the Laplace–Beltrami operator. Then the finite-volume method is applied to discretize the weak formulation. Weak forms of elliptic operators are expressed through surface gradients. In our numerical scheme we use a piece-wise linear approximation of a solution in space and the backward Euler time discretization. Furthermore, we extend a linear diffusion on surface to the regularized surface Perona–Malik model. It represents a nonlinear diffusion equation, which at the same time reduces noise and preserves main edges and other details important for a correct interpretation of the real data. We present four numerical experiments. The first one has an illustrative character showing how an additive noise is filtered out from an artificial function defined on a sphere. Other three examples deal with the real geodetic data on the Earth’s surface, namely (i) we reduce a stripping noise from the GOCE satellite only geopotential model up to degree 240, (ii) we filter noise from the real GOCE measurements (the component $T_{zz})$ , and (iii) we reduce a stripping noise from the satellite only mean dynamic topography at oceans. In all experiments we focus on a comparison of the results obtained by both the linear and nonlinear models presenting advantages of the nonlinear diffusion.     

New Solution for the Earth’s Free Wobble and Its Geophysical Implications

In this paper, the theory of the free wobble of the triaxial Earth is developed and new conclusions are drawn: the Euler period should be actually expressed by the first kind of complete elliptic integral; the trace of the free polar motion is elliptic and the orientations of its semi-minor and major axes are approximately parallel to the Earth’s principal axes A and B, respectively. In addition, the present theory shows that there is a mechanism of frequency-amplitude modulation in the Chandler wobble, whi...     

Error analysis of the NGS’ surface gravity database

Are the National Geodetic Survey’s surface gravity data sufficient for supporting the computation of a 1 cm-accurate geoid? This paper attempts to answer this question by deriving a few measures of accuracy for this data and estimating their effects on the US geoid. We use a data set which comprises ${\sim }1.4$ million gravity observations collected in 1,489 surveys. Comparisons to GRACE-derived gravity and geoid are made to estimate the long-wavelength errors. Crossover analysis and $K$ -nearest neighbor predictions are used for estimating local gravity biases and high-frequency gravity errors, and the corresponding geoid biases and high-frequency geoid errors are evaluated. Results indicate that 244 of all 1,489 surface gravity surveys have significant biases ${>}2$  mGal, with geoid implications that reach 20 cm. Some of the biased surveys are large enough in horizontal extent to be reliably corrected by satellite-derived gravity models, but many others are not. In addition, the results suggest that the data are contaminated by high-frequency errors with an RMS of ${\sim }2.2$  mGal. This causes high-frequency geoid errors of a few centimeters in and to the west of the Rocky Mountains and in the Appalachians and a few millimeters or less everywhere else. Finally, long-wavelength ( ${>}3^{\circ }$ ) surface gravity errors on the sub-mGal level but with large horizontal extent are found. All of the south and southeast of the USA is biased by +0.3 to +0.8 mGal and the Rocky Mountains by $-0.1$ to $-0.3$  mGal. These small but extensive gravity errors lead to long-wavelength geoid errors that reach 60 cm in the interior of the USA.     

Observation of the Earth’s nutation by the VLBI: how accurate is the geophysical signal

We compare nutation time series determined by several International VLBI Service for geodesy and astrometry (IVS) analysis centers. These series were made available through the International Earth Rotation and Reference Systems Service (IERS). We adjust the amplitudes of the main nutations, including the free motion associated with the free core nutation (FCN). Then, we discuss the results in terms of physics of the Earth’s interior. We find consistent FCN signals in all of the time series, and we provide corrections to IAU 2000A series for a number of nutation terms with realistic errors. It appears that the analysis configuration or the software packages used by each analysis center introduce an error comparable to the amplitude of the prominent corrections. We show that the inconsistencies between series have significant consequences on our understanding of the Earth’s deep interior, especially for the free inner core resonance: they induce an uncertainty on the FCN period of about 0.5 day, and on the free inner core nutation (FICN) period of more than 1000 days, comparable to the estimated period itself. Though the FCN parameters are not so much affected, a 100 % error shows up for the FICN parameters and prevents from geophysical conclusions.     

Using DORIS measurements for modeling the vertical total electron content of the Earth’s ionosphere

The Doppler orbitography and radiopositioning integrated by satellite (DORIS) system was originally developed for precise orbit determination of low Earth orbiting (LEO) satellites. Beyond that, it is highly qualified for modeling the distribution of electrons within the Earth’s ionosphere. It measures with two frequencies in L-band with a relative frequency ratio close to 5. Since the terrestrial ground beacons are distributed quite homogeneously and several LEOs are equipped with modern receivers, a good applicability for global vertical total electron content (VTEC) modeling can be expected. This paper investigates the capability of DORIS dual-frequency phase observations for deriving VTEC and the contribution of these data to global VTEC modeling. The DORIS preprocessing is performed similar to commonly used global navigation satellite systems (GNSS) preprocessing. However, the absolute DORIS VTEC level is taken from global ionospheric maps (GIM) provided by the International GNSS Service (IGS) as the DORIS data contain no absolute information. DORIS-derived VTEC values show good consistency with IGS GIMs with a RMS between 2 and 3 total electron content units (TECU) depending on solar activity which can be reduced to less than 2 TECU when using only observations with elevation angles higher than \(50^\circ \) . The combination of DORIS VTEC with data from other space-geodetic measurement techniques improves the accuracy of global VTEC models significantly. If DORIS VTEC data is used to update IGS GIMs, an improvement of up to 12  % can be achieved. The accuracy directly beneath the DORIS satellites’ ground-tracks ranges between 1.5 and 3.5 TECU assuming a precision of 2.5 TECU for altimeter-derived VTEC values which have been used for validation purposes.     

Effects of topographic–isostatic masses on gravitational functionals at the Earth’s surface and at airborne and satellite altitudes

Topographic–isostatic masses represent an important source of gravity field information, especially in the high-frequency band, even if the detailed mass-density distribution inside the topographic masses is unknown. If this information is used within a remove-restore procedure, then the instability problems in downward continuation of gravity observations from aircraft or satellite altitudes can be reduced. In this article, integral formulae are derived for determination of gravitational effects of topographic–isostatic masses on the first- and second-order derivatives of the gravitational potential for three topographic–isostatic models. The application of these formulas is useful for airborne gravimetry/gradiometry and satellite gravity gradiometry. The formulas are presented in spherical approximation by separating the 3D integration in an analytical integration in the radial direction and 2D integration over the mean sphere. Therefore, spherical volume elements can be considered as being approximated by mass-lines located at the centre of the discretization compartments (the mass of the tesseroid is condensed mathematically along its vertical axis). The errors of this approximation are investigated for the second-order derivatives of the topographic–isostatic gravitational potential in the vicinity of the Earth’s surface. The formulas are then applied to various scenarios of airborne gravimetry/gradiometry and satellite gradiometry. The components of the gravitational vector at aircraft altitudes of 4 and 10 km have been determined, as well as the gravitational tensor components at a satellite altitude of 250 km envisaged for the forthcoming GOCE (gravity field and steady-state ocean-circulation explorer) mission. The numerical computations are based on digital elevation models with a 5-arc-minute resolution for satellite gravity gradiometry and 1-arc-minute resolution for airborne gravity/gradiometry.     

On the influence of the ground track on the gravity field recovery from high–low satellite-to-satellite tracking missions: CHAMP monthly gravity field recovery using the energy balance approach revisited

In this paper, the influence of the ground track coverage on the quality of a monthly gravity field solution is investigated for the scenario of a high–low satellite- to-satellite tracking mission. Data from the CHAllenging Minisatellite Payload (champ) mission collected in the period April 2002 to February 2004 has been used to recover the gravity field to degree and order 70 on a monthly basis. The quality is primarily restricted by the accuracy of the instruments. Besides, champ passed through a 31/2 repeat mode three times during the period of interest resulting in an insufficient spatial sampling and a degraded solution. Contrary to the rule of thumb by Colombo (The global mapping of gravity with two satellites, Publications on Geodesy, vol 7(3), Netherlands Geodetic Commission, The Netherlands, 263 pp, 1984), see also Wagner (J Geod 80(2): 94–103, 2006), we found that the monthly solutions themselves could be recovered to about degree 30, not 15. In order to improve the monthly gravity solutions, two strategies have been developed: the restriction to a low degree, and the densification of the sampling by the introduction of additional sensitive measurements from contemporaneous satellite missions. The latter method is tested by combining the champ measurements with data from the Gravity Recovery And Climate Experiment (grace). Note that the two grace satellites are considered independent here, i.e. no use is made of the K-band ranging data. This way, we are able to almost entirely remove the influence of the ground track leaving the accuracy of the instruments as the primary restriction on the quality of a monthly solution. These findings are especially interesting for the upcoming swarm-mission since it will consist of a similar configuration as the combined champ and (grace) missions.     

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.     

Modeling of earth’s gravity fields visualization based on Quad Tree

The problems of the earth’s gravity fields’ visualization are both focus and puzzle currently. Aiming at multiresolution rendering, modeling of the Earth’s gravity fields’ data is discussed in the paper by using LOD algorithm based on Quad Tree. First, this paper employed the method of LOD based on Quad Tree to divide up the regional gravity anomaly data, introduced the combined node evaluation system that was composed of viewpoint related and roughness related systems, and then eliminated the T-cracks that appeared among the gravity anomaly data grids with different resolutions. The test results demonstrated that the gravity anomaly data grids’ rendering effects were living, and the computational power was low. Therefore, the proposed algorithm was a suitable method for modeling the gravity anomaly data and has potential applications in visualization of the earth’s gravity fields.     

Development of the second-order derivatives of the Earth’s potential in the local north-oriented reference frame in orthogonal series of modified spherical harmonics

The second-order derivatives of the Earth’s potential in the local north-oriented reference frame are expanded in series of modified spherical harmonics. Linear relations are derived between the spectral coefficients of these series and the spectrum of the geopotential. On the basis of these relations, recurrence procedures are developed for evaluating the geopotential coefficients from the spectrum of each derivative and, inversely, for simulating the latter from a known geopotential model. Very simple structure of the derived expressions for the derivatives is convenient for estimating the geopotential coefficients by the least-squares procedure, at a certain step of processing satellite gradiometry data. Due to the orthogonality of the new series, the quadrature formula approach can be also applied, which allows avoidance of aliasing errors caused by the series truncation. The spectral coefficients of the derivatives are evaluated on the basis of the derived relations from the geopotential models EGM96 and EIGEN-CG01C at a mean orbital sphere of the GOCE satellite. Various characteristics of the spectra are studied corresponding to the EGM96 model. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Theory and Methods of the Earth＇s Gravity Field Model Recovery from GOCE Data

研究GOCE卫星测量恢复地球重力场模型的理论与方法。论文的主要工作和创新点有： （1）建立扰动重力梯度张量各分量没有奇异性的详细计算模型,解决重力梯度张量Txx分量在两极地区计算的奇异性难题。