Towards Dependence of Tropical Cyclone Intensity on Sea Surface Temperature and Its Response in a Warming World

Tropical Cyclone (TC) systems affect global ocean heat transport due to mixing of the upper ocean and impact climate dynamics. A higher Sea Surface Temperature (SST), other influencing factors remaining supportive, fuels TC genesis and intensification. The atmospheric thermodynamic profile, especially the sea-air temperature contrast (SAT), also contributes due to heat transfer and affects TC’s maximum surface wind speed (Vmax) explained by enthalpy exchange processes. Studies have shown that SST can approximately be used as a proxy for SAT. As a part of an ongoing effort in this work, we simplistically explored the connection between SST and Vmax from a climatological perspective. Subsequently, estimated Vmax is applied to compute Power Dissipation Index (an upper limit on TC’s destructive potential). The model is developed using long-term observational SST reconstructions employed on three independent SST datasets and validated against an established model. This simple approach excluded physical parameters, such as mixing ratio and atmospheric profile, however, renders it generally suitable to compute potential intensity associated with TCs spatially and weakly temporally and performs well for stronger storms. A futuristic prediction by the HadCM3 climate model under doubled CO2 indicates stronger storm surface wind speeds and rising SST, especially in the Northern Hemisphere.


Introduction
In recent years there has been a growing interest in understanding the response of tropical cyclone (TC) activity to climate change [1][2][3].Most work on TC involves its track prediction that primarily needs knowledge of large-scale environmental circulation, available for most of the tropics.However, TC intensity (TCI) determination requires consideration of fine sub-grid scale processes and hence poses a challenge to the scientific community.The theoretical upper bound on a TC's intensity is termed as its Potential Intensity (PI) which can be calculated using Sea Surface Temperature (SST), Sea Level Pressure (SLP) and profiles of atmospheric temperature and mixing ratio [1,2].
There are three major techniques developed to date to explore this subject (cf., [4]).Most of the published works on the subject mention examining the intensity, based on atmospheric thermodynamics and equations governing atmospheric flows in association with TC [5][6][7][8].Another approach to determine cyclone intensity is based on a time-dependent numerical modeling technique.This method considers a weak TC seeded into a known environment.The model is then integrated with the progressing time until the cyclone's maximum intensity is accomplished [9,10].A third technique makes use of observational data to statistically find TC intensity as a function of observed oceanic-atmospheric variables, such as SST [11][12][13].For example, a previous study using ~45 year statistical regressions has shown that SST has a dominating influence in high hurricane seasons, independent of all known other factors and that sea warming explained 61% ˘34% of the anomalous Atlantic hurricane activity in 1995 to 95% [14].
All approaches have their relative merits for the intended applications and their potential areas of improvements.Generally speaking, evolved models using thermodynamic equations (physical models) are expected to yield better characterization of tropical cyclones (e.g., [5][6][7][8]), however, if these models consider cyclones to be in isolation, then the effects of the environment will be ignored.Also, they require a range of accurate input parameters that may not always be precisely available.Numerical models also require multiple parameters but most users likely use the default values (however, exact initial values are desirable) and the model must inherently be accurately parameterized.The use of observational data is based on statistical relationships (stochastic models), which certainly has the advantage of ease-of-use when data is available but may yield results which are grossly over-or underestimated and will be precise only for those points where the condition matches with the input to the regression model.It is understandably insufficient to attempt to characterize TC with a simple set of coefficients, however, they are still effective as long as the retrieval errors are in a permissible range.This third methodology is extended in this work and its results are compared with those of a physical model.The assumption, therefore, is that SST is one of the main driving forces for TC intensification.
Hazard mitigation planning and management requires a good estimate of tropical cyclone intensity and their identification, in a prompt and timely manner.As stated earlier, the genesis of thunderstorms, its intensification and conversion to cyclonic systems and dissipation is a synergistic process and depends on a number of physical parameters, such as SST, sea-air temperature contrast (SAT), SLP and mixing ratio.This paper documents the first step of an ongoing investigation where we simplistically explored the connection between SST (as the source of energy to fuel a TC system and as a proxy for SAT) and TC's maximum surface wind speed (V max ) using climatological data (not instantaneous).The aim is to see if this simplistic approach captures the major part of the variability of the maximum wind speed, as compared to a more evolved model.Provided that such simplistic approaches match developed methods at the least initially, another goal is to provide a simple empirical relation to determine cyclone intensity using statistical analysis, which incorporates sensitivity of storm strength on the ocean temperature.In future studies of this ongoing effort, we will extend our model to capture the remaining inconsistencies by including more model variables, keeping in mind the ease of application as a key requirement.Another focus of this study is to estimate changes in tropical cyclone intensity patterns in the warming world, including under a doubled CO 2 scenario, using a model simulation.
The paper is organized as follows: Section 2 briefly describes the theoretical background.The methods employed and datasets used are given in Section 3. Section 4 comprises the major results and discussions concerning: (a) development of a linear model for maximum achievable cyclone winds (V max ); (b) the relation between SST and sustained wind speed and understanding its role through enthalpy and spatial distribution of V max ; (c) application of linear model outputs to the HadCM3 Model.Temporal analyses using wavelet spectrum and time series plots are given in Section 5. Section 6 summarizes the paper and gives a future outlook.

Theoretical Background
The total power dissipation of a TC is proportional to the cube of its wind speeds [3,6,15,16].The Potential Intensity (PI) of a TC can be expressed in terms of the surface drag coefficient, the surface air density and the magnitude of the surface winds, integrated over a TC's radius ( [3,6]).Historical data does not provide information of a storm's dimension and therefore it is difficult to use historical datasets in order to determine their PI.Further studies shed light on the radial wind contours [17] and characterization of storm's dimension [18] and the drag coefficient [19].With this background and further simplifications, the Power Dissipation Index (PDI) of a TC is given by [3]: Here, V max represents the maximum sustained wind speeds at 10 meters above the surface.A TC's PDI is a better indicator of a storm's threat than a cyclone's frequency or its intensity alone [3].PDI is also of interest in reference to its influence on thermohaline circulation [20] and upper ocean mixing [3].Owing to the simplicity of practically implementing Equation (1) and its relevance to define destructiveness in association with a tropical cyclone, we have made use of this index instead of more traditional indices, such as the hurricane severity index [21], the cyclone intensity scale [22] and the accumulated cyclone energy index [23], all of which require a range of input parameters.
It is also well known that SST plays a vital role in cyclone genesis and intensification [20,24].SST in its own may not always be sufficient to fully determine TC intensification.For instance, the Atlantic storm Earl in 1992 attained a maximum speed of 28 m/s at an SST of about 27 ˝C, while the Atlantic hurricane Bonnie, during the same season, gained a maximum speed of 49 m/s with an ocean temperature of about 25 ˝C [11].However, SST is a critical parameter necessary to fuel intensification of thunderstorms that may or may not grow into a cyclone, depending on other parameters.In this study we present an empirical model to estimate cyclone intensity using SST.

