Multi-Criteria Decision Analysis Coupled with GIS and Remote Sensing Techniques for Delineating Suitable Artificial Aquifer Recharge Sites in Tafilalet Plain (Morocco)

Article history: Received: 12 October, 2020 Accepted: 17 November, 2020 Online: 14 December, 2020 Despite that groundwater is an important and vital water resource, it is not well managed; depletion of aquifers around the world due to overexploitation is of serious concern especially in arid regions where the situation is much more alarming. Tafilalet plain in Morocco which belongs to this type of environment is certainly no exception and is viewing its groundwater disappearing. Artificial aquifer recharge (AAR) is found to be appropriate to such an urgent issue. Thus, the objective of this paper is the exploration of suitable sites to process Artificial aquifer recharge in Tafilalet plain by the joint use of remote sensing (RS), geographic information systems (GIS), and the analytic hierarchy process (AHP) method. For doing so, eight parameters were considered as groundwater influencing parameters such as slope, soil, geology, land Cover/land Use, depth to the water table, aquifer transmissivity, electrical conductivity and drainage density. The laters were integrated and processed in a GIS, their thematic layers were created, and their relative weights were generated by The AHP method based on their significance in recharging the aquifer. Afterward, thematic layers were reclassified and assigned their weights, so the GIS overlay tool was used for inferring artificial groundwater recharge potential regions in the study area. 47 % were identified as suitable while only 12 % were identified as unsuitable. Such studies facilitate groundwater management for stakeholders and water managers because important decisions may be taken in a record time which will preserve water resources and prevent them from being in an alarming situation.


