Using Artificial Intelligence to Retrieve the Optimal Parameters and Structures of Adaptive Network-Based Fuzzy Inference System for Typhoon Precipitation Forecast Modeling

This study aims to construct a typhoon precipitation forecast model providing forecasts one to six hours in advance using optimal model parameters and structures retrieved from a combination of the adaptive network-based fuzzy inference system (ANFIS) and artificial intelligence. To enhance the accuracy of the precipitation forecast, two structures were then used to establish the precipitation forecast model for a specific lead-time: a single-model structure and a dual-model hybrid structure where the forecast models of higher and lower precipitation were integrated. In order to rapidly, automatically, and accurately retrieve the optimal parameters and structures of the ANFIS-based precipitation forecast model, a tabu search was applied to identify the adjacent radius in subtractive clustering when constructing the ANFIS structure. The coupled structure was also employed to establish a precipitation forecast model across short and long lead-times in order to improve the accuracy of long-term precipitation forecasts. The study area is the Shimen Reservoir, and the analyzed period is from 2001 to 2009. Results showed that the optimal initial ANFIS parameters selected by the tabu search, combined with the dual-model hybrid method and the coupled structure, provided the favors in computation efficiency and high-reliability predictions in typhoon precipitation forecasts regarding short to long lead-time forecasting horizons.


