Search NASA⌕ Search

SEARCH · Search NASA

Results for “magmatic”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High pressure experimental investigation of clinopyroxene dissolution in a K-basaltic melt

Dissolution of clinopyroxene (cpx) in a K-basaltic melt from the Campi Flegrei Volcanic District (Italy) has been investigated through dissolution and dissolution-crystallization experiments at pressure of 0.8 GPa, superliquidus temperature of 1350 °C, and dwell times between 0.5 and 1 h. The obtained dissolution rates range from 7.9·10 –6 cm s –1 to 6.1·10 –6 cm s –1 as a function of dwell time. In the dissolution-crystallization experiment (1300 °C; 0.8 GPa; 2 h), the formation of overgrowth rims accompanied by new cpx crystals suggests that the injection of recharge magmas in basaltic reservoir may lead to inverse or oscillatory zonation. The interaction between cpx and K-basaltic melt at ~1570 °C was studied by in situ radiography using synchrotron radiation combined with the Paris-Edinburgh press. Here, this resulted in cpx resorption to occur depending on the temperature conditions with respect to the liquidus temperature of the cpx (T cpxL ). The calculated cpx dissolution rates are ~5·10 –3 cm s –1 at T ≤ T cpxL and ~ 3·10 –2 cm s –1 at T ≥ T cpxL . The role of crystal dissolution in the estimation of magma residence times has been also tested for a natural magmatic system by interpolating the dissolution rates (~10 –5 –10 –6 cm s –1 ) with the textural data of cpx phenocrysts from the Agnano-Monte Spina pyroclastic deposit at Campi Flegrei caldera (Campanian region, Italy). Results from calculations indicate that the time required for partial or complete resorption of phenocrysts varies from ~0.5 to ~40 h, and that the effect of crystal dissolution may be relevant to estimate magma residence times whether significant dissolution occurs during magma mixing processes.

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

U-Pb zircon age constraints on the earliest eruptions of the Deccan Large Igneous Province, Malwa Plateau, India

Climate instability driven by emission of volatiles during emplacement of large igneous provinces (LIPs) is frequently invoked as a potential cause of mass extinctions. However, documenting this process in the geologic record requires a holistic understanding of eruption rates, the location of eruptive centers, and potential sources of climate-changing volatiles. We present new chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb zircon geochronology from Malwa Plateau basalts on the northern margin of the Deccan LIP, India. These basalts have been previously interpreted as either an extension of the province's main volcanic stratigraphy or as an independent eruptive center active up to millions of years prior to the main eruptive phase. Our data instead demonstrate that the lower Malwa Plateau basalts are temporally correlative with the first pulse of Deccan volcanism and provide new constraints on its initiation and duration. Paleomagnetic data further indicate that upper Malwa Plateau basalts may be age-equivalent to the second, third, and fourth pulses of Deccan volcanism. The relative thicknesses of age-equivalent packages of basalt are consistent with eruption of the Deccan LIP from a southward-migrating eruptive center. The first eruptive pulse is coeval with a ~200 kyr Late Maastrichtian warming event preserved globally in contemporaneous stratigraphic sections. We propose that the first pulse of Deccan magmatism was more voluminous in the north, where it erupted through organic-rich sedimentary rocks of the Narmada-Tapti rift basin. Thermal metamorphism of these sediments could have been a source of sufficient CO 2 to drive the Late Maastrichtian warming event, which if true would reconcile the apparent dampened warming signals associated with later Deccan eruptive pulses.

58 GEOSCIENCES↗
DOE OSTI2022

Sound velocity and compressibility of melts along the hedenbergite (CaFeSi2O6)-diopside (CaMgSi2O6) join at high pressure : Implications for stability and seismic signature of Fe-rich melts in the mantle

