Next Article in Journal
Future Wave Climate-Driven Longshore Sediment Transport and Shoreline Evolution along the Southwestern Black Sea
Previous Article in Journal
Integrated Struvite Precipitation and Fenton Oxidation for Nutrient Recovery and Refractory Organic Removal in Palm Oil Mill Effluent
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rainwater Harvesting Site Selection for Drought-Prone Areas in Somali and Borena Zones, Oromia Regional State, Ethiopia: A Geospatial and Multi-Criteria Decision Analysis

by
Amanuel Kumsa Bojer
1,*,
Desta Jula Bekalo
1,
Taye Girma Debelee
1,2,
Saralees Nadarajah
3 and
Ayad M. Fadhil Al-Quraishi
4
1
Ethiopian Artificial Intelligence Institute, Addis Ababa P.O. Box 40782, Ethiopia
2
Department of Electrical and Computer Engineering, Addis Ababa Science and Technology University, Addis Ababa P.O. Box 16417, Ethiopia
3
Department of Mathematics, University of Manchester, Manchester M13 9PL, UK
4
Petroleum and Mining Engineering Department, Tishk International University, Erbil 44001, Kurdistan Region, Iraq
*
Author to whom correspondence should be addressed.
Water 2024, 16(13), 1789; https://doi.org/10.3390/w16131789
Submission received: 24 April 2024 / Revised: 24 May 2024 / Accepted: 27 May 2024 / Published: 25 June 2024
(This article belongs to the Section Urban Water Management)

Abstract

:
Rainwater collection systems play a crucial role in enhancing water availability in regions with low precipitation. This study focused on identifying potential rainwater harvesting (RWH) locations in drought-prone areas, specifically the Borena zone of Oromia and the regional states of Somali. This research leveraged geospatial techniques and a multi-criteria decision analysis (MCDA) to assess feasible RWH sites. The dataset comprises essential factors such as rainfall, drainage density, slope, soil texture type, and land use/land cover. These thematic layers serve as inputs for analysis, with each factor being weighted using the Analytical Hierarchy Process (AHP) method based on its significance. Reclassifying factors into subclasses facilitates suitability analysis. The weighted linear combination (WLC) technique is applied to identify and prioritize potential rainwater harvesting (PRWH) locations based on four suitability classes: highly suitable, moderately suitable, low suitability, and unsuitable. Our findings reveal that 1% of the study area, covering approximately 3288 km2, is highly suitable for RWH. Areas with moderate suitability constitute approximately 12% (37,498 km2), while regions with low suitability, representing the majority, encompass about 75% (242,170 km2). Additionally, 13% (41,000 km2) of the study area is deemed unsuitable for RWH. The proposed technique for identifying suitable RWH sites is adaptable to other low-precipitation regions. However, before implementing RWH structures, further research is imperative. This study proposed the exploration of socioeconomic variables in future research and urged for an in-depth examination of various aspects of environmental sustainability. Our research paves the way for adapting rainwater harvesting systems to align with community needs and life cycles while also exploring the socio-economic and environmental dimensions of sustainability for future study. The insights offer promising solutions to address the urgent issues associated with water scarcity. This should include comprehensive site depictions, an exploration of social and economic activities, and the meticulous preparation of a cost-benefit analysis.

1. Introduction

Rainwater harvesting involves collecting and maintaining rainwater before it is dispersed as surface runoff and is a key water intervention system. Effective water management, especially in arid and semi-arid regions, is a critical concern for human existence [1,2]. Addressing this global challenge requires the development of innovative solutions to ensure sustainable water management. Out of the total amount of water in the world, only 1% of freshwater is available for human consumption, and resources are also threatened by a multiplicity of problems, including climate change, overdevelopment, and water pollution, despite their significance to natural existence [3,4,5]. The challenge of water scarcity has been intensified by climate change, global warming, and population expansion, leading to critical shortages of water resources worldwide [6,7]. In response to this, rainwater harvesting (RWH), an age-old practice involving the collection, storage, and reuse of rainwater, has regained prominence as a viable solution for water supply in various regions across the globe [8]. Water is the most significant component of all food availability and nutrition determinants, such as agricultural production, energy, poverty reduction, economic growth, and survival [9,10]. Water resources and rainwater, in particular, are renewable sources and are some of the most important resources for the development and livelihood of societies with limited water capital in low precipitation areas. Rainwater also increases groundwater potential if properly harvested [11,12].
Ethiopia has a variety of topographies, with climatic differences between regions, unequal temperatures, and precipitation. High rainfall and humidity are prevalent in the south and southwest of the country, whereas the northeast, eastern, and southern lowlands are characterized by arid and semi-arid climates with very little rainfall. The Bale Mountain in the south of the country is characterized by Afro-alpine vegetation. The central and northern highlands of the country experience cooler climates [13]. Despite there being several major perennial rivers and lakes in Ethiopia, uneven distribution and inconsistencies in rainfall have affected the country on several occasions, particularly in areas below 1500 m [14]. There is substantial interannual variability in rainfall across the country [15,16]. Low precipitation has affected rural communities in Ethiopia’s lower land area for a decade. Despite persistent low precipitation, water remains vital for communities in this region, which depend heavily on livestock for their livelihoods. These communities not only depend on cattle and other farm animals for income through sales, but also for essential resources such as milk and meat for their dietary needs.
From 1980 to the present, rainfall data for southeast Ethiopia show an overall decline in rainfall between March and September. Across the country, this decrease has led to intense and frequent droughts [17]. In the study region, recurrent drought events, both historical and recent, have inflicted considerable damage, leading to livestock fatalities and crop losses. As a result, drought-prone areas of eastern Somali and the Borena zone of the Oromia region have been severely affected by droughts since 2020, and a low rainfall record for the third consecutive year has resulted in a high number of livestock depletions and food insecurity, particularly in the Somali and Borena zones. Droughts have serious implications for the lives and livelihoods of approximately 1.8 million people [18]. Areas with an average altitude of less than 1500 m have a warm, hot climate where rainfall is deficient (almost all of the Somali and Borena zones of the Oromia region). In these areas, situated on the peripheries, most of the land is arid and semi-arid, and most people living in the region are either pastoralists or semi-pastoralists. These areas constitute 10 and 12% of the country’s population [19]. RWH techniques are increasingly being recognized as valuable strategies to mitigate water scarcity in drought-prone regions [20,21]. This is particularly relevant for arid and semi-arid areas like Ethiopia’s Borena zone and Somali region, where rainfall variability and limited water resources pose significant challenges [22]. Insufficient rainfall during the growing season has resulted in crop failure across numerous agricultural lands. The scarcity of water has posed challenges for local villagers, particularly during the summer of 2022, when water shortages significantly affected animal husbandry.
RWH presents a compelling approach to managing agricultural water resources in regions prone to drought [23]. It enables farmers to capture and store rainwater during periods of high precipitation, making it available during critical dry spells [24,25]. Dams can hold water for later use, and their auxiliary reservoirs can protect against extreme weather events that may be detrimental to livestock, irrigation, fodder, or domestic consumption [26]. The best solutions to ensure improved access to and availability of water from rainfall for communities’ livelihoods and to improve the quality of life within regions are selecting rainwater harvesting sites and building structures [27]. RWH encompasses all techniques for collecting rainwater from rooftops or ground surfaces. This collected water can then be used for various purposes, including agriculture, domestic needs, or even drinking. However, identifying suitable locations for RWH structures presents a complex challenge [28]. The selection of suitable sites for rainwater collection necessitates the consideration of specific criteria that account for the geoenvironmental characteristics of the area [29]. The integrated RWH technique and storage facilities within landscapes in a planned and systematic manner are crucial. Therefore, suitable areas should be established for harvesting water and constructing structures to reduce vulnerability to droughts and seasonal fluctuations in rainfall [30]. Rainwater harvesting is crucial for reducing runoff, protecting soil erosion, and increasing groundwater resources through infiltration into the ground, ensuring a continuous supply throughout the year [31].
Rainwater harvesting holds immense significance, particularly in drought-prone areas of Ethiopia, where water scarcity poses substantial challenges. The implementation of rainwater harvesting techniques serves as a sustainable solution to augmenting water resources, addressing the pressing issue of water availability during extended dry periods. In Ethiopia, characterized by erratic rainfall patterns and recurring droughts, rainwater harvesting emerges as a critical strategy to enhance water resilience and alleviate the impacts of water scarcity on agriculture, livelihoods, and ecosystems. Studies, such as those conducted in [32], emphasize the positive impact of rainwater harvesting on water security, crop production, and community well-being in drought-prone regions of Ethiopia. The adoption of these techniques not only contributes to local water self-sufficiency but also aligns with broader water management strategies, fostering sustainable development in the face of climatic uncertainties [33].
There is a need to introduce an innovative and cost-efficient water-saving alternative to mitigate the problem of drought caused by low precipitation and the uneven distribution of rainfall. The assessment of RWH potential and the identification of suitable sites is a demanding and time-consuming task for water resource managers and planners, especially when dealing with large-scale applications. The selection of a potential site for water harvesting structures requires a multi-factor disciplinary approach, in which the use of remote sensing data with geospatial technique applications has become common practice in recent times [34,35]. However, advancements in geographic information system (GIS) and remote sensing (RS) technologies have streamlined and expedited the site selection process for RWH structures. These technologies have not only reduced the number of recommended sites but have also facilitated the identification of the most optimal locations. Different methodologies have been developed to select appropriate RWH sites. Hydrological modeling techniques enable the assessment of water availability and flow patterns, while field surveys provide on-the-ground insights into local conditions and community needs [36]. Also, MCDA methods, such as the AHP or Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS), can be employed to prioritize potential rainwater harvesting sites based on multiple criteria, including water availability, land suitability, and socioeconomic factors [37]. Currently, remote sensing (RS) and geographic information systems (GISs) stand out as highly valuable tools for the management of ecosystems and natural resources [38]. Also, participatory approaches involving stakeholder engagement and community consultations contribute to the identification of culturally appropriate and socially acceptable RWH sites [39].
Furthermore, the application of MCDA plays a pivotal role in the determination of appropriate zones for RWH [40]. The fusion of MCDA with GIS involving the integration of spatial data layers has been extensively employed in the RWH process [41]. The preferred method for choosing suitable sites for RWH in small geographical areas is field surveys, while another method also plays a role in determining which locations are appropriate. Such a method may include the main elements that play essential roles in finding suitable areas for RWH: rainfall, land use/land cover (LULC), a Digital Elevation Model (DEM), topography, soil texture, and environmental factors [42]. RWH projects can be significantly enhanced by employing a multi-criteria analysis (MCA) approach for site selection [43,44]. This technique involves systematically evaluating potential sites based on a defined set of criteria that influence the success and sustainability of the RWH system.
Below-average and uneven rainfall distributions have repeatedly affected the study area, resulting in drought and floods. However, the RWH site selection technique as a foundation for building structures must be considered by analyzing factors that have never been applied. This research aims to develop a robust methodology for creating a suitability zone map for rainwater harvesting (RWH) in drought-prone regions, specifically focusing on the Borena zone and Somali region. The proposed approach will be customized to align with the prevailing environmental, economic, and social conditions of the target region. Moreover, this study aims to identify optimal locations for installing various RWH systems. These systems encompass a variety of water management structures, such as waterholes, potential rainwater harvesting zones, and ponds designated for agricultural and livestock utilization. To achieve this, this study will leverage advanced technologies, including RS, GIS, and MCDA techniques. These tools will be employed to systematically identify potential RWH sites in the study area. Finally, this research aims to employ geospatial techniques and a multi-criteria decision analysis to discern viable water harvesting sites in drought-prone areas of the Borena zone in Oromia and Somali regional states.

2. Materials and Methods

2.1. Study Area

