Analyse multicritère couplée au radiocarbone en contexte de variabilité climatique pour l’identification des zones potentielles de recharge des aquifères fracturés du Sud-Est de la région de la Nawa (Côte d’Ivoire)

Présentation de la zone d’étude

Située entre les longitudes 6°13’ et 6°39’ Ouest et les latitudes 5°29’ et 5°50’ Nord, la zone d’étude (Figure 1) est une portion du Sud-Est de la région de la Nawa, région située au Sud-Ouest de la Côte d’Ivoire. Cette région est limitée au Nord par les régions du Guémon et du Haut-Sassandra, au Sud par les régions du Gbôklè et de San-pédro, à l’Ouest par la région du Cavally et à l’Est par les régions du Gôh et du Lôh Djiboua. La zone d’étude s’étend sur une superficie de 9193 km².

La Nawa est composée de quatre départements: Soubré (chef-lieu), Méagui, Guéyo, Buyo. Elle compte également six communes et dix sous-préfectures. La population de la région est estimée à 1 053 084 habitants dont 464 554 habitants pour le département de Soubré (RGPH, 2014RGPH (Recensement Général de la Population et de l’Habitat) (2014). Répertoire des localités. Région de la NAWA, 29 pp.). Cette population est composée d’autochtones (Bété, Bakwé, Kouzié et Godié), d’allochtones (Baoulé et Malinké) et d’allogènes (Burkinabè, Maliens et libanais). Plus de 70 % de la population vit en milieu rural. Les aptitudes culturales des sols de la région sont propices tant aux cultures de rente (café, cacao, palmier à huile et hévéa), qu’à une gamme variée de cultures vivrières (riz, manioc, maïs, banane plantain, igname) et maraîchères généralement exploitées dans les bas-fonds (PEMED-CI, 2015PEMED-CI (2015). Etudes monographiques et économiques des Districts de Côte d’Ivoire, District du Bas-Sassandra. Note de synthèse, 69 pp.). La Nawa est la principale région productrice de cacao de la Côte d’Ivoire fournissant plus de 40% (500 000 T) de la production nationale. L’agriculture est la principale source de revenu représentant 75% des activités économiques, et contribue à 94% des revenus (Tahoux, 1993 in Kouamé et al., 2016Kouamé, C.; Diomande S.; Kouassi, J.L.; Yeboi, B. & Diby, L. (2016). Dynamiques d’utilisation des terres dans la Région de la Nawa. Atelier sur l’étude analytique de l’occupation du sol en vue de l’intégration du mécanisme REDD+ dans le Schéma Régional d’Aménagement et de Développement du Territoire (SRADT) de la Région de la Nawa, Soubre, Côte d’Ivoire. http://doi.org/10.13140/RG.2.1.3766.5523 ). 93 % et 52,46% des paysans vivent respectivement du cacao et du café contre 31,2% pour les cultures vivrières (Kouadjo et al., 2002 in Kouamé et al., 2016Kouamé, C.; Diomande S.; Kouassi, J.L.; Yeboi, B. & Diby, L. (2016). Dynamiques d’utilisation des terres dans la Région de la Nawa. Atelier sur l’étude analytique de l’occupation du sol en vue de l’intégration du mécanisme REDD+ dans le Schéma Régional d’Aménagement et de Développement du Territoire (SRADT) de la Région de la Nawa, Soubre, Côte d’Ivoire. http://doi.org/10.13140/RG.2.1.3766.5523 ).

La géologie de la zone d’étude (Figure 2), qui s’intègre à celle du domaine SASCA (Sassandra-Cavally), est constituée de formations archéennes (migmatites et gneiss de Grabo) et éburnéennes (granitoïdes éburnéens concordants et syénites du mont Troquoi) (Papon & Lemarchand, 1973Papon, A. & Lemarchand, R. (1973). Géologie et minéralisation du sud-ouest de la Côte d’Ivoire. Synthèse des travaux de l’opération SASCA (19621968). Mémoire B.R.G.M., 80: 284 pp.).

Les migmatites font partie de l’ensemble granito-gneissique cartographié par Bolgarsky (1950)Bolgarsky, M. (1950). Étude géologique et description pétrographique du Sud-Ouest de la Côte d’Ivoire. Bulletin de la Direction des mines, 9, Dakar 55 pp..

Elles forment une auréole autour de la granodiorite d’anatexie et s’étend du Sud-Ouest (Kragui) jusqu’au-delà de la limite Nord-Est du département de Soubré. Les roches appartenant à ces formations sont des gneiss métamorphisés dans la catazone (niveau structural des gneiss inférieurs ou base du faciès amphibolite) et affectés par une anatexie différentielle.

Les gneiss de Grabo sont intrudés par les massifs de granites à deux micas éburnéens rencontrés au Sud-Est (zone d’Okrouyo). Ses affleurements se présentent sous forme de bancs avec une alternance de lits phylliteux et de lits quartzeux.

Les granitoïdes éburnéens concordants cartographiés sous le nom de granodiorite d’anatexie de type « rivière Bao », affleurent dans cette rivière non loin de la voie joignant Soubré à Okrouyo. Ce massif est entièrement pris dans les migmatites. Il se distingue de son encaissant par sa grande richesse en feldspath potassique et par sa texture granoblastique encore plus engrenée.

Les syénites du mont Troiquoi forment six mas sifs d’inégale importance, encaissés dans les gneiss à biotite ou dans les migmatites. Trois de ces massifs correspondent aux reliefs regroupés sous le nom de Mont Troiquoi. L’un d’eux affleure au Nord de Koziyao et à l’Est d’Okrouyo.

Du point de vue géomorphologique, la zone d’étude fait partie d’une vaste pénéplaine légèrement ondulée, dont l’altitude varie de 70 m à plus de 270 m. Les vallées (70 m à 170 m), occupent la majeure partie de la zone d’étude. Les bas plateaux (170 m-250 m) et les hauts plateaux (250 m-270 m) sont surmontés par endroit de quelques élévations comme le mont Troiquoi (444 m) à Okrouyo (PEMED-CI, 2015PEMED-CI (2015). Etudes monographiques et économiques des Districts de Côte d’Ivoire, District du Bas-Sassandra. Note de synthèse, 69 pp.). D’autres élévations existent dans la région de la Nawa comme celle du mont Zatro à Grand Zattry développés sur des roches éruptives anciennes (Yao, 2009Yao, K.T. (2009). Hydrodynamisme de l’eau souterraine dans les aquifères de socle cristallin et cristallophyllien du Sud-Ouest de la Côte d’Ivoire: cas du département de Soubré. Apports de la télédétection, de la géomorphologie et de l’hydrogéochimie. Thèse de Doctorat, Conservatoire National des Arts et Métiers, Paris, France, 284 pp.).

La région est drainée par le fleuve Sassandra et ses nombreux affluents qui participent à la recharge de la nappe d’eau souterraine. La puissance du fleuve a engendré la construction du barrage de Soubré en 2017, qui produit de l’électricité en soutien à la capacité énergétique déjà disponible en Côte d’Ivoire. Il permet d’irriguer une grande partie de la Nawa, au bénéfice de l’agriculture et apporte d’importantes ressources halieutiques.

Du point de vue climatique, la zone d’étude est soumise à un climat subéquatorial ou climat attiéen (Eldin, 1971Eldin, M. (1971). Le climat. Mémoire ORSTOM, 50: 77108.). La température a une valeur moyenne d’environ 25°C avec une amplitude annuelle de l’ordre de 3°C. La végétation de la Nawa, initialement marquée par la forêt sempervirente intermédiaire, dense et humide, s’est réduite aujourd’hui au profit des lambeaux de forêts et d’immenses plantations de cultures pérennes traditionnelles ou industrielles.

Sur le plan pédologique, la quasi-totalité de la zone d’étude est constituée de sols ferrallitiques fortement désaturés, à l’exception d’une infime portion qui comporte des sols ferrallitiques moyennement et/ou faiblement désaturés, notamment au Sud-Est de Soubré et au Sud-Ouest de Ottawa.

Données et matériel

Plusieurs données ont été utilisées dans le cadre de cette étude. Il s’agit de:

Les différentes données ont fait l’objet de traitement à l’aide d’un certain nombre de logiciels. Le logiciel Khronostat 1.01, mis au point par l’Institut de Recherche et de Développement (IRD) (ex-ORSTOM), a permis de traiter les données pluviométriques en vue de caractériser la variabilité climatique. Le tableur Excel™ de Microsoft Office 2010 a facilité la réalisation les différents tableaux et graphiques. Le logiciel Envi 5.3 a été utilisé pour traiter les images satellitaires en vue de l’extraction des différentes structures linéamentaires. Les linéaments structuraux ont été traités avec le logiciel Linwin 2.0 pour la réalisation des rosaces directionnelles. Le logiciel Mapinfo 7.5 a permis de numériser les différentes cartes. ArcGis 10.5 a facilité le croisement de cartes thématiques pour aboutir à une carte de synthèse.

Caractérisation de la variabilité climatique et détermination des paramètres hydrologiques

Les changements dans le régime du climat de la zone d’étude sont relevés par des dates charnières (UICN-BRAO et al., 2003UICN-BRAO, GWP-WAWP & CILSS (2003). Eau, changement climatique et désertification en Afrique de l’Ouest: Stratégie régionale de préparation et d’adaptation. Ouagadougou, Burkina-Faso, 57 pp.) qui marquent les ruptures à l’intérieur de la série climatique. Le test de Pettitt et la segmentation de Hubert, réputés pour leur robustesse (Lubès-Niel et al., 1998Lubès-Niel, H.; Masson, J.M.; Paturel, J.E. & Servat, E. (1998). Variabilités climatiques et statistiques. Etude par simulation de la robustesse de quelques tests utilisés pour vérifier l’homogénéité de chroniques. Revue des Sciences de l’Eau, 3: 383-408. https://doi.org/10.7202/705313ar ; Kingumbi et al., 2000Kingumbi, A.; Bergaoui, Z.; Bourges, J.; Hubert, P. & Kallel, R. (2000). Etude de l’évolution pluviométrique de la Tunisie Centrale. Documents Techniques en Hydrologie, 51: 341-345. ; Bigot et al., 2002Bigot S.; Brou, T.Y.; Bonnardot, V. & Servat, E. (2002). Stabilité inter-annuelle des précipitations en Côte d’Ivoire de 1950 à 1996. Hydrological Sciences Journal, 274: 507-514.; Ardoin-Bardin, 2004Ardoin-Bardin, S. (2004). Variabilité hydroclimatique et impacts sur les ressources en eau de grands bassins hydrographiques en zone soudano-sahélienne. Thèse de Doctorat, Université de Montpellier II, France, 437 pp.; Berthelot et al., 2004Berthelot, M.; Genin, A. & Gourdon, D. (2004). Variabilité interannuelle des précipitations dans les Alpes du Sud pour la période 1951/2000: Analyse statistique des séries. Université François Rabelais, Tours, France, 15 pp. http://alexandre.letort.pagesperso-orange.fr/berthelot/alpessud.htm (consultation le 22 novembre 2021).) sont utilisés pour la mise en évidence de ruptures au sein de la série pluviométrique annuelle de 1951 à 2017 de Soubré. Le logiciel Khronostat 1.01 (Boyer et al., 1998Boyer, J.F.; Lubès-Niel, H.; Paturel, J.E. & Servat E. (1998). Khronostats-Logiciel de détection des fluctuations dans une série chronologique par application de tests statistiques. IRD, UMR Hydrosciences.) a servi à la mise en œuvre de ces méthodes de détection de ruptures.

Le déficit pluviométrique de part et d’autre de l’année de rupture, permet de comprendre le comportement de l’évolution de la pluviométrie entre les périodes antérieure et postérieure à la rupture. Il est évalué à partir de la formule 1:

Où D représente le déficit pluviométrique, ${\stackrel{-}{x}}_J$ la moyenne de la pluviométrie après la rupture et, ${\stackrel{-}{x}}_i$ la moyenne de la pluviométrie avant la rupture.

La fluctuation interannuelle de la pluviométrie permettant d’identifier les grandes tendances climatiques (période sèche, période normale et période humide) du département de Soubré, est déterminée à partir de l’indice pluviométrique qui correspond à une « variable centrée et réduite ». Cet indice est calculé à partir de la formule de Nicholson (1980)Nicholson S.E. (1980). The nature of rainfall fluctuations in subtropical West Africa. Monthly Weather Review, 108: 473-487. équation (2):

Où Ip représente l’indice pluviométrique, x_i la hauteur de pluie totale pour une station donnée pendant une année i, $\stackrel{-}{x}$ la moyenne interannuelle de la pluie à la station pendant la durée entière de l’enregistrement (période d’étude), et $\sigma$ l’écart-type de la pluviométrie annuelle.

L’estimation de la quantité d’eau infiltrée dans le sol permettant d’obtenir le taux de recharge des aquifères de la zone d’étude passe par la détermination de son bilan hydrologique. Les paramètres hydrologiques (Evapotranspiration Potentielle, Evapotranspiration réelle, Excédent, Déficit, etc.) sont déterminés en suivant la méthode de Thornthwaite (1944)Thornthwaite, C.W. (1944). The measurement of potential evapotranspiration. J. M. Mather, Sea-Brook, New Jersey, 225 pp..

La lame d’eau infiltrée pour recharger les aquifères de 1951 à 2017 est appréciée à partir des fluctuations de l’infiltration (I) (équation 4). Celle-ci est déduite à partir de l’équation 3 du bilan hydrologique (équation 3). Elle est calculée en soustrayant de la pluviométrie (P), la lame ruisselée (R) et évapotranspirée (ETR).

Où P est la pluie, ETR l’évapotranspiration réelle, S l’excédent ou le surplus, R le ruissellement et I l’Infiltration.

L’approche de la représentation spatiale de la recharge consiste à combiner les informations spatiales relatives à différents critères susceptibles d’influencer l’infiltration.

Modélisation spatiale des zones de recharge de l’aquifère par analyse multicritère et technique de validation

Chaque donnée de l’étude, soumise à un traitement approprié (saisie manuelle, tri, interpolation, numérisation, digitalisation, maillage, traitement numérique), aboutit à l’obtention d’une carte thématique qui représente un critère de décision intervenant dans la réalisation de la recharge. Les critères disponibles intervenant dans la recharge des aquifères fracturés de la zone d’étude, sont au nombre de 6 en référence aux travaux de Assoma et al. (2012)Assoma, T.V.; Adiaffi, B. & Koudou, A. (2012). Contribution of Remote Sensing, Multi-Criteria Analysis and GIS for Mapping Recharge Areas of the Coastal Aquifers in Southeast Ivory Coast. American Journal of Scientific Research, 80: 112-125., Oularé et al. (2017)Oularé, S.; Adon, G.C.; Akpa, Y.L.; Saley, M.B.; Kouamé, K.F. & Therrien, R. (2017). Identification Des Zones Potentielles De Recharge Des Aquifères Fracturés Du Bassin Versant Du N’zo (Ouest De La Côte d’Ivoire): Contribution Du SIG Et De La Télédétection. European Scientific Journal, 13: 192-217. https://doi.org/10.19044/esj.2017.v13n36p192 , Aké et al. (2018)Aké, G.E.; Kouamé, K.J.; Koffi, A.B. & Jourda, J.P. (2018). Cartographie des zones potentielles de recharge de la nappe de Bonoua (sud-est de la Côte d’Ivoire). Revue des Sciences de l’Eau, 31: 129-144. https://doi.org/10.7202/1051696ar et Koudou et al. (2021)Koudou, A.; Adiaffi, B. & Kra, C.C. (2021). Cartographie des zones de recharge des aquifères de fractures du département de Tanda par analyse multicritère (Nord-Est de la Côte d’Ivoire). Bulletin de l’Institut Scientifique, Rabat, Section Sciences de la Terre, 43: 27-42..

La lame d’eau infiltrée, qui correspond au gain d’eau à la recharge des aquifères fracturées, est évaluée à partir des données pluviométriques par le calcul du bilan hydrique.

La lithologie se rapporte à la nature des roches issues des formations géologiques de la zone d’étude. L’aptitude à la recharge varie d’une formation à une autre.

L’occupation du sol décrit la couverture biophy sique de la surface des terres émergées; c’est-à-dire ce qui recouvre le sol (FAO, 1998FAO (1998). Les conséquences de l’accord sur l’agriculture du cycle d’Uruguay pour les pays en développement.). Les différentes catégories biophysiques, les zones de végétation (arbres, buissons, champs, pelouses), les sols nus (même s’il s’agit d’un manque de couverture), les surfaces dures (roches, bâtiments), les surfaces humides et les plans d’eau (nappes et cours d’eau, zones inondables), présentes à la surface du sol, influencent différemment (favorablement ou défavorablement) la recharge des aquifères. L’occupation du sol de la zone d’étude est élaborée à partir du logiciel Envi 5.3. Les images OLI de Landsat 8 ont fait l’objet d’un prétraitement (extraction de la zone d’étude) suivi d’un traitement numérique proprement dit (compositions colorées, ACP, etc.). Plusieurs images transformées de la zone d’étude permettent une bonne discrimination spectrale des différents types d’occupation du sol en général, et de la couverture végétale en particulier (Girard & Girard, 1999Girard M.-C. & Girard, C.M. (1999). Traitement des données de télédétection. Interprétation physique des données, les comportements spectraux, 529 pp.). Une classification, réalisée conjointement à l’interprétation visuelle de la composition colorée des bandes 8, 6 et 5 de qualité optimale, se base sur la connaissance de l’occupation des sols acquise sur la carte de Yao (2009)Yao, K.T. (2009). Hydrodynamisme de l’eau souterraine dans les aquifères de socle cristallin et cristallophyllien du Sud-Ouest de la Côte d’Ivoire: cas du département de Soubré. Apports de la télédétection, de la géomorphologie et de l’hydrogéochimie. Thèse de Doctorat, Conservatoire National des Arts et Métiers, Paris, France, 284 pp.. La classification, supervisée par le choix des parcelles d’entraînement, amène à la définition de 5 classes radiométriquement différentes. Cette classification est évaluée par différents tests de performances (la matrice de confusion, la précision globale et le coefficient kappa). Un filtre médian 3×3 est appliqué à la carte d’occupation du sol pour réduire l’hétérogénéité intra-classes.

La densité de fracturation ou l’intensité de la fracturation s’exprime soit en nombre de fractures, soit en longueurs cumulées de fractures par unité de surface. Sa connaissance dans un milieu traduit le degré de fracturation des roches. Les zones hautement fracturées sont les endroits perméables à travers lesquels l’eau s’infiltre. Les différentes techniques de traitement des images OLI de Landsat 8 à partir de Envi 5.3, notamment l’Analyse en composantes principales, les rapports de bandes, les compositions colorées, le filtrage directionnel de Sobel (7×7) ont facilité l’interprétation visuelle et le relevé des linéaments structuraux. Ceux-ci ont été validés à partir des débits des forages d’eau productifs, du réseau hydrographique de la zone d’étude. La carte qui en résulte est soumise à un maillage de taille régulière de 5×5 km² à partir de Linwin 2.0. Les longueurs cumulées de linéaments structuraux par maille sont interpolées selon la technique IDW pour réaliser la carte thématique de densité de fracturation.

La densité de drainage résulte du réseau hydrographique de la zone d’étude qui est numérisé à partir de la carte topographique. Elle est réalisée à partir de l’interpolation selon IDW des longueurs cumulées des drains par maille.

La pente est destinée à faciliter l’écoulement des eaux. Les pentes faibles sont favorables à l’infiltration donc à la recharge des aquifères, alors que les pentes fortes favorisent le ruissellement, et donc constituent un caractère défavorable à la recharge d’un aquifère. La réalisation de la carte thématique est faite automatiquement sous ArcGis 10.5 à partir du Modèle Numérique d’Elévation (MNE). Celui-ci est obtenu à partir des données SRTM (Shuttle Radar Topography Mission) avec une résolution de 30 m.

Les différents critères de décision sont structurés au sein d’une base de données géospatiales. Les plans d’informations qui en résultent sont mis en conformité géométrique dans le système de référence cartographique longitude / latitude WGS 1984.

Chaque critère considéré est subdivisé en classes qui représente soit un milieu particulier, soit un intervalle de valeurs. Les différentes classes sont standardisés pour une bonne analyse multicritère. A cet effet, ces classes sont attribuées d’une note variant de 1 à 10 en fonction de leur influence particulière (défavorable ou favorable) sur la recharge des aquifères. Chaque critère considéré est aussi pondéré sur la base de son importance relative sur la recharge des aquifères. La procédure mise en œuvre pour déterminer les poids à attribuer à chaque critère se décompose en quatre (4) étapes:

Les poids attribués à chaque critère varient entre 0 et 1 et leur somme est égale à 1. Le développement théorique de la détermination des poids de chaque critère peut être consulté auprès de Koudou (2013)Koudou, A.; Adiaffi, B.; Assoma, T.V.; Sombo, A.P.; Amani, E.M.E. & Biémi, J. (2013). Conception d’un outil d’aide à la décision pour la prospection des eaux souterraines en zone de socle du Sud-Est de la Côte d’Ivoire. Revue internationale de géologie, de géographie et d’écologie tropicales, 37: 211-226.. Le tableau 3 présente les classes des critères influençant la recharge qui sont classifiés, standardisés et pondérés.

Les différents critères classés, standardisés et pondérés sont aussi agrégés.

La méthode d’agrégation complète émet l’hypothèse que tous les critères considérés sont parfaitement comparables, et le résultat de la fonction objective qu’elle conçoit et qui intègre tous les critères considérés, peut être cartographié. Ainsi, l’agrégation complète des critères par pondération (Joerin, 1995Joerin, F. (1995). Méthode multicritère d’aide à la décision et SIG pour la recherche d’un site. Revue Internationale de Géomatique, 1: 37-51.; Martin et al., 2004Martin, J.L.; Henry, E.; Boulemia, C. & Masson, F.X. (2004). Proposition d’outils de gestion et d’aide à la décision dans une moyenne collectivité locale: application à l’exploitation de la ressource en eau souterraine. Actes de la journée d’études « Les territoires de l’eau », Université d’Artois, Arras: 74-82.) est appliquée dans la zone d’étude et consiste à la sommation des valeurs standardisées et pondérées de chaque critère intervenant dans la réalisation de la recharge. Elle est résumée par l’équation 4:

Où S est le résultat, W_i est le poids du critère i, et, X_i est la valeur standardisée du critère i.

Cette procédure aboutit à l’obtention d’un indice d’aptitude situé sur une échelle de 0 à 10. Cet indice est reporté dans l’espace et reclassifié en quatre classes qualifiées de Faible, Moyenne, Forte et Très forte, pour une meilleure lisibilité et une bonne interprétation de la carte de synthèse.

La validation de la carte des zones de recharge est effectuée à partir de la mesure de l’activité en carbone 14. L’activité en carbone 14 est un outil très efficace pour mettre en évidence les zones de recharge. Lorsque l’activité en carbone 14 est élevée (autour de 100% pMC), la recharge de la nappe au point d’échantillonnage est récente. En revanche, lorsqu’elle est basse (bien loin de 100% pMC dans le sens de la diminution), la recharge est dite ancienne (Oga, 1998Oga, Y.M.S. (1998). Ressources en eau souterraines dans la region du grand Abidjan (Côte d’Ivoire): Approches hydrochimique et isotopique. Thèse de Doctorat en Sciences, Université de Paris-Sud, Orsay, France, 290 pp.; Adiaffi, 2008Adiaffi, B. (2008). Apport de la géochimie isotopique, de l’hydrochimie et de la télédétection à la connaissance des aquifères de la zone de contact ‘’socle-bassin sédimentaire” du Sud-est de la Côte d’Ivoire. Thèse de Doctorat en Sciences, Université de Paris-Sud, Orsay, 230 pp.; Adiaffi et al., 2009Adiaffi, B.; Marlin, C.; Oga, Y.M.S.; Massault, M.; Nauret, A. & Biémi, J. (2009). Palaeoclimatic and deforestation effect on the coastal fresh groundwater resources of SE Ivory Coast from isotopic and chemical evidence. Journal of Hydrology, 369: 130-141. https://doi.org/10.1016/j.jhydrol.2009.02.046 ; Adiaffi et al., 2021Adiaffi, B.; Marlin, C.; Koudou, A.; Assoma, V.T.; Guédé, K.; Oga, Y.M.S.; Massault, M.; Noret, A. & Monvoisin, G. (2021). Geochemical processes in the groundwater of Soubré, Ivory Coast, West Africa: evidence from ¹⁴C, ¹³C, ¹⁸O, ²H and major ions. Journal of African Earth Sciences, 184: 1-14. https://doi.org/10.1016/j.jafrearsci.2021.104382 ). Les âges des eaux sont déterminés à partir des activités en carbone 14.

C’est une méthode fréquemment utilisée en géochimie des eaux pour apprécier le taux de renouvellement des aquifères étudiés (Girard, 1993Girard, P. (1993). Techniques isotopiques appliquées à l’étude des nappes des altérites et du socle fracturés de l’ouest africain. Etude de cas: l’ouest du Niger. Thèse de Doctorat en ressources minérales, Université du Québec à Chicoutimi, 138 pp. https://doi.org/10.1522/1486981 ; Oga, 1998Oga, Y.M.S. (1998). Ressources en eau souterraines dans la region du grand Abidjan (Côte d’Ivoire): Approches hydrochimique et isotopique. Thèse de Doctorat en Sciences, Université de Paris-Sud, Orsay, France, 290 pp.; Adiaffi, 2008Adiaffi, B. (2008). Apport de la géochimie isotopique, de l’hydrochimie et de la télédétection à la connaissance des aquifères de la zone de contact ‘’socle-bassin sédimentaire” du Sud-est de la Côte d’Ivoire. Thèse de Doctorat en Sciences, Université de Paris-Sud, Orsay, 230 pp.). Dans cette étude, le mode de validation consiste à superposer l’ensemble des forages où les données des A¹⁴C sont connues à la carte des zones de recharge, puis à établir une relation entre les classes de recharge et les A¹⁴C des eaux de forages.

La figure 3 illustre la méthodologie élaborée pour la cartographie des zones potentielles de recharge par analyse multicritère.

Résultats

Caractérisation de la variabilité climatique et estimation des paramètres du bilan hydrologique

 ⌅

Deux ruptures successives (Figure 4 et Tableau 4), l’une primaire et très significative au seuil de confiance de 99%, a été identifiée en 1970 et l’autre secondaire en 1992, s’observent dans la série pluviométrique annuelle (1951-2017) de Soubré.

Figure 4.  Courbe de la série pluviométrique de la station de Soubré aprés la variable U du test de Pettitt (dates de rupture: 1970 et 1992).

Tableau 4.  Test de segmentation d’Hubert

Début	Fin	Moyenne (mm)	Ecart-type (mm)
1951	1970	1635,605	254,489
1971	1992	1263,514	191,711
1993	2017	1520,256	214,228

Le déficit pluviométrique calculé de part et d’autre de la rupture primaire de 1970 est de 14%, confirmant ainsi une baisse significative de la moyenne des pluies au cours du XX^ème siècle à la station de Soubré (Tableau 5).

Tableau 5.  Déficit pluviométrique à la station de Soubré (1951-2017)

		Avant rupture		Après rupture
Station	Référence	Moyenne (mm)	Ecart-type (mm)	Moyenne (mm)	Ecart-type (mm)	Déficit (%)
Soubré	1970	1635,6	254,5	1400,1	239,8	14

L’indice pluviométrique de Nicholson (Figure 5) permet de confirmer les grandes tendances dans la série chronologique identifiée par les tests de rupture.

Figure 5.  Evolution interannuelle de la pluviométrie de Soubré de 1951 à 2017.

Soubré est caractérisé par 3 périodes:

• une période humide entre 1951 et 1970 avec une moyenne de 1635,6 mm et un écart-type de 254,5 mm;

• une période sèche de 1971 à 1992, avec une moyenne de 1263,5 mm et un écart-type de 191,7 mm;

• une période humide de 1993 à 2017 avec une moyenne de 1520,3 mm et un écart-type de 214,2 mm.

La variation mensuelle et saisonnière des précipitations avant et après la rupture de 1970 met en évidence 4 saisons à la station de Soubré (Figure 6). Une grande saison sèche de novembre à février pour la période avant la rupture, et qui se prolonge jusqu’à mars pour la période après rupture. Une petite saison sèche qui se résume au mois d’août avant la rupture, et qui s’étend sur 2 mois (juillet et août) après la rupture. Une grande saison pluvieuse qui s’étend de mars à juillet (5 mois) avant la rupture, et qui passe d’avril à juin, soit 3 mois après la rupture. La petite saison pluvieuse se résume à 2 mois (septembre et octobre) avant et après rupture. La pluviométrie de chaque mois avant la rupture reste sensiblement supérieure à celle après rupture, hormis les mois de février, mai et août durant lesquels la pluviométrie est presqu’identique.

Figure 6.  Variation mensuelle et saisonnière des précipitations avant et après la rupture de 1970 à Soubré.

Le bilan hydrologique à la station de Soubré sur la période 1951-2017 est consigné dans le Tableau 6.

Tableau 6.  Bilan hydrologique de la station de Soubré (1951-2017) selon la méthode de Thornthwaite

	J	F	M	A	M	J	J	A	S	O	N	D	Total
P(mm)	25	47	102	159	297	515	172	32	71	177	153	78	1828
ETP(mm)	140	144	166	160	155	130	114	102	110	132	141	143	1638
ETR(mm)	25	47	102	159	155	130	114	102	101	132	141	135	1446
RFU(mm)	0	0	0	0	100	100	100	30	0	45	57	0	432
DRFU(mm)	0	0	0	0	100	0	0	-70	-30	45	12	-57	00
P-ETR(mm)	0	0	0	0	142	385	58	-70	-30	45	12	57	485
Surplus(exc)°	0	0	0	0	42	385	58	0	0	0	0	0	485
Ecoulement	2	1	1	1	21	203	131	65	32	12	8	4	485

Pour une pluviométrie moyenne annuelle (19512017) de 1828 mm, la lame d’eau évapotranspirée est de 1638 mm. La lame d’eau ruisselée est de 242 mm. Celle destinée à recharger les aquifères du Sud-Est de la Nawa est évaluée à I = (485-242) = 243 mm.

Cartes thématiques du Sud-Est de la région de la Nawa

 ⌅

Cinq cartes thématiques ont été élaborées dans le Sud-Est de la région de la Nawa. Il s’agit de la carte d’occupation du sol, de la carte de densité de fracturation, de la carte de densité de drainage, de la carte des pentes, et de celle des zones de recharge validées du Sud-Est de la région de la Nawa et leur proportion.

La carte des zones de recharge validées du Sud-Est de la région de la Nawa et leur proportion, est illustrée à la Figure 7.

Figure 7.  Carte des zones de recharge validées du Sud-Est de la région de la Nawa et leur proportion.

Cette carte de synthèse est caractérisée par une variabilité spatiale des zones à potentialité de recharge des aquifères. Quatre (4) classes de zones de recharge, de proportion variable, se rencontrent dans la zone d’étude.

Les zones de recharge très forte et forte se localisent au Sud-Est et s’étendent jusqu’au Sud-Ouest de la zone d’étude. Ces zones représentent respectivement 40% et 35% de la superficie globale, et caractérisent les zones où le taux d’infiltration pourrait être élevé en cas de pluies. Elles désigneraient également les zones de transmission rapide des contaminants de la surface vers l’aquifère, donc elles doivent être protégées. À l’opposé, les zones de recharge faible sont localisées majoritairement à l’Ouest et au Nord de la zone d’étude. Quelques portions disséminées apparaissent au Sud et à l’Est. Ces endroits représentent environ 14% de la superficie totale de la zone d’étude. Enfin, les zones de recharge moyenne, qui occupent 11% de la superficie totale, apparaissent dans le lit du fleuve Sassandra dans le Sud-Est de la Nawa. Les fortes teneurs en carbone 14 se superposent aux zones de recharge forte et très forte.

La distribution spatiale des teneurs en carbone 14 est tout à fait en concordance avec les zones de recharge élaborées par analyse multicritère.

Discussion

La cartographie des zones de recharge des aquifères fracturés en contexte de variabilité climatique par analyse multicritère a certes donné des résultats satisfaisants, mais soulève des constatations, des incertitudes et quelques difficultés.