More than unfamiliar environmental connection to super typhoon climatology

More than unfamiliar environmental connection to super typhoon climatology | Scientific Reports

Abstract

This study employs a refined geometric variability model to look at the environmental relationship to super typhoon climatology, which is one of the major concerns about climate change and disasters. It is noted that adding only several recent years leads to a remarkable weakening of the environmental explanatory power on super typhoon climatology. Looking into the annual covariance elements, we find that the recent observations showing a group of outlying events with a particular drift are more than unfamiliar compared to the former stable relationship from 1985 through 2012. Greater uncertainty thereby amplifies concerns about the looming climate crisis.

Introduction

The Earth’s climate is configured by air, water, ice, rocks, and living things all of which change over time 1. If the forcings external to the climate system remain static, there is only internal variability. In the real world the climate system is stimulated by external forcings such as the Earth’s orbit, solar activity, volcanic activity, and aerosol concentrations, and so on 2. During the recent century, increasing concentration of greenhouse gases, despite their contribution to urbanizing society, is exerting a force on the earth’s climate 3. As a response, record-breaking extremes are continually and everywhere being reported creating a climate crisis 4, 5. The larger responses, even at the extremes, do not necessarily mean the climate is out of balance (Supplementary Fig. S2). In the study, “balance” denotes the stable relationship among the climate conditions, and potentially more dire is the possibility that an imbalance could appear as an unusual climatic response 6, 7, 8.

Tropical cyclones (TCs) create enormous socio-economic costs 6, 9, 10, 11 but also contribute to the earth’s energy distribution 12, 13, 14, 15. Because of this role, the relationship between TC activity and the environment is supposed to show a certain aspect of the functional balance. This study employs a refined geometric variability model to look at the environmental relationship to super typhoon activity, which is one of the major concerns about climate change and disasters. Super typhoon climatology refers to a set of observations including total genesis frequency, rate of super typhoon occurrences, overall TC activity, and the efficiency of intensification 16. While the climatology has shown a strong relationship with the large scale forcings of El Niño-Southern Oscillation (ENSO) and global ocean warmth 17, we find more than unfamiliar connections beginning in the early 2010s.

A geometric variability model. The climate connection between the two variability planes is configured by rotation (\(\theta _{1}\)), scaling (r) in the maximum covariance direction (\(\theta _{2}\)), and tilting (\(\theta _{3}\)) of the response variability plane. The geometric model gives a presentation of the response variability plane onto the explanatory variability plane in the three-dimension variability space. All variables use standardized values. The center of the planes represents the mean value, while the boundary indicates 1.0. The response variability plane and the explanatory variability plane are constructed from the TC variables and the environmental variables, respectively.

Reduction in the explanatory power of the environmental variables. Comparison of the environmental relationships to TC variables between two periods. Variables are averaged over June to November. The former 30-year (1985–2014) TC-climate connection to the environment is shown in blue, while the recent update (1985–2020) is shown in red. The positions of the TC variables show the best explanatory variability and its explanatory power. A large reduction in the explanatory power together with the changes in the corresponding directions is noted when including data from the years since 2013.

Sudden and still ongoing drift of the climate connection. The direction and magnitude of the annual maximum covariance elements (\(mX_{t} \cdot mY_{t}\)) distributed on (A) a time domain and (B) an XY plane, and (C) the interruption (\(mX_{t} \cdot mY_{t} -mX_{t} \cdot mX_{t}\)). All variables are annually averaged from June through November. Since the variability direction is cyclic, the distribution is tiled as a mosaic. It is noted that \(\theta _{1}\) had been stable around the value of − 51.9\(^{\circ }\) until 2012 (sky-blue lines), and the sudden and still ongoing drift of the climate connection began in 2013. Simultaneously, the interruption started a distinct fluctuation in 2014 (red line).

Time series of the atmospheric suppression and the efficiency of TC intensity in the western North Pacific. The atmospheric suppression is indicated by geopotential height at 500 hPa. The efficiency of TC intensity is the out-of-phase principal component of annual TC intensity and frequency. All values are annually averaged over JJASON in the tropical region of the western North Pacific (0\(^\circ\)–30\(^\circ\)N, 100\(^\circ\)E–180\(^\circ\)). Green vertical line indicates the year 2012. For demonstration, annual values are standardized by the values during the 28 years (1985–2012).

