Suitable Site Selection for Rainwater Harvesting and Storage Case Study Using Dohuk Governorate

: The Middle East is an inherently dry zone. It has experienced severe drought for the last seven years, and climate change has made the situation worse. The Dohuk governorate has been su ﬀ ering from an appalling water crisis. One possible way of relieving this water crisis is by properly harvesting the rainwater. Rainwater harvesting is a widely used method of storing rainwater in the countries presenting with drought characteristics. Several pieces of research have derived and developed di ﬀ erent criteria and techniques to select suitable sites for harvesting rainwater. The main aim of this research was to identify and select suitable sites for the potential erection of dams, as well as to derive a model builder in ArcMap 10.4.1. The model combined several parameters, such as slope, runo ﬀ potential, land cover / use, stream order, soil quality, and hydrology to determine the suitability of the site for harvesting rainwater. To compute the land use / cover categories, the study depended on Landsat image data from 2018. Supervised classiﬁcation was applied using the ENVI 5 software, while the slope mapping and drainage order were extracted using a digital elevation model. Inverse distance weighting (IDW) was used for the spatial interpolation of the rain data. The results demonstrated that suitable areas for water harvesting, are located in the middle and northern part of the research area, and in intensively cultivated zones. The main soil texture in these suitable sites was loam, while the rainfall rate amounted to 750 to 900 mm. This research shows that 15% and 13% of the area studied can be categorized as having excellent and good suitability for water harvesting, respectively. Furthermore, 21% and 27% of the area studied were of moderate and poor suitability, while the remaining 24% were not suitable at all.


Introduction
Water scarcity is a major issue in several countries in the world, particularly in developing countries [1]. These adverse conditions are further exacerbated by an unsustainable usage of water, such as in overexploitation. Overexploitation causes groundwater levels and water quality to decline [2]. Nevertheless, not only groundwater and surface-water hydrology are affected, but the alteration of the soil properties has led among others to deforestation and soil erosion [3]. It is recognized that changing the natural structure of vegetation and surface coarseness affects the balance between precipitation, evaporation, and runoff responses [4].

Study Area
The study area covers the Dohuk Province in the far north-west of Iraq, situated at 36° 54' 27.72" latitude, 43° 3' 47.52" longitude, and 585 m above the sea level [29] (Figure 1). The Dohuk province was selected because of its importance, especially in terms of history and geography. The effects and sculptures discovered in the hills and covers demonstrate the significance of its geographical location, at the border of Iraq, Turkey, and Syria. The province of Dohuk is characterized by diverse terrains of high mountains and complexes. In addition, the plains are spacious and rich in agricultural resources. These are located in the south of the province. While the surrounding mountains have relatively high altitudes, the climate of the study area is related to the Mediterranean zone. In this climate, summer is characterized by a dry and hot climate with low moisture. In contrast, winter is characterized by lower temperatures, higher humidity, and reasonable rainfall as well as snowfall in the mountains. The spring seasons are characterized by cold temperatures as well as significant rainfall rates. The occasional dry season which lasts for longer periods of time leads to the undernourishment of the surface.

Data Collection
In the present study, remote sensing data (Landsat satellite imagery), digital evaluation models (DEM), soil type maps, and rainfall data were adapted to be used for runoff evaluation and water harvesting. Landsat imaging was used to detect the land use/cover maps. The Landsat 8 Operational Land Imager (OLI) data were obtained from the United States Geological Survey (USGS) website with geo-reference to UTM zone 38, WGS 84, and was taken in October 2018 with a 30 m resolution. The drainage layer, slope, and topographic maps were generated using DEM. The soil map of the study area was obtained from the Ministry of Agriculture and Water Resources of the Kurdistan Regional Government, Iraq. The soil map is defined by the texture, color, and depth of the soil in the study area. The monthly rainfall data were delivered by the Ministry of Agriculture and Water Resources of the Kurdistan Region for the years 2001 to 2017 as measured by the 15 meteorological stations in the study area.

Figure 1.
An illustration of the location of the study area. Light blue area is the Kurdistan Region and blue is study area.

Data Collection
In the present study, remote sensing data (Landsat satellite imagery), digital evaluation models (DEM), soil type maps, and rainfall data were adapted to be used for runoff evaluation and water harvesting. Landsat imaging was used to detect the land use/cover maps. The Landsat 8 Operational Land Imager (OLI) data were obtained from the United States Geological Survey (USGS) website with geo-reference to UTM zone 38, WGS 84, and was taken in October 2018 with a 30 m resolution. The drainage layer, slope, and topographic maps were generated using DEM. The soil map of the study area was obtained from the Ministry of Agriculture and Water Resources of the Kurdistan Regional Government, Iraq. The soil map is defined by the texture, color, and depth of the soil in the study area. The monthly rainfall data were delivered by the Ministry of Agriculture and Water Resources of the Kurdistan Region for the years 2001 to 2017 as measured by the 15 meteorological stations in the study area.

