Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Groundwater Flux”

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

The importance of explicitly representing the streambed in watershed models

Abstract The streambed is the critical interface between the aquatic and terrestrial systems and hosts important biogeochemical hot spots within river corridors. Although the streambed characteristics are significantly different from those of its surrounding soil, the streambed itself has not been explicitly represented in watershed models. Here, we explicitly incorporated a streambed layer into an integrated hydrologic model through model parameterization and discretization. We examined the hydrological effects of streambed characteristics, including hydraulic conductivity ( K ), layer thickness, and resolution, on the exchange fluxes across the streambed as well as the streamflow at the watershed outlet. The numerical experiments were performed in the American River Watershed, a headwater, mountainous watershed within the Yakima River Basin in central Washington. Despite having a negligible effect on the watershed streamflow, an explicit representation of the streambed with distinctive properties dramatically changed the magnitude and variability of the exchange flux. In general, a larger streambed K along with a thicker streambed layer induced larger exchange fluxes. The exchange flux was most sensitive to the streambed resolution. A finer streambed resolution increased exchange fluxes per unit area while reducing the overall exchange volumes across the entire streambed. The amount of baseflow decreased by 6% as the streambed resolution increased from 250 to 50 m. This finding is important because these hydrological changes may, in turn, affect the exchange of nutrients and contaminants between surface water and groundwater and the associated biogeochemical processes. Our work demonstrated the importance of representing streambeds in fully distributed, process‐based watershed models to better capture the exchange flow dynamics in river corridors.

54 ENVIRONMENTAL SCIENCES↗

Flow Directions and Ages of Subsurface Water in a Salt Marsh System Constrained by Isotope Tracing

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Due to the counteracting terrestrial and marine forcings that control these environments, we do not comprehensively understand water fluxes in these complex coastal systems. To understand the water sources, flow directions, and velocities in salt marsh porewater, we employed a combination of geochemical tracers and analytical models across a hillslope-to-salt marsh continuum in a salt marsh experiencing daily inundation of estuarine surface water (SW) from tides and mixing of fresh seasonal groundwater. We used tritium ( 3 H) as a hydrologic tracer to assess porewater ages and stable water isotope ($\delta$ 2 H and $\delta$ 18 O) analyses to separate isotopically distinct estuarine and terrestrial groundwater across different depths and landscape positions in the study transect. We employed electrical conductivity to constrain the role of source mixing and evapotranspiration in salt marsh hydrology. Salinity and stable isotopes revealed that transpiration, rather than evaporation, increased subsurface water salinity to concentrations above estuarine SW during summer. Elevated salinity at depth indicated that salt marsh subsurface water is recharged during the dry growing season. Seasonal recharge patterns drive long-term deep subsurface water dynamics across the salt marsh, with 3 H ages of 3–7 years, and daily tidal cycles drive short-term shallow porewater dynamics with 3 H ages of 0 ± 3.6 years. Furthermore, our conceptual understanding of the spatiotemporal changes in SW-subsurface water interactions at the terrestrial-marine interface quantifies the hydrological constraints we are missing to improve our understanding of biogeochemical cycles within the salt marsh.

54 ENVIRONMENTAL SCIENCES↗

Water Level Data from Wells PLM1 and PLM6 for the East River Watershed, Colorado

This dataset (Williams et al., 2020) contains the original un-QA/QC-ed water level data for PLM1 and PLM6 and has been obsoleted. The data contained within this dataset is not to be used. Refer to Faybishenko et al., 2022 (DOI: 10.15485/1866836) for the latest QA/QC-ed data available via ESS-DIVE.This data set contains water level data for the PLM1 and PLM6 wells. PLM1 and PLM6 are location identifiers used by the Watershed Function SFA project for two groundwater monitoring wells along an elevation gradient located along the lower montane life zone of a hillslope near the Pumphouse location. These wells used to monitor subsurface water and carbon inventories and fluxes at the East River Watershed, Colorado, USA. Complete metadata information on the PLM1 and PLM6 wells are available in the related data package reference Varadharajan C, et al (2020). https://doi.org/10.15485/1660962.Data are reported in .csv files per well. The latitude and longitude of each location are given in a file called locations.csv. These data are used for determining the seasonally dependent flow of groundwater under the PLM hillslope. The downslope flow of groundwater in combination with data on groundwater chemistry can be used to estimate rates of solute export from the hillslope to the floodplain and river.These data products are part of the Watershed Function Scientific Focus Area collection effort to further scientific understanding of biogeochemical dynamics from genome to watershed scales.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Subsurface Flow and Solute Transport in a Snowmelt‐Recharged Hillslope With Multiyear Water Balance

