Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrologic exchange flows”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗
DOE OSTI2024

Hydrological controls of a riparian wetland based on stable isotope data and model simulations

Isotopic evidence of groundwater and stream water is frequently used to investigate water exchanges with groundwater. Monthly sampling of rain, stream water, and groundwater was conducted at Tims Branch watershed in South Carolina for the oxygen and hydrogen stable isotope (δ 2 H and δ 18 O) measurement, as well as pH and oxidation–reduction potential (ORP). Together with a mass balance perspective, it was determined that it takes a few weeks to one month for groundwater in the hyporheic zone to fully exchange with stream water. From hydrodynamic modelling, we show that substantial (up to 70 %) groundwater exchange occurs at gaining and losing sites. Groundwater exfiltration, i.e. inflow into stream water, contributes up to 4 % to stream water, with the remainder from upstream exfiltration. A 2–4 % per day renewal rate of adjacent groundwater would indirectly indicate a groundwater residence time in the order of half a month to a full month (assuming either a well-mixed case or large dispersion rate in pulse flow case), in agreement with a greatly reduced variability of δ 2 H and δ 18 O of groundwater compared to stream water and rain. This reduced variability of stable isotope signal from groundwater confirms our hypothesis that riparian groundwater mixing at Tims Branch is more of a mixed type rather than a pulse flow type. As a result, a monthly time scale is sufficient for groundwater to become anoxic at exit points into stream water resulting in the episodic production of natural organic matter- and iron-rich flocs upon oxidation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hydrologic connectivity and dynamics of solute transport in a mountain stream: Insights from a long-term tracer test and multiscale transport modeling informed by machine learning

The movement of solutes in a watershed is a complex process with multiple interactions and feedbacks across spatial and temporal scales. Modeling the dynamics of solute transport along diverse hydrologic pathways within watersheds – from hillslopes to stream channels and in and out of the hyporheic zones – is challenging but critically important, as these processes integrate and contribute to the biogeochemical functioning of the river corridor up to the river network scale. Here we use results from a long-term network-scale tracer test at the H.J. Andrews experimental forest in western Cascade Mountains, Oregon, USA to inform a multiscale framework for transport in stream corridors. The framework uses a Lagrangian-based subgrid model to represent the effects of hyporheic exchange flow and advective transport at stream network scales. The spatially and temporally resolved stream discharge needed for the transport model is imputed across the river system by an entity-aware long short-term memory network. Modeled concentrations show good agreements with the observations and exhibit power scaling laws indicative of a very wide range of timescales over which hyporheic exchange flow occurs. Our results demonstrate a data-informed modeling framework that links dynamical processes occurring at small scales to a network context to help understand how changes at reach scale cascade into network-scale effects, providing a useful tool for sustainable river basin management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impacts of Sea‐Level Rise on Coastal Groundwater Table Simulated by an Earth System Model With a Land‐Ocean Coupling Scheme

Abstract Sea‐level rise (SLR) poses a severe threat to the coastal environment through seawater intrusion into freshwater aquifers. The rising groundwater table also exacerbates the risk of pluvial, fluvial, and groundwater flooding in coastal regions. However, current Earth system models (ESMs) commonly ignore the exchanges of water at the land‐ocean interface. To address this gap, we developed a novel land‐ocean hydrologic coupling scheme in a state‐of‐the‐science ESM, the Energy Exascale Earth System Model version 2 (E3SMv2). The new scheme includes the lateral exchange between seawater and groundwater and the vertical infiltration of seawater driven by the SLR‐induced inundation. Simulations were performed with the updated E3SMv2 for the global land‐ocean interface to assess the impacts of SLR on coastal groundwater under a high CO 2 emission scenario. By the middle of this century, seawater infiltration on the inundated areas will be the dominant component in the land‐ocean coupling process, while the lateral subsurface flow exchange will be much smaller. The SLR‐induced seawater infiltration will raise the groundwater levels, enhance evapotranspiration, and increase runoff with distinct spatial patterns globally in the future. Although the coupling process is induced by SLR, we found topography and warming temperature have more control on the coupling impacts, probably due to the relatively modest magnitude of SLR during the selected future period. Overall, our study suggests significant groundwater and seawater exchange at the land‐ocean interface, which needs to be considered in ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Toward more mechanistic representations of biogeochemical processes in stream and river networks: Implementation and demonstration of a multiscale model : Modeling Archive