Introduction
Taiwan is located in the path of typhoons as they move in from the Western Pacific, and as a result, three to five typhoons hit Taiwan annually.Conventional methods used in disaster prevention mainly include engineered structures, but these methods have been proven to be mostly palliative, as typhoons are capable of damaging the hydraulic structures if their associated rainfall and flood discharge exceed the design limit.Therefore, effective methods of ameliorating typhoonrelated disasters need to include nonengineered disaster prevention programs, such as effective disaster forewarning and associated response mechanisms, which include the ability to identify the disaster before it occurs.One such method is the regional Quantitative Precipitation Forecast (QPF) information system.In the output of QPF, an accurate forecast of the duration of rainfall and its location is equally as important as the forecasted amount of precipitation.Information related to forecast precipitation can provide the public (and meteorological or water conservation departments) an estimate of expected precipitation in particular catchment areas, in addition to the ability to forewarn of potential debris flow hazards or of the need to employ flood control operations in reservoirs.Therefore, the QPF is crucial for use in disaster forewarning and associated responses.
Due to the powerful capability of artificial neural network (ANNs) to model nonlinear systems without the need to make any assumptions, in recent years ANNs have been used increasingly in applications for modeling hydrological processes, including precipitation forecasts (Govindaraju, [1]; 2 Advances in Meteorology Lin and Wu, [2]).The advantages of using ANNs include the ability to derive accurate short-term forecasts with low building costs, while the disadvantages are (1) unsatisfactory long-term precipitation forecasts with a significant time delay and (2) lower forecast accuracy in relation to higher amounts of precipitation.
Research related to the use of ANNs is summarized as follows.Hsu et al. [3] developed a modified counter propagation artificial neural network for transforming satellite infrared images to precipitation over a watershed, and Kuligowski and Barros [4] developed a feedforward neural network (FFNN) architecture with upper atmospheric wind direction and antecedent observed rainfall data to forecast short-term (0-6 h) precipitation for a specific location.Luk et al. [5] employed backpropagation neural networks (BPNN) in a short-term precipitation forecast in urban catchment areas, focusing on the effects of temporal and spatial information on the precipitation forecast.This study found that there was high forecast accuracy along with a reduced time delay when the input of the neural network included data from a certain number of adjacent rainfall stations.Grecu and Krajewski [6] applied BPNN and a statistical method to radar echo data in QPF and discussed the differences in results arising from the different spatial resolutions and echo sizes.Toth et al. [7] used the linear autoregressive model, BPNN, and nonparametric nearest-neighbors method for rainfall prediction by taking precipitation as the only input of the model; the results showed a significant improvement in long-term precipitation forecast using the ANN.Luk et al. [8] studied the performance of multilayer feedforward neural networks, partial recurrent neural networks, and time delay neural networks in precipitation forecasts within urban catchment areas.Results showed that the three models were all effective in precipitation forecasts and further indicated that precipitation time series were characterized with a shortterm memory and that the predicted precipitation presented time delays under large lead-times.Liu et al. [9] proposed a scheme of instantaneous rainfall estimation based on the radial basis function neural network trained with WSR-88D Doppler radar observations.Lin and Chen [10] applied the ANN to typhoon rainfall forecast, and their results showed that precipitation was significantly underestimated when the model input only included typhoon characteristic parameters but that the forecast accuracy was improved by adding appropriate data from adjacent rainfall stations.Valverde Ramírez [11] developed two precipitation forecast models in São Paulo State, Brazil.Their results showed that a better performance can be achieved with the ANN model compared to the multiple linear regressions (MLR) model, particularly when using nonlinear phenomenon such as precipitation forecasting.Hong [12] developed a hybrid model of RNNs (recurrent artificial neural networks) and SVMs (support vector machines), namely, RSVR, to forecast typhoon precipitation, where a chaotic particle swarm optimization algorithm (CPSO) was employed to choose the parameters of the SVR model.The results showed that the proposed model yielded a good forecasting performance.Lin and Wu [2] combined the self-organizing map (SOM) and the multilayer perceptron network (MLPN) to develop a hybrid typhoon precipitation forecast model.In the proposed model, a data analysis technique was developed based on the SOM, which is able to perform cluster analysis and discrimination analysis in one step, and the MLPN is used as a nonlinear regression technique to construct the relationship between the input and output data.The results showed the proposed model to deliver a more precise forecast than the model developed using a conventional neural network approach.
Based on this previous research, there is an evident lack of studies related to instantaneous typhoon precipitation forecasts using ANN techniques.However, the accurate forecast of rainfall during typhoons is extremely difficult because of the variable tracks of typhoons and the complicated interaction between typhoon circulations and diverse terrain (Huang et al., [13]).As a result, the accurate and effective application of ANN techniques (within the field of artificial intelligence) in precipitation forecasts for multiple lead-times during typhoons is a subject that urgently requires development and scientific breakthrough.Typhoon precipitation is characterized by a fuzzy, chaotic, and highly uncertain nonlinear system.The adaptive network-based fuzzy inference system (ANFIS) which is widely used as solution for various real world problems contains fuzzy logic reasoning and learning, and is therefore considered to be a promising candidate technique for use in typhoon precipitation forecasts.Current available literature related to ANNs shows that the selection of network architecture (both input selection parameters and network architecture) and its efficient training are very time consuming, and major obstacles are inherent in its day to day applications (Nasseri et al., [14]).The most commonly used method to retrieve the parameters and structures of ANNs is a random trial-anderror method, which is time consuming and proves difficult in retrieving the optimal parameters, and there is only limited previous research performed on the optimization of parameters and structures of ANNs.Of these researches, El-Fallahi et al. [15] employed the generalized reduced gradient (GRG) method to optimally select initial weights and bias values of BPNN.Results showed that the GRG method was able to effectively search the optimal weights and bias values within a reasonable computing time.In addition, Nasseri et al. [14] established the precipitation forecast model based on FFNN and used the back propagation learning algorithm and the genetic algorithm in network training and optimization, respectively.Results showed the superior performance of the genetic algorithm in precipitation forecasting by selecting optimal network inputs, compared with randomly selected inputs.Furthermore, Leahy et al. [16] applied a simulated annealing algorithm and tabu search to adjust the weights and the number of nodes of the network; results showed that the optimal ANN structures and parameters selected by the simulated annealing algorithm, along with a tabu search, significantly reduced the error and network complexity.In addition, this technique yielded promising forecast results.
However, these above-mentioned studies have not addressed the problems regarding optimization of the parameters and structures of ANFIS.The objective of this study, therefore, is to develop multiple mechanisms and methods for building a typhoon precipitation forecast model based on ANFIS and BPNN, and to compare their respective advantages and disadvantages.Through the study of various methods, which include the optimal selection of the ANFIS parameters and structures based on a merging of tabu search and subtractive clustering, a precipitation forecast for a specific lead-time using a single-model structure and dual-model hybrid structure, and a multiple leadtime precipitation forecast using coupled and noncoupled structures, we aim to improve the modeling efficiency for ANFIS, ameliorate the forecasted ability for high intensity precipitation, and enhance the forecasted accuracy for longer lead-time, respectively.

Development of Methodology
2.1.Procedures.The procedures used in this study are divided into four steps as shown in Figure 1.The detailed procedures are thoroughly described as follows.
Step 1.The optimal inputs for the typhoon precipitation forecast model were selected.A nonparametric statistical correlation analysis was applied to data related to meteorological and precipitation factors during previous typhoons.
After determining threshold values, the variables with higher correlation coefficients were selected as the model inputs.
Step 2. The optimal construction mechanism of ANFIS was determined (Section 2.3).Fuzzy subtractive clustering approach (Section 2.3.2) was used to determine the number of membership functions in the fuzzy rule base as well as the function parameters during the model construction.The initial adjacent radius in fuzzy subtractive clustering was selected using two methods: the conventional random trialand-error method and the tabu search algorithm, in order to rapidly, automatically, and accurately retrieve the optimal parameters and structures of ANFIS (Section 2.3.4).
Step 3. The optimal forecast model structure for a specific lead-time was retrieved.Precipitation forecast models based on ANFIS and BPNN were separately constructed using a dual-model hybrid structure for a future period from 1 to 6 hours.The dual-model hybrid structure was created by constructing and combining the higher and lower precipitation models.The single-model structure and dual-model hybrid structure were then compared to determine which model was superior.
Step 4. The optimal structure for the multiple lead-time precipitation forecast model was retrieved.Precipitation forecast models at different lead-times were coupled and compared with the noncoupled structure.

