Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VOLCANICS”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Methanogens acquire and bioaccumulate nickel during reductive dissolution of nickelian pyrite

ABSTRACT Nickel (Ni) is a key component of the active site metallocofactors of numerous enzymes required for methanogenesis, including [NiFe]-hydrogenase, carbon monoxide dehydrogenase, and methyl CoM reductase, leading to a high demand for Ni among methanogens. However, methanogens often inhabit euxinic environments that favor the sequestration of nickel as metal-sulfide minerals, such as nickelian pyrite [(Ni,Fe)S 2 ], that have low solubilities and that are not considered bioavailable. Recently, however, several different model methanogens ( Methanosarcina barkeri , Methanococcus voltae , Methanococcus maripaludis ) were shown to reductively dissolve pyrite (FeS 2 ) and to utilize dissolution products to meet iron and sulfur biosynthetic demands. Here, using M. barkeri Fusaro, and laboratory-synthesized (Ni,Fe)S 2 that was physically isolated from cells using dialysis membranes, we show that trace nickel (<20 nM) abiotically solubilized from the mineral can support methanogenesis and limited growth, roughly fivefold less than the minimum concentration known to support methanogenesis. Furthermore, when provided direct contact with (Ni,Fe)S 2 , M. barkeri promoted the reductive dissolution of (Ni,Fe)S 2 and assimilated solubilized nickel, iron, and sulfur as its sole source of these elements. Cells that reductively dissolved (Ni,Fe)S 2 bioaccumulated approximately fourfold more nickel than those grown with soluble nickel and sulfide but had similar metabolic coupling efficiencies. While the mechanism for Ni uptake in archaeal methanogens is not known, homologs of the bacterial Nik uptake system were shown to be ubiquitous across methanogen genomes. Collectively, these observations indicate that (Ni,Fe)S 2 is bioavailable in anoxic environments and that methanogens can convert this mineral into nickel-, iron-, and sulfur-containing metalloenzymes to support methanogenesis and growth. IMPORTANCE Nickel is an essential metal, and its availability has changed dramatically over Earth history due to shifts in the predominant type of volcanism in the late Archean that limited its availability and an increase in euxinic conditions in the early Proterozoic that favored its precipitation as nickel sulfide minerals. Observations presented herein indicate that the methanogen, Methanosarcina barkeri , can acquire nickel at low concentration (<20 nM) from soluble and mineral sources. Furthermore, M. barkeri was shown to actively reduce nickelian pyrite; use dissolution products to meet their iron, sulfur, and nickel demands; and bioaccumulate nickel. These data help to explain how M. barkeri (and possibly other methanogens and anaerobes) can acquire nickel in contemporary and past anoxic or euxinic environments.

(NiFe)S2↗
DOE OSTI2021

Ecological Dichotomies Arise in Microbial Communities Due to Mixing of Deep Hydrothermal Waters and Atmospheric Gas in a Circumneutral Hot Spring

Understanding the source and availability of energy capable of supporting life in hydrothermal environments is central to predicting the ecology of microbial life on early Earth when volcanic activity was more widespread. Little is known of the substrates supporting microbial life in circumneutral to alkaline springs, despite their relevance to early Earth habitats.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Continental rift evolution and drainage reorganization along the Dead Sea rift since the Miocene

