Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling sediment diagenesis processes on riverbed to better quantify aquatic carbon fluxes and stocks in a small watershed of the Mid-Atlantic region

Background: Despite the widely recognized importance of aquatic processes for bridging gaps in the global carbon cycle, there is still a lack of understanding of the role of riverbed processes for carbon flows and stocks in aquatic environments. Here, we added a sediment diagenesis and sediment carbon (C) resuspension module into the SWAT-C model and tested it for simulating both particulate organic C (POC) and dissolved organic C (DOC) fluxes using 4 years of monthly observations (2014–2017) in the Tuckahoe watershed (TW) in the U.S. Mid-Atlantic region. Results: Sensitivity analyses show that parameters that regulate POC deposition in river networks are more sensitive than those that determine C resuspension from sediments. Further analyses indicate that allochthonous contributions to POC and DOC are about 36.6 and 46 kgC ha –1 year –1 , respectively, while autochthonous contributions are less than 0.72 kgC ha –1 year –1 for both POC and DOC (less than 2% of allochthonous sources). The net deposition of POC on the riverbed (i.e., 11.4 kgC ha –1 year –1 ) retained ca. 31% of terrestrial inputs of POC. In addition, average annual buried C was 0.34 kgC ha –1 year –1 , accounting for only 1% of terrestrial POC inputs or 3% of net POC deposition. The results indicate that about 79% of deposited organic C was converted to inorganic C (CH 4 and CO 2 ) in the sediment and eventually released into the overlying water column. Conclusion: This study serves as an exploratory study on estimation of C fluxes from terrestrial to aquatic environments at the watershed scale. We demonstrated capabilities of the SWAT-C model to simulate C cycling from uplands to riverine ecosystems and estimated C sinks and sources in aquatic environments. Overall, the results highlight the importance of including carbon cycle dynamics within the riverbed in order to accurately estimate aquatic carbon fluxes and stocks. The new capabilities of SWAT-C are expected to serve as a useful tool to account for those processes in watershed C balance assessment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modelling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California

ABSTRACT Coastal wetlands, some of the most productive ecosystems on Earth, provide critical ecosystem services, including support of biodiversity, carbon sequestration and flood protection. In recent decades, these ecosystems have experienced extensive coastal wetland loss. Coastal wetland restoration provides a beacon of hope, offering a chance to reclaim these important habitats. However, even with billions of dollars invested worldwide in restoring coastal wetlands, we still lack comprehensive knowledge about the effectiveness of these restoration efforts in recovering wetland ecosystem functions and how future climate change may affect these efforts. The ability to evaluate how these ecosystems will function in the future is vital for examining current investments and developing future protection and management plans. We selected Elkhorn Slough, a tidal estuary, in California, to investigate the impact of wetland restoration and sea level rise (SLR) on coastal hydrology using the process‐based coastal hydrologic model, Advanced Terrestrial Simulator (ATS), informed by site‐specific data. We designed a novel modelling workflow for incorporating wetland restoration features into land cover and soil properties for the model parameterization. The validation results demonstrate a strong agreement between modelled and observed data. We studied the characteristics of coastal watershed hydrology, then focused on the surface water dynamics at two wetland sites within Elkhorn Slough, a reference site and a restored site. Our simulation results indicate that the restored site successfully maintains surface elevation, resulting in reduced surface inundation. We also examined the impact of wetland restoration under expected SLR over the next few decades. The low‐lying Yampah Marsh, the reference site, is likely to be inundated due to future SLR when highest tides arrive, while a higher percentage of Hester Marsh, the restored site, would retain marsh vegetation in coming decades, regardless of tidal conditions. Our study provides important information for examining the outcome of restoration practices that include surface elevation in tidal wetlands under climate changes.

advanced terrestrial simulator↗
DOE OSTI2023

Cycles-L: A Coupled, 3-D, Land Surface, Hydrologic, and Agroecosystem Landscape Model

