Spin crossover and iron-rich silicate melt in the Earth's deep mantle

Trace element partitioning in a deep magma ocean and the origin of the Hf-Nd mantle array

Crystallization in Earth's deep magma ocean could have caused trace element fractionation in the lower mantle that might be inherited to the isotopic compositions of the present-day mantle. However, the trace element partitioning has been experimentally investigated only up to the uppermost lower-mantle pressures. Here, we determined the bridgmanite/melt partition coefficients D of La, Nd, Sm, Lu, and Hf from 24 to 115 gigapascals, covering the wide pressure range of the lower mantle. Results demonstrate substantial reductions in DLu and DHf from >1 to ≪1 with increasing pressure to 91 gigapascals. We also found DLu/DHf > 1 and DSm/DNd < 1 under deep lower-mantle conditions, evolving melts toward low Lu/Hf and high Sm/Nd ratios by crystallizing bridgmanite. If residual melts form a dense hidden reservoir in the lowermost mantle, the complementary accessible mantle has the Hf and Nd isotopic compositions matching the observed terrestrial mantle array that deviates from the bulk silicate Earth reference.

Earth's core-mantle boundary shaped by the crystallization of a hydrous terrestrial magma ocean

Enigmatic anomalous structures of Earth's lowermost mantle may have been incubated at the hydrous magma ocean of the Hadean eon.

Moon-forming impactor as a source of Earth's basal mantle anomalies

Seismic images of Earth's interior have revealed two continent-sized anomalies with low seismic velocities, known as the large low-velocity provinces (LLVPs), in the lowermost mantle1. The LLVPs are often interpreted as intrinsically dense heterogeneities that are compositionally distinct from the surrounding mantle2. Here we show that LLVPs may represent buried relics of Theia mantle material (TMM) that was preserved in proto-Earth's mantle after the Moon-forming giant impact3. Our canonical giant-impact simulations show that a fraction of Theia's mantle could have been delivered to proto-Earth's solid lower mantle. We find that TMM is intrinsically 2.0-3.5% denser than proto-Earth's mantle based on models of Theia's mantle and the observed higher FeO content of the Moon. Our mantle convection models show that dense TMM blobs with a size of tens of kilometres after the impact can later sink and accumulate into LLVP-like thermochemical piles atop Earth's core and survive to the present day. The LLVPs may, thus, be a natural consequence of the Moon-forming giant impact. Because giant impacts are common at the end stages of planet accretion, similar mantle heterogeneities caused by impacts may also exist in the interiors of other planetary bodies.

Ultrafast x-ray detection of low-spin iron in molten silicate under deep planetary interior conditions

The spin state of Fe can alter the key physical properties of silicate melts, affecting the early differentiation and the dynamic stability of the melts in the deep rocky planets. The low-spin state of Fe can increase the affinity of Fe for the melt over the solid phases and the electrical conductivity of melt at high pressures. However, the spin state of Fe has never been measured in dense silicate melts due to experimental challenges. We report detection of dominantly low-spin Fe in dynamically compressed olivine melt at 150 to 256 gigapascals and 3000 to 6000 kelvin using laser-driven shock wave compression combined with femtosecond x-ray diffraction and x-ray emission spectroscopy using an x-ray free electron laser. The observation of dominantly low-spin Fe supports gravitationally stable melt in the deep mantle and generation of a dynamo from the silicate melt portion of rocky planets.

Globally distributed subducted materials along the Earth's core-mantle boundary: Implications for ultralow velocity zones

Ultralow velocity zones (ULVZs) are the most anomalous structures within the Earth's interior; however, given the wide range of associated characteristics (thickness and composition) reported by previous studies, the origins of ULVZs have been debated for decades. Using a recently developed seismic analysis approach, we find widespread, variable ULVZs along the core-mantle boundary (CMB) beneath a largely unsampled portion of the Southern Hemisphere. Our study region is not beneath current or recent subduction zones, but our mantle convection simulations demonstrate how heterogeneous accumulations of previously subducted materials could form on the CMB and explain our seismic observations. We further show that subducted materials can be globally distributed throughout the lowermost mantle with variable concentrations. These subducted materials, advected along the CMB, can provide an explanation for the distribution and range of reported ULVZ properties.

