火星磁场和行星发电机理论

火星磁场是火星主要观测物理场之一,火星磁场研究对火星探索具有重要的科学意义.本文侧重介绍火星磁场探测所取得的主要成果,在介绍火星行星磁场起源的行星发电机理论的基础上,重点讨论了动力学机制、起始时间、停止的原因等关键性问题,并指出了研究中存在的一些问题.     

A multi-layered kinematic dynamo model: implications of a stratified upper layer in the Earth's core

It has been suggested that there exists a stably stratified electrically conducting layer at the top of the Earth's outer fluid core and that lateral temperature gradients in the lower mantle is capable of a driving thermal-wind-type flow near the core–mantle boundary. We investigate how such a flow in a stable layer could influence the geomagnetic field and the geodynamo using a very simple two-dimensional kinematic dynamo model in Cartesian geometry. The dynamo has four layers representing the inner core, convecting lower outer core, stable upper core, and insulating mantle. An α² dynamo operates in the convecting outer core and a horizontal shear flow is imposed in the stable layer. Exact dynamo solutions are obtained for a range of parameters, including different conductivities for the stable layer and inner core. This allows us to connect our solutions with known, simpler solutions of a single-layer α² dynamo, and thereby assess the effects of the extra layers. We confirm earlier results that a stable, static layer can enhance dynamo action. We find that shear flows produce dynamo wave solutions with a different spatial structure from the steady α² dynamos solutions. The stable layer controls the behavior of the dynamo system through the interface conditions, providing a new means whereby lateral variations on the boundary can influence the geomagnetic field.     

Inner-core conductivity in numerical dynamo simulations

Speculation about its possible super-rotation has drawn the attention of many geophysical researchers to the Earth’s inner core. An issue of special interest for geodynamo modelling is the influence of the inner-core conductivity. It has been suggested that the finite magnetic diffusivity of the inner core prevents more frequent reversals of the Earth’s magnetic field. We explore the possible influence of the inner-core conductivity by comparing convection-driven 3D dynamo simulations with insulating or conducting inner cores (CIC) at various parameters. The influence on the field structure in the outer core is only marginal. The time behaviour of dipole-dominated non-reversing dynamos is also little affected. Concerning reversing dynamos, the inner-core conductivity reduces the number of short dipole-polarity intervals with a typical length of a few thousand years. Reversals are always correlated with low dipole strength and these short intervals are found in periods where the dipole moment stays low. Polarity intervals longer than about 10,000 years, where the dipole moment has time recover in strength, are equally likely in insulating and CIC models. Since these latter intervals are of more geophysical relevance, we conclude that the influence of the inner-core conductivity on Earth-like reversal sequences is insignificant for the dynamo model employed here.     

Fast dynamo action in the Ponomarenko dynamo

Abstract

An analysis of small-scale magnetic fields shows that the Ponomarenko dynamo is a fast dynamo; the maximum growth rate remains of order unity in the limit of large magnetic Reynolds number. Magnetic fields are regenerated by a “stretch-diffuse” mechanism. General smooth axisymmetric velocity fields are also analysed; these give slow dynamo action by the same mechanism.     

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.     

The flux ejection dynamo effect

Abstract

The mean-field effects of cyclonic convection become increasingly complex when the cyclonic rotation exceeds ½-π. Net helicity is not required, with negative turbulent diffusion, for instance, appearing in mirror symmetric turbulence. This paper points out a new dynamo effect arising in convective cells with strong asymmetry in the rotation of updrafts as against downdrafts. The creation of new magnetic flux arises from the ejection of reserve flux through the open boundary of the dynamo region. It is unlike the familiar α-effect in that individual components of the field may be amplified independently. Several formal examples are provided to illustrate the effect. Occurrence in nature depends upon the existence of fluid rotations of the order of π in the convective updrafts. The flux ejection dynamo may possibly contribute to the generation of field in the convective core of Earth and in the convective zone of the sun and other stars.     

A dynamo model with double diffusive convection for Mercury's core

A recent dynamo model for Mercury assumes that the upper part of the planet's fluid core is thermally stably stratified because the temperature gradient at the core–mantle boundary is subadiabatic. Vigorous convection driven by a superadiabatic temperature gradient at the boundary of a growing solid inner core and by the associated release of light constituents takes place in a deep sub-layer and powers a dynamo. These models have been successful at explaining the observed weak global magnetic field at Mercury's surface. They have been based on the concept of codensity, which combines thermal and compositional sources of buoyancy into a single variable by assuming the same diffusivity for both components. Actual diffusivities in planetary cores differ by a large factor. To overcome the limitation of the codensity model, we solve two separate transport equations with different diffusivities in a double diffusive dynamo model for Mercury. When temperature and composition contribute comparable amounts to the buoyancy force, we find significant differences to the codensity model. In the double diffusive case convection penetrates the upper layer with a net stable density stratification in the form of finger convection. Compared to the codensity model, this enhances the poloidal magnetic field in the nominally stable layer and outside the core, where it becomes too strong compared to observation. Intense azimuthal flow in the stable layer generates a strong axisymmetric toroidal field. We find in double diffusive models a surface magnetic field of the observed strength when compositional buoyancy plays an inferior role for driving the dynamo, which is the case when the sulphur concentration in Mercury's core is only a fraction of a percent.     

