Thin elastic shells with variable thickness for lithospheric flexure of one-plate planets

Planetary topography can either be modelled as a load supported by the lithosphere, or as a dynamic effect due to lithospheric flexure caused by mantle convection. In both cases the response of the lithosphere to external forces can be calculated with the theory of thin elastic plates or shells. On one-plate planets the spherical geometry of the lithospheric shell plays an important role in the flexure mechanism. So far the equations governing the deformations and stresses of a spherical shell have only been derived under the assumption of a shell of constant thickness. However, local studies of gravity and topography data suggest large variations in the thickness of the lithosphere. In this paper, we obtain the scalar flexure equations governing the deformations of a thin spherical shell with variable thickness or variable Young's modulus. The resulting equations can be solved in succession, except for a system of two simultaneous equations, the solutions of which are the transverse deflection and an associated stress function. In order to include bottom loading generated by mantle convection, we extend the method of stress functions to include loads with a toroidal tangential component. We further show that toroidal tangential displacement always occurs if the shell thickness varies, even in the absence of toroidal loads. We finally prove that the degree-one harmonic components of the transverse deflection and of the toroidal tangential displacement are independent of the elastic properties of the shell and are associated with translational and rotational freedom. While being constrained by the static assumption, degree-one loads can deform the shell and generate stresses. The flexure equations for a shell of variable thickness are useful not only for the prediction of the gravity signal in local admittance studies, but also for the construction of stress maps in tectonic analysis.     

Global anisotropic phase velocity maps for higher mode Love and Rayleigh waves

It is well established that the Earth's uppermost mantle is anisotropic, but there are no clear observations of anisotropy in the deeper parts of the mantle. Surface waves are well suited to observe anisotropy since they carry information about both radial and azimuthal anisotropy. Fundamental mode surface waves, for commonly used periods up to 200 s, are sensitive to structure in the first few hundred kilometres, and therefore, do not provide information on anisotropy below. Higher mode surface waves have sensitivities that extend to and beyond the transition zone, and should thus give insight about azimuthal anisotropy at greater depths. We have measured higher mode Love and Rayleigh phase velocities using a model space search approach, which provides us with consistent relative uncertainties from measurement to measurement and from mode to mode. From these phase velocity measurements, we constructed global anisotropic phase velocity maps. Prior to inversion, we determine the optimum relative weighting for anisotropy. We present global azimuthal phase velocity maps for higher mode Rayleigh waves (up to the sixth higher mode) and Love waves (up to the fifth higher mode) with corresponding average model uncertainties. The anisotropy we derive is robust within the uncertainties for all modes. Given the ray theoretical sensitivity kernels of Rayleigh and Love wave modes, the source of anisotropy is complex, but mainly located in the asthenosphere and deeper. Our models show a good correspondence with other studies for the fundamental mode, but we have been able to achieve higher resolution.     

Identification of symmetry planes in weakly anisotropic elastic media

A procedure is proposed to obtain symmetry properties of weakly anisotropic (WA) elastic media by giving 15 WA parameters of qP wave in an arbitrary Cartesian coordinate system. The 15 WA parameters, which linearly depend on 21 elastic elements, form a complete set to determine the symmetry planes in WA materials. The procedure is based on the eigenvalue problems of two matrices. One of the matrices consists of the Voigt and dilatational matrices, and the other is an acoustic tensor defined by an irreducible, completely symmetric and traceless, fourth-rank tensor resulting from decomposition of the fourth-rank elastic tensor. If the eigenvectors are taken as the axes of a new coordinate system (called symmetry Cartesian coordinate system), the transformation of WA parameters from an arbitrary Cartesian coordinate system to a symmetry Cartesian coordinate system can reduce the number of distinct WA parameters of elastic materials except in triclinic medium.     

Wide-band coupling of Earth's normal modes due to anisotropic inner core structure

