Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hyporheic”

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 37 records · Page 2

Hyporheic Exchange Due to Cobbles on Sandy Beds

Hyporheic exchange induced by riverbed topography and roughness provides important ecosystem functions. We investigated how cobble clusters embedded in a finer-grained sand bed affect the near-bed channel flow and the exchange of surface and subsurface water. We tested how the spacing and embeddedness of cobbles altered hyporheic exchange through a three-dimensional fully coupled surface-subsurface model. The 3D modeling framework captured the full physics of the exchange process locating the lateral position of upwelling zones on the side of the cobble. As the cobbles protrude more into the channel, eddies appear downstream, with smaller vortices where cobbles are closer together. In our simulations, hyporheic exchange increased with the spacing and protrusion ratio. The travel time in the hyporheic zone also increased with the protrusion ratio because the hyporheic flow paths are longer when the cobble is more protruded. With different spacings at the highest protrusion ratio, the average slope of early time breakthrough curve was steeper for large spacing settings: the travel time increases as the spacing decreases.

54 ENVIRONMENTAL SCIENCES↗

Hyporheic‐Zone Processes and Stream Oxygen Dynamics: Insights From a Multiscale Reactive Transport Model

Aquatic ecosystem metabolism encapsulates the daily fixation (gross primary production, GPP d ) and mineralization (ecosystem respiration, ER d ) of organic carbon. In fluvial systems, these are commonly estimated by inverse solutions to field observations using a model that describes oxygen concentrations varying in the water column in response to metabolic fluxes and air‐water gas exchange controlled by a rate coefficient (K 600 ). The most common conceptual model is the single‐station metabolism (SSM) model. The simplicity and flexibility of this conceptualization make it attractive; however, it implicitly assumes that all the processes that consume oxygen in fluvial systems can be lumped into a bulk estimate of respiration with poorly understood consequences for estimates of GPP d , ER d , and K 600 . Here, we focus on the implications of using SSM conceptualization when estimating metabolic fluxes from oxygen dynamics in channels where hyporheic exchange occurs. We use a new multiscale numerical model for reactive transport in streams that represents hyporheic exchange and streambed heterotrophic respiration. Nondimensionalization of this model reveals dimensionless groups that collectively control oxygen dynamics. Numerical experiments offer a mechanistic understanding of the impacts of hyporheic exchange on diel oxygen dynamics revealing that potential biases arise from neglecting mass transfer limitations. Specifically, we found that hyporheic exchange significantly affects diel oxygen dynamics, even for nonreactive streambed sediments. Moreover, while the SSM performs well in many situations, we find conditions where significant bias is produced by hyporheic exchange, even when oxygen data are well‐fitted. These situations pose a major challenge in the interpretation of metabolism assessment estimates.

Gomez‐Velez, Jesus D. [Oak Ridge National Laborato↗

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗

Bed form-induced hyporheic exchange and geochemical hotspots

Small-scale bed form topographies control hyporheic exchange and biogeochemical processes within aquatic sediments, which ultimately affect water quality and nutrient cycling at the watershed scale. The impact of three-dimensional and small-scale bed form topographies on hyporheic exchange and solute mixing is investigated in the present work. The effect of bed form morphologies on the development of zones of enhanced reaction rates (i.e., hotspots) is also studied. A computational fluid dynamics model to simulate river flow over bed forms is combined with a subsurface flow and multicomponent reactive solute transport model. A wide variety of bed form topographies are generated using geometric models by varying parameters controlling curvature as well as bed form wavelength and amplitude. The results in this research suggest that out-of-phase bed forms generate more complex hyporheic flow patterns which reduce the efficiency of solute transformations. Higher phase shifts in bed form shapes result in overall higher average velocity, larger zones of enhanced pressure and reaction rates, and higher amounts of solute exchange. Moreover, the bed form shapes control the reaction process for a wide range of sediment conductivities. Here this study advances the understanding of the effects of complex and small scale morphological features on hyporheic exchange processes including the rate and spatio-temporal distribution of reaction hotspots.

54 ENVIRONMENTAL SCIENCES↗

Release Timing and Duration Control the Fate of Photolytic Compounds in Stream-Hyporheic Systems

Predicting environmental fate requires an understanding of the underlying, spatiotemporally variable interaction of transport and transformation processes. Photolytic compounds, for example, interact with both time-variable photolysis and the perennially dark hyporheic zone, generating potentially unexpected dynamics that arise from time-variable reactivity. This interaction has been found to significantly impact environmental fate but is commonly oversimplified in predictive models. Our primary objective was to explore how time-variable photolysis and hyporheic storage interact across a range of photolysis rates to control the fate and transport of photolytic solutes in stream-hyporheic systems. In this study, we released a photolytic compound in a natural system at different times of day and simulated variable release timing and durations of photolytic compounds spanning half-lives ranging from about 2 min to 900 hr. To contextualize these results, we simulated five photolytic compounds with varied half-lives to systematically study the interaction between release timing, duration, and reactivity. We found that the environmental fate and transport of photolytic compounds are highly variable as a function of release timing, which controls when, where, and for how long solute is stored in the hyporheic zone or exposed to in-channel photolysis. Finally, this knowledge can be used to improve predictions for photolytic compounds or assess potential impacts for an anticipated discharge or treatment.

54 ENVIRONMENTAL SCIENCES↗

Simultaneous attenuation of trace organics and change in organic matter composition in the hyporheic zone of urban streams

Trace organic compounds (TrOCs) enter rivers with discharge of treated wastewater. These effluents can contain high loads of dissolved organic matter (DOM). In a 48 h field study, we investigated changes in molecular composition of seven DOM compound classes (FTICR-MS) and attenuation of 17 polar TrOCs in a small urban stream receiving treated wastewater. Correlations between TrOCs and DOM were used to identify simultaneous changes in surface water and the hyporheic zone. Changes in TrOC concentrations in surface water ranged between a decrease of 29.2% for methylbenzotriazole and an increase of 152.2% for the transformation product gabapentin-lactam. In the hyporheic zone, only decreasing TrOC concentrations were observed, ranging from 4.9% for primidone to 93.8% for venlafaxine . TrOC attenuation coincided with a decline of molecular diversity of easily biodegradable DOM compound classes while molecular diversity of poorly biodegradable DOM compound classes increased. This concurrence indicates similar or linked attenuation pathways for biodegradable DOM and TrOCs. Strong correlations between TrOCs and DOM compound classes as well as high attenuation of TrOCs primarily occurred in the hyporheic zone. This suggests high potential for DOM turnover and TrOC mitigation in rivers if hyporheic exchange is sufficient.

54 ENVIRONMENTAL SCIENCES↗

Maximum respiration rates in hyporheic zone sediments are primarily constrained by organic carbon concentration and secondarily by organic matter chemistry

Abstract. River corridors are fundamental components of the Earth system, and their biogeochemistry can be heavily influenced by processes in subsurface zones immediately below the riverbed, referred to as the hyporheic zone. Within the hyporheic zone, organic matter (OM) fuels microbial respiration, and OM chemistry heavily influences aerobic and anaerobic biogeochemical processes. The link between OM chemistry and respiration has been hypothesized to be mediated by OM molecular diversity, whereby respiration is predicted to decrease with increasing diversity. Here we test the specific prediction that aerobic respiration rates will decrease with increases in the number of unique organic molecules (i.e., OM molecular richness, as a measure of diversity). We use publicly available data across the United States from crowdsourced samples taken by the Worldwide Hydrobiogeochemical Observation Network for Dynamic River Systems (WHONDRS) consortium. Our continental-scale analyses rejected the hypothesis of a direct limitation of respiration by OM molecular richness. In turn, we found that organic carbon (OC) concentration imposes a primary constraint over hyporheic zone respiration, with additional potential influences of OM richness. We specifically observed respiration rates to decrease nonlinearly with the ratio of OM richness to OC concentration. This relationship took the form of a constraint space with respiration rates in most systems falling below the constraint boundary. A similar, but slightly weaker, constraint boundary was observed when relating respiration rate to the inverse of OC concentration. These results indicate that maximum respiration rates may be governed primarily by OC concentration, with secondary influences from OM richness. Our results also show that other variables often suppress respiration rates below the maximum associated with the richness-to-concentration ratio. An important focus of future research will identify physical (e.g., sediment grain size), chemical (e.g., nutrient concentrations), and/or biological (e.g., microbial biomass) factors that suppress hyporheic zone respiration below the constraint boundaries observed here.

58 GEOSCIENCES↗

Combined Effects of Stream Hydrology and Land Use on Basin‐Scale Hyporheic Zone Denitrification in the Columbia River Basin

Abstract Denitrification in the hyporheic zone (HZ) of river corridors is crucial to removing excess nitrogen in rivers from anthropogenic activities. However, previous modeling studies of the effectiveness of river corridors in removing excess nitrogen via denitrification were often limited to the reach‐scale and low‐order stream watersheds. We developed a basin‐scale river corridor model for the Columbia River Basin with random forest models to identify the dominant factors associated with the spatial variation of HZ denitrification. Our modeling results suggest that the combined effects of hydrologic variability in reaches and substrate availability influenced by land use are associated with the spatial variability of modeled HZ denitrification at the basin scale. Hyporheic exchange flux can explain most of spatial variation of denitrification amounts in reaches of different sizes, while among the reaches affected by different land uses, the combination of hyporheic exchange flux and stream dissolved organic carbon (DOC) concentration can explain the denitrification differences. Also, we can generalize that the most influential watershed and channel variables controlling denitrification variation are channel morphology parameters (median grain size (D50), stream slope), climate (annual precipitation and evapotranspiration), and stream DOC‐related parameters (percent of shrub area). The modeling framework in our study can serve as a valuable tool to identify the limiting factors in removing excess nitrogen pollution in large river basins where direct measurement is often infeasible.

54 ENVIRONMENTAL SCIENCES↗

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↗

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↗

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↗

Testing Hidden Assumptions of Representativeness in Reach-Scale Studies of Hyporheic Exchange

Field studies of hyporheic exchange in mountain systems are often conducted using short study reaches and a limited number of observations. It is common practice to assume these study reaches represent hyporheic exchange at larger scales or different sites and to infer general relationships among potential causal mechanisms from the limited number of observations. However, these assumptions of representativeness are rarely tested. In this study, we develop numerical models from four segments of mountain streams in different geomorphologic settings and extract shorter reaches to test how representative exchange metrics are in shorter reaches compared to their reference segments. We also map the locations of the representative reaches to determine if a pattern exists based on location. Finally, we compare variance of these shorter within-site reaches to 29 additional reaches across the same basin to understand the impacts of inferring causal mechanisms, for example, the expectation that wide and narrow valley bottoms will yield different hyporheic exchange patterns. Our results show that the location and length strategy of the study reach must be considered before assuming an exchange metric to be representative of anything other than the exact segment studied. Further, it is necessary to quantify within and between site variations before making causal inferences based on observable characteristics, such as valley width or stream morphology. Our findings have implications for future field practices and how those practices are translated into models.

54 ENVIRONMENTAL SCIENCES↗

Toward a Universal Model of Hyporheic Exchange and Nutrient Cycling in Streams

In this paper we demonstrate that several ubiquitous hyporheic exchange mechanisms can be represented simply as a one-dimensional diffusion process, where the diffusivity decays exponentially with depth into the streambed. Based on a meta-analysis of 106 previously published laboratory measurements of hyporheic exchange (capturing a range of bed morphologies, hydraulic conditions, streambed properties, and experimental approaches) we find that the reference diffusivity and mixing length-scale are functions of the permeability Reynolds Number and Schmidt Number. These dimensionless numbers, in turn, can be estimated for a particular stream from the median grain size of the streambed and the stream's depth, slope, and temperature. Application of these results to a seminal study of nitrate removal in 72 headwater streams across the United States, reveals: (a) streams draining urban and agricultural landscapes have a diminished capacity for in-stream and in-bed mixing along with smaller subsurface storage zones compared to streams draining reference landscapes; (b) under steady-state conditions nitrate uptake in the streambed is primarily biologically controlled; and (c) median reaction timescales for nitrate removal in the hyporheic zone are 0.5 and 20 hr for uptake by assimilation and denitrification, respectively. While further research is needed, the simplicity and extensibility of the framework described here should facilitate cross-disciplinary discussions and inform reach-scale studies of pollutant fate and transport and their scale-up to watersheds and beyond.

54 ENVIRONMENTAL SCIENCES↗

Production of hydrogen peroxide in an intra-meander hyporheic zone at East River, Colorado.

The traditionally held assumption that photo-dependent processes are the predominant source of H 2 O 2 in natural waters has been recently questioned by an increasing body of evidence showing the ubiquitiousness of H 2 O 2 in dark water bodies and in groundwater. In this study, we conducted field measurement of H 2 O 2 in an intra-meander hyporheic zone and in surface water at East River, CO. On-site detection using a sensitive chemiluminescence method suggests H 2 O 2 concentrations in groundwater ranging from 6 nM (at the most reduced region) to ~ 80 nM (in a locally oxygen-rich area) along the intra-meander transect with a maxima of 186 nM detected in the surface water in an early afternoon, lagging the maximum solar irradiance by ~1.5 h. Our results suggest that the dark profile of H 2 O 2 in the hyporheic zone is closely correlated to local redox gradients, indicating that interactions between various redox sensitive elements could play an essential role. Due to its transient nature, the widespread presence of H 2 O 2 in the hyporheic zone indicates the existence of a sustained balance between H 2 O 2 production and consumption, which potentially involves a relatively rapid succession of various biogeochemically important processes (such as organic matter turnover, metal cycling and contaminant mobilization). More importantly, this study confirmed the occurrence of reactive oxygen species at a subsurface redox transition zone and further support our understanding of redox boundaries on reactive oxygen species generation and as key locations of biogeochemical activity.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Continental-Scale Controls on Hyporheic Respiration Revealed by Knowledge-Guided Machine Learning

Hyporheic zone sediments regulate organic matter turnover and in-stream respiration, yet controls on sediment respiration remain poorly constrained across heterogeneous river networks, limiting prediction of stream metabolism and carbon processing at continental scales. Here, we integrate observations from ~90 river corridors across the United States in the WHONDRS consortium with a knowledge-guided machine learning (KGML) framework that couples thermodynamic rate theory with machine learning to identify dominant controls on hyporheic respiration. Diagnostic analyses show that organic matter concentration and thermodynamic favorability define an upper bound on respiration potential, whereas biological catalytic capacity and physical accessibility jointly govern realized respiration rates through interaction effects. To represent unmeasurable accessibility constraints, we use the mechanistic model as a scaffold for KGML, allowing machine learning to target residual structure not explained by process theory. This hybrid framework improves predictive skill relative to both the mechanistic model alone and fully data-driven models while preserving interpretability. These results indicate that variability in hyporheic respiration is largely mechanistically structured and demonstrate how integrating process theory with explainable AI enhances predictive performance while enabling scalable synthesis of river corridor observations.

Zheng, Jianqiu↗

Dissolved oxygen sensor in an automated hyporheic sampling system reveals biogeochemical dynamics

Many river corridor systems frequently experience rapid variations in river stage height, hydraulic head gradients, and residence times. The integrated hydrology and biogeochemistry of such systems is challenging to study, particularly in their associated hyporheic zones. Here we present an automated system to facilitate 4-dimensional study of dynamic hyporheic zones. It is based on combining real-time in-situ and ex-situ measurements from sensor/sampling locations distributed in 3-dimensions. A novel dissolved oxygen (DO) sensor was integrated into the system during a small scale study. We measured several biogeochemical and hydrologic parameters at three subsurface depths in the riverbed of the Columbia River in Washington State, USA, a dynamic hydropeaked river corridor system. During the study, episodes of significant DO variations (~+/- 4 mg/l) were observed, with minor variation in other parameters (e.g., <~+/-0.15 mg/l NO 3 ). DO concentrations were related to hydraulic head gradients, showing both hysteretic and non-hysteretic relationships with abrupt (hours) transitions between the two types of relationships. The observed relationships provide a number of hypotheses related to the integrated hydrology and biogeochemistry of dynamic hyporheic zones. We suggest that preliminary high-frequency monitoring is advantageous in guiding the design of long term monitoring campaigns. The study also demonstrated the importance of measuring multiple parameters in parallel, where the DO sensor provided the key signal for identifying/detecting transient phenomena.

54 ENVIRONMENTAL SCIENCES↗

River and hyporheic zone water geochemical data from OsmoSamplers, East River Pumphouse, Colorado, Nov 2017-Sept 2018

This data package contains geochemical and river water temperature data from OsmoSamplers deployed at East River, Colorado pump house. There are three data files. The first one contains sulfate (SO4) and chloride (Cl) concentrations from riverwater and hyporheic zone (sediment water interface (SWI), 4 cm below seafloor (cmbsf), and 12 (cmbsf)). The second file contains methane concentrations, methane stable carbon isotope ratios and water isotopes from hyporheic zone porewaters. And the third is from a temperature logger during the deployment. Files are provided as .xls or .csv. This research was performed as part of the SBR funded project “Quantifying Subsurface Biogeochemical Variability in a High Altitude Watershed During Winter Isolation” DE-SC0018328. The novel part of this dataset is that we quantify these geochemical parameters throughout the winter isolation period, which is difficult to measure, as well as within the hyporheic zone. The data files were updated 3/1/2021 because the column label for the water isotopes was incorrect. Also, the csv files were also updated since they were not formatted correctly.

54 ENVIRONMENTAL SCIENCES↗

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↗