Isostatic response of the lithosphere beneath the Mozambique Ridge (SW Indian Ocean) and geodynamic implications

International audienceS U M M A R Y The SW Indian Ocean is characterized by the presence of several aseismic features. The Mozambique Ridge, an elongated feature lying roughly parallel to the SE coast of Africa, is by far the least known of those structures, mainly due to the scarcity of marine data. To date, the crustal nature and the origin of the ridge are still controversial points. Since knowledge of the origin of the Mozambique Ridge is important for a better understanding of the evolution of the SW Indian Ocean, the isostatic response of the lithosphere beneath the ridge is analysed in order to characterize its effective elastic thickness and the emplacement process of the feature. Two different approaches are applied, the direct computation of the geoid anomaly over the ridge, by means of a 2.5-dimensional method, and the computation of the admittances between the bathymetry and both the geoid and gravity anomalies. Both approaches point to a local isostatic response of the lithosphere. The crustal thickness beneath the Mozambique Ridge ranges from 17 to 30 km and the average density, from 2.78 to 2.7 X lo3 kg m-3, varying with locality

Mantle lateral variations and elastogravitational deformations – I. Numerical modelling

International audienceThe Earth response (deformation and gravity) to tides or to surface loads is traditionally computed assuming radial symmetry in stratified earth models, at the hydrostatic equilibrium. The present study aims at providing a new earth elastogravitational deformation model which accounts for the whole complexity of a more realistic earth. The model is based on a dynamically consistent equilibrium state which includes lateral variations in density and elastic parameters, and interface topographies. The deviation from the hydrostatic equilibrium has been taken into account as a first-order perturbation. We use a finite element method (spectral element method) and solve numerically the gravitoelasticity equations. As a validation application, we investigate the deformation of the Earth to surface loads. We first evaluate the classical loading Love numbers with a relative precision of about 0.3 per cent for PREM earth model. Then we assume an ellipsoidal homogeneous incom-pressible earth with hydrostatic pre-stresses. We investigate the impact of ellipticity on loading Love numbers analytically and numerically. We validate and discuss our numerical model. At periods greater than 1 hr, the solid earth is mainly deformed by luni-solar tides and by surface loads induced by different external fluid layers (ocean, atmosphere, continental hydrology, ice volumes). This work is devoted to the analytical and numerical development to compute the response of the Earth to such forcing. The body tides have been investigated since the 19th century. In 1862, Lord Kelvin (Sir William Thomson) made the first calculation of the elastic deformation of a homogeneous incompressible earth under the action of the tidal gravitational potential (Thomson 1862). Some years later, Love (1911) studied a compressible homogeneous earth model and showed that the tidal effects could be represented by a set of dimensionless numbers, the so-called Love numbers. Takeuchi (1950) obtained a first estimation of the Love numbers by a numerical integration of the equations using a reference earth model deduced from seismology. These results have been later extended (Smith 1974; Wahr 1981) to an ellipsoidal, rotating Earth with hydrostatic pre-stresses and a liquid core, and finally the effects of mantle anelasticity have been included (Wahr & Bergen 1986; Dehant 1987). In addition to tidal forces, mass changes in the atmosphere cause deformation and mass redistribution inside the planet. The Earth's response to such forcing involves both local and global surface motions and variations in the gravity field, which may be observed in geodetic experiments. These hydrological, atmospheric or oceanic effects on the Earth's gravity field are usually modelled for a spherical Earth with hydrostatic pre-stress (e.g. Farrell 1972; Wahr et al. 1998), generally identified to the preliminary reference earth model (PREM) developed by Dziewonski & Anderson (1981). However, the internal structure of the Earth is more complex than in a spherical non-rotating elastic isotropic (SNREI) earth model like PREM. Seismology and fluid dynamic studies show that the mantle presents heterogeneous structure induced by a thermochemical convection (Davaille 1999; Gu et al. 2001; Forte & Mitrovica 2001) and a bias from hydrostatic state. Large lateral heterogeneities have taken place on a million year timescale (Courtillot et al. 2003), like the two supposed superplumes under the Pacific and South Africa superswells, or like descending slabs. These aspects of the mantle structure are classically not taken into account in the deformation models. The elastogravitational deformations are presently observed with very high accuracy. The accuracy of superconducting gravimeter and of positioning techniques (GPS, VLBI) has seen a large improvement in the last decade. Moreover, the global gravity field will be of interest in the next 10 yr with the launch of the missions GRACE (in 2002) and GOCE (in 2007), which are dedicated to gravimetry and gradiometry 106

