Search NASA⌕ Search

SEARCH · Search NASA

Results for “Total dissolved nitrogen”

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.

83 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Climate and atmospheric deposition drive the inter-annual variability and long-term trend of dissolved organic carbon flux in the conterminous United States

The lateral flux of dissolved organic carbon (DOC) from soils to inland waters and ultimately to the ocean represents a fundamental component of the global carbon cycle. To estimate the DOC flux, we developed an empirical terrestrial-aquatic DOC fluxes model (TAF-DOC). TAF-DOC incorporates various environmental factors (e.g., meteorology, sulfur, and nitrogen deposition) that to-date have not been comprehensively considered or well-represented in existing modeling frameworks. TAF-DOC was applied to estimate spatial-temporal patterns of DOC flux and potential fates across the conterminous United States during the 1985 to 2018 time period. Our results suggest that TAF-DOC successfully characterized spatial-temporal of DOC flux. As expected, the interannual pattern of DOC flux was strongly regulated by precipitation, but the long-term trend was significantly influenced by the rate of atmospheric wet sulfur deposition. From 1985 to 2018, TAF-DOC estimated DOC loading from terrestrial to aquatic ecosystems in the conterminous United States to be 33.5 ± 2.2 TgC per year, which was roughly 0.39–0.49% of total soil organic carbon stock estimates. The dominant fate of terrestrially-derived DOC was delivery to the coastal ocean in riverine export (41%), with another 21% buried in sediment and the remaining 12.8 ± 0.4 TgC per year (38%) returned to the atmosphere through outgassing from inland waters. Finally, assuming the quantities of DOC sediment burial and export to the ocean as an annual sink of terrestrially-derived carbon, budget inventories and models that do not account for DOC flux in the conterminous United States will underestimate the net annual carbon sink by as much as 5.5–6.4%.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Manganese effects on plant residue decomposition and carbon distribution in soil fractions depend on soil nitrogen availability

Recent studies have highlighted the critical role of manganese (Mn) in plant litter decomposition and soil organic carbon (C) cycling in forest ecosystems. Long term nitrogen (N) deposition and N fertilization can increase soil acidity and mobilize bioavailable Mn (Mn 2+ ) in soil. However, no studies have examined the interactive effect of N and Mn fertilization on litter decomposition and carbon distribution in agricultural soils, despite agroecosystems being subject to both N and Mn management. We hypothesized that increased soil N and Mn availability would accelerate plant residue decomposition and transfer of its C to mineral-associated organic matter (MAOM), and that the combined effect of Mn and N enrichment would be greater than the individual effect. Here, we conducted a laboratory incubation experiment by adding 13 C-labeled residue of perennial grass Glyceria striata (Lam.) to agricultural soils that had received 225 kg N ha –1 yr –1 for 27 years (N 1 ) and comparable soils that received no N (N 0 ). Before the experiment, these soils also received three levels of dissolved Mn 2+ , designated M 0 (no additional Mn), M 1 (50 mg kg –1 ), or M 2 (250 mg kg –1 ). We measured total CO 2 production as well as distribution of 13 C from the residue into CO 2 , particulate organic matter (POM), MAOM, and dissolved organic carbon (DOC) over a 1-year period. Manganese amendments significantly increased CO 2 production from residue decomposition in the N 1 soil, but no such effect was observed in the N 0 soil. Manganese also accelerated the loss of residue-derived C from POM and DOC, but increased its recovery in MAOM. However, the positive effect of added Mn in decomposition and recovery in MAOM in the presence of N fertilization occurred only during the initial 30-day decomposition period, where M 2 showed a 12% increase in cumulative CO 2 production from residue, 8% increase in POM loss, and 43% increase in recovery of residue C in MAOM compared to M 0 . For M 1 , only CO 2 emission from residue was significantly higher than Mo during this period. At 365 days M 2 showed 8% increase in CO 2 production, 1% increase in POM loss, and 16% increase in recovery of residue C in MAOM compared to M 0 , but none of these were statistically significant (p < 0.05). This study adds to the growing evidence that increasing Mn availability enhances plant litter decomposition. However, the occurrence and magnitude of Mn-induced stimulation of decomposition is context specific. Further investigation with greater temporal resolution, involving a multitude of litter and soil types and including microbial compositional and functional characterization, is recommended to fully elucidate the interactive role of Mn and N on C cycling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Degradation of Poly- and Perfluoroalkyl Substances (PFAS) in Water via High Power, Energy-Efficient Electron Beam Accelerator

