Under anthropogenic climate change (CC), the global hydrological cycle intensifies at a rate known as hydrologic sensitivity (HS). Global climate models (GCMs) exhibit substantial uncertainty in HS. Past work suggests that another form of HS, derived from internal climate variability (IV), is useful for constraining this uncertainty. However, these two forms of HS are weakly related. Here we show that decomposing HS under both CC and IV, based on the global energy budget, provides insight into the likely range of future HS. We find that sensible heat exchange between the atmosphere and ocean is not accounted for in the atmospheric energy budget under IV, masking the connection between HS under IV and CC. Removing this term, a closer relationship emerges. We use observations in conjunction with this relationship to suggest an upward shift in the likely range of future HS (66% confidence interval: 2.00–2.36 W m⁻² K⁻¹).

Global climate models (GCMs) are unable to produce detailed runoff conditions at the basin scale. Assumptions are commonly made that dynamical downscaling can resolve this issue. However, given the large magnitude of the biases in downscaled GCMs, it is unclear whether such projections are credible. Here, we use an ensemble of dynamically downscaled GCMs to evaluate this question in the Sierra-Cascade mountain range of the western US. Future projections across this region are characterized by earlier seasonal shifts in peak flow, but with substantial inter-model uncertainty (−25 ± 34.75 days, 95% confidence interval (CI)). We apply the emergent constraint (EC) method for the first time to dynamically downscaled projections, leading to a 39% (−28.25 ± 20.75 days, 95% CI) uncertainty reduction in future peak flow timing. While the constrained results can differ from bias corrected projections, the EC is based on GCM biases in historical peak flow timing and has a strong physical underpinning.

Localized projections of 21st-century hydroclimate variables obtained from downscaling Global Climate Model (GCM) output are central to informing regional impact assessments and infrastructure planning. Regional GCM biases can be significant and, for dynamical downscaling, can be addressed either before (a priori) or after (a posteriori) downscaling. However, a priori bias correction (APBC) has generally unexplored effects on climate change signals. Here we analyze dynamically downscaled solutions of CMIP6 GCMs over the Western U.S., with and without APBC, and quantify APBC's impact on climate change signals relative to other irreducible uncertainty sources. For temperature and precipitation, the uncertainty introduced by APBC is negligible compared to that arising from GCM choice or internal variability. Furthermore, APBC greatly reduces regional models' unrealistically high snow-water-equivalent (SWE) biases that result directly from GCM errors. We leverage this finding to encourage the dynamical downscaling community to adopt APBC as a standard operating procedure.

Polarization surrounding bias correction (BC) in creating climate projections arises from its lack of physicality. Here, we perform and analyze 18 dynamical downscaling simulations (with and without BC) to better understand the physical impacts of BC, applied before downscaling, on regional climate output across the western United States. Without BC, downscaled precipitation is systematically and unrealistically wet biased compared to a hierarchy of observationally based datasets over the 1980–2014 period due to cascading mean-state Global Climate Model (GCM) biases: (a) overly strong lower-tropospheric lapse rates (5 K/km), (b) overly cold (2 K) tropospheric temperatures, and (c) anomalous mid-tropospheric cyclonic vorticity advection. With BC, downscaled precipitation (snow) biases are virtually eliminated (halved). Identified GCM biases are common to the broader Coupled Model Intercomparison Project ensemble. Physical effects of BC on the quality of the regionalized projections, pending an evaluation of BC's distortion of the downscaled climate response, may motivate its broader application by dynamical downscalers.

Predicting future climate change over a region of complex terrain, such as the western United States (US), remains challenging due to the low resolution of global climate models (GCMs). Yet the climate extremes of recent years in this region, such as floods, wildfires, and drought, are likely to intensify further as climate warms, underscoring the need for high-quality and high-resolution predictions. Here, we present an ensemble of dynamically downscaled simulations over the western US from 1980–2100 at 9 km grid spacing, driven by 16 latest-generation GCMs. This dataset is titled the Western US Dynamically Downscaled Dataset (WUS-D3).

We describe the challenges of producing WUS-D3, including GCM selection and technical issues, and we evaluate the simulations' realism by comparing historical results to temperature and precipitation observations. The future downscaled climate change signals are shaped in physically credible ways by the regional model's more realistic coastlines and topography. (1) The mean warming signals are heavily influenced by more realistic snowpack. (2) Mean precipitation changes are often consistent with wetting on the windward side of mountain complexes, as warmer, moister air masses are uplifted orographically during precipitation events. (3) There are large fractional precipitation increases on the lee side of mountain complexes, leading to potentially significant changes in water resources and ecology in these arid landscapes. (4) Increases in precipitation extremes are generally larger than in the GCMs, driven by locally intensified atmospheric updrafts tied to sharper, more realistic gradients in topography. (5) Changes in temperature extremes are different from what is expected by a shift in mean temperature and are shaped by local atmospheric dynamics and land surface feedbacks. Because of its high resolution, comprehensiveness, and representation of relevant physical processes, this dataset presents a unique opportunity to evaluate societally relevant future changes in western US climate.

