1. Introduction
The extent of flood damage caused by locally concentrated heavy rainfall has been increasing owing to extreme weather events such as heavy rainfall and snowfall, which have become more frequent because of climate change. The risk of flood damage is also rising because of increased runoff due to continuing urbanization and limited pervious spaces for absorbing and holding stormwater. Flood damage reduction measures for urban areas can be categorized into structural and non-structural methods. Typical structural measures include the maintenance of streams and embankments and maintenance and extension of drainage facilities including pipelines.
Non-structural measures are means of actively preparing for and responding to anticipated weather events by methods such as decision-making systems, for example, flash flood forecasting and warning systems [
1].
Given that most riverside areas in urban watersheds are used as roads or residential areas, it is easier to apply non-structural measures such as flash flood forecasting than short-term structural measures. In particular, Korea is characterized by great differences between its maximum and minimum precipitation amounts, making it prone to floods and droughts. The precipitation in summer (June to September) accounts for approximately 70% of the annual precipitation and heavy rainfall events during this period cause increasingly frequent flood/inundation damage. It is essential to predict such events accurately and to provide detailed forecasting. That is, a decision support system based on accurate prediction and forecasting is necessary to minimize the damage caused by locally concentrated and sporadic heavy rainfall events.
Most previous researches have focused on investigating the causes of floods using rainfall-runoff models. In some studies, physical models have been developed for flood forecasting by considering regional characteristics or specific parameters. These investigations involved rainfall-runoff simulations taking into account climate change [
2,
3]. Attempts have been made to provide results with high reliability for decision making by accurately predicting the changes in runoff intensity and frequency [
4,
5]. Zahmatkesh et al. [
6] proposed a rainfall-runoff model based on K-means clustering to minimize the uncertainties of the parameters constituting the model and applied it to the Bronx River watershed, New York City. Using this model, the runoff based on the predicted rainfall was simulated to predict the future runoff for the watershed. The runoff simulation using the clustering scheme was found to be up to 50% more accurate than the simulations performed using individual models. Zhang and Yang [
7] used genetic algorithms to derive a unified value from eight different rainfall-runoff models. They applied the multi-model ensemble method to the Yellow River Basin and demonstrated that it outperformed individual rainfall-runoff models in terms of prediction reliability. Mehr and Kahya [
8] and Mehr and Nourani [
9] also proposed hybrid models using genetic algorithms on rainfall-runoff modeling results. In all of these studies, the extents of the parameter uncertainties were assessed at each stage of the rainfall-runoff analysis and their impacts on the runoff analysis were evaluated.
The current research trend is to integrate new techniques into rainfall-runoff models. For example, artificial neural network (ANN) modeling has been used for flood forecasting based on the water levels of rivers/streams [
10,
11]. Wu et al. [
12] presented an ANN-based method to predict the variations in watershed runoff and stream water level and proposed a regression equation for water level prediction. Aichouri et al. [
13] compared various ANN-based models for rainfall-runoff analysis. They pointed out that the models compared in the study mostly dealt with linear (regular) rainfall-runoff patterns and highlighted the need to consider nonlinear (irregular) rainfall-runoff. They additionally proposed an ANN-based model designed to overcome this limitation and listed the improved features of the developed model for irregular rainfall-runoff prediction based on the results obtained by applying the model to catchments in the Mediterranean region. Zhang et al. [
14] used an Elman neural network to improve the prediction accuracy of the time-series runoff input data by taking into account the nonlinear characteristics of the data. They proposed a water level prediction model for warning purposes by considering the hidden layer for lead time (t − 1) along with the originally predicted time (t). Studies on ANN-based water level prediction focus on improving the reliability of the predicted values based on the input data (input layer) selection. Stated differently, there is a lack of research on the relative structural suitability of the input, hidden, and output layers and a comprehensive model based on an optimal structure is required. Therefore, flash flood forecasting based on simple prediction using ANN results alone has room for improvement.
With the improved computational capabilities resulting from advanced computer technology, studies on hybrid models combining geographic information system (GIS) technology and distributed runoff models using spatially distributed data have recently received increasing interest [
15,
16,
17]. Liu et al. [
18] focused on runoff behavior (responses) during heavy rainfall events.
They noted that runoff behavior plays an important role in runoff routing and presented the results of diffusive transport using a GIS-based model by applying a grid mesh to the Attert catchment in Luxembourg. For model verification, they compared the results of runoff routing from each cell with the data collected over the previous 30 months, which yielded similar hydrographs for the predicted and observed values. Khalfallah and Saidi [
19] predicted runoff probabilities according to GIS-based rainfall frequency mapping and presented the flood simulation results obtained by applying the prediction results to the Hydrologic Engineering Center’s River Analysis System (HEC-RAS). Cohen et al. [
20] regarded the longitudinal river slope as a key parameter that determines flood discharge variability. They proposed a methodology for calculating the river slope using the Global River-Slope dataset and used the modified river slope for modeling. They also presented regression analysis results depicting the correlations between the catchment parameters and river slope, which are necessary for flood discharge estimation. Prediction using a distributed model is very important because it takes into account the non-uniformity of the geomorphological spatial distribution. However, this prediction process is time consuming because its rainfall-runoff simulation requires hydrodynamic numerical analysis. Furthermore, studies on flood forecasting and warning systems using distributed models are limited to small sub-watershed areas or low resolution.
The UK and US provide user-centered customized flood forecasting based on flood warning, preparedness, and response systems in an effort to reduce flood damage. The UK provides flood risk and flood zone maps covering all riverine areas across the country and operates a three-step flood warning system [
21]. In the US, flood warnings are issued based on 75% of the threshold water level using the streamflow observations received from over 6800 stream gauges in the Automated Flood Warning System network distributed across the country [
22].
As described above, flood risk in urban watersheds is a topic of extensive research, and numerous studies have been conducted on flood warning and response systems. However, it remains necessary to improve the knowledge and application of the prediction methods. Most importantly, the results obtained from various studies should support flood prediction with adequate lead time. Locally concentrated, sudden, excessive rainfall causes rapid watershed response, which makes city centers vulnerable to flood damage. Therefore, flash flood forecasting is very important for ensuring sufficient lead time for disaster prevention.
The primary goal of flash flood forecasting is to prevent disasters that cause loss of human life and property. This objective can be achieved only by rapid and accurate transmission of reliable flood forecasting and warning information. Speed is a very important aspect of flood forecasting and warning. An adequate lead time ensuring timely response directly impacts the application of the flood forecasting and warning system. Models that combine parametric optimization, distributed models, and hydro-geomorphological features are inherently restricted from providing sufficient lead time for flood forecasting and warning. Research methods based on rainfall-runoff analysis are conclusively considered to be the most suitable for ensuring sufficient lead time and are highly effective in terms of practical aspects based on speed and accuracy.
The main purpose of this study was to develop a flood forecasting and warning system that supports rapid decision making. Specifically, the objective was to provide reliable information necessary for flood forecasting and warning by estimating flood and water levels with high accuracy based on a one-dimensional flood analysis model. The high-water marks of a stream were used to develop a stable flash flood warning operation plan. High-water marks are customarily used independently to collect local rain observation records or establish criteria for situation assessment. However, when used independently, high-water marks cannot function as reliable measures for flood warning because they can only help in assessing the current situation or indicate disaster situations that occurred in the past. Furthermore, damage by heavy rainfall should be reduced by providing forecasts with sufficient lead time to ensure appropriate preparedness and response by considering the spatial movement and temporal evolution of rainfall, which conventional water level and rain gauge networks cannot provide.
In summary, it is necessary to establish a rainfall and flood prediction system that considers the upstream and downstream relationships between the currently used water level points to ensure sufficient lead time. Considering the growing need to establish a flash flood forecasting and warning system based on a reliable prediction system, the objective of this study was to develop a flood forecasting and warning method capable of integrated consideration of the past and current states of the upstream and downstream points.
Most urban watersheds are significantly influenced by the backwater conditions of their mainstream channels. Therefore, the accuracy of flood forecasting and warning for urban streams can be ensured only through meticulous analysis of the relationships between rainfall and runoff and between flow rate and water level along with their specific manifestations. In this study, we prepared an analysis system for flood forecasting and warning, which can be directly applied to heavy rainfall events, by establishing a method for hydrodynamic analysis of urban streams exposed to such backwater effects and by creating a database for rainfall rates depending on rainfall duration. The objective of the model is to enable flood forecasting and warning based on the threshold water level at a specified warning point along the stream.
In this research, a flood-prone area located in the typical urban area of Seoul, Korea, which is at risk of flood damage, was selected as the study area. The aim was to support decision makers in operating flood forecasting and warning systems by analyzing the conditions according to the temporal distributions of various rainfall types based on rainfall predictions and allowing sufficient lead time for evacuation. The ultimate goal is to minimize the loss of human life and property in flood events.
4. Predicted Rainfall Verification and Validity Analysis
This section describes the process of checking the results presented in
Section 3 by comparison with actual water level observations to test the validity of the flash flood warning issuance criteria proposed in this report. In general, urban watershed runoff is characterized by sensitivity to rainfall characteristics. Furthermore, the runoff characteristics that undergo various changes under the influence of the backwater of the main stream at the confluence point affect the selection of the flash flood warning time and warning point. Therefore, the optimal time for flash flood warning issuance should be identified by applying a consistent correlation analysis method throughout the analysis process, from selecting the warning point to obtaining the observed rainfall and water level curves of the given stream from the minimum to the maximum duration. In other words, the identification and estimation of the time for flash flood warning issuance should be determined by considering the efficiency of the decision support related to the flash flood warning issuance customized for all rainfall conditions in line with the purpose of constructing a flood forecasting and warning system.
4.1. Identification of Optimal Flash Flood Warning Time
We selected the optimal flash flood warning time by segmenting the water level and rainfall data at regular intervals based on the corresponding real-time observations of the Dorim Stream collected at the Sindorim Bridge point. The collected real-time observations were divided into four segments at intervals of at least 10 min (10, 20, 40, and 60 min), depending on the data characteristics. The maximum interval of 60 min was applied to consider the relatively rapid watershed response to runoff generation according to the actual characteristics of an urban stream.
To test the validity of the flash flood warning time for each duration interval considering the actual water level and rainfall data observed in the watershed, we determined the flash flood warning time using the coefficient of determination (
R2). In general, the coefficient of determination is used when analyzing the degree of correlation between two variables (
x,
y). In this study, the correlations between two variables associated with various rainfall types were tested based on the observed rainfall curve representing the rainfall during each interval and the corresponding observed water level curve. Accordingly,
xn, and
yn in Equation (1) denote the rainfall and water level, respectively, during each observed interval, and
x and
y denote their respective mean observed values. The coefficients of determination for the duration intervals were calculated using the standard deviations of each item,
xi and
yi, and the number of observed data points
n:
As test rainfall events, we used 10 datasets based on the rainfall event observations from July to September 2013 for validation. The rainfall amounts estimated using the proposed flash flood warning method were verified against four rainfall events that exceeded the threshold for issuing a flash flood warning (caution level) out of the 10 real-time observations.
Table 8 describes the four rainfall event datasets used for comparison.
We calculated the coefficients of determination based on the observed rainfall and water level data by segmenting the duration at intervals between 10 and 60 min to select the duration yielding high values for both sets of coefficients of determination as the optimal flash flood warning time. We chose to use the coefficient of determination because of the wide variation in the standard deviation of the unit rainfall applied for each duration. Given that the behavior diverged from that of the actual observed hydrographs, we considered it necessary to analyze the behavior best reflecting the tendency in the context of efficient operation of flash flood warning.
Table 9 presents the calculation results, revealing that the accuracy of the calculated coefficients of determination increases by up to 70% compared with the existing coefficients of determination as the interval increases from 10 min to 60 min. In effect, in contrast to the 26% average reliability with a 10 min interval, over 70% reliability was achieved with the 60 min interval, except for the rainfall event on 8 July 2013, which yielded up to 90.3% reliability, demonstrating that the predicted values well reflect the water level fluctuations. Consequently, we chose 60 min as the optimal interval for flash flood issuance.
4.2. Comparison with Real-Time Water Level Observations and Improvements
An urban watershed, which responds sensitively to runoff trends depending on the rainfall characteristics, is affected by the backwater of the mainstream channel, which affects the selection of the flash flood criteria. To estimate the rainfall amount depending on the optimal water level prediction, the discharge outlet opening and closing times should be determined by applying a consistent correlation analysis method throughout the entire process for flash flood warning, from selecting the warning point to obtaining the observed rainfall and water level curves of the given stream from the minimum to maximum duration.
We selected the optimal prediction time in the regression model using the coefficient of determination (R2) based on the observed water level curve corresponding to the observed rainfall curve. Considering that the water level observation interval of the Dorim Stream is 10 min, we segmented the duration into four intervals (10, 20, 40, and 60 min) for prediction time estimation.
In the estimation results, the highest correlation coefficient was obtained for the longest duration (60 min), which was set as the duration for validity analysis. For the rainfall amount, however, we considered the periods before and after 30 min from the flash flood warning issuance time (t = 0) to ensure safety in consideration of the time gap of 20–30 min between the peak rainfall and peak water level. Thus, a set of three time series (t = −40 ± 20 min, t = −30 ± 30 min, and t = −20 ± 40 min) was considered for warning issuance instead of only one. A flash flood warning is triggered for the relevant watershed when the rainfall amount corresponding to one of these time series exceeds the rainfall threshold estimated in this study.
We tested the validity of the proposed flash flood warning issuance criteria by comparing the predicted values with the corresponding observed values using the datasets for all four rainfall events listed in
Table 9 (
Table 10,
Table 11,
Table 12 and
Table 13).
Out of the 10 major rainfall events that occurred from July to September 2013, we selected four to obtain the validity analysis data based on the observation records for rainfall events exceeding 50 mm/h.
Table 10,
Table 11,
Table 12 and
Table 13 summarize the application results of each rainfall event.
The real-time observations were used as follows: ① corresponds to the predicted rainfall amount at 60 min; ②, ③, and ④ are the aforementioned set of three time series for flash flood warning issuance depending on the water level prediction considering the safety aspect (① ≤ ② or ③ or ④ leads to flash flood warning issuance); and ⑤ is the time of actual issuance of a flash flood warning for the Dorim Stream, which is used for validity testing by checking the predicted value against it.
In
Table 8, the real-time observations used for the validity analysis have the following characteristics: ① corresponds to the 60 min rainfall threshold; and ②, ③, and ④are the aforementioned set of three time series for decision making on flash flood warning issuance considering the safety aspect (① ≤ ② or ③ or ④ meets the conditions for a flash flood warning issuance). In fact, the observed water level reached 6.83 m at 11:40, and the 60 min rainfall threshold for flash flood warning issuance according to the corresponding known WS is 13.26 mm. Although this value is 0.72 m lower than the threshold water level for issuing a caution-level flood warning (7.55 m, as presented in
Table 3), the observed water level will reach 7.65 m at 12:40, thus exceeding the 60 min water level threshold. This finding suggests that the flash flood warning should be issued stepwise starting from 11:40 based on the proposed methodology, well ahead of the water level rise to the threshold level. In effect, the rise in water level in the Dorim Stream starting from 11:40 was such that the water level reached the rainfall threshold for flash flood warning issuance pre-calculated for each duration interval at the given point. Moreover, given that observations at 10 min intervals are used stepwise for actual decision making regarding flash flood warning issuance, the methodology proposed in this report is considered suitable for setting the optimal criteria for issuing flash flood warnings and, consequently, capable of supporting decision-making regarding flash flood warning issuance.
Looking at the four heavy rainfall events used for validity analysis, encroachment occurred as per the 60 min rainfall predicted in each of the cases ②, ③, and ④, based on the Dorim Stream levee encroachment water level, 7.55 m. Furthermore, only one of the four heavy rainfall events showed a time difference of 10 min between the predicted and observed levee encroachment times, with the remaining three heavy rainfall events triggering the pumping operation due to encroachment. Finally, based on the comparison of the predicted rainfall amounts with the corresponding observations at 10 min intervals, it can be concluded that stepwise application of the proposed methodology is valid for actual decision making on flash flood warning issuance.
5. Conclusions
To establish an efficient water damage prevention plan for flood damage minimization in urban watersheds, it is essential to prepare an integrated flash flood forecasting and warning system considering the meteorological, hydraulic, and hydrologic characteristics of streams. The implementation of an integrated flash flood forecasting and warning system must be preceded by the development of a system customized to the user demands, which vary depending on the local and disaster characteristics. In view of most flood-induced damage cases occurring in urban watersheds, there is a compelling need to provide appropriate response strategies.
Most studies conducted and systems established in relation to flood forecasting and warning have been focused on large streams and rivers. However, recent years have seen increased flooding events in small streams and related flood damage. In particular, various cultural and leisure spaces created along urban streams are prone to inundation damage, in addition to overflow onto stream embankments. Despite the urgency of countering such flood damage, there are hardly any systematic analysis results regarding urban streams and only a few flash flood forecasting and warning systems are in place.
We established an analysis process for appropriate flash flood forecasting and warning for an urban stream and presented a method for deriving the rainfall thresholds for issuing flash flood warnings. To this end, we developed a rainfall threshold estimation approach enabling short-term prediction based on the real-time stream water level at a flood warning point, as well as the backwater effect of the main stream at the confluence point, and performed analysis to predict the optimal flood warning time depending on the changes in the known downstream WS conditions in the Dorim Stream. We also highlighted the need to consider the effects of water level fluctuations based on the observed upstream and downstream water levels to ensure efficient operation of a flash flood warning system. Based on the results derived, we aim to develop analysis techniques and use them as elementary methods for studying the adequacy of flood forecasting and warning system operation based on reliable validity testing because such techniques are ultimately determined by the criteria for ensuring efficient resident evacuation. The flash flood warning issuance criteria for urban streams proposed in this report are expected to facilitate the issuance of adequate flash flood warnings and setting up of related response strategies considering the duration-dependent rainfall amounts estimated based on scenarios for various backwater conditions. This work is also expected to contribute to the establishment of evacuation-related disaster prevention strategies. However, the short-term rainfall amount estimation proposed in this study is based on the rainfall threshold for estimating the time for flash flood warning issuance for the given stream and is not linked to analysis using radar-based rainfall information and the like. Therefore, it is necessary to develop an integrated model, such as a short-term prediction system, in follow-up research. Furthermore, the most important aspect in flood warning-related research and the selection of warning issuance time is the collection and storage of high-water mark points and high-quality observation data. These elements are still lacking, and continuous observations and validity testing as well as correction of forecasting and warning models are thought necessary. In this regard, the flash flood warning technique and analysis results presented in this report will prove their value only through continuous observation and validation.