On the theory of core-mantle coupling

This article commences by surveying the basic dynamics of Earth's core and their impact on various mechanisms of core-mantle coupling. The physics governing core convection and magnetic field production in the Earth is briefly reviewed. Convection is taken to be a small perturbation from a hydrostatic, “adiabatic reference state” of uniform composition and specific entropy, in which thermodynamic variables depend only on the gravitational potential. The four principal processes coupling the rotation of the mantle to the rotations of the inner and outer cores are analyzed: viscosity, topography, gravity and magnetic field. The gravitational potential of density anomalies in the mantle and inner core creates density differences in the fluid core that greatly exceed those associated with convection. The implications of the resulting “adiabatic torques” on topographic and gravitational coupling are considered. A new approach to the gravitational interaction between the inner core and the mantle, and the associated gravitational oscillations, is presented. Magnetic coupling through torsional waves is studied. A fresh analysis of torsional waves identifies new terms previously overlooked. The magnetic boundary layer on the core-mantle boundary is studied and shown to attenuate the waves significantly. It also hosts relatively high speed flows that influence the angular momentum budget. The magnetic coupling of the solid core to fluid in the tangent cylinder is investigated. Four technical appendices derive, and present solutions of, the torsional wave equation, analyze the associated magnetic boundary layers at the top and bottom of the fluid core, and consider gravitational and magnetic coupling from a more general standpoint. A fifth presents a simple model of the adiabatic reference state.     

On asymptotics of low-frequency oscillations of the liquid core: 2. Frequency dependence of the inertial core-mantle coupling

New, unique information on the inertial and dissipative coupling of the liquid core and the mantle has been retrieved from modern high-precision (radiointerferometer and GPS) data on tidal variations in the rotation velocity and nutation of the Earth. Comparison of theoretical and observed data provided new estimates for the dynamic flattening of the outer liquid and the inner solid cores, mantle quality factor, viscosity of the liquid core, and electromagnetic coupling of the liquid core and the mantle [Molodensky, 2004, 2006]. As was shown in the first part of the paper [Molodensky, 2008] (further referred to as [I]), generation of eddy flows in Proudman-Taylor columns, whose orientation is controlled by the topography of the liquid core-mantle boundary, should be taken into account for correct estimation of the inertial coupling (see formulas (8) and (34) in [I]). The range of periods within which this effect plays a significant role is determined by the decay time of these flows. This time is estimated in the paper for the case where dissipation is related to viscous friction at the core-mantle boundary or with the electromagnetic coupling of the liquid core and the mantle. Because of significant uncertainties in modern data on the viscosity of the liquid core, the magnetic field intensity at the core-mantle boundary, and the electrical conductivity of the lower mantle, the dissipative coupling of the liquid core and the mantle cannot be calculated as yet. However, as shown in the paper, the decay time of eddy flows is connected with the attenuation time of subdiurnal free nutation and with the liquid core viscosity. This enables the estimation of the frequency dependence of the dissipative coupling in a fairly wide range. It is shown that the range of periods for which relations (8) and (34) in [I] are valid encompasses the best-studied length-of-day variations and, therefore, these relations are applicable to analysis of the majority of modern data.     

核幔耦合对地球自由核章动的激发影响

地球自由核章动（FCN）是地幔与液核相互作用的重要动力学现象,其激发机制涉及地表流体层、地幔和地核等圈层之间的耦合,此前研究多利用地表流体层角动量数据单独研究其对FCN的激发,对核幔耦合的影响考虑不足.本文基于角动量守恒理论分析了核幔耦合对FCN周期及振幅的影响,并结合多个大气及海洋角动量函数时间序列首次估算了核幔耦合在FCN激发过程中的贡献.结果表明核幔耦合对FCN周期产生的固定和时变影响对FCN激发的作用均不可忽视,尤其时变影响可达几十个微角秒,对于进一步解释FCN时变特征非常重要;核幔耦合对FCN振幅的直接影响是地表流体层的激发与实测FCN不相符的主要原因,黏滞、电磁和地形等耗散耦合的存在对地表流体的激发振幅有67%左右的减弱效果.     

On the effect of oceanic tides on the inertial coupling between the liquid core and the mantle