Methods
One multi-objective and multi-criteria faced problem is discovering locations for rainwater harvesting. Different criteria were used to detect potential areas for the construction of small dams.
According to the Food and Agriculture Organization (FAO), six criteria were selected for evaluating the suitability of the areas for water harvesting which included the runoff potential, slope (topography), land cover/use, drainage, rainwater data, soil quality, and hydrological soil groups as presented in Figures 2 and 3. ArcGIS 10.4.1 (ESRI, Redlands, CA, USA) and ENVI 5.1 (Exelis Visual Information Solutions, Boulder, CO, USA) were used to combine and analyze the data as a part of the procedure for the site of potential runoff areas and appropriate locations for runoff harvesting.  The study depends on the Landsat 8 image satellite to extract land use and land cover data with a 30 m spatial resolution. Supervised classification with a maximum likelihood logarithm was used to assess the land cover/use map. Figure 2d demonstrates the land cover classification in the study area. This research classified the study area into five main types of land use and land cover classes, namely water body, forest area, farmland, bare soil, and urban area. Table 1 summarizes the statistical details of the land cover classification: water body (0.77%), forest land (16.33%), farmland (34.49%), bare soil (46.25%), and urban area (2.15%). The water body, forest, and farmland are a more suitable class of land cover for rainwater harvesting, whereas infiltration and evapotranspiration depend on land use patterns.  Figure 3. Flowchart for generation rain water harvesting mapping, Q v is the volume of surface runoff (mm 3 ), Q is the runoff depth (mm) and A is the basin area (km 2 ).

Slope Map
The slope is a significant parameter for identifying the best site for water harvesting in order to acquire the extreme storing volume in the channel [30]. This parameter is determined by the terrain ratio, which denotes the ratio of the height difference between two points divided by the horizontal straight line distance between them [18,31]. The slope parameter has an effect on the runoff, sedimentation ratio, recharge, surface water velocity, and the number of materials essential to creating dams [25]. The DEM with a 30 m resolution was used to generate the slope map, and the area was then divided into four classes ( Figure 2c). Critchley et al. [32] indicated that the areas with slopes of 5% were not appropriate for water harvesting due to their vulnerability to high erosion degrees.

Stream Order Map
The hydrological parameters from the spatial analysis tools were derived to produce the stream order map based on the DEM data. The main surface water is watershed. In the dry season, water gathering is used for human life, livestock, and other purposes. The arrangement of the stream is based on the linking of tributaries. The stream order refers to the hierarchical linking between the flow sections and allows the drainage basins to be classified according to their size. In addition, for mapping rainwater conservation streams, the order analysis is significant due to higher stream orders having lower permeability and infiltration.

Land Cover and Land Use
One of the key criteria for selecting suitable sites for water harvesting is the land cover and land use. The hydrological response of the river basins is affected by the land cover/use change and rainfall.
Hydrologists have been concerned about the impact of land use and land cover change in the runoff dynamics of the river basin. This is a significant parameter for runoff evaluation. Due to the decrease of the initial abstraction, the runoff rate is increased. It is therefore essential to identify alterations in the runoff feature of the watershed under the influence of the land cover change. Land cover is associated with high runoff, generated by rainwater in a certain area due to higher infiltration, while low runoffs are associated with vegetation areas [33].
The study depends on the Landsat 8 image satellite to extract land use and land cover data with a 30 m spatial resolution. Supervised classification with a maximum likelihood logarithm was used to assess the land cover/use map. Figure 2d demonstrates the land cover classification in the study area. This research classified the study area into five main types of land use and land cover classes, namely water body, forest area, farmland, bare soil, and urban area. Table 1 summarizes the statistical details of the land cover classification: water body (0.77%), forest land (16.33%), farmland (34.49%), bare soil (46.25%), and urban area (2.15%). The water body, forest, and farmland are a more suitable class of land cover for rainwater harvesting, whereas infiltration and evapotranspiration depend on land use patterns.