The increasing prevalence of low snow conditions in a warming climate has attracted substantial attention in recent years, but a focus exclusively on low snow leaves high snow years relatively underexplored. However, these large snow years are hydrologically and economically important in regions where snow is critical for water resources. Here, we introduce the term “snow deluge” and use anomalously high snowpack in California’s Sierra Nevada during the 2023 water year as a case study. Snow monitoring sites across the state had a median 41 y return interval for April 1 snow water equivalent (SWE). Similarly, a process-based snow model showed a 54 y return interval for statewide April 1 SWE (90% CI: 38 to 109 y). While snow droughts can result from either warm or dry conditions, snow deluges require both cool and wet conditions. Relative to the last century, cool-season temperature and precipitation during California’s 2023 snow deluge were both moderately anomalous, while temperature was highly anomalous relative to recent climatology. Downscaled climate models in the Shared Socioeconomic Pathway-370 scenario indicate that California snow deluges—which we define as the 20 y April 1 SWE event—are projected to decline with climate change (58% decline by late century), although less so than median snow years (73% decline by late century). This pattern occurs across the western United States. Changes to snow deluge, and discrepancies between snow deluge and median snow year changes, could impact water resources and ecosystems. Understanding these changes is therefore critical to appropriate climate adaptation.

The serial occurrence of atmospheric rivers (ARs) along the US West Coast can lead to prolonged and exacerbated hydrologic impacts, threatening flood-control and water-supply infrastructure due to soil saturation and diminished recovery time between storms. Here a statistical approach for quantifying subseasonal temporal clustering among extreme events is applied to a 41-year (1979–2019) wintertime AR catalog across the western United States (US). Observed AR occurrence, compared against a randomly distributed AR timeseries with the same average event density, reveals temporal clustering at a greater-than-random rate across the western US with a distinct geographical pattern. Compared to the Pacific Northwest, significant AR clusters over the northern Coastal Range of California and Sierra Nevada are more frequent and occur over longer time periods. Clusters along the California Coastal Range typically persist for 2 weeks, are composed of 4–5 ARs per cluster, and account for over 85% of total AR occurrence. Across the northwest Coast-Cascade Ranges, clusters account for ∼50% of total AR occurrence, typically last 8–10 days, and contain 3–4 individual AR events. Based on precipitation data from a high-resolution dynamical downscaling of reanalysis, the fractions of total and extreme hourly precipitation attributable to AR clusters are largest along the northern California coast and in the Sierra Nevada. Interannual variability among clusters highlights their importance for determining whether a particular water year is anomalously wet or dry. The mechanisms behind this unusual clustering are unclear and require further research.

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₂ 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⁻¹). However, the magnitude of this sensitivity is reduced to 6.8% °C⁻¹ when considering vegetation response to historical CO₂. For present-day conditions, this translates to runoff reductions of 10.3% due to anthropogenic increases in both temperature and CO₂ 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.

Downslope wind-driven fires have resulted in many of the wildfire disasters in the western United States and represent a unique hazard to infrastructure and human life. We analyze the co-occurrence of wildfires and downslope winds across the western United States (US) during 1992–2020. Downslope wind-driven fires accounted for 13.4% of the wildfires and 11.9% of the burned area in the western US yet accounted for the majority of local burned area in portions of southern California, central Washington, and the front range of the Rockies. These fires were predominantly ignited by humans, occurred closer to population centers, and resulted in outsized impacts on human lives and infrastructure. Since 1999, downslope wind-driven fires have accounted for 60.1% of structures and 52.4% of human lives lost in wildfires in the western US. Downslope wind-driven fires occurred under anomalously dry fuels and exhibited a seasonality distinct from other fires—occurring primarily in the spring and fall. Over 1992–2020, we document a 25% increase in the annual number of downslope wind-driven fires and a 140% increase in their respective annual burned area, which partially reflects trends toward drier fuels. These results advance our understanding of the importance of downslope winds in driving disastrous wildfires that threaten populated regions adjacent to mountain ranges in the western US. The unique characteristics of downslope wind-driven fires require increased fire prevention and adaptation strategies to minimize losses and incorporation of changing human-ignitions, fuel availability and dryness, and downslope wind occurrence to elucidate future fire risk.

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.

