Search NASA⌕ Search

SEARCH · Search NASA

Results for “MAGMA”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The ultralow viscosity of volatile-rich kimberlite magma: Implications for the water content of primitive kimberlite melts

The eruption of deeply sourced kimberlite magma offers the fastest route to bring deep-seated volatiles back to the Earth’s surface. However, the viscosity of kimberlite magma, a factor governing its migration and eruption dynamics within Earth, remains poorly constrained. We conducted synchrotron in situ falling sphere viscometry experiments to examine kimberlite magma with different volatile contents (0 to 5 wt % H 2 O and 2 to 8 wt % CO 2 ) under high pressure-temperature conditions. The results reveal that the viscosity of volatile-rich kimberlite magma is ~1 to 2 orders lower than that of mid-ocean ridge basalt (MORB) and comparable to the ultramobile pure carbonate melt. Using the measured viscosity values, we simulated the ascent and eruption process of kimberlite magma. We found that a minimum content of ~0.5 wt % water in the primitive magma is necessary to allow the ultrafast eruption process of kimberlite, thereby enabling the preservation of diamonds and high-pressure mineral inclusions transported by the magma.

58 GEOSCIENCES↗
DOE OSTI2025

A sticky situation: The heterogeneity and implications of magma mush unlocking

Magma, partially molten rock below Earth's surface, can be made up of both fluid (melt) and particles (crystals). The relative fraction of the fluid and particles matters for the magma's mobility and ability to reach Earth's surface. Here we conduct physical laboratory experiments to examine the mechanisms through which a magma becomes mobile and the speed at which pure melt and infiltrate a mixture of fluids and particles. The results have implications for understanding magma mobility and the mechanisms that move magma to Earth's surface. Poster presented at the 2025 International Association of Volcanology and Chemistry of Earth's Interior Meeting (IAVCEI).

multiphase flows↗
DOE OSTI2026

Electrical conductivities of (Mg,Fe)O at extreme pressures and implications for planetary magma oceans

During planet formation, planets undergo many impacts that can generate magma oceans. When these crystallize, part of the magma densifies via iron enrichment and migrates to the core–mantle boundary, forming an iron-rich basal magma ocean (BMO). The BMO could generate a dynamo in early Earth and super-Earths if the electrical conductivity of the BMO, which is thought to be sensitive to its Fe content, is sufficiently high. To test this hypothesis, here we conduct laser-driven shock experiments on ferropericlase (Mg x ,Fe 1−x )O (0.95 ≤ x ≤ 1) as an Fe-rich BMO analogue, perform density functional theory molecular dynamics simulations on MgO and calculate the long-term evolution of super-Earths. We find that the d.c. conductivities of MgO and (Mg,Fe)O are indistinguishable between 467 GPa and 1,400 GPa, despite previous predictions. Here, we predict that super-Earths larger than 3–6 Earth masses can produce BMO-driven dynamos that are almost one order of magnitude stronger than core-driven dynamos for several billion years.

Exoplanets↗
DOE OSTI2025

Completion of lunar magma ocean solidification at 4.43 Ga

Crystallization of the lunar magma ocean yielded a chemically unique liquid residuum named KREEP. This component is expressed as a large patch on the near side of the Moon and a possible smaller patch in the northwest portion of the Moon’s South Pole-Aitken basin on the far side. Thermal models estimate that the crystallization of the lunar magma ocean (LMO) could have spanned from 10 and 200 My, while studies of radioactive decay systems have yielded inconsistent ages for the completion of LMO crystallization covering over 160 My. Here, we show that the Moon achieved >99% crystallization at 4,429 ± 76 Ma, indicating a lunar formation age of ~4,450 Ma or possibly older. Using the 176 Lu– 176 Hf decay system (t 1/2 = 37 Gy), we found that the initial 176 Hf/ 177 Hf ratios of lunar zircons with varied U–Pb ages are consistent with their crystallization from a KREEP-rich reservoir with a consistently low 176 Lu/ 177 Hf ratio of 0.0167 that emerged ~140 My after solar system formation. The previously proposed younger model age of ~4.33 Ga for the source of mare basalts (240 My after solar system formation) might reflect the timing of a large impact. Our results demonstrate that lunar magma ocean crystallization took place while the Moon was still battered by planetary embryos and planetesimals leftover from the main stage of planetary accretion. The study of Lu–Hf model ages for samples brought back from the South Pole-Aitken basin will help to assess the lateral continuity of KREEP and further understand its significance in the early history of the Moon.

