Application of a Novel Hybrid Method for Flood Susceptibility Mapping with Satellite Images: A Case Study of Seoul, Korea

: This paper proposes a novel hybrid method for ﬂood susceptibility mapping using a geographic information system (ArcGIS) and satellite images based on the analytical hierarchy process (AHP). Here, the following nine multisource environmental controlling factors inﬂuencing ﬂood susceptibility were considered for relative weight estimation in AHP: elevation, land use, slope, topographic wetness index, curvature, river distance, ﬂow accumulation, drainage density, and rainfall. The weight for each factor was determined from AHP and analyzed to investigate critical regions that are more vulnerable to ﬂoods using the overlay weighted sum technique to integrate the nine layers. As a case study, the ArcGIS-based framework was applied in Seoul to obtain a ﬂood susceptibility map, which was categorized into six regions (very high risk, high risk, medium risk, low risk, very low risk, and out of risk). Finally, the ﬂood map was veriﬁed using real ﬂood maps from the previous ﬁve years to test the model’s effectiveness. The ﬂood map indicated that 40% of the area shows high ﬂood risk and thus requires urgent attention, which was conﬁrmed by the validation results. Planners and regulatory bodies can use ﬂood maps to control and mitigate ﬂood incidents along rivers. Even though the methodology used in this study is simple, it has a high level of accuracy and can be applied for ﬂood mapping in most regions where the required datasets are available. This is the ﬁrst study to apply high-resolution basic maps (12.5 m) to extract the nine controlling factors using only satellite images and ArcGIS to produce a suitable ﬂood map in Seoul for better management in the near future.


Introduction
Globally, flooding and storms as natural disasters have become a serious problem currently, accounting for 67% of all the natural disasters worldwide [1]. Thus, flood susceptibility mapping has recently received significant attention [2][3][4]. In addition, flood damage, which is a significant factor in extreme flood events, is expected to increase in the near future because of climate change [5]. Meanwhile, population increase is expected to severely affect the increase in floods, as it leads to a loss of vegetation owing to construction and increased exposure to flood events. Thus, comprehensive flood management is essential to reduce the effects of floods on human lives. To achieve this aim, a flood susceptibility map can be an effective method, offering information on the probability of flood occurrence, magnitude of the disaster, and location of inundation for flood control.
Flood occurrence is an unpredictable event; therefore, estimating the time of occurrence of a flood and its location is one of the biggest challenges being faced recently. Traditionally, floods are measured using stream gauges installed in water bodies located near populated areas. These gauges constantly monitor the water levels and collect data. Although this is a direct method for determining flood locations, it is costly and timeconsuming. Thus, the integration of remote sensing (RS) and geographic information system (ArcGIS) technology has received substantial attention in the field of flood susceptibility mapping. RS can cover different kinds of images at large spatial scales that are critical Peninsula between the latitudes of 37 • 42 to 37 • 70 N and longitudes of 126 • 76 to 127 • 18 E. It can experience high volume of rain each year which can cause flooding, even though some parts of this city are surrounded by densely forested mountains.
Previous studies have used a combination of directly extracted man-made region maps and ArcGIS maps to develop final flood maps; however, in this study, only high-resolution (12.5 m) satellite images were used to develop flood susceptibility maps. Thus, this study can be considered novel as it uses only maps extracted from ArcGIS and satellite images without direct intervention. Overall, the main objectives of this study were (1) to identify flood potential zones using ArcGIS and AHP and (2) to validate the flood susceptibility map by using a real flood map.