Iron-rich silicate melts play an important role in the magmatic history of the Earth and the 16 Moon. However, their elastic properties at high pressures, especially the sound velocity, are poorly 17 understood. Here we determined the ultrasonic sound velocity for the first time of a hedenbergite 18 (Hd, CaFeSi2O6) melt and a melt mixture of 50 mol% hedenbergite + 50 mol% diopside (Hd50Di50) 19 at high pressure and temperature conditions up to 6 GPa and 2329 K, using high-pressure ultrasonic 20 technique combined with synchrotron radiation in a multi-anvil apparatus. Our results show that 21 Fe can significantly reduce the sound velocity while increase the density of silicate melts. 22 Comparing the Di, Hd, and Hd50Di50 melts, we find that the sound velocity does not mix linearly 23 2 for melts in the Hd-Di join, whereas the density for Hd-Di melts at high pressures can be well-24 described by linear mixing. Combined with melt geometry and melt compaction models, we 25 applied our results to study the stability and seismic signature of Fe-rich silicate melts in the Earth’s 26 upper mantle. For the low-velocity zone (LVZ) in mantle asthenosphere, although the degree of 27 seismic velocity reduction can be explained by the presence of a small amount of partial melt 28 distributed in film/band geometry along grain boundaries, Fe-rich melts formed at this depth are 29 unlikely to be gravitationally stable, but may be dynamically unextractable if the melt supply is 30 continuous. For the low-velocity layer (LVL) above the mantle transition zone, the presence of Fe-31 rich melts (with FeO>~10 wt%) distributed in textural equilibrium with the ambient mantle is a 32 plausible explanation.

diopside↗
DOE OSTI2021

The gallium isotopic composition of the Moon

Here, in this study, we present new Ga isotope data from a suite of 28 mare basalts and lunar highland rocks. The δ 71 Ga values of these samples range from -0.10 to +0.66‰ (where δ 71 Ga is the relative difference between the 71 Ga/ 69 Ga ratio of a sample and the Ga-IPGP standard), which is an order of magnitude more heterogeneous than δ 71 Ga values in terrestrial magmatic rocks. The cause of this isotopic heterogeneity must be established to estimate the bulk δ 71 Ga value of the Moon. In general, low-Ti basalts and ferroan anorthosite suite (FAS) rocks have δ 71 Ga values that are lower than high-Ti basalts and KREEP-rich rocks. The observation that rocks derived from later forming LMO cumulates have higher δ 71 Ga values suggests that Ga isotopes are fractionated by processes that operate within the chemically evolving LMO, rather than localized degassing or volatile redistribution. Correlations between indices of plagioclase removal from the LMO (e.g. Eu/Eu*) with Ga isotope ratios suggest that a Δ 71 Ga plagioclase-melt of -0.3‰, (where Δ 71 Ga plagioclase-melt is the isotopic fractionation associated with crystallization of plagioclase from a melt), could drive the observed isotopic fractionation in high-Ti mare basalts and KREEP-rich rocks. This would be consistent with the observation that FAS rocks have δ 71 Ga values that are lower than mare basalts. However, the addition of KREEP-like material into the mare basalt source regions would not contribute enough Ga to perturb the isotopic composition outside of analytical uncertainty. Thus, basalts derived from early formed LMO cumulates such as those from Apollo 15, would preserve light Ga isotopic compositions despite containing modest amounts of urKREEP. We estimate that the δ 71 Ga value of the LMO was ~0.14‰ prior to the onset of plagioclase crystallization and extraction. Whether this δ 71 Ga value is representative of the initial BSM cannot be ascertained from the current dataset. It remains plausible that the Moon accreted with a heavier Ga isotopic composition than the Earth. Alternatively, the Moon and Earth could have accreted with similar isotopic compositions (BSE = 0.00 ± 0.06‰, Kato et al., 2017) and volatile loss drove the LMO to higher δ 71 Ga values prior to formation of the lunar crust.

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

Quartz-hosted inclusions and embayments reveal storage, fluxing, and ascent of the Mesa Falls Tuff, Yellowstone