The goal of the 2-year workplan was to see if electron beam (EB) could be used to break down a sub-set of the larger chemical family of per and polyfluoroalkylated substances (PFAS) in an energy efficient and economical manner when compared to conventional water treatment technologies. Year one (Y1) work focused on sample EB treatment work in the Fermi National Accelerator Laboratory’s (FNALs) Accelerator Applications Demonstration and Development (A2D2) EB accelerator. While there are reportedly thousands of types of PFAS, for the point of most of the work herein, a small subset was examined, typically perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA). PFOA and PFOS are two of the most well studied PFAS and are studied for baseline evaluations and are considered most useful. The work from Y1 provided information about the optimal operating parameters and additives to use when treating PFOS and PFOA via EB. The data were then used to see where in a water treatment system an EB accelerator would be best suited to treat PFAS. A conventional water treatment technology, GAC, was then compared to e-beam treatment technology with respect to energy and costs for treatment. In year two (Y2), several conventional e-beam accelerator designs, and FNAL’s developmental compact SRF accelerator design, were evaluated for their suitability in PFAS treatment, from an energy efficiency and cost standpoint. Several EB parameters were evaluated and optimized for the removal of PFOA and PFOS from water at normal pressure and temperature, measured as total PFAS removal. Under the optimized test conditions both PFOA showed complete destruction to inorganic fluoride, and PFOS to inorganic fluoride and sulfate, with mass balance. The effect on PFAS removal relative to solution pH, total EB dose, EB dose rate, dissolved oxygen concentration (DO), temperature, and initial PFAS concentration were evaluated. In general, PFOA was easier to destroy than PFOS. Degradation products, typically observed under less-than-optimal EB conditions, provided insight to degradation mechanisms. Products were identified to rule out possible deleterious biproduct formation. The water radiolysis radical reaction kinetics with PFOS and PFOA were not dependent on the initial concentration over 5-orders of magnitude from 2 μg/L to 20 mg/L. This is thought to be because there was an overabundance of the reactive water radiolysis radicals relative to PFAS molecules and largely attributed to aqueous electrons. The reaction rates appeared to be diffusion limited. Testing at higher concentrations (100-200 mg/L) showed a decrease in removal efficiency, suggesting alternative kinetics, possibly second order rates, at higher concentrations. In all, we successfully defined a set of optimal EB parameters to treat PFOA and PFOS at concentrations of 20 mg/L in water with destruction efficiencies near 100%. We further tested the optimized EB parameters with other types of PFAS, including shorter and longer fluorocarbon chain homologs of PFOA and PFOS, and PFAS with alternative functional groups such as sulfonamides. Based on our results EB can be optimized as an effective destructive technology for removing PFAS from water. The conditions optimized for PFOA and PFOS were less effective with ultra-short fluorocarbon compounds like TFMS, PFES, PFPS and PFBS, and likely require re-optimization of parameters to them. In all, it was determined that from a cost and energy efficiency standpoint, EB would be best applied to waste streams with relatively high concentrations of PFOS and PFOA and is not as cost effective as GAC treatment for removing low concentrations of PFAS from water. Higher concentrations of PFAS can be found in the wastewater of conventional treatment processes such as RO and IE and therefore EB may be used to supplement such treatment technologies. Some real-world IE regeneration wash water and RO reject water containing higher concentrations of PFAS and obtained from pilot scale industrial wastewater treatment system at a fluorochemical manufacturing facility, showed that EB could remove PFAS from such types of wastewaters. The IE regenerant wash water appeared to be the most efficient of the two types of wastewaters tested. However, some further optimization of the EB parameters for the specific PFAS types present in those wastewaters may be required. Also, the effects of co-present TOC and mineral salts should be considered during such optimization efforts. From the experimental Y1 results it was seen that the aqueous electron drives degradation of the PFAS. In a hypothetical water treatment skid using EB for PFAS destruction the parameters of the system should be optimized to promote aqueous electron production. Before EB treatment, the PFAS should be preconcentrated when possible, the pH should be raised to pH 10 or higher to enhance aqueous electron production, and the water should be nitrogen purged to remove dissolved oxygen to minimize aqueous electron scavenging. An excel spreadsheet was created that calculates optimal conditions based on inlet PFAS concentration and desired outlet concentration, by optimizing the accelerator power, dose rate, water treatment rate, pH and dissolved oxygen levels to reach the desired endpoint. Given this information on accelerator operating conditions five different EB accelerator systems were compared. One EB system was a continuous-wave, linear superconducting accelerator being designed at Fermilab. Three other EB systems (IMPELA at 5% and 25% duty factor and the ILU-14) were normal conducting pulsed linear accelerators. The fifth system was an IBA Rhodotron which is a normal conducting, circular, continuous-wave accelerator. The accelerator efficiency (% of the incoming power that is used in water treatment) was the dominating factor in accelerator choice. The radio frequency (RF) power supply and the accelerator design (superconducting versus warm technology) drive the accelerator efficiency. The IBA Rhodotron was seen to be the most energy efficient commercially available technology with a wall-plug (total) power efficiency of 43% at 400 kW. The Fermilab design, with a prototype for a different application currently being fabricated, was the most energy efficient at 55% when driven by a Klystron RF power supply and as high as 77% when powered by a magnetron. As the Fermilab design was the most energy efficient by approximately 10-30%, further design work was done on the accelerator and beam delivery system specific to the destruction of PFAS in water. The Fermilab design is unique from industrial accelerators in that is superconducting. Superconducting technology allows for the acceleration of electrons without losses. The accelerator must be cooled to below the point where it is superconducting and is operated around 4 degrees Kelvin. The bulk of the design work for the accelerator is on making the accelerator as energy efficient as possible so that it does not require liquid helium and can be cooled with conduction cooling via cryocoolers. Final design work resulted in an EB accelerator that would operate at minimally 200 kW and 10 MeV. Prototype construction would cost $\$ $7.8 million dollars when driven by a Klystron power supply. A second version of the same accelerator would cost $\$ $5.5 million dollars when driven by a magnetron that is still under development. The commercially available 300 kW IBA Rhodotron cost was estimated at approximately $\$ $9 million. While it is hard to directly compare, an operational GAC system used by 3M for groundwater treatment capital cost (2022 dollars) was estimated to cost $\$ $3.3 million. While the capital expense of the EB accelerator systems was higher than GAC, the accelerator EB treatment would result in destruction of the PFAS and not just sequestration of PFAS to form a new waste stream that requires further treatment or disposal. The operating cost to destroy the PFAS via 400 kw EB system was less than $\$ $1000/kg of PFAS destroyed when treating at a 20 mg/L PFAS concentration, compared to GAC with operating costs that calculated at $\$ $27,530 per kg of PFAS sequestered when treating 100 μg/L PFOA and PFOS combined concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Effects of nitrogen fertilization and bioenergy crop type on topsoil organic carbon and total Nitrogen contents in middle Tennessee USA