Abstract Quantifying flow and transport from hillslopes is vital for understanding water quantity and quality in rivers, but remains obscure because of limited subsurface measurements. Using measured hydraulic conductivity K profiles and water balance over a single year to calibrate a transmissivity feedback model for a hillslope in the East River watershed (Colorado) proved unsatisfactory for predicting flow over the subsequent years. Well‐constrained field‐scale K were obtained by optimizing subsurface flux predictions over years having large differences in recharge, and by including estimates of interannual transfer of excess snowmelt recharge. Water and solute exports during high snowmelt recharge occur predominantly via shallow groundwater flow through weathered rock and soil because of their enlarged transmissivities under saturated conditions. Conversely, these shallow pathways are less active in snow drought years when the water table remains deeper within the weathering zone. Hillslope soil water monitoring showed that rainfall does not infiltrate deeply during summer and fall months, and revealed water losses consistent with model ET predictions. By combining water table‐dependent fluxes with pore water chemistry in different zones, time‐dependent rates of solute exports become predictable. As an example, calibrated K were combined with dissolved nitrogen concentrations in pore waters to show the snowmelt‐dependence of reactive nitrogen exported from the hillslope, further supporting the recent finding that the weathering zone is the dominant source of reactive nitrogen at this site. Subsurface export predictions can now be obtained for wide ranges of recharge based on measurements of water table elevation and profiles of pore water chemistry.

54 ENVIRONMENTAL SCIENCES↗

Water balance model (WBM) v.1.0.0: a scalable gridded global hydrologic model with water-tracking functionality

Abstract. This paper describes the University of New Hampshire Water Balance Model, WBM, a process-based gridded global hydrologic model that simulates the land surface components of the global water cycle and includes water extraction for use in agriculture and domestic sectors. The WBM was first published in 1989; here, we describe the first fully open-source WBM version (v.1.0.0). Earlier descriptions of WBM methods provide the foundation for the most recent model version that is detailed here. We present an overview of the model functionality, utility, and evaluation of simulated global river discharge and irrigation water use. This new version adds a novel suite of water source tracking modules that enable the analysis of flow-path histories on water supply. A key feature of WBM v.1.0.0 is the ability to identify the partitioning of sources for each stock or flux within the model. Three different categories of tracking are available: (1) primary inputs of water to the surface of the terrestrial hydrologic cycle (liquid precipitation, snowmelt, glacier melt, and unsustainable groundwater); (2) water that has been extracted for human use and returned to the terrestrial hydrologic system; and (3) runoff originating from user-defined spatial land units. Such component tracking provides a more fully transparent model in that users can identify the underlying mechanisms generating the simulated behavior. We find that WBM v.1.0.0 simulates global river discharge and irrigation water withdrawals well, even with default parameter settings, and for the first time, we are able to show how the simulation arrives at these fluxes by using the novel tracking functions.

54 ENVIRONMENTAL SCIENCES↗

Toward an Improved Understanding of the Global Fresh Water Budget

The major components of the global fresh water cycle include the evaporation from the land and ocean surfaces, precipitation onto the Ocean and land surfaces, the net atmospheric transport of water from oceanic areas over land, and the return flow of water from the land back into the ocean. The additional components of oceanic water transport are few, principally, the mixing of fresh water through the oceanic boundary layer, transport by ocean currents, and sea ice processes. On land the situation is considerably more complex, and includes the deposition of rain and snow on land; water flow in runoff; infiltration of water into the soil and groundwater; storage of water in soil, lakes and streams, and groundwater; polar and glacial ice; and use of water in vegetation and human activities. Knowledge of the key terms in the fresh water flux budget is poor. Some components of the budget, e.g. precipitation, runoff, storage, are measured with variable accuracy across the globe. We are just now obtaining precise measurements of the major components of global fresh water storage in global ice and ground water. The easily accessible fresh water sources in rivers, lakes and snow runoff are only adequately measured in the more affluent portions of the world. presents proposals are suggesting methods of making global measurements of these quantities from space. At the same time, knowledge of the global fresh water resources under the effects of climate change is of increasing importance and the human population grows. This paper provides an overview of the state of knowledge of the global fresh water budget, evaluating the accuracy of various global water budget measuring and modeling techniques. We review the measurement capabilities of satellite instruments as compared with field validation studies and modeling approaches. Based on these analyses, and on the goal of improved knowledge of the global fresh water budget under the effects of climate change, we suggest priorities for future improvements in global fresh water budget monitoring. The priorities are based on the potential of new approaches to provide improved measurement and modeling systems, and on the need to measure and understand the potential for a speed-up of the global water cycle under the effects of climate change.