Calcium dissolution in bridgmanite in the Earth's deep mantle

Accurate knowledge of the mineralogy is essential for understanding the lower mantle, which represents more than half of Earth's volume. CaSiO₃ perovskite is believed to be the third-most-abundant mineral throughout the lower mantle, following bridgmanite and ferropericlase^1-3. Here we experimentally show that the calcium solubility in bridgmanite increases steeply at about 2,300 kelvin and above 40 gigapascals to a level sufficient for a complete dissolution of all CaSiO₃ component in pyrolite into bridgmanite, resulting in the disappearance of CaSiO₃ perovskite at depths greater than about 1,800 kilometres along the geotherm^4,5. Hence we propose a change from a two-perovskite domain (TPD; bridgmanite plus CaSiO₃ perovskite) at the shallower lower mantle to a single-perovskite domain (SPD; calcium-rich bridgmanite) at the deeper lower mantle. Iron seems to have a key role in increasing the calcium solubility in bridgmanite. The temperature-driven nature can cause large lateral variations in the depth of the TPD-to-SPD change in response to temperature variations (by more than 500 kilometres). Furthermore, the SPD should have been thicker in the past when the mantle was warmer. Our finding requires revision of the deep-mantle mineralogy models and will have an impact on our understanding of the composition, structure, dynamics and evolution of the region.

Primitive noble gases sampled from ocean island basalts cannot be from the Earth's core

Noble gas isotopes in plumes require a source of primitive volatiles largely isolated in the Earth for 4.5 Gyrs. Among the proposed reservoirs, the core is gaining interest in the absence of robust geochemical and geophysical evidence for a mantle source. This is supported by partitioning data showing that sufficient He and Ne could have been incorporated into the core to source plumes today. Here we perform ab initio calculations on the partitioning of He, Ne, Ar, Kr and Xe between liquid iron and silicate melt under core forming conditions. For He our results are consistent with previous studies allowing for substantial amounts of He in the core. In contrast, the partition coefficient for Ne is three orders of magnitude lower than He. This very low partition coefficient would result in a ³He/²²Ne ratio of ~10³ in the core, far higher than observed in ocean island basalts (OIBs). We conclude that the core is not the source of noble gases in OIBs.

Kilometer-scale structure on the core-mantle boundary near Hawaii

The lowermost mantle right above the core-mantle boundary is highly heterogeneous containing multiple poorly understood seismic features. The smallest but most extreme heterogeneities yet observed are 'Ultra-Low Velocity Zones' (ULVZ). We exploit seismic shear waves that diffract along the core-mantle boundary to provide new insight into these enigmatic structures. We measure a rare core-diffracted signal refracted by a ULVZ at the base of the Hawaiian mantle plume at unprecedentedly high frequencies. This signal shows remarkably longer time delays at higher compared to lower frequencies, indicating a pronounced internal variability inside the ULVZ. Utilizing the latest computational advances in 3D waveform modeling, here we show that we are able to model this high-frequency signal and constrain high-resolution ULVZ structure on the scale of kilometers, for the first time. This new observation suggests a chemically distinct ULVZ with increasing iron content towards the core-mantle boundary, which has implications for Earth's early evolutionary history and core-mantle interaction.

Slab control on the mega-sized North Pacific ultra-low velocity zone

Ultra-low velocity zones (ULVZs) are localized small-scale patches with extreme physical properties at the core-mantle boundary that often gather at the margins of Large Low Velocity Provinces (LLVPs). Recent studies have discovered several mega-sized ULVZs with a lateral dimension of ~900 km. However, the detailed structures and physical properties of these ULVZs and their relationship to LLVP edges are not well constrained and their formation mechanisms are poorly understood. Here, we break the degeneracy between the size and velocity perturbation of a ULVZ using two orthogonal seismic ray paths, and thereby discover a mega-sized ULVZ at the northern edge of the Pacific LLVP. The ULVZ is almost double the size of a previously imaged ULVZ in this region, but with half of the shear velocity reduction. This mega-sized ULVZ has accumulated due to stable mantle flow converging at the LLVP edge driven by slab-debris in the lower mantle. Such flow also develops the subvertical north-tilting edge of the Pacific LLVP.

Magma Ocean Evolution of the TRAPPIST-1 Planets