This Modeling Archive supports the ORNL-SFA and IDEAS-Waterhsed publication submitted to Environmental Modeling and Software. A recently introduced multiscale model for representing the combined effect of hyporheic exchange flows and small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The ADELS (Advection Dispersion Equation with Lagrangian Subgrid) model uses a one-dimensional advection-dispersion-reaction equation for the channel network and couples that equation at each channel location to a one-dimensional advection-reaction subgrid model representing an ensemble of streamlines that are diverted into the hyporheic zone before returning to the channel. The subgrid model is written in a computationally advantageous Lagrangian form with hyporheic age replacing the hyporheic travel distance. In the paper, we summarized implementation in the integrated surface/subsurface hydrology modeling system Advanced Terrestrial Simulator (ATS).

East Fork Poplar Creek↗
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↗
NASA NTRS2011

An Integrated Hydrological and Water Management Study of the Entire Nile River System - Lake Victoria to Nile Delta

The Nile basin River system spans 3 million km(exp 2) distributed over ten nations. The eight upstream riparian nations, Ethiopia, Eretria, Uganda, Rwanda, Burundi, Congo, Tanzania and Kenya are the source of approximately 86% of the water inputs to the Nile, while the two downstream riparian countries Sudan and Egypt, presently rely on the river's flow for most of the their needs. Both climate and agriculture contribute to the complicated nature of Nile River management: precipitation in the headwaters regions of Ethiopia and Lake Victoria is variable on a seasonal and inter-annual basis, while demand for irrigation water in the arid downstream region is consistently high. The Nile is, perhaps, one of the most difficult trans-boundary water issue in the world, and this study would be the first initiative to combine NASA satellite observations with the hydrologic models study the overall water balance in a to comprehensive manner. The cornerstone application of NASA's Earth Science Research Results under this project are the NASA Land Data Assimilation System (LDAS) and the USDA Atmosphere-land Exchange Inverse (ALEXI) model. These two complementary research results are methodologically independent methods for using NASA observations to support water resource analysis in data poor regions. Where an LDAS uses multiple sources of satellite data to inform prognostic simulations of hydrological process, ALEXI diagnoses evapotranspiration and water stress on the basis of thermal infrared satellite imagery. Specifically, this work integrates NASA Land Data Assimilation systems into the water management decision support systems that member countries of the Nile Basin Initiative (NBI) and Regional Center for Mapping of Resources for Development (RCMRD, located in Nairobi, Kenya) use in water resource analysis, agricultural planning, and acute drought response to support sustainable development of Nile Basin water resources. The project is motivated by the recognition that accurate, frequent, and spatially distributed estimates of the water balance are necessary for effective water management. This creates a challenge for watersheds that are large, include data poor regions, and/or span multiple nations. All of these descriptors apply to the Nile River basin, yet successful management of the Nile is critical for development and political stability in the region. For this reason, improved hydrological data to support cooperative water management in the Nile basin is a priority for USAID, the US State Department, the World Bank and other international organizations. In this project, the U.S. based research team is working with partners at RCMRD, Nile Basin Initiative (NBI), and their member national-level agencies to develop satellite-based land cover maps, satellite-derived evapotranspiration estimates (using the ALEXI algorithm), and NASA's Land Data Assimilation System (LDAS) customized to match identified information needs. The cornerstone applied sciences product of the project is the development of a customized "Nile LDAS" that will produce optimal estimates of hydrological states and fluxes, as vetted against the in situ observations of NBI and RCMRD member organizations and independent satellite-derived hydrological estimates. Nile LDAS will be applied to improve the reliability of emerging Decision Support Systems in applications that include drought monitoring, reservoir management, and irrigation planning. The end-users such as RCMRD, NBI, Ethiopian and Kenya Meteorological and Famine Early Warning System Network (FEWSNet) will be the eventual benefactors of this work. There will be a capacity building process involving the above end-user organizations and transfer the models and the results for these organizations to execute for future use. The team has already initiated this study and the early results of first years' work are shown. The plan is to complete this work by late 2013.

Habib, Shahid↗
DOE OSTI2026

Xanthos-Lake Model Source Code

