Dynamos with weakly convecting outer layers: implications for core-mantle boundary interaction

Convection in the Earth's core is driven much harder at the bottom than the top. This is partly because the adiabatic gradient steepens towards the top, partly because the spherical geometry means the area involved increases towards the top, and partly because compositional convection is driven by light material released at the lower boundary and remixed uniformly throughout the outer core, providing a volumetric sink of buoyancy. We have therefore investigated dynamo action of thermal convection in a Boussinesq fluid contained within a rotating spherical shell driven by a combination of bottom and internal heating or cooling. We first apply a homogeneous temperature on the outer boundary in order to explore the effects of heat sinks on dynamo action; we then impose an inhomogeneous temperature proportional to a single spherical harmonic Y ₂² in order to explore core-mantle interactions. With homogeneous boundary conditions and moderate Rayleigh numbers, a heat sink reduces the generated magnetic field appreciably; the magnetic Reynolds number remains high because the dominant toroidal component of flow is not reduced significantly. The dipolar structure of the field becomes more pronounced as found by other authors. Increasing the Rayleigh number yields a regime in which convection inside the tangent cylinder is strongly affected by the magnetic field. With inhomogeneous boundary conditions, a heat sink promotes boundary effects and locking of the magnetic field to boundary anomalies. We show that boundary locking is inhibited by advection of heat in the outer regions. With uniform heating, the boundary effects are only significant at low Rayleigh numbers, when dynamo action is only possible for artificially low magnetic diffusivity. With heat sinks, the boundary effects remain significant at higher Rayleigh numbers provided the convection remains weak or the fluid is stably stratified at the top. Dynamo action is driven by vigorous convection at depth while boundary thermal anomalies dominate in the upper regions. This is a likely regime for the Earth's core.     

Transitions to chaotic thermal convection in a rapidly rotating spherical fluid shell

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell heated from below and within have been carried out with a nonlinear, three-dimensional, time-dependent pseudospectral code. The investigated phenomena include the sequence of transitions to chaos and the differential mean zonal rotation. At the fixed Taylor number T _a =10⁶ and Prandtl number Pr=1 and with increasing Rayleigh number R, convection undergoes a series of bifurcations from onset of steadily propagating motions SP at R=R _c = 13050, to a periodic state P, and thence to a quasi-periodic state QP and a non-periodic or chaotic state NP. Examples of SP, P, QP, and NP solutions are obtained at R = 1.3R _c , R = 1.7 R _c , R = 2R _c , and R = 5 R _c , respectively. In the SP state, convection rolls propagate at a constant longitudinal phase velocity that is slower than that obtained from the linear calculation at the onset of instability. The P state, characterized by a single frequency and its harmonics, has a two-layer cellular structure in radius. Convection rolls near the upper and lower surfaces of the spherical shell both propagate in a prograde sense with respect to the rotation of the reference frame. The outer convection rolls propagate faster than those near the inner shell. The physical mechanism responsible for the time-periodic oscillations is the differential shear of the convection cells due to the mean zonal flow. Meridional transport of zonal momentum by the convection cells in turn supports the mean zonal differential rotation. In the QP state, the longitudinal wave number m of the convection pattern oscillates among m = 3,4,5, and 6; the convection pattern near the outer shell has larger m than that near the inner shell. Radial motions are very weak in the polar regions. The convection pattern also shifts in m for the NP state at R = 5R _c , whose power spectrum is characterized by broadened peaks and broadband background noise. The convection pattern near the outer shell propagates prograde, while the pattern near the inner shell propagates retrograde with respect to the basic rotation. Convection cells exist in polar regions. There is a large variation in the vigor of individual convection cells. An example of a more vigorously convecting chaotic state is obtained at R = 50R _c . At this Rayleigh number some of the convection rolls have axes perpendicular to the axis of the basic rotation, indicating a partial relaxation of the rotational constraint. There are strong convective motions in the polar regions. The longitudinally averaged mean zonal flow has an equatorial superrotation and a high latitude subrotation for all cases except R = 50R _c , at this highest Rayleigh number, the mean zonal flow pattern is completely reversed, opposite to the solar differential rotation pattern.     

