Gravitational viscoelastic relaxation of eccentrically nested spheres

We present a semi-analytical solution to the 2-D forward modelling of viscoelastic relaxation in a heterogeneous model consisting of eccentrically nested spheres. Several numerical methods for 2-D and 3-D viscoelastic relaxation modelling have been applied recently, including finite-element and spectral-finite-difference schemes. The present semi-analytical approach provides a model response against which more general numerical algorithms can be validated. The eccentrically nested sphere solution has been tested by comparing it with the analytical solutions for viscoelastic relaxation in a homogeneous sphere and in two concentrically nested spheres, and good agreement was obtained.     

Postglacial rebound and lateral viscosity variations: a semi-analytical approach based on a spherical model with Maxwell rheology

A spectral method is employed to study the response to surface loads of a Maxwell earth including lateral viscosity variations. In particular, we focus on the effects of lithospheric cratons on the long-wavelength time-dependent displacement field for simple earth models. The viscosity contrast of the craton with respect to the surrounding mantle is kept fixed, whereas its thickness is allowed to vary. We show that the long-wavelength vertical displacement is not greatly affected by the presence of a lithospheric craton, while the tangential displacement is severely modified for the case of a homogeneous mantle. With increasing harmonic degree and thickness of the craton, the load-deformation coefficients deviate from those pertaining to a homogeneous mantle with a viscosity of 10²¹ Pa s. These deviations are particularly enhanced on timescales larger than a few hundred years. These findings indicate that the interpretation of the viscosity structure of the mantle inferred from postglacial rebound signatures based on radially stratified models is affected by the presence of lateral viscosity variations.     

Analytical approach for the toroidal relaxation of viscoelastic earth

This paper is concerned with post-seismic toroidal deformation in a spherically symmetric, non-rotating, linear-viscoelastic, isotropic Maxwell earth model. Analytical expressions for characteristic relaxation times and relaxation strengths are found for viscoelastic toroidal deformation, associated with surface tangential stress, when there are two to five layers between the core–mantle boundary and Earth's surface. The multilayered models can include lithosphere, asthenosphere, upper and lower mantles and even low-viscosity ductile layer in the lithosphere. The analytical approach is self-consistent in that the Heaviside isostatic solution agrees with fluid limit. The analytical solution can be used for high-precision simulation of the toroidal relaxation in five-layer earths and the results can also be considered as a benchmark for numerical methods. Analytical solution gives only stable decaying modes—unstable mode, conjugate complex mode and modes of relevant poles with orders larger than 1, are all excluded, and the total number of modes is found to be just the number of viscoelastic layers between the core–mantle boundary and Earth's surface—however, any elastic layer between two viscoelastic layers is also counted. This confirms previous finding where numerical method (i.e. propagator matrix method) is used. We have studied the relaxation times of a lot of models and found the propagator matrix method to agree very well with those from analytical results. In addition, the asthenosphere and lithospheric ductile layer are found to have large effects on the amplitude of post-seismic deformation. This also confirms the findings of previous works.     

Some remarks about the degree-one deformation of the Earth

The degree-one deformation of the Earth (and the induced discrepancy between the figure centre and the mass centre of the Earth) is computed using a theoretical approach (Love numbers formalism) at short timescales (where the Earth has an elastic behaviour) as well as at long timescales (where the Earth has a viscoelastic or quasi-fluid behaviour). For a Maxwell model of rheology, the degree-one relaxation modes associated with the viscoelastic Love numbers have been investigated: the M_o mode does not exist and there is only one transition mode (instead of two) generated by a viscosity discontinuity.
The translations at each interface of the incompressible layers of the earth model [surface, 670 km depth discontinuity, core-mantle boundary (CMB) and inner-core boundary (ICB)] are computed. They are elastic with an order of magnitude of about 1 mm when the excitation source is the atmospheric continental loading or a magnetic pressure acting at the CMB. They are viscoelastic when the earth is submitted to Pleistocene deglaciation, with an order of magnitude of about 1 m. In a quasi-fluid approximation (Newtonian fluid) because of the mantle density heterogeneity their order of magnitude is about 100 m (except for the ICB, which is in quasi-hydrostatic equilibrium at this timescale).     

First-order asymptotics for the eigenfrequencies of the Earth and application to the retrieval of large-scale lateral variations of structure

Summary . We present a variety of examples, showing systematic fluctuations as a function of angular order of measured eigenfrequencies for given normal modes of the Earth. The data are single station measurements from the GEOSCOPE network. Such fluctuations are attributed to departures from the lowest order asymptotic expression of the geometrical optics approximation. We derive first-order asymptotic expressions for the location parameter for all three components of the Earth's motion, by a method based on the stationary phase approximation and geometric relations on the unit sphere.
We illustrate the sensitivity of the fluctuations to the different parameters involved (source parameters, epicentral distance, laterally heterogeneous earth model) with synthetic examples corresponding to GEOSCOPE observations. Finally, we show the results of first attempts at inversion, which indicate that, when the fluctuations are taken into account, more accurate estimates of the great circle average eigenfrequencies can be obtained, and additional constraints put on the structure in the neighbourhood of the great circle.     

About the influence of pre-stress upon adiabatic perturbations of the Earth

Summary. In this paper we examine the influence of the state of stress in the equilibrium configuration of the Earth (i.e. the pre-stress) upon its adiabatic perturbations. The equations governing these perturbations to the first order (Woodhouse & Dahlen; Dahlen) are re-derived using a Lagrangian approach. Different expressions of the sesquilinear form associated to the elastic-gravitational operator are given. One of these provides a way to extend to hydrostatically pre-stressed solids the criterion of local stability given by Friedman & Schutz for uniformly rotating fluids. Then the propagation in the Earth of seismic wavefronts is considered. It is shown that the nature of these different wavefronts is entirely determined by the quadratic coefficients of the development of the specific internal energy variation, corresponding to isentropic evolution, with respect to the Lagrangian finite deformation tensor. Expressions for the velocities of the various waves are given as functions of incidence angle and pre-stress for orthotropic elastic material. In the particular case where the elastic parameters depend only on one coordinate of a curvilinear system and the axis of orthotropy of the material coincides with the corresponding natural base vector, the elastodynamic equations are reduced to a simple system for a displacement stress vector, using surface operators. In particular for spherical geometry, equations are obtained which generalize to orthotropic pre-stress those given by Alterman et al. and Takeuchi & Saito.