Search NASA⌕ Search

SEARCH · Search NASA

Results for “Streamflow”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Changing Snow Cover and Stream Discharge in the Western United States - Wind River Range, Wyoming

Earlier onset of springtime weather has been documented in the western United States over at least the last 50 years. Because the majority (>70%) of the water supply in the western U.S. comes from snowmelt, analysis of the declining spring snowpack has important implications for the management of water resources. We studied ten years of Moderate-Resolution Imaging Spectroradiometer (MODIS) snow-cover products, 40 years of stream discharge and meteorological station data and 30 years of snow-water equivalent (SWE) SNOw Telemetry (SNOTEL) data in the Wind River Range (WRR), Wyoming. Results show increasing air temperatures for.the 40-year study period. Discharge from streams in WRR drainage basins show lower annual discharge and earlier snowmelt in the decade of the 2000s than in the previous three decades. Changes in streamflow may be related to increasing air temperatures which are probably contributing to a reduction in snow cover, although no trend of either increasingly lower streamflow or earlier snowmelt was observed within the decade of the 2000s. And SWE on 1 April does not show an expected downward trend from 1980 to 2009. The extent of snow cover derived from the lowest-elevation zone of the WRR study area is strongly correlated (r=0.91) with stream discharge on 1 May during the decade of the 2000s. The strong relationship between snow cover and streamflow indicates that MODIS snow-cover maps can be used to improve management of water resources in the drought-prone western U.S.

Hall, Dorothy K.↗
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↗
NASA NTRS2015

Hazard Assessment from Storm Tides and Rainfall on a Tidal River Estuary

Here, we report on methods and results for a model-based flood hazard assessment we have conducted for the Hudson River from New York City to Troy/Albany at the head of tide. Our recent work showed that neglecting freshwater flows leads to underestimation of peak water levels at up-river sites and neglecting stratification (typical with two-dimensional modeling) leads to underestimation all along the Hudson. As a result, we use a three-dimensional hydrodynamic model and merge streamflows and storm tides from tropical and extratropical cyclones (TCs, ETCs), as well as wet extratropical cyclone (WETC) floods (e.g. freshets, rain-on-snow events). We validate the modeled flood levels and quantify error with comparisons to 76 historical events. A Bayesian statistical method is developed for tropical cyclone streamflows using historical data and consisting in the evaluation of (1) the peak discharge and its pdf as a function of TC characteristics, and (2) the temporal trend of the hydrograph as a function of temporal evolution of the cyclone track, its intensity and the response characteristics of the specific basin. A k-nearest-neighbors method is employed to determine the hydrograph shape. Out of sample validation tests demonstrate the effectiveness of the method. Thus, the combined effects of storm surge and runoff produced by tropical cyclones hitting the New York area can be included in flood hazard assessment. Results for the upper Hudson (Albany) suggest a dominance of WETCs, for the lower Hudson (at New York Harbor) a case where ETCs are dominant for shorter return periods and TCs are more important for longer return periods (over 150 years), and for the middle-Hudson (Poughkeepsie) a mix of all three flood events types is important. However, a possible low-bias for TC flood levels is inferred from a lower importance in the assessment results, versus historical event top-20 lists, and this will be further evaluated as these preliminary methods and results are finalized. Future funded work will quantify the influences of sea level rise and flood adaptation plans (e.g. surge barriers). It would also be valuable to examine how streamflows from tropical cyclones and wet cool-season storms will change, as this factor will dominate at upriver locations.

Hazard assessment↗
NASA NTRS2016

Multi-Model Assessment of Global Hydropower and Cooling Water Discharge Potential Under Climate Change

