Next Article in Journal
Modeling Wave Overtopping on a Seawall with XBeach, IH2VOF, and Mase Formulas
Next Article in Special Issue
Integrated Technology for Evaluation and Assessment of Multi-Scale Hydrological Systems in Managing Nonpoint Source Pollution
Previous Article in Journal
Effects of Combined High Temperature and Waterlogging Stress at Booting Stage on Root Anatomy of Rice (Oryza sativa L.)
Previous Article in Special Issue
Sea Level Rise Effect on Groundwater Rise and Stormwater Retention Pond Reliability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Analytical Hierarchy Process and Multi-Influencing Factors to Map Groundwater Recharge Zones in a Semi-Arid Mediterranean Coastal Aquifer

1
LR01ES06 Laboratory of Geological Resources and Environment, Department of Geology, Faculty of Sciences of Tunis, University of Tunis El Manar, Tunis–El Manar 2092, Tunisia
2
Department of Biosystems and Agricultural Engineering, Oklahoma State University, 111 Agricultural Hall, Stillwater, OK 74078, USA
3
Department of Geology, Faculty of Sciences of Tunis, University of Tunis El Manar, Tunis–El Manar 2092, Tunisia
4
National Institute of Agronomy, University of Carthage, Tunis 1082, Tunisia
5
AGHYLE, SFR Condorcet FR CNRS 3417, Polytechnic Institute UniLaSalle Beauvais, 19 Pierre Waguet Street, 60026 Beauvais, France
6
HSM, Univ. Montpellier, IRD, CNRS, 34090 Montpellier, France
7
Georesources Laboratory, Water Research and Technology Centre, Borj Cedria Ecopark 8, P.O. Box 273, Soliman 8020, Tunisia
*
Author to whom correspondence should be addressed.
Water 2020, 12(9), 2525; https://doi.org/10.3390/w12092525
Submission received: 28 July 2020 / Revised: 1 September 2020 / Accepted: 1 September 2020 / Published: 10 September 2020

Abstract

:
Mapping groundwater recharge zones (GWRZs) is essential for planning artificial recharge programs to mitigate groundwater decline and saltwater intrusion into coastal aquifers. We applied two multi-criteria decision-making approaches, namely the analytical hierarchy process (AHP) and the multi-influencing factors (MIF), to map GWRZs in the Korba aquifer in northeastern Tunisia. GWRZ results from the AHP indicate that the majority (69%) of the area can be classified as very good and good for groundwater recharge. The MIF results suggest larger (80.7%) very good and good GWRZs. The GWRZ maps improve groundwater balance calculations by providing estimates of recharge-precipitation ratios to quantify percolation. Lithology, land use/cover and slope were the most sensitive parameters followed by geomorphology, lineament density, rainfall, drainage density and soil type. The AHP approach produced relatively more accurate results than the MIF technique based on correlation of the obtained GWRZs with groundwater well discharge data from 20 wells across the study area. The accuracy of the approaches ultimately depends on the classification criteria, mean rating score and weights assigned to the thematic layers. Nonetheless, the GWRZ maps suggest that there is ample opportunity to implement aquifer recharge programs to reduce groundwater stress in the Korba aquifer.

1. Introduction

Excessive groundwater pumping is a global problem, which particularly affects arid/semi-arid regions where surface water is scarce [1,2,3,4]. Demand for fresh water resources is noticeably increasing as a result of growing agricultural water demand, population growth and rapid industrialization. Therefore, groundwater abstraction has become an integral component of many water management plans worldwide, including in rural areas of the Western Mediterranean region [4,5,6,7,8].
The Korba coastal area in northeastern Tunisia is faced with groundwater over-pumping, causing seawater intrusion and degradation of groundwater quality. In this region, groundwater is the primary water source for domestic and agricultural purposes for about 90% of the rural population [2,5,6,7,8,9] and 50% of urban residents [6,9,10], putting pressure on the aquifer to cope with the rising demand [2,5,9,10]. Groundwater overdraft was identified as a problem in the early 1980s when the Tunisian government called for a reduction of water pumping by encouraging farmers to adopt drip irrigation systems by covering up to 60% of the installation cost [11]. Since then, groundwater levels have continued to decline and seawater intrusion is steadily progressing in most of the coastal areas with the expansion of agricultural lands despite efforts to implement water conservation policies, build recharge dams, transfer water from the north, stop drilling new wells and encourage the re-use of treated wastewater [4,8]. Most of the deficit is observed in the central part of the Korba area where 53% of croplands are irrigated [2,5,6,7,11].
Delineation of groundwater recharge zones (GWRZ) is essential for planning artificial recharge programs to mitigate groundwater decline [12,13,14,15]. Groundwater recharge is affected by human activities and natural possesses on the ground surface and within the aquifer system. Many factors affect the occurrence and movement of groundwater in a region, including rainfall, lineament density, slope, drainage density, land use/land cover, lithology, geomorphology and soil types [12,14,16,17,18,19].
Geographic information system (GIS) and Remote Sensing (RS) techniques facilitate estimating surface and subsurface water over large areas [15,18,20,21,22,23]. A variety of methods have been implemented for groundwater recharge zone mapping worldwide [24,25,26,27] such as frequency ratio (FR), certainty factor (CF), weights-of-evidence (WOE), fuzzy logic index models, analytical hierarchy process (AHP) and multi-influencing factors (MIF). In particular, the AHP and MIF techniques have gained attention as useful approaches for groundwater prediction for rapid, precise and cost-effective assessment of groundwater recharge potential [28,29,30,31,32,33] by significantly reducing the mathematical complexity of decision-making based on systematic expert judgment. AHP is a GIS-based multi-criteria decision-making (MCDM) technique for pairwise comparison of spatial parameters by assigning weights based on expert opinion [34,35,36,37]. The MIF is an MCDM technique that computes spatial relationships between a dependent variable and independent variables according to scores that are assigned based on major and minor influential factors affecting the GWRZ [38,39]. These two approaches have gained popularity because they are relatively simple and useful for practical applications prior to launching costly field explorations [19,20,23].
The objective of this paper is to use RS data and multi-criteria analysis using AHP and MIF techniques to delineate GWRZs in the Korba aquifer and evaluate the results of these approaches based on field observations. The accuracy of the AHP and MIF techniques was evaluated using a receiver operating characteristic (ROC) curve based on the area under the ROC curve (AUC), which includes all the possible decision thresholds from field observations of well yield data and water level depth. Although many investigations have assessed the groundwater resources of the Korba area [2,5,6,7,8,10,11] none of them have used remotely derived data to demarcate GWRZs in this region. The GWRZ maps facilitate groundwater management through better water balance calculations based on improved estimates of recharge-precipitation ratios to quantify percolation.

2. Study Area

The Korba aquifer is part of the Cap-Bon Peninsula located off the Mediterranean Coast of northeastern Tunisia (Figure 1). The main geomorphological features of the study area are shallow flood plain and shallow buried pediplain. The spatial extent of the surface of the shallow flood plain is 56 km2 (12.76% of the study area), which enhances the recharge potential. This coastal area is dominated by shallow buried pediplain (denudational origin) which covers an area of 290.5 km2 (66.32%). This feature has Good groundwater storage potential and it is located around rivers and the beach, which was filled with alluvial material deposited by wadis and its effluents.
The region is one of the most productive agricultural areas in Tunisia, contributing about 15% of the country’s total agricultural production [9]. Major crops grown in the study area are vegetables (potatoes, tomatoes, strawberries, etc.), vineyard plantations and citrus [8,11]. The region is characterized by a semi-arid climate (average annual rainfall: 420 mm) with Mediterranean influence. Average monthly high temperature occurs in July (37 °C) and the lowest average monthly low is 6.1 °C in January. There are ephemeral streams known as wadis, which originate from the Sidi Abderrahmen anticline and receive seasonal runoff. The most important wadis in the region (i.e., Chiba, M’laabi and Lebna) are obstructed by upstream dams. Consequently, groundwater has become the primary water resource in downstream areas.
Three main geological formations constitute the Korba aquifer system (Figure 2a) [2,5,6,10,40]:
  • The Tyrrhenian Quaternary made of arenitic limestones overlying conglomeratic units with a total thickness ranging 10–50 m, which forms a 1.2 km wide strip parallel to the Mediterranean coastline;
  • The Pliocene age corresponds to marine sediments deposited in the Dakhla syncline in the Northern part of the Korba City. The subsurface lithological details collected from some private wells show that vertical variations consist of yellow sands with alternating clay and sandstone layers;
  • Upper Miocene (the Somâa sands) laying on Medium Miocene (sandstones and marls of the Saouaf Formation), which constitutes the base of the aquifer. This formation contains thick fine sand layers of continental origin, including conglomeratic layers.
The main compressive structures that bend in the northeast-southwest direction (Figure 2a) were generated during the Atlasique phase of the Upper Miocene [5,6,10,40]. These structures include: (1) Djebel Sidi Abderrahmene anticlinal, an elongated antiform parallel to the atlasic folds along the northeast–southwest direction; (2) Dakha synclinal, which contains the Korba coastal aquifer; (3) Menzel Temime anticlinal, which stretches along the Northeast-Southwest direction, beginning in the north of Nabel due to seismicity and ending at Korba. This anticlinal is influenced by inverse faults; and (4) Diare el Hajjej anticlinal and Taffeloun synclinal.
A distensible regime dominated the last phase of the Upper Miocene and a majority of Pliocene causing subsidence of the entire South-West region (Nabeul and Hammamet). This subsidence continued during the Upper Pliocene and Quaternary when the sinking trenches were formed. Normal faults developed during Tortonian synsedimentary in the Somâa Formation along the northwest-southeast and west/northwest-east/southeast directions. These lineaments dip 40–60° in the southeast direction.
The Korba unconfined aquifer system, assigned to the Plio-Quaternary unit, occupies the coastal plain but might also be locally semi-confined due to less permeable lenses of large lateral extent (Figure 2b). The saturated thickness of the aquifer varies from 3–10 m in the east toward the coastal area to a maximum of 50 m in the west near Somâa village with a Miocene marl Formation. The western part of the aquifer may contain sandstone lenses with different thicknesses and lateral variations, creating a third Miocene reservoir with considerable water resources [2,10,11].
The unconfined aquifer system receives infiltration from rainfall, that is, the main natural recharge source. The recharge is estimated to be more than 10% of the 420 mm/year average annual rainfall but also from lateral inflow from the surrounding relief, dunes and wadis during the rainy season [7,10,11]. Other important sources of recharge are irrigation return flow and lateral recharge by leakage (or through faults) from the underlying Miocene sandstone, which are not yet evaluated [2].
Since 1960s, groundwater exploitation in the Cap-Bon area has increased from 270 wells pumping 4 mm3 in 1962 [10] to more than 13,000 wells pumping 68 mm3 in 2014 [11]. There is a congestion of wells in the coastal part of the plain, which draws landward the freshwater-saltwater interface [8]. The spatial distribution of the piezometric heads measured in 1962, 1977, 2014 and 2018 are presented in Figure 3. The piezometric maps of 1962 and 1977 (Figure 2b and Figure 3a) show a smooth potentiometric surface located above sea level before the major perturbations of surface water regime. The growth of groundwater abstraction during the 1970s caused the hydraulic heads in the Korba aquifer to approach zero piezometric level. The groundwater table decline formed a concentric depression in the region of Diar El-Hojjaj, mainly due to intensive exploitation coupled with the reduction of recharge by damming the wadis (i.e., Chiba (1963) and M’laabi (1964)).
The piezometric map of 2014 (Figure 3c) shows a groundwater cone of depression reaching 12 m below sea level observed 3 Km from the Mediterranean shoreline. In addition, a localized piezometric depression of up to 5 m below sea level is observed in the western part of the study area [8]. The latter reveals (Figure 3d) a general water table decline (at least 10 m) during the 56-year period, creating local cones of depression that increase in size from year to year. The greatest piezometric depression (up to 12 m below sea level) is observed in the western part, several kilometers from the Mediterranean shoreline in a region that is prone to vertical saltwater up-coning. Nearer to the shoreline, the piezometric level contours become increasingly negative, reversing the hydraulic gradients toward the central part of the aquifer, which accelerates saltwater intrusion from the Mediterranean Sea.

