What caused the increase of tropical cyclones in the western North Pacific during the period of 2011–2020?

Based on satellite data after 1979, we find that the tropical cyclone (TC) variations in the Western North Pacific (WNP) can be divided into three-periods: a high-frequency period from 1979 to 1997 (P1), a low-frequency period from 1998 to 2010 (P2), and a high-frequency period from 2011 to 2020 (P3). Previous studies have focused on WNP TC activity during P1 and P2. Here we use observational data to study the WNP TC variation and its possible mechanisms during P3. Compared with P2, more TCs during P3 are due to the large-scale atmospheric favorable conditions of vertical velocity, relative vorticity and relative humidity. Warm sea surface temperature (SST) anomalies are found during P3 and migrate from east to west, which is also favorable for TC genesis. The correlation between the WNP TC frequency and SST shows a significant positive correlation around the equator and a significant negative correlation around 36°N, which is similar to the warm phase of the Pacific Decadal Oscillation (PDO). The correlation coefficient between the PDO and TC frequency is 0.71, above the 95% confidence level. The results indicate that the increase of the WNP TC frequency during 2011–2020 is associated with the PDO and warm SST anomalies.


Introduction
Tropical cyclone (TC) is one of the most intense weather systems in the world. A TC is called a typhoon when its wind speed exceeds 34 knots in the Western North Pacific (WNP). The WNP basin (0º-30ºN, 100ºE-180º) exhibits a high level of TC activity and is the most active region globally. About 26 named TCs occurred over the WNP basin each year on a climatological average, leading to huge loss of life and property (Zhang et al., 2009, Wang et al., 2021. Many studies point out an increasing pattern of TC-induced damages over Dr. Chunzai Wang cwang@scsio.ac.cn 1 atmospheric changes and remote oceanic SST changes. The rest of this paper is organized as follows. Section 2 describes the data and method used in this study. The TC increase is illustrated in Sect. 3. In Sect. 4, we explore the possible influences of large-scale atmospheric environment and oceanic conditions associated with TC increase. The discussion and conclusions are presented in Sect. 5.

TC data
The TC best-track data, including the TC intensity, latitude and longitude, are from the United States Joint Typhoon Warning Center (JTWC). The selected period for the present study is the period after the beginning of the satellite era from 1979 to 2020 (https://www.ncdc.noaa.gov/ibtracs/ index.php?name=ib-v4-access). In this study, we focus on the TC frequency throughout the year (from January to December). The TC frequency in this study is calculated over the whole WNP region. The WNP domain (0°-30°N, 100°E-180°) is divided into 5° × 5° resolution grid boxes. TC genesis positions are counted for each grid box. The first position of a TC recorded in JTWC is defined as the TC genesis location. Student's t test is performed to test statistical significance.

Atmospheric and SST data
The monthly mean atmospheric fields (e.g., winds) are obtained from the United States National Centers for Environmental Prediction and National Center for Atmospheric Research (NCEP/NCAR) II monthly reanalysis dataset with a 2.5° × 2.5° grid (Kalnay et al. 1996) (https://psl.noaa.gov/ data/gridded/). The monthly mean SST is taken from the United States National Oceanic and Atmospheric Administration (NOAA) Extended Reconstruction SST version 5 (ERSSTv5) with a resolution of 2° × 2° (Huang et al. 2017).

Selection of ENSO events
The Niño 3.4 SST index is calculated from the HadISST1 dataset (Rayner et al. 2003) (https://psl.noaa.gov/gcos_ wgsp/Timeseries/Nino34/). The Oceanic Niño Index (ONI) from the Climate Prediction Center is used to distinguish El Niño and La Niña years based on a 0.5 °C threshold. (https://origin.cpc.ncep.noaa.gov/products/analysis_monitoring). The years with SST anomalies more positive (more negative) than 0.5 °C (− 0.5 °C) during October to December are defined as El Niño (La Niña) years. The other years are considered as neutral years. Detailed classifications are Mei and Xie 2016). The Madden-Julian Oscillation (MJO) with a 30-90 day period can regulate WNP TC activity (Li and Zhou, 2013;Liu et al., 2021;Wang et al., 2018;Zhao and Wang, 2019). In the WNP basin, the active phase of the MJO enhances the local convection, which in turn leads to the increase in TC activity in the basin. The quasi-biweekly oscillation is another intra-seasonal mode, but with a more localized effect, that also affects WNP TC activity (Li and Zhou 2013). The El Niño-Southern Oscillation (ENSO) is known as the main phenomenon occurring at interannual time scales (Trenberth and Caron 2000;Timmermann et al. 2018;Kim et al. 2020). ENSO has significant impacts on WNP TCs (Kim et al., 2011). Wang (2016, 2019) examine a stronger interannual relationship between ENSO and TCs in the WNP during the boreal summer from 1998 to 2015. The correlation is 0.60 at the 95% confidence level. At decadal time scales, many factors affecting the TC frequency have already been suggested. Researchers mainly focus on the following factors: (a) the strengthened western North Pacific subtropical high (Liu and Chan, 2020); (b) the anomalous westward WNP SST gradient (Choi et al., 2015); (c) phase shifts of the Pacific Decadal oscillation (PDO) (Zhao and Wang, 2016); (d) the interdecadal Pacific oscillation (He et al., 2015;Zhao et al., 2018), and also (e) the Atlantic Multi-decadal Oscillation (AMO) (Zhang et al., 2020;Zhang et al., 2017).
Most studies suggest that the frequency of WNP TCs decreased from 1979 to 2010 (Tu et al., 2009;Yokoi and Takayabu, 2013;Hsu et al., 2014). Few studies have examined the change of TC activity from 2010 to 2020. Most of the previous analyses on TCs focus on the first or the second half of the year, or a certain period of the year, and the number of TCs throughout the year is less covered in previous studies. With global warming and the emergence of various extreme weather phenomena, how did the TC frequency per year change during 2011-2020 compared to the pre-2010 decades? What caused those changes? To our knowledge, these issues have not been explained in previous studies. Here we investigate the possible reasons associated with the changes in the WNP TC frequency during 2011-2020 based on observations and reanalysis data.
The aim of this paper is to explore the change of the WNP TC frequency, and the corresponding local large-scale Table 1 The list of El Niño years, La Niña years and neutral years during P2 (1998-2010) and P3 (2011-2020) P2 (

Effective degrees of freedom
In order to consider the reduction in the degrees of freedom of the filtered time series, the effective degrees of freedom, which was used by Chen et al. (2013) and Li and Zhou (2014), are used to evaluate the significance between the filtered variables. The effective degrees of freedom are obtained as n/T, where n is the sample size and T = K t=0 c xx (t) c yy (t). c xx and c yy are the autocorrelation coefficients of the two variables with a time lag of t, the maximum of integer K corresponds to n/2.

TC frequency changes
The present work focuses on the total number of WNP TCs each year (from January to December), and all analyses use the annual mean data. The time series of the WNP TC number is displayed in Fig. 1. Based on the variations, we divide 1979-2020 into three epochs, with turning points in 1998 and 2010, respectively. The number of TCs in the WNP during these two years of turning points is the least. Moreover, 1998 is approximately the PDO phase transition year from the positive to negative (warm to cold) phase (Zhao and Wang, 2016). Compared with P2, there is a significant and abrupt increase in the WNP TC frequency during P3 (Fig. 1). Since the TC frequency between P1 and P2 has been previously studied (e.g., Wang, 2016, 2019;Choi and Kim, 2019;Zhao et al., 2021), this study explores shown in Table 1. P1, P2 and P3 represent the periods of 1979-1997, 1998-2010 and 2011-2020, respectively.

Genesis potential index
To determine the factors that influence changes in TC genesis, we used a modified version of genesis potential index (GPI) (Murakami et al. 2011): where η is the absolute vorticity (s − 1 ) at 850 hPa, RH is the relative humidity (%) at 600 hPa, V s is the magnitude of the vertical wind shear (ms − 1 ) between 850 and 200 hPa, and ω is the vertical pressure wind velocity (Pa s − 1 ) at 500 hPa, V pot is the maximum potential intensity (PI) (ms − 1 ): source code for V pot is available online at http://texmex.mit.edu/pub/ emanuel/TCMAX/pcmin.m), C k is the exchange coefficient for enthalpy, T s is SST (K), C D is the drag coefficient, T 0 is the mean outflow temperature (K), CAP E * is the value of convective available potential energy of the air lifted from saturation at sea level, and the variable CAP E b is the boundary layer air. To assess the individual importance of these five variables, we recalculate the GPI using the longterm climatology of four variables but varying the annual values for one variable. We carried out this procedure for each of the five variables, followed by recalculating GPI differences between P2 and P3. In neutral years, the average TC counts in P2 (P3) are about 23.5 (28.3). In El Niño years, 25.25 (26.33) TCs are observed during P2 (P3). The TC frequency difference in El Niño years, La Niña years and neutral years is significant lower than the 95% confidence level. Figure 2 also shows that there is more TC generation in La Niña years and neutral years during P3 compared to P2. Furthermore, most TCs form in neutral years than in La Niña and El Niño the reasons for the WNP TC frequency rise during P3 compared to P2. Figure 2 shows that more TCs can be found during P3 than during P2. The difference between P2 and P3 is 3.3, which is significant at the 95% confidence level. It is shown that the major difference between P2 and P3 is due to the differences during neutral and La Niña years. In La Niña years, the average TC number in P2 (P3) is about 20 (24).

Fig. 3
Differences in WNP TC genesis distribution (×10) between P2 and P3 (P3 minus P2) in each 5° x 5° grid box. The dashed rectangle indicates the region used for TC numbers difference amalgamating in the west (6°-14°N,100°-130°E) and east (3°-14°N,140°E-180°) regions. The slash in the west box indicates significance at the 95% confidence level. The significance in the east box is lower than the 95% confidence level  (1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010) and P3 (2011-2020). The TC frequency differences between P2 and P3 are statistically significant at the 95% confidence level for all years, but not for La Niña years, neutral years, and El Niño years. During P2 and P3, the frequency difference between El Niño years and La Niña years is statistically lower than the 95% confidence level during P3. The genesis of TCs is generally influenced by favorable environmental conditions, such as the increase of the relative vorticity and the decrease of the vertical wind shear, high SST and strengthened vertical velocity (omega). The increase in the relative vorticity at the 850 hPa level can concentrate more convective heating, which is beneficial to the TC development. The vertical wind shear is usually calculated as the difference between wind vectors at 200 hPa and 850 hPa. The decreased vertical wind shear can reduce the advection of heat and moisture from the TC center. Figure 4a shows that GPI differences between P2 and P3 are broadly consistent with the corresponding TC genesis differences over the South China Sea (Fig. 3). There are significant positive areas over the South China Sea. However, some differences also exist in Figs. 3 and 4a, particularly years in P3, while most TCs are in El Niño years during P2. It is shown that there is no significant change of the TC frequency between P2 and P3 in El Niño years. During P2 and P3, the frequency difference between El Niño years and La Niña years is statistically significant lower than the 95% confidence level. Figure 3 shows that TCs in the northwest, southeast, and central areas increase significantly during P3 compared to P2. There is a slight reduction in the areas 6°-12°N and 137°-143°E, 12°-18°N and 122°-128°E. The TCs move while they develop, and stochastic factors will also be an issue for TCs. Physical forcing in one region may lead to TC changes in a different region. So here we try to amalgamate TC numbers and average quantities over large regions. The amalgamating TC number difference in the west (6°-14°N,100°-130°E) and east (3°-14°N,140°E-180°) regions is shown in Fig. 3. The TC frequency in the west box increases (1.63) significantly (at the 95% confidence level). However, the TC frequency change in the east box (1.86) is lower than the 95% significant confidence level. TC frequency change is more significant in the western WNP than in the eastern part. More TCs can be observed over the South China Sea while relative vorticity and vertical pressure velocity contribute positively to the GPI. The other three terms occupy a large proportion of the equation for calculating GPI, which in the eastern box. According to the contribution of the five variables to the GPI in the eastern region, wind shear, relative humidity and PI contribute negatively to the GPI,  Fig. 3. GPI denotes the GPI change between P2 and P3; Shear denotes the GPI change caused by the vertical wind shear; Vor denotes the GPI change caused by absolute vorticity at 850 hPa; Hum denotes the GPI change caused by the relative humidity at 600 hPa; and W denotes the GPI change caused by the vertical pressure velocity at 500 hPa vorticity is positive and has an important influence on the TC genesis.

Differences in GPI between P2 and P3
To investigate the differences in environmental variables among El Niño and La Niña events, Fig. 6 displays the anomalies of five environmental variables constituting the GPI, including PI, Hum, Vor, W, and Shear. Figure 6a-e shows the contribution of each term of Eq. (1) in individual La Niña year during P2. In contrast to P2, the GPI in La Niña years of P3 is all positive (Fig. 6f-i), which is favorable for TC genesis. It is consistent with the result shown (Fig. 2) that the TC frequency in La Niña year of P3 is more than that in P2. Figure 6f-i also demonstrates that the vertical pressure velocity at 500 hPa positively contributes to the TC genesis. However, Fig. 6j-p shows that different environmental factors that are beneficial to TC activities in different El Niño year, as previous studies have shown that El Niño does not significantly impact WNP TC number (Li and Zhou, 2013). Figure 6q shows that the TC number generated in La Niña year is less than the number generated in El Niño year in the domain 5°-15°N, 105°E-180°. In La Niña years, the TC number of P3 is greater than that of P2; and in El Niño years, the TC number between P2 and P3 is slightly different.
It is worth noting that the GPI (Fig. 6a-p) has some differences with the TC number. One reason may be that the TC genesis position is not in the selected domain. There are TCs generated outside the selected domain in El Niño or La Niña years. The selected domain has little influence on the generation of these TCs, but has a very important influence on the subsequent development of TCs.Another reason maybe the change of environmental variable that affects weakens the influence of vertical pressure velocity and relative vorticity. Only small areas of the WNP have significant differences in vertical wind shear (Fig. 4b). Figure 4c shows that the variable PI does not contribute to the TC increase in P3, while the relative vorticity greatly contributes to the increase of TCs frequency in P3 (Fig. 4d). Although the vorticity is not significant in the east of 150°E, the vorticity appears to impact TC formation (Fig. 3). Figure 4e and f indicate that the relative humidity (west box) and vertical pressure velocity are favorable for the TC frequency increase in P3, especially the vertical pressure velocity. These findings suggest that the vorticity, relative humidity and vertical pressure velocity cause more TCs in the western box in P3. Additionally, larger vorticity and vertical pressure velocity play an important role in increasing the TC frequency in the eastern box. Figure 5 shows the contributions of each term of Eq.
(1) to TC genesis. Figure 5a indicates favorable conditions for TC frequency increase in the domain of 5°-15°N, 105°E-180°. The larger vorticity, relative humidity and vertical pressure velocity influence the TC genesis in the domain. The whole GPI is positive, which is in favor of TC genesis. Figure 5b displays the contribution of the five variables to the TC genesis in the west box of Fig. 3. The whole GPI is positive, and wind shear, relative vorticity, relative humidity and vertical pressure velocity contribute positively to TC genesis, except PI. Figure 5c indicates that the whole GPI is negative, the relative vorticity is positive and contribute to the TC genesis. According to the contribution of the five variables to the TC genesis shown in Fig. 5, the relative Fig. 7 Differences of (a) relative humidity (%) at 600 hPa, (b) omega (Pa s − 1 ) at 500 hPa, (c) relative vorticity (10 − 6 s − 1 ) at 850 hPa, and (d) vertical wind shear (m s − 1 ) between P2 and P3 (P3 minus P2). The black crosses indicate that the statistical difference in this region is significant at the 95% confidence level Figure 7 shows that there are significant positive relative humidity differences in the western box, and a significant enhancement in vertical velocity (omega) in the southeastern part (123°E-180°, 0°-13°N) of the WNP between P2 and P3. An increase in relative vorticity can also be found in 123°E-180°, 0°-18°N area. The increase in TC frequency in the east of the WNP (east of 140˚E) is mainly influenced by TC activity during El Niño events is complex. Song et al. (2020b) examined the impacts of the three types of El Niño events on TCs over the WNP. They pointed out that different types of El Niño correspond to different environmental factors that are beneficial to TC activities. The TC frequency over western part of the WNP increases when the SST anomalies migrate from east to west. Zhang et al. (2015) reported that the occurrence of TCs in the WNP basin is different during the different kinds of El Niño years.   Fig. 7, but for La Niña years that conditions in the above areas are conducive to more TC precursor disturbances that can spread to the west or north to develop into TCs in P3 than P2. Corresponding to the genesis position of TCs, we can find that the occurrence of relative vorticity and vertical velocity. Figure 7d also shows that the vertical wind shear in the WNP during P3 is smaller than that during P2. However, only the 165°E-180°,0°-12°N area is significant at the 95% confidence level. We speculate  respectively. Although these correlations are not significant, it is interesting to see the sign change of the correlations. In P1, the correlation is the largest and negative. The larger number of TCs occurs during La Niña events. But during the phase transition of the PDO around 1997, the increase in TC corresponds to the increase in Niño index, and the larger number of TCs appears during El Niño events.

Differences in atmospheric conditions between P2 and P3
To emphasize the decadal variations, we perform the 7-year running mean on the data and repeat the correlation calculations (Fig. 11c). The patterns of correlation distribution in Fig. 11a are much more similar to Fig. 11c than that to Fig. 11b, indicating that the correlation pattern in Fig. 11a (at both interannual and decadal time scales) is dominated by decadal variability. A significant negative correlation around 36°N can also be seen in Fig. 11a. The correlation patterns in the Pacific are similar to those of the PDO. Therefore, we compare the TC frequency variation at decadal time scales with the PDO index (Fig. 11e). The correlation between the PDO and WNP TC frequency during 1982-2017 is 0.71 which is above the 95% confidence level. Figure 11e also shows that the decadal relationships of WNP TCs with the PDO are opposite before and after the middle 1990s. The correlations during 1982-1995 and 1996-2017 are − 0.84 and 0.90, respectively. The TC frequency is negatively (positively) correlated with PDO index before (after) 1995. The TC frequency in the positive phase of the PDO is larger than that in the negative phase. When the PDO is in the positive phase, the corresponding SST can affect the large-scale environmental factors that are more favorable for TC development.These suggest that in the Pacific, the PDO has a stronger relationship with WNP TCs than ENSO.

Compared with P1
The TC frequency during the late season (October-December) or from June to August in P1 has previously been studied Wang, 2016, 2019). In this section, the changes in the factors influencing the TC frequency between P1 and P3 are discussed. Wang (2016, 2019) previously point out that the warm phase (P1) is related to the increase in the TC frequency and the decrease of interannual correlation between ENSO and TCs, while the cold phase (P2) is related to the decrease in the TC frequency and the increase in interannual correlation between ENSO and TCs. The changes of relative vorticity and vertical wind shear in the WNP play a major role in the change of the TC frequency during P1 compared to P2. The change in SST in the tropical Indo-Pacific is closely associated with the correlation between ENSO and TCs. According to the analysis of the P3 period, we find that the WNP TC frequency also increases, which are consistent with the result of P1 (Zhao and Wang, 2016, TCs in the western box is mainly affected by relative humidity. The increase in the TC frequency in the eastern box is mainly due to the relative vorticity and vertical velocity. The increase in the TC frequency is mainly related to the relative humidity, the enhanced vertical velocity and the low-level relative vorticity, which is also consistent with the histogram of GPI (Fig. 5). Figures 8 and 9 display differences in atmospheric conditions between P2 and P3 in La Niña years and El Niño years, respectively. Figure 8 shows that in La Niña years, the vertical velocity (130°E-180°,0°-12°N), the relative vorticity (120°E-180°,0°-18°N), the relative humidity (100°-135°E, 0°-3°N) and the vertical wind shear (135°E-180°,0°-12°N) are conducive to the development of TCs. The large-scale climate differences for El Niño years are shown in Fig. 9. Compared with P2, the favorable factors for the development of TCs in P3 are not significant, although vertical velocity (145°E-180°,2.5°-10°N) and relative vorticity (145°-153°E,12°-15°N; 105°-120°E,0°-8°N) play a role. This also corresponds to the fact that the frequency changes of TCs in El Niño years in Fig. 2 are not significant.

Relations to SST variations
As suggested by previous studies (Zhan et al., 2011;Hsu et al., 2014;Yamaguchi et al., 2020), the different SST anomaly distributions can impact atmospheric circulation in the WNP. The SST changes in the WNP region are warm between P2 and P3 (Fig. 10), the warm SST anomalies migrate from east to west. We speculate that the higher SST can influence large-scale environmental factors related to TCs. We calculate the correlation of the WNP TC frequency throughout the year (January-December) with SST over the three ocean basins during 1979-2020 (Fig. 11a). A significant negative correlation appears in the northern Indian Ocean, due to the remote SST warming effect on WNP TCs (Wang et al. 2021). The significant negative correlation in the Atlantic Ocean is related to the AMO (Zhang et al. 2018). Figure 11b shows the interannual relationship between the TCs frequency and SST. The interannual variable is calculated by subtracting the decadal variable which is represented by a 7-year running mean, from the total time series. It shows that a negative relationship appears in the equatorial central and eastern Pacific. The distribution is most likely the ENSO pattern. Figure 11d compares the TC frequency throughout the year (January-December) with the Niño 3.4 index. The correlation during 1982-2017 is -0.16, indicating that ENSO is not significantly correlated with the TC number in the whole WNP (e.g., Wang and Chan 2002). The correlations between the Niño 3.4 index and TC frequency during P1, P2 and P3 are − 0.35, 0.23 and − 0.11, years are the lowest. The reason may be that the larger value of the TC frequency appears in the negative phase of ENSO in P1, while the larger value of the TC frequency appears in the positive phase of ENSO in P3 (Figs. 6q and 11d). It also may be due to the warm SST during the P3 period compared with P2 period (Fig. 10). Li and Zhou (2014) point out that the zonal SST gradient between the cooling WNP SST and the warming northern Indian Ocean SST is not favorable for TC genesis. The warming WNP SST during P3 reduces the zonal SST gradient, resulting in a favorable environment for TC genesis (Liu and Chan, 2020).
Our research also suggests that the changes of dynamical factors like low-level relative vorticity, relative humidity and vertical velocity play an essential role in TC genesis during P3. Vertical velocity, relative vorticity and relative humidity play a major role in La Niña years. For El Niño years, compared with P2, the favorable factors for the development of TCs in P3 are not significant. This also corresponds to the fact that the frequency changes of TCs in El Niño years in Fig. 2 are not significant. The WNP SST has an important influence on TC activity (Hsu et al. 2014;Song et al. 2020). The SST changes in the WNP are characterized by positive SST anomalies in P3 that migrate westward in the North Pacific. We study the correlation between the WNP TC frequency and SST. The result shows that the WNP TC frequency and SST have a significant positive correlation region and a significant negative correlation region. The correlation distributions in the WNP are related to the PDO pattern at decadal time scales (Fig. 11). To further confirm the influence of the PDO, we compare the time series of the PDO and TC frequency. The calculation indicates that the PDO is positively correlated to the frequency of TC. The correlation coefficient is 0.71, which is significant at the 95% confidence level.
SST in P3 is warmer than that in P2, and the positive SST differences are moving towards the western Pacific. This may be one of the reasons why the average TC frequency increases in P3. Warm SST anomalies during the warm phase of the PDO may cause atmospheric circulation anomalies in the WNP which in turn affect TC activity. During the period 2011-2020, the PDO and the warm tropical Pacific SST anomalies are favorable for increasing the TC frequency in the WNP. The large-scale atmospheric environmental conditions like relative vorticity, relative humidity and vertical velocity play a significant role in TC genesis during P3. 2019): the PDO influences the WNP TC frequency. The TC frequency is the highest in neutral years and the smallest in La Niña years ( Fig. 2), which is different from P1 described in Zhao and Wang (2016): La Niña years are with the highest TC frequency and El Niño years are with the lowest. The TC frequency in El Niño years is no longer the least during P3. This may be caused by the PDO: the correlation between the TC frequency per year and PDO is different during P1 and P3 (Fig. 11e). In addition, the difference between P1 and P3 may be caused by the warmer SST (Fig. 10). The SST gradient between the northern Indian Ocean and WNP influences the South China Sea TC frequency (Li and Zhou, 2014). The WNP SST during the P3 period is higher than that during P2, which is different from P1. Warm SST anomalies in the WNP are observed during P3, which can induce equatorial wind anomalies and thereby further change atmospheric circulations. During the La Niña years of P3, the WNP TC frequency is mainly affected by relative humidity and relative vorticity, as well as by vertical velocity (Figs. 8 and 9).

Discussion and conclusions
In this study, the increase in the WNP TC frequency during 2011-2020 is analyzed based on large-scale atmospheric and ocean conditions. During 2011-2020, approximately 26 TCs occurred on average every year in the WNP and caused catastrophic damages to the Philippines and other Asian countries (Fig. 3). In recent years, several studies found a decrease in TC counts during 1979-2010 (Hsu et al., 2014;Wang, 2016, 2019). It is associated with ENSO and PDO changes. Previous studies mainly focused on the TC frequency before 2015 and did not investigate the change of the whole year TC frequency after 2015. Investigations of the whole year TC frequency are necessary and useful to improve our understanding of TC activity in the WNP.
Based on the JTWC best-track data, this study examines the variations in the WNP TC frequency between January and December of 1979-2020. The WNP TC count shows a significant decline in 1998 and a rise since 2010, resulting in a low cycle (P2: 1998-2010) and two high cycles (P1: 1979-1997 and P3: 2011-2020). Previous studies have analyzed the decrease between P1 and P2. In this study, the analysis is mainly focused on the increased counts of WNP TCs in P3. A difference in the TC number between P2 and P3 can be found (Figs. 1 and 2). Figure 2 shows that the all year difference between P2 and P3 is significant at the 95% confidence level; however the differences in La Niña years, El Niño years and neutral years are lower than the 95% level. P3 is different from P1 described in Zhao and Wang (2016) where La Niña years are the highest and El Niño