The forced nutations of an elliptical, rotating, elastic and oceanless earth

Summary. We compute the luni-solar forced nutations of an elliptical, rotating, self-gravitating, elastic, hydrostatically prestressed and oceanless earth. Several recent structural models are considered, each possessing a fluid outer core and solid inner core. Complete results are given for the nutation of the 'axis of figure for the Tisserand mean surface' which best represents the observational effects of the Earth's nutational motion. Differences between results for different structural models are observationally insignificant. Differences between our results and Molodensky's are as large as ∼ 0.002 arcsec at six month and at 18.6 yr.     

A diurnal resonance in the ocean tide and in the Earth's load response due to the resonant free 'core nutation'

Summary. The luni-solar forced nutations and body tide are believed to be resonant at frequencies near (1 + 1/460) cycle sidereal day⁻¹ as seen from the rotating Earth. This resonance is due to the Earth's rotating, elliptical fluid core. We show here that tides in the open ocean and the Earth's response to those tides must also be resonant at (1 + 1/460) cycle day⁻¹. We examine these resonant oceanic effects on the Earth's nutational motion and on the body tide. Effects on the forced nutations might be as large as 0.002 arcsec at 18.6 yr. The effects on the observed resonance in the body tide are more important. For tidal gravity, for example, the difference between K ₁ and 0 ¹ which is usually used to determine the resonance, can be perturbed by 30 per cent or more due to the oceanic resonance effects.     

A survey of Earth strain tides in Great Britain

Summary . We report 35 measurements of Earth strain tides at 16 near-surface sites in Great Britain. This is the first widespread survey undertaken specifically to examine the problem of the inhomogeneity of elastic strain fields near the Earth's surface. Some sites were instrumented intensively in order to examine variations of tidal admittance over distances of several hundred metres, while measurements from single instruments at other sites were compared with theoretically predicted strain tides. After allowing for cavity and topographic effects, our data show variations of up to 50 per cent in tidal admittance. We interpret such large anomalies as being due to variations in the regional elastic parameters of wavelength about a 100 m or less. The data indicate that strain measurements from single instrument sites must be interpreted with caution.     

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.     

Could the energy near the FCN and the FICN be explained by luni-solar or atmospheric forcing?

The observed time-series of precession/nutation show residuals with respect to an empirical model based on the rigid Earth theoretical nutations and a frequency dependent transfer function with resonances to the Earth's normal modes. These residuals display energy mainly in the frequency domain around 430 and 500 days in the inertial frame. In this frequency band, the energy is possibly related to two normalmode frequencies: the free core nutation (FCN) and the free inner core nutation (FICN). In this paper, we examine the possibility of obtaining this energy from the resonance effect induced by a luni-solar (or planetary) forcing, or by an atmospheric forcing at a frequency very close to these Earth free nutations. The amplification factor due to the resonance is computed from an analytical formula expressed in the case of a simplified three-layer ellipsoidal rotating earth (with an elastic inner core, a liquid outer core and an elastic mantle), as well as the empirical formula based on the analysis of VLBI observations. For the tidal forcing, the theoretical results do not show any resonance at the level of precision we have examined but it is still possible to find one frequency near the FCN or FICN frequencies which could be excited. In contrast, for the atmospheric pressure the level of energy needed could be obtained from the diurnal pressure, depending on the noise level of the Earth's global pressure. We also show that the combination of three waves can explain the observed decrease of energy with time. While the tidal potential amplitudes are too small, a pressure noise level of 0.5 Pa would be sufficient to excite these waves.     

Topographic, geologic and cavity effects on the harmonic content of tidal strain

Summary. The combined effects of topography, geology and cavities upon the measurements of tidal strains are considered. A theoretical model is proposed in which the topography and geology are modelled with finite elements, and the cavity (assumed to be elliptical) modelled using the analytical theory. The three effects are then superimposed to form a total model. Previous experimental results obtained for the solid tidal strain in Eastern Australia (Blair & Sydenham) are compared with the present theoretical model and reasonable agreement is found assuming values of the Love numbers for the Gutenberg—Bullen A earth model in which h = 0.61 and l = 0.08 (Farrell). Further experimental results are presented for two strain meters, only one of which is not significantly influenced by the cavity. The relative difference, in harmonic content between the recordings of both strain meters is compared with that predicted by the theoretical model.     

