Glacial isostatic adjustment in Fennoscandia for a laterally heterogeneous earth

Glaciation and deglaciation in Fennoscandia during the last glacial cycles has significantly perturbed the Earth's equilibrium figure. Changes in the Earth's solid and geoidal surfaces due to external and internal mass redistributions are recorded in sequences of ancient coastlines, now either submerged or uplifted, and are still visible in observations of present‐day motions of the surface and glacially induced anomalies in the Earth's gravitational field. These observations become increasingly sophisticated with the availability of GPS measurements and new satellite gravity missions.
Observational evidence of the mass changes is widely used to constrain the radial viscosity structure of the Earth's mantle. However, lateral changes in earth model properties are usually not taken into account, as most global models of glacial isostatic adjustment assume radial symmetry for the earth model. This simplifying assumption contrasts with seismological evidence of significant lateral variations in the Earth's crust and upper mantle throughout the Fennoscandian region.
We compare predictions of glacial isostatic adjustment based on an ice model over the Fennoscandian region for the last glacial cycle for both radially symmetric and fully 3‐D earth models. Our results clearly reveal the importance of lateral variations in lithospheric thickness and asthenospheric viscosity for glacially induced model predictions. Relative sea‐level predictions can differ by up to 10–20 m, uplift rate predictions by 1–3 mm yr⁻¹ and free‐air gravity anomaly predictions by 2–4 mGal when a realistic 3‐D earth structure as proposed by seismic modelling is taken into account.     

Inference of mantle viscosity from GRACE and relative sea level data

Gravity Recovery And Climate Experiment (GRACE) satellite observations of secular changes in gravity near Hudson Bay, and geological measurements of relative sea level (RSL) changes over the last 10 000 yr in the same region, are used in a Monte Carlo inversion to infer-mantle viscosity structure. The GRACE secular change in gravity shows a significant positive anomaly over a broad region (>3000 km) near Hudson Bay with a maximum of ∼2.5 μGal yr⁻¹ slightly west of Hudson Bay. The pattern of this anomaly is remarkably consistent with that predicted for postglacial rebound using the ICE-5G deglaciation history, strongly suggesting a postglacial rebound origin for the gravity change. We find that the GRACE and RSL data are insensitive to mantle viscosity below 1800 km depth, a conclusion similar to that from previous studies that used only RSL data. For a mantle with homogeneous viscosity, the GRACE and RSL data require a viscosity between  1.4 × 10²¹  and  2.3 × 10²¹  Pa s. An inversion for two mantle viscosity layers separated at a depth of 670 km, shows an ensemble of viscosity structures compatible with the data. While the lowest misfit occurs for upper- and lower-mantle viscosities of  5.3 × 10²⁰  and  2.3 × 10²¹  Pa s, respectively, a weaker upper mantle may be compensated by a stronger lower mantle, such that there exist other models that also provide a reasonable fit to the data. We find that the GRACE and RSL data used in this study cannot resolve more than two layers in the upper 1800 km of the mantle.     

The reduction of aliasing in gravity anomalies and geoid heights using digital terrain data

Observations of gravity can be aliased by virtue of the logistics involved in collecting these data in the field. For instance, gravity measurements are often made in more accessible lowland areas where there are roads and tracks, thus omitting areas of higher relief in between. The gravimetric determination of the geoid requires mean terrain-corrected free-air anomalies; however, anomalies based only on the observations in lowland regions are not necessarily representative of the true mean value over the topography. A five-stage approach is taken that uses a digital elevation model, which provides a more accurate representation of the topography than the gravity observation elevations, to reduce the unrepresentative sampling in the gravity observations. When using this approach with the Australian digital elevation model, the terrain-corrected free-air anomalies generated from the Australian gravity data base change by between 77.075 and −84.335 mgal (−0.193 mgal mean and 2.687 mgal standard deviation). Subsequent gravimetric geoid computations are used to illustrate the effect of aliasing in the Australian gravity data upon the geoid. The difference between 'aliased' and 'non-aliased' gravimetric geoid solutions varies by between 0.732 and −1.816 m (−0.058 m mean and 0.122 m standard deviation). Based on these conceptual arguments and numerical results, it is recommended that supplementary digital elevation information be included during the estimation of mean gravity anomalies prior to the computation of a gravimetric geoid model.     

An Analytic Model For A Mantle Plume Fed By A Boundary Layer