Quartz-hosted glasses from the Mesa Falls Tuff provide a geochemical window into the pre-eruptive magmatic system from one of Yellowstone's largest-volume caldera-forming eruptions. H 2 O and CO 2 concentrations, along with major and trace elements, were measured in both fully enclosed glass inclusions and partially enclosed embayments in the same quartz crystals. Major elements are largely consistent between the inclusions and embayments, except for K 2 O and Na 2 O. Of note, K 2 O is enriched by ~1 wt.% in embayment interiors relative to inclusions. Most trace elements are also enriched in the embayment interiors compared to inclusions from the same crystals. Fractionation trends of trace elements are consistent with ~30-60% crystallization. Quartz-hosted glass inclusions preserve 3.1±0.9 wt.% H 2 O and 493±227 ppm CO 2 whereas embayment interiors have 0.9±0.1 wt.% H 2 O and 399±229 ppm CO 2 . The CO 2 is roughly similar, but the distinct ~2 wt.% discrepancy between inclusion and embayment interior H 2 O contents may have been produced by CO 2 fluxing sourced from underplated Yellowstone basalts. H 2 O gradients within embayments are flat in their interiors and modified by sharp positive gradients near embayment exteriors which were produced by post-eruptive rehydration. CO 2 gradients occur as gently sloping concentration gradients that extend inward 150 to 250 μm from the embayment exterior. Finally, finite-difference 1D and 2D diffusion modeling indicates the distribution of H 2 O and CO 2 in embayments was produced by slow, fluid-saturated decompression that preceded rapid ascent during the caldera-forming eruption.

2D diffusion modeling↗
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 OSTI2025

Probing Cretaceous-Paleogene crustal thickness in southern Tibet using quartz-zircon chronobarometry

Knowledge of the crustal thickness history of southern Tibet during the India-Asia collision is key to understanding what deformation mechanisms accommodated northward propagating crustal shortening. Thermoisotopic models suggest that a relatively thin margin of southern Asia persisted from ca. 200 to 45 Ma whereas trace element paleodepth proxies are interpreted to indicate a thinning phase from ca. 100 to 65 Ma when the crust reached ∼30 km-thickness. As neither of these methods directly measures crustal thickness, resolution of this conflict awaits development of a method that can. Here, in this study, we place bounds on the thickening history of southern Tibet using coupled quartz-zircon thermobarometry and U-Pb geochronology of granitoid plutons in the Gangdese batholith. We find that crustal thicknesses were as high (or higher) as 60–75 km at 65 Ma, or 10–15 Ma prior to the onset of hard continental collision. Magmatic inflation in the lower crust likely contributed to crustal thickening and thermal weakening of the lower crust, suggesting that pure shear in the lower crust was likely the primary accommodation mechanism for N-S shortening during collision. These data are in partial agreement with results of thermoisotopic models of crustal thickness but contradict thickness histories derived from empirical trace element proxies.

Chronobarometry↗
DOE OSTI2021

Chemical characteristics of iron meteorite parent bodies

The projected relative abundances of the highly siderophile elements (HSE; Re, Os, Ir, Ru, Pt, and Pd) for bulk parent bodies of 10 magmatic iron meteorite groups/grouplet (IC, IIAB, IIC, IID, IIF, IIIAB, IIIF, IVA, IVB, and South Byron Trio) are broadly similar and show no resolvable differences between noncarbonaceous (NC) and carbonaceous (CC) genetic heritage. The processes driving genetic isotopic heterogeneity in the early Solar System, therefore, evidently did not leave discernable chemical fingerprints with respect to HSE relative abundances on the bulk planetesimal scale. By contrast, the absolute abundances of HSE projected for parent body cores, which reflect core size, are more variable and, on average, higher in CC bodies compared to NC bodies. Overall, bulk core chemical compositions, as well as core size, are linked to the distribution of Fe within a parent body, which is controlled by its oxidation state. The CC parent bodies are constrained to have formed under heterogeneous oxidizing conditions which were, on average, more oxidizing than those of the NC environment.

Hilton, Connor D.↗
DOE OSTI2022

Nickel-rich, volatile depleted iron meteorites: Relationships and formation processes

