Next Article in Journal
Sugar Beet Processing Wastewater Treatment by Microalgae through Biosorption
Next Article in Special Issue
Generating Continuous Rainfall Time Series with High Temporal Resolution by Using a Stochastic Rainfall Generator with a Copula and Modified Huff Rainfall Curves
Previous Article in Journal
Multi–Proxy Reconstruction of Drought Variability in China during the Past Two Millennia
Previous Article in Special Issue
Using a Self-Organizing Map to Explore Local Weather Features for Smart Urban Agriculture in Northern Taiwan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Threshold Rainfall in Ungauged Areas Using Machine Learning

Department of Urban and Environmental and Disaster Management, Graduate School of Disaster Prevention, Kangwon National University, Samcheok 25913, Korea
*
Author to whom correspondence should be addressed.
Water 2022, 14(6), 859; https://doi.org/10.3390/w14060859
Submission received: 22 January 2022 / Revised: 4 March 2022 / Accepted: 6 March 2022 / Published: 10 March 2022

Abstract

:
In recent years, Korea has seen abnormal changes in precipitation and temperature driven by climate change. These changes highlight the increased risks of climate disasters and rainfall damage. Even with weather forecasts providing quantitative rainfall estimates, it is still difficult to estimate the damage caused by rainfall. Damaged by rainfalls differently for inch watershed, but there is a limit to the analysis coherent to the characteristic factors of the inch watershed. It is time-consuming to analyze rainfall and runoff using hydrological models every time it rains. Therefore, in fact, many analyses rely on simple rainfall data, and in coastal basins, hydrological analysis and physical model analysis are often difficult. To address the issue in this study, watershed characteristic factors such as drainage area ( A ), mean drainage elevation ( H ), mean drainage slope ( S ), drainage density ( D ), runoff curve number ( C N ), watershed parameter ( L p ), and form factor ( R s ) etc. and hydrologic factors were collected and calculated as independent variables, and the threshold rainfall calculated by the Ministry of Land, Infrastructure and Transport (MOLIT) was calculated as a dependent variable and used in the machine learning technique. As for machine learning techniques, this study uses the support vector machine method (SVM), the random forest method, and eXtreme Gradient Boosting (XGBoost). As a result, XGBoost showed good results in performance evaluation with RMSE 20, MAE 14, and RMSLE 0.28, and the threshold rainfall of the ungauged watersheds was calculated using the XGBoost technique and verified through past rainfall events and damage cases. As a result of the verification, it was confirmed that there were cases of damage in the basin where the threshold rainfall was low. If the application results of this study are used, it is judged that it is possible to accurately predict flooding-induced rainfall by calculating the threshold rainfall in the ungauged watersheds where rainfall-outflow analysis is difficult, and through this result, it is possible to prepare for areas vulnerable to flooding.

1. Introduction

