Memoire Online - Structure des peuplements de la subéraie tunisienne:situation actuelle et devenir d'un écosystème.

Structure des peuplements de la subéraie tunisienne : situation actuelle et devenir d'un écosystème.

TABLE DE MATIERES

TABLES DE MATIERES..............................................................................I

RESUME...................................................................................................V

SUMMARY .............................................................................................VI

ÇáÎáÇÕÉ ....................................................................................................VII

LISTES DES TABLEAUX.........................................................................VIII

LISTES DES FIGURES................................................................................X

LISTE DES ABREVIATIONS .....................................................................XII

Introduction................................................................................................1

Chapitre I : Données bibliographique ...............................................................3

1.1. Structure du peuplement.............................................................................3

1.1.1. Etude du peuplement...............................................................................3

1.1.2. Etude de la structure...............................................................................3

1.1.2.1. La structure élémentaire.......................................................4

1.1.2.2. La composition floristique.....................................................6

1.1.2.3. La densité........................................................................6

1.1.2.4. Structures verticales, horizontales et diamétriques des tiges..............6

1.1.3. Régime et stade de développement du peuplement........................................7

1.2. Devenir d'un écosystème.......................................................................9

1.2.1. Evolution de la structure du peuplement forestier........................................9

1.2.2. Facteurs d'évolution de la structure du peuplement forestier........................10

1.2.2.1. Facteurs topographiques......................................................10

1.2.2.2. Le Dépérissement.............................................................11

1.2.2.3. Le Surpâturage.................................................................11

1.2.2.4. La Pression démographique............................................. ...12

1.3. Données générales sur le chêne-liège........................................................13

1.3.1. Répartition dans le monde......................................................................13

1.3.2. Répartition en Tunisie...........................................................................14

1.3.3. Sylviculture du chêne-liège.....................................................................16

Chapitre 2 : Le cadre naturel de l'étude.........................................................17

2.1. Présentation de la région.......................................................................17

2.2. Climat.............................................................................................17

2.2.1. Pluviométrie....................................................................................... 17

2.2.2. Température........................................................................................18

2.2.3. Vents..................................................................................................19

2.2.4. Neige.................................................................................................19

2.2.5. Humidité relative de l'air.......................................................................19

2.2.6. Synthèse bioclimatique..........................................................................19

2.3. Formations géologiques........................................................................21

2.4. Végétation........................................................................................21

Chapitre 3 : Matériel et Méthodes...............................................................23

3.1. Matériel végétal................................................................................. 23

3.1.1. Exigences écologiques du chêne-liège.................................................... ..23

3.1.2. Caractéristiques botaniques...................................................................24

3.2. Méthodologie....................................................................................24

3.2.1. Cartographie.......................................................................................24

3.2.2. Inventaire de l'état actuel de la subéraie en Kroumirie................................25

3.2.2.1. Pré-échantillonnage et calcul statistique...................................25

3.2.2.2 Echantillonnage................................................................27

3.2.2.3. Description des stations......................................................27

3.2.2.4. Description de la structure actuelle de la subéraie en Kroumirie.........30

3.2.2.5. Effet des facteurs du milieu sur la densité et la structure diamétrique du chêne-liège en Kroumirie..............................................................30

3.2.3. Dynamique de l'écosystème chêne-liège en Kroumirie ............... ................31

3.2.3.1. Dynamique de la subéraie en Kroumirie en fonction des facteurs du milieu ....................................................................................31

3.3. Analyses statistiques...........................................................................32

Chapitre 4 : Résultats..............................................................................33

4.1. Constat de la structure diamétrique actuelle de la subéraie en Kroumirie..............33

4.1.1. Description de la structure diamétrique actuelle .......................................33

4.1.2. Impact de certains facteurs du milieu sur la structure diamétrique actuelle.......................................................................................................35

4.1.2.1. Effet de l'altitude sur la densité et la structure diamétrique

du chêne-liège...........................................................................35

4.1.2.2. Effet de l'exposition sur la densité et la structure diamétrique

du chêne-liège...........................................................................37

4.1.2.3. Effet de la pente sur la densité et la structure diamétrique

du chêne-liège...........................................................................39

