Gary A. Glatzmaier

A three-dimensional convective dynamo solution with rotating and finitely conducting inner core and mantle

We present the first three-dimensional (3D), time-dependent, self-consistent numerical solution of the magneto- hydrodynamic (MHD) equations that describe thermal convection and magnetic field generation in a rapidly rotating spherical fluid shell with a solid conducting inner core. This solution, which serves as a crude analog for the geodynamo, is a self-sustaining supercritical dynamo that has maintained a magnetic field for three magnetic diffusion times, roughly 40 000 years. Fluid velocity in the outer core reaches a maximum of 0.4 cm s-1, and at times the magnetic field can be as large as 560 gauss. Magnetic energy is usually about 4000 times greater than the kinetic energy of the convection that maintains it. Viscous and magnetic coupling to both the inner core below and the mantle above cause time-dependent variations in their respective rotation rates; the inner core usually rotates faster than the mantle and decadal variations in the length of the day of the mantle are similar to those observed for the Earth. The pattern and amplitude of the radial magnetic field at the core-mantle boundary (CMB) and its secular variation are also similar to the Earths. The maximum amplitudes of the longitudinally averaged temperature gradient, shear flow, helicity, and magnetic field oscillate between the northern and southern hemispheres on a time scale of a few thousand years. However, only once in many attempts does the field succeed in reversing its polarity because the field in the inner core, which has the opposite polarity to the field in most of the outer core, usually does not have enough time to reverse before the field in the outer core changes again. One successful magnetic field reversal occurs near the end of our simulation.

The role of the Earth's mantle in controlling the frequency of geomagnetic reversals

A series of computer simulations of the Earths dynamo illustrates how the thermal structure of the lowermost mantle might affect convection and magnetic-field generation in the fluid core. Eight different patterns of heat flux from the core to the mantle are imposed over the core–mantle boundary. Spontaneous magnetic dipole reversals and excursions occur in seven of these cases, although sometimes the field only reverses in the outer part of the core, and then quickly reverses back. The results suggest correlations among the frequency of reversals, the duration over which the reversals occur, the magnetic-field intensity and the secular variation. The case with uniform heat flux at the core–mantle boundary appears most ‘Earth-like’. This result suggests that variations in heat flux at the core–mantle boundary of the Earth are smaller than previously thought, possibly because seismic velocity anomalies in the lowermost mantle might have more of a compositional rather than thermal origin, or because of enhanced heat flux in the mantles zones of ultra-low seismic velocity.

Numerical simulations of stellar convective dynamos. I. the model and method

Abstract A numerical model used to simulate global convection and magnetic field generation in stars is described. Nonlinear, three-dimensional, time-dependent solutions of the anelastic magnetohydrodynamic equations are presented for a stratified, rotating, spherical, fluid shell heated from below. The velocity, magnetic field, and thermodynamic perturbations are expanded in spherical harmonics to resolve their horizontal structure and in Chebyshev polynomials to resolve their radial structure. An explicit Adams-Bashforth time integration scheme is used with an implicit Crank-Nicolson treatment of the diffusion terms. Nonlinear terms are computed in physical space; and spatial derivatives are computed in spectral space. The resulting second-order differential equations are solved with a Chebyshev collocation method. Preliminary solutions for a Solar-like model are briefly discussed. Convective motions driven in the upper, superadiabatic part of the zone penetrate into the lower, subadiabatic part of the zone. Differential rotation, induced by the interaction of convection and rotation, is an equatorial acceleration at the surface as observed on the Sun; below the surface angular velocity decreases with depth. The meridional circulation and the equator-pole temperature excess are within Solar observational constraints; however, the giant-cell velocities at the surface are larger than observed. The large-scale magnetic field, induced by the differential rotation and helical motions, peaks in the subadiabatic region below the convection zone; and, near the top of the convection zone, it is concentrated over the downdrafts of giant cells. A systematic drift in latitude of the magnetic field or a field reversal has not yet been seen.

Numerical modeling of the geodynamo: Mechanisms of field generation and equilibration

Numerical calculations of fluid dynamos powered by thermal convection in a rotating, electrically conducting spherical shell are analyzed. We find two regimes of nonreversing, strong field dynamos at Ekman number 10 -4 and Rayleigh numbers up to 11 times critical. In the strongly columnar regime, convection occurs only in the fluid exterior to the inner core tangent cylinder, in the form of narrow columnar vortices elongated parallel to the spin axis. Columnar convection contains large amounts of negative helicity in the northern hemisphere and positive helicity in the southern hemisphere and results in dynamo action above a certain Rayleigh number, through a macroscopic α 2 mechanism. These dynamos equilibrate by generating concentrated magnetic flux bundles that limit the kinetic energy of the convection columns. The dipole-dominated external field is formed by superposition of several flux bundles at middle and high latitudes. At low latitudes a pattern of reversed flux patches propagates in the retrograde direction, resulting in an apparent westward drift of the field in the equatorial region. At higher Rayleigh number we find a fully developed regime with convection inside the tangent cylinder consisting of polar upwelling and azimuthal thermal wind flows. These motions modify the dynamo by expelling poloidal flux from the poles and generating intense toroidal fields in the polar regions near the inner core. Convective dynamos in the fully developed regime exhibit characteristics that can be compared with the geomagnetic field, including concentrated flux bundles on the core-mantle boundary, polar minima in field intensity, and episodes of westward drift.

