SEARCH · Engineering Papers

Results for “integrated hydrologic models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Disentangling atmospheric, hydrological, and coupling uncertainties in compound flood modeling within a coupled Earth system model

Compound riverine and coastal flooding is usually driven by complex interactions among meteorological, hydrological, and ocean extremes. However, existing efforts to model this phenomenon often do not integrate hydrological processes across atmosphere–land–river–ocean systems, leading to substantial uncertainties that have not been fully examined. To bridge this gap, we leverage the new capabilities of the Energy Exascale Earth System Model (E3SM) that enable a multi-component framework that integrates coastally refined atmospheric, terrestrial, and oceanic components. We evaluate compound uncertainties arising from two-way land–river–ocean coupling in E3SM and track the cascading meteorological and hydrological uncertainties through ensemble simulations over the Delaware River basin and estuary during Hurricane Irene (2011). Our findings highlight the importance of two-way river–ocean coupling to compound flood modeling and demonstrate E3SM's capability in capturing compound flood extent near the coast, with a hit rate over 0.75. Our study shows the growing uncertainties that transition from atmospheric forcings to flood distribution and severity. Furthermore, an analysis based on artificial neural networks is used to assess the roles of hydrological drivers, such as infiltration and soil moisture, in the generation of compound flooding. The response of compound floods to tropical cyclones (TCs) is found to be susceptible to these often overlooked drivers. For instance, the flooded area could increase more than 2-fold (∼2.4) if Hurricane Irene were preceded by an extreme antecedent soil moisture condition (AMC). The results not only support the use of a multi-component framework for interactive flooding processes, but also underscore the necessity of broader definitions of compound flooding that encompass the simultaneous occurrence of intense precipitation, storm surge, and high AMC during TCs.

