Corps de l’article

1. Introduction

L’utilisation des eaux usées brutes en agriculture est une pratique répandue dans les régions du monde affectées par la pénurie en eau (EL KETTANI et al., 2006). Au Maroc, elle est réalisée d’une manière courante et très ancienne à la périphérie des grandes villes continentales. Elle s’est développée ces dernières années, autour de certaines agglomérations pourvues d’un réseau d’assainissement (JEMALI et KEFATI, 1994). Cette pratique peut être très bénéfique pour la production agricole, car elle est économique, permet d’irriguer les cultures, d’épurer l’eau et de recycler des éléments nutritifs (DERE et al., 2006). Toutefois, elle peut également entraîner des problématiques importantes de pollution de l’eau et des sols, qui peuvent engendrer des risques pour la santé humaine (LAAMARI et al., 2004).

La région de Settat au Maroc, comme plusieurs régions arides en Afrique, fait face à cette problématique depuis plusieurs années. L’oued Boumoussa qui est le principal cours d’eau de la région reçoit des eaux usées domestiques et industrielles de la ville de Settat qui sont évacuées sans traitement préalable dans l’oued Boumoussa. Le long de son passage, les riverains utilisent l'oued pour l’irrigation des cultures céréalières et fourragères. Les polluants contenus dans ces eaux peuvent donc s’accumuler dans les sols et l'aquifère au niveau de la zone au nord de la ville de Settat. Compte tenu de ceci, il est clair que les risques de dégradation des composantes de l’écosystème et de la santé humaine sont élevés.

Cette problématique a été étudiée dans le cadre du projet Écosystème et santé humaine (financé par le Centre de recherche pour le développement international du Canada et la Fondation Ford des États-Unis). Cette étude (2002-2005) a été menée sur les eaux souterraines et les eaux superficielles de la communauté de Mzamza (plaine de Berrechid) dans le but de déceler l’impact de la réutilisation des eaux usées pour l'irrigation sur l’état de santé de la population. Étant donné l’impact négatif sur la santé de la population mis en évidence dans ce projet (LAAMARI et al., 2004), une station d’épuration et de traitement des eaux usées domestiques par lagunage naturel (STEP) a été mise en place par la Régie Autonome de Distribution d’Eau et d’Électricité de la Chaouia (RADEEC) en partenariat avec d’autres organismes nationaux et internationaux et est opérationnelle depuis mai 2007.

En 2009, un projet de recherche entre la Faculté des Sciences et Techniques de Settat et l’Institut de l’Eau de l’Université de Grenade (financé par l’Agence espagnole de coopération internationale pour le développement) a porté sur la contamination des sols et des eaux souterraines par les eaux d'irrigation chargées en contaminants industriels. En 2010, un échantillonnage a été réalisé au niveau des puits, de l'oued Boumoussa et des sources situées sur le plateau de Settat (en amont sur le bassin de l’oued Boumoussa) pour comprendre l’origine de l'eau souterraine qui s'écoule à Settat.

L'ensemble des études réalisées depuis 2002 a permis d'obtenir une meilleure compréhension du fonctionnement hydrodynamique de l'aquifère de Settat, mais aucune nouvelle mesure de gestion de l’eau n’a encore été mise en place. L’objectif de ce travail était de réaliser une synthèse des données physicochimiques et isotopiques afin de proposer un modèle conceptuel de circulation des eaux au niveau de l’aquifère plioquaternaire superficiel de Settat au Maroc, notamment pour déceler l’origine des apports à l’aquifère provenant de la surface et de l’aquifère sous-jacent. Cette étude pourra guider les gestionnaires vers une meilleure gestion de l’eau dans la région.

2. Description de la zone d’étude et méthodologie

2.1 Région d'étude

