Taiwan, being located on a path in the west Pacific Ocean where typhoons often strike, is often affected by typhoons. The accompanying strong winds and torrential rains make typhoons particularly damaging in Taiwan. Therefore, we aimed to establish an accurate wind speed prediction model for future typhoons, allowing for better preparation to mitigate a typhoon’s toll on life and property. For more accurate wind speed predictions during a typhoon episode, we used cutting-edge machine learning techniques to construct a wind speed prediction model. To ensure model accuracy, we used, as variable input, simulated values from the Weather Research and Forecasting model of the numerical weather prediction system in addition to adopting deeper neural networks that can deepen neural network structures in the construction of estimation models. Our deeper neural networks comprise multilayer perceptron (MLP), deep recurrent neural networks (DRNNs), and stacked long short-term memory (LSTM). These three model-structure types differ by their memory capacity: MLPs are model networks with no memory capacity, whereas DRNNs and stacked LSTM are model networks with memory capacity. A model structure with memory capacity can analyze time-series data and continue memorizing and learning along the time axis. The study area is northeastern Taiwan. Results showed that MLP, DRNN, and stacked LSTM prediction error rates increased with prediction time (1–6 hours). Comparing the three models revealed that model networks with memory capacity (DRNN and stacked LSTM) were more accurate than those without memory capacity. A further comparison of model networks with memory capacity revealed that stacked LSTM yielded slightly more accurate results than did DRNN. Additionally, we determined that in the construction of the wind speed prediction model, the use of numerically simulated values reduced the error rate approximately by 30%. These results indicate that the inclusion of numerically simulated values in wind speed prediction models enhanced their prediction accuracy.
Ministry of Science and Technology, TaiwanMOST107-2111-M-019-0031. Introduction
A typhoon is a severe natural disaster that affects tropical and subtropical coastal countries, and it occurs most frequently in the northwestern Pacific Ocean. Taiwan is to the east of the Eurasian Continent and at the western side of the Pacific Ocean; its climate is intermediate between tropical and subtropical climate. Thus, typhoons frequently occur in Taiwan, generally in summer and fall. Typhoons affecting Taiwan typically develop at the sea surface southeast of Taiwan, and most typhoons are accompanied by torrential rains and strong winds [1]. Such rain and wind add to the damage from typhoons, posing a great threat to the transportation, economic, agricultural, and fishery activities in Taiwan even if the typhoon itself does not hit Taiwan. Therefore, typhoons constitute a serious natural disaster in Taiwan.
For example, the 2015 Typhoon Soudelor was the most destructive typhoon that occurred in Taiwan in recent history, with gust intensity exceeding 12 on the Beaufort wind force scale (32.7 m/s). Its strong gusts caused widespread damage to infrastructure, affecting gas supply, power and utilities, transportation and communication, and weather radar stations. Electricity was cut off in approximately 4.5 million households simultaneously during Typhoon Soudelor—the greatest recorded number in recent history. The economic loss from the typhoon was estimated to be as high as US$76 million [2].
Therefore, we aimed to establish an accurate wind speed prediction model for future typhoons, allowing for better preparation to mitigate a typhoon’s toll on life and property. In this study, cutting-edge machine learning (ML) techniques were used to improve predictive accuracy. In general, ML algorithms learn from a huge dataset, improving their ability to identify patterns in the data. Specifically, ML involves creating algorithms to make prediction from sets of unknown data. Given their ability to perform parameter adjustment and achieve optimization through self-learning, neural network algorithms in ML are particularly powerful [3]. Such algorithms have been extensively used in recent wind speed prediction models, and these models are increasingly data-driven due to developments in ML [4–11]. In the development of neural networks, multilayer perceptron (MLP) networks are a classic approach that are often used and compared with other neural network models. For example, Wei [12] compared the accuracy of MLP with that of adaptive network-based fuzzy-inference-system neural networks in the construction of typhoon wind speed prediction models.
Deep learning has become possible due to the exponential increase in computing power in recent years. This approach is the further derivation of multiple neural layers from the original neural layers of a model. Such derivation improves an algorithm’s ability to learn, better approximating the complex neural network structure of a human being. For example, Hu et al. [13] formulated multilayer deep neural networks that were trained using data from data-rich wind plants. These networks extracted wind speed patterns, and the mapping was finely tuned using data from newly constructed wind plants. Tiancheng et al. [14] proposed a sandstorm prediction method that considered both the effect of atmospheric movement and ground factors on sandstorm occurrence, called improved naive Bayesian convolutional neural network classification algorithm.
In the field of neural networks, recurrent neural networks (RNNs), which can analyze sequential (or time-series) data, have recently been developed [15–20]. RNNs are connectionist models with the ability to selectively pass information across sequence steps while processing sequential data one element at a time [21]. Therefore, RNNs are important, especially in the analysis of sequential data. A particular type of RNNs is long short-term memory (LSTM), a class of artificial neural networks where connections between units form a directed cycle [22] introduced the LSTM primarily to overcome the problem of diminishing gradients. LSTM creates an internal network state that allows it to exhibit dynamic temporal behavior. Unlike feedforward neural networks, RNNs can use their internal memory to process arbitrary sequences of input [23]. To model the time series of wind speed data, Byeon et al. [24] developed the LSTM for the prediction of typhoon wind speeds. Hence, feature enhancement from RNNs has been explored in wind prediction.
The widespread application of ensembles in numerical weather prediction (NWP) has helped researchers improve weather forecasts [25–30]. Numerical models can be used to calculate all climate parameters through atmospheric dynamics and numerical methods, allowing for simulated results to be explained in terms of physical relationships. NWP models simulate the atmosphere on user-defined grid scales as a moving fluid. Through several types of parameterizations, NWP models also account for the influence of subgrid physical processes on grid-resolved motions [31, 32]. NWP models, such as the Weather Research and Forecasting model (WRF), have been increasingly popular as a low-cost alternative source of data for such assessments of climate parameters [33]. WRF offers a wide variety of physical and dynamical elements to choose from; these elements must be put together to form model configurations, with which the model can be run [34]. However, because of imperfect models and uncertain initial boundary atmospheric conditions, errors exist in the NWP output [35, 36].
Recent studies have applied the NWP model to typhoons and tropical cyclones [37–39]. However, the prediction of severe meteorological phenomena (such as typhoons), which result from a multiplicity of mutually interacting multiscale processes, remains a major challenge for NWP systems [25, 40, 41]. These models are typically unable to predict wind intensity with satisfactory accuracy, even at short forecast times and a high horizontal resolution [42, 43]. Some studies have tried combining NWP with ML models to develop an integrated climate prediction model. For example, Zhao et al. [44] evaluated the performance and enhanced accuracy of a day-ahead wind-power forecasting system. The system comprised artificial neural networks and an NWP model.
Furthermore, to enhance the predictive accuracy of constructed ML models, other researchers have used numerically simulated results as input data for the construction of ML models. For example, Zhao et al. [45] developed the ARIMAX model, where wind speed results from the WRF simulation were chosen as an exogenous input variable. However, to the best of our knowledge, in the literature on short-term wind speed prediction for typhoons (or tropical cyclones), numerically simulated values have seldom been used as an input variable for ML prediction models. Therefore, in relation to the construction of a ML-based wind speed prediction model, we evaluated the improvements to predictive accuracy afforded by the use of numerically simulated values (by comparing between its use and nonuse).
Therefore, our study has two primary aims: (1) develop an ML- and neural network-based wind speed prediction model and compare the predictive accuracy of various neural network-based algorithms and (2) evaluate the improvements to predictive accuracy afforded by the use of numerical solutions obtained from NWP models (by comparing between its use and nonuse) in a typhoon-surface wind-speed prediction model.
2. Methodology and Algorithms
Figure 1 illustrates the flow of the construction, involving NWP numerical solutions, of our neural network-based typhoon wind-velocity prediction model. In the first stage, we collected data on the typhoon characteristics and ground wind speed of the research area’s historical typhoon events. In the second step, wind speed solutions were obtained from a wind field simulation of typhoons, involving an NWP numerical model. In our study, we also employed a WRF numerical model to simulate circulation distribution, thus obtaining the wind speed values in the research area. Two datasets could be built, one comprising typhoon data and the other comprising wind simulation results from the NWP model. The datasets were split into testing and training-validation subsets. The training-validation set was used for the learning of several ML-based wind velocity prediction models, and the testing set was used for the identification of the optimal prediction model among these models.
Flowchart of neural network-based typhoon wind velocity prediction model using NWP numerical solutions.
The best model among a set of ML neural network-based models, involving MLP, DRNN, and stacked LSTM, was determined. According to [46], real-time dynamics constitute the most challenging aspect of wind speed forecasting. We determined the stacked LSTM model to be the best for wind speed forecasting due to its appropriate handling of long- and short-term time dependency. In the final stage, the forecast accuracy of the stacked LSTM model was evaluated against other neural networks. We also evaluated whether forecast accuracy in ML-based models improved when NWP numerical solutions were used as input.
2.1. Frameworks Underlying the Proposed Neural Networks
In this section, we describe the neural network-based architectures, which used the MLP, DRNN, and stacked LSTM algorithms, that were adopted for model construction. As illustrated in Figure 2(a), the MLP is a typical type of feedforward backpropagation neural network that uses processing units placed in the input, hidden, and output layers [47–49]. Each unit (with an associated weight) in a layer is connected to the units in adjacent layers [50, 51]. In the study, to enhance learning efficacy—and by implication, approximation, and prediction accuracy—we added hidden layers to a simple MLP neural network; the MLP was trained through backpropagation. Generally, the weight updates between layers are calculated in terms of the stochastic gradient descent [52]. Specifically,(1)Δwijt+1=βΔwijt+η∂E∂wij,where wijt is the weight set connecting layers i and j at time t, Δw is the weight correction, η is the learning rate, β is a momentum coefficient, and E is a cost function that indicates the difference between the target and predicted values. In particular, η and β are hyperparameters for adjusting the spacing of weight correction.
Architecture of (a) MLP, (b) DRNN, and (c) stacked LSTM neural networks.
As illustrated in Figure 2(b), the multilayers of the RNN structure comprise an input layer, multiple recurrent layers, and an output layer. When the length of recurrent layer is 1, the framework is a simple RNN. Here, the recurrent network is based on the networks developed by [53]. In the RNN, the hidden units are connected to context units; in the successive time step, the units feed back into the hidden units. The hidden state at any time step can contain information from an (almost) arbitrarily long context window [21]. The DRNN model framework has multiple recurrent layers before the forwarding to a dropout layer and output layer at the final output. In this paper, the dropout layer excludes 10% of the neurons to avoid overfitting.
A stacked LSTM architecture is defined as an LSTM model that comprises multiple LSTM layers (Figure 2(c)). The stacked LSTM, also known as deep LSTM, was first formulated by [54] and was applied to speech recognition problems. Similar to the framework underpinning the DRNN model, the stacked LSTM model uses multiple LSTM layers that are stacked before the forwarding to a dropout layer and output layer at the final output. In a stacked LSTM, the first LSTM layer produces sequence vectors used as the input of the subsequent LSTM layer. Moreover, the LSTM layer receives feedback from its previous time step, thus allowing for the capturing of data patterns. The dropout layer also excludes 10% of the neurons to avoid overfitting.
The basic structure of LSTM, as illustrated in the LSTM layer in Figure 2(c), comprises an input gate it, output gate ot, forget gate ft, and memory cell ct. A single LSTM layer has a second-order RNN architecture that excels at storing sequential short-term memories and retrieving them at many time steps later [55]. An LSTM network is identical to a standard RNN, with the exception of the summation units in the hidden layer being replaced by memory blocks [56]. Equations (2)–(6) describe how output values are updated at each step [22, 46, 57]. Specifically,(2)ft=σWf⋅xt+Uf⋅ht−1+bf,(3)it=σWi⋅xt+Ui⋅ht−1+bi,(4)ot=σWo⋅xt+Uo⋅ht−1+bo,(5)ct=ft⋅ct−1+it⋅σWc⋅xt+Uc⋅ht−1+bc,(6)ht=ot⋅σct,where xt is the input vector; σ is the activation function; Wf, Wi, Wo, Wc, Uf, Ui, Uo, and Uc are the weight vector terms; bf, bi, bo, and bc are the corresponding bias terms; and ht and h−1 are the current and previous hidden vectors, respectively.
In this study, because the function converges quickly and has no problems with a vanishing gradient, we used the ReLU activation function in the middle layers (the hidden, recurrent, and LSTM layers) of the aforementioned three models. The ReLU is defined as the positive part of its argument; that is, f(x) = max (0, x), where x is the input to a neuron. The ReLU function has been recently used to replace the sigmoid function in neural networks, resulting in good performance and fast training times [58]. Moreover, in these proposed models, the input layers receive the observed and simulated meteorological values as their input.
3. Study Area and Data
Taiwan is located at the western Pacific Ocean, a region frequently in the path of typhoons. Typhoons in Taiwan are often accompanied by torrential rains, which severely endanger life and property. Thus, more accurate regional wind speed prediction is needed to improve typhoon protection. In this study, we chose northeastern Taiwan as the study area (Figure 3), wherein Taipei City and Yilan City are the two major cities.
Map of Taiwan’s northeastern region and the meteorological stations therein.
3.1. Data from Gauge and Typhoon Weather Stations
According to typhoon data from the Central Weather Bureau (CWB), between 2000 and 2018, 29 typhoons have directly struck Taiwan (Table 1; 2003, 2011, and 2018 were excluded from the table because typhoons did not directly affect Taiwan in those years). The tracks of these historical typhoons are illustrated in Figure 4. These typhoons can be classified according to the Saffir–Simpson wind scale depending on the intensity of their maximum sustained winds. A tropical storm has the wind speed range of 18–33 m/s, and category 1, 2, and 3 typhoons have the wind speed ranges of 33–43 m/s, 43–50 m/s, and 50–58 m/s, respectively. In these collected typhoons, the numbers for the category 1, 2, and 3 typhoons were 9, 5, and 8, respectively, and the number for the tropical storm was 7.
Typhoons affecting the study area from 2000 to 2017 (29 incidents).
Year
Typhoon
Period (UTC)
2000
Kai-Tak
6–10 Jul
2001
Toraji
28–31 Jul
2001
Nari
10–20 Sep
2002
Nakri
9–11 Jul
2004
Mindulle
29 Jun–3 Jul
2004
Aere
23–26 Aug
2004
Haima
12-13 Sep
2004
Nock-Ten
24–26 Oct
2005
Haitang
17–19 Jul
2005
Longwang
1-2 Oct
2006
Bilis
12–14 Jul
2007
Sepat
16–19 Aug
2007
Krosa
5–8 Oct
2008
Kalmaegi
16–19 Jul
2008
Fung-Wong
27–29 Jul
2008
Sinlaku
10–16 Sep
2008
Jangmi
27–30 Sep
2009
Morakot
6–9 Aug
2010
Namtheun
30-31 Aug
2012
Saola
30 Jul–3 Aug
2013
Soulik
10–14 Jul
2013
Trami
16–22 Aug
2013
Kong-Rey
27–30 Aug
2014
Matmo
21–23 Jul
2014
Fung-Wong
19–22 Sep
2015
Soudelor
6–9 Aug
2015
Dujuan
25–29 Sep
2016
Megi
25–28 Sep
2017
Nesat
26–30 Jul
Tracks of historical typhoon incidents that affected Taiwan.
The data included two types of datasets: one on typhoon dynamic characteristics and another (comprising 2621 hourly records) on surface wind speed observations. Surface wind speed data were released by the CWB and measured from eight gauge stations: Anbu, Banqiao, Keelung, Pengjiayu, Su-ao, Taipei, Tamsui, and Yilan (Figure 3). In addition, data on typhoon dynamics were released by the CWB. The data comprised six variables: pressure at the typhoon center, latitude of the typhoon center, longitude of the typhoon center, typhoon radius (i.e., the distance from typhoon center for points with winds greater than 15.5 m/s), moving speed of the typhoon, and maximum wind speed at the typhoon center.
3.2. Data from Numerical Solutions
WRF model simulations using typhoon data were used to generate wind speed values at each meteorological station. To set the initial field and boundary conditions, we used data from the Final Operational Global Analysis dataset. The dataset is a part of the Global Data Assimilation System of the US Government’s National Center for Environmental Protection. The grid was set to be a two-domain nested grid. The horizontal grid spacing for the coarser domain (118°E–126°E, 21°N–29°N) was 15 km, and that of the finer domain (121°E–123°E, 24°N–26°N) was 3 km. Both domains had 32 vertical levels. This study used the WRF physical parameters recommended by [59, 60]. They studied the tracks and rainfall of typhoons that have invaded Taiwan in addition to other physical parameters, and they conducted wind forecasts using the WRF model. They determined the following physical parameters to be suitable for wind speed forecasts: with regard to the planetary layer, those from the Yonsei University (YSU) scheme; with regard to microphysics, those from the WRF Single-Moment 5-Class (WSM5) scheme; for cumulus parameterization, those from the Kain–Fritsch scheme; and with regard to longwave radiation, those from the Rapid Radiative Transfer Model scheme.
After our simulation of a large number of typhoon events, the WRF model generated wind velocity simulations at eight ground gauge stations. In the verification of subsequent WRF wind outcomes, the output objective for the WRF model was the simulated wind value at an altitude of 10 m, because the ground gauge-station wind speed reached the observation value at an altitude of 10 m. Figure 5 displays the scatter plot of the WRF model simulation values and observed gauge-station values. According to Figure 5, several observation stations (Keelung, Anbu, Yilan, and Su-ao) have consistently higher wind speeds. These higher wind speeds are caused by their surrounding topographies—which cause them to face windwards—in conjunction with the anticlockwise circulatory flow characteristic of typhoons. Pengjiayu Station is similar to the four stations. However, being located on an island on the sea around northeastern Taiwan, the station is directly affected only by a typhoon’s circulatory flow and not the station’s surrounding topography; this station also records higher wind speed values. The Taipei, Banqiao, and Tamsui Stations are located within the Tamsui Basin and are thus fully enclosed by mountains. Because the circulatory flow of typhoons can be easily disrupted by the surrounding topography, the wind speed values observed at these stations are slightly lower.
Scatterplots depicting observations vs. WRF simulations at eight stations: (a) Pengjiayu, (b) Keelung, (c) Yilan, (d) Su-ao, (e) Banqiao, (f) Taipei, (g) Tamsui, and (h) Anbu.
For the four windward observation stations, wind speed simulation values were lower and approximately equal to actual values for high and low wind speeds, respectively. For the three observation stations in the Tamsui Basin, wind speed simulation values were approximately equal to the actual values, for all wind speeds.
To evaluate the simulation outcomes, this study used the mean absolute error (MAE) and root mean squared error (RMSE) that are defined as follows:(7)MAE=1n∑i=1nOi−Yi,RMSE=1n∑i=1nOi−Yi2,where Yi is the estimated value of record i, Oi is the observation of record i, and n is the total number of records.
The performance of the estimation indicators (MAE and RMSE) in terms of errors is presented in Figure 6. The estimation indicators for the Pengjiayu Station had the highest error values (MAE = 3.547 m/s and RMSE = 5.139 m/s), followed by those of the four observation stations (Keelung, Anbu, Yilan, and Su-ao) facing windward—their MAE and RMSE values had errors in the ranges 1.623–2.020 m/s and 2.405–2.728 m/s, respectively. The estimation indicators for the three observation stations (Taipei, Banqiao, and Tamsui) in the Tamsui Basin had the lowest error values; their MAE and RMSE values had errors in the ranges of 1.092–1.362 m/s and 1.422–1.654 m/s, respectively.
Simulation performance in terms of prediction errors.
4. Modeling
We used ML-based models to construct our wind speed prediction model for our study area. Observation stations at Taipei and Yilan were chosen as the test location. When constructing the model, the attribute data entered into the model included data on typhoon dynamics, data from ground meteorological observation stations, and data obtained from the aforementioned simulation. In this study, we performed data splitting for all typhoon episodes. Our training-validation set comprised data on 23 typhoon episodes between 2000 and 2013 (2093 records in total); our testing set comprised data on six typhoon episodes between 2014 and 2017 (528 records in total). Model training and validation were performed through 10-fold cross-validation, in which the training set was divided into 10 subsamples, one of which was retained for model verification, and the other nine were used for model training. In the verification process, each subsample must be verified.
As mentioned in Section 2.1, multilayer neural networks were used to construct our ML neural network-based models. We used the adaptive moment estimation optimization algorithm (also known as the Adam optimizer) to optimize the momentum and learning rate [61]. Generally, the Adam optimizer is more broadly applied in neural networks [62–64]. The Adam optimizer can be used instead of the classical stochastic gradient descent procedure to iteratively update network weights based on training data. We also calibrated the hyperparameters, specifically, the number of neurons in a middle layer and the length of the middle layers in the MLP, DRNN, and stacked LSTM models.
These hyperparameters were evaluated through trial and error. First, the number of neurons in a middle layer was adjusted from 10 to 100 until the curves of the RMSE of the RMSE errors were approximated. The middle layers of the MLP, DRNN, and stacked LSTM models are the hidden, recurrent, and LSTM layers, respectively. Subsequently, the prediction errors corresponding to the minimum RMSE were obtained in the form of the optimal solution for the number of neurons in the middle layer. When the lead time was 1 h, for the Taipei Station, the optimal numbers of neurons were 30, 50, and 50 for the MLP, DRNN, and LSTM models, respectively; for the Yilan Station, those numbers for the three models were 40, 30, and 30, respectively (Figures 7(a) and 8(a)). Subsequently, we calibrated the length of the middle layers in the networks, adjusting them to be between 1 and 10 layers. For the Taipei Station, the optimal lengths of the middle layers were 7, 5, and 4 for the MLP, DRNN, and LSTM models, respectively; for the Yilan Station, those for the three models were 7, 6, and 5, respectively (Figures 7(b) and 8(b)).
Sensitivity of model parameters on the Taipei gauge: (a) number of neurons in a middle layer and (b) length of middle layers.
Sensitivity of model parameters on the Yilan gauge: (a) number of neurons in a middle layer and (b) length of middle layers.
Using the aforementioned method, we conducted parameter testing for the prediction models. Tests were conducted using the training-validation set at six hourly lead times between 1 and 6 h. The RMSE performance values (for all lead times) for the three models are presented in Figure 9. The stacked LSTM model outperformed the MLP and DRNN models for all lead times and for both the Taipei and Yilan Stations. We further tested and fine-tuned these models using the testing set to confirm the accuracy and feasibility of each model.
Performance levels of neural network-based predictions for lead times between 1 and 6 h using a training-validation set for the (a) Taipei Station and (b) Yilan Station.
5. Evaluation5.1. Predictions of Forecasting Horizons
We tested and evaluated the MLP, DRNN, and stacked LSTM models using the testing set. The testing set comprised data (measured at the Yilin and Taipei Stations) on six typhoon episodes that occurred between 2014 and 2017. Figures 10 and 11 illustrate the time-series charts of the simulated and observed wind speed values for these typhoon episodes. The six typhoon episodes were Matmo in 2014, Fung-Wong in 2014, Soudelor in 2015, Dujuan in 2015, Megi in 2016, and Nesat in 2017. According to Figure 10(a), the maximum observed wind speeds for the Taipei Station were between 7.8 m/s (Fung-Wong) and 14.9 m/s (Soudelor). According to Figure 11(a), those for the Yilan Station were between 9.9 m/s (Fung-Wong) and 26.8 m/s (Soudelor). As mentioned, the windward Yilan Station is more susceptible to the circular flow of the typhoon, unlike the Taipei Station, which is shielded by the mountains of the Tamsui Basin. Therefore, wind speed was consistently higher at the Yilan station. Figures 10(a)–10(c) illustrate the prediction results for the Taipei Station at the lead times of 1, 3, and 6 h, respectively; for all prediction models, prediction accuracy was inversely related to the lead time. The data for the Yilan Station exhibited a similar trend (Figures 11(a)−11(c)). Such a situation is not atypical for prediction models. Intuitively, the further a model predicts into the future, the harder it is to obtain useful, real-time features for prediction. Therefore, to evaluate predictive accuracy, we used the error evaluation indicators to quantify each model’s rate of error in predicting the subsequent hour’s wind speed.
Predicted results of typhoons (2014–2017) in the testing set on the Taipei gauge for the lead times of (a) 1 h, (b) 3 h, and (c) 6 h.
Predicted results of testing typhoons (2014–2017) on the Yilan gauge for the lead times of (a) 1 h, (b) 3 h, and (c) 6 h.
To compute term-by-term comparisons of the relative error in the prediction with respect to the actual value of the variable, the mean absolute percentage error (MAPE) and root mean square percentage error (RMSPE) were calculated. MAPE and RMSPE can be calculated using the following formulae:(8)MAPE=1n∑i=1nOi−YiOi,RMSPE=1n∑i=1nOi−YiOi2.
Generally, MAPE is used to express the MAE as a percentage of the observations. MAPE is an unbiased statistic for measuring the predictive capability of a model [65]. RMSPE has the same properties as the RMSE but is expressed as a percentage [66].
Figures 12 and 13 illustrate the performance of each model in terms of the four error indicators (MAE, RMSE, MAPE, and RMSPE) for all lead times and for both Taipei and Yilan Stations. For the MLP model, the error rate was steeply and positively related to the lead time. However, for the DRNN and stacked LSTM models, this relation was positive but slighter. To explain the results using the RMSPE error as an example, the RMSPE error increased as the lead time increased. For the Taipei Station and for any given lead time, the MLP model exhibited the steepest increase in error, followed by the DRNN and RMSPE models. Specifically, when the lead time increased from 1 to 6 h, RMSPE increased from 1.494 to 2.616, from 1.386 to 2.181, and from 1.205 to 2.030 for the MLP, DRNN, and stacked LSTM models, respectively. Similar results were obtained for the Yilan Station. When the lead time increased from 1 to 6 h, RMSPE increased from 3.391 to 4.142, from 2.409 to 3.613, and from 2.127 to 3.457 for the MLP, DRNN, and stacked LSTM models, respectively.
Performance measures at lead times from 1 to 6 h on the Taipei gauge: (a) MAE, (b) RMSE, (c) MAPE, and (d) RMSPE.
Performance measures at lead times from 1 to 6 h on the Yilan gauge: (a) MAE, (b) RMSE, (c) MAPE, and (d) RMSPE.
Table 2 presents the average measurement values of the four error indicators for all lead times. In terms of the two unbiased statistical indicators, MAPE and RMSPE, the stacked LSTM and MLP models exhibited the most and least favorable performance, respectively. Stacked LSTM’s superiority is attributable to its model structure. Specifically, LSTM is a neural network that contains LSTM blocks. LSTM blocks can be described as a type of smart network unit that can memorize numerical values of different time lengths (to determine the quantity of useful information). Additionally, a gate in the LSTM blocks can help determine if the input is important enough to be remembered and if the input can be exported as outputs. Because the data involved in wind speed prediction are sequential, the decision on how far back the retention of memory data should go becomes consequential for predictive accuracy. By contrast, because the MLP model has no memory capacity, its predictions are based on limited and currently available information, thus decreasing its predictive accuracy. For the DRNN model, although it can receive memory information from long ago, the lack of a gate filter system can cause the receiving of too much (i.e., overly long) memory information, potentially undermining predictive accuracy.
Average performance measures of absolute error terms (MAE and RMSE) and relative error terms (MAPE and RMSPE) for 1–6 h predictions.
Station
Measure
MLP
DRNN
Stacked LSTM
Taipei
MAE (m/s)
1.241
1.070
0.928
RMSE (m/s)
1.666
1.439
1.228
MAPE
0.730
0.599
0.560
RMSPE
2.129
1.764
1.593
Yilan
MAE (m/s)
1.863
1.714
1.528
RMSE (m/s)
2.680
2.489
2.185
MAPE
1.151
0.974
0.927
RMSPE
3.643
3.030
2.847
5.2. Evaluation of Forecast Efficiency with and without Numerical Solutions
We evaluated whether the use of WRF simulation values as input affects the model’s accuracy in predicting wind speed. We focused on the stacked LSTM model because it was the best performing model. Figures 14 and 15 illustrate the performance of the stacked LSTM model in terms of the four error indicators (MAE, RMSE, MAPE, and RMSPE) for all lead times. We compared forecast efficiency with and without WRF simulated values. Data for the Taipei and Yilan stations are presented in Figures 14 and 15, respectively.
Performance measures of predictive accuracy for the use and nonuse of WRF simulation values on the Taipei gauge: (a) MAE, (b) RMSE, (c) MAPE, and (d) RMSPE.
Performance measures of predictive accuracy for the use and nonuse of WRF simulation values on the Yilan gauge: (a) MAE, (b) RMSE, (c) MAPE, and (d) RMSPE.
We determined that the use of WRF simulation values as input increased the model’s predictive accuracy. Therefore, the use of numerically simulated values as a part of the input data aids in the reduction of predictive error.
We define the improvement rates for the MAE and RMSE (denoted IRMAE and IRRMSE, respectively) as follows:(9)IRMAE=MAEwith−MAEwithoutMAEwithout,IRRMSE=RMSEwith−RMSEwithoutRMSEwithout,where MAEwith and MAEwithout are the MAE results on the use and nonuse of WRF simulation values, respectively, and RMSEwith and RMSEwithout are the RMSE results on the use and nonuse of the WRF simulation values, respectively.
Figure 16 details the improvement rates for the Taipei and Yilan Stations for all lead times. For the Taipei Station, IRMAE ranged between 25.5% and 29.7%, and IRRMSE ranged between 27.0% and 30.8%. For the Yilan Station, IRMAE ranged between 26.4% and 36.3%, and IRRMSE ranged between 28.4% and 35.7%. Generally, the average IRMAE (for all lead times) was 27.3% and 30.3% for the Taipei and Yilan Stations, respectively, and the average IRRMSE (for all lead times) was 28.7% and 31.1% for the Taipei and Yilan Stations, respectively. For both the Taipei and Yilan Stations, improvement was demonstrated in predictive accuracy, although the improvement rate was higher for the Yilan Station. In summary, although wind speed was difficult to predict accurately at the Yilan Station, the use of WRF simulation values increased the predictive accuracy of the model at the Yilan Station.
Improvement rate of predictive accuracy in terms of MAE and RMSE, comparing the use and nonuse of WRF simulation values for the Taipei and Yilan Stations.
6. Conclusions
To accurately predict the wind speed of future typhoons, we constructed a typhoon wind speed prediction model using cutting-edge ML techniques. RNNs have been recently developed as a type of neural network that can analyze sequential data. The structure of such networks facilitates the effective processing of wind speed-relevant climate data over an extended period. That is, the structure imbues RNN models with long-term memory capacity. Therefore, such networks are suitable for predicting typhoon wind speeds. LSTM is a type of RNN that allows the user to decide the memory time’s length. Additionally, LSTM gives users the option to filter output results, increasing LSTM’s predictive accuracy. According to current developments in deep learning, learning performance is enhanced when the layers of neural networks are deepened. Therefore, we used deep learning neural networks in this study. Additionally, we compared the performance of three types of RNNs—MLP, DRNN, and stacked LSTM—in predicting wind speed values. These three types of model structure differ by their memory capacity: MLPs are model networks with no memory capacity, whereas DRNNs and stacked LSTM are model networks with memory capacity.
We chose northeastern Taiwan as the study area, and the observation stations at Taipei and Yilan were selected as the study subjects. The results indicated that for both the Taipei and Yilan Stations and for the MLP, DRNN, and stacked LSTM models, prediction error is positively related to the prediction lead time (of which there were six, one for each hour between 1 and 6 h). In other words, the lead time is inversely related to predictive accuracy. Model networks with memory capacity (DRNN and stacked LSTM) are more accurate than those without memory capacity (MLP). Stacked LSTM was also discovered to be more accurate than DRNN. Stacked LSTM is better than the traditional RNN method. Because stacking the LSTM’s hidden layers deepens the model, model accuracy is increased; LSTM is closer to being a deep learning technique than its counterparts. Thus, stacked LSTMs is a stable technique for overcoming problems in wind velocity prediction.
Finally, we examined if model accuracy is increased by the use of WRF simulation values as an input variable in the wind speed prediction model. The stacked LSTM model was used in the evaluation, given its best performance in most evaluation tests of the current study. Evaluating predictive accuracy in terms of the four error indicators, namely, MAE, RMSE, MAPE, and RMSPE, verified that the inclusion of the numerically simulated input greatly increased model accuracy. Specifically, for all lead times, the use of numerically simulated values reduced the MAE and RMSE error rates of the Taipei Station by 27.3% and 28.7%, respectively. Similarly, the MAE and RMSE rates for the Yilan Station were reduced by 30.3% and 31.1%, respectively. Therefore, the inclusion of numerically simulated values in wind speed prediction models enhanced their prediction accuracy.
Data Availability
The Final Operational Global Analysis data were obtained from the Global Data Assimilation System of the US Government’s National Center for Environmental Protection, which is available at https://rda.ucar.edu/datasets/ds083.2/. The typhoon data were obtained from the Central Weather Bureau of Taiwan, which are available at https://rdc28.cwb.gov.tw/and https://e-service.cwb.gov.tw/wdps/.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This study was supported by the Ministry of Science and Technology, Taiwan, under Grant no. MOST107-2111-M-019-003. The author acknowledges the data provided by the Central Weather Bureau of Taiwan and the Research Data Archive (RDA) at NCAR.
WeiC.-C.Soft computing techniques in ensemble precipitation nowcast201313279380510.1016/j.asoc.2012.10.0062-s2.0-84881652537ChangC. H.WangU. T.FuH. C.LinU. C.2015Taipei, TaiwanNational Science and Technology Center for Disaster Reduction, Investigation Report148WeiC. C.PengP. C.TsaiC. H.HuangC. L.Regional forecasting of wind speeds during typhoon landfall in Taiwan: a case study of westward-moving typhoons2018910.3390/atmos90401412-s2.0-85045296850CadenasE.RiveraW.Campos-AmezcuaR.HeardC.Wind speed prediction using a univariate ARIMA model and a multivariate NARX model2016910.3390/en90201092-s2.0-84959329983KimM.ParkM.ImJ.ParkS.LeeM. I.Machine learning approaches for detecting tropical cyclone formation using satellite data20191110.3390/rs111011952-s2.0-85066743716KozarM. E.MisraV.PowellM. D.Hindcasts of integrated kinetic energy in Atlantic tropical cyclones: a neural network prediction scheme2016144124591460310.1175/mwr-d-16-0030.12-s2.0-84997018030MadhiarasanM.DeepaS. N.A novel criterion to select hidden neuron numbers in improved back propagation networks for wind speed forecasting201644487889310.1007/s10489-015-0737-z2-s2.0-84949495989NaikJ.SatapathyP.DashP. K.Short-term wind speed and wind power prediction using hybrid empirical mode decomposition and kernel ridge regression2018701167118810.1016/j.asoc.2017.12.0102-s2.0-85039452101PanapakidisI. P.MichailidesC.AngelidesD. C.A data-driven short-term forecasting model for offshore wind speed prediction based on computational intelligence2019810.3390/electronics80404202-s2.0-85066812474WeiC.-C.Conceptual weather environmental forecasting system for identifying potential failure of under-construction structures during typhoons2017168485910.1016/j.jweia.2017.05.0102-s2.0-85019935817PangG.HeJ.HuangY.ZhangL.A binary logistic regression model for severe convective weather with numerical model data2019201915612728110.1155/2019/6127281WeiC.-C.Surface wind nowcasting in the Penghu Islands based on classified typhoon tracks and the effects of the Central Mountain Range of Taiwan20142961425145010.1175/waf-d-14-00027.12-s2.0-84919682240HuQ.ZhangR.ZhouY.Transfer learning for short-term wind speed prediction with deep neural networks201685839510.1016/j.renene.2015.06.0342-s2.0-84934887601TianchengL.Qing-dao-er-jiR.YingQ.Application of improved naive Bayesian-CNN classification algorithm in sandstorm prediction in inner Mongolia2019201913517657610.1155/2019/5176576CortezB.CarreraB.KimY.-J.JungJ.-Y.An architecture for emergency event prediction using LSTM recurrent neural networks20189731532410.1016/j.eswa.2017.12.0372-s2.0-85039755181GuoL.LiN.JiaF.LeiY.LinJ.A recurrent neural network based health indicator for remaining useful life prediction of bearings20172409810910.1016/j.neucom.2017.02.0452-s2.0-85014137440LiuT.YuS.XuB.YinH.Recurrent networks with attention and convolutional networks for sentence representation and classification201848103797380610.1007/s10489-018-1176-42-s2.0-85046008176QureshiA. S.KhanA.ZameerA.UsmanA.Wind power prediction using deep neural network based meta regression and transfer learning20175874275510.1016/j.asoc.2017.05.0312-s2.0-85020035518YuanX.ChenC.JiangM.YuanY.Prediction interval of wind power using parameter optimized Beta distribution based LSTM model20198210.1016/j.asoc.2019.1055502-s2.0-85066951367ZhangY.WangX.TangH.An improved Elman neural network with piecewise weighted gradient for time series prediction201935919920810.1016/j.neucom.2019.06.0012-s2.0-85067024011LiptonZ. C.BerkowitzJ.ElkanC.A critical review of recurrent neural networks for sequence learning2015http://arxiv.org/abs/1506.00019HochreiterS.SchmidhuberJ.Long short-term memory1997981735178010.1162/neco.1997.9.8.17352-s2.0-0031573117ZhangQ.WangH.DongJ.ZhongG.SunX.Prediction of sea surface temperature using long short-term memory201714101745174910.1109/lgrs.2017.27335482-s2.0-85028961894ByeonW.LiwickiM.BreuelT. M.Scene analysis by mid-level attribute learning using 2D LSTM networks and an application to web-image tagging201563232910.1016/j.patrec.2015.06.0032-s2.0-84962765341FioriE.ParodiA.SiccardiF.Uncertainty in prediction of deep moist convective processes: turbulence parameterizations, microphysics and grid-scale effects2011100444745610.1016/j.atmosres.2010.10.0032-s2.0-79955664615MolteniF.BuizzaR.PalmerT. N.PetroliagisT.The ECMWF ensemble prediction system: methodology and validation19961225297311910.1002/qj.49712252905MurphyA. H.What is a good forecast? An essay on the nature of goodness in weather forecasting199382281293ShenX.JiangB.CaoQ.LinW.ZhangL.Sensitivity of precipitation and structure of typhoon Hato to bulk and explicit spectral bin microphysics schemes201920192536413StormB.BasuS.The WRF model forecast-derived low-level wind shear climatology over the United States Great Plains20103225827610.3390/en30202582-s2.0-77949410402TrachteK.Atmospheric moisture pathways to the highlands of the tropical Andes: analyzing the effects of spectral nudging on different driving fields for regional climate modeling201891145610.3390/atmos91104562-s2.0-85056766534CarvalhoD.RochaA.Gómez-GesteiraM.Silva SantosC.WRF wind simulation and wind energy production estimates forced by different reanalyses: comparison with observed data for Portugal201411711612610.1016/j.apenergy.2013.12.0012-s2.0-84890929317ChadeeX. T.SeegobinN. R.ClarkeR. M.Optimizing the weather research and forecasting (WRF) model for mapping the near-surface wind resources over the southernmost caribbean islands of Trinidad and Tobago20171010.3390/en100709312-s2.0-85025161443AwanN. K.TruhetzH.GobietA.Parameterization-induced error characteristics of MM5 and WRF operated in climate mode over the Alpine region: an ensemble-based analysis201124123107312310.1175/2011jcli3674.12-s2.0-79959327126DzebreD. E. K.AcheampongA. A.AmpofoJ.AdaramolaM. S.A sensitivity study of surface wind simulations over coastal Ghana to selected time control and nudging options in the weather research and forecasting model2019510.1016/j.heliyon.2019.e013852-s2.0-85063411797KalneyE.2003Cambridge, UKCambridge University PressChuY.LiC.WangY.LiJ.LiJ.A long-term wind speed ensemble forecasting system with weather adapted correction201691189410.3390/en91108942-s2.0-85020514684DiZ.GongW.GanY.ShenC.DuanQ.Combinatorial optimization for WRF physical parameterization schemes: a case study of three-day typhoon simulations over the Northwest Pacific Ocean201910523310.3390/atmos100502332-s2.0-85066300995WeiC.-C.Wavelet support vector machines for forecasting precipitation in tropical cyclones: comparisons with GSVM, regression, and MM5201227243845010.1175/waf-d-11-00004.12-s2.0-84863003352WuZ.JiangC.DengB.ChenJ.LiuX.Sensitivity of WRF simulated typhoon track and intensity over the South China Sea to horizontal and vertical resolutions2019387748310.1007/s13131-019-1459-z2-s2.0-85067523731CassolaF.FerrariF.MazzinoA.Numerical simulations of Mediterranean heavy precipitation events with the WRF model: a verification exercise using different approaches2015164-16521022510.1016/j.atmosres.2015.05.0102-s2.0-84936124965DavolioS.MastrangeloD.MigliettaM. M.DrofaO.BuzziA.MalguzziP.High resolution simulations of a flash flood near Venice2009951671167810.5194/nhess-9-1671-20092-s2.0-77953739238WeismanM. L.DavisC.WangW.ManningK. W.KlempJ. B.Experiences with 0-36-h explicit convective forecasts with the WRF-ARW model200823340743710.1175/2007waf2007005.12-s2.0-48649108162MigliettaM. M.RotunnoR.Application of theory to simulations of observed cases of orographically forced convective rainfall201214093039305310.1175/mwr-d-11-00253.12-s2.0-84867929965ZhaoP.WangJ.XiaJ.DaiY.ShengY.YueJ.Performance evaluation and accuracy enhancement of a day-ahead wind power forecasting system in China20124323424110.1016/j.renene.2011.11.0512-s2.0-84856532432ZhaoE.ZhaoJ.LiuL.SuZ.AnN.Hybrid wind speed prediction based on a self-adaptive ARIMAX model with an exogenous WRF simulation20169LiangS.NguyenL.JinF.A multi-variable stacked long-short term memory network for wind speed forecastingProceedings of the 2018 IEEE International Conference on Big DataDecember 2018Seattle, WA, USA4561456410.1109/bigdata.2018.86223322-s2.0-85062600212RostamiA.BaghbanA.MohammadiA. H.Hemmati-SarapardehA.HabibzadehS.Rigorous prognostication of permeability of heterogeneous carbonate oil reservoirs: smart modeling and correlation development201923611012310.1016/j.fuel.2018.08.1362-s2.0-85052912167SaghafiH. R.RostamiA.ArablooM.Evolving new strategies to estimate reservoir oil formation volume factor: smart modeling and correlation development201918110618010.1016/j.petrol.2019.06.0442-s2.0-85068263322RostamiA.ArablooM.KamariA.MohammadiA. H.Modeling of CO2 solubility in crude oil during carbon dioxide enhanced oil recovery using gene expression programming201721076878210.1016/j.fuel.2017.08.1102-s2.0-85029420723BehzadiS. S.PrakasvudhisarnC.KlockerJ.WolschannP.ViernsteinH.Comparison between two types of artificial neural networks used for validation of pharmaceutical processes2009195215015710.1016/j.powtec.2009.05.0252-s2.0-67650558726WeiC.-C.Comparing lazy and eager learning models for water level forecasting in river-reservoir basins of inundation regions20156313715510.1016/j.envsoft.2014.09.0262-s2.0-84908689236LoD.-C.WeiC.-C.TsaiE.-P.Parameter automatic calibration approach for neural-network-based cyclonic precipitation forecast models20157123963397710.3390/w70739632-s2.0-84940421508ElmanJ. L.Finding structure in time199014217921110.1207/s15516709cog1402_1GravesA.Generating sequences with recurrent neural networks2013http://arxiv.org/abs/1308.0850WeningerF.GeigerJ.WöllmerM.SchullerB.RigollG.Feature enhancement by deep LSTM networks for ASR in reverberant multisource environments201428488890210.1016/j.csl.2014.01.0012-s2.0-84900534601GravesA.2012Berlin, GermanySpringerGersF. A.SchmidhuberJ.CumminsF.Learning to forget: continual prediction with LSTM200012102451247110.1162/0899766003000150152-s2.0-0034293152NairV.HintonG.Rectified linear units improve restricted Boltzmann machinesProceedings of Machine Learning ResearchMay 2010Sardinia, Italy807814YehT. C.ChenS. H.HongJ. S.2008Taipei, TaiwanFunded by Ministry of Science and Technology of TaiwanNo. NSC96-2625-Z-052-003 in ChineseYehT. C.TerngC. T.LeeC. S.YangM. J.The forecast technique development studies on the typhoon track, rainfall and winds forecast over taiwan area: a study on the implementation of WRF typhoon forecasting component in the operational environment of CWB (II)2009Funded by Ministry of Science and Technology of TaiwanNo. NSC97-2625-M-052-002 (in Chinese)KingmaD. P.BaJ. L.Adam: a method for stochastic optimizationProceedings of the International Conference on Learning RepresentationsMay 2015San Diego, CA, USA113HanF.JiangJ.LingQ.-H.SuB.-Y.A survey on metaheuristic optimization for random single-hidden layer feedforward neural network201933526127310.1016/j.neucom.2018.07.0802-s2.0-85056303500SharmaA.Guided stochastic gradient descent algorithm for inconsistent datasets2018731068108010.1016/j.asoc.2018.09.0382-s2.0-85055211186YaoC.CaiD.BuJ.ChenG.Pre-training the deep generative models with adaptive hyperparameter optimization201724714415510.1016/j.neucom.2017.03.0582-s2.0-85017377432WangW.-C.ChauK.-W.ChengC.-T.QiuL.A comparison of performance of several artificial intelligence methods for forecasting monthly discharge time series20093743-429430610.1016/j.jhydrol.2009.06.0192-s2.0-68349105875SwansonD. A.TaymanJ.BryanT. M.MAPE-R: a rescaled measure of accuracy for cross-sectional subnational population forecasts2011282-322524310.1007/s12546-011-9054-52-s2.0-79957952788