The Economic Loss Prediction of Flooding Based on Machine Learning and the Input ‐ Output Model

: As climate change becomes increasingly widespread, rapid, and intense, the frequency of heavy rainfall and floods continues to increase. This article establishes a prediction system using feature sets with multiple data dimensions, including meteorological data and socio ‐ economic data. Based on data of historical floods in 31 provinces and municipalities in China from 2006 to 2018, five machine learning methods are compared to predict the direct economic losses. Among them, GBR performs the best with a goodness ‐ of ‐ fit of 90%. Combined with the input ‐ output (IO) model, the indirect economic losses of agriculture to other sectors are calculated, and the total economic losses caused by floods can be predicted effectively by using the GBR ‐ IO model. The model has a strong generalization ability with a minimum requirement of 80 pieces of data. The results of the data show that in China, provinces heavily reliant on agriculture suffered the most with the pro ‐ portion of direct economic losses to provincial GDP exceeding 1‰. Therefore, some policy implica ‐ tions are provided to assist the government to take timely pre ‐ disaster preventive measures and conduct post ‐ disaster risk management, thereby reducing the economic losses caused by floods.


Introduction
The acceleration of climate change has intensified the water cycle, affected rainfall patterns, and caused rising sea levels, triggering more frequent heavy rainfall and worsening flood situations [1].Due to global warming, research shows that the death toll and economic damage of floods have risen by around 75% and 200%, respectively [2].According to a comprehensive analysis published in Geneva, 23 July 2021 by the World Meteorological Organization (WMO), of the top 10 meteorological disasters, floods have inflicted the top three largest economic tolls and human losses around the world over the past halfcentury.Fatal floods, incurred by heavy rains, occur with a low frequency, a wide range of influence, and can cause serious economic losses [3].The economic loss from such natural disasters should be evaluated as soon as possible and effective countermeasures should be taken to manage natural hazards [4].Hence, if potential economic losses of floods can be accurately predicted prior to or just after the beginning of heavy rainfall, more time can be provided for people to take precautions and put in place measures to avoid and alleviate unnecessary losses and secondary hazards.
In response, relevant studies of economic loss assessment have been conducted.Generally, economic losses, which can be measured with monetary value, are classified into two sorts: direct economic losses and indirect economic losses [2].Direct economic losses refer to property losses of residences, asset losses of enterprises, loss of infrastructure, loss of natural resources, etc., while indirect economic losses arise from the suspension and reduction of production, investment premiums, material shortages, overstock, and losses from related industries [5].
For direct economic losses, assessments are normally conducted by using a flood stage-damage calculation [6], loss rate approach, comprehensive loss value method, Category-Unit Loss Functions [7], approach based on RS and GIS [8], and Regression Analysis Rapid Evaluation Model [9], etc.By far, most of the direct economic assessments are conducted post-disaster, while pre-disaster assessment is currently at an initial stage.Yet with the boom in artificial intelligence, advanced approaches, such as backpropagation neural networks (BPNN) [10], random forest, support vector regression (SVR) [11], and gradient boosted regression tree (GBR) [12], etc., are implemented to predict the direct economic losses of natural disasters.The most popular approach among these is BPNN.Although BPNN better fits the complex nonlinear function relationship and has higher assessment accuracy, it requires a huge amount of data, which is lacking in extreme weather data as the number of flood occurrences across history are limited [10].Another commonly used method, SVR-the extension of support vector machine (SVM)-is suitable for solving high-dimensional regression problems using a small number of samples, which is in accordance with the features of natural disasters [13].Additionally, ensemble learning-a commonly used algorithm-integrates several machine learning methods to reduce deviation and improve prediction accuracy [14].Extreme Gradient Boosting Regressor (XGB) has been demonstrated to be a good fit for evaluating economic losses of natural disasters, but the generalization ability is unsatisfactory [12].Sun et al. used GBR in the field of direct economic loss evaluation caused by storm surge disasters, improving the prediction accuracy and reducing model overfitting caused by small datasets [11].However, there has not yet been any report of research applying GBR to the disaster assessment of floods.
For indirect economic loss assessment, empirical analysis methods, Computable General Equilibrium (CGE) models, and Input-Output (IO) models are commonly used.The empirical analysis method assumes a certain proportional relationship between the indirect economic losses caused by a flood to different sectors and industries and its direct economic losses in the inundated area, but the result is relatively ineffective [15].The dynamic CGE model can realize a comprehensive assessment of the indirect economic impact of flood disasters through parameter simulation, flow restriction simulation, and variable shock simulation [16].However, CGE models are based on more restrictive assumptions, which typically assume optimizing behavior and equilibrium economy.Such assumptions are easily violated under real-world economic conditions [17].In contrast, the input-output approach, an improvement in methodologies, is more suitable for assessing economic losses caused by exogenous shocks [18].Zhang et al. implemented the IO model to analyze the indirect economic losses of floods in Hunan Province, China, in 1998 [2].It is proven that the model is an effective and flexible approach to assess indirect economic losses as it can select the number of sectors according to the retrieved data.
Therefore, in this paper, we establish a comprehensive prediction system to assess the direct and indirect economic losses caused by floods in 31 regions of China.There are two main contributions in this paper.Firstly, we designed a timely and effective pre-disaster prediction system with the GBR-IO model by using indicators which can be collected pre-disaster or at the beginning of a flooding episode and found the minimum amount of data required for the prediction system to be able to provide results with high effectiveness.Secondly, based on the meteorological features of natural disasters, we combined machine learning methods, using advanced regression prediction approaches, with the input-output model, a traditional economic method, resulting in what can have a profound impact on inter-disciplinary research.Thirdly, the method we used has a high generalization ability, meaning that it can be applied to other countries and regions which experience flooding, especially those with small datasets.

