Next Article in Journal
Prediction of Mean Wave Overtopping Discharge Using Gradient Boosting Decision Trees
Previous Article in Journal
Effect of Cobalt, Cadmium and Manganese on Nitrogen Removal Capacity of Arthrobacter arilaitensis Y-10
Previous Article in Special Issue
Improved Landsat-Based Water and Snow Indices for Extracting Lake and Snow Cover/Glacier in the Tibetan Plateau
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integration of GIS-Based Multicriteria Decision Analysis and Analytic Hierarchy Process to Assess Flood Hazard on the Al-Shamal Train Pathway in Al-Qurayyat Region, Kingdom of Saudi Arabia

1
Research Center, Ministry of Housing, Riyadh 11461, Saudi Arabia
2
Geography Department, College of Arts, Princess Nourah bint Abdulrahman University, Riyadh 84428, Saudi Arabia
3
The High Institute of Tourism, Hotels and Computer (H.I.T.H.C.), The Ministry of Higher Education and Scientific Research (MHESR), Alexandria 21500, Egypt
4
Department of Geography, Faculty of Art, Monofiya University, Shbeen El Koom 32511, Egypt
*
Author to whom correspondence should be addressed.
Water 2020, 12(6), 1702; https://doi.org/10.3390/w12061702
Submission received: 13 May 2020 / Revised: 9 June 2020 / Accepted: 11 June 2020 / Published: 14 June 2020
(This article belongs to the Special Issue Applications of Remote Sensing and GIS in Hydrology II)

Abstract

:
Understanding the dynamics of floods in dry environments and predicting an accurate flood hazard map considering multiple standards and conflicting objectives is of great political and planning importance in the Kingdom of Saudi Arabia’s vision for the year 2030, in order to reduce losses in lives, property, and infrastructure. The objectives of this study are (1) to develop a flood vulnerability map identifying flood-prone areas along the Al-Shamal train railway pathway; (2) to forecast the vulnerability of urban areas, agricultural land, and infrastructure to possible future floods hazard; and (3) to introduce strategic solutions and recommendations to mitigate and protect such areas from the negative impacts of floods. In order to achieve these objectives, multicriteria decision analysis based on geographic information systems (GIS-MCDA) is used to build a flood hazard map of the study area. The analytic hierarchy process (AHP) is applied to extract the weights of eight criteria which affect the areas which are prone to flooding hazards, including flow accumulation, distance from the wadi network, slope, rainfall density, drainage density, and rainfall speed. Furthermore, the receiver operating characteristic (ROC Curve) method is used to validate the presented flood hazard model. The results of the study reveal that there are five degrees of flooding hazard along the Al-Shamal train path, ranging from very high to very low. The high and very high hazard zones comprise 19.2 km along the path, which constitutes about 26.45% of the total path length, and are concentrated at the intersections of the Al-Shamal train pathway with the Bayer and Al-Makhrouk wadis. Moderate, low, and very low flood severity areas constitute nearly 53.39 km, representing 73.55% of the total length (72.59 km) of the track. These areas are concentrated at the intersection of the Al-Shamal train track with the Haseidah Al-Gharbiyeh and Hsaidah Umm Al-Nakhleh wadis. Urban and agricultural areas that are vulnerable to high and very high flooding hazards are shown to have areas of 29.23 km2 (22.12%) and 59.87 km2 (46.39%), respectively.

1. Introduction

Floods are considered to be the most common and frequent natural hazards [1]. In recent decades, flood hazards have increased throughout the world. The severity of a flood is directly related to several factors, including population growth in urban areas and expansion in economic and social infrastructure. Furthermore, changes in land-use play an important role in the hydrological behavior of drainage basins, the hydrological cycle, and the morphology of wadis, which lead to an increase in flood hazard [2]. A wadi is a natural watercourse depression on the earth’s surface, which is dry except during periods of rainfall.
The occurrence of extreme climate changes around the world have a major impact on the frequency of occurrence of floods, as well as their dynamics, rapidness, and destructivity [3]. Climate changes affect the amount of water falling into the wadi network and, consequently, the amount of runoff. Other geographical factors also affect the severity of the floods, including the topography, geomorphology, and morphometry.
In recent years, flood hazard maps have represented a major and strategic component of flood mitigation. Flood hazard maps aim to provide residents with information on probable damages and potential disaster prevention [4]. Flood hazard maps are considered nonstructural representations for the mitigation of the potential impacts of floods. These maps can be produced in various forms, such as maps for flood tracking, emergency action plans, flood information, and flood hazard reduction, as well as for various other purposes [5].
Mapping of natural hazards and hazard analysis along road pathways and railway tracks involves several criteria [6,7]. Geographic information system (GIS) and remote sensing (RS) technologies have made significant contributions to the analysis of natural hazards [8,9,10]. Over the past few decades, researchers have developed various methods and models for natural hazard mapping, including GIS-RS [11,12,13], frequency rate [14,15], analytical hierarchy process [16], fuzzy logic [17], logistic regression [18], artificial neural networks [19,20,21], weights-of-evidence [22], multicriteria decision [23], support vector machine models [24], adaptive neuro-fuzzy inference system (ANFIS) [25,26], biogeography-based optimization [27], decision tree [28], multivariate adaptive regression splines [29], and hydraulic modeling techniques, which are considered essential tools for flood hazard management and mitigation [30,31,32].
This study combines the multicriteria analysis based on geographic information systems (GIS-MCDA) approach and the analytic hierarchy process (AHP). Multicriteria analysis (MCA) is considered an important tool in the analysis of complex decision problems, which often involve noncomparable criteria or data [33]. The AHP is the most common method for multicriteria decision-making [34,35].
Olga et al. [36], Nerantzis et al. [37], and Somaiyeh and Mehran [38] discussed the mapping of flood hazard using multicriteria analysis integrated with GIS, RS, and the AHP, respectively. These various studies identified six factors which affect the occurrence of floods, including the cumulative runoff, distance from the water drainage network, elevation, land-uses, rainfall density, and geology.
The main concern of this study is to shed light on the flood hazards in the study area, concentrating on the impacts of five main wadis (from west to east and from north to south), namely: Wadi Al-Makhrouk, which influences Hadithah village to the north of Qurayyat city and strikes the railway in the north; Wadi Haseidah Al-Gharbiyeh, which affects the north of Al-Qurayyat city; Wadi Haseidah, Umm Nakhleh, which influences the southern part of Al-Qurayyat city; Wadi Bayer, which threatens the villages of Bayer and Ghati, and covers the railway track part facing them; and, finally, Wadi Sarmada, which impacts on the villages of Al-Nasifah, Skullim, and Al-Rudaifa.
The flow of these wadis have severe negative impacts, resulting in the occurrence of frequent annual floods. Table 1 presents the records of historical flooding events which occurred in the study area. The floods that occurred on 10 November 2018 were the most severe, which resulted in halting traffic completely, the flooding of houses and schools, and the partial destruction of Al-Qurayyat hospital. The Ghati, Bayer, and Hadithah villages were flooded and their people were displaced, as can be seen in Figure 1.
The topographical nature of the drainage basins of the study area worsens the flood impacts. The study area is characterized by its varied topographic features and slopes, which intersect with the dense networks of dry wadis. These characteristics negatively affect the study area and make it permanent vulnerable to flood hazards, in addition to the increased unplanned urban growth centralized in the pathways and estuaries of dry wadis, which are considered to be the most vulnerable to flood hazards. The construction of the train pathway in 2017 and its crossing near the cities and villages of Al-Qurayyat led to changes in the land-use, the morphology of the wadis, and increased the flood hazard in the region. The importance of this study lies in the scarcity of applied studies which address flood hazards in the drainage basins of the northern region in the Kingdom of Saudi Arabia (KSA). Furthermore, this study determines the negative impacts of floods in Al-Qurayyat city and in the railway of the Al-Shamal train pathway facing it. Although railways are of strategic and economic importance in KSA, most of these railways are subject to geomorphological, hydrological, and flood hazards due to the varied climates and topographies that they pass through. Only two studies have focused on railways in KSA, where the first study addressed the relationship between the effect of land-use changes and the increase of flood hazard due to the Riyadh–Dammam train in the Eastern region in the time period from 2011–2017 [2]; meanwhile, the second study concentrated on the Al-Shamal train pathway in Riyadh and examined the impacts of some of the morphological characteristics of the region on it [39]. The main aims of this study can be summarized in five items: (1) evaluate the flood hazard on the Al-Shamal train pathway in the Al-Qurayyat region using multicriteria decision analysis based on geographic information systems (GIS-MCDA) and hierarchical analysis (AHP); (2) identify areas vulnerable to flood hazard along the northern railway track and the hazard levels for each part thereof; (3) predict the urban areas, agricultural lands, and infrastructure prone to possible future flood hazards and the degree of hazard in each; (4) reduce the flood hazard and adapt to serve future planning and decision-making in the study area; and (5) propose suitable strategic solutions to mitigate and protect against the negative impacts which are expected due to flood hazards in the study area.

2. The Study Area