Methodology
The local SST tends to influence the atmospheric state above it and affects the Potential Intensity of a cyclone [2,25,26].Temperature anomalies of the troposphere also depend upon variation in tropical SST [27].
Three datasets have been used in our study for formulation and verification purposes [28][29][30][31]; one for model formulation and the other two for independent dataset verification.The first part of this work makes use of observation-based reanalysis of datasets, ERSST.V3b [30], for constructing a simple statistical model to determine (tropical and) local maximum sustained wind speeds, V max , from (tropical and) regional SSTs.The output of Emanuel's model (with ERSST.V3b) provides TC potential intensity measured in terms of V max and P min , and provides input V max to the regression model.This simple model is then tested on two additional and independent SST input datasets, Optimally Interpolated SST (OISST.V2) [29] and COBE [31], ranging from 1800 [30] and 1891 [31] onwards, making the model suitable for computing the cyclone's strength and also for retrospective studies during the pre-satellite era.For purposes of comparison and validation, our model outputs are then compared against those of an established model [5] that is also implemented in this work (Section 3.1.1).The second part of this work focuses on exploring temporal and spatial variation of V max derived using the relationship built in the first part.After understanding how V max values change over time, in the remaining part of this work we also explore how it responds to a warmer world with doubled CO 2 .The outcome of the HadCM3 climate model is employed for this purpose.
Our simple approach, which holds only for a climatological snapshot of average values (not necessarily for instantaneous values), is based on two assumptions.First, TC's maximum surface wind speed V max values and sea-air temperature contrast (SAT) show a distinct correlation in their variabilities (cf., [32]).Second, SST and air temperature values are linearly related to a first order approximation [33,34] and thus SST and SAT (SST minus air temperature) too are correlated.Based on these assumptions, we developed a linear relationship between SST and V max whose reliability is tested further using two long-term SST reanalysis datasets.A test of hypothesis was further applied to estimate the strength of this correlation.

Computation of Sustained Wind Speeds
Sustained wind speeds, V max , are derived using monthly SST (averaged over 792 months from 1948 to 2014), sounding data, sea level pressure and specific humidity from reanalysis datasets as input to the algorithm available at ftp://texmex.mit.edu/pub/emanuel/TCMAX/pcmin_2013.f[5].In this work, we first calculate V max using the established model referenced above, using 3 different SST datasets.This algorithm also requires additional parameters (see Section 3.2.1).Then we define V maxNew as TC's sustained wind speeds computed using only SST as input to our model, as explained later in Section 4. The results of these two methods are statistically compared, showing a reasonable agreement (see Section 4.1).

Warmer World Scenario with Increased CO 2
In order to explore changes in SST and TCI when the CO 2 is doubled in the Earth's climate system, a suite of coupled oceanic-atmospheric models was employed.The model runs to get SST are obtained from the HadCM3 based on present day CO 2 levels and then in a doubled CO 2 scenario.We use Ensemble Perturbed Physics Ensembles (PPE) by estimating the mean response of 53 different model versions with changes in physical parameters involved in the model.
There are two primary reasons for exploiting the PPE technique here.The first reason is the capability of the method to sample mean climate states in relatively broad bounds [35].Secondly, a wide range of climate forcing and feedback are efficiently handled under increased levels of greenhouse gases (like CO 2 ).The demand for better computation of uncertainties in climate forecast estimates gave rise to the Perturbed Physics Technique [36].In this approach, we consider a single model structure and perturb the magnitudes of uncertain variables within an allowed range.The range is decided with the help of the people involved in parameterization development and by rigorous inspection of the modeling literature.Experimental runs with single models can then be developed in "ensemble mode", depending upon the availability of computation power and time.
Considering a changed climate situation of doubled CO 2 , SST values thus obtained from the model runs are then plugged into our statistical model to obtain V maxNew for the case.Then newly obtained V maxNew values over a grid in the tropics for the present day and the climate change scenario are subtracted to determine the net change in V max .

Reanalysis Data
A monthly mean reanalysis dataset provided by the National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) is used as input to determine tropical cyclone strength using an algorithm based on the thermodynamic model for cyclone intensity.The algorithm [15] takes SST, SLP, atmospheric temperature and mixing ratio as input variables and yields maximum achievable cyclone winds (V max ) and theoretically minimum attainable TC pressure (P min ).Unlike SST [29][30][31] and SLP [37,38] which are considered at surface level, atmospheric temperature and mixing ratio are considered in the vertical atmosphere.Atmospheric soundings are available up to the lower stratosphere for 16 pressure levels (plus the surface): 1000, 925, 850, 700, 600, 500, 400, 300, 200, 150, 100, 70, 50, 30, 20 and 10 mb.SLP, atmospheric soundings and mixing ratio datasets, used in developing the statistical model, are considered from 1948 onwards.(SST is expressed in ˝C, SLP in mb, atmospheric temperature in ˝C and the mixing ratio in g/kg.)The value of the mixing ratio ranges between 0 g/kg for cold continents to 20 g/kg for warm tropical areas.
Three climatological SSTs are used here.The first set of SST data, ERSST.V3b [30], is the National Oceanic and Atmospheric Administration (NOAA) extended reconstruction made using the most recently available Comprehensive Ocean-Atmosphere Data Sets (COADS) and improved, statistical methods that allow stable re-constructions using sparse data.(ERSST.V3b is an improved, extended reconstruction over the second version.)The second set of SST dataset, COBE SST [31], is a monthly mean of global fields created during June, 2011 at NOAA's Physical Sciences Division (PSD) of the Earth System Research Laboratory (ESRL) using the gridded data from the Japanese Reanalysis Project (JRA).The third SST dataset employed is NOAA's second version optimum interpolation OISST.V2 [29].These optimally interpolated blended data are an in situ and satellite SST combined analysis for climate [28].The satellite SSTs ingested are the Advanced Very High Resolution Radiometer (AVHRR) infrared retrievals onboard NOAA -11 through -19 polar-orbiters, with only one platform used at a time and the data are anchored to in situ drifters (approximating a bulk depth of ~20 cm).It is worth mentioning that in recent years, a substantial effort has been invested to generate blended satellite SSTs from multiple platforms, many of which also ingest in situ data in their optimal interpolation scheme [39,40].We have chosen OISST.V2 because of its extended range going back in time.The datasets, other than SST used here, were available at 2.5 ˝ˆ2.5 ˝(144 ˆ73 cells, globally), OISST.V2 and COBE SST at 1 ˝ˆ1 ˝and ERSST.V3b at 2 ˝ˆ2 ˝grids.These datasets were aggregated to a common grid of 1 ˝ˆ1 ˝employing bilinear interpolation for analyses.

Climate Model Data
In order to explore the influence of projected change of SST on cyclone intensity with enhanced levels of green house gases, we make use of climate model ensembles as derived from the third version of UK Met Office Hadley Centre Climate Model, HadCM3 [41].Considered climate model ensembles taken from the physics ensemble are retrieved from HadCM3.
HadCM3 is a coupled ocean-atmospheric general circulation model developed at the Hadley Centre, UK and is one of the primary models used for IPCC third assessment report (2001) [42].The model possesses a high resolution over the oceans and thus does not require flux correction (artificial heat and fresh water fluxes over the ocean floor) [43].The model has a grid resolution of 3.75 ˝ˆ2.5 ˝(longitude ˆlatitude; Arakawa-B type grid), which corresponds to approximately 300 km spacing between the points and is comparable to T42 truncation in a particular spectral model.The model has a time step of 30 min which includes sub-time steps, per time step of half an hour.