The amplitudes and phases of forced nutation and diurnal earth tides depend significantly on the moment of forces between the liquid core and mantle of the Earth, resulting from the differential rotation of the core. The solution to the dynamic problem of rotation of an imperfectly elastic mantle with an imperfectly liquid core and an ocean indicates that the predominant role is played by the so-called core-mantle inertial coupling (related to the effect of hydrodynamic pressure in the liquid core on the ellipsoidal core-mantle boundary). The magnitude of this coupling depends significantly not only on the dynamic flattening of the liquid core but also on the elastic and inelastic properties of the mantle, as well as on the amplitudes and phases of oceanic tides. In this paper, the effects of oceanic tides on the magnitude of inertial coupling between the liquid core and the mantle and on the period and damping decrement of free nearly diurnal nutation are estimated.     

Magnetic probing of planetary interiors

The following general question is addressed: what can be learned about a planetary interior from measurements of the global planetary magnetic field at (or near) its surface? The discussion is placed in the context of Earth, for clarity, but the considerations apply to terrestrial planets in general (so long as the observed magnetism is either predominantly of internal origin, or else external source effects can be successfully filtered out of the observations). Attention is given to the idealized but typical situation of a rotating but spherically symmetric planet containing a highly conducting uniform fluid core surrounded by a nearly insulating rigid mantle, whose conductivity, a function of at most radius only, falls monotonically from its largest value at the base of the planetary mantle to zero at the planetary surface; the largest value of mantle conductivity as well as the mean value for the whole mantle and the mantle conductance are assumed small compared to the corresponding values of the core. Exterior to the planet is vacuum in the sense of an electrically uncharged insulator. The core fluid is inviscid, Boussinesq and gravitationally driven.Complete and perfect observations of either the instantaneous internal vector magnetic field together with its secular variation at a single epoch, or more realistically, the instantaneous internal vector magnetic field alone at two separated epochs are presumed available; the time separation between measurement epochs is long compared the Ohmic diffusion time of the planetary mantle, but small compared to that of the liquid core.Under such circumstances we describe how information about each of the following planetary properties can, in principle (though not without practical difficulty) be retrieved from the observations: (1) depth of the core-mantle boundary (a result of Hide); (2) depth to the current and motion sources responsible for the planetary dynamo; (3) presence or absence of small-scale turbulence in the upper reaches of the core; (4) large-scale horizontal fluid motion at the top of the core; (5) strength of horizontal currents, zonal magnetic fields, Coriolis and Lorentz forces at the top of the core; and (6) current system in the mantle and strength of electromagnetic core-mantle coupling.     

Axial poloidal electromagnetic core-mantle coupling torque: A re-examination for different conductivity and satellite supported geomagnetic field models

We investigate the temporal behaviour of the axial component of the electromagnetic core-mantle coupling torque that is associated with the poloidal part of the geomagnetic field observable at the Earth surface. For its computation, we use different models of the geomagnetic field, expanded into spherical harmonics (Wardinski and Holme, 2006; Sabaka et al., 2004), and the mantle conductivity. The geomagnetic field, which we have to know at the core-mantle boundary for the associated computations, will be inferred from the field at the Earth surface by the non-harmonic field continuation through a conducting mantle shell. The aims of this investigation are (i) to check how sensitive is the computation of the torque with respect to the different geomagnetic field models, (ii) to check its dependence on the spherical harmonic degree n, and (iii) to determine the difference between the mechanical torque derived from the observed length-of-day variations (atmospheric influence subtracted) and the poloidal electromagnetic torque in dependence on the assumed conductivity. To use the non-harmonic field continuation for the torque calculation and to obtain an insight into the influence of the different geomagnetic field models on the EM torques are the major aspects of this paper. grm@gfz-potsdam.de     

Core-mantle coupling and viscoelastic deformations

Writing the angular momentum theorem for the Earth and for its fluid core, we show that there are couplings between the core and the mantle induced by viscomagnetic torque, by external active torque, by topographic torque acting at the core-mantle boundary (CMB) but also by viscoelastic deformations of the CMB which may perturb the axial rotations of the Earth and of the core. We compute these deformations at the CMB induced by the Pleistocenic deglaciation. The time-dependence of inertia tensor perturbations, i.e. the rheology of the mantle, is very important in the calculation of the coupling. Taking into account the passive viscomagnetic torque of tangential traction acting at the CMB, we investigate, for different values and various temporal evolutions of the topographic torque, the perturbations in the rotations of the Earth and of the core induced by the deglaciation, by the constant torque of tidal friction and by the 18.6 year tidal potential. We show that, for these excitation sources, the existence of a constant topographic torque involves the core oscillating with respect to the mantle and thus forbids any large drift of the core with respect to the mantle. However, it seems theoretically possible to have an excitation source with enough energy which involves a shift of the core with respect to the mantle. If the pressure within the fluid core varies with time, the motion of the core with respect to the mantle could be drastically different.     

