Search NASA⌕ Search

SEARCH · Search NASA

Results for “EARTH MANTLE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Origin of 182 W Anomalies in Ocean Island Basalts

Ocean island basalts (OIB) show variable 182 W deficits that have been attributed to either early differentiation of the mantle or core-mantle interaction. However, 182 W variations may also reflect nucleosynthetic isotope heterogeneity inherited from Earth's building material, which would be evident from correlated 182 W and 183 W anomalies. Some datasets for OIB indeed show hints for such correlated variations, meaning that a nucleosynthetic origin of W isotope anomalies in OIB cannot be excluded. We report high-precision W isotope data for OIB from Samoa and Hawaii, which confirm previously reported 182 W deficits for these samples, but also demonstrate that none of these samples have resolvable 183 W anomalies. These data therefore rule out a nucleosynthetic origin of the 182 W deficits in OIB, which most likely reflect the entrainment of either core material or an overabundance of late-accreted materials within OIB mantle sources. If these processes occurred over Earth's history, they may have also been responsible for shifting the 182 W composition of the bulk mantle to its modern-day value. We also report Mo isotope data for some Hawaiian OIB, which reveal no resolved nucleosynthetic Mo isotopic anomalies. This is consistent with inheritance of 182 W deficits in OIB from the addition of either core or late-accreted material, but only if these materials have a non-carbonaceous (NC) meteorite-like heritage. As such, these data rule out significant contributions of carbonaceous chondrite (CC)-like materials to either Earth's core or late accretion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The stability of subducted glaucophane with the Earth’s secular cooling

The blueschist to eclogite transition is one of the major geochemical–metamorphic processes typifying the subduction zone, which releases fluids triggering earthquakes and arc volcanism. Although glaucophane is an index hydrous mineral for the blueschist facies, its stability at mantle depths in diverse subduction regimes of contemporary and early Earth has not been experimentally determined. Here, we show that the maximum depth of glaucophane stability increases with decreasing thermal gradients of the subduction system. Along cold subduction geotherm, glaucophane remains stable down ca. 240 km depth, whereas it dehydrates and breaks down at as shallow as ca. 40 km depth under warm subduction geotherm or the Proterozoic tectonic setting. Our results imply that secular cooling of the Earth has extended the stability of glaucophane and consequently enabled the transportation of water into deeper interior of the Earth, suppressing arc magmatism, volcanism, and seismic activities along subduction zones.

58 GEOSCIENCES↗
DOE OSTI2026

Improved Regional Moment Tensor Inversion for Moderately Large Earthquakes in the Western United States Using a 3D Earth Model Based on Full Waveform Tomography

The nature of seismic sources for moderately large (moment magnitude, M w 5.0–6.5) events are commonly characterized by their moment tensor (MT) solutions and obtained by inversion of regional distance (200–1600 km) long‐period (20–50 s) waveforms. Regional MT estimates are often calculated from average plane‐layered, one‐dimensional (1D) velocity models. However, 1D model calculations can produce misfits in the arrival times and waveform shapes that introduce errors, particularly at longer distances or for shorter periods, which are necessary for analyzing lower magnitude events. Approximate Earth models (e.g., 1D) representing broad areas may be inadequate, particularly in the crust and uppermost mantle of tectonically complex regions. In this study, we show how a three‐dimensional (3D) Earth model obtained from full waveform inversion tomography can improve waveform fits and decrease phase errors. We developed a platform and workflow to perform routine 3D MT inversions and inverted MTs for 25 earthquakes in the western United States and seven nuclear explosions using an average 1D and a recent 3D Earth model, WUS256 (Rodgers et al., 2022). Using the 3D model improves waveform fits (variance reduction and phase time shifts) compared with the 1D model, and the 3D MT solutions are stable across large distances. This study shows that 3D models obtained from full waveform tomography can improve MTs and source characterization especially at far regional distances (>800 km).

Geosciences↗
DOE OSTI2020

Comprehensive geoneutrino analysis with Borexino

