Development of Statistical Typhoon Intensity Prediction: Application to Satellite-Observed Surface Evaporation and Rain Rate (STIPER)

Transcription

1 240 W E A H E R A N D F O R E C A S I N G VOLUME 27 Development of Statistical yphoon Intensity Prediction: Application to Satellite-Observed Surface Evaporation and Rain Rate (SIPER) SI GAO AND LONG S. CHIU Department of Atmospheric, Oceanic and Earth Sciences, College of Science, George Mason University, Fairfax, Virginia (Manuscript received 8 March 2011, in final form 10 September 2011) ABSRAC A statistical dynamical model has been used for operational guidance for tropical cyclone (C) intensity prediction. In this study, several multiple linear regression models and neural network (NN) models are developed for the intensity prediction of western North Pacific Cs at 24-, 48-, and 72-h intervals. he multiple linear regression models include a model of climatology and persistence (CLIPER), a model based on the Statistical yphoon Intensity Prediction System (SIPS), which serves as the base regression model (BASE), and a model of SIPS with additional satellite estimates of surface evaporation (SLHF) and innercore rain rate (IRR, SIPER model). A revised equation for the C maximum potential intensity is derived using ropical Rainfall Measuring Mission Microwave Imager optimally interpolated sea surface temperature data, which have higher temporal and spatial resolutions. Analyses of the resulting models show the marginal improvement of SIPER over BASE. However, IRR and SLHF are found to be significant predictors in the predictor pool. Neural network models using the same predictors as SIPER show reductions of the mean absolute errors of 7%, 11%, and 16% relative to SIPER for 24-, 48-, and 72-h forecasts, respectively. he largest improvement is found for the intensity forecasts of the rapidly intensifying and rapidly decaying Cs. 1. Introduction ropical cyclones (Cs) cause huge economic losses and human casualties worldwide. Hurricane Katrina in 2005 and the associated flood resulted in over 1800 lives being lost, and damage estimates were well in excess of $100 billion in the aftermath of the actual hurricane and the subsequent flooding. Cyclone Nargis in 2008 in Myanmar claimed at least fatalities and its estimated death toll was more than Damage was estimated at over $10 billion, making Nargis the most damaging cyclone ever recorded in the North Indian Ocean Basin. Accurate forecasting of C tracks and intensities is therefore essential to minimizing economic losses and human casualties caused by Cs. An evaluation of the National Hurricane Center and Joint yphoon Warning Center (JWC) operational tropical cyclone intensity forecasts for the three major Northern Hemisphere tropical cyclone basins (Atlantic, eastern North Pacific, and western North Pacific) for the past two Corresponding author address: Long S. Chiu, Dept. of Atmospheric, Oceanic and Earth Sciences, College of Science, George Mason University, 4400 University Dr., Fairfax, VA lchiu@gmu.edu decades only shows marginal improvement in intensity forecasting while the track forecasts have seen steady improvement (DeMaria et al. 2007), indicating that there are challenges in C intensity forecast remaining. he Statistical yphoon Intensity Prediction Scheme (SIPS) for the western North Pacific developed by Knaff et al. (2005) has been implemented into operations at the Joint yphoon Warning Center (JWC) since 2003, and it has produced skillful forecasts through 4 days. SIPS contains static predictors related to climatology and persistence and time-dependent predictors related to environmental and sea surface temperature (SS) conditions. With the advent of remote sensing technology, a number of remote sensing products are now available in real time or near real time. he inclusion of satellite information, such as brightness temperatures from infrared imagery and ocean heat contents (OHCs) from satellitebased altimetry (DeMaria et al. 2005; Mainelli et al. 2008), as well as predictors relating to the inner-core (within 110 km of the storm s center) precipitation and convective characteristics of Cs derived from passive microwave imagery (Jones et al. 2006), have shown improvements of up to 8% in the Statistical Hurricane Intensity Prediction Scheme (SHIPS) for the eastern North Pacific and Atlantic (DeMaria and Kaplan 1994a; DOI: /WAF-D Ó 2012 American Meteorological Society

2 FEBRUARY 2012 G A O A N D C H I U 241 ABLE 1. Datasets used in this study, as well as the temporal and spatial coverages and resolutions. Dataset Spatial resolution emporal resolution Spatial coverage emporal coverage RSMC okyo best track hourly Western North Pacific and the South China Sea NCEP GFS analysis 18 6 hourly Global 2000 present MPA rain rate hourly 408N 408S 1998 present OAFlux SLHF 18 Daily Global MI OI SS Daily 408N 408S 1998 present DeMaria and Kaplan 1999). However, the value of satellitebased information has yet to be fully exploited for the operational prediction of C intensity. Latent heat transfer at the air sea interface, which can be quantified by SLHF, and latent heat release occurring within the inner-core region in the atmosphere, which could be indicated by IRR, are two major heat sources for C intensification. Gao and Chiu (2010) used satellitebased SLHF and rain-rate data to show that high initial SLHFs and IRRs are usually associated with C rapid intensification in 24 h over the western North Pacific, suggesting that SLHF and IRR have the potential to be useful new predictors for C intensity forecasting. he goal of this study is to examine the impacts of SLHF and IRR on the intensity prediction of Cs over the western North Pacific. In addition to multiple linear regression models, artificial neural network (NN) models are also constructed for intensity prediction. NN models have been shown to work well in C intensity forecasts (Baik and Paek 2000; Jin et al. 2008). Section 2 describes the datasets. he predictors and the development and evaluation of multiple linear regression models are presented in section 3. NN models are developed and assessed in section 4. Finally, a summary and discussion are given in section Datasets he intensity and location information for each C is taken from best-track data produced by the Japan Meteorological Agency s Regional Specialized Meteorological Center okyo (RSMC okyo). his postanalysis best-track dataset contains 6-h the location, minimum central pressure, and 10-min maximum sustained wind speed (MWS) of all of the Cs over the western North Pacific including the South China Sea. Environmental data are derived from the National Centers for Environmental Prediction s (NCEP) Global Forecasting System (GFS) Final (FNL) gridded analysis (Yang et al. 2006) at and 6-h resolution. he environmental data include the wind, air temperature, and relative humidity at 200, 250, 300, 350, 400, 450, 500, 700, 750, 800, and 850 hpa. he divergence, relative eddy flux convergence, and relative vorticity at each grid are calculated using wind field information and the central difference method. All the environmental predictors are derived by averaging corresponding data within some specific radius; detail calculations of these parameters will be presented in the next section. he rainfall data are taken from the ropical Rainfall Measuring Mission (RMM) Multisatellite Precipitation Analysis (MPA; Huffman et al. 2007). MPA provides combined precipitation estimates from multiple satellites, as well as gauge analyses where possible. he spatial and temporal resolutions are and 3 hourly, respectively. Daily gridded SLHF data at resolution are obtained from the third release of the Objectively Analyzed Air Sea Fluxes dataset (OAFlux; Yu et al. 2008). OAFlux is an objectively blended dataset; data sources are from satellite and numerical weather prediction (NWP) model outputs while in situ observations are used to assign the weights. RMM Microwave Imager (MI) optimally interpolated (OI) daily SSs (Wentz et al. 2000) at resolution produced by Remote Sensing Systems, together with best-track data, are utilized to estimate the maximum potential intensity (MPI). able 1 shows the availability of the datasets, including their temporal and spatial coverages and resolutions. All the datasets are collected over the period from 2000 to his period represents an intersection of all available data used in this study. RSMC okyo best-track data, NCEP GFS FNL environmental data, and MPA rainrate data are collected twice daily at 0000 and 1200 UC. Only overwater C samples are considered, land effects on C intensity change are not considered in this work, and the OAFlux SLHF and MI SS data are available over the ocean only. Except for SS and SLHF, the datasets are available for all storms at all forecast times. At 24-, 48-, and 72-h forecast times, SS and SLHF are concurrently available for 75%, 68%, and 62% of all the cases, respectively. 3. Linear regression models Multiple linear regressionmodelsareusedtopredict how the overwater 24-, 48-, and 72-h intensities