3. Material and Methods

We identified suitable sites for groundwater recharge through a knowledge-based analysis of eight factors, namely rainfall, lineament density, slope, drainage density, land use/land cover, lithology, geomorphology and soil type layers [41,42,43,44]. These factors were examined independently for groundwater recharge zoning to estimate recharge rates. Figure 4 illustrates the methodology used to delineate suitable sites for groundwater recharge.

3.1. Preparation of Input Database

A range of field data and remotely sensed data were collected from various government agencies to map suitable artificial recharge areas. Table 1 summarizes the data types, sources and application. Primary data include digital elevation model (DEM) available from the Shuttle Radar Topography Mission (SRTM) and Sentinel-2A (S2A) multispectral imagery. Rainfall maps were generated by linear kriging interpolation method on a Surfer platform using a long-term average point rainfall for six pluviometry stations in the Cap-Bon over 37 years (1980–2017) [45].
The slope factor was calculated using the 30 m DEM based on the maximum rate of change in value from each cell to neighboring cells [46]. The data related to drainage density were generated indirectly from the slope data. Drainage density (Dd) is defined as the closeness of spacing of stream channels [47] estimated as the measure of the total length of the stream segment of all orders per unit area using Equation (1) [48]. The drainage frequency value of the area ranges from 0 to >4 and the map indicates high moderate, medium and low drainage frequency.
D d =   i = 1 n ( D i A ) ,
where, ΣDi is the total length of all streams in stream order i (km) and A is the watershed area (km2).
The lineament density of the Korba region was prepared from a combination of Sentinel-2A (S2A) multispectral imagery, DEM and black and white photographs using computer-aided supervised extraction. To achieve optimum interpretation with reference to frequency, the “sliding window” method was applied and the number of lineaments in each unit area or pixel was counted [49,50]. The thematic layer for lineament density (Ld) can be defined as the total length of all recorded lineaments divided by the area of the catchment under consideration as follows [51]:
L d =   i = 1 n ( L i A ) ,
where, Li is the length of the ith lineament, ΣLi is the total length of all lineaments in km and A is the area of the grid (km2).
Lineament raster was created using false color composite in spatial analysis tool in ArcGIS 10.1 followed by additional editing of the boundary divide line and road features to ensure the quality of the extracted lineaments. The frequency value ranges from 0 to over 2 km/km2, providing a basis for evaluating the lineament frequency and their distribution. A lineament frequency contour map was prepared along with a synoptic classification indicating very low, low, moderate, high and very high frequency value.
The spatial soil groups and geomorphological map were derived from the soil map of La Carte Agricole of Nabeul governorate at 1:50,000 for the year 2005, which is an official spatial data source. The vector layer was converted to a grid and then reclassified according to recharge weighted rating [51,52]. The land-use/land cover (LU/LC) raster (30 m resolution) was obtained by processing satellite images from the Landsat 8 Sensor on a scale of 1/500,000 by ENVI software. The lithological map for the Korba area was obtained using a subsurface geology map and several well drilling profiles of the Cap-Bon region [53]. The Korba area is mainly covered by six types of lithological units that were ranked based on available porosity and permeability values of aquifer rocks. Sedimentary features, such as sand and sandstone, generally demonstrate very high porosity (20–31%) and permeability (500–7000 mD) and they are capable of accumulating and storing large volumes of groundwater [54,55,56].

3.2. Analytic Hierarchy Process (AHP)

The AHP [57,58] facilitates decision making by organizing perceptions and judgments into a multi-level hierarchical structure that accounts for the forces that control a decision through assessing multiple factors [59]. The AHP is used as a complex multi-criteria decision-solving method to determine the weights assigned to different thematic layers and their respective features. This technique breaks down the multi-criteria decision problem into a hierarchy based on a pair-wise comparison of the importance of different criteria and sub-criteria within the judgment matrix [60,61,62]. The hierarchy allows analysis to consider each of several properties separately [54] for the selection of potential recharge zones from competing sets of parameters. The method was implemented in four steps to delineate potential recharge zones in the Korba study area: (1) selecting factors that influence groundwater recharge zones (2) developing pairwise comparison matrix, (3) estimating relative weights and (4) assessing matrix consistency.

3.2.1. Selection of Factors Influencing Groundwater Recharge Zones

In the first step of the AHP, each factor that influences recharge was given a score between 1 and 9, depending on its significance compared to other factors in pairwise comparisons [58]. For this, a standard Saaty’s 1–9 scale was used (Table 2) to describe the relative influence of parameters, where score 1 denotes equal influence of parameters and score 9 denotes extreme influence of a parameter on groundwater recharge compared to the other parameters [58].

3.2.2. Pairwise Comparison Matrix

The AHP method integrates and transforms spatial data (input) into decision (output), where qualitative information of individual thematic layers and features are converted into quantitative scores based on Saaty’s scale. Then, a pairwise comparison matrix (PCM) [19,26,29] is constructed (Equation (3)) using Saaty’s scores obtained in the previous step. In the PCM, the matrix column is constructed based on a descending order of parameter influence on recharge. The first element is assigned a score of 1 when compared to itself (see Table 3). Other elements of the rows are filled using the actual Saaty’s scores when a more influential parameter is compared with a less influential parameter or the reciprocal of the Saaty’s scores score when a less influential parameter is compared to a more influential parameter.
A = [ X 11 X 12   X 1 n X 21 X 22   X 2 n       X n 1 X n 2   X nn ] ,
where A is a pairwise comparison matrix where element Xnn denotes relative significance of a parameter for recharge compared to another parameter.
Table 3 provides the PCM for the parameters examined in this study. Lithology was selected as the first parameter of the matrix because has a higher influence on recharge potential compared to the other factors. Thus, lithology was assigned the value 8. Land use/land cover was selected as the second most important parameter influencing recharge followed by slope, geomorphology, lineaments, rainfall, drainage and finally soil parameter in a descending order of influence. Each parameter in the selected set was assigned a Saaty’s score based on its influence on recharge relative to lithology.

3.2.3. Estimating Relative Weights

Weights were assigned to the variables based on ‘expert’ opinion to estimate the relative importance of variables compared to other variables and to quantify the relative influence of each variable on recharge [26,60,61]. The layers were assigned weights derived by normalizing the pair comparison matrix (NPCM) [62,63]. The NPCM elements were computed by dividing thematic element values by their corresponding total column values from the PCM (Equation (4)) (see Table 4):
X i j = C i j L i j ,
where Xij is normalized pair-wise matrix value at ith row and jth column, Cij is the value assigned to each criteria at ith row and jth column and Lij is the total values in each column of the pair-wise matrix.
Subsequently, a standard weight was calculated for variable i by dividing each normalized pairwise matrix elements by criterion number (N) (Equation (5)) [64,65].
W i =   X i j N ,
where Wi is standard weight.
Then, Eigen vector and eigenvalue calculations help determine the percentage of effect of the thematic layers and the classification of the constraints (Table 5). The eigenvector was calculated by dividing column elements by the column sum in Table 4. The principal Eigen vector was obtained by averaging across the rows to quantify relative weights of each parameter [57,59,66]. A consistency vector was obtained by multiplying two different matrix values from selected thematic layers (Equation (6)), namely, pair-wise comparison matrix value and normalized pair-wise matrix value [3,4,67].
    λ =   ( C i j X i j )   ,
where, λ is the consistency vector.
The sum of eigenvalues called principal eigenvalue (λmax) is a measure of matrix deviation from consistency [66]. According to Saaty [57], for a pairwise comparison matrix to be consistent it must have a principal eigenvalue (λmax) greater than or equal to the number of the parameters considered (n). The principal eigenvalue of 8.62 was obtained for the 8 × 8 matrix (see Table 5), which was used for the calculation of consistency index.

3.2.4. Assessing Matrix Consistency

Consistency is checked based on Consistency Index (CI) and Consistency Ratio (CR) [48,50,53,55,56]:
    C I = λ m a x   n n 1  
C R = C I R I ,
where, CI: consistency index, λmax: highest matrix Eigenvalue, n: number of variables (thematic layers), CR: consistency ratio and RI: Random Index value based on the number of variables.
A perfectly consistent decision maker should always yield CI = 0 but small values of inconsistency may be tolerated if the CI < 0.1 [57,68]. We obtained an acceptable CI value of 0.08. Also, if the CR is greater than 0.1, the pairwise comparison judgments must be re-evaluated. With a matrix of eight variables, the RI is 1.41 (see Table 6). The applied weighting yielded a CR of 0.05, which shows that the weights (Table 5) assigned to GIS thematic layers of parameters are consistent [58,59,60].

3.3. Multi-Influencing Factors (MIF)

