A GIS-Based Framework for Real-Time Debris-Flow Hazard Assessment for Expressways in Korea

Debris flows caused by heavy rainfall in mountain areas near expressways lead to severe social and economic losses and sometimes result in casualties. Therefore, the development of a real-time system for debris-flow hazard assessment is necessary to provide preliminary information for rapid decision making about evacuations or restoration measures, as well as to prevent secondary disasters caused by debris flows. Recently, various map-based approaches have been proposed using multi-attribute criteria and assessment methods for debris-flow susceptibilities. For the macrozonation of debris-flow hazard at a national scale, a simplified method such as the Korea Expressway Corporation (KEC) debris-flow hazard assessment method can be applied for systematic analysis based on geographic information systems (GIS) and monitoring networks. In this study, a GIS-based framework of real-time debris-flow hazard assessment for expressway sections is proposed based on the KEC debris-flow hazard assessment method. First, the KEC-based method was standardized in a systematic fashion using ArcGIS, enabling the objective and quantitative acquisition of various attribute datasets. The quantification of rainfall criteria also was considered. A safety management system for debris-flow hazard was developed based on the GIS platform. Finally, the method was applied and verified on three expressway sections in Korea. The grading standard for each individual influencing attribute was subsequently modified to more accurately assess the debris-flow hazards.


Introduction
Landslides, and in particular debris flows, are one of the most damaging natural hazards in mountainous terrain that experiences heavy torrential rainfall. Rainfall-triggered landslides are a recurring problem in Korea due to a combination of mountainous terrain, a shallow layer of alluvial soil, and associated weather conditions (Hungr et al. 2001;Park et al. 2013). ''Flowage'' or ''flow'' is the term used to describe the downslope movement of unconsolidated material in which particles move about and mix within the mass, such as earthflows, debris flows, or avalanches. To avoid confusion resulting from the interchangeable use of terminology, the term ''debris flow'' in this study refers to the fluid mixture of rocks, sand, mud, and water that is intermediate between a landslide and a water flood (Feusto et al. 1999;Keller and DeVecchio 2008).
Debris-flow damage includes loss of human life, destruction of facilities, and damage to roads, pipelines, and vehicles (Jakob and Hungr 2005). Until now, only post-event repairs and protective infrastructure improvements have generally been executed after debris-flow occurrences. Recently, an increase has occurred in Korea in both the number of occurrences of and the costs for countermeasures against debris flows. In order to sufficiently manage and protect expressway sections and facilities that are exposed to debris-flow occurrences, a method to assess the hazard of debris flows during certain rainfall events at a regional scale is needed. Infiltration of prolonged, intense rainfall, which causes soil saturation and a temporary increase in pore water pressure, is the mechanism by which most shallow landslides, and more specifically debris flows, are generated during rainstorms (Iverson 2000;Keller and DeVecchio 2008). Because Korea is a region with heavy, concentrated, torrential rainfall in the summer, most of the debris flows occur during the rainy summer periods in mountainous regions. Due to climate change and other environmental factors, the intensity and frequency of rainfall events are increasing. This increase results in a greater number of debris-flow occurrences and higher restoration costs than in the past. Consequently, more debris-flow countermeasures are now necessary for the safety and preservation of human lives and infrastructure.
Assessment of landslides, including debris flows, has been carried out with geographic information systems (GIS) techniques combined with statistical analyses and physical-based approaches by various researchers such as Carrara et al. (1999), Dai and Lee (2002), Ohlmacher and Davis (2003), Ayalew et al. (2004), Wang et al. (2008), and Kritikos and Davies (2015). In spite of the importance of GIS approaches, there are few techniques for debris-flow zonation and geodatabase analysis, and there is no practical and standardized debris-flow mapping method using GIS (Paulín et al. 2014). Through comparison and review of existing studies on the factors that influence debris flow and appropriate assessment methods, the Korea Expressway Corporation (KEC) debris-flow hazard assessment method was developed as a fundamental assessment tool (Korea Expressway Corporation 2009). This eclectic method focuses on the probable occurrence of road structure hazards in a highway system. The method can be quantitatively and objectively performed in a simple manner by using documents such as numerical maps and expressway design files. By drawing on existing databases, the method minimizes the need for tiresome field investigations in countless potential debris-flow occurrence regions in vast areas. Above all, in order to immediately assess rainfall-induced debris flows, rainfall criteria were utilized to determine fundamental trigger values and indices for the severity level of debris flows.
According to the KEC method, zonation of debris-flow hazard is evaluated through two indices: susceptibility value and vulnerability value. The susceptibility value indicates the likelihood of debris-flow occurrence, and can be estimated with the topography information available on target locations. The vulnerability value is used to represent whether an occurred debris flow will damage or have an impact on specific expressway sections, and can be assessed from the capacity of the drainage facility and the roadside margin area for sedimentation of debris-flow materials before they reach expressway structures. Determination of these influencing factors is achieved through the use of digital maps and expressway design files. The calculated susceptibility value and vulnerability value are used to indicate a single hazard class (Fig. 1), which represents the likelihood of damage by debris-flow events under a given rainfall intensity.
In this study, the GIS-based framework of real-time debris-flow hazard assessment for expressway sections was proposed based on the KEC debris-flow hazard assessment method. To standardize the KEC method, a systematic sequence for the acquisition of attribute values was proposed using Arcgis 10.2 (Esri 2006). Then to consider the real-time link with the Automatic Weather Station (AWS) network in the Korean Meteorological Administration, 1 the rainfall recurrence periods were quantified as accumulated rainfall using the trigger values of past debris-flow occurrence events. Following the optimization of the KEC debris-flow hazard assessment method, a safety management system for debris-flow hazards was developed based on the ArcGIS platform. This system consists of a database (DB) and four systematic submodules: the input module, the debris-flow hazard assessment module, decision-making module of debris-flow countermeasure, and the real-time prediction module of debris-flow hazard, which is linked to the rainfall monitoring network. Finally, an applicability evaluation of the framework was carried out on three expressway sections of the Korean highway system that have suffered damage due to debris-flow occurrences during heavy rainfall events in 2005 and 2006: the Pyeongchang area of the Yeongdong Expressway, the Deogyu Mountain area of the Daejun-Jinju Expressway, and the Juksan-Geochang area of the 88 Expressway. The reliability of the assessment method was investigated by comparing the actual debris-flow occurrence and nonoccurrence cases.

