Hydrogeology and Groundwater Modeling of the Boudenib Senonian Aquifer (Southeastern Morocco)

Volume 1 - Number 1 - March 2017

EWASH & TI Journal, 2019 Volume 3 Issue 4, Page 221-230 Environmental and Water Sciences, Public Health & Territorial Intelligence Env.Wat. Sci. pub. H. Ter. Int. J. ISSN Electronic Edition : 2509 - 1069 Acces on line : http://revues.imist.ma/?journal=ewash-ti/

Section : Environmental and Water Sciences Publication type : Full Paper Hydrogeology and groundwater modeling of the Boudenib senonian aquifer (Southeastern Morocco) Hydrogéologie et modélisation hydrodynamique de la nappe sénonienne de Boudenib (Sud-Est Marocain) Received 27 Jul. 2019 JIHANE TLEMCANI1, MOHAMED HILALI1, ABDERRAHMAN MAHBOUB2, Accepted 25 Aug. 2019 HALA MOUDDEN2, ABDERAHMAN OUTALEB2 & SOUAD BOUMESHOUL2 On line 31 Dec. 2019 (1) Equipe de recherche de Génie des Mines, de l’Eau et de l’Environnement / Faculté des Sciences et Techniques / Université Moulay Ismail 52000, BP 509, Boutalamine, Errachidia, Maroc Email : [email protected] ; [email protected] (2) Agence du Bassin Hydraulique du Guir-Ziz-Rhéris Errachidia, Maroc Email : [email protected], [email protected], [email protected], [email protected]

KEY WORDS Abstract The region of Boudenib, located in the Guir basin in Southeastern Boudenib senonian aquifer, Morocco, contains, in terms of agricultural land, important hydraulic Hydrogeology, Modeling. potentialities and water resources more or less abundant by place. These resources make it possible to develop collective lands in the form of agricultural farms. The Guir basin contains two main aquifers: the Plio-Quaternaire alluvionnaire (water table) and the sandy and sandstone Sénonien (deep water table). The latter is considered to be the most important water table in arid areas in South-East Morocco. Nevertheless, it faces major natural and anthropogenic pressures. Due to its good hydrodynamic characteristics and accessibility, The aquifer of Senonia is of particular interest for its exploitation in the form of pumped samples. Numerous boreholes have been used to test this aquifer at varying depths. The aquifer reservoir of Sénonien has a diversified lithology, consisting of fine sands and red clays more or less sandy or gypsy, most often in sandstone beds with depths generally varying from 10 to 70 m. The hydrodynamic modelling of the Senegalese aquifer was carried out in order to follow the responses of this aquifer in relation to the scenarios of pumped water withdrawals for irrigation of large farms (palm groves). The aim is to have elements of decision-making support for the monitoring and accompaniment of agricultural investments at the level in the region of Boudenib.

MOTS CLES Résumé La région de Boudenib, située dans le bassin du Guir au Sud-Est Hydrogéologie, marocain renferme, en termes de terres agricoles, des potentialités hydrauliques Modélisation, Nappe importantes et des ressources en eau plus ou moins abondantes par endroit. Ces sénonienne de Boudenib. ressources permettent de mettre en valeur des terrains collectifs, sous forme de fermes agricoles. Le bassin du Guir renferme deux principaux aquifères : le Plio- Quaternaire alluvionnaire (nappe phréatique) et le Sénonien sableux et gréseux (nappe profonde). Cette dernière est considérée comme la plus importante nappe des zones arides située dans le Sud-Est marocain. Néanmoins, elle est confrontée à des pressions naturelles et anthropiques majeures.De par ses bonnes caractéristiques hydrodynamiques et son accessibilité, l’aquifère du Sénonien présente un intérêt particulier pour son exploitation sous forme de prélèvements par pompages.

J. TLEMCANI & al 221