Recent observations of the potentially habitable planets TRAPPIST-1 e, f, and g suggest that they possess large water mass fractions of possibly several tens of weight percent of water, even though the host star's activity should drive rapid atmospheric escape. These processes can photolyze water, generating free oxygen and possibly desiccating the planet. After the planets formed, their mantles were likely completely molten with volatiles dissolving and exsolving from the melt. To understand these planets and prepare for future observations, the magma ocean phase of these worlds must be understood. To simulate these planets, we have combined existing models of stellar evolution, atmospheric escape, tidal heating, radiogenic heating, magma-ocean cooling, planetary radiation, and water-oxygen-iron geochemistry. We present MagmOc, a versatile magma-ocean evolution model, validated against the rocky super-Earth GJ 1132b and early Earth. We simulate the coupled magma-ocean atmospheric evolution of TRAPPIST-1 e, f, and g for a range of tidal and radiogenic heating rates, as well as initial water contents between 1 and 100 Earth oceans. We also reanalyze the structures of these planets and find they have water mass fractions of 0-0.23, 0.01-0.21, and 0.11-0.24 for planets e, f, and g, respectively. Our model does not make a strong prediction about the water and oxygen content of the atmosphere of TRAPPIST-1 e at the time of mantle solidification. In contrast, the model predicts that TRAPPIST-1 f and g would have a thick steam atmosphere with a small amount of oxygen at that stage. For all planets that we investigated, we find that only 3-5% of the initial water will be locked in the mantle after the magma ocean solidified.

Investigating Magma Ocean Solidification on Earth Through Laser-Heated Diamond Anvil Cell Experiments

We carried out a series of silicate fractional crystallization experiments at lower mantle pressures using the laser-heated diamond anvil cell. Phase relations and the compositional evolution of the cotectic melt and equilibrium solids along the liquid line of descent were determined and used to assemble the melting phase diagram. In a pyrolitic magma ocean, the first mineral to crystallize in the deep mantle is iron-depleted calcium-bearing bridgmanite. From the phase diagram, we estimate that the initial 33%-36% of the magma ocean will crystallize to form such a buoyant bridgmanite. Substantial calcium solubility in bridgmanite is observed up to 129 GPa, and significantly delays the crystallization of the calcium silicate perovskite phase during magma ocean solidification. Residual melts are strongly iron-enriched as crystallization proceeds, making them denser than any of the coexisting solids at deep mantle conditions, thus supporting the terrestrial basal magma ocean hypothesis (Labrosse et al., 2007).

Reconciling metal-silicate partitioning and late accretion in the Earth

Highly siderophile elements (HSE), including platinum, provide powerful geochemical tools for studying planet formation. Late accretion of chondritic components to Earth after core formation has been invoked as the main source of mantle HSE. However, core formation could also have contributed to the mantle's HSE content. Here we present measurements of platinum metal-silicate partitioning coefficients, obtained from laser-heated diamond anvil cell experiments, which demonstrate that platinum partitioning into metal is lower at high pressures and temperatures. Consequently, the mantle was likely enriched in platinum immediately following core-mantle differentiation. Core formation models that incorporate these results and simultaneously account for collateral geochemical constraints, lead to excess platinum in the mantle. A subsequent process such as iron exsolution or sulfide segregation is therefore required to remove excess platinum and to explain the mantle's modern HSE signature. A vestige of this platinum-enriched mantle can potentially account for ¹⁸⁶Os-enriched ocean island basalt lavas.

Formation of large low shear velocity provinces through the decomposition of oxidized mantle

Large Low Shear Velocity Provinces (LLSVPs) in the lowermost mantle are key to understanding the chemical composition and thermal structure of the deep Earth, but their origins have long been debated. Bridgmanite, the most abundant lower-mantle mineral, can incorporate extensive amounts of iron (Fe) with effects on various geophysical properties. Here our high-pressure experiments and ab initio calculations reveal that a ferric-iron-rich bridgmanite coexists with an Fe-poor bridgmanite in the 90 mol% MgSiO₃-10 mol% Fe₂O₃ system, rather than forming a homogeneous single phase. The Fe³⁺-rich bridgmanite has substantially lower velocities and a higher V_P/V_S ratio than MgSiO₃ bridgmanite under lowermost-mantle conditions. Our modeling shows that the enrichment of Fe³⁺-rich bridgmanite in a pyrolitic composition can explain the observed features of the LLSVPs. The presence of Fe³⁺-rich materials within LLSVPs may have profound effects on the deep reservoirs of redox-sensitive elements and their isotopes.