The dynamo basis of solar cycle precursor schemes

We investigate the dynamo underpinning of solar cycle precursor schemes based on direct or indirect measures of the solar surface magnetic field. We do so for various types of mean-field-like kinematic axisymmetric dynamo models, where amplitude fluctuations are driven by zero-mean stochastic forcing of the dynamo number controlling the strength of the poloidal source term. In all stochastically forced models considered, the surface poloidal magnetic field is found to have precursor value only if it feeds back into the dynamo loop, which suggests that accurate determination of the magnetic flux budget of the solar polar fields may hold the key to dynamo model-based cycle forecasting.     

A shell model for a large-scale turbulent dynamo

This article addresses the interesting and important problem of large-scale magnetic field generation in turbulent flows, using a self-consistent dynamo model recently developed. The main idea of this model is to consider the induction equation for the large-scale magnetic field, integrated consistently with the turbulent dynamics at smaller scales described by a magnetohydrodynamic shell model. The questions of dynamo action threshold, magnetic field saturation, magnetic field reversals, nature of the dynamo transition and the changes of small-scale turbulence as a consequence of the dynamo onset are discussed. In particular, the stability curve obtained by the model integration is shown in a very wide range of values of the magnetic Prandtl number not yet accessible by direct numerical simulation but more realistic for natural dynamos. Moreover, from our analysis it is shown that the large-scale dynamo transition displays a hysteretic behaviour and therefore a subcritical nature. The model successfully reproduces magnetic polarity reversals, showing the capability to generate persistence times which are increasing for decreasing magnetic diffusivity. Moreover, when the system reaches a statistically stationary dynamo state, where the large-scale magnetic field can abruptly reverse its polarity (magnetic reversal state) or not, keeping the same polarity (steady state), it shows an unmistakable tendency towards the energy equipartition for the turbulence at small scale.     

Periodic and aperiodic dynamo waves

Abstract

In order to show that aperiodic magnetic cycles, with Maunder minima, can occur naturally in nonlinear hydromagnetic dynamos, we have investigated a simple nonlinear model of an oscillatory stellar dynamo. The parametrized mean field equations in plane geometry have a Hopf bifurcation when the dynamo number D=1, leading to Parker's dynamo waves. Including the nonlinear interaction between the magnetic field and the velocity shear results in a system of seven coupled nonlinear differential equations. For D>1 there is an exact nonlinear solution, corresponding to periodic dynamo waves. In the regime described by a fifth order system of equations this solution remains stable for all D and the velocity shear is progressively reduced by the Lorentz force. In a regime described by a sixth order system, the solution becomes unstable and successive transitions lead to chaotic behaviour. Oscillations are aperiodic and modulated to give episodes of reduced activity.     

Nonlinearities in dynamo

By the example of the dynamo model in the rotating plane layer heated from below, the effects are examined that lead to the stabilization of an exponentially growing magnetic field in the magnetostrophic convection in passing from the kinematic dynamo mode to the nonlinear mode. The estimates of the energy redistribution in the spectrum are given, and the mechanisms of suppression of helicity are presented. Equalization of the field of velocity and the magnetic field is analyzed. The modes examined are close to those utilized in the up-to-date models of the planetary dynamo in the cores of planets.     

Implementation of a high-order combined compact difference scheme in problems of thermally driven convection and dynamo in rotating spherical shells

We present an improved solution method for modeling thermally driven convection and dynamo in a rotating spherical shell. In this method, we introduce a high-order three-point combined compact difference scheme (CCDS) on non-uniform grid points in radius, while spherical harmonic expansion is conventionally performed in the angular direction. The governing equations in the spectral form are time-stepped together with the implicit CCDS up to the second derivative. To improve stability of the scheme, a boundary closure scheme is developed on non-uniform mesh. Numerical comparison with a published benchmark solution at moderate Ekman and Rayleigh numbers demonstrates that accuracy and convergence of the CCDS is fairly good and superior to the existing finite difference scheme using more stencil. With this scheme, we could more accurately solve problems of convection and also dynamo action in planetary core with less grid points.     

The initial temporal evolution of a feedback dynamo for Mercury

