Next Article in Journal
Influence of Rainfall Pattern on Wetness Index for Infinite Slope Stability Analysis
Next Article in Special Issue
Managed Aquifer Recharge: A Key to Sustainability
Previous Article in Journal
Watershed Hydrological Responses to Land Cover Changes at Muger Watershed, Upper Blue Nile River Basin, Ethiopia
Previous Article in Special Issue
Integration of Managed Aquifer Recharge into the Water Supply System in the Algarve Region, Portugal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Fuzzy AHP-TOPSIS Model for Assessing Managed Aquifer Recharge Potential in a Hot Dry Region: A Case Study of Djibouti at a Country Scale

by
Rachid Mohamed Mouhoumed
1,2,*,
Ömer Ekmekcioğlu
3,
Eyyup Ensar Başakın
1 and
Mehmet Özger
1
1
Hydraulics Division, Civil Engineering Department, Istanbul Technical University, 34469 Istanbul, Turkey
2
Energy and Environment Research Center, Faculty of Engineering, University of Djibouti, Djibouti 1904, Djibouti
3
Disaster and Emergency Management Department, Disaster Management Institute, Istanbul Technical University, 34469 Istanbul, Turkey
*
Author to whom correspondence should be addressed.
Water 2023, 15(14), 2534; https://doi.org/10.3390/w15142534
Submission received: 5 June 2023 / Revised: 1 July 2023 / Accepted: 3 July 2023 / Published: 10 July 2023
(This article belongs to the Special Issue Managed Aquifer Recharge: A key to Sustainability)

Abstract

:
Given the prevailing arid climate and rapid population growth, groundwater resources face unprecedented challenges globally, including depletion, seawater intrusion, and contamination. Managed aquifer recharge (MAR) technologies have emerged as valuable solutions to address these pressing issues. However, identifying suitable regions for MAR activities is a complex task, particularly at the country level. Therefore, in this study, we propose a robust approach that combines the fuzzy analytical hierarchy process (AHP) and the technique for order of preference by similarity to ideal solution (TOPSIS) to delineate suitable sites for MAR structures. The proposed model was applied to Djibouti, a hot, dry, and water-stressed country. We identified a set of nine decision criteria and conducted a pairwise comparison survey to determine their relative importance. Additionally, the TOPSIS method was employed to integrate the decision layers and prioritize the study area. The results highlight the significance of rainfall, the slope, and the NDVI as the most influential decision parameters, while the drainage density has the least impact. A suitability analysis reveals that 16.38%, 17.96%, and 30.41% of the country have a very high, high, and moderate potential for MAR activities, respectively. Furthermore, a sensitivity analysis demonstrates the stability of the proposed model, affirming the usefulness of the generated suitability map.

1. Introduction

The escalating demand for groundwater due to population growth, coupled with a drying climate that limits recharge, is exacerbating the water situation and intensifying the deterioration of subsurface resources [1]. In the Horn of Africa (Djibouti, Somalia, Ethiopia, and Eritrea), water scarcity and climate variability have already manifested in many challenges threatening millions of humans and living species. For instance, Djibouti’s population heavily relies on groundwater resources, as permanent surface water sources, such as rivers and lakes, are nonexistent. However, the heavy abstraction of underground resources for several decades, coupled with highly variable rainfall in the region, has caused significant water table declines and seawater intrusion, forcing more rural dwellers and herders to seek refuge in urban areas after losing their livelihood sources [2]. The situation in Somalia is much worse, as recurrent droughts in the region caused by four consecutive failed rainy seasons have resulted in unsustainable means of groundwater extraction, including the creation of more wells and boreholes, further aggravating environmental degradation [3]. Consequently, about 6.1 million Somalis have been impacted by the consequences of these droughts, leading to massive internal population displacement [3].
The need to implement sustainable practices to mitigate the challenges faced by aquifer systems has become increasingly pressing. In this regard, managed aquifer recharge (MAR) has emerged as a promising technology that can address the risks confronting subsurface water resources while simultaneously promoting groundwater recharge [4]. MAR is also an efficient sustainable urban water management technique [5]. The technology involves the intentional injection of reclaimed water into an aquifer for subsequent retrieval or the achievement of environmental objectives [6]. Various built recharge structures are used to facilitate the accelerated replenishment of natural underground reservoirs [7]. An example of a widely embraced MAR technology in arid and semi-arid regions is rainwater harvesting (RWH). RWH involves the collection and storage of storm runoff to augment the water supply for domestic and agricultural purposes while also replenishing local aquifers [8]. Infiltration basins serve as another effective method of MAR, relying on the infiltration process whereby water is spread over wide areas using basins, dams, and channels to recharge an aquifer [9,10]. In cases where confined layers are present, deep injection wells are commonly used to convey water to the aquifer system [11]. On the other hand, vadose zone wells (i.e., drywells) offer advantages over other MAR techniques in recharging shallow urban aquifers due to their minimal land requirement and cost-effectiveness [12,13,14]. In general, Bouwer [15] asserted that the utilization of MAR systems offers numerous advantages compared to traditional water storage methods such as large dams. The artificial recharge of underground water via MAR technology not only minimizes water loss due to evaporation but also enables a smaller spatial footprint requirement and enhanced long-term storage capacity with lower costs [16].
In the context of MAR technologies, identifying suitable locations for recharge structures during the design stage of MAR projects is a challenging task, given the multiple factors involved. Consequently, the research community has developed several approaches to tackle this issue. For example, Anbazhagan and Ramasamy [17] utilized geophysical methods (particularly resistivity surveys) to identify appropriate sites for artificial recharge. Similarly, Christy and Lakshmanan [18] employed electrical resistivity and ground-penetrating radar methods to delineate feasible locations for MAR structures in Chennai, India, aiming to address saltwater intrusion. Brown et al. [19], on the other hand, proposed a statistical analysis-based methodology to determine suitable sites for Aquifer Storage and Recovery (ASR) projects. In a different vein, Zaidi et al. [20], Tiwari et al. [21], and Ahirwar et al. [22] utilized an integrated geographic information system (GIS) and remote sensing methods to pinpoint promising regions for MAR activities. Nevertheless, the application of coupled GIS and multi-criteria decision analysis (GIS-MCDA) techniques can be regarded as one of the most embraced approaches for delineating sites favorable for MAR applications [23,24,25,26,27,28]. Hence, the ability of a GIS to store and process spatial data [1] and the robustness of MCDA methods in resolving multi-tiered decision-making problems [29] make their integration a promising alternative for identifying regions promising for MAR implementation.
Sallwey et al. [30] noted that the most commonly adopted GIS-MCDA approaches in the relevant literature for determining the importance levels of decision parameters and combining thematic layers are the analytic hierarchy process (AHP) and the weighted linear combination (WLC), respectively. To exemplify this, Kazakis [25] applied AHP-WLC to identify potential MAR sites as a means to mitigate and prevent saltwater intrusion in a coastal aquifer. Similarly, Zhang et al. [31] employed the AHP to assess the significance levels of criteria and employed the WLC to combine thematic layers, leading to the production of a final MAR suitability map for the water-stressed West Coast region of South Africa. The authors claimed that the combination of the AHP with WLC offers an effective approach to address spatial decision-making problems. In line with this, Shadmehri Toosi et al. [32] applied the AHP to derive the weights of the utilized decision parameters (namely, slope, rainfall, soil type, soil depth, land use, and drainage density) to delineate suitable sites for RWH structures. However, it should be noted that some scholars have highlighted the limitations of the standard AHP, particularly in incorporating the inherent fuzziness of human judgment [33]. Consequently, hybridizing fuzzy set theory and the AHP technique has been widely acknowledged by various researchers in recent attempts [4,10,34]. In addition, various MCDA-based prioritization techniques, such as the evaluation based on distance from average solution (EDAS), simple additive weighting (SAW), and the technique for order of preference by similarity to ideal solution (TOPSIS), have been utilized in the pertinent literature [35,36,37,38].
From a spatial aspect, the assessment of MAR potential has been extensively conducted at the watershed or district scale. However, it is imperative to emphasize the significance of examining MAR potential at a countrywide level, particularly in regions where MAR is a relatively new concept. In this regard, Mahmoud et al. [39] developed an RWH site suitability map for Egypt at a national scale to promote the sustainability of subsurface resources and address water scarcity across the country. Likewise, Mahmoud and Tang [40] generated an RWH potential map for the entire United Kingdom considering divergent decision attributes, including rainfall, slope, curve number, land cover/use, and soil texture. Bonilla Valverde et al. [24] mapped the potential of MAR technology (notably surface-spreading techniques) in Costa Rica at the country scale based on four decision factors, namely, the hydrogeological aptitude, slope, soil texture, and drainage network density. Kadhem and Zubari [41] prioritized the Kingdom of Bahrein in terms of MAR suitability as a means to address the groundwater table decline and saltwater intrusion in the country. Notably, Mati et al. [42] conducted the largest-scale assessment of MAR potential, evaluating RWH technology for ten African countries, namely, Botswana, Ethiopia, Kenya, Malawi, Mozambique, Rwanda, Tanzania, Uganda, Zambia, and Zimbabwe.
This study entailed the assessment of MAR potential in a hot and arid region through the utilization of a novel hybrid fuzzy AHP and TOPSIS approach. The proposed framework was applied to the water-stressed country of Djibouti, where ensuring the long-term sustainability of groundwater resources is of paramount importance for the local population’s survival. The research objectives can be summarized as follows: (i) conduct a comprehensive assessment of MAR potential in Djibouti for the first time, (ii) identify a comprehensive set of surface, environmental, and subsurface decision factors based on an extensive review of the literature to determine regions suitable for MAR activities, (iii) initiate a pairwise comparison questionnaire involving experts from various backgrounds and ensure the consistency of the preferences of the participants, (iv) employ the fuzzy AHP algorithm to determine the importance levels of the decision criteria while utilizing the TOPSIS method to prioritize the study area based on different degrees of suitability, (v) perform a two-stage sensitivity analysis to assess the stability and reliability of the proposed model, which represents an initial attempt in the existing literature, and (vi) utilize the web-based INOWAS tool to select feasible MAR techniques for the study region and assess their suitability in terms of their implementation in Djibouti.

2. Materials and Methods

The present research aimed to establish a robust decision framework encompassing the adoption of different MCDA techniques. Hence, a thorough literature survey was conducted first, leading to the extraction of decision criteria that influence the determination of (non-) suitable MAR regions. To quantitatively evaluate the corresponding decision criteria, this study introduced the fuzzy AHP algorithm with the triangular fuzzy membership function. It is important to note that, along with dealing with interpersonal uncertainty by adopting fuzzy set theory, the reliability of the subjective judgments of the experts and the stability of the entire criteria weighting process were ensured by consistency ratio checks and a sensitivity analysis, respectively. Given the necessity to mimic real-world conditions in this analysis to extrapolate the results that can potentially be applied in regions of interest, this research further integrated the criteria weights with the current values of the decision criteria in the TOPSIS analysis in order to designate the MAR-suitable regions in Djibouti. The research steps considered in this study are depicted in Figure 1.

2.1. Study Area