This paper presents a comprehensive geoneutrino measurement using the Borexino detector, located at Laboratori Nazionali del Gran Sasso (LNGS) in Italy. The analysis is the result of 3262.74 days of data between December 2007 and April 2019. The paper describes improved analysis techniques and optimized data selection, which includes enlarged fiducial volume and sophisticated cosmogenic veto. The reported exposure of ( 1.29 ± 0.05 ) × 10 32 protons × year represents an increase by a factor of two over a previous Borexino analysis reported in 2015. By observing 52.6 - 8.6 + 9.4 ( stat ) - 2.1 + 2.7 ( sys ) geoneutrinos (68% interval) from U 238 and Th 232 , a geoneutrino signal of 47.0 - 7.7 + 8.4 ( stat ) - 1.9 + 2.4 ( sys ) TNU with - 17.2 + 18.3 % total precision was obtained. This result assumes the same Th/U mass ratio as found in chondritic CI meteorites but compatible results were found when contributions from U 238 and Th 232 were both fit as free parameters. Antineutrino background from reactors is fit unconstrained and found compatible with the expectations. The null-hypothesis of observing a geoneutrino signal from the mantle is excluded at a 99.0% C.L. when exploiting detailed knowledge of the local crust near the experimental site. Measured mantle signal of 21.2 - 9.0 + 9.5 ( stat ) - 0.9 + 1.1 ( sys ) TNU corresponds to the production of a radiogenic heat of 24.6 - 10.4 + 11.1 TW (68% interval) from U 238 and Th 232 in the mantle. Assuming 18% contribution of K 40 in the mantle and 8.1 - 1.4 + 1.9 TW of total radiogenic heat of the lithosphere, the Borexino estimate of the total radiogenic heat of the Earth is 38.2 - 12.7 + 13.6 TW , which corresponds to the convective Urey ratio of 0.78 - 0.28 + 0.41 . These values are compatible with different geological predictions, however there is a ~ 2.4 σ tension with those Earth models which predict the lowest concentration of heat-producing elements in the mantle. In addition, by constraining the number of expected reactor antineutrino events, the existence of a hypothetical georeactor at the center of the Earth having power greater than 2.4 TW is excluded at 95% C.L. Particular attention is given to the description of all analysis details which should be of interest for the next generation of geoneutrino measurements using liquid scintillator detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

No detectable redox exchange between sulfur and iron during rapid cooling of basalts

The valence states of iron and sulfur in mantle-derived melts influence and respond to the composition of the mantle as well as the subsequent evolution of those melts in Earth's crust. Characterizing the valence states of iron and sulfur in silicate melts that have quenched to form glasses is therefore critical to understanding the composition of the mantle and the magmatic processes that take place before and during volcanic eruptions. Glasses quenched from silicate melts offer the opportunity to precisely measure sulfur and iron valence states, but whether the transition from melt to glass affects the valence state of sulfur and iron is uncertain. Here, we use the glassy margins of two pillow basalts, one from a mid-ocean ridge and one from a back-arc basin, as a natural experiment to test the effect of quench rate on sulfur and iron valence states. We report micro-X-ray Absorption Near Edge Structure (XANES) measurements of S 6+ /ΣS (S 6+ /[S 6+ +S 2- ]) and Fe 3+ /ΣFe (Fe 3+ /[Fe 2+ +Fe 3+ ]) along transects from the rapidly quenched rims of the pillows down into the slowly cooled crystal-rich interiors. The range of average quench rates estimated for our selected analysis areas range from 1.2 x 10 5 °C/sec at the pillow margin to 1.1 °C/sec in the interior and span quench rates experienced by natural and experimental glasses. On average, S6+/ΣS = 0.09 ± 0.01 and Fe 3+ /ΣFe = 0.147 ± 0.002 in the mid-ocean ridge sample and S 6+ /ΣS = 0.21 ± 0.01and Fe 3+ /ΣFe = 0.153 ± 0.004 in the back-arc sample, where the quoted uncertainty is the 1σ standard deviation of n= 9-16 analyses. We find that the S 6+ /ΣS and Fe 3+ /ΣFe values along transects in each sample show no systematic changes with depth, and by proxy, no systematic changes with quench rate. In each pillow, all measured S 6+ /ΣS and Fe 3+ /ΣFe ratios are within the absolute uncertainty of each technique. We propose that these observations are consistent with either sluggish reaction kinetics between sulfur and iron, or an insignificant temperature effect on sulfur-iron redox exchange equilibria. These observations demonstrate that sulfur and iron valence states in glasses with varying quench rates can be directly compared, and that changes in sulfur and iron valence during quench are below the limits of detection with XANES.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

