Recent studies have investigated trends and interannual variability in the potential intensity (PI) of tropical cyclones (TCs), but relatively few have examined TC PI seasonality or its controlling factors. Potential intensity is a function of environmental conditions that influence thermodynamic atmosphere–ocean disequilibrium and the TC thermodynamic efficiency—primarily sea surface temperatures and the TC outflow temperatures—and therefore varies spatially across ocean basins with different ambient conditions. This study analyzes the seasonal cycles of TC PI in each main development region using reanalysis data from 1980 to 2013. TC outflow in the western North Pacific (WNP) region is found above the tropopause throughout the seasonal cycle. Consequently, WNP TC PI is strongly influenced by the seasonal cycle of lower-stratospheric temperatures, which act to damp its seasonal variability and thereby permit powerful TCs any time during the year. In contrast, the other main development regions (such as the North Atlantic) exhibit outflow levels in the troposphere through much of the year, except during their peak seasons. Mathematical decomposition of the TC PI metric shows that outflow temperatures damp WNP TC PI seasonality through thermodynamic efficiency by a quarter to a third, whereas disequilibrium between SSTs and the troposphere drives 72%–85% of the seasonal amplitude in the other ocean basins. Strong linkages between disequilibrium and TC PI seasonality in these basins result in thermodynamic support for powerful TCs only during their peak seasons. Decomposition also shows that the stratospheric influence on outflow temperatures in the WNP delays the peak month of TC PI by a month.
It is important to understand the intensity of tropical cyclones (TCs) because coastal societies are vulnerable to their dangerous and costly impacts (e.g., Pielke and Pielke 1997). Studies over the last few decades have established a theoretical framework for investigating the potential intensity (PI) that a TC may reach given local environmental conditions (e.g., Emanuel 1985; Holland 1997; Bister and Emanuel 1998, 2002, hereafter BE02). Evidence suggests that variations in actual tropical cyclone intensities scale with variability in potential intensity on multiple time scales (Emanuel 2000; Wing et al. 2007). While it is recognized that warming sea surface temperatures (SSTs) have been and will continue to be a key driver for TC PI trends under anthropogenic climate forcing (e.g., Henderson-Sellers et al. 1998; Emanuel 2005; Knutson et al. 2010; Holland and Bruyère 2014; Strazzo et al. 2015; Walsh et al. 2016; and references therein), recent work has also focused on how lower-stratospheric temperatures influence TC PI (e.g., Emanuel et al. 2013; Ramsay 2013; Wang et al. 2014). Some studies have explored stratospheric linkages with TC PI trends over the past three decades (Emanuel et al. 2013; Wing et al. 2015), in the future under anthropogenic climate forcing (Vecchi et al. 2014; Sobel et al. 2016; Walsh et al. 2016), and in a Montreal Protocol world–avoided scenario (Polvani et al. 2016). Less attention has been given to interannual variability in TC PI and its connections with the middle atmosphere (Wing et al. 2015), and only a few studies have considered the seasonality of potential intensity (Free et al. 2004; Tonkin et al. 2000). In this study, we calculate the seasonal cycles of tropical cyclone potential intensity in the main TC development regions with 34 years of reanalysis data and use a decomposition method to determine the main factors that drive TC PI seasonality in these regions.
Tonkin et al. (2000) investigated the seasonal cycles of TC potential intensity with a set of radiosonde stations in the northwestern Pacific, Australian southwestern Pacific, and North Atlantic regions. Potential intensity calculations using the early Emanuel (1985, 1991) method and the Holland (1997) method were compared to one another and to observed TC intensities. Tonkin et al. (2000) found that both methods generally performed well at predicting observed seasonal intensities in the tropics (though each was overly sensitive to certain environmental conditions). However, their study did not directly attribute the seasonal cycles to local thermodynamic contributions beyond SSTs, and their data were spatially limited. Free et al. (2004) also explored potential intensity seasonal cycles. Using data from 14 radiosonde stations, they showed that there were seasonal differences in TC PI across the tropics. In the Caribbean these differences were attributed to each station’s distance from the equator (with calculations showing larger seasonal cycle amplitudes at stations farther from the equator), whereas no explanation was offered for the Pacific basin differences. Further, no analysis was performed to determine the thermodynamic conditions contributing to the seasonal cycles in each ocean basin. In this study, we will show that spatial differences in the seasonal cycles of TC PI are linked to the seasonality of each ocean basin’s sea surface temperatures, TC outflow levels, and outflow temperatures. Note that we do not explore the seasonality of other factors known to be important for the genesis or intensity of tropical cyclones, such as vertical wind shear (e.g., Frank and Ritchie 2001; Emanuel 2006).
The goals of this study are to characterize the seasonal cycles of TC PI pertaining to the main TC development regions in each ocean basin and to identify the key factors contributing to these cycles. In section 2, we describe the TC PI calculation and decomposition, along with the reanalysis data used to compute TC PI. Section 3 shows the resulting TC PI seasonal cycles and quantifies the thermodynamic contributions to these cycles. The study is summarized in section 4.
2. Data and methods
Tropical cyclone potential intensity is defined following Bister and Emanuel (1998):
where Ck and CD are the enthalpy and momentum surface exchange coefficients, respectively; Ts is the sea surface temperature; To is the mean outflow temperature; is the saturation moist static energy at the sea surface; and is the moist static energy of the free troposphere. The ratio is an uncertain quantity and determining its value is an area of active research (e.g., Emanuel 2003, and references therein), but it is typically taken to be a constant. Here we set following the precedent of previous studies (Wang et al. 2014; Wing et al. 2015), but we have tested a range of reasonable ratios (0.4–1.0) and found our results to be qualitatively insensitive to the choice. The term is the thermodynamic efficiency of the potential tropical cyclone. It is readily determined that kelvin for kelvin the efficiency term is more sensitive to than . However, seasonal amplitude is generally larger than , and and may be related to one another through the outflow level, with key implications for the seasonal cycles of TC PI (see section 3). The final term is the thermodynamic disequilibrium between the sea surface and the free troposphere, which is principally controlled by sea surface temperatures (e.g., Emanuel 2007).
We compute the potential intensity using the algorithm of BE02. The BE02 scheme is more accurate for TC PI calculations in the tropics than directly calculating and summing each individual term on the right-hand side of Eq. (1) because of its incorporation of water loading and virtual temperature (Garner 2015; Wing et al. 2015). Besides its increased accuracy over direct computation of Eq. (1), the BE02 algorithm is advantageous because it has been frequently employed in previous studies of PI (e.g., Emanuel 2007; Sobel and Camargo 2011; Emanuel et al. 2013; Wang et al. 2014; Wing et al. 2015; Polvani et al. 2016), allowing straightforward comparisons between those studies and this work. We computed TC PI directly with Eq. (1) and found that the thermodynamic disequilibrium term is larger than that inferred from the BE02 algorithm by 10%–25%, but results are qualitatively similar to those in section 3.
The BE02 algorithm requires profiles of temperature and humidity, along with sea surface temperatures and mean sea level pressures. The algorithm is configured to output , the outflow temperature, and the outflow temperature level (OTL). The OTL is the level of neutral buoyancy for an air parcel saturated at sea level pressure, corresponding to the outflow temperature. Outflow temperatures allow direct calculation of the efficiency term [second term on the right-hand side of Eq. (1)] at each temporal and spatial location. Given the output potential intensity and efficiency, the thermodynamic disequilibrium can be derived from Eq. (1). In the BE02 algorithm, parcels are lifted assuming reversible adiabatic ascent. Assuming pseudoadiabatic ascent instead leads to qualitatively similar results, except that OTLs are typically found at higher altitudes, allowing more stratospheric influence on TC PI in all regions (see section 3). We decompose Eq. (1) by taking the natural logarithm of both sides:
This log-additive model is advantageous because it quantifies thermodynamic disequilibrium and thermodynamic efficiency contributions to TC PI at every spatial and temporal location.
Temperature, water vapor, sea surface temperature, and mean sea level pressure data from 1980 to 2013 are taken from the second Modern-Era Retrospective Analysis for Research and Applications (MERRA-2; Bosilovich et al. 2016; NASA 2016) and the European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-Interim, hereafter ERA-I; Dee et al. 2011; ECMWF 2016). The 6-hourly data from both reanalyses are regridded onto a monthly 2.5° × 2.5° grid. As we average over tropical cyclone main development regions (see below), the horizontal resolution is not expected to play a significant role in the results. Previous work by Kossin (2015) validated observed TC cloud-top brightness temperatures against MERRA/ERA-I outflow temperatures and found that they vary consistently with one another on multiple time scales, suggesting that seasonal evaluation of TC outflow temperatures and potential intensity with these reanalyses is appropriate. For clarity and brevity we show only the MERRA-2 results herein; ERA-I results are largely consistent with MERRA-2 and are shown in the supplemental material (Figs. S1–S4).
Our methodology is designed to evaluate the seasonal cycles of potential intensity dictated by seasonal atmospheric and oceanic states. We determine the monthly means of each state variable from reanalyses and compute TC PI with the BE02 algorithm at every ocean grid point. Monthly means of state variables and algorithm outputs are then time averaged over the 34-yr period and area averaged over the tropical cyclone main development regions in each ocean basin (Table 1) following Emanuel (2005) and Wing et al. (2015), except that we retain the full seasonal cycle in each case (rather than averaging over the peak months) and use a somewhat different NA main development region definition than Wing et al. (2015). Our NA region focuses specifically on TC PI in the far-western Atlantic and Gulf of Mexico, where tropical cyclones frequently develop (e.g., Gray 1984); results for the NA are qualitatively insensitive to whether we use this definition or that of Wing et al. (2015). Furthermore, results for each main development region in Table 1 are qualitatively robust to modest (~5°) poleward shifts (see Fig. 2, described in greater detail below).
We do not account for interannual variability in our results [e.g., ENSO, quasi-biennial oscillation (QBO), and volcanoes], but the 34-yr period averages are expected to smooth such variability and demonstrate typical TC PI seasonal cycles in each ocean basin. To test the sampling uncertainty in the mean seasonal cycles we employ an empirical bootstrap method (e.g., Efron and Gong 1983). Bootstrapping is used to construct 1000 simulations of 34-yr subsets that are resampled from 1980 to 2013 MERRA-2 data with replacement. For each simulation we recalculate the mean seasonal cycles of SSTs and algorithm outputs (i.e., TC PI, outflow temperature, and OTL) and determine their 95% confidence intervals in each region and month. Results from this exercise show that each of these seasonal cycles averaged over the main development regions is statistically robust (see Fig. S7 and Table S1).
a. Seasonal cycle overview
We first provide an overview of the seasonal cycles of TC PI and associated atmospheric conditions in each main development region. Figure 1 shows the seasonal cycles of sea surface temperature, outflow temperature, the outflow temperature pressure level, and TC PI for each main development region from MERRA-2 data. The seasonal cycles calculated with ERA-I data are similar to those of MERRA-2 (see Figs. S1 and S2). Figure 2 shows the horizontal structures of the monthly peak-to-peak seasonal cycle amplitudes (hereafter amplitudes) of SSTs, outflow temperature, and TC PI from MERRA-2. Positive amplitudes in Fig. 2 indicate that the seasonal cycle maximizes in the boreal summer and negative amplitudes indicate that the seasonal cycle maximizes in the boreal winter; boxes outlining each main development region are also included (except NI, see below).
SST variability is driven by complex interactions in the coupled ocean–atmosphere system (e.g., Deser et al. 2010). SST seasonal cycles in the NA, ENP, WNP, and SH regions (Fig. 1a) each have a single annual maxima and a single annual minima, which coincide with their hemispheric summer or winter, respectively. The phase of the SST seasonal cycle is delayed by several months relative to the annual solar cycle owing to atmosphere–ocean thermal inertia (e.g., Cronin and Emanuel 2013). In the NI region, the monsoonal circulation drives a semiannual seasonal cycle in SSTs (Schott and McCreary 2001; Schott et al. 2009; and references therein). Because our amplitude definition does not account for a semiannual component, we do not consider the NI region in our seasonal amplitude analyses (e.g., Fig. 2; Table 2).
In the WNP, SSTs have limited seasonality (Fig. 1a) and are above the canonical 26.5°C threshold for TC development year-round (e.g., Palmén 1948; although this requisite has been recently reexamined, and its specific value is expected to increase with climate change; e.g., Vecchi and Soden 2007; Johnson and Xie 2010; McTaggart-Cowan et al. 2015). The deep mixed layer in the Pacific warm pool region results in a weak SST seasonal cycle (mean amplitude of 1.5 K) at the lower tropical latitudes (Fig. 2a; Wyrkti 1961; Schneider and Zhu 1998). The ENP region also has weak seasonality (amplitude of 1.3 K), but SSTs are cooler by 1–1.5 K (Fig. 1a) and have a less homogenous seasonal cycle across the region (Fig. 2a). The NA and SH regions both have larger seasonal amplitudes in SST (3.8 and 2.3 K, respectively) that are spatially inhomogeneous and associated with seasonally varying mixed layer depths (de Boyer Montegut et al. 2004).
Figure 3 shows the annual mean OTLs in the tropics (Fig. 3a), along with the seasonal cycles of the OTL pressure minus the lapse rate tropopause pressure averaged over each main development region (Fig. 3b). Here the lapse rate tropopause is defined as the lowest height where the lapse rate falls below 2 K km−1 (World Meteorological Organization 1957) and was determined with MERRA-2 data following the methodology of Reichler et al. (2003). A positive difference in Fig. 3b indicates that the monthly mean OTL is at a higher pressure (lower altitude) than the monthly mean tropopause. Note that the height of the tropical tropopause also has a consistent seasonal cycle: it is found at minimum pressures in the boreal winter and maximum pressures in the boreal summer (e.g., Kim and Son 2012). For reference, Fig. 4 shows the seasonal anomalies of air temperature (with the annual mean at each pressure level removed) in tropical troposphere and lower stratosphere, the lapse rate tropopause pressure, the outflow temperature levels, and the seasonal cycles of SSTs in each main development region.
The seasonal cycles in outflow temperature (Fig. 1b) are tightly linked to outflow temperature levels (Figs. 1c and 3). While the OTL remains in the troposphere, lower pressures are associated with cooler outflow temperatures (as temperatures decrease with height up to the cold point). Although there are temperature seasonal cycles in the tropical troposphere, their amplitudes on any tropospheric pressure level are small (<4 K; Fig. 4; cf. Donohoe and Battisti 2013), and instead the lapse rate and OTL pressure dominates outflow temperature seasonality. Recall that the OTL is the level of neutral buoyancy for a saturated air parcel lifted from sea level and is hence thermodynamically coupled to local SSTs (section 2; Fig. 4): warmer SSTs will be associated with higher (and therefore cooler) outflow levels, whereas cooler SSTs with be associated with lower (and warmer) outflow levels. Outflow temperatures and SSTs are therefore generally opposite phased when outflow is found in the troposphere (e.g., in the NA region from November to May; Fig. 3b). This inverse relationship has important implications for TC PI seasonality because the difference between SSTs and outflow temperatures determines the thermodynamic efficiency of the TC PI [Eq. (1)]. When the OTL migrates across the tropopause into the stratosphere, however, outflow temperatures are regulated by the lower-stratospheric temperature seasonal cycle. In each region OTLs penetrate the stratosphere at some point in the seasonal cycle, typically during their periods of largest PI (Figs. 3b–d and 4).
Annual mean OTLs in the WNP (Fig. 3a) are largely homogeneous over the region and are always found in the stratosphere (<95 hPa; Figs. 3b and 4c). This is consistent with the annual presence of the western Pacific warm pool, which “anchors” the ITCZ perennial deep convection in the west Pacific (Wytki 1961, 1989; Mitchell and Wallace 1992). In contrast, annual mean OTLs over the NA vary considerably from ~90 hPa in the southwest part of the region to ~150 hPa in the northeast. Consistent with the observed seasonal shifts of the ITCZ (Mitchell and Wallace 1992), the NA OTL only penetrates the stratosphere climatologically during the boreal summer and fall months (Figs. 3b and 4a; i.e., during the peak NA hurricane season). The other main development regions (ENP, NI, and SH) exhibit similar annual OTL inhomogeneity and seasonal variance, with their OTLs penetrating the stratosphere primarily during their peak-TC seasons (Fig. 3b). Notably, observed tropical cyclones are important and disproportionately large contributors to annual overshooting convection in the tropics, especially in the WNP (Romps and Kuang 2009). Our results suggest that TCs attaining their potential intensity would exhibit similar overshooting behavior with regionally dependent seasonal variance (Fig. 3b).
The seasonality of NA OTLs results in a large outflow temperature seasonal cycle amplitude (10.9 K) that maximizes in the boreal winter and minimizes in the boreal fall (Fig. 1b), with a considerable gradient across the basin (Fig. 2b). The ENP and SH both have seasonal cycles of outflow temperature that are inversely related to their seasonal SSTs when their outflow is found in the troposphere, and like the NA these effects are not spatially homogenous. The semiannual cycles of NI SSTs also tracks inversely with regional average outflow temperatures when its outflow is found in the troposphere. In each region, the outflow temperature seasonal amplitude is dominated by the vertical migrations of their OTLs in the troposphere.
The perennial penetration of the WNP OTL into the stratosphere contrasts with the other regions and has important implications for WNP TC PI seasonality. By “seeing” the lower stratosphere year-round, WNP seasonal outflow temperatures are set by the factors controlling tropical lower-stratospheric temperatures. The seasonal cycle of tropical lower-stratospheric temperatures is annual rather than semiannual (Reed and Vlcek 1969), maximizing in the boreal summer and minimizing in the boreal winter, and is driven primarily by planetary wave forcing from the troposphere (e.g., Reid 1994; Yulaeva et al. 1994; Fueglistaler et al. 2009, 2011). Enhanced wave driving in the Northern Hemisphere compared with the Southern Hemisphere (SH) results in an annual cycle of upwelling, adiabatically cooling lower-stratospheric temperatures in the boreal winter, and warming them in the boreal summer. Recent studies have shown that this seasonal cycle is modified radiatively by the correlated upwelling of ozone, which amplifies the annual cycle by ~30% (Folkins et al. 2006; Chae and Sherwood 2007; Fueglistaler et al. 2011; Gilford and Solomon 2017). The net result is an annual cycle across the tropical lower stratosphere with an amplitude peaking at ~8 K (e.g., Reid 1994; Randel and Wu 2015). Although this seasonal cycle does not vary much in the zonal average between the two hemispheres (Fig. S5), the largest seasonal amplitudes of lower-stratospheric temperatures are found over the western Pacific and Caribbean regions, with important implications for tropical cyclones (Fig. S6; Fig. 4; see section 3b). Because the outflow of WNP TCs is always within the stratosphere (Fig. 3b), the outflow temperature seasonal cycle follows that of the lower stratosphere, maximizing in the boreal summer and minimizing in the boreal winter (Figs. 1b and 4c).
Like the WNP, the other basins show signs of lower-stratospheric influences on their monthly mean outflow temperatures in those months when their OTLs penetrate into the stratosphere (Fig. 4). Outflow temperatures in the Northern Hemisphere main development regions show similar curvature to the WNP outflow temperatures from August through November, when outflow temperatures are falling from ~198 to ~194 K following the stratospheric seasonal cycle (Figs. 1b and 4). The NA region being more poleward than the other main development regions (Table 1) results in warmer average outflow temperatures during the boreal summer and fall, delaying the NA outflow temperature minima until October (Fig. 1b). NI outflow reaches the stratosphere in both the boreal spring and fall, incidentally missing the seasonal maximum in lower-stratospheric temperatures (Fig. 4d) and resulting in an inverse relationship between SSTs and outflow temperatures in these seasons (Fig. 1b). In the months that SH OTLs are found in the stratosphere, from November through May (Figs. 3b and 4e), the magnitude and shape of the SH outflow temperatures resemble those in the WNP (Fig. 1b) and follow the lower-stratospheric seasonal cycle (Fig. 4e).
TC PI seasonality reflects regional differences in sea surface temperatures and outflow temperatures (Fig. 1d). The NA TC PI seasonal cycle has a large amplitude (~22 m s−1) and maximizes in September/October, coinciding with a large amplitude seasonal cycle in SSTs. The TC PI seasonal cycles in the ENP are spatially inconsistent within the region (Fig. 2c) but on average show similar structure to their SSTs (cf. Figs. 1a,d). The SH seasonal cycle amplitude is similar to the NA (~23 m s−1), maximizing in the austral fall. The NI seasonality highlights the important influence of the monsoonal circulation in setting TC potential intensity in this region. In the WNP, TC PI is strong year-round (always >72 m s−1) and has a small seasonal amplitude (~6 m s−1) relative to the other regions. This small amplitude indicates that either the disequilibrium term has similar weak seasonality, the thermodynamic efficiency term is damping TC PI because of stratospheric influences, or both. The relative contributions from these terms are explicitly investigated in section 3b.
In contrast to the proposal of Free et al. (2004), Fig. 2c shows that TC PI seasonal cycles are not simply zonal in nature. For instance, although the ENP and WNP have similar latitudinal ranges (Table 1), they have distinct seasonal cycle structures in part because of zonal asymmetries in SSTs and outflow temperatures (Figs. 2a,b).
The TC PI seasonal cycles in Fig. 1d compare generally well with those found by Tonkin et al. (2000). In the NA, both studies find large amplitude TC PI seasonality, with values maximizing during the peak-season months (Fig. 1d). Their TC PI seasonal cycles in the low-latitude Pacific stations (equatorward of 15°N) also agree with our finding of a weak WNP TC PI seasonal cycle. Tonkin et al. (2000) also show good agreement between observed TC intensities and TC PI, and taken together our results are illustrative of why powerful TCs have been observed at all times of year in the lower latitudes of the western Pacific basin (particularly shown by Guam station data; 14°N, 145°E; their Fig. 5).
b. Potential intensity decomposition
We have shown that the seasonal cycle of WNP outflow temperatures is perennially linked to lower-stratospheric temperatures, whereas outflow temperatures in the other regions are found in the troposphere except during their peak seasons. These differences are broadly consistent with the derived TC PI in each main development region (Fig. 1d). We now employ the decomposition in Eq. (2) to specifically quantify to what extent these differences affect TC PI seasonality. Figure 5 shows the seasonal anomalies (annual means removed) of the decomposed terms averaged over the main development regions. The monthly variations in disequilibrium and efficiency terms show their specific monthly contributions to the TC PI magnitude, and these seasonal cycles are statistically robust to sampling uncertainty (Table S2). Note that is a constant that does not contribute to seasonal variability. Table 2 shows the TC PI seasonal cycle amplitudes (from Fig. 1d) in each region excluding the NI region, along with the percent changes in the TC PI term amplitude if the efficiency or disequilibrium term seasonal cycles are subtracted from the TC PI term seasonal cycle. A negative percent change indicates that the term climatologically amplifies the TC PI seasonal cycle, whereas a positive percent change indicates that the term climatologically damps the TC PI seasonal cycle. Percent changes do not add to 100% because of seasonal phase differences between the terms (discussed below). We see from Table 2 that in every region the disequilibrium term is the largest overall driver of TC PI seasonality and hence always amplifies TC PI seasonal cycles.
The disequilibrium and efficiency terms in the NA (Fig. 5a) both maximize in the boreal summer/fall and minimize in the boreal winter/spring. SSTs and outflow temperatures both overall amplify the efficiency seasonality in the NA (Figs. 1a,b; Table 2), but the disequilibrium term—driven principally by SSTs—is a much larger contributor to TC PI (Table 2). Consequently, NA TC PI has strong seasonality (with an amplitude that is 35% of its annual average value) linked to the SST seasonal cycle, whereas stratospheric temperatures play a secondary role. The main influences of the NA stratospheric temperatures (when OTLs reach the stratosphere from June to October; Fig. 3b) are to slightly reduce the magnitude of the efficiency term (Fig. 5a) and slightly delay the peak month of the TC PI seasonal cycle (cf. the peak month between the NA disequilibrium term and the TC PI term). Overall the NA environment is thermodynamically most conducive for intense TCs in the boreal summer/fall, consistent with the observed NA hurricane season.
Results for the SH and ENP regions generally agree with the NA: SSTs dominate the seasonal cycle through the disequilibrium term (Figs. 5b,e), and the efficiency term is a minor amplifier (Table 2). Note that because the seasonal cycle of lower-stratospheric temperatures is similar across the tropics in both hemispheres (Fig. 4; Fig. S5), the seasonal cycle of stratospheric outflow temperatures in the SH (over November–May; Fig. 3b) amplifies the efficiency seasonal cycle (and hence TC PI) in that region rather than damping it (Table 2; Fig. 5e). TC PI in the NI region is strongly coupled to the disequilibrium term, and the efficiency term only slightly amplifies its seasonal variability (Fig. 5d).
WNP disequilibrium and efficiency terms have opposing signs and more comparable magnitudes (Fig. 5c; Table 2). The phase of the efficiency term seasonal cycle, maximizing in January–February and minimizing in August–September, is indicative that this term is principally controlled by lower-stratospheric outflow temperatures. Contributions from the disequilibrium term maximize in September–October and minimize in February, following regional SSTs (Fig. 1a). The weak seasonality of WNP SSTs (amplitude of 1.5 K) coupled with the relatively strong seasonality in lower-stratospheric outflow temperatures (amplitude of 6.4 K) allows efficiency impact to TC PI seasonality on the order of thermodynamic disequilibrium’s impact (Table 2). The efficiency term damps the amplitude of the WNP TC PI seasonal cycle by about 26%. Because efficiency and TC PI are seasonally out of phase, the efficiency term’s influence on PI is in fact larger than the amplitude metric shows. If, for example, the efficiency term was phase shifted by two months so that the minima in efficiency aligned with the maxima in TC PI, the WNP efficiency term would damp the TC PI seasonal cycle by 39%.
The phase offsets between the efficiency and disequilibrium terms in the WNP depress TC PI in the boreal summer. Because SSTs lagging the seasonal insolation (e.g., Braconnot et al. 2000) are offset from the asymmetric wave forcing on lower-stratospheric temperatures (e.g., Yulaeva et al. 1994), the WNP peak season occurs later than would be expected if SSTs were the only driver of WNP TC PI seasonality (shifting its TC PI maxima from September into October; Fig. 5c). The positive contribution of the NA efficiency term also delays the TC PI peak in that region, resulting in nearly equal monthly averages of NA TC PI in September and October.
To gauge the relative importance of each term in Eq. (2) between the main development regions, we examine the ratios of their amplitudes (see Table S2). If TC PI seasonality is influenced by the same factors in two regions, contributions from each term should scale consistently between those regions. Here we take the ratio of the NA and WNP term amplitudes, as they show the most dramatic contrasts in how the stratosphere impacts TC PI through the efficiency term (the NA region is characteristic of the behavior in the ENP, NI, and SH regions). We find that the amplitude ratio of NA/WNP TC PI terms is 3.8, whereas the amplitude ratios of the disequilibrium and efficiency terms are 2.2 and −2.6, respectively. Although the absolute magnitudes of their influences are different (e.g., vertical axes of Fig. 5), the differences between their proportionalities make clear that the NA seasonality (and by extension the ENP, NI, and SH regional seasonality) is strongly driven by SST seasonality through both thermodynamic disequilibrium and efficiency, whereas the WNP seasonal cycle is effectively damped by stratospheric influences on thermodynamic efficiency.
The seasonal cycles of tropical cyclone (TC) potential intensity (PI) have been calculated with 34 years of reanalysis data across five TC main development regions. Potential intensity depends not only on sea surface temperatures driving thermodynamic disequilibrium with the atmosphere but also on the difference between the sea surface temperatures and TC outflow temperatures (e.g., Emanuel 1985, 2003). Whereas previous studies have focused primarily on the long-term trends and interannual variability in TC PI, in this study we have elucidated the seasonality of PI and determined the relative seasonal contributions of thermodynamic disequilibrium and thermodynamic efficiency with a decomposition method.
Environmental conditions in the western North Pacific imply that the region has perennial outflow levels that penetrate into the stratosphere and outflow temperatures that are stratospherically controlled. In contrast, outflow temperatures in the North Atlantic main development region are generally found at lower and warmer tropospheric altitudes except during the boreal summer and fall seasons (from June through October). The net result is a large-amplitude (~22 m s−1 from peak to peak) TC PI seasonal cycle in the North Atlantic that maximizes in the Atlantic hurricane season, compared with the weaker-amplitude but perpetually powerful (> 70 m s−1) TC PI seasonal cycle in the western North Pacific. Although SSTs can be thought of as dictating the seasonal cycles in all basins, in the western North Pacific they do so in part by allowing TC outflow to be influenced by the stratosphere year-round, whereas in the other main developments regions SSTs are linked more directly to tropospheric outflow temperatures through seasonal variations in the outflow level.
Decomposing the TC PI metric reveals that thermodynamic disequilibrium is the main driver of TC PI seasonal cycles, contributing 72%–85% of the seasonal variability in each region except the western North Pacific. Thermodynamic efficiency contributions in these regions are smaller and amplify the seasonal cycle of TC PI because they are tropospherically controlled (through SSTs) during most of the year. In the western North Pacific, however, efficiency damps TC PI by between a quarter and a third because TC outflow temperatures there are linked to the seasonal cycles of lower-stratospheric temperatures, which enhance TC PI in the boreal winter and decrease it in the boreal summer. Furthermore, the stratospheric influence on western North Pacific outflow temperatures delays the peak month of TC PI by about a month relative to the peak implied by SST seasonality alone.
Our results imply that thermodynamic conditions in the western North Pacific are ripe for intense TCs throughout the year, whereas ambient conditions in the North Atlantic region typically support powerful TC intensities only during the peak hurricane season (late boreal summer and fall). In this study we have only assessed the thermodynamic conditions for intense tropical cyclones. Although powerful typhoons have been observed year-round in the western Pacific, there remains clear seasonality in their frequency (Sopko et al. 2014). This distinction is likely owing to other factors such as wind shear and midtropospheric humidity, which are not accounted for in our methodology.
The lead author thanks Allison Wing, Paul O’Gorman, Qiang Fu, and Tom Beucler for helpful discussions. We thank the anonymous reviewers whose suggestions improved this work. DG was supported by NASA Headquarters under the NASA Earth and Space Science Fellowship Program by Grant NNX14AK83H. SS, KE, and DG were supported by NSF Grant AGS-1461517.
Supplemental information related to this paper is available at the Journals Online website: http://dx.doi.org/10.1175/JCLI-D-16-0827.s1.