Solar dynamo and geomagnetic activity

The correlation between geomagnetic activity and the sunspot number in the 11-year solar cycle exhibits long-term variations due to the varying time lag between the sunspot-related and non-sunspot related geomagnetic activity, and the varying relative amplitude of the respective geomagnetic activity peaks. As the sunspot-related and non-sunspot related geomagnetic activity peaks are caused by different solar agents, related to the solar toroidal and poloidal fields, respectively, we use their variations to derive the parameters of the solar dynamo transforming the poloidal field into toroidal field and back. We find that in the last 12 cycles the solar surface meridional circulation varied between 5 and 20 m/s (averaged over latitude and over the sunspot cycle), the deep circulation varied between 2.5 and 5.5 m/s, and the diffusivity in the whole of the convection zone was ～10⁸ m²/s. In the last 12 cycles solar dynamo has been operating in moderately diffusion dominated regime in the bulk of the convection zone. This means that a part of the poloidal field generated at the surface is advected by the meridional circulation all the way to the poles, down to the tachocline and equatorward to sunspot latitudes, while another part is diffused directly to the tachocline at midlatitudes, “short-circuiting” the meridional circulation. The sunspot maximum is the superposition of the two surges of toroidal field generated by these two parts of the poloidal field, which is the explanation of the double peaks and the Gnevyshev gap in sunspot maximum. Near the tachocline, dynamo has been operating in diffusion dominated regime in which diffusion is more important than advection, so with increasing speed of the deep circulation the time for diffusive decay of the poloidal field decreases, and more toroidal field is generated leading to a higher sunspot maximum. During the Maunder minimum the dynamo was operating in advection dominated regime near the tachocline, with the transition from diffusion dominated to advection dominated regime caused by a sharp drop in the surface meridional circulation which is in general the most important factor modulating the amplitude of the sunspot cycle.     

Solar dynamo theory: a new look at the origin of small-scale magnetic fields

Fausto Cattaneo and David W Hughes delve beneath the surface of the Sun with numerical models of turbulent convection.
Although magnetic dynamo action is traditionally associated with rotation, fast dynamo theory shows that chaotic flows, even without rotation, can act as efficient small-scale dynamos. Indeed, numerical simulations suggest that granular and supergranular convection may generate locally a substantial part of the field in the quiet photosphere.     

Toroidal flux oscillation as possible cause of geomagnetic excursions and reversals

It is proposed that convection driven dynamos operating in planetary cores could be oscillatory even when the oscillations are not directly noticeable from the outside. Examples of dynamo simulations are pointed out that exhibit oscillations in the structure of the azimuthally averaged toroidal magnetic flux while the mean poloidal field shows only variations in its amplitude. In the case of the geomagnetic field, global excursions may be associated with these oscillations. Long period dynamo simulations indicate that the oscillations may cause reversals once in a while. No special attempt has been made to use most realistic parameter values. Nevertheless some similarities between the simulations and the paleomagnetic record can be pointed out.     

Convective dynamos with intermediate and strong fields

Abstract

The weak-field Benard-type dynamo treated by Soward is considered here at higher levels of the induced magnetic field. Two sources of instability are found to occur in the intermediate field regime M ～ T ^1/12, where M and T are the Hartmann and Taylor numbers. On the time scale of magnetic diffusion, solutions may blow up in finite time owing to destabilization of the convection by the magnetic field. On a faster time scale a dynamic instability related to MAC-wave instability can also occur. It is therefore concluded that the asymptotic structure of this dynamo is unstable to virtual increases in the magnetic field energy.

In an attempt to model stabilization of the dynamo in a strong-field regime we consider two approximations. In the first, a truncated expansion in three-dimensional plane waves is studied numerically. A second approach utilizes an ad hoc set of ordinary differential equations which contains many of the features of convection dynamos at all field energies. Both of these models exhibit temporal intermittency of the dynamo effect.     

Distribution of magnetic energy in αΩ-dynamos,I: The method

Abstract

In this paper a method for solving the equation for the mean magnetic energy <BB> of a solar type dynamo with an axisymmetric convection zone geometry is developed and the main features of the method are described. This method is referred to as the finite magnetic energy method since it is based on the idea that the real magnetic field B of the dynamo remains finite only if <BB> remains finite. Ensemble averaging is used, which implies that fields of all spatial scales are included, small-scale as well as large-scale fields. The method yields an energy balance for the mean energy density ε ≡ B ²/8π of the dynamo, from which the relative energy production rates by the different dynamo processes can be inferred. An estimate for the r.m.s. field strength at the surface and at the base of the convection zone can be found by comparing the magnetic energy density and the outgoing flux at the surface with the observed values. We neglect resistive effects and present arguments indicating that this is a fair assumption for the solar convection zone. The model considerations and examples presented indicate that (1) the energy loss at the solar surface is almost instantaneous; (2) the convection in the convection zone takes place in the form of giant cells; (3) the r.m.s. field strength at the base of the solar convection zone is no more than a few hundred gauss; (4) the turbulent diffusion coefficient within the bulk of the convection zone is about 10¹⁴cm²s^?1, which is an order of magnitude larger than usually adopted in solar mean field models.     

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.     