大地测量在近代发展的内部地球物理学中的作用

我们的目的是说明大地测量学的近代发展如何为解决地球动力学的基本问题提供新的工具。例如:1.全球规模的重力场的知识将给出关于地幔中过渡带位置的信息,或关于地核-地幔交界面的地形的资料;2.在地球表面上的高精度定位将允许测量构造板块的位移,同时将对地震预报作出贡献;3.用新技术所确定的极移的分析,将能阐明其阻尼问题,维持这一运动的能源问题以及地核-地幔耦合等问题。     

Torque balance and energy budget for the precessionally driven dynamo

The feasibility of a precessionally driven dynamo is investigated. The relative orientation of the angular-velocity vectors of the mantle and core and the precession vector of the earth are determined from a torque balance. The core and mantle are acted upon by separate gravitational torques and mutual interaction torques resulting from pressure, viscous and magnetic stresses at the core-mantle interface. The viscous and magnetic torques are determined using the results of a detailed analysis of the Ekman-Hartmann and magnetic diffusion layers generated at the core-mantle interface by the misalignment of the mantle and core angular-velocity vectors. The dissipative torques are found to be weaker by a factor of 10^?4 than those estimated by Malkus (1968) and Stacey (1973), resulting in only 3.5 · 10⁷ W being extracted from the rotational kinetic energy of the earth by these mechanisms. Furthermore, it is found that all of this energy is dissipated in the boundary layers at the core-mantle interface and none is available to drive the geodynamo.     

Visualization of whole mantle propagation of seismic shear energy using normal mode summation

Correct representation of seismic waveforms propagating through the mantle from a 600 km deep earthquake is presented using graphic interpolation between synthetic seismograms computed across a grid of mantle depths and distances. All torsional normal modes with periods above 12 s are summed to create 72,846 seismograms at depths between the surface and the core-mantle boundary. The resulting time snapshots show the manner by which seismic shear energy propagates around the core away from the source.     

Dependence of the Love numbers upon the inner structure of the Earth and comparison of theoretical models with results of measurements

Summary The question this paper is examining is the following: to what extent are the Love numbers dependent on certain characteristics of the inner structure of the Earth? It has been proven — on the basis of calculations carried out by the author-that these quantities are only in a small degree dependent on the density values measured on the surface of the Earth and on the selection of the density function in the mantle of the Earth. On the other hand the value of Love numbersh, k andl is considerably influenced by the assumptions made about the core of the Earth, namely by the position of the boundary between the core and the mantle and by the magnitude of the rigidity coefficient presumed in the core in the vicinity of the core-mantle boundary.The results of the calculations are compared with those mean values of Love numbers obtained from the data of stations operating at different places of the Earth. By reason of this it can be assumed that the core of the Earth has, in the vicinity of the core-mantle boundary, a coefficient of effective rigidity of the order of 10¹⁰ dyn/cm², if the core-mantle boundary is placed at the relative Earth radius of 0.545 from the centre of the Earth.     

Inhomogeneities near the core-mantle boundary evidenced from scattered waves: A review

Statistical properties of small-scale inhomogeneities (wavelengths between 20 and 70 km) near the core-mantle boundary are inferred from scattered core waves. Observations of scattered core waves at large seismic arrays and worldwide networks indicate that the inhomogeneities have a global nature with similar characteristics. However, there may exist a few regions having markedly stronger or weaker strengths. Scattering by volumetric inhomogeneities of about 1% inP-wave velocity in the lower mantle or by about 300 m of topographic relief of the core-mantle boundary can explain the observations. At present it is not possible to rule out either of these two alternatives, or a combination of both.     

Research on the scatter features of the PKiKP/PcP amplitude ratio and the inner core boundary density contrasts beneath Northeast Asia