Buoyancy and Structure of Volatile-Rich Silicate Melts

The early Earth was marked by at least one global magma ocean. Melt buoyancy played a major role for its evolution. Here we model the composition of the magma ocean using a six-component pyrolite melt, to which we add volatiles in the form of carbon as molecular CO or CO2 and hydrogen as molecular H2O or through substitution for magnesium. We compute the density relations from first-principles molecular dynamics simulations. We find that the addition of volatiles renders all the melts more buoyant compared to the reference volatile-free pyrolite. The effect is pressure dependent, largest at the surface, decreasing to about 20 GPa, and remaining roughly constant to 135 GPa. The increased buoyancy would have enhanced convection and turbulence, and thus promoted the chemical exchanges of the magma ocean with the early atmosphere. We determine the partial molar volume of both H2O and CO2 throughout the magma ocean conditions. We find a pronounced dependence with temperature at low pressures, whereas at megabar pressures the partial molar volumes are independent of temperature. At all pressures, the polymerization of the silicate melt is strongly affected by the amount of oxygen added to the system while being only weakly affected by the specific type of volatile added.

200-mm segmented cylindrical figured crystal for von Hamos x-ray spectrometer

A von Hamos Bragg crystal spectrometer at 1C beamline of Pohang Accelerator Laboratory for x-ray emission spectroscopy (XES) is described. Diced Si crystals of different orientations ([111], [110], [100], and [311]) are glued onto a planoconcave glass substrate having 250/500 mm radius of curvature. To enhance the spectrometer efficiency, the length of the crystal analyzer is kept 200 mm. The emission spectra of Cu foil and Fe foil and elastic scattering from Al foil are measured using the von Hamos geometry in which curved crystals disperse the x-rays. Spectrometer efficiency and energy resolution are measured at various x-ray photon energies. X-rays are incident at 6.54 keV, 9.00 keV, 9.205 keV, and 11.51 keV for Si(440), Si(444), Si(800), and Si(933) crystal analyzers, respectively. The cylindrical figured analyzer is placed near 80° with respect to the sample, which gives better energy resolution. The spectrometer efficiency of the Si(444) crystal analyzer increases by ∼2 times when the length of the analyzer is increased from 100 mm to 200 mm. Furthermore, to measure Fe Kα₁, Kα₂, and K_β simultaneously, we made a mixed crystal analyzer in which alternative strips of Si[111] and Si[110] are glued onto one preshaped cylindrical substrate. The enhanced efficiency and simultaneous measurement of K_α and K_β emission lines will give an edge over in situ and time-resolved x-ray emission spectroscopy studies. The information extracted with a high efficiency spectrometer from low intensity XES emission lines will be useful for the in situ elemental characterization in catalytic reactions.

An improved laboratory-based x-ray absorption fine structure and x-ray emission spectrometer for analytical applications in materials chemistry research

X-ray absorption fine structure (XAFS) and x-ray emission spectroscopy (XES) are advanced x-ray spectroscopies that impact a wide range of disciplines. However, unlike the majority of other spectroscopic methods, XAFS and XES are accompanied by an unusual access model, wherein the dominant use of the technique is for premier research studies at world-class facilities, i.e., synchrotron x-ray light sources. In this paper, we report the design and performance of an improved XAFS and XES spectrometer based on the general conceptual design of Seidler et al. [Rev. Sci. Instrum. 85, 113906 (2014)]. New developments include reduced mechanical degrees of freedom, much-increased flux, and a wider Bragg angle range to enable extended x-ray absorption fine structure (EXAFS) measurement and analysis for the first time with this type of modern laboratory XAFS configuration. This instrument enables a new class of routine applications that are incompatible with the mission and access model of the synchrotron light sources. To illustrate this, we provide numerous examples of x-ray absorption near edge structure (XANES), EXAFS, and XES results for a variety of problems and energy ranges. Highlights include XAFS and XES measurements of battery electrode materials, EXAFS of Ni with full modeling of results to validate monochromator performance, valence-to-core XES for 3d transition metal compounds, and uranium XANES and XES for different oxidation states. Taken en masse, these results further support the growing perspective that modern laboratory-based XAFS and XES have the potential to develop a new branch of analytical chemistry.