La région d’étude couvre une superficie de 80 km2 au sud de la ville de Casablanca au Maroc (Figure 1). À la sortie de la ville de Settat, l’oued Boumoussa traverse, sur une distance de 8 km, deux localités de la communauté de Mzamza, soit Oulad Boukhalou et Dladla. L’oued Boumoussa est une rivière intermittente et de faible débit, il prend sa naissance en amont dans le plateau de Settat et vient se dessécher à son aval dans la plaine de Berrechid. Lorsqu'il est actif, il s’écoule sur un trajet de 25 km du sud vers le nord, traversant la ville de Settat en longeant la route principale Casablanca-Marrakech.

Figure 1

Situation géographique de la zone d’étude et localisation des points de mesure

Geographic location of the study area and the sampling points

Situation géographique de la zone d’étude et localisation des points de mesure

-> Voir la liste des figures

Le climat est semi-aride méso thermique (25 °C en moyenne) et reçoit une pluviométrie de l’ordre de 325 mm∙a-1 (EL ASSAOUI, 2017). L’évapotranspiration potentielle moyenne est de l’ordre 1 182 mm∙a-1 (LAAMARI et al., 2004). La saison des pluies s’étend de novembre à mai.

La région est une zone charnière entre le plateau de Settat au sud et la plaine de Berrechid au nord. Cette limite plateau de Settat-plaine de Berrechid témoigne d’une tectonique plioquaternaire par l’existence de failles importantes. Les deux ensembles qui constituent le plateau de Settat et la plaine de Berrechid sont formés de roches sédimentaires sur lesquelles vient se superposer la couverture quaternaire récente. La plaine de Berrechid se caractérise par des affleurements de limons argileux quaternaires ou de grès calcaires pliocènes. Le plateau de Settat se situe dans une cuvette d’âge cénomanien. Les formations marno-calcaires du Cénomanien, affleurant largement sur le plateau de Settat, contiennent une nappe d'eau de faible extension participant en partie à l’alimentation de la nappe plioquaternaire de Berrechid. Leur rôle hydrogéologique est important à l'amont de la plaine où des circulations de type karstique peuvent constituer un relais entre les formations profondes du plateau de Settat et l'aquifère plioquaternaire de la plaine de Berrechid.

La nappe de Berrechid se caractérise au fond de la cuvette de Settat par des niveaux piézométriques élevés qui contribuent à la vulnérabilité de cette ressource aux différentes pollutions de surface. Les limons quaternaires forment la quasi-totalité du recouvrement de la plaine. Peu perméables, ils jouent un rôle hydrogéologique important du fait qu'ils rendent captive la nappe des formations plioquaternaires et qu’ils sont localement à l'origine de petites nappes perchées. À petite échelle, l'absence de flaques d’eau de grande dimension (daya) sur la plaine montre toutefois que la conductivité hydraulique de ces limons est suffisamment importante pour permettre une infiltration des eaux de pluie et donc une alimentation de la nappe à partir de la surface (EL BOUQDAOUI, 1995).

Il ressort de l’analyse de la piézométrie relative à la période de septembre 2003 à octobre 2004 sur l’ensemble des puits que les variations de niveaux piézométriques sont faibles durant cette période (l’écart moyen varie entre 0,20 et 0,68 m). L’alimentation de la nappe semble conditionnée par un apport relativement constant à partir des précipitations (en hiver) et à partir du retour d’irrigation (en été). L’analyse de la piézométrie pour le mois de janvier 2004 (période intermédiaire entre les hautes eaux et l’étiage) (Figure 2) montre que la direction générale d’écoulement de la nappe est sensiblement SE-NO dans le sud de la zone d’étude avec deux dépressions piézométriques au niveau des localités Oulad Boukhalou et Dladla, qui témoignent de l’effet des pompages pour l’irrigation des cultures (EL AMRANI PAAZA et al., 2008).

Figure 2

Carte piézométrique de janvier 2004 (EL AMRANI et al., 2008). Carte en coordonnées Lambert; les lignes noires sont des courbes isopièzes (m snm)

Piezometric map of January 2004 (EL AMRANI et al., 2008). Map in Lambert coordinates; black lines correspond to piezometric levels (m asl)

Carte piézométrique de janvier 2004 (EL AMRANI et al., 2008). Carte en coordonnées Lambert; les lignes noires sont des courbes isopièzes (m snm)

