Hydrostatic figure and related properties of the Earth

Summary. The second-order internal theory of hydrostatic equilibrium initially developed by Kopal is extended to form equations equivalent to those of de Sitter. This new development is used together with the density distribution of Dziewonski, Hales & Lapwood to compute flattening, precessional constant and moments of inertia of the hydrostatic earth. These values are compared with those of the real Earth.     

On the density distribution within the Earth

The distribution of density as a function of position within the Earth is much less well constrained than the seismic velocities. The primary information comes from the mass and moment of inertia of the Earth and this information alone requires that there be a concentration of mass towards the centre of the globe. Additional information is to be found in the frequencies of the graver normal modes of the Earth which are sensitive to density through self-gravitation effects induced in deformation.
  The present generation of density models has been constructed using linearized inversion techniques from earlier models, which ultimately relate back to models developed by Bullen and based in large part on physical arguments. A number of experiments in non-linear inversion have been conducted using the PREM reference model, with fixed velocity and attenuation, but with the density model constrained to lie within fixed bounds on both density and density gradient. A set of models is constructed from a uniform probability density within the bound and slope constraints. Each of the resultant density models is tested against the mass and moment of inertia of the Earth, and for successful models a comparison is made with observed normal mode frequencies. From the misfit properties of the ensemble of models the robustness of the density profile in different portions of the Earth can be assessed, which can help with the design of parametrization for future reference models. In both the lower mantle and the outer core it would be desirable to allow a more flexible representation than the single cubic polynomial employed in PREM.     

The moment of inertia of Mars and the existence of a core

Summary. A relation is obtained between the true value of the moment of inertia of a planet and the value calculated from the dynamical form factor, J ₂, on the assumption of hydrostatic equilibrium. The result is applied to Mars and it is shown that the difference between the true and calculated moments of inertia is probably insignificant in considering models of the interior of Mars and in particular does not affect an argument for a core based on models calculated by Lyttleton.     

Implications of a non-adiabatic density gradient for the Earth's viscoelastic response to surface loading

The response of a viscoelastic Earth to the melting of the Late Pleistocene ice sheets has been the subject of a number of investigations employing PREM. In PREM, a non-adiabatic density gradient (NADG) exists in the upper mantle, and to understand the implications of this model it is thus important to examine the effects of this NADG on the Earth's response to surface loads. This paper is based on the assumption that the contribution to the depth dependence of the density that is not due to self-compression is due to compositional change. This contribution is referred to as 'non-adiabatic'. We evaluate the effects of a non-adiabatic density jump (NADJ) for the 670  km discontinuity and the NADG in the upper mantle by adopting a compressible earth model with both a compositional density gradient and a density jump. Numerical calculations based on these models indicate that the magnitude of the Earth's response associated with the NADG is much smaller than that associated with the NADJ at 670  km depth. It is also confirmed that the higher modes associated with the NADJ and the NADG are much more sensitive to the existence of an elastic lithosphere than the fundamental modes associated with the density jumps at the surface and core–mantle boundary.     

On crustal corrections in surface wave tomography

Mantle models from surface waves rely on good crustal corrections. We investigated how far ray theoretical and finite frequency approximations can predict crustal corrections for fundamental mode surface waves. Using a spectral element method, we calculated synthetic seismograms in transversely isotropic PREM and in the 3-D crustal model Crust2.0 on top of PREM, and measured the corresponding time-shifts as a function of period. We then applied phase corrections to the PREM seismograms using ray theory and finite frequency theory with exact local phase velocity perturbations from Crust2.0 and looked at the residual time-shifts. After crustal corrections, residuals fall within the uncertainty of measured phase velocities for periods longer than 60 and 80 s for Rayleigh and Love waves, respectively. Rayleigh and Love waves are affected in a highly non-linear way by the crustal type. Oceanic crust affects Love waves stronger, while Rayleigh waves change most in continental crust. As a consequence, we find that the imperfect crustal corrections could have a large impact on our inferences of radial anisotropy. If we want to map anisotropy correctly, we should invert simultaneously for mantle and crust. The latter can only be achieved by using perturbation theory from a good 3-D starting model, or implementing full non-linearity from a 1-D starting model.     

Coseismic rotation changes from the 2004 Sumatra earthquake: the effects of Earth's compressibility versus earthquake induced topography

