Search

Moving-base gravity gradiometry has proven to be a convenient method to determine the Earth’s gravity field. The ESA mission GOCE (Gravity field and steady-state Ocean Circulation Explorer) has enabled to map the Earth gravity field and its gradients with a resolution of 80 km, leading to significant advances in physical oceanography and solid Earth physics. At smaller scales, airborne gravity gradiometry has been increasingly used during the past decade in mineral and hydrocarbon exploration. In both cases the sensitivity of gradiometers to the short wavelengths of the gravity field is of crucial interest. Here, we quantify and characterize the error on the gravity gradients estimated from measurements performed with a new instrument concept, called GREMLIT, for typical airborne conditions. GREMLIT is an ultra-sensitive planar gravitational gradiometer which consists in a planar acceleration gradiometer together with 3 gyroscopes. To conduct this error analysis, a simulation of a realistic airborne survey with GREMLIT is carried out. We first simulate realistic GREMLIT synthetic data, taking into account the acceleration gradiometer and gyroscope noises and biases and the variation of orientation of the measurement reference frame. Then, we estimate the gravity gradients from these data. Special attention is paid to the processing of the gyroscopes measurements whose accuracy is not commensurate with the ultra-sensitive gradiometer. We propose a method to calibrate the gyroscopes biases with a precision of the order $$10^{-8}$$ rad/s. In order to transform the tensor from the measurement frame to the local geodetic frame, we estimate the error induced when replacing the non-measured elements of the gravity gradient tensor by an a priori model. With the appropriate smoothing, we show that it is possible to achieve a precision better than 2E for an along-track spatial resolution of 2 km.

The estimation of the zero-height geopotential level of a local vertical datum (LVD) is a key task towards the connection of isolated physical height frames and their unification into a common vertical reference system. Such an estimate resolves, in principle, the ‘ambiguity’ of a traditional crust-fixed LVD by linking it with a particular equipotential surface of Earth’s gravity field under the presence of an external geopotential model. The aim of this paper is to study the estimation scheme that can be followed for solving the aforementioned problem based on the joint inversion of co-located GPS and leveling heights in conjunction with a fixed Earth gravity field model. Several case studies with real data are also presented that provide, for the first time, precise estimates of the LVD offsets for a number of Hellenic islands across the Aegean and Ionian Sea.

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.

Since its launch in 2002, the Gravity Recovery and Climate Experiment (GRACE) mission has been providing measurements of the time-varying Earth gravity field. The GRACE mission architecture includes two satellites in near-circular, near-polar orbits separated in the along-track direction by approximately 220 km (e.g. collinear). A microwave ranging instrument measures changes in the distance between the spacecraft, while accelerometers on each spacecraft are used to measure changes in distance due to non-gravitational forces. The fact that the satellites are in near-polar orbits coupled with the fact that the inter-satellite range measurements are directed in the along-track direction, contributes to longitudinal striping in the estimated gravity fields. This paper examines four candidate mission architectures for a future gravity recovery satellite mission to assess their potential in measuring the gravity field more accurately than GRACE. All satellites were assumed to have an improved measurement system, with an inter-satellite laser ranging instrument and a drag-free system for removal of non-gravitational accelerations. Four formations were studied: a two-satellite collinear pair similar to GRACE; a four-satellite architecture with two collinear pairs; a two-satellite cartwheel formation; and a four-satellite cartwheel formation. A cartwheel formation consists of satellites performing in-plane, relative elliptical motion about their geometric center, so that inter-satellite measurements are, at times, directed radially (e.g. parallel to the direction towards the center of the Earth) rather than along-track. Radial measurements, unlike along-track measurements, have equal sensitivity to mass distribution in all directions along the Earth’s surface and can lead to higher spatial resolution in the derived gravity field. The ability of each architecture to recover the gravity field was evaluated using numerical simulations performed with JPL’s GIPSY-OASIS software package. Thirty days of data were used to estimate gravity fields complete to degree and order 60. Evaluations were done for 250 and 400 km nominal orbit altitudes. The sensitivity of the recovered gravity field to under-sampled effects was assessed using simulated errors in atmospheric/ocean dealiasing (AOD) models. Results showed the gravity field errors associated with the four-satellite cartwheel formation were approximately one order of magnitude lower than the collinear satellite pair when only measurement system errors were included. When short-period AOD model errors were introduced, the gravity field errors for each formation were approximately the same. The cartwheel formations eliminated most of the longitudinal striping seen in the gravity field errors. A covariance analysis showed the error spectrum of the cartwheel formations to be lower and more isotropic than that of the collinear formations.