De nombreux forages ont permis de tester cet aquifère à des profondeurs variables. Le réservoir aquifère du Sénonien présente une lithologie diversifiée, constituée de sables fins et d’argiles rouges plus ou moins sableuses ou gypseuses, le plus souvent à intercalations de bancs gréseux avec des profondeurs variant généralement de 10 à 70 m. La modélisation hydrodynamique de la nappe sénonienne a été effectuée afin de suivre les réponses de cette nappe par rapport aux scénarios de prélèvements d’eau par pompages pour l’irrigation des grandes exploitations agricoles (palmeraies). Le but est d’avoir des éléments d’aide à la prise de décision pour le suivi et l’accompagnement des investissements agricoles au niveau de la région de Boudenib.

1. Introduction Ce sont généralement des vents EW, SW et SSE. La vitesse La région de Boudenib, faisant partie des zones arides où le moyenne des vents atteint près de 3 m/s. problème de l’eau, souci séculaire, se pose avec le plus d’ampleur. Dans cette région, l’agriculture constitue la principale activité économique. Dans le cadre du Plan Maroc Vert, et à fin d’accompagner le processus de développement de la région par le biais des études de planification et de gestion des ressources en eau, la présente étude permet d’effectuer une modélisation de la nappe de Boudenib dans l’objectif de ressortir des éléments de base pour la prise de décision, par la réalisation d’un modèle hydrodynamique de la nappe de Boudenib. Présentation de la zone d’étude La zone d’étude, située dans le Sud-Est marocain, est limitée au Nord par le Haut Atlas, à l’Ouest par l’oued Guir, au Sud par les Hamadas du Guir et à l’Est par la palmeraie de Bouanane (Cf. Fig.1). La principale activité économique de la région est l’agriculture *1]. Fig. 1 : Situation géographique de la zone d'étude Source : Travail personnel Les précipitations moyennes interannuelles sont caractérisées par une variation aussi bien dans le temps que Le principal cours d’eau traversant la région de Boudenib est dans l’espace et ne dépassent guère 130 mm (Cf. Fig. 2). l’oued Guir. Les apports des eaux de surface proviennent Elles diminuent généralement du Nord vers le Sud. Les principalement des bassins haut-atlasiques sous forme de périodes pluvieuses s’enregistrent en octobre ou novembre, crues violentes et aléatoires [3]. Les apports observés depuis janvier-février et avril ou mai, tandis que la période sèche les années 60 à la station hydrologique de Tazouguert sont s’étale généralement sur 5 à 8 mois. La pluie mensuelle se caractérisés par une grande irrégularité interannuelle des caractérise par un régime pluviométrique très variable d’une débits reflétant le régime hydrologique pluvial de la zone année à l’autre. d’étude (Cf. Fig.3). Le débit moyen est de l’ordre de 1.2 m3/s, avec un débit minimal de 0.02 m3/s enregistré au cours des Les températures mensuelles interannuelles sont mois de juillet et aout, et un débit maximal de 40 m3/s relativement élevées et confirment l’aridité du climat, elles observé pendant le mois d’octobre. présentent d’importantes variations saisonnières avec un été très chaud et un hiver très froid. Les températures Le régime annuel des apports se caractérise par deux saisons moyennes annuelles sont de l’ordre de 21°C. L’humidité de hautes eaux (sous forme de crues), automne et atmosphérique est relativement faible toute l’année, la printemps séparées par deux saisons d’étiage, l’hiver et l’été moyenne mensuelle est de l’ordre de 30% [2]. L’évaporation pendant lesquelles, le débit est très faible et il est capté moyenne annuelle est de 3500 mm avec des valeurs totalement en amont des bassins versants. extrêmes en décembre et juillet. Les vents violents chauds soufflent essentiellement pendant la période s’étalant de mars à septembre. Les plus violents arrivent entre les mois d’avril et d’août.

J. TLEMCANI & al 222

300 250 200 150 100

50 Pluie en mm en Pluie