-> Voir la liste des figures

À hauteur de la ville de Settat, l’oued Boumoussa est canalisé sur 3 km. Sur cette distance, il collecte les eaux d’orages, les eaux usées de la ville (actuellement canalisées en partie, et évacuées vers la STEP) et celles de la zone industrielle. Les unités industrielles doivent toutefois se doter de systèmes de prétraitement pour que leurs rejets aient des teneurs en éléments chimiques pouvant être traités à la STEP et avoir l’autorisation de se brancher dans le réseau d’assainissement municipal. Les établissements industriels de Settat contribuent actuellement au débit total des eaux usées par temps sec de la ville à hauteur de 25 %, ce qui correspond à une proportion importante de la charge polluante globale générée (RADEEC, 2007). Les débits jaugés à l’oued varient de 20 L∙s-1 en basses eaux jusqu'à près de 600 L∙s-1 en période de crues et en temps de rejets. Cette évolution des débits est accompagnée par une hausse de la salinité de l’eau de rivière; cette salinité des eaux étant fonction du type de rejets, soit domestiques ou industriels. La STEP permet de traiter un volume d’eau estimé à 9 632 m3∙j-1 (RADEEC, 2007). La mise en oeuvre du lagunage naturel en aval de la ville de Settat a permis de mettre en conformité les effluents de la ville par rapport aux normes de rejets (OMS, 1989), et éventuellement les réutiliser pour les besoins de l’agriculture.

En amont, le parc industriel de Settat comprend essentiellement des usines de textile, une tannerie, des abattoirs et des verreries-cristallerie. On trouve également des usines de fabrication de produits céramiques, de marbre, des plastiques, des équipements électriques et électroniques, des minoteries, des usines de production de margarine, un abattoir de volaille, des boyauderies, une usine de fabrication de clous et une imprimerie (RADDEC, 2006). Les eaux usées industrielles ne peuvent être traitées à la STEP du fait de leur salinité très élevée (conductivité électrique des effluents varie entre 9 000 et 23 000 µS∙cm-1) qui pourrait endommager les installations et inhiber la fonction biologique des bactéries responsables de l’épuration. De même, des teneurs élevées en éléments traces métalliques (Cu, Cr, PB, ZN, Ni, Fe, Cd et Al) émanant de ces activités industrielles ont été détectés au niveau des eaux usées, eaux souterraines et sols (EL AMRANI PAAZA et al., 2010). Le détournement des eaux usées domestiques vers la STEP a entraîné une diminution significative du débit de l’oued Boumoussa qui s’assèche maintenant à 3 km à l’aval de la ville. Seule la localité d'Oulad Boukhalou peut maintenant utiliser l'eau de l'oued pour l’irrigation. Plus en aval de la zone d’étude, la localité de Dladla doit maintenant utiliser l'eau souterraine pour irriguer; ce recyclage des eaux souterraines chargées en nitrates induit une contamination plus accrue de l'aquifère plioquaternaire par les nitrates.

2.2 Échantillonnage et analyses

