Search NASA⌕ Search

SEARCH · Search NASA

Results for “Watershed modeling”

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

Represent precipitation-induced geological hazards in Earth system models using artificial intelligence

Precipitation-induced geological hazards, such as debris flow, landslides, mudflow and rockfalls (hereafter referred to as landslides), pose serious threats to public safety in many areas through the world. As residential properties and infrastructure in the US have increasingly expanded into landslide-prone areas and the drivers of landslides (e.g., wildfires and hurricanes) are predicted to intensify under climate warming, losses and fatalities from landslides are likely to increase in the future. Although our understanding of geoenvironmental factors and mechanisms contributing to landslides has greatly improved, only moderate progress has been made in predicting landslides out of well-studied watersheds. Furthermore, explicitly representing various landslide-related processes in Earth system models (ESMs), from the buckling of local bearing elements in granular materials, to frictional sliding between grains, formation of microcracks in the soil matrix, rupture of capillary bridges, or breakage of plant roots3 is still very unlikely within the next decade, even with the help of exascale computers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

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

Advancing a Watershed Hydro-Biogeochemical Theory: Linking Water Travel Time and Reaction Rates Under Changing Climate

The major goal of the project is to answer the question how and to what degree does a warming climate alter hydrological functioning and (biogeochemical and chemical weathering) reaction rates at the watershed scale? In particular, how does hydrological functioning (water storage, flow paths, and travel time) change responding to a warming climate? What are the corresponding alterations in biogeochemical and chemical weathering rates at the watershed scale? What is the general watershed-scale reaction rate law that causally links reaction rates to metrics of hydrological functioning that are determined by external hydroclimatic conditions and internal watershed structure characteristics? To answer these questions, we have used field measurements and reactive transport modeling. We have found that: 1) stream chemistry has altered tremendously in warmer years; 2) In this water-limited environment, hydrology is the dominant driver for soil respiration and carbon lateral export at the watershed scale; 3) deeper carbon respiration plays an essential role in regulating the patterns of stream carbon in low and high flow conditions, and deeper groundwater will become increasingly important as warming intensifies. This indicates in mountain streams, a warming climate will modify not only the timing and magnitude of streamflow as expected, but also the chemistry and quality of stream water.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A hybrid data–model approach to map soil thickness in mountain hillslopes

Abstract. Soil thickness plays a central role in the interactions between vegetation, soils, and topography, where it controls the retention and release of water, carbon, nitrogen, and metals. However, mapping soil thickness, here defined as the mobile regolith layer, at high spatial resolution remains challenging. Here, we develop a hybrid model that combines a process-based model and empirical relationships to estimate the spatial heterogeneity of soil thickness with fine spatial resolution (0.5 m). We apply this model to two aspects of hillslopes (southwest- and northeast-facing, respectively) in the East River watershed in Colorado. Two independent measurement methods – auger and cone penetrometer – are used to sample soil thickness at 78 locations to calibrate the local value of unconstrained parameters within the hybrid model. Sensitivity analysis using the hybrid model reveals that the diffusion coefficient used in hillslope diffusion modeling has the largest sensitivity among all input parameters. In addition, our results from both sampling and modeling show that, in general, the northeast-facing hillslope has a deeper soil layer than the southwest-facing hillslope. By comparing the soil thickness estimated between a machine-learning approach and this hybrid model, the hybrid model provides higher accuracy and requires less sampling data. Modeling results further reveal that the southwest-facing hillslope has a slightly faster surface soil erosion rate and soil production rate than the northeast-facing hillslope, which suggests that the relatively less dense vegetation cover and drier surface soils on the southwest-facing slopes influence soil properties. With seven parameters in total for calibration, this hybrid model can provide a realistic soil thickness map with a relatively small amount of sampling dataset comparing to machine-learning approach. Integrating process-based modeling and statistical analysis not only provides a thorough understanding of the fundamental mechanisms for soil thickness prediction but also integrates the strengths of both statistical approaches and process-based modeling approaches.

58 GEOSCIENCES↗
DOE OSTI2021

Modeling Terrestrial Dissolved Organic Carbon Loading to Western Arctic Rivers