1983/1984 1998/1999 1971/1972 1974/1975 1977/1978 1980/1981 1986/1987 1989/1990 1992/1993 1995/1996 2001/2002 2004/2005 2007/2008 2010/2011 2013/2014 2016/2017 Années

Fig. 2 : Précipitations moyennes interannuelles observées au niveau de la station de Tazouguert Source : Travail personnel Fig. 4 : Esquisse géologique simplifiée de la zone d'étude Source : Travail personnel 9 8 La coupe ci-après (Cf. Fig.5) montre une succession 7 6 lithologique des formations essentiellement crétacées, avec 5 un approfondissement du Nord vers le Sud. 4 3

Débit en m3/s en Débit 2 1

1961/1962 1963/1964 1967/1968 1969/1970 1973/1974 1975/1976 1979/1980 1981/1982 1985/1986 1987/1988 1991/1992 1993/1994 1997/1998 1999/2000 2003/2004 2005/2006 2009/2010 2011/2012 2015/2016 1965/1966 1971/1972 1977/1978 1983/1984 1989/1990 1995/1996 2001/2002 2007/2008 2013/2014 Année

Fig. 3 : Evolution des débits moyens interannuels del’ouedGuir, enregistrés à la station de Tazouguert Source : Travail personnel

Géologie de la zone d’étude La zone d’étude fait partie du bassin crétacé de Boudnib, et correspond à une succession de plaines et de bassins qui Fig. 2 : Coupe géologique synthétique de la région de sépare le Haut Atlas au Nord de l’Anti-Atlas au Sud et ses Boudenib prolongements [4]. Source : [6] Sur le plan stratigraphique, le Plio-Quaternaire couvre la Caractérisation hydrogéologique de la région de partie centrale et sud de la région de Boudnib [5]. Il est Boudenib formé de limons, sable et alluvions avec une épaisseur variable ne dépassant guère les 25 m (Cf. Fig.4). La zone de Boudnib, faisant partie du bassin crétacé d’Errachidia, est caractérisée par la présence de deux Le Tertiaire continental couvre de grandes superficies. principales nappes, la nappe quaternaire et la nappe Il s’étend vers l’Est et le Sud pour former la Hamada du Guir. sénonienne. Il est représenté par des grès, des conglomérats, des argiles, des siltstones et des calcaires lacustres. Il comprend - Nappe du Plio-Quaternaire : quelques passées de gypse et d’anhydrite. L’aquifère est formé de conglomérats plio-quaternaires plus Le Sénonien affleure au Nord et au Sud de Boudenib. Il est ou moins cimentés et d’alluvions. L’épaisseur du constitué essentiellement d’argile relativement sableuse ou recouvrement plio-quaternaire varie entre 15 et 20 m. Le gréseuse et même gypsifère par endroits. substratum de la nappe de Boudenib est constitué de sables et d’argiles du Sénonien. Le niveau d’eau par rapport au sol Le Cénomanien-Turonien affleure dans la partie nord de la se situe généralement entre 3 et 5 m. Les transmissivités région et il est constitué de formations calcaires avec des sont de l’ordre de 10-2 m²/s. passages marneux et devient gréseux à sa base.

J. TLEMCANI & al 223