Selection of Model Inputs
2.2.1.Candidate Predictor.Three types of predictors were used in this study.
(2) Rainfall station factor: rainfall at the forecasting ground station and rainfall from adjacent stations.
(3) Factors calculated from the typhoon factors using fundamental meteorological formulae as follows.
(1) The distance between the typhoon and the forecasting station (()): this distance can be obtained using a conversion formula from longitude/latitude to distance: where lat  () and lat fos () are the latitudes of the typhoon center and the forecasting station at time  and lon  () and lon fos () are longitudes of the typhoon center and the forecasting station at time .
(2) The upper wind speed at the rainfall station (V()): the observed instantaneous grade 7 and grade 10 typhoon speed and radius were used in the present study to calculate the circulation precipitation in typhoon model.The relationship between wind speed and the distance from typhoon center is described by where  is an empirical value that varies with location;  is a constant number; and () is the distance between the typhoon and the forecasting station at time .The values of  and  can be solved from the known speed and radius values of grade 7 and grade 10 typhoons.The upper wind speed at the rainfall station V() is then obtained by substituting the distance between the typhoon and the station into (2).(3) Typhoon characteristic parameter (()): the typhoon characteristic parameter is calculated from the radius and radial speed of grade 7 and grade 10 typhoons using the following formula: where  7 () and  10 () are grade 7 and grade 10 typhoon radii, respectively, at time  and  7 () and  10 () are radial wind speeds of grade 7 and grade 10 typhoon, respectively, at time .

Selection of Model Inputs.
In this study, the favorable inputs used in the forecast model use correlation analysis for decision making, and Spearman's rank correlation coefficient (Spearman, [17]) is adopted as an analysis index.The analysis mechanism used for the correlation depends on the rank relationship of the time series of two variables, and hence this analysis can determine the correlation and suitability of input, regardless of the kind of relationship that exists between the candidate input and output; that is,  and  bench were used to verify the suitability of the solution selected by MAE: where Ŷ() is the forecasted value at time ; () is the measured value at time t;  bench () is the benchmark time series; and  is the number of data.Smaller values of MAE and RMSE imply a higher accuracy of the forecast model, and larger CC values indicate a closer coupling between the forecasted and measured series.The  bench value enables the effectiveness of the model to be more realistic and for the model to be formularized generally (which is particularly useful for hydrological forecast models).Precipitation is characterized by a continuous time system, in which precipitation at time  can be shifted along  time units to form a benchmark value at time  + , and in this way a benchmark series  bench can be constructed.Errors of the forecasted series and the benchmark series were both calculated from the measured values using (7).A negative value of  bench indicates a poorer performance of the forecast model compared with the benchmark series;  bench = 0 shows an equivalent performance of the forecast model and the benchmark series; and a positive  bench value means that the forecast model is superior to the benchmark series.

Establishment of the Optimal ANFIS Construction Mechanism.
The most commonly used method for retrieving ANN parameters and structure is merging the random trialand-error method with conventional parametric searching algorithm, which is used to identify the most appropriate initial structural parameters and train the best data-driven parameters under the random set up structure.However, the random trial-and-error method is time consuming, and the retrieval of the optimal parameters and structure becomes difficult if the initial solutions are set improperly and if the modeled system is enormous.Therefore, the optimization method was employed in the present study to retrieve the optimal parameters and structures of ANFIS.Specifically, the tabu search was first applied to optimize the adjacent radius in the subtractive clustering algorithm, and subtractive clustering was then used to determine the fuzzy membership function parameters and network structures of ANFIS.Finally,  the parameters of the consequent layer were estimated by the least squares estimator to complete the ANFIS model.The model identified was tested at a lead-time of +1 in the present study, as shown in Figure 2. ANFIS, subtractive clustering, and the tabu search used in the tests are described in detail below.