A geometric variability model

By employing a ‘three-dimensional variability space’ 17, we refined a geometric variability model to quantify how the response variability space harmonizes with explanatory variability space and how that harmony is changing (Fig. 1). For this, the two variability planes are configured and investigated. The model includes components describing the shape of the response variability plane (\(P1_{Y}\) \(\sim\) \(P2_{Y}\) \(\sim\) − \(P1_{Y}\) \(\sim\) − \(P2_{Y}\) \(\sim\) \(P1_{Y}\)) juxtaposed above the components describing the explanatory variability plane (\(P1_{X}\) \(\sim\) \(P2_{X}\) \(\sim\) -\(P1_{X}\) \(\sim\) -\(P2_{X}\) \(\sim\) \(P1_{X}\)). The principal component analysis produces the same number of principal components as the number of input variables. Here, \(P1_{Y}\) and \(P2_{Y}\) are the two principal components from Y1 and Y2. \(P1_{Y}\) and \(P2_{Y}\) respectively show the in-phase and out-of-phase modes between Y1 and Y2. The variability space is constructed with perpendicular axes for illustration but that does not imply the variables Y1 and Y2 are orthogonal. The same for \(P1_{X}\) and \(P2_{X}\) from X1 and X2.

The overall shape of the response variability plane is constrained by (1) rotation (\(\theta _{1}\)), (2) scaling (r) in the maximum covariance direction (\(\theta _{2}\)), and (3) tilting (\(\theta _{3}\)) of the response variability plane. mX and mY denote each maximum covariance direction for the X-plane and Y-plane. In this variability framework, a correlation is the projection length between the variables. All variables are standardized by subtracting the mean and dividing by the standard deviation using observations over the entire period. The center and the boundary in the screen show the mean of the annual values and a correlation coefficient of one, respectively. The larger the angle between the directions defined by the variables, the weaker the correlation between them. The final projection of the response variability plane onto the explanatory variability plane shows a bird’s-eye view of how the environment is statistically connected to the TC climate, which effectively integrates all regression coefficients from individual ordinary-least-squares methods (Y \(\sim\) X) (Supplementary Fig. S5).

Two variability planes

To identify features of the TC response to environmental climate conditions, we focus on super typhoons in the western North Pacific during the period June through November (JJASON) because they are the most energetic. Considering the reliability of the best-track data 18, 19, 20, we use observations over the period 1985–2020, inclusive (36 years). A TC has a lifetime-maximum intensity (LMI) exceeding 17 m \(\hbox {s}^{-1}\). A super typhoon has an LMI of at least 65 m \(\hbox {s}^{-1}\) and we use intensities estimated by the US Joint Typhoon Warning Center.

First, super typhoon climatology refers to a set of observations including total genesis frequency, rate of super typhoon occurrences, overall TC activity, and the efficiency of intensification 16. Here, we configure a TC variability plane using the annual frequency defined by the annual number of TC occurrences and the annual TC intensity defined by the annual proportion of TCs that become super typhoons S6 A). Other choices for annual TC intensity include the annual mean LMI and 70th percentile intensity 21 but the number of Category 3, 4 or 5 typhoons as an inner product of the total number of TCs and the proportion of strong typhoons (i.e., the frequency and the intensity) is statistically comparable to other metrics of TC activity such as the Accumulated Cyclone Energy 22 and the Power Dissipation Index 23. With the intensity and the frequency denoted as Y1 and Y2, respectively, the two principal components of in-phase mode (\(P1_{Y}\)) and out-of-phase mode (\(P2_{Y}\)) are obtained. Here, the subscript “Y” in the variable names denotes they are the response variables to the environment. The principal components \(P1_{Y}\) and \(P2_{Y}\) are interpreted as the annual TC activity and the efficiency of intensification, respectively. The efficiency of intensification is the intensification at the expense of TC frequency 17, 24, and thus, the larger \(P2_{Y}\) implies a larger proportion of super typhoons even with fewer TCs. The variability plane is helpful to locate and compare any available TC variables in a continuous variability space.

