Évaluation de la productivité actuelle et potentielle de la pessière à épinette noire dans la forêt boréale de la ceinture d’argile à l’aide de différentes approches de cartographie à haute résolution

La forêt boréale canadienne a une importance écologique et économique considérable. Toutefois, les forêts d’épinettes noires situées dans la ceinture d’argile, une région boréale de l’est de l’Amérique du Nord, sont sujettes à la paludification. Ce phénomène est un processus naturel par lequel une couche organique s’accumule sur le sol forestier conduisant à une diminution importante de la productivité de ces forêts. Théoriquement, il existe deux types de paludification, à savoir la paludification permanente et réversible. La paludification permanente se produit dans des endroits où les conditions d’humidité du sol sont élevées (ex., reliefs plats, dépressions topographiques); alors que la paludification réversible intervient dans des sites à pente faible ou moyennement forte au fil du temps en réponse à une perturbation telle qu’un feu peu sévère. L’épaisseur de la couche organique (ECO) et la topographie constituent des paramètres clefs de la présence de la paludification dans cette région et y affectent négativement la productivité. La recherche proposée dans cette thèse consiste à approfondir la compréhension et la détection du phénomène de paludification dans les forêts d’épinette noire dans une perspective de maintien ou d’accroissement de la productivité des arbres. Cette thèse a pour objectif principal de déterminer et de sélectionner les variables permanentes des site permettant d’expliquer les écarts de productivité de la pessière à épinette noire observés dans la ceinture d’argile à l’aide des méthodes à haute résolution et des données recueillies sur le terrain. L’expression de ces variables a été ensuite utilisée pour prédire la productivité actuelle et potentielle des sites soumis à la paludification. Les objectifs spécifiques de cette thèse étaient de (1) détecter et identifier d’une manière continue l’interface sol minéral/couche organique afin de cartographier la topographie du sol minéral à l’échelle des sites paludifiés; (2) étudier d’une façon quantitative les relations entre l’ECO et la topographie (au niveau du sol minéral et de la surface) afin de caractériser ces relations à l’échelle du paysage, notamment la distribution et la variabilité spatiale de l’ECO; (3) identifier les variables topographiques permettant de distinguer et de cartographier la paludification réversible et la paludification permanente à l’échelle du paysage; et (4) évaluer l’effet de l’ECO et des variables topographiques, exprimées à différentes résolutions spatiales, sur la productivité forestière des forêts paludifiées de la ceinture d’argile. Le but ultime est d’améliorer notre compréhension de la façon dont ces variables ainsi que leurs résolutions influencent la productivité des arbres dans les forêts d’épinette noire. Les résultats du premier chapitre de la thèse ont démontré que la méthode géophysique géoradar, ayant une bonne corrélation de ses résultats avec les données du terrain (r = 0,93; P < 0,001), a permis d’obtenir une cartographie précise, continue et fiable de l’interface couche organique/sol minéral dans des sites faiblement à modérément paludifiés. Cependant, en dépit de son incapacité à cartographier l’interface couche organique/sol minéral dans les sites hautement paludifiés, le recours au géoradar s’est révélé pertinent dans la mise en évidence de l’interface horizon fibrique/couche organique et de sa continuité spatiale. Cela rend le géoradar particulièrement intéressant dans la détection des niveaux d’entourbement constituant ainsi une méthode de détection indirecte prometteuse pour l’aménagement des forêts paludifiées. Le deuxième objectif a été abordé dans deux différents chapitres (II et III) à l’aide d’une approche quantitative de modélisation de l’ECO par arbre de régression. Différentes variables topographiques (élévation, pente, exposition, indice topographique d’humidité (TWI), courbure totale, courbure transversale, et courbure horizontale) ont été utilisées dans les modèles sélectionnés des deux chapitres. D’une façon générale, nous avons démontré que les topographies de surface et du sol minéral influencent l'accumulation de la couche organique à l'échelle du paysage dans la ceinture d’argile. Les résultats du deuxième chapitre ont permis de délimiter les principaux patrons de l’ECO et d’élucider trois relations spatiales entre l’ECO et les variables explicatives: (i) les zones avec une couche organique épaisse (62 cm) avaient des pentes douces (≤ 1,8%); (ii) les zones avec pentes plus raides (> 3,2 %) ont été associées à une couche organique peu profonde (27 cm); et (iii) les résultats les plus significatifs ont été obtenus avec des résolutions 10 et 20 m en comparaison au 1 et 5 m. Le troisième chapitre a permis de mettre en évidence les différentes relations entre la topographie du sol minéral et l’ECO à l’échelle du paysage. La construction d’un modèle numérique d'élévation au niveau du sol minéral à l’échelle du paysage est un élément central de notre démarche. Les modelés développés nous permettent d’affirmer que : (i) la pente du sol minéral, la composition du sol minéral (argile, till et régolithe), le TWI et l’exposition sont les quatre principales variables influençant l’accumulation de la couche organique; (ii) les valeurs seuils de pente du sol minéral > 3,5% et ≤ 2% permettent respectivement de distinguer les zones les plus prometteuses et les plus vulnérables pour l’aménagement forestier; (iii) les zones avec une exposition nord et est étaient associées à une couche organique plus profonde par rapport à celles exposées vers le sud et l’ouest; et (iv) la distinction entre les zones paludifiées et non paludifiées sur la base d’une valeur seuil de la pente du sol minéral de l’ordre de 3,5% constitue un des apports majeurs de cette étude. Afin de répondre au troisième objectif, une approche semi-automatique de subdivision sous SIG du territoire à l’étude en des entités du paysage distinctes a été réalisée en combinant des données topographiques, notamment l’indice topographique de position (TPI), l’indice topographique d’humidité (TWI) ainsi que la pente de surface. Cette approche s’est révélée efficace, car elle a permis de délimiter des entités possédant des caractéristiques géomorphologiques semblables, notamment en terme de susceptibilité à l’accumulation de la couche organique, et par conséquent ont été assignées à l’un ou l’autre type de paludification, soit réversible ou permanente. Un apport majeur de cette approche semi-automatisée est la mise en évidence de deux sous-entités statistiquement différentes (le test HSD de Tukey, P < 0,001), à savoir des dépressions ouvertes préférentiellement drainées (paludification réversible) et des dépressions fermées potentiellement engorgées (paludification permanente) du fait de leurs positions topographiques et conditions d’humidité. Cela rend l'outil développé particulièrement utile pour la mise en oeuvre des stratégies d’aménagement durable dans les forêts paludifiées. Pour atteindre le quatrième objectif, deux modèles ont été explorés pour la modélisation de la productivité (exprimée par l’indice de qualité de station (IQS) dans notre cas) en utilisant une approche par arbre de régression. Le premier modèle contient les variables topographiques et l’ECO, alors que le deuxième modèle inclut seulement les variables topographiques issues des données LiDAR. Les résultats de cette modélisation ont démontré que l’ECO, l'exposition et la pente sont les trois variables les plus importantes pour expliquer la productivité forestière à l’échelle du paysage; et pour déterminer des seuils d’ECO et des variables topographiques qui permettent de caractériser, à la fois, des zones productives et improductives. En effet, les zones avec une productivité élevée étaient associées à une couche organique peu profonde ( 85 cm), favorisant l’invasion de mousses et de sphaignes. Du point de vue de l’échelle (résolutions), le premier modèle semble relativement indépendant de l’échelle, alors que la réponse du deuxième modèle augmentait significativement avec la taille du pixel. Ces résultats pourraient donc être appliqués à des échelles opérationnelles et là où des informations sur l’ECO sont disponibles afin de prédire la productivité. Des cartes thématiques prédictives de la distribution spatiale de la productivité ont été réalisées avec nos deux modèles. Les résultats de cette thèse ont permis d'approfondir les connaissances sur les variables permettant d’expliquer les écarts de productivité dans la pessière à épinette noire ainsi que sur l’importance des variables topographiques dans la modélisation de l’ECO et la productivité à l’échelle du paysage. De plus, cette étude a permis de caractériser la distribution spatiale des deux types de paludification (permanente et réversible) à l’échelle du paysage. Quoique cette étude s’intéresse plus particulièrement à la pessière à épinette noire, les connaissances acquises pourraient être applicables à d’autres territoires paludifiés de la forêt boréale

