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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

De-risking superhot EGS development through 3D play fairway analysis: Methodology development and application at Newberry Volcano, Oregon, USA

The DErisking Exploration for multiple geothermal Plays in magmatic ENvironments (DEEPEN) project aims to de-risk exploration of geothermal plays in magmatic systems, with a focus on superhot and supercritical geothermal systems. This article presents advancements in geoscientific assessment of superhot enhanced geothermal system (EGS) resources as part of the DEEPEN project, particularly in exploration and identification of these systems. While there exists a general consensus on key components of conventional hydrothermal systems, there has not been a prior consensus on such components for EGS, let alone superhot EGS. The DEEPEN project identified these key components for superhot EGS plays, and used them to modify the traditional geothermal play fairway analysis (PFA) approach to identify favorable areas in superhot EGS plays in 3D. This was done through modeling both favorability and uncertainty separately for each key component of a superhot EGS play, incorporating diverse 3D geoscientific datasets, including geologic features, models produced from direct observations, single inversions, and joint inversions. The PFA is applied to Newberry Volcano as a form of validation. Finally, the results appear to validate the methodology, aligning with expectations derived from conceptual modeling, and highlighting the area targeted for EGS stimulation well NWG 55-29 as favorable, and suggesting potential additional areas worthy of exploration at Newberry.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

DEEPEN 3D PFA Index Models for Exploration Datasets at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the development of the DEEPEN 3D play fairway analysis (PFA) methodology for magmatic plays (conventional hydrothermal, superhot EGS, and supercritical), index models needed to be developed to map values in geoscientific exploration datasets to favorability index values. This GDR submission includes those index models. Index models were created by binning values in exploration datasets into chunks based on their favorability, and then applying a number between 0 and 5 to each chunk, where 0 represents very unfavorable data values and 5 represents very favorable data values. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were created using histograms of the distributions of each exploration dataset in combination with literature and input from experts about what combinations of geophysical, geological, and geochemical signatures are considered favorable at Newberry. This is in attempt to create similar sized bins based on the current understanding of how different anomalies map to favorable areas for the different types of geothermal plays (i.e., conventional hydrothermal, superhot EGS, and supercritical). For example, an area of partial melt would likely appear as an area of low density, high conductivity, low vp, and high vp/vs. This means that these target anomalies would be given high (4 or 5) index values for the purpose of imaging the heat source. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were produced for the following datasets: - Geologic model - Alteration model - vp/vs - vp - vs - Temperature model - Seismicity (density*magnitude) - Density - Resistivity - Fault distance - Earthquake cutoff depth model

15 GEOTHERMAL ENERGY↗
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 OSTI2020

Composition of Magma and Characteristics of the Hydrothermal System of Newberry Volcano, Oregon, From Magnetotellurics

Abstract We use 3‐D magnetotellurics to improve our understanding of the structure and magma composition of Newberry Volcano in Oregon, USA. Newberry is a broad shield volcano with a summit caldera and strongly bimodal magmatism. Newberry has long been the subject of geothermal exploration research, but that work has focused on the volcano's west flank, leaving the caldera largely unstudied with geophysical methods until recently. Our modeling shows a relatively resistive magma reservoir of approximately 50 Ωm. Our work builds upon recent seismic models and petrological analysis to interpret Newberry's magma reservoir as a dry rhyolite with 8–11% partial melt, which matches the seismically determined melt fraction. Finding the conditions within the magma reservoir that allow the resistivity and seismic analysis to come to the same melt fraction helps us narrow down the magma temperature and composition. From this we infer a dry rhyolitic magma at 850 °C. We also image a prominent vertical conductive anomaly along the south rim below the vent that produced the most recent eruption. The anomaly extends from magma reservoir depths of 3 km to 1 km below the surface where it disperses into the caldera fill. We interpret this as the main conduit for magmatic fluids to reach Newberry's hydrothermal system. Other features that our model shows include higher conductivity along the caldera rim and a resistive older pluton on the west flank.

Bowles‐Martinez, E.↗
DOE OSTI2021

Volatile Element Evidence of Local MORB Mantle Heterogeneity Beneath the Southwest Indian Ridge, 48°–51°E

