SEARCH · Engineering Papers

Results for “river discharge”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A multi-model assessment of global freshwater temperature and thermoelectric power supply under climate change

Water temperature is a key abiotic factor influencing aquatic ecosystem health and the services provided to both nature and humans. Global water temperature models offer possibilities to improve our understanding of water temperature regimes, which is increasingly important against the backdrop of climate change. Yet, existing studies have predominantly relied on a single model, which can lead to an incomplete representation of uncertainty and potential biases, in addition to limited insight into the range of possible future conditions, which ultimately reduces the robustness of climate impact assessments. Here, we provide a comprehensive assessment of surface freshwater temperature changes from various river and lake models for both past conditions and under future scenarios of climate change. Global models consistently simulate that surface water temperatures are currently 0.5 °C–0.8 °C higher than at the turn of the century (i.e. 1981–2000), and that warming will extend and intensify with future global change throughout the 21st century. While the strength of warming is highly sensitive to the different water temperature models, emissions scenarios and global climate models, our multi-model ensemble shows a global average annual water temperature rise of between +1.3 °C and +4.1 °C by the end of the century. To illustrate a potential societal impact of our results, we evaluate how future changes in discharge and water temperature may affect existing thermoelectric power plants, estimating average annual reductions of 1.5%–6% in global usable capacity by the end of the century. However, with river water temperatures projected to exhibit more pronounced seasonal patterns in the future—especially under the more extreme climate change scenarios and during summer months in the Northern Hemisphere—intra-annual reductions in usable capacity can be much more severe. Given the challenges associated with (large-scale) adaptation to control water temperature regimes, strong climate change mitigation is crucial for minimising water temperature rises and its associated negative impacts on humankind and ecosystems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2025

Assessing the Role of Hydrodynamics in Enhancing Height-Above-the-Nearest-Drainage Derived Synthetic Rating Curves: A Comparative Study in the Wu River Basin, Taiwan

The conventional approach to generating synthetic rating curves (SRC) using the Height-Above-the-Nearest-Drainage (HAND) method typically relies on the assumption of uniform flow, such as Manning's equation, to establish stage-discharge ratings. The zero-physics application of the uniform flow equation is insufficient for capturing detailed hydraulic features (e.g., backwater effect) and neglects the hydraulic effects from adjacent channels. This lack of hydrodynamic computation can impact the accuracy and effectiveness of riverine flood risk estimation and management. To reduce this foreseeable error, we introduce the HAND-hd workflow, which integrates sophisticated hydrodynamic computations in the production of HAND-based SRC with hydrodynamic features (SRC hd ). The results indicate that SRC hd demonstrates consistent agreement with both gauge observations and benchmark solutions. Additionally, the comparative analysis suggests that SRC hd provides notable improvements in stage-discharge ratings over conventional HAND-based SRCs, particularly in channels with mild bed gradients, where it reduces water stage prediction errors and percent biases. In steeper channel segments, SRC hd maintains comparable accuracy to conventional methods. The comprehensive evaluation in this study emphasizes the potential discrepancies and inaccuracies associated with the adoption of the uniform flow assumption in the conventional HAND-SRCs and addresses the necessity of including hydrodynamic physics in the application of HAND-based SRC (e.g., inundation map) in channels with mild gradients.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Structuring Nutrient Yields throughout Mississippi/Atchafalaya River Basin Using Machine Learning Approaches