Study Area
According to government statistics, Seoul is located in the northwestern part of the country with a total area of 605.25 km 2 [27] and is composed of 25 administrative districts, called GU [28]. The population of Seoul Metropolitan City, which is the world's tenth largest capital (with more than a quarter of the entire population of South Korea), was approximately 10 million in 2020 [28]. The Han River (Hangang River) is a major river in South Korea and the fourth longest on the Korean Peninsula, after the Amnok (Yalu), Tuman (Tumen), and Nakdong rivers, and crosses Seoul from east to west [29]. The river starts as two smaller rivers in the eastern mountains of the Korean Peninsula, merging near Seoul.
The Han River, located at 126 • 59 40 E and 37 • 33 59 N, divides the city area into northern and southern portions. It flows through Seoul and then merges with the Rimjin River shortly before flowing into the Yellow Sea. The total length of the Han River is approximately 494 km [30]. Seoul was selected as the case study because more than 80% of rainfall occurs during the monsoon months (July-August) with high potential risks of flood events (Figure 1). Particularly in 2010 and 2011, extreme rainfall wreaked havoc on the city's southern area, and caused extensive damage to property affected by the flood with worth approximately $65 million [31]. A torrential downpour of 259.5 mm per day on 21 September 2010, caused significant damage to the city center, especially around Gwangwhamun and Cheonggyecheon. Furthermore, on the morning of 27 July 2011, heavy rainfall in the central part of Korea resulted in a rainfall of 164 mm, exceeding the frequency value of over 100 years [28].
The study area's climate is characterized by dry winter (called "Dwa" in the Köppen climate classification) with four distinct seasons; however, temperature is extremely different between the hottest part of summer and the coldest part of winter. Winters are cold and dry, with the average low temperatures ranging from −6.1 • C (21 • F) to −2.8 • C (27 • F), and temperatures often drop to −17.8 • C (−0 • F). The spring and autumn seasons are warm and dry, with calm weather. During the year, the annual sunshine is 2066 h and the average humidity is 65%. The annual precipitation rate in this area is 1370 mm, with summer accounting for more than 70% of the city's rainfall.
Due to the high population density and heavy rainfall, the flood occurrence of a large city such as Seoul Metropolitan City becomes a more serious issue; therefore, efficient management to prevent flood events is essential to reduce the considerable cost of disaster recovery. The determination of the main factors influencing flood vulnerability was the first step in preparing the dataset for this study.

Satellite Images
In this study, the ALOS PALSAR (Phased Array type L-band Synthetic Aperture Radar) sensor was used to derive digital elevation model (DEM) with 12.5 resolution. Regarding this issue, the DEM was downloaded from the website https://search.asf.alaska.edu (accessed date: 20 June 2020). On the other hand, Landsat 8 images (i.e., thematic mapper), with the lowest percentage of cloud cover, were downloaded to generate the land use classification map in ArcGIS using data management tools. By using the DEM as a basic layer, other basic thematic maps (slope, drainage density, topographic wetness index (TWI), flow accumulation, distance from river, and plan curvature) were extracted using ArcGIS. Finally, for interpolation of Landsat 8 images, Erdas Imagine and ArcGIS were used to prepare a land use category map of Seoul, considering google earth images to identify different categories of land use. All above mentioned controlling factors were converted to a raster grid with 12.5 m × 12.5 m grid cells for application of the AHP model. The flood susceptibility map was prepared after overlaying thematic layers and validating the consistency of the AHP model. A flowchart of the methodology is presented in Figure 2.

Delineation of Thematic Layers
Determining the relevant flood controlling factors is a necessary step in preparing the flood susceptibility map. In this study, nine thematic layers were selected for flood mapping such as DEM, land use, slope, topographic wetness index, plan curvature, river distance, flow accumulation, drainage density, and mean rainfall ( Figure 3). As can be seen in this Figure, all thematic layers were extracted from DEM using spatial analyst tools in ArcGIS except a land use map which was prepared from Landsat 8 images with integration of ArcGIS. Taking advantage of the fact that most of the area of Seoul has become an urban area indicated from the new land use map prepared in this study for the year of 2020, the soil map of the area was ignored because of relatively low impact on potential floods. After these layers were prepared, they were converted to a raster format of 12.5-m resolution. Then, the AHP model was used for reclassification using ArcGIS.