Introduction of ANFIS.
ANFIS was proposed by Jang (1993) [18] and is based on a fuzzy inference system constructed by combining the self-organization characteristics of a neural network.Hence, ANFIS integrates two algorithms to improve its accuracy and solves for the best parameters by employing capabilities of learning and selfadaption.ANFIS is composed of an input layer, a rule layer, a normalization layer, a consequent layer, and an output layer, as shown in Figure 3.The modeling tool can transform the fuzzy-complex process and phenomenon into artificial logic language that therefore is a potential approach for typhoon precipitation forecast.The computation and transmission of each layer is described as follows.
(1) Input Layer.This layer projects input to a group of fuzzy sets and estimates the values of a group of membership functions.We adopt a group of Gaussian functions as the membership functions, which can be expressed as follows: where   (  ) is the membership function,   and   are the antecedent parameters,  is the number of inputs, and   is the number of the fuzzy membership functions of input i.
(2) Rule Layer.This layer precedes the antecedent match of the fuzzy logic rule between variables and then applies a Tnorm product operation to obtain the weighted value of each rule; that is, where   is the weighted value and  is the number of rules.
(3) Normalization Layer.The node of this layer computes the output ratio between the node and all other nodes; that is, (4) Consequent Layer.The output of the consequent layer node is the product of the outputs of the normalization layer and the Sugeno fuzzy model (Takagi and Sugeno, [19]); that is, where   represents the consequent parameters and  0 is equal to 1.
(5) Output Layer.This layer sums the outputs of the previous layer to compute the model output; that is, ANFIS is a feedforward neural network and is constructed by supervised learning.The network parameters can be divided into antecedent parameters (nonlinear parameters:   ,   ) and consequent parameters (linear parameters:   ), and the model structure is determined by setting the number of membership functions in the input layer and the number of nodes of the rule layer.In tradition, the model structure and initial parameters are determined by random trial-and-error method, and the calibrated parameters after data-driven training are retrieved by least square estimator (LSE) and gradient steepest descent method.
To decrease the time for model construction in obtaining the best network structures and parameters, this study constructs ANFIS using hybrid algorithms including subtractive clustering (SC) and a least square estimator (LSE).The input and output vectors were first classified by subtractive clustering before training the model.The number of clusters obtained from the classification was set as the number of membership functions for node fuzzification at the various input layers and the number of nodes of the rule layers, and the th membership function node of each input layer node was connected to the corresponding th rule layer node.After determining the network structures, the center point and standard deviation of each cluster were taken as the initial parameters of the input layer membership functions (Gaussian function).The training data were then fed into the network with the consequent linear parameter set and the antecedent nonlinear parameter set solved by the least squares estimator and the gradient steepest descent method, respectively.The corresponding algorithm flowchart of the model construction is shown in Figure 4; this composite construction mechanism only requires the setting of adjacent radius in subtractive clustering and the subsequent training of network structures and parameters are automatically completed.The network structure significantly reduces the time required to retrieve the optimal number of fuzzy membership functions, number of rules, and network parameters, and the optimal network structure and parameters can be obtained after simply setting the adjacent radius between 0 and 1 (Jang, [18]; Jang et al., [20]).

Introduction of Subtractive
Clustering.Subtractive clustering was employed in the present study to construct fuzzy if-then rules in order to reduce the number of parameters of the fuzzy membership function in the ANFIS model.This was performed to establish a suitable rule base in the fuzzy inference system.Subtractive clustering was proposed by Chiu [21], in which every data point is treated as the candidate of the cluster center.Subtractive clustering is a fast and independent clustering method: the computational complexity is proportional to the number of data and is independent from the system dimension.For example,   ( = 1, 2, . . ., ) are  sets of data in an -dimensional space and the corresponding density measures  are defined as where the adjacent radius   is a positive number representing the distance near the center, and the data points outside the radius have minimum impact on the density measure.The density measure is calculated for each data point (  ), and the one with the highest density ( 1 ) is selected as the first cluster center ( 1 ).The definition of density measure is then modified to select the next cluster center.Assuming that   is the cluster center selected at the th round, and the corresponding density measure is   , the modified formula is as follows: where radius   has the same definition as   and is usually set as 1.5  so that the selected center will not be too close to that of the previous one.The above procedure of cluster center selection is repeated until a termination condition is reached or there are a sufficient number of cluster centers.

Establishment of the Optimal
A smaller value of the objective function ( value) indicates a higher accuracy of the forecast model.The value of the adjacent radius is confined to between 0 and 1, and the calculations across the layers of ANFIS that all belong to constraints are shown as follows: for  = 1, . . .,  ,  ∈ where  is the adjacent radius.In order to improve the forecast accuracy, single forecast models were constructed for each lead-time with no connections between them and were classified as noncoupled structures, wherein the singlemodel structure and the dual-model hybrid structure were first compared to determine the subsequent research topic regarding forecast model structure for a specific lead-time (the multiple lead-time precipitation forecast models based on coupled and noncoupled structures are described in Section 2.5).

Introduction of Tabu Search.
The tabu search was proposed by Glover [22] and Glover and Laguna [23] and guides the search direction and region using different types of memory.During the search, a search direction or region can be favored or prohibited according to the memory and rules.Additionally, the search can exit at a local optimum region and avoid repeated searches through the definition of a tabu list that includes the type and length of the search variables and the associated objective function value, where every iteration only searches to find the best candidate solution.
Hence, this search mechanism is not only able to significantly improve the search efficiency and accuracy but also able to obtain the best global solution.The tabu search was thus applied to optimize the adjacent radius in subtractive clustering, as required in the ANFIS model construction.A diagram of the method is illustrated in Figure 2.