Managing landscapes to increase agricultural productivity and environmental stewardship can be informed by spatially-distributed models that operate at spatial and temporal scales that are intervention-relevant. This paper presents Cycles-L, a landscape-scale agroecosystem and hydrologic modeling system, using as a test case a watershed in Pennsylvania. Cycles-L emerges from melding the landscape and hydrology structure of Flux-PIHM, a 3-D land surface hydrologic model, with the agroecosystem processes in the Cycles model. Consequently, Cycles-L can simulate processes affected by topography, soil heterogeneity, and management practices, owing to its physically-based hydrology that can simulate horizontal and vertical transport of solutes with water. The model was tested at a 730-ha experimental watershed within the Mahantango Creek watershed. Cycles-L simulated well stream water and mineral nitrogen discharge (Nash-Sutcliffe coefficient 0.55 and 0.60, respectively) and grain yield (root mean square error 1.2 Mg ha -1 ). Cycles-L outputs are as good or better than those obtained with the uncoupled Flux-PIHM (water discharge) and Cycles (grain yield) models. Modeled spatial patterns of nitrogen fluxes like denitrification illustrate the combined control of crop management and topography. For example, denitrification is almost twice as high when simulated with Cycles-L than when simulated with Cycles 1-D. Due to its spatial and temporal resolution, Cycles-L fills a gap in the availability of models that operate at a scale relevant to evaluate interventions in the landscape. Cycles-L can become a central component in tools for climate change scenario analysis, precision agriculture, precision conservation, and artificial intelligence-based decision support systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

WHONDRS-GUI: a web application for global survey of surface water metabolites

Background The Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) is a consortium that aims to understand complex hydrologic, biogeochemical, and microbial connections within river corridors experiencing perturbations such as dam operations, floods, and droughts. For one ongoing WHONDRS sampling campaign, surface water metabolite and microbiome samples are collected through a global survey to generate knowledge across diverse river corridors. Metabolomics analysis and a suite of geochemical analyses have been performed for collected samples through the Environmental Molecular Sciences Laboratory (EMSL). The obtained knowledge and data package inform mechanistic and data-driven models to enhance predictions of outcomes of hydrologic perturbations and watershed function, one of the most critical components in model-data integration. To support efforts of the multi-domain integration and make the ever-growing data package more accessible for researchers across the world, a Shiny/R Graphical User Interface (GUI) called WHONDRS-GUI was created. Results The web application can be run on any modern web browser without any programming or operational system requirements, thus providing an open, well-structured, discoverable dataset for WHONDRS. Together with a context-aware dynamic user interface, the WHONDRS-GUI has functionality for searching, compiling, integrating, visualizing and exporting different data types that can easily be used by the community. The web application and data package are available at https://data.ess-dive.lbl.gov/view/doi:10.15485/1484811 , which enables users to simultaneously obtain access to the data and code and to subsequently run the web app locally. The WHONDRS-GUI is also available for online use at Shiny Server ( https://xmlin.shinyapps.io/whondrs/ ).

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1981

Remote sensing applications to resource problems in South Dakota

The procedures used as well as the results obtained and conclusions derived are described for the following applications of remote sensing in South Dakota: (1) sage grouse management; (2) censusing Canada geese; (3) monitoring grasshopper infestation in rangeland; (4) detecting Dutch elm disease in an urban environment; (5) determining water usage from the Belle Fourche River; (6) resource management of the Lower James River; and (7) the National Model Implantation Program: Lake Herman watershed.

Myers, V. I.↗
NASA NTRS1982

The application of remote sensing to resource management and environmental quality programs in Kansas

Capabilities to process data for state agencies in Kansas were upgraded through the vehicle of a low cost data processing system. Short term projects in which agencies identified areas of immediate needs, and longer terms projects, continued from previous years, are described including studies of Arkansas River irrigation; evaluation of rangeland in the Cimarron National Grassland; a model of the Walnut Creek watershed groundwater; selection of a pronghorn antelope release site; the establishment of a geographical data base for tax reassessment in Finney County; a land use/land cover inventory and hazards assessment; and applied R & D in agricultural remote sensing. The topics discussed at a short course in remote sensing and publications are listed.

Barr, B. G.↗
NASA NTRS2008

Assimilation of GRACE Terrestrial Water Storage Data into a Land Surface Model