International audience; DORIS (Détermination d'Orbite et Radiopositionnement Intégrés par Satellite) is a system used for precise orbit determination (POD) and ground-station positioning. It has been implemented on-board various satellites: the SPOT (Système pour l'Observation de la Terre) remote sensing satellites SPOT-2, SPOT-3, SPOT-4, SPOT-5, TOPEX/Poseidon and more recently on its successors Jason-1 and ENVISAT. DORIS is also a terrestrial positioning system that has found many applications in geophysics and geodesy; in particular, it contributes to the realization of the International Terrestrial Reference Frame, ITRF2000 and the forthcoming ITRF2005. Although not its primary objective, DORIS can bring information on Earth orientation monitoring, mainly polar motion and length of day (LOD) variations that complement other astrogeodetic techniques. In this paper, we have analyzed various recent polar motion solutions derived from independent analysis centers using different software packages and applying various analysis strategies. Comparisons of these solutions to the International Earth Rotation and Reference Systems Service (IERS) C04 solution are performed. Depending on the solutions, the accuracy of DORIS polar components are in the range of 0.5 1 mas corresponding to a few centimeters on the Earth's surface. This is approximately ten times larger than results derived from GPS, which are typically 0.06 mas in both components. This does not allow DORIS results to be taken into account in the IERS EOP combinations. A gain in the precision could come from technical improvements to the DORIS system, in addition to improvement of the orbit, tropospheric, ionospheric and Earth gravity field modeling.

A synthetic [simulated] Earth gravity model (SEGM) of the geoid, gravity and topography has been constructed over Australia specifically for validating regional gravimetric geoid determination theories, techniques and computer software. This regional high-resolution (1-arc-min by 1-arc-min) Australian SEGM (AusSEGM) is a combined source and effect model. The long-wavelength effect part (up to and including spherical harmonic degree and order 360) is taken from an assumed errorless EGM96 global geopotential model. Using forward modelling via numerical Newtonian integration, the short-wavelength source part is computed from a high-resolution (3-arc-sec by 3-arc-sec) synthetic digital elevation model (SDEM), which is a fractal surface based on the GLOBE v1 DEM. All topographic masses are modelled with a constant mass-density of 2,670 kg/m3. Based on these input data, gravity values on the synthetic topography (on a grid and at arbitrarily distributed discrete points) and consistent geoidal heights at regular 1-arc-min geographical grid nodes have been computed. The precision of the synthetic gravity and geoid data (after a first iteration) is estimated to be better than 30 μ Gal and 3 mm, respectively, which reduces to 1 μ Gal and 1 mm after a second iteration. The second iteration accounts for the changes in the geoid due to the superposed synthetic topographic mass distribution. The first iteration of AusSEGM is compared with Australian gravity and GPS-levelling data to verify that it gives a realistic representation of the Earth’s gravity field. As a by-product of this comparison, AusSEGM gives further evidence of the north–south-trending error in the Australian Height Datum. The freely available AusSEGM-derived gravity and SDEM data, included as Electronic Supplementary Material (ESM) with this paper, can be used to compute a geoid model that, if correct, will agree to in 3 mm with the AusSEGM geoidal heights, thus offering independent verification of theories and numerical techniques used for regional geoid modelling.

The remove-compute-restore (RCR) technique is the most well known method for regional gravimetric geoid determination today. Its basic theory is the first-order approximation of either Molodensky’s method for quasi-geoid determination or the classical geoid modelling by Helmert’s second method of condensing the topography onto the geoid. Although the basic approximate formulae do not meet today’s demands for a 1-cm geoid, it is sometimes assumed that the removal of the less precise long-wavelength terrestrial gravity anomaly field from Stokes’s integral by utilising a higher-order reference field represented by a more precise Earth gravity model (EGM) and the restoration of the EGM as a low-degree geoid contribution will produce a geoid model of the desired accuracy. Further improvement is achieved also by removing and restoring a residual topographic effect, which favourably smoothes the gravity anomaly to be integrated in Stokes’s formula. However, it is shown here that the RCR technique fails to tune down the long-wavelength gravity signal from the terrestrial data, and the EGM actually only reduces, in a non-optimised way, the truncation error committed by limiting the Stokes integration to a small region around the computation point. Hence, in order to take full advantage of a precise EGM, especially one from new dedicated satellite gravimetry, Stokes’s kernel must be modified in a suitable way to match the errors of terrestrial gravity, EGM and truncation. In addition, topographic, atmospheric and ellipsoidal effects must be carefully applied.