This repository contains the source code for Xanthos-Lake, a lake-modeling extension of the Xanthos framework that introduces a coupled lake component comprising the Xanthos-Lake Snow and Ice Model (xLSIM) and the Xanthos-Lake Water Balance Model (xLWBM). xLSIM is a basin-aware machine-learning model for lake snow, ice, and thermal conditions. It predicts monthly lake ice thickness, snow depth, snow-cover fraction, mixing-layer temperature, and lake ice fraction from meteorological forcing and lake surface-area information. It uses sequence-based deep-learning architectures, including Transformer and hybrid Long Short-Term Memory–Transformer (LSTM–Transformer) models, together with seasonal encoding, multi-lake learning, physical masking, and basin-level cryospheric and non-cryospheric classification. The training workflow uses Ray for scalable execution and includes optional Ray Tune hyperparameter optimization. Model predictions, observations, diagnostics, and feature-importance outputs are written in NetCDF. xLWBM is the water-balance component of the new lake framework. It simulates monthly lake storage, surface area, evaporation, inflow, outflow, and lake–groundwater exchange. It combines physical water-balance equations with calibrated bathymetric relationships, weir-based outlet flow, modified Penman open-water evaporation, groundwater head relaxation, Penman–Monteith snow and ice sublimation, and snow, ice, and thermal conditions supplied by xLSIM. The model calibrates lake parameters against satellite-derived surface-area data, using evaporation-based calibration where surface-area data are unavailable, and supports small, medium, and large lake classes. For large lakes, xLWBM is integrated with the managed-routing workflow so that lake storage and outflow interact directly with downstream river routing and reservoir operations. Together, xLSIM and xLWBM provide Xanthos with a coupled lake-modeling capability. xLSIM supplies the snow, ice, and thermal conditions that affect lake evaporation and snow- and ice-related water exchanges, while xLWBM translates those conditions into dynamic lake storage, surface area, evaporation, and discharge. In return, xLWBM supplies evolving lake surface area to xLSIM. This coupling enables Xanthos to represent lakes as active hydrologic components within basin-scale water-availability and routing simulations.

Machine Learning↗
DOE OSTI2023

Dissolved Oxygen and Temperature Data from the Hyporheic Zone of the East River Watershed July 2017 to October 2018