The NASA Gravity Recovery and Climate Experiment (GRACE) system of satellites provides observations of large-scale, monthly terrestrial water storage (TWS) changes. In. this presentation we describe a land data assimilation system that ingests GRACE observations and show that the assimilation improves estimates of water storage and fluxes, as evaluated against independent measurements. The ensemble-based land data assimilation system uses a Kalman smoother approach along with the NASA Catchment Land Surface Model (CLSM). We assimilated GRACE-derived TWS anomalies for each of the four major sub-basins of the Mississippi into the Catchment Land Surface Model (CLSM). Compared with the open-loop (no assimilation) CLSM simulation, assimilation estimates of groundwater variability exhibited enhanced skill with respect to measured groundwater. Assimilation also significantly increased the correlation between simulated TWS and gauged river flow for all four sub-basins and for the Mississippi River basin itself. In addition, model performance was evaluated for watersheds smaller than the scale of GRACE observations, in the majority of cases, GRACE assimilation led to increased correlation between TWS estimates and gauged river flow, indicating that data assimilation has considerable potential to downscale GRACE data for hydrological applications. We will also describe how the output from the GRACE land data assimilation system is now being prepared for use in the North American Drought Monitor.

Reichle, Rolf H.↗
DOE OSTI2026

Water Observations of Flow/No-Flow for the East-Taylor Watershed, Colorado (June-July 2025 and 2026)

This dataset provides multi-year, ground-truth visual observations of surface water flow/no-flow conditions within the East-Taylor Watershed, Colorado, collected during June and July of 2025 and 2026. In June and July 2025, on-the-ground visual observations of flow/no-flow were collected as part of the Watershed Function Scientific Focus Area (SFA) and Rocky Mountain Biological Laboratory (RMBL) Colorado Headwaters Ecological Spectroscopy Study (CHESS) campaign (further details are provided within the CHESS Project Description). We obtained 377 water observations of flow/no-flow within the East-Taylor Watershed, Colorado. These ground-truth observations were collected to validate classification maps from remote sensing data and model results within the East-Taylor Watershed. In 2025, flow/no-flow measurements were collected using a field-based app for the CHESS Campaign (Zerion iForm). Within the field app, a water observation form was created to collect coordinates and metadata about the observation. Information collected for the water observation points included information about visually-assessed streamflow presence/absence (standard question obtained from Colorado State University’s StreamTracker project), flow estimate, stream or ponded area width, canopy cover, manganese films, iron seeps, and beaver activity. For 2025 water observations, this dataset contains: (1) a data file with the water observations and coordinates (2025_Water_Observations.csv); (2) a Keyhole Markup Language Zipped (KMZ) with the water observation locations and metadata (2025_Water_Observations_Locations.kmz); (3) photos (.jpg and .jpeg) of the water observation points, organized by location, contained within 2025_Water_Observations_FieldPhotographs.zip file; and (4) water observation protocols and figures (2025_Water_Observation_Protocols.pdf). In June and July 2026, on-the-ground visual observations of flow/no-flow were collected as part of the Watershed Function SFA project. We obtained 365 water observations of flow/no-flow within the East-Taylor Watershed, Colorado. The 2026 observations focused on collecting repeat measurements at the 2025 flow/no-flow observation locations conducted as part of the CHESS campaign. These ground-truth observations were collected to understand differences in flow/no-flow in 2026, given the unprecedented 2026 drought in Colorado. In 2026, flow/no-flow measurements were collected using ArcGIS (Geographic Information System) Survey123. Within the field app, a water observation form was created to collect coordinates and metadata about the observation. Information collected for the water observation points included repeat information from the 2025 water observation effort, including visually-assessed streamflow presence/absence (standard question obtained from Colorado State University’s StreamTracker project), flow estimate, stream or ponded area width, canopy cover, manganese films, iron seeps, beaver activity, and a new metadata component of estimated stream depth (for select locations). For 2026 water observations, this dataset contains: (1) a data file with the water observations and coordinates (2026_Water_Observations.csv); (2) a Keyhole Markup Language Zipped (KMZ) with the water observation locations and metadata (2026_Water_Observations_Locations.kmz); (3) photos (.jpg) of the water observation points, organized by location, contained within 2026_Water_Observations_FieldPhotographs.zip file; and (4) water observation protocols and figures (2026_Water_Observation_Protocols.pdf). For 2025 and 2026 water observations, this dataset contains: (1) a location metadata file (locations.csv); (6) a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and (7) a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. 2026-09-02: This dataset was updated to include 2026 water observation measurements. The 2025 observation files were also updated to ensure a consistent file naming convention across water observation years.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗
DOE OSTI2026

Model Data Archive Associated with Manuscript "Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon"