Djibouti, situated in the Horn of Africa, is the smallest country in the region, with a land area of approximately 23,200 km2 (Figure 2). The country has a current estimated population of one million and shares terrestrial borders with Eritrea, Ethiopia, and Somalia, as well as maritime borders with Yemen. The climate of Djibouti is characterized by hot and arid conditions, with an annual mean rainfall of 150 mm and temperatures ranging from 20 °C to 30 °C in winter (October to April) and 30 °C to 45 °C in summer (May to September) [43]. In Djibouti, the principal groundwater reserves consist of fractured volcanic aquifers, primarily represented by three geological formations: the Dalha basalts, the Stratoid basalts, and the Mabla rhyolites [44,45]. These geological formations serve as the main sources of groundwater supply, characterized by a fractured nature and high permeability, allowing for substantial quantities of groundwater to be extracted. Furthermore, sedimentary aquifers with a transmissivity between 5.6 × 10−6 and 1.3 × 10−3 m2/s [44] are present along some watersheds and the eastern coastal part of the study region.
The ability to meet the water demand of the country heavily relies on groundwater due to the absence of permanent surface water resources. This overdependence has led to the widespread use of abstraction wells without proper planning or sustainable management practices [46]. Consequently, this excessive groundwater extraction, coupled with reduced recharge rates due to climate variability, poses substantial challenges to subsurface resources (i.e., depletion, contamination, and saltwater intrusion). In light of these challenging circumstances, investigating the feasibility of MAR in Djibouti assumes great importance, as it presents a viable solution to mitigate water shortage and achieve long-term sustainability.

2.2. Data Source and Layer Processing

After performing an extensive review of the existing literature and consulting with experts in the field, a comprehensive analysis resulted in the identification of nine pertinent parameters for delineating regions suitable for MAR activities. These selected criteria were subsequently categorized into three primary clusters, each encompassing three criteria, in order to establish an improved hierarchical framework (Table 1). The surface cluster consisted of the slope, soil texture, and curve number, while rainfall, the normalized difference vegetation index (NDVI), and drainage density were incorporated into the decision framework under the umbrella of the environmental cluster. Lastly, the subsurface cluster comprised depth to groundwater, geology, and groundwater quality (specifically the electrical conductivity in this investigation). Additionally, a concise description of the influence of each criterion in the MAR potential mapping decision framework can also be found in Table 1.
The generation of a spatial thematic layer for each factor involved in the decision framework represents an important step in MAR mapping. In line with this, a raster layer was generated for each criterion based on various data sources. For instance, slope and drainage density maps (depicted in Figure 3a and Figure 3f, respectively) were derived from the ALOS World 3D (AW3D) digital elevation model with a 5 m resolution [47]. The available slope analysis tools in QGIS were used to generate the slope layer of the study region in percent, while the drainage density raster layer was obtained by applying the line density interpolation algorithm. Soil texture was obtained from the Harmonized World Soil Database version 2.0 (HWSD v2.0) [48], and the resulting map is given in Figure 3b. The curve number map (Figure 3c) of the study region was clipped from the global curve number (GCN250) dataset provided by Jaafar and Ahmad [49]. The GCN250 dataset has a resolution of 250 m and can be downloaded in raster format. To establish the rainfall distribution throughout the country, the average annual rainfall at various stations was acquired from Dabar et al. [50], and the precipitation raster layer (Figure 3d) was generated using the inverse distance weighting (IDW) interpolation method. Regarding the calculation of the NDVI illustrated in Figure 3e, Sentinel 2 imagery was employed, and the following formula was applied:
NDVI = NIR RED NIR + RED
where NIR and RED represent the near-infrared and red bands, respectively. The acquisition of Sentinel imagery from Google Earth Engine and the necessary computations to obtain the NDVI were performed using the Geemap Python package [51]. The geological representation of the country was digitized from a technical report [46], while the layers depicting depth to groundwater and electrical conductivity were generated using IDW interpolation based on raw data obtained from the Ministry of Agriculture, Water, Fisheries and Livestock, in charge of Fishery Resources (MAEPE-RH) of Djibouti. Hence, Figure 3g–i illustrate the depth to groundwater, geology, and electrical conductivity maps, respectively. It is important to highlight that all geoprocessing tasks in this research were conducted using QGIS 3.28 [52].
During this stage, the need for standardizing the decision layers becomes imperative due to their diverse nature in this study. The existing literature encompasses various standardization techniques, including fuzzy membership functions, stepwise functions, AHP, and the utilization of rates/scores for feature classes [23,24,31,32,53,54,55]. Consequently, this study employed the rating technique to reduce the complexity of the proposed framework, whereby a score of 5 was assigned to cells exhibiting high suitability for MAR potential, while a rating of 1 was attributed to the least suitable cells (Table 1).
Table 1. List of selected decision criteria, accompanied by their corresponding descriptions, along with the ratings assigned to the criteria classes.
Table 1. List of selected decision criteria, accompanied by their corresponding descriptions, along with the ratings assigned to the criteria classes.
ClustersCriteriaIDClassesRatingSource Reference
SurfaceSlope a (%)S0–25AW3D [47]
2–54
5–103
10–252
>251
Soil texture bSTEutric Fluvisols (FLeu)5HWSD v2.0 [48]
Petric Gypsisols (GYp)4
Eutric Leptosols (LPeu)2
Lithic Leptosols (LPli)2
Haptic Solonchaks (SCh)1
Curve number cCN75–855Jaafar and Ahmad [49]
67–754
62–673
55–622
40–551
EnvironmentRainfall d (mm)Rl153–2205Dabar et al. [50]
134–1534
118–1343
95–1182
43–951
Normalized difference eNDVI0.10–0.245Sentinel 2 [56]
vegetation index 0.05–0.104
0.02–0.053
0–0.021
Drainage density f (km/km2)DD0–4.135AW3D [47]
4.13–7.104
7.10–10.173
10.17–14.122
14.12–24.471
SubsurfaceDepth to groundwater g (m)DG0.7–131MAEPE-RH
13–215
21–353
35–672
>671
Geology hGSedimentary formations5MAEPE-RH [46]
Recent basalt (1 Ma-actual)3
Stratoid basalt (3.4–1 Ma)2
Golf basalt (3.4–1 Ma)2
Somali basalt (9–3.4 Ma)2
Dalha basalt (9–3.4 Ma)2
Adolei basalt (25Ma)2
Stratoid rhyolite1
Mabla rhyolite1
Cretaceous–Jurassic base1
Groundwater quality i (µS/cm)EC200–10005MAEPE-RH
1000–21004
2100–26003
2600–37002
>37001
Notes: a Slope serves as a highly utilized decision criterion in the context of MAR mapping due to its significant influence on the convergence and divergence of runoff water, ultimately affecting the infiltration capacity [24]. b Rajasekhar et al. [57] noted that the soil type plays a critical role in regulating both infiltration rates and the potential generation of runoff. Specifically, soils with a high clay content tend to display diminished infiltration rates and increased runoff, whereas soils with a high sand content tend to exhibit enhanced infiltration rates and reduced runoff. c As a dimensionless index that characterizes the soil’s ability to absorb water, the curve number can be used to indirectly estimate the volume of runoff that can be harvested in a particular area for MAR usage [58]. d Rainfall plays a pivotal role in MAR mapping since rainwater represents the main source of water for MAR projects worldwide [59]. e According to Ansems et al. [60], high-NDVI regions could be indicative of the temporal availability of water and thereby have the potential to reclaim large volumes of water. f Drainage density is inversely proportional to permeability such that areas with a high drainage density indicate the presence of low-permeability rock, whereas a low drainage density suggests the presence of more permeable rock [61]. g The depth to groundwater influences the feasibility of the MAR project as well as the recharge rates [20]. h The rock type prevalent in a given region plays a vital role in governing the movement and distribution of groundwater [62]. i Injecting reclaimed water into a poor-groundwater-quality region could jeopardize the MAR benefit; therefore, it is important to include groundwater quality parameters in the decision framework [58].

2.3. MCDA Rationale

MCDA is a technique used in complex decision-making processes. MCDA methods are generally employed to evaluate alternative options by considering multiple criteria or objectives, and rather than providing an ideal solution, they enable the extraction of a compromise solution [63]. With the MCDA rationale, the first step involves identifying criteria to consider in the decision-making process, which represent the metrics used to evaluate the alternatives. Subsequently, weighting is applied to assign importance degrees to the identified criteria, reflecting their relative significance. Following that, each alternative is subjected to a performance evaluation based on the predefined criteria involving the objective evaluation of the alternatives’ performance. Finally, by utilizing the corresponding outcomes, an aggregation method is employed to determine the overall performance of the alternatives, aiding in ranking the alternatives and identifying the optimal solution(s). MCDA is applied in various domains, including business [64], public policy [65], environmental management [66], and healthcare [67] projects. Different techniques and tools are available within the MCDA framework, including the AHP, analytic network process (ANP), TOPSIS, and preference ranking organization method for enrichment evaluations (PROMETHEE). However, the choice of technique depends on the specific characteristics of the decision problem, the available data, and the preferences of the decision-makers. Overall, MCDA provides a structured and systematic approach to decision making, enabling stakeholders to consider multiple criteria and objectives in a comprehensive manner. It supports more informed and rational decision making, especially in complex situations where there are trade-offs and conflicting preferences among different criteria.

2.3.1. Description of the Fuzzy AHP Algorithm