The Somali and Borena zones in the Oromia region are situated in the southern and southeastern parts of the country (Figure 1). They lie between latitude 3.40° to 9.11° N and longitude 36.64° to 47.94° E. The study area is approximately 323,955 km2, accounting for nearly 1/3 of the country’s total area. Dry and partly arid climatic conditions dominate it; most of the population’s livelihoods are nomadic and based on livestock production. The Somali regional state alone is the second-largest region, following the Oromia region of Ethiopia. The Wabe-Shebeli and Juba rivers are the only perennial rivers in the region. 90% of their flow originates within the study area and flows to Somalia in East Africa. The altitudes of the area extend from 163 to 2459 m above sea level, with the highest altitude being recorded in the northern parts of Somalia and the areas of Borena in the Oromia region. Most areas (approximately 77%) are below 1000 m in altitude.
According to [45], Ethiopia is mainly divided into five traditional agroecological zones, of which three are located in the area based on the altitude range: “Bereha” (hot-arid lowlands), which is situated at less than a 500 m altitude, characterized by very limited crop production or not at all; “Kolla” (hot lowlands), which is located between 500 and 1500 m in altitude and experiences crop production such as sorghum; And “Woina dega” (midlands), which is situated between 1500 and 2300 m in altitude above sea level. The altitude range is good for crop production, such as teff, wheat, and barley, but a very small proportion of the study area participates in growing such crops [46]. When looking into the soil textural properties of loam, loamy sand is the dominant coverage, followed by sandy loam and sand, and a small proportion of clay and clay loam exists.

2.2. Data Source

The central data sources for this study are satellite images, meteorological data, and soil characteristics. Satellite images were used to produce slope and drainage density data using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) DEM and Shuttle Radar Topographic Mission (SRTM) [47,48] at a 30 m resolution. The data were downloaded from an open-source geo-database of the United States Geological Survey (www.usgs.gov). Rainfall records that were converted to an annual average precipitation of ten years from 2012 to 2021 were collected from the Ethiopian metrological service based on 15 meteorological stations. LULC data were downloaded from the open-source dataset “Ethiopian Sentinel 2: LULC 2016” [49] and trimmed for the research area. Soil characteristic data for the study were adopted from [50,51].

2.3. Methodology

Identifying suitable areas for PRWH involves multiple criteria for a multi-factor problem. This study adopted a GIS-based MCDA model. Thematic input databases, such as rainfall, drainage density, slope, land use, cover, and soil texture, were prepared for suitability analysis. The ArcGIS 10.8 platform is ideal for locating the appropriate areas and suitable levels for potential RWH site selection. This study employed a multi-step approach for RWH site suitability mapping in the study area. First, relevant factors influencing RWH site selection were identified. Suitability levels were then assigned to each factor based on established criteria. Then, thematic layers were generated using spatial data for these factors. Weights were assigned to each factor based on its relative importance for RWH, followed by a weighted overlay analysis. Finally, the integrated weighted thematic layers were synthesized to create an RWH suitability map, delineating the most favorable locations for implementing RWH systems (Figure 2).

Parameter Selection

The identification of thematic layers for RWH site selection is determined by several factors: a literature review of earlier studies on RWH, climatic conditions, physical characteristics, the socioeconomic state of the area, and the availability of data [52]. In previous related studies, a large percentage of the parameters were identified, and most of the studies focused on climate and topographic characteristics, such as rainfall, slope, soil type, and LULC [53,54]. Hameed [55] reviewed 48 studies; in most of them, the parameters indicated were the factors that are frequently utilized to identify probable RWH sites.
This study utilizes geospatial MCDA methods to select optimal RWH sites. Environmental criteria that are crucial to RWH suitability, such as rainfall patterns, soil texture, slope, LULC, and drainage density, are identified and assessed. Thematic layers depicting the spatial distribution of each criterion are developed using geospatial data and RS techniques. Additionally, weights are allocated to these criteria according to their significance in RWH site selection. Through a weighted overlay analysis, the thematic layers are integrated to produce a comprehensive RWH suitability map. This method identifies areas with the greatest potential for successful RWH implementation by collectively evaluating the influence of these environmental factors. This method facilitates the identification of optimal locations for RWH implementation, promoting efficient water resource utilization and mitigating drought impacts.

2.4. Data Processing

2.4.1. Rainfall

Rainfall is the key parameter that affects RWH storage and receives the maximum weighting factor when deciding where to store it [56,57]. RWH site selection largely depends on rainfall variability and the volume of water collected in a suitable space. Annual rainfall data for ten years (2012–2021) were provided by Ethiopian Metrology Services for 14 meteorological stations inside and outside the study location. The average annual rainfall was again averaged based on the location of the meteorological stations, as shown in Figure 3, and latitude and longitude points were created using ArcGIS 10.8, ArcPro 2.9, and the QGIS 3.4 environment to assign the average rainfall data to a rainfall distribution map (Figure 3). This process aimed to assess and quantify the amount of rainfall in the study area, as this information constitutes a critical meteorological parameter. Annual rainfall data collected from the Ethiopian Meteorological Institute (EMI) were processed in Microsoft Excel for cleaning and analysis. Subsequently, the data were loaded into ArcMap 10.8 to generate a rainfall map, aiding in understanding spatial variability and informing RWH site selection with sufficient precipitation. The geostatistical method known as Inverse Distance Weighted (IDW), a widely accepted approach for interpolating rainfall variables, was employed to estimate the spatial distribution of rainfall across the region. This method relies on the geographical direction and distance of existing data points to approximate rainfall values for unobserved locations [58]. The careful utilization of these tools and methods enhances the accuracy and reliability of the rainfall data analysis, contributing valuable insights to the study’s findings. Spatial variability in rainfall is a key factor in the selection of RWH sites. Conducting an analysis of meteorological station data allows us to visualize this variability through rainfall distribution maps. These maps are crucial for identifying areas with consistent precipitation, a prerequisite for successful RWH implementation (Figure 3).

2.4.2. Slope

The slope has an essential influence on runoff generation, water infiltration, the modulation of sedimentation, and the velocity of water flow. This is a valuable parameter for the selection of RWH sites. A slope greater than 8% may not be suitable for RWH, whereas areas with low slopes are suitable for RWH site identification [59,60]. ArcGIS was used as a spatial tool for the processing of DEMs at a 30 m resolution. A DEM map of the study area is shown in Figure 4. The slope tool in ArcGIS 10.8 was used to create the slope map. To determine the suitability level for RWH site identification, the slopes derived in the study area were subdivided into percentage categories.

2.4.3. Soil Texture

Soil texture plays a critical role in RWH site selection effectiveness due to its influence on the infiltration rate and water-holding capacity. This characteristic is determined by the relative proportions of mineral particles (sand, silt, and clay) within the soil. Finer-textured soils with higher clay content exhibit greater water-holding capacity, making them more suitable for RWH site selection [61,62]. Conversely, coarser-textured soils with higher sand content may allow for faster infiltration but have a lower capacity to retain water, potentially limiting their suitability for RWH applications. When the ratio of clays is increased, the pore space is reduced, and thus, there are also reductions in infiltration. This limits water movement through the soil, improves drainage, and is useful for RWH [63,64].
The soil map utilized in this study was sourced from the Food and Agriculture Organization (FAO), and the global soil gridded information is available at https://www.isric.org/explore/soilgrids accessed on 20 April 2024 [65]. The country-level soil texture map was georeferenced and converted to the GCS_WGS_1984 Coordinate Reference System (CRS), which was also clipped to the study area and digitized on-screen using the GIS platform for spatial analysis (Table 1). According to [66,67], the soil texture types were reclassified into four soil texture categories based on similar properties. This reclassification was also in line with the classification of the hydrologic soil category of soil conservation services [68]. The SCS was also separated into one of the four HSGs, as shown in Table 1. Soil groups can be used to estimate rain overflow precipitation generation and identify rainwater storage sites [69].

2.4.4. LULC

Land use types and coverage influence the generation of overflow, volume, and velocity of water that streams to the lesser slope zones [70]. Areas with higher plant coverage tended to have a higher penetration rate owing to obstruction by dense vegetation coverage. However, areas with thin vegetation cover increased the overflow volume [71]. The Ethiopian Sentinel-2 LULC 2016 dataset initially consisted of six land use classes. These classes were then reclassified into four groups based on their suitability for rainwater harvesting RWH. The reclassified groups include shrubland, bare land, grassland, woodland, cropland, forest, and water bodies. The LULC analysis in the study area involved cross-referencing data through the use of Google Earth maps, randomly selected points, and on-site ground observations. This approach was employed to precisely define the different land use classes within the satellite imagery that had been classified. This reclassification was performed using spatial analysis tools, considering the suitability of land use/cover for RWH.

2.4.5. Drainage Density

Drainage density is defined as the total stream length per unit catchment [72]. Regions characterized by an elevated drainage density typically exhibit increased runoff suitable for harvesting [73]. Sites with high drainage density, characterized by lower-order streams, are particularly suitable due to their high infiltration and permeability rates [74]. This is because lower-order streams often have less channelization, allowing water to infiltrate the ground more effectively. Groundwater availability, the capacity of overflow loss due to penetration, and most hydrological processes are mainly influenced by drainage density; the lower the drainage density, the lesser the possibilities for RWH, while a higher drainage density has a higher probability of RWH [75,76]. RWH proves to be more effective in regions characterized by higher drainage densities, as it establishes a framework that facilitates the swift capture and retention of runoff [77]. The drainage density was derived from the DEM using ArcGIS 10.8 environment hydrological tools. Drainage density layers were generated through a series of steps involving DEM processing. These steps included filling depressions in the DEM, determining flow direction, delineating streams, and calculating line density.

2.5. AHP and Weight for Different Factors