Deux périodes d'échantillonnage ont été considérées dans la présente étude, celles de juin 2009 et de novembre 2010. Ces campagnes d’échantillonnage ont permis d’échantillonner 21 stations, c.-à-d. cinq points dans l’oued Boumoussa, 12 puits privés utilisés pour l’approvisionnement et l’irrigation (huit dans le Plioquaternaire et quatre dans le Cénomanien), trois sources (Oum Karma, Hlilifa et Baida) et l’oued Zokch (voir figure 1 pour la localisation des stations). L'échantillonnage de juin 2009 a consisté à prélever, sur la plaine de Berrechid, 14 échantillons d’eau (dix au niveau des puits, trois au niveau de l’oued Boumoussa et la source Oum Karma). L’échantillonnage de novembre 2010 a consisté à prélever, sur le plateau de Settat, sept échantillons (deux dans des puits, deux dans l’oued Boumoussa, un dans l’oued Zokch et deux dans des sources (Hlilifa et Baida).

Les forages ont été purgés jusqu’à ce que les paramètres physicochimiques (pH et conductivité électrique) soient stabilisés. Les sources, l’oued Boumoussa et l’oued Zokch ont été échantillonnés par prélèvement direct dans le cours d’eau. Les échantillons mis dans des bouteilles en polyéthylène ont été conservés à 4 °C jusqu’au moment de leur analyse. Les paramètres physicochimiques (conductivité électrique, pH et température) ont été mesurés in situ. La profondeur de l’eau par rapport au sol a été mesurée pour l’ensemble des puits.

Les analyses chimiques (ions majeurs, ions azotés et éléments traces métalliques) des échantillons prélevés en 2009 ont été réalisées à la division de l’Unité d’Appui technique à la Recherche scientifique (UATRS) du Centre National de la Recherche Scientifique à Rabat. Les analyses chimiques des échantillons de 2010 ont été réalisées à l’Institut de l’Eau de l’Université de Grenade (Espagne) et ont concerné les ions majeurs et les nitrates. La composition isotopique de l'eau a été analysée en 2009 à l’Alfred-Wegener Institut for Polar and Marine Research, Isotope Laboratory (Research Unit Potsdam, Allemagne), et en 2010 au laboratoire IDES de l’Université Paris-Sud (France), dans les deux cas par équilibration puis mesures par spectrométrie gazeuse. Les teneurs en 2H et 18O sont exprimées selon la notation conventionnelle δ (‰) par rapport au V-SMOW. Les erreurs analytiques sont de ±1‰ et ±0,1‰ vs V-SMOW respectivement pour le δ2H et le δ18O.

3. Résultats et discussion

3.1 Géochimie et types d'eau

La pollution chimique des eaux de la nappe de l’aquifère de Settat par les écoulements superficiels de l’oued Boumoussa et par le recyclage des eaux souterraines pour l’irrigation a été décelée à partir des résultats de la campagne d’analyse physicochimique de 2004. Les eaux souterraines présentaient des niveaux de salinité élevés, plus prononcés dans les périmètres irrigués par les eaux souterraines. Les teneurs en nitrates de l'eau souterraine dépassaient largement la norme de 50 mg NO3∙L-1 préconisée par le Maroc (NORMES MAROCAINES, 2002). Les teneurs en nitrites étaient plus élevées au niveau de la zone d’épandage des eaux usées au nord de la ville de Settat et diminuent vers l’aval. Le plomb, le chrome, le zinc et l’aluminium présentaient des teneurs importantes dans les puits ainsi qu’au niveau des eaux usées, mettant en évidence l’apport des rejets industriels (EL ASLOUJ et al., 2007); ce fait a été observé sur l’oued Moulouya témoignant de l’impact des rejets urbains de la ville de Taourirt au Maroc (EL FOUNTI, 2000, cité par BOUABDLI et al., 2005). Ces valeurs en métaux non tolérées par les normes marocaines pour l’irrigation font que ces eaux soient inaptes à l’irrigation et peuvent constituer une menace pour la santé publique et l’environnement. Aussi, elles peuvent causer une contamination importante de l’horizon de surface du sol (DERE et al., 2006).

Les cartes de distributions spatiales des conductivités électriques de 2004 (Figure 3a) et 2009 (Figure 3b), correspondant à la période avant et après la mise en place de la STEP (janvier 2004 et juin 2009, respectivement), montrent que la salinité des eaux souterraines a augmenté au sud de la zone étudiée (Oulad Boukhalou). Ceci s’explique par la présence d’irrigation par les eaux acheminées par l’oued contaminé par les effluents industriels qui ne sont plus dilués par les rejets domestiques acheminés vers la STEP depuis 2007.

Figure 3

Cartes d’isoconductivité pour a) janvier 2004 et b) juin 2009 (EL AMRANI et al., 2011). Cartes en coordonnées Lambert; les lignes noires sont des courbes d’isoconductivité (µS∙cm-1)

