Body tides on an elliptical, rotating, elastic and oceanless earth

Summary. The Earth's deformation caused by the luni-solar tidal force is defined as the 'body tide'. We compute the effects of the Earth's rotation and elliptical stratification on the body tide for a number of modern elastic structural models. Rotation and ellipticity within the mantle are found to affect tidal observations by about 1 per cent. A consequence is an improved estimate for the fluid core resonance in the diurnal tidal band. Agreement between results for the different structural models is very good. As a result, the results computed here can be used to model the tidal effects of a globally averaged, oceanless, rotating, elliptical and elastic earth to an accuracy of at least one part in 300.     

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.     

On the upper limit of undertone periods

Summary. The severity of the effect of the Earth's rotation on the transverse motions within the fluid core of the Earth is examined. It is shown that, for any model of the fluid core, the formalism adapted from free oscillation theory is applicable only for periods less than 12 hr. In the limiting case, as the frequency is decreased to 2 cycle day^-1, all modes of response in the fluid are excited to the same order.     

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.     

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.     

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.     

Asymptotic behaviour of decay rates of internal waves in a rotating stratified spherical shell

Summary. Boundary layer techniques are used to examine the dissipative decay of an internal oscillation that is a member of the inviscid spectrum of normal modes for a rotating fluid shell stratified under a radially directed gravitational field. A formula is derived for the decay factor on the so-called homogeneous spin-down time-scale. Estimates are obtained for the size of the decay factor as a function of wavelength, a function of the frequency and a function of a parameter A which measures the ratio of the stratification strength to the rotation strength. It is shown that all modes decay on the spin-down time-scale. The results are interpreted in the context of a model for the Earth's fluid core. It is observed that the presence of regions of unstable stratification may increase the decay rate for oscillations at frequencies less than twice the rotation frequency.     

Internal waves in a rotating stratified spherical shell: asymptotic solutions

Summary. Small amplitude oscillations of a rotating, density-stratified fluid bounded by a spherical shell are examined. No restrictions are placed on the thickness of the shell. The internal mode spectrum is examined in the complete rotation-stratification parameter range including the regime that is appropriate for a plausible stratification distribution in the Earth's fluid core. A mathematical model is derived in terms of an eigenvalue PDE of mixed type. The existence of oscillatory solutions is exhibited in the limits of no rotation and no stratification. The frequency spectrum is extended asymptotically away from these limiting cases. A reduction in the complexity of the PDE for modes oscillating at the inertial frequency is exploited. A variational formulation is constructed in which the stratification parameter is treated as an eigenvalue of the system for fixed wave frequency. The spectral information is again extended asymptotically away from these 'accessible' points. Although the PDE reduces to Laplace's tidal equations (LTE) only under stringent parameter restrictions, it is observed that aspects of the behaviour of low frequency LTE modes are reproduced in the general model.     

The period and Q of the Chandler wobble

Summary. We have extended our calculation of the theoretical period of the Chandler wobble to account for the non-hydrostatic portion of the Earth's equatorial bulge and the effect of the fluid core upon the lengthening of the period due to the pole tide. We find the theoretical period of a realistic perfectly elastic Earth with an equilibrium pole tide to be 426.7 sidereal days, which is 8.5 day shorter than the observed period of 435.2 day. Using Rayleigh's principle for a rotating Earth, we exploit this discrepancy together with the observed Chandler Q to place constraints on the frequency dependence of mantle anelasticity. If Qμ in the mantle varies with frequency σ as σα between 30 s and 14 months and if Qμ in the lower mantle is of order 225 at 30 s, we find that 0.04 ρ^α≤ 0.11; if instead Qμ in the lower mantle is of order 350 near 200 s, we find that 0.11 ≤α≤ 0.19. In all cases these limits arise from exceeding the 68 per cent confidence limits of ± 2.6 day in the observed period. Since slight departures from an equilibrium pole tide affect the Q much more strongly than the period we believe these limits to be robust.     