To minimize the eutrophication pressure along the Gulf of Mexico or reduce the size of the hypoxic zone in the Gulf of Mexico, it is important to understand the underlying temporal and spatial variations and correlations in excess nutrient loads, which are strongly associated with the formation of hypoxia. This study’s objective was to reveal and visualize structures in high-dimensional datasets of nutrient yield distributions throughout the Mississippi/Atchafalaya River Basin (MARB). For this purpose, the annual mean nutrient concentrations were collected from thirty-three US Geological Survey (USGS) water stations scattered in the upper and lower MARB from 1996 to 2020. Eight surface water quality indicators were selected to make comparisons among water stations along the MARB over the past two decades. Principal component analysis (PCA) was used to comprehensively evaluate the nutrient yields across thirty-three USGS monitoring stations and identify the major contributing nutrient loads. The results showed that all samples could be analyzed using two main components, which accounted for 81.6% of the total variance. The PCA results showed that yields of orthophosphate (OP), silica (SI), nitrate–nitrites (NO 3 -NO 2 ), and total suspended sediment (TSS) are major contributors to nutrient yields. It also showed that land-planted crops, density of population, domestic and industrial discharges, and precipitation are fundamental causes of excess nutrient loads in MARB. These factors are of great significance for the excess nutrient load management and pollution control of the Mississippi River. It was found that the average nutrient yields were stable within the sub-MARB area, but the large nitrogen yields in the upper MARB and the large phosphorus yields in the lower MARB were of great concern. t-distributed stochastic neighbor embedding (t-SNE) revealed interesting nonlinear and local structures in nutrient yield distributions. Clustering analysis (CA) showed the detailed development of similarities in the nutrient yield distribution. Moreover, PCA, t-SNE, and CA showed consistent clustering results. This study demonstrated that the integration of dimension reduction techniques, PCA, and t-SNE with CA techniques in machine learning are effective tools for the visualization of the structures of the correlations in high-dimensional datasets of nutrient yields and provide a comprehensive understanding of the correlations in the distributions of nutrient loads across the MARB.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Knowledge-guided graph machine learning for spatially distributed prediction of daily discharge and nitrogen export dynamics

Spatially distributed prediction of streamflow and nitrogen export dynamics is essential for precision management of agricultural watersheds. While temporal deep learning models such as Long Short-Term Memory (LSTM) have shown strong performance at basin scales, their ability to generalize spatially is limited by insufficient representation of spatial dependencies and flow paths, particularly under data-scarce conditions. To address this gap, we propose HydroGraphNet, a knowledge-guided graph machine learning framework that integrates process-based knowledge and explicit spatial learning into temporal modeling. This framework incorporates directed graph topology to encode watershed connectivity and upstream inflows, with mass balance constraints to improve physical consistency. To enhance generalization in sparsely monitored regions, HydroGraphNet is pretrained on synthetic data generated by the SWAT+ (Soil and Water Assessment Tool Plus) model. We evaluated HydroGraphNet in the Upper Sangamon River Basin (44 HUC-12 subwatersheds, 2001–2020) against two LSTM baselines: a lumped basin-level model and a distributed variant. When benchmarked on SWAT+ simulations in pretraining, HydroGraphNet improved test NSEs by 8.9% (discharge) and 13.7% (NO₃–N load) in temporal extrapolation, and by 27.1% and 34.7% in spatial extrapolation, relative to the Lumped LSTM baseline. After fine-tuning with USGS monitoring data, the model achieved mean test NSE (KGE) scores of 0.768 (0.861) for discharge and 0.626 (0.664) for NO₃–N load, substantially outperforming baselines. Attribution analysis further highlighted the importance of upstream inflow representation and graph-based spatial learning in capturing cross-subwatershed dependencies. The model also reproduced seasonal hydrological and biogeochemical patterns consistent with known processes, demonstrating its robustness and process fidelity for spatially distributed prediction. Altogether, HydroGraphNet advances the integration of physical knowledge and spatially explicit learning in hydrological modeling, offering a generalizable framework for distributed modeling to support spatially targeted water quality management in data-scarce watersheds.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A Look At Dissolved Organic Carbon In Streams and Its Effect On DGT Hg Collection