Quantifying the density contrasts of the Earth's inner core boundary(ICB) is crucial to understand core-mantle coupling and the generation of the geodynamo. The PKiKP/PcP amplitude ratio is commonly used to obtain the density contrast at the ICB, but its applications are limited by scattered observed data. In this study, we selected the PKiKP and PcP phases reflected at the same region of inner-core and core-mantle boundaries beneath Northeast Asia from different earthquakes for the first time, and the observations suggested that the PKiKP/PcP amplitude ratio is widely scattered. We also compared the PKiKP and PcP amplitudes, which demonstrated that the scatter cannot be attributed only to ICB anomalies but might also arise from raypath differences and heterogeneities throughout the crust and mantle. By fitting the observed PKiKP/PcP amplitude ratio, we obtained a density contrast of approximately 0.65 g cm~(-3) and a compressional velocity contrast of approximately 0.87 km s~(-1) at the ICB beneath Northeast Asia. The larger contrast values indicate the possible occurrence of local crystallization occurring at the inner core surface.     

Detecting thermal boundary control in surface flows from numerical dynamos

The geomagnetic field and secular variation exhibit asymmetrical spatial features which are possibly originating from an heterogeneous thermal control of the Earth's lower mantle on the core. The identification of this control in magnetic data is subject to several difficulties, some of which can be alleviated by the use of core surface flow models. Using numerical dynamos driven by heterogeneous boundary heat flux, we confirm that within the parameter space accessible to simulations, time average surface flows obey a simple thermal wind equilibrium between the Coriolis and buoyancy forces, the Lorentz, inertial and viscous forces playing only a secondary role, even for Elsasser numbers significantly larger than 1. Furthermore, we average the models over the duration of three vortex turnovers, and correlate them with a longer time average which fully reveals the signature of boundary heterogeneity. This allows us to quantify the possibility of observing mantle control in core surface flows averaged over a short time period. A scaling analysis is performed in order to apply the results to the Earth's core. We find that three vortex turnovers could represent between 100 and 360 years of Earth time, and that the heat flux heterogeneity at the core-mantle boundary could be large enough to yield an observable signature of thermal mantle control in a time average core surface flow within reach of the available geomagnetic data.     

Mass heterogeneities and convection in the earth's mantle inferred from gravity and core-mantle boundary irregularities

Mass heterogeneities in the earth's mantle are retrieved from the gravity data and the topography of the core-mantle boundary as well as the topography of the earth's surface. A mantle circulation induced by the heterogeneities is modelled by solving the Stokes problem for incompressible Newtonian fluid. The derived models of mantle motions correlate well with the plate tectonics and point at a close relation between the surface tectonic activity and the processes in the vicinity of the core-mantle boundary.     

The amplitude ratioPcP/P and the core-mantle boundary

Summary Using the Haskell matrix formulation, theoretical reflection coefficient curves have been calculated for a multi-layered core-mantle boundary for comparison with observational data. Two cases are considered, first when the shear velocity in the core is equal to zero and second when the core has a finite rigidity. If the velocity contrast is large between the imbedded layer and the mantle, the reflection coefficient curves for the multi-layered medium are irregular in shape as compared to those for two half-spaces, representing the core and the mantle, respectively. The reflection coefficient curves show an oscillatory character if the imbedded layer is thick and has a high velocity contrast.The observational data consist of short-period vertical-component seismograph records ofP andPcP from nuclear explosions in the Aleutian chain, Nevada, Novaya Zemlya, Kazakh and Sahara. Attenuation and geometrical spreading are taken into consideration. Four different models for the quality factorQ are applied to the observational data. The data are found to be much affected by theQ-model used for the corrections.Based on proposedQ-values, a model for the core-mantle boundary is found, characterized by two low-velocity layers at the bottom of the mantle. The thicknesses are 16.10 km (outer layer) and 19.96 km (inner layer), the compressional wave velocities 12.17 km/sec and 10.94 km/sec and the shear wave velocities are 6.29 km/sec and 5.33 km/sec, respectively. A better fit to this model is found when in addition the shear velocity in the outer core is 2.20 km/sec and the density ratio at the core-mantle boundary is 1.07. In other words, the observations favour a layer of finite rigidity in the outer core rather than a fluid one.     

On the coherence between changes in biota and geomagnetic reversals in the Phanerozoic

