SEARCH · Search NASA

Results for “elevation gradient”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2026

NGEE Arctic Integrated Modeling (IM2): Improved subgrid hillslope hydrologic connectivity

This data product represents the integration of new code capability for arctic tundra hillslope hydrologic processes into the Energy Exascale Earth System Model (E3SM), through the E3SM Land Model (ELM) component. This code integration is the result of collaborative effort between the NGEE Arctic project and the E3SM project. The current ELM represents water movement primarily through vertical processes, such as precipitation, canopy interception, evaporation, infiltration, and soil water movement. Lateral water movement—such as surface runoff, subsurface flow, and river transport—plays a significant role in the hydrological cycle, especially in regions with varied topography. While E3SM includes a runoff routing component representing water transport in the river network, the lateral transport of water at the subgrid scale within the land model has previously not been taken into account. With the recent development of topographic units within the ELM subgrid data structure, there is an opportunity to simulate hillslope hydrologic connectivity by introducing water transport along topographic gradients. We expect that more realistic representation of hillslope hydrologic processes will lead to improved predictions of both soil water content and river network flows. Lateral transport of water at and near the surface is represented as a sub-grid process in this new code development. Water is tracked as it moves from higher to lower elevations within a gridcell. This capability uses the nested hierarchical sub-grid scheme within ELM to connect water fluxes from sub-grid elements with higher elevation to those with lower elevation. This data record consists of a single document (pdf format) that describes the theoretical basis for the hillslope hydrology processes added to ELM, and describes the modifications made to the ELM code. The Methods section of this metadata record includes a link to the public E3SM code repository where the exact code modifications as integrated in E3SM can be accessed. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

Thornton, Peter E [ORNL] (ORCID:0000000247595158)↗

DOE OSTI2020

Mechanistic Models as Framework for Understanding Biomarker Disposition: Prediction of Creatinine‐Drug Interactions

Creatinine is widely used as a biomarker of glomerular filtration, and, hence, renal function. However, transporter‐mediated secretion also contributes to its renal clearance, albeit to a lesser degree. Inhibition of these transporters causes transient serum creatinine elevation, which can be mistaken as impaired renal function. The current study developed mechanistic models of creatinine kinetics within physiologically based framework accounting for multiple transporters involved in creatinine renal elimination, assuming either unidirectional or bidirectional‐OCT2 transport (driven by electrochemical gradient). Robustness of creatinine models was assessed by predicting creatinine‐drug interactions with 10 perpetrators; performance evaluation accounted for 5% intra‐individual variability in serum creatinine. Models showed comparable predictive performances of the maximum steady‐state effect regardless of OCT2 directionality assumptions. However, only the bidirectional‐OCT2 model successfully predicted the minimal effect of ranitidine. The dynamic nature of models provides clear advantage to static approaches and most advanced framework for evaluating interplay between multiple processes in creatinine renal disposition.

Scotcher, Daniel↗

DOE OSTI2023

Pleistocene Colorado River terraces in the canyonlands region (Utah, USA) record unsteady, transient incision and growth of the Cataract Canyon knickzone by salt tectonics

Colorado River terraces in the canyonlands region of southeastern Utah can be dated and analyzed to address the controls of incision and nature of the Colorado Plateau's largest channel-steepness anomaly, Cataract Canyon. Field correlations supported by luminescence and cosmogenic-nuclide ages on strath terraces along Meander Canyon, upstream of Cataract Canyon, reveal a complex record of unsteady incision over the past ~340 k.y. at an average rate of ~0.4 mm/yr. Both an upstream progression of rapid incision and a unique sigmoidal long-profile pattern of terraces indicate incision in response to episodic baselevel fall. Also, terraces converge downstream with the anomalously low-gradient modern channel above Cataract Canyon. Here we interpret these results as indicating that growth of the Cataract Canyon knickzone is due to an erosion-salt tectonics feedback since at least the Mid Pleistocene, which has imparted unsteady, and currently elevated, local baselevel. More broadly, the canyonlands region is marked by rapid and unsteady incision that is complicated by local geologic controls even while being absent of any regional or mantle-driven uplift.

58 GEOSCIENCES↗

DOE OSTI2023

Uniform momentum and temperature zones in unstably stratified turbulent flows