Development of a Linear V max Model and Its Validation against an Established Model
The assumptions behind our approach to formulate a simplistic approach to compute V max only from cyclone-prone SST values and regions are two-fold: (a) TC's V max values and sea-air temperature contrast (SAT) are correlated; (b) SST and air temperature (AT) values are approximately related, therefore are sea-air temperature contrast and SST.On the first assumption, the time-evolution of V max and SAT has been shown to be identical in their directions (not amplitude) during phases of intense storms (cf., [32]) due to enhanced vertical transport.(An offline analysis also showed significant spatial linear correlation with sea-air temperature contrast.).Secondly, mean AT calculated employing the mean value theorem, cf., equation A18 in [33], has been subsequently shown to be related to surface temperature (T s ) in a linear way.Various researchers have used this formulation, where the form of the equation remains the same but the coefficients vary depending on the area of their study.For example, AT = 55.8 + 0.77 ˆTs is reported in ([34], also see Figure 3 therein) for their study area.We have not used these values but made the assumption that SAT and SST are correlated with a scatter permissible for our purposes and developed a linear relationship between SST and V max , whose reliability is tested further using two independent inputs.
Based on the above assumptions that V max and SST for cyclone-prone conditions (all SST values in the tropics) are correlated, at least to a first order approximation, we have performed a linear regression (r 2 = 0.95) for the averaged datasets described in Section 3.1.1.Subsequently, the model is applied for climatological averaged values and with a time series analysis also performed in Section 5. Here, we denote our modeled V max as V maxNew .
V maxNew " 6.3 ˆSST ´102 (2) V maxNew, as obtained from our statistical model (Equation 2), also shows a significantly high correlation with that of V max obtained from Emanuel's model (r 2 = 0.95).In order to validate these modeled TC wind speeds, both against an established method [15] and across different SST inputs, we calculate both V max and V maxNew for each of the three SST data as inputs.Results using COBE [29][30][31] are shown in Figure 1 (COBE SST is independent of the input used for determining the model parameters) and using ERSST.V3b and OISST.V2 are supplementary, shown in Figures S1 and S2, respectively.We acknowledge that this simplistic linear model can be improved by including more parameters to reduce the scatter, however as we show here, this model works well for V max in the tropics, especially for the higher values.In addition, our goal is to be able to use and provide other users with a simple, easy-to-use model based on parameters that are readily available and in a timely manner.Figure 1 provides a visual comparison of V max from two different outputs, from Emanuel's model and from our simplistic model, using COBE SST as input.VmaxNew, as obtained from our statistical model (Equation 2), also shows a significantly high correlation with that of Vmax obtained from Emanuel's model (r 2 = 0.95).In order to validate these modeled TC wind speeds, both against an established method [15] and across different SST inputs, we calculate both Vmax and VmaxNew for each of the three SST data as inputs.Results using COBE [29][30][31] are shown in Figure 1 (COBE SST is independent of the input used for determining the model parameters) and using ERSST.V3b and OISST.V2 are supplementary, shown in Figure S1 and S2, respectively.We acknowledge that this simplistic linear model can be improved by including more parameters to reduce the scatter, however as we show here, this model works well for Vmax in the tropics, especially for the higher values.In addition, our goal is to be able to use and provide other users with a simple, easy-to-use model based on parameters that are readily available and in a timely manner.Figure 1 provides a visual comparison of Vmax from two different outputs, from Emanuel's model and from our simplistic model, using COBE SST as input.The structural resemblance between Figure 1 (row 1) and Figure 1 (row 2) implies that VmaxNew is capable of capturing most of the Vmax structures in the tropics, which is encouraging given our simplistic approach.However on an average, our model shows somewhat higher Vmax values of 2.89 V maxNew employing our approach (Equation 2) with COBE SST; Row-3: Spatial distribution of V maxNew minus V max to visually identify areas of disagreement, considering V max as a benchmark (reference); Row-4: Probability density function (PDF) of "∆V max = V maxNew ´Vmax ".Statistical parameters (number of V max matches, minimum, maximum, mean, standard deviation, root mean square difference), a Gaussian fit, X~N (mean, Std Dev) and number of grids beyond "Mean ˘4ˆStd Deviation", henceforth called 'extreme grids' are also annotated on the PDF.
The structural resemblance between Figure 1 (row 1) and Figure 1 (row 2) implies that V maxNew is capable of capturing most of the V max structures in the tropics, which is encouraging given our simplistic approach.However on an average, our model shows somewhat higher V max values of 2.89 m/s mean difference as seen in Figure 1 (row 4).These 'higher differences' are relatively evenly distributed to a smaller magnitude throughout the tropics and distinctly off the coasts of Angola, Namibia, Senegal and Peru to higher magnitudes.These higher values are seen on the histogram in Figure 1 (row 4) as ~0.45% points beyond "Mean ˘4 ˆStd Dev".
Supplementary figures show analyses as in Figure 1 but using ERSST.V3b (Figure S1) and OISST.V2 (Figure S2) as SST inputs.For Emanuel's model, additional data for other required parameters have been used.The results of Figures S1 and S2 are consistent with the results observed from Figure 1, with a slight change of statistical parameters and hence are excluded from re-discussion for brevity.
High correlations between SST and V max simulations for the three datasets (as mentioned later in this section) correspond to a linear relation between the two variables, thus simplifying the calculation of PDI, at least to a first order approximation.As shown in Table 1, this simplistic approach yields results approximately similar to those from the more elaborate approach that, however, requires a range of inputs.It is also important to mention here that this approach is valid only for time-integrated, low-spatial resolution datasets (e.g., monthly mean averages at 1 ˝ˆ1 ˝latitude-longitude grids) and for SST conditions potentially suitable for TC genesis (SSTs>25 ˝C).By no means, however, does it imply that this relationship between SST and V max will hold globally for all SST values and in short time-scales (for example on a daily scale and in daytime low wind-speed conditions with reduced mechanical mixing, warming will be higher.)Put differently, on short time scales and high spatial resolutions, SST only cannot predict instantaneous 10 m maximum wind-speeds.
The minor dissimilarities across outputs of various SST inputs, e.g., OISST.V2 and COBE can be further attributed to minor variations in the distribution of SST across the datasets.Discrepancy between the considered SST products can be attributed to difference in the statistical reconstructions used to find SST values in the regions where data was not available.A second reason could be in association with the techniques applied in collecting various SST datasets, for instance, going from the "Bucket" to the "Hull" method.In the Bucket approach, SST is considered to be the water temperature of buckets of water taken from oceans.In the "hulls of ships" method, SST is measured using sensors placed at the hull of a ship.The data products use both in situ and satellite data to obtain SST covering all the ocean basins on earth.(Note that the cause of differences in SST is not the focus of this study; readers are referred to [28][29][30][31]39,40] for more information on blended analysis SST fields.) We have also studied the effect of differences between various SST analysis fields on the model output (figures not shown for brevity).For example, SAT and SST as obtained from [37] and COBE SST respectively, show significant correlation (r 2 = 0.98) and so do SAT and V max (r 2 = 0.93).These strong connections are then verified using the bootstrapping method [44].Modeled TC wind speeds, V maxNew , is appreciably connected with a linear correlation of r 2 = 0.90 to the storm winds V max , as computed from reanalysis datasets using ERSST.V3b (figures not shown).Another set of SST values (OISST.V2) employed to compute V maxNew also give high correlation with V max (r 2 = 0.89).These results thus validate our statistical model and further support it.The statistical significance of these correlations was further confirmed by using null hypothesis where p is a matrix of p-values to test the hypothesis of no correlation.In our case, the p-value was found to be close to zero, implying that the correlation is significant.
Since reliable atmospheric temperature profiles before the satellite-era at remote regions are unavailable, we make use of a relationship between PDI and SST which can be used to determine cyclone power in remote regions.It is known [3] that maximum TC wind speed is better correlated with ocean-atmosphere temperature difference than with SST.That is, though warmer SSTs contribute towards a highly supportive environment for TC strengthening, its formation and intensification are strongly influenced by the contrast between ocean temperature and the vertical thermodynamic profile of the atmosphere above it.However, this result suggests that when vertical moisture content and SLP are not available, SST values alone can be used by simply employing the linear regression model suggested here, to approximate TC wind speeds, V maxNew , and hence the intensity of a cyclone as reflected by the PDI.This result showing indirect association between SST and sustained wind speed can be explained by considering the contribution of enthalpy exchange between the ocean and the air above during cyclone genesis.A tropical cyclone draws heat energy from the warm ocean water and uses it to evolve and move across the underlying surface against friction.Thus, a TC system can be regarded as a heat engine.At equilibrium, the input energy is equal to the energy going out of the system or is used in a mechanical motion against the surface friction.Using simplifications as described by previous work (cf., [5]) and considering the bulk aerodynamic formula for evaporation [45] and further simplifications as shown earlier in [5], the rate of total energy generated per unit surface area is given by: where, T s and T 0 refer to temperature at the surface of the ocean (SST) and the atmosphere above it (AT) respectively; C k represents the enthalpy exchange coefficient; ρ denotes the density of near surface air in kg/m 3 ; v represents wind speeds close to the ground in m/s; ∆k is the difference of enthalpy from the sea surface and the air above it and C d is the surface drag coefficient.Equation (3) provides a sound thermodynamic approximation of the behavior of SST and the atmospheric temperature profile during TC genesis in relation to changes in enthalpy (∆k), though deriving this equation involves various assumptions and hence has its own limitations.