The AHP, introduced by Saaty [68], is frequently employed to construct a hierarchical framework that captures intricate problem structures, aligning with the subjective assessments of experts. This methodological framework enables a detailed and granular analysis of complex problems, allowing for the incorporation of the subjective judgments of experts while maintaining scientific rigor. The AHP employs a systematic approach to hierarchically structure a problem, beginning with a primary goal and progressively delving into criteria, sub-criteria, and subsequent levels, ultimately encompassing a comprehensive set of alternatives [69]. This hierarchical representation offers experts a comprehensive understanding of the intricate relationships inherent in the context, enabling a holistic view of the problem. Furthermore, it aids in evaluating the comparability of elements at the same level, allowing for meaningful comparisons and assessments [33]. Subsequently, elements are subjected to pairwise comparisons on a nine-level scale to derive their relative weights, facilitating [70,71] the quantification of their respective importance in the decision-making process. Although the primary objective of the AHP is to incorporate expert knowledge, the classical AHP method falls short in addressing the non-numeric uncertainties inherent in human thinking processes. Consequently, the implementation of the fuzzy AHP to provide more robust solutions has gained much more interest in the recent decade. The fuzzy AHP method integrates the principles of fuzzy sets (initially introduced by Zadeh [72]) to address the intricate and ambiguous nature of the environment and experts’ preferences. The fuzzy AHP tackles the challenges posed by vagueness and imprecision in decision making by enabling fuzzy linguistic assessments and the representation of interpersonal uncertainties through fuzzy numbers. By capturing and quantifying the inherent ambiguity, the fuzzy AHP offers a more reliable and comprehensive approach to handling complex problems and supports more informed decision-making processes. When decision-makers evaluate criteria and alternatives, they can utilize natural language expressions alongside precise numerical values. Therefore, the fuzzy AHP method closely resembles human thinking and perceptions. As a result, it has been consistently employed by numerous researchers across various fields [70,71,73].
In the initial step of utilizing the fuzzy AHP, the decision hierarchy structure is established by incorporating all the criteria, including primary criteria, sub-criteria, goals, and alternatives. To determine the order, each pair of items is compared using a nine-point significance scale. The steps of the fuzzy AHP approach are commonly arranged in the following manner:
  • After the hierarchical structure is revealed, decision-makers construct binary comparison matrices in accordance with their perspectives. These matrices encompass the relative evaluations and favored choices among components at each level of the hierarchy. The reciprocals of linguistic variables regarding the importance degrees of the criteria are incorporated into the preferences of the experts who attended the surveys. Thus, l i j , m i j , and u i j , indicating the lower, mean, and upper widths of the pairwise judgments of the experts for criterion i compared to criterion j, respectively, are determined (Table 2).
  • In the fuzzy AHP approach, an additional step is implemented to verify the consistency of experts’ pairwise comparisons. This is achieved by calculating the consistency ratio (CR), where CR values exceeding 0.1 indicate inconsistent judgments made by respondents, while CR values below the threshold indicate a more consistent set of expert preferences. The following expression can be utilized to determine the CR values, thus assessing the level of consistency in the decision-making process.
    C R = λ m a x n n 1 R I
    where λ m a x and n denote the maximum eigenvalue of the matrix and the number of criteria in the matrix, respectively. In addition, R I is the random index, which was introduced by Saaty [74] and is based on the size of the matrix (Table 3).
  • In this step, the fuzzy equivalents of each linguistic variable are calculated. Equation (3) outlines the method for determining the lower ( l i j k ), mean ( m i j k ), and upper ( u i j k ) widths of the fuzzy equivalents using the triangular membership function.
    l i j = ( k = 1 K l i j k ) 1 / K ;   m i j = ( k = 1 K m i j k ) 1 / K ;   u i j = ( k = 1 K u i j k ) 1 / K
    in which K is the total number of respondents.
  • To address the inherent vagueness and uncertainty in experts’ judgments (Table 4), Chang’s [75] extent analysis was employed. In this approach, crisp mathematical notations were utilized to obtain fuzzy quantities. The object set, represented by X = x 1 ,   x 2 , x 3 , , x n , and the goal set, denoted by U = u 1 ,   u 2 , u 3 , , u n , were considered in the extent analysis. For each goal, denoted by u i , extent analysis values represented by m are obtained for each object.
    l i j = ( k = 1 K l i j k ) 1 / K ;   m i j = ( k = 1 K m i j k ) 1 / K ;   u i j = ( k = 1 K u i j k ) 1 / K
    To calculate M g i j , the fuzzy extent analysis M value addition operation is performed on the matrix. This operation involves adding each triangular fuzzy number (TFN) in each row of the matrix using the addition operation, as described in Equation (5).
    j = 1 m M g i j = j = 1 m l j , j = 1 m m j , j = 1 m u j
    with i = 1 ,   2 ,   , n . The score j = 1 n j = 1 m M g i j 1 is obtained by calculating the sum of the entire triangular fuzzy number set M g i j j = 1 ,   2 ,   , m .
    j = 1 n j = 1 m M g i j = j = 1 n j = 1 m l j , j = 1 n j = 1 m m j , j = 1 n j = 1 m u j
    The inverse of the initial equation can be computed using the formula presented in Equation (7).
    j = 1 n j = 1 m M g i j 1 = 1 i = 1 n u 1 , 1 i = 1 n m 1 , 1 i = 1 n l 1
    A comparative calculation is performed to assess the level of possibility between fuzzy numbers. This comparison is utilized to determine the weight value assigned to each criterion. When comparing two triangular fuzzy numbers M 1 = l 1 ,   m 1 ,   u 1 and M 2 = l 2 ,   m 2 ,   u 2 , where the probability level S2S1, a definition can be established.
    M 2 M 1 = 1 , i f   m 1 m 2 0 , i f   l 1 l 2 l 1 u 2 m 2 u 2 m 1 l 1 , f o r   o t h e r s
    To compare M 1 and M 2 , it is necessary to calculate the values V ( M 1 M 2 ) and V ( M 2 M 1 ). Once the fuzzy synthetic values have been compared, the minimum value is determined using Equation (9).
    D A i = m i n V S i S k
    For each k value ranging from 1 to n, where k i, the weight vector is calculated to facilitate the interpretation of the defined criteria.
    W = d A 1 , d A 2 , , d A n   T
    where A i i = 1 , 2 , n is n elements, and d’( A i ) is the score describing each decision attribute of the compared options.
Subsequently, the weights are normalized using Equation (11) to transform the values in the weight vector into analog weights. This normalization process ensures that the weights consist of non-fuzzy numbers.
d A i = d A i i = 1 n d A i                                 f o r   i = 1 , 2 , , n  
5.
The last step is considered crucial in determining the degrees of importance for the considered criteria [76]. Therefore, a sensitivity analysis was conducted to examine the variations in criteria importance based on different degrees of fuzziness. The initial degree of fuzziness in the adopted FAHP method was set to 1, determined by the distances between l, m, and u values (Table 2). Additionally, five additional fuzziness degrees (1.2, 1.4, 1.6, 1.8, and 2.0) were explored in the current study. Consequently, if the order of importance remains relatively unchanged, it can be concluded that the decision analysis framework yields reliable results and is not significantly influenced by changes in fuzziness degrees [77].

2.3.2. Description of the TOPSIS Algorithm

The TOPSIS method is a well-established decision-making technique that has been widely used in various fields. According to Hwang and Yoon [78], TOPSIS is based on the concept of ranking alternatives by their relative proximity to the ideal solution. On the one hand, the positive ideal solution represents the alternative that maximizes benefit criteria and minimizes cost criteria. On the other hand, the negative ideal solution minimizes benefit criteria and maximizes cost criteria. The selection of the most suitable alternative is performed based on its proximity to the positive ideal solution and distance from the negative ideal solution. The TOPSIS method provides a systematic approach to decision making, considering multiple criteria simultaneously. Overall, the TOPSIS method has proven to be effective in evaluating and ranking alternatives, and its application has been documented in numerous scientific studies [79,80,81].
The TOPSIS essentially contains six steps:
  • Defining the Decision Matrix: A decision matrix is formulated, encompassing all available alternatives along with their corresponding performance values on various criteria. The decision matrix is typically represented as an m × n matrix, where m is the number of alternatives and n is the number of criteria. The decision matrix according to the TOPSIS method is shown in Equation (12).
    A i j = a 11 a 12 a 1 n a 21 a 22 a 2 n . . . . . . . . a m 1 a m 2 a m n  
    In this study, the generation of the decision matrix involved several steps. Initially, a rectangular grid layer with a spatial resolution of 500 m × 500 m was created and subsequently clipped with the study area vector layer to confine the analysis within the defined study boundaries. Then, the zonal statistics tool was employed to calculate the mean value for each pixel, considering the nine decision layers as input raster layers. Consequently, the outcome of this process yielded a decision matrix consisting of 90,177 rows and 9 columns.
  • Normalizing the Decision Matrix: The decision matrix is normalized to eliminate any scale differences among the criteria. This step ensures that all criteria are given equal weightage in the decision-making process. Various normalization methods can be used, such as min–max normalization or vector normalization. The normalization of the decision matrix can be calculated through the formula depicted in Equation (13).
    r i j = a i j k = 1 m a k j 2                       i = 1 ,   2 ,   , m ; j = 1 ,   2 ,   ,   n  
  • Assigning Weights to the Criteria: The relative importance or weights of the criteria are determined. The weights reflect the significance of each criterion in the decision-making process. The determination of weights can be subjective, based on expert judgment, or derived using mathematical techniques, such as the analytic hierarchy process (AHP) or Entropy Weight Method. At this stage, the weighted decision matrix is obtained by multiplying the data in the normalized decision matrix, obtained in the second step of the TOPSIS method, by the weight values determined through the previously conducted weighting method. The sum of the data obtained from the weighting method must be equal to 1.
  • Determining the Positive and Negative Ideal Solutions: The positive ideal solution and the negative ideal solution are identified based on the maximum and minimum values, respectively, for each criterion. The positive ideal solution represents the alternative that achieves the maximum benefit and the minimum cost, while the negative ideal solution represents the alternative that minimizes the benefit and maximizes the cost. The TOPSIS method assumes that each criterion exhibits a monotonically increasing or decreasing trend. To determine the ideal solution set, the maximum value of the column in the weighted decision matrix is selected. If the criterion is cost-oriented or in a minimization direction, the smallest criterion is chosen. The relevant formula for the ideal solution set is shown in Equation (14).
    A + = m a x v i j j   J , m i n v i j j   J
    In the negative ideal solution set, the smallest values of the data in the columns containing criterion values in the weighted decision matrix are examined. The formula for the negative ideal solution set is shown in Equation (15).
    A = m i n v i j j   J , m a x v i j j   J
  • Calculating Euclidean Distances: The Euclidean distance between each alternative and the positive and negative ideal solutions is calculated. The Euclidean distance represents the overall proximity or distance of each alternative to the ideal solutions in the multi-dimensional criteria space. As a result of this process, the deviation values for the alternatives are defined as the ideal separation ( S 1 + ) and negative ideal separation ( S 1 ) measures. The formulas for ideal separation and negative ideal separation are shown in Equations (16) and (17).
    S 1 + = j = 1 n v i j v j + 2  
    S 1 = j = 1 n v i j v j 2  
  • Calculating the Proximity to Ideal Solutions: The relative proximity of each alternative to the positive and negative ideal solutions is determined. This can be achieved by calculating the relative closeness coefficient, which is the ratio of the distance from the negative ideal solution to the sum of the distances from the positive and negative ideal solutions. The calculation of the relative closeness to the ideal solution is shown in Equation (18).
    C C i = S 1 S 1 + S 1 +  
    The obtained C C i value takes a value in the range of 0 ≤ C C i ≤ 1. C C i = 1 indicates the absolute proximity of the alternative to the ideal solution, while C C i = 0 indicates the absolute proximity of the alternative to the negative ideal solution. The ranking of alternatives is determined by sorting the obtained C C i values in descending order, indicating their level of importance.

3. Results

3.1. Criteria Weighting

This research has explored the role of different criteria in assessing the MAR mapping based on three dimensions, i.e., surface, environment, and subsurface. The fuzzy AHP technique was employed not only to reveal the contribution of each dimension but also to demonstrate the impact of nine different criteria on the MAR suitability assessments. Hence, the present study utilized the outcomes of the judgments of the expert group (Table 4) to compute the weights of both the main clusters and the corresponding criteria. To ensure the reliability of the obtained results, consistency ratios were initially investigated. It is especially worth mentioning that obtaining CR values smaller than 10% indicates that the decision framework is consistent, and further computations, including criteria weighting, can be pursued. Figure 4 reveals that all the experts were quite consistent in pairwise comparisons. In qualitative assessments, the experts may be confused about their results due to the complex nature of the pairwise comparison matrix, leading to inconsistent evaluations. However, this research refined the list of criteria (3 × 3) and the main clusters (1 × 3), and the experts yielded consistent results in the first round of the surveys.
Table 5 highlights that, among the three main clusters, the environment dimension (with a weight of 48.85%) was determined as the most significant, followed by the surface (30.36%) and subsurface (20.79%) dimensions, respectively. Table 5 also includes the local and global weights of the criteria considered. Accordingly, in the surface cluster, the slope was the most determinant criterion, as it gained the highest weight with 43.92%, while the curve number and soil texture followed the slope with 31.27% and 24.81%. In addition, the rainfall criterion (62.80%) outperformed the NDVI (24.33%) and drainage density (12.86%) with respect to the criteria weights calculated through the fuzzy AHP. The local criteria weights regarding the subsurface dimension depict a more homogeneous distribution among the three considered criteria, i.e., groundwater quality (37.27%), geology (33.41%), and depth to groundwater (29.32%), indicating that each criterion has no distinctive impact on the determination of suitable MAR sites.
Furthermore, global weights, which directly indicate the contribution of each criterion to the final suitability assessment of the study region, are also presented in Table 5. According to the results obtained, rainfall, having a global weight of 30.68%, outperformed its counterparts. In addition, rainfall was followed by the slope (with a global weight of 13.33%) and NDVI (with a global weight of 11.89%). On the other hand, among the nine decision criteria, the fuzzy AHP assessments underestimated the role of depth to groundwater, drainage density, and geology, as the lowest weights were assigned to these factors.

