Preliminary Test of ERTS-1 Imagery for Improving Definition of Natural Streamflow
Preliminary test of ERTS-1 imagery for improving definition of natural streamflow, using two basins in Chesapeake Bay area
Hollyday, E. F.↗
SEARCH · Search NASA
Results for “naturalized 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.
Quantifying Dam-Induced Fluctuations in Streamflow Frequencies Across the Colorado River Basin
Periodic fluctuations in natural streamflow are a major driver of river ecosystem dynamics and water resource management. However, most U.S. rivers are impacted both by long-term hydroclimatic trends and dams that alter flow variability. Despite these impacts, it remains largely unexplored how dams affect the dominant frequencies of natural streamflow over a highly regulated river network. We investigated the entire Colorado River Basin (CRB) to understand how the annual (10–14 months) and multi-annual (24–60 months) frequencies in natural flow regimes have been progressively altered by dams. Given the significant alteration over the CRB, we captured changes in streamflow frequencies between naturalized and observed monthly flows via wavelet analysis. Based on the similarity of changes in streamflow frequencies (annual and multi-annual) over the last 30 years, sections of the riverine network were classified into four groups. The annual frequency was relatively well preserved downstream of Hoover Dam, while showing a systematic trend of alteration downstream of Glen Canyon Dam until Hoover Dam. Meanwhile, the multi-annual frequency component was highly altered for the entire Lower Colorado main stem (i.e., downstream of Glen Canyon). We also identified dams with significant impacts on streamflow frequency by comparing wavelet coherence estimates. This study advances the notion that dams fundamentally alter river flow regimes across multiple frequencies and with varying amplitudes over time and space, with alteration propagating – or being dampened – by both hydroclimatic fluctuations and water resource management.
54 ENVIRONMENTAL SCIENCES↗
Estimating Future Surface Water Availability Through an Integrated Climate‐Hydrology‐Management Modeling Framework at a Basin Scale Under CMIP6 Scenarios
Abstract Climate change and increasing water demand due to population growth pose serious threats to surface water availability. The biggest challenge in addressing these threats is the gap between climate science and water management practices. Local water planning often lacks the integration of climate change information, especially with regard to its impacts on surface water storage and evaporation as well as the associated uncertainties. Using Texas as an example, state and regional water planning relies on the use of reservoir “Firm Yield” (FY)—an important metric that quantifies surface water availability. However, this existing planning methodology does not account for the impacts of climate change on future inflows and on reservoir evaporation. To bridge this knowledge gap, an integrated climate‐hydrology‐management (CHM) modeling framework was developed, which is generally applicable to river basins with geographical, hydrological, and water right settings similar to those in Texas. The framework leverages the advantages of two modeling approaches—the Distributed Hydrology Soil Vegetation Model (DHSVM) and Water Availability Modeling (WAM). Additionally, the Double Bias Correction Constructed Analogues method is utilized to downscale and incorporate Coupled Model Intercomparison Project Phase 6 GCMs. Finally, the DHSVM simulated naturalized streamflow and reservoir evaporation rate are input to WAM to simulate reservoir FY. A new term—“Ratio of Firm Yield” (RFY)—is created to compare how much FY changes under different climate scenarios. The results indicate that climate change has a significant impact on surface water availability by increasing reservoir evaporation, altering the seasonal pattern of naturalized streamflow, and reducing FY.
54 ENVIRONMENTAL SCIENCES↗
Insights From Dayflow: A Historical Streamflow Reanalysis Dataset for the Conterminous United States
Abstract Reconstructed historical streamflow time series can supplement limited streamflow gauge observations. However, there are common challenges of typical modeling approaches: process‐based hydrologic models can be data/computation‐intensive, and statistics‐based models can be region/stream‐specific. Here we present a nationally scalable modeling framework integrating the simulated runoff from the Variable Infiltration Capacity (VIC) model with the Routing Application for Parallel computatIon of Discharge (RAPID) routing model leveraging high‐performance computing. We demonstrate an efficient method of assimilating streamflow at US Geological Survey (USGS) streamflow monitoring sites using a simple hierarchical approach in the VIC‐RAPID framework. The result is a reconstructed 36‐year (1980–2015) daily and monthly streamflow dataset (Dayflow) at ∼2.7 million NHDPlusV2 stream reaches in the conterminous US (CONUS). We perform a comprehensive evaluation at 7,526 USGS sites and characterize their error statistics. The results demonstrate that 49% of the USGS sites demonstrate Kling–Gupta Efficiency (KGE) > 0.5 and 58% of the sites show percentage bias within ±20% for the daily naturalized streamflow. Streamflow data assimilation across CONUS shows an overall improvement over naturalized streamflow, notably in the western semiarid‐to‐arid regions. Comparison to other national and global streamflow reanalysis datasets such as the National Water Model and Global Reach‐scale A priori Discharge Estimates for SWOT demonstrates improved KGE, reduced bias, and directions for Dayflow improvements. Investigations of error statistics with key hydrologic, hydroclimatic, and geomorphologic basin characteristics reveal region‐specific patterns which may help improve future framework applications. Overall, Dayflow may enable a better understanding of hydrologic conditions in a changing environment, especially in locations currently not represented by streamflow monitoring networks.
54 ENVIRONMENTAL SCIENCES↗
Achieving Realistic Runoff in the Western United States with a Land Surface Model Forced by Dynamically Downscaled Meteorology
Abstract In this study, we calibrate a regional climate model’s (RCM) underlying land surface model (LSM). In addition to providing a realistic representation of runoff across the hydroclimatically diverse western United States, this is done to take advantage of the RCM’s ability to physically resolve meteorological forcing data in ungauged regions, and to prepare the calibrated hydrologic model for tight coupling, or the ability to represent land surface–atmosphere interactions, with the RCM. Specifically, we use a 9-km resolution meteorological forcing dataset across the western United States, from the fifth generation ECMWF Reanalysis (ERA5) downscaled by the Weather Research Forecasting (WRF) regional climate model, as an offline forcing for Noah-Multiparameterization (Noah-MP). We detail the steps involved in producing an LSM capable of accurately representing runoff, including physical parameterization selection, parameter calibration, and regionalization to ungauged basins. Based on our model evaluation from 1954 to 2021 for 586 basins with daily natural streamflow, the streamflow bias is reduced from 24.2% to 4.4%, and the median daily Nash–Sutcliffe efficiency (NSE) is improved from 0.12 to 0.36. When validating against basins with monthly natural streamflow data, we obtain a similar reduction in bias and a median monthly NSE improvement from 0.18 to 0.56. In this study, we also discover the optimal setup when using a donor-basin method to regionalize parameters to ungauged basins, which can vary by 0.06 NSE for unique designs of this regionalization method. Significance Statement This study provides useful guidance for improving a land surface model to accurately represent runoff across a spatially extensive and hydroclimatically diverse region (the western United States). The land surface model is updated to represent runoff more accurately at gauged basins, and then additionally updated for basins without observational data using a mathematical approach called the donor-basin method. We make use of a regional climate model’s reanalysis-derived meteorological data and its underlying land surface model to achieve realistic runoff. The calibrated land surface model can thus be tightly coupled in subsequent studies in a manner that should more accurately reflect runoff conditions. Findings from this study will serve as a useful reference for the atmospheric (and hydrologic) modeling communities and their ability to represent large-scale hydrology accurately.
Meteorology & Atmospheric Sciences↗
Dayflow-PR: High-Resolution Streamflow Reanalysis for Puerto Rico, Version 1.0
This dataset presents a high-resolution historical streamflow reanalysis for NHDPlusV2 stream reaches across Puerto Rico (PR) spanning 1950 - 2019. The reanalysis is generated using the calibrated VIC-RAPID hydrologic modeling framework at the Hydrologic Unit Code Sub-basin (HUC08) scale, forced with sub-daily and daily meteorological forcings from Daymet. Runoff is simulated on 1- and 6-km grids, and the resulting total runoff is routed through the NHDPlusV2 river network using the RAPID routing model to produce Naturalized Streamflow Reanalysis. Where complete observational records are available over 1980 - 2019, streamflows are assimilated (substituted) and subsequently routed downstream through the river network to produce Assimilated Streamflow Reanalysis. The dataset includes streamflow outputs from eight distinct hydrologic modeling configurations along with key performanc evaluation metrics at daily and monthly scales, supporting a wide range of water resource applications. This dataset is derived to support the Non-Powered Dam Assessment, as well as 9505 Secure Water Assessment projects for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details, refer to Ghimire et al. (2023), Kao et al. (2024), and Ghimire et al. (2025).
13 HYDRO ENERGY↗
Prediction of Hydrological Drought: What Can We Learn From Continental-Scale Offline Simulations?
Land surface model experiments are used to quantify, across the coterminous United States, the contributions (isolated and combined) of soil moisture and snowpack initialization to the skill of seasonal streamflow forecasts at multiple leads and for different start dates. Forecasted streamflows are compared to naturalized streamflow observations where available and to synthetic (model-generated) streamflow data elsewhere. We find that snow initialization has a major impact on skill in the mountainous western U.S. and in a portion of the northern Great Plains; a mid-winter (January 1) initialization of snow in these areas leads to significant skill in the spring melting season. Soil moisture initialization also contributes to skill, and although the maximum contributions are not as large as those seen for snow initialization, the soil moisture contributions extend across a much broader geographical area. Soil moisture initialization can contribute to skill at long leads (up to 5 or 6 months), particularly for forecasts issued during winter.
Koster, Randal↗
Dayflow: CONUS Daily Streamflow Reanalysis, Version 2 (DayflowV2)
The DayflowV2 dataset provides multiple meteorologic forcings driven hourly streamflow information for approximately 2.7 million NHDPlusV2 stream reaches in the conterminous US (CONUS). DaymetV4, Stage-IV, and Analysis of Period of Record for Calibration (AORC) forcings and their corresponding hybrids drive a nationally scalable modeling framework integrating the simulated runoff from the Variable Infiltration Capacity (VIC) model with the Routing Application for Parallel computatIon of Discharge (RAPID) routing model. Streamflow with (Assimilated) and without (Naturalized) streamflow assimilation at US Geological Survey (USGS) streamflow monitoring sites are included in DayflowV2. A comprehensive evaluation of streamflow at 7,526 USGS gauges is performed for both streamflow types. The resulting key evaluation metrics are also included in the Dayflow dataset. The reanalysis data are available for variable periods; 36 years (1980-2015) for DaymetV4 (DayflowV1), 18 years (2002-2019) for Stage-IV and its hybrids, and 40 years (1980-2019) for AORC and its hybrids.
13 HYDRO ENERGY↗
Assessing the complex influences of water management on hydrological drought characteristics in Texas
The state of Texas in the United States is highly susceptible to drought. Its major rivers are subject to extensive water management (WM) activities in order to sustain multisectoral water demands, particularly during drought conditions. However, the impact of WM on the propagation dynamics and characteristics of hydrological drought (HD) in Texas remains unclear. To fill this gap, this study quantifies the influence of WM across 32 streamflow gauges along the mainstems of seven major rivers in Texas by comparing a variety of drought metrics under natural and managed conditions. Notably, we leveraged an extensive, naturalized streamflow dataset constructed by the Texas Commission on Environmental Quality, paired with gauge observations of managed conditions. Results indicate that at the multi-decadal scale, WM significantly reduced HD frequency across all seven rivers and at 81% of the gauges analyzed. Additionally, it increased the response timescale of HD across Texas’ major rivers by a median of 2.5 months. Conversely, the average-event duration and severity increased in most locations. Temporal analysis reveals that the WM impact on HD varied seasonally, with attenuation effects during mid-summer and early fall and intensification effects during late winter and spring. Additionally, WM was found to greatly increase the spatial variability of HD characteristics across the region. These findings emphasize the complexity of WM effects on HD and the necessity for nuanced strategies in managing HD under WM influences.
54 ENVIRONMENTAL SCIENCES↗
The Contribution of Soil Moisture Information to Forecast Skill: Two Studies
This talk briefly describes two recent studies on the impact of soil moisture information on hydrological and meteorological prediction. While the studies utilize soil moisture derived from the integration of large-scale land surface models with observations-based meteorological data, the results directly illustrate the potential usefulness of satellite-derived soil moisture information (e.g., from SMOS and SMAP) for applications in prediction. The first study, the GEWEX- and ClIVAR-sponsored GLACE-2 project, quantifies the contribution of realistic soil moisture initialization to skill in subseasonal forecasts of precipitation and air temperature (out to two months). The multi-model study shows that soil moisture information does indeed contribute skill to the forecasts, particularly for air temperature, and particularly when the initial local soil moisture anomaly is large. Furthermore, the skill contributions tend to be larger where the soil moisture initialization is more accurate, as measured by the density of the observational network contributing to the initialization. The second study focuses on streamflow prediction. The relative contributions of snow and soil moisture initialization to skill in streamflow prediction at seasonal lead, in the absence of knowledge of meteorological anomalies during the forecast period, were quantified with several land surface models using uniquely designed numerical experiments and naturalized streamflow data covering mUltiple decades over the western United States. In several basins, accurate soil moisture initialization is found to contribute significant levels of predictive skill. Depending on the date of forecast issue, the contributions can be significant out to leads of six months. Both studies suggest that improvements in soil moisture initialization would lead to increases in predictive skill. The relevance of SMOS and SMAP satellite-based soil moisture information to prediction are discussed in the context of these studies.
Koster, Randal↗
Physical Controls on Irrigation Return Flow Contributions to Stream Flow in Irrigated Alluvial Valleys
Irrigation can be a significant source of groundwater recharge in many agricultural regions, particularly in arid and semi-arid climates. Once infiltrated, irrigation recharge can travel via subsurface flowpaths that return to the river system in a lagged manner, supplementing natural streamflow weeks, months, or even years from when the irrigation was applied. In regions that experience low flows during summer and early fall, return flows can be a significant source of supplementary streamflow. Many water planning and operations models either ignore return flows or roughly approximate them with analytical solutions. Thus, return flows represent an important but often overlooked component of the hydrological exchange and overall water balance in agricultural regions. This study uses groundwater models to explore a wide range of factors that control irrigation return flow timing in irrigated alluvial valleys. A sensitivity analysis approach is used to assess how factors such as the extent of irrigated land adjacent to a stream, irrigation recharge rate, aquifer hydraulic conductivity, aquifer thickness, water table configuration, and seasonal fluctuations in stream stage control the timing of subsurface return flows. Modeling is conducted using MODFLOW models representing an irrigated alluvial valley adjacent to a stream. While a simplification of the full complexity in real systems, the models are a significant advancement from the analytical solution and provide new insight into the timescales of return flows over a broad range of possible conditions. To contextualize our modeling results, they are compared to an analytical solution commonly used for approximating return flows to evaluate its performance. Our findings show what factors and conditions influence return flow timing and control whether they contribute to streamflow over short term (months) or longer term (seasonal) time scales.
Ferencz, Stephen B.↗
Nature-based solutions in mountain catchments reduce impact of anthropogenic climate change on drought streamflow
Quantifying how well Nature-based Solutions can offset anthropogenic climate change impacts is important for adaptation planning, but has rarely been done. Here we show that a widely-applied Nature-based Solution in South Africa – invasive alien tree clearing – reduces the impact of anthropogenic climate change on drought streamflow. Using a multi-model joint-attribution of climate and landscape-vegetation states during the 2015–2017 Cape Town “Day Zero” drought, we find that anthropogenic climate change reduced streamflow by 12–29% relative to a counterfactual world with anthropogenic emissions removed. This impact on streamflow was larger than corresponding reductions in rainfall (7–15%) and reference evapotranspiration (1.7–2%). Clearing invasive alien trees could have ameliorated streamflow reductions by 3–16% points for moderate invasions levels. Preventing further invasive alien tree spread avoided potential additional reductions of 10–27% points. Total clearing could not have offset the anthropogenic climate change impact completely. Invasive alien tree clearing is an important form of catchment restoration for managing changing hydroclimatic risk, but will need to be combined with other adaptation options as climate change accelerates.
54 ENVIRONMENTAL SCIENCES↗
Natural Variability Has Concealed Increases in Western US Flood Hazard Since the 1970s
Flood hazard across the western United States (US) has generally shown decreasing trends in recent decades. This region's extreme streamflow is highly influenced by natural variability, which could either mask or amplify anthropogenic streamflow trends. Here, in this study, we utilize a technique known as dynamical adjustment to assess historical (1970–2020) annual maximum 1-day streamflow (Qx1d) from unregulated basins across the western US with and without the impact of natural variability. After removing natural variability, the fraction of basins with a positive (>5%) trend in Qx1d shifts from 25% to 53%. Basins with increasing (decreasing) Qx1d trends after dynamical adjustment exhibit weak (strong) drying, and furthermore are associated with intensifying precipitation extremes and/or large decreases in snowpack. Increasing flood hazard will likely emerge for such basins as the current phase of natural decadal variability shifts, and anthropogenic signals continue to intensify.
54 ENVIRONMENTAL SCIENCES↗
Predicting Hydrological Drought: Relative Contributions of Soil Moisture and Snow Information to Seasonal Streamflow Prediction Skill
in this study we examine how knowledge of mid-winter snow accumulation and soil moisture conditions contribute to our ability to predict streamflow months in advance. A first "synthetic truth" analysis focuses on a series of numerical experiments with multiple sophisticated land surface models driven with a dataset of observations-based meteorological forcing spanning multiple decades and covering the continental United States. Snowpack information by itself obviously contributes to the skill attained in streamflow prediction, particularly in the mountainous west. The isolated contribution of soil moisture information, however, is found to be large and significant in many areas, particularly in the west but also in region surrounding the Great Lakes. The results are supported by a supplemental, observations-based analysis using (naturalized) March-July streamflow measurements covering much of the western U.S. Additional forecast experiments using start dates that span the year indicate a strong seasonality in the skill contributions; soil moisture information, for example, contributes to kill at much longer leads for forecasts issued in winter than for those issued in summer.
Koster, R.↗
Hydropower potential derived from streamflow extremes for Alaska, USA
Alaska is an expansive region known for its abundant natural resources, including thousands of miles of streams and rivers. These rivers represent potential opportunities for future hydropower development that could provide reliable energy supply for local communities. There is limited long-term high temporal resolution streamflow data available for the region, making data-driven estimates of potential hydropower and its variability across the state challenging. This study provides a novel data-driven approach for hydropower capacity estimation across Alaska. We use supervised machine learning to develop a relationship between the daily and peak flow duration curves in order to augment the size of our dataset from 44 sites to 67 sites. We perform a stochastic hydropower estimation across the 67 sites and identify approximately 1000 MW of total potential hydropower capacity distributed across these sites. Our study provides the first step towards more comprehensive hydropower estimation for this critical region, highlighting the need for future work integrating high-resolution spatial data, community needs, and economic constraints in estimates of potential hydropower development in Alaska.
Hydropower↗
Streamflow simulation in data-scarce basins using Bayesian and physics-informed machine learning models
Hydrologic predictions at rural watersheds are important but also challenging due to data shortage. Long short-term memory (LSTM) networks are a promising machine learning approach and have demonstrated good performance in streamflow predictions. However, due to its data-hungry nature, most LSTM applications focus on well-monitored catchments with abundant and high-quality observations. In this work, we investigate predictive capabilities of LSTM in poorly monitored watersheds with short observation records. To address three main challenges of LSTM applications in data-scarce locations, i.e., overfitting, uncertainty quantification (UQ), and out-of-distribution prediction, we evaluate different regularization techniques to prevent overfitting, apply a Bayesian LSTM for UQ, and introduce a physics-informed hybrid LSTM to enhance out-of-distribution prediction. Through case studies in two diverse sets of catchments with and without snow influence, we demonstrate that 1) when hydrologic variability in the prediction period is similar to the calibration period, LSTM models can reasonably predict daily streamflow with Nash–Sutcliffe efficiency above 0.8, even with only 2 years of calibration data; 2) when the hydrologic variability in the prediction and calibration periods is dramatically different, LSTM alone does not predict well, but the hybrid model can improve the out-of-distribution prediction with acceptable generalization accuracy; 3) L2 norm penalty and dropout can mitigate overfitting, and Bayesian and hybrid LSTM have no overfitting; and 4) Bayesian LSTM provides useful uncertainty information to improve prediction understanding and credibility. In conclusion, these insights have vital implications for streamflow simulation in watersheds where data quality and availability are a critical issue.
54 ENVIRONMENTAL SCIENCES↗
Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"
The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.
54 ENVIRONMENTAL SCIENCES↗
Data From: "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater"
This repository contains the data and code associated with the paper titled "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater," published in Nature Geoscience, 2026. This study seeks to answer how various ages of groundwater interact with mountainous streamflow in mountainous headwaters such as the East River. It includes various model-data processing scripts, primarily for ParFlow-CLM analysis of simulated water years 2015-2021, and two numerical warming experiments (+2.5 and +4.0 degrees C), including run scripts, forcing scripts, and post-processing, as well as comparison to observation datasets, detailed below. This data requires the use of R (.r, .rmd), Python (.py), Jupyter Notebook or Jupyter Lab (.ipynb), ParFLOW-CLM, EcoSLIM. Further information on the use of all file formats mentioned below (e.g. .tff. .nc) are provided within the associated scripts and directory where the files are located. Contents & Usage ASO/: Contains the bash and python scripts used to convert airborne snow observatory (ASO) data (ASO, 2023) in various data formats (georeferenced tiff file, NetCDF, UTM, and to latitude/longitude) then regrided to the ParFlow equivalent grid. Output data are in regrid_regll_data.zip and subsequently visualized and analyzed in plot_and_compare.py for Supplementary Figures A14 and A15. The wksht_ASO_comparison.xlsx spreadsheet is used to calculate the data for Supplementary Figure A16. EcoSLIM/: Contains the scripts and input files to run the EcoSLIM particle tracking simulations (/run_scripts) and the post-processing python script (/plot_scripts/eco_agedist_plots.ipynb). Jasechko et al./: Contains the jupyter notebook (Extract_Elevation.ipynb) to determine the outlet elevations of the 260 watersheds used in Jasechko et al. (2016), and the corresponding table, Table_S1_Watersheds_alt.csv. Used to create Supplementary Information Figure A2. PLM_Wells/: Contains the QA/QC-ed groundwater level time series of the PLM-1 and PLM-6 Monitoring Wells from Faybishenko et al. (2023), reformatted to water years used for Supplementary Figures A19 and and A20. ParFlow/: Contains the input files and run scripts to run ParFlow-CLM (/run_scripts), the python and tool command language (Tcl) scripts to create and distribute the ParFlow forcing simulation files (/forcing), and various scripts and intermediary files to analyze the model outputs (/post_process). SQUIRE/: Contains the processing scripts and intermediary files for the Surface QUantitatIve pRecipitation Estimation (SQUIRE) data (Grover, 2023) used to generate Supplementary Figure A18. USGS_Streamflow/: Contains the raw and gap-filled United States Geological Survey streamflow data (U.S. Geological Survey, 2026) used at the Almont station (site number 09112500). Gap-filling is performed in the R script with data from the Taylor station (site number 09110000). (/USGS_09112500_EAST_RIVER_AT_ALMONT_GAP_FILLED/code_almont_streamflow_gap_fill.Rmd). discharge/: Contains the gap-filled discharge data at the Watershed Function SFA East River pumphouse site (Newcomer et al., 2022) used to generate Supplementary Figure A13 and to compute hourly Nash-Sutcliffe model efficiency coefficients (NSE) in Table A4. snotel_and_flux_tower/: Contains the snow telemetry data (U.S. Department of Agriculture, 2024) from the Butte (site ID 380) and Schofield (site ID 737) stations, reformatted by water year, accessed with the snotelr R package. Used to create Supplementary Figure A17. Also contains the flux tower observational data (FluxTower_Pumphouse_ESS-DIVE.ET_only.h.txt) from Ryken et al. (2022) and sap flux transpiration data (MaxB_Transpiration_5Sites.daily_sums.h.txt) from Ryken (2021), used to create Supplementary Figures A22 and A23, respectively. Raw EcoSLIM model outputs are in excess of 24TB, and are stored on National Energy Research Scientific Computing Center (NERSC) and publicly available via the external link provided in the paper.
atmospheric warming↗