Worldwide, 98% of total electricity is currently produced by thermoelectric power and hydropower. Climate change is expected to directly impact electricity supply, in terms of both water availability for hydropower generation and cooling water usage for thermoelectric power. Improved understanding of how climate change may impact the availability and temperature of water resources is therefore of major importance. Here we use a multi-model ensemble to show the potential impacts of climate change on global hydropower and cooling water discharge potential. For the first time, combined projections of streamflow and water temperature were produced with three global hydrological models (GHMs) to account for uncertainties in the structure and parametrization of these GHMs in both water availability and water temperature. The GHMs were forced with bias-corrected output of five general circulation models (GCMs) for both the lowest and highest representative concentration pathways (RCP2.6 and RCP8.5). The ensemble projections of streamflow and water temperature were then used to quantify impacts on gross hydropower potential and cooling water discharge capacity of rivers worldwide. We show that global gross hydropower potential is expected to increase between +2.4% (GCM-GHM ensemble mean for RCP 2.6) and +6.3% (RCP 8.5) for the 2080s compared to 1971-2000. The strongest increases in hydropower potential are expected for Central Africa, India, central Asia and the northern high-latitudes, with 18-33% of the world population living in these areas by the 2080s. Global mean cooling water discharge capacity is projected to decrease by 4.5-15% (2080s). The largest reductions are found for the United States, Europe, eastern Asia, and southern parts of South America, Africa and Australia, where strong water temperature increases are projected combined with reductions in mean annual streamflow. These regions are expected to affect 11-14% (for RCP2.6 and the shared socioeconomic pathway (SSP)1, SSP2, SSP4) and 41-51% (RCP8.5-SSP3, SSP5) of the world population by the 2080s.

global hydrological models↗
NASA NTRS2022

Assimilation of NASA’s Airborne Snow Observatory Snow Measurements for Improved Hydrological Modeling: A Case Study Enabled by the Coupled LIS/WRF-Hydro System

he NASA LIS/WRF-Hydro system is a coupled modeling framework that combines the modeling and data assimilation (DA) capabilities of the NASA Land Information System (LIS) with the multi-scale surface hydrological modeling capabilities of the WRF-Hydro model, both of which are widely used in both operations and research. This coupled modeling framework builds on the linkage between land surface models (LSMs), which simulate surface boundary conditions in atmospheric models, and distributed hydrologic models, which simulate horizontal surface and sub-surface flow, adding new land DA capabilities. In the present study, we employ this modeling framework in the Tuolumne River basin in central California. We demonstrate the added value of the assimilation of NASA Airborne Snow Observatory (ASO) snow water equivalent (SWE) estimates in the Tuolumne basin. This analysis is performed in both LIS as an LSM column model and LIS/WRF-Hydro, with hydrologic routing. Results demonstrate that ASO DA in the basin reduced snow bias by as much as 30% from an open-loop (OL) simulation compared to three independent datasets. It also reduces downstream streamflow runoff biases by as much as 40%, and improves streamflow skill scores in both wet and dry years. Analysis of soil moisture and evapotranspiration (ET) also reveals the impacts of hydrologic routing from WRF-Hydro in the simulations, which would otherwise not be resolved in an LSM column model. By demonstrating the beneficial impact of SWE DA on the improving streamflow forecasts, the article outlines the importance of such observational inputs for reservoir operations and related water management applications.

Hydrologic modeling↗
NASA NTRS2022

Assimilation of NASA's Airborne Snow Observatory Snow Measurements for Improved Hydrological Modeling: A Case Study Enabled by the Coupled LIS/WRF-Hydro System