Large uncertainty exists in hydrologic sensitivity (HS), the global-mean precipitation increase per degree of warming, across global climate model (GCM) ensembles. Meanwhile, the global circulation and hence global precipitation are sensitive to variations of surface temperature under internal variability. El Niño–Southern Oscillation (ENSO) is the most dominant mode of global temperature variability and hence of precipitation variability. Here we show in phase 6 of the Coupled Model Intercomparison Project (CMIP6) that the strength of HS under ENSO is predictive of HS in the climate change context (r = 0.56). This correlation increases to 0.62 when only central Pacific ENSO events are considered, suggesting that they are a better proxy for HS under future warming than east Pacific ENSO events. GCMs with greater HS are associated with greater weakening of the Walker circulation and expansion of the Hadley circulation under ENSO. Observations of HS under ENSO suggest that it is significantly underestimated by the GCMs, with the lower bound of observational uncertainty almost double even the highest-HS GCMs. The ENSO-related transformation of the tropical circulation holds clues into how the GCMs may be improved in order to more reliably simulate future hydrological cycle intensification.

Throughout the world, the hydrologic cycle is projected to become more variable due to climate change, posing challenges in semi-arid regions with high water resource vulnerability. Precipitation whiplash results from hydrologic variability, and refers to interannual shifts between wet (⩾80th historical percentile) and dry (⩽20th historical percentile) years. Using five model large ensembles, we show that whiplash is projected to increase in frequency (25%–60%) and intensity (30%–100%) by 2100 across several semi-arid regions of the globe, including Western North America and the Mediterranean. These changes can be driven by increases in the frequency of wet years or dry years, or both, depending on the region. Moisture budget calculations in these regions illuminate the physical mechanisms behind increased whiplash. Thermodynamic changes generally dominate, with modulations by dynamics, evaporation, and eddies on regional or global scales. These findings highlight increasingly volatile hydrology in semi-arid regions as the 21st Century progresses.

Future projections of global meteorological drought are evaluated in the Multi-Model Large Ensemble Archive, including an evaluation of the atmospheric moisture budget, conditioned on drought years. Drought is defined as 5-year running-mean annual precipitation below some threshold, for example, 10th percentile. Drought increases in frequency over the subtropics, in addition to certain tropical regions, consistent with previous studies. The moisture-budget decomposition allows drought to be defined as mean-flow, eddy, or feedback droughts, depending on which term in the equation contributes the largest negative interannual anomaly. In the historical climate, mean-flow droughts constitute most droughts at low latitudes; eddy droughts are equally common at higher latitudes; feedback droughts (i.e., droughts exacerbated by land–atmosphere feedbacks) constitute almost all droughts in water-limited subtropical/Mediterranean regions. The future drought increases are predominantly due to increases in feedback droughts in regions where these droughts are common historically but also over the Amazon. However, over most Mediterranean-type regions mean-flow droughts are also large contributors, resulting from dynamics. Eddy droughts also contribute to future increases along the equatorward flanks of historical eddy-driven jets, likely reflecting poleward shifts therein. Model uncertainty is particularly large over the Amazon and Australia, a reflection of model diversity in processes associated with land-atmosphere interaction. Based on these results, an availability of 3-D atmospheric data from a wider swath of global climate model large ensembles could help constrain global drought projections based on the representation of drought mechanisms in the historical climate.

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. 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.

A key indicator of climate change is the greater frequency and intensity of precipitation extremes across much of the globe. In fact, several studies have already documented increased regional precipitation extremes over recent decades. Future projections of these changes, however, vary widely across climate models. Using two generations of models, here we demonstrate an emergent relationship between the future increased occurrence of precipitation extremes aggregated over the globe and the observable change in their frequency over recent decades. This relationship is robust in constraining frequency changes to precipitation extremes in two separate ensembles, and under two future emissions pathways (reducing intermodel spread by 20-40%). Moreover, this relationship is also apparent when the analysis is limited to near-global land regions. These constraints suggest that historical global precipitation extremes will occur roughly 32 ± 8% more often than present by 2100 under a medium-emissions pathway (and 55 ± 13% under high-emissions).

Southern California is a biodiversity hotspot and home to over 23 million people. Over recent decades the annual wildfire area in the coastal southern California region has not significantly changed. Yet how fire regime will respond to future anthropogenic climate change remains an important question. Here, we estimate wildfire probability in southern California at station scale and daily resolution using random forest algorithms and downscaled earth system model simulations. We project that large fire days will increase from 36 days/year during 1970–1999 to 58 days/year under moderate greenhouse gas emission scenario (RCP4.5) and 71 days/year by 2070–2099 under a high emission scenario (RCP8.5). The large fire season will be more intense and have an earlier onset and delayed end. Our findings suggest that despite the lack of a contemporary trend in fire regime, projected greenhouse gas emissions will substantially increase the fire danger in southern California by 2099.