The administrative boundary of the study area is in the north of the Al-Jawf region, which is located in the north of KSA. The main drainage basins that affect Al-Qurayyat city and the Al-Shamal train pathway facing it are located on the international border between the KSA and the Hashemite Kingdom of Jordan. Al-Qurayyat city is about 1200 km from the capital of Riyadh, about 310 km from Sakaka city (the headquarters of the Emirate of Al-Jawf region), and about 350 km from Arar city (the headquarters of the Emirate of the northern region). The study area is bounded to the north and west by the Hashemite Kingdom of Jordan, to the east by Tarif Governorate, and to the south by the Al-Jawf region, as shown in Figure 2.
The study area is considered an ideal drainage environment for the drainage basins located in the southeast of the Hashemite Kingdom of Jordan, especially the main drainage basins mentioned in this study. The drainage basins (wadis Al-Makhrouk, Haseidah Al-Gharbiyeh, Haseidah Umm Nakhleh, Bayer, and Sarmada) affecting Al-Qurayyat city and the Al-Shamal train pathway, from the north to the south, are located between the latitudes of 30°16′9.75” and 31°31′33.01” N, and the longitudes of 36°30′45.71” and 37°48′52.18” E. The study area is shown in Figure 3.

3. Methodology and Data Processing

The key goal of using MCDA is to “study the selection of a number of multiple criteria and conflicting objectives”. The MCDA method permits the assessment of a region based on multiple objectives and criteria and supports decision-making in the flood hazard assessment process. Figure 4 presents the methodology steps which are followed to obtain the final flood hazard map.

3.1. Data Collection

Multiple data types and sources were used to conduct the multicriteria analysis, in order to determine the degree of flood hazard. All data were standardized for integration within the GIS environment. Table 2 shows the collected data used for the study and presents their characteristics, including the spatial resolution and sources. A digital elevation model (DEM) with spatial resolution of 12.5 m was downloaded freely from http://vertex-retired.daac.asf.alaska.edu, in addition to topographic maps (of scale 1:50,000) obtained from the Survey Public Authority. Landsat 8 OLI/TIRS satellite images with 30 m spatial resolution were acquired in January 2020, downloaded freely from http://earthexplorer.usgs.gov. A geological map (of scale 1:250,000) was obtained from the Saudi Geological Survey. Rainfall rate (mm day−1) data were downloaded from the NASA Godard Earth Sciences Data and Information Service Center (GES DISC GSFC).

3.2. Derivation of the Study Criteria Layers

The collected data were used to extract the flow accumulation, stream order, and drainage basins through hydrological analysis. The slope was derived through surface analysis of the DEM. Land cover was extracted from Landsat 8 OLI/TIRS images.
The hydrological soil groups in the study area were determined, according to the soil classification groups in the Soil Conservation Service (SCS) method, by analyzing the geological map. The runoff speed was derived from the DEM using a hydraulic modeling program (HEC-RAS). Figure 5 and Figure 6 show the MCDA criteria used in the study.

3.3. Processing and Classifying of Criteria Layers

The processing and classification of criteria layers used in this study are summarized below:

3.3.1. Standardization of Criteria Layers

All derived raster layers, including flow accumulation, slopes, rainfall density, and land cover were standardized using a raster calculator in map algebra in ArcGIS, where the Euclidean distance was implemented, using the wadi network vector layer, in order to determine the distance from the wadi network. The other vector layers, including drainage density, runoff speed, and the hydrological group, were converted to raster mode in ArcGIS.

3.3.2. Reclassification of Criteria Layers

The eight extracted criteria layers were reclassified in ArcGIS, in order to determine the degree of severity and to create the hazard maps. The reclassify tool (in the spatial analysis toolbox in ArcGIS) was applied to determine and define the hazard classes and ranks for each criterion. Hazard classes and ranks were assigned for each criterion using the experience of the authors and extensive previous studies. For example, the distance from wadis criterion was classified as: a distance between 0 to 0.28 km from the wadis had the highest hazard (rank = 10), while a distance between 3.37 to 5.51 km from the wadis had the lowest hazard (rank = 1). Furthermore, with respect to the flow accumulation criterion, the class with the lowest flow had the lowest hazard and was assigned rank = 1, and vice versa. Table 3 shows the hazard ranks for each criterion.

3.3.3. Application of AHP

AHP is an important tool in multicriteria decision-making. AHP is a mathematical theory of measurement developed by Thomas Saaty and which has been proven successful in many different fields, including management and economics. AHP is the best hierarchical framework to derive the importance of each criterion, in relation to the other criteria corresponding to it. Table 4 illustrates the importance of AHP criteria, as presented by Thomas Saaty.
The production of a flood hazard map using AHP can be summarized in five stages, as presented below:
i. First Stage: Calculate the Importance Values of the Criteria
In this stage, the importance (priority) values of each criterion relative to the other criteria are established, according to the previous table. The assignment of the importance values in this study involved evaluating each criterion against all other criteria at a hierarchical level. Table 5 shows the comparison between the eight criteria in this study and the assigned importance values.
ii. Second Stage: Calculate the Percentage of Importance Values
The percentage of importance values are calculated between every two criteria (criterion in row and criterion in a column) using Equation (1):
a ¯ j k = a j k 1 = I m a I k ,
where:   a ¯ j k is the percentage of importance value between two criteria,   a j k is the importance value between two criteria (one a row and the other a column), and 1 = I m a I k is the total of the standardized columns. The relative weight value of each of the row criteria is determined by Equation (2):
w j   = 1 = I m a j l m ,
where:   w j is the value of the relative weight of the standardized rows, 1 = I m a j l is the sum of percentages of importance values for a criterion row, and m is the final value of 1 = I m a j l for all rows.
Table 6 shows the percentages of importance values using the AHP.
iii. Third Stage: Create the Weight Values Matrix
Table 7 and Table 8 show the matrices of the weight values, according to the priority values between the criteria.
where:   ʎ M a x is the largest eigenvalue of the matrix
iv. Fourth Stage: Calculate Consistency index (Consistency Verification)
A matrix is said to be consistent if, for each two values that are reciprocal to each other, the multiplication of the two values is equal to one, and if the elements in the column are the reciprocal of the elements in the corresponding row. The consistency index is calculated using Equation (3):
C I = ʎ M a x n n 1 ,
where:   ʎ M a x is the largest eigenvalue of the matrix (as seen in the previous table), which is always greater than or equal to the number of rows or columns; and n is the number of criteria.
C I = 8.218 8 8 1 = 0.03
The “CI” value is more acceptable if it is closer to zero. The closer the CI value is to zero, the more confidence we can have in the consistency index. To the contrary, the farther the CI value from zero, the more inconsistent. The consistency ratio is calculated using Equation (5):
c o n s i s t e n c y   r a t i o = C I R ,
where: R is the random index, according to the number of criteria; its value is determined from Table 9.
From the previous table, the value of the random index in this case was (R) = 1.4, as we used eight criteria in this study. The consistency ratio was   0.03 1.4 = 2% = 0.02, which is an acceptable percentage, as the consistency ratio was within 0.1 (10%). A consistency ratio exceeding 0.1 indicates an inconsistency (contradiction). Table 10 demonstrates the weights of the criteria using the AHP.
v. Final Stage: Production of the Flood Hazard Map
This stage is the prefinal part of the study, which involves the creation of the structural model constructed within the ArcGIS environment to extract all the needed criteria automatically by adding the required ArcGIS modules to get the flood hazard map [40]. The model included the following: (1) the distance maps were generated to produce the proximity distance between different criteria by using the ArcCatalog, the Spatial Analysis toolbox, and the Euclidean Distance. (2) The slope criterion was used to create a slope-surface map. (3) All criteria maps were reclassified to define the categories of spatial suitability for each criterion. (4) As a subsequent step, all criteria layers were clipped using the ‘Extract by mask’ tool to match the boundary of the study area. (5) Then, the Raster Calculator was applied to overlay all the criteria layers according to the obtained AHP weights specified for each criterion (weighted overlay). (6) The final step was to run the structure model to get the flood hazard map in five classes: very high, high, moderate, low, and very low. The structural model of flood hazard on the Al-Shamal train pathway created by integrating GIS-MCDA and AHP is presented in Figure 7.

4. Results and Discussion

4.1. Analysis of Criteria Affecting the Flood Hazard of the Al-Shamal Train Pathway

The eight criteria affecting the occurrence of floods along the Al-Shamal train pathway in the study area were distance from the wadi network, slope, land cover, hydrological soil groups, flow accumulation, rainfall density, drainage density, and runoff speed. By studying and analyzing these criteria, the extent of the impact and the spatial distribution of floods were determined. More details of the analysis of these criteria are presented below.

4.1.1. Distance from the Wadi Network