Iso-conductivity maps for a) January 2004 and b) June 2009 (EL AMRANI et al., 2011). Maps in Lambert coordinates; black lines correspond to iso-conductivity levels (µS∙cm-1)

Cartes d’isoconductivité pour a) janvier 2004 et b) juin 2009 (EL AMRANI et al., 2011). Cartes en coordonnées Lambert; les lignes noires sont des courbes d’isoconductivité (µS∙cm-1)

-> Voir la liste des figures

Le diagramme de Piper élaboré à partir des analyses physicochimiques de 2009 et 2010 montre trois nuages de points (Figure 4, Tableau 1). Un premier groupement (I) à faciès sodique bicarbonaté sulfaté et chloruré est composé des eaux de l’oued Boumoussa (eaux usées [1, 4 et 5] et de la source [3]). Un deuxième groupement (II) à faciès sodique chloruré, est composé de puits captant le plioquaternaire (2, 6, 7, 8, 9, 10 et 12) où probablement il y a échange cationique calcium/sodium ou le calcium est retenu lors du passage des eaux au travers d’une formation argileuse (les argiles étant de bons échangeurs d’ions) (BANTON et BANGOY, 1997; TAUPIN et al., 2009; FEHDI et al., 2009). Le troisième groupement (III) à faciès calcique magnésien sodique et chloruré correspond aux puits situés les plus à l’aval sur la zone d’étude et les plus profonds. Comme ils présentent le même faciès que les puits (17 et 18) et les sources (16 et 20) situés sur le plateau de Settat, ces puits (11 et 13) captent probablement le cénomanien marno-calcaire et aussi ces eaux présentent un rapport Ca/SO4 supérieur à 1 ce qui montre qu’on a une dissolution des minéraux carbonatés (environnement calcaire). Le mélange entre les eaux de surface et les eaux souterraines est observé au niveau des diagrammes anionique et cationique, dans le sens de la flèche, par un enrichissement des eaux souterraines en chlorures et une diminution en calcium (puits 2, 9, 10, 6, 7, 8). De même une évaporation avant infiltration dans les zones irriguées peut induire un enrichissement en NaCl, ce phénomène a été observé par STRAUCH et al. (2006) dans les eaux superficielles et souterraines dans l’aquifère fracturé du bassin versant de l’Elqui dans la région des Andes au Chili.

Figure 4

Diagramme de Piper distinguant les types de faciès obtenus lors de la campagne de 2009- 2010. Les flèches indiquent le mélange entre les eaux de surface et les eaux souterraines. Les cercles indiquent les groupements de points à faciès hydrochimique distinct : I) faciès sodique bicarbonaté sulfaté, II) faciès sodique chloruré et III) facies calcique magnésien chloruré

Piper diagram distinguishing the types of facies obtained during the 2009-2010 campaign. Arrows indicate the mixing of surface water and groundwater. The circles indicate the groups of points with distinct hydrochemical facies (I, II and III)

Diagramme de Piper distinguant les types de faciès obtenus lors de la campagne de 2009- 2010. Les flèches indiquent le mélange entre les eaux de surface et les eaux souterraines. Les cercles indiquent les groupements de points à faciès hydrochimique distinct : I) faciès sodique bicarbonaté sulfaté, II) faciès sodique chloruré et III) facies calcique magnésien chloruré

-> Voir la liste des figures

Tableau 1

Résultats analytiques des campagnes d’échantillonnage de juin 2009 et novembre 2010

Analytical results of the June 2009 and November 2010 sampling campaigns

Résultats analytiques des campagnes d’échantillonnage de juin 2009 et novembre 2010

a OB : oued Boumoussa; OK: source Oum Karma; AH : source Aïn Hlilifa; AB : source Aïn Baida; OZ : oued Zokch

b Mètres au-dessus du niveau de la mer

-> Voir la liste des tableaux

