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.     

Slow differential rotation of the Earth's inner core indicated by temporal changes in scattering

The finding that the Earth's inner core might be rotating faster than the mantle has important implications for our understanding of core processes, including the generation of the Earth's magnetic field. But the reported signal is subtle--a change of about 0.01 s per year in the separation of two seismic waves with differing paths through the core. Subsequent studies of such data have generally supported the conclusion that differential rotation exists, but the difficulty of accurately locating historic earthquakes and possible biases induced by strong lateral variations in structure near the core-mantle boundary have raised doubt regarding the proposed inner-core motion. Also, a study of free oscillations constrained the motion to be relatively small compared to previous estimates and it has been proposed that the interaction of inner-core boundary topography and mantle heterogeneity might lock the inner core to the mantle. The recent detection of seismic waves scattered in the inner core suggests a simple test of inner-core motion. Here we compare scattered waves recorded in Montana, USA, from two closely located nuclear tests at Novaya Zemlya, USSR, in 1971 and 1974. The data show small but coherent changes in scattering which point toward an inner-core differential rotation rate of 0.15 degrees per year--consistent with constraints imposed by the free-oscillation data.     

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.     

An impact-driven dynamo for the early Moon

The origin of lunar magnetic anomalies remains unresolved after their discovery more than four decades ago. A commonly invoked hypothesis is that the Moon might once have possessed a thermally driven core dynamo, but this theory is problematical given the small size of the core and the required surface magnetic field strengths. An alternative hypothesis is that impact events might have amplified ambient fields near the antipodes of the largest basins, but many magnetic anomalies exist that are not associated with basin antipodes. Here we propose a new model for magnetic field generation, in which dynamo action comes from impact-induced changes in the Moon's rotation rate. Basin-forming impact events are energetic enough to have unlocked the Moon from synchronous rotation, and we demonstrate that the subsequent large-scale fluid flows in the core, excited by the tidal distortion of the core-mantle boundary, could have powered a lunar dynamo. Predicted surface magnetic field strengths are on the order of several microteslas, consistent with palaeomagnetic measurements, and the duration of these fields is sufficient to explain the central magnetic anomalies associated with several large impact basins.     

Chemical interaction of Fe and Al(2)O3 as a source of heterogeneity at the Earth's core-mantle boundary

Seismological studies have revealed that a complex texture or heterogeneity exists in the Earth's inner core and at the boundary between core and mantle. These studies highlight the importance of understanding the properties of iron when modelling the composition and dynamics of the core and the interaction of the core with the lowermost mantle. One of the main problems in inferring the composition of the lowermost mantle is our lack of knowledge of the high-pressure and high-temperature chemical reactions that occur between iron and the complex Mg-Fe-Si-Al-oxides which are thought to form the bulk of the Earth's lower mantle. A number of studies have demonstrated that iron can react with MgSiO3-perovskite at high pressures and high temperatures, and it was proposed that the chemical nature of this process involves the reduction of silicon by the more electropositive iron. Here we present a study of the interaction between iron and corundum (Al(2)O3) in electrically- and laser-heated diamond anvil cells at 2,000-2,200 K and pressures up to 70 GPa, simulating conditions in the Earth's deep interior. We found that at pressures above 60 GPa and temperatures of 2,200 K, iron and corundum react to form iron oxide and an iron-aluminium alloy. Our results demonstrate that iron is able to reduce aluminium out of oxides at core-mantle boundary conditions, which could provide an additional source of light elements in the Earth's core and produce significant heterogeneity at the core-mantle boundary.     

Magnetic reversals and mass extinctions

Previous analyses of the time distribution of reversals of the Earth's magnetic field have yielded mixed results. Some authors have claimed significant periodicities of order 10(7) yr whereas others have reported failure to reject null hypothesis of random spacing at that scale. Because of repeated suggestions that field reversal is linked to biological extinction, further analysis of the magnetic times series is appropriate. I present here the results of a study of the reversal record of the past 165 Myr. A stationary periodicity of 30 Myr emerges (superimposed on the non-stationarities already established by others), which predicts pulses of increased reversal activity centered at 10, 40, 70, ... Myr BP.     

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.     

Geochemistry: does U-Pb date Earth's core formation?

