Spin and valence dependence of iron partitioning in Earth's deep mantle

Pressure stabilizes ferrous iron in bridgmanite under hydrous deep lower mantle conditions

Earth's lower mantle is a potential water reservoir. The physical and chemical properties of the region are in part controlled by the Fe3+/ΣFe ratio and total iron content in bridgmanite. However, the water effect on the chemistry of bridgmanite remains unclear. We carry out laser-heated diamond anvil cell experiments under hydrous conditions and observe dominant Fe2+ in bridgmanite (Mg, Fe)SiO3 above 105 GPa under the normal geotherm conditions corresponding to depth > 2300 km, whereas Fe3+-rich bridgmanite is obtained at lower pressures. We further observe FeO in coexistence with hydrous NiAs-type SiO2 under similar conditions, indicating that the stability of ferrous iron is a combined result of H2O effect and high pressure. The stability of ferrous iron in bridgmanite under hydrous conditions would provide an explanation for the nature of the low-shear-velocity anomalies in the deep lower mantle. In addition, entrainment from a hydrous dense layer may influence mantle plume dynamics and contribute to variations in the redox conditions of the mantle.

Iron silicate perovskite and postperovskite in the deep lower mantle

Ferromagnesian silicates are the dominant constituents of the Earth's mantle, which comprise more than 80% of our planet by volume. To interpret the low shear-velocity anomalies in the lower mantle, we need to construct a reliable transformation diagram of ferromagnesian silicates over a wide pressure-temperature (P-T) range. While MgSiO3 in the perovskite structure has been extensively studied due to its dominance on Earth, phase transformations of iron silicates under the lower mantle conditions remain unresolved. In this study, we have obtained an iron silicate phase in the perovskite (Pv) structure using synthetic fayalite (Fe2SiO4) as the starting material under P-T conditions of the lower mantle. Chemical analyses revealed an unexpectedly high Fe/Si ratio of 1.72(3) for the Pv phase in coexistence with metallic iron particles, indicating incorporation of about 25 mol% Fe2O3 in the Pv phase with an approximate chemical formula (Fe2+0.75Fe3+0.25)(Fe3+0.25Si0.75)O3. We further obtained an iron silicate phase in the postperovskite (PPv) structure above 95 GPa. The calculated curves of compressional (VP) and shear velocity (VS) of iron silicate Pv and PPv as a function of pressure are nearly parallel to those of MgSiO3, respectively. To the best of our knowledge, the iron silicate Pv and PPv are the densest phases among all the reported silicates stable at P-T conditions of the lower mantle. The high ferric iron content in the silicate phase and the spin-crossover of ferric iron at the Si-site above ~55 GPa should be taken into account in order to interpret the seismic observations. Our results would provide crucial information for constraining the geophysical and geochemical models of the lower mantle.

espm: A Python library for the simulation of STEM-EDXS datasets

We present two open-source Python packages: "electron spectro-microscopy" (espm) and "electron microscopy tables" (emtables). The espm software enables the simulation of scanning transmission electron microscopy energy-dispersive X-ray spectroscopy datacubes, based on user-defined chemical compositions and spatial abundance maps of constituent phases. The simulation process uses X-ray emission cross-sections generated via state-of-the-art calculations made with emtables. These tables are designed to be easily modifiable, either manually or using espm. The simulation framework is designed to test the application of decomposition algorithms for the analysis of STEM-EDX spectrum images with access to a known ground truth. We validate our approach using the case of a complex geology-related sample, comparing raw simulated and experimental datasets and the outputs of their non-negative matrix factorization. In addition to testing machine learning algorithms, our packages will also help experimental design, for instance, predicting dataset characteristics or establishing minimum counts needed to measure nanoscale features.

Compositional and thermal state of the lower mantle from joint 3D inversion with seismic tomography and mineral elasticity

The compositional and thermal state of Earth's mantle provides critical constraints on the origin, evolution, and dynamics of Earth. However, the chemical composition and thermal structure of the lower mantle are still poorly understood. Particularly, the nature and origin of the two large low-shear-velocity provinces (LLSVPs) in the lowermost mantle observed from seismological studies are still debated. In this study, we inverted for the 3D chemical composition and thermal state of the lower mantle based on seismic tomography and mineral elasticity data by employing a Markov chain Monte Carlo framework. The results show a silica-enriched lower mantle with a Mg/Si ratio less than ~1.16, lower than that of the pyrolitic upper mantle (Mg/Si = 1.3). The lateral temperature distributions can be described by a Gaussian distribution with a standard deviation (SD) of 120 to 140 K at 800 to 1,600 km and the SD increases to 250 K at 2,200 km depth. However, the lateral distribution in the lowermost mantle does not follow the Gaussian distribution. We found that the velocity heterogeneities in the upper lower mantle mainly result from thermal anomalies, while those in the lowermost mantle mainly result from compositional or phase variations. The LLSVPs have higher density at the base and lower density above the depth of ~2,700 km than the ambient mantle, respectively. The LLSVPs are found to have ~500 K higher temperature, higher Bridgmanite and iron content than the ambient mantle, supporting the hypothesis that the LLSVPs may originate from an ancient basal magma ocean formed in Earth's early history.