The Dead Sea fault is a section of the Arabian-African plate boundary. Widespread field relations indicate that three major drainage systems (stages) occupied the landscape west of the Dead Sea fault since its initiation at ca. 20 Ma. Specifically, (1) an early to middle Miocene drainage system, only minorly reconfigured by the fault (all sediments of this system belong to the Hazeva Formation); (2) a late Miocene to early Pleistocene fault-parallel drainage system named Paran-Neqarot (all sediments of this system belong to the Arava and Zehiha Formations; and (3) the early Pleistocene to present drainage configuration. The temporal and spatial frameworks of drainage stage 1 are generally constrained by radiometric dating of interfingering volcanic units, and the onset and temporal and spatial frameworks of drainage stage 3 are well constrained by cosmogenic 10 Be surface exposure ages. The timing and longevity of the rift-parallel drainage (stage 2) have until now been evasive to direct dating. The overall time gap between stages 1 and 3 is ~12–13 million years. Thus, an early age (within this time gap) of stage 2 would imply an immediate response of drainage reorganization to rift tectonics, while a later age of this drainage system would imply a delayed response. We present 11 10 Be- 26 Al cosmogenic burial ages of alluvial and colluvial units related to the fault-parallel drainage system (stage 2), which collectively constrain the time of deposition of the Arava Formation sediments in the central Negev to ca. 8 Ma. The general lack of stratigraphic order, together with the large dispersion of ages both across and within the sampling sites, attests to significant recycling of sediments from drainage stage 1 into the Arava Formation deposits. The termination of Arava and Zehiha Formation sediment deposition at ca. 1.8 Ma was determined previously using cosmogenic exposure ages of desert pavements that cover the formations. Combining the previously published data with our new data, we established the longevity and character of the Paran-Neqarot drainage system. In conlcusion, this framework highlights the temporal aspect of drainage system build-up and collapse as it responded to transform and extensional plate boundary tectonics during the Neogene.

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 OSTI2021

The spatial and temporal evolution of the Portland and Tualatin forearc basins, Oregon, USA

The Portland and Tualatin basins are part of the Salish-Puget-Willamette Lowland, a 900-km-long, forearc depression lying between the volcanic arc and the Coast Ranges of the Cascadia convergent margin. Such inland seaways are characteristic of warm, young slab subduction. We analyzed the basins to better understand their evolution and relation to Coast Range history and to provide an improved tectonic framework for the Portland metropolitan area. We model three key horizons in the basins: (1) the top of the Columbia River Basalt Group (CRBG), (2) the bottom of the CRBG, and (3) the top of Eocene basement. Isochore maps constrain basin depocenters during (1) Pleistocene to mid-Miocene time (0–15 Ma), (2) CRBG (15.5–16.5 Ma), and (3) early Miocene to late Eocene (ca. 17–35 Ma) time. Results show that the Portland and Tualatin basins have distinct mid-Miocene to Quaternary depocenters but were one continuous basin from the Eocene until mid-Miocene time. A NW-striking gravity low coincident with the NW-striking, fault-bounded Portland Hills anticline is interpreted as an older graben coincident with observed thickening of CRBG flows and underlying sedimentary rocks. Neogene transpression in the forearc structurally inverted the Sylvan-Oatfield and Portland Hills normal faults as high-angle dextral-reverse faults, separating the Portland and Tualatin basins. An eastward shift of the forearc basin depocenter and ten-fold decrease in accommodation space provide temporal constraints on the emergence of the Coast Range to the west. Clockwise rotation and northward transport of the forearc is deforming the basins and producing local earthquakes beneath the metropolitan area.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Tracking the spatial and temporal links between late Cenozoic extension and magmatism in the Death Valley region (California, USA)

The relationship between late Cenozoic magmatism and extension in the central Basin and Range province (western United States) is complex, necessitating high-precision geochronology to understand its spatiotemporal connections. In the Death Valley region (California), the lack of high-precision U-Pb zircon ages has limited our understanding of the timing of pluton formation and its links to regional extension. We present new high-precision chemical abrasion-isotope dilution-thermal ionization mass spectrometry 206 Pb/ 238 U zircon ages and trace element analyses for eight Death Valley plutons. Our findings reveal three distinct phases of intrusive magmatism: (1) emplacement of shallow rapakivi granites at 13.2 Ma, (2) construction of the mid-crustal Black Mountains intrusive complex at 11.3 Ma, and (3) late emplacement of shallow, compositionally diverse intrusions at 8.2 Ma. A gap in zircon crystallization between 10 Ma and 8.2 Ma coincides with exhumation of the Black Mountains and a transition from sill to dike emplacement. The dominance of rapakivi granites in the Death Valley region, which is rare among Cenozoic granitoids, is likely a result of rapid crustal extension that induces adiabatic decompression. A comparison of the timing of volcanism, plutonism, and tectonic events in Death Valley reveals that intrusive magmatism closely tracks the locus of extension, underscoring the plutonic record as a vital link for understanding regional tectonics and changes in plate boundary dynamics during this period.

58 GEOSCIENCES↗
DOE OSTI2022

Earth, Wind, Fire, and Pollution: Aerosol Nutrient Sources and Impacts on Ocean Biogeochemistry

A key Earth system science question is the role of atmospheric deposition in supplying vital nutrients to the phytoplankton that form the base of marine food webs. Industrial and vehicular pollution, wildfires, volcanoes, biogenic debris, and desert dust all carry nutrients within their plumes throughout the globe. In remote ocean ecosystems, aerosol deposition represents an essential new source of nutrients for primary production. The large spatiotemporal variability in aerosols from myriad sources combined with the differential responses of marine biota to changing fluxes makes it crucially important to understand where, when, and how much nutrients from the atmosphere enter marine ecosystems. This review brings together existing literature, experimental evidence of impacts, and new atmospheric nutrient observations that can be compared with atmospheric and ocean biogeochemistry modeling. We evaluate the contribution and spatiotemporal variability of nutrient-bearing aerosols from desert dust, wildfire, volcanic, and anthropogenic sources, including the organic component, deposition fluxes, and oceanic impacts.

Marine & Freshwater Biology↗
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 OSTI2021

A Frequency-Domain-Based Algorithm for Detecting Microseismicity Using Dense Surface Seismic Arrays

We propose a new frequency-domain-based algorithm for detecting small-magnitude seismic events using dense surface seismic arrays. Our proposed method takes advantage of the high energy carried by S waves, and approximate known source locations, which are used to rotate the horizontal components to obtain the maximum amplitude. By surrounding the known source area with surface geophones, we achieve a favorable geometry for locating the detected seismic events with the backprojection method. To test our new detection method, we used a dense circular array, consisting of 151 5 Hz three-component geophones, over a 5 km aperture that was in operation at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) in southcentral Utah. We apply the new detection method during a small-scale test injection phase at FORGE, and during an aftershock sequence of an Mw 4.1 earthquake located ~30 km north of the geophone array, within the Black Rock volcanic field. We are able to detect and locate microseismic events (Mw<0) during injections, despite the high level of anthropogenic activity, and several aftershocks that are missing from the regional catalog. By comparing our method with known algorithms that operate both in the time and frequency domain, we show that our proposed method performs better in the case of the FORGE injection monitoring, and equally well for the off-array aftershock sequence. Our new method has the potential to improve microseismic event detections even in extremely noisy environments, and the proposed location scheme serves as a direct discriminant between true and false detections.