Summary. An approximate analytical solution for flow in a mantle plume of constant radius, viscosity, and density contrast is obtained in cylindrical coordinates. the differential equations for vertical velocity of the mantle surrounding the plume and for topography are homologous to the equation for flexure of an elastic plate. Although the model is too simple to be fully applicable to the Earth, one can conclude that the vertical velocity in the mantle changes significantly away from plumes, that the viscosity of the plume is important for controlling flow rate, and that the long-wavelength geoid anomalies are sensitive to the viscosity of the surrounding mantle. the first induced upwelling away from a plume is quite weak and unlikely to control the spacing of plumes.     

Impact of the formation and ablation of Antarctic ice sheet on the global geoid and sea level

１ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅＡｎｔａｒｃｔｉｃｉｃｅｓｈｅｔｎｅａｒｌｙｏｃｕｐｉｅｓ９０％ｏｆｔｈｅｇｌｏｂａｌｏｎｅｓ,ｔｈｅｆｏｒｍａｔｉｏｎａｎｄａｂｌａｔｉｏｎｏｆｗｈｉｃｈｈａｖｅａｓｔｒｏｎｇｉｍｐａｃｔｏｎｔｈｅｇｌｏｂａｌｇｅｏ...     

Detailed 1×1° gravimetric Indian Ocean geoid and comparison with GEOS-3 radar altimeter geoid profiles

Summary. A new set of 1×1° mean free-air anomalies in the Indian Ocean is determined on the basis of previously published free-air anomaly maps (Talwani & Kahle) and the most recent Lamont surface ship gravity measurements. The data are then used to compute a (total) 1×1° gravimetric Indian Ocean geoid. The computation is carried out by combining the Goddard Space Flight Center (GSFC) GEM-6 geoid and a difference geoid that corresponds to the differences between the set of 1×1° surface gravity values and the GEM-6 gravity anomalies. The difference geoid is highest over the Madagascar Ridge (+ 20 m) and lowest over the Timor Trough (-30 m). The total geoid is compared with GEOS-3 radar altimeter derived geoid profiles and geophysical implications are discussed.     

Viscous gravitational relaxation

Summary This paper is concerned with a detailed examination of the response of Maxwell models of the planet to surface mass loads. Particular attention is devoted to an examination of the factors which determine the isostatic response since the understanding of this response is crucial in a number of different geodynamic problems. One particular example which we discuss in detail is concerned with the prediction of free air gravity anomalies produced by large-scale deglaciation events. Using the methods developed here we are able to provide the first direct assessment of the importance of initial isostatic disequilibrium on the observed relative sea-level variations and free air gravity anomalies forced by the melting of the Laurentide ice sheet. We are therefore able to estimate the extent to which such initial disequilibrium might influence the inference of mantle viscosity from isostatic adjustment data. Our calculations establish that free air gravity data, although they are sensitive to the degree of initial disequilibrium, provide an extremely high quality constraint upon the viscosity of the lower mantle.     

On the gravitational potential and field anomalies due to thin mass layers

Summary. The gravitational potential and field anomalies for thin mass layers are derived using the technique of matched asymptotic expansions. An inner solution is obtained using an expansion in powers of the thickness and it is shown that the outer solution is given by a surface distribution of mass sources and dipoles. Coefficients are evaluated by matching the inner expansion of the outer solution with the outer expansion of the inner solution. The leading term in the inner expansion for the normal gravitational field gives the Bouguer formula. The leading term in the expansion for the gravitational potential gives an expression for the perturbation to the geoid. The predictions given by this term are compared with measurements by satellite altimetry. The second-order terms in the expansion for the gravitational field are required to predict the gravity anomaly at a continental margin. The results are compared with observations.     

Secular variation of the poloidal magnetic field at the core–mantle boundary

A region of enhanced conductivity at the base of the mantle is modelled by an infinitesimally thin sheet of uniform effective conductance adjacent to the core–mantle boundary. Currents induced in this sheet by the temporally varying magnetic field produced by the geodynamo give rise to a discontinuity in the horizontal components of the poloidal magnetic field on crossing the sheet, while the radial component is continuous across the sheet. Treating the rest of the mantle as an insulator, the horizontal components of the poloidal magnetic field and their secular variation at the top of the core are determined from geomagnetic field, secular variation and secular acceleration models. It is seen that for an assumed effective conductance of the sheet of 10⁸  S, which may be not unrealistic, the changes produced in the horizontal components of the poloidal field at the top of the core are usually ≤10 per cent, but corrections to the secular variation in these components at the top of the core are typically 40 per cent, which is greater than the differences that exist between different secular variation models for the same epoch. Given the assumption that all the conductivity of the mantle is concentrated into a thin shell, the present method is not restricted to a weakly conducting mantle. Results obtained are compared with perturbation solutions.     