This data package supports the publication “Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon” by Li et al. (2026). The package contains processed model inputs, configuration files, restart files, simulation outputs, scripts, and visualization products used to evaluate post-fire dissolved organic carbon (DOC) dynamics in the Naches River Watershed, Washington, USA, following the 2021 Schneider Springs Fire. The modeling workflow couples ELM-BGC, the biogeochemistry-enabled Energy Exascale Earth System Model Land Model; ATS, the Advanced Terrestrial Simulator for integrated surface-subsurface hydrology; and PFLOTRAN, a reactive transport model for multicomponent aqueous geochemistry. Together, these models simulate how wildfire-induced changes in vegetation, litter, coarse woody debris, and soil organic matter influence DOC production, transport, and reaction from burned hillslopes to stream networks. The archive includes preprocessed meteorological, geospatial, hydrologic, and biogeochemical forcing data; ELM-BGC-derived DOC source terms; ATS mesh files; PFLOTRAN reactive-transport inputs; model configuration files; spin-up and transient restart files; watershed-scale diagnostic outputs; stream concentration time series; and figures or visualization files used to inspect and reproduce key results. File types include Hierarchical Data Format 5 (HDF5) files for gridded forcing and model-coupling data, model input and configuration files for ELM-BGC, ATS, and PFLOTRAN, restart and simulation-output files generated by the modeling workflow, tabular or time-series diagnostic outputs, scripts for post-processing and figure generation, and image or visualization products associated with the manuscript. Use of the package depends on the intended task. Re-running the simulations requires the relevant modeling software, including ELM-BGC, ATS, and PFLOTRAN as ATS's geochemical engine. Inspecting outputs and reproducing figures requires Python with scientific plotting libraries such as Matplotlib, and three-dimensional model outputs may be viewed with ParaView. Geographic information system files or maps may be inspected with ArcGIS Pro or comparable GIS software. The data package is intended to enable traceability, reuse, and partial reproduction of the coupled land-to-watershed hydro-biogeochemical modeling workflow used to test how wildfire disturbance affects terrestrial carbon pools and downstream DOC dynamics.

ATS↗
DOE OSTI2022

rabpro: global watershed boundaries, river elevation profiles, and catchment statistics

River and Basin Profiler (rabpro) is a Python package to delineate watersheds, extract river flowlines and elevation profiles, and compute watershed statistics for any location on the Earth’s surface. As fundamental hydrologically-relevant units of surface area, watersheds are areas of land that drain via aboveground pathways to the same location, or outlet. Delineations of watershed boundaries are typically performed on digital elevation models (DEMs) that represent surface elevations as gridded rasters. Depending on the resolution of the DEM and the size of the watershed, delineation may be very computationally expensive. With this in mind, we designed rabpro to provide user-friendly workflows to manage the complexity and computational expense of watershed calculations given an arbitrary coordinate pair. In addition to basic watershed delineation, rabpro will extract the elevation profile for a watershed’s mainchannel flowline. This enables the computation of river slope, which is a critical parameter in many hydrologic and geomorphologic models. Finally, rabpro provides a user-friendly wrapper around Google Earth Engine’s (GEE) Python API to enable cloud-computing of zonal watershed statistics and/or time-varying forcing data from hundreds of available datasets. Altogether, rabpro provides the ability to automate or semi-automate complex watershed analysis workflows across broad spatial extents.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics Within Water‐Tagging Enabled Hydrologic Models

Determining the age distribution of water exiting a catchment is important for understanding groundwater storage and mixing. New water-tagging capabilities within models track precipitation events as they move through simulated storages, yet forward modeling of individual events may not systematically capture the full transit time distribution (TTD). Here, we present a “sequential precipitation input tagging” (SPIT) framework to tag all input precipitation at regular intervals during extended model simulations. Monthly tags over 7 years were applied at six National Ecological Observatory Network sites to calculate TTDs and derive mean virtual tracer age, $\overline{T_{V}}$, fractions of young water, F yw , and hydrologic tracer concentrations (water isotopes δ 18 O and δ 2 H) within a tagging enabled version of the Weather Research and Forecast hydrologic model (WRF-Hydro). Throughout seven simulation years, the fraction of simulated discharge derived from tagged events, F tag , increased each year, with the final year's F tag ranging from 66% to 100% and highlights the need to apply SPIT over many years to understand TTDs. When the F tag was >75%, simulated $\overline{T_{V}}$ ranged 179–923 days and F yw 0.6%–23.9%, with daily values exhibiting a power-law relationship with precipitation, discharge, and groundwater. Through implementation of SPIT, we find this hydrologic model configuration performs poorly in estimation of $\overline{T_{V}}$ and F yw (root mean squared error of 469 days and 14.4% respectively), suggesting it misrepresents subsurface mixing. Thus, the SPIT framework provides a reproducible approach to calculate watershed transit times within tagging enabled models and thereby assess and improve representation of hydrologic processes.