Polar motion is modelled for the large 2004 Sumatra earthquake via dislocation theory for an incompressible elastic earth model, where inertia perturbations are due to earthquake-triggered topography of density–contrast interfaces, and for a compressible model, where inertia perturbation due to compression-dilatation of Earth's material is included; density and elastic parameters are based on a multilayered reference Earth. Both models are based on analytical Green's functions, propagated from the centre to the Earth's surface. Preliminary and updated seismological solutions are considered in elucidating the effects of improving earthquake parameters on polar motion. The large Sumatra thrust earthquake was particularly efficient in driving polar motion since it was responsible for large material displacements occurring orthogonally to the strike of the earthquake and to the Earth's surface, as imaged by GRACE gravity anomalies over the earthquake area. The effects of earthquake-induced topography are four times larger than the effects of Earth's compressibility, for l = 2 geopotential components. For varying compressional Earth properties and seismic solution, modelled polar motion ranges from 8.6 to 9.4 cm in amplitude and between 117° and 130° east longitude in direction. The close relationship between polar motion direction, earthquake longitude and thrust nature of the event, are established in terms of basic physical concepts.     

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.     

On the amplitude of the Chandler wobble

The change in the inertia tensor of the Earth, due to the mass shift following a seismic event, has been computed by several authors for non-rotating earth models. Rotation is taken into account in the present paper, and the additional change in the inertia tensor is computed for an equivalent earth model, in which the axis of geometrical symmetry becomes tilted instead of the axis of greatest inertia. Rotation is thus seen to produce an increase by a factor 1.4 in the amplitude variation of the Chandler wobble, with respect to the non-rotating case, which, when added to the 1.4 amplitude increase due to the precessional re-adjustment of the equatorial bulge, gives a factor of 2 increase of the Chandler wobble amplitude with respect to the case of a rigid earth model.     

The effects of post-seismic motions on the moment of inertia of a stratified viscoelastic earth with an asthenosphere

Summary. We give the analytical formulation for calculating the transient displacement of fields produced by earthquakes in a stratified, selfgravitating, incompressible, viscoelastic earth. We have evaluated the potential of viscous creep in the asthenosphere in exciting the Chandler wobble by a four-layer model consisting of an elastic lithosphere, a two-layer Maxwell viscoelastic mantle, and an inviscid core. The seismic source is modelled as an inhomogeneous boundary condition, which involves a jump condition of the displacement fields across the fault in the lithosphere. The response fields are derived from the solution of a two-point boundary value problem, using analytical propagator matrices in the Laplace-transformed domain. Transient flows produced by post-seismic rebound are found to be confined within the asthenosphere for local viscosity values less than 10²⁰P. The viscosity of the mantle below the low-viscosity channel is kept at 10²²P. For low-viscosity zones with widths greater than about 100 km and asthenospheric viscosities less than 10¹⁸P, we find that viscoelasticity can amplify the perturbations in the moment of inertia by a factor of 4–5 above the elastic contribution within the time span of the wobble period. We have carried out a comparative study on the changes of the inertia tensor from forcings due to surface loading and to faulting. In general the global responses from faulting are found to be much more sensitive to the viscosity structure of the asthenosphere than those produced from surface loading.     

Numerical simulation of salinity anomaly propagation in the Nordic seas and the Arctic Ocean

We investigate the response of the Nordic seas-Arctic Ocean system to surface freshwater flux anomalies that we regard as typical for long-term atmospheric variability. We employ response experiments with a coupled sea ice-ocean model where we introduce a surface freshwater flux anomaly (A) over the Norwegian Sea and (B) in the Laptev Sea. Case A offers an explanation for the intermediate depth salinity changes observed in the Amundsen Basin. The signal observed there belongs to an original perturbation that, according to the model, occurred around a decade earlier. Salinity fluctuations in the Laptev Sea could play a role in changes in the near surface salinity in the Amundsen Basin.     

Regional tomographic inversion of the amplitude and phase of Rayleigh waves with 2-D sensitivity kernels