Ungrouped iron meteorites Tishomingo, Willow Grove, and Chinga, and group IVB iron meteorites, are Ni-rich. Similarities include enrichments of 10–100 × CI for some refractory siderophile elements, and equivalent depletions in more volatile siderophile elements. Superimposed on the overall enrichment/depletion trend, certain siderophile elements (P, W, Fe, Mo) are depleted relative to elements of similar volatility. All three ungrouped irons derive from parent bodies formed in the early Solar System. Willow Grove and Chinga are characterized by cosmic ray exposure corrected 182 W/ 184 W consistent with metal-silicate segregation on their parent bodies within 1–3 Myr of Solar System formation, within the age range determined for segregation of magmatic iron meteorite parent bodies, including group IVB irons. Tishomingo is characterized by a younger model age 4–5 Myr subsequent to Solar System formation, reflecting either late stage melting resulting from 26 Al decay, or an impact resetting. The discovery of stishovite in Tishomingo, indicating exposure to a minimum shock pressure of 8–9 GPa, is consistent with the latter. Stishovite in Tishomingo and chromite included in troilite-daubréelite in IVB irons allows oxygen isotopic composition comparison between these meteorites. Different mass independent oxygen isotopic compositions of IVB irons and Tishomingo indicate genetically distinct parent bodies. By contrast, mass independent Mo isotopic compositions overlap within analytical uncertainties, indicating a similar, carbonaceous chondrite (CC) type genetic heritage. Molybdenum and 183 W isotopic data for Chinga and Willow Grove indicate derivation from CC type parent bodies. Willow Grove shares Mo and 183 W isotopic compositions with the Ni-rich South Byron Trio (SBT) grouplet and the Milton pallasite. These Ni-rich meteorites likely formed in the same general nebular environment as other CC planetesimals, likely the outer Solar System. Highly siderophile element (HSE) abundances of Willow Grove and Tishomingo are similar to some IVB meteorites, consistent with formation by moderate degrees of fractional crystallization from initial metallic melts with low S and P, and modestly fractionated HSE. The comparable HSE abundances of Tishomingo and Willow Grove to some IVB irons, yet substantially higher Ni concentrations, indicate formation on parent bodies with lower bulk HSE abundances or HSE concentration in proportionally smaller volumes of metal. HSE abundances in Chinga are considerably lower than in IVB irons, highly fractionated, and processes responsible for these remain elusive. For IVB irons and these ungrouped irons, high temperature condensation likely dominated the enrichment and depletion of the refractory and volatile siderophile elements, respectively. Parent body degassing may have also played a role. Relative depletion of volatile siderophile elements is not, however, a universal feature of high-Ni meteorites. The SBT and Milton pallasite are Ni-rich, but less depleted in the more volatile siderophile elements. Nickel enrichment was likely driven by oxidation of Fe metal during parent body accretion or core segregation. Oxidation of the Tishomingo and Willow Grove parent bodies may have occurred at ~ IW + 1, indicated by relative Mo and W depletions due to metal/water reaction during differentiation. Late-stage reduction, indicated by the presence of Cr-bearing sulfides in Tishomingo and IVB irons, may have resulted from exhaustion of the oxidant.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗
DOE OSTI2023

The effect of composition, temperature and pressure on the oxidation state and coordination environment of copper in silicate melts