Kinematic stationary geodynamo models with separated toroidal and poloidal motions

Abstract

This paper is concerned with a three-dimensional spherical model of a stationary dynamo that consists of a convective layer with a simple poloidal flow of the S^2c ₂ kind between a rotating inner body core and solid outer shell. The rotation of the inner core and the outer shell means that there are regions of concentrated shear or differential rotation at the convective layer boundaries. The induction equation for the inside of the convective layer was solved numerically by the Bullard-Gellman method, the eigenvalue of the problem being the magnetic Reynolds number of the poloidal flow (R _M2) and it was assumed that the magnetic Reynolds number of the core (R _M1) and of the shell (R _M3) were prescribed parameters. Hence R _M2 was studied as a function of R _M1 and R _M3, along with the orientation of the rotation axis, the radial dependence of the poloidal velocity and the relative thickness of the layers for the three different situations, (i) the core alone rotating, (ii) the shell alone rotating and (iii) the core and the shell rotating together. In all three cases it was found that, at definite orientations of the rotation axis, there is a good convergence of both the eigenvalues and the eigenfunctions of the problem as the number of spherical harmonics used to represent the problem increases. For R _M1 =R _M3= 10³, corresponding to the westward drift velocity and the parameters of the Earth's core, the critical values of R _M2 are found to be three orders of magnitude lower than R _M1, R _M3 so that the poloidal flow velocity sufficient for maintaining the dynamo process is 10-20 m/yr. With only the core or the shell rotating, the velocity field generally differs little from the axially symmetric case. However, for R _M2 (or R _M3) lying in the range 10² to 10⁵, the self-excitation condition is found to be of the form R _M2˙R ^½ _M1=constant (or R _M2˙R^½ _M3=constant) and the solution does not possess the properties of the Braginsky near-axisymmetric dynamo. We should expect this, in particular, in the Braginsky limit R _M2˙R^?½; _M1=constant.

An analysis of known three-dimensional dynamo models indicates the importance of the absence of mirror symmetry planes for the efficient generation of magnetic fields.     

Generation of magnetic fields by convection in a rotating sphere,I

The magnetohydrodynamic dynamo problem is solved for an electrically conducting spherical fluid shell with spherically symmetric distributions of gravity and heat sources. The dynamics of motions generated by thermal buoyancy are dominated by the effects of rotation of the fluid shell. Dynamos are found for low and intermediate values of the Taylor number, T ? 10⁵, if the scale of the nonaxisymmetric component of the velocity field is sufficiently small. The generation of magnetic fields of quadrupolar symmetry is preferred at Rayleigh numbers close to the critical value R_c for onset of convection. As the Rayleigh number increases, the generation of dipolar magnetic fields becomes preferred.     

Finite prandtl number convection in spherical shells

Abstract

Finite amplitude convection in spherical shells with spherically symmetric gravity and heat source distribution is considered. The nonlinear problem of three-dimensional convection in shells with stress-free and isothermal boundaries is solved by expanding the dependent variables in terms of powers of the amplitude of convection. The preferred mode of convection is determined by a stability analysis in which arbitrary infinitesimal disturbances are superimposed on the steady solutions. The shell is assumed to be thick and only shells for which the ratio ζ of outer radius to inner radius is 2 or 3 are considered. Three cases, two of which lead to a self adjoint problem, are treated in this paper. The stable solutions are found to be l=2 modes for ζ=3 where l is the degree of the spherical harmonics and an l=3 non-axisymmetric mode which exhibits the symmetry of a tetrahedron for ζ=2. These stable solutions transport the maximum amount of heat. The Prandtl number dependence of the heat transport is computed for the various solutions analyzed in the paper.     