Feng, Dongyu [Pacific Northwest National Laborator↗

NASA NTRS2007

Assimilation of GRACE Terrestrial Water Storage Data into a Land Surface Model: Results for the Mississippi River Basin

NASA's GRACE mission has the potential to be extremely valuable for water resources applications and global water cycle research. What makes GRACE unique among Earth Science satellite systems is that it is able to monitor variations in water stored in all forms, from snow and surface water to soil moisture to groundwater in the deepest aquifers. However, the space and time resolutions of GRACE observations are coarse. GRACE typically resolves water storage changes over regions the size of Nebraska on a monthly basis, while city-scale, daily observations would be more useful for water management, agriculture, and weather prediction. High resolution numerical (computer) hydrology models have been developed, which predict the fates of water and energy after they strike the land surface as precipitation and sunlight. These are similar to weather and climate forecast models, which simulate atmospheric processes. We integrated the GRACE observations into a hydrology model using an advanced technique called data assimilation. The results were new estimates of groundwater, soil moisture, and snow variations, which combined the veracity of GRACE with the high resolution of the model. We tested the technique over the Mississippi River basin, but it will be even more valuable in parts of the world which lack reliable data on water availability.

Zaitchik, Benjamin F.↗

DOE OSTI2024

Watching the Grand Ethiopian Renaissance Dam from a distance: Implications for sustainable water management of the Nile water

Increased demands for sustainable water and energy resources in densely populated basins have led to the construction of dams, which impound waters in artificial reservoirs. In many cases, scarce field data led to the development of models that underestimated the seepage losses from reservoirs and ignored the role of extensive fault networks as preferred pathways for groundwater flow. We adopt an integrated approach (remote sensing, hydrologic modeling, and field observations) to assess the magnitude and nature of seepage from such systems using the Grand Ethiopian Renaissance Dam (GERD), Africa's largest hydropower project, as a test site. The dam was constructed on the Blue Nile within steep, highly fractured, and weathered terrain in the western Ethiopian Highlands. The GERD Gravity Recovery and Climate Experiment Terrestrial Water Storage (GRACETWS), seasonal peak difference product, reveals significant mass accumulation (43 ± 5 BCM) in the reservoir and seepage in its surroundings with progressive south-southwest mass migration along mapped structures between 2019 and 2022. Seepage, but not a decrease in inflow or increase in outflow, could explain, at least in part, the observed drop in the reservoir's water level and volume following each of the three fillings. Using mass balance calculations and GRACETWS observations, we estimate significant seepage (19.8 ± 6 BCM) comparable to the reservoir's impounded waters (19.9 ± 1.2 BCM). Investigating and addressing the seepage from the GERD will ensure sustainable development and promote regional cooperation; overlooking the seepage would compromise hydrological modeling efforts on the Nile Basin and misinform ongoing negotiations on the Nile water management.

GRACE and GRACE-FO↗

NASA NTRS2023

A Multi-Satellite Framework to Rapidly Evaluate Extreme Biosphere Cascades: The Western US 2021 Drought and Heatwave

The increasing frequency and intensity of climate extremes and complex ecosystem responses motivate the need for integrated observational studies at low-latency to determine biosphere responses and carbon-climate feedbacks. Here, we develop a satellite-based rapid attribution workflow and demonstrate its use at a 1–2-month latency to attribute drivers of the carbon cycle feedbacks during the 2020-2021 Western US drought and heatwave. In the first half of 2021, concurrent negative photosynthesis anomalies and large positive column CO 2 anomalies were detected with satellites. Using a simple atmospheric mass balance approach, we estimate a surface carbon efflux anomaly of 132 TgC in June 2021, a magnitude corroborated 28 independently with a dynamic global vegetation model. Integrated satellite observations of hydrologic processes, representing the soil-plant-atmosphere continuum (SPAC), show that these surface carbon flux anomalies are largely due to substantial reductions in photosynthesis because of a spatially widespread moisture-deficit propagation through the SPAC between 2020 and 2021. A causal model indicates deep soil moisture stores partially drove photosynthesis, maintaining its values in 2020 and driving its declines throughout 2021. The causal model also suggests legacy effects may have amplified photosynthesis deficits in 2021 beyond the direct effects of environmental forcing. The integrated, observation framework presented here provides a valuable first assessment of a biosphere extreme response and an independent testbed for improving drought propagation and mechanisms in models. The rapid identification of extreme carbon anomalies and hotspots can also aid mitigation and adaptation decisions.

Causal model↗

DOE OSTI2021

DPC Disposal Thermal Scoping Analysis

This is a progress report on thermal modeling for dual-purpose canister (DPCs) direct disposal that covers several available calculation methods and addresses creep and temperature-dependent properties in a salt repository. Three modeling approaches are demonstrated: A semi-analytical calculation method that uses linear solutions with superposition and imaging, to represent a central waste package in a larger array; A finite difference model of coupled thermal creep, implemented in FLAC2D; and An integrated finite difference thermal-hydrologic modeling approach for repositories in different generic host media, implemented in PFLOTRAN. These approaches are at different levels of maturity, and future work is expected to add refinements and establish the best applications for each.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Data and scripts associated with the manuscript "A Hydrogeophysical Framework to Assess Infiltration During a Simulated Ecosystem-Scale Flooding Experiment”

This package contains the data and scripts used in "A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment" (Adebayo et al., 2023). The Geophysics folder contains two subfolders containing Ground Penetrating Radar (GPR) and Electrical Resistivity Imaging (ERI) datasets acquired between August 23 and August 26, 2021. The GPR subfolder contains raw data (Tempest_Project.GPZ) and processed results using the Sensors and Software EKKO Project software. The ERI subfolder contains two additional subfolders named Static and Timelapse. The Static subfolder contains background resistivity datasets acquired before infiltration, both raw and inverted results using the ResIpy code. The Timelapse subfolder contains repeated resistivity datasets acquired during infiltration on two profiles, along with raw and inverted results using the ResIpy code. The Flow model folder contains ATS model simulation files. The Geophysics-informed subfolder contains simulation datasets for two profiles generated using information from the geophysics data. The National_database-informed subfolder contains simulation datasets for two profiles generated using information from a national soil database. Each subfolder contains initial simulation results, multiple ensembles, and final simulation files. The Sensitivity_analysis folder contains the observation files exported from the flow model, notebooks used in running the code, and sensitivity analysis results showing the sensitivity of different parameters over a 48-hour period. We used mutual information-based global sensitivity analysis for the sensitivity analysis, following Jiang et al., 2022. This study shows a novel approach that integrates geophysical data with hydrological modeling to understand subsurface flow during coastal flooding.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Decision support for United States—Canada energy integration is impaired by fragmentary environmental and electricity system modeling capacity

The renewable energy transition is leading to increased electricity trade between the United States and Canada, with Canadian hydropower providing firm lower-carbon power and buffering variability of wind and solar generation in the U.S. However, long-term power purchase agreements and transborder transmission projects are controversial, with two of four proposed transmission lines between Quebec, Canada and the northeast U.S. cancelled since 2018. Here, we argue that controversies are exacerbated by a lack of open-source data and tools to understand tradeoffs of new hydropower generation and transmission infrastructure in comparison to alternatives. This gap includes impacts that incremental transmission and generation projects have on the economics of the entire system, for example, how new transmission projects affect exports to existing markets or incentivize new generation. We identify priority areas for data synthesis and model development, such as integrating linked hydropower and hydrologic interactions in energy system models and openly releasing (by utilities) or back-calculating (by researchers) hydropower generation and operational parameters. Publicly available environmental (e.g. streamflow, precipitation) and techno-economic (e.g. costs, reservoir size,) data can be used to parameterize freely usable and extensible models. Existing models have been calibrated with operational data from Canadian utilities that are not publicly available, limiting the range of scientific and commercial questions these tools have been used to answer and the range of parties that have been involved. Studies conducted using highly resolved, national-scale public data exist in other countries, notably, the United States, and demonstrate how greater transparency and extensibility can drive industry action. Improved data availability in Canada could facilitate approaches that (1) increase participation in decarbonization planning by a broader range of actors; (2) allow independent characterizations of environmental, health, and economic outcomes of interest to the public; and (3) identify decarbonization pathways consistent with community values.

13 HYDRO ENERGY↗

NASA NTRS2010

Utilization of Hydrologic Remote Sensing Data in Land Surface Modeling and Data Assimilation: Current Status and Challenges

Recent advances in remote sensing technologies have enabled the monitoring and measurement of the Earth's land surface at an unprecedented scale and frequency. The myriad of these land surface observations must be integrated with the state-of-the-art land surface model forecasts using data assimilation to generate spatially and temporally coherent estimates of environmental conditions. These analyses are of critical importance to real-world applications such as agricultural production, water resources management and flood, drought, weather and climate prediction. This need motivated the development of NASA Land Information System (LIS), which is an expert system encapsulating a suite of modeling, computational and data assimilation tools required to address challenging hydrological problems. LIS integrates the use of several community land surface models, use of ground and satellite based observations, data assimilation and uncertainty estimation techniques and high performance computing and data management tools to enable the assessment and prediction of hydrologic conditions at various spatial and temporal scales of interest. This presentation will focus on describing the results, challenges and lessons learned from the use of remote sensing data for improving land surface modeling, within LIS. More specifically, studies related to the improved estimation of soil moisture, snow and land surface temperature conditions through data assimilation will be discussed. The presentation will also address the characterization of uncertainty in the modeling process through Bayesian remote sensing and computational methods.

Kumar, Sujay V.↗

DOE OSTI2026

Flood Estimation under Snowmelt and Rain-on-Snow Processes in Alaska: A Military Installation Perspective

Accurate flood estimation in snow-dominated and high-latitude regions remains challenging due to complex interactions among rainfall, snowmelt, and rain-on-snow (ROS) processes, which are not captured in conventional precipitation-based intensity-duration-frequency (PREC-IDF) curves. This study evaluates the Next-Generation IDF (NG-IDF) framework, an extension of PREC-IDF that incorporates total water available for runoff (precipitation minus changes in snow water equivalent), in two contrasting Alaskan watersheds of Department of Defense (DoD) significance: the Little Chena River Basin (LCRB) in interior Alaska, a tributary of the Chena River that flows through Fort Wainwright and near Eielson Air Force Base, and the Upper Ship Creek Basin (USCB), which drains Joint Base Elmendorf-Richardson near Anchorage. Using long-term SNOTEL observations and event-based rainfall-runoff modeling, NG-IDF and PREC-IDF flood estimates were compared against observation-based flood frequency analyses. Results show that NG-IDF consistently reduces flood-estimation bias by 15–20% relative to PREC-IDF, particularly for snowmelt- and ROS-dominated events. In the interior LCRB, permafrost conditions can substantially amplify flood responses during snowmelt events, an effect not explicitly represented in standard design tools. These findings demonstrate that NG-IDF provides a more physically consistent and transferable framework for flood estimation in cold regions, with potential relevance to mission-critical DoD installation resilience. Projected increases in ROS frequency and permafrost degradation across Alaska further emphasize the need to integrate physics-based hydrologic models that explicitly represent snow and permafrost processes to enhance design resilience and operational readiness for military and other critical infrastructure.

Yan, Hongxiang (ORCID:000000022387403X)↗

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↗

DOE OSTI2020

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Advanced Terrestrial Simulator (ATS) evaluation dataset at 7 catchments across the continental United States

This dataset comprises of the input files and other files required for Advanced Terrestrial Simulator (ATS) simulations at 7 catchments across the continental United States. ATS is an integrated surface-subsurface hydrology model. We include Jupyter notebooks (within scripts folder) for individual catchments showing information (including data sources, river network, soil, geology, landuse types etc.) on preparing the machine readable input files. ATS observation output files are provided in the output folder. Figures and analyses (.xlsx sheets) are also provided. The catchments include, Taylor River Upstream (Colorado); (b) Cossatot River (Arkansas); (c) Panther Creek (Alabama); (d) Little Tennessee River (North Carolina and Georgia); (e) Mayo River (Virginia); (f) Flat Brook (New Jersey); (g) Neversink River headwaters (New York). Readme files are provided inside the directories providing more details. Files types include: .xml, .h5, .xlsx, .png, .ipynb, .py, .nc, .txt. All of the files types can be accessed by open source software, details on software requirements are following: .xml (any text editors including notepad and textedit), .h5 (in python using hdf libraries), .xlsx (WPS Office Spreadsheets, OpenOffice Calc, LibreOffice Calc, Microsoft Office etc.), .png (any image viewer), .ipynb (Jupyter notebook), .py (any text editors including notepad and textedit), .nc (using python or other open source software).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Discrete global grid system-based flow routing datasets in the Amazon and Yukon basins

Abstract. Discrete global grid systems (DGGS) are emerging spatial data structures widely used to organize geospatial datasets across scales. While DGGS have found applications in various scientific disciplines, including atmospheric science and ecology, their integration into physically based hydrological models and Earth system models (ESMs) has been hindered by the lack of flow routing datasets based on DGGS. In response to this gap, this study pioneers the development of new flow routing datasets using icosahedral Snyder equal-area (ISEA) DGGS and a novel mesh-independent flow direction model. We present flow routing datasets for two large basins, the tropical Amazon River basin and the Arctic Yukon River basin. These datasets (1) facilitate the adoption of DGGS for hydrological models and (2) provide flow routing inputs for evaluation of DGGS-based flow routing in the Amazon and Yukon river basins. The data are available at https://doi.org/10.5281/zenodo.8377765 (Liao, 2023).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Estuarine hurricane wind can intensify surge-dominated extreme water level in shallow and converging coastal systems

Based on the projected increase in hurricane landfall frequency on the middle to lower US east coast, we examined the crucial role of the estuarine wind field in exacerbating coastal flooding. A regionally refined atmospheric and two high-resolution hydrology and ocean models are integrated to provide plausible and physically consistent ensembles of hurricane events and the associated flooding inside the Delaware Bay and River, a US mid-Atlantic estuary. Model results show that the hurricane propagation direction, estuarine geometry, remote surge from the open ocean, and direct nearshore upwind stress could magnify the flood magnitude. More specifically, inland-bound tracks that make landfall before reaching the mid-Atlantic coast produce a more significant surge within Delaware Bay than the shore-parallel tracks, where the estuarine wind direction plays the primary role in surge amplification. Ultimately, this study emphasized the need for integrated models to capture the nonlinear dynamics and interactions in flood hazard modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Ambient field seismology in critical zone hydrological sciences

Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.

58 GEOSCIENCES↗

DOE OSTI2026

A coupled hydrologic-agroeconomic modeling framework to evaluate adaptive irrigation strategies under groundwater withdrawal restrictions

Growing groundwater scarcity requires integrated tools to capture interactions among hydrology, agricultural production, markets, and land use. This study presents an iterative modeling framework that couples hydrologic, crop-yield, and economic models to capture two-way feedback among water availability, agricultural production, and market responses under groundwater constraints. The primary goal of this paper is to describe the methodological development of the coupled framework and demonstrate the significance of iterative model interaction. Applied to the western United States, we evaluated adaptive responses to restricting groundwater use beyond recharge levels, represented through changes in irrigation management and expansion or shrinkage of crop markets through land reallocation. Results demonstrate that the iterative coupling converges to stable equilibrium responses within 10 iterations. At equilibrium, deficit irrigation emerges as the dominant adaptation strategy in California, with irrigation levels stabilizing at approximately 70% of full irrigation demand, while Arizona and New Mexico experience stronger yield sensitivities. Early iterations produce commodity price increases of up to 10% for fruit and vegetable crops; however, these responses moderate as land allocation and production patterns adjust across regions. Deficit irrigation and spatial reallocation of irrigated land partially offset production losses, with variability observed across different states: California maintains yields primarily via deficit irrigation, whereas Arizona and New Mexico will rely mainly on reducing irrigated area to absorb the shock. By capturing feedback between biophysical and economic processes, this approach highlights how irrigation strategies and land-use decisions evolve under water stress and provides a transferable platform for evaluating water management policies.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2022

Tropical peatland hydrology simulated with a global land surface model

Tropical peatlands are among the most carbon-dense ecosystems on Earth, and their water storage dynamics strongly control these carbon stocks. The hydrological functioning of tropical peatlands differs from that of northern peatlands, which has not yet been accounted for in global land surface models (LSMs). Here, we integrated tropical peat-specific hydrology modules into a global LSM for the first time, by utilizing the peatland-specific model structure adaptation (PEATCLSM) of the NASA Catchment Land Surface Model (CLSM). We developed literature-based parameter sets for natural (PEATCLSM(Trop,Nat)) and drained (PEATCLSM(Trop,Drain)) tropical peatlands. Simulations with PEATCLSM(Trop,Nat) were compared against those with the default CLSM version and the northern version of PEATCLSM (PEATCLSM(North,Nat)) with tropical vegetation input. All simulations were forced with global meteorological reanalysis input data for the major tropical peatland regions in Central and South America, the Congo Basin, and Southeast Asia. The evaluation against a unique and extensive data set of in situ water level and eddy covariance-derived evapotranspiration showed an overall improvement in bias and correlation compared to the default CLSM version. Over Southeast Asia, an additional simulation with PEATCLSM(Trop,Drain) was run to address the large fraction of drained tropical peatlands in this region. PEATCLSM(Trop,Drain) outperformed CLSM, PEATCLSM(North,Nat) and PEATCLSM(Trop,Nat) over drained sites. Despite the overall improvements of PEATCLSM(Trop,Nat) over CLSM, there are strong differences in performance between the three study regions. We attribute these performance differences to regional differences in accuracy of meteorological forcing data, and differences in peatland hydrologic response that are not yet captured by our model.

Peatland↗

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↗