Literature Review of Debris-Flow Hazard Assessment
A review of methods for debris-flow hazard assessment, including those by Dai and Lee (2002), Lin et al. (2002), Ayalew et al. (2004), Lee and Pradhan (2007), Blahut et al. (2010), Yune et al. (2010), and Lee et al. (2012), permits a comparative analysis of the influential attributes and hazard assessment methods of debris flows. Debris-flow hazards are dependent on a specific set of factors and processes that are usually investigated by various experts (such as hydrologists, geologists, or civil engineers). Prediction concerns where and/or when debris flow will occur, and depends on the type of movement and the scope of the forecasting. Many attributes (geography, rainfall, geology, vegetation, wildfire history, and conditions of existing structures) are related to the mechanism/initiation of debris flows (Table 1). To obtain information on the influencing factors other than topographical properties (elevation, slope, valley shape, and watershed size) and rainfall data, field surveys should be thoroughly conducted throughout entire expressway facility sections. Because the assessment of debris-flow hazards in this study is to be applied on a regional scale for management of the entire national expressway system, the method needs to be simple and applicable to the macrozonation of debris-flow hazards. Other factors, such as the size and shape of a valley along with the variations in slope direction, properties of the subsoil, geological properties, and vegetation cover also have an influence on the movement of debris flows. To simplify the method, only slope information was considered for debris-flow movement possibilities. Because obtaining and processing all the attributes stated above in the prediction stage is a time-consuming and difficult task, we decided that only easily accessible document data (such as digital maps, geological maps, and so on) were to be used in the assessment process. For example, the spatial correlations (for site-specific microscale areas) between influence of debris flow and geological or vegetation characteristics were not well investigated in Korea. Since most debris flows regardless of scale are induced by regional torrential rains, available and reliable influencing attribute data such as rainfall and digital map is necessary to immediately evaluate the debris-flow hazard for a wide coverage area. Consequently, the KEC debris-flow hazard assessment method has become a fundamental assessment tool owing to its simplicity. As a precautionary measure, the KEC method assesses the hazard of debris flows at a regional scale using a limited number of datasets. Only digital elevation models (DEMs) and expressway design files of the target area were used. The DEMs that were used in the assessment process were those provided by the National Geographic Information Institute (NGII) of Korea, whereas the expressway design files were provided by the Korea Expressway Corporation. Through the application and use of only easily accessible datasets, the KEC hazard assessment method minimizes the need for time-consuming field investigation, and allows rapid hazard assessment of debris flows in large area (administrative sections or nationwide). Accordingly, the DEM for elevation and slope was used to evaluate the watershed information containing topographical route condition of rainfall-induced debris flow. Two indices-the susceptibility value and vulnerability value-were used to evaluate debris-flow hazard.
The susceptibility value indicates the likelihood of debris-flow occurrence in a target area and is assessed using a total of four attributes. The mean watershed slope and area percentage of watershed with slopes over 35°are used to assess debris-flow initiation. Mean valley slope and length percentage of valley with slopes over 15°are used for the assessment of debris-flow movement. Other factors such as the size and shape of the valley along with variations in slope direction, properties of the subsoil, and vegetation also have an influence on the initiation and movement of debris flows. To simplify the method, only the slope information derived from the DEMs is considered for debris-flow possibilities. Each of the influencing factors are assigned points from 0 to 5 based on the grading standard set (Table 2) considering past debris-flow occurrence cases for 10 years (from 2000 to 2009), when the annual precipitation rate amounted to 3000 mm and adds up to a maximum total susceptibility value of 20 points.  For the weight considerations of the four attributes, logistic regression was carried out through the Statistical Package for Social Science. The occurrence of debris flow is ''one'' when it happened and ''zero'' when there was no occurrence. Logistic regression has been used frequently for multivariable analysis (Ohlmacher and Davis 2003;Ayalew and Yamagishi 2005;Lee and Pradhan 2007). In this study logistic regression is used to estimate the correlation between the dependent variable (occurrence of debris flow) and independent variables (influencing attributes) by evaluating probabilities using a logistic function. Results of logistic regression analyses using datasets from 30 debrisflow cases show that the four susceptibility value attributes had weights of 0.27, 0.24, 0.26, and 0.23, respectively.
Since the weights showed no significant difference, the attributes were considered to have identical weights. The vulnerability value indicates whether a debris flow will actually damage or have an impact on expressway sections. The vulnerability value is assessed by two attributes: the volume of margin area available in which to deposit debris-flow materials before reaching expressway structures; and the size of drainage facilities running through the expressway. For the acquisition of the attribute values, expressway design files provided by the Korea Expressway Corporation were used. Each of the attributes is given points ranging from 0 to 5 based on a grading standard (Table 2), and these points are added up to provide the total vulnerability value of 10 points.
By combining the calculated susceptibility value and vulnerability value (Fig. 1), a hazard class is given for a target expressway section. In the table of severity rating shown in the chart of hazard class ( Fig. 1), the x-axis and y-axis indicate the vulnerability value and susceptibility value, respectively. Through investigations on past debrisflow occurrences, the hazard classes were categorized according to the rainfall recurrence period for expressway design purposes. Hazard class S indicates a likelihood of debris-flow occurrences in areas with rainfall recurrence periods of 2-5 years. Hazard classes A, B, C, and D have rainfall recurrence periods of 5-20, 20-50, 50-100 years, and over 100 years, respectively. Hazard class E indicates an area with a very low likelihood of debris-flow damage (Korea Expressway Corporation 2009).

