1. Introduction
The projected damage arising from tropical cyclones (TCs) in a future climate, particularly as anthropogenic global warming continues, is a topic of scientific and societal interest, and will be influenced by changes in storm intensity, population, and vulnerability (e.g., Mendelsohn et al. 2012; Peduzzi et al. 2012). The most intense hurricanes (category 3–5) are responsible for the vast majority of the tropical cyclone damage in the United States, even though they represent only one-fourth of the overall landfalling tropical cyclone activity (e.g., Pielke et al. 2008).
Theoretical considerations (e.g., Emanuel 1987; Holland 1997) and high-resolution modeling studies (e.g., Knutson and Tuleya 2004; Oouchi et al. 2006; Emanuel et al. 2008; Bender et al. 2010; Knutson et al. 2010) generally suggest an increase in the intensity of tropical cyclones in a warming climate. High-resolution models can represent the most intense storms directly, but the required computational expense generally limits them to single-model runs and/or time slice experiments. At the present time, it is unclear what outcome we would get by running these high-resolution models in a multimodel fashion over the entire twenty-first century.
An alternative approach to counting the number of the most intense storms is to employ the seasonally integrated power dissipation index (PDI; Emanuel 2005, 2007), which is a metric that convolves storm duration, frequency, and intensity. Storm intensity is accounted for by taking the third power of the wind speed. Emanuel (2005) found that there was a high correlation between PDI and tropical Atlantic sea surface temperature (SST). Swanson (2008) obtained a higher correlation using the difference between tropical Atlantic SSTAtl and tropical-mean SSTTrop. Vecchi et al. (2008b) showed that projections of PDI based on relative SST (the difference between SSTAtl and SSTTrop) are in better agreement with results from dynamical models than using SSTAtl alone. Recently, Villarini and Vecchi (2012a) developed a statistical model in which SSTAtl and SSTTrop are used as predictors (see section 2). We have also recently shown that this statistical model can be used to make retrospective skillful forecasts of seasonally integrated North Atlantic PDI from November of the previous season, allowing the skillful forecast of the upcoming season as the current one is still coming to an end (Villarini and Vecchi 2013).
Here, we apply the model of Villarini and Vecchi (2012a) to outputs from 17 global climate models (GCMs) produced under phase 5 of the Coupled Model Intercomparison Project (CMIP5; Taylor et al. 2012; Table 1) to address questions related to future changes in North Atlantic tropical cyclone intensity.
Summary of the 17 GCMs used in this study. For all of them, data for the RCP 2.6, RCP 4.5, and RCP 8.5 are available. The same holds true for the 2 × CO2 runs with the exception of MIROC-ESM-CHEM. The final column indicates the models for which the “greenhouse only” historical experiments were available (see Fig. 5). Here, Y indicates available.
2. Methodology
This parsimonious model can describe well the interannual and decadal variability and change of the PDI record over the period 1949–2008 (cf. red and blue lines in Fig. 1; Villarini and Vecchi 2012a) and allows for century-scale reconstructions of PDI (orange line Fig. 1).
The statistical frameworks modeling Atlantic hurricane activity using SSTAtl relative to SSTTrop, rather than SSTAtl alone, are supported by both modeling and empirical results (e.g., Latif et al. 2007; Vecchi and Soden 2007a; Swanson 2008; Bender et al. 2010; Zhao et al. 2010; Ramsay and Sobel 2011; Vecchi et al. 2008b, 2011; Villarini et al. 2010; 2011, 2012) and are the basis of skillful seasonal forecasts (Vecchi et al. 2011; Zhao et al. 2010; Villarini and Vecchi 2013). In the model of Villarini and Vecchi (2012a), the positive coefficient on Atlantic SSTs is larger than the negative coefficient on tropical-mean SSTs [Eq. (2)], and a similar difference in the magnitude of the Atlantic and tropical-mean coefficients to the fit of PDI was found by Swanson (2008). This indicates that uniform warming (cooling) of the tropics should lead to an increase (decrease) in PDI. In addition, warming (cooling) of the Atlantic relative to the tropical average will also lead to an increase (decrease) in PDI.
In this study, we examine projected changes in PDI by applying the statistical model of Villarini and Vecchi (2012a) to outputs from 17 GCMs (see Table 1 for a list) to address questions related to future changes in North Atlantic tropical cyclone intensity, using 1986–2005 as our reference period to compute anomalies, and the median of the gamma distribution as reference value. We have recently analyzed these GCMs for possible changes in tropical cyclone activity (Villarini and Vecchi 2012b). We showed that, while both tropical Atlantic and tropical-mean SSTs are projected to increase over the twenty-first century, there is a significant radiatively forced increase in North Atlantic tropical storm frequency only over the first half of the century, which is not driven by CO2 but is likely the results of anthropogenic aerosol changes (with a particular influence of a projected clearing of aerosols over the North Atlantic). This increase, however, does not extend over the entire twenty-first century, for which the sign of the trend is uncertain. Differences in the behavior of PDI and tropical cyclone frequency indicate changes in tropical cyclone intensity and duration at the strongest intensities (Villarini and Vecchi 2012a).
3. Results
Figure 1 shows the time series of PDI anomalies for three representative concentration pathways (RCPs; each RCP is labeled to reflect the radiative forcing change at the end of the twenty-first century, in watts per meter squared). The observations over the period 1949–2008 are within the GCMs ensemble spread. The projected PDI values tend to be larger than the PDI values over the last 130 years. Regardless of the RCP, there is a tendency toward increases in PDI over the twenty-first century. There are statistically significant trends in six out of the 17 models over the 2006–99 period for RCP 2.6 (five positive and one negative trends). This number increases to nine and 13 for RCP 4.5 and 8.5, respectively, all of which are increasing trends. The magnitude of these increases depends on the RCP, with RCP 2.6 showing smaller increases relative to RCP 8.5. This magnitude increase is coupled with an increase in the ensemble spread (Figs. 1 and 2). The projected mean over 2016–35 is, on average, 20% larger than the corresponding values over 1986–2005 (Fig. 2). Averaged over 2046–65, the multimodel mean PDI values are between 50% and 75% larger than the values for the reference period. At the end of the twenty-first century, the magnitude of the changes shows a positive dependence on the strength of projected greenhouse gas (GHG) forcing.
With the exception of RCP 2.6, which has a midcentury maximum in GHG forcing, the largest projected increases in PDI are toward the end of the twenty-first century. For RCP 4.5 (approximately a CO2 doubling at the end of the twenty-first century) the PDI values over 2080–99 are on average 50% larger than the reference period. Meanwhile, the PDI values for RCP 8.5 (approximately a quadrupling of CO2 at the end of the twenty-first century) are, on average, 100% larger than over 1986–2005, with a large increase in intermodel variability.
Storm duration, frequency, and intensity are used to compute the PDI. Could it be that the increases in PDI are reflecting an increase in tropical cyclone frequency? A recent analysis of the same GCMs (Villarini and Vecchi 2012b) showed that the ensemble-mean North Atlantic tropical cyclone frequency is not projected to change significantly over the entire twenty-first century, regardless of the RCP. In further contrast to the PDI projections, frequency increases that were found were the largest in the first half of the twenty-first century (driven by aerosol changes) and showed no relation to GHG forcing. The increase in PDI, therefore, indicates a projection for an increase in North Atlantic tropical cyclone intensity or the duration over which these storms are at their strongest intensities. A scaling argument in Villarini and Vecchi (2012a) suggests that these projections include an increase in the time during which North Atlantic tropical cyclones are at their strongest intensities; that is, all other things equal, the probability of a major hurricane occupying a place at any given time is projected to increase. This is related to the rather uniform SST increase in the Atlantic Ocean and the tropics (Villarini and Vecchi 2012b), resulting in small changes in the frequency of tropical storms and in an increase in PDI (Villarini and Vecchi 2012a). In indicating that North Atlantic tropical cyclone frequency and intensity can behave disparately, these results using the CMIP5 GCMs are qualitatively consistent with previous analyses using models from phase 3 of the Coupled Model Intercomparison Project (CMIP3) (Emanuel et al. 2008; Bender et al. 2010).
Multiple forcing agents (GHG, aerosols, ozone, etc.) are changing in these RCPs. An idealized suite of experiments in which CO2 concentrations are doubled over a 70-yr period isolates the impact of CO2 in the projections (Fig. 3). An increase in CO2 results in an average slight increase in PDI, driven by the overall warming of the tropics (Villarini and Vecchi 2012a), in contrast to a CO2-driven decrease found for tropical storm frequency (Villarini and Vecchi 2012b). However the magnitude of the PDI sensitivity to CO2 in these models is not large enough to explain the increase in PDI for the three RCPs (e.g., compare the response at year 70 of Fig. 3 with 2100 in RCP 8.5 of Fig. 1), indicating that other forcing agents also contribute to the projected intensity increases.
To assess the role of aerosols in the PDI projections, we explore a partial perturbation experiment using a GCM from the CMIP5 suite [Geophysical Fluid Dynamics Laboratory Climate Model, version 3 (GFDL CM3)], in which aerosol precursors in RCP 4.5 are not allowed to change after 2005 and compare that with the response of the same model to the full RCP 4.5 projected forcing (Fig. 4). This set of experiments indicates that the projected increase in PDI in the projections from GFDL CM3 is driven both by GHG increases and aerosol changes (largely through a projected decrease in Atlantic and global aerosols; Villarini and Vecchi 2012b). The aerosol-driven increase in North Atlantic TC intensity in the GFDL CM3 projections arises both from the tropical-mean warming driven by aerosol optical depth decreases and from a warming of the Atlantic relative to the tropics driven by a more rapid reduction of aerosols over the Atlantic than over the global tropics in this RCP scenario. Unfortunately, this experiment is not available for the full CMIP5 model suite; because of the large potential role of aerosols in climate change, future coordinated experiments should include idealized experiments like these to allow an assessment of the relative contributions of aerosols and GHGs to projected climate changes.
An analysis of a couplet of historical CMIP5 experiments over the 1880–2005 period, one using “all forcings” (changing GHGs, aerosols, natural forcing, etc.) and another using past GHG forcing only (in which aerosols, natural forcings, etc. are kept at preindustrial values) indicates that the projected influence of aerosol changes on PDI may have begun in the 1990s (Fig. 5). These historical perturbation experiments are only available for a subset of the CMIP5 models (Table 1). In particular, the impact of GHGs alone leads to an increase in PDI over the twentieth century in these models, while the non-GHG forcing leads to a decrease between the 1960s and 1980s, and a rebound following that. Unfortunately again, experiments isolating the role of aerosols are available from few CMIP5 models at this time; however, the timing of the non-GHG decrease suggests that increasing aerosol loading in the Atlantic was a key countervailing force against a GHG-induced increase of PDI over the twentieth century, but over the recent decades (and in the projections of the twenty-first century—see Fig. 4) GHG and aerosol forcing both act to increase PDI in these models. In particular, the timing of the observed local minimum in PDI over the 1970s and 1980s is well captured by the radiatively forced models (comparing the dashed line with the “all forcing” CMIP5 line in Fig. 5), suggesting that this minimum was due, at least in part, to past changes in radiative forcing—likely anthropogenic aerosols. So the recent (since the mid-1980s) radiatively forced increase in our PDI index in the CMIP5 models is largely a rebound from a radiatively forced decrease over the 1960s and 1970s, likely because of anthropogenic aerosols, and is not representative of the sensitivity of PDI to GHG forcing. The amplitude of the observed changes in PDI is larger than that of the multimodel ensemble; we speculate that this is in part because the PDI is based on a statistical model that by construction will only capture part of the variance, and in part it is an indication that the observed changes included a nonradiative component. However, we cannot exclude the possibility that it may also indicate that this statistical downscaling of GCMs may underestimate the magnitude of forced changes in PDI, either because of problems with the statistical model, the GCMs, or the radiative forcing. Future work should explore these issues.
4. Discussion and conclusions
In this study, we used output from a new suite of coupled climate simulations (CMIP5; Taylor et al. 2012) and a recently developed statistical model (Villarini and Vecchi 2012a) to project changes in North Atlantic PDI over the twenty-first century. These analyses are based on 17 GCMs and explore three future radiative forcing scenarios (or RCPs). Comparison of the PDI projections with projections of North Atlantic TC frequency (Villarini and Vecchi 2012b) allows us to interpret the changes in terms of tropical cyclone intensity. Our results suggest that the North Atlantic PDI, driven primarily by changes to tropical cyclone intensity and the duration of TCs at highest intensity, is projected to increase over the current century in all three RCPs. By the end of the twenty-first century, the magnitude of the projected increase depends on the projected GHG forcing. The projected intensification of North Atlantic tropical cyclone is in response to GHG increases and aerosol changes.
The results of this study are based on the statistical model described in Villarini and Vecchi (2012a), in which the dominance of certain physical processes is implicit. In particular, the model assumes that tropical tropospheric warming will follow something close to a “moist adiabat” (warming will be about twice as large in the upper troposphere than at the surface) and that the “weak temperature gradient” (WTG) approximation holds (Sobel et al. 2002; the WTG approximation reflects the tendency of tropospheric temperature anomalies to be relatively spatially homogeneous in the tropics). Therefore, we do not account for the impact of direct radiative heating on free-atmospheric temperatures or other factors that would make tropospheric temperature changes deviate from a moist adiabat (e.g., Emanuel 2010; Emanuel et al. 2013; Vecchi et al. 2013), which would require other relevant predictors that are currently not included in the model. With these caveats in mind, these results point to a substantial increase in North Atlantic tropical cyclone intensity and growing probability of extreme hurricane seasons over this century (Fig. 6).
The CMIP5 coupled GCM experiment suite leads to a projection of increases in North Atlantic PDI over the twenty-first century in response to projected increases in GHGs and changes in atmospheric aerosols (largely reductions in Atlantic aerosol loading). The projections for increased PDI reflect a projection of increase tropical cyclone intensity and duration at the highest intensities, rather than an increase in frequency. These projected changes in PDI are large, indicating substantially increased probability of years as or more active than 2005 (Fig. 6), which may have been the most active year on record. However, these same models do not indicate that we should have seen an increase over the past century—nor do reconstructions of PDI from SST (Fig. 1; Villarini and Vecchi 2012a). The lack of an expectation for increasing PDI over the past century in these GCMs appears to arise in part because of large internal variability (Fig. 1), but also because the slight GHG-driven increase over the past century has been masked by a sharp non-GHG-driven decrease around the 1960s–80s—the timing of which suggests a role for aerosols (Fig. 5), similar to a potential masking of a GHG-induced weakening of oceanic circulation (e.g., Delworth and Dixon 2006). Only in the recent decades has the GHG and non-GHG PDI response in these models been in the same direction (Fig. 5), a constructive influence that is projected to continue over the next few decades (Fig. 4).
These results add to the growing body of work (e.g., Rotstayn and Lohmann 2002; Mann and Emanuel 2006; Evan et al. 2009; Chang et al. 2011; Booth et al. 2012; Villarini and Vecchi 2012b) suggesting that the observed multidecadal variability in the North Atlantic, and its related impacts (e.g., the inactive Atlantic hurricane period between the late 1960s and early 1990s), may include a component driven by changes in atmospheric aerosols (primarily through an increase in the post–World War II era, and a decrease in the 1980s–90s, of aerosol optical depth over the Atlantic). However, the CMIP5 historical experiments only explain a fraction of the recently observed multidecadal swing in PDI (Fig. 1), indicating that factors such as internal variability (e.g., Zhang and Delworth 2005, 2006, 2009; Robson et al. 2012) may also have contributed.
In contrast to projections of surface warming, which have already been observed and attributed in part to increasing GHGs (Solomon et al. 2007), for PDI we are in an uncomfortable position where the GCMs are projecting potentially dramatic and societally relevant changes, while at the same time indicating no detectable changes should be present in the record. Therefore, tests of these projections must be indirect and are intimately tied to our confidence in the fidelity of the GCMs and the projected radiative forcing. In particular, because of the role of aerosol changes in the historical simulations and projections of PDI, and because there are currently substantial uncertainties in the role of aerosols in past climate variations (e.g., Booth et al. 2012; Zhang et al. 2013), efforts should continue to improve our understanding and modeling capability for the role of aerosols in regional and global climate change. More generally, the mechanisms behind patterns of SST must be better understood (e.g., Leloup and Clement 2009; Clement et al. 2010; Xie et al. 2010), as should the character of past changes in regional SST (e.g., Vecchi and Soden 2007b; Vecchi et al. 2008a; Deser et al. 2010), to develop confident projections and assessment of past causes for changes in Atlantic hurricane intensity.
Acknowledgments
Gabriele Villarini acknowledges financial support by the Iowa Flood Center, IIHR-Hydroscience and Engineering. We thank Tom Knutson and Xiasong Yang for comments and suggestions and Larry Horowitz for access to the GFDL CM3 perturbation experiments. We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 1) 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 thank Dr. Stasinopoulos, Dr. Rigby, and Dr. Akantziliotou for making the GAMLSS (Stasinopoulos et al. 2007) package freely available in R (R Development Core Team 2008).
REFERENCES
Bender, M. A., T. R. Knutson, R. E. Tuleya, J. J. Sirutis, G. A. Vecchi, S. T. Garner, and I. M. Held, 2010: Model impact of anthropogenic warming on the frequency of intense Atlantic hurricanes. Science, 327, 454–458.
Bister, M., and K. A. Emanuel, 1998: Dissipative heating and hurricane intensity. Meteor. Atmos. Phys., 65, 233–240.
Booth, B. B., N. J. Dunstone, P. R. Halloran, T. Andrews, and N. Bellouin, 2012: Aerosols implicated as a prime driver of twentieth-century North Atlantic climate variability. Nature, 484, 228–232.
Chang, C.-Y., J. C. H. Chiang, M. F. Wehner, A. Friedman, and R. Ruedy, 2011: Sulfate aerosol control of tropical Atlantic climate over the twentieth century. J. Climate, 24, 2540–2555.
Clement, A. C., A. C. Baker, and J. Leloup, 2010: Climate change: Patterns of tropical warming. Nat. Geosci., 3, 8–9, doi:10.1038/ngeo728.
Delworth, T. L., and K. W. Dixon, 2006: Have anthropogenic aerosols delayed a greenhouse gas-induced weakening of the North Atlantic thermohaline circulation? Geophys. Res. Lett., 33, L02606, doi:10.1029/2005GL024980.
Deser, C., A. S. Phillips, and M. A. Alexander, 2010: Twentieth century tropical sea surface temperature trends revisited. Geophys. Res. Lett., 37, L10701, doi:10.1029/2010GL043321.
Emanuel, K. A., 1987: The dependence of hurricane intensity on climate. Nature, 326, 483–485.
Emanuel, K. A., 2005: Increasing destructiveness of tropical cyclones over the past 30 years. Nature, 436, 686–688.
Emanuel, K. A., 2007: Environmental factors affecting tropical cyclone power dissipation. J. Climate, 20, 5497–5509.
Emanuel, K. A., 2010: Tropical cyclone activity downscaled from NOAA-CIRES Reanalysis, 1908-1958. J. Adv. Model. Earth Syst., 2, 1–12.
Emanuel, K. A., R. Sundararajan, and J. Williams, 2008: Hurricanes and global warming—Results from downscaling IPCC AR4 simulations. Bull. Amer. Meteor. Soc., 89, 347–367.
Emanuel, K. A., S. Solomon, D. Folini, S. Davis, and C. Cagnazzo, 2013: Influence of tropical tropopause layer cooling on Atlantic hurricane activity. J. Climate, 26, 2288–2301.
Evan, A. T., D. J. Vimont, A. K. Heidinger, J. P. Kossin, and R. Bennartz, 2009: The role of aerosols in the evolution of tropical North Atlantic Ocean temperature anomalies. Science, 324, 778–781.
Holland, G. J., 1997: The maximum potential intensity of tropical cyclones. J. Atmos. Sci., 54, 2519–2541.
Knutson, T. R., and R. E. Tuleya, 2004: Impact of CO2-induced warming on simulated hurricane intensity and precipitation: Sensitivity to the choice of climate model and convective parameterization. J. Climate, 17, 3477–3495.
Knutson, T. R., and Coauthors, 2010: Tropical cyclones and climate change. Nat. Geosci., 3, 157–163.
Landsea, C. W., 1993: A climatology of intense (or major) Atlantic hurricanes. Mon. Wea. Rev., 121, 1703–1713.
Latif, M., N. Keenlyside, and J. Bader, 2007: Tropical sea surface temperature, vertical wind shear, and hurricane development. Geophys. Res. Lett., 34, L01710, doi:10.1029/2006GL027969.
Leloup, J., and A. Clement, 2009: Why is there a minimum in projected warming in the tropical North Atlantic Ocean? Geophys. Res. Lett., 36, L14802, doi:10.1029/2009GL038609.
Mann, M. E., and K. A. Emanuel, 2006: Atlantic hurricane trends linked to climate change. Eos, Trans. Amer. Geophys. Union, 87, 233–241, doi:10.1029/2006EO240001.
Mendelsohn, R., K. Emanuel, S. Chonabayashi, and L. Bakkensen, 2012: The impact of climate change on global tropical cyclone damage. Nat. Climate Change, 2, 205–209.
Oouchi, K., J. Yoshimura, H. Yoshimura, R. Mizuta, S. Kusunoki, and A. Noda, 2006: Tropical cyclone climatology in a global warming climate as simulated in a 20-km-mesh global atmospheric model: Frequency and wind intensity analysis. J. Meteor. Soc. Japan, 84, 259–276.
Peduzzi, P., B. Chatenou, H. Dao, A. De Bono, C. Herold, J. Kossin, F. Mouton, and O. Nordbeck, 2012: Global trends in tropical cyclone risk. Nat. Climate Change, 2, 289–294.
Pielke, R. A., J. Gratz, C. W. Landsea, D. Collins, M. A. Saunders, and R. Musulin, 2008: Normalized hurricane damages in the United States: 1900–2005. Nat. Hazards Rev., 9, 29–42.
R Development Core Team, cited 2008: R: A language and environment for statistical computing. R Foundation for Statistical Computing. [Available online at http://www.R-project.org.]
Ramsay, H. A., and A. H. Sobel, 2011: Effects of relative and absolute sea surface temperature on tropical cyclone potential intensity using a single-column model. J. Climate, 24, 183–193.
Robson, J., R. Sutton, K. Lohmann, D. Smith, and M. Palmer, 2012: Causes of the rapid warming of the North Atlantic Ocean in the mid-1990s. J. Climate, 25, 4116–4134.
Rotstayn, L. D., and U. Lohmann, 2002: Tropical rainfall trends and the indirect aerosol effect. J. Climate, 15, 2103–2116.
Smith, T. M., R. W. Reynolds, T. C. Peterson, and J. Lawrimore, 2008: Improvements to NOAA’s historical merged land–ocean surface temperature analysis (1880–2006). J. Climate, 21, 2283–2296.
Sobel, A. H., I. M. Held, and C. S. Bretherton, 2002: The ENSO signal in tropical tropospheric temperature. J. Climate, 15, 2702–2706.
Solomon, S., D. Qin, M. Manning, M. Marquis, K. Averyt, M. M. B. Tignor, H. L. Miller Jr., and Z. Chen, Eds., 2007: Climate Change 2007: The Physical Science Basis. Cambridge University Press, 996 pp.
Stasinopoulos, D. M., R. A. Rigby, and C. Akantziliotou, cited 2007: GAMLSS: Generalized Additive Models for Location, Scale and Shape. R Package Version 1.6-0. [Available online at http://cran.r-project.org/web/packages/gamlss/index.html.]
Swanson, K. L., 2008: Nonlocality of Atlantic tropical cyclone intensities. Geochem. Geophys. Geosyst., 9, Q04V01, doi:10.1029/2007GC001844.
Taylor, K. E., R. J. Stouffer, and G. A. Meehl, 2012: An overview of CMIP5 and the experiment design. Bull. Amer. Meteor. Soc., 93, 485–498.
Vecchi, G. A., and B. J. Soden, 2007a: Effect of remote sea surface temperature change on tropical cyclone potential intensity. Nature, 450, 1066–1071.
Vecchi, G. A., and B. J. Soden, 2007b: Global warming and the weakening of the tropical circulation. J. Climate, 20, 4316–4340.
Vecchi, G. A., A. Clement, and B. J. Soden, 2008a: Examining the tropical Pacific’s response to global warming. Eos, Trans. Amer. Geophys. Union, 89, 81–83, doi:10.1029/2008EO090002.
Vecchi, G. A., K. L. Swanson, and B. J. Soden, 2008b: Whither hurricane activity? Science, 322, 687–689.
Vecchi, G. A., M. Zhao, H. Wang, G. Villarini, A. Rosati, A. Kumar, I. M. Held, and R. Gudgel, 2011: Statistical–dynamical predictions of seasonal North Atlantic hurricane activity. Mon. Wea. Rev., 139, 1070–1082.
Vecchi, G. A., S. Fueglistaler, I. M. Held, T. R. Knutson, and M. Zhao, 2013:Impacts of atmospheric temperature trends on tropical cyclone activity. J. Climate, in press.
Villarini, G., and G. A. Vecchi, 2012a: North Atlantic power dissipation index (PDI) and accumulated cyclone energy (ACE): Statistical modeling and sensitivity to sea surface temperature changes. J. Climate, 25, 625–637.
Villarini, G., and G. A. Vecchi, 2012b: Twenty-first-century projections of North Atlantic tropical storms from CMIP5 models. Nat. Climate Change, 2, 604–607.
Villarini, G., and G. A. Vecchi, 2013: Multiseason lead forecast of the North Atlantic power dissipation index (PDI) and accumulated cyclone energy (ACE). J. Climate, in press.
Villarini, G., G. A. Vecchi, and J. A. Smith, 2010: Modeling of the dependence of tropical storm counts in the North Atlantic basin on climate indices. Mon. Wea. Rev., 138, 2681–2705.
Villarini, G., G. A. Vecchi, T. R. Knutson, M. Zhao, and J. A. Smith, 2011: North Atlantic tropical storm frequency response to anthropogenic forcing: Projections and sources of uncertainty. J. Climate, 24, 3224–3238.
Villarini, G., G. A. Vecchi, and J. A. Smith, 2012: U.S. landfalling and North Atlantic hurricanes: Statistical modeling of their frequencies and ratios. Mon. Wea. Rev., 140, 44–65.
Xie, S. P., C. Deser, G. A. Vecchi, J. Ma, H. Teng, and A. T. Wittenberg, 2010: Global warming pattern formation: Sea surface temperature and rainfall. J. Climate, 23, 966–986.
Zhang, R., and T. L. Delworth, 2005: Simulated tropical response to a substantial weakening of the Atlantic thermohaline circulation. J. Climate, 18, 1853–1860.
Zhang, R., and T. L. Delworth, 2006: Impact of Atlantic multidecadal oscillations on India/Sahel rainfall and Atlantic hurricanes. Geophys. Res. Lett., 33, L17712, doi:10.1029/2006GL026267.
Zhang, R., and T. L. Delworth, 2009: A new method for attributing climate variations over the Atlantic hurricane basin’s main development region. Geophys. Res. Lett., 36, L06701, doi:10.1029/2009GL037260.
Zhang, R., and Coauthors, 2013: Have aerosols caused the observed Atlantic multidecadal variability? J. Atmos. Sci., 70, 1135–1144.
Zhao, M., I. M. Held, and G. A. Vecchi, 2010: Retrospective forecasts of the hurricane season using a global atmospheric model assuming persistence of SST anomalies. Mon. Wea. Rev., 138, 3858–3868.