Abstract The mantle source beneath the Southwest Indian Ridge (SWIR) reflects a complex history of contamination. Magmatic volatile contents are vital tracers of these kinds of heterogeneities, which may fractionate otherwise constant volatile/non‐volatile elemental ratios, such as the H 2 O/Ce ratio. Although several studies have recently used trace element and isotopic data to address mantle source heterogeneity and magmatic processes at SWIR 48°–51°E region, volatile element constraints provide a valuable test of models for the origins of mantle heterogeneities in this region. Here, we present new data for nine rare basaltic glass samples from the 48°–51°E region, which enable careful assessment of the effects of primary versus secondary processes on the glass volatile contents. These samples are strongly affected by variable extents of carbon degassing, and shallow assimilation of Cl‐rich fluid, but also reveal consistently high H 2 O/Ce ratios (458.8 ± 14.9), among the highest in MORBs, that cannot be explained by late‐stage secondary processes, crustal assimilation, or simple melting of peridotite mantle at variable depths. Instead, the high H 2 O/Ce ratios are features of the mantle source composition. The 48°–51°E region is notably more depleted in highly incompatible trace elements relative to other regions of the SWIR, although this depletion is not apparent in H 2 O, which is similarly abundant throughout the SWIR. We link the high H 2 O/Ce ratios in these glasses with other trace element characteristics diagnostic of subduction and fluid addition, suggesting that the mantle source reflects signatures of a refractory mantle residue that previously melted within a subduction zone.

Wang, Wei↗
DOE OSTI2025

Eocene-Oligocene Metamorphism, Fluid Flow and Deformation in the Ruby Mountains-East Humboldt Range Metamorphic Core Complex

Metamorphic core complexes are ubiquitous in collapsed orogens globally and play a primary role in crustal exhumation. Here, we investigate the metamorphic history of the Ruby Mountains-East Humboldt Range metamorphic core complex, Nevada, using petrochronology to understand how magmatism, metamorphism and deformation interact to modulate crustal rheology, and the timing of exhumation within the Sevier orogenic belt. Field study, microstructural analysis, thermobarometry and laser ablation split-stream monazite, titanite and allanite petrochronology of the mylonitic footwall were integrated to elucidate the pressure-temperature-time (P-T-t) evolution. Major-, trace-element and quartz-in-garnet thermobarometry show peak metamorphism occurred at 5.5–6 kbar and 600–650°C across the study area. Monazite and titanite U-(Th)-Pb petrochronology constrain this episode of metamorphism to 88–81 Ma with the structurally deepest samples yielding dates down to ca. 71 Ma. The preservation of peak metamorphic mineral assemblages indicates these rocks remained at near-peak temperature conditions into the Cenozoic. Cenozoic metamorphism occurred during a punctuated episode contemporaneous with the Great Basin ignimbrite flare up. All samples show a distinct 39–37 Ma population of monazite, titanite and allanite dates that correlate with the emplacement of quartz diorite and gabbro intrusions throughout the footwall. Accessory phases contain zoning textures in backscattered electron images that are consistent with fluid-mediated dissolution-reprecipitation reactions. The U-(Th)-Pb analyses from these zones yield 39–27 Ma dates interpreted as the timing of fluid flow. The results of this study elucidate the integrated importance of magmatism, melting and fluid flow in driving the coupled evolution of Late Cretaceous metamorphism within the Sevier hinterland and the subsequent Cenozoic metamorphic core complex development that exhumed the middle crust.

58 GEOSCIENCES↗
DOE OSTI2020

The internal structure and geodynamics of Mars inferred from a 4.2-Gyr zircon record

Combining U–Pb ages with Lu–Hf data in zircon provides insights into the magmatic history of rocky planets. The Northwest Africa (NWA) 7034/7533 meteorites are samples of the southern highlands of Mars containing zircon with ages as old as 4476.3 ± 0.9 Ma, interpreted to reflect reworking of the primordial Martian crust by impacts. We extracted a statistically significant zircon population (n= 57) from NWA 7533 that defines a temporal record spanning 4.2 Gyr. Ancient zircons record ages from 4485.5 ± 2.2 Ma to 4331.0 ± 1.4 Ma, defining a bimodal distribution with groupings at 4474 ± 10 Ma and 4442 ± 17 Ma. We interpret these to represent intense bombardment episodes at the planet’s surface, possibly triggered by the early migration of gas giant planets. The unradiogenic initial Hf-isotope composition of these zircons establishes that Mars’s igneous activity prior to ~4.3 Ga was limited to impact-related reworking of a chemically enriched, primordial crust. A group of younger detrital zircons record ages from 1548.0 ± 8.8 Ma to 299.5 ± 0.6 Ma. The only plausible sources for these grains are the temporally associated Elysium and Tharsis volcanic provinces that are the expressions of deep-seated mantle plumes. The chondritic-like Hf-isotope compositions of these zircons require the existence of a primitive and convecting mantle reservoir, indicating that Mars has been in a stagnant-lid tectonic regime for most of its history. Our results imply that zircon is ubiquitous on the Martian surface, providing a faithful record of the planet’s magmatic history.