3.2. Sensitivity Analysis

In this study, a sensitivity analysis with respect to different fuzziness degrees was implemented to ensure the stability of the obtained results as well as underscore the transparency of the criteria weighting attempts. Hence, in addition to dealing with interpersonal uncertainty by utilizing the fuzzy AHP rationale [82], the robustness of the criteria assessments is ensured via the sensitivity analysis conducted through different fuzziness degrees (i.e., 1.2, 1.4, 1.6, 1.8, and 2.0). This method was embraced due to its previous implementations [83]. Fuzziness degrees capture the level of ambiguity or uncertainty associated with linguistic judgments or criteria weights in the AHP. A sensitivity analysis helps understand the impact of variations in fuzziness degrees on the final results, enabling better management and the interpretation of uncertainty (Figure 5). The most significant aspect of assessing the sensitivity plot is focusing on variations in the criteria importance order. The figure reveals that, with regard to the alterations in fuzziness degrees, the fuzzy AHP results yielded consistent outcomes, as there are no considerable variations in criteria rankings. Therefore, one can conclude that the factors having a significant impact on model outcomes account for the similar importance considering various fuzziness degrees. The level of confidence/agreement in the decision process is ensured with the sensitivity analysis, and it is highlighted that the quality of the experts’ preferences is aligned with the consistent decision context (Figure 5).
In this study, an additional sensitivity analysis approach was conducted. To accomplish this aim, the weight of each criterion was shifted to another criterion, and the obtained results were compared with the baseline scenario (i.e., originally obtained criteria weights) with respect to the root mean square error (RMSE) computations (Table 6). The results showed that the vast majority of the scenarios yielded considerably lower RMSE values, except when Rl (i.e., rainfall) was exchanged with its counterparts. The corresponding results can be explained by the fact that rainfall was originally overestimated by the experts, illustrating its high influence on the final MAR suitability decision.

3.3. MAR Potential Mapping

The current research analyzed the suitability of MAR potential in Djibouti using the TOPSIS approach. To accomplish this aim, the criteria weights obtained through the fuzzy AHP analysis were integrated with the actual values of the corresponding criteria in the field and then subjected to the TOPSIS algorithm, aiding in exploring the prioritization of suitable regions. Figure 6 encompasses two different sub-figures, in which (a) the first represents the map generated based on the closeness coefficients, whereas (b) the latter provides the reclassification of the suitability assessments based on the closeness coefficients using the Jenks natural break method. In addition, the higher the closeness coefficient, the higher the suitability of MAR implementation. Hence, one can conclude that the northern and western parts of the country (with orange and red colors) are less suitable compared to other regions (Figure 6a). Conversely, the middle- and south-eastern parts of the country, mainly covering the vicinity of the Gulf of Tadjourah, are more suitable for MAR applications. In addition to the western part of the country, south Djibouti demonstrated very high suitability for MAR implementation (Figure 6b). The corresponding results can further be explained by the decision criteria; for example, highly suitable parts are significantly correlated to high rainfall values (Figure 3d), which was found to be the most determinant decision criterion.
Furthermore, Figure 7 illustrates the distribution of the different suitability clusters based on the corresponding regions’ surface areas. According to the figure, one can conclude that nearly one-third of the country is moderately suitable for MAR implementation. Likewise, one-third of the country is suitable (high and very high suitability) for the corresponding sustainable solutions, whereas the remaining parts are not (low and very low suitability). The obtained results, on the one hand, showed that approximately 16.38% of the country demonstrates very high suitability for MAR implementation, and these regions are mainly located in the middle and southern parts. On the other hand, nearly 15.76% of the entire country is mainly located in the north-eastern part, which implies very low suitability for MAR solutions.

4. Discussion

4.1. Assessment of the Decision Criteria

The current research found that rainfall was the most influential decision criterion in assessing the MAR suitability of a region. Here, the corresponding result is consistent with nature, as it is well established that not only the amount of rainfall but also the frequency of rainfall events is significantly influential in terms of the recharge potential of the area of interest [84,85]. Accordingly, regions characterized by higher rainfall rates can be regarded as more favorable for the implementation of MAR techniques, in contrast to regions experiencing limited amounts of rainfall. Additionally, alongside the rainfall amount, the intensity of rainfall events holds paramount importance, as it directly affects the efficiency of recharge. For example, intense rainfall events can result in runoff, reducing the effectiveness of recharge efforts, whereas gentle and steady rainfall facilitates better infiltration and a higher recharge efficiency [86]. Given all these facts, taking the amount and characteristics of rainfall into account is essential when determining suitable regions, allowing for the better planning, implementation, optimization, and long-term sustainability of MAR initiatives. Table 5 also illustrates that the slope is the second most determinant factor in assessing the MAR suitability potential, with a global weight of 13.33%. Here, the existing literature further recognizes the critical role played by the slope factor such that not only hydro-morphological variabilities (such as infiltration rates, surface water flow, sedimentation, etc.) [87,88] but also the feasibility of the MAR techniques (engineering considerations, compatibility with land use, etc.) [89] are ascertained based on the topographical conditions. For instance, Fathi et al. [90] highlighted that steep slopes generally result in faster runoff, which limits the amount of water available for recharge, and gentle slopes, on the other hand, promote the penetration of water, allowing for higher recharge rates and the improved effectiveness of MAR techniques. Compatibility with real-world conditions is also critical, particularly considering that steep slopes can require additional efforts (stabilization, erosion control structures, etc.) [91] to ensure the stability and longevity of recharge facilities. Likewise, the slope impact may differentiate the MAR type that is planned to be installed. For instance, steeper slopes may be less suitable for certain recharge methods, such as infiltration basins [87,92], due to limited space and the potential for erosion. Gentle slopes provide more opportunities for siting and integrating recharge facilities within the existing landscape, considering factors such as land availability, land use restrictions, and environmental compatibility.
According to the criteria evaluations, the NDVI was found to be the third most influential factor, with a global weight of 11.89%. As it reflects the vegetation condition of the study region, the NDVI provides valuable insights into the suitability of MAR selection mechanisms. Since the NDVI measures the amount of live green vegetation in an area [93], it plays a vital role in enhancing infiltration and reducing runoff, making areas with higher NDVI values more conducive to recharge efforts. In addition, monitoring NDVI variabilities helps estimate water loss due to evapotranspiration, aiding in quantifying the water demand in regions. Consequently, one can assess the potential availability of water for recharge purposes. Focusing on NDVI values further aids in taking the infiltration capacity of soils [94], ecological considerations [95], and land use patterns [96] into account and thereby gaining a better overview regarding the suitability of MAR attempts. The fuzzy AHP analysis conducted based on pairwise comparisons underestimated the role of geology (global weight of 6.95%), drainage density (global weight of 6.28%), and depth to groundwater (global weight of 6.10%), as they ranked at seventh, eighth, and ninth, respectively. It is especially worth noting that these factors influence the selection of suitable sites for potential MAR implementations; however, their impact was found to be limited compared to other criteria incorporated into the decision framework.

4.2. Assessment of the Adopted Decision Framework and Its Limitations

This research introduces a holistic hybrid MCDA framework aimed at identifying favorable regions for MAR activities in the arid country of Djibouti as a response to the groundwater challenges faced by the nation. To determine the relative importance of individual decision parameters, the fuzzy AHP method was employed. Additionally, the TOPSIS technique was utilized to prioritize the study region and generate a countrywide MAR potential map. Although previous studies have separately applied the fuzzy AHP and TOPSIS methods [10,38], the integration of these approaches has been largely overlooked in the existing literature, with Mouhoumed et al. [4] being the only previous authors to attempt to adopt a similar combination. However, their framework was limited to a city-scale assessment covering an area of 218.53 km2 and, more importantly, was focused on a predefined MAR technology, namely, drywells. In contrast, this research focused on an assessment at the country scale, and rather than concentrating on a specific approach, all MAR techniques are targeted.
Refining the list of decision parameters poses a challenge in multi-tiered decision problems, and this challenge is also evident in the context of MAR potential mapping. As highlighted by Sallwey et al. [30], the choice of criteria for site delineation is influenced by the problem definition. The average number of selected decision factors in the relevant literature tends to be around eight [30]. However, the corresponding number of decision parameters tends to be lower in country-level cases, likely due to data availability and accessibility limitations in such large areas. Mahmoud et al. [39], Mahmoud and Tang [40], and Mati et al. [42], for instance, used a set of five criteria to map MAR at a country scale, while Bonilla Valverde et al. [24] selected four decision variables. However, it is worth noting that in a recent study conducted by Kadhem and Zubari [41] to identify optimal locations for MAR in Bahrain, eight decision criteria were utilized. The present research incorporated the largest criteria set (nine factors) into a country-level MAR potential assessment.
The validation of the proposed models further represents an additional burden in MAR site suitability studies [53]. Apart from the meticulous selection of the number and nature of parameters involved in the decision framework, assigning appropriate weights to individual parameters can also introduce ambiguity since they highly rely on the subjective judgments of the pairwise survey participants. The scarcity of reported successful MAR projects in many countries, particularly in regions where MAR is a relatively novel topic, renders commonly used validation techniques less applicable. Nevertheless, some scholars argue that a combined consistency check and sensitivity analysis are enough to validate proposed models [25,97]. In line with this, several sensitivity analysis approaches have been proposed so far, ranging from the adoption of different weighting schemes [53], the addition or removal of relationships in the case of the multi-influence factor (MIF) method [24], or changes in the linguistic quantifier when applying the ordered weighted averaging (OWA) technique [97,98]. While the fuzzy AHP algorithm incorporates the inherent ambiguity of human decision making, we also tested the robustness of our model by adjusting the degree of fuzziness within the fuzzy AHP framework. This application of sensitivity analysis by modifying the degree of fuzziness represents a rare endeavor in the literature, further enhancing the novelty of the proposed framework.
This research further compared the current attempt with previously applied studies, as depicted in Table 7. The table illustrates that most of the efforts regarding MAR implementation were devoted to a watershed-scale analysis, while problem statements, i.e., major objectives, were diverse among the research community. In addition, a diligent review of the pertinent literature also demonstrated that the vast majority of studies have utilized traditional MCDA techniques for both criteria weighting and alternative prioritization, such as the AHP and WLC, respectively. As mentioned before, the MCDA rationale is implemented based on the preferences of experts having divergent backgrounds regarding the topic of interest. Despite this fact, only two of the past studies, i.e., Sandoval and Tiburan [27] and Shadmehri Toosi et al. [32], provided details of the experts who attended pairwise comparison surveys. From a different aspect, although Kharazi et al. [36] and Itani et al. [99] used the preferences of local experts for their criteria weight assessment, the researchers presented limited information regarding their fields of expertise. Given the importance of providing transparency and clarity to MCDA attempts, this study further provided the job description, role, and experience in the relevant field of the experts who were incorporated into the decision-making framework. Still, such a clear demonstration of experts’ backgrounds has not been acknowledged in the pertinent literature. Along with the experts’ profiles, the number of experts contributing to the analysis is of critical importance in managing decision processes. Including a limited number of experts may restrict the generation of creative ideas, whereas a high number of experts may endanger the diversification of the criteria weights, especially in fuzzy AHP analysis. Hence, the present research benefited from the preferences of a total of eight experts, who all provided reliable outcomes based on consistency checks. Finally, one can also consider performing a sensitivity analysis in order to ensure the stability of the criteria weighting attempts. In this regard, only three of the past studies, i.e., Kazakis [25], Fuentes and Vervoort [53], and Itani et al. [99], conducted sensitivity analyses to ensure the stability of MCDA outcomes. At this point, it is especially worth noting that this study is also differentiated from its counterparts as it covers two different types of sensitivity analyses. First, this research sought to explore the impact of fuzziness degrees on the criteria weights, and secondly, it applied 36 additional MCDA analyses (i.e., 9 × 8/2) based on shifting the weight of one criterion to another. Such a comprehensive evaluation of reliability control regarding criteria weighting and alternative prioritization can be considered a first attempt in the pertinent literature. It is also significant to mention that changing the fuzziness degrees has little impact on criteria weighting, and a similar conclusion can be drawn for the subsequent sensitivity checks, as the computed RMSE values are considerably low.
Despite the valuable insights offered by this research in assessing MAR potential at a countrywide level, as well as the type of MAR technologies feasible in Djibouti, there are still some limitations that need to be covered in follow-up attempts. The model could be challenged at a watershed scale, resulting in much lower resolution MAR site suitability maps. The number of decision factors could also be increased in future research by incorporating different dimensions that encompass divergent criteria (such as socio-economic parameters) in MAR potential mapping. Other MCDA methods could be introduced for similar objectives regarding both the criteria weight assessment process (e.g., ANP) and the alternative prioritization phase (e.g., VIKOR and PROMETHEE, as well as their fuzzy variants). In addition, despite the TOPSIS method possessing a more advanced mathematical foundation for prioritizing alternatives, in contrast to the commonly employed GIS overlay techniques found in the relevant literature, it is imperative to substantiate its superior accuracy compared to the latter techniques in effectively delineating locations suitable for MAR implementation. This can be achieved by comparing TOPSIS with WLC, WOA, and other MCDA-based prioritization techniques that address some of the limitations of TOPSIS, serving as a guideline for future research. Furthermore, along with comparing the performance of the proposed decision-making process with the traditional approaches (such as BWM or AHP), evaluating it based on more recent techniques, e.g., RANCOM [101], would not only strengthen the provided scheme but also ensure valuable insights in follow-up studies. In addition, future research can also incorporate probabilistic approaches (e.g., Monte Carlo simulations) into both criteria weighting and alternative prioritization in order to deal with intrapersonal uncertainties that exist in such decision frameworks. Finally, the current study performed two types of sensitivity analyses and employed similarity checks based on RMSE computations. However, different metrics (such as the WS rank similarity coefficient and weighted Spearman correlation coefficient) [101] can also be adopted for the extraction of ranking similarities and the performance evaluation of the sensitivity checks.