Spatial Distribution of SST and Comparison with that of V max
It is observed in Figure 2 that higher temperatures responsible for tropical cyclone intensification are prevalently present along the equator and primarily accumulated in the Northern Hemisphere.In general, the Northern Hemisphere warms more than the Southern Hemisphere, mostly due to more landmass and less ocean surface in the north and the fact that sea-water heats up relatively slowly.Owing to a wider area occupied by the Pacific Ocean and warmer SSTs in the Pacific in comparison to other oceans, a higher amount of heat energy is available to an evolving TC.However, vertical atmospheric sounding also contributes to cyclogenesis (Section 4.1).SST and theoretically maximum sustained winds speeds, Vmax, show roughly similar structures (cf., Figure 1 and 2).

Application of Linear Model Outputs to HadCM3 Model
We examined changes in SST and hence VmaxNew and PDI projected from the third version of the Hadley Centre Climate Model, HadCM3, where coupled oceanic-atmospheric model outputs were considered in the present day and for a doubled CO2 scenario.We further applied the resulting ensemble outputs generated by perturbed physics outcomes from the HadCM3 climate model to our statistical model of determining TC winds.This was done to investigate the change in maximum sustained winds and hence cyclone strength ( max 3 ) when CO2 levels are doubled in our atmosphere.
It should be noted that the climate model outputs used here are inefficient at resolving the details of tropical cyclone intensity.Both reanalysis and model projections of cyclone wind speeds are found to be highly correlated to SAT, which is in turn are strongly linked to SST.Modeled VmaxNew derived from SST again show high correlation between the reanalysis and HadCM3 model outputs for present day CO2 level.Using estimated values of SST from the climate model outputs employing both present day and doubled CO2 cases in our statistical model (Equation ( 2)), we computed VmaxNew for the two cases mentioned before.The difference in the wind speeds obtained with enhanced CO2 levels show a significant increase in the maximum sustained winds as shown in Figure 3 and hence cyclone strength.This can be interpreted as stronger storms in the warming world with higher levels of CO2.Increase in maximum sustained wind speeds is heterogeneous in both hemispheres, with higher increases in the Northern Hemisphere.
The Northeast Pacific shows the highest rise in storm winds [35].The North Atlantic Ocean comes next in the list of regions.The Indian Ocean followed by the West Pacific region show little influence due to increased levels of CO2.This change is reversed in the Southern Hemisphere with the West Pacific and Indian Oceans representing regions most prone to destruction by TCs.However, possibly significant statistical differences of the maximum wind speed under a double CO2 climate appear in the Southeast Pacific.To summarize, it is some sort of dipole effect with regions in the North showing the highest increase, reflecting the lowest increase in the South.Owing to a wider area occupied by the Pacific Ocean and warmer SSTs in the Pacific in comparison to other oceans, a higher amount of heat energy is available to an evolving TC.However, vertical atmospheric sounding also contributes to cyclogenesis (Section 4.1).SST and theoretically maximum sustained winds speeds, V max , show roughly similar structures (cf., Figures 1 and 2).