L’évolution piézométrique de la nappe met en évidence une - Nappe du Sénonien : stabilité relative (fluctuations faibles) des niveaux d’eau ; des Le réservoir aquifère du Sénonien présente une lithologie baisses qui demeurent faibles, et ne dépassent pas 1.5 m diversifiée, constituée de sables fins et d’argiles rouges plus (cas par exemple du piézomètre numéro 969/48) (Cf. Fig. 6). ou moins sableuses ou gypseuses, le plus souvent à La recharge de la nappe s'effectue principalement par intercalations de bancs gréseux. L’épaisseur du Sénonien est l’infiltration directe des précipitations, la réinfiltration des variable ; les faibles épaisseurs se situent au Sud de l’aire eaux d'irrigation et les eaux de crues. Les sorties de la nappe d’extension du Sénonien, et les épaisseurs maximales sont essentiellement représentées par les résurgences de (dépassant les 700 m) sont rencontrées au centre du sources, les khettaras (galeries souterraines captant la bassin.Lestransmissivités de l’aquifère varient entre 2 10-5 et nappe) et les pompages. 6.8 10-2 m²/s avec une moyenne de 1.5 10-2 m²/s. Les débits mis au jour sont variables avec une moyenne de près de La qualité chimique des eaux de la nappe est relativement 15 l/s [7]. La nappe du Sénonien est artésienne dans le bonne avec un résidu sec n’excédant pas 1.5 g/l. secteur situé entre Boudenib et Bouanane. La partie supérieure du Sénonien est un aquifère relativement perméable et en contact direct avec les aquifères plio- quaternaires. L’écoulement d’eau souterraine se fait essentiellement du NW vers le SE (Cf. Fig.7).

Fig. 3 : Exemple d’évolution piézométrique de la nappe plioquaternaire de Boudnib Source : Travail personnel

Fig. 4 : Piézométrie de la nappe sénonienne (Année 2018) Source : Travail personnel

J. TLEMCANI & al 224

Les fluctuations saisonnières des niveaux piézométriques du Dans ce cas on a choisi un maillage rectangulaire ayant Sénonien sont peu importantes ; elles varient généralement quatre nœuds qui permettent de définir une fonction H entre 0.5 et 1 m. La hausse des niveaux d’eau par endroit est (Potentiel Pondéré) en fonction des 4 nœuds constituant en relation avec une bonne recharge de la nappe à partir des chaque élément (Cf. Fig. 9). formations jurassiques perméables du haut Atlas (Cf. Fig. 8). Le principe est d’optimiser le résidu de l’équation de Poisson La qualité des eaux de la nappe sénonienne est sur tout le domaine de l’étude par intégration par partie en généralement bonne, à l’exception de certains secteurs utilisant la méthode de Galerkine [10], [11], [12], [13]. dominés par la présence des évaporites qui peuvent Le maillage utilisé pour cette étude est un maillage carré impacter négativement la qualité de ces eaux en particulier (500 m×500 m) (Cf. Fig. 10). les zones situées plus à l’Est de Boudenib en allant vers Bouanane. Ces zones ont une eau généralement assez minéralisée ; la conductivité électrique de l’eau est comprise entre 630 et 8600 µS/cm [8]. L’alimentation de la nappe du Sénonien se fait par abouchement latérale à partir du Jurassique du haut atlas situé au Nord de Boudenib. Elle peut se faire également par infiltration pluviale au niveau des zones d’affleurement *9+. Les prélèvements par pompage et le drainage souterrain constituent l’essentiel des sorties de l’aquifère.

Les ressources en eau souterraines renouvelables au niveau de la zone d’étude sont reparties comme suit : Fig. 6 : Maillage carré utilisé - Les réserves renouvelables au niveau de la nappe Source : [13] quaternaire de Boudnib sont estimées à 2.5 Mm3/an [3] ; - Les ressources renouvelables au niveau de la nappe crétacée sont estimées à 58.2 Mm3/an. Les ressources exploitées actuellement au niveau de la zone de Boudenib sont celles des nappes quaternaires et sénonienne. Par la suite, on traite le système des nappes quaternaire et sénonienne comme un aquifère bi-couche, avec des caractéristiques équivalentes, communément appelé dans notre cas la nappe de Boudenib.

Fig. 7 : Transformation de coordonnées d’un élément quadrilatéral(e), deux par deux en un élément carré Source : [13] La correction des frontières du maillage qui ne sont pas régulières a été faite de la manière suivante : Les équations utilisées pour la transformation des coordonnées d’un élément quadrilatéral (e) en un élément carré sont les suivantes : Fig. 5 : Exemple d’évolution piézométrique de la nappe sénonienne de Boudenib