fraction of young water↗
DOE OSTI2022

Improving tide-estuary representation in MPAS-Ocean (Final Report)

The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.

58 GEOSCIENCES↗
DOE OSTI2024

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

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

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2013

Application and Evaluation of a Snowmelt Runoff Model in the Tamor River Basin, Eastern Himalaya Using a Markov Chain Monte Carlo (MCMC) Data Assimilation Approach

Previous studies have drawn attention to substantial hydrological changes taking place in mountainous watersheds where hydrology is dominated by cryospheric processes. Modelling is an important tool for understanding these changes but is particularly challenging in mountainous terrain owing to scarcity of ground observations and uncertainty of model parameters across space and time. This study utilizes a Markov Chain Monte Carlo data assimilation approach to examine and evaluate the performance of a conceptual, degree-day snowmelt runoff model applied in the Tamor River basin in the eastern Nepalese Himalaya. The snowmelt runoff model is calibrated using daily streamflow from 2002 to 2006 with fairly high accuracy (average Nash-Sutcliffe metric approx. 0.84, annual volume bias <3%). The Markov Chain Monte Carlo approach constrains the parameters to which the model is most sensitive (e.g. lapse rate and recession coefficient) and maximizes model fit and performance. Model simulated streamflow using an interpolated precipitation data set decreases the fractional contribution from rainfall compared with simulations using observed station precipitation. The average snowmelt contribution to total runoff in the Tamor River basin for the 2002-2006 period is estimated to be 29.7+/-2.9% (which includes 4.2+/-0.9% from snowfall that promptly melts), whereas 70.3+/-2.6% is attributed to contributions from rainfall. On average, the elevation zone in the 4000-5500m range contributes the most to basin runoff, averaging 56.9+/-3.6% of all snowmelt input and 28.9+/-1.1% of all rainfall input to runoff. Model simulated streamflow using an interpolated precipitation data set decreases the fractional contribution from rainfall versus snowmelt compared with simulations using observed station precipitation. Model experiments indicate that the hydrograph itself does not constrain estimates of snowmelt versus rainfall contributions to total outflow but that this derives from the degree-day melting model. Lastly, we demonstrate that the data assimilation approach is useful for quantifying and reducing uncertainty related to model parameters and thus provides uncertainty bounds on snowmelt and rainfall contributions in such mountainous watersheds.

runoff↗
DOE OSTI2022

Impacts of climate change on future hurricane induced rainfall and flooding in a coastal watershed: A case study on Hurricane Harvey

The warming climate is likely to increase hurricane-associated extreme rainfall and lead to sea-level rise (SLR). Thus, how the floods induced by intense hurricanes respond to these potential changes is of great concern. This study investigates the future warmer climate impacts on hurricane-induced extreme rainfall, and—more importantly—the subsequent compound flooding at the watershed scale (from an event-based analysis perspective). To this goal, a modeling framework is designed based on the Distributed Hydrology Soil Vegetation Model (DHSVM), the Two-Dimensional fvand the Regional Community Earth System Model (R-CESM). The framework was applied to Hurricane Harvey (2017) at the Clear Creek watershed (a coastal watershed in the southern Houston) as a case study. The results show that the projected maximum rainfall totals over the watershed would be exacerbated by 17.7 % and 49.7 % in the 2050s and 2090s (respectively) under Representative Concentration Pathway 8.5 (RCP 8.5). This means a 16.1 % increase in Harvey rainfall over the watershed per degree Celsius increase in Mean Surface Temperature over the Gulf of Mexico region (18°~31° N, 77°~98° W). Meanwhile, the increases in maximum inundation extent would be 11.0 % (2050s) and 19.5 % (2090s). Furthermore, considerable increases in maximum inundation depth and duration in regions along the middle and downstream of Clear Creek (and also those around Clear Lake) are expected. The projected SLR will have little effect on the maximum inundation depth and extent if storm surge changes are not taken into account; meanwhile, it will influence the inundation duration at downstream locations. In conclusion, this modeling framework can be also applied at other coastal watersheds to evaluate the projected climate change impacts on the compound flooding induced by extreme climate events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

