Search NASA⌕ Search

SEARCH · Search NASA

Results for “Runoff”

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 145 records · Page 8

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗

Significant Rural Contributions to Severe Coastal Urban Flooding Through Flow Connectivity

Coastal cities are increasingly threatened by hurricane-induced compound flooding. However, the relative contributions of rainfall-surge interactions and other flood drivers, such as runoff from surrounding rural and local impervious areas, in shaping urban flooding remain unclear. Here, we use a multi-scale Earth system modeling framework to disentangle the contributions of diverse flooding drivers. The results show that urban flooding is primarily affected by topographic factors (low elevation, flat slopes, and localized flow constrictions), rather than increased runoff from their impervious surfaces. Importantly, runoff from rural areas contributes more to severely flooded urban areas with inundation depth exceeding 1.7 m than local runoff. Although rainfall–surge interactions and sea-level rise exacerbate inundation along the coastline, over 90% of the compounding effects are absorbed by wetlands. By resolving flooding processes across scales and components, this study emphasizes the role of runoff from surrounding rural areas in shaping zones of severe urban flooding.

58 GEOSCIENCES↗

The impact of tillage row orientation on physical and chemical sediment enrichment

Abstract This study aimed to better understand how tillage row orientation with respect to dominant flow‐pathway along hillslope impacts runoff and the transport of different sediment size fractions. Experimental plots were constructed in contour ridge till (CRT) and parallel ridge till (PRT) sites to monitor runoff and sediment fluxes. Particle size fractions of the sediment, along with organic C and N contents were measured to quantify physical and chemical enrichment ratios for the two tillage orientations. In the CRT plot, tillage produced large oriented roughness elements along the contours, which acted as little check dams, while in the PRT site, the roughness elements helped confine and concentrate the runoff. In the CRT, the “check dams” resulted in a runoff coefficient of .03, while in the PRT, the flow confinement between rows produced a runoff coefficient of .82. Moreover, the erosion rates at the CRT site were 97% less than those in the PRT plot. Large contours produced finer sediment fractions due to the selective sorting in ponded furrows. More aggregate sediment fractions were present in the PRT site, which was dominated by rill erosion. Physical enrichment revealed selective entrainment of finer sediment particles during erosion. Finer‐grain particles with higher specific surface areas, attached more organic C resulting in chemical enrichment. Physical and chemical enrichment methods were found to be in good agreement. These findings suggest transport models that can simulate size fraction updates to the soil active layer can be used to estimate SOC redistribution and hence more accurate C budgets.

Wacha, Kenneth M.↗

Aridification of Colorado River Basin's Snowpack Regions Has Driven Water Losses Despite Ameliorating Effects of Vegetation

Abstract The Colorado River Basin is an important natural resource for the semi‐arid southwestern United States (US), where it provides water to more than 40 million people. While nearly 1.5°C of anthropogenic warming has occurred across this region from the 1880s to 2021, climate models show little agreement in the precipitation change during the same historical period, with no trend in the mean of the latest (sixth) generation of Global Climate Models. As such, here we focus on how the CO 2 increase and associated anthropogenic warming over the historical period has impacted runoff across the Colorado Basin. We find that the Colorado Basin's runoff over the historical period has decreased by 8.1% per degree Celsius of warming (°C −1 ). However, the magnitude of this sensitivity is reduced to 6.8% °C −1 when considering vegetation response to historical CO 2 . For present‐day conditions, this translates to runoff reductions of 10.3% due to anthropogenic increases in both temperature and CO 2 since 1880. We demonstrate that Colorado Basin's natural flow has been decreased by roughly the storage of Lake Mead during the 2000–2021 megadrought due to this long term anthropogenic influence, suggesting the basin's first shortage in 2021 would likely not have occurred without anthropogenic warming. We further show warming has led to disproportionate aridification in snowpack regions, causing runoff to decline at double the rate relative to non‐snowpack regions. Thus, despite only making up ∼30% of the basin's drainage area, 86% of runoff decreases in the Colorado Basin is driven by water loss in snowpack regions.

54 ENVIRONMENTAL SCIENCES↗

Divergent responses of historic rain-on-snow flood extremes to a warmer climate