GIS-Based Framework for Debris-Flow Assessment
Based on the processing of attributes included in the KEC method, a real-time debris-flow hazard assessment framework for expressways was newly developed using the ArcGIS 10.2 software (Fig. 1). Various ArcGIS tools, such as the Spatial Analyst tools and the Analysis tools, were used for a quantitative and objective assessment of the attributes using rainfall datasets on a real-time basis.

Attribute Processing for the Susceptibility Value
In order to process watershed slope and valley slope datasets, DEMs provided by the NGII of Korea were used. Numerical maps with the highest resolution were those of a scale of 1:1000. But 1:1000 scale numerical maps were only provided for major urban areas. Because numerical maps of the highest resolution provided for the entire Korean Peninsula were those of 1:5000 scale, numerical maps with a scale of 1:5000 were employed in attribute processing to obtain the susceptibility value.
Of the entities within the DEMs, only the polyline and point entities having elevation value were extracted from the numerical map and used for construction of the DEMs. Because the system focuses on the debris-flow hazard assessment of expressway facilities, the expressway layers were selected. For the processing of slopes in the surrounding area of expressways, the elevation layers were also selected (Fig. 2a). With the elevation layers of the DEMs, the elevation and slope raster images with the smallest cell sizes possible were obtained, as shown in Fig. 2b, c. Because the minimum cell size that could be considered with 1:5000 DEMs was 5 m, the raster with cell sizes of 5 m by 5 m was processed. Based on the elevation raster, the flow direction datasets were computed. The Flow Direction tool creates a raster of flow direction from each cell to its steepest downslope neighbor (Olivera et al. 2002a). From the flow direction raster, the flow accumulation datasets were obtained (Fig. 2d). The Flow Accumulation tool creates a raster of accumulated flow into each cell. With a flow accumulation grid as the accumulated weight of all cells flowing into each downslope cell, valleys can be defined through the use of the flow accumulation value (Olivera et al. 2002b) (Fig. 2e). For a more accurate visualization of valley areas, the properties of the flow accumulation grids were altered in various ways. Through trial and error and by comparing with the actual site survey (or the valley layer of the digital map) and the flow accumulation grids, the accumulation grids (white lines in Fig. 2e) having 0.1 STD of valley shape (geometry of valley based on the DEMs) were applied to define the boundary of valleys.
After setting a pour point (output point) on the route of the assessed expressway, the flow direction and pour point were taken into consideration to obtain the watershed. The Watershed tool creates the drainage areas that contribute to the flow from the land surface to the water system (Esri 2002) (Fig. 2f). Through the Extract by Mask tool, the slopes of the cell in the watershed area were obtained (Fig. 2g). Through the histogram in the raster properties, the values for the attributes of the mean watershed slope and area percentage of the watershed with slopes over 35°w ere acquired. Because the slope raster can indicate the steepest slope with regard to the surrounding pixels and the flow direction, but not the slope in the valley direction, the previously obtained slope raster cannot be used to derive the attributes related to the whole valley. Other means should be used for the acquisition of the valley slope data. An approach that assesses the slope through the elevation raster in the valley direction was proposed and applied.
From the flow accumulation layer, the valley shapes in the watershed were obtained. The valley paths were manually plotted on the elevation layer. The elevations of the cells in the path of the valley were obtained through the Extract by Mask tool of ArcGIS using the plotted valley path and elevation layer. With the length of the plotted valley path, and the total elevation difference between the expressway and the highest point in the valley path, the mean valley slope was simply calculated as follows: DH is the elevation difference between the expressway and the highest point in valley path, while L path is the length of the plotted valley.
In the estimation and prediction stage, the valley in which a debris flow will initiate is unknown. Therefore, for watersheds with more than one valley, all valleys in the watershed were considered in the process. The mean valley slopes for each valley were calculated, and were averaged using the valley lengths as weight factors: showing abrupt valley direction changes at right angles (Fig. 3a). In order to appropriately assess the slopes in the actual flow directions of the valleys for the calculation of the length percentage of the valley with slopes over 15°, additional ArcGIS tools were used. The extracted valley path elevations of cells were converted to points through the Raster to Point tool, and were assigned to the map coordinate system (Fig. 3b). With the Buffer tool, a buffer of 5 m was set around the manually plotted valley path (Fig. 3c). In order to obtain the valley path points that are positioned inside the 5 m buffer zones of the plotted valley path, the Intersect tool was used. Through the process, only the valley path points that were in the vicinity of the actual valley path were obtained (Fig. 3d). Using the shape value (as file type of ArcGIS) of extracted points, the slope between the DEMs in the direction of the valley travel path was calculated using Microsoft Excel. Through the calculation results of the type of shape files, the distances between cells with slopes of over 15°were obtained, which in turn allowed the calculation of the length percentage of a valley with slopes over 15°. According to the mean valley slope calculation process in watersheds with more than one valley, the length percentage of valley with slopes over 15°w as calculated using a similar process: Overall length percentage of valley with slopes over 15 (L 15 ) i is the length of valley with slopes over 15°for the individual valley.