Source : Travail personnel

Avec : 2. Matériel et méthodes 2.1. Modélisation hydrodynamique de la nappe Approche méthodologique de la modélisation La méthode choisie est celle des éléments finis, en fait c’est la méthode la plus adaptée à notre cas, car les frontières de la nappe de Boudnib ne sont pas régulières, le milieu ni homogène ni isotrope.

J. TLEMCANI & al 225

Et : BL : Recharge telle que : Avec A : aire de la maille, Q=débit sortant (-) ou entrant (+). Avec : Nous avons élaboré un code de calcul pour résoudre cette équation couplée avec VBA et SIG (version 10.1.4). 2.2. Discrétisation du domaine d’écoulement dans le système aquifère de Boudenib La nappe de Boudenib est généralement captive, mais elle Donc : devient libre par endroit. On peut considérer l’écoulement comme un écoulement plan, en considérant un écoulement verticalement hydrostatique quasi plan, l’équation de masse L’équation de l’hydrodynamique souterraine est la suivante : et la loi de Darcy peuvent être intégrées verticalement pour obtenir l’équation de Dupuit-Boussinesq en écoulement plan (x, y) puisque ∆H<

Avec : nappe (épaisseur supérieure à 200 m). Les variables d’état Tx : Transmissivité suivant X ; deviennent : la charge hydraulique moyenne H [m], et le flux Ty : Transmissivité suivant Y ; spécifique en 2D. La seconde variable provient de S : Coefficient d’emmagasinement ; l’intégration de la vitesse de Darcy. R : Recharge (source et puits). Pour discrétiser le système on tient compte de la densité et de la variabilité des données disponibles, selon les mailles Pour résoudre l’équation sur tout le domaine, on a utilisé la carrées régulières de l’ordre de 500 m /500 m disposées sur méthode de Galerkine qui est une classe de méthodes 46 lignes et 101 colonnes (Cf. Fig. 11). permettant de transformer un problème continu (par exemple une équation différentielle) en un problème discret. Cette méthode est couramment utilisée dans la méthode des éléments finis. Elle consiste à utiliser un maillage du domaine d'étude, et considérer la restriction de la fonction solution sur chacune des mailles.

Les équations utilisées sont : : Fonction d’interpolation :

Pour le cas particulier d’un l’écoulement unidimensionnel :

peut être résolue par intégration par Fig. 8 : Maillage de la nappe de Boudenib partie comme suit : Source : Travail personnel Posant : et 2.3. Conditions aux limites La synthèse hydrogéologique de la zone d’étude a permis de Alors : déterminer la nature des limites du système aquifère afin de L’application de l’écoulement bidimensionnel les intégrer au sein du modèle. Ces conditions aux limites se donne : répartissent comme suit : - Du côté Nord on a un flux souterrain provenant du Haut Atlas au Nord de la nappe en plus de celui du barrage Kadoussa sur la frontière Nord-Ouest ; - Le flux aux autres côtés est nul ; - Le flux entrant au niveau de la partie Nord de la zone d’étude est estimé à 45 Mm3/an ; F : fonction flux (conditions aux limites) - Un flux entrant à partir du barrage Kadoussa au niveau de Donc : la première maille du modèle de la partie Nord-Ouest de la zone d’étude est calculé en utilisant le coefficient Finalement : d’infiltration de 2.5% et le volume des eaux des lâchers à 3 Donc : partir du barrage Kadoussa est estimé à 26 Mm /an.

J. TLEMCANI & al 226