Similarly, we configure an environmental variability plane using the global mean sea surface temperature (GMSST) to indicate the forcing by global ocean warmth, and the Southern Oscillation Index (SOI) to indicate forcing by oceanic-atmospheric internal variability through ENSO (Supplementary Fig. S6 B). These environmental variables are denoted as X1 and X2, respectively. Since a positive value of SOI reflects the La Niña pattern of ENSO, we use negative SOI as the positive sign on X2 to indicate ENSO variability. The environmental variability plane is defined by the two principal components of in-phase mode (\(P1_{X}\)) and out-of-phase mode (\(P2_{X}\)). The subscript “X” represents the explanatory variables. \(P1_{X}\) indicates the in-phase variability direction of a warmer El Niño condition, while \(P2_{X}\) indicates the out-of-phase variability direction of a warmer La Niña condition 17, 24. Similarly, a colder La Niña condition and a colder El Niño condition are in the opposite directions of \(P1_{X}\) and \(P2_{X}\), respectively. With the help of this variability space, the directions of the environmental variables are identified by their projection onto the variability plane. The X2 appears orthogonal to X1 (\(r= +0.10; [-0.24, 0.41]\), 95 % confidence interval), and other variables found in the literature for describing the TC environment project onto the plane in directions similar to their counterparts.

The possibility exists that a longer-term internal variability might seem to have some trend in a shorter time window. For instance, periodical solar activity can modulate Pacific meridional mode (PMM) variability 25, and subsequent responses could occur in the tropical SST through the wind-evaporation-SST process 26. PMM variability could also be influenced by external factors such as greenhouse warming as well as the volcanic eruptions 27. Regarding global warming, Niño indices such as Niño 1+2, Niño 3, Niño 3\(\cdot\)4, and Niño 4, as defined by tropical SSTs, cannot be free from environmental changes in the tropics 28. Changes to the tropical SST pattern like those across the central Pacific SST indicated by the Niño 4 index can excite a PMM response through the Aleutian low variability exhibiting the Pacific decadal oscillation pattern 29. Consequently, an environmental variable between X1 and X2 implies the merged reflection of global ocean warmth and ENSO on interannual timescales. To the extent that such a variable is associated with ocean dynamics, its response to the forcing would be slow due to thermal inertia of the ocean 30. This would weaken the correlation between the longer-term and interannual variability. The orthogonal variability plane on interannual timescales effectively spans the environmental variability directions onto which the TC variables can be projected.

Remarkable changes to the recent climate connection

The environmental connection to TC activity is described by associating the two corresponding variability planes. The correlation wavelets shown in Fig. 2 are understood as the expression of the directions and the lengths from this geometric projection (Supplementary Fig. S5). The thirty-year (1985–2014) connection is revisited 17 as a prior reference (Fig. 2 blue line). We define the connection by the direction of the closest environmental variability and the explanatory power (\(r^{2}\)) to each TC variability. Red and blue dots in Fig. 31. Physically, TC genesis frequency (Y2) in a warmer environment is hindered by active moist adiabatic warming in the troposphere 32, 33. Suppression (− Y2) is most effective when La Niña conditions occur during a warmer year (\(P2_{X}\)) which is likely to be accompanied by decreasing upward mass flux 34, 35, increasing saturation deficit 36, 37, and anomalous highs in the middle and upper troposphere 17, 24. The environmental connection to TC climate implies that the continuous warming of the environment is likely to lead to record-breaking TC intensity during El Niño as well as a falling frequency during La Niña

