Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed hydrologic model”

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 145 records · Page 8

Real-time deep-learning inversion of seismic full waveform data for CO 2 saturation and uncertainty in geological carbon storage monitoring

Deep-learning inversion has recently drawn attention in geological carbon storage research due to its potential of imaging and monitoring carbon storage in real time, significantly improving efficiency and safety of carbon storage operations. We present a deep-learning full waveform inversion method that after the neural network has been trained can image CO 2 saturation and its uncertainty in real time. Our deep-learning inversion method is based on the U-Net architecture with the neural network trained on pairs of synthetic seismic data and CO 2 saturation models. Accordingly, our training establishes a mapping relationship between seismic data and CO 2 saturation models and once fully trained directly estimates CO 2 saturation as a function of subsurface location. We further quantify uncertainties of CO 2 saturation estimates using the Monte Carlo dropout method and a bootstrap aggregating method. For this proof-of-concept study, the CO 2 training models and data are derived from the Kimberlina 1.2 model, a hypothetical 3D geological carbon storage model that is constructed based on various geological and hydrological data from the Southern San Joaquin Basin, California. We perform deep-learning inversion experiments using noise-free and noisy training and test data sets and compare the results. Our modelling experiments show that (1) the deep-learning inversion can estimate 2D distributions of CO 2 fairly well even in the presence of Gaussian random noise and (2) both CO 2 saturation imaging and uncertainty quantification can be done in real time. Our results suggest that the deep-learning inversion method can serve as a robust real-time monitoring tool for geological carbon storage and/or other time-varying reservoir/aquifer properties that result from injection, extraction, and/or other subsurface transport phenomena.

58 GEOSCIENCES↗

River Dynamics Control Transit Time Distributions and Biogeochemical Reactions in a Dam-Regulated River Corridor

Residence Time Distributions (RTDs) exert an important control on biogeochemical translation in the hydrological systems. Previous tracer studies have revealed that RTDs often followed time-invariant exponential, lognormal, or heavy-tailed RTDs that have power-law behavior for long tails in headwater or low-order streams. However, there is a recognition that RTDs can be more complicated and time-variable in response to dynamic hydrological forcing. Here, we use particle tracking to estimate RTDs along the Hanford Reach of the Columbia River and to quantify the influences of river stage fluctuations. Particle tracking is conducted using the velocity field from high-resolution 3D groundwater flow simulations. The effects of hydrological forcing on the residence time distribution were evaluated by varying river flow boundary conditions and releasing particles in different time windows. Our results revealed that dynamic stage fluctuations created rapidly changing losing-gaining conditions in the river and led to highly transient RTDs, which contributes to multiple modes of RTDs. Dam-induced high-frequency (sub-daily) flow variation contributes to the short-time (sub-daily) component of the RTDs. Deviation of the reactant consumption under the single-mode assumption compared to the multimodal RTDs is relatively small (~5%) and appears when the Damköhler number was close to one.More specifically, high-frequency (daily) dam-induced stage variations have higher impacts on biogeochemical reactions with faster reaction rate, and can potentially increase the reactant consumption by 26.68%. Based on these findings, we suggest that current river basin models could be improved by including bank storage and more complex RTDs influenced by both short and long term river stage fluctuations.

54 ENVIRONMENTAL SCIENCES↗

WHOLESCALE - Water & Hole Observations Leverage Effective Stress Calculations And Lessen Expenses (Final Technical Report 2020 - 2024)