In this study, we test the adequacy of 2-D sensitivity kernels for fundamental-mode Rayleigh waves based on the single-scattering (Born) approximation to account for the effects of heterogeneous structure on the wavefield in a regional surface wave study. The calculated phase and amplitude data using the 2-D sensitivity kernels are compared to phase and amplitude data obtained from seismic waveforms synthesized by the pseudo-spectral method for plane Rayleigh waves propagating through heterogeneous structure. We find that the kernels can accurately predict the perturbation of the wavefield even when the size of anomaly is larger than one wavelength. The only exception is a systematic bias in the amplitude within the anomaly itself due to a site response.
An inversion method of surface wave tomography based on the sensitivity kernels is developed and applied to synthesized data obtained from a numerical simulation modelling Rayleigh wave propagation over checkerboard structure. By comparing recovered images to input structure, we illustrate that the method can almost completely recover anomalies within an array of stations when the size of the anomalies is larger than or close to one wavelength of the surface waves. Surface wave amplitude contains important information about Earth structure and should be inverted together with phase data in surface wave tomography.     

Effect of the viscoelastic lithosphere on polar wander speed caused by the Late Pleistocene glacial cycles

Previous studies of the wander of the rotation pole associated with the Late Pleistocene glacial cycles indicate that the predicted polar wander speed is sensitive to the density jump at the 670 km discontinuity, the thickness of the elastic lithosphere, and the lower mantle viscosity. In particular, the M1 mode related to the density jump at 670 km depth has been shown to contribute a dominant portion of predicted polar wander speed for sufficiently small lower mantle viscosities. In this study, we examine the sensitivity of polar wander to variations in the viscosity of the viscoelastic lithosphere using simplified compressible Maxwell viscoelastic earth models. Model calculations for earth models with a viscoelastic lithosphere of finite viscosity indicate that the contribution of the M1 mode is similar to those associated with the density discontinuity at the core–mantle boundary (C0 mode) and the lithosphere (L0 mode). We speculate that this is due to the interaction between the M1 mode and the transient mode associated with the viscoelastic lithosphere, which reduces the magnitude of polar wander rates. Therefore, the M1 mode does not contribute a dominant portion of the predicted polar wander speed for earth models with a viscoelastic lithosphere of finite viscosity. In this case, predictions of polar wander speed as a function of lower mantle viscosity exhibit the qualitative form of an 'inverted parabola', as predicted for the J ˙₂ curve. We caution, however, that these results are obtained for simplified earth models, and the results for seismological earth models such as PREM may be complicated by the interaction between the M1 mode and the large set of transient modes.     

Observation and modelling of fault-zone fracture seismic anisotropy – II. P-wave polarization anomalies

Summary. In Part I of this paper we modelled shear-wave splitting observed in crystalline rock bordering an active, normal fault-zone at Oroville, California, with Červený's ray-tracing system applied to anisotropic heterogeneous media using Hudson's formulation of elastic constants for a medium containing aligned cracks. In Part II we use the ray-tracing results of Part I to quantitatively interpret P -wave polarization anomalies observed in the three-component seismograms recorded in the Oroville fault zone. We show that the eigenvectors of the first-order Christoffel tensor defined by the ray-tracing slowness vector and Hudson's first-order anisotropic corrections to the isotropic elastic tensor correctly account for P -wave first motion that deviates from the ray vector.     

Effects of a mushy transition zone at the inner core boundary on the Slichter modes

This article investigates the effects of a mushy inner core boundary on the eigenperiods of the Slichter modes for a simple, but realistic, earth model (rotating, spherical configuration, elastic inner core and mantle, neutrally stratified, inviscid, compressible liquid core). It is found that the influence of the mushy boundary layer is substantial compared with some other effects, such as those from elasticity of the mantle, non-neutral stratification of the liquid outer core and ellipticity of the Earth and centrifugal potential. The results obtained here may set a lower bound on the eigenperiods of the Slichter modes for a realistic earth model. For example, for a PREM model, the lower bound of the central period of the Slichter modes should be about 5.3 hr.     

Iceocean mass balance during the Late Pleistocene glacial cycles in view of CHAMP and GRACE satellite missions

