1Les perturbations climatiques observées ces dernières décennies, ne demeurent pas sans conséquences sur les terres agricoles de même que sur les ressources eaux. En effet, la variabilité du climat à travers ses principaux paramètres, notamment les précipitations et les températures, constitue aujourd’hui une menace potentiellement majeure pour l’environnement selon Totin et al, (2016). Les incidences d’événements climatiques extrêmes mettent en évidence la grande vulnérabilité et le degré élevé d’exposition de certains écosystèmes et de nombreux systèmes humains à la variabilité actuelle du climat (GIEC, 2014). Les communautés pauvres seront les plus vulnérables du fait de leurs capacités d’adaptation limitées et leur grande dépendance aux ressources à forte sensibilité climatique telles que l’eau et les systèmes de production agricole. Ces changements sont une réalité à l’échelle des terres émergées et des océans où la température moyenne est à la hausse et le réchauffement de 0,85 [0,65 à 1,06] °C (GIEC, 2013).
2Au Bénin, comme à l’instar des autres pays du monde, les ressources en eau et en sol sont sérieusement menacées par l’érosion (Mazour et Roose, 2002 ; Touaibia, 2010). L’eau et le sol, ressources vitales et de production, connaissent actuellement une dégradation en termes quantitatif et de qualitatif. Les déficits de pluviométrie, la fréquence et l’intensité des sécheresses et les pluies torrentielles accélèrent de plus en plus le stress hydrique et le phénomène érosif. L’avenir et la pérennité des lacs et des territoires ruraux est en jeu. L’érosion hydrique demeure une préoccupation majeure pour les pays en développement. Environ six millions d’hectares sont exposés aujourd’hui à une érosion active (Heddadj, 1997).
3Les pertes annuelles d’eaux, dues à l’envasement dans les barrages, sont estimées à environ 20 millions de m3 (Remini, 2000). Avec une érosion spécifique annuelle moyenne variant entre 2000 et 4 000 t/km‑2, le Bénin est classé parmi les pays où les terres agricole sont les plus exposées au phénomène érosif(Eétné, 2010). Les conséquences ne cessent de s’amplifier et se traduisent par un épuisement des ressources en eau et en sol et par l’aggravation des problèmes environnementaux.
4Sur le plateau de Zagnanado, les décennies 1970 et 1980 ont été marquées par une baisse considérable des précipitations (Idani et al., 2013) et une sécheresse sévère. Ces différentes crises climatiques poussent les populations à prendre d’assaut les berges des cours d’eau pour des activités agricoles. Ces pratiques ne demeurent pas sans conséquence sur la vie de ces écosystèmes fragiles.
5La présente étude vise à analyser les effets de la dynamique climatique sur les terres agricoles autour des petits lacs du plateau d’Abomey-Zagnanado au Bénin.
6Le secteur d’étude est constitué de l’espace des petits lacs (Azli et Sélé)du plateau de Zagnanado au Bénin (figure 1). Il est compris entre la latitude 6°53’ et 7°28 Nord et la longitude 1°58’ et 2°29’ Est. Le secteur des petits lacs du plateau de Zagnanado sont fortement influencée par le climat subéquatorial dont la moyenne pluviométrique annuelle est d’environ 1 137,40 mm avec un régime bimodal. Les pluies sont réparties sur deux saisons dont les maxima culminent à 186,2 mm et 147,4 mm respectivement en juin et octobre ; les températures moyennes mensuelles varient entre 29,5 °C au mois de mars et 26,9 °C au mois d’août (Atchadé, 2014). Durant la période 1952-2018, le nombre de jours de pluie a constamment baissé dans l’ordre de 11 à 25 % alors que la température a augmenté de 1 °C. Le régime hydrologique des deux lacs est fortement influencé par celui du fleuve Ouémé ; ainsi, les hautes eaux ont lieu en septembre et octobre, en liaison avec les apports du cours supérieur du fleuve, situé en climat soudanien au nord.
Figure 1. Situation géographiques des lacs Azli et Sélé
Sources : IGN, 2015 ; Atlas Sud-Bénin, 1992.
7La population du milieu d’étude est ancienne, avec une civilisation agraire exploitant alternativement les ressources du plateau et de la plaine d’inondation au cours de l’année (Houndagba et al., 2007). L’habitat y est groupé et concentré sur les versants de plateau et les bourrelets de berge. Les densités de population (versant et plaine d’inondation) dépassent très souvent 150 habitants au km², ce qui pose le problème de l’adéquation entre les ressources disponibles et les besoins de la population.
8Les précipitations et les températures mensuelles de la station synoptique de Bohicon et les hauteurs de pluie du poste pluviométrique de Zagnanado, collectées à l’ASECNA-Cotonou (Bénin) sur la période 1952-2018 ont été utilisées.
9Les données obtenues sur ces deux stations ont été régionalisées sur le secteur d’étude qui ne dispose pas d’un réseau d’observation bien fourni.
10A partir de l’hypothèse de l’homogénéité spatiale des précipitations mensuelles de Brunet-Moret (1971), les données manquantes ont été comblées par la méthode des doubles cumuls des stations du même domaine climatique. Soit pour une année i, les précipitations mensuelles xi et yi à deux stations x et y peuvent s’écrire :
11X et Y sont les moyennes inter-mensuelles sur une longue série des stations retenues. La variable εi est négligeable et nulle. Ensuite, une reconstitution de ces valeurs manquantes a été faite grâce à la régression linéaire multiple à partir de trois à quatre stations du domaine climatique les plus proches.
12Les données cartographiques concernent les images suivantes :
-
les images landsat TM (classifiées) de 2005 couvrant les lacs Azili et Sélé ;
-
les images landsat 8 (classifiées) de 2018 couvrant les lacs Azili et Sélé.
13La précipitation ou la température centrée réduite, (notée PCR) est calculée à l’aide de la formule :
avec ; xi pluviométrie ou température de l’année i ; X et σ respectivement la moyenne et l’écart-type de la pluviométrie ou la température moyenne interannuelle sur la période de référence.
14Ces indices de pluie centrée réduite (PCR), correspondant aux indices standardisés de précipitation de McKee et al. (1993), ont été utilisés pour déterminer les années excédentaires ou déficitaires. Cinq régimes de pluie ont été ainsi retenus en fonction des classes de PCR. Il s’agit de : régime très humide (PCR maximale), régime humide (1 < PCR < 1,49), régime normal (PCR = 0), régime sec (-1 < PCR < -1,49) et régime très sec (PCR minimale).
15Par ailleurs, le taux de variation (Tv) des totaux pluviométriques annuels des régimes extrêmes par rapport au régime normal est calculé par la formule :
avec Mextrême moyenne de pluie d’un quelconque régime extrême rapportée à un régime normal (Mnormale).
16Les indices calculés de température (TCR) ont permis de déterminer l’évolution thermométrique sur la période 1961-2018 et l’effet du pouvoir évaporant sur les ressources en eau.
17Elle permet de mettre en relief, le taux moyen annuel d’expansion spatiale, le taux de conversion et l’intensité des changements. Ces taux ont été calculés sur la base de la matrice de transition, ce qui a permis de comprendre les différentes mutations intervenues dans le temps autour des lacs Azili et Sélé entre 2005 et 2018.
18- Matrice de transition
C’est un tableau synthétique qui récapitule en termes de superficies, les différentes transformations subies par les formations végétales et les autres unités d’occupation des terres entre deux dates (Arouna, 2012). Elle a permis de décrire les changements d’état des unités d’occupation des terres entre 2005 et 2018.
19- Taux moyen annuel d’expansion spatiale
Le taux moyen annuel d’expansion spatiale exprime la proportion de chaque unité de végétation naturelle qui change annuellement. Ce taux annuel est calculé à partir de la formule suivante :
Avec s1 la superficie d’une unité de végétation à la date t1 et s2 la superficie de la même unité de végétation à la date t2 et t le nombre d’années d’évolution et ln logarithme népérien ; e la base des logarithmes népériens (e = 2,71828).
20- Taux de conversion
Le taux de conversion d’une classe de végétation correspond au degré de transformation subie par cette classe de végétation en se convertissant vers d’autres classes (Arouna, 2012). C’est donc la quantité de changements observés au niveau d’une formation végétale entre deux dates t0 et t1. Il permet ainsi de mesurer le degré de conversion d’une formation végétale en d’autres unités d’occupation des terres. Il est calculé à partir de la formule suivante :
Avec : TC : taux de conversion ; sit : Superficie de l’unité d’occupation des terres i à la date initiale t ; sis : Superficie de la même catégorie de cette unité i demeurée stable à la date t1.
21- Analyse de l’intensité des changements
Deux programmes d’analyse “PontiusMatrix41” et “Intensity Analysis03.xlms” de Aldwaik et Pontius (2012), ont permis de mesurer (en %) les intensités des changements selon les intervalles de temps, les unités d’occupation des terres et les transitions entre unités d’occupation des terres.
22- PontiusMatrix41
Le programme “PontiusMatrix41” fournit une variété d’analyses et les sommaires statistiques d’une matrice de tableau-croisé (www.clarku.edu/~rpontius). Ce programme “PontiusMatrix41” s’est basé sur les matrices de transition de 2005 à 2018 pour générer deux graphiques montrant lesdites intensités selon l’intervalle. Le premier graphique présente les intensités des changements survenus au sein de chaque catégorie, tandis que le deuxième graphique présente les intensités de pertes et de gains par catégorie. A ces niveaux, l’état de rapidité des changements est déterminé grâce à la ligne verticale en tireté bleu, appelée ligne de zone uniforme. Si le graphe est à gauche de cette ligne, le changement est lent ou dormant. Mais si celui-ci est à droite de cette ligne, le changement est rapide ou actif.
23- Intensity Analysis03.xlms
Le programme “Intensity Analysis03.xlms” grâce à la matrice de transition a permis de générer des statistiques pour les transitions selon les intervalles de temps, entre chaque catégorie d’occupation des terres et les autres. Il en est de même pour les pertes et gains survenus lors des transitions entre unités.
24La technique utilisée pour la collecte de ces données est celle du piquet d’érosion. Le piquet est constitué d’une règle en bois de forme parallélépipédique d’une longueur de 50 cm. Il comprend deux parties délimitées par une corde dont :
-
la partie inférieure du piquet à bout pointu permet de le fixer au sol sur une profondeur de 25 cm ;
-
la partie supérieure graduée en pas de 5 cm par des rayures tout autour de la règle.
25Pour faciliter l’identification des piquets sur le terrain, la partie supérieure est peinte en jaune sur une longueur de 5 cm. Le piquet est enfoncé dans le sol par le bout pointu à l’aide d’un marteau ou d’un maillet jusqu’à la première marque rouge. Les piquets sont installés en ligne droite suivant les facettes topographiques, mais avec un décalage de 4 m par rapport au piège à sédiments.
26La variabilité du climat autour des petits lacs du plateau de Zagnanado est marquée par une perturbation du régime pluviométrique et thermométrique dont dépend la disponibilité des ressources en eau.
27Les figures 2 et 3 montrent l’évolution de la pluviométrie sur la période 1952-2018 autour des lacs Azli et Sélé. La variabilité pluviométrique dans le secteur d’étude est marquée par des séquences de période humide (1952-1970), sèche (1971-1990) et d’instabilité (2011-2018) pour la station de Bohicon. Par contre, à la station de Zagnanado, les séquences humides couvrent la période (1952-1969), la période sèche (1970-1990) et une période de reprise (1991-2018). La fin de la période humide est en phase avec la rupture de stationnarité, située entre la fin de la décennie 1960 et le début de la décennie 1970 en Afrique occidentale et centrale, y compris le Sahel. Les conditions climatiques de cette période ont rendu difficile la production agricole par le manque d’eau pour la croissance des plantes.
Figure 2. Variabilité interannuelle de la pluviométrie à Bohicon sur la période 1952-2018
Figure 3. Variabilité interannuelle de la pluviométrie à Zagnanado sur la période 1952-2018
28La récurrence d’années sèches (anomalies négatives de -0,1 à -1,4), 53 % de la série d’étude, contrarie la proportion de 47 % d’années de pluie excédentaire (+0,1 à +4) et de 4 % seulement d’années normales pour les deux stations retenues pour la présente recherche.
29Ce résultat corrobore celui d’Atchade (2014) qui évoque une tendance à la persistance de la sécheresse progressive dans la partie centrale du Bénin aux confins du début du domaine de la transition climatique entre le subéquatorial et le soudanien. Les ménages agricoles (53 %) évoquent un bouleversement de la durée des saisons agricole dans ce domaine. Ce qui les oblige à une adaptation du système agricole à la nouvelle donne climatique.
30La dominance d’années déficitaires et le nombre très faible d’années normales (anomalies sensiblement égales à 0) permet d’apprécier ce qui se passe au niveau du régime de pluie. Les figures 4 et 5 présentent le régime de pluie autour des lacs Azli et Sélé sur le plateau de Zagnanado.
Figure 4. Régime de pluie à Bohicon sur la période (1952-2018)
Figure 5. Régime de pluie à Zagnanado sur la période (1952-2018)
31L’analyse des figures 4 et 5 montre que la hauteur de pluie annuelle moyenne est d’environ 1 328,96 mm en régime humide. A l’opposé, la pluviométrie annuelle moyenne est de 976,40 mm en régime sec. Au regard des conditions climatiques, le plateau de Zagnado connaît un régime pluviométrique subéquatorial (pluie > 1 200 mm).
32L’instabilité du régime pluviométrique compromet, dans une certaine mesure, l’intégrité des écosystèmes, source nécessaire pour la production agricole. Les cumuls annuels de pluie des régimes humide et sec sont très variables par rapport à celui du régime normal. Les taux sont de +55,5 % en régime humide à -44,5 % en régime sec, comparé à un régime climatique normal.
33A l’échelle saisonnière, les lames d’eau précipitées peuvent varier de 97 à 100 %. Ainsi, de la saison pluvieuse à la saison sèche, les hauteurs de pluie peuvent passer de 1 328,96 mm à 54,45 mm en régime humide, de 1 137,40 mm à 11,2 mm en régime normal, de 976,40 mm à 8,3 mm en régime sec. Cela dénote de la variabilité pluviométrique du milieu subéquatorial du Bénin et parfois un fort contraste saisonnier des précipitations.
34La combinaison de la baisse de la pluviométrie et de sa forte variabilité a un effet plus que proportionnel sur les systèmes de production, surtout autour des petits lacs du plateau de Zagnanado.
35L’analyse des figures 4 et 5 montre que la hauteur de pluie annuelle moyenne est d’environ 1 328,96 mm en régime humide. A l’opposé, la pluviométrie annuelle moyenne est de 976,40 mm en régime sec. Au regard des conditions climatiques, le plateau de Zagnado connaît un régime pluviométrique subéquatorial (pluie > 1 200 mm).
36L’instabilité du régime pluviométrique compromet, dans une certaine mesure, l’intégrité des écosystèmes, source nécessaire pour la production agricole. Les cumuls annuels de pluie des régimes humide et sec sont très variables par rapport à celui du régime normal. Les taux sont de +55,5 % en régime humide à -44,5 % en régime sec, comparé à un régime climatique normal.
37A l’échelle saisonnière, les lames d’eau précipitées peuvent varier de 97 à 100 %. Ainsi, de la saison pluvieuse à la saison sèche, les hauteurs de pluie peuvent passer de 1 328,96 mm à 54,45 mm en régime humide, de 1 137,40 mm à 11,2 mm en régime normal, de 976,40 mm à 8,3 mm en régime sec. Cela dénote de la variabilité pluviométrique du milieu subéquatorial du Bénin et parfois un fort contraste saisonnier des précipitations. La combinaison de la baisse de la pluviométrie et de sa forte variabilité a un effet plus que proportionnel sur les systèmes de production, surtout autour des petits lacs du plateau de Zagnanado.
38La tendance de la température annuelle moyenne est à la hausse (figure 6), de 26,35 à 28,75 °C entre 1952 et 2018 ; soit une augmentation de 0,03 °C/an (1,1 °C sur la période d’étude).
39A l’échelle mensuelle, les températures varient de 22,57 °C (décembre) à 23,54 °C (avril) pour les minima, de 29,45 °C (juillet) à 35,84 °C (février) pour les maxima et 25,59 °C (août) à 30,3 °C (février) pour les moyennes. En général, le mois de février est le plus chaud de l’année et constitue dans le milieu subéquatorial, celui de « présaison » avec des valeurs élevées d’évapotranspiration.
40Le renforcement du pouvoir évaporant de l’air en février correspond à une pression climatique sur les ressources en eau du milieu d’étude. Cette influence du climat suite à la sécheresse absolue des mois de janvier-février induit une intensification des activités agricoles (65 % des enquêtés) surtout celle des cultures de contre saison.
Figure 6. Variabilité interannuelle des anomalies de température (1952-2018).
41L’observation des cartes d’occupation des terres en 2005 (figure 7 et 8) montre que les formations végétales naturelles dominantes sont les savanes arborées et arbustives qui s’observent principalement dans les parties nord, nord-est et nord-ouest du lac Azli. Cette même formation végétale naturelle est observée également vers le nord-ouest et le sud-ouest du lac Sélé. Ces formations végétales sont suivies des formations marécageuses tout autour du lac Azli et les cultures et jachères sous palmiers et les formations marécageuses surtout dans la partie sud et ouest du lac Sélé. Les minces forêts galeries du bassin sont situées le long des différents cours d’eau dont le plus important est l’Ouémé. Les agglomérations sont très faiblement représentées au sud du lac Azli et à l’ouest et sud du lac Sélé.
Figure 7. Dynamique des unités d’occupation des terres autour du lac Azili entre 2005 et 2018
Figure 8. Dynamique des unités d’occupation des terres autour du lac Azili entre 2005 et 2018
42En 2018 (figure 7), les savanes arborées et arbustives ont connu une régression très significative (0,95 %) autour des deux lacs. On remarque la présence des plantations au niveau du lacs Azli et la présence des forêts claires et savanes boisées de même que la forêt galerie au niveau du lac Sélé.
43Les savanes arborées et arbustives sont devenues les formations naturelles les plus répandues autour de deux lacs. Par contre, les mosaïques de cultures et jachères ont connu une forte extension (65 %) elles se rencontrent un peu partout autour des lacs. Les sols dénudés et les plantations se sont étendus.
44Le tableau 1, présente la matrice de transition des unités d’occupation des terres autour des lacs Azli et Sélé entre 2005 et 2018.
Tableau 1 - Matrice de transition des unités d’occupation des terres autour du lac Azili
Forêt galarie (FD), Forêt dense (FD), Forêt claire et savane boisée (FCSB), Savane arborée et arbustive (SAA), Plantation (P), Formation marécageuse (FM), Champs et jachère (CJ), Champs et jachères sous-palmeraie (CJSP), Plan d’eau (PE) et Agglomération (A).
45Entre 2005 et 2018, les FG, FD, FCSB et CJP n’ont pas connu de mutation autour du lac Azli.
46Les SASA qui occupaient une superficie de 4 165,52 ha ont connu une régression de 542,07 ha au profit des CJ ; 106,90 ha au profit des FM et 6,89 ha au profit des HA avec 3 509,66 ha de superficie restée stable entre 2005 et 2018. Le taux de conversion et le taux moyen annuel d’expansion spatiale des SASA sont respectivement 15,74 % et – 0,60 %. Les PL ont connu une évolution régressive (1,97) de leur superficie au profit des CJ entre 2005 et 2018 avec un taux de conversion de 11,66 % et un taux moyen annuel d’expansion spatiale de 3,44 %. Les CJ ont connu une évolution régressive de leur superficie au profit des SASA (342,69 ha), des PL (11,48 ha) et des HA (5,90 ha) avec 1 336,44 ha de superficie resté stable. Le taux de conversion des CJ et le taux moyen annuel d’expansion spatiale sont respectivement 21,22 % et 0,79 %. Les FM qui occupaient une superficie de 1 529,46 ha en 2005 ont connu une évolution régressive (6,07 ha) de leur superficie au profit des PE avec un taux de conversion de 0,40 % et un taux moyen annuel d’expansion spatiale de 0,49 %. Les PE occupaient une superficie de 429,91 ha en 2005, ont connu une évolution progressive de leur superficie et occupent 435,98 ha en 2018 avec un taux moyen annuel d’expansion spatiale de 0,11 %.
47Entre 2005 et 2018, les HA sont passées de 26, 73 ha à 39,52 ha avec un taux moyen annuel d’expansion spatiale de 3,01 %.
48Le tableau 2, montre la transition qui avait eu lieu dans les unités d’occupation des terres autour du lac Sélé entre 2005 et 2018.
Tableau 2. Matrice de transition des unités d’occupation des terres autour du lac Sélé
Forêt galarie (FD), Forêt dense (FD), Forêt claire et savane boisée (FCSB), Savane arborée et arbustive (SAA), Plantation (P), Formation marécageuse (FM), Champs et jachère (CJ), Champs et jachères sous-palmeraie (CJSP), Plan d’eau (PE) et Agglomération (A).
49Entre 2005 et 2018, les unités d’occupation des terres telles que les FG, FD, FCSB, CJP et PL n’ont pas connu de mutation autour du lac Sélé.
50Les SASA par contre ont connu une régression de 327 ha et de 6,41 ha respectivement au profit des CJ et HA avec 2 406,02 ha de superficie restée stable. Le taux de conversion et le taux moyen annuel d’expansion spatiale des SASA est respectivement de 12,17 % et de – 0,87 %. Les CJ ont connu une mutation de leur superficie au profit des SASA (39,10 ha), des PE (5,50 ha) et des HA (45,21 ha) avec 2 856,73 ha de superficie resté stable. Le taux de conversion de cette unité et le taux moyen annuel d’expansion spatiale sont respectivement 3,05 % et 0,73 %. Les FM qui occupaient une superficie de 953,06 ha en 2005 ont connu une évolution régressive de leur superficie au profit des CJ (40,02 ha) et PE (1,07) avec un taux de conversion de 4,31 % et un taux moyen annuel d’expansion spatiale de -0,21 %. En 2005, les PE occupaient une superficie de 376,49 ha. 15,88 ha se sont convertis en FM et 15,43 ha en CJ avec 345,18 ha de superficie restée stable. Le taux de conversion et le taux moyen annuel d’expansion spatiale sont respectivement 8,32 % et -0,52 %.
51Entre 2005 et 2018, les HA sont passés de 65,22 ha à 116,84 ha avec un taux moyen annuel d’expansion spatiale de 4,48 %.
52Les figures 9 et 10 présentent l’intensité des changements des unités d’occupation des terres entre 2005 et 2018 autour des Lacs Azli et Sélé.
Figure 9. Intensité des changements des unités d’occupation des terres entre 2005 et 2018 autour du Lac Azli
Figure 10. Intensité des changements des unités d’occupation des terres entre 2005 et 2018 autour du lac Sélé
53L’observation de la figure 9 indique que les FG, FD, FCSB et CJP n’ont enregistré aucun gain ni perte entre 2005 et 2018 autour du lac Azili. Par contre, les SASA sont les unités qui ont le plus connu de mutation en termes de perte de superficie alors que les CJ ont enregistré le plus de gain de superficie. L’intensité des changements en gain et perte de superficie au niveau des autres unités d’occupation des terres est faible autour du lac Azli( perte entre 2005 et 2018). Par contre, les SASA sont les unités qui ont le plus connu de mutation en termes de perte de superficie au profit des CJ, PE et FM. L’intensité des changements en gain et perte de superficie au niveau des autres unités d’occupation des terres est faible autour du lac Sélé.
54Les figures 11a et 11b met en relief l’intensité et la vitesse des changements d’occupation des terres entre 2005 et 2018 autour des petits lacs du plateau de Zagnanado.
Figure 11a. Intensité et vitesse des changements des unités d’occupation des terres entre 2005 et 2018 autour du lac Azli.
Figure 11b. Intensité et vitesse des changements des unités d’occupation des terres entre 2005 et 2018 autour du lac Sélé.
55Au niveau de la figure 11a, il est montré que des unités d’occupation FG, FD, FCSB, CJP, FM et PE les gains et pertes sont dormants. Par contre au niveau des PL, HA et CJ l’intensité et la vitesse des changements en termes de gains sont rapides. On observe la même tendance pour les pertes aux seins des SASA et des CJ autour du lac Azli.
56L’examen de la figure 11b montre qu’au niveau des FG, FD, FCSB, CJP et PE les gains et pertes sont dormants. Au niveau des SASA, PL et CJ l’intensité et la vitesse des changements en termes de gain et perte de superficie sont rapides. Par contre, on observe des gains rapides aux seins HA et des FM autour du lac Sélé.
57Au total, il est retenu que l’analyse de la dynamique de l’occupation des terres entre 2005 et 2018 est marquée par la régression de toutes les formations végétales naturelles au profit des unités d’occupation anthropiques. Entre ces deux dates, les champs et jachères ont connu la plus grande expansion (72 %). Ce qui contribue énormément à la dégradation des terres agricoles.
58La dynamique érosive est observée le long des transects établis dans la zone d’étude. Les quantités des matières en mouvement des sommets vers les bas de pentes varient selon la position topographique, la nature des roches et la force de l’agent érosif. Les tableaux 3 et 4 renseignent sur la quantité des terres érodées le long de deux écosystèmes lacustres.
59Orienté NE-SW du lac Sélé, le transect établi est long de 900 m. Le transect présente 3 versants dont les pentes sont respectivement 7,4 et 2 %, un sommet avec une pente de 1 % et un bas de pente de 0,5 % (Tableau 3).
Tableau 3. Relevé des piquets d’érosion à Sélé
Source : Travaux de terrain, juillet, 2019.
60Il résulte des données du tableau 3 que les versants sont marqués par un phénomène érosif très important (12 cm en 2018 et 15,26 cm en 2019) (Photo 1).
Photo 1. Phénomène de rigole et du sapement dans un champ de maïs sur versant du lac Sélé
Prise de vue : Etene, juillet 2020.
61Ce départ se traduit surtout par l’importance de la pluviométrie observée ces années. Par contre, l’érosion est moins importante dans le bas de pente et le sommet de l’interfluve, ce qui permet d’expliquer que la position topographique intervient en grande partie dans le départ des sédiments. Le ruissellement étant la principale cause de l’érodibilité du sol, les forts mouvements des matériaux sont généralement observés au niveau des piquets où la végétation est faible ou quasi nulle.
62Situé presque en amont du lac Azli et orienté NE-SW ce transect à une longueur 500 m (Tableau 4).
Tableau 4. Relevé des piquets d’érosion à Azli
Source : Travaux de terrain, juillet, 2019
63Ledit transect présente dans l’ensemble un sommet d’interfluve suivi d’un long versant légèrement convexe qui draine les eaux pluviales vers le lac. Les pentes au niveau des différentes facettes topographiques sont comprises entre 0,5 à 1 % sur le sommet, 1 à 6,4 % sur le versant et 0,1 à 0,3 % dans le bas de pente. L’occupation humaine en termes d’espace agricole de ce sous bassin est estimée à 70 %. Les plus importants mouvements de particules sont observés au niveau des piquets 2 et 3 qui sont installés sur les versants. Cet important départ des sédiments s’explique par la présence d’une forte pente, d’une surface imperméabilisée et par l’absence du couvert végétal.
64Toutefois, il est important de rappeler que pendant la saison pluvieuse, les deux lacs bénéficient énormément des apports hydriques et sédimentologiques substantiel de l’Ouémé, puis du Zou. Ces apports peuvent être bénéfiques mais aussi dramatiques pour les champs et les cultures (photo 2).
Photo 2. Inondation d’un champ de maïs du lac Sélé
Prise de vue : Etene, juillet 2020.
65La connaissance du comportement possible du climat d’une région donnée, en cas de perturbations climatiques, est un élément de base pour toute évaluation de la vulnérabilité de cette localité. La variabilité climatique dans le secteur d’étude est marquée par des séquences humides (1952-1970), séquences sèches (1971-1990) et la période d’instabilité (2011-2018) sur le plateau de Zagnanado où sont localisés les lacs Azli et Sélé.
66Les pluies annuelles en régime normal/moyen sont de 1 157 mm à 1 262 mm et confirment l’appartenance de cette région à un domaine climatique subéquatorial comme l’ont montré Boko (1988) et Houndénou (1999). L’instabilité du régime pluviométrique compromet, dans une certaine mesure, l’intégrité des écosystèmes, substrat pour les activités agricoles. Les cumuls annuels de pluie des régimes irréguliers (sec et humide) sont très variables par rapport à celui du régime normal.
67En effet, le plateau de Zagnanado a connu entre les années 1971 et 1990 une phase de sécheresse qui peut être interprétée comme une manifestation de la grande sécheresse qui a frappé l’Afrique de l’Ouest à partir de la fin des années 1960 (Bricquet et al., 1997). Ce déficit pluviométrique a été également mis en évidence, dans le Bénin méridional et central (Ogouwalé, 2006) et dans le bassin du Zou (Atchade, 2014).
68La période des années 70, qui est caractérisée par la séquence sèche, coïncide également avec l’époque de la révolution au Bénin, qui a fragilisé les réglementations traditionnelles en matière de gestion et de protection des écosystèmes. Ainsi, le point de départ de la régression rapide des lacs devrait être situé au début des années 70 avec la conjonction de la sécheresse ayant entraîné la destruction de nombreuses forêts sacrées(Houndagba, 2007).
69Ce résultat est conforme au résultat de Mâche et Olivry (1995), en étudiant les variations des précipitations et des écoulements entre 1951 et 1989 en Afrique de l’ouest et du centre dans 34 bassins hydrographiques dont ceux de l’Ouémé et du Mono au Bénin. Ces auteurs ont, en effet, montré que le déficit pluviométrique persistant des années 70 et 80 a engendré un tarissement progressif des fleuves, particulièrement prononcé au cours des années 80 en Afrique de l’ouest, et qui s’est traduit par une diminution des coefficients d’écoulement. Au total, le tarissement des lacs Azliet Sélé, mis en évidence par les données de télédétection, doit être considéré comme un phénomène à caractère régional, amplifié par les facteurs locaux (Houndagba et al, 2007).
70Les cartes d’occupation des terres réalisées ont permis de mettre en évidence la tendance générale de la dynamique de l’occupation des terres. Cette tendance générale révèle une augmentation de la superficie des formations anthropiques (cultures et jachères, agglomérations, plantations) au détriment des formations végétales naturelles (forêts, savanes boisées). Atchade (2014) et Houndagba et al, (2007) étaient déjà parvenus à la même conclusion concernant le bassin versant de la rivière Zou et l’étude de la dynamique des petits bassins lacustres du Sud-Bénin. Plus inquiétant, cette tendance de dégradation de la couverture végétale n’épargne même pas les terres agricoles dans le milieu d’étude.
71L’analyse de l’étude de la dégradation des terres a montré que les versants sont marqués par un phénomène érosif très important (12 cm en 2018 et 15,26 cm en 2019) au niveau du lac Sélé. Par contre, au niveau du lac Azli, l’occupation humaine en termes d’espace agricole est estimée à 70 %. Les plus importants mouvements de particules sont observés au niveau des versants. Ce résultat peut être comparé à celui de Boutkhil et al., (2013) et Morsliet al., (2005)qui ont montré que la dégradation est le résultat des incendies et des défrichements qui se sont succédés tout au long de l’histoire. Suite à cette dégradation, le sol a subi à son tour une modification. Situé sur pente forte et occupé par une végétation clairsemée, le sol est très affecté par l’érosion en nappe. Il en résulte un décapage des horizons de surface, une charge colluviale importante de 30 %, un tassement du sol et une diminution d’activité biologique. Le décapage du sol est si important que l’érosion a fini par affleurer la roche‑mère gênant ainsi fortement l’infiltration des eaux.
72La présente recherche a permis de caractériser la dynamique climatique et l’évolution de l’occupation des terres autour des petits lacs sur le plateau de Zagnanado au Bénin. Sur le plan pluviométrique, à partir de la chronique des précipitations annuelles (1952-2018), on note que le secteur d’étude est caractérisé par une succession de phases humides (47 %) et sèches (53 %). La tendance de la température moyenne est à la hausse de 26,35 à 28,75 °C entre 1952 et 2018 ; soit une augmentation de 1,1 °C sur la période d’étude.
73Le déficit pluviométrique lié à la fréquence des années sèches, a entraîné l’accentuation des pressions humaines sur les écosystèmes des zones humides, l’aggravation des phénomènes érosifs, dont les conséquences sont la dégradation des terres agricoles.
74L’analyse de la dynamique de l’occupation des terres entre 2005 et 2018, à partir des cartes d’occupation des terres, a montré que cette dynamique est marquée par la régression de toutes les formations végétales naturelles (forêt galerie, et savane boisée) au profit des unités d’occupation anthropiques (mosaïque de champs et jachères, plantation, agglomération). Ce qui ne demeure pas sans conséquence sur les terres agricoles.