The water runoff flow in the wadis constitutes a major factor in the occurrence of a floods, where floods occur due to the immersion of the wadi network in a drainage basin. Such inundation increases in areas near the wadis (which are more vulnerable to floods) and decreases in areas far from them (which are less vulnerable to floods) [36,41,42,43]. In this study, the distance from the wadi network was divided into ten classes, where the closest class to the wadi network was less than 0.28 km. This class was the most vulnerable and concentrated around the five wadi networks in the study area. The least vulnerable areas were at a distance of more than 3.37 km, which were concentrated in the estuaries of the Bayer and Al-Makhrouk wadis, as shown in Table 11 and Figure 8 and Figure 9.

4.1.2. Slope

There is a strong positive correlation between the slope of the surface and the runoff speed. Regions with low inclines are exposed to high runoff speeds, causing the occurrence of severe floods [36,42,43,44]. The slope criterion was divided into ten classes. Areas with a slope less than 1.07° (the lowest slope) were considered to have the highest hazard, which were concentrated in the northern and western parts of the study area, in the Sarmada, Haseidah Al-Gharbiyeh, and Haseidah Umm Nakhleh drainage basins; while areas with a slope greater than 20.34° (the highest slope) were considered to be the least vulnerable to floods, which were concentrated in the center of the study area, around the wadi networks of the Bayer and Al-Makhrouk drainage basins, as shown in Table 12 and Figure 8.

4.1.3. Land Cover

The land cover determines the areas which are vulnerable to flood hazards. Urban areas and agricultural lands are considered less vulnerable to floods (the lowest hazard), as they obstruct and consequently decrease the runoff speed, in addition to the presence of branched irrigation and drainage channels in such areas. In contrast, wadis, swamps, and desert areas are the most vulnerable to floods and have the highest hazard [36,42,43,44].
The land cover criterion was divided into five classes. The urban areas and agricultural lands concentrated in Al-Qurayyat city to the north of the study area were considered the least vulnerable to floods (lowest hazard). Wadis, swamps, and desert areas, representing the majority of the study area, had the highest flood hazard, as shown in Table 13 and Figure 8.

4.1.4. Hydrological Soil Groups

The hydrological soil group signifies the ability of specific soil types to reduce flood hazards. The hydrological group criterion in this study was divided into four groups: A, B, C, and D. Group A represents deep sandy soils with very high intrusion rates, which minimize the flood hazard. The hydrological group (B) includes relatively fine-grained soils with moderate intrusion rates. The hydrological group (C) includes fine-grained soils with low intrusion rates. The hydrological group (D) indicates cohesive soil with very fine grains and very low intrusion rates, which increase the flood hazard [45,46]. Groups A and B constitute low flood hazard and covered a large part of the study area. Groups C and D represent high flood hazard and were mostly situated in the northeast and in scattered areas in the center of the study area. The areas of the four hydrological soil groups are shown in Table 14 and Figure 8.

4.1.5. Flow Accumulation

Flow accumulation is one of the most important criteria in determining areas vulnerable to flood hazard. High values of flow accumulation indicate potential high flood hazard, whereas low values indicate low flood hazard [36,37]. The flow accumulation map of the study area was divided into five categories: (1) the first flow accumulation class (greater than 700,033 m3) had the highest flood hazard and was concentrated in the estuaries of the large wadis located in Bayer and Sarmada; (2) the second flow accumulation class (390,018–700,033 m3) was situated in the estuaries of the Bayer, Sarmada, Al-Makhrouk, and Haseidah Al-Gharbiyeh wadis; (3) the third flow accumulation class (180,008–390,018 m3) was located in some tributaries of all wadis, except for that of Umm Nakhla; (4) the fourth flow accumulation class (45,002–180,008 m3) was concentrated in some tributaries of all wadis; and (5) the fifth flow accumulation class (less than 45,002 m3) represented the land between the wadis. The areas of the flow accumulation categories are shown in Table 15 and Figure 10.

4.1.6. Rainfall Density

Rain is the main element for flood occurrence, where areas with high rainfall density are the most vulnerable to floods and areas with low rainfall density are the least susceptible to floods [35,36,47].
The rainfall density criterion was divided into ten classes. The class corresponding to rainfall greater than 14.3 mm had the highest flood hazard, which was concentrated in the west of the study area. The rainfall density class with less than 8.96 mm had the lowest flood hazard, which was concentrated in the southeast of the study area, as shown in Table 16 and Figure 10.

4.1.7. Drainage Density

Drainage density is one of the factors affecting the flood risk, as areas with high drainage density are more susceptible to floods and, therefore, have higher flood hazard than areas with less drainage density [45,46]. The drainage density criterion was divided into five classes, where areas with a drainage density greater than 3.91 km/km2 had the highest flood hazard, which was represented in the Al-Makhrouk basin in the west of the study area. Areas with a drainage density of less than 3.8 km/km2 represented the lowest flood hazard, which was located in Haseidah Umm Nakhleh basin (whose estuary ends at Al-Qurayyat city), as shown in Table 17 and Figure 10.

4.1.8. Runoff Speed

Runoff speed is one of the factors influencing the severity of a flood, as areas with high runoff speed are the most susceptible to flood hazard, while areas with low runoff speed are the least susceptible to flood hazard [45]. The runoff speed criterion was divided into five classes, where runoff with a speed greater than 1.21 km/h indicated the highest flood hazard class, which was situated in the Bayer basin and in the center and west of the study area. The runoff speed class with a speed less than 1.12 km/h had the lowest flood hazard, which was located in the Haseidah Umm Nakhleh basin (whose estuary ends in Al-Qurayyat city) in the northeast of the study area, as presented in Table 18 and Figure 10.
Figure 11 shows the correlation relationships between each of the eight criteria which affect the occurrence of floods and the degree of flood hazard. There was a positive relationship between the degree of flood hazard and some criteria affecting the occurrence of floods. These criteria were the flow accumulation, rainfall density, drainage density, and runoff speed, where the degree of hazard increased with an increase in the value of these criteria. Meanwhile, there were negative relationships between the degree of flood hazard and some other criteria; namely, slope and distance from the wadi network, where the degree of severity of the flood hazard decreased with an increase in the value of these criteria. Overall, the wadi network and desert areas had high hazard and the hydrological group (D) had the highest hazard.

4.2. Analysis of the Flood Hazard

Five classes of flood hazard were identified for the Al-Shamal train pathway in Al-Qurayyat city in 2020, as illustrated in Figure 12 and Table 19.

4.2.1. Very High Hazard Areas

The very high flood hazard areas were concentrated around the dry wadis of the Bayer basin, which mediates the main drainage basin and the estuary of which ends at the agricultural lands to the southeast of Al-Qurayyat city; the Al-Shamal train railway; and the Al-Makhrouk basin, in the west of the main basin (its estuary ends at the north railway track). The percentage of hazard in these areas ranged between 68.2% and 91.2%, and the area of the high flood hazard zone was about 1744.22 km2; that is, 15.84% of the total area of the basin.

4.2.2. High Hazard Areas

The high flood hazard areas were situated around the dry wadis in the Hsaidah Al-Gharbiah basin, whose estuary ends at the north railway. The hazard percentage in these areas was between 61.6% and 68.2%, with a total area of nearly 3172.4 km2; that is, 28.81% of the total area of the basin.

4.2.3. Moderate Hazard Areas

The areas exposed to moderate flood hazards were located around the dry wadis in the Sarmada basin, which ends at the northern railway. The percentage of hazard in these areas ranged between 55.6% and 61.6%, with an area reaching 3184.52 km2; that is, 28.92% of the total area of the basin.

4.2.4. Low Hazard Areas

The low hazard areas were concentrated around the dry wadis in the Umm Nakhleh basin, whose estuary ends at Al-Qurayyat city and the northern railway. The total area was 2050.34 km2; that is, 18.62% of the total area.

4.2.5. Very Low Hazard Areas

The very low flood hazard areas were concentrated in the lands between the wadi networks of the Hsaidah Um Nakhla and Sarmada basins, the hazard percentage in these areas being between 24% and 48.5%, with an area of about 860 km2; that is, 7.81% of the total area.

4.3. Flood Hazards in the Urban and Agricultural Areas

The flood hazard map in Figure 13 reveals that the urban areas exposed to very high, high, medium, low, and very low flood hazards were 4.44, 24.79, 56.3, 32.65, and 13.97 km2, respectively, with percentages of 3.36%, 18.76%, 42.61%, 24.71%, and 10.57%, respectively, of the total urban area in the study area (which was nearly 132.14 km2). The higher flood hazards in the urban areas were concentrated in Ghati village, which is part of the major estuary of the Bayer wadi; the middle of the basin of Hadithah village (the major estuary of Al-Makhrouq wadi in the west of the basin Al-Qurayyat city). Al-Qurayyat city was subject to potential flood hazards with a moderate degree of severity, especially as it is considered part of the major estuary for the Haseidah wadi, in addition to the high hazard from the urban expansion between the center of the city and the Al-Qurayyat airport, which is considered part of the main estuary of the Hassidah Umm Nakhleh wadi.
Figure 13 and Table 20 show that the agricultural lands in the study area exposed to very high, high, moderate, low, and very low flood hazard had areas of 18.38, 41.49, 37.29, 19.1, and 12.8 km2, respectively, comprising 14.24%, 32.15%, 28.89%, 14.8%, and 9.92%, respectively, of the total agricultural land in the study area (for a total of 129.06 km2). The very high flood hazards in the agricultural lands were concentrated around the village of Ghati, which is considered part of the main estuary of the Bayer basin. Very high flood hazards also existed in the south of the village of Ghati and around Hadithah village (major estuary of the Al-Makhrouk wadi in the west of the drainage basin).
Figure 14 illustrates the degree of flood hazard in urban areas and agricultural lands, in addition to the displacement of residents of the villages of Bayer and Ghati.