The NASA LIS/WRF-Hydro system is a coupled modeling framework that combines the modeling and data assimilation (DA) capabilities of the NASA Land Information System (LIS) with the multi-scale surface hydrological modeling capabilities of the WRF-Hydro model, both of which are widely used in both operations and research. This coupled modeling framework builds on the linkage between land surface models (LSMs), which simulate surface boundary conditions in atmospheric models, and distributed hydrologic models, which simulate horizontal surface and sub-surface flow, adding new land DA capabilities. In the present study, we employ this modeling framework in the Tuolumne River basin in central California. We demonstrate the added value of the assimilation of NASA Airborne Snow Observatory (ASO) snow water equivalent (SWE) estimates in the Tuolumne basin. This analysis is performed in both LIS as an LSM column model and LIS/WRF-Hydro, with hydrologic routing. Results demonstrate that ASO DA in the basin reduced snow bias by as much as 30% from an open-loop (OL) simulation compared to three independent datasets. It also reduces downstream streamflow runoff biases by as much as 40%, and improves streamflow skill scores in both wet and dry years. Analysis of soil moisture and evapotranspiration (ET) also reveals the impacts of hydrologic routing from WRF-Hydro in the simulations, which would otherwise not be resolved in an LSM column model. By demonstrating the beneficial impact of SWE DA on the improving streamflow forecasts, the article outlines the importance of such observational inputs for reservoir operations and related water management applications.

Timothy M. Lahmers↗
NASA NTRSDate not supplied

Economic Impact Assessments (EIA) of application of GEOGLOWS in Ecuador: Data Gaps, Limitations and Recommendations

In 2022, the United Nations launched the Early Warnings for All (EW4ALL) Program to establish global early warning systems by 2027. To assess the impact of the substantial $3.1 billion annual investment over five years, EW4ALL will consider factors that will require national coordination for the data needed for these assessments. In 2023, Ecuador was identified as one of the world's most climate-vulnerable countries, emphasizing the need to enhance its early warning systems. In 2020, the SERVIR Amazonia hub implemented the GEOGLOWS streamflow forecast service in collaboration with Ecuador's national meteorological agency (INAMHI). GEOGloWS provides 15-day ensemble forecasts and 80 years of historical streamflow data for every river worldwide through a free web service. The World Meteorological Organization has recognized this initiative as essential in contributing to the UN's call to ensure an 'Early Warning for All' by 2027. In 2023, as part of NASA's continuous efforts to fund research for Policy-Relevant Implementations, an economic impact assessment (EIA) was performed to understand the potential socioeconomic benefits of Early streamflow predictions in Ecuador using the GEOGLOWS service. Preliminary findings highlighted that gaps remain in effectively integrating socioeconomic and Earth observation (EO) data to capture the total value of these predictions. Implementing GEOGLOWS has led to valuable hydrological forecasts; however, the total economic benefits have yet to be documented. This study addresses the gaps and makes recommendations for future work that should focus on capturing the socio-economic benefits and costs associated with these forecasts, including their impact on decision-making at national and local levels. Despite the daily use of GEOGLOWS by key figures, including the President of Ecuador, the need for comprehensive recommendations and assessments is urgent.

Reetwika Basu↗
DOE OSTI2024

JAMES BUTTLE REVIEW: Interflow, subsurface stormflow and throughflow: A synthesis of field work and modelling

Interflow, throughflow and subsurface stormflow are interchangeable terms that refer to the lateral subsurface flow above a restricting layer of lower hydraulic conductivity that occurs during and following storm events. Interflow (used here) is a more dominant process in steeper catchments with high infiltration capacity soils overlying a more impermeable soil or geologic layer. Interflow as a runoff process was first recognised in the early 1900s, yet hydrologists still struggle to predict its occurrence, persistence, importance, interaction with other streamflow generation processes, and potential to connect to valleys and streams during and following storms. We review the history of interflow research and address some of the challenges in understanding its role in runoff production. We argue that characterising the controls on interflow initiation and occurrence relies on detailed field observations of subsurface properties, which exist only in limited experimental settings. This data shortcoming contributes to our inability to predict interflow or determine its contribution to streamflow more broadly. There remain many opportunities to advance our understanding of interflow that include both modelling and experimental or observational approaches in hydrology.

hillslope hydrology↗
DOE OSTI2024

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↗
DOE OSTI2025

Data and scripts associated with “When do Riverine Systems 'Feel the Burn'? Simulating How Burn Extent and Severity Modulate Hydrologic Controls on Biogeochemical Export” (v2)

