1. Introduction
The western North Pacific (WNP) is the most active tropical cyclone (TC) basin on Earth. About 30% of global TCs develop in the WNP and the lowest central pressure of TCs is observed in this area as well. Every year, many TCs in the WNP threaten coastal countries, resulting in enormous socioeconomic damages and loss of lives. Thus, it is important to know when the TC season will start to try to prevent possible damages from the TCs.
Normally, TC activity is modulated by large-scale environmental parameters such as sea surface temperature (SST), vertical wind shear, low-level vorticity, and relative humidity (Gray 1979; McBride and Zehr 1981). These environmental parameters vary seasonally and thus can induce changes in the seasonality of TC activity including the start, end, and length of the TC season. Previous studies on changes in the TC season have largely focused on the Atlantic basin, which has an apparent active hurricane season from June to November. Kossin (2008) examined the Atlantic hurricane season using quantile regression analysis and suggested that this season is possibly getting longer in association with an increase in the SST in the main region where the hurricanes develop. In addition, Karloski and Evans (2016) suggested that low-level vorticity, El Niño–Southern Oscillation (ENSO), and the Pacific decadal oscillation are related to the long-term trends of the hurricane season.
Climatologically, it is known that TCs in the WNP, unlike Atlantic hurricanes, form throughout the year, so the variation in the TC season in the WNP has been of less concern. However, the major period that TCs are active in the WNP has varied year to year. In particular, the start of the TC season in the WNP has shown large interannual variation. For example, no TC formed before July in 2016, whereas nine named TCs occurred before July in 2015. Thus, it is worthwhile to examine the variations of the TC season in the WNP and the large-scale environmental parameters affecting them. In this paper, the analysis of the characteristics of the variation in the TC season in the WNP is described in section 2 and the investigation concerning large-scale environmental parameters affecting the start of the TC season is presented in section 3. The data and methodologies used are described in each section. Finally, a discussion regarding the results of this study occurs in section 4.
2. Variation of the TC season in the WNP
We analyzed the TC season in the WNP based on the annual TC formation dates. The observed TC data were obtained from the best-track data of the Regional Specialized Meteorological Center (RSMC) Tokyo–Typhoon Center. These data include 6-hourly location and intensity information of TCs over the WNP collected since 1951. In this study, the TC formation date is defined as the day when the TC was first labeled as a tropical storm, which is defined as having a maximum sustained wind speed of more than 17 m s−1. In the WNP, the monthly mean number of TCs is at its annual minimum in February (Kim et al. 2014). Thus, we defined a “TC year” in the WNP as the period from 1 February in the corresponding year to 31 January the following year.
Figure 1 shows the box-and-whisker plots for annual TC formation dates from 1951 to 2015 of the TC year. The boxes go from lower to upper quartiles of the TC formation dates in each TC year and whiskers are extended to their 5th and 95th percentiles. Thus the whiskers cover the period that during which around 90% of annual TCs are active. The first and last TC formation dates are also marked with dots. In the calculation of the percentiles, the first and last TC formation dates are regarded as the 0th and 100th percentiles, respectively. On average, the first date of TC formation was 6 April, the 5th percentile was 18 May, the 95th percentile was 27 November, and the last date was 27 December. The 5th to 95th percentile of the TC formation dates properly occurred from June to November, which is known as the active TC season in the WNP. Thus, we defined the TC season in the WNP as the period from the 5th percentile to the 95th percentile of TC formation dates in each TC year. This definition, based on the percentiles, can remove the effect of annual TC counts in the estimation of the TC season. The correlation coefficients between the TC count in each TC year and the start, end, and length of the TC season were −0.07, 0.09, and 0.12, respectively (Table 1), which are statistically insignificant, confirming that the TC season defined in this study is independent of the annual TC counts.

Box-and-whisker plots for the western North Pacific (WNP) tropical cyclone (TC) formation dates for 1951–2015. The boxes cover from the 25th to 75th percentiles (lower to upper quartiles, respectively) and the whiskers extend to the 5th and 95th percentiles, of the annual TC formation dates. The first and last TC formation dates are marked with dots. The lower and upper dashed lines denote the linear trends for the 5th and 95th percentiles, respectively, for the 1951–2015 period.
Citation: Journal of Climate 30, 9; 10.1175/JCLI-D-16-0888.1

