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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Magmatic water content controls the pre-eruptive depth of arc magmas

We report that vanguard efforts in forecasting volcanic eruptions are turning to physics-based models, which require quantitative estimates of magma conditions during pre-eruptive storage. Below active arc volcanoes, observed magma storage depths vary widely (~0 to 20 kilometers) and are commonly assumed to represent levels of neutral buoyancy. Here we show that geophysically observed magma depths (6 ± 3 kilometers) are greater than depths of neutral buoyancy, ruling out this commonly assumed control. Observed depths are instead consistent with predicted depths of water degassing. Intrinsically wetter magmas degas water and crystallize deeper than dry magmas, resulting in viscosity increases that lead to deeper stalling of ascending magma. The water–depth relationship provides a critical constraint for forecasting models by connecting depth of eruption initiation to its volatile fuel.

58 GEOSCIENCES↗
DOE OSTI2026

Mobility of heavy rare earth elements in magmatic-hydrothermal systems: Experimental determination of DyPO 4 solubility in supercritical fluids

Dysprosium (Dy) is part of the heavy rare earth elements (REE) and a critical component in the production of high-performance super magnets used in wind turbines, electric vehicles, and other green technologies. The HREE phosphate, xenotime, is commonly associated with hydrothermal systems, where the mobility, fractionation, and deposition of REE are controlled by the stability of aqueous complexes and solubility of REE minerals which depend on temperature, pressure, and fluid composition (salinity, pH, and ligand availability). Here, we conducted a series of synthetic DyPO 4 solubility experiments in NaCl-HCl-bearing aqueous solutions (pH of 2 at 25 °C) from 500 to 700 °C and 1.5 kbar using externally heated pressure vessels. DyPO 4 displays a strong prograde solubility with temperature across all studied fluid salinities (0–0.5 m NaCl). The DyPO₄ solubility isotherms at 500 °C, 600 °C, and 700 °C show a complex relationship with salinity where DyPO₄ solubility is highest (39.2–1552 ppm at 500–700 °C) in 0 m NaCl solutions, decreases strongly (21–331 ppm at 500–700 °C) in the 0.06 m NaCl solutions and remains relatively constant in the 0.25 m (42.8–248 ppm) and 0.5 m NaCl (40.2–124 ppm) solutions. Thermodynamic modeling suggests that pH increases as a function of temperature and salinity, and controls solubility due to the stability of chloride complexes (i.e., DyCl +2 and DyCl 2 + ) at 500–600 °C, transitioning to hydroxyl complexes (i.e., Dy(OH) 3 0 ) dominating at 700 °C and relatively alkaline pH. Comparison with NdPO 4 solubility data indicates that the light REE in monazite is significantly more soluble than the heavy REE in xenotime in saline fluids at high temperature (500–700 °C), implying efficient light/heavy REE fractionation in supercritical fluids. The lower solubility of DyPO 4 relative to NdPO 4 suggests that at high temperatures, acidic (pH 2–4), low-salinity (∼5 wt% NaCl) fluids will favor the precipitation of xenotime over monazite as it will cool from 700 °C to below 500 °C. These observations align with paragenetic mineral assemblages in NYF pegmatites (e.g., Baveno, Italy), where xenotime commonly immobilizes heavy REE during the high-temperature hydrothermal stages.

58 GEOSCIENCES↗
DOE OSTI2021

Uplift and Seismicity Driven by Magmatic Inflation at Sierra Negra Volcano, Galápagos Islands