Nitrogen (N) fertilization affects bioenergy crop growth and productivity and consequently carbon (C) and N contents in soil, it however remains unclear whether N fertilization and crop type individually or interactively influence soil organic carbon (SOC) and total N (TN). In a three-year long fertilization experiment in switchgrass (SG: Panicum virgatum L.) and gamagrass (GG: Tripsacum dactyloides L.) croplands in Middle Tennessee USA, soil samples (0–15cm) were collected in plots with no N input (NN), low N input (LN: 84 kg N ha -1 yr -1 in urea) and high N input (HN: 168 kg N ha -1 yr -1 in urea). Besides SOC and TN, the aboveground plant biomass was also quantified. In addition to a summary of published root morphology data based on a separated mesocosm experiment, the root leachable dissolved organic matter (DOM) of both crops was also measured using archived samples. Results showed no significant interaction of N fertilization and crop type on SOC, TN or plant aboveground biomass (ABG). Relative to NN, HN (not LN) significantly increased SOC and TN in both crops. Though SG showed a 15–68% significantly higher ABG than GG, GG showed a 9.3–12% significantly higher SOC and TN than SG. The positive linear relationships of SOC or TN with ABG were identified for SG. However, GG showed structurally more complex and less readily decomposed root DOM, a larger root volume, total root length and surface area than SG. Collectively, these suggested that intensive N fertilization could increase C and N stocks in bioenergy cropland soils but these effects may be more likely mediated by the aboveground biomass in SG and root chemistry and morphology in GG. Future studies are expected to examine the root characteristics in different bioenergy croplands under the field fertilization experiment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Spatial Microbial Respiration Variations in the Hyporheic Zones Within the Columbia River Basin