This data package is associated with the publication “When do Riverine Systems 'Feel the Burn'? Simulating How Burn Extent and Severity Modulate Hydrologic Controls on Biogeochemical Export” published in Water Resources Research (Wampler et al. 2025; preprint: https://doi.org/10.22541/essoar.174438106.63564767/v1). This study used the Soil and Water Assessment Tool (SWAT), a processed based model to explore the impacts of area burned and burn severity on streamflow, nitrate, and dissolved organic carbon (DOC) in two test basins: a semi-arid, mixed land use basin and a humid, primarily forested basin. We developed 1800 wildfire scenarios that we ran in each basin: 20 different burn extents (5 to 100% by 5%), 3 different burn severities (low, moderate, and high), and 30 different post-fire precipitation scenarios. We also ran an additional 30 scenarios associated with no wildfire for the 30 post-fire precipitation scenarios. For each scenario we were interested in the change in runoff ratio (streamflow) and average concentration and annual loads (nitrate and DOC) across the wildfire scenarios. This data package contains the data and scripts required to build SWAT models for the two test basins, create and run the wildfire scenarios, and generate the data summaries and figures used in the associated manuscript. This data package was originally published in March 2025. It was updated in January 2026 (v2; new and modified files) to include the final files after the manuscript went through reviews. See the change history section below for more details. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2018

Ground and Satellite Based Observation Datasets for the Lower Mekong River Basin

In ‘Satellite observations and modeling to understand the Lower Mekong River Basin streamflow variability’ [1] hydrological fluxes, meteorological variables, land cover land use maps, and soil characteristics and parameters data were compiled and processed for the Lower Mekong River Basin. In this work, daily streamflow time series data at nine gauges located at five different countries in the Mekong region (Thailand, Laos People׳s Democratic Republic (PDR), Myanmar, Cambodia, and Viet Nam) is presented. Satellite-based daily precipitation and air temperature (minimum & maximum) data is processed and provided over the entire basin as part of the dataset provided in this work. Moreover, land cover land use raster data that contains 18 classes that cover agriculture, urban, range and forests land cover land use classes for the basin is offered. In addition, a soil data that contains physical and chemical characteristics needed by physically based hydrological models to simulate the cycling of water and air is also provided.

Streamflow↗
NASA NTRS2023

Towards an Indian Land Data Assimilation System (ILDAS): A Coupled Hydrologic-Hydraulic System for Water Balance Assessments

Effective management of water resources requires reliable estimates of land surface states and fluxes, including water balance components. But most land surface models run in uncoupled mode and do not produce river discharge at catchment scales to be useful for water resources management applications. Such integrated systems are also rare over India where hydrometeorological extremes have wreaked havoc on the economy and people. So, an Indian Land Data Assimilation System (ILDAS) with a coupled land surface and a hydrodynamic model has been developed and driven by multiple meteorological forcings (0.1°, daily) to estimate land surface states, channel discharge, and floodplain inundation. ILDAS benefits from an integrated framework as well as the largest suite of observation records collected over India and has been used to produce a reanalysis product for 1981–2021 using four forcing datasets, namely, Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), ECMWF’s ERA-5, and Indian Meteorological Department (IMD) gridded precipitation. We assessed the uncertainty and bias in these precipitation datasets and validated all major components of the terrestrial water balance, i.e., surface runoff, soil moisture, terrestrial water storage anomalies, evapotranspiration, and streamflow, against a combination of satellite and in situ observation datasets. Our assessment shows that ILDAS can represent the hydrological processes reasonably well over the Indian landmass with IMD precipitation showing the best relative performance. Evaluation against ESA-CCI soil moisture shows that MERRA-2 based estimates outperform the others, whereas ERA-5 performs best in simulating evapotranspiration when evaluated against MODIS ET. Evaluations against observed records show that CHIRPS-based estimates have the highest performance in reconstructing surface runoff and streamflow. Once operational, this system will be useful for supporting transboundary water management decision making in the region.