Soil Types
The soil permeability is a significant parameter which controls the infiltration ratio and storage of water in the layers of the soil [34]. The soil textural class is determined through the percentages of sand, silt, and clay. The fine and medium textured soils are more suitable for rainwater harvesting, due to their ability to hold a higher content of water. The suitable soils for the water harvesting process have a higher capacity for water retention [35]. The soils that have a high clay content have a relatively higher runoff depth and do not permit more infiltration, whereas sand soils have lower runoff and higher infiltration. The soil map contains nine different classes in the study area such as lithosolic soils, sand and gypsum, shallow to moderately deep, grave to very grave, and lime-rich silt loam (Figure 2g). Areas in Figure 2e are simplified and represented by the numbers 32 to 39, excluding 37. Lithosolic soil in limestone soils is the first, as shown in zone 32. Brown soils and the medium shallow phase of the Bakhtiary gravel is referred to by number 33. The zone of Gypsum is 34. The deep phases of brown soil are represented in zone 35. Chestnut soils, shallow, stiny and sloping phases are represented in zone 36. As for 38 and 39, they represent rough broken stony land and rough mountainous alpine phases, respectively.

SCS-CN Model and Proposed Method
The design of the major of hydraulic structures widely utilizes the SCS-CN runoff model, which is an event-based rainfall-runoff model. The estimation of runoff depth is determined through the SCS model using the CN values as input parameters [36]. CN evaluated the study area on a per pixel base, employed land cover/use type, antecedent soil texture, and hydrological situation. The hydrological soil group's map was produced by reclassifying the soils using information indicating that infiltration rates depend on soil properties. According to the United State Department of Agriculture, soil classes were grouped under three hydrological categories-B, C, and D-depending on the soil textural classes; matching the soil layer and land use cover helped prepare the index file that was referred to by the ArcGIS runoff tool to assign the curve number to the intersected land-soil layers. Table 2 explains the soils, and their textures determine their hydrological soil groups (HSG). The spatial distributions of the HSGs are presented in Figure 2e. A, B, and C represent the three main types of hydrologic soil groups found in the study area. Soil maps of the northern governorates assisted in determining the soil type's moderate, high-to-moderate, and high potential runoff, as found in the study area and are portrayed in Table 2. HSG, land use, land treatment, and hydrological conditions are used to determine the dimensionless runoff index, CN. This method effectively presents a runoff to reflect the changes in land use. The CN is capable of reversing the influence on alterations in land use based on runoff [20,37], and its range from 1 to 100; low values from CN indicate the low runoff [2,18,31,38].

Runoff Depth Estimation
Runoff mapping is a significant principle for identifying an appropriate location for rainwater harvesting. The water source during the runoff has been calculated through the runoff depth. To evaluate the runoff, the depth Soil Conservation Service (SCS) and Curve Number (CN) were derived. Low runoff depth potentials which start from 391 mm have been located in the south and west part of the study area, while high runoff depths of up to 1014 mm have been located in the northeast part of the study area ( Figure 4).

Runoff Depth Estimation
Runoff mapping is a significant principle for identifying an appropriate location for rainwater harvesting. The water source during the runoff has been calculated through the runoff depth. To evaluate the runoff, the depth Soil Conservation Service (SCS) and Curve Number (CN) were derived. Low runoff depth potentials which start from 391 mm have been located in the south and west part of the study area, while high runoff depths of up to 1014 mm have been located in the northeast part of the study area ( Figure 4).   The Soil Conservation Service Curve Number (SCS-CN) modeling was established by the United States Department of Agriculture to be utilized in outside cities [16,40,41]. The most widely used runoff equation is the SCS-CN since it is versatile, simple, and flexible [42,43]. The SCS-CN model is a crucial method for computing the surface runoff from the rainwater [44][45][46]. The SCS-CN model was used to estimate each pixel's spatial variation runoff depth. The estimation of runoff depth used the rainfall distribution layer as the main parameter. The interpolation determined pixel values with an Inverse Distance Weighting creating a set of station points for the pixel values of the rainfall that when combined, became linear. The estimation of each runoff pixel was achieved through the application of the SCS-CN model, using the rainfall depth and its specific CN value raster layers as an input. Figure 3 shows the runoff volumes based on their land cover categories. This was achieved through intersecting the runoff rater layers with the land use/covers. The model depends on the equation of the water balance and two essential hypotheses; firstly, the equation of the actual rate of the direct surface runoff depth (Q) to the overall precipitation (P) or the maximum probability surface runoff of the actual infiltration rate (F) to the maximum possible retaining volume, and secondly, the preliminary abstraction (Ia) to the potential maximum retention (S). The main equations for the SCS-CN process are expressed below: The experiential relative was developed for the term Ia and it is provided by the empirical relationship as: where Q is the runoff depth in mm, P is the rainfall in mm, S is the potential maximum retention after the onset of runoff in mm, Ia is an initial abstraction in mm that contains all losses before the runoff commences, infiltration evaporation, and water interception by vegetation [47]. The potential maximum retention after runoff commences in mm (S) can be designed by the Curve Number (CN) as follows: By replacing Ia and S from Equations (1) in Equations (2) and (3), an equation with only two parameters was obtained:

Site Selection Criteria Using WLC
All processes for the creation of a suitable rainwater harvesting map were applied in a convenience model developed in ArcGIS 10.4.1. The convenience model creates rainwater harvesting compatibility maps by merging different criteria using a weighted linear combination process [48,49]. The Weight liner combination is a widely used multi-criteria evaluation (MCE) process for suitability site analysis [50]. This model includes the standardization of appropriateness maps, the weighting of the comparative significance of appropriateness maps, and the merging of weights and uniformity maps to achieve a suitability value [51]. The model is based on the concept of weighted averages, and it is prepared by multiplying the weight of the factor by its consistent membership score (Equation (5)): where S = suitable site, W i = weighting of factor i, and X i = the membership value of criteria i.

Suitable Sites for Water Harvesting Structures
The research applies appropriate criteria for rainwater harvesting based on a review of the literature and is mindful of soil groups, rainfall depths, slope drainage areas, and a land use and land cover map. The suitable location for rainwater harvesting structures was developed with integrates of remote sensing and GIS-based on an Analytical Hierarchy Process (AHP) analysis, taking into account five layers. Rainwater harvesting suitability sites were determined using multi-criteria evaluation (MCE). Image satellite and DEM were derived to determine the land cover/use map, the slope map, the drainage map, and the soil map. A hydrological soil group map and four land use classes were prepared as well as all the criteria for the selection of a suitable location for rainwater harvesting. Five comparable units are used as indicators for potential sites for water harvesting (according to their suitability): excellent suitability, good suitability, moderate suitability, poor, and unsuitable. Figure 5 indicates that excellent suitability sites for water harvesting were concentrated in the northern part of the study area. The outcome indicated that 15% (2204 km 2 ) and 13% (1377 km 2 ) of the study area exhibited excellent and good suitability for water harvesting, while 21% (1506 km 2 ), 27% (2842 km 2 ), and 24% (2486 km 2 ) of the site were classified as moderate, poor, and unsuitable according to Figure 6. It is clear that areas with a large excess of rainfall will have high appropriateness ranking; where this excess points to the availability of runoff to be gathered. The slope is a crucial limiting aspect to rainwater harvesting. The sites with a slighter slope were more suitable for rainwater harvesting; a small slope is required for enhanced runoff harvesting. Critchley and Siegert [32] suggested that an area with a slope of larger than 5% is unsuitable for water harvesting because of the rough runoff distribution and the great volume of earthwork required. The slope of the study area (Figure 2c) has been classified into four classes, as nearly level (0-3%), gentle slope (3-5%), moderate slope (5-8%), and strong slope (>10%). The northern part of the Dohuk governorate including Amedi, Zaxo, and Akre was identified as suitable areas for water harvesting, as it has been shown that these areas exhibit maximum rainfall rates (949 mm), as well as downstream sites of the watershed and dense hydrological networks with slopes less than 3. This led to the conclusion that the area is a suitable location for water harvesting ( Figure 6). slope (3-5%), moderate slope (5-8%), and strong slope (>10%). The northern part of the Dohuk governorate including Amedi, Zaxo, and Akre was identified as suitable areas for water harvesting, as it has been shown that these areas exhibit maximum rainfall rates (949 mm), as well as downstream sites of the watershed and dense hydrological networks with slopes less than 3. This led to the conclusion that the area is a suitable location for water harvesting (Figure 6).

Sites Proposed for Small and Medium Dams
The study derived a Triangulated Irregular Network (TIN), drainage layer, and contour line to construct the dams and harvest rainwater in the study area. The height of the dams was identified slope (3-5%), moderate slope (5-8%), and strong slope (>10%). The northern part of the Dohuk governorate including Amedi, Zaxo, and Akre was identified as suitable areas for water harvesting, as it has been shown that these areas exhibit maximum rainfall rates (949 mm), as well as downstream sites of the watershed and dense hydrological networks with slopes less than 3. This led to the conclusion that the area is a suitable location for water harvesting (Figure 6).