Climate change has increased rainfall in Korea, resulting in various natural disasters that cause rapidly increasing social and economic loss [1]. However, Korean weather forecasts only provide rainfall information in absolute terms, and the same heavy rain warnings and special reports apply to all areas in Korea, which means a failure to reflect regional differences. For this reason, even with accurate forecasts, the forecast system fails to provide specific information on how different areas are affected and damaged by weather events. Forecasts focused on physical aspects of weather events do not provide sufficient information on how people’s properties and safety are affected by them.
It is for this reason that the World Meteorological Organization (WMO) emphasizes the need for ‘impact forecasts’ that consider the socioeconomic effects that may be caused by weather events [2]. In Korea, different organizations provide different definitions of impact forecasting. However, they can be summarized as follows: forecast that scientifically estimates the socioeconomic impact of weather at different times and places and delivers the estimates along with detailed weather information [3,4]. Outside of Korea, the WMO defines impact forecast as forecast that provides information on expected risks along with weather forecasts when disaster-causing high impact weather is expected. According to the Met Office of the United Kingdom, it is defined as a forecast that estimates the socioeconomic impact of a climate disaster at the time and place of its occurrence by considering meteorological disasters, level of exposure to disasters, and regional vulnerabilities. The National Weather Service of the United States defines it as a service aimed at providing the people with information on the social, economic, and environmental impact of weather, hydrological, and climate events [5]. Leading countries in the field of meteorology already provide information on socioeconomic impact of weather events along with high-resolution weather information. In the United Kingdom, the Flood Forecasting Centre (FFC) provides Flood Guidance Statements (FGS) that assess the risks of all flood types over five days and publish the findings daily [6,7]. The FFC uses the information to publish a table of flood risks which divides flood impact into four stages.
Impact forecasting requires threshold rainfall. Threshold rainfall means the rainfall amount that causes inundation. Accurate impact assessment requires calculation of the precise inundation-causing rainfall in each area. However, in Korea, research on threshold rainfall has been lacking. Most researchers use simplified analysis methods rather than refined hydraulic and hydrological analyses. Hydrological analyses of coastal areas are too complex to conduct properly.
As for previous literature on threshold rainfall calculation, ref [8] developed a flash food monitoring and prediction (F2MAP) model to calculate the flash flood-threshold runoff from rainfall. Ref [9] analyzed the relationship between flash flood index and runoff number characteristics to develop an equation between the two. Ref [10] proposed a threshold runoff calculation method using the flash flood guidance (FFG) model, which is more suitable for Korea rather than those used in the United States. The researchers presented the method as a way to acquire basic data for a flash flood forecast system. Ref [11] analyzed runoff in Jeju using the SWAT-K model that combines DEM, landcover, soil map methods, and developed a threshold runoff simulation method (TRSM) specifically for the island. Ref [12] used ArcGIS and HEC-GEOHMS to divide the Nakdonggang River watershed into 2268 sectors, drew rainfall-peak flow curves for different initial loss scenarios and antecedent moisture conditions, and calculated the threshold rainfalls. Ref [13] estimated threshold rainfalls for different durations using events with damage caused by past rainfall in urban areas and others without such damage. Ref [14] stressed the need for impact forecast and estimated threshold rainfalls using the SWMM model. Ref [7] also linked the grid base inundation analysis model (GIAM) for grid-based inundation analysis. Using the Huff distribution [14], the researcher converted the data into time-series rainfall data to simulate inundation depths, and inversely estimated the threshold rainfall based on the inundation depths. Ref [15] collected data on rainfall and typhoon damage over the last five years where inundation was caused, analyzed the relationship between rainfall and the damages, and developed an equation for threshold rainfall ( y = a x b ). As can be seen from the literature cited above, Korean studies on threshold rainfall mostly used hydrological models. Few researchers studied threshold rainfall by considering hydrological characteristics.
More recently, a number of researchers used machine learning to improve the accuracy of threshold rainfall analysis [16,17,18]. Additionally growing is the body of literature that study rainfall-runoff, rainfall damage, and flood estimation with machine learning and deep learning rather than hydrological models [19,20,21,22,23]. However, few studies were identified in Korea that used machine learning to calculate threshold rainfall. Ref [19] sought to predict river water levels using observation data and deep learning algorithms. To that end, the researchers used tensor flow to predict water levels at the Okcheon Observatory location along the upper stream sectionof the Daecheong Dam within the Geumgang watershed and used TensorFlow to develop a multiple regression model and a long short-term memory (LSTM) artificial neural network model. Ref [24] used three machine learning techniques (support vector machine, decision-making tree, and random forest) to develop a function for predicting rainfall damage in the Seoul Metropolitan Area (SMA) and found that support vector machine analysis using meteorological observation data from two days before yields the highest prediction performance. Ref [25] used the machine learning method on Gyeonggi-do, the province that suffers the worst rainfall damage each year. Choi used the data on rainfall damage of public facilities from the 2006–2015 Disaster Yearbooks published by the Ministry of the Interior and Safety (MIST) as the dependent variable. Ref [26] used machine learning methods such as ESN and DeepESN to predict rainfall using rainfall, pressure, and humidity from 2004 to 2014 as mediating variables. The correlation factors calculated using DeepESN yielded better results. Ref [27] performed hydrological rainfall adjustment using Light GBM and XGBoost. They found clear adjustment effects across all rainfall events after Light GBM and XGBoost learning, despite the fact that rainfall is adjusted 5 to 20 mm less.
Much of the literature cited above only used a single hydrological model. However, in this study, two models were coupled and used to calculate the marginal rainfall [8] and it is considered an advantage of this paper to apply the results to machine learning. Figure 1 shows the flow chart of this study. It was analyzed through machine learning using threshold rainfall and topographic factors of standard watershed units. The threshold rainfall was used as a dependent variable, and the topographic factor of the standard watershed unit was designated as an independent variable. In addition, the model with the smallest error was selected using error performance analysis to calculate the threshold rainfall in the ungauged basin where hydrological analysis was difficult. The ungauged basin means a coastal area where it is difficult to calculate the threshold rainfall.

2. Theoretical Background

2.1. Definition of Threshold Rainfall

In this study, threshold rainfall is calculated using the method used by the MOLIT in 2007, according to which threshold rainfall means the three-hour rainfall causing inundation depth at which the flow overflows the river embankment [8]. Threshold rainfall can be calculated by determining the rainfall of the rainfall-runoff curve corresponding to the threshold rainfall. In general, the runoff calculation equation for a rainfall-runoff model can be expressed as follows [15,28].
R t = R i + R p
where R t is the total runoff, R i is the runoff at the impermeable layer, and R p is the flow at the permeable area.
R t =   FFG × I + f ( FFG ) × ( I 1 )
In a rainfall-runoff model, rainfall and soil moisture constitute the inputs. However, the opposite is true with the flash flood threshold; calculation of flash flood threshold requires current soil moisture and required flow as inputs. As such, the equation on upper stream water and small-sized rivers is converted for FFG using the repetitive calculation method, as shown in Figure 2, to calculate the rainfall that causes threshold runoff. The FFG is the rainfall corresponding to the threshold runoff in the relationship of the rainfall and runoff curve. If there is no impervious area, the relationship between R and FFG can be expressed in Figure 2 and Equation (2). R means the threshold runoff (mm), FFG means the flash flow guidance (mm), and f ( ) means the fall-runoff process. Moreover, I means rainfall intensity.