3 242 W E A H E R A N D F O R E C A S I N G VOLUME 27 ABLE 2. Potential climatological, environmental, and satellitebased predictors. he predictors that are evaluated at the beginning of the forecast period are static (S), and the predictors that are averaged along a storm track from the initial time to the forecast time are time dependent (). Predictor Description Static (S) or time dependent () Climatology and persistence MWS0 Initial maximum sustained wind S speed (MWS) DMWS Change in MWS during the past 12 h S JDAY Absolute value of (Julian day 2 248) S SPD Storm translational speed S LA Latitude of storm center S LON Longitude of storm center S Environmental PO Maximum potential intensity based on Eq. (1) minus initial maximum wind speed RHLO Area-averaged ( km) RH at hpa RHHI Area-averaged ( km) RH at hpa U200 Area-averaged ( km) zonal wind at 200 hpa 200 Area-averaged ( km) temperature at 200 hpa d200 Area-averaged ( km) divergence at 200 hpa REFC Relative eddy flux convergence within 600 km at 200 hpa SHR Area-averaged ( km) hPa wind shear USHR Area-averaged ( km) hPa zonal wind shear z850 Area-averaged ( km) 850-hPa relative vorticity Satellite based SLHF Area-averaged ( box) S surface latent heat flux IRR Area-averaged (0 100 km) inner-core rain rate S (i.e., MWS) change (predictand, DELV) from the initial forecast time. a. Model formulation able 2 summarizes the predictors (including 16 original and 2 new satellite-based predictors: SLHF and IRR) used in this study. he computation of the traditional climatological and environmental predictors follows the approach in Knaff et al. (2005). All of the environmental predictors are obtained using a perfect prog approach (Kalnay 2003). Both the NCEP GFS FNL analysis and the actual C best track (from RSMC okyo) are used to develop the models. he predictors that are evaluated at the beginning of the forecast period are static, such as those predictors related to climatology and persistence; predictors that are averaged along the track of the storm from the initial observation to the forecast time are time dependent, providing the mean conditions for the storm, such as those predictors related to SS, moisture, and wind fields. he MPI is defined as the upper bound of the C intensity for a set of given atmospheric and oceanic thermal conditions (Camp and Montgomery 2001). MPI can be estimated theoretically (e.g., Miller 1958; Emanuel 1988; Holland 1997) or empirically (e.g., Merrill 1987; DeMaria and Kaplan 1994b; Whitney and Hobgood 1997; Knaff et al. 2005; Zeng et al. 2007). heoretically derived MPIs are often represented by minimum central pressures. he relationship between C maximum wind and minimum pressure relationship is complex and is determined by factors such as latitude, C size, environment pressure, and intensification rate (Knaff and Zehr 2007). o be consistent with the RMM-derived rainfall, we use the SS derived from MI and hence the empirical approach of determining the MPI as an exponential function of SS is used. he interpolated monthly SS climatology (Levitus 1982) or weekly SS analysis (Reynolds et al. 2002) used in previous studies is replaced by the high-resolution daily SS data retrieved from MI in this study. Daily SSs are expected to provide more precise and up-to-date thermal information on the ocean. MI SSs were stratified into SS bins with midpoints from to 32.58C at 0.58C intervals, and each observation is assigned to the nearest SS midpoint. he empirical exponential MPI function described by Eq. (1) is derived with an SS cutoff of 28.58C since the flattening or decrease of maximum winds over the warmest waters occurs at this SS. he resulting coefficients are given by A kt (or m s 21,where1kt m s 21 ), B kt, C C 21,and C. Figure 1 shows this MPI function together with the data used for its development. he highest MPI is given as 140 kt: MPI 5 A 1 Be C(2 0 ). (1) he coefficients in the MPI SS relation are different from those of Knaff et al. (2005). Knaff et al. used 1-min sustained maximum wind speed in JWC best-track data whereas the RMSC okyo 10-min sustained maximum wind speed is used here. SLHF and IRR predictors are computed as the average within a box of and within a radius of 100 km centered at the C position, respectively. he averaging areas optimize their correlation coefficients with intensity and intensity changes, as has also been suggested in previous studies (Chang et al. 1997; Rodgers et al. 1994).