High‐Pressure Phase Stability and Thermoelastic Properties of Iron Carbonitrides and Nitrogen in the Deep Earth

Abstract Iron‐dominant metallic phases are likely the primary hosts for nitrogen in the reduced deep Earth, hence the storage of nitrogen in the lower mantle and the core is governed by the behavior of the Fe‐N‐C system at high temperatures and pressures. In this study, phase transitions and thermoelastic properties of iron carbonitrides were investigated at high pressure‐temperature conditions by diamond anvil cell experiments and first‐principles calculations. Experimental data revealed no phase transition in ε ‐type Fe 4 (N 0.6 C 0.4 ) or Fe 7 (N 0.75 C 0.25 ) 3 up to 60 GPa at room temperature. At high temperature, Fe 7 (N 0.75 C 0.25 ) 3 transforms into the Fe 3 C‐type phase at ∼27 GPa, and then into the Fe 7 C 3 ‐type phase at ∼45 GPa, which is also corroborated by our theoretical calculations. We found that the phase stability of iron carbonitrides mainly depends on the N/C ratio, and the elastic properties of iron carbonitrides are dominantly affected by the Fe/(N+C) ratio. Iron carbonitrides with diverse structures may be the main host for nitrogen in the deep mantle. Some iron carbonitride inclusions in lower mantle diamonds could be the residue of the primordial mantle or originate from subducted nitrogen‐bearing materials, rather than iron‐enriched phases of the outer core. In addition, our experiments confirmed the existence of Fe 7 C 3 ‐type Fe 7 C 3 ‐Fe 7 N 3 solid solutions above 40 GPa. Fe 7 C 3 ‐type Fe 7 (C, N) 3 has comparable density and thermoelastic properties to its isostructural endmembers and may be a promising candidate constituent of the Earth's inner core.

Huang, Shengxuan↗
DOE OSTI2025

Sound Velocity of Hematite up to 13 GPa

We performed acoustic measurements on polycrystalline hematite using ultrasonic interferometric techniques in a Kawai-type multi-anvil apparatus. The elastic moduli and their pressure derivatives were determined by fitting the experimental data to third-order finite strain equations, yielding: 𝐾 𝑆,0 = 214.8 (20) GPa, 𝐺 0 = 92.0 (2) GPa, 𝐾$^{′}_{S,0}$ = 5.70 (1), and 𝐺$^{′}_{0}$ = 0.88 (4). Modeling of the velocity profiles of hematite along hot and cold slab geotherms shows that compressional (V P ) and shear (V S ) wave velocities are 22%–24% and 53%–55% lower, respectively, compared to the Preliminary Reference Earth Model at the top of the lower mantle. Furthermore, these findings likely offer a potential explanation for the low-velocity anomalies associated with slab graveyards at the base of the lower mantle.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High Sodium Solubility in Magnesiowüstite in Iron‐Rich Deep Lower Mantle

