Large-scale convective dynamos in a stratified rotating plane layer

We study the effect of stratification on large-scale dynamo action in convecting fluids in the presence of background rotation. The fluid is confined between two horizontal planes and both boundaries are impermeable, stress-free and perfectly conducting. An asymptotic analysis is performed in the limit of rapid rotation (τ???1 where τ is the Taylor number). We analyse asymptotic magnetic dynamo solutions in rapidly rotating systems generalising the results of Soward [A convection-driven dynamo I. The weak field case. Philos. Trans. R. Soc. Lond. A 1974, 275, 611–651] to include the effects of compressibility. We find that in general the presence of stratification delays the efficiency of large-scale dynamo action in this regime, leading to a reduction of the onset of dynamo action and in the nonlinear regime a diminution of the large-scale magnetic energy for flows with the same kinetic energy.     

Models for convectively driven hydromagnetic dynamos

Abstract

Models of a convectively driven hydromagnetic dynamo are constructed using a truncated modal expansion. The resulting nonlinear partial differential equations are integrated numerically. The results confirm that rotation is a necessary condition for effective dynamo action, and suggest that equipartition of kinetic and magnetic energies is qualitatively valid, and that toroidal field energies can be much larger than poloidal.     

The rotational quenching of the rotation-induced kinetic alpha-effect

Abstract

A theory of the non-diffusive anisotropic kinetic alpha-effect (“Γ-effect”) for densitystratified rotating turbulent fluids is developed. No limitations on the rotation rate are imposed and the fully nonlinear dependence of the Γ-effect on the angular velocity is studied. When the Coriolis number, ω? = 2τ ω, is small the dimensionless “dynamo number”, Cτ, characterising the power of the Γ-effect, grows with ω?. The dependence, however, reaches a maximum for ω? ～ 2. For still higher rotation rates C_Λ decreases as 1/ω?. In opposition, the corresponding number, C_x, of the hydromagnetic α² -dynamo problems remains finite for very large ω?. Hence, for fast rotation the hydrodynamic Γ-effect is small while the hydromagnetic α-effect remains large. In consequence, the large-scale magnetic and velocity structures are expected to be generated with roughly equal power in slowly rotating objects. In the rapid rotators, however, generation of the large-scale flows is problematic.     

Maximally-efficient-generation approach in the dynamo theory

Abstract

We propose a method of derivation of global asymptotic solutions of the hydromagnetic dynamo problem at large magnetic Reynolds number. The procedure reduces to matching the local asymptotic forms for the magnetic field generated near individual extrema of generation strength. The basis of the proposed method, named here the Maximally-Efficient-Generation Approach (MEGA), is the assertion that properties of global asymptotic solutions of the kinematic dynamo are determined by the distribution of the generation strength near its leading extrema and by the number and distribution of the extrema.

The general method is illustrated by the global asymptotic solution of the α²-dynamo problem in a slab. The nature of oscillatory solutions revealed earlier in numerical simulations and the reasons for the dominance of even magnetic modes in slab geometry are clarified.

Applicability of the asymptotic solutions at moderate values of the asymptotic parameter is also discussed. We confirm this applicability using comparisons with complementary asymptotic expansions and numerical simulations. In particular, this justifies application of the MEGA solutions to estimation of the generation threshold.     

On the modified Taylor constraint

A dynamo driven by motions unaffected by viscous forces is termed magnetostrophic. Although such a model might describe magnetic field generation in Earth’s core well, a magnetostrophic dynamo has not yet been found even though Taylor [Proc. R. Soc. Lond. A 1963, 274, 274–283] devised an apparently viable method of finding one. His method for determining the fluid velocity from the magnetic field and the energy source involved only the evaluation of integrals along lines parallel to the Earth’s axis of rotation and the solution of a second-order ordinary differential equation. It is demonstrated below that an approximate solution of this equation for a broad family of magnetic fields is immediate. Furthermore inertia, which was neglected in Taylor’s theory, is restored here, so that the modified theory includes torsional waves, whose existence in the Earth’s core has been inferred from observations of the length of day. Their theory is reconsidered.     

Note on the scalar dynamo model

Abstract