Altered chemistry of oxygen and iron under deep Earth conditions

A drastically altered chemistry was recently discovered in the Fe-O-H system under deep Earth conditions, involving the formation of iron superoxide (FeO₂Hx with x = 0 to 1), but the puzzling crystal chemistry of this system at high pressures is largely unknown. Here we present evidence that despite the high O/Fe ratio in FeO₂Hx, iron remains in the ferrous, spin-paired and non-magnetic state at 60-133 GPa, while the presence of hydrogen has minimal effects on the valence of iron. The reduced iron is accompanied by oxidized oxygen due to oxygen-oxygen interactions. The valence of oxygen is not -2 as in all other major mantle minerals, instead it varies around -1. This result indicates that like iron, oxygen may have multiple valence states in our planet's interior. Our study suggests a possible change in the chemical paradigm of how oxygen, iron, and hydrogen behave under deep Earth conditions.

Metastable silica high pressure polymorphs as structural proxies of deep Earth silicate melts

Modelling of processes involving deep Earth liquids requires information on their structures and compression mechanisms. However, knowledge of the local structures of silicates and silica (SiO2) melts at deep mantle conditions and of their densification mechanisms is still limited. Here we report the synthesis and characterization of metastable high-pressure silica phases, coesite-IV and coesite-V, using in situ single-crystal X-ray diffraction and ab initio simulations. Their crystal structures are drastically different from any previously considered models, but explain well features of pair-distribution functions of highly densified silica glass and molten basalt at high pressure. Built of four, five-, and six-coordinated silicon, coesite-IV and coesite-V contain SiO6 octahedra, which, at odds with 3rd Pauling's rule, are connected through common faces. Our results suggest that possible silicate liquids in Earth's lower mantle may have complex structures making them more compressible than previously supposed.

Stability and nature of the volume collapse of ε-Fe2O3 under extreme conditions

Iron oxides are among the major constituents of the deep Earth's interior. Among them, the epsilon phase of Fe₂O₃ is one of the less studied polymorphs and there is a lack of information about its structural, electronic and magnetic transformations at extreme conditions. Here we report the precise determination of its equation of state and a deep analysis of the evolution of the polyhedral units under compression, thanks to the agreement between our experiments and ab-initio simulations. Our results indicate that this material, with remarkable magnetic properties, is stable at pressures up to 27 GPa. Above 27 GPa, a volume collapse has been observed and ascribed to a change of the local environment of the tetrahedrally coordinated iron towards an octahedral coordination, finding evidence for a different iron oxide polymorph.

Pressure Sensor with a Color Change at Room Temperature Based on Spin-Crossover Behavior

Two new iron(II) complexes with 1D chain and 2D network structures have been successfully synthesized and characterized. One of the complexes exhibits a pressure-induced spin-crossover property with a reversible color change from white to purple at room temperature. The special property makes it a suitable candidate as a pressure sensor.

Hadean silicate differentiation preserved by anomalous 142Nd/144Nd ratios in the Réunion hotspot source

Active volcanic hotspots can tap into domains in Earth's deep interior that were formed more than two billion years ago. High-precision data on variability in tungsten isotopes have shown that some of these domains resulted from differentiation events that occurred within the first fifty million years of Earth history. However, it has not proved easy to resolve analogous variability in neodymium isotope compositions that would track regions of Earth's interior whose composition was established by events occurring within roughly the first five hundred million years of Earth history. Here we report ¹⁴²Nd/¹⁴⁴Nd ratios for Réunion Island igneous rocks, some of which are resolvably either higher or lower than the ratios in modern upper-mantle domains. We also find that Réunion ¹⁴²Nd/¹⁴⁴Nd ratios correlate with helium-isotope ratios (³He/⁴He), suggesting parallel behaviour of these isotopic systems during very early silicate differentiation, perhaps as early as 4.39 billion years ago. The range of ¹⁴²Nd/¹⁴⁴Nd ratios in Réunion basalts is inconsistent with a single-stage differentiation process, and instead requires mixing of a conjugate melt and residue formed in at least one melting event during the Hadean eon, 4.56 billion to 4 billion years ago. Efficient post-Hadean mixing nearly erased the ancient, anomalous ¹⁴²Nd/¹⁴⁴Nd signatures, and produced the relatively homogeneous ¹⁴³Nd/¹⁴⁴Nd composition that is characteristic of Réunion basalts. Our results show that Réunion magmas tap into a particularly ancient, primitive source compared with other volcanic hotspots, offering insight into the formation and preservation of ancient heterogeneities in Earth's interior.

