Small-scale structure of the geodynamo inferred from Oersted and Magsat satellite data

The 'geodynamo' in the Earth's liquid outer core produces a magnetic field that dominates the large and medium length scales of the magnetic field observed at the Earth's surface. Here we use data from the currently operating Danish Oersted satellite, and from the US Magsat satellite that operated in 1979/80, to identify and interpret variations in the magnetic field over the past 20 years, down to length scales previously inaccessible. Projected down to the surface of the Earth's core, we found these variations to be small below the Pacific Ocean, and large at polar latitudes and in a region centred below southern Africa. The flow pattern at the surface of the core that we calculate to account for these changes is characterized by a westward flow concentrated in retrograde polar vortices and an asymmetric ring where prograde vortices are correlated with highs (and retrograde vortices with lows) in the historical (400-year average) magnetic field. This pattern is analogous to those seen in a large class of numerical dynamo simulations, except for its longitudinal asymmetry. If this asymmetric state was reached often in the past, it might account for several persistent patterns observed in the palaeomagnetic field. We postulate that it might also be a state in which the geodynamo operates before reversing.     

Power requirement of the geodynamo from ohmic losses in numerical and laboratory dynamos

In the Earth's fluid outer core, a dynamo process converts thermal and gravitational energy into magnetic energy. The power needed to sustain the geomagnetic field is set by the ohmic losses (dissipation due to electrical resistance). Recent estimates of ohmic losses cover a wide range, from 0.1 to 3.5 TW, or roughly 0.3-10% of the Earth's surface heat flow. The energy requirement of the dynamo puts constraints on the thermal budget and evolution of the core through Earth's history. Here we use a set of numerical dynamo models to derive scaling relations between the core's characteristic dissipation time and the core's magnetic and hydrodynamic Reynolds numbers--dimensionless numbers that measure the ratio of advective transport to magnetic and viscous diffusion, respectively. The ohmic dissipation of the Karlsruhe dynamo experiment supports a simple dependence on the magnetic Reynolds number alone, indicating that flow turbulence in the experiment and in the Earth's core has little influence on its characteristic dissipation time. We use these results to predict moderate ohmic dissipation in the range of 0.2-0.5 TW, which removes the need for strong radioactive heating in the core and allows the age of the solid inner core to exceed 2.5 billion years.     

A crystallizing dense magma ocean at the base of the Earth's mantle

The distribution of geochemical species in the Earth's interior is largely controlled by fractional melting and crystallization processes that are intimately linked to the thermal state and evolution of the mantle. The existence of patches of dense partial melt at the base of the Earth's mantle, together with estimates of melting temperatures for deep mantle phases and the amount of cooling of the underlying core required to maintain a geodynamo throughout much of the Earth's history, suggest that more extensive deep melting occurred in the past. Here we show that a stable layer of dense melt formed at the base of the mantle early in the Earth's history would have undergone slow fractional crystallization, and would be an ideal candidate for an unsampled geochemical reservoir hosting a variety of incompatible species (most notably the missing budget of heat-producing elements) for an initial basal magma ocean thickness of about 1,000 km. Differences in 142Nd/144Nd ratios between chondrites and terrestrial rocks can be explained by fractional crystallization with a decay timescale of the order of 1 Gyr. These combined constraints yield thermal evolution models in which radiogenic heat production and latent heat exchange prevent early cooling of the core and possibly delay the onset of the geodynamo to 3.4-4 Gyr ago.     

Structural and temporal requirements for geomagnetic field reversal deduced from lava flows

Reversals of the Earth's magnetic field reflect changes in the geodynamo--flow within the outer core--that generates the field. Constraining core processes or mantle properties that induce or modulate reversals requires knowing the timing and morphology of field changes that precede and accompany these reversals. But the short duration of transitional field states and fragmentary nature of even the best palaeomagnetic records make it difficult to provide a timeline for the reversal process. 40Ar/39Ar dating of lavas on Tahiti, long thought to record the primary part of the most recent 'Matuyama-Brunhes' reversal, gives an age of 795 +/- 7 kyr, indistinguishable from that of lavas in Chile and La Palma that record a transition in the Earth's magnetic field, but older than the accepted age for the reversal. Only the 'transitional' lavas on Maui and one from La Palma (dated at 776 +/- 2 kyr), agree with the astronomical age for the reversal. Here we propose that the older lavas record the onset of a geodynamo process, which only on occasion would result in polarity change. This initial instability, associated with the first of two decreases in field intensity, began approximately 18 kyr before the actual polarity switch. These data support the claim that complete reversals require a significant period for magnetic flux to escape from the solid inner core and sufficiently weaken its stabilizing effect.     