Bayly (1993) introduced and investigated the equation (? _t + v·▽-η ▽²)S = RS as a scalar analogue of the magnetic induction equation. Here, S(r, t) is a scalar function and the flow field v(r, t) and “stretching” function R(r, t) are given independently. This equation is much easier to handle than the corresponding vector equation and, although not of much relevance to the (vector) kinematic dynamo problem, it helps to study some features of the fast dynamo problem. In this note the scalar equation is considered for linear flow and a harmonic potential as stretching function. The steady equation separates into one-dimensional equations, which can be completely solved and therefore allow one to monitor the behaviour of the spectrum in the limit of vanishing diffusivity. For more general homogeneous flows a scaling argument is given which ensures fast dynamo action for certain powers of the harmonic potential. Our results stress the singular behaviour of eigenfunctions in the limit of vanishing diffusivity and the importance of stagnation points in the flow for fast dynamo action.     

Analytic base of geodynamo-like scaling laws in the planets,geomagnetic periodicities and inversions

Scaling laws for hydromagnetic dynamo in planets initially express the characteristic strength of the magnetic field through the primary values, such as the size of the conductive core of the planet, the angular rotation rate, electrical conductivity and energy flows. Most of the earlier proposed scaling laws based only on observations and assumptions about force balances. Recent and my new approaches to fully take into account the energy and induction balance has additionally expressed here in terms of primary values such important characteristics as forces, magnitudes, energies, scales and orientations of hydromagnetic fields. The direct numerical simulation of the hydromagnetic dynamo and modeling ability in a fairly wide range of parameters for the first time allowed direct test such laws. The obtained numerical geodynamo-like results for the Earth, Jupiter and partially Saturn postulated previously not identified analytically simplest law that predicts the field strength is only depended on the specific energy density of convection and the size of the dynamo area. This simplest and already widely used law was original way analytically grounded here along with other previously known and new laws. This analytic identifies the physics determining geomagnetic periodicities for jerk, secular variations and inversions. Mean period between the inversions is found to be roughly proportional to the intensity of the geomagnetic field that is confirmed by some paleomagnetic researches. Possible dynamos in Mercury, Ganymede, Uranus and Neptune are also discussed.     

Hydromagnetic dynamo model with spiral convection

Summary The present paper deals with a hydromagnetic dynamo model of the generation mechanism of the Earth's magnetic field. An attempt has been made at selecting a flow-velocity field in the Earth's core which would satisfy the condition 0 for regenerating the field according to [2], and which would yield a velocity field pattern on the core surface as given in the papers by Kahle et al. [9]. These conditions are satisfied by the velocityv=V ₁+U ₂ ^c –V ₂ ^c and, geometrically, this velocity field is represented in space by a spiral convective motion. On the core surface two downflows and two upflows with the corresponding rotating cells may then be found. Only the axisymmetric harmonic component regeneration of the magnetic field has been considered. Adequate regeneration equations have been obtained by means of Braginski's method of quantity estimates in order of magnitude.     

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.     

Convective Instability of a Mushy Layer - I: Uniform Permeability

Mushy layers arise and are significant in a number of geophysical contexts, including freezing of sea ice, solidification of magma chambers and inner-core solidification. A mushy layer is a region of solid and liquid in phase equilibrium which commonly forms between the liquid and solid regions of a solidifying system composed of two or more constituents. We consider the convective instability of a plane mushy layer which advances steadily upwards as heat is withdrawn at a uniform rate from the bottom of a eutectic binary alloy. The solid which forms is assumed to be composed entirely of the denser constituent, making the residual liquid within the mush compositionally buoyant and thus prone to convective motion. In this article we focus on the large-scale mush mode of instability, arguing that the 'boundary-layer' mode is not amenable to the standard stability analysis, because convective motions occur on that scale for any non-zero value of the Rayleigh number. We quantify the minimum critical Rayleigh number and determine the structure of the convective modes of motion within the mush and the associated deflections of the mush-melt and mush-solid boundaries. This study of convective perturbations differs from previous analyses in two ways; the inhibition of motion and deformation of the mush-melt interface by the stable stratification of the overlying melt is properly quantified and deformation of the mush-solid interface is permitted and quantified. We find that the mush-melt interface is almost unaffected by convection while significant deformation of the mush-solid interface occurs. We show that each of these effects causes significant (unit-order) changes in the predicted critical Rayleigh number. The marginal modes depend on three dimensionless parameters: a scaled eutectic temperature, τ _e (which characterizes the eutectic temperature relative to the depression of the liquidus), a scaled superheat, τ_∞ (which measures the amount by which the temperature of the incoming melt exceeds the liquidus temperature) and the Stefan number, S (which measures the latent heat of crystallization). To survey parameter space, we focus on seven cases, a standard case having S = τ_∞ = τ _e = 1, and six others in which one of the parameters is either large or small compared with unity: a nearly pure case (τ _e = 100; having little of the light constituent), the large superheat limit (τ_∞→ ∞), a case of large latent heat (S = 100), the near eutectic limit (τ _e → 0), a case of small superheat (τ_∞ = 0.01) and the case of zero latent heat (S = 0). The critical Rayleigh number and the associated wavelength of the convection pattern are determined in each case. The eigenvector for each case is presented in terms of the streamlines and the isolines of the perturbation temperature and solid fraction.     