Abstract (Mg,Fe)O ferropericlase‐magnesiowüstite has been proposed to host the majority of Earth's sodium, but the mechanism and capacity for incorporating the alkali cation remain unclear. In this work, experiments in the laser‐heated diamond anvil cell and first‐principles calculations determine the solubility of sodium and favorability of sodium incorporation in iron‐rich magnesiowüstite relative to (Mg,Fe)SiO 3 bridgmanite. Reaction of Mg/(Mg + Fe) (Mg#) 55 and 28 olivine with NaCl at 33–128 GPa and 1600–3000 K produces iron‐rich magnesiowüstite containing several percent sodium, while iron‐rich bridgmanite contains little to no detectable sodium. In sodium‐saturated magnesiowüstite, sodium number [Na/(Na + Mg + Fe)] is 2–5 atomic percent at pressures below 60 GPa and drastically increases to 10–20 atomic percent at deep lower mantle pressures. For these two compositions, there is no significant dependence of the results on Mg#. Our calculations not only show consistent results with experiments but further indicate that such an increase in solubility and partitioning of Na into magnesiowüstite is driven by the spin transition in iron. These results provide fundamental constraints on the crystal chemistry of sodium at lower‐mantle conditions. If the sodium capacity of (Mg,Fe)O is not strongly dependent on Mg#, (Mg,Fe)O in the lower mantle may have the capacity to store the entire sodium budget of the Earth.

Geochemistry & Geophysics↗
DOE OSTI2021

Enhanced Visibility of Subduction Slabs by the Formation of Dense Hydrous Phase A

Abstract Phase A (Mg 7 Si 2 O 8 (OH) 6 ) is one of the important dense hydrous magnesium silicates in subducting slab, because it forms after the breakdown of antigorite serpentine and could be the dominant hydrous phase in the upper‐mantle deep slab for the water transportation into deep Earth. In this study, the compressional ( P ) and shear ( S ) wave velocities of phase A were measured at simultaneous pressure and temperature conditions up to 11 GPa and 1073 K. Combined with elastic properties of olivine and pyroxenes, we calculate the hydration effect on the velocities of harzburgite lithology in cold subduction zones throughout the depth ranges where phase A is thermodynamically stable. Our calculations suggest that the hydration increases both P ‐ and S ‐wave velocities of harzburgite; at ∼5 wt% hydration, its seismic detectability is enhanced by 1%–1.5% in velocity contrasts relative to its anhydrous counterpart.

Cai, Nao↗
DOE OSTI2020

Redox state of Earth’s magma ocean and its Venus-like early atmosphere

Exchange between a magma ocean and vapor produced Earth’s earliest atmosphere. Its speciation depends on the oxygen fugacity (fO 2 ) set by the Fe 3+ /Fe 2+ ratio of the magma ocean at its surface. Here, we establish the relationship between fO 2 and Fe 3+ /Fe 2+ in quenched liquids of silicate Earth-like composition at 2173 K and 1 bar. Mantle-derived rocks have Fe 3+ /(Fe 3+ +Fe 2+ ) = 0.037 ± 0.005, at which the magma ocean defines an fO 2 0.5 log units above the iron-wüstite buffer. At this fO 2 , the solubilities of H-C-N-O species in the magma ocean produce a CO-rich atmosphere. Cooling and condensation of H 2 O would have led to a prebiotic terrestrial atmosphere composed of CO 2 -N 2 , in proportions and at pressures akin to those observed on Venus. Present-day differences between Earth’s atmosphere and those of her planetary neighbors result from Earth’s heliocentric location and mass, which allowed geologically long-lived oceans, in-turn facilitating CO 2 drawdown and, eventually, the development of life.

58 GEOSCIENCES↗
DOE OSTI2024

Probing the state of hydrogen in δ − AlOOH at mantle conditions with machine learning potential

Hydrous and nominally anhydrous minerals are a fundamental class of solids of enormous significance to geophysics. They are the water carriers in the deep geological water cycle and impact structural, elastic, plastic, and thermodynamic properties and phase relations in Earth's forming aggregates (rocks). They play a critical role in the geochemical and geophysical processes that shape the planet. Their complexity has prevented predictive calculations of their properties, but progress in materials simulations ushered by machine-learning potentials is transforming this state of affairs. Here, we adopt a hybrid approach that combines deep learning potentials (DPs) with the strongly constrained and appropriately normed meta-generalized gradient approximation functional to simulate a prototypical hydrous system. We illustrate the success of this approach to simulate δ − AlOOH ( δ ) , a phase capable of transporting water down to near the core-mantle boundary of the Earth ( ∼ 2900 km depth and ∼ 135 GPa ) in subducting slabs. A high-throughput sampling of phase space using molecular dynamics simulations with DPs sheds light on the hydrogen-bond behavior and proton diffusion at geophysical conditions. These simulations provide a pathway for a deeper understanding of these crucial components that shape Earth's internal state. Published by the American Physical Society 2024

