Outer-core compositional stratification from observed core wave speed profiles

Extensive iron-water exchange at Earth's core-mantle boundary can explain seismic anomalies

Seismological observations indicate the presence of chemical heterogeneities at the lowermost mantle, just above the core-mantle boundary (CMB), sparking debate over their origins. A plausible explanation for the enigmatic seismic wave velocities observed in ultra-low-velocity zones (ULVZs) is the process of iron enrichment from the core to the silicate mantle. However, traditional models based on diffusion of atoms and penetration of molten iron fail to account for the significant iron enrichment observed in ULVZs. Here, we show that the chemical reaction between silicate bridgmanite and iron under hydrous conditions leads to profound iron enrichment within silicate, a process not seen in anhydrous conditions. Our findings suggest that the interaction between the core and mantle facilitates deep iron enrichment over a few kilometres at the bottom of the mantle when water is present. We propose that the seismic signatures observed in ULVZs indicate whole mantle convection, accompanied by deep water cycles from the crust to the core through Earth's history.

Seismic low-velocity equatorial torus in the Earth's outer core: Evidence from the late-coda correlation wavefield

Thermochemical inhomogeneities in the Earth's outer core that enhance our understanding of the geodynamo have been elusive. Seismic constraints on such inhomogeneities would provide clues on the amount and distribution of light elements in the core apart from iron and nickel. Here, we present evidence for a low-velocity volume within the outer core via the global coda correlation wavefield. Several key correlogram features with a unique sensitivity to the liquid core show variations with wave paths remarkably slower in the equatorial than polar planes. We constrain a torus structure at low latitudes with ~2% lower velocity than the surrounding liquid outer core via waveform modeling. We propose a thermochemical origin for such a low-velocity torus, providing important constraints on the dynamical processes of the Earth's outer core.

Evidence for a liquid silicate layer atop the Martian core

Seismic recordings made during the InSight mission1 suggested that Mars's liquid core would need to be approximately 27% lighter than pure liquid iron2,3, implying a considerable complement of light elements. Core compositions based on seismic and bulk geophysical constraints, however, require larger quantities of the volatile elements hydrogen, carbon and sulfur than those that were cosmochemically available in the likely building blocks of Mars4. Here we show that multiply diffracted P waves along a stratified core-mantle boundary region of Mars in combination with first-principles computations of the thermoelastic properties of liquid iron-rich alloys3 require the presence of a fully molten silicate layer overlying a smaller, denser liquid core. Inverting differential body wave travel time data with particular sensitivity to the core-mantle boundary region suggests a decreased core radius of 1,675 ± 30 km associated with an increased density of 6.65 ± 0.1 g cm-3, relative to previous models2,4-8, while the thickness and density of the molten silicate layer are 150 ± 15 km and 4.05 ± 0.05 g cm-3, respectively. The core properties inferred here reconcile bulk geophysical and cosmochemical requirements, consistent with a core containing 85-91 wt% iron-nickel and 9-15 wt% light elements, chiefly sulfur, carbon, oxygen and hydrogen. The chemical characteristics of a molten silicate layer above the core may be revealed by products of Martian magmatism.

Satellite magnetic data reveal interannual waves in Earth's core

SignificanceThe physics responsible for most of the interannual geomagnetic field changes, continually recorded by satellites for 20 years, is a long-standing open issue. By analyzing magnetic data, we detect Magneto-Coriolis waves in the Earth's outer core that account for a significant part of this signal. We further propose theoretical advances in the physical characterization of these waves, enabling a deeper understanding of the dynamics behind the geomagnetic signal. It should allow one to better sketch the heterogeneous magnetic field deep within the core, shedding further light on the mechanisms that sustain the geodynamo. Our interpretation does not require the presence of a stratified layer at the top of the core, with potent consequences regarding the Earth's thermal history.

Stratification in planetary cores by liquid immiscibility in Fe-S-H

Liquid-liquid immiscibility has been widely observed in iron alloy systems at ambient pressure and is important for the structure and dynamics in iron cores of rocky planets. While such previously known liquid immiscibility has been demonstrated to disappear at relatively low pressures, here we report immiscible S(±Si,O)-rich liquid and H(±C)-rich liquid above ~20 GPa, corresponding to conditions of the Martian core. Mars’ cosmochemically estimated core composition is likely in the miscibility gap, and the separation of two immiscible liquids could have driven core convection and stable stratification, which explains the formation and termination of the Martian planetary magnetic field. In addition, we observed liquid immiscibility in Fe-S-H(±Si,O,C) at least to 118 GPa, suggesting that it can occur in the Earth’s topmost outer core and form a low-velocity layer below the core-mantle boundary.