Attribute Processing for the Vulnerability Value
From the expressway design files (Korea Expressway Corporation 2009), the volume of area available for sedimentation was calculated by simplifying the area as a triangular pyramid bounded by the valley and expressway embankment (Fig. 4). The volume of the margin area of sedimentation is calculated using Eq. 4. The drainage facility information of each watershed was obtained through the expressway design files, and the attribute points were determined according to the grading standard set by the KEC (Table 2).

Quantification of Rainfall Criteria for Hazard Class
According to the KEC method, the rainfall recurrence periods (ranging from 2 years to more than 500 years) were used as trigger values of the landslide hazard. But to immediately determine the rainfall event by linking the network sever of meteorological observatories and realtime estimates of the landslide hazard, it is impossible to directly apply the rainfall recurrence periods for the proposed framework. Thus, rainfall recurrence periods were quantified as representative accumulated rainfall criteria: 1 h rainfall (Yune et al. 2010), 6 h rainfall (Oh and Park 2013;Ham and Hwang 2014), and 3 day rainfall (Yoo et al. 2012). The baseline rainfalls were set at the lower bounds of 1, 6 h, and 3 day rainfall data from AWS and at the major interchanges at 18 stations on the nine expressway sections in which the 163 debris-flow events (the highest number of debris-flow events from 2000 to 2009) occurred in Korea. Also, various sets of these rainfall baselines using accumulated rainfall data recorded from AWS were cross-validated with the actual damage state for early warning. These referred to the recurrence intervals of rainfall corresponding to the hazard class in the KEC-based debris-flow assessment (Choi et al. 2015). From the correlations between recurrence period and accumulated rainfall data regarding to the hazard class, the range and lower bound of 1, 6 h, and 3 day accumulated rainfall datasets were determined (Table 3). Thus, the rainfall criteria were set with the lower bound values of accumulated rainfall ranges for five hazard classes based on conservative hazard management. The severity levels of debris flow were defined as three levels: safe, caution, and warning (Fig. 1). If the measured rainfall datasets do not correspond to all three cases of accumulated rainfall datasets, the target expressway section is considered as safe from debris-flow hazard. If more than one indicator exceeds the three baseline rainfalls with regard to the five hazard classes, the debris-flow hazard needs to be forecast as a level of caution or warning. The criteria are appropriate for conservative risk management because of the cross validation between rainfall datasets and hazard class in the case of actual damage status on the representative nine expressway sections (see Sect. 5.2). The quantified rainfall criteria for the hazard class of debris flow create difficulties for real-time management of debrisflow hazard on the national expressway, although it is desirable to apply the probability of the possible rainfall recurrence period (Choi et al. 2015). Therefore, simulations using relationships between recurrence period and accumulated rainfall data at locations where debris flows have occurred is required to establish the appropriate rainfall criteria for debris-flow risk assessment.

