Simulated reduction in Atlantic hurricane

Using projected boundary conditions for the end of the twenty-first century, the frequency of Atlantic tropical cyclones and hurricanes in a regional climate model of the Atlantic basin is reduced compared with observed boundary conditions at the end of the twentieth century. This is inconsistent with the idea that higher levels of atmospheric greenhouse gases will result in increased Atlantic hurricane activity. Increasing sea surface temperatures in the tropical Atlantic Ocean and measures of Atlantic hurricane activity have been reported to be strongly correlated since at least 1950 (refs 1, 2, 3, 4, 5), raising concerns that future greenhouse-gas-induced warming6 could lead to pronounced increases in hurricane activity. Models that explicitly simulate hurricanes are needed to study the influence of warming ocean temperatures on Atlantic hurricane activity, complementing empirical approaches. Our regional climate model of the Atlantic basin reproduces the observed rise in hurricane counts between 1980 and 2006, along with much of the interannual variability, when forced with observed sea surface temperatures and atmospheric conditions7. Here we assess, in our model system7, the changes in large-scale climate that are projected to occur by the end of the twenty-first century by an ensemble of global climate models8, and find that Atlantic hurricane and tropical storm frequencies are reduced. At the same time, near-storm rainfall rates increase substantially. Our results do not support the notion of large increasing trends in either tropical storm or hurricane frequency driven by increases in atmospheric greenhouse-gas concentrations.