4.1.2.4. Effet du facteur `distances du peuplement aux douars les plus proches sur la densité et la structure diamétrique du chêne-liège...........................41

4.1.3. Variation de la densité du chêne-liège, toutes classes de diamètres confondues en fonction de certains facteurs du milieu pour chaque station..............................43

4.1.3.1. Répartition de la densité du chêne-liège, toutes classes de diamètres confondues, en fonction d'altitude...................................................43

4.1.3.2. Répartition de la densité du chêne-liège en fonction de l'exposition toutes classes de diamètres confondues, et par frange d'altitude.................44

4.1.3.3. Répartition de la densité du chêne-liège en fonction de la pente, toutes

d'exposition .............................................................................48

4.1.3.4. Répartition de la densité du chêne-liège, toutes classes de diamètres confondues, en fonction des distances aux douars les plus proches..............53

4.2. Dynamique de l'écosystème du chêne-liège en Kroumirie................................ 55

4.2.1. Dynamique de la subéraie de la Kroumirie en fonction des facteurs

du milieu ........................................................................................................55

4.2.1.1. Evolution de la densité et de la structure diamétrique de la subéraie

en fonction de l'altitude.............................................................55

en fonction de l'exposition...........................................................58

en fonction de la pente...............................................................60

4.2.1.4. Evolution de la densité et de la structure diamétrique de la subéraie

en fonction des distances aux douars les plus proches...........................62

4..2. Etude comparative entre la structure actuelle et les structures antérieures........65

CHAPITRE 5 : Discussions.......................................................................67

Conclusion générale................................................................................73

Références bibliographiques ..........................................................................74

Annexes...............................................................................................82

RESUME

Résumé : Structure des peuplements de la subéraie tunisienne : situation actuelle et devenir d'un écosystème.

Dans le présent travail nous avons essayé d'inventorier l'état actuel de la subéraie en Kroumirie, représenter cette structure actuelle et la comparer avec celles étudiées par les autres inventaires. De même, nous avons étudié l'impact des facteurs du milieu sur la densité et la structure diamétrique du chêne-liège. Enfin, nous avons évalué l'évolution de cette structure durant la période (1983 - 2003) en fonction des facteurs de perturbations.

Les résultats du dépouillement de la fiche placette (239 placettes inventoriées) ont montré que la structure des peuplements de la subéraie en Kroumirie est essentiellement représentée par les classes de semis et de semenciers. Les classes intermédiaires sont moins représentées. L'étude comparative a révélé que la structure actuelle est comparable à celle qui prévalait présentée en 1983.

Les analyses de variance de l'impact des facteurs du milieu ont montré que l'altitude, l'exposition, la pente et la distance par rapport aux douars les plus proches ont un effet sur la densité et la structure diamétrique du chêne-liège en Kroumirie. Les peuplements les plus équilibrés sont situés dans les stations alticoles, à pentes faibles ou modérés, d'expositions plus fraîches (Nord) et les plus éloignées des douars.

Les analyses de variance de l'effet des facteurs du milieu sur la répartition de la densité de chêne-liège, toutes classes de diamètres confondues, pour chaque station ont montré une grande variation au sein d'une même classe ou d'une classe à une autre d'altitude, d'exposition, de pente, des distances par rapport aux douars les plus proches. En effet, les fortes densités sont enregistrées dans les stations alticoles, à pentes faibles, d'expositions (Nord) et plus éloignées des douars.

L'étude de la dynamique de la subéraie en Kroumirie en fonction de ces facteurs durant la période (1983-2008) a été réalisée par la méthode d'inventaires périodiques. Les résultats de la comparaison ont montré que la régression de la subéraie a touché tous les stades des peuplements.

SUMMARY

Summary: Structure of the stands of Oak Tunisia: current situation and future of an ecosystem.

In this work we tried to identify the current state of Oak in Kroumirie represents the current structure and compare with those studied by other inventories. Similarly, we studied the impact of environmental factors on the density and diameter structure of cork oak. Finally, we evaluated the evolution of this structure during the period (1983 - 2003) in terms of disruption.

The results of the examination of the record plot (239 plots inventoried) have the structure that stands in the Oak in Kroumirie is mainly represented by the classes of seedlings and seed. The intermediate classes are less represented. The comparative study revealed that the current structure is comparable to that prevailing presented in 1983.

Analysis of variance of the impact of environmental factors showed that the altitude, exposure, slope and distance to nearest douars have an effect on the density and structure diameter of the cork oak in Kroumirie. The settlements are the most balanced in alticola stations at low or moderate slopes, exhibitions cooler (northern) and the most remote douars.

Analysis of variance of the effect of environmental factors on the distribution of the density of cork oak, all diameter classes combined, for each station showed a wide variation within the same class or a class another altitude, exposure, slope, relative distances to the nearest douars. Indeed, high densities were recorded at stations alticola, low slope, exposure (North) and more distant from the douars.

The study of the dynamics of Oak in Kroumirie according to these factors during the period (1983-2008) was carried out by the of periodic inventories. The results of the comparison showed that the regression of has affected all stages of growth.

ÇáãáÎÕ :ÚÑíÞÉ åíßáÉ ÇáãÏÑÌÇÊ ÛÇÈÇÊ Çááíä ÇáÊæäÓí : ÇáÍÇáÉ ÇáÑÇåäÉ æãÓÊÞÈá ÇáäÙÇã ÇáÅíßæáæÌí.

ÈÛíÉ ÇáÊÚÑ Úáì ÇáæÖÚíÉ ÇáÍÇáíÉ áåíßáÉ ÛÇÈÇÊ Çááíä í ÎãíÑ æ ÊÔÎíÕåÇ æãÞÇÑäÊåÇ ãÚ Êáß ÇáÊí ÏÑÓÊåÇ ÌÑÏ ÎÑì æ ÏÑÇÓÉ ÊËíÑ ÇáÚæÇãá ÇáÈíÆíÉ Úáì ßËÇÉ æ åíßáÉ ÈáæØ Çááíäí æ ÊÞíã ÊØæÑ åÐå ÇáåíßáÉ ÎáÇá ÇáÊÑÉ (1983 - 2003) ÊÍÊ ÊËíÑ åÐå ÇáÚæÇãá ÞãäÇ ÈÅÚÏÇÏ ÐÇ ÇáÚãá.

ËÈÊÊ äÊÇÆÌ ÍÕ 239 ÞØÚÉ Êã ÌÑÏåÇ ä åíßáÉ ÛÇÈÇÊ Çááíä í ÎãíÑ ÊãËáå ÓÇÓÇ ÆÇÊ ÇáÔÊáÇÊ æÇáÈÐæÑ. ßãÇ ßÔÊ ÚãáíÉ ÇáãÞÇÑäÉ ä ÇáåíßáÉ ÇáÍÇáíÉ åí ããÇËáÉ áÊáß ÇáÓÇÆÏÉ í ÚÇã 1983.

ÙåÑÊ ÇáÊÍÇáíá ä ÇáÇÑÊÇÚ æÇáÊÚÑÖ æÇáãäÍÏÑ æÇáãÓÇÉ Åáì ÞÑÈ ÏæÇÑ áåÇ ÊËíÑ Úáì ßËÇÉ æÇáåíßáÉ ÇáÞØÑíÉ áÕíáÉ ÇáÈáæØ Çááíäí ÈÎãíÑ.

äáÇÍÙ ä ÇáãÓÊæØäÇÊ ÇáßËÑ ÊæÇÒäÇ åí ÇáãæÌæÏÉ í ÇáãÑÊÚÇÊ æ í ãäÍÏÑÇÊ ÇáãäÎÖÉ æ ÇáãÊæÓØÉ æÇáãÚÇÑÖÉ ááÈÑæÏÉ (ÔãÇá) æÇáÈÚíÏÉ Úä ÇáÊÌãÚÇÊ ÇáÓßäíÉ (ÇáÏæÇÑ).

ßãÇ ÙåÑÊ ÇáÊÍÇáíá ä ÇáÚæÇãá ÇáÈíÆíÉ áåÇ ÊËíÑ Úáì ÊæÒíÚ ßËÇÉ ÇáÈáæØ Çááíäí í ßá ãÍØÉ (ßá ÇáØÈÞÇÊ ÇáÞØÑíÉ ãÚÇ) ÅÐ äáÇÍÙ ä åäÇáß ÊÇæÊÇ æÇÓÚÇ í ÏÇÎá äÓ ÏÑÌÉ ÇÑÊÇÚ æÇáÊÚÑÖ æÇáÇäÍÏÇÑ æ ÞÑÈ ÇáãÓÇÉ ãä ßá ÏæÇÑ. áÞÏ Êã ÊÓÌíá ßËÇÉ ÚÇáíÉ í ÇáãÍØÇÊ ÇáÚÇáíÉ ÇáÇÑÊÇÚ æÇáãäÎÖÉ ÇáÇäÍÏÇÑ æÇáãÚÇÑÖÉ ááÈÑæÏÉ (ÔãÇá) æÇáÈÚíÏÉ Úä ÇáÏæÇÑ.

ÊãÊ ÏÑÇÓÉ ÏíäÇãíÇÊ Çááíä í ÎãíÑ æÞÇ áåÐå ÇáÚæÇãá ÎáÇá ÇáÊÑÉ (1983-2008) ÈØÑíÞÉ ÓáæÈ ÇáÌÑÏ ÇáÏæÑí. äÊÇÆÌ ÇáãÞÇÑäÉ ÙåÑÊ ä åÐÇ ÇáÊÑÇÌÚ ËÑ Úáì ÌãíÚ ãÑÇÍá Çáäãæ ÇáÈáæØ Çááíäí.

ßáãÇÊ ãÇÊíÍ : åíßáÇááíä ÇáÊæäÓí ÎãíÑ æÇáÍÇÖÑ æÇáãÓÊÞÈá. LISTE DES TABLEAUX

Tableau 1 : L'évolution démographique entre 1956 et 1994, dans la zone de Kroumirie

et gouvernorat de Jendouba (I.N.S)....................................................................12

Tableau 2 : Répartition des subéraies selon les superficies........................................14

Tableau 3 : Régression de la superficie de chêne-liège en Tunisie (1950-2005)............... 14

Tableau 4 : Répartition de la superficie de chêne-liège par gouvernorat........................15

et Tabarka (1900-2005) .................................................................................17

Tableau 6. Répartition saisonnière des précipitations des pluies (1900-2005)....................18

et Ain Draham (1900-2005).............................................................................18

Tableau 8 : Tableau récapitulatif des placettes de pré échantillonnage.............................26

Tableau 9 : Tableau récapitulatif du calcul statistique de pré échantillonnage..................27

Tableau 10 : Caractéristiques écologiques des stations.............................................27

Tableau 11 : Classes de diamètres retenus pour présenter la structure de la subéraie.........30

Tableau 12 : Les différentes variables des analyses statistiques..................................31

Tableau 13 : Résultats de l'analyse de variance des effets de l'altitude sur la répartition

des classes de diamètres.................................................................................36

Tableau 14 : Résultats de l'analyse de variance des effets de l'exposition sur la répartition des classes de diamètres.................................................................................38

Tableau 15 : Résultats de l'analyse de variance des effets de la pente sur la répartition

des classes de diamètres................................................................................40

Tableau 16 : Résultats de l'analyse de variance de l'effet de la nature du sol sur la répartition des classes de diamètres.................................................................................41

Tableau 17 : Résultats de l'analyse de variance des effets de l'altitude sur la densité du chêne-liège (toutes classes d'altitudes confondues)..................................................43

Tableau 18 : Résultats de l'analyse de variance des effets de l'exposition sur la densité du chêne-liège, à basse altitude et toutes classes de diamètres confondues .........................45

Tableau 19 : Résultats de l'analyse de variance des effets de l'exposition sur la densité du chêne-liège, à moyenne altitude (toutes classes de diamètres confondues).......................46

Tableau 20 : Résultats de l'analyse de variance des effets de l'exposition sur la densité du chêne-liège, à haute altitude (toutes classes de diamètres confondues) ...........................47

Tableau 21 : Résultats de l'analyse de variance des effets de la pente sur la densité

en exposition Nord........................................................................................48

Tableau 22 : Résultats de l'analyse de variance des effets de la pente sur la densité

en exposition Sud........................................................................................51

Tableau 23 : Résultats de l'analyse de variance des effets de la distance des douars les plus proches sur la densité du chêne-liège (toutes classes de diamètres confondues).....................53

Tableau 24 : Résultats de l'analyse de variance des effets de l'altitude sur la densité

et la structure diamétrique (en 1983)...................................................................56

Tableau 25 : Résultats de l'évolution de la densité et la structure diamétrique du chêne-liège en Kroumirie en fonction de l'altitude, durant la période (1983-2008)...........................56

Tableau 26 : Résultats de l'analyse de variance des effets de l'exposition sur la densité et la structure diamétrique (en 1983)........................................................................58

Tableau 27 : Résultats de l'évolution de la densité et la structure diamétrique du chêne-liège en Kroumirie en fonction de l'exposition durant la période (1983-2008)........................59

Tableau 28 : Résultats de l'analyse de variance des effets de la pente sur la densité

et la structure diamétrique (en 1983)..................................................................60

Tableau 29 : Résultats de l'évolution de la densité et la structure diamétrique du chêne-liège en Kroumirie en fonction des classes de pente durant la période (1983-2008)..................61

Tableau 30 : Résultats de l'analyse de variance des effets de la distance des peuplements les plus proches par rapport aux douars sur la répartition des classes de diamètres (en 1983)....62

Tableau 31 : Résultats de l'évolution de la densité et la structure diamétrique du chêne-liège en Kroumirie en fonction de la distance des douars les plus proches

durant la période (1983-2008).............................................................................................63

LISTE DES FIGURES

Figure 1 : Représentation de la structure élémentaire régulière....................................4

Figure 2 : Représentation de la structure élémentaire jardinée....................................5

Figure 3 : Représentation de la structure élémentaire irrégulière.................................5

Figure 4 : Répartition mondiale du chêne-liège......................................................13

Figure 5: Proportion des espèces feuillues en Tunisie..............................................15

de Tabarka (1900-2005).................................................................................20

Ain Draham (1900-2005)................................................................................20

Figure 8: Localisation des placettes d'échantillonnages............................................29

Figure 9: Représentation de la structure diamétrique actuelle pour chaque station

(en 2008)..................................................................................................34

Figure 10: Structure diamétrique globale de la subéraie en Kroumirie (en 2008)..............35

Figure 11 : Effet de l'altitude sur les densités (effectifs par hectare) et la structure

diamétrique de peuplements de chêne-liège au niveau de l'ensemble de la zone d'étude.......36

Figure 12 : Effet de l'exposition sur les densités (effectifs par hectare) et la structure diamétrique des peuplements de chêne-liège au niveau de l'ensemble de la zone d'étude.....38

Figure 13 : Effet de la pente sur les densités (effectifs par hectare) et la structure

diamétrique des peuplements de chêne-liège..........................................................40

Figure 14 : Effet de la distance des peuplements aux douars les plus proches sur les densités (effectifs par hectare) et la structure diamétrique des peuplements de chêne-liège..............42

Figure 15 : Effet de l'altitude sur la variation de la densité de chêne-liège, toutes

classes de diamètres confondues au niveau des stations d'études.................................44

Figure 16 : Effet de l'exposition sur la densité du chêne-liège, toutes classes de diamètres confondues au niveau des stations de basses altitudes............................................. 45

Figure 17 : Effet de l'exposition sur la densité du chêne-liège, toutes classes de diamètres confondues au niveau des stations de moyennes altitudes..........................................46

Figure 18 : Effet de l'exposition sur la densité du chêne-liège, toutes classes de diamètres confondues, au niveau des stations alticoles.........................................................47

Figure 19 : Effet de la pente sur la densité de chêne-liège, toutes classes de diamètres et d'altitudes confondues et en expositions Nord.................................................................50

Figure 20 : Effet de la pente sur la densité de chêne-liège, à toutes classes de diamètres et d'altitudes confondues et en exposition Sud..........................................................53

Figure 21 : Effet de la distance des douars les plus proches sur la densité

du chêne-liège, toutes classes de diamètres confondues, pour l'ensemble des stations d'étude..................................................................................................... 54

Figure 22 : Evolution de la densité et de la structure de chêne-liège en fonction de l'altitude

(, et C3)............................................................................................58

Figure 23 : Evolution de la densité et de la structure diamétrique du chêne-liège

en fonction de l'exposition..............................................................................60

Figure 24 : Evolution de la densité et de la structure diamétrique du chêne-liège

en fonction de la pente...................................................................................62

Figure 25 : Evolution de la densité et de la structure diamétrique du chêne-liège en fonction de la distance aux douars les plus proches.............................................................64

Figure 26 : Structure diamétrique de la subéraie en Kroumirie durant la période 1995........65

Figure 27 : Structure diamétrique de la subéraie en Kroumirie durant la période 2005........65

Figure 28: Structure diamétrique globale de la subéraie en Kroumirie (1983)...................66

Figure 29 : Structures diamétriques de la subéraie en Kroumirie durant les années

1983, 1995, 2005 et 2008...............................................................................66

LISTE DES ABREVIATIONS

J : Janvier ; F : Février, M : Mars ; Av : Avril ; Jn : Juin ; Jt : Juillet ; At : Août ; S : Septembre ; O : Octobre ; N : Novembre ; D : Décembre ;

*** : paramètre très hautement significatif (P < 0,001) ; ** : paramètre hautement significatif (P < 0,01) ; * : paramètre significatif (P < 0,05) ; ns : paramètre non significatif.

Dédicaces

A' tous ceux qui m'ont soutenu moralement pendant la rédaction de ce mémoire. Plus particulièrement à mon père, ma mère, Asma, Rayen ;

Remerciements

Je tiens à dresser mes vifs remerciements, à Monsieur Brahim Hasnaoui mon directeur de mémoire, professeur à l'Institut Sylvo - Pastoral de Tabarka pour son encadrement, ses conseils enrichissants, ses orientations pertinentes, et sa patience enthousiasmée.

Mes remerciements s'adressent aussi à Madame Samira SMITI, professeur à la faculté des sciences de Tunis, pour sa patience, sa gentillesse, sincère reconnaissance et qui nous a fait l'honneur de présider notre jury de mémoire.

Merci à Monsieur Mohamed Néjib REJEB, Directeur de recherche Général de l'Institut National de recherche en Génie Rural, Eaux et Forêts, pour sa rigueur scientifique et pour avoir accepté de juger ce travail.

Je tiens à exprimer ma profonde gratitude à Monsieur Toumi LAMJED, Maître de conférences et Directeur général de l'Institut Sylvo - Pastoral de Tabarka.

Je remercie également Monsieur Chaâbane ABBES, Maître assistant à l'Institut Sylvo - Pastoral de Tabarka, pour ces conseilles.

J'adresse ma reconnaissance à Monsieur Foued Hasnaoui, pour son aide généreuse, ses conseils, ses encouragements tout au long de travail et pour m'avoir aidé aux analyses statistiques de mes données.

Mes remerciements s'adressent également à Monsieur Naceur BOUSSAIDI, pour leur soutien moral.

Je tiens à remercier Monsieur Ali ALOUI, Maître assistant à l'Institut Sylvo - Pastoral de Tabarka pour ses remarques.

Merci aux enseignants de l'Institut Sylvo - Pastoral : Rabeh NSIBI, Foued ALOUI, Youssef SAIDI et Abdellah RJEIBI, pour leurs efforts sincères et leurs conseilles.

Je salue également l'équipe m'ayant accompagné lors de la récolte des données : Mounir Mlaouhia, Sarra SAIDI, Salwa ALOUI

Merci aussi à Jdidi Mouhamed, Fethi SELMI, Souhir MLEIKI, Latifa JDAIDI, Imen MLEIKI, Ikbel ZOUAOUI, Sameh MHAMDI, Hanen GHAZGAZI, Emna KRIMI.

Je vaudrais aussi remercier ma famille pour son soutien, sa patience, et sa compréhension dans les moments les plus difficiles de cette expédition. Merci à mon oncle Chérif.

Introduction

Les forêts couvrent environ 30 % de la surface du globe et contiennent une grande partie des espèces qui peuplent les écosystèmes terrestres. La diversité des espèces des forêts est positivement corrélée avec l'augmentation de la complexité des structures et l'hétérogénéité spatiale à diverses échelles écologiques (Douglas et al, 2007).

Selon Goreaud et al (2005), ces forêts sont des écosystèmes complexes, elles abritent de nombreuses espèces animales et végétales en interaction, et participent ainsi à maintenir une grande biodiversité. Ces forêts peuvent jouer également un rôle important dans les grands cycles biogéochimiques, notamment au niveau du climat.

Courbeau et al, (1993) ont montré que l'étude de la structure du peuplement forestier et leur évolution revêt une importance particulière. D'une part, ces structures ont présenté une grande diversité notamment, tout un ensemble de structures régulières ou irrégulières mal définies est à prévoir ; d'autre part, elles ont déterminé la capacité de la forêt à remplir ses multiples fonctions : production de bois, mais aussi rôle paysager et écologique en fonction de la protection humaine.

Sokopn et Biaou (2002) ont prouvé que la structure des peuplements et en particulier les distributions de diamètres de l'espèce est une variable importante qui doit être examinée dans le développement durable des systèmes de gestion et une base pour suivre leur développement après la récolte des ressources.

Naggar (2000) a suggéré que la forêt des pays du Maghreb joue un rôle stratégique tant sur les plans socio-économique et pastoral qu'environnemental. Elle a constitué un patrimoine par la diversité des systèmes écologiques et de biodiversité qu'elle a intégrés et par l'importance de son étude sur environ 13,5 millions d'hectares.

La forêt tunisienne constitue un patrimoine national et joue un rôle aussi bien sur le plan écologique, environnemental qu'économique et social, elle occupe une superficie totale de 654.220 ha, la superficie totale des espèces dans la région Nord-Ouest est de 368.363 ha soit 44 % de celle enregistrée à l'échelle nationale. Elle est caractérisée par 22 % de feuillus, 55 % de résineux et 23 % de maquis et de garigue (O.D.N.O, 2004).

La composition, la structure et la densité des peuplements forestiers en Tunisie varient avec les principaux étages bioclimatiques, les forêts les plus denses et les plus riches en composition floristique s'étalent dans la Kroumirie.

La forêt de chêne-liège occupe une place importante dans le patrimoine forestier tunisien, elle représente 3 % de la superficie subéricole du bassin méditerranéen occidental et de la côte

Elle présente une situation alarmante qui se traduit par une régression temporelle et spatiale continue dans la majorité de son aire de répartition.

Selon Hasnaoui (1998), la superficie de cette espèce a régressé de 40.000 ha en 40 ans (entre 1952 et 1992), soit 1.000 ha / an ou environ 1% / an. Actuellement, il ne reste que 69.870 ha dont 49.142 ha au niveau de la Kroumirie, soit 70 % de celle enregistrée à l'échelle nationale (D.G.F, 2005).

Cette forêt, soumise durant ces dernières années à un pâturage intense, défrichement remarquable pour gagner des terrains agricoles, des coupes abusives des maquis et même des grands arbres, ainsi qu'à une grande pression naturelle (sécheresse, chaleur estivale, facteurs climatiques et orographiques). Cette situation provoque des modifications progressives au niveau de la structure, de la diversité et de la densité des peuplements au niveau des écosystèmes de chêne-liège.

L'action combinée de ces facteurs provoque une élimination progressive des strates basses et des semis de chêne-liège et en amène à des difficultés majeures de la régénération naturelle et à des modifications des structures des peuplements (Hasnaoui, 1992).

Les travaux de recherche sur la structure et le devenir de l'écosystème de chêne-liège sont rares en Tunisie. En 1992, Hasnaoui a donné une idée sur l'état de la structure diamétrique et les difficultés qui rencontrent nos subéraies à se reconstituer naturellement par semis. En 2004, Hasnaoui et al. ont présenté, par deux méthodes (structure et architecture de la végétation), l'échec de la reconstitution naturelle de ces semis et la prévision de leur avenir.

Le présent travail a pour objectifs d'inventorier l'état actuel de la subéraie en Kroumirie et essayer de dégager les relations possibles entre densité et structure diamétrique en fonction des facteurs du milieu, et de mieux comprendre enfin la dynamique du chêne-liège en fonction de ces mêmes facteurs durant une période d'un quart de siècle qui s'étend de 1983 à 2008 et de comparer sa structure diamétrique actuelle avec des structures antérieures.

Au premier chapitre, nous avons passé en revue les données bibliographiques se rapportent à la structure des subéraies, leur devenir, leur répartition dans le monde et en Tunisie ainsi qu'aux techniques sylvicoles de cette essence.

Le troisième chapitre décrit la méthodologie suivie et le matériel utilisé.

Et enfin un cinquième chapitre réservé à la discussion et la conclusion générales.

Chapitre 1 : Données bibliographiques

1.1. Structure du peuplement

1.1.1. Etude d'un peuplement

Un peuplement est un ensemble d'arbres ayant une uniformité jugée suffisante quant à sa composition, sa structure, son âge, sa répartition, son état sanitaire, etc., pour se distinguer des peuplements voisins, et pouvant ainsi former une unité élémentaire sylvicole ou d'aménagement (D.G.F, 2005).

Un peuplement constitué par une seule essence est dit pur. Il peut être pur feuillu (chêne-liège, chêne zeen, etc.) ou pur résineux (Pin d'Alep, Pin Maritime, etc.).

Un peuplement de chêne-liège est dit mélangé si la densité en chêne-liège est insuffisante ou si la dynamique naturelle des autres essences entrave les possibilités de régénération du chêne-liège (Gaudin., 1996).

Les peuplements réguliers sont composés d'arbres d'une même espèce et de même âge. Ils peuvent être issus de plantation ou de régénération naturelle menée sur l'ensemble d'une parcelle et constituant la méthode sylvicole de référence pour les essences sociales en plaine (Jarret, 2004).

Pour Turckheim et al. (2005), les peuplements irréguliers sont composés d'arbres d'âges et de tailles différents et ont une dynamique plus complexe que les peuplements purs réguliers et sont plus difficiles à gérer.

D'après Goureaud et al. (2005), ces peuplements très fréÿÿents dans les forêts de montagne, moins faciles d'accès et moins rentables, ont joué surtout un rôle de protection contre l'érosion. Dans ces peuplements, la régénération naturelle est exprimée de façon hétérogène, notamment en fonction de la lumière disponible au sol et entraîne une grande hétérogénéité des structures forestières.

1.1.2. Etude de la structure

La structure d'un peuplement, par définition, est une caractéristique complexe associant le régime et la distribution verticale. Elle est déterminée indépendamment du mode de traitement appliqué par l'administration forestière (D.G.F, 2005).

D'après Walter (1979), cette structure explique le mode d'organisation d'un écosystème ou de ses composants.

Bourlière et Lamote (1978) ont suggéré que la structure permet la connaissance, dans l'espace et dans le temps de la position des individus de toutes tailles et tous âges constituant la communauté, ce qui conduit à la compréhension des interactions entre ses différents composants et le contrôle du fonctionnement de l'écosystème.

D'après ces mêmes auteurs, cette structure peut être élémentaire ou globale. Elle est caractérisée par la diversité floristique, la densité, la surface terrière, les répartitions horizontale, verticale et diamétrique des tiges.

1.1.2.1. Structure élémentaire

Selon Franclet (1972), cette structure est représentée par une courbe en cloche (courbe de Gauss), elle caractérise la structure des peuplements équiennes (âges voisins et réguliers), au niveau de laquelle tous les arbres ont tous des âges voisins et ont donc des diamètres peu différents (figure 1).

Figure 1 : Répartition de la structure élémentaire régulière (Gaudin, 1996)

- Structure élémentaire jardinée Selon Gaudin (1996), cette structure est représentée sous forme d'une fonction sensiblement exponentielle, au niveau de laquelle tous les âges ou tous les diamètres sont représentés. On a une décroissance régulière du nombre de tiges quand on passe d'une classe de diamètre à la suivante (figure 2).

Figure 2 : Représentation de la structure élémentaire jardinée (Gaudin, 1996)

- Structure élémentaire irrégulière : Franclet (1972) a montré que cette structure n'est pas reconnue comme une structure naturelle, elle est en effet très artificielle et provient d'exploitation ou de conditions très particulières.

Letocart (2002) a révélé que la structure irrégulière offre de nombreux avantages d'ordre écologique et économique et elle a rendu les peuplements plus résistants aux divers dangers; elle a garanti une couverture permanente du sol et favorise la régénération.

Figure 3 : Représentation de la structure élémentaire irrégulière (Gaudin, 1996)

1.1.2.2. Composition floristique

La composition floristique a permis de décrire la forêt sous un aspect écologique et fonctionnel. Sa liaison éventuelle avec la dynamique du peuplement, la stabilité de l'écosystème et sa productivité, influe évidemment sur la dynamique moyenne du peuplement et elle peut être indicatrice de l'intensité de perturbations (Vincent et al, 1998). Le peuplement de chêne-liège a montré une composition beaucoup plus homogène du point de vue floristique, car il évite les stations trop humides, les endroits trop froids et les sols non acides (le chêne-liège est intimement lié aux formations gréseuses de l'Oligocène Numidien et au Trias non calcaire) (Saoudi, 1994).

Selon Amandier (2002), le feuillage du chêne-liège est assez léger, transparent et la subéraie, dans son état optimal, est un peuplement forestier relativement clair où les arbres ne couvrent que 60% du terrain. De même, Saoudi (1994) a signalé que, en altitude et le long des vallons frais, le chêne-liège est en contact avec le chêne zeen. Sur les dunes, le chêne-liège a beaucoup de peine à supplanter le chêne kermès.

1.1.2.3. La densité

La densité du peuplement est l'effectif par unité de surface, c'est une caractéristique intensive discrétisée en classe de densité de peuplement (D.G.F, 2005). Selon Franclet (1972), le nombre de tiges augmente pendant la période de construction du peuplement, puis il diminue très rapidement et d'une façon continue sous l'effet de plusieurs facteurs tels que l'influence de l'homme, la concurrence vitale entre les espèces et les sujets d'une même espèce et les facteurs accidentels comme la neige collante, les vents violents et les invasions d'insectes.

1.1.2.4. Structures verticale, horizontale et diamétrique des tiges

- Structure verticale : Elle exprime le type d'occupation de l'espace aérien par les végétaux et l'agencement des différentes strates constituées par des espèces différentes ou par une même espèce mais à différents stades de croissance (Hasnaoui, 1992).

Pour Collinet (1997), la structure verticale représente la distribution des individus par classes de hauteur, elle offre l'intérêt de pouvoir fournir un indicateur de richesse du site.

Selon M'Hirit (1982), la structure verticale de la végétation correspond au recouvrement du sol par les différentes strates arborescente, arbustive, sous arbustive et herbacée, elle peut renseigner sur les potentialités de la station et sur le dynamisme du peuplement forestier.

- Structure horizontale : Elle rend compte de la répartition des individus dans le plan. L'étude de la structure horizontale est importante car cette répartition résulte des mécanismes de fonctionnement des espèces (mode de régénération, sensibilité aux facteurs du milieu, compétition, etc.).

On distingue souvent l'étude des répartitions individu - individu et celle des répartitions individu - milieu (Dérout, 1994).

Dauglas et al (2007) ont signalé que la variabilité horizontale dans la structure du peuplement a été proposée comme un mécanisme qui contribue à la diversité des espèces dans les forêts anciennes.

- Structure diamétrique : Selon Vincent et al (1988) la structure diamétrique est la répartition des tiges par classes de diamètres, elle peut être établie en prenant en compte tous les individus et toutes espèces confondues. Elle est porteuse d'informations sur la stabilité du peuplement, elle peut être également établie par espèce et l'on parle de structure spécifique.

1.1.3. Régime et stade de développement du peuplement

On a distingué trois régimes caractérisés principalement par le mode de régénération.

Selon Bourdru (1989), le peuplement est un taillis et est soumis au régime du taillis lorsqu'il est rajeuni par des rejets de souche ou par des drageons.

Pour Schütz (1989), le peuplement est un taillis sous futaie et est soumis au régime du taillis sous futaie quand les arbres se reproduisant par graines surplombent d'autres végétaux rajeunis par rejets ou drageons. Il comprend, au-dessus, la futaie ou réserve et au-dessous, le taillis.

Selon Franclet (1972) le régime de la futaie est basé sur la régénération à l'aide de semences. Cette régénération se poursuit, sans discontinuité, pendant des générations dans les forêts cultivées et modifie la structure et la composition botanique des peuplements primitifs.

Goreaud et al. (2005) distinguent, quant à eux, Futaie régulière (équienne) et futaie jardinée (inéquienne). Une futaie est régulière lorsqu'elle est composée de peuplements réguliers élémentaires, d'âges convenablement gradués pour qu'elle puisse être parcourue à intervalles réguliers par des coupes principalement dans les peuplements les plus vieux. Une futaie est inéquienne lorsqu'elle est constituée d'arbres d'âges divers ou multiples, dont les dimensions en grosseur et en hauteur sont assez différentes pour qu'elles ne puissent se régulariser avec le temps.

Certains auteurs comme Boudru (1973), Schütz (1989) et Franclet (1972) ont classé la futaie jardinée en plusieurs classes :

La futaie est dite jardinée par pieds d'arbres et les surfaces élémentaires du jardinage par pied couvrent l'étendue de toute la forêt. Dans ce type de peuplement, les arbres de tous les âges et de toutes les dimensions sont intimement mélangés.

Graphiquement, les fréquences décroissantes par catégorie de grosseur sont représentées par une courbe à concavité tournée vers le haut.

La futaie est dite par groupes ou par bouquets, lorsqu'elle comprend des groupes étendus, variables, confusément mélangés.

Les âges et les dimensions d'arbres varient de groupe en groupe de telle sorte que, sur une surface suffisante on rencontre tous les âges ou toutes les dimensions. Chaque groupe est considéré comme un peuplement régulier.

La futaie jardinée par quartier est un type intermédiaire, par sa structure, son traitement et même son aménagement. La forêt se répartit en plusieurs quartiers, trois par exemple, respectivement occupés par les vieilles tiges, par les bois d'âge moyen et par les jeunes plantes. Ces peuplements sont d'âges multiples, chaque quartier est formé de groupes équiennes.

Il est préférable de délimiter les états de développement par les dimensions des arbres, plus précisément par leur diamètre. Certains chercheurs comme Bourdru (1973) et Schütz (1989) ont utilisé la notion de diamètre dominant (Do).

Les états de développement des peuplements se repartissent en classes de diamètre dominant de 10 cm.

Ces mêmes auteurs, ont distingué les stades de développements de taillis et de futaie suivants :

- Le jeune recru : Après la coupe, les rejets apparaissent sensiblement à la même époque. Les cépées de rejets sont séparées par des espaces qui se garnissent d'un tapis herbacé et de quelques bois - morts.

- Le jeune taillis : Les cépées se développent en hauteur, elles s'étalent en largeur, les branches s'entrecroisent.

- Le vieux taillis : Rappelle un peu du bas - perchis passant au haut - perchis de la futaie régulière.

- L'état de semis : il est constitué par de jeunes sujets issus de graines (hauteur inférieure à 40 cm). Il s'agit des tiges non ou peu ramifiées.

- L'état de fourré : il est formé par des tiges ligneuses issues de graines, ramifiées (40 ou 50 cm à1 à 3 m de haut).

- L'état gaulis : comprend de jeunes tiges de faibles hauteurs (3 à 8 m) et de diamètres (1 à 5 cm).

- Le bas - perchis : comprend les arbres de diamètre 10 à 20 cm (circonférence de 30 à 60 cm) et de hauteur supérieure à 8 m.

- Le haut - perchis : comprend les arbres de diamètre 20 à 30 cm (circonférence 60 à 90 cm) et de hauteur 15 m.

- L'état de futaie : caractérisé par des arbres de diamètres supérieurs à 30 cm et de hauteur entre 15 à 20 m.

1.2. Devenir d'un écosystème

1.2.1. Evolution de la structure d'un peuplement forestier

La structure de la forêt naturelle est considérée, sur une grande surface, comme une forêt jardinée avec un certain excès de gros bois. Elle peut varier d'allure de place en place mais reste généralement étagée. La régénération est due au hasard, l'élimination des arbres, le plus souvent les plus gros, se fait isolément (sénilité, foudre, champignons, insectes), par groupes plus ou moins étendus, ou encore par grands blocs (accidents météorologiques) (Boudru, 1989).

Ce même auteur suggérait que les semis et les arbres nés dans ces trouées (ou grands vides) forment des peuplements de structures très diverses depuis les groupes mélangés jardinés par pieds jusqu'aux parcelles équiennes mixtes ou pures.

Zukrigl et al (1963) ont montré que dans la forêt naturelle existe un cycle d'évolution de longue durée avec des stades différents et successifs (phases) répartis localement : l'ensemble de la forêt conserve toujours sa composition botanique générale et sa structure jardinée.

D'après Greaud et al (2005), pour comprendre la dynamique des peuplements forestiers et prédire leur évolution pour mieux les gérer, il est important de prendre en compte les processus de régénération à l'échelle du petit massif forestier, car ce sont eux qui sont à l'origine de la colonisation de ces espaces, et qui continuent à modifier la composition spécifique des peuplements, dans le cadre d'une succession écologique des espèces. Ces mêmes auteurs, ont montré que la croissance et la mortalité sont toutes deux fonctions de la compétition entre arbres (qu'ils soient ou non de la même espèce, et quels que soient leur âge et leurs dimensions).

Zukrigl et al (1963) ont montré également que le cycle d'évolution de la structure des peuplements forestiers présente des stades successifs qui sont :

? Une phase optimale (A) ou le peuplement fermé, régularisé, de hauteur uniformisée, très vigoureux, renferme peu d'arbres dépérissants et de rares jeunes bois. Il est constitué par les espèces de grandes tailles du milieu, d'autres étant ravalées dans les étages subordonnés.

? Une phase de vieillissement (B), ou de sénilité, à structure étagée avec arbres géants. La jeunesse naît timidement dans quelques trouées formées par disparition d'individus surâgés.

? Une phase de destruction (C) temporaire de la forêt très vieille. Les sujets disparus, morts ou dépérissants, sont nombreux. De multiples trouées se garnissent abondamment de recrus jeunes.

? Une phase jardinée initiale (D) ou phase de régénération avec jeunesse abondante et persistante des derniers vieux bois de l'ancienne forêt dont les cadavres jonchent le sol.

? Une jardinée typique (O) composée d'individus d'essences dominantes de toutes dimensions. Les espèces accessoires sont également présentes à tous les niveaux. Les chablis de vieillesse sont absents ou exceptionnels.

1.2.2. Facteurs d'évolution de la structure du peuplement

De nombreux facteurs, tels que la productivité des semences et la dispersion, la concurrence, les facteurs environnementaux, les agents pathogènes et la perturbation anthropique, pourraient influer sur la structure du peuplement (Block et Treter, 2001).

1.2.2.1. Facteurs topographiques

- Altitude : les facteurs environnementaux, par exemple la température et l'humidité du sol, changent généralement avec un gradient d'altitude (Block & Treter, 2001).

Carles et al. (2001) ont pu montrer que le climat varie essentiellement selon l'altitude et selon l'éloignement de la mer. Les températures baissent quand l'altitude augmente. Le nombre de jours de gel varie entre Novembre et Mars selon l'altitude. De même, la pluviosité annuelle augmente avec la proximité de la mer. Le déficit hydrique estival, fortement marqué en zone basse, s'atténue avec l'altitude et avec l'éloignement de la mer.

Ben Brahim et al (2004) ont noté l'action du vent sur le relief et sur la croissance végétale.

- L'exposition : Hasnaoui (2008) a souligné que la Kroumirie est caractérisée par des chaînons de montagne et des collines qui sont orientés de Sud-Ouest vers Nord-Est, ce qui leur confère une position perpendiculaire au sens des vents dominants (Nord-Ouest vers Sud - Est), qui amènent les pluies, d'où une humidité nettement plus élevée en exposition Nord qu'en celle Sud.

Selon Aussenac (2000), Les expositions fortement ensoleillées, les sols minéraux bruts, peu évolués et superficiels et les terrains pentus accentuent les conditions hydriques. Ce stress réduit la croissance en diamètre et en hauteur et affaiblit la capacité des arbres à résister aux stress.

- La pente : Les bas de versant qui ne sont pas trop abrupts et les replats sont constitués de sols d'accumulation profonds et riches.

Les versants exposés au Sud sont plus secs que ceux exposés au Nord. Les incendies y sont plus fréquents et la végétation a plus de mal à se réinstaller. L'érosion est donc plus forte et les sols sont généralement superficiels (Carles et al, 2001). D'après ces mêmes auteurs, sur les versants exposés au Nord les sols sont généralement profonds et frais, sur schistes notamment, la profondeur du sol est également fortement influencée par les différences de pentes : les sols sont superficiels et pierreux sur les pentes fortes, plus frais et meubles sur les pentes faibles.

Bergers (2004) a prouvé que la position de la plantation sur les pentes semble également exercer une influence sur la hauteur et le diamètre, soit une meilleure croissance de bas en haut.

1.2.2.2. Le dépérissement

Hasnaoui (2008) a montré que le dépérissement a touché environ 20 % des individus des différentes subéraies tunisiennes. Les subéraies de basse altitude, sur des versants exposés au Sud et sur des pentes fortes sont les plus exposées au dépérissement. L'anthropisation de la subéraie augmente encore ce phénomène.

Abdendi (2003) a révélé que l'impact d'origine anthropique de plus en plus croissant (pâturage excessif, écimage des arbres, coupes, cueillettes), qui se traduit dans la plupart des cas, par une dégradation très avancée des structures des groupements des végétaux et des sols, a rendu les arbres très vulnérables aux processus de dépérissement.

Aussenac (2000) a souligné que la sécheresse semble être le facteur le plus associé au dépérissement des forêts.

1.2.2.3 : Le surpâturage

Le surpâturage contribue à la dégradation des potentialités pastorales, la disparition de la couverture végétale, l'absence de la régénération naturelle par semis des forêts de chêne, l'apparition de l'érosion hydrique sous ces différentes formes et par conséquent le déséquilibre de tout l'écosystème Sylvo - pastoral (Hasnaoui, 2008).

D'après Boussaidi (2005), la dégradation des pâturages, d'une façon générale, est due plutôt à leur surexploitation lorsque la charge animale est excessive (la quantité de matière végétale prélevée chaque année devient supérieure à la quantité de matière consommable produite et prise sur le capital végétal). Beltrain (2002) a suggéré que les autres actions humaines qui provoquent la dégradation des subéraies sont : le tourisme (excursions, zones très visitées (les Week-end), l'urbanisation et la pollution (pluies acides, rejets toxiques)).

Dans le même sens, Bendaanoun (1998) a trouvé que les facteurs de dégradation et de destruction de la subéraie sont nombreux et leurs impacts néfastes sur le recul et la régression des écosystèmes de chêne-liège.

Ben Tiba (1980) montrait que ces défrichements sont le résultat de deux facteurs différents, mais agissant de pair, qui sont l'explosion démographique et la perte de fertilité des sols les plus soumis à l'érosion.

1.2.2.4. La pression démographique

Les densités de populations de la Kroumirie sont parmi les plus élevées pour une région rurale en Tunisie. En effet, elles ont atteint 92 habitants / Km²en 1994 contre 45 habitants / Km² en 1956, pour une région forestière dont les terres agricoles sont limitées aux clairières forestières résultant de la conversion de superficies couvertes par la végétation naturelle en terres arables par défrichement (Yacoubi, 2000).

La dynamique démographique de cette région est passée par quatre phases (Yacoubi, 2000) :

- Une phase de croissance assez lente entre 1921 et 1956 avec des taux annuels compris entre 1,16 et 1,17 % (Makhlouf, 1969) contre une moyenne nationale de 1,71 % ;

- Une phase de croissance démographique rapide entre 1956 et 1966 (2,68 % / an contre une moyenne nationale de 2,35 %) ;

- Une phase de ralentissement progressif entre 1966 et 1975 (2,11 % / an) qui s'est confirmée entre 1975 et 1984 (avec 1,97 % / an) ;

- Après 1984, ce taux est descendu à 1 % / an voire moins dans certaines zones où on commence à assister à une diminution de la population.

Tableau 1 : L'évolution démographique entre 1956 et 1994, dans la zone de Kroumirie et gouvernorat de Jendouba (I.N.S).

Selon Yacoubi (2000), les populations de la Kroumirie ont généralement des liens de parenté très forts et occupent généralement des clairières forestières. Le « Douar » compte en moyenne 25 ménages, leur ressource principale est l'exploitation du milieu environnant en pratiquant l'agriculture vivrière et l'arboriculture dans les limites de leur clairière et en exerçant leur droit d'usage dans la forêt par le biais du pastoralisme et le prélèvement du bois (pour différents usages) et d'autres menus produits.

1.3. Données générales sur le Chêne-liège

Le chêne-liège est désigné par le nom de suveraie ou subéraies données aux peuplements. C'est une essence endémique du domaine méditerranéo - atlantique du bassin méditerranéen. Il est présent en Méditerranée occidentale depuis plus de 60 millions d'années. Son aire de répartition s'est fortement réduite suite à de fortes variations climatiques et surtout suite à l'action anthropique de cette longue période (Benadid, 1989).

1.3.1. Répartition dans le monde :

Le chêne-liège est localisé presque exclusivement dans le bassin méditerranéen occidental. Sur la façade atlantique, trois taches principales : l'une en France, l'autre au Portugal et en Espagne et enfin, dans le Nord-ouest du Maroc.

Les trois autres taches le long des côtes de la Méditerranée occidentale : en France (les Pyrénées Orientales, le Var et la Corse), en Italie

(le littoral Ouest, la Sicile et la Sardaigne), en Afrique du Nord (la bordure Nord du continent en Tunisie, en Algérie et jusqu'au Rif marocain, à la rencontre des influences atlantiques).

La Subéraie mondiale serait d'environ 2.687.000 hectares repartis dans sept pays (tableau 2).

Tableau 2 : Répartition des subéraies selon la superficie (Benadid, 1989).

Pays	Superficie (ha)	%	Production de liège (T)	%
Portugal Espagne France Italie	650.000 500.000 100.000 100.000	28,5 22,0 04,4 04,4	162.000 85.000 10.000 27.000	50 26 03 08
Total Europe	1.350.000	59,2	284.000	87
Algérie Maroc Tunisie	480.000 350.000 100.000	21,0 15,3 4,4	20.000 12.000 9.000	06 04 03
Total Maghreb	930.000	40,8	4ÿÿ000	13
Total	2.280.000	100,000	325.000	100

1.3.2. Répartition en Tunisie

L'aire du chêne liège en Tunisie est localisée (confinée) dans le Nord Ouest du pays qui correspond à une unité géographique remarquable constituant principalement la Kroumirie et les Mogods (EL Afsa, 1978).

Une partie de Nefza et les forêts de Mogods constituent la subéraie orientale, celle-ci présente un faciès généralement pauvre et dégradé, par contre la Kroumirie et l'autre partie de Nefza constituent la subéraie occidentale qui présente un faciès plus productif et relativement en bon état. (DGF, 2000). Dans cette aire, les forêts de chêne-liège occupent une superficie d'environ 69.870 ha (D.G.F, 2005).

Hasnaoui (1991) signalait le chiffre de 90.000 ha .De même, Boudy (1950) avançait le chiffre 127.000 ha, répartis du NO au NE en trois grand massifs distincts. Le tableau (3) résume cette régression.

Tableau 3 : Régression de la superficie du chêne-liège en Tunisie (1950-2005)

Tableau 4 : Répartition de la superficie de peuplements purs de chêne-liège par gouvernorat (D.G.F, 2005).

Par rapport aux feuillus, le chêne-liège en Tunisie occupe plus de 53 % de la superficie des espèces forestières présentes (Figure 8). L'autre superficie est occupée par l'Eucalyptus (14 %), Acacia (5 %), Chêne zeen (7 %), feuillus mélangés (5 %) et divers feuillus (16 %).

1.3.3. Sylviculture du chêne-liège

Le traitement sylvicole est l'ensemble des opérations spéciales que l'homme fait subir à un peuplement où à un groupe de peuplements dont le but est de satisfaire les besoins de la société, d'une manière continue et soutenue (Boudru, 1989).

Bergers (2004) a signalé que le passage d'un traitement à un autre est graduel à l'intérieur d'un même régime. Des formes de transition peuvent exister entre le taillis et le taillis sous futaie, entre le taillis sous futaie et la futaie jardinée, entre celle-ci et la futaie régulière. La forêt passe par trois grands stades de développements : la jeunesse, l'âge mûr et la vieillesse.

La futaie régulière se conduit assez facilement. Les éclaircies ont pour objectif de diminuer la densité au profit des arbres les plus productifs, autour d'une à trois classes de diamètres, tout en conservant un couvert suffisant. Elles sont réalisées à chaque récolte. La densité finale préconisée est de 350 à 400 tiges par hectares.

La futaie régulière de chêne-liège présente généralement un sous bois plus propre. La régénération de ce type de peuplement est cependant plus compliquée. Les coupes d'ensemencement entraînent une explosion du maquis peu profitable au semis. Cette régénération peut être assistée (débroussaillement de trouées, crochetage).

Elle doit se faire par bouquet afin d'éviter les coupes rases trop importantes tout en tenant compte de la durée de survie des sujets. Ce type de gestion assure un revenu qui augmente en fonction de l'âge du peuplement. Cela sous-entend un trou de production de la régénération à l'âge d'exploitation (diamètre > 25 cm ; environ 30-40 ans) (I.M.L, 2005).

La futaie irrégulière demande une attention particulière. Son principe repose sur une régénération continue. Les éclaircies se font dans toutes les classes de diamètres afin de conserver un équilibre entre les jeunes sujets improductifs (diamètre < 25 cm) et les arbres productifs. La régénération se fait dans les trouées laissées lors de l'élimination des arbres improductifs. Cette gestion assure un revenu régulier, sans trou de production, grâce à un renouvellement continu du peuplement. Elle nécessite cependant des interventions sylvicoles régulières pour maintenir l'équilibre (I.M.L, 2005).

La sylviculture du chêne-liège présente deux objectifs : L'amélioration de la production du liège par la réalisation des coupes de régénération et essentiellement des coupes d'éclaircies dans les peuplements les plus denses. La protection de ces peuplements contre les maladies par l'intermédiaire des coupes sanitaires.

Chapitre 2 : Cadre naturel de l'étude

2.1 : Présentation de la région

Notre zone d'étude constitue un centre remarquable de la biodiversité terrestre notamment les formations forestières et arbustives qui se distinguent par les forêts de chêne-liège et de chêne zen. La Kroumirie est constituée par une bande étroite qui occupe presque la partie septentrionale et qui est délimitée au Nord par la Méditerranée, la frontière algérienne à l'Ouest et la ligne Ghardimaou - Fernana - Nefza - Cap Négro, au Sud et à l'Est (Ben Jamaa et al, 2006).

2.2 : Climat

La Kroumirie est bioclimatiquement la région la plus privilégiée de la Tunisie (Emberger, 1942, 1954,1955). Cette région est caractérisée par un gradient pluviométrique croissant d'Est en Ouest (de Bizerte vers Tabarka) et de bas en haut (de Tabarka vers Aïn Draham), elle est ainsi incluse dans un bioclimat humide inférieur tandis qu'Ain Draham, Ain Zena et El Feidja se situent dans l'humide supérieur (Hasnaoui, 1992).

2.2.1 : Pluviométrie

La zone d'étude est caractérisée par la présence des vents les plus fréquents de direction Nord - Ouest, perpendiculaire à l'orientation générale du relief, qui se chargent d'eau qu'ils déversent en grande partie sur ses sommets relativement élevés, d'où les effets très marqués de versant et d'exposition.

D'après les résultats du tableau ci-dessous, les précipitations sont supérieures à 1521 mm/an à Aïn Draham à 750 m d'altitude, elles sont en moyenne de 989 mm/an à Tabarka sur la côte. Selon Selmi (1985), ces précipitations varient avec l'altitude et sont mal réparties tout au long de l'année.

Tableau 5 : Pluviométrie moyenne mensuelle des deux stations Ain Draham et Tabarka (1900-2005)

Mois

Mai

Moy

ann

Tabarka

152

118,4

82,6

73,5

40,65

16,5

3,6

10,1

56,9

105,6

152,5

176,4

989

Aïn Draham

240,8

176,1

151,6

144,6

76,1

25,5

5,75

15,5

70,35

143,15

208,5

262

1521

Tableau 6 : Répartition saisonnière des précipitations des pluies (1900-2005)

D'après ce tableau, la répartition saisonnière des pluies fait apparaître que 44% des pluies qui tombent annuellement, sont concentrés pendant les trois mois d'hiver. Alors que, pendant la période sèche on en a 3% pour les deux stations.

Le régime saisonnier est du type H.A.P.E (Hiver, Automne, Printemps, Eté). Les pluies sont mal réparties tout le long de l'année, l'hiver étant la saison la plus arrosée avec 40 à 50% du total et avec un maximum au mois de Janvier et l'été la saison la plus sèche avec un minimum d'eau au mois d'Août. A ces mauvaises répartitions, s'ajoutent des fluctuations périodiques annuelles importantes (Selmi, 1985).

La comparaison des précipitations mensuelles annuelles des deux zones, durant les deux périodes (1900-1980 et (1981-2005), a montré que les quantités moyennes de pluies enregistrées à Ain Draham au cours des deux périodes sont sensiblement les mêmes, alors que la station de Tabarka a connu une légère baisse au cours de la deuxième période ; le déficit s'élève à -65 mm / an.

2.2.2 : Température

La température a une implication assez forte sur la phénologie des espèces végétales qui se trouvent dans la région (Hasnaoui, 1992).

Tableau 7° : Les températures mensuelles moyennes des deux stations Tabarka et Aïn Draham (1900-2005)

Mois

Mai

Moy

ann.

Tabarka

11,3

11,55

13,4

15,35

18,9

22,65

25,25

25,95

23,95

20,3

15,8

12,55

18,1°C

Aïn Draham

7,05

7,65

9,85

12,1

16,55

20,5

23,9

24,75

21,6

17,15

11,9

8,4

15,15C

Selon les résultats du tableau ci-dessus, la température moyenne annuelle décroît avec l'altitude 18,1°C à Tabarka sur la côte et 15,15°C à Aïn Draham à 750 m d'altitude.

Les températures maximales absolues de 47 °C à Tabarka et 43°C à Aïn Draham au cours du mois d'Août.

2.2.3 : Vents

Les vents dominants en hiver et au printemps soufflent du Nord-Ouest. En été, ce sont ceux du Sud-Est et du Sud-Ouest qui sont redoutés par les forestiers à cause de leur action asséchante, ils sont appelés Siroco (Hasnaoui, 1992).

2.2.4 : Neige

En Kroumirie, les chutes annuelles de neige ont lieu régulièrement chaque année entre janvier et février et souvent aux mois de décembre et mars dans les stations d'altitudes supérieures à 400 m et peuvent dépasser souvent un mètre d'épaisseur.

2.2.5 : Humidité relative de l'air

Ce facteur va conditionner l'évapotranspiration des plantes et donc influencer grandement leur croissance (Gerkens, 1976). L'humidité relative de l'air est très forte pendant la période d'octobre à février (85 % en janvier), moins importante en août (60 %), et très souvent, tout au long de l'année l'air est proche de la saturation (Ben Tiba, 1980).

2.2.6 : Synthèse bioclimatique :

Pour avoir une idée synthétique des paramètres climatiques, nous avons choisi, la représentation en diagrammes ombrothermiques de BAGNOULS et GAUSSEN.

Emberger a défini, à partir du quotient Q2, les étages bioclimatiques suivants : humide, subhumide, semi - aride, aride, saharien. Ces étages ont été subdivisés en inférieur et supérieur. En outre, Emberger, a fait ressortir des variantes thermiques au sein de chacun de ces étages, en se basant sur les valeurs m :

Le calcul des valeurs de Q₂:

où P est la pluviométrie moyenne annuelle en mm, M est la moyenne des maxima du mois le plus chaud et m moyenne des minima du mois le plus froid exprimée en °K), nous permet de situer les stations de Tabarka et Aïn Draham comme suit :

Tabarka se situe dans l'étage bioclimatique humide, sous étage inférieur à hiver chaud (Q₂ = 146,9 et m = 7,45°C).

Aïn Draham se situe dans l'étage bioclimatique humide sous étage supérieur à hiver tempéré (Q₂ = 202,62 et m = 4,05°C).

La lecture des diagrammes ombrothermiques (Figure 6 et 7) des deux stations montre l'existence d'une sécheresse estivale plus longue à Tabarka (4 mois) qu'à Ain Draham (2,5 mois).

Figure 6 : Diagramme ombrothermique de Bagnouls et Gaussen de la région de Tabarka (1900-2005).

Figure 7 : Diagramme ombrothermique de Bagnouls et Gaussen de la région de

2.3. Formations géologiques

La Kroumirie ainsi que les Mogods et Nefza font partie de la chaîne des Atlas septentrionaux représentée au Maroc par le Rif et le pré - Rif et en Algérie par le Tell.

Rouvier (1977) a souligné que cette zone correspond au domaine du flysch dont l'unité numidienne est la plus importante et qui s'étend de l'Oligocène au Miocène inférieur, le Numidien est un flysch argileux, gréseux comportant :

-Une série argileuse à intercalations de bancs gréseux (oligocène inférieur) les argiles sont grises, rouges, noires et parfois vertes ; la série constitue la partie basale de l'Oligocène.

-Un ensemble gréseux formé par des conglomérats successifs c'est la partie médiane (Aquitanien) et une assise marneuse au sommet (le Burdigalien).

Outre les formations argileuses et gypseuses, des formations gypseuses du Trias apparaissent à moyennes altitudes à l'ouest d'Aïn Draham et dans la plaine de Nefza (Ben Tiba, 1980).

La Kroumirie est connue en Tunisie pour l'importance de ses formations forestières liées à un substrat géologique acide. Elle se caractérise par ses reliefs, la nature des roches qui la constituent et par la mise en place de ses unités structurales.

Elle montre un relief montagneux, dont les sommets s'abaissent du Sud-Ouest au Nord-Est avec une inflexion vers le Nord, par rapport à la direction générale du tell Algérien Ouest - Est dont il est le terme final en Tunisie (Solignac, 1927). Les altitudes les plus élevées se trouvent en Kroumirie occidentale tandis que la Kroumirie orientale, à l'Est d'Aïn Draham et de Tabarka présente des altitudes plus basses, la Kroumirie centrale constitue une région aux altitudes intermédiaires (Selmi ,1985).

2.4 : Végétation

Grâce à la pluviométrie globale annuelle reçue, la Kroumirie est une région bien arrosée appartenant à l'étage bioclimatique humide, tempéré et frais. Ces conditions climatiques favorables auxquelles s'ajoutent la nature du sol et les facteurs altitude et exposition (facteurs orographiques) ont une influence sur la répartition et l'épanouissement d'une végétation naturelle et riche en taxons arborescents (Quercus suber, Quercus canariensis ; Quercus coccifera ; Pinus pinaster ; Olea europea, etc.). Ces forêts sont constituées sur plus de la moitié de leur superficie par des peuplements purs de chêne-liège et de chêne - zen et sur le reste des parties par des mélanges chêne - liège et chêne - zen. Il est aussi possible de rencontrer un mélange stable et équilibré de chêne - liège et de pin maritime.

Les études palynologiques effectuées par Ben Tiba (1980) dans la région centrale de la Kroumirie semblent apporter des éléments nouveaux à la dynamique et à la répartition de la végétation de cette région.

Ces études révèlent que la Kroumirie a connu deux types de végétations : La première est une chênaie à Quercus canariensis (chêne - zen) ayant existé avant le Pléistocène. Cette association est relayée depuis l'Holocène par une chênaie à Quercus suber (chêne-liège) qui coïncide, en outre, avec l'apparition des Ericacées notamment Erica arborea, Erica scoparia et Arbutus Unedo.

Chapitre 3 : Matériel et Méthodes

3.1. Matériel végétal

L'étude a touché à la subéraie de la Kroumirie. Nous nous sommes limités aux peuplements purs de chêne-liège.

3.1.1. Exigences écologiques du chêne-liège

Le chêne-liège est un arbre caractéristique de la région méditerranéenne occidentale. Il se développe sous un climat typiquement méditerranéen qu'Emberger (1930) définit comme étant un climat de la zone extratropicale à photopériodisme quotidien et saisonnier et à pluviométrie concentrée sur la période froide de l'année, l'été, saison la plus chaude, étant sec.

Tous les chênes sont inclus dans les variantes à hivers tempéré, doux et chaud de l'étage bioclimatique humide d'Emberger et débordent dans l'étage subhumide supérieur à hivers tempéré à doux. Le chêne - liège se comporte dans le bassin méditerranéen comme une espèce relativement plastique. Il apparaît en Tunisie et en Algérie jusqu'à 1300 m d'altitude et au Maroc jusqu'au 3200 m dans le grand Atlas (Boudy, 1952). La plasticité écologique du chêne-liège par rapport à la pluviosité est grande. En Tunisie, il a besoin d'une tranche pluviométrique de 550 mm à 600 mm/an au minimum. L'humidité de l'air notamment, qui atteint des valeurs élevées dans les régions proches du littoral, permet à cette espèce pionnière d'occuper largement les sols sablonneux ou gréseux et les terres rouges (Bellkhadar, 2003).

Il est également exigeant en ce qui concerne la chaleur, l'espèce est dite thermophile, espèce assez frileuse (assez sensible au froid elle ne supporte pas de froid prolongé inférieurs à 5°C ou de froid de courte durée inférieur à -9°C) par contre il supporte bien les grandes chaleurs, sous réserve qu'elle soit accompagnée d'une humidité suffisante (Ben Seghir, 1996). Le préférendum thermique du chêne - liège oscille dans un intervalle de 13 à 18°C. Ce même auteur a suggéré que le chêne-liège est une espèce héliophile à l'âge adulte, les semis ont besoin d'ombre pour éviter les brûlures de l'apex c'est pourquoi les jeunes glands poussent souvent dans les touffes du maquis (Et tobi, 1996). Le sol des subéraies est largement acide ou voisin de la neutralité entre 6 et 6,5 de pH. La répartition limitée du chêne-liège est aussi en relation avec ses exigences édaphiques, il est calcifuge, il ne supporte ni le sous sol calcaire ni les chlorures (Natividade, 1956). Yessad (2000) affirme que le chêne-liège prospère sur les terrains siliceux à horizon supérieur sableux qui lui est vraisemblablement un caractère général et d'une nécessité édaphique fondamentale.

Les forêts de chêne-liège offrent généralement une structure pluri - strates : Le sous - bois est dense et très varié et la strate herbacée est bien développée.

Selon l'altitude et la nature du sol, on peut dégager plusieurs types structuraux (physionomiques) et groupement floristiques (Schoenenberger, 1962). La structure verticale et horizontale du peuplement de chêne-liège introduit une grande hétérogénéité dans l'éclairement du sol. Ainsi sous la subéraie peut se développer un sous bois assez riche ; arbuste et petits ligneux du maquis, nombreuses espèces herbacées, composant à la fois des espèces siophiles (d'ombre) et des espèces héliophiles (de lumière) (Amandier, 2002).

3.1.2. Caractéristiques botaniques

Le chêne-liège se présente sous la forme d'un arbre de taille moyenne pouvant atteindre 20 à 25 mètre de hauteur, mais le plus souvant ne dépassant pas 15 mètre à l'âge adulte et vivre jusqu'à 500 ans à l'état naturel. Par contre lorsqu'il est exploité, la durée de vie de l'arbre varie entre 150 et 200 ans. Son feuillage est peu abondant et son couvert est léger facilitant le développement d'un sous bois de type maquis très dense (Cantat & Piazzetta, 2005).

Ces mêmes auteurs ajoutent que le chêne-liège est un arbre à feuilles persistantes (2 à 3 ans) bombées, de forme ovale, au limbe quelque peu denté. Hasnaoui (2002) et Bellakhder (2003) évoquent le fait que le chêne-liège est une espèce monoïque, les fleurs mâles sont regroupées en grappes, appelées (Châtons), au bout des pousses de l'année précédente, les fleurs femelles poussent isolément à la base des feuilles de la pousse de l'année. La floraison a lieu de la fin d'avril jusqu'à la fin de mai. Le fruit, connu sous le nom de gland est une capsule ovoïde enchâssée à la base dans une cupule hémisphérique et contenant une amande.

Ksontini (1996) souligne que le chêne-liège possède un enracinement fort grâce à une racine pivotante, qui se développe en profondeur ainsi que, de nombreuses racines latérales, ce qui donne à l'arbre une bonne alimentation en eau.

Le chêne-liège se multiplie par ensemencement naturel, semis ou plantation, rejets de souches et drageonnement (Hasnaoui, 1991).

3.2. Méthodologie

3.2.1. Cartographie

A partir des cartes de peuplements contenues dans les procès verbaux d'aménagements de chaque zone, nous avons dégagé les parcelles qui représentent les peuplements de chêne-liège purs couvrant la partie centrale de la Kroumirie (forêts de Tabarka et Ain Draham).

Cette phase de cartographie nous a permis de dresser les différentes cartes nécessaires pour ce travail.

3.2.2. Inventaire de l'état actuel de la subéraie en Kroumirie

3.2.2.1. Pré.échantillonnage et calcul statistique

Afin de déterminer le nombre de placettes à retenir dans notre zone d'étude, à partir d'un échantillonnage aléatoire, nous avons réalisé une phase de pré échantillonnage. Elle a concerné les stations suivantes : Hamdia I, Mekna IV, Mekna V, Mekna II, Ain Saida, Ain Zena III, Tbeinia, Souiniet, Ain Bacouch et Amdoun I en essayant de couvrir la majorité de la zone d'étude.

Nous avons fixé trois placettes par station et pour une dizaine de stations afin de tenir compte, au maximum de l'hétérogénéité des peuplements.

Chaque placette ayant 500 m² de superficie, de forme circulaire. Pour déterminer le nombre des placettes dans cette étude, nous nous sommes basés sur le nombre d'arbres à l'hectare (ou densité) de chaque placette.

Avec e= erreur =6% et 1,96 = la constante de la loi normale fixée auparavant pour un pré échantillonnage de 30 placettes.

Station	Nombre de placette	Densité à l'hectare de chaque placette
Mekna IV	Placette 1	360
	Placette 2	420
	Placette 3	400
Mekna II	Placette 1	220
	Placette 2	240
	Placette 3	200
Ain Saida	Placette 1	300
	Placette 2	400
	Placette 3	460
Ain Zena III	Placette 1	760
	Placette 2	700
	Placette 3	540
Souiniet	Placette 1	560
	Placette 2	520
	Placette 3	480
Tbeinia	Placette 1	660
	Placette 2	560
	Placette 3	600
Ain Bacouch	Placette 1	180
	Placette 2	220
	Placette 3	200
Amdoun I	Placette 1	360
	Placette 2	340
	Placette 3	280
Hamdia I	Placette 1	220
	Placette 2	160
	Placette 3	180
Mekna V	Placette 1	380
	Placette 2	480
	Placette 3	460

Tableau 9 : Tableau récapitulatif du calcul statistique de pré.échantillonnage

A partir de ce calcul, nous constatons que nous pouvons prendre un nombre de 239 placettes, qui vont servir comme unité d'échantillonnage représentant notre zone d'étude.

3.2.2.2. Echantillonnage

En tenant compte de la superficie de la région de la Kroumirie et de la nature du peuplement (peuplement pur de chêne-liège), nous avons opté pour un échantillonnage aléatoire. Nous avons matérialisé de manière aléatoire, des placettes de forme circulaire et de 500 m² de superficie. Les cordonnées du centre de chaque placette ainsi que l'altitude et l'orientation ont été relevés par le GPS. Ces positions géographiques ont été, par la suite, portées sur une carte au 1/50.000. Le nombre de placettes pris en compte et déterminé par le pré.échantillonnage est de 239.

3.2.2.3. Description des stations

Pour décrire nos placettes, nous avons élaboré des fiches que nous avons remplies lors de nos visites à ces placettes (voir annexes). Dans le tableau 10, figurent les caractéristiques écologiques sommaires de nos stations.

Station	Nombre de placette	Caractéristiques écologiques
Mekna IV	15	Exposée Est à Nord Est, elle évolue de 200 à 300 m d'altitude. Elle relève de l'étage bioclimatique humide inférieur à hiver doux. L'essence principale est le chêne-liège avec présence d'espèces de maquis comme Pistacia lentiscus, Myrtus communis et Erica arborea.
Mekna II	21	D'exposition générale Est et de 150 à 250 m d'altitude, elle a les mêmes caractéristiques bioclimatiques que la station précédente. La végétation est composée essentiellement de Quercus suber à l'état pur avec un cortège d'espèces de maquis dont : Erica arborea, Arbutus unedo, Olea europea et Myrtus communis.
Ain Bacouch	27	Exposition générale Nord à Nord Ouest. Elle est située de 100 à 250 m d'altitude. Les caractéristiques climatiques sont celles enregistrées à Tabarka. Le chêne-liège existe à l'état pur accompagné de certaines espèces de maquis comme Myrtus communis, Erica arborea et Smilax aspera.
Hamdia I	44	D'exposition Nord à Nord Ouest et de 100 à 200 m d'altitude, elle appartient à l'étage bioclimatique humide inférieur. Le peuplement de chêne-liège est pur avec maquis de différentes espèces comme Erica arborea, Pistacia lentiscus et Arbutus unedo.
Amdoun I	28	Exposition Est à Nord Est. Elle se trouve de 300 à 500 m d'altitude et appartient à l'étage bioclimatique humide à hiver doux. Le chêne-liège est la seule essence arborescente.
Souiniet	24	Située sur un plateau de 600 à 700 m d'altitude. D'exposition générale Nord à Nord Est, elle révèle de l'étage bioclimatique humide supérieur à hiver tempéré. La végétation est composée de Quercus suber avec un maquis riche en Erica arborea, Smilax aspera, Myrtus communis et Arbutus unedo.
Tbeinia	20	Exposée vers le Sud, elle se situe entre 600 et 800 m d'altitude. Elle possède les mêmes caractéristiques climatiques qu'Ain Draham. Le peuplement de chêne-liège est pur associé à un maquis très dense composé essentiellement de Pistacia lentiscus, Erica arborea, Arbutus unedo et Smilax aspera
Ain Saida	42	D'exposition Nord Ouest, elle est située de 350 à 580 m d'altitude. Elle appartient à l'étage bioclimatique humide supérieur à hivers tempérés. Le chêne-liège pur est accompagné de : Erica arborea, Myrtus communis et Arbutus unedo.
Ain Zena III	08	Située à une altitude de 900 à 1000 m, elle a les mêmes caractéristiques bioclimatiques que la station d'Ain Draham. Le chene-liege pur est accompagné d'un maquis composé essentiellement de Pistacia lentiscus et Erica arborea.
Mekna V	10	Située à environ 380 à 450 m d'altitude et d'exposition générale Nord Est. Peuplement pur de chêne-liège avec maquis composé essentiellement de Myrtus communis et Erica arboreaa.

3.2.2.4. Description de la structure actuelle de la subéraie en Kroumirie

Au cours de nos investigations sur le terrain, nous avons déterminé pour chaque placette les caractéristiques dendrométriques individuelle de l'espèce étudiée (Diamètre à 1,3 m en cm et Hauteur en m). Au niveau de chaque placette, nous avons calculé le diamètre moyen en cm, la hauteur moyenne en m et la densité à l'hectare.

Pour présenter la structure actuelle de la subéraie, nous avons choisi les classes de diamètre réparties comme suit :

Tableau 11 : Classes de diamètres retenues pour représenter la structure de la subéraie

3.2.2.5. Effet des facteurs du milieu sur la densité et la structure diamétrique du chêne-liège en Kroumirie

Pour étudier l'effet des facteurs du milieu sur la répartition des différentes classes de diamètres pour chaque station, nous nous sommes appuyés sur les données de la fiche d'observation des placettes (voir annexes).

Cette fiche regroupe des facteurs comme les caractéristiques géographiques de la station (exposition, altitude, pente et distance des peuplements aux douars les proches); les caractéristiques sylvicoles de la station (origine du peuplement, structure du peuplement, stade de croissance, essences principale et secondaire avec leur proportion), les caractéristiques dendrométriques des arbres échantillonnés, densité de chaque placette (nombre d'arbre par 500 m²), le nombre de plants de l'espèce principale par placette (régénération), l'intensité de pâturage et la présence des incendies. De même, nous avons estimé les distances séparant les placettes des douars environnants.

Les facteurs édaphiques n'ont pas été pris en considération parce qu'ils sont supposés assez homogènes dans nos différentes stations et que leurs impacts dans les variations spatio - temporelles des structures des peuplements intra et inter - stationnelles sont minimes.

En ce qui concerne les facteurs climatiques, nous avons utilisé les données météorologiques des stations les plus proches des nôtres (Tabarka et Ain Draham).

3.2.3. Dynamique de l'écosystème chêne-liège en Kroumirie

3.2.3.1. Dynamique de la subéraie en Kroumirie en fonction des facteurs du milieu.

Pour étudier cette dynamique, du point de vue structure des peuplements, plusieurs méthodes à suivre sont possibles.

Selon Foundi et al.(2001), ces méthodes peuvent être subdivisées en trois catégories ; les inventaires forestiers, les méthodes diachroniques et la méthode de la FAO.

Les approches basés sur les inventaires périodiques sont considérées comme étant les méthodes les plus précises pour évaluer la dynamique. Cependant, elles exigent une standarisation des normes et des méthodes de prélèvement des données.

Les méthodes diachroniques consistent à suivre l'évolution de plusieurs aspects d'une parcelle au cours du temps par la comparaison de deux plans cadastraux ou deux couvertures photographiques de divers époques. Le principe de la méthode est basé sur l'exploitation des données planimétriques (la densité ou le couvert), soit d'une façon systématique ou par échantillonnage (méthode de grilles) et des outils statistiques tels que les matrices de succession. L'analyse des résultats des comparaisons permet d'exprimer par rapport à la période étudiée, un rythme moyen d'évolution qui se traduit par des pertes ou des gains de surfaces boisés d'une certaine densité ou d'un certain niveau de recouvrement.

En outre, cette méthode exige également une standarisation des critères de stratifications des peuplements à différentes dates de comparaison.

La méthode de FAO, a été utilisée pour évaluer les ressources forestières tropicales. Elle consiste à élaborer un modèle de déforestation qui rattache le changement de la couverture forestière au temps et d'autres variables telles que la densité de la population et le taux de croissance démographique pour la période correspondante.

Dans ce travail, nous nous proposons d'étudier l'évolution de la densité des peuplements purs de chêne-liège en Kroumirie à différentes dates et en fonction des facteurs de perturbations (orographiques et anthropiques) selon la méthode des inventaires périodiques, quoiqu'il ne s'agit pas, pour les dates de départ (1983) et de fin (2008), de résultats d'inventaire proprement dits, mais de résultats de recherches.

Ainsi, et pour l'année 1983, nous nous sommes appuyés sur les données dendrométriques et les résultats d'analyses statistiques des différentes variables du milieu (facteurs biotiques et abiotiques) présents dans les travaux de Hasnaoui, 1992. Pour l'année 2008, nous avons utilisé nos propres données.

Les données intermédiaires sont tirées des inventaires nationaux forestiers de 1993 et 2005.

3.3. Analyses statistiques

Nos données ont été traitées en utilisant des calculs statistiques simples (Ecart - type, Covariance, Moyenne, etc.) qui ont permis de décrire la situation actuelle.

En ce qui concerne l'effet des facteurs du milieu sur la structure, le calcul de la variance a été effectué en utilisant le logiciel «SAS». La comparaison des moyennes a été effectuée en utilisant le test (SNK). Les graphiques et les courbes ont été réalisés à l'aide du logiciel « Excel ».

Chapitre 4 : Résultats

Dans ce chapitre, et en un premier temps, nous avons tenu, à travers nos mesures de l'année 2008, à donner un constat de la situation actuelle de la structure de la subéraie tout en le rapportant à certains facteurs topographiques et anthropiques du milieu.

La topographie (altitude, exposition et pente) traduit des facteurs écologiques notamment la température et la pluviométrie et d'autres édaphiques, exprimés essentiellement par la pente, comme la profondeur du sol et son bilan hydrique.

Les facteurs anthropiques, représentés ici par la présence humaine et le cheptel (éloignement de la forêt par rapport aux douars), s'expriment surtout par le surpâturage et son impact sur la régénération naturelle et sur les nouvelles pousses d'une manière générale, ainsi que par les délits éventuels touchant notamment aux stades gaulis, perchis et adulte de l'essence étudiée.

Dans une deuxième étape, nous avons procédé à la comparaison de nos données avec celles d'études antérieures à savoir les données du début des années 1980 (Hasnaoui, 1992) et celles des inventaires Nationaux Forestiers de 1995 et 2005.

Cette comparaison vise la détermination de l'évolution de la structure diamétrique de ces subéraies en vue d'essayer d'en dégager une idée sur leur avenir.

Il est à souligner et comme nous l'avons mentionné au chapitre 3, que les facteurs édaphiques (pH, structure et texture du sol, etc.) ont été neutralisés puisque nous avons travaillé dans le même type de sol (moder) et que notre objectif principal est d'apprécier la dynamique de la subéraie. Concernant les facteurs climatiques, les stations de références sont celles de Tabarka et Ain Draham, et étant donnée l'absence de stations en milieu forestier, nous avons procédé à des corrélations entre les placettes en fonction de l'altitude et l'exposition (La pluviométrie augmente, dans la région, de 0,7 mm / m d'altitude, l'exposition Nord est plus arrosée que l'exposition Sud, et la température diminue avec l'altitude).

4.1. Constat de la structure diamétrique actuelle de la subéraie en Kroumirie

4.1.1. Description de la structure diamétrique actuelle

La structure des stations forestières constitue un élément d'information primordial de leur écologie et de la dynamique des diverses essences forestières. Elle traduit le mode d'organisation de l'écosystème forestier et de ses composants (Hasnaoui, 1992).

Dans cet esprit, nous avons présenté la répartition des différentes classes diamètriques (structure diamétrique) des peuplements de chêne-liège purs en Kroumirie dans le cadre de nos stations d'étude. Ces classes diamètriques correspondent, conventionnellement, aux stades de semis, de gaulis, de bas et de haut perchis et d'adultes.

Les résultats des densités des individus par classe de diamètre et par station sont consignés dans la figure 9.

La lecture de la figure montre que les semenciers ou classe de diamètre D5 sont mieux représentés dans les stations alticoles Mekna V (87 individus/ha), Ain Zena III (88,75 individus/ha), Souiniet (95 individus/ha) et Tbeinia (103 individus/ha). Par contre, les densités faibles sont enregistrées dans les stations basses de Hamdia I (45 individus/ha), Mekna II (33 individus/ha) et Ain Bacouch (65 individus/ha).

La densité de régénération est très élevée dans les stations alticoles de Mekna V, Ain Zena III, Tbeinia et Souiniet respectivement 275, 430, 325 et 315 individus/ha. Les densités assez faibles sont présentées dans les stations basses : Hamdia I (150 individus/ha), Mekna II (178 individus/ha) et Ain Bacouch (143 individus/ha).

En ce qui concerne les individus relevant des classes intermédiaires D2, D3 et D4, les densités sont très rapprochées pour la totalité des stations d'étude. Les fortes densités sont enregistrées dans la station de Souiniet, soient 23 individus/ha (D2), 27 individus/ha (D3) et 35 individus/ha (D4), alors que, les densités les plus faibles sont rencontrées dans la station de Hamdia I, soient 9 individus/ha (D2), 10 individus/ha (D3) et 17 individus/ha (D4).

Les données recueillies de l'ensemble des placettes et portées dans la figure 10 montrent que les classes les mieux représentées dans cette subéraie sont respectivement D1, avec un effectif de 250 individus/ha, soit 65,78 % de l'effectif total, et D5 avec un taux de 18,42 %, soit 70 individus/ha. Par contre, les classes D2, D3 et D4 sont très peu représentées (respectivement 4 %, 5 % et 6,80 %).

Les densités faibles des classes intermédiaires : gaulis, bas et haut perchis sont inquiétants et traduisent les difficultés rencontrées par la subéraie en Kroumirie à se reconstituer naturellement par semis.

Figure 10 : Structure diamétrique globale (toutes stations confondues) de la subéraie en Kroumirie (en 2008)

4.1.2. Impacts de certains facteurs sur la structure diamétrique actuelle

4.1.2.1. Effet de l'altitude sur la densité et la structure diamétrique du chêne-liège

Ce constat revêt, en premier lieu, un aspect général où nous avons travaillé sur des moyennes de données se rapportant à l'ensemble de la zone d'étude, toutes stations confondues et, en second lieu, nous nous sommes orientés vers plus de détails, en traitant séparément les placettes de chacune des stations, les paramètres observés et / ou mesurés étant les mêmes.

Les représentations graphiques ne portent que sur les paramètres dont les effets sur la structure sont significatifs d'après les résultats de l'analyse de variance figurant sur les tableaux correspondants.

Les analyses statistiques (analyse de variance) appliquées à nos données dendrométriques, récoltées en 2008, ont abouti aux résultats figurant dans le tableau 13.

Tableau 13 : Résultats de l'analyse de variance des effets de l'altitude sur la répartition des classes de diamètre (2008)

En effet, la variance est très hautement significative (P< 0,001) pour les diamètres les plus bas (D1 et D2), alors qu'elle est seulement significative pour les troisième et quatrième classes (D3 et D4) et elle est non significative pour la classe des semenciers (D5) (tableau 13).

La lecture de la figure 11 montre ainsi que la classe D1, correspondant aux plantules ayant des diamètres inférieurs à 5 cm et représentant la régénération naturelle de l'espèce, est présente avec des effectifs importants (autour de 295 individus / ha), dans les placettes appartenant à la troisième frange d'altitude (> 600 m).

Entre 300 et 600 m d'altitude la densité des individus de D1 tourne autour de 160 / ha, loin derrière la troisième classe d'altitude.

L'effectif le plus faible est rencontré au niveau des placettes appartenant à la frange d'altitude la plus faible.

Concernant la classe D2 (diamètres compris entre 5 et 10 cm correspondant aux stades fourrés et gaulis), elle reste assez bien représentée dans les placettes alticoles, avec environ 17 individus / ha, suivies des placettes de moyennes altitudes et enfin de celles les moins élevées.

Les individus relevant des classes D3 et D4 ont des effectifs semblables mais assez faibles, avec les mêmes tendances que pour les classes précédentes quant à l'effet de l'altitude.

Il ressort de cette figure, que les semis s'installent mieux dans les stations alticoles qui sont les mieux arrosées, mais il est fort probable que le facteur eau n'est pas le seul responsable de cette distribution et que la prédation participe, pour beaucoup, dans ce constat (Hasnaoui, 1992, 1995, 1998, Hasnaoui, 2002). La même remarque peut être faite pour les classes de diamètres D2, D3 et D4.

Quant aux effectifs des semenciers (D5), nous constatons qu'ils n'ont pas été affectés significativement par l'altitude.

4.1.2.2. Effet de l'exposition sur la densité et la structure diamétrique du chêne-liège

Le relief de la Kroumirie est constitué de chaînons de montagnes et de collines, faisant partie de l'Atlas septentrional qui débute au Maroc et se termine au Nord de la Tunisie, du côté de Bizerte. La direction générale de cet ensemble est de Sud - Ouest à Nord - Est.

Les vents dominants de la région viennent du Nord - Ouest et se dirigent vers le Sud - Est, soit perpendiculairement à la direction du relief.

Les nuages amenés par ces vents, et notamment les nuages bas qui peuvent donner des précipitations, sont donc arrêtés par cette barrière naturelle, d'où la chute de pluie.

Il en résulte que l'exposition Nord, d'une manière générale et par notre région d'étude, est plus arrosée que celle Sud, d'où l'effet probable sur la végétation naturelle et sur l'écosystème dans son ensemble.

Tenant compte de cette situation, nous avons noté, lors de la récolte des données, l'exposition de chaque placette (nous avons considéré uniquement les expositions Nord et Sud), afin de la corréler avec les paramètres observés, à savoir la classe de diamètres et la densité du peuplement.

Les résultats de l'analyse de variance des effets de l'exposition sur la répartition des classes de diamètres de l'essence en question, sont cosignés dans le tableau 14.

Tableau 14 : Résultats de l'analyse de la variance des effets de l'exposition sur la répartition des classes de diamètre

Ce tableau montre l'existence d'un effet très hautement significatif (P<0,001) entre l'exposition et le nombre d'individus / ha pour les classes de diamètres D1 et D2.

Dans les stations exposées au Nord, les semis de chêne-liège (D1) sont assez denses, puisqu'on compte 295 individus / ha, par contre dans les stations exposées au Sud, ces effectifs sont beaucoup plus réduits, puisqu'on ne rencontre que 91 individus / ha, soit environ le 1/3 de ceux rencontrés dans l'exposition opposée (figure 12).

Quant aux individus relevant de la classe (D2) de diamètres, ils présentent des effectifs de loin inférieurs à ceux de la classe (D1).

Ainsi, pour l'exposition Nord, cette classe est représentée par 19 individus / ha, soit 6 % de l'effectif des individus appartenant à la classe D1, pour la même exposition, et par 4 individus / ha pour l'exposition Sud, soit 4 % de ceux de D1.

Si on compare les densités des individus de cette classe dans les deux expositions en question on constate comme pour D1, que les individus exposés au Sud, ne représentent que 45% de ceux exposés au Nord.

Les trois classes les plus grandes (D3, D4 et D5) ne subissent aucun effet significatif de l'exposition.

4.1.2.3. Effet de la pente sur la densité et la structure diamétrique du chêne-liège

Le degré d'inclinaison d'un terrain ou la pente est le troisième paramètre orographique d'une station donnée, qui peut avoir des effets sur l'installation et la croissance des végétaux.

L'action de ce paramètre peut s'exprimer par le biais de deux types de facteurs écologiques du milieu d'une influence capitale sur l'existence des plantes naturelles et leur répartition :

- La première catégorie de ces facteurs sont les facteurs édaphiques qui portent sur la

En effet, le degré d'inclinaison du terrain, associé au taux de couverture végétale de ce terrain (couvertures muscinale, herbacée, arbustive et arborescente), vont commander la mobilité et le dépôt des substrats solides (transport du sol) sous l'effet de l'eau et de la pesanteur.

Plus la pente est forte et plus le terrain est découvert, plus le transport est favorisé et donc le sol qui reste est plus mince ou peu profond et plus pauvre.

Par contre, lorsque la pente est faible ou nulle, le dépôt est favorisé et le sol est plus profond et plus riche.

Ce phénomène est d'autant plus important que les précipitations sont de fortes intensités, ce qui constitue une caractéristique du climat méditerranéen, sous lequel existe notre écosystème.

- Le deuxième facteur écologique, fortement lié aux précédents est le bilan hydrique. Ce bilan est d'autant plus positif et favorable que la pente est faible et le terrain est couvert.

Dans ce cas, l'infiltration de l'eau est favorisée (nature du sol aidant), par contre, plus le terrain est incliné et nu plus le ruissellement l'emporte, causant des départs importants de substrats.

L'analyse statistique de l'effet de la pente sur la répartition des classes de diamètres des arbres de chêne-liège a été réalisée et elle a donné le tableau 15.

Tableau 15 : Résultats de l'analyse de variance des effets de la pente sur la répartition des classes de diamètre

D'après le tableau (15), il existe un effet hautement significatif (P < 0,01) pour D2 et très hautement significatif (P < 0,001) entre la pente et le nombre d'individus/ha des individus de classe D1.

La figure 13 montre que les effectifs des semis les plus élevés sont situés au niveau des placettes de faibles pentes (297 individus/ha).

Dans les reliefs de pentes moyennes, la fréquence des plantules de chêne-liège est de 178 par hectare, soit 55 % de ceux des terrains plats (P1).

Au niveau des endroits où les pentes sont fortes, la densité des semis est très faible (18 individus/ha), elle est nulle au niveau des terrains à très fortes pentes.

En ce qui concerne les individus de classes de diamètres D2, les fortes densités sont enregistrées au sein des placettes à pentes faibles ou modérées, respectivement 23 et 16 individus/ha.

Par ailleurs, la densité est minimale dans les terrains à pentes fortes (8 individus/ha) et est nulle aux endroits à pentes très fortes.

4.1.2.4. Effet du facteur « distances du peuplement aux douars » les plus proches sur la densité et la structure diamétrique du chêne-liège

Le pâturage en forêt entraîne une dégradation progressive de la végétation arborée. Tous les animaux se nourrissent, essentiellement, de jeunes pousses, empêchant ainsi toute régénération. On constate l'existence de plusieurs causes de déforestation qui peuvent agir simultanément ou chacune à part.

Selon Hasnaoui (1995), la prédation, en Kroumirie, peut avoir lieu au niveau du gland de chêne-liège comme elle peut l'être au niveau de la plantule. Concernant le gland, plusieurs agents peuvent intervenir dont notamment : l'Homme (ramassage des glands que ce soit par gaulage ou après leur chute naturelle), les animaux domestiques, le sanglier, les oiseaux et les rongeurs.

Cette destruction des forêts de chêne-liège provoque des perturbations très graves des écosystèmes naturels, induisant ainsi l'augmentation de l'impact de certains désastres réputés naturels comme la sécheresse.

Dans ce contexte il est important de rappeler que la totalité des placettes inventoriées dans notre zone d'étude n'est pas protégée.

A partir des données recueillies au niveau des placettes, nous avons dressé l'effet des facteurs anthropiques sur les peuplements de chêne-liège en fonction de la distance aux douars les plus proches. Le tableau 16 donne les résultats de l'analyse de variance des effets de ce facteur.

Tableau 16 : Résultats de l'analyse de variance des effets de la distance des peuplements par rapport aux douars les plus proches sur la répartition des classes de diamètres

D'après le tableau (16), il existe un effet très hautement significatif (P < 0,001) de l'effet de la distance des douars les plus proches sur la densité des individus relevant des classes D1 et D2.

De même, il existe une interaction significative (P < 0,05) pour les arbres de classes de diamètres D3, D4 et D5.

La figure 14 montre une domination complète, du point de vue effectifs, des jeunes plantules de chênes-lièges (D1) dans les placettes les plus éloignées des douars (299 individus/ha) par rapport à ceux moyennement éloignées (130 individus/ha) et les plus proches (25 individus/ha).

Cette domination des effectifs dans les lieux les plus éloignés se confirme aussi pour les individus du stade gaulis (D2) où on compte respectivement, 20 pieds/ha contre 6 dans les peuplements les plus proches des douars.

En consultant la même figure, nous constatons, d'une manière générale, que les effectifs des individus relevant des autres classes de diamètres (D3, D4 et D5) rencontrés dans les lieux les plus éloignés, sont beaucoup plus élevés que ceux trouvés aux alentours des habitations.

En conclusion, on peut dire que le facteur anthropique joue un rôle majeur dans la perturbation du milieu, les stations les plus éloignées des douars sont les plus denses, par contre, les stations les plus proches des douars sont les plus pâturées et les moins denses.

4.1.3. Variation de la densité du chêne-liège en fonction de certains facteurs du milieu pour chaque station.

4.1.3.1. Répartition de la densité des individus, toutes classes de diamètres confondues, en fonction de l'altitude.

L'altitude a une influence directe sur la pluviométrie et la température qui, à leur tour, ont un impact direct sur la répartition et la distribution des végétaux.

Les stations de notre étude sont réparties sur trois classes d'altitudes (basse, moyenne et haute altitude).

Les analyses statistiques (analyse de variance) appliquées à nos données, ont abouti aux résultats figurant dans le tableau 17.

Tableau 17 : Résultats de l'analyse de variance des effets de l'altitude sur la densité du chêne-liège (toutes classes d'altitudes confondues)

	Station	Fobs	P>F	Sign
Basses altitudes	H I	2,21	0,05	*
	AB	2,87	0,05	*
	M II	2,89	0,05	*
altitudes moyennes	M IV	3,79	0,01	**
	AS	3,28	0,05	*
	AM I	3,11	0,05	*
	M V	4,26	0,001	***
Hautes altitudes	AZ III	6,54	0,001	***
	TB	5,27	0,001	***
	SO	5,56	0,001	***

Le tableau 17 montre l'existence de relations significatives (5 %) pour les stations Hamdia I (H I), Mekna II (M II), Ain Baccouch (AB), Amdoun I (AM I) et Ain Saida (AS), hautement significatives (1%) pour la station M IV et très hautement significatives (0,1 %) pour les stations Mekna V (M V), Ain Zena III (AZ III), Tbeinia (TB) et Souiniet (SO).

La figure 15 montre que la densité de chêne-liège varie entre un minimum de 235 individus/ha pour la station de Hamdia I et un maximum de 588 individus/ha pour la station d'Ain Zena III.

Pour les stations de basse altitude, la densité varie d'une station à une autre. Elle atteint 235 individus/ha pour Hamdia I, 253 et 246 individus/ha, respectivement pour les stations de Mekna II et Ain Bacouch.

Dans les stations de la deuxième frange d'altitude, les effectifs sont élevés et varient entre 426 individus/ha pour la station de Mekna V et 267 individus/ha dans la station d'Amdoun I.

Les stations alticoles sont les plus denses, ainsi, la densité moyenne passe de 490 à 494 et à 588 individus/ha respectivement pour les stations de Tbeinia, Souiniet et Ain Zena III.

A partir de ces informations, on remarque une nette augmentation de la densité du chêne-liège avec l'altitude.

La lecture de ces résultats révèle que la densité de chêne-liège varie à l'intérieur d'une même classe d'altitude, cette variation peut être le résultat de l'effet d'autres facteurs du milieu.

4.1.3.2. Répartition de la densité en fonction d'exposition, toutes classes de diamètres confondues et par frange d'altitude.

En Tunisie, on constate l'existence d'une différence considérable entre le versant Nord, relativement froid et humide et le versant Sud relativement chaud et sec, conséquence directe du relief (altitude et exposition).

Les résultats des analyses de variance de l'effet de l'exposition sur la densité des peuplements de chêne-liège, toutes classes de diamètres confondues et à basse altitude sont consignés dans le tableau 18.

Tableau 18 : Résultats de l'analyse de variance des effets de l'exposition sur la densité du chêne-liège, à basse altitude et toutes classes de diamètres confondues

Selon le tableau 18, il existe une relation hautement significative (1%) pour les placettes d'exposition (Nord) et significative (5%) pour les placettes d'exposition (Sud) des deux stations Hamdia I (H I) et Ain Bacouch (AB).

La lecture de la figure 16 montre que les fortes densités sont enregistrées dans les placettes d'exposition Nord des deux stations H I (245 individus/ha) et A B (235 individus/ha).

Par ailleurs, on constate des densités relativement faibles sur les versants Sud, soit 85 individus/ha (H I) et 135 individus/ha (AB).

Le tableau 19 porte les résultats de l'analyse de variance de l'effet de l'exposition au niveau des stations de moyennes altitudes.

Tableau 19 : Résultats de l'analyse de variance des effets de l'exposition sur la densité du chêne-liège, à moyenne altitude (toutes classes de diamètres confondues)

	Station	Fobs	P>F	Sign
Exposition (Nord)	M IV	1,08	0,44	NS
	AS	4,53	0,001	***
	AM I	3,53	0,01	**
	M V	1,29	0,27	NS
Exposition (Sud)	M IV	0,71	0,55	NS
	AS	2,99	0,05	*
	AM I	2,53	0,05	*
	M V	0,96	0,48	NS

Il existe des interactions très hautement significatives (P<0,001) pour Ain Saida (AS), hautement significative (P<0,01) pour Amdoun I (AM I) pour une exposition Nord et significative (P<0,05) pour les deux stations à une exposition Sud (tableau 19).

L'exposition semble avoir une influence marquée sur la densité de chêne-liège, les fortes densités abondent dans les versants Nord des stations (AMI) et (AS), respectivement 295 et 339 individus/ha.

Par contre, dans les placettes d'expositions Sud, la densité des peuplements est faible, elle atteint 256 individus/ha (AS) et 180 individus/ha (AM I) (figure 17).

Le tableau 20 porte les résultats de l'analyse de variance des effets de l'exposition sur la densité du chêne-liège, toutes classes de diamètres confondues au niveau des stations alticoles.

Tableau 20 : Résultats de l'analyse de variance des effets de l'exposition sur la densité du chêne-liège, à haute altitude (toutes classes de diamètres confondues)

Il existe des relations très hautement significatives (0,1%) pour les deux stations de Tbeinia (TB) et Souiniet (SO) en exposition Nord, hautement significatives (1%) pour (TB) et significatives (5%) pour (SO) en exposition Sud (tableau 20).

Les individus de chêne-liège abondent dans les expositions Nord avec un effectif moyen de 535 individus/ha à Souiniet (SO) et 523 individus/ha à Tbeinia (TB).

Le minimum est enregistré dans les placettes à exposition Sud, avec un effectif moyen de 291 individus/ha dans la station de Souiniet et 275 individus/ha pour la station de Tbeinia (figure 18).

D'une manière générale, les fortes densités sont rencontrées dans les stations exposées au Nord et de hautes altitudes, alors que, les densités les plus faibles des peuplements de chêne-liège sont localisées dans les terrains à exposition Sud et à basses altitudes.

4.1.3.3. Répartition de la densité en fonction de la pente, toutes classes de diamètres et d'altitudes confondues, en fonction de l'exposition.

La pente relativement importante, par sa raideur, détermine la stabilité du sol et la rétention de l'eau. Elle facilite le ruissellement de l'eau, avec amoindrissement des sols, et accumulation des dépôts dans les fonds et peut jouer un rôle déterminant dans la pédogenèse (Akrimi (2000)).

Les résultats des analyses de variance de l'effet de la pente sur la densité du chêne-liège en expositions Nord sont consignés dans le tableau 21.

Tableau 21 : Résultats de l'analyse de variance des effets de la pente sur la densité du chêne-liège (toutes classes de diamètres et d'altitudes confondues) en exposition Nord

	Station	Fobs	P>F	Sign
Pente faible (P1)	H I	3,08	0,05	*
	AB	3,53	0,01	**
	M II	1,72	0,16	NS
	M IV	3,88	0,01	**
	AS	3,85	0,01	**
	AM I	5,66	0,001	***
	M V	1,62	0,15	NS
	TB	5,87	0,001	***
	SO	6,54	0,001	***
	AZ III	1,10	0,32	NS
Pente moyenne (P2)	H I	2,63	0,05	*
	AB	2,85	0,05	*
	M II	0,98	0,48	NS
	M IV	2,89	0,05	*
	AS	2,50	0,05	*
	AM I	3,28	0,05	*
	M V	1,05	0,37	NS
	TB	4,26	0,001	***
	SO	5,17	0,001	***
	AZ III	0,78	0,65	NS
Pente forte (P3)	H I	2,21	0,05	*
	AB	2,58	0,05	*
	M II	0,15	0,92	NS
	M IV	2,55	0,05	*
	AS	2,09	0,10	NS
	AM I	2,75	0,05	*
	M V	0,71	0,19	NS
	TB	2,89	0,05	*
	SO	3,12	0,05	*
	AZ III	0,51	0,67	NS
Pente très forte (P4)	H I	1,54	0,36	NS
	AB	2,04	0,10	NS
	M II	0,42	0,73	NS
	M IV	1,88	0,14	NS
	AS	1,91	0,12	NS
	AM I	2,49	0,05	*
	M V	0,17	0,91	NS
	TB	2,55	0,05	*
	SO	2,35	0,05	*
	AZ III	0,15	0,92	NS

La lecture de la figure 19 nous permet de constater, d'une manière générale la variation de la densité du chêne-liège en fonction des classes de pentes.

Les effectifs élevés sont distribués sur les terrains à pentes faibles des stations de Tbeinia (TB), Souiniet (SO), Amdoun I (AM I) et Mekna IV (M IV), respectivement 582, 566, 345 et 426 individus/ha. Par contre, les densités les plus faibles sont enregistrées au sein des endroits à pentes fortes des stations de Hamdia I (H I), Ain Bacouch (AB) et Ain Saida (AS), respectivement 224, 241 et 312 individus/ha. Dans les terrains à pentes moyennes (P2), la densité des individus de chêne-liège tourne autour de 198/ha pour la station de Hamdia I et 512/ha pour la station de Tbeinia.

Concernant la densité de la troisième frange de pente (P3), elle reste assez bien représentée dans les stations d'AM I, SO et TB, respectivement 210, 235 et 280 individus/ha.

Les densités les plus faibles sont rencontrées dans les terrains à pentes très fortes, soient 175 individus/ha (SO), 156 individus/ha (TB) et 125 individus/ha (AM I).

Nos résultats montrent que le nombre d'individus/ha de chêne-liège par station est généralement plus important en terrains à pentes faibles ou modérés qu'en terrains à pentes fortes et très fortes pour une exposition Nord et à toutes altitudes.

Les résultats des analyses de variance de l'effet de la densité du chêne-liège, toutes franges d'altitudes confondues et à expositions Sud sont consignés dans le tableau 22.

Tableau 22 : Résultats de l'analyse de variance des effets de la pente sur la densité du chêne-liège (toutes classes de diamètres et d'altitudes confondues) et à exposition Sud

Classes de pentes	Station	Fobs	P>F	Sign
Pente faible (P1)	H I	2,58	0,05	*
	AB	3,10	0,05	*
	M II	0,95	0,45	NS
	M IV	3,28	0,05	*
	AS	2,85	0,05	*
	AM I	3,66	0,01	**
	M V	1,38	0,24	NS
	TB	3,11	0,05	*
	SO	3,74	0,01	**
	AZ III	1,20	0,31	NS
Pente moyenne (P2)	H I	1,20	0,31	NS
	AB	2,25	0,05	*
	M II	0,89	0,44	NS
	M IV	2,67	0,05	*
	AS	1,32	0,27	NS
	AM I	2,28	0,05	*
	M V	0,96	0,45	NS
	TB	3,55	0,01	**
	SO	3,17	0,05	*
	AZ III	1,12	0,33	NS
Pente forte (P3)	H I	2,45	0,05	*
	AB	2,03	0,10	NS
	M II	0,93	0,39	NS
	M IV	1,99	0,11	NS
	AS	2,09	0,10	NS
	AM I	2,25	0,05	*
	M V	0,62	0,59	NS
	TB	2,55	0,05	*
	SO	2,62	0,05	*
	AZ III	0,22	0,88	NS
Pente très forte (P4)	H I	0,42	0,73	NS
	AB	1,99	0,11	NS
	M II	0,28	0,84	NS
	M IV	0,91	0,40	NS
	AS	1,01	0,36	NS
	AM I	0,89	0,43	NS
	M V	0,47	0,59	NS
	TB	1,95	0,12	NS
	SO	0,79	0,50	NS
	AZ III	1,58	0,19	NS

La figure 20 nous montre que, en exposition Sud, le nombre d'individus par hectare du chêne-liège diminue avec l'accroissement du degré d'inclinaison des terrains, les densités maximales sont enregistrées dans les lieux à pentes faibles ou modérées.

Pour une pente faible, la densité des individus varie d'une station à une autre, les effectifs les plus élevés sont enregistrés au niveau des stations de Souiniet (SO) et Tbeinia (TB) (321 et 341 individus/ha), tandis que ceux les faibles sont rencontrés au sein des stations de Hamdia I (H I) et Ain Bacouch (AB) (160 et 125 individus/ha).

A partir de la même figure, on constate que les placettes à pentes modérées sont caractérisées par des effectifs moyens, soit 95 individus/ha (AB), 156 individus/ha (Mekna IV), 305 individus/ha (TB) et 287 individus/ha (SO).

La densité minimale du chêne-liège est située au niveau des terrains de fortes pentes, soit 88 individus/ha (H I), 93 individus/ha (AM I), 231 individus/ha (TB) et 209 individus/ha (SO).

Les résultats de l'effet des pentes très fortes sur la densité du chêne-liège sont non significatifs pour la totalité des stations d'études.

D'une manière générale, les résultats de l'effet de la pente sur la densité des peuplements du chêne-liège montrent que les effectifs les plus élevés sont enregistrés dans les endroits à pentes faibles ou moyennes, alors que, les minimales sont enregistrées dans les terrains à pentes fortes ou très fortes, quelques soient la nature de l'exposition et l'altitude.

4.1.3.4. Répartition de la densité du chêne-liège, toutes classes de diamètres confondues, en fonction des distances aux douars les plus proches

En Kroumirie, la pression humaine a entraîné le surpâturage dans les écosystèmes forestiers. Actuellement, les parcours sont, en général, incapables de répondre aux besoins nutritionnels du cheptel.

Le tableau 23 porte les résultats des analyses statistiques de l'effet des distances des douars aux peuplements sur les densités de ces derniers.

Tableau 23 : Résultats de l'analyse de variance des effets de la distance des douars les plus proches sur la densité du chêne-liège (toutes classes de diamètres confondues)

Station	Fobs	P>F	Sign
H I	5,54	0,01	***
AB	3,79	0,01	**
M II	2,53	0,05	*
M IV	4,26	0,001	***
AS	3,11	0,05	*
AM I	3,99	0,01	**
M V	5,17	0,001	***
AZ III	3,63	0,01	**
TB	5,16	0,001	***
SO	3,55	0,01	**

La corrélation entre ces deux variables est très hautement significative (P < 0,001) pour les stations de H I, M IV, M V et TB, hautement significative (P < 0,01) pour AB, AM I, AZ III et SO et significative (P < 0,05) pour M II et AS (tableau 23).

La figure 21 représente la relation entre la distance des peuplements des douars les plus proches et la densité de chêne-liège au niveau de chaque station.

Les densités les élevées sont celles comprises entre 660 et 665 individus/ha au sein des placettes les plus éloignées des douars des deux stations de Souiniet (SO) et Ain Zena III

(AZ III). Par contre, les densités des peuplements les plus faibles sont rencontrées au niveau des placettes les plus proches des douars pour ces deux stations, respectivement 245 individus/ha (Ain Zena III) et 256 individus/ha (Souiniet).

La lecture de la même figure montre que les densités les plus faibles sont enregistrées au niveau des stations de Hamdia I (H I), Ain Bacouch (AB) et Mekna II (M II), respectivement 243, 304 et 264 individu/ha.

On constate ici que l'action anthropique est très importante à cause de l'extension des parcours et l'évolution de l'effectif du cheptel. Ces trois stations sont les plus fréquentées par le cheptel, avec présence de traces de délits de coupes d'arbres et de maquis.

La position géographique des placettes par rapport aux douars influence donc, sérieusement, la densité du peuplement ; plus la distance est grande, plus la densité est forte.

4.2. Dynamique de l'écosystème chêne-liège en Kroumirie

Introduction

La dynamique de peuplements forestiers est liée à l'état pédologique (structure, PH, humidité, nature du sol, etc.), climatiques (température, pluviométrie), orographiques (altitude, pente et exposition), anthropiques (Homme, animaux domestiques et sauvages).

Selon Hasnaoui (2004), les facteurs de perturbations entraînent des modifications profondes du micro - climat au niveau de la trouée (rayonnement important, hygrométrie plus faible, etc.).

Cette nouvelle situation provoque l'installation d'une nouvelle dynamique de végétation.

Le manque de données écologiques pour les périodes 1995 et 2005, ainsi que la grande diversité entre les données dendrométriques actuelles et celles existant dans les deux inventaires, rendent l'établissement d'un modèle d'évolution écologique de l'écosystème de chêne-liège en Kroumirie (évolution de chaque classe de diamètre en fonction des facteurs de perturbations et en fonction du temps) difficile.

Le présent travail nous a permis de dresser l'évolution de la subéraie en Kroumirie durant la période 1983 - 2008, en nous basant sur des données dendrométriques et écologiques.

4.2.1. Dynamique de la subéraie de la Kroumirie en fonction des facteurs du milieu

Dans cette partie, nous avons dressé l'évolution de la densité et de la structure diamétrique de la subéraie en fonction des facteurs du milieu, en utilisant les moyennes de chaque classe de diamètre trouvées à partir des données de Hasnaoui (1992) comme base de référence pour l'état 1983 et de nos données (2008).

Les données des Inventaires ne sont pas prises en considération à cause de l'absence de la classe D1 (régénération) et l'absence des données concernant les facteurs du milieu.

Nous avons déterminé le degré d'évolution (régression ou progression) et le rythme moyen annuel de croissance de chaque classe de diamètres par la méthode des inventaires périodiques.

4.2.1.1. Evolution de la densité et de la structure diamétrique de la subéraie en fonction de l'altitude.

Les résultats des analyses de variance de l'effet de l'altitude sur la densité et la structure diamétrique (1983) sont consignés dans le tableau 24.

Tableau 24 : Résultats de l'analyse de la variance des effets de l'altitude sur la densité et la structure diamétrique (1983)

D'après le tableau 24, il existe un effet très hautement significatif (P<0,001) pour les classes de diamètres (D1 et D2) et un autre significatif (P<0,05) pour les arbres de classes de diamètres D3. Les résultats de variance sont significatifs pour les autres classes D4 et D5.

Pour étudier l'évolution de la densité et la structure diamétrique de la subéraie en fonction des classes d'altitudes, nous avons effectués une étude comparative des moyennes entre les différentes classes de diamètres durant la période 1983 et 2008, les résultats de l'évolution sont rapportés par le tableau 25 et la figure 22.

Tableau 25: Résultats de l'évolution de la densité et de la structure diamétrique

	Classe de diamètre	1983	2008	Différence	Rythme moy / an
	Classe de diamètre	N / ha	N / ha	N / ha	N / ha
Basse altitude	D1	453	75	- 378	- 15,75
	D2	09	05	- 4	- 0,16
	D3	07	03	-4	- 0,16
Moyenne altitude	D1	525	160	- 365	- 15,20
	D2	15	13	- 2	- 0,08
	D3	16	16	0	0
Haute altitude	D1	650	295	- 355	- 14,79
	D2	27	17	- 10	- 0,41
	D3	23	21	- 2	- 0,08

La lecture du tableau 25 montre que le stade semis (D1) a subi une régression importante avec les différentes classes d'altitudes, soit 378 individus/ha (basse altitude), 365 individus/ha (moyenne altitude) et 355 individus/ha (haute altitude).

Leur rythme moyen annuel varie de 15,75 individus/ha/an (basse altitude) à 15,20 individus/ha/an (moyenne altitude) et à 14,79 individus/ha/an (haute altitude).

La régression des classes D2 et D3 est très faible par rapport à la classe précédente, leurs rythmes moyens annuels de régressions varient avec les classes d'altitudes (respectivement 0,16 individus/ha/an (D2 et D3) à une basse altitude, 0,08 et 0 individus/ha/an (D2 et D3) à une moyenne altitude et 0,41 et 0,08 individus/ha/an (D2 et D3) à haute altitude).

La figure 22, représentant l'évolution de la densité de chêne-liège en fonction de l'altitude durant la période 1983 à 2008 montre:

- que la classe D1, correspondant aux semis ou régénération du chêne-liège est toujours mieux représentée que les classes D2 et D3, et ce quelles que soient la frange d'altitude et l'année d'observation.

- que l'effectif des semis rencontré en 1983 est, de loin, supérieur à celui comptabilisé en 2008.

- que le passage de la classe D1 à la classe D2 de diamètres est caractérisé par une chute brutale d'effectifs qui s'approchent de zéro et ce quelles que soient la frange d'altitude et l'année d'observation.

- que les effectifs des classes D2 et D3 sont très rapprochés en 1983 et 2008 et à toutes altitudes.

Il ressort de ces constats que la germination et l'installation des semis sont parfaitement possibles quoiqu'avec des effectifs hétérogènes d'une période à une autre et d'une frange d'altitude à une autre, avec un écart maximum entre les semis des stations alticoles (C3) de 1983 et ceux des stations de basses altitudes en 2008 (respectivement 650 individus/ha et 80 individus/ha).

Les classes D2 et D3 sont très peu représentées dénotant de difficultés de croissance rencontrés par les plantes et dont les causes sont multiples (Hasnaoui, 1992 et 1998).

Figure 22: Evolution de la densité et de la structure diamétrique de chêne-liège en fonction de l'altitude (C1, et C3)

4.2.1.2. Evolution de la densité et la structure diamétrique de la subéraie en fonction de l'exposition.

Les résultats de l'analyse de variance de l'effet de l'exposition sur la densité et la structure diamétrique (1983) sont portés sur le tableau 26.

Tableau 26 : Résultats de l'analyse de variance des effets de l'exposition sur la densité et la structure diamétrique (en 1983)

Selon le tableau 26, il existe une relation très hautement significative (P<0,001) pour les deux clases de diamètres D1 et D2. Les résultats des analyses de variance sont non significatifs pour les classes D3, D4 et D5.

L'évolution de la densité et de la structure diamétrique de chêne-liège durant la période 1983 et 2008 en fonction de l'exposition sont consignées dans le tableau 27 et la figure 23.

Tableau 27 : Résultats de l'évolution de la densité et de la structure diamétrique

du chêne-liège en fonction de l'exposition durant la période (1983-2008)

Durant cette période et pour une même exposition, la densité des semis a subi une régression, de 400 individus/ha, avec un rythme moyen annuel de régression 16,68 individus/ha/an dans les expositions (Nord).

Dans les stations chaudes (Sud), la régression est de 442 individus/ha, leur rythme annuel de régression est 18,41 individus/ha/an.

L'action combinée de la chaleur et de la sécheresse estivale d'une part, et de la prédation d'autre part, sont responsables de la régression constatée lors du passage du stade semis (D1) au stade gaulis (D2).

En revanche, les effectifs des individus de la classe D2 apparaissent très stables au fil du temps et en relation avec l'exposition, leur rythme moyen annuel de régression est de 0,08 individus/ha/an (Nord) et de 0,16 individus/ha/an (Sud).

En ce qui concerne l'évolution des autres classes de diamètres (D3, D4 et D5), les résultats des analyses de variance sont non significatifs pour les deux dates.

Pour une même frange d'exposition et pour la même année d'observation, la lecture de la figure 23 montre que l'évolution de la densité du chêne-liège, durant la période 1983 - 2008 des stades semis (D1) est mieux représentée que les stades gaulis (D2).

En passant d'une classe de diamètres à l'autre, l'effectif des semis (D1) connaît une chute brutale, et ce quelles que soient l'exposition et l'année d'observation.

Les effectifs des classes D2 sont très rapprochés durant cette période de vingt quatre ans et à toutes expositions.

Figure 23: Evolution de la densité et de la structure diamétrique du chêne-liège en fonction de l'exposition

4.2.1.3. Evolution de la densité et la structure diamétrique de la subéraie en fonction de la pente.

Les résultats de l'analyse de variance de l'effet de la pente sur la densité et la structure diamétrique de chêne-liège (1983) figurent dans le tableau 28.

Tableau 28 : Résultats de l'analyse de variance des effets de la pente sur la densité et la structure diamétrique (en 1983)

D'après le tableau 28, il existe une relation très hautement significative pour D1 et une autre significative pour D2.

Au terme de l'évolution de la densité et la structure diamétrique du chêne-liège en fonction de la pente, les principaux résultats d'évolution sont donnés par le tableau 29 et figure 24.

Tableau 29 : Résultats de l'évolution de la densité du chêne-liège et la structure diamétrique en fonction des classes de pente durant la période (1983-2008)

	Classe de diamètre	1983	2008	Différence	Rythme moy / an
	Classe de diamètre	N / ha	N / ha	N / ha	N / ha
Faible pente (P1)	D1	585	297	- 288	- 12
Faible pente (P1)	D2	28	23	- 5	- 0,20
Moyenne pente (P2)	D1	475	178	- 297	- 12,37
Moyenne pente (P2)	D2	23	16	- 7	- 0,29
Forte pente (P3)	D1	325	18	- 307	- 12,79
Forte pente (P3)	D2	15	8	- 7	- 0,29
Très forte pente (P4)	D1	0	0	0	0
	D2	9	0	- 9	- 0,37

D'après le tableau 29, les faibles différences d'effectifs sont observées dans les terrains à pentes faibles (P1) ou modérées (P2) pour le stade des semis, respectivement 288 et 297 individus/ha ; leur rythme moyen annuel de régression sont de 12 et 12,37 individus/ha/an.

Les écarts maximums sont rencontrés sur les terrains à fortes pentes (307 individus/ha), avec un rythme annuel de régression de 12,79 individus/ha/an.

En effet, on a une absence totale de l'installation des semis sur les terrains à pentes très fortes (P4).

En ce qui concerne les stades gaulis (D2), les fortes régressions sont constatées dans les endroits à pentes très fortes (9 individus/ha), leur rythme moyen annuel de régression est de 0,37 individus/ha/an.

La pente n'a pas d'effet significatif sur la classe de diamètres pour D3, D4 et D5.

La figure 24 montre que les densités maximales sont mieux représentées pour la classe de semis (D1) que la classe de gaulis (D2).

On constate une chute brutale des effectifs en passant de la classe de semis à celle de gaulis en relation avec les classes de pentes et l'année d'observation.

En consultant les mêmes résultats de cette figure, on ne remarque pas une grande différence pour les densités du gaulis (D2) durant cette période et toutes franges de pentes confondues.

Figure 24 : Evolution de la densité et de la structure diamétrique de chêne-liège en fonction de la pente

4.2.1.4. Evolution de la densité et de la structure diamétrique de la subéraie en fonction des distances aux douars les plus proches.

Les résultats de l'analyse de variance de l'effet de la distance des peuplements par rapport aux

douars les plus proches sur la densité et la structure diamétrique (1983) sont consignés dans le

Tableau 30 : Résultats de l'analyse de variance des effets de la distance des peuplements les plus proches par rapport aux douars sur la répartition des classes de diamètres (1983)

Selon le tableau 30, il existe une relation très hautement significative pour D1et D2 et une autre significative pour D3 et D5.

Les résultats de l'évolution de la densité et de la structure diamétrique du chêne-liège en fonction des distances aux douars les plus proches durant la période 1983-2008 sont portés sur le tableau 31 et figure 25.

Tableau 31 : Résultats de l'évolution de la densité et la structure diamétrique

du chêne-liège en fonction de la distance des douars les plus proches durant la période (1983-2008)

	Classe de diamètre	1983	2008	Différence	Rythme moy / an
	Classe de diamètre	N / ha	N / ha	N / ha	N / ha
Douars très proches	D1	480	25	- 455	- 18,95
	D2	11	06	- 6	- 0,16
	D3	14	05	-9	- 0,37
	D5	30	12	- 18	- 0,75
Assez éloignés	D1	525	130	- 395	- 16,45
	D2	13	10	- 2	- 0,12
	D3	13	08	- 5	- 0,20
	D5	45	33	- 10	- 0,41
Très éloignés	D1	595	299	- 296	- 12,33
	D2	22	20	- 2	- 0,08
	D3	21	18	- 3	- 0,12
	D5	78	75	- 3	- 0,12

Selon le tableau 31 que ce soit en 1983 ou en 2008, les effectifs des jeunes semis varient d'une classe de distance à une autre et d'une période à une autre, il apparaît que les peuplements les plus proches et assez éloignés des douars ont subi de fortes régressions (respectivement 455 et 395 individus/ha).

Par ailleurs, les fortes densités sont enregistrées dans les lieux les plus éloignés des douars, avec un écart maximum de 296 individus/ha entre 1983 et 2008. Le rythme moyen annuel de régression des jeunes plantules diminue en s'éloignant des douars les plus proches, soient 18,95 individus/ha/an (douars plus proches), 16,45 individus/ha/an (douars assez éloignés) et 12,33 individus/ha/an (douars plus éloignés).

En revanche, les effectifs des autres classes de diamètres D2, D3 et D5 présentent certaines ressemblances en passant d'une période à une autre, avec une différence maximale entre les semenciers (D5) des placettes les plus éloignées (T3), assez éloignées (T2) et les plus proches des douars (T1) (respectivement de 18, 10 et 3 individus/ha).

Le rythme moyen annuel des semenciers (D5) diminue en effectifs en s'éloignant des douars, soient 0,75 individus/ha/an (douars plus proches), 0,41 individus/ha/an (douars assez éloignés) et 0,12 individus/ha/an (douars très éloignés). Le rythme moyen annuel de régression des classes D2 et D3 déminue en passant des endroits les plus proches des douars (T1) vers les endroits les plus éloignés (T3) (respectivement 0,16 et 0, 37 individus/ha/an (T1) et 0,08 et 0,12 individus/ha/an (T3)).

Cette forte régression des différentes classes de diamètres est liée principalement à l'évolution des effectifs du cheptel (ovin, bovin et caprin), l'explosion démographique, l'évolution des pratiques délictuelles (coupes de bois de chauffage, carbonisation, terres agricoles, etc.).

Figure 25 : Evolution de la densité et de la structure diamétrique de chêne-liège en fonction des distances aux douars

La lecture de la figure 25 montre la variation de l'évolution des peuplements de chêne-liège en Kroumirie en fonction des distances aux douars les proches.

Ces résultats mettent en exigence l'importance de la fréquence des semis quelles que soient la distance aux douars les plus proches et l'année d'observation.

Le passage de la classe D1 à la classe D2 de diamètre est caractérisé comme pour les autres paramètres par une chute brutale d'effectifs en fonction des distances par rapport aux douars les plus proches et l'année d'observation.

Les densités des autres classes de diamètres (D2, D3 et D5) sont très rapprochées en 1983 et 2008 et à toutes classes de distances.

4.2.2. Etude comparative entre la structure actuelle et les structures antérieures

Les différentes informations acquises sur la structure des peuplements de la subéraie en Kroumirie depuis les années 1980, permettent d'avoir plusieurs repères de comparaison sur les états successifs de ces peuplements.

D'après les résultats des deux inventaires nationaux 1995 et 2005, la structure diamétrique de la subéraie en Kroumirie est régulière.

Elle définit un aspect plus ou moins normal avec une dissymétrie à gauche caractéristique des peuplements pas trop âgés. C'est une forme de courbe gaussienne (figures 26 et 27).

Le dépouillement des données similaires relatives à la subéraie des données de Hasnaoui (1983) a montré que la distribution diamétrique est représentée par des taux élevés des semis (D1) et des semenciers (D5). Les classes intermédiaires sont les moins représentées (figure 28).

Figure 28 : Structure diamétrique globale de la subéraie en Kroumirie (1983)

L'étude comparative entre les différentes structures diamétriques des peuplements de la subéraie, montre que les résultats de 2008 sont semblables à ceux de 1983. Par contre, ils sont divergents avec ceux figurant dans les deux inventaires 1995 et 2005 (figure 29).

Figure 29 : Structures diamétriques de la subéraie en Kroumirie durant les années 1983, 1995, 2005 et 2008

Chapitre 5 : Discussions

Le présent travail nous a permis de représenter la structure actuelle de la subéraie en Kroumirie, de dégager les principaux facteurs qui influencent la répartition des individus de peuplements purs de chêne-liège en classes de diamètres et de donner une idée sur l'évolution de chaque classe de diamètres durant la période (1983-2008).

L'étude de la structure diamétrique actuelle des peuplements de chêne-liège en Kroumirie nous permet de conclure à une grande diversité au sein d'une même station et d'une station à une autre.

Selon les résultats obtenus, la distribution diamétrique actuelle du peuplement de la subéraie est marquée par une très forte proportion dans les classes de diamètres D1 soient 65,78 %.

Par contre, elle est définie par des densités faibles au niveau des classes D2, D3 et D4, respectivement 4 %, 5 % et 6,80 %. La classe D5, représente, quant à elle, 18,42 %.

L'étude comparative entre la structure actuelle de la subéraie en Kroumirie par rapport aux autres structures définies par les deux inventaires périodiques (1995 et 2005) et les données de Hasnaoui, 1983, montre que la structure des peuplements (1983) est semblable à celle du peuplement actuel. Les deux inventaires 1995 et 2005 présentent des structures différentes.

Les résultats du premier inventaire (1995) montrent que, la distribution diamétrique est marquée par une très forte densité dans les classes de diamètres D3 et D4, soient 44,07 % et 24,50 %, les deux classes D2 et D5 étant moins représentées, soient 14,82 % et 16,60 %.

L'étude du deuxième inventaire (2005) montre que, la distribution diamétrique du peuplement de la subéraie est définie par des effectifs élevés pour la classe de diamètres D3 (59,12 % de l'ensemble) et assez faibles pour D4 (16,78 %). Les classes D2 et D5 présentent, respectivement, des taux d'effectifs de 13,87 % et 10,21 %. Dans les deux inventaires, les jeunes semis, D1, ne sont pas représentés.

Nos résultats conformes à ceux trouvés par Hasnaoui (1992) qui a souligné que la subéraie tunisienne est caractérisée par deux strates ; une strate arborescente comprenant les individus adultes ou semenciers et la strate des plantules. Entre les deux on ne rencontre que de rares sujets appartenant aux stades intermédiaires (fourrées, gaulis, perchis).

Ce même auteur a montré que l'étude de la structure diamétrique des populations des différents chênes présents en Kroumirie, nous permet de conclure que celle-ci est loin d'être satisfaisante car elle montre, en particulier, une rupture nette de la croissance des jeunes plantules issues de semis, ce qui permet de prévoir un avenir plutôt inquiétant de ces chênaies. Les causes de cette situation peuvent être multiples : propres à l'essence étudiée dans certains milieux, liées à la nature et à la composition du groupement végétal, et enfin, marquées par l'impact humain direct ou indirect.

La variation de la structure diamétrique au sein d'une même station ou d'une station à une autre est le résultat de l'interaction des facteurs d'ordres orographiques (pente, altitude, exposition) et anthropiques (distance par rapport aux douars les proches).

Les résultats des analyses statistiques de l'effet de l'altitude sur la densité et la structure diamétrique montrent que les zones alticoles (supérieurs à 250 m) contiennent des effectifs des jeunes plantules et des semenciers plus importantes que dans les sites de basses altitudes.

De même, la densité des différentes classes de diamètres varie aussi avec l'exposition, nos résultats montrent que les versants Nord (humides) sont les plus denses.

Sur les versants Sud, où la chaleur et la sécheresse sont plus accentuées, la densité des semis et des arbres de chêne-liège est faible.

Selon les analyses statistiques de l'effet de la pente sur la densité des différentes classes de diamètres, les fortes densités sont enregistrées sur les terrains de faibles ou moyennes pentes.

Les pentes très fortes ne permettent l'installation que d'un effectif très faible de semis, parce que les glands relativement lourds, sous l'action de la pesanteur et de l'eau, sont charriés vers le bas où ils vont trouver des conditions meilleures de sol et d'humidité pour germer.

Nos résultats confirment ceux trouvés par Khanfouci (2005), l'altitude rend le milieu plus humide ou plus sec par une augmentation ou diminution des précipitations et des températures. Selon ce même auteur, les zones alticoles assurent les meilleures conditions favorables à la germination et le développement des plantes.

Selon Block et Treter (2001), les facteurs environnementaux influent sur la régénération des plantes et sur leurs structures diamétriques.

Nos résultats confirment aussi ceux trouvés par d'autres auteurs : Nsibi et al 2006, Ezzahari et al 2000, Ben Brahim et al 2004, Durand et al 2002, Hasnaoui (1998), Arnold et al 1989 et Khanfouci (2005).

Nsibi et al (2006) ont montré qu'un versant orienté vers le sud reçoit, à surface égale, plus de radiation qu'un versant exposé au nord. Par suite de cette forte insolation, la lumière est plus intense, la chaleur est plus grande, l'évaporation plus considérable et l'air plus sec. Inversement, le versant exposé au nord est moins éclairé, moins chaud et l'air y est plus humide.

Selon Ben Brahim et al (2004), l'exposition et la forme des versants ont également une influence importante : à une altitude donnée les caractéristiques climatiques peuvent être très différentes entre un versant exposé vers le nord et un autre exposé vers le sud.

Ezzahiri et al (2000) ont montré que les stations exposées à l'Ouest et au Nord-ouest offrent des conditions favorables à la germination des graines, les expositions Sud sont généralement les plus chaudes et la germination des graines est précoce. Dans cette situation, la chaleur estivale, trop élevée, fait disparaître les jeunes plantules les plus exposées au soleil.

Messaoudène et al (1998) ont montré que les expositions chaudes et ensoleillées et à basses altitudes apparaissent plus favorables à l'installation des semis de chêne-liège que les stations à expositions froides et à hautes altitudes.

Hasnaoui (1998) a noté que les pentes fortes sont défavorables à l'installation des glands et même s'ils arrivent à germer il leur sera difficile de survivre à cause du manque de la rétention de l'eau par le sol, et à cause de la pauvreté de ce sol. Dans les bas fonds, le sol est généralement riche, profond, meuble et humide.

D'après Durand et al (2002), les terrains à drainage vertical se montrent plus favorables que les terrains à drainage latéral superficiel pour la répartition des classes diamètres.

Selon Arnold et al (1989), plus la pente est rapide, plus l'eau érodera le sol. L'érosion hydrique augmente aussi avec la longueur de la pente à cause de l'augmentation du ruissellement. Les répercussions de l'érosion des sols vont au-delà de la perte du sol de la couche arable. La levée des plantes, leur croissance et leur rendement sont directement affectés par la perte d'éléments nutritifs du sol et des engrais éventuels qui y sont appliqués. La semence et les plantes peuvent être dérangées ou totalement transportées hors du lieu où se produit l'érosion. Le sol érodé, déposé au bas des pentes, peut empêcher ou retarder l'émergence de la semence, enterrer les jeunes pousses.

Nsibi et al (2006) ont montré que les fortes pentes ont une influence marquée sur la répartition des semis à cause de l'érosion active. Les glands, en tombant ne trouvent pas les conditions favorables pour germer et s'exposent aux risques de destruction (déprédation et desséchement). En revanche, dans les classes à pentes faibles, les semis naturels de chêne-liège bénéficient d'un supplément d'eau pluviale et d'élément nutritifs provenant du ruissellement de l'amont.

Ainsi, Khanfouci (2005), a indiqué que les terrains en pente, exposés vers le Sud et à basses altitudes sont les stations les moins arrosées et qui sont extrêmement défavorables au maintien du développement des jeunes plantules.

En effet, les stations à expositions Nord, à pentes modérées ou faibles et à hautes altitudes sont les stations les plus fertiles, plus denses et les plus arrosées.

L'action des distances des peuplements par rapport aux douars les plus proches (facteurs anthropiques), a un effet très grave sur la structure de la subéraie et sa stabilité.

En effet, la densité des semis est très élevée au niveau des stations les plus éloignées des douars.

Les analyses statistiques (analyse de variance) de l'effet des distances des douars les plus proches sur la densité des peuplements de chêne-liège (toutes classes de diamètres confondues), pour chaque station, montrent que les effectifs élevés sont enregistrés sur les placettes les plus éloignées des douars.

D'après ces résultats, la densité régresse en se rapprochant des douars, Or les placettes proches des habitations sont très fréquentées par les animaux domestiques (bovins, ovins, caprins, etc.).

L'action des facteurs anthropiques est très forte dans les stations de Hamdia I, Mekna II et Ain Bacouch. Ces trois stations sont caractérisées par une grande extension des parcours à cause des effectifs élevés du cheptel.

Au niveau des stations de Tbeinia, Souiniet et Ain Zena III, nous avons souligné que l'intensité des coupes de bois de chauffage et de carbonisation surtout dans les peuplements les plus proches.

La pauvreté des habitants des clairières forestières et des zones limitant des forêts a un impact direct sur la structure de la subéraie. En effet, la satisfaction de leurs besoins nécessite l'extension de leur terre agricole au déterminent de la végétation naturelle (défrichement) et entraine une pression de plus en plus forte sur l'écosystème forestier en pratiquant des prélèvements de bois, de glands, de feuillage et en coupant les arbres pour la carbonisation. Tous ces facteurs rendent les possibilités de sauvegarde et de reconstitution de la subéraie en Kroumirie très difficile.

Des résultats similaires furent rapportés par Hasnaoui (1998), Beltrain (2002), Boussaidi (2005), Ben Jemâa et al (2006), Bendaanoun (1998), Benchekroun (1998), Yacoubi (2000) et Ben Brahim et al 2004 .

Selon Hasnaoui (1998), l'homme agit par l'intermédiaire de ses animaux domestiques ou en perturbant le milieu forestier par la coupe de maquis et même d'arbres pour ses différents usages ; bois de chauffage, carbonisation, confection de clôture.

Beltrain (2002) a démontré que les ébranchages ont été abusifs. Les «élagages excessifs » sont néfastes pour le développement de l'arbre, l'affaiblissant et le rendant plus sensible à la sécheresse, et aux maladies. De plus, la forte éclaircie des houppiers entraîne une invasion du sous étage par des espèces héliophiles qui dégradent la subéraie.

Pour Bendaanoun (1998), les troupeaux, par leur action sélective, ont tendance à éliminer la flore la plus palatable au profit des taxons toxiques et très peu appétés, ce qui constitue un déséquilibre au niveau des structures de la végétation à chêne-liège et un frein pour la régénération.

Selon Boussaidi (2005), à long terme le surpâturage provoque la dégradation de la pâture mais à court terme, un pâturage excessif diminue la production en réduisant la repousse.

Ben Brahim et al 2004 confirme que l'action anthropique a ouvert les massifs montagneux et collinaires rifains et pré rifains à la dégradation de l'environnement en délimitant les espaces forestiers.

D'après Ben Jemâa et al (2006), la poussée démographiques et sans doute le principal facteur de la dégradation de la subéraie en Kroumirie. Les densités des populations sont parmi les plus élevées des pays : 100 habitants au km² et par fois plus.

Selon Yacoubi (2000), la poussée démographique et sans doute le principal facteur de dégradation de la subéraie en Kroumirie. Les densités des populations sont très élevées : 96 habitants au Km² en 1994, leur ressource principale est l'exploitation du milieu environnant en pratiquant l'agriculture vivrière et l'arboriculture dans les limites de leurs clairière et en expliquant leur droit d'usage dans la forêt par le biais du pastoralisme et le prélèvement du bois (par différents usages) et d'autres menu produits.

Benchekroun (1998) a confirmé que l'explosion de la population et son besoin annuel en bois de feu, la surexploitation et le problème de disponibilité sont les causes de la régression des subéraies.

De même, Ezzahiri et al 2001 ont montré que les actions de la mutilation, l'absence de la régénération et le dessèchement localisé et même généralisé dans certains endroits sont probablement dus à des sécheresses prolongés et répétées.

La sécheresse, la chaleur, le froid, la sénescence et les défoliations répétées sont des facteurs qui affaiblissent les arbres et favorisent l'installation des insectes xylophages et des pathogènes conduisant au dépérissement des arbres et à leur mortalité (El Abidine (2003)).

D'après Hasnaoui (2008), dans la forêt de la Kroumirie la période de sécheresse estivale est assez longue (3 à 5 mois), la perte d'eau par ruissellement due à la texture légère est accentuée par le dommage du sol provoqué par le surpâturage.

L'affaiblissement des peuplements, évoluant vers un dépérissement grave au sein des arbres, résulte de la diminution des réserves en eau du sol suite à la sécheresse prolongée.

Pour la subéraie pure de la Kroumirie, le taux de dépérissement le plus fort est enregistré dans zones de basses altitudes (27 %). Dans les zones alticoles (à partir de 250 m d'altitude), le dépérissement est moins représenté et va en diminuant (taux de dépérissement moyen de 20 à 14 %). L'altitude a un effet sur l'état phytosanitaire des arbres (par le bais de l'humidité et de la température). Cet état est de plus en plus dégradé pour les peuplements à basse altitude.

Ce même auteur montre que, les placettes en pentes exposées au Sud (Sud, Sud - Est et Sud - Ouest) présentent des taux de dépérissement importants, respectivement 25, 26 et 31 %.

Les autres expositions (Nord - Ouest, Nord - Est, Nord, Ouest et Est) présentent des taux de dépérissement plus faibles.

Concernant l'effet de la pente, Hasnaoui (2008) montre que l'état sanitaire des arbres de chêne-liège s'affaiblit pour les terrains en pente forte (supérieur à 30 %) avec un taux de dépérissement moyen de 26 %. Cet état s'améliore lorsque la pente devient modérée à faible avec un taux de dépérissement moyen (inférieur à 20 %).

Conclusion générale

Au cours de ce travail, nous avons utilisé la structure diamétrique des peuplements purs de chêne-liège en Kroumirie comme indicateur de l'état actuel de ces peuplements et de leur évolution.

La comparaison des données de 1983 avec celles de 2008 montre que la structure diamétrique présente le même aspect, mais avec une baisse considérable, en valeur absolue, des effectifs des semis, ce qui dénote d'une pression de plus en plus forte, des facteurs du milieu, y compris les facteurs anthropiques, sur l'écosystème .

L'absence quasi-totale, des stades intermédiaires issues de semis, laisse présager le pire quant à l'avenir de nos subéraie.

Il est à souligner que nos résultats sont différents de ceux des inventaires nationaux forestiers de 1995 et 2005.

Le devenir de cet écosystème reste incertain si nous ne prenons pas les mesures nécessaires pour la sauvegarde.

Parmi ces mesures vient en premier la régénération naturelle par rejets de souches. Ce type de régénération, qui peut être entrepris dans le cadre d'une sylviculture appropriée du chêne-liège, peut donner des résultats satisfaisantes moyennement des précautions à prendre d'une manière stricte, comme la date de coupe, la mise en défens et la conduite des rejets (Hasnaoui, 1992).

Par ailleurs, et étant donnée l'infrastructure déjà existante en matière de pépinières modernes, il est impératif de procéder à des reboisements à grande échelle par les espèces autochtones, dont notamment le chêne-liège.

Ce type de reboisement est d'autant plus recommandé, que les superficies déboisées dans les régions humides et subhumides sont considérables.

Références bibliographiques

Abdendi Z. A., 2003. Forest decline in Morroco: Causes and control strateg. Science et changement planétaires/sécheresse. Volume 14 ; pp : 9-18.

Akrimi N., 2000. Cours d'écologie dispensé aux étudiants du cycle ingénieur « formation continue ».

Alder D., 1980. Estimation des volumes et accroissement des peuplements forestiers avec référence particulière aux forêts tropicales. Étude et prévisions de la production. Vol 2. pp 23.

Aloui A., 2007. Analyse et diagnostic de l'état actuel de la subéraie tunisienne et proposition d'une stratégie de développement durable. Etude stratégique pour le développement durable de la subéraie tunisienne. Rapport de mission de la consultation nationale en aménagements forestiers et sylvicoles. Pp : 12.

Aloui A., 1982. Recherches dendroclimatologiques en Kroumirie (Tunisie). Thèse de doctorat présentée l'université d'Aix-Marseille III. 203 p.

Amandier L., 2002. La subéraie, biodiversité et paysage. Colloque 2002. pp 16.

André P., Reginster P et Reginster J., 1983. Evolution théorique d'une pessière en fonction de l'intensité de l'éclaircie. Bull. Soc. R. For.Belg. pp 180-187.

Arnold, J. B, Wall G., Moore N., Baldwin C. S et Shelton I. J., 1989. L'érosion du sol - causes et effets. 89-064. I.S.S.N 1198 6 7138.

Ashton, L. S et Hall P., 1992. Comparaison of structure amony mixed dipterocarps forests of nother-western Borneo J Ecol. pp: 459-481.

Aussenac, G., 2000. Introductions betwen forest stands and microclimate: Ecophysiological aspects and consequences for silviculture. Ann.For.Sci.2000. pp : 287-301.

Bedel F., Durrieu L., Dupy B., Vincent F., Maître H., Borhen A et Narboni P., 1989. Dynamique de croissance dans des peuplements exploités et éclaircies de forêt dense africaine. Le dispositif de M'Baiki en république centre africaine (1982-1995). Pp : 32.

Bellkhdar J., 2003. Les Maghreb à travers ces plantes ; plantes, productions végétales et traditions au maghreb. Ethomobotanique, histoire et société : Les espèces Emblématiques de la tradition locale. Pp : 4.

Beltrain R. S., 2002. Subéraies, biodiversité et production de liège. Colloque 2002. pp : 10.

Benabid A., 2002. Le Rif et le Moyen atlas (Maroc) : Biodiversité, Menaces, Préservation, pp : 9.

Benabid A., 1998. Ecologie, conservation et restauration des subéraies. Formation continue à l'E.N.F.I. 20 p.

Ben Brahim K, Ismail M et Tribak A., 2004. Problèmes de dégradation de l'environnement par la désertification : impact du phénomène au Maroc. Science et changement planétaires /sécheresse volume 15, Numéro 4, 307 p.

Benchakroun F., 1998. La nécessité d'intégration des aspects socio-économiques dans les approches techniques de la régénération. Annales de l'I.N.G.R.E.F 1998.

Bendaanoum M., 1998. Contribution à l'étude des facteurs écologiques, de l'impact de la régénération des subéraie du rif, du Moyen Atlas oriental et de la Mamora. (Maroc septentritionnal) Actes du séminaire Méditerranéen sur la régénération des forêts de chêne-liège Tabarka. Annales de l'I.N.G.R.E.F 1998.

Ben Jamaa M. L, Abid H et Nouioui M., 2006. Impact de la gestion du liège sur l'économie rurale et les incendies des forêts en Tunisie. Collègue international vivexpo. L'homme et le liège 2006, pp : 1-10.

Ben Jamaa M.L., Hasnaoui B., 1996. Le dépérissement du chêne-liège (Quercus Suber L.) en Tunisie. In call Nat sur le dépérissement des forêts du Maroc, CNRF Rabat -Maroc 28 et 29 février 1996.

Ben Tiba B., 1980. Contribution pollenalytique à l'histoire Holocène de la végétation de la Kroumirie (Tunisie septentrionale). Thèse doctorat ingénieur univers. D'Aix-Marseille III, 76 p.

Ben Seghir L. A., 1996. Mastère en science forestière institut méditerranéen du liége.

Bergérs L., 2004. Rôles des coupes de la stratification verticale et du mode de traitement sur la biodiversité. In cosselin M, Laroussinie, O (Eds). Gestion Forestière et Biodiversité : connaître pour préserver synthèse bibliographique. CEMAGREF Edition, Antony, pp : 149-215.

Block J et Treter V., 2001. Les facteurs limitants à la partie supérieure et inférieure des forêts dans les limites de montagne -bois. 18 p.

Bortoli, L., Gounot M., et Jacquiot J CL., 1969. Climatologie et bioclimatologie de la Tunisie septentrionale. Annales de l'I.N.R.A. Tunisie Vol 42. Fasc. 1, pp : 9-104.

Boudru M., 1973. La futaie et ses traitements. Ann. Gembloux. 79(3), pp : 155-164.

Boudru M., 1989. Forêt et sylviculture traitement des forêts. Les presses agronomiques de Gembloux, A.S.B. pp : 15.

Boudy P., 1950. Economie forestière Nord-africaine I, II. Monographie et traitement des essences forestières. Ed. La rose. 525 p.

Boudy P., 1952. Guide forestier en Afrique du Nord. La Maison Rustique, 26, rue Jacob -Paris 6, 505 p.

Bourgenot L., 1968. En France. Nouvelles tendances en matière d'aménagement. Rev. For, Franç. 27 (3), pp : 179-348.

Bourgouin D., 1998. La subéraie Française état, réhabilitation et programmes de recherche. Séminaire méditerranéen sur la régénération des forêts de chêne-liège. Annales de l'I.N.G.R.E.F 1998.

Bourliere P, Lamotte M., 1978. La notion d'écosystème. Un problème d'écologie. Ecosystèmes terrestres, Ed, Masson. Paris, pp : 10.

Bousaidi N., 2005. Parcours en forêt et risque de dégradation des potentialités pastorales dans la IV éme série forestière de Mekna (Tabarka - Tunisie). Mémoire du troisième cycle universitaire Tunis - cartage (I.N.A.T), pp : 15.

Busby J.A., 1965. Sttudies on the stability of conifer stands. Scottish Forestry. 15 : 86-102.

Carles S., Lauriac A., Lecontre B., et Mariton B., 2001. Forêts privées des Albères et de la cote rocheuse orientations de gestion, orientations régionales de production. Schéma régional de gestion sylvicole. Tome 2, centre régionale de la propriété forestière langue doc - Roussilon. 32 p.

Christina L., Valerie M., 2001. Introduction and evolution of possible indices of stand structural diversity. Pp: 10.

Collinet F., 1997. Essai de regroupement des principales essences structurales d'une forêt dense humide d'après l'analyse de leur répartition spatiale (Forêt de Paracon - Guyane). Thèse de doctorat de l'université Claude Bernard Lyon I/ENGREF. 203 p + annexes.

Cantat R., Piazzetta R., 2005. Le levé du liège, ce qu'il faut savoir sur l'exploitation du chêne-liège, Guide technique et de vulgarisation. Institut méditerranéen du liège. 45 p.

Courbeau B., Houlier H., Rupe C., 1993. Un modèle de croissance en hauteur des arbres en pessière pure irrégulière de montagne. Ann. Sci. For (1993) 50, pp : 337-351. Elsevier/INRA.

Courbeau B., 1991. Introduction à une modélisation de la croissance en forêt irrégulière de montagne. DEA de biologie forestière, université Nancy I / E.N.G.R.E.F / C.E.M.A.G.R.E.F. 13 p.

Danniel D., 2008. L'effet du patron de répartition des coupes et des variabilités du milieu sur les pertes par chablis dans les lisières. Cas de la pessière à mousses de l'Est. Mémoire d'obtention du grade de maître en science. Faculté de foresterie et de géomatique université Laval Québec.

Dauglas A., Charles B., David L., 2007. Changes in forest structure following variable retention harvests in Dauglas-fidominoted forests. Pp: 9.

Débierre F., 1922. Le chêne-liège en Tunisie. Imp. Centrale. Tunis, 55 p.

De Lemps F., 1958. Qu'est-ce que la futaie par bouquets ? Rev. For. Franç : 10 (3), pp : 183-194.

Derout R., 1994. Méthode d'étude de la structure des populations en forêts tropicales humides. DEA de l'université Lyon I, 51 p.

Direco Geral dos Reursos Florestais, 2006. Séminaire vitale des peuplements de chênes-lièges et chênes verts : situation actuelle, état des connaissances et actions à entre prendre. (25-26 octobre 2006), pp : 33.

Direction Générales des forêts, 1995. Résultats du premier inventaire forestier et pastoral : Ministère d'agriculture ; Tunisie, pp : 9-14.

Direction Générale des Forêts, 2000 : Résultats de l'inventaire forestier National.

Direction Générale des forêts, 2005. Inventaire Forestier National ; Tunisie .88 p.

Durant J.M., Hosgood N et Toutchkov M., 2002. Etude de l'influence du sol et du substrat géologique sur la structure d'un peuplement forestier Guyanais, pp : 13-15.

Dupieux N., 2004. Une proposition de protocole commun pour la description et le suivi des forêts alluviales du bassin de la loire. Programme Loire nature, mission scientifique, 41 p.

Dyer M.E et Bailey R.L., 1987. A test of six methods of estimating tree heights from stem analysis data. Forest.Sci.p33, pp : 3-13.

El Afsa M., 1978. Ecologie, phytosociologie, régénération et production des subéraies tunisiennes Fa.Sci et Tech. ST Jérôme Aix - Marseille. III, pp : 15.

Emberger L., 1930. La végétation de la région méditerranéenne. Essaie d'une classification des groupements végétaux. Rev. Genbot 42, pp : 641-662.

Emberger L., 1954. Projet d'une classification biogéographique des climats, les divisions écologiques du monde C.N.R.S Paris.

Et toubi M., 1996. Contribution à l'étude de la dynamique et du dépérissement du chêne-liège au Mamora. Mémoire de 3éme cycle pour l'obtention de diplôme d'ingénieur des eaux et forêts Salé Maroc.

Ezzahiri M., Belghazi B et Ouchtarmoune M., 2001. Evaluation de la dynamique de cedraie du haut atlas oriental par les méthodes diachroniques. (Cas de la cedraie de Midkane). Intrenationnal Meeting on silviculture of cork oak (Quercus suber L.) and cedra (Cedrus atlantica (Endl) M). Rabat - Maroc, 2001, pp : 53.

Ezzahiri M., Belghazi B., 2000. Synthèse de quelques résultats sur la régénération naturelle du cèdre de l'Atlas au Moyen - Atlas (Maroc). Science et changement planétaires/sécheresse. Volume 11, Numéro 2, pp : 79-84, Juin 2000.

Ezzahiri M., Bellghazi B et Bahmed M., 1994. Bilan de la régénération naturelle de la cédraie dans les parcelles culture du Moyen - Atlas, Maroc. Ann. Rech. For, Maroc, pp 13-16.

Foundi. M., Ezzahiri M et Belghazi B., 2001. Evaluation de la dégradation des subéraies du rif central en vue d'aménagement adéquat (cas de la subéraie de Jbel Asri à Chechoouen). Intrenationnal Meeting on silviculture of cork oak (Quercus suber L.) and cedra (Cedrus atlantica (Endl) M). Rabat - Maroc, 2001, pp : 22.

Franclet R., 1972. Cours de sylviculture. Projet FAO (Education et formation forestière).

Gader G., Daly H., 2000. L'étude prospective du secteur forestier en Afrique (FOSA) Tunisie.

Gaudin S., 1996. Dendrométrie des peuplements. BTSA Gestion forestière Module D 42, pp : 42.

Goreaud F., De Coligny F., Courbeaud B., DHote J.F., Dreyfus P.H., Pérot T., 2005. La modélisation : un outil pour la gestion et l'aménagement en forêt. La revue en sciences de l'environnement. Volume 6 n : 2 Septembre 2005, 13 p.

Gounot M., 1969. Méthodes d'étude quantitative de la végétation. Masson, Paris, 314 p.

Hasnaoui B., 1991. Régénération naturelle par rejet de souche et par drageonnement d'une subéraie dans le Nord-Ouest de la Tunisie. XVII (Ecologia Méditerranaea, 79-87). Pp : 13-21.

Hasnaoui B., 1992. Chênaies du Nord de la Tunisie, Ecologie et régénération. Doctorat d'état des sciences naturelles, Univ de Provence Aix-Marseille I, 186 p.

Hasnaoui B., 1998. Régénération naturelle du chêne-liège : Difficulté et propositions de solutions. Annales de l'I.N.G.R.E.F. 1989.

Hasnaoui B., 2002. Le chêne - afarès : menaces d'extinction et possibilité de sauvegardes. Annales de l'I.N.G.R.E.F.

Hasnaoui B., Abbes C., Elghzel A et Yaccoubi W., 2004. Structure et architecture des chênaies du nord ouest tunisien : présent et devenir. Annales de l'I.N.G.R.E.F. (2004), pp : 10-15.

Hasnaoui F., 2008. Le dépérissement des chênaies du Nord - Ouest Tunisien : Diagnostic, causes et conséquences. Thèse de doctorat en science agronomique. Institut National Agronomique de Tunis. 203 p.

Hasnaoui F., Rjeibi N., Chaabane A., Yacoubi W., Hasnaoui B., 2005. Contribution à l'étude du dépérissement du chêne-liège de la forêt de Tabarka ; relations possibles entre nutriments et état sanitaire de l'arbre. I O B C /W P R S. Bulletin Vol 28 (8) 2005.

Harper J. L., 1977. Population Biology of plants. Academic Press. London, 892 p.

Institut Méditerranéenne du Liège, 2005. Une valorisation allant de la sylviculture à transformation du bois, pp : 2.

Jarret P., 2004. Guide sylviculture de la chênaie atlantique : O.N.F. 335 p.

Jian D., Wang T., Liangb Y., Renb H., Jian N et Keping M., 2003. La structure par âge de l'Epicéa des forêts le long d'un gradient d'altitude dans la partie centrale des montagnes, le nord ouest de la Chine, 13 p.

Johnson D., et Miller R., 2006. Structure and developpment of expanding westren Juniper wood londs as influenced by two topographic variables, 35 p.

Khanfouci M.S., 2005. Contribution à l'étude de la frictification et de la régénération du cédre de l'Atlas. Mémoire de mastère en sciences forestières. Université Hadj Lakhdar- Batna. Mémoire online 250 pages.

Kelen G.V., Lessard G., 2004. Entre les coupes progressives irrégulières synthèse de littérature. C.E.R.F.O Réf 3457. G l - am - 13 / 09 / 04, pp: 15-22.

Ksontini M., 1996. Etude éco physiologique des réponses à la contrainte hydrique du chêne-liège (Quercus suber L.), dans le Nord - Tunisie : comparaison avec le chêne kermès (Q. cocciféra) et le chêne zen (Q. faginea). Thèse université Paris XII, 159 p.

Lamhamdi M., Chbouki M., 1994. Les principaux) facteurs influençant la régénération naturelle du cèdre de l'Atlas (Cedrus atlantica) .Ann. Rech.For.Maroc, (1994). Pp 205-217.

Lauvaerts T., 1998. Faisabilité d'une mesure de la hauteur des arbres en forêt tropicale humide, Mémoire E.N.G.R.E.F / C.I.R.A.D - Forêt, 66 p.

Messaoudenen M., Metna B., Djouher N., 1996. La régénération naturelle de Quercus suber L. dans la forêt dominale des Béni - Ghori (Algérie). Actes du séminaire Méditerranéen sur la régénération des forêts de chêne-liège Tabarka, du 22 au 24 Octobre 1996. Annales de l'I.N.G.R.E.F 1998.

M'Hirit O., 1994. Croissance et productivité du cèdre : Approche multidimensionnelle de l'étude des liaisons stations productions. Ann. Rech. For. Maroc (1994), 11 p.

M'Hirit O., 1982. Etude écologique et forestière des cédraies du Rif Marocain. Essai sur une approche multidimensionnelle de la phytoécologie et de la productivité du cèdre (Cedrus atlantica manetti). Thèse de doctorat des sciences, uni - d'Aix - Marseille III ; 436 p+ annexes.

Nager M., 2000. Eléments de base d'une stratégie sylvopastoralisme en Afrique du Nord. C.I.H.E.A.M. Options Méditerranées.

Natividade J.V., 1950. Subériculture : Ecole national des eaux et forêts Nancy Edit Française de l `ouvrage Portugais à Subériculture, 303 p.