A model for energy dissipation at the mantle—core boundary

Summary. An existing experimentally verified model for energy dissipation in a processing spherical cavity filled with liquid assumed to be in a semirigidized state except for a viscous Ekman boundary layer is applied to the Earth's liquid core to assess energy dissipation magnitudes. Application of the model to the best available Earth data occurs at the derived energy dissipation maximum for the model. Other existing research showing that the Earth's atmosphere appears to adjust to a state of maximum dissipation led to generic models for systems of maximum dissipation. The maximum dissipation mantle—core model with core motion driven by Earth precession alone, coupled to the mantle only by viscous shear stresses, and with a spherical mantle—core boundary leads to energy dissipation rates on the order of 10⁴ times those necessary for an Earth dynamo. The maximum dissipation model also leads to excessive magnetic field drift rates and to excessive retardation of the Earth's rotation rate. Effects of the mantle—core ellipticity and of magnetic field coupling are briefly discussed and are used to help develop a less than maximum dissipation model also driven by precession alone but using the additional coupling to yield a model more consistent with observed phenomena.     

Tidal friction in historical time and in the remote past

Summary. Section 1 is a commentary on Harold Jeffreys's extension of his treatment of tidal friction in historical time. Sections 2 and 3 are amendments of the theory given in The Earth concerning the Moon's orbit and the Earth's rotation in the remote past.     

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.     

Dislocation-mediated melting of iron and the temperature of the Earth's core

Summary . Dislocation theories of melting provide a possibility of calculating the melting temperature, from first principles, as the temperature at which the free energy of a crystal saturated with dislocations becomes equal to that of the dislocation-free crystal. After a brief review of the physical bases of the dislocation melting theories, Ninomiya's theory is used to determine the melting temperature as well as the volume and entropy of melting and the slope of the melting curve for iron at atmospheric pressure and under conditions prevailing at the Earth's inner core boundary. The necessary parameters (elastic moduli, Grüneisen parameter) are drawn from seismological earth models. We find a melting temperature of the material of the inner core of about 6150 K, independent of shock-wave experiments but in good agreement with them and with extrapolations using Lindemann's law. With usually accepted values of the melting point depression due to light elements in solution, the temperature at the inner core boundary is found to be T _ICB≅ 5000 K. This temperature is compatible with a temperature of the outer core at the core-mantle boundary T _CMB≅ 3800 K. Dislocation melting theories can thus help constrain the temperature profile in the Earth's core.     

Experiments on precessing flows in the Earth's liquid core

Experiments simulating flow in the Earth's liquid core induced by luni-solar precession of the solid mantle indicate, to a first approximation, that the core behaves like a rigidized fluid sphere spinning slower than the mantle and with its spin axis lagging the mantle spin axis in precession. Secondary flow patterns are always present. At low precession rates the fluid sphere is subdivided into a set of cylinders coaxial with the fluid spin axis, the cylinders rotating alternately at slightly faster and slower rates relative to the net retrograde motion of the fluid as a whole. Slow non-axisymmetric columnar wave patterns develop between the differentially rotating cylinders. Axial flows between the spheroidal cavity boundary and the interior are observed. Fluid motion becomes turbulent only at precession rates large enough to cause the fluid spin axis to align nearly with the precession axis. There is no evidence that the Earth's liquid spin axis direction departs more than a fraction of a degree from geographic north. Our observations suggest precession induces a complex variety of laminar flows, including slowly varying and/or periodic patterns, in the Earth's liquid core.     

Wide-scale detection of earthquake waveform doublets and further evidence for inner core super-rotation

