We compare historical and end‐of‐century temperature and precipitation patterns over California from one dynamically downscaled simulation using the Weather Research and Forecast (WRF) model and two simulations statistically downscaled using Localized Constructed Analogs (LOCA). We uniquely separate causes of differences between dynamically and statistically based future climate projections into differences in historical climate (gridded observations versus regional climate model output) and differences in how these downscaling techniques explicitly handle future climate changes (numerical modeling versus analogs). In these methods, solutions between different downscaling techniques differ more in the future compared to the historical period. Changes projected by LOCA are insensitive to the choice of driving data. Only through dynamical downscaling can we simulate physically consistent regional springtime warming patterns across the Sierra Nevada, while the statistical simulations inherit an unphysical signal from their parent Global Climate Model (GCM) or gridded data. The results of our study clarify why these different techniques produce different outcomes and may also provide guidance on which downscaled products to use for certain impact analyses in California and perhaps other Mediterranean regimes.

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. This huge increase in runoff during the most extreme AR events could present major flood control challenges for the region.

Precipitation extremes will likely intensify under climate change. However, much uncertainty surrounds intensification of high-magnitude events that are often inadequately resolved by global climate models. In this analysis, we develop a framework involving targeted dynamical downscaling of historical and future extreme precipitation events produced by a large ensemble of a global climate model. This framework is applied to extreme “atmospheric river” storms in California. We find a substantial (10 to 40%) increase in total accumulated precipitation, with the largest relative increases in valleys and mountain lee-side areas. We also report even higher and more spatially uniform increases in hourly maximum precipitation intensity, which exceed Clausius-Clapeyron expectations. Up to 85% of this increase arises from thermodynamically driven increases in water vapor, with a smaller contribution by increased zonal wind strength. These findings imply substantial challenges for water and flood management in California, given future increases in intense atmospheric river-induced precipitation extremes.

Atmospheric rivers (ARs) are characterized by intense moisture transport, which, on landfall, produce precipitation which can be both beneficial and destructive. ARs in California, for example, are known to have ended drought conditions but also to have caused substantial socio-economic damage from landslides and flooding linked to extreme precipitation. Understanding how AR characteristics will respond to a warming climate is, therefore, vital to the resilience of communities affected by them, such as the western USA, Europe, East Asia and South Africa. In this Review, we use a theoretical framework to synthesize understanding of the dynamic and thermodynamic responses of ARs to anthropogenic warming and connect them to observed and projected changes and impacts revealed by observations and complex models. Evidence suggests that increased atmospheric moisture (governed by Clausius–Clapeyron scaling) will enhance the intensity of AR-related precipitation — and related hydrological extremes — but with changes that are ultimately linked to topographic barriers. However, due to their dependency on both weather and climate-scale processes, which themselves are often poorly constrained, projections are uncertain. To build confidence and improve resilience, future work must focus efforts on characterizing the multiscale development of ARs and in obtaining observations from understudied regions, including the West Pacific, South Pacific and South Atlantic.

Atmospheric rivers (ARs) are responsible for a majority of extreme precipitation and flood events along the U.S. West Coast. To better understand the present‐day characteristics of AR‐related precipitation extremes, a selection of nine most intense historical AR events during 1980–2017 is simulated using a dynamical downscaling modeling framework based on the Weather Research and Forecasting Model. We find that the chosen framework and Weather Research and Forecasting Model configuration reproduces both large‐scale atmospheric features—including parent synoptic‐scale cyclones—as well as the filamentary corridors of integrated vapor transport associated with the ARs themselves. The accuracy of simulated extreme precipitation maxima, relative to in situ and interpolated gridded observations, improves notably with increasing model resolution, with improvements as large as 40–60% for fine scale (3 km) relative to coarse‐scale (27 km) simulations. A separate set of simulations using smoothed topography suggests that much of these gains stem from the improved representation of complex terrain. Additionally, using the 12 December 1995 storm in Northern California as an example, we demonstrate that only the highest‐resolution simulations resolve important fine‐scale features—such as localized orographically forced vertical motion and powerful near hurricane‐force boundary layer winds. Given the demonstrated ability of a targeted dynamical downscaling framework to capture both local extreme precipitation and key fine‐scale characteristics of the most intense ARs in the historical record, we argue that such a configuration may be highly conducive to understanding AR‐related extremes and associated changes in a warming climate.