58 GEOSCIENCES↗
DOE OSTI2025

A tale of two planets: Disparate evolutionary models for Mars inferred from radiogenic isotope compositions of Martian meteorites

The radiogenic isotopic compositions of basaltic Martian meteorites (shergottites) and clinopyroxene/olivine cumulate meteorites (nakhlite/chassignites) are used to define the global evolution of Mars. However, the two main groups of meteorites demonstrate that their sources underwent divergent styles of magmatic evolution. The shergottites portray a planet that differentiated ~4.52 billion years ago via solidification of a magma ocean, producing incompatible element-depleted and -enriched reservoirs that remained isolated until melt production. In contrast, the reservoir from which the nakhlite/chassignites derive may have formed earlier, produced melts that fractionated Sm/Nd and Hf/W differently, was compositionally less variable, and experienced a significantly more complex history following primordial differentiation than the shergottite sources. The disparate histories recorded by these two groups of meteorites elucidate important questions that could be addressed by acquiring additional samples. Obtaining samples that shared the isotopic systematics of the shergottites would provide confidence that extrapolating the primordial differentiation history of Mars from shergottite radiogenic isotope systematics is reasonable. Returned samples from Mars will also constrain the physical locations of the meteorite source regions, providing insights into the general structure of the Martian mantle. In addition, they will help constrain the phases present in the martian mantle during melting and the conditions under which they are stable. Finally, identifying an evolved lithology that satisfies the geochemical and isotopic constraints placed on the incompatible element-enriched endmember observed in the shergottites would define the nature of magmatic evolution on Mars and whether it is more akin to processes on the Earth or the Moon.

Mars↗
DOE OSTI2023

Validation of Play Fairway Analysis of the geothermal potential of Camas Prairie, south-central Idaho, by an exploration well

Abstract Play Fairway Analysis (PFA) methodology was adapted for geothermal exploration at Camas Prairie, Idaho. Geophysical data, structural and geologic mapping, volcanic rock ages and vent locations, and the distribution of thermal springs and wells all indicated a relatively high geothermal potential along the southern margin of the Prairie. An exploration well (USU Camas-1) was drilled to a depth of 618.3 m to validate the PFA. A permeable zone was encountered at ~ 357.5 m with a maximum measured temperature of ~ 80 °C, which was suppressed following the injection of cold water. A moderate transmissivity of ~ 0.25–1 cm 2 /s estimated from an injection test as well a seasonal artesian flow at ~ 0.7 L/s corroborate the presence of a permeable zone. The existence of a lacustrine clay seal was confirmed near the bottom of the basin-fill sediment occupying the upper 314 m of the well. Geothermometers suggest the USU Camas-1 well water equilibrated at a reservoir temperature of ~ 120 °C. Based on the locations of both thermal and cold wells, geothermal fluids appear to be flowing upward along one or both of two fault systems. The presence of young basalts and elevated helium isotope ratios suggest that the heat source of Camas Prairie is magmatic. However, the faults may be acting as a conduit for geothermal fluids to rise from great depth without a shallow magmatic source being present. Camas Prairie is a promising area for geothermal development, but the relatively low reservoir temperatures indicate this resource may not be suitable for electric generation. Perhaps the best use would be for heating.

Lachmar, Thomas E. (ORCID:0000000249640119)↗
DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗
DOE OSTI2023

Andrieslombaardite, RhSbS, a new platinum-group mineral from the platiniferous Onverwacht Pipe, Republic of South Africa