Copper is a redox variable element that may occur as Cu 0 , Cu + , and Cu 2+ in the Earth’s crust. The oxidation state will affect its partitioning between coexisting minerals, melts and fluids and hence its behaviour in magmatic processes. Copper bearing silicate glasses were quenched from melts with 23 synthetic compositions (19 CaO-Na 2 O-MgO-Al 2 O 3 -SiO 2 (CNMAS), two “granites” containing K 2 O ± H 2 O, and Fe-bearing “MORB” and “andesite”) equilibrated at oxygen fugacities (fO 2 ), expressed in log units relative to the fayalite-magnetite-quartz (FMQ) buffer, ranging from −0.7 to 14, temperatures from 900 to 1500 °C and pressures from 0 to 2.5 GPa. Cu K-edge X-ray absorption near edge structure (XANES) spectra were recorded from the glasses and a pre-edge feature in the XANES spectra was found to scale with the proportion of Cu + . Cu + /ΣCu (where ΣCu = Cu + + Cu 2+ ) was quantified by fitting the intensity of the pre-edge feature as a function of fO 2 to the thermodynamically expected relationship. Cu + /ΣCu was found to only weakly depend on melt composition, with more basic melts (e.g., basalts rather than granites) preferentially stabilising Cu 2+ . Increasing temperature stabilises Cu + , while increasing pressure had little effect on Cu + /ΣCu in CNMAS melts but preferentially stabilised Cu 2+ in granite melts. Cu + /ΣCu can be predicted in silicate melts by the empirical equation: log(Cu 2+ /Cu + ) = 0.25(ΔFMQ + 8.58 − 25050/T + 940P/T – 0.02P) − 4.73 + 5400/T + 1.99Λ + (280P – 90P2) where T is temperature in K, P is pressure in GPa and is the optical basicity of the composition. The effects of fO 2 , melt composition, temperature and pressure on Cu + /ΣCu indicate that Cu + will be the dominant oxidation state in terrestrial silicate melts (e.g., Cu + /ΣCu = 99% in an andesitic melt at 900 °C, 1 GPa and ΔFMQ = 1). The electron exchange reaction Cu 2+ + Fe 2+ → Cu + + Fe 3+ occurs on cooling and, given the abundance of Fe in natural melts, the oxidation state of Cu in natural glasses is unlikely to correspond to that in the original melt.

Copper↗
DOE OSTI2024

Hydrothermal solution calorimetry in acidic aqueous solutions and revisiting the standard partial molal thermodynamic properties of Nd 3+ from 25 to 300 °C

The mobility of rare earth elements (REE) can be predicted in aqueous fluids using geochemical modeling but the accuracy of these models strongly depends on the availability of robust thermodynamic properties for the REE aqueous species. The REE 3+ aqua ions are important in the derivation of the formation constants of all the major REE complexes including the chloride, sulfate, and fluoride species which predominate in many hydrothermal-magmatic systems. However, the thermodynamic properties of the REE 3+ aqua ions are still commonly derived from the Helgeson-Kirkham-Flowers (HKF) equation of state parameters tabulated several decades ago. The standard state thermodynamic properties at reference conditions (25 °C and 1bar) and their extrapolations to high temperature need to be verified, if not revised, based on hydrothermal experiments. In this study, the enthalpy of solution was measured for synthetic Nd hydroxide from 25 to 150 ºC to retrieve the standard partial molal thermodynamic properties of Nd 3+ as a function of temperature. The experiments were conducted in aqueous perchloric acid based solutions with starting pH of 2 and varying ionic strength (0.01 to 0.09 mol/kg NaClO 4 ). The standard partial molal enthalpy of formation (Δ f H°) of Nd 3+ derived from the experimental study displays differences of up to 10 kJ/mol compared to the enthalpy values derived from the HKF equation of state in the studied temperature range. These inaccuracies are resolved by adjusting the standard partial molal Gibbs energy of formation (Δ f G°) of Nd 3+ at 25 °C and 1 bar from -672.0 to -679.7 kJ/mol. The heat capacity function (C p °) derived between 25 and 150 ºC can be described by: C p ° = a 0 + a 1 ·T + a 2 ·T -2 , with a 0 = 1256, a 1 = -2.68, a 2 = -55.56·10 6 and T in Kelvin. A set of recommended thermodynamic properties is provided for the Nd 3+ aqua ions and corrections are provided for the chloride and fluoride species to remain internally consistent with the experimentally derived properties. These results allow predicting accurately the solubility of monazite between 25 and 300 ºC. Before these corrections, the properties for the Nd 3+ aqua ions derived from the HKF parameters resulted in up to ~1.5 orders of magnitude lower monazite solubility than determined experimentally. Here, a revision of the REE +3 aqua ions properties is necessary to accurately predict the mobility of REE in hydrothermal acidic solutions.