Helicity–vorticity turbulent pumping of magnetic fields in the solar convection zone

We study the effect of turbulent drift of a large-scale magnetic field that results from the interaction of helical convective motions and differential rotation in the solar convection zone. The principal direction of the drift corresponds to the direction of the large-scale vorticity vector. Thus, the effect produces a latitudinal transport of the large-scale magnetic field in the convective zone wherever the angular velocity has a strong radial gradient. The direction of the drift depends on the sign of helicity and it is defined by the Parker–Yoshimura rule. The analytic calculations are done within the framework of mean-field magnetohydrodynamics using the minimal τ-approximation. We estimate the magnitude of the drift velocity and find that it can be a few m/s near the base of the solar convection zone. The implications of this effect for the solar dynamo are illustrated on the basis of an axisymmetric mean-field dynamo model with a subsurface shear layer. The model shows that near the bottom of the convection zone the helicity–vorticity pumping results mostly from the kinetic helicity contributions. We find that the magnetic helicity contributions to the pumping effect are dominant at the subsurface shear layer. There the magnitude of the drift velocity is found to be a few cm/s. We find that the helicity–vorticity pumping effect can have an influence on the features of the sunspot time–latitude diagram, producing a fast drift of the sunspot activity maximum at the rise phase of the cycle and a slow drift at the decay phase of the cycle.     

Numerical simulations of MHD dynamos

The generation of magnetic fields in space plasmas and in astrophysics is usually described within the framework of magnetohydrodynamics. Turbulent helical flows produce magnetic fields very efficiently, with correlation length scales larger than those characterizing the flow. Within the context of the solar magnetic cycle, a turbulent dynamo is responsible for the so-called alpha effect, while the Omega effect is associated to the differential rotation of the Sun.We present direct numerical simulations of turbulent magnetohydrodynamic dynamos including two-fluid effects such as the Hall current. More specifically, we study the evolution of an initially weak and small-scale magnetic field in a system maintained in a stationary regime of hydrodynamic turbulence, and explore the conditions for exponential growth of the magnetic energy. In all the cases considered, we find that the dynamo saturates at the equipartition level between kinetic and magnetic energy, and the total energy reaches a Kolmogorov power spectrum.     

The weak taylor state in an αω-dynamo

Abstract

A spherical αω-dynamo is studied for small values of the viscous coupling parameter ε ～ v^1/2, paying attention particularly to large dynamo numbers. The present study is a follow-up of the work by Hollerbach et al. (1992) with their choice of α-effect and Archimedean wind including also the constraint of magnetic field symmetry (or antisymmetry) due to equatorial plane. The magnetic field scaled by ε^1/2 is independent of ε in the solutions for dynamo numbers smaller than a certain value of D _b (the Ekman state) which are represented by dynamo waves running from pole to equator or vice-versa. However, for dynamo numbers larger than D _b the solution bifurcates and subsequently becomes dependent on ε. The bifurcation is a consequence of a crucial role of the meridional convection in the mechanism of magnetic field generation. Calculations suggest that the bifurcation appears near dynamo number about 33500 and the solutions for larger dynamo numbers and ε = 0 become unstable and fail, while the solutions for small but non-zero ε are characterized by cylindrical layers of local maximum of magnetic field and sharp changes of geostrophic velocity. Our theoretical analysis allows us to conclude that our solution does not take the form of the usual Taylor state, where the Taylor constraint should be satisfied due to the special structure of magnetic field. We rather obtained the solution in the form of a “weak” Taylor state, where the Taylor constraint is satisfied partly due to the amplitude of the magnetic field and partly due to its structure. Calculations suggest that the roles of amplitude and structure are roughly fifty-fifty in our “weak” Taylor state solution and thus they can be called a Semi-Taylor state. Simple estimates show that also Ekman state solutions can be applicable in the geodynamo context.     

Magnetic tracers,a probe of the solar convective layers

Abstract

A meridional circulation of sunspots has been measured through the digital analysis of the Meudon spectroheliograms from 1978 to 1983. Old and young sunspots follow a zonal meridional circulation, in several bands of latitude, in which two adjacent bands have opposite motions. This meridional circulation pattern is time-dependent. Using the H _α filaments as magnetic field tracers, a large-scale magnetic pattern has been found that was also obtained independently by direct measurement of the magnetic field (Hoeksema, 1988).

The coincidence of a large-scale magnetic pattern with a zonal meridional circulation suggests the existence of azimuthal rolls below the surface, and these azimuthal rolls can explain a number of properties of the solar cycle. New rolls occur with increasing proximity to the Equator, thereby indicating the direction of propagation of the dynamo wave. The occurrence of rolls is very favorable to the emergence of the magnetic regions. The rolls also influence the magnetic complexity of the active regions. They modulate the surface rotation through the Coriolis force, which accelerates or decelerates the fluid particles. They therefore offer a plausible explanation of the torsional oscillation pattern.