4.3. Feasible MAR Technologies in Djibouti and Practical Utilization of the Proposed Framework

MAR potential maps serve as effective tools for evaluating the artificial recharge feasibility of large areas. However, one limitation of such broad-scale MAR assessments is their inability to provide insights into the feasibility of specific MAR technologies within a focused region. For instance, in the case of Djibouti, where permanent rivers are absent, bank filtration technology is deemed unfeasible. Fortunately, several tools have been developed to enhance the practicality of MAR and address the complexities associated with its implementation. Among these tools, the web-based INOWAS platform stands out, providing a diverse array of functionalities aimed at facilitating MAR initiatives (such as groundwater flow/transport numerical and analytical simulation models, MAR databases, data-driven tools, etc.). Specifically, tool no. 6 within the platform enables the selection of appropriate MAR techniques based on specific characteristics and conditions of the area of interest, including factors such as the availability of water sources and other pertinent considerations. Notably, the development of this tool by the INOWAS research group is based on extensive investigations encompassing several hundred MAR studies conducted worldwide [102]. Therefore, the data-driven tool was used to assess the type of MAR techniques feasible in Djibouti, and the result is depicted in Figure 8. The findings revealed that Aquifer Storage and Recovery (ASR) and Aquifer Storage Transfer and Recovery (ASTR) technologies could be adopted to store reclaimed water in the aquifer for future recovery or to address environmental concerns such as saltwater intrusion. These techniques could potentially utilize the treated water from the wastewater treatment plant in Douda, having a daily capacity to treat 2000 m3 of water [103], and from the Ambouli dam, storing a considerable amount of water that can be diverted for MAR activities alongside its flood mitigation purpose. Additionally, surface-spreading techniques (such as infiltration ponds, ditches, and furrows, as well as barriers, bunds, etc.) requiring rainwater as target sources were found to be feasible for recharging local aquifers. Furthermore, drywells and flooding MAR techniques show promise in recharging shallow aquifers in the country, aiming to achieve groundwater sustainability and address recurring droughts.
As mentioned earlier, this research contains a comprehensive MCDA analysis containing the implementation of fuzzy AHP and TOPSIS algorithms to identify the suitability of MAR techniques in Djibouti. To provide valuable insights to both policymakers and decision-makers, the practical utilization of such a framework is needed. Figure 9 provides an example of how the corresponding decision-making framework is used in conjunction with real-world implementation scenarios. As depicted in the figure, these analyses start by gathering relevant data and are followed by pairwise comparison surveys to explore the importance of differences in the criteria considering potential application regions. It is therefore important to identify key stakeholders carefully, since the contribution of the criteria taken into account may vary based on the focus region. Once the judgments of the experts are collected, criteria weights are computed. Subsequently, the actual data, which can be obtained from either model outcomes, in situ measurements, or remote sensing products, were incorporated into criteria weights to determine the site suitability map for MAR implementation. The resultant maps illustrate the regions suitable for the application of MAR technologies (from very high to very low suitability). The final step of the analysis covers the determination of the type of MAR technique, which is beyond the scope of the current research. Once the decision is made and the selected MAR technique is implemented, the outcomes based on various aspects are evaluated to review the impact of the corresponding solutions. Finally, based on the performance evaluation of the employed MAR technique(s), the entire framework is repeated to ensure the continuous improvement of the decision-making scheme.

5. Conclusions

The existing literature recognizes the promising role played by the integration of GIS and MCDA techniques for conducting MAR potential investigations. Thus, establishing a robust GIS-MCDA scheme enables the assessment of large areas while mitigating the computational burden and financial constraints. Consequently, in this study, a coupled model combining the fuzzy AHP and TOPSIS was proposed to evaluate the MAR potential at a country level, considering Djibouti as the specific case study. A comprehensive set of criteria pertaining to the surface, environmental, and subsurface dimensions were identified, and their relative importance was determined using the fuzzy AHP algorithm. Subsequently, the TOPSIS technique was employed to prioritize the study area based on the weighted decision layers. A sensitivity analysis was also employed to ensure the robustness and stability of the decision framework.
The results of the analysis highlighted the significance of rainfall, the slope, and the NDVI as the most influential decision criteria in identifying regions suitable for MAR implementation. Conversely, the fuzzy AHP determined that depth to groundwater and the drainage density were comparatively less influential in the decision-making process. In terms of the prioritization of the study area, it was observed that approximately 10.63%, 23.20%, and 31.06% of the country, corresponding to approximately 2466 km2, 5382 km2, and 7206 km2, respectively, exhibited very high, high, and moderate suitability for hosting MAR activities. Furthermore, the sensitivity analyses conducted to evaluate the stability of the framework indicated its robustness, as there were no significant changes in the ranks of the criteria with respect to the various degree-of-fuzziness values, and considerably less variation was observed in RMSE values computed based on closeness coefficients.
Overall, the comprehensive approach proposed in this study empowers policymakers and stakeholders to identify and prioritize areas with promising MAR potential, facilitating informed decision making and the efficient allocation of resources for sustainable water management practices.

Author Contributions

