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.

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.

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.

Biases in climate model simulations introduce biases in subsequent impact simulations. Therefore, bias correction methods are operationally used to post-process regional climate projections. However, many problems have been identified, and some researchers question the very basis of the approach. Here we demonstrate that a typical cross-validation is unable to identify improper use of bias correction. Several examples show the limited ability of bias correction to correct and to downscale variability, and demonstrate that bias correction can cause implausible climate change signals. Bias correction cannot overcome major model errors, and naive application might result in ill-informed adaptation decisions. We conclude with a list of recommendations and suggestions for future research to reduce, post-process, and cope with climate model biases.

Using hybrid dynamical–statistical downscaling, 3-km-resolution end-of-twenty-first-century runoff timing changes over California’s Sierra Nevada for all available global climate models (GCMs) from phase 5 of the Coupled Model Intercomparison Project (CMIP5) are projected. All four representative concentration pathways (RCPs) adopted by the Intergovernmental Panel on Climate Change’s Fifth Assessment Report are examined. These multimodel, multiscenario projections allow for quantification of ensemble-mean runoff timing changes and an associated range of possible outcomes due to both intermodel variability and choice of forcing scenario. Under a “business as usual” forcing scenario (RCP8.5), warming leads to a shift toward much earlier snowmelt-driven surface runoff in 2091–2100 compared to 1991–2000, with advances of as much as 80 days projected in the 35-model ensemble mean. For a realistic “mitigation” scenario (RCP4.5), the ensemble-mean change is smaller but still large (up to 30 days). For all plausible forcing scenarios and all GCMs, the simulated changes are statistically significant, so that a detectable change in runoff timing is inevitable. Even for the mitigation scenario, the ensemble-mean change is approximately equivalent to one standard deviation of the natural variability at most elevations. Thus, even when greenhouse gas emissions are curtailed, the runoff change is climatically significant. For the business-as-usual scenario, the ensemble-mean change is approximately two standard deviations of the natural variability at most elevations, portending a truly dramatic change in surface hydrology by the century’s end if greenhouse gas emissions continue unabated.

High-resolution gridded datasets are in high demand because they are spatially complete and include important finescale details. Previous assessments have been limited to two to three gridded datasets or analyzed the datasets only at the station locations. Here, eight high-resolution gridded temperature datasets are assessed two ways: at the stations, by comparing with Global Historical Climatology Network–Daily data; and away from the stations, using physical principles. This assessment includes six station-based datasets, one interpolated reanalysis, and one dynamically downscaled reanalysis. California is used as a test domain because of its complex terrain and coastlines, features known to differentiate gridded datasets. As expected, climatologies of station-based datasets agree closely with station data. However, away from stations, spread in climatologies can exceed 6°C. Some station-based datasets are very likely biased near the coast and in complex terrain, due to inaccurate lapse rates. Many station-based datasets have large unphysical trends (>1°C decade⁻¹) due to unhomogenized or missing station data—an issue that has been fixed in some datasets by using homogenization algorithms. Meanwhile, reanalysis-based gridded datasets have systematic biases relative to station data. Dynamically downscaled reanalysis has smaller biases than interpolated reanalysis, and has more realistic variability and trends. Dynamical downscaling also captures snow–albedo feedback, which station-based datasets miss. Overall, these results indicate that 1) gridded dataset choice can be a substantial source of uncertainty, and 2) some datasets are better suited for certain applications.

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.

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.