The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States. To reach this goal, the WHOLESCALE team has developed a methodology to incorporate and interpret data from four methods of measurement into a multi-physics model that couples thermal, hydrological, and mechanical (T H-M) processes. The WHOLESCALE team has applied this methodology at the San Emidio geothermal field, located ~100 km north of Reno, Nevada in the northwestern Basin and Range province. The WHOLESCALE team includes 30 individuals working at two universities, two national laboratories, and one industry partner. Two master-degree students and five post-doctoral researchers have gained professional experience and earned partial financial support via the WHOLESCALE project. The WHOLESCALE team has taken advantage of the perturbations created by changes in pumping operations during planned shutdowns in 2016, 2021, and 2022 to infer temporal changes in the state of stress in the geothermal system at San Emidio, Nevada, U.S. The WHOLESCALE results support the working hypothesis that increasing pore-fluid pressure reduces the effective normal stress acting across fault zones. During normal operations, pumping in deep production wells decreases fluid pressures and thus increases the effective normal stresses on faults, reducing microseismicity. During planned shutdowns, the cessation of production increases pore-fluid pressure and reduces effective normal stress. The WHOLESCALE products generated during the 4-year period between 2020 and 2024 include: three articles published in the open-access, peer-reviewed scientific literature, two master’s theses, 20 presentations or papers at scientific conferences, and 17 data sets available on public repositories. The WHOLESCALE project has been completed in two phases that included three performance periods separated by two Go/No-go Stage Gate Reviews. Tasks were classified by data type (i.e., Geologic Structure, Borehole, Geodesy, Hydrology, Seismology, and Modeling). The first phase of the project started July 31, 2020 and included ongoing project coordination (Task 1), a project kickoff (Task 2), analysis of existing data (Task 3), development of the initial stress model & deployment design (Task 4), and Go/No-go Decision Point #1 (Task 5). Phase II began with implementing the 2022 deployment (Task 6), followed by Go/No-go Decision Point #2 (Task 7) The remainder of Phase II consisted of analyzing data collected during deployment (Task 8), calibration of the stress model on all observations (Task 9), and the Final Review (August 23, 2024) & Reporting (Task 10).

15 GEOTHERMAL ENERGY↗

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗

Multi-scale planning model for robust urban drought response

Increasingly severe droughts are straining municipal water resources and jeopardizing urban water security, but uncertainty in their duration, frequency, and intensity challenges drought planning and response. We develop the Drought Resilient Interscale Portfolio Planning model (DRIPP) to generate optimal planning responses to urban drought. DRIPP is a generalizable multi-scale framework for optimizing dynamic planning strategies of long-term infrastructure deployment and short-term drought response. It integrates climate and hydrological variability with high-fidelity representations of urban water distribution, available technology options, and demand reduction measures to yield robust and cost-effective water supply portfolios that are location-specific. We apply DRIPP in Santa Barbara, California to assess how least cost water supply portfolios vary under different drought scenarios and identify portfolios that are robust across drought scenarios. In Santa Barbara, we find that drought intensity, not duration or frequency, drives cost increases, reliability risk, and regret of overbuilding infrastructure. Under uncertain drought conditions, a diversified technology portfolio that includes both rapidly deployable, decentralized technologies alongside larger centralized technologies minimizes water supply cost while maintaining high robustness to climate uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Solving high-dimensional inverse problems using amortized likelihood-free inference with noisy and incomplete data

Here, we present a likelihood-free probabilistic inversion method based on normalizing flows for high-dimensional inverse problems. The proposed method is composed of two complementary networks: a summary network for data compression and an inference network for parameter estimation. The summary network encodes raw observations into a fixed-size vector of summary features, while the inference network generates samples of the approximate posterior distribution of the model parameters based on these summary features. The posterior samples are produced in a deep generative fashion by sampling from a latent Gaussian distribution and passing these samples through an invertible transformation. We construct this invertible transformation by sequentially alternating conditional invertible neural network and conditional neural spline flow layers. The summary and inference networks are trained simultaneously. We apply the proposed method to an inversion problem in groundwater hydrology to estimate the posterior distribution of the log-conductivity field conditioned on spatially sparse time-series observations of the system’s hydraulic head responses. The conductivity field is represented with 706 degrees of freedom in the considered problem. Comparison with the likelihood-based iterative ensemble smoother PEST-IES method demonstrates that the proposed method accurately estimates the parameter posterior distribution and the observations’ predictive posterior distribution at a fraction of the inference time of PEST-IES.

conditional invertible neural network↗

Representing plant diversity in land models: An evolutionary approach to make “Functional Types” more functional

Plants are critical mediators of terrestrial mass and energy fluxes, and their structural and functional traits have profound impacts on local and global climate, biogeochemistry, biodiversity, and hydrology. Yet, Earth System Models (ESMs), our most powerful tools for predicting the effects of humans on the coupled biosphere-atmosphere system, simplify the incredible diversity of land plants into a handful of coarse categories of "Plant Functional Types" (PFTs) that often fail to capture ecological dynamics such as biome distributions. The inclusion of more realistic functional diversity is a recognized goal for ESMs, yet there is currently no consistent, widely accepted way to add diversity to models, that is, to determine what new PFTs to add and with what data to constrain their parameters. Here we review approaches to representing plant diversity in ESMs and draw on recent ecological and evolutionary findings to present an evolution-based functional type approach for further disaggregating functional diversity. Specifically, the prevalence of niche conservatism, or the tendency of closely related taxa to retain similar ecological and functional attributes through evolutionary time, reveals that evolutionary relatedness is a powerful framework for summarizing functional similarities and differences among plant types. We advocate that Plant Functional Types based on dominant evolutionary lineages ("Lineage Functional Types") will provide an ecologically defensible, tractable, and scalable framework for representing plant diversity in next-generation ESMs, with the potential to improve parameterization, process representation, and model benchmarking. We highlight how the importance of evolutionary history for plant function can unify the work of disparate fields to improve predictive modeling of the Earth system.