4 FEBRUARY 2012 G A O A N D C H I U 243 FIG. 1. he empirical relationship between MPI (kt) and SS (8C). he relationship is derived from 9 yr ( ) of data and the individual data points used for its development are also shown. For operational purposes, SLHF and IRR are static predictors, since no coupled atmosphere ocean models are used for prognostic purposes. A linear stepwise regression procedure is used to select parameters from the potential predictor pool. A 99% statistical significance level based on an F test (e.g., Wilks 2006) is the threshold for an individual predictor to be added initially into the model. Once selected, a predictor can only be removed if its significance level becomes less than 98% after the addition removal of another predictor. hree regression models are developed. he first serves as a control. he stepwise procedure is applied on the original 16 SIPS predictors to select significant predictors and create a base regression model (hereafter referred to as BASE). he second model is developed by conducting a stepwise procedure on the original predictor pool plus two more satellite-based parameters SLHF and IRR (hereafter referred to as SIPER). he stepwise procedure is also employed on the six predictors related to climatology and persistence to create another regression model called CLIPER, which is generally a baseline for evaluating the skill of the operational models. A model can be considered to produce a skillful intensity forecast if it has a smaller error than CLIPER. Assuming the independence of annual statistics, the samples in one year are used for verification and the samples in the other years are used for model development. As a result, for each CLIPER, BASE, and SIPER model, there are nine regression equations (one for each verification year), which may contain different sets of significant predictors due to different training samples. For consistency and a fair comparison, the final predictor set for each model (CLIPER, BASE, and SIPER) is the union of all the significant predictors identified in the nine regression equations. Each model is rerun using the same set of predictors for the nine verification years. he predictand as well as the predictors are normalized by subtracting their means and dividing by their standard deviations before regression; the resulting coefficients can be used to compare the relative contribution of each predictor directly. b. Model interpretation 1) 24-h MODELS able 3 lists the normalized coefficients associated with each predictor for each SIPER 24-h forecast equation. he numbers of samples used to develop the regression equations are shown in parentheses at the top of the table. he 99% statistically significant predictors for each verification years are indicated in boldfaced italics in able 3. here are around 1100 samples for training in each model. he 24-h SIPER models contain 10 statistically significant predictors: initial intensity (MWS0), previous 12-h intensity change (DMWS), ABLE 3. Normalized regression coefficients in the 24-h SIPER forecast model during different verification years. he predictors are listed in the first column and the verification years are listed in the top row. Here, N (shown in parentheses) is the number of samples used to develop the equation. Coefficients significant above the 99% statistical significance level from an F test are indicated in boldfaced italics. Year (N) Predictor 2000 (1126) 2001 (1102) 2002 (1065) 2003 (1097) 2004 (1039) 2005 (1108) 2006 (1111) 2007 (1144) 2008 (1176) 1) MWS ) DMWS ) JDAY ) LA ) PO ) RHLO ) RHHI ) SHR ) SLHF ) IRR