Global warming is altering flood risks induced by rain-on-snow events. However, decision-makers lack guidance on how rain-on-snow induced extreme floods could be altered with warming. Here, storyline analyses using a kilometer-scale land surface model reveal diverse responses of four historically-impactful, decision-relevant rain-on-snow induced extreme flood events over the contiguous U.S. to warming, due to alterations in their water budgets. For the 2017-Feb California floods, runoff first increases and then decreases with warming, peaking under the +3 K scenario, while runoff of the 2017-Jan California floods increases monotonically by ~53%/K. Contrastingly, runoff of the 1996-Jan Mid-Atlantic floods decreases gradually with warming. Despite these differences, warming generally shifts flood-generating regimes along elevation profiles. High elevations could experience notably increased runoff, while low elevations encounter a shift from rain-on-snow-driven to rainfall-dominated runoff. These findings underscore the need for flood control planning to quantify region- and elevation-specific changes in rain-on-snow events in a warmer climate.

54 ENVIRONMENTAL SCIENCES↗

Water balance response of permafrost-affected watersheds to changes in air temperatures

Abstract Observations show increases in river discharge to the Arctic Ocean especially in winter over the last decades but the physical mechanisms driving these changes are not yet fully understood. We hypothesize that even in the absence of a precipitation increase, permafrost degradation alone can lead to increased annual river runoff. To test this hypothesis we perform 12 millennium-long simulations over an idealized hypothetical watershed (1 km 2 ) using a distributed, physically based water balance model (Water flow and Balance Simulation Model, WaSiM). The model is forced by both a hypothetical warming defined by an air temperature increase of 7.5 ∘ C over 100 years, and a corresponding cooling scenario. To assess model sensitivity we vary soil saturated hydraulic conductivity and lateral subsurface flow configuration. Under the warming scenario, changes in subsurface water transport due to ground temperature changes result in a 7%–14% increase in annual runoff accompanied by a 6%–20% decrease in evapotranspiration. The increase in runoff is most pronounced in winter. Hence, the simulations demonstrate that changes in permafrost characteristics due to climate warming and associated changes in evapotranspiration provide a plausible mechanism for the observed runoff increases in Arctic watersheds. In addition, our experiments show that when lateral subsurface moisture transport is not included, as commonly done in global-scale Earth System Models, the equilibrium water balance in response to the warming or cooling is similar to the water balance in simulations where lateral subsurface transport is included. However, the transient changes in water balance components prior to reaching equilibrium differ greatly between the two. For example, for high saturated hydraulic conductivity only when lateral subsurface transport is considered, a period of decreased runoff occurs immediately after the warming. This period is characterized by a positive change in soil moisture storage caused by the soil moisture deficit developed during prior cooling.

Debolskiy, Matvey V. (ORCID:0000000296343627)↗

To What Extent Can Vegetation Mitigate Greenhouse Warming? A Modeling Approach

Climate models participating in the IPCC Fourth Assessment Report indicate that under a 2xCO2 environment, runoff would increase faster than precipitation overland. However, observations over large U.S watersheds indicate otherwise. This inconsistency suggests that there may be important feedbacks between climate and land surface unaccounted for in the present generation of models. We postulate that the increase in precipitation associated with the increase in CO2 is also increasing vegetation density, which may already be feeding back onto climate. Including this feedback in a climate model simulation resulted in precipitation and runoff trends consistent with observations and reduced the warming by 0.6OC overland. This unaccounted for missing water may be linked to about 10% of the missing land carbon sink. A recent compilation of outputs from 19 coupled atmosphere-ocean general circulation models used in the IPCC Fourth Assessment Report (AR4) shows projected increases in air temperature, precipitation and river discharge for 24 major rivers in the world in response to doubling CO2 by the end of the century (1). The ensemble mean from these models also indicates that, compared to their respective baselines overland, the global mean of the runoff change would increase faster (8.9% per year) than that of the precipitation (5% per year). We analyze century-scale observed annual runoff time-series (1901-2002) over 9 hydrological units covering large regions of the Eastern United States (Fig.1) compiled by the United States Geological Survey (USGS)(2). These regions were selected because they are the most forested; the least water-limited and are not under extensive irrigation. We compare these time-series to similar time-series of observed annual precipitation anomalies spanning the period 1900-1995 (3). Both time-series exhibit a positive longterm trend (Fig. 2); however, in contrast to the analysis of (I), these historic data records show that the rate of precipitation increase is 5.5 % per year, roughly double the rate of runoff increase of 3.1 % per year.

Bounoua, L.↗

Century Scale Evaporation Trend: An Observational Study

