Erica R. Siirila-Woodburn

Atmospheric warming is reducing snowpack, with uncertain effects on mountainous streamflow, a crucial water resource. Despite limited historical observations of groundwater–streamflow interactions above 2,500 m, new measurements in the Upper Colorado River headwaters indicate declining groundwater storage that is dated decades to millennia old. Here we use integrated hydrologic modelling spanning water years 2015–2021 to determine whether the loss of old-age groundwater buffers streamflow during low-snow years and whether that loss is exacerbated with warming. Results show that old-groundwater contributions to streams remain relatively steady through time, unlike the more variable contributions from young groundwater. Numerical experiments of increased surface air temperatures (+2.5 °C and +4 °C) increase rain–snow fractions and evapotranspiration and decrease runoff ratio by 2–3% per degree Celsius increase. As streamflow declines with warming, the age of groundwater supporting it gets older, in part owing to intermediate-aged (1–3 year) groundwater declining twice as fast. Simulations show that water table depths at higher elevations (>3,700 m) decline disproportionately and fail to recover even during wet years. These findings suggest altered groundwater–streamflow interactions with warming and snow loss, with implications for water resources.

Climate change will impact mountain watershed streamflow both directly—with changing precipitation amounts and variability—and indirectly—through temperature shifts altering snowpack, melt, and evapotranspiration. To understand how these complex processes will affect ecosystem functioning and water resources, we need tools to distinguish connections between water sources (rain/snowmelt), groundwater storage, and exit fluxes (streamflow/evapotranspiration), and to determine how these connections change seasonally and as climate shifts. Here, we develop novel watershed‐scale approaches to understand water source, storage, and exit flux connections using a dynamic‐flux particle tracking model (EcoSLIM) applied in California's Cosumnes Watershed, which connects the Sierra Nevada and Central Valley. This work develops new visualizations and applications to provide mechanistic understanding that underpins the interpretation of isotopic field data at watershed scales to distinguish sources, flow paths, residence times, and storage selection. In our simulations, streamflow comes primarily from snow‐derived water while evapotranspiration generally comes from rain. Most streamflow starts above 1,000 m while evapotranspiration is sourced relatively evenly across the watershed and is generally younger than streamflow. Modeled streamflow consists primarily of water sourced from precipitation in the previous 5 years but before the current water year, while ET consists primarily of water from precipitation in the current water year. ET, and to a lesser extent streamflow, are both younger than water in groundwater storage. However, snowmelt‐derived streamflow preferentially discharges older water from snow‐derived storage. Dynamic‐flux particle tracking and new approaches presented here enable novel model‐tracer comparisons in large‐scale watersheds to better understand watershed behavior in a changing climate.

Groundwater interactions with mountain streams are often simplified in model projections, potentially leading to inaccurate estimates of streamflow response to climate change. Here, using a high-resolution, integrated hydrological model extending 400 m into the subsurface, we find groundwater an important and stable source of historical streamflow in a mountainous watershed of the Colorado River. In a warmer climate, increased forest water use is predicted to reduce groundwater recharge resulting in groundwater storage loss. Losses are expected to be most severe during dry years and cannot recover to historical levels even during simulated wet periods. Groundwater depletion substantially reduces annual streamflow with intermittent conditions predicted when precipitation is low. Expanding results across the region suggests groundwater declines will be highest in the Colorado Headwater and Gunnison basins. Our research highlights the tight coupling of vegetation and groundwater dynamics and that excluding explicit groundwater response to warming may underestimate future reductions in mountain streamflow.

Groundwater storage and release in mountain systems is often ignored compounding uncertainty in how these systems may respond to climate change. To assess the importance of groundwater to mountain streamflow, we combine an extensive data network with a high-resolution, integrated hydrologic model stretching into the deep subsurface for a Colorado River headwater basin. In a warmer climate, increased forest water use reduces recharge, thereby destabilizing groundwater storage such that it cannot recover to historical levels even during simulated wet periods. Groundwater depletion significantly reduces average annual streamflow and forces the basin toward non-perennial conditions during dry water years. Groundwater declines are amplified in high-relief, conifer-dominated catchments. Upscaling model results across the region indicates concentrated risk of groundwater decline in the Colorado Headwater and Gunnison Basins. Research highlights the tight coupling of topography, vegetation and groundwater dynamics and that excluding groundwater response to warming will underestimate future reductions in mountain streamflow.

Wildfires naturally occur in many landscapes, however they are undergoing rapid regime shifts. Despite the emphasis in the literature on the most severe hydrological responses to wildfire, there remains a knowledge gap on the thresholds of wildfire (i.e., burned area/drainage area ratio, BAR) required to initiate hydrological responses. We investigated hydrological changes in the Russian River Watershed (RRW) in California, a coastal, Mediterranean, drought‐prone, wildfire‐adapted ecosystem, following ten wildfires that burned 30% of the watershed. Our findings suggest that sub‐watersheds of the RRW have not burned beyond an intrinsic, unknown, threshold required to initiate change. Using paired watersheds, we examined spatiotemporal patterns of pre‐and‐post wildfire hydrology with a rainfall‐runoff hydrological model. Even though these successive wildfires burned 1%–50% of each sub‐watershed (1%–30% at moderate/high severity), we found little evidence of wildfire‐related shifts in hydrology. As a function of BAR, wildfire imposed limited effects on runoff ratios (runoff/precipitation) and runoff residuals (observations—model simulations). Our findings that post‐wildfire runoff enhancements asymptote beyond 30% burn indicate that when a watershed is burned beyond a certain threshold, the magnitude of the hydrologic response no longer increases. Drought and storm conditions explained much of the variability observed in streamflow, whereas wildfire explained only moderate variability in streamflow even when wildfire accounted for >45% BAR. While the BAR in the RRW was sufficiently beyond previously reported minimum disturbance thresholds (>20% burned forest), the lack of hydrological response is attributed to buffering effects of wildfire adaptation and drought factors that are unique to Mediterranean ecoregions.

In the coastal regions of the western United States, atmospheric rivers (ARs) are associated with the largest precipitation generating storms and contribute up to half of annual precipitation, but the impact of ARs on the integrated hydrologic cycle, specifically on groundwater storage and hydrodynamics, is largely unknown. To better explore the hydrologic behavior of AR versus non‐AR event precipitation, we present a novel combination of two water tracking methods (one in the atmosphere and one in the subsurface) to explicitly track the full lifecycle of water parcels generated by ARs. Simulations of northern California's Cosumnes River watershed during the record wet 2017 water year are performed via the coupling of a high‐resolution regional climate model and a land surface‐groundwater model accounting for lateral groundwater flow. Despite ARs contributing more precipitation than non‐AR storms, we find less AR water is preferentially stored in aquifers by year end. Fractionally, ARs result in 300% less snow derived groundwater‐recharged compared to non‐AR precipitation. Rain‐on‐snow (RoS) plays an important role in AR‐driven discharge, where over 50% of total discharge from ARs snow is from RoS events. Finally, despite record‐breaking annual precipitation, simulated groundwater depletion occurs by year end due to estimates of groundwater pumping activities. The results from these simulations serve as a partial analogue of future hydrologic conditions where ARs are expected to intensify and provide a greater fraction of annual precipitation due to climate change.

Increasing wildfire and declining snowpacks in mountain regions threaten water availability. We combine satellite‐based fire detections with snow seasonality classifications to examine fire activity in California's seasonal and ephemeral snow zones. We find a nearly tenfold increase in fire activity during 2020–2021 versus 2001–2019. Accumulation season broadband snow albedo declined 25%–71% at two burned sites (2021 and 2022) according to in‐situ data relative to un‐burned conditions, with greater declines associated with increased burn severity. By enhancing snowpack susceptibility to melt, both decreased snow albedo and canopy drove midwinter melt during a multi‐week dry spell in 2022. Despite similar meteorological conditions in December–February 2013 and 2022–linked to persistent high pressure weather regimes–minimal melt occurred in 2013. Post‐fire snowpack differences are confirmed with satellite measurements. With growing geographical overlap between wildfire and snow, our findings suggest California's snowpack is increasingly vulnerable to the compounding effects of dry spells and wildfire.