Study Area
Statistics show that Asia is the continent where floods are, by far, the most frequent and devastating natural disasters around the world [19].Among Asian countries, Chinawith a monsoon climate and major rivers of the world-is the most frequently affected country.Around the world, the frequency of floods and economic losses caused by floods in China ranked first, and the casualties and death toll ranked second during the last decade [20].Recently, in mid to late July 2021, rainstorms and devastating floods battered large portions of North China and the Huanghuai Region, especially Zhengzhou, the state capital of Henan province, causing the deaths of 302 people, 50 people to be missing, and 114.3 billion RMB yuan (US $17.7 billion) in direct economic losses [21].For studying floods or the economic losses caused by floods, China is ideal in terms of sample size and practical significance.Thus, in this research we focused on 31 provinces in China (excluding the Hong Kong, Macao, and Taiwan regions) with the aim having the ability to extend the prediction system to other countries affected by flood threats.

Data
The process of establishing the pre-disaster prediction system is data-driven.We primarily used two sorts of data: first, meteorological data; and second, socio-economic data.To predict the direct economic losses in terms of geographical information and timescale, while considering the integrity, continuity, and variety of data, we selected a time-series dataset of 31 provinces and municipalities in China composed of a period of 13 years, from 2006 to 2018.With the aim of ensuring the timeliness of the prediction system, we chose variables based on significance and whether they can be retrieved pre-disaster or at the beginning of a flooding episode.Therefore, 23 independent variables and a dependent variable, the latter being direct economic loss, were used in the study, as shown in Table 1.There were 403 pieces of data in total showing the meteorological and socio-economic conditions.To predict the indirect economic losses, which are based on the direct economic losses of the agricultural sector, we further utilized the direct economic loss data as the input of the IO model.With the combination of other socio-economic data, we generated the indirect economic losses caused by the demand reduction in the agricultural sector.

Meteorological Data
Meteorological data were obtained from "Daily meteorological dataset of basic meteorological elements of the China National Surface Weather Station (V3.0)" including China's national basic weather stations, reference climatological station, and general weather stations, with a total of 2474 stations of the China Meteorological Administration.All of the research data were retrieved at an annual provincial level.The data included daily maximum precipitation, precipitation anomaly percentage, precipitation anomaly percentage in spring, precipitation anomaly percentage in summer, precipitation anomaly percentage in autumn, precipitation days, moderate rainy days, heavy rainy days, torrential rainy days, and maximum continuous precipitation days.

Socio-Economic Data
To predict the direct losses, the socio-economic data included: casualties, death toll, direct economic losses, sown area with 10% reduced production, sown area with 30% reduced production, sown area with 80% reduced production, railway disruption, road disruption, number of reservoirs, capacity of reservoirs, area with flood prevention measures, area of soil erosion under control, and length of city sewage pipes.The above data were retrieved from the Bulletin of flood and drought disasters in China (2006China ( -2018))

