1. Introduction
Global warming-initiated, extreme weather events receive great attention worldwide. South Korea, in particular, has faced such events, including increasing temperature and rainfall and a growing number of heavy rain days, for the recent 100 years [
1], which has led to natural disasters such as localized heavy rainfall, wind and waves, droughts, and heavy snows. Notably, the summer season from June to September shows a tendency of having an increased number of debris flows [
2]. Debris flows are a type of natural disaster that occurs by a complex interaction between flooding from heavy rainfall and ground soil, as well as by a wide range of other factors such as thawing during spring, indiscriminate logging, and forest fire. They are also, commonly, secondary damage from typhoons and localized heavy rains, with the latter being their main cause because of how heavy rainfall brings an increase in flow speed, soil loss, and large-scale movement of rocks that lead to huge disasters [
3]. In South Korea, damage from debris flows has been reported frequently nationwide, with examples such as Inje County and Pyeongchang County of Gangwon Province in 2006; Seoul, Chuncheon City, and Pocheon City in 2011; Samcheok City in 2012; Busan Metropolitan city in 2014; and Cheongju City and Cheonan City in 2017. For this study, debris flows are seen as mainly from localized heavy rains. In this regard, it requires a thorough understanding of the characteristics of rainfall events that cause debris flows, when establishing an early-warning system for debris flow damage and related planning, maintaining, or managing disaster prevention facilities.
In South Korea, studies on forecasting of debris flows and landslides are mainly about using the related standards provided by the Korea Forest Service and the Korea Meteorological Administration to review their relevance with an analysis of rainfall events that cause debris flows and landslides or to quantitatively calculate the standards. However, studies on debris flow forecasting based on rainfall events have not been actively conducted [
4,
5,
6,
7,
8,
9,
10].
Table 1 and
Table 2 show the forecasting standards for landslides and rainfall, provided by the Korea Forest Service and the Korea Meteorological Administration, respectively. Such standards mainly defined rainfall and accumulated rainfall separately.
South Korea forecasts landslides and debris flows by analyzing rainfall and basin characteristics and using models to calculate the triggering factors. With the current advancement in radar technologies, studies are continuously conducted for forecasting using radar data [
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21]. Therefore, the study attempted to establish a method that considers rainfall intensity and accumulated rainfall not as an independent factor but a function. To this end, it modified the RTI calculation method developed by Jan and Lee [
22] to support the forecasting of debris flows potentially caused by rainfall.
The study used past rainfall data from 80 stations located at the areas that experienced damage from debris flows from 2012 to 2013 for rainfall intensity and accumulated rainfall for each rainfall duration. Based on this, it classified debris flow damages to estimate the rainfall triggering index (RTI). In addition, it calculated the average intensity of the rainfall that causes debris flows. For debris flow forecasting, the study classified the forecasting standards for accumulated rainfall into ALERT (RTI from 10 to 50%), WARNING (RTI from 50 to 70%), and EMERGENCY (RTI from 70% or higher). The 10%, 50%, and 70% RTIs were divided by the average rainfall intensity to estimate the critical accumulated rainfall () and its curve by duration. The calculated was applied to the actual cases of Umyeon Mountain of Seoul, Chuncheon of Gangwon County, and Cheongju City of Chungcheongbuk Province, where damage actually occurred, to make the debris flow forecasting for 24 h of the rainfall triggering such, which aims to determine its applicability for debris flow forecasting.
2. Materials and Methods
To analyze the influence from the interlinkage between accumulated rainfall and rainfall intensity, the study collected the rainfall data of 80 areas that experienced debris flow damage in Gangwon Province from 2012 to 2013 and used the rainfall amount by duration with a maximum of 24 h in which debris flows occurred. Based on this, the RTI, an index for accumulated rainfall and rainfall intensity was calculated for 6 h, 12 h, and 24 h, respectively. Furthermore, the study estimated an average rainfall intensity at the time of debris flow occurrence before using the RTI equation to calculate
for each duration (6 h, 12 h, and 24 h). The
of 10%, 50%, and 70% was then used with the occurrence probability to define the three risk levels. In addition, based on actual damage cases, the study developed a nomogram for continuous rainfall to verify its applicability for debris flow forecasting (
Figure 1).
3. Theoretical Background
3.1. Debris Flow
A debris flow refers to the dynamic phenomenon where soil, rocks, and floating substances flow down a slope by gravity with changes in their shape and sizes. Sharpe [
23] differentiated debris flows from debris avalanches in his United States–based studies, with the former as a movement of soil and rocks saturated with water at a water channel with a steep slope, and the latter as a phenomenon where fragmented soil of an upper layer at a steep slope flow fast, similar to a snow avalanche. As shown in
Figure 2, the path of debris flows comprises three zones: initiation, transportation, and deposition [
24]. Since debris flows have pressure 4–5 times higher than that of flooding water, given that they are mixed with soil and rocks, their external force is 10 times higher than that of flooding water when conflicting with facilities [
25].
Major factors that have influence on the occurrence of debris flows include topographic factors (slope angle, slope impact, and facilities to reduce the flow of pumice stones and soil), geographical factors (depth of soil layers and characteristics of top soil), and hydrological factors (amount of rainfall). Among such factors, rainfall increases pore water pressure and soil weight and leads to erosion and scour of the surface. The analysis of scales and accumulated rainfall indicates that an area with 200 mm or higher of rainfall and 20 mm/h of rainfall intensity will face severe damage with increasing frequency (
Figure 3). This result suggests that areas with low vulnerability may experience a higher probability of debris flow occurrence, in a rainfall episode with a certain level and intensity. Therefore, for the rainfall that triggers a debris flow, it is standard to consider both accumulated rainfall and rainfall intensity observed at the time of its occurrence.
3.2. Estimation of Critical Accumulated Rainfall Using RTI
The RTI model developed by Jan and Lee [
22] was designed to forecast debris flows triggered by rainfall in real time. For the RTI calculation, rainfall intensity (
I) and accumulated rainfall (
) are used as follows.
In the equation above,
I indicates rainfall intensity (mm/h) and
is the accumulated rainfall (mm) observed shortly before the occurrence of debris flows. Among of the rainfall episodes for up to seven days, the one that continues for 24 h with a direct influence on debris flows is considered as antecedent rainfall. The study used the rainfall accumulated for duration 6 h, 12 h, and 24 h to estimate the RTI. Since rainfall has a direct impact on the occurrence of debris flows, especially its accumulation and intensity, the existing system for forecasting landslides uses forecasting for accumulated rainfall and rainfall intensity and daily rainfall, whereas the RTI is calculated with accumulated rainfall and rainfall intensity to consider both the amount and intensity. However, the RTI can be difficult to understand for communities where debris flow–related damage is expected, as it is only a combination of rainfall intensity and accumulated rainfall and does not directly deliver the information about a risk level of debris flow. Therefore, the RTI was converted into critical accumulated rainfall (
) to aid understanding in the provided forecasting. Since the RTI focused on damage in Taiwan during the country’s developed stage, it showed a gap for the rainfall and intensity of South Korea. Therefore, the study changed the level to 10%, 50%, and 70%, taking into consideration the flood forecasting standards provided by the Han River Flood Control Office [
26].
Figure 4 shows the definition of RTI and
.
4. Result and Discussion
4.1. Analysis of Debris Flow-Triggering Rainfall Data
In South Korea, mountainous areas account for 60% of its territory. Since most of them are concentrated in Gangwon Province, debris flow damage is frequently reported for the province. In this regard, the study collected data on the debris flow–triggering rainfall from 80 stations for 2012 to 2013 in Gangwon Province, where debris flows easily occur, and calculated accumulated rainfall and rainfall intensity at the site of damage occurrence (
Table 3).
Figure 5 shows the points of debris flows and the current status of rainfall monitoring stations.
Figure 6 and
Figure 7 show dispersion of the maximum accumulated rainfall and rainfall intensity for 6 h, 12 h, and 24 h at the 80 stations in the damaged areas.
4.2. Development of Nomogram for Debris Flow Forecasting, Using RTI and
The study used rainfall information from the 80 stations mentioned above to calculate the RTIs by phase (ALERT, WARNING, and EMERGENCY) for each rainfall duration (6 h, 12 h, and 24 h). The RTIs were estimated as 600 (10%), 1350 (50%), and 2321 (70%) for the 6 continuous hours; 494 (10%), 1496 (50%), and 1900 (70%) for the 12 hours; and 570 (10%), 950 (50%), 1442 (70%) for the 24 hours.
Table 4 summarizes the calculated RTIs and accumulated rainfall and rainfall intensity for each duration.
Prior to forecasting debris flow, related standards should be established. In South Korea, flood forecasting is made, wherein flood levels are standardized with 50 to 70% of the project flood water levels, in general, applied for the warning and alerting. As explained above, the study referred to the flood forecasting standards of the Han River Flood Control Office [
26], with the following set for each level: 10 to 50% of the occurrence possibility for ALERT, 50 to 70% for WARNING, and 70% or higher for EMERGENCY. Furthermore, the study classified three forecasting levels for the durations of 6 h, 12 h, and 24 h.
Figure 8 shows events of the 80 stations in relation with RTIs, whereas
Figure 9,
Figure 10 and
Figure 11 show graphs of the RTI estimations.
RTIs are not information obtained directly from rainfall. Moreover, most people find RTIs difficult to understand and use. Therefore, the study converted RTIs to
to aid understanding. To estimate values, average rainfall intensity was used for each duration.
that corresponds to the average rainfall intensity is shown in
Table 5.
Figure 12 shows RTIs for 10%, 50%, and 70% calculated from
Figure 9,
Figure 10 and
Figure 11 and
estimation graphs.
The study developed a nomogram for debris flow forecasting by rainfall duration, using the critical accumulated rainfall (
) for each occurrence possibility (10%, 50%, and 70%) and duration (6 h, 12 h, and 24 h). As shown in
Figure 13, a nomogram is a graph of debris flow forecasting levels for the rainfall accumulated from the start to 24 h of the duration. For each duration, the debris flow forecasting levels (ALERT, WARNING, and EMERGENCY) are classified with different colors to aid the visual expression of each level by duration of accumulated rainfall.
4.3. Review on Applicability of Debris Flow Nomogram with Actual Cases
To review applicability of the debris flow nomogram that the study developed, it applied the nomogram to cases of damage caused in the past by debris flows. The representative cases include Umyeon Mountain of Seoul Seoul Metropolitan City, the capital of Korea in 2011, Chuncheon City of Gangwon Province in 2011, and Cheongju City of Chungcheongbuk Province in 2017. The study estimated the response time before the damage occurrence by forecasting debris flows with the actual rainfall data for the cases. The case of Umyeon Mountain where debris flows occurred at 10:00 in 27 July 2011, resulted in 18 deaths and the evacuation of 400 people. In 2011, Chuncheon City of Gangwon Province experienced debris flows that occurred at 24:00 and caused 13 deaths and 26 injuries. The case of Cheongju City of Chungcheongbuk Province occurred at 11:00 on 16 July 2017, causing two deaths.
Figure 14 shows the photos of damaged areas taken at those times.
The results of debris flow forecasting with the nomogram the study developed are as follows.
4.3.1. Case 1: Umyeon Mountain, Seoul
For the case of Umyeon Mountain of Seoul, the capital of Korea, it started raining at 17:00 on 26 July and recorded the maximum accumulated rainfall 307 mm (
Figure 15) until 16:00 on 27 July with damage occurring at 9:00 on 27 July. The debris flow forecasting results were ALERT for 18:00 on 26 July, WARNING for 19:00 of the same day, and EMERGENCY for 5:00 on 27 July (
Figure 16). Based on this, it can be assumed that damage occurs after the EMERGENCY level. Therefore, it is estimated that 4 h of response time is secured prior to damage occurrence. When forecasting is made additionally for the WARNING level, the response time that can be secured is estimated as 7 h.
Regarding the comparison analysis with the alerting standards of the Korea Forest Service and the Korea Meteorological Administration, the former provided the same level of risk; however, it produced the Alarm level for 18:00 of 26 July and 2:00 of 27 July, which are some hours before the damage occurrence, with its response time delayed for an hour. On the other hand, the latter provided the Alarm level for 19:00 of 26 July, which is some hours before the damage occurrence, and issued the alert for 24:00, which is 3 h passed the actual damage occurrence (
Table 6).
4.3.2. Case 2: Chuncheon, Gangwon Province
In the case of Chuncheon City, located in Gangwon Province, the rainfall started at 1:00 of 27 July, and the 230 mm of maximum accumulated rainfall was recorded until 24:00 of the same day (
Figure 17). The damage occurred at 24:00 of 27 July, and forecasting for debris flows was made on 4:00 for ALERT, 19:00 for WARNING, and 21:00 for EMERGENCY (
Figure 18). With the application of the EMERGENCY level, it was found that 4 h of response time was secured prior to the damage occurrence, and with the additional forecasting for the WARNING level, a total of 6 h of the time was secured. Regarding the comparison analysis with the alerting standards of the Korea Forest Service and the Korea Meteorological Administration, the former provided the same risk level; however, it produced the Alarm level for 19:00 of 27 July and for 21:00 of 27 July again before the actual damage occurrence. The standards of the latter issued Alarm from 1:00 of 27 July, which is some hours before the damage occurrence. This is a level lower, compared to the actual risk level of Warning at the time of damage occurrence (
Table 7).
4.3.3. Case 3: Cheongju, Chungcheongbuk Province
In the case of Chuncheon City, located in Chungcheongbuk Province, the rainfall started at 1:00 of 16 July, and the 290 mm of maximum accumulated rainfall was recorded until 14:00 of the same day (
Figure 19). The damage occurred at 11:00 of 16 July, and forecasting for debris flows was made at 8:00 for WARNING and 9:00 for EMERGENCY (
Figure 20). With the application of the EMERGENCY level, it was found that 2 h of response time was secured prior to the damage occurrence, and with the additional forecasting for the WARNING level, a total of 3 h of the time was secured. Regarding the comparison analysis with the alerting standards of the Korea Forest Service and the Korea Meteorological Administration, similar tendency risk levels were shown for all three alerting standards (
Table 8).
5. Conclusions
In this study, we collected rainfall data targeting the areas that experienced damage from debris flows from 2012 to 2013, and developed the debris flow nomogram that reflects both accumulated rainfall and rainfall intensity. It used the two elements observed shortly before the occurrence of debris flows to estimate RTIs and set the three levels according to the possibility of debris flow occurrence: 10 to 50% for ALERT, 50 to 70% for WARNING, and 70% or higher for EMERGENCY. In addition, to help the understanding of the residents in the areas where debris flows can occur, the study converted RTIs to actual accumulated rainfall values () for use in forecasting. In this study, the debris flow nomogram was developed for each duration (6 h, 12 h, and 24 h) and applied to actual cases of debris flow damage for Umyeon Mountain of Seoul, Inje County of Gangwon Province, and Cheongju City of Chungcheongbuk Province.
As a result, the use of the nomogram for debris flow forecasting that the study developed could secure sufficient response time for the cases of Umyeon Mountain of Seoul and Chuncheon of Gangwon Province, where rainfall continues for long durations, and the case of Cheongju of Chungcheongbuk Province where heavy rain is localized. Results for each case are summarized as follows.
Case 1: In the case of Umyeon Mountain of Seoul, 280 mm of the rain that continued for 17 h caused the occurrence of debris flows. The results of using the nomogram in forecasting debris flows for the EMERGENCY level showed that it could secure 4 h of the response time. When the forecasting was made additionally for the WARNING level, a total of 7 h of the response time could be secured to ensure reactive actions.
Case 2: In the case of Chuncheon of Gangwon Province, 260 mm of the rain for about 24 h caused the occurrence of debris flows. The results of using the nomogram in forecasting debris flows showed that it could secure 4 h of the response time. In addition to the forecasting for the WARNING level, a total of 6 h of the response time could be secured.
Case 3: In the case of Cheongju of Chungcheongbuk Province, 290 mm of the rain for about 11 h caused the occurrence of debris flows. The results of using the nomogram in forecasting debris flows showed that it could secure 2 h of the response time. With addition to the forecasting for the WARNING level, a total of 3 h of the response time could be secured.
The results above suggest that the debris flow forecasting nomogram provided by the study is applicable for the actual forecasting on debris flow damages that can be caused by the long-term increase in rainfall and short-term, localized heavy rain. Meanwhile, in the cases of Seoul and Chuncheon, the forecasting standards of the Korean Meteorological Administration and the Korean Forest Service led to the indiscriminate issuance of alerts at the starting point of rainfall. However, the forecasting with the nomogram of the study is expected to support the understanding of rainfall value by general users with a visual representation of the risk level, and allow a proper forecasting or response system to the situation.
Due to the diverse causes of debris flows, rainfall-related factors are not enough in determining debris flow occurrence. Therefore, it is crucial to provide the standards that general people can use to make decisions even without expert knowledge. As rainfall is considered the most common factor that causes debris flows, it is expected that the forecasting on debris flows using the nomogram can support the easier interpretation of general users for debris flows. In addition, the forecasting that uses the nomogram the study developed and radar rainfall information can prevent debris flow damage in real time.