DEM
Elevation plays a dominant role in flood susceptibility studies and is used as a basic map for slope, flow accumulation, TWI, plan curvature, and drainage density. As mentioned earlier, a DEM with 12.5-m resolution was obtained from the ALOS PALSAR sensor and was processed in ArcGIS. According to this type of map, floods are less likely to occur in areas far away from the water bodies because of their high elevation. In Figure 3a, high-altitude areas are marked in green and low-altitude areas are marked in dark pink. As areas along the Han River at the central part of Seoul have lower elevation, they are more susceptible to flooding.

Slope
The topographic slope is a common constraint on water velocity and flood strength and is always presented as a percentage or degree. Thus, areas with flat topography and a very low regional slope are more susceptible to flooding when a large volume of water becomes stagnant. As shown in Figure 3b, the southwestern part of Seoul has a slope ranging from 1 • to 6 • , which is more prone to floods, whereas the northern and southern parts of the city with high slopes have more protection against floods.

Distance from River
The distance from river is another main controlling factor in flood susceptibility mapping which was prepared using the DEM and streamflow. It is an extracted layer from flow accumulation layer using spatial analyst tools in ArcGIS map was prepared using a buffer tool in ArcGIS and classified into six classes in raster format (Figure 3c). According to this type of map, floods are less likely to occur in areas far away from the water.

Mean Rainfall
Among the eight mentioned factors, rainfall seems to have a direct effect on floods. Thus, it is considered to be one of the most crucial factors for flood susceptibility mapping. To prepare the mean rainfall map using ArcGIS, the daily rainfall data for the past 33 years (1987-2019) were downloaded from https://data.kma.go.kr (accessed date: 20 June 2020) for two stations, named Daegokgyo (Latitude 37.47 • and Longitude 127.12 • ) and Seoul (Latitude 37.57 • and Longitude 126.97 • ). This map was extracted using the average monthly rainfall data in ArcGIS. In this study, Inverse Distance Weighting (IDW) was used to extract average monthly rainfall of 33 years. In Figure 3d, dark and light pink areas are more prone to flooding and are located in the center and west part of Seoul.

Flow Accumulation
According to Zhou [32], one of the most significant factors that has been commonly measured and used in various studies is flow accumulation, which represents the total amount of water that has flown from the upstream areas to a particular point within the catchment. This map was extracted using spatial analyst tools in ArcGIS. A higher flow accumulation value indicates a greater likelihood of water accumulation and thus a greater likelihood of flooding. In Figure 3e, the central part of Seoul (near the Han River) shows a high amount of flow accumulation.

TWI
This is a secondary topographical attribute that depicts the spatial distribution of wetness conditions and can be described as [33]: where α represents the upslope area per unit contour length through a point and tan β is the slope angle at that point. In this study, TWI was extracted using map algebra tool in ArcGIS and divided into four classes by using three basic maps: DEM, TWI, and flow accumulation ( Figure 3f). TWI has a direct relation with flood events, implying that areas with a higher TWI are more vulnerable to flooding.

Plan Curvature
The plan curvature is the curvature in a horizontal plane and is also known as the curvature of a hypothetical contour line that passes through a particular cell. Cells with concave contours have a positive plan curvature, whereas those with convex contours have a negative plan curvature. The plan curvature can be analyzed to extract useful geomorphological details [34]. The concave, flat, and convex planes are the three types of plan curvatures used in this study (Figure 3g). A negative value indicates a concave shape, a positive value indicates a convex shape, and a value close to zero indicates a flat field. This map was created using 3D analyst tools in ArcGIS. The negative value has the highest weight in the heightening curvature map because flood events are more common.

Drainage Density
This factor is calculated by dividing the total length of all streams in a drainage basin by the drainage basin's total area [35] using spatial analyst tools in ArcGIS. The stream layer extracted from the DEM is the basic map to prepare a drainage map. According to Das and Pardeshi [36], a high drainage density indicates considerable surface runoff and low-groundwater-potential zone. In this study, the drainage density varied from 0.77 to 3.14 km/km 2 , as shown in Figure 3h, and was classified into the following four classes: <0.78, 0.78-1.57, 1.57-2.36, and >3.14 km/km 2 .