2.2. Machine Learning Method

Machine learning is an area of artificial intelligence where numerical models, algorithms, and programs are used to have a machine learn from given data as humans do, and new information is derived, or decisions are made based on what it learns [29]. In other words, machine learning means a system that uses accumulated empirical data to build models and improve performance. The amount of data matters in machine learning, and higher-quality data leads to higher-performing results. As for machine learning methods, this study used random forest, support vector machine, and XGBoost.
(1) Random Forest
The random forest method uses bootstraps to create several samples and applies them to a decision tree model to compile the results [30]. A decision-making trees produces estimates by creating and learning one-time training data from a given dataset. On the other hand, a random forest creates multiple training data from a given data set and creates and combines multiple decision-making threes for improved prediction [31]. The observations not used by individual decision-making trees are out-of-bagging (OOB) data and are used for estimating prediction probability and identifying variables. The prediction probability of OOB observations for each observation k within the x i category (0 or 10) [32].
This study used Python and the random forest method to calculate threshold rainfalls. Figure 3 shows the conceptual diagram of a random forest.
p k ^ ( x i ) = j O O B i I y ^ ( x i , t j ) = k | O O B i |   ,   for   k = 0 ,   1
(2) Support Vector Machine
Support vector machine (SVM) is a supervised learning algorithm used for both linear and non-linear classification issues. The purpose of SVM algorithms is to determine the lines or boundaries dividing an n-dimensional space into separate groups, so that they can be classified as their proper categories when new data are given. There may be multiple lines or boundaries for dividing an n-dimensional space into classes. However, the optimal boundary should be identified to determine categories. This optimal determination boundary is called the hyperplane. A support vector is the vector closest to the hyperplane and affects its position. The support vector machine is an algorithm that determines the optimal hyperplane that maximizes the margin, which means the distance between different data points.
The support vector regression (SVR) model has a small number of support vectors, and thus is known to be less sensitive to outliers. Ref [33] developed support vector regression that adopts a ε -insensitive loss function into the support vector machine. Support vector regression is estimated using the function shown in Equation (4) [32].
f ( x ) = ω t x + b
Equation (5) shows the constraints for calculating the optimal hyperplane function while calculating the error that minimizes Equation (4). 1 2 ω 2 describes the degree of flattening of the function. If the data cannot be completely linearly separated, a slack variable ξ ( i = 1 , I ) is introduced to process it. ξ means the distance between the margin and the data outside the boundary between the margins. The main superparameters of the support vector regression are C (cost) and γ , and C adjusts the complexity of the estimation model and the degree of error tolerance. An increase in C means imposing a high penalty on errors within the margin. ϵ is not considered in the calculation process if the error is less than ϵ due to the maximum deviation between the actual value and the estimated value.
In this study, SVM was converted into SVR to predict arbitrary real values and used, and a Gaussian kernel (RBF) known for its excellent performance was applied. Figure 4 shows the conceptual diagram of the support vector machine.
Minimuze : 1 2 w 2 + C i = 1 I ( ξ i + ξ i * ) Subject   to : y i w x i b ϵ + ξ i   w x i + b y i ϵ + ξ i * ξ i ξ i * 0
(3) eXtreme Gradient Boosting
Similarly to random forests, XGBoost is an ensemble algorithm that addresses the errors of multiple decision trees [35]. It offers improved prediction performance over gradient boosting machine (GBM) through distribution and parallel processing. In general, it is ten times faster than GBM. The efficiency and scalability of this method has been validated in multiple previous studies [36,37]. This boosting method lowers errors by grouping multiple classification and regression trees (CARTs).
y ^ i = k = 1 k f k ( x i ) , f k F
Equation (6) shows an ensemble model of trees, where K is the number of trees and F represents the set of CARTs. f k corresponds to the weight of each independent tree and leaf. The scores of the leaves are summed up and compared for final prediction.
O b j = i n = l ( y i , y i ^ ) + k = 1 K Ω ( f k )
Equation (7) represents an XGBoost model. The first l ( y i , y i ^ ) is a loss function that represents the difference between a prediction and an actual observation. The second Ω ( f k ) is the normalization term that controls the complexity of the model to prevent overfitting.

2.3. Performance Assessment Using K-Fold Cross Validation