Pressure-induced structural change in MgSiO3 glass at pressures near the Earth's core-mantle boundary

Knowledge of the structure and properties of silicate magma under extreme pressure plays an important role in understanding the nature and evolution of Earth's deep interior. Here we report the structure of MgSiO₃ glass, considered an analog of silicate melts, up to 111 GPa. The first (r1) and second (r2) neighbor distances in the pair distribution function change rapidly, with r1 increasing and r2 decreasing with pressure. At 53-62 GPa, the observed r1 and r2 distances are similar to the Si-O and Si-Si distances, respectively, of crystalline MgSiO₃ akimotoite with edge-sharing SiO₆ structural motifs. Above 62 GPa, r1 decreases, and r2 remains constant, with increasing pressure until 88 GPa. Above this pressure, r1 remains more or less constant, and r2 begins decreasing again. These observations suggest an ultrahigh-pressure structural change around 88 GPa. The structure above 88 GPa is interpreted as having the closest edge-shared SiO₆ structural motifs similar to those of the crystalline postperovskite, with densely packed oxygen atoms. The pressure of the structural change is broadly consistent with or slightly lower than that of the bridgmanite-to-postperovskite transition in crystalline MgSiO₃ These results suggest that a structural change may occur in MgSiO₃ melt under pressure conditions corresponding to the deep lower mantle.

SiO_{2} Glass Density to Lower-Mantle Pressures

The convection or settling of matter in the deep Earth's interior is mostly constrained by density variations between the different reservoirs. Knowledge of the density contrast between solid and molten silicates is thus of prime importance to understand and model the dynamic behavior of the past and present Earth. SiO_{2} is the main constituent of Earth's mantle and is the reference model system for the behavior of silicate melts at high pressure. Here, we apply our recently developed x-ray absorption technique to the density of SiO_{2} glass up to 110 GPa, doubling the pressure range for such measurements. Our density data validate recent molecular dynamics simulations and are in good agreement with previous experimental studies conducted at lower pressure. Silica glass rapidly densifies up to 40 GPa, but the density trend then flattens to become asymptotic to the density of SiO_{2} minerals above 60 GPa. The density data present two discontinuities at ∼17 and ∼60  GPa that can be related to a silicon coordination increase from 4 to a mixed 5/6 coordination and from 5/6 to sixfold, respectively. SiO_{2} glass becomes denser than MgSiO_{3} glass at ∼40  GPa, and its density becomes identical to that of MgSiO_{3} glass above 80 GPa. Our results on SiO_{2} glass may suggest that a variation of SiO_{2} content in a basaltic or pyrolitic melt with pressure has at most a minor effect on the final melt density, and iron partitioning between the melts and residual solids is the predominant factor that controls melt buoyancy in the lowermost mantle.

Coulomb Correlations Intertwined with Spin and Orbital Excitations in LaCoO_{3}

We carried out temperature-dependent (20-550 K) measurements of resonant inelastic x-ray scattering on LaCoO_{3} to investigate the evolution of its electronic structure across the spin-state crossover. In combination with charge-transfer multiplet calculations, we accurately quantified the renomalized crystal-field excitation energies and spin-state populations. We show that the screening of the effective on-site Coulomb interaction of 3d electrons is orbital selective and coupled to the spin-state crossover in LaCoO_{3}. The results establish that the gradual spin-state crossover is associated with a relative change of Coulomb energy versus bandwidth, leading to a Mott-type insulator-to-metal transition.

Compositionally-distinct ultra-low velocity zones on Earth's core-mantle boundary