Dissolved oxygen (DO) is critical for aquatic ecosystems. Our focus is on the long-term DO dynamics in hyporheic zone of rivers, which are a function of both transport (hydrologic exchange between river and hyporheic zone) and uptake by biogeochemical reactions or respiration. The study site is the alpine East River watershed in Colorado, USA, meander A downstream from the pump house. Opti O2 probes were deployed in the water column and directly within the river-bed at 10, 20, and 35 cm depth (38°55'23.38"N, 106°57'4.21"W, 9046.88m) to monitor DO and temperature. A continuous data stream from July 24, 2017 to Oct 24, 2018 was autonomously telemetered to the cloud. This 14-month DO and temperature time series were obtained without any servicing for maintenance or data downloads; additionally the ability to remotely verify probe performance during field deployment was essential to confirm data validity during winter freeze-in and hydrological/weather events, such as spring melt and summer monsoons. We investigate the variations in dissolved oxygen dynamics of this snow-pack dominated watershed during a comparatively low flow water year (2018) and a relatively normal water year (2017), enabled by distinctive, in-situ, high frequency (∆t = 5min) sensors that provided a continuous time-series from the undisturbed study site over multiple seasons.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Spatial Study 2022: Water Column, Sediment, and Total Ecosystem Respiration Rates across the Yakima River Basin, Washington, USA (v2)

This dataset supports a broader study examining the drivers of spatial variability in sediment respiration rates in the Yakima River Basin and is associated with the manuscript “Sediment-associated processes account for most of the spatial variation in ecosystem respiration in the Yakima River basin” submitted to Nature Communications Earth & Environment (Garayburu-Caruso et al., in review). The dataset provides ecosystem metabolism estimates generated from streamMetabolizer (Appling et al.; 2018) using data collected during the same five-week period at 48 sites within multiple rivers throughout the Yakima River Basin in Washington, USA. Additionally, it includes the scripts used for the analysis and producing the figures in the manuscript. The contents include streamMetabolizer inputs and outputs and additional relevant data needed to generate the main manuscript results. The data included are: total ecosystem respiration, water respiration, calculated sediment-associated respiration, gross primary production outputs from the river corridor model for the Yakima River Basin, median grain size (d50), depth, dissolved oxygen, water temperature, pressure, and annual oxygen consumption. The associated GitHub repository can be found at https://github.com/river-corridors-sfa/SSS_metabolism. Samples collected during this study were labeled as “Second Spatial Study” or “SSS.” Raw time series sensor data, total suspended solids, and depth data from SSS were published at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1969566. A subset of data from the SSS samples were published in the contiguous United States (CONUS)-Scale Model-Sample (CM) study data package available at https://data.ess-dive.lbl.gov/view/doi:10.15485/1923689 that presents data from across the CONUS. They include dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC), total nitrogen (TN), grain size, aerobic sediment respiration, dissolved oxygen (DO), and temperature. Parent IDs and Site IDs are consistent between the SSS and CM data packages, and they can be mapped directly so data across packages can be used together. Field metadata for the samples in this da This dataset is comprised of one main data folder with four subfolders. The main data folder contains of (1) file-level metadata; (2) data dictionary; (3) total/water column/sediment respiration; (4) gross primary production (GPP); (5) median grain size (d50); and (6) annual oxygen consumption. The “Figures” subfolder contains the figures used in the paper and all intermediate files (including geospatial files). The “Published_Data” contains a readme directing the user to download the public data to reproduce analyses and figures. The “Scripts” folder contains all scripts used in the analyses that were not part of running StreamMetabolizer. Lastly, the “Stream_Metabolizer” folder contains all files associated with running StreamMetabolizer including (1) model input files, (2) model output files, (3) processing scripts, (4) histogram plots of the outputs, and (5) an R project. All files are .csv, .pdf, .R, .Rmd, .Rproj, .html, .png, .txt, .qgz, .cpg, .dbf, .prj, .shp, .shp.ea.iso.xml, .shp.iso.xml, .shx, .sbn. ta package can be found at either link. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hydrologic windows into the crystalline basement and their controls on groundwater flow patterns across the Paradox Basin, western USA

Conceptual models of sedimentary basin groundwater flow systems typically assume that the crystalline basement acts as an impermeable boundary and can be neglected. In this study, we use hydrologic models constrained by isotopic and geochemical datasets to argue that the La Sal Mountains, Utah, USA, act as a hydrologic window into the Paradox Basin’s lower aquifer system and underlying crystalline basement. We conducted a sensitivity study in which we varied crystalline basement/laccolith permeability as well as fault zone connectivity along a cross-sectional transect from the La Sal Mountains to Lisbon Valley. When the crystalline basement/laccolith units are set at relatively permeable levels (10–14 m 2 ), simulated tracers that include total dissolved solids, oxygen isotopic composition of pore fluids (δ 18 O), and groundwater residence times are in closest agreement with field measurements. Model results indicate that pore fluids in the basal aquifer system underlying the Paradox Formation confining unit are a mixture of relatively young meteoric fluids and older Paradox Formation brines. The presence of faults did not significantly modify fluid exchange between the upper and lower aquifer systems. This was due, in part, to underpressuring within the Paradox Formation. Our study concludes that the Paradox Basin represents a regional recharge area for the Colorado Plateau, with groundwater discharge occurring along the Colorado River within the Grand Canyon some 375 km away to the southwest. This is only possible with a permeable crystalline basement. Our findings help explain the genesis of Mississippi Valley-type ore deposits of the US Midcontinent, where the presence of a permeable basement may be useful in addressing issues related to solute mass and energy balance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Interactions between molecular-scale biogeochemical processes and hyporheic exchange for understanding coupled Fe-S-C cycling in iron-rich riparian wetlands (Final Technical Report)

Riparian wetlands are dynamic interfaces that exert strong control over water quality, contaminant mobility, and greenhouse gas emissions. These systems are characterized by hyporheic exchange between oxic surface water and anoxic groundwater, which generates steep redox gradients and promotes spatially and temporally variable microbial activity. Despite growing recognition of tightly coupled iron (Fe), sulfur (S), and carbon (C) cycling in these environments, the mechanisms governing these interactions—particularly under low-sulfate freshwater conditions—remain poorly constrained. This project developed a mechanistic understanding of how hydrologic variability and microbial processes interact to control Fe–S–C cycling in iron-rich riparian wetlands. Using a multi-scale and multi-method approach integrating field observations, geochemical and spectroscopic analyses, metagenomics, and reactive transport modeling, we demonstrate that “cryptic” sulfur cycling—rapid sulfur transformations involving intermediate-valence species—plays a dominant and previously underrecognized role in freshwater wetlands. These processes persist despite low sulfate concentrations and significantly influence iron reduction, carbon mineralization, and methane dynamics. The results show that cryptic sulfur cycling enhances dissolved Fe 2+ production, regulates methane concentrations, and is strongly controlled by climate-driven hydrologic fluxes. By linking hydroclimate, subsurface flow, and biogeochemical reactions, this work provides a predictive framework for understanding how wetland systems respond to environmental change, with direct implications for water quality and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data and script associated with “Shifts in Rain-Snow Partitioning Drive Faster Water Transit Times in the US Pacific Northwest”

This data package contains the data and code to use and run the Water Tracer enabled version of the Weather Research and Forecasting Hydrologic model (WT-WRF-Hydro) with the Sequential Precipitation Input Tagging (SPIT) framework. It is associated with the publication “Shifts in Rain-Snow Partitioning Drive Faster Water Transit Times in the US Pacific Northwest” published in Scientific Reports (Butler et al., 2026; https://doi.org/10.1038/s41598-026-46539-1). We use the Continental U.S. (CONUSII; Rasmussen et al., 2021) dataset to force the model with an historical climate (2006–2013) and a future climate (2086–2093) with a representative carbon pathway (RCP) 8.5 scenario. We use the model to calculate water transit times in five headwater catchments within the U.S. Pacific Northwest. We also show key hydrologic and environmental variables that affect water transit times and changes in the future. Finally, we use observed data to validate the model such as stream water isotopes, snowpack characteristics, and stream discharge. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. The data package consists of 11 folders: (1) "Figures" contains the exported figures used in the manuscript; (2) "Model_Isotope_Date" contains the WT-WRF-Hydro isotope date used in model validation; (3) “Model_Outputs_Future” contains the WT-WRF-Hydro future climate outputs; (4) “Model_Outputs_Historical” contains the WT-WRF-Hydro historical climate outputs; (5) “Model_Outputs_Weights_Areas” contains the WT-WRF-Hydro weights per catchment used to calculate water transit times and isotopes in stream water; (6) “MODIS_data_scripts” contains data used to validate snow conditions in the study area; (7) “Observed_Flow_Data” contains the observed streamflow data used in model validation; (8) “Observed_Isotope_Data” contains the observed stream water isotope data used in model validation; (9) “Scripts” contains the Python scripts used to general results and the figures; (10) “Statistic_Outputs” contains the water transit time statistical outputs reported in this manuscript; (11) “Validation_SNOTEL” contains the SNOTEL data used in model validation. The files in this data package have the following file extensions: .tif, .txt, .csv, .pdf, .py, .jpg, and .png.

American River↗
DOE OSTI2023

Detecting Permafrost Active Layer Thickness Change From Nonlinear Baseflow Recession

We report permafrost underlies about one fifth of the global land area and affects ground stability, freshwater runoff, soil chemistry, and surface-atmosphere gas exchange. The depth of thawed ground overlying permafrost (active layer thickness) has broadly increased across the Arctic in recent decades, coincident with a period of increased streamflow, especially the lowest flows (baseflow). Mechanistic links between active layer thickness and baseflow have recently been explored using linear reservoir theory, but most watersheds behave as nonlinear reservoirs. We derive theoretical nonlinear relationships between long-term average saturated soil thickness $\bar{η}$(proxy for active layer thickness) and long-term average baseflow. When applied to 38 years of daily streamflow data for the Kuparuk River basin on the North Slope of Alaska, the theory predicts $\bar{η}$ increased 0.17 ± 0.22[2σ] cm a -1 between 1983 and 2020 (6.4 ± cm total). The rate of increase nearly doubled to 0.29 ± 0.31 cm a -1 between 1990 and 2020, during which time local field measurements from Circumpolar Active Layer Monitoring sites indicate the active layer increased 0.31 ± 0.22 cm a -1 . The predicted rate of increase more than doubled again between 2002 and 2020, outpacing a near doubling of observed active layer thickening, consistent with trends in terrestrial water storage inferred from Gravity Recovery and Climate Experiment satellite gravimetry and Modern-Era Retrospective Analysis for Research and Applications climate reanalysis. Overall, hydrologic change is accelerating in the Kuparuk River basin, and we provide a theoretical framework for estimating basin-scale changes in active layer water storage from streamflow measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Conservative solute transport processes and associated transient storage mechanisms: Comparing streams with contrasting channel morphologies, land use and land cover

Abstract Land use within a watershed impacts stream channel morphology and hydrology and, therefore, in‐stream solute transport processes and associated transient storage mechanisms. This study evaluated transport processes in two contrasting stream sites where channel morphology was influenced by the surrounding land use, land cover, climate and geologic controls: Como Creek, CO, a relatively undisturbed, high gradient, forested stream with a gravel bed and complex channel morphology, and Clear Creek, IA, an incised, low‐gradient stream with low‐permeability substrate draining an agricultural landscape. We performed conservative stream tracer injections at these sites to address the following questions: (1) How does solute transport vary between streams with differing morphologies? and (2) How does solute transport at each stream site change as a function of discharge? We analysed in‐stream tracer time series data and compared results quantifying solute attenuation in surface and subsurface transient storage zones. Significant trends were observed in these metrics with varying discharge conditions at the forested site but not at the agricultural site. There was a broad range of transport mechanisms and evidence of substantial exchange with both surface and hyporheic transient storage in the relatively undisturbed, forested stream. Changing discharge conditions activated or deactivated different solute transport mechanisms in the forested site and greatly impacted advective travel time. Conversely in the simplified agricultural stream, there was a narrow range of solute transport behaviour across flows and predominantly surface transient storage at all measured discharge conditions. These results demonstrate how channel simplification inhibits available solute transport mechanisms across varying discharge conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

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

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

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

On the Representation of Hyporheic Exchange in Models for Reactive Transport in Stream and River Corridors

Efforts to include more detailed representations of biogeochemical processes in basin-scale water quality simulation tools face the challenge of how to tractably represent mass exchange between the flowing channels of streams and rivers and biogeochemical hotspots in the hyporheic zones. Multiscale models that use relatively coarse representations of the channel network with subgrid models for mass exchange and reactions in the hyporheic zone have started to emerge to address that challenge. Two such multiscale models are considered here, one based on a stochastic Lagrangian travel time representation of advective pumping and one on multirate diffusive exchange. The two models are formally equivalent to well-established integrodifferential representations for transport of non-reacting tracers in steady stream flow, which have been very successful in reproducing stream tracer tests. Despite that equivalence, the two models are based on very different model structures and produce significantly different results in reactive transport. In a simple denitrification example, denitrification is two to three times greater for the advection-based model because the multirate diffusive model has direct connections between the stream channel and transient storage zones and an assumption of mixing in the transient storage zones that prevent oxygen levels from dropping to the point where denitrification can progress uninhibited. By contrast, the advection-based model produces distinct redox zonation, allowing for denitrification to proceed uninhibited on part of the hyporheic flowpaths. These results demonstrate that conservative tracer tests alone are inadequate for constraining representation of mass transfer in models for reactive transport in streams and rivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Data, figures, animations, and scripts associated with the manuscript "Integrated Effects of Site Hydrology and Vegetation on Exchange Fluxes and Carbon Cycling at the Coastal Terrestrial Aquatic Interface"

This package contains the data, figures, animations, and scripts used in "Integrated Effects of Site Hydrology and Vegetation on Exchange Fluxes and Carbon Cycling at the Coastal Terrestrial Aquatic Interface" (Li et al., 2023). This study examines the interactions between soil, vegetation and hydrologic conditions in coastal areas, with a focus on the Chesapeake Bay region. The researchers used a mechanistic model called ATS-PFLOTRAN to explore how these interactions affect exchange process in different regions and the carbon/nitrogen cycle across the terrestrial-aquatic interface (TAI). Simulation scenarios isolate the effects of control factors and reaction constants derived from laboratory experiments. The results show a carbon cycle "hot zone" in coastal wetlands and in the transition zone between wetlands and uplands. Transpiration enhances fluxes between the surface and subsurface domains and increases dissolved oxygen in TAI. The decomposition of leaf-derived organic carbon provides an additional source of carbon for aerobic respiration and denitrification in the TAI. Microbial activity plays a key role in controlling redox conditions and their variability. This modeling study improves the understanding of complex TAI interactions and facilitates the representation of coastal ecosystems in larger-scale Earth system models.Several files can be found from this data package.1. Readme,md: This file describe the Title, Target Journal, Target submission date, Co-author, Science Questions, Hypotheses, Key words, Key message, Model and Data, Repo Structure. The user can read this file first and then go to details.2. mesh.zip: This file contains the mesh file for simulation cases. 3. Figure.zip: This file contains figures used in the manuscript . 4. animation.zip: This file contains animations used in the manuscript . 5. Data.zip: This file contains all required input data, such as concentration and flow boundaries. Also, it contains the DEM and processes results on area fraction. 6. Simulation_setup_and_results.zip: This file contains setup of cases used in the manuscript . The result files are too large, if you need that data, please contact to the Author. 7. notebooks.zip: This file contains the Jupyter notebooks for performing sensitivity analysis and other pre-process and post-process analyses.

54 ENVIRONMENTAL SCIENCES↗