Abstract While the hyporheic zone (HZ) accounts for a significant portion of whole stream CO 2 concentrations, HZ respiration modeling studies are lacking in quantifying their contributions to the total CO 2 at large watershed/basin scales. Quantifying the contribution of anaerobic respiration is also underappreciated. This study used a carbon‐nitrogen‐coupled river corridor model to quantify HZ aerobic and anaerobic respiration and determined the key factors controlling their spatial variability within the Columbia River Basin (CRB). The modeled respiration patterns showed high spatial variability. Among the nine sub‐basins composing the CRB, the Lower Columbia and the Willamette, which receive higher precipitation, had higher respiration. Medium‐sized rivers (fourth to sixth orders) produced the highest aerobic and anaerobic respiration among reaches of different sizes. At the basin scale, aerobic respiration is dominant, representing approximately 98.7% of the total respiration across the CRB. While most of the reaches were dominant with aerobic respiration, reaches in agricultural land showed a relatively higher anaerobic respiration (18%) ratio. A variable importance analysis showed that hyporheic exchange flux controlled most of the spatial variability of HZ respiration, dominating over other physical variables such as residence time, stream dissolved organic carbon (DOC), nitrate, and dissolved oxygen (DO). The influence of substrate concentration (DOC and DO) is larger in modeling anaerobic respiration than aerobic respiration. Future efforts will focus on improving the estimation of the HZ exchange flux and the implementation of spatially explicit parameterizations for the reactions of interest to reduce model uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Data from: Leaching losses of dissolved organic carbon and nitrogen from agricultural soils in the upper US Midwest

Leaching losses of dissolved organic carbon (DOC) and nitrogen (DON) from agricultural systems are important to water quality and carbon and nutrient balances but are rarely reported; the few available studies suggest linkages to litter production (DOC) and nitrogen fertilization (DON). In this study we examine the leaching of DOC, DON, NO3-, and NH4+ from no-till corn (maize) and perennial bioenergy crops (switchgrass, miscanthus, native grasses, restored prairie, and poplar) grown between 2009 and 2016 in a replicated field experiment in the upper Midwest U.S. Leaching was estimated from concentrations in soil water and modeled drainage (percolation) rates. DOC leaching rates (kg ha-1 yr-1) and volume-weighted mean concentrations (mg L-1) among cropping systems averaged 15.4 and 4.6, respectively; N fertilization had no effect and poplar lost the most DOC (21.8 and 6.9, respectively). DON leaching rates (kg ha-1 yr-1) and volume-weighted mean concentrations (mg L-1) under corn (the most heavily N-fertilized crop) averaged 4.5 and 1.0, respectively, which was higher than perennial grasses (mean: 1.5 and 0.5, respectively) and poplar (1.6 and 0.5, respectively). NO3- comprised the majority of total N leaching in all systems (59-92%). Average NO3- leaching (kg N ha-1 yr-1) under corn (35.3) was higher than perennial grasses (5.9) and poplar (7.2). NH4+ concentrations in soil water from all cropping systems were relatively low (<0.07 mg N L-1). Perennial crops leached more NO3- in the first few years after planting, and markedly less after. Among the fertilized crops, the leached N represented 14-38% of the added N over the study period; poplar lost the greatest proportion (38%) and corn was intermediate (23%). Requiring only one third or less of the N fertilization compared to corn, perennial bioenergy crops can substantially reduce N leaching and consequent movement into aquifers and surface waters.