Abstract Although episodes of surface uplift and elevated seismicity precede many volcanic eruptions, their temporal evolution is often complex, and apparently in contradiction to simple trends predicted by mechanical deformation models. Here, we use continuous global positioning system and seismic data recorded at Sierra Negra volcano, Galápagos Islands, to show how the edifice responded to stress changes driven by magma accumulation in a shallow sill. The rate of uplift varied during the 13 years and 6.5 m of inflation before the 2018 eruption. The number of earthquakes per unit of uplift increased exponentially with total uplift as the differential stress increased. Accordingly, the temporal seismicity rate varied in time as a function of both the total uplift and the uplift rate. The Gutenberg‐Richter b ‐value decreased as a function of total uplift. In the final six months before the eruption, a sequence of large ( M > 4) earthquakes regulated the state of stress on the fault, each being followed by 2–3 days of postseismic quiescence, and retarding the increase in seismicity rate. These earthquakes did not affect the overall uplift rate. Subsidence of 8.5 m accompanied the 2‐month eruption. On resumption of uplift, the number of earthquakes per unit of uplift was very low, and the b ‐value high, reflecting the relaxed stress state of the fault system. These observations show that crustal deformation becomes increasingly brittle at higher stress states, and supports theoretical models based on elastic‐brittle mechanics. They suggest that joint interpretation of deformation and seismicity is key for forecasting future eruptions in similar volcanic settings.

Bell, Andrew Forbes↗
DOE OSTI2023

Assessing the Behavior of Moderately Volatile Elements on the Moon in Order to Constrain Processes of Magmatic Evolution and Planetary Accretion

Heterogeneous isotopic compositions of Zn and Ga in lunar basalts and highlands rocks indicate that secondary processing affected the stable isotopic composition of major lunar reservoirs. It follows that these fractionating processes must be identified and characterized in order to make robust estimates for the isotopic composition of the mare basalt source regions or the bulk silicate Moon (BSM). Contrary to previous assertions (e.g., Kato & Moynier, 2017), MVEs such as Zn and Ga do not behave similarly on the Moon, in line with differences in their chemical behavior and volatility. There is clear evidence that impact reworking fractionates zinc isotopes but, in addition, zinc isotopes are also highly sensitive to the formation of Ti-bearing oxides such as ilmenite. In contrast, gallium is far less susceptible to volatile reworking, but its isotopic composition is fractionated during incorporation into plagioclase (Wimpenny et al., 2022; Render et al., 2023). Based on these findings we conclude that current isotopic compositions of Zn and Ga in lunar rocks, such as those measured for the mare basalt suite, are unlikely to be representative of the bulk Moon when it accreted. As such, estimates for the MVE isotopic composition of the BSM derived from mare basalts may not be accurate, which may have implications for the incorporation of these data into dynamical models of the Moon’s formation. In the following, we discuss the main findings from this project in more detail, focusing on the behaviour of zinc and gallium isotopic systematics in lunar rocks and the processes that control their isotopic fractionation on the Moon.

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 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 OSTI2020

The Role of Premagmatic Rifting in Shaping a Volcanic Continental Margin: An Example From the Eastern North American Margin

Both magmatic and tectonic processes contribute to the formation of volcanic continental margins. Such margins are thought to undergo extension across a narrow zone of lithospheric thinning (~100 km). New observations based on existing and reprocessed data from the Eastern North American Margin contradict this hypothesis. With ~64,000 km of 2-D seismic data tied to 40 wells combined with published refraction, deep reflection, receiver function, and onshore drilling efforts, we quantified along-strike variations in the distribution of rift structures, magmatism, crustal thickness, and early post-rift sedimentation under the shelf of Baltimore Canyon Trough (BCT), Long Island Platform, and Georges Bank Basin (GBB). Results indicate that BCT is narrow (80–120 km) with a sharp basement hinge and few rift basins. The seaward dipping reflectors (SDR) there extend ~50 km seaward of the hinge line. In contrast, the GBB is wide (~200 km), has many syn-rift structures, and the SDR there extend ~200 km seaward of the hinge line. Early post-rift depocenters at the GBB coincide with thinner crust suggesting “uniform” thinning of the entire lithosphere. Models for the formation of volcanic margins do not explain the wide structure of the GBB. We argue that crustal thinning of the BCT was closely associated with late syn-rift magmatism, whereas the broad thinning of the GBB segment predated magmatism. Correlation of these variations to crustal terranes of different compositions suggests that the inherited rheology determined the premagmatic response of the lithosphere to extension.