R.M.M.: Conceptualization, Methodology, Software, Formal Analysis, Writing—Original Draft, and Visualization. Ö.E.: Conceptualization, Methodology, Software, Formal analysis, Writing—Original Draft, and Visualization. E.E.B.: Conceptualization, Methodology, Software, Formal Analysis, Writing—Original Draft, and Visualization. M.Ö.: Resources, Writing—Review and Editing, and Supervision. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by the African Center of Excellence in Logistics and Transport (CEALT).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Singh, L.K.; Jha, M.K.; Chowdary, V.M. Multi-criteria analysis and GIS modeling for identifying prospective water harvesting and artificial recharge sites for sustainable water supply. J. Clean. Prod. 2017, 142, 1436–1456. [Google Scholar] [CrossRef]
  2. World Bank Group. Climate Risk Country Profile: Djibouti; The World Bank Group: Washington, WA, USA, 2021; p. 20433. [Google Scholar]
  3. Broek, E.; Hodder, C.M. Towards an Integrated Approach to Climate Security and Peacebuilding in Somalia; Stockholm International Peace Research Institute: Solna, Sweden, 2022. [Google Scholar]
  4. Mouhoumed, R.M.; Ekmekcioglu, O.; Ozger, M. A hybrid MCDA approach for delineating sites suitable for artificial groundwater recharge using drywells. J. Hydrol. 2023, 620, 129387. [Google Scholar] [CrossRef]
  5. Page, D.; Bekele, E.; Vanderzalm, J.; Sidhu, J. Managed aquifer recharge (MAR) in sustainable urban water management. Water 2018, 10, 239. [Google Scholar] [CrossRef] [Green Version]
  6. IAH. IAH Commission on Managed Aquifer Recharge. International Association of Hydrogeologists. 2022. Available online: https://recharge.iah.org/ (accessed on 4 June 2023).
  7. World Bank Group. What the Future Has in Store: A New Paradigm for Water Storage; World Bank Group: Washington, WA, USA, 2023. [Google Scholar]
  8. Ammar, A.; Riksen, M.; Ouessar, M.; Ritsema, C. Identification of suitable sites for rainwater harvesting structures in arid and semi-arid regions: A review. Int. Soil Water Conserv. Res. 2016, 4, 108–120. [Google Scholar] [CrossRef] [Green Version]
  9. Alam, S.; Borthakur, A.; Ravi, S.; Gebremichael, M.; Mohanty, S.K. Managed aquifer recharge implementation criteria to achieve water sustainability. Sci. Total. Environ. 2021, 768, 144992. [Google Scholar] [CrossRef]
  10. Papadopoulos, C.; Spiliotis, M.; Pliakas, F.; Gkiougkis, I.; Kazakis, N.; Papadopoulos, B. Hybrid Fuzzy Multi-Criteria Analysis for Selecting Discrete Preferable Groundwater Recharge Sites. Water 2022, 14, 107. [Google Scholar] [CrossRef]
  11. Brown, C.J.; Weiss, R.; Verrastro, R.; Schubert, S. Development of an Aquifer, Storage and Recovery (ASR) Site Selection Suitability Index in Support of the Comprehensive Everglades Restoration Project. J. Environ. Hydrol. 2005, 13, 1–13. [Google Scholar]
  12. Edwards, E.C.; Harter, T.; Fogg, G.E.; Washburn, B.; Hamad, H. Assessing the effectiveness of drywells as tools for stormwater management and aquifer recharge and their groundwater contamination potential. J. Hydrol. 2016, 539, 539–553. [Google Scholar] [CrossRef] [Green Version]
  13. Händel, F.; Liu, G.; Fank, J.; Friedl, F.; Liedl, R.; Dietrich, P. Assessment of small-diameter shallow wells for managed aquifer recharge at a site in southern Styria, Austria. Hydrogeol. J. 2016, 24, 2079–2091. [Google Scholar] [CrossRef]
  14. Sasidharan, S.; Bradford, S.A.; Šimůnek, J.; Kraemer, S.R. Comparison of recharge from drywells and infiltration basins: A modeling study. J. Hydrol. 2020, 594, 125720. [Google Scholar] [CrossRef]
  15. Bouwer, H. Artificial recharge of groundwater: Hydrogeology and engineering. Hydrogeol. J. 2002, 10, 121–142. [Google Scholar] [CrossRef] [Green Version]
  16. Dillon, P.; Stuyfzand, P.; Grischek, T.; Lluria, M.; Pyne, R.D.G.; Jain, R.C.; Bear, J.; Schwarz, J.; Wang, W.; Fernandez, E.; et al. Sixty years of global progress in managed aquifer recharge. Hydrogeol. J. 2018, 27, 1–30. [Google Scholar] [CrossRef] [Green Version]
  17. Anbazhagan, S.; Ramasamy, S.M. Geophysical resistivity survey and potential site selection for artificial recharge in central Tamil Nadu, India. In Engineering Geology and the Environment, 2nd ed.; Taylor & Francis: Milton Park, UK, 1997; pp. 1169–1173. [Google Scholar]
  18. Christy, R.M.; Lakshmanan, E. Percolation pond as a method of managed aquifer recharge in a coastal saline aquifer: A case study on the criteria for site selection and its impacts. J. Earth Syst. Sci. 2017, 126, 66. [Google Scholar] [CrossRef] [Green Version]
  19. Brown, C.J.; Ward, J.; Mirecki, J. A Revised Brackish Water Aquifer Storage and Recovery (ASR) Site Selection Index for Water Resources Management. Water Resour. Manag. 2016, 30, 2465–2481. [Google Scholar] [CrossRef]
  20. Zaidi, F.K.; Nazzal, Y.; Ahmed, I.; Naeem, M.; Jafri, M.K. Identification of potential artificial groundwater recharge zones in Northwestern Saudi Arabia using GIS and Boolean logic. J. Afr. Earth Sci. 2015, 111, 156–169. [Google Scholar] [CrossRef]
  21. Tiwari, A.K.; Lavy, M.; De Maio, M.; Singh, P.K.; Mahato, M.K.; Amanzio, G. Identification of artificial groundwater recharging zone using a GIS-based fuzzy logic approach: A case study in a coal mine area of the Damodar Valley, India. Appl. Water Sci. 2017, 7, 4513–4524. [Google Scholar] [CrossRef]
  22. Ahirwar, S.; Malik, M.S.; Ahirwar, R.; Shukla, J. Identification of suitable sites and structures for artificial groundwater recharge for sustainable groundwater resource development and management. Groundw. Sustain. Dev. 2020, 11, 100388. [Google Scholar] [CrossRef]
  23. Aju, C.; Achu, A.; Raicy, M.; Reghunath, R. Identification of suitable sites and structures for artificial groundwater recharge for sustainable water resources management in Vamanapuram River Basin, South India. Hydroresearch 2021, 4, 24–37. [Google Scholar] [CrossRef]
  24. José, P.; Bonilla, V.; Clemens, B.; Mario, R.; Lisa, S.; Catalin, S. Application of a GIS Multi-Criteria Decision Analysis for the Identification of Intrinsic Suitable Sites in Costa Rica for the Application of Managed Aquifer Recharge (MAR) through Spreading Methods. Water 2016, 8, 391. [Google Scholar] [CrossRef] [Green Version]
  25. Kazakis, N. Delineation of Suitable Zones for the Application of Managed Aquifer Recharge (MAR) in Coastal Aquifers Using Quantitative Parameters and the Analytical Hierarchy Process. Water 2018, 10, 804. [Google Scholar] [CrossRef] [Green Version]
  26. Rajasekhar, M.; Ajaykumar, K.; Bhagat, V. Identification of artificial groundwater recharge zones in semi-arid region of southern India using geospatial and integrated decision-making approaches. Environ. Chall. 2021, 5, 100278. [Google Scholar] [CrossRef]
  27. Sandoval, J.A.; Tiburan, C.L. Identification of potential artificial groundwater recharge sites in Mount Makiling Forest Reserve, Philippines using GIS and Analytical Hierarchy Process. Appl. Geogr. 2019, 105, 73–85. [Google Scholar] [CrossRef]
  28. Tsangaratos, P.; Kallioras, A.; Pizpikis, T.; Vasileiou, E.; Ilia, I.; Pliakas, F. Multi-criteria Decision Support System (DSS) for optimal locations of Soil Aquifer Treatment (SAT) facilities. Sci. Total. Environ. 2017, 603–604, 472–486. [Google Scholar] [CrossRef] [PubMed]
  29. Malczewski, J.; Rinner, C. Multicriteria Decision Analysis in Geographic Information Science; Springer: New York, NY, USA, 2015; Volume 1. [Google Scholar]
  30. Sallwey, J.; Valverde, J.P.B.; López, F.V.; Junghanns, R.; Stefan, C. Suitability maps for managed aquifer recharge: A review of multi-criteria decision analysis studies. Environ. Rev. 2019, 27, 138–150. [Google Scholar] [CrossRef]
  31. Zhang, H.; Xu, Y.; Kanyerere, T. Site Assessment for MAR through GIS and Modeling in West Coast, South Africa. Water 2019, 11, 1646. [Google Scholar] [CrossRef] [Green Version]
  32. Toosi, A.S.; Tousi, E.G.; Ghassemi, S.A.; Cheshomi, A.; Alaghmand, S. A multi-criteria decision analysis approach towards efficient rainwater harvesting. J. Hydrol. 2019, 582, 124501. [Google Scholar] [CrossRef]
  33. Liu, Y.; Eckert, C.M.; Earl, C. A review of fuzzy AHP methods for decision-making with subjective judgements. Expert Syst. Appl. 2020, 161, 113738. [Google Scholar] [CrossRef]
  34. 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]
  35. Kamangar, M.; Katorani, S.; Tekyekhah, J.; Sohrabnejad, C.; Haderi, F.G. A novel hybrid MCDM model to select a suitable location for implement groundwater recharge *. Plant Arch. 2019, 19, 87–98. [Google Scholar]
  36. Kharazi, P.; Yazdani, M.R.; Khazealpour, P. Suitable identification of underground dam locations, using decision-making methods in a semi-arid region of Iranian Semnan Plain. Groundw. Sustain. Dev. 2019, 9, 100240. [Google Scholar] [CrossRef]
  37. Phankamolsil, Y.; Rittima, A.; Rantasewee, S.; Talaluxmana, Y.; Surakit, K.; Tabucanon, A.S.; Sawangphol, W.; Kraisangka, J. Analysis of Potential Site for Managed Aquifer Recharge Scheme in the Upper Greater Mae Klong Irrigation Project, Thailand. Appl. Environ. Res. 2022, 44, 80–94. [Google Scholar] [CrossRef]
  38. Tahvili, Z.; Khosravi, H.; Malekian, A.; Sigaroodi, S.K.; 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, 10. [Google Scholar] [CrossRef]
  39. 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, 4822. [Google Scholar] [CrossRef] [PubMed]
  40. Mahmoud, S.H.; Tang, X. Monitoring prospective sites for rainwater harvesting and stormwater management in the United Kingdom using a GIS-based decision support system. Environ. Earth Sci. 2015, 73, 8621–8638. [Google Scholar] [CrossRef]
  41. Kadhem, G.M.; Zubari, W.K. Identifying Optimal Locations for Artificial Groundwater Recharge by Rainfall in the Kingdom of Bahrain. Earth Syst. Environ. 2020, 4, 551–566. [Google Scholar] [CrossRef]
  42. Mati, B.; De Bock, T.; Malesu, M.; Khaka, E.; Oduor, A.; Meshack, M.; Oduor, V. Mapping the Potential of Rainwater Harvesting Technologies in Africa: A GIS Overview on Development Domains for the Continent and Ten Selected Countries; Technical Manual No. 6; World Agroforestry Centre (ICRAF): Nairobi, Kenya, 2006. [Google Scholar]
  43. Dabar, O.A.; Camberlin, P.; Pohl, B.; Waberi, M.M.; Awaleh, M.O.; Silah-Eddine, S. Spatial and temporal variability of rainfall over the Republic of Djibouti from 1946 to 2017. Int. J. Clim. 2021, 41, 2729–2748. [Google Scholar] [CrossRef]
  44. Awaleh, M.O.; Baudron, P.; Soubaneh, Y.D.; Boschetti, T.; Hoch, F.B.; Egueh, N.M.; Mohamed, J.; Dabar, O.A.; Masse-Dufresne, J.; Gassani, J. Recharge, groundwater flow pattern and contamination processes in an arid volcanic area: Insights from isotopic and geochemical tracers (Bara aquifer system, Republic of Djibouti). J. Geochem. Explor. 2017, 175, 82–98. [Google Scholar] [CrossRef]
  45. Jalludin, M.; Razack, M. Assessment of hydraulic properties of sedimentary and volcanic aquifer systems under arid conditions in the Republic of Djibouti (Horn of Africa). Hydrogeol. J. 2004, 12, 159–170. [Google Scholar] [CrossRef]
  46. MAEPE-RH, “Elaboration du programme national d’approvisionnement en eau potable et d’assainissement en milieu rural a l’horizon 2030, Djibouti,” Ministère de l’agriculture, de l’eau, de la pêche, de l’élevage et des ressources halieutiques (MAEPE-RH), 2015. Available online: https://www.afdb.org/fileadmin/uploads/afdb/Documents/Project-and-Operations/Djibouti - Activités dans le secteur de l’eau potable et de l’assainissement - mai 2013.pdf (accessed on 28 May 2023).
  47. NTT DATA Corp. AW3D Standard Product Guide: Product Specification Document. 2016. Available online: https://www.engesat.com.br/wp-content/uploads/AW3D-st-product-guide_201702.pdf (accessed on 28 May 2023).
  48. FAO; IIASA. Harmonized World Soil Database Version 2.0; FAO: Rome, Italy; IIASA: Laxenburg, Austria, 2023. [Google Scholar] [CrossRef]
  49. Jaafar, H.; Ahmad, F. GCN250, global curve number datasets for hydrologic modeling and design. Figshare Dataset 2019, 6, 145. [Google Scholar] [CrossRef]
  50. Dabar, O.A.; Adan, A.-B.I.; Ahmed, M.M.; Awaleh, M.O.; Waberi, M.M.; Camberlin, P.; Pohl, B.; Mohamed, J. Evolution and Trends of Meteorological Drought and Wet Events over the Republic of Djibouti from 1961 to 2021. Climate 2022, 10, 148. [Google Scholar] [CrossRef]
  51. Wu, Q. geemap: A Python package for interactive mapping with Google Earth Engine. J. Open Source Softw. 2020, 5, 2305. [Google Scholar] [CrossRef]
  52. QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation Project. 2022. Available online: https://qgis.org/en/site/ (accessed on 28 May 2023).
  53. Fuentes, I.; Vervoort, R.W. Site suitability and water availability for a managed aquifer recharge project in the Namoi basin, Australia. J. Hydrol. Reg. Stud. 2019, 27, 100657. [Google Scholar] [CrossRef]
  54. Ghazavi, R.; Babaei, S.; Erfanian, M. Recharge Wells Site Selection for Artificial Groundwater Recharge in an Urban Area Using Fuzzy Logic Technique. Water Resour. Manag. 2018, 32, 3821–3834. [Google Scholar] [CrossRef]
  55. Hussaini, M.S.; Farahmand, A.; Shrestha, S.; Neupane, S.; Abrunhosa, M. Site selection for managed aquifer recharge in the city of Kabul, Afghanistan, using a multi-criteria decision analysis and geographic information system. Hydrogeol. J. 2021, 30, 59–78. [Google Scholar] [CrossRef]
  56. European Space Agency. SENTINEL-2 User Handbook; European Space Agency: Cologne, Germany, 2015. [Google Scholar]
  57. 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 approach. Arab. J. Geosci. 2020, 13, 24. [Google Scholar] [CrossRef]
  58. Amineh, Z.B.A.; Hashemian, S.J.A.-D.; Magholi, A. Integrating Spatial Multi Criteria Decision Making (SMCDM) with Geographic Information Systems (GIS) for delineation of the most suitable areas for aquifer storage and recovery (ASR). J. Hydrol. 2017, 551, 577–595. [Google Scholar] [CrossRef]
  59. Sallwey, J.; Schlick, R.; Valverde, J.P.B.; Junghanns, R.; López, F.V.; Stefan, C. Suitability Mapping for Managed Aquifer Recharge: Development of Web-Tools. Water 2019, 11, 2254. [Google Scholar] [CrossRef] [Green Version]
  60. Ansems, N.; Visser, R.; Oord, A.; Mwango, F. Regional Mapping of the Potential of Managed Aquifer Recharge: A landscape- Based Approach. 2015. Available online: https://www.researchgate.net/publication/342390244 (accessed on 28 May 2023).
  61. Arshad, M.; Sarah, S.; Chatterjee, A.; Venkatarao, A.; Ahmed, S. Integrated approach to delineate sites for implementation of managed aquifer recharge (MAR) structures in fluoridated crystalline aquifer of south India. J. Earth Syst. Sci. 2022, 131, 67. [Google Scholar] [CrossRef]
  62. Selvarani, A.G.; Maheswaran, G.; Elangovan, K. Identification of Artificial Recharge Sites for Noyyal River Basin Using GIS and Remote Sensing. J. Indian Soc. Remote Sens. 2016, 45, 67–77. [Google Scholar] [CrossRef]
  63. Campello, B.S.C.; Duarte, L.T.; Romano, J.M.T. Dealing with multi-criteria decision analysis in time-evolving approach using a probabilistic prediction method. Eng. Appl. Artif. Intell. 2022, 116, 105462. [Google Scholar] [CrossRef]
  64. Yalcin, A.S.; Kilic, H.S.; Delen, D. The use of multi-criteria decision-making methods in business analytics: A comprehensive literature review. Technol. Forecast. Soc. Chang. 2021, 174, 121193. [Google Scholar] [CrossRef]
  65. Anastasiadou, K.; Gavanas, N. Enhancing urban public space through appropriate sustainable mobility policies. A multi-criteria analysis approach. Land Use Policy 2023, 132, 106765. [Google Scholar] [CrossRef]
  66. Koc, K.; Ekmekcioğlu, Ö.; Özger, M. An integrated framework for the comprehensive evaluation of low impact development strategies. J. Environ. Manag. 2021, 294, 113023. [Google Scholar] [CrossRef] [PubMed]
  67. Zhao, Y.; Zhou, Y. Identification of the critical hospitals in the urban post-disaster healthcare system based on the network modeling and multi-criteria decision-making. Int. J. Disaster Risk Reduct. 2023, 93, 103795. [Google Scholar] [CrossRef]
  68. Saaty, T. A scaling method for priorities in hierarchical structures. J. Math. Psychol. 1977, 15, 234–281. [Google Scholar] [CrossRef]
  69. Saaty, T.L. How to make a decision: The analytic hierarchy process. Eur. J. Oper. Res. 1990, 48, 9–26. [Google Scholar] [CrossRef]
  70. Goodarzi, M.R.; Niknam, A.R.R.; Barzkar, A.; Niazkar, M.; Mehrjerdi, Y.Z.; Abedi, M.J.; Pour, M.H. Water Quality Index Estimations Using Machine Learning Algorithms: A Case Study of Yazd-Ardakan Plain, Iran. Water 2023, 15, 1876. [Google Scholar] [CrossRef]
  71. Mallick, J.; Khan, R.A.; Ahmed, M.; Alqadhi, S.D.; Alsubih, M.; Falqi, I.; Hasan, M.A. Modeling Groundwater Potential Zone in a Semi-Arid Region of Aseer Using Fuzzy-AHP and Geoinformation Techniques. Water 2019, 11, 2656. [Google Scholar] [CrossRef] [Green Version]
  72. Zadeh, L.A. Fuzzy sets. Inf. Control 1965, 8, 338–353. [Google Scholar] [CrossRef] [Green Version]
  73. Ma, L.; Xu, Y.; Ngo, I.; Wang, Y.; Zhai, J.; Hou, L. Prediction of Water-Blocking Capability of Water-Seepage-Resistance Strata Based on AHP-Fuzzy Comprehensive Evaluation Method—A Case Study. Water 2022, 14, 2517. [Google Scholar] [CrossRef]
  74. Saaty, T.L. Decision making—the Analytic Hierarchy and Network Processes (AHP/ANP). J. Syst. Sci. Syst. Eng. 2004, 13, 1–35. [Google Scholar] [CrossRef]
  75. Chang, D.-Y. Applications of the extent analysis method on fuzzy AHP. Eur. J. Oper. Res. 1996, 95, 649–655. [Google Scholar] [CrossRef]
  76. Ekmekcioğlu, Ö.; Koc, K.; Özger, M. District based flood risk assessment in Istanbul using fuzzy analytical hierarchy process. Stoch. Environ. Res. Risk Assess. 2020, 35, 617–637. [Google Scholar] [CrossRef]
  77. Ishizaka, A.; Labib, A. Review of the main developments in the analytic hierarchy process. Expert Syst. Appl. 2011, 38, 14336–14345. [Google Scholar] [CrossRef] [Green Version]
  78. Hwang, C.I.; Yoon, K. Multiple Attribute Decision Making; Springer: Berlin, Germany, 1981. [Google Scholar]
  79. Albulescu, A.-C.; Minea, I.; Boicu, D.; Larion, D. Comparative Multi-Criteria Assessment of Hydrological Vulnerability—Case Study: Drainage Basins in the Northeast Region of Romania. Water 2022, 14, 1302. [Google Scholar] [CrossRef]
  80. Chen, C.-H. A New Multi-Criteria Assessment Model Combining GRA Techniques with Intuitionistic Fuzzy Entropy-Based TOPSIS Method for Sustainable Building Materials Supplier Selection. Sustainability 2019, 11, 2265. [Google Scholar] [CrossRef] [Green Version]
  81. Nguyen, H.X.; Nguyen, A.T.; Ngo, A.T.; Phan, V.T.; Nguyen, T.D.; Do, V.T.; Dao, D.C.; Dang, D.T.; Nguyen, A.T.; Nguyen, T.K.; et al. A Hybrid Approach Using GIS-Based Fuzzy AHP–TOPSIS Assessing Flood Hazards along the South-Central Coast of Vietnam. Appl. Sci. 2020, 10, 7142. [Google Scholar] [CrossRef]
  82. Koc, K.; Ekmekcioğlu, Ö.; Işık, Z. Developing a probabilistic decision-making model for reinforced sustainable supplier selection. Int. J. Prod. Econ. 2023, 259, 108820. [Google Scholar] [CrossRef]
  83. Ekmekcioğlu, Ö.; Koc, K.; Dabanli, I.; Deniz, A. Prioritizing urban water scarcity mitigation strategies based on hybrid multi-criteria decision approach under fuzzy environment. Sustain. Cities Soc. 2022, 87, 104195. [Google Scholar] [CrossRef]
  84. Nasta, P.; Adane, Z.; Lock, N.; Houston, A.; Gates, J.B. Links between episodic groundwater recharge rates and rainfall events classified according to stratiform-convective storm scoring: A plot-scale study in eastern Nebraska. Agric. For. Meteorol. 2018, 259, 154–161. [Google Scholar] [CrossRef]
  85. Tashie, A.M.; Mirus, B.B.; Pavelsky, T.M. Identifying long-term empirical relationships between storm characteristics and episodic groundwater recharge. Water Resour. Res. 2016, 52, 21–35. [Google Scholar] [CrossRef] [Green Version]
  86. Dunne, T.; Zhang, W.; Aubry, B.F. Effects of Rainfall, Vegetation, and Microtopography on Infiltration and Runoff. Water Resour. Res. 1991, 27, 2271–2285. [Google Scholar] [CrossRef]
  87. Ajjur, S.B.; Mogheir, Y.K. Identification of intrinsic suitable sites in Gaza Strip for the application of artificial groundwater recharge using a geographic information system multicriteria decision analysis. J. Multi-Criteria Decis. Anal. 2019, 27, 255–265. [Google Scholar] [CrossRef]
  88. Jourgholami, M.; Karami, S.; Tavankar, F.; Monaco, A.L.; Picchio, R. Effects of Slope Gradient on Runoff and Sediment Yield on Machine-Induced Compacted Soil in Temperate Forests. Forests 2020, 12, 49. [Google Scholar] [CrossRef]
  89. Martos-Rosillo, S.; Ruiz-Constán, A.; González-Ramón, A.; Mediavilla, R.; Martín-Civantos, J.; Martínez-Moreno, F.; Jódar, J.; Marín-Lechado, C.; Medialdea, A.; Galindo-Zaldívar, J.; et al. The oldest managed aquifer recharge system in Europe: New insights from the Espino recharge channel (Sierra Nevada, southern Spain). J. Hydrol. 2019, 578, 124047. [Google Scholar] [CrossRef]
  90. Fathi, S.; Hagen, J.S.; Haidari, A.H. Synthesizing existing frameworks to identify the potential for Managed Aquifer Recharge in a karstic and semi-arid region using GIS Multi Criteria Decision Analysis. Groundw. Sustain. Dev. 2020, 11, 100390. [Google Scholar] [CrossRef]
  91. Ben Meftah, M.; Mossa, M. New Approach to Predicting Local Scour Downstream of Grade-Control Structure. J. Hydraul. Eng. 2020, 146, 04019058. [Google Scholar] [CrossRef]
  92. Kallali, H.; Anane, M.; Jellali, S.; Tarhouni, J. GIS-based multi-criteria analysis for potential wastewater aquifer recharge sites. Desalination 2007, 215, 111–119. [Google Scholar] [CrossRef]
  93. Helbich, M. Spatiotemporal Contextual Uncertainties in Green Space Exposure Measures: Exploring a Time Series of the Normalized Difference Vegetation Indices. Int. J. Environ. Res. Public Heal. 2019, 16, 852. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Farswan, S.; Vishwakarma, C.A.; Mina, U.; Kumar, V.; Mukherjee, S. Assessment of rainwater harvesting sites in a part of North-West Delhi, India using geomatic tools. Environ. Earth Sci. 2019, 78, 329. [Google Scholar] [CrossRef]
  95. Rohde, M.M.; Stella, J.C.; Roberts, D.A.; Singer, M.B. Groundwater dependence of riparian woodlands and the disrupting effect of anthropogenically altered streamflow. Proc. Natl. Acad. Sci. USA 2021, 118, e20264531182021. [Google Scholar] [CrossRef]
  96. Ekmekcioğlu, Ö.; Koc, K. Explainable step-wise binary classification for the susceptibility assessment of geo-hydrological hazards. Catena 2022, 216, 106379. [Google Scholar] [CrossRef]
  97. Rahman, M.A.; Rusteberg, B.; Uddin, M.S.; Lutz, A.; Abu Saada, M.; Sauter, M. An integrated study of spatial multicriteria analysis and mathematical modelling for managed aquifer recharge site suitability mapping and site ranking at Northern Gaza coastal aquifer. J. Environ. Manag. 2013, 124, 25–39. [Google Scholar] [CrossRef] [PubMed]
  98. Soliman, K.; Sallam, O.M.; Schüth, C. Delineating MAR Sites Using GIS-MCDA for Nuweiba Alluvial Fan Aquifer, Sinai, Egypt. Water 2022, 14, 475. [Google Scholar] [CrossRef]
  99. Itani, N.; Harik, G.; Alameddine, I.; El-Fadel, M. Managed aquifer recharge in karstic systems:Site suitability mapping by coupling multi-criteria decision analysis with remote sensing and hydrologic modeling. J. Environ. Manag. 2022, 322, 116162. [Google Scholar] [CrossRef]
  100. Ezzeldin, M.; Konstantinovich, S.E.; Igorevich, G.I. Determining the suitability of rainwater harvesting for the achievement of sustainable development goals in Wadi Watir, Egypt using GIS techniques. J. Environ. Manag. 2022, 313, 114990. [Google Scholar] [CrossRef]
  101. Więckowski, J.; Kizielewicz, B.; Shekhovtsov, A.; Sałabun, W. RANCOM: A novel approach to identifying criteria relevance based on inaccuracy expert judgments. Eng. Appl. Artif. Intell. 2023, 122, 106114. [Google Scholar] [CrossRef]
  102. INOWAS. T06. MAR Method Selection-INOWAS. 2018. Available online: https://inowas.com/t06-mar-method-selection/ (accessed on 6 January 2023).
  103. ONEAD. Station d’Epuration de Douda. 2023. Available online: https://www.onead.dj/station-depuration-de-douda/ (accessed on 28 May 2023).