One important application of an Earth Gravity Model (EGM) is to determine the geoid. Since an EGM is represented by an external-type series of spherical harmonics, a biased geoid model is obtained when the EGM is applied inside the masses in continental regions. In order to convert the downward-continued height anomaly to the corresponding geoid undulation, a correction has to be applied for the analytical continuation bias of the geoid height. This technique is here called the geoid bias method. A correction for the geoid bias can also be utilised when an EGM is combined with terrestrial gravity data, using the combined approach to topographic corrections. The geoid bias can be computed either by a strict integral formula, or by means of one or more terms in a binomial expansion. The accuracy of the lowest binomial terms is studied numerically. It is concluded that the first term (of power H2) can be used with high accuracy up to degree 360 everywhere on Earth. If very high mountains are disregarded, then the use of the H2 term can be extended up to maximum degrees as high as 1800. It is also shown that the geoid bias method is practically equal to the technique applied by Rapp, which utilises the quasigeoid-to-geoid separation. Another objective is to carefully consider how the combined approach to topographic corrections should be interpreted. This includes investigations of how the above-mentioned H2 term should be computed, as well as how it can be improved by a correction for the residual geoid bias. It is concluded that the computation of the combined topographic effect is efficient in the case that the residual geoid bias can be neglected, since the computation of the latter is very time consuming. It is nevertheless important to be able to compute the residual bias for individual stations. For reasonable maximum degrees, this can be used to check the quality of the H2 approximation in different situations.

A new isostatic model for the Earth's gravity field is presented based on a simple hypothesis of layers approximating constant density contrasts. The spherical layer distribution used to describe the hydrostatic equilibrium of the Earth's masses leads to a new set of spherical harmonic coefficients for the gravitational potential. First attempts to quantify the information content of these coefficients led to the outcome that they seem to explain the observed gravity field for a certain wavelength band, while they are insufficient for short and very long wavelengths. A synthesis of the derived coefficients over specific degree ranges provided a computation of band-limited geoid undulations on a global scale. The association of these potential quantities with known tectonic structures, such as the topography of the core-mantle boundary, strengthens the belief that the interpretation of Earth gravity models, especially those arising from global digital elevation models, should be considered in close relation with deep-Earth structure.

Today the combination of Stokes’ formula and an Earth gravity model (EGM) for geoid determination is a standard procedure. However, the method of modifying Stokes’ formula varies from author to author, and numerous methods of modification exist. Most methods modify Stokes’ kernel, but the most widely applied method, the remove compute restore technique, removes the EGM from the gravity anomaly to attain a residual gravity anomaly under Stokes’ integral, and at least one known method modifies both Stokes’ kernel and the gravity anomaly. A general model for modifying Stokes’ formula is presented; it includes most of the well-known techniques of modification as special cases. By assuming that the error spectra of the gravity anomalies and the EGM are known, the optimum model of modification is derived based on the least-squares principle. This solution minimizes the expected mean square error (MSE) of all possible solutions of the general geoid model. A practical formula for estimating the MSE is also presented. The power of the optimum method is demonstrated in two special cases.

In a modern application of Stokes’ formula for geoid determination, regional terrestrial gravity is combined with long-wavelength gravity information supplied by an Earth gravity model. Usually, several corrections must be added to gravity to be consistent with Stokes’ formula. In contrast, here all such corrections are applied directly to the approximate geoid height determined from the surface gravity anomalies. In this way, a more efficient workload is obtained. As an example, in applications of the direct and first and second indirect topographic effects significant long-wavelength contributions must be considered, all of which are time consuming to compute. By adding all three effects to produce a combined geoid effect, these long-wavelength features largely cancel. The computational scheme, including two least squares modifications of Stokes’ formula, is outlined, and the specific advantages of this technique, compared to traditional gravity reduction prior to Stokes’ integration, are summarised in the conclusions and final remarks.

In precise geoid determination by Stokes formula, direct and primary and secondary indirect terrain effects are applied for removing and restoring the terrain masses. We use Helmert's second condensation method to derive the sum of these effects, together called the total terrain effect for geoid. We develop the total terrain effect to third power of elevation H in the original Stokes formula, Earth gravity model and modified Stokes formula. It is shown that the original Stokes formula, Earth gravity model and modified Stokes formula all theoretically experience different total terrain effects. Numerical results indicate that the total terrain effect is very significant for moderate topographies and mountainous regions. Absolute global mean values of 5–10 cm can be reached for harmonic expansions of the terrain to degree and order 360. In another experiment, we conclude that the most important part of the total terrain effect is the contribution from the second power of H, while the contribution from the third power term is within 9 cm.

A variety of mathematical models exist with which to combine long-wavelength geoid information from an Earth gravity model (EGM) with terrestrial gravity data. As new data become available, such as an extended or refined EGM and/or some variance–covariance information from terrestrial gravity data, this new information can be combined with a previous geoid model through a so-called hybrid estimator. The hybrid estimator is defined as a two-step estimator, which provides a solution that is a compromise between computational labour and matching of the entire data set. Some commonly used geoid modifications are compared, and several hybrid estimators are introduced. It is concluded that, in general, hybrid estimators linking the new information to old data by the least squares principle should be preferred to deterministic combinations.