A deep dynamo generating Mercury's magnetic field

Mercury has a global magnetic field of internal origin and it is thought that a dynamo operating in the fluid part of Mercury's large iron core is the most probable cause. However, the low intensity of Mercury's magnetic field--about 1% the strength of the Earth's field--cannot be reconciled with an Earth-like dynamo. With the common assumption that Coriolis and Lorentz forces balance in planetary dynamos, a field thirty times stronger is expected. Here I present a numerical model of a dynamo driven by thermo-compositional convection associated with inner core solidification. The thermal gradient at the core-mantle boundary is subadiabatic, and hence the outer region of the liquid core is stably stratified with the dynamo operating only at depth, where a strong field is generated. Because of the planet's slow rotation the resulting magnetic field is dominated by small-scale components that fluctuate rapidly with time. The dynamo field diffuses through the stable conducting region, where rapidly varying parts are strongly attenuated by the skin effect, while the slowly varying dipole and quadrupole components pass to some degree. The model explains the observed structure and strength of Mercury's surface magnetic field and makes predictions that are testable with space missions both presently flying and planned.     

Thermochemical flows couple the Earth's inner core growth to mantle heterogeneity

Seismic waves sampling the top 100 km of the Earth's inner core reveal that the eastern hemisphere (40 degrees E-180 degrees E) is seismically faster, more isotropic and more attenuating than the western hemisphere. The origin of this hemispherical dichotomy is a challenging problem for our understanding of the Earth as a system of dynamically coupled layers. Previously, laboratory experiments have established that thermal control from the lower mantle can drastically affect fluid flow in the outer core, which in turn can induce textural heterogeneity on the inner core solidification front. The resulting texture should be consistent with other expected manifestations of thermal mantle control on the geodynamo, specifically magnetic flux concentrations in the time-average palaeomagnetic field over the past 5 Myr, and preferred eddy locations in flows imaged below the core-mantle boundary by the analysis of historical geomagnetic secular variation. Here we show that a single model of thermochemical convection and dynamo action can account for all these effects by producing a large-scale, long-term outer core flow that couples the heterogeneity of the inner core with that of the lower mantle. The main feature of this thermochemical 'wind' is a cyclonic circulation below Asia, which concentrates magnetic field on the core-mantle boundary at the observed location and locally agrees with core flow images. This wind also causes anomalously high rates of light element release in the eastern hemisphere of the inner core boundary, suggesting that lateral seismic anomalies at the top of the inner core result from mantle-induced variations in its freezing rate.     

Intense equatorial flux spots on the surface of the Earth's core

A large number of high-accuracy vector measurements of the Earth's magnetic field have recently become available from the satellite Oersted, complementing previous vector data from the satellite Magsat, which operated in 1979/80. These data can be used to infer the morphology of the magnetic field at the surface of the fluid core, approximately 2,900 km below the Earth's surface. Here I apply a new methodology to these data to calculate maps of the magnetic field at the core surface which show intense flux spots in equatorial regions. The intensity of these features is unusually large--some have intensities comparable to high-latitude flux patches near the poles, previously identified as the major component of the dynamo field. The tendency for pairing of some of these spots to the north and south of the geographical equator suggests they might be associated with the tops of equatorially symmetric columnar structures in the fluid, or their antisymmetric equivalents. The drift of the equatorial features may represent material flow or could represent wave motion; discrimination of these two effects based on future data could provide new information on the strength of the hidden toroidal magnetic field of the Earth.     

Dependence of the duration of geomagnetic polarity reversals on site latitude

An important constraint on the processes governing the geodynamo--the flow in the outer core responsible for generating Earth's magnetic field--is the duration of geomagnetic polarity reversals; that is, how long it takes for Earth's magnetic field to reverse. It is generally accepted that Earth's magnetic field strength drops to low levels during polarity reversals, and the field direction progresses through a 180 degrees change while the field is weak. The time it takes for this process to happen, however, remains uncertain, with estimates ranging from a few thousand up to 28,000 years. Here I present an analysis of the available sediment records of the four most recent polarity reversals. These records yield an average estimate of about 7,000 years for the time it takes for the directional change to occur. The variation about this mean duration is not random, but instead varies with site latitude, with shorter durations observed at low-latitude sites, and longer durations observed at mid- to high-latitude sites. Such variation of duration with site latitude is predicted by simple geometrical reversal models, in which non-dipole fields are allowed to persist while the axial dipole decays through zero and then builds in the opposite direction, and provides a constraint on numerical dynamo models.     