ParFlow-CLM model simulation files for Maina et al., Journal of Hydrology, 2022

This dataset contains the files to run a ParFlow-CLM integrated hydrologic model simulation for Maina et al., HESS, 2022. It also contains the associated daily pressure output of those simulations. Simulations are for a hillslope-similarity approach based on seasonal groundwater changes in the East River, Colorado. Comparisons are made to different clustering, or functional zonation approaches. We assess the ability of these clustering approaches to identify and categorize hillslopes with similar static characteristics, hydroclimate, land surface processes, and subsurface dynamics in a mountainous watershed. ParFlow binary files (.pfb) are model specific, and can be modified/viewed with conversion to .sa, .silo, or .vtk file formats. ParFlow input scripts are in the .tcl file format. Associated post-processing scripts are included as .f90 files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Location Identifiers, Metadata, and Map for Field Measurements at the East-Taylor Watershed Community Observatory, Colorado, USA (Version 3.3)

This dataset contains identifiers, metadata, and a map of the locations where field measurements have been conducted at the East-Taylor Watershed Community Observatory located in the Upper Colorado River Basin, United States. This is version 3.3 of the dataset and replaces the prior version 3.2 (see below for details on changes between the versions). Dataset description: The East River-Taylor Watershed is the primary field site of the Watershed Function Scientific Focus Area (WFSFA) and the Rocky Mountain Biological Laboratory. Researchers from several institutions generate highly diverse hydrological, biogeochemical, climate, vegetation, geological, remote sensing, and model data at the East-Taylor Watershed in collaboration with the WFSFA. Thus, the purpose of this dataset is to maintain an inventory of the field locations and instrumentation to provide information on the field activities in the East-Taylor Watershed and coordinate data collected across different locations, researchers, and institutions. The dataset contains (1) a README file with information on the various files, (2) three csv files describing the metadata collected for each surface point location, plot and region registered with the WFSFA, (3) csv files with metadata and contact information for each surface point location registered with the WFSFA, (4) a csv file with with metadata and contact information for plots, (5) a csv file with metadata for geographic regions and sub-regions within the watershed, (6) a compiled xlsx file with all the data and metadata which can be opened in Microsoft Excel, (7) a kml map of the locations plotted in the watershed which can be opened in Google Earth, (8) a jpg image of the kml map which can be viewed in any photo viewer, and (9) a zipped file with the registration templates used by the SFA team to collect location metadata. The zipped template file contains two csv files with the blank templates (point and plot), two csv files with instructions for filling out the location templates, and one compiled xlsx file with the instructions and blank templates together. Additionally, the templates in the xlsx include drop down validation for any controlled metadata fields. Persistent location identifiers (Location_ID) are determined by the WFSFA data management team and are used to track data and samples across locations. Dataset uses: This location metadata is used to update the Watershed SFA’s publicly accessible Field Information Portal (an interactive field sampling metadata exploration tool; https://wfsfa-data.lbl.gov/watershed/), the kml map file included in this dataset, and other data management tools internal to the Watershed SFA team. Version Information: The latest version of this dataset publication is version 3.3. This version contains 167 new point locations, 1 new plot, and 2 new geographic regions. Overall, there are a total of 1439 point locations, 75 plots, and 54 geographic regions. Additionally, the kml map of locations and image now includes two boundaries (Upper Ohio Creek (UO) and Carbon Creek (CA)) outside of the East River watershed (USGS HUC-10) and accompanying stream network that represents areas of focus. Refer to methods for further details on the version history. This dataset will be updated on a periodic basis with new measurement location information. Researchers interested in having their East-Taylor Watershed measurement locations added to this list should reach out to the WFSFA data management team at wfsfa-data@googlegroups.com. Acknowledgments: Please cite this dataset if using any of the location metadata in other publications or derived products. If using the location metadata for the 2018 NEON hyperspectral campaign, additionally cite Chadwick et al. (2020). doi:10.15485/1618130. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. Part of this work was performed at SLAC Accelerator Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-76SF00515.

2018 NEON and 2025 CHESS Campaigns↗