REVIEWS OF GEOPHYSICS, SUPPLEMENT, PAGES 443-450, JULY 1995 U.S. NATIONAL REPORT TO INTERNATIONAL UNION OF GEODESY AND GEOPHYSICS 1991-1994 Dynamics of the core, geodynamo Paul H. Roberts Institute of Geophysics and Planetary Physics, University of California. Los Angeles 1. Introduction standing affiliations with US Institutions. US theoreti­ cians have played a major role in elucidating fast dy­ T h e mechanism for generating the geomagnetic field namos, which amplify fields on the same time scales as remains one of the central unsolved problems in geo­ the flows. In contrast, slow dynamos act on a diffusive science." So states the report on the National Geo­ timescale, based on the magnetic diffusivity, n, of the magnetic Initiative (NGI) prepared by the U.S. Geo- conductor. Interestingly, McFadden and Merrill [1993] dynamics Committee, et al [1993], with advice from have recently derived 17 ~ 1 m 2 s"1 from the paleomag­ the NGI Workshop held in Washington D.C. in March netic data, but we shall take n = 3 m2 s""1. The diffu­ 1992. All analyses of the geomagnetic data point to the sive time scale of the core is therefore about 104 years, core as containing the source of the field and "The ba­ and the geodynamo problem is to understand how the sic premise that virtually everyone accepts is that the geomagneticfieldis maintained over times that are sub­ Earth's magnetism is created by a self-sustaining dy­ stantially longer than this. Fast dynamo theory has no namo driven by fluid motions in Earth's core" (NGI, obvious contributions to make, and will therefore not p.135). Dynamical questions at once arise, such as be considered here. The core is a slow dynamo. "What is the energy source driving those motions?" Ja­ cobs [1953] proposed that the solid inner core (SIC) is 2. C o r e D y n a m i c s ( M i c r o s c a l e ) the result of the freezing of the fluid outer core (FOC). Verhoogen [1961] noticed that the release of latent heat at the inner core boundary (ICB) during freezing would There has been increasing interest in the role played help drive thermal convection in the FOC, and Bragin­ by small-scale motions in transporting and mixing ther­ sky [1963] pointed out that the release of the light al­ mal and chemical inhomogeneities in the core. When loying elements during fractionation at the ICB would the heavy constituents (mainly Fe) of core fluid freeze provide compositional buoyancy. These two sources suf­ onto the ICB, latent heat and the light constituents are fice to supply the geodynamo with energy throughout released. Moffatt [1989] proposed that this light, hot geological time, even in the absence of dissolved radioac­ fluid would congregate into "blobs" at the ICB that tivity in the core [Braginsky and Roberts, 1994a; Kuang would, when large enough, break away from the ICB et al, 1994]. Stevenson [1991] argues that potential dif­ and rise through the core, stirring it as it did so, and ferences on the core-mantle boundary (CMB) of elec­ possibly retaining their identity until they reach the trochemical origin may be partially responsible for the CMB, where some may remain to form a light layer geomagnetic field. (see §4 below). This idea has been pursued by Ruan Whatever the driving mechanism, it is clear that the and Loper [1993], Loper and Moffatt [1993] and Moffatt magnetohydrodynamics (MHD) of the core must be un­ derstood. This has proved to be a challenging task; and Loper [1994]; see also Bush et al [1992,1994], Loper et al [1994] progress has been slow. The directional property of Braginsky and Meytlis [1990] argued that core tur­ the magnetic compass needle demonstrates that Corio- bulence is totally unlike classical turbulence of Kol- lis forces play an essential role. Because the molecular mogoroff type, in which energy 'cascades' from large diffusivities of heat and composition are so small, these to small eddies. Nor is it related to classical MHD tur­ sources of buoyancy must instead be transported across bulence, where a reverse cascade may create large-scale the core by turbulence. The phenomena of interest arise magnetic fields by turbulent dynamo action, possibly from slight deviations in the FOC from a well-mixed through a turbulent a-effect. [The a-effect is the cre­ adiabatic state, and theory must consistently disentan­ ation of a mean electromotive force (emf) parallel to the gle these from the large "background." In short, it is far mean magnetic field.] They argue that the microscale from obvious what equations best describe large scale fields are so tiny that they produce neither a significant core MHD [Braginsky and Roberts, 1994a]. turbulent or—effect, nor an enhancement in the mean We shall concentrate below on the research of U S field diffusivity. Nevertheless, the turbulent diffusivi­ scientists, even when it was carried out abroad. We ties of the mean thermal and chemical inhomogeneities shall add the work of foreign scientists who have long are enormously greater than their molecular counter­ parts, at least in some directions. Because of the strong influence of Coriolis and Lorentz forces on motions of Copyright 1995 by the American Geophysical Union. all scales, the turbulence is highly anisotropic, forming 'plate-like' eddies that have their long dimensions par­ Paper number 95RG00735. allel to the rotation axis (Oz) and to the mainly zonal 8755-1209/95/95RG-00735$15.00 (^—)direction of the prevailing toroidalfield;the short 443 444 ROBERTS: DYNAMICS OF THE CORE, GEODYNAMO dimension is in the s—direction, away from the rotation larity reversal mechanism. Some MC instabilities are of axis. Turbulent diffusion is represented by one tensor, short time scale, of the order of 103 years. These are the the same for both heat and composition. The elements so-called ideal instabilities, where "ideal" refers to the of that tensor corresponding to diffusion in the z and fact that, unlike the so-called resistive instabilities, they directions are large, of the same order as the molec­ do not rely on the finite resisitivity of the fluid. Ideal ular magnetic diffusivity, 17; turbulent diffusion in the and resistive instabilities are the counterparts of sim­ s—direction is comparatively weak. Braginsky-Meytlis ilar instabilities that arise in laboratory plasmas but, theory has recently been taken further by Braginsky and because of the importance of Coriolis forces in the core, Roberts [1994a]; it still contains ad hoc elements. they evolve there on longer time scales (see above). It is the highly dispersive character of rotating flu­ London [1992a, b] examined MC waves and instabili­ ids that led to this unusual picture of core turbulence. ties, supposing that the prevailing magnetic field, B, is The effect of that dispersion on the Moffatt mecha­ zonal and proportional in strength to distance, s, from nism has recently been studied by St Pierre [1994a, the rotation axis, Oz. Assuming a geostrophic dynam­ b] who argues, on the basis of computer simulations, ical balance of the kind used in atmospheric dynamics, that a blob released at the ICB will be stretched, lami­ he developed a uniform approximation for waves that nated in plates, and absorbed into its surroundings be­ have a short wavelength in the s—direction, i.e. away fore it can rise far into the FOC; see also St Pierre from the rotation axis. He showed [London, 1992b] that and Roberts [1994]. It is difficult in the laboratory to these waves propagate in a westward direction. In a mimic core conditions in which the effective diffusivities later work [London, 1993], he generalized to other zonal of heat and composition are (see above) the same. Ex­ fields; see also London [1994]. periments have however been performed by Cardin and Stablefieldconfigurations may become unstable when Olson [1992]. the fluid is top heavy, and there is much interest there­ fore in MAC waves and instabilities, where the added 3. Core Dynamics (Macroscale) 'A* stands for Mrchimedean' (i.e. buoyancy) forces. Simple examples of MC and MAC instabilities have Nearly all existing models of the geodynamo are ax­ isymmetric. Since the axisymmetric part of the ge­ omagnetic field cannot, according to Cowling's theo­ rem, be self-maintained, the emf created by the nonaxisymmetric components of fluid flow and magnetic field must be retained. For simplicity, it is usually parameterized by an a—effect. Since turbulent induc­ tion is unimportant (see §2), the a—effect is created by asymmetric waves/instabilities of global scale. Recent work on axisymmetric geodynamo models is reported in §5 below; here we describe studies of waves/instabilities. Viewed from the inertial frame, a contained rotat­ ing fluid is filled with vortex lines parallel to the ro­ tation axis that impart the "elasticity" that inertial waves require. In a non-rotating electrically conducting fluid, the "elasticity" of the lines of force of the pre­ vailing magnetic field is responsible for Alfven waves. In a rapidly rotating conductor, the Alfven and inertial waves are replaced by 'fast* and 'slow' waves. The for­ mer resemble the inertial waves and have a timescale of the order of a day; the latter act on timescales of order r8 = 2QjAopR2/B%1 where Q « 7 x 10"5s""1 is the angu­ lar velocity of Earth, fio « 4w x 10"7H m"1 is the mag­ been analyzed by Kuang and Roberts [1991,1992], Lan, Kuang and Roberts [1993]. Fearn and Kuang [1994] and Kuang [1994] stress the importance of the conductivity of the boundaries on the instabilities. Bergman and Madden [1993] studied core convection, paying particular attention to the steady mean poloidal circulation in the core, for which they argued equato­ rial upwelling would occur. Such a circulation has a profound effect on the functioning of an ocu—dynamo; see §5. Compositional buoyancy is important for driv­ ing core convection in the considerations of Kuang et al [1994], and Bergman et al [1994]. The surface of the inner core is probably constitutionally supercooled, so that a mushy layer exists there. If no magnetic field is present, chimneys form in such a layer through which the light fluid, released during fractionation inside the layer, is ejected into the FOC. Bergman et al investi­ gate how this mechanism is affected by the prevailing magnetic field. For simplicity, many investigations of core dynamics and the geodynamo ignore the SIC, by assuming that the entire core is fluid. At first sight, this unrealism seems not too serious: the SIC is only 4% of the vol­ ume of the core and 5% of its mass. Nevertheless, the netic permeability, p « 104kg m~3 is the core density, SIC may have a disproportionate effect on core flows R « 3.5 x 106m is the radius of the core and B is a char­ and field generation. It has long been known [Stewart- acteristic strength of the (zonal) field. If B at the CMB that creates an electromagnetic the T C and its adjustment layer may be even greater for torque that best fits the length of day data. They make the asymmetricfields;see also Hollerbach [1994]. Glatz­ three demands which they find cannot be simultane­ maier and Olson [1993] studied non-magnetic convec­ ously met: (a) Br(R) does not exceed the upper limit tion in a rotating sphere and showed that the amplitude provided by dynamo theory (see also Levy and Pearce of the convective motions is greatest outside the TC; see [1991] who argue that B T ( R ) is less than lOmT, and also Cardin and Olson [1994a]. In contrast, for the cor­ is probably less than lmT), (b) the poloidal electric responding MHD situation, where a zonal magnetic field currents which generate that toroidal field and which was imposed, Olson and Glatzmaier [1993, 1994] and leak into the mantle do not exceed bounds on the elec­ Glatzmaier and Olson [1994] found that the convection tricfieldinferred from measurements at Earth's surface was strongest inside the TC, the Taylor columns outside [Lanzerotti et al, 1992, 1993, 1994], (c) the ohmic dis­ that cylinder being suppressed by the Lorentz force; see sipation in the mantle caused by those currents does also Jones et al [1994], Cardixt and Olson [1994b] and not exceed the heat flux from the Earth. They con­ §5 below. clude that magnetic stresses cannot be the main factor in core-mantle coupling. Their treatment offluxdiffu­ 4. Effects of t h e M a n t l e o n C o r e M H D sion in the analysis leading to their conclusion merits further study. Love and Bloxham [1994b] have recently proposed a second application of their idea. The core is only one component of the coupled core- Diffusion offluxplays an important role in the study mantle system. Each component profoundly affects, of Braginsky and Le Mouel [1993], who are particularly and is affected by, the other. Strictly, the core cannot interested in the inductive effects of high shears in a be considered in isolation from the mantle but, when it "A—layer" at the top of the core. Kuang and Bloxham is, the mantle is replaced by a set of conditions on the [1993] analyze how magnetic field in the upper core af­ CMB. The resulting theoretical simplification is enor­ fects topographic core-mantle coupling. Theyfindthat mous, but sometimes is an over-simplification. In par­ the field may change the strength of the coupling by ticular, to suppose that the form, and the physical state, of the CMB are uniform in space and unvarying in time is simplistic. Larson and Olson [1991] argue that vari­ several orders of magnitude but, for parameters appro­ priate to the core, the stress is of order 10~1N m~2, ations in the convective regime in the mantle, and in which is adequate to account for the decadal variations particular the changing configuration of mantle plumes, in Earth's rotation. Angular momentum exchange be­ control the rate of geomagnetic field reversals. It has tween core and mantle is also discussed by Bloxham and recently been realized that lateral variations in the tem­ Kuang [1994]. perature of the CMB will have a strong effect on core Malkus has long urged that precessional driving of motions and therefore on core-mantle coupling; see Sun the core is an important and, possibly, the dominant et al [1994]. A new type of geodynamo is also possible; source of energy for core motions and the geodynamo; see §5. In this section we shall ignore lateral variations e.g. see Malkus [1994]. Interest in this idea has been on the CMB, apart from topography. revitalized by the discovery that flows with elliptical The assumption (§1) of an adiabatic well-mixed core streamlines, somewhat similar to flows driven by the becomes suspect near the CMB, and several authors luni-solar precession, are unstable. Malkus has pro­ have argued that a layer of comparatively light fluid vided experimental demonstrations of the instability in exists adjacent to the CMB. Braginsky [1993] has chris­ an elliptically distorted cylinder of fluid. Experiments tened this "the hidden ocean of the core" and has ar­ have also been performed by Vanyo [1991], Vanyo et gued that waves propagating in this stable layer may be al [1992, 1994b] and Wilde and V anyo [1994]; see also partially responsible for the short period geomagnetic Vanyo et al [1994a], Vanyo and Lods [1994]. So far, secular variation. It may also strongly affect core-man­ all studies have been non-magnetic, but it is hoped tle coupling, particularly topographic coupling. Waves that they will provide stepping stones to the corre­ in a stratified layer at the top of the core have been sponding MHD situations. The a—effect created by studied by Bergman [1993], who developed a theory of precessionally-driven flows has already been estimated magnetic Rossby waves based on a generalization of by Barenghi et al [1994]. The effects of the SIC on the Laplace's tidal equation in which the Lorentz force is forced nutation of the Earth have been studied theoret­ included and the induction equation is added. He also ically, and the results have been compared with obser­ developed /?—plane solutions analytically and showed vational data by Mathews et al [1991a, b] and Herring that the magneticfieldcan release equatorially trapped et al [1991]. Cognate issues are analyzed by de Vries Rossby waves. and Wahr [1991]. Love and Bloxham [1994a] have recently investigated The exchange of the z—component of angular mo­ a new idea which may lead to the abandonment of mag­ mentum between core and mantle is accomplished via netic core-mantle coupling in comparison with other, geostrophic motions in the core; these are zonal flows 446 ROBERTS: DYNAMICS OF THE CORE, GEODYNAMO that depend only on distance s from the rotation axis, therefore be strong. It used to be said that an appeal to Oz, and on time t. By analyzing the field extrapolated an invisible zonal field is a return to armchair science, downwards to the CMB, Jault and Le Mouel [1988] esti­ but galaxies are transparent to observation. Their fields mated the geostrophicflowin the recent past and could are predominantly toroidal and seem to be produced by therefore monitor the angular momentum of the core an aw—mechanism; see §7 of Krause et al [1993]. It is as a function of time. They showed that its variations in principle possible to detect a zonal field in Earth's are roughly equal but opposite to those of the angular core through the electric fields it creates outside the momentum of the mantle over the same period, as deter­ core, in particular the potential difference between the mined by changes in the length of day; the net angular two ends of a trans-oceanic cable [Lanzerotti et al, 1992, momentum of the core-mantle system as a whole is con­ 1993, 1994]. In practice, the obscuring, long period, in­ stant. Jackson ei al [1993] have developed this theme ductive effects of ocean currents have so far prevented a and have used their analysis of the core geostrophic flow convincing demonstration [Runcorn and Winch, 1991). to predict, with encouraging results, variations in the The axisymmetric force balance in an MHD dynamo length of day. is not easily understood. Many studies of 2D "inter­ mediate" dynamo models have been launched to elu­ 5. G e o d y n a m o Modeling cidate it. These are so named because, while they do not address the full MHD problem, they take a Dynamo models that solve the induction equation alone are called 'kinematic', and several such geodynamos models have been integrated. The main chal­ lenge today is to solve the fully dynamic dynamo prob­ lem, sometimes also called 'the MHD dynamo problem' or 'the fully self-consistent dynamo problem', in which the induction equation is solved and the equation of motion for the fluid. This nonlinear problem raises formidable difficulties. Because of Cowling's theorem, a true MHD dynamo model should be 3D, but a super­ step beyond kinematic models. An a—source is in­ voked to maintain the axisymmetric field. Whether an a2—, a w - or a2w—dynamo results depends on the dy­ namical balance. In an a2—dynamo the primary bal­ ance is geostrophic, i.e. between Coriolis and pressure forces; the magneticfieldstrength, B, is determined by a secondary balance, e.g. between the Lorentz and vis­ cous forces, which gives B ~ (fiopi/rj)1/2/R) where v is core viscosity. In a strong field dynamo, the mag­ netic field plays a role in the primary balance and computer is then required to integrate it numerically. therefore B = ^(B^BM) ~ (2f2f7/i0/?)1/2 ~ 2.4mT, Axisymmetric (2D) models can be solved on worksta­ a relation confirmed by Benton [1992]. Benton ar­ tions, but the magnetic field will decay to zero unless gued further that a typical zonal flow velocity would the emf generated by the omitted asymmetricfieldsand be (2Qn2/R)1/3 - 7 x 10"*4m s""1, which is similar to flows is reintroduced in some way, through an a—effect. the speed with which some features of the geomagnetic Axisymmetric a—effect models are of two extreme field drift westward. types, a2— and aw—dynamos, together with a range St. Pierre [1993b] demonstrated that, when a strong of a2w—models between them. In an a2—model, zonal field branch exists, weak field solutions are likely to be field creates meridional field and vice versa; in an aw- nonlinearly unstable. Although his plane layer model is model, the a-effect creates meridional field from zonal geometrically too simple to represent the geodynamo, field, but the zonalfieldis created by an w—effect, i.e. by it is a convective MHD dynamo of similar physical the inductive effects of zonal shearing motions. Some­ type. It clearly demonstrates the existence of a strong times aw—models are called "strong field dynamos" field branch, one that also operates subcritically, i.e. at since the zonal field, which is locked inside the con­ smaller thermal forcing than that at which kinematic ductor, is large compared with the observed merid­ dynamo action is first possible. The model of St. Pierre ional field, in contrast to the "weak field" dynamos is fully 3D, as is the spherical 3D model of Glatzmaier of a2—type where the strength, BM , of the meridional and Roberts [1994] described below. field is characteristic of the strength of the entire field. The axisymmetric force balance in intermediate mod­ A strong field dynamo functions only if the product of els is so dominated by magnetic and Coriolis forces that a and w, as measured by the so-called "dynamo num­ inertial forces are often omitted. One of two extreme ber" , D, exceeds in magnitude a certain marginal value, scenarios arise, or perhaps some intermediate scenario. D M ; for an a2—model to function the a—effect magnetic At one extreme is the model—Z state [Braginsky, 1975, Reynolds number, Ra, a dimensionless measure of a, 1991,1994] which relies on the coupling of core to man­ must be large enough. While a2—dynamos are usually tle and in which the geostrophic motions in the core steady, aw—dynamos tend to be oscillatory, but they are large. At the other extreme is the Taylor state too may become steady, and more efficient (as judged [Taylor, 1963] in which core-mantle coupling is insignif­ by a smaller value of D M ) , when a sufficiently strong icant, but in which a certain integral demand (the Tay­ meridional flow is present. Meridionalflowis produced lor constraint) must be satisfied. Sometimes models of by Lorentz forces or by core-mantle coupling; see also either type can exist under the same conditions of exci­ Bergman and Madden [1993]. tation. Model—Z is energetically the more expensive to Zonal shearing motion is comparatively easily excited run, because of core-mantle friction, and the external in rotatingfluids,for example by pole-equator temper­ fields it produces therefore tend to be smaller than in ature differences; the w—effect and the zonal field may the corresponding Taylor-type model. There are there- ROBERTS: DYNAMICS OF THE CORE, GEODYNAMO 447 fore two contenders for the geodynamo, a strong field of anisotropic aw—models; see also Kono and Roberts (model—Z) mechanism and a very strong field (Taylor- [1994]. A general method for solving weakfieldMHD type) mechanism. In trying to decide between these, models when conditions for kinematic dynamo action it is usually supposed for simplicity that core-mantle are only marginally exceeded, has been explored by coupling is viscous — it is the existence of this cou­ Kono and Roberts [1991, 1992]. pling rather than its precise nature that is significant. Lateral variations in the temperature of the CMB St. Pierre [1993a] has examined the stability of Taylor bring about concomitant changes in the electrical con­ states. ductivity of the lower mantle so that new current paths This then is the background against which much of are allowed and old ones forbidden. The axisymmetry the recent work on intermediate geodynamos may be assumed in most geodynamo modeling is destroyed and viewed. Hollerbach and Ierley [1991] analyzed an in­ with it the applicability of Cowling's theorem [Busse, termediate dynamo of a2—type and showed that, as 1992]]. A zonal shear can readily create zonal magnetic Ra exceeds its marginal value, Ram > the solution is at field from meridional field through the w—effect, but first viscously controlled. As Ra is further increased, a it is incapable, in an axisymmetric system, of creating second critical value, iZaT, is reached at which Taylor meridional field from zonal field. This, however, is no states appear. When Hollerbach et al [1992] carried longer true when longitudinal inhomogeneities destroy out a parallel study for an aw—model, they uncovered the axial symmetry. A zonal shear can then produce a more complex situation. Despite very simple choices zonal field from meridional field and vice versa. This of a and ut they found that, as the dynamo number D fact enabled Busse and Wicht [1991] and Wicht and increases beyond D m , the solution is at first viscously Busse [1993] to construct new, simple models of dy­ controlled but that, as D increases through a second namo action that make use of the broken symmetry critical value DXi oscillations arise in which the Taylor and which work through zonal shear alone. balance is struck during part of each cycle but in which An unusual approach to the geodynamo problem was viscous coupling is essential during the remainder. initiated by Ruzmaikin et al [1993]. They divide the Braginsky and Roberts [1994b] continued earlier in­ fluid domain into fixed cells, each of which randomly vestigations of one particular model. They observed a amplifies or destroys field by dynamo action; nonlin- transition from Taylor-type behavior to model—Z—type earity, diffusion and correlations between cells are then behavior as D increases. An au;—dynamo model inte­ added. An initially smooth field becomes intermittent, grated by Glatzmaier and Roberts [1993] developed an the field concentrating mainly in a few cells, the loca­ interesting bifurcation as D was increased. Against the tion of which changes with time, a phenomenon they background of an approximately steady dipole compo­ liken to the motion of geomagnetic field anomalies. nent, an oscillatory quadrupole field causes the merid­ Roberts [1992] introduced a "mapping method" that has ional field lines to bunch up alternately in one hemi­ been successfully tested against axisymmetric [Naka- sphere and the other. The role of the dipole and jima and Roberts, 1994a] and asymmetric [Nakajima quadrupole families of solutions of the geodynamo equa­ and Roberts, [1994b] dynamo models; see also Nakajima tion in geomagnetic field reversals is an oft recurring et al [1993a]. theme, and was raised again in a novel way by Hoffman The resources of the NSF Pittsburg Supercomputing [1991]. Questions of parity coupling in a2—models also Center were enlisted to generate the first 3D time- arose in the work of Hollerbach [1991]. dependent, fully self-consistent numerical solution of In integrating a kinematic ao;—geodynamo model, the MHD equations that describes thermal convection Braginsky [1964] found that the fields induced outside and magnetic field generation in a low-viscosity rap­ the TC differed substantially from those generated in­ idly-rotating spherical shell with a solid conducting in­ side it. Dynamic (intermediate) models have recently ner core. The resulting solution, reported by Glatz­ been studied by Hollerbach and Jones [1993a]. They maier and Roberts [1994], serves as a crude simulation find that most dynamo action takes place outside the of the geodynamo, crude because because the trunca­ TC, a conclusion that may depend on their, choices of a tion was too severe and because geophysicially realis­ and o> since it was not confirmed by the recent 3D inte­ tic values for some parameters were not numerically grations of Glatzmaier and Roberts [1994]. The model accessible (e.g. v was several orders of magnitude too of Hollerbach and Jones was used to benchmark that large, though still apparently not very influential). The of Glatzmaier and Roberts [1993]; the agreement was heat flux from the core was taken to be 4 x 1013W and nearly perfect. Hollerbach and Jones [1993b, 1994] ar­ the integration was continued over approximately three gued that the SIC plays a potent role in the reversal magnetic diffusion times, during which the field showed mechanism; its electromagnetic inertia diminishes chaos no signs of disappearing. Field generation takes place in the FOC. Glatzmaier and Roberts [1994] agreed. The mainly within and near the T C . The pattern and am­ effects of conducting boundaries were investigated by plitude of the radial magnetic field at the CMB is qual­ Hirsching and Busse [1993]. itatively similar to that of the Earth. The toroidal field Although the emphasis of the subject has moved to­ energy is rather larger than the poloidalfieldenergy but wards MHD models, kinematic geodynamos are still be­ the maximum amplitudes attained by the two fields are ing profitably studied. In particular, Hagee and Ol­comparable (~ 0.05T); the maximum fluid velocity is son [1991] have suggested an interesting connection be­ of order 4 x 10~3m s - 1 . An irregular exchange of field tween the observed secular variation and certain types between hemispheres takes place, similar to that found 448 ROBERTS: DYNAMICS OF THE CORE, GEODYNAMO in the 2D model of Glatzmaier and Roberts [1993]. Its an inviscid drop in a bounded rotating fluid, Phys. Fluids timescale is about 10% of the magnetic diffusion time of the FOC. Excitingly, the dynamo sometimes reverses its polarity spontaneously. Preliminary to doing so, the poloidal field in the SIC has to reverse; if it does not do so, the reversal is aborted (as in geomagnetic excur­ sions). It is hard not to be excited by such similarities A4, 1142-1147, 1992. Bush, J.W.M., H.A. Stone and J. Bloxham, Axial drop mo­ tion in rotating fluids, J. Fluid Mech., t o appear, 1994. Busse, F.H., Theory of the geodynamo and core-mantle coupling, in Chaotic Processes in the Geological Sciences, edited by D.A. Yuen, pp. 281-292, Springer, Heidelberg, 1992. between the computed model and the real Earth. Per­ Busse, F.H., and J. Wicht, A simple dynamo caused by con­ haps the answer to the challenging sentence that opened ductivity variations, Geophys. Astrophys. Fluid Dynam., this review is at last within sight? 64, 135-144, 1992. Cardin, P., and P. Olson, An experimental approach to ther- mochemical convection in the Earth's core, Geophys. Res. References Lttrs., 19, 1995-1998, 1992. Cardin, P., and P. Olson, Chaotic thermal convection in a rapidly rotating spherical shell: consequences for flow Barenghi, C F . , R. Hollerbach and R.R. Kerswell, Alphaeffects produced by precessionally induced flows, Abstract, 4th SEDI Symposium, Whistler, Session 3: Core Dynam­ ics and Thermodynamics, 1994. Benton, E.R., Hydromagnetic scale analysis at Earth's coremantle boundary, Geophys. Astrophys. Fluid Dynam., 67, 259-272, 1992. Bergman, I., Magnetic Rossby waves in a stably stratified layer near t h e surface of the Earth's inner core, Geo­ phys. Astrophys. Fluid Dynam., 68, 151-176, 1993. Bergman, M.I., D.R. Fearn, P. Shannon and J. Bloxham, The effect of a magnetic field on chimney convection in a solidifying, metallic alloy melt, Abstract, 4th SEDI Sym­ posium, Whistler, Session 3: Core Dynamics and Ther­ modynamics, 1994. Bergman, M.I., and T.R. Madden, Steady, finite amplitude, rotating magnetoconvection and the mean poloidal circu­ lation of t h e Earth's core, Phys. Earth planet. Int., 77, 267-284, 1993. in t h e outer core, Phys. Earth planet. Int., 82, 235-259, 1994a. Cardin, P., and P. Olson, The influence of toroidal mag­ netic field o n thermal convection in t h e core, Abstract, 4th SEDI Symposium, Whistler, Session 3: Core Dynam­ ics and Thermodynamics, 1994b. de Vries, D., and J.M. Wahr, The effects of the solid in­ ner core on the Earth's forced nutations and Earth tides, /. Geophys. Res., 96, 8275-8293, 1991. Fearn, D.E,, and W. Kuang, Resistive instabilities in a rap­ idly rotating fluid without critical layers, Geophys. Astr­ ophys. Fluid Dynam., 94, 181-206, 1994. Glatzmaier, G.A., and P. Olson, Highly supercritical ther­ mal convection in a rotating spherical shell: Centrifugal v s radial gravity, Geophys. Astrophys. Fluid Dynam., 70, 113-136, 1993. Glatzmaier, G.A., and P. Olson, A study of the struc­ ture of convection in t h e Earth's liquid core, Abstract S4.I8, IASPEI 27th General Assembly, Wellington, New Bloxham, J., and W. Kuang, Core-mantle coupling and the Zealand, August 1994. exchange of angular momentum, Abstract, 4 SEDI Sym­ Glatzmaier, G.A., and P.H. Roberts, Intermediate dynamo posium, Whistler, Session 2: Core Mantle Interactions and models, J. Geomag. Geoelec, 45, 1605-1616, 1993. the CMB, 1994. Glatzmaier, G.A., and P.H. Roberts, A three-dimensional Braginsky, S.L, Structure of the F-layer and reasons for convective dynamo solution with rotating and finitely con­ convection in the Earth's core, Doklady Akad. Nauk SSSR, ducting inner core and outer mantle, Abstract, 4th SEDI 149, 1311-1314, 1963. Symposium, Whistler, Session 4: T h e Geodynamo and Braginsky, S.L, Kinematic models of the Earth's hydromag­ Observational Constraints, 1994. netic dynamo, Geomagn. & Aeron., 4* 5 7 2 - 5 8 3 , 1964. Hagee, V.L., and P. Olson, Dynamo models with permanent Braginsky, S.L, Nearly axisymmetric model of the hydro- dipole fields and secular variation, / . Geophys. Res., 96, magnetic dynamo of t h e Earth, Geomag. & Aeron., 15, 11673-11687, 1991. 122-128, 1975. Herring, T.A., B.A. Buffett, P.M. Mathews and LI. Shapiro, Braginsky, S.L, Towards a realistic theory of the geodynamo, Geophys. Astrophys. Fluid Dynam., 60, 8 9 - 1 3 4 , 1991. Braginsky, S.L, MAC oscillations and the hidden ocean of the core, J. Geomag. Geoelec., 45, 1517-1538, 1993. Braginsky, S.L, The nonlinear dynamo and model-Z, in Lectures on Solar and Planetary Dynamos, edited by M.R.E. Proctor and A.D. Gilbert, pp. 267-304, University Press, Cambridge U.K., 1994. Braginsky, S.L, and J.L. Le Mouel, Two-scale model of a geomagnetic field variation, Geophys. J. Int., 112, 1 4 7 158, 1993. Braginsky, S.L, and V.P. Meytlis, Local turbulence in the Earth's core, Geophys. Astrophys. Fluid Dynam., 55, 7 1 87, 1990. Braginsky, S.L, and P.H. Roberts, Equations governing con­ Forced nutations of the Earth: influence of inner core dy­ namics. 3. Very long baseline interferometry data analy­ sis, / . Geophys. Res., 96, 8259-8273, 1991. Hirsching, W., and F.H. Busse, A numerical study of dy­ namos in spherical shells with conducting boundaries, in Solar and Planetary Dynamos, edited by M.R.E. Proctor, P.C. Matthews and A.M. Rucklidge, pp. 137-144, Univ­ ersity Press, Cambridge U.K., 1993. Hoffman, K.A., Long-lived transitional states of the geo­ magnetic field and the two dynamo families, Nature 354, 273-277, 1991. Hollerbach, R., Parity coupling in a 2 - d y n a m o s , Geophys. As­ trophys. Fluid Dynam., 60, 245-260, 1991. Hollerbach, R., Magnetohydrodynamic Ekman and Stew- vection in Earth's core and the geodynamo, Abstract, 4th artson layers in a rotating spherical shell, Proc. R. Soc. SEDI Symposium, Whistler, Session 4: T h e Geodynamo Lond., A444, 333-346, 1993. and Observational Constraints, 1994a. Hollerbach, R., Imposing a magnetic field across a non- Braginsky, S.L, and P.H. Roberts, From Taylor state to axisyrametric shear layer in a rotating spherical shell, m o d e l - Z , Abstract, 4th SEDI Symposium, Whistler, Ses­ Phys. Fluids, t o appear, 1994. sion 4: The Geodynamo and Observational Constraints, Hollerbach, R., C F . Barenghi and C. A. Jones, Taylor's 1994b. constraint in a spherical aw-dynamo, Geophys. Astroph­ Bush, J.W.M., H.A. Stone and J. Bloxham, The motion of ys. Fluid Dynam., 67, 3-25, 1992. ROBERTS; DYNAMICS OF THE CORE, GEODYNAMO 449 Hollerbach, R. and G.R. Ierley, A modal or2—dynamo in the limit of asymptotically small viscosity, Geophys. Astrophys. Fluid Dynam., 56, 133-158, 1991. Hollerbach, R., and C. A. Jones, A geodynamo model in­ corporating a finitely conducting inner core, Phys. Earth planet. Int., 75, 317-327, 1993a. Hollerbach, R., and C. A. Jones, Influence of the Earth's inner core on geomagnetic fluctuations and reversals, Nature, 865, 541-543, 1993b. Hollerbach, R., and C. A. Jones, On the magnetically stabi­ lizing role of Earth's inner core, Phys. Earth planet. Int.9 to appear, 1994. Hollerbach, R., and M.R.E. Proctor, Non-axisymmetric sh­ ear layers in a rotating spherical shell, in Solar and Plane­ tary Dynamos, edited by M.R.E. Proctor, P.C. Matthews and A.M. Rucklidge, pp. 145-152, University Press, Cam­ bridge U.K., 1993. Jackson, A., J. Bloxham and D. Gubbins, Time dependent flow at the core surface and conservation of angular mo­ mentum in the coupled core-mantle system, in Dynamics of Earth's Deep Interior and Earth Rotation, edited by J.-L. Le Mouel, D.E. Smylie and T. Herring, pp. 97-107. AGU Monograph 72, American Geophysical Union, Wash­ ington DC., 1993. Jacobs, J.A., T h e Earth's inner core, Nature 172, 297-298, 1953. Jault, D. and J.-L. Le Mouel, Westward drift, core motions and exchanges of angular momentum between core and mantle, Nature SSS, 353-356, 1988. Jones, C.A., A.W. Longbottoms, K. Zhang and R. Holler­ bach, Convection in the outer core of the Earth, Abstract, 4th SEDI Symposium, Whistler, Session 4: T h e Geody­ namo and Observational Constraints, 1994. Kono, M., and P.H. Roberts, Small amplitude solutions of the dynamo problem. 1. The adjoint system and its solu­ tions, / . Geomag. Geoelec., 42, 839-862, 1991. Kono, M., and P.H. Roberts, Small amplitude solutions of the dynamo problem. 2. The case of the a2—dynamos, Geophys. Astrophys. Fluid Dynam., 67, 6 5 - 8 5 , 1992. Kono, M., and P.H. Roberts, Kinematic dynamos with anisotropic or—effect, Geophys. Astrophys. Fluid Dynam., to appear, 1994. Krause, F., K.-H. Radler and G. Rudiger (editors), The Cos­ mic Dynamo, IAU Symposium 157, Kluwer, Dordrecht, 1993. Kuang, W., Effect of boundary conditions on resistive insta­ bilities in rapidly rotating fluids, Geophys. Astrophys. Fluid Dynam., to appear, 1994. Kuang, W., and J. Bloxham, On the effect of the boundary topography on flow in the Earth's core, Geophys. Astro­ phys. Fluid Dynam.f 72, 161-195, 1993. Kuang, W., J. Bloxham and D.E. Fearn, Thermal and com­ positional convection in Earth's fluid outer core, Abstract, 4th SEDI Symposium, Whistler, Session 3: Core Dynam­ ics and Thermodynamics, 1994. Kuang, W., and P.H. Roberts, Resistive instabilities in rapidly rotating fluids: linear theory of the g—modes, Geophys. Astrophys. Fluid Dynam., 60, 295-333, 1991. Kuang, W., and P.H. Roberts, Resistive instabilities in rapidly rotating fluids: linear theory of the convective modes, Geophys. Astrophys. Fluid Dynam., 67, 129-162, 1992. Lan, S., W. Kuang and P.H. Roberts, Ideal instabilities in rapidly rotating MHD systems that have critical layers, Geophys. Astrophys. Fluid Dynam., 69, 131-160, 1993. Lanzerotti, L.J., How strong is the invisible component of the magnetic field in Earth's core? Geophys. Res. Lttrs., 21, 2339-2340, 1994. Lanzerotti, L.J., C.H. Sayres, L.V. Medford, J.S. Kran and C.G. Maclennan, Earth potential over 4000km between Hawaii and California, Geophys. Res. Lttrs., 19, 11771180, 1992. Lanzerotti, L.J., R.A. Langel and A.D. Crane, Geoelectromagnetism, in Encyclopedia of Applied Physics 7, edited by G.L. Trigg, pp. 109-123, VCH Publishers, New York, 1993. Larson, R.L., and P. Olson, Mantle plumes control magnetic reversal frequency, Earth planet. Science Lttrs., 107, 4 3 7 447, 1991. Levy, E.H., and S.J. Pearce, Steady state toroidal magnetic field at Earth's core-mantle boundary, / . Geophys. Res., 96, 3935 - 3942, 1991. London S.D., Hydromagnetic waves in a rotating sphere, Geophys. Astrophys. Fluid Dynam., 63, 9 1 - 1 0 9 , 1992a. London S.D., A uniform asymptotic approximation for un­ stable waves in a rotating sphere, Geophys. Astrophys. Fluid Dynam., 67, 185-198, 1992b. London S.D., Stability and trapping of slow hydromagnetic waves in rotating cylindrical geometry, J. Geomag. Geo­ elec., 45, 1553-1561, 1993. London S.D., Stability and trapping of slow hydromagnetic waves in rotating spherical geometry, Abstract, 4TH SEDI Symposium, Whistler, Session 4: T h e Geodynamo and Observational Constraints, 1994. Loper, D.E., and H.K. Moffatt, Small-scale hydromagnetic flow in the Earth's core: Rise of a vertical buoyant plume, Geophys. Astrophys. Fluid Dynam., 68, 177-202, 1993. Loper, D.E., H.K. Moffatt and H. Shimizu, Magnetostrophic dynamics of localized upwellings near the C M B , Abstract, 4TH SEDI Symposium, Whistler, Session 2: Core Mantle Interactions and the CMB, 1994. Love, J.J., and J. Bloxham, Electromagnetic coupling and the toroidal magnetic field at the core-mantle boundary, Geophys. J. Int., 117, 235-256, 1994a. Love, J.J., and J. Bloxham, Dynamos by inverse theory, Abstract, 4th SEDI Symposium, Whistler, Session 4: T h e Geodynamo and Observational Constraints, 1994b. Malkus, W.V.R., Energy sources for planetary dynamos, in Lectures on Solar and Planetary Dynamos, edited by M.R.E. Proctor and A.D. Gilbert, pp. 161-179, University Press, Cambridge U.K., 1994. Mathews, P.M., B.A. Buffett, T.A. Herring and LI. Shapiro, Forced nutations of the Earth: influence of inner core dynamics. 1. Theory, / . Geophys. Res., 96, 8219-8242, 1991a. Mathews, P.M., B.A. Buffett, T.A. Herring and I.t. Shapiro, Forced nutations of the Earth: influence of inner core dy­ namics. 2. Numerical results and comparisons, J. Geo­ phys. Res., 96, 8243-8257, 1991b. McFadden, P.L., and R.T. Merrill, Inhibition and geomag­ netic field reversals, / . Geophys. Res., 98, 6189-6199, 1993. Moffatt, H.K., Liquid metal MHD and the geodynamo, in Liquid Metal Magneto-hydrodynamics, edited by J. Lielpeteris and R. Moreau, pp. 403-412. Kluwer, Dordrecht, 1989. Moffatt, H.K., and D . E . Loper, T h e magnetostatic rise of a buoyant parcel in the Earth's core, Geophys. J. Int., 117, 394-402, 1994. Nakajima, T., and P.H. Roberts, A mapping method for solving dynamo equations, Proc. R. Soc. Lond., to ap­ pear, 1994a, Nakajima, T., and P.H. Roberts, An application of the map­ ping method t o the asymmetric dynamo equations, Ab­ stract, 4TH SEDI Symposium, Whistler, Session 4: T h e Geodynamo and Observational Constraints, 1994b. Nakajima, T., P.H. Roberts and M. Kono, An application of a mapping method t o the geodynamo problem. Abstract 01.01.02, 7th IAGA Scientific Assembly, Buenos Aires, 1993. 450 ROBERTS: DYNAMICS OF THE CORE,GEODYNAMO Olson, P., and G.A. Gla.tzmaier, Magnetoconvection in a Stewartson, K., On almost rigid rotations, 1. Fluid Mech., rotating sphere and the strncture of How in the liquid core, EOS, 74, Supp., 79, 1993 . Olson, P., and G.A. Gla.tzmaier, Magnetoconvection in a 16, 131-144, 1966. Sun, Z.-P., G. Schubert and G.A. Gla.tzmaier, Numerical simulations of thermal convection in a rapidly rotating rotating sphere and the structure of How in the liquid core, Ph"•. Earth planet. Int., to appear, 1994. Proctor, M .R.E., P.C. Matthews and A.M. Rucklidge (ed­ itors), Solar and PianetarrJ DflRa mo., University Press, Cambridge U.K., 1993. Proctor, M.R.E., and A.D. Gilbert (editors), Lecture, on So­ lar and PianetarrJ DflRamo" University Press,Cambridge U.K.,1994 . Roberts, P.H., A mapping method for integration o f the spherical shell cooled inhomogeneously from above, Geo­ ph"•. Altroph",. Fluid D"nam., to appear, 1994. Taylor, J . B .,The m agneto-hydrodynamics of a rotating fluid and the Earth'. dynamo problem, Proc. R. Soc. Lond., A!74, 274-283, 1963. U.S. Geodynamics Committee, Board on Earth Sciences and Resources, Commission on Geosciences and Resources, and National Research Council, The Geomagnetic Initia­ title, National Academy Pres,s Washington D.C.,1993. geodynamo equations, Abstract, 3rd SEDI Symposium, Vanyo, J.P., A geodynamo powered by Iuni-solar precesosi n, Mizusawa, 1992. Geophfl.l A.troph",. Fluid DflRam., 59, 2 09-234, 1991. Ruan,K. and D. Loper, On smal-l scale hydromagnetic flow in the Earth's core: motion of a rigid cylinder parallel to its axis, 1. Geomag. Geoelec., 45, 1539-1552,1993. Runcorn, S.K., and D.E. Winch,The Pacific cable experi­ ment, IAGA Program andAb,tractl, Vienna, p. 607,1991. Ruzmaikin, A., On the role of rotation of the inner core Vanyo, J.P., and D. Lods, MC B relative velocities in preces­ sionally drjven Hows, EOS, 75 Su pp., 120, 1994. Vanyo, J., core and Dso. liLdodcosraenidntPer.aW ctiioldnes,,PArbelcte ras ct,gsin 4tmh aSnEtlDe,IliSqyumid­ po.ium, Whi,Uer, Session 2: Core Mantle Interactions and the CMB, 1994b. rela.tive to the mantle, in Solar and PianetarrJ Dynamo" edited by M.R.E. Proctor, P.C. Matthews and A.M. Ruck­ lidge, pp . 265-270, University Press, Cambridge U.K.,19- Vanyo,J., P. Wilde, P. Cardin and P. Olson, Flow strnctures induced in a liquid core by precession, EOS, 79 Supp., 63, 1992. 9 3. Vanyo, J., P. Wilde, P. Cardin and P. Olson, Experim­ Ruzmaikin, A.A.,P.C. Liewer and J. Feynman,Random cell dynamo, Geophy•. A,trophy,. Fluid Dynam., 79, 163 - ents on precessing flows in the Earth's liquid core, Geo­ ph y•. 1. Int., to appear, 199480. 117, 1993. St. Pierre, M.G., The stability of the magnetostrophic ap­ proximation I: Taylor state solutions,Geophy,. A,trophy,. Fluid Dynam., 67, 99-115, 199380. Verhoogen, J.,Heat balance of the Ear�h's core,Geophy,. J., 4, 2 7 6-281, 1961. Wicht, J., and F .H . Busse,Dynamo action induced by lat­ eral variation of electrical conductivity,in Solar and Plan­ St . Pierre, M.G.,The strong field .branch of the Childres­s Soward dynamo,in Solar and Planetary Dynamo., edited etary Dynamo" edited by M.R.E. Proctor,P.C. Matthews and A.M. Rucklidge,pp. 329-337,University Press,Cam­ by M.R.E. Proctor,P.C . Matthews and A.M . Rucklidge, pp . 295-302, University Pres,s Cambridge U.K.,1993b. St. Pierre, M.G., On the local nature of fluid flow in the Earth's outer core,submitted,199480. bridge U.K., 1993. Wuildpet,oPfu.Dlly.,daenvdelJo.pPe.dVparneycoe,sFsiolonw,Aobb,sterravcatt,io4nths from spin­ SED] Sym­ po,ium, Whi,tler, Sesis on 2: Core Mantle Interactions and Sto. uPtieerrrceo,rM e,.AGb.,,tTrahcet,st4athbiSlitEyDoIfSbyumopyoa,nitum st,ruWcthuir,etslerin, Stehse­ the CMB,1994. sion 3: Core Dynamics and Thermodynamics, 1994b. St . Pierre,M.G., and P.H . local turbulence in the REoarbtehr'tss,cDoriere,cAt bsi,mtrualcatt, io4nths of the SEDI Sympo.ium, Whi,tler, Session 2: Core Mantle Interactions Paul H. Roberts, Institute of Geophysics and Planetary Physics,University of California,Los Angeles,CA90095. and the CMB, 1994. Stevenson,D .J., The core-mantle electrochemical cell as an explanation for part of the geomagnetic field, EOS, 7!, Supp., 133, 1991. (Received June 13, 1994; revised January 18, 1995; accepted January 23, 1995.)