The Earth’s lowermost mantle large low velocity provinces are accompanied by small-scale ultralow velocity zones in localized regions on the core-mantle boundary. Large low velocity provinces are hypothesized to be caused by large-scale compositional heterogeneity (i.e., thermochemical piles). The origin of ultralow velocity zones, however, remains elusive. Here we perform three-dimensional geodynamical calculations to show that the current locations and shapes of ultralow velocity zones are related to their cause. We find that the hottest lowermost mantle regions are commonly located well within the interiors of thermochemical piles. In contrast, accumulations of ultradense compositionally distinct material occur as discontinuous patches along the margins of thermochemical piles and have asymmetrical cross-sectional shape. Furthermore, the lateral morphology of these patches provides insight into mantle flow directions and long-term stability. The global distribution and large variations of morphology of ultralow velocity zones validate a compositionally distinct origin for most ultralow velocity zones.

Ultralow velocity zones are detected on the core-mantle boundary, but their origin is enigmatic. Here, the authors find that the global distribution and large variations of morphology of ultralow velocity zones are consistent with most having a compositionally-distinct origin.

Solid-liquid density and spin crossovers in (Mg, Fe)O system at deep mantle conditions

The low/ultralow-velocity zones in the Earth's mantle can be explained by the presence of partial melting, critically depending on density contrast between the melt and surrounding solid mantle. Here, first-principles molecular dynamics simulations of (Mg, Fe) O ferropericlase in the solid and liquid states show that their densities increasingly approach each other as pressure increases. The isochemical density difference between them diminishes from 0.78 (±0.7) g/cm3 at zero pressure (3000 K) to 0.16 (±0.04) g/cm3 at 135 GPa (4000 K) for pure and alloyed compositions containing up to 25% iron. The simulations also predict a high-spin to low-spin transition of iron in the liquid ferropericlase gradually occurring over a pressure interval centered at 55 GPa (4000 K) accompanied by a density increase of 0.14 (±0.02) g/cm3. Temperature tends to widen the transition to higher pressure. The estimated iron partition coefficient between the solid and liquid ferropericlase varies from 0.3 to 0.6 over the pressure range of 23 to 135 GPa. Based on these results, an excess of as low as 5% iron dissolved in the liquid could cause the solid-liquid density crossover at conditions of the lowermost mantle.

Origins of ultralow velocity zones through slab-derived metallic melt

Understanding the ultralow velocity zones (ULVZs) places constraints on the chemical composition and thermal structure of deep Earth and provides critical information on the dynamics of large-scale mantle convection, but their origin has remained enigmatic for decades. Recent studies suggest that metallic iron and carbon are produced in subducted slabs when they sink beyond a depth of 250 km. Here we show that the eutectic melting curve of the iron-carbon system crosses the current geotherm near Earth's core-mantle boundary, suggesting that dense metallic melt may form in the lowermost mantle. If concentrated into isolated patches, such melt could produce the seismically observed density and velocity features of ULVZs. Depending on the wetting behavior of the metallic melt, the resultant ULVZs may be short-lived domains that are replenished or regenerated through subduction, or long-lasting regions containing both metallic and silicate melts. Slab-derived metallic melt may produce another type of ULVZ that escapes core sequestration by reacting with the mantle to form iron-rich postbridgmanite or ferropericlase. The hypotheses connect peculiar features near Earth's core-mantle boundary to subduction of the oceanic lithosphere through the deep carbon cycle.

Fate of MgSiO3 melts at core-mantle boundary conditions

One key for understanding the stratification in the deep mantle lies in the determination of the density and structure of matter at high pressures, as well as the density contrast between solid and liquid silicate phases. Indeed, the density contrast is the main control on the entrainment or settlement of matter and is of fundamental importance for understanding the past and present dynamic behavior of the deepest part of the Earth's mantle. Here, we adapted the X-ray absorption method to the small dimensions of the diamond anvil cell, enabling density measurements of amorphous materials to unprecedented conditions of pressure. Our density data for MgSiO3 glass up to 127 GPa are considerably higher than those previously derived from Brillouin spectroscopy but validate recent ab initio molecular dynamics simulations. A fourth-order Birch-Murnaghan equation of state reproduces our experimental data over the entire pressure regime of the mantle. At the core-mantle boundary (CMB) pressure, the density of MgSiO3 glass is 5.48 ± 0.18 g/cm(3), which is only 1.6% lower than that of MgSiO3 bridgmanite at 5.57 g/cm(3), i.e., they are the same within the uncertainty. Taking into account the partitioning of iron into the melt, we conclude that melts are denser than the surrounding solid phases in the lowermost mantle and that melts will be trapped above the CMB.