Introduction
Groundwater is a valuable natural resource and constitutes a significant freshwater reservoir as it accounts for 30 % of the earth's freshwater while surface water accounts for just 0.3 % [1,2]. Besides, groundwater offers a multitude of benefits in different forms. In fact, groundwater meets the need for drinking water to relatively half of the world's population [3,4] and around 38 % of the world's irrigated lands are supplied by groundwater [5], not to mention its contribution to the industry sector and its paramount role for the environment. Nevertheless, there is growing evidence that groundwater is under threat of numerous difficulties such as climate change impacts, contamination, pollution in addition to depletion [6][7][8][9][10]. Actually, many aquifers throughout the world suffer from groundwater depletion [11][12][13][14][15][16][17] and this situation is much more serious in arid regions since groundwater constitutes their main source of water supply. Mismanagement and namely over-exploitation are considered as its major cause alongside many other factors. In fact, rapid population growth, agricultural and industrial activities that require a considerable availability of water together with climatic conditions notably of arid environments consisting of high evapotranspiration rates, low precipitation, and low vegetation cover are placing greater pressure on limited aquifer resources which has led to their intensive extraction in many parts of the world [9,12,[17][18][19][20][21][22][23][24][25]. Moreover, If this abstraction is greater than the natural groundwater recharge, it is referred to as an overdraft; so, groundwater levels will drop sharply over time, triggering depletion [8,9,14,17,[25][26][27][28][29][30], which can make the situation complicated as the replenishment of the aquifers ASTESJ ISSN: 2415-6698 to the original levels takes a lot of time [3] and it is more complicated in cases of fossil or compacting aquifers because of their irreversible and non-renewable character [14,31]. In such a way, countless negative impacts and many political, social, environmental, and economic issues, such as conflicts, land subsidence, degradation of groundwater dependent ecosystems, deterioration of water quality, aquifers drying and obviously water scarcity will take place particularly in countries which are groundwater dependent [10, 11, 14, 23-25, 27, 28, 32-39].
The above situation has been a major point of concern of water researchers, water managers, and water decision-makers, seeing that in many parts of the world a lot of measures have been taken to properly manage their groundwater resources and to prevent the worsening of the situation. The implementation of monitoring networks to control groundwater levels, the establishment of groundwater laws, and the setting up of exploitation rules for regulating its use, in addition to the organization of awarenessraising companies to urge on the importance of water and also education on water conservation are a good illustration of the undertaken measures [9,37,38,40,41], but even so, the decline of groundwater resources levels has been noticed and is still ongoing in many places around the globe. This is especially due to the violations of the established rules and the intensive mining of underground water notably in case of droughts. Added to that, the natural recharge of aquifers by streams, precipitation, rivers, pond seepage, reservoirs, irrigation canals, and lakes has become overwhelmed by the intensity of groundwater withdrawals [25] and has become difficult to occur especially in urbanized areas [42], because the presence of impervious surfaces such as roads, sidewalks, driveways and parking lots has led to many modifications to natural processes like impeding the infiltration of water. Thereby, the natural recharge of aquifers has proved to be insufficient to fill the quantity of water coming out of the aquifers in many parts of the world. Given these trends, water researchers, water managers, and water decision-makers have been prompted to think about solutions and innovative scientific analyses in order to bring better and sound management strategies that may address this groundwater issue and guarantee its sustainable use. If decreasing discharge and regulating abstraction has not brought important results, then this depletion can only be alleviated by increasing and boosting recharge.
Thus, artificial or managed aquifer recharge is a promising solution to this problem. Actually, artificial aquifer recharge (AAR) or managed aquifer recharge (MAR) is defined as an intentional human intervention arranged to enhance groundwater resources through processes implemented to heighten and improve the replenishment of depleted aquifers or depleting aquifers [22,[43][44][45][46][47]. The (AAR) has many gains and benefits, the major ones are the restoration of groundwater balance and the revival of the overexploited aquifers [11,19,25,39,48,49], the other ones are enhancement and raising the water table which will cut down its operating costs, improving groundwater quality and prevention of seawater intrusion in the case of coastal aquifers [28,46,49,50], alleviating floods adverse impacts given that during strong storms aquifers can be managed to play the role of a water storage reservoir which will mitigate the magnitude of flooding [22,45,[51][52][53], besides (AAR) can limit the loss of land covered with surface water storage structures and limit the loss of water by evaporation in dams, ponds, and other surface reservoirs by storing it in aquifers [12,39]. In order to artificially recharge aquifers, treated wastewater, floodwater, or water pumped from another aquifer can be used and redirected across the land surface. So, a variety of ways and systems have been developed such as infiltration basins, canals, ponds, ditches, sprinkler systems, rills, and irrigation furrows where the surface water is put on and subsequently infiltrates into the soil and moves downward to recharge aquifers or simply artificial recharge can be achieved by injecting water directly into the aquifer through injection wells especially in case of confined aquifers and each of these ways is implemented according to specific conditions. Applying the (AAR) needs some criteria to be met notably, the presence of highly permeable geological formations allowing percolation of recharge water, endowed with high capacity of water storage (unsaturated zone thickness should be deep enough to increase the piezometric level) and have high transmissivity. Moreover, the presence of a water source ensuring volumes of water needed for the recharge such as dams, stormwater runoff, wastewater treatment station, and desalination station besides, the water intended for artificial groundwater recharge as well as those of receiving aquifer must be of good quality and meet a certain number of standards and also good knowledge of the climatic characteristics of the region is essential in (AAR) projects (temperature, evaporation, rainfall, frequency of floods, periods of drought, snowmelt, etc.). Consequently, performing the (AAR) requires a suitable place whose selection depends on solid scientific analysis. Actually, selecting optimal sites for (AAR) involves integrating several parameters proved to have an influence on the infiltration of water to aquifers. To do this, a combination of geospatial technologies and multi-criteria decision analysis (MCDA) is used in this study.
Tafilalet plain was chosen as the study area in the present study as it is one of the Moroccan regions calling for urgent intervention and requiring important measures to be taken at all levels, including the environmental one; it suffers from enormous environmental difficulties mainly scarcity and poor water quality, due essentially to long periods of drought that has affected the region in recent years, combined to severe nature of the climate which makes groundwater resources the main source of water supply in low-flow periods. In 1967, the Hassan Addakhil dam was constructed in the north of Errachidia city and located at a distance of about 80 km from the Tafilalet plain, it was built on wadi Ziz in order to regulate floodwaters and, simultaneously, it serves a variety of goals, ranging from irrigation to water supply for domestic and industrial uses and remarkably it helps to establish a certain balance of groundwater of Tafilalet aquifer both quantitatively and qualitatively, but in front of considerable population growth and an agrarian-based economy resulting in increasing the number of wells notably prived ones, the exploitation of groundwater resources is proceeding at a very fast pace which has lead to their overdraft, so the implementation of artificial aquifer recharge in this region is essential to recover its overexploited groundwater.
Thus, the prime focus of the current paper is to delineate the artificial groundwater recharge potential regions in Tafilalet plain (Morocco), using the analytic hierarchy process (AHP) method which is a forme of Multiple-criteria decision analysis (MCDA) and geospatial technologies represented here by remote sensing (RS) and geographic information systems (GIS). After this introduction, the following parts are organized in five sections, the first one is an overview of the study area, the second one is a presentation of the materials and definition of the methods used, the third one is a description of the major results, the fourth one is a discussion about the outcomes and the last one is a conclusion of the main ideas of this study.

Study Area
Tafilalet plain belongs to the pre-Saharan region located in the extreme southeast of the kingdom of Morocco, particularly in the Drâa-Tafilalet region, at an altitude of 700 to 800 m, it lies in between the coordinates of projection system Lambert zone Morocco (Zone I) X = 580000 m to X = 630000 m and Y = 50000 m to Y = 120000 m (Figure 1), covers an area of about 700 km², and the majority of its population are rural. Besides, Tafilalet plain is part of the quaternary basin of Tafilalet; it is an erosional depression arising out of the secondary and tertiary coverage clearing during the quaternary era, the bottom of this depression was, afterward, covered by the alluvium of wadis Ziz and Rheris. In fact, it constitutes the zone of gathering and accumulation of alluvial deposits, conveyed by surges of the Ziz and Rheris wadis and their tributaries. Moreover, Tafilalet plain extends over three different geological formations, it is enclosed by the High Atlas in the north, it crosses the pre-Africain furrow and is surrounded at the south by the old massif of the Anti-Atlas domain as well as overlying an impermeable Palaeozoic and pre-Cambrian substratum.
Concerning the climate, Tafilalet plain is dominated by an arid to semi-arid climate, typified by a great rainfall scarcity as the average annual precipitation barely reaches 70 mm all over the area, this is mainly due to the presence of Atlas reliefs in the north which constitute a barrier to oceanic influences. The climate is also characterized by a great irregularity of the temperature, rising to very high values in summer (about 40 degree C ° especially in July) and dropping to very low values in winter (surroundings 0 degree C ° especially in January), this great difference in temperature brings out two seasons: a hot season from May to October and a cold season from November to April so, in the Tafilalet plain, The dry periods last almost all the year with extremely hot summers and highly cold winters which is mainly due to its openness on the Sahara. In addition, The mean annual evapotranspiration is close to 500 mm which is higher than the mean annual rainfall, this means that Tafilalet plain suffers from water deficit all over the year. Adding to that, the dominant winds are commonly known as Chergui, they are characterized by their hotness and dryness and they blow from the southeast to the northeast.
Regarding the water resources, Tafilalet plain is, on the one hand, underlain by confined and unconfined aquifers, the north of the region encloses the infra Cenomanian aquifer which is a totally confined aquifer, whereas the entire area is almost occupied by an unconfined, shallow and very extensive aquifer named as "Tafilalet aquifer" which is the aquifer targeted in this study. In fact, this aquifer is not very thick and the water below the ground surface is 2 m to more than 30 m in some places, its shape matches the slightly convex shape of the plain and its slope is very weak. Consequently, this makes the water of this aquifer sensitive to hydrological and atmospheric conditions such as pollution and evapotranspiration. The spatial distribution of Tafilalet aquifer transmissivities indicates dispersed values ranging from 10 -2 m²/s to 10 -4 m²/s and groundwater flows mainly from the north-east to south-west, south, and south-east characterized by an average hydraulic gradient that declines from upstream to downstream [54]. From a lithological point of view, the base of the aquifer is constituted by conglomerates and lacustrine limestones including gravelly levels and at the top of the aquifer, there are sands, gravel, and pebbles, covered with unevenly distributed powerful silts. Moreover, the main sources of groundwater recharge of the Tafilalet aquifer are essentially from the High Atlas mountains surrounding the plain in the north, from surface runoff including rivers, mainly Ziz and Rheris wadis, as well as irrigation water which infiltrates to the aquifer, and without forgetting precipitation and dams contribution. Conversely, the Tafilalet aquifer discharges mostly by evapotranspiration, by extraction from wells, and by Khettarats and springs. On the other hand, Tafilalet plain is crossed by many ephemeral wadis, the major ones are wadi Ziz in the center, wadi Rheris in the west. Besides, the region holds some springs resulting from the emergence of groundwater resources to the surface, such as Aïn Dolla-Attrouss and Aïn Chorfa.
Furthermore, the heart of the Tafilalet plain is made up of a very large palm grove, which occupies most of its surface so, a great part of the plain is devoted to agriculture, especially irrigated agriculture, making it the main economic activity in this region and is the main source of income for a large population. Because of climate aridity, the majority of crops are water-demanding which implies a strong mobilization of water notably during the months of June and July and due to the scarcity of surface waters in the region, growers are turning to groundwater exploitation to fill this water gap and meet their water needs. This has made the groundwater balance in aquifers of Tafilalet plain quickly in deficit.

Materials and Methods
The overall process followed in this study of the artificial aquifer recharge site selection has been schematically presented in (Figure 2).

Influencing factors Selection reason
The prediction of suitable sites for artificial aquifer recharge depends on a number of parameters that have an influence in either a direct or indirect way on the process of aquifer recharge. The selection of these factors was based firstly on a deep review of studies of artificial recharge carried out in other areas in the world, secondly on a complete understanding of the characteristics and specificity of the study area "Tafilalet Plain" and thirdly on the point of view of experts. Thus, eight factors were revealed to have control over the recharge of the Tafilalet aquifer, including slope, soil, geology, land cover/ land use, depth to the water table, aquifer transmissivity, electrical conductivity and drainage density as mentioned in (Table 1).
• Slope: The slope is one of the essential factors influencing the artificial aquifer recharge because the slope directly affects the penetration of water into the soil. Actually, the flatter the ground, the more time the water has to sink in, unlike the steep slopes where the water runs over the ground rather than entering it [25,[55][56][57]. However, the penetration of water is not guaranteed even if the slope is suitable because it depends on the soil's nature.
• Soil: Soil is a factor of significant importance in aquifer recharge because soil's nature determines the percolation ability of the soil as well as the water infiltration quantity and eventually governs the recharging capacity of the bedrock [58]. Moreover, soil acts as a purificator of water [29] that is why it should contain high rates of loam and sand to avoid groundwater contamination and to ensure a high infiltration rate. • Geology: Taking into consideration geology in artificial aquifer recharge site selection is of major importance because geological formations which must be selected should enjoy high permeability so that water can flow seamlessly downward towards the groundwater besides, they should have adequate transmissivity and high water storage capacity. • Depth to the water table: Depth to the water table or unsaturated zone thickness is the distance from the ground surface to the groundwater, it is another major factor to consider in artificial groundwater recharge site selection since, the deeper the water table, the greater the amount of groundwater storage, and the more efficient the artificial aquifer recharge will be Moreover, the deeper the water table, the more purified the water will be, because it will pass through a very large layer of soil, which is called soil aquifer treatment (SAT) [19,29]. Acceptable depth is about 3 m consequently, sites, where the depth is less than 3 m, should not be considered [59]. • Transmissivity: transmissivity is an important hydrodynamic parameter of an aquifer which makes its consideration paramount in the selection of suitable sites for artificial groundwater recharge. Actually, transmissivity determines the rate of water flowing per unit width of the fully saturated thickness of an aquifer under a unit hydraulic gradient [25,[60][61][62] so, to avoid the excessive formation of groundwater mounds that interfere with the infiltration process, aquifers should be unconfined and sufficiently transmissive to allow the lateral flow of infiltrated water well away from the recharge zone [49]. • Drainage density: Drainage density is defined as the total length of streams of all orders per unit area and provides an indication of the watershed draining state. Drainage density was heeded as another factor of artificial aquifer recharge site selection because it indirectly signals the suitability of a site. In fact, high drainage density allows less penetration and more surface runoff and areas with great infiltration capacity are endowed with low drainage density [25,39,51,58,[63][64][65][66][67][68] hence appropriate to artificial aquifer recharge. • Electrical Conductivity: Groundwater quality is a paramount parameter that should be investigated in the artificial aquifer recharge projects because it is a determinant factor of the future usage of that water. Usually, the water intended for artificial aquifer recharge is of good quality but when it penetrates across a geological formation and mingles with different dissolved materials in it, its quality may get altered. Electrical conductivity (EC) is considered in this study because its variations have similar trends to those of total dissolved solids (TDS) as shown in the following formula : In fact, the higher the electrical conductivity, the poorer the water quality and vice versa [69]. Thus, the quality of groundwater will be assessed through the EC in the artificial aquifer recharge site selection.
• Land Cover / Land Use (LCLU) : (LCLU) is one of the important aquifer recharge influencing factors, it affects both the infiltration capacity and the water restraining of the land. Actually, areas with an important density of vegetation cover are appropriate for undertaking artificial aquifer recharge projects because the denser the vegetation cover, the more the runoff is slowed down and the more water infiltration takes place and also the less evapotranspiration and soil surface erosion will occur. In addition, organisms present in the soil and plant roots form openings in the subsoil, allowing water to infiltrate easily, unlike in areas without vegetation and in built-up areas, where the soil is mostly impervious, which favors runoff and prevents infiltration.

Data extraction, preparation, and GIS processing:
Given the interdisciplinary nature of groundwater, using remote sensing (RS) and geographic information systems (GIS) techniques are required because RS is a technology that provides important and information-rich geospatial data, especially for unreachable and very large areas. Moreover, GIS is commonly known for its usefulness to handle, analyze, integrate, and process various thematic geospatial data. Therefore, in order to get the geographic data related to the above-mentioned factors, various data sources were used: FAO database, NASA database, USGS database, and other data were given by the hydraulic basin agency of Guir-Ziz-Rheris as detailed in Table 1. Furthermore, in order to process and to analyze those data, the geographic data processing functions of the ArcGIS software were used. In fact, as indicated in Table 1, slope and hydrographic network were extracted from ASTER DEM at 30 m, then the drainage density was formed by applying the Line Density function to the hydrographic network, the soil was extracted from the digital soil map of the world taken from the FAO database, geology was digitized from the geological map of Tafilalet-Taouz (1/200000 scale), electrical conductivity, transmissivity and depth to the water table were taken from field surveys carried out with the hydraulic basin agency of Guir-Ziz-Rheris and land cover/land use (LCLU) was classified from satellite image (Landsat 8) drawn from USGS database. Thus, a geospatial database was constructed and filled with the aforementioned data. Afterward, using ArcGIS software a thematic layer of each factor was created ( Figure 3) and each thematic layer has been converted to raster format and reclassified into five classes pertaining to its importance for artificial groundwater recharge (Figure 4), each class is translated to a category and each category has assigned a score as it is detailed in Table 2. Actually, the five categories adopted are very good, good, moderate, poor, and very poor in accordance with a specific evaluation scale (Table 3).

AHP method application
After reclassification of all the thematic layers, appropriate weights were generated by the analytic hierarchy process method (AHP).
Indeed, AHP is a multi-criteria decision analysis method and a mathematically-based technique [70], the mathematical underpinnings and the philosophy of which were developed in the seventies of the last century by the mathematician Thomas L. Saaty [71], AHP involves breaking down a convoluted multi-criteria decision problem into a hierarchy also called decision modeling in order to organize and better understand the problem and analyze the decision to be made.  After structuring the problem in such a way, the second step is to evaluate, compare and prioritize in a pairwise comparison matrix, based on decision-makers' judgments each criterion against the others pertaining to their importance in the problem decision making. To represent the magnitude of judgments, Saaty's fundamental scale should be used (Table 4), in doing so, relative weights of each criterion is deduced using Saaty's eigenvector technique and the generated weights must be validated because judgments are subjective and can typically be made with errors, so Saaty defined an index for measuring the degree of inconsistency which is the consistency ratio C. R. and it should be computed. Actually, C. R. is a way of measuring the number of "errors" created when making the judgments, the rule-of-thumb is that if the C.R. is below 0.1, the errors are relatively small, thus the final estimate is acceptable [70,[72][73][74][75][76].
The above-explained concept of AHP could be explained mathematically as follows: Suppose c 1 , c 2 , …, c n are n criteria influencing a decision with respect to a given problem, the pairwise comparisons matrix B=[b ij ] of that problem is an n×n matrix and it may be represented as follows : comparing the n criteria in pairs should comply with the following conditions: for all i and j which means that B should be a positive and a reciprocal matrix. Besides, the values that can be attributed to b ij for all i and j must be on the scale of intensities defined by Saaty (Table 4). Moderate importance one criterion is slightly favored over another 5 Strong importance one criterion is strongly favored over another 7 Very strong importance one criterion is highly favored over another 9 Extreme importance one criterion is extremely favored over another 2,4,6,8 Intermediate values when compromise is needed Reciprocals if criterion i has one of the above values assigned to it when compared with criterion j , then j has the reciprocal value when compared with i .
Deriving the weights w 1 ,w 2 , …, w n relative to the n criteria c 1 , c 2 , …, c n can be done by Saaty's eigenvector technique when judgments are inconsistent. Saaty's eigenvector technique stipulate that weights can be computed as the principal right eigenvector of the matrix B : where λ max is the maximum eigenvalue of the matrix B , or λ max for all i = 1, 2, …, n () After generating the weights, the consistency of the matrix should be measured by C.R. as explained above but first consistency index C.I. should be computed by the following formula: R.I. is the random index and is defined for each size of the matrix which is illustrated in Table 5. C.R. must be less than 0.1 to consider the result as acceptable, larger values require the decision-maker to revise judgments in order to reduce the inconsistencies.
Therefore, for the current study, application of the AHP gives firstly the hierarchy structure illustrated in Figure 5, secondly, the pairwise matrix (Table 6) was developed by comparing pairs of factors with the help of expert point of views and using the fundamental scale of Saaty (Table 4) to translate judgments into values. Comparisons were done with regard to the relative influence of each factor on aquifer recharge and on implementing the (AAR) technique in Tafilalet plain. R.I. used is 1.41 because the dimension of the pairwise matrix of the current study is 8×8 (Table 6) so, after computation, the obtained C.R. is equal to 0.041 which is less than 0.1 and according to the AHP method this result is satisfactory and the generated weights are acceptable.
Based on the generated weights, geology is the highest eigenvector while drainage density is the lowest ( Table 6). In fact, geology, soil type, depth to the water table, and transmissivity are the factors that were given the highest weights. This was expected because, according to the analysis that was made during pairwise comparisons, the first thing to look out for in a project for artificial recharge of an aquifer suffering from overexploitation is the geological formations, which must be sufficiently permeable and transmissive, the soil, which must have a high infiltration capacity and the groundwater which must be deep enough to provide sufficient space for water storage to allow the water table to regenerate its level.

AAR index generation
After generating the relative weights for each factor by the AHP method, the reclassified thematic layers of the factors were multiplied by their weights and finally added up using the weighted sum function of ArcGIS software.
Thus, the potential sites for the artificial aquifer recharge map was produced by calculating the AAR index according to the following formula: where: -I AAR is the AAR index -n is the number of factors influencing (AAR) which is eight in the current study -F i stands for factor -w i stands for weight

Result
The map of potential sites for artificial aquifer recharge which is illustrated in Figure 6 was produced as explained above and then reclassified according to five categories which are: very poorly suitable, poorly suitable, moderately suitable, suitable, and very suitable. It indicates that the most suitable areas for the artificial groundwater recharge of Tafilalet aquifer are located notably in the northeast and in the west of the plain. Statistically, About 47 % of the study area is suitable for artificial aquifer recharge and about 40 % are moderately suitable. Thus, those results can be used to formulate an efficient groundwater management plan for sustainable utilization of these limited resources because identified areas will be prioritized in the implementation of artificial aquifer recharge projects.
Based on the produced map of potential sites for artificial recharge of Tafilalet aquifer (Figure 6), most suitable zones are underlain by high permeable geological formations that mostly belong to the late quaternary and covered with soils with high infiltration capacity most often constituted by alluvium and basic aeolian reliefs rich in meteoric materials and located in areas where the water table is very deep more than 20 m deep in addition the slope of the plain is generally not too steep and the abovementioned description of the Tafilalet aquifer revealed that it is an alluvial and unconfined aquifer so, the artificial aquifer recharge technique most conducive with regard to the aforementioned conditions is infiltration basins nevertheless, the climatic characteristics of the region known by very high evapotranspiration rates will constitute a big challenge because more water will be lost by evaporation so less water will infiltrate unless the artificial aquifer recharge is conducted just during wet seasons when evapotranspiration rate decreases. Thus, several other factors must be taken into account when deciding whether a specific (AAR) system is suitable to be implemented, even if the appropriate locations are identified.

Discussion
The present-day situation of Tafilalet aquifer should urge water managers to work hard to avoid its worsening so, great efforts in risk mitigation need to be greatly accelerated and they should be centered on a range of scientific and engineering-based solutions. The current study presents one of those solutions which is the (AAR) and how to identify, in a scientific way, the most suitable locations for its implementation. In fact, A GIS-based multicriteria analysis was performed, using eight thematic maps in order to produce an appropriate map for the identification of sites for (AAR) in Tafilalet plain. The resulting map shows that the most suitable areas are identified in the northeast and in the west of the region while inappropriate locations are limited to small portions in the northern, southern, and eastern parts of the study area. This result is extremely important because it will help taking appropriate decisions for mitigating the overexploitation problem of the Tafilalet aquifer that is why the scientific approach should be taken into consideration in such important projects because it gives nearly accurate results and it brings time savings. The (AAR) will certainly sort out the problem of depletion of groundwaters in the Tafilalet aquifer as it was proved in some other aquifers in Morocco and in other countries. In fact, the (AAR) has already initiated and applied to some aquifers suffering from the same problem in Morocco such as the Charf El Akab aquifer in the region of Tangier, the Haouz Mejjate aquifer in the region of Marrakech, the Souss aquifer in the region of Agadir and the Berrechid aquifer located in the south of Casablanca, those aquifers have already started to reap its benefits. Even if the (AAR) practices in Morocco are still very much in their infancy, the technology is beginning to be recognized at the national level as having a value in tackling groundwater depletion issues.
Some attempts have been made in Tafilalet plain by the hydraulic basin agency of Guir-Ziz-Rheris in this regard and artificial recharging techniques have already begun to emerge on the plain, such as the construction of dikes or sills in the bed of the wadi Ziz in order to slow down the speed of the water and increase the wetted areas along the bed of the wadi, however, if the working field is extended to other locations that must be well chosen, these locations must be identified by means of scientific decisionmaking like GIS and multi-criteria analysis taking into consideration several criteria such as geological, pedological, hydrogeological, and many other factors to make sure that they are suitable for artificial recharge, in doing so, other artificial recharge techniques can be performed and can better and quickly restore the groundwater level situation of Tafilalet aquifer to its normal stage.

Conclusion
The aridity of the climate together with the excessive need for water is putting pressure on groundwater resources of Tafilalet plain so, the water level of the Tafilalet aquifer is declining to a threateningly low level which may undermine its productivity and causing its depletion, problems which are more likely to have irreversible adverse effects may occur like imperiling groundwater-dependent ecosystems, human activities as well as economic development of the region. Thus, the core objective of this study is to investigate suitable places for artificially recharge the Tafilalet aquifer to restore its levels to normal and fill the gap caused by overexploitation. For this purpose, thematic layers of factors that control groundwater recharge like depth to the water table, drainage density, electrical conductivity, slope, aquifer transmissivity, geology, soil, LCLU were prepared mapped, and reclassified in a GIS (ArcGIS software), after that, AHP method was applied to compute weights of each factor and the map of appropriate sites of artificial groundwater recharge was generated by the application of weighted sum tool of ArcGIS software. The result has revealed that around 47% of the study area has been identified as a high potential zone for groundwater recharge which is a very important outcoming. However, identifying suitable places is essential but not sufficient for implementing (AAR) because as mentioned above even if infiltration basins are estimated as the appropriate (AAR) system to implement in Tafilalet plain, there are other factors to consider like climatic conditions of the region, need for periodic maintenance of the basin in order to sustain infiltration rates over the long term because organic matter build-up, silting and/or clogging can make the surface of the basin impervious and the capacity of infiltration will decrease over time in addition to the existence of a good quality water source nearby that will ensure the supply of water to be put into the basin and many other conditions that if not met will render the project impractical. Lastly, the present study has shown an example of the benefits of coupling multi-criteria analysis with GIS and remote sensing techniques in taking relevant decisions to better manage and protect groundwater resources.