Recherche par télédétection géologique de sites potentiels d'hydrocarbures dans l'anti-atlas marocain

Les bassins de Ouarzazate et de Tindouf, respectivement au nord et au sud de l'Anti-Atlas, présentent des indices d'hydrocarbures dans les formations paléozoïques. Il s'agit de réservoirs carbonates géorgiens (Cambrien inférieur) et gréseux géorgiens, acadiens (Cambrien moyen) et ordovicien. Les mêmes formations paléozoïques constituent une couverture partielle du socle dans l'Anti-Atlas et dessinent un grand pli complexe à l'ouest de Zagora et au sud de Ouarzazate. L'image Landsat-TM du 08 avril 1994 de cette région de l'Anti-Atlas central a été utilisée pour préciser et raffiner, par télédétection géologique, la carte géologique existante tant au point de vue lithologique que structural. Après une correction radiométrique pour la lumière diffuse, la division des bandes TM1, TM3, TM4, TM5 et TM7 par la bande TM2 minimise les effets topographiques sur la réflectance et permet d'obtenir une discrimination lithologique maximale; la bande TM2 était celle qui présente le moins de variation de réflectance d'une lithologie à l'autre. Une image en fausses couleurs constituée à partir des deux premières composantes principales des 5 néo-canaux obtenus et du néo-canal 3/2 a permis la distinction de 27 unités stratiformes dans la couvertures infra-cambrienne et paléozoïque. Certaines recettes de bandes, déjà utilisées sur des terrains arides à lithologies semblables, sont applicables dans la zone d'étude. Le plissement principal forme un grand pli en S à axe généralement est-ouest dans la demie nord de la région. D'autres plis secondaires transverses sont visibles dans la demie sud et des ensellements des charnières des plis principaux dessinent des dômes et des bassins. Les principales failles sont identifiées par inspection de la CP1 des canaux originaux. C'est dans la partie sud de la région étudiée, montrant une stratigraphie complète, que se concentrent les diverses structures d'intérêt pétrolier identifiées, essentiellement des structures anticlinales en association avec les roches réservoirs potentielles encore coiffées de leur couverture imperméable respective de la zone d'étud