Arctic sea ice has decreased substantially over recent decades, a trend projected to continue. Shrinking ice reduces surface albedo, leading to greater surface solar absorption, thus amplifying warming and driving further melt. This sea-ice albedo feedback (SIAF) is a key driver of Arctic climate change and an important uncertainty source in climate model projections. Using an ensemble of models, we demonstrate an emergent relationship between future SIAF and an observable version of SIAF in the current climate’s seasonal cycle. This relationship is robust in constraining SIAF over the coming decades (Pearson’s r = 0.76), and then it degrades. The degradation occurs because some models begin producing ice-free conditions, signalling a transition to a new ice regime. The relationship is strengthened when models with unrealistically thin historical ice are excluded. Because of this tight relationship, reducing model errors in the current climate’s seasonal SIAF and ice thickness can narrow SIAF spread under climate change.

Purpose of Review

Highlight significant developments that have recently been made to enhance our understanding of how snow responds to climate forcing and the role that snow plays in the climate system.

Recent Findings

Widespread snow loss has occurred in recent decades, with the largest decreases in spring. These changes are primarily driven by temperature and precipitation, but changes in vegetation, light-absorbing impurities, and sea ice also contribute to variability. Changes in snow cover can also affect climate through the snow albedo feedback (SAF). Recently, considerable progress has been made in better understanding the processes contributing to SAF. We also highlight advances in knowledge of how snow variability is linked to large-scale atmospheric changes. Lastly, large-scale snow losses are expected to continue under climate change in all but the coldest climates. These projected changes to snow raise considerable concerns over future freshwater availability in snow-dominated watersheds.

Summary

The results discussed here demonstrate the widespread implications that changes to snow have on the climate system and anthropogenic activity at large.

Earth system models (ESMs) are key tools for providing climate projections under different scenarios of human-induced forcing. ESMs include a large number of additional processes and feedbacks such as biogeochemical cycles that traditional physical climate models do not consider. Yet, some processes such as cloud dynamics and ecosystem functional response still have fairly high uncertainties. In this article, we present an overview of climate feedbacks for Earth system components currently included in state-of-the-art ESMs and discuss the challenges to evaluate and quantify them. Uncertainties in feedback quantification arise from the interdependencies of biogeochemical matter fluxes and physical properties, the spatial and temporal heterogeneity of processes, and the lack of long-term continuous observational data to constrain them. We present an outlook for promising approaches that can help to quantify and to constrain the large number of feedbacks in ESMs in the future. The target group for this article includes generalists with a background in natural sciences and an interest in climate change as well as experts working in interdisciplinary climate research (researchers, lecturers, and students). This study updates and significantly expands upon the last comprehensive overview of climate feedbacks in ESMs, which was produced 15 years ago (NRC, 2003).

This study uses dynamical and statistical methods to understand end‐of‐century mean changes to Sierra Nevada snowpack. Dynamical results reveal mid‐elevation watersheds experience considerably more rain than snow during winter, leading to substantial snowpack declines by spring. Despite some high‐elevation watersheds receiving slightly more snow in January and February, the warming signal still dominates across the wet‐season and leads to notable declines by springtime. A statistical model is created to mimic dynamical results for April 1 snowpack, allowing for an efficient downscaling of all available General Circulation Models (GCMs) from the Coupled Model Intercomparison Project Phase 5. For all GCMs and emissions scenarios, dramatic April 1 snowpack loss occurs at elevations below 2500 meters, despite increased precipitation in many GCMs. Only 36% (±12%) of historical April 1 total snow water equivalent volume remains at the century's end under a “business‐as‐usual” emissions scenario, with 70% (±12%) remaining under a realistic “mitigation” scenario.

In recent years, an evaluation technique for Earth System Models (ESMs) has arisen—emergent constraints (ECs)—which rely on strong statistical relationships between aspects of current climate and future change across an ESM ensemble. Combining the EC relationship with observations could reduce uncertainty surrounding future change. Here, we articulate a framework to assess ECs, and provide indicators whereby a proposed EC may move from a strong statistical relationship to confirmation. The primary indicators are verified mechanisms and out-of-sample testing. Confirmed ECs have the potential to improve ESMs by focusing attention on the variables most relevant to climate projections. Looking forward, there may be undiscovered ECs for extremes and teleconnections, and ECs may help identify climate system tipping points.

Differences among climate models in equilibrium climate sensitivity (ECS; the equilibrium surface temperature response to a doubling of atmospheric CO₂) remain a significant barrier to the accurate assessment of societally important impacts of climate change. Relationships between ECS and observable metrics of the current climate in model ensembles, so-called emergent constraints, have been used to constrain ECS. Here a statistical method (including a backward selection process) is employed to achieve a better statistical understanding of the connections between four recently proposed emergent constraint metrics and individual feedbacks influencing ECS. The relationship between each metric and ECS is largely attributable to a statistical connection with shortwave low cloud feedback, the leading cause of intermodel ECS spread. This result bolsters confidence in some of the metrics, which had assumed such a connection in the first place. Additional analysis is conducted with a few thousand artificial metrics that are randomly generated but are well correlated with ECS. The relationships between the contrived metrics and ECS can also be linked statistically to shortwave cloud feedback. Thus, any proposed or forthcoming ECS constraint based on the current generation of climate models should be viewed as a potential constraint on shortwave cloud feedback, and physical links with that feedback should be investigated to verify that the constraint is real. In addition, any proposed ECS constraint should not be taken at face value since other factors influencing ECS besides shortwave cloud feedback could be systematically biased in the models.

Mediterranean climate regimes are particularly susceptible to rapid shifts between drought and flood—of which, California’s rapid transition from record multi-year dryness between 2012 and 2016 to extreme wetness during the 2016–2017 winter provides a dramatic example. Projected future changes in such dry-to-wet events, however, remain inadequately quantified, which we investigate here using the Community Earth System Model Large Ensemble of climate model simulations. Anthropogenic forcing is found to yield large twenty-first-century increases in the frequency of wet extremes, including a more than threefold increase in sub-seasonal events comparable to California’s ‘Great Flood of 1862’. Smaller but statistically robust increases in dry extremes are also apparent. As a consequence, a 25% to 100% increase in extreme dry-to-wet precipitation events is projected, despite only modest changes in mean precipitation. Such hydrological cycle intensification would seriously challenge California’s existing water storage, conveyance and flood control infrastructure.

This study investigates temperature impacts to snowpack and runoff‐driven flood risk over the Sierra Nevada during the extremely wet year of 2016–2017, which followed the extraordinary California drought of 2011–2015. By perturbing near‐surface temperatures from a 9‐km dynamically downscaled simulation, a series of offline land surface model experiments explore how Sierra Nevada hydrology has already been impacted by historical anthropogenic warming and how these impacts evolve under future warming scenarios. Results show that historical warming reduced 2016–2017 Sierra Nevada snow water equivalent by 20% while increasing early‐season runoff by 30%. An additional one third to two thirds loss of snowpack is projected by the end of the century, depending on the emission scenario, with middle elevations experiencing the most significant declines. Notably, the number of days in the future with runoff exceeding 20 mm nearly doubles under a mitigation emission scenarios and triples under a business‐as‐usual scenario. A smaller snow‐to‐rain ratio, as opposed to increased snowmelt, is found to be the primary mechanism of temperature impacts to Sierra snowpack and runoff. These findings are consequential to the prevalence of early‐season floods in the Sierra Nevada. In the Feather River Watershed, historical warming increased runoff by over one third during the period of heaviest precipitation in February 2017. This suggests that historical anthropogenic warming may have exacerbated runoff conditions underlying the Oroville Dam spillway overflow that occurred in this month. As warming continues in the future, the potential for runoff‐based flood risk may rise even higher.

Snow albedo feedback (SAF) behaves similarly in the current and future climate contexts; thus, constraining the large intermodel variance in SAF will likely reduce uncertainty in climate projections. To better understand this intermodel spread, structural and parametric biases contributing to SAF variability are investigated. We find that structurally varying snowpack, vegetation, and albedo parameterizations drive most of the spread, while differences arising from model parameters are generally smaller. Models with the largest SAF biases exhibit clear structural or parametric errors. Additionally, despite widespread intermodel similarities, model interdependency has little impact on the strength of the relationship between SAF in the current and future climate contexts. Furthermore, many models now feature a more realistic SAF than in the prior generation, but shortcomings from two models limit the reduction in ensemble spread. Lastly, preliminary signs from ongoing model development are positive and suggest a likely reduction in SAF spread among upcoming models.