09 BIOMASS FUELS↗
DOE OSTI2022

Insights into Controls of Mineralogy and Pore Structure on the Density of Methane Adsorption Phase in Shales under Supercritical Conditions

With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.

58 GEOSCIENCES↗
DOE OSTI2022

Assessing the Potential Impact of River Chemistry on Arctic Coastal Production

The Arctic coastal margin receives a disproportionately large fraction of the global river discharge. The bio-geochemistry of the river water as it empties into the marine environment reflects inputs and processes that occur as the water travels from its headwaters. Climate-induced changes to Arctic vegetation and permafrost melt may impact river chemistry. Understanding the impact of river nutrients on coastal marine production, and how this may change in the future, are important for resource managers and community members who monitor and rely on coastal food resources. Using the Energy Exascale Earth System Model we explore the impact of timing and river nutrient concentrations on primary production in each coastal Arctic region and then assess how this influences secondary production and particle fluxes supporting the benthic food web. Our results indicate that while the concentration of Arctic river nitrogen can have a significant impact on annual average nitrogen and primary production in the coastal Arctic, with production increases of up to 20% in the river influenced interior Seas, the timing of the river nutrient inputs into the marine environment appears less important. Bloom timing and partitioning between small and large phytoplankton were minimally impacted by both river nutrient concentration and timing, suggesting that in general, coastal Arctic ecosystem dynamics will continue to be primarily driven by light availability, rather than nutrients. Under a doubling river nutrient scenario, the percentage increase in the POC flux to the benthos on river influenced Arctic coastal shelves was 2-4 times the percentage increase in primary production, suggesting changes to the river nutrient concentration has the potential to modify the Arctic food web structure and dynamics. Generally, the nutrient-induced changes to primary production were smaller than changes previously simulated in response to ice reduction and temperature increase. However, in the Laptev Sea, the production increase resulting from a doubling of river nutrients exceeded the production increase simulated with an atmospheric warming scenario. Dissolved organic carbon is presently poorly represented in the model so its impact on production is hard to simulate. Applying established relationships between modeled DOC, total DOC, and light absorption we illustrate that DOC could play a very important role in modulating production. Our findings highlight the importance of developing more realistic river nutrient and discharge forcing for Earth System Models such that their impact on the critical Arctic coastal domain can be more adequately resolved.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Comparison of Extended and Accelerated VOG Tests