2.4. Recharge par les précipitations Les sorties d’eau de la nappe sont comme suit : Le potentiel pluviométrique de la zone d’étude est établi sur - Le débit moyen des khettaras existantes dans la zone la base des moyennes interannuelles des séries actualisées d’étude qui est de 140 l/s ; des deux stations proches (Tazouguert et Bouanane). Nous - Le débit moyen de la source de Meski situées à l’Ouest de avons procédé au découpage de la zone d’étude en deux Boudenib qui est de 60 l/s ; parties, les précipitations dans la partie Nord dans de la zone - L’évaporation et l’évapotranspiration au niveau de la zone d’étude sont de 130 mm et les précipitations dans la partie sont calculées par les formules de Turc et Thornthwaite, à Sud sont de 120 mm. Le volume de la pluie infiltrée partir des données mensuelles des précipitations et des directement dans la zone d’étude est calculé moyennant un températures. coefficient d’infiltration de 2.5%. 3. Résultats et discussion 2.5. Les prélèvements 3.1. Calage et validation du modèle Les prélèvements par pompage au niveau de la zone d’étude Le calage du modèle a été réalisé sur la base des données jusqu’en 2013 sont de l’ordre de 16.7 Mm3/an. Pour les piézométriques observées en 2009. Ce calage a permis prélèvements pendant la durée allant de 2014 à 2017, ils d’ajuster la transmissivité à 1.3 10-2 m2/s. sont de 15.3 Mm3/an, déduits des autorisations de prélèvements délivrées par les services de l’Office Régional La validation a été faite en comparant les résultats du de Mise en Valeur Agricole de Tafilalet sur la base des débits modèle à ceux observés en 2017, on constate que les fixés après des essais de pompages réalisés par les résultats obtenus par le modèle sont satisfaisants, ce qui pétitionnaires concernés. permet de l’utiliser sans réserves pour simuler différents scénario (Cf. Fig. 12).

Fig. 9 : Carte de validation du modèle par les données piézométriques observées Source : Travail personnel

3.2. Application du modèle 3ème scénario correspondant à la superficie autorisée et non Maintenant que le modèle est validé, on va l’utiliser pour encore aménagée ni plantée qui est de l’ordre de 24000 ha ; simuler les scénarios suivants : 4ème scénario est celui qui prévoit la baisse des flux entrants 1er scénario correspondant à une superficie irriguée de à partir du Haut Atlas et de l’infiltration des eaux des lâchers 2500 ha ; du futur barrage de Kadoussa pour une superficie irriguée de 24 000 ha, cette baisse est respectivement de 10%, 20% et 2ème scénario correspondant à la superficie autorisée et 40%. aménagée et non encore plantée qui est de l’ordre de 4800 ha ; J. TLEMCANI & al 227

La hauteur d’eau est calculée par rapport au substratum de - 2ème scénario : la nappe qui se situe à 250 m par rapport au sol. Les cartes ci-dessous (Cf. Fig. 15.16) représentent les cartes piézométriques au niveau de la nappe de Boudnib pour une Les résultats des différents scénarios sont présentés dans ce surface irriguée de 4800 ha pendant des durées d’une et de qui suit : cinquante (50) années respectivement. - 1er scénario : On constate que la baisse maximale du niveau d’eau pour Les cartes ci-dessous représentent une comparaison de la ces périodes est presque un mètre. variation du niveau piézométrique de la nappe de Boudnib pour une surface irriguée de 2500 ha pour des durées d’une et cinquante (50) années respectivement (Cf. Fig. 13.14). Pour ce cas la baisse maximale du niveau piézométrique pour ces périodes n’est pas significative, elle est de l’ordre de quelques centimètres. Il est à noter que l’apparition du dôme piézométrique positive sur les cartes est due au flux entrant du lâcher du barrage Kadoussa.

Fig. 12 : Piézométrie de la nappe de Boudenib pour une surface irriguée de 4800 ha pendant une durée d’un an Source : Travail personnel

Fig. 10 : Piézométrie de la nappe de Boudenib pour unesurface irriguée de 2500 ha pendant une durée d’une Fig. 13 : Piézométrie de la nappe de Boudenib pour une année surface irriguée de 4800 ha pendant une durée de 50 ans Source : Travail personnel Source : Travail personnel