Data Processing
Step 1 processing the missing data: As mentioned in the previous section, all data used in the study were collected online and offline from official datasets, official reports, and reference books.For some missing data, if the data were proved to be below the statistical standard, we substituted it with 0. For the data where the true value was missing, the data were substituted with the average figure in order to mitigate the impact on the prediction model.
Step 2 normalization: Each sample of the original dataset had 23 features (independent variables) to reflect the flood disasters from a specific aspect or information related to the local province.If components are with different magnitudes, they will not make equal contributions to fit the model and are likely to cause a bias.Due to this, we used the normalization process of Min-Max normalization to make the indicators comparable.By doing so, all features were transformed into the range [0, 1], meaning that the minimum value of a feature becomes 0, while the maximum value of that becomes 1.The formula is as follows: where x denotes the input of independent variables.

Gradient Boosting Regression Trees (GBR)
A new machine learning method, Gradient Boosting Regression Trees (GBRT or GBR) was adopted in our research to efficiently predict the direct economic losses of floods.It is a modification of the gradient boosting (GB) algorithm and Classification and Regression Trees (CARTs) by using a regression tree of fixed size as weak learners [23].
The gradient boosting algorithm, proposed by Friedman in 1999 [24], integrates various machine learning and statistical methods, such as gradient algorithm, boosting algorithm and tree algorithm.The rationale of the gradient boosting algorithm is to create new base learners which are maximally correlated with a negative gradient of the loss function.Compared with traditional boosting algorithms, along the direction of the gradient, every new model is built with the aim of reducing the loss function of the previous model.By using the gradient boosting algorithm during the training process, regression was achieved by continuously reducing the residuals [25].The process of gradient boosting with multiple iterations is shown in Figure 1.To elaborate on the training process, in each iteration, a weak classifier is generated by the model, which goes through further training according to the residuals of the previous classifier.Based on the performance of a classifier, the weight is generated.The poorer the performance, the more weight will be given.Finally, the ensemble model is achieved by the weighted summarization of all weak classifiers.With a variance reduction, the prediction accuracy of a classifier is improved.Both continuous and discrete values can be dealt with gradient boosting, and this algorithm can mitigate the drawbacks of overfitting.In the process, a base classifier, also known as a weak classifier, is generally CARTs, a series of decision tree regression models.The weight is generated based on the performance of weak classifiers, and the process of learning consists of inputting the residuals from the previous iteration as the object function of the next base classifier.Classification and regression trees (CARTs), proposed by Breiman et al. in 1984 [26], can be used for both classification and regression models [27][28][29].The two types of trees used in these two models are decision trees and regression trees.Compared with other artificial intelligence models, CARTs have better performance in terms of prediction as they can obtain nonlinear relationships without requiring prior information about the probability distribution of variables.
The GBR algorithm combines weak learners by iteratively concentrating on the errors resulting at each step until a sum of the successive weak learners can a create suitable strong learner.
Given the Dataset  ,  (i.e., historical flooding dataset), the loss function is used to evaluate a set of weights to seek a minimized error with the aim of optimizing the model.Let Xi denote a set of explanatory variables (i.e., disaster-inducing factors, disasteraffected bodies and the disaster-prevention capabilities) and Y be dependent variables (i.e., direct economic losses).There are six main steps of the GBR, which can be expressed as follows: 1. Initialize the parameters in the learning machine with the following equation: where   is the parameter set. is the parametric variable that minimizes the loss function, and L  ,  is the square error loss function in our model.The negative gradient of the loss function is used in the current model as an approximation of the residual.The calculation process of the residual is shown as follows: For iterations m = 1 to M: where   refers to the objective function.
2. A regression tree is generated with J leaf nodes, described as follows: where  refers to a regression tree with J leaf nodes, TNT refers to terminal node tree.
3. Estimating the value of the leaf nodes in the regression tree.The value can be estimated by the following equation: where  refers to the value of the leaf nodes in the regression tree.
4. The learning machine of this iteration can be obtained, as shown in the following: 5. After iterations, the final regression model can be shown as follows: , where Numerous researchers have developed an index system-based assessment method to gain a better understanding of the relationship between factors and natural disaster loss.Sun et al. classified the indicators into several different categories, including disasterinducing factors, disaster-affected factors and disaster-prevention capabilities [30].In our research, we also designed an index system for independent variables as the input of GBR, as shown in Table 1.To predict the direct economic losses of floods pre-disaster or at the beginning of a flooding episode is to establish a regression model.Since the dataset of flooding, a disaster occurring with a low frequency, is discrete with a small scale and large time span, five machine learning methods are generally used, including Bayesian Ridge, Line Linear, Elastic Net, XGB, and GBR [31].Despite the advantages of the other four methods, we compared the five different models to verify the effectiveness of GBR.The comparative experiment is designed as follows in Figure 2: To further reduce the autocorrelation between the variables, principal component analysis (PCA) is implemented.In this study, we used a randomized truncated singular value decomposition (SVD) by the method Halko et al. 2009 [32].
To compare the validation of the performance of the five models, four indicators are used, including Root-Mean-Square Error (RMSE), Mean Absolute Error (MAE), Explained Variance (EV), and R-squared score(R 2 ).The formulas are as follows: )