Several climate models with different complexity indicate that under increased CO2 forcing, runoff would increase faster than precipitation overland. However, observations over large U.S watersheds indicate otherwise. This inconsistency between models and observations suggests that there may be important feedbacks between climate and land surface unaccounted for in the present generation of models. We have analyzed century-scale observed annual runoff and precipitation time-series over several United States Geological Survey hydrological units covering large forested regions of the Eastern United States not affected by irrigation. Both time-series exhibit a positive long-term trend; however, in contrast to model results, these historic data records show that the rate of precipitation increases at roughly double the rate of runoff increase. We considered several hydrological processes to close the water budget and found that none of these processes acting alone could account for the total water excess generated by the observed difference between precipitation and runoff. We conclude that evaporation has increased over the period of observations and show that the increasing trend in precipitation minus runoff is correlated to observed increase in vegetation density based on the longest available global satellite record. The increase in vegetation density has important implications for climate; it slows but does not alleviate the projected warming associated with greenhouse gases emission.

Bounoui, Lahouari↗

Oceanic Transport of Surface Meltwater from the Southern Greenland Ice Sheet

The Greenland ice sheet has undergone accelerating mass losses during recent decades. Freshwater runoff from ice melt can influence fjord circulation and dynamic1 and the delivery of bioavailable micronutrients to the ocean. It can also have climate implications, because stratification in the adjacent Labrador Sea may influence deep convection and the strength of the Atlantic meridional overturning circulation. Yet, the fate of the meltwater in the ocean remains unclear. Here, we use a high-resolution ocean model to show that only 1-15% of the surface meltwater runoff originating from southwest Greenland is transported westwards. In contrast, up to 50-60% of the meltwater runoff originating from southeast Greenland is transported westwards into the northern Labrador Sea, leading to significant salinity and stratification anomalies far from the coast. Doubling meltwater runoff, as predicted in future climate scenarios, results in a more-than-double increase in anomalies offshore that persists further into the winter. Interannual variability in offshore export of meltwater is tightly related to variability in wind forcing. The new insight that meltwaters originating from the west and east coasts have different fates indicates that future changes in mass loss rates and surface runoff will probably impact the ocean differently, depending on their Greenland origins.

Physical oceanography↗

The Potential for Snow to Supply Human Water Demand in the Present and Future

Runoff from snowmelt is regarded as a vital water source for people and ecosystems throughout the Northern Hemisphere (NH). Numerous studies point to the threat global warming poses to the timing and magnitude of snow accumulation and melt. But analyses focused on snow supply do not show where changes to snowmelt runoff are likely to present the most pressing adaptation challenges, given sub-annual patterns of human water consumption and water availability from rainfall. We identify the NH basins where present spring and summer snowmelt has the greatest potential to supply the human water demand that would otherwise be unmet by instantaneous rainfall runoff. Using a multi-model ensemble of climate change projections, we find that these basins - which together have a present population of approx. 2 billion people - are exposed to a 67% risk of decreased snow supply this coming century. Further, in the multi-model mean, 68 basins (with a present population of more than 300 million people) transition from having sufficient rainfall runoff to meet all present human water demand to having insufficient rainfall runoff. However, internal climate variability creates irreducible uncertainty in the projected future trends in snow resource potential, with about 90% of snow-sensitive basins showing potential for either increases or decreases over the near-term decades. Our results emphasize the importance of snow for fulfilling human water demand in many NH basins, and highlight the need to account for the full range of internal climate variability in developing robust climate risk management decisions.

Mankin, Justin S.↗

Phosphorus and cover crop management practices affect phosphorus speciation in soils and eroded sediments

Abstract Agricultural runoff often contains P in dissolved and sediment‐bound forms, decreasing surface water quality. No‐till and cover cropping conservation practices have been recommended for reducing erosion and nutrient loss from cropping systems. The overall aims of this study were to characterize and evaluate the effects of fertilizer (placement and source) and cover crop management on P speciation in surface runoff sediments and source soil. In 2014, a field‐scale experiment was established in a no‐till, corn (Zea maysL.)–soybean (Glycine maxL.) cropping system with two cover crop treatments (with and without a winter crop; winter wheat [Triticum aestivumL.], rapeseed [Brassica napusL.], hairy vetch [Vicia villosaRoth], winter triticale [×Triticosecale Wittm.], and cereal rye [Secale cerealeL.]) and three P fertilizer management treatments (no P, fall broadcast diammonium phosphate, and spring subsurface injected ammonium polyphosphate). Phosphorus fractionation in the source soil collected in the fall of 2019 and sediment samples collected throughout 2020 were analyzed using a modified sequential P extraction method to evaluate the cumulative effects of imposing the treatment factors over 5 years. The direct P speciation was done using X‐ray absorption near edge structure spectroscopy. The indirect P speciation (fractionation) results showed that the management practices affected the exchangeable, organic matter‐associated, and Fe‐bound P fractions in sediments and the exchangeable and residual fractions in source soil. Direct P speciation results showed a depletion of Fe‐associated P in soil and sediment from cover crop treatment, suggesting that Fe‐associated P species were affected by cover crops. Changes in soil and runoff sediment P speciation would change the proportions and forms of soluble and particulate P in runoff sediments and may influence P bioavailability in aquatic ecosystems. Developing P fertilizer and cropping system management options with an understanding of soil P transformations helps maintain environmental sustainability.