- 3ème scénario : Les cartes ci-dessous (Cf. Fig. 17.18) représentent les cartes piézométriques au niveau de la nappe de Boudnib pour une surface irriguée de 24000 ha pendant des durées d’une et cinquante (50) années respectivement. On révèle que la baisse maximale du niveau d’eau pour ces périodes est de l’ordre 5m.

Fig. 11 : Piézométrie de la nappe de Boudenib pour une surface irriguée de 2500 ha pendant une durée de 50 ans Source : Travail personnel

J. TLEMCANI & al 228

Fig. 16 : Piézométrie de la nappe de Boudenib pour une surface irriguée de 24000 ha avec une baisse du flux de 10%

pendant une durée de 50 ans Fig. 14 : Piézométrie de la nappe de Boudenib pour une Source : Travail personnel surface irriguée de 24000 ha pendant une durée d’un an Source : Travail personnel Baisse du flux de 20% : On constate que la baisse maximale du niveau d’eau pour une période de 50 ans est de l’ordre de 15 m.

Fig. 17 : Car Piézométrie de la nappe de Boudenib pour une Fig. 15 : Piézométrie de la nappe de Boudenib pour une surface irriguée de 24000 ha avec une baisse du flux de 20% surface irriguée de 24000 ha pendant une durée de 50 ans pendant une durée de 50 ans Source : Travail personnel Source : Travail personnel ème - 4 scénario : Baisse du flux de 40% : Le quatrième scénario est celui qui prévoit la baisse des flux La baisse maximale du niveau d’eau pour une période de 50 entrants à partir du Haut Atlas et de l’infiltration des eaux ans est de l’ordre de 25 m. des lâchers du barrage de Kadoussa, cette baisse est respectivement de 10%, 20% et 40%. Les cartes isopièzes L’apparition du dôme piézométrique négatif sur les cartes ci- représentatives de chaque cas sont présentées ci-dessous dessus est due aux prélèvements au niveau de la région de (Cf. Fig. 19. 20. 21). Boudenib. Baisse du flux de 10% : On observe que la baisse maximale du niveau d’eau pour une période de 50 ans est de l’ordre de 5.5 m.

Fig. 18 : Piézométrie de la nappe de Boudenib pour une surface irriguée de 24000 ha avec une baisse du flux de 40% pendant une durée de 50 ans Source : Travail personnel J. TLEMCANI & al 229