Constraining the timing of the formation of Earth's core, which defines the birth of our planet, is essential for understanding the early evolution of Earth-like planets. Wood and Halliday and Halliday discuss the apparent discrepancy between the U-Pb (60-80 Myr) and Hf-W clocks (30 Myr) in determining the timescale of Earth's accretion and core formation. We find that the information the authors present is at times contradictory (for example, compare Fig. 1 in ref. 1 with Fig. 1 in ref. 2) and confusing and could suggest that the U-Pb clock constrains core formation better than the Hf-W system. Here we point out the limitations of the U-Pb system and show that the U-Pb age cannot be used to argue for protracted accretion and/or core formation (>50 Myr) because this clock only records the processes that occurred during the last 1% of Earth's accretion and core formation in the Wood and Halliday mechanism.     

Subducted banded iron formations as a source of ultralow-velocity zones at the core-mantle boundary

Ultralow-velocity zones (ULVZs) are regions of the Earth's core-mantle boundary about 1-10 kilometres thick exhibiting seismic velocities that are lower than radial-Earth reference models by about 10-20 per cent for compressional waves and 10-30 per cent for shear waves. It is also thought that such regions have an increased density of about 0-20 per cent (ref. 1). A number of origins for ULVZs have been proposed, such as ponding of dense silicate melt, core-mantle reaction zones or underside sedimentation from the core. Here we suggest that ULVZs might instead be relics of banded iron formations subducted to the core-mantle boundary between 2.8 and 1.8 billion years ago. Consisting mainly of interbedded iron oxides and silica, such banded iron formations were deposited in the world's oceans during the late Archaean and early Proterozoic eras. We argue that these layers, as part of the ocean floor, would be recycled into the Earth's interior by subduction, sink to the bottom of the mantle and may explain all of the observed features of ULVZs.     

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.     

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.     

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.     

Tidal dissipation and the strength of the Earth's internal magnetic field

Magnetic fields at the Earth's surface represent only a fraction of the field inside the core. The strength and structure of the internal field are poorly known, yet the details are important for our understanding of the geodynamo. Here I obtain an indirect estimate for the field strength from measurements of tidal dissipation. Tidally driven flow in the Earth's liquid core develops internal shear layers, which distort the internal magnetic field and generate electric currents. Ohmic losses damp the tidal motions and produce detectable signatures in the Earth's nutations. Previously reported evidence of anomalous dissipation in nutations can be explained with a core-averaged field of 2.5?mT, eliminating the need for high fluid viscosity or a stronger magnetic field at the inner-core boundary. Estimates for the internal field constrain the power required for the geodynamo.     

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.     

Experimental evidence that potassium is a substantial radioactive heat source in planetary cores

The hypothesis that (40)K may be a significant radioactive heat source in the Earth's core was proposed on theoretical grounds over three decades ago, but experiments have provided only ambiguous and contradictory evidence for the solubility of potassium in iron-rich alloys. The existence of such radioactive heat in the core would have important implications for our understanding of the thermal evolution of the Earth and global processes such as the generation of the geomagnetic field, the core-mantle boundary heat flux and the time of formation of the inner core. Here we provide experimental evidence to show that the ambiguous results obtained from earlier experiments are probably due to previously unrecognized experimental and analytical difficulties. The high-pressure, high-temperature data presented here show conclusively that potassium enters iron sulphide melts in a strongly temperature-dependent fashion and that (40)K can serve as a substantial heat source in the cores of the Earth and Mars.     

Iron-silica interaction at extreme conditions and the electrically conducting layer at the base of Earth's mantle

The boundary between the Earth's metallic core and its silicate mantle is characterized by strong lateral heterogeneity and sharp changes in density, seismic wave velocities, electrical conductivity and chemical composition. To investigate the composition and properties of the lowermost mantle, an understanding of the chemical reactions that take place between liquid iron and the complex Mg-Fe-Si-Al-oxides of the Earth's lower mantle is first required. Here we present a study of the interaction between iron and silica (SiO2) in electrically and laser-heated diamond anvil cells. In a multianvil apparatus at pressures up to 140 GPa and temperatures over 3,800 K we simulate conditions down to the core-mantle boundary. At high temperature and pressures below 40 GPa, iron and silica react to form iron oxide and an iron-silicon alloy, with up to 5 wt% silicon. At pressures of 85-140 GPa, however, iron and SiO2 do not react and iron-silicon alloys dissociate into almost pure iron and a CsCl-structured (B2) FeSi compound. Our experiments suggest that a metallic silicon-rich B2 phase, produced at the core-mantle boundary (owing to reactions between iron and silicate), could accumulate at the boundary between the mantle and core and explain the anomalously high electrical conductivity of this region.     