Dynamical downscaling remains a powerful tool for studying regional climate processes, and the genesis of high-resolution historical and future climate data. This technique is particularly important over areas of complex terrain, such as the western United States (WUS), where global models are especially limited in representing regional climate. After identifying a suite of WRF options that best simulate snow and precipitation for an average water year (2010) over the WUS, we evaluate the performance of the dynamically downscaled European Centre for Medium-range Weather Forecasting's fifth Reanalysis (ERA5) from 1980 to 2020 on 45-km, 9-km, and two 3-km grids. We find that by decreasing the horizontal grid spacing within WRF, improvements to Sierra Nevada and Northern Rocky Mountain snow, Santa Ana and Diablo winds, and coastal meteorology occur. For landfalling atmospheric rivers (ARs), the downscaled reanalysis simulates greater upstream integrated vapor transport (IVT) than ERA5. However, WRF skillfully simulates the positioning of the IVT and the timing and magnitude of AR precipitation. This potential IVT bias, in conjunction with increasing resolution, leads to a wet precipitation bias across the Sierra Nevada in the 3-km experiment. This conclusion is supported by streamflow analysis, although we note that the bias in the 3-km experiment can also be explained by in situ undercatch issues. Meanwhile, the 9-km experiment is more biased than the 3-km experiment across the Northern Rocky Mountains compared to in situ measured SWE and precipitation, indicating a geographic sensitivity to biases.

Quantifying the responses of forest disturbances to climate warming is critical to our understanding of carbon cycles and energy balances of the Earth system. The impact of warming on bark beetle outbreaks is complex as multiple drivers of these events may respond differently to warming. Using a novel model of bark beetle biology and host tree interactions, we assessed how contemporary warming affected western pine beetle (Dendroctonus brevicomis) populations and mortality of its host, ponderosa pine (Pinus ponderosa), during an extreme drought in the Sierra Nevada, California, United States. When compared with the field data, our model captured the western pine beetle flight timing and rates of ponderosa pine mortality observed during the drought. In assessing the influence of temperature on western pine beetles, we found that contemporary warming increased the development rate of the western pine beetle and decreased the overwinter mortality rate of western pine beetle larvae leading to increased population growth during periods of lowered tree defense. We attribute a 29.9% (95% CI: 29.4%–30.2%) increase in ponderosa pine mortality during drought directly to increases in western pine beetle voltinism (i.e., associated with increased development rates of western pine beetle) and, to a much lesser extent, reductions in overwintering mortality. These findings, along with other studies, suggest each degree (°C) increase in temperature may have increased the number of ponderosa pine killed by upwards of 35%–40% °C⁻¹ if the effects of compromised tree defenses (15%–20%) and increased western pine beetle populations (20%) are additive. Due to the warming ability to considerably increase mortality through the mechanism of bark beetle populations, models need to consider climate's influence on both host tree stress and the bark beetle population dynamics when determining future levels of tree mortality.

Previous studies have identified a recent increase in wildfire activity in the western United States (WUS). However, the extent to which this trend is due to weather pattern changes dominated by natural variability versus anthropogenic warming has been unclear. Using an ensemble constructed flow analogue approach, we have employed observations to estimate vapor pressure deficit (VPD), the leading meteorological variable that controls wildfires, associated with different atmospheric circulation patterns. Our results show that for the period 1979 to 2020, variation in the atmospheric circulation explains, on average, only 32% of the observed VPD trend of 0.48 ± 0.25 hPa/decade (95% CI) over the WUS during the warm season (May to September). The remaining 68% of the upward VPD trend is likely due to anthropogenic warming. The ensemble simulations of climate models participating in the sixth phase of the Coupled Model Intercomparison Project suggest that anthropogenic forcing explains an even larger fraction of the observed VPD trend (88%) for the same period and region. These models and observational estimates likely provide a lower and an upper bound on the true impact of anthropogenic warming on the VPD trend over the WUS. During August 2020, when the August Complex “Gigafire” occurred in the WUS, anthropogenic warming likely explains 50% of the unprecedented high VPD anomalies.

Overestimation of precipitation frequency and duration while underestimating intensity, that is, the “drizzling” bias, has been a long-standing problem of global climate models. Here we explore this issue from the perspective of precipitation partitioning. We found that most models in the Climate Model Intercomparison Project Phase 5 (CMIP5) have high convective-to-total precipitation (PC/PR) ratios in low latitudes. Convective precipitation has higher frequency and longer duration but lower intensity than non-convective precipitation in many models. As a result, the high PC/PR ratio contributes to the “drizzling” bias over low latitudes. The PC/PR ratio and associated “drizzling” bias increase as model resolution coarsens from 0.5° to 2.0°, but the resolution's effect weakens as the grid spacing increases from 2.0° to 3.0°. Some of the CMIP6 models show reduced “drizzling” bias associated with decreased PC/PR ratio. Thus, more reasonable precipitation partitioning, along with finer model resolution should alleviate the “drizzling” bias within current climate models.