Earth's core and the geodynamo

Earth's magnetic field is generated by fluid motion in the liquid iron core. Details of how this occurs are now emerging from numerical simulations that achieve a self-sustaining magnetic field. Early results predict a dominant dipole field outside the core, and some models even reproduce magnetic reversals. The simulations also show how different patterns of flow can produce similar external fields. Efforts to distinguish between the various possibilities appeal to observations of the time-dependent behavior of the field. Important constraints will come from geological records of the magnetic field in the past.     

Elastic anisotropy of Earth's inner core

Earth's solid-iron inner core is elastically anisotropic. Sound waves propagate faster along Earth's spin axis than in the equatorial plane. This anisotropy has previously been explained by a preferred orientation of the iron alloy hexagonal crystals. However, hexagonal iron becomes increasingly isotropic on increasing temperature at pressures of the inner core and is therefore unlikely to cause the anisotropy. An alternative explanation, supported by diamond anvil cell experiments, is that iron adopts a body-centered cubic form in the inner core. We show, by molecular dynamics simulations, that the body-centered cubic iron phase is extremely anisotropic to sound waves despite its high symmetry. Direct simulations of seismic wave propagation reveal an anisotropy of 12%, a value adequate to explain the anisotropy of the inner core.     

Body-centered cubic iron-nickel alloy in Earth's core

Cosmochemical, geochemical, and geophysical studies provide evidence that Earth's core contains iron with substantial (5 to 15%) amounts of nickel. The iron-nickel alloy Fe(0.9)Ni(0.1) has been studied in situ by means of angle-dispersive x-ray diffraction in internally heated diamond anvil cells (DACs), and its resistance has been measured as a function of pressure and temperature. At pressures above 225 gigapascals and temperatures over 3400 kelvin, Fe(0.9)Ni(0.1) adopts a body-centered cubic structure. Our experimental and theoretical results not only support the interpretation of shockwave data on pure iron as showing a solid-solid phase transition above about 200 gigapascals, but also suggest that iron alloys with geochemically reasonable compositions (that is, with substantial nickel, sulfur, or silicon content) adopt the bcc structure in Earth's inner core.     

The structure of iron in Earth's inner core

Earth's solid inner core is mainly composed of iron (Fe). Because the relevant ultrahigh pressure and temperature conditions are difficult to produce experimentally, the preferred crystal structure of Fe at the inner core remains uncertain. Static compression experiments showed that the hexagonal close-packed (hcp) structure of Fe is stable up to 377 gigapascals and 5700 kelvin, corresponding to inner core conditions. The observed weak temperature dependence of the c/a axial ratio suggests that hcp Fe is elastically anisotropic at core temperatures. Preferred orientation of the hcp phase may explain previously observed inner core seismic anisotropy.     

Sediments at the top of Earth's core

Unusual physical properties at the core-mantle boundary have been inferred from seismic and geodetic observations in recent years. We show how both types of observations can be explained by a layer of silicate sediments, which accumulate at the top of the core as Earth cools. Compaction of the sediments expels most of the liquid iron but leaves behind a small amount of core material, which is entrained in mantle convection and may account for the isotopic signatures of core material in some hot spot plumes. Extraction of light elements from the liquid core also enhances the vigor of convection in the core and may increase the power available to the geodynamo.     

Physics of iron at Earth's core conditions

The bulk properties of iron at the pressure and temperature conditions of Earth's core were determined by a method that combines first-principles and classical molecular dynamic simulations. The theory indicates that (i) the iron melting temperature at inner-core boundary (ICB) pressure (330 gigapascals) is 5400 (+/-400) kelvin; (ii) liquid iron at ICB conditions is about 6% denser than Earth's outer core; and (iii) the shear modulus of solid iron close to its melting line is 140 gigapascals, consistent with the seismic value for the inner core. These results reconcile melting temperature estimates based on sound velocity shock wave data with those based on diamond anvil cell experiments.     

Spin transition zone in Earth's lower mantle

Mineral properties in Earth's lower mantle are affected by iron electronic states, but representative pressures and temperatures have not yet been probed. Spin states of iron in lower-mantle ferropericlase have been measured up to 95 gigapascals and 2000 kelvin with x-ray emission in a laser-heated diamond cell. A gradual spin transition of iron occurs over a pressure-temperature range extending from about 1000 kilometers in depth and 1900 kelvin to 2200 kilometers and 2300 kelvin in the lower mantle. Because low-spin ferropericlase exhibits higher density and faster sound velocities relative to the high-spin ferropericlase, the observed increase in low-spin (Mg,Fe)O at mid-lower mantle conditions would manifest seismically as a lower-mantle spin transition zone characterized by a steeper-than-normal density gradient.     

On the presence of liquid in Earth's inner core

Seismological studies indicate that the inner core of Earth is anisotropic for compressional waves (P waves), and has low shear wave (S wave) velocity, and high seismic attenuation. Using an effective medium theory for composite materials, we show that the presence of a volume fraction of 3 to 10% liquid in the form of oblate spheroidal inclusions aligned in the equatorial plane between iron crystals is sufficient to explain the aforementioned seismic phenomena. Variation of S-wave velocity between the polar axis and equatorial plane is more sensitive to the addition of liquid than that of P waves. The liquid could arise from the presence of dendrites or a mixture of elements other than iron that exist in liquid form under inner-core conditions.     

Isotopic evidence of Cr partitioning into Earth's core