Nsibi R., Souayha N., Khouia L.M., Bouzid S., 2006. La régénération naturelle par semis de la subéraie de Tabarka - Ain Draham face aux facteurs écologiques et anthropiques. Geo - Eco - Trop, 2006, 30.1 : 35-48, 16 pages.

Office de Développement du Nord-Ouest (O.D.N.O), 2004. Valorisation des produits forestiers 2004 au Nord-Ouest. 67 pages.

O.N.F., 1969. Manuel d'aménagement 2éme édition. Directive concernant l'utilisation des méthodes d'aménagement dans les forêts soumises au régime forestier, 120 p.

Poskin A., 1949. Traité de sylviculture. Encyclopédie agronomique et vétérinaire. Gembloux, Duculot, 530 p.

Rose E., De Noni G., 2004. Recherches sur l'érosion hydrique en Afrique. Revue et perspectives sciences et changements planétaires /sécheresse. Volume 15, 42 p.

Rouvier H., 1977. Géologie de l'extrême Nord Tunisien tectonique et paléogéographies composées.

Saoudi A., 1994. Les bases phytosociologiques de l'aménagement et protection des chênaies tunisiennes. Thèse d'Etat. Université libre de Bruxelles (Laboratoire de botanique systématique et physiologique), 173 p.

Schütz J. P., 1997. Sylviculture 2. La gestion des forêts irrégulières et mélangées. Press polytechniques et universitaire Romandes, 177 p.