Land Use
According to Yalcin [37], this factor is crucial in determining several hydrological processes, such as the rate of soil infiltration, evapotranspiration, and surface runoff. As mentioned before, Landsat 8 images with the lowest percentage of cloud cover were obtained from the United States Geological Survey official website, and the ArcGIS and Erdas Imagine were used to process satellite imageries for land use change assessment of this study area. In order to derive supervised land use classification, maximum likelihood classification algorithm was used in this study. The pixel color tone for each land use class was identified using ArcGIS and verified using google earth. For the final stage, Erdas Imagine was used for accuracy assessment of supervised land use classification. As Figure 3i shows, the major land uses in the study area are: (i) urban areas, (ii) several types of vegetation, (iii) water bodies, (iv) forests, and (v) bare land with a kappa accuracy of 0.92.

AHP
After preparing all nine layers, including DEM, mean rainfall, land use, slope, distance from river, flow accumulation, TWI, drainage density, and plan curvature (Figure 3), the AHP model was used to assign different weights to each thematic layer. AHP is the MCDM method implemented by ArcGIS, which was developed by Saaty [38], and this approach involves determining the weights and ranks of different flooding parameters. The total weight of each map of the final integrated layer is computed as where FH is the flood hazard index, w represents the weight of each criterion, and r is the rating of each criterion with the following names: DEM (E), mean rainfall (R), land use (L), slope (S), distance from river (D), flow accumulation (F), TWI (T), drainage density (DR), and plan curvature (C). The FH model was used to assess the vulnerability of flood incidence, and finally, the flood map was categorized into six regions using the weights of each sub-criterion in ArcGIS: very high risk, high risk, medium risk, low risk, very low risk, and out of risk. In this study, pairwise comparison matrices were developed using the AHP extension in the ArcGIS, which involved all nine criteria to determine which of them has higher priority. In this regard, Saaty [38] suggested a scale of 1 to 9 for each criterion compared to the other elements (Table 1). After weighting all elements, the values in the pairwise comparison are normalized and the average random consistency indices ( Table 2) according to the number of criteria (here, the number of criteria is nine) are chosen to express the consistency ratio with the following equation: where CR is the consistency ratio, CI represents the consistency index, RI is the random consistency index, λmax represents the principal eigenvalue of the matrix, and n is the number of components or criteria in the matrix.

Nine Controlling Factors with Their Relative Weights
In this study, nine thematic layers were considered for flood susceptibility mapping and their weights were determined from the AHP model using pairwise comparison matrix. It should be noted that, the novelty of this study lies in the preparation of maps from satellite images without direct intervention. Thus, all nine maps were extracted, using collected DEM from satellite images, to prepare the flood susceptibility map. The details of the pairwise comparison matrix is listed in Table 3. According to Saaty [38], the CR is acceptable if the value is <0.10, which indicates a reasonable level for the pairwise comparison matrix. In the ArcGIS program, when the sentence "The matrix is considered consistent enough" is displayed, it means that the result is accepted. Subsequently, all parameters are assessed individually and are classified into different subcategories. Table 4 shows the different ranks for each sub-category using the AHP extension in the ArcGIS program.  According to experts' opinion, elevation plays the most effective role in flood mapping [40,41], and thus receives the highest weight in the method developed by Saaty [38], as shown in Table 3. This means that regions located at higher elevations have lower potential to be flooded and vice versa for lower-altitude areas. Thus, the highest altitudes receive the lowest weight, as shown in Table 4. According to the study conducted by Das and Pardeshi [36] on hydrology, slope also plays a crucial role in directing water to downstream areas and is an important factor after elevation. Regarding this issue, as the slope increases, the flow rate also increases, and the infiltration process is reduced. Thus, high amounts of water flowing in streams often leads to flooding. Distance from the river, as a geomorphic factor, is another important factor in flood occurrence. Regions far from the river have lower vulnerability to flood occurrence. The drainage network and flow accumulation had almost the same weight. Flow accumulation refers to the concentration of flow and is directly related to flood occurrence [42]. Meanwhile, a higher density of drainage network simplifies the draining runoff, and thus reduces the risk of flooding.
Land use of area is the type of natural and man-made land. In this study, the following five categories of land use were adopted: water bodies, urban areas, bare land, forest, and vegetation. Some activities such as deforestation and urbanization have serious effects on flood occurrence. In highly built-up regions, flood situations occur frequently due to an increase in population. As shown in Figure 3i, the urban area in red has the largest area and exhibits the highest effect on flood occurrence. TWI represents the effect of topography on runoff generation and flow accumulation in river catchments [41]. It has a direct effect on flood occurrences, which means that regions with higher TWI have a higher vulnerability to flooding. Among the nine factors, plan curvature has a minor effect on floods, which cannot be neglected. In this study, the following three shapes were considered: concave, flat, and convex. Negative values represent concave shapes, positive values show convex shapes, and near-zero indicates flat areas. In the heightening curvature map, the negative value has the highest weight because it plays a more important role in flood occurrence. Rainfall, as the final layer, plays a significant role and has a direct relation with the Han River discharge; thus, it can control the flood occurrence.
After determining weights for each of the criteria and sub-criteria, the AHP extension in ArcGIS was used to find the optimum CR for the highest accuracy. When the statement "The matrix is regarded to be consistent enough" appears in the produced table, it implies that the result has been accepted and this weight is suitable for use in reclassification and overlay operations. In reclassification, all criteria were divided into five sub-criteria by considering the weights developed in AHP.