58 GEOSCIENCES↗
DOE OSTI2025

NdPO 4 solubility and aqueous Neodymium speciation in supercritical fluids: An experimental study at 500–700 °C and 1.7 kbar

A key aspect in the formation of rare earth elements (REE) deposits is the role of REE transport as aqueous REE complexes in supercritical hydrothermal solutions, where the nature of the aqueous complex is controlled by solution composition, temperature and pressure. Despite chloride being considered as one of the most abundant transporting ligands in magmatic-hydrothermal fluids, experimental investigations on the stability of aqueous REE chloride complexes are scarce above 300 °C. In this study, synthetic NdPO 4 crystals were reacted with non-saline and saline (0, 0.05 and 0.5 mNaCl), acidic (0.01 mHCl) aqueous solutions in a series of solubility experiments conducted at 500–700 °C and 1.7 kbar, where the solubilities were determined using a stable Nd isotope ( 145 Nd isotope spike) dilution technique. NdPO 4 solubility ranges between 28 ppm and 10,858 ppm, where solubility increases with both temperature and salinity. At 500 °C, log mNdPO 4 increases from –3.93 to –1.60 and there is a strong correlation between NdPO 4 solubility and NaCl concentrations (slope of 1.2 ± 0.3), indicating stabilization of the Nd chloride aqueous complexes with a stoichiometry corresponding to NdCl 2+ . At 600 °C, this correlation is weaker (slope of 0.4, log mNdPO 4 increases from –2.63 to –1.88) indicating the stabilization of both Nd chloride and hydroxyl species controlling solubility. At 700 °C, NdPO 4 solubility is largely independent of NaCl concentration indicating that solubility is controlled by Nd hydroxyl complexes, where stoichiometry suggests the neutral Nd(OH) 3 0 species is dominant. The solubility product (Ksp) of NdPO4 is derived from experimental data with the relation: log K sp = -41.81 – 0.057T – 20987/T, with T temperature in Kelvin. Comparison of the measured Nd phosphate solubility to thermodynamic predictions using the available Helgeson-Kirkham-Flowers equation of state parameters for aqueous Nd complexes indicate that predictions are up to three orders of magnitude lower compared to experimental observations. This discrepancy is most pronounced in saline solutions, suggesting that thermodynamic properties of the REE chloride species in supercritical fluids require revision. Numerical simulations of fluid-rock interaction between acidic, saline fluids and a Strange Lake felsic mineral assemblage demonstrates that NdPO 4 solubility predictions from models are four to six orders of magnitude lower than those calculated based on empirical fits from experiments, which suggests that acidic, saline fluids may play an important role in mobilizing large amounts of light REE from 450 to 700 °C.

58 GEOSCIENCES↗
DOE OSTI2020

Seismic tomography of compressional wave velocity and attenuation structure for Makushin Volcano, Alaska

Located in the central portion of the Alaska-Aleutian subduction zone, Makushin Volcano is among the most active volcanoes in the United States and has been classified as a very high threat by the National Volcano Early Warning System based on its eruptive history and proximity to the City of Unalaska and international air routes. In 2015, five standalone seismic stations and three mini seismic arrays of 15 stations each were deployed on Unalaska Island to supplement the Alaska Volcano Observatory permanent seismic network. This temporary array was operational for one year. Taking advantage of the increased azimuthal coverage and the array's increased earthquake detection capability, we developed body-wave Vp and Vp/Vs seismic images of the velocity structure beneath the volcano. Our results show a complex structure with the upper 4 km of the crust. The shallow high-Vp features possibly delineate remnant magma pathways or conduits. Low-Vp regions are found east of the caldera at approximately 4–8 km depth. This is in agreement with previous seismic tomographic work and InSAR-based geodetic models, which had identified this region as a possible long-term source of magma. A three-dimensional compressional wave attenuation model complemented the velocity models. In general, we observe regions of corresponding low P-wave velocity and high attenuation east of the summit caldera between 5 and 7 km depth, supporting the presence of a melt-rich zone. The absolute Vp values (5.5–6.0 km/s) in this region are consistent with the presence of a magma body with a mafic signature. Here, we suggest that a triggering mechanism that involves mafic injections into a more evolved shallow magma chamber is representative of the current magmatic system below the volcano, similarly to what petrological studies have inferred for past eruptions.