Savannah River Site (SRS) and Oak Ridge National Laboratory (ORNL) were both major material production sites for the U.S.'s nuclear program during the 20. century. Both sites discharged mercury (Hg) to streams. Oak Ridge streams have accumulated significant levels of Hg from early 'historical' discharges. Diffusive Gradient in Thin Films (DGT) samplers may be a labor and material saving method to monitor the Hg contamination of streams at both sites. Humic and Fulvic acids are released during the degradation of plant and animal residue. They are part of the Dissolved Organic Matter that colors the black water streams of the southern U.S.. Humic acid and Fulvic acids have a tendency to sequester heavy metals, like Hg, and may interfere with the reactions needed to draw the Hg through the thin films of the DGT samplers. Project Objective: Use a UV-Vis spectrophotometer to estimate the levels of Dissolved Organic Carbon (DOC) in the streams selected to test the effect of DOC on DGT samplers. Conclusion: UV-Vis is a convenient way to estimate the amount of DOC in streams and may be used to provide more context in DGT Hg sampling. efforts. Future Work: Analyze data from SRS and ORNL DGT samples; Compare the results of the DGT samplers with more traditional speciating methods from ORNL; Use a TOC Spectrophotometer to gain a more accurate picture of stream DOC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Scaling hyporheic nitrogen cycling in large river alluvial aquifers

Research Element 1: Denitrification is a difficult flux to measure when fluxes are low and when the ecosystem is difficult to sample. These attributes describe Nyack floodplain, a large gravel-bedded aquifer adjacent to the Middle Fork Flathead River. Given that denitrification is simply respiration with nitrate as a terminal electron acceptor, we reasoned that N 2 production should negatively covary with oxygen respiration so long that oxygen is available. We sampled 16 wells 8 times throughout a year and analyzed these water samples on a membrane-inlet mass spectrometer for N 2 :Ar and O 2 . We fit multilevel models to N 2 :Ar as a function of dissolved oxygen. Within wells, there was no relationship between O 2 and N 2 :Ar. Among wells N 2 :Ar increased as O 2 declined. Based on the slope of this relationship, ~ -0.07, and estimated oxygen respiration in Nyack derived from prior studies, a possible range for denitrification is 0-8 mmol m -2 d -1 corresponding to a turnover rate for nitrate of 0-0.2/d. This estimated turnover rate is in the lower 20% of rates measured in the LINX 2 study, yet shows that even in a low nitrate (0.005 mmol/L) aquifer denitrification can be a substantial sink for dissolved N. Research Element 3 and 4: We developed a new theoretical approach to describe the relationship between water inflow, water outflow, and water age in porous media control volumes under dynamic flow conditions. Our approach had three distinct steps. First we described the application of well-established equations borrowed from Chemical Engineering to describe water age distributions in the hyporheic zone under steady state conditions. We published the associated results in the scientific journal PLOS ONE. Second, we developed a new theory ("flow spans") to describe how predicted steady-state water age distributions will respond to time varying imbalances in hyporheic discharge and recharge. Finally, we developed a proof-of-concept simulation model based on the new theory and used the model to examine the effects of flow regulation on a simulated river with characteristics representative of western U.S. montane rivers with expansive hyporheic zones.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Technology Development for Dry Storage of Aluminum-Clad Spent Nuclear Fuel - 20490