Reclassification for Each Thematic Layer
In this study, nine criteria were converted into the raster format to identify flood sensitive areas. AHP was applied for the weighting, ranking, and reclassifying of these maps in ArcGIS. These criteria can provide a flood hazard location in a specific region when combined with the AHP technique. According to Saaty [38], a pairwise comparison matrix and the amount of CI and CR were established using ArcGIS. The results showed that CR is acceptable because it has a value of less than 0.10. This demonstrates that pairwise comparison has a sensible consistency dimension. After obtaining this reasonable result, all parameters were assessed individually and classified into five subcategories using ArcGIS; for instance, high-elevation areas have a lower risk of flooding, whereas low-altitude areas have a higher risk. Thus, higher elevations have a high value, whereas lower elevations have a low value (Figure 4a). The classification was repeated for all thematic layers, based on their effects on floods and their detailed information was summarized in Figure 4. To be more descriptive, ratings were provided to classes in each thematic layer based on its relative importance for flood susceptibility, and all thematic layers were reclassified to five different categories regarding to the Table 4. According to all nine reclassified thematic layers (Figure 4), central part to the western regions are more vulnerable to flood than the northern and southern regions because of low elevation (purple color with rank 1 in Figure 4a), low slope (dark brown color with rank 1 in Figure 4b), and high amount of rainfall (light blue color with rank 1 in Figure 4d).

