Search NASA⌕ Search

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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A network approach for multiscale catchment classification using traits

Abstract. The classification of river catchments into groups with similar biophysical characteristics is useful to understand and predict their hydrological behavior. The increasing availability of remote sensing and other large-scale geospatial datasets has enabled the use of advanced data-driven approaches to classify catchments using traits such as topography, geology, climate, land cover, land use, and human influence. Unsupervised clustering algorithms based on the Euclidean distance are commonly used for trait-based classification but are not suitable for highly dimensional data. In this study we present a new network-based method for multi-scale catchment classification, which can be applied to large datasets and used to determine the traits associated with different catchment groups. In this framework, two networks are analyzed in parallel: the first being where the nodes are traits and the second being where the nodes are catchments. In both cases, edges represent pairwise similarity, and a network cluster detection algorithm is used for the classification. The trait network is used to investigate redundancy in the trait data and to condense this information into a small number of interpretable categories. The catchments network is used to classify the catchments into clusters and to identify representative catchments for the different groups using the degree centrality metric. We apply this method to classify 9067 river catchments across the contiguous United States at both regional and continental scales using 274 non-categorical traits. At the continental scale, we identify 25 interpretable trait categories and 34 catchment clusters of sizes greater than 50. We find that catchments with similar trait categories are typically located in the same region, with different spatial patterns emerging among clusters dominated by natural and anthropogenic traits. We also find that the catchment clusters exhibit distinct hydrological behavior based on an analysis of streamflow indices. This network approach provides several advantages over traditional means of classification, including better separation of clusters, the use of alternate similarity metrics that are more suitable for highly dimensional data, and reducing redundancy in the trait information. The paired catchment–trait networks enable analysis of hydrological behavior using the dominant trait categories for each catchment cluster. The approach can be used at multiple spatial scales since the network topologies adjust automatically to reflect the trait patterns at the scale of investigation. Finally, the representative catchments identified as hub nodes in the network can be used to guide transferable observational and modeling strategies. The method is broadly applicable beyond hydrology for classification of other complex systems that utilize different types of trait datasets.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Using Historical Precipitation, Temperature, and Runoff Observations to Evaluate Evaporation Formulations in Land Surface Models

Key to translating soil moisture memory into subseasonal precipitation and air temperature forecast skill is a realistic treatment of evaporation in the forecast system used - in particular, a realistic treatment of how evaporation responds to variations in soil moisture. The inherent soil moisture-evaporation relationships used in today's land surface models (LSMs), however, arguably reflect little more than guesswork given the lack of evaporation and soil moisture data at the spatial scales represented by regional and global models. Here we present a new approach for evaluating this critical aspect of LSMs. Seasonally averaged precipitation is used as a proxy for seasonally-averaged soil moisture, and seasonally-averaged air temperature is used as a proxy for seasonally-averaged evaporation (e.g., more evaporative cooling leads to cooler temperatures)~ the relationship between historical precipitation and temperature measurements accordingly mimics in certain important ways nature's relationship between soil moisture and evaporation. Additional information on the relationship is gleaned from joint analysis of precipitation and streamflow measurements. An experimental framework that utilizes these ideas to guide the development of an improved soil moisture-evaporation relationship is described and demonstrated.

Koster, Randal D.↗
DOE OSTI2022

Technology Pathways Could Help Drive the U.S. West Coast Grid's Exposure to Hydrometeorological Uncertainty