The Canary Islands swell: a coherence analysis of bathymetry and gravity

The Canary Archipelago is an intraplate volcanic chain, located near the West African continental margin, emplaced on old oceanic lithosphere of Jurassic age, with an extended volcanic activity since Middle Miocene. The adjacent seafloor does not show the broad oceanic swell usually observed in hotspot-generated oceanic islands. However, the observation of a noticeable depth anomaly in the basement west of the Canaries might indicate that the swell is masked by a thick sedimentary cover and the influence of the Canarian volcanism. We use a spectral approach, based on coherence analysis, to determine the swell and its compensation mechanism. The coherence between gravity and topography indicates that the swell is caused by a subsurface load correlated with the surface volcanic load. The residual gravity/geoid anomaly indicates that the subsurface load extends 600 km SSW and 800 km N and NNE of the islands. We used computed depth anomalies from available deep seismic profiles to constrain the extent and amplitude of the basement uplift caused by a relatively low-density anomaly within the lithospheric mantle, and coherence analysis to constrain the elastic thickness of the lithosphere ( T_e ) and the compensation depth of the swell. Depth anomalies and coherence are well simulated with T_e =28–36 km, compensation depth of 40–65 km, and a negative density contrast within the lithosphere of ∼33 kg m⁻³. The density contrast corresponds to a temperature increment of ∼325°C, which we interpret to be partially maintained by a low-viscosity convective layer in the lowermost lithosphere, and which probably involves the shallower parts of the asthenosphere. This interpretation does not require a significant rejuvenation of the mechanical properties of the lithosphere.     

Mantle convection under simulated plates: effects of heating modes and ridge and trench migration, and implications for the core—mantle boundary, bathymetry, the geoid and Benioff zones

Summary. Numerical convection models are presented in which plates are simulated by imposing piecewise constant horizontal velocities on the upper boundary. A 4 × 1 box of constant viscosity fluid and two-dimensional (2-D) flow is assumed. Four heating modes are compared: the four combinations of internal or bottom heating and prescribed bottom temperature or heat flux. The case with internal heating and an isothermal base is relevant to lower mantle or whole mantle convection, and it yields a lower thermal boundary layer which is laterally variable and can be locally reversed, corresponding to heat flowing back into the core locally. When scaled to the whole mantle, the surface deflections and gravity and geoid perturbations calculated from the models are comparable to those observed at the Earth's surface. For models with migrating ridges and trenches, the flow structure lags well behind the changing surface 'plate'configurations. This may help to explain the poor correlation between the main geoid features and plate boundaries. Trench migration substantially affects the dip of the cool descending fluid because of induced horizontal shear in the vicinity of the trench. Such shear is small for whole mantle convection, but is large for upper mantle convection, and would probably result in the Tonga Benioff zone dipping to the SE, opposite to the observed dip, for the case of upper mantle convection.     

The onset of Pleistocene glaciation in the Barents Sea: implications for glacial isostatic adjustment

In this paper the effect of a delayed onset of glaciation in the Barents Sea on glacial isostatic adjustment is investigated. The model calculations solve the sea-level equation governing the total mass redistributions associated with the last glaciation cycle on a spherically symmetric, linear, Maxwell viscoelastic earth for two different scenarios for the growth phase of the Barents Sea ice sheet. In the first ice model a linear growing history is used for the Barents Sea ice sheet, which closely relates its development to the build-up of other major Late Pleistocene ice sheets. In the second ice model the accumulation of the Barents Sea ice sheet is restricted to the last 6 ka prior to the last glacial maximum.
The calculations predict relative sea levels, present-day radial velocities, and gravity anomalies for the area formerly covered by the Weichselian ice sheet. The results show that observed relative sea levels in the Barents Sea are appropriate for distinguishing between the different glaciation histories. In particular, present-day observables such as the free-air gravity anomaly over the Barents Sea, and the present-day radial velocities are sensitive to changes in the glaciation history on this scale.
A palaeobathymetry derived from relative sea-level predictions before the last glacial maximum based on the second ice model essentially agrees with a palaeobathymetry derived by Lambeck (1995). The additional emerged areas provide centres for the build-up of an ice sheet and thus support the theory of Hald, Danielsen & Lorentzen (1990) and Mangerud et al. (1992) that the Barents Sea was an essentially marine environment shortly before the last glacial maximum.     