4.4. Degrees of Floods Hazards along the Al-Shamal Train Pathway

Inspection of the flood hazard map showing the degree of severity on the Al-Shamal train pathway in the study area reveals that (1) the lengths of the sections exposed to very high, high, medium, low, and very low flood hazards corresponded to 6.96, 12.24, 6.54, 15.24, and 31.64 km, respectively; a percentage of 9.59%, 16.86%, 9.01%, 20.99%, and 43.55%, respectively, of the total pathway length; (2) the sections exposed to very high and high flood hazards were concentrated in front of Wadis Bayer, Al-Makhrouk, and Haseidah Umm Nakhleh at the villages of Ghati, Bayer, Hadithah, and Al-Nasifah; (3) the sections exposed to low and very low flood hazards were concentrated in the areas between the estuaries of the wadis, particularly in the section between Al-Qurayyat city and Al-Qurayyat airport in the northwest of the city and the city of Ghati to the southeast of the city. Table 21 and Figure 15 illustrate the degrees of severity on the Al-Shamal train pathway in the study area. Figure 16 shows a field visit to inspect the water drainage installations at the bottom of the Al-Shamal train pathway in the study area.

4.5. Validation of Flood Hazard Model

Validation of the accuracy and reliability of the results for decision-making is essential and crucial, as the flood hazard map is a key element for optimal decision-making to mitigate the negative impacts of potential flood hazards and for future planning decisions and strategies. There are many scientific methods for the validation of the results of the flood hazard map. The receiver operating characteristic (ROC) curve and the area under the curve (AUC) are considered to be the most popular in flood hazard studies [48,49].
The ROC curve is a mathematical and statistical method which can be applied to examine the accuracy and reliability of the results of a flood hazard map [50,51,52]. The ROC curve is a graph with two axes (x,y), where the false-positive rate (FPR) is displayed along the x-axis to express the number of samples in the flood hazard map which are unlikely to occur and, in fact, where floods have not occurred. The true positive rate (TPR) is displayed along the y-axis, which expresses the number of samples in the flood hazard map which are likely to occur and, in fact, where floods have occurred.
The reliability of the prediction rate can be achieved by calculating the area under the curve (AUC). If the AUC value is less than 0.5, this means that the flood hazard map is unreliable and that its results are misleading and incorrect. However, if the AUC value is 1, then the results of the flood hazard map are reliable. Table 22 shows the values of the AUC and the degree of validation of the results.
The ROC curve in Figure 17 was drawn using the SPSS program, considering the flood hazards value from the samples collected from the field study and the values extracted from the flood hazard map. The AUC value (using the SPSS program results) was 0.971 (AUC = 0.971) with a 95% confidence level. This AUC value indicates that the reliability of the flood hazard map produced in this study is excellent (around 97.1%).

5. Conclusions

Understanding the dynamics of floods and obtaining an accurate flood hazard map comprise major political, planning, and popular demands in KSA. Such a map, in turn, reduces and mitigates flood hazards and assists urban and regional planning authorities in developing optimal protection measures. The study area considered in this paper is one of the largest dry environments in the world, characterized by the sudden occurrence of floods; almost every year due to the presence of five main wadis. Flood hazards in the study area involve damage to infrastructure, urban areas, and transportation, such as roads and railways.
The main aim of this scientific research was to carefully study and inspect the flood hazard in the study area after the construction of the Al-Shamal train pathway in 2017 near the cities and villages of Al-Qurayyat, considering that there were several historical flooding events in the study area which flooded urban communities close to the train pathway, producing significant damage to lives, property, and infrastructure. Thus, the flood hazard map presented in this study can help in developing an appropriate strategic solution by the responsible authorities and policymakers, shedding light on the hot-spot flood areas in order to mitigate and reduce the hazard.
Integrating GIS-MCDA and AHP in flood hazard studies is of great importance in the case of multiple criteria and conflicting objectives. The results of this study revealed that the very high and high flood hazard zones constituted about 44.65% of the study area, concentrated in the Bayer wadi and Al-Makhrouq wadi drainage basins. These zones represent a severe hazard for the urban areas, agricultural lands, and the Al-Shamal train pathway. The moderate, low, and very low flood hazard areas constituted about 55.35% of the study area, which were concentrated in the drainage basins of wadis Haseidah Al-Gharbiyeh, Haseidah, Umm Nakhleh, and Sarmada. About 22.12% of the urban areas were exposed to very high and high flood hazards, as concentrated in Ghati, Bayer, and Hadithah villages. Meanwhile, the urban areas in the center and the northern border of Al-Qurayyat city were exposed to about 77.88% of average, low, and very low hazards. About 26.45% of the Al-Shamal train pathway was exposed to high and very high flood hazards, while about 73.55% of the train pathway was exposed to moderate, low, and very low flood hazards.

6. Recommendations

Flood hazard prevention projects, such as drainage basins, barriers, channels, dams, collecting basins, water harvesting, and so on, should be planned in the upper or middle sectors of the wadis, and not in the areas targeted by the flood hazards. The projects that are currently established to mitigate the frequent floods every year are not commensurate with the nature of the floods in the Al-Qurayyat region. The presence of five main wadis, originating from west to east and carrying a huge amount of floodwater, exposes urban and agricultural areas to loss every year without taking advantage of them, in addition to the destruction of infrastructure.
It was noticed, during the field visits, that there are few implemented projects to protect the study area from flood hazards. However, these projects should be focused on the Haseidah Al-Gharbiyeh, Haseidah Umm Nakhleh, Bayer, and Sarmada wadi areas, in order to prevent floods in Al-Qurayyat, its villages, and the Al-Shamal train pathway.
The study area urgently needs to develop integrated strategies to mitigate flood hazards. These strategies should involve studying the wadis affecting the study area, including wadis Al-Makhrouk, Haseidah Al-Gharbiyeh Umm Nakhleh, Bayer, and Sarmada. Moreover, hydrological studies must include the maximum floodwater flow in these wadis in addition to the estimated floodwater amount and arrival time, based on the amount of rain falling on these basins during the different return periods every 50 and 100 years. Furthermore, setting up floodwater drainage facilities is required, according to these recent studies, in order to establish suitable protection measures.
We recommend the implementation of urgent plans to mitigate the hazards of floods. Priority of these plans has to be given to areas of high and very high flood hazards, which are represented by the Wadi Bayer and Wadi Al-Makhrouk drainage basins; specifically the point at which these wadis intersect the Al-Shamal train pathway. The areas subject to frequent flood hazards which have exposed the residents of Bayer, Ghati, and Hadithah villages to full displacement in recent years require special interest.
We recommend the establishment of early warning systems (EWS) through the building of permanent stations to measure and monitor the water flow in the five wadis in the study area. EWS enable the monitoring of floods automatically, sending warnings to the monitoring stations when the water level reaches a dangerous level.
We recommend directing future urban and infrastructure planning in the study area away from the high and very high hazard areas in the floods hazard map obtained in this study. Furthermore, either construction in the high and very high hazard areas should be prohibited, or the encroachment of agricultural areas should only be allowed after the provision of mechanisms to mitigate flood hazards which are commensurate with the nature of the wadis.
We also recommend developing long-term plans to deal with the critical areas in the study area, especially in the Faydat Al-Rashrashyia region in the north of Al-Qurayyat city, the northwest region of Al-Qurayyat city along the new international road, and the villages of Aqeelah, Daabousiya, Bayer, Ghati, and Ain Hawas.
Finally, we recommend that the Saar company (the owner of the Al-Shamal train pathway project) should perform a site suitability analysis, in order to study and evaluate the current location of the Al-Shamal train pathway (2750 km) from Riyadh city in the south to Hadithah village in the north; in addition to intensifying studies on such topics as sand migration, floods, and geomorphological hazards, in order to ensure the safety of life and property.

Author Contributions

A.A. designed the study, developed the research idea and planned the research activities, designed the methodology, replied to the reviewer’s comments; wrote the manuscript, and professionally edited the manuscript; S.S.A.-A. performed climatic and statistical criteria analysis; H.M.A. contributed the Geological and hydrological criteria; S.A.M. developed the research methodology and wrote and reviewed the manuscript; I.I.A. and I.Y.I. carried out statistical and spatial analyses for monitoring the changes in land-use and provided valuable comments in writing this paper; and A.A., S.S.A.-A., H.M.A., S.A.M., I.I.A., and I.Y.I. carried out the research, including collecting the input data, preparing the manuscript, and carrying out a statistical analysis of the obtained results. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University, through the Research Groups Program Grant No. (RGP-1440-0029).

