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 181 records · Page 10

Simulating water dynamics related to pedogenesis across space and time: Implications for four-dimensional digital soil mapping

Digital soil mapping (DSM) relies on machine-learning and geostatistics to represent soil property observations across space. DSM techniques are powerful but often empirical, being limited to the quality and density of point samples. Water dynamics are closely related to soil variability, and the physics that govern water movement are well known. Hydrological properties can hence be simulated by physical models through space and time, unveiling key characteristics about soils. We propose the use of hydrologic models to map soils across the surface (2D), depth (1D), and time (1D)–which provides a 4D approach to digital soil mapping (4DSM). The Distributed Hydrology Soil Vegetation Model (DHSVM) was applied to a watershed currently under pasture. Moisture sensors and wells were installed at different depths in the watershed on summit, sideslope and toeslope positions to validate the model. DHSVM simulations of soil moisture distribution and depth to saturation were performed during the hydrological year (October 2008-September 2009). Clusters of similar pixels based on soil moisture values were determined using Dynamic Time Warping (DTW) to align temporal data and K-means. Clustering was performed both seasonally and for the entire year. Temporal patterns simulated by DHSVM matched measurements given by moisture sensors and wells. Seasonal clusters differed from the annual cluster. Distinct clusters were observed for each season and with depth, showing that spatiotemporal soil variability is lost when statically assessing soils. Spatiotemporal clusters corroborated field observations of fragipan occurrence not explicitly spatially mapped by Soil Survey Geographic Database (SSURGO). If a connection can be made between water and soils, static and dynamic soil variability can be predicted using physically based hydrologic models. Hydrologic models can benefit soil mapping by enabling reliable 4D simulation of water dynamics, which are fundamental to soil variability and soil classification and directly relate to biological, physical and chemical soil processes not captured by typical soil sampling protocols.

54 ENVIRONMENTAL SCIENCES↗

Latest Pleistocene glacial chronology and paleoclimate reconstruction for the East River watershed, Colorado, USA

Reconstructing Pleistocene glaciation timing and extent is vital for understanding paleoclimate. Whereas late Pleistocene glaciation has been studied extensively in western North American mountain ranges, the glacial history of the western Elk Range in Colorado remains understudied, particularly in the East River watershed, a site of intense scientific focus. Here we use cosmogenic nuclide exposure and depth–profile dating methods to determine the timing of glaciation in the East River watershed. We use glacier modeling to reconstruct paleoglacier extents and quantify past climate conditions. Our findings indicate that the East River glacier retreated from its maximum position approximately 17–18 ka, moving to recessional positions between 13 and 15 ka, before experiencing more substantial retreat to high-elevation cirques around 13 ka. Glacier modeling suggests that the maximum ice extents at 17–18 ka could have been sustained by temperature depressions of approximately −6.5°C compared to modern conditions, assuming consistent precipitation. Additionally, the ice position at 13–15 ka could have been supported by temperature depressions of around −4.0°C. These results offer insights into the deglaciation timeline in the East River watershed and broader western Elk Range as well as paleoclimate conditions during the late Pleistocene, which may aid future research on critical zone evolution in the East River watershed.

54 ENVIRONMENTAL SCIENCES↗

Hydrologic land use classification of the Patuxent River watershed using remotely sensed data

The Patuxent River Watershed is located in central Maryland between Baltimore and Washington, D.C. and is approximately 2330 sq km in area and 175 km long. This region is now at a critical point because of major concerns such as water management and quality, flooding and land use within the watershed. Data from the NASA-directed LANDSAT and Earth Resources Aircraft Programs were used to provide a new dimension in information collection and processing for the management of watersheds. Digital data from LANDSAT-1 were analyzed along with selected IR photography from U-2 flight number 74-060B taken 28 April 1974, which was digitized in three channels. Processing of the data was accomplished using a multispectral analysis system. Land use themes consisting of surface water, wetlands, forest, residential, cropland/pasture, urban, and extractive were developed and delineated through the watershed. Area measurements of watershed themes were obtained and will serve as a calibration input to a deterministic hydrologic model on a sub-watershed. Using the derived residential and urban theme areas from LANDSAT an estimated basin imperviousness was also calculated. Thematic maps were produced at 1:62,500 scale. Floodprone areas were also classified and delineated at a scale of 1:24,000. Comparison with standard floodprone area maps at the same scale have indicated a few areas of discrepancy. Such information can be used for updating or checking floodprone area boundaries as well as monitoring changes in floodplain areas.

