Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydrological design”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Representing the Unrepresented Impact of River Ice on Hydrology, Biogeochemistry, Vegetation, and Geomorphology: A Hybrid Physics-Machine Learning Approach

Focal areas include: Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system composed of a hierarchy of models (e.g., AI driven model/component/parameterization selection). Insight gleaned from complex data (both observed and simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Baseflow Identification via Explainable AI With Kolmogorov‐Arnold Networks

Abstract Hydrological models often involve constitutive laws that may not be optimal in every application. We propose to replace such laws with the Kolmogorov‐Arnold networks (KANs), a class of neural networks designed to identify symbolic expressions. We demonstrate KAN's potential on the problem of baseflow identification, a notoriously challenging task plagued by significant uncertainty. KAN‐derived functional dependencies of the baseflow components on the aridity index outperform their original counterparts; they demonstrate that water availability, rather than potential evapotranspiration, drives baseflow by constraining actual evapotranspiration under arid conditions. On a test set, they increase the Nash‐Sutcliffe efficiency (NSE) by 65%, decrease the root mean squared error by 29%, and increase the Kling‐Gupta efficiency by 34%. This superior performance is achieved while reducing the number of fitting parameters from three to two. Next, we use data from 378 catchments across the continental United States to refine the water‐balance equation at the mean‐annual scale. The KAN‐derived equations based on the refined water balance outperform both the current aridity index model, with up to a 105% increase in NSE, and the KAN‐derived equations based on the original water balance. While the performance of our model and tree‐based machine learning methods is similar, KANs offer the advantage of simplicity and transparency and require no specific software or computational tools. This case study focuses on the aridity index formulation, but the approach is flexible and transferable to other hydrological processes. Plain Language Summary Equations used in hydrologic model are often suboptimal, resulting in reduced prediction accuracy and efficiency. We implemented Kolmogorov‐Arnold networks (KAN), a machine learning algorithm for deriving symbolic formulations, to estimate groundwater recharge and showed that it outperforms an existing state‐of‐the‐art semi‐empirical formulation. In hydrology, Nash‐Sutcliffe efficiency (NSE), root mean squared error (RMSE), and Kling‐Gupta efficiency (KGE) are commonly used to evaluate model performance. Higher NSE and KGE values indicate better performance, while lower RMSE values are preferable. Our results show that NSE increased by 71%, RMSE decreased by 32%, and KGE improved by 25%. In addition, KAN identifies an optimal functional form and can be used to derive new analytical formulas using the prior knowledge. The KAN‐inspired equation outperformed the original formulation and reduced the fitting parameters. Furthermore, we refined the water‐balance equation at the mean‐annual scale and showed that, based on the new water‐balance equation, KAN can derive new formulations that are superior to the original aridity index formulations (up to 105% increase in NSE) and KAN‐derived equations based on the original water balance. These findings highlight the significant potential of KAN to advance the scientific understanding of a wide range of hydrologic processes. Key Points Kolmogorov‐Arnold networks (KANs) enhance interpretability of machine‐learned hydrological models KAN‐derived symbolic formulations outperform state‐of‐the‐art semi‐empirical aridity indices KAN‐identified functional form yields an analytical index with fewer fitting parameters and improved performance

baseflow↗
DOE OSTI2020

Y-12 Groundwater Protection Program Monitoring Well Inspection and Maintenance Plan

This plan describes the systematic approach for: inspecting the physical condition of monitoring wells at Y-12, determining maintenance needs that extend the life of a well, and identifying those wells that no longer meet acceptable monitoring well design or well construction standards and require plugging and abandonment. The inspection and maintenance of groundwater monitoring wells is one of the primary management strategies of the Y-12 Groundwater Protection Program (GWPP) Management Plan, that is, the “proactive stewardship of the extensive monitoring well network at Y-12" (Consolidated Nuclear Security, L.L.C. [CNS], 2018). Effective stewardship, and a program of routine inspections of the physical condition of each monitoring well, ensures that representative water-quality samples and hydrologic data are obtained from the well network and protects the subsurface environment. In accordance with the Y-12 GWPP Monitoring Optimization Plan (MOP) for Groundwater Monitoring Wells at the Y-12 National Security Complex, Oak Ridge, Tennessee (CNS, 2017), the status designation (active or inactive) for each well determines the scope and extent of well inspections and maintenance activities. This plan, in conjunction with the above document, formalizes the GWPP approach to focus available resources on monitoring wells which provide the most useful data, and for that reason the GWPP inspects and performs maintenance on the wells sampled by the GWPP. This plan applies to groundwater monitoring wells installed at Y-12 and the related waste management facilities located within the three hydrogeologic regimes: (1) the Bear Creek Hydrogeologic Regime (Bear Creek Regime), (2) the Upper East Fork Poplar Creek Hydrogeologic Regime (East Fork Regime), and (3) the Chestnut Ridge Hydrogeologic Regime (Chestnut Ridge Regime). The Bear Creek Regime encompasses the section of Bear Creek Valley (BCV) immediately west of Y-12. The East Fork Regime encompasses most of the Y-12 process, operations, and support facilities in BCV west of Scarboro Road. The Chestnut Ridge Regime is directly south of Y-12 and encompasses a section of Chestnut Ridge that is bounded to the west by a surface drainage feature (Dunaway Branch, located immediately west of Industrial Landfill II) and by Scarboro Road to the east. The GWPP maintains an extensive database of geographic and construction details and related information for the monitoring wells in each hydrogeologic regime in the Updated Subsurface Database for Bear Creek Valley, Chestnut Ridge, and Parts of Bethel Valley on the U.S. DOE Oak Ridge Reservation (CNS, 2019). A detailed description of the hydrogeologic framework at Y-12 can be found in the GWPP Management Plan (CNS, 2018).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Monitoring Potential Transport of Radioactive Contaminants in Shallow Ephemeral Channels: FY2019

Desert Research Institute (DRI) conducted a field assessment of the potential for contaminated soil to be transported from the Smoky Site Contamination Area (CA) because of storm runoff. This activity supported U.S. Department of Energy (DOE) Environmental Management Nevada Program (EM NV) efforts to establish post-closure monitoring plans for the Smoky Site Soils Corrective Action Unit (CAU) 550. The work was intended to confirm the likely mechanism of transport and determine the meteorological conditions that might cause the movement of contaminated soils, as well as determine the particle size fraction most closely associated with transported radionuclide-contaminated soils. These data will facilitate the design of the appropriate post-closure monitoring program. In 2011, DRI installed a meteorological monitoring station on the west side of the Smoky Site CA and a hydrologic (runoff) monitoring station within the CA, near the east side. The meteorological station collected air temperature, wind speed, wind direction, relative humidity, precipitation, solar radiation, barometric pressure, soil temperature, and soil water content data. The maximum, minimum, and average or total values (as appropriate) for each of these parameters were recorded for each 10-minute interval. The maximum, minimum, and average water depth in the flume installed at the hydrologic station were also recorded for each 10-minute interval. This report presents the data collected from these stations during fiscal year (FY) 2019. During the FY2019 reporting period, the warmest months were June, July, and August and the coldest were December and February. The highest monthly solar radiation (i.e., sunshine) values were recorded in June, July, and August, whereas December through February had the lowest solar radiation values. Monthly mean wind speeds were highest in the spring (April) and early summer (June). Winds were predominantly from the north, except in February, July, August, and September when winds were mostly from the west. The monthly average relative humidity ranged from 17 percent to approximately 60 percent. Humidity was lowest in the summer and highest in the winter. Monthly total precipitation ranged from 0.0 during August and September to 2.22 inches (in) (56.39 millimeters [mm]) in February. Total precipitation for FY2019 was 7.83 in (198.88 mm). During the reporting period, a single, nonzero, flow event was recorded at the flume. This occurred over two days between March 6 and 7, 2019, in response to 1.16 in (29.46 mm) of precipitation falling within 9.3 hours on March 6, 2019. This event produced a peak water depth of 0.39 in (0.99 centimeters [cm]) in the flume, and runoff was measured for approximately 23.5 hours. However, the maximum recorded flow depth was less than the limit for flow discharge and velocity that can be calculated for the six-inch (15 cm) throat of a Parshall flume. The runoff event produced insufficient flow velocities through the natural channel to cause local erosion and transport of bedload materials. No bedload samples were collected for radiological analysis. Observing meteorological and environmental conditions that lead to storm runoff will help identify the parameters and threshold conditions that should be incorporated into post-closure monitoring plans for this and similar Soils CAU sites. Routine monitoring of meteorological and hydrologic parameters at the Smoky Site CA has been discontinued at the request of EM NV, effective at the conclusion of FY2019. Environmental data acquisition was completed and the meteorological and runoff monitoring stations were removed from the Smoky Site on October 15, 2019.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Sensitivities of subgrid-scale physics schemes, meteorological forcing, and topographic radiation in atmosphere-through-bedrock integrated process models: a case study in the Upper Colorado River basin

Abstract. Mountain hydrology is controlled by interacting processes extending from the atmosphere through the bedrock. Integrated process models (IPMs), one of the main tools needed to interpret observations and refine conceptual models of the mountainous water cycle, require meteorological forcing that simulates the atmospheric process to predict hydroclimate then subsequently impacts surface–subsurface hydrology. Complex terrain and extreme spatial heterogeneity in mountainous environments drive uncertainty in several key considerations in IPM configurations and require further quantification and sensitivity analyses. Here, we present an IPM using the Weather Research and Forecasting (WRF) model which forces an integrated hydrologic model, ParFlow-CLM, implemented over a domain centered over the East River watershed (ERW), located in the Upper Colorado River basin (UCRB). The ERW is a heavily instrumented 300 km2 region in the headwaters of the UCRB near Crested Butte, CO, with a growing atmosphere-through-bedrock observation network. Through a series of experiments in the water year 2019 (WY19), we use four meteorological forcings derived from commonly used reanalysis datasets, three subgrid-scale physics scheme configurations in WRF, and two terrain shading options within WRF to test the relative importance of these experimental design choices for key hydrometeorological metrics including precipitation and snowpack, as well as evapotranspiration, groundwater storage, and discharge simulated by the ParFlow-CLM. Our hypothesis is that uncertainty from synoptic-scale forcings produces a much larger spread in surface–subsurface hydrologic fields than subgrid-scale physics scheme choice. Results reveal that the WRF subgrid-scale physics configuration leads to larger spatiotemporal variance in simulated hydrometeorological conditions, whereas variance across meteorological forcing with common subgrid-scale physics configurations is more spatiotemporally constrained. Despite reasonably simulating precipitation, a delay in simulated discharge peak is due to a systematic cold bias across WRF simulations, suggesting the need for bias correction. Discharge shows greater variance in response to the WRF simulations across subgrid-scale physics schemes (26 %) rather than meteorological forcing (6 %). The topographic radiation option has minor effects on the watershed-average hydrometeorological processes but adds profound spatial heterogeneity to local energy budgets (±30 W m−2 in shortwave radiation and 1 K air temperature differences in late summer). This is the first presentation of sensitivity analyses that provide support to help guide the scientific community to develop observational constraints on atmosphere-through-bedrock processes and their interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

rabpro: global watershed boundaries, river elevation profiles, and catchment statistics

River and Basin Profiler (rabpro) is a Python package to delineate watersheds, extract river flowlines and elevation profiles, and compute watershed statistics for any location on the Earth’s surface. As fundamental hydrologically-relevant units of surface area, watersheds are areas of land that drain via aboveground pathways to the same location, or outlet. Delineations of watershed boundaries are typically performed on digital elevation models (DEMs) that represent surface elevations as gridded rasters. Depending on the resolution of the DEM and the size of the watershed, delineation may be very computationally expensive. With this in mind, we designed rabpro to provide user-friendly workflows to manage the complexity and computational expense of watershed calculations given an arbitrary coordinate pair. In addition to basic watershed delineation, rabpro will extract the elevation profile for a watershed’s mainchannel flowline. This enables the computation of river slope, which is a critical parameter in many hydrologic and geomorphologic models. Finally, rabpro provides a user-friendly wrapper around Google Earth Engine’s (GEE) Python API to enable cloud-computing of zonal watershed statistics and/or time-varying forcing data from hundreds of available datasets. Altogether, rabpro provides the ability to automate or semi-automate complex watershed analysis workflows across broad spatial extents.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Integrated modeling for assessing climate change impacts on water resources and hydropower potential in the Himalayas

Regional hydroclimatic variability and change can affect water resources and hydropower generation. It is essential to assess hydropower potential under current and future climatic conditions to inform the design and operation of hydropower infrastructures. Here, we employ an integrated modeling framework to assess the impact of projected hydroclimatic conditions on water resource systems and hydropower generation. The integrated framework samples climate model outputs under different scenarios to force a hydrologic model and produces streamflow projections. The projected streamflows are inputs for the future hydropower potential assessment. We implement the framework in the central Himalayan river basin. Our results demonstrate substantial spatiotemporal variability in different water balance components (precipitation, evapotranspiration, and water yield) under current and future climatic conditions. For the Himalayan Tila river basin, the annual average energy production is expected to increase under future hydroclimatic conditions (up to 39% in Tila-2 hydropower project, suggested by ensemble mean). Furthermore, this increase in energy is driven mainly by the increased streamflow projections, particularly during the dry season and in the late century. Our results highlight the impacts of hydroclimatic variability in hydropower productions and are of practical use to provide decision-relevant information for designing and operating hydropower infrastructures. The integrated modeling framework presented here is region-specific; however, the approach is reproducible, and the overall insights are generalizable across the Himalayan region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Utilizing Temperature and Brine Inflow Measurements to Constrain Reservoir Parameters During a Salt Heater Test

Brine availability in salt has multiple implications for the safety and design of a nuclear waste storage facility. Brine availability includes both the distribution and transport of brine through a damaged zone around boreholes or drifts excavated into the salt. Coupled thermal, hydrological, mechanical, and chemical processes taking place within heated bedded salt are complex; as part of DECOVALEX 2023 Task E this study takes a parsimonious modeling approach utilizing analytical and numerical one-dimensional simulations to match field measurements of temperature and brine inflow around a heater. The one-dimensional modeling results presented arrive at best-fit thermal conductivity of intact salt, and the permeability and porosity of damaged salt of 5.74 W/m∙K, 10 –17 m 2 , and ≈ 0.02, respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

WHOLESCALE: Seismic Survey Metadata from San Emidio Nevada 2021

This is a collection of metadata from the 2021 seismic survey at the San Emidio geothermal field in Nevada. In April and May 2021, 37 tri-axial short period seismographs were deployed in a 1.8km diameter cluster centered on 40.367278 deg N, 119.409019 deg W. The first data record started at 2021-04-06T07:09:10Z UTC and the last record ended 2021-05-11T02:58:52Z UTC. The 37 stations included 29 SmartSolo IGU-16HR 3C all-in-one 5 Hz seismographs and 8 DataCube seismographs with 4.5 Hz HGS HG-6(B coil) tri-axial geophones. The raw format (level 0) data includes 353 GB of 500 sps SmartSolo data in native DLD format, 113 GB of 400 sps DataCube data in native DataCube format, and 3.4 GB of GPS data collected during the RTK GPS survey. The SAC data (level 1) files include 564 GB of hourly SAC files. The experiment was designed to monitor seismic activity before, during, and after the planned three-day plant maintenance shutdown April 19-21, 2021. The pumping stopped at 2021/04/19 12:51:45 UTC and resumed about 2021/04/21 21:00:00 UTC. The raw and processed data are in an associated GDR submission, linked below. The metadata here includes files containing experiment details, station locations, seismic data logger specifications, instrumentation specifications, and descriptions of data. Also included are data and metadata from a 2016 seismic survey at the same site.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Model Inputs, Outputs, and Scripts associated with: “Spatial microbial respiration variations in the hyporheic zones within the Columbia River Basin”

This data package is associated with the publication “Spatial microbial respiration variations in the hyporheic zones within the Columbia River Basin” published in the Journal of Geophysical Research: Biogeosciences (Son et al. 2022) available at doi: 10.1029/2021JG006654. This data package includes the key model inputs/outputs of the river corridor model for the Columbia River Basin (CRB) and the model source codes, which were used in the manuscript. The model is a carbon-nitrogen-coupled river corridor model (RCM), and the model is used to quantify hyporheic zone (HZ) aerobic and anaerobic respiration at the NHDPLUS stream reach scales. The RCM used in this study combines empirical substrate models derived from observations and three microbially driven reactions to compute respiration of the HZ for each National Hydrography Dataset (NHD) reach within the CRB. The reactions in HZs of each NHD reach include anaerobic respiration and two-step anaerobic respiration via denitrification. Our HZ respiration estimates are limited to the lotic (or flowing) stream/river systems, and do not account for the respiration process in water column. Note that the RCM only simulates the HZ’s contribution to the dissolved carbon dioxide (CO2) concentrations in the streams, and the CO2 emissions to the atmosphere are not modelled. The model computes at hourly timesteps because of the fast reaction rates. The key input data of the model are exchange flux, residence time, and stream solute (dissolved organic carbon (DOC), dissolved oxygen (DO), and nitrate concentrations). These inputs are constant over time and represent long-term averaged values.This modeling framework successfully quantified HZ respiration components over multiple scales. It revealed key mechanisms driving the spatial variation of HZ aerobic and anaerobic respiration in reaches with varying hydrologic and substrate conditions. Thus, this modeling study offers a testing hypothesis in different river system (e.g., climate and biomes) for the HZ respiration processes, and can be used as a sampling design tool for large-scale HZ experimental studies.This dataset contains five folders: (1) model_inputs, (2) model_outputs, (3) Rscripts, (4) figures, and (5) model_codes. It also contains a readme, file level metadata (FLMD), and data dictionary (dd). Please see the FLMD for a list of all the files contained in this data package and descriptions for each. The model_inputs folder contains the model inputs used to drive the model simulations. The model_outputs folder contains key model output files from the river corridor model. The Rscripts folder contains the Rscripts for pre- and post- processing model results. The figures folder contains the raw figures associated with the manuscript. The model_codes folder includes key model source codes/input files. All files are .jpg, .jpeg, .out, .e, .od, .dat, .sub, .F90, .0, .R, .sbx, .cpg, .sbn, .shx, .shp, .dbf, .prj, .tfw, .tif, .xml, .pdf, or .csv.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Influence of local river hydraulics on Kootenai River white sturgeon ( Acipenser transmontanus ) habitat selection during four spawning years, 2017–2020

Understanding fine-scale habitat selection of endangered Kootenai River white sturgeon ( Acipenser transmontanus) is an important component for monitoring and recovery efforts. Fine-scale habitat selection and quantifying temporal changes in suitable habitat contributes to the work of addressing recruitment failure within the Kootenai River population. Habitat suitability indices were developed using over 96 000 acoustic telemetry sturgeon detections and two-dimensional hydrodynamic model simulations near Bonners Ferry, Idaho, USA. The selected habitat was assessed to develop habitat suitability indices for sturgeon; females undergoing spawn migrations and non-spawners. The most frequented locations were 8–9 m deep and water velocities of 0.3–0.7 m·s −1 . These observations suggest sturgeon with different spawning capabilities selected similar habitat. Weighted usable area was calculated to understand temporal variability in habitat quality, which showed a positive relationship with increases in flow. Results help understand the habitat limiting factors in regulated hydrologic regimes; provide biologists insight for monitoring efforts in discrete habitat conditions; guidance for water managers and the regulation of upstream water resources; and guidance to restoration practitioners for in-stream structure designs.

Dudunake, Taylor J. (ORCID:0000000176502419)↗
DOE OSTI2023

Designing Iowa Agricultural Landscapes to Improve Environmental Co-Benefits of Bioenergy Production

Cellulosic bioenergy feedstocks are needed to improve carbon (C) management while provisioning biomass for bioproducts and biofuel. The transition to increased cellulosic biomass production can be guided by land management plans designed to improve economic, environmental, and ecological performance. We constructed a sustainability model to compare landscape designs for biofuel production from corn (Zea mays L.) stover and switchgrass (Panicum virgatum L.) in central Iowa, USA. We used the model to compare environmental and socioeconomic outcomes associated with four landscape management strategies, with and without cellulosic biomass markets. We evaluated (1) a fuelshed area containing over 1.2 million ha (3 million acres) of corn and soybean (Glycine max (L.) Merr.) within 80 km (50 miles) of a commercial-scale cellulosic biorefinery in Nevada, Iowa, and (2) the South Fork watershed containing over 72,000 ha (178,000 acres) of these row crops within eight north central Iowa HUC-12 (hydrologic unit code) watersheds. At both landscape scales, we found that it is possible to achieve multiple environmental and socioeconomic benefits concomitantly with cellulosic biomass production by strategically collecting corn stover and converting the 10% of the lowest-profitability row crop land to perennial switchgrass. Potential benefits from landscape design include increased biodiversity, soil and water quality improvements, increased soil carbon sequestration for climate change mitigation, and reduced fertilizer use and cost. Our model results showed that increasing benefits can accrue when complementary conservation practices (e.g., reduced tillage, use of a rye cover crop) are combined and integrated throughout a fuelshed or watershed area. We conclude that ecologically based landscape designs offer valuable insights about costs and benefits of land management alternatives, with relevance for achieving stakeholder goals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Energy Exascale Earth System Model v2.1.0

First release of version 2.1 of the Energy Exascale Earth System Model. [ATM] The atmosphere component remains EAM. There are no major changes in the default configuration since 2.0. New features include: A semi-lagrangian tracer transport for theta-l dycore, a new algorithm for finding the tropopause, new RRM mesh configurations. Add and update SSP370 and SSP585 cases. Restore the FIDEAL case. [LAND] The land component is ELM. There are no major changes in the default configuration since 2.0. Several option features have been added including: implementation of topography-based subgrid structure (topounits) and accompanying parameterizations and atmospheric forcing downscaling methods; a new plant hydraulics scheme; two-way land-river hydrological coupling through the infiltration of floodplain water; an implementation of perennial crops; updates to the SNICAR-AD snow radiative transfer model; and implementation of soil erosion and sediment yield in ELM-Erosion. Each of these new changes is modular in design and can be turned on or off as the user specifies; they are currently being tested in different “BGC” configurations. [OCEAN] The ocean component remains MPAS-Ocean. Major change since version 2.0 include the addition of the Fox-Kemper et al. 2011 parameterization for submesocale eddies, a correction for barotropic thickness consistency that reduces divergence noise, and the addition of an ocean carbon conservation analysis member. [SEAICE] The sea-ice component remains MPAS-Seaice. Major changes since version 2.0 include: A correction to how shortwave parameters are interpolated in the snicar-ad 5-band radiation scheme, the addition of a sea ice carbon conservation analysis member, updates to the default sea ice biogeochemistry namelist parameters to be consistent with version 2.0 improvements to nitrogen cycling and a correction in the ice-ocean dissolved organic nitrogen coupling. [LAND ICE] The land-ice component remains MPAS-Albany-landIce (MALI). Major changes since 2.0 include an update to the MALI version and the Greenland mesh.[RIVER] The river model is MOSART. There are no major changes in the default configuration since 2.0. A major new optional feature is two-way river-ocean hydrological coupling between MOSART and MPAS-O. This change can be turned on or off as the user specifies, and is being tested in different configurations. [COUPLER] The coupler remains cpl7/MCT. Major changes since version 2.0 include: Carbon budget calculated when heat/water budgets active. Fix a bug in land-atm fluxes for tri-grid configurations. [OTHER] a small bug in the zenith angle calculation was fixed in the data models.

ECP↗
DOE OSTI2020

A Model Tree Generator (MTG) Framework for Simulating Hydrologic Systems: Application to Reservoir Routing

Data-driven algorithms have been widely used as effective tools to mimic hydrologic systems. Unlike black-box models, decision tree algorithms offer transparent representations of systems and reveal useful information about the underlying process. A popular class of decision tree models is model tree (MT), which is designed for predicting continuous variables. Most MT algorithms employ an exhaustive search mechanism and a pre-defined splitting criterion to generate a piecewise linear model. However, this approach is computationally intensive, and the selection of the splitting criterion can significantly affect the performance of the generated model. These drawbacks can limit the application of MTs to large datasets. To overcome these shortcomings, a new flexible Model Tree Generator (MTG) framework is introduced here. MTG is equipped with several modules to provide a flexible, efficient, and effective tool for generating MTs. The application of the algorithm is demonstrated through simulation of controlled discharge from several reservoirs across the Contiguous United States (CONUS).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

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↗
DOE OSTI2020

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in deep clay/shale/argillaceous rock. International collaboration activities such as heater tests and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests has as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum vs. discrete) to tackle issues related to flow and transport at various scales of the host-rock and EBS design concept. Consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the EBS design. This is particularly important for backfilled repository concepts where temperature plays a key role in the EBS behavior and long-term performance. This report describes multiple R&D efforts on disposal in argillaceous geologic media through development and application of coupled THMC process models, experimental studies on clay/metal/cement barrier and host-rock (argillite) material interactions, molecular dynamic (MD) simulations of water transport during (swelling) clay dehydration, first-principles studies of metaschoepite (UO 2 corrosion product) stability, and advances in thermodynamic plus surface complexation database development. Drift-scale URL experiments provides key data for testing hydrological-chemical (HC) model involving strong couplings of fluid mixing and barrier material chemical interactions. The THM modeling focuses on heater test experiments in argillite rock and gas migration in bentonite as part of international collaboration activities at underground research laboratories (URLs). In addition, field testing at an URL involves in situ analysis of fault slip behavior and fault permeability. Pore-scale modeling of gas bubble migration is also being investigated within the gas migration modeling effort. Interaction experiments on bentonite samples from heater test under ambient and elevated temperatures permit the evaluation of ion exchange, phase stability, and mineral transformation changes that could impact clay swelling. Advances in the development, testing, and implementation of a spent nuclear fuel (SNF) degradation model coupled with canister corrosion focus on the effects of hydrogen gas generation and its integration with Geologic Disposal Safety Assessment (GDSA). GDSA integration activities includes evaluation of groundwater chemistries in shale formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Informing climate adaptation strategies using ecological simulation models and spatial decision support tools