Effect of the fluid core on changes in the length of day due to long period tides

Summary. The long period luni-solar tidal potential is known to cause periodic changes in the Earth's rotation rate. We find that the effect of a dissipationless fluid outer core is to reduce the amplitudes of these tidal perturbations by ∼ 11 per cent. When the fluid core effect is added to Agnew & Farrell's estimate of the effect of an equilibrium ocean, the result is in accord with observation.
The effects of dissipative processes within the fluid core are also examined. We find out-of-phase perturbations which could be as large as ∼ 10ms at 18.6 yr. We conclude, however, that the poorly understood decade fluctuations in the Earth's rotation rate will prohibit observation of this effect.     

Quasi-static analysis of strike fault growth in layered media

We study the effects of structural inhomogeneity on the quasi-static growth of strike-slip faults. A layered medium is considered, made up of an upper layer bounded by a free surface and welded to a lower half-space with different elastic property. Mode III crack is employed as a mathematical model of strike-slip fault, which is nucleated in the lower half-space and then propagates towards the interface. We adopt FEM-β, newly proposed analysis method for failure, to simulate the quasi-statistic crack growth governed by the stress distribution in layered media. Our results show that along planar traces across interfaces a compliant upper layer has significant effects on promoting/suppressing crack growth before/after its extension into the layer and vice versa for a rigid one. This proposes a possibility that surface breaks due to strike-slip faulting could be arrested by deposit layers at the topmost part of the Earth's crust.     

Tidal gravity and ocean tide loading in Europe

b
The results are presented from tidal gravity measurements at five sites in Europe using LaCoste and Romberg ET gravimeters. Improvements that we have made to the accuracies of these gravimeters are discussed. It is shown that the 'standard' calibration of the International Center for Earth Tides, used for worldwide tidal gravity profiles, is 1.2 per cent too high. The M₂ and O₁ observations are compared with model calculations of the Earth's body tide and ocean tide loading and it is shown that there is a very significant improvement in the agreement between observations and models compared to that obtained with previous tidal gravity measurements. For O₁, where the ocean tide loading and attraction in central Europe is only 0.4 per cent of the body tide, our measurements verify that the Dehant-Wahr anelastic body tide model gravimetric factor is accurate to 0.2 per cent. It is also shown that the effects of lateral heterogeneities in Earth structure on tidal gravity are too small to explain the large anomalies in previously published tidal gravity amplitudes. The observations clearly show the importance of conserving tidal mass in the Schwiderski ocean tide model. For sites in central Europe, the M₂ and O₁ observations and the models are in agreement at the 0.1 μgal (10⁻⁹ m s⁻²) level and tidal corrections to this accuracy can now be made to absolute gravity measurements.     

Tides and the evolution of the Earth—Moon system

Summary. A model of the tides in a hemispherical ocean is used to investigate the effect of changes in the Earth's rotation rate on the power dissipated by the ocean tides. The results obtained are then used in an idealized astronomical model to investigate how they affect the history of the Earth—Moon system.
Using the tidal model it is found that at rotation rates higher than that of the present Earth, the power dissipated by the semi-diurnal tides in the ocean drops off rapidly as a result of the increased tidal frequency. Thus if the Earth's rotation rate is doubled from its present value, then the rate of energy dissipation in the ocean is reduced to approximately one-third of its present value and the tidal torque is reduced by a factor of about 6.
The present value for secular acceleration of the Moon, calculated from the results of the tidal model is -30.5 arcsec century^-2. Using this value in the astronomical model, which has the Moon and Sun in circular orbits above the equator, and assuming that the tidal torque is independent of the tidal frequency, the Gerstenkorn event is predicted to have occurred 1.3 × 10⁹ yr ago.
When the astronomical model is run with a torque determined at all times from the tidal model, the reduction in the energy dissipated early in the history of the system, leads to a Gerstenkorn date of 5.3 × 10⁹ yr ago. However, dissipation within the solid earth is found to be important early in the history of the system and when this effect is included it gives a date for the Gerstenkorn event of 3.9 × 10⁹ yr ago.     

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.     