A new approach to computing accurate gravity time variations for a realistic earth model with lateral heterogeneities

International audienceWe have developed a new elasto-gravitational earth model able to take into account lateral variations, deviatoric pre-stresses and topographies. As a first application, we assume an el-lipsoidal earth with hydrostatic pre-stresses, and validate and discuss our numerical model by comparison with previous studies on the M 2 body tide. We then study the response of the ellipsoidal earth to zonal atmospheric loads, and find that global lateral variations within the Earth, such as ellipticity, have a weak impact (about 1 per cent) on the elasto-gravitational deformations induced by atmospheric loading. At low frequencies, the Earth is deformed mainly by luni-solar tides and by surface loads, including ocean, atmosphere, ice volumes and post-glacial rebound. In this work, we focus our attention on the Earth's body tides and atmospheric loadings. The most accepted Earth body-tide models presently deal with an ellipsoidal, rotating earth, containing a liquid core and an anelastic mantle with hydrostatic pre-stresses (Wahr 1981; Wahr & Bergen 1986). The Earth, however, is not an exact ellipsoid, but presents lateral variations and deviatoric pre-stresses: there are long-wavelength density anomalies within the mantle, as shown by geoid anomalies and tomography studies (e.g. Romanowicz & Gung 2002). Wang (1994) and Dehant et al. (1999) studied the influence of lateral heterogeneities on Earth tides and showed that this effect is small but not necessarily negligible. They did not, however, take into account possible deviatoric pre-stresses: these effects on the Earth's body tides are totally unknown. In addition to tidal forces, mass changes in the atmosphere also cause deformation and mass redistribution inside the planet, involving both local and global surface motions and variations in the gravity field, which may be observed in geodetic experiments. For several decades, satellite geodesy has provided information on the temporal variation of the Earth's geopotential, and especially on the low-degree zonal harmonics (J 2 , J 3. . .) (Gegout & Cazenave 1993), which are essentially controlled by surface loads. These hydrological , atmospheric or oceanic effects on the Earth's gravity field are usually modelled assuming a spherical earth with hydro-static pre-stress (e.g. Farrell 1972; Wahr et al. 1998). With the advent of the new generation of gravity measurements, one of the challenges of the coming decade will be to provide more realistic earth models that show the variation of gravity with time. In particular, global studies based on gravity data from satellites such as GRACE, GOCE, and future GRACE/GOCE follow-on ones require accurate body-tide deformation models. More realistic gravity variation models are also needed for local and ground measurements, particularly for the very accurate superconducting gravimeters and the associated gravimetric observatory network such as the Global Geodynamic Project (Crossley et al. 1999). The formalism developed to compute this elasto-gravitational model is usually based on spherical harmonic analysis. The addition of lateral variations leads to couplings between spherical harmonics , i.e. to a more complex formalism that requires a large numerical effort (e.g. Wang 1994; Plag et al. 1996). We develop here a new approach for a non-radially symmetrical earth model using a finite-element method known as the spectral element method. The efficiency of this method is less dependent on the shape of the lateral heterogeneities than the spherical harmonic method. Our method is therefore well adapted to studying the impact of global and local lateral variations on the Earth deformation. We solve the elasto-gravitational equations taking into consideration the lateral variations within the Earth by using a first-order perturbation theory (Smith 1974; Dahlen & Tromp 1998). This new model allows us to take into account lateral variations of density and rheological parameters, deviatoric pre-stresses and interface topography. In order to validate our calculations, we tackle a well-known problem: the impact of the hydrostatic ellipticity on the Earth body tides. An analytical solution for this problem can be derived for a simple model in which the earth is assumed to be homogeneous and incompressible. The gravitational potential and the vertical displacement on the surface of the deformed ellipsoid were first derived by Love (1911) and then corrected by Wang (1994). We have recently extended these analytical results to the tangential surface displacement (Greff-Lefftz et al. 2005). We first validate our model with our analytical solutions, and then compare our results wit