Flood Susceptibility Map
In Table 4, the weight of the relevant factors and their relationships with flood occurrence are shown for each thematic layer. Based on the values in this table, among the different land use types, water bodies, bare land, and urban area had the highest weights, which indicated a positive influence on flooding, whereas forest area had a positive influence on decreasing floods, and thus its corresponding weight was lesser than that of the other categories. In the case of the rainfall map, the area with a high amount of rain has a higher weight, and this value decreases for areas with low levels of rainfall. This classification is the same for slope, TWI, and flow accumulation. However, in the case of DEM, distance from the river, drainage network, and plan curvature, the weighting is the opposite; thus, areas with higher values receive lower weights. Based on this relation, after reclassification for each criterion, the flood susceptibility map was extracted and categorized into the following six regions using the natural breaking method in ArcGIS ( Figure 5): very high risk, high risk, medium risk, low risk, very low risk, and out of risk. For example, the area with the level of very high risk has the highest elevation, slope, rainfall, TWI and flow accumulation, as well as the lowest distance from river and lowest drainage network. Moreover, the very high risk areas have a negative value for plan curvature and located in urban area and bare lands. Figure 5 shows that the high-floodpotential location is to the western area of Seoul, near the Han River. To better understand this type of flood susceptibility map, it is reasonable to refer to the previous basic maps (Figure 3) because these maps are directly relevant to the final map. DEM and slope play crucial roles in flood occurrence. As shown in Figure 3, the green areas have the lowest slope and are situated in the central portion of the basin, which is approximately 28%. In addition, the central part of Seoul from east to west shows the lowest amount of elevation, which is indicated in dark pink in Figure 3a. That is, elevation and slope can contribute to make the Han River overflow, along with monsoon rainfall in the region. In the study area, the roles of forests and vegetation would be ignored in the central part of Seoul because this location is replaced with urban areas and there are no natural places, except for some parks and man-made areas. Zaharia [43] indicated that places with slopes exceeding 15 • do not favor water accumulation, whereas those with a large cover of forest and dense vegetation and a flat surface favors the retention of excess surface water during a flood.
The results revealed that the area was categorized into the following six flood potential zones: very high risk (4%), high risk (36%), medium risk (37%), low risk (16%), very low risk (6%), and out of risk (1%). The area with a very high risk level is concentrated in the central-western part of Seoul, which is characterized by a low elevation and low slope, with high levels of rainfall. Moreover, the small part of the northern and southern parts of Seoul either exhibit very low risk or are out of risk, as the elevation is very high, and the slope exceeds 40 • . This area, with a dense cover of forest and vegetation, is at risk of flood occurrence.
The flood susceptibility map was validated from comparison results with a real flood map ( Figure A1) obtained from Seoul Metropolitan government for the years of 2001, 2010, 2011, 2012, and 2013. At this point, the total flooded areas ( Figure A1) were merged into Figure 5 and the number of pixels associated with each polygon was counted to evaluate the accuracy and validation of the model. Table 5 and Figure 6 describe the number of cells assigned to each of the six categories and the number of corresponding pixels according to the flood susceptibility map with resolution of 12.5 m in Figure 5, which was extracted using spatial analyst toolbar in ArcGIS for validation purposes. According to Table 5, in total, 94.64% of areas are located in the very high risk, high risk and medium risk regions, and 5.37% is located in low risk, very low risk, and out of risk regions. It means that the proposed method can be useful in identifying flood-prone areas with a priority for mitigation and adaptation approaches for water resources management.  On the other hand, the comparison indicates that the flood susceptibility map obtained from this study is consistent with a real flood hazard and that the western part of Seoul has more potential for flooding ( Figure 5). The accuracy of the AHP technique depends on the assigned weights, which vary for each topic. For instance, elevation [44] and slope [45] as the most important factors in the topic of flood prone areas receive the highest weights. These results indicate that flooding occurs most often along riverbanks and only occasionally in areas far from rivers.