Figure 1. Flow of the research steps.
Figure 1. Flow of the research steps.
Water 15 02534 g001
Figure 2. Study area.
Figure 2. Study area.
Water 15 02534 g002
Figure 3. Thematic maps of the relevant decision criteria for MAR potential mapping: (a) slope, (b) soil texture, (c) curve number, (d) rainfall, (e) normalized difference vegetation index, (f) drainage density, (g) depth to groundwater, (h) geology, and (i) electrical conductivity.
Figure 3. Thematic maps of the relevant decision criteria for MAR potential mapping: (a) slope, (b) soil texture, (c) curve number, (d) rainfall, (e) normalized difference vegetation index, (f) drainage density, (g) depth to groundwater, (h) geology, and (i) electrical conductivity.
Water 15 02534 g003
Figure 4. Consistency ratios of the experts who completed the pairwise comparison questionnaire.
Figure 4. Consistency ratios of the experts who completed the pairwise comparison questionnaire.
Water 15 02534 g004
Figure 5. Sensitivity analysis with respect to different fuzziness degrees.
Figure 5. Sensitivity analysis with respect to different fuzziness degrees.
Water 15 02534 g005
Figure 6. MAR suitability potential of Djibouti based on the integrated FAHP-TOPSIS approach: (a) the closeness coefficient (CCi) distribution and (b) the reclassified version of it.
Figure 6. MAR suitability potential of Djibouti based on the integrated FAHP-TOPSIS approach: (a) the closeness coefficient (CCi) distribution and (b) the reclassified version of it.
Water 15 02534 g006
Figure 7. Percentage distributions of the different suitability classes in terms of surface area.
Figure 7. Percentage distributions of the different suitability classes in terms of surface area.
Water 15 02534 g007
Figure 8. Graphical representation of the feasible MAR techniques in Djibouti.
Figure 8. Graphical representation of the feasible MAR techniques in Djibouti.
Water 15 02534 g008
Figure 9. Practical utilization of the constructed decision-making scheme.
Figure 9. Practical utilization of the constructed decision-making scheme.
Water 15 02534 g009
Table 2. Linguistic scales and triangular fuzzy reciprocals of AHP and fuzzy AHP.
Table 2. Linguistic scales and triangular fuzzy reciprocals of AHP and fuzzy AHP.
Linguistic VariablesAHPFuzzy AHP
ImportanceValue for Reciprocals Triangular   Fuzzy   Numbers   l i j , m i j , u i j Triangular   Fuzzy   Reciprocals   1 / u i j , 1 / m i j , 1 / l i j
Equally important1(1/1)(1,1,1)(1,1,1)
Intermediate value2(1/2)(1,2,3)(1/3,1/2,1)
Moderately important3(1/3)(2,3,4)(1/4,1/3,1/2)
Intermediate value4(1/4)(3,4,5)(1/5,1/4,1/3)
Important5(1/5)(4,5,6)(1/6,1/5,1/4)
Intermediate value6(1/6)(5,6,7)(1/7,1/6,1/5)
Very important7(1/7)(6,7,8)(1/8,1/7,1/6)
Intermediate value8(1/8)(7,8,9)(1/9,1/8,1/7)
Extremely important9(1/9)(9,9,9)(1/9,1/9,1/9)
Table 3. Random index.
Table 3. Random index.
n 12345678910
Random Index000.520.891.111.251.351.401.451.49
Table 4. Profiles of the experts who performed pairwise comparisons.
Table 4. Profiles of the experts who performed pairwise comparisons.
IDSectorJob DescriptionBackgroundExperience (Years)
Expert 1AcademiaProfessorCivil and Environmental Engineering16
Expert 2AcademiaAssociate ProfessorEnvironmental Engineering11
Expert 3MunicipalityHead of DepartmentCivil Engineering (PhD)20
Expert 4MunicipalityPlanning EngineerArchitecture (MSc)6
Expert 5Water administrationUnit ManagerCivil Engineering8
Expert 6Water administrationTechnical Office EngineerGeological Engineering (MSc)5
Expert 7Private sectorGeneral ManagerEnvironmental Engineering (MSc)15
Expert 8Private sectorModeling and Design EngineerGeological Engineering6
Table 5. Decision criteria weights and their respective ranking.
Table 5. Decision criteria weights and their respective ranking.
ClusterWeightCriteriaWeightRank
LocalGlobalLocalGlobal
Surface30.36%Slope43.92%13.33%12
Soil texture24.81%7.53%36
Curve number31.27%9.50%24
Environment48.85%Rainfall62.80%30.68%11
Normalized difference vegetation index24.33%11.89%23
Drainage density12.86%6.28%38
Subsurface20.79%Depth to groundwater29.32%6.10%39
Geology33.41%6.95%27
Groundwater quality (EC)37.27%7.75%15
Table 6. Sensitivity analysis results based on shifting criteria weights.
Table 6. Sensitivity analysis results based on shifting criteria weights.
Scenario IDShiftRMSEScenario IDShiftRMSEScenario IDShiftRMSE
Scenario 1S–ST0.0209Scenario 13ST–DG0.0062Scenario 25Rl–G0.1932
Scenario 2S–CN0.0259Scenario 14ST–G0.0014Scenario 26Rl–EC0.1519
Scenario 3S–Rl0.1416Scenario 15ST–EC0.0008Scenario 27NDVI–DD0.0242
Scenario 4S–NDVI0.0075Scenario 16CN–Rl0.2059Scenario 28NDVI–DG0.0327
Scenario 5S–DD0.0366Scenario 17CN–NDVI0.0190Scenario 29NDVI–G0.0213
Scenario 6S–DG0.0384Scenario 18CN–DD0.0130Scenario 30NDVI–EC0.0182
Scenario 7S–G0.0239Scenario 19CN–DG0.0095Scenario 31DD–DG0.0005
Scenario 8S–EC0.0246Scenario 20CN–G0.0154Scenario 32DD–G0.0028
Scenario 9ST–CN0.0126Scenario 21CN–EC0.0067Scenario 33DD–EC0.0037
Scenario 10ST–Rl0.1687Scenario 22Rl–NDVI0.1234Scenario 34DG–G0.0036
Scenario 11ST–NDVI0.0141Scenario 23Rl–DD0.1633Scenario 35DG–EC0.0046
Scenario 12ST–DD0.0049Scenario 24Rl–DG0.1938Scenario 36G–EC0.0030
Table 7. Comprehensive evaluation of the present study compared with past attempts.
Table 7. Comprehensive evaluation of the present study compared with past attempts.
ReferenceCountryScaleProblemIIIIIIIVVVI
Kazakis [25]GreeceWatershedSaltwater intrusion10AHPWLC
Fuentes and Vervoort [53]AustraliaWatershedWater table decline9AHPWLC
Sandoval and Tiburan [27]PhilippinesWatershedGroundwater depletion10AHPWOA7
Kharazi et al. [36]IranWatershedWater scarcity16N/AAHP, TOPSIS, and EDAS7
Itani et al. [99]LebanonWatershedSaltwater intrusion9AHPWLC4
Hussaini et al. [55]AfghanistanCityWater table decline7AHP and ANPFL and WOA
Papadopoulos et al. [10]GreeceWatershedExcess water storage9Fuzzy AHPFIS
Zhang et al. [31]South AfricaWatershedWater scarcity12AHPWLC
Arshad et al. [61]IndiaWatershedChemical contamination7AHPWOA
Shadmehri Toosi et al. [32]IranWatershedWater scarcity6AHPWLC
Ezzeldin et al. [100]EgyptWatershedWater scarcity11AHPWLC
This StudyDjiboutiCountryWater scarcity and saltwater intrusion9Fuzzy AHPTOPSIS8
Notes: I: Number of decision criteria; II: criteria weighting technique; III: alternative prioritization techniques; IV: experts’ details; V: number of experts; VI: sensitivity analysis; WOA: weighted overlay analysis; FL: fuzzy logic; FIS: fuzzy inference system.
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