Retrieval of the Optimal Forecast Model Structure for a
Specific Lead-Time.During typhoons, an accurate forecast of higher amounts of precipitation is considered to be much more important than amounts of lower precipitation.
According to previous research, precipitation has often been underestimated in higher rainfall, indicating different underlying mechanisms involved in higher and lower precipitation (Luk et al. [5]; Lin and Chen [10]; Lin and Wu [2]).To address this issue, the present study optimized the forecast of higher precipitation to reduce the forecast error.In addition, it was expected that the forecast error of higher precipitation would be lower than that obtained from the single-model structure.
According to the Central Weather Bureau in Taiwan, a heavy rainfall is defined as 24-hour accumulated rainfall of 50 mm or above, during which time there is at least one hour where precipitation amounts to 15 mm or above.In this study therefore, high precipitation was defined as an amount of 15 mm/hr or above.The dual-model hybrid method proposed in the present study constructs both lower and higher precipitation forecast models that can be subsequently selected within instantaneous forecasting.The lower precipitation model is usually initially applied, and if the forecasted precipitation is less than 15 mm/hr, the forecast value is directly accepted.However, if the forecasted precipitation is greater than 15 mm/hr, the higher precipitation model is used instead to obtain the corresponding forecast value.In this section, the single-model method and the dual-model hybrid method are first compared within noncoupled structures in order to determine the model construction for a specific lead-time in the combined multiple lead-time forecast model, based on noncoupled and coupled structures (see Section 2.5).

Combined Multiple Lead-Time Precipitation Forecast
Model Based on Coupled and Noncoupled Structures.Due to the complexity and uncertainty of the relationship between the atmosphere and physiography, the accuracy of any precipitation forecast decreases with any increase in the forecast lead-time.In this present study, two structures are proposed to improve the accuracy of long-term precipitation forecasts by constructing a multiple lead-time precipitation forecast model.The first structure is a noncoupled model, where the inputs of the forecast models for different leadtimes are all instantaneous observed data, and the outputs of forecast models at each lead-time are not coupled with the model inputs at the next lead-time.The second structure is a coupled model, which better characterizes the relationship between series and provides enhanced forecasting information.The coupled structure is therefore expected to provide a more accurate precipitation forecast for a longer lead-time.Observed information includes the typhoon factors that are input and the meteorological and physical factors, whereas the precipitation factor couples both observed and forecast information.Taking the present study as an example, forecast precipitation at  + 1 is required for the  + 2 coupled model, where the forecast precipitation at  + 1 replaces the input precipitation at  in the  + 2 model, the forecast precipitation at  replaces the input precipitation at  − 1 in the  + 2 model, and so on.The noncoupled and coupled structures are shown in Figure 5.After model training, it will then be intended to investigate the performance of the two structures in practical applications using model verification and assessment indices from the testing stage.

Introduction of BPNN.
The most commonly used BPNN was employed in the present study to construct models based on a noncoupled single-model and dual-model hybrid methods, as well as a coupled structure.The BPNN-based forecast performance was then compared with that of the newly developed ANFIS model.BPNN, developed by Rosenblatt [24], is a single-layer perceptron with no hidden layers and is therefore incapable of solving certain problems.The conception of a multilayer neural network was not realized until the well-known back propagation learning method was proposed by Rumelhart and McClelland [25], which developed an effective strategy for solving large amount of node weights.The BPNN structure is composed of multilayer perceptron (MLP) and generally employs the error back propagation (EBP, also abbreviated as back propagation (BP)) algorithm as the learning algorithm.The combination of MLP and EBP is known as BPNN.The BP algorithm is classified as a multilayer FFNN, and the nonlinear mapping between inputs and outputs is handled by supervised learning.The commonly adopted BPNN structure is shown in Figure 6 and includes an input layer, a hidden layer, and an output layer.The input and output layers represent the input and output values with the number of nodes determined by the type of problem to be solved.The number of hidden nodes is generally retrieved using the trial-and-error method.In addition, the number of hidden layers can be increased from one to multiple layers based on the complexity of the problem.Input layer Hidden layer Output layer The nodes of different layers are connected by related weights.The input values are directly transmitted to the hidden layer through the input layer and converted into the outputs by the activation functions after weighted summation.The output values are fed back into the input layer in the same manner.
In BPNN, the output of the jth node in the th layer is the conversion function of the node outputs in the  − 1th layer: where    represents the jth node output in the th layer, which is the input for the first layer. is the activation function and net   is the weighted summation of the outputs from the  − 1th layer: where    is the node weight connecting the jth node in the th layer and the th node in the  − 1th layer and    is the bias value of the jth node in the th layer.In the present study, the activation functions at the BPNN hidden layer and output layer were set as the hyperbolic function and the linear activation function, respectively.As a supervised learning network, the BPNN aims to reduce the difference between network outputs and objective outputs.The objective function (i.e., cost function) E is defined as follows: where   is the target output value of the kth node and   is the network output value of the kth node.Therefore, the network learning process is a process that minimizes the objective function .In the present study, the optimization of BPNN parameters was handled by the conjugate gradient method to search for the minimum value of the objective function E.  7).