Dallam, W. C.↗

Denudation, solute export, landscape evolution modeling, and geographic information system data for the East River watershed, Colorado, USA (2020-2024)

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of lithologic controls on denudation, solute export, carbon-scaling relationships, and transient landscape evolution in the East River watershed near Crested Butte, Colorado, USA. The package includes GIS layers used to produce the Figure 2 map, including drainage, hillshade, lithology, sample locations, and basin polygons, together with comma-separated value (CSV) tables and matching CSV data dictionaries. One group of tables reports sample-level and catchment-level information for river-sediment samples analyzed for in situ-produced cosmogenic beryllium-10 (10Be), including sample names, outlet elevations, geographic coordinates, upstream drainage area, rock-type classes, production-rate scaling scheme, analyzed nuclide, catchment-averaged denudation rates, and associated lower and upper analytical uncertainties. Sample and catchment attributes provide the basis for comparing denudation rates across intrusive, shale, sedimentary, and mixed-lithology settings. A second group of tables reports supporting information for landscape-evolution modeling and the mapped geologic framework of the study area. Included files list parameter values and definitions for the two-phase landscape-evolution simulations, summarize full-domain model erosion fluxes and topographic metrics for different simulation configurations, provide a fixed-area carbon-model scaling table, and summarize mapped geologic units within the East River study domain, including geologic code, formation name, lithologic description, mapped area, and lithologic class grouping. Model outputs and geologic summaries support interpretation of transient landscape behavior and its relation to the mapped distribution of shale, intrusive, sedimentary, and surficial units. A third group of tables reports hydrologic and hydrochemical information used to quantify dissolved export from the watershed. Included files provide site-level values for drainage area, mean annual solute export, standard error of annual export, area-normalized solute yield, and equivalent weathering rate for five East River monitoring sites, along with metadata describing the number, sampling cadence, and date range of discharge records and partial and full total dissolved solids observations used in the solute-yield analyses. The package also contains a supplementary daily ion-load time series with daily mean discharge, discharge observation counts, dissolved concentrations, and daily loads for calcium, magnesium, sodium, potassium, chloride, sulfate, nitrate, fluoride, dissolved silica, charge-balance bicarbonate, and total dissolved solids. The package contains GIS files, comma-separated value files (.csv), CSV data dictionaries, a file-level metadata table, a package-tree text file, and a readme text file.

10Be↗

A path to intelligent watersheds: coordinating the data to decision pipeline

Operations of multi-reservoir systems are challenged in-part by the interplay of complex physical processes functioning within the watershed. The employment of intelligent systems can be of aid by linking environmental sensing, information technology, data analytics, simulation and decision support to achieve a data-to-decision flow of information. A further challenge is that watershed resources are managed for multiple purposes requiring some level of coordination among numerous resource managers, asset operators and users. System intelligence in this context relies on shared community platforms (data portals, community models), and coordinated communication between decision makers. Opportunities to enrich watershed intelligence has been the subject of a roadmapping exercise for the Department of Energy’s Water Power Technologies Office which has relied on broad stakeholder engagement. Initial phases of engagement involved personal interviews and a series of virtual group meetings, which focused on identifying opportunities to improve the intelligence of the physical infrastructure within our watersheds—examples of feedback include improved sensing of snowpack and runoff, data standards for facilitated data sharing, and better forecasting tools. The latter phase of engagement involved the conduct of a case study in the Upper Colorado River basin where key stakeholders were interviewed to map how their decisions are informed by intelligence from other basin stakeholders. Our presentation will highlight the interdisciplinary flow of information in complex watershed systems and identify physical and institutional opportunities toward the strategic operation of water infrastructure.

Colorado River↗

Association of small vessel disease with tau pathology