Various possibilities are currently under discussion to explain the observed weakness of the intrinsic magnetic field of planet Mercury. One of the possible dynamo scenarios is a dynamo with feedback from the magnetosphere. Due to its weak magnetic field, Mercury exhibits a small magnetosphere whose subsolar magnetopause distance is only about 1.7 Hermean radii. We consider the magnetic field due to magnetopause currents in the dynamo region. Since the external field of magnetospheric origin is antiparallel to the dipole component of the dynamo field, a negative feedback results. For an αΩ-dynamo, two stationary solutions of such a feedback dynamo emerge: one with a weak and the other with a strong magnetic field. The question, however, is how these solutions can be realized. To address this problem, we discuss various scenarios for a simple dynamo model and the conditions under which a steady weak magnetic field can be reached. We find that the feedback mechanism quenches the overall field to a low value of about 100–150 nT if the dynamo is not driven too strongly.     

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.     

Properties of nonlinear dynamo waves

Abstract

Dynamo theory offers the most promising explanation of the generation of the sun's magnetic cycle. Mean field electrodynamics has provided the platform for linear and nonlinear models of solar dynamos. However the nonlinearities included arc (necessarily) arbitrarily imposed in these models. This paper conducts a systematic survey of the role of nonlinearities in the dynamo process, by cousidering the behaviour of dynamo waves in the nonlinear regime. It is demonstrated that only by considering realistic nonlinearities that are non-local in space and time can modulation of the basic dynamo wave be achieved. Moreover this modulation is greatest when there is a large separation of timescales provided by including a low magnetic Prandtl number in the equation for the velocity perturbations.     

Kinematic dynamo induced by helical waves

We investigate numerically the kinematic dynamo induced by the superposition of two helical waves in a periodic box as a simplified model to understand the dynamo action in astronomical bodies. The effects of magnetic Reynolds number, wavenumber and wave frequency on the dynamo action are studied. It is found that this helical-wave dynamo is a slow dynamo. There exists an optimal wavenumber for the dynamo growth rate. A lower wave frequency facilitates the dynamo action and the oscillations of magnetic energy emerge at some particular wave frequencies.     

Self-consistent dynamo models driven by hydromagnetic instabilities

The dynamics of the Earth's core are dominated by a balance between Lorentz and Coriolis forces. Previous studies of possible (magnetostrophic) hydromagnetic instabilities in this regime have been confined to geophysically unrealistic flows and fields. In recent papers we have treated rather general fields and flows in a spherical geometry and in a computationally simple plane-layer model. These studies have highlighted the importance of differential rotation in determining the spatial structure of the instability. Here we have proceeded to use these results to construct a self-consistent dynamo model of the geomagnetic field. An iterative procedure is employed in which an α-effect is calculated from the form of the instability and is then used in a mean field dynamo model. The mean zonal field calculated there is then input back into the hydromagnetic stability problem and a new α-effect calculated. The whole procedure is repeated until the input and output zonal fields are the same to some tolerance.     

A nonlinear dynamo driven by rapidly rotating convection

In Kim et al. (Kim, E., Hughes, D.W. and Soward, A.M., “An investigation into high conductivity dynamo action driven by rotating convection”, Geophys. Astrophys. Fluid Dynam. 91, 303–332 ().) we investigated kinematic dynamo action driven by rapidly rotating convection in a cylindrical annulus. Here we extend this work to consider self-consistent nonlinear dynamo action in which the back-reaction of the Lorentz force on the flow is taken into account. In particular, we investigate, as a function of magnetic Prandtl number, the evolution of an initially weak magnetic field in two different types of convective flow – one chaotic and the other integrable. On saturation, the latter shows a systematic dependence on the magnetic Prandtl number whereas the former appears not to. In addition, we show how, in keeping with the findings of Cattaneo et al. (Cattaneo, F., Hughes, D.W. and Kim, E., “Suppression of chaos in a simplified nonlinear dynamo model”, Phys. Rev. Lett. 76, 2057–2060 ().), saturation of the growth of the magnetic field is brought about, for the originally chaotic flow, by a strong suppression of chaos.     

Nonlinear dynamics of a stellar dynamo in a spherical shell

Abstract

A simple mean-field model of a nonlinear stellar dynamo is considered, in which dynamo action is supposed to occur in a spherical shell, and where the only nonlinearity retained is the influence of the Lorentz forces on the zonal flow field. The equations are simplified by truncating in the radial direction, while full latitudinal dependence is retained. The resulting nonlinear p.d.e.'s in latitude and time are solved numerically, and it is found that while regular dynamo wave type solutions are stable when the dynamo number D is sufficiently close to its critical value, there is a wide variety of stable solutions at larger values of D. Furthermore, two different types of dynamo can coexist at the same parameter values. Implications for fields in late-type stars are discussed.     

Inverse dynamo problem in a cylinder

An inverse dynamo problem is presented in which we search for either kinematic dynamos which produce the same external magnetic fields or an invisible dynamo. The existence of flows which produce the same external magnetic fields is proved. However, we have not found general conditions necessary for such kind of dynamos. An “invisible dynamo” operates in an electrically conducting region surrounded by vacuum and generates a magnetic field trapped in the electrically conducting region so that no magnetic field exists in the vacuum. Invisible magnetic decay modes exist in cylinders, but no invisible growing field supported by the dynamo mechanism has been found.