Low core-mantle boundary temperature inferred from the solidus of pyrolite

The melting temperature of Earth's mantle provides key constraints on the thermal structures of both the mantle and the core. Through high-pressure experiments and three-dimensional x-ray microtomographic imaging, we showed that the solidus temperature of a primitive (pyrolitic) mantle is as low as 3570 ± 200 kelvin at pressures expected near the boundary between the mantle and the outer core. Because the lowermost mantle is not globally molten, this provides an upper bound of the temperature at the core-mantle boundary (T(CMB)). Such remarkably low T(CMB) implies that the post-perovskite phase is present in wide areas of the lowermost mantle. The low T(CMB) also requires that the melting temperature of the outer core is depressed largely by impurities such as hydrogen.

An X-ray Raman spectrometer for EXAFS studies on minerals: bent Laue spectrometer with 20 keV X-rays

An X-ray Raman spectrometer for studies of local structures in minerals is discussed. Contrary to widely adopted back-scattering spectrometers using ≤10 keV X-rays, a spectrometer utilizing ~20 keV X-rays and a bent Laue analyzer is proposed. The 20 keV photons penetrate mineral samples much more deeply than 10 keV photons, so that high intensity is obtained owing to an enhancement of the scattering volume. Furthermore, a bent Laue analyzer provides a wide band-pass and a high reflectivity, leading to a much enhanced integrated intensity. A prototype spectrometer has been constructed and performance tests carried out. The oxygen K-edge in SiO(2) glass and crystal (α-quartz) has been measured with energy resolutions of 4 eV (EXAFS mode) and 1.3 eV (XANES mode). Unlike methods previously adopted, it is proposed to determine the pre-edge curve based on a theoretical Compton profile and a Monte Carlo multiple-scattering simulation before extracting EXAFS features. It is shown that the obtained EXAFS features are reproduced fairly well by a cluster model with a minimal set of fitting parameters. The spectrometer and the data processing proposed here are readily applicable to high-pressure studies.

Solid-liquid iron partitioning in Earth's deep mantle

Melting processes in the deep mantle have important implications for the origin of the deep-derived plumes believed to feed hotspot volcanoes such as those in Hawaii. They also provide insight into how the mantle has evolved, geochemically and dynamically, since the formation of Earth. Melt production in the shallow mantle is quite well understood, but deeper melting near the core-mantle boundary remains controversial. Modelling the dynamic behaviour of deep, partially molten mantle requires knowledge of the density contrast between solid and melt fractions. Although both positive and negative melt buoyancies can produce major chemical segregation between different geochemical reservoirs, each type of buoyancy yields drastically different geodynamical models. Ascent or descent of liquids in a partially molten deep mantle should contribute to surface volcanism or production of a deep magma ocean, respectively. We investigated phase relations in a partially molten chondritic-type material under deep-mantle conditions. Here we show that the iron partition coefficient between aluminium-bearing (Mg,Fe)SiO(3) perovskite and liquid is between 0.45 and 0.6, so iron is not as incompatible with deep-mantle minerals as has been reported previously. Calculated solid and melt density contrasts suggest that melt generated at the core-mantle boundary should be buoyant, and hence should segregate upwards. In the framework of the magma oceans induced by large meteoritic impacts on early Earth, our results imply that the magma crystallization should push the liquids towards the surface and form a deep solid residue depleted in incompatible elements.

Experimental and theoretical evidence for pressure-induced metallization in FeO with rocksalt-type structure

Electrical conductivity of FeO was measured up to 141 GPa and 2480 K in a laser-heated diamond-anvil cell. The results show that rock-salt (B1) type structured FeO metallizes at around 70 GPa and 1900 K without any structural phase transition. We computed fully self-consistently the electronic structure and the electrical conductivity of B1 FeO as a function of pressure and temperature, and found that although insulating as expected at ambient condition, B1 FeO metallizes at high temperatures, consistent with experiments. The observed metallization is related to spin crossover.