5 244 W E A H E R A N D F O R E C A S I N G VOLUME 27 Year (N) Predictor 2000 (759) 2001 (736) ABLE 4. As in able 3, but for 48-h forecasts (702) 2003 (726) 1) MWS ) DMWS ) JDAY ) LA ) PO ) RHLO ) RHHI ) REFC ) SHR ) USHR ) SLHF ) IRR (684) 2005 (738) 2006 (744) 2007 (770) 2008 (837) absolute value of Julian day minus 248 (JDAY), initial latitude of the storm (LA), potential (PO), hPa average relative humidity (RHLO), hPa average relative humidity (RHHI), hPa wind shear (SHR), SLHF, and IRR. DMWS, JDAY, PO, and SHR are significant for all the verification years. IRR is significant for the verification cases in the years 2001, 2002, 2003, 2005, 2007, and 2008 (six out of nine), and SHLF is significant for the 2001, 2004, 2006, 2007, and 2008 (five out of nine) verification cases. Correspondingly, MWS0, DMWS, JDAY, and LA are used to develop the CLIPER models and MWS0, DMWS, JDAY, LA, PO, RHLO, RHHI, and SHR are used in the BASE models. In all of these models, the four most important predictors are PO, DMWS, SHR, and MWS0. As expected, the persistence term, DMWS, is associated with a positive regression coefficient, since storms that have intensified in the previous 12 h tend to intensify in the next 24 h (Knaff et al. 2005). Intensity change is negatively correlated with MWS0 because weak storms are further from their MPIs and hence have more potential to intensify. Vertical wind shear has a negative impact on the intensification of Cs. One explanation for this is that the heat and moisture at upper levels are advected in a different direction relative to the low-level cyclonic circulation and therefore the ventilation of heat away from the circulation inhibits the development of the storm (Gray 1968). DeMaria (1996) proposed an alternate explanation: the tilt of the upper- and lower-level potential vorticities due to vertical wind shear produces a midlevel temperature increase near the vortex center. his midlevel warming is hypothesized to reduce the convective activity and thus inhibit storm development. he coefficient of JDAY is negative since this variable indicates the number of days from the peak of the typhoon season. LA is negatively correlated with intensity change since the SS generally decreases toward the north in the western North Pacific Basin. RHLO and RHHI can affect C intensification rates because high relative humidity in the middle atmosphere reduces the entrainment of dry air into the cumulus convection, which is a direct source of C energy. Nearly the same positive coefficients associated with the two new satellitebased predictors suggest that latent heat release in the atmosphere and latent heat transfer at the ocean-atmosphere interface have comparable effects on C development in 24 h. 2) 48-h MODELS he statistically significant predictors in SIPER selected from the potential predictor pool for 48-h intensity forecasts include MWS0, DMWS, JDAY, LA, PO, RHLO, RHHI, 200-hPa relative eddy flux convergence (REFC), SHR, the zonal component of the hPa wind shear (USHR), SLHF, and IRR (able 4). DMWS, PO, SHR, and IRR are significant for all of the verification years. SHLF is significant only for the years 2001 and here are around 700 samples for training of each 48-h model. he 48-h CLIPER models have the same predictors as in the 24-h version, and the 48-h BASE models contain 12 predictors. he additional predictors in the 24-h models are REFC and USHR. he predictors related to SS and vertical wind shear (PO and SHR) are most important. he contribution of the persistence term DMWS is weaker than in the 24-h SIPER model; persistence is more important for shorter-term forecasts. Relative to the 24-h SIPER model, the initial SLHF has a lower impact on the 48-h intensity change; however, IRR is significant in all of the verification years and has a larger normalized regression coefficient, suggesting the increasing role of initial innercore latent heating in longer-term forecasts.

6 FEBRUARY 2012 G A O A N D C H I U 245 Year (N) Predictor 2000 (522) 2001 (505) ABLE 5. As in able 3, but for 72-h forecasts (470) 2003 (489) 1) MWS ) DMWS ) JDAY ) LA ) LON ) PO ) RHHI ) REFC l ) SHR ) USHR ) SLHF ) IRR (464) 2005 (509) 2006 (509) 2007 (536) 2008 (572) he positive regression coefficients associated with USHR indicate that westerly shear (and westerly 200-hPa winds) is favorable for C intensification. his positive relationship is consistent with the finding in SIPS, but different from the negative relationships found in the Eastern Pacific basin and in the Atlantic basin (DeMaria and Kaplan 1994a; DeMaria and Kaplan 1999). here are two possible explanations for this. 1) he accompanying westerly winds from tropical upper-tropospheric troughs (Us; Sadler 1976, 1978) or midlatitude troughs on the north sides of the Cs weaken the easterly shear normally observed over the storm center and result in weak westerlies within the km annulus where the vertical wind shear is computed. his mechanism is likely more effective in the western North Pacific Basin, where most of the tropical cyclogenesis is associated with monsoon troughs (Zehr 1992; Briegel and Frank 1997; Ritchie and Holland 1999), which often have intense upper-level easterlies (Wang and Xu 1997) that inhibit storm intensification. 2) Some observations suggest that typhoon peak intensity often occurs at or near recurvature (Riehl 1972; Evans and McKinley 1998; Knaff 2009). he positive relationships between the momentum flux predictor REFC and the 48-h intensity change are expected, since REFC tends to be large when a C is moving toward an upper-level trough in the midlatitude westerlies or interacting with upper-level cold lows at low latitudes, and the large REFC that makes the upper-level circulation more cyclonic could result in an increase in the storm intensification rate (Holland and Merrill 1984; Molinari and Vollaro 1989; DeMaria et al. 1993). Note that REFC is not a significant predictor for the 24-h BASE models. 3) 72-h MODELS welve predictors [MWS0, DMWS, JDAY, LA, initial longitude of the storm (LON), PO, RHHI, REFC, SHR, and USHR] are selected as statistically significant predictors in SIPER from the potential predictor pool for 72-h intensity forecasts (able 5). DMWS, PO, SHR, and IRR are significant for all of the verification cases. SHLF is significant only for the 2001 and 2006 verification cases. Around 500 samples are available for each 72-h model development. he 72-h CLIPER models have one more predictor (LON) compared to the 24- and 48-h CLIPER models. he 72-h BASE and 72-h SIPER models have two additional predictors (REFC and USHR) compared to the 24-h version of SIPER, while RHLO, which is included in the 48-h model, is eliminated in the 72-h models. he predictors related to SS and vertical wind shear (PO and SHR), and the initial intensity (MWS0), are the most important, while the persistence term DMWS is less important in the 24- and 48-h versions of SIPER. As the 48-h intensity changes, the initial SLHF has a smaller impact on the 72-h intensity change than in the 24-h case. he other satellite-based predictor (IRR) has a comparable normalized regression coefficient as in the 48-h models, suggesting that the importance of the initial inner-core latent heating extends to 72-h forecasts. he positive regression coefficients associated with USHR and REFC indicate that westerly shear (and westerly 200-hPa winds) and large upper-level momentum flux are also favorable for C intensification in 72 h. he positive regression coefficient of LON indicates that the intensification rate would be higher if the storms are at higher longitudes. A possible explanation is that storms situated farther west are at their initial stages and, hence, have a greater potential to intensify. Furthermore, SSs are higher to the west over the western North Pacific. 4. Neural network models he nonlinear response of the C intensities suggests the use of nonlinear models. Models with nonlinear terms, such as those developed by Knaff et al. (2005),

7 246 W E A H E R A N D F O R E C A S I N G VOLUME 27 ABLE 6. MAEs (kt) with adjusted explained variances (R 2,%) in parentheses for CLIPER, BASE, and SIPER 24-h forecasts during different verification years with the number of samples (N). he best model for each verification year is indicated by boldfaced italics. FIG. 2. Architectural graph of an NN with one hidden layer. were developed and tested. Improvements of up to 4% are found for individual verification years. However, the NN method seems to work best. Back-propagation threelayer NN models for 24-, 48-, and 72-h intensity predictions (NN24, NN48, and NN72) are developed using the same predictors from the linear regression model SIPER. he NN architecture is indicated in Fig. 2. he input layers of the NN24, NN48, and NN72 models contain 10, 12, and 12 neurons, in that order, which correspond to the predictors used in the linear regression forecasting models for the same time period. All of the NN models have seven neurons in the hidden layers and one neuron in the output layers. he neuron in the output layers corresponds to the predictand (24-, 48-, or 72-h intensity change). he log-sigmoid transfer function is selected from the input layer to the hidden layer and the linear transfer function is used from the hidden layer to the output layer in all NN models. Similar to linear regression models, the NN models are developed using 8-yr data samples and the remaining 1-yr data samples are used for model verification. he samples in the development set of NN models are separated into three sets: a training set, a verification set, and a test set (Demuth et al. 2007). he numerical solution of the model is one that minimizes mean absolute error (MAE) while trying to avoid model overfitting. Overfitting can result in some unreasonable forecast models. 5. Model evaluation he linear regression models and NN models developed above are evaluated in this section. ables 6 8 show MAEs with adjusted explained variances (R 2 )for the CLIPER, BASE, and SIPER models at 24-, 48-, and 72-h forecasts for different verification years. he best model for each verification year is indicated in boldfaced Year N CLIPER BASE SIPER MAE R 2 MAE R 2 MAE R italics. he numbers of samples used for verification are also shown. able 9 indicates the verification statistics for the NN24, NN48, and NN72 models. For 24-h models, there are as few as 70 samples for verification year 2008 and up to 207 samples for verification year he BASE, SIPER, and NN24 models each can produce better forecasts than CLIPER. Among nine verification years, SIPER is the best model in six of the cases, and SIPER provides a 1% 3% MAE reduction relative to BASE for those six cases. he average MAEs for all verification years of CLIPER, BASE, and SIPER are 7.58, 7.18, and 7.16 kt, respectively, indicating that SIPER is only slightly better than BASE for the 24-h forecast. he MAEs of the NN24 models range from 5.97 kt in 2001 to 8.04 kt in 2007 (able 9). Figure 3 shows a direct comparison of three linear regression models (CLIPER, BASE, and SIPER) and the NN models at 24-h forecast interval for the nine verification years. Compared to 24-h SIPER, NN24 provides a 6% 10% MAE reduction for the nine verification years. he average MAE of NN24 for all nine verification years is 6.63 kt, or 7% improvement relative to that of 24-h SIPER (7.16 kt). Year ABLE 7. As in able 6, but for 48-h forecasts. N CLIPER BASE SIPER MAE R 2 MAE R 2 MAE R