The climate warming effects of accelerated urbanization along with projected global climate change raise an urgent need for sustainable mitigation and adaptation strategies to cool urban climates. Our modeling results show that historical urbanization in the Los Angeles and San Diego metropolitan areas has increased daytime urban air temperature by 1.3 °C, in part due to a weakening of the onshore sea breeze circulation. We find that metropolis-wide adoption of cool roofs can meaningfully offset this daytime warming, reducing temperatures by 0.9 °C relative to a case without cool roofs. Residential cool roofs were responsible for 67% of the cooling. Nocturnal temperature increases of 3.1 °C from urbanization were larger than daytime warming, while nocturnal temperature reductions from cool roofs of 0.5 °C were weaker than corresponding daytime reductions. We further show that cool roof deployment could partially counter the local impacts of global climate change in the Los Angeles metropolitan area. Assuming a scenario in which there are dramatic decreases in greenhouse gas emissions in the 21st century (RCP2.6), mid- and end-of-century temperature increases from global change relative to current climate are similarly reduced by cool roofs from 1.4 °C to 0.6 °C. Assuming a scenario with continued emissions increases throughout the century (RCP8.5), mid-century warming is significantly reduced by cool roofs from 2.0 °C to 1.0 °C. The end-century warming, however, is significantly offset only in small localized areas containing mostly industrial/commercial buildings where cool roofs with the highest albedo are adopted. We conclude that metropolis-wide adoption of cool roofs can play an important role in mitigating the urban heat island effect, and offsetting near-term local warming from global climate change. Global-scale reductions in greenhouse gas emissions are the only way of avoiding long-term warming, however. We further suggest that both climate mitigation and adaptation can be pursued simultaneously using 'cool photovoltaics'.

Future snowfall and snowpack changes over the mountains of Southern California are projected using a new hybrid dynamical–statistical framework. Output from all general circulation models (GCMs) in phase 5 of the Coupled Model Intercomparison Project archive is downscaled to 2-km resolution over the region. Variables pertaining to snow are analyzed for the middle (2041–60) and end (2081–2100) of the twenty-first century under two representative concentration pathway (RCP) scenarios: RCP8.5 (business as usual) and RCP2.6 (mitigation). These four sets of projections are compared with a baseline reconstruction of climate from 1981 to 2000. For both future time slices and scenarios, ensemble-mean total winter snowfall loss is widespread. By the mid-twenty-first century under RCP8.5, ensemble-mean winter snowfall is about 70% of baseline, whereas the corresponding value for RCP2.6 is somewhat higher (about 80% of baseline). By the end of the century, however, the two scenarios diverge significantly. Under RCP8.5, snowfall sees a dramatic further decline; 2081–2100 totals are only about half of baseline totals. Under RCP2.6, only a negligible further reduction from midcentury snowfall totals is seen. Because of the spread in the GCM climate projections, these figures are all associated with large intermodel uncertainty. Snowpack on the ground, as represented by 1 April snow water equivalent is also assessed. Because of enhanced snowmelt, the loss seen in snowpack is generally 50% greater than that seen in winter snowfall. By midcentury under RCP8.5, warming-accelerated spring snowmelt leads to snow-free dates that are about 1–3 weeks earlier than in the baseline period.

The area burned by Southern California wildfires has increased in recent decades, with implications for human health, infrastructure, and ecosystem management. Meteorology and fuel structure are universally recognized controllers of wildfire, but their relative importance, and hence the efficacy of abatement and suppression efforts, remains controversial. Southern California's wildfires can be partitioned by meteorology: fires typically occur either during Santa Ana winds (SA fires) in October through April, or warm and dry periods in June through September (non-SA fires). Previous work has not quantitatively distinguished between these fire regimes when assessing economic impacts or climate change influence. Here we separate five decades of fire perimeters into those coinciding with and without SA winds. The two fire types contributed almost equally to burned area, yet SA fires were responsible for 80% of cumulative 1990–2009 economic losses ($3.1 Billion). The damage disparity was driven by fire characteristics: SA fires spread three times faster, occurred closer to urban areas, and burned into areas with greater housing values. Non-SA fires were comparatively more sensitive to age-dependent fuels, often occurred in higher elevation forests, lasted for extended periods, and accounted for 70% of total suppression costs. An improved distinction of fire type has implications for future projections and management. The area burned in non-SA fires is projected to increase 77% (±43%) by the mid-21st century with warmer and drier summers, and the SA area burned is projected to increase 64% (±76%), underscoring the need to evaluate the allocation and effectiveness of suppression investments.

Regional information on climate change is urgently needed but often deemed unreliable. To achieve credible regional climate projections, it is essential to understand underlying physical processes, reduce model biases and evaluate their impact on projections, and adequately account for internal variability. In the tropics, where atmospheric internal variability is small compared with the forced change, advancing our understanding of the coupling between long-term changes in upper-ocean temperature and the atmospheric circulation will help most to narrow the uncertainty. In the extratropics, relatively large internal variability introduces substantial uncertainty, while exacerbating risks associated with extreme events. Large ensemble simulations are essential to estimate the probabilistic distribution of climate change on regional scales. Regional models inherit atmospheric circulation uncertainty from global models and do not automatically solve the problem of regional climate change. We conclude that the current priority is to understand and reduce uncertainties on scales greater than 100 km to aid assessments at finer scales.

Changes to mean and extreme wet season precipitation over California on interannual time scales are analyzed using twenty-first-century precipitation data from 34 global climate models. Models disagree on the sign of projected changes in mean precipitation, although in most models the change is very small compared to historical and simulated levels of interannual variability. For the 2020/21–2059/60 period, there is no projected increase in the frequency of extremely dry wet seasons in the ensemble mean. Wet extremes are found to increase to around 2 times the historical frequency, which is statistically significant at the 95% level. Stronger signals emerge in the 2060/61–2099/2100 period. Across all models, extremely dry wet seasons are roughly 1.5 to 2 times more common, and wet extremes generally triple in their historical frequency (statistically significant). Large increases in precipitation variability in most models account for the modest increases to dry extremes. Increases in the frequency of wet extremes can be ascribed to equal contributions from increased variability and increases to the mean. These increases in the frequency of interannual precipitation extremes will create severe water management problems in a region where coping with large interannual variability in precipitation is already a challenge. Evidence from models and observations is examined to understand the causes of the low precipitation associated with the 2013/14 drought in California. These lines of evidence all strongly indicate that the low 2013/14 wet season precipitation total can be very likely attributed to natural variability, in spite of the projected future changes in extremes.

Using the hybrid downscaling technique developed in part I of this study, temperature changes relative to a baseline period (1981–2000) in the greater Los Angeles region are downscaled for two future time slices: midcentury (2041–60) and end of century (2081–2100). Two representative concentration pathways (RCPs) are considered, corresponding to greenhouse gas emission reductions over coming decades (RCP2.6) and to continued twenty-first-century emissions increases (RCP8.5). All available global climate models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) are downscaled to provide likelihood and uncertainty estimates. By the end of century under RCP8.5, a distinctly new regional climate state emerges: average temperatures will almost certainly be outside the interannual variability range seen in the baseline. Except for the highest elevations and a narrow swath very near the coast, land locations will likely see 60–90 additional extremely hot days per year, effectively adding a new season of extreme heat. In mountainous areas, a majority of the many baseline days with freezing nighttime temperatures will most likely not occur. According to a similarity metric that measures daily temperature variability and the climate change signal, the RCP8.5 end-of-century climate will most likely be only about 50% similar to the baseline. For midcentury under RCP2.6 and RCP8.5 and end of century under RCP2.6, these same measures also indicate a detectable though less significant climatic shift. Therefore, while measures reducing global emissions would not prevent climate change at this regional scale in the coming decades, their impact would be dramatic by the end of the twenty-first century.

In this study (Part I), the mid-twenty-first-century surface air temperature increase in the entire CMIP5 ensemble is downscaled to very high resolution (2 km) over the Los Angeles region, using a new hybrid dynamical–statistical technique. This technique combines the ability of dynamical downscaling to capture finescale dynamics with the computational savings of a statistical model to downscale multiple GCMs. First, dynamical downscaling is applied to five GCMs. Guided by an understanding of the underlying local dynamics, a simple statistical model is built relating the GCM input and the dynamically downscaled output. This statistical model is used to approximate the warming patterns of the remaining GCMs, as if they had been dynamically downscaled. The full 32-member ensemble allows for robust estimates of the most likely warming and uncertainty resulting from intermodel differences. The warming averaged over the region has an ensemble mean of 2.3°C, with a 95% confidence interval ranging from 1.0° to 3.6°C. Inland and high elevation areas warm more than coastal areas year round, and by as much as 60% in the summer months. A comparison to other common statistical downscaling techniques shows that the hybrid method produces similar regional-mean warming outcomes but demonstrates considerable improvement in capturing the spatial details. Additionally, this hybrid technique incorporates an understanding of the physical mechanisms shaping the region’s warming patterns, enhancing the credibility of the final results.

A new hybrid statistical–dynamical downscaling technique is described to project mid- and end-of-twenty-first-century local precipitation changes associated with 36 global climate models (GCMs) in phase 5 of the Coupled Model Intercomparison Project archive over the greater Los Angeles region. Land-averaged precipitation changes, ensemble-mean changes, and the spread of those changes for both time slices are presented. It is demonstrated that the results are similar to what would be produced if expensive dynamical downscaling techniques were instead applied to all GCMs. Changes in land-averaged ensemble-mean precipitation are near zero for both time slices, reflecting the region’s typical position in the models at the node of oppositely signed large-scale precipitation changes. For both time slices, the intermodel spread of changes is only about 0.2–0.4 times as large as natural interannual variability in the baseline period. A caveat to these conclusions is that interannual variability in the tropical Pacific is generally regarded as a weakness of the GCMs. As a result, there is some chance the GCM responses in the tropical Pacific to a changing climate and associated impacts on Southern California precipitation are not credible. It is subjectively judged that this GCM weakness increases the uncertainty of regional precipitation change, perhaps by as much as 25%. Thus, it cannot be excluded that the possibility that significant regional adaptation challenges related to either a precipitation increase or decrease would arise. However, the most likely downscaled outcome is a small change in local mean precipitation compared to natural variability, with large uncertainty on the sign of the change.

Techniques to downscale global climate model (GCM) output and produce high-resolution climate change projections have emerged over the past two decades. GCM projections of future climate change, with typical resolutions of about 100 km, are now routinely downscaled to resolutions as high as hundreds of meters. Pressure to use these techniques to produce policy-relevant information is enormous. To prevent bad decisions, the climate science community must identify downscaling's strengths and limitations and develop best practices. A starting point for this discussion is to acknowledge that downscaled climate signals arising from warming are more credible than those arising from circulation changes.

Changes in wintertime 10 m winds due to the El Niño-Southern Oscillation are examined using a 6 km resolution climate simulation of Southern California covering the period from 1959 through 2001. Wind speed statistics based on regional averages reveal a general signal of increased mean wind speeds and wind speed variability during El Niño across the region. An opposite and nearly as strong signal of decreased wind speed variability during La Niña is also found. These signals are generally more significant than the better-known signals in precipitation. In spite of these regional-scale generalizations, there are significant sub-regional mesoscale structures in the wind speed impacts. In some cases, impacts on mean winds and wind variability at the sub-regional scale are opposite to those of the region as a whole. All of these signals can be interpreted in terms of shifts in occurrences of the region’s main wind regimes due to the El Niño phenomenon. The results of this study can be used to understand how interannual wind speed variations in regions of Southern California are influenced by the El Niño phenomenon.

Projections of possible precipitation change in California under global warming have been subject to considerable uncertainty because California lies between the region anticipated to undergo increases in precipitation at mid-to-high latitudes and regions of anticipated decrease in the subtropics. Evaluation of the large-scale model experiments for phase 5 of the Coupled Model Intercomparison Project (CMIP5) suggests a greater degree of agreement on the sign of the winter (December–February) precipitation change than in the previous such intercomparison, indicating a greater portion of California falling within the increased precipitation zone. While the resolution of global models should not be relied on for accurate depiction of topographic rainfall distribution within California, the precipitation changes depend substantially on large-scale shifts in the storm tracks arriving at the coast. Significant precipitation increases in the region arriving at the California coast are associated with an eastward extension of the region of strong Pacific jet stream, which appears to be a robust feature of the large-scale simulated changes. This suggests that effects of this jet extension in steering storm tracks toward the California coast constitute an important factor that should be assessed for impacts on incoming storm properties for high-resolution regional model assessments.

The frequency and character of Southern California’s Santa Ana wind events are investigated within a 12-km-resolution downscaling of late-20th and mid-21st century time periods of the National Center for Atmospheric Research Community Climate System Model global climate change scenario run. The number of Santa Ana days per winter season is approximately 20% fewer in the mid 21st century compared to the late 20th century. Since the only systematic and sustained difference between these two periods is the level of anthropogenic forcing, this effect is anthropogenic in origin. In both time periods, Santa Ana winds are partly katabatically-driven by a temperature difference between the cold wintertime air pooling in the desert against coastal mountains and the adjacent warm air over the ocean. However, this katabatic mechanism is significantly weaker during the mid 21st century time period. This occurs because of the well-documented differential warming associated with transient climate change, with more warming in the desert interior than over the ocean. Thus the mechanism responsible for the decrease in Santa Ana frequency originates from a well-known aspect of the climate response to increasing greenhouse gases, but cannot be understood or simulated without mesoscale atmospheric dynamics. In addition to the change in Santa Ana frequency, we investigate changes during Santa Anas in two other meteorological variables known to be relevant to fire weather conditions—relative humidity and temperature. We find a decrease in the relative humidity and an increase in temperature. Both these changes would favor fire. A fire behavior model accounting for changes in wind, temperature, and relative humidity simultaneously is necessary to draw firm conclusions about future fire risk and growth associated with Santa Ana events. While our results are somewhat limited by a relatively small sample size, they illustrate an observed and explainable regional change in climate due to plausible mesoscale processes.

The depth and timing of snowpack in the Sierra Nevada Mountains are of fundamental importance to California water resource availability, and recent studies indicate a shift toward earlier snowmelt consistent with projected impacts of anthropogenic climate change. In order for future studies to assess snowpack variability on seasonal to centennial time scales, physically based models of snowpack evolution at high spatial resolution must be improved. Here we evaluate modeled snowpack accuracy for the central Sierra Nevada in the Weather Research and Forecasting regional climate model coupled to the Noah land surface model. A simulation with nested domains at 27, 9, and 3 km grid spacings is presented for November 2001 to July 2002. Model outputs are compared with daily snowpack observations at 41 locations, air temperature at 31 locations, and precipitation at 10 locations. Comparison of snowpack at different resolutions suggests that 27 km simulations substantially underestimate snowpack, while 9 and 3 km simulations are closer to observations. Regional snowpack accumulation is accurately simulated at these high resolutions, but model snowmelt occurs an average of 22–25 days early. Some error can be traced to differences in elevation and observation scale between point‐based measurements and model grid cells, but these factors cannot explain the persistent bias toward early snowmelt. A high correlation between snowmelt and error in modeled surface air temperature is found, with melt coinciding systematically with excessively cold air temperatures. One possible source of bias is an imbalance in turbulent heat fluxes, erroneously warming the snowpack while cooling the surface atmosphere.

A study of the California Sierra Nevada snowpack has been conducted using snow station observations and reanalysis surface temperature data. Monthly snow water equivalent (SWE) measurements were combined from two datasets to provide sufficient data from 1930 to 2008. The monthly snapshots are used to calculate peak snow mass timing for each snow season. Since 1930, there has been an overall trend toward earlier snow mass peak timing by 0.6 days per decade. The trend toward earlier timing also occurs at nearly all individual stations. Even stations showing an increase in 1 April SWE exhibit the trend toward earlier timing, indicating that enhanced melting is occurring at nearly all stations. Analysis of individual years and stations reveals that warm daily maximum temperatures averaged over March and April are associated with earlier snow mass peak timing for all spatial and temporal scales included in the dataset. The influence is particularly pronounced for low accumulation years indicating the potential importance of albedo feedback for the melting of shallow snow. The robustness of the early spring temperature influence on peak timing suggests the trend toward earlier peak timing is attributable to the simultaneous warming trend (0.1°C decade⁻¹ since 1930, with an acceleration in warming in later time periods). Given future scenarios of warming in California, one can expect acceleration in the trend toward earlier peak timing; this will reduce the warm season storage capacity of the California snowpack.

The Arctic climate is changing rapidly¹. From 1979 to 2006, September sea-ice extent decreased by almost 25% or about 100,000 km² per year (ref. 2). In September 2007, Arctic sea-ice extent reached its lowest level since satellite observations began³and in September 2008, sea-ice cover was still low. This development has raised concerns that the Arctic Ocean could be ice-free in late summer in only a few decades, with important economic and geopolitical implications. Unfortunately, most current climate models underestimate significantly the observed trend in Arctic sea-ice decline⁴, leading to doubts regarding their projections for the timing of ice-free conditions. Here we analyse the simulated trends in past sea-ice cover in 18 state-of-art-climate models and find a direct relationship between the simulated evolution of September sea-ice cover over the twenty-first century and the magnitude of past trends in sea-ice cover. Using this relationship together with observed trends, we project the evolution of September sea-ice cover over the twenty-first century. We find that under a scenario with medium future greenhouse-gas emissions, the Arctic Ocean will probably be ice-free in September before the end of the twenty-first century.