This study uses MAE, RMSE, and RMSLE to compare the performance of different models. Most studies use the above three indicators a lot for data comparison [38,39,40]. They are widely used to objectively assess the accuracy of a regression equation by analyzing differences between observations and estimates. MAE and RMSE are statistical indicators for confirming the degree of errors included in an estimate calculated using an equation, when compared with an observation. A value closer to 0 represents better fit. RMSLE represents the average ratio of observations to predictions.
MAE = 1 N i = 1 N | y i y i ^ |
RMSE = 1 n i = 1 n ( y i y i ^ ) 2
RMSLE = 1 n i = 1 n ( log ( y i + 1 ) log ( y i ^ + 1 ) ) 2
Ideally, these errors need to be tested by applying them to actual ungauged watersheds. However, due to data and time constraint, the prediction models were validated using five-fold cross validation. K-fold cross validation is a model assessment method that uses a part of the overall data as a validation set. It ensures that all data are used as dataset at least once. Figure 5 shows dividing the data into five datasets and validating the models with a different dataset each time. An average cross-validation uses five datasets. This study selects the optimal parameters following cross validation to calculate threshold rainfalls.

3. Selection of Target Watersheds and Variables

3.1. Selection of Target Watersheds

This study chose the Han River watershed as its target, as the area includes the highest number of standard watersheds according to the water resource unit map. 290 of Korea’s 850 standard watersheds are included in the Han River watershed. 237 of the 290 watersheds are inland, and the other 53 are coastal watersheds, as shown in Figure 6.
The (a) section of Figure 7 shows the learning watersheds; 80% of the learning watersheds were used for machine learning, and the other 20% were used for validation. High-performing models were selected with (a), and predictions were performed for the watersheds highlighted yellow in (b). The data set of the basin used for machine learning was randomly selected and proceeded.

3.2. Dependent and Independent Variables

This study used the threshold rainfalls calculated using the MOLIT method [8,15,41] as dependent variables. Figure 8 shows the calculated threshold rainfalls on the map.
Characteristic factors of the watersheds were used as independent variables. The analysis only considered topographical factors and hydrological factors. The watershed characteristic factors used in were collected from the Water Resources Management Information System (www.wamis.go.kr, accessed on 31 December 2011) and the geographic information system (GIS). Data on 15 characteristic factors were collected, including: drainage area ( km 2 ), mean drainage elevation ( m ), mean drainage slope (%), highest drainage elevation ( m ), drainage density, runoff curve number, river length ( km ), drainage perimeter ( km ), form factor, circularity ratio, stream frequency, channel maintenance constant, relative relief, number of reliefs, and river length ratio.
The drain area refers to the area on the plane of the basin and refers to the plane area within the closed curve, which is usually made up of a watershed. The basin average elevation is calculated by arithmetically averaging the elevation values corresponding to each cell of the DEM (Digital Elevation Model). The mean drain slope is calculated by arithmetically averaging the slope corresponding to each cell of the DEM in degrees. Highest drain evaluation means the highest elevation in the basin, and drain density means the length of rivers per unit area. It means that the degree of outflow of the basin is quantified by the Soil Conservation Service (SCS) using the runoff curve number land use and soil map. River length is the total length of all rivers in a given drainage basin. The drain perimeter is defined as the length measured along the boundary of the watershed of a given order projected on the horizontal plane of the map, and the form factor is defined as the ratio of the main river length of the watershed to the diameter of the circle having the same area as the watershed area. The circularity ratio is a dimensionless parameter defined as the ratio of the basin area to the area of a circle with the same length as the basin circumference. Stream frequency is defined as the ratio of river water in the basin to the basin area, and the channel maintenance constant is the reciprocal of the aqueous density. Relative relief is defined as the ratio of watershed undulations to watershed circumference, number of reliefs is defined as the product of watershed undulations and water density, and river length ratio is defined as the ratio of river length w to average river length w − 1. Table 1 shows a summary of the watershed characteristics factor.
A correlation analysis was performed to select statistically correlated independent variables, as independent variables not correlated to dependent variables may lower the prediction performance. The correlation analysis was performed as shown in Figure 9. Given the fact that the threshold runoffs required for calculating threshold rainfalls were calculated from peak flood volumes and overflowing runoffs, the following variables were determined to be significantly correlated: drainage area, river length, drainage perimeter, relative relief, and river length ratio. Among those factors, river length, river length ratio, relative relief, and drainage perimeter were determined to be more highly correlated with the independent variables. The correlation coefficient was 0.65 for threshold rainfall and drainage area, 0.64 for drainage perimeter, and 0.31 for river length. As such, drainage area, drainage perimeter, and river length were finally selected as independent variables. In most data analyses, principal component analysis (PCA) should be used to reduce initial independent variables [42], but in this study, a principal component analysis was omitted because the amount of data for each independent variable was not large.

4. Machine Learning Application and Results