Yokoo et al. find the liquid immiscibility between H-rich and S-rich liquids Fe above 20 GPa. The separation of immiscible liquids could explain the disappearance of Mars’ magnetic field and the formation of low-velocity layer atop the Earth’s core.

Thermal conductivity of Fe-Si alloys and thermal stratification in Earth's core

Light elements in Earth's core play a key role in driving convection and influencing geodynamics, both of which are crucial to the geodynamo. However, the thermal transport properties of iron alloys at high-pressure and -temperature conditions remain uncertain. Here we investigate the transport properties of solid hexagonal close-packed and liquid Fe-Si alloys with 4.3 and 9.0 wt % Si at high pressure and temperature using laser-heated diamond anvil cell experiments and first-principles molecular dynamics and dynamical mean field theory calculations. In contrast to the case of Fe, Si impurity scattering gradually dominates the total scattering in Fe-Si alloys with increasing Si concentration, leading to temperature independence of the resistivity and less electron-electron contribution to the conductivity in Fe-9Si. Our results show a thermal conductivity of ∼100 to 110 W⋅m^-1⋅K^-1 for liquid Fe-9Si near the topmost outer core. If Earth's core consists of a large amount of silicon (e.g., > 4.3 wt %) with such a high thermal conductivity, a subadiabatic heat flow across the core-mantle boundary is likely, leaving a 400- to 500-km-deep thermally stratified layer below the core-mantle boundary, and challenges proposed thermal convection in Fe-Si liquid outer core.

Core Eigenmodes and their Impact on the Earth's Rotation

Changes in the Earth's rotation are deeply connected to fluid dynamical processes in the outer core. This connection can be explored by studying the associated Earth eigenmodes with periods ranging from nearly diurnal to multi-decadal. It is essential to understand how the rotational and fluid core eigenmodes mutually interact, as well as their dependence on a host of diverse factors, such as magnetic effects, density stratification, fluid instabilities or turbulence. It is feasible to build detailed models including many of these features, and doing so will in turn allow us to extract more (indirect) information about the Earth's interior. In this article, we present a review of some of the current models, the numerical techniques, their advantages and limitations and the challenges on the road ahead.

Gravity Variations and Ground Deformations Resulting from Core Dynamics

ABSTRACT

Fluid motion within the Earth's liquid outer core leads to internal mass redistribution. This occurs through the advection of density anomalies within the volume of the liquid core and by deformation of the solid boundaries of the mantle and inner core which feature density contrasts. It also occurs through torques acting on the inner core reorienting its non-spherical shape. These in situ mass changes lead to global gravity variations, and global deformations (inducing additional gravity variations) occur in order to maintain the mechanical equilibrium of the whole Earth. Changes in Earth's rotation vector (and thus of the global centrifugal potential) induced by core flows are an additional source of global deformations and associated gravity changes originating from core dynamics. Here, we review how each of these different core processes operates, how gravity changes and ground deformations from each could be reconstructed, as well as ways to estimate their amplitudes. Based on our current understanding of core dynamics, we show that, at spherical harmonic degree 2, core processes contribute to gravity variations and ground deformations that are approximately a factor 10 smaller than those observed and caused by dynamical processes within the fluid layers at the Earth's surface. The larger the harmonic degree, the smaller is the contribution from the core. Extracting a signal of core origin requires the accurate removal of all contributions from surface processes, which remains a challenge.

ARTICLE HIGHLIGHTS

Dynamical processes in Earth's fluid core lead to global gravity variations and surface ground deformationsWe review how these processes operate, how signals of core origin can be reconstructed and estimate their amplitudesCore signals are a factor 10 smaller than the observed signals; extracting a signal of core origin remains a challenge.

Rapid geomagnetic changes inferred from Earth observations and numerical simulations

Extreme variations in the direction of Earth’s magnetic field contain important information regarding the operation of the geodynamo. Paleomagnetic studies have reported rapid directional changes reaching 1° yr⁻¹, although the observations are controversial and their relation to physical processes in Earth’s core unknown. Here we show excellent agreement between amplitudes and latitude ranges of extreme directional changes in a suite of geodynamo simulations and a recent observational field model spanning the past 100 kyrs. Remarkably, maximum rates of directional change reach  ~10° yr⁻¹, typically during times of decreasing field strength, almost 100 times faster than current changes. Detailed analysis of the simulations and a simple analogue model indicate that extreme directional changes are associated with movement of reversed flux across the core surface. Our results demonstrate that such rapid variations are compatible with the physics of the dynamo process and suggest that future searches for rapid directional changes should focus on low latitudes.