Abstract Previous studies investigating deep decarbonization of bulk electric power systems and wholesale electricity markets have not sufficiently explored how future grid pathways could affect the grid's vulnerability to hydrometeorological uncertainty on multiple timescales. Here, we employ a grid operations model and a large synthetic weather ensemble to “stress test” a range of future grid pathways for the U.S. West Coast developed by ReEDS, a well‐known capacity planning model. Our results show that gradual changes in the underlying capacity mix from 2020 to 2050 can cause significant “re‐ranking” of weather years in terms of annual wholesale electricity prices (with “good” years becoming bad, and vice versa). Nonetheless, we find the highest and lowest ranking price years in terms of average electricity price remain mostly tied to extremes in hydropower availability (streamflow) and load (summer temperatures), with the strongest sensitivities related to drought. Seasonal dynamics seen today involving spring snowmelt and hot, dry summers remain well‐defined out to 2050. In California, future supply shortfalls in our model are concentrated in the evening and occur mostly during periods of high temperature anomalies in late summer months and in late winter; in the Pacific Northwest, supply shortfalls are much more strongly tied to negative streamflow anomalies. Under our more robust sampling of stationary hydrometeorological uncertainty, we also find that the ratio of dis‐patchable thermal (i.e., natural gas) capacity to wind and solar required to ensure grid reliability can differ significantly from values reported by ReEDS.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Evaluating 3 decades of precipitation in the Upper Colorado River basin from a high-resolution regional climate model

Abstract. Convection-permitting regional climate models (RCMs) have recently become tractable for applications at multi-decadal timescales. These types of models have tremendous utility for water resource studies, but better characterization of precipitation biases is needed, particularly for water-resource-critical mountain regions, where precipitation is highly variable in space, observations are sparse, and the societal water need is great. This study examines 34 years (1987–2020) of RCM precipitation from the Weather Research and Forecasting model (WRF; v3.8.1), using the Climate Forecast System Reanalysis (CFS; CFSv2) initial and lateral boundary conditions and a 1 km × 1 km innermost grid spacing. The RCM is centered over the Upper Colorado River basin, with a focus on the high-elevation, 750 km2 East River watershed (ERW), where a variety of high-impact scientific activities are currently ongoing. Precipitation is compared against point observations (Natural Resources Conservation Service Snow Telemetry or SNOTEL), gridded climate datasets (Newman, Livneh, and PRISM), and Bayesian reconstructions of watershed mean precipitation conditioned on streamflow and high-resolution snow remote-sensing products. We find that the cool-season precipitation percent error between WRF and 23 SNOTEL gauges has a low overall bias (x^ = 0.25 %, s = 13.63 %) and that WRF has a higher percent error during the warm season (x^ = 10.37 %, s = 12.79 %). Warm-season bias manifests as a high number of low-precipitation days, though the low-resolution or SNOTEL gauges limit some of the conclusions that can be drawn. Regional comparisons between WRF precipitation accumulation and three different gridded datasets show differences on the order of ± 20 %, particularly at the highest elevations and in keeping with findings from other studies. We find that WRF agrees slightly better with the Bayesian reconstruction of precipitation in the ERW compared to the gridded precipitation datasets, particularly when changing SNOTEL densities are taken into account. The conclusions are that the RCM reasonably captures orographic precipitation in this region and demonstrates that leveraging additional hydrologic information (streamflow and snow remote-sensing data) improves the ability to characterize biases in RCM precipitation fields. Error characteristics reported in this study are essential for leveraging the RCM model outputs for studies of past and future climates and water resource applications. The methods developed in this study can be applied to other watersheds and model configurations. Hourly 1 km × 1 km precipitation and other meteorological outputs from this dataset are publicly available and suitable for a wide variety of applications.

Geology↗
DOE OSTI2023

Pore-resolved investigation of turbulent open channel flow over a randomly packed permeable sediment bed

Pore-resolved direct numerical simulations are performed to investigate the interactions between streamflow turbulence and groundwater flow through a randomly packed porous sediment bed for three permeability Reynolds numbers, Re K = 2.56 , 5.17 and 8.94, representative of natural stream or river systems. Time–space averaging is used to quantify the Reynolds stress, form-induced stress, mean flow and shear penetration depths, and mixing length at the sediment–water interface (SWI). Here, the mean flow and shear penetration depths increase with Re K and are found to be nonlinear functions of non-dimensional permeability. The peaks and significant values of the Reynolds stresses, form-induced stresses, and pressure variations are shown to occur in the top layer of the bed, which is also confirmed by conducting simulations of just the top layer as roughness elements over an impermeable wall. The probability distribution functions (p.d.f.s) of normalized local bed stress are found to collapse for all Reynolds numbers, and their root-mean-square fluctuations are assumed to follow logarithmic correlations. The fluctuations in local bed stress and resultant drag and lift forces on sediment grains are mainly a result of the top layer; their p.d.f.s are symmetric with heavy tails, and can be well represented by a non-Gaussian model fit. The bed stress statistics and the pressure data at the SWI potentially can be used in providing better boundary conditions in modelling of incipient motion and reach-scale transport in the hyporheic zone.