Box-and-whisker plots for the western North Pacific (WNP) tropical cyclone (TC) formation dates for 1951–2015. The boxes cover from the 25th to 75th percentiles (lower to upper quartiles, respectively) and the whiskers extend to the 5th and 95th percentiles, of the annual TC formation dates. The first and last TC formation dates are marked with dots. The lower and upper dashed lines denote the linear trends for the 5th and 95th percentiles, respectively, for the 1951–2015 period.
Citation: Journal of Climate 30, 9; 10.1175/JCLI-D-16-0888.1
Box-and-whisker plots for the western North Pacific (WNP) tropical cyclone (TC) formation dates for 1951–2015. The boxes cover from the 25th to 75th percentiles (lower to upper quartiles, respectively) and the whiskers extend to the 5th and 95th percentiles, of the annual TC formation dates. The first and last TC formation dates are marked with dots. The lower and upper dashed lines denote the linear trends for the 5th and 95th percentiles, respectively, for the 1951–2015 period.
Citation: Journal of Climate 30, 9; 10.1175/JCLI-D-16-0888.1
Statistics for the start, end, and length of the TC season over the WNP. The ordinal date is the count of days from 1 February. The asterisk (*) denotes the correlation coefficient is significant at the 99% confidence level. The season start and end are 5th and 95th percentiles of TC formation dates in the TC year and the season length is the period between the season start and end.


The TC season in the WNP showed apparent interannual variability, as seen in Fig. 1. The standard deviation σ and correlation coefficient r for the interannual variation in the parameters for the TC season are given in Table 1. The start and length of the TC season had large interannual variations of more than one month (σstart = 35.1 and σlength = 35.1 days), whereas the variation in the end of the TC season was relatively small (σend = 18.6 days). The start and length of the TC season were highly correlated with each other (rstart,length = −0.86), whereas the end date was less correlated with the start and length of the TC season (rend,start = 0.26 and rend,length = 0.27). These results suggest that the variation in the length of the TC season is mainly associated with the start of the TC season rather than its end. Meanwhile, the start and end of the TC season in the WNP had long-term variations as well as the interannual ones. The linear trends for the start and end of the TC season are plotted on Fig. 1. The start of the TC season was delayed (0.26 day yr−1), whereas the end occurred earlier (−0.23 day yr−1), indicating that the length of the TC season in the WNP is getting shorter (−0.49 day yr−1) (Table 1). However, these long-term trends are not statistically significant at the 90% confidence level. Because of the weak signal in the long-term trends in this study, we paid attention to interannual variation of TC season rather than the long-term variations. In addition, further analysis focused on the variation of the start of TC season, which is a main factor modulating the variation of TC season.
3. Large-scale environmental parameters related to the variation of TC season start
We analyzed large-scale environmental parameters that modulate the start of TC season using the extended reconstructed sea surface temperature version 4 (Huang et al. 2015) and the National Centers for Environmental Prediction (NCEP)–National Center for Atmospheric Research (NCAR) reanalysis dataset (Kalnay et al. 1996). Figure 2 shows the correlation maps between the ordinal date of the TC season start and the large-scale environmental parameters (SST, 850-hPa geopotential height, relative vorticity, and horizontal wind). In the figures, the positive (negative) correlations indicate that the TC season starts later (earlier) when the environments’ anomalies are positive (negative). To examine the factors affecting TC season start, we considered the environmental parameters during the preceding winter (December–February) and spring (March–May) of each TC season.

The correlation coefficients between the start day of TC season and (a),(b) SST, (c),(d) 850-hPa geopotential height, and (e),(f) 850-hPa vorticity and horizontal winds in (left) December–February and (right) March–May. The solid (dashed) line of contour represents positive (negative) correlation coefficients. The contour interval is 0.1 and the zero line is omitted. The light (dark) shading regions indicate statistical significance at the 90% (95%) confidence level. The correlation coefficients with the zonal and meridional winds are illustrated as vectors. The correlation vectors are drawn only if the correlation coefficients with either the zonal or meridional winds are significant at the 95% confidence level.
Citation: Journal of Climate 30, 9; 10.1175/JCLI-D-16-0888.1