Induced magnetic field of the equatorial electrojet at the surface of the Earth

Summary. The northward component of the induced magnetic field due to the equatorial electrojet at the Earth's surface is calculated using a more realistic local time variaton of the external field due to the electrojet than is provided for by models of the electrojet currently used in induction calculations. It is seen that appreciable induction effects can be expected about an hour before local noon for the kind of local time variation considered. Our results are in qualitative agreement with direct observations of Earth currents in the equatorial region in Nigeria. At local times when observable induction effects are present, the magnetic field due to the electrojet is necessarily three-dimensional; hence in order to obtain the internal part directly from the observed total field due to the electrojet at the Earth's surface, a three-dimensional formulation is required.     

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.     

Tidal tilt at Llanrwst, North Wales: tidal loading and Earth structure

Summary. Tidal tilt measurements in a disused lead mine at Llanrwst, North Wales, are used to examine the magnitudes of the local strain—tilt coupling perturbations that can be expected in mine and tunnel installations. From the spatial variability of the tilt signal and comparisons with tidal loading calculations on radially stratified Earth models, it is concluded that typical strain—tilt coupling perturbations are at the ± 10 per cent level. The general problem of using tidal loading tilt and strain measurements to investigate the elastic parameters of layered models of the crust and upper mantle is studied by examining the sensitivities of the load Green's functions to these parameters. It is shown that the load tilt essentially gives information on Young's modulus whereas the load strain, although more sensitive to the structure, depends on both Young's modulus and Poisson's ratio. These investigations show that, if load tilt or load strain can be measured to an accuracy of 1 or 2 per cent, then they can be used to improve models of the crust and upper mantle. It is concluded that due to the local strain—tilt coupling perturbations this accuracy cannot be achieved in mines or tunnels and that future measurements should be made in boreholes or near surface trenches. The tidal loading models are also used to examine the spatial variation of the M² loading tilt in Europe.     

Surface wave tomography for azimuthal anisotropy in a strongly reduced parameter space

Large scale seismic anisotropy in the Earth's mantle is likely dynamically supported by the mantle's deformation; therefore, tomographic imaging of 3-D anisotropic mantle seismic velocity structure is an important tool to understand the dynamics of the mantle. While many previous studies have focused on special cases of symmetry of the elastic properties, it would be desirable for evaluation of dynamic models to allow more general axis orientation. In this study, we derive 3-D finite-frequency surface wave sensitivity kernels based on the Born approximation using a general expression for a hexagonal medium with an arbitrarily oriented symmetry axis. This results in kernels for two isotropic elastic coefficients, three coefficients that define the strength of anisotropy, and two angles that define the symmetry axis. The particular parametrization is chosen to allow for a physically meaningful method for reducing the number of parameters considered in an inversion, while allowing for straightforward integration with existing approaches for modelling body wave splitting intensity measurements. Example kernels calculated with this method reveal physical interpretations of how surface waveforms are affected by 3-D velocity perturbations, while also demonstrating the non-linearity of the problem as a function of symmetry axis orientation. The expressions are numerically validated using the spectral element method. While challenges remain in determining the best inversion scheme to appropriately handle the non-linearity, the approach derived here has great promise in allowing large scale models with resolution of both the strength and orientation of anisotropy.     

Strain-tide spectroscopy

Summary. We have used two years of strain-tide data to study the response of the Earth to the diurnal and semidiurnal tidal excitations. Our results show that there is significant structure in the response of the Earth to tidal excitations near one cycle/sidereal day. This structure agrees with the resonance behaviour predicted from calculations of the forced elastic-gravitational response of an elliptical, rotating earth with a liquid outer core. The data can also be used to test for possible preferred frames and spatial anisotropies. We find that upper bounds on the parameterized post-Newtonian (PPN) parameters which characterize these effects are α₂≤ 0.007 and ζ_w≤ 0.005. We also infer that there is a significant frequency dependence to the ocean load near one cycle/day and that the coherence between strain and barometric pressure fluctuations are significant at periods longer than a few days.     