The mobilization and land-to-ocean transfer of dissolved organic carbon (DOC) in Arctic watersheds is intricately linked with the region's climate and water cycle, and furthermore at risk of changes from climate warming and associated impacts. This study quantifies model-simulated estimates of runoff, surface and active layer leachate DOC concentrations and loadings to western Arctic rivers, specifically for basins that drain into coastal waters between and including the Yukon and Mackenzie Rivers. Model validation leverages data from other field measurements, synthesis studies, and modeling efforts. The simulations effectively quantify DOC leaching in surface and subsurface runoff and broadly capture the seasonal cycle in DOC concentration and mass loadings reported from other studies that use river-based measurements. A marked east-west gradient in simulated spring and summer DOC concentrations of 24 drainage basins on the North Slope of Alaska is captured by the modeling, consistent with independent data derived from river sampling. Simulated loadings for the Mackenzie and Yukon show reasonable agreement with estimates of DOC export for annual totals and four of the six seasonal comparisons. Nearly equivalent loading occurs to rivers which drain north to the Beaufort Sea and west to the Bering and Chukchi Seas. The modeling framework provides a basis for understanding carbon export to coastal waters and for assessing impacts of hydrological cycle intensification and permafrost thaw with ongoing warming in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Role of underlying surface, rainstorm and antecedent wetness condition on flood responses in small and medium sized watersheds in the Yangtze River Delta region, China

Floods are becoming more frequent and are hard to control due to the shortage of water conservancy projects in small and medium sized basins, especially in developing regions. Understanding the hydrologic responses and estimated flood characteristics to storm events can help to predict flood disasters and form better effective mitigation and adaptation strategies. Here, using observations of several representative watersheds and Artificial Neural Networks (ANN) and Principal Component Regression (PCR) models, we conducted data-driven analyses to examine the flood responding characteristics. We quantified the relative contributions of the influencing factors to the variation of each flood characteristic. Statistical analysis of the observations shows that the drainage area plays a key role in determining the distribution of lag time and peak discharge. Rainstorm variability has direct influence on floods, and typhoon-induced rainstorms with high total rainfall and rainfall intensity generate higher lag time, flood peak, unit discharge and runoff depth, but lower runoff coefficient. The ANN and PCR models accurately predicted the variations of flood features using the driving factors including physical geographical characteristics, rainstorm features, and antecedent condition. Physical geographic characteristics are key influential factors of lag time, flood peak and runoff coefficient, while the rainstorm features control the magnitude of unit discharge and runoff depth. These results indicate that floods are mainly affected by rainstorm features and physical geographic characteristics in the Yangtze River delta, and it might become more damaging with the increasing rainfall extremes and sprawling impervious surfaces in a changing environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Diagnosing dominant controls on carbon exchanges in high altitude, western U.S. headwaters. Final report

An observational and modeling research approach was implemented in the East River CO, DOE Watershed Function SFA site. An eddy-covariance flux tower was emplaced located within the basin to observe movement of moisture and carbon. These new observations provided critical information for understanding subsurface controls on terrestrial- aquatic exchanges of carbon and water.

58 GEOSCIENCES↗
DOE OSTI2023

Artificial intelligence models, photos, and data associated with the manuscript “Quantifying Streambed Grain Size, Uncertainty, and Hydrobiogeochemical Parameters Using Machine Learning Model YOLO” (v2)