THREE-DIMENSIONAL SPHERICAL SIMULATIONS OF SOLAR CONVECTION. I. DIFFERENTIAL ROTATION AND PATTERN EVOLUTION ACHIEVED WITH LAMINAR AND TURBULENT STATES

Rotationally constrained convection possesses velocity correlations that transport momentum and drive mean —ows such as diUerential rotation. The nature of this transport can be very complex in turbu- lent —ow regimes, where large-scale, coherent vorticity structures and mean —ows can be established by smaller scale turbulence through inverse cascades. The dynamics of the highly turbulent solar convection zone therefore may be quite diUerent than in early global-scale numerical models, which were limited by computational resources to nearly laminar —ows. Recent progress in high-performance computing tech- nology and ongoing helioseismic investigations of the dynamics of the solar interior have motivated us to develop more sophisticated numerical models of global-scale solar convection. Here we report three- dimensional simulations of compressible, penetrative convection in rotating spherical shells in both laminar and turbulent parameter regimes. The convective structure in the laminar case is dominated by ii banana cells,ˇˇ but the turbulent case is much more complex, with an intricate, rapidly evolving down- —ow network in the upper convection zone and an intermittent, plume-dominated structure in the lower convection zone and overshoot region. Convective patterns generally propagate prograde at low lati- tudes and retrograde at high latitudes relative to the local rotation. The diUerential rotation pro—les show some similarity with helioseismic determinations of the solar rotation but still exhibit signi—cantly more cylindrical alignment. Strong, intermittent, vortical down—ow lanes and plumes play an important dynamical role in turbulent —ow regimes and are responsible for signi—cant diUerences relative to laminar —ows with regard to momentum and energy transport and to the structure of the overshoot region at the base of the convection zone. Subject headings: convectionhydrodynamicsstars: interiorsSun: interiorSun: rotation ¨ turbulence

Rotation and Magnetism of Earth's Inner Core

Three-dimensional numerical simulations of the geodynamo suggest that a superrotation of Earths solid inner core relative to the mantle is maintained by magnetic coupling between the inner core and an eastward thermal wind in the fluid outer core. This mechanism, which is analogous to a synchronous motor, also plays a fundamental role in the generation of Earths magnetic field.

A numerical dynamo benchmark

We present the results of a benchmark study for a convection-driven magnetohydrodynamic dynamo problem in a rotating spherical shell. The solutions are stationary aside from azimuthal drift. One case of non-magnetic convection and two dynamos that differ in the assumptions concerning the inner core are studied. Six groups contributed numerical solutions which show good agreement. This provides an accurate reference standard with high confidence.

Effects of multiple phase transitions in a three‐dimensional spherical model of convection in Earth's mantle

Numerical models of mantle convection that incorporate the major mantle phase changes of the transition zone reveal an inherently three-dimensional flow pattern, with cylindrical features and linear features that behave differently in their ability to penetrate the 670-km discontinuity. The dynamics are dominated by accumulation of cold linear downwellings into rounded pools above the endothermic phase change at 670 km depth, resulting in frequent “avalanches” of upper mantle material into the lower mantle. The effect of the exothermic phase transition at 400 km depth is to reduce the overall degree of layering by pushing material through the 670-km phase change, resulting in smaller and more frequent avalanches, and a wider range of morphologies. Large quantities of avalanched cold material accumulate above the coremantle boundary (CMB), resulting in a region of strongly depressed mean temperature at the base of the mantle. The 670-km phase change has a strong effect on the temperature field, with three-distinct regions being visible: (1) the upper mantle, containing linear downwellings and pools of cold material in the transition zone and characterized by a high amplitude long wavelength spectrum; (2) the midmantle, containing quasi-cylindrical avalanche conduits and characterized by a low amplitude, broad spectrum; and (3) the deep mantle, containing large pools of cold, avalanched material and characterized by a high amplitude, ultra-red (i.e., long wavelength) spectrum. The effect on the velocity field is very different. Flow penetration across the 670-km phase change is strongly wavelength-dependent, with easy penetration at long wavelengths but strong inhibition at short wavelengths. Thus, when comparing numerical models with long wavelength seismic tomography, diagnostics based on the density field, such as the radial correlation function, are much more sensitive to the effects of phase transitions than those based on the velocity field. The amplitude of the geoid is not significantly affected by the partial layering, because the contribution from the strong heterogeneity in the transition zone is almost perfectly balanced by the contribution from deflection of the 670-km discontinuity. Avalanches are associated with geoid lows. However, a more complex viscosity structure is required to correctly match the sign of the geoid over slabs in Earth.

Three-Dimensional Spherical Models of Convection in the Earth's Mantle

Three-dimensional, spherical models of mantle convection in the earth reveal that upwelling cylindrical plumes and downwelling planar sheets are the primary features of mantle circulation. Thus, subduction zones and descending sheetlike slabs in the mantle are fundamental characteristics of thermal convection in a spherical shell and are not merely the consequences of the rigidity of the slabs, which are cooler than the surrounding mantle. Cylindrical mantle plumes that cause hotspots such as Hawaii are probably the only form of active upwelling and are therefore not just secondary convective currents separate from the large-scale mantle circulation. Active sheetlike upwellings that could be associated with mid-ocean ridges did not develop in the model simulations, a result that is in agreement with evidence suggesting that ridges are passive phenomena resulting from the tearing of surface plates by the pull of descending slabs.

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 Earths 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 a...