Acknowledgments

The authors highly appreciate the great support from the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University, Grant No. (RGP-1440-0029).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Termeh, S.V.R.; Kornejady, A.; Pourghasemi, H.R.; Keesstra, S. Flood susceptibility mapping using novel ensembles of adaptive neuro fuzzy inference system and metaheuristic algorithms. Sci. Total Environ. 2018, 615, 438–451. [Google Scholar] [CrossRef] [PubMed]
  2. Abdelkarim, A.; Gaber, A.; Alkadi, I.; Alogayell, H. Integrating remote sensing and hydrologic modeling to assess the impact of land-use changes on the increase of flood risk: A case study of the Riyadh–Dammam train track, Saudi Arabia. Sustainability 2019, 11, 6003. [Google Scholar] [CrossRef] [Green Version]
  3. Brunner, M.I.; Sikorska-Senoner, A.E.; Seibert, J. Bivariate analysis of floods in climate impact assessments. Sci. Total Environ. 2018, 1392–1403. [Google Scholar] [CrossRef] [PubMed]
  4. Abdel Karim, A. Assessment of the expected flood hazards of the Jizan-Abha highway, Kingdom of Saudi Arabia by integrating spatial-based hydrologic and hydrodynamic modeling. Glob. J. Res. Eng. Glob. J. 2019, 19, 1024. Available online: https://engineeringresearch.org/index.php/GJRE/article/view/1948 (accessed on 8 April 2020). [CrossRef] [Green Version]
  5. Abdel Karim, A.; Awawdeh, M.M.; Alogayell, H.M. Integration of remote sensing, hydrological and hydraulic modeling for flood risk assessment and mitigation in the coastal city of Al-Lith, Saudi Arabia. Int. J. GEOMATE 2020, 18, 252–280. [Google Scholar] [CrossRef]
  6. Xu, C.; Chen, Y.; Chen, Y.; Zhao, R.; Ding, H. Responses of surface runoff to climate change and human activities in the arid region of Central Asia: A case study in the Tarim River Basin, China. Environ. Manag. 2013, 51, 926–938. [Google Scholar] [CrossRef]
  7. Poussin, J.K.; Botzen, W.W.; Aerts, J.C. Factors of influence on flood damage mitigation behaviour by households. Environ. Sci. Policy 2014, 40, 69–77. [Google Scholar] [CrossRef]
  8. Patel, D.; Srivastava, P.K. Flood Hazards Mitigation Analysis Using Remote Sensing and GIS: Correspondence with Town planning scheme. Water Resour. Manag. 2013, 27, 2353–2368. [Google Scholar] [CrossRef]
  9. De Moel, H.; Van Vliet, M.; Aerts, J.C.J.H. Evaluating the effect of flood damage-reducing measures: A case study of the unembanked area of Rotterdam, the Netherlands. Reg. Environ. Chang. 2013, 14, 895–908. [Google Scholar] [CrossRef]
  10. Mohamed, S.A.; El-Raey, M.E. Vulnerability assessment for flash floods using GIS spatial modeling and remotely sensed data in El-Arish City, North Sinai, Egypt. Nat. Hazards 2019, 102, 707–728. [Google Scholar] [CrossRef]
  11. Mohamed, S.A. Application of satellite image processing and GIS-Spatial modeling for mapping urban areas prone to flash floods in Qena governorate, Egypt. J. Afr. Earth Sci. 2019, 158, 103507. [Google Scholar] [CrossRef]
  12. Althuwaynee, O.F.; Pradhan, B.; Park, H.-J.; Lee, J.H. A novel ensemble bivariate statistical evidential belief function with knowledge-based analytical hierarchy process and multivariate statistical logistic regression for landslide susceptibility mapping. Catena 2014, 114, 21–36. [Google Scholar] [CrossRef]
  13. Van Westen, C.; Rengers, N.; Soeters, R. Use of geomorphological information in indirect landslide susceptibility assessment. Nat. Hazards 2003, 30, 399–419. [Google Scholar] [CrossRef]
  14. Lee, M.-J.; Kang, J.-E.; Jeon, S. Application of frequency ratio model and validation for predictive flooded area susceptibility mapping using GIS. In Proceedings of the Geoscience and Remote Sensing Symposium (IGARSS), Munich, Germany, 22–27 July 2012; pp. 895–898. [Google Scholar] [CrossRef]
  15. Tehrany, M.S.; Pradhan, B.; Mansor, S.; Ahmad, N. Flood susceptibility assessment using GIS-based support vector machine model with different kernel types. Catena 2015, 125, 91–101. [Google Scholar] [CrossRef]
  16. Stefanidis, S.; Stathis, D. Assessment of flood hazard based on natural and anthropogenic factors using analytic hierarchy process (AHP). Nat. Hazards 2013, 68, 569–585. [Google Scholar] [CrossRef]
  17. Pradhan, B. Use of GIS-based fuzzy logic relations and its cross application to produce landslide susceptibility maps in three test areas in Malaysia. Environ. Earth Sci. 2010, 63, 329–349. [Google Scholar] [CrossRef]
  18. Pradhan, B. Flood susceptible mapping and risk area delineation using logistic regression, GIS and remote sensing. J. Spat. Hydrol. 2010, 9, 1–18. [Google Scholar]
  19. Kia, M.B.; Pirasteh, S.; Pradhan, B.; Mahmud, A.R.; Sulaiman, W.N.A.; Moradi, A. An artificial neural network model for flood simulation using GIS: Johor River Basin, Malaysia. Environ. Earth Sci. 2011, 67, 251–264. [Google Scholar] [CrossRef]
  20. Lohani, A.K.; Goel, N.; Bhatia, K. Improving real time flood forecasting using fuzzy inference system. J. Hydrol. 2014, 509, 25–41. [Google Scholar] [CrossRef]
  21. Samanta, S.; Pal, D.K.; Lohar, D.; Pal, B. Interpolation of climate variables and temperature modeling. Theor. Appl. Clim. 2011, 107, 35–45. [Google Scholar] [CrossRef]
  22. Tehrany, M.S.; Pradhan, B.; Jebur, M.N. Flood susceptibility mapping using a novel ensemble weights-of-evidence and support vector machine models in GIS. J. Hydrol. 2014, 512, 332–343. [Google Scholar] [CrossRef]
  23. Malczewski, J. GIS-based multicriteria decision analysis: A survey of the literature. Int. J. Geogr. Inf. Sci. 2006, 20, 703–726. [Google Scholar] [CrossRef]
  24. Tehrany, M.S.; Pradhan, B.; Jebur, M.N. Flood susceptibility analysis and its verification using a novel ensemble support vector machine and frequency ratio method. Stoch. Environ. Res. Risk Assess. 2015, 29, 1149–1165. [Google Scholar] [CrossRef]
  25. Bui, D.T.; Panahi, M.; Dou, J.; Singh, V.P.; Shirzadi, A.; Chapi, K.; Khosravi, K.; Chen, W.; Panahi, S.; Li, S.; et al. Novel hybrid evolutionary algorithms for spatial prediction of floods. Sci. Rep. 2018, 8, 15364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Bui, D.T.; Khosravi, K.; Li, S.; Dou, J.; Panahi, M.; Singh, V.P.; Chapi, K.; Shirzadi, A.; Panahi, S.; Chen, W.; et al. New hybrids of ANFIS with several optimization algorithms for flood susceptibility modeling. Water 2018, 10, 1210. [Google Scholar] [CrossRef] [Green Version]
  27. Ahmadlou, M.; Karimi, M.; Alizadeh, S.; Shirzadi, A.; Parvinnejhad, D.; Dou, J.; Panahi, M.; Parvinnejad, D. Flood susceptibility assessment using integration of adaptive network-based fuzzy inference system (ANFIS) and biogeography-based optimization (BBO) and BAT algorithms (BA). Geocarto Int. 2018, 34, 1252–1272. [Google Scholar] [CrossRef]
  28. Khosravi, K.; Pham, B.T.; Chapi, K.; Shirzadi, A.; Dou, J.; Revhaug, I.; Prakash, I.; Bui, D.T. A comparative assessment of decision trees algorithms for flash flood susceptibility modeling at Haraz watershed, northern Iran. Sci. Total Environ. 2018, 627, 744–755. [Google Scholar] [CrossRef]
  29. Shafizadeh-Moghadam, H.; Valavi, R.; Dou, J.; Chapi, K.; Shirzadi, A. Novel forecasting approaches using combination of machine learning and statistical models for flood susceptibility mapping. J. Environ. Manag. 2018, 217, 1–11. [Google Scholar] [CrossRef] [Green Version]
  30. Abdelkarim, A.; Gaber, A.; Youssef, A.; Pradhan, B.; Youssef, A. Flood hazard assessment of the urban area of Tabuk City, Kingdom of Saudi Arabia by Integrating spatial-based hydrologic and hydrodynamic modeling. Sensors 2019, 19, 1024. [Google Scholar] [CrossRef] [Green Version]
  31. Abdelkarim, A.; Gaber, A. Flood risk assessment of the Wadi Nu’man Basin, Mecca, Saudi Arabia (During the Period, 1988–2019) based on the integration of geomatics and hydraulic modeling: A case study. Water 2019, 11, 1887. [Google Scholar] [CrossRef] [Green Version]
  32. Bradbrook, K. JFLOW: A multiscale two-dimensional dynamic flood model. Water Environ. J. 2006, 20, 79–86. [Google Scholar] [CrossRef]
  33. Pourghasemi, H.R.; Pradhan, B.; Gokceoglu, C. Application of fuzzy logic and analytical hierarchy process (AHP) to landslide susceptibility mapping at Haraz watershed, Iran. Nat. Hazards 2012, 63, 965–996. [Google Scholar] [CrossRef]
  34. Mohamed, S.A. Application of geo-spatial analytical hierarchy process and multi-criteria analysis for site suitability of the desalination solar stations in Egypt. J. Afr. Earth Sci. 2020, 164, 103767. [Google Scholar] [CrossRef]
  35. Rozos, D.; Bathrellos, G.; Skillodimou, H.D. Comparison of the implementation of rock engineering system and analytic hierarchy process methods, upon landslide susceptibility mapping, using GIS: A case study from the Eastern Achaia County of Peloponnesus, Greece. Environ. Earth Sci. 2010, 63, 49–63. [Google Scholar] [CrossRef]
  36. Patrikaki, O.; Kazakis, N.; Kougias, I.; Patsialis, T.; Theodossiou, N.; Voudouris, K. Assessing flood hazard at river basin scale with an index-based approach: The case of Mouriki, Greece. Geosciences 2018, 8, 50. [Google Scholar] [CrossRef] [Green Version]
  37. Kazakis, N.; Kougias, I.; Patsialis, T. Assessment of flood hazard areas at a regional scale using an index-based approach and Analytical Hierarchy Process: Application in Rhodope–Evros region, Greece. Sci. Total Environ. 2015, 538, 555–563. [Google Scholar] [CrossRef]
  38. Somaiyeh, K.; Mehran, M. Assessment of flood hazard zonation in a mountainous area based on gis and analytical hierarchy process. Carpathian J. Earth Environ. Sci. 2017, 12, 311–322. [Google Scholar]
  39. Altuwaijri, H.A. Correction to: Morphometric network drainage analysis for railway location: Case study of Saudi Railway Company’s project. Arab. J. Geosci. 2019, 12, 1–13. [Google Scholar] [CrossRef] [Green Version]
  40. Marull, J.; Mallarach, J.M. A GIS methodology for assessing ecological connectivity: Application to the Barcelona Metropolitan Area. Landsc. Urban Plan. 2005, 71, 243–262. [Google Scholar] [CrossRef]
  41. Afifi, Z.; Chu, H.-J.; Kuo, Y.-L.; Hsu, Y.-C.; Wong, H.-K.; Ali, M.Z. Residential flood loss assessment and risk mapping from high-resolution simulation. Water 2019, 11, 751. [Google Scholar] [CrossRef] [Green Version]
  42. Gigovic, L.; Pamucar, D.; Bajić, Z.; Drobnjak, S. Application of GIS-interval rough AHP methodology for flood hazard mapping in Urban Areas. Water 2017, 9, 360. [Google Scholar] [CrossRef] [Green Version]
  43. Samanta, S.; Koloa, C.; Pal, D.K.; Palsamanta, B. Flood risk analysis in lower part of markham river based on Multi-Criteria Decision Approach (MCDA). Hydrology 2016, 3, 29. [Google Scholar] [CrossRef] [Green Version]
  44. Cao, C.; Xu, P.; Wang, Y.; Chen, J.; Zheng, L.; Niu, C. Flash flood hazard susceptibility mapping using frequency ratio and statistical index methods in Coalmine Subsidence Areas. Sustainability 2016, 8, 948. [Google Scholar] [CrossRef] [Green Version]
  45. Patra, S.; Mishra, P.; Mahapatra, S.C. Delineation of groundwater potential zone for sustainable development: A case study from Ganga Alluvial Plain covering Hooghly district of India using remote sensing, geographic information system and analytic hierarchy process. J. Clean. Prod. 2018, 172, 2485–2502. [Google Scholar] [CrossRef]
  46. Valverde, J.P.B.; Blank, C.; Roidt, M.; Schneider, L.; Stefan, C. Application of a GIS multi-criteria decision analysis for the identification of intrinsic suitable sites in costa rica for the application of Managed Aquifer Recharge (MAR) through spreading methods. Water 2016, 8, 391. [Google Scholar] [CrossRef]
  47. Hong, H.; Panahi, M.; Shirzadi, A.; Ma, T.; Liu, J.; Zhu, A.-X.; Chen, W.; Kougias, I.; Kazakis, N. Flood susceptibility assessment in Hengfeng area coupling adaptive neuro-fuzzy inference system with genetic algorithm and differential evolution. Sci. Total Environ. 2018, 621, 1124–1141. [Google Scholar] [CrossRef]
  48. Khosravi, K.; Nohani, E.; Maroufinia, E.; Pourghasemi, H.R. A GIS-based flood susceptibility assessment and its mapping in Iran: A comparison between frequency ratio and weights-of-evidence bivariate statistical models with multi-criteria decision-making technique. Nat. Hazards 2016, 83, 947–987. [Google Scholar] [CrossRef]
  49. Rimba, A.; Setiawati, M.D.; Sambah, A.B.; Miura, F. Physical flood vulnerability mapping applying geospatial techniques in Okazaki City, Aichi Prefecture, Japan. Urban Sci. 2017, 1, 7. [Google Scholar] [CrossRef] [Green Version]
  50. Nandi, A.; Mandal, A.; Wilson, M.; Smith, D. Flood hazard mapping in Jamaica using principal component analysis and logistic regression. Environ. Earth Sci. 2016, 75, 1–16. [Google Scholar] [CrossRef]
  51. Lin, K.; Chen, H.; Xu, C.-Y.; Yan, P.; Lan, T.; Liu, Z.; Dong, C. Assessment of flash flood risk based on improved analytic hierarchy process method and integrated maximum likelihood clustering algorithm. J. Hydrol. 2020, 584, 124696. [Google Scholar] [CrossRef]
  52. Sarkar, D.; Mondal, P. Flood vulnerability mapping using frequency ratio (FR) model: A case study on Kulik river basin, Indo-Bangladesh Barind region. Appl. Water Sci. 2019, 10, 17. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Change in the natural appearance of land-use due to the presence of the Al-Shamal train pathway, which impedes the natural wadi runoff; (b) detaining of Wadi Bayer floodwater in front of the Al-Shamal train pathway on 10 November 2018; and (c) flooding of Bayer village due to the detaining of floodwater in front of the Al-Shamal train pathway on 10 November 2018.
Figure 1. (a) Change in the natural appearance of land-use due to the presence of the Al-Shamal train pathway, which impedes the natural wadi runoff; (b) detaining of Wadi Bayer floodwater in front of the Al-Shamal train pathway on 10 November 2018; and (c) flooding of Bayer village due to the detaining of floodwater in front of the Al-Shamal train pathway on 10 November 2018.
Water 12 01702 g001
Figure 2. Location of the Al-Shamal train pathway relative to KSA in 2020.
Figure 2. Location of the Al-Shamal train pathway relative to KSA in 2020.
Water 12 01702 g002
Figure 3. Drainage basins affecting the Al-Shamal train pathway in Al-Qurayyat city in 2020.
Figure 3. Drainage basins affecting the Al-Shamal train pathway in Al-Qurayyat city in 2020.
Water 12 01702 g003
Figure 4. Methodology and data processing.
Figure 4. Methodology and data processing.
Water 12 01702 g004
Figure 5. The criteria used in the study: (A) slope; (B) distance from wadis; (C) hydrological soil groups; and (D) land cover affecting the Al-Shamal train pathway in Al-Qurayyat city.
Figure 5. The criteria used in the study: (A) slope; (B) distance from wadis; (C) hydrological soil groups; and (D) land cover affecting the Al-Shamal train pathway in Al-Qurayyat city.
Water 12 01702 g005
Figure 6. Criteria affecting the Al-Shamal train pathway to flood hazard in Al-Qurayyat city: (A) rainfall density; (B) flow accumulation; (C) runoff speed; and (D) drainage density.
Figure 6. Criteria affecting the Al-Shamal train pathway to flood hazard in Al-Qurayyat city: (A) rainfall density; (B) flow accumulation; (C) runoff speed; and (D) drainage density.
Water 12 01702 g006
Figure 7. Structural model of flood hazard on the Al-Shamal train pathway by integrating multicriteria analysis based on geographic information systems (GIS-MCDA) and AHP.
Figure 7. Structural model of flood hazard on the Al-Shamal train pathway by integrating multicriteria analysis based on geographic information systems (GIS-MCDA) and AHP.
Water 12 01702 g007
Figure 8. The degrees of hazard of the criteria affecting the Al-Shamal train pathway in Qurayyat city in 2020: (A) slope; (B) distance to wadis; (C) hydrological soil groups; and (D) land cover.
Figure 8. The degrees of hazard of the criteria affecting the Al-Shamal train pathway in Qurayyat city in 2020: (A) slope; (B) distance to wadis; (C) hydrological soil groups; and (D) land cover.
Water 12 01702 g008
Figure 9. (a) Field visit to the area below the Al-Shamal train pathway in Al-Qurayyat city; and (b,c) Wadi Bayer during the floods on 10 November 2018.
Figure 9. (a) Field visit to the area below the Al-Shamal train pathway in Al-Qurayyat city; and (b,c) Wadi Bayer during the floods on 10 November 2018.
Water 12 01702 g009
Figure 10. The degrees of floods hazard affecting the Al-Shamal train pathway in Qurayyat city in 2020 in the study: (A) rainfall density; (B) flow accumulation; (C) runoff speed; (D) drainage density.
Figure 10. The degrees of floods hazard affecting the Al-Shamal train pathway in Qurayyat city in 2020 in the study: (A) rainfall density; (B) flow accumulation; (C) runoff speed; (D) drainage density.
Water 12 01702 g010
Figure 11. Correlations between criteria affecting flood hazards and the severity of floods in 2020.
Figure 11. Correlations between criteria affecting flood hazards and the severity of floods in 2020.
Water 12 01702 g011
Figure 12. Degrees of flood hazard for drainage basins affecting the Al-Shamal train pathway in Al-Qurayyat city in 2020.
Figure 12. Degrees of flood hazard for drainage basins affecting the Al-Shamal train pathway in Al-Qurayyat city in 2020.
Water 12 01702 g012
Figure 13. Degrees of flood hazard for urban and agricultural areas in the study area in 2020.
Figure 13. Degrees of flood hazard for urban and agricultural areas in the study area in 2020.
Water 12 01702 g013
Figure 14. (a) Displacement of the inhabitants of the Bayer wadi as a result of flooding without the presence of floodwater drainage facilities; and (b) flooding of urban communities near the Al-Shamal train pathway.
Figure 14. (a) Displacement of the inhabitants of the Bayer wadi as a result of flooding without the presence of floodwater drainage facilities; and (b) flooding of urban communities near the Al-Shamal train pathway.
Water 12 01702 g014
Figure 15. Degrees of flood hazard on the Al-Shamal train pathway in the study area in 2020.
Figure 15. Degrees of flood hazard on the Al-Shamal train pathway in the study area in 2020.
Water 12 01702 g015
Figure 16. (a) Part of the Al-Shamal train pathway; (b) the existing structures under the train pathway which are blocked by Wadi Bayer floodwater sediments in front of the train pathway; and (c,d) the process of removing sediment from under the Al-Shamal train pathway.
Figure 16. (a) Part of the Al-Shamal train pathway; (b) the existing structures under the train pathway which are blocked by Wadi Bayer floodwater sediments in front of the train pathway; and (c,d) the process of removing sediment from under the Al-Shamal train pathway.
Water 12 01702 g016
Figure 17. Receiver operating characteristic (ROC) curve between the flood hazard values collected from the field study and the values extracted from the flood hazard map.
Figure 17. Receiver operating characteristic (ROC) curve between the flood hazard values collected from the field study and the values extracted from the flood hazard map.
Water 12 01702 g017
Table 1. Records of historical flood events that occurred in the study area during the period 1988–2018.
Table 1. Records of historical flood events that occurred in the study area during the period 1988–2018.
Floods DateThe Most Prominent Damages
1/4/1988 to 10/4/19881—Traffic halted partially.
2—Flooding of streets, houses, schools, and Al-Qurayyat hospital.
3—Partial destruction of the northern railway track and damage to its infrastructure.
4—Partial destruction of some villages and districts, such as Ghati, Bayer, and Hadithah, and the displacement of their people.
5—Flooding of whole districts (such as Al-Faisaliyah district) in 2018.
13/2/1993
4/11/1994 to 8/11/1994
10/11/1997 to 17/11/1997
14/4/2000
5/4/2001
27/10/2015
10/11/2018
Table 2. Description of the data sets used in this study.
Table 2. Description of the data sets used in this study.
CriteriaData Used to Extract CriteriaAuthority/Website Used to Download DataSpatial Resolution of Criteria
SourceScale/Spatial Resolution
Wadis networkTopographic maps1:50,000Survey Public Authority12.5 m
DEM12.5 mhttp://vertex-retired.daac.asf.alaska.edu
SlopeDEM
Land coverLandsat 830 mhttp://earthexplorer.usgs.gov
Hydrological soil groupGeological map1:250,000Saudi Geological Survey
Flow accumulationTopographic maps1:50,000Survey Public Authority
DEM12.5 mhttp://vertex-retired.daac.asf.alaska.edu
Rainfall densityRainfall rate0.25° × 0.25°http://daac.gsfc.nasa.gov
Drainage densityWadis network12.5 mhttp://vertex-retired.daac.asf.alaska.edu
Runoff speedDEM
Table 3. The classes of flood hazard for criteria affecting the Al-Shamal train pathway in Al-Qurayyat city in 2020.
Table 3. The classes of flood hazard for criteria affecting the Al-Shamal train pathway in Al-Qurayyat city in 2020.
CriteriaClassesRanks of Hazard
Flow accumulation (m3)0–45,0022
45,002–180,0084
180,008–390,0186
390,018–700,0338
700,033–1,275,06010
Distance from wadis network (km)0–0.2810
0.28–0.69
0.6–0.918
0.91–1.217
1.21–1.516
1.51–1.835
1.83–2.24
2.2–2.653
2.65–3.372
3.37–5.511
Slope (%) 0–1.71
1.7–2.782
2.78–4.073
4.07–5.574
5.57–7.495
7.49–9.636
9.63–11.997
11.99–15.28
15.2–20.349
20.34–54.610
Land coverUrban and built-up areas2
Agricultural areas4
Swamps6
Desert areas8
Wadis10
Hydrological soil groupA (deep sandy soils with very high intrusion rates)2
B (relatively fine grains soil with moderate intrusion rates)4
C (fine grains soil with low intrusion rates)6
D (cohesive soil with very fine grains and very low intrusion rates)5
Rainfall density (mm)8.96–9.891
9.89–10.482
10.48–11.023
11.02–11.494
11.49–11.955
11.95–12.446
12.44–12.967
12.96–13.598
13.59–14.39
14.3–15.2210
Drainage density (km/km2)3.8–3.832
3.83–3.864
3.86–3.886
3.88–3.918
3.91–3.9410
Runoff speed (km/h)1.12–1.132
1.13–1.164
1.16–1.176
1.17–1.218
1.21–1.2710
Table 4. Relative importance of criteria, according to Thomas Saaty.
Table 4. Relative importance of criteria, according to Thomas Saaty.
Scale/Degree of ImportanceExplanation
1Equal importance
3One of the criteria is of moderate importance with respect to the other
5One of the criteria is of high importance with respect to the other
7One of the criteria is of very high importance with respect to the other
9One of the criteria is extremely important with respect to the other
2–4–6–8Intermediate values used between the previous weights in numerical comparison
Table 5. Pairwise comparison matrix for factor criteria.
Table 5. Pairwise comparison matrix for factor criteria.
CriteriaDWFASRDDDRSLCHG
DW12234567
FA0.51123456
S0.51123456
RD0.330.50.512345
DD0.250.330.30.51234
RS0.20.250.250.330.5123
LC0.170.20.20.250.330.512
HG0.140.170.170.20.250.330.51
Total3.095.455.459.2814.0819.8326.534
DW: distance from wadis, FA: flow accumulation, S: slope, RD: rainfall density, DD: drainage density, RS: runoff speed, LC: land cover, and HG: hydrological groups.
Table 6. Percentage importance values using the analytic hierarchy process (AHP).
Table 6. Percentage importance values using the analytic hierarchy process (AHP).
CriteriaDWFASRDDDRSLCHG
DW0.3240.3670.3670.3230.2840.2520.2260.206
FA0.1620.1830.1830.2160.2130.2020.1890.176
S0.1620.1830.1830.2160.2130.2020.1890.176
RD1.070.0920.0920.1080.1420.1510.1510.147
DD0.0810.0610.0610.0540.0710.1010.1130.118
RS0.0650.0460.0460.0360.0360.050.0750.088
LC0.0550.0370.0370.0270.0230.0250.0380.059
HG0.0450.0310.0310.0220.0180.0170.0190.029
Total11111111
Table 7. Relative weight matrix, according to importance values between criteria.
Table 7. Relative weight matrix, according to importance values between criteria.
CriteriaDWFASRDDDRSLCHG
DW0.2940.3800.3800.3720.3280.2750.2280.189
FA0.1470.1900.1900.2480.2460.2200.1900.162
S0.1470.1900.1900.2480.2460.2200.1900.162
RD0.0970.0950.0950.1240.1640.1650.1520.135
DD0.0740.0630.0630.0620.0820.1100.1140108
RS0.0590.0480.0480.0410.0410.0550.0760.081
LC0.0500.0380.0380.0310.0270.0280.0380.054
HG0.0410.0320.0320.0250.0210.0180.0190.027
Total0.2940.1900.1900.1240.0820.0550.0380.027
Table 8. Weight values of each criteria using AHP.
Table 8. Weight values of each criteria using AHP.
CriteriaValue WeightRelative WeightTotal of Rows
DW2.4460.2948.320
FA1.5930.1908.384
S1.5930.1908.384
RD1.0270.1248.282
DD0.6760.0828.244
RS0.4490.0558.164
LC0.3040.0388.000
HG0.2150.0277.963
Average 8.218 = ʎMax
Table 9. The order of the matrix (N) and the equivalent random index (R).
Table 9. The order of the matrix (N) and the equivalent random index (R).
N12345678910
R000.520.891.111.251.31.41.451.49
Table 10. Weights of the eight criteria of the study using the AHP.
Table 10. Weights of the eight criteria of the study using the AHP.
CriteriaTotal of RowsRelative WeightWeight (in %)
DW2.3490.29429.4
FA1.5240.19019
S1.5240.19019
RD0.990.12412.4
DD0.660.0828.2
RS0.4420.0555.5
LU0.3010.0383.8
HG0.2120.0272.7
Total8.0021100
Table 11. Analysis of the distance from the wadi network in the study area in 2020.
Table 11. Analysis of the distance from the wadi network in the study area in 2020.
Class (km)Area (km2)Percentage (%)
0–0.282556.8623.22
0.28–0.62458.8622.33
0.6–0.911871.9517
0.91–1.211445.8113.13
1.21–1.511044.999.49
1.51–1.83720.156.54
1.83–2.2471.294.28
2.2–2.65282.992.57
2.65–3.37129.941.18
3.37–5.5128.630.26
Total11,011.47100
Table 12. Slope analysis results for Al-Shamal train pathway in Al-Qurayyat city in 2020.
Table 12. Slope analysis results for Al-Shamal train pathway in Al-Qurayyat city in 2020.
Class (Degrees)Area (km2)Percentage (%)
0–1.72164.8619.66
1.7–2.784491.5740.79
2.78–4.072131.8219.36
4.07–5.571445.8113.13
5.57–7.49562.695.11
7.49–9.63138.741.26
9.63–11.9950.650.46
11.99–15.218.720.17
15.2–20.345.510.05
20.34–54.61.110.01
Total11,011.47100
Table 13. Land cover analysis for Al-Shamal train pathway in Al-Qurayyat city in 2020.
Table 13. Land cover analysis for Al-Shamal train pathway in Al-Qurayyat city in 2020.
ClassesArea (km2)Percentage (%)
Urban areas132.141.2
Agricultural lands129.061.17
Wadis2.630.02
Swamps8950.5781.28
Desert areas1797.0716.32
Total11,011.47100
Table 14. Hydrological group analysis for Al-Shamal train pathway in Al-Qurayyat city in 2020.
Table 14. Hydrological group analysis for Al-Shamal train pathway in Al-Qurayyat city in 2020.
ClassArea (km2)Percentage (%)
A27.530.25
B10,825.3898.31
C48.450.44
D110.111
Total11,011.47100
Table 15. Flow accumulation analysis for the Al-Shamal train pathway in Al-Qurayyat city in 2020.
Table 15. Flow accumulation analysis for the Al-Shamal train pathway in Al-Qurayyat city in 2020.
Class (m3)Area (km2)Percentage (%)
0–45,00210,956.4199.50
45,002–180,00829.730.27
180,008–390,01814.310.13
390,018–700,0335.510.05
700,033–1,275,0605.510.05
Total11,011.47100
Table 16. Rainfall density analysis for Al-Shamal train pathway in Al-Qurayyat city in 2020.
Table 16. Rainfall density analysis for Al-Shamal train pathway in Al-Qurayyat city in 2020.
Class (mm)Area (km2)Percentage (%)
8.96–9.89789.527.17
9.89–10.48948.098.61
10.48–11.021435.913.04
11.02–11.491509.6713.71
11.49–11.951760.7315.99
11.95–12.442233.1320.28
12.44–12.96895.238.13
12.96–13.59770.87
13.59–14.3528.554.8
14.3–15.22139.851.27
Total11011.47100
Table 17. Drainage density analysis for Al-Shamal train pathway in Al-Qurayyat city in 2020.
Table 17. Drainage density analysis for Al-Shamal train pathway in Al-Qurayyat city in 2020.
Class (km/km2)Area (km2)Percentage (%)
3.8–3.83799.437.26
3.83–3.865983.754.34
3.86–3.8800
3.88–3.912900.726.34
3.91–3.941327.6412.06
Total11,011.47100
Table 18. Runoff speed for Al-Shamal train pathway in Al-Qurayyat city in 2020.
Table 18. Runoff speed for Al-Shamal train pathway in Al-Qurayyat city in 2020.
Class (km/h)Area (km2)Percentage (%)
1.12–1.13799.437.26
1.13–1.1600
1.16–1.171271.311.54
1.17–1.214712.442.8
1.21–1.274228.3438.4
Total11,011.47100
Table 19. Degrees of flood hazard in the study area in 2020.
Table 19. Degrees of flood hazard in the study area in 2020.
Degree of HazardPercentage of HazardFlood-Prone Areas
Area (km2)Percentage (%)
Very high68.2–91.21744.2215.84
High61.6–68.23172.428.81
Moderate55.6–61.63184.5228.92
Low48.5–55.62050.3418.62
Very low24–48.58607.81
Total11,011.47100
Table 20. Degrees of flood hazard for urban and agricultural areas exposed to flood hazards in the study area in 2020.
Table 20. Degrees of flood hazard for urban and agricultural areas exposed to flood hazards in the study area in 2020.
Degree of HazardPercentage of HazardUrban AreasAgricultural Areas
Area (km2)Area (%)Area (km2)Area (%)
Very high68.2–91.24.443.3618.3814.24
High61.6–68.224.7918.7641.4932.15
Moderate55.6–61.656.342.6137.2928.89
Low48.5–55.632.6524.7119.114.8
Very low24–48.513.9710.5712.89.92
Total132.14100129.06100
Table 21. Flood hazard degrees on the Al-Shamal train pathway in the study area in 2020.
Table 21. Flood hazard degrees on the Al-Shamal train pathway in the study area in 2020.
Degree of HazardPercentage of HazardLength (km)Length (%)
Very high68.2–91.26.969.59
High61.6–68.212.2416.86
Moderate55.6–61.66.549.01
Low48.5–55.615.2420.99
Very low24–48.531.6143.55
Total72.59100
Table 22. Area under the curve (AUC) values and the degree of validation of the results.
Table 22. Area under the curve (AUC) values and the degree of validation of the results.
The Degree of Validation of the ResultsAUC
The Results Are Unreliable and Cannot Be Relied uponLess than 0.5
Poor0.5–0.6
Fair0.6–0.7
Good0.7–0.8
Very good0.8–0.9
Excellent0.9–1