De même, les puits à niveaux piézométriques élevés présentent les teneurs en nitrates les plus faibles, alors que les puits à niveaux piézométriques bas enregistrent des teneurs en nitrates les plus élevées (supérieures à 50 mg∙L-1) (Figure 5), ces derniers sont situés dans la zone irriguée par les eaux souterraines (Dladla). La présence des ions nitrates dans les niveaux aquifères profonds confirme l’infiltration des eaux depuis les niveaux aquifères superficiels vers des niveaux plus profonds (TRABELSI et al., 2005); aussi, un recyclage des eaux d’irrigation provenant des eaux souterraines (puits) se fait plus patent à ce niveau de la nappe. Les concentrations plus faibles dans les puits plus superficiels pourraient s'expliquer par des écoulements préférentiels rapides entre la zone non saturée et la zone saturée.

Figure 5

Relation entre les niveaux piézométriques et les concentrations en nitrates

Relationship between piezometric levels and nitrate concentrations

Relation entre les niveaux piézométriques et les concentrations en nitrates

-> Voir la liste des figures

3.2 Isotopes stables de l'eau

La composition de l'eau en isotopes stables permet de mieux comprendre le fonctionnement des aquifères, facilitant ainsi la gestion de la ressource en eau. Généralement, l'eau conserve sa signature isotopique sauf si elle est diluée ou mélangée avec des eaux d'une autre composition isotopique (FONTES 1980). Cependant, ce type d’investigation nécessite au préalable une bonne connaissance du signal isotopique véhiculé par la pluie (CELLE-JEANTONA et al., 2001). La composition isotopique des précipitations locales n’étant pas disponible, la droite des eaux météoriques mondiales est utilisée.

Les stations représentatives du plateau de Settat s'alignent sur une droite d’évaporation de pente 6,81 sous la droite météorique mondiale (GMWL) (Figure 6). Les stations de la plaine de Berrechid se répartissent sur une droite d’évaporation de pente 2,60, ce qui témoigne que ces eaux ont subi une évaporation avant d’atteindre les niveaux aquifères par infiltration (CLARK et FRITZ, 1997). Des résultats similaires ont été observés au niveau de l’aquifère de Takesla en Tunisie (TRABELSI et al., 2005), de Souss (BOUCHAOU et al., 2008) et de Tadla (BOUCHAOU et al., 2009) au Maroc. Il est intéressant de noter que l’enrichissement des eaux souterraines en δ18O est plus remarquable pour les eaux à salinité élevée, ce qui reflète l’effet de l’évaporation. Les puits ayant les niveaux piézométriques les plus élevés présentent également les signatures en δ18O les plus enrichies, ce qui reflète l'importance de l'évaporation à partir de la zone non saturée. Aussi, faut-il noter que ces eaux voient leur salinité s’accroitre à cause de l’effet anthropique de l’infiltration des eaux usées industrielles à salinité élevée déversées directement dans l’oued Boumoussa.

Figure 6

Composition en isotopes stables de l'eau des échantillons étudiés

Stable isotope composition in the water samples studied

Composition en isotopes stables de l'eau des échantillons étudiés

-> Voir la liste des figures

Le déplacement des points relatifs aux eaux souterraines par rapport à la GMWL offre une réelle estimation de la recharge par les précipitations. Cela a été possible grâce aux analyses isotopiques en appliquant la méthode d’ALLISON et al. (1984) basée sur celle de CRAIG (1961). La méthode d’ALLISON et al. (1984) consiste à estimer la recharge (R en mm∙a-1) à partir de la déviation des eaux souterraines par rapport à la droite météorique par l’une des formules empiriques soit par le δ18O = 3/√R ou le δ2H = 22/√R. En appliquant cette méthode, la recharge par les précipitations au niveau de la plaine de Berrechid a été estimée à 3,36 mm∙a-1. Dans notre cas, la recharge par les précipitations est faible, ce qui est observé aussi par l’éloignement des points par rapport à la GMWL (Figure 6), mais on remarque aussi que les points relatifs à l’oued Boumoussa (15, 19) au niveau du plateau de Settat sont les plus proches de la GWML, ce qui reflète l’apport important par les précipitations à ce niveau; cela a été confirmé par la valeur estimée de la recharge par les précipitations de 20 mm∙a-1 au niveau du plateau de Settat.