58 GEOSCIENCES↗
DOE OSTI2021

Oxygen isotope ratios in zircon and garnet: A record of assimilation and fractional crystallization in the Dinkey Dome peraluminous granite, Sierra Nevada, California

Abstract The 119 Ma Dinkey Dome pluton in the central Sierra Nevada Batholith is a peraluminous granite and contains magmatic garnet and zircon that are complexly zoned with respect to oxygen isotope ratios. Intracrystalline SIMS analysis tests the relative importance of magmatic differentiation processes vs. partial melting of metasedimentary rocks. Whereas δ18O values of bulk zircon concentrates are uniform across the entire pluton (7.7‰ VSMOW), zircon crystals are zoned in δ18O by up to 1.8‰, and when compared to late garnet, show evidence of changing magma chemistry during multiple interactions of the magma with wall rock during crustal transit. The evolution from an early high-δ18O magma [δ18O(WR) = 9.8‰] toward lower values is shown by high-δ18O zircon cores (7.8‰) and lower δ 18O rims (6.8‰). Garnets from the northwest side of the pluton show a final increase in δ18O with rims reaching 8.1‰. In situ REE measurements show zircon is magmatic and grew before garnets. Additionally, δ18O in garnets from the western side of the pluton are consistently higher (avg = 7.3‰) relative to the west (avg = 5.9‰). These δ18O variations in zircon and garnet record different stages of assimilation and fractional crystallization whereby an initially high-δ18O magma partially melted low-δ18O wallrock and was subsequently contaminated near the current level of emplacement by higher δ18O melts. Collectively, the comparison of δ18O zoning in garnet and zircon shows how a peraluminous pluton can be constructed from multiple batches of variably contaminated melts, especially in early stages of arc magmatism where magmas encounter significant heterogeneity of wall-rock assemblages. Collectively, peraluminous magmas in the Sierran arc are limited to small <100 km2 plutons that are intimately associated with metasedimentary wall rocks and often surrounded by later and larger metaluminous tonalite and granodiorite plutons. The general associations suggest that early-stage arc magmas sample crustal heterogeneities in small melt batches, but that with progressive invigoration of the arc, such compositions are more effectively blended with mantle melts in source regions. Thus, peraluminous magmas provide important details of the nascent Sierran arc and pre-batholithic crustal structure.

Geochemistry & Geophysics↗
DOE OSTI2023

Chronology of alkali anorthosite 14304 clast “b” records basin forming impact at ~3.95 Ga

Alkali-suite rocks constitute one of three major suites of lunar crustal rocks. As such, constraining their formation timescales and petrogenesis is important for understanding the earliest magmatic history of the Moon. However, the magmatic history of alkali-suite rocks is partly obscured by superimposed effects of major basin-forming impact events on the lunar nearside. Consequently, unambiguous crystallization ages of samples from this suite of rocks have not been determined. Here, the aim of this study is to better understand the petrogenetic history of the alkali-suite and the potential superimposed effects of impact metamorphism by determining Sm-Nd, Rb-Sr, and 40 Ar/ 39 Ar ages for an alkali anorthosite clast from Apollo 14 lithic breccia 14304 termed clast “b”. The new chronologic measurements of clast “b” yield concordant Sm-Nd, Rb-Sr, and 40 Ar/ 39 Ar ages of 3947 ± 13 Ma, 3975 ± 34 Ma, and 3937 ± 37 Ma respectively, resulting in a weighted mean age of 3949 ± 11 Ma. This age is not interpreted to date an igneous event related to production of the lunar highlands crust and instead the chronology and petrography of clast “b” are most readily explained by an impact event at ~3.95 Ga that caused near-complete isotopic re-equilibration of the Sm-Nd, Rb-Sr, and 40 Ar/ 39 Ar chronometers. The weighted mean age of 3949 ± 11 Ma of clast “b” is several hundred Ma younger than 4.3–4.4 Ga ages typically determined for samples of other crustal rock suites but in very good agreement with independent estimates for the formation of Imbrium basin ejecta and other marginally older impact events which are thought to have been sampled at the Apollo 14 landing site. Thus, although petrologic and geochemical examination suggest that clast “b” is a pristine igneous clast, its age likely records an impact event at the Apollo 14 landing site. Whereas the various determined ages of clast “b” do not reflect the timescales of alkali-suite magmatism, the relatively low initial Sr and Nd isotopic compositions of clast “b” indicate that its protolith evolved with very low 147 Sm/ 144 Nd and 87 Rb/ 86 Sr that are distinct from estimates for urKREEP but similar to that of lunar plagioclase. This implies that the igneous protolith of clast “b” derived from a plagioclase-dominated KREEP-rich source that must have formed after the formation of the urKREEP source at ~4.35 Ga but well before the impact event recorded by clast “b” at ~3.95 Ga.