A candidate disposition pathway for the > 13 MTHM of aluminum-clad spent nuclear fuel (ASNF), owned and managed by the U.S. Department of Energy, is the drying and placement of the SNF into sealed-canister dry storage, with the ASNF-in- canisters 'road-ready' for transportation to and final direct disposal in a repository waste package. Technical information gaps in fuel drying, and fuel dry storage behavior, have previously challenged the declaration of technology readiness for drying and placement of this fuel into the DOE Standard Canister design for > 50 years of safe dry storage. The principal technical information gaps included: i) characterization and thermal dehydration behavior of aluminum (oxy)hydroxide films attendant on the cladding due to film formation during reactor operation and during post-discharge up to long-term wet storage histories; and ii) G-value data to enable estimation of the radiolytic gas generation from the cladding with its (oxy)hydroxide films. Thus, the oxide films on the ASNF challenged the safety of a sealed storage canister with thermal and radiolytic decomposition of the waters on the films that can lead to corrosion, pressurization, and flammability issues. These gaps in the technical information base have largely been closed. This paper discusses the investigations at the Idaho National Laboratory (INL) and the Savannah River National laboratory (SRNL), and outlines the pending technology development work for input to an engineering design to enable a road-ready dry storage system for ASNF. The ASNF inventory considered for road-ready dry storage is stored at the Savannah River Site (SRS) and at the Idaho Nuclear Technology and Engineering Center (INTEC) at the INL. The ASNF inventory in the SRS L Basin is from foreign and domestic research reactors (FRR and DRR), and is diverse in terms of design, irradiation, and post-reactor-discharge storage conditions; these factors yield a range of characteristics of cladding oxide films on ASNF. Mixed aluminum (oxy)hydroxide (boehmite and bayerite/gibbsite) films, non-uniform in thicknesses up to a maximum local thickness bounded by 25 μm, were observed on ASNF materials removed from wet storage in the L Basin and in non-sealed dry storage at the SRS following reactor service and a long-term (up to 40+ years) interim storage history. The ASNF inventory at INTEC, in both wet and vented dry storage, is predominantly from the Advanced Test Reactor (ATR), but it also includes DRR and FRR fuel. To address a profound behavior of these films, radiolysis testing of aluminum specimens with mixed type boehmite/bayerite oxide films was performed to develop basic data on G-values for production of radiolytic hydrogen under dry storage conditions with nominal relative humidity, temperature, and cover gas. Modeling and simulation of canister internal environments with postulated inventories of oxide films provides estimation of the evolution of the conditions of the canister loaded with ASNF. Simultaneous Thermal Analysis (TGA/DSC) of hydrated oxide powders, and laboratory-scale aluminum specimens with a bayerite film (∼10 μm) using TGA methods, inform drying time/temperature conditions to be used for the ASNF. A demonstration project is recommended for Verification and Validation of the drying and storage of the ASNF. Remaining major tasks leading up to the hot demonstration include scale-up radiolysis testing and scale-up drying testing. Engineering design with the information from the technology program will establish the safety basis and enable long-term (> 50 years) dry storage compatibility with ASNF in the DOE Standard Canister pending its transportation to and disposal in a repository. This full capability would show the ASNF-in-canister storage system to be road-ready. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Remote Sensing of Multitemporal Functional Lake‐To‐Channel Connectivity and Implications for Water Movement Through the Mackenzie River Delta, Canada

The Mackenzie River Delta in Canada is a mediator of hydrological transport between the expansive Mackenzie River watershed and the Beaufort Sea. Within the delta, lakes frequently act as water and sediment traps, limiting or delaying the movement of material to the coastal ocean. The degree to which this filtering takes place depends on the ease with which sediment-laden water is transported from distributary channels into deltaic lakes, referred to as functional lake-to-channel connectivity, which varies both spatially and temporally. Tracking of connectivity has previously been limited to either small regions of the delta or has focused on a snapshot of connectivity at a single instance in time. Here we describe an algorithm that uses Landsat imagery to track summertime functional lake-to-channel connectivity of 10,362 lakes between 1984 and 2022 on an image-by-image basis. We calculate a total average connected lake area of 1400.7 km 2 during the 2 weeks after peak discharge, 763.6 km 2 higher than previous estimates, suggesting a larger influence of connected lakes on water movement through the delta than previously estimated. We also identify water level thresholds that lead to the initiation of high sediment river water movement into 5,989 lakes (908 lakes with uncertainty ≤±0.5 m), and identify an additional 2899 lakes whose connectivity does not vary at all. As the Arctic hydrological cycle responds to climate change, this work lays a foundation for tracking the movement of water, and the matter it carries, from the Mackenzie River watershed to the Beaufort Sea.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Rheology and Flow Evaluation of Neutralized Sodium Reactor Experiment Fuel with Manganous Nitrate