This data package is associated with the manuscript “Quantifying Streambed Grain Size, Uncertainty, and Hydrobiogeochemical Parameters Using Machine Learning Model YOLO” published in Water Resources Research (Chen et al., 2024). This data package includes the training, validation, testing, and prediction data used by the artificial intelligence (AI) model for automated grain size and hydro-biogeochemistry quantification using streambed photos. The grain size data are extracted for each photo using You Look Only Once (YOLO), a pre-trained object detection model. This data package was originally published in October 2023. It was updated August 2025 (v2; new and modified files). File and folder names were not revised to indicate changes. See the change history section in the readme for more details. Please see flmd.csv for a list of all files contained in this data package and descriptions for each. Please see dd.csv for a data dictionary that defines the column headers of .csv files in the data package. This dataset is comprised of one data folder containing (1) file-level metadata; (2) data dictionary; (3) readme; and (4) six subfolders. Subfolders 1 to 4 include the training, validation, testing, and prediction data. Subfolder 5_Summary includes the summary results of different combinations of training, validation, testing, and prediction data. Subfolder 6_SupplementalData includes additional data downloaded from public sources (Kaufman et al., 2023a; Kaufman et al., 2023b; Garefalakis et al., 2023; Mair et al., 2024; https://github.com/river-corridors-sfa/Geospatial_variables). In total, the data package includes 110 folders and 44,283 files. These files include 9,047 .jpg photos, 1 .png photo, 3 .tif photos; 26,639 photo labels and individual grain sizes and probability from AI (.txt); 8,447 grain size distribution data (.dat); and 126 CSV files for results summary, and 14 required metadata files (.xlsx). The summary CSV files contain 68 columns and approximately 2,200 rows that represent photo names, site locations, recording time, GPS coordinates, grains sizes (D10, D50, D60, and D84), number of grains, and additional hydro-biogeochemical data such as water depth, flow velocity, Manning’s coefficient, friction factor, hydraulic conductivity, permeability, streambed interstitial velocity magnitude, mass transfer rate, and nitrate uptake velocity. The photos were obtained from 75 sites in the Yakima River Basin and the Columbia River shorelines, and other associated data from samples and sensors obtained when the photos were taken are publicly available (Fulton et al. 2022; Grieger et al. 2023). All files are .csv, .txt, .dat, .jpg, or .pdf. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected some of 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 OSTI2025

Steel Creek, Pen Branch, and D-Area Watershed Stream Gauging Stations

A network of stream gauging systems were installed in the Steel Creek, Pen Branch, D-Area Discharge Canal, and the D006 Stream in support of the groundwater modeling efforts for the P-Area Groundwater Operable Unit (OU); Chemical, Metals, and Pesticides (CMP) Pits OU; and the D-Area Watershed, respectively. Each location is monitored by a MACE Floseries3 FloPro data logger and a MACE doppler ultrasonic area/velocity sensor. Each stream gauging system is powered by an internal 12-volt battery supplied by a solar panel with a trickle charger. Information collected by each data logger is logged internally and telecommunicated via a cellular network to an online server for real time analysis and monitoring. The MACE doppler ultrasonic area/velocity sensor can measure stream depth and velocities to give output values of flow rates, total flow, net flow, and volumes. The water depth is measured by a ceramic pressure transducer located on the top of the sensor. The velocity is measured by a continuous wave doppler sensor to give an average velocity across the whole stream profile. This report discusses the equipment and methods used to install continuous stream gauging stations and provides a summary of data collected through FY2025.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data and scripts associated with a manuscript analyzing ELM-FATES parameter sensitivity under pre-fire and postfire scenarios using machine learning

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “Fire Severity-Dependent Shifts in Vegetation Parameter Sensitivity: A Pre- and Post-Fire Analysis Using ELM-FATES and Explainable AI” submitted to Journal of Advances in Modeling Earth Systems (Zahura et al. 2026). The study examines vegetation physiological parameters controlling pre-fire and post-fire vegetation dynamics. To support this analysis, 73 vegetation parameters in Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (Fisher et al., 2018) , which is coupled with E3SM (Energy Exascale Earth System Model) land model (ELM, ELM-FATES), were perturbed using a Sobol sequence to generate 1,024 ensemble members for two plant functional types: needleleaf evergreen extratropical trees (NEET) and C3 grass. Simulations were conducted for the pre-fire period (2016) and post-fire period (2018–2023). Burn severity was represented by modifying the Nesterov index in FATES to 75,000, 150,000, and 300,000 for low, moderate, and high severity, respectively. A no-fire scenario was also included. Simulations were performed for 16 grid cells in the American River Watershed across different burn severities and plant functional types. XGBoost (eXtreme Gradient Boosting) models were trained using the parameter ensembles and ELM-FATES-simulated outputs, including leaf area index (LAI), gross primary productivity (GPP), aboveground biomass, vegetation evaporation, transpiration, and soil evaporation. Models were trained separately for each year and burn severity, followed by SHAP (SHapley Additive exPlanations) analysis to identify changes in dominant parameters after fire disturbance. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package contains the ELM-FATES simulation data. The scripts and data related to the analysis will be added later. The inputs and outputs from ELM-FATES are inside the “FATES” folder. “FATES_domain_surface” contains the domain and surface netcdfs that were used to run ELM-FATES in the study area. “FATES_parameters” contains the 1024 ensembles that were generated using Sobol sequence. “FATES_outputs” folder contains ELM-FATES simulated variables. All files are .csv and .nc (NetCDF).

Aboveground biomass↗
DOE OSTI2023

Estimation of hydraulic conductivity in a watershed using sparse multi-source data via Gaussian process regression and Bayesian experimental design

Enhanced water management systems depend on accurate estimation of subsurface hydraulic properties. However, geologic formations can vary significantly, so information from a single source (e.g., widely spaced boreholes) is insufficient in characterizing subsurface aquifer properties. Therefore, multiple sources of information are needed to complement the hydrogeology understanding of a region. Here, this study presents a numerical framework in which information from different measurement sources is combined to characterize the 3D random field in a multi-fidelity prediction model. Coupled with the model, a Bayesian experimental design was used to determine the best future sampling locations. The Upper Sangamon watershed in east-central Illinois was selected as the case study site, where the multi-fidelity Gaussian process model was used to estimate the hydraulic conductivity in the region of interest. Multi-source observation data were obtained from electrical resistivity and borehole pumping tests. The accuracy of the model prediction is dependent on the locations and the distribution of both high- and low-fidelity data. Furthermore, the multi-fidelity model was compared with the single-fidelity model. The uncertainties and confidence in the measurements and parameter estimates were quantified and used to design future cycles of data collection to further improve the confidence intervals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Post-fire time series photos from five sites across the Oak Creek watershed, Washington

This dataset supports a broader study examining wildfire impacts on hydrologic connectivity across 5 sites within the Oak Creek watershed and the resulting biogeochemical impacts. Sites were selected using the Advanced Terrestrial Simulator (ATS) hydrologic model to identify locations with varying groundwater contributions and hydrologic responses across different burn severity scenarios. The Retreat Fire burned from July 23 to August 2, 2024, affecting all five sites. This dataset provides time series game camera photos, while the broader study includes continuous water quality monitoring, biogeochemical sampling of water and soils, precipitation data, and organic matter analysis. The other data types and additional metadata (include site environmental information) can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3018020. Because this study is ongoing, this data package will be updated regularly to include newly collected photos. 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. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) readme; (5) field protocol; and (6) folders of game camera photos. The game camera photos are organized by site with subfolders by month of collection. The field metadata contains a subset of the information collected that is most relevant to photo-processing. The full set of field metadata can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3018020. All files are .csv, .pdf, or .jpg.