4. Conclusion le bassin sud-atlasique : Cas du Sénonien de la plaine de Boudenib (Sud du Maroc), 2011. La nappe sénonienne de Boudenib est renfermée dans des [3] Agence du Bassin Hydraulique du Guir-Ziz-Rhéris, sables et grès, la transmissivité moyenne est de l’ordre de Actualisation du plan directeur d'aménagement intégré 1.3 10-2 m²/s. L’écoulement général de la nappe se fait du des ressources en eau des bassins de Guir - Ziz - Rhéris NW vers le SE ; dans tous les cas l’oued Guir draine la nappe. et Maider. Rapport de synthèse, 2011. L’évolution piézométrique de la nappe montre une relative [4] J.Hindermayer, Observations géologiques dans les stabilité du niveau d’eau dans le cas où le flux ne change Hamadas au Sud de Boudenib. Tome 3, Notes et pas. Mémoires du Service Géologique du Maroc, 1950, n°76, P. 105-139. Les résultats de la modélisation hydrodynamique ont [5] A.Jossen, J.Filali, Le versant sud du bassin d’Errachidia : montré une baisse des niveaux d’eau au niveau des points Synthèse stratigraphique et structurale ; Contribution à de prélèvements de la nappe de Boudenib pendant une l’étude des aquifères profonds. Rapport de la Direction période de 50 ans qui est de l’ordre de 25 m, dans le cas de la Recherche et de la Planification de l’Eau-Rabat, d’une superficie irriguée de 24000 ha et une baisse du flux 1989. du 40%. Il est à noter que les baisses maximales du niveau [6] Campagnie Africaine de Géophysique, Prospection d’eau constatées dans tous les scenarios étudiés se situent géophysique par méthode électrique dans le bassin au niveau de la zone sollicitée par les prélèvements par d’Errachidia-Boudnib. Rapport de synthèse, 1988. pompages. [7] M. Hilali, A.Mahboub, A.Kettani, M.Boualoul, H. Sahbi, T.Bahaj, I.Kacimi, Bassin crétacé au Nord d’Erfoud : les Si on comptabilise les eaux de surface au niveau du bassin nouveaux acquis des prospections géophysiques et du Guir qui sont de l’ordre de 59 Mm3/an en plus des eaux hydrogéologiques. Notes et Mémoires du Service souterraines, les investissements agricoles peuvent Géologique du Maroc n°577. 2014. Editions du Service atteindre plus de 20 000 ha sans porter préjudice à la nappe Géologique du Maroc-Rabat, pp : 147-154. sénonienne de Boudnib. [8] A. El Ouali, Contribution du Haut Atlas central au Sud de Midelt à l’alimentation des aquifères profonds du Ces résultats ont été confirmés par l’observation des bassin crétacé d’Errachidia. Thèse de doctorat es mesures du niveau d’eau dans les piézomètres contrôlant la sciences de la terre. Université de Franche-Comté. nappe sénonienne de Boudenib. Cette dernière enregistre France, 1992. une hausse du niveau piézométrique qui prouve que la [9] T. Bahaj , I.Kacimi , M .Hilali, N.Kassou & A .Mahboub, nappe est bien alimentée par les eaux du l’aquifère Preliminarystudy of the groundwater geochemistry in jurassique du Haut Atlas situé en amont hydraulique du the sub-desert area in Morocco: Case of the Ziz-Rheris bassin de Boudenib et que les prélèvements au profit des basins. Volume 7 Sciverse Science Direct. Procedia palmeraies sont renouvelés. Earth and Planetary Science. Elsevier, ISSN 1878-5220. Contribution des auteurs Editors: Roland Hellmann and Helmut Pitsch, 2013 p. 44-48. Le travail est réalisé par une équipe de l’Agence du Bassin [10] CA. Brebbia, J.J. Connor, Finite Element Techniques for Hydraulique du Guir-Ziz-Rhéris sous la direction de fluid flow. London, Boston : New nes-Butterworths, Monsieur Abderrhaman Mahboub, et un groupe de 1977, 310p,. doctorants sous la direction de Monsieur Mohamed Hilali, [11] C.R.Faust, J.W.Mercer, Ground water modeling: Recent affilié à l’Equipe de Recherche de Génie des Mines, de l’Eau developments. Ground Water, 1980, 18 (6): 569-577. et de l’Envrionnement de l’université Moulay Ismail. [12] P.W.Huntoon, Finite Difference methods as applied to Remerciement the solution of groundwater flow problems. Laramie : Wyoming water resources research institute, 1974, Les auteurs souhaitent remercier les services de l’Agence du University of Wyoming, 108p. Bassin Hydraulique du Guir-Ziz-Rhéris d’Errachidia pour leur [13] H.F.Wang, M.P.Anderson, Introduction to groundwater aide fructueuse et disponibilité de données. modeling; Finite Difference and Finite Element Références bibliographiques methods. University of Wisconsin, Madison. Academic [1] Agence du Bassin Hydraulique du Guir-Ziz-Rhéris, Etude Press Inc. California, 1982, 246p. de synthèse hydrogéologique du bassin crétacé d’Errachidia. Rapport de la mission 1 : Collecte, actualisation des données et élaboration d’un SIG, 2008. [2] M. Hilali, M.Boualoul, H.Sahbi, et T Bahaj, Potentialités hydrogéologiques des réservoirs aquifères situés dans

J. TLEMCANI & al 230