Mapping of Flood-Prone Areas Utilizing GIS Techniques and Remote Sensing: A Case Study of Duhok, Kurdistan Region of Iraq

: One of the most common types of natural disaster, ﬂoods can happen anywhere on Earth, except in the polar regions. The severity of the damage caused by ﬂooding can be reduced by putting proper management and protocols into place. Using remote sensing and a geospatial methodology, this study attempts to identify ﬂood-vulnerable areas of the central district of Duhok, Iraq. The analytical hierarchy process (AHP) technique was used to give relative weights to 12 contributing parameters, including elevation, slope, distance from the river, rainfall, land use land cover, soil, lithology, topographic roughness index, topographic wetness index, aspect, the sediment transport index, and the stream power index in order to calculate the Flood Hazard Index (FHI). The relative importance of each criterion was revealed by a sensitivity analysis of the parameter values. This research developed a ﬁnal ﬂood susceptibility map and identiﬁed high-susceptible zones. This was classiﬁed anywhere from very low to very high classiﬁcations for its potential ﬂood hazard. The generated map indicates that 44.72 km 2 of the total land area of the study area in Duhok city has a very high susceptibility to ﬂooding, and that these areas require signiﬁcant attention from government authorities in order to reduce ﬂood vulnerability.


Introduction
Flooding is a regular natural disaster that causes significant damage to people and property [1].Droughts and floods have historically been the most deadly natural disasters [2].For instance, drought has been a challenge on both the Indian Peninsula and in California [3,4], and floods have also affected countries in southern and southeast Asia, especially Pakistan and Thailand.According to Kron [5], floods are produced by severe rainfall and the melting of snow, which causes rivers to overflow and become temporarily stagnant.This in turn causes damage to the communities that are located along the river.Flooding can occur as a result of a river overflowing, as well as other factors like abrupt cloudbursts in semi-arid locations, resulting in flash floods [6].The incidence of floods can be considerably influenced not only by climatic factors, but also by human interventions such as reservoir storage and operations [7].Because flooding is such a complex natural event, it has long attracted researchers from all over the world who seek to better understand and analyze the factors that can lead to improved management and avoidance of Duhok city, which is situated in the northern region of Iraq.Mustafa, et al. [37] investigated the role that extreme rainfall played in the occurrence of urban flash floods in Erbil.The purpose of this work is to make use of remote sensing and geospatial analysis in order to investigate a variety of contributing factors that are accountable for the occurrence of floods in the city of Duhok in Iraq.Another aim of this study was to identify flood-prone areas in order to develop solutions for the assessment and management of floods, as well as to establish which area merits the most involvement in the development of risk reduction measures or mitigation procedures.Particularly in Duhok, no one has ever attempted a study for the purpose of conducting an analysis of flood susceptibility by comparing historic flood areas.

Study Area
The city of Duhok is located in the Kurdistan Region of Iraq.In terms of population, it is one of the least populous cities in the region.Duhok is the region's economic hub.The study area covers an area of 180.9 km 2 and is located between 36 • 46 N and 36 • 56 N latitudes and 42 • 48 E and 43 • 08 E longitudes (Figure 1).The district of Duhok is situated in a low-lying terrain that is sandwiched between the mountain ranges of Bekhair to the north and northeast and Zawa to the southeast [38].The climate in Duhok is influenced by that of the Mediterranean; therefore, the summers are typically dry and hot, while the winters are cold and wet [39].There are four distinct seasons in Duhok, each with its own climate, which range from generally cold in the winter (December-February), to hot in the summer (June-August), to moderate in the fall (September-November) and spring (March-May) [40].The average annual rainfall is 700 mm and the temperature varies between 21 to 29 • C. The elevation of Duhok district varies from 334 to 1212 m above sea level.The research area is frequently flooded near the river due to strong orographic rainfall during the spring and winter.Many tiny communities are located in such areas, and they are severely impacted both physically and economically.

Data sources
For the purpose of flood susceptibility mapping, open-source spatial data for the research area, as well as secondary information from a variety of sources including local authorities, were collected (Table 1).

Data Type
Description Source Downloaded from mostwiedzy.pl

Data Sources
For the purpose of flood susceptibility mapping, open-source spatial data for the research area, as well as secondary information from a variety of sources including local authorities, were collected (Table 1).

Methods for Preparing Flood Thematic Layers
The creation of the vulnerable spatial flood model, which is highly intricate and detailed, necessitates a large amount of topographical and hydrological data.The flood susceptibility maps are verified with the help of measurements that have been established through scientific research.As a result, one of the most important tasks is to identify the variables that cause flooding.The twelve flood-inducing parameters in the current study area were chosen based on the available flood susceptibility literature, including elevation, slope, distance from the river, rainfall, Land Use Land Cover (LULC), soil, lithology, Topographic Roughness Index (TRI), Topographic Wetness Index (TWI), aspect, the Sediment Transport Index (STI), and the Stream Power Index (SPI) [42][43][44].All of the variables that had an effect on the outcome were put through a rasterization process, and the resulting raster had a spatial resolution of 30 m.
For flood modeling, topographic parameters are critical because they influence the hydrological properties of the studied region both directly and indirectly [45,46].In the first phase, a Digital Elevation Model (DEM) for the research area was used in the ArcGIS 10.7 environment using the SRTM with a 30 m spatial resolution.We used DEM to compute slope, aspect, curvature, TWI, STI, SPI, and TRI.

Elevation
When it comes to flood modeling, elevation is the most important factor [27].Flooding is inversely proportional to height.Flooding is less likely to occur at higher elevations, and vice versa [47].As a direct consequence of this, floods are more likely to occur in locations that are either completely flat or have lower elevations.The elevation map of the research area is displayed in Figure 2a below.

Aspect
Another factor that influences the direction of flooded areas is aspect.Water flows, which helps to keep the soil moist [48].As a result, the aspect has an indirect impact on floods.For example, the section of a shady slope where soil humidity is quite high and with a lot of runoff.Figure 2b shows the aspect of the study area.

Downloaded from mostwiedzy.pl
The SPI has a significant influence on the fluvial system.The SPI is calculated using the equation below.

=
(2) The exact catchment area is represented by and the slope gradient is denoted by (radians) [52].The stream power index map is shown in Figure 2f.where each pixel of the slope is denoted by , and the upstream region is denoted by .The STI value ranges from 0 to >1 (Figure 3a).

Slope
Another important factor that has a significant impact on a flood is the slope, which determines how quickly the water will move when it floods [49].When the rate of infiltration is low, there is less risk that water will become stagnant.On the other hand, when the flow velocity is great, the higher the slope angle.As a result, lower and flat areas are more prone to floods.Figure 2c illustrates the slope of the study area.

Downloaded from mostwiedzy.pl
Remote Sens. 2023, 15, 1102 6 of 20 2.3.4.Topographic Roughness Index (TRI) One of the key contributing variables for flood episodes is TRI.It is determined by the study basin's local topography.The lower the TRI number, the more likely it is that a flood will occur.The TRI map was preserved in stretch format in this investigation, with values ranging from 0.11-0.88(Figure 2d).

Topographic Wetness Index (TWI)
There is a strong relationship between TWI and the occurrence of floods [50], due to the fact that it makes it possible to spatially visualize the differences in a basin's level of wetness [51].This index depicts the quantity of water present in each pixel of the area, and it is calculated using the following equation.
A s and β denote the catchment area (m 2 m −1 ) and slope gradient (in degrees) for the particular catchment area (m 2 m −1 ), respectively.Floods and high TWI levels are inextricably related.In the study area, the TWI value ranges from 2.86 to 22.22 (Figure 2e).

Stream Power Index (SPI)
The SPI has a significant influence on the fluvial system.The SPI is calculated using the equation below.
The exact catchment area is represented by A s and the slope gradient is denoted by β (radians) [52].The stream power index map is shown in Figure 2f.

Sediment Transport Index (STI)
Another component that causes flooding is the STI, which can upturn the frequency of floods, resulting in foundation damage.The STI is derived from DEM using the following equation [53].

STI =
A s 22.13 0.6 Sinβ 0.0896 1.3 (3) where each pixel of the slope is denoted by β, and the upstream region is denoted by A s .The STI value ranges from 0 to >1 (Figure 3a).

Land Use Land Cover (LULC)
Land use land cover has an effect on surface runoff as well as the transport of sediment, which in turn has a direct effect on the number of floods that occur [54] because land use has a direct impact on surface runoff generation and infiltration.Because metropolitan areas do not allow water to enter and generate surface water, flooding is more prevalent.The LULC map was produced using remote sensing data.When cloud cover was less than 5%, a satellite image from Sentinel-2A (2017-present) was downloaded from the European Space Agency's Copernicus open access hub platform (https://scihub.copernicus.eu/)(accessed on 22 September 2022).The image was captured by a satellite on 18 July 2022.In this work, we used the Arc GIS image classification approach to create a land use map from the Sentinel-2 satellite imageries.Extensive preprocessing was performed, including georeferencing, layer stacking, extraction, and image improvement.The land use map was divided into four categories: agriculture land, forest and vegetation, built-up area, and bare land.Using ArcMap tools, supervised classification was performed.The user of the software creates spectral signatures of known classes, such as urban, water, and agricultural, and then the software assigns each pixel to the cover type to which its signature is most similar (Figure 3b).Downloaded from mostwiedzy.pl orded maximum average rainfall (747 mm) in the same period.The average annual rainfall in the area of study is 588 mm.For this study, the average yearly rainfall values were derived from the average monthly rainfall.To compute the spatial distribution of average rainfall for the whole study area, the Inverse Distance Weighted (IDW) method was used within the ArcGIS environment.The raster layer of rainfall spatial distribution was created at 30-m resolution.The amount of rainfall in the study area increased from southwest to northeast, as shown in the figure (Figure 3f).

Analytical Hierarchy Process (AHP) Model
After completing all of the thematic layers, the analytical hierarchy process model was utilized to assign different weights to each of the parameters investigated.The analytical hierarchy method is a versatile and well-structured technique for interpreting complex decision scenarios involving several factors [55].AHP is a well-established approach for solving complex decision-making problems.The advantage of using this method lies in the allotment of user-specified weights to a hierarchical structure of criteria and sub criteria.Therefore, AHP assists the user to evaluate the different components that are more important according to the requirements, and comprehend the problem with deeper understanding.In various problems such as flood susceptible mapping, landslide mapping, and groundwater potential zones assessment, AHP has been found to be a powerful and a beneficial tool with high accuracy [29,[56][57][58][59][60][61].

Distance from the River
The majority of communities that are impacted by floods are typically situated close to the river's edge.Distance from the river is a critical conditioning feature for defining the areas of a basin that are prone to flooding, because it governs the incidence of flooding and the river flow to river factor.The more away from the river a location is, the less probable it is that it will be flooded.In this work, we determined the distance to the river by delineation of the watershed in the studied area, using the hydrology toolbox in ArcMap.The distance from the river (m) chart is depicted in Figure 3c.

Soil
In the process of rainfall and runoff, the different types of soil are one of the most important influencing factors.However, the qualities of the soil directly limit water infiltration, which in turn determines the generation of rainfall-runoff.Other factors, such as the local meteorological conditions and erosion, also play a role in regulating rainfall-runoff production.In the current investigation, the soil map of the studied area was extracted using a basin boundary from the prepared digital soil map of Iraq by the Food and Agriculture Organization (FAO).Figure 3d shows the soil map of the study area.

Lithology
The lithological properties of a region influence porosity and permeability.Such factors/properties of rock can impact flood intensity.In the present study, the geology map was prepared using Geological map of Iraq data (1:1,000,000 scale).Most of the study area is covered by Clastic rock, followed by marl limestone and gypsum, carbonates, etc. [41] (Figure 3e).

Rainfall
It has been determined that rainfall is one of the most important factors in influencing the occurrence of floods, as a large amount of rain in a short period of time might produce floods.The Ministry of Agriculture and Water Resources of the Kurdistan Regional Government was the source of the collected rainfall data.The rainfall data of the study area were extracted from five meteorological stations; three stations are within the study area and two of them are around the study area (Semel, Malta Nasara, Duhok, Qarqarava, and Zawita stations).Spatial points distribution of rainfall stations (monthly rainfall totals are available for every year) were used to interpolate the whole rainfall of the study area.Semel station recorded minimum average rainfall (430 mm), while Zawita station recorded maximum average rainfall (747 mm) in the same period.The average annual rainfall in the area of study is 588 mm.For this study, the average yearly rainfall values were derived from the average monthly rainfall.To compute the spatial distribution of average rainfall for the whole study area, the Inverse Distance Weighted (IDW) method was used within the ArcGIS environment.The raster layer of rainfall spatial distribution was created at 30-m resolution.The amount of rainfall in the study area increased from southwest to northeast, as shown in the figure (Figure 3f).

Analytical Hierarchy Process (AHP) Model
After completing all of the thematic layers, the analytical hierarchy process model was utilized to assign different weights to each of the parameters investigated.The analytical hierarchy method is a versatile and well-structured technique for interpreting complex decision scenarios involving several factors [55].AHP is a well-established approach for solving complex decision-making problems.The advantage of using this method lies in the allotment of user-specified weights to a hierarchical structure of criteria and sub criteria.Therefore, AHP assists the user to evaluate the different components that are more important according to the requirements, and comprehend the problem with deeper understanding.In various problems such as flood susceptible mapping, landslide mapping, and groundwater potential zones assessment, AHP has been found to be a powerful and a beneficial tool with high accuracy [29,[56][57][58][59][60][61].
The AHP methodology is divided into three steps.The first step of AHP is the preparation of a pairwise comparison matrix of all of the parameters affecting a system (Table 2).This judgement is purely based on the knowledge or the expertise of the user.The second step deals with the evaluation of different factors or parameters.The third and final step is the evaluation of the consistency of the AHP model to ascertain whether all the concerned parameters are consistent with each other or not.A well-established workflow for AHP model is shown in Figure 4. Downloaded from mostwiedzy.plThe AHP model's widespread appeal can be attributed to the fact that it is able to manage both quantitative and qualitative data in the process of evaluating natural disasters, such as floods.Generally, for flood-related studies, research scholars generally utilizes eight to nine factors [17,62,63]; however, the present study considered as many as twelve flood conditioning factors.For investigation problems regarding flood modeling or susceptibility studies, the AHP model is used to estimate the weights and rankings of various factors.This paper's analytical hierarchy approach is divided into two sections.The first component is the primary categorization system for all of the parameters, in which values are assigned to each parameter based on its relevance, and weights are calculated.The second portion is built by breaking down all of the criteria into subcategories (Table 3).In this case study, each of the twelve parameters utilized for hazard area mapping are individually evaluated.In addition, we divided each factor into a number of classes and gave a The AHP model's widespread appeal can be attributed to the fact that it is able to manage both quantitative and qualitative data in the process of evaluating natural disasters, such as floods.Generally, for flood-related studies, research scholars generally utilizes eight to nine factors [17,62,63]; however, the present study considered as many as twelve flood conditioning factors.For investigation problems regarding flood modeling or susceptibility studies, the AHP model is used to estimate the weights and rankings of various factors.This paper's analytical hierarchy approach is divided into two sections.The first component is the primary categorization system for all of the parameters, in which values are assigned to each parameter based on its relevance, and weights are calculated.The second portion is built by breaking down all of the criteria into subcategories (Table 3).In this case study, each of the twelve parameters utilized for hazard area mapping are individually evaluated.In addition, we divided each factor into a number of classes and gave a weight to every class.The maximum and minimum weights for every parameter class ranged between 1 and 5.The class weights utilized were proportional to the impact Downloaded from mostwiedzy.pl of each class on the hazard of flooding.All factors used the same classification, although their contributions to amplitudes varied [64].The normalized weightage was multiplied with the rank of all individual parameter sub classes to obtain the overall weightage of all thematic layers.Finally, all of the weightage values were attributed in an Arc GIS file to obtain the final flood susceptibility map.The consistency of the following method is carried out by the following equation [55]: where, CR denotes the consistency ratio, CI is the consistency index, RI denotes the random index, λ max denotes the matrix's principal eigenvalue, and n denotes the number of components or factors.In Saaty [55], a significant number of components were used to obtain RI values.It has been discovered that RI is influenced by a variety of circumstances.According to Saaty [55], the comparable value of RI was 1.45 when there were nine components.The consistency ratio (CR) is acceptable in this study if it is less than 0.10, and suggests a fair level for the pairwise comparison matrix.Eventually, the calculated consistency ratio was 0.091677, which was lower than the threshold 0.10; the weights' consistency was affirmed in the present study.

Flood Hazard Index (FHI)
The model that was implemented is capable of conducting a multi-criteria analysis, and incorporates the FHI.The FHI aims to be of assistance in the localization of high-risk areas for flooding.In the beginning, the GIS is used to process the information that has been collected from a variety of data sources.To proceed with the evaluation of the flooding probability rate, the FHI is computed using the following equation [65]: where W i is the weight of each parameter, R i is the rating of the factors, and n corresponds to the number of parameters.Therefore, Equation ( 2) is written in the next form: In order to assess the effect of each parameter on the vulnerability index, a final single parameter sensitivity analysis is performed.This test of sensitivity is carried out as a result of the arbitrary nature of the numerical values that are assigned to these parameters [66].The arbitrary weights that were used in the beginning by AHP have been replaced with effective weights [64,65], and the results are derived using Equation (8).

Downloaded from mostwiedzy.pl
where P r and P w stand for the ratings and weights of the parameter P, respectively, and V is the vulnerability index that was generated using Equation (7).After that, the effective weights, which can be found in Table 4, are utilized in the process of calculating the revised Flood Hazard Index that comes from the sensitivity analysis (FHIS).The FHI index and the FHIS index both conduct research on the same factors and class rating, but the FHIS index uses a different weighting system (the average effective weight of the sensitivity analysis is transformed in the scale of 5).

Results
In the GIS framework, the twelve maps that were created after the AHP and weighting approach were integrated and superposed with each other.The AHP weightage for all the thematic layers, namely elevation, slope, distance from the river, rainfall, LULC, soil, lithology, TRI, TWI, aspect, STI, and SPI, were 0.207, 0.174, 0.174, 0.134, 0.085, 0.085, 0.0378, 0.037, 0.0229, 0.0174, 0.0116, and 0.0116 respectively.These weightage values were multiplied with all parameter factors and the FHI was prepared.
After constructing the final map of the flood-prone area, we divided it into five primary classes, ranging from very low to very high in terms of the vulnerability of flooding occurring there.The final FHIS classes ranged from 0.951 to 38.06, and they were classified into five classes, namely very low, low, moderate, high, and very high.These classes occupied the areas of 2.10, 27.72, 38.97, 69.40, and 44.72, respectively.The majority of the area was covered by high and very high flood susceptibility zones (Figure 5).Clearly, the majority of investigated areas are susceptible to flooding.According to the categorization, out of the total area of 180 km 2 , approximately 114.12 km 2 was classified as having a high or a very high flooding hazard.This represents approximately 63.4% of the investigated total area.The high and very high flood susceptibility zones were derived from low-lying areas (<1000 m), flat areas with a slope less than 10 • , a distance from the river less than 500 m, the highest stream power and sediment transport index region, and finally, hard rock terrain and a low permeability soil type region.Furthermore, areas in the high and very high flood susceptibility zones were located in built-up areas.This tells us that the lives and properties of residents are at danger in any heavy rainfall event since built-up areas are more likely to experience flooding.Urban surface flooding is defined by the Federal Emergency Management Agency (FEMA) of the United States as the inundation of assets or any property in a developed environment, particularly in areas where the population is denser and an impervious surface dominates with poor drainage or sewer systems.Urban surface flooding can also include water that enters houses through sewer pipes [67].The total area of the zones that had a very low possibility of flooding was only 2.10 km 2 , which was equivalent to the 1.20% area under consideration.Very low flood hazard areas are generally found in areas that are both extremely elevated (at the top of mountains) and extremely steep and rocky.The terrain between mountains and high hills, as well as hills Downloaded from mostwiedzy.pl in the northern and southwest parts of the research area, contained the majority of the moderately susceptible areas.Only 20% of the total research area was included in the zone of moderate susceptibility, which was 38.97 km 2 in area.
assets or any property in a developed environment, particularly in areas where the population is denser and an impervious surface dominates with poor drainage or sewer systems.Urban surface flooding can also include water that enters houses through sewer pipes [67].The total area of the zones that had a very low possibility of flooding was only 2.10 km 2 , which was equivalent to the 1.20% area under consideration.Very low flood hazard areas are generally found in areas that are both extremely elevated (at the top of mountains) and extremely steep and rocky.The terrain between mountains and high hills, as well as hills in the northern and southwest parts of the research area, contained the majority of the moderately susceptible areas.Only 20% of the total research area was included in the zone of moderate susceptibility, which was 38.97 km 2 in area.The very high and high flood susceptibility zones in the research area are generally distributed throughout the geography of the area and are not limited to a specific locality.This could result in more severe flooding in places that are already at a moderate level.It is more probable for flash floods to occur in regions that have a dry climate and impervious terrain.This is because a lack of pervious surface or vegetation enables heavy rains to flow overland rather than seeping into the earth, which increases the risk of flooding.The results further display that most of Duhok city or the nearby area comes under very high flood-susceptible zones, and substantial damage could occur in the future in these areas.In addition, the majority of the neighborhoods or streets that are located close to the Hishakru and Duhok rivers have a high to very high flood susceptibility.This demonstrates the presence of a river that transports a significant volume of water during the wet season and occasionally floods areas that are located nearby; this will be discussed more in the validation section.Flood susceptibility maps make it abundantly clear that certain areas in the northwest and southwest directions contain high and very high flood-prone areas; however, due to the fact that these areas primarily consist of agricultural and bare land, it is expected that the impact of the flood will be less severe in these areas compared to those that contain built-up areas.As a result of the heavy rains and flooding that occurred in March 2020 in the province of Duhok, more than 400 houses, 61 vehicles, and hundreds of shops sustained damage (Figure 6) [68].These locations can be found inside the map of areas that are highly susceptible to flooding.Post-rain flooding is a recurrent problem across the country, largely due to chronically inadequate infrastructure, such as ineffective or poorly maintained urban rainwater management and sewage systems.

Downloaded from mostwiedzy.pl
to those that contain built-up areas.As a result of the heavy rains and flooding that oc curred in March 2020 in the province of Duhok, more than 400 houses, 61 vehicles, and hundreds of shops sustained damage (Figure 6) [68].These locations can be found inside the map of areas that are highly susceptible to flooding.Post-rain flooding is a recurren problem across the country, largely due to chronically inadequate infrastructure, such a ineffective or poorly maintained urban rainwater management and sewage systems.

Sensitivity Analysis and Validation
Sensitivity analysis and validation are two phases that are crucial aspects of flood susceptibility modeling.Overall, the reliability of the suggested approach for the assess ment of the flood hazardous areas was greatly improved as a result of the validation o the weights of the FHI index.Table 4 shows the effective weight that was determined by the sensitivity analysis.

Sensitivity Analysis and Validation
Sensitivity analysis and validation are two phases that are crucial aspects of flood susceptibility modeling.Overall, the reliability of the suggested approach for the assessment of the flood hazardous areas was greatly improved as a result of the validation of the weights of the FHI index.Table 4 shows the effective weight that was determined by the sensitivity analysis.
To validate the flood susceptibility mapping, a success rate was required.Using 41 training flood areas, the success rate was computed.The historic flood areas were compiled using information from news reports, reports from the directorate of water and sewage in Duhok, and before and after satellite pictures (https://earthexplorer.usgs.gov/)(accessed on 25 October 2022).There are several areas in Duhok city that are prone to flooding, as can be observed, and this is because of the previously described variables that lead to the generation of surface floods.These causes include terrain, poor drainage systems, and changes in land cover.As a result of this, it has been determined that the neighborhoods of Baroshke, Mahabad, Diyari, Serhildan, Zirka, KRO, Bntika, and Malta are the most vulnerable places inside the city of Duhok.The upcoming flood event might occur in high and extremely high vulnerability zones [69,70].In addition, basic statistical methods were applied in order to check and confirm the results of the susceptibility maps.These maps demonstrate the degree to which the model accurately predicts the extent of flooding compared to historical flood locations.The validation of flood susceptibility mapping shows thirty-six, four, and one flood location in the area of very high, high, and moderate flood susceptibility zones (Figure 7).According to the findings of the validation, 97.5% of the places that have been flooded in the past are areas that are classified as having a high or very high vulnerability in this study.Downloaded from mostwiedzy.pl

Discussion
The methodology that has been presented for the estimate of flood hazard zones has the potential to be an effective strategy for the mitigation of the devastation that is caused by floods.In addition, the validation method that was used, which also considered floods that occurred in the past, led to the computation of the modified FHIS index, which is able to provide support for the study.Within the scope of this investigation, the FHIS index has shed light on the fact that distance from rivers and elevation are both significant factors in determining the extent of flooding, indicating that the consideration of these factors is essential for developing effective flood prevention strategies.
Particularly, the modified index (FHIS) reveals that lowland areas close to rivers are even more susceptible to flooding.The sensitivity analysis specifically considers the vicinity of rivers in the category of very highly prone locations.The FHIS findings are supported by historical flood records, which also show that settlements near rivers are frequently flooded.Additionally, Figure 7 shows how FHIS analysis categorizes locations with a lot of recorded flood episodes as highly susceptible areas, another sign of its accuracy.
According to the effective weights that were utilized in the sensitivity analysis (Table 4), the FHI assumption inaccurately underestimated elevation, distance from the river, LULC, and soil.Simultaneously, the FHI assumption overestimated slope, rainfall, and TWI.Initially, FHI considered that elevation was the most important factor to consider.However, the results of the sensitivity analysis showed that the factors with the greatest impact on the investigated area were distance from the river, followed by elevation, soil, and slope.This perspective was supported as well in [71].
Obviously, there is a correlation between the zones with the highest vulnerability to flooding and urban areas.Consequently, this index gained a greater weight.However, this effect is particularly crucial to flooding and has a substantial impact on the resulting flood susceptibility map.It is also important to note that four factors, distance from river, elevation, soil type, and slope, were given considerable weights, as shown.The comparison of the FHI and FHIS indices highlighted the influence and weight of each parameter in flood hazard assessment.

Conclusions
The AHP and geospatial techniques were utilized to demarcate flood-prone regions in Iraq's Duhok district, resulting in a time-consuming and cost-effective methodology.Elevation, slope, distance from the river, topographic wetness index, rainfall, the topographic roughness index, the stream power index, the sediment transport index, land use, geology, aspect, and soil map were all included in the ArcGIS platform.As a consequence, a flood susceptibility map was developed, revealing that about 48.16 km 2 of the research region was extremely vulnerable to floods.These regions are mainly found in the middle section of the Duhok district, where the slope grade rapidly changes and numerous minor streams join the main river.The validation of the data further reveals that 97.5% of the historic flood locations fall into high and very high flood susceptible zones.The approach utilized in this study has the potential to be used to provide flood control recommendations in the study area.The availability of temporal datasets for these factors will determine future temporal analysis.Furthermore, the current study's future scope includes vulnerability and damage assessments of various subdivisions based on demographic criteria by assessing the vulnerability index.These models may be enhanced further by conducting sensitivity analysis on many key factors.Flood prevention techniques have recently received more attention, particularly in Duhok, where flash floods occur on an annual basis.Other basins and sites, on the other hand, have yet to undergo flood control testing.As a consequence, the current study was conducted in the Duhok district, employing AHP to obtain relevant information on ways that may be utilized to aid local governments and other parties in the implementation of successful flash flood mitigation measures and land-use policy planning.Downloaded from mostwiedzy.pl

Figure 2 .
Figure 2. Thematic layers: (a) Elevation, (b) Aspect, (c) Slope, (d) Topographic Roughness Index, (e) Topographic Wetness Index, (f) Stream Power Index.2.3.7.Sediment Transport Index (STI)Another component that causes flooding is the STI, which can upturn the frequency of floods, resulting in foundation damage.The STI is derived from DEM using the following equation[53].

Figure 4 .
Figure 4. Workflow of AHP methodology to prepare flood susceptibility map for part of Duhok district.

Figure 4 .
Figure 4. Workflow of AHP methodology to prepare flood susceptibility map for part of Duhok district.

Table 2 .
Each parameter's comparison matrix and relative score.

Table 3 .
Sub criteria of each parameter and their weights.