59 BASIC BIOLOGICAL SCIENCES↗

Temporal, Spatial, and Temperature Controls on Organic Carbon Mineralization and Methanogenesis in Arctic High-Centered Polygon Soils

Warming temperatures in continuous permafrost zones of the Arctic will alter both hydrological and geochemical soil conditions, which are strongly linked with heterotrophic microbial carbon (C) cycling. Heterogeneous permafrost landscapes are often dominated by polygonal features formed by expanding ice wedges: water accumulates in low centered polygons (LCPs), and water drains outward to surrounding troughs in high centered polygons (HCPs). These geospatial differences in hydrology cause gradients in biogeochemistry, soil C storage potential, and thermal properties. Presently, data quantifying carbon dioxide (CO 2 ) and methane (CH 4 ) release from HCP soils are needed to support modeling and evaluation of warming-induced CO 2 and CH 4 fluxes from tundra soils. This study quantifies the distribution of microbial CO 2 and CH 4 release in HCPs over a range of temperatures and draws comparisons to previous LCP studies. Arctic tundra soils were initially characterized for geochemical and hydraulic properties. Laboratory incubations at −2, +4, and +8°C were used to quantify temporal trends in CO 2 and CH 4 production from homogenized active layer organic and mineral soils in HCP centers and troughs, and methanogen abundance was estimated from mcrA gene measurements. Results showed that soil water availability, organic C, and redox conditions influence temporal dynamics and magnitude of gas production from HCP active layer soils during warming. At early incubation times (2–9 days), higher CO 2 emissions were observed from HCP trough soils than from HCP center soils, but increased CO 2 production occurred in center soils at later times (>20 days). HCP center soils did not support methanogenesis, but CH 4 -producing trough soils did indicate methanogen presence. Consistent with previous LCP studies, HCP organic soils showed increased CO 2 and CH 4 production with elevated water content, but HCP trough mineral soils produced more CH 4 than LCP mineral soils. HCP mineral soils also released substantial CO 2 but did not show a strong trend in CO 2 and CH 4 release with water content. Knowledge of temporal and spatial variability in microbial C mineralization rates of Arctic soils in response to warming are key to constraining uncertainties in predictive climate models.

59 BASIC BIOLOGICAL SCIENCES↗

Foliar element determination from field survey in association with the National Ecological Observatory Network Airborne Observation Platform survey, East River, Colorado 2018

The purpose of this dataset is to support research aimed at understanding the coupling between hydrologic and biogeochemical processes at watershed scale, particularly the relationship between aboveground vegetation characteristics and subsurface soil properties. These data are intended to inform and calibrate models of catchment-scale biogeochemical fluxes, including rock-derived nutrient cycling, and they were procured to address the following questions: (1) What is the distribution of vegetation characteristics across the study catchments? (2) Are foliar concentrations of rock-derived nutrients related to underlying lithology and soil availability, or are these signals masked by biotic nutrient cycling and retention processes?This data package contains foliar elemental data collected during the 2018 National Ecological Observatory Networks (NEON) Airborne Observation Platform (AOP) imaging spectroscopy and lidar surveys in Gunnison County, Colorado. Folair samples were collected across the East River, Washington Gulch, Slate River, and Coal Creek watersheds and contain a mixture of vegetation including meadow, shrub, and tree foliar samples. The samples were processed using aqua regia digestion and analyzed for elemental determination on inductively coupled plasma optical emission spectrometry (ICP-OES).The data package includes: (1) raw foliar elemental data files in CSV and PDF formats, (2) quality control certificates in PDF format, and (3) an aggregated CSV file containing all elemental measurements compiled across samples. No specialized software is required to access or use these files.

2018 National Ecological Observatory Network Campa↗

Investigating the Impact of Irrigation on Malaria Vector Larval Habitats and Transmission Using a Hydrology‐Based Model