Mouhoumed, R.M.; Ekmekcioğlu, Ö.; Başakın, E.E.; Özger, M. Integrated Fuzzy AHP-TOPSIS Model for Assessing Managed Aquifer Recharge Potential in a Hot Dry Region: A Case Study of Djibouti at a Country Scale. Water 2023, 15, 2534. https://doi.org/10.3390/w15142534

AMA Style

Mouhoumed RM, Ekmekcioğlu Ö, Başakın EE, Özger M. Integrated Fuzzy AHP-TOPSIS Model for Assessing Managed Aquifer Recharge Potential in a Hot Dry Region: A Case Study of Djibouti at a Country Scale. Water. 2023; 15(14):2534. https://doi.org/10.3390/w15142534

Chicago/Turabian Style

Mouhoumed, Rachid Mohamed, Ömer Ekmekcioğlu, Eyyup Ensar Başakın, and Mehmet Özger. 2023. "Integrated Fuzzy AHP-TOPSIS Model for Assessing Managed Aquifer Recharge Potential in a Hot Dry Region: A Case Study of Djibouti at a Country Scale" Water 15, no. 14: 2534. https://doi.org/10.3390/w15142534

APA Style

Mouhoumed, R. M., Ekmekcioğlu, Ö., Başakın, E. E., & Özger, M. (2023). Integrated Fuzzy AHP-TOPSIS Model for Assessing Managed Aquifer Recharge Potential in a Hot Dry Region: A Case Study of Djibouti at a Country Scale. Water, 15(14), 2534. https://doi.org/10.3390/w15142534

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