Surface vertical displacements, potential perturbations andgravity changes of a viscoelastic earth model induced by internal point dislocations

Using the viscoelastic correspondence principle, we utilize the surface coseismic spheroidal deformation fields (i.e. vertical displacements, potential perturbations and gravity changes) of SNREI earth models caused by four typical types of point dislocation, derived by Sun & Okubo (1993 ), to deduce the fundamental formulas for spheroidal fields relevant to viscoelastic earth models. In computations, we employ a strike-slip dislocation on a vertical plane buried at the bottom of the lithosphere to estimate the maximal viscous relaxation responses to this kind of source that possibly exist on the surface of the earth. We take the seismic moment as 10²²  N  m, which is characteristic of an average large earthquake. The numerical results demonstrate that, if we take the viscosity as 10¹⁹  Pa  s in the asthenosphere, and 10²¹  Pa  s in the other mantle layers, the rates of surface vertical displacements and gravity changes within about 2.5° for the 10 postseismic years are respectively 1.5–8.1  cm  yr⁻¹ and 4.0–14.9  μgal  yr⁻¹ : the viscous relaxation for this mantle viscosity profile proceeds much faster than for a constant mantle viscosity of 10²¹  Pa  s.     

A technique for computing geoid height anomalies over two-dimensional earth models

Summary A technique is presented for calculating geoid height anomalies over two-dimensional models of Earth structure. The method consists of convolving gravity anomalies over the structure with filters which take into account the finite size of the structure in the third dimension and the curvature of the Earth. Similar filters are also developed for a flat earth case. The method is applied to a sea-surface gravity profile crossing the Tonga-Kermadec trench and is found to give good agreement with a Geos-3 radar altimetry profile in the same region. The example demonstrates that introducing arbitrary offsets in computing gravity anomalies can result in spurious long-wavelength effects in the computed geoid. Comparison of the results obtained using flat earth and spherical earth filters suggests that the effects of the curvature of the Earth only become significant for wavelengths in the gravity field greater than about 1000 km.     

Structural effects of the crust on the geoid modelled using deep seismic sounding interpretations

According to the theory of isostasy, the Earth has a tendency to deform its surface in order to reach an equilibrium state. The land-uplift phenomenon in the area of the Fennoscandian Shield is thought to be a process of this kind. The geoid, as an equipotential surface of the Earth's gravity field, contains information on how much the Earth's surface departs from the equilibrium state. In order to study the isostatic process through geoidal undulations, the structural effects of the crust on the geoid have to be investigated.
  The structure of the crust of the Fennoscandian Shield has been extensively explored by means of deep seismic sounding (DSS). The data obtained from DSS are used to construct a 3-D seismic-velocity structure model of the area's crust. The velocity model is converted to a 3-D density model using the empirical relationship that holds between seismic velocities and crustal mass densities. Structural effects are then estimated from the 3-D density model.
  The structural effects computed from the crustal model show that the mass deficiency of the crust in Fennoscandia has caused a geoidal depression twice as deep as that observed from the gravimetric geoid. It proves again that the crust has been isostatically compensated by the upper mantle. In other words, an anomalously high-density upper mantle must exist beneath Fennoscandia.     

Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity

The GRACE satellite mission has been measuring the Earth's gravity field and its temporal variations since 2002 April. Although these variations are mainly due to mass transfer within the geofluid envelops, they also result from mass displacements associated with phenomena including glacial isostatic adjustment and earthquakes. However, these last contributions are difficult to isolate because of the presence of noise and of geofluid signals, and because of GRACE's coarse spatial resolution (>400 km half-wavelength). In this paper, we show that a wavelet analysis on the sphere helps to retrieve earthquake signatures from GRACE geoid products. Using a wavelet analysis of GRACE geoids products, we show that the geoid variations caused by the 2004 December ( M _w= 9.2) and 2005 March ( M _w= 8.7) Sumatra earthquakes can be detected. At GRACE resolution, the 2004 December earthquake produced a strong coseismic decrease of the gravity field in the Andaman Sea, followed by relaxation in the area affected by both the Andaman 2004 and the Nias 2005 earthquakes. We find two characteristic timescales for the relaxation, with a fast variation occurring in the vicinity of the Central Andaman ridge. We discuss our coseismic observations in terms of density changes of crustal and upper-mantle rocks, and of the vertical displacements in the Andaman Sea. We interpret the post-seismic signal in terms of the viscoelastic response of the Earth's mantle. The transient component of the relaxation may indicate the presence of hot, viscous material beneath the active Central Andaman Basin.     