Further Discussion
In flood susceptibility analysis, the major aim is to find areas that may be affected by floods in the future. Thus, it is irrelevant which integration method is used (e.g., AHP as a simple method or AI as a complex method), but the validation of this flood susceptibility map is the most important in terms of determining the unknown future flood events [46]. As the first step, creating high-resolution basic maps helps increase the accuracy of the final map; thus, choosing satellite images of the case study area of Seoul with the lowest cloud cover would be effective in determining flood-prone areas. The use of a real flood map is the second step for validation, which helps decision-makers and administrative bodies improve the accuracy of flood-mapping in management and planning by comparing the mapping results with those obtained using real maps. In Figure 5, the locations that are very high risk and high risk are colored in red and pink, respectively; however, even though moderately sensitive regions are categorized in the third classification, they should not be ignored in terms of high priority for flood management. On the other hand, regarding to the Figure 3i, the vegetation cover has decreased compared to other previous years (not included in this study), which could be one of the factors associated with increasing floods. Thus, the pink and red areas need serious attention from administrative bodies to prevent flood situations in the near future.
The major challenge of this study is the opinion of experts and Saaty's [38] method for preparing the appropriate weights of each parameter in extracting the flood map, because factors such as expert tiredness, lack of sufficient expertise, and lack of appropriate measuring devices greatly influence the weighting of the controlling variables. However, a soil texture map [47] and a geology map [48] were not considered for flood susceptibility mapping because of inaccurate satellite images obtained for the year of 2020, but the new land use map was extracted for the year of 2020 to specify the percentage of increasing the urban areas. At the final stage, the results for the year of 2020 were compared with a map that belongs to the years of 2001, 2010, 2011, 2012, and 2013 ( Figure A1). As it has been shown, little difference was found between the two maps because this study extracted new maps to divide areas in terms of flood susceptibility.
Flood susceptibility mapping through the AHP model is expected to be used in management studies related to risk assessment and mitigation strategies for flooding, as it determines the most vulnerable areas based on physical characteristics that influence flooding risk. Overall, future research will focus on issues related to data consolidation, such as resolution, generalization, and integration. Furthermore, the AHP methodology will be compared to other approaches for evaluating flood susceptibility, in order to determine the best alternative, including more robust validation with data collection.

Conclusions
Floods are the most destructive natural disasters worldwide. Therefore, to achieve sustainable growth, flood susceptibility mapping is required for integrated watershed management. In this study, flood prone areas were identified by integrating the AHP method and geospatial techniques, and for validating the results, a flood inventory map containing flood regions was obtained from the Seoul Metropolitan government [49]. First, nine thematic layers, due to the most effective role in the occurrence of floods, were derived using satellite images and ArcGIS. Then, a flood susceptibility map was extracted using appropriate weights and CR. Approximately 40% of the total area was found to belong to the very high risk and high risk zones, which are located in the western and central parts of Seoul, near the Han River, with a flat slope and low elevation. It is noted that the methodology presented in this study can be a useful and easier approach for creating a flood susceptibility map.
Satellite images were used to detect flood hazard-zones. In particular, Landsat 8 was used to extract a new map of land use/landscape in 2020. Previous studies only used the old version of the land use map and considered the basic maps that were extracted in different years. Therefore, using all maps from a particular year (i.e., 2020 in this study) can be effective to obtain suitable results; this was the focus of this study. For instance, in the case study of Seoul, the growth of urban areas is a significant factor in the preparation of land use maps, and it should not be ignored for countries like Korea. Furthermore, a DEM was obtained from the ALOS PALSAR sensor with a high resolution of 12.5 m, which is an advantage of the presented method because this map helps extract the seven basic maps including elevation, slope, flow accumulation, TWI, drainage network, plan curvature, and distance from the river with high resolution for additional use in subsequent map development and analysis. Another novelty of this study is the use of basic maps without human intervention, which are appropriate for ungauged areas. In addition, each map was reclassified into five categories to determine the direct impact of each part on the flood-prone area, which is another advantage of the proposed method. Finally, the number of pixels related to each polygon extracted from the real flood map was counted to assess the accuracy and validation of the flood susceptibility map obtained from this study. The results indicate that the proposed method shows the high performance of identifying the flood prone area. The region was classified into six classes in terms of flood susceptibility by integration of ArcGIS and the study area.
This study is a new approach for obtaining flood susceptibility map, which focused on: (1) preparation of basic maps with high resolution of satellite images; (2) extraction of a new land use map for the year of 2020 to clarify the growth in urban area with all updated maps for the same year of 2020 without direct human intervention; (3) validation of the extracted flood susceptibility map with the real flood map by counting the number of pixels in flooded areas. Therefore, it should be noted that the proposed method just considers satellite images, ArcGIS, and AHP for obtaining a flood susceptibility map, which can be used in any other geographical area with available satellite images. It can help decision-makers and planners to perform suitable flood management in the future, especially in the identification and analysis stages of risk assessment projects.

Data Availability Statement:
The data presented in this study are available on request from the first author.

Conflicts of Interest:
The authors declare no conflict of interest.