AHP procedures were used to calculate the weight of each factor and identify an appropriate RWH location, followed by a GIS-based multi-criteria decision-making (MCDM) procedure. The primary objective of the AHP is to quantify the influence level of each criterion and show their relationship. An AHP analysis was conducted for each parameter and subsequent ranking, and the determined percentage of influence was then incorporated as input for the weighted overlay analysis within the ArcGIS software [78]. AHP prioritizes options based on biophysical (rainfall, topography, and soil), socio-economic (water needs and land ownership), and technical (technology type and storage capacity) factors. Pairwise comparisons conducted by experts help to determine criterion weights, ensuring that water needs, infrastructure, and sustainability are considered for optimal RWH site selection. The AHP technique is a structured system for forming and examining multifaceted decision making to govern suitable areas [79]. Thomas L. Saaty developed this technique from 1971 to 1975 [80,81]. For the location of proper RWH site selection, this method was employed to determine the significance of each factor in percentages. The AHP procedure was used to compare the factors identified using inverse matrices. Gradings from 1 to 9 were assigned to the criteria based on their importance for rainwater harvesting suitability. Pairwise comparisons ascertain the comparative significance of two factors and aid in assessing their appropriateness for a specific purpose. The evaluation involves a continuous 9-point scale, wherein the two factors are compared and assigned values (Table 2).
The lowest rank is equally important, whereas the highest rank is excellent for rainwater harvesting site selection. To assess the reliability of weights assigned to different factors in the matrix, the Random Index (RI) is determined based on the size of the comparison matrix (n), and its values are derived from simulations involving random matrices of the same size. The Consistency Ratio (CR) is calculated to assess the level of consistency in the expert judgments [82]. The CR is derived from the following formula:
CR = C I R I
where CI is the Consistency Index and RI is the Random Index.
The Consistency Index (CI) is calculated using the following formula:
C I = ( λ m a x n ) ( n 1 )
where λmax is the largest eigenvalue of the comparison matrix and n is the size of the comparison matrix. Once the pairwise comparison matrix is constructed, eigenvalues are calculated. An eigenvalue (λ) of a matrix is a special scalar value associated with a non-zero eigenvector (v) of the matrix. An eigenvalue represents a scaling factor that stretches or shrinks the eigenvector when multiplied by the matrix [82]. The equation for calculating eigenvalues can be written as follows:
A v = λ v  
where A represents the pairwise comparison matrix, v represents the eigenvector, and λ represents the eigenvalue. Solving this equation for λ requires specialized mathematical methods, often handled by software packages designed for AHP analysis. The RI serves as a standardized metric established by [83]. The RI is a benchmark for comparing the Consistency Index (CI) against random matrices. The RI value depends on the size of the matrix (n). It is derived from a large number of randomly generated pairwise comparison matrices and is given in a tabulated form for matrices of different sizes. This involves summing up all of the CIs and dividing by the number of matrices generated. The formula used to calculate RI is
RI n =     C I   r a n d o m N
where RIn is the Random Index for a matrix of order n, ∑CI random represents the sum of the consistency indices obtained from randomly generated pairwise comparison matrices of order n, and N is the number of randomly generated matrices used to compute the average.
The RI serves as a benchmark for evaluating the consistency of pairwise comparisons made by experts or stakeholders. It provides a standardized measure of consistency across different pairwise comparison matrices. In the calculation of the Consistency Index (CI), the RI is used as the denominator, serving to assess the reliability of the comparisons. The formula for both the RI and CI involves the term (n − 1), where n represents the size of the comparison matrix. This standardization ensures that the consistency of comparisons can be objectively assessed regardless of the matrix size. The variability of RI values is contingent upon the number of parameters delineated in Table 3. The importance of these criteria in identifying appropriate locations may vary according to their respective weights, and the decision-making process heavily relies on these criterion weights. In this study, the Analytical Hierarchy Process Consistency Ratio (AHP-CR) was employed to evaluate the consistency of the assigned weights across various layers. The Random Index (RI) value, a reference point for calculating the CR, can be found in Table 3 based on the number of compared parameters. By rigorously applying the AHP-CR, this study ensured that the assigned weights across various criteria were consistent and reliable. This methodological rigor strengthens the validity of the multi-criteria decision analysis, ultimately leading to more accurate and credible site selection for rainwater harvesting in the drought-prone Somali and Borena zones. Consequently, this approach supports better-informed and more effective interventions for water security in these regions.
The AHP is a decision-making framework used to simplify complex problems by organizing them into hierarchical structures comprising criteria, sub-criteria, and alternatives. It prioritizes these elements through pairwise comparisons. However, due to the potential for inconsistencies in such comparisons, the AHP incorporates a consistency ratio to verify the reliability of judgments. AHP pairwise comparisons assigned weights to RWH suitability factors. A ratio lower than a predefined threshold, typically 0.1, signifies acceptable consistency. According to Saaty [83], the value of CR should be less than 0.1 or less than 10%; otherwise, the weights should be re-evaluated to maintain uniformity. In our study, we identified 5 parameters, and according to Saaty’s CR threshold, comparisons are deemed sufficiently consistent when CR is less than 0.1. We proceeded to calculate the RI based on these 5 parameters. Specifically, for order 5, the calculation yielded CI = 1.11 − 5/5 − 1 = −3.89/4 = −0.9725. In this study, the AHP pairwise matrix for the factors indicated in Table 4 was calculated according to the importance of the factor assigned. The process involves systematically comparing each factor against every other factor to determine their pairwise importance. Participants in the decision-making process assign numerical values, typically on a scale of 1 to 7, to express the relative importance or preference of one factor over another. These judgments are then organized into a pairwise matrix, quantifying the importance ratio between each factor pair. This structured approach facilitates decision making by revealing the hierarchical relationships among factors.
Rankings from 1 to 7 were assigned to the criteria for the study, and for each criterion, the percentage influence was calculated, as shown in Table 5 Several factors were considered when ranking criteria for the RWH site selection in this study. Firstly, criteria directly impacting rainwater harvesting, such as the average annual rainfall or drainage density, were given higher rankings due to their direct influence on the amount of rainwater collected. Secondly, criteria affecting water storage and usability, like soil texture or slope, were ranked based on their ability to retain and make harvested rainwater accessible. Lastly, criteria influencing long-term sustainability, such as LULC, were ranked with a focus on minimizing issues like land degradation or competition for water resources, ensuring the viability of the site selection over the long term. Finally, weights were assigned to each criterion using CR, RI, and CI calculations, ensuring a systematic and comprehensive evaluation process. In addition, fractional values in the table may arise due to various reasons. Participants may detect slight variations in factor importance, resulting in fractional assignments. Equally important factors might receive numerical values reflecting their shared significance. Factors such as soil texture and slope play crucial roles in rainwater harvesting, each with a unique level of importance. Employing fractional assignments in the pairwise comparison matrix allows participants to precisely articulate these distinctions. In our study, we found that drainage density holds a moderate level of importance compared to LULC but exhibits greater significance than soil texture, and they were assigned values of 3 (moderate importance) and 5 (strong importance), respectively, reflecting these differences. To quantify the weight and influence of each factor and to standardize the comparison matrix, each column’s elements are divided by the sum of that column. This normalization process ensures that the values are appropriately scaled, facilitating fair comparisons across factors. For instance, consider the drainage density column in Table 5: the sum of its values is 1 + 1/7 + 1/5 + 1/3 + 1 = 1.796. Normalizing involves dividing each value in the drainage density column by 1.796. Following normalization, a consistency check is conducted to ensure the internal reliability of the pairwise comparisons. This check entails computing the CI and the CR. The CI is determined from the largest eigenvalue (λ_max) of the comparison matrix, while the CR is computed by dividing the CI by the RI. A CR below 0.1 indicates that the matrix is internally consistent and reliable. Finally, to normalize the weights, we summed the average weights and adjusted them to ensure they add up to 1. The sum of the drainage density column is calculated as follows: 1 + 1/7 + 1/5 + 1/3 + 1 = 1 + 0.1429 + 0.2 + 0.3333 + 1 = 2.67621 + 1/7 + 1/5 + 1/3 + 1 = 1 + 0.1429 + 0.2 + 0.3333 + 1 = 2.6762. We normalized the values by dividing each element of the column by the value of 2.6762. Drainage density: 1/2.6762 = 0.37371/2.6762 = 0.3737. Land use/land cover: 1/7/2.6762 = 0.00541/7/2.6762 = 0.0054. Soil texture: 1/5/2.6762 = 0.07471/5/2.6762 = 0.0747. Slope: 1/3/2.6762 = 0.12461/3/2.6762 = 0.1246. MAAR: 1/2.6762 = 0.37371/2.6762 = 0.3737.
The eigenvalue (λ) used in the calculation of the CR is derived from the CI and the RI. The RI is a reference value based on the order of the pairwise comparison matrix, which, in this case, is 5, as there are 5 criteria (drainage density, land use/land cover, soil texture, slope, and MAAR). The RI value for a matrix of order 5 is empirically determined to be 1.12. Considering a CI of 0.024 and an RI of 1.12, the pairwise comparison matrix is deemed consistent if the CR falls below 0.1. Since the calculated CR of 0.024 indicates a consistent matrix, this allows us to proceed with calculating the eigenvalue (λ) using the following formula:
CR = CI/RI. Given that CR = 0.0216, we rearrange the formula to solve for CI: CI = CR × RI, CI = 0.0216 × 1.12 and CI = 0.024192. Now that we have the CI value, we can use it to find the eigenvalue (λ) using the formula CI = λ − 1/n – 1, where λ is the eigenvalue and n is the order of the matrix (number of criteria). Given that n = 5 (since there are 5 criteria) and CI = 0.024192, we rearrange the formula to solve for λ:
λ = CI × (n − 1) +n, λ = 0.024192 × (5 − 1) +5, λ = 0.024192 × 4 + 5, λ = 0.096768 + 5 and λ = 5.096768. So, the eigenvalue (λ) derived from the given data is approximately 5.097.
The AHP model employs the pairwise comparison matrix from Table 4 to compute the weights for each criterion in Table 5. Each entry in Table 4 signifies the relative significance of the criterion in its respective row compared to the criterion in its corresponding column. For instance, to determine the weight for drainage density, we sum the values within its column in Table 4, resulting in a total of 1.796. Subsequently, we normalize the matrix by dividing each entry in the drainage density row by this sum, generating the values presented in Table 5. Then, to ascertain the weighted average for drainage density, we multiply each normalized value by its corresponding weight in Table 5 and sum the products. For example, the calculation proceeds as follows: weighted average for drainage density = (0.3737 × 0.3442) + (0.0054 × 0.3442) + (0.0747 × 0.3442) + (0.1246 × 0.3442) + (0.3737 × 0.3442) = 0.127 + 0.00186 + 0.02567 + 0.04273 + 0.127 = 0.32496. This procedure is replicated for each criterion to establish its weighted average, subsequently recorded in the “Average (Weight)” column of Table 5. Moving on to Table 6, the percent influences for each criterion are determined by dividing the weight of each criterion by the sum of all weights and then multiplying by 100 to express the result as a percentage. For instance, to compute the percentage of influence for drainage density, we perform the following calculation: percent influence for drainage density = (0.3442/1) × 100 = 34.42%. This process is iterated for each criterion, culminating in the values showcased in the “% of Influences” column. Rankings were assigned based on the relative priority of each factor Table 6. To calculate the weight for each factor, we can use the drainage density as an example: the sum of the drainage density row was 17, and the total column sum of all factors was 49.37. Therefore, the weight for drainage density was calculated as 17/49.37 = 0.3442. The value of CR = 0.0216, which is less than 10%; therefore, the assignment of weight above was consistent and acceptable (Table 6). After the weight of the criteria was determined, the potential RWH site was acquired using the weighted overlay process (WOP), intersecting all of the weighted criteria using the weighted linear combination (WLC) technique using the ArcGIS tools. This technique applies to suitability analysis and has been used in many previous studies [83,84]. A weighted linear combination is an analysis technique used when working on multi-criteria decision making (MCDM) or when considering more than one criterion [85,86]. A WLC suitability map was generated using the following equation:
S =     W i S L i  
where S is the suitability location and Wi is the relative importance weight of the criteria (rainfall, slope, soil texture, land use, land cover, and drainage density) for the RWH site selection procedure and SLi suitability levels for input factors.

2.6. Addressing Subjectivity in AHP for RWH Site Selection

This study addresses the inherent subjectivity of the AHP when used for rainwater harvesting site selection in drought-prone areas. Firstly, we advocate for a multi-stakeholder approach involving experts in hydrology, environmental science, community development, and local communities [84,85,86,87,88,89]. This diverse team can contribute specialized knowledge and potentially reduce bias by including various perspectives. Secondly, experts not involved in the AHP process can evaluate the final decision model and its outcomes, objectively assessing its validity and ensuring credibility [90]. Thirdly, we emphasize the importance of comprehensive transparency and documentation. This includes clear justifications for criteria selection, assigned weights, and final site selection results [91]. Finally, to reduce bias and ensure transparency, we propose establishing pre-agreed criteria weighting. This can be achieved through stakeholder consensus or expert judgment at the project’s outset [92]. Furthermore, we support the concept employing diverse criteria and sub-criteria to enhance the robustness of the AHP results. As [93] suggests, analyzing how changes in criteria weights affect the final decision strengthens the overall credibility and reliability of the process. Similarly, a sensitivity analysis, where weights and pairwise comparisons are slightly altered to observe the impact on site selection, can identify potential biases in the AHP judgments [94].

2.7. Validation

The validation process is a crucial step in ensuring the accuracy of the suitability map generated for diverse RWH site selection. This evaluation assesses the effectiveness of the methods and techniques employed. The determination of RWH suitability in arid and semi-arid regions, specifically within the Borena zone and Somali region, primarily relies on a combination of field studies, existing data, high-resolution satellite imagery from Google Earth, and ground truth information collection using optical satellite maps. Validation was conducted using ground truth data obtained through visits to selected RWH sites identified geospatially by MCDA. This validation process encompassed various aspects, including measuring actual rainfall, assessing soil characteristics, and evaluating slope conditions in comparison to the data utilized in the analysis. Furthermore, the presence of suitable locations for water storage and accessibility for the local population was examined on-site.
Additionally, potential limitations not accounted for in the geospatial data, such as land ownership issues or social factors, were assessed during field validation. This approach, as outlined by [95], involves leveraging geospatial technology. Moreover, ancillary data from various institutions are incorporated to validate and compare the study results against a classified and analyzed map of potential RWH suitability. The RWH suitability zone map is meticulously compared to the functional RWH structures within the study area. To achieve this, existing RWH structure points are sourced from reputable entities such as the Ministry of Water and Energy, the National Disaster Risk Management Commission, and the central statistical agency. By analyzing the spatial distribution of existing RWH structures relative to the suitability classes, this study evaluated the level of agreement between the model’s predictions and the on-ground reality. A high degree of overlap between existing RWH structures and the highly suitable zones on the map strengthens the model’s credibility and the potential effectiveness of the proposed RWH sites.

3. Results