The geometric model depicts how the overall TC-climate relationship changes over time. It is surprising that the TC-climate connection after adding only six more years of observations (2015–2020) is remarkably weakened in explanatory power (Fig. 2 red line). Annual TC activity (\(P1_{Y}\)) and intensity (Y1) retain some of the earlier connection, but much of the TC variability indicates weaker connections. The environmental direction of \(P1_{Y}\) indicating a rotation factor (\(\theta _{1}\)) shows comparatively small change. Maximum covariance direction (\(\theta _{2}\)), the angle of the peak correlation and \(P1_{Y}\), is shifted by the narrower correlation peak. Scaling factor (r) in \(\theta _{2}\) reduced from 0.74 to 0.70 by the larger interruptions. Reduction of the correlation coefficients in most environmental variability directions implies a steeper tilting (\(\theta _{3}\)). \(\theta _{3}\) from 46.5\(^{\circ }\) to 73.1\(^{\circ }\) proves that the recent climatic events were unevenly distributed over the variability plane compared to the former events. The largest change occurs in the variability direction around the response to global ocean warmth (X1) implying that recent additional warming produces a response that is now different (Supplementary Fig. S9). Did the past observations contribute to the appearance of a stable climate connection by chance or are the recent observations evidence of an emerging imbalance in the connection among climate conditions? Since the current correlation approach is only an average of the annual covariance elements, the correlation itself cannot answer how the climatic relationship changes over time.

Sudden and still ongoing drift

We look into the annual covariance elements to see how the noticeable change in the geometric variability model occurred, which is considered a novel approach to validating a statistical model. Statistically, the correlation coefficient (r) between the two standardized variables of X and Y is calculated by \(\sum _{t=1}^{N}X_{t} \cdot Y_{t}/(N-1)\), where N is the number of years. This study defines \(X_{t} \cdot Y_{t}\) as the annual covariance element. Let’s think of an annual covariance field at certain paired variability directions (Supplementary Fig. S12). We focus on the functionality of the maximum covariance element. Because any directional covariance element (\(X_{t} \cdot Y_{t}\)) in year (t) can be functionally determined by the maximum TC variability (\(mX_{t}\)) and the maximum environmental variability (\(mY_{t}\)), the direction and the magnitude of annual maximum covariance element (\(mX_{t} \cdot mY_{t}\)) are the most important determinants of the TC-climate connection in a variability model. The prefixed “m” denotes annual “maximum value”.

Firstly, we examine the annual variation of the rotation factor (\(\theta _{1}\)). The position of \(mX_{t} \cdot mY_{t}\) identifies the environmental forcing and the TC response at the same time. The angle between the variability directions of \(mX_{t}\) and \(mY_{t}\) implies \(\theta _{1}\) which shapes overall climate relationship in a certain year (t). Fig. 3 A shows annual variation of \(\theta _{1}\) indicating how TC variability plane is positioned on the environmental variability plane. Since the variability direction is cyclic, the distribution is tiled as a mosaic. One of the key findings is that random annual \(mX_{t}\) and \(mY_{t}\) events are in reality constrained by a climate relationship defined by \(\theta _{1}\). \(\theta _{1}\) in Fig.

Secondly, the annual variation of the scaling factor (r) in the maximum covariance direction (\(\theta _{2}\)) is examined. Here, the annual maximum covariance element (\(mX_{t} \cdot mY_{t}\)) is usefully exploited again. Though r cannot be estimated by a single \(mX_{t} \cdot mY_{t}\), its difference from each given environmental condition (\(mX_{t} \cdot mY_{t}-mX_{t} \cdot mX_{t}\)) clearly shows how the environmental connection is interrupted by unknown variables (Fig. 3 C). Since the horizontal dashed line at zero means no interruption, the farther values from the line imply the larger interruptions and the greater reductions in the explanatory power of the environment. Distinct unprecedented interruptions occurred in 2014 and 2015 and still fluctuate. A notable point is that, as we showed above, during 2014 and 2015 it looked as if \(\theta _{1}\) temporarily returned to the former status during the drift but in reality, the TC-climate connection during these two years suffered from huge interruptions.

Lastly, the tilting factor (\(\theta _{3}\)) is examined by the group of recent observations. The smaller \(\theta _{3}\) means the longer projection length of the TC variability onto the environmental variability plane, which occurs when the maximum covariance elements are evenly distributed over the environmental variability directions. While the points over the earlier period 1985–2012 are evenly distributed on the plot, the points over the later period are irregularly clustered with gaps indicating a large difference between the maximum and minimum covariance (see Fig. 3 B) confirming a contribution to the larger \(\theta _{3}\) as shown by a loss of explanatory power (Fig. 2).