Data Used in Model
Construction.This study used instantaneous observed information, such as rainfall and typhoon characteristics, to forecast precipitation for six leadtimes from  + 1 to  + 6.The output variables were taken as the precipitation forecast of the two rainfall stations of  + 1 to  + 6 in the Shimen Reservoir catchment area: Yu-Feng and Hsia-Yun stations.A total of 25 typhoons were selected in the present study based on the criterion that a typhoon land warning was issued and the corresponding cumulated precipitation was greater than 50 mm.The model construction included three stages, namely, a training, validation, and testing stage (Table 1).

Model Inputs Selection. Correlation analysis was
applied in the present study to assess the correlation coefficient between each input factor and the precipitation at each lead-time for each rainfall station.The selected model inputs and corresponding correlation coefficients for Yu-Feng station and Hsia-Yun station are shown in Tables 2 and 3, respectively.Research conducted by Lin and Chen [10] revealed that excessive model inputs could introduce additional noise into the model, therefore 10 input factors were selected based on the correlation coefficients and subsequent tests were used to avoid excessive inputs while retaining the forecasting information necessary for the model.In addition, we aimed to include certain typhoon factors as the model inputs in order to supplement information (Lin and Chen [10]; Lin and Wu [2]) that could not be represented by the rainfall stations, and therefore three typhoon factors and seven rainfall station factors were selected as the inputs.
The results showed that the forecast factor mostly used was precipitation information, and that typhoon factors and calculated physical factors were also indispensable.These parameters are considered to be essential for use in eliminating the time lag and for characterizing rainfall.

Determination of the Optimal ANFIS Construction
Mechanism.The present study combined a tabu search and subtractive clustering to determine the optimal structures and parameters of ANFIS.The adjacent radius in subtractive clustering was optimized using two methods: the traditional random trial-and-error method and the tabu search.The MAE values of the forecasted precipitation at  + 1 during validation were used to compare the two methods and to identify the superior model.The 10 inputs and one output (11 variables in total) were included in the proposed model, and, therefore, the number of decision variables for the adjacent radius was set as 11.The random trial-and-error method generated 11 random numbers between 0 and 1 each time for the value of adjacent radius in subtractive clustering, followed by model parameter training using LSE and the gradient steepest descent method.A total of 1000 models were generated and trained, and the one with the smallest MAE value was selected as the optimal model in the random trial-and-error method.The tabu search first generated 50 sets of random solutions for training and then selected the one with the smallest MAE value as the initial solution of the tabu search.The optimal solution was then obtained by conducting an additional 20 iterations of the tabu search, with a step shift of the decision variables set at 0.05.
The adjacent radius in subtractive clustering was found using the tabu search, and validation of the MAE values of the precipitation forecast model at  + 1 at Yu-Feng station and Hsia-Yun station is illustrated in Figure 8.These figures show that the MAE values of the constructed forecast models significantly decreased with an increase in the number of iterations.In addition, the optimal MAE values were markedly lower than those from the random trial-and-error method.Furthermore, the tabu search required only 490 iterations, which was less than half of the computation amount and time required by the random trial-and-error method.The testing results shown in Table 4 also indicate that the tabu search was able to retrieve the optimal ANFIS parameters and structures more efficiently than the random trial-anderror method.The lower testing values of MAE and RMSE and higher testing value of CC and  bench received from the tabu search approach represent that the constructed ANFIS-based precipitation forecast model using the tabu search approach is more accurate, stable, and effective than the random trial-and-error method.Therefore, the tabu search was selected for determining the optimal adjacent radius.