Small-scale instabilities below the cooling oceanic lithosphere

Interpretation of satellite altimetry data as well as ship bathymetry data revealed strongly elongated anomalies roughly perpendicular to the mid-ocean ridges in the Indian and east Pacific oceans. A spectral analysis of gravity altimetry data along profiles parallel to the East Pacific Rise indicated wavelengths of about 150–180  km close to the ridge and about 250  km further away. A simple model of Rayleigh–Taylor instabilities developing at the base of the cooling lithosphere is discussed and applied to the data. By considering thermal diffusion and comparing Rayleigh–Taylor growth rates to the velocity of the thermal front in the cooling lithosphere, we are able to explain the observed anomalies by instabilities developing below the lithosphere in a layer with a viscosity of about 10¹⁹  Pa  s above an asthenospheric layer with a viscosity reduction of 2–3 orders of magnitude.     

Gravity prediction using collocation and taking known mass density anomalies into account

Summary. The anomalous (gravitational) potential of the Earth, T , is split in two parts, T= T ^C + T _M. Here T _M is a harmonic function generated by known mass density anomalies and T ^C =T-T _M. This function will also be a harmonic function, which therefore may be approximated using the method of collocation, based on known gravity anomalies or altimeter derived geoid undulations, for example. Gravity anomalies can then be predicted using the known linearized relationship between T and Δ g . This procedure may give a 40–50 per cent increase in the precision of the prediction results as compared to a procedure where mass density anomalies are not taken into account.     

The Free Air Geoid for Australia

Summary The free air geoid, which is the co-geoid obtained by the use of free air anomalies in Stokes' integral, is computed for Australia from available gravity data. The set of anomalies used to represent the outer zones had been obtained previously using a combined solution from satellite data and terrestrial gravimetry. The solutions so obtained for the free air geoid are compared with the astrogeodetic determination of the geoid on the Australian Geodetic Datum by Fischer and Slutsky and the accuracy of the comparisons is estimated.     

Gravity anomalies, subsidence history and the tectonic evolution of the Malay and Penyu Basins (offshore Peninsular Malaysia)

The tectonic subsidence and gravity anomalies in the Malay and Penyu Basins, offshore Peninsular Malaysia, were analysed to determine the isostatic compensation mechanism in order to investigate their origin. These continental extensional basins contain up to 14  km of sediment fill which implies that the crust had been thinned significantly during basin development. Our results suggest, however, that the tectonic subsidence in the basins cannot be explained simply by crustal thinning and Airy isostatic compensation.
The Malay and Penyu Basins are characterized by broad negative free-air gravity anomalies of between −20 and −30  mGal. To determine the cause of the anomaly, we modelled four gravity profiles across the basins using a method that combines two-dimensional flexural backstripping and gravity modelling techniques. We assumed a model of uniform lithospheric stretching and Airy isostasy in the analysis of tectonic subsidence. Our study shows that the basins are probably underlain by relatively thinned crust, indicating that some form of crustal stretching was involved. To explain the observed gravity anomalies, however, the Moho depth that we calculated based on the free-air gravity data is about 25% deeper than the Moho predicted by assuming Airy isostasy (Backstrip Moho). This suggests that the Airy model overestimates the compensation and that the basins are probably undercompensated isostatically. In other words, there is an extra amount of tectonic subsidence that is not compensated by crustal thinning, which has resulted in the discrepancy between the gravity-derived Moho and the Backstrip Moho. We attribute this uncompensated or anomalous tectonic subsidence to thin-skinned crustal extension that did not involve the mantle lithosphere. The Malay and Penyu Basins are interpreted therefore as basins that formed by a combination of whole-lithosphere stretching and thin-skinned crustal extension.