The regional aspect of global warming could be looked into by the regional responses to GMSST. It has been known that the global warming indicated by GMSST significantly increases the efficiency of TC intensification by the dominant high pressure anomalies with the more unstable atmosphere in the western North Pacific 24. Even with the drifting environmental connection to TC climatology, we find that this region has experienced the stronger environmental responses to global warming in recent years (Supplementary Fig. S11). The correlation implies that the tropical region where most TC activity occurs has become even warmer. High pressure anomaly is also shown much stronger over the same region. On the other hand, Fig. 4 shows the time series of geopotential height at 500 hPa compared to the efficiency of TC intensity which is considered to regulate TC intensity and frequency at a given level of TC activity. All values are annually averaged over JJASON in the tropical region of the western North Pacific (0\(^\circ\)–30\(^\circ\) N, 100\(^\circ\) E–180\(^\circ\)). The correlation during 28 years (1985–2012) drops dramatically from 0.72 {[0.48, 0.86], 95 % confidence interval (CI)} to 0.42 ([0.11, 0.66], 95 % CI) when the 36-year period (1985–2020) is applied, which implies that the recent TC intensity and frequency now cannot be effectively explained by atmospheric suppression due to global warming. The recent split of the two variables in the time series also proves that the weaker correlation does not simply mean some smaller magnitude of the environmental connection but a mixture of unexpected observations in recent years.

This analysis shows conclusively that the former environmental relationship to TC climate are far less valid, and they are drifting to a new state. The analysis is done on western North Pacific super typhoons since collectively they are the most energetic but similar changes to the TC-climate environment relationships are found for North Atlantic and eastern North Pacific hurricanes and their environments (Supplementary Fig. S13) and collectively for all basins (Supplementary Fig. S14) indicative of a global imbalance in TC climatology. This study deals with TC climatology which is only related to intensity and frequency. For a better understanding of TC climatology, the scope of future research may include TC track climatology.

Conclusions and discussion

This study investigates the balance of the climate system by decomposing the environmental connection to super typhoon activity which is one of the central issues regarding climate change and natural disasters. A model for the environmental connection is refined using geometric factors of the TC variability plane (rotation, maximum covariance direction, scaling, tilting), which depicts the variability directions and the explanatory powers in a three-dimension variability space. The geometric variability model effectively shows the functional connection between TC climate and the environment. The framework provides a bird’s-eye view of where the individual variables are located and how they are experiencing changes on interannual timescales. It is noted that only adding several recent years leads to remarkable reduction in the explanatory power of the environmental forcings on the super typhoon climatology. This study first time shows how a deformation of the former climatic connection among variables, in a changing climate, is printed on annual covariance elements. We find that the recent observations showing a group of outlying events with a certain drift are more than unfamiliar compared to the stable connection from 1985 through 2012.

The future research could focus on how it happens. The functional imbalance may be caused by a feedback process involving TCs themselves or by changes to the larger climate conditions. On the other hand, the sudden and persistent drift suggests to us that the super typhoon climatology might be entering a “tipping point” over a threshold 38. The environmental conditions over a tipping point may turn out to be the emergence of new environmental relationships in earth’s climate system. For example, changes to the environmental conditions induced by sea-ice loss and accelerated glacier’s melting due to global warming 39 could be disturbing the former TC-climate connection. Increasing freshwater by the melting of Greenland’s ice sheet can weaken the thermohaline circulation 40, and the weakened circulation would influence the SST pattern in the North Atlantic 41. Subsequently, the modified North Atlantic SST pattern may regulate the storm activities in the tropical Pacific 42, 43, which describes a connection of the Atlantic meridional overturning circulation as a potential tipping element 44 for super typhoon climatology. A drifting climate will affect the performance of numerical approaches as well. Despite the reliability of numerical models in simulating the future TC climate with higher resolutions 45, 46, performance is limited to the warming environment 47 as well as to ENSO variability 48. To the extent that the historical environmental and TC observations are used to correct and verify the models, a destabilized climate with a new TC-climate relationship to the environment will make these approaches less reliable in their predictive performance. Consequently, this study suggests that a climate-change impact could take the form of uncertain diagnostics and prognostics about TC activity and thereby magnifying concerns about how to cope with the climate crisis.