Safety Management System for Debris-Flow Hazard
The GIS-based safety management system for debris-flow hazard consists of a database (DB) and four systematic modules (Fig. 5). The database contains all field data and processed data in the system. The submodules execute specific functions that manage and utilize information in the database. These include the input module, the debrisflow hazard assessment module, the decision-making module of debris-flow countermeasures, and the real-time, debris-flow hazard assessment module linked to the rainfall monitoring network. The framework, which includes all these functions, focuses on user-friendliness and real-time applications.

System Database
A database is the backbone of the developed debris-flow framework. It stores not only collected data such as the digital numerical map, expressway, and real-time based rainfall monitoring data (which is standardized using AWS, transmitted from meteorological observatory server), but also contains processed data obtained from the debris-flow hazard assessment and real-time prediction modules (Fig. 6). The DB contains all data as alphanumeric values according to standard data formats, which represent the outcome of data classification and standardization with spatial information. The data stored in the database can be easily utilized in the framework. The digital numerical map can be used as the basic topographical information of the system because it offers an easy way to construct topographical information for a target area. The expressway information is composed of organization categories for managing Korean expressway route data, and spatial datasets designed for various coordinates systems (longitude and latitude, GRS80, Google coordinates). The data format of rainfall events (belonging to the AWS network) was arranged previously for an expressway in Korea. The general information about rainfall observatory and time-series rainfall monitoring information were standardized considering real-time connection with AWS network. Furthermore, the results of hazard analysis and alert information were built into the database at the same time as the framework of debris-flow assessment.