Geochemistry & Geophysics↗
DOE OSTI2024

Joint Inversion of Regional Waveform, First-Motion Polarity, and Synthetic Aperture Radar Surface Displacement for the Fourth and Sixth North Korean Declared Nuclear Explosions

Here, this study analyzed the Democratic People’s Republic of Korea’s (DPRK) fourth (DPRK4, 6 January 2016 M w 4.49) and sixth (DPRK6, 7 September 2017 M w 5.2) declared nuclear tests, employing a joint seismic and Interferometric Synthetic Aperture Radar (InSAR) inversion to improve understanding of these events and enhance moment tensor (MT) inversion capabilities. The recent efforts have focused on employing seismic waveform and InSAR geodetic deformation data separately to analyze these and the previous nuclear tests (e.g., Chiang et al., 2018; Myers et al., 2018; Wang et al., 2018). Building upon our previous work (Chi-Durán et al., 2021), we performed a joint regional waveform, first-motion (FM) polarity, and surface displacement inversion, which demonstrated improved source-type discrimination, a revised MT solution with reduced scalar moment uncertainty, and an independently constrained location. In this article, we build on the previous results for DPRK6 by including an analysis using a four-layered velocity model with free-surface topography to compute the near-source static deformation Green’s functions. The model consists of a 50 m basalt layer (⁠V P = 2.07 km/s, V S = 1.2 km/s⁠), a 250 m stratified volcanic deposit layer (⁠V P = 1.73 km/s, V S = 1.0 km/s⁠), a 700 m weathered granodiorite layer (⁠V P = 2.5 km/s, V S = 1.3 km/s⁠⁠), and a granodiorite half-space (⁠V P = 5.35 km/s, V S = 3.09 km/s⁠⁠). The half-space shares the velocity of the regional MDJ2 velocity model (Ford et al., 2010), which has proven effective for waveform inversion in the region. This model considers the range of reported values for various lithologies and weathering effects. Our findings show that using the layered velocity model enhances the recovery of source location and depth for both the explosions by improving fits and reducing uncertainties. The joint inversion also improves source-type discrimination and better constrains the scalar seismic moment necessary for downstream yield estimation.