A hundred years after the discovery of the Merensky Reef in 1924, it is appropriate to present the new mineral andrieslombaardite in honour of Andries Frederik Lombaard who was instrumental in its discovery. Andrieslombaardite, RhSbS, was first described as an unknown mineral from placer deposits associated with the Tulameen Alaskan-Uralian type complex, British Colombia, Canada (Raicevic and Cabri, 1976) but has since been reported from several other deposits including the platiniferous Driekop, Mooihoek, and Onverwacht pipes in the eastern Bushveld Complex, South Africa. The mineral and the name were approved by the Commission on New Minerals Nomenclature and Classification (CNMNC) of the International Mineralogical Association (IMA no. 2022-076) based on data in the co-type samples from Onverwacht and a co-type sample from the Yubdo stream, Birbir River, Ethiopia. Andrieslombaardite in the Onverwacht sample is a single 8 x 20 μm grain attached to laurite in a matrix of altered silicate and Fe-oxyhydroxide minerals. In the Yubdo samples, there are many grains of pale brownish gray andrieslombaardite up to 25 x 55 μm in size, included in Pt-Fe alloys, some associated with erlichmanite, and others attached to bornite and chalcopyrite. The reflectance values (R%) measured in air and in oil at the COM wavelengths are 48.3 and 33.0 (470 nm), 49.3 and 34.0 (546 nm), 51.0 and 35.9 (589 nm), and 51.8 and 36.7 (650 nm). The colour values x, y, Y, λd, and Pe in air are 0.317, 0.322, 50.3, 580, and 3.2, and in oil are 0.319, 0.324, 35.6, 579, and 4.5. The composition of andrieslombaardite is ideally RhSbS, but it contains variable amounts of Fe, Pt, Pd, and Ir that may substitute for Rh. The mineral is cubic with unit-cell dimensions of a = 6.0278(4) Å, V = 219.01(6) Å 3 and Z = 4. It was synthesised at 400 and 550°C using stoichiometric elemental amounts. It is a member of the cobaltite group. The mineralisation of the intrusive dunite pipes was probably introduced at high temperatures, under magmatic conditions. The primary assemblages were to a certain degree overprinted and redistributed by low-temperature hydrothermal fluids. The Pt-Fe alloys from Yubdo containing PGM inclusions such as andrieslombaardite in the Yubdo-Alaskan-type complex were formed at some post-magmatic stage owing to PGE remobilisation during hydrothermal or metamorphic episodes.

58 GEOSCIENCES↗
DOE OSTI2024

Utah FORGE: Geochemical Data for Cold Groundwaters and Produced Geothermal Fluids

Geochemical data for cold groundwaters and produced geothermal fluids around the Utah FORGE site. The data is compiled into four tables in the attached Excel File. Table 1 is a compilation of compositions (anions, cations, weak acids, oxygen, hydrogen, and carbon isotopes) for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 2 is a compilation of noble gas (He, Ne, Ar) and He and Ne isotopic compositions for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 3 provides values for calculated advective and diffusive fluxes of helium. Table 4 provides values of calculated subsurface stored heat between the Opal Mound fault and the Utah FORGE site, which are related to volumes of recently solidified magmatic heat sources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Transport of coexisting Ni-Cu sulfide liquid and silicate melt in partially molten peridotite

Transport of coexisting sulfide and silicate melts in partially molten peridotite contributes to the redistribution of chalcophile elements within the upper mantle as well as the genesis of magmatic Ni-Cu sulfide deposits, but has not been investigated systematically. Using laboratory experiments, theoretical calculations, and X-ray synchrotron microtomography, this study documents the topology and considers controls on the extraction of two immiscible liquids during partial melting of mantle peridotite. Under hydrostatic conditions, the measured dihedral angle at silicate melt-mineral-mineral junctions is 13.7–21.3°. Silicate melt is distributed along grain edges forming incompletely interconnected melt channels that disconnected by some dead ends at melt fractions ~7–9 vol%. Application of theoretically predicted permeability (k ~10 -14 -10 -16 m 2 ) permits estimation of the extraction velocity of silicate melt of 0.7-11.1 μm/day within a single interconnected melt channel. In the absence of silicate melt, isolated sulfide droplets (3.77 vol%) show a sulfide-olivine-olivine dihedral angle of 91.5–101.3°. However, in the presence of silicate melt, sulfide droplets (average size ~2.53±2.14μm, 1σ) are partially surrounded by silicate melt and stranded in triple junctions or melt pockets due to the limitation of the smallest dimension (0.3μm) of melt channels. Thus, the extraction of sulfide liquid is highly restricted by these dead ends and the smallest dimension of melt channels during porous flow of silicate melt. In contrast, during large-strain shear deformation (strain ~1.6–2.5), initially stranded sulfide droplets were elongated and extracted with silicate melt into liquid-rich sheets with a length of several hundred microns, constantly oriented at 14.3±4.5°to the shear plane and antithetic to the shear direction. The angle is lower than that (18-30°) of those sheets containing sulfide liquid only, indicating that silicate melt dominates liquid-rich sheets. Driven by stress, silicate melt-dominated liquid-rich sheets open the appropriately oriented grain boundaries between silicate minerals, thus providing an efficient pathway for the extraction of sulfide liquid during deformation. When such sheet-like channels remain open, sulfide droplets (> millimeter-scale) can be potentially mobile through high strain domains of the upper mantle, contributing to the addition of chalcophile elements and the fertilization of the lithospheric mantle.