A highly uncertain aspect of anthropogenic climate change is the rate at which the global hydrologic cycle intensifies. The future change in global‐mean precipitation per degree warming, or hydrologic sensitivity, exhibits a threefold spread (1–3%/K) in current global climate models. In this study, we find that the intermodel spread in this value is associated with a significant portion of variability in future projections of extreme precipitation in the tropics, extending also into subtropical atmospheric river corridors. Additionally, there is a very tight intermodel relationship between changes in extreme and nonextreme precipitation, whereby models compensate for increasing extreme precipitation events by decreasing weak‐moderate events. Another factor linked to changes in precipitation extremes is model resolution, with higher resolution models showing a larger increase in heavy extremes. These results highlight ways various aspects of hydrologic cycle intensification are linked in models and shed new light on the task of constraining precipitation extremes.

Emergent constraints use relationships between future and current climate states to constrain projections of climate response. Here we introduce a statistical, hierarchical emergent constraint (HEC) framework in order to link future and current climates with observations. Under Gaussian assumptions, the mean and variance of the future state are shown analytically to be a function of the signal‐to‐noise ratio between current climate uncertainty and observation error and the correlation between future and current climate states. We apply the HEC to the climate change, snow‐albedo feedback, which is related to the seasonal cycle in the Northern Hemisphere. We obtain a snow‐albedo feedback prediction interval of (−1.25,−0.58)%/K. The critical dependence on signal‐to‐noise ratio and correlation shows that neglecting these terms can lead to bias and underestimated uncertainty in constrained projections. The flexibility of using HEC under general assumptions throughout the Earth system is discussed.

This paper describes ESM-SnowMIP, an international coordinated modelling effort to evaluate current snow schemes, including snow schemes that are included in Earth system models, in a wide variety of settings against local and global observations. The project aims to identify crucial processes and characteristics that need to be improved in snow models in the context of local- and global-scale modelling. A further objective of ESM-SnowMIP is to better quantify snow-related feedbacks in the Earth system. Although it is not part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6), ESM-SnowMIP is tightly linked to the CMIP6-endorsed Land Surface, Snow and Soil Moisture Model Intercomparison (LS3MIP).

California’s Sierra Nevada is a high-elevation mountain range with significant seasonal snow cover. Under anthropogenic climate change, amplification of the warming is expected to occur at elevations near snow margins due to snow albedo feedback. However, climate change projections for the Sierra Nevada made by global climate models (GCMs) and statistical downscaling methods miss this key process. Dynamical downscaling simulates the additional warming due to snow albedo feedback. Ideally, dynamical downscaling would be applied to a large ensemble of 30 or more GCMs to project ensemble-mean outcomes and intermodel spread, but this is far too computationally expensive. To approximate the results that would occur if the entire GCM ensemble were dynamically downscaled, a hybrid dynamical–statistical downscaling approach is used. First, dynamical downscaling is used to reconstruct the historical climate of the 1981–2000 period and then to project the future climate of the 2081–2100 period based on climate changes from five GCMs. Next, a statistical model is built to emulate the dynamically downscaled warming and snow cover changes for any GCM. This statistical model is used to produce warming and snow cover loss projections for all available CMIP5 GCMs. These projections incorporate snow albedo feedback, so they capture the local warming enhancement (up to 3°C) from snow cover loss that other statistical methods miss. Capturing these details may be important for accurately projecting impacts on surface hydrology, water resources, and ecosystems.

Sierra Nevada climate and snowpack is simulated during the period of extreme drought from 2011 to 2015 and compared to an identical simulation except for the removal of the twentieth century anthropogenic warming. Anthropogenic warming reduced average snowpack levels by 25%, with middle‐to‐low elevations experiencing reductions between 26 and 43%. In terms of event frequency, return periods associated with anomalies in 4 year 1 April snow water equivalent are estimated to have doubled, and possibly quadrupled, due to past warming. We also estimate effects of future anthropogenic warmth on snowpack during a drought similar to that of 2011–2015. Further snowpack declines of 60–85% are expected, depending on emissions scenario. The return periods associated with future snowpack levels are estimated to range from millennia to much longer. Therefore, past human emissions of greenhouse gases are already negatively impacting statewide water resources during drought, and much more severe impacts are likely to be inevitable.