MIF is a widely used MCDM technique for environmental management [12,34,36]. The method is useful for mapping groundwater recharge sites by assigning an appropriate weight to various factors and feature classes based on the influence on groundwater flow and storage (Figure 5).
The ranks and weights are assigned to each factor and different classes based on their relative contribution to groundwater recharge. We specified major and minor influential relationships among the factors (Figure 5). Each relationship was weighted according to its influence on recharge. A large weight indicates that the factor has a large influence on groundwater recharge. The representative weight of a factor influencing recharge potential is the sum of all weights from each factor. Major effects were assigned a weight of 1.0 whereas minor influential relationships were given a weight of 0.5 [55,56,57]. Furthermore, the classes with no effect on groundwater recharge were assigned a weight of zero. The weights for each factor were summed to obtain the recharge potential factor weight [34,35]. The estimated weight for each influencing factor was obtained as a percentage using Equation (9):
S c o r e = ( M a j o r   e f f e c t + M i n o r   e f f e c t ) ( M a j o r   e f f e c t + M i n o r   e f f e c t ) × 100 .
The results are shown in Table 7.
Figure 6 illustrates the GIS analysis procedure to map the spatial distribution of groundwater recharge using the weighted overlay tool of Weighted Overlay analysis tool of ArcGIS 10.1 based on the eight thematic layers and their corresponding percentage influence on recharge [69,70]. Values in the input raster layers were reclassified into a common evaluation scale of 1 (very good), 2 (good), 3 (moderate), 4 (poor) and 5 (very poor). This was done by multiplying the cell values of each factor class by the factor weight and summing the resulting cell values to produce a map of potential recharge zones [71,72], as summarized in Equation (10):
GWRZ = i = 0 n W i × R i = ( LI c LI w +   LULC c LULC w +   SLO c SLO w +   GM c GM w + LD c LD w , +   RF c RF w +   DD c DD w +   SL c SL w ) ,
where GWRZ is the groundwater recharge zone, Wi is the weight of each thematic layer, Ri is the rating of each class of each thematic layer, LI is the lithology of the aquifer, LU/LC is the land-use/land-cover, SLO is the slope, GM is the geomorphology, LD is the lineament, RF is the rainfall, DD is the drainage density and SL is the soil cover. The subscripts c and w refer, respectively, to the factor class of a thematic layer and its percent influence on recharge [46,56]. It means that each layer class rank (each individual class has a rank) was multiplied by the weight from Saaty’s assumption of each layer (each layer has a unique weight) to obtain the position of layers in the overlay analysis to map GWRZ.

3.4. Sensitivity Analysis

The GWRZ index is sensitive to the factor weight and the rating value for each class of each layer [73,74]. Sensitivity analysis reveals how much each thematic layer and its weight influence the output map. A map removal sensitivity analysis (MRSA) and a single parameter sensitivity analysis (SPSA) were conducted to examine the effect of each of the eight factors’ weights on the final GWRZ.

3.4.1. Map Removal Sensitivity Analysis (MRSA)

MRSA was conducted to determine the most or the least influential thematic layer for groundwater recharge zones [30]. The sensitivity index, S, was computed using Equation (11):
S = [ | GWRZ N GWRZ n | GWRZ   ]   × 100 ,
where GWRZ and GWRZ′ are the output of groundwater recharge map index of all the thematic layers and when one of the thematic layers is removed, respectively. N is the number of the full thematic layers used to compute the GWRZ and n is the number of thematic layers used to compute GWRZ’.

3.4.2. Single Parameter Sensitivity Analysis (SPSA)

The single-parameter sensitivity analysis was performed to evaluate the effect of each factor in the groundwater recharge potential index. This sensitivity test was conducted because the numerical scores assigned to the parameters, which are the basis of parameter weightings, are essentially arbitrary [11]. The initial arbitrary weights are replaced with effective weights [22,75] and the sensitivity index was calculated according to the Equation (12):
S = [ W i R i GWRZ ]   × 100 ,
where; Wi is the weight for each layer, Ri is the rating value for each class of each layer, GWRZ is the groundwater recharge zone.

3.5. Verification/Validation of Groundwater Recharge Zone

An important step for any model is validation to assess the performance by establishing the relationship between groundwater recharge level class and GWRZ map [76,77,78,79]. Several methods can be used to validate groundwater recharge maps, such as receiver performance analysis, curve area, groundwater yield estimation of wells during field visits and a comparative study between the water level profile and groundwater recharge zones, among others [78,80,81,82]. In this study, two methods were used to validate the obtained groundwater recharge zones, namely: Receiver Operating Characteristics (ROC) curve and groundwater level monitoring and assessment at 25 shallow wells. The groundwater yield data were collected during several field investigations or obtained from local management agencies [83].
The ROC curve was applied as a mathematical technique to determine the map accuracy [78,80,81]. Using SPSS-statistical software, ROC curves were obtained by considering cumulative percentage of area (on the x axis) and the cumulative percentage of number of wells (on the y axis) [84]. The area under the curve (AUC) was calculated to quantify the ability to predict the occurrence or non-occurrence of pre-defined ‘‘events’’ [83,84,85,86]. The resultant AUC values range from 0.5 (i.e., a random prediction) to 1 (i.e., an excellent prediction) [81]. The quantitative–qualitative relationship between the AUC value and prediction accuracy can be excellent (0.9–1), very good (0.8–0.9), good (0.7–0.8), average (0.6–0.7) and poor (0.5–0.6) [85]. To illustrate the extent to which the occurrence or non-occurrence of pre-defined events was correctly predicted, the AUC method was verified through frequency percentage for AHP and MIF techniques. A comparative study was carried out between the water level profile and groundwater recharge zones through two water level profiles for two sections in the study area AA’ and BB’ (see Figure 8), for AHP and MIF methods, respectively.

4. Results and Discussion

4.1. Evaluation of Factors Governing Groundwater Recharge Zones

Table 8 reports the final weights of the factors governing GWRZs, which are described below.
Lithology (LI): LI governs the porosity and permeability of aquifer rocks [18], which in turn influences the occurrence and distribution of groundwater recharge through physio-mechanical properties that control the ability of the aquifer materials to transmit water and the rate at which groundwater flows [27]. The lithology of the study area consists of quaternary alluvium (12% of the surface area; 52.56 km2), sandy clay (22%; 96.36 km2), yellow sand with clay (6%; 26.28 km2), tyrrhenian sandstones (3%; 13.14 km2), sandy sandstone (3%; 13.14 km2) and silty sand (54%; 236.52 km2) with ratings varying from 1 to 6 (Table 8 and Figure 7a). The quaternary alluvial sediments are identified along the wadis across the study area. The lithology is predominantly silty sand of Dakhla syncline and sandy clay with sandstone intercalations of Pliocene Formation. Laterally, the facies changes southward to more differentiated yellow sand with clay, sand and sandstone layers. The higher ratings (6 and 5), apply to 60% of the area and are assigned to the sediments of sandstone, gravel and sand located at the “Sable of the Somâa Formation” and in the central part of the study area around Diarr El Hojjaj and Menzel Horr villages. These delineated zones are characterized by high permeability index, which implies that rainwater can easily infiltrate [11], providing high hydraulic conductivities (e.g., ranging from 10−6 to 10−3 m/s) while hydraulic conductivities at the Somâa Formation are smaller [10]. The Tyrrhenian sandstones form the most permeable unit, followed by the central part of the Pliocene (alluvial deposits), the Pliocene, the Somâa sand and finally, the early Miocene Formations [2].
Land use/Land cover (LU/LC): LU/LC has a major influence on the occurrence and development of groundwater in a terrain [26,87]. Agricultural lands typically allow more infiltration due to pore spaces in the soil, which trap and hold the water in the roots, providing a pathway for water to percolate into the surface by loosening up the rock and soil. By contrast, built-up and barren lands reduce infiltration due to loss of permeable surface and increased runoff potential. Therefore, the areas with agriculture and water bodies are considered as good sites for groundwater recharge, while settlements and barren lands have poor groundwater recharge potential [73]. A major portion of the study area is under agriculture (69.42%; 304.06 km2) followed by barren land (16.58%; 72.62 km2), rural/urban settlements (8.65%; 37.89 km2) and river (5.35%; 24.43 km2) (Figure 7b).
Slope (SLO): SLO is a characteristic of local and regional relief, which is an important factor because it influences the water retention, intensity of infiltration [49,52], aquifer recharge and groundwater movement. The higher the ground slope, the smaller the infiltration rate will be due to larger runoff potential [12,27,58]. The study area has five slope classes (Figure 7c), namely very low (0–5%), low (5–10%), moderate (10–15%), high (15–20%) and very high (>20%). The topography layer displayed a gentle slope (0–5%) over most of the study area (87%; 381.06 km2) which has been assigned a weight score of 6 with Very Good potential for artificial recharge. The slope increases from east to west due to the presence of the Abderrahman Mountain range. Slopes (5–10%) are considered Good for groundwater storage. The third category (i.e., Poor) includes two slope ranges of 10–15% and 15–20%. Gently sloping Very Poor recharge areas cover 8.76 km2 (2%) of the study area, while the area under steep slope is negligible. A higher priority was given to a nearly level to gentle slope categories followed by moderate and steep slopes (Table 8).
Geomorphology (GM): GM of an area reflects various landforms, which includes their description, species and physical processes that help assess groundwater recharge potential and evaluate possible groundwater areas [20,23,25,50]. Geomorphic features of the study area were categorized into five units namely: shallow flood plain and beach, shallow buried pediplain, moderately buried pedipalin, deep buried pediplain and sedimentary high ground (Figure 7d).
Moderately buried pedipalin covered 9.01% (39.5 km2) of the study area. It comprises gently undulating plains covered with sand layers of continental origin including conglomeratic layers and clay lenses that are moderately favorable for groundwater recharge. Deep buried pediplain and sedimentary high ground cover 4.32% (18.92 km2) and 7.59% (33.24 km2) of the study area, respectively and are characterized by highly sloping topography, smaller amounts of infiltration and high surface runoff characteristic of a poor recharge zone.
Lineament (LD): LD density networks such as fractures, joints and faults increase porosity and therefore play an important role in groundwater movement and high groundwater recharge potential [26,27,33]. The LD density, extracted from satellite images and remote sensing, is an indicator of the degree of fracturing of the rock, meaning that groundwater recharge potential is high near lineaments. While high lineaments frequency indicates very high recharge potential due to the presence of recharge pathways, low frequency does not necessarily translate into very low recharge potential. In this study, the lineament density was subdivided and ranked into five classes: Lineament with densities ranging from 0.0–5.0 km/km2 dominated the study area (57.22%; 250.62 km2). Lineament density ranging from 0.5–1 km/km2 and higher than 1 km/km2 covers an area of 45.10 km2 (32.09%) and 33.78 km2 (10.69%), respectively. The rose diagram shows the dominance of NE–SW and NW-SE directions of lineaments (Figure 7e). These conjugate lineament directions are characteristic of Tunisian Atlas domain and they increase reservoir permeability through creating interconnected lineaments/fractures.
Rainfall (RN): RN is typically the primary source of groundwater recharge in which the water infiltrates into subsurface through fractures and soils. It governs the amount of runoff that would be available to capture in recharge basins to increase infiltration. The mean annual rainfall ranges from 400 to 450 mm in the lowland areas along the coastline, while it ranges from 450 and 500 mm in the mountainous areas. An increase of annual rainfall with altitude is detected. The resulting rainfall map was classified into two major classes (Table 8 and Figure 7f): 450–500 mm/yr and 400–450 mm/yr. Generally, 75% of the 420 mm/yr average annual rainfall occurs during the wet season (from September to March), while summers are usually dry. From the map of rainfall and piezometric levels (see Figure 3), higher-altitude areas have greater recharge potential than lower altitudes. In this study, the rainfall factor was assigned a weight of 0.8 and 0.5 for MIF and AHP, respectively, in the total groundwater recharge potential.
Drainage density (DD): DD is another factor that affects the movement of water and groundwater recharge [44,48,55]. It is a measure of how a watershed is drained by stream channels. DD is influenced by numerous factors, including resistance to erosion, infiltration capacity, vegetation cover, surface roughness and runoff intensity index and climatic conditions [16,18,19]. Areas with high drainage density have less recharge rate, whereas low drainage density areas have a high recharge rate and can directly influence the groundwater recharge. Five drainage density categories were identified in the study area (Figure 7g), namely ‘very good’ (3–4 km/km2 or more), ‘good’ (2–3 km/km2), ‘moderate’ (1–2 km/km2) and ‘poor’ (0.1–1 km/km2). High drainage was recorded in the high relief near to Somâa village in the southwestern part. A higher priority was given to a low drainage density category followed by medium and high drainage density categories.
Soil (SL): soil texture controls the percolation and infiltration of surface water into the aquifers, influencing groundwater recharge through properties such as porosity, structure, adhesion and consistency [26,44]. Based on soil texture data, the study area has three main soil types, namely gravel and sand, loam-sandy loam and clay-clay loam (Figure 7h). Sandy loam is the dominant soil textural class in the study area (69.72%; 305.37 km2). Gravel and sandy soils cover a major portion of the northern and central parts (27.18%; 119.05 km2). Nearly 3% of the clay loam soil partially covers north eastern and southwestern portions of the study area. The ranks of soils were assigned according to their degree of infiltration [21,22,46]. The gravel and sandy soil has a high degree of infiltration and, therefore, has a higher influence, while the clay and clay loam soil has the lowest degree of infiltration and, therefore, is given the lowest influence (Table 8).

4.2. Groundwater Recharge Maps

4.2.1. Analytical Hierarchy Process

The identification of GWRZ was accomplished based on the rates and weights of the eight thematic layers according to Equation (13):
GWRZ = 0.33   ×   L I + 0.23   × L U L C + 0.16   × S L O + 0.11   × G M + 0.07   × L D + 0.05   × R N + 0.03   × D D + 0.02   × S L .
Results of the GWRZ mapping are shown in Figure 8a and Figure 9. The Korba area has 70 km2 of very good groundwater recharge zones (about 16% of the area). The GWRZ map indicates the high groundwater recharge potential zone situated mostly near the shallow flood pediplain zone. This is mainly attributed to the availability of the high quaternary alluvial sediment in plain terrain along with the rivers course where there are high to very high lineament density (1–>2 km/km2) with gravel and sand soil type. The GWRZ designated as good is mostly found in the north and the central part of the Korba area, which covers 53% (123 km2). The good potential zone is distributed across irrigated agricultural lands (e.g., vineyards, citrus, strawberries, potatoes, tomatoes, etc.) and grassland on gentle slopes (0–5%) and flat buried pediplains with sandy loam soil type. The moderate zone covers 11.5% (50.37 km2) of the area, scattered along the north-west direction where the lineament density is high (>2 km/km2). The poor 15.8% (69.2 km2) and very poor 2.9% (0.37 km2) groundwater recharge areas are found especially in the south-eastern portion of the study area. This may be due to the fact that the area exhibits high sedimentary ground of mountainous geomorphology, relatively less water bearing lithological Formations (Somâa sand with clay) and steep topography (>10% slope category).

4.2.2. Multi-Influencing Factors Techniques (MIF)

The GWRZ index using the MIF technique was calculated using Equation (14):
GWRZ = 0.22   ×   L I + 0.19   × L U L C + 0.13   × S L O + 0.13   × G M   + 0.11   × D L   + 0.08   × R N + 0.08   × D D + 0.06 × S L .
Very high GWRZ (134.47 km2 or 30.7%) is identified especially in the center of the study area, along the Mediterranean coastline and along the Chiba and Lebna river channels (Figure 8b and Figure 9). The geomorphology of this region consists of shallow pediplain with shallow quaternary alluvium deposit, which allows water to infiltrate. However, the river channels increase the drainage density, routing the surface runoff downstream, which reduces groundwater recharge. Good GWRZ is about 219 km2 (50%) covering the north, central and southwestern parts. The lithological Formations with the highest groundwater recharge potential are the alluvial deposits and silt-sand. Furthermore, the good GWRZ is attributed to flat buried pediplains with sandy loam soil type and the relatively high amounts of rainfall on low slopes. Moderate GWRZ (42.05 km2 or 9.6%) is mainly concentrated in the western part and scattered patches throughout the study area. The low potential zones (poor and very poor) are identified in the north-western and southwestern parts, respectively, due to higher slope and unfavorable geology (Miocene age deposition) and geomorphological condition (high ground) where hard rock terrains are present. The areas covered by poor and very poor GWRZ are about 35.04 Km2 (8.07%) and 6.92 km2 (1.58%), respectively.

4.3. Results of the Sensitivity Analysis

4.3.1. Map Removal Sensitivity Analysis (MRSA)

In Map Removal Sensitivity Analysis (MRSA) for AHP and MIF methods, the thematic layers were removed from the mapping process one at a time, using the remaining layers to map GWRZ. The computed sensitivity index (S) values are reported in Table 9, indicating that the most sensitive layers for GWRZ map are LI, with the mean variations index (SMVI) of 1.62% and 1.54%, LU/LC, with SMVI of 1.24% and 1.26% and SLO with SMVI of 1.04% and 1.01%, for AHP and MIF, respectively. This is essentially due to the relatively high theoretical weight assigned to these layers (Table 8). In addition to the relatively high weight of LI (33% for AHP and 22% for MIF), this layer has high rating values in almost every sub-area and it is thus considered the most influential thematic layer for mapping the GWRZ in the Korba area.
Likewise, the removal of LU/LC (AHP weight: 23%; MIF weight: 19%) and SLO (AHP weight: 16%; MIF weight: 13%) has a great impact on the GWRZ map.
The fourth and fifth layers are GM (SMVI for AHP: 0.74%; SMVI for MIF: 0.71%) and LD (SMVI for AHP: 0.64%; SMVI for MIF: 0.69%). The least influential factors for the output map determination were the RN thematic map (SMVI = 0.48% and 0.29% for AHP and MIF, respectively) and DD thematic map (SMVI = 0.31% and 0.29% for AHP and MIF, respectively).
The lowest MRSA values, in AHP and MIF, were 0.27% and 0.21%, respectively, for the SL layer. Thus, the sensitivity order of the thematic layers for the AHP method is LI > LU/LC > SLO > GM > LD > RN > DD > SL, which is consistent with the magnitudes of the theoretical weights in the AHP method. However, the sensitivity order of variables is slightly different in MIF technique (LI > LU/LC > SLO > GM > LD > RN = DD > SL). The variation of the sensitivity index is governed by the collective impact of the large data range, the mean rating score and especially the weight assigned to the thematic layer. Overall, the MRSA indicates that almost all the factors are necessary for GWRZ mapping in the Korba aquifer (Table 9).

4.3.2. Single Parameter Sensitivity Analysis (SPSA)

Single Parameter Sensitivity Analysis (SPSA) analysis confirms that most of the study area has good to very good groundwater recharge potential and all the parameters have a significant impact on recharge. The analysis was performed to better understand the contribution of each factor to the groundwater recharge map [88]. In this section, the “effective” weight of each factor was compared with the theoretical weight assigned to it by the AHP and MIF methods (Table 10) in order to determine the “most effective impact” factors. As the factor’s effective weight increases, the higher its sensitivity compared to other factors with a lower effective weight.
For AHP, groundwater recharge tends to be most sensitive to LU/LC, LI, SLO and GM layers (Table 10) because their mean effective weights, 86%, 66%, 44% and 39%, respectively, are higher than their corresponding theoretical weights, 23%, 33%, 16% and 11% (see Table 7). On the other hand, LD, RF, DD and SL are the least effective, with a mean effective weight of 25%, 25%, 10% and 8% respectively, which are still higher than their corresponding theoretical weights 7%, 5%, 3% and 2%. For MIF, too, the results (Table 10), indicate that LU/LC is the most effective factor, with a mean effective weight of 71%, compared to its theoretical weight (19%). The GM (46%), LI (44%), RN (40%), LD (39%) and SLO (36%) also show a higher effective weight compared to their theoretical weights, 13%, 22%, 8%, 11% and 13%, respectively. The rest of the Factors such as DD (28%) and SL (24%) exhibit less effective weights than others.

4.4. Verification/Validation of Groundwater Recharge Zones

The hydrogeological data in combination with the hydrological data provide new information to advance understanding of key hydrological processes controlling groundwater flow and recharge in the shallow Korba aquifer. Lithology is the most important component in determining groundwater recharge potential in the study area due to the nature of the geological formations and how they influence infiltration rates in the surficial materials. Pumping rates are expected to be correlated with parameters such as transmissivity (high transmissivity allows high discharge) or depth of the groundwater table (shallow depth facilitates groundwater abstraction). Ten years of with 10 years of measured well discharge data from 25 monitoring wells were used to cross-validate the GWRZs. Locations of the monitoring wells were overlain on the map of groundwater recharge potential to evaluate the results from the AHP and MIF techniques (Figure 8a,b). The maximum and minimum depths of the monitoring wells are 28 m and 5 m, respectively and the average depth is 9 m. The evaluation approach is deemed adequate given the generally shallow depths of groundwater wells.
The accuracy of the AHP method to map groundwater recharge potential in the Korba area was about 73% whereas the MIF method’s accuracy was 67%, indicating good agreement between the two approaches. In relation to topography, the good and very good groundwater recharge zones coincide with shallow and moderate water level depth locations, which can be explained by the low thickness of the unsaturated zone (ranging from 3.5 to 10 m) composed of thick fine sand layers of continental origin. Areas with small depth to groundwater table are located where there is high infiltration. The good recharge zone is mostly found in wet regions located in the central part of the study area with intensive irrigated cultivations on highly permeable soil where the number of wells is proportional to the number of farms. These nearly flat alluvial areas with large hydraulic conductivities (ranging from 10−6 to 10−3 m/s) and flood pediplain have good groundwater recharge potential. The poor and very poor recharge zones are found where water level is deep to very deep (Figure 8). These areas exhibit diverse geomorphological variations (numerous structural hills) where the western most areas are characterized by very steep escarpment and low soil porosity, resulting in low recharge potential in the upland plain.
Eight very high discharge wells (20–30 m3/h) were located in the very good recharge zones obtained from both AHP and MIF methods. These wellfields are located primarily in the sandstone and quaternary alluvium in the central part of the study area, which has the highest transmissivity zones of about 10−2 m2/s. Three of five high discharge (10–20 m3/h) wells fell within good recharge zones and two were in moderate zones. Further, three medium discharge (5–10 m3/h) wells were located in the moderate recharge zone except one that fell within the good zone according to the MIF generated map. For AHP method, five of six low discharge (1–5 m3/h) wells exist in the poor zone and one well in the moderate zone. For MIF technique, two of the six low discharge wells occur in the moderate recharge zone. All three wells with very low discharge were found in the very poor zone except one well in moderate zone. Finally, the three very low discharge wells (<1 m3/h) fell within very poor recharge zones except one that was located in the moderate zone based on the MIF method.
The groundwater recharge maps were also validated based on AUC. The area under the ROC curve (AUC) is an indicator of model quality. Generally, a good model has an AUC value of 0.7–0.9, while an excellent model has values over 0.9. According to Figure 10, the AHP method outperformed MIF (AHP AUC: 75.6%; MIF AUC: 70.4%), although both methods were reasonably accurate for mapping groundwater recharge in the Korba area. Overall, the validation results increase confidence in the applied methodology as a useful framework for rapid assessment of groundwater recharge to inform siting artificial recharge structures and other groundwater management efforts.