Wall-bounded turbulent flows exhibit a zonal arrangement, in which streamwise velocity organizes into uniform momentum zones (UMZs), separated by thin layers of elevated interfacial shear. While significant research efforts have focused on these structural features in neutrally stratified flows, the effects of unstable thermal stratification on UMZs and on analogous uniform temperature zones (UTZs) have not been considered previously. In this article, statistical properties of UMZs and UTZs are investigated using a suite of large eddy simulations of unstably stratified turbulent channel flow spanning weakly to highly convective conditions. When normalized by the friction velocity and stability-dependent mixing length, the mean velocity gradient based on UMZ interfacial velocity jumps and the vorticity thickness exhibits good collapse for all stabilities, establishing a link between UMZ properties and scaling predictions from Monin–Obukhov similarity theory. A similar relationship is found between UTZ properties and surface-layer scaling of the mean temperature gradient. In the mixed layer, mean UMZ depth is quasi-constant with wall-normal distance, while the deepest UTZs are found in the centre of the boundary layer. These instantaneous structures are found to be linked to the well-mixed velocity and temperature profiles in the convective mixed layer. Conditional averaging indicates that both UMZ and UTZ interfaces are associated with ejections of momentum and warm updrafts below the interface and sweeps of momentum and cool downdrafts above the interface. These results demonstrate a tangible connection between instantaneous structural features, mean properties and scaling laws in unstably stratified flows.

Mechanics↗

DOE OSTI2021

Seawater exposure causes hydraulic damage in dying Sitka-spruce trees

Sea-level rise is one of the most critical challenges facing coastal ecosystems under climate change. Observations of elevated tree mortality in global coastal forests are increasing, but the mechanism of salinity-stress induced tree mortality is currently unknown. We monitored progressive mortality and associated gas exchange and hydraulic shifts in Sitka-spruce (Picea sitchensis) trees located within a salinity gradient under an ecosystem-scale change of seawater exposure in Washington State, USA. Percentage of live foliated crown (PLFC) decreased and tree mortality increased with increasing soil salinity during the study period. A strong reduction in gas exchange and xylem hydraulic conductivity (Ks) occurred during tree death, with an increase of the percentage loss of xylem conductivity (PLC) and turgor loss point (ptlp). Hydraulic and osmotic shifts reflected that hydraulic function declined from seawater exposure and dying trees were unable to support osmotic adjustment. Constrained gas exchange was strongly related to hydraulic damage at both stem and leaf levels. Significant correlations between foliar sodium (Na+) concentration and gas exchange and key hydraulic parameters (Ks, PLC and ptlp ) suggest that cellular injury related to the toxic effects of ion accumulation impacted the physiology of these dying trees. This study provides evidence of toxic effects on cellular function that become manifest in all aspects of plant functioning, leading to unfavourable osmotic and hydraulic conditions.

Zhang, Hongxia↗

DOE OSTI2021

Interpretation of electromagnetic modes in the sub-TAE frequency range in JET plasmas with elevated monotonic q-profiles

Recent JET deuterium experiments with an advanced tokamak scenario using an internal transport barrier (ITB) have been performed to clearly observe destabilized toroidicity-induced Alfvén eigenmodes (TAEs) by fast ions; interestingly, these also exhibit unstable electromagnetic (EM) perturbations in the sub-TAE frequency range. We identify such EM perturbations to be beta-induced ion temperature gradient (BTG) eigenmodes and not beta-induced Alfvén eigenmodes nor beta-induced Alfvén acoustic eigenmodes, which are often unstable in such high-beta plasmas with high-power neutral beam injection (NBI). The BTG modes are the most unstable modes due to the high thermal ion temperature gradient related to the ITB, high thermal ion temperature compared to thermal electron temperature (high T i /T⁠ e ), and a high ion beta regime. BTG mode experimental characteristics match the analytical theory, i.e., location in the vicinity of a rational magnetic surface with a low magnetic shear, mode frequency scaling with the ion diamagnetic frequency (⁠ω$^*_i$⁠), and a coupling among Alfvén and drift waves. We also perform linear gyrokinetic simulations with validated plasma profiles and equilibrium, and find a mode kinetically driven by thermal ions with similar characteristics as the experimental BTG modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Soil water percolation and nutrient fluxes as a function of topographical, seasonal and soil texture variation in Central Amazonia, Brazil

Abstract Understanding soil water dynamics and transport of nutrients is challenging in tropical rainforests due to the uniqueness of several properties related to soils, vegetation and seasonality that make relating patterns found in temperate environments to tropical sites difficult. We address the need for better edaphic characterization in tropical environments by investigating soil water percolation rates and chemistry across topographic, soil texture and seasonal gradients in a mature tropical rainforest in Central Amazonia, Brazil. We utilized a passive wick flux meter (e.g., drainage lysimeter) to directly measure real‐time percolation fluxes at 60‐cm depth, and to sample a suite of chemical species across plateau, slope and valley topographic positions. We found percolation flux volume and chemical exports generally increase with decreasing elevation and clay content, which was lowest in the valley. Daily percolation flux was observed to be 2.39 ± 0.44 in plateau, 3.01 ± 0.50 in slope and 6.16 ± 0.83 mm in valley. Most solutes were present in small amounts of <1 mg L −1 , such as PO₄ 3− , Fe 2+ /Fe 3+ and Mn 2+ ; however, NO 3 − concentrations were >20 mg L −1 , even exceeding 100 mg L −1 in the valley. Based on additional isotopic analysis, we speculate high NO 3 − concentrations are partially an artefact of root decomposition following installation of the flux meters. The empirical relationships we show among percolation volume and nutrient exports under varying topographies and soil textures can improve Earth System Model performance by better constraining ecohydrological relationships to nutrient fluxes, which can in‐turn better illuminate the important factors that govern their behaviour.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Data from: "Snowmelt Timing Regulates Community Composition, Phenology, and Physiological Performance of Alpine Plants"

This archive contains data that were used to support conclusions drawn in “Snowmelt Timing Regulates Community Composition, Phenology, and Physiological Performance of Alpine Plants”, by Winkler et al., 2018. Data were collected throughout the 2009 growing season on Niwot Ridge, Colorado, before the site became part of the Alpine Treeline Warming Experiment (ATWE). Geospatial files are included in this archive to provide additional locational context. The files in this data package consist of five comma-separated-values (.csv) files, one keyhole markup language (.kml) file, and two ESRI shapefiles (.shp). The .csv files can be opened by Microsoft Excel, R, or any simple text-editor software, such as TextEdit (MacOS) or Notepad (Windows). The .kml files can be opened by Google Maps or Google Earth, and the .shp files are compatible with GIS softwares such as ESRI’s ArcGIS suite, and QGIS.------------------------------------------------------------------------------------------------------------------------------------------------------------------------------We asked how plant community composition, phenology, plant water relations, and photosynthetic gas exchange of alpine-restricted and wide-ranging species differ in their responses to a ca. 40-day snowmelt gradient in the Colorado Rocky mountains (Lewisia pygmaea, Sibbaldia procumbens, and Hymenoxys grandiflora were alpine-restricted and Artemisia scopulorum, Carex rupestris, and Geum rossii were wide-ranging species). To do this, we measured percent cover and flowering initiation across 20 plots varying in snowmelt timing and measured net photosynthesis and stomatal conductance in multiple individuals of each target species in these plots in 2009.As hypothesized, species richness and foliar cover increased with earlier snowmelt, due to a greater abundance of wide-ranging species present in earlier melting plots. Flowering initiation occurred earlier with earlier snowmelt for 12 out of 19 species analyzed, while flowering duration was shortened with later snowmelt for six species (all but one were wide ranging species). We observed >50% declines in net photosynthesis from July to September as soil moisture and plant water potentials declined. Early-season stomatal conductance was higher in wide-ranging species, indicating a more competitive strategy for water acquisition when soil moisture is high. Even so, there were no associated differences in photosynthesis or transpiration, suggesting no strong differences between these groups in physiology.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Mercury methylation potential and bioavailability in the sediments of two distinct aquatic systems

Here, this study explored mercury (Hg) methylation potential in two distinct aquatic systems. Fourmile Creek (FMC) was historically polluted with Hg effluents from groundwater as it is a typical gaining stream, where organic matter and microorganisms in streambed are continuously winnowed. The H02 constructed wetland only receives atmospheric Hg and is rich in organic matter and microorganisms. Both systems receive Hg from atmospheric deposition now. Surface sediments were collected from FMC and H02, spiked with inorganic Hg, and cultivated in an anaerobic chamber to stimulate microbial Hg methylation reactions. Total mercury (THg) and methylmercury (MeHg) concentrations were measured at each spiking stage. Mercury methylation potential (MMP, %MeHg in THg) and Hg bioavailability were assessed with the deployment of diffusive gradients in thin films (DGTs). During the methylation process and at the same incubation stage, FMC sediment showed faster increasing rates of %MeHg and higher MeHg concentrations than H02, demonstrating a stronger MMP in the FMC sediment. Similarly, higher Hg bioavailability was observed in FMC sediment compared to the H02 as indicated by DGT-Hg concentrations. In conclusion, the H02 wetland with high levels of organic matter and microorganisms presented low MMP. But the Fourmile Creek as a gaining stream and a historical site of Hg pollution showed strong MMP and high Hg bioavailability. A related study on microbial community activities characterized the microorganisms between FMC and H02, which is attributed to be the main reason for their different methylation capabilities. Our study further brought up the considerations on remediated sites from Hg contamination: Hg bioaccumulation and biomagnification can still be elevated and higher than the surrounding environment due to lagged changes in microbial community structures. This study supported the sustainable ecological modifications of legacy Hg contamination and raised the necessity of long-term monitoring actions even after executing a remediation plan.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Re-evaluation of the TSL for Yttrium Hydride

Yttrium hydride (YH x ) is of interest as a high-temperature moderator material in advanced nuclear reactor systems because of its superior ability to retain hydrogen at elevated temperatures. Thermal neutron scattering laws (TSL) for hydrogen bound in yttrium hydride (H-YH 2 ) and yttrium bound in yttrium hydride (Y-YH 2 ) were previously evaluated by Naval Nuclear Laboratory using the ab initio approach and released in ENDF/B-VIII.0. In that work, density functional theory, incorporating the generalized gradient approximation (GGA) for the exchange-correlation energy, was used to simulate the face-centered cubic structure of YH 2 and calculate the interatomic Hellmann-Feynman forces for a 2×2×2 supercell containing 96 atoms. Lattice dynamics calculations using PHONON were used to determine the phonon density of states. The calculated phonon density of states for H and Y in YH 2 were then used to prepare H-YH 2 and Y-YH 2 TSL evaluations, in the incoherent approximation, using the LEAPR module of NJOY. In addition, elastic scattering was assumed to be incoherent for both H and Y. While the incoherent elastic scattering approximation is appropriate for H-YH 2 , it introduces an undesirable approximation for Y-YH 2 . In this work, we re-evaluate the TSL for Y-YH 2 using FLASSH (Full Law Analysis Scattering System Hub). Y-YH 2 is evaluated using the FLASSH generalized coherent elastic scattering capability in order to capture the Bragg peaks associated with the YH 2 crystal structure which were neglected in the prior NJOY-based evaluation due to limitations in LEAPR. An experimental approach to validate the Y-YH 2 TSL using neutron transmission measurements is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss)

Abstract Background and Aims New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae. Methods We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally. Key Results Using sliding-window population genetic analyses we find that differentiation is concentrated within ‘islands’ of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap. Conclusions We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth’s climate rapidly changes.

58 GEOSCIENCES↗

DOE OSTI2025

A radar view of ice microphysics and turbulence in Arctic cloud systems

Ice microphysical processes are inherently complex because of their sensitivity to temperature and humidity, the diversity of ice crystal habits, and their interaction with supercooled liquid water (SCL) and turbulence. Long-term surface-based radar observations have been systematically used to unravel the different processes that affect ice particle growth. In this study, we present a statistical analysis of 6.5 years of Ka-band radar observations in Arctic cloud systems, combined with thermodynamic profiles derived from radiosonde measurements. For the first time, ice particle growth and sublimation – diagnosed from vertical gradients of radar reflectivity and mean Doppler velocity – are systematically mapped across a broad range of temperature and moisture conditions. These vertical gradients correspond closely to saturation levels relative to ice and exhibit a strong temperature dependence in supersaturated regions. Notably, distinct signatures near −15 °C are indicative of dendritic growth. Turbulence, quantified via the eddy dissipation rate (EDR), is most frequently observed in regions containing SCL. The co-occurrence of SCL and elevated turbulence results in significantly enhanced ice particle growth compared to conditions in which either is present alone. This work provides new observational constraints that are critical for improving the representation of ice microphysics in atmospheric models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

A research agenda for nonvascular photoautotrophs under climate change