YIN et al. (2011) ont appliqué cette même méthode d’ALLISON et al. (1984) et la méthode du bilan de masse en chlorures pour évaluer la recharge au niveau de l’aquifère crétacé dans le plateau d’Ordos (Chine) et ils ont trouvé que la méthode d’ALLISON et al. (1984) donne une valeur de la recharge beaucoup plus faible comparativement avec la valeur évaluée par la méthode de bilan de masse en chlorures. Il est à signaler que la méthode d’ALLISON et al. (1984) est une formule empirique développée pour l’aquifère en sable du Bao’s groupe (Chine) dans des conditions hydroclimatiques différentes; nonobstant, elle donne une idée sur la dimension de la recharge.

Au niveau de la plaine de Berrechid, la recharge par les eaux d’irrigation est significativement plus importante; cette zone étant irriguée par les eaux souterraines et par les eaux de la rivière contenant les effluents liquides industriels, ce qui est montré par le déplacement des points par rapport à la GMWL qui ne peut être qu’un reflet du mélange des différents types de recharge à la nappe (LAROCQUE et PHARAND, 2010; ZOUARI et al., 2011). Ce processus de mélange est aussi observable et indiqué par les flèches sur les triangles anionique et cationique du diagramme de Piper de la figure 4 ou les eaux souterraines s’enrichissent en chlorures et en sodium et s’appauvrissent en calcium dans le sens de l’écoulement des eaux souterraines.

Comme il a été signalé, l’enrichissement des eaux souterraines en δ18O est plus remarquable pour les eaux à salinité élevée. Trois groupements de points se distinguent lorsque les conductivités électriques sont comparées aux compositions en δ18O (Figure 7a). Le premier groupement (G1) concerne les échantillons pris dans les puits profonds (Cénomanien) qui sont moins enrichis en δ18O, le deuxième groupement (G2) concerne les échantillons pris dans les puits superficiels (Plioquaternaire) plus enrichis en δ18O(évaporation et mélange) et le troisième groupement (G3) est relatif aux échantillons du Plioquaternaire, plus enrichis en δ18O et aux valeurs de conductivités électriques les plus élevées. Les puits où les niveaux piézométriques sont les plus élevés présentent les δ18O les plus importants (Figure 7b), ce qui reflète une évaporation à partir de la zone non saturée beaucoup plus importante et que l’évaporation est d’autant plus importante que les eaux sont prélevées près de la surface (ADORNI-BRACCISI et al., 2001), ce qui a été rencontré dans de nombreuses études par un enrichissement isotopique affectant les eaux souterraines quand le niveau piézométrique se rapproche de la surface (ALLISON, 1995; IAEA, 1983).

Figure 7

Relations entre les compositions en δ18O et a) les conductivités électriques et b) les niveaux piézométriques. G1 : échantillons pris dans les puits profonds (Cénomanien) qui sont moins enrichis en δ18O; G2 : échantillons pris dans les puits superficiels (Plioquaternaire) plus enrichis en δ18O (évaporation et mélange) et G3 : échantillons du Plioquaternaire, plus enrichis en δ18O et aux valeurs de conductivité électrique les plus élevées

Relationship between compositions in δ18O and a) electrical conductivity and b) piezometric levels. G1: samples taken in deep wells (Cenomanian) which are less enriched in δ18O; G2: samples taken from superficial wells (Plioquaternary) enriched in δ18O (evaporation and mixing) and G3: Plioquaternary samples, more enriched in δ18O and with the highest electrical conductivity

Relations entre les compositions en δ18O et a) les conductivités électriques et b) les niveaux piézométriques. G1 : échantillons pris dans les puits profonds (Cénomanien) qui sont moins enrichis en δ18O; G2 : échantillons pris dans les puits superficiels (Plioquaternaire) plus enrichis en δ18O (évaporation et mélange) et G3 : échantillons du Plioquaternaire, plus enrichis en δ18O et aux valeurs de conductivité électrique les plus élevées

-> Voir la liste des figures

3.3 Modèle conceptuel de l’écoulement