During the last glacial cycles, global sea level dropped several times by about 120 m and large ice sheets covered North America, northern Europe and Antarctica during the glacial stages. The changes in the iceocean mass balance have displaced mantle material mainly via viscous flow, and the perturbation of the equilibrium figure of the Earth by glacial isostatic adjustment is still observable today in timedependent changes of gravitational and rotational observations. Contemporary iceocean mass balance from volume changes of polar ice caps also contributes to secular variations of the Earth's gravitational field.
In the near future, several satellite gravity missions will significantly improve the accuracy of the observed timedependent gravitational field. In view of the expected improvements in the observations, we predict glacially induced perturbations of the gravitational field, induced by Late Pleistocene and contemporary ice volume changes, for a variety of radial mantle viscosity profiles. We assess the degree of uncertainty for the glacially induced contributions to gravitational and rotational parameters, both in the spectral and the spatial domain.
Predictions of power spectra for the glacially induced freeair gravity and geoid anomalies are about one order of magnitude lower than the observed values, and uncertainties arising from different plausible viscosity profiles are around 0.150.4 mGal and 0.21.5 m, respectively. Uncertainties from different ice models are of secondary importance for the predicted power spectra. Predicted secular changes in geoid anomalies in formerly glaciated areas are mainly controlled by the viscosity profile and contemporary ice volume changes. We also show that the simple threelayer viscosity profiles currently employed for the majority of postglacial rebound studies represent a limited subset for model predictions of the timedependent gravitational field.     

Full three-dimensional tomography: a comparison between the scattering-integral and adjoint-wavefield methods

This paper analyses the computational issues of full 3-D tomography, in which the starting model as well as the model perturbation is 3-D and the sensitivity (Fréchet) kernels are calculated using the full physics of 3-D wave propagation. We compare two formulations of the structural inverse problem: the adjoint-wavefield (AW) method, which back-propagates the data from the receivers to image structure, and the scattering-integral (SI) method, which sets up the inverse problem by calculating and storing the Fréchet kernels for each data functional. The two inverse methods are closely related, but which one is more efficient depends on the overall problem geometry, particularly on the ratio of sources to receivers, as well as trade-offs in computational resources, such as the relative costs of compute cycles to data storage. We find that the SI method is computationally more efficient than the AW method in regional waveform tomography using large sets of natural sources, although it requires more storage.     

A search for Chandler and nearly diurnal free wobbles using Liouville equations

Summary. The basic equations describing the dynamical effects of the Earth's fluid core (Liouville, Navier-Stokes and elasticity equations) are derived for an ellipsoidal earth model without axial symmetry but with an homogeneous and deformable fluid core and elastic mantle.
We develop the balance of moment of momentum up to the second order and use Love numbers to describe the inertia tensor's variations. The inertial torque takes into account the ellipticity and the volume change of the liquid core. On the core—mantle boundary we locate dissipative, magnetic and viscous torques. In this way we obtain quite a complete formulation for the Liouville equations.
These equations are restricted in order to obtain the usual Chandler and nearly diurnal eigenfrequencies.
Then we propose a method for calculating the perturbations of these eigenfrequencies when considering additional terms in the Liouville equations.     

Random walk of the Earth's pole related to the Chandler wobble excitation

Summary. The equation governing the polar motion shows that the polar secular drift and the Chandler wobble amplitude are related to each other. In particular, a drift of the mean pole position comes out as a consequence of the maintenance of the Chandler wobble by possible step perturbations of the Earth's inertia tensor.
The minimum excitation functions necessary to explain the Chandler wobble amplitude variations for the period 1901–84 are derived from the Chandler term, with the hypothesis that the excitations follow a uniform random distribution in time. It is shown that they have the statistical properties of the steps of a two-dimensional random walk. These functions are then used to derive, from a statistical simulation, a lower limit of the secular drift which may result from the excitation of the Chandler wobble.
The drift generated by the random walk is of the same order of magnitude as the observed secular drift for the period 1901–84, but their time dependence is different. This indicates that the observed secular drift cannot be explained as the consequence of an excitation of the Chandler wobble by random steps of the Earth's inertia tensor. However, the possible contribution of the Chandler wobble excitation to the polar drift has to be taken into account when other mechanisms, such as lithospheric rebound related to deglaciation, are proposed.     

Whole mantle SH velocity model constrained by waveform inversion based on three-dimensional Born kernels

A whole mantle SH velocity model is obtained by using a unique data set and techniques. Body and surface waveforms including major and multi-orbit phases are used as a data set and are inverted by using 3-D Born kernels. The resultant model, SH18CE, reveals the different natures of the two major upwelling systems: the strong low velocity anomalies beneath Africa extend for more than 1000 km from the core–mantle boundary (CMB), whereas those beneath the Pacific are restricted to 300–400 km from the CMB. The results also show the variable natures of stagnant slabs on the 670 discontinuity around Japan: the depths of the strongest high velocity anomalies within the stagnant slabs are different region by region, which is consistent with the detailed delay time tomography model in this area.     

Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.