This study used SVM, random forest, and XGBoost. Excel 2010 and Python ver. 3.6 were used to record and statistically analyze the collected data and generate graphs. This study also used the model packages provided by Python-based Scikit-learn.
Effective machine learning requires pre-processing of the data to be used. The independent variable data went through data scaling and missing values were removed. As data scales vary depending on the variable, the data were standardized to render them more suitable for machine learning. Independent variables were analyzed using RobustScaler, which is less affected by outliers. A higher accuracy can be expected by removing outliers. However, the small number of inputs in this study means possible overfitting. Therefore, this study addressed outliers through pre-processing rather than outlier removal.

4.1. Validation of Prediction Models

The optimal parameters for each model were selected through K-fold cross-validation. An error closer to 0 indicates a better result. Table 2 shows the MAE, RMSLE, and RMSLE values of each of the five datasets created by dividing the datasets through the k-fold cross-validation. All three performance assessments found that XGBoost produces the results closest to actual observations compared with the other models.
Parameters were applied to increase the accuracy of machine learning, and n_estimator represents variables that adjust the number of trees to generate. max_depth means the number of tree depths. min_samples_split represents the minimum number of sample data to split nodes, and min_samples_leaf means the minimum number of sample data required for a leaf node. learning_rate means a parameter that, in machine learning and statistics, moves toward the min loss function and determines the size of each stage of repetition. The calculated parameter values area n_estimators: 100, learning_rate = 0.04, min_samples_leaf = 3, min_samples_split = 2, max_depth = 4.
Figure 10 compares the existing threshold rainfalls with those calculated using XGBoost. Most threshold rainfalls are distributed between 40 and 60 mm and between 60 and 80 mm and are close to actual observations.
Figure 11 is a map representing threshold rainfall values. The watersheds with low threshold rainfalls in (a) are reflected in (b) as well.

4.2. Calculation of Threshold Rainfalls in Ungauged Watersheds

This study used XGBoost, which produced good results in error performance assessment, to calculate the threshold rainfalls of ungauged watersheds. Figure 12 shows the distribution of the threshold rainfalls calculated for the ungauged basins. The majority of watersheds show threshold rainfalls between 40 mm and 80 mm. Figure 13 is a map showing the threshold rainfalls of the ungauged basins other than the those in the inland areas.

4.3. Validation Using Real World Cases and Assessment

As shown in Figure 14, based on the rainfall events in 2017, 2020, and 2021 of each affected watershed, among the ungauged watersheds outside the purple lines around the inland part of the Han River, Yongin, Cheonan, Samcheok, Gangneung, and Sokcho watersheds were found to be vulnerable against heavy rain. An application to actual rainfall events showed that damage was caused when the rainfall exceeds the specified rainfall in the legends. However, the researchers’ ability to verify damages in other areas was restricted by the fact that damage was verified from news articles and social network posts.

5. Conclusions

Damage caused by localized heavy rain continues to increase in Korea. However, research on inundation-causing threshold rainfall is still largely absent in the country. More researchers need to study technologies for predicting and responding to inundation in advance. This study can be summarized as follows.
The purpose of this study was to identify threshold rainfalls in areas not readily available for hydrological analysis, using the calculation method and characteristic factors used by the MOLIT. Three machine learning methods (SVM, random forest, and XGBoost) were compared in terms of accuracy using MAE, RMSE, and RMSLE, and XGBoost was selected as the best-performing method. Watershed characteristics, hydrological factors, and XGBoost were used to calculate the threshold rainfalls of the ungauged coastal watersheds. In this study, it is judged that what can reflect actual topographic and hydrological factors can be differentiated from other machine learning and marginal rainfall papers. In addition, distinct from conventional simple data, data using physical models were used in machine learning techniques, so high accuracy could be secured through a small number of data, and anyone could use it by using widely known machine learning techniques.
However, this study has its limitations. First, outliers were found while calculating the threshold runoffs of the hydrological models. More sophisticated hydrological models and more accurate data may be needed for analysis. In addition, threshold rainfall calculation based on the runoff-rainfall curve simply used polynomials. However, higher accuracy may be achieved by applying a machine learning method to threshold runoff and runoff-rainfall curve calculation.
This study compared the calculated threshold rainfalls with real world cases identified from news reports and social network posts, which was found to pose limitations to quantitative assessment.
The researcher plans to conduct a similar study nation-wide. Watershed data with more diverse hydrological models and outliers will improve the accuracy of the findings. Although not included in this study, quantitative validation using real world events will yield meaningful results. The current water forecast system provides only quantitative figures without considering the damage caused, which some regard as insufficient for supporting effective decision-making to prevent and prepare for damage caused by natural disasters. Threshold rainfall prediction suggested in this study may, if implemented on a continued basis, provide accurate information on rainfall damage in advance and help decisionmakers make better decisions on disaster control.

Author Contributions