Science & Technology - Other Topics↗
DOE OSTI2025

Shock compression of FeOOH and implications for iron-water interactions in super-earth magma oceans

Iron(Fe)-water reactions in a magma ocean can influence water storage and density of planets. These reactions can form Fe-O-H phases, whose density, melting, and electronic properties at planetary interior conditions are important for informing planetary models. Here, we study natural goethite (α-FeOOH) that is shock-compressed along its principal Hugoniot. Analysis of our velocity interferometer system for any reflector (VISAR) results extends the equation of state to over 800 GPa. X-ray diffraction and VISAR reflectivity results indicate the onset of melting occurs at ~95 GPa with complete melting by 166 GPa, which may be relevant to low seismic velocity anomalies observed above the core-mantle boundary. Analysis of X-ray emission spectroscopy results up to 285 GPa shows the spin crossover of Fe, with dominantly low spin Fe above ~265 GPa in the melt, supporting formation of dense basal magma oceans in terrestrial planets. Using our measured FeOOH densities, we model planetary interiors up to 10 Earth masses. Assuming FeOOH forms via iron-water reactions, the radius decreases by up to 28%, while the density increases by up to 165% compared to the unreacted cases, providing an avenue to investigate water storage and evolution in super-Earths and sub-Neptunes.

Zhang, Yanyao [Stanford Univ., CA (United States)]↗
DOE OSTI2025

Building wet planets through high-pressure magma–hydrogen reactions

Close-in transiting sub-Neptunes are abundant in our Galaxy. Planetary interior models based on their observed radius–mass relationship suggest that sub-Neptunes contain a discernible amount of either hydrogen (dry planets) or water (wet planets) blanketing a core composed of rocks and metal. Water-rich sub-Neptunes have been believed to form farther from the star and then migrate inwards to their present orbits. Here we report experimental evidence of reactions between warm, dense hydrogen fluid and silicate melt that release silicon from the magma to form alloys and hydrides at high pressures. We found that oxygen liberated from the silicate melt reacts with hydrogen, producing an appreciable amount of water up to a few tens of weight per cent, which is much greater than previously predicted based on low-pressure ideal gas extrapolation. Consequently, these reactions can generate a spectrum of water contents in hydrogen-rich planets, with the potential to reach water-rich compositions for some sub-Neptunes, implying an evolutionary relationship between hydrogen-rich and water-rich planets. Therefore, detection of a large amount of water in exoplanet atmospheres may not be the optimal evidence for planet migration in the protoplanetary disk, calling into question the assumed link between composition and planet formation location.

