TROPICULTURA, 201 7, 35, 4, 31 2-324

Modélisation de la dynamique du paysage forestier de la Réserve

Spéciale d’Ambatovaky (Nord-Est de Madagascar)
L.O. Rakotondrasoa 1 * , F. Malaisse 2 & J. Bogaert2
Keywords : Landscape dynamic- Forest- Ambatovaky- Madagascar
Résumé Summary
Madagascar abrite une biodiversité unique qui doit Modelling of the Forest Landscape Dynamics
faire face à de nombreuses menaces, dont of Réserve Spéciale d’Ambatovaky
l’anthropisation, entraînant la déforestation. Cette
étude a pour objectif de déterminer la dynamique (Northeast of Madagascar)
paysagère des forêts de la Réserve Spéciale Madagascar shelters a unique biodiversity that faces
d’Ambatovaky de 1996 à 2050 à partir d’une chaîne numerous threats, such as human activities, which
de Markov de premier ordre. Pour ce faire, des cause deforestation. The objective of this study was
images Landsat datées de 1996, 2004 et 2014 ont to identify the landscape forest’s dynamics of the
été utilisées pour l’élaboration des cartes Réserve Spéciale d’Ambatovaky (RSA) from 1996 to
d’occupations de sol multi-dates. Une matrice de 2050 based upon a first order Markov chain model.
transition a ensuite permis de simuler la composition For this purpose, LANDSAT images of 1996, 2004
du paysage jusqu’en 2050. Les résultats ont montré and 2014 were used to compose the land cover
que la superficie forestière de la Réserve diminue au maps. A transition matrix made it possible to simulate
profit d’une formation secondaire appelée « savoka» the composition of the landscape until 2050. The
qui apparaît quelques années après la mise en results showed that the forest area of RSA decreased
culture temporaire. La culture sur brûlis, appelée and is substituted by a secondary formation called
localement «tavy», une pratique commune des “savoka” that appears several years after the end of
paysans de la côte malgache, constitue la principale agricultural activities. Slash-and-burn agriculture,
raison de la régression de la superficie forestière. called locally "tavy", is a common practice of farmers
in the coastal areas of Madagascar. Tavy is the main
reason for the decline of the forest area.

1 Universitéd’Antananarivo, Ecole Supérieure des Sciences Agronomiques, Département des Eaux et Forêts, Antananarivo, Madagascar.
2 Universitéde Liège, Gembloux Agro-Bio Tech, Gembloux, Belgique.
*Auteur correspondant: Email: [email protected]
31 2
Reçu le 23.03.1 7 et accepté pour publication le 1 9.06.1 7
 TROPICULTURA, 201 7, 35, 4, 31 2-324

Introduction Notre recherche se situe dans ce cadre; elle a pour