Emerging evidence suggests that small vessel disease (SVD) is a risk factor for clinical dementia and may contribute to AD neuropathological changes. Watershed brain regions are located at the most distal areas between arterial territories, making them vulnerable to SVD-related changes. We examined the association of pathologic markers of SVD, specifically arteriolosclerosis in watershed brain regions, with AD pathologic changes. Participants (N= 982; mean age-at-death = 90; 69% women) were enrolled as part of one of two cohort studies of aging and dementia. At autopsy, neuropathological evaluation included semi-quantitative grading of arteriolosclerosis pathology from 2 cortical watershed regions: the anterior watershed (AWS) and posterior watershed (PWS), densities for cortical β-amyloid and tau-tangle pathology, and other common age-related pathologies. Linear regression models examined the association of watershed arteriolosclerosis pathology with β-amyloid and tau-tangle burden. In follow-up analyses, available ex-vivo MRI and proteomics data in a subset of decedents were leveraged to examine the association of whole brain measure of WMH, as a presumed MRI marker of SVD, with β-amyloid and tau-tangle burden, as well as to examine the association of watershed arteriolosclerosis with proteomic tau. Watershed arteriolosclerosis was common, with 45% of older persons having moderate-to-severe arteriolosclerosis pathology in the AWS region, and 35% in the PWS. In fully adjusted models that controlled for demographics and common age-related pathologies, an increase in severity of PWS arteriolosclerosis was associated with a higher burden of tau-tangle burden, specifically neocortical tau burden, but not with β-amyloid. AWS arteriolosclerosis was not associated with β-amyloid or tau pathology. Ex-vivo WMH was associated with greater tau-tangle pathology burden but not β-amyloid. Furthermore, PWS arteriolosclerosis was associated with higher abundance of tau phosphopeptides, that promote formation of tau aggregates. These data provide compelling evidence that SVD, specifically posterior watershed arteriolosclerosis pathology, is linked with tau pathological changes in the aging brain.

60 APPLIED LIFE SCIENCES↗

Model Scripts for "Old-Aged Groundwater Contributes to Mountain Hillslope Hydrologic Dynamics"