Var y y EV y y
Var y ) where i y is the actual measurement, ˆi y is the predicted value In this article, Anaconda Navigator 2.0.3, python 3.8.8,Jupyter notebook 6.3.0 were implemented.

Input-Output (IO) Model
Natural disasters, especially floods, have detrimental impacts on the agricultural industry and disrupt public transportation.The agricultural industry and the transportation industry are closely related to other industries, such as real estate, construction, warehousing and retail, accommodation, and catering, etc. Floods also cause indirect losses in upstream and downstream industries that are not directly related, such as the financial, mining, and other industries.Therefore, losses of directly affected industries were predicted by using the machine-learning prediction system, and then as the input, the predicted direct losses were used to evaluate the indirect losses of other industries using the Input-Output (IO) model.
The IO model has been implemented to assess the economic effect of natural disasters since the 1970s.Results have shown that the model can assess related economic losses effectively.In this article, indirect economic losses among the industries are evaluated by using a static IO model [18].
The correlations among the industries in the IO table can be expressed as: That is, where  is the direct consumption coefficient,  is the total output of sector i, and  is the final demand for sector i.The above formula can then be transformed as: where  is the identity matrix and   is the inverse matrix of Leontief.Taking the sectional direct economic losses as losses in final products,   ,  , . . .,  .The total product loss is then: where  denotes the total economic losses, and  denotes the direct economic losses Thus, the loss of indirect input is expressed by the reduction of intermediate input, as  .
To improve the accuracy of the indirect loss assessment of various departments, this paper uses the complete consumption coefficient for analysis.Let B be a complete consumption coefficient matrix obtained by transforming the direct consumption coefficient matrix A, then B =   .Therefore, the total loss of the product can be further expressed as: Assume that  is the economic loss in the sector i caused by floods and the final use of other sectors has no change.The total output of the entire economic system then becomes: where  ,  1,2, … ,  are the complete consumption coefficients.The total output loss of sector i is then: where  is the direct economic loss of sector i and   is the indirect economic loss of sector i.The total product losses of other sectors are:

The Pre-Disaster Prediction System
The overall prediction system process of the GBR-IO model is shown, as follows, in Figure 3.

Direct Economic Loss Prediction
To analyze which features make the most contribution in causing direct economic losses, we firstly retrieved feature importance indicating the usefulness and value of each feature in the construction of the boosted regression trees within the model.The more an attribute is used to make key decisions with decision trees, the higher the F score is.From Figure 4, below, it can be seen that reservoir loss, sown area with 10% reduced production, road disruption, annual anomaly percentage and daily maximum precipitation are the top five features of importance.Further, we developed a heat map of the correlation coefficient between the top 10 features and the model.The heat map is used to directly illustrate the correlation between the indicators and economic losses of flooding.The color shade of the heat map shows the degree of the correlation; that is, the greater the correlation, the lighter the color.The correlation between the indicators and direct economic losses, ranging between −1 to 1, shows whether they are positively correlated or negatively correlated.From the heat map, shown in Figure 5, it is notable that reservoir loss (with the correlation of above 0.8), sown area with 10% reduced production (correlation of 0.6-0.8),casualties (correlation of 0.6), road disruption, number of reservoirs, annual anomaly percentage, summer anomaly percentage and daily maximum precipitation are positively correlated with direct economic losses.In contrast, length of city sewage pipes (correlation of −0.2-0), reservoir capacity (correlation of −0.3) and annual rainy days are negatively correlated (correlation of below −0.2).From the results, it can be inferred that: (1) the agricultural industry is the most directly affected sector; (2) there is a time lag in constructing reservoirs for the mitigation of economic impacts caused by floods; (3) since a region with high annual rainy days is less likely to suffer from economic losses of floods, a possible explanation could be that such regions have a balanced precipitation trend and are unlikely to suffer from floods caused by short-term torrential or heavy rains; (4) the improvement of disaster prevention capabilities is effective in reducing economic impacts; and (5) reservoir loss contributes the most to direct economic losses, among other factors.To further reduce the autocorrelation between variables, principal component analysis (PCA) was implemented.PCA is used to decompose a multivariate dataset in a set of successive orthogonal components that explain a maximum amount of variance.The PCA screen plot in Figure 6 shows the explained variance of all of the variables used in the model.Based on the ranking of variance that is explained, we selected the top 16 principal components as new features, i.e., F1 to F16 in Figure 6.Such new features retained 90.75% of the original information.Due to the characteristics of decision tree algorithms, they do not address the value of variables.Therefore, there is no need for feature normalization and PCA.The tree models (GBR and XGB) in the study were not processed with PCA.The above new features are input into the prediction system to be trained with five different machine learning algorithms.During the training process, K-fold cross-validation (k=5 in this prediction system) is performed.Firstly, the original data set is split into five sets.We then use the four of the five sets as the training set to train the model without repetition.Such an approach can obtain a substantial amount of information from limited data.The cross-validation results are shown in Figure 7.Among the five machine learning methods, Bayesian Ridges and GBR performed better than the other three methods, with a cross score of around 0.7.To further determine which method is better, four regression model evaluation metrics, including RMSE, MAE, EV, and R 2 , were used to compare the prediction effect of the model.RMSE, MAE, EV, and R 2 were calculated for the gradient boosting regression technique, as well as the other four machine learning methods.The lower value of RMSE and MAE and the higher value of EV indicate the more accurate prediction result, and the higher value of R 2 indicates a heightened match between the analytical and predicted values.The advancement of gradient boosting regression, compared to other approaches, is shown in Table 2. Observing the regression metrics, it is notable that, although Bayesian Ridge, Linear Regression, and Elastic Net perform better than GBR in some error metrics, the prediction of a specific dataset can substantially deviate from the true value, indicating overfitting.GBR has the highest R 2 and low RMSE, MAE and high EV in comparison to the predicted and the actual value, indicating that GBR has the best fit of values and a better fit than the other four methods.Therefore, combining the cross-validation results and regression metrics, we obtained GBR as the best machine learning regression model in predicting the direct economic losses of floods.The learning curve is generally used in machine learning to evaluate both the performance of the training and the validation of datasets to diagnose whether the model is underfit, overfit, or well-fit [33].Training samples on the horizontal axis show the size of datasets used in the learning process, while the score on the vertical axis represents R 2 , which is used to evaluate the overall performance.The learning curve of the GBR model is shown in Figure 8.With the expansion of the training set, the cross-validation score approaches a desirable level, and it converges with the training score.It can be seen from the figure that when the training data size reaches about 80 pieces, the over-fitting and under-fitting of the model have been significantly reduced, proving that the GBR model is well-fit.Adding more training data can decrease the variance and bias of the GBR model.When the data size reaches about 150 pieces, the variance and deviation of the prediction results can be further mitigated.Therefore, when the model is generalized to analyze similar problems in other countries or regions, the required data size must reach at least 80 pieces, and the prediction accuracy will be improved with the expansion of data size.Integrated learning is considered to have high predictive precision, especially in terms of those algorithms which use a decision tree as the base learner.To validate GBR as the optimal combined model, we compared the regression result with Bayesian Ridge, Linear Regression, Elastic Net, XGB, and GBR, as is shown in Figure 9.It can be seen that GBR is the best fit for the true value.To verify the effectiveness of the prediction model, we did descriptive analysis and empirical analysis based on the dataset in China from 2006 to 2018.
According to the descriptive statistical analysis for the years 2006 to 2018, the direct economic losses caused by floods in China were overwhelmingly higher in the years 2010 and 2016 than of that in other years and showed an overall increasing trend, as shown in Figure 10.Among the floods in China over this period, floods in 2010 caused by torrential rainfall in southern China were the most devastating to occur since 1998, leading to the worst economic conditions in the past half century [34].Global warming had caused a surge of moisture in the atmosphere.With the combination of a strong El Niño effect and climate change, severe flooding occurred along the Yangtze River in the summer of 2016 [35].The floods that broke out in these two years were the worst in the recent decades resulting in the direct economic losses being much more severe than in the other years analyzed.As for the empirical analysis, by splitting the samples into training sets (data from 2006 to 2017) and the validation set (data in 2018), we gained the regression result of direct economic loss with GBR, as shown in Table 3.The comparison between the true value and predicted value is shown in Figure 11.It is notable that the deviation between the predicted direct economic losses and the true economic losses is relatively negligible.From Table 3, we can see that among 31 provinces in China, Heilongjiang, Shandong, Inner Mongolia, and Gansu provinces, which are heavily reliant on agriculture, suffered the most, with the proportions of the direct economic losses to provincial GDP at 1.73‰, 0.94‰, 0.86‰, and 0.75‰ respectively.From the charts above, it can be seen that, for provinces suffering mild losses from natural disasters, the GBR model can provide relatively accurate prediction results in terms of the direct economic losses, while for the losses caused by extreme weather, despite the suddenness and uncertainty, the prediction results can still partially reflect the characteristics; reflecting the robustness and effectiveness of the model.It is also notable that the provincial direct economic losses caused by floods are geographically unevenly distributed.Gansu, Guangdong, Shandong, Inner Mongolia, and Yunnan suffered more direct economic losses, while Ningxia, Shanxi, and Tianjin suffered less.Such distribution can be attributed to the following reasons: (1) coastal regions are more susceptible to floods caused by typhoons and monsoon climates; (2) provinces with a high reliance on agriculture, such as Henan, Shandong, and Gansu, are more likely to be economically affected by floods as the agriculture sector is directly related to natural disasters; (3) provinces with more flood prevention facilities and infrastructure, such as Beijing, Tianjin, and Shanghai have strong resistance to floods.

Indirect Economic Loss Prediction
In this paper, the IO model was introduced to predict the indirect economic losses caused by floods in 42 sectors of 31 provinces in China.The precondition of the IO model is to apply the economic loss data of directly affected sectors.Since agriculture is the most directly affected sector among all of the 42 sectors, according to the Chinese National Bureau of Statistics, we took it as the directly affected sector to calculate the indirect losses of all of the sectors by using the industry linkage generated from the IO table [36].
For the direct economic loss of the agricultural sector, according to the descriptive statistical analysis (as shown in Figure 11), most of the floods in Northern China happened in autumn, while those in Southern China happened in summer.The sown areas in both northern China in autumn and southern China in summer were used for planting rice.According to China's Ministry of Agriculture and Rural Affairs of the People's Republic of China, the price of rice and the average yield per hectare are published in each season and can be denoted with Pi.To verify the model validation, we obtained the official average price in 2018, i.e., 2.56 RMB/kg, and the average yield per hectare, which was 2347.6667kg/hectare.Thus, the average production per hectare is 2.56 * 2347.6667= 6010.0267RMB/hectare.Since the national standard for affected sown areas includes the sown areas with 10%, 30% and 80% reduced production, according to the bulletin of flood and drought disasters in China (2006-2019), we did a scenario analysis to reflect an average scenario of affected production with the assumption that 40% of the production was reduced due to flood disasters [37].In other extreme situations, related calculations can be also conducted according to the specific situation.The estimated direct economic losses of the agriculture sector in 31 provinces in 2018 is shown, as follows, in Figure 12.Further, to estimate the indirect economic losses of floods of other industries related to the agricultural industry, we analyzed the IO model.Although the 42 sectors are integrated and related, in the Input-Output table, they are split into different sectors without overlapping.Considering the analysis method of the input-output relationship from direct and indirect economic losses, we chose the direct economic loss of the directly affected sectors as the final product loss, where ΔY and ΔX − ΔY are defined the same as above.The complete consumption coefficient was introduced to obtain an accurate assessment of the indirect input loss of various sectors.The complete consumption coefficient is the quantity of the products in sector i that need to be consumed directly and indirectly to produce final products per unit of sector j.As defined above, B was used to represent the complete consumption coefficient matrix.Accordingly, B =    is the relationship between the complete consumption coefficient and the direct consumption coefficient.Thus, B is obtained by the direct consumption coefficient matrix of each one given by the 42 sectors in the IO table of nationwide provinces in 2017.In this paper, we investigated the most directly hit industry-the agricultural sector-which has an indirect effect on all industries.Here we only demonstrated the top six most affected industries' complete consumption coefficients of agriculture out of 42 industries, as shown in Table 4.The figure indicates the correlation between other sectors and the agricultural sector, which can be treated as the input for the following calculations.Based on the above matrix, taking the computation of indirect economic loss of agriculture, forestry, animal husbandry and fishery to other sectors, the total indirect economic loss of the 31 provinces was calculated, as shown in Figure 13.It can be seen that Shandong, Henan, Heilongjiang, and Inner Mongolia, all of which are provinces that are reliant on agriculture, suffered the most indirect economic losses among all provinces, with the predicted indirect economic losses reaching 6235.55 million RMB, 1770.49 million RMB, 2227.45 million RMB, and 1389.13 million RMB, respectively.Thus, to prevent unnecessary economic losses and secondary disasters, these provinces should be given more attention.