objectif d’étudier la dynamique spatio-temporelle
Madagascar est connu mondialement pour la d’une aire protégée de Madagascar, à savoir la
richesse de sa biodiversité terrestre et marine unique Réserve Spéciale Ambatovaky (RSA), qui constitue
et exceptionnelle. Plus de 1 2.000 espèces végétales une des plus grandes réserves spéciales malgaches
y sont répertoriées, dont le degré d’endémisme a été ayant une vaste forêt de basse altitude. Notre étude
estimé à des valeurs comprises entre 82 et 90% (1 4, vise également à dégager la tendance évolutive de la
1 5, 21 ). En outre, le taux d’endémisme de la faune RSA à l’horizon 2050, en termes d’occupation du sol.
est également élevé, étant de l’ordre de 80% (6). A cette fin, nous avons retenu comme hypothèse de
Mais, parallèlement à cette richesse, il faut prendre travail que les pratiques culturales, notamment
conscience d’une poussée démographique l’agriculture itinérante sur brûlis ou «tavy», favorisent
préoccupante avec un taux d’accroissement de la régression du couvert forestier du paysage
2,62% en 201 4 (34). À l’échelle mondiale, la d’Ambatovaky, malgré son statut d’aire protégée. Le
déforestation a plus ou moins suivi le rythme de la « tavy » est un terme malgache pour désigner le
croissance démographique (1 2). En effet, les causes déboisement de la forêt primaire par abattage suivi de
de la déforestation sont en relation avec brûlis en vue de la culture du sol (1 3). Une croissance
l'augmentation de la croissance démographique, démographique non accompagnée d’une croissance
entre autres par la nécessité de plus de terres pour la économique pourrait être fatale à la situation de la
production agricole. La poussée démographique est population, menaçant la pérennisation des ressources
également la cause principale de la déforestation en naturelles. En effet, pour maximiser les productions
Afrique et certainement à Madagascar (1 0). La agricoles, en réponse à l'accroissement
déforestation de la Grande Ile résulte de démographique très rapide, l'agriculture sur brûlis
défrichements et de la pratique de cultures itinérantes reste la principale technique de production vivrière
sur brûlis et de l’exploitation forestière en vue de dans la région. Après quelques années (2 à 5 ans) de
fournir les besoins en bois d’œuvre et de production, la parcelle est laissée à l’abandon et la
construction, ainsi que la production de bois d’énergie population va défricher une autre partie de la forêt. La
de la population locale (4, 6, 1 0, 1 9, 28). Le taux végétation secondaire qui se développe sur le
annuel de déforestation de la Grande Ile était de « tavy » abandonné, pour quelques années, est
0,3% de 1 950 à 1 970; il a augmenté et a atteint 1 ,7% dénommée «savoka» (1 3). Ce dernier peut être aussi
pour la période 1 970 à 1 990 (1 6). l’objet d’un nouveau défrichement par la
Face à la menace qui pèse sur la biodiversité de population locale.
Madagascar, le gouvernement malgache a établi en
1 990 la première version de la Charte de Matériels et méthodes
l’Environnement qui a pour finalité d’enrayer la spirale
de dégradation de l’environnement (26). La création Zone d’étude
des aires protégées figure parmi les solutions Ambatovaky est située dans la partie Nord-Est de
entreprises. Les aires protégées, gérées par Madagascar (Figure 1 ), elle appartient au district de
Madagascar National Parks (MNP), jouent plusieurs Sonierana Ivongo et à la Région administrative
rôles à la fois: lieux de conservation de la d’Analanjirofo. La zone est soumise à un climat
biodiversité; lieux d’éducation relative à tropical humide de type perhumide chaud. Ce type de
l’environnement; lieux de recréation et de découverte; climat se caractérise par l'absence de saison sèche.
et lieux de recherche et pôles de développement Les précipitations moyennes annuelles s’élèvent à
économique. Les stratégies de conservation de MNP 2000 mm réparties sur 1 80 jours, quant à la
sont basées et développées de manière température, ses valeurs moyennes mensuelles sont
systématique et scientifique selon un Plan de Gestion comprises entre 21 °C 87 (février) et 24 °C
et de Conservation (PGC). Le PGC évalue les (septembre) (1 8).
menaces sur les cibles de conservation, les moyens La RSA a été créée en 1 958. Cette Aire Protégée,
et capacités de gestion. Ces informations sont gérée par Madagascar Nationals Parks (MNP) est
ensuite utilisées pour développer la stratégie de classée dans la catégorie IV selon l’UICN, c'est-à-
conservation, dont le système de surveillance et de dire, aire protégée principalement à des fins de
contrôle forme l’outil principal de la conservation. conservation, avec intervention au niveau de la
Suite à la mise en place des aires protégées, le taux gestion. Elle est localisée entre les points 1 6°43’21 ’’
annuel de déforestation malgache a diminué et a été et 1 6°43’58’’ de latitude Sud et 49°01 ’07’’ et 49°1 5’06’’
de 0,41 % entre 1 990 et 201 0 (1 3); spécifiquement, il de longitude Est, et s’étend sur environ 78 000 ha
est de 1 ,7% entre 2005 à 2011 pour la forêt humide (1 8). La zone présente un relief accidenté, avec des
du Nord-Est (32). Cependant, ces valeurs sont pentes abruptes localement. L’altitude est comprise
considérables en regard du taux mondial entre 31 4 m au niveau de Sandragato et
qui est de 0,1 4% entre 1 990 et 201 0 (1 3). 11 75 m àl’ouest, au mont Ambatovaky.

31 3
 TROPICULTURA, 201 7, 35, 4, 31 2-324

Figure 1 : Localisation de la Réserve Spéciale d’Ambatovaky (Projection WGS 84).

Le réseau hydrographique de la réserve est Méthodologie