Submodules for Rainfall and Safety Management
Input function provides an effective way to store and arrange all collected data, including electric or nonelectric documents, general information, digital numerical maps, expressway data, rainfall monitoring data, and analyzed data, according to a standardized data format in the DB. Before applying the rainfall criteria for severity levels of debris flows, the target expressway sections are grouped within the spatial range of influence of AWS observatory (5 km radius), around the nearest AWS in Korea (Fig. 7). The influencing areas vary in different localities of a watershed, but we group the nearest expressway sections in accordance with the average distance (5 km). In order to link the site-specific rainfall datasets with expressway sections, additional investigation about influencing areas of AWS observatories considering the hydrological characteristics of mountain areas is important, but this study did not consider this complexity. The reclassified AWS datasets for rainfall criteria (1, 6 h, and 3 day accumulated rainfall) are linked to the grouped expressway sections using the system database.
Using the real-time-based rainfall transmission module, as the rainfall accumulation records are transmitted from the AWS in real time, the rainfall event protocol is retained to consider the event trigger time and transferred event file. In Fig. 7, as soon as a rainfall event is detected, the rainfall event time is transmitted from the AWS to the developed system. Then, the generated event file (stored in a oracle server) is transferred in company with the completion message, which is the computational message to inform the transmission or reception of event file and enable the DB, to operate the safety management system. At the developed system server, an event file is calibrated based on the rainfall criteria, according to the file transfer protocol. The rainfall accumulation datasets, including completion message of event time and file, are stored in system database.
The alerted rainfall datasets are compared with rainfall criteria for hazard class (Table 3), and the severity level can be determined in 30 s after the occurrence of the rainfall event.
The real-time debris-flow assessment framework has four functional phases within the database based on the proposed schematic sequence of the debris-flow hazard assessment module (Fig. 7). In the first phase, linked with the digital numerical map and the DEM, the watershed DEM and valley layer are extracted using the ArcGIS desktop program and input system DB. In the second phase, the susceptibility value and vulnerability value for the target route are constructed into DB and combined with geospatial information. In the third phase, reliable rainfall monitoring data for the target route is transmitted from the widely distributed AWS server on a real-time basis, the routes of the completed site investigation for debris-flow hazards are grouped into the same datasets and focused on the adjacent rainfall station in target areas. The rainfall values for debris-flow hazard assessment are automatically computed based on the monitoring criteria with reference to recurrence periods of rainfall (Korea Expressway Corporation 2009 when the rainfall monitoring data is input to the DB. Once the rainfall threshold level for debris-flow hazard is reached (from the KEC method), the severity levels (safe, caution, and warning) are determined using map symbols (blue, yellow, and red symbols, respectively, in Fig. 7) at the target route in real time. In addition, sound signals and message windows are alarmed to notify the expressway administrator of the hazard status.
A decision-making module of debris-flow countermeasures was developed by standardization of site investigation reports as part of a nation-wide expressway survey in mountainous areas of Korea. As soon as the severity levels based on the rainfall monitoring data in real time are decided, it is possible that optimum countermeasures for debris-flow minimization can be rapidly put in place for a threatened target route depending on whether a caution or warning severity level is activated. Facilities (ring net, debris barrier, and so on) for debris-flow prevention can be preposition beforehand based on the site investigation for target routes having potential of debris flow. The method and cost of construction can be studied and visualized in a report, corresponding to the anticipated severity level in near real time.
The debris-flow hazard prediction function shows all the attributive information in the database by using tables and graphics according to its characteristics, on screen or as a document. Also, all data in the DB can be exported as a chart or graph. The graphic functions simultaneously display interpolated data and field data over a designated area. All of the charts and graphs can then be printed. Especially, the debris-flow hazard can be visualized and forecast as 2D maps overlaid with satellite images. Also, the severity level of the target route can be determined using zonation criteria in real time.
In this proposed framework, the computer-based method for the real-time assessment of spatial debris-flow hazard was embedded based on a stand-alone system developed using Microsoft Visual BASIC and the Esri ArcGIS developer tool (Lee and Wong 2001;Esri 2006). The ArcGIS developer tool was mainly used for the development of the database, evaluation of the results, and spatial visualization. Several assumptions and precedent assessments are needed to estimate the possible debris-flow hazard for a target site in real time, at the point at which the debris flow occurs. Especially, the precedent procedures consist of the building of the database, DEM construction, and determination of susceptibility value and vulnerability value in order to consider the overall debris flow potential in the target area, prior to the occurrence of rainfall and debris flow. As the debris flow occurs near the target site, the possible severity level can be estimated in real-time by using the rainfall data monitored from the AWS server.

Application of the Assessment Method on Selected Sites
The developed real-time GIS-based debris-flow hazard assessment framework is applied to three debris flowdamaged expressway areas in Korea to verify its validity. The standard set of the original KEC method was subsequently revised based on the case study.

Condition of Application
The Pyeongchang area of the Yeongdong Expressway, the Deogyu Mountain area of the Daejeon-Jinju Expressway, and the Juksan-Geochang area of the 88 Expressway in Korea were selected for investigation of their debris-flow hazards (Fig. 8). Many debris flows occurred in the Pyeongchang area of the Yeongdong Expressway during a heavy rainfall event (244.0 mm/day, 66.0 mm/h) in the summer of 2006 (Table 4). Debris flows occurred in the Deogyu Mountain area of the Daejeon-Jinju Expressway in the summer of 2005. The rainfall intensity in the region at that time was 312.0 mm/day and 54.5 mm/h. The Juksan-Geochang area faced several debris flows during a heavy rainfall event in the summer of 2006. The rainfall intensity at that time was 121.0 mm/day, and 31.5 mm/h. When comparing all target areas, the Pyeongchang and Deogyu Mountain areas had similar rainfall intensities, indicating the same rainfall recurrence periods, whereas the Juksan-Geochang area had the lowest daily rainfall. For each area, a test bed area was set for the application of the assessment method. A test bed with a length of 11 km along the Yeongdong Expressway was chosen for the Pyeongchang area. The lengths of the test beds for the Daejeon-Jinju Expressway and 88 Expressway were 3 and 2 km respectively, presented as black dotted lines in Fig. 8.

Verification of the Framework for Debris-Flow Assessment
All existing watersheds in the selected expressway test beds were analyzed. Of all the watersheds in the selected regions, the areas with expressways positioned on bridges and tunnels, or near vast areas of fields were excluded from the analysis due to their very low likelihood of damage.
Since reported debris flows were based on the damage made to road structures, those areas without any damage were not reported. No debris-flow damage was reported in areas where the space for possible sedimentation deposit is vast, and consequently these areas were not considered in the analysis process. After the exclusion of low-risk sites, watersheds were then classified according to whether or not debris-flow damage had been reported. Areas with debrisflow damages reports are referred to as ''occurrences,'' and those without damage reports are referred to as ''nonoccurrences.'' As a result, 18 debris-flow occurrences and 14 nonoccurrences were analyzed for the Pyeongchang area (divided into 32 subsections) based on the proposed framework for debris-flow assessment; 12 debris-flow occurrences and 8 nonoccurrences were analyzed for the Deogyu Mountain area (divided into 20 subsections); and 8 debris-flow occurrences (evaluated as 1 warning and 7 caution level according to the developed system) and 8 nonoccurrences were analyzed for the Juksan-Geochang area (divided into 16 subsections) as shown in Fig. 9. The developed safety management system for debris-flow hazard constructed the debris-flow hazard map according to hazard class for watershed and forecasted severity level linked with AWS (hourly max rainfall datasets) on realtime basis at the Jusan-Geochang area (Fig. 9).
Applications of the KEC method show results that roughly coincide with the actual debris-flow occurrences and nonoccurrences are presented in the table of hazard value (varying from 0 to 30) and hazard class (Fig. 10). Occurrence cases are roughly positioned in the upper righthand side, which indicate higher susceptibility and vulnerability values, whereas nonoccurrence cases are located on the lower left side, with relatively lower susceptibility and vulnerability values. Although this tendency may seem correct to some extent, it does not always show flawless results. In the C and D hazard classes, both debris-flow occurrence and nonoccurrence cases are mixed up, not always indicating a result in which occurrences have higher hazard classes, and nonoccurrences fall into lower classes. Debris flows occurred even in areas with an E hazard class, which was initially expected to indicate a very low likelihood of debris flow. The total sum of susceptibility value and vulnerability value even showed a greater number of nonoccurrences than the occurrences due to the low hazard values of occurrence cases.
Although the sites showed results that generally presented higher susceptibility values and vulnerability values  in occurrence cases, they did not perfectly represent the total hazard value differences. But the effectiveness of the proposed framework was validated by the three cases because the potential of debris-flow hazard is generally high for the actual occurrence cases. The average of the hazard value of occurrence cases is 1.13 times that of nonoccurrence cases (Table 5). In order to appropriately represent the differences between occurrences and nonoccurrences of debris flows, modifications can be made on the grading standard of the KEC method using more validation tests based on the proposed framework, because the KEC standard set was established using logistic regression with a small number of samples of debris-flow damaged expressway sections. Specifically, grading standards for attributes need to be revised by taking into consideration the attribute values of both occurrences and nonoccurrences, maximizing the distinction between the two instances.

Modifications to the Framework
The existing KEC method assigned points to each attribute according to a grading standard set based on a few case studies of debris-flow occurrences. But it did not take the nonoccurrence cases into consideration (Kim 2012). In addition, the KEC grading standard set used all data samples collected for the study. In order to appropriately represent the differences between debris-flow occurrences and nonoccurrences, modifications can be made to the grading standard using additional tests. The three applications in this study included a large number of debris-flow occurrence sites and can be used for this purpose. To modify the grading standards, the average values of attributes in the three target areas were calculated. When calculating these average values, the mismatched cases (between estimated severity level and actual damage) were eliminated from the calculation. Likewise, occurrence cases that showed considerably high hazard values, such as 25 and higher, and nonoccurrence cases with total values of 11 and lower, were considered as outliers for statistical normalization and were not considered when calculating the average values. After calculating the average values of each attribute for the occurrence and nonoccurrence cases of the three target areas, the highest attribute values of the occurrence and nonoccurrence cases were respectively defined as the upper and lower boundaries of each scoring  index. The range of each scoring index of the four susceptibility value attributes was divided into six classes using equal intervals. Logistic regression was performed to reestablish the grades of scoring indices in the range of 0-5. Results of logistic regression analysis for the four susceptibility value attributes indicate that the six grades have weights of 0.29, 0.24, 0.19, 0.14, 0.10, and 0.05, respectively. Accordingly, the redefined grades of the scoring indices were assigned points from 0 to 5. The modified susceptibility grading standard is presented in Table 6. The vulnerability grading standard was also modified based on the debris-flow occurrence and nonoccurrence cases of the three target areas (Table 6).
Although the grading standard of the original KEC method was established based on 30 field applications (KEC 2009), the scale of damage by debris flow during the considered period was very small (Kim 2014). Considering the fact that only a very small number of the considered sites had deposit areas with volumes exceeding 2000 m 3 , the grading standard was modified to reflect this range, that is, the highest class of debris storage was altered from above 5000 m 3 to above 2000 m 3 . The other standards were also modified accordingly.
Based on the same condition of application for the three target areas, the hazard value and hazard class were determined using the revised grading standard. On the scoring chart, the occurrence cases are distributed in hazard classes S, A, and B, which have a high potential of debris flow, as shown in Fig. 11. In contrast, the points of nonoccurrence cases are distributed in hazard class C, D, and E, which have comparatively lower potential for debris flow. In addition, the average of the hazard value of occurrence cases is 1.58 times that of nonoccurrence cases (Table 7); this is significantly higher than the result (1.13) using the original KEC standard. Moreover, the method has potential to be improved with more data accumulation and more case histories. Also attributes other than those regarding slope should be considered, such as watershed size and valley curvature. According to Kim (2014), with larger watershed sizes both the debris-flow initiation risk and occurrence risk increase. With larger valley bending ratios, more debrisflow materials are subjected to sedimentation, lowering the possibility of damage on road structures.

Conclusions
In this study, to provide preliminary information for rapid decision making of evacuations or restoration measures, and to prevent secondary disasters caused by debris flows, a GIS-based real-time framework of debris-flow hazard assessment for expressway sections was proposed based on the KEC debris-flow hazard assessment method. First, the KEC-based method was standardized using Esri ArcGIS, which enabled various attribute datasets to be acquired in an objective and quantitative manner with a fixed data acquisition sequence. Second, considering the real-time link with the AWS network, the rainfall recurrence periods were quantified as accumulated rainfall using trigger values of actual debris-flow events. Based on the optimized KEC debris-flow hazard assessment method, a safety management system for debris-flow hazards was developed based on the GIS platform. For field applicability evaluation, the framework was applied on three expressway sections in Korea: the Pyeongchang area of the Yeongdong Expressway, the Deogyu Mountain area of the Daejeon-Jinju Expressway, and the Juksan-Geochang area of the 88 Expressway. The reliability of the assessment method was investigated by comparing with actual debris-flow occurrence and nonoccurrence cases. A revised grading standard set was proposed and validated in a case study of damaged expressways in comparison with the KEC method. In the future, additional attributes including watershed size and bending of valley should be considered in the hazard assessment.