58 GEOSCIENCES↗
DOE OSTI2024

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure of the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

Here we present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. The new model enables fully 3D simulations of complete three-component waveforms and provides improved fits that were not possible with previous models. We inverted over 32,000 waveforms from 192 earthquakes recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region with highly uneven coverage. Inversion iterations proceeded from the period band 50–100 s in six stages and 54 total iterations reducing the minimum period to 30 s. Our final model, MESWA, improves waveform fits compared to the starting and other models for both the data used in the inversion and an independent validation set of 66 events. Restitution tests indicate that the model resolves features in the central part of the model to depths of about 150 km. The new model reveals tectonic features imaged by other studies and methods but in a new holistic model of anisotropic shear and compressional wavespeeds (V S and V P , respectively) covering a larger domain with smaller scale length and amplified features. Examples include low crustal V S in the Tethyan belt and low mantle VS following divergent (Gulf of Aden, Red Sea) and transform (Dead Sea fault) margins of the Arabian plate. Low V S is imaged below Cenozoic volcanic centers of the Mecca–Madina–Nafud Line, Arabian Peninsula, and the Türkiye–Iran border region. Elevated V S tracks Makran subduction under southeast Iran with near vertical dip. MESWA could be used as a starting model for further improvements, say, using waveforms from in-country seismic networks that are not currently openly available and/or smaller-scale studies targeting a shorter period. The model could be used to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗
DOE OSTI2021

Humming Trains in Seismology: An Opportune Source for Probing the Shallow Crust

Abstract Seismologists are eagerly seeking new and preferably low-cost ways to map and track changes in the complex structure of the top few kilometers of the crust. By understanding it better, they can build on what is known regarding important, practical issues. These include telling us whether imminent earthquakes and volcanic eruptions are generating telltale underground signs of hazard, about mitigation of induced seismicity such as from deep injection of wastewater, how the Earth and its atmosphere couple, and where accessible natural resources are. Passive seismic imaging usually relies on blind correlations within extended recordings of Earth’s ceaseless “hum” or coda of well-mixed, small vibrations. In this article, we propose a complementary approach. It is seismic interferometry using opportune sources—specifically ones not stationary in time and moving in a well-understood configuration. Its interpretation relies on an accurate understanding of how these sources radiate seismic waves, precise timing, careful placement of pairs of listening stations, and seismic phase differentiation (surface and body waves). Massive freight trains were only recently recognized as such a persistent, powerful cultural (human activity-caused) seismic source. One train passage may generate a tremor with an energy output of a magnitude 1 earthquake and be detectable for up to 100 km from the track. We discuss the source mechanisms of train tremors and review the basic theory on sources. Finally, we present case studies of body- and surface-wave retrieval as an aid to mineral exploration in Canada and to monitoring of a southern California fault zone. We believe noise recovery from this new signal source, together with dense data acquisition technologies such as nodes or distributed acoustic sensing, will deeply transform our ability to monitor activity in the shallow crust at sharpened resolution in time and space.