Un modèle conceptuel a été développé à partir de l’ensemble des résultats, et intègre les aspects hydrogéochimiques (Figure 8). Les résultats de l’étude physicochimique montrent que les eaux souterraines présentent des conductivités électriques élevées et que celles-ci sont plus élevées dans les périmètres irrigués par les eaux de la rivière Boumoussa, ce qui indique l’infiltration vers l’aquifère d’une eau de surface formée principalement de rejets industriels très chargés en contaminants. La présence de nitrates dans les puits à niveau piézométrique bas confirme la possibilité d’une recharge de l'aquifère régional par une eau ayant percolé à travers le sol agricole fertilisé. Aussi la lithologie des formations traversées contribue à donner un faciès particulier aux eaux souterraines. Trois types de faciès sont distingués :

  • groupe (I) à faciès sodique bicarbonaté sulfaté et chloruré, est composé des eaux de l’oued Boumoussa et de la source Oum Karma au niveau de la plaine;

  • groupe (II) à faciès sodique chloruré, est composé de puits captant le Plioquaternaire;

  • groupe (III) à faciès calcique magnésien sodique et chloruré correspond aux puits situés les plus à l’aval sur la plaine et les plus profonds et les puits et sources situés sur le plateau de Settat qui captent le Cénomanien marno-calcaire.

Figure 8

Coupe schématique du modèle conceptuel de circulation hydrogéologique de la zone d’étude

Schematic cross-section representing the conceptual hydrogeological model of the study area

Coupe schématique du modèle conceptuel de circulation hydrogéologique de la zone d’étude

-> Voir la liste des figures

De même, le modèle proposé peut aider à une meilleure gestion de la ressource en eau dans cette région par une revalorisation des eaux épurées de la STEP dans les besoins de l’industrie et de l’agriculture et par une utilisation des eaux souterraines dans les besoins en approvisionnement en eau potable.

4. Conclusion

La croissance économique et démographique de la région de Settat a été accompagnée par une augmentation des débits des effluents liquides urbains et industriels qui sont directement évacués dans l’oued Boumoussa et qui ont induit une détérioration de la qualité de ces eaux superficielles et souterraines. L’objectif de ce travail était de réaliser une synthèse des données physicochimiques et isotopiques afin de proposer un modèle conceptuel de circulation des eaux au niveau de l’aquifère plioquaternaire superficiel de Settat au Maroc, notamment pour déceler l’origine des apports à l’aquifère provenant de la surface et de l’aquifère sous-jacent.

Les résultats confirment le comportement hydrodynamique et hydrogéochimique des eaux souterraines. L’enrichissement des eaux souterraines en δ18O est plus remarquable pour les eaux à salinité élevée, ce qui reflète l’effet de l’évaporation. Les puits ayant les niveaux piézométriques les plus élevés présentent également les signatures en δ18O les plus enrichies, ce qui reflète l'importance de l'évaporation à partir de la zone non saturée. Aussi faut-il noter que ces eaux voient aussi leur salinité s’accroitre à cause de l’effet anthropique de l’infiltration des eaux usées industrielles à salinité élevée déversées directement dans l’oued Boumoussa. Au niveau de la nappe de Settat, la recharge naturelle par les précipitations est faible, ce qui se traduit par l'alignement des compositions isotopiques de l'eau souterraine sur une droite de pente inférieure à celle de la droite météorique mondiale. Le déplacement des points par rapport à celle-ci peut aussi être un reflet du mélange de la recharge naturelle avec la recharge par les eaux d'irrigation (provenant de la rivière Boumoussa et des eaux souterraines surtout à l’aval de la zone d’étude au niveau de la localité Dladla) et par les eaux provenant de l'aquifère cénomanien du plateau de Settat.

Le modèle proposé peut aider à une meilleure gestion de la ressource en eau par une revalorisation des eaux épurées de la STEP dans les besoins de l’industrie et de l’agriculture et par une utilisation des eaux souterraines provenant du Cenomanien, vu leur meilleure qualité, dans les besoins en eau potable.