important. La réserve est principalement caractérisée L’image satellite est une source d’informations
par une Forêt Dense Humide Sempervirente (FDHS); abondantes et précises sur l’occupation du sol et sur
elle constitue une excellente représentation de les phénomènes qui se déroulent à la surface de la
l’écorégion de la forêt humide de l’Est. Cette dernière planète (8). Trois images satellitaires, datées de
abrite une biodiversité originale très élevée et 1 996, 2004 et 201 4, ont été utilisées dans le cadre de
possède un taux d’endémicité important. 291 cette étude (tableau 1 ).
espèces de plantes y ont été recensées, dont 221 Ces images LANDSAT ont été téléchargées via le
sont endémiques à Madagascar (1 8). Ces 291 portail Earth Explorer (http://earthexplorer.usgs.gov/)
espèces sont réparties dans 1 59 genres et 59 de l’USGS (United States Geological Survey). Elles
familles. Notons l’existence de 5 familles endémiques se trouvaient en LPGS (Level 1 Product Generation
de Madagascar: Asteropeiaceae, Melanophyllaceae, System) et contiennent encore des erreurs
Physenaceae, Sarcolaenaceae et géométriques et radiométriques. Par conséquent, il
Sphaerosepalaceae. La réserve abrite 31 espèces de leur a été appliqué des corrections radiométriques
mammifères dont 11 lémuriens, 1 3 insectivores, 3 systématiques et d’amélioration des précisions
carnivores et 4 rongeurs; 11 0 espèces d’oiseaux, 11 3 géométriques en utilisant des Ground Control Points
espèces Herpetofauniques puis 34 espèces de provenant de GLS 2000 (Global Land Survey 2000) et
poissons (1 8). d’ajustement topographique par le biais de l’utilisation
d’un DEM (SRTM) (33).

31 4
 TROPICULTURA, 201 7, 35, 4, 31 2-324

Tableau 1
Liste des images utilisées et leurs caractéristiques

Des prétraitements ont été opérés afin de réajuster Sur le tableau de la matrice de confusion, les erreurs
les valeurs radiométriques et spectrales et de de classifications sont totalisées. Il existe deux types
minimiser les variabilités atmosphériques tenant d’erreurs: les erreurs d'omission dans le cas où des
compte la nature multi date des images utilisées. Les points de référence d'une catégorie sont attribués à
calibrations ont été ainsi appliquées sur les bandes une autre catégorie; et les erreurs de commission
multispectrales afin de transformer les valeurs DN dans le cas contraire où les points d'une catégorie lui
(Digital Number) en des unités radiance et ToA (Top sont attribués par erreur. En complément de ces deux
of Atmospheric) réflectance. Le DOS (Dark Object types d’erreur, la valeur du coefficient kappa (κ) a été
Subtraction) a été par la suite implémenté afin de calculée (25). Ce coefficient mesure l’intensité de la
ramener les valeurs spectrales à zéro ou en d’autres concordance réelle entre deux jugements catégoriels
termes de réduire les effets de diffusion appariés; il est proche de 0 lorsque la concordance
atmosphérique à un seuil minimal (9). est peu significative et de 1 lorsque la
Par la suite, une classification supervisée a été concordance est très forte.
effectuée avec l’image datée de 201 4, en utilisant le Les caractéristiques spectrales des classes
l’algorithme du maximum de vraisemblance via le d’occupation du sol obtenues suite à la classification
logiciel de traitement d’image ENVI. C'est une finale de l’image Landsat 201 4 ont été utilisées
classification par zones d'entraînement comme base pour une classification supervisée par
représentatives formées par les points relevés sur maximum de vraisemblance des images datées de
terrain des principaux types d’occupation du sol, 2004 et de 1 996. Sur base des cartes multi-dates
respectivement 1 22, 87 et 69 observations pour les obtenues, la simulation de l’évolution future des
classes «forêt», «savoka» et «champs de culture». classes d’occupation de la zone d’étude a été réalisée
En effet, dans cette analyse de l’évolution globale de à l’aide de l’élaboration d’une matrice de transition (5,
l’occupation du sol, trois classes thématiques ont été 27) et de l’utilisation d’une chaîne de Markov (29, 30).
retenues. Il s’agit: d’une classe «forêt» (Photo 1 ) La matrice de transition a servi à construire une
regroupant les forêts naturelles et les forêts pour les matrice de probabilités annuelles qui représente le
cultures de sous-bois, notamment pour la vanille; changement d’état d’une classe sur l’horizon d’un an
d’une classe «savoka» (2) (Photo 2) associant les (30). Les taux des superficies des classes du
végétations secondaires de quelques années et paysage en 1 996 ont été transcrits en un vecteur. En
d’une classe «champ de culture», constituée par des considérant l’occupation du sol comme un processus
rizières, des champs de manioc et patate douce entre stochastique (22) et les différents changements
autres et les jachères d’un à deux ans (Photo 3). comme des états de cette chaîne (31 ), le produit de la
Les jeunes jachères ont été incluses dans la classe matrice de probabilités annuelles et du vecteur
«champs de culture» car elles comportent représentant la composition du paysage d'une année
essentiellement des plantes herbacées, entre autres particulière en utilisant une chaîne de Markov de
Lantana camara, Aframomum angustifolium, Rubus premier ordre permet de simuler la superficie des
sp. (1 7), et par conséquent, elles ont été souvent classes du paysage à une date ultérieure. La
confondues avec les champs de culture au cours de simulation obtenue a été vérifiée et validée en
la classification de l’image Landsat. comparant la composition paysagère simulée avec
La carte obtenue a été vérifiée au cours d’une celle observée respectivement sur les images 2004 et
mission de validation sur le terrain en mois d’avril 201 4 par l’intermédiaire des tests χ2. L’année 2050 a
201 5 permettant d’élaborer une matrice de confusion été choisie comme année de projection dans la
et de calculer l’indice Kappa, pour estimer l’exactitude simulation. Ce choix a été basé sur les objectifs
de la classification. La matrice de confusion est un d’Aïchi (7), qui est un plan stratégique en faveur de la
tableau à double entrée dans lequel sont confrontés biodiversité stipulant que «d’ici 2050, la diversité
la réalité sur le terrain et le résultat de la biologique est valorisée, conservée, restaurée et
classification. Il permet d'évaluer l'intensité de la utilisée avec sagesse, en assurant le maintien des
liaison entre des données de référence et le résultat services fournis par les écosystèmes, en maintenant
de la classification. la planète en bonne santé et en procurant des
avantages essentiels à tous les peuples». 31 5
 TROPICULTURA, 201 7, 35, 4, 31 2-324

(Rakotondrasoa, 201 5)
Photo 1 : Forêt Dense Humide Sempervirente de la Réserve Spéciale
d’Ambatovaky.

(Rakotondrasoa, 201 5)
Photo 2 : «Savoka» à Ambavala (formations secondaires apparaissant
quelques années après la mise en culture temporaire).

31 6
 TROPICULTURA, 201 7, 35, 4, 31 2-324
pas significativement différentes entre les années
Résultats 1 996 et 2004 (χ2= 28,67; p>0,05). Par contre, la
Matrice de confusion et Indice de Kappa relative à différence des distributions est significative entre les
la classification de l’image de 201 4 années 2004 et 201 4 (χ2=1 32,29; p<0,05): il y a eu
L’évaluation de la classification de l’image de 201 4 donc un changement significatif pour la deuxième
est représentée par la matrice de confusion (Tableau période. Le changement majeur d’occupation du sol
2). Les valeurs en diagonale indiquent le nombre de de la Réserve au cours de la période 2004-201 4 a été
pixels bien classés (la classification effectuée sur une augmentation d’environ 47% de la superficie des
ordinateur correspond à la réalité sur terrain) et ceux «savoka» ainsi qu’une augmentation de 43% de la
hors diagonale le nombre de pixels mal classés pour classe «champs de culture » qui sont respectivement
chaque unité d'occupation du sol. Le rapport de la passés de 4489 ha à 6592 ha et de 1 053 ha à 1 501
somme des pixels bien classés sur le total de pixels ha.
utilisés dans la classification nous donne le Pour la totalité de la période d’étude (1 996-201 4), les
pourcentage de classification global qui est de l'ordre distributions de la proportion des classes d’occupation
de 84,1 %. En effet, sur les 1 26 points GPS collectés du sol présentent une différence significative (χ2=
sur terrain, 1 06 correspondent à la classification. 31 ,92; p<0,05) : on assiste à une nette régression de
L’analyse de cette matrice permet aussi d’avancer la superficie forestière (20753 ha en 1 996 et 1 71 31
que le niveau de précision de la classe «savoka» a ha en 201 4, soit une perte de 1 7,5% correspondant à
été le plus faible: sur les 57 points GPS collectés un taux annuel de déforestation de 0,9%). Si on
dans cette catégorie, 45 correspondent à la raisonne selon deux périodes, les taux de
classification de l’image satellite, soit 79%; le reste a déforestation annuels de la RSA sont de 0,6% entre
été confus aux deux autres types d’occupation du sol. 1 996 et 2004 et de 1 ,3% entre les années 2004 et
Ce niveau de précision serait dû à l'hétérogénéité 201 4.
structurale et spécifique de la physionomie du Dynamique du paysage
«savoka» en fonction de son âge. En effet, on La figure 3 illustre les changements observés de 1 996
distingue trois types de savoka selon leur âge, à à 201 4 dans une partie de la zone d’études.
savoir: (i) savoka jeune de deux à cinq ans qui peut Les changements d’occupation de sol survenus dans
atteindre deux à trois mètres de hauteur avec RSA pendant les deux périodes 1 996-2004 et 2004-
quelques arbustes; (ii) savoka d’âge moyen de cinq à 201 4 sont présentés dans les tableaux 4 et 5. Chaque
dix ans à l’intérieur duquel les arbustes commencent valeur des tableaux correspond à la proportion de
à dominer en éliminant peu à peu les plantes l’aire convertie de la classe indiquée sur la ligne vers
herbacées; et (iii) savoka vieux de plus de 1 0 ans la classe en tête de colonne; et les valeurs en gras
avec une dominance des arbustes et de Ravenala indiquent les taux de permanences des classes. A
madagascariensis et d’autres espèces pionnières titre d’exemple (pour la première ligne du tableau 4):
telles que Harungana madagascariensis et Albizzia sur 82,3% de la surface du paysage qu’occupait la
lebbeck (1 7). forêt en 1 996, 75,5% sont restés intacts; 5,5% ont été
L'indice de Kappa est un indice qui permet de dégradés et convertis en «savoka» et 1 ,3% convertis
«retirer» la portion de hasard ou de subjectivité de en champs de culture.
l'accord entre les techniques. Après calcul, l’indice Selon les tableaux 4 et 5, la forêt s’est révélée être la
Kappa est égal à 82%. D’après Pontius, une étude de classe qui a subi la plus grande conversion, plus
l’occupation du sol peut être validée si κ> 50% (25); précisément des régressions. En effet, sur les 82,3%
cela permet de conclure que le résultat de la qu’occupait la forêt en 1 996; 5,5% ont été convertis
classification de l’image 201 4 est statistiquement en «savoka» et 1 ,3% en «champs de culture» en
acceptable. 2004. Ainsi, en 2004, la proportion restante de la
Cartographie d’occupation du sol par époque «forêt» était de 78,0%. Avec la dégradation
Les classifications des images Landsat ont permis continuelle des forêts, cette proportion n’a pas cessé
d’obtenir trois cartes d’occupation du sol datées de de diminuer et la forêt n’occupait plus que 67,9% de
1 996, 2004 et 201 4 (Figure 1 ) avec chacune trois la RSA;
classes. Concernant la classe «champs de culture», sa
Les proportions des classes d’occupation du sol du superficie reste plutôt stable (5,5% en 1 996, 4,2% en
paysage de la RSA pour les trois années d’études 2004 et 6,0% en 201 4). Par contre, si on examine
sont présentées dans le tableau 3. On en déduit que bien les tableaux 4 et 5, on peut voir que le taux de
le paysage de la RSA est dominé par la classe permanence de cette classe (sur les diagonaux, en
«forêt», mais avec une régression notable (d’environ gras) est faible. Autrement dit, la population exploite
82 à 68 %). Les 1 8% à 32% restants sont formés par les terrains seulement pendant quelques années et
les « savoka » et les «champs de culture». Les tests change de terrains à cultiver au fur et à mesure du
χ2 indiquent que les distributions de la proportion des temps.
différents types d’occupation du sol de la RSA ne sont
31 7
 TROPICULTURA, 201 7, 35, 4, 31 2-324

(Rakotondrasoa, 201 5)
Photo 3 : Champs de culture à Ambinaniroa (culture de riz pluvial
sur «tavy» ou culture sur brûlis).

Tableau 2
Matrice de confusion de la classification de l’image de 201 4

Tableau 3
Proportions des classes d’occupation du sol du paysage de la Réserve Spéciale
d’Ambatovaky pour les trois années d’étude

(Les données sont issues des cartes élaborées à partir des images Landsat datées de 1 996, 2004 et 201 4)

31 8
 TROPICULTURA, 201 7, 35, 4, 31 2-324

(Obtenue après classification des images LANDSAT , projection WGS84)

Figure 2 : Cartes d’occupation du sol de la Réserve

Spéciale d’Ambatovaky en 1 996, 2004 et 201 4.

31 9
 TROPICULTURA, 201 7, 35, 4, 31 2-324

(Obtenue sur base des cartes d’occupation du sol issues des classifications des images
LANDSAT datées de 1 996, 2004 et 201 4 de la Réserve Spéciale d’Ambatovaky)
Figure 3 : Illustration des changements observés au niveau de la
Réserve Spéciale d’Ambatovaky entre 1 996 et 201 4.

Tableau 4
Matrice de transition illustrant, en pourcentage (%) la
surface des types d’occupation du sol du paysage de
la RSA [données de 1 996 (lignes) et 2004
(colonnes)].

Tableau 5
Matrice de transition illustrant, en pourcentage (%) la
surface des types d’occupation du sol du paysage de
la RSA [données de 2004 (lignes) et 201 4
(colonnes)].

320
 TROPICULTURA, 201 7, 35, 4, 31 2-324

(Obtenu à partir de l'utilisation de la Chaine de Markov sur base des cartes d’occupation du sol datées de
1 996, 2004 et 201 4 de la Réserve Spéciale d’Ambatovaky)
Figure 4: Dynamique de l’occupation du sol de la Réserve Spéciale
d’Ambatovaky de 1 996 à 2050.

(Rakotondrasoa, 201 2)

Photo 4: «Tavy» à Ambavala (destruction de la forêt primaire par

abattage suivi de brûlis en vue de la culture du sol).

321
 TROPICULTURA, 201 7, 35, 4, 31 2-324

Les nouveaux terrains de culture sont issus à la fois Discussion

du défrichement de la forêt et des «savoka», comme
indiquent les tableaux 4 et 5 à travers les proportions Sur le plan méthodologique, l’étude cartographique et
de ces deux types d’occupation du sol converties en la production des cartes à partir du traitement
champs de culture pour les deux périodes d’étude. d’images satellites par télédétection ont permis de
En conséquence, pour les deux périodes, une détecter l’évolution spatiale et temporelle de
augmentation considérable de la superficie des l’occupation du sol. La validation de la classification
«savoka» a été notée. En 1 996, le «savoka» occupait effectuée sur ordinateur grâce aux relevés de points
1 2,2% de la superficie de la RSA; cette proportion a GPS sur terrain a permis de confirmer les cartes
augmenté et a atteint 1 7,8% en 2004 et 26,1 % en produites. En effet, la matrice de confusion a permis
201 4 (soit une augmentation de 45,9% entre 1 996 et de calculer la précision globale de la classification et
2004 et de 46,6% entre 2004 et 201 4). Cette l’indice Kappa, qui sont respectivement 84,1 % et
augmentation est liée à la conquête périodique de 78,9%. L’élaboration des cartes multidates a ressorti
nouveau terrain de culture par la population locale. où ont eu lieu les changements d’occupation du sol
En effet, les terrains de culture laissés évoluent en de la RSA et quelle a été la superficie d’une classe
«savoka», qui sont des formations secondaires convertie en une autre classe.
composées par des espèces à croissance rapide. Et Par ailleurs, la chaîne Markov a été choisie pour
cette conquête de terrain est une des principales effectuer la modélisation de la dynamique du paysage
causes de la déforestation; la population défriche la forestier de la RSA. En effet, pour maîtriser l’évolution
forêt pour la convertir en champs de culture. de ces processus de transformation du paysage, il est
La simulation de l’évolution future de l’occupation du important de s’appuyer sur des outils d’études
sol de 1 996 à l’horizon 2050 est présentée dans la cinématiques ou prospectifs permettant l’anticipation
figure 4. Les tests de comparaison de la composition de l’évolution du risque et parmi ces méthodes, la
paysagère simulée avec celle observée sur les cartes modélisation apparaît comme un réel outil d’aide à la
ont affichées les valeurs suivantes χ2 = 29,6; p > 0,05 décision (3). La chaîne de Markov a permis
et χ2 = 1 7,2; p > 0,05 respectivement pour l’année de d’appréhender les changements de catégorie au
vérification 2004 et pour l’année de validation 201 4. niveau des classes d’occupation du sol pour les deux
Autrement dit, il n’y a pas de différence significative périodes d’études 1 996-2004 et 2004-201 4. Elle a
entre la proportion des types d’occupations du sol l’avantage de fournir une méthode de complexité
issue de la simulation et ceux des cartes élaborées limitée par laquelle un système dynamique pourrait
pour les deux années 2004 et 201 4. En être détaillé et examiné (3). En plus de mesurer les
conséquence, la simulation obtenue est vérifiée, changements d’utilisation du territoire qui se sont
validée et jugée recevable. produits, cette approche permet de faire ressortir les
La simulation obtenue montre que la forêt et le tendances qui étaient non-évidentes dans le jeu de
«savoka» occuperont presque la même superficie en données de départ (23).
2050. En effet, la couverture forestière diminuera Pour la première période d’étude, c’est-à-dire entre
progressivement. Elle occupait 82,3% du paysage 1 996 et 2004, la couverture forestière de la RSA a
total en 1 996 et serait de 47,4% en 2050. La régressée avec un taux annuel de déforestation égal
réduction de cette classe s’effectuerait principalement à 0,6%. Ce taux est bas comparé à celui que
en faveur de l’augmentation de la superficie des Zaehringer a ressorti pour toute la forêt humide du
«savoka», qui passeront de 1 2,3% en 1 996 à 43,6% Nord-Est de Madagascar qui est de 1 % pour la
en 2050. Concernant les champs de culture, leurs période de 1 995-2005 (32). Cependant, pour presque
superficies vont connaître une augmentation: 5,5% les mêmes périodes, ce taux est élevé par rapport au
en 1 996 et 9,0% en 2050. Il est à signaler que ces taux moyen annuel de déforestation des aires
tendances partent de l’hypothèse selon laquelle les protégées malgaches qui est de 0,3% entre 1 990 et
facteurs potentiels et les variables qui influenceront la 2000 (24). En outre, entre les années 2004 et 201 4,
dynamique du paysage garderont la même emprise la déforestation de la RSA a eu plus d’ampleur, une
que celle de la période de 1 996 à 201 4; cette perte annuelle de forêt d’environ 1 ,3% a été
hypothèse correspond au principe de la «stationarité» enregistrée. Ce taux rejoint celui trouvé par Allnutt et
(30). al. pour le parc national Masoala entre 201 0 et 2011
qui est de 1 ,27% (1 ). Cependant, pour d’autres forêts
de l’Est, les taux de déforestations restent stables,
c’est le cas du corridor Manompana avec un taux de
1 ,07% entre 1 991 et 2004 et 1 ,09 entre 2004 et 2009
(11 ). Ainsi, on peut conclure que la couverture
forestière de la RSA a diminuée malgré son statut

322
 TROPICULTURA, 201 7, 35, 4, 31 2-324

d’aire protégée. Les actions de conservation ne se de la production agricole ne peut plus suivre les
limitent pas à la seule protection des ressources évolutions démographiques du milieu rural (20). En
naturelles, mais elles sont inséparables des actions effet, les populations dictées par leurs besoins accrus
de développement économique et social. Une sont contraintes à conquérir de nouveaux terrains
croissance démographique non accompagnée d’une pour l’agriculture de subsistance dans certaines
croissance économique pourrait être fatale à la régions à Madagascar (1 0).
pérennisation des ressources naturelles, car, dans ce Dans le cas d’Ambatovaky, la population riveraine
contexte, la population devra avoir recours aux défriche la forêt pour la convertir en champs de
ressources naturelles. En effet, la zone d’étude culture qu’elle abandonne après quelques années
présente un taux de croissance démographique pour défricher une autre partie de la forêt. Cette
annuel élevé d’environ 2,8% entre 2011 et 201 4 (35). pratique tend à réduire annuellement la superficie des
Cet accroissement démographique engendre de différents types de forêts humides à Madagascar. En
sérieux problèmes d’insuffisance de parcelles effet, la déforestation des forêts primaires du versant
culturales et de migration. Cette difficulté est Est peut en grande partie être attribuée à la culture
ressentie progressivement par les familles paysannes sur brûlis (20)
qui exploitent les mêmes parcelles dont elles ont
héritées de leurs parents et grand-parents. La Conclusions
migration se manifeste par le déplacement vers les Dans cette étude, nous avons déterminé la
régions plus fertiles et surtout à la périphérie des dynamique spatio-temporelle du paysage forestier de
forêts primaires. la Réserve Spéciale d’Ambatovaky. Des
Les tendances obtenues suite à la simulation à partir classifications des images Landsat datées de 1 996,
de la chaîne de Markov de premier ordre, ont révélé 2004 et de 201 4 ont permis de produire des cartes
une diminution de la superficie forestière de la RSA. d’occupation du sol de la RSA. Ces cartes ont
Ces tendances partent de l’hypothèse selon laquelle constitué une base pour la modélisation de l’évolution
les facteurs potentiels et les variables qui de la dynamique du paysage forestier de la RSA. Il en
influenceront la dynamique du paysage garderont la ressort que la couverture forestière de la zone d’étude
même emprise que celle de la période de 1 996 à régressera au fur et à mesure de temps malgré son
201 4. Il en découlerait que la couverture forestière va statut d’aire protégée. La perte de la forêt se
diminuer au profit principalement des «savoka», des manifestera au profit des champs de culture, qui sera
formations secondaires post-culturales (1 ). laissés après quelques années d’exploitation et se
L’augmentation de la superficie des «savoka» convertira en « savoka ». Cette situation peut être en
découle du développement et de l’extension du relation avec la pression démographique élevée
«tavy» (photo 4). Notre hypothèsee stipulant que «les entraînant une insuffisance des terres agricoles
pratiques culturales, notamment l’agriculture disponibles. Notre recherche a mis en évidence que
itinérante sur brûlis ou «tavy», favorisent la les actions de conservation ne doivent pas se limiter à
régression du couvert forestier du paysage la seule protection des ressources naturelles, mais
d’Ambatovaky, malgré son statut d’aire protégée» est qu’elles sont inséparables des actions de
ainsi acceptée. La limitation spatiale incite les développement économique et social de la population
paysans à cultiver le riz sur «tavy». Le fond du riveraine.
problème réside donc dans le fait que l'intensification

Références bibliographiques
1. Allnutt, T. F., Asner G. P., Golden C. D. & Powell G. V. N., 5. Bell E. J., 1 974, Markov analysis of land use change: an
201 3, Mapping recent deforestation and forest disturbance in application of stochastic processes to remotely sensed data,
northeastern Madagascar, Trop. Conserv. Sci., 6, 1 -1 5. Socio-Econ. Plan. Sci., 8 , 311 -31 6.
2. Aubert S., Razafiarison S. & Bertrand A., 2003, Déforestation 6. Blanc-Pamard C. & Rakoto Ramiarantsoa H., 2003,
et systèmes agraires à Madagascar. Les dynamiques des Madagascar: Les enjeux environnementaux. In: Lesourd M.,
tavy sur la côte orientale. CIRAD, CITE, FOFIFA, Montpellier, L’Afrique. Vulnérabilité et défis. Nantes: Editions du temps,
Antananarivo. 354-376.
3. Barima Y.S.S., Egnankou M.W., N’doumé C.T.A., 7. Bodin B., Ravilious C., Mant R. & Bastianelli C., 201 4, Les
Kouamé F.N. & Bogaert J., 201 0, Modélisation de la synergies entre la REDD+ et les objectifs d’Aichi de la
dynamique du paysage forestier dans la région de transition Convention sur la Diversité Biologique en Afrique Centrale -
forêt-savane à l’est de la côte d’Ivoire, Rev. Télédétect., L’apport de l’analyse spatiale pour la planification conjointe
9 , 1 29-1 38. de deux engagements internationaux sur les forêts. UNEP-
4. Barrett C.B., 1 999, Stochastic food prices and slash-and- WCMC, Cambridge, UK.
burn agriculture, Environ. Dev. Econ., 4, 1 61 -1 76. 323
 TROPICULTURA, 201 7, 35, 4, 31 2-324

8. Caloz R. & Collet C., 2001 , Traitements numériques d'images 22. Moore A.D., 1 990, The Semi-Markov process - a useful tool
de télédétection. Quebec, Canada: Presses de l'Université du in the analysis of vegetation dynamics for management,
Québec. J. Environ. Manage. 30 , 111 -1 30.
9. Chavez J.P.S., 1 988, An Improved Dark-Object Subtraction 23. Muller M.R. & Middleton J., 1 994, A Markov model of land-
Technique for Atmospheric Scattering Correction of use change dynamics in the Niagara region, Ontario, Canada,
Multispectral Data , Remote Sens. Environ., 24, 459 - 479. Landscape Ecol., 9 , 1 51 -1 57.
1 0. Clark M., 201 2, Deforestation in Madagascar: Consequences 24. ONE, DGF, FTM, MNP & CI, 201 3, Evolution de la couverture
of Population Growth and Unsustainable Agricultural des forêts naturelles à Madagascar 2005-2010. Antananarivo,
Processes, Global Majority E-Journal, 3 , 61 -71 . Madagascar.
11 . Eckert S., Ratsimba H. ., Rakotondrasoa L.O., 25. Pontius R.G. Jr., 2000, Quantification error versus location in
Rajoelison L.G. & Ehrensperger A., 2011 , Deforestation and comparison of categorical maps, Photogram. Eng, Remote
forest degradation monitoring and assessment of biomass Sens., 66, 1 011 -1 01 6.
and carbon stock of lowland rainforest in the Analanjirofo 26. République de Madagascar, 1 990, Loi n° 90-033 du 21
region, Madagascar, For. Ecol. Manage. , 262 , 1 996–2007. décembre 1990 portant Charte de l'Environnement malagasy,
1 2. FAO, 201 2, Situation des forêts 2012. Food and Agriculture Madagascar.
Organization of the United Nations, Rome. 27. Schlaepfer R., 2002, Analyse de la dynamique du paysage.
1 3. FAO, 201 3, FAO Statistical Yearbook 201 3. Food and Fiche d’enseignement 4.2, Laboratoire de Gestion des
Agriculture Organization of the United Nations, Rome. Ecosystèmes, Ecole Polytechnique de Lausanne, Suisse.
1 4. Gautier L. & Goodman S.M., 2003, Introduction to the flora of 28. Styger E, Rakotondramasy H.M., Pfeffer M.J., Fernandes E.
Madagascar .In The Natural History of Madagascar, ed. S.M. C.M. & Bates D.M., 2007, Influence of slash-and-burn
Goodman & J.P. Benstead, 229-250. Chicago, IL, USA: farming practices on fallow succession and land degradation
University of Chicago Press. in the rainforest region of Madagascar, Agric.,
1 5. Groupe des Spécialistes des Plantes de Madagascar, 2011 , Ecosyst. & Environ., 11 9 , 257–269.
Liste rouge des plantes vasculaires endémiques de 29. Toyi M.S., Barima Y.S.S., Mama A., André M., Bastin J -F.,
Madagascar. De Cannière C., Sinsin B. & Bogaert J., 201 3. Tree Plantation
1 6. Harper G.J., Steininger M.K., Tucker C.J., Juhn D. & Will Not Compensate Natural Woody Vegetation Cover Loss
Hawkins F., 2007, Fifty years of deforestation and forest in the Atlantic Department of Southern Benin, Tropicultura,
fragmentation in Madagascar, Environ. Conserv., 34, 31 , 62-70.
325–333. 30. Urban D.L. & Wallin D.O., 2002, Models of landscape
1 7. Kiener A., 1 963, Le “tavy” à Madagascar : ses différentes change. In: Learning landscape ecology: a practical guide to
forms et denominations, bilan du tavy et problèmes humains, concepts and techniques, Eds S.E. Gergel & T.M.G,
moyens de lutte, Bois et For. Trop., 90 , 9-1 6. Springer, New York, pp. 35-48.
1 8. Madagascar National Parks, 201 0, Plan de sauvegarde 31 . Weng Q., 2002, Land use change analysis in the Zhujiang
sociale et environnementale de la réserve spéciale Delta of China using satellite remote sensing, GIS and
Ambatovaky. stochastic modelling, J. Environ. Manage., 64, 273-284.
1 9. McConnell, 2002, Madagascar: Emerald isle or paradise 32. Zaehringer J.G., Eckert S. & Peter M., 201 5,
lost?, Environ. , 44, 1 0-1 4. Revealing Regional Deforestation Dynamics in North-Eastern
20. Messerli, P. 2003, Alternatives à la culture sur brûlis sur la Madagascar—Insights from Multi-Temporal Land Cover
Falaise Est de Madagascar: Stratégies en vue d'une gestion Change Analysis, Land, 4, 454-474.
plus durable des terres, African Studies Series, A1 7, 33. http://landsat.usgs.gov/Landsat_Processing_Details.php
Geographica Bernensia, Berne. 34. http://www.indexmundi.com/fr/madagascar/croissance_
21 . Minten B. & Moser C., 2003, F orêts: usages et menaces sur demographique_%28taux_de_croissance%29.html
une ressource. In : Agriculture, pauvreté rurale et politiques 35. https://www.instat.mg/madagascar-en-chiffre/
économiques à Madagascar, ed. B. Minten, J.C.
Randrianarisoa & L. Randrianarison, 86-89, Cornell
University, Ithaca.

L.O. Rakotondrasoa, Malgache, PhD, Enseignant chercheur, Université d’Antananarivo, Ecole Supérieure des Sciences Agronomiques, Département des Eaux
et Forêts, Antananarivo, Madagascar.

F. Malaisse, Belge, PhD, Professeur honoraire, Université de Liège, Gembloux Agro-Bio Tech, Gembloux, Belgique; Membre ARSOM.

J. Bogaert, Belge, PhD, Professeur ordinaire, Université de Liège, Gembloux Agro-Bio Tech, Gembloux, Belgique; Professeur ERAIFT, Knshasa, R.D. Congo;
Membre ARSOM.

324