58 GEOSCIENCES↗
DOE OSTI2023

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

The timing of lunar solidification and mantle overturn recorded in ferroan anorthosite 62237

Ferroan anorthosite suite (FAS) rocks are widely interpreted to represent primordial lunar crust. Despite their importance in pinpointing the timing of lunar crust formation, robust chronological investigations for this rock type are scarce. Here, we report the Ar-Ar, Rb-Sr, and Sm-Nd isotopic systematics for the FAS troctolitic anorthosite 62237. The Ar-Ar isotopic system has been reset by a thermal event at 3710 ± 48 Ma, and the Rb-Sr isotopic systematics has been disturbed such that a Rb-Sr isochron age cannot be determined. However, an internal isochron for the Sm-Nd isotopic system has yielded an age of 4350 ± 73 Ma (MSWD = 2.0) with an initial ε 143 Nd CHUR of -0.53 ± 0.26. The mineral and whole-rock fractions of 62237 plot on the same internal isochron as FAS sample 60025. The combined datasets define an age of 4372 ± 35 Ma (MSWD = 4.0) with an initial ε 143 Nd CHUR of -0.17 ± 0.22. Literature Sm-Nd data for FAS and Mg-suite whole-rocks also plot on the 60025-62237 isochron. The coherence of data from both FAS and Mg-suite rocks examined thus far suggests that both rock suites formed contemporaneously from identical, or nearly identical, sources. In addition, the concordance of FAS and Mg-suite ages suggests that primordial crust solidification either involved both magmatic suites, or that Mg-suite magmatism was contemporaneous with FAS magmatism within resolution of the Sm-Nd chronometer. The ages for FAS and Mg-suite also coincide with the formation ages of the mare basalt source regions and urKREEP. Ferroan anorthosite suite rocks and urKREEP are thought to represent primordial LMO solidification products, whereas Mg-suite and the mare basalt source regions are argued to represent mixtures of various LMO crystallization products that were formed during density-driven overturn of the LMO. The concordance of ages implies that the 4372 ± 35 Ma Sm-Nd isochron records the age of mantle overturn, and that overturn occurred during, or shortly after, solidification of the LMO.

58 GEOSCIENCES↗
DOE OSTI2020

The formation and evolution of the Moon’s crust inferred from the Sm-Nd isotopic systematics of highlands rocks

Ages determined for magnesian and ferroan anorthosite crustal rock suites overlap, suggesting they formed contemporaneously about 4.3–4.5 Ga. A notable exception is the Sm-Nd age previously determined on Mg-suite gabbronorite 67667 which is at least 100 Ma younger than the youngest ferroan anorthosite. New chronologic measurements of 67667 presented here yield concordant Sm-Nd and Rb-Sr mineral isochron ages of 4349 ± 31 Ma and 4368 ± 67 Ma, suggesting the sample is older than previous estimates. Furthermore, a whole rock Sm-Nd isochron of Mg-suite rocks from the Apollo 14, 15, 16, and 17 landing sites yields an age of 4348 ± 25 Ma, indicating that Mg-suite magmatism was widespread and roughly contemporaneous on the lunar nearside. Here, analysis of Sm-Nd internal isochron ages confirms that Mg-suite magmatism was restricted to a period between about 4.33 and 4.35 Ga at the Apollo 14, 15, 16, and 17 landing sites and was synchronous with magmatism at the Apollo 16 site associated with the ferroan anorthosite suite between 4.35 and 4.37 Ga. Magnesian- and ferroan anorthosite suite rocks with ages younger than ~4.33 Ga appear to have experienced slow cooling in the deep lunar interior, so that the ages record when the samples cooled below the closure temperature of the Sm-Nd isotopic system and not the time they crystallized.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Complex subsurface hydrothermal fluid mixing at a submarine arc volcano supports distinct and highly diverse microbial communities