The reprocessing of used nuclear fuel would release volatile radionuclides into the off-gas streams of a processing plant, including ¹²⁹I. The dissolver off-gas and more dilute vessel off-gas streams (VOG) must both be targeted when mitigating ¹²⁹I environmental release because each contains some amount of iodine. Iodine-129 will likely be found as elemental iodine (I₂) and methyl iodide (CH₃I) in these off-gas streams. Reduced silver-exchanged mordenite (AgZ or Ag⁰Z) has been investigated as a potential sorbent for iodine abatement and is studied here under prototypic VOG conditions. Because of the relatively low iodine abundance and high flow rates of the VOG, total iodine concentrations are expected to be in the parts per billion range. Thus, VOG experiments need to run for extended durations at low concentrations to reach sorbent saturation. Because the sorbent will be exposed to oxidizing gas streams for extended periods of time, aging effects and sorbent degradation need to be considered when designing a sorbent-based abatement system. Three sets of experiments were performed with the aim of determining how the adsorption of iodine by AgZ is affected by differing test durations and gas compositions. The first tested the capacity of sorbent aged for 8 months under a humid air stream. The second tested differences in sorbent capacity and mass transfer zone (MTZ) length during high-concentration (1200 ppb v ) CH₃I loading in a humid nitrogen gas (N₂) stream and a humid air gas stream over 28 days. The third tested sorbent capacity and MTZ length during low concentration (< 200 ppb v ) CH₃I and I₂ loading in a humid air stream over 9 months. The results of these tests suggest that AgZ capacity drops by ~50% after 1 month of aging, by ~60% after 2 months of aging, and then does not continue to decrease significantly at up to 8 months of aging. One-month aging tests that compared N₂ and air as the gas stream diluent resulted in similar maximum loading capacities and overall loading curves, indicating that the effects of aging cannot be mitigated by avoiding air as the balance gas. The 9-month extended VOG tests did not result in clear sorbent saturation at the inlet, but it can be inferred using an assumption of a maximum capacity of 45 mg I/g AgZ. Applying this assumption, then the MTZ of CH₃I is 12.8 cm, and the MTZ of I₂ is 12.4 cm. These results are similar to previous tests. Comparisons to tests completed at Idaho National Laboratory suggest that the CH₃I MTZ may be dependent on concentration. Future work should re-evaluate the design of the VOG iodine capture system with this updated data and should seek to understand fundamental characteristics of CH₃I adsorption by AgZ, such as the effect of concentration, gas velocity, and gas composition.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Model Inputs, Outputs, and Scripts associated with: “Combined effects of stream hydrology and land use on basin-scale hyporheic zone denitrification in the Columbia River Basin”

This data package is associated with the publication “Combined effects of stream hydrology and land use on basin‐scale hyporheic zone denitrification in the Columbia River Basin”, published in Water Resource Research (Son et al.2022) available at https://doi.org/10.1029/2021WR031131. 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 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) denitrification at the NHDPLUS stream reach scales. The RCM used in this study combines empirical substrate models derived from observations and three microbially driven reactions, including two-step denitrification and aerobic respiration, are considered within the HZ. 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 study uses the RCM to explore the spatial patterns of HZ denitrification across reaches with different sizes and land use in the CRB. Our main objective is to use the RCM as a virtual reality model, and the machine-learning models as surrogates that encapsulate the complexities of the physics-based model while identifying the importance of different variables that are not evident in the model conceptualization. We do not include a direct comparison of the modeled HZ denitrification and measurements; however, the RCM can capture the overall spatial patterns of the HZ denitrification because the model inputs and its reaction networks are based on well-established theory and a physical-based model. The combination of the model-based predictions and a machine-learning approach (e.g., random forest) is used to improve our understanding of what variables of the model are associated with spatial patterns of the modeled denitrification across reaches with different sizes and land uses, and to develop a proxy model using measurable variables to reproduce the simulated patterns.This dataset contains five folders: (1) model_inputs, (2) model_outputs, (3) Rscripts, (4) figures, and (5) model_codes. It also contains a readme, file level metadata (FLMD), and data dictionary (dd). Please see the FLMD for a list of all the files contained in this data package and descriptions for each. The model_inputs folder contains the model inputs used to drive the model simulations. The model_outputs folder contains key model output files from the river corridor model. The Rscripts folder contains the Rscripts for pre- and post- processing model results. The figures folder contains the raw figures associated with the manuscript. The model_codes folder includes key model source codes/input files. All files are .jpg, .jpeg, .out, .e, .od, .dat, .sub, .F90, .0, .R, .sbx, .cpg, .sbn, .shx, .shp, .dbf, .prj, .tfw, .tif, .xml, .pdf, or .csv.

54 ENVIRONMENTAL SCIENCES↗
DOE 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↗