Abstract

The Northwest Pacific tropical cyclone (TC) intensification is classified into rapid intensification (RI), normal intensification (NI), and slow intensification (SI) categories. The initial location and intensity, the preceding intensity change, the motion direction, the occurrence month, and the intensification duration time are all found to differ for RI cases compared with NI and SI cases. The dependence of RI, NI, and SI on environmental conditions is further examined statistically by using the intensification rates of named TCs, for the 21-year period 1995–2015, obtained from JTWC best track data, and the environmental conditions derived from the ERA-Interim reanalysis data and GODAS high-resolution global ocean analysis data. It was found that deep-layer and upper-mid vertical wind shear (VWS), upper-level outflow, sea surface temperature (SST), and ocean heat content (OHC) are statistically different among RI, NI, and SI both before and during intensification. RI is enhanced by weaker and decreasing VWS, warmer oceans, and stronger and increasing outflow. In contrast, SI typically occurs with larger and increasing VWS, cooler oceans, and weaker, static outflow. The impacts of low-level VWS and net moisture inflow are only significantly different between RI and SI and between NI and SI, but not between RI and NI. Another key finding is that increased upper-level outflow and decreased VWS are important precursors and hence are possible predictors, of RI onset. The direction of upper-level outflow affects TC intensification, with NW and NE outflow being more favorable for TC RI than SE and SW outflow.

1. Introduction

The operational prediction of tropical cyclone (TC) intensity change remains a major challenge. Operational forecasts often are required to distinguish rapid intensification (RI) and slow intensification (SI) from normal intensification (NI) [13]. However, during RI, the forecast intensification rate (IR) typically is far lower than the observed IR, whereas it is much higher than the observed IR in the SI cases [3, 4]. In some operational TC forecast centers (Figure 1), when forecasters are unsure of how rapidly TCs will intensify, the default forecast issued is that the TC will undergo NI [3]. The difficulty of distinguishing between RI, NI, and SI during the forecast is due partly to the challenge of quantifying the relationship between the predictors and the IR. Most previous studies focused solely on comparing environmental conditions between RI and non-RI cases. Some of these studies revealed statistically significant differences in the environmental conditions between RI and non-RI cases [5, 6]. However, in these studies, the non-RI category does not distinguish between NI, SI, and weakening storms. Hendricks et al. [7] classified the intensification into two categories: rapidly intensifying and intensifying and showed that the environmental conditions, including sea surface temperature (SST), and other environmental parameters associated with these two categories are similar. The present study further explores the relationship between environmental conditions and TC IR and, in contrast to previous studies, also classifies the IR into the three categories of RI, NI, and SI. The main motivation for this study is toward addressing the need for, and challenges of, distinguishing between the various intensification phases present during the operational forecast period.

Most of the earlier studies have focused on hurricanes in the North Atlantic or the East Pacific basins [5, 6, 8]. However, this study concentrates on TCs over the Northwest Pacific Basin (NWP). The NWP is controlled by the monsoon background, with about 70% of all TC genesis occurring in the monsoon trough [9, 10]. This is very different from other ocean basins. For example, in the East Pacific and North Atlantic basins, the main genesis source is African Easterly Waves [11]. Hence, the relative importance of the individual predictors for TC intensification is likely to be different in various basins [6, 12]. Rios-Berrios and Torn [13] also suggested it is problematic to consider all basins jointly, because the relative importance of environmental conditions for TC intensity change depends on the ocean basin.

The present study seeks to address an important, but as yet unanswered, question: “Do significant differences in environmental factors, including atmospheric conditions and ocean heat conditions, exist between rapidly, normally, and slowly intensifying TCs in the NWP?” In this study, a range of environmental factors is examined, including vertical wind shear (deep-layer, upper-mid, mid-low, and low-level), net moisture inflow, upper-level outflow (strength and direction), and ocean heat conditions (both SST and ocean heat content (OHC)). Earlier, limited studies of the NWP TC basin either examined different TC categories [14] or focused only on one environmental factor [12].

Specifically, the following environmental factors and their relationships with RI, NI, and SI are assessed. Environmental vertical wind shear (VWS) is one of the most important factors influencing TC intensity change, and it is widely accepted that strong VWS generally inhibits TC intensification [5, 7, 12, 1518]. VWS also provides probability estimates of RI in dynamical-statistical models [5, 6, 16, 19] and can assist in determining if, and when, RI will occur [2022]. The deep layer, 200–850 hPa VWS, is commonly used, whereas other studies showed that the VWS at different levels can also impact TC intensity change [12, 2325]. Hence, it is too simplistic to assess the impact of VWS on TC intensity change by using a VWS calculated from just two levels [2629]. In this study, the VWS between different vertical levels are used to analyze the impact of VWS on the three TC intensification rates, RI, NI, and SI.

Environmental moisture is also suggested as a possible important factor for TC intensification. Both theoretical and modeling studies showed that high environmental moisture favors TC intensification [30, 31], whereas dry air intrusion can cause TC decay [21, 32]. The environmental relative humidity (RH) of the lower and middle troposphere is higher for RI cases than non-RI cases [57]. However, other studies suggested that moisture might have a negative or no significant impact on TC intensification [3335]. Therefore, the impact of environmental moisture on TC IR requires further, more detailed, examination. In this study, the impact of environmental moisture on the TC IR is studied by using the low-level to mid-level net moisture inflow, defined as the sum of net moisture inflow from all sides of a box surrounding the TC center and extending from the surface to 500 hPa. This quantity is used to represent the total low-level to mid-level moisture transport from the environment into the TC core.

Due to the asymmetric structure [36, 37], low Richardson number [38, 39], and weaker inertial instability [4043] of the TC upper-level outflow, the outflow can easily interact with the upper-level larger-scale environment and the inner-core of a TC, thereby playing a mediating role between the environment and the storm core [41]. As such, it can affect the secondary circulation and, consequently, the TC intensity [37, 41, 4349]. The upper-level radial outflow is much stronger [6, 37, 43] and usually is more concentrated close to the TC center for intensifying cases than for nonintensifying cases [43]. In addition, TC intensification is associated with the multiple upper-level outflow channels, and TCs usually weaken when one of the efficient upper-level outflow channels is cut off [44, 5053]. Kaplan and DeMaria [5] and Shu et al. [14] noted that upper-tropospheric flows for RI cases are more easterly than for non-RI cases. Rappin et al. [41] showed that the ventilation of outflow in north and northeast directions minimizes the hurricane energy expenditure, leading to larger intensification rates. Hence, the direction of TC outflow possibly modifies the effect of outflow on TC intensity change. In this study, the significant differences in the impacts of both the strength and direction of outflow are examined, for all three different TC intensification categories. As far as the authors are aware, this relationship between outflow direction and TC IR has not previously been statistically analyzed.

The effect of the ocean heat conditions on TC intensity change has been addressed for decades [5457]. SST is widely known to affect TC intensification [5, 6, 58]. Numerous supporting case studies exist [21, 5961]. Statistical analyses [8, 12] showed that TCs typically intensify (weaken) over warmer (colder) ocean eddy regions with higher (lower) SST. Some studies have shown that the OHC, measured from the surface to the depth of the 26°C isotherm, is important for TC intensity change [59], and that SST alone is not always a good measure of the total heat content that affects TC intensity change [62]. This hypothesis is further supported by studies showing that the use of OHC can help reduce TC intensity forecast errors in the statistical model, SHIPS [63, 64]. In the present study, both SST and OHC are used to determine the impact of ocean heat conditions on RI, NI, and SI.

In summary, the main question that needs to be addressed in this study is whether there are significant differences in large-scale atmospheric environmental conditions, such as VWS at different vertical levels, net environmental moisture inflow, strength and direction of upper-level outflow, and ocean heat conditions (SST and OHC) between RI, NI, and SI, for NWP TCs. Additionally, a unique aspect of this study is the use of a relatively long period of data to derive statistically robust results, as all named NWP TCs during the 21-year period from 1995 to 2015 are included. The findings from this study are relevant to improving TC forecast skill and to increasing our understanding of the TC intensification processes.

Section 2 describes the data and provides definitions of the environmental conditions. Section 3 defines the RI, NI, and SI categories and shows the distribution of TC intensity changes in the NWP basin, and Section 4 describes and discusses the main results. Finally, Section 5 summarizes the findings of this study.

2. Data and Methodology

2.1. Data

The Joint Typhoon Warning Center (JTWC) best track data, including the TC center position and maximum 1 min mean sustained 10 m wind speed (Vmax), at 6 h intervals for all named TCs over the NWP basin from 1995 to 2015, were used to calculate TC intensity change and to define the TC environment.

The 6-hourly wind and specific humidity fields and SST data from 1995 to 2015 are from the European Center for Medium-Range Weather Forecasts (ECMWF) interim reanalysis website (https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era-interim) [65]. The same ECMWF reanalysis data were widely used in published statistical studies to investigate the relationship between the environmental conditions and the tropical cyclone intensity change (e.g., [25, 66]). Schenkel and Hart show that for the ECMWF interim reanalysis, most of the TC position errors are less than 150 km in the Northwest Pacific basin [67]. The study of Hodges et al. showed that for the ECMWF interim reanalysis [68], the majority of the TCs in the reanalysis have a mean position difference from the best track that is less than 2° (∼220 km). They also showed that the peak of the position error distribution is at about 110 km. Hence, most of the TC position errors in the ECMWF reanalysis are less than 110 km [68]. In the studies of Schenkel and Hart [67] and Hodges et al. [68], the earlier ECMWF interim reanalysis was evaluated. In our study, we used the most recent 20-year ECMWF reanalysis (1995–2015). Because this period is more recent and therefore the quality of the ECMWF interim reanalysis data should be improved, the TC position errors of the ECMWF reanalysis in our study should be smaller than that of Schenkel and Hart [67] and Hodges et al. [68]. In addition, the horizontal grid spacing of the ECMWF interim reanalysis is 0.5° longitude × 0.5° latitude. In other words, the smallest resolvable scale is about ∼100 km, comparable to the peak position error in Hodges et al.’s study [68]. Moreover, the computational domain which is used to calculate the environmental conditions in this study is relatively large compared to the TC position error (∼110 km). Based on all of the above, the TC position error in the ECMWF interim analysis is not expected to impact our results. So the reanalysis data are used to calculate different level VWS, upper-level outflow, net moisture inflow, and SST around the TC center, during the intensification period.

To study the impact of OHC, the high-resolution global ocean analysis sea temperature data (34 levels up to 5000 m, 12 km grid spacing) were provided by U.S. GODAE (U.S. Global Ocean Data Assimilation Experiment) (http://usgodae.org/). GODAE was shown by early studies to provide high-quality ocean analysis data [69] and was used to analyze the upper-ocean heat condition (e.g., [70]). These data are used to calculate the OHC around the TC center in this study. Although the GODAE data have high spatial resolution, they are only available after 2007. Therefore, only a 9-year worth of data (2007–2015) were used. Given the statistical relationship between intensification rates and other environmental variables was studied using 20-year data, caution should be taken to intercompare the OHC impact with other environmental variables.

2.2. Definitions of TC Environmental Conditions
2.2.1. Vertical Wind Shear

The environmental vertical wind shear, VWS, is defined as the wind vector difference between two levels as shown in the following equation:where the and variables are annular-averaged horizontal wind vector at level 1 and level 2. In this calculation, the annular is between 200 and 800 km from the storm center. Therefore, the inner-core region of the storm is removed or filtered when the environmental shear is calculated. The subscripts “lev1”and “lev2” indicate two pressure levels. This method of defining and calculating the vertical wind shear follows the studies [14, 19]. We calculate 4 types of wind shear , , , and as follows. The deep-layer shear is calculated. It is commonly used to diagnose the relationship between TCs intensity change and environmental vertical wind shear (e.g., [18, 19]). 850 hPa, 500 hPa, and 200 hPa are commonly used to analyze the low-level, mid-level, and upper-level circulations of tropical cyclone. Therefore, in addition to , and are also used because early studies suggest these shears can also impact TC intensity change [19, 71]. Moreover, both the statistical and idealized simulation studies found that the low-level shear could have a negative effect on TC intensification [12, 24]. In particular, the study by Wang et al. [12] showed that the low-level shear between 850 and 1000 hPa can have a larger negative impact on TC intensity change than VWS at other levels during active seasons. So the low-level shear is also calculated in this study.

2.2.2. Net Moisture Inflow

The net inflow of environmental moisture is used to assess the impact of environmental moisture on TC intensity change. The net moisture inflow flux is defined aswhere are the net moisture inflows from the eastern, southern, western, and northern boundaries, respectively, of a 8° × 8° latitude-longitude box around the TC center. Hence, Quv is the total net moisture transport to the TC core from the surface to 500 hPa. The units of Quv are kg m−2 s−1. In equation (2),where q is the specific humidity, u and are the zonal wind and meridional wind components, respectively, and is 9.8 m s−2.

2.2.3. Upper-Level Outflow

The strongest TC outflow is, on average, concentrated at about 200 hPa and around 500 km from the TC center [38, 43]. In this study,where from Green’s law, the upper-level outflow strength is defined as the area-averaged divergence over the 8° × 8° latitude-longitude box around the TC center, at 200 hPa.

For simplicity, the impact of outflow direction on IR is investigated only for cases with a single-outflow channel. A single-outflow channel is defined as occurring when the quadrant-mean wind within 500 km from the storm center is outward and greater than 6 m s−1. In the Northern Hemisphere, the anticyclonic outflow of TCs mainly has four radial directions: northwest (NW), northeast (NE), southwest (SW), and southeast (SE). Therefore, in this study, only the NW, NE, SW, and SE quadrants are considered.

2.2.4. SST and OHC

The OHC definition follows Buckley et al. [72] except that we choose the 26°C isotherm depth as the integral depth during the vertical integration. This replacement was according to early studies (e.g., [73]) which showed that the TC genesis and/or TC intensity change can be sensitive to the depth of the 26°C isotherm. Accordingly, the areal density of the OHC is defined aswhere is the seawater density, is the specific heat of sea water, is the depth of the 26°C sea temperature contour, and T(z) is the sea temperature at the depth z. The unit of the areal density of ocean heat content is kJ cm−2. Hence, OHC is the total ocean heat content per square centimeter, from the depth of the 26°C contour sea temperature to the sea surface. In this study, the area-averaged SST and OHC over the 2° × 2° square box around the TC center is used to calculate the ocean heat conditions.

2.3. Significant Difference Tests

The testing of significant differences between mean values was carried out using a standard bootstrapping method [74, 75]. This method was used to obtain confidence levels for the differences in the sample means of the various parameters in this study. As an example, it was used to compare the mean intensification rates for RI and SI in deep-layer VWS environments. Here, confidence levels of 90% and 95% are used, and the number of resamples in this study is 10,000. The bootstrap technique was used in early studies by Wang and Bishop [75] and Hamill et al. [76].

3. Definitions and Distributions of RI, NI, and SI in the NWP Basin

3.1. Definitions of RI, NI, and SI

The frequency and cumulative frequency distributions of 12 h intensity changes for all 499 named TCs over the NWP basin during 1995–2015 are shown in Figure 1. As this study focuses on TC intensification, only intensifying cases are chosen. The IR is classified into 3 categories: rapid intensification (RI), normal intensification (NI), and slow intensification (SI). These are formally defined in Table 1.

Note that in earlier studies, RI is usually defined as the 95th percentile of all 24 h TC intensity changes, i.e., approximately 30 kt or greater [5, 12, 14]. These previous definitions did not include the 6 h TC intensity change. As shown later, the average intensification duration time of the continuous RI process found in this study is 18.8 hr. Therefore, using a 24 h period to define the intensification period might also include a weakening phase. Hence, in this study, a 12 h period is used to define RI. In addition, to ensure the selection of continuously intensifying cases, where no weakening phases are included, the definition of intensification category is also constrained by the 6 h intensity change. Specifically as shown in Figure 1, RI is defined as a 12 h intensity change ≥20 kt (i.e., the 95th percentage of TC intensity change for 12 hr) and a 6 h intensity change ≥5 kt. NI is defined as a 12 h intensity change between 10 kt and 20 kt (i.e., between 85th and 95th percentage of TC intensity change for 12 hr), and a 6 h intensity change is between 5 and 10 kt. SI requires that the 12 h intensity change must lie between 0 kt and 10 kt (i.e., between 60th and 85th percentage of TC intensity change for 12 hr), and the 6 h intensity change between 0 kt and 5 kt. Finally, the total numbers of RI, NI, and SI cases examined in this study are 763, 1345, and 1218, respectively.

3.2. Distributions of RI, NI, and SI

In this section, the statistical distribution of RI, NI, and SI with respect to initial location, initial intensity, moving direction, occurrence month, and intensification duration time are discussed. The initial TC center locations and intensities of the three intensity categories are shown in Figure 2. Hereafter, TC intensities are classified into six categories, tropical depression (TD, 10.8 m s−1 ≤ Vmax ≤ 17.1 m s−1), tropical storm (TS, 17.2 m s−1 ≤ Vmax ≤ 24.4 m s−1), severe tropical storm (STS, 24.5 m s−1 ≤ Vmax ≤ 32.6 m s−1), typhoon (TY, 32.7 m s−1 ≤ Vmax ≤ 41.4 m s−1), strong typhoon (STY, 41.5 m s−1 ≤ Vmax ≤ 50.9 m s−1), and super typhoon (Super TY, Vmax ≥ 51 m s−1) based on the CMA (China Meteorological Administration) criteria of 2006 [77]. Only named TCs are chosen in this study and TDs are not included.

From Figure 2, the TC data reveal that most RI cases occur at latitudes between 10°N and 25°N, whereas most NI cases occur between 5°N and 30°N, and SI cases occur in an even larger area between 5°N and 35°N. RI cases mainly are concentrated east of the Philippines and in the northern South China Sea. For each intensification category, the distribution of initial intensity as a function of TC intensity grade is shown in Figure 3. For RI cases, the initial intensity most frequently falls in the TY grade. For NI and SI cases, the initial intensity most commonly falls in the TS intensity category. The combined frequency of TY, STY, and super TY categories is about 64% for RI but only around 33% for both NI and SI. Overall, the initial intensity of the RI cases generally is much greater than that of NI and SI. However, the RI distribution over the TY, STY, and super TY intensities in Figure 3 shows that fewer RI cases are initiated from the STY and super TY intensities than from TYs. This possibly is because TYs are farther from their MPI (maximum potential hurricane intensity [78]) than STYs and super TYs and hence have greater potential to intensify [79].

The frequency distributions of the TC moving directions for the three intensification rate categories (SI, NI, and RI) are shown in Figure 4. Overall, most of the intensification cases (about 75%) move northwest (including WNW and NNW), whereas only about 16% and 7% of the intensification cases move northeast (including NNE and ENE) and west-southwest, and very few cases move southward (including SSE and SSW) or east-southeastward. It is also found that the percentage of west-northwestward motion is higher for RI cases than NI and SI cases, whereas the percentage of north-northwestward moving TCs is lower for RI cases than NI and SI cases. That is, RI cases tend to move WNW more than NI and SI cases.

Most of the RI cases (about 61%) occur in August, September, and October, with a peak in September (Figure 5), but about 56% of NI and SI cases occur mainly in July, August, and September, with an earlier peak in August. As expected, very few TC intensification cases occur in January, February, and March.

As shown in Figure 6, the mean intensification duration times for a continuous RI, NI, and SI process are 18.8 hours, 18.3 hours, and 14.6 hours, respectively. About 70% of the RI and NI cases intensify for 12–18 hours and only about 30% of the RI and NI cases continue for 24 hours or longer (Figure 6). Very few TCs can intensify rapidly or normally for longer than 36 hours. Most SI cases (about 71%) intensify only for 12 hours. Overall, RI has the longest intensification period and SI has the shortest. As discussed earlier, these results form the basis to define different intensification categories using the 12 h period rather than the 24 h period.

4. Results of Dependence of IR on Environmental Factors

Section 3.2 shows the statistical distributions of initial location, initial intensity, moving direction, occurrence month, and the intensification duration time, which are different among the RI, NI and SI cases. This finding warrants an assessment of their dependence on environmental conditions.

As discussed in Section 1, the selection of the environmental conditions is based on factors that are known to strongly influence TC intensity change. The variables selected are deep-layer shear , the upper-mid shear and mid-low- shear , and low-level shear , the net environmental moisture inflow, the upper-level outflow, and ocean heat condition (including both SST and OHC).

4.1. Environmental Vertical Wind Shear

Figure 7 shows the VWS before, during, and after the intensification for the three intensification categories. The significant differences of the mean magnitudes between the three intensification categories are tested using a bootstrap technique. Overall, the deep-layer shear and the upper-mid shear are significantly (greater than 95% confidence level) lower before, during, and after the RI process than the NI and SI processes. This is not the case for the low-level shears. Notably, though not statistically significantly different, the average low-level shear is slightly higher for RI than NI. For mid-low VWS500–850, consistent significant differences exist between RI and SI, but the differences between RI and NI are significant only starting 12 hours before the onset of the intensification and during the intensification (greater than 90% confidence level). These results suggest that RI is distinguished from NI and SI primarily by deep-layer, upper-level shear, and mid-low VWS rather than the low-level VWS. In addition, , , VWS500–850, and are also significantly (greater than 90% confidence level) lower during the NI processes than the SI processes, except during the onset period for VWS500–850. These results show that SI can be distinguished from NI by all VWS.

On average, there are clear decreasing trends of VWS before the onset of RI, except for (Figure 7). The decreasing trends commence at least 24 hours prior to the onset of RI. For example, VWS200–850 and VWS200–500 reach their minima at the onset of RI (Figures 7(a) and 7(b)). decreases to a minimum at 12 hours prior to the onset of RI (Figure 7(c)). Further diagnostics (not shown) reveal that the decreasing trends are accompanied by a reduction in vortex tilt.

For NI and SI, the VWS usually exhibits slight changes before the onset of RI. In particular, VWS500–850 shows a slight increasing trend for NI and SI before the onset of RI. In other words, in addition to the magnitude, the variation of VWS differs for the three TC intensification categories and therefore can serve as one of the preceding signals for predicting the three intensification categories.

Figure 8 shows the frequency distributions of the three intensification categories as a function of VWS. For VWS850–1000, the percentage peak for all three intensification categories corresponds to the same shear value of 1-2 m s−1. For shear at other levels, VWS200–850, VWS200–500, and VWS500–850, the percentage peak of RI corresponds to a smaller VWS shear than NI and SI. Specifically, the percentage peak of RI is at 2–6 m s−1 for VWS200–850, 2–4 m s−1 for VWS200–500, and about 2-3 m s−1 for VWS500–850. In comparison, the percentage peak for NI and SI is about 2-3 m s−1 larger than RI. NI and SI are also further compared. For VWS200–850 and VWS200–500, the percentage peak for NI is at a smaller shear than SI, whereas for VWS500–850 and VWS850–1000, they both peak at the same shear. Figure 8 also shows that the shape of the frequency distribution as a function of VWS is different among all three intensification categories. The distribution of RI is sharper and has a higher peak value than NI and SI for VWS200–850, VWS200–500, and VWS500–850, followed by NI. However, for VWS850–1000, the distributions among the three categories are in general similar. And for VWS500–850, the percentage peaks of NI and SI are both at 2-3 m s−1; the differences just are that that the percentage and peak value is slightly higher for NI than SI when VWS500–850 is less than 4 m s−1. These results further suggest that deep-layer and upper-mid VWS of the three intensification categories are more significantly different and play more important roles in distinguishing the three intensification categories.

Figure 9 shows the relative occurrence probability of the three intensification categories in each VWS bin. Evidently, the RI occurrence probability decreases linearly with increasing deep-layer VWS () and upper-mid VWS (). It decreases only marginally with increasing mid-low VWS (). RI probability distribution is almost flat when the low-level shear () is less than 3 m s−1 and decreases with increasing low-level shear when is greater than 3 m s−1. This result suggests that the RI occurrence probability is more sensitive to the deep layer and midupper VWS. The NI mean probability decreases only with increasing and with the increasing if the is higher than 4 m s−1. It is almost flat across the bin. However, it increases slightly with the increasing . These results indicate that the NI occurrence probability is more sensitive to the upper-mid and low-level VWS than VWS at other levels. The mean occurrence probability of SI increases especially with an increasing deep-layer and low-level VWS, suggesting its occurrence is more sensitive to the variation of the shear at these levels.

Notably, the error bar range is much larger when VWS is stronger, especially for , , and . This result shows that the occurrence probabilities of SI, NI, and RI at large shear are with high uncertainty, suggesting that the probabilities of SI, NI, and RI are more difficult to estimate under strong environmental VWS. These results are consistent with the study of Zhang and Tao [20] that suggests that the larger the VWS, the larger the uncertainty of the predictability of TC intensity. Although the RI occurrence probability shows a decreasing trend for larger shear values for all levels, and the error bars of occurrence probabilities are much larger in low-level than other levels, especially when low-level shear is larger than 4 ms−1. This result together with those from Figures 7 and 8 show that RI, NI, and SI are overall better delineated by the deep-layer and upper-mid shears than the mid-low and low-level shears. The role of mid-low and low-level shears on distinguishing the three intensification categories is mixed. These results are in general consistent with earlier studies showing the VWS at different levels can impact TC intensity change differently [12, 2325].

4.2. Environmental Moisture Condition

Earlier studies on the impact of environmental moisture on the TC intensification produce mixed results [6, 7, 33, 35]. This study reexamines the impact of the net environmental moisture in distinguishing RI, NI, and SI. The net moisture inflow (Quv) which measures the moisture transport from the environment to the TC core is used. For each intensification category (SI, NI, and RI), the net moisture inflow changes before and during the intensification are shown in Figure 10(a). The net moisture inflow apparently increases before, during, and after the RI and NI processes. The increasing rate of the net moisture inflow is much smaller for SI. On average, the net moisture inflow is slightly higher for RI followed by NI and then SI. The difference in the net moisture inflow between SI and RI (NI) is significant at the 90% level just after the onset of the intensification. However, the difference in net moisture inflow between RI and NI is not statistically significant.

Figure 10(b) shows the percentage distribution of the three intensification categories as a function of net moisture inflow. The percentage peak of RI corresponds to a net moisture inflow of 20–30 kg m−2 s−1, about 10 kg m−2 s−1 larger than that of SI. As expected, the distribution of NI is in between showing a near flat top between 10 and 30 kg m−2 s−1. Results from Figure 10 are in general consistent with early studies which show that the environmental moisture of the lower and middle troposhere is higher for RI cases than non-RI cases [6, 7]. However, when the intensification cases are further delineated among RI, SI, and NI, this study finds the net moisture inflow before, during, and after the intensification is not significantly higher for RI than NI and the moisture inflow of RI and NI is only significantly higher than SI after the intensification.

4.3. Upper-Level Outflow
4.3.1. Strength of Upper-Level Outflow

The upper-level outflow strength is defined as the area-averaged divergence over a (8° × 8°) box around the TC center, at 200 hPa. Figure 11(a) shows the variation of the mean upper-level outflow strength 24 hours before and after the onset of the intensification for the three intensification categories. On average, the upper-level outflow strength increases prior to the RI and NI onset and continues strengthening thereafter. In contrast, there is little variation of outflow strength for SI cases. Moreover, the outflow is strongest in RI, followed by NI. The outflow strength difference between RI and NI is statistically significant at the 90% level. SI shows the weakest outflow, with SI being statistically significant weaker than RI and NI at the 95% level.

Figure 11(b) shows the frequency distribution of the three intensification categories as a function of outflow strength. The frequency peak of RI occurs at a larger outflow strength value than NI and SI. Specifically, the frequency peak of RI is at 0.9 × 10−5 s−1 while those of NI and SI are at about 0.6 × 10−5 s−1. Although NI and SI have similar peak outflow strength, the frequency distribution suggests that SI occurs more frequently than NI when the outflow strength is smaller than 1.2 × 10−5 s−1 and vice versa. Consistent with Figure 11(a), these results further suggest that the occurrence of RI is associated with stronger outflow strength than both NI and SI, and SI usually occurs under the environment with the weaker outflow. In other words, the three intensification categories occur under significantly different outflow strengths. These results are consistent with the previous studies that the upper-level radial outflow is much stronger for intensifying cases than nonintensifying cases [43]. But our results focus only on the intensification cases without including weakening cases and examine the association of the occurrence of difference intensification rates with the strength of the outflow.

4.3.2. Direction of Upper-Level Outflow

For each intensification category, the distributions of the occurrence probabilities in the four outflow directions (NW, SW, NE, and SE) are shown in Figure 12(a). In brief, different outflow directions favor different intensification rates. RI occurs more frequently under NW and NE outflow environments than under SW and SE outflow environments. In contrast, NI occurs more frequently under NW and SW outflow environments and occurs least frequently under a SE outflow environment. SI most likely occurs within SE outflow followed by the SW, NE, and NW outflow.

Figure 12(b) shows the variation of the mean intensification rate in the four outflow directions, which further underlines the sensitivity of the intensification rate to both the outflow direction and strength. For outflow strength greater than 8–10 m s−1, the intensification rate associated with the NW and NE outflows is higher than that with the SW and SE outflows. The intensification rate with SW outflow is the second highest, and the lowest intensification rates are under SE outflow. For the weak outflow strength (6–8 m s−1), intensification rates associated with the NW, SW, and NE outflows are similar and higher than that of the SE outflow. In short, the larger intensification rate is associated with the NW and NE outflow and a slower intensification rate is associated with SE outflow. Additional calculation (not shown) also shows that the NE and NW outflow environment is often associated with a much lower mean VWS than the other two outflow directions. The SE outflow is accompanied with the largest mean VWS. In addition, however, the mean VWS decreases prior to the onset of intensification in the NE and NW outflow environment, and the mean VWS for the SE outflow increases. In brief, significant outflow differences exist between the three intensification categories. For the same outflow strength, faster intensification is most favored in the NW and NE outflows and least favored in the SE outflow. As far as the authors are aware, this study is the first published work that provides the statistic relationship between direction of outflow and intensification categories.

The relative eddy momentum flux convergence (REFC) has been used as a diagnostic to identify and/or measure the upper-level TC-environment interaction [8083]. In this study, the REFC was calculated at 200 hPa over 300–600 km radially from the TC center, for all three TC IR categories. However, consistent with Shu et al. [14], no significant differences were found between the categories (not shown).

4.4. Ocean Conditions

This section studies the statistical relationship between ocean conditions and intensification rates. Note as stated at the end of Section 2.1, the OHC is based on data for a shorter 9-year period (2007–2015) due to the availability of the GOADE analysis data. Figure 13 shows the variation of SST and OHC 24 hours before and after the onset of intensification. For all three intensification categories, both SST and OHC decrease during the 48-hour period. As shown in Section 3.2, most intensification cases move from warm pool toward northwest. The ocean is generally warmer in the tropical eastern and southern areas compared with the tropical western and northern NWP basin. Figure 13 also shows that SST and OHC are the highest for RI, followed by NI and then SI. Their differences are statistically significant at least at the 90% confidence interval.

Figure 14 shows the frequency distributions of the three different intensification categories, as functions of SST and OHC. For SST, the frequencies of the three intensification categories all peak at 29.5°C (Figure 14(a)), but the frequenc peaks of RI and NI are higher than SI. Notably, the integrated frequency of RI is higher than NI and SI, and the integrated frequency of SI is the lowest when the SST is above 29°C. In contrast, the integrated frequency of RI is lower than NI and SI, and the frequency of SI is the highest when the SST is below 28°C. This result further suggests that faster intensification usually occurs over the warmer SST. In addition, Figure 14(a) shows that only a small percentage (about 5%) of RI cases occurs with SST below 28°C. This value is far higher than the threshold (26.5°C) for TC genesis and intensification [84].

For OHC, the frequency distributions of the three intensification categories are also different (Figure 14(b)). The frequency of SI is the highest followed by NI and then RI when OHC is less than 40 kJ cm−2. However, the ranking is reversed when the OHC is above 70 kJ cm−2. More than 50% of the RI cases occur over the deeper and warmer ocean where the OHC is above 70 kJ cm−2, but less than 25% SI cases occur there. The distribution of NI lies between SI and RI. This result further suggests that faster intensification usually occurs more frequently under deeper, warmer ocean conditions.

The occurrence probabilities of the three intensification categories under the same ocean heat condition are shown in Figure 15. First, the occurrence of RI is rarer than NI and SI for all SST. The occurrence probability of RI increases when the SST is above 27°C (Figure 15(a)). The occurrence probability of SI decreases with an increasing SST, whereas the occurrence probability of NI does not vary with the SST. The error bars are larger when the SST is lower. Therefore, when the ocean is relatively cold, the occurrence of the three different intensification cases is more uncertain than when the ocean is warmer. The occurrence probability as a function of OHC is in general similar to SST (Figure 15(b)). But the occurrence probability of RI increases more significantly with increase of OHC when OHC is less than 90 kJ cm−2 and decreases slightly when OHC is more than 90 kJ cm−2.

Therefore, the warmer the ocean is, the higher occurrence probability of RI is and the lower occurrence probability of SI is. To some extent, this result is consistent with Xu and Wang that shows that the maximum potential intensification rate increases with increasing SST [79]. This study however emphasizes the dependence of RI, NI, and SI on ocean conditions.

5. Conclusions

Using statistical approach, this study aims to address the following question: “Do significant differences in environmental factors, including atmospheric conditions and ocean heat conditions, exist between rapidly, normally, and slowly intensifying TCs in the North West Pacific Basin?” Based on the JTWC best track data, the intensification rates of all named TCs over the NWP for the 21-year period 1995–2005 were categorized into three intensity categories, RI, NI, and SI. The initial location and intensity, the preceding intensity change, the moving direction, the occurrence month, and the intensification duration time of TC intensification are all found to differ for RI cases compared with NI and SI cases. Most RI cases occur east of the Philippines and over the northern South China Sea. The mean initial intensity of RI cases is stronger than that of both NI and SI. The preceding intensity change is much faster for RI than for NI and SI, and RI cases tend to move WNW more than the NI and SI TCs. The mean intensification duration time of both RI and NI is found to be about 18 hr, whereas it is only about 14 hr for SI. The peak active month for RI is found to be September but is August for NI and SI. Using the ERA-Interim reanalysis data and the high-resolution global ocean analysis data, the differing impacts of the large-scale environmental atmospheric and ocean heat conditions between RI, NI, and SI were examined. Overall, there are statistically significant differences between RI and NI or RI and SI responses to deep-layer , upper-mid , and mid-low , upper-level outflow, SST, and OHC. The difference between RI and SI is significant for low-level . For net moisture inflow, the difference between RI and SI is significant only after intensification onset. There is no statistically significant difference between RI and NI in the net moisture inflow and the low-level . The differences between NI and SI significantly exist in deep-layer , upper-mid , and low-level , upper-level outflow, and SST/OHC before, during, and after the intensification but only exist in mid-low and net moisture inflow after intensification onset. These results are summarized in Table 2.

In summary, RI usually occurs under lower and decreasing VWS, enhanced and higher net moisture inflow, deeper warm ocean conditions, and stronger and increasing upper-level outflow. In contrast, SI often occurs in higher and increasing VWS, slightly increasing and lower net moisture inflow, colder ocean conditions, and weaker and static upper-level outflow. NI easily occurs in the environmental conditions that lie in between those of SI and RI. The deep-layer VWS and upper-mid layer VWS are more impactful in distinguishing the three intensification rate categories than mid-low and low-level VWS. The analysis also reveals that the direction, strength, and the increasing trend of outflow all can impact the rate of TC intensification. For example, both the NW and NE outflows favor RI, but the SE outflow is often associated with slower intensification. In addition, the differing impacts of OHC were found to be highly significant and more clearly distinguish the three different intensification categories than with only the SST.

This study is only limited to examine the statistical relationship between RI, NI, and SI and the individual large-scale environmental factors. It is possible that favorable conditions sometimes can cancel the negative impact of unfavorable conditions on TC intensification. For example, in the situation of weak vertical wind shear and the high upper-ocean heat content, TC may not intensify if other environmental conditions are not favorable. In other words, the TC intensity change largely depends on the integrated impact of all the environmental factors. The results of this study indeed verify that VWS, outflow, ocean heat condition, and moisture flux can all impact TC intensification. In other words, the result of the current study also suggests that one should consider all these environmental factors when determining and predicting the TC intensification categories. Although the focus of this paper is not to directly improve SHIPS (e.g., [64]), the results of the study suggest SHIPS may benefit from using additional predictors such as the trend of the VWS, the trend of the upper-level outflow strength, and the direction of outflow during the TC intensity forecast.

It is also found that both the strength and direction of upper-level outflow can impact TC intensification rates. Precisely how the outflow direction and how the interaction of outflow and the wind shear affect intensity intensification rates are still not clear, which will require further modeling studies. Difference of moisture conditions between the RI and NI categories is not significant. However, this does not imply that the moisture conditions are unimportant for TC intensification. To some extent, the monsoon background over NWP basin usually provides favorable moisture condition for TC intensification during an active TC season. It is also found that the occurrence probability of both RI and the mean intensification rate decreases with an increase of Quv when Quv is higher than 40 kg m−2 s−1 (not shown), suggesting that excessive moisture is not favorable for TC RI. In addition, satellite observations showed that the impact of moisture on TC intensity change depends on its quadrant distribution or relative location to the TC motion [35, 85]. This study uses statistical approach to associate the three intensification categories with the environmental factors. Additional modeling studies should be planned to understand the dynamical and thermodynamical processes that explain the relationship between the large-scale environment and the three TC intensification rates.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The research documented in this paper was supported by ONR grants N00014-14-1–0125 and N000141712111.