Hildebrand, Peter H.↗

Permafrost Promotes Shallow Groundwater Flow and Warmer Headwater Streams

Abstract The presence of permafrost influences the flow paths of water through Arctic landscapes and thereby has the potential to impact stream discharge and thermal regimes. Observations from 11 headwater streams in Alaska showed that July water temperatures were higher in catchments with more near‐surface permafrost. We apply a fully coupled cryohydrology model to investigate if the impact of permafrost on flow path depth could cause the same pattern in temperatures of groundwater discharging from hillslopes to streams. The model simulates surface energy and water balances, snow, and subsurface water and energy balances for two‐dimensional hillslope model cases with varying permafrost extent. We find that hillslopes with continuous permafrost have more shallow flow paths and twice as high rates of evapotranspiration, compared to hillslopes with no permafrost. For our simulated cases, 6.7% of the horizontal water flux moves through the top organic soil layers when there is continuous permafrost, while only 0.5% moves through organic layers without permafrost. The deeper flow paths in permafrost‐free simulations buffer seasonal temperature extremes, so that summer groundwater discharge temperatures are highest with continuous permafrost. Our results suggest that permafrost thawing alters groundwater flow paths and can lead to decreases in summer stream temperatures and reductions in evapotranspiration in headwater catchments. These changes are of potential importance for stream biotic components of ecosystems, however, the full impact remains unknown.

54 ENVIRONMENTAL SCIENCES↗

Small, Untethered, Mobile Robots for Inspecting Gas Pipes

Small, untethered mobile robots denoted gas-pipe explorers (GPEXs) have been proposed for inspecting the interiors of pipes used in the local distribution natural gas. The United States has network of gas-distribution pipes with a total length of approximately 109 m. These pipes are often made of iron and steel and some are more than 100 years old. As this network ages, there is a need to locate weaknesses that necessitate repair and/or preventive maintenance. The most common weaknesses are leaks and reductions in thickness, which are caused mostly by chemical reactions between the iron in the pipes and various substances in soil and groundwater. At present, mobile robots called pigs are used to inspect and clean the interiors of gas-transmission pipelines. Some carry magnetic-flux-leakage (MFL) sensors for measuring average wall thicknesses, some capture images, and some measure sizes and physical conditions. The operating ranges of pigs are limited to fairly straight sections of wide transmission- type (as distinguished from distribution- type) pipes: pigs are too large to negotiate such obstacles as bends with radii comparable to or smaller than pipe diameters, intrusions of other pipes at branch connections, and reductions in diameter at valves and meters. The GPEXs would be smaller and would be able to negotiate sharp bends and other obstacles that typically occur in gas-distribution pipes.

Wilcox, Brian↗

Revised Kimberlina Groundwater Simulations

Kimberlina Groundwater Qaulity Simulations (Version 3: 2015 October) by Lawrence Livermore National Labratory. Submitted to the EDX team at the National Energy Technology Laboratory for research use. This dataset contains 562 updated groundwater quality simulations based on the Kimberlina site. It is a six-layer model including the aquifers and aquitards inbetween the ground surface and the reservoir. The input files and associated leakage information for these simulations are included in the zip file "nrap_kim_well31_fluxes and inputs.tar.gz". These 562 simulations are the successful results out of 1000 monte carlo simulations, with time length greater than 5 years. Each simulation consists 6 output files, therefore, there are 3372 files (6*562) in total. One single example simulation of full time length (200 year) is also provided in the zip file "single_simulation_200yr_sim0012.tar.gz". These 6 groundwater quality parameters are: - pH (use absolute value) - TDS (mg/kg) - Pressure (Pa) - CO2 Liquid Saturation (fraction range 0.0-1.0) - CO2 Gas Concentration (mol/L) - CO2 Liquid Concentration (mol/kg) Their correponding file names, take the first simulation for example, are: ph0001.ntab (for pH) , tds0001.ntab (for TDS), usnt0001.p.ntab (for pressure), usnt0001.s_liquid.ntab (for liquid saturation), usnt0001.c_co2_gas.ntab (for CO2 gas concentration) and usnt0001.c_co2_liquid.ntab (for CO2 liquid concentration). Each output file contains at least 19 columns (variables). The total number of columns depends on the simulation time length. The column data include the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the jth grid of y-axis 4. k: the kth grid of z-axis 5. element_ref: element type reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The rest (14, 15, 16...) variables are the simulated parameter values over time, take pH as an example, are: 14. ph000y: initial pH value. 15. ph001y: simulated pH value at the end of the 1st year. 16. ph002y: simulated pH value at the end of the 2nd year. ... (and so on)... The zip input files have the input fluxes of brine and CO2 (in kg/s) allocated to model layers 1,2,3 and 5 (layers 4 and 6 are assumed to be no-flow aquitards) through the following nodes: (src1 WELL 7 7 25 25 10 32) #Etchegoin (src2 WELL 7 7 25 25 33 72) #Macoma Chanac (src3 WELL 7 7 25 25 73 83) #Santa Margarita Mclure (src4 WELL 7 7 25 25 84 89) #Fruitvale Round Mountain (aquitard) (src5 WELL 7 7 25 25 90 96) #Olcese Sand (thief zone) (src6 WELL 7 7 25 25 97 99) #tremblor Freeman (aquitard) ex: brn1 is assigned to src 1 from i=(7,7), j=(25,25), k=(10,32) Edited by Ya-Mei Yang at National Energy Technology Laboratory

brine leakage↗

Quantifying Parameter Sensitivity, Interaction and Transferability in Hydrologically Enhanced Versions of Noah-LSM over Transition Zones

We use sensitivity analysis to identify the parameters that are most responsible for shaping land surface model (LSM) simulations and to understand the complex interactions in three versions of the Noah LSM: the standard version (STD), a version enhanced with a simple groundwater module (GW), and version augmented by a dynamic phenology module (DV). We use warm season, high-frequency, near-surface states and turbulent fluxes collected over nine sites in the US Southern Great Plains. We quantify changes in the pattern of sensitive parameters, the amount and nature of the interaction between parameters, and the covariance structure of the distribution of behavioral parameter sets. Using Sobol s total and first-order sensitivity indexes, we show that very few parameters directly control the variance of the model output. Significant parameter interaction occurs so that not only the optimal parameter values differ between models, but the relationships between parameters change. GW decreases parameter interaction and appears to improve model realism, especially at wetter sites. DV increases parameter interaction and decreases identifiability, implying it is overparameterized and/or underconstrained. A case study at a wet site shows GW has two functional modes: one that mimics STD and a second in which GW improves model function by decoupling direct evaporation and baseflow. Unsupervised classification of the posterior distributions of behavioral parameter sets cannot group similar sites based solely on soil or vegetation type, helping to explain why transferability between sites and models is not straightforward. This evidence suggests a priori assignment of parameters should also consider climatic differences.

Rosero, Enrique↗

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗

Hydrological Controls on Riverbed Methane Emissions: A Numerical Investigation of Hydrodynamic and Ebullitive Mechanisms from Site to Basin Scales

Methane production and emission from riverbed sediments constitute a significant yet underexplored component of the river methane budget. Despite their importance, a mechanistic and quantitative understanding of these processes under dynamic flow conditions remains elusive, particularly at the basin scale. In this study, we investigate the multiscale (site-to-basin) mechanisms governing methane emissions from riverbeds, focusing on the interplay between biogeochemical processes and hydrological dynamics. We developed a reactive transport model that integrates methane production with hydrodynamic and ebullitive pathways and coupled it with a basin-scale groundwater flow model to simulate methane emissions across scales. By incorporating key factors such as vertical hydrological exchange flow (VHEF), particulate organic carbon availability, sediment properties, and temperature sensitivity, our model captures the spatiotemporal variability of methane ebullition at both the site and basin scales. Our results reveal that methane ebullition is strongly modulated by VHEF, with lower flow rates promoting methane accumulation and subsequent ebullition. Sensitivity analyses underscore nonlinear responses of methane fluxes to environmental drivers. Furthermore, these findings emphasize the critical roles of sediment composition, hydrological exchange, and environmental conditions in shaping methane dynamics. Overall, this work advances our understanding of riverine methane emissions and offers insights for guiding climate models and mitigation strategies.

Fluxes↗

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗

Effects of groundwater pumping on ground surface temperature: A regional modeling study in the North China Plain

Over-exploitation of groundwater (GW) in the North China Plain (NCP) since the 1960s has many environmental consequences. However, understanding of the dominant mechanisms remains limited, particularly at regional scale. In this study, the coupled ParFlow.CLM model representing subsurface and land-surface processes and their interactions was applied in the NCP at high spatio-temporal resolutions. The model was validated using the water and energy fluxes reported in previous studies and from the JRA-55 reanalysis. Numerical experiments were designed to examine the relative impacts of GW pumping and irrigation on the ground surface temperature (GST). Results show significant effects of GW pumping on GST in the NCP. Generally, the subsurface acts as a buffer to temporal variations in heat fluxes at the land-surface, but long-term pumping can gradually weaken this buffer, resulting in increases in the spatio-temporal variability of GST, as exemplified by hotter summers and colder winters. Considering that changes of water table depth (WTD) can significantly affect land surface heat fluxes when WTD ranges roughly between 0.01–10 m, the 0.5 m/year increase of WTD simulated by the model due to pumping can continue to increase the average WTD and hence, disturb GST, for about 20 years from the pre-pumping WTD in the NCP, before the WTD exceeds 10 m. The variations of GST are expected to be faster initially and gradually slow down due to the nonlinear behaviors of GST with WTD. The findings from this study in the NCP may also have implications for other regions with GW depletion.

Yang, Chen↗

Prototype Hanford Barrier Performance Monitoring Report: Fiscal Year 2024

This report provides an annual update on performance monitoring of the Prototype Hanford Barrier (PHB) at the Hanford Site. The PHB is part of a long-term study that serves as the scientific basis for many of the evapotranspiration-capillary barrier designs currently used globally and for future engineered barrier designs planned for remedial actions over waste and demolition sites at the Hanford Site. The PHB allows us to identify potential future issues and develop better monitoring techniques before these engineered barriers are constructed. Surface barriers like the PHB are essential for preventing water infiltration and curbing the spread of contaminants to groundwater. Effective monitoring of soil moisture levels above and below these barriers is crucial given that performance metrics could span up to 1,000 years. From July 2023 to June 2024, the water flux (as measured by tipping buckets) through the 2-m-thick silt loam layer of the PHB – a fine silty material that stores water under high capillary tension – remained well below the 0.5-mm-per-year performance threshold, demonstrating its effectiveness in preventing water penetration. In 2024, degraded tipping bucket gauges were replaced to ensure the accuracy of future data. Additionally, neutron probes revealed that the wetting front from the rainy season only penetrated to a maximum depth of 1.2 m into the silt loam. This further demonstrated the barrier's efficiency, as the wetting front did not fully penetrate the 2-m-thick capillary barrier. The western gravel slope of the PHB, composed of Hanford Site pit gravel of varying sizes, demonstrated low flux rates. In contrast, the eastern riprap slope exhibited higher flux rates. Zhang (2017) found that the riprap side slope of the PHB had the highest drainage rates in January and the lowest in late summer or early fall, indicating significant summer evaporation. Drainage rates increased significantly under enhanced precipitation conditions, far exceeding the design criterion, which could lead to water infiltration into the waste zone. The study introduced the “edge effect,” where elevated drainage rates from the riprap may migrate laterally beneath the barrier, compromising its ability to isolate waste, and recommended expanding the barrier and conducting further research to mitigate this issue. This report provides a review on surface barrier edge effects and their potential impact to subsurface contaminant migration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Can Simple Machine Learning Tools Extend and Improve Temperature-Based Methods to Infer Streambed Flux?

Temperature-based methods have been developed to infer 1D vertical exchange flux between a stream and the subsurface. Current analyses rely on fitting physically based analytical and numerical models to temperature time series measured at multiple depths to infer daily average flux. These methods have seen wide use in hydrologic science despite strong simplifying assumptions including a lack of consideration of model structural error or the impacts of multidimensional flow or the impacts of transient streambed hydraulic properties. We performed a “perfect-model experiment” investigation to examine whether regression trees, with and without gradient boosting, can extract sufficient information from model-generated subsurface temperature time series, with and without added measurement error, to infer the corresponding exchange flux time series at the streambed surface. Using model-generated, synthetic data allowed us to assess the basic limitations to the use of machine learning; further examination of real data is only warranted if the method can be shown to perform well under these ideal conditions. We also examined whether the inherent feature importance analyses of tree-based machine learning methods can be used to optimize monitoring networks for exchange flux inference.

54 ENVIRONMENTAL SCIENCES↗