Retrieval of Optimal Forecast Model Structure for a Specific Lead-Time (1) ANFIS-Based Single and Hybrid Precipitation Forecast
Model for a Specific Lead-Time.In this section, the noncoupled single-model structure and the dual-model hybrid structure were compared for each lead-time.The construction method used for the optimal model was the same as that used for constructing the coupled model for each lead-time.
The forecast performance of the single and hybrid models at Yu-Feng and Hsia-Yun stations are compared in Figures 9 and 10.Based on the lower MAE value of assessment index of dual-model hybrid structure than single-model, the results show that the hybrid model effectively improves        the forecast accuracy especially for high intensity precipitation by between  + 1 and  + 5.The validated and testing results of Yu-Feng station and Hsia-Yun station all can be verified.However, lower forecast accuracy at Yu-Feng station at +6 was observed because the forecast at a greater lead-time was prone to a time lag, and when using a combination of the higher and lower precipitation models, the precipitation forecast by the lower precipitation model was actually forecast by the higher precipitation model, which increased the forecast error.Based on test results, the dual-model hybrid structure of the ANFIS forecast model outperformed that of the single-model structure and was therefore selected to construct the coupled forecast model for multiple lead-times.
(             the simulated mechanism of ANFIS owns the higher ability of tolerance and adaption for the transmitted error based on the fuzzy-adaptive routing structure, the coupled structure of ANFIS can perform better forecasted accuracy, stability, and effectiveness than BPNN.The number of nodes in hidden layer of the optimal BPNN-noncoupled precipitation forecast model for Yu-Feng station from lead-time 1 hour to 6 hours is 8, 8, 9, 9, 9, and 9, respectively, and, for Hsia-Yun station, is 8, 8, 9, 9, 10, and 10, respectively.In BPNN-noncoupled forecast model, the optimal activation function of hidden layer for shorter lead-time 1 hour to 3 hours is hyperbolic function, but, for longer lead-time 4 hours to 6 hours, the optimal activation function is linear function.It represents that the real-time observed meteorological rainfall information can retrieve the future precipitation mechanism, and linear activation function can better solve the problem of underestimating future precipitation in longer lead-time.Moreover, the number of rules derived from the optimal ANFIS-coupled precipitation forecast model for Yu-Feng station from lead-time 1 hour to 6 hours is 6, 7, 7, 8, 8, and 8, respectively, and for Hsia-Yun station is 7, 8, 8, 8, 9 and 9, respectively.The optimal model structure no matter ANFIS or BPNN all would be getting more complex along the increase of forecasted lead-time both for Yu-Feng station and for Hsia-Yun station.It represents that because the rainfall uncertainty of longer forecasted lead-time is too high, the forecast model including ANFIS and BPNN needs more rules and network connections, respectively, to retrieve the typhoon-rainfall mechanism.Furthermore, because the location of Hsia-Yun station is more near at marginal basin than Yu-Feng station, under the effect of terrain and lower representative information offered by adjacent observation station, the forecast model for Hsia-Yun station including ANFIS and BPNN needs more complex rules and network structures than Yu-Feng station, respectively, to retrieve the typhoon-rainfall mechanism.
From the perspective of model construction time, the ANFIS model was found to be significantly better than the BPNN model.Based on the computer computation in the present study, it took BPNN approximately 40000 seconds to construct a forecast model for a single lead-time, while the ANFIS only required approximately 1800 seconds.This means that without parallel computation using multiple computers, BPNN will require three days to complete the model construction from  + 1 to  + 6, and ANFIS will only require three hours.As a result, the ANFIS-based precipitation forecast model is considered to be superior based on its construction speed and convenience that can solve the time-consuming problem.Additionally, the overall performance of the models at each forecast lead-time presented little difference, and most of the validation and testing results of the ANFIS-coupled model were close to optimum.

Conclusion
Typhoon precipitation is characterized as a chaotic, fuzzy, highly uncertain, and nonlinear system.The adaptive network-based fuzzy inference system (ANFIS) has the ability to reason and learn using fuzzy logic.Therefore, the present study applied ANFIS, combined with multiple artificial intelligence techniques, to develop the optimal   construction mechanism using model parameters and structures for a typhoon precipitation forecast model of oneto six-hour lead-time.The proposed system was evaluated by a comparison with the most commonly used back propagation neural networks (BPNN).The present study first employed nonparametric correlation analysis to assess the correlation between multiple observed and forecasted variables (typhoon factors, calculated physical factors, and precipitation information from adjacent rainfall stations) and the precipitation of the forecasting station, so that the most appropriate input variables at each lead-time were identified.In order to effectively and accurately retrieve the optimal parameters and structure of the ANFIS-based precipitation forecast model, the tabu search was applied to identify the adjacent radius in subtractive clustering when constructing the ANFIS structure.The results of the tabu search were further compared with the random trial-anderror method to obtain a fast and automatic solution.To improve the forecast accuracy of higher precipitation, the single-model method and the dual-model hybrid method were used in the present study to forecast higher and lower precipitation for a specific lead-time.In addition, coupled structures were adopted for the construction of multiple leadtime precipitation forecast model to improve the accuracy of long-term precipitation forecast.Lastly, four forecast models for multiple lead-times were compared, including ANFIScoupled, ANFIS-noncoupled, BPNN-coupled, and BPNNnoncoupled.
The proposed method was applied to the Taiwan Shimen Reservoir catchment area with a study period from 2001 to 2009.The results showed that (1) subtractive clustering was able to effectively simplify the construction and algorithm parameters of ANFIS and the tabu search was able to effectively optimize the parameters in subtractive clustering and hence retrieve the optimal ANFIS parameters and structures.Compared with the random trial-and-error method, the tabu search significantly provided the favors in computation efficiency and high-reliability predictions in typhoon precipitation forecasts.(2) Due to different precipitation mechanisms occurring regarding short to long leadtime forecasting horizons, the present study combined the forecast results from higher and lower precipitation models and predicted typhoon precipitation for a specific lead-time based on the current precipitation level.The results showed that the dual-model hybrid method is superior to the singlemodel structure in overall MAE and forecasting accuracy during maximum precipitation.(3) The ANFIS model not only outperformed the BPNN model in model construction efficiency, but its forecast accuracy was also equal to that of the BPNN model.The ANFIS model showed a favorable performance in characterizing the highly complex nonlinear relationship between atmosphere and precipitation during typhoons.Therefore, the use of coupling between the inputs and outputs of the ANFIS forecast model at different leadtimes was able to effectively improve the accuracy of longterm precipitation forecasts.The forecast performance barely changed between the BPNN-coupled and BPNN-noncoupled structures for multiple lead-times forecasting, indicating that the BPNN model was relatively insensitive to updated precipitation values, and that the error tolerance of the BPNN algorithm mechanism was relatively low.

Figure 2 :
Figure 2: Flowchart of the developed methodology for determining the optimal parameter and structure of ANFIS.
Model.In this study, the tabu search was used in the present study to optimize the adjacent radius in subtractive clustering to obtain the optimal Advances in Meteorology ANFIS parameters and structures.The objective function minimizing the forecast error is as follows: Min  = ∑  =1      Ŷ () −  ()       .

Figure 4 :
Figure 4: Flowchart of training the optimal parameters and structure of ANFIS by combining subtractive clustering, least squares estimator, and gradient steepest descent method.

Figure 5 :
Figure 5: Noncoupled and coupled structures for precipitation forecast regarding short to long lead-time forecasting horizons.

Figure 7 :
Figure 7: Study area and the distribution of rainfall gauges.

Figure 8 :
Figure 8: Comparison of MAE of validation between tabu search approach and random trial-and-error method for construction of ANFISbased precipitation forecast model: (a) Yu-Feng station; (b) Hsia-Yun station.

Figure 9 :
Figure 9: Comparison of the forecast performance of the single and hybrid ANFIS-based models at Yu-Feng station: (a) validation; (b) testing.

Figure 10 :
Figure 10: Comparison of the forecast performance of the single and hybrid ANFIS-based models at Hsia-Yun station: (a) validation; (b) testing.

Figure 11 :
Figure 11: Comparison of the forecast performance of the single and hybrid BPNN-based models at Yu-Feng station: (a) validation; (b) testing.

Figure 12 :
Figure 12: Comparison of the forecast performance of the single and hybrid BPNN-based models at Hsia-Yun station: (a) validation; (b) testing.

Figure 13 :
Figure 13: The validation MAE values of the models at multiple forecast lead-times: (a)Yu-Feng station; (b) Hsia-Yun station.

Figure 15 :
Figure 15: The testing index values of the coupled and noncoupled models at Hsia-Yun station: (a) MAE and RMSE; (b) CC and  bench .

Figure 16 :Figure 17 :
Figure 16: The validation and testing precipitation forecast results at  + 1 of the ANFIS-based hybrid dual-model at Yu-Feng station: (a) validation; (b) testing.

Figure 18 :
Figure 18: The validation and testing precipitation forecast results at  + 2 of the ANFIS-based coupled hybrid dual-model at Yu-Feng station: (a) validation; (b) testing.

Figure 19 :
Figure 19: The validation and testing precipitation forecast results at +2 of the ANFIS-based coupled hybrid dual-model at Hsia-Yun station: (a) validation; (b) testing.
and Rank   and Rank   are the sort values of   and   in their individual time series of the variable, respectively.

Table 1 :
The adopted typhoon events among training, validation, and testing stage.The methodology proposed in the present study was applied to the Shimen Reservoir catchment area, which measures approximately 763.4 km 2 .The main stream within this area is the Dahan Creek, which is the upper stream of the Tamsui River.Rainfall stations in the Shimen Reservoir catchment area used in the present study were set up by the Northern Region Water Resources Office at the Water Resources Agency of Ministry of Economic Affairs and include locations in Shimen, Hsia-Yun, Gao-Yi, Ba-Ling, Galahe, Yu-Feng, Bai-Shi, Zhenxibao, Siciouish Mountain, and Chi-Duan (Figure 3.1.Study Area.

Table 2 :
The selected model inputs and corresponding correlation coefficients for Yu-Feng station.

Table 3 :
The selected model inputs and corresponding correlation coefficients for Hsia-Yun station.

Table 4 :
Comparison of the value of assessment indices by applying tabu search and trial-and-error method for retrieving the optimal parameter and structure of ANFIS on testing stage.