Geochemistry & Geophysics↗
DOE OSTI2025

AK112: Full Waveform Inversion Tomography of Alaska Improves Waveform Fits While Imaging Crustal, Mantle, and Slab Structure

We report a full waveform inversion tomography model of Alaska and the surrounding regions, inferring radially anisotropic shear and isotropic compressional wavespeeds by fitting complete waveforms from 120 regional earthquakes. Our multiscale approach inverted time–frequency phase misfits (maximum period of 100 s), starting with a minimum period of 40 s and ending at 20 s in 7 stages and 112 total iterations. The model (AK112) was evaluated by computing the misfits for 36 independent validation events. We find that misfit reductions were large and equal (∼55%) for both the inversion and validation data sets, providing confidence in the model. AK112 also provides much better waveform fits compared to other reported models for the region, including an isotropic version of itself, highlighting the importance of anisotropy. The model resolves known crustal, upper mantle, and slab structure to depths of 100 km with new detail: sedimentary basins in the Alaskan Shelf, Cook Inlet, and Colville basins, among others; discontinuous lithospheric structure across major terrane boundaries; and subducting slab geometry and back‐arc volcanic sources. In addition to tectonic interpretations, the model enables full waveform simulations for long‐period earthquake ground motions and source characterization (e.g., moment tensor and finite‐fault inversion).

