A Multidisciplinary Approach for Groundwater Potential Mapping in a Fractured Semi-Arid Terrain (Kerdous Inlier, Western Anti-Atlas, Morocco)

<jats:p>This study is focused on developing an approach for spatial mapping of groundwater by considering four types of factors (geological, topographical, hydrological, and climatic factors), and by using different bivariate statistical models, such as frequency ratio (FR) and Shannon’s entropy (SE). The developed approach was applied in a fractured aquifer basin (Ameln Basin, Western Anti-Atlas, Morocco), to map the spatial variation of groundwater potential. Fifteen factors (15) influencing groundwater were considered in this study, including slope degree, slope aspect, elevation, topographic wetness index (TWI), slope length (LS), topographic position index (TPI), plane curvature, profile curvature, drainage density, lineament density, distance to rivers and fault network, normalized difference vegetation index (NDVI), lithology, and land surface temperature (LST). The potential maps produced were then classified into five classes to illustrate the spatial view of each potential class obtained. The predictive capacity of the frequency ratio and Shannon’s entropy models was determined using two different methods, the first one based on the use of flow data from 49 boreholes drilled in the study area, to test and statistically calibrate the predictive capacity of each model. The results show that the percentage of positive water points corresponds to the most productive areas (high water flow) (42.86% and 30.61% for the FR and SE models, respectively). On the other hand, the low water flows are consistent with the predicted unfavorable areas for hydrogeological prospecting (4.08% for the FR model and 6.12% for the SE model). Additionally, the second validation method involves the integration of 7200 Hz apparent resistivity data to identify conductive zones that are groundwater circulation zones. The interpretation of the geophysical results shows that the high-potential zones match with low apparent resistivity zones, and therefore promising targets for hydrogeological investigation. The FR and SE models have proved very efficient for hydrogeological mapping at a fractured basement area and suggest that the northern and southern part of the study area, specifically the two major fault zones (Ameln Valley in the north, and the Tighmi-Tifermit Valley in the south) has an adequate availability of groundwater, whereas the central part, covering the localities of Tarçouat, Boutabi, Tililan, and Ighalen, presents a scarcity of groundwater. The trend histogram of the evolution of positive water points according to each potentiality class obtained suggests that the FR model was more accurate than the SE model in predicting the potential groundwater areas. The results suggest that the proposed approach is very important for hydrogeological mapping of fractured aquifers, and the resulting maps can be helpful to managers and planners to generate groundwater development plans and attenuate the consequences of future drought.</jats:p>