Application of Linear Model Outputs to HadCM3 Model
We examined changes in SST and hence V maxNew and PDI projected from the third version of the Hadley Centre Climate Model, HadCM3, where coupled oceanic-atmospheric model outputs were considered in the present day and for a doubled CO 2 scenario.We further applied the resulting ensemble outputs generated by perturbed physics outcomes from the HadCM3 climate model to our statistical model of determining TC winds.This was done to investigate the change in maximum sustained winds and hence cyclone strength (V 3 max ) when CO 2 levels are doubled in our atmosphere.It should be noted that the climate model outputs used here are inefficient at resolving the details of tropical cyclone intensity.Both reanalysis and model projections of cyclone wind speeds are found to be highly correlated to SAT, which is in turn are strongly linked to SST.Modeled V maxNew derived from SST again show high correlation between the reanalysis and HadCM3 model outputs for present day CO 2 level.Using estimated values of SST from the climate model outputs employing both present day and doubled CO 2 cases in our statistical model (Equation ( 2)), we computed V maxNew for the two cases mentioned before.The difference in the wind speeds obtained with enhanced CO 2 levels show a significant increase in the maximum sustained winds as shown in Figure 3 and hence cyclone strength.This can be interpreted as stronger storms in the warming world with higher levels of CO 2 .Increase in maximum sustained wind speeds is heterogeneous in both hemispheres, with higher increases in the Northern Hemisphere.
The Northeast Pacific shows the highest rise in storm winds [35].The North Atlantic Ocean comes next in the list of regions.The Indian Ocean followed by the West Pacific region show little influence due to increased levels of CO 2 .This change is reversed in the Southern Hemisphere with the West Pacific and Indian Oceans representing regions most prone to destruction by TCs.However, possibly significant statistical differences of the maximum wind speed under a double CO 2 climate appear in the Southeast Pacific.To summarize, it is some sort of dipole effect with regions in the North showing the highest increase, reflecting the lowest increase in the South.Using Equation ( 2), we temporally determined a proxy index for Vmax, the VmaxNew and then computed the same parameter (VmaxNew) from temporal reanalysis outputs.The linear correlation between the two time series comes out to be only 0.86 in the tropical Atlantic Ocean and about 0.74 in tropical parts of the West Pacific Ocean.Thus, we cannot efficiently employ the spatially derived Vmax model to determine the parameter varying temporally.This could be due to the non-stationary behavior of maximum sustained winds.However, as we showed earlier the spatially derived proxy index (VmaxNew), averaged over time, successfully determines maximum sustained wind speeds in the tropics and the cyclone prone ocean basins.
It should be noted that TCs rarely occur within 5° of the Equator due to the negligible value of planetary angular momentum [46].However, TC wind speed Vmax, derived using [1,15] as shown in Figure 1 does show high values due to the absence of this concept in the model and has been excluded from interpretation, although shown on the plot for contouring purposes.

Model Application on Temporal Resolution
In order to analyze if this simplistic approach holds adequately in the temporal domain, a time series analysis is performed using wavelet spectrum and standard time series plots.The conclusion that is drawn on spatial domain that this approach holds is also verified for longer time series, implying that the simplistic approach may be used as a fallback for TC studies when reliable atmospheric temperature profiles before the satellite-era, at remote regions, are unavailable for evolved models.
As a first step, we have used wavelet analyses on our output using ERSST.V3.b, individually for all the basins that provided us with the frequency of the VmaxNew signals and their associated timeperiods.The major known trends were found to be identified from this simplistic approach, which is encouraging.Whereas the outputs of wavelet spectrum analysis show the frequency vs. time pattern of VmaxNew, it does not provide a clear picture of the temporal behavior.Therefore, as a second step, time series variation of VmaxNew is shown for the corresponding basins.

Patterns in the VmaxNew
Continuous wavelet power spectra for VmaxNew in different basins are plotted here for further analyses.The black contours on Figures 4 through to 8 represent the 5% significance level against the red noise and the cone of influence demarcates the region where edge effects come into effect.Using Equation ( 2), we temporally determined a proxy index for V max , the V maxNew and then computed the same parameter (V maxNew ) from temporal reanalysis outputs.The linear correlation between the two time series comes out to be only 0.86 in the tropical Atlantic Ocean and about 0.74 in tropical parts of the West Pacific Ocean.Thus, we cannot efficiently employ the spatially derived V max model to determine the parameter varying temporally.This could be due to the non-stationary behavior of maximum sustained winds.However, as we showed earlier the spatially derived proxy index (V maxNew ), averaged over time, successfully determines maximum sustained wind speeds in the tropics and the cyclone prone ocean basins.

North Indian Ocean
It should be noted that TCs rarely occur within 5 ˝of the Equator due to the negligible value of planetary angular momentum [46].However, TC wind speed V max , derived using [1,15] as shown in Figure 1 does show high values due to the absence of this concept in the model and has been excluded from interpretation, although shown on the plot for contouring purposes.

Model Application on Temporal Resolution
In order to analyze if this simplistic approach holds adequately in the temporal domain, a time series analysis is performed using wavelet spectrum and standard time series plots.The conclusion that is drawn on spatial domain that this approach holds is also verified for longer time series, implying that the simplistic approach may be used as a fallback for TC studies when reliable atmospheric temperature profiles before the satellite-era, at remote regions, are unavailable for evolved models.
As a first step, we have used wavelet analyses on our output using ERSST.V3.b, individually for all the basins that provided us with the frequency of the V maxNew signals and their associated time-periods.The major known trends were found to be identified from this simplistic approach, which is encouraging.Whereas the outputs of wavelet spectrum analysis show the frequency vs. time pattern of V maxNew , it does not provide a clear picture of the temporal behavior.Therefore, as a second step, time series variation of V maxNew is shown for the corresponding basins.

Patterns in the V maxNew
Continuous wavelet power spectra for V maxNew in different basins are plotted here for further analyses.The black contours on Figures 4-8 represent the 5% significance level against the red noise and the cone of influence demarcates the region where edge effects come into effect.

North Atlantic Ocean
The continuous wavelet power spectrum in the North Atlantic is illustrated in Figure 6.A significant peak in the 5 years band around 1934 is apparent.The wavelet spectrum also shows high power with the periodicity of 2-3 years which is widespread and show significant but intermittent occurrences at a difference of about 8-10 years.

North Atlantic Ocean
The continuous wavelet power spectrum in the North Atlantic is illustrated in Figure 6.A significant peak in the 5 years band around 1934 is apparent.The wavelet spectrum also shows high power with the periodicity of 2-3 years which is widespread and show significant but intermittent occurrences at a difference of about 8-10 years.

North Atlantic Ocean
The continuous wavelet power spectrum in the North Atlantic is illustrated in Figure 6.A significant peak in the 5 years band around 1934 is apparent.The wavelet spectrum also shows high power with the periodicity of 2-3 years which is widespread and show significant but intermittent occurrences at a difference of about 8-10 years.It should be noted that the findings from this part of the study shows various significant cycles in the considered ocean basins that may be associated with ENSO, PDO, North Atlantic Oscillation and IOD, as inferred in the description above.However, more detailed analysis is recommended in each of the ocean basins in relation with such natural cycles.(2) The South Indian Ocean: Trends in VmaxNew in the South Indian Ocean (Figure 10) show a lull period between 1854 and 1940 followed by a year of sharp rise and then fall.VmaxNew since 1970 onwards illustrates a sharp rise from 59 m/s to 63 m/s.It should be noted that the findings from this part of the study shows various significant cycles in the considered ocean basins that may be associated with ENSO, PDO, North Atlantic Oscillation and IOD, as inferred in the description above.However, more detailed analysis is recommended in each of the ocean basins in relation with such natural cycles.It should be noted that the findings from this part of the study shows various significant cycles in the considered ocean basins that may be associated with ENSO, PDO, North Atlantic Oscillation and IOD, as inferred in the description above.However, more detailed analysis is recommended in each of the ocean basins in relation with such natural cycles.(2) The South Indian Ocean: Trends in VmaxNew in the South Indian Ocean (Figure 10) show a lull period between 1854 and 1940 followed by a year of sharp rise and then fall.VmaxNew since 1970 onwards illustrates a sharp rise from 59 m/s to 63 m/s.(2) The South Indian Ocean: Trends in V maxNew in the South Indian Ocean (Figure 10) show a lull period between 1854 and 1940 followed by a year of sharp rise and then fall.V maxNew since 1970 onwards illustrates a sharp rise from 59 m/s to 63 m/s.(4) The North Indian Ocean: The North Indian Ocean looks like an active region accommodating various intense tropical storms since early 1990 (Figure 12).Although there was a comparatively dormant period during 1895 to 1990, a dramatic rise is observed thereafter.(5) The North Atlantic Ocean: The North Atlantic Ocean seems to drift from moderately active to a dramatically active TC region, in terms of the associated destructive potential (Figure 13).Year 2005 and 2010 marked years of severe tropical storms.The period prior to 1910 was of low cyclone activity.

