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

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↗

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↗

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↗

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↗

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↗

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)↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Multicontinuum Flow Models for Assessing Two-Phase Flow in Containment Science

We present a new pre-processor tool written in Python that creates multicontinuum meshes for PFLOTRAN to simulate two-phase flow and transport in both the fracture and matrix continua. We discuss the multicontinuum modeling approach to simulate potentially mobile water and gas in the fractured volcanic tuffs at Aqueduct Mesa, at the Nevada National Security Site.

58 GEOSCIENCES↗

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↗