section_content: Water is a precious and vital resource around the world. It is a resource that is becoming increasingly scarce in many parts of the world. Moreover, the request for drinking water to satisfy the water needs of the population is increasingly strong. Since Morocco is one of the North African countries that is mainly characterized by its arid to semi-arid bioclimatic stage, the country is heading for a water crisis mainly caused by its inefficient water resources management policy, spatiotemporal water irregularity, and drought caused by climate change [1] [2] [3] [4] [5] [6] [7] [8] [9]. Groundwater resources in the basement zone represent a major resource for rural populations in the mountainous areas of Morocco. Since surface water is not perennial, groundwater conditions the food security of these populations [10] [11] [12] [13] [14]. Basement aquifers are generally characterized by a geographically irregularly distributed water resource. Access to this groundwater for domestic use is achieved by drilling wells or boreholes capable of providing quality water and sustainability over time. These wells should be located in aquifers with significant storage and/or recharge. The scarce groundwater in the Kerdous crystalline basement (Ameln Basin, Western Anti-Atlas, Morocco) developed in a fractured environment and is distributed illegally throughout the basin. Their sustainable management represents a major socio-economic development challenge for the region [6, [15] [16] [17]. This basin is located in a complex geological setting where the Kerdous crystalline basement formations are dominant. In such environments, groundwater resources can hardly supply the increasing demand induced by the demographic evolution of rural populations [18] [19] [20]. Therefore, the identification of an appropriate method for groundwater management and prediction at regional and local levels is necessary, using efficient, accurate, advanced, and cost-effective modeling methods such as remote sensing and geographic information systems [21] [22] [23] because, generally, conventional approaches (geophysical, hydrochemical, etc. ) are expensive, less accurate, and take too much time [24, 25]. Remote sensing (RS) and geographic information systems (GIS) are widely applied in water resources research. GIS is a useful tool that has some advantages for acquiring, storing, analyzing, modeling, archiving, and sharing geographic information to manage a massive quantity of spatial data used in the decision-making process [26] ; on the other hand, remote sensing facilitates a wide range of studies and a synoptic view of spatio-temporal distribution observations in a better and faster assessment [27]. Exploration using RS and GIS integration has recently attracted attention as a cost-effective and efficient method [28, 29]. In this work, groundwater potential zones were delineated by integrating various thematic maps (NDVI, LST, and lineament density) generated from remote-sensing data (Landsat 8 Oli); on the other hand, other criteria such as slope, slope aspect, elevation, topographic wetness index (TWI), slope length (LS), topographic position index (TPI), plane curvature, profile curvature, drainage density, distance to rivers, and fault network were plotted in a GIS environment. There are several hydrogeological mapping methods, such as the analytical hierarchy process (AHP) [6, [30] [31] [32], analytical network process (ANP) [33], multi-influence factor analysis (MIF) [34, 35], boosted regression tree [36], advanced machine learning algorithms [37], weighted overlay analysis [38], frequency ratio model [39] [40] [41] [42], and Shannon's entropy models [42] [43] [44] were used to identify groundwater exploration areas and artificial recharge sites. The application of the frequency ratio and Shannon's entropy model has not yet been used in Morocco, especially for hydrogeological mapping in an environment that is characterized by fracture porosity [45] [46] [47]. FR is a simple bivariate statistical tool to calculate the probabilistic relationship between independent and dependent variables; it is a very simple and flexible technique to delineate potential areas of discontinuous aquifers [48, 49]. In contrast, Shannon's entropy model represents the average unpredictability of a random variable. Groundwater potential entropy refers to the extent to which various controlling groundwater occurrences influence groundwater productivity. References [41, 44] conducted a comparative study between FR and entropy index, for mapping the potential of qanats in the Moghan watershed, Iran. They concluded that both techniques are effective for delineating areas with similar potential and climate as compared to our case study. The main objective of this study is to delineate the groundwater availability of the aquifer system in the study area using the FR and SE models. The resulting groundwater potential maps will be useful for water resource managers to efficiently guide future drilling and, consequently, sustainably manage groundwater resources.
section_name: Introduction

section_content: The Ameln Basin is located in the northern part of the Tiznit district and covers an area of 712 km 2 with an average altitude of 1324 m (Figure 1 ). From a climatic point of view, the study area is characterized by a semi-arid climate, where the minimum and the maximum temperatures are 2 and 41 • C, respectively, and the average annual precipitation is about 287 mm (Figure 2 ). The most current land uses in this region are the building areas, almond trees, and argan trees, with an abundance of skeletal soils and superficial rock. They are classified as soils of lesser quality; consequently, the agricultural production potential of these soils is low and the forest cover remains low throughout the territory. The hydrographic network in the study area is characterized by a strong seasonal and interannual irregularity of the flows. Floods are often violent and rapid. The rivers are not perennial and dry most of the year because of the aridity of the region's climate, which makes the exploitation of these water resources difficult. From a geologic point of view, the study area is part of the Anti-Atlas (AA) belt of southern Morocco that is located at the northern boundary of the West African Craton (WAC) (Figure 3a ). The AA belt corresponds to a broad mountain chain following the direction ENE-WSW and it is parallel to the High-Atlas belt. The contact between the aforementioned belts is marked by the southern Atlasic Fault. The Precambrian terrain (old units) outcrops as "inliers" that are distributed along two major faults zones (the southern Atlasic Fault and the Anti-Atlas major fault) (Figure 3b ) [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61]. The Kerdous inlier, where the study area is located, is the largest one in the western Anti-Atlas. The age of the units forming this terrain range from Paleoproterozoic to Cambrian, where they are distinguished from each other by tectonics unconformities [53, [62] [63] [64] [65]. The Paleoproterozoic of the Ameln Basin includes a weakly metamorphic schist series: green-schists, sericite-schists, mica-schists, basic volcanic rocks (gabbro and dolerites), and granites that crosscut them all, forming a metamorphism contact. The tectonic directions recorded in these old formations are E-W and ESE-WNW [62, [64] [65] [66]. The Upper Neoproterozoic formations include thick quartzite mostly corresponding to the major ridges (Jbel Lkst of Kerdous), basic volcanic rocks (gabbros and dolerites), and schists with conglomerates and basic vulcanites. After successive volcanic events, the Anti-Atlasic domain recognized a great transgression from the west forming the Lower Paleozoic cover that consists of a transgressive limestone or dolomitic level [67]. The Quaternary deposits consist of alluvial deposits that are forming the small generalized aquifers along the major valleys in the study area (Figure 4 ).
section_name: Study Area

section_name: Materials and Methods

section_content: In this study, to delineate the groundwater potential areas in the Ameln Basin, several GIS-based techniques were used, including numerical, cartographic, and geostatistical techniques. The methodology proposed for modeling aquifer productivity involves four steps: (1) preparation of data recorded from boreholes in the study area, (2) influencing factors data collection and spatial database preparation, (3) productive zones evaluation via the relationship between borehole data and influencing factors, and (4) results from validation and interpretation. The geospatial data were achieved using ArcGIS 10. 4. 1 software and the weights generation for the groundwater exploration factors and the associated functions were computed using Microsoft Excel. The validation of results was prepared as a histogram showing the productivity of each obtained groundwater potential class. The final map of the groundwater potentiality was achieved by multiplying each factor with its weight using a raster calculation tool in ArcGIS 10. 4. 1 software. Figure 5 shows the methodology applied in this study. 

Water 2022, 14, x FOR PEER REVIEW 6 of 26
section_name: GIS Applications

section_content: In this study, to delineate the groundwater potential areas in the Ameln Basin, several GIS-based techniques were used, including numerical, cartographic, and geostatistical techniques. The methodology proposed for modeling aquifer productivity involves four steps: (1) preparation of data recorded from boreholes in the study area, (2) influencing factors data collection and spatial database preparation, (3) productive zones evaluation via the relationship between borehole data and influencing factors, and (4) results from validation and interpretation. The geospatial data were achieved using ArcGIS 10. 4. 1 software and the weights generation for the groundwater exploration factors and the associated functions were computed using Microsoft Excel. The validation of results was prepared as a histogram showing the productivity of each obtained groundwater potential class. The final map of the groundwater potentiality was achieved by multiplying each factor with its weight using a raster calculation tool in ArcGIS 10. 4. 1 software. Figure 5 shows the methodology applied in this study.
section_name: GIS Applications

section_content: In general, mapping groundwater potentiality depends on a variety of conditioning factors, which depend on data available in the study area [68]. In the present study, 15 groundwater influencing factors were used, which are grouped into four categories: topographic, hydrological, geological, and climatologic factors. The integration of these factors in a GIS-based model is anticipated to be an effective tool for predicting groundwater resource potentiality in the whole study area.
section_name: Groundwater Occurrence Factors

section_content: Topographic factors, including slope degree, slope aspect, elevation, topographic wetness index (TWI), slope length (LS), topographic position index (TPI), plane curvature, and profile curvature, are considered essential to the spatial variation of surface water infiltration conditions. In this study, the topographic factors were prepared, in the ArcGIS
section_name: Topographic Factors

section_content: Topographic factors, including slope degree, slope aspect, elevation, topographic wetness index (TWI), slope length (LS), topographic position index (TPI), plane curvature, and profile curvature, are considered essential to the spatial variation of surface water infiltration conditions. In this study, the topographic factors were prepared, in the ArcGIS platform, as individual maps from a digital elevation model (DEM) of 12. 5 m of spatial resolution that was collected free of charge from the Alaska satellite facility site. Among the topographic factors, the slope is considered an important one for identifying suitable locations for surface water infiltration, where low slope areas correspond to an area with low surface runoff and high percolation rates, while high slope areas indicate areas with high and rapid runoff [69]. The slope aspect identifies the downslope direction of the maximum rate of change in value from each pixel to its neighbors [70]. The slope aspect strongly controls the hydrological parameters, like evapotranspiration, precipitation direction, and weather process, especially in dry and arid areas without vegetation cover [71]. The elevation (altitude) is the height of an object above sea level. The elevation impacts the climate conditions that affect soil infiltration degree, vegetation type, and groundwater occurrence [49, 72]. The high altitudes indicate areas with rapid and high runoff, while the lower altitudes indicate high surface water recharge and then high water infiltration [73]. The topographic wetness index (TWI) is an index used to describe spatial humidity models and explain the effects of the topographic conditions on them [74]. TWI affects runoff movement and accumulation [75]. TWI can be calculated using the following equation:

where As is the cumulative surface of the ascending slope and β is the slope gradient. Slope length (LS) is the combination of slope (S) and slope length (L). LS can be derived from the universal soil loss equation (USLE), which is mostly used to estimate the water erosion and transported sediments rates [75]. Moore and Burch [76] (2) where Bs is the specific catchment area (m 2 ) and β is the local slope gradient (degree). 

The topographic position index (TPI) is an index used to analyze the topographic landscape units in the upper, middle, and lower portions of the watershed [77]. TPI evaluates the elevation of a pixel and the average elevation of pixels around it using the following Equation (3):

where M0 is the elevation of the point (pixel) under evaluation, Mn is the elevation of grid, and n is the total number of surrounding pixels used for evaluation. Curvature refers to the morphology of the topography and consists of three parameters: profile, plane, and a combination of both. Profile curvature and plane curvature mainly affect the acceleration and deceleration of flows, as well as the convergence and divergence of flows [44].
section_name: Topographic Factors

section_content: In general, to estimate surface water flow and to determine groundwater infiltration, hydrological factors like drainage density and distance from the river are to be taken into account. The drainage is controlled by the nature, structure, and lithology of the bedrock, the soil's capacity to infiltrate, type of vegetation, and slope angle [68, 73]. The drainage density is the ratio of the total length of streams per unit of watershed area [78], as expressed in the following equation:

where Dd is drainage density (km -1 ), ∑Li is the total length of streams per unit (km), and A is an area (km 2 ). 

Water 2022, 14, 1553 8 of 26
section_name: Hydrological Factors

section_content: Lithology is considered a crucial factor controlling aquifer porosity and permeability [68, 79]. In addition, lineaments are considered an important geological factor and are linear or curvilinear structures often associated with various geomorphic features or tectonic structures. The orientation is not the only parameter characterizing a lineament network, but its density and connectivity are essential [80]. Lineament density is often used as a hydrodynamic indicator for groundwater resources functionality and determination of productive areas [6, 68]. Furthermore, the distance from lineaments is considered an important factor to determine the hydrogeological interest areas, which are generally very close to these structures [6]. In this study, the lithology map of the study area was prepared from several geological maps (scale 1:50,000) through digitalization in the ArcGIS platform. On the other hand, the lineaments map was prepared by applying Sobel directional filtering on satellite images acquired by Landsat 8 OLI on 29 October 2020. They used satellite data collected free of charge from the US Geological Survey (USGS). First, several processing steps were applied to the collected data in image-processing software (ENVI 5) to extract the lineaments, including radiometric calibration and atmospheric correction using fast line-of-sight atmospheric analysis of spectral hypercubes (FLAASH) model, and transformation by applying principal component analysis (PCA). Afterward, the automatic lineaments extraction was achieved by applying different combinations of values for each parameter of the LINE module in Geomatica software. Then, the obtained map was verified using a faults map that was extracted from geological maps of the study area.
section_name: Geological Factors

section_content: Climatological factors are most important to identify the areas of hydrogeological importance [81, 82]. A normalized difference vegetation index (NDVI) is usually used to estimate the spatiotemporal changes in vegetation distribution [81, 82]. In general, the presence of vegetation cover reflects the presence of near-surface aquifers [41, 83]. In this study, an NDVI map was prepared, from the same Landsat 8 OLI scene used for lineaments extraction, using the following equation:

where ρNIR and ρR are spectral reflectance in the red and near-infrared, respectively. Furthermore, the land surface temperature (LST) was mapped using the following equation:

where ε and BT are emissivity and brightness temperature conversion, respectively, extracted from the thermal infrared bands (10 and 11) of the same Landsat scene imagery.
section_name: Climatological Factors

section_content: This study is based on a combination of remote sensing (RS) and geographic information systems (GIS) to prepare the parameters controlling the availability of groundwater resources. Finally, weights were assigned for different thematic layers and their classes using statistical analysis and to calculate the groundwater potential index using the FR and SE models; the results were interpreted and validated by exploiting the hydrogeological and geophysical data to choose the most adequate model for hydrogeological exploration. Figure 6 shows the role of all the techniques incorporated in this study.
section_name: Modeling Techniques

section_content: The frequency ratio is a quantitative method popularly used to estimate the probability of occurrence and non-occurrence of a given attribute [84]. In this study, the frequency ratio model was used to determine the relationship between each borehole location and a class of the groundwater occurrence factor using Equation ( 7 ):

where A is an area of a class of the factor influencing groundwater, B is a total area of the factor, C is a number of pixels in the class area of the factor, and D is a total number of pixels of the total study area. Moreover, the groundwater productivity based on the FR model was calculated using the equation in following [39, 41, 85] :

where FRi is the frequency ratio of the factor and n is a total number of the used factors.
section_name: Frequency Ratio Model

section_name: Shannon's Entropy Model

section_content: The frequency ratio is a quantitative method popularly used to estimate the probability of occurrence and non-occurrence of a given attribute [84]. In this study, the frequency ratio model was used to determine the relationship between each borehole location and a class of the groundwater occurrence factor using Equation ( 7 ):

where A is an area of a class of the factor influencing groundwater, B is a total area of the factor, C is a number of pixels in the class area of the factor, and D is a total number of pixels of the total study area. Moreover, the groundwater productivity based on the FR model was calculated using the equation in following [39, 41, 85] :

where FRi is the frequency ratio of the factor and n is a total number of the used factors. 

Water 2022, 14, 1553 10 of 26
section_name: Frequency Ratio Model

section_content: Entropy provides a measure of the uncertainty of a random variable system [86], where the groups are equally distributed within the system [87]. Precisely, the entropy model provides the extent of instability, disorder, disequilibrium, and uncertainty in a system [88]. Herein, the application of the entropy model consists of determining the occurrence factors impacting groundwater productivity. Therefore, the entropy value can be used to calculate the system weights index [85]. The information coefficient (Vj) is the weight value of each occurrence-controlling factor [43, 89]. The information coefficient (Vj) can be calculated using the following Equations ( 9 )-( 14 ):

H jmax = log 2 Sj ( 12 )

V j = I j P ij (14) The groundwater potential index based on the SE model (GWPIvj) was calculated using the equation below [43] :

where FR is frequency ratio, Pij is a probability density, Hj and Hjmax are entropy values, Sj is a number of classes, Ij is information coefficient, Vj is the weight value of the factor, i is a number of particular factors maps, Z is a number of classes in the factor map of the whole classes, mi is a number of classes in the particular factor map, C is the class value after secondary classification, and Vj is the weight value of the factor.
section_name: Shannon's Entropy Model

section_name: Results and Discussion

section_content: The slope and altitude maps show that the low slope areas correspond to the central and southern parts of the study area, and to the valley draining the Ameln Basin, while the high slope areas correspond to the quartzites of Jbel Lkest, sandstones of Tizoughrane, and granitic arenas of Tafraout (Figure 7a, b ). In general, due to slow runoff flows in the Ameln Valley, south of Ait Ouafqa and north of the Tarçouat areas located in the north of the study area are more suitable for groundwater recharge. The slope aspect map shows that the valleys (like Ameln and Oumaghouz) characterized by downward-direction slopes oriented to the western side are considered areas of hydrogeological interest (Figure 8a ). In general, the boreholes were installed along these valleys because of the descending underflows from the Precambrian massifs. The topographic wetness index (TWI) map shows that the TWI high values correspond to the valley areas (Ameln Valley and Assif Oumaghouz) and to the hydrographic network areas that have intense ramifications (Figure 8b ). These results prove that the TWI factor is an indicator of possible water drainage. The analysis of the obtained slope length (LS) map reveals that the LS high values correspond to the areas of high altitude (reliefs, like Jbel Lkest, Jbel Azourozour, and Jbel Ayehgd, while the low values of LS correspond to the low altitudes (Figure 9a ). In general, the LS low values are more exposed to water recharge. The high TPI values correspond to the northern part of the Ameln Basin (Jbel Lkest quartzites), while the low TPI values correspond to the low altitude areas (Ameln and Oumaghouz valleys) and TPI values close to zero correspond to the flatlands (Figure 9b ). The topographic wetness index (TWI) map shows that the TWI high values correspond to the valley areas (Ameln Valley and Assif Oumaghouz) and to the hydrographic network areas that have intense ramifications (Figure 8b ). These results prove that the TWI factor is an indicator of possible water drainage. The analysis of the obtained slope length (LS) map reveals that the LS high values correspond to the areas of high altitude (reliefs, like Jbel Lkest, Jbel Azourozour, and Jbel Ayehgd, while the low values of LS correspond to the low altitudes (Figure 9a ). In general, the LS low values are more exposed to water recharge. The high TPI values correspond to the northern part of the Ameln Basin (Jbel Lkest quartzites), while the low TPI values correspond to the low altitude areas (Ameln and Oumaghouz valleys) and TPI values close to zero correspond to the flatlands (Figure 9b ). The topographic wetness index (TWI) map shows that the TWI high values correspond to the valley areas (Ameln Valley and Assif Oumaghouz) and to the hydrographic network areas that have intense ramifications (Figure 8b ). These results prove that the TWI factor is an indicator of possible water drainage. The analysis of the obtained slope length (LS) map reveals that the LS high values correspond to the areas of high altitude (reliefs, like Jbel Lkest, Jbel Azourozour, and Jbel Ayehgd, while the low values of LS correspond to the low altitudes (Figure 9a ). In general, the LS low values are more exposed to water recharge. The high TPI values correspond to the northern part of the Ameln Basin (Jbel Lkest quartzites), while the low TPI values correspond to the low altitude areas (Ameln and Oumaghouz valleys) and TPI values close to zero correspond to the flatlands (Figure 9b ). The idea behind using the profile curvature and plane curvature factors is to high light different aspects of the ground surface shape. In general, profile curvature has par alleled the slope and indicates the direction of the maximum slope. The profile curvature affects the acceleration and deceleration of the water flow at the surface (Figure 9c ). In deed, the profile curvature negative values indicate that the surface is convex, causing dispersion of water. In contrast, the positive profile curvature values indicate that the sur face is concave upward, favoring water stagnation and, then, increasing the water infil tration rate. On the other hand, the plane curvature refers to the convergence and diver gence of water flow at the surface (Figure 9d ). The plane curvature positive values indicate that the surface is a convex causing divergence of surface flows and a consequent loss o the surface water. On the contrary, the plane curvature negative values indicate that the surface is laterally concave causing convergence of the surface flows, which favors water reassembly toward zones favorable to its stagnation. 

In our study area, the suitable areas for water stagnation correspond to the contac area between the great reliefs (quartzite formations of Jbel Lkest and granitic outcrops o Precambrian age) and great valleys of Quaternary fillings. In general, the high drainage The idea behind using the profile curvature and plane curvature factors is to highlight different aspects of the ground surface shape. In general, profile curvature has paralleled the slope and indicates the direction of the maximum slope. The profile curvature affects the acceleration and deceleration of the water flow at the surface (Figure 9c ). Indeed, the profile curvature negative values indicate that the surface is convex, causing dispersion of water. In contrast, the positive profile curvature values indicate that the surface is concave upward, favoring water stagnation and, then, increasing the water infiltration rate. On the other hand, the plane curvature refers to the convergence and divergence of water flow at the surface (Figure 9d ). The plane curvature positive values indicate that the surface is a convex causing divergence of surface flows and a consequent loss of the surface water. On the contrary, the plane curvature negative values indicate that the surface is laterally concave causing convergence of the surface flows, which favors water reassembly toward zones favorable to its stagnation. 

In our study area, the suitable areas for water stagnation correspond to the contact area between the great reliefs (quartzite formations of Jbel Lkest and granitic outcrops of Precambrian age) and great valleys of Quaternary fillings. In general, the high drainage density decreases infiltration and increases surface runoff and is, therefore, not suitable for groundwater development [90], which means that drainage density is an inverse function of permeability. In our study area, the high values of runoff correspond to the Ameln Valley, Assif Idawsmlal, and Assif Oumaghouz (Figure 10a ). Moreover, the criterion distance from river networks had been considered an important factor in hydrogeological research, especially for arid areas, because the local alluvial aquifers usually prospect in areas near the rivers [6]. The obtained results reveal that the areas close to the rivers (200 m) characterized by a thick alluvial level (very permeable formations) are the most promising for water infiltration, while it seems difficult to find water resources in the areas far away from the rivers (distance > 500 m) (Figure 10b ).
section_name: Conditioning Factors

section_content: Valley, Assif Idawsmlal, and Assif Oumaghouz (Figure 10a ). Moreover, the criterion distance from river networks had been considered an important factor in hydrogeological research, especially for arid areas, because the local alluvial aquifers usually prospect in areas near the rivers [6]. The obtained results reveal that the areas close to the rivers (200 m) characterized by a thick alluvial level (very permeable formations) are the most promising for water infiltration, while it seems difficult to find water resources in the areas far away from the rivers (distance > 500 m) (Figure 10b ). Based on the lithology type and geologic age of the formations, the study area was divided into 48 units. These units were classified according to their characteristics and permeabilities: sedimentary rock type (A), sedimentary rock type (B), and igneous rocks (Figure 11a ). The sedimentary type (A) rocks are very permeable of the Quaternary formations, and consist of silts, screes, and pebbles, while the sedimentary type (B) rocks are moderately permeable formations of Adoudounian and Georgian ages, and consist of fractured limestones, dolomites, conglomerates, and shales. On the contrary, the sedimentary type (C) rock consists of igneous rocks that are impermeable to water filtration. The lineaments extraction shows that the NE-SW direction is the major direction detected by the 45° filter. In addition, the lineaments are dense in the north of Tafraout and in Tarçouat (Afra granites) where the geological formations are strongly tectonized and correspond to granitic arenas (Figure 11b ). Furthermore, the areas close to the lineaments (300 m) are the most promising for effective infiltration; however, the effect of this parameter decreases at distances greater than 700 m (Figure 11c ). Based on the lithology type and geologic age of the formations, the study area was divided into 48 units. These units were classified according to their characteristics and permeabilities: sedimentary rock type (A), sedimentary rock type (B), and igneous rocks (Figure 11a ). The sedimentary type (A) rocks are very permeable of the Quaternary formations, and consist of silts, screes, and pebbles, while the sedimentary type (B) rocks are moderately permeable formations of Adoudounian and Georgian ages, and consist of fractured limestones, dolomites, conglomerates, and shales. On the contrary, the sedimentary type (C) rock consists of igneous rocks that are impermeable to water filtration. The lineaments extraction shows that the NE-SW direction is the major direction detected by the 45 • filter. In addition, the lineaments are dense in the north of Tafraout and in Tarçouat (Afra granites) where the geological formations are strongly tectonized and correspond to granitic arenas (Figure 11b ). Furthermore, the areas close to the lineaments (300 m) are the most promising for effective infiltration; however, the effect of this parameter decreases at distances greater than 700 m (Figure 11c ). The NDVI map of the Ameln Basin reveals five interval classes varying between -0. 12 and 0. 38 (Figure 12a ). The high NDVI values correspond to the various valleys and topographic depressions, indicating temperate and humid areas and the increasing possibility of positive groundwater prospection. On the contrary, the NDVI negative values correspond to arid areas of low groundwater potential. On the other hand, the land surface temperature (LST) map shows values ranging from 16 to 39 °C with a mean of 29 °C and a standard deviation of 3. 5 (Figure 12b ). The NDVI map of the Ameln Basin reveals five interval classes varying between -0. 12 and 0. 38 (Figure 12a ). The high NDVI values correspond to the various valleys and topographic depressions, indicating temperate and humid areas and the increasing possibility of positive groundwater prospection. On the contrary, the NDVI negative values correspond to arid areas of low groundwater potential. On the other hand, the land surface temperature (LST) map shows values ranging from 16 to 39 • C with a mean of 29 • C and a standard deviation of 3. 5 (Figure 12b ).
section_name: R REVIEW 13 of 26

section_content: The results of the spatial relationship between the water boreholes and conditioning factors are tabulated in Table 1. Commonly, the ratio value of 1 indicates a medium correlation; on the other hand, the ratio values greater than 1 indicate a high correlation while the values less than 1 indicate a lower correlation e. g., [23, 24, 30, 33]. For elevations between 800 and 1000 m, the FR had a value of 2. 59, which indicates a very high probability of groundwater availability in these altitudes. Meanwhile, the elevation classes more than 1000 m had an FR value of less than 1, which indicates a low probability of groundwater availability, promoting a significant water runoff surface. In contrast, water boreholes are concentrated in areas that have a topographic wetness index (TWI) of more than 8 (FR value 2. 48). In addition, the slope classes 0-5 and 5-15 showed high FR values, 2. 82 and 1. 13, respectively, indicating high groundwater availability. The low FR value (0. 25) corresponds to the slope degree class >30. Of note is that, in this study, the FR decreases where the gradient of the slope increases. The high FR value of 7. 81 corresponds to the slope length ranges of between 10-15. In terms of slope aspect, the water boreholes are dominated by slopes that have a planar aspect and orient northeast and southeast, while being less dominated by west-orienting slopes.
section_name: Frequency Ratio and Shannon Entropy Result

section_content: The results of the spatial relationship between the water boreholes and conditioning factors are tabulated in Table 1. Commonly, the ratio value of 1 indicates a medium correlation; on the other hand, the ratio values greater than 1 indicate a high correlation, while the values less than 1 indicate a lower correlation e. g., [23, 24, 30, 33]. For elevations between 800 and 1000 m, the FR had a value of 2. 59, which indicates a very high probability of groundwater availability in these altitudes. Meanwhile, the elevation classes more than 1000 m had an FR value of less than 1, which indicates a low probability of groundwater availability, promoting a significant water runoff surface. In contrast, water boreholes are concentrated in areas that have a topographic wetness index (TWI) of more than 8 (FR value 2. 48). In addition, the slope classes 0-5 and 5-15 showed high FR values, 2. 82 and 1. 13, respectively, indicating high groundwater availability. The low FR value (0. 25) corresponds to the slope degree class >30. Of note is that, in this study, the FR decreases where the gradient of the slope increases. The high FR value of 7. 81 corresponds to the slope length ranges of between 10-15. In terms of slope aspect, the water boreholes are dominated by slopes that have a planar aspect and orient northeast and southeast, while being less dominated by west-orienting slopes. For plan and profile curvature, the results show that the flat form had the highest FR value (1. 62), while the convex form had the lowest FR value (0. 83), which means that the flat and concave curvature control the groundwater availability in the study area. Furthermore, geologically, the highest values of FR (3. 38) correspond to the highly permeable groups of Quaternary formations (alluvial terrace, ancient Quaternary, and scree), indicating suitable areas of water accumulation and therefore enhancing groundwater recharge. Moreover, the classes of distance from hydrographic networks 0-100 and 100-200 showed the highest FR values of 1. 27 and 0. 61, respectively, which means the areas corresponding to those classes and are effectively zones of the occurrence of groundwater resources. The areas of drainage density of more than 4 km/km 2 have higher values of FR (1. 35 and 2. 10), indicating soil characterized by high infiltration, which means percolation of water to the surface and then intense recharge of the aquifer, which correspond to the alluvium sediments of the Ameln Valley in the study area. In addition, the classes 250-500 m, 500-750 m, and 1000-1250 m of the distance to the lineaments have higher values of FR, 1. 29, 1. 25, and 1. 95, respectively, which mean the areas closer to the lineaments are of importance for groundwater exploration. In our study area, the lineaments density class (0-1, 5) had an FR value of 1, 62 corresponding to the highly tectonized geological formations of Tafraout granite, Tarçouat, and the quartzite of Jbel Lkest. The topographic position index shows that flat ground surfaces (medium slope zones) have higher FR values (0. 95 and 1. 72), which means these areas are suitable for surface water stagnation. 

The NDVI class (0. 1-0. 2) shows a high FR value (2. 45), indicating temperate and humid areas, consequently corresponding to areas of hydrogeological potential. In addition, the higher values of FR, 0. 12 and 1. 20, correspond to the surface temperature class's 21-27 • C and 27-32 • C, respectively. For the Shannon entropy model, the hydrological factors, such as slope degree, slope aspect, elevation, plan curvature, profile curvature, TWI, LS, and TPI, show descendent weights 0. 33, 0. 29, 0. 24, 0. 05, 0. 17, 0. 24, and 0. 06, respectively, which indicate that the slope degree is more suitable for hydrogeological exploration (Table 1 ). In addition, drainage density may be considered an important factor too, where it had a high weight value of 0. 30, and also a distance from the river that shows a value weight of 0. 20 (Table 1 ). On the other hand, lithology, lineaments, and the distance to lineaments factors have weight values of 0. 45, 0. 25, and 0. 16, respectively, indicating that lithology and lineaments are important indicators for groundwater reservoir mapping. Similarly, the NDVI and LST correspond to high weight values (0. 53 and 0. 24, respectively). 

Figures 13 and 14 show the FR and SE hydrogeologic potential maps, respectively, which are generated using Equations ( 16 ) and (17) based on all criteria found by the used models:

The FR groundwater potential map reveals that the area of very low, low, moderate, high, and very high groundwater potential area covers about 12. 6%, 22. 7%, 32. 3%, 21. 7%, and 10. 45% of the total study area, respectively (Figure 13, Table 2 ). While the SE groundwater potential map indicates that the areas of very high, high, moderate, low, and very low groundwater potential cover about 5. 73%, 13. 34%, 27. 09%, 33. 9%, and 19. 84% of the total study area, respectively (Figure 14, Table 2 ). Consequently, the favorable areas for groundwater resources exploration are double for the FR model compared to the SE model (10. 45% and 5. 73% for FR and SE, respectively), while the unfavorable areas are much more for the SE model than the FR model: 19. 84% and 12. 6%, respectively. 

The FR groundwater potential map reveals that the area of very low, low, moderate, high, and very high groundwater potential area covers about 12. 6%, 22. 7%, 32. 3%, 21. 7%, and 10. 45% of the total study area, respectively (Figure 13, Table 2 ). While the SE groundwater potential map indicates that the areas of very high, high, moderate, low, and very low groundwater potential cover about 5. 73%, 13. 34%, 27. 09%, 33. 9%, and 19. 84% of the total study area, respectively (Figure 14, Table 2 ). Consequently, the favorable areas for groundwater resources exploration are double for the FR model compared to the SE model (10. 45% and 5. 73% for FR and SE, respectively), while the unfavorable areas are much more for the SE model than the FR model: 19. 84% and 12. 6%, respectively. In this study, the hydrogeological mapping using the frequency ratio and Shannon's entropy models highlighted the importance of using the following factors as indicators of groundwater resources availability: slope, altitudes, lithology of geological formations, lineaments (faults), and distance from them. For future water drilling in the study area, both models confirmed the suitability of several zones, such as the Ameln, Oumaghouz, and Idikel valleys that drain the northern slope of Tikmirt and the valley of Tkatrt in the south of the Ameln Basin. The aforementioned valleys are generally small alluvial deposits of the Quaternary consisting of silt, scree, and gravel. In addition, the sedimentary outcrop formations of the Adoudounian and Georgian ages, essentially consisting of fractured limestones, dolomites, conglomerates, and shale, are moderately productive. Nevertheless, the shales are generally impermeable but can be a potential hydrogeological zone for the thick altered layers forming moderate slopes in depression In this study, the hydrogeological mapping using the frequency ratio and Shannon's entropy models highlighted the importance of using the following factors as indicators of groundwater resources availability: slope, altitudes, lithology of geological formations, lineaments (faults), and distance from them. For future water drilling in the study area, both models confirmed the suitability of several zones, such as the Ameln, Oumaghouz, and Idikel valleys that drain the northern slope of Tikmirt and the valley of Tkatrt in the south of the Ameln Basin. The aforementioned valleys are generally small alluvial deposits of the Quaternary consisting of silt, scree, and gravel. In addition, the sedimentary outcrop formations of the Adoudounian and Georgian ages, essentially consisting of fractured limestones, dolomites, conglomerates, and shale, are moderately productive. Nevertheless, the shales are generally impermeable but can be a potential hydrogeological zone for the thick altered layers forming moderate slopes in depression areas. In our study area, the shale plays the role of reservoir for the underflow of descending rivers of the Precambrian Water 2022, 14, 1553 20 of 26 massifs and boreholes along the valleys. In general, igneous rocks consist of rhyolites, andesites, granites, migmatites, and impermeable quartzite formations. Also, the eruptive rocks such as rhyolites, ignimbrites, and volcanic tuffs are impermeable but perhaps the groundwater potential reservoirs are found to be cracked. In our study area, all the water sources in these formations were found in fractured zones. The field survey confirms that the water sources, always of low flow, are distributed along with the contact of two different eruptive rocks (e. g., rhyolite and tuffs). In addition, these formations are marked by the presence of springs but they are always of low flow. On the other hand, at the surface of depression zones, the granites and migmatites formations are very altered, while they are not much altered in elevated areas (reliefs) where the network circulations are remarkable. In the hollows, the granites are covered by permeable arenas indicating aquifers. In other places, where the arenas are of low thickness, there are no generalized aquifers, but the alluvial aquifers that are drained by the rivers may be exploitable. Despite the slope steepness of the quartzite series of the Jbel Lkest group, it may be considered a suitable groundwater reservoir because it is very thick (more than 1000 m), is tectonized (fractured), and has fissure permeability. In addition, in the study area, the dolerites that cross the quartzites are of hydrogeological importance. In fact, at the low elevation points, the water sources were found in the quartzite-dolerite contacts forming a depression zone. In general, several springs identified in quartzites are of important flow (up to 20 L/s in the valley of Ameln at the foot of Jbel Lkest, near Tafraout city).
section_name: Frequency Ratio and Shannon Entropy Result

section_name: Validation of Groundwater Potential Map

section_content: Commonly, in the context of groundwater potential map validation, the used criterion in the evaluation must reflect the reality of the field, the criterion independence principles, and the conformity of the criterion [91, 92]. In this study, the visual interpretation of the overlapped water flow points (recorded from well points over the study area) on the obtained groundwater potential maps reveals that the percentage of the positive wells increases progressively toward the most productive areas (42. 86% for the FR model and 30. 61% for the SE model). In addition, we noted that few of these wells of low water flow are concentrated in unfavorable areas of hydrogeological importance (4. 08% for the FR model and 6. 12 for the SE model) (Figure 15 ). These findings confirm that Shannon's entropy model has a slightly better performance than the frequency ratio model in mapping groundwater potential areas. 

areas. In our study area, the shale plays the role of reservoir for the underflow of descending rivers of the Precambrian massifs and boreholes along the valleys. In general, igneous rocks consist of rhyolites, andesites, granites, migmatites, and impermeable quartzite formations. Also, the eruptive rocks such as rhyolites, ignimbrites, and volcanic tuffs are impermeable but perhaps the groundwater potential reservoirs are found to be cracked. In our study area, all the water sources in these formations were found in fractured zones. The field survey confirms that the water sources, always of low flow, are distributed along with the contact of two different eruptive rocks (e. g., rhyolite and tuffs). In addition, these formations are marked by the presence of springs but they are always of low flow. On the other hand, at the surface of depression zones, the granites and migmatites formations are very altered, while they are not much altered in elevated areas (reliefs) where the network circulations are remarkable. In the hollows, the granites are covered by permeable arenas indicating aquifers. In other places, where the arenas are of low thickness, there are no generalized aquifers, but the alluvial aquifers that are drained by the rivers may be exploitable. Despite the slope steepness of the quartzite series of the Jbel Lkest group, it may be considered a suitable groundwater reservoir because it is very thick (more than 1000 m), is tectonized (fractured), and has fissure permeability. In addition, in the study area, the dolerites that cross the quartzites are of hydrogeological importance. In fact, at the low elevation points, the water sources were found in the quartzite-dolerite contacts forming a depression zone. In general, several springs identified in quartzites are of important flow (up to 20 L/s in the valley of Ameln at the foot of Jbel Lkest, near Tafraout city).
section_name: Validation Using Wells Data

section_content: Commonly, in the context of groundwater potential map validation, the used criterion in the evaluation must reflect the reality of the field, the criterion independence principles, and the conformity of the criterion [91, 92]. In this study, the visual interpretation of the overlapped water flow points (recorded from well points over the study area) on the obtained groundwater potential maps reveals that the percentage of the positive wells increases progressively toward the most productive areas (42. 86% for the FR model and 30. 61% for the SE model). In addition, we noted that few of these wells of low water flow are concentrated in unfavorable areas of hydrogeological importance (4. 08% for the FR model and 6. 12 for the SE model) (Figure 15 ). These findings confirm that Shannon's entropy model has a slightly better performance than the frequency ratio model in mapping groundwater potential areas.
section_name: Validation Using Wells Data

section_content: The visual interpretation of the 7200 Hz coplanar apparent resistivity map of the geological formations of about 60 m in depth indicates the presence of four areas of hydrogeological importance (Figure 16 ). The first one corresponds to the alluvium and to the altered and fractured shales of the Ameln Valley, where the registered conductivity is less than 120 Ohm. m. This area is of high flow, where the water flow registered at the wells may be up to 10 L/s. On the contrary, the area in the north of the Ameln Valley shows more resistivity (>5000 Ohm. m), indicating an unfavorable area for hydrogeological exploitation. This area corresponds to the quartzite formations of Jbel Lkest, to shaleresisted alteration, or a very compact basic vein intrusion. In addition, the map shows that the second area of hydrogeological interest is located in the southwest part of Tafraout city, where a conductivity break appears. This area may correspond to a deep NE-SW fault located between two granitic bedrocks that have a resistivity of more than 5500 Ohm. m and show an average flow of 3 L/s. The third area corresponds to the Quaternary depression located in the north of Tikmirt. In this area, the apparent resistivity of the sediments covering resistant conglomerate of about 60 m of thickness is less than 100 Ohm. m. Further, due to the major and deep fault of Tighmi-Tizi Oumanouz oriented E-W (the southern equivalent of the Ameln fault), and due to the intense alteration of fractured shales, the geological formations located in the south of the Ameln Basin show a resistivity of less than 150 Ohm. m, indicating an area of hydrogeological potential that is favorable exploitation. This was confirmed by the flow rate of the aquifer productivity, which can reach 12 L/s. Overall, these findings highlight the importance of using the apparent resistivity data to confirm the developed hydrogeological maps of both models. Of note is that the areas of high conductivity indicate potential drinking water targets.
section_name: Validation Using Apparent Resistivity Data

section_content: The visual interpretation of the 7200 Hz coplanar apparent resistivity map of the geological formations of about 60 m in depth indicates the presence of four areas of hydrogeological importance (Figure 16 ). The first one corresponds to the alluvium and to the altered and fractured shales of the Ameln Valley, where the registered conductivity is less than 120 Ohm. m. This area is of high flow, where the water flow registered at the wells may be up to 10 L/s. On the contrary, the area in the north of the Ameln Valley shows more resistivity (>5000 Ohm. m), indicating an unfavorable area for hydrogeological exploitation. This area corresponds to the quartzite formations of Jbel Lkest, to shale-resisted alteration, or a very compact basic vein intrusion. In addition, the map shows that the second area of hydrogeological interest is located in the southwest part of Tafraout city, where a conductivity break appears. This area may correspond to a deep NE-SW fault located between two granitic bedrocks that have a resistivity of more than 5500 Ohm. m and show an average flow of 3 L/s. The third area corresponds to the Quaternary depression located in the north of Tikmirt. In this area, the apparent resistivity of the sediments covering resistant conglomerate of about 60 m of thickness is less than 100 Ohm. m. Further, due to the major and deep fault of Tighmi-Tizi Oumanouz oriented E-W (the southern equivalent of the Ameln fault), and due to the intense alteration of fractured shales, the geological formations located in the south of the Ameln Basin show a resistivity of less than 150 Ohm. m, indicating an area of hydrogeological potential that is favorable for exploitation. This was confirmed by the flow rate of the aquifer productivity, which can reach 12 L/s. Overall, these findings highlight the importance of using the apparent resistivity data to confirm the developed hydrogeological maps of both models. Of note is that the areas of high conductivity indicate potential drinking water targets.
section_name: Validation Using Apparent Resistivity Data

section_content: In Morocco, for many years, the Water Basin Agencies, responsible for water resources planning, management, and protection, tried to find an adequate approach to identify the most favorable zones for water exploitation to serve the population needs, especially in rural areas. The present study represents an approach anticipated to map groundwater potential areas in a fractured semi-arid area (Ameln Basin, Morocco). In this study, two groundwater potential area maps were prepared using the FR and SE models based on the geologic, topographic, hydrologic, and climatic data of the study area. The obtained maps were verified with wells and apparent resistivity data. The SE model reveals slope length, lithology, and slope degree that control the availability of groundwater in the study area more than the other factors. In contrast, for the FR model, the controlling factors are slope aspect, altitude, curvatures, and distance from lineament networks. In general, both models show promoting results in terms of identification of favorable areas for groundwater exploration, with a preference for the FR model that shows more percentage of surface area than the SE model (10. 45% and 5. 73% of the total study area for FR and SE, respectively). The favorable areas correspond to the Ameln and Oumaghouz valleys, altered surface depressions in the north of Tafraout, and the Quaternary formations located in the south of Tarçouat, in the north of Aoussift and Tikmirt, and in the southern part of the Ameln Basin. In addition, 42. 86% and 30. 61% of the total surface area are concordant to the positive wells (high water flow) for the FR and SE models, respectively, indicating a slightly better performance of the FR model in mapping groundwater potential areas than the SE model. Eventually, despite the limited number of wells in our study area used in the verification step, the obtained results allow us to conclude that both methods may be considered suitable tools for facilitating and orienting hydrogeological exploration strategies (water resources management, development, and exploitation). The conclusion is that both models can produce groundwater prospecting areas with very good accuracy. Future work should highlight the integration of mathematical approaches such as the boosted regression tree or advanced machine learning. The efficacy of the remote sensing methodology in groundwater delineation could be further improved by considering a sufficient quantity of information and data that have direct or indirect control over groundwater storage and occurrence. Indeed, to better appreciate the accuracy of these methods and to improve them, specific surveys are required in the most reliable and productive areas. In addition, the integration of aeromagnetic, gravimetric, and radar data will be helpful to enhance the lineaments map and then improve the accuracy of the used methods.
section_name: Conclusions

section_name: Data Availability Statement:

section_content: The data presented in this study are available on request from the corresponding authors.
section_name: Data Availability Statement:

section_name: Conflicts of Interest:

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: