1. Introduction
Elsberry et al. (2020) [
1] have developed ensemble-based tropical cyclone (TC) formation guidance products based on the European Centre for Medium-range Weather Forecasts ensemble (ECEPS) and the National Centers for Environmental Prediction Global Ensemble Forecast System (GEFS). Elsberry et al. [
1] used the Marchok (2002, 2020) [
2,
3] vortex tracker outputs from the ECEPS and GEFS to provide precise (nearest 6-h time) Time-to-Formation (T2F) timings and positions along Weighted-Mean Vector Motion (WMVM) track forecasts in which the largest (smallest) weight is given to the ensemble member motion vector that is closest to (farthest from) the most recent 12-h WMVM vector. The T2F (25 kt) timing is defined to be when the weighted-mean (same weighting factors as for the track) intensities of the various ensemble members exceed the 25 kt threshold (>22.5 kt for the ECEPS). For the T2F (35 kt) timing, two Marchok vortex tracker [
2,
3] “genesis parameters” along the ensemble member tracks are utilized to calculate the weighted-mean (again, same weighting factors as for the track forecast) Warm Core Magnitude (WCM) along the WMVM track forecast. To reflect the (1/p) factor in the hypsometric equation, the Upper (600–300 mb)-layer warm core is weighted two times the Lower (900–600 mb)-layer warm core.
Elsberry et al. [
1] used examples of pre-TC circulations in the western North Pacific during the latter half of the 2019 season to calibrate the weighted-mean WCM thresholds that are associated with the T2F (35 kt) for both the ECEPS and the GEFS. Elsberry et al. [
1] provided detailed validations for Typhoons Lingling and Bualoi (and one hurricane in the eastern North Pacific and one hurricane in the Atlantic). The T2F (35 kt) timing and position errors were often smaller and more consistent than the T2F (25 kt) timing and position errors. The likely explanation for the more erratic T2F (25 kt) timing errors is the initial intensities of the pre-TC circulations in the ECEPS and (especially) GEFS were often already 25 kt.
An example of the ECEPS weighted-mean WCM evolutions for a series of forecasts from 00 UTC 27 August (082700) through 12 UTC 1 September (090112) for Typhoon (TY) Lingling (2019) is provided in
Figure 1b. The corresponding WMVM track forecasts in
Figure 1a illustrate that Lingling originated east of the Philippines and had a poleward track that is typically more difficult to predict than a westward-mover at lower latitudes or a recurver storm. According to the Joint Typhoon Warning Center (JTWC) Working Best Track (WBT), Lingling became a Tropical Storm at 06 UTC 2 September 2019 near 16.8° N, 125.1° E. The 00 UTC 29 August forecast with 23 ensemble members (inset,
Figure 1b) had the earliest detection (−4.5 days) of the pre-Lingling circulation before the T2F (35 kt), and yet the WCM threshold corresponding to that T2F (35 kt) was predicted very accurately (−0.25 days; Elsberry et al. [
1], Figure 16, column 13). Indeed, all nine of these ECEPS predictions of when the WCM evolution surpassed the threshold for the T2F (35 kt) were accurate within 0.75 days.
The actual Lingling intensity is over-plotted in
Figure 1b (kt, dash-dot line, right-axis). While the WCM = 25 threshold is for a tropical storm, the actual intensity of Lingling was only 25 kt at the time the ECEPS-predicted WCM threshold was exceeded. Whereas the ECEPS WCM= 45 threshold for a typhoon is surpassed by all nine of these forecasts, eight of the nine forecasts in
Figure 1b achieve that typhoon threshold 24–36 h after the actual time that Lingling became a typhoon according to the JTWC WBT (00 UTC 3 September). Thus, the ECEPS does predict WCM evolutions that exceeded the WCM threshold value corresponding to a Likely Storm Category (LSC) of Typhoon, the timing of when that typhoon intensity (surface winds > 64 kt) is achieved as not well-predicted by the ECEPS in this Lingling case. Since this poor intensification guidance inferred just from the WCM evolution was found with all typhoons (and the two hurricanes examined), Elsberry et al. [
1] provided intensity forecasts following the T2F from the Tsai and Elsberry (2019) [
4] combined three-stage 7-day Weighted Analog Intensity Pacific (WAIP) prediction technique (see
Section 2 for details of the WAIP). While these WAIP intensity predictions were successful in predicting the timing and magnitude of the intensification from the T2F to the Typhoon stage, the WAIP intensity would then level off and under-predict the peak intensity—especially when rapid intensification (RI; here ≥ 30 kt/24 h) occurred. Furthermore, it was difficult to estimate the timing of the peak intensity within that extended period of level intensities.
Predicting RI events is a difficult forecasting problem in all TC basins (see comprehensive reviews in Courtney et al., 2019a [
5], Courtney et al., 2019b [
6]), and many research studies have focused on RI (see summaries in Hendricks et al., 2019 [
7]; Braun et al., 2018 [
8]; and Vigh et al., 2018 [
9]). When the environmental conditions over the western North Pacific (WPAC) are favorable for TC formation, a RI event will frequently follow formation. The intensity evolutions for the WPAC cases studied by Elsberry et al. [
1], and including the TY Lekima case studied by Tsai et al. (2020) [
10], are displayed in
Figure 2. These cases form the database for this study that focuses on the importance of the T2F accuracy and of a modification of the combined WAIP to better predict the peak intensity when at least one RI event follows the T2F.
Based on the JTWC WBT, the TY Hagibus (20 W) case is an example of the extreme RI events that can occur in WPAC, and that extreme RI began only 18 h after the T2F (35 kt), which was only 36 h after JTWC WBT began for 20 W (
Figure 2). TY Halong (24 W) was the most intense among this sample of cases from the 2019 WPAC season, and the intensification from 50 kt just 18 h after T2F (35 kt) to 155 kt only 54 h later is obviously an extremely difficult intensity forecast challenge. The TY Bualoi (22 W) rapidly intensified almost linearly from 20 kt at 12 h before the T2F (35 kt) to 100 kt only 48 h after T2F (35 kt). However, other typhoons such as Lekima (10 W) [Neoguri (21 W)] did not immediately rapidly intensify after T2F (35 kt), but starting from 45 kt [50 kt] at T2F (35 kt) + 36 h [+30 h] then rapidly intensified to 130 kt [95 kt] in just 60 h [24 h]. Although Tapah (18 W) only reached 60 kt according to the JTWC WBT (
Figure 2), it is included within this sample to emphasize that an intensity forecast technique must also address slower intensification cases, and of course non-developers.
The objectives of this study are to: (i) demonstrate the ECEPS performance in early detection of pre-TC circulations in WPAC; (ii) predict the T2F (25 kt) via weighted-mean model intensity changes and the T2F (35 kt) via the predicted WCM evolution; and (iii) utilize the Likely Storm Categories based on the predicted WCM to better select the analogs in the combined WAIP technique to predict the magnitude and timing of the peak intensity of the storm within 7 days after the T2F.
Section 2 will describe the changes in methodology compared to the Elsberry et al. [
1] and the Tsai et al. [
10] studies of these same storms during the 2019 WPAC season. Examples of the ECEPS performance and some statistical results in terms of Mean Absolute Errors (MAEs) and Mean Biases (MBs) will be presented in
Section 3. A summary and some concluding remarks will be given in
Section 4. Although Elsberry et al. [
1] also described the performance of the GEFS, a new FV3 version of the GEFS was implemented in late September 2020, so it is uncertain whether that performance evaluation of those original GEFS forecasts during the 2019 season is still relevant.
2. Methodology
The Tsai and Elsberry [
4] combined three-stage 7-day WAIP technique was developed to provide the intensity evolution and the intensity uncertainty that is most likely to occur given the track forecast. The basic hypothesis is that the predominant factor in the intensity forecast beyond 72 h is the TC track, and thus in the selection of analogs more weight is given to the portion of the track beyond 72 h (Tsai and Elsberry 2014) [
11]. In each step in the development of the combined WAIP technique, the JTWC best-track was utilized rather than the official track forecast, which allowed Tsai and Elsberry (2015) [
12] to extend the WAIP technique to 7 days even though JTWC did not have official 7-day track forecasts. Thus, the evaluations of the combined WAIP technique in Tsai and Elsberry [
4] are labeled as an “optimum performance demonstration.”
The three stages in the combined WAIP intensity predictions are: (i) Pre-formation stage from the first detection of a pre-TC circulation in the ECEPS until the Time-to-Formation (T2F), where formation may be defined as either Tropical Depression (25 kt) or Tropical Storm (35 kt); (ii) Intensification stage immediately following T2F, in which 16 analog storm intensities are weighted according to the same weighting factors that were assigned to the analog tracks to generate the weighted-mean forecast intensities each 12 h (rather than a simple mean intensity that assumes each analog intensity is equally likely); (iii) Ending storm stage, which may be due to landfall, extratropical transition, or to non-development within the 7-day forecast interval. A highly accurate ECEPS track forecast is obviously very important for determining landfall timing, and also for anticipating intensity changes associated with extratropical transition. Each of these three stages is defined in terms of the weighted-mean warm core magnitude (WCM) evolution along the WMVM forecast track of the ensemble storm predicted by the ECEPS starting from the first detection of a pre-TC circulation.
The preparatory steps for the modified combined WAIP intensity predictions in this study are summarized at the top of
Figure 3. As indicated in
Section 1, the ECEPS 51-member ensemble WMVM track forecasts are the basic input to the combined WAIP intensity prediction technique, so that the potential real-time performance can be demonstrated. Three important benefits from these ECMWF ensemble forecasts are that they often detect a pre-TC circulation 2–3 days before the JTWC typically begins issuing storm warnings, the 7-day track forecasts after TC formation are usually highly accurate, and the up to 15 day total length of the tracks provides a longer-range outlook as to what regions/countries may be affected.
Elsberry et al. [
1] describe the methodology for calculating the weighted-mean WCM from the two Marchok genesis parameters, which originally are from the Hart (2003) [
13] Cyclone Phase Space variables along the WMVM track forecast, and then how the T2F (35 kt) timing was estimated when the WCM exceeded a threshold value as illustrated in
Figure 1b for pre-TY Lingling. While Elsberry et al. [
1] attempted to distinguish an early and a later T2F (35 kt) timing with two WCM threshold values, these two threshold values were typically surpassed at the same time or within six hours for the rapidly intensifying storms. Consequently, only the WCM threshold for the early T2F (35 kt) described in Elsberry et al. [
1] will be utilized in this study.
Elsberry et al. [
1] also describe how the maximum positive weighted-mean WCM for the ECEPS forecasts could be used to define the Likely Storm Category (LSC): Maximum WCM < 25.0 for Tropical Depression; 25.0 ≤ WCM < 46.0 for Tropical Storm; and WCM ≥ 46.0 for TY. Note that this WCM value for a TY is also displayed on the left ordinate of
Figure 1b, and eight of the ECEPS forecasts surpassed that maximum WCM threshold for a Likely TY. Whereas six of these ECEPS forecasts had maximum positive WCM values just beyond the 46.0 threshold for a TY to magnitudes ~70, the ECEPS forecast from 00 UTC 29 August (082900) with a maximum WCM ~140 and the forecast from 12 UTC 30 August with a maximum WCM ~110 far exceed the threshold value such that the JTWC forecaster could have confidence that the pre-Lingling circulations in those forecasts would become a typhoon.
Whereas Elsberry et al. [
1] had calibrated the WCM threshold for the T2F (35 kt) timing along the WMVM track, they found the likely timing of when TY intensity would be achieved could not be inferred from just the predicted WCM evolution. Consequently, Elsberry et al. [
1] calculated the intensity forecasts after the T2F from the Tsai and Elsberry [
4] combined three-stage 7-day WAIP. These variables are then the preparatory steps for the modified combined WAIP intensity forecast (
Figure 3, top row).
The pre-formation stage for the modified WAIP utilized in this study is not changed (
Figure 3, middle). The bifurcation version of the Tsai and Elsberry [
4, Appendix] combined WAIP is utilized for the intensification stage after formation (
Figure 3, lower flowchart). If the weighted-mean intensity spread over the 7-day forecast period exceeds a specified threshold value, a hierarchical cluster analysis (Wilks 2011) [
14] separates the 16 analog intensities into two clusters. If there are at least three analogs in each cluster, the WAIP technique is separately applied to these two clusters to produce weighted-mean intensities (and weighted-mean intensity spreads) each 12 h over the 7-day forecast period. The intensity Cluster-1 is always that cluster with the larger peak intensities, and Cluster-2 is the alternate solution with the lower peak intensity. These two cluster WAIP forecasts are considered to represent a substantial intensity bifurcation if a 15 kt intensity difference threshold is exceeded for at least 25% of the 12-h forecast intervals within the 7-day forecast period.
The new procedure in this modified WAIP is to now check the LSC based on the maximum positive WCM as described above. If there is a strong LSC indication of an early and long-lasting WCM increase indicative of a likely TY (e.g., almost all of the storms in
Figure 1b), then the bifurcation WAIP Cluster 1 is selected. As indicated in
Section 1, a simple mean of the analog peak intensities that occur at different times results in an extended time period of level intensities. Therefore, a special weighting of the analog intensities with the largest weight given to the analog with the largest peak intensity is added in order to increase the Cluster-1 peak intensity and to define a specific timing of the peak intensity (see
Appendix). This special weighting of the peak intensities results in a larger peak WAIP intensity at a better-defined peak intensity time estimate compared to an average of all Cluster-1 analog intensities.
If the WCM evolution indicates a more moderate intensification rate, the original 16-analog WAIP intensification scale is considered to be a more appropriate choice when there is no clear indication of a rapid intensification based on the predicted WCM evolution (
Figure 3, middle). If there is only a weak (or no) indication of a substantial intensification following the T2F (which may be for 25 kt), then the bifurcation WAIP Cluster-2 is selected, and the weighted peak intensity and peak intensity time procedure is applied. In the latter case, the objective is to identify in the ECEPS forecasts those situations in which the environmental constraints (e.g., vertical wind shear and sea-surface temperature) are likely to inhibit intensification beyond the T2F (35 kt).
In some ECEPS forecasts there may be no bifurcation indicated by the analog intensity spread according to the rules in the Appendix of Tsai and Elsberry [
4], but the maximum positive WCM value clearly indicates a likely TY (
Figure 3, bottom branch of the flow chart). This situation may occur over the most favorable TC formation areas of the WPAC (e.g., east of the Philippines) when nearly all of the selected analogs have similar large intensification rates such that the potential Cluster-1 and potential Cluster-2 do not meet the required 15 kt intensity difference threshold for at least 25% of the 12-h forecast intervals. Rather than averaging a large number of similar analogs with peak intensities at different times, which again would lead to an extended period with level intensities, a weighted-mean peak intensity and peak intensity time are defined from only the top three analog peak intensities. Finally, if no bifurcation exists and the maximum positive WCM value does not indicate a likely TY, the original WAIP that is the weighted-mean intensities of the 16 analogs becomes the final combined WAIP intensity forecast (
Figure 3, bottom decision in the flow chart).
In summary of this extension of the Elsberry et al. [
1] article, the ECEPS forecasts during the 2019 WPAC season from the first detections will be examined with the objective of obtaining the earliest possible predictions that the pre-TC circulation will have a precise formation time (T2F) and will likely have a RI such that it is highly likely a typhoon will subsequently occur. As indicated in
Figure 1b, Elsberry et al. [
1] had calibrated the T2F (35 kt) timing in terms of a WCM threshold representative of a Tropical Storm, and had defined a second WCM threshold that would indicate a Likely Storm Category of Typhoon would occur at some time during the life cycle. Whereas Elsberry et al. [
1] had applied the Tsai and Elsberry [
4] combined WAIP technique to predict the intensification after the T2F (35), they had included all 16 analog intensities rather than selecting only analog intensities that included an RI when a typhoon was being predicted. Consequently, the WAIP tended to have an extended period of level intensities rather than a well-defined peak intensity since the different timings of the peak intensities were not accounted for. The modified combined WAIP in
Figure 3 that will be applied after the T2F (35 kt) in this study utilizes the bifurcation version and then also uses three categories of the maximum WCM-based Likely Storm Categories to be highly selective as to which analog intensities should be used. Selecting only bifurcation WAIP Cluster-1 intensities that have RI events leading to the largest peak intensities will be utilized when the WCM-based Likely Storm Category is TY. Another important modification will be to implement a weighted peak intensity technique that takes into account the different times (and magnitudes) of those more appropriately selected analog intensities for the bifurcation WAIP cluster intensity predictions, with the objective of predicting larger peak WAIP intensities associated with RI events at a better-defined time of peak intensity. A description of this weighted peak intensity procedure is given in first section of the
Appendix A.
3. Examples of the Modified Combined WAIP Intensity Predictions
A detailed description of the modified combined WAIP intensity predictions based on both the ECEPS T2F (25 kt) and the T2F (35 kt) timing forecasts for pre-TY Lingling (
Figure 4) will first be given, and then for the other examples only the WAIP predictions based on the ECEPS T2F (35 kt) forecasts will be described. These WAIP intensity predictions will be compared with the JTWC WBT intensities as the post-season best-track files were not available. However, when the ECEPS has predicted a pre-TC circulation existed before the JTWC had begun a WBT, it was necessary to extrapolate the WBT intensities back in time to the beginning of the pre-TC circulation in the ECEPS, which is assumed to have an intensity V
max ≥ 15 kt. If the WBT initial intensity (V
0) is already 15 kt, the intensity records are simply backward-extended as 15 kt to the beginning of the pre-TC circulation (
Figure 2, lower left). If V
0 is 20 kt, the timing of 15 kt is estimated by using the decreasing trend between the 20 kt timing and 25 kt timing of that WBT. If the initial V
0 in the WBT is 25 kt, the timing of the 15 kt is backward-extended based on the trend between the 25 kt timing and the 35 kt. As indicated by the blue-circled crosses in
Figure 2, these backward-extended intensities/timings provide reasonable values for validation of the predicted WAIP intensities during the pre-formation stage.
This modified WAIP intensity prediction in
Figure 4a based on the ECEPS T2F (25 kt) has been included for pre-TY Lingling from the ECEPS because seven of the nine predictions are outstanding in terms of the sustained rapid intensification and peak intensity predictions. As indicated in
Figure 1a, these pre-Lingling circulation track forecasts began just east of the Philippines (note that the dots indicating the T2F are for 35 kt rather than 25 kt), and then had a long track poleward over warm ocean water, which would be favorable for a sustained RI event. Two-thirds of these ECEPS forecasts also predicted well the timing of the rapid decay following the peak intensity. Consequently, the intensity Mean Absolute Errors (MAEs;
Figure 4c) are incredibly small for the 7-day WAIP intensity forecast period following the T2F (25 kt). While the inclusion of the first ECEPS forecast (inset, 082700) and last ECEPS forecast (090112) no doubt contributed to an intensity MAE = 20 kt at tau = 126 h, most of the other forecast intervals have MAEs < 15 kt. The intensity Mean Bias (MB;
Figure 4e) is plotted relative to the actual T2F (25 kt, tau = 0 h) to illustrate that the pre-T2F period of the combined WAIP has very small errors, which may be attributed to most of the ECEPS forecasts not predicting the pre-Lingling circulation very early prior to the T2F (25 kt). As is also evident in
Figure 4a, the last (090112) ECEPS forecast prior to the T2F (25 kt) had a large under-forecast bias, and the first (082700) forecast had an increasingly larger under-forecast bias from tau = +48 h to tau = +90 h. These two under-forecast bias cases distorted the MB at the end of the forecast period because the remaining ECEPS forecasts of the T2F (25 kt) had only small positive biases until tau = +72 h.
For the pre-Lingling case, the modified WAIP intensity predictions in
Figure 4b based on the ECEPS T2F (35 kt) and WCM evolution are not as accurate as the T2F (25 kt)-based WAIP intensity predictions in
Figure 4a. The pre-formation period WAIP intensity predictions in
Figure 4b are slightly better because the T2F (35 kt) is better predicted, except for the first (082700) ECEPS forecast. Consequently, the MAEs (
Figure 4d) during the first 72 h and the MB during the pre-formation period are slightly smaller. After the T2F (35 kt), almost all of the modified WAIP intensification stage predictions were about one day late to typhoon stage and missed the extremely rapid intensification on 4 September.
The two exceptions were the 082812 ECEPS forecast, which detected the pre-Lingling circulation of the 082906 [note that this is four days before T2F (35 kt)] and the 083012 ECEPS forecast, which both predicted well the time-to-typhoon and then achieved a peak intensity of ~105 kt (albeit more than one day late). The first exception in predicting the RI of Lingling was from the first ECEPS forecast that had a very early detection of the pre-Lingling circulation and also had an accurate T2F prediction, but the ECEPS then under-predicted the vortex spin-up. It was not unusual for the first detection of a pre-TC circulation in the ECEPS to have an under-prediction of the intensification, and especially as in this case when that intensification is occurring ~5 days after the first detection. We can only speculate that the more-accurate next ECEPS forecast 12 h later had the benefit from initial conditions that had a better definition of the pre-Lingling circulation. [Note: ECMWF does not bogus the vortex based on warning center TC Vitals.] The second exception in predicting the RI of Lingling was from the ECEPS forecast that was initiated just 6 h before Lingling became a Tropical Storm according to the JTWC. Without any history of the earlier stage of the pre-Lingling vortex, and again ECMWF does not bogus the vortex, the T2F WAIP diagram in
Figure 4a indicates the Lingling vortex was only slowly spun-up, which means the WCM only slowly increased in the model. Consequently, the WCM threshold for a Typhoon was not met (i.e., only Cluster 1 analog intensities were not selected), and therefore all 16 analog intensities were selected for the T2F WAIP (
Figure 4b). With those 16 analog intensities, and a decay stage following the first peak intensity that was not well-defined, our peak intensity modification was not applied.
Note the five WAIP intensity forecasts (
Figure 4b) with peak intensities between 65 kt and 75 kt, and that the timings of those peak intensities are not well-defined as in
Figure 4a. For these five WAIP forecasts the weighted peak intensity modification technique was not applied because the bifurcation WAIP Cluster 1 versus Cluster 2 intensity difference (15 kt over at least 25 percent of the 12-hourly intensities) was not satisfied. That is, too few Cluster 1 peak intensity analogs were identified, presumably because of the favorable formation conditions that exist north and east of the Philippines where the pre-Lingling circulation began (first set of dots in
Figure 1a). Consequently, the MAEs (
Figure 4d) increase to 30 kt at the time of peak intensity of Lingling, but then decrease to near 10 kt during the decay stage because all but two decayers of the ECEPS T2F (35 kt)-based WAIP intensity predictions are relatively accurate during the decay stage. The MAEs are inversely related to the MB errors, which are negative (under-forecast) and grow slowly during the first 48 h after the T2F (35 kt), but then become large negative during the extreme RI period of TY Lingling.
While the ECEPS T2F (35 kt)-based WAIP intensity predictions are not as accurate as the T2F (25 kt)-based predictions, they do include highly accurate formation (35 kt) timing forecasts and all but the first (082700) ECEPS-based WAIP prediction indicate that the pre-Lingling circulations will intensify to at least typhoon intensity within 2–3 days after that T2F (35 kt). In conjunction with the ECEPS WMVM track forecast (
Figure 1a) that now also provides the positions at which formation (35 kt) and the TY intensity will likely occur, these WAIP intensity predictions would be helpful guidance for many TC disaster preparation activities even if the actual peak intensity of 120 kt is not predicted.
5. Summary and Discussion
When the environmental conditions over the western North Pacific are favorable for tropical cyclone formation, a RI event will frequently follow formation. Thus, an accurate prediction of that RI event will require an accurate prediction of formation. Elsberry et al. [
1] have used the Marchok [
2,
3] vortex tracker outputs from the ECEPS to provide precise (nearest 6-h synoptic time) T2F timings and positions along WMVM track forecasts. The T2F (25 kt) timing, which is simply the weighted-mean of the times that the various ECEPS ensemble storm members exceed 22.5 kt, is occasionally quite accurate, and when that accurate T2F (25 kt) leads to analog selections in the WAIP that had RI events following their formation, then an accurate WAIP intensification stage was predicted. Even though the ECMWF often predicts the existence and track of a pre-TC circulation 3–4 days before the T2F (25 kt), the predicted intensities of those circulations tend to be irregular with a T2F (25 kt) that is biased early, which leads to too-early predictions of RI.
For the T2F (35 kt) timing, Elsberry et al. [
1] utilized two Marchok vortex tracker “genesis parameters” to calculate the WCM along the WMVM track forecast, and defined weighted-mean WCM thresholds that defined the T2F (35 kt) for pre-TC circulations in the western North Pacific during the latter half of the 2019 typhoon season. Similarly, Elsberry et al. [
1] defined a weighted-mean WCM threshold that for the ECEPS was indicative of a Typhoon intensity. Whereas the ECEPS does frequently predict WCM evolutions that exceeded the WCM threshold value corresponding to a Likely Storm Catetgory (LSC) of Typhoon, the timing of when that typhoon intensity is achieved tends to be biased late. Therefore, Elsberry et al. [
1] also provided 7-day WAIP intensity predictions following T2F (35 kt) that generally improved the timing of when the typhoon stage would be achieved. However, the WAIP intensity evolutions tended to level-off after reaching typhoon intensity such that the peak intensity timing was not well-defined and the magnitude was under-predicted (e.g.,
Figure A1a).
A summary of the ECEPS-predicted T2F and modified WAIP intensity predictions for the western North Pacific TCs during the late 2019 typhoon season as displayed in
Figure 2 is given in
Table 1. As the first objective of the T2F-WAIP technique is to provide earlier guidance as to the intensity during the pre-formation stage of WAIP, the first column in
Table 1 is the time before formation that the pre-TC circulation is detected in the ECEPS forecasts. In all but one TC, the circulation was detected at least 66 h prior to T2F (35 kt), and in the pre-TY Lekima case the first detection was 132 h in advance (
Figure 5). The exception is for pre-TY Hagibis, which was detected in only three ECEPS forecasts that originated from only—6 h to—36 h in advance of T2F (35 kt). Such a short detection time is a likely indicator that the ECEPS did not have a good handle on this pre-TC circulation, and that is reflected in the “very poor” in
Table 1, columns 3 and 4, respectively, for pre-TY Hagibis.
Tsai and Elsberry [
4] had demonstrated the potential benefit for highly accurate WAIP pre-formation stage intensity predictions, IF an accurate T2F timing is provided. In the Tsai and Elsberry [
4] “optimum performance” test, that timing was from the best-track file. In
Table 1, column 2, the T2F timing accuracy of these ECEPS forecasts during the late 2019 typhoon season is evaluated. The most common entry is “Good-slightly early,” which means within 12 h of the actual T2F. Two entries are “Excellent,” which means perfect timing or within 6 h, but in the Hagibis case those two excellent predictions were for very short detection times (
Table 1, column 1). The “Very delayed” T2F entry for pre-Bualoi is noteworthy because the WAIP intensification stage (
Table 1, column 3) and the peak intensity (
Table 1, column 4) were also very delayed, and consequently the intensity MAEs and MBs were quite large (
Figure 7b).
The first objective of the modified WAIP prediction in the flowchart in
Figure 3 was to improve the bifurcation WAIP intensification rate following the T2F (35 kt) by being more selective when the WCM evolution indicated a Likely Storm Category of Typhoon. Specifically, by selecting only the Cluster 1 analogs with the larger peak intensities in these situations, the WAIP intensification rate will be based on RI events following T2F (35 kt). However, the first “Excellent” entry for pre-TY Lingling in
Table 1, column 3, is actually for T2F (25 kt), because it was used in
Figure A1 to illustrate the modified peak intensity timing and magnitude. A majority of the storms in
Table 1, column 3, have at least a “Good” WAIP intensification rate, but some also have qualifying notations such as “mostly slow” or “but early”. As mentioned above for pre-Bualoi, the “very delayed” T2F (35 kt) leads to a very delayed intensification stage. The noteworthy Hagibis and Halong cases with “Very poor” performance ratings are examples in which extreme RI occurred, but the modified WAIP analog selection condition limiting those analogs to RI event cases was not appliable. Consequently, large intensity MAEs and large negative MBs are found for these cases in
Figure 8b and
Figure A3b, respectively.
The second objective of the modified WAIP predictions in the flowchart in
Figure 3 has been to improve the peak intensity magnitude and timing by fitting a cubic spline curve starting from a point on the intensification side through the weighted-mean of the peak intensities and then to a point on the decay side (
Figure A1b). In addition to that “Excellent” example of a modified peak intensity magnitude and timing for the first pre-Lingling case (
Table 1, column 4), the other “Excellent” result was for the pre-Tapah case in
Figure 6a that had an “Excellent” WAIP intensification stage prediction all the way to peak intensity. For the pre-Bualoi case (
Figure 7a), the modified peak intensity magnitudes were good, but they were delayed by two days due to the delayed T2F (35 kt). In four other cases in which the peak intensity was poor, the modified peak intensity conditions were not applied because the intensification stage or the decay stage did not meet the present criteria for a cubic spline fit. Consequently, further testing is required to improve the criteria for modification, although it is not expected that a statistical analog technique such as the modified WAIP will be able to match the extreme RIs that occurred in Hagibis (
Figure 8a) or Halong (
Figure A3a). Nevertheless, the ECEPS WCM evolutions and the modified WAIP predictions are successfully indicating that these pre-TC circulations detected as much as 3–4 days before T2F (35 kt) are going to be at least Typhoon intensity. Even if the very high peak intensities of Super Typhoons like Halong are not being matched, guidance that a typhoon threat exists along a 15-day ECEPS WMVM track could be used for improved risk management decisions.
One of the requirements for this T2F-WAIP guidance product is that it be objective as the forecasters do not have time to examine individual ensemble member forecasts among the 51 ECEPS members. Thus, the T2F (25 kt) and T2F (35 kt) timing (nearest 6-h synoptic time) and locations along the ECEPS and GEFS during the 2019 season WMVM track forecast are provided, and those variables are the inputs to the combined three-stage, 7-day WAIP intensity predictions following the T2F. Because of the larger track uncertainty in the TC pre-formation stage than in the mature stage, there is considerable uncertainty in the predicted variables that go into the T2F and thus into the WAIP predictions. Thus, these T2F-WAIP predictions must be considered by the forecasters in conjunction with other formation guidance (e.g., climatology, Madden-Julian Oscillation, deterministic global and regional models, etc.). Because T2Fs and 7-day WAIP intensity predictions are provided along the 15-day ECEPS WMVM track forecasts, the forecaster is able to consider the regions that might (or might not) be at risk of a typhoon. A good example is the pre-TY Lekima formation predictions by the ECEPS described by Tsai et al. [
10], in which there was an opportunity to provide early (10 days in advance) and consistent warnings of the threat to the Zhejiang, China area.
As an example of the critical need for forecaster interpretation of this T2F-WAIP guidance, the ECEPS forecasts of TS Peipah (2019) are presented in
Figure 9. As described by Elsberry et al. ([
1], Section 3.1.2), the first ECEPS track forecast from 00 UTC 13 September for pre-Peipah had a maximum of 22 of a possible 51 ensemble members, and although there was considerable scatter the WMVM track was northwestward toward Japan (
Figure 9a). This early detection of pre-Peipah was 4.75 days before the T2F (25 kt) of 18 UTC 14 September, and the ECEPS perfectly predicted this T2F (25 kt) (
Figure 9c). The WAIP intensity prediction based on that T2F (25 kt) timing and WMVM track was a relatively slow increase to a peak intensity of 80 kt after 120 h (
Figure 9c). Note in
Figure 9b that the WCM was near zero until 12 UTC 14 September and only slowly increased to the WCM = 25 threshold for a Tropical Storm at 00 UTC 17 September.
It is noteworthy that with the T2F (35 kt) estimate and the WMVM track forecast in
Figure 9a, the WAIP intensity prediction was that the Peipah circulation would intensify to 75 kt at 96 h and maintain that intensity through 168 h after the T2F (35 kt) (
Figure 9d). Given that this ECEPS forecast with a maximum of 22 ensemble members is predicting a track toward Japan as well as crossing a major maritime shipping lane, and with the WAIP intensity predictions of a typhoon based on both the T2F (25 kt) and the T2F (35 kt), presumably the forecaster should issue at least an alert of a Typhoon threat.
What complicates the situation is there was then a 3.5 day gap between the first ECEPS forecast and the next ECEPS forecast from 12 UTC 13 September, which had 45 of a possible 51 members and had pre-Peipah beginning at 13.2° N, 155.7° E, and later had a track brushing the south coast of Japan (
Figure 9a). The next ECEPS forecast 12 h later starting from 00 UTC 14 September had an even longer track passing along the southern coast of Japan so the circulation would be moving over the Kuroshio ocean current. Both of these ECEPS forecasts had T2F (25 kt) times that were slightly early, and then the modified WAIP intensity predictions were ~65 kt and ~110 kt, respectively (
Figure 9c). By contrast, neither of these forecasts had WCMs that surpassed the WCM = 25 threshold indicating a T2F (35 kt) time (
Figure 9b). Consequently, no WAIP intensity predictions were triggered for these two ECEPS forecasts, which would be a conundrum for the forecaster when the T2F (25 kt) had led to such strong WAIP intensity forecasts.
Elsberry et al. [
1], their Figure 7c also provide the GEFS forecasts from 12 UTC 13 September, and not only was the WCM threshold for a TS not met, the T2F (25 kt) was not met until 00 UTC 20 September, which was too late to initiate a 7-day WAIP intensity forecast. Because Elsberry et al. (2020) had repeatedly documented that the model physical processes in the GEFS were more active than in the ECEPS, the fact that the GEFS forecasts did not indicate a Tropical Storm would have given the forecaster a basis to not re-initiate an alert after the 3.5 day gap, even though the first ECEPS had seemed to strongly indicate the threat of a typhoon.
Finally, the last two ECEPS from 12 UTC 14 September and 00 UTC 15 September had accurate track forecasts (
Figure 9a) of a recurvature rather than a landfall, neither had weighted-mean intensity exceeding 22.5 kt so there were no WAIP intensity forecasts in
Figure 9c for those forecast times. While both of these forecasts had initial WCM values above the threshold for a Tropical Storm (
Figure 9b), those WCM values then quickly decreased to zero, which indicated extratropical transition so no WAIP forecasts were generated (
Figure 9d).
The purpose of this discussion of pre-Peipah forecasts has been to emphasize the need for forecaster involvement for interpreting these T2F-WAIP predictions in conjunction with other guidance. Secondly, this case indicates that FV3 GEFS-based forecasts should also be tested for their performance so that an alternative to the ECEPS is available for confirmation. We are confident that the forecasters will quickly learn the characteristics of each version and how to access the level of threat of a typhoon, and now be able to issue alerts or warnings earlier with WAIP analogs that include (or do not include) the second modification to optimize the RI events when the Likely Storm Category is for (is not for) a Typhoon, and with the WAIP peak intensity timing and magnitude, the forecasters will have an opportunity to provide better typhoon risk management guidance.