K.-S.C. and C.-H.O., J.-R.C. carried out the survey of previous study and wrote the graph of the data. B.-S.K. suggested idea of study and contributed to the writing of the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Korea Meteorological Administration Research and Development Program under Grant KMI [2021—00312]. This work was financially supported by Ministry of the Interior and Safety as Human Resource Development Project in Disaster Management(C2001777-01-01).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, J.Y.; Choi, C.H.; Kang, D.S.; Kim, B.S.; Kim, T.W. Estimating Design Floods at Ungauged Watersheds in South Korea Using Machine Learning Models. Water 2020, 12, 3022. [Google Scholar] [CrossRef]
  2. WMO. Guideline on Multi-Hazard Impact-Based Forecast and Warning System; WMO-No. 1150; World Meteorological Organization: Geneva, Switzerland, 2015; p. 22. [Google Scholar]
  3. Jeong, K.Y. Impact Forecast Vision and Direction. Meteorol. Technol. Policy 2016, 9, 6–22. [Google Scholar]
  4. Go, Y.H. Development and Activation of Impact Forecasting Service. Available online: https://www.kma.go.kr/down/t_policy/t_policy_201706.pdf (accessed on 2 June 2021).
  5. Lee, D.H.; Kim, G.Y. Weather Disaster Impact Forecast. KISTEP Technol. Trend Brief 2019, 10, p1–p4. [Google Scholar]
  6. FFC (Flood Forecasting Centre). Flood Guidance Statement User Guide (Version 4); Flood Forecasting Centre: Exeter, UK, 2017; p. 6. [Google Scholar]
  7. Lee, B.J. Analysis on Inundation Characteristics for Flood Impact Forecasting in Gangnam Drainage Basin. Korean Meteorol. Soc. 2017, 27, 189–197. [Google Scholar]
  8. Kim, B.S.; Hong, J.B.; Kim, H.S.; Yoon, S.Y. Development of Flash Flood model Using Digital Terrain Analysis Model and Rainfall RADAR: I. Methodology and Model Development. J. Korean Soc. Civ. Eng. B 2007, 27, 151–159. [Google Scholar]
  9. Kim, B.S.; Kim, H.S. Estimation of the Flash Flood Severity using Runoff hydrograph and Flash flood index. J. Korea Water Resour. Assoc. 2008, 41, 185–196. [Google Scholar] [CrossRef] [Green Version]
  10. Cho, B.G.; Ji, H.S.; Bae, D.H. Development and Evaluation of Computational Method for Korean Threshold Runoff. J. Korea Water Resour. Assoc. 2011, 44, 875–887. [Google Scholar] [CrossRef] [Green Version]
  11. Chung, I.M.; Lee, C.W.; Kim, J.T.; Na, H.N.; Kim, N.W. Development of Threshold Runoff Simulation Method for Runoff Analysis of Jeju Island. J. Environ. Sci. Int. 2011, 20, 1347–1355. [Google Scholar] [CrossRef]
  12. Park, J.B.; Kang, D.G.; Park, M.J. Rainfall Thresholds Estimation to Develop Flood Forecasting and Warning System for Nakdong Small River Basins. J. Korean Soc. Hazard Mitig. 2013, 13, 311–317. [Google Scholar] [CrossRef] [Green Version]
  13. Kang, H.S.; Choi, C.W.; Lee, H.S.; Hwang, J.G.; Moon, H.J. Development of an ANN-Based Urban Flood Alert Criteria Prediction Model and the Impact of Training Data Augmentation. J. Korean Soc. Hazard Mitig. 2021, 21, 257–264. [Google Scholar] [CrossRef]
  14. Huff, F.A. Time Distribution of Rainfall in Heavy Storms. Water Resour. Res. 1967, 3, 1007–1019. [Google Scholar] [CrossRef]
  15. Lee, S.H.; Kang, D.H.; Kim, B.S. A Study on the Method of Calculating the Threshold Rainfall for Rainfall Impact Forecasting. J. Korean Soc. Hazard Mitig. 2018, 18, 93–102. [Google Scholar] [CrossRef] [Green Version]
  16. Zagorecki, A.T.; Johnson, D.E.; Ristvej, J. Data mining and machine learning in the context of disaster and crisis management. Int. J. Emerg. Manag. 2013, 9, 351–365. [Google Scholar] [CrossRef]
  17. Chang, F.-J.; Guo, S. Advances in hydrologic forecasts and water resources management. Water 2020, 12, 1819. [Google Scholar] [CrossRef]
  18. Park, J.; Park, J.-H.; Choi, J.-S.; Joo, J.C.; Park, K.; Yoon, H.C.; Park, C.Y.; Lee, W.H.; Heo, T.-Y. Ensemble Model Development for the Prediction of a Disaster Index in Water Treatment Systems. Water 2020, 12, 3195. [Google Scholar] [CrossRef]
  19. Jung, S.H.; Lee, D.E.; Lee, K.S. Prediction of River Water Level Using Deep-Learning Open Library. J. Korean Soc. Hazard Mitig. 2018, 18, 1–11. [Google Scholar] [CrossRef]
  20. Kao, I.F.; Zhou, Y.; Chang, L.C.; Chang, F.J. Exploring a Long Short-Term Memory based Encoder-Decoder framework for multi-step-ahead flood forecasting. J. Hydrol. 2020, 583, 124631. [Google Scholar] [CrossRef]
  21. Chang, L.C.; Chang, F.J.; Yang, S.N.; Tsai, F.H.; Chang, T.H.; Herricks, E.E. Self-organizing maps of typhoon tracks allow for flood forecasts up to two days in advance. Nat. Commun. 2020, 11, 1983. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Chang, F.J.; Tsai, M.J. A nonlinear spatio-temporal lumping of radar rainfall for modeling multi-step-ahead inflow forecasts by data-driven techniques. J. Hydrol. 2016, 535, 256–269. [Google Scholar] [CrossRef]
  23. Zhou, Y.; Guo, S.; Chang, F.J. Explore an evolutionary recurrent ANFIS for modelling multi-step-ahead flood forecasts. J. Hydrol. 2019, 570, 343–355. [Google Scholar] [CrossRef]
  24. Choi, C.H.; Kim, J.S.; Kim, D.H.; Lee, J.H.; Kim, D.H.; Kim, H.S. Development of Heavy Rain Damage Prediction Functions in the Seoul Capital Area Using Machine Learning Techniques. J. Korean Soc. Hazard Mitig. 2018, 18, 435–447. [Google Scholar] [CrossRef] [Green Version]
  25. Choi, C.H. Development of Combined Heavy Rain Damage Prediction Models Using Machine Learning and Effectiveness of Disaster Prevention Projects. Ph.D. Thesis, Inha University, Incheon, Korea, 2019. [Google Scholar]
  26. Yen, M.H.; Liu, D.W.; Hsin, Y.C.; Lin, C.E.; Chen, C.C. Application of the deep learning for the prediction of rainfall in Southern Taiwan. Sci. Rep. 2019, 9, 12774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Lee, Y.M.; Ko, C.M.; Shin, S.C.; Kim, B.S. The Development of a Rainfall Correction Technique based on Machine Learning for Hydrological Applications. J. Environ. Sci. Int. 2019, 28, 125–135. [Google Scholar] [CrossRef]
  28. Yu, B.I. Study on Flash Floods and Debris Flow Guidance of Gangwon-do. Master’s Thesis, Kangwon University, Chuncheon-si, Korea, 2016. [Google Scholar]
  29. Kim, D.W.; Hong, S.J.; Kim, J.W.; Han, D.G.; Hong, I.P.; Kim, H.S. Water Quality Analysis of Hongcheon River Basin under Climate Change. J. Wetl. Res. 2015, 17, 348–358. [Google Scholar] [CrossRef]
  30. Suh, J.D. Foreign Exchange Rate Forecasting Using the GARCH extended Random Forest Model. J. Ind. Econ. Bus. 2016, 29, 1607–1628. [Google Scholar]
  31. Choi, C.H.; Park, K.H.; Park, H.K.; Lee, M.J.; Kim, J.S.; Kim, H.S. Development of Heavy Rain Damage Prediction Function for Public Facility Using Machine Learning. J. Korean Soc. Hazard Mitig. 2017, 17, 443–450. [Google Scholar] [CrossRef]
  32. Jang, D.R. A Study on the Art Price Prediction Model Using Machine Learning. Master’s Thesis, Hongik University, Seoul, Korea, 2020. [Google Scholar]
  33. Vapnik, V. Support-Vector Networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar]
  34. Kleynhans, T.; Montanaro, M.; Gerace, A.; Kanan, C. Predicting Top-of-Atmosphere Thermal Radiance Using MERRA-2 Atmospheric Data with Deep Learning. Remote Sens. 2017, 9, 1133. [Google Scholar] [CrossRef] [Green Version]
  35. Kim, I.H.; Lee, K.S. Tree based ensemble model for developing and evaluating automated valuation models: The case of Seoul residential apartment. J. Korean Data Inf. Sci. Soc. 2020, 31, 375–389. [Google Scholar] [CrossRef]
  36. Chen, T.; Guestrin, C. XGBoost: A scalable tree boosting system. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016; pp. 785–794. [Google Scholar]
  37. Ke, G.; Meng, Q.; Finley, T.; Wang, T.; Chen, W.; Ma, W.; Ye, Q.; Tie, Y.L. LightGBM: A Highly efficient gradient boosting decision tree. In Proceedings of the 31st Conference on Neural Information Processing Systems, Long Beach, CA, USA, 4–9 December 2017; pp. 3149–3157. [Google Scholar]
  38. Ioannou, K.; Karampatzakis, D.; Amanatidis, P.; Aggelopoulos, V.; Karmiris, I. Low-Cost Automatic Weather Stations in the Internet of Things. Information 2021, 12, 146. [Google Scholar] [CrossRef]
  39. Stefanidis, S.; Dafis, S.; Stathis, D. Evaluation of Regional Climate Models (RCMs) Performance in Simulating Seasonal Precipitation over Mountainous Central Pindus (Greece). Water 2020, 12, 2750. [Google Scholar] [CrossRef]
  40. Ioannou, K.; Myronidis, D.; Lefakis, P.; Stathis, D. The use of artificial neural networks (anns) for the forecast of precipitation levels of lake doirani(N. greece). Fresenius Environ. Bull. 2010, 19, 1921–1927. [Google Scholar]
  41. Choo, K.S.; Kang, D.H.; Kim, B.S. A Study on the Estimation of the Threshold Rainfall in Standard Watershed Units. J. Korean Soc. Disaster Secur. 2021, 14, 1–11. [Google Scholar]
  42. Myronidis, D.; Ivanova, E. Generating Regional Models for Estimating the Peak Flows and Environmental Flows Magnitude for the Bulgarian-Greek Rhodope Mountain Range Torrential Watersheds. Water 2020, 12, 784. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flowchart of study.