Abstract A combination of accelerated population growth and severe droughts has created pressure on food security and driven the development of irrigation schemes across sub‐Saharan Africa. Irrigation has been associated with increased malaria risk, but risk prediction remains difficult due to the heterogeneity of irrigation and the environment. While investigating transmission dynamics is helpful, malaria models cannot be applied directly in irrigated regions as they typically rely only on rainfall as a source of water to quantify larval habitats. By coupling a hydrologic model with an agent‐based malaria model for a sugarcane plantation site in Arjo, Ethiopia, we demonstrated how incorporating hydrologic processes to estimate larval habitats can affect malaria transmission. Using the coupled model, we then examined the impact of an existing irrigation scheme on malaria transmission dynamics. The inclusion of hydrologic processes increased the variability of larval habitat area by around two‐fold and resulted in reduction in malaria transmission by 60%. In addition, irrigation increased all habitat types in the dry season by up to 7.4 times. It converted temporary and semi‐permanent habitats to permanent habitats during the rainy season, which grew by about 24%. Consequently, malaria transmission was sustained all‐year round and intensified during the main transmission season, with the peak shifted forward by around 1 month. Lastly, we evaluated the spatiotemporal distribution of adult vectors under the effect of irrigation by resolving habitat heterogeneity. These findings could help larval source management by identifying transmission hotspots and prioritizing resources for malaria elimination planning.

60 APPLIED LIFE SCIENCES↗

On the transferability of residence time distributions in two 10-km long river sections with similar hydromorphic units

Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface and their residence time distributions (RTDs) in the subsurface are important for managing the water quality and ecosystem health in dynamic river corridors. However, direct simulating high-spatial resolution HEFs and RTDs can be time-consuming, especially for watershed-scale modeling. Efficient surrogate models linking RTDs to hydromorphic units (HUs) can be alternatives for simulating RTDs in large-scale models. A common concern of these surrogate models, though, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this issue, this work evaluates the HEFs and resulting RTD-HU relationships for two 10-km long river corridors along the Columbia River leveraging a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework we previously developed. Applying such a framework at the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. Finally, our comparison shows that the similarity and transferability of the RTD-HU relationship is very low for the two investigated river sections, which suggests that devising a general algorithm to estimate RTDs based solely on surface water hydrodynamics and short-distance river channel topography data, as well as HU classification, might be nearly impossible.

54 ENVIRONMENTAL SCIENCES↗

The Kimberlina synthetic multiphysics dataset for CO 2 monitoring investigations

Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

58 GEOSCIENCES↗

Data and scripts associated with a manuscript on residence time distribution simulation in two 10-kilometer long river sections

This data package is associated with the publication “On the Transferability of Residence Time Distributions in Two 10-km Long River Sections with Similar Hydromorphic Units” submitted to the Journal of Hydrology (Bao et al. 2024).Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface, along with their residence time distributions (RTDs) in the subsurface, is crucial for managing water quality and ecosystem health in dynamic river corridors. However, directly simulating high-spatial resolution HEFs and RTDs can be a time-consuming process, particularly for watershed-scale modeling. Efficient surrogate models that link RTDs to hydromorphic units (HUs) may serve as alternatives for simulating RTDs in large-scale models. One common concern with these surrogate models, however, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this, we evaluated the HEFs and the resulting RTD-HU relationships for two 10-kilometer-long river corridors along the Columbia River, using a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework that we previously developed. Applying this framework to the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. This data package includes the model inputs files and the simulation results data. This data package contains 10 folders. The modeling simulation results data are in the folders 100H_pt_data and 300area_pt_data, for the study domain Hanford 100H and 300 area respectively. The remaining eight folders contain the scripts and data to generate the manuscript figures. The file-level metadata file (Bao_2024_Residence_Time_Distribution _flmd.csv) includes a list of all files contained in this data package and descriptions for each. The data dictionary file (Bao_2024_Residence_Time_Distribution _dd.csv) includes column header definitions and units of all tabular files.

54 ENVIRONMENTAL SCIENCES↗

Mitigation of spatial nonstationarity with vision transformers