4.5. Groundwater Recharge Estimation

Quantifying groundwater recharge rate, which is difficult to measure directly [66,70,72,75], improves understanding of watershed scale hydrologic processes. For the Korba case study, groundwater recharge has been estimated to be up to 10% of the 420 mm/year mean average rainfall [7,10]. Ennabli [10] provides an excellent mathematical description of the problems related to artificial recharge value from precipitation (ranging between 15% and 25%). Furthermore, a 3D numerical model allowed Kerrou et al. [2] to estimate the areal recharge in the range of 8% and 30% of the annual rainfall (33.6–126 mm/yr), which is enough to counteract the intrusion of about 7 mm3/yr of seawater into the aquifer. However, groundwater overdraft reduced the submarine groundwater discharge of 16 mm3/yr, causing saltwater front to advance inland. Based on the works of Zghibi et al. [8,11], the precipitation recharge contributes 43.5% of the inflow to the aquifer system, depending especially on the frequency and the timing and intensity of the rain events, soil type and geological conditions.
The GWRZ maps were used to derive recharge-precipitation ratios for the two calculation methods (Table 11). A simple estimation of water volume (W) in the subsurface media was performed for the AHP and MIF methods by calculating the volume of annual precipitation (P) × recharge ratio ×% of area, where P = 175.56 106 m3/yr (420 mm/yr) and S (surface) = 438 km2:
  • For AHP method:
    W AHP   =   175.75   ×   10 6   [ ( 0.475   ×   0.16 )   +   ( 0.325   ×   0.53 )   +   ( 0.15   ×   0.115 ) +   ( 0.075   ×   0.0158 )   +   ( 0.025   ×   0.029 ) ]   =   48.69   10 6   m 3 / yr
  • For MIF technique:
    W MIF   = 175.75 × 10 6   [ ( 0.475 × 0.307 ) + ( 0.325 × 0.50 ) + ( 0.15 × 0.096 ) + ( 0.075 × 0.0807 ) + ( 0.025 × 0.0158 ) ] = 57.60   10 6   m 3 / yr .
Thus, based on the results of the AHP method 27.73% (116.5 mm/yr) of the annual rainfall recharges the Korba aquifer whereas the MIF method indicates that recharge is about 32.60% (137.8 mm/yr) of the annual rainfall. These estimated recharge values can be further verified using tracer techniques and water-table fluctuations method to inform groundwater conservation planning.

5. Conclusions

Two MCDM methods, namely analytical hierarchy process (AHP) and multi-influencing factor (MIF) were applied using a geospatial analysis framework to delineate groundwater recharge zones (GWRZ) in the Korba unconfined aquifer in northeastern Tunisia. The factors that were considered include lithology, land use/land cover, slope, geomorphology, lineament density, rainfall, drainage density and, soil type. The AHP method classified the GWRZs as very good (16%), good (53%), moderate (11.5%), poor (15.8%) and very poor (2.9%). The classification based on the MIF method suggests that 30.7% of the study area has very good recharge potential while 50% is considered good. The remainder of the Korba area has moderate (9.6%), poor (8.1%) and very poor potential (1.6%). The GWRZ maps were used to estimate recharge-precipitation ratios to quantify percolation (i.e., 27.7% (116.5 mm/yr) based on AHP and 32.6% (137.8 mm/yr) based on MIF). The groundwater recharge maps were verified using field yield data and water level depth from 25 shallow wells. While the accuracy of the GWRZs generated by both MCDM methods was satisfactory, the validation phase indicates that the AHP outperformed the MIF in the Korba region (AHP accuracy: 75.6% compared to MIF accuracy: 70.4%). The results are valuable for groundwater management and planning artificial recharge to mitigate groundwater tale decline and salt water intrusion in this coastal aquifer. The presented case study illustrates how groundwater recharge potential can be assessed using AHP and MIF to support sustainable groundwater management.

Author Contributions

All authors collaborated in the research presented in this publication by making the following contributions: research conceptualization, A.Z., M.H.M., A.M. (Amira Merzougui) and A.C.; methodology, A.Z., A.M. (Ali Mirchi), L.Z., J.-D.T., A.C. and I.C.; formal analysis, A.Z.; data curation, A.Z., M.H.M. and A.M. (Amira Merzougui); writing—original draft preparation, A.Z. and A.M. (Ali Mirchi); writing—review and editing, A.Z., A.M. (Ali Mirchi), L.Z., J.-D.T., A.C., I.C. and J.T.; supervision, A.M. (Ali Mirchi), L.Z., J.-D.T. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the United States Agency for International Development (USAID) through Partnerships for Enhanced Engagement in Research (PEER; Grant number 7-444).

Acknowledgments