Each parameter was subdivided into different importance levels for the RWH to produce the final potential suitability map. The importance level (suitability level) was divided into four categories: unsuitable, low suitability, medium suitability, and high suitability. In the results of this study, the basis for assigning limits to each category is primarily derived from a combination of expert knowledge, relevant studies in the literature, and an empirical data analysis. Existing research findings were carefully reviewed for each parameter, and consultation with domain experts was undertaken. Statistical analyses were conducted to establish thresholds that differentiate between various suitability levels. These thresholds are determined based on specific criteria relevant to each parameter and aim to capture meaningful distinctions in the suitability of the study area to carry out RWH. Furthermore, a systematic approach was employed to ensure consistency and objectivity in the categorization process. While the exact limits for each category may vary depending on the parameter and the specific characteristics of the study area, the aim is to provide a clear and scientifically justified framework for assessing suitability.

3.1. Rainfall

This study assesses rainfall levels in the area using Ethiopia’s agroclimatic zone classification system, focusing on the drought-prone Somali and Borena zones. The regional states of the Somali and Borena zones of the Oromia regional states of Ethiopia are more often characterized by semi-arid and arid climatic conditions. Prioritizing rainfall data is essential for RWH site selection in these regions. Annual rainfall in the area ranges from 294 to 601 mm. However, some areas in the north and southwest of the study area are characterized as sub-moist areas, where the annual RF is between 350 and 566. The annual average rainfall for the given stations for 10 years was produced using the Inverse Distance Weight (IDW) spatial interpolation method and clipped to the study area. The RWH site identification was determined based on the relative availability of rainfall. High-response areas showed the highest opportunities for RWH. Rainfall variability in the area ranged from 294 to 619 mm. A relatively high rainfall variability was observed in the north, southwest, and some points in the center of the study area, whereas the eastern and most of the central parts scored lower rainfall amounts.
A mean annual average rainfall > 500 mm was considered highly suitable for RWH, offering the potential for larger storage capacities and more frequent harvests (referred to as dry kolla). Sites with medium suitability, receiving 500–600 mm of rainfall annually (labeled as Bereha), can still benefit from RWH, albeit with potentially lower storage capacities and harvest frequencies [96,97]. Conversely, regions with less suitability, receiving 400–500 mm of rainfall annually, may face challenges due to limited water availability (Figure 5). However, a further analysis, such as an analysis examining the drainage density, could identify micro-locations suitable for RWH implementation within these zones. Regions receiving less than 400 mm annually are unsuitable for RWH due to insufficient rainfall and being highly prone to drought. The reclassification of rainfall showed that a 79,284 km2 (24.47%) area was classified as being unsuitable, a 73,110 km2 (22.57%) area was classified as having low suitability, a 104,666 km2 (32.31%) area was classified as having medium suitability, and a 66,896 km2 (20.65%) area was classified as having high suitability for RWH site selection.

3.2. Slope

A slope map was derived from a DEM with a 30 m resolution for the entire area. The GIS environment in ArcGIS was used to calculate the slope percentage or degree of slope. For this study, the percentage of slope (the percentage ratio of elevation change with horizontal distance) was calculated and used for the suitability analysis. Table 7 shows the percentage range of the slope calculated and the suitability levels of the slope for the RWH site selection procedure.
Figure 6 shows the distribution of the slopes and topography of the study area. The consensus among researchers worldwide underscores the significance of slope factors in determining suitable sites for RWH, a notion affirmed by our study’s classification system. Moderate slopes, typically falling within the 5–10% range, emerge as optimal locations for RWH because they facilitate runoff while mitigating soil erosion risks [98]. Gentle slopes, spanning 0–5%, are highly suitable for RWH, offering advantages such as efficient water collection, reduced erosion risks, enhanced groundwater recharge, and simplified construction processes. Conversely, steeper slopes, ranging from 10 to 15%, pose challenges, including heightened erosion risks and construction complexities. Such slopes are prone to soil erosion, potentially leading to sedimentation in RWH storage structures and diminishing their efficacy. Extremely steep slopes exceeding 15% are deemed unsuitable for RWH sites due to substantial erosion hazards and construction difficulties. Most of the study area, which is predominantly flat and gentle, covered approximately 98% of the entire area, and has medium and high appropriateness for RWH site identification, and less than 2% of the area has steep slopes. A smaller slope and flat area will slow down water movement and increase the chance of it staying on the soil’s surface. Most of the area is suitable for rainwater storage, but the porosity of the soil governs the duration of water holding in the area. Incorporating slope considerations into MCDA for RWH site selection in the Somali and Borena zones is imperative for supporting the effectiveness and sustainability of RWH interventions in alleviating drought impacts.

3.3. Soil Texture

Our study aligns with global research highlighting the importance of soil characteristics in RWH site selection through a soil data analysis. The soil texture within the study area was reclassified into four categories: clay, clay loam, loam, and loamy sand. Soil data are commonly utilized in site selection studies, with soils possessing good infiltration rates, such as sandy loam or sandy clay loam, generally considered more suitable due to their capacity to facilitate water percolation and groundwater recharge capabilities. Our findings present a unique case. Clay soil emerged as the most suitable type in this specific area due to its ideal balance between infiltration and water-holding capacity, aligning with previous studies [99]. However, clay loam, while moderately suitable, may require additional measures for optimal water capture due to its lower infiltration rates. Conversely, loam and loamy sand soils exhibit lower suitability due to very low infiltration rates, potentially causing surface runoff and limiting water storage within the soil profile. Sandy loam and sand are entirely unsuitable for RWH due to their rapid drainage and minimal water retention despite having high infiltration rates. Figure 7 depicts the spatial distribution of soil textures, with sandy loam and sand covering the majority of the study area (59%). Loam and loamy sand are less prevalent, followed by the most suitable clay soil type, which accounts for approximately 11% of the total area. Table 8 summarizes the soil types, their corresponding suitability classes, and their respective area coverage.

3.4. LULC

Ethiopian Sentinel-2 land use/land cover 2016 was the major input to drive LULC data, and reclassification was performed while considering the type of LULC potential for RWH. This study revealed that areas with minimal human intervention, such as shrublands or open woodlands, are well suited for RWH due to their sparse vegetation cover, which reduces interference, erosion, and contamination risks. Bare land, which is commonly found in drought-prone regions, was classified as moderately suitable for RWH due to its flat, non-vegetated nature, facilitating potential water collection, albeit with attention regarding erosion risks. Conversely, areas characterized by dense settlements, industrial zones, closed woodlands, or dense grasslands were considered less suitable for RWH due to their high interception, potential contamination, and increased runoff, limiting water availability [100]. Furthermore, croplands, forests, and water bodies were considered unsuitable due to concerns regarding agrochemical pollution and dense vegetation, which can interrupt rainfall. A comprehensive understanding of LULC dynamics is crucial for accurate RWH site suitability assessments in drought-prone areas. The integration of information on vegetation cover, impervious surfaces, agricultural practices, and land degradation into MCDA frameworks can greatly enhance the effectiveness of RWH site selection strategies, contributing to sustainable water resource management in the Somali and Borena zones of Ethiopia’s Oromia regional state. As shown in Table 9, shrubland constitutes the largest spatial coverage at about 44% (142,929 km2), followed by grassland and woodland with about 37% (121,030 km2); cropland, forest, and water bodies with 16% (53,193 km2); and bare land with the least spatial coverage, constituting only 2% (6803 km2) of the entire area.
In terms of suitability for RWH site selection, 44% of the entire coverage in the study area was determined to be highly appropriate because the area was covered by shrubland. Approximately 2% of the area covered by bare land was determined to have medium suitability. Grassland and woodland covered 37% of the area and was determined to have low suitability, and cropland, forest, and water bodies contributed 16% of the area and were determined to be unsuitable for RWH site selection (Figure 8).

3.5. Drainage Density

The drainage density of the area extended from 0 to 13. Table 10 lists the deranged density intervals and suitability classes assigned. The higher the drainage density, the more likely the RWH site is to be located, which means that when the area is dominated by rainwater accumulation in several streams, it is more likely to occur in an area with high drainage density and vice versa. The analysis of drainage density revealed its importance in RWH site selection for drought-prone regions. Areas with high drainage density (8–13 km/km2) were classified as highly suitable due to their potentially higher infiltration rates and higher surface runoff, leading to greater water storage within the soil profile and increased capture by RWH structures [101]. Conversely, areas with less drainage density (3–4 km/km2) were classified as being less suitable due to rapid surface water runoff and limited infiltration. Moderate drainage density (5–7 km/km2) areas might still be suitable with additional measures to improve infiltration, while very low drainage density (below 2 km/km2) areas were deemed unsuitable for RWH.
Figure 9 shows the drainage density and the reclassification map for suitability. Less than 2% of the area had the highest stream concentration, indicating higher suitability. The RWH covered 5864 km2 of the study area. Approximately 13% (42,017 km2) of the area was determined to have medium suitability, 39% (126,342 km2) was determined to have low suitability, and the remaining 46% (149,700 km2) was unsuitable for RWH site selection.

3.6. Final Potential RWH Suitability

The final RWH suitability map was produced using the weighted linear combination (WLC) method and weighted overlay process in ArcGIS. Weight and suitability rankings were assigned using a scale of 1 to 9, capturing the relative importance of factors in RWH site selection. These weights were derived from a comprehensive consideration of previous research, expert perceptions from hydrologists, environmental scientists, and community development specialists, as well as local knowledge facilitated through a multi-stakeholder approach. Further enhancement of the weightings was achieved through a pairwise comparison matrix and the AHP, with the CR ensuring their reliability. The thematic layers’ suitability levels were also determined on a scale of 1 to 9 using a weighted linear combination tool in ArcGIS. The suitability levels were determined using a weighted linear combination equation, assigning values 1, 3, 5, and 7 based on the importance scale. A suitability level of 1 indicated areas that were considered unsuitable for RWH site identification due to their lack of significance, and they were considered restricted zones. On the other hand, a suitability level of 7 represented areas that were highly suitable for RWH implementation. The importance scale delineated two distinct categories: areas considered practically unsuitable (rated 1 and 3) and areas deemed suitable (rated 5 and 7). This classification ensured clear differentiation between zones considered viable for RWH and those deemed unsuitable. Each factor was then rated based on the calculated weight, and the factor sub-levels were prioritized, as shown in Table 11.
The suitability scale denoted by 1 considers a restricted area for suitability analysis because of its insignificance for RWH site identification; a suitability scale of 7 denotes the most appropriate area for RWH. Before the suitability analysis, all thematic layers were converted to a similar spatial reference system, and a final suitability map for PRWH with four suitability categories was generated based on multi-criteria decision (MCD): high suitability, medium suitability, low suitability, and unsuitable categories. Figure 10 depicts the distribution of suitability for PRWH.
This study employed the AHP with a pairwise comparison matrix to establish weights for the factors influencing RWH suitability. These factors included rainfall, slope, soil texture, LULC, and drainage density. For each potential RWH location, specific criteria were evaluated for each factor. The rainfall amount, slope angle, soil texture type, LULC category, and drainage density values were identified. A score was assigned to each factor value based on predefined importance scales. For instance, rainfall between 501 and 619 mm received a score of 7 (highly suitable), while areas with rainfall between 294 and 400 mm received a score of 1 (unsuitable). Similar importance scales were established for other factors to determine suitability classes. The suitability score (S) for a location is calculated by multiplying the weight of each factor by its corresponding score and summing these values across all factors. Suitability classes were determined using two common approaches: Quantile Classification, which divides the range of S values into equal intervals (quartiles), and Natural Breaks Classification, which identifies significant changes in the S value distribution. Suitability classes are then assigned based on these breaks. Following this classification, suitability limits were established for each class (unsuitable, low suitability, medium suitability, and highly suitable), and a tabulated distribution of suitability values along with a statistical analysis was provided. The analysis results in Table 12 show that among the entire study area, 1% of the area was highly appropriate for RWH, covering approximately 3288 km2. Approximately 13% (37498 km2) of the area had medium suitability; most of the study area was represented by regions with low suitability, constituting 75% or 242,170 km2 of the area. A region with 13%, or 41,000 km2, was regarded as unsuitable for RWH. Finally, adopting sustainable water management strategies, such as integrated water resource management that considers groundwater and other resources, is essential to combat water scarcity.