Horn, H. W. [Arizona State University, Tempe, AZ (↗
DOE OSTI2025

Temporal relationships among lunar crustal rocks

Temporal relationships among the three most common suites of lunar crustal rocks have been investigated by obtaining new high precision ages on Felsic/Alkali-suite Quartz monzodiorite Clast B from breccia 15405 and Magnesian-suite norite 78235/6/8/55/56 and comparing them to previously dated ferroan anorthosite sample 60025. The weighted average age of 4337.19 ± 0.49 Ma of 15405 Clast B is defined by zircon U-Pb and Pb-Pb ages as well as mineral isochron Sm-Nd and Nd-Nd ages. It is identical to the weighted average age for Apollo 17 norite 78235/6/8/55/56 of 4334.1 ± 3.5 Ma which is defined by Pb-Pb ages measured on baddeleyites in this investigation and less precise Pb-Pb and Sm-Nd ages reported in the literature. Both ages are ∼ 25 Ma younger than the weighted average of Sm-Nd and Pb-Pb ages reported in the literature on ferroan anorthosite 60025 of 4359.3 ± 2.3 Ma. The fact that ages of all three samples are defined by multiple U-Pb, Pb-Pb, Sm-Nd, and 142 Nd- 143 Nd chronometers provide confidence that they record the igneous crystallization history of the samples and do not represent disturbances or mixing lines with no temporal significance. Here, the extent to which these three ages represent broader scale magmatism is difficult to evaluate. Nevertheless, the age defined for 15405 Clast B, 78235/6/8/55/56, and 60025 are contemporaneous with the peak of ages observed in detrital zircons from the Apollo 12, 14, 15, and 17 landing sites (4340 ± 20 Ma), a Mg-suite Sm-Nd whole rock isochron defined by samples from Apollo 14, 15, 16, and 17 landing sites (4348 ± 25 Ma), and a Ferroan Anorthosite-suite Sm-Nd whole rock isochron defined by samples from the Apollo 15 and 16 landing sites (4354 ± 29 Ma). This implies that Ferroan Anorthosite-suite magmatism is temporally distinct and earlier than magmatism associated with the Mg-suite and the Felsic/Alkali-suite, as predicted by the lunar magma ocean model of lunar differentiation. The short 35 ± 10 Ma interval between primary ferroan anorthosite magmatism and secondary magmatism suggests that the lunar crust formed over a limited period of time. Although heat from decay of long-lived isotopes, large impacts, tidal heating associated with interactions between the Earth and Moon, and density driven overturn of the magma ocean have all been invoked to explain production of ancient secondary crustal magmatism, only tidal heating and cumulate overturn are consistent with the apparent short duration of secondary crustal magmatism and the great depth of crystallization implied for some Mg-suite samples. The initial ε 143 Nd values derived from the 15405 Clast B and 78238 Mg-suite norite isochrons, as well as a Mg-suite whole rock isochron are −0.23 ± 0.11, −0.27 ± 0.74, and −0.25 ± 0.09, respectively. They are identical within uncertainty indicating that Mg-suite and Felsic/Alkali-suite magmas were derived from materials that had the same time averaged Sm/Nd ratios since the formation of the solar system. This, combined with the contemporaneous nature of 15405 Clast B and 78235/6/8/55/56 Mg-suite norite, is consistent with evolution of both samples, and likely both magma suites, from a common source through closed system fractional crystallization or partial melting processes.

Chronology↗
DOE OSTI2025

Apatite geochemistry as a tool for understanding the petrogenesis of layered mafic-ultramafic rocks in the Bushveld Complex, South Africa

The sources of the magmas that formed the Rustenburg Layered Suite of the Bushveld Complex in South Africa remain debated, despite decades of research. Vertical and lateral variation in bulk rock and mineral separate Sr-Nd isotopic compositions, which generally indicate enriched sources, demonstrate that the layered sequence was formed by the emplacement of multiple batches of magma, crucially resulting in episodes of PGE-Cr-V mineralisation. The Lu-Hf isotope compositions of zircon are, however, at odds with the bulk rock Sr-Nd isotopic heterogeneity as they show near homogeneous compositions throughout the layered sequence (εHf (2.06 Ga) =−8). This lack of variation in Hf isotope composition has been attributed to deep, continental lithospheric mantle-related and/or crustal contamination of plume-derived Bushveld magmas. In this study, we analysed the major, trace element and Sr-Nd isotope geochemistry of apatite in the Rustenburg Layered Suite. Apatite occurs as an intercumulus mineral in the lowermost regions and a cumulus mineral in the uppermost regions of the layered sequence and can therefore be used to test existing models for the isotopic disequilibrium between bulk rock Sr-Nd and zircon Hf isotopic compositions. Apatite is largely chlorapatite in the lowermost regions and fluorapatite in the uppermost regions of the layered sequence. The Merensky Reef is unusual in that it comprises both chlorapatite and fluorapatite. Apatite throughout the layered sequence is generally unzoned and shows no evidence of late-stage alteration. Trace element data show that apatite is enriched in L/HREE, with common negative Eu-Sr anomalies. These trace element signatures are consistent with a magmatic origin for the apatite grains, with prior, or concurrent, plagioclase crystallization from the same melt. Variability in in situ Sr and Nd isotope compositions of apatite is recorded throughout the layered sequence with εNd (2.06 Ga) compositions varying between −2.5 and − 10.2 and initial 87 Sr/ 86 Sr compositions varying between 0.7079 and 0.7103 (for the Marikana dikes only). The variability in Sr-Nd isotope compositions of apatite is consistent with the bulk rock (and mineral separate) variation in Sr-Nd isotope compositions, suggesting apatite preserves primary magmatic compositions in the Rustenburg Layered Suite.

Apatite↗
DOE OSTI2024

Formation of a Large Cold Groundwater Mantle Helium Anomaly and High Temperature Geothermal Resources in Response to Bimodal Magmatism Near Roosevelt Hot Springs and Utah FORGE, Milford Valley, Southwest Utah

Abstract A large mantle helium anomaly and separate domains of high heat flow are the predominant manifestations of bimodal magmatic activity in the Milford valley. The mantle helium anomaly (1.9–2.6 R/Ra) covers 270 km 2 and is subdivided into two separated domains: a cold shallow groundwater regime and high temperature hydrothermal activity. The zone of anomalous heat flow covers >100 km 2 and is also subdivided into two adjacent domains, comprising hydrothermal activity at Roosevelt Hot Springs (RHS) (3–7 W/m 2 ) and conductive heat flow (100–180 mW/m 2 ). While the transfer of heat and mantle helium at RHS are coupled, heat and helium transfer are decoupled in the adjacent cold groundwater regime to the west. Both the mantle helium and geothermal anomalies are attributed to recent mafic‐felsic magmatic intrusions of >400 km 3 , however, the absence of volcanic eruptions <500,000 years indicates magmas stall before rising to shallow crustal level <10 km depth. Deep level magmatism produces a felsic composition melt, which is inferred to be responsible for the widespread and near uniform range of diluted mantle helium values. A thick and impermeable mass of crystalline granitic basement rock at the mid‐crustal level divides the ascent of mantle helium into separate flow paths. It may also impede the rise of buoyant magma trapping thermal energy that facilitates partial melting, slow cooling, and development of a thick thermal aureole. Partitioning of convective and conductive thermal regimes and independent flow paths supplying deeply derived helium characterize the development of a large long‐lived magma‐related geothermal system.

Geochemistry & Geophysics↗
DOE OSTI2024

Formation and interior evolution of Mars

The formation and earliest evolution of Mars involved core formation followed by magma ocean solidification and primordial crust formation. Here, in this study, we review the formation and early interior evolution of Mars as understood through meteorite studies. Martian meteorites exhibit ages that span the entire history of the Solar System, but most derive from compositionally distinct sources that formed early in the history of Mars. Geochemical and isotopic studies of Martian meteorites have provided critical information on Mars’ building blocks, bulk composition, and its differentiation into a core, mantle and crust. Radiogenic isotope systems provide vital insights into the petrogenesis of Martian meteorite source regions and into the nature and timescales of primordial differentiation on Mars. Application of short-lived isotope chronometers to Martian meteorites indicates that the Martian core formed early within <10 million years (Ma), while magma ocean differentiation and crust formation occurred soon thereafter, within the first tens of Ma of Solar System history. Bulk silicate Mars is variably depleted in moderately volatile elements relative to chondrites, but to a lesser extent than the Earth. The nature and evolution of more volatile elements and volatile species on Mars (e.g., H 2 O, C, S, halogens) remains an active area of investigation, but overall suggest that Mars is more enriched in volatile elements than the Moon but less than the Earth. While geochemical exploration of Mars coupled with new discoveries of Martian meteorites has provided a wealth of information, the compositional and lithologic diversity on Mars significantly exceeds that seen in the suites of Martian meteorites currently studied. This implies that our current perspective on the planet's evolution is very likely biased, highlighting the need for future sample return missions to Mars.

Kruijer, Thomas S.↗
DOE OSTI2024

Temporal covariation of island arc Sr isotopes and seawater chemistry over the past 2 billion years

The chemical compositions of island arc basalts (IAB) reflect contributions from the mantle as well as fluids and melts from the subducting slab. Addition of radiogenic seawater Sr to oceanic crust through hydrothermal alteration and subsequent subduction is often invoked to explain elevated 87 Sr/ 86 Sr signatures in modern IAB. However, changes in the 87 Sr/ 86 Sr of island arc magmatic rocks through time has not been investigated, limiting our understanding of the factors influencing the Sr budgets of arcs throughout Earth’s history. To address this, we compiled 87 Sr/ 86 Sr values from island arc magmatic rocks ranging in age from modern to Paleoproterozoic, only including data from island arc localities that best preserve initial magmatic 87 Sr/ 86 Sr. Median initial 87 Sr/ 86 Sr values are consistently elevated compared to depleted mantle 87 Sr/ 86 Sr over this period, indicating persistent enrichment in radiogenic Sr in island arcs. Moreover, the elevation in island arc 87 Sr/ 86 Sr relative to the depleted mantle is variable. A notable rise in island arc 87 Sr/ 86 Sr during the late Neoproterozoic coincides with a steep increase in seawater 87 Sr/ 86 Sr and Sr concentration. To investigate this potential connectivity, we modeled the 87 Sr/ 86 Sr of island arc magmas between 0 and 830 Ma with inputs of depleted mantle 87 Sr/ 86 Sr, seawater 87 Sr/ 86 Sr, and seawater Sr concentration. The model reproduces the overall trajectory of the compiled data. We interpret the observed temporal variation in island arc 87 Sr/ 86 Sr values and its close association with fluctuations in seawater chemistry as evidence that changes in marine geochemistry have strongly influenced the Sr isotopic record of island arc magmas over time.

Science & Technology - Other Topics↗
DOE OSTI2024

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The structure of high-Mg alkali-bearing aluminosilicate glasses investigated in situ at ambient and high pressure by multi-angle energy dispersive X-ray diffraction and infrared microspectroscopy

The structural properties of synthetic high-Mg alkali-bearing aluminosilicate glasses analogues of natural picritic-to-komatiitic magmas were investigated in situ by multiangle energy dispersive X-ray diffraction at 2.1 GPa and ambient pressure and by Fourier Transform infrared spectroscopy up to 5.4 GPa in a cycle of compression and decompression experiments. Our results show that the intermediate range ordering of the glass structure at 2.1 GPa is 3.14 Å, increasing to 3.19 Å when decompressed. The local structure shows T-O lengths of 1.66 Å (2.1 GPa) and 1.65 Å (ambient pressure), T-T distances of 3.19 Å at high pressure, which lengthen to 3.21 Å at ambient pressure, causing the T-O-T angle of 147° determined at 2.1 GPa to widen to 154° upon decompression. The deconvoluted infrared spectra result in the presence of Q 1 , Q 2 , Q 3 populations in the aluminosilicate spectral region, whose proportions remain relatively unchanged up to 5.4 GPa. The structural response of the investigated glasses to cold-compression does not involve changes in polymerization, but rather a shrinking and compaction of the structure as evidenced by the Qn species shifting to higher wavenumbers as a function of pressure. The structural properties determined from X-ray diffraction for this glass composition are discussed together with those of glasses emerging from previous studies to highlight a compositional dependence mainly dictated by the amount of SiO 2 and Al 2 O 3 .

glass structure↗
DOE OSTI2026

Chronology of ungrouped achondrites Northwest Africa 11119 and 7325/8486: implications for early evolution of a heterogeneous crust on a differentiated planetesimal

The chronology of ungrouped achondrites provides key insights into the timeline of igneous activity in the early Solar System. Two ungrouped achondrites – Northwest Africa (NWA) 11119 (andesite-dacite) and NWA 8486 and its pair NWA 7325 (olivine gabbro) – are the focus of this study, where their chronologies are investigated using multiple high-resolution techniques. Here we report the lead-lead (207Pb-206Pb) and manganese-chromium (53Mn-53Cr) systematics of NWA 11119, as well as the 207Pb-206Pb systematics of NWA 8486 alongside aluminum-magnesium (26Al-26Mg) systematics for NWA 7325. The U-corrected 207Pb-206Pb ages of NWA 11119 and the combined ages of NWA 7325/8486 are 4566.4 ± 0.8 Ma and 4563.8 ± 1.9 Ma, respectively. Additionally, we report the 53Mn-53Cr age of NWA 11119 to be 4564.4 ± 2.5 Ma and the 26Al-26Mg age of NWA 7325/8486 to be 4563.1 ± 0.3 Ma. The formation of NWA 11119 requires partial melting of a (likely chondritic) source reservoir leading to eruption of Si rich, alkali depleted magmas, while NWA 7325/8486 likely formed from a chemically fractionated reservoir with superchondritic Al/Mg. The clear geochemical and isotopic similarities of these achondrites, combined with the chronology reported here, is suggestive of formation of these two ungrouped achondrites on a common parent body which likely formed in the inner Solar System and experienced early differentiation under reducing conditions. If these achondrites did share a parent body, it would suggest that primary asteroids in the early Solar System commonly produced mineralogically and geochemically heterogeneous crusts. While mineralogical and geochemical heterogeneity is known to exist in the (mostly mafic) crust of asteroid Vesta, our findings show that even more significant crustal heterogeneity (representing felsic and mafic compositions) may exist on other asteroids.

Cartwright, Julia [University of Leicester, UK]↗
DOE OSTI2024

The Global Biogeochemical Cycle of Rhenium

Here, this paper is the first comprehensive synthesis of what is currently known about the different natural and anthropogenic fluxes of rhenium (Re) on Earth's surface. We highlight the significant role of anthropogenic mobilization of Re, which is an important consideration in utilizing Re in the context of a biogeochemical tracer or proxy. The largest natural flux of Re derives from chemical weathering and riverine transport to the ocean (dissolved = 62 × 10 6 g yr -1 and particulate = 5 × 10 6 g yr -1 ). This review reports a new global average [Re] of 16 ± 2 pmol L -1 , or 10 ± 1 pmol L -1 for the inferred pre-anthropogenic concentration without human impact, for rivers draining to the ocean. Human activity via mining (including secondary mobilization), coal combustion, and petroleum combustion mobilize approximately 560 × 10 6 g yr -1 Re, which is more than any natural flux of Re. There are several poorly constrained fluxes of Re that merit further research, including: submarine groundwater discharge, precipitation (terrestrial and oceanic), magma degassing, and hydrothermal activity. The mechanisms and the main host phases responsible for releasing (sources) or sequestrating (sinks) these fluxes remain poorly understood. This study also highlights the use of dissolved [Re] concentrations as a tracer of oxidation of petrogenic organic carbon, and stable Re isotopes as proxies for changes in global redox conditions.

58 GEOSCIENCES↗
DOE OSTI2024

Oxidized Sulfur Species in Slab Fluids as a Source of Enriched Sulfur Isotope Signatures in Arcs

Recycling of oxidized sulfur from subducting slabs to the mantle wedge provides simultaneous explanations for the elevated oxygen fugacity (fO 2 ) in subduction zones, their high hydrothermal and magmatic sulfur outputs, and the enriched sulfur isotopic signatures (i.e., δ 34 S > 0‰) of these outputs. However, a quantitative understanding of the abundance and speciation of sulfur in slab fluids consistent with high pressure experiments is lacking. Here we analyze published experimental data for anhydrite solubility in H 2 O-NaCl solutions to calibrate a high-pressure aqueous speciation model of sulfur within the framework of the deep earth water model. We characterize aqueous complexes, required to account for the high experimental anhydrite solubilities. We then use this framework to predict the speciation and solubility of sulfur in chemically complex fluids in equilibrium with model subducting mafic and ultramafic lithologies, from 2 to 3 GPa and 400 to 800°C at log fO 2 from FMQ-2 to FMQ+4. We show that sulfate complexes of calcium and sodium markedly enhance the stability of sulfate in moderately oxidized fluids in equilibrium with pyrite at fO 2 conditions of FMQ+1 to +2, causing large sulfur isotope fractionations up to 10‰ in the fluid relative to the slab. Such fluids could impart oxidized, sulfur-rich and high δ 34 S signatures to the mantle wedge that are ultimately transferred to arc magmas, without the need to invoke 34 S-rich subducted lithologies.

58 GEOSCIENCES↗
DOE OSTI2025

Sound Velocity and Grüneisen Parameter in Shock‐Melted Silica at Deep Earth Pressures

Silica is a primary component of rocky planet interiors and its melt properties are important for understanding planetary formation and differentiation, magma oceans, and the deep mantle. Although well understood in the solid state, the high-pressure behavior of liquid silica is poorly constrained at lower mantle pressures. Using laser interferometry to measure shock wave profiles, we report measured stress-density states and longitudinal sound speeds in shock-synthesized stishovite, from fused silica staring material, and across the solid-liquid phase boundary up to 154 GPa. Our results constrain completion of melt at 80 GPa and show that at pressures relevant to the deep mantles of Earth-sized, rocky planets, the Grüneisen parameter for liquid silica increases with compression. This finding is consistent with a continuous increase in Si-O coordination above six for liquid silica at core-mantle boundary relevant pressures and temperatures.

Ocampo, I. K. [Princeton University, NJ (United St↗