The authors here use numerical simulations combined with a time-dependent model of Earth’s magnetic field spanning the last 100 kyrs. They identify field directional changes to be 10 times faster than previously thought.

Anomalous SmKS induced by postcritical reflection and refraction at the core-mantle boundary

Earth's outer core is generally thought to be a well-mixed liquid consisting mostly of iron and a small amount of lighter elements. Recent seismic studies using SmKS waves show that the top a few hundred kilometers of the outer core possess a P-wave velocity slightly lower than the PREM model, which cannot be explained by self-compression of a chemically homogeneous outer core. We investigated the SmKS waveforms of a deep earthquake occurring beneath South America recorded by a large and dense seismic array in China, and measured the differential arrival times of the SmKS pairs. We found significant waveform distortion of the SmKS caused by postcritical refraction and reflection at the core-mantle boundary. This waveform distortion can introduce significant bias to the measured differential times, leading to incorrect estimate of P-wave velocity of the outer core. Whether stable stratification is occurring in outer core or not requires further seismic investigations.

Evidence for Fe-Si-O liquid immiscibility at deep Earth pressures

Seismic observations suggest that the uppermost region of Earth's liquid outer core is buoyant, with slower velocities than the bulk outer core. One possible mechanism for the formation of a stably stratified layer is immiscibility in molten iron alloy systems, which has yet to be demonstrated at core pressures. We find immiscibility between liquid Fe-Si and Fe-Si-O persisting to at least 140 GPa through a combination of laser-heated diamond-anvil cell experiments and first-principles molecular dynamics simulations. High-pressure immiscibility in the Fe-Si-O system may explain a stratified layer atop the outer core, complicate differentiation and evolution of the deep Earth, and affect the structure and intensity of Earth's magnetic field. Our results support silicon and oxygen as coexisting light elements in the core and suggest that [Formula: see text] does not crystallize out of molten Fe-Si-O at the core-mantle boundary.

Seismically determined elastic parameters for Earth's outer core

Turbulent convection of the liquid iron alloy outer core generates Earth's magnetic field and supplies heat to the mantle. The exact composition of the iron alloy is fundamentally linked to the processes powering the convection and can be constrained by its seismic properties. Discrepancies between seismic models determined using body waves and normal modes show that these properties are not yet fully agreed upon. In addition, technical challenges in experimentally measuring the equation-of-state (EoS) parameters of liquid iron alloys at high pressures and temperatures further complicate compositional inferences. We directly infer EoS parameters describing Earth's outer core from normal mode center frequency observations and present the resulting Elastic Parameters of the Outer Core (EPOC) seismic model. Unlike alternative seismic models, ours requires only three parameters and guarantees physically realistic behavior with increasing pressure for a well-mixed homogeneous material along an isentrope, consistent with the outer core's condition. We show that EPOC predicts available normal mode frequencies better than the Preliminary Reference Earth Model (PREM) while also being more consistent with body wave-derived models, eliminating a long-standing discrepancy. The velocity at the top of the outer core is lower, and increases with depth more steeply, in EPOC than in PREM, while the density in EPOC is higher than that in PREM across the outer core. The steeper profiles and higher density imply that the outer core comprises a lighter but more compressible alloy than that inferred for PREM. Furthermore, EPOC's steeper velocity gradient explains differential SmKS body wave travel times better than previous one-dimensional global models, without requiring an anomalously slow ~90- to 450-km-thick layer at the top of the outer core.

Geomagnetic spikes on the core-mantle boundary

Extreme variations of Earth's magnetic field occurred in the Levant region around 1000 BC, when the field intensity rapidly rose and fell by a factor of 2. No coherent link currently exists between this intensity spike and the global field produced by the core geodynamo. Here we show that the Levantine spike must span >60° longitude at Earth's surface if it originates from the core-mantle boundary (CMB). Several low intensity data are incompatible with this geometric bound, though age uncertainties suggest these data could have sampled the field before the spike emerged. Models that best satisfy energetic and geometric constraints produce CMB spikes 8-22° wide, peaking at O(100) mT. We suggest that the Levantine spike reflects an intense CMB flux patch that grew in place before migrating northwest, contributing to growth of the dipole field. Estimates of Ohmic heating suggest that diffusive processes likely govern the ultimate decay of geomagnetic spikes.

Highly siderophile elements were stripped from Earth's mantle by iron sulfide segregation