4. Discussion

This study aimed to identify an appropriate site for harvesting rainwater for communities lacking rainfall, especially in semi-arid and arid areas. Several studies have confirmed the effectiveness of RWH systems in providing water for diverse purposes in drought-prone regions [102,103,104,105,106,107]. In the study area, low rainfall patterns are characterized by infrequent yet intense precipitation events, resulting in significant surface runoff and flooding. To mitigate water scarcity during extended dry seasons, the implementation of large-scale RWH structures is imperative. Despite the irregular nature of rainfall, data spanning from 1997 to 2015 indicate an adequate number of rainy days annually (up to 49), supporting the feasibility of RWH through strategic site selection for large-scale structures [108]. Identifying appropriateness using a suitability analysis for any plan often requires several possible options that can be assessed based on multiple criteria. Creating suitability index maps through GIS and MCDA models can be challenging due to the influence of numerous factors [109]. This study employed the WLC method to assign weights to various thematic layers within the model. This approach offers significant flexibility in weight allocation, allowing for a more nuanced consideration of these factors [110]. This study established the factors used to locate suitable areas for RWH based on their appropriateness and applicability to potentially suitable sites. Based on the leveling assignment, the weighted factors and their subgroup (internally ranked factor) intersections produced the PRWH site. The areas identified were technically in line with similar studies, and the results depict how the importance of weight determines the selection of suitable areas. Most sites identified as highly suitable for PRWH are situated in areas with clay soil, high rainfall, gentle slopes, shrublands, bare land, and high drainage density.
According to [111], clay soil has the highest runoff potential with a very low infiltration rate and is highly suitable for water harvesting, whereas sandy loam and sandy soils have the least potential for water storage because they have the lowest water-holding capacity due to their high infiltration rate. The same is true for sparsely vegetated land, such as bare land and shrubland, which are types of land cover that are more often recommended for water harvesting site selection [112]. The flat and gentle regions in the study area also had the highest share of the PRWH site. This aligns with findings from another study [113], steep slopes are unsuitable for rainwater storage; however, they significantly contribute to overflow generation and the speed of water flow to the harvesting area on flat and gentle slopes.
A total of 13% of the area was an excellent and reasonably proper zone for RWH; therefore, attention should be given to this area for the construction of RWH structures. These areas are best suited for rainfall, slope, land use, land cover, soil texture, and drainage density. However, in most of the areas studied, approximately 87% were unsuitable for RWH. By naming the suitability location by administrative unit, the following districts (smaller administrative units) and zones (higher administrative units) in the study area were identified. These were the Dugda Dawa, Yabelo, and Melka Soda districts and some areas in the Borena zone of Oromia. The Hudet District in the Liben zone; the Bare, Chereti, and West Imi districts in the Afder zone; and East Imi in the Shebelle zone in the Somali regional state were identified as suitable for RWH. The exact locations can be determined by extracting the geographic coordinates (X and Y) from the potential water harvesting (PRWH) suitability map.
The final PRWH suitability map, shown in Figure 9, reflects the suitability levels of the RWH site. For instance, the most suitable sites for RWH are located in areas with higher rainfall (>450 mm), with slopes extending from 0 to 10, and in large drainage density areas (8–13) because rainfall and drainage density had the highest weight factors, followed by slope. These results agree with the assumption of factor suitability level allocation. The results also reveal that areas of high and medium suitability for RWH were mostly located in areas where shrubland, bare land, clay, and clay loam were predominant. These LULC and soil types are appropriate for rainwater storage, as has been explained in several similar studies [114,115].
Generally, the analysis results show that significant conclusions were asserted from factors and methods used to develop a suitability result, including criteria levels of suitability and relative importance weights. All factors were combined in a weighted overly analysis using the weighted linear combination technique to produce suitable rainwater collection sites. According to a study by [116], RWH requires a significantly flat topography to harvest overflow water and easily construct the structure. When considering the study area, predominantly flat topography accounts for approximately 94% of the area (<5% slope), which makes it easy to harvest overflow water and reduces earthwork costs; however, rainfall variability and drainage density were the most important factors given the greater weights out of all factors, which were governed for the identification of a highly suitable location. Since the area has often experienced very low rainfall amounts, resulting in continued drought and the depletion of livestock in recent years, the implementation of RWH projects in the Borena zone in the Oromia and Somali regional states will be very important for crop production, water for livestock fodder, and even water for human consumption.

5. Conclusions

This study identified a suitable site for RWH using geospatial techniques and an MCDA approach for the Borena zone of the Oromia and Somali regional states. ArcGIS tools were the main instruments used to combine different thematic layers, and it was convenient to analyze the data using a weighted overly process with weighted linear combination (WLC) to identify suitable locations. The criteria/factor importance weights were extracted using the AHP method. The area is more often affected by drought due to low precipitation over a longer period and is usually called a drought-prone area. This area is more dominantly characterized by arid and semi-arid climatic conditions that attract attention to RWH options to mitigate drought problems as a result of low precipitation.
RWH is an important method for overcoming water deficiency during longer drought periods. This study highlights the potential of RWH for drought mitigation in drought-prone areas of the Somali and Borena zones. However, for long-term water security, an integrated approach is crucial. Combining RWH with existing groundwater resources and other water management strategies can maximize water availability. The result of this study leveraged multiple factors influencing RWH suitability, including elevation, slope, precipitation, land use, soil type, and drainage density. The AHP was employed to assess the relative importance of these factors. Through pairwise comparisons and weight assignment, we established a robust method for evaluating potential RWH sites. This approach yielded a suitability map with clearly defined criteria levels and weightings, providing valuable insights for informed RWH site selection. Drainage density and rainfall factors are given higher weights than other factors to increase their usefulness for RWH. The suitability model generated four suitability classes: high suitability, medium suitability, low suitability, and unsuitable. Eighteen suitable sites were identified in the study area, of which thirteen were highly suitable, and five were considered moderately suitable. Almost all highly suitable areas were surrounded by moderately suitable areas, meaning that these two suitable areas were adjacent to each other. This adjacency provides flexible opportunities for constructing RWH structures.
The findings of this study will guide regional state authorities, concerned federal state bodies, and decision makers in planning future development approaches for handling water insufficiency within the selected study area. Additionally, this study guides the design of cost-effective interventions aimed at maximizing rainwater utilization and alleviating water scarcity in the Ethiopian regions under study. The technique used to identify appropriate RWH sites can also be adapted for other low-precipitation areas. However, lack of site visits, ground validation, and complete socio-economic data were some of the limitations in this study, so additional investigation is needed to execute the RWH system, such as comprehensive site depictions, social and economic viability characterization, in-depth assessments of environmental sustainability, and detailed cost–benefit analyses.

6. Future Work

Further studies should be conducted focusing on refining the decision-making framework by incorporating additional criteria or datasets to improve the accuracy of RWH site selection. This could involve integrating socioeconomic factors, such as population density or land tenure arrangements, to better capture the sociocultural context of the study area. Additionally, exploring the integration of climate change scenarios into the MCDA framework may improve the robustness of RWH interventions. Moreover, exploring advanced modeling techniques, such as machine learning algorithms or spatial optimization methods, could enhance the predictive capabilities of the RWH suitability analysis and facilitate more informed decision making. Additionally, longitudinal studies tracking the performance of selected RWH sites over time could provide valuable insights into their effectiveness and inform adaptive management strategies for sustainable water resource management in drought-prone regions.

Author Contributions

A.K.B. contributed to writing, the development of remote sensing models, dataset preprocessing, and correction. D.J.B. contributed to writing, data collection, methodology, and the development of remote sensing models. A.M.F.A.-Q. contributed to editing, enhancing the entire text, and providing supervision. T.G.D. provided supervision, contributed to the methodology, and participated in reviewing. S.N. provided support in structuring the manuscript, offering guidance, and conducting quality checks. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All data supporting the findings of this study are available in the manuscript.

Acknowledgments