The correlation coefficients between the start day of TC season and (a),(b) SST, (c),(d) 850-hPa geopotential height, and (e),(f) 850-hPa vorticity and horizontal winds in (left) December–February and (right) March–May. The solid (dashed) line of contour represents positive (negative) correlation coefficients. The contour interval is 0.1 and the zero line is omitted. The light (dark) shading regions indicate statistical significance at the 90% (95%) confidence level. The correlation coefficients with the zonal and meridional winds are illustrated as vectors. The correlation vectors are drawn only if the correlation coefficients with either the zonal or meridional winds are significant at the 95% confidence level.
Citation: Journal of Climate 30, 9; 10.1175/JCLI-D-16-0888.1
The correlation coefficients between the start day of TC season and (a),(b) SST, (c),(d) 850-hPa geopotential height, and (e),(f) 850-hPa vorticity and horizontal winds in (left) December–February and (right) March–May. The solid (dashed) line of contour represents positive (negative) correlation coefficients. The contour interval is 0.1 and the zero line is omitted. The light (dark) shading regions indicate statistical significance at the 90% (95%) confidence level. The correlation coefficients with the zonal and meridional winds are illustrated as vectors. The correlation vectors are drawn only if the correlation coefficients with either the zonal or meridional winds are significant at the 95% confidence level.
Citation: Journal of Climate 30, 9; 10.1175/JCLI-D-16-0888.1
Figures 2a and 2b show that the ordinal date of the TC season start was positively correlated with the tropical SSTs in the Indian Ocean (IO), South and East China Seas, and tropical central-eastern Pacific (CEP) but negatively correlated with those in the southeastern Philippines Sea. The correlation pattern was strong in the preceding winter (Fig. 2a). The positive correlations in the IO and tropical CEP remained until spring (Fig. 2b). The correlation patterns between the TC season start and SST reflect the variability associated with ENSO (Harrison and Larkin 1998). Thus, we examined the start of the TC season following strong ENSO phases that are selected when the Niño-3.4 anomaly from December to February exceeds one standard deviation from the mean. The selected TC years following strong ENSO events in winter and their start dates of the TC season are shown in Table 2. The mean start date of the TC season following a strong El Niño was 21 June whereas for one following a strong La Niña it was 5 May. It is noted that because of the increasing trend in the TC season start, the TC start dates in early years of El Niño are slightly earlier than those in late years. Despite of the trend, the mean start date of the TC season following El Niño was delayed more than one and a half months from the mean start date after La Niña and the difference is statistically significant at the 99% confidence level based on both the Student’s t test and Mann–Whitney U test. This suggests that the TC season in the WNP is significantly modulated by El Niño–related SST variations.
Start dates of TC season following a strong ENSO event in the preceding winter. The ordinal date is the count of days from 1 February.


Along with the ENSO-like SST variation, atmospheric environmental parameters closely varied with the start of the TC season. Figures 2c–f illustrate the correlation coefficients between the low-level atmospheric environmental parameters in the WNP and the ordinal date of the TC season start. The ordinal date of the TC season start was positively correlated with the 850-hPa geopotential height in the western Pacific and eastern IO during winter (Fig. 2c). In spring, the positive correlation area was smaller but the strong positive correlations were sustained over the WNP TC basin (Fig. 2d). In addition, the correlation coefficients with 850-hPa vorticity and horizontal wind show that anomalous anticyclonic circulation prevailed over the WNP TC basin in spring in association with the late-starting TC season (Fig. 2f). The anomalous high geopotential height and low-level anticyclonic circulation over the WNP provided unfavorable conditions for TC formation in spring, resulting in a delay of the TC season start over the WNP.
Many previous studies have pointed out that the ENSO-related SST anomalies over the tropics affect the anomalous circulation over the WNP basin. Anticyclonic flow develops in the Philippine Sea during the boreal winter of El Niño and is sustained until the following spring by a Rossby wave response to El Niño–induced subsidence over the Maritime Continent and an increase in westerlies over tropical western Pacific, which is forced by the SST warming over the CEP (e.g., Wang et al. 2000; Wang and Zhang 2002; Wu et al. 2010a). In addition to this process in the tropical Pacific sector, the El Niño–related warm SST anomalies over the IO also intensify and sustain the anomalous Philippine Sea anticyclone by a Kelvin wave response that forces surface divergence and descent over the WNP (e.g., Watanabe and Jin 2002; Lau and Nath 2003; Annamalai et al. 2005; Yang et al. 2007; Wu et al. 2010b; Zhan et al. 2011a,b; Zhan et al. 2014; Yuan et al. 2012). Thus, the correlation patterns shown in Fig. 2 can be attributed to a similar physical process as the aforementioned ENSO-related variations. In summary, when the warm SST anomalies occurred in the CEP and IO in the preceding winter, anomalous high pressure was induced over the vicinity of the Philippines Sea by the atmosphere–ocean coupled processes. The anomalous high pressure over the Philippines Sea was maintained until spring, suppressing TC formation in the spring and causing a late start of the TC season.
4. Summary and discussion
In this study, the variation of TC season in the WNP was analyzed based on the percentile analysis on the TC formation days in each TC year, which is defined from 1 February in the corresponding year to 31 January the following year. The TC season is defined as the period from the 5th to the 95th percentiles of annual TC formation days. The results show that the length of TC season was highly dependent on its start, which has large interannual variation. The start of the TC season was significantly associated with the ENSO-related variations in oceanic and atmospheric environmental parameters during the preceding TC season. When the El Niño–like tropical SST anomalies (i.e., anomalous warming over the IO and CEP) occurred in the preceding winter and anomalous high pressures and anticyclonic circulations were induced and maintained around the WNP TC basin until the spring, TC formations over the WNP were suppressed in the spring, resulting in the late start of the TC season.
The statistics for the start of the TC season following a strong El Niño and La Niña also supported the relationship between the ENSO-related SST variation in the preceding winter and the start of the TC season (Table 2). It is notable that the start of the 2016 TC season in the WNP was extremely late. A named TC was not reported until 3 July by the RSMC-Tokyo Typhoon Center. Meanwhile, the El Niño of 2015–16 had the strongest El Niño events observed since 1950 (L’Heureux et al. 2017). Thus, the relationship between the strong El Niño and late start of the TC season analyzed in this study was valid for the 2016 TC season of the WNP.
As mentioned in section 2, the long-term trends show that the length of TC season has been slightly shortening as the start of the season is occurring later and the end is occurring earlier. Although the trends are very weak and statistically insignificant, it is notable that the reduction in the TC season length in the WNP is different from the positive trends in the Atlantic hurricane season length reported by Kossin (2008). They interpreted that the increases in local SST warming in the hurricane-developing region might affect the longer hurricane season. However, the simple linkage between the long-term trends in the local SST and length of TC season may not be valid for the WNP. In fact, the local SSTs in the WNP TC basin have also increased in the last several decades, whereas the length of the TC season has shortened slightly. From the results of this study, we speculate that the long-term variation of the TC season over the WNP might be associated with the remote SST variations over the IO and CEP and related atmospheric conditions of the WNP TC basin. Unfortunately, because of the large uncertainty in the insignificant long-term variations of the TC season over the WNP, the environmental factors affecting the possible long-term trends of the WNP TC season remain unknown.
Acknowledgments
This work was funded by the Korea Meteorological Administration Research and Development Program under Grant KMIPA 2015-1100.
REFERENCES
Annamalai, H., P. Liu, and S.-P. Xie, 2005: Southwest Indian Ocean SST variability: Its local effect and remote influence on Asian monsoons. J. Climate, 18, 4150–4167, doi:10.1175/JCLI3533.1.
Gray, W. M., 1979: Hurricanes: Their formation, structure and likely role in the tropical circulation. Meteorology over Tropical Oceans, D. B. Shaw, Ed., Royal Meteorological Society, 155–218.
Harrison, D. E., and N. K. Larkin, 1998: El Niño–Southern Oscillation sea surface temperature and wind anomalies 1946–1993. Rev. Geophys., 36, 353–399, doi:10.1029/98RG00715.
Huang, B., and Coauthors, 2015: Extended Reconstructed Sea Surface Temperature version 4 (ERSST.v4): Part I. Upgrades and intercomparisons. J. Climate, 28, 911–930, doi:10.1175/JCLI-D-14-00006.1.
Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437–471, doi:10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.
Karloski, J. M., and C. Evans, 2016: Seasonal influences upon and long-term trends in the length of the Atlantic hurricane season. J. Climate, 29, 273–292, doi:10.1175/JCLI-D-15-0324.1.
Kim, H.-S., G. A. Vecchi, T. R. Knutson, W. G. Anderson, T. L. Delworth, A. Rosati, F. Zeng, and M. Zhao, 2014: Tropical cyclone simulation and response to CO2 doubling in the GFDL CM2.5 high-resolution coupled climate model. J. Climate, 27, 8034–8054, doi:10.1175/JCLI-D-13-00475.1.
Kossin, J. P., 2008: Is the North Atlantic hurricane season getting longer? Geophys. Res. Lett., 35, L23705, doi:10.1029/2008GL036012.
Lau, N.-C., and M. J. Nath, 2003: Atmosphere–ocean variations in the Indo-Pacific sector during ENSO episodes. J. Climate, 16, 3–20, doi:10.1175/1520-0442(2003)016<0003:AOVITI>2.0.CO;2.
L’Heureux, M. L., and Coauthors, 2017: Observing and predicting the 2015–16 El Niño. Bull. Amer. Meteor. Soc., doi:10.1175/BAMS-D-16-0009.1, in press.
McBride, J. L., and R. Zehr, 1981: Observational analysis of tropical cyclone formation. Part II: Comparison of non-developing versus developing systems. J. Atmos. Sci., 38, 1132–1151, doi:10.1175/1520-0469(1981)038<1132:OAOTCF>2.0.CO;2.
Wang, B., and Q. Zhang, 2002: Pacific–East Asian teleconnection. Part II: How the Philippine Sea anomalous anticyclone is established during El Niño development. J. Climate, 15, 3252–3265, doi:10.1175/1520-0442(2002)015<3252:PEATPI>2.0.CO;2.
Wang, B., R. Wu, and X. Fu, 2000: Pacific–East Asian teleconnection: How does ENSO affect East Asian climate? J. Climate, 13, 1517–1536, doi:10.1175/1520-0442(2000)013<1517:PEATHD>2.0.CO;2.
Watanabe, M., and F.-F. Jin, 2002: Role of Indian Ocean warming in the development of Philippine Sea anticyclone during ENSO. Geophys. Res. Lett., 29, 1478, doi:10.1029/2001GL014318.
Wu, B., T. Li, and T. Zhou, 2010a: Asymmetry of atmospheric circulation anomalies over the western North Pacific between El Niño and La Niña. J. Climate, 23, 4807–4822, doi:10.1175/2010JCLI3222.1.
Wu, B., T. Li, and T. Zhou, 2010b: Relative contributions of the Indian Ocean and local SST anomalies to the maintenance of the western North Pacific anomalous anticyclone during the El Niño decaying summer. J. Climate, 23, 2974–2986, doi:10.1175/2010JCLI3300.1.
Yang, J., Q. Liu, S.-P. Xie, Z. Liu, and L. Wu, 2007: Impact of the Indian Ocean SST basin mode on the Asian summer monsoon. Geophys. Res. Lett., 34, L02708, doi:10.1029/2006GL028571.
Yuan, Y., S. Yang, and Z. Zhang, 2012: Different evolutions of the Philippine Sea anticyclone between the eastern and central Pacific El Niño: Possible effects of Indian Ocean SST. J. Climate, 25, 7867–7883, doi:10.1175/JCLI-D-12-00004.1.
Zhan, R., Y. Wang, and X. Lei, 2011a: Contributions of ENSO and East Indian Ocean SSTA to the interannual variability of northwest Pacific tropical cyclone frequency. J. Climate, 24, 509–521, doi:10.1175/2010JCLI3808.1.
Zhan, R., Y. Wang, and C.-C. Wu, 2011b: Impact of SSTA in East Indian Ocean on the frequency of northwest Pacific tropical cyclones: A regional atmospheric model study. J. Climate, 24, 6227–6242, doi:10.1175/JCLI-D-10-05014.1.
Zhan, R., Y. Wang, and L. Tao, 2014: Intensified impact of east Indian Ocean SST anomaly on tropical cyclone genesis frequency over the western north Pacific. J. Climate, 27, 8724–8738, doi:10.1175/JCLI-D-14-00119.1.