The partitioning of water inputs between deep and shallow groundwater flow paths is a fundamental processes, yet is challenging to observe. Numerical models provide a valuable tool to further develop insights on these groundwater mixing processes. This repository contains scripts to run the ParFlow-CLM and EcoSLIM integrated hydrologic models along the Pumphouse Hillslope in the East River Watershed and associated python scripts to process model outputs. Files includes input decks for the associated models, python scripts, and text and css files to support the model runs and processing. The 2-D hillslope model simulates from 2000 to 2021 using transient forcing conditions from the publicly available NLDAS-2 dataset. We develop an ensemble of models with variable hydrogeologic and soil parameters. Model outputs from the ensemble of runs are compared to to water level observations at the PLM1 and PLM6 wells published at (https://data.ess-dive.lbl.gov/view/doi:10.15485/1866836) and groundwater mean ages published at (https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1960042. The model results are used to evaluate hillslope water mass-balance transience and mixing dynamics between groundwater flow paths with young (<10 year) and old (>10 year) ages.

54 ENVIRONMENTAL SCIENCES↗

Hydrology Outweighs Temperature in Driving Production and Export of Dissolved Carbon in a Snowy Mountain Catchment

Abstract Terrestrial production and export of dissolved organic and inorganic carbon (DOC and DIC) to streams depends on water flow and biogeochemical processes in and beneath soils. Yet, understanding of these processes in a rapidly changing climate is limited. Using the watershed‐scale reactive‐transport model BioRT‐HBV and stream data from a snow‐dominated catchment in the Rockies, we show deeper groundwater flow averaged about 20% of annual discharge, rising to ∼35% in drier years. DOC and DIC production and export peaked during snowmelt and wet years, driven more by hydrology than temperature. DOC was primarily produced in shallow soils (1.94 ± 1.45 gC/m 2 /year), stored via sorption, and flushed out during snowmelt. Some DOC was recharged to and further consumed in the deeper subsurface via respiration (−0.27 ± 0.02 gC/m 2 /year), therefore reducing concentrations in deeper groundwater and stream DOC concentrations at low discharge. Consequently, DOC was primarily exported from the shallow zone (1.62 ± 0.96 gC/m 2 /year, compared to 0.12 ± 0.02 gC/m 2 /year from the deeper zone). DIC was produced in both zones but at higher rates in shallow soils (1.34 ± 1.00 gC/m 2 /year) than in the deep subsurface (0.36 ± 0.02 gC/m 2 /year). Deep respiration elevated DIC concentrations in the deep zone and stream DIC concentrations at low discharge. In other words, deep respiration is responsible for the commonly‐observed increasing DOC concentrations (flushing) and decreasing DIC concentrations (dilution) with increasing discharge. DIC export from the shallow zone was ~66% of annual export but can drop to ∼53% in drier years. Numerical experiments suggest lower carbon production and export in a warmer, drier future, and a higher proportion from deeper flow and respiration processes. These results underscore the often‐overlooked but growing importance of deeper processes in a warming climate.

54 ENVIRONMENTAL SCIENCES↗

ESS-DIVE Reporting Format for Amplicon Abundance Table

While standardized sequencing data is available in public repositories and efforts such as MIxS for common sample collection and processing metadata are well established, the lack of common bioinformatic processing metadata has hindered the ability to do large-scale metaanalyses and the potential for data re-use by non-experts such as ecosystem, watershed, or earth system modelers. To address this need for Department of Energy researchers, we have developed an amplicon reporting format which captures both sample preparation and bioinformatic processing metadata and stores processed amplicon data as a paired abundance table and sequencing file to maximize the potential for re-use of these data. To aid in the adoption of accessible and reproducible analysis workflows, this reporting format was developed in concert with amplicon functionality within the Department of Energy’s Systems Biology Knowledgebase (KBase) to ensure common data and metadata requirements and facilitate seamless transfer between these platforms.This dataset contains support documentation for the amplicon reporting format (README.md and instructions.md), templates for both bioinformatic and sequencing metadata (amplicon_bioinformatic_metadata_template_2021_10_03.csv and amplicon_sequencing_metadata_template_2021_10_03.csv), a crosswalk indicating how this reporting format relates to the current MIxS format (ESSDIVE-MIxS_crosswalk.csv), a list of available instrument terms (amplicon_seq_instrument_terms_2021_10_03.csv), a map between QIIME2 parameter settings and metadata fields (amplicon_qiime2_plugin_metadata_map.csv), a data dictionary (amplicon_CSV_dd.csv), and file-level metadata (amplicon_FLMD.csv).

54 ENVIRONMENTAL SCIENCES↗

Carbon dioxide, water vapor and methane soil efflux (soil respiration) in a Pinus palustris root exclusion in Georgetown, SC

This dataset contains processed data from a combination of survey flux chambers and long-term automated flux chambers. Soil flux measurements were conducted from June 2023 through December 2025 in a mature longleaf pine forest in Georgetown, SC. Soil respiration measurements were conducted approximately biweekly for two and a half years, before and after a root exclusion that took place on May 5, 2024. Processed, QAQC’d data for the treatment (root exclusion) and control (roots intact) before and after the root exclusion can be found in the file: 1_DATA_ESS_DOE_HR_RS_HB2_QAQC_Survey_Data_20260223.csv. Two multiday deployments were also conducted prior to the root exclusion using long-term automated chambers to continuously monitor greenhouse gas soil efflux. Processed, QAQC’d data for both long-term deployments can be found in the file: 2_DATA_ESS_DOE_HR_RS_HB2_QAQC_Longterm_Data_20260209.csv. Raw and working data files (.json, .81x, & .82z format) from LI-COR equipment are included for reference and can be accessed using SoilFluxPro software. CSV metadata files describe the raw data and modifications made using SoilFluxPro v5 and Matlab R2024b, as well as formatting and units for processed CSVs. Matlab code is included for reading in the processed CSVs, with sample figures comparing treatment and control. This research was performed as part of the project: “Improving models of stand and watershed carbon and water fluxes with more accurate representations of soil-plant-water dynamics in southern pine ecosystems”, which examines in part the effects hydraulic redistribution on soil efflux of carbon dioxide, water vapor and methane, as well as soil moisture and temperature in a southern pine ecosystem with sandy soils and high water table.

CARBON DIOXIDE FLUX↗

Carbon dioxide, water vapor and methane soil efflux (soil respiration) in a Pinus palustris restoration site in Georgetown, SC

This dataset contains processed data from a combination of survey flux chambers and long-term automated flux chambers. Biweekly soil flux measurements were conducted from June 2023 through December 2025 at a longleaf pine restoration site in Georgetown, SC. Processed, QAQC’d data can be found in the file: 1_DATA_ESS_DOE_HR_RS_HB3_QAQC_Survey_Data_20260223.csv. Raw and working data files (.json, & .81x format) from LI-COR equipment are included for reference and can be accessed using SoilFluxPro software. CSV metadata files describe the raw data and modifications made using SoilFluxPro v5 and Matlab R2024b, as well as formatting and units for processed CSVs. Matlab code is included for reading in the processed CSVs. This research was performed as part of the project: “Improving models of stand and watershed carbon and water fluxes with more accurate representations of soil-plant-water dynamics in southern pine ecosystems”, which examines in part the effects hydraulic redistribution on soil efflux of carbon dioxide, water vapor and methane, as well as soil moisture and temperature in a southern pine ecosystem with sandy soils and high water table.

CARBON DIOXIDE FLUX↗

Modeling the Impact of Riparian Hollows on River Corridor Nitrogen Exports

Recent studies in snowmelt-dominated catchments have documented changes in nitrogen (N) retention over time, such as declines in watershed exports of N, though there is a limited understanding of the controlling processes driving these trends. Working in the mountainous headwater East River Colorado watershed, we explored the effects of riparian hollows as N-cycling hotspots and as important small-scale controls on observed watershed trends. Using a modeling-based approach informed by remote sensing and in situ observations, we simulated the N-retention capacity of riparian hollows with seasonal and yearly hydrobiogeochemical perturbations imposed as drivers. We then implemented a scaling approach to quantify the relative contribution of riparian hollows to the total river corridor N budget. We found that riparian hollows primarily serve as N sinks, with N-transformation rates significantly limited by periods of enhanced groundwater upwelling and promoted at the onset of rainfall events. Given these observed hydrologic controls, we expect that the nitrate ( NO 3 - ) sink capacity of riparian hollows will increase in magnitude with future climatic perturbations, specifically the shift to more frequent rainfall events and fewer snowmelt events, as projected for many mountainous headwater catchments. Our current estimates suggest that while riparian hollows provision ~5–20% of NO 3 - to the river network, they functionally act as inhibitors to upland NO 3 - reaching the stream. Our work linking transient hydrological conditions to numerical biogeochemical simulations is an important step in assessing N-retaining features relative to the watershed N budget and better understanding the role of small-scale features within watersheds.

54 ENVIRONMENTAL SCIENCES↗

The application of remote sensing to the development and formulation of hydrologic planning models

Regional hydrologic planning models built upon remote sensing capabilities and suited for ungaged watersheds are developed. The effectiveness of such models is determined along with which parameters impact most the minimization of errors associated with the prediction of peak flow events (floods). Emphasis is placed on peak flood prediction because of its significance to users for the purpose of planning, sizing, and designing waterworks.

Castruccio, P. A.↗

Bridging Hydrological Ensemble Simulation and Learning Using Deep Neural Operators

Ensemble-based simulation and learning (ESnL) has long been used in hydrology for parameter inference, but computational demands of process-based ESnL can be quite high. To address this issue, we propose a deep neural operator learning approach. Neural operators are generic machine learning algorithms that can learn functional mappings between infinite-dimensional spaces, providing a highly flexible tool for scientific machine learning. Our approach is built upon DeepONet, a specific deep neural operator, and is designed to address several common problems in hydrology, namely, model parameter estimation, prediction at ungaged locations, and uncertainty quantification. Here we demonstrate the effectiveness of our DeepONet-based workflow using an existing large model ensemble created for an eastern U.S. watershed that is instrumented with 10 streamflow gages. Results suggest DeepONet achieves high efficiency in learning an ML surrogate model from the model ensemble, with the modified Kling-Gupta Efficiency exceeding 0.9 on holdout test sets. Parameter inference, carried out using the trained DeepONet surrogate model and genetic algorithm, also yields robust results. Additionally, we formulate and train a separate DeepONet model for physics-informed, seq-to-seq streamflow forecasting, which further reduces biases in the pre-trained DeepONet surrogate model. While this study focuses primarily on a single watershed, our approach is general and may be extended to enable learning from model ensembles across multiple basins or models. Thus, this research represents a significant contribution to the application of hybrid machine learning in hydrology.

54 ENVIRONMENTAL SCIENCES↗

Optimizing process-based models to predict current and future soil organic carbon stocks at high-resolution

From hillslope to small catchment scales (< 50 km 2 ), soil carbon management and mitigation policies rely on estimates and projections of soil organic carbon (SOC) stocks. Here we apply a process-based modeling approach that parameterizes the MIcrobial-MIneral Carbon Stabilization (MIMICS) model with SOC measurements and remotely sensed environmental data from the Reynolds Creek Experimental Watershed in SW Idaho, USA. Calibrating model parameters reduced error between simulated and observed SOC stocks by 25%, relative to the initial parameter estimates and better captured local gradients in climate and productivity. The calibrated parameter ensemble was used to produce spatially continuous, high-resolution (10 m 2 ) estimates of stocks and associated uncertainties of litter, microbial biomass, particulate, and protected SOC pools across the complex landscape. Here, subsequent projections of SOC response to idealized environmental disturbances illustrate the spatial complexity of potential SOC vulnerabilities across the watershed. Parametric uncertainty generated physicochemically protected soil C stocks that varied by a mean factor of 4.4 × across individual locations in the watershed and a – 14.9 to + 20.4% range in potential SOC stock response to idealized disturbances, illustrating the need for additional measurements of soil carbon fractions and their turnover time to improve confidence in the MIMICS simulations of SOC dynamics.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Effectiveness of Soil Profile Rehabilitation for Pluvial Flood Mitigation Through Two-Dimensional Hydrodynamic Modeling

Pluvial flooding, driven by increasingly impervious surfaces and intense storm events, presents a growing challenge for urban areas worldwide. In Baltimore City, MD, USA, climate change, rapid urbanization, and aging stormwater infrastructure are exacerbating flooding impacts, resulting in significant socio-economic consequences. This study evaluated the effectiveness of a soil profile rehabilitation scenario using a 2D hydrodynamic modeling approach for the Tiffany Run watershed, Baltimore City. This study utilized different extreme storm events, a high-resolution (1 m) LiDAR Digital Terrain Model (DTM), building footprints, and hydrological soil data. These datasets were integrated into a fully coupled 2D hydrodynamic model, the City Catchment Analysis Tool (CityCAT), to simulate urban flood dynamics. The pre-soil rehabilitation simulation revealed a maximum water depth of 3.00 m in most areas, with hydrologic soil groups C and D, especially downstream of the study area. The post-soil rehabilitation simulation was targeted at vacant lots and public parcels, accounting for 33.20% of the total area of the watershed. This resulted in a reduced water depth of 2.50 m. Additionally, the baseline runoff coefficient of 0.49 decreased to 0.47 following the rehabilitation, and the model consistently recorded a peak runoff reduction rate of 4.10 across varying rainfall intensities. The validation using a contingency matrix demonstrated true-positive rates of 0.75, 0.50, 0.64, and 0 for the selected events, confirming the model’s capability at capturing real-world flood occurrences.

Baltimore City↗

Remote sensing techniques for prediction of watershed runoff

Hydrologic parameters of watersheds for use in mathematical models and as design criteria for flood detention structures are sometimes difficult to quantify using conventional measuring systems. The advent of remote sensing devices developed in the past decade offers the possibility that watershed characteristics such as vegetative cover, soils, soil moisture, etc., may be quantified rapidly and economically. Experiments with visible and near infrared data from the LANDSAT-1 multispectral scanner indicate a simple technique for calibration of runoff equation coefficients is feasible. The technique was tested on 10 watersheds in the Chickasha area and test results show more accurate runoff coefficients were obtained than with conventional methods. The technique worked equally as well using a dry fall scene. The runoff equation coefficients were then predicted for 22 subwatersheds with flood detention structures. Predicted values were again more accurate than coefficients produced by conventional methods.

Blanchard, B. J.↗

Modeling Riverine Dissolved and Particulate Organic Carbon Fluxes from Two Small Watersheds in the Northeastern United States

The coupled carbon (C) cycle across terrestrial and aquatic environments at the watershed scale has been identified as an important, but poorly constrained component of the global carbon budget. Here, we extended Soil and Water Assessment Tool (SWAT) with coupled riverine particulate organic carbon (POC) and dissolved organic carbon (DOC) modules (referred to as SWAT-C hereafter). Results show that SWAT-C reproduced daily POC and DOC fluxes well in two watersheds in the Northeastern United States. We found that SWAT-C tended to underestimate high flows and peak DOC and POC fluxes. Uncertainty analysis indicated flux uncertainties associated with POC and DOC simulation were larger than those for flow simulation. Sensitive parameters controlling POC and DOC biogeochemical processes were identified along with how these parameters influence mechanisms underlying C cycling. We anticipate that the tool developed and applied here will inform C related ecosystem services in watershed assessment and planning.

Dissolved organic carbon, Particulate organic carb↗