Spatial nonstationarity, the location variance of features’ statistical distributions, is ubiquitous in many natural settings. For example, in geological reservoirs rock matrix porosity varies vertically due to geomechanical compaction trends, in mineral deposits grades vary due to sedimentation and concentration processes, in hydrology rainfall varies due to the atmosphere and topography interactions, and in metallurgy crystalline structures vary due to differential cooling. Conventional geostatistical modeling workflows rely on the assumption of stationarity to be able to model spatial features for geostatistical inference. Nevertheless, this is often not a realistic assumption when dealing with nonstationary spatial data and this has motivated a variety of nonstationary spatial modeling workflows such as trend and residual decomposition, cosimulation with secondary features, and spatial segmentation and independent modeling over stationary subdomains. The advent of deep learning technologies has enabled new workflows for modeling spatial relationships. However, there is a paucity of demonstrated best practice and general guidance on mitigation of spatial nonstationarity with deep learning in the geospatial context. We demonstrate the impact of two common types of geostatistical spatial nonstationarity on deep learning model prediction performance and propose the mitigation of such impacts using self-attention (vision transformer) models. We demonstrate the utility of vision transformers for the mitigation of nonstationarity with relative errors as low as 10%, exceeding the performance of alternative deep learning methods such as convolutional neural networks. We establish best practice by demonstrating the ability of self-attention networks for modeling large-scale spatial relationships in the presence of commonly observed geospatial nonstationarity.

58 GEOSCIENCES↗

High–Resolution Maps of Near–Surface Permafrost for Three Watersheds on the Seward Peninsula, Alaska Derived From Machine Learning

Permafrost soils are a critical component of the global carbon cycle and are locally important because they regulate the hydrologic flux from uplands to rivers. Furthermore, degradation of permafrost soils causes land surface subsidence, damaging infrastructure that is crucial for local communities. Regional and hemispherical maps of permafrost are too coarse to resolve distributions at a scale relevant to assessments of infrastructure stability or to illuminate geomorphic impacts of permafrost thaw. Here we train machine learning models to generate meter–scale maps of near–surface permafrost for three watersheds in the discontinuous permafrost region. The models were trained using ground truth determinations of near–surface permafrost presence from measurements of soil temperature and electrical resistivity. We trained three classifiers: extremely randomized trees (ERTr), support vector machines (SVM), and an artificial neural network (ANN). Model uncertainty was determined using k–fold cross validation, and the modeled extents of near–surface permafrost were compared to the observed extents at each site. At–a–site near–surface permafrost distributions predicted by the ERTr produced the highest accuracy (70%–90%). However, the transferability of the ERTr to the sites outside of the training data set was poor, with accuracies ranging from 50% to 77%. The SVM and ANN models had lower accuracies for at–a–site prediction (70%–83%), yet they had greater accuracy when transferred to the non–training site (62%–78%). These models demonstrate the potential for integrating high–resolution spatial data and machine learning models to develop maps of near–surface permafrost extent at resolutions fine enough to assess infrastructure vulnerability and landscape morphology influenced by permafrost thaw.

54 ENVIRONMENTAL SCIENCES↗

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗

Deep Direct-Use Geothermal Feasibility Study for Hawthorne NV

The Hawthorne Nevada, deep direct-use geothermal study is a two-year effort funded by the U.S. Department of Energy to determine the techno-economic feasibility of implementing a large-scale, direct-use facility for the Hawthorne Army Depot (HAD) and the public facilities of the city of Hawthorne and Mineral County. The approach links a production side analysis (PSA) and a demand side analysis (DSA) into a whole-system analysis (WSA) to provide an integrated assessment of the resource and the probability of delivering economically viable direct-use energy to Hawthorne. Hawthorne, Nevada is in the western part of the Basin and Range province and has been the focus of geothermal investigations for over 40 years. Over the last 15 years, several studies completed by the U.S. Navy Geothermal Program Office (GPO) in conjunction with industry professionals quantified the existence of several low temperature geothermal prospects, the most promising of which is called Prospect A. The promise of Prospect A is based on drilling and flow testing that produced ~100 °C water at flow rates of up to 31 l/s (500 gallons per minute). Measured productivity indexes range from 40-85 l/s/MPa, suggesting a warm and productive heat source. Despite the promise of the resource, uncertainties in its spatial extent and long-term sustainability mean that techno-economic analyses must include probabilities of the sustainability of the resource under different operating scenarios. Here, the PSA is conducted by integrating a wide range of disparate data to estimate lognormal P90, P50, and P10 resource capacities. These capacities are used as input to a thermal-hydrologic (T-H) model to estimate thermal drawdown for each capacity estimate for several different DSA scenarios. Using a systems-based approach, the WSA links the dynamic T-H simulations of the PSA/DSA combinations with the techno-economic model GEOPHIRES to account for both the temporal dynamics and uncertainties in the system to produce probabilistic distributions of several performance metrics including the levelized cost of heat (LCOH) and the return on investment (ROI). This report is the final delivery for the project and documents the study's activities and results.

15 GEOTHERMAL ENERGY↗