Conclusions
This paper proposes an effective prediction model consisting of a cutting-edge machine learning regression model and a traditional economic model, i.e., the GBR-IO model.The performance of the GBR-IO model used in the paper is superior to other models because it improves the prediction effectiveness, reduces the issue of overfitting and combines inter-disciplinary methods.The prediction model obtains outstanding results of direct economic losses on datasets in 31 provinces (excluding Hong Kong, Macau, and Taiwan) in China from 2006 to 2018, with a goodness-of-fit of 90%.Based on the predicted direct economic losses, we obtained the indirect economic losses by using the Input-Output model.Compared with previous studies, the GBR-IO model can predict regional direct economic losses and indirect economic losses pre-disaster or at the beginning of a flooding episode with effectiveness and efficiency.Further, the GBR-IO model has a high generalization ability, which can be applied to other countries, especially to those with small datasets at the minimum requirement of 80 pieces.
However, there is still room for improvement of the pre-disaster prediction model.First, as flood forecasting systems and the numerical simulation technology advance in meteorology, the variables used in the model could be estimated more accurately before a disaster happens.Second, the limitation of the datasets results in a failure to categorize floods.If floods are categorized into river floods, flash floods, and drainage problem floods, the accuracy of the model can be further improved.Third, since geographical conditions are not included in the system, global climate models (GCMs) and regional climate models (RCMs) can be further integrated into the system to boost performance.With the above improvements, the model performance can be further enhanced.
According to the prediction model, it can be concluded that positively correlated indicators include reservoir loss, sown area with 10% reduced production, casualties, road disruption, annual anomaly percentage, and daily maximum precipitation while the length of city sewage pipes, reservoir capacity and annual rainy days are negatively correlated with economic losses.From the empirical case based on datasets in China, we can also conclude that from the proportion of direct economic losses to provincial GDP, it is notable that the impact of flooding on the economy is relatively significant.Among the 31 provinces in China, especially Gansu, Inner Mongolia, Shandong, and Heilongjiang, the direct economic losses even approached or exceeded one-thousandth of the province's GDP.
Thus, some policy implications can be drawn from the results.Since the GBR-IO model has a strong generalization ability, regional flooding databases are encouraged to be established in order to improve the prediction accuracy with more data.Although there is a time lag of reservoir capacity, giving priority to reservoir construction and city drainage capacity can mitigate the economic impact of flooding.Policymakers could also pay more attention to those regions that are heavily reliant on agriculture as they are more vulnerable to flood disasters.Further, if the accuracy of meteorological forecasts can be improved, more effective measures can be taken in advance.

Figure 2 .
Figure 2. The flowchart of model comparison.

Figure 4 .
Figure 4. Features of importance, measured with the F score.

Figure 5 .
Figure 5. Heat map of the correlation coefficient between features and the model.

Figure 7 .
Figure 7. Cross validation of machine learning models.

Figure 8 .
Figure 8. Learning curve of GBR model.

Figure 9 .
Figure 9.The fitting and prediction result of compared models.

Figure 11 .
Figure 11.Predicted direct economic loss versus true value in 2018.

Figure 12 .
Figure 12.Direct economic losses of agriculture sector in 2018.
, published by the Ministry of Water Resources of the People's Republic of China, China statistical Yearbook (2006-2018).To predict the indirect economic losses, the Input and Output table with 42 sectors of 31 provinces from the National Bureau of Statistics (2017) was used to calculate the industry linkage [22].

Table 1 .
Floods disasters direct economic loss prediction index system.

Table 3 .
Predicted direct economic losses.

Table 4 .
Matrix of cumulative input coefficients of agriculture, forestry, animal husbandry and fishery sectors to other sectors.