Tropical cyclone genesis across palaeoclimates

Introduction Conclusions References

. However, recent analyses of future simulations performed as part of the Coupled Model Intercomparison Project Phase 5 (CMIP5) appear equivocal: statistical downscaling indicates an increase in both cyclone intensity and genesis (Emanuel, 2013); dynamical downscaling indicates an increase in intensity combined with a reduction in frequency (Knutson et al., 2013); tracking algorithms of global coupled models do likewise (Camargo, 2013); large-scale cyclogenesis indices have shown both frequency increases (Emanuel, 2013) and decreases (Camargo, 2013).
Understanding past climates provides a means for scientists to contextualise future climate change impacts.Palaeoclimates with altered climate forcings, such as the elevated levels of carbon dioxide during the Pliocene period, may provide clues about how the trend of cyclone genesis would respond to ongoing anthropogenic emissions of greenhouse gases.
The mid-Piacenzian warm portion of the mid-Pliocene (around 3 million years ago, henceforth "Pliocene") was a recent episode in Earth's geological history where mean global temperatures were warmer by 2-3 • C compared to modern times (Haywood et al., 2013), but the warming was not constant across the globe.Sea surface temperatures (SST) anomalies were more pronounced at the higher latitudes (up to 20 • C in the high Arctic; Ballantyne et al., 2009), while the lower latitudes exhibited minimal change in places (Dowsett et al., 2010).The geography of the continents and oceans were relatively similar to earth's current configuration (Haywood et al., 2011).Carbon dioxide concentrations were at or near present day levels during the mid-Pliocene (Pagani et al., 2009).There is potential for using the Pliocene to learn about the equilibrium state of the earth's warm climate following anthropogenic greenhouse gas influence (Haywood et al., 2009).Meanwhile, the icehouse climate of the Last Glacial Maximum (LGM) at 21ka serves as a contrast to the current greenhouse climate.Proxy estimates by Annan and Hargreaves (2013) suggest that LGM tropical SST was around 1.6 • C lower than preindustrial, while global surface air temperatures were 3.1-4.LGM and pre-industrial intervals, the Palaeoclimate Model Intercomparison Project (PMIP)'s focus on these times help facilitate studies that examine the effect of carbon dioxide concentration changes on the tropical climate (Table 1).
On the other hand, simulations for the mid-Holocene at 6ka differ from pre-industrial conditions mainly in the orbital parameters that result in an increased insolation in the high latitudes.The tropical region during the mid-Holocene interval might have encountered slightly elevated sea-surface temperatures (SST) of around 1 • C (Gagan et al., 1998), although recent studies indicate some uncertainty in terms of negative SST anomaly for regions such as the western Indian Ocean (Kuhnert et al., 2014).Despite the limited proxy record agreement on whether tropical oceans may have warmed (Koutavas et al., 2002;Rimbu et al., 2004;Stott et al., 2004), prior PMIP simulations suggest SST in the Northern Hemisphere was generally warmer by less than 1 • C in the mid-Holocene period compared to the pre-industrial era, and the Southern Hemisphere might have been slightly cooler (Braconnot et al., 2007).
Given the lack of data on tropical cyclone frequency for these palaeoclimates, model simulation studies do not seek to verify model response on cyclone formation, but rather aim to describe tropical cyclone trends with the assumption that signals would be detectable by using indicators such as cyclogenesis potential.Using PMIP Phase 2 (PMIP2) data, studies have been conducted to investigate TC genesis activity during the LGM and mid-Holocene intervals.For the LGM, Korty et al. (2012a)  ensemble (PlioMIP) is included to investigate whether there is a robust response to carbon dioxide concentrations.A further objective is to explore how TC activity in these palaeoclimates, with their equilibrium conditions, relates to that in future simulations.
The various model simulations in this study are described in Sect. 2. The calculation of genesis potential index (GPI) that underpins this study will be presented in Sect. 3 of this paper.Section 4 consolidates the results from the GPI analysis of the various palaeoclimates derived from the GCM ensembles.A discussion of how the climate during the Pliocene, LGM and mid-Holocene may affect TC genesis relative to the pre-industrial will be covered in Sect.5, and the effects of elevated carbon dioxide concentration on TC genesis will also be investigated.Section 6 will summarise the paper's key findings.

Climate simulations
The Pliocene Model Intercomparison Project (PlioMIP), which complements the LGM and the mid-Holocene aspects of the PMIP Phase 3 (PMIP 3), coordinates the efforts of various international climate modelling teams to quantify uncertainties in model outputs using the average interglacial conditions of the mid-Piacenzian (hereafter known as Pliocene) climate boundary conditions between 3.29 and 2.97 Ma (Haywood et al., 2011).
Nine coupled climate models participated in PlioMIP (Haywood et al., 2013), although only five are analysed here.The GCM selection for this study is largely dependent on data availability for the large-scale climatic variables such as the atmospheric temperature and humidity profile from the PlioMIP project for the Pliocene Epoch.PMIP3 data for the LGM, mid-Holocene and pre-industrial are taken from the same GCM that is used in the Pliocene simulation, subject to data availability.Where a different GCM from the same model family has to be selected, a pre-industrial control from that particular GCM will also be used for the comparisons.Data for the representative concentration pathway culminating in a radiative forcing of 8. 8.5) is likewise analysed as an example future elevated carbon dioxide concentration scenario.The GCMs that have been included for this study are outlined in Table 2.
In situations where a pre-computed monthly climatology of a particular simulation is not available on the Earth System Federation Grid, a 50 year time-slice from the end of the interval of interest is used to generate the monthly climatology data so as to minimise stochastic effects, model drift and internal variability.To estimate the ranges associated with internal variability, ten decadal time-slices are taken from a hundred year dataset of the pre-industrial dataset of each model.The standard deviation (SD) is found to be within 1-3 % of the pre-industrial (PI) TC genesis annual frequencies simulated across the five GCMs (Table 2).