Analysis of the Characteristics of Low-Latitude GPS Amplitude Scintillation Measured During Solar Maximum Conditions and Implications for Receiver Performance

Ionospheric scintillations are fluctuations in the phase and/or amplitude of trans-ionospheric radio signals caused by electron density irregularities in the ionosphere. A better understanding of the scintillation pattern is important to make a better assessment of GPS receiver performance, for instance. Additionally, scintillation can be used as a tool for remote sensing of ionospheric irregularities. Therefore, the study of ionospheric scintillation has both scientific as well as technological implications. In the past few years, there has been a significant advance in the methods for analysis of scintillation and in our understanding of the impact of scintillation on GPS receiver performance. In this work, we revisit some of the existing methods of analysis of scintillation, propose an alternative approach, and apply these techniques in a comprehensive study of the characteristics of amplitude scintillation. This comprehensive study made use of 32?days of high-rate (50?Hz) measurements made by a GPS-based scintillation monitor located in S?o José dos Campos, Brazil (23.2°S, 45.9°W, ?17.5° dip latitude) near the Equatorial Anomaly during high solar flux conditions. The variability of the decorrelation time (τ₀) of scintillation patterns is presented as a function of scintillation severity index (S ₄). We found that the values of τ₀ tend to decrease with the increase of S ₄, confirming the results of previous studies. In addition, we found that, at least for the measurements made during this campaign, averaged values of τ₀ (for fixed S ₄ index values) did not vary much as a function of local time. Our results also indicate a significant impact of τ₀ in the GPS carrier loop performance for S ₄?≥?0.7. An alternative way to compute the probability of cycle slip that takes into account the fading duration time is also presented. The results of this approach show a 38% probability of cycle slips during strong scintillation scenarios (S ₄ close to 1 and τ₀ near 0.2?s). Finally, we present results of an analysis of the individual amplitude fades observed in our set of measurements. This analysis suggests that users operating GPS receivers with C/N ₀ thresholds around 30?dB-Hz and above can be affected significantly by low-latitude scintillation.     

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.     

The DRESDYN project: liquid metal experiments on dynamo action and magnetorotational instability

ABSTRACT

Magnetic fields of planets, stars and galaxies are generated by self-excitation in moving electrically conducting fluids. Once produced, magnetic fields can play an active role in cosmic structure formation by destabilising rotational flows that would be otherwise hydrodynamically stable. For a long time, both hydromagnetic dynamo action as well as magnetically triggered flow instabilities had been the subject of purely theoretical research. Meanwhile, however, the dynamo effect has been observed in large-scale liquid sodium experiments in Riga, Karlsruhe and Cadarache. In this paper, we summarise the results of liquid metal experiments devoted to the dynamo effect and various magnetic instabilities such as the helical and the azimuthal magnetorotational instability and the Tayler instability. We discuss in detail our plans for a precession-driven dynamo experiment and a large-scale Tayler–Couette experiment using liquid sodium, and on the prospects to observe magnetically triggered instabilities of flows with positive shear.     

Earthquake response spectra taking account of number of response cycles

Proposed is a new definition of earthquake response spectra, which takes account of the number of response cycles N. The Nth largest amplitude of absolute acceleration response of a linear oscilator with natural period T and damping ratio h, which is subjected to ground motion at its base, is defined as S_A(T, h, N). By defining a reduction factor η(T, h, N) as S_A(T, h, N)/S_A(T, h, 1), characteristics of η(T, h, N) were investigated based on 394 components of strong motion records obtained in Japan. Two practical empirical formulae to assess the reduction factor η(T, h, N) are proposed.     