Environmental Sciences & Ecology↗

Multi‐year hydroclimatic droughts and pluvials across the conterminous United States

Abstract Time series of water‐year runoff for 2,109 hydrologic units (HUs) across the conterminous United States (CONUS) for the 1900 through 2014 period were used to identify drought and pluvial (i.e., wet) periods. Characteristics of the drought and pluvial events including frequency, duration, and severity were examined and compared. Additionally, a similar analysis was performed using gridded tree‐ring reconstructions of the Palmer Drought Severity Index (PDSI) for the period 1475 through 2005 to place the drought and pluvial characteristics determined using water‐year runoff for 1900 through 2014 in the context of multi‐century climate variability. The temporal and spatial variability of droughts and pluvials determined using runoff for the 1900 through 2014 period indicated that most drought events in the CONUS occurred before about 1970, whereas most pluvial periods occurred after about 1970. This change in the frequencies of drought and pluvial events around 1970 was largely related to an increase in fall (October through December) precipitation across much of the central United States. Also, the duration and severity of droughts and pluvials identified using runoff for the 1900 through 2014 period generally were not significantly different from the drought and pluvial characteristics identified using the PDSI for the 1475 through 2005 period.

McCabe, Gregory J.↗

Role of underlying surface, rainstorm and antecedent wetness condition on flood responses in small and medium sized watersheds in the Yangtze River Delta region, China

Floods are becoming more frequent and are hard to control due to the shortage of water conservancy projects in small and medium sized basins, especially in developing regions. Understanding the hydrologic responses and estimated flood characteristics to storm events can help to predict flood disasters and form better effective mitigation and adaptation strategies. Here, using observations of several representative watersheds and Artificial Neural Networks (ANN) and Principal Component Regression (PCR) models, we conducted data-driven analyses to examine the flood responding characteristics. We quantified the relative contributions of the influencing factors to the variation of each flood characteristic. Statistical analysis of the observations shows that the drainage area plays a key role in determining the distribution of lag time and peak discharge. Rainstorm variability has direct influence on floods, and typhoon-induced rainstorms with high total rainfall and rainfall intensity generate higher lag time, flood peak, unit discharge and runoff depth, but lower runoff coefficient. The ANN and PCR models accurately predicted the variations of flood features using the driving factors including physical geographical characteristics, rainstorm features, and antecedent condition. Physical geographic characteristics are key influential factors of lag time, flood peak and runoff coefficient, while the rainstorm features control the magnitude of unit discharge and runoff depth. These results indicate that floods are mainly affected by rainstorm features and physical geographic characteristics in the Yangtze River delta, and it might become more damaging with the increasing rainfall extremes and sprawling impervious surfaces in a changing environment.

54 ENVIRONMENTAL SCIENCES↗

Coupled effects of future rainfall and land use on urban stormwater drainage system in Tampa, Florida (USA)

Future rainfall and land use are two important factors for flood management since they both directly and indirectly affect the functionality and performance of urban stormwater drainage infrastructures. Although some studies have already examined the coupled effects of future rainfall and land use conditions, they have not concluded how such effects on both surface runoff and flood areas would change with regional surface elevation variation and infrastructure conditions. This paper analyzed the coupled effects of future rainfall and land use on urban drainage systems in terms of surface runoff quantity and flood area changes using EPA SWMM. Future downscaled and bias-corrected precipitation projections for 2040–2060 and 2080–2099 from “Downscaled CMIP3 and CMIP5 Climate and Hydrology Projections” Archive and future land use for 2050 s and 2090 s from EPA ICLUS project were used for the City of Tampa in Florida (USA). It was shown that potential runoff volume and flood area changes ranged from –40% to 160%, and –40% to 400% due to rainfall change, while potential changes induced by land use change ranged from 0% to 3.5%, and 0% to 18% respectively. Additionally, this study examined the coupled effects of future rainfall and land use changes on flood area variations, considering the capacity of drainage infrastructure and elevation variation of study site. Results showed that flood area variations are not solely influenced by runoff quantity changes due to the two external drivers, but also by other factors. Specifically, in low-relief areas with inadequate drainage infrastructure and less elevation variation, these two drivers have no additive effects on the percentage of flood area changes. In contrast, in areas with adequate drainage infrastructure and greater surface elevation variation, there are additive effects from rainfall and land use changes on flood area expansions. These findings have important implications for stakeholders involved in city planning and flooding control, particularly for urban areas facing future precipitation and urbanization challenges.