H-Canyon is preparing the Sodium Reactor Experiment (SRE) solutions in Tanks 16.3 and 16.4 for discard to the Savannah River Site (SRS) High Level Waste (HLW) Tanks into Sludge Batch (SB) 10. To meet HLW's criticality requirements, manganese will be added to the SRE solutions. The addition of manganese to the existing thorium and uranium in the SRE solution raised concerns with the flow of this neutralized material as it is discharged from H-Canyon through the gravity drain system to the H-Area Pump Pit (HPP). A neutralized Mn adjusted SRE stream will produce primarily Mn, Th, and U solids that can deter transfer. A 2012 rheology study was completed as part of the flowsheet development.1 However, the initial rheology study did not include the addition of manganous nitrate as a poison. Therefore, H-Canyon Engineering requested Savannah River National Laboratory (SRNL) to determine if the neutralized, Mn-adjusted SRE solutions will flow through the waste header to the HPP. The H-Canyon Technical Task Request (TTR) specified a target of 80 to 1 to bound the uncertainty in Mn target. Parallel studies were being performed to ensure that freshly precipitated Mn did not have a solubility that would result in challenging the DWPF WAC requirements of 70 to 1. This task was requested via a TTR and is governed by a Task Technical and Quality Assurance Plan (TTQAP).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Quantifying microbial roles in environmental iron oxidation via an integrated kinetics, `omics and metabolic modeling study (Final Report)

Iron oxyhydroxides are extremely reactive components of environmental systems, and therefore exert a strong influence on biogeochemical cycles. These oxyhydroxides strongly adsorb many biologically-relevant elements, including organic carbon and phosphate, as well as a wide range of metals including uranium and actinide species. Thus, the formation mechanism of iron oxyhydroxides are key to understanding both nutrient and contaminant cycling. Microorganisms can catalyze iron oxidation and promote the formation of Fe biominerals and thus are increasingly recognized as important players in biogeochemical cycling. However, it is completely unknown how much of environmental iron oxidation is biologically mediated versus abiotic, and various challenges in studying microbial iron oxidation have hindered accurate incorporation into hydrobiogeochemical models. The overarching goal of our work was to quantify and constrain microbial iron oxidation rates and use ‘omics to gain insight into the controls on this process, while developing tools to enable integration of biotic iron oxidation into hydrobiogeochemical models. Our work focused on the Savannah River Site (SRS) in South Carolina, where extensive microbial iron oxidation has been observed. At Tims Branch, part of the Argonne National Laboratory Wetland Hydrobiogeochemistry Science Focus Area (Argonne SFA), where groundwater discharges into a stream, iron-oxidizing microbial mats form and appear to be a major sink of uranium. In the wetlands that surround Tims Branch, there are wide swaths of iron microbial mats and flocs (mobilized mat). We measured biotic and abiotic iron oxidation rates using mats and water sampled from these sites and found that iron oxidation is primarily carried out by chemolithotrophic microorganisms. The resulting rate constants can be incorporated into models. These mats were characterized by metagenomics and metatranscriptomics, which showed that aerobic chemolithotrophs were the dominant iron-oxidizing bacteria (FeOB), and these included Gallionellaceae and Leptothrix, and possibly Rhodoferax, which is known as an Fe-reducer but may also oxidize Fe(II). This demonstrated that diverse FeOB can coexist and suggests that there are a range of niches and therefore drivers of chemolithotrophic iron oxidation. Analysis of reconstructed genomes strongly suggests that a major factor in diversity is carbon source, as genomes contained varied pathways for autotrophy and heterotrophy. We performed an in-depth analysis of Leptothrix ochracea genomes, since this sheath-former is one of the primary mat builders, yet its physiology remained unresolved. A combination of genomics, transcriptomics, and metabolic modeling suggest that L. ochracea grows mixotrophically using a combination of Fe(II) and organics for energy and both inorganic and organic carbon to create biomass. This contrasts with the largely autotrophic Gallionellaceae (Gallionella, Sideroxydans, and Ferriphaselus) also present in the mats and flocs. Remarkably, multiple FeOB, both Leptothrix and Gallionellaceae, showed activity in response to Fe(II) in live mat incubations. We tracked the gene expression of individual MAGs to Fe(II) and found that various autotrophic and heterotrophic FeOB responded to Fe(II), increasing expression of both carbon fixation and organic utilization genes. The results of the integrated field, kinetics, and omics studies give detailed insight into 1) the taxa that oxidize Fe, and 2) how they connect Fe, C, and N cycles. Towards the goal of connecting omics data to hydrobiogeochemical models, we worked with the KBase team to create a template metabolic model for chemolithotrophic iron oxidation. We initially modeled the well-characterized isolate Gallionellaceae Sideroxydans lithotrophicus, and also applied the model to the mixotroph L. ochracea. In all, we have characterized diverse FeOB in a representative wetland system and solved key problems that enable better incorporation of iron-oxidizing microbes into hydrobiogeochemical models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗

DOE OSTI2022

Hyporheic hydraulic geometry: Conceptualizing relationships among hyporheic exchange, storage, and water age

Hyporheic exchange is now widely acknowledged as a key driver of ecosystem processes in many streams. Yet stream ecologists have been slow to adopt nuanced hydrologic frameworks developed and applied by engineers and hydrologists to describe the relationship between water storage, water age, and water balance in finite hydrosystems such as hyporheic zones. Here, in the context of hyporheic hydrology, we summarize a well-established mathematical framework useful for describing hyporheic hydrology, while also applying the framework heuristically to visualize the relationships between water age, rates of hyporheic exchange, and water volume within hyporheic zones. Building on this heuristic application, we discuss how improved accuracy in the conceptualization of hyporheic exchange can yield a deeper understanding of the role of the hyporheic zone in stream ecosystems. Although the equations presented here have been well-described for decades, our aim is to make the mathematical basis as accessible as possible and to encourage broader understanding among aquatic ecologists of the implications of tailed age distributions commonly observed in water discharged from and stored within hyporheic zones. Our quantitative description of “hyporheic hydraulic geometry,” associated visualizations, and discussion offer a nuanced and realistic understanding of hyporheic hydrology to aid in considering hyporheic exchange in the context of river and stream ecosystem science and management.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Augmentation of WRF-Hydro to simulate overland-flow- and streamflow-generated debris flow susceptibility in burn scars

In steep wildfire-burned terrains, intense rainfall can produce large runoff that can trigger highly destructive debris flows. However, the ability to accurately characterize and forecast debris flow susceptibility in burned terrains using physics-based tools remains limited. Here, we augment the Weather Research and Forecasting Hydrological modeling system (WRF-Hydro) to simulate both overland and channelized flows and assess postfire debris flow susceptibility over a regional domain. We perform hindcast simulations using high-resolution weather-radar-derived precipitation and reanalysis data to drive non-burned baseline and burn scar sensitivity experiments. Our simulations focus on January 2021 when an atmospheric river triggered numerous debris flows within a wildfire burn scar in Big Sur – one of which destroyed California's famous Highway 1. Compared to the baseline, our burn scar simulation yields dramatic increases in total and peak discharge and shorter lags between rainfall onset and peak discharge, consistent with streamflow observations at nearby US Geological Survey (USGS) streamflow gage sites. For the 404 catchments located in the simulated burn scar area, median catchment-area-normalized peak discharge increases by ~ 450 % compared to the baseline. Catchments with anomalously high catchment-area-normalized peak discharge correspond well with post-event field-based and remotely sensed debris flow observations. We suggest that our regional postfire debris flow susceptibility analysis demonstrates WRF-Hydro as a compelling new physics-based tool whose utility could be further extended via coupling to sediment erosion and transport models and/or ensemble-based operational weather forecasts. Given the high-fidelity performance of our augmented version of WRF-Hydro, as well as its potential usage in probabilistic hazard forecasts, we argue for its continued development and application in postfire hydrologic and natural hazard assessments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Hanford Site Composite Analysis: Recharge and Inventory Sensitivity Analyses for the A Trenches, BC Cribs and Trenches, and PUREX Area Vadose Zone Models

