Minimizing the Impacts of Desertification in an Arid Region: A Case Study of the West Desert of Iraq

University of Anbar, College of Engineering, Department of Dams and Water Resources Engineering, University Campus, Ramadi, Anbar, Iraq Faculty of Engineering and Built Environment, Universiti Kebangsaan Malaysia, 43600 UKM Bangi, Selangor Darul Ehsan, Malaysia Faculty of Civil Engineering, Universiti Teknologi MARA, 40450 Shah Alam, Selangor Darul Ehsan, Malaysia Civil, Environmental and Natural Resources Engineering, Lulea University of Technology, 97187 Lulea, Sweden


Introduction
With the rising expectations and needs of local citizens, developing countries are turning to greater resource development to create sustainable projects. However, desertification is a soil degradation issue that represents tremendous environmental challenges and negative impacts, especially in term of future agricultural and socioeconomic development in semiarid and arid regions, many scholars [1][2][3] define the desertification as the land degradation resulting from different factors, containing climatic difference and human activities in semiarid, arid, and dry subhumid areas. Climatic and physical conditions have produced widespread natural losses and cumulative soil erosion by water and wind which result in a decrease of the biological productivity of the land [4,5]. Desertification seems once land degradation develops irreparable or when forfeiture of total productivity extents 50% to 66% [6]. Currently, desertification is a major problem in the west desert of Iraq and leads to population displacement. A total of 4,263 families were displaced due to desertification between December 2007 and June 2009 [7]. Additionally, desertification has caused the deterioration of plants and soil covers that has harmfully affected nearly 70% of the dry land. e total area that is exposed to desertification in Iraq is estimated to be 167,000 km 2 , which is amounting to 40% of its entire area of 437,500 km 2 [8,9]. erefore, the desertification phenomenon, especially in the west desert of Iraq in the near future is irrefutable, unless effective methods and research are conducted to develop a database regarding this environmental threat. e severe nature, remoteness, and size of the desert make it difficult and expensive to understand this environmental threat. Parallel to this effort, mitigation plans can also be implemented to minimize the risk of desertification. However, this process would be limited or hindered by a lack of information.
erefore, an effective solution to mitigate the problem of desertification is needed.
One of the proposed methods for reducing the risk and impact of desertification in this study is rainwater harvesting (RWH). RWH is a significant component in any water resource system for ensuring a continuous supply, in order to meet the demand for domestic uses, livestock, agriculture, and small industries [10][11][12][13]. Moreover, the proposed rainwater harvesting system should be able to maintain a sustainable water resource and minimize the risk of further degradation of resources in the area. Four procedures are taken into consideration for sustainable RWH consisting of the selection of appropriate criteria, rating of the appropriateness of the ranking for every criterion, choosing of the locations and generation of appropriateness maps intended for the nominated sites for RWH. e process of sustainable RWH may be time-consuming and difficult when the whole of these procedures was considered particularly in large watersheds with restricted data. Geographic information systems (GISs) and geospatial data techniques are integrated to find baseline data related to droughts, desertification, and degradation, and to ease the process of attaining sustainable RWH which may directly affect agricultural development and activities in arid regions [14][15][16]. e combination of geospatial data with GIS approaches has guided much attention towards identifying suitable sites for RWH. e adaption of geospatial data together with GIS techniques has made it effortless to produce database concerning the hydrological potential of an area [17,18]. Several methods and criteria have been exercised to identify appropriate sites for RWH. Many researchers [19][20][21][22][23] have applied geospatial data with GIS for identifying suitable sites for RWH structures. Since water deficiency is the major problem in semiarid and arid areas, such regions would benefit from the outcomes of water managing research. Inquiries related to selecting suitable locations for RWH found that taking into consideration both the physical and socio-economic criteria is important for RWH. e past studies discussed earlier required considerable data availability and accessibility to achieve their objectives. Considering the limited data availability and severity of desertification in the study area of the west desert of Iraq, the main objective of this study is to propose a comprehensive and cost-effective method for identifying and minimizing the impacts of desertification through the combination of geospatial data with GIS techniques. Rapid development in satellite image technology and GIS means that these components are easily accessible at reasonable costs. is factor is highly important, especially in areas where data are scarce or unavailable. e aim of this study is to use the Land Degradation Index (LDI) and Normalized Differential Vegetation Index (NDVI) to estimate the potential risk of desertification, identify the most vulnerable areas, and determine the most appropriate sites for RWH. e proposed method may help decision-makers and planners for a quick reference and determination of the potential risk of desertification, as well as prioritize for determining suitable sites for RWH to minimize the risk of desertification.