Time Series Trends
Of all these discussed ocean basins, the Northwest Pacific and the North Indian Ocean represent the regions of high intensity storms, followed by the North Atlantic.The South Indian Ocean accommodates the least active storms of all the considered oceans.However on an average, the TC intensity in the South Indian Ocean shows gradual growth since 1960.The Southwest Pacific also follows a rising trend in the TC intensity since 1960 onwards.(4) The North Indian Ocean: The North Indian Ocean looks like an active region accommodating various intense tropical storms since early 1990 (Figure 12).Although there was a comparatively dormant period during 1895 to 1990, a dramatic rise is observed thereafter.(5) The North Atlantic Ocean: The North Atlantic Ocean seems to drift from moderately active to a dramatically active TC region, in terms of the associated destructive potential (Figure 13).Year 2005 and 2010 marked years of severe tropical storms.The period prior to 1910 was of low cyclone activity.
Of all these discussed ocean basins, the Northwest Pacific and the North Indian Ocean represent the regions of high intensity storms, followed by the North Atlantic.The South Indian Ocean accommodates the least active storms of all the considered oceans.However on an average, the TC intensity in the South Indian Ocean shows gradual growth since 1960.The Southwest Pacific also follows a rising trend in the TC intensity since 1960 onwards.(4) The North Indian Ocean: The North Indian Ocean looks like an active region accommodating various intense tropical storms since early 1990 (Figure 12).Although there was a comparatively dormant period during 1895 to 1990, a dramatic rise is observed thereafter.(5) The North Atlantic Ocean: The North Atlantic Ocean seems to drift from moderately active to a dramatically active TC region, in terms of the associated destructive potential (Figure 13).Year 2005 and 2010 marked years of severe tropical storms.The period prior to 1910 was of low cyclone activity.
Of all these discussed ocean basins, the Northwest Pacific and the North Indian Ocean represent the regions of high intensity storms, followed by the North Atlantic.The South Indian Ocean accommodates the least active storms of all the considered oceans.However on an average, the TC intensity in the South Indian Ocean shows gradual growth since 1960.The Southwest Pacific also follows a rising trend in the TC intensity since 1960 onwards.The Northern Hemispheric oceans can thus be regarded as the more active region for TC intensification, as compared to the Southern Hemisphere.
Decadal Trends: Decadal fluctuations are most prominent in the Southwest Pacific and the North Atlantic.Other oceans show either small cycles or a very long cycle as in the Northwest Pacific.All the TC prone basins however show a dramatic rise since 1970.It is worth mentioning that here we only briefly describe the observed overall tendency.However, trend detection is a separate subject of investigation where it is necessary to isolate the different components from the time series (seasonal, error in data and the trend) and quantify the trend for a chosen time period e.g., decadal or centennial.
To conclude this section, similar to the steep rise predicted using a full calculation model [47,48], statistically simulated wind speeds show a sharp rise post the satellite era, after 1970s.

Summary and Future Work
The three considered SST datasets support the notion that in tropical cyclone-prone areas, oceanic-atmospheric temperature contrast and, thus, SST strongly contributes to maximum surface wind speeds of TCs.This high correlation is explained by the enthalpy exchange process between the  The Northern Hemispheric oceans can thus be regarded as the more active region for TC intensification, as compared to the Southern Hemisphere.
Decadal Trends: Decadal fluctuations are most prominent in the Southwest Pacific and the North Atlantic.Other oceans show either small cycles or a very long cycle as in the Northwest Pacific.All the TC prone basins however show a dramatic rise since 1970.It is worth mentioning that here we only briefly describe the observed overall tendency.However, trend detection is a separate subject of investigation where it is necessary to isolate the different components from the time series (seasonal, error in data and the trend) and quantify the trend for a chosen time period e.g., decadal or centennial.
To conclude this section, similar to the steep rise predicted using a full calculation model [47,48], statistically simulated wind speeds show a sharp rise post the satellite era, after 1970s.

Summary and Future Work
The three considered SST datasets support the notion that in tropical cyclone-prone areas, oceanic-atmospheric temperature contrast and, thus, SST strongly contributes to maximum surface wind speeds of TCs.This high correlation is explained by the enthalpy exchange process between the The Northern Hemispheric oceans can thus be regarded as the more active region for TC intensification, as compared to the Southern Hemisphere.
Decadal Trends: Decadal fluctuations are most prominent in the Southwest Pacific and the North Atlantic.Other oceans show either small cycles or a very long cycle as in the Northwest Pacific.All the TC prone basins however show a dramatic rise since 1970.It is worth mentioning that here we only briefly describe the observed overall tendency.However, trend detection is a separate subject of investigation where it is necessary to isolate the different components from the time series (seasonal, error in data and the trend) and quantify the trend for a chosen time period e.g., decadal or centennial.
To conclude this section, similar to the steep rise predicted using a full calculation model [47,48], statistically simulated wind speeds show a sharp rise post the satellite era, after 1970s.

Summary and Future Work
The three considered SST datasets support the notion that in tropical cyclone-prone areas, oceanic-atmospheric temperature contrast and, thus, SST strongly contributes to maximum surface wind speeds of TCs.This high correlation is explained by the enthalpy exchange process between the sea surface and the atmospheric profile above it, considering the TC system as a natural Carnot engine.The higher this contrast is, the greater is the energy transferred upwards during cyclone genesis.According to the simple linear model developed based on significant correlations, it is possible to determine TC intensity primarily using only SST values at the corresponding time scales, with somewhat suboptimal but acceptable performance, when other data (SLP, atmospheric profile) are missing.In this paper, we explored this possibility using a simplistic approach, where the results indicate that they match another evolved models to a first degree of approximation for climatologically averaged datasets.This assumes importance in light of the fact that it is challenging to obtain parameters like mixing ratio and other atmospheric temperature profile variables at remote, less accessible locations or for retrospective studies during the pre-satellite era.In our future work, we will first extend the method to reduce the scatter using other parameters and then test its performance on different categories of tropical cyclones, using available data for known major cyclones.
Tropical cyclone intensities vary in the different ocean basins.Changes in the simulated surface wind speeds under 2ˆand 1ˆCO 2 levels show higher values for the Northeast Pacific and the North Atlantic Oceans and comparatively lower values for the North Indian and the Northwest Pacific Oceans (Figure 3).This situation reverses in the Southern Hemisphere with higher wind speeds in the Southwest Pacific and the South Indian Oceans than that in the South Atlantic and the East Pacific.Also, statistically significant differences of the maximum wind speed under a double CO 2 climate appear in the Southeast Pacific.Spatial distribution of simulated wind speeds in the Pacific Ocean shows patterns similar to the El Niño Southern Oscillation (ENSO) warm pool and cold tongue.This pattern is unambiguously visible in reanalysis and model derived wind speeds as well.Variability of these wind speed structures indicates a possible connection of ENSO and maximum sustained winds.This aspect needs more rigorous investigation and is beyond the scope of this work.
Interestingly, correlation between oceanic atmospheric temperature differences and simulated surface wind speeds is stronger for major cyclones (belonging to category 3 or more).That is, stronger storms follow the linear relationship obtained through this model, better than weaker storms.One possible explanation of this could be the role of WISHE (Wind-Induced Surface Heat Exchange) in deciding the weaker or stronger storms (in terms of V max ) and hence their relation with the oceanic atmospheric temperature difference.
Model analysis shows that high surface wind speeds associated with stronger TCs are located in relatively lower, up to approximately 20 ˝latitudes, in both the hemispheres.Also, the Northern Hemisphere warms more and usually experiences higher cyclone winds than the Southern Hemisphere.
The North Atlantic and the East Pacific Ocean basins are the regions more prone to high category storms than the West Pacific and the Indian Ocean.The peculiarity of a vast land region, north of the Indian Ocean, could account for the relatively diminished heat energy available to feed a TC, in comparison to other oceans.
In a futuristic prediction by a climate model, both tropical SST and a storm's surface wind speeds increase.The computed difference between the two V max values (with and without when using a doubled CO 2 scenario) in the Indian Ocean, the West Pacific and the North Atlantic Ocean indicate a rise in PDI when doubling CO 2 in our environment.In addition, significant differences of maximum wind speeds under a double CO 2 climate also appear in the Southeast Pacific.
The following is a summary of the key observations: (a) Correlation between the ocean temperature and V max is stronger for higher values of the latter (figure not shown); (b) Northern Hemisphere shows comparatively stronger storms than the Southern Hemisphere, owing to a higher ocean warming in the Northern Hemisphere; (c) The North Atlantic and the Pacific Oceans are prone to stronger storms; (d) In a futuristic scenario of doubled CO 2 , both SST and V max increases.The oceans more prone to such intensification are the North Indian, the West Pacific and the Southeast Pacific basins; (e) Tropical cyclone intensity is highly dependent on ocean temperature.SST can be considered as one of the main driving forces for tropical cyclone intensification.Hence, in the absence of other parameters, we can use SST to simulate tropical cyclone intensity to a reasonable extent for global tropical oceans as has also been demonstrated by previous case studies (cf., [49]).
The following are the main limitations of this work: (a) This work employed the WISHE model approach which is also debated in some research results [49].Thus we suggest a further investigation using a model closer to the non-axisymmetric spin-up state of an ideal tropical cyclone; (b) The developed model is based on SST alone which has an advantage while studying tropical cyclone intensity at locations where all the needed atmospheric profiles are missing and SST is available.However, any error in SST will be propagated to the estimated wind speeds as there is no mechanism to prevent error propagation in such an overly simplistic approach; (c) The current validation is performed as follows: Verification using two independent SSTs, other than the one used to derive the equation and comparison against Emanuel's model output.However, our results are not compared against actual observations, which remain the subject for future studies.Nevertheless, it should be noted that no matter what we validate against (product minus reference), the statistical numbers (e.g., mean, standard deviation) have contributions from both the "target product" and the "reference product" and are therefore subject to some error.To alleviate this issue and characterize true product-specific random errors, a triple-collocation method (cf., [50] employed on near-surface wind fields) needs to be explored in an extension of this work.

Figure 1 .
Figure 1.Comparison of cyclone wind speed Vmax obtained from Emanuel's model and our simplistic empirical approach with only SST as input.Row-1: Vmax employing Emanuel's model with reanalysis datasets as inputs: SLP, specific humidity, atmospheric temperature profile and COBE SST; Row-2: VmaxNew employing our approach (Equation 2) with COBE SST; Row-3: Spatial distribution of VmaxNew minus Vmax to visually identify areas of disagreement, considering Vmax as a benchmark (reference); Row-4: Probability density function (PDF) of "∆Vmax = VmaxNew − Vmax".Statistical parameters (number of Vmax matches, minimum, maximum, mean, standard deviation, root mean square difference), a Gaussian fit, X~N (mean, Std Dev) and number of grids beyond "Mean ± 4×Std Deviation", henceforth called 'extreme grids' are also annotated on the PDF.

Figure 1 .
Figure 1.Comparison of cyclone wind speed V max obtained from Emanuel's model and our simplistic empirical approach with only SST as input.Row-1: V max employing Emanuel's model with reanalysis datasets as inputs: SLP, specific humidity, atmospheric temperature profile and COBE SST; Row-2: V maxNew employing our approach (Equation2) with COBE SST; Row-3: Spatial distribution of V maxNew minus V max to visually identify areas of disagreement, considering V max as a benchmark (reference); Row-4: Probability density function (PDF) of "∆V max = V maxNew ´Vmax ".Statistical parameters (number of V max matches, minimum, maximum, mean, standard deviation, root mean square difference), a Gaussian fit, X~N (mean, Std Dev) and number of grids beyond "Mean ˘4ˆStd Deviation", henceforth called 'extreme grids' are also annotated on the PDF.
less ocean surface in the north and the fact that sea-water heats up relatively slowly.

Figure 4 19 Figure 4
Figure 4 plots the power spectrum in the North Indian Ocean that shows high power in ~6 year band during 1877 and a 5-year cycle in 1901.The two elongated bands (contoured in black) range between ~2-6 years.After these two bands, we can observe a lull period and then a cycle between 3-4 years during 1943, followed by a similar significant cycle around 1984 and again during 1988.All these bands show highest power around 5 years and 3 years that may be associated with the El Niño-Southern Oscillation (ENSO) in the Northwest Pacific and dormant Indian Ocean Dipole cycle respectively.Climate 2016, 4, 30 11 of 19 Figure 4 plots the power spectrum in the North Indian Ocean that shows high power in ~6 year band during 1877 and a 5-year cycle in 1901.The two elongated bands (contoured in black) range between ~2-6 years.After these two bands, we can observe a lull period and then a cycle between 3-4 years during 1943, followed by a similar significant cycle around 1984 and again during 1988.All these bands show highest power around 5 years and 3 years that may be associated with the El Niño-Southern Oscillation (ENSO) in the Northwest Pacific and dormant Indian Ocean Dipole cycle respectively.

Figure 4 .
Figure 4. Wavelet power spectrum of VmaxNew in the North Indian Ocean highlighting temporal patterns of occurrence of maximum wind speeds.5.1.2.Southwest Pacific Ocean The continuous wavelet transform of VmaxNew in the Southwest Pacific is shown in Figure 5. Cyclic patterns in the wavelet power are clearly observed in the VmaxNew series.A significant peak in the 2-9 years cycle is apparent around 1954.This cycle of ~9 years may coincide with the Pacific Decadal Oscillations (PDO) which spans 8-12 years period.About a 4 year cycle is also seen around 1972.Another significant cycle possesses a period of 3-7 years during 1910.These oscillations may coincide with the ENSO cycle of 3-7 years (needs further investigation).

Figure 5 .
Figure 5. Wavelet power spectrum of VmaxNew in the Southwest Pacific Ocean highlighting temporal patterns of occurrence of maximum wind speeds.

Figure 4 . 19 Figure 4
Figure 4. Wavelet power spectrum of V maxNew in the North Indian Ocean highlighting temporal patterns of occurrence of maximum wind speeds.5.1.2.Southwest Pacific OceanThe continuous wavelet transform of V maxNew in the Southwest Pacific is shown in Figure5.Cyclic patterns in the wavelet power are clearly observed in the V maxNew series.A significant peak in the 2-9 years cycle is apparent around 1954.This cycle of ~9 years may coincide with the Pacific Decadal Oscillations (PDO) which spans 8-12 years period.About a 4 year cycle is also seen around 1972.Another significant cycle possesses a period of 3-7 years during 1910.These oscillations may coincide with the ENSO cycle of 3-7 years (needs further investigation).

Figure 4 .
Figure 4. Wavelet power spectrum of VmaxNew in the North Indian Ocean highlighting temporal patterns of occurrence of maximum wind speeds.5.1.2.Southwest Pacific Ocean The continuous wavelet transform of VmaxNew in the Southwest Pacific is shown in Figure 5. Cyclic patterns in the wavelet power are clearly observed in the VmaxNew series.A significant peak in the 2-9 years cycle is apparent around 1954.This cycle of ~9 years may coincide with the Pacific Decadal Oscillations (PDO) which spans 8-12 years period.About a 4 year cycle is also seen around 1972.Another significant cycle possesses a period of 3-7 years during 1910.These oscillations may coincide with the ENSO cycle of 3-7 years (needs further investigation).

Figure 5 .
Figure 5. Wavelet power spectrum of VmaxNew in the Southwest Pacific Ocean highlighting temporal patterns of occurrence of maximum wind speeds.

Figure 5 .
Figure 5. Wavelet power spectrum of V maxNew in the Southwest Pacific Ocean highlighting temporal patterns of occurrence of maximum wind speeds.

Figure 6 .
Figure 6.Wavelet power spectrum of VmaxNew in the North Atlantic Ocean.

Figure 7 .
Figure 7. Wavelet power spectrum of VmaxNew in the Northwest Pacific.

Figure 6 .
Figure 6.Wavelet power spectrum of V maxNew in the North Atlantic Ocean.5.1.4.Northwest Pacific Figure 7 presents the continuous wavelet power spectrum of V maxNew in the Northwest Pacific which shows a significant cycle of ~8 years around 1950 and may be related to the Pacific Decadal Oscillations (PDO) that show up around 1916.Cycles of 2-3 years are widespread during the considered period of 161 years from 1854 to 2014.

Figure 6 .
Figure 6.Wavelet power spectrum of VmaxNew in the North Atlantic Ocean.

Figure 7 .
Figure 7. Wavelet power spectrum of VmaxNew in the Northwest Pacific.

Figure 7 .
Figure 7. Wavelet power spectrum of V maxNew in the Northwest Pacific.5.1.5.Southern Indian Ocean Figure 8 exhibits the continuous wavelet power spectrum of V maxNew in the South Indian Ocean.The figure suggests the occurrence of very few cycles.Most of the power of the signal appears to be

Figure 8 .
Figure 8. Wavelet power spectrum of VmaxNew in the Southern Indian Ocean.

( 1 )
The Southwest Pacific: Simulated TC wind speed VmaxNew reveals a steep increase since 1910 until 2014 in the Southwest Pacific (Figure 9).A period of 161 years between 1854 to 2014 shows steady fluctuations between 67 to 70 m/s wind speed.

Figure 9 .
Figure 9.Time series of Simulated TC wind speed VmaxNew in the Southwest Pacific.

Figure 8 .
Figure 8. Wavelet power spectrum of V maxNew in the Southern Indian Ocean.

( 1 )
The Southwest Pacific: Simulated TC wind speed V maxNew reveals a steep increase since 1910 until 2014 in the Southwest Pacific (Figure 9).A period of 161 years between 1854 to 2014 shows steady fluctuations between 67 to 70 m/s wind speed.Climate 2016, 4, 30 13 of 19here and detailed analysis is suggested for this region.The most dominant cycle is of ~3 years which occurs around 1991, 1971 and 1891.

Figure 8 .
Figure 8. Wavelet power spectrum of VmaxNew in the Southern Indian Ocean.

( 1 )
The Southwest Pacific: Simulated TC wind speed VmaxNew reveals a steep increase since 1910 until 2014 in the Southwest Pacific (Figure 9).A period of 161 years between 1854 to 2014 shows steady fluctuations between 67 to 70 m/s wind speed.

Figure 9 .
Figure 9.Time series of Simulated TC wind speed VmaxNew in the Southwest Pacific.

Figure 9 .
Figure 9.Time series of Simulated TC wind speed V maxNew in the Southwest Pacific.

Figure 10 .
Figure 10.Time series of Simulated TC wind speed VmaxNew in the South Indian Ocean.

Figure 11 .
Figure 11.Time series of Simulated TC wind speed VmaxNew in the Northwest Pacific.

Figure 10 .
Figure 10.Time series of Simulated TC wind speed V maxNew in the South Indian Ocean.

( 3 ) 19 Figure 10 .
Figure 10.Time series of Simulated TC wind speed VmaxNew in the South Indian Ocean.

Figure 11 .
Figure 11.Time series of Simulated TC wind speed VmaxNew in the Northwest Pacific.

Figure 11 .
Figure 11.Time series of Simulated TC wind speed V maxNew in the Northwest Pacific.

Climate 2016, 4 , 30 15 of 19 Figure 12 .
Figure 12.Time series of Simulated TC wind speed VmaxNew in the North Indian Ocean.

Figure 13 .
Figure 13.Time series of Simulated TC wind speed VmaxNew in North Atlantic Ocean.

Figure 12 . 19 Figure 12 .
Figure 12.Time series of Simulated TC wind speed V maxNew in the North Indian Ocean.

Figure 13 .
Figure 13.Time series of Simulated TC wind speed VmaxNew in North Atlantic Ocean.

Figure 13 .
Figure 13.Time series of Simulated TC wind speed V maxNew in North Atlantic Ocean.

Table 1 .
Summary and intercomparison of tropical cyclone wind speed obtained from Emanuel's model and our approach with only SST as input.Comparisons are performed for three different SST inputs: COBE, ERSST.V3b and OISST.V2.Other parameters (SLP, specific humidity, atmospheric temperature profile) required for implementing Emanuel's model are from NCEP reanalysis datasets.