54 ENVIRONMENTAL SCIENCES↗

Using SWMM for emergency response planning: A case study evaluating biological agent transport under various rainfall scenarios and urban surfaces

To assist in emergency preparedness for a biological agent terrorist attack or accidental pathogen release, potential contaminant levels and migration pathways of spores spread by urban stormwater were evaluated using a Storm Water Management Model (SWMM) of U.S. Coast Guard Base Elizabeth City, North Carolina. The high temporal-spatial resolution SWMM model was built using spore concentrations in stormwater runoff from asphalt, grass, and concrete collected from a point-scale field study. The subsequent modeled contamination scenarios included a notional plume release and point releases mimicking the field study under three rainfall conditions. The rainfall scenarios included a 6-hour natural rainfall event on Dec. 8, 2021 and two design storms (2-year and 100-year events). The observed spore concentrations from asphalt and concrete from the actual field experiment were applied to calibrate the washoff parameters in the SWMM model, using an exponential washoff function. The calibrated washoff coefficient (c 1 ) and exponent (c 2 ) were 0.01 and 1.00 for asphalt, 0.05 and 1.45 for grass, and 2.45 and 1.00 for concrete, respectively. The calibrated SWMM model simulated spore concentrations in runoff at times and magnitudes similar to the field study data. In the point release modeled scenario, the concrete surface generated 55.6% higher average spore concentrations than asphalt. Similarly, in the field experiment, a 175% (p < 0.05) higher average spore concentration in surface runoff was observed from concrete than from asphalt. Here, this study demonstrates how SWMM may be used to evaluate spore washoff from urban surfaces under different precipitation amounts, intensities, and durations, and how visualized spatial migration pathways in stormwater runoff may be used for emergency planning and remediation.

54 ENVIRONMENTAL SCIENCES↗

Future Warming and Intensification of Precipitation Extremes: A “Double Whammy” Leading to Increasing Flood Risk in California

This study focuses on quantifying future anthropogenic changes in surface runoff associated with extreme precipitation in California's Sierra Nevada. The method involves driving a land surface model with output from a high resolution regional atmospheric simulation of the most extreme atmospheric rivers (ARs). AR events were selected from an ensemble of global climate model simulations of historical and late 21st century climate under the “high-emission” RCP8.5 scenario. Average precipitation during the future ARs increases by ~25% but a much lower proportion falls as snow. The resulting future runoff increase is dramatic—nearly 50%, reflecting both the precipitation increase and simultaneous conversion of snow to rain. The “double whammy” impact on runoff is largest in the 2,000–2,500 m elevation band, where the snowfall loss and precipitation increase are both especially large. Overall, this huge increase in runoff during the most extreme AR events could present major flood control challenges for the region.

54 ENVIRONMENTAL SCIENCES↗

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in situ and remote sensing data sets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

Cheng, Yanyan↗

Causes of Missing Snowmelt Following Drought

Abstract Water management in snowy mountainous regions hinges on forecasting snowmelt runoff. However, droughts are altering snowpack‐runoff relationships with ongoing debate about the driving mechanisms. For example, in 2021 in California, less than half of predicted streamflow arrived. Mechanisms proposed for this “missing” streamflow included changes in evapotranspiration ( ET ), rainfall, and subsurface moisture conditions. Here, we demonstrate that ET in drought years generates dry subsurface conditions that reduce runoff in subsequent years. A model including this legacy of depleted moisture storage reduced median error in 2021 forecasts from 60% to 20% at 15 minimally disturbed basins and from 18% to 2% at 6 water supply basins in the Sierra Nevada (basins range in area from 5 to 23,051 km 2 and mean annual precipitation from 814 to 1,549 mm). Our findings indicate that the relationship between snowpack and runoff will evolve as plant ecosystems respond to climate change and alter subsurface water storage dynamics.

Lapides, Dana A.↗