58 GEOSCIENCES↗
DOE OSTI2022

Volcanology, Geochemistry, and Petrology Perspectives on Integrated, Coordinated, Open, Networked (ICON) Science

This article is composed of a commentary about the state of Integrated, Coordinated, Open, and Networked (ICON) principles in Volcanology, Geochemistry, and Petrology (VGP), and discussion on the opportunities and challenges of adopting them. VGP encompasses a broad field that addresses volcanic, magmatic, hydrothermal, geomicrobial systems; process investigations that span the physical, geochemical and biological realms, including planetary geology; and one that is extensively supported by state-of-the-art research facilities. We suggest that an open, inclusive, collaborative and evolving model of an international coordinated network is critical to answering the most pressing challenges in VGP. In this commentary piece, we begin to discuss the elements of, challenges to, and path forward in developing such a model. For this team, ICON means collaboration, equitable access to data for the entire scientific community, and forging of partnerships that potentially contribute to more innovative ways of coordinating and sharing research. It also means bringing more equity to science, by implementing effective measures which consider access to funding, analytical equipment, resources, and mentors. More importantly, ICON to us means having important conversations around what we value in the advancement of science, perhaps exploring outside the idea of meritocracy and evaluating what individual traits can contribute to science outside what has traditionally been considered the norm.

58 GEOSCIENCES↗
DOE OSTI2022

Strong Upper-Plate Heterogeneity at the Hikurangi Subduction Margin (North Island, New Zealand) Imaged by Adjoint Tomography

We use earthquake-based adjoint tomography to invert for three-dimensional structure of the North Island, New Zealand, and the adjacent Hikurangi subduction zone. The study area, having a shallow depth to the plate interface below the North Island, offers a rare opportunity for imaging material properties at an active subduction zone using land-based measurements. Starting from an initial model derived using ray tomography, we perform iterative model updates using spectral element and adjoint simulations to fit waveforms with periods ranging from 4–30 s. We perform 28 model updates using an L-BFGS optimization algorithm, improving data fit and introducing P- and S-wave velocity changes of up to ±30%. Resolution analysis using point spread functions show that our measurements are most sensitive to heterogeneities in the upper 30 km. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone. Lateral velocity structures in the upper 5 km correlate well with New Zealand geology. Here, the inversion reveals increased along-strike heterogeneity on the margin. In Cook Strait we observe a low-velocity zone interpreted as deep sedimentary basins. In the central North Island, low-velocity anomalies are linked to surface geology, and we relate velocity structures at depth to crustal magmatic activity below the Taupō Volcanic Zone.

58 GEOSCIENCES↗
DOE OSTI2023

Episodic Magma Hammers for the 15 January 2022 Cataclysmic Eruption of Hunga Tonga‐Hunga Ha'apai

Understanding the forces and magma system dynamics on timescales of seconds to minutes remains challenging. In the January 2022 phreatoplinian Hunga Tonga-Hunga Ha'apai eruption, four remarkably similar seismic subevents within a 5-min interval occurred during the intensifying early eruptive phase. The subevents are similar in waveforms and durations (~25 s each). Each subevent begins with an unusual negative P-wave polarity which is inferred, using full-wave seismic modeling, to be caused by an upward single-force mechanism at the volcano created by a magma hammer likely in response to magma flow blockage/constriction during the early part of the eruption as discharge rapidly increased over orders of magnitude with concomitant conduit geometry evolution and instability. Our proposed episodic magma hammer model is consistent with thermodynamic and phase properties of the magmatic mixture, and yields an estimate of conduit mass flow in agreement with vent discharge rates derived from satellite imagery of plume heights.

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↗