The authors would like to thank the staff of the Regional Commissariat for Agricultural Development of Nabeul, its local representatives and the Tunisian Union of Agriculture and Fisheries for their support and interest in this work. We appreciate logistical support from the National Academies of Sciences, Engineering and Medicine (NASEM).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ashraf, B.; AghaKouchak, A.; Alizadeh, A.; Mousavi-Baygi, M.; Moftakhari, H.R.; Mirchi, A.; Anjileli, H.; Madani, K. Quantifying Anthropogenic Stress on Groundwater Resources. Sci. Rep. 2017, 7, 12910. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kerrou, J.; Renard, P.; Tarhouni, J. Status of the Korba groundwater resources (Tunisia): Observations and three-dimensional modelling of seawater intrusion. Hydrogeol. J. 2010, 18, 1173–1190. [Google Scholar] [CrossRef]
  3. Murmu, P.; Kumar, M.; Lal, D.; Sonker, I.; Singh, S.K. Delineation of groundwater potential zones using geospatial techniques and analytical hierarchy process in Dumka district, Jharkhand, India. Groundw. Sustain. Dev. 2019, 9, 100239. [Google Scholar] [CrossRef]
  4. Nithya, C.N.; Srinivas, Y.; Magesh, N.; Seenipandi, K.; Nair, N.C. Assessment of groundwater potential zones in Chittar basin, Southern India using GIS based AHP technique. Remote Sens. Appl. Soc. Environ. 2019, 15, 100248. [Google Scholar] [CrossRef]
  5. Chekirbane, A.; Tsujimura, M.; Lachaal, F.; Khadhar, S.; Mlayah, A.; Kawachi, A.; Isoda, H.; Tarhouni, J.; Benalaya, A. Quantification of Groundwater—Saline Surface Water Interaction in a Small Coastal Plain in North-East Tunisia using Multivariate Statistical Analysis and Geophysical Method. Water Environ. Res. 2016, 88, 2292–2308. [Google Scholar] [CrossRef]
  6. Kouzana, L.; Benassi, R. Geophysical and geochemical study of the seawater intrusion in Mediterranean semi arid zones. Case of the Korba coastal aquifer (Cap-Bon, Tunisia). J. Afr. Earth Sci. 2010, 58, 242–254. [Google Scholar]
  7. Paniconi, C.; Khlaifi, I.; Lecca, G.; Giacomelli, A.; Tarhouni, J. A modeling study of sea water intrusion in the Korba coastal plain, Tunisia. Phys. Chem. Earth 2001, 26, 345–351. [Google Scholar] [CrossRef]
  8. Zghibi, A.; Mirchi, A.; Zouhri, L.; Taupin, J.-D.; Chekirbane, A.; Tarhouni, J. Implications of groundwater development and seawater intrusion for sustainability of a Mediterranean coastal aquifer in Tunisia. Environ. Monit. Assess. 2019, 191, 696. [Google Scholar] [CrossRef]
  9. CCI Cap-Bon. Chambre de Commerce et d’Industrie du Cap-Bon, Tunisia. 2018. Available online: http://www.ccicapbon.org.tn (accessed on 24 March 2018).
  10. Ennabli, M. Etude Hydrogélogique des Aquifères du Nord-Est de la Tunisie Pour Une Gestion Intégrée des Ressources en Eau. Ph.D. Thesis, Nice University, Nice, France, 1980. [Google Scholar]
  11. Zghibi, A.; Merzougui, A.; Chenini, I.; Ergaieg, K.; Zouhri, L.; Tarhouni, J. Groundwater vulnerability analysis of Tunisian coastal aquifer: An application of DRASTIC index method in GIS environment. Groundw. Sustain. Dev. 2016, 2, 169–181. [Google Scholar] [CrossRef]
  12. Etikala, B.; Golla, V.; Li, P.; Renati, S. Deciphering groundwater potential zones using MIF technique and GIS: A study from Tirupati area, Chittoor District, Andhra Pradesh, India. HydroResearch 2019, 1, 1–7. [Google Scholar] [CrossRef]
  13. Oikonomidis, D.; Dimogianni, S.; Kazakis, N.; Voudouris, K. A GIS/Remote Sensing-based methodology for groundwater potentiality assessment in Tirnavos area, Greece. J. Hydrol. 2015, 525, 197–208. [Google Scholar] [CrossRef]
  14. Rajasekhar, M.; Raju, G.S.; Raju, R.S.; Basha, U.I. Data on artificial recharge sites identified by geospatial tools in semi-arid region of Anantapur District, Andhra Pradesh, India. Data Brief 2018, 19, 462–474. [Google Scholar] [CrossRef] [PubMed]
  15. Yousif, M. Hydrogeological inferences from remote sensing data and geoinformatic applications to assess the groundwater conditions: El-Kubanyia basin, Western Desert, Egypt. J. Afr. Earth Sci. 2019, 152, 197–214. [Google Scholar] [CrossRef]
  16. Mishra, A.K.; Upadhyay, A.; Srivastava, A.; Rai, S.C. Probabilistic groundwater recharge zonation in hard rock terrain using geospatial techniques in Veniar watershed, South India. Ecohydrol. Hydrobiol. 2019. [CrossRef]
  17. Mogaji, K.A.; Omobude, O.B. Modeling of geoelectric parameters for assessing groundwater potentiality in a multifaceted geologic terrain, Ipinsa Southwest, Nigeria—A GIS-based GODT approach. NRIAG J. Astron. Geophys. 2017, 6, 434–451. [Google Scholar] [CrossRef]
  18. Rahmati, O.; Melesse, A.M. Application of Dempster–Shafer theory, spatial analysis and remote sensing for groundwater potentiality and nitrate pollution analysis in the semi-arid region of Khuzestan, Iran. Sci. Total Environ. 2016, 568, 1110–1123. [Google Scholar] [CrossRef]
  19. 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]
  20. Kumar, V.A.; Mondal, N.C.; Ahmed, S. Identification of Groundwater Potential Zones Using RS, GIS and AHP Techniques: A Case Study in a Part of Deccan Volcanic Province (DVP), Maharashtra, India. J. Indian Soc. Remote Sens. 2020, 48, 497–511. [Google Scholar] [CrossRef]
  21. Nhu, V.-H.; Rahmati, O.; Falah, F.; Shojaei, S.; Al-Ansari, N.; Shahabi, H.; Shirzadi, A.; Górski, K.; Nguyen, H.; Ahmed, B. Mapping of Groundwater Spring Potential in Karst Aquifer System Using Novel Ensemble Bivariate and Multivariate Models. Water 2020, 12, 985. [Google Scholar] [CrossRef] [Green Version]
  22. Saidi, S.; Bouri, S.; Ben Dhia, H. Sensitivity analysis in groundwater vulnerability assessment based on GIS in the Mahdia-Ksour Essaf aquifer, Tunisia: A validation study. Hydrol. Sci. J. 2011, 56, 288–304. [Google Scholar] [CrossRef]
  23. Wu, Q.; Si, B.C.; He, H.; Wu, P. Determining Regional-Scale Groundwater Recharge with GRACE and GLDAS. Remote Sens. 2019, 11, 154. [Google Scholar] [CrossRef] [Green Version]
  24. Arshad, A.; Zhang, Z.; Zhang, W.; Dilawar, A. Mapping favorable groundwater potential recharge zones using a GIS-based analytical hierarchical process and probability frequency ratio model: A case study from an agro-urban region of Pakistan. Geosci. Front. 2020. [Google Scholar] [CrossRef]
  25. Bhattacharya, R.K.; Das Chatterjee, N.; Das, K. An integrated GIS approach to analyze the impact of land use change and land cover alteration on ground water potential level: A study in Kangsabati Basin, India. Groundw. Sustain. Dev. 2020, 11, 100399. [Google Scholar] [CrossRef]
  26. Lentswe, G.B.; Molwalefhe, L. Delineation of potential groundwater recharge zones using analytic hierarchy process-guided GIS in the semi-arid Motloutse watershed, eastern Botswana. J. Hydrol. Reg. Stud. 2020, 28, 100674. [Google Scholar] [CrossRef]
  27. Rajasekhar, M.; Raju, G.S.; Sreenivasulu, Y.; Raju, R.S. Delineation of groundwater potential zones in semi-arid region of Jilledubanderu river basin, Anantapur District, Andhra Pradesh, India using fuzzy logic, AHP and integrated fuzzy-AHP approaches. HydroResearch 2019, 2, 97–108. [Google Scholar] [CrossRef]
  28. Abijith, D.; Saravanan, S.; Singh, L.; Jennifer, J.J.; Saranya, T.; Parthasarathy, K. GIS-based multi-criteria analysis for identification of potential groundwater recharge zones—A case study from Ponnaniyaru watershed, Tamil Nadu, India. HydroResearch 2020, 3, 1–14. [Google Scholar] [CrossRef]
  29. Abrams, W.; Ghoneim, E.; Shew, R.; LaMaskin, T.; Al-Bloushi, K.; Hussein, S.; Abubakr, M.; Al-Mulla, E.; Al-Awar, M.; El-Baz, F. Delineation of groundwater potential (GWP) in the northern United Arab Emirates and Oman using geospatial technologies in conjunction with Simple Additive Weight (SAW), Analytical Hierarchy Process (AHP), and Probabilistic Frequency Ratio (PFR) techniques. J. Arid. Environ. 2018, 157, 77–96. [Google Scholar] [CrossRef]
  30. Arefin, R. Groundwater potential zone identification at Plio-Pleistocene elevated tract, Bangladesh: AHP-GIS and remote sensing approach. Groundw. Sustain. Dev. 2020, 10, 100340. [Google Scholar] [CrossRef]
  31. Das, B.; Pal, S.C. Combination of GIS and fuzzy-AHP for delineating groundwater recharge potential zones in the critical Goghat-II block of West Bengal, India. HydroResearch 2019, 2, 21–30. [Google Scholar] [CrossRef]
  32. Hammami, S.; Zouhri, L.; Souissi, D.; Souei, A.; Zghibi, A.; Marzougui, A.; Dlala, M. Application of the GIS based multi-criteria decision analysis and analytical hierarchy process (AHP) in the flood susceptibility mapping (Tunisia). Arab. J. Geosci. 2019, 12, 653. [Google Scholar] [CrossRef]
  33. 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]
  34. Magesh, N.; Chandrasekar, N.; Soundranayagam, J.P. Delineation of groundwater potential zones in Theni district, Tamil Nadu, using remote sensing, GIS and MIF techniques. Geosci. Front. 2012, 3, 189–196. [Google Scholar] [CrossRef] [Green Version]
  35. Thapa, R.; Gupta, S.; Guin, S.; Kaur, H. Assessment of groundwater potential zones using multi-influencing factor (MIF) and GIS: A case study from Birbhum district, West Bengal. Appl. Water Sci. 2017, 7, 4117–4131. [Google Scholar] [CrossRef]
  36. Al-Ruzouq, R.; Shanableh, A.; Merabtene, T.; Siddique, M.; Khalil, M.A.; Idris, A.; Almulla, E. Potential groundwater zone mapping based on geo-hydrological considerations and multi-criteria spatial analysis: North UAE. Catena 2019, 173, 511–524. [Google Scholar] [CrossRef]
  37. Jothibasu, A.; Anbazhagan, S. Modeling groundwater probability index in Ponnaiyar River basin of South India using analytic hierarchy process. Model. Earth Syst. Environ. 2016, 2, 109. [Google Scholar] [CrossRef] [Green Version]
  38. Meena, S.R.; Mishra, B.; Piralilou, S.T. A Hybrid Spatial Multi-Criteria Evaluation Method for Mapping Landslide Susceptible Areas in Kullu Valley, Himalayas. Geosciences 2019, 9, 156. [Google Scholar] [CrossRef] [Green Version]
  39. Milevski, I.; Dragićević, S. Landslides susceptibility zonation of the territory of north macedonia using analytical hierarchy process approach. Math. Biotech. Sci. 2019, 40, 115–126. [Google Scholar] [CrossRef] [Green Version]
  40. Horriche, F.J.; Benabdallah, S. Assessing Aquifer Water Level and Salinity for a Managed Artificial Recharge Site Using Reclaimed Water. Water 2020, 12, 341. [Google Scholar] [CrossRef] [Green Version]
  41. Al-Shabeeb, A.A.-R.; Al-Adamat, R.; Al-Fugara, A.; Al Amoush, H.; Alayyash, S. Delineating groundwater potential zones within the Azraq Basin of Central Jordan using multi-criteria GIS analysis. Groundw. Sustain. Dev. 2018, 7, 82–90. [Google Scholar] [CrossRef]
  42. Hamdani, N.; Baali, A. Height Above Nearest Drainage (HAND) model coupled with lineament mapping for delineating groundwater potential areas (GPA). Groundw. Sustain. Dev. 2019, 9, 100256. [Google Scholar] [CrossRef]
  43. Sophiya, M.S.; Syed, T.H. Assessment of vulnerability to seawater intrusion and potential remediation measures for coastal aquifers: A case study from eastern India. Environ. Earth Sci. 2013, 70, 1197–1209. [Google Scholar] [CrossRef]
  44. Souissi, D.; Msaddek, M.H.; Zouhri, L.; Chenini, I.; May, M.; Dlala, M. Mapping groundwater recharge potential zones in arid region using GIS and Landsat approaches, southeast Tunisia. Hydrol. Sci. J. 2018, 63, 251–268. [Google Scholar] [CrossRef]
  45. Institut National de la Météorologie. 2018 Archive INM Pour la Période de 1985–2015; Institut National de la Météorologie-Tableaux Climatiques Mensuels: Nabeul, Tunisia, 2018. [Google Scholar]
  46. Yeh, H.-F.; Cheng, Y.-S.; Lin, H.-I.; Lee, C.-H. Mapping groundwater recharge potential zone using a GIS approach in Hualian River, Taiwan. Sustain. Environ. Res. 2016, 26, 33–43. [Google Scholar] [CrossRef] [Green Version]
  47. Sar, N.; Khan, A.; Chatterjee, S.; Das, A. Hydrologic delineation of ground water potential zones using geospatial technique for Keleghai river basin, India. Model. Earth Syst. Environ. 2015, 1, 25. [Google Scholar] [CrossRef] [Green Version]
  48. Rahmati, O.; Samani, A.A.N.; Mahdavi, M.; Pourghasemi, H.R.; Zeinivand, H. Groundwater potential mapping at Kurdistan region of Iran using analytic hierarchy process and GIS. Arab. J. Geosci. 2014, 8, 7059–7071. [Google Scholar] [CrossRef]
  49. Adham, M.I.; Jahan, C.S.; Mazumder, Q.H.; Hossain, M.M.A.; Haque, A.-M. Study on groundwater recharge potentiality of Barind Tract, Rajshahi District, Bangladesh using GIS and Remote Sensing technique. J. Geol. Soc. India 2010, 75, 432–438. [Google Scholar] [CrossRef]
  50. Cai, Z.; Ofterdinger, U.S. Analysis of groundwater-level response to rainfall and estimation of annual recharge in fractured hard rock aquifers, NW Ireland. J. Hydrol. 2016, 535, 71–84. [Google Scholar] [CrossRef] [Green Version]
  51. Edet, A.E.; Okereke, C.S.; Teme, S.C.; Esu, E.O. Application of remote sensing data to groundwater exploration: A case study of the Cross River State, southeastern Nigeria. Hydrogeol. J. 1998, 6, 393–404. [Google Scholar] [CrossRef]
  52. Su, Z. Remote sensing of land use and vegetation for mesoscale hydrological studies. Int. J. Remote Sens. 2000, 21, 213–233. [Google Scholar] [CrossRef]
  53. Shaban, A.; Khawlie, M.; Abdallah, C. Use of remote sensing and GIS to determine recharge potential zones: The case of Occidental Lebanon. Hydrogeol. J. 2005, 14, 433–443. [Google Scholar] [CrossRef]
  54. CRDA (Commissariat Régionale au Développement Agricole). Carte des Périmètres Irrigués de la Côte Orientale de Cap-Bon; Rapport D’activités de Ministère de L’agriculture et des Ressources Hydrauliques; CRDA: Nabeul, Tunisia, 2012; p. 64. [Google Scholar]
  55. Ismail, E.; El-Sayed, E.; Sakr, S.; Youssef, E. Characteristic of groundwater potentialities in West Nile Valley South, Minia Governorate, Egypt. Arab. J. Geosci. 2017, 10, 521. [Google Scholar] [CrossRef]
  56. Mageshkumar, P.; Subbaiyan, A.; Lakshmanan, E.; Thirumoorthy, P. Application of geospatial techniques in delineating groundwater potential zones: A case study from South India. Arab. J. Geosci. 2019, 12, 151. [Google Scholar] [CrossRef]
  57. Yeh, H.-F.; Lee, C.-H.; Hsu, K.-C.; Chang, P.-H. GIS for the assessment of the groundwater recharge potential zone. Environ. Earth Sci. 2008, 58, 185–195. [Google Scholar] [CrossRef]
  58. Saaty, T.L. The Analytic Hierarchy Process: Planning, Priority Setting, Resource Allocation; McGraw: New York, NY, USA, 1980; p. 281. [Google Scholar]
  59. Saaty, T.L. Decision Making for Leaders: The Analytic Hierarchy Process for Decisions in a Complex World; RWS Publications: Pittsburgh, PA, USA, 1990. [Google Scholar]
  60. Saaty, T.L. Fundamentals of the Analytic Network Process. In Proceedings of the International Symposium of the Analytic Hierarchy Process (ISAHP), Kobe, Japan, 12–14 August 1999. [Google Scholar]
  61. Saaty, T.L. Fundamentals of the analytic network process—Multiple networks with benefits, costs, opportunities and risks. J. Syst. Sci. Syst. Eng. 2004, 13, 348–379. [Google Scholar] [CrossRef]
  62. Saaty, T.L. Decision making with the analytic hierarchy process. Int. J. Serv. Sci. 2008, 1, 83. [Google Scholar] [CrossRef] [Green Version]
  63. Rezaei-Moghaddam, K.; Karami, E. A multiple criteria evaluation of sustainable agricultural development models using AHP. Environ. Dev. Sustain. 2007, 10, 407–426. [Google Scholar] [CrossRef]
  64. Sahoo, S.; Dhar, A.; Kar, A.; Ram, P. Grey analytic hierarchy process applied to effectiveness evaluation for groundwater potential zone delineation. Geocarto Int. 2016, 32, 1188–1205. [Google Scholar] [CrossRef]
  65. Pang, Z.; Rahaman, F.; Arefin, R.; Ali, S.; Mazumder, Q.H. Delineation of groundwater potential zones of Atrai–Sib river basin in north-west Bangladesh using remote sensing and GIS techniques. Sustain. Water Resour. Manag. 2018, 5, 689–702. [Google Scholar] [CrossRef]
  66. Valverde, J.P.B.; Blank, C.; Roidt, M.; Schneider, L.; Stefan, C. 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]
  67. Brunnelli, M. Introduction to the Analytic Hierarchy Process; Springer: Amsterdam, The Netherlands, 2015. [Google Scholar]
  68. Pilevar, A.R.; Matinfar, H.R.; Sohrabi, A.; Sarmadian, F. Integrated fuzzy, AHP and GIS techniques for land suitability assessment in semi-arid regions for wheat and maize farming. Ecol. Indic. 2020, 110, 105887. [Google Scholar] [CrossRef]
  69. Leake, C.; Malczewski, J. GIS and Multicriteria Decision Analysis. J. Oper. Res. Soc. 2000, 51, 247. [Google Scholar] [CrossRef]
  70. Chatterjee, R.; Pranjal, P.; Jally, S.; Kumar, B.; Dadhwal, V.K.; Srivastav, S.; Kumar, D. Potential groundwater recharge in north-western India vs spaceborne GRACE gravity anomaly based monsoonal groundwater storage change for evaluation of groundwater potential and sustainability. Groundw. Sustain. Dev. 2020, 10, 100307. [Google Scholar] [CrossRef]
  71. Malczewski, J. GIS-based multicriteria decision analysis: A survey of the literature. Int. J. Geogr. Inf. Sci. 2006, 20, 703–726. [Google Scholar] [CrossRef]
  72. Senanayake, I.P.; Dissanayake, D.; Mayadunna, B.; Weerasekera, W. An approach to delineate groundwater recharge potential sites in Ambalantota, Sri Lanka using GIS techniques. Geosci. Front. 2016, 7, 115–124. [Google Scholar] [CrossRef] [Green Version]
  73. Rao, B.V.; Briz-Kishore, B. A methodology for locating potential aquifers in a typical semi-arid region in India using resistivity and hydrogeologic parameters. Geoexploration 1991, 27, 55–64. [Google Scholar] [CrossRef]
  74. Berhanu, K.G.; Hatiye, S.D. Identification of Groundwater Potential Zones Using Proxy Data: Case study of Megech Watershed, Ethiopia. J. Hydrol. Reg. Stud. 2020, 28, 100676. [Google Scholar] [CrossRef]
  75. Patra, S.; Mishra, P.; Mahapatra, S.C. Delineation of groundwater potential zone for sustainable development: A case study from Ganga Alluvial Plain covering Hooghly district of India using remote sensing, geographic information system and analytic hierarchy process. J. Clean. Prod. 2018, 172, 2485–2502. [Google Scholar] [CrossRef]
  76. Lodwick, W.A.; Monson, W.; Svoboda, L. Attribute error and sensitivity analysis of map operations in geographical informations systems: Suitability analysis. Int. J. Geogr. Inf. Syst. 1990, 4, 413–428. [Google Scholar] [CrossRef]
  77. Rahman, A. A GIS based DRASTIC model for assessing groundwater vulnerability in shallow aquifer in Aligarh, India. Appl. Geogr. 2008, 28, 32–53. [Google Scholar] [CrossRef]
  78. Maity, D.K.; Mandal, S. Identification of groundwater potential zones of the Kumari river basin, India: An RS & GIS based semi-quantitative approach. Environ. Dev. Sustain. 2017, 21, 1013–1034. [Google Scholar] [CrossRef]
  79. Janizadeh, S.; Avand, M.; Jaafari, A.; Phong, T.V.; Bayat, M.; Ahmadisharaf, E.; Prakash, I.; Pham, B.T.; Lee, S. Prediction Success of Machine Learning Methods for Flash Flood usceptibility Mapping in the Tafresh Watershed, Iran. Sustainability 2019, 11, 5426. [Google Scholar] [CrossRef] [Green Version]
  80. Pal, S.; Kundu, S.; Mahato, S. Groundwater potential zones for sustainable management plans in a river basin of India and Bangladesh. J. Clean. Prod. 2020, 257, 120311. [Google Scholar] [CrossRef]
  81. Malczewski, J.; Rinner, C. Multicriteria Decision Analysis in Geographic Information Science; Springer Science and Business Media LLC.: Berlin, Germany, 2015. [Google Scholar]
  82. Nguyen, V.-T.; Tran, T.H.; Ha, N.A.; Tsangaratos, P.; Al-Ansari, N.; Van Phong, T.; Nguyen, D.H.; Malek, M.; Amini, A.; Prakash, I.; et al. GIS Based Novel Hybrid Computational Intelligence Models for Mapping Landslide Susceptibility: A Case Study at Da Lat City, Vietnam. Sustainability 2019, 11, 7118. [Google Scholar] [CrossRef] [Green Version]
  83. Anderson, A.; Khorram, M.A. Exposure to D.D.T. BMJ 1948, 1, 1132–1134. [Google Scholar] [CrossRef] [PubMed]
  84. DGRE (Direction Générale des Ressources en Eau). Annuaires Pluviométriques (1998–2017); DGRE: Tunis, Tunisia, 2018. [Google Scholar]
  85. Naghibi, S.A.; Pourghasemi, H.R.; Dixon, B. GIS-based groundwater potential mapping using boosted regression tree, classification and regression tree, and random forest machine learning models in Iran. Environ. Monit. Assess. 2015, 188, 1–27. [Google Scholar] [CrossRef]
  86. Jaafari, A.; Pourghasemi, H.R. Factors Influencing Regional-Scale Wildfire Probability in Iran. In Spatial Modeling in GIS and R for Earth and Environmental Sciences; Elsevier: Amsterdam, The Netherlands, 2019; pp. 607–619. [Google Scholar]
  87. Das, S. Comparison among influencing factor, frequency ratio, and analytical hierarchy process techniques for groundwater potential zonation in Vaitarna basin, Maharashtra, India. Groundw. Sustain. Dev. 2019, 8, 617–629. [Google Scholar] [CrossRef]
  88. Napolitano, P.; Fabbri, A.G. Single-parameter sensitivity analysis for aquifer vulnerability assessment using DRASTIC and SINTACS. In Application of Geographic Information Systems in Hydrology and Water Resources Management; IAHS: Vienne, Austria, 1996; pp. 559–566. [Google Scholar]
  89. United Nations. Hydrogeologic Map of Lebanon; Carte Hydrogelogique du Liban au 1/100000 Me; United Nations: Beyrouth, Lebanon, 1967. [Google Scholar]