The distribution of chemical elements in primitive meteorites (chondrites), as building blocks of terrestrial planets, provides insight into the formation and early differentiation of Earth. The processes that resulted in the depletion of some elements [such as chromium (Cr)] in the bulk silicate Earth relative to chondrites, however, remain debated between leading candidate causes: volatility versus core partitioning. We show through high-precision measurements of Cr stable isotopes in a range of meteorites, which deviate by up to ~0.4 per mil from those of the bulk silicate Earth, that Cr depletion resulted from its partitioning into Earth's core, with a preferential enrichment in light isotopes. Ab initio calculations suggest that the isotopic signature was established at mid-mantle magma ocean depth as Earth accreted planetary embryos and progressively became more oxidized.     

Localized temporal change of the Earth's inner core boundary

Compressional waves of an earthquake doublet (two events occurring in the South Sandwich Islands on 1 December 1993 and 6 September 2003), recorded at three seismic stations in Russia and Kyrgyzstan and reflected off Earth's inner core boundary, arrived at least from 39 to 70 milliseconds earlier in the 2003 event than in the 1993 event. Such changes indicate that Earth's inner core radius enlarged locally beneath middle Africa by 0.98 to 1.75 kilometers between the times of these two events. Changes of the inner core radius may be explained by either a differential motion of the inner core, assuming that irregularities are present at the inner core boundary and fixed to the inner core, or a rapid growth of the inner core by this amount.     

Origin of the low rigidity of the Earth's inner core

Earth's solid-iron inner core has a low rigidity that manifests itself in the anomalously low velocities of shear waves as compared to shear wave velocities measured in iron alloys. Normally, when estimating the elastic properties of a polycrystal, one calculates an average over different orientations of a single crystal. This approach does not take into account the grain boundaries and defects that are likely to be abundant at high temperatures relevant for the inner core conditions. By using molecular dynamics simulations, we show that, if defects are considered, the calculated shear modulus and shear wave velocity decrease dramatically as compared to those estimates obtained from the averaged single-crystal values. Thus, the low shear wave velocity in the inner core is explained.     

Equatorially dominated magnetic field change at the surface of Earth's core

Slow temporal variations in Earth's magnetic field originate in the liquid outer core. We analyzed the evolution of nonaxisymmetric magnetic flux at the core surface over the past 400 years. We found that the most robust feature is westward motion at 17 kilometers per year, in a belt concentrated around the equator beneath the Atlantic hemisphere. Surprisingly, this motion is dominated by a single wavenumber and persists throughout the observation period. This phenomenon could be produced by an equatorial jet of core fluid, by hydromagnetic wave propagation, or by a combination of both. Discrimination between these mechanisms would provide useful constraints on the dynamics of Earth's core.     

Temperatures in Earth's Core Based on Melting and Phase Transformation Experiments on Iron

Experiments on melting and phase transformations on iron in a laser-heated, diamond-anvil cell to a pressure of 150 gigapascals (approximately 1.5 million atmospheres) show that iron melts at the central core pressure of 363.85 gigapascals at 6350 +/- 350 kelvin. The central core temperature corresponding to the upper temperature of iron melting is 6150 kelvin. The pressure dependence of iron melting temperature is such that a simple model can be used to explain the inner solid core and the outer liquid core. The inner core is nearly isothermal (6150 kelvin at the center to 6130 kelvin at the inner core-outer core boundary), is made of hexagonal closest-packed iron, and is about 1 percent solid (MgSiO(3) + MgO). By the inclusion of less than 2 percent of solid impurities with iron, the outer core densities along a thermal gradient (6130 kelvin at the base of the outer core and 4000 kelvin at the top) can be matched with the average seismic densities of the core.     

In Situ Determination of the NiAs Phase of FeO at High Pressure and Temperature

In situ synchrotron x-ray diffraction measurements of FeO at high pressures and high temperatures revealed that the high-pressure phase of FeO has the NiAs structure (B8). The lattice parameters of this NiAs phase at 96 gigapascals and 800 kelvin are a = 2.574(2) angstroms and c = 5.172(4) angstroms (the number in parentheses is the error in the last digit). Metallic behavior of the high-pressure phase is consistent with a covalently and metallically bonded NiAs structure of FeO. Transition to the NiAs structure of FeO would enhance oxygen solubility in molten iron. This transition thus provides a physiochemical basis for the incorporation of oxygen into the Earth's core.     

Experimental Evidence for a New Iron Phase and Implications for Earth's Core

Iron is known to occur in four different crystal structural forms. One of these, the densest form (epsilon phase, hexagonal close-packed) is considered to have formed Earth's core. Theoretical arguments based on available high-temperature and high-pressure iron data indicate the possibility of a fifth less dense iron phase forming the core. Study of iron phase transition conducted between pressures of 20 to 100 gigapascals and 1000 to 2200 Kelvin provides an experimental confirmation of the existence of this new phase. Thee epsilon iron phase transforms to this lower density phase before melting. The new phase may form a large part of Earth's core.     

Raman spectroscopy of iron to 152 gigapascals: implications for Earth's inner core

Raman spectra of hexagonal close-packed iron (varepsilon-Fe) have been measured from 15 to 152 gigapascals by using diamond-anvil cells with ultrapure synthetic diamond anvils. The results give a Gruneisen parameter gamma(0) = 1.68 (+/-0.20) and q = 0.7 (+/-0.5). Phenomenological modeling shows that the Raman-active mode can be approximately correlated with an acoustic phonon and thus provides direct information about the high-pressure elastic properties of iron, which have been controversial. In particular, the C(44) elastic modulus is found to be lower than previous determinations. This leads to changes of about 35% at core pressures for shear wave anisotropies.