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

WUS324: Multiscale Full Waveform Inversion Approaching Convergence Improves Waveform Fits While Imaging Seismic Structure of the Western United States

Abstract We report a new model of radially anisotropic crustal and upper mantle structure of the western United States (WUS324) obtained from full waveform inversion of earthquake data. We ran three multiscale inversion stages beyond model WUS256 (Rodgers et al., 2022, https://doi.org/10.1029/2022jb024549 ) allowing them to approach convergence to fit a larger data set to a shorter minimum period of 16 s. WUS324 is based on 324 total iterations from its starting model, significantly more (16 times) than previous studies. Waveform misfit reductions are 66%–70% for both the inversion data and an independent validation data set providing confidence in the predictive power of the model. WUS324 provides much better fits and reveals shear wavespeed, v S , structure of this large region with more detail than previous waveform tomography models. We show representative images demonstrating the resolution of diverse seismic structure across this highly heterogeneous region including oceanic lithosphere, subducting slabs and continental magmatism.

58 GEOSCIENCES↗
DOE OSTI2025

Structure of the Ecuadorian Upper Plate From a Joint Seismic‐Gravity Inversion

The Ecuadorian portion of the South American subduction zone presents an interesting case study in the structure and complex evolution of an upper plate. There are outstanding questions about its tectonic history, composition, and magmatic processes. While previous studies have employed ambient noise tomography to image the Ecuadorian upper plate, surface wave inversions alone often lack sensitivity at relevant shallow depths. This limitation can be overcome with an independent, complementary data set, such as gravity. We have jointly inverted Rayleigh wave phase velocities and Bouguer gravity anomalies to provide a more detailed seismic velocity model of the Ecuadorian upper plate. Our joint inversion has yielded several key improvements from previous models. First, we observe much shallower slow velocities beneath major basins (the Manabí, Progreso, and Gulf of Guayaquil), better aligning with expected basin structure. Second, we identify a high-velocity block beneath the entire forearc, corresponding to the Piñon Terrane, with velocities suggesting the presence of ultramafic material. Third, we highlight a new narrow swath of slow velocities beneath the Ecuadorian Andes, which closely follows the active volcanoes along the Eastern Cordillera. The extent of these slow velocities coincides with the termination of active arc volcanism and the predicted location of the subducted Carnegie Ridge. The predicted compositions for the mid to lower crust in the region preclude a purely compositional explanation for these velocities, suggesting that some level of partial melt is necessary.

Birkey, Andrew [Univ. of Delaware, Newark, DE (Uni↗
DOE OSTI2021

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

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

58 GEOSCIENCES↗
DOE OSTI2021

Oxidized sulfur-rich arc magmas formed porphyry Cu deposits by 1.88 Ga

Most known porphyry Cu deposits formed in the Phanerozoic and are exclusively associated with moderately oxidized, sulfur-rich, hydrous arc-related magmas derived from partial melting of the asthenospheric mantle metasomatized by slab-derived fluids. Yet, whether similar metallogenic processes also operated in the Precambrian remains obscure. Here we address the issue by investigating the origin, fO 2 , and S contents of calc-alkaline plutonic rocks associated with the Haib porphyry Cu deposit in the Paleoproterozoic Richtersveld Magmatic Arc (southern Namibia), an interpreted mature island-arc setting. We show that the ca. 1886–1881 Ma ore-forming magmas, originated from a mantle-dominated source with minor crustal contributions, were relatively oxidized (1–2 log units above the fayalite-magnetite-quartz redox buffer) and sulfur-rich. These results indicate that moderately oxidized, sulfur-rich arc magma associated with porphyry Cu mineralization already existed in the late Paleoproterozoic, probably as a result of recycling of sulfate-rich seawater or sediments from the subducted oceanic lithosphere at that time.

58 GEOSCIENCES↗
DOE OSTI2025

U redox state tracked in mineralized hydrothermal carbonate with implications for U-Pb geochronology

U-Pb carbonate geochronology can directly constrain the timing and rates of important geological processes. However, the mechanisms and controls on U incorporation, distribution, and retention in carbonate minerals remain unclear, limiting geological interpretations. Here X-ray absorption spectroscopy (µXAS) and in-situ U-Pb carbonate geochronology are combined to temporally track U distribution and redox state in a porphyry-epithermal system. In this setting, multiple generations of carbonate minerals record fluid conditions and processes which control the solubility and deposition of metals, including U. This novel approach provides the first evidence of both oxidized UO 2 2+ and reduced U 4+ species in temporally distinct generations of carbonate within a single sample. Preservation of two different U oxidation states during discrete precipitation events requires U retentivity within older domains, demonstrating that the U-Pb carbonate geochronometer is robust under hydrothermal conditions. Furthermore, crystal zones with abundant fluid/vapour inclusions linked to boiling processes coincide with relatively high levels of U and favourable U/Pb. Targeting carbonate domains with these textures may therefore increase success in U-Pb geochronology. U-Pb carbonate dating combined with µXAS can track the temporal evolution of processes critical for metal deposition in long-lived and multistage hydrothermal-magmatic ore deposit settings.

58 GEOSCIENCES↗
DOE OSTI2020

Hydrous olivine alteration on Mars and Earth

In this report hydrous alteration of olivine macrocrysts in a Martian olivine phyric basalt, NWA 10416, and a terrestrial basalt from southern Colorado are examined using SEM, EPMA, TEM, and µXRD techniques. The olivines in the meteorite contain linear nanotubes of hydrous material, amorphous areas, and fluid dissolution textures quite distinct from alteration identified in other Martian meteorites. Instead, they bear resemblance to terrestrial deuteric alteration features. The presence of the hydrous alteration phase Mg‐laihunite within the olivines has been confirmed by µXRD analysis. The cores of the olivines in both Martian and terrestrial samples are overgrown by unaltered rims whose compositions match those of a separate population of groundmass olivines, suggesting that the core olivines are xenocrysts whose alteration preceded crystallization of the groundmass. The terrestrial sample is linked to deep crustal metasomatism and the “ignimbrite flare‐up” of the Oligocene epoch. The comparison of the two samples suggests the existence of an analogous relatively water‐rich magmatic reservoir on Mars.

58 GEOSCIENCES↗
DOE OSTI2024

Rehydrated glass embayments record the cooling of a Yellowstone ignimbrite

Abstract Hydration fronts penetrate 50–135 μm into glassy rhyolite embayments hosted in quartz crystals from the Mesa Falls Tuff in the Yellowstone Plateau volcanic field. The hydration fronts occur as steep enrichments that reach 2.4 ± 0.6 wt% H2O at the embayment opening, representing much higher values than interior concentrations of 0.9 ± 0.2 wt% H2O. Molecular water accounts for most of the water enrichment. Water speciation indicates the hydration fronts comprise absorbed meteoric water that modified the original magmatic composition of the rhyolitic glass. We used finite difference diffusion models to demonstrate that glass rehydration was likely produced over a few decades as the ignimbrite cooled. Such temperatures and time scales are consistent with rare firsthand observations of decadal hydrothermal systems associated with cooling ignimbrites at Mount Pinatubo (Philippines) and the Valley of Ten Thousand Smokes (Alaska).

Geology↗
DOE OSTI2025

Zinc isotope constraints on the cycling of carbon in the Bermuda mantle source

Volatile recycling and storage in the mantle transition zone (MTZ) is important for the refertilization of the upper mantle and is associated with the generation of high-µ (HIMU, where µ is 238 U/ 204 Pb) mantle. One way to probe the MTZ and the processes associated with mantle convection is to sample lavas that originate from the shallow mantle and were contaminated by upwelling from the MTZ, such as at the previously proposed shallow plume of Bermuda. Here we present the first δ 66 Zn isotopic compositions of Bermuda silica-undersaturated and silica-saturated lavas to explore the origin of the carbon-rich lithologies and the genesis of the large seamount found in the western North Atlantic Ocean. Contrasting with global δ 66 Zn data sets, our results (δ 66 Zn between 0.24 ± 0.04 and 0.41 ± 0.04) do not support direct sampling of recycled marine carbonates in the Bermuda HIMU mantle. Instead, we show that δ 66 Zn fractionation toward higher values is associated with magmatic processes and incorporation of carbon sourced from deep fluids associated with the formation of carbonatites. These carbon-rich fluids are likely sourced from the metasomatic reactions between the subducted cold slab of the Iapetus oceanic lithosphere ca. 500 Ma and the thickened continental lithospheric mantle of Pangea. Melting of this metasomatized mantle was triggered by the arrival of the Farallon slab to the eastern North American margin in the late Cenozoic via shallow convection.

Mazza, Sarah E. [Smith College, Northampton, MA (U↗
DOE OSTI2021

Interpretation of hydrothermal conditions, production-injection induced effects, and evidence for enhanced geothermal system-type heat exchange in response to >30 years of production at Roosevelt Hot Springs, Utah, USA

The Roosevelt Hot Springs hydrothermal system is located at the base of the Mineral Mountains in southwestern Utah on the eastern side of the Basin and Range. Hydrothermal activity is related to relatively recent bimodal magmatism, and the system is hosted in coarsely crystalline rock made of Oligocene–Miocene granitoids and Precambrian gneiss. The hydrothermal plume covers ~5 km 2 , with a maximum temperature of 268 °C at ~750 m depth, and a vertically extensive fault-fracture mesh east of the Opal Mound fault controls the upflow of hydrothermal fluids. Power generation (currently 38 MWe gross) began in 1984, and up through 2016, four wells were used for fluid production, and three wells were used for edge-field injection. Chemical analyses of produced fluids show that modern reservoir fluid compositions are similar to but more concentrated than those at the start of production, having near-neutral pH, total dissolved solids of 7000–10,000 mg/kg, and ionic ratios of Cl/HCO 3 ~50–100, Cl/SO4 ~50–100, and Na/K ~4–5. Chemical geothermometers indicate equilibration temperatures that mainly range between 240° and 300 °C. Early production induced a steep drop in pressure (~3.0–3.5 MPa), which was accompanied by a 250–300 m lowering of piezometric levels in wells and development of a shallow steam zone across the system. Hydrothermal fluid compositions evolved continuously in response to production-related steam-loss and injection breakthrough, which is reflected by gradual increases in chloride of up to 35% and stable isotope ratios of up to ~2‰δ 18 O and ~10‰δD. Simple mixing model calculations suggest that there has been a significant amount, ~10–20 MW th , of sustained multi-decadal heat mining and enhanced geothermal system (EGS)–type heat transfer by the injectate as it returns to the production zone. Overall, the two factors that have sustained long-term power production (currently 38 MWe gross) are the increased upflow of deep chloride water and, to a lesser extent, the mining of heat at <1 km depth.

58 GEOSCIENCES↗
DOE OSTI2020

29Si solid state NMR and Ti K-edge XAFS pre-edge spectroscopy reveal complex behavior of Ti in silicate melts

Abstract An understanding of the mechanisms of Ti is incorporation into silicate glasses and melts is critical for the field of petrology. Trace-element thermobarometry, high-field-strength element partitioning, and the physical properties of magmas are all be influenced by Ti incorporation into glasses and changes therein in response to changes in composition and temperature. In this study, we combine 29 Si solid state NMR and Ti K -edge XAFS spectroscopy to investigate how Ti is incorporated into quenched Na-silicate glasses, and the influence of Ti on the structure of silicate species in these glasses. 29 Si NMR shows that in both Ti-bearing Na 2 O•4SiO 2 (NS4) and Na 2 O•8SiO 2 (NS8) glasses, increasing the amount of Ti in the melt results in a shift of Si Q 4 peak in the 29 Si NMR spectra reflecting Ti nearest neighbors for Si in Q 4 speciation. The Ti XAFS results from NS8 glass indicate that Ti is primarily incorporated in [5]-fold coordination. At higher Ti content, there is a shift of the XAFS pre-edge feature suggesting mixing of [4]-fold Ti into the spectra. Combined, the 29 Si NMR and XAFS pre-edge data are consistent with Ti incorporation as isolated [5] Ti atoms and the formation of [5] Ti clusters at relatively low Ti concentrations, with no evidence for Ti–Na interactions as suggested by previous studies. As the Ti content increases, the Ti atoms begin to occupy 4-fold coordinated sites interacting primarily with Si in Q 4 speciation (no significant Na– [4] Ti bonding). The internal consistency of these two techniques provides a uniquely complete snapshot of the complexity of Ti incorporation in silicate melts and underlies the importance of understanding Ti incorporation mechanisms in natural magmatic systems.

Ackerson, Michael R.↗
DOE OSTI2023

Gas-mediated trace element incorporation into rhyolite-hosted topaz: A synchrotron microbeam XAS study

Magmatic gas exsolving during late-stage cooling of shallow magmas has been considered an important facilitator of low-pressure alteration and metal transport. However, the chemical properties of such gas, particularly its metal transport mechanisms and capacity, remain elusive. Trace elements in minerals produced by gas-mediated surface reaction or precipitation from gas capture details of gas composition and reaction pathways. However, interpretation of mineral trace element contents is dependent on understanding crystallographic controls on gas/mineral partitioning. Here, this work investigates the structural accommodation of As, Mn, Ga, Ge, Fe, and Ti in vapor-deposited topaz of vesicular topaz rhyolite from the Thomas Range, Utah, through single-crystal synchrotron microbeam X-ray techniques on picogram quantities of those trace elements. X-ray absorption near edge structure (XANES) data indicates that these elements are incorporated into topaz as As 5+ , Fe 3+ , Mn 3+ , Ti 4+ , Ga 3+ , and Ge 4+ . Extended X-ray absorption fine structure (EXAFS) analysis for these trace elements, compared to EXAFS of structural Al and Si, reveals that As 5+ and Ge 4+ are incorporated directly into the tetrahedral site of the topaz structure, with the octahedral site accommodating Mn 3+ , Fe 3+ , Ga 3+ , and Ti 4+ . For As 5+ and Fe 3+ , the structural impact of substitution extends to at least second neighbors (other elements were only resolvable to first neighbors). Further interpretation of the EXAFS results suggests that the substitution of Ti 4+ results in increased distortion of the octahedral site, while the other trace elements induce more uniform expansion correlating in magnitude to their ionic radius. Comparison of quantified X-ray fluorescence (XRF) data for two topaz crystals from this rhyolite reveals variable trace element concentrations for As 5+ , Fe 3+ , Ga 3+ , and Ti 4+ , reflective of a source gas undersaturated in these trace elements changing in concentration over the period of topaz deposition. The identical Ge 4+ content of the two topaz crystals suggests that Ge 4+ in the gas was buffered by the growth of another Ge 4+ -bearing phase, such as quartz. The very low Mn 3+ content in the topaz crystals does not reflect the abundance of Mn 3+ in the gas (saturation of Mn is evidenced by coexisting bixbyite). Instead, it suggests a strong Jahn-Teller inhibitory effect to the substitution of Mn 3+ for Al 3+ in the distorted octahedral site of topaz. It is proposed that exsolution of an HF-enriched gas from cooling rhyolitic magma led to local scouring of Al, Si, and trace metals from the magma. Once topaz crystals nucleated, self-catalyzed reactions that recycle HF led to continued growth of topaz.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Geothermal Fault Zone and Fluid Imaging through Joint Airborne ZTEM and Ground MT Data Inversion Analysis

This project has aimed to achieve detailed electrical resistivity resolution at geothermal reservoir scales by combining airborne natural electromagnetic (EM) field surveying (ZTEM) with ground magnetotelluric (MT) measurements to approximate an airborne MT geophysical method. MT alone is relatively expensive and may have permitting challenges in sensitive areas. Airborne ZTEM field data contains only the magnetic field, requires a background assumption, and has been limited to relatively high frequencies, thus suffering uniqueness problems. Based on proto-type 2D simulations, ZTEM ambiguities may be reduced through formal incorporation with possibly sparse ground MT soundings, which we pursued in full 3D for this project. The methodology was tested at the high-temperature Roosevelt Hot Springs geothermal system, Utah, which was considered advantageous given the near total exposure of crystalline reservoir rocks across the project area. ZTEM and ground MT survey data were acquired in 2017, subcontracted to outside parties with which we have worked in the past. These included 80 remote-referenced tensor MT soundings over the Mineral Mountains and adjacent Roosevelt Hot Spring producing geothermal system. These MT stations abut later coverage of a similar number of MT stations taken for the Utah FORGE project providing excellent total data aperture to re-solve structure beneath both project areas better than either set alone. The airborne ZTEM survey covered 704 line kilometers in E-W flight lines with a 250 m line spacing. Although this survey was timed during a maintenance-related shutdown of power production at the Roosevelt Hot Springs, other noise sources difficult to identify but including two high-voltage state-scale transmission lines compromised the ZTEM survey badly leading to unusable responses. Thus, with DOE management concurrence, the project proceeded to emphasize inversion and interpretation of the joint SubTER-FORGE MT data sets with regard to the Roosevelt Hot Springs reservoir recharge and to deep heat sources for both it and the Utah FORGE EGS project area. We also investigated the joint ZTEM-MT sampling concept with data sets from the Eleven Mile Canyon prospect area donated by the U.S. Navy (A. Sabin, PoC). Inversion of the SubTER-FORGE MT data using the HexMT 3D finite element algorithm reveals a large, low-resistivity anomaly extending sub-vertically through the depth range of the crust beneath the western Mineral Mountains. The steep conductive zone connects in the lower crust to a more tabular conductor characteristic of much of the Great Basin that generally is ascribed to current mafic magmatic underplating, hybridization and fluid release. The location of the resolved anomaly relative to the recent (0.5-0.8 Ma) eruptive centers of the Mineral Mountains implicates it as remnants of the magma body which fed these centers. This structure appears to be currently feeding heat and fluids upward into the Roosevelt Hot Springs hydrothermal system, as well as heat laterally to the FORGE project area. Separate and joint inversion models were carried out for the donated Eleven Mile Canyon MT-ZTEM data set to demonstrate concept. ZTEM only inversion showed two main alteration zones in the western portion of the project area known from geological mapping. Joint inversion including an E-W profile of MT soundings sharpened these features considerably. It also resolved in much greater detail the graben related normal faulting structure of the central project area which lies at depths exceeding the sensitivity of ZTEM alone. The sparse number of MT da-ta relative to the ZTEM required upweighting the former by a factor of several, but an exact procedure awaits future research. Our final impression is that sparse MT data can improve resolution of the subsurface over that of ZTEM alone. However, well sampled MT data are to be preferred and offer the simplicity of interpreting just one data type, and possess the superior resolution capability coming with the electric field everywhere, and from their high bandwidth.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Innovative Subsurface Learning and Hawaiian Exploration using Advanced Tomography (ISLAND HEAT) Phase 1 (Final Report)

The Innovative Subsurface Learning and Hawaiian Exploration using Advanced Tomography (ISLAND HEAT) project aims to leverage and enhance existing Hawai'i geothermal play fairway analysis (PFA) results and validate a conceptual model-driven, optimized, least-cost exploration and geophysical suite at a proven geothermal field with known permeability, the Puna Geothermal Venture, located in the Lower East Rift Zone. The abundance of geophysical scrutinization before, during, and after the 2018 Kilauea eruption at and around the broader Puna system offers a rare opportunity to elucidate geophysical expressions within a dynamic magmatic rift setting. Development of the methodology continues with its application at a second prospective site identified by the PFA, Mauna Kea, which serves as confirmation of the integrated approach.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

The geology, geochemistry, and geographic distribution of southern Idaho obsidian

The Snake River Plain Volcanic Province of Southern Idaho is host to an abundance of pyroclastic obsidian deposits that constitute the principal lithic raw material utilized by Native Americans of the region throughout the Precontact Period. Archaeologists have long looked to obsidian as a source of information on past patterns of raw material selection, landscape use, and mobility through geochemical analysis of obsidian artifacts and lithic raw materials. Formation of the Snake River Plain through cataclysmic “super-eruptions” of the Yellowstone hotspot, however, presents one of the most complex and poorly understood obsidian source-scapes of North America. To clarify the geographic distribution of Southern Idaho obsidian source groups, the INL CRMO has compiled a comprehensive reference library of geologic obsidian from 138 source locales distributed across the southern half of the state. In this report, we define 28 geochemically distinct obsidian source groups that occur at these deposits through X-ray Fluorescence spectrometry and contextualize the geographic distribution of each source group relative to the history of silicic magmatism on the Snake River Plain. This research facilitates ongoing investigations at the INL CRMO into patterns of obsidian procurement, conveyance, and exchange in Southern Idaho throughout the Precontact Period.

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

Linear Correlations of Gibbs Free Energy of REE Phosphates (Monazite, Xenotime, and Rhabdophane) and Internally Consistent Binary Mixing Properties

Rare Earth Elements (REE) phosphates (monazite, xenotime, and rhabdophane) are critical REE-bearing minerals typically formed in hydrothermal and magmatic ore deposits. The thermodynamic properties of those REE minerals are crucial to understanding the solubility, speciation, and transport of REE complexes. However, the standard-state Gibbs free energy of formation (∆G°f) values reported for these minerals in the literature vary by up to 25 kJ mol−1. Here, we present linear free energy relationships that allow the evaluation and estimation of the ∆G°f values at 25 °C and 1 bar for the three minerals from the ionic radius (rREE3+) and the non-solvation Gibbs free energy contribution to the REE3+ aqua ion (∆G°n, REE3+): ∆G°f,monazite − 399.71 rREE3+ = 1.0059 ∆G°n,REE3+ − 2522.51; ∆G°f,xenotime − 344.08 rREE3+ = 0.9909 ∆G°n,REE3+ − 2451.53; and ∆G°f,rhabdophane − 416.17 rREE3+ = 1.0067 ∆G°n, REE3+ − 2688.86. Moreover, based on the new dataset derived for REE end-members, we re-fitted the binary Margules parameter (W) from previous theoretical calculations into linear correlations: W + 0.00204 ∆G°’n,monazite = 39.3549 ∆V + 0.0641; W + 0.00255 ∆G°’n,xenotime = 25.4885 ∆V − 0.0062. The internally consistent thermodynamic properties of these REE phosphates are incorporated into the computer program Supcrtbl, which is available online at Zhu’s research website.

Pan, Ruiguang (ORCID:0000000324896545)↗
DOE OSTI2026

3D Play Fairway Analysis for Examining of Superhot Reservoir Production Scenarios

The DEEPEN (DE-risking Exploration for geothermal Plays in magmatic ENvironments) project was a multi-laboratory, international effort to reduce uncertainty and improve resource characterization in superhot geothermal systems. Building on this foundation, this work advances open-source tools designed to lower the exploration risk and cost of superhot geothermal projects while promoting transparency, reproducibility, and efficiency in exploration workflows. These tools are being tested at two key sites: (1) the Nesjavellir Geothermal Area in Iceland, where the Icelandic Deep Drilling Project (IDDP) will drill its third well, and (2) Newberry Volcano in Oregon, USA, where Mazama Energy will pilot the first superhot enhanced geothermal system (EGS). A major outcome is the creation of a modular, open-source Python framework for play fairway analysis (PFA) in 2D and 3D, called geoPFA. The PFA workflow has been expanded to produce pseudo conceptual models, and will soon be refined to assess reservoir components through integration with the thermo-hydraulic-mechanical-chemical (THMC) simulator TReactMech, to enable iterative coupling between PFA and THMC models, improving characterization of superhot systems. All three of the Icelandic Deep Drilling Project's production scenarios were analyzed via this framework: (1) a superhot deep injection well paired with conventional production wells at Nesjavellir, (2) a superhot deep production well at Nesjavellir, and (3) superhot enhanced geothermal system at Newberry Volcano. This analysis provides useful insights around conceptual modeling of these production scenarios, helping to inform decisions around which scenario is best suited for which types of environments.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Combined Structural Analysis of Core and Image Log of TGH 76-31 South East of Mt Baker, Washington State

Despite active volcanism, few geothermal energy resources have been developed in the Cascades Range. Temperature Gradient Hole 76-31 was drilled to ~440 m measured depth to probe for zones where fractures provide fluid conduits that transport deep volcanic heat to shallow depths that could support baseload, carbon neutral electrical generation. These zones were predicted by a Play Fairway assessment (PFA) of resource potential along a zone 11 km west-southwest of the summit of Mount Baker Volcano. Rock core, temperature logs, and an acoustic image log were obtained. By comparison to outcrops, the core has been initially interpreted as the Chilliwack group comprised of partially metamorphosed basaltic to andesitic volcanics. Core mapping reveals complex, steeply dipping networks of fractures and brecciation along slickensided strike slip faults; clay alteration is common in many of these structures. The majority of fractures are thoroughly healed by layers of chlorite and calcite; preservation of pyrite indicates anoxic conditions. The majority of fracture porosity resides in very dense fractures a few centimeters or less in length. The image log provides good insights into attitude of fractures that fully transect the core, but generally underestimates fracture density. The combination of complex, non-planar fracture zones containing many short fractures and healing promote misinterpretation of natural fracture attitude and density in the image log. The healing and anoxic minerals are consistent with the conductive temperature gradient measured in the well below a shallow isothermal zone, although, several fractures are open or only partially healed and resulted in fluid entries into the well. Here, the equilibrated measured temperature gradient of 64ºC/km and calculated heat flow of 145 mW/m2 is more than twice the regional average, indicating local influence of the Mt Baker magmatic system at the Little Park Creek TGH site.

Cascades↗