Figure 1. Location map of Tunisia and the Korba region located in the Cap-Bon Peninsula (projection system: TUN UTM; datum: WGS_1984; zone: 32).
Figure 1. Location map of Tunisia and the Korba region located in the Cap-Bon Peninsula (projection system: TUN UTM; datum: WGS_1984; zone: 32).
Water 12 02525 g001
Figure 2. (a) Geologic and structures map of the Korba aquifer: ((1) Miocene sands, (2) Pliocene sands, (3) Tyrrhenian sandstones, (4) Quaternary alluvial deposits, (5) Marine Quaternary) and (b) Schematic cross section of the aquifer showing the subsurface hydro-geological and groundwater flow movements [2].
Figure 2. (a) Geologic and structures map of the Korba aquifer: ((1) Miocene sands, (2) Pliocene sands, (3) Tyrrhenian sandstones, (4) Quaternary alluvial deposits, (5) Marine Quaternary) and (b) Schematic cross section of the aquifer showing the subsurface hydro-geological and groundwater flow movements [2].
Water 12 02525 g002
Figure 3. Potentiometric map evolution of the Korba aquifer: (a) 1965 [10]; (b) 1977 [10]; (c) 2014 [8] and (d) 2018 (Chiba, M’Laabi and Lebna dams were built in 1963, 1964 and 1986, respectively).
Figure 3. Potentiometric map evolution of the Korba aquifer: (a) 1965 [10]; (b) 1977 [10]; (c) 2014 [8] and (d) 2018 (Chiba, M’Laabi and Lebna dams were built in 1963, 1964 and 1986, respectively).
Water 12 02525 g003
Figure 4. Framework to delineate groundwater recharge zones.
Figure 4. Framework to delineate groundwater recharge zones.
Water 12 02525 g004
Figure 5. Interrelationships between multiple factors influencing groundwater recharge.
Figure 5. Interrelationships between multiple factors influencing groundwater recharge.
Water 12 02525 g005
Figure 6. A schematic of thematic layer overlay analysis using geographic information system (GIS) to demarcate groundwater recharge potential zones.
Figure 6. A schematic of thematic layer overlay analysis using geographic information system (GIS) to demarcate groundwater recharge potential zones.
Water 12 02525 g006
Figure 7. The eight thematic layers used in the multi-criteria decision making (MCDM) frameworks to map groundwater recharge potential: (a) Lithology, (b) Land use/Land cover, (c) Slope, (d) Geomorphology, (e) Lineament density map, (f) Rainfall, (g) Drainage density and (h) soil types.
Figure 7. The eight thematic layers used in the multi-criteria decision making (MCDM) frameworks to map groundwater recharge potential: (a) Lithology, (b) Land use/Land cover, (c) Slope, (d) Geomorphology, (e) Lineament density map, (f) Rainfall, (g) Drainage density and (h) soil types.
Water 12 02525 g007
Figure 8. Groundwater recharge zones of the Korba area: ((a) using AHP method and (b) MIF technique) and validation of groundwater recharge zones through comparison with water level profile in different parts of the study area: (c) AA’ for AHP and (d) BB’ for MIF, respectively, indicated the W-E and NW-SE directions.
Figure 8. Groundwater recharge zones of the Korba area: ((a) using AHP method and (b) MIF technique) and validation of groundwater recharge zones through comparison with water level profile in different parts of the study area: (c) AA’ for AHP and (d) BB’ for MIF, respectively, indicated the W-E and NW-SE directions.
Water 12 02525 g008
Figure 9. Distribution of groundwater recharge classes for AHT and MIF in the study area.
Figure 9. Distribution of groundwater recharge classes for AHT and MIF in the study area.
Water 12 02525 g009
Figure 10. Average of Receiver Operating Characteristic (ROC) curves for AHP and MIF methods.
Figure 10. Average of Receiver Operating Characteristic (ROC) curves for AHP and MIF methods.
Water 12 02525 g010
Table 1. Summary of data, sources and purpose of use.
Table 1. Summary of data, sources and purpose of use.
Data TypeSourceFormatProduct
Digital elevation model (DEM)Shuttle Radar Topography Mission (SRTM) data (USGS), Resolution: 30 mDigital
Drainage densityShuttle Radar Topography Mission (SRTM) data (USGS), Resolution: 30 mRasterDD
SlopeDeveloped from SRTM data (resolution: 30 m) in Arc GISDigitalSLO
RainfallAnnual rainfall data (1980–2017) from Institut National de la MétéorologieTableRN
Land Use/CoverPrepared from Landsat 8, NIR Band with resolution 30 mRasterLU/LC
SoilSoil texture map La Carte Agricole of Nabeul governorate (CRDA Nabeul)DigitalSL
GeomorphologySoil texture map La Carte Agricole of Nabeul governorate (CRDA Nabeul)DigitalGM
Lineament densitySentinel 2 data (European Space Agency) and SRTM (DEM) data (U.S. Geological Survey)RasterLD
LithologyGeological map La Carte Agricole of Nabeul governorate (CRDA Nabeul)DigitalLI
Table 2. Analytic Hierarchy Process (AHP) relative class rate scale according to Saaty [58].
Table 2. Analytic Hierarchy Process (AHP) relative class rate scale according to Saaty [58].
ImportanceEqualWeak ModerateModerate PlusStrongStrong PlusVery StrongVery, Very StrongExtreme
Scale123456789
1/91/81/71/61/51/41/31/2123456789
Less important
Water 12 02525 i001
More important
Water 12 02525 i002
Table 3. Pair-wise comparison matrix for all parameters.
Table 3. Pair-wise comparison matrix for all parameters.
ParameterLILU/LCSLOGMLDRNDDSL
Lithology (LI)1.002.003.004.005.006.007.008.00
Land use/Land cover (LU/LC)1/21.002.003.004.0056.007.00
Slope (SLO)1/31/21.002.003.0045.006.00
Geomorphology (GM)1/41/31/21.002.0034.005.00
Lineaments (LD)1/51/41/31/21.0023.004.00
Rainfall (RN)1/61/51/41/31/212.003.00
Drainage (DD)1/71/61/51/41/31/21.002.00
Soil (SL)1/81/71/61/51/41/31/21.00
TOTALS2.724.597.4511.2816.0821.8328.5036.00
Table 4. Standardized pairwise comparison matrix and weight factors influencing recharge.
Table 4. Standardized pairwise comparison matrix and weight factors influencing recharge.
ParameterLILU/LCSLOGMLDRNDDSLEigen VectorInfluence (%)AHP Weight
Lithology (LI)0.370.440.40.350.310.270.250.220.440.070.33
Land use/Land cover (LULC)0.180.220.270.270.250.230.210.190.1919.30.23
Slope (SLO)0.120.120.130.180.190.180.160.170.1211.450.16
Geomorphology (GM)0.100.070.070.090.120.140.140.140.19.510.11
Lineaments (LD)0.070.050.040.040.060.090.120.110.066.280.07
Rainfall (RN)0.060.040.030.030.030.050.070.080.055.190.05
Drainage (DD)0.050.040.030.020.020.020.040.060.044.10.03
Soil (SL)0.050.030.020.020.020.020.020.030.044.10.02
Total effect1.001.001.001.001.001.001.001.001.00100.001.00
Table 5. Calculation of the principal Eigenvalue to rank parameter influence.
Table 5. Calculation of the principal Eigenvalue to rank parameter influence.
Thematic Map(1) Total Relative Weight of Each Factor (from Table 3)(2) Eigenvector Value of Each Factor (from Table 4)Eigenvalues (1) × (2)
Lithology2.720.41.09
Land use/Land cover4.590.190.87
Slope7.450.120.89
Geomorphology11.280.11.13
Lineaments16.080.060.96
Rainfall21.830.051.09
Drainage28.500.041.14
Soil36.000.041.44
Principal Eigenvalue (λmax) 8.62
Table 6. Random inconsistency indices [58].
Table 6. Random inconsistency indices [58].
n23 456789
RI00.580.91.121.241.321.411.49
Table 7. The values of major, minor, relative effects and proposed weight of influencing factors.
Table 7. The values of major, minor, relative effects and proposed weight of influencing factors.
FactorMajor EffectMinor EffectRelative EffectProposed Weight
Lithology1 + 1 + 1 +10422
Land use/Land cover1 + 10.5 + 0.5 + 0.53.519
Slope1 + 10.52.513
Geomorphology1 + 10.52.513
Lineaments1 + 10211
Rainfall10.51.58
Drainage10.51.58
Soil1016
Total∑18.5∑100
Table 8. Criteria-assigned ratings and weights for AHP and Multi-Influencing Factors (MIF) methods.
Table 8. Criteria-assigned ratings and weights for AHP and Multi-Influencing Factors (MIF) methods.
FactorDomain of EffectRankMIFAHP
WeightWeighted RatingTotal WeightWeightWeighted RatingTotal Weight
LithologyTyrrhenian sandstone60.221.325.940.331.988.91
Sandy sandstone61.321.98
Silty sand51.11.65
Sand clay–sand40.881.32
Quaternary alluvial deposits40.881.32
Yellow sand with clay20.440.66
Land use/Land coverWadis (ephemeral rivers)60.191.142.850.231.383.45
Grassland + Agricultural land50.951.15
Barren land30.570.69
Buildings10.190.23
Slope0–5%60.130.781.820.160.962.24
5–10%40.520.64
10–15%20.260.32
15–20%10.130.16
>20%10.130.16
GeomorphologyShallow Flood Plain; Beach sand60.130.782.340.110.661.98
Shallow buried pediplain50.650.55
Moderately buried pediplain40.520.44
Deep buried pediplain20.260.22
Sedimentary high ground10.130.11
Lineament density (km/km2)>260.110.661.980.070.421.26
1–250.550.35
1–1.540.440.28
0.5–120.220.14
0–0.510.110.07
Rainfall450–500 mm60.080.480.80.050.30.5
400–450 mm40.320.2
Drainage density (km/km2)>460.080.481.440.030.180.54
3–450.40.15
2–340.320.12
1–220.160.06
0–110.080.03
SoilGravel and sand60.060.360.720.020.120.24
Loam + sandy loam40.240.08
Clay + clay loam20.120.04
Table 9. Statistics of the sensitivity analysis index numbers.
Table 9. Statistics of the sensitivity analysis index numbers.
Layer RemovedSensitivity Analysis Index, S (%)
MinimumMaximumMean Variation Index (SMVI))Standard Deviation (Sd)
AHPMIFAHPMIFAHPMIFAHPMIF
Slope0.210.161.651.561.041.010.280.20
Rainfall0.380.201.711.070.480.290.210.30
Soil0.080.060.390.320.270.210.130.20
Lineament density0.060.051.361.380.640.690.090.08
Land use/Land cover0.720.651.992.021.241.260.200.24
Lithology0.210.192.352.121.621.540.160.14
Geomorphology0.710.661.571.620.740.710.070.09
Drainage density0.130.090.760.800.310.290.140.17
Table 10. Statistics of the single parameter sensitivity analysis.
Table 10. Statistics of the single parameter sensitivity analysis.
FactorEffective Weighting, %Assigned Weights
MinimumMaximumMeanStandard Deviation (Sd)
AHPMIFAHPMIFAHPMIFAHPMIFAHPMIF
Slope161396784436282216%13%
Rainfall2032304825407115%8%
Soil41212368244122%6%
Lineament density7114266253914227%11%
Land use/Land cover23191381148671504223%19%
Lithology6644198132664460.252.433%22%
Geomorphology111366783946222611%13%
Drainage density38184810286163%8%
Table 11. Groundwater recharge potential zones and categorization.
Table 11. Groundwater recharge potential zones and categorization.
Delineated SurfaceGroundwater Recharge Zones
Very GoodGoodModeratePoorVery Poor
AHPMIFAHPMIFAHPMIFAHPMIFAHPMIF
Area (%)16.0030.7053.0050.0011.509.6015.808.072.901.58
Area (km2)66.88128.32221.54209.0048.0740.1366.0433.7312.126.60
Average (%) [89]47.5032.5015.007.502.50

