Core precession: flow structures and energy

Experiments using a precessing liquid-filled oblate spheroid with ellipticity ( a − b )/ a =1/400 extend and clarify earlier research. They yield flow data useful for estimating flows in the Earth's liquid core. Observed flows illustrate and confirm a nearly rigid liquid sphere with retrograde drift and lagging a cavity (mantle) axis in precession. The similarities of the observed lag angle with that computed for a rigid sphere, and earlier energy dissipation research both support the use of a rigid sphere analytical model to predict energy dissipation and first-order flow within the core–mantle boundary (CMB). Second-order boundary layer and interior cylindrical flow structures also are photographed and measured. Interior flows are never turbulent or unstable at near-Earth parameters, although complex and transient flow patterns are observed within the boundary layer. Other mechanisms proposed to explain net heat loss from the Earth and maintenance of the geodynamo typically require acceptance of some critical but unproven premise. Precession and CMB configuration are known with certainty and precision. Analytical difficulties have been the obstacle. Experiments illustrate the consequences of precession and ellipticity, provide criteria for validating analytical and numerical models, and may yield direct knowledge of the Earth's deep interior with careful scaling.     

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.     

Magnetic and viscous coupling at the core–mantle boundary: inferences from observations of the Earth's nutations

Dissipative core–mantle coupling is evident in observations of the Earth's nutations, although the source of this coupling is uncertain. Magnetic coupling occurs when conducting materials on either side of the boundary move through a magnetic field. In order to explain the nutation observations with magnetic coupling, we must assume a high (metallic) conductivity on the mantle side of the boundary and a rms radial field of 0.69 mT. Much of this field occurs at short wavelengths, which cannot be observed directly at the surface. High levels of short-wavelength field impose demands on the power needed to regenerate the field through dynamo action in the core. We use a numerical dynamo model from the study of Christensen & Aubert (2006) to assess whether the required short-wavelength field is physically plausible. By scaling the numerical solution to a model with sufficient short-wavelength field, we obtain a total ohmic dissipation of 0.7–1 TW, which is within current uncertainties. Viscous coupling is another possible explanation for the nutation observations, although the effective viscosity required for this is 0.03 m² s⁻¹ or higher. Such high viscosities are commonly interpreted as an eddy viscosity. However, physical considerations and laboratory experiments limit the eddy viscosity to 10⁻⁴ m² s⁻¹, which suggests that viscous coupling can only explain a few percent of the dissipative torque between the core and the mantle.     

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.     

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.     

Electromagnetic and Viscous Damping of Core Oscillations

Summary. A detailed investigation of the interaction of core oscillations with the main magnetic field is made. Selection rules and interaction coefficients are obtained. Although the strengths of the interactions vary over nearly three orders of magnitude, all are extremely weak. Ohmic and viscous dissipation are considered in parallel calculations. It is found that core oscillations are almost unaffected by either damping mechanism. Most of the Ohmic dissipation takes place within a few skin depths of the boundary. Viscous dissipation is estimated only for the body of the liquid core. Even though it is expected to be much larger near the boundaries, it would seem unlikely to be sufficiently large to be of any significance. Our conclusion is that core oscillations are virtually free of damping. In a subadiabatic core, gravitational oscillations might be persistent enough to sustain an oscillatory geodynamo.     

Geomagnetic field analysis - I. Stochastic inversion

Summary. Most of the Earth's magnetic field and its secular change originate in the core. Provided the mantle can be treated as an electrical insulator, stochastic inversion enables surface observations to be analysed for the core field. A priori information about the variation of the field at the core boundary leads to very stringent conditions at the Earth's surface. The field models are identical with those derived from the method of harmonic splines (Shure, Parker & Backus) provided the a priori information is specified appropriately.
The method is applied to secular variation data from 106 magnetic observatories. Model predictions for fields at the Earth's surface have error estimates associated with them that appear realistic. For plausible choices of a priori information the error of the field at the core is unbounded, but integrals over patches of the core surface can have finite errors. The hypothesis that magnetic fields are frozen to the core fluid implies that certain integrals of the secular variation vanish. This idea is tested by computing the integrals and their standard and maximum errors. Most of the integrals are within one standard deviation of zero, but those over the large patches to the north and south of the magnetic equator are many times their standard error, because of the dominating influence of the decaying dipole. All integrals are well within their maximum error, indicating that it will be possible to construct core fields, consistent with frozen flux, that satisfy the observations.     

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.     

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.     