Indian Land Data Assimilation System (ILDAS)↗
DOE OSTI2024

Interactions Between Climate and Species Drive Future Forest Carbon and Water Balances

Global change is altering forest carbon and water balances; however, the extent to which tree species shape ecosystem‐scale responses to climate, particularly in biodiverse forests, remains unclear. To address this, we simulated the effects of an envelope of future climate conditions on watershed carbon and water balances and quantified the contributions of tree species based on their xylem anatomy. We accomplished this by incorporating species‐level transpiration calculations into a landscape‐scale ecosystem process model. Our revised model linked the effects of forest succession, species composition, and climate change on water and carbon. Calibration of forest water fluxes using sap flux measurements and catchment water balances captured variability in species transpiration and interannual ET in biodiverse, humid temperate forest catchments in the southern Blue Ridge Mountains, USA. Across wet and dry future climate projections, ET increased, and streamflow and net carbon uptake decreased, particularly under a scenario of increasing drought. Despite accounting for just 30% of current biomass, diffuse‐porous tree species were the main driver of carbon and water flux responses now and in the future, thus intensifying the increase in ET and decline in streamflow. As diffuse‐porous biomass continues to increase, these forests will be increasingly sensitive to drought, amplifying losses of carbon sequestration and freshwater delivery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Improving the prediction of daily reservoir releases over the CONUS using conditioned LSTM

Reservoirs play a vital role in regulating streamflow timing and variability for hydroelectricity, flood control, water supply, irrigation, and recreation. Despite their importance, many reservoirs lack comprehensive operational guidelines, making their management complex due to conflicting operational objectives. Hence traditional policy-based reservoir models often fail to capture real-world conditions accurately and they depend on perfect streamflow predictions, which are not always available. In contrast, data-driven models like Long Short-Term Memory (LSTM) networks offer a robust alternative. This study introduces an approach that integrates reservoir characteristics—such as main use, climate, and maximum capacity—into the LSTM model to enhance reservoir release predictions. Using data from nearly 200 reservoirs in the contiguous United States (CONUS), our conditioned LSTM model (LSTM_cond) was compared with both the vanila LSTM and a traditional policy-based approach. Furthermore, our results show that while both LSTM_cond and LSTM perfoms better than the policy-based approach, LSTM_cond consistently outperforms LSTM for hydroelectric, water supply, irrigation, and recreation reservoirs. The KGE median values for LSTM_cond for out-sample reservoirs are 0.764, 0.565, 0.821, and 0.779, respectively, for the aforementioned reservoir types, which are consistently higher that the corresponding KGE values of 0.737, 0.413, 0.775, and 0.713 of LSTM, demonstrating its advantages in improving generalizability.

CONUS↗
DOE OSTI2025

Advancing stream temperature prediction with a generalizable large-sample framework across CONUS river reaches

Accurately predicting stream temperature in ungauged basins remains a critical challenge for water resource management, thermoelectric power plant cooling, and ecosystem conservation. Large-sample machine learning models trained on hundreds of well-monitored river basins have shown remarkable performance; however, such models have yet to be developed solely using forcing data that can be readily extracted to simulate stream temperatures anywhere in the contiguous United States (CONUS). In this study, we present a scalable, large-sample deep learning framework using Long Short-Term Memory (LSTM) networks to simulate daily stream temperatures in ungauged basins across the CONUS. The framework leverages both modeled reanalysis of meteorological and streamflow inputs as well as static attributes available for all 2.7 million CONUS river reaches in the National Hydrography Dataset Plus (NHDPlusV2). By generating dynamical inputs from predefined thermally relevant upstream contributing areas, rather than the entire upstream basin, the model also offers improvements in very large basins where full-basin averaging can dilute the most important influences on stream temperature. Evaluated across 300 basins, the model achieves a median Mean Absolute Error (MAE) of 1.1 °C and a Nash-Sutcliffe Efficiency (NSE) of 0.95 on temporally and spatially distinct test folds—comparable to models trained exclusively using meteorological and streamflow observational data. The flexible, high-performing framework generalizes to any unmonitored river reach without significant regulation or unnatural thermal input immediately upstream, substantially expanding predictive capabilities in data-scarce regions.