A three-dimensional,iterative mapping procedure for the implementation of an ionosphere-magnetosphere anisotropic Ohm’s law boundary condition in global magnetohydrodynamic simulations

The mathematical formulation of an iterative procedure for the numerical implementation of an ionosphere-magnetosphere (IM) anisotropic Ohm’s law boundary condition is presented. The procedure may be used in global magnetohydrodynamic (MHD) simulations of the magnetosphere. The basic form of the boundary condition is well known, but a well-defined, simple, explicit method for implementing it in an MHD code has not been presented previously. The boundary condition relates the ionospheric electric field to the magnetic field-aligned current density driven through the ionosphere by the magnetospheric convection electric field, which is orthogonal to the magnetic field B, and maps down into the ionosphere along equipotential magnetic field lines. The source of this electric field is the flow of the solar wind orthogonal to B. The electric field and current density in the ionosphere are connected through an anisotropic conductivity tensor which involves the Hall, Pedersen, and parallel conductivities. Only the height-integrated Hall and Pedersen conductivities (conductances) appear in the final form of the boundary condition, and are assumed to be known functions of position on the spherical surface R=R₁ representing the boundary between the ionosphere and magnetosphere. The implementation presented consists of an iterative mapping of the electrostatic potential , the gradient of which gives the electric field, and the field-aligned current density between the IM boundary at R=R₁ and the inner boundary of an MHD code which is taken to be at R₂>R₁. Given the field-aligned current density on R=R₂, as computed by the MHD simulation, it is mapped down to R=R₁ where it is used to compute by solving the equation that is the IM Ohm’s law boundary condition. Then is mapped out to R=R₂, where it is used to update the electric field and the component of velocity perpendicular to B. The updated electric field and perpendicular velocity serve as new boundary conditions for the MHD simulation which is then used to compute a new field-aligned current density. This process is iterated at each time step. The required Hall and Pedersen conductances may be determined by any method of choice, and may be specified anew at each time step. In this sense the coupling between the ionosphere and magnetosphere may be taken into account in a self-consistent manner.     

On the onset of thermal convection in slowly rotating fluid shells

Abstract

The problem of the removal of the degeneracy of the patterns of convective motion in a spherically symmetric fluid shell by the effects of rotation is considered. It is shown that the axisymmetric solution is preferred in sufficiently thick shells where the minimum Rayleigh number corresponds to degree l = 1 of the spherical harmonics. In all cases with l > 1 the solution described by sectional spherical harmonics Y^l _l (θ,φ) is preferred.     

The Influence of a Homogeneous Magnetic Field on the Ekman and Stewartson Layers

This contribution is aimed at a comparison of two different methods of how to deal with the solid inner core in geodynamo models. The first method, based on a direct application of the non-slip boundary conditions, was frequently used in the past. The second one, developed by the authors of the present paper, is based on an advanced analytical solution within the boundary layers and consequent formulation of new boundary conditions on the flow in the volume of the outer core. As an example we have used the results obtained by Hollerbach (1997) in the study of the influence of an imposed axial magnetic field on the fluid flow in a differentially rotating spherical shell. In the case of a weak imposed magnetic field, our solutions are very similar to those of Hollerbach. This non-trivial correspondence confirms the correctness of both methods, which are different not only in the formulation of boundary conditions, but also in the numerical methods: whereas Hollerbach used spectral methods, our computer code is based on finite differences. The influence of the conductivity of the inner core on the fluid flow was also studied.     

横向黏度变化的全地幔对流应力场初步研究

将地幔地震波速度异常转换为地幔横向黏度变化(达到3个数量级),在球坐标系下计算了瑞雷数为106、上边界为刚性、下边界为应力自由等温边界条件下的岩石层底部的地幔对流极型和环型应力场.结果表明,地幔对流极型应力场与地表大尺度构造具有良好的对应关系:俯冲带和碰撞带的应力呈现挤压状态,而洋中脊处的应力则呈现拉张状态.地幔对流环...     

Shear flow instabilities in the earth’s plasma sheet region