Contaminant Transport Modeling for Technology Evaluation and Long-Term Monitoring in the Tims Branch Testbed, SC - 20343

Studies conducted at U.S. DOE sites have shown the presence of heavy metals, radionuclides, and volatile organic compounds in surface water, groundwater, and soil as a result of nuclear activity in the Cold War era. Since the 1990's, innovative cleanup methods have been implemented in the Tims Branch watershed at Savannah River Site (SRS) to limit the contaminant flux to the stream that have reduced the contaminant concentrations to acceptable regulatory levels in the dissolved phase. A tin-based treatment which effectively eliminated all local anthropogenic mercury inputs to this ecosystem resulted in a known step function addition of inert tin oxide particles which now serve as a potential tracer for sedimentation and particle transport processes in the stream. The long-term effectiveness of this and other remediation techniques and the potential for remobilization of adsorbed contaminant in sediment during extreme hydrologic conditions however remains unclear. It is therefore important to understand not only the fate and transport of dissolved contaminants, but also the movement of sediment and the relevant interactions with dissolved contaminant. To narrow this knowledge gap, a study is being conducted using the Tims Branch watershed as a stream-scale ecosystem test-bed to identify the primary transport processes of major contaminants of concern (such as mercury, nickel and uranium) with an emphasis on interactions with sediment transport. This involves the development of a fully distributed hydrologic watershed model of the Tims Branch watershed to predict streamflow under extreme weather conditions, as well as the development of a comprehensive contaminant transport model that can properly account for coupled contaminant and sediment transport. Review of relevant research reports and peer-reviewed journals revealed that aside from advection-dispersion transport of dissolved contaminants, adsorption and desorption with suspended solids and bed sediment also play an important role in the transport of those contaminants of concern. To develop the fully distributed hydrologic watershed model, the MIKE SHE 2-dimensional (2D) land surface/3D groundwater model that simulates surface/subsurface hydrologic processes (such as overland flow, evapotranspiration, and infiltration) was coupled with a 1D streamflow model that accounts for stream water hydraulics (such as hydraulic structures, cross-sections, and network). To model contaminant transport, the MIKE 11 streamflow component was coupled with the MIKE 11 AD module that simulates solute transport through advection and dispersion, and MIKE ECO Lab module that accounts for both sediment transport and interactions with dissolved contaminant. At this stage, the development and optimization of the fully distributed hydrologic model has been completed, achieving satisfactory statistical results between observed and predicted discharge as indicated by a root mean square error (RMSE) of 0.039 cms and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.764. The ongoing development of the contaminant transport model has also yielded realistic results from preliminary tests. Results from this study are a key to evaluating the effectiveness of tin (II)-based mercury treatment of wetlands at the SRS site, and are also relevant to evaluating the potential of using this type of novel remediation technology in other mercury-contaminated stream systems. Knowledge acquired from this research will also support interpretation of historical data on the trends of contaminant concentration distribution in Tims Branch, particularly considering the effect of extreme hydrological events on the stream flow and pollutant transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Scale-dependent spatial variabilities of hydrological exchange flows and transit time in a large regulated river

Hydrological exchange flows (HEF) across the river-aquifer interface and the associated residence time of river water in the aquifer have important implications for contaminant plume migration and biogeochemical processes in the river corridor. HEFs and residence time are influenced by both subsurface physical features and hydrologic forcing related to the transport process, which can exhibit complex spatial and temporal variations. In this study, we used a massively parallel subsurface flow model and a particle-tracking model to study the influences of different control factors on spatial variability of HEFs and residence time distributions (RTD) in the Hanford Reach of the Columbia River in Washington State. A total number of 100M particles were randomly injected in time and space and then tracked in a model domain that covers a 51-km 2 area (15.1M model cells). We used hourly river stages and groundwater levels to drive the model to provide dynamic velocity fields for the particle tracking in the simulation period that was longer than 2 years. The groundwater flow simulation and particle-tracking results provide the first comprehensive assessment of the spatial distribution of HEFs and residence time in large complex river corridors. Overall, our results show that the aquifer hydrogeological structure has the strongest correlation with the extent and magnitude of exchange flux. The residence time exhibits complex patterns that are impacted by all the river geomorphologic, hydrodynamic, and hydrogeologic factors and are strongly correlated with the downwelling ratio of exchange flux. The new insights gained through this study can be used to support the development of reduced-order models of HEFs and RTDs for large complex river systems.

54 ENVIRONMENTAL SCIENCES↗