Hydrology↗
DOE OSTI2026

Systematic Evaluation of Atmospheric Forcing, Surface Datasets, and Mesh Effects on Kilometer-Scale Land Surface and River Modeling

Earth system models are advancing toward kilometer-scale resolution to capture local climate impacts and extremes. High-resolution land and river modeling depends on multiple factors, including mesh, surface datasets, and atmospheric forcing, but their relative effects at kilometer scales remain unquantified. We evaluated five Energy Exascale Earth System Model land and river configurations over the Mid-Atlantic region using two mesh (1/8° structured versus variable-resolution unstructured mesh), two surface datasets (default versus newly developed), and three atmospheric forcings (NLDAS2, MSWX, GSWP). Evaluation against satellite, reanalysis, and in situ benchmarks across water, energy, and carbon cycles quantifies how these factors affect model performance. Forcing selection produces the largest bias reductions (12-99% across variables), followed by surface datasets (7-75%) and mesh (up to 21%). Forcing effects vary by variable, with MSWX reducing biases for snow water equivalent, evapotranspiration, albedo, temperature, and gross primary productivity, GSWP for snow cover and runoff, and NLDAS for soil moisture and streamflow. The use of newly developed surface datasets improves gross primary productivity (58% bias reduction) and evapotranspiration but increase soil moisture and albedo biases due to current modeling limitations. Variable-resolution unstructured mesh improves the simulation of small-basin streamflow through better capturing drainage networks, though mesh minimally affects other land variables. These findings provide important guidance for high-resolution modeling development and actionable science.

Land and River modeling↗
DOE OSTI2025

The Role of Snowmelt and Subsurface Heterogeneity in Headwater Hydrology of a Mountainous Catchment in Colorado: A Model‐Data Integration Approach

Mountainous headwater streams are sustained by both snowmelt‐driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model‐data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface‐subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed data sets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater levels. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process‐based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.

Wang, Lijing [University of Connecticut, Storrs, C↗
DOE OSTI2026

When Do Riverine Systems “Feel the Burn”? Simulating How Burn Extent and Severity Modulate Hydrologic Controls on Biogeochemical Export

Wildfires impact terrestrial landscapes and downstream river corridors through shifts in vegetation and soil properties leading to downstream hydrologic and water quality impacts. The magnitude of these impacts depend on a complex and interconnected set of wildfire, landscape, and aquatic processes. Here, we isolate the impact of post-fire hydrologic changes on streamflow, nitrate, and dissolved organic carbon using the Soil and Water Assessment Tool (SWAT) model. We explore how responses differ across burn severity and area burned in two test basins: a humid forested basin and a semi-arid mixed land use basin. We ran 1830 wildfire simulations testing impacts of area burned, burn severity, and post-fire precipitation on streamflow, nitrate, and dissolved organic carbon. Our work suggests that area burned thresholds differ with burn severity and analyte. Additionally, post-fire transport of dissolved organic carbon was sensitive to both area burned and severity, while nitrate was primarily sensitive to area burned. Despite a muted (−9.5 to 5.7 mm yr −1 change) hydrologic response in the semi-arid basin, the model predicted large (7%–288% increase) shifts in dissolved organic carbon, suggesting that post-fire shifts in flow pathways and soil properties are key in its response. The limited shifts in nitrate responses in the simulations highlight that terrestrial post-fire transformations, rather than hydrologic changes, may control the increases in stream nitrate often observed post-fire. As wildfire regimes are shifting, improving understanding of post-fire nutrient export responses is critical to protect freshwater resources and aquatic ecosystems.

Wampler, Katherine A. [Pacific Northwest National ↗