How strong is the magnetic field in the Earth's liquid core?

A crucial step in the investigation of the energetics of motions in the Earth's core and the generation of the geomagnetic field by the hydromagnetic dynamo process is the estimation of the average strength $\text{B}$ of the magnetic field B = B_p + B_T in the core. Owing to the probability that the toroidal field B_T in the core, which has no radial component, is a good deal stronger than the poloidal field B_p, direct downward extrapolation of the surface field to the core-mantle interface gives no more than an extreme lower limit to $\text{B}$. This paper outlines the indirect methods by which $\text{B}$ can be estimated, arguing that $\text{B}$ is probably about 10^?2 T (100 Γ) but might be as low as 10^?3 T (10 Γ) or as high as 5 × 10^?2 T (500 Γ).     

The distribution of 210Pb and 210Po in the surface waters of the Pacific Ocean

By modelling the observed distribution of²¹⁰Pb and²¹⁰Po in surface waters of the Pacific, residence times relative to particulate removal are determined. For the center of the North Pacific gyre these are τ_Po = 0.6years andτ_Pb = 1.7years. The surface ocean τ_Pb is determined by particulate transport rather than plankton settling. The fact that it is about two orders of magnitude smaller than τ_Pb for the deep ocean implies a sharp change in the adsorptive quality of particles during descent through the water column.     

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.     

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.     

Magnetic instabilities in rapidly rotating spherical geometries I. from cylinders to spheres

Abstract

The solution of the full nonlinear hydromagnetic dynamo problem is a major numerical undertaking. While efforts continue, supplementary studies into various aspects of the dynamo process can greatly improve our understanding of the mechanisms involved. In the present study, the linear stability of an electrically conducting fluid in a rigid, electrically insulating spherical container in the presence of a toroidal magnetic field B_o(r,θ)lø and toroidal velocity field U_o(r,θ)lø, [where (r,θ,ø) are polar coordinates] is investigated. The system, a model for the Earth's fluid core, is rapidly rotating, the magnetostrophic approximation is used and thermal effects are excluded. Earlier studies have adopted a cylindrical geometry in order to simplify the numerical analysis. Although the cylindrical geometry retains the fundamental physics, a spherical geometry is a more appropriate model for the Earth. Here, we use the results which have been found for cylindrical systems as guidelines for the more realistic spherical case. This is achieved by restricting attention to basic states depending only on the distance from the rotation axis and by concentrating on the field gradient instability. We then find that our calculations for the sphere are in very good qualitative agreement both with a local analysis and with the predictions from the results of the cylindrical geometry. We have thus established the existence of field gradient modes in a realistic (spherical) model and found a sound basis for the study of various other, more complicated, classes of magnetically driven instabilities which will be comprehensively investigated in future work.     

Kinematic dynamo action in a sphere: Effects of periodic time-dependent flows on solutions with axial dipole symmetry

Choosing a simple class of flows, with characteristics that may be present in the Earth's core, we study the ability to generate a magnetic field when the flow is permitted to oscillate periodically in time. The flow characteristics are parameterised by D, representing a differential rotation, M, a meridional circulation, and C, a component characterising convective rolls. The dynamo action of all solutions with fixed parameters (steady flows) is known from earlier studies. Dynamo action is sensitive to these flow parameters and fails spectacularly for much of the parameter space where magnetic flux is concentrated into small regions, leading to high diffusion. In addition, steady flows generate only steady or regularly reversing oscillatory fields and cannot therefore reproduce irregular geomagnetic-type reversal behaviour. Oscillations of the flow are introduced by varying the flow parameters in time, defining a closed orbit in the space ( D,?M ). When the frequency of the oscillation is small, the net growth rate of the magnetic field over one period approaches the average of the growth rates for steady flows along the orbit. At increased frequency time-dependence appears to smooth out flux concentrations, often enhancing dynamo action. Dynamo action can be impaired, however, when flux concentrations of opposite signs occur close together as smoothing destroys the flux by cancellation. It is possible to produce geomagnetic-type reversals by making the orbit stray into a region where the steady flows generate oscillatory fields. In this case, however, dynamo action was not found to be enhanced by the time-dependence. A novel approach is being taken to solve the time-dependent eigenvalue problem where, by combining Floquet theory with a matrix-free Krylov-subspace method, we can avoid large memory requirements for storing the matrix required by the standard approach.