Burn severity↗
DOE OSTI2020

Conceptual Model Update of Mercury Sources and Flux at Y-12 and Upper East Fork Poplar Creek, Oak Ridge,Tennessee

The interactions between the key processes that control the fate and transport of mercury (Hg) in the sediments, groundwater, and stream water within the Y-12 National Security Complex (Y-12) in Oak Ridge, Tennessee are complex, and many aspects are not well understood. In contaminated environments such as Y-12, conceptual models are frequently developed to aid in visualizing and understanding the dynamic nature of the hydrologic, geochemical, and physical environment. These conceptual models integrate data in an internally consistent manner to understand processes that control the fate and transport of contaminants. Over the past few decades of environmental investigation at Y-12, a number of conceptual models have been developed to identify and define various technical processes at various scales. The main purpose of the activities described in this report is to update the 2011 conceptual model for Upper East Fork Poplar Creek using the most recent Hg concentration and flux data and informed scientific interpretation. Emphasis was placed on improving the previous model by implementing a more detailed spatial and temporal approach to visualize transport pathways in the watershed and trends in flux and concentration over time, and to compare baseflow and stormflow system dynamics. Detailed descriptions of historical and current Hg sources and transport pathways are also provided for the east and west ends of the facility. This conceptual model will allow the US Department of Energy to evaluate past and present remedial activities and provide a strong technical basis for prioritizing and optimizing remedial responses in a cost-effective and efficient manner. Twenty years ago, Hg flux from Outfall (OF) 200 at the headwaters of East Fork Poplar Creek represented approximately 20% of the overall flux leaving Y-12. By the time of the 2011 conceptual model report, that percentage had increased to 70%–80%. Flux estimates compiled for this report using data from 2009 to 2018 show that OF200 accounted for roughly 60% of the flux leaving Y-12. The relative role of the four storm drain conduits to OF200 appeared similar to 2011, with OF163 being the greatest contributor to downstream flux. An increase in Hg concentrations and flux occurred throughout and downstream of west end storm drains in 2011 in response to a storm drain cleanout. Leading up to the 2011 storm drain cleanout, annual baseflow flux was higher at OF200A6 than at Station 17. From 2011 to 2018, the trend reversed with concentrations and flux being higher at Station 17 than at OF200A6. A spike in concentrations occurred in July 2018, potentially due to a one-time influx of Hg that occurred during COLEX decontamination and decommissioning activities at Y-12. This report includes recommendations for improving our conceptual understanding of Hg sources, transport pathways, and flux at Y-12. Important recommendations include further evaluation of transport pathways, concurrent measurement of Hg concentration and flux under baseflow and stormflow conditions at multiple sites, and research to understand the connections among periphyton and methylmercury concentration, bioavailability, and bioaccumulation. Ongoing facility decontamination and decommissioning efforts, site characterization, remedial actions, and research are facilitating the collection of new Hg data at Y-12. The conceptual model update provided herein takes advantage of the extensive recent Hg sampling to provide an updated assessment of Hg mass balance and support a refined understanding of Hg behavior at or near Y-12. This assessment should assist in future environmental management decisions and in mitigating the impacts of Hg on the surrounding environment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Evaluation of Flow Routing on the Unstructured Voronoi Meshes in Earth System Modeling