Schütz J. P., 1989. Sylviculture 1. Principe d'éducation des forêts. Press Polytechnique et universitaire Romandes, 230 p.

Schoenembeger A., 1962. Ecologie forestière Bull Station de Rech. Forest. N : 20, pp : 97. Tunis.

Selmi M., 1985. Différenciation des sols et fonctionnement des écosystèmes sur grés numidiens de Kroumirie (Tunisie). Ecologie de la subéraie zenaie. Thèse de doctorat des Sciences Naturelles Université de Nancy I, 200 p.

Sokpon N., S. Biaou. 2002. L'utilisation de diamètre des distribution en gestion, utilisation des vestiges des forêts au Bénin : cas de Bassila réserve forestière dans le Nord Bénin, 15 p.

Solignac., 1927. Etude Géologique de la Tunisie septentrionale. Direction des travaux publiés Tunis.

Teller H.L., 1967. Influence de la forêt sur les crues et Inondation, 122 p.

Turchhien B., Bruciamucchie M., 2005. La futaie irrégulière - Théorie et pratique de la sylviculture irrégulière, continue et proche de la nature Edisud : 286 p.

Vincent F., Sylvie G., Hélène D., 1998. Parcelles permanentes de recherches en forêt dense tropicale humide. Eléments pour une méthodologie d'analyse des données. C.I.R.A.D. Forêt CAMPUS International de Baillar gent B.P 5035.-.Monplier Cedex 1 France, pp : 28.