Properties of iron at the Earth's core conditions

Summary. The phase diagram of iron up to 330 GPa is solved using the experimental data of static high pressure (up to 11 GPa) and the experimental data of shock wave data (up to 250 GPa). A solution for the highest triple point is found ( P = 280 GPa and T = 5760 K) by imposing the thermodynamic constraints of triple points. This pressure of the triple point is less than the pressure of the inner core–outer core boundary of the Earth. These results indicate that the density of iron at the inner core–outer core boundary pressure is close to 13 g cm⁻³, which lies close to the seismic solutions of the Earth at that pressure. It is thus concluded that the Earth's inner core is very likely to be virtually pure iron in its hexagonal close packed (hcp) phase.
It is shown that four properties of the Earth's inner core determined from seismology are close in value to the corresponding properties of hcp iron at inner core conditions: density, bulk modulus, longitudinal velocity, and Poisson's ratio. The density–pressure profile of hcp iron at inner core conditions matches the density–pressure profile of the inner core as determined by seismic methods, within the spread of values given by recent seismic models.
This indicates that the Earth is slowly cooling, the Earth's inner core is growing by crystallization, and the impurities of the core are concentrated in the outer core. The calculated temperature at the Earth's centre is 6450 K.     

Geomagnetic effects of the Earth's ellipticity

We analyse the external field generated by a uniform distribution of magnetic susceptibility contained in an oblate spheroidal shell when it is magnetized by an internal magnetic field of arbitrary complexity. The situation is more relevant to the Earth than that of a spherical shell considered by Runcorn (1975a ) (in the context of lunar magnetism), because of the larger flattening of the Earth than that of the Moon. We find that, to first order in the susceptibility, each internal harmonic in a spheroidal harmonic expansion of the magnetic potential generates just one non-vanishing external field coefficient, unlike in the spherical case when all harmonics vanish identically. The field generated is proportional to the susceptibility, thickness of the shell and square of the Earth's eccentricity, and hence it appears that this field amplification mechanism will be very ineffective for the Earth.     

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.     

Spectral-finite-element approach to two-dimensional electromagnetic induction in a spherical earth

We present a spectral-finite-element approach to the 2-D forward problem for electromagnetic induction in a spherical earth. It represents an alternative to a variety of numerical methods for 2-D global electromagnetic modelling introduced recently (e.g. the perturbation expansion approach, the finite difference scheme). It may be used to estimate the effect of a possible axisymmetric structure of electrical conductivity of the mantle on surface observations, or it may serve as a tool for testing methods and codes for 3-D global electromagnetic modelling. The ultimate goal of these electromagnetic studies is to learn about the Earth's 3-D electrical structure.
Since the spectral-finite-element approach comes from the variational formulation, we formulate the 2-D electromagnetic induction problem in a variational sense. The boundary data used in this formulation consist of the horizontal components of the total magnetic intensity measured on the Earth's surface. In this the variational approach differs from other methods, which usually use spherical harmonic coefficients of external magnetic sources as input data. We verify the assumptions of the Lax-Milgram theorem and show that the variational solution exists and is unique. The spectral-finite-element approach then means that the problem is parametrized by spherical harmonics in the angular direction, whereas finite elements span the radial direction. The solution is searched for by the Galerkin method, which leads to the solving of a system of linear algebraic equations. The method and code have been tested for Everett & Schultz's (1995) model of two eccentrically nested spheres, and good agreement has been obtained.     

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.     