Flow routing is a fundamental process of Earth System Models' (ESMs) river component. Traditional flow routing models rely on Cartesian rectangular meshes, which exhibit limitations, particularly when coupled with unstructured mesh-based ocean components. They also lack the support for regionally refined models. While previous studies have highlighted the potential benefits of unstructured meshes for flow routing, their widespread application and comprehensive evaluation within ESMs remain limited. This study extends the river component of the Energy Exascale Earth System Model to unstructured Voronoi meshes. We evaluated the model's performance in simulating river discharge and water depth across three watersheds spanning the Arctic, temperate, and tropical regions. The results show that while providing several benefits, unstructured mesh-based flow routing can achieve comparable performance to structured mesh-based routing, and their difference is often less than 10%. Although the unstructured mesh-based method could address several existing limitations, this research also shows that additional improvements in the numerical method are needed to fully exploit the advantages of unstructured mesh for hydrologic and ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The impact of hurricane disturbances on a tropical forest: implementing a palm plant functional type and hurricane disturbance module in ED2-HuDi V1.0

Hurricanes commonly disturb and damage tropical forests. Hurricane frequency and intensity are predicted to change under the changing climate. The short-term impacts of hurricane disturbances to tropical forests have been widely studied, but the long-term impacts are rarely investigated. Modeling is critical to investigate the potential response of forests to future disturbances, particularly if the nature of the disturbances is changing with climate. Unfortunately, existing models of forest dynamics are not presently able to account for hurricane disturbances. Therefore, we implement the Hurricane Disturbance in the Ecosystem Demography model (ED2) (ED2-HuDi). The hurricane disturbance includes hurricane-induced immediate mortality and subsequent recovery modules. The parameterizations are based on observations at the Bisley Experimental Watersheds (BEW) in the Luquillo Experimental Forest in Puerto Rico. We add one new plant functional type (PFT) to the model – Palm, as palms cannot be categorized into one of the current existing PFTs and are known to be an abundant component of tropical forests worldwide. The model is calibrated with observations at BEW using the generalized likelihood uncertainty estimation (GLUE) approach. The optimal simulation obtained from GLUE has a mean relative error of –21 %, –12 %, and –15 % for stem density, basal area, and aboveground biomass, respectively. The optimal simulation also agrees well with the observation in terms of PFT composition (+1 %, –8 %, –2 %, and +9 % differences in the percentages of “Early”, “Mid”, “Late”, and “Palm” PFTs, respectively) and size structure of the forest (+0.8 % differences in the percentage of large stems). Lastly, using the optimal parameter set, we study the impact of forest initial condition on the recovery of the forest from a single hurricane disturbance. The results indicate that, compared to a no-hurricane scenario, a single hurricane disturbance has little impact on forest structure (+1 % change in the percentage of large stems) and composition (<1 % change in the percentage of each of the four PFTs) but leads to 5 % higher aboveground biomass after 80 years of succession. The assumption of a less severe hurricane disturbance leads to a 4 % increase in aboveground biomass.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Inductive predictions of hydrologic events using a Long Short-Term Memory network and the Soil and Water Assessment Tool

We present machine learning methods to predict hydrologic features such as streamflow and soil moisture from spatially and temporally varying hydrological and meteorological data. Here, we used a temporal reduction technique to reduce computation and memory requirements and trained a Long Short-Term Memory (LSTM) network to predict soil moisture and streamflow over multiple watersheds. We show LSTM networks can be trained in a fraction of the time required by complex process-based and attention-based models such as Soil and Water Assessment Tool (SWAT) and GeoMAN without sacrificing accuracy. We also demonstrate that outside data - sourced from a watershed other than the target - can be used to train LSTM to comparable or even superior prediction accuracy. The success of LSTM in such spatially-inductive settings shows hydrologic features can be predicted with minimal prior knowledge of the watershed in question. Finally, we make all methodologies of this work publicly available as an end-to-end software pipeline that facilitates rapid prototyping of hydrologic learners.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Deep groundwater↗
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↗