The authors extend their gratitude to the Ethiopian Artificial Intelligence Institute (EAI) Research and Development Sector, the Ethiopian Space Science and Geospatial Institute (ESSGI), and the Ministry of Water and Energy (MoE) for generously providing the dataset that was essential for this study. Their support was instrumental in the successful completion of this research.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Food and Agricultural Organization (FAO). Strengthening Agricultural Water Efficiency and Productivity on the African and Global Level: Status, Performance, and Scope Assessment of Water Harvesting in Uganda, Burkina Faso, and Morocco; FAO: Rome, Italy, 2016. [Google Scholar]
  2. Al-Quraishi, A.M.; Gaznayee, H.A.; Crespi, M. Drought trend analysis in a semi-arid area of Iraq based on normalized difference vegetation index, normalized difference water index, and standardized precipitation index. J. Arid. Land 2021, 13, 413–430. [Google Scholar] [CrossRef]
  3. Cassardo, C.; Jones, J.A.A. Managing water in a changing world. Water 2011, 3, 618–628. [Google Scholar] [CrossRef]
  4. Sayl, K.N.; Mohammed, A.S.; Ahmed, A.D. GIS-based approach for rainwater harvesting site selection. IOP Conf. Ser. Mater. Sci. Eng. 2020, 737, 012246. [Google Scholar] [CrossRef]
  5. Tiwari, K.; Goyal, R.; Sarkar, A. GIS-based methodology for identification of suitable locations for rainwater harvesting structures. Water Resour. Manag. 2018, 32, 1811–1825. [Google Scholar] [CrossRef]
  6. Salem, H.S.; Pudza, M.Y.; Yihdego, Y. Water strategies and water–food Nexus: Challenges and opportunities towards sustainable development in various regions of the World. Sustain. Water Resour. Manag. 2022, 8, 114. [Google Scholar] [CrossRef] [PubMed]
  7. Matomela, N.; Li, T.; Ikhumhen, H.O. Siting of rainwater harvesting potential sites in arid or semi-arid watersheds using GIS-based techniques. Environ. Process. 2020, 7, 631–652. [Google Scholar] [CrossRef]
  8. Waghaye, A.M.; Singh, D.K.; Sarangi, A.; Sena, D.R.; Sahoo, R.N.; Sarkar, S.K. Identification of suitable zones and sites for rainwater harvesting using GIS and multicriteria decision analysis. Environ. Monit. Assess. 2023, 195, 279. [Google Scholar] [CrossRef] [PubMed]
  9. Gaznayee, H.A.A.; Al-Quraishi, A.M.F.; Mahdi, K.; Ritsema, C. A geospatial approach for analysis of drought impacts on vegetation cover and land surface temperature in the Kurdistan Region of Iraq. Water 2022, 14, 927. [Google Scholar] [CrossRef]
  10. Almamalachy, Y.S.; Al-Quraishi, A.M.F.; Moradkhani, H. Agricultural drought monitoring over Iraq utilizing MODIS products. In Environmental Remote Sensing and GIS in Iraq; Springer: Berlin/Heidelberg, Germany, 2020; pp. 253–278. [Google Scholar] [CrossRef]
  11. Kumar, M.G.; Agarwal, A.K.; Bali, R. Delineation of potential sites for water harvesting structures using remote sensing and GIS. J. Indian Soc. Remote Sens. 2008, 36, 323–334. [Google Scholar] [CrossRef]
  12. Santosh, N.; Vinay, G. Rainwater harvesting for domestic use: A systematic review and outlook from the utility policy and management perspectives. Util. Policy 2022, 77, 101383. [Google Scholar] [CrossRef]
  13. Gwenzi, W.; Nyamadzawo, G. Hydrological impacts of urbanization and urban roof water harvesting in water-limited catchments: A review. Environ. Process. 2014, 1, 573–593. [Google Scholar] [CrossRef]
  14. World Bank Group. Climate Change Knowledge Portal. For Development Practitioners and Policymakers. Available online: https://climateknowledgeportal.worldbank.org/ (accessed on 11 March 2021).
  15. Kebede, S.; Hailu, A.; Crane, E.; Ó Dochartaigh, B.É.; Bellwood-Howard, I. Africa Groundwater Atlas: Hydrogeology of Ethiopia. Br. Geol. Surv. 2018, 8, 36. [Google Scholar] [CrossRef]
  16. Alemayehu, Z.Y.; Minale, A.S.; Legesse, S.A. Spatiotemporal rainfall and temperature variability in Suha Watershed, Upper Blue Nile Basin, Northwest Ethiopia. Environ. Monit. Assess. 2022, 194, 538. [Google Scholar] [CrossRef] [PubMed]
  17. Alene, A.; Yibeltal, M.; Abera, A.; Andualem, T.G.; Lee, S.S. Identifying rainwater harvesting sites using integrated GIS and a multi-criteria evaluation approach in semi-arid areas of Ethiopia. Appl. Water Sci. 2022, 12, 238. [Google Scholar] [CrossRef]
  18. World Meteorological Organization (WMO). Weather. Climate. Water. Available online: https://public.wmo.int/en/our-Mandate/climate/wmo-statement-state-of-global-climate (accessed on 21 April 2023).
  19. Food and Agricultural Organization (FAO). Rapid Drought Assessment Report for Borena and Dawa Zones 2021; FAO: Rome, Italy, 2021. [Google Scholar]
  20. Ding, W.; Wang, F.; Dong, Y.; Jin, K.; Cong, C.; Han, J.; Ge, W. Effects of rainwater harvesting system on soil moisture in rain-fed orchards on the Chinese Loess Plateau. Agric. Water Manag. 2021, 243, 106496. [Google Scholar] [CrossRef]
  21. Oweis, T.; Hachum, A. Water harvesting and supplemental irrigation for improved water productivity of dry farming systems in West Asia and North Africa. Agric. Water Manag. 2006, 80, 57–73. [Google Scholar] [CrossRef]
  22. Gebremichael, S.G.; Yismaw, E.; Tsegaw, B.D.; Shibeshi, A.D. Determinants of water source use, quality of water, sanitation and hygiene perceptions among urban households in North-West Ethiopia: A cross-sectional study. PLoS ONE 2021, 16, e0239502. [Google Scholar] [CrossRef] [PubMed]
  23. Worku, M.A.; Feyisa, G.L.; Beketie, K.T.; Garbolino, E. Rainfall variability and trends in the Borana Zone of Southern Ethiopia. J. Water Clim. Chang. 2022, 13, 3132–3151. [Google Scholar] [CrossRef]
  24. Hatibu, N.; Mutabazi, K.; Senkondo, E.M.; Msangi, A.S.K. Economics of rainwater harvesting for crop enterprises in semi-arid areas of East Africa. Agric. Water Manag. 2006, 80, 74–86. [Google Scholar] [CrossRef]
  25. Zhao, L.; Hou, R.; Wu, F. Rainwater harvesting capacity of soils subjected to reservoir tillage during rainfall on the Loess Plateau of China. Agric. Water Manag. 2019, 217, 193–200. [Google Scholar] [CrossRef]
  26. Tolossa, T.T.; Abebe, F.B.; Girma, A.A. Rainwater harvesting technology practices and implication of climate change characteristics in Eastern Ethiopia. Cogent Food Agric. 2020, 6, 1724354. [Google Scholar] [CrossRef]
  27. Liersch, S.; Fournet, S.; Koch, H.; Djibo, A.G.; Reinhardt, J.; Kortlandt, J.; Van Weert, F.; Seidou, O.; Klop, E.; Baker, C.; et al. Water resources planning in the Upper Niger River basin: Are there gaps between water demand and supply? J. Hydrol. Reg. Stud. 2019, 21, 176–194. [Google Scholar] [CrossRef]
  28. Ahmed, S.; Jesson, M.; Sharifi, S. Selection frameworks for potential rainwater harvesting sites in arid and semi-arid regions: A systematic literature review. Water 2023, 15, 2782. [Google Scholar] [CrossRef]
  29. Bera, A.; Mukhopadhyay, B.P. Identification of suitable sites for surface rainwater harvesting in the drought-prone Kumari River Basin, India in the context of irrigation water management. J. Hydrol. 2023, 621, 129655. [Google Scholar] [CrossRef]
  30. Owusu, S.; Asante, R. Rainwater harvesting and primary uses among rural communities in Ghana. J. Water Sanit. Hyg. Dev. 2020, 10, 502–511. [Google Scholar] [CrossRef]
  31. Raimondi, A.; Quinn, R.; Abhijith, G.R.; Becciu, G.; Ostfeld, A. Rainwater harvesting and treatment: State of the art and perspectives. Water 2023, 15, 1518. [Google Scholar] [CrossRef]
  32. Ramakrishnan, D.; Durga Rao, K.H.V.; Tiwari, K.C. Delineation of potential sites for water harvesting structures through remote sensing and GIS techniques: A case study of Kali watershed, Gujarat, India. Geocarto Int. 2008, 23, 95–108. [Google Scholar] [CrossRef]
  33. Araya, A.; Stroosnijder, L. Assessing drought risk and irrigation need in northern Ethiopia. Agric. For. Meteorol. 2011, 151, 425–436. [Google Scholar] [CrossRef]
  34. Nyamwanza, A.M.; Kujinga, K.K. Climate change, sustainable water management and institutional adaptation in rural sub-Saharan Africa. Environ. Dev. Sustain. 2017, 19, 693–706. [Google Scholar] [CrossRef]
  35. Mahmoud, S.H.; Alazba, A.A.; Adamowski, J.; El-Gindy, A.M. GIS methods for sustainable stormwater harvesting and storage using remote sensing for land cover data-location assessment. Environ. Monit. Assess. 2015, 187, 598. [Google Scholar] [CrossRef]
  36. Jha, M.K.; Chowdary, V.M.; Kulkarni, Y.; Mal, B.C. Rainwater harvesting planning using geospatial techniques and multicriteria decision analysis. Resour. Conserv. Recycl. 2014, 83, 96–111. [Google Scholar] [CrossRef]
  37. Desta, S.; Block, P. Application of multi-criteria analysis in prioritizing rainwater harvesting technologies for development interventions: A case study in the Upper Blue Nile Basin, Ethiopia. Agric. Water Manag. 2013, 123, 45–54. [Google Scholar] [CrossRef]
  38. Modak, S.; Das, D. Delineation of Suitable Zone for Rainwater Harvesting in Upper Catchment of Kumari River Basin, West Bengal and Jharkhand, India: Using AHP and Geospatial Techniques. Paideuma 2022, 14, 20–43. [Google Scholar]
  39. O’Donovan, J.; Thompson, A.; Stiles, C.; Opintan, J.A.; Kabali, K.; Willis, I.; Mutreja, A. Participatory approaches, local stakeholders, and cultural relevance facilitate an impactful community-based project in Uganda. Health Promot. Int. 2020, 35, 1353–1368. [Google Scholar] [CrossRef] [PubMed]
  40. Ramasamy, S.M.; Anbazhagan, S. Criteria and techniques of detecting site-specific mechanisms for artificial recharge case study from Ayyar basin, India. J.-Geol. Soc. India 1997, 50, 449–456. [Google Scholar]
  41. Ngangom, B.; Das, A.; Lal, R.; Idapuganti, R.G.; Layek, J.; Basavaraj, S.; Babu, S.; Yadav, G.S.; Ghosh, P.K. Double mulching improves soil properties and productivity of maize-based cropping systems in the eastern Indian Himalayas. Int. Soil Water Conserv. Res. 2020, 8, 308–320. [Google Scholar] [CrossRef]
  42. Shadeed, S. Developing a GIS-based suitability map for rainwater harvesting in the West Bank, Palestine. In Proceedings of the International Conference on Environmental Education for Sustainable Development: Plugging the Hole, Birzeit, Palestine, 16–17 November 2011; Birzeit University: Birzeit, Palestine, 2011; p. 13. [Google Scholar]
  43. Rejani, R.; Rao, K.V.; Srinivasa Rao, C.H.; Osman, M.; Sammi Reddy, K.; George, B.; Pratyusha Kranthi, G.S.; Chary, G.R.; Swamy, M.V.; Rao, P.J. Identification of potential rainwater-harvesting sites for the sustainable management of a semi-arid watershed. Irrig. Drain. 2017, 66, 227–237. [Google Scholar] [CrossRef]
  44. Pauw, E.; Oweis, T.; Youssef, J. Integrating Expert Knowledge in GIS to Locate Biophysical Potential for Water Harvesting: Methodology and a Case Study for Syria. Aleppo, Syrian Arab Republic: International Center for Agricultural Research in the Dry Areas (ICARDA). 2008. Available online: https://hdl.handle.net/20.500.11766/8855 (accessed on 10 May 2024).
  45. Adham, A.; Riksen, M.; Ouessar, M.; Ritsema, C.J. A Methodology to assess and evaluate rainwater harvesting techniques in (Semi-) arid regions. Water 2016, 8, 198. [Google Scholar] [CrossRef]
  46. Ministry of Agriculture (MoA). Agro-Ecological Zone of Ethiopia, Natural Resources Management, and Regulatory Department with the support of the German Agency for Technical Cooperation (GTZ); Ministry of Agriculture: Addis Ababa, Ethiopia, 1998.
  47. Bogale, G.A.; Erena, Z.B. Drought vulnerability and impacts of climate change on livestock production and productivity in Different Agro-Ecological Zones of Ethiopia. J. Appl. Anim. Res. 2022, 50, 471–489. [Google Scholar] [CrossRef]
  48. NASA. ASTER Digital Elevation Model AST14DEM Relative. Available online: https://asterweb.jpl.nasa.gov/content/03_data/01_Data_Products/release_DEM_relative.htm (accessed on 10 May 2024).
  49. Das, S.; Patel, P.P.; Sengupta, S. Evaluation of different digital elevation models for analyzing drainage morphometric parameters in mountainous terrain: A case study of the Supin–Upper Tons Basin, Indian Himalayas. SpringerPlus 2016, 5, 1544. [Google Scholar] [CrossRef]
  50. RCMRD Geo-Portal. Ethiopia Sentinel2 Land Use Land Cover 2016-GeoNode. Available online: https://opendata.rcmrd.org/datasets/rcmrd::ethiopia-sentinel2-lulc2016/about (accessed on 10 May 2024).
  51. Belete, E.; Ayalew, A.; Ahmed, S. Associations of biophysical factors with Faba Bean Root Rot (Fusarium solani) epidemics in the Northeastern Highlands of Ethiopia. Crop Prot. 2013, 52, 39–46. [Google Scholar] [CrossRef]
  52. Assefa, W.W.; Abebe, W.B. GIS modeling of potentially suitable sites for aquaculture development in the Lake Tana Basin, Northwest Ethiopia. Agric. Food Secur. 2018, 7, 72. [Google Scholar] [CrossRef]
  53. Wable, P.S.; Jha, M.K.; Chowdary, V.M.; Mahapatra, S. GIS-based multi-criteria decision analysis for identifying rainwater harvesting structures sites in a semiarid river basin. In Impacts of Urbanization on Hydrological Systems in India; Springer International Publishing: Cham, Switzerland, 2023; pp. 1–23. [Google Scholar] [CrossRef]
  54. Islam, M.R.; Aziz, M.T.; Alauddin, M.; Kader, Z.; Islam, M.R. Site suitability assessment for solar power plants in Bangladesh: A GIS-based Analytical Hierarchy Process (AHP) and Multi-Criteria Decision Analysis (MCDA) approach. Renew. Energy 2024, 220, 119595. [Google Scholar] [CrossRef]
  55. Hameed, H. Water Harvesting in Erbil Governorate, Kurdistan Region, Iraq: Detection of Suitable Sites Using Geographic Information System and Remote Sensing. Master’s Thesis, Lund University, Lund, Sweden, 2013. Student Thesis Series INES. [Google Scholar]
  56. Islam, M.R.; Ahmad, M.R.; Alam, M.N.; Marma, M.S.; Gafur, A. Micro-watershed delineation and potential site selection for runoff water harvesting using remote sensing and GIS in a Hilly Area of Bangladesh. Am. J. Water Resour. 2020, 8, 134–144. [Google Scholar]
  57. Mahmoud, S.H.; Alazba, A.A. Integrated Remote sensing and GIS-based approach for deciphering groundwater potential zones in the Central Region of Saudi Arabia. Environ. Earth Sci. 2016, 75, 344. [Google Scholar] [CrossRef]
  58. Food and Agricultural Organization (FAO). Water Harvesting. In A Manual for the Design and Construction of Water Harvesting Schemes for Plant Production; FAO: Rome, Italy, 1991. [Google Scholar]
  59. Mugo, G.M.; Odera, P.A. Site selection for rainwater harvesting structures in Kiambu County-Kenya. Egypt. J. Remote Sens. Space Sci. 2019, 22, 155–164. [Google Scholar] [CrossRef]
  60. Rajasekhar, M.; Gadhiraju, S.R.; Kadam, A.; Bhagat, V. Identification of groundwater recharge-based potential rainwater harvesting sites for sustainable development of a semiarid region of Southern India using geospatial, AHP, and SCS-CN approaches. Arab. J. Geosci. 2020, 13, 24. [Google Scholar] [CrossRef]
  61. Endalkachew, Y.S.; Ejegu, M.A.; Tolossa, A.T.; Teka, A.H.; Andualem, T.G.; Tegegne, M.A.; Walle, W.M.; Shibeshie, S.E.; Dirar, T.M. Geospatial and AHP approach rainwater harvesting site identification in drought-prone areas, South Gonder Zone, Northwest Ethiopia. J. Indian Soc. Remote Sens. 2022, 50, 1321–1331. [Google Scholar] [CrossRef]
  62. Yang, R.; Xing, B. A comparison of the performance of different interpolation methods in replicating rainfall magnitudes under different climatic conditions in Chongqing Province (China). Atmosphere 2021, 12, 1318. [Google Scholar] [CrossRef]
  63. Adham, A.; Sayl, K.N.; Abed, R.; Abdeladhim, M.A.; Wesseling, J.G.; Riksen, M.; Fleskens, L.; Karim, U.; Ritsema, C.J. A GIS-based approach for identifying potential sites for harvesting rainwater in the western desert of Iraq. Int. Soil Water Conserv. Res. 2018, 6, 297–304. [Google Scholar] [CrossRef]
  64. Chowdhury, M.; Paul, P.K. Identification of suitable sites for rainwater harvesting using fuzzy AHP and fuzzy gamma operator: A case study. Arab. J. Geosci. 2021, 14, 585. [Google Scholar] [CrossRef]
  65. Soil Grids 250m v2.0 API. Available online: https://gee-community-catalog.org/projects/isric/ (accessed on 9 March 2023).
  66. Al-Ali, M. Soil Water Conservation and Water Balance Model for Micro-Catchment Water Harvesting System. Doctoral Thesis, Loughborough University, Loughborough, UK, 2012. Available online: https://hdl.handle.net/2134/10941 (accessed on 16 November 2012).
  67. Ibrahim Faqe, G.R.; Rasul, A.; Hamid, A.A.; Ali, Z.F.; Dewana, A.A. Suitable site selection for rainwater harvesting and storage case Study Using Dohuk Governorate. Water 2019, 11, 864. [Google Scholar] [CrossRef]
  68. Rana, V.K.; Suryanarayana, T.M.V. GIS-Based multi-criteria decision-making method to identify potential runoff storage zones within the watershed. Ann. GIS 2020, 26, 149–168. [Google Scholar] [CrossRef]
  69. Camera, C.B.; Zomeni, Z.; Noller, J.S.; Zissimos, A.M.; Christoforou, I.C.; Bruggeman, A. A high-resolution map of soil types and physical properties for Cyprus: A digital soil mapping optimization. Geoderma 2017, 285, 35–49. [Google Scholar] [CrossRef]
  70. Al Marsumi, K.J.; Al Shamma, A.M. Selection of suitable sites for water harvesting structures in a flood-prone area using remote sensing and GIS–case study. J. Environ. Earth Sci. 2017, 7, 91–100. [Google Scholar] [CrossRef]
  71. Zhang, Z.; Zhou, Y.; Liu, H.; Gao, H. In-Situ water level measurement using a NIR-imaging video camera. Flow Meas. Instrum. 2019, 67, 95–106. [Google Scholar] [CrossRef]
  72. US Department of Agriculture, Soil Conservation Service (USDA-SCS). National Engineering Handbook; USDA: Washington, DC, USA, 1975; Section 4.
  73. Syarifuddin, M.; Oktaviantoro, D. A GIS-based approach to determine the priority area for rainwater harvest in Kupang. In IOP Conference Series: Earth and Environmental Science; IOP Publishing: 2024; Volume 1311, Issue 1, p. 012032. Available online: https://iopscience.iop.org/article/10.1088/1755-1315/1311/1/012032/meta (accessed on 20 April 2024).
  74. Ashraf, M.U.; Masood, M.; Iqbal, M.; Ahsan, M.F.; Mustafa, A.U. Ascertainment of the most suitable water harvesting sites in a data-scarce region via an integrated PROMETHEE-II–AHP method in GIS environment. Water Pract. Technol. 2024, 19, 796–811. [Google Scholar] [CrossRef]
  75. Shadmehri, A.; Ghasemi, E.; Ghassemi, A.; Cheshomi, A.; Alaghmand, S. A multi-criteria decision analysis approach towards efficient rainwater harvesting. J. Hydrol. 2020, 582, 124501. [Google Scholar] [CrossRef]
  76. Alaghmand, S.; Beecham, S.; Hassanli, A. Impacts of vegetation cover on surface-groundwater flows and solute interactions in a semi-arid saline floodplain: A case study of the Lower Murray River, Australia. Environ. Process. 2014, 1, 59–71. [Google Scholar] [CrossRef]
  77. Al-Ghobari, H.; Dewidar, A.Z. Integrating GIS-Based MCDA Techniques and the SCS-CN method for identifying potential zones for rainwater harvesting in a semi-arid area. Water 2021, 13, 704. [Google Scholar] [CrossRef]
  78. Gebremedhin, A.Y.; Getahun, Y.S.; Moges, A.S.; Tesfay, F. Identification of suitable rainwater harvesting sites using geospatial techniques with AHP in Chacha Watershed, Jemma Sub-Basin Upper Blue Nile, Ethiopia. Air Soil Water Res. 2023, 16. [Google Scholar] [CrossRef]
  79. Aziz, S.F.; Abdulrahman, K.Z.; Ali, S.S.; Karakouzian, M. Water harvesting in the Garmian region (Kurdistan, Iraq) Using GIS and remote sensing. Water 2023, 15, 507. [Google Scholar] [CrossRef]
  80. Gavhane, K.P.; Mishra, A.K.; Sarangi, A.; Singh, D.K.; Sudhishri, S. Targeting of rainwater harvesting structures using geospatial tools and Analytical Hierarchy Process (AHP) in the semi-arid region of Rajasthan (India). Environ. Sci. Pollut. Res. 2023, 30, 61682–61709. [Google Scholar] [CrossRef]
  81. Wijesinghe, D.C.; Mishra, P.K.; Withanage, N.C.; Abdelrahman, K.; Mishra, V.; Tripathi, S.; Fnais, M.S. Application of GIS, Multi-Criteria Decision-Making Techniques for Mapping Groundwater Potential Zones: A Case Study of Thalawa Division, Sri Lanka. Water 2023, 15, 3462. [Google Scholar] [CrossRef]
  82. Saaty, R.W. The analytic hierarchy process-what it is and how it is used. Math Model. 1987, 9, 161–176. [Google Scholar] [CrossRef]
  83. Saaty, T.L. The Analytic Hierarchy Process; McGraw-Hill: New York, NY, USA, 1980. [Google Scholar]
  84. Fagbohun, B.J. Integrating GIS and multi-influencing factor technique for delineation of potential groundwater recharge zones in parts of Ilesha Schist Belt, Southwestern Nigeria. Environ. Earth Sci. 2018, 77, 69. [Google Scholar] [CrossRef]
  85. Saha, A.; Patil, M.; Karwariya, S.; Pingale, S.M.; Azmi, S.; Goyal, V.C.; Rathore, D.S. Identification of potential sites for water harvesting structures using geospatial techniques and multi-criteria decision analysis. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2018, 42, 329–334. [Google Scholar] [CrossRef]
  86. Bojer, A.K.; Woldesilassie, F.F.; Debelee, T.G.; Kebede, S.R.; Esubalew, S.Z. AHP and machine learning-based military strategic site selection: A case study of Adea District East Shewa Zone, Ethiopia. J. Sens. 2023, 2023, 6651486. [Google Scholar] [CrossRef]
  87. Preeti, P.; Shendryk, Y.; Rahman, A. Identification of suitable sites using gis for rainwater harvesting structures to meet irrigation demand. Water 2022, 14, 3480. [Google Scholar] [CrossRef]
  88. Golub, G.H.; Van Loan, C.F. Matrix Computations; JHU Press: St, Baltimore, MD, USA, 2013. [Google Scholar]
  89. Belton, V.; Stewart, T. Multiple Criteria Decision Analysis: An Integrated Approach; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  90. Ramanathan, R. Data envelopment analysis for weight derivation and aggregation in the analytic hierarchy process. Comput. Oper. Res. 2006, 33, 1289–1307. [Google Scholar] [CrossRef]
  91. Jozaghi, A.; Alizadeh, B.; Hatami, M.; Flood, I.; Khorrami, M.; Khodaei, N.; Ghasemi Tousi, E. A comparative study of the AHP and TOPSIS techniques for dam site selection using GIS: A case study of Sistan and Baluchestan Province, Iran. Geosciences 2018, 8, 494. [Google Scholar] [CrossRef]
  92. Vaidya, O.S.; Kumar, S. Analytic hierarchy process: An overview of applications. Eur. J. Oper. Res. 2006, 169, 1–29. [Google Scholar] [CrossRef]
  93. Montibeller, G.; Patel, P.; del Rio Vilas, V.J. A critical analysis of multi-criteria models for the prioritization of health threats. Eur. J. Oper. Res. 2020, 281, 87–99. [Google Scholar] [CrossRef]
  94. Drobne, S.; Lisec, A. Multi-attribute decision analysis in GIS: Weighted linear combination and ordered weighted averaging. Informatica 2009, 33, 4. [Google Scholar]
  95. Singh, V.; Kumar, V.; Singh, V.B. A hybrid novel fuzzy AHP-Topsis technique for selecting parameter-influencing testing in software development. Decis. Anal. J. 2023, 6, 100159. [Google Scholar] [CrossRef]
  96. Taye, M.; Simane, B.; Zaitchik, B.F.; Selassie, Y.G.; Setegn, S. Rainfall variability across the agro-climatic zones of a tropical highland: The case of the Jema watershed, northwestern Ethiopia. Environments 2019, 6, 118. [Google Scholar] [CrossRef]
  97. Yona, Y.; Matewos, T.; Sime, G. Analysis of rainfall and temperature variabilities in Sidama regional state, Ethiopia. Heliyon 2024, 10, e28184. [Google Scholar] [CrossRef]
  98. Karakuş, C.B.; Yıldız, S. GIS-multi criteria decision analysis-based land suitability assessment for dam site selection. Int. J. Environ. Sci. Technol. 2022, 19, 12561–12580. [Google Scholar] [CrossRef]
  99. Gashaw, D.Y.; Suryabhagavan, K.V.; Nedaw, D.; Gummadi, S. Rainwater harvesting in Modjo Watershed, Upper Awash River Basin, Ethiopia through remote sensing and Fuzzy AHP. Geocarto Int. 2022, 37, 14785–14810. [Google Scholar] [CrossRef]
  100. Sarkar, S.; Biswas, S. Application of integrated AHP-Entropy model in suitable site selection for rainwater harvesting structures: A case study of Upper Kangsabati Basin, India. Arab. J. Geosci. 2022, 15, 1684. [Google Scholar] [CrossRef]
  101. De Winnaar, G.; Jewitt, G.P.W.; Horan, M. A GIS-Based approach for identifying potential runoff harvesting sites in the Thukela River Basin, South Africa. Phys. Chem. Earth Parts A/B/C 2007, 32, 1058–1067. [Google Scholar] [CrossRef]
  102. Hashim, H.Q.; Sayl, K.N. Detection of Suitable Sites for Rainwater Harvesting Planning in an Arid Region Using Geographic Information System. Appl. Geomat. 2021, 13, 235–248. [Google Scholar] [CrossRef]
  103. Kapli, F.W.A.; Azis, F.A.; Suhaimi, H.; Shamsuddin, N.; Abas, P.E. Feasibility studies of rainwater harvesting system for ablution purposes. Water 2023, 15, 1686. [Google Scholar] [CrossRef]
  104. Rahman, A.; Keane, J.; Imteaz, M.A. Rainwater harvesting in Greater Sydney: Water savings, reliability and economic benefits. Resour. Conserv. Recycl. 2012, 61, 16–21. [Google Scholar] [CrossRef]
  105. Abas, P.E.; Mahlia, T.M.I. Techno-economic and sensitivity analysis of rainwater harvesting system as alternative water source. Sustainability 2019, 11, 2365. [Google Scholar] [CrossRef]
  106. Staddon, C.; Rogers, J.; Warriner, C.; Ward, S.; Powell, W. Why Doesn’t every family practice rainwater harvest? Factors that affect the decision to adopt rainwater harvesting as a household water security strategy in Central Uganda. Water Int. 2018, 43, 1114–1135. [Google Scholar] [CrossRef]
  107. Al-Batsh, N.; Al-Khatib, I.A.; Ghannam, S.; Anayah, F.; Jodeh, S.; Hanbali, G.; van der Valk, M. Assessment of rainwater harvesting systems in poor rural communities: A case study from Yatta Area, Palestine. Water 2019, 11, 585. [Google Scholar] [CrossRef]
  108. Tekalign, S. Rainfall trends and variability in selected areas of Ethiopian Somali Regional State, Southeastern Ethiopia. Ethiop. J. Environ. Stud. Manag. 2011, 10, 162–175. [Google Scholar] [CrossRef]
  109. Aklan, M.; Al-Komaim, M.; de Fraiture, C. Site suitability analysis of indigenous rainwater harvesting systems in arid and data-poor environments: A case study of Sana’a Basin, Yemen. Environ. Dev. Sustain. 2023, 25, 8319–8342. [Google Scholar] [CrossRef]
  110. Ahmad, B.A.; Salar, S.G.; Shareef, A.J. An integrated new approach for optimizing rainwater harvesting system with dams site selection in the Dewana Watershed, Kurdistan Region, Iraq. Heliyon 2024, 10, e27273. [Google Scholar] [CrossRef]
  111. Al-Hasani, B.; Abdellatif, M.; Carnacina, I.; Harris, C.; Al-Quraishi, A.; Maaroof, B.F.; Zubaidi, S.L. Integrated geospatial approach for adaptive rainwater harvesting site selection under the impact of climate change. Stoch. Environ. Res. Risk Assess. 2024, 38, 1009–1033. [Google Scholar] [CrossRef]
  112. Bahrani, S.; Ebadi, T.; Ehsani, H.; Yousefi, H.; Maknoon, R. Modeling landfill site selection by multi-criteria decision making and fuzzy functions in GIS, Case Study: Shabestar, Iran. Environ. Earth Sci. 2016, 75, 337. [Google Scholar] [CrossRef]
  113. Taherdoost, H.; Madanchian, M. Multi-Criteria Decision-Making (MCDM) methods and concepts. Encyclopedia 2023, 3, 77–87. [Google Scholar] [CrossRef]
  114. Tahvili, Z.; Khosravi, H.; Malekian, A.; Khalighi Sigaroodi, S.; Pishyar, S.; Singh, V.P.; Ghodsi, M. Locating Suitable Sites for Rainwater Harvesting (RWH) in the Central Arid Region of Iran. Sustain. Water Resour. Manag. 2021, 7, 1–11. [Google Scholar] [CrossRef]
  115. Maina, C.W.; Raude, J.M. Assessing land suitability for rainwater harvesting using geospatial techniques: A case study of Njoro Catchment, Kenya. Appl. Environ. Soil Sci. 2016, 2016, 4676435. [Google Scholar] [CrossRef]
  116. Gao, Z.; Zhang, Q.; Wang, Y.; Jv, X.; Dzakpasu, M.; Wang, X.C. Evolution of water quality in rainwater harvesting systems during long-term storage in non-rainy seasons. Sci. Total Environ. 2024, 912, 168784. [Google Scholar] [CrossRef]