Introduction: Forest landscapes offer resources and ecosystem services that are vital to the social, economic, and cultural well-being of human communities, but managing for these provisions can require socially and ecologically relevant trade-offs. We designed a spatial decision support model to reveal trade-offs and synergies between ecosystem services in a large eastern Cascade Mountain landscape in Washington State, USA. Methods: We used process-based forest landscape (LANDIS-II) and hydrology (DHSVM) models to compare outcomes associated with 100 years of simulated forest and wildfire dynamics for two management scenarios, Wildfire only and Wildfire + Treatments. We then examined the strength and spatial distribution of potential treatment effects and trends in a set of resources and ecosystem services over the simulation period. Results: We found that wildfire area burned increased over time, but some impacts could be mitigated by adaptation treatments. Treatment benefits were not limited to treated areas. Interestingly, we observed neighborhood benefits where fire spread and severity were reduced not only in treated patches but in adjacent patches and landscapes as well, creating potential synergies among some resource benefits and services. Ordinations provided further evidence for two main kinds of outcomes. Positive ecological effects of treatments were greatest in upper elevation moist and cold forests, while positive benefits to human communities were aligned with drier, low- and mid-elevation forests closer to the wildland urban interface. Discussion: Our results contribute to improved understanding of synergies and tradeoffs linked to adaptation and restoration efforts in fire-prone forests and can be used to inform management aimed at rebuilding resilient, climate-adapted landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