The main objectives of the composite analysis (CA) compliance case vadose zone (VZ) facet are to simulate transport of radionuclide releases from the surface to the water table, and to provide radionuclide transfer rates to groundwater for inclusion into the Plateau-to-River (P2R) saturated zone (SZ) model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). The CA compliance case approach to address the scale and distribution of contaminant sources in the VZ facet is to subdivide the Hanford Site Central Plateau into 26 models that contain contaminant sources and liquid discharges likely to commingle during migration through the VZ to the water table

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Train, Inform, Borrow, or Combine? Approaches to Process–Guided Deep Learning for Groundwater–Influenced Stream Temperature Prediction

Although groundwater discharge is a critical stream temperature control process, it is not explicitly represented in many stream temperature models, an omission that may reduce predictive accuracy, hinder management of aquatic habitat, and decrease user confidence. We assessed the performance of a previously-described process-guided deep learning model of stream temperature in the Delaware River Basin (USA). We found lower accuracy (root mean square error [RMSE] of 1.71 versus 1.35°C) and stronger seasonal bias (absolute mean monthly bias of 1.06 vs. 0.68°C) for reaches primarily influenced by deep groundwater as compared to atmospheric conditions. We then tested four approaches for improving groundwater process representation: (a) a custom loss function leveraging the unique patterns of air and water temperature coupling characteristic of different temperature drivers, (b) inclusion of additional groundwater-relevant catchment attributes, (c) incorporation of additional process model outputs, and (d) a composite model. The custom loss function and the additional attributes significantly improved the predictive accuracy in groundwater-dominated reaches (RMSE of 1.37 and 1.26°C) and reduced the seasonal bias (absolute mean monthly bias of 0.44 and 0.48°C), but neither approach could identify holdout groundwater reaches. Variable importance analysis indicates the custom loss function nudges the model to use the existing inputs more efficiently, whereas with the added features the model relies on a broader suite of inputs. This analysis is a substantial step toward more accurately representing groundwater discharge processes in stream temperature models and will improve predictive accuracy and inform habitat management.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Ascorbic Acid Stability and pH Testing to Support RWM-018 Pilot Scale Wellhead Treatment System Feasibility Determination

A multitude of groundwater remediation techniques for chlorinated volatile organic compounds (cVOCs) have been conducted at the A/M-Areas of the Savannah River Site (SRS). Historical in-situ chemical oxidation (ISCO) injections of potassium permanganate and sodium persulfate conducted in 2018 and 2020, have left residual oxidant concentrations in the A/M-Area groundwater system. ISCO events targeted dense non-aqueous phase liquids (DNAPLs) through the injection of strong oxidants within plume zones upgradient of current groundwater recovery wells RWM-008 and RWM-018. These residual oxidants could significantly impede the operation of the mercury removal system needed to meet National Pollutant Discharge Elimination System (NPDES) permit requirements for the ongoing pump-and-treat system that hydraulically controls potions of the A/M-Area plume. Effluent treated groundwater from this system is discharged to surface water at the receiving outfall of the M-1 Air Stripper. To prevent interference of the M-1 Air Stripper system by residual oxidants, groundwater recovery has been ceased at well RWM-018. Modeling of oxidant transport in A/M-Area from ISCO estimated that the potential oxidant load contributed by RWM-018 could be as high as 120 mg/L. It is also possible that oxidants may be contributed via pumping at well RWM-008 in the near future as residual oxidants within the aquifer migrate toward its zone of influence. Savannah River National Laboratory (SRNL) has conducted research to investigate potential pre-treatment chemicals to neutralize these residual oxidants, before reaching the M-1 Air Stripper, to allow for resumed groundwater recovery at RWM-018 and continued pumping at RWM-008. Deployment of ascorbic acid for neutralization of residual oxidants at RWM-018 (and possibly RWM-008) is recommended based on results from previous reductant testing (SRNL, 2021). For a maximum estimated oxidant load at RWM-0018, a well-head treatment system would need to dose about 400 gallons of ascorbic acid each month to ensure complete neutralization. A second system with additional reagent should also be considered for deployment at RWM-008.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