Highly siderophile elements (HSEs) are strongly depleted in the bulk silicate Earth (BSE) but are present in near-chondritic relative abundances. The conventional explanation is that the HSEs were stripped from the mantle by the segregation of metal during core formation but were added back in near-chondritic proportions by late accretion, after core formation had ceased. Here we show that metal-silicate equilibration and segregation during Earth's core formation actually increased HSE mantle concentrations because HSE partition coefficients are relatively low at the high pressures of core formation within Earth. The pervasive exsolution and segregation of iron sulfide liquid from silicate liquid (the "Hadean matte") stripped magma oceans of HSEs during cooling and crystallization, before late accretion, and resulted in slightly suprachondritic palladium/iridium and ruthenium/iridium ratios.

Experimental constraints on light elements in the Earth's outer core

Earth's outer core is liquid and dominantly composed of iron and nickel (~5-10 wt%). Its density, however, is ~8% lower than that of liquid iron, and requires the presence of a significant amount of light element(s). A good way to specify the light element(s) is a direct comparison of density and sound velocity measurements between seismological data and those of possible candidate compositions at the core conditions. We report the sound velocity measurements of a model core composition in the Fe-Ni-Si system at the outer core conditions by shock-wave experiments. Combining with the previous studies, we found that the best estimate for the outer core's light elements is ~6 wt% Si, ~2 wt% S, and possible ~1-2.5 wt% O. This composition satisfies the requirements imposed by seismology, geochemistry, and some models of the early core formation. This finding may help us to further constrain the thermal structure of the Earth and the models of Earth's core formation.

Gyre-driven decay of the Earth's magnetic dipole

Direct observations indicate that the magnitude of the Earth's magnetic axial dipole has decreased over the past 175 years; it is now 9% weaker than it was in 1840. Here we show how the rate of dipole decay may be controlled by a planetary-scale gyre in the liquid metal outer core. The gyre's meridional limbs on average transport normal polarity magnetic flux equatorward and reverse polarity flux poleward. Asymmetry in the geomagnetic field, due to the South Atlantic Anomaly, is essential to the proposed mechanism. We find that meridional flux advection accounts for the majority of the dipole decay since 1840, especially during times of rapid decline, with magnetic diffusion making an almost steady contribution generally of smaller magnitude. Based on the morphology of the present field, and the persistent nature of the gyre, the current episode of dipole decay looks set to continue, at least for the next few decades.

The magnitude of the Earth's magnetic dipole has decreased by 9% over the past 175 years. Here, the authors suggest that the rate of dipole decay is controlled by a huge gyre in the liquid metal outer core acting on a field asymmetry, and that decay is set to continue for the next few decades.

Seismological evidence for a non-monotonic velocity gradient in the topmost outer core

The Earth's core is mostly an Fe-Ni alloy with a fraction of light elements (~10 wt%, mainly O, S and Si). Accumulation of these light elements under the core-mantle boundary (CMB) may lead to chemical stratification. Seismic observations have been presented both for and against the stratification in the topmost region of the outer core. Here we investigate the structure under the CMB using differential travel times between SKKS and S3KS waves. We obtain 606 high-quality S3KS-SKKS differential travel times with global path coverage. Result from a Bayesian inversion of these differential times indicates that the seismic velocity in the top 800 km of the outer core is ~0.07% on average lower than that in model PREM. The depth-dependent velocity profile, in particular a low-velocity zone of up to ~0.25% lower than PREM at ~80 km below the CMB, strongly favors the existence of stratification at the top of the outer core.

Electrical resistivity and thermal conductivity of liquid Fe alloys at high P and T, and heat flux in Earth's core

Earth's magnetic field is sustained by magnetohydrodynamic convection within the metallic liquid core. In a thermally advecting core, the fraction of heat available to drive the geodynamo is reduced by heat conducted along the core geotherm, which depends sensitively on the thermal conductivity of liquid iron and its alloys with candidate light elements. The thermal conductivity for Earth's core is very poorly constrained, with current estimates based on a set of scaling relations that were not previously tested at high pressures. We perform first-principles electronic structure computations to determine the thermal conductivity and electrical resistivity for Fe, Fe-Si, and Fe-O liquid alloys. Computed resistivity agrees very well with existing shock compression measurements and shows strong dependence on light element concentration and type. Thermal conductivity at pressure and temperature conditions characteristic of Earth's core is higher than previous extrapolations. Conductive heat flux near the core-mantle boundary is comparable to estimates of the total heat flux from the core but decreases with depth, so that thermally driven flow would be constrained to greater depths in the absence of an inner core.