Across vast, agriculturally intensive regions of the United States, the spread in predictions of summer temperature and soil moisture under global warming is curiously elevated in current climate models. Some models show modest warming of 2–3C° and little drying or slight moistening by the 22nd century, while at the other extreme are simulations with warming as large as 7–8C° and 20–40% reductions in soil moisture. We show this region of large spread arises from differences in simulations of snow albedo feedback. During winter and early spring, models with strong snow albedo feedback exhibit large reductions in snowpack and hence water storage. This water deficit persists in summer soil moisture, with reduced evapotranspiration yielding warmer temperatures. Comparison of simulated feedback strength to observations of the feedback from the current climate's seasonal cycle suggests the inter‐model differences are excessive. At the same time, the multi‐model mean feedback strength agrees reasonably well with the observed value. We estimate that if the next generation of models were brought into line with observations of snow albedo feedback, the unusually wide divergence in simulations of summer warming and drying over the US would shrink by roughly one third to one half.

We examine the climatological diurnal cycle of surface air temperature in a 6 km resolution atmospheric simulation of Southern California from 1995 to the present. We find its amplitude and phase both have significant geographical structure. This is most likely due to diurnally-varying flows back and forth across the coastline and elevation isolines resulting from the large daily warming and cooling over land. Because the region’s atmosphere is generally stably stratified, these flow patterns result in air of lower (higher) potential temperature being advected upslope (downslope) during daytime (nighttime). This suppresses the temperature diurnal cycle amplitude at mountaintops where diurnal flows converge (diverge) during the day (night). The nighttime land breeze also advects air of higher potential temperature downslope toward the coast. This raises minimum temperatures in land areas adjacent to the coast in a manner analogous to the daytime suppression of maximum temperature by the cool sea breeze in these same areas. Because stratification is greater in the coastal zone than in the desert interior, these thermal effects of the diurnal winds are not uniform, generating spatial structures in the phase and shape of the temperature diurnal cycle as well as its amplitude. We confirm that the simulated characteristics of the temperature diurnal cycle as well as those of the associated diurnal winds are also found in a network of 30 observation stations in the region. This gives confidence in the simulation’s realism and our study’s findings. Diurnal flows are probably mainly responsible for the geographical structures in the temperature diurnal cycle in other regions of significant topography and surface heterogeneity, their importance depending partly on the degree of atmospheric stratification.

The primary regimes of local atmospheric variability are examined in a 6-km regional atmospheric model of the southern third of California, an area of significant land surface heterogeneity, intense topography, and climate diversity. The model was forced by reanalysis boundary conditions over the period 1995–2003. The region is approximately the same size as a typical grid box of the current generation of general circulation models used for global climate prediction and reanalysis product generation, and so can be thought of as a laboratory for the study of climate at spatial scales smaller than those resolved by global simulations and reanalysis products. It is found that the simulated circulation during the October–March wet season, when variability is most significant, can be understood through an objective classification technique in terms of three wind regimes. The composite surface wind patterns associated with these regimes exhibit significant spatial structure within the model domain, consistent with the complex topography of the region. These regimes also correspond nearly perfectly with the simulation’s highly structured patterns of variability in hydrology and temperature, and therefore are the main contributors to the local climate variability. The regimes are approximately equally likely to occur regardless of the phase of the classical large-scale modes of atmospheric variability prevailing in the Pacific–North American sector. The high degree of spatial structure of the local regimes and their tightly associated climate impacts, as well as their ambiguous relationship with the primary modes of large-scale variability, demonstrate that the local perspective offered by the high-resolution model is necessary to understand and predict the climate variations of the region.

This chapter assesses the scientific literature on projected changes in major climate phenomena and more specifically their relevance for future change in regional climates, contingent on global mean temperatures continue to rise.