Shear flow instability arising from the velocity shear between the inner and the outer central plasma sheet regions is studied by treating the plasma as compressible. Based on the linearized MHD equations, dispersion relations for the surface wave modes occurring at the boundary of the inner central plasma sheet (ICPS) and the outer central plasma sheet (OCPS) are derived. The growth rates and the eigenmode frequencies are obtained numerically. Three data sets consisting of parameters relevant to the earth’s magnetotail are considered. The plasma sheet region is found to be stable for constant plasma flows unless M_A>9.6, where M_A is the Alfvén Mach number in the ICPS. However, for a continuously varying flow velocity profile in the ICPS, the instability is excited for M_A\geq1.4. The excited modes have oscillation periods of 2–10 min and 1.5–6 s, and typical transverse wavelengths of 30–100 R_E and 0.5–6 R_E for data sets 1 and 2 (i.e., case of no neutral sheet) respectively. For the data set 3, which corresponds to a neutral sheet at the center of the plasma sheet, the excited oscillations have periods of 2 s-1 min with transverse wavelengths of 0.02–1 R_E.     

Evolution of non-linear α 2-dynamos and taylor's constraint

A key non-linear mechanism in a strong-field geodynamo is that a finite amplitude magnetic field drives a flow through the Lorentz force in the momentum equation and this flow feeds back on the field-generation process in the magnetic induction equation, equilibrating the field. We make use of a simpler non-linear?α?²-dynamo to investigate this mechanism in a rapidly rotating fluid spherical shell. Neglecting inertia, we use a pseudo-spectral time-stepping procedure to solve the induction equation and the momentum equation with no-slip velocity boundary conditions for a finitely conducting inner core and an insulating mantle. We present calculations for Ekman numbers (E) in the range 2.5× 10^?3 to 5.0× 10^?5, for?α?=α ₀cos?θ?sin?π?(r?r_i ) (which vanishes on both inner and outer boundaries). Solutions are steady except at lower E and higher values of?α?₀. Then they are periodic with a reversing field and a characteristic rapid increase then equally rapid decrease in magnetic energy. We have investigated the mechanism for this and shown the influence of Taylor's constraint. We comment on the application of our findings to numerical hydrodynamic dynamos.     

Chaotic axisymmetrical spherical convection and large-scale mantle circulation

This paper presents a study of high Rayleigh number (up to 200 times supercritical) axisymmetrical convection in a spherical shell with an aspect ratio relevant for the Earth's lower mantle. Both bottom-heated and internal heated cases have been considered. Computations have been carried out for an infinite Prandtl number isoviscous fluid with free slip isothermal boundary conditions. The first part of the paper is devoted to the influence of the resolution on the accuracy of the numerical results. It is shown that the resolution strongly influences the onset of time dependence. Recent methods of non-linear physics have been used to prove that the time dependence and the chaotic behaviors of the solutions are real ones. From these results we can confirm that convection is chaotic, in this particular geometry, even for Rayleigh numbers 200 times critical. Aperiodic boundary layer instabilities are found to be incapable of breaking up the large-scale flow, owing to the shear of the global circulation. Spectral analysis of the power associated with the thermal anomalies shows that there is an upward cascade of energy, due to small-scale chaotic instabilities, from l = 2 to l = 4–6 at the bottom boundary, in agreement with new seismic observations at the core-mantle boundary [1–3].     

Calculation of diffusivities of passive fields in turbulent media

Abstract