Luo, Chenxing (ORCID:0000000341166851)↗
DOE OSTI2020

Key problems of the four-dimensional Earth system

Compelling evidence indicates that the solid Earth consists of two physicochemically distinct zones separated radially in the middle of the lower mantle at ~1800 km depth. The inner zone is governed by pressure-induced physics and chemistry dramatically different from the conventional behavior in the outer zone. These differences generate large physical and chemical potentials between the two zones that provide fundamental driving forces for triggering major events in Earth’s history. One of the main chemical carriers between the two zones is H 2 O in hydrous minerals that subducts into the inner zone, releases hydrogen, and leaves oxygen to create superoxides and form oxygen-rich piles at the core–mantle boundary, resulting in localized net oxygen gain in the inner zone. Accumulation of oxygen-rich piles at the base of the mantle could eventually reach a supercritical level that triggers eruptions, injecting materials that cause chemical mantle convection, superplumes, large igneous provinces, extreme climate changes, atmospheric oxygen fluctuations, and mass extinctions. Interdisciplinary research will be the key for advancing a unified theory of the four-dimensional Earth system.

58 GEOSCIENCES↗
DOE OSTI2023

Equation of State and Spin Crossover of (Al, Fe)–Phase H

The transport of hydrogen into Earth's deep interior may have an impact on lower mantle dynamics as well as on the seismic signature of subducted material. Due to the stability of the hydrous phases δ-AlOOH (delta phase), MgSiO 2 (OH) 2 (phase H), and ε-FeOOH at high temperatures and pressures, their solid solutions may transport significant amounts of hydrogen as deep as the core-mantle boundary. Here we have constrained the equation of state, including the effects of a spin crossover in the Fe 3+ atoms, of (Al, Fe)-phase H: Al 0.84 Fe 3+ 0.07 Mg 0.02 Si 0.06 OOH, using powder X-ray diffraction measurements to 125 GPa, supported by synchrotron Mössbauer spectroscopy measurements on (Al, Fe)-phase H and δ-(Al, Fe)OOH. The changes in spin state of Fe 3+ in (Al, Fe)-phase H results in a significant decrease in bulk sound velocity and occurs over a different pressure range (48–62 GPa) compared with δ-(Al, Fe)OOH (32–40 GPa). Changes in axial compressibilities indicate a decrease in the compressibility of hydrogen bonds in (Al, Fe)-phase H near 30 GPa, which may be associated with hydrogen bond symmetrization. The formation of (Al, Fe)-phase H in subducted oceanic crust may contribute to scattering of seismic waves in the mid-lower mantle (~1,100–1,550 km). Accumulation of 1–4 wt.% (Al, Fe)-phase H could reproduce some of the seismic signatures of large, low seismic-velocity provinces. Our results suggest that changes in the electronic structure of phases in the (δ-AlOOH)-(MgSiO 2 (OH) 2 )-(ε-FeOOH) solid solution are sensitive to composition and that the presence of these phases in subducted oceanic crust could be seismically detectable throughout the lower mantle.

58 GEOSCIENCES↗
DOE OSTI2024

Xenon compatibility in magmatic processes: Hadean to current contexts