Rodgers, Arthur [Lawrence Livermore National Labor↗
DOE OSTI2020

IC Project: w19_hossfault “Modelling of stick-slip behavior in sheared granular fault gouge & nonlinear elasticity behavior in cracked solid”

Low frequency earthquakes, non-volcanic tremor, and acoustic emissions are examples of weak seismic signals that may help detect major seismic events, i.e., earthquakes. In a laboratory setting, in which stick-slip events are simulated, acoustic emissions are detected far from the stick-slip events. In both, field and laboratory scale cases, the acoustic emission signals are sourced or detected in a volume remote from the volume that spawns the earthquake. Therefore, it is relevant to establish the causal relationship between signals detected on passive, remote monitors and the dynamics of the elastic structures that launch important seismic events. The work conducted under this IC allocation allowed us to utilize a numerical model that let us follow this causality, i.e., examine and connect the dynamics in a granular system (fault gouge), to signals detected on passive remote monitors. It was demonstrated that stress chains are key in the dynamics of granular systems and that their evolution is the source of the acoustic emission.

58 GEOSCIENCES↗
DOE OSTI2020

Advancing Measurements and Understanding of the Rate and Structure of Atmospheric Warming

The rate and geographic pattern of atmospheric warming are key indicators of historical climate change and play a role in modulating future changes in Earth’s climate. The past and future evolution of atmospheric temperature is simulated using global climate models while observational estimates of past warming are derived from satellite microwave measurements. Individual climate models simulate widely varying rates of past and future atmospheric warming and most model simulations exhibit greater tropospheric (lowest ~10 km of atmosphere) temperature change than satellite observations between 1979 and 2020. This project examined intermodel differences in the pattern of atmospheric warming and how these differences influence climate feedbacks that can amplify or damp the rate of global surface warming. A key result is that climate model representation of the current climate can influence climate feedbacks and the simulation of future changes in climate. Another focus of this project was to analyze how models respond to different input datasets, such as volcanic aerosols or sea surface temperature (in atmosphere-only simulations). We find that different prescribed inputs can affect simulated changes in atmospheric temperature, even though “input uncertainty” is often unconsidered in model-observational comparisons. A final focus of this research project was to consider the influence of natural internal climate variability on satellite era changes in climate. Although the Earth is warming substantially due to anthropogenic emissions of greenhouse gases, the observed warming rate can be modulated by natural variations in the Earth’s climate. We find that natural climate variability has slowed the rate of tropical tropospheric warming, which explains model-satellite differences in the rate of tropospheric warming.

58 GEOSCIENCES↗
DOE OSTI2020

Acid and heat stable enzymes for economical fuel production from biomass

CinderBio’s vision is to harness three billion years of biology to transform many industrial practices into sustainable processes with superior performance. Environments such as volcanic hot springs have evolved extremophilic microbes and enzymes that are naturally heat and acid stable. CinderBio was founded based on technologies developed at Lawrence Berkeley National Laboratory. Our core technologies are the molecular tools and know how to efficiently express extremophilic genes in an extremophilic host such that they retain their natural heat and acid stability. During this project, we tested our enzymes at the Advanced Biofuels and Bioproducts Process Development Unit (ABPDU), automated our fermentation and downstream processing, and built a comprehensive technoeconomic assessment (TEA) model for our manufacturing costs.

09 BIOMASS FUELS↗
DOE OSTI2021

Shocked Transport: Experiments to Study Dynamic Particle Behavior under Varying Volume Fraction Conditions

Shocks interacting with densely-packed porous media is common in volcanic eruptions or in engineering applications such as ordinance. Large-scale, exploratory experiments have been performed by Frost and collaborators examining the behavior of a bed of particles subjected to a detonation wave. The detonation replaces the central explosive with a high-density, high-temperature gas that propagates outward. A shock wave travels through the particle bed, followed by a contact discontinuity between the shocked air and the explosive products. In the early times, the detonation wave causes many of the particles to collide and compact with each other. As the particles disperse, the experiments demonstrate a series of large-scale jets or fingers that form in the advancing particle front.

42 ENGINEERING↗
DOE OSTI2021

Discovering Hidden Geothermal Signatures using Unsupervised Machine Learning

Discovering hidden geothermal resources is a very challenging task. It requires the mining of large datasets, including various diverse data attributes representing subsurface hydrogeological and geothermal conditions. The commonly used Play Fairway Analysis (PFA) typically relies on subject-matter expertise to analyze site or regional data to estimate geothermal conditions and prospectivity. Here, we demonstrate an alternative approach based on machine learning (ML) to process a geothermal dataset of Southwest New Mexico (SWNM). The study region includes low- and medium-temperature hydrothermal systems. However, most of these systems are poorly characterized because of insufficient existing data and limited past explorative studies. This study aims to discover hidden patterns and relationships in the SWNM geothermal dataset to better understand regional hydrothermal conditions. This is achieved by applying an unsupervised machine learning algorithm based on non-negative matrix factorization coupled with customized k-means clustering (NMFk). NMFk can automatically identify (1) hidden (latent) signatures characterizing datasets, (2) the optimal number of these signatures, (3) dominant data attributes associated with each signature, and (4) spatial distribution of the extracted signatures. Here, NMFk is applied to analyze 18 geological, geophysical, hydrogeological, geothermal attributes at 44 locations in SWNM. NMFk successfully finds data patterns and identifies the spatial associations of hydrothermal signatures with the four physiographic provinces in SWNM (Colorado Plateau, Volcanic Field, Basin and Range, and the Rio Grande rift). The algorithm identified up to 5 hydrothermal signatures in the SWNM datasets that differentiate between low- and medium-temperature hydrothermal systems in different provinces. Also, the algorithm identifies two medium-temperature hydrothermal systems in SWNM that require further exploration for geothermal resource development. Based on our analyses, 12 of the attributes are important to identify medium-temperature hydrothermal systems, and the remaining six attributes are critical to characterize low-temperature hydrothermal systems. Based on the obtained results, we identify potential physiographic provinces for further exploration to characterize them as geothermal resources. The resulting NMFk model can be applied to predict geothermal conditions and their uncertainties at new SWNM locations based on limited data from unexplored areas.

58 GEOSCIENCES↗