Figure 1. Flowchart of study.
Water 14 00859 g001
Figure 2. Concept of flash flood concept.
Figure 2. Concept of flash flood concept.
Water 14 00859 g002
Figure 3. Concept of Random Forest.
Figure 3. Concept of Random Forest.
Water 14 00859 g003
Figure 4. Concept of Support Vector Machine [34].
Figure 4. Concept of Support Vector Machine [34].
Water 14 00859 g004
Figure 5. Concept of cross vaildation.
Figure 5. Concept of cross vaildation.
Water 14 00859 g005
Figure 6. Study Areas.
Figure 6. Study Areas.
Water 14 00859 g006
Figure 7. Training and Prediction Watersheds for Machine Learning.
Figure 7. Training and Prediction Watersheds for Machine Learning.
Water 14 00859 g007
Figure 8. Standard watershed Threshold rainfall.
Figure 8. Standard watershed Threshold rainfall.
Water 14 00859 g008
Figure 9. Variable correlation analysis.
Figure 9. Variable correlation analysis.
Water 14 00859 g009
Figure 10. Threshold rainfall distribution.
Figure 10. Threshold rainfall distribution.
Water 14 00859 g010
Figure 11. Comparison of training data against actual values.
Figure 11. Comparison of training data against actual values.
Water 14 00859 g011
Figure 12. Ungauged watershed threshold rainfall distribution.
Figure 12. Ungauged watershed threshold rainfall distribution.
Water 14 00859 g012
Figure 13. Threshold rainfall values for ungauged watersheds.
Figure 13. Threshold rainfall values for ungauged watersheds.
Water 14 00859 g013
Figure 14. Threshold rainfall values for ungauged watersheds.
Figure 14. Threshold rainfall values for ungauged watersheds.
Water 14 00859 g014
Table 1. Summary of independent variables.
Table 1. Summary of independent variables.
A H S E m H C N L w L p R s R c C f C R p R n R L
Count290290290290290290290290290290290290290290290
mean144.6324.035.4253.31.758.712.967.31.00.42.40.713.313761.8
Max571.6930.365.1302.74.087.963.3262.33.60.712.69.536.83921.44.4
min39.04.94.0103.70.133.70.032.70.00.00.10.30.932.80.7
Table 2. Comparison of model performance evaluation.
Table 2. Comparison of model performance evaluation.
ModelMAERMSERMSLE
Support VectorFold 115190.26
Fold 223380.4
Fold 319260.29
Fold 421260.47
Fold 528380.34
Random ForestFold 112190.28
Fold 220320.46
Fold 316200.34
Fold 422270.47
Fold 521260.45
XGBoostFold 114200.28
Fold 220330.38
Fold 316200.29
Fold 421270.46
Fold 525370.35
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chu, K.-S.; Oh, C.-H.; Choi, J.-R.; Kim, B.-S. Estimation of Threshold Rainfall in Ungauged Areas Using Machine Learning. Water 2022, 14, 859. https://doi.org/10.3390/w14060859

AMA Style

Chu K-S, Oh C-H, Choi J-R, Kim B-S. Estimation of Threshold Rainfall in Ungauged Areas Using Machine Learning. Water. 2022; 14(6):859. https://doi.org/10.3390/w14060859

Chicago/Turabian Style

Chu, Kyung-Su, Cheong-Hyeon Oh, Jung-Ryel Choi, and Byung-Sik Kim. 2022. "Estimation of Threshold Rainfall in Ungauged Areas Using Machine Learning" Water 14, no. 6: 859. https://doi.org/10.3390/w14060859

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop