In part III of a three-part study on North American climate in phase 5 of the Coupled Model Intercomparison Project (CMIP5) models, the authors examine projections of twenty-first-century climate in the representative concentration pathway 8.5 (RCP8.5) emission experiments. This paper summarizes and synthesizes results from several coordinated studies by the authors. Aspects of North American climate change that are examined include changes in continental-scale temperature and the hydrologic cycle, extremes events, and storm tracks, as well as regional manifestations of these climate variables. The authors also examine changes in the eastern North Pacific and North Atlantic tropical cyclone activity and North American intraseasonal to decadal variability, including changes in teleconnections to other regions of the globe. Projected changes are generally consistent with those previously published for CMIP3, although CMIP5 model projections differ importantly from those of CMIP3 in some aspects, including CMIP5 model agreement on increased central California precipitation. The paper also highlights uncertainties and limitations based on current results as priorities for further research. Although many projected changes in North American climate are consistent across CMIP5 models, substantial intermodel disagreement exists in other aspects. Areas of disagreement include projections of changes in snow water equivalent on a regional basis, summer Arctic sea ice extent, the magnitude and sign of regional precipitation changes, extreme heat events across the northern United States, and Atlantic and east Pacific tropical cyclone activity.
The twenty-first-century projections generated by phase 5 of the Coupled Model Intercomparison Project (CMIP5; Taylor et al. 2012) are analyzed here to assess climate change in North America (NA). This study accompanies two companion papers (Sheffield et al. 2013a, hereafter Part I; Sheffield et al. 2013b, hereafter Part II) that assess the CMIP5 models’ potential to accurately simulate regional climate in the twentieth century. Additionally, it provides an overview and is a first step toward integrating the understanding of climate projection results from the individual papers in the Journal of Climate special collection entitled “North American climate in CMIP5 experiments.” This paper first examines the changes in the continent-wide distribution of seasonal precipitation and temperature in simulations making use of representative concentration pathway 8.5 (RCP8.5; Taylor et al. 2012). It then focuses on a select set of regional climate features. These changes are considered in the context of the ability of models to accurately simulate current climate, discussed in the two companion papers (Part I and Part II), which is generally comparable to that of CMIP3 models, with some improvement noted for individual models.
Previous projections of NA climate change (e.g., CMIP3) have been evaluated as part of earlier climate assessments (Solomon et al. 2007). The CMIP3 consensus projection indicated that, by 2080–99, annual mean temperature increases are very likely across NA with the greatest changes in northern Canada and Alaska, where 10°C mean wintertime temperature increases are projected in some scenarios (Solomon et al. 2007). Increases in annual mean precipitation are projected for the northern tier of the United States, northward into Canada, with projected decreases for the southwest United States, east Pacific warm pool, Caribbean, and adjacent land areas (e.g., Neelin et al. 2006; Seager et al. 2007; Seager and Vecchi 2010).
Beyond mean state changes, CMIP3 models predict a general increase in precipitation intensity (e.g., Diffenbaugh et al. 2005; Mahajan et al. 2012), particularly in the northern tier of the United States and Canada (Tebaldi et al. 2006). Increases in the duration and severity of drought are projected in regions such as Central America and midlatitude NA (e.g., Sheffield and Wood 2008), of which increased temperatures and evapotranspiration are major components (Francina et al. 2010; Gutzler and Robbins 2011; Wehner et al. 2011). A general increase in heat waves, decrease in cold extremes, decrease in frost days, and increase in length of the growing season have been projected across large portions of NA (Meehl and Tebaldi 2004; Diffenbaugh et al. 2005; Biasutti et al. 2012; Christiansen et al. 2011; Diffenbaugh and Scherer 2011; Duffy and Tebaldi 2012; Lau and Nath 2012), projected trends that are generally consistent with observed trends in such quantities over the last century (Alexander et al. 2006). Decreases in the duration of the snowpack have been projected for many regions, in particular low altitude areas of the Pacific Northwest and Rockies (e.g., Brown and Mote 2009; Elsner et al. 2010). Such changes are likely to lead to earlier spring snowmelt in many areas of the west (e.g., Hay et al. 2011). While model agreement is good on projected overall snow water equivalent declines in many areas by the end of the twenty-first century, some models show increases in snowpack along the Arctic Rim by 2100 (e.g., Brown and Mote 2009), particularly at the height of the winter season, even though the length of the snow season shortens (e.g., Räisänen 2008).
The projected response of NA climate in future emission scenarios is often more nuanced on the regional and local scales than for the continental-scale features, especially when considering the evolution during the seasonal cycle. For example, Rauscher et al. (2008) noted an earlier onset of the midsummer drought in Mexico and Central America in model projections for the end of the twenty-first century. Previous studies project a redistribution of precipitation in monsoon regimes such as the southwest United States with reduced spring rainfall and increased late rainy season rainfall (Seth et al. 2010, 2011; Biasutti and Sobel 2009). Ruiz-Barradas and Nigam (2010) showed that projections for the twenty-first century indicate a wetter north-central United States during spring (increase in number of extreme springs) and a drier southwest United States but little consistency in summer rainfall tendencies among models in these same regions. The uncertainty in projected summer precipitation extends to adjacent land areas of the Gulf of Mexico (Biasutti et al. 2012). Studies using CMIP3 projections suggest that, while the total number of North Atlantic tropical cyclones (TCs) will decrease and the number of intense hurricanes will increase, changes in North Atlantic TC activity remain uncertain. This is likely because climate models produce differing patterns of tropical SST change and different representations of tropical Atlantic SST relative to the tropical mean SST, which has been suggested to be a strong regulator of Atlantic TC activity (e.g., Latif et al. 2007; Swanson 2008; Vecchi et al. 2008; Wang and Lee 2008; Knutson et al. 2010; Vecchi et al. 2011; Zhao and Held 2012). Past studies using CMIP3-class models have generally indicated that climate projections for the twenty-first century at the local and regional levels remain a substantial challenge.
The present study provides a summary of projected twenty-first-century NA climate change in the updated state-of-the-art climate and Earth system models used in CMIP5. The results contained herein are contributed by members of the CMIP5 Task Force of the National Oceanographic and Atmospheric Administration (NOAA) Modeling, Analysis, Predictions and Projections Program (MAPP). Where appropriate, we make reference to individual papers submitted in parallel with this comprehensive study to the Journal of Climate special collection entitled “North American climate in CMIP5 experiments.” These individual contributions provide further depth to and physical interpretation of the findings summarized here. The current paper is one of three papers (with Part I and Part II) that synthesize the results and form the core of the special collection, and they represent an initial step toward integrating our understanding of CMIP5 evaluations and projections for North America. We largely focus on RCP8.5 in a core set of 17 CMIP5 models.
Section 2 provides a brief introduction to CMIP5 as well as the primary climate change experiment (RCP8.5). Section 3 presents an assessment of continental climate changes over the twenty-first century, and section 4 assesses regional climate changes. How intraseasonal variability will change in the twenty-first century is assessed in section 5. Changes in Atlantic and east Pacific TC activity are examined in section 6. Multidecadal trends in interannual to decadal hydroclimate variability are analyzed in section 7. Conclusions and a discussion are presented in section 8.
2. CMIP5 models and experiments
We use CMIP5 multimodel datasets of historical climate and climate change experiments (Taylor et al. 2012). These are long-term century-scale projections of climate based on coupled simulations that include a representation of future atmospheric composition from the RCPs (Meinshausen et al. 2011). Table 1 summarizes information on the models used in this study. As noted in Taylor et al. (2012), in addition to physical improvements made in many models, one advantage provided by the CMIP5 experiments versus the CMIP3 effort is that the horizontal resolution of the atmospheric components of the coupled models has significantly increased. About one-third of the models have atmospheric resolution of approximately 1.5° latitude or less, an improvement over CMIP3 where only about 10% of models met this criterion. This higher resolution is of some help in discerning the regional structure of hydroclimate variables over NA. However, in regions of complex topography and coastlines, the resolution of CMIP5 models remains insufficient for resolving important dynamic and thermodynamic features. Where appropriate in the text, we provide contrasts between the current CMIP5 results and those previously derived from CMIP3.
Results based on RCP8.5 will be highlighted here, as it represents one of the core concentration pathways used for the CMIP5 project (Taylor et al. 2012). This experiment represents a high concentration pathway in which radiative forcing due to anthropogenic factors reaches 8.5 W m−2 by 2100 (e.g., Meinshausen et al. 2011, Fig. 4) and continues to grow thereafter. Selected analyses also provide a comparison to a more moderate mitigation pathway (RCP4.5) in which stabilization at 4.5 W m−2 occurs around 2050, and then forcing remains fixed. In terms of the time evolution and value of globally averaged radiative forcing at 2100, RCP8.5 and RCP4.5 most closely resemble the A2 and B1 scenarios for CMIP3 used in the International Panel on Climate Change (IPCC) Fourth Assessment Report (Solomon et al. 2007, Fig. 10.26), respectively. The projection experiments are compared to historical runs of the same models forced by observed trace gases, natural and anthropogenic aerosols, solar forcing, and other agents from the mid-nineteenth century onward (Taylor et al. 2012). A more comprehensive analysis of model performance in the historical runs is provided in the companion papers (Part I and Part II), which provide an additional baseline for comparison with the results shown here. No downscaling or bias correction is used before presentation of results. The exception to lack of downscaling is contained within the supplementary material, where a high-resolution model is used to assess future changes in tropical cyclone activity. Further, the use of downscaling in past and potential future studies is referenced at certain other points in the manuscript.
Multimodel ensemble mean (MEM) differences are highlighted for most of the analyses, as the MEM produces demonstrably superior results in historical climate assessment to those from an individual model (e.g., Gleckler et al. 2008; Pierce et al. 2009). We also use intermodel variability about the MEM to assess model consensus, including the likelihood of specific climate changes, with aspects demonstrating lack of model consensus summarized in the conclusions. In places, the methodology used is more diverse than a simple MEM analysis, given that this paper represents a synthesis of ideas from individual papers in the special collection. Ideally, we would like to compare common future and historical base periods among analyses. Unfortunately, from a practical standpoint this was not always possible. In the analysis of historical simulations in Part I and Part II, the base periods were often determined by the availability of the observational data to assess the models, which were application specific. Further, we assess projected changes in phenomena that have different time scales that range from synoptic to decadal. For example, assessment of interannual to decadal variability requires a longer record than assessment of tropical intraseasonal variability to assess statistical significance.
For consistency with Part I and Part II, our analysis concentrates on the core set of CMIP5 model highlighted by asterisks in Table 1. Part I discussed the selection criteria for these models, which meet the need to include contributions from a large and diverse set of modeling groups and model types. The number of models used in a particular analysis shown below is often limited by availability of data at the time of this study or local storage space, although we try to be as comprehensive as possible. For example, for many of the analyses requiring high-resolution data, a smaller subset of models was used because of the lack of RCP8.5 data availability. Further, because of the large number of contributors from different institutions, overall coordination and unified model selection were not always possible. For some analyses, the number of models used was significantly lower than that of the core set, or the RCP4.5 scenario was used rather than RCP8.5. In these cases, while the results are still enlightening, we have placed the details of these analyses into the supplementary material. These include analyses of moisture transport and diurnal temperature range changes, as well as an analysis of tropical cyclone activity change using a downscaling technique with a high-resolution model. We also occasionally include a more expansive set than the core models in an individual analysis, although we comment on how results would differ if a smaller subset including only core models were used.
3. Continental climate
a. Temperature and precipitation
We first examine projected changes on the continental scale at the end of the twenty-first century relative to the twentieth century climate. Part I noted that CMIP5 models have success in capturing the broad-scale features of NA surface temperature and precipitation in current climate, although with some regional-scale biases. Figure 1 shows the 17-member MEM December–February (DJF) and June–August (JJA) precipitation changes during 2070–99 for RCP8.5 relative to a 1961–90 base period. For models that have more than one run with the same forcing, the average is taken over all runs prior to forming the MEM. A two-sided t test comparing the MEM change to a standard error associated with interannual variability is shown at the 95% confidence level. Note that this tests only the sampling error associated with interannual variability in forming the MEM. Figure 2 shows the model agreement for the precipitation changes, along with two additional criteria for evaluation of significance that are described in the figure caption.
Figures 1 and 2 indicate increases in MEM wintertime precipitation along the west coast of NA from California to Alaska, as well as along the NA east coast from the Mid-Atlantic states northward. Model agreement for these changes is high north of about 40°N, where all but one or two models agree on the sign of changes for all locations. Comparison to a similar ensemble of 16 CMIP3 models indicates that, while the large-scale pattern of precipitation increases at middle to high latitudes and precipitation decreases in the subtropics are similar between the two intercomparisons, one notable difference is that the boundary between these changes has shifted slightly south. This yields projected precipitation increases over parts of California in CMIP5, passing the binomial test for agreement on sign at levels exceeding 95% at points seen in Fig. 2. High interannual variance over coastal land points prevents these from passing the stricter Neelin et al. (2006) criterion. Area averages pass significance tests on the model ensemble (Neelin et al. 2013), which points out a relationship between the extension of storm-track-associated precipitation in this region and the regional manifestation of jet stream increases at the steering level. For such differences at the boundaries of precipitation features, it remains an open question whether the CMIP5 ensemble should be given any additional weight in assessment relative to the CMIP3 ensemble.
Summertime MEM precipitation changes are characterized by higher precipitation amounts in Alaska and the Yukon, where the models all agree on the sign of the changes. All models also suggest precipitation increases along the Arctic coast across the entire length of NA. The MEM indicates reduced summertime precipitation in the east Pacific warm pool and the Caribbean, with agreement of all models in the vicinity of major Caribbean islands, the Yucatan, and in southwestern Mexico adjacent to the east Pacific warm pool. The agreement on these changes for the Caribbean and Mexico was high for CMIP3 models (e.g., Neelin et al. 2006) and is reinforced as a region of even higher intermodel agreement in CMIP5. Because of the model disagreement in projections of future tropical cyclone activity for the Atlantic and east Pacific (shown in section 6 below), it is unlikely that changes in tropical cyclone activity are responsible for these precipitation decreases given the strong model agreement in the precipitation change. Figure S1 (in the supplemental material) provides maps of the MEM percentage precipitation change and its corresponding multimodel standard deviation for the core models used here.
Figure 3 shows the comparable MEM changes for surface air temperature (2-m level) during JJA and DJF. As expected, warming is projected across all regions of NA, with the greatest warming concentrated during wintertime at high latitude regions, where the MEM temperature increase peaks near 15°C in the vicinity of Hudson Bay. Land regions warm more than ocean regions, associated in part with ocean heat storage causing a lag relative to the ongoing greenhouse gas increase, as in the CMIP3-based assessments (Meehl et al. 2007). Over the lower 48 United States and most of Canada, the MEM warming exceeds 5°C in JJA, with an intermodel standard deviation slightly over 1°C in the United States increasing to 2°C in northern Canada. In DJF, both the ensemble mean warming and the intermodel standard deviation have a strong poleward gradient, with warming around 4°C (standard deviation of 1°C) in the southern United States, increasing to over 12°C (standard deviation of over 3°C) in far northern Canada.
While we cannot provide extensive detail on all aspects of projected regional precipitation and temperature change here, Tables S1–S10 (in the supplemental material) provide the MEM and multimodel standard deviation of precipitation and temperature changes as a function of season for the regions defined in Fig. 4. Comparing the individual models or the standard deviation of the model ensemble provides additional information regarding how much confidence should be placed in the MEM mean for a particular region. For instance, in central North America in JJA (Table S9), the intermodel standard deviation is roughly double the MEM value, indicating low confidence that this should be interpreted as significantly different from zero. This is consistent with the summary statistics in Fig. 2, in which only a small portion of this region passes either the binomial test for agreement on sign or the Neelin et al. (2006) criterion that includes a model-by-model t test for significance with respect to interannual variability.
b. Evapotranspiration and runoff
Future changes in precipitation and how it is partitioned into evaporation and runoff have implications for water availability and the occurrence of extreme hydrological events such as floods and droughts. CMIP3 projections indicated that changes in precipitation coupled with increased potential for evaporation from higher atmospheric moisture holding capacity leads to the drying of subtropical regions including Central America and the southwestern United States (Wang 2005; Seager et al. 2007; Sheffield and Wood 2008) and an increased potential for flooding because of more frequent and intense precipitation events, changes in snow accumulation and melt timing, and changes in antecedent soil moisture conditions (e.g., Hamlet and Lettenmaier 2007; Das et al. 2011). Here we analyze changes in the terrestrial water budget on a regional basis by calculating 30-yr averages in annual mean precipitation, evapotranspiration, and runoff between 1971–2000 and 2071–2100 for 15 core models for six regions of NA (Fig. 4). Changes in water storage (i.e., surface water, soil moisture, groundwater, etc.) over the 30-yr periods are assumed to be small compared to the other terms. Thus, precipitation should equal the sum of evapotranspiration and runoff.
For western and eastern NA regions, an increase in MEM annual precipitation is projected, consistent with Fig. 1. The precipitation increases in these regions are apportioned more to evapotranspiration than to runoff (which is important for understanding changes in the availability of water), although models tend to overestimate evapotranspiration in historical simulations (Part I). In the central region, annual mean precipitation increases are more modest. In high latitudes (Alaska–northwest Canada and northeast Canada) the MEM precipitation is projected to increase, consistent with Fig. 1, and is mostly partitioned into increases in runoff, rather than increases in evapotranspiration. This is likely because 1) higher temperatures will increase the proportion of rainfall to snowfall and will melt the snowpack earlier and faster in the spring and 2) the increased precipitation will come in more intense events. In Central America, precipitation is projected to decrease with most of the decrease manifesting in decreasing runoff.
Reductions in snow cover extent and amount are expected in the future as a result of increasing temperatures modified by changes in precipitation and their seasonal interactions. This has important implications for water supply, hydropower generation, and ecosystems and feedbacks with the underlying soil and permafrost (Lawrence and Slater 2010) and to the climate system through changes in albedo (Qu and Hall 2007). CMIP3 projections (Räisänen 2008) indicated that warming reduces the snow season length from both ends across NA, but midwinter snow water equivalent (SWE) is expected to increase in high latitude colder regions because of increased winter snowfall but decrease in regions to the south, where temperature effects on precipitation phase and melting dominate any changes in precipitation amount.
Changes in snow are calculated using the CMIP5 model SWE values. Figure 5 (top) shows the seasonal cycle of changes between 2070–99 and 1971–2000 in monthly mean SWE averaged over NA for 15 core models. All models project a decrease in SWE throughout the year with maximum changes during the peak of the snow season in January–April. The MEM decrease in SWE averaged over NA is about −30 mm (with models ranging from −80 to −10 mm) in the spring and about −10 mm in the summer (ranging from 0 to −65 mm). Spatially, the majority of NA (south of 70°N) is projected to experience a decline in SWE where increasing temperatures have a dominant effect by reducing the ratio of snowfall to rainfall and accelerating melting (Fig. 5, bottom). These reductions are concentrated in the Rocky Mountains to southern Alaska, in the eastern provinces of Canada, and with lower magnitude in the Canadian Prairies. North of 70°N, SWE is projected to increase in places because of increasing precipitation, which outweighs the effects of increasing temperature. Uncertainties across models are likely associated with differences in the temperature projections, to which modeled snow is highly sensitive (Räisänen 2008). At higher latitudes, the sign of the change is uncertain in transitional regions because of the competing effects of increasing snowfall and warming temperatures and in regions of increasing SWE where the magnitude of the precipitation increase is also quite uncertain.
d. Growing season length
Projected warming will likely impact temperature-sensitive flora and fauna, as well as agriculture. We calculate changes in biophysically relevant temperature thresholds including the last spring freeze, the first autumn freeze, and the growing season length at 2071–2100. The growing season length is defined as the number of days between the last spring freeze and the first autumn freeze in the same year. A hard freeze occurs when the daily maximum temperature drops below −2°C. (Schwartz et al. 2006). An analysis of 14 core models (Fig. 6) indicates that the growing season will increase across the NA continent by the end of the century, although substantial variability in the magnitude of these changes exists on a regional basis. All changes are statistically significant at the 95% level relative to interannual variability in the observations (see Part I) with implications for impacts on agriculture and natural vegetation. The largest changes occur over the western United States and northern Mexico, where MEM increases of 40 days or more are projected. It should be noted that these same regions have complex terrain and are characterized by some of the largest negative biases in historical simulations (Part I), as well as the largest multimodel standard deviation in growing season length change of up to 8 days. In the central United States and Canada, increases of about 3–5 weeks are projected. The lengthening of the growing season is caused by both last spring freezes that are earlier and first autumn freezes that are later (not shown), but the change in the former is generally larger. A complementary analysis detailing changes in frost days is shown in the supplementary material.
e. Extreme events
To assess the projected changes in the frequency of occurrence (FOC) of extreme persistent dry–wet events over NA, we use the eight core models that contain three or more runs for both the historical and RCP8.5 experiments (see caption). Because these events are rare, three runs from each model are needed in order to produce enough events for meaningful statistics. The methodology used to define extreme events is the same as in Part I and repeated in the Fig. 7 caption. Note that because these calculations account for only liquid precipitation, results in the coldest northern regions are questionable and not discussed.
We first calculated the difference in FOC of extreme precipitation events using each experiment’s own climatology as a baseline for that experiment. Therefore, the historical climatology is from 1850 to 2005 and the RCP8.5 climatology is from 2006 to 2100. Each model shows little difference in the FOCs for both positive (wet) and negative (dry) events (not shown). The MEM difference between the positive and negative events (Fig. 7j) also shows no robust change between the model projections and historical data when the different climatologies are used. When the historical climatology is used as a baseline for the RCP8.5 experiment instead of the RCP8.5 climatology, all models but HadGEM2-ES show a decrease in the number of positive events in Mexico and the southwest United States and an increase in such events in the northeastern United States. The opposite is true for negative (dry) events. This shows the impact of the changing climatology between the experiments. These results are illustrated in Figs. 7a–h for individual models and in Fig. 7i for the MEM. The patterns among the models are similar, but they differ in magnitude and the precise location of the zero line. For example, IPSL-CM5A-LR shows an increase in dry events over the southern United States and Mexico, but the CCSM4 indicates that only Mexico is impacted. A metric showing the difference between area averages of the northeast quadrant (94°–75°W, 35°–48°N) and the southwest quadrant (123°–95°W, 15°–35°N) of NA is also given (Table 2). The metric values also have a small range between 0.31 and 0.48 when HadGEM2-ES is removed. With the spread taken into consideration, we conclude that more droughts are projected over Mexico and more persistent wet spells are projected over the northeast United States (Fig. 7i).
The projected changes of extreme surface temperature during 2081–2100 relative to 1981–2000 are shown in Fig. 8, calculated using one ensemble member from each of the 11 models described in the caption (10 are core models). Daily maximum surface temperatures Tmax are used to compute the number of days per year that exceed 90° (Fig. 8a) and 100°F (Fig. 8b), respectively. The MEM projections show an increase of 60%–300% (50–80 days) annually with Tmax warmer than 90°F in the southern United States and northern Mexico. In the southeastern United States, southern Texas, and northern Mexico, the number of days with Tmax warmer than 90°F is projected to increase to nearly 80. The MEM projections also show that the number of days with Tmax warmer than 100°F will increase 80%–400% (40–80 days) in parts of the south-central and southwestern United States. Across the southern United States and northern Mexico, the change in frequency of extreme surface temperatures are robust, suggested by the MEM projections having greater difference than the intermodel spread of the changes. However, greater uncertainty exists in other areas where MEM increases have the same magnitude as the standard deviation, in particular increases of 90°F days in the northeastern United States and northern Rockies and 100°F day changes across the northern half of the United States.
4. Regional climate
a. North Pacific and North Atlantic storm tracks
The projected change in Northern Hemisphere storm-track (ST) activity is examined based on 6-hourly data provided by 15 of the 17 core models (Fig. 9). Here, ST activity is defined based on meridional wind and SLP variance statistics computed using a 24-h difference filter (Wallace et al. 1988) that highlights the synoptic time scale (1.2–6 days).
Near the tropopause (250 hPa), the models project a strengthening of ST activity on the poleward flank of the historical ST peak and a slight weakening on the equatorward flank during winter. In summer, the models project a significant decrease in upper-tropospheric ST activity south of the ST peak and weak increase north of it. These results are consistent with previous studies based on CMIP3 (e.g., Yin 2005; Teng et al. 2008) that indicate a poleward shift of the ST under global warming. In the midtroposphere (500 hPa) and near the surface, the models generally project a significant decrease in ST activity extending from the Pacific across North America into the Atlantic in both seasons. This contrasting upper-level and near-surface change in winter can be related to contrasting projected change in the meridional temperature gradient. Near the surface, the temperature gradient is projected to significantly decrease because of greater warming at high latitudes, while the temperature gradient near the tropopause is projected to increase because of warming in the tropical upper troposphere and cooling in the polar lower stratosphere due to increases in greenhouse gases. The general decrease in ST activity in summer can be related to a projected decrease in mean available potential energy because of a decrease in the midlatitude temperature gradient and an increase in static stability.
As discussed in Chang et al. (2012), the significant near-surface ST activity decrease over NA represents one of the largest differences between CMIP5 and CMIP3 ST projections over the globe. As seen in Table 3, during winter (DJF), CMIP5 models project a −9.9% ± 3.6% MEM change in sea level pressure (SLP) variance in winter over the region roughly covering the contiguous United States and −19.8% ± 6.9% change in summer (JJA), with 14 out of 15 models projecting a decrease in winter and all 15 models projecting a decrease in summer. On the other hand, 11 CMIP3 models project −0.4% ± 4.0% change in the same quantity in winter [based on the Special Report on Emissions Scenarios (SRES) A2 scenario] and −9.2% ± 6.3% change in summer, with 7 out of 11 models projecting a decrease in winter and 10 of 11 models projecting a decrease in summer. We have also examined 8 other CMIP5 models not listed in Table 3, and all of them showed a decrease for both seasons. More details are given in Chang (2013), who showed that models projecting a larger decrease in ST activity over North America also project a more northward intrusion of the decrease in subtropical precipitation into southern United States.
We now provide a complementary analysis to that above using the Hodges (1994, 1995) cyclone-tracking scheme on 6-hourly mean SLP data to assess changes in extratropical cyclone activity along the U.S. East Coast. Part I presented the historical (1979–2004) predictions of western Atlantic extratropical cyclones during the cool season (November–March), which show substantial skill at simulating the distribution of cyclone activity, although with modest underprediction of amplitude. Colle et al. (2013) highlighted the details of this historical cyclone analysis and the twenty-first-century predictions in this region using these 15 CMIP5 models. Figures 10b–d show the MEM difference in cool season cyclone-track density for each of the three separate 30-yr future periods in RCP8.5 (2009–38, 2038–69, and 2069–98) and the historical period (1979–2004; Fig. 10a). Only a slight decrease in cyclone activity is projected over parts of the western Atlantic storm track for 2009–38 (Fig. 10b); however, Colle et al. (2013) show that this reduction may be more widespread if only the highest-resolution CMIP5 models are considered. The MEM reduction in cyclone density is more apparent for the 2038–69 period, with a reduction of 5%–15%, primarily along the southern half of the cyclone storm track, which is near the Gulf Stream boundary. Meanwhile, a slight increase in cyclone density is projected to the north over parts of northern New England and Nova Scotia, enhanced in the highest-resolution models (not shown). Future changes relative to the historical period continue to increase in size and amplitude (10%–20%) for the 2069–98 period. The results above are very similar if only the 12 of 15 models contained in the core model list (Table 1) are used. Figure S5 provides a commentary analysis of changes in cyclone intensity.
b. Northeast United States and western Atlantic precipitation
In Part I, 14 of the 17 core CMIP5 models in Table 1 were evaluated to determine how well they can simulate precipitation over the northeast United States and western Atlantic during the cool season (November–March) in the historical period (1979–2004). Figure 11 shows the MEM CMIP5 precipitation for this historical period, as well as the precipitation difference (in mm season−1 and percent change) for the early twenty-first century (2009–38) minus historical and late twenty-first century (2069–98) minus historical for these same models. For the early twenty-first century (Fig. 11b), the precipitation increases 5%–10% (10–30 mm) over the northeast United States. Less than a 5% increase occurs over the western Atlantic associated with the midlatitude storm track, while the largest percentage increase is over northeastern Canada (10%–20%). By the late twenty-first century (Fig. 11c), the largest increase of 35%–80% (40–100 mm season−1) occurs in eastern Canada. Over the northeast United States, the mean precipitation increases by 15%–25% by the late twenty-first century. The number of relatively heavy precipitation events (>25 mm day−1) over the northeast United States increases by 50% by the early twenty-first century and increases by 4–5 times by the late twenty-first century (Fig. 11d). These results suggest that the potential exists for a dramatic increase in the number of extreme rainfall events over the northeast United States during the next 50–75 yr.
c. Western water
Changes in snow are important in western regions of the United States because of the implications for water resources and winter tourism. Figure 12 shows changes in April SWE over the western United States from 15 core CMIP5 models for RCP8.5. April SWE is an indicator of total snow accumulation over the winter and the potential water resource availability for the coming year, which is especially important in California and the upper Colorado River basin (Fassnacht 2006). A MEM decrease in April SWE of up to 80 mm or greater is projected for the central Rockies and Canadian Rockies, with general decreases of smaller amplitude indicated elsewhere. The spatial resolution of the models is generally too coarse to represent snow-related processes throughout western North America because of the smoothed topography. For example, the details of the Sierra Nevada range on the California–Nevada border are absent. Nevertheless, this broad decline projected by the models is supported by high-resolution hydrological changes using downscaled projections (Hayhoe et al. 2004). The decrease in SWE is driven by higher temperatures that increase the ratio of rainfall to snowfall (wintertime precipitation is projected to remain the same or increase slightly; see Fig. 1) and increased melt efficiency, therefore moving the spring melt earlier in the season. The shift in snowmelt timing may also have consequences because water rights can be month dependent (Hayhoe et al. 2004).
d. North American monsoon
In Part I, the seasonal cycle of precipitation in 21 CMIP5 models over the historical period (1979–2005) was evaluated to identify models that have a reasonable precipitation climatology over the core NA monsoon (NAM) region (24°–29°N, 105°–109°W). The results of this analysis indicate that 9 of the 21 models have small (lag = 0) phase errors with respect to the observations. Here we present the projected behavior at 2070–99 for this subset of models that best simulate the historical precipitation climatology in this region.
The RMS difference of the future minus historical monthly rainfall climatology and the annual mean rainfall percent differences over the core NAM region are shown in Table 4 for these nine models. The results suggest that even for models correctly capturing the timing of the seasonal cycle of precipitation in the region, large differences exist in what these models project for the change in the monthly mean magnitude (range of 0.4–0.8 mm day−1) and the relative change in the overall annual mean (range from −34% to 3.7%) in precipitation for the monsoon region. However, seven out of nine models project that conditions for the NAM will be drier in the future under the RCP8.5 warming scenario suggesting some consistency in the sign of the change. Using a larger set of 16 core models that provided daily RCP8.5 rainfall, which include those with nonzero phase errors in the historical period, the change in the monthly mean magnitude of rainfall in the core monsoon region is somewhat greater, with a mean of 0.72 mm day−1 (range of 0.40–1.35 mm day−1), and the change in the annual mean rainfall indicates more drying, with a mean of −22.2% (range from −72.3% to 3.7%) compared to the −15.4% for the better performing models shown in Table 4.
In Fig. 13, MEM monthly precipitation from 16 core models are analyzed for a longitudinal belt (102.5°–115°W) from 7.5° to 35°N to assess the seasonal migration of precipitation in the NAM and ITCZ regions. These models capture the northward migration of precipitation between the ITCZ and the NAM region (north of 20°N) during the warm season, although the migration within the NAM region is less evident. The models’ historical precipitation (Fig. 13b) begins later and is stronger than the observed (Fig. 13a) in the NAM region but is weaker than the observed south of this region. Model projections from RCP8.5 (Fig. 13d) are consistent with the CMIP3 results (see introduction) and show reduced precipitation from 10° to 25°N through the cold season and extending into June and July, with increased precipitation in September and October. The monthly mean precipitation response from the individual models indicates strong model agreement on reduced December–July rainfall and a weaker consensus on increased late summer rainfall (Fig. 13c). These reductions in precipitation from 10° to 25°N in June and July are also consistent with reduced Mexico and Caribbean precipitation in this latitude band seen in Figs. 1 and 2 for JJA.
e. Great Plains low-level jet
The Great Plains (GP) low-level jet (LLJ) is a basic component of the warm season circulation in NA that provides a moisture source for GP precipitation. It emerges in April, strengthens and peaks in June and July, and vanishes in late September (e.g., Part I). Figure 14a reviews the ability of 16 of the core models (excluding HadCM3) to capture the intensity and seasonal cycle of meridional flow in the LLJ region. The MEM is best able to capture the amplitude of the LLJ during boreal summer, and vertical cross sections and horizontal maps of the jet verify this behavior (not shown). Figures 14b–e show how the jet strength evolves for the core models in the MEM and ±1 standard deviation limits for RCP4.5 and RCP8.5. The MEM projection indicates a strengthening of the LLJ by about 10%–20% during boreal summer by 2071–2100 in RCP8.5, with modest increases during 2035–64 and in RCP4.5. Notably, the lower one standard deviation bound is separable from zero in RCP8.5 for 2071–2100 during March–July. When the same analysis is applied to only models that produce the best simulation of the LLJ in current climate (Part I), the results are generally consistent.
The development of warm season precipitation anomalies or extreme wet and dry conditions in the GP and the Midwest is largely a result of the weather/precipitation systems that develop and depend strongly on the dynamic stability of the large-scale circulation (e.g., Moore et al. 2012; Veres and Hu 2013) and is not only a function of the strength of the LLJ. In fact, Klein et al. (2006) show that models that produce a reasonable LLJ may have trouble reproducing an accurate GP mean summertime rainfall distribution if they cannot correctly simulate eastward-propagating convective systems that develop over terrain during the diurnal cycle, since such systems produce 50% of summertime GP rainfall. This notion might explain the apparent discrepancy between the predicted strengthening LLJ shown in Fig. 14 and the decrease in summer precipitation in the GP shown in Figs. 1 and 2. This discrepancy suggests changes in the future summer circulation regime that would produce less precipitation in a moister environment with possibly a different intensity distribution of precipitation events in the GP, similar to what is observed in the northeast United States presented in section 4b.
f. Arctic sea ice
Significant reductions in the extent and thickness of the Arctic sea ice cover have occurred during the past several decades (e.g., Stroeve et al. 2012a,b; Kwok and Rothrock 2009), with wide-ranging impacts on marine ecosystems, coastal communities, prospects for resource extraction, and weather conditions in the Arctic and beyond (e.g., Francis and Vavrus 2012). While negative trends are observed in all calendar months, the largest reductions are observed in September at the end of the summer melt season. Through 2013, the rate of decline in September sea ice extent since 1979 has been −14.0% per decade, representing a reduction of more than 40% in the amount of sea ice covering the Arctic Ocean.
All models participating in CMIP5 correctly simulate a declining Arctic sea ice cover over the period of observations (Stroeve et al. 2012a). However, trends from most models remain smaller than observed. Figure 15 shows the CMIP5 MEM and ±1 standard deviation of the historical (gray line and shading) and future evolution of the sea ice cover under RCP4.5 (blue line and shading) and RCP8.5 (red line and shading). Observations based on a combination of satellite data and ship and aircraft observations (Meier et al. 2012) are shown in black. In constructing the MEM, the same criteria as used in Stroeve et al. (2012a) was applied, such that models that had more than 75% of their extents outside the observed range from 1953 to 1995 were excluded. Based on the extent criteria, seven models were excluded (CanCM4, CanESM2, CSIRO Mk3.6.0, EC-EARTH, FGOALS-s2, GISS-E2-R, and HadGEM2-AO). This resulted in 19 models retained for RCP4.5 and 18 models for RCP8.5. If we include only the core models in Table 1 without using the selection criterion applied above, conclusions do not significantly change from those described below.
Under RCP4.5, the MEM does not reach ice-free conditions by 2100, though the −1 standard deviation bound reaches nearly ice-free conditions (defined as less than 1 × 106 km2) around 2050. In contrast, the MEM ice extent drops below 1 × 106 km2 around 2060 for the RCP8.5 emission scenario, with the −1 standard deviation bound dropping below 1 × 106 km2 in 2030 and showing completely ice-free conditions by 2050. Note however that a large spread exists among CMIP5 models as to when a seasonally ice-free Arctic state may be realized. CMIP3 models were found to be very conservative with regards to Arctic sea ice loss (Stroeve et al. 2007, 2012a) despite having a realistic representation of the seasonal cycle in the ice cover. However, despite better overall agreement in CMIP5 with the historical observations, the spread in projected extent through the twenty-first century remains about the same between both modeling efforts, with an equivalent number of models showing summer ice-free conditions at the end of the century. While constraining future projections based on model performance does not significantly change the date when seasonally ice-free conditions will be reached, it does help reduce the spread in the projections.
Spatially, significant reductions in sea ice concentration off the coast of North America are expected during summer, even under the RCP4.5 emission scenario (Fig. 16). MEM mean sea ice concentrations for three decades are shown (2000–09, 2040–49, and 2090–99) together with a measure of the spread of the MEM, defined as the percentage of models that had more than 15% sea ice concentration for each grid cell. This measure conveys model agreement and, while the selection of a 15% sea ice concentration threshold is somewhat arbitrary, it is a useful metric for marine shipping and is generally used as a cutoff for sea ice extent calculations.
As early as the middle of the century, MEM sea ice concentrations decrease to less than 10% compared to somewhat higher sea ice concentrations of 10%–20% during the most recent decade. More importantly is that fewer models suggest sea ice concentrations in excess of 15% in the region. By the end of the century, the region with less than 10% sea ice concentration has grown and even fewer models suggest there will be ice with concentrations in excess of 15%. North of Greenland and the Canadian Archipelago, the last refuges of old thick ice, MEM ice concentrations in September are less than 50% by the end of the century. In contrast, during winter (e.g., March), no change in sea ice concentrations along coastal North America is expected under the RCP4.5 emission scenario (Fig. 17). The same is true for most of the Arctic, though the number of models indicating at least 15% sea ice concentration in the North Atlantic declines throughout the century.
g. Warming hole in the eastern United States
During the second half of the twentieth century, the central-eastern United States experienced cooling trends while global mean temperatures warmed. We refer to this cooling region as a “warming hole” (WH), following Pan et al. (2004). A number of studies have attributed the mechanisms for this abnormal trend to large-scale decadal oscillations such as the Pacific decadal oscillation (PDO) and Atlantic multidecadal oscillation (AMO; Robinson et al. 2002; Kunkel et al. 2006; Wang et al. 2009; Meehl et al. 2012a). Other studies indicate that regional-scale processes such as hydrological cycle (Pan et al. 2004) and land surface interaction (Liang et al. 2006) may contribute to the WH. Leibensperger et al. (2012) have attributed the WH to anthropogenic aerosol forcing. While the atmosphere–ocean coupled simulations do not allow for distinguishing different types of forcings, neither CMIP3 nor CMIP5 twentieth-century simulations show the WH as an externally forced response signal (Kunkel et al. 2006; Kumar et al. 2013; Pan et al. 2013). Kumar et al. (2013) found that observed temperature trend variability in the eastern United States is significantly correlated with the AMO. In Kumar et al. (2013) and Part II, it was shown that the 95% uncertainty range of historical CMIP5 simulations brackets the observed negative temperature trend, although the MEM time series has limited skill at reproducing the WH.
Here we assess whether evidence of an eastern United States WH exists in twenty-first-century projections. Figure 18 shows the MEM annual temperature trends during 2015–50 and 2051–98 for RCP4.5 and RCP8.5 in the core models (Table 1; one ensemble member each). The MEM shows a warming trend in all regions of NA. In the first half of the twenty-first century, the warming rate is 28% higher in RCP8.5 (0.55°C decade−1) than in RCP4.5 (0.43°C decade−1). In the second half of the twenty-first century (2051–98), the warming rate is more than 4 times higher in RCP8.5 (0.73°C decade−1) than in RCP4.5 (0.17°C decade−1). The 50% reduction in the late-twenty-first-century warming rate compared to the first half of the twenty-first century is consistent with stabilization of CO2 emissions in RCP4.5 (Moss et al. 2010; Meehl et al. 2012b). In both RCP4.5 and RCP8.5 the projected warming rate is significantly higher than the twentieth-century warming rate in CMIP5 simulations (0.07°C decade−1; Kumar et al. 2013).
Figures 19a and 19b show time series of 30-yr running trends of eastern United States annual temperature for RCP4.5 and RCP8.5 simulations from the core models including all available ensemble members (a total of 55 ensemble members in RCP4.5 simulations, and 46 ensemble members in RCP8.5 simulations). The multimodel ensemble median of RCP8.5 simulations indicates a continued increase in warming rate from 0.4°C decade−1 at the start of the twenty-first century to 0.7°C decade−1 toward the end of the twenty-first century, whereas RCP4.5 simulations indicates a decline in warming rate from 0.3° to 0.1°C decade−1 over the twenty-first century. The entire 95% uncertainty range in RCP8.5, as well as the majority of the 95% uncertainty range (>90%) in RCP4.5, is above the zero line for the 30-yr running annual temperature trend in eastern United States (Figs. 19a,b). Kumar et al. (2013) found similar results for summer (JJA) temperature trends. Hence, the negative temperature trend in eastern United States is shown to disappear under RCP4.5 projections with 90% probability and to disappear under RCP8.5 projections with 100% probability. We did not find a significant difference in trends between the eastern and western United States (defined in Fig. 18) in the twenty-first-century climate projections (not shown).
5. Tropical intraseasonal variability
a. Midsummer drought
For most of Central America and southern Mexico, climatological precipitation has a maximum in June and September, bracketing a period of reduced rainfall during July–August known as the midsummer drought (MSD; Portig 1961; Magaña et al. 1999). Hence, this variability in the annual cycle represents a climatological intraseasonal oscillation, akin to those found in other regions of the globe (e.g., Wang and Xu 1997). A previous assessment of CMIP3 model performance at simulating the MSD and future projections (Rauscher et al. 2008) suggested that many CMIP3 models are capable of simulating the MSD despite an overall dry bias and that the MSD is projected to become stronger with an earlier onset. An updated evaluation of this feature (Part II) indicates that many CMIP5 models are also able to capture both the spatial and temporal aspects of the MSD. This success is promising, as accurately simulating or including all of the air–sea interaction processes relevant to the MSD is difficult.
The CMIP5 MEM projection for precipitation change during each of the summer months (June through September) is shown in Fig. 20 for the core models in Table 1. During each of the summer months, the east Pacific ITCZ is projected to shift southward in concert with a drying over the east Pacific warm pool (EPWP), Central America/southern Mexico, and the Caribbean with enhanced drying over the Caribbean islands of Cuba, Hispaniola, and Jamaica. The strongest drying is projected to occur during July and August, which are the months when the MSD occurs in many regions throughout the inter-Americas region. The precipitation changes are in general consistent with the JJA average precipitation patterns shown in Fig. 1, although western Mexico is projected to experience wetter conditions in September, consistent with the analysis of Fig. 13.
Next, a simple algorithm for quantifying the MSD strength is used that does not assume a priori which months are climatological maxima and which months constitute the MSD (Karnauskas et al. 2013). Consistent with the month-by-month rainfall projections, the CMIP5 core MEM provides a very robust projection of a stronger MSD for most regions that experience an MSD today (Fig. 21). The maximum MSD increases from ~2.5 mm day−1 to ~(3–4) mm day−1 in the RCP4.5 forcing experiment. The peak amplitude increases slightly more in the RCP8.5 forcing experiment, and areas of amplitude greater than 3 mm day−1 substantially increase in spatial extent. The projection in each of the CMIP5 models that best replicates the observed MSD is qualitatively consistent with the MEM projection (Part II). The stronger MSD is a result of early and midsummer rainfall being reduced relative to the late summer peak (see also Fig. 13). The extension of the MSD northward along the Gulf coast of Mexico and into the United States is projected to strengthen in both RCP4.5 and RCP8.5 experiments. Some regions that did not previously exhibit a MSD (e.g., Panama) develop a moderate MSD under both forcing experiments.
b. Transient intraseasonal oscillations
Some have argued that the leading mode of intraseasonal variability (ISV) in the east Pacific warm pool is a regional manifestation of the Madden–Julian oscillation (Maloney et al. 2008). Jiang et al. (2013) and Part II documented the ability of CMIP5 models to capture this leading mode of 30–90-day precipitation variability over the east Pacific and adjacent regions of Mexico, Central America, and the Caribbean during June–September. This analysis was done using complex empirical orthogonal function (CEOF; Barnett 1983; Horel 1984) analysis that was previously used to document the leading mode of ISV over the east Pacific during boreal summer (Maloney et al. 2008). Based on a pattern correlation and assessment of amplitude of the leading CEOF as compared to TRMM observations, eight models were judged to produce realistic ISV.
Figure 22 shows a seven-model MEM amplitude of the leading CEOF mode for 1981–2005, the 2076–2100 projection in RCP8.5, and the difference in amplitude relative to present day. As indicated in the caption, these models include six core models that were assessed by Part II to produce realistic ISV, along with one additional core model that produces poor ISV (CCSM4), although results do not differ if only good performing models are used. Stippling on the difference plot (Fig. 22c) shows where amplitude increases are statistically significant from zero at the 95% confidence level. Robust changes include significantly increased MEM amplitude of the leading CEOF on the southern fringe of the amplitude maximum, with inconsistent changes in amplitude elsewhere. A plot of intraseasonal precipitation variance supports these results (Jiang et al. 2013). These regions of significantly increased variance coincide with areas of mean precipitation increase shown in Figs. 1 and 20, a tendency noted in other studies that have examined projected ISV increases in a warming climate (Maloney and Xie 2013).
6. East Pacific and Atlantic tropical storm-track and cyclone activity
a. Tropical storm track
Traveling tropical synoptic-scale disturbances found in tropical storm tracks (e.g., Thorncroft and Hodges 2001; Serra et al. 2008, 2010) serve as the precursors to a majority of TCs in the Atlantic and eastern North Pacific. Their frequency at 850 hPa over Africa and the eastern Atlantic has been shown to be positively correlated with Atlantic hurricane activity (Thorncroft and Hodges 2001). Projected changes in the density and strength of these disturbances may therefore contribute to our understanding of how the statistics of TCs may change in the future, with the advantage that these synoptic systems are better resolved by low-resolution global models than the TCs themselves (see section 6b). As described in Part II, storm-track statistics are calculated from smoothed, 6-hourly 850-hPa relative vorticity following the method of Hodges (1994, 1995) for nine CMIP5 models. This method tracks 6-hourly 850-hPa positive vorticity centers with a minimum threshold of 0.5 × 10−6 s−1 that persist for at least 2 days and have tracks of at least 1000 km in length. This method primarily identifies westward moving disturbances such as easterly waves (e.g., Serra et al. 2010), although more intense storms that could potentially reach hurricane intensity are not excluded. Track density is calculated as the number of tracks in a 5° spherical cap per month while mean track strength is the mean vorticity of those tracks at each location. Comparisons of the CMIP5 track statistics with the European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-Interim) for these same models were presented in Part II, which found that the best performing models within this subset are the HadGEM2-ES and CNRM-CM5.1 while the worst performing models are the CCSM4, BCC_CSM1.1, and GFDL-ESM2M.
Figure 23 shows the RCP8.5 nine-model MEM and standard deviation of the mean (SDOM) for track density and mean strength for the May–November of 2070–2100 period, as well as the future minus historical period differences in these quantities. Overall, the models indicate a southward shift of the main storm track as well as an increase in track density for the future projections, consistent with results found using CMIP3 models (Bengtsson et al. 2006, 2007; Colbert et al. 2013). This southward shift is also observed in the pattern of mean precipitation change seen in the JJA analysis in Figs. 1 and 20. Notably, the two best performing models for the historical period, HadGEM2-ES and CNRM-CM5.1, do not show strong agreement on the magnitude of the projected future track density change (not shown). Changes in the mean track strength are less significant, with a small region of decreased strength seen in the west Atlantic and a small region of increased strength seen in the east Pacific.
b. Tropical cyclones
Using the Camargo and Zebiak (2002) algorithm, we analyzed changes in TC-like vortices in RCP4.5 and RCP8.5 in 14 core CMIP5 models, as shown in the figures of the supplementary material section (e). As discussed in Part II, a subset of these models was shown to produce the best simulation of TC activity in the current climate (GFDL CM3, GFDL-ESM2M, MPI-ESM-LR, MRI-CGCM3, and MIROC5), although all of these models still severely underestimate the number of TCs per year in the Atlantic and east Pacific. The definition of the model TCs includes thresholds for low-level vorticity and wind speed (model and resolution dependent), as well as imposes a constraint of warm core and that the vortex is an SLP minimum. Details of this analysis can be found in Camargo (2013).
Figure 24 shows the mean track density in the North Atlantic for the historical simulations for the period 1951–2000, as well as the difference of the track density between the period 2051–2100 and the historical track density for RCP4.5 and RCP8.5. While the MRI models shows a slight increase in TC track density for the Atlantic north of 25°N, the MIROC model is characterized by substantial decreases along the U.S. coastline and the GFDL CM3 and the MPI models show a small northward shift in the track density. To be more quantitative, Fig. 25 shows the number of TCs in the Atlantic and eastern north Pacific for these five models in the historical runs, RCP4.5, and RCP8.5. For the North Atlantic, the MRI model exhibits an increase in Atlantic TC numbers in RCP8.5 relative to the historical period, whereas the MIROC5 model exhibits a significant decrease in Atlantic TC numbers in both RCP4.5 and RCP8.5. For the east Pacific, results are also mixed, with significant increases in TC numbers in future climate in the MPI model in both RCP4.5 and RCP8.5 and significant decreases in the MIROC5 model. This inconsistency in future TC changes among models is supported by use of a high-resolution model as a downscaling tool to examine projected changes in hurricane activity, as described in the supplementary material.
In the context of intensity, the CMIP5 models are not able to simulate the most intense TCs due to model resolution. For the seven models that have the highest global number of TCs, a statistically significant difference is found in the maximum wind speed distributions of the RCP4.5 and RCP8.5 simulations compared with the historical simulation, based on the Kolmogorov–Smirnov test at the 99% significance level. For these seven models (CSIRO Mk3.6.0, GFDL CM3, GFDL-ESM2M, IPSL-CM5A-LR, MIROC5, MPI-ESM-LR, and MRI-CGCM3), increases in the values of the maximum wind speed occur for various percentiles of the wind speed distribution (50th, 75th, and 90th percentiles), with six (five) of the models having increases in the maximum wind speed values in the 95th (99th) percentile.
c. Analysis of Atlantic wind shear and relative SST changes
We provide further diagnosis here to help explain some of the intermodel inconsistency in Atlantic TC activity in the twenty-first century. Zhao and Held (2012) found that most of the intermodel spread in the North Atlantic TC frequency response among the CMIP3 models can be explained by a simple relative SST (RSST) index defined as the tropical Atlantic SST minus tropical mean SST. Under global warming scenarios the SST difference between the Atlantic main development region (MDR; 85°–15°W, 10°–20°N) and the other tropical ocean basins varies from model to model with implications for TC activity (Latif et al. 2007; Swanson 2008; Vecchi et al. 2008; Wang and Lee 2008; Xie et al. 2010). Other studies suggest the strong influence of Atlantic vertical shear (VS) variations in future climate on TC activity (e.g., Vecchi and Soden 2007). Here we assess the trends over the twenty-first century in Atlantic MDR RSST and VS in CMIP5 models for RCP8.5.
A scatterplot of the MDR RSST and VS trends in the twenty-first-century CMIP5 simulations is shown in Fig. 26. The conclusions are similar if using only a subset including the core models. The individual models show different responses of the RSST and VS trends, and hence suggest inconsistent changes in TC activity, consistent with the findings in section 6b. As shown in the supplementary material, similar uncertainty in the sign and magnitude of VS and RSST change is seen in the east Pacific. The linear fit of all models (the solid line in Fig. 26) shows that the VS trend decreases with the RSST trend. That is, if the MDR SST trend under global warming is smaller than that of the global tropical ocean, the VS trend in the MDR is increased and vice versa. This is consistent with a recent modeling study (Lee et al. 2011), showing that a slower warming in the tropical North Atlantic compared with the tropical Indo-Pacific Oceans increases the VS and static stability in the MDR for Atlantic TCs. The preferential warming of the tropical Indo-Pacific Oceans induces a warming of the tropical troposphere over the tropical North Atlantic, via a tropical teleconnection mechanism, and thus increases atmospheric static stability and decreases convection over the tropical North Atlantic. The slower warming (or the relative cooling) of the tropical North Atlantic induces an increase of the VS in the Atlantic MDR region (Wang et al. 2008). Because of the increased VS and static stability in the MDR, the number of Atlantic TCs in the twenty-first century is projected to decrease in the lower RSST and higher VS conditions predicted by some models.
We note, however, that the projected MIROC5, MRI-CGCM3, GFDL, and MPI-ESM-LR changes shown in Fig. 25 derived from a direct TC-tracking method do not scale as expected from Fig. 26. This result may highlight the limitations of the direct tracking approach described in section 6b. Indeed, in the supplementary material we describe results using a high-resolution model as a downscaling tool in which the projected changes in TC activity in the GFDL model come more in line with those suggested by Fig. 26. It is also described in the supplementary material that, at very high greenhouse gas concentrations, the larger direct effect of the atmospheric greenhouse gas concentrations can suppress global and regional TC frequency and thus weaken the relationship between RSST and TC activity (Held and Zhao 2011).
7. Interannual to decadal hydroclimate
a. Interannual variability
The large internal variability of the El Niño–Southern Oscillation (ENSO), a major source of interannual NA hydroclimate variability, makes it difficult to detect changes in ENSO properties over periods of a couple of centuries (Wittenberg 2009; Stevenson et al. 2012). Even without significant changes in ENSO itself, the long-term trend in tropical SST and the associated changes in the midlatitude basic state may result in substantial changes in ENSO-related teleconnections (Meehl et al. 2006; Meehl and Teng 2007; Lau et al. 2008; Kug et al. 2010; Kumar et al. 2010; Stevenson et al. 2012). In this analysis, we evaluate whether the NA seasonal temperature and precipitation patterns associated with ENSO are projected to change significantly in the twenty-first century in CMIP5 models in RCP4.5 and RCP8.5. Performance in accurately simulating ENSO teleconnections in current climate is relatively strong in several of the CMIP5 models and particularly in the MEM (Part II). We do not distinguish different types of ENSO episodes, such as central Pacific from east Pacific El Niño episodes that may change in proportion with climate change (Yeh et al. 2009), but rather consider a single broad class of ENSO events identified by the Niño-3.4 SST index. The Niño-3.4 SST index is defined as the average SST in the region 5°N–5°S, 120°–170°W (Trenberth 1997). We identify ENSO episodes in the same way as described in the historical analysis (Part II). Because all anomalies are detrended for this analysis, we focus on interannual variability superimposed upon the long-term trend. On the basis of these calculations, we find that the frequency of ENSO episodes in the RCP4.5 and RCP8.5 simulations is approximately the same as in the historical simulations. A hint of increased ENSO amplitude exists, as the ensemble and episode mean peak Niño-3.4 SST amplitude of El Niño episodes for the historical, RCP4.5, and RCP8.5 episodes is 1.51°, 1.51°, and 1.59°C, respectively; for La Niña episodes, the mean peak amplitude is 1.39°, 1.46°, and 1.57°C, respectively. However, large intermodel variability exists, so these changes are not robust.
We calculate MEM 300-hPa geopotential height, surface air temperature (SAT), precipitation, and tropical surface temperature composites for RCP8.5 in the same way as in the historical analysis (Part II) and then show the RCP8.5 (2006–2100) minus historical (1850–2005) composite differences in Fig. 27 (see Table 5 for the list of 21 models used that encompass the core set). Comparison of Fig. 27 with the historical composites for El Niño (Fig. 9 in Part II) indicates a projected strengthening and slight northeastward shift of the El Niño and La Niña teleconnection patterns in the twenty-first century. These projected changes include a strengthening of the negative southeastern United States temperature anomalies during El Niño (Fig. 27c), warming over much of NA during La Niña (Fig. 27d), and a strengthening of the precipitation patterns on the west and southeastern coastal regions of the United States (Figs. 27e,f). These strengthened patterns are consistent with significant increases in the tropical convective forcing (Figs. 27i,j). The changes in ENSO SST anomalies are modest (Figs. 27g,h), and so the increase in the convection anomalies may relate more strongly to the warming of the climatological SSTs, particularly in the eastern equatorial Pacific.
The projected northeastward shift of the Aleutian anomaly during El Niño has been noted in previous studies (Meehl and Teng 2007; Kug et al. 2010), but not all models exhibit this change. As noted in Stevenson (2012), some models project this northeastward shift whereas others tend to indicate a strengthening in place. This uncertainty is reflected in high ensemble standard deviations of 300-hPa height over the North Pacific (Fig. 27a), which then contributes to the high ensemble standard deviations in projected temperature changes over northwest NA (Fig. 27c). For the La Niña teleconnections, the projected changes in the geopotential height field (Fig. 27b) resemble the response to Indo–western Pacific warming shown in Lau et al. (2008), which suggests that the changes may be a direct response to the increased Indo–western Pacific convection anomalies (Fig. 27j). For both the El Niño and La Niña teleconnection changes, the highest intermodel variance generally corresponds with the regions of highest precipitation change (Figs. 27e,f).
Table 5 provides centered pattern correlations and RMS differences between the historical and twenty-first-century RCP8.5 seasonal composites of SAT and precipitation over the NA region (6°–75°N, 180°–50°W) for one run of each CMIP5 model. In addition, we indicate whether the composite patterns are significantly different based on a “false discovery rate” (FDR) (Benjamini and Hochberg 1995; Wilks 2006) field significance test, as discussed in Table 5. Table 5 indicates a large degree of intermodel spread for each seasonal composite, but the total ensemble pattern correlations are high for each composite and yet some of the projected changes are statistically significant. This finding reinforces the tendency for projected strengthening of the MEM teleconnections with little change in pattern shape.
b. Decadal variability
Similar to our analysis for ENSO, we evaluate whether NA temperature and precipitation patterns associated with the PDO are projected to change in the CMIP5 models. We define the PDO as the leading empirical orthogonal function of extended winter (November–April) monthly mean sea surface temperature (SST) anomalies in the North Pacific poleward of 20°N. We follow the same procedures for calculating SST anomalies, performing the EOF analysis, and calculating the PDO index as described in the historical analysis (Part II), except that we now base the EOF analysis on the period from 2006 to 2090. As in the preceding analysis of ENSO teleconnection changes, we examine changes in decadal variability superimposed on the long-term trend by removing the global mean SST for the PDO EOF calculations and by detrending the NA temperature and precipitation anomalies prior to the regressions. Part II showed that CMIP5 models have success in reproducing the PDO temperature teleconnection pattern but mixed success at capturing the teleconnection pattern in precipitation.
Overall, we find that SAT and precipitation pattern changes for individual models are generally not significant, but some robust changes are apparent in the ensemble mean, particularly for the temperature regressions (Table 6). Table 6 indicates that, although the ensemble-mean pattern correlations between the historical and RCP8.5 twenty-first-century temperature regressions are quite high, the projected pattern changes are statistically significant at least at the 5% significance level for all seasons. Inspection of the regression pattern changes reveals that the PDO-related temperature patterns are projected to weaken in the twenty-first century. Figure 28 shows the projected changes in the RCP8.5 boreal spring SAT regressions, a season that features robust changes under both emissions pathways. This figure indicates a clear weakening of the PDO-related SAT variability in the warmer climate, with largest weakening over western North America. Similar changes are noted in all other seasons and under the RCP4.5 emissions pathway (not shown), although the RCP4.5 changes are weaker in magnitude and statistically significant only in the spring. We find that the PDO-related North Pacific SST variability does not show a similar weakening, however.
8. Conclusions and discussion
We have examined twenty-first-century projections of NA climate in CMIP5 models under RCP8.5, including RCP4.5 in some cases. In terms of robust projected changes relative to historical runs, models indicate that warming will occur everywhere in NA, with the greatest warming occurring in northern Canada, where MEM temperature increases of 15°C can be found during boreal winter at the end of the twenty-first century. Models suggest a very likely disappearance of the negative temperature trend (warming hole) in the eastern United States in the twenty-first century. Models agree that more precipitation will fall north of 40°N in winter, with the largest increases along the west coast and northeast United States and Canada. The potential exists for a large increase in the number of extreme rainfall events over the northeast United States, since many of the CMIP5 models suggest a 4–5-fold increase in heavy precipitation events by the late twenty-first century. Model agreement is good on reduced summertime precipitation in the Caribbean and southern Mexico. Projected precipitation increases are robust across models for Alaska and the Arctic rim during boreal summer.
Further robust projected changes include a decrease in snow water equivalent throughout the year, with the greatest decreases in January–April. Especially large decreases are noted in the Rockies and northwest Canada. The growing season is projected to lengthen across NA by the end of the twenty-first century, with the largest changes in the western United States and northern Mexico, where increases of 40 days or more are projected. In the agriculturally important regions of the central United States and Canada, increases of 3–4 weeks are projected. The projections show an increase in the frequency of extended wet periods for the eastern United States and extended dry periods for Mexico. Warm temperature extremes are projected to increase across NA. Particularly notable is an increase of 40–80 100°F days in parts of the south-central and southwest United States relative to the historical climatology. Models project an increase in Great Plains lower-level jet strength in the twenty-first century. A poleward shift of the Atlantic storm track is indicated in the upper troposphere during all seasons. Near-surface storm-track and cyclone activity are projected to decrease over NA, except along the United States and Canadian east coast. Models agree on a drier early summer (June and July) and wetter fall (September and October) from the Pacific ITCZ north to the NAM region, accompanied by robust strengthening of the midsummer drought in Central America and the greater Caribbean region.
Although many projected changes in NA climate are robust across CMIP5 models, substantial disagreement exists in other areas. The sign of mean precipitation changes across the southern United States is inconsistent among the models, as is the annual mean precipitation change in the core NA monsoon region. Models also disagree on snow water equivalent changes on a regional basis, especially in transitional regions where competing effects occur because of greater snowfall and warming temperatures. The western United States is characterized by large intermodel variability in changes in the number of frost days, where MEM decreases in frost days (greater than 40 days) are also largest. Substantial intermodel spead exists for projections of how ENSO teleconnection changes will affect precipitation and temperature variability in western NA. Projected changes in seasonal mean Atlantic and east Pacific TC activity are inconsistent among models, which disagree on the sign and amplitude of changes in environmental factors that modulate TC activity. For example, changes in Atlantic and east Pacific SSTs relative to the tropical mean differ among models. Xie et al. (2010) imply that differences in model physics that regulate how ocean heat transport and wind–evaporation feedbacks may change in future climate may be partially responsible for different patterns of SST change and hence differing relative SST changes among models. Understanding how and why the pattern of SST change varies among models may provide a way to constrain future projections of east Pacific and Atlantic TC activity. It should be noted that even areas of substantial agreement among models may not imply more confidence that projections are correct, as common errors or deficiencies in model parameterizations may provide false confidence in the robustness of future projections.
Model biases in some areas decrease confidence in projections, including changes of the timing of NA monsoon precipitation that exhibits projected changes in the same sense as historical biases. Growing season length is projected to increase most strongly along the West Coast, where models tend to display large negative biases in historical runs. Models have substantial difficulties in simulating the historical distribution of persistent drought and wet spells, which might produce less confidence in the pattern of extreme precipitation event changes in future climate. However, Pierce et al. (2009) found that model success in producing historical climate has little bearing on regional projections. Further, the MEM tends to provide superior performance when compared to an individual model. Hence, MEM statistics were presented for most of the analyses shown here and may provide the best estimate of projected climate changes.
Finally, this paper represents an overview and only first attempt at integrating results for CMIP5 projections of North American climate, and many gaps remain in characterizing NA climate change. Further characterization of changes in projections from CMIP3 to CMIP5 is necessary, including how differences in concentration pathways versus model physical differences dictate projection differences between the two archives. More precise characterization and attribution of regional precipitation changes is also needed, including changes in extremes, as such projections are of particular importance for agricultural and domestic water supplies, hydropower generation, ecosystems, and recreational interests. For this task, it is likely that regional downscaling will provide a useful complement to these CMIP5 results. Further, extended process-oriented diagnosis of climate models is necessary to provide confidence that models are simulating modern climate for the correct reasons, which provides greater confidence in future projections.
We acknowledge three reviewers whose insightful comments led to significant improvements in the manuscript. We acknowledge the World Climate Research Programme’s Working Group on Coupled Modeling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 1 of this paper) for producing and making available their model output. For CMIP, the U.S. Department of Energy’s Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. The authors acknowledge the support of NOAA Climate Program Office Modeling, Analysis, Predictions and Projections (MAPP) Program as part of the CMIP5 Task Force. EDM was also supported by the Climate and Large-Scale Dynamics Program of the National Science Foundation under Grant AGS-1025584.
Supplemental information related to this paper is available at the Journals Online website: http://dx.doi.org/10.1175/JCLI-D-13-00273.s1.
This article is included in the North American Climate in CMIP5 Experiments special collection.