We report on more than 100 earthquake waveform doublets in five subduction zones, including an earthquake nest in Bucaramanga, Colombia. Each doublet is presumed to be a pair of earthquakes that repeat at essentially the same location. These doublets are important for studying earthquake physics, as well as temporal changes of the inner core. Particularly, our observation from one South Sandwich Islands (SSI) doublet recorded at station INK in Canada shows an inner core traveltime change of ∼0.1 s over ∼6 yr, confirming the inner-core differential motion occurring beneath Central America. Observations from one Aleutian Islands doublet, recorded at station BOSA in South Africa, and from one Kuril Islands doublet, recorded at station BDFB in Brazil, show an apparent inner core traveltime change of ∼0.1 s over ∼7 yr and ∼6 yr, respectively, providing evidence for the temporal change of inner core properties beneath Central Asia and Canada, respectively. On the other hand, observations from one Tonga–Fiji–Solomon Islands doublet, recorded at station PTGA in Brazil, and from one Bucaramanga doublet, recorded at station WRAB in Australia and station CHTO in Thailand, show no/little temporal change (no more than 0.005 s yr⁻¹, if any) of inner core traveltimes for the three corresponding ray paths for which the path in the inner core is nearly parallel to the equatorial plane. Such a pattern of observations showing both presence and possible absence of inner-core traveltime change can be explained by the geometry and relative directions of ray path, lateral velocity gradient and inner-core particle motion due to an eastward super-rotation of a few tenths of a degree per year.     

A variational principle for the subseismic wave equation

Summary. A variational principle is developed for the subseismic wave equation governing the normal modes of the outer liquid core with frequencies below seismic frequencies (>300 μHz). The calculation of these modes is important both in determining the core contribution to the Earth's dynamical response to tidal and other forces and because their detection at the surface could provide valuable new insight into the density structure of the core, critical to theories of the geomagnetic dynamo. Included as a special case is a variational principle for the Poincaré equation governing inertial oscillations studied in the laboratory by Aldridge and others. This opens up the possibility of 'tracking' laboratory results over to the real Earth.     

Electromagnetic Core-mantle Coupling

i
The work of Bullard (1950) and Rochester (1960) on the geomagnetic westward drift and its effects on the Earth's rotation is extended to investigate the effects of assuming various distributions of electrical conductivity in the mantle. By a proper choice of conductivities, one is able to increase the theoretical value for the tightness of the coupling by a factor of at least six over that afforded by Rochester's model, without sacrificing agreement with observations on the rapidity with which changes in the secular variation are established at the Earth's surface. It is shown that it is reasonable to attribute the observed random changes in the length of the day to perturbations in the electromagnetic coupling.     

Global heat budget, plate tectonics and climatic change

For the past 2000 Ma, the temperature of the Earth's surface has fluctuated around a mean similar to that of today, although individual locations have undergone long-term changes of ∼30°C at different times in different places. Water bodies absorb at least five times as much solar radiation as land surfaces, and ocean currents transport the excess heat absorbed in the tropics towards the poles. Changes in the distribution of land and sea due to plate tectonics explain the major temperature fluctuations (>25°C) around the globe in the last 350 Ma, and are first-order controls. Large-scale changes in ocean currents and thermohaline circulations are probably second-order controls (15–25°C). The Milankovitch orbital cycles are third-order controls producing variations in air temperature of the order of 10°C, while massive volcanic eruptions and changes in carbon dioxide are amongst the fourth-order controls producing minor perturbations (<5°C). The major climatic fluctuations are continuous but regional in effect and not global. Extraterrestrial factors may not cause major changes in climate when viewed from a geological perspective.     

Limitations of Waveform Modelling of Long-Period Seismograms

In recent years the use of synthetic seismograms calculated for radially stratified models has gained increasing popularity as a means of placing further constraints on the velocity structure of the Earth's mantle. Such synthetics do. however, have a number of limitations. At short periods (∼1 s) the amplitudes as well as the wave shapes of travel-time branches are affected by seismograph siting, the structure immediately beneath the seismograph and any laterally heterogeneous structure in the mantle. Later arrivals can also be masked by signal-generated noise and by extended source times. At longer periods (∼-20 s) the larger wavelengths reduce the sensitivity of amplitudes and waveforms to contaminating effects. As a result the use of long-period synthetics can only lead to the resolution of the gross features of the Earth's interior.