The Straightening of a River Meander Leads to Extensive Losses in Flow Complexity and Ecosystem Services

To assist river restoration efforts we need to slow the rate of river degradation. This study provides a detailed explanation of the hydraulic complexity loss when a meandering river is straightened in order to motivate the protection of river channel curvature. We used computational fluid dynamics (CFD) modeling to document the difference in flow dynamics in nine simulations with channel curvature (C) degrading from a well-established tight meander bend (C = 0.77) to a straight channel without curvature (C = 0). To control for covariates and slow the rate of loss to hydraulic complexity, each of the nine-channel realizations had equivalent bedform topography. The analyzed hydraulic variables included the flow surface elevation, streamwise and transverse unit discharge, flow velocity at streamwise, transverse, and vertical directions, bed shear stress, stream function, and the vertical hyporheic flux rates at the channel bed. The loss of hydraulic complexity occurred gradually when initially straightening the channel from C = 0.77 to C = 0.33 (i.e., the radius of the channel is three-times the channel width), and additional straightening incurred rapid losses to hydraulic complexity. Other studies have shown hydraulic complexity provides important riverine habitat and is positively correlated with biodiversity. This study demonstrates how hydraulic complexity can be gradually and then rapidly lost when unwinding a river, and hopefully will serve as a cautionary tale.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Uncertainty quantification of bank vegetation impacts on the flood flow field in the American River, California, using large‐eddy simulations

Bank vegetation plays a key role in both hydrodynamics and morphodynamics of natural rivers; however, these effects are often unaccounted for in the computational flow dynamics of natural waterways. Recent studies using the large‐eddy simulation (LES), however, have attempted to gain insights into the impacts of bank vegetation on the mean flow field of the natural rivers using a vegetation model, which applies a sink term to the momentum equations of motion. This approach accounts for the effects of the vegetation and provides a practical approach to account for the complex patches of bank vegetation in large‐scale rivers. To implement the vegetation model, a drag coefficient reflecting the overall resistance of vegetal structures to the flow is needed, but due to the lack of calibrated data and range of size, density and type of vegetation, this parameter can be a significant source of uncertainty in the model results. Here, in this study, we use uncertainty quantification (UQ) to investigate the hydrodynamics and bed shear results when a bank vegetation is incorporated in an LES model. To this end, we used the polynomial chaos expansion and Monte Carlo sampling techniques to determine the uncertainties associated with the drag coefficient in the vegetation model and from uncertainties in the bed roughness and inflow discharge. The UQ analysis provided spatially varying confidence levels for the spanwise and vertical distribution of velocity magnitude and for the bed shear stress distributions. In addition, Sobol indices were computed to indicate the relative influence that each parameter had on the overall uncertainty. In general, it was found that uncertainty in flow discharge was the dominant source of uncertainty; however, the drag coefficient in the vegetation model and the bed roughness parameter also made significant contribution to the uncertainty near the banks and bed, respectively.

54 ENVIRONMENTAL SCIENCES↗