We investigate the importance of wide-band coupling of normal modes due to inner core anisotropy. We compare four different seismic models of inner core anisotropy, which were obtained by others using the splitting of Earth's normal modes. These models have been developed using a self-coupling (SC) approximation, which assumes that coupling between nearby modes through anisotropic inner core structure is negligible. We test the SC approximation by comparing the frequencies and quality factors of 90 inner core sensitive modes, computed for these models using either the SC approximation or full-coupling (FC) among large groups of modes. We find significant shifts in the quality factors and frequencies for some modes. Groups of modes which significantly couple together are constructed for six target modes. These groups are model dependent and in some cases contain large numbers of modes. Synthetic seismograms are calculated to show that the difference between SC and FC is observable on the scale of seismograms and of the same order of magnitude as the difference between synthetic and observed seismograms. Thus, future models of inner core anisotropy should take cross-coupling between large groups of modes into account.     

Propagation of harmonic plane waves in a general anisotropic porous solid

Wave propagation is studied in a general anisotropic poroelastic solid. The presence of dissipation due to fluid-viscosity as well as hydraulic anisotropy of pore permeability are also considered. Biot's theory is used to derive a system of modified Christoffel equations for the propagation of plane harmonic waves in porous media. A non-trivial solution of this system is ensured by a determinantal equation. This equation is separated into two different polynomial equations. One is the quartic equation whose roots represent the complex velocities of four attenuating waves in the medium. The other is a eighth-degree polynomial whose roots represent the vertical slowness values for the four waves propagating upward and downward in a finite porous medium. Procedure is explained to associate the numerically obtained roots with the waves propagating in the medium. The slowness surfaces of waves reflected at the boundary of the medium are computed for a realistic numerical model. The behaviours of phase velocity surfaces are analysed with the help of numerical examples.     

Existence of transverse waves in anisotropic poroelastic media

Out of the four waves in an anisotropic poroelastic medium, two are termed as quasi-transverse waves. The prefix 'quasi' refers to their polarizations being nearly, but not exactly, perpendicular to direction of propagation. In this composite medium, unlike perfectly elastic medium, the propagation of a longitudinal wave along a phase direction may not be accompanied by transverse waves. The existence of a transverse wave in anisotropic poroelastic media is ensured by the two equations restricting the choice of elastic coefficients of porous aggregate as well as fluid–solid coupling. Necessary and sufficient conditions for the existence of transverse waves along the coordinate axes and in the coordinate planes for general anisotropy are discussed. The discussion is extended to the case of orthotropic materials and existence for few specific phase directions is also explored. The conditions for the transverse waves decided on the basis of their apparent polarizations, that is, particle motion being perpendicular to ray direction, are also discussed. For a particular numerical model, the existence of these apparent transverse waves is solved numerically for phase directions in coordinate planes. For general directions of phase propagation, the existence of these transverse waves is checked graphically for the chosen numerical model.     

Gravitational torque on the inner core and decadal polar motion

A decadal polar motion with an amplitude of approximately 25 milliarcsecs (mas) is observed over the last century, a motion known as the Markowitz wobble. The origin of this motion remains unknown. In this paper, we investigate the possibility that a time-dependent axial misalignment between the density structures of the inner core and mantle can explain this signal. The longitudinal displacement of the inner core density structure leads to a change in the global moment of inertia of the Earth. In addition, as a result of the density misalignment, a gravitational equatorial torque leads to a tilt of the oblate geometric figure of the inner core, causing a further change in the global moment of inertia. To conserve angular momentum, an adjustment of the rotation vector must occur, leading to a polar motion. We develop theoretical expressions for the change in the moment of inertia and the gravitational torque in terms of the angle of longitudinal misalignment and the density structure of the mantle. A model to compute the polar motion in response to time-dependent axial inner core rotations is also presented. We show that the polar motion produced by this mechanism can be polarized about a longitudinal axis and is expected to have decadal periodicities, two general characteristics of the Markowitz wobble. The amplitude of the polar motion depends primarily on the Y ¹₂ spherical harmonic component of mantle density, on the longitudinal misalignment between the inner core and mantle, and on the bulk viscosity of the inner core. We establish constraints on the first two of these quantities from considerations of the axial component of this gravitational torque and from observed changes in length of day. These constraints suggest that the maximum polar motion from this mechanism is smaller than 1 mas, and too small to explain the Markowitz wobble.