Yakoubi Hasnaoui W., 2000. Impacts des projets du fonds national de solidarité 26-26 sur les populations cibles et sur le milieu naturel environnant à Tabarka et Ain Draham. Mémoire de fin d'études du cycle de formation continue Ingénieur GREF. INAT - ISP Tabarka.58 p+ annexe.

Yessad S. A., 2000. Le chêne-liège dans les pays Méditerranée occidentale 190p.

Walter M.N., 1976. Arbres et forêts alluviales du Rhin, Bull, soc, Hist. Nar. Colomar, 55, pp: 37-88

Wetzel S., Burges D., 2001. Understory environnement and vegetation response after partial cutting and site preparation in pinus strobus L, pp: 35 - 45.

Zukrigl K., Eckhart G., Nather J et Roller M., 1963. Standrotskundliche and waldderosterreiche untersuchungen in urwaldersten der neiderosterreischissihen kalpen. Mitt. Der forest .Bundes - versnchsanstalt, 62. 244 p.

Date de l'observation :.....................Région :..................Subdivision :.....................

Triage :.......................................Série :.........................................................

Parcelle :......................................Numéro Placette :...........................................

Longitude du centre placette :............ Latitude du centre placette :...............................

Caractéristiques géographiques
Exposition :................................................................................................... Altitude :...................................................................................................... Pente (%) :....................................................................................................
Information sur le peuplement
Origine : naturelle artificielle Traitement : Futaie régulière Futaie jardinée Taillis Taillis sous futaie Stade de l'évolution : Gaulis Bas perchis Haut perchis Futaie Vielle futaie Etage principale : fermé ouvert très ouvert Essence principale :.............................................Proportion :.............................. Essences secondaires :..........................................Proportion :.............................. Densité (nbre d'arbres par 500 m²) :...................................................................... Régénération : Densité (nbre de plants par placette) :.................................................. Distance des peuplements par rapport ou douars les plus proches : très proches moyennement proches très éloignés Pâturage : nul modéré intense Incendies : oui non année :...................
Les arbres échantillons
N° arbre	C à 1,30 m	Htot	N° arbre	C à 1,30 m	Htot
1 2 3 4 5 6 7 8 8 10 11 12 13 14 15 16 17 18 19 20	................. .................. ................. ................. ................. .................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. .................. ..................	................. .................. ................. ................. ................. .................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. .................. ..................	21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40	................. .................. ................. ................. ................. .................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. .................. ..................	................. .................. ................. ................. ................. .................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. ................. .................. ..................