Melting of the Earth's inner core

The Earth's magnetic field is generated by a dynamo in the liquid iron core, which convects in response to cooling of the overlying rocky mantle. The core freezes from the innermost surface outward, growing the solid inner core and releasing light elements that drive compositional convection. Mantle convection extracts heat from the core at a rate that has enormous lateral variations. Here we use geodynamo simulations to show that these variations are transferred to the inner-core boundary and can be large enough to cause heat to flow into the inner core. If this were to occur in the Earth, it would cause localized melting. Melting releases heavy liquid that could form the variable-composition layer suggested by an anomaly in seismic velocity in the 150 kilometres immediately above the inner-core boundary. This provides a very simple explanation of the existence of this layer, which otherwise requires additional assumptions such as locking of the inner core to the mantle, translation from its geopotential centre or convection with temperature equal to the solidus but with composition varying from the outer to the inner core. The predominantly narrow downwellings associated with freezing and broad upwellings associated with melting mean that the area of melting could be quite large despite the average dominance of freezing necessary to keep the dynamo going. Localized melting and freezing also provides a strong mechanism for creating seismic anomalies in the inner core itself, much stronger than the effects of variations in heat flow so far considered.     

Formation of current coils in geodynamo simulations

Computer simulations have been playing an important role in the development of our understanding of the geodynamo, but direct numerical simulation of the geodynamo with a realistic parameter regime is still beyond the power of today's supercomputers. Difficulties in simulating the geodynamo arise from the extreme conditions of the core, which are characterized by very large or very small values of the non-dimensional parameters of the system. Among them, the Ekman number, E, has been adopted as a barometer of the distance of simulations from real core conditions, in which E is of the order of 10(-15). Following the initial computer simulations of the geodynamo, the Ekman number achieved has been steadily decreasing, with recent geodynamo simulations performed with E of the order of 10(-6). Here we present a geodynamo simulation with an Ekman number of the order of 10(-7)-the highest-resolution simulation yet achieved, making use of 4,096 processors of the Earth Simulator. We have found that both the convection flow and magnetic field structures are qualitatively different from those found in larger-Ekman-number dynamos. The convection takes the form of sheet plumes or radial sheet jets, rather than the columnar cell structures that are usually found. We have found that this sheet plume convection is an effective dynamo and the generated current is organized as a set of coils in the shape of helical springs or at times as a torus.     