The exact numerical and approximate analytical solutions of the simplest nonlinear integral equation with second order nonlinearity for the averaged Green function are presented. It is assumed that the turbulence is stationary, homogeneous, isotropic and incompressible. Numerous examples of turbulent spectra are considered (peak-like spectrum, spectra of Kolmogorov's type with different forms of “pumping” regions, stepwise spectra etc.). Special emphasis is given to investigating the case of so called “frozen” turbulence when the parameter ξ =u ₀τ/R→∞ where uτ₀,R ₀ are characteristic velocity, lifetime and space scale of turbulent pulsations, respectively. It is shown that these solutions allow us to calculate the turbulent diffusivities accurately for arbitrary spectra with any values of the parameter ξ. The results take into account the possible helicity of turbulence concerned only with scalar passive fields (number density and temperature).     

The patterns of high-degree thermal free convection and its features in a spherical shell

Introduction Mantle convection is thought to be the most effective way of heat transportation in the earth and the source of driving lithospheric plates (Elasser, 1971). The velocity field of plate motion can be split into poloidal and toroidal parts, which are corresponding to the vertical (i.e. radial) and horizontal motions, respectively, in global model. The toroidal component is manifested in the existences of transform faults and the poloidal one in the presences of convergence and diver…     

A Grid-Spectral Method of the Solution of the 3D Kinematic Geodynamo with the Inner Core

A 3D kinematic geodynamo model in a sphere with the conductive solid inner core is considered. The 3D magnetic field and velocity field are resolved in the physical space for r- and -coordinates, whereas the sin- and cos-decomposition is applied to the -coordinate. The additional boundary conditions for the case of non-zero velocity field on the boundaries of the liquid spherical shell and for different magnetic diffusivities of the inner and outer core are applied. The computer code was tested by free decay mode solutions and comparisons were made also with results reported by other authors. This work is a part of a project to study 3D inviscid geodynamo models.     

TIME DOMAIN ELECTROMAGNETIC RESPONSE OF A SHIELDED CONDUCTOR*

Velekin and Bulgakov (1967) in an interesting model experiment while studying the transient electromagnetic response of a conductive sphere placed below a thin conductive sheet found that at the earlier stages of the transience, the composite system response corresponded to the response due to the overlying sheet alone and at the later stages, it corresponded to that of the sphere alone. To examine whether such a separation of responses due to individual components can be analytically studied and applied to other source configurations, we have analyzed an idealized model consisting of two spherical shells. We find that in corroboration with the above results, the general nature of the curve consists of two humps representing the responses dominated by the outer shell and the inner shell respectively. In addition, however, we find that the two humps gradually disappear to yield a smooth decay curve for increasing values of the ratio σd₁b/σ₂d₂a (where σ₁, σ₂ are the conductivities, d₁, d₂ are the thicknesses of the outer and inner shells respectively, and b and a are their respective distances from the centre) and the effect of inner shell on the composite system response is considerably reduced.     

Necessary conditions for the magnetohydrodynamic dynamo

Abstract

A magnetohydrodynamic, dynamo driven by convection in a rotating spherical shell is supposed to have averages that are independent of time. Two cases are considered, one driven by a fixed temperature difference R and the other by a given internal heating rate Q. It is found that when q, the ratio of thermal conductivity to magnetic diffusivity, is small, R must be of order q ^?4/3 and Q of order q ^?2 for dynamo action to be possible; q is small in the Earth's core, so it is hoped that the criteria will prove useful in practical as well as theoretical studies of dynamic dynamos. The criteria can be further strengthened when the ohmic dissipation of the field is significant in the energy balance. The development includes the derivation of two necessary conditions for dynamo action, both based on the viscous dissipation rate of the velocity field that drives the dynamo.     

Effect of the oceanic lithosphere velocity on free convection in the asthenosphere beneath mid-ocean ridges

The paper presents results obtained in experiments on a horizontal layer heated from below in its central part and cooled from above; the layer models the oceanic asthenosphere. Flow velocity and temperature profiles are measured and the flow structure under boundary layer conditions is determined (at Rayleigh numbers Ra > 5 × 10⁵). The flow in the core of a plane horizontal layer heated laterally and cooled from above develops under conditions of a constant temperature gradient averaged over the layer thickness. The flow core is modeled by a horizontal layer with a moving upper boundary and with adiabatic bounding surfaces under conditions of a constant horizontal gradient of temperature. Exact solutions of free convection equations are found for this model in the Boussinesq approximation. Model results are compared with experimental data. Temperature and flow velocity ranges are determined for the boundary layer regime. Based on the experimental flow velocity profiles, an expression is found for the flow velocity profile in a horizontal layer with a mobile upper boundary heated laterally and cooled from above. Free convection velocity profiles are obtained for the asthenosphere beneath a mid-ocean ridge (MOR) with a mobile lithosphere. An expression is obtained for the tangential stress at the top of the asthenosphere beneath an MOR and the total friction force produced by the asthenospheric flow at the asthenosphere-lithosphere boundary is determined.     

Penetrative convection in the upper ocean due to surface cooling

Abstract

A new non-linear model of mixing and convection based on a modelling of two buoyant interacting fluids is applied to penetrative convection in the upper ocean due to surface cooling. In view of simple algebra, the model is one-dimensional. Dissipation is included, but no mean shear is present. A non-similar analytical solution is found in the case of a well-mixed layer bounded below by a sharp thermocline treated as a boundary layer. This solution is valid if the Richardson number, R _i , defined as the ratio of the total mixed-layer buoyancy to a characteristic rms vertical velocity, is much greater than unity. The model predicts a deepening rate proportional to R _i ^?3/4. The thermocline remains of constant thickness, and the ratio thermocline thickness to mixed-layer depth decreases as R _i ^?3/4 as the mixed layer deepens. If the surface flux is constant, the mixed-layer depth increases with time as t ^½. The vertical structure throughout the mixed layer and thermocline is given by the analytical solution, and vertical profiles of mean temperature and vertical fluxes are plotted. Computed profiles and available laboratory data agree remarkably well. Moreover, the accuracy of the simple analytical results presented here is comparable to that of sophisticated turbulence numerical models.     

Numerical simulations of mantle convection: Time-dependent,three-dimensional,compressible, spherical shell

Abstract

We describe nonlinear time-dependent numerical simulations of whole mantle convection for a Newtonian, infinite Prandtl number, anelastic fluid in a three-dimensional spherical shell for conditions that approximate the Earth's mantle. Each dependent variable is expanded in a series of 4,096 spherical harmonics to resolve its horizontal structure and in 61 Chebyshev polynomials to resolve its radial structure. A semiimplicit time-integration scheme is used with a spectral transform method. In grid space there are 61 unequally-spaced Chebyshev radial levels, 96 Legendre colatitudinal levels, and 192 Fourier longitudinal levels. For this preliminary study we consider four scenarios, all having the same radially-dependent reference state and no internal heating. They differ by their radially-dependent linear viscous and thermal diffusivities and by the specified temperatures on their isothermal, impermeable, stress-free boundaries. We have found that the structure of convection changes dramatically as the Rayleigh number increases from 10⁵ to 10⁶ to 10⁷. The differences also depend on how the Rayleigh number is increased. That is, increasing the superadiabatic temperature drop, δT, across the mantle produces a greater effect than decreasing the diffusivities. The simulation with a Rayleigh number of 10⁷ is approximately 10,000 times critical, close to estimates of that for the Earth's mantle. However, although the velocity structure for this highest Rayleigh number scenario may be adequately resolved, its thermodynamic structure requires greater horizontal resolution. The velocity and thermodynamic structures of the scenarios at Rayleigh numbers of 10⁵ and 10⁶ appear to be adequately resolved. The 10⁵ Rayleigh number solution has a small number of broad regions of warm upflow embedded in a network of narrow cold downflow regions; whereas, the higher Rayleigh number solutions (with large δT) have a large number of small hot upflow plumes embedded in a broad weak background of downflow. In addition, as would be expected, these higher Rayleigh number solutions have thinner thermal boundary layers and larger convective velocities, temperatures perturbations, and heat fluxes. These differences emphasize the importance of developing even more realistic models at realistic Rayleigh numbers if one wishes to investigate by numerical simulation the type of convection that occurs in the Earth's mantle.