Structural transition and re-emergence of iron's total electron spin in (Mg,Fe)O at ultrahigh pressure

Fe-bearing MgO [(Mg_1−xFe_x)O] is considered a major constituent of terrestrial exoplanets. Crystallizing in the B1 structure in the Earth’s lower mantle, (Mg_1−xFe_x)O undergoes a high-spin (S = 2) to low-spin (S = 0) transition at ∼45 GPa, accompanied by anomalous changes of this mineral’s physical properties, while the intermediate-spin (S = 1) state has not been observed. In this work, we investigate (Mg_1−xFe_x)O (x ≤ 0.25) up to 1.8 TPa via first-principles calculations. Our calculations indicate that (Mg_1−xFe_x)O undergoes a simultaneous structural and spin transition at ∼0.6 TPa, from the B1 phase low-spin state to the B2 phase intermediate-spin state, with Fe’s total electron spin S re-emerging from 0 to 1 at ultrahigh pressure. Upon further compression, an intermediate-to-low spin transition occurs in the B2 phase. Depending on the Fe concentration (x), metal–insulator transition and rhombohedral distortions can also occur in the B2 phase. These results suggest that Fe and spin transition may affect planetary interiors over a vast pressure range.

Iron spin transition occurs at ultrahigh pressure. The total electron spin increases from 0 to 1 as the structural transition of (Mg,Fe)O occurs (~0.6 TPa) and drops back to 0 at higher pressure. Its effects on exoplanet interiors are anticipated.

Seismological expression of the iron spin crossover in ferropericlase in the Earth's lower mantle

The two most abundant minerals in the Earth’s lower mantle are bridgmanite and ferropericlase. The bulk modulus of ferropericlase (Fp) softens as iron d-electrons transition from a high-spin to low-spin state, affecting the seismic compressional velocity but not the shear velocity. Here, we identify a seismological expression of the iron spin crossover in fast regions associated with cold Fp-rich subducted oceanic lithosphere: the relative abundance of fast velocities in P- and S-wave tomography models diverges in the ~1,400-2,000 km depth range. This is consistent with a reduced temperature sensitivity of P-waves throughout the iron spin crossover. A similar signal is also found in seismically slow regions below ~1,800 km, consistent with broadening and deepening of the crossover at higher temperatures. The corresponding inflection in P-wave velocity is not yet observed in 1-D seismic profiles, suggesting that the lower mantle is composed of non-uniformly distributed thermochemical heterogeneities which dampen the global signature of the Fp spin crossover.

This study identifies the predicted seismic expression of the high-to-low iron spin crossover in the deep Earth mineral ferropericlase. A depth-dependent signal is detected in the fastest and slowest regions, related to lateral temperature variations, of several global seismic tomography models.

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).

Valence and spin states of iron are invisible in Earth's lower mantle

Heterogeneity in Earth's mantle is a record of chemical and dynamic processes over Earth's history. The geophysical signatures of heterogeneity can only be interpreted with quantitative constraints on effects of major elements such as iron on physical properties including density, compressibility, and electrical conductivity. However, deconvolution of the effects of multiple valence and spin states of iron in bridgmanite (Bdg), the most abundant mineral in the lower mantle, has been challenging. Here we show through a study of a ferric-iron-only (Mg_0.46Fe³⁺_0.53)(Si_0.49Fe³⁺_0.51)O₃ Bdg that Fe³⁺ in the octahedral site undergoes a spin transition between 43 and 53 GPa at 300 K. The resolved effects of the spin transition on density, bulk sound velocity, and electrical conductivity are smaller than previous estimations, consistent with the smooth depth profiles from geophysical observations. For likely mantle compositions, the valence state of iron has minor effects on density and sound velocities relative to major cation composition.

Stability of ferrous-iron-rich bridgmanite under reducing midmantle conditions

Our current understanding of the electronic state of iron in lower-mantle minerals leads to a considerable disagreement in bulk sound speed with seismic measurements if the lower mantle has the same composition as the upper mantle (pyrolite). In the modeling studies, the content and oxidation state of Fe in the minerals have been assumed to be constant throughout the lower mantle. Here, we report high-pressure experimental results in which Fe becomes dominantly Fe²⁺ in bridgmanite synthesized at 40-70 GPa and 2,000 K, while it is in mixed oxidation state (Fe³⁺/∑Fe = 60%) in the samples synthesized below and above the pressure range. Little Fe³⁺ in bridgmanite combined with the strong partitioning of Fe²⁺ into ferropericlase will alter the Fe content for these minerals at 1,100- to 1,700-km depths. Our calculations show that the change in iron content harmonizes the bulk sound speed of pyrolite with the seismic values in this region. Our experiments support no significant changes in bulk composition for most of the mantle, but possible changes in physical properties and processes (such as viscosity and mantle flow patterns) in the midmantle.