Efficacy of the post-perovskite phase as an explanation for lowermost-mantle seismic properties

Constraining the chemical, rheological and electromagnetic properties of the lowermost mantle (D') is important to understand the formation and dynamics of the Earth's mantle and core. To explain the origin of the variety of characteristics of this layer observed with seismology, a number of theories have been proposed, including core-mantle interaction, the presence of remnants of subducted material and that D' is the site of a mineral phase transformation. This final possibility has been rejuvenated by recent evidence for a phase change in MgSiO3 perovskite (thought to be the most prevalent phase in the lower mantle) at near core-mantle boundary temperature and pressure conditions. Here we explore the efficacy of this 'post-perovskite' phase to explain the seismic properties of the lowermost mantle through coupled ab initio and seismic modelling of perovskite and post-perovskite polymorphs of MgSiO3, performed at lowermost-mantle temperatures and pressures. We show that a post-perovskite model can explain the topography and location of the D' discontinuity, apparent differences in compressional- and shear-wave models and the observation of a deeper, weaker discontinuity. Furthermore, our calculations show that the regional variations in lower-mantle shear-wave anisotropy are consistent with the proposed phase change in MgSiO3 perovskite.     

Fine-scale heterogeneity in the Earth's inner core

The seismological properties of the Earth's inner core have become of particular interest as we understand more about its composition and thermal state. Observations of anisotropy and velocity heterogeneity in the inner core are beginning to reveal how it has grown and whether it convects. The attenuation of seismic waves in the inner core is strong, and studies of seismic body waves have found that this high attenuation is consistent with either scattering or intrinsic attenuation. The outermost portion of the inner core has been inferred to possess layering and to be less anisotropic than at greater depths. Here we present observations of seismic waves scattered in the inner core which follow the expected arrival time of the body-wave reflection from the inner-core boundary. The amplitude of these scattered waves can be explained by stiffness variations of 1.2% with a scale length of 2 kilometres across the outermost 300 km of the inner core. These variations might be caused by variations in composition, by pods of partial melt in a mostly solid matrix or by variations in the orientation or strength of seismic anisotropy.     

Seismological constraints on a possible plume root at the core-mantle boundary

Recent seismological discoveries have indicated that the Earth's core-mantle boundary is far more complex than a simple boundary between the molten outer core and the silicate mantle. Instead, its structural complexities probably rival those of the Earth's crust. Some regions of the lowermost mantle have been observed to have seismic wave speed reductions of at least 10 per cent, which appear not to be global in extent. Here we present robust evidence for an 8.5-km-thick and approximately 50-km-wide pocket of dense, partially molten material at the core-mantle boundary east of Australia. Array analyses of an anomalous precursor to the reflected seismic wave ScP reveal compressional and shear-wave velocity reductions of 8 and 25 per cent, respectively, and a 10 per cent increase in density of the partially molten aggregate. Seismological data are incompatible with a basal layer composed of pure melt, and thus require a mechanism to prevent downward percolation of dense melt within the layer. This may be possible by trapping of melt by cumulus crystal growth following melt drainage from an anomalously hot overlying region of the lowermost mantle. This magmatic evolution and the resulting cumulate structure seem to be associated with overlying thermal instabilities, and thus may mark a root zone of an upwelling plume.     

Hemispherical variations in seismic velocity at the top of the Earth's inner core

Knowledge of the seismic velocity structure at the top of the Earth's inner core is important for deciphering the physical processes responsible for inner-core growth. Previous global seismic studies have focused on structures found 100 km or deeper within the inner core, with results for the uppermost 100 km available for only isolated regions. Here we present constraints on seismic velocity variations just beneath the inner-core boundary, determined from the difference in travel time between waves reflected at the inner-core boundary and those transmitted through the inner core. We found that these travel-time residuals-observed on both global seismograph stations and several regional seismic networks-are systematically larger, by about 0.8 s, for waves that sample the 'eastern hemisphere' of the inner core (40 degrees E to 180 degrees E) compared to those that sample the 'western hemisphere' (180 degrees W to 40 degrees E). These residuals show no correlation with the angle at which the waves traverse the inner core; this indicates that seismic anisotropy is not strong in this region and that the isotropic seismic velocity of the eastern hemisphere is about 0.8% higher than that of the western hemisphere.