turbulence simulation↗
DOE OSTI2022

A hydrological simulation dataset of the Upper Colorado River Basin from 1983 to 2019

Abstract This article presents a hydrological reconstruction of the Upper Colorado River Basin with an hourly temporal resolution, and 1-km spatial resolution from October 1982 to September 2019. The validated dataset includes a suite of hydrologic variables including streamflow, water table depth, snow water equivalent (SWE) and evapotranspiration (ET) simulated by an integrated hydrological model, ParFlow-CLM. The dataset was validated over the period with a combination of point observations and remotely sensed products. These datasets provide a long-term, natural-flow, simulation for one of the most over-allocated basins in the world.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2014

"Efficiency Space" - A Framework for Evaluating Joint Evaporation and Runoff Behavior

At the land surface, higher soil moisture levels generally lead to both increased evaporation for a given amount of incoming radiation (increased evaporation efficiency) and increased runoff for a given amount of precipitation (increased runoff efficiency). Evaporation efficiency and runoff efficiency can thus be said to vary with each other, motivating the development of a unique hydroclimatic analysis framework. Using a simple water balance model fitted, in different experiments, with a wide variety of functional forms for evaporation and runoff efficiency, we transform net radiation and precipitation fields into fields of streamflow that can be directly evaluated against observations. The optimal combination of the functional forms the combination that produces the most skillful stream-flow simulations provides an indication for how evaporation and runoff efficiencies vary with each other in nature, a relationship that can be said to define the overall character of land surface hydrological processes, at least to first order. The inferred optimal relationship is represented herein as a curve in efficiency space and should be valuable for the evaluation and development of GCM-based land surface models, which by this measure are often found to be suboptimal.

GCM↗
NASA NTRS2005

The Hurricane-Flood-Landslide Continuum