Xenon (Xe) behaviour in petrological processes, albeit essential to constrain mantle ingassing and degassing models, is elusive due to its volatile nature, and lack of direct investigation at the pressures (P) and temperatures (T) relevant to magma formation and crystallisation at depth. Xenon stands out amongst noble gases due to its unique reactivity with silicates of the lower crust and upper mantle, which could at least partially explain that published mineral/melt partitioning coefficients span up to six orders of magnitude. We report partition coefficients of Xe using in situ X-ray fluorescence at high P and T, and mass spectrometry analyses. Xenon is found to be moderately incompatible in anorthite-clinopyroxene mix in equilibrium with basalt (partition coefficient value of 0.16 ± 0.06), and compatible in olivine in equilibrium with basalt (partition coefficients in the range 88 ± 22 to 302 ± 46). While Xe is, thus, concentrated in basaltic melts coexisting with crystallising pyroxenes and feldspars, it is strongly retained in olivine at depth. Consequently, Xe originally contained in solid Earth has been preferentially retained at depth throughout Earth’s history, from the magma ocean stages to present day partial mantle melting processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Melting Curve and Phase Relations of Fe-Ni Alloys: Implications for the Earth's Core Composition

We discuss that nickel is the second most abundant element in the Earth's core. However, the properties of Fe-Ni alloys are still poorly constrained under planetary cores conditions, in particular concerning the effect of Ni on the melting curve of Fe. Here we show that Ni alloying up to 36 wt% does not affect the melting curve of Fe up to 100 GPa. However, Ni strongly modifies the hexagonal-closed-packed/face-centered-cubic (hcp/fcc) phase boundary, pushing the hcp/fcc/liquid triple point of Fe-20wt%Ni to higher pressures and temperatures. Our results allow constraining the triple point for Fe-10wt%Ni, a composition relevant for the Earth interior, and point out a decrease of the melting temperature at core-mantle boundary by 400 K with respect to pure Fe. A lower amount of light elements than previously predicted is thus required to reduce the crystallization temperature of core materials below that of a peridotitic lower mantle, in better agreement with geochemical observations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Capturing the interactions between ice sheets, sea level and the solid Earth on a range of timescales: a new “time window” algorithm

Retreat and advance of ice sheets perturb the gravitational field, solid surface and rotation of the Earth, leading to spatially variable sea-level changes over a range of timescales $\textit{O}$ (10 0–6 years), which in turn feed back onto ice-sheet dynamics. Coupled ice-sheet–sea-level models have been developed to capture the interactive processes between ice sheets, sea level and the solid Earth, but it is computationally challenging to capture short-term interactions $\textit{O}$ (10 0–2 years) precisely within longer $\textit{O}$ (10 3–6 years) simulations. The standard forward sea-level modelling algorithm assigns a uniform temporal resolution in the sea-level model, causing a quadratic increase in total CPU time with the total number of input ice history steps, which increases with either the length or temporal resolution of the simulation. In this study, we introduce a new “time window” algorithm for 1D pseudo-spectral sea-level models based on the normal mode method that enables users to define the temporal resolution at which the ice loading history is captured during different time intervals before the current simulation time. Utilizing the time window, we assign a fine temporal resolution $\textit{O}$ (10 0–2 years) for the period of ongoing and recent history of surface ice and ocean loading changes and a coarser temporal resolution $\textit{O}$ (10 3–6 years) for earlier periods in the simulation. This reduces the total CPU time and memory required per model time step while maintaining the precision of the model results. We explore the sensitivity of sea-level model results to the model temporal resolution and show how this sensitivity feeds back onto ice-sheet dynamics in coupled modelling. We apply the new algorithm to simulate sea-level changes in response to global ice-sheet evolution over two glacial cycles and the rapid collapse of marine sectors of the West Antarctic Ice Sheet in the coming centuries and provide appropriate time window profiles for each application. The time window algorithm reduces the total CPU time by ~ 50 % in each of these examples and changes the trend of the total CPU time increase from quadratic to linear. This improvement would increase with longer simulations than those considered here. Our algorithm also allows for coupling time intervals of annual temporal scale for coupled ice-sheet–sea-level modelling of regions such as West Antarctica that are characterized by rapid solid Earth response to ice changes due to the thin lithosphere and low mantle viscosities.

59 BASIC BIOLOGICAL SCIENCES↗