Study Area
Wadi Al-Gahdaf is situated in the western desert of Iraq.
is study area is situated to the east of the Euphrates River, between 32°20′ 14″ to 33°08′ 00″ N and 40°01′ 00″ to 43°30′ 00″ E and has a catchment area of 8513 km 2 ( Figure 1). e catchment experiences an arid environment with cool winters and dry summers. It is classified as a desert due to the low mean annual precipitation of between 100 and 150 mm. About 48% of the precipitation occurs during the winter, 16% in autumn, and 36% in spring [24]. e average annual evaporation is 3300 mm, and the monthly mean evaporation fluctuates significantly with the seasons. Considerable parts of the study area are covered with very fertile soil, in which different types of vegetation grow, such as tulips, lilies, colocynth, and numerous types of flowers grow during rainy seasons. Many drought-resistant shrubs and spiny plants are grown all over the year also covering large areas, especially in wide valleys and depressions, such as colocynth and thyme. Exposed rock outcrops of hard limestone provide a sound base for construction of barriers or dams. e study area is considered to be a very important part of the western Iraqi desert because it represents one of main valleys that receives a significant amount of rainfall during the rainy season, and it also faces issues related to desertification. So, an effective solution to mitigate the problem of desertification is needed. An Earth dam is one technique to conserve rainwater during the rainy periods for use in dry periods, and this is considered an effective way for developing a sustainable water supply in Iraq.

Data Requirements
e data used for this study are minimal and can be divided into three major categories, that is, satellite images, meteorological data, and soil characteristics. Satellite images that include the operational land imager sensor (Landsat 8), thematic mapper sensor, and 30-meter shuttle radar topography mission (SRTM) digital elevation model were provided by the United States Geological Survey (USGS). e OLI and ETM + provide data at a spatial resolution of 30 m for visible, near infrared (NIR), and short-wave infrared, and of 15 m for panchromatic bands [25,26]. e study area coincided with the extent of Landsat images scene (path/row: 169/037, 169/038, 170/037, 170/038, 171/037, and 171/038). A total of 6 of OLI images and 6 of Landsat 7 ETM + over the study area were downloaded from United 2 Advances in Civil Engineering States Geological Survey (USGS) (https://earthexplorer.usgs. gov/). e images were acquired in dry and wet seasons. e images were geometrically and atmospherically corrected. According to UTM-WGS84 Zone 38 N projection, the coordinate systems of all the data used were unified. Longterm and reliable daily and monthly precipitation data of 34 years from 1980 to 2013 are also used in this study. Twentythree samples were collected throughout the study area, and the soil characteristics were determined in the laboratory. e main sources of data used in this research are presented in Table 1.

Methodology
e research methodology can be divided into two major sections: (1) an analysis of the change in the LDI and the NDVI using geospatial data in 1999 and 2014 to define the existing land degradation and (2) the determination of suitable sites for RWH to minimize the impact of desertification via Model Builder within ArcGIS 10.2. is includes the selection of suitable sites, classification of the appropriateness for every criterion, and GIS analysis for generation of suitability maps based on the combination of the weighted linear combination (WLC) and Boolean operators to order the alternative zones according to ranking with respect to evaluation criteria. e overall methodology is presented in Figure 2, while the specifics of each major section are given in the subsequent sections.

Determination of the Desertification Extent.
e NDVI and LDI were used as indicators of desertification phenomenon in the study area. ematic maps were developed for the period of 1999 to 2013, and detailed methodologies are given in the following paragraphs.

Normalized Difference Vegetation Index (NDVI).
e application of NDVI has included the impact of various atmospheric conditions on vegetation [27][28][29], drought monitoring [30], identifying crop production and yield [31], and studies on desertification [32,33]. Most of these studies have been limited to quantifying the degree of desertification using NDVI. is index is an effective tool for quantifying the risk and monitoring changes of desertification as it estimates the change in vegetation according to the spectral response of the land cover [34,35]. e ERDAS IMAGINE software was used to determine the vegetation density changes and zones for the period of 1999 to 2013. e NDVI was estimated using the following equations [36].
From the thematic mapper sensor: where B 3 is the spectral reflectance value of band 3 (0.63-0.69) μm, and B 4 is the spectral reflectance of band 4 (0.77-0.90) μm. From the operational land imager sensor: where B 4 is the spectral reflectance of band 4 (0.64-0.67) μm, and B 5 is the spectral reflectance of band 5 (0.85-0.88) μm. e NDVI values fluctuated between −1 and +1. Values greater than 0 represent existing vegetation. More specifically, values between 0.0 and 0.2 represent low vegetation areas, and values of more than 0.2 represent dense vegetation areas. NDVI values between −1 and 0 represent nonvegetated areas, that is, abandoned land and urban land. e density slicing method was applied to classify the NDVI range level to estimate the effectiveness of vegetation variations between 1999 and 2013.

Land Degradation Index (LDI).
e Land Degradation Index (LDI) reflects the condition of booths degraded land and land cover, and it has a clear physical and biological meaning [37]. e LDI is representative of the total degree of difficulty in re-educating degraded land in a specified region [38]. e LDI can be used to further enhance and demonstrate the desertification phenomenon in the study area. Many researchers [3,39] have applied LDI to assess and analyze the hazards of the extension of soil degradation. e LDI was used by Meng et al. [40] to determine the degradation of soil based on band 2 and band 3 of the thematic Advances in Civil Engineering mapper sensor and band 3 and band 4 of the operational land imager sensor. Each of these two bands was used because they give high reflectance for soil degradation. Meng et al. [40] suggested four classes of soil degradation, namely <30 (extreme), (30-60) strong, (60-90) moderate, and >90 light. e estimation of LDI using the thematic mapper sensor is given as where B 2 is the spectral reflectance of band 2 (0.52-0.60) μm and B 3 is the spectral reflectance of band 3 (0.63-0.69) μm. e estimation of LDI using the operational land imager sensor (Landsat 8) is given as where B 3 is the spectral reflectance value of band 3 (0.53-0.59) μm and B 4 is the spectral reflectance of band 4 (0.64-0.67) μm.

Identification of the Most Appropriate Site for RWH
Structures. Different criteria and constraints directly affect the process of decision-making regarding site selection. is study offers a robust method to facilitate this process combining the WLC with the Boolean operators to order the alternative zones according to favorites with respect to evaluation criteria. ese methods are the most mainly used decision judgment in GIS [41][42][43].

Weighted Linear Combination
Operation. e WLC integrates the maps by assigning a normalized score to every category and gives weight to the criteria themselves [44]. e WLC technique has been done in two steps: (i) multiplying the weight of each parameter by rates of category of the same parameter and (ii) adding all the criteria layers in the raster calculator.
e WLC method was adopted because of its flexibility in choosing suitable sites. Several researchers have adopted this method including [45][46][47][48]. Six criteria were considered in this study in the identification of the most suitable sites for rainwater harvesting, namely slope, rainfall and runoff, land use/cover, soil texture, remoteness from irrigated land, and remoteness from roads. e criteria are measured using a variety of scales and measurements. Each criterion was first classified individually using pixel values from 0 to 10, as given in detail in Table 2. e least convenient were categorized as 1, while the most suitable areas were categorized as 10. Weights, classification, appropriateness scales, and scores are based on previous studies performed in study areas with the same nature [45,46,[49][50][51].
For each criterion, the scaled map was produced with pixel amounts ranging between 0 and 10. e combined appropriateness map was created by merging these criteria layers using the raster calculator. e appropriateness amounts were then categorized into five categories namely very low, low medium, high, and very high suitability. e six criteria that are important for selecting the potential locations for RWH are discussed in the following ( Figure 3): (1) Runoff depth: it is an essential criterion for selecting a suitable site for RWH and is utilized to evaluate the possible water supply through runoff [20,52]. e curve number (CN), which is expected from the consequence of land cover and soil type on rainfall/ runoff relationship, was used to estimate the runoff depth [52][53][54]. ese data were manipulated within GIS using the spatial analysis model to generate a digital map of runoff depth in the study area based on the following equation: where Q represents the runoff depth (mm); P is the rainfall (mm); S represents the possible maximum detention after the onset of runoff (mm); and I a represents a primary abstraction (mm) that contains all losses prior the start of infiltration, runoff, evaporation, and water interception by vegetation and is specified as I a � 0.2S from analysis of rainfall information for small watershed [55]. Equation (5) can therefore be expressed as where S can be calculated using CN as CN lies between 0 and 100, and this value signifies the runoff constrained to a specified precipitation amount. High CN values show that a considerable proportion of the rain will develop surface runoff [56]. e CN map and rainfall data for the study area were used to infer the layer of runoff depth in terms of the hourly extreme rainfall of the whole research area (Figure 3(a)). e minimum and maximum values of the runoff depth vary between 7 and 35 mm, with runoff increasing in the downstream part. e first category of <20 mm, represents 17%, the second category of 20-25 mm represents 22%, the third category of 25-30 mm represents 32%, the fourth category of 30-35 mm represents 16%, and the final category of >35 mm represents 13% of the entire area.
(2) Slope: the generation of runoff also greatly depends on the slope which affects the velocity of flow, the amount of sedimentation, and the quantities of construction material needed to erect the RWH structure. Critchley et al. [57] recommended that an RWH structure is not suitable for the slopes exceeding 5% because erosion rates are higher and required larger earthworks [46]. A DEM generated from SRTM data for the research area was used to generate the slope map which was categorized into five categories (Figure 3(b)): the first category of <2%, represents 9%, the second category of 2% to 3% represents 9%, the third category of 3% to 4% represents 10%, the fourth category of 4% to 5% represents 10%, and the final category of >5% represent 62% of the entire area. (3) Soil texture: the soil texture category is determined by the percentage of clay, silt, and sand, which affects the rate of infiltration and surface runoff. Mediumand fine-textured soils are typically more suitable for RWH due to their improved water retention. Since the low permeability of clay and its ability to retain the harvested water, sites with clay soil proved to be the best for water storage [58]. erefore, the soil texture is a critical factor for selecting sites for RWH. A soil map was generated based on soil samples that were taken from 23 sites in the study area using an auger, with depths ranging from 20 to 40 cm. e particle size distributions were used to classify the soil samples according to the USGS method of texture classification. e spatial analysis tool with Advances in Civil Engineering  GIS platform was used to develop the soil map (Figure 3(c)). A large portion of the study area is covered by sandy clay loam, which represents fertile soil when water is available. As shown in Figure 3(c), the soil map was categorized into five classes, that is, the first category (loamy sand) represents 30%, the second category (sandy clay loam) represents 59%, the third category (silty clay loam) represents 5%, the fourth category (clay loam) represents 3%, and the final category (loam) represents 3% of the entire area. (4) Land use land cover (LULC): Kahinda et al. [59] found that denser vegetation represents the higher rates of infiltration and lower runoff. e land use/ cover map was extracted from Landsat 8 satellite image taken in June 2019. e land use/cover, as a result of supervised classification, categorized the study area into three main categories consisting of: grass and vegetation (17%), residential (2%), and barren land (81%) (Figure 3(d)). (5) Distance from irrigated land: it is crucial that the RWH structure sites be close to the cultivable zone [21]. erefore, the distance from irrigated lands was adopted in this study. e NDVI was used to identify the green area (irrigated/activity growing vegetation). A vegetation cover map was derived from Landsat 8 satellite image taken in April 2019. e vegetated cover represents 25% of the entire study area (Figure 3(e)). (6) Distance from roads: roads have a considerable socio-economic significance for the area's population in the study region by providing access to water and grass for their livestock. e existing roads close to the suggested site contribute to reducing transportation costs. People, there may drive their trucks and tankers across these roads from one location to another. e distance from the road was categorized into four classes 250 to 500, 500 to 1,000, 1,000 to 2,000, and >2,000 m (Figure 3(f )). Less than 250 m did not recommend prohibiting any future inconsistency between RWH structure and roads [50]. A map of this parameter was obtained, based on the buffer areas adopted in Table 1, from Landsat 8 images taken in June 2019.

Boolean Overlay Operation.
e Boolean overlay method separates the operation of site selection based on using either (the OR) or (the AND); it is crispy (false or true), and it is limited to a small discrete location. Many researchers [15,60] have adopted this approach. e Boolean overlay approach is an important method for eliminating a specific area identified by the WLC method. e stream order and distance from faults are adopted in this method. e details of the two criteria are discussed as follows ( Figure 4): (1) Distance from faults: faults are a major hindrance when determining an RWH system. e location of faults is determined by the use of an obtainable geological map of the research area. e exclusion of the fault area from the adopted sites is critical when identifying a site for RWH [46]. Consequently, a value of 0 is given to areas that are at a distance of 1,000 m or less from the fault area, whilst a value of 1 is given to areas more than 1,000 m away (Figure 4(a)). (2) Stream order: the stream order signifies the hierarchical linking between stream sections and allows the categorization of drainage watershed via their size. e suitability of RWH depends strongly on the wadi density, and highly dense regions are the most appropriate. Analysis of stream order is very important for site selection for RWH; a lower stream order has a higher infiltration and permeability. A stream order larger than a fourth order is given the value 1, while a stream order smaller than a third is given the value 0 in this study [52]. e thematic map of the stream order is classified as very low (2), low (3), moderate (4), high (5), and very high (6). Higher stream order (sixth) was showed in the north east part with a total length of 214.7 km; the total length of the fifth order is 152.2 km; the total length of the fourth order is 296.4 km; the total length of the third order is 643 km (Figure 4(b)).

Development of Potential Sites for Rainwater Harvesting.
To identify the possible sites for RWH, the vector and raster GIS database was establish using Model Builder in ArcGIS 10.2. After that process was completed, suitability maps were developed. e appropriate areas for RWH were determined by rescaling all layers of WLC criteria, that is, slope, rainfall and runoff, land use/cover, soil texture, distance from irrigated land, and distance from roads and combining them by spatial analyst tool with the platform of ArcGIS 10.2. Spatial analyst tools with raster calculators were used in the ArcGIS platform for the Boolean criteria, that is, stream order and distance to the faults to obtain the sites with values of 1 which were suitable for RWH. Following this, the land suitability of the combined WLC and Boolean methods was categorize based on the NDVI and LDI to determine the areas that should be considered first and that were most appropriate for RWH structures.

Results and Discussion
e application of NDVI is an effective tool for quantifying the risk and monitoring changes of desertification, which include the impact of various atmospheric conditions on vegetation, and drought monitoring. e density slicing method was used to classify the NDVI change from 1999 to 2013. e results show that the area with a significant change of NDVI is 6,981 km 2 , which represents 82% of the total study area, while the area with a low change of NDVI is 1,532 km 2 , which represents 18% of total study area. e rate of deterioration of vegetation cover reached 18% in less than 15 years. erefore, the percentage of land classified as nonvegetation increased, and it is expected to accelerate as a Advances in Civil Engineering result of the desertification phenomenon in the study region ( Figure 5).
Other than NDVI, LDI could also be used to further enhance and prove the desertification phenomenon in the study area. e findings of the changes in LDI were mostly classified as extreme degradation (766 km 2 , represents 9%), strong degradation (5363 km 2 , represents 63%), followed by moderate degradation (2299 km 2 , represents 27%), and none degradation (85 km 2 , represents 1%) for the period of 1999 to 2013 ( Figure 6). e change in land degradation increased as a result of the desertification phenomenon in the study region.
As a result of the NDVI and LDI analysis, the percentage of land classified as nonvegetation increased, and it is expected to accelerate as a result of the desertification phenomenon in the study region. erefore, RWH is a significant component for minimizing the risk of desertification.
is study aimed to simplify the process of determining suitable sites based on a combination of the WLC method and the Boolean overlay to order the alternate zone according to favorites with respect to valuation criteria. After determining the watershed properties in the stage of GIS, the suitability sites for RWH were determined using multicriteria evaluation for six criteria, that is, the runoff depth, soil texture, slope, land use (LULC), distance from roads, and distance from irrigated lands, adopted in the WLC step. All layers of the criteria were combined in the raster calculator to gain a definitive map of the WLC criteria (Figure 7). e results map show that the maximum value calculated through the WLC way is 46, while the minimum value is 16. e differences between the maximum and minimum value of the WLC method can be divided into three classes. For this work in study area, the values that ranged between 36 and 46 represents suitable area (16%), the values between 26 and 36 represents moderately suitable areas (60%), and the values between 16 and 26 represents unsuitable areas (24%) for the entire study area (Figure 7).
e map of WLC criteria shows that the most downstream zone of the study area is suitable for RWH. e main soil texture in the suitable regions was clay, clay loam, and silty clay loam, and was intensively cultivated. Moreover, the regions suitable for RWH, had a runoff depth diversity of between 30 and 40 mm and slopes varied between 2% and 4%.
e results map is in agreement with the studies of Adham et al. and Mbilinyi et al. [20,58].
ese studies demonstrated that areas with fine to moderate slopes incorporated with clay and clay loam were appropriate for RWH.
In the western desert of Iraq, the scarcity of surface water resources, and the high cost of groundwater investment because of the 200 m depth in some region, require special effective water resources management [9,21]. An Earth dam is one method to conserve rainwater during rainy periods for usage in dry periods, and it is considered an effective way for developing a sustainable water supply in Iraq. erefore, the Boolean overlay method was used to limit the potential sites for RWH throughout the seasonal drainage network. Based on the faults distributed in the entire study area and stream order, a finding map of the Boolean overlay was generates as shown in Figure 8. is is a significant method for excluding certain areas identified by the WLC process. e findings of Boolean overlay, for which the value 1 represents 35% and the value 0 represents 65%, show that the very high potential for RWH sites represents the highest stream (sixth order), which was showed in the outlet of the catchment of the study area.
e findings map of WLC criteria, and Boolean overly criteria were multiplied to yield the concluding raster map (Figure 9). e resulting map included three units that have been used as markers for possible RWH sites due to their suitability: suitable, moderate, and unsuitable. ese units meant that only 2% of the catchment is suitable, 24% of the catchment is moderately suitable, and 74% of the catchment is not appropriate.
According to the analysis of the potential sites for RWH in the catchment area with the NDVI and LDI changes, it was demonstrated that 65% of highly suitable sites are in areas with a large change of NDVI, and 35% of highly suitable sites are in areas with a small change of NDVI, while 90% of moderately suitable sites are in areas with large change of NDVI, and 10% of moderately suitable sites are in areas with a small change of NDVI, 85% of highly suitable sites are in areas with a moderate change of LDI, 12% of highly suitable sites are in areas with a large change, and only 3% of highly suitable sites are in areas with no change.
In comparison with the literature, the adaption of together the WLC and Boolean methods within the GIS environment and the analysis of NDVI with LDI changes can enable hydrologists, decision-makers, and planners to quickly determine and minimize the risk of desertification and prioritize for determining suitable sites for RWH. Moreover, this methodology can help decision-makers and planners to estimate the potential risk of desertification, identify the most vulnerable areas, and determine the most proper sites for RWH.

Conclusion
It is very hard to revitalize a region that has a high probability of suffering from partial or total desertification. erefore, it is more efficient and economical to protect the land before its degradation and minimize the negative impacts of desertification. is can be done in several ways, including by maintaining water resources, exporting these resources, and using technologically advanced methods to protect them from total loss. is study has concentrated on the usage of RS with GIS to minimize desertification in the western desert of Iraq, with a consideration of the existent soil degradation condition and estimated possibility sites for RWH. e outcomes show that the changes in LDI with extreme, strong, and moderate degradations increased by 9%, 63%, and 27%, respectively, for the period 1999 to 2013. In addition, significant change of NDVI covered 6,981 km 2 , that is, 82% of the total study area, while the areas with a marginal change of NDVI was 1,532 km 2 (18% of the whole study area). e rate of deterioration in the vegetation cover reached 18% in less than 15 years. e results of the combined rainwater harvesting suitability map and the maps of vegetation and degradation changes found that 65% of highly suitable sites are in areas with a large change of NDVI, and 35% of highly suitable sites are in areas with a small change of NDVI, while 90% of moderately suitable sites are in areas with large change of NDVI and 10% of moderately suitable sites are in Advances in Civil Engineering areas with small change of NDVI. Moreover, 85% of highly suitable sites are in areas with a moderate change LDI, 12% of highly suitable sites are in areas with a strong change, and only 3% of highly suitable sites are in areas with no change. e appropriateness map is beneficial to decisionmakers, and hydrologists for a quick reference and determination of the potential risk of desertification, as well as the identification of a proper sites for RWH to minimize the risk of desertification. e proposed method is effective in terms of time and cost, particularly for extent study areas where data availability is scarce. Validation of this approach was proved using the nature and characteristics of the study area.