Share and Cite

MDPI and ACS Style

Abdelkarim, A.; Al-Alola, S.S.; Alogayell, H.M.; Mohamed, S.A.; Alkadi, I.I.; Ismail, I.Y. Integration of GIS-Based Multicriteria Decision Analysis and Analytic Hierarchy Process to Assess Flood Hazard on the Al-Shamal Train Pathway in Al-Qurayyat Region, Kingdom of Saudi Arabia. Water 2020, 12, 1702. https://doi.org/10.3390/w12061702

AMA Style

Abdelkarim A, Al-Alola SS, Alogayell HM, Mohamed SA, Alkadi II, Ismail IY. Integration of GIS-Based Multicriteria Decision Analysis and Analytic Hierarchy Process to Assess Flood Hazard on the Al-Shamal Train Pathway in Al-Qurayyat Region, Kingdom of Saudi Arabia. Water. 2020; 12(6):1702. https://doi.org/10.3390/w12061702

Chicago/Turabian Style

Abdelkarim, Ashraf, Seham S. Al-Alola, Haya M. Alogayell, Soha A. Mohamed, Ibtesam I. Alkadi, and Ismail Y. Ismail. 2020. "Integration of GIS-Based Multicriteria Decision Analysis and Analytic Hierarchy Process to Assess Flood Hazard on the Al-Shamal Train Pathway in Al-Qurayyat Region, Kingdom of Saudi Arabia" Water 12, no. 6: 1702. https://doi.org/10.3390/w12061702

APA Style

Abdelkarim, A., Al-Alola, S. S., Alogayell, H. M., Mohamed, S. A., Alkadi, I. I., & Ismail, I. Y. (2020). Integration of GIS-Based Multicriteria Decision Analysis and Analytic Hierarchy Process to Assess Flood Hazard on the Al-Shamal Train Pathway in Al-Qurayyat Region, Kingdom of Saudi Arabia. Water, 12(6), 1702. https://doi.org/10.3390/w12061702

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

Article Metrics

Back to TopTop