For the present-day simulations, the model uses observed SSTs, and the interior solution is nudged towards observed atmospheric reanalyes 16 on large spatial scales (see the Methods section) for each August-October season from 1980 to 2006. Hurricanes develop from weak disturbances that are generated spontaneously by the model, in ways that are dependent on the large-scale environment. To demonstrate the model's hurricane simulation skill, Fig. 1a shows simulations obtained using two slightly different versions of the model (referred to here as Model1 and Model2) 7 . Model2 is our control for the climate warming experiments described below. The year-to-year variability in Atlantic hurricane counts  is well reproduced by the model ensemble mean (r = 0.84), as well as by Model1 and Model2 individually 7 . El Niño-related interannual variations in hurricane activity are also well reproduced 7 . The model simulates an increase in hurricane counts , although the simulated trend is ∼40% larger than the observed increase (Fig. 1a).
A limitation of the model is that it does not simulate hurricanes as intense as observed (Fig. 1b, Supplementary Information,  Fig. S1). Minimum central pressures of ∼937 hPa and maximum surface winds of ∼47 m s −1 are the most intense simulated in the control, compared with 882 hPa and ∼85 m s −1 maximum intensities in observations. Thus, model simulations of indices that are strongly influenced by intense hurricanes (such as PDI) should be treated with caution. The ensemble model produces an increasing trend in PDI   7 , although it is ∼25% smaller than the observed increase. For major hurricanes (using the central pressure criterion of <965 hPa (ref. 7)), both the mean number and increasing trend  in the model are about half as large as observed, although the model's variability correlates well (r = 0.7) with observations 7 .
Storm trajectories for contrasting active and inactive seasons ( Fig. 2a-d) further illustrate the control model's performance. The truncated tracks near the northern boundary of the model domain are an artefact of the boundary nudging. Animated fields from sample seasons are available in the Supplementary Information.
We use Model2 to explore the influence of future climate warming on Atlantic hurricane activity. We re-run each August-October season , keeping the daily to multidecadal variations unchanged, but altering the August-October mean atmospheric state (used by the interior nudging) and SSTs according to late twenty-first-century changes   Fig. S2a), enhanced upper-tropospheric warming relative to the surface (see Supplementary Information, Fig. S2b) and enhanced vertical wind shear and reduced low-to mid-tropospheric relative humidity across the Caribbean 17 .
In the warm-climate runs, our model generates substantially fewer tropical storms (−27%) and hurricanes (−18%) (Fig. 3, Supplementary Information, Fig. S3). Comparing 27 warm-climate seasons with the control, 22 years have a decrease in tropical storms, 3 have an increase and 2 show no change. Major hurricane counts decrease by 8%. Thus, the reduction in numbers is most pronounced for weaker storms. Storm tracks for the warming and control climates, as well as observed, are shown in Fig. 2 and Supplementary Information, Fig. S4 (the latter for each individual year). Dynamical reasons for the reduced storm frequency remain under investigation. Preliminary sensitivity experiments indicate that the projected combined change in tropical SST and atmospheric lapse rate is not the primary driver.
Environmental changes in circulation and/or moisture are more likely the dominant factors in reducing storm frequency.
In agreement with previous studies [18][19][20][21] , the frequency and intensity of the strongest hurricanes simulated by the model are increased for both surface wind speed (Fig. 1b,c) and central pressure (see Supplementary Information, Fig. S5). More than twice as many hurricanes occur with wind speeds exceeding 45 m s −1 in the warm-climate runs (twelve versus five). The normalized frequency distributions ( Fig. 1c) show, for the warm climate, an increased relative occurrence rate for all storm intensity categories above 35 m s −1 . The mean maximum wind speed increases slightly, (+2.9%) for tropical storms and hurricanes combined and by 1.7% for hurricanes alone. However, simulated PDI decreases by about 25% (see Supplementary  Information, Fig. S3f), owing to the decrease in tropical storm and hurricane frequency. As both the frequency of major hurricanes (see Supplementary Information, Fig. S3g) and the maximum intensities of the strongest hurricanes ( Fig. 1b) are under-predicted in the control model 7 , results related to the most intense simulated storms should be regarded as only suggestive.
Near-hurricane rainfall rates (see Supplementary  centre. This is in reasonable quantitative agreement with a previous hurricane modelling study 22 . A caveat is that there is some evidence (see Supplementary Information, Fig. S6) that the control model's composite hurricane precipitation rates are higher than observed 23 . As in previous studies, there is no indication that the region of Atlantic tropical storm formation expands with greenhouse warming (see Supplementary Information, Fig. S7). Storms typically develop in the control simulation, and in observations, over SSTs in excess of a 'threshold' temperature of about 26 • C. In the warmer climate, given the relatively homogeneous tropical warming, the 'threshold' temperature for storm formation increases in the model by the same amount as the mean tropical SSTs (Fig. 4d,e). Consequently, no pronounced expansion of the tropical storm formation region is simulated. This 'threshold' temperature for cyclone formation is clearly climate dependent, and is closely related to the temperature needed for a saturated parcel near the surface to be sufficiently buoyant to rise to the tropical tropopause.
The decrease in hurricane and tropical storm occurrence is slightly more evident in the western than eastern half of the basin (see Supplementary Information, Fig. S7). This is qualitatively similar to occurrence trend patterns since the late 1800s in a recent tropical storm reconstruction 15 . Consequently, the simulated number of landfalling hurricanes decreases by more than the hurricane count itself (about 30% versus 18%; Supplementary  Information, Fig. S3). We note that the band of increased vertical shear and decreased mid-tropospheric relative humidity stretching across the Caribbean in the multimodel ensemble projections 17 should suppress activity more in the western half of the basin than in the east. As the detailed spatial pattern of shear and humidity change differs among the CMIP3 models 17 , this result may not be robust.
Uncertainties in the large-scale climate projections will translate into uncertainties in the downscaled tropical storm projections in ways that have not yet been quantified. Here, we use only the timemean changes projected by the CMIP3 climate models. We assume that changes in the time-mean state are the dominant control over tropical cyclone formation and that the CMIP3 ensemble mean provides the best available projection of these changes. Any simulated changes in atmospheric variability from the CMIP3 models are neglected, but our confidence in such changes obtained using current models is low.
Further simulations are underway with our downscaling model using mean climate changes generated by individual CMIP3 models rather than the ensemble mean across the models. Preliminary results suggest that the spread in hurricane activity projections obtained using individual models will be substantial (as also suggested in a recent study using a different modelling strategy 24 ). As an initial example, we used the mean climate change from one model (GFDL CM2.1) that we expected would be more favourable than the ensemble mean for increased hurricane activity, on the basis of changes in large-scale shear over the Atlantic. A sample of 13 simulated seasons shows relatively little change in tropical storm (−8%) or hurricane counts (+7%), but a relatively large increase (∼+70%) in major hurricane counts. This sensitivity of the hurricane response to details of the climate model projections highlights the need to better constrain regional climate responses to increased CO 2 .
A key limitation of our regional downscaling model is that it does not simulate the most intense hurricanes or the strong dependence of upper-limit intensities on SST observed in the present climate (see Supplementary Information, Fig. S1). Owing LETTERS to these limitations, it is plausible that the model's quantitative projections of increased intensity and increased numbers of the most intense storms with climate warming are underestimates. Follow-on studies are underway with models capable of simulating higher-intensity hurricanes.
Our results using the ensemble-mean global model projections (Fig. 4) are inconsistent with the notion of large, upward trends in tropical storm and hurricane frequency over the twentieth century, driven by greenhouse warming 3,5 . The twentieth-century warming attributed to greenhouse gases by global models is sufficiently similar in spatial structure to the twenty-first-century projections used here that we are confident that our model's downscaled storm response to twentieth-century greenhouse warming would also be qualitatively similar to the results presented here. However, the general issue of twentieth-century attribution is made more complex by other radiative forcing changes (for example, aerosols, stratospheric ozone) that have no direct counterparts in the twentyfirst-century A1B scenario.
These results do not conflict with the hypothesis that low-frequency variability of Atlantic hurricane activity is largely controlled by SSTs, but they do disagree with the idea that these variations are primarily determined by local tropical Atlantic SSTs alone. When the North Atlantic warms more rapidly than other ocean basins, as it has since 1980, the changing SST gradients favour northward displacement of the Atlantic Intertropical Convergence Zone, and-in the main development region for Atlantic storms-reduced vertical shear 25,26 , reduced stability of the tropospheric thermodynamic profile 27,28 and increased potential intensity 29 . Although the relative importance of these factors is unclear, they all favour tropical storm genesis, and probably produce the positive correlation between Atlantic SST and hurricane frequency in both observations and our model (Fig. 4a,b). The observed correlation between the Atlantic basin PDI and the difference between local tropical Atlantic SSTs and tropical mean SSTs 30 also supports this interpretation. However, during the twenty-first century, tropical North Atlantic SSTs are not projected to warm more rapidly than other tropical regions according to the ensemble of CMIP3 models 29 , implying conditions that do not favour increased tropical storm genesis.
In summary, our model supports the hypothesis that the primary driver of the recent increase in Atlantic hurricane numbers was the warming of the tropical Atlantic relative to the other tropical basins, rather than local increases in tropical Atlantic SSTs alone. The decrease in Atlantic tropical storm and hurricane frequency simulated by our model-in response to a pronounced twenty-first-century greenhouse warming-is markedly different from the large increase in future hurricane activity that might be projected using statistical analysis based only on the recent covariability of hurricanes and local SSTs.

METHODS
This study uses a regional atmospheric modelling framework designed specifically for dynamically downscaling Atlantic hurricane activity 7 . The non-hydrostatic model has a grid spacing of 18 km and is run without a subgrid-scale moist convective parameterization. The model is forced at the lower boundary by observed SSTs, and is relaxed at its horizontal boundaries, on a timescale of 3 h, to agree with the time-evolving atmospheric state estimated by the National Centers for Environmental Prediction/National Center for Atmospheric Research reanalysis 16 . We also Fourier analyse the time-evolving reanalysis in the domain interior and relax the model (36 h timescale) on the largest scales (basin wavenumbers 0-2) to these observed atmospheric conditions, including the winds, temperatures and relative humidity at all levels. The strength of the interior large-scale relaxation is the sole explicit tuning parameter in the model. This tuning, along with other details of the model design, is described in ref. 7. This parameter has the same value in every seasonal simulation and in both the control and warm-climate integrations.
Our evaluations of the model 7 thus far have focused mostly on hurricane metrics, including hurricane precipitation. Many other metrics, such as the general space-time characteristics of the model's tropical convection, apart from seasonal-mean precipitation 7 , have not been evaluated in detail.
For our climate change experiment, we re-run each August-October season , relaxing towards a new solution consisting of the original observed conditions altered by a time-invariant three-dimensional climate change perturbation. This procedure keeps unchanged the largescale interannual to multidecadal variations in the interior as well as the high-frequency weather variability imposed at the model boundaries. The time-invariant climate change is the average change simulated by a multimodel ensemble of CMIP3 models 8 used in the IPCC Fourth Assessment Report 6 (years 2080-2099, A1B emission scenario). This method avoids some difficulties in the direct use of climate model simulations, which have known biases that can distort tropical storm simulations. We use the multimodel ensemble rather than separately downscaling each climate model, the latter being computationally prohibitive. The method assumes that the average downscaled hurricane response of the models can be approximated through downscaling the average climate change from the models. It further assumes that the atmospheric variability in the large (interior) scales and boundary conditions of our control simulation is also representative of conditions under the global warming scenario.
Using a multimodel ensemble to estimate the climate change perturbation also helps distinguish the models' radiatively forced response from model internal variability, which can be particularly pronounced in the Atlantic basin. In summary, our method simultaneously reduces possible contamination by internal climate variability, while also leveraging the benefits of ensemble averaging across climate models with different physics, tunings and parameterizations.
An important research question is whether hurricanes are able to maintain their intensities longer as they travel into higher latitudes in a warmer climate. Unfortunately, the placement of the model's northern boundary prevents us from addressing this issue.
In addition, as our simulations are restricted to the core of the Atlantic hurricane season (August-October), we do not address whether climate warming will change the season length by altering pre-August or post-October storms.
Although we originally intended to use the 18-model ensemble described in ref. 17 as our multimodel ensemble, owing to technical problems, we used a 14-model average to form the relative humidity perturbation, a 17-model average for atmospheric temperatures and winds and all 18 models of ref. 17 for SSTs. The difference between the full 18-model ensemble and the ensemble mean fields we used in our climate change runs is a very small fraction of the standard deviation of the 18-model estimates. It is extremely unlikely that a set of runs using means over all 18 models could be statistically distinguished from those reported here. The primary source of uncertainty in our climate perturbations is the uncertainty in the climate model projections themselves, rather than the precise set of models used to define our multimodel ensemble climate perturbation.