Share and Cite

MDPI and ACS Style

Zghibi, A.; Mirchi, A.; Msaddek, M.H.; Merzougui, A.; Zouhri, L.; Taupin, J.-D.; Chekirbane, A.; Chenini, I.; Tarhouni, J. Using Analytical Hierarchy Process and Multi-Influencing Factors to Map Groundwater Recharge Zones in a Semi-Arid Mediterranean Coastal Aquifer. Water 2020, 12, 2525. https://doi.org/10.3390/w12092525

AMA Style

Zghibi A, Mirchi A, Msaddek MH, Merzougui A, Zouhri L, Taupin J-D, Chekirbane A, Chenini I, Tarhouni J. Using Analytical Hierarchy Process and Multi-Influencing Factors to Map Groundwater Recharge Zones in a Semi-Arid Mediterranean Coastal Aquifer. Water. 2020; 12(9):2525. https://doi.org/10.3390/w12092525

Chicago/Turabian Style

Zghibi, Adel, Ali Mirchi, Mohamed Haythem Msaddek, Amira Merzougui, Lahcen Zouhri, Jean-Denis Taupin, Anis Chekirbane, Ismail Chenini, and Jamila Tarhouni. 2020. "Using Analytical Hierarchy Process and Multi-Influencing Factors to Map Groundwater Recharge Zones in a Semi-Arid Mediterranean Coastal Aquifer" Water 12, no. 9: 2525. https://doi.org/10.3390/w12092525

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