Methods

First the response variability plane and the explanatory variability plane is constructed from the TC variables and the environmental variables, respectively. Then, the geometric association of the two variability planes is modeled in the three-dimensional space. The geometric factors are (1) rotation (\(\theta _{1}\)), (2) scaling (r) in the maximum covariance direction (\(\theta _{2}\)), and (3) tilting (\(\theta _{3}\)) of the response variability plane (see Fig. 3). The annual contribution to geometric factors and the TC-climate connection to the environment is examined by the annual maximum covariance elements.

Data and code availability

The TC variables are calculated from the best-track data provided by the US Joint typhoon warning center ( https://www.metoc.navy.mil/jtwc/jtwc.html?best-tracks). GMATS, SOI, PMM, and PDO from the National Oceanic and Atmospheric Administration (NOAA)/Climate Prediction Center are the environmental variables. The global mean sea surface temperature (SST) is calculated from the Extended Reconstructed Sea Surface Temperature version 5 (ERSSTv5) of the NOAA/National Centers for Environmental Prediction (NCEP) reanalysis 49. The same for the Niño indices such as Niño 1+2, Niño 3, Niño 3\(\cdot\)4, and Niño 4 for the following regions: (1) 10° S–0\(^\circ\) and 90\(^\circ\) W–80\(^\circ\)W for Niño 1+2, (2) 5\(^\circ\) S–5\(^\circ\) N and 150\(^\circ\) W–90\(^\circ\) W for Niño 3, (3) 5\(^\circ\) S–5\(^\circ\) N and 170\(^\circ\) W-120\(^\circ\) W for Niño 3\(\cdot\)4, and (4) 5\(^\circ\) S–5\(^\circ\) N and 160\(^\circ\) E–50\(^\circ\) W for for Niño 4. All SST-based variables calculated from ERSSTv5 are area-weighted values. All of the statistics and figures are made using the software R ( https://www.r-project.org) and code is available online ( https://rpubs.com/Namyoung/P2022a).
Alimonti, G., Mariani, L., Prodi, F. & Ricci, R. A. A critical assessment of extreme events trends in times of global warming. Eur. Phys. J. Plus 137(1), 1–20 (2022).
Von Der Heydt, A. S. et al. Quantification and interpretation of the climate variability record. Global Planet. Change 197, 103399.
https://doi.org/10.1016/j.gloplacha.2020.103399(2021).
IPCC, 2018: Summary for Policymakers. In: Global Warming of \(1.5^{\circ }\text{C}\). An IPCC Special Report on the impacts of global warming of \(1.5^{\circ }\text{ C }\) above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty. Cambridge University Press, Cambridge, UK and New York, NY, USA, pp. 3–24.
https://doi.org/10.1017/9781009157940.001(2018).
Ripple, W. J., Wolf, C., Newsome, T. M., Barnard, P. & Moomaw, W. R. Corrigendum: World scientists’ warning of a climate emergency. Bioscience 70(1), 100 (2020).
Ripple, W. J. et al. World scientists’ warning of a climate emergency 2021. Bioscience 71(9), 894–898 (2021).
Cai, Y., Lenton, T. M. & Lontzek, T. S. Risk of multiple interacting tipping points should encourage rapid Co\(_2\) emission reduction. Nat. Clim. Change 6(5), 520–525 (2016).
Lenton, T. M. et al. Climate tipping points-too risky to bet against. Nature 575, 592–595 (2019).
Boers, N., & Rypdal, M. Critical slowing down suggests that the western Greenland ice sheet is close to a tipping point. Proc. Natl. Acad. Sci. 118(21), e2024192118 (2021).
Mendelsohn, R., Emanuel, K., Chonabayashi, S. & Bakkensen, L. The impact of climate change on global tropical cyclone damage. Nat. Clim. Change 2, 205–209 (2012).
Koks, E. E. et al. A global multi-hazard risk analysis of road and railway infrastructure assets. Nat. Commun. 10(1), 1–11 (2019).
Pinnegar, J. K., Engelhard, G. H., Norris, N. J., Theophille, D. & Sebastien, R. D. Assessing vulnerability and adaptive capacity of the fisheries sector in dominica: Long-term climate change and catastrophic hurricanes. ICES J. Mar. Sci. 76(5), 1353–1367 (2019).
Trenberth, K. E., Caron, J. M. & Stepaniak, D. P. The atmospheric energy budget and implications for surface fluxes and ocean heat transports. Clim. Dyn. 17, 259–276 (2001).
Emanuel, K. Contribution of tropical cyclones to meridional heat transport by the oceans. J. Geophys. Res. Atmos. 106(D14), 14771–14781 (2001).
Trenberth, K. E. & Stepaniak, D. P. The flow of energy through the earth’s climate system. Quart. J. Roy. Meteor. Soc. 130(603), 2677–2701 (2004).
Trenberth, K. E., & Fasullo, J. Water and energy budgets of hurricanes and implications for climate change. J. Geophys. Res. Atmos. 112, D23107 (2007).
Kang, N., Kim, D. & Elsner, J. B. The contribution of super typhoons to tropical cyclone activity in response to ENSO. Sci. Rep. 9, 5046.
https://doi.org/10.1038/s41598-019-41561-y(2019).
Kang, N. & Elsner, J. B. Trade-off between intensity and frequency of global tropical cyclones. Nat. Clim. Change 5, 661–664 (2015).
Chu, J. H., Sampson, C. R., Levine, A. S., & Fukada, E. The Joint Typhoon Warning Center tropical cyclone best-tracks, 1945–2000. Joint Typhoon Warning Center, 22 pp (2002).
Klotzbach, P. J. Trends in global tropical cyclone activity over the past twenty years (1986–2005). Geophys. Res. Lett. 33, 10805 (2006).
Kang, N. & Elsner, J. B. Consensus on climate trends in western North Pacific tropical cyclones. J. Clim. 25, 7564–7573 (2012).
Elsner, J. B., Kossin, J. P. & Jagger, T. H. The increasing intensity of the strongest tropical cyclones. Nature 455, 92–95.
https://doi.org/10.1038/nature07234(2008).
Bell, G. D. et al. Climate assessment for 1999. Bull. Am. Meteor. Soc. 81(6), 1–50.
https://doi.org/10.1175/1520-0477(2000)81[s1:CAF]2.0.CO;2(2000).
Emanuel, K. A. Increasing destructiveness of tropical cyclones over the past 30 years. Nature 436, 686–688.
https://doi.org/10.1038/nature03906(2005).
Kang, N. & Elsner, J. B. Climate mechanism for stronger typhoons in a warmer world. J. Clim. 29, 1051–1057.
https://doi.org/10.1175/JCLI-D-15-0585.1(2016).
Lin, Y.-F., Yu, J.-Y., Wu, C.-R. & Zheng, F. The footprint of the 11-year solar cycle in Northeastern Pacific SSTs and its influence on the Central Pacific el niño. Geophys. Res. Lett. 48(5), 2020–091369 (2021).
Amaya, D. J. et al. The north pacific pacemaker effect on historical ENSO and its mechanisms. J. Clim. 32(22), 7643–7661 (2019).
Sanchez, S. C., Amaya, D. J., Miller, A. J., Xie, S.-P. & Charles, C. D. The pacific meridional mode over the last millennium. Clim. Dyn. 53(5), 3547–3560 (2019).
Kang, N. & Elsner, J. B. The changing validity of tropical cyclone warnings under global warming. NPJ Clim. Atmos. Sci. 1, 36.
https://doi.org/10.1038/s41612-018-0047-3(2018).
Stuecker, M. F. Revisiting the pacific meridional mode. Sci. Rep. 8(1), 1–9 (2018).
Peixoto, J. P. & Oort, A. H. Physics of Climate (American Institute of Physics, New York, 1992).
Elsner, J. B. Continued increases in the intensity of strong tropical cyclones. Bull. Am. Meteor. Soc. 101(8), 1301–1303 (2020).
Williams, I. N., Pierrehumbert, R. T., & Huber, M. Global warming, convective threshold and false thermostats. Geophys. Res. Lett. 36, L21805.
https://doi.org/10.1029/2009GL039849(2009)
Johnson, N. C. & Xie, S.-P. Changes in the sea surface temperature threshold for tropical convection. Nat. Geosci. 3(12), 842–845 (2010).
Sugi, M., Murakami, H. & Yoshimura, J. On the mechanism of tropical cyclone frequency changes due to global warming. J. Meteorol. Soc. Jpn. Ser. II 90, 397–408 (2012).
Sugi, M., Yoshida, K. & Murakami, H. More tropical cyclones in a cooler climate?. Geophys. Res. Lett. 42(16), 6780–6784 (2015).
Tang, B. & Emanuel, K. A ventilation index for tropical cyclones. Bull. Am. Meteor. Soc. 93(12), 1901–1912 (2012).
Chavas, D. R. A simple derivation of tropical cyclone ventilation theory and its application to capped surface entropy fluxes. J. Atmos. Sci. 74(9), 2989–2996 (2017).
Gladwell, M. The Tipping Point: How Little Things Can Make a Big Difference (Little Brown, New York, 2000).
Hugonnet, R. et al. Accelerated global glacier mass loss in the early twenty-first century. Nature 592(7856), 726–731 (2021).
Caesar, L., Rahmstorf, S., Robinson, A., Feulner, G. & Saba, V. Observed fingerprint of a weakening Atlantic ocean overturning circulation. Nature 556(7700), 191–196 (2018).
Goldenberg, S. B., Landsea, C. W., Mestas-Nuñez, A. M. & Gray, W. M. The recent increase in Atlantic hurricane activity: Causes and implications. Science 293(5529), 474–479 (2001).
Dima, M., Lohmann, G. & Rimbu, N. Possible north Atlantic origin for changes in ENSO properties during the 1970s. Clim. Dyn. 44(3), 925–935 (2015).
Park, J.-H., Li, T., Yeh, S.-W. & Kim, H. Effect of recent Atlantic warming in strengthening Atlantic-Pacific teleconnection on interannual timescale via enhanced connection with the Pacific meridional mode. Clim. Dyn. 53(1), 371–387 (2019).
Lenton, T. M. et al. Tipping elements in the earth’s climate system. Proc. Natl. Acad. Sci. 105(6), 1786–1793 (2008).
Satoh, M. et al. Global cloud-resolving models. Curr. Clim. Change Rep. 5(3), 172–184 (2019).
Stevens, B. et al. Dyamond: The dynamics of the atmospheric general circulation modeled on non-hydrostatic domains. Prog. Earth Planet. Sci. 6(1), 1–17 (2019).
Jerez, S. et al. Impact of evolving greenhouse gas forcing on the warming signal in regional climate model experiments. Nat. Commun. 9(1), 1–7 (2018).
Bellenger, H., Guilyardi, É., Leloup, J., Lengaigne, M. & Vialard, J. Enso representation in climate models: From cmip3 to cmip5. Clim. Dyn. 42(7), 1999–2018 (2014).
Huang, B. et al. Extended reconstructed sea surface temperature, version 5 (ersstv5): upgrades, validations, and intercomparisons. J. Clim. 30(20), 8179–8205 (2017).

Acknowledgements

This work was supported by the Korea Meteorological Administration (KMI2021-01511), and the Korea Institute of Ocean Science and Technology (PEA0112).

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

About this article

Cite this article

Kang, N., Jang, C.J. & Elsner, J.B. More than unfamiliar environmental connection to super typhoon climatology. Sci Rep 13, 6372 (2023). https://doi.org/10.1038/s41598-023-33104-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-33104-3

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.