Sites Proposed for Small and Medium Dams
The study derived a Triangulated Irregular Network (TIN), drainage layer, and contour line to construct the dams and harvest rainwater in the study area. The height of the dams was identified

Sites Proposed for Small and Medium Dams
The study derived a Triangulated Irregular Network (TIN), drainage layer, and contour line to construct the dams and harvest rainwater in the study area. The height of the dams was identified based on the contour lines. A 5 m interval, as well as TIN, was used to calculate the height and volume of the dams with the tools of the ArcGIS 10.4.1 software (Figure 7). Sites of potential dams were then selected depending on other criteria serving to identify the suitable sites to construct the dams, as well as evaluating the runoff rate that could be stored in the dams. The preferred locations for the dams were placed where the drainage was narrow and straight to acquire the lowest cost for the constructions. The construction of the dam should be built to regulate the flow of water. based on the contour lines. A 5 m interval, as well as TIN, was used to calculate the height and volume of the dams with the tools of the ArcGIS 10.4.1 software (Figure 7). Sites of potential dams were then selected depending on other criteria serving to identify the suitable sites to construct the dams, as well as evaluating the runoff rate that could be stored in the dams. The preferred locations for the dams were placed where the drainage was narrow and straight to acquire the lowest cost for the constructions. The construction of the dam should be built to regulate the flow of water. The northern part of the study site is the most complex in topography with deep drainage networking compared to the southern part, which contains plains or has moderate slopes with depthless drainage networks. The study identified nine potential sites for new dams. Four potential dams (No.1,No.2,No.3,and No.4) were evaluated by computed features from the TIN including a profile of the dam, the storage capacity, and the height of the dams (Figure 7). The volume of the dam depends on the storing and height of the dam. In addition, to demonstrate the correlation between the height of the dams and the storage capacity for harvesting water, four sites were considered as presented in Figure 8. The volume of water storage in the dams at each level was computed by the water level that reaches the top of the dam. These outcomes illustrate a strong relationship between the height of dams and the storing capacity; Figure 7 and Table 3 Table 2). The profiles of the dams are based on the method described throughout this paper. Nevertheless, some of the proposed dam heights may not be feasible in practice because of some other factors such as the movement of material, or the dam's width [25]. In addition to the five criteria used in this study, there are other criteria not considered in this study such as geology and socio-economic aspects [18].  The northern part of the study site is the most complex in topography with deep drainage networking compared to the southern part, which contains plains or has moderate slopes with depthless drainage networks. The study identified nine potential sites for new dams. Four potential dams (No. 1,No.2,No.3,and No.4) were evaluated by computed features from the TIN including a profile of the dam, the storage capacity, and the height of the dams (Figure 7). The volume of the dam depends on the storing and height of the dam. In addition, to demonstrate the correlation between the height of the dams and the storage capacity for harvesting water, four sites were considered as presented in Figure 8. The volume of water storage in the dams at each level was computed by the water level that reaches the top of the dam. These outcomes illustrate a strong relationship between the height of dams and the storing capacity; Figure 7 and Table 3 Figure 8. Location of the dams in the study area, Elevation is shown in meter.

Conclusions
The Dohuk governorate is challenged by problems related to scarce water sources. Some of the reasons for the water crisis are low rainfall, high runoff, and evaporation. In this paper, remote sensing data were applied to monitor changes in land use/cover in the Dohuk governorate as well as the derived GIS-based suitability model to determine potential sites for rainwater conservation and the erection of nine dams in suitable areas. The study attempted to select sites for harvesting rainwater. MCE was used to assess the plans. MCE was obtained by using a weighted linear combination of a dependable method built on an appropriate map for water harvesting. The study considered various parameters such as slope, runoff potential, land use/cover, soil types, and stream order. The applied techniques used in this paper were very effective in attaining the subject of this work. The outcome of the identified site selection showed that the northern part of the study area included mostly hill and mountain areas with deep drainage networking compared to the southern part, which contained plains or had moderate slopes with depthless drainage networks. The results demonstrated that excellent suitable sites cover 15% of the area, while the moderate suitable sites covered 21%. Meanwhile, the rest of the zones, 51%, spread in the southern part of the study area had poor suitability or were entirely unsuitable for rainwater harvesting. In addition, the paper applied the model of soil preservation for evaluating the runoff depth, using the average rainfall measurements from the years 2000 to 2017, which had been registered in Dohuk. It results that the southern parts of the study area have minimum runoff while the northern parts have maximum runoff rates.