Hydrothermally active submarine volcanoes are mineral-rich biological oases contributing significantly to chemical fluxes in the deep sea, yet little is known about the microbial communities inhabiting these systems. Here we investigate the diversity of microbial life in hydrothermal deposits and their metagenomics-inferred physiology in light of the geological history and resulting hydrothermal fluid paths in the subsurface of Brothers submarine volcano north of New Zealand on the southern Kermadec arc. From metagenome-assembled genomes we identified over 90 putative bacterial and archaeal genomic families and nearly 300 previously unknown genera, many potentially endemic to this submarine volcanic environment. While magmatically influenced hydrothermal systems on the volcanic resurgent cones of Brothers volcano harbor communities of thermoacidophiles and diverse members of the superphylum “DPANN,” two distinct communities are associated with the caldera wall, likely shaped by two different types of hydrothermal circulation. The communities whose phylogenetic diversity primarily aligns with that of the cone sites and magmatically influenced hydrothermal systems elsewhere are characterized predominately by anaerobic metabolisms. These populations are probably maintained by fluids with greater magmatic inputs that have interacted with different (deeper) previously altered mineral assemblages. However, proximal (a few meters distant) communities with gene-inferred aerobic, microaerophilic, and anaerobic metabolisms are likely supported by shallower seawater-dominated circulation. Furthermore, mixing of fluids from these two distinct hydrothermal circulation systems may have an underlying imprint on the high microbial phylogenomic diversity. Collectively our results highlight the importance of considering geologic evolution and history of subsurface processes in studying microbial colonization and community dynamics in volcanic environments.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2025

Seismic imaging of the Ecuadorian forearc and arc from joint ambient noise, local, and teleseismic tomography: catching the Nazca slab in the act of flattening

SUMMARY The Ecuadorian Andes are a complex region characterized by accreted oceanic terranes driven by the ongoing subduction of the oceanic Nazca plate beneath South America. Present-day tectonics in Ecuador are linked to the downgoing plate geometry featuring the subduction of the aseismic, oceanic Carnegie Ridge, which is currently entering the trench. Using seismic tomography, we jointly invert arrival times of P and S waves from local and teleseismic earthquakes with surface wave dispersion curves to image the structure of the forearc and magmatic arc of the Ecuadorian Andes. Our data set includes > 100 000 traveltimes recorded at 294 stations across Ecuador. Our images show the basement of the central forearc is composed of accreted oceanic terranes with high elastic wave speeds. Inboard of the Carnegie Ridge, the westernmost forearc and coastal cordilleras display relatively low Vp and Vs and high Vp/Vs values, which we attribute to the increased hydration and fracturing of the overriding plate due to the subduction of the thick oceanic crust of the Carnegie Ridge. We additionally image across-arc differences in magmatic architecture. The frontal volcanic arc overlies accreted terranes and is characterized by low velocities and high Vp/Vs indicative of partial melt reservoirs which are limited to the upper crust. In contrast, the main arc displays regions of partial melt across a wider range of depths. The Subandean zone of Ecuador has two active volcanoes built on continental crust suggesting the arc is expanding eastwards. The mid to lower crust does not show indications of being modified from the magmatic process. We infer that the slab is in the process of flattening as a consequence of early-stage subduction of the buoyant Carnegie Ridge.

58 GEOSCIENCES↗