An official website of the United States government
Abstract Three-dimensional models of Earth’s seismic structure can be used to identify temperature-dependent phenomena, including mineralogical phase and spin transformations, that are obscured in 1-D spherical averages. Full-waveform tomography maps seismic wave-speeds inside the Earth in three dimensions, at a higher resolution than classical methods. By providing absolute wave speeds (rather than perturbations) and simultaneously constraining bulk and shear wave speeds over the same frequency range, it becomes feasible to distinguish variations in temperature from changes in composition or spin state. We present a quantitative joint interpretation of bulk and shear wave speeds in the lower mantle, using a recently published full-waveform tomography model. At all depths the diversity of wave speeds cannot be explained by an isochemical mantle. Between 1000 and 2500 km depth, hypothetical mantle models containing an electronic spin crossover in ferropericlase provide a significantly better fit to the wave-speed distributions, as well as more realistic temperatures and silica contents, than models without a spin crossover. Below 2500 km, wave speed distributions are explained by an enrichment in silica towards the core-mantle boundary. This silica enrichment may represent the fractionated remains of an ancient basal magma ocean.
more »
« less
Seismological expression of the iron spin crossover in ferropericlase in the Earth’s lower mantle
Abstract 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.
more »
« less
- Award ID(s):
- 1918126
- PAR ID:
- 10308202
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present LDA + U sc calculations of high-spin (HS) and low-spin (LS) states in ferropericlase (fp) with an iron concentration of 18.75%. The Hubbard parameter U is determined self-consistently with structures optimized at arbitrary pressures. We confirm a strong dependence of U on the pressure and spin state. Static calculations confirm that the antiferromagnetic configuration is more stable than the ferromagnetic one in the HS state, consistent with low-temperature measurements. Phonon calculations guarantee the dynamical stability of HS and LS states throughout the pressure range of the Earth mantle. Compression curves for HS and LS states agree well with experiments. Using a non-ideal mixing model for the HS to LS states solid solution, we obtain a crossover starting at ∼45 GPa at room temperature and considerably broader than previous results. The spin-crossover phase diagram is calculated, including vibrational, magnetic, electronic, and non-ideal HS–LS entropic contributions. Our results suggest the mixed-spin state predominates in fp in most of the lower mantle.more » « less
-
null (Ed.)Abstract The high-pressure phases of oxyhydroxides (δ-AlOOH, ε-FeOOH, and their solid solution), candidate components of subducted slabs, have wide stability fields, thus potentially influencing volatile circulation and dynamics in the Earth’s lower mantle. Here, we report the elastic wave velocities of δ-(Al,Fe)OOH (Fe/(Al + Fe) = 0.13, δ-Fe13) to 79 GPa, determined by nuclear resonant inelastic X-ray scattering. At pressures below 20 GPa, a softening of the phonon spectra is observed. With increasing pressure up to the Fe 3+ spin crossover (~ 45 GPa), the Debye sound velocity ( v D ) increases. At higher pressures, the low spin δ-Fe13 is characterized by a pressure-invariant v D . Using the equation of state for the same sample, the shear-, compressional-, and bulk-velocities ( v S , v P , and v Φ ) are calculated and extrapolated to deep mantle conditions. The obtained velocity data show that δ-(Al,Fe)OOH may cause low- v Φ and low- v P anomalies in the shallow lower mantle. At deeper depths, we find that this hydrous phase reproduces the anti-correlation between v S and v Φ reported for the large low seismic velocity provinces, thus serving as a potential seismic signature of hydrous circulation in the lower mantle.more » « less
-
Abstract Seismic tomography has demonstrated that the shear‐wave velocity is relatively high over a 3,000‐km wide region in the lowermost mantle beneath southern and eastern Asia. This seismic anomaly demarcates the current position of slab remnants that may have subducted in the Cretaceous. To further characterize the seismic structure at smaller scales, we measure 929 residual travel time differences (δt) between the phasesScSandSusing recordings of eight earthquakes beneath the Indian Ocean at stations from the Chinese Digital Seismic Network. We interpret variations of δtup to 10 s as due to horizontal shear‐velocity variations in D″ beneath northern India, Nepal, and southwestern China. The shear velocity can vary by as much as 7% over distances shorter than 300 km. Our observations provide additional observational evidence that compositional heterogeneity and possibly melt contribute to the seismic structure of the lower mantle characterized by long‐term subduction and mantle downwelling.more » « less
-
Abstract The dynamics of Earth's D″ layer at the base of the mantle plays an essential role in Earth's thermal and chemical evolution. Mantle convection in D″ is thought to result in seismic anisotropy; therefore, observations of anisotropy may be used to infer lowermost mantle flow. However, the connections between mantle flow and seismic anisotropy in D″ remain ambiguous. Here, we calculate the present‐day mantle flow field in D″ using 3D global geodynamic models. We then compute strain, a measure of deformation, outside the two large‐low velocity provinces (LLVPs) and compare the distribution of strain with previous observations of anisotropy. We find that, on a global scale, D″ materials are advected toward the LLVPs. The strains of D″ materials generally increase with time along their paths toward the LLVPs and toward deeper depths, but regions far from LLVPs may develop relative high strain as well. Materials in D″ outside the LLVPs mostly undergo lateral stretching, with the stretching direction often aligning with mantle flow direction, especially in fast flow regions. In most models, the depth‐averaged strain in D″ is >0.5 outside the LLVPs, consistent with widespread observations of seismic anisotropy. Flow directions inferred from anisotropy observations often (but not always) align with predictions from geodynamic modeling calculations.more » « less