The data on geomagnetic reversals, organic changes, and lower-mantle plume magmatism in the Phanerozoic are collected and correlated. No direct relationship is revealed between the geomagnetic reversals, plumes, and biozones. However, the frequency of geomagnetic reversals is found to correlate to the frequency of biozonal alterations. We relate this inconsistency to the coupling of the two processes, which are mutually independent, with the long-term changes in the Earth’s rotation. The plumes are formed at the core-mantle boundary and, thus, the reversals should have a different source. We hypothesize that the change in the geomagnetic polarity is due to the nonuniform rotation of the inner core relative to the mantle in combination with the changes in the axial tilt of the Earth’s rotation.     

Core cooling by subsolidus mantle convection

Although vigorous mantle convection early in the thermal history of the Earth is shown to be capable of removing several times the latent heat content of the core, we are able to construct a thermal evolution model of the Earth in which the core does not solidify. The large amount of energy removed from the model Earth's core by mantle convection is supplied by the internal energy of the core which is assumed to cool from an initial high temperature given by the silicate melting temperature at the core-mantle boundary. For the smaller terrestrial planets, the iron and silicate melting temperatures at the core-mantle boundaries are more comparable than for the Earth, and the cores of these planets may not possess enough internal energy to prevent core solidification by mantle convection. Our models incorporate temperature-dependent mantle viscosity and radiogenic heat sources in the mantle. The Earth models are constrained by the present surface heat flux and mantle viscosity. Internal heat sources produce only about 55% of the Earth model's present surface heat flow.     

Axial and equatorial rotations of the Earth's cores associated with the Quaternary ice age

The differential axial and equatorial rotations of both cores associated with the Quaternary glacial cycles were evaluated based on a realistic earth model in density and elastic structures. The rheological model is composed of compressible Maxwell viscoelastic mantle, inviscid outer core and incompressible Maxwell viscoelastic inner core. The present study is, however, preliminary because I assume a rigid rotation for the fluid outer core. In models with no frictional torques at the boundaries of the outer core, the maximum magnitude of the predicted axial rotations of the outer and inner cores amounts to ∼2° year⁻¹ and ∼1° year⁻¹, respectively, but that for the secular equatorial rotations of both cores is ∼0.0001° at most. However, oscillating parts with a period of ∼225 years are predicted in the equatorial rotations for both cores. Then, I evaluated the differential rotations by adopting a time-dependent electromagnetic (EM) torque as a possible coupling mechanism at the core-mantle boundary (CMB) and inner core boundary (ICB). In a realistic radial magnetic field at the CMB estimated from surface magnetic field, the axial and equatorial rotations couple through frictional torques at the CMB, although these rotations decouple for dipole magnetic field model. The differential rotations were evaluated for conductivity models with a conductance of 10⁸ S of the lowermost mantle inferred from studies of nutation and precession of the Earth and decadal variations of length of day (LOD). The secular parts of equatorial rotations are less sensitive to these parameters, but the magnitude for the axial rotations is much smaller than for frictionless model. These models, however, produce oscillating parts in the equatorial rotations of both cores and also in the axial rotations of the whole Earth and outer and inner cores. These oscillations are sensitive to both the magnitude of radial magnetic field at the CMB and the conductivity structure. No sharp isolated spectral peaks are predicted for models with a thin conductive layer (∼200 m) at the bottom of the mantle. In models with a conductive layer of ∼100 km thickness, however, sharp spectral peaks are predicted at periods of ∼225 and ∼25 years for equatorial and axial rotations, respectively, although these depend on the strength of radial magnetic field at the CMB. While the present study is preliminary in modelling the fluid outer core and coupling mechanism at the CMB, the predicted axial rotations of the whole Earth may be important in explaining the observed LOD through interaction between the equatorial and axial rotations.     

Mercury's magnetic field: a thermoelectric dynamo?

Permanent magnetism and conventional dynamo theory are possible but problematic explanations for the magnitude of the Mercurian magnetic field. A new model is proposed in which thermoelectric currents driven by temperature differences at a bumpy core-mantle boundary are responsible for the (unobserved) toroidal field, and the helicity of convective motions in a thin outer core (thickness 10² km) induces the observed poloidal field from the toroidal field. The observed field of 3 × 10⁻⁷ T can be reproduced provided the electrical conductivity of Mercury's semiconducting mantle approaches 10³ Ω⁻¹ m⁻¹. This model may be testable by future missions to Mercury because it predicts a more complicated field geometry than conventional dynamo theories. However, it is argued that polar wander may cause the core-mantle topography to migrate so that some aspects of the rotational symmetry may be reflected in the observed field.