Precise location of unsurveyed seamounts in the Austral archipelago area using SEASAT data

SEASAT altimetric geoid data are used to detect uncharted seamounts in the Austral archipelago area. The various physical parameters which affect the geoid signature of a seamount are inspected to analyse their influence on the precision of the location. If the correct elastic thickness is assumed, the precision on the location is order 15 km. Ten previously unsurveyed seamounts have been located in the Austral archipelago. It appears that they are emplaced along two well-defined azimuths and that two parallel distinct volcanic chains form the Austral archipelag

Extracting low frequency climate signal from GRACE data

International audienceFor more than four years, the GRACE pair of satellites have been orbiting the Earth, monitoring the time variable mass distribution for scales ranging from regional to global. The GRACE data have been released for a broad scientific community and sets of gravity fields are available. This paper shows that there are evidences at interrannual time scales for the presence of ENSO signal in the data, strongly correlated with the hydrological mass distribution, and also similar to the expected hydrological signature associated with the ENSO cycle. This signal dominates, at global scale, the one associated with geodynamic sources

Characterization of geological boundaries using 1‐D wavelet transform on gravity data: Theory and application to the Himalayas

International audienceWe investigate the use of the continuous wavelet transform for gravity inversion. The wavelet transform operator has recently been introduced in the domain of potential fields both as a filtering and a source-analysis tool. Here we develop an inverse scheme in the wavelet domain , designed to recover the geometric characteristics of density heterogeneities described by simple-shaped sources. The 1-D analyzing wavelet we use associates the upward continuation operator and linear combinations of derivatives of any order. In the gravity case, we first demonstrate how to localize causative sources using simple geometric constructions. Both the upper part of the source and the whole source can be studied when considering low or high altitudes, respectively. The ho-mogeneity degree of the source is deduced without prior information and allows us to infer its shape. Introducing complex wavelets, we derive analytically the scaling behavior of the wavelet coefficients for the dyke and the step sources. The modulus term is used in an inversion procedure to recover the thickness of the source. The phase term provides its dip. This analysis is performed on gravity data we measured along a profile across the Himalayas in Nepal. Good agreement of our results with well-documented thrusting structures demonstrates the applicability of the method to real data. Also, deeper, less constrained structures are characterized

A broken plate beneath the North Baikal rift zone revealed by gravity modelling

International audienceWe modelled a 1200 km long gravimetric profile in the North Baikal rift to assess the mechanical behaviour of the lithosphere, using a numerical model that accounts for realistic brittle-elasto-ductile rheology. We use published seismicity and re-fraction data, a new 5'x7.5' free-air/Bouguer gravity and topography data set, and a detailed map of faults obtained from high resolution SPOT imagery. Analysis of the gravity field over the North Baikal rift zone indicates significant asymmetry of the mechanical processes governing the deformation of the diverging sides of the rift. These anomalies cmmot be explained by a conventional continuous plate undergoing extension beneath the rift zone, whereas a strong mechanical discontinuity (wedge shaped detachment zone beneath the rift axis) is able to reproduce observations. Such a discontinuous model provides a good fit to the gravity and crustal thickness data and explains the deep seismicity reported there

A comparison of surface fitting algorithms for geophysical data

Cet article présente les résultats d'une comparaison de différents algorithmes d'approximation de surface. Pour chacun de ces algorithmes (approximation polynomiale, combinaison spline-laplace, krigeage, approximation aux moindres carrés, méthode des éléments finis) la pertinence pour différents ensembles de données et les limites d'application sont discutée