Experimental earth tidal models in considering nearly diurnal free wobble of the Earth''''s liquid core

Based on the 28 series of the high precision and high minute sampling tidal gravity observations at 20 stations in Global Geodynamics Project (GGP) network, the resonant parameters of the Earth's nearly diurnal free wobble (including the eigenperiods, resonant strengths and quality factots) are precisely determined. The discrepancy of the eigenperiod between observed and theoretical values is studied, the important conclusion that the real dynamic ellipticity of the liquid core is about 5% larger than the one under the static equilibrium assumption is approved by using our gravity technique. The experimental Earth's tidal gravity models with considering the nearly diurnal free wobble of the Earth's liquid core are constructed in this study. The numerical results show that the difference among three experimental models is less than 0.1%, and the largest discrepancy compared to those widely used nowdays given by Dehant (1999) and Mathews (2001) is only about 0.4%. It can provide with the most recent real     

Experimental earth tidal models in considering nearly diurnal free wobble of the Earth＇s liquid core

Based on the 28 series of the high precision and high minute sampling tidal gravity observations at 20 stations in Global Geodynamics Project (GGP) network, the resonant parameters of the Earth‘s nearly diurnal free wobble (including the eigenperiods, resonant strengths and quality factots) are precisely determined. The discrepancy of the eigenperiod between observed and theoretical values is studied, the important conclusion that the real dynamic ellipticity of the liquid core is about 5% larger than the one under the static equilibrium assumption is approved by using our gravity technique. The experimental Earth‘s tidal gravity models with considering the nearly diurnal free wobble of the Earth‘s liquid core are constructed in this study. The numerical results show that the difference among three experimental models is less than 0.1%, and the largest discrepancy compared to those widely used nowdays given by Dehant (1999) and Mathews (2001) is only about 0.4%. It can provide with the most recent real experimental tidal gravity models for the global study of the Earth‘s tides, geodesy and space techniques and so on.     

Evidence for enhanced mixing over rough topography in the abyssal ocean

The overturning circulation of the ocean plays an important role in modulating the Earth's climate. But whereas the mechanisms for the vertical transport of water into the deep ocean--deep water formation at high latitudes--and horizontal transport in ocean currents have been largely identified, it is not clear how the compensating vertical transport of water from the depths to the surface is accomplished. Turbulent mixing across surfaces of constant density is the only viable mechanism for reducing the density of the water and enabling it to rise. However, measurements of the internal wave field, the main source of energy for mixing, and of turbulent dissipation rates, have typically implied diffusivities across surfaces of equal density of only approximately 0.1 cm2 s(-1), too small to account for the return flow. Here we report measurements of tracer dispersion and turbulent energy dissipation in the Brazil basin that reveal diffusivities of 2-4 cm2 s(-1) at a depth of 500 m above abyssal hills on the flank of the Mid-Atlantic Ridge, and approximately 10 cm2 s(-1) nearer the bottom. This amount of mixing, probably driven by breaking internal waves that are generated by tidal currents flowing over the rough bathymetry, may be large enough to close the buoyancy budget for the Brazil basin and suggests a mechanism for closing the global overturning circulation.     

Evidence for an oxygen-depleted liquid outer core of the Earth

On the basis of geophysical observations, cosmochemical constraints, and high-pressure experimental data, the Earth's liquid outer core consists of mainly liquid iron alloyed with about ten per cent (by weight) of light elements. Although the concentrations of the light elements are small, they nevertheless affect the Earth's core: its rate of cooling, the growth of the inner core, the dynamics of core convection, and the evolution of the geodynamo. Several light elements-including sulphur, oxygen, silicon, carbon and hydrogen-have been suggested, but the precise identity of the light elements in the Earth's core is still unclear. Oxygen has been proposed as a major light element in the core on the basis of cosmochemical arguments and chemical reactions during accretion. Its presence in the core has direct implications for Earth accretion conditions of oxidation state, pressure and temperature. Here we report new shockwave data in the Fe-S-O system that are directly applicable to the outer core. The data include both density and sound velocity measurements, which we compare with the observed density and velocity profiles of the liquid outer core. The results show that we can rule out oxygen as a major light element in the liquid outer core because adding oxygen into liquid iron would not reproduce simultaneously the observed density and sound velocity profiles of the outer core. An oxygen-depleted core would imply a more reduced environment during early Earth accretion.     

Thermal and electrical conductivity of iron at Earth's core conditions

The Earth acts as a gigantic heat engine driven by the decay of radiogenic isotopes and slow cooling, which gives rise to plate tectonics, volcanoes and mountain building. Another key product is the geomagnetic field, generated in the liquid iron core by a dynamo running on heat released by cooling and freezing (as the solid inner core grows), and on chemical convection (due to light elements expelled from the liquid on freezing). The power supplied to the geodynamo, measured by the heat flux across the core-mantle boundary (CMB), places constraints on Earth's evolution. Estimates of CMB heat flux depend on properties of iron mixtures under the extreme pressure and temperature conditions in the core, most critically on the thermal and electrical conductivities. These quantities remain poorly known because of inherent experimental and theoretical difficulties. Here we use density functional theory to compute these conductivities in liquid iron mixtures at core conditions from first principles--unlike previous estimates, which relied on extrapolations. The mixtures of iron, oxygen, sulphur and silicon are taken from earlier work and fit the seismologically determined core density and inner-core boundary density jump. We find both conductivities to be two to three times higher than estimates in current use. The changes are so large that core thermal histories and power requirements need to be reassessed. New estimates indicate that the adiabatic heat flux is 15 to 16 terawatts at the CMB, higher than present estimates of CMB heat flux based on mantle convection; the top of the core must be thermally stratified and any convection in the upper core must be driven by chemical convection against the adverse thermal buoyancy or lateral variations in CMB heat flow. Power for the geodynamo is greatly restricted, and future models of mantle evolution will need to incorporate a high CMB heat flux and explain the recent formation of the inner core.     

Geomagnetic intensity variations over the past 780 kyr obtained from near-seafloor magnetic anomalies

Knowledge of past variations in the intensity of the Earth's magnetic field provides an important constraint on models of the geodynamo. A record of absolute palaeointensity for the past 50 kyr has been compiled from archaeomagnetic and volcanic materials, and relative palaeointensities over the past 800 kyr have been obtained from sedimentary sequences. But a long-term record of geomagnetic intensity should also be carried by the thermoremanence of the oceanic crust Here we show that near-seafloor magnetic anomalies recorded over the southern East Pacific Rise are well correlated with independent estimates of geomagnetic intensity during the past 780 kyr. Moreover, the pattern of absolute palaeointensity of seafloor glass samples from the same area agrees with the well-documented dipole intensity pattern for the past 50 kyr. A comparison of palaeointensities derived from seafloor glass samples with global intensity variations thus allows us to estimate the ages of surficial lava flows in this region. The record of geomagnetic intensity preserved in the oceanic crust should provide a higher-time-resolution record of crustal accretion processes at mid-ocean ridges than has previously been obtainable.     

The origin of geomagnetic jerks

Geomagnetic jerks, which in the second half of the twentieth century occurred in 1969 (refs 1, 2), 1978 (refs 3, 4), 1991 (ref. 5) and 1999 (ref. 6), are abrupt changes in the second time-derivative (secular acceleration) of the Earth's magnetic field. Jerks separate periods of almost steady secular acceleration, so that the first time-derivative (secular variation) appears as a series of straight-line segments separated by geomagnetic jerks. The fact that they represent a reorganization of the secular variation implies that they are of internal origin (as has been established through spherical harmonic analysis), and their short timescale implies that they are due to a change in the fluid flow at the surface of the Earth's core (as has also been established through mapping the time-varying flow at the core surface). However, little is understood of their physical origin. Here we show that geomagnetic jerks can be explained by the combination of a steady flow and a simple time-varying, axisymmetric, equatorially symmetric, toroidal zonal flow. Such a flow is consistent with torsional oscillations in the Earth's core, which are simple oscillatory flows in the core that are expected on theoretical grounds, and observed in both core flow models and numerical dynamo models.     

Proterozoic low orbital obliquity and axial-dipolar geomagnetic field from evaporite palaeolatitudes

Palaeomagnetism of climatically sensitive sedimentary rock types, such as glacial deposits and evaporites, can test the uniformitarianism of ancient geomagnetic fields and palaeoclimate zones. Proterozoic glacial deposits laid down in near-equatorial palaeomagnetic latitudes can be explained by 'snowball Earth' episodes, high orbital obliquity or markedly non-uniformitarian geomagnetic fields. Here I present a global palaeomagnetic compilation of the Earth's entire basin-scale evaporite record. Magnetic inclinations are consistent with low orbital obliquity and a geocentric-axial-dipole magnetic field for most of the past two billion years, and the snowball Earth hypothesis accordingly remains the most viable model for low-latitude Proterozoic ice ages. Efforts to reconstruct Proterozoic supercontinents are strengthened by this demonstration of a consistently axial and dipolar geomagnetic reference frame, which itself implies stability of geodynamo processes on billion-year timescales.     

Constraints on the composition of the Earth's core from ab initio calculations

Knowledge of the composition of the Earth's core is important for understanding its melting point and therefore the temperature at the inner-core boundary and the temperature profile of the core and mantle. In addition, the partitioning of light elements between solid and liquid, as the outer core freezes at the inner-core boundary, is believed to drive compositional convection, which in turn generates the Earth's magnetic field. It is generally accepted that the liquid outer core and the solid inner core consist mainly of iron. The outer core, however, is also thought to contain a significant fraction of light elements, because its density--as deduced from seismological data and other measurements--is 6-10 per cent less than that estimated for pure liquid iron. Similar evidence indicates a smaller but still appreciable fraction of light elements in the inner core. The leading candidates for the light elements present in the core are sulphur, oxygen and silicon. Here we obtain a constraint on core composition derived from ab initio calculation of the chemical potentials of light elements dissolved in solid and liquid iron. We present results for the case of sulphur, which provide strong evidence against the proposal that the outer core is close to being a binary iron-sulphur mixture.