In August 2004, representatives from NOAA, NASA, the USGS, and other government agencies convened in San Juan, Puerto Rim for a workshop to discuss a proposed research project called the Hurricane-Flood-Landslide Continuum (HFLC). The essence of the HFLC is to develop and integrate tools across disciplines to enable the issuance of regional guidance products for floods and landslides associated with major tropical rain systems, with sufficient lead time that local emergency managers can protect vulnerable populations and infrastructure. All three lead agencies are independently developing precipitation-flood-debris flow forecasting technologies, and all have a history of work on natural hazards both domestically and overseas. NOM has the capability to provide tracking and prediction of storm rainfall, trajectory and landfall and is developing flood probability and magnTtude capabilities. The USGS has the capability to evaluate the ambient stability of natural and man-made landforms, to assess landslide susceptibilities for those landforms, and to establish probabilities for initiation of landslides and debris flows. Additionally, the USGS has well-developed operational capacity for real-time monitoring and reporting of streamflow across distributed networks of automated gaging stations (http://water.usgs.gov/waterwatch/). NASA has the capability to provide sophisticated algorithms for satellite remote sensing of precipitation, land use, and in the future, soil moisture. The Workshop sought to initiate discussion among three agencies regarding their specific and highly complimentary capabilities. The fundamental goal of the Workshop was to establish a framework that will leverage the strengths of each agency. Once a prototype system is developed for example, in relatively data-rich Puerto Rim, it could be adapted for use in data-poor, low-infrastructure regions such as the Dominican Republic or Haiti. This paper provides an overview of the Workshop s goals, presentations and recommendations with respect to the development of the HFLC.

Negri, Andrew J.↗
DOE OSTI2022

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗
DOE OSTI2023

Data for Hadjimichael et al. - Multi-actor, multi-impact scenario discovery

Naming convention for each file is S[scenario_id]_[realization_id].parquet All realizations are generated as per the accompanying paper: Hadjimichael, A., Reed, P.M., Quinn, J.D, Vernon, C.R., Thurber, T., Multi-actor, multi-impact scenario discovery of consequential narrative storylines for human-natural systems planning. Earth's Future (In Review) Each .parquet file contains the following columns: 'structure_id' : StateMod structure ID data corresponds to 'year': simulation year 'month': simulation month 'demand': water demand in acre-feet 'shortage': water shortage in acre-feet 'river_outflow': streamflow in diversion location in acre-feet 'control_location': other structure or reason for shortage

Hadjimichael, Antonia↗
DOE OSTI2022

A graph neural network (GNN) approach to basin-scale river network learning: the role of physics-based connectivity and data fusion

Abstract. Rivers and river habitats around the world are under sustained pressure from human activities and the changing global environment. Our ability to quantify and manage the river states in a timely manner is critical for protecting the public safety and natural resources. In recent years, vector-based river network models have enabled modeling of large river basins at increasingly fine resolutions, but are computationally demanding. This work presents a multistage, physics-guided, graph neural network (GNN) approach for basin-scale river network learning and streamflow forecasting. During training, we train a GNN model to approximate outputs of a high-resolution vector-based river network model; we then fine-tune the pretrained GNN model with streamflow observations. We further apply a graph-based, data-fusion step to correct prediction biases. The GNN-based framework is first demonstrated over a snow-dominated watershed in the western United States. A series of experiments are performed to test different training and imputation strategies. Results show that the trained GNN model can effectively serve as a surrogate of the process-based model with high accuracy, with median Kling–Gupta efficiency (KGE) greater than 0.97. Application of the graph-based data fusion further reduces mismatch between the GNN model and observations, with as much as 50 % KGE improvement over some cross-validation gages. To improve scalability, a graph-coarsening procedure is introduced and is demonstrated over a much larger basin. Results show that graph coarsening achieves comparable prediction skills at only a fraction of training cost, thus providing important insights into the degree of physical realism needed for developing large-scale GNN-based river network models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2022

The Simulation and Subseasonal Forecasting of Hydrological Variables: Insights from a Simple Water Balance Model

Past work has shown that a land surface model’s (LSM’s) implicit (not explicitly coded) relationships between soil moisture and both evapotranspiration (ET) and runoff largely determine the LSM’s hydrological behavior. Here we estimate the relationships that appear to be operating in the real world and compare them to those of the LSM component of a state-of-the-art Earth system model (ESM). The two sets of relationships are determined by calibrating them within a simple water balance model (WBM): once using stream gauge observations from small, unregulated rivers over the eastern half of the U.S., and once using the runoffs generated by the LSM as part of a state-of-the-art atmospheric reanalysis. Hydrological simulations and subseasonal hydrological forecasts performed with the two calibrated versions of the WBM provide two key results. First, the version calibrated to the LSM-generated runoffs does successfully reproduce, to first order, the hydrological behavior of the full LSM within its ESM environment. Second, of the two WBM versions, the one calibrated to the observations reproduces more accurately a broad collection of fully independent streamflow observations as well as a similarly broad collection of in-situ soil moisture measurements. Taken together, the two results suggest that the observations-calibrated ET and runoff efficiency functions do successfully represent, at least to some degree, soil moisture controls over hydrological variability in Nature and can serve as potentially useful targets for further LSM development.

Water Balance Model↗
DOE OSTI2025

An Integrated Hydroclimatic Assessment of Future Reservoir and Hydropower Operations in the U.S.

The engineering of rivers by dams is a formative feature of human-nature systems and the interconnectivity of water, energy, and the climate. Sufficient and broad-based representations of dams in large-scale hydrological models prove essential to mapping their extensive regulation of river flow and biogeochemistry and gauging climate-linked provisions, including freshwater supply and hydropower. We present an integrated modeling framework to investigate future streamflow and hydropower generation in the Contiguous U.S. (1990–2075), leveraging an ensemble of six downscaled and bias-corrected General Circulation Models (GCMs) from the high-end SSP585 scenario of the CMIP6. To achieve this, we develop a reservoir operations and parameterization scheme for 1,384 dams in a high-resolution river network, including simulated hydropower generation for 326 dams. For the GCM ensemble mean, we simulate a widespread increase in regulated streamflow into the late-century (11% annual and 17% in winter for the dam median) with region-specific changes in summer streamflow that feature prominent declines in the Northwest (−7%). Mediation by reservoirs is shown to dampen intra-annual streamflow changes, delivering additional summer releases that partially mitigate declining flows. Total hydropower generation is projected to increase modestly (+3%), with boosted generation in the winter (+9%) and spring (+5%) offsetting declined summer generation (−3.4%), suggesting strong adaptation potential for hydropower in the future energy portfolio. Further analysis reveals that the choice of GCM, particularly in western regions, has significant bearing on projected streamflow and hydropower changes.

13 HYDRO ENERGY↗
DOE OSTI2025

Ensemble‐Based Spatially Distributed CLM5 Hydrological Parameter Estimation for the Continental United States

Abstract One of the major challenges in large‐domain hydrological modeling efforts lies in the estimation of spatially distributed hydrological parameters while simultaneously accounting for their associated uncertainties. Addressing this challenge is particularly difficult in ungauged locations. With growing societal demands for large‐scale streamflow projections to inform water resource management and long‐term planning, evaluating and constraining hydrological parameter uncertainty is increasingly vital. This study introduces a hybrid regionalization approach to enhance hydrological predictions of the Community Land Model version 5 (CLM5) across the Continental United States (CONUS), with a total of 50,629 1/8° grid cells. This hybrid method combines the strengths of two existing techniques: parameter regionalization and streamflow signature regionalization. It identifies ensemble behavioral parameters for each 1/8° grid cell across the CONUS domain, tailored to three distinct streamflow signatures focused on low flows, high flows, and annual water balance. Evaluating this hybrid method for 464 CAMELS (Catchment Attributes and Meteorology for Large‐sample Studies) basins demonstrates a significant improvement in CLM5 hydrological predictions, even in challenging arid regions. In CONUS applications, the derived spatially distributed parameter sets capture both spatial continuity and variation of parameters, highlighting their heterogeneous nature within specific regions. Overall, this hybrid regionalization approach offers a promising solution to the complex task of improving hydrological modeling over large domains for important hydrological applications.

CLM5↗
NASA NTRS2015

The Pattern Across the Continental United States of Evapotranspiration Variability Associated with Water Availability

The spatial pattern across the continental United States of the interannual variance of warm season water-dependent evapotranspiration, a pattern of relevance to land-atmosphere feedback, cannot be measured directly. Alternative and indirect approaches to estimating the pattern, however, do exist, and given the uncertainty of each, we use several such approaches here. We first quantify the water dependent evapotranspiration variance pattern inherent in two derived evapotranspiration datasets available from the literature. We then search for the pattern in proxy geophysical variables (air temperature, stream flow, and NDVI) known to have strong ties to evapotranspiration. The variances inherent in all of the different (and mostly independent) data sources show some differences but are generally strongly consistent they all show a large variance signal down the center of the U.S., with lower variances toward the east and (for the most part) toward the west. The robustness of the pattern across the datasets suggests that it indeed represents the pattern operating in nature. Using Budykos hydroclimatic framework, we show that the pattern can largely be explained by the relative strength of water and energy controls on evapotranspiration across the continent.

Air Temperature↗
NASA NTRS2019

Assimilation of Remotely Sensed Leaf Area Index into the Noah-MP Land Surface Model: Impacts on Water and Carbon Fluxes and States over the Continental U.S.

Accurate representation of vegetation states is required for the modeling of terrestrial water-energy-carbon exchanges and the characterization of the impacts of natural and anthropogenic vegetation changes on the land surface. This study presents a comprehensive evaluation of the impact of assimilating remote sensing-based Leaf Area Index (LAI) retrievals over the Continental U.S. in the Noah-MP land surface model, during a time period of 2000 to 2017. The results demonstrate that the assimilation has a beneficial impact on the simulation of key water budget terms such as soil moisture, evapotranspiration, snow depth, terrestrial water storage and streamflow, when compared with a large suite of reference datasets. In addition, the assimilation of LAI is also found to improve the carbon fluxes of Gross Primary Production (GPP) and Net Ecosystem Exchange (NEE). Most prominent improvements in the water and carbon variables are observed over the agricultural areas of the U.S., where assimilation improves the representation of vegetation seasonality impacted by cropping schedules. The systematic, added improvements from assimilation in a configuration that employs high quality boundary conditions highlight the significant utility of LAI data assimilation in capturing the impacts of vegetation changes.

anthropogenic vegetation↗
DOE OSTI2025

A Scale‐Adaptive Urban Hydrologic Framework: Incorporating Network‐Level Storm Drainage Pipes Representation

Abstract Below‐ground urban stormwater networks (BUSNs) significantly influence urban flood dynamics, yet their representation at the watershed or larger scales remains challenging. We introduce a scalable urban hydrologic framework that centers on a novel network‐level BUSN representation, balancing the needs for physical basis, parameter parsimony, and computational efficiency. Our framework conceptualizes an urban watershed into four interacting zones: hillslopes (natural), storm‐sewersheds (urban), a sub‐network channel (tributaries), and a main channel. We develop an innovative Graph Theory‐based algorithm to derive network‐level BUSN parameters from publicly available datasets, enabling efficient, scalable parameterization. We demonstrate this framework's applicability at nine representative watersheds in the Houston metropolitan region, USA, with urban imperviousness ranging from 0% to 64% and drainage areas ranging from 24 to 302 . Our model achieves satisfying computational efficiency, completing hourly time step simulations for 18 years in less than 5 sec per watershed on a standard PC. Validation against observed daily streamflow confirms that the model can capture small‐to‐large flood peaks and seasonal and annual water balance over these watersheds. Comparisons with the National Water Model show better performance in predicting flood peaks and overall water balance, underscoring the promises of our new framework for urban hydrologic modeling at large scales. Furthermore, analysis reveals nonlinear relationships between BUSNs' designed capacities and flood reduction effects. Our approach bridges the gap between detailed hydraulic and large‐scale hydrologic models, providing a valuable tool for urban flood prediction and management across broader spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Exploring Water System Vulnerabilities in California's Central Valley Under the Late Renaissance Megadrought and Climate Change

Abstract California faces cycles of drought and flooding that are projected to intensify, but these extremes may impact water users across the state differently due to the region's natural hydroclimate variability and complex institutional framework governing water deliveries. To assess these risks, this study introduces a novel exploratory modeling framework informed by paleo and climate‐change based scenarios to better understand how impacts propagate through the Central Valley's complex water system. A stochastic weather generator, conditioned on tree‐ring data, produces a large ensemble of daily weather sequences conditioned on drought and flood conditions under the Late Renaissance Megadrought period (1550–1580 CE). Regional climate changes are applied to this weather data and drive hydrologic projections for the Sacramento, San Joaquin, and Tulare Basins. The resulting streamflow ensembles are used in an exploratory stress test using the California Food‐Energy‐Water System model, a highly resolved, daily model of water storage and conveyance throughout California's Central Valley. Results show that megadrought conditions lead to unprecedented reductions in inflows and storage at major California reservoirs. Both junior and senior water rights holders experience multi‐year periods of curtailed water deliveries and complete drawdowns of groundwater assets. When megadrought dynamics are combined with climate change, risks for unprecedented depletion of reservoir storage and sustained curtailment of water deliveries across multiple years increase. Asymmetries in risk emerge depending on water source, rights, and access to groundwater banks.

Gupta, Rohini S. [School of Civil and Environmenta↗