8 FEBRUARY 2012 G A O A N D C H I U 247 Year ABLE 8. As in able 6, but for 72-h forecasts. N CLIPER BASE SIPER MAE R 2 MAE R 2 MAE R o examine the 24-h model s performance at specific stages of Cs, MAEs of the BASE, SIPER, and NN24 models are stratified by initial C intensity and 24-h intensity change for all nine verification years as indicated in Fig. 4. Figure 4 (top) shows MAEs as a function of initial intensity classified into 5-kt bins. he SIPER model outperforms the BASE model for typhoons with intensity between 80 and 105 kt except for the kt bins. he two models have very similar levels of performance for bins with lesser intensity. he NN24 model outperforms the SIPER and BASE models in all C stages except for the kt bin, where the BASE model produces the best forecasts. Figure 4 (bottom) indicates MAEs binned at every 5-kt change as a function of 24-h intensity change. Most of the samples are located within an intensity range of 620 kt. he SIPER model outperforms the BASE model for nearly all intensity change bins especially for the change between 35 and 40 kt. he NN24 model outperforms the SIPER and BASE models for all intensity change bins. he largest MAE reduction of the NN24 model occurs for the rapidly intensifying (24-h intensity change $ 35 kt) or rapidly decaying (24-h intensity change # 230 kt) storms. For 48-h models, there are 40 samples for verification year 2008 and up to 153 samples for verification year Each of the BASE, SIPER, and NN48 models can provide better forecasts than CLIPER except for SIPER in verification year Among the nine verification years of the three linear regression models, SIPER produces the best forecasts except for verification years 2001 and 2005, and MAE reductions of SIPER are 1% 3% compared to BASE for the seven years. he average MAE for all nine verification years of the 48-h CLIPER, BASE, and SIPER simulations are 10.88, 9.90, and 9.79 kt, respectively, indicating that SIPER is only slightly better than BASE. he MAEs of the NN48 models range from 7.81 kt in 2003 to kt in 2007 (able 9). Figure 5 shows the direct comparison of three linear regression models (CLIPER, BASE, and SIPER) and the NN models at ABLE 9. MAEs (kt) and adjusted explained variances (R 2,%) of neural network modes (NN24, NN48, and NN72) for 24-, 48-, and 72-h intensity forecasts during different verification years. NN24 NN48 NN72 Year MAE R 2 MAE R 2 MAE R h forecast intervals for the nine verification years. NN48 provides a 4% 17% MAE reduction relative to the 48-h SIPER simulation for the nine verification cases, with an average MAE of 8.68 kt, which shows an 11% improvement relative to that of 48-h SIPER (9.79 kt). For 72-h models, there are 26 samples for verification year 2008 and up to 108 samples for verification year Among the nine verification years, five SIPER models, three CLIPER models, and one BASE model provide the best forecasts, and MAEs of the SIPER models are 1% 6% smaller than those of the corresponding BASE models for the five cases. he 72-h models are less stable than shorter-term forecast models for 24- and 48-h intensity prediction. his is probably due to the decreasing sample size with 72-h models compared to 24- and 48-h models, as indicated in ables 3 5. he average MAE for all nine verification years of the 72-h CLIPER, BASE, and SIPER simulations are 11.63, 10.61, and kt, respectively, again, indicating that SIPER is only slightly better than BASE. FIG. 3. he MAEs (kt) of four regression models (CLIPER, BASE, SIPER, and NN) in different verification years at 24-h forecast intervals.

9 248 W E A H E R A N D F O R E C A S I N G VOLUME 27 FIG. 5. As in Fig. 3, but for the 48-h forecast interval. some predictors related to CLIPER, which have a greater impact on the intensity changes over shorter time periods. 6. Summary and discussion FIG. 4. (top) he 24-h BASE, SIPER, and NN24 model MAEs stratified by best-track initial intensity (MWS0) in 5-kt bins and (bottom) 24-h intensity change (DELV) in 5-kt bins for all nine verification years. Lower values of MAEs represent better forecasts. Dashed dotted lines represent the numbers of valid observations within a particular bin. he MAEs of the NN72 models range from 7.79 kt in 2003 to kt in 2004 (able 9). Figure 6 shows a direct comparison of the CLIPER, BASE, SIPER, and NN models at 72-h forecast intervals for the nine verification years. NN48 provides a 6% 28% improvement relative to 72-h SIPER for the nine verification years. he average MAE of NN72 for all nine verification years is 8.92 kt, which shows a 16% improvement relative to that (10.58 kt) of 72-h SIPER. Student s t tests show that, for all of the 9-yr forecasts, the differences between the NN models and each of the linear regression models for the three forecast intervals are statistically significant at a 95% confidence level, indicating that the NN models significantly outperform the linear regression models for the typhoon intensity forecast. he increase of the average MAEs with the longer forecast time intervals is primarily due to the inclusion of A new maximum potential intensity equation is derived using MI OI sea surface temperature data with high temporal and spatial resolutions (daily and 1 /48). his equation together with environmental information obtained from NCEP GFS FNL analysis, best tracks taken from RSMC okyo, and SLHF and IRR results derivedfromsatellitedataarethenutilizedtodevelop multiple linear regression models and NN models for western North Pacific C intensity forecasting at 24-, 48-, and 72-h intervals. Compared to the multiple linear regression models (BASE) with climatology, persistence, and environmental predictors only, the linear regression models that include additional satellite-based SLHF and IRR (SIPER) data provides a 1% 3% improvement in 24-h forecasting (for six out of nine cases) and 48-h forecasting (for seven FIG. 6. As in Fig. 3, but for the 72-h forecast interval.

10 FEBRUARY 2012 G A O A N D C H I U 249 out of nine cases) and a 1% 6% improvement in 72-h forecasting (for five out of nine cases). In terms of the average MAE for all nine of the verification years, SIPER is only slightly better, and the largest improvement of the satellite-enhanced SIPER model occurs for those rapidly intensifying storms with intensity change between 35 and 40 kt. he NN models developed using the same predictors as those used in the SIPER models outperform those of the SIPER models. he improvements are up to 10%, 17%, and 28% at 24-, 48-, and 72-h forecast times for individual verification cases, respectively. he average MAEs of all nine of the verification years show 7%, 11%, and 16% reductions relative to SIPER for the 24-, 48-, and 72-h NN models, respectively. he error reduction exists at almost all of the C stages and for all intensity change categories, and forecasts of the rapidly intensifying and rapidly decaying storms show the largest improvement. he NN model outperforms CLIPER at all forecast ranges (24, 48, and 72 h), indicating skill in all ranges. he other three models (BASE, SIPER, and NN) show improvements over CLIPER for the 24-h forecasts (see Fig. 3), indicating forecast skill at near range. However, for extended-range (72-h) forecasts, SIPER and BASE outperform CLIPER only six out of nine verification years (able 8), while CLIPER performs better than BASE and SIPER in three out of nine cases. It would be interesting to add IRR and SLHF to the static parameters in the CLIPER model and to compare their use in extended forecasts. he predictor IRR is significant for six out of the nine verification years at 24 h, and for all 48- and 72-h forecast intervals, indicating the importance of IRR. However, the impacts of the SLHF predictor are only significant for some of the verification years. he SLHF used in this study is based on a merged analysis and has relatively coarse temporal and spatial resolutions (daily and 18). Higher temporal and spatial resolution versions of SLHF may provide a more detail depiction of C ocean interaction. he progress on the third version of the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center Satellite-based Surface Fluxes (GSSF3; Shie et al. 2009) SLHF data with resolutions of 12 h and 1 /48, which is retrieved from a number of satellite datasets, may be potentially useful for further improvements to our ability to forecast C intensity. A study by Knaff and DeMaria (2005) found 5% 10% improvements by using a NN version of SHIPS for the Atlantic but they were not able to duplicate the NN improvement in operational forecasts. In fact, the NN results were inferior to the original multiple regression version of the SHIPS model when both versions were run. hey examined 2 yr s worth of independent cases using real-time GFS forecast fields and National Hurricane Center (NHC) forecast tracks. he use of GFS forecasts injected forecast errors into their C intensity model while perfect prognoses are used in our model development. hus, in order to ascertain the value of IRR and SLHF in C intensity model forecasts, our models need to be run in an operational setting to determine if the NN models with IRR and SLHF can actually improve C intensity forecasts operationally. Acknowledgments. his work constitutes part of a Ph.D.dissertationsubmittedbySGtotheInstituteof Space and Earth Information Science, Chinese University of Hong Kong, Hong Kong, China. Support from the CEMAPS project IS/058/09FP is acknowledged. he authors thank Dr. Roongroj Chokngamwong for helpful discussions and the two anonymous reviewers for their valuable comments. LSC thanks the NASA PMM/RMM program for its support. REFERENCES Baik, J.-J., and J.-S. Paek, 2000: A neural network model for predicting typhoon intensity. J. Meteor. Soc. Japan, 78, Briegel, L. M., and W. M. Frank, 1997: Large-scale influences on tropical cyclogenesis in the western North Pacific. Mon. Wea. Rev., 125, Camp, J. P., and M.. Montgomery, 2001: Hurricane maximum intensity: Past and present. Mon. Wea. Rev., 129, Chang, A.. C., L. S. Chiu, G. R. Liu, and K. H. Wang, 1997: Analysis of 1994 typhoons in the aiwan region using satellite data. Proceedings of COSPAR Colloquium on Space Remote Sensing of Subtropical Oceans (SRSSO), C.-. Liu, Ed., COSPAR Colloquia Series, Vol. 8, Elsevier Science, DeMaria, M., 1996: he effect of vertical shear on tropical cyclone intensity change. J. Atmos. Sci., 53, , and J. Kaplan, 1994a: A Statistical Hurricane Intensity Prediction Scheme (SHIPS) for the Atlantic basin. Wea. Forecasting, 9, , and, 1994b: Sea surface temperature and the maximum intensity of Atlantic tropical cyclones. J. Climate, 7, , and, 1999: An updated Statistical Hurricane Intensity Prediction Scheme (SHIPS) for the Atlantic and eastern North Pacific basins. Wea. Forecasting, 14, , J.-J. Baik, and J. Kaplan, 1993: Upper-level eddy angular momentum fluxes and tropical cyclone intensity change. J. Atmos. Sci., 50, , M. Mainelli, L. K. Shay, J. Knaff, and J. Kaplan, 2005: Further improvements to the updated Statistical Hurricane Intensity Prediction Scheme (SHIPS). Wea. Forecasting, 20, , J. A. Knaff, and C. Sampson, 2007: Evaluation of long-term trends in tropical cyclone intensity forecasts. Meteor. Atmos. Phys., 97, Demuth, H., M. Beale, and M. Hagan, 2007: Neural Network oolbox 5 users guide. MathWorks, Inc. [Available online at Emanuel, K. A., 1988: he maximum intensity of hurricanes. J. Atmos. Sci., 45,

11 250 W E A H E R A N D F O R E C A S I N G VOLUME 27 Evans, J. L., and K. McKinley, 1998: Relative timing of tropical storm lifetime maximum intensity and track recurvature. Meteor. Atmos. Phys., 65, Gao, S., and L. S. Chiu, 2010: Surface latent heat flux and rainfall associated with rapidly intensifying tropical cyclones over the western North Pacific. Int. J. Remote Sens., 31, Gray, M. W., 1968: Global view of the origin of tropical disturbances and storms. Mon. Wea. Rev., 96, Holland, G. J., 1997: he maximum potential intensity of tropical cyclones. J. Atmos. Sci., 54, , and R.. Merrill, 1984: On the dynamics of tropical cyclone structure changes. Quart. J. Roy. Meteor. Soc., 110, Huffman, G. J., and Coauthors, 2007: he RMM Multisatellite Precipitation Analysis (MPA): Quasi-global, multiyear, combined-sensor precipitation estimates at fine scales. J. Hydrometeor., 8, Jin, L., C. Yao, and X. Y. Huang, 2008: A nonlinear artificial intelligence ensemble prediction model for typhoon intensity. Mon. Wea. Rev., 136, Jones,. A., D. Cecil, and M. DeMaria, 2006: Passive-microwave enhanced Statistical Hurricane Intensity Prediction Scheme. Wea. Forecasting, 21, Kalnay, E., 2003: Atmospheric Modeling, Data Assimilation and Predictability. Cambridge University Press, 341 pp. Knaff, J. A., 2009: Revisiting the maximum intensity of recurving tropical cyclones. Int. J. Climatol., 29, , and M. DeMaria, 2005: Improvements in deterministic and probabilistic tropical cyclone surface wind predictions. U.S. Weather Research Program, Joint NOAA/Navy/NASA Hurricane estbed Final Rep., 7 pp., and R. M. Zehr, 2007: Reexamination of tropical cyclone pressure wind relationships. Wea. Forecasting, 22, , C. R. Sampson, and M. DeMaria, 2005: An operational statistical typhoon intensity prediction scheme for the western North Pacific. Wea. Forecasting, 20, Levitus, S., 1982: Climatological Atlas of the World Ocean. NOAA Prof. Paper 13, 173 pp. and 17 microfiche. Mainelli,M.,M.DeMaria,L.K.Shay,andG.Goni,2008:Application of oceanic heat content estimation to operational forecasting of recent Atlantic category 5 hurricanes. Wea. Forecasting, 23, Merrill, R.., 1987: An experiment in statistical prediction of tropical cyclone intensity change. NOAA ech. Memo. NWS NHC 34, 33 pp. Miller, B. I., 1958: On the maximum intensity of hurricanes. J. Meteor., 15, Molinari, J., and D. Vollaro, 1989: External influences on hurricane intensity. Part I: Outflow layer angular momentum fluxes. J. Atmos. Sci., 46, Reynolds, R. W., N. A. Rayner,. M. Smith, D. C. Stokes, and W. Wang, 2002: An improved in situ and satellite SS analysis for climate. J. Climate, 15, Riehl, H., 1972: Intensity of recurved typhoons. J. Appl. Meteor., 11, Ritchie, E. A., and G. J. Holland, 1999: Large-scale patterns associated with tropical cyclogenesis in the western Pacific. Mon. Wea. Rev., 127, Rodgers, E. B., S. W. Chang, and H. F. Pierce, 1994: A satellite observational and numerical study of the precipitation characteristics in western North Atlantic tropical cyclones. J. Appl. Meteor., 33, Sadler, J. C., 1976: A role of the tropical upper tropospheric trough in early season typhoon development. Mon. Wea. Rev., 104, , 1978: Mid-season typhoon development and intensity changes and the tropical upper tropospheric trough. Mon. Wea. Rev., 106, Shie, C.-L., and Coauthors, 2009: A note on reviving the Goddard Satellite-Based Surface urbulent Fluxes (GSSF) dataset. Adv. Atmos. Sci., 26, Wang, B., and X. Xu, 1997: Northern Hemispheric summer monsoon singularities and climatological intraseasonal oscillation. J. Climate, 10, Wentz, F. J., C. L. Gentemann, D. K. Smith, and D. B. Chelton, 2000: Satellite measurements of sea-surface temperature through cloud. Science, 288, Whitney, L. D., and J. S. Hobgood, 1997: he relationship between sea surface temperatures and maximum intensities of tropical cyclones in the eastern North Pacific Ocean. J. Climate, 10, Wilks, D. S., 2006: Statistical Methods in the Atmospheric Sciences. 2nd ed. International Geophysics Series, Vol. 59, Academic Press, 627 pp. Yang, F., H.-L. Pan, S. K. Krueger, S. Moorthi, and S. J. Lord, 2006: Evaluation of the NCEP Global Forecast System at the ARM SGP site. Mon. Wea. Rev., 134, Yu, L., X. Jin, and R. A. Weller, 2008: Multidecade Global Flux Datasets from the Objectively Analyzed Air Sea Fluxes (OAFlux) project: Latent and sensible heat fluxes, ocean evaporation, and related surface meteorological variables. Woods Hole Oceanographic Institution, OAFlux Project ech. Rep. OA , 64 pp. Zehr, R. M., 1992: ropical cyclogenesis in the western North Pacific. NOAA ech. Rep. NESDIS 61, 181 pp. Zeng, Z., Y. Wang, and C.-C. Wu, 2007: Environmental dynamical control of tropical cyclone intensity An observational study. Mon. Wea. Rev., 135,