high P-T experiment↗
DOE OSTI2023

DEEPEN Leapfrog Geodata Model Cleaned and Reformatted Exploration Datasets from Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the DEEPEN 3D play fairway analysis (PFA) conducted at Newberry Volcano for multiple play types (conventional hydrothermal, superhot EGS, and supercritical), existing geoscientific exploration datasets needed to be acquired, cleaned, reformatted, and assembled in Leapfrog Geothermal. This GDR submission includes all of the cleaned and reformatted (X (m), Y (m), elevation (m), processed data values) datasets used to build the Leapfrog Geodata model. Existing datasets were acquired from the GDR, from AltaRock, and from other sources. This yielded the following datasets: - Digital elevation model produced from LiDAR data by Ramsey and Bard, 2016 - MT surveys from 2006, 2011, 2014, and 2017 (including single inversions) - Gravity surveys from 2006, 2007, and 2011 (including single) - Earthquake catalogs from PNSN, LLNL, and the Newberry EGS Demonstration project - Seismic velocity model from Templeton et al., 2014 - The Frone, 2015 temperature model and a new one produced through extrapolating downhole temperature measurements and the SMU temperature at depth maps. Two versions of the new model are provided: 250 m spacing and 500 m spacing - EarthVision geologic model with alteration from Moser et al., 2016 - Well data from EGS well 55-29, deep geothermal wells, coreholes (GEO N-2 through 5) and several thermal gradient holes - "Newberry Well Data:" Location, simple lithology, directional survey data, and temperature data for the 34 wells and coreholes used in the Newberry PFA Although there are additional 2D datasets available in the area, such as aeromagnetic surveys, these were not included in the analysis. While it may be possible to project these datasets into three dimensions by assuming the surface measurements do not vary with depth, this method is associated with high uncertainty. Preexisting inversions of these data were unavailable, and inverting additional geophysical datasets is outside the scope of this project.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

DEEPEN: Newberry Volcano MT and Gravity Data - 2022 and 2023 Acquisition and Processing

As part of DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments), a 3D play fairway analysis (PFA) was conducted at Newberry Volcano in Central Oregon for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). For use in this PFA, combined full tensor broadband magnetotelluric (MT) and gravity data were acquired, processed, and inverted by Enthalpion Energy LLC (Enthalpion) with support from NREL staff. The data collection efforts took place from June 19th to July 24. Data were collected with the goal of gaining an improved understanding of the South Flank and the extent of the magma chamber. This GDR submission includes the raw data, single inversions, joint inversions, resolution matrices, and a report on these processes. More detailed information about the folder structure of the datasets included in this submission may be found on pages 94-97 of Magnetotelluric and Gravity Survey Report.pdf below.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

The solubility of titanite in silicate melt determined from growth and dissolution experiments

The solubility of titanite (CaTiSiO 5 ) in Si-rich melts was measured experimentally through growth experiments at 800–1000 °C, 0.5–1.0 GPa, log fO 2 ~ CCO–0.8, t = 72–168 h, and H 2 O = 0 to 4 wt.%, and in dissolution experiments at 925–1300 °C, 0.8 GPa, t = 18–118 h, and H 2 O = 1–10 wt.% in a piston-cylinder apparatus. Run product glasses in growth experiments were homogeneous, and iron loss suppressed ilmenite crystallization. Saturation concentrations in dissolution experiments were estimated by fitting measured diffusion profiles. Titanite solubility increases with increasing temperature and melt composition parameter \(M = {\text{ molar }}\left( {{\text{Na }} + {\text{ K }} + {\text{ 2Ca}}} \right)/\left( {{\text{Al }} \times {\text{ Si}}} \right)\) . Multiple linear regression of glass composition data from growth and dissolution experiments ( n = 29) plus 39 experiments from the LEPR database (Hirschmann et al. 2008 ) yielded the titanite solubility equation (adj. r 2 = 0.95): \(\left( {TiO_{2} } \right)^{melt} \left( {wt.\% } \right) = 0.978 \times M + 0.0048 \times T\left( K \right){-}5.90\) . This model correctly predicted undersaturation in 95% of 2344 experiments from the LEPR database that did not have titanite or rutile. Application to natural rocks yields saturation temperatures that are similar to independent temperature estimates. Finally, this equation should be useful for constraining the temperatures of titanite-saturated melts, for determining whether titanite saturation in magmatic source regions is likely, and for determining when titanite can crystallize and begin to exert an influence on melt geochemistry.

58 GEOSCIENCES↗
DOE OSTI2025

Isolated vesicles in submarine pumice: insights from the 2019 Volcano F eruption, Kingdom of Tonga

No reliable method can distinguish between subaerial and submarine pumice, but recent work shows that highly vesicular (total porosity, φ > 65%) submarine pumice from the 2012 Havre and 2021 Fukutoku Oka-no-Ba eruptions contains abundant isolated porosity (connectivity, C < 0.6). This differs from textural measurements of subaerial pumice with similarly high vesicularities, where high connectivities (C > 0.8) are often measured. To investigate the implications of this under-studied texture for bubble nucleation, growth, and coalescence dynamics at submarine eruptions, we analyze the textural characteristics of rafted pumice clasts from the 2019 submarine eruption of Volcano F, located along the Tonga-Kermadec Arc. We examined clasts collected while floating near the vent (50–200 km) and distally on the shores of Fiji (900 km). Using helium pycnometry and X-ray tomography (XRT), we quantified porosity and connectivity in 45 pumice lapilli, classifying them by vesicle macrotexture and edge morphology. Microvesicular clasts, the majority of which have breadcrust and cauliflower textures, exhibit the lowest connectivities (C = 0.61 ± 0.2), while macrovesicular clasts have high connectivities with C ~ 1. We propose that isolated porosity abundance varies according to clast location in the eruption column post-fragmentation. Microvesicular clasts erupted on the edges of the column and experienced high levels of seawater ingestion and subsequent quenching of the entire clast. Macrovesicular clasts, in contrast, erupted in the center of the eruption column and were thus thermally insulated, allowing magmatic bubbles in the molten pumice to nucleate, grow, and coalesce for longer. Our results and interpretation imply that low connectivity, microvesicular clasts preserve pre and syn-eruptive bubble processes at the 2019 Volcano F eruption.

Permeability↗
DOE OSTI2025

Initial Conditions of Planet Formation: Time Constraints from Small Bodies and the Lifetime of Reservoirs in the Solar Protoplanetary Disk

This review explores the timescales of the initial phase of planet formation, from nebular dust (CAIs and chondrules) to planetesimal accretion and differentiation, using evidence from meteorite research. Aluminium-Mg systematics of CAIs indicate either an extended period of CAI formation (∼0.3 Ma) or an initial 26 Al heterogeneity, with evidence supporting a homogeneous 26 Al abundance in the protoplanetary disk. Thermal and aqueous alteration on the parent body can disturb the U-Pb and Al-Mg chronometers in chondrules. Focusing on relatively robust isochron data from plagioclase of pristine (types ≤3.05) chondrites indicates a shift in chondrule formation locations, moving from the inner to the outer disk over time. Ages of basaltic achondrites show that silicate differentiation on small bodies was well underway within the first few million years (Ma) of our solar system. Their age record, however, reveals inconsistencies between different chronometers, partly caused by (i) secondary disturbances due to thermal metamorphism, aqueous alteration, or impacts, (ii) the presence of xenolithic minerals, and (iii) potentially variable initial 26 Al abundances due to disturbances at the mineral scale. Nucleosynthetic isotope data indicate that parent bodies of iron and stony meteorites formed in two distinct regions within the protoplanetary disk: the inner, non-carbonaceous (NC) and the outer, carbonaceous (CC) region. Based on Hf-W chronometry it has been demonstrated that NC and CC parent bodies of magmatic iron meteorites segregated their cores within ∼1–3 Ma after CAI formation, implying that parent body accretion occurred within <1 Ma in both reservoirs. Combining accretion ages with nucleosynthetic data further reveals that, at first order, NC and CC reservoirs in the solar protoplanetary disk were established within 1 Ma and existed over several Ma with limited exchange between them. In the CR chondrite accretion region of the disk, planetary bodies formed over at least 3 Ma, while in most other regions, formation spanned at least 1 Ma, with minimal changes in nucleosynthetic isotope compositions. Aerodynamical size sorting of dust likely introduced or amplified some of these variations.

26Al distribution↗