Nonvascular photoautotrophs (NVP), including bryophytes, lichens, terrestrial algae, and cyanobacteria, are increasingly recognized as being essential to ecosystem functioning in many regions of the world. Current research suggests that climate change may pose a substantial threat to NVP, but the extent to which this will affect the associated ecosystem functions and services is highly uncertain. Here, we propose a research agenda to address this urgent question, focusing on physiological and ecological processes that link NVP to ecosystem functions while also taking into account the substantial taxonomic diversity across multiple ecosystem types. Accordingly, we developed a new categorization scheme, based on microclimatic gradients, which simplifies the high physiological and morphological diversity of NVP and world-wide distribution with respect to several broad habitat types. We found that habitat-specific ecosystem functions of NVP will likely be substantially affected by climate change, and more quantitative process understanding is required on (1) potential for acclimation, (2) response to elevated CO 2 , (3) role of the microbiome, and (4) feedback to (micro)climate. We suggest an integrative approach of innovative, multimethod laboratory and field experiments and ecophysiological modelling, for which sustained scientific collaboration on NVP research will be essential.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Tethered Balloon Observations of Vertically Resolved Aerosol Size Distributions during the U.S. DOE ARM CoURAGE and BNF Campaigns

This dataset contains vertically resolved aerosol particle number size distributions measured using a miniature Scanning Electrical Mobility Spectrometer (mSEMS) deployed aboard the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Tethered Balloon System (TBS) during the 2025 Coast-Urban-Rural Atmospheric Gradient Experiment (CoURAGE) and Bankhead National Forest (BNF) field campaigns. The mSEMS measured aerosol size distributions over approximately 10–300 nm during repeated TBS vertical profiles through the atmospheric boundary layer. The dataset provides measurements of particle number size distributions as a function of particle diameter, altitude, and time and supports investigation of aerosol vertical structure, elevated ultrafine and Aitken-mode particle layers, boundary-layer evolution, and residual-layer processes.

ARM↗

DOE OSTI2020

Plausibility of the Formose Reaction in Alkaline Hydrothermal Vent Environments

Prebiotic processes required a reliable source of free energy and complex chemical mixtures that may have included sugars. The formose reaction is a potential source of those sugars. At moderate to elevated temperature and pH ranges, these sugars rapidly decay. Here it is shown that CaCO 3 -based chemical gardens catalyze the formose reaction to produce glucose, ribose, and other monosaccharides. These thin inorganic membranes are explored as analogs of hydrothermal vent materials—a possible place for the origin of life—and similarly exposed to very steep pH gradients. Supported by simulations of a simple reaction-diffusion model, this study shows that such gradients allow for the dynamic accumulation of sugars in specific layers of the thin membrane, effectively protecting formose sugar yields. Furthermore, the formose reaction may be a plausible prebiotic reaction in alkaline hydrothermal vent environments, possibly setting the stage for an RNA world.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗

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

Groundwater level elevation and temperature data, Oct 2018-Dec 2021, Slate River Floodplain, Crested Butte, CO

This data package includes a time series of water level and temperature measurements from October 2018 to December 2021 in groundwater and surface water from the Slate River floodplain outside Crested Butte, CO, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? The data were recorded by pressure transducers were installed in four types of piezometers: 1) a network of gravel bed ("GB") drive-point piezometers with a 6-inch screen interval installed all at ~330 cm below ground surface. 2) a network of piezometers screened across the water table, used to measure evapotranspiration ("ET") using the White method. Each of these piezometers is screened along almost its entire length.3) a suite nested piezometers used the measure the vertical hydraulic gradient ("VHG") across the fines-cobble interface. Each of these piezometers uses a 6-inch screen length.4) a group of piezometers scattered across the boundaries of the floodplain, used to monitor boundary conditions ("BC") flowing into and out of the floodplain. With the exception of "SR-BD-WT", each of these piezometers is screened along its entire length. Within the data package, "FLMD.csv" describes file-level metadata and "dd.csv" defines column headers and universal terms across the dataset. The data package includes 11 "*data.csv" files, one for each piezometer type for each calendar year. Because piezometers have been added over time, not every sensor has data dating back to Oct 2018. Each "*data.csv" file has a corresponding "*_InstallationMethods.csv" file that describes the location, elevation, screen depth, sediment type and sensor metadata for each piezometer and pressure transducer.

54 ENVIRONMENTAL SCIENCES↗