Genesis potential index
The use of a metric related to genesis potential is particularly useful to describe cyclogenesis as the grid resolution of most GCMs is not sufficiently refined to simulate mesoscale processes taking place in TCs.Many studies have used genesis potential indices as a less computationally intensive approach to describe TC genesis (Bruyère et al., 2012;Camargo et al., 2007;Emanuel and Nolan, 2004;Korty et al., 2012a, b;Menkes et al., 2012;Tippett et al., 2010).Gray (1975) pioneered work on the genesis potential index (GPI) by demonstrating the use of selected diagnostics such as mid-troposphere humidity, vertical shear of the horizontal winds between the high and low level troposphere, low level relative vorticity, and thermal parameters related to SST to characterise climatic conditions that are favourable for cyclone genesis.The subsequent GPI improved by Emanuel and Nolan (2004) is considered state-of-the-art (Tippett et al., 2010) and incorporates the potential intensity theory (Emanuel, 1988;Holland, 1997) that evaluates the maximum wind speed that may be attainable using the available thermodynamic energy imparted from the atmospheric environment and the sea surface to the cyclone.The GPI proposed by Emanuel and Nolan (2004) serves to synergise the thermodynamic and kinematic factors affecting TC genesis into a single index.With the aim of facilitating comparison with previous investigations into palaeoclimate cyclone genesis, the "clipped vorticity" version of the GPI employed by Korty et al. (2012a, b) has likewise been adopted for this study: (1) Here, η represents the absolute vorticity computed at the 850 hPa level (Nolan and Rappin, 2008), V shear is the 200-850 hPa wind shear value, χ m is the moist entropy deficit.PI is the maximum potential intensity a TC can theoretically achieve (a potential intensity of less than 35 m s −1 would not be classified as a tropical storm).Due to the inherent biases in convection schemes and parameterisations employed by GCMs, the global annual total TC genesis has to be calibrated (Emanuel et al., 2008b).b is therefore an empirically derived normalisation factor that calibrates the GPI to achieve a pre-industrial control annual cyclone genesis frequencies on the order of the roughly eighty storms observed per year in the modern period.Following Korty et al. (2012a, b), a value of 5.55 × 10 −5 is used throughout for the value of b.
Wind shear and absolute vorticity are the two kinematic factors of the GPI, while potential intensity and moist entropy deficit are both thermodynamic factors (Korty et al., 2012a).Wind shear, which is the vertical shear of the horizontal winds between the upper and lower troposphere, causes asymmetries in the developing cyclone which results in the ventilation of the upper level warm core through the flushing of relatively cooler and drier air from the top (Frank and Ritchie, 2001).Stronger wind shear therefore influences inflow dynamics and weakens cyclone formation (Riemer et al., 2013).While noting caveats where such two-level vector differentials may be inadequate to describe the resultant wind shear in some scenarios (Velden and Sears, 2014) study defines the wind shear as the difference between the 200 and 850 hPa winds given its ease of computation.Meanwhile, the vorticity serves as a spin-up mechanism that initiates cyclone formation in a recirculating flow that is quasi-closed in the lower troposphere.Taking the analogy of a protective pouch, the quasi-closed streamlines surround the enhanced vorticity while nurturing the thermodynamic and convective processes that favour TC development (Tory et al., 2012).Tippett et al. (2010) observed that vorticity has a greater influence on cyclone formation at lower latitudes, and other factors play a greater role at higher latitudes.Korty et al. (2012a, b) subsequently incorporated a "clipped vorticity" diagnosis in place of absolute vorticity in the GPI, so as to moderate its response in over-estimating TC activity for the sub-tropics.Potentially, the clipping threshold may have varied in the past through large-scale changes in the atmosphere circulation.The sensitivity of different clipping thresholds was investigated in the case of the Pliocene period using the MIROC model.Changes in the clipping threshold appear to have little substantive impact on the resulting change in GPI for this study.
The non-dimensional term (χ m ) measures the moist entropy difference between the mid-troposphere and the boundary layer that is derived from asymmetric cyclone models (Emanuel, 1995b), as shown below: s m , s b and s * o represent the moist entropies of the mid-troposphere layer, boundary layer, and the sea surface saturation entropy respectively.Taking the assumption that the lapse rate of the tropical atmosphere is largely moist adiabatic (Emanuel et al., 2008b), s * which is the saturation entropy above the boundary layer, is assumed to be constant throughout the atmospheric column.This allows the numerator term in Eq. ( 2) to be evaluated at 600 hPa which is taken to represent the mid-troposphere as defined by Emanuel (1994).needed for an initial perturbation to moisten the middle troposphere before intensification occurs (Emanuel et al., 2008b).
Taking on the analogy of a cyclone's evolution process as equivalent to Nature's Carnot engine (Emanuel, 1988(Emanuel, , 1991)), the potential intensity diagnostic derived by Bister and Emanuel (1998Emanuel ( , 2002) ) that takes into account the effects of dissipative heating is: C k and C d are the surface exchange coefficients for enthalpy and momentum.Its ratio could range between 0.1 to 1.3 (Montgomery et al., 2010) and is likely between 0.75 and 1.5 for naturally occurring cyclones (Emanuel, 1995a).In this study, a ratio of 1 is taken to allow for ease of comparison with previous work that used a similar assumption (Korty et al., 2012a).T o is an entropy-weighted mean temperature of the outflow.The convective available potential energy (CAPE * ) describes an air parcel of maximum wind intensity that has been earlier saturated at the sea surface, while CAPE b describes a boundary layer air parcel which has been isothermally lowered from an equivalent air parcel of maximum wind intensity.Climate variables that are required for the potential intensity calculation include SST and pressure of the sea surface, as well as the humidity and temperature profile of the atmospheric column.The calculation of potential intensity for this study is facilitated by the use of an previously applied algorithm (Emanuel et al., 2008a).Figures

Back Close
Full

Potential intensity
In the tropical region, the Pliocene saw higher SSTs by about 2 • C relative to the preindustrial control (and the mid-Holocene), while SSTs were lower by about 2 • C at the LGM (Fig. 1).Korty et al. (2012a) suggested that high values of potential intensity, typically higher than 55 m s −1 , are needed to induce deep tropospheric convection in TC genesis.Interestingly the locations of the 55 m s −1 potential intensity contour appears to be relatively insensitive to these wholesale SST changes.For example, the contour in the North Pacific is associated with SSTs ranging from 26 • C during the Pliocene to 22 • C at the LGM.
During the Pliocene, there is a reduction in potential intensity over the North Atlantic, despite an SST increase in the same region (Fig. 1b).This supports research showing that absolute SST by itself can be an inadequate indicator of storm strength (Vecchi et al., 2008).Whilst, this may appear to depart from early understanding of threshold SST values (e.g. 26 • C) in influencing cyclone genesis (Palmen, 1948), it rather underscores the importance of other factors, such as atmospheric humidity and upper troposphere outflow temperature relative to the SST, that jointly determine the magnitude of energy available to a tropical cyclone (Emanuel, 1998).

Pre-industrial
The pre-industrial serves as a useful reference climate as it is before Earth's environment came under substantial anthropogenic influence, especially over the tropical oceans (Lewis and Maslin, 2014).period of January, February, March, April (JFMA).Monthly storm genesis will be discussed in Sect.4.6.
The GPI distribution of the various GCMs compares favourably with the outcomes from similar model analysis by Camargo (2013) for the pre-industrial period, despite the use of different genesis indices.All models simulate cyclone genesis in the eastern and western Pacific in the Northern Hemisphere during JASO, as well as the eastern Pacific near the South Pacific Convergence Zone (SPCZ) during JFMA.Stronger cyclone genesis in the southern Indian Ocean is found during JFMA, with limited genesis in the northern Indian Ocean during JASO apart from some areas such as the northern Bay of Bengal.The North Atlantic features some genesis potential at the deep and sub-tropics, but the South Atlantic shows almost negligible potential for TC genesis.These features are all shown in observations of tropical cyclone genesis (Knapp et al., 2010).
However the various models do show some biases.CCSM4, HadCM3 and IPSL-CM5A-LR exhibit greater intensity of cyclone genesis over the North Pacific and at the SPCZ; as compared to FGOALS-G2 and MIROC-ESM which appear to have a more extensive area of cyclone genesis but with lower intensity.While IPSL-CM5A-LR suggests that the central-western Pacific would experience the greatest intensity of cyclone genesis, CCSM4 and HadCM3 show the greatest genesis intensity arising in the northeastern Pacific.MIROC-ESM simulates its strongest genesis intensity in the maritime continent region.FGOALS-G2 shows a relatively uniform strength of genesis potential across all the oceans, apart from an area of increased intensity in the eastern North Pacific and Philippine Sea.The genesis potential also stretches across a greater area in FGOALS-G2 relative to the other models.
The ensemble mean (Fig. the highest intensity of genesis potential is located between the 10-20 • latitude belts of the respective peak storm periods of both hemispheres, and both are of comparable cumulative genesis magnitude of between 3-5 occurrences m −2 month −1 (not shown).
The pre-industrial climate thus exhibits consistency in cyclogenesis between the PMIP3 and PMIP2 simulations (only HadCM3 occurs in both ensembles).

Mid-Holocene
The key difference between the mid-Holocene and pre-industrial climate lies in the changes in solar insolation arising from different angular precession (Table 1).As a result, the Northern Hemisphere receives proportionally greater insolation during its storm season compared to the Southern Hemisphere.The summer and annual mean insolation for the high latitudes in both hemispheres is also increased (Braconnot et al., 2007).These insolation changes drive responses in diagnosed cyclogenesis across the five models (Fig. 3).CCSM4 shows the strongest response while the response from MIROC is relatively muted.HadCM3 exhibits a reduction of genesis potential in the Northern Hemisphere compensated for by an increase in the Southern Hemisphere.The response of FGOALS and CCSM bear similarities to each other in that their bands of TC genesis in the North Pacific become more zonal.IPSL and MIROC show little change from the pre-industrial.
The ensemble genesis potential for the mid-Holocene (Fig. 4a) shares largely similar distribution as the pre-industrial period (Fig. 2), although some localised differences in genesis strength are observed (Fig. 4b).The Southern Hemisphere exhibits a weak change from mid-Holocene over pre-industrial, except for pockets around Northern Australia that show a stronger increase in genesis magnitude.There is a good agreement across the ensemble on the sign of the mid-Holocene change in most areas amongst the five GCMs (Fig. 4c).There is a general decrease in cyclone genesis in the Northern Hemisphere, and an increase in cyclone genesis as one moves polewards in the Southern Hemisphere.Although several regions show strong agreement in increased cyclone genesis, such as the South Pacific and South Atlantic, the fact that the increased cyclone genesis are not depicted in Fig. 4b indicates only a very small degree of cyclone genesis increase in these areas.
The results for the mid-Holocene using these PMIP3 models compare well with findings from Korty et al. (2012b) that detail cyclone genesis using an ensemble from ten GCMs from PMIP2.The magnitude and distribution of genesis potential share similar patterns across all oceans.Nonetheless this study simulates a slightly weaker genesis potential for the western South Indian Ocean and the South Atlantic, as well as a slightly weaker increase in genesis potential for mid-Holocene over pre-industrial in both hemispheres.The model agreement (Fig. 4c) shares remarkable similarity with Korty et al. (2012b) on the sign of change in genesis potential of the mid-Holocene era compared to the pre-industrial period, where the Northern Hemisphere shows a decrease in genesis potential while the Southern Hemisphere shows general increase in genesis potential, especially in the higher tropical latitudes.The ensemble genesis potential for the LGM (Fig. 6a) shares again, at a first glance, a similar distribution with the pre-industrial.However, it exhibits greater intensity of genesis potential in the central North Pacific and near the SPCZ (Fig. 6b).The centraleastern South Indian Ocean shows decrease in genesis potential along 10

Last Glacial
• S, and slight decreases of genesis potential are observed in the North Atlantic, as well as the eastern and western North Pacific There is some model agreement (Fig. 6c) focussed around the largest changes in genesis potential in the LGM period for most oceans relative to pre-industrial.The North Atlantic exhibits a very robust decrease in genesis potential that spreads over Central America into the eastern North Pacific.This is likely a response to the imposition of the Laurentide ice sheet and its impact on the regional circulation.There appears to be a dipole pattern in the Indian Ocean (most noticeable in Fig. 6c), although it is not as robust.This is likely an expression of the alteration in Walker Circulation (DiNezio et al., 2011), whose fidelity varies across models depending on their parameterisations and boundary conditions (DiNezio and Tierney, 2013).These patterns of the model agreement are qualitatively similar to those seen in the PMIP2 experiments (Korty et al., 2012a).

Pliocene
The Pliocene is a warmer climate compared to pre-industrial (Dowsett et al., 2010;Haywood et al., 2013), with the area-averaged tropical SST from the five GCMs in this study over the peak storm season being 1.7 • C warmer.In terms of the tropical cyclogenesis difference from pre-industrial (Fig. 7 The key areas of Pliocene cyclone genesis in the ensemble mean can be found in the eastern and western areas of the North Pacific, the SPCZ and central region of the South Pacific, as well as the north-western corner of the South Indian Ocean (Fig. 8a).In terms of the difference in genesis potential between the Pliocene and pre-industrial periods (Fig. 8b), the North Atlantic, North Pacific, and South Indian oceans and the SPCZ region would largely experience a decline in genesis during the Pliocene.Nonetheless, there may be localised areas where genesis potential could have increased, such as to the east of New Guinea in the South Pacific.
This large-scale pattern appears to be robust as most models suggest a general decrease in genesis potential for the Pliocene relative to the pre-industrial for most oceans (Fig. 8c), although the magnitude of change might be small in areas such as the South Atlantic and eastern South Pacific.There appears to be weaker model agreement on the sign of change for the subtropical latitudes for the Pacific and Indian oceans in both hemispheres, although a slight increase in genesis potential may be expected.

Genesis frequency
Figure 9 illustrates the annual TC genesis frequency generated from the five GCMs across the various palaeoclimates as a percentage of the pre-industrial.The ensemblemean annual genesis frequency for the Pliocene, LGM and mid-Holocene are determined to be 91, 97 and 100 % of the pre-industrial respectively.
Estimating the natural variability of an ensemble mean number is problematic.As a conservative measure, we take that from the model with the highest internal variability (HadCM3) that has a SD (σ) of 2.9 %.Given that the ensemble annual frequencies are generally within the standard measure of 2σ (Haywood et al., 2013), the annual genesis frequencies for both the LGM and mid-Holocene is considered to have not deviated significantly from the pre-industrial era.Whilst the ensemble mean value for the Pliocene is statistically significant by this metric, in fact the magnitude of the reduction is driven solely by the HadCM3 member (the ensemble average without it is 100 % of 195 Figures

Back Close
Full the pre-industrial).The assumption of a Gaussian distribution inherent in this metric of significance is clearly not valid for this ensemble.It is therefore not clear we can consider the reduction seen in Pliocene ensemble as robust feature.In Fig. 10, the Northern (Southern) Hemisphere peak in JASO (JFMA) appears consistent across the various intervals, justifying this study's choice of the peak storm seasons for the respective hemisphere.Previous work from Korty et al. (2012a, b) using PMIP2 data showed a stronger peak from the Southern Hemisphere relative to the north, while this study suggests a stronger Northern Hemisphere peak.This suggests that the PMIP3 simulations may have improved accuracy in describing present day trends of a more conducive Northern Hemisphere for cyclone genesis (Gray, 1968;Klotzbach, 2006;Webster et al., 2005).
The observation of a slight increase in LGM cyclone genesis relative to the preindustrial period is consistent with earlier work by Korty et al. (2012a).The reduced TC genesis associated with the warm Pliocene conforms to the Knutson et al. (2010) view of future behaviour.It does differ from the sole prior Pliocene TC study (Fedorov et al., 2010).However that work, aimed at the early Pliocene and using a statistical downscaling approach, involved the imposition of a substantially different SST field.
For the mid-Holocene epoch, a salient increase in October activity is observed by Korty et al. (2012b), which has been attributed to a delayed SST response from the TOA insolation forcing, resulting in a shift of the Northern Hemisphere storm season.However, this trend is not observed in this study.Annual SST changes are found to have varied minimally relative to the pre-industrial (Fig. 1), suggesting that the ocean component during the mid-Holocene may play a lesser role in comparison to the Pliocene and LGM epochs where more substantial SST changes are observed.observed to be relatively insensitive to these changes (Fig. 1).Meanwhile, the annual number of tropical cyclogenesis is likewise found to be fairly consistent across the various palaeoclimates.Despite disagreement about the change of global annual TC frequency (Fig. 9), there is some model consensus on the spatial patterns of tropical cyclogenesis change.These changes may be attributable to changes in large-scale atmospheric properties such as carbon dioxide levels, altered topography and orbital forcing.
The Pliocene represents a sustained interval of elevated carbon dioxide concentration resulting in a warmer climate relative to the pre-industrial period, while the LGM era experienced an opposite cooling effect arising from lower carbon dioxide levels present at that time.Carbon dioxide, being a well-mixed greenhouse gas, causes globally coherent temperature changes in contrast to orbital forcing.In response to the warming during the Pliocene, a general decrease in genesis potential in the Northern Hemisphere and sub-tropical Southern Hemisphere is observed (Fig. 7).The Pliocene simulations have a weaker Hadley and Walker circulation that results in a broadening of the Inter-tropical Convergence Zone (ITCZ; Contoux et al., 2012).Kamae et al. (2011) show that Equatorial specific humidity during the Pliocene increased in the lower troposphere and decreased in the mid-troposphere arising from the weakened ascent of the Walker circulation.Convective processes were curtailed and an increase in moist entropy deficit might have weakened the ITCZ.The ITCZ is associated with favourable conditions for TC genesis (Gray, 1968), and the weakened ITCZ is in line with signs of general decrease in TC activity during Pliocene.
The key difference in forcing between the mid-Holocene and pre-industrial lies in the orbital parameters (Table 1).Solar insolation received in the Northern Hemisphere is enhanced relative to the Southern Hemisphere as a result of the altered precession (Braconnot et al., 2007).There is a slight tropospheric warming in the Northern Hemisphere for the middle and high latitudes as a consequence of this, while general tropospheric cooling is found in the tropical region and the Southern Hemisphere.Increased TC genesis is observed during the mid-Holocene in the Southern Hemisphere, Figures along with slight reduction in the Northern Hemisphere (Fig. 4c).This is associated with higher entropy deficit in the Northern Hemisphere hindering cyclone genesis compared to the Southern Hemisphere (not shown).The potential intensity increases slightly at all latitudes.

Relationship to future projections
Records do not currently exist to either confirm or refute the tropical cyclogenesis changes simulated by this ensemble.They probably never will.Yet the Earth will shortly experience carbon dioxide concentrations beyond those of the Pliocene period.Therefore, it is interesting to consider how the results above correspond to future projections of tropical cyclogenesis.One further motivation to do this is that the palaeoclimate simulations are all equilibrium experiments, whilst the future projections are transient.It is therefore anticipated that the climate change signal will be easier to detect in the palaeoclimate simulations.In transient simulations, large-scale forcings may not fully account for the observed variability (Menkes et al., 2012), as stochastic effects may potentially account for up to half of the observed variability (Jourdain et al., 2010).
The RCP8.5 scenario is used to project how TC genesis may develop in future.It is chosen as it is the most extreme scenario and so should have the biggest signal.In this scenario, carbon dioxide concentrations reach 936 ppmv by 2100 (Collins et al., 2013), more than double the level in the Pliocene simulations.
The GCMs selected in this study all show future changes in TC genesis (as calculated by the GPI) under the RCP8.5 transient scenario (Fig. 11).Yet these trends are not consistent between the models.Note that HadCM3 has not contributed results for RCP8.5, so a later generation of the model (HadGEM2) has been substituted.Two models suggest an increase in TC genesis, while three models suggest a decrease, resulting in an ensemble mean with a trend of slightly reduced TC genesis by 2095.
The future response is also seemingly inconsistent with the palaeoclimate responses in the same GCM.For example, MIROC shows a decrease in cyclone activity in the Figures

Back Close
Full warm Pliocene and an increase during the LGM: counter-intuitively it also shows an increase under RCP8.5.Interestingly, the multi-model mean GPI difference between the future RCP8.5 (2071-2100) scenario and historical  simulation from Camargo (2013) shows an opposite response to the difference for the equilibrium Pliocene-control difference in Fig. 8b of this study.The transient RCP8.5 GPI difference in Camargo (2013) suggests a global increase, except for a small area in the central South Pacific where a decrease is expected.Meanwhile the equilibrium Pliocene-pre-industrial difference in this study shows a general decrease, and a small region of the western South Pacific has an increase in GPI.The stark difference in cyclogenesis trends between the RCP8.5 and Pliocene therefore throws some questions on the suitability for the choice of the Pliocene as a projection of modern day greenhouse climate, at least in terms of cyclogenesis projections.Emanuel (2013) downscaled six CMIP5 GCMs for the RCP8.5 projection, and concluded that an increase in future global tropical cyclone activity might be expected.The same paper also acknowledged that other modelling groups obtained contrasting results where modest decreases (Knutson et al., 2010) and no robust change (Camargo, 2013) in future tropical cyclone activity had been detected.Two models (CCSM4 and HADGEM2-ES) that Emanuel (2013) used for his RCP8.5 scenario conclusions are also incorporated in this study, but an increasing trend is not detected for the two particular models.Potential reasons that could account for the difference include the use of a modified "clipped" vorticity GPI in this study, and a different choice of 250-850 hPa tropospheric wind shear in Emanuel (2013).The striking difference in TC activity outcomes, despite a similar GCM choice, suggests that the GPI may be highly sensitive to adjustments in diagnostic definition.Kossin et al. (2014) showed that the lifetime-maximum intensity of tropical cyclones is migrating polewards at a rate of about one degree of latitude per decade, similar to the rate of expansion of the tropics (Lucas et al., 2014).No coherent message about poleward expansion of TC genesis was found in this ensemble and changes in TC Figures

Back Close
Full genesis are found largely in the 10-20 • regions of both hemispheres, with minimum adjustment in the sub-tropics.

Conclusions
The global frequency of tropical cyclone formation is estimated to have been relatively constant over the range of past climates using a genesis potential index (GPI).
This range encompasses both greenhouse (Pliocene) and icehouse (Last Glacial Maximum) climates and changing orbital forcing.These conditions are thought to represent the extremes of climates Earth has experienced in the past three million years.At times the spread of model estimates of the global, total number of tropical cyclones can be large (Fig. 9), something that is especially noticeable for the Pliocene simulations.
The ensemble shows much higher levels of consistency on the regional scale.All five models agree on a reduction in the genesis of the North Atlantic hurricanes at the Last Glacial Maximum; compensated for by an increase in potential cyclogenesis in the central North Pacific, to a greater or lesser degree.Presumably this is a circulation response to the existence of a large ice-sheet over North America.A qualitatively similar feature has been seen previously (Korty et al., 2012a), but would obviously not be expected in the reverse in future projections.The mid-Holocene ensemble shows alterations of GPI associated with shifts in the intertropical convergence zone driven by the altered incoming solar distribution.Again the results from this ensemble are qualitatively similar to those from prior model ensembles (Korty et al., 2012b).
The tropical cyclone distribution for the mid-Pliocene has not previously been considered.Fedorov et al. (2010) anticipated an expansion of cyclone location, strength and magnitude in response to the observed weakened tropical temperature gradients of the early Pliocene (Brierley et al., 2009).The PlioMIP ensemble of the mid-Pliocene does not show those weakened gradients (Haywood et al., 2013;Brierley, 2014)  One motivation for studying the past climate tropical cyclone response is to investigate its relationship to future projections.The genesis potential under the RCP8.5 scenario was computed and contrasted with the palaeoclimate response.There is no simple relationship that emerges between global tropical cyclogenesis and say CO 2 level or hemispheric temperature.This result means that analysis of the changes in total frequency of tropical cyclone remains much less robust than the changes in the regional response noted above.The conclusion is further strengthened by the apparent sensitivity of projected future global frequency to the precise genesis potential index used -with our analysis not fully supporting either the results of Emanuel (2013) nor the opposing results of Camargo (2013) TC) constitute one of the most powerful forces of nature and can cause severe destruction to human life and property.How TC genesis may change in the face of climate change is thus an area of strong interest.Past studies using high resolution general circulation models (GCMs) have generally suggested that cyclone intensity would strengthen, yet cyclone genesis would decline in a warming climate Discussion Paper | Discussion Paper | Discussion Paper | 7 • C cooler.Given the relatively similar orbital parameters controlling Earth's solar insolation during the Pliocene, Discussion Paper | Discussion Paper | Discussion Paper | observed higher annual cyclone genesis relative to the pre-industrial era.For the mid-Holocene interval,Korty et al. (2012b) demonstrated that the difference in distribution of the topof-atmosphere (TOA) radiation in comparison to the pre-industrial control altered the seasonal cycle of potential intensity in the Northern Hemisphere.There was mixed response in TC genesis for the mid-Holocene relative to the pre-industrial period, with slightly lower TC activity in the Northern Hemisphere (and slightly higher in the Southern Hemisphere).This study aims to investigate whether similar behaviours are seen in the subsequent generation of PMIP; namely the PMIP3 model ensemble.The related Pliocene Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | s b and s * o are calculated at 925 hPa for the boundary layer and at the sea surface respectively.Physically, a larger χ m signifies a longer duration Discussion Paper | Discussion Paper | Discussion Paper | Figure 2 illustrates the genesis potential simulated from the various GCMs in their pre-industrial simulations.The Northern Hemisphere shows cyclone genesis potential for the peak storm periods of July, August, September and October (JASO), while the Southern Hemisphere corresponds to the peak storm Figures Back Close Full Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 2f) is obtained by interpolating the individual model GPI onto the HadCM3 grid which is comparatively lower in resolution.The ensemble also reveals highly similar distribution of genesis intensity for regions such as the North Atlantic, Pacific and Indian oceans in comparison with the 0 ka genesis potential from Korty et al. (2012a) simulated using PMIP2 data from seven GCMs.In both instances, Discussion Paper | Discussion Paper | Discussion Paper | Magnitude changes are more salient in the Northern Hemisphere.A slight decrease in genesis potential is seen in the North Atlantic.Changes in the North Pacific occur with dipole patterns at the periphery of the East and West Pacific indicating shifts in the location of cyclone genesisDiscussion Paper | Discussion Paper | Discussion Paper | Maximum (LGM) During the LGM, the tropics experienced cooling of −5 to −2 • C over land, while most of the tropical surface ocean did not encounter cooling beyond −2 • C especially in the Southern Hemisphere (Waelbrook et al., 2009).The LGM mean tropical SST from the five GCMs in this study during the peak storm period is 2.0 • C cooler than preindustrial.Simulated genesis potential responses for the LGM appear mixed (Fig. 5).CCSM, HadCM3 and MIROC show generally stronger genesis, while FGOALS and IPSL show a weakening in genesis potential relative to pre-industrial.All of the models show some form of compensation, indicative of shifts in the relative dominance of the TC formation locales.Discussion Paper | Discussion Paper | Discussion Paper | ), most models suggest a mixed response in the direction of change for various oceans, apart from MIROC that shows only a limited change.The majority of models indicate a decrease in genesis potential for the North Atlantic and South Indian oceans, and relatively no change in the South Atlantic Ocean.In the North Pacific Ocean, majority of models suggest a decrease in genesis potential in the eastern region, but appear to have mixed responses for the western region and the SPCZ.Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | and LGM, changes in carbon dioxide lead to sea surface temperature (SST) changes throughout the tropics, yet the potential intensity of TCs are Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | nor their associated TC genesis pattern.Rather there are large swathes of the ocean that show a robust reduction in genesis potential index.Discussion Paper | Discussion Paper | Discussion Paper | despite all three using the same simulations.Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Figure 9 .Figure 10 .Figure 11 .
Figure 9. Model and ensemble annual cyclone genesis frequency as percentage of preindustrial control frequency.

Table 1 .
Trace gases and Earth's orbital parameters recommended for PMIP.

Table 2 .
List of GCMs used in this study along with their storm frequency and variability seen in the pre-industrial simulation.