Numerical models of mantle convection with secular cooling

A 2-D time-dependent finite-difference numerical model is used to investigate the thermal character and evolution of a convecting layer which is cooling as it convects. Two basic cooling modes are considered: in the first, both upper and lower boundaries are cooled at the same rate, while maintaining the same temperature difference across the layer; in the second, the lower boundary temperature decreases with time while the upper boundary temperature is fixed at 0°C. The first cooling mode simulates the effects of internal heating while the second simulates planetary cooling as mantle convection extracts heat from, and thereby cools, the Earth's core. The mathematical analogue between the effects of cooling and internal heating is verified for finite-amplitude convection. It is found that after an initial transient period the central core of a steady but vigorous convection cell cools at a constant rate which is governed by the rate of cooling of the boundaries and the viscosity structure of the layer. For upper-mantle models the transient stage lasts for about 30 per cent of the age of the Earth, while for the whole mantle it lasts for longer than the age of the Earth. Consequently, in our models the bulk cooling of the mantle lags behind the cooling of the core-mantle boundary. Models with temperature-dependent viscosity are found to cool in the same manner as models with depth-dependent viscosity; the rate of cooling is controlled primarily by the horizontally averaged variation of viscosity with depth. If the Earth's mantle cools in a similar fashion, secular cooling of the planet may be insensitive to lateral variations of viscosity.     

A Hamiltonian approach to dissipative phenomena between the Earth's mantle and core, and effects on free nutations

The Hamiltonian formalism was recently applied by Getino (1995a,b) for the study of the rotation of a non-rigid earth with a heterogeneous and stratified liquid core. That earth model is generalized here by including the effect of the dissipation arising from the mantle-core interaction, using a model similar to that of Sasao, Okubo & Saito (1980), which includes both viscous and electromagnetic coupling. First, a solution for the free nutations is obtained following a classical approach, which in our opinion is more familiar to most of the readers than the Hamiltonian treatment. This solution provides a theoretical basis clear enough to study both the qualitative and quantitative effects of the dissipations considered in the hypotheses. The main qualitative features are, besides the delays, that the free core nutation (FCN) suffers an exponential damping, while the chandler wobble (CW) is not damped at first order, by the dissipation considered. The numerical values obtained for the complex compliances agree with the most recent experimental computations.
Next, the problem is studied under a Hamiltonian formalism, and a solution equivalent to the above is obtained. Besides its interest from a theoretical point of view, this formalism is necessary in order to apply canonical perturbation methods in order to obtain analytical nutation series.     

Maximum entropy regularization of time-dependent geomagnetic field models

We incorporate a maximum entropy image reconstruction technique into the process of modelling the time-dependent geomagnetic field at the core–mantle boundary (CMB). In order to deal with unconstrained small lengthscales in the process of inverting the data, some core field models are regularized using a priori quadratic norms in both space and time. This artificial damping leads to the underestimation of power at large wavenumbers, and to a loss of contrast in the reconstructed picture of the field at the CMB. The entropy norm, recently introduced to regularize magnetic field maps, provides models with better contrast, and involves a minimum of a priori information about the field structure. However, this technique was developed to build only snapshots of the magnetic field. Previously described in the spatial domain, we show here how to implement this technique in the spherical harmonic domain, and we extend it to the time-dependent problem where both spatial and temporal regularizations are required. We apply our method to model the field over the interval 1840–1990 from a compilation of historical observations. Applying the maximum entropy method in space—for a fit to the data similar to that obtained with a quadratic regularization—effectively reorganizes the magnetic field lines in order to have a map with better contrast. This is associated with a less rapidly decaying spectrum at large wavenumbers. Applying the maximum entropy method in time permits us to model sharper temporal changes, associated with larger spatial gradients in the secular variation, without producing spurious fluctuations on short timescales. This method avoids the smearing back in time of field features that are not constrained by the data. Perspectives concerning future applications of the method are also discussed.     

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.     

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.