An excitation mechanism for the free 'core nutation'

Summary. The Earth is believed to possess a free nutational mode due to its rotating, elliptical, fluid core, with an eigenfrequency of approximately (1 + 1/460) cycle sidereal day⁻¹ as seen from the sidereally rotating Earth. This free 'core nutation' has not yet been undisputably observed. Furthermore, there has been considerable doubt that any known mechanism could excite this mode to an observable level. We show here that diurnal atmospheric and oceanic loading of the Earth's surface provide an efficient excitation mechanism which depends critically on the physical damping of the mode. Possible effects of the mode on geodetic measurements are discussed. We also consider the effects of 'wobble' and 'nutation' on astrometric observations.     

The effect of polar wander on the tides of a hemispherical ocean

Summary. A numerical model is constructed of the tides in a hemispherical ocean driven by the forces corresponding to the Y₂^–2 equilibrium tide. The model is used to study how tidal dissipation is affected by changes in the position of the ocean relative to the Earth's rotational axis and to test a hypothesis concerning the Gerstenkorn event.
As the position of the Earth's axis is varied with respect to the ocean, the model shows changes in the dissipation rate due to the changing position and importance of individual resonances of the ocean. However, a cooperative effect is also observed which results, for an ocean of depth 4400 m, in broad frequency bands near 10 rad day⁻¹ and-6 rad day⁻¹ in which the dissipation rate remains high.
The cooperative effect is found to arise from the existence, in an unbounded ocean, of resonances at these frequencies which match the tidal forces. When ocean boundaries are introduced, the new resonances near these frequencies contain a large component of the underlying resonance and as a result are themselves a good match to the driving forces.
For the real ocean, these findings imply that changes in the position of the pole, and also possibly changes in the shape of the ocean, will on average have little effect on the energy dissipated by the tides. However in the past changes in the mean depth and area of the ocean or the increased rotation rate of the Earth may have resulted in a smaller dissipation rate.     

On the Stability of a Spherical Gravitating Compressible Liquid Planet without Spin

During the past ten years or so there has been considerable discussion in the literature regarding the author's 1963 contention that (neglecting temperature effects and spin) the Earth's liquid core cannot be stable unless the Adams-Williamson condition relating density distribution and compressibility holds there.
The present paper throws light on this question by showing mathematically that a sphere of gravitating compressible liquid cannot be internally stable unless this condition is fulfilled. Physical reasons for the necessity of this condition, which implies that particles of the liquid are in neutral equilibrium, are also discussed. By internal stability is meant stability of the density distribution while the spherical shape is maintained.
The question of shape stability is not treated here, since it may be assumed that the Earth's mantle is sufficiently rigid to keep the core essentially spherical.
The liquid is assumed to be a perfect fluid, elastic, and in the discussion only small strains are considered from an equilibrium configuration of initial hydrostatic stress. Furthermore thermodynamic effects are neglected and there is no spin.     

Surface-wave polarization data and global anisotropic structure

In the past few years, seismic tomography has begun to provide detailed images of seismic velocity in the Earth's interior which, for the first time, give direct observational constraints on the mechanisms of heat and mass transfer. The study of surface waves has led to quite detailed maps of upper-mantle structure, and the current global models agree reasonably well down to wavelengths of approximately 2000 km. Usually, the models contain only elastic isotropic structure, which provides an excellent fit to the data in most cases. For example, the variance reduction for minor and major arc phase data in the frequency range 7–15 mHz is typically 65–92 per cent and the data are fit to within 1–2 standard deviations. The fit to great-circle phase data, which are not subject to bias from unknown source or instrument effects, is even better. However, there is clear evidence for seismic anisotropy in various places on the globe. This study demonstrates how much (or little) the fit to the data is improved by including anisotropy in the modelling process. It also illuminates some of the trade-offs between isotropic and anisotropic structure and gives an estimate of how much bias is introduced by neglecting anisotropy. Finally, we show that the addition of polarization data has the potential for improving recovery of anisotropic structure by diminishing the trade-offs between isotropic and anisotropic effects.