Figure 1. Location map of the study area.
Figure 1. Location map of the study area.
Water 16 01789 g001
Figure 2. A conceptual model for the selection of potential rainwater harvesting (PRWH) sites.
Figure 2. A conceptual model for the selection of potential rainwater harvesting (PRWH) sites.
Water 16 01789 g002
Figure 3. A map of meteorological stations in the study area.
Figure 3. A map of meteorological stations in the study area.
Water 16 01789 g003
Figure 4. A DEM map of the study area.
Figure 4. A DEM map of the study area.
Water 16 01789 g004
Figure 5. Average annual rainfall (2012–2021) distribution and reclassified rainfall map.
Figure 5. Average annual rainfall (2012–2021) distribution and reclassified rainfall map.
Water 16 01789 g005
Figure 6. Slope distribution and suitability class map.
Figure 6. Slope distribution and suitability class map.
Water 16 01789 g006
Figure 7. Soil texture class and reclassification suitability map.
Figure 7. Soil texture class and reclassification suitability map.
Water 16 01789 g007
Figure 8. LULC distribution and suitability map.
Figure 8. LULC distribution and suitability map.
Water 16 01789 g008
Figure 9. Drainage density distribution and drainage suitability class map.
Figure 9. Drainage density distribution and drainage suitability class map.
Water 16 01789 g009
Figure 10. Potential rainwater harvesting (PRWH) suitability map.
Figure 10. Potential rainwater harvesting (PRWH) suitability map.
Water 16 01789 g010
Table 1. Hydrologic soil group (HSG).
Table 1. Hydrologic soil group (HSG).
Soil GroupDescriptionTexture Class
ALowest overflow possibilities because of high penetration rateSand and loamy sand
BModerately low overflow possibilities and moderate penetration rateSandy loam and loam
CModerately high overflow possibilities due to slow infiltration rateSilt loam, silt, and sandy clay loam
DHighest overflow possibilities with very low infiltration rateClay loam, silt clay, sandy clay, and clay
Table 2. Foundational scale used for assessing relative significance of criteria in Analytic Hierarchy Process (AHP).
Table 2. Foundational scale used for assessing relative significance of criteria in Analytic Hierarchy Process (AHP).
Relative ImportanceIntensitiesDescription
1Equally importantBoth activities contribute equally to the objective.
3Moderately importantOne activity is slightly preferred over another.
5Strongly importantOne activity is significantly preferred over another.
7Very strongly importantOne activity is exceptionally preferred, dominating in practice.
9Extremely importantThe evidence supporting one activity surpasses all others with the highest level of confirmation.
2, 4, 6, 8Values between adjacent judgmentsAdditional subdivision or compromise when required.
Table 3. Random Index values.
Table 3. Random Index values.
Order12345678910
RI000.520.891.111.251.351.41.451.49
Table 4. AHP pair-wise matrix for factors.
Table 4. AHP pair-wise matrix for factors.
Drainage DensityLand Use/Land CoverSoil TextureSlopeMAAR
Drainage density17531
Land use/land cover1/711/31/51/7
Soil texture1/5311/21/5
Slope1/35311/3
Monthly average annual rainfall (MAAR)17531
Table 5. Weights of each criterion for RWH calculated by AHP model.
Table 5. Weights of each criterion for RWH calculated by AHP model.
CriteriaDrainage DensityLand Use/Land CoverSoil TextureSlopeMAARAverage (Weight)
Drainage Density0.37370.77780.6250.64290.37370.5586
Land Use/Land Cover0.00540.11110.04170.10710.00540.0541
Soil Texture0.07470.33330.1250.21430.07470.1644
Slope0.12460.55560.3750.42860.12460.3217
MAAR0.37370.77780.6250.64290.37370.5586
Table 6. Percentage of influences CI, RI, and CR for each criterion.
Table 6. Percentage of influences CI, RI, and CR for each criterion.
No.CriteriaWeight% of Influences
1Drainage density0.344234.42
2Land use/land cover0.03683.68
3Soil texture0.09929.92
4Slope0.175517.55
5MAAR0.344234.42
SUM1100
Note: Eigenvalue (λ) = 5.09797, CI = 0.024, RI = 1.11, CR = CI/RI, CR = 0.0216.
Table 7. Slope percentage categories, suitability level, and area coverage.
Table 7. Slope percentage categories, suitability level, and area coverage.
FactorSlope RangeSuitability LevelCoverage (%)Area (Km2)
Slope>15Not suitable0.411328.22
10−15Low suitability1.063433.92
5−10Medium suitability4.7515,387.86
<5High suitability93.78303,805
Total 100323,955
Table 8. Soil type suitability level and coverage.
Table 8. Soil type suitability level and coverage.
Factor TypeSuitability ClassCoverage (%)Coverage (Km2)
Soil typeClayHigh suitability10.834,987
Clay loamMedium suitability0.973142
Loam and loamy sandLow suitability29.6596,053
Sandy loam and sandNot suitable58.58189,773
Total 100323,955
Table 9. LULC, suitability class, and area coverage.
Table 9. LULC, suitability class, and area coverage.
FactorTypeSuitability ClassCoverage (%)Coverage (Km2)
Land use/coverShrublandHigh suitability44.12142,929
Bare landMedium suitability2.16803
Grassland and woodlandLow suitability37.36121,030
Cropland, forest, and water bodyNot suitable16.4253,193
Total 100323,955
Table 10. Drainage density intervals, suitability class, and coverage.
Table 10. Drainage density intervals, suitability class, and coverage.
FactorIntervalSuitability ClassCoverage (%)Coverage (km2)
Drainage density0–2Not suitable46.21149,700
3–5Low suitability39126,342
5–7Medium suitability12.9742,017
8–13High suitability1.815864
Total 100323,923
Table 11. Importance scale for suitability classes by factor.
Table 11. Importance scale for suitability classes by factor.
Importance Scale7531
Rainfall501–619451–500401–450294–400
Slope<55–1010–15>15
Soil textureClayClay loamLoam and loamy sandSandy loamy and sand
Land use/coverShrub landBare landGrassland and woodlandCropland, forest, and water body
Drainage density8–135–73–40–2
Table 12. Suitability levels and coverage area for PRWH site.
Table 12. Suitability levels and coverage area for PRWH site.
Suitability LevelCoverage (%)Coverage (km2)
Not suitable12.6641,000
Low suitability74.75242,170
Medium suitability11.5837,498
High suitability1%3288
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bojer, A.K.; Bekalo, D.J.; Debelee, T.G.; Nadarajah, S.; Al-Quraishi, A.M.F. Rainwater Harvesting Site Selection for Drought-Prone Areas in Somali and Borena Zones, Oromia Regional State, Ethiopia: A Geospatial and Multi-Criteria Decision Analysis. Water 2024, 16, 1789. https://doi.org/10.3390/w16131789