There are a number of problems raised by such an unexpected circulation pattern: for example, the coexistence of axisymmeric rolls with hypothetical giant cells, the location of the dynamo source below or within the convective zone, and the coupling of the radiative interior and the convective layers. To resolve these important issues, continuous observational studies are needed of the manifestation of solar activity, as well as of radius and luminosity variations. So, we have aimed our paper at an audience of theoreticians in the hope that they take up the challenges we describe.     

Magnetic field generation by convective flows in a plane layer: the dependence on the Prandtl numbers

Investigation of magnetic field generation by convective flows is carried out for three values of kinematic Prandtl number: P = 0.3, 1 and 6.8. We consider Rayleigh–Bénard convection in Boussinesq approximation assuming stress-free boundary conditions on horizontal boundaries and periodicity with the same period in the x and y directions. Convective attractors are modelled for increasing Rayleigh numbers for each value of the kinematic Prandtl number. Linear and non-linear dynamo action of these attractors is studied for magnetic Prandtl numbers P _m ≤ 100. Flows, which can act as magnetic dynamos, have been found for all the three considered values of P, if the Rayleigh number R is large enough. The minimal R, for which of magnetic field generation occurs, increases with P. The minimum (over R) of critical P_m for magnetic field generation in the kinematic regime is admitted for P = 0.3. Thus, our study indicates that smaller values of P are beneficial for magnetic field generation.     

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.     

Thermal convection and the dynamo during rapid rotation

A dynamo model for an incompressible liquid in a rectangular rapidly rotating box is considered. The model is based on a pseudospectral method adapted for multiprocessor technology. The effect of rotation in Boussinesq models on some of the basic characteristics of flows is studied: the spatial form of fields, helicity generation, and spectra of fields. A mechanism stabilizing the magnetic field rise associated with the transition of the system to a nonlinear regime is analyzed separately. The results obtained in the paper provide deeper insights into the processes of magnetoconvection in the cores of planets.     

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.     

Evolution of the dipole geomagnetic field. Observations and models

The works on paleomagnetic observations of the dipole geomagnetic field, its variations, and reversals in the last 3.5 billion years have been reviewed. It was noted that characteristic field variations are related to the evolution of the convection processes in the liquid core due to the effect of magnetic convection and solid core growth. Works on the geochemistry and energy budget of the Earth’s core, the effect of the solid core on convection and the generation of the magnetic field, dynamo models are also considered. We consider how core growth affects the magnetic dipole generation and variations, as well as the possibility of magnetic field generation up to the appearance of the solid core. We also pay attention to the fact that not only the magnetic field but also its configuration and time variations, which are caused by the convection evolution in the core on geological timescales, are important factors for the biosphere.     

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.     

Electrostatic radiation of plasma in the upper ionosphere in the inhomogeneous geomagnetic field

It has been indicated that the spectrum of electrostatic waves in the ionospheric plasma depends on the geophysical conditions and solar wind parameters. The wave field measurements in the frequency band 0.1–10 MHz in the topside ionosphere were used to analyze the electrostatic instabilities of the plasma electron content (the APEX satellite experiment). A change of the sign of one magnetic field component at the geomagnetic equator can result in the formation of the large-scale irregular plasma structure with a decay of the natural electrostatic oscillations and vortices in unstable plasma. The plasma particle polarization drift from the region of decay of electrostatic oscillations and vortices can cause large plasma density and temperature gradients across the geomagnetic field. New vortices can originate at these gradients. This mechanism of plasma vortex formation and decay can be important for mass and energy convection in the topside ionosphere.     

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.     

From stable dipolar towards reversing numerical dynamos

We are using a three-dimensional convection-driven numerical dynamo model without hyperdiffusivity to study the characteristic structure and time variability of the magnetic field in dependence of the Rayleigh number (Ra) for values up to 40 times supercritical. We also compare a variety of ways to drive the convection and basically find two dynamo regimes. At low Ra, the magnetic field at the surface of the model is dominated by the non-reversing axial dipole component. At high Ra, the dipole part becomes small in comparison to higher multipole components. At transitional values of Ra, the dynamo vacillates between the dipole-dominated and the multipolar regime, which includes excursions and reversals of the dipole axis. We discuss, in particular, one model of chemically driven convection, where for a suitable value of Ra, the mean dipole moment and the temporal evolution of the magnetic field resemble the known properties of the Earth’s field from paleomagnetic data.     

Baroclinically-driven flows and dynamo action in rotating spherical fluid shells

The dynamics of stably stratified stellar radiative zones is of considerable interest due to the availability of increasingly detailed observations of Solar and stellar interiors. This article reports the first non-axisymmetric and time-dependent simulations of flows of anelastic fluids driven by baroclinic torques in stably stratified rotating spherical shells – a system serving as an elemental model of a stellar radiative zone. With increasing baroclinicity a sequence of bifurcations from simpler to more complex flows is found in which some of the available symmetries of the problem are broken subsequently. The poloidal component of the flow grows relative to the dominant toroidal component with increasing baroclinicity. The possibility of magnetic field generation thus arises and this paper proceeds to provide some indications for self-sustained dynamo action in baroclinically-driven flows. We speculate that magnetic fields in stably stratified stellar interiors are thus not necessarily of fossil origin as it is often assumed.