AMA Style

Bojer AK, Bekalo DJ, Debelee TG, Nadarajah S, Al-Quraishi AMF. Rainwater Harvesting Site Selection for Drought-Prone Areas in Somali and Borena Zones, Oromia Regional State, Ethiopia: A Geospatial and Multi-Criteria Decision Analysis. Water. 2024; 16(13):1789. https://doi.org/10.3390/w16131789

Chicago/Turabian Style

Bojer, Amanuel Kumsa, Desta Jula Bekalo, Taye Girma Debelee, Saralees Nadarajah, and Ayad M. Fadhil Al-Quraishi. 2024. "Rainwater Harvesting Site Selection for Drought-Prone Areas in Somali and Borena Zones, Oromia Regional State, Ethiopia: A Geospatial and Multi-Criteria Decision Analysis" Water 16, no. 13: 1789. https://doi.org/10.3390/w16131789

APA Style

Bojer, A. K., Bekalo, D. J., Debelee, T. G., Nadarajah, S., & Al-Quraishi, A. M. F. (2024). Rainwater Harvesting Site Selection for Drought-Prone Areas in Somali and Borena Zones, Oromia Regional State, Ethiopia: A Geospatial and Multi-Criteria Decision Analysis. Water, 16(13), 1789. https://doi.org/10.3390/w16131789

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

Article Metrics

Back to TopTop