Abstracts
Résumé
Dans la région du lac Guillaume-Delisle, au Québec subarctique, des cônes de déjection d’âge et d’altitude variables se sont déposés tout au cours de l’Holocène sur des plages soulevées, offrant ainsi la possibilité d’étudier une chronoséquence holocène de sols. Le cadre temporel a été établi par la datation radiocarbone de 38 horizons organiques enfouis et la mise à jour de la courbe d’émersion des terres pour la région. L’émersion se serait effectuée selon un taux moyen de 5,6 m/siècle entre 8500 et 4500 années cal. BP et de 1,2 m/siècle après 4500 années cal. BP. Les analyses morphologiques et chimiques ont montré que les sols et paléosols subarctiques sableux, bien drainés et soustraits à l’influence des combes à neige, sont formés principalement par l’accumulation dans l’horizon B de matière organique humifiée combinée à de l’aluminium et du fer (podzolisation). L’enrichissement de l’horizon B en fer et en aluminium se fait lentement en suivant une fonction linéaire. La durée de pédogenèse requise pour obtenir le seuil minimal (Fep+Alp) de 0,4 % pour l’horizon Bf fixé par le Système canadien de classification des sols serait d’environ 12 400 années. Bien que théorique, cette durée indique à quel point les conditions nécessaires à la formation des podzols sont difficiles à réunir dans les sols subarctiques où la végétation demeure basse et le couvert nival peu épais.
Abstract
In the Lake Guillaume-Delisle area, alluvial fans of various ages and elevations were deposited throughout the Holocene over raised beach deposits, which allowed the studying of a subarctic soil chronosequence. The temporal framework is based on 38 14C dates obtained from buried organic soil horizons, and on an update of the Lake Guillaume-Delisle emersion curve. Postglacial uplift due to glacio-isostatic rebound occurred at a mean rate of 5.6 m/century from 8500 to 4500 cal. yr BP, and of 1.2 m/century after 4500 cal. yr BP. Soil morphology and chemistry showed that well-drained sandy soil and paleosoil development, when not influenced by late-lying snowpatches and snowmelt, is characterized by B horizons having a dominant accumulation product of amorphous material composed mainly of humified organic matter combined with Al and Fe (podzolisation). However, pedogenic development is weak and would take about 12 400 years to reach the minimum values of (Fep+Alp) ≥0,4% required for meeting the criteria of a Bf horizon in the Canadian system of soil classification, according to a linear function. Although theoretical, this value indicates that conditions required for Podzol formation are hard to reach in subarctic soils with low vegetation and thin snow cover.
Article body
Introduction
La podzolisation fait partie des mécanismes de transformation des sols bien drainés des environnements nordiques. Son intensité diminue avec la latitude (Ugolini, 1986 ; Ugolini et al., 1987). Dans l’Est du Canada, Moore (1978) a montré que des Podzols répondant aux critères du Système canadien de classification des sols (Groupe de travail sur la classification des sols, 2002) se développaient préférentiellement dans la forêt boréale et rarement dans la toundra forestière située plus au nord. Des Podzols répondant aux critères du Système canadien de classification des sols ont tout de même été observés dans le domaine de la toundra forestière dans des sols basiques à acides (voir Payette, 1973 ; Protz et al., 1984).
Bien que la podzolisation ait été identifiée comme un processus actif en milieu subarctique, peu d’études ont permis d’en mesurer l’évolution temporelle. Les chronoséquences de sols subarctiques publiées à ce jour portent principalement sur les stades initiaux de la pédogenèse (<250 ans) et sur l’évolution de la concentration en nutriments (Stork, 1963 ; Viereck, 1966, 1970). À plus long terme, le taux de podzolisation — tel qu’indiqué par l’accumulation de fer, d’aluminium et de carbone organique dans l’horizon B (Groupe de travail sur la classification des sols, 2002) — demeure peu connu pour les sols de la toundra forestière et de la toundra arbustive. Selon Protz (1984), la formation de Podzols répondant aux critères du Système canadien de classification des sols de la toundra forestière nécessiterait une durée de pédogenèse minimale de 2300 ans.
Les taux de pédogenèse sont difficiles à évaluer, car ils reposent sur l’établissement de chronoséquences de sols auxquelles plusieurs pré-requis sont associés. D’abord, l’âge des surfaces initiales doit être établi avec précision afin de permettre de mesurer quantitativement la pédogenèse à l’aide de chronofonctions (Catt, 1991). Ensuite, l’établissement d’une chronoséquence de sols repose sur le postulat selon lequel, à l’exception du temps, tous les facteurs de pédogenèse (climat, matériel parental, topographie, organismes) sont identiques entre les profils de sols étudiés (Jenny, 1941).
Dans la région du lac Guillaume-Delisle, au Québec subarctique (fig. 1A-B), des cônes de déjection d’origine torrentielle se sont déposés au cours de l’Holocène sur des plages soulevées ayant émergé en réponse au relèvement isostatique postglaciaire (Lafortune et al., 2006). Ces cônes de déjection d’âge variable permettent l’étude d’une chronoséquence de sols à partir de la datation au radiocarbone de nombreux horizons organiques enfouis. L’homogénéité des conditions édaphiques des dépôts de cônes de déjection, et notamment en ce qui concerne le matériel parental (granulométrie et minéralogie), suggère la possibilité de mesurer quantitativement l’évolution de la pédogenèse à l’aide de chronofonctions.
Localisation de la région du lac Guillaume-Delisle et des sites d’étude. Les noms de lieux mentionnés dans le texte sont également indiqués.
Location of the Lake Guillaume-Delisle area and of the two study sites. Place names mentioned in the text are also indicated.
Les sols enfouis sous les cônes de déjection, en fournissant un âge minimal de mise en place des surfaces, offrent aussi la possibilité de préciser la courbe d’émersion des terres de la région du lac Guillaume-Delisle (Hillaire-Marcel, 1976) en renseignant sur l’âge approximatif des surfaces émergées. Les courbes d’émersion disponibles pour le Québec nordique reposent surtout sur la datation de coquillages marins et peu sur des datations de sols enfouis (Gray et al., 1993). Les environnements littoraux sableux récents du lac Guillaume-Delisle sont cependant pauvres en coquillages marins et la courbe d’émersion des terres est basée principalement sur des datations de coquillages, toutes plus anciennes que 4500 années cal. BP (Hillaire-Marcel, 1976). La datation d’horizons organiques enfouis offre une avenue intéressante pour la période postérieure à 4500 années cal. BP dans ces environnements sableux.
Cette étude a pour objectifs (1) de caractériser le développement des sols subarctiques sableux bien drainés et soustraits à l’influence des combes à neige à partir de l’analyse d’une chronoséquence holocène et de sa relation avec les autres facteurs de pédogenèse et (2) d’établir le cadre temporel de mise en place des plages soulevées en mettant à jour la courbe d’émersion des terres de la région du lac Guillaume-Delisle.
Région et sites d’étude
Les sols enfouis au sein de six cônes de déjection localisés dans les vallées de l’Inverbule et de l’Amphithéâtre ont été étudiés (fig. 1C et 2). Les vallées sont incisées dans des cuestas d’âge protérozoïque qui reposent en discordance sur un socle de granite d’âge archéen et qui sont formées, de la base au sommet, de séquences d’arkose rose, de dolomie et de basalte (Woodcock, 1960 ; Chandler et Schwarz, 1980). Les revers des cuestas s’inclinent vers la baie d’Hudson et leur front, tourné vers l’est, atteint une altitude maximale de 420 m.
Après une brève période de stagnation du front glaciaire vers 8300 années cal. BP, les vallées ont rapidement été déglacées (Lajeunesse et Allard, 2003a). Les eaux de la mer de Tyrrell ont alors inondé les vallées de l’Amphithéâtre et de l’Inverbule jusqu’à une altitude respective de 254 m et 270 m au-dessus du niveau marin actuel (Hillaire-Marcel, 1976). La mer s’est ensuite retirée en raison d’un relèvement isostatique rapide dont le taux fut de l’ordre de 6,5 m/siècle jusque vers 6000 années 14C BP, puis de 1 m/siècle jusqu’à nos jours (Hillaire-Marcel, 1976). De nouvelles données présentées dans ce travail viendront préciser ces taux d’émersion pour la période récente. Le retrait des eaux de la mer de Tyrrell a laissé une succession de plages soulevées dans plusieurs vallées de la région, notamment dans la vallée de l’Amphithéâtre. Les sédiments quaternaires qui recouvrent les vallées sont essentiellement d’origine marine et littorale avec, par endroits, un remaniement éolien ou nivéo-éolien (Filion et Morisset, 1983).
Au cours de l’Holocène, des cônes de déjection ont recouvert, à l’occasion d’événements torrentiels, les sédiments marins et littoraux des vallées de l’Amphithéâtre et de l’Inverbule. La vallée de l’Amphithéâtre compte cinq cônes de déjection, alors que la vallée de l’Inverbule n’en compte qu’un seul (fig. 2 ; tabl. I). Quatre des cinq cônes de la vallée de l’Amphithéâtre (CF1 à CF4) et celui de la vallée de l’Inverbule ont une forme conique simple résultant de la superposition de dépôts torrentiels. Le cône de déjection de la vallée de l’Amphithéâtre est imbriqué, c’est-à-dire qu’il est formé d’une succession de cônes simples emboîtés qui ont été progressivement abandonnés suite à l’encaissement permanent du chenal principal (fig. 2A et 3 ; Lafortune et al., 2006). Ces cônes de déjection résultent d’écoulements torrentiels en nappe, d’écoulements hyperconcentrés ou d’écoulements intermédiaires entre les deux premiers types (Lafortune et al., 2006). En règle générale, les cônes formés par ces types de dépôts se forment rapidement lors d’orages violents (Harvey, 1986 ; Blair, 1987 ; Wells et Harvey, 1987 ; Ballantyne et Whittington, 1999 ; Blair, 1999).
Blocs diagrammes montrant les diverses étapes de la formation des cônes de déjection dans la vallée de l’Amphithéâtre en relation avec le relèvement isostatique postglaciaire.
3D sketch showing various stages of alluvial fan development in the Amphithéâtre valley in relationship with postglacial land emersion.
Les données provenant de la station météorologique de Kuujjuarapik, située à une centaine de kilomètres au sud du lac Guillaume-Delisle (fig. 1B), indiquent que la température moyenne annuelle et les précipitations totales annuelles sont respectivement de ‑4,4 °C et de 649 mm (Environnement Canada, 2004). Les précipitations tombent sous forme de neige dans une proportion de 36 %. Les vents de tempête et les vents anticycloniques du sud-ouest, de l’ouest et du nord-ouest soufflent la neige et la redistribuent dans les endroits protégés et les dépressions topographiques (Filion et Payette, 1976).
Dans les vallées à l’étude, des combes à neige d’une superficie de 209 à 3 045 m2 sont présentes dans les bassins versants du cône de déjection de la vallée de l’Inverbule, du cône CF1 et du cône imbriqué de la vallée de l’Amphithéâtre. Une combe à neige unique, adossée à une paroi rocheuse, constitue l’essentiel des bassins versants des cônes CF2, CF3 et CF4. Les congères accumulées dans les combes à neige constituent un apport hydrique important favorisant la pédogenèse au moment de la fonte qui peut se produire jusque tard dans l’été (Payette et Filion, 1993b).
La zone d’étude se situe dans la zone de pergélisol discontinu et dispersé, c’est-à-dire là où le pergélisol couvre moins de 50 % des surfaces (Allard et Séguin, 1987). Sur les hautes terres, le pergélisol peut atteindre une profondeur de 100 m, alors qu’il est absent des sites bien drainés des basses terres (Allard et Séguin, 1987).
Les vallées étudiées font partie du domaine de la toundra forestière maritime (Payette, 1976), caractérisé par la présence de l’épinette blanche [Picea glauca (Moench.) Voss]. Adaptée au climat maritime (Payette, 1975 ; Payette et Filion, 1975), cette espèce prédomine dans les sites de basse altitude. Plus haut en altitude ou plus au nord, l’épinette noire [Picea mariana (Mill.) B.S.P.] sous forme prostrée remplace graduellement l’épinette blanche et la toundra forestière cède le pas à la toundra arbustive (Payette et Filion, 1975). Il n’y a actuellement aucun arbre à la surface actuelle des cônes de déjection étudiés. La majeure partie des cônes de déjection est couverte de lichénaies. Le bouleau glanduleux (Betula glandulosa Michx) est parfois présent dans le chenal principal, dans les chenaux secondaires et au pourtour des combes à neige.
Méthodologie
Échantillonnage sur le terrain
Dans la vallée de l’Amphithéâtre, 14 profils de sol ont été dégagés dans les dépôts du cône de déjection imbriqué (N1 à N6 et S1 à S8) et un dans chacun des cônes situés au pied de la cuesta (CF1 à CF4 ; fig. 2A-B). Six profils ont été dégagés dans le cône de déjection de la vallée de l’Inverbule (Inv1 à Inv6 ; fig. 2C). L’altitude de chaque profil a été déterminée à l’aide d’un altimètre numérique (précision de 0,1 m) lors de courtes séquences de mesures et en partant de la ligne des hautes eaux.
Les propriétés morphologiques (profondeur/limite des horizons et structure des sédiments) des sols de surface et des paléosols ont été décrites pour chaque profil, selon la nomenclature du Système canadien de classification des sols (Groupe de travail sur la classification des sols, 2002). Les horizons organiques ont été échantillonnés pour la datation au radiocarbone et les horizons minéraux à des fins d’analyses en laboratoire, à l’exception des horizons trop minces pour être isolés, notamment les horizons Ae.
Analyses en laboratoire
Les échantillons minéraux ont été tamisés et seule la fraction inférieure à 2 mm a été retenue pour effectuer les analyses. La granulométrie de la matrice a été déterminée suivant la méthode décrite par Head (1992) et la granulométrie moyenne suivant la méthode de Folk et Ward (1957), exprimée selon la classification de Folk (1974). La minéralogie a été évaluée sommairement à la loupe binoculaire. La couleur de l’échantillon a été déterminée en laboratoire sur la fraction tamisée, à l’état sec et à l’état humide, à l’aide d’une charte Munsell.
Le pH CaCl2 (1 : 2) a été déterminé pour la majorité des échantillons. La teneur en carbone organique (Corg) a été mesurée pour tous les échantillons par la méthode de Walkley-Black (Walkley et Black, 1934). Les teneurs en fer et en aluminium ont été obtenues par extraction au pyrophosphate de sodium (Fep, Alp), à l’oxalate d’ammonium acide (Feo, Alo) et au dithionite-citrate-bicarbonate (DCB) (Fed). Toutes les analyses ont été effectuées selon le protocole de McKeague (1978). Le pyrophosphate de sodium extrait principalement les complexes organo-métalliques (McKeague et al., 1971). L’oxalate d’ammonium acide extrait une portion des complexes organo-métalliques de fer et d’aluminium, mais aussi les allophanes, les ferrihydrites, l’imogolite et la proto-imogolite (Wang, 1990). Comme l’oxalate d’ammonium acide dissout la magnétite du matériel parental (Baril et Bitton, 1969 ; Walker, 1983), cette dernière a préalablement été extraite à l’aide d’un aimant rotatif. Le DCB extrait les oxydes de fer cristallins en plus des fractions déjà extraites par l’oxalate (McKeague et Day, 1966).
Datations au radiocarbone
Les horizons organiques enfouis des profils de sols ont été datés au radiocarbone (tabl. II) et les datations 14C obtenues ont été étalonnées à l’aide du Radiocarbon calibration program rev5.0.1 (Stuiver et al., 2005). L’étalonnage a été effectué à partir des datations conventionnelles, c’est-à-dire corrigées pour le fractionnement isotopique. Le temps de séjour de la matière organique dans un sol (c’est-à-dire l’âge moyen de la matière organique en décomposition constituant un horizon H) doit être considéré, les humus livrant un âge apparent qui diffère de l’âge réel du sol (Paul et al., 1964 ; Matthews, 1980 ; Kinnard et Lewkowicz, 2006). En raison de l’incertitude liée au temps de séjour et aux différents prétraitements des échantillons (Martin et Johnson, 1995 ; Wang et Hsieh, 2002), nous présumons que les datations 14C obtenues sur les paléosols fournissent un âge minimal pour l’émersion des terres et un âge maximal pour l’enfouissement d’une surface donnée.
Pour l’établissement de la courbe d’émersion, des échantillons provenant du domaine marin (coquillages) et datés au 14C ont aussi été utilisés. En plus de l’étalonnage des datations au radiocarbone à partir des datations conventionnelles, la différence (∆R) entre l’âge apparent local et celui de l’océan a été prise en compte lors de l’utilisation du logiciel de Stuiver et al. (2005). Le ∆R utilisé provient de coquillages échantillonnés dans la région et il correspond à la datation moyenne de 135 ± 50 années obtenue pour trois sites sur la côte est de la Baie d’Hudson (181 ± 40, 131 ± 50 et 92 ± 50 années ; Barber et al., 1999 ; Reimer, 2005).
Établissement de la chronoséquence de sols et des chronofonctions
La chronoséquence de sols et les chronofonctions correspondantes sont de type post-incisives (Vreeken, 1975), c’est-à-dire qu’elles portent sur des sols qui ont commencé à se former à différentes périodes dans le passé et qui sont toujours exposés. La chronoséquence repose donc sur l’évolution de propriétés de l’horizon B (ou C pour les régosols) des sols de surface en fonction du temps.
La chronoséquence analysée dans cette étude porte sur neuf des 24 profils de sols échantillonnés. Ces neuf sols sont tous bien datés et situés sur des surfaces bien drainées qui ne sont pas localisées dans une combe à neige. Ce n’est pas le cas des 15 autres profils qui ont été exclus de la chronoséquence en raison de la présence d’une combe à neige (S8, Inv6 et CF4) ou de problèmes concernant l’établissement de la durée de la pédogenèse liés à l’absence de datations au radiocarbone (S1, S2, S3, N4, N5, N6, Inv1, Inv2 et Inv3) ou à des inversions de dates (CF2, N2 et S7). Les sols qui n’ont pas été utilisés pour la chronoséquence ont tout de même été analysés afin d’évaluer la pédogenèse au sein de chaque profil (par exemple, tendances à l’enrichissement ou à la déplétion des propriétés de sol dans l’horizon B).
Des régressions linéaires multiples pas à pas ont été utilisées afin d’obtenir des équations prédictives (chronofonctions) pour les indicateurs pédogénétiques dans l’horizon B (ou C pour les régosols) à l’aide de la procédure REG du logiciel SAS 9.1 (SAS Institute inc.).
Les variations dans les conditions édaphiques, soit la topographie, le matériel parental, le climat et les organismes (Jenny, 1941), ont été prises en compte lors des régressions multiples pas à pas.
Pour la topographie, l’altitude médiane prévalant lors de la formation du sol [(altitude lors de la mise en place des surfaces + altitude lors de l’enfouissement ou altitude actuelle) / 2] a été retenue comme variable indépendante. La topographie locale n’a pas été utilisée puisqu’elle était similaire entre les profils, les profils S8, Inv6 et CF4 situés dans des combes à neige ayant déjà été exclus de la chronoséquence.
Pour le matériel parental, la granulométrie est similaire entre les profils avec une proportion de particules <63 µm inférieure à 5 % (tabl. III), ce qui, dans la majorité des cas, se situe sous le seuil de détectabilité de la méthode utilisée pour déterminer la granulométrie des sédiments fins. Les variations dans le contenu en particules fines au sein des profils n’ont donc pu être utilisées. La granulométrie moyenne de la matrice (<2 mm) a cependant été retenue puisque ce facteur varie légèrement entre les profils (tabl. III). La composition minéralogique était la suivante : quartz, feldspaths potassiques et plagioclases, magnétite, biotite — parfois chloritisée dans les horizons Bm, muscovite et amphiboles et/ou pyroxènes. Elle était similaire d’un profil à l’autre. Elle n’a pas été prise en compte pour l’établissement de chronofonctions.
Pour ce qui est du climat et de la végétation (organismes), ces facteurs ont connu des variations importantes au cours de l’Holocène (voir Payette et al., 1985 ; Allard et Séguin, 1987 ; Bégin et Payette, 1988 ; Filion et al., 1991 ; Payette et Filion, 1993a), mais ces changements peuvent difficilement être quantifiés et n’ont pu être inclus dans l’analyse statistique.
Établissement de la courbe d’émersion
Onze datations 14C ont été obtenues sur l’horizon Hb de plages enfouies situées entre 10 et 196 m d’altitude. Ces datations 14C du domaine terrestre ont été utilisées pour mettre à jour la courbe d’émersion des terres, de même que neuf datations 14C du domaine marin, soit huit datations obtenues de coquillages et déjà publiées (tabl. II ; Blake, 1988 ; Lajeunesse et Allard, 2003b) et une datation inédite obtenue sur des coquillages échantillonnés dans une plage soulevée, à 218 m d’altitude, dans la partie amont de la vallée de l’Inverbule. Il s’agit du site le plus élevé à avoir été daté à ce jour dans cette région.
Le début de l’émersion des terres a été estimé à partir de la datation de transgression marine la plus ancienne pour la région (I‑9005 : 8370 (+290/‑310) années cal. BP) obtenue sur des coquilles de Macoma calcarea récoltées dans les premiers 50 cm d’argile franche dans la région de Kuujjuarapik (tabl. II ; fig. 1A ; Hillaire-Marcel, 1976). Cette datation a été associée à l’altitude des plus hautes plages dans la région à l’étude (270 m) et reportée sur la courbe d’émersion. La remontée eustatique postglaciaire n’a pas été prise en compte.
La courbe d’émersion obtenue fut comparée à celle de Hillaire-Marcel (1976) qui est basée sur le postulat d’une périodicité moyenne de 45 ans pour la formation des cordons de plages établie à partir de quatre datations 14C provenant directement de la région du lac Guillaume-Delisle. Pour cette courbe, aucune transformation n’a été effectuée puisque l’échelle de temps était considérée par l’auteur comme étant en années sidérales. Notre courbe d’émersion fut aussi comparée à la courbe du détroit de Manitounuk (fig. 1B) de Allard et Tremblay (1983a) qui a été construite à partir de nombreuses datations dont les marges d’erreur altitudinales et temporelles ont été prises en compte. Les datations conventionnelles de Allard et Tremblay (1983a) ont été étalonnées de sorte que le tracé de la courbe a été légèrement modifié par rapport à celui de la publication originale.
Résultats
Séquence sédimentaire et caractéristiques morphologiques
Dans les trois sites d’étude, tous les profils ont été excavés dans des dépôts grossiers allant des sables moyens aux sables grossiers avec graviers. À l’exception des sols S3 et N5 qui ne présentent pas de paléosols, tous les profils étudiés montrent, du bas vers le haut, la séquence sédimentaire suivante : un dépôt de plage, un mince horizon organique enfoui (avec ou sans charbons de bois) ou une discordance d’érosion, un dépôt de cône de déjection de type torrentiel ou de combe à neige dans lequel le sol actuel s’est développé (tabl. III). Plusieurs profils montrent une succession de dépôts de cône séparés par un horizon organique, correspondant à autant d’événements sédimentaires suivis d’une stabilisation par le couvert végétal. Les sols S3 et N5 se sont formés uniquement dans des dépôts de cônes de déjection.
Les caractères morphologiques des sols sont faiblement développés, mais présentent tout de même une faible coloration de l’horizon B et la présence d’un horizon éluvié Aej dans plusieurs profils (tabl. III).
Caractéristiques chimiques et classification des sols et des paléosols
Les sols sont acides avec un pH CaCl2 passant d’un minimum de 3,5 dans l’horizon B à un maximum de 5,4 dans l’horizon C. Les concentrations en Corg, Fep, Alp, Feo et Alo sont généralement plus élevées dans l’horizon B, mais elles demeurent faibles (tabl. IV et V). Les concentrations en (Feo-Fep) et (Alo-Ap) (correspondant aux allophanes et ferrihydrites pour Fe et à l’imogolite et proto-imogolite pour Al) ne montrent pas de tendance générale à l’enrichissement dans l’horizon Bm (tabl. IV). En conséquence, les concentrations en Feo et Alo n’apportent pas d’informations supplémentaires dans les sols étudiés puisqu’elles dépendent principalement des variations dans la concentration en complexes organo-métalliques (Fep et Alp) et non en allophanes, ferrihydrites, imogolite et proto-imogolite ((Feo-Fep) et (Alo-Ap)). De même, la concentration en Fed ne serait pas un bon indicateur pédogénétique puisque l’analyse des profils (tabl. IV) ne permet pas de confirmer la tendance à l’enrichissement de l’horizon B indiqué par les moyennes pour Fed (tabl. V). Les concentrations en Fed, (Feo-Fep) et (Alo-Ap) seraient liées aux variations dans le matériel parental. Ainsi, seuls le pH CaCl2 et les concentrations en Corg, Fep et Alp caractérisent la pédogenèse des sols et des paléosols étudiés. Ce sont donc ces variables qui seront utilisées pour l’établissement des chronofonctions.
Les sols et paléosols étudiés font partie des brunisols dystriques orthiques ou éluviés (Groupe de travail sur la classification des sols, 2002). Les profils de base correspondant aux plages enfouies dans les profils S1, S2 et N3, de même que le sol de surface du cône de déjection du profil S5, ne montrent aucun indice de développement pédogénétique et sont des régosols orthiques.
Certains sols et paléosols pourraient faire partie des cryosols statiques dystriques brunisoliques ou des cryosols statiques régosoliques si le pergélisol était présent à moins de 1 m de profondeur par le passé (Groupe de travail sur la classification des sols, 2002), comme le suggère la présence de fentes de gel fossiles dans les profils N1, S1 et S5 (Groupe de travail sur la classification des sols, 2002).
Des sols podzoliques à plaques (Payette et Filion, 1993b) caractérisent certains profils dont la matrice contient peu ou pas de graviers (profils S7, N6, CF2, CF1 et Inv1). Les dépôts grossiers ne semblent pas permettre le développement de sols à plaques.
Temps et pédogenèse
La chronoséquence de sols montre que les premiers signes morphologiques de podzolisation apparaissent à partir de 570 ans avec la formation d’un horizon Aej (profil S5 ; tabl. III). Les indices chimiques de podzolisation suivent de près, montrant un enrichissement à partir de 650 ans (profil S6).
La procédure pas à pas des régressions multiples effectuée sur les indicateurs pédogénétiques (pH CaCl2, Corg, Fep Alp et (Fep+Alp) ; tabl. VI) a montré que seule la concentration en Alp montre une relation significative avec le temps (p = 0,0144), les deux autres variables (granulométrie moyenne de la matrice et altitude médiane) n’étant pas significatives. La chronofonction obtenue est la suivante :
Pour le pH CaCl2 et les concentrations en Corg et Fep, aucune variable n’a été retenue par la procédure pas à pas. Des régressions linéaires simples montrent cependant que les concentrations en Corg et Fep augmentent avec la durée de la pédogenèse (fig. 4), bien que les relations ne soient pas significatives (p < 0,05). C’est aussi le cas de la relation entre (Fep+Alp) — une variable utilisée par le Système canadien de classification des sols (Groupe de travail sur la classification des sols, 2002) pour identifier les Podzols — et la durée de la pédogenèse. La procédure pas à pas n’a retenu aucune variable permettant d’expliquer les variations dans les concentrations en (Fep+Alp), mais elle montre une relation marginalement significative avec le temps (p = 0,06 ; fig. 4). Les équations obtenues lors des régressions simples sont les suivantes :
D’après l’équation 5, la durée de pédogenèse requise pour obtenir le seuil minimal (Fep+Alp) de 0,4 % pour l’horizon Bf fixé par le Système canadien de classification des sols (Groupe de travail sur la classification des sols, 2002) serait de 12 398 années, toutes autres choses étant égales par ailleurs. Cette valeur demeure théorique, notamment en raison du fait que la relation entre (Fep+Alp) et la durée de pédogenèse n’est que marginalement significative.
Durée de la pédogenèse et propriétés chimiques des sols sableux bien drainés et exempts de l’influence de combes à neige dans les vallées de l’Amphithéâtre et de l’Inverbule. (A) Corg, (B) Fep, (C) Alp, (D) Fep+Alp et (E) pH CaCl2. L’échelle de l’ordonnée diffère en (A) et (E).
Time and soil chemical properties of well-drained soils not located in late-lying snowpatches environments in the Amphithéâtre and Inverbule valleys. (A) Corg, (B) Fep, (C) Alp, (D) Fep+Alp and (E) pH CaCl2. The ordinate scale differs in (A) and (E).
Courbe d’émersion des terres
La courbe d’émersion des terres élaborée dans ce travail pour la région du lac Guillaume-Delisle précise la portion postérieure à 4500 années cal. BP. L’émersion s’est effectuée selon un taux moyen de 5,6 m/siècle entre 8500 et 4500 années cal. BP et de 1,2 m/siècle après 4500 années cal. BP (fig. 5).
Les courbes d’émersion de la région du lac Guillaume-Delise, c’est-à-dire celle de Hillaire-Marcel (1976) et la courbe d’émersion mise à jour ici, sont similaires et se trouvent dans l’intervalle établi pour le détroit de Manitounuk (Allard et Tremblay, 1983b) pour la portion postérieure à 4500 années cal. BP (fig. 5). Une de nos datations du domaine terrestre obtenue pour une plage enfouie à 23 m d’altitude (S6) a toutefois donné un âge jugé trop ancien (TO‑9124 : 3620 (+290/‑260) années cal. BP). Comme la courbe de Allard et Tremblay (1983a) a été construite à partir de nombreuses datations 14C dont les marges d’erreur altitudinales et temporelles ont été prises en compte, cette datation n’a pas été retenue pour la confection de la courbe d’émersion. Il est possible que le lit de matière organique échantillonné provienne en fait d’une laisse de mer.
Les datations Inv4, Inv5 et Inv6 situées entre 175 et 200 m d’altitude se trouvent considérablement en marge de la courbe (fig. 5). Les sols de ces plages enfouies montrent un développement avancé avec la présence d’un horizon Aej, suggérant l’existence d’un délai important entre l’émersion et l’enfouissement des surfaces. Ces datations 14C n’ont donc pas été prises en compte pour la confection de la courbe d’émersion.
Les quatre datations 14C obtenues à partir d’échantillons de coquilles de Mytilus edulis, un coquillage de la zone intertidale reconnu pour donner une bonne approximation de l’altitude des surfaces lors de l’émersion (Allard et Tremblay, 1983a ; Mitrovica et al., 2000), proviennent de plages situées à une altitude comprise entre 77 et 170 m (tabl. II ; fig. 5), mais elles ne sont pas significativement différentes (t = 4.31 ; p > 0,05). Le tracé de la courbe prend en compte la datation de Mytilus edulis obtenue sur la surface la plus élevée des quatre sites (170 m). Les trois autres sites localisés à 77, 122 et 134 m sont considérés ici comme montrant un âge plus ancien que celui de l’émersion des terres en raison d’un remaniement des dépôts de plage par des processus littoraux (vagues de tempête et poussées glacielles).
Discussion
Processus pédogénétiques
La migration vers le bas du profil de complexes organo-métalliques avec les eaux de percolation serait le principal mode d’enrichissement en fer et en aluminium de l’horizon B dans les sols sableux et bien drainés de la région du lac Guillaume-Delisle, considérant l’augmentation de la concentration en Corg, Fep et Alp dans l’horizon B (tabl. IV). L’augmentation des concentrations en Feo et Alo de la majorité des horizons B pourrait suggérer que la podzolisation se ferait aussi par l’enrichissement de l’horizon B sous forme de fer et d’aluminium amorphes ou au réseau cristallin peu structuré (ferrihydrites, proto-imogolite, imogolite et allophanes ; Wang, 1990 ; Childs, 1992). Cependant, l’augmentation en Feo et Alo reflète principalement l’enrichissement en complexes organo-métalliques, car les concentrations de (Feo-Fep) et (Alo-Ap) ne montrent pas de tendance à l’enrichissement dans l’horizon B.
Durée de la pédogenèse
L’étude d’une chronoséquence de sols repose sur le postulat selon lequel tous les facteurs de développement pédogénétique (matériel parental, climat, organismes et topographie) sont constants, sauf le temps (Jenny, 1941) qui devient alors le seul facteur responsable de l’évolution d’une propriété du sol. Cependant, les propriétés d’un sol sont plus fréquemment le résultat de l’action de plusieurs facteurs (Moore, 1978 ; Bockheim, 1980). Dans cette étude, les régressions linéaires multiples pas à pas qui portaient sur les variables caractérisant la pédogenèse des sols étudiés, soit le pH CaCl2 et les concentrations en Fep, Alp et Corg, ont toutefois montré que seul le temps permettait d’expliquer significativement une partie des variations observées et ce, uniquement pour la concentration en Alp.
La relation marginalement significative entre la concentration en Fep+Alp et la durée de la pédogenèse a tout de même permis de proposer une première approximation de la durée minimale pour l’obtention d’un horizon B podzolique répondant aux critères du Système canadien de classification des sols (Groupe de travail sur la classification des sols, 2002), soit environ 12 400 ans. Cette valeur demeure théorique, d’abord d’un point de vue statistique (p = 0,06), mais aussi parce qu’elle constitue une extrapolation au-delà des données échantillonnées (>4120 années cal. BP) et même au-delà de l’âge même des surfaces émergées (>8300 années cal. BP). D’ailleurs, la durée de pédogenèse obtenue est légèrement supérieure à celles qui ont été répertoriées jusqu’ici, les podzols répondant aux critères des différentes classifications nationales et internationales prenant habituellement entre 1000 et 10 000 années à se former (Muhs, 1984). Cependant, la durée estimée de 12 400 ans permet de souligner que les conditions requises pour la formation d’un horizon Bf n’ont pu, à toutes fins pratiques, être atteintes dans les sols subarctiques sableux, bien drainés et où la végétation reste basse et le couvert nival peu épais.
La durée minimale de 2300 années estimée par Protz et al. (1984) pour la formation d’un Podzol subarctique (côte ouest de la Baie d’Hudson) répondant aux critères du Système canadien de classification des sols (Groupe de travail sur la classification des sols, 2002) est nettement plus rapide qu’au lac Guillaume-Delisle. Les sols s’étaient aussi développés dans des sables et graviers (altitude variant de 4,2 à 104,7 m), mais sous une pessière à lichen où la quantité de neige au sol est plus abondante que sur les surfaces de toundra de la région du lac Guillaume-Delisle, ce qui pourrait expliquer la différence.
Payette (1973) a d’ailleurs montré que, sur la côte est de la Baie d’Hudson, les Podzols trouvent leur limite nordique un peu au nord de Kuujjuarapik (fig. 1) et sont confinés aux sites à fort enneigement et bien alimentés en eau. Payette et Filion (1993) ont aussi montré que la pédogenèse des sols podzoliques à plaques était dépendante de la présence d’un couvert nival associé à la présence de combes à neige ou d’une végétation haute captant la neige et favorisant la podzolisation lors de la fonte printanière. Les sols de la région du lac Guillaume-Delisle qui ont été retenus pour l’établissement de la chronoséquence sont caractérisés par une topographie relativement plane et par la présence d’une végétation basse qui ne contribue pas à la rétention de la neige. Lors de la fonte printanière, l’apport en eau est faible et ne permet pas une migration rapide des oxydes de fer et d’aluminium vers le bas du profil et, par conséquent, la formation d’un horizon Bf. Le fait que le drainage des sols soit rapide diminue encore davantage le temps de percolation de l’eau au sein du profil.
Le temps requis pour atteindre les seuils minimums exigés par les différentes classifications des sols, notamment celle du Système canadien (Groupe de travail sur la classification des sols, 2002) peut donc varier considérablement en fonction des conditions édaphiques (Mokma et al., 2004), notamment la température de l’air et du sol pendant la saison de croissance, les précipitations, de même que la topographie et la hauteur du couvert végétal qui contrôlent le régime d’enneigement (Payette et Filion, 1993b). Dans la région du lac Guillaume-Delisle, les conditions édaphiques prévalant à la formation des sols subarctiques bien drainés n’ont pas permis, depuis l’émersion des terres, d’obtenir de véritables Podzols.
Courbe d’émersion
Les datations 14C obtenues du domaine terrestre livrent un âge minimum pour l’émersion des terres et celles du domaine marin, un âge maximum. Les sols enfouis sont particulièrement utiles puisqu’ils représentent un période de stabilité géomorphologique permettant la colonisation végétale des dépôts de plages émergées. Les paléosols étudiés ne montrent pas de signes de remaniement par des vagues de tempêtes ou des poussées glacielles (à l’exception du sol S6). L’incertitude associée à l’altitude de coquillages qui ne sont pas échantillonnés en position de vie (Allard et Tremblay, 1983a) ne s’applique pas aux sols enfouis. Cependant, l’incertitude temporelle ne peut être exclue, ni pour le domaine marin (âge apparent local ou effet réservoir), ni pour le domaine terrestre (temps de séjour).
Les valeurs de pédogenèse des sols ont montré une relation significative avec la durée de la pédogenèse (temps) pour Alp. La plupart des paléosols peu développés ont été enfouis rapidement après l’émersion, de sorte qu’ils offrent d’excellentes possibilités pour l’établissement d’une courbe d’émersion des terres. En ce sens, les plages enfouies des profils S1 et N3 sont particulièrement importantes : les paléosols de ces plages ne montrent pratiquement aucun développement pédogénétique (régosols orthiques), ce qui implique un enfouissement rapide des plages suite à l’émersion. Le tracé de la courbe d’émersion proposée ici, et notamment dans son inflexion après ca. 4500 années cal. BP, prend en compte ces données pédogénétiques en ne laissant qu’une courte période de temps entre l’émersion et l’enfouissement des plages. À l’inverse, les plages des profils Inv4, Inv5 et Inv6 ont été enfouies longtemps après l’émersion et montrent des valeurs de pédogenèse plus élevées. Elles présentent donc un intérêt limité pour l’établissement de la courbe d’émersion.
Les deux courbes obtenues pour la région du lac Guillaume-Delisle sont similaires même si elles ont été obtenues par des méthodes différentes, c’est-à-dire une extrapolation à partir d’une périodicité moyenne de 45 ans pour la formation des cordons de plage (Hillaire-Marcel, 1976) et l’utilisation de datations 14C étalonnées provenant du domaine terrestre et du domaine marin (cette étude). Bien que seules des datations non corrigées aient été utilisées par Hillaire-Marcel (1976), le fait de mesurer le nombre de cordons entre des niveaux datés relativement rapprochés dans le temps a permis de limiter les effets des variations holocènes dans le flux de 14C atmosphérique sur la précision des datations au radiocarbone (Stuiver et al., 1998). Entre 5000 et 2500 années cal. BP, la courbe d’émersion révisée indique cependant une inflexion plus prononcée que ce qui avait été proposé par Hillaire-Marcel (1976), mais elle se situe tout de même dans l’intervalle établi pour le détroit de Manitounuk (Allard et Tremblay, 1983b). La courbe d’émersion établie à partir de la périodicité des plages repose sur quatre datations 14C provenant directement de la région du lac Guillaume-Delisle, toutes plus anciennes que 6000 années 14C BP, ce qui pourrait expliquer les différences avec notre courbe : les nouvelles données présentées ici permettent donc de préciser la courbe existante de l’Holocène moyen à aujourd’hui, par l’utilisation de la datation 14C de sols enfouis et de leur développement pédogénétique.
Conclusion
L’ensemble des datations 14C obtenues de sols et des paléosols a permis de montrer que, dans la région du lac Guillaume-Delisle, les sols sableux bien drainés se forment principalement par la migration de complexes organo-métalliques vers le bas du profil. La chronoséquence de sols subarctiques formés dans des cônes de déjection montre que des signes de podzolisation peuvent apparaître dès 570 ans (horizon Aej) après la mise en place des surfaces. La durée de pédogenèse requise pour obtenir le seuil minimal (Fep+Alp) de 0,4 % pour l’horizon Bf fixé par le Système canadien de classification des sols (Groupe de travail sur la classification des sols, 2002) est estimée à environ 12 400 ans, ce qui est plus élevé que les valeurs obtenues par Muhs (1984) et Protz et al. (1984) dans des environnements subarctiques. La datation d’horizons organiques enfouis et l’établissement de chronofonctions ont aussi permis de mettre à jour et de préciser la courbe d’émersion pour la côte ouest du lac Guillaume-Delisle pour la période allant de ca. 4500 ans à nos jours, ce qui constitue une contribution importante dans un contexte où les sédiments littoraux du secteur sont pauvres en coquillages.
Appendices
Remerciements
Les auteurs tiennent à remercier Stéphane Babin, Richard Dubuc, Yann Gilbert, Marianne Lafond, Michel Lajoie, Éva Stephani, Emmanuel Tousignant et le regretté Stéphane Beaudoin pour leur assistance sur le terrain et en laboratoire. Les commentaires de Michel Caillier, François Courchesne, Michel Allard, Hugo Asselin et Caroline Lavoie sur une version antérieure du manuscrit ont aussi été grandement appréciés, ainsi que ceux de Michel Nolin, Claude Hillaire-Marcel et d’un lecteur anonyme qui ont fait une lecture critique du manuscrit. Les analyses statistiques multivariées ont été réalisées par Hugo Asselin. Ce projet de recherche a été financé par le Fonds Québécois de la Recherche sur la Nature et les Technologies (FQRNT) et le Conseil de Recherche en Sciences Naturelles et en Génie du Canada (CRSNG). Les travaux de terrain ont été réalisés en partie grâce au soutien du Ministère des Affaires indiennes et du Nord du Canada.
Références
- Allard, M. et Séguin, M.K., 1987. The Holocene evolution of permafrost near the tree line, on the eastern coast of Hudson Bay (northern Quebec). Canadian Journal of Earth Sciences, 24 : 2206-2222.
- Allard, M. et Tremblay, G., 1983a. La dynamique littorale des îles Manitounuk durant l’Holocène. Zeitschrift für Geomorphologie, Supplement bänd 47 : 61-95.
- Allard, M. et Tremblay, G., 1983b. Les processus d’érosion littorale périglaciaire de la région de Poste-de-la-Baleine et des îles Manitounuk sur la côte est de la mer d’Hudson, Canada. Zeitschrift für Geomorphologie, Supplement bänd 47 : 27-60.
- Ballantyne, C.K. et Whittington, G., 1999. Late Holocene floodplain incision and alluvial fan formation in the central Grampian Highlands, Scotland : chronology, environment and implications. Journal of Quaternary Science, 14 : 651-671.
- Barber, D.C., Dyke, A., Hillaire-Marcel, C., Jennings, A.E., Andrews, J.T., Kerwin, M.W., Bilodeau, G., McNeely, R., Southon, J., Morehead, M.D. et Gagnon, J.-M., 1999. Forcing of the cold event of 8,200 years ago by catastrophic drainage of Laurentide lakes. Nature, 400 : 344-348.
- Baril, R. et Bitton, G., 1969. Teneurs élevées de fer libre et identification taxonomique de certains sols du Québec contenant de la magnétite. Canadian Journal of Soil Science, 49 : 1-9.
- Bégin, Y. et Payette, S., 1988. Dendroecological evidence of lake-level changes during the last three centuries in subarctic Québec. Quaternary Research, 30 : 210-220.
- Blair, T.C., 1987. Sedimentary processes, vertical stratification sequences, and geomorphology of the Roaring River alluvial fan, Rocky Mountain National Park, Colorado. Journal of Sedimentary Petrology, 57 : 1-18.
- Blair, T.C., 1999. Sedimentary processes and facies of the waterlaid Anvil Spring Canyon alluvial fan, Death Valley, California. Sedimentology, 46 : 913-940.
- Blake, W., Jr., 1988. Geological survey of Canada radiocarbon dates. Geological Survey of Canada, Ottawa, Paper 87-7, 100 p.
- Bockheim, J.G., 1980. Solution and use of chronofunctions in studying soil development. Geoderma, 24 : 71-85.
- Catt, J.A., 1991. Soils as indicators of Quaternary climatic change in mid-latitude regions. Geoderma, 51 : 167-187.
- Chandler, F.W. et Schwarz, E.J., 1980. Tectonics of the Richmond Gulf area, northern Quebec — A hypothesis. Geological Survey of Canada, Ottawa, Paper 80-1C, p. 59-68.
- Childs, C.W., 1992. Ferrihydrite : a review of structure, properties and occurence in relation to soils. Zeitschrift für Pflanzenernährung und Bodenkunde, 155 : 441-448.
- Environnement Canada, 2004. Normales et moyennes climatiques au Canada 1971-2000, Service météorologique du Canada. Disponible en ligne à http://www.climat.meteo.ec.gc.ca/climate_normals/index_f.html, page consultée le 15 décembre 2006.
- Filion, L. et Morisset, P., 1983. Eolian landforms along the eastern coast of Hudson Bay, Northern Québec. Nordicana, 47 : 73-94.
- Filion, L. et Payette, S., 1976. La dynamique de l’enneigement en région hémi-arctique, Poste-de-la-Baleine, Nouveau Québec. Cahiers de géographie de Québec, 20 : 275-302.
- Filion, L., Saint-Laurent, D., Desponts, M. et Payette, S., 1991. The late Holocene record of aeolian and fire activity in northern Québec, Canada. The Holocene, 1 : 201-208.
- Folk, R.L., 1974. Petrology of sedimentary rocks. Hemphill publishing company, Austin, 182 p.
- Folk, R.L. et Ward, W.C., 1957. Brazos river bar : a study in the significance of grain size parameters. Journal of Sedimentary Petrology, 27 : 3-26.
- Gray, J., Lauriol, B., Bruneau, D. et Ricard, J., 1993. Postglacial emergence of Ungava Peninsula, and its relationship to glacial history. Canadian Journal of Earth Sciences, 30 : 1676-1696.
- Groupe de travail sur la classification des sols, 2002. Le système canadien de classification des sols, 3e édition. Direction générale de la recherche, Ministère de l’Agriculture et Agro-Alimentaire du Canada, Ottawa, 196 p.
- Harvey, A.M., 1986. Geomorphic effects of a 100 year storm in the Howgill Fells, Northwest England. Zeitscrift für Geomorphologie, 30 : 71-91.
- Head, K.H., 1992. Manual of soil laboratory testing, 2e édition. Pentech Press, London, 388 p.
- Hillaire-Marcel, C., 1976. La déglaciation et le relèvement isostatique sur la côte est de la baie d’Hudson. Cahiers de géographie de Québec, 20 : 185-220.
- Hughen, K.A., Baillie, M.G.L., Bard, E., Bayliss, A., Beck, J.W., Bertrand, C.J.H., Blackwell, P.G., Buck, C.E., Burr, G.S., Cutler, K.B., Damon, P.E., Edwards, R.L., Fairbanks, R.G., Friedrich, M., Guilderson, T.P., Kromer, B., McCormac, F.G., Manning, S.W., Bronk Ramsey, C., Reimer, P.J., Reimer, R.W., Remmele, S., Southon, J.R., Stuiver, M., Talamo, S., Taylor, F.W., van der Plicht, J. et Weyhenmeyer, C.E., 2004. Marine04 Marine radiocarbon age calibration, 26-0 ka BP. Radiocarbon, 46 : 1059-1086.
- Jenny, H., 1941. Factors of soil formation and system of quantitative pedology. McGraw-Hill, London, 281 p.
- Kinnard, C. et Lewkowicz, A.G., 2006. Frontal advance of turf-banked solifluction lobes, Kluane Range, Yukon Territory, Canada. Geomorphology, 73 : 261-276.
- Lafortune, V., Filion, L. et Hétu, B., 2006. Impacts of Holocene climatic variations on alluvial fan activity below snowpatches in subarctic Québec. Geomorphology, 76 : 375-391.
- Lajeunesse, P. et Allard, M., 2003a. The Nastapoka drift belt, eastern Hudson Bay : implications of a stillstand of the Quebec-Labrador ice margin in the Tyrrell Sea at 8 ka BP. Canadian Journal of Earth Sciences, 40 : 65-76.
- Lajeunesse, P. et Allard, M., 2003b. Late Quaternary deglaciation, glaciomarine sedimentation and glacioisostatic recovery in the Rivière Nastapoka area, eastern Hudson Bay, northern Québec. Géographie physique et Quaternaire, 57 : 65-83.
- Martin, C.W. et Johnson, W.C., 1995. Variation in radiocarbon ages of soil organic matter fractions from Late Quaternary buried soils. Quaternary Research, 43 : 232-237.
- Matthews, J.A., 1980. Some problems and implications of 14C dates from a podzol buried beneath an end moraine at Haugabreen, Southern Norway. Geografiska Annaler, 62A : 185-208.
- McKeague, J.A., 1978. Manuel de méthodes d’échantillonnage et d’analyse des sols, 2e édition. Canadian Society of Soil Science, Ottawa, 250 p.
- McKeague, J.A., Brydon, J.E. et Miles, N.M., 1971. Differentiation of forms of extractable iron and aluminum in soils. Soil Science Society of America Proceedings, 35 : 33-38.
- McKeague, J.A. et Day, J.H., 1966. Dithionite- and oxalate-extractable Fe and Al as aids in differentiating various classes of soils. Canadian Journal of Soil Science, 46 : 13-22.
- Mitrovica, J.X., Forte, A.M. et Simons, M., 2000. A reappraisal of postglacial decay times from Richmond Gulf and James Bay, Canada. Geophysical Journal International, 142 : 783-800.
- Mokma, D., Yli-Halla, M. et Lindqvist, K., 2004. Podzol formation in sandy soils of Finland. Geoderma, 120 : 259-272.
- Moore, T.R., 1978. Soil formation in Northeastern Canada. Annals of the Association of American Geographers, 68 : 518-534.
- Muhs, D.R., 1984. Intrinsic thresholds in soil systems. Physical Geography, 5 : 99-110.
- Paul, E.A., Campbell, C.A., Rennie, D.A. et McCallum, K.J., 1964. Investigations of the dynamics of soil humus utilizing carbon dating techniques. Transactions of the 8th International Congress of Soil Science, 3 : 201-208.
- Payette, S., 1973. Contribution à la pédologie de la zone hémi-arctique : région de Poste-de-la-baleine, Nouveau-Québec. Le Naturaliste canadien, 100 : 123-163.
- Payette, S., 1975. La limite septentrionale des forêts sur la côte orientale de la baie d’Hudson, Nouveau Québec. Le Naturaliste canadien, 102 : 317-329.
- Payette, S., 1976. Les limites écologiques de la zone hémi-arctique entre la mer d’Hudson et la baie d’Ungava, Nouveau-Québec. Cahiers de géographie du Québec, 20 : 347-364.
- Payette, S. et Filion, L., 1975. Écologie de la limite septentrionale des forêts maritimes, Baie d’Hudson, Nouveau-Québec. Le Naturaliste canadien, 102 : 783-802.
- Payette, S. et Filion, L., 1993a. Holocene water-level fluctuations of a subarctic lake at the tree line in northern Québec. Boreas, 22 : 7-14.
- Payette, S. et Filion, L., 1993b. Origin and significance of subarctic patchy podzolic soils and paleosols. Arctic and Alpine Research, 25 : 267-276.
- Payette, S., Filion, L., Gauthier, L. et Boutin, Y., 1985. Secular climate change in old-growth tree-line vegetation of northern Quebec. Nature, 315 : 135-138.
- Protz, R., Ross, G. J., Martini, I. P. et Terasmae, J., 1984. Rate of podzolic soil formation near Hudson Bay, Ontario. Canadian Journal of Soil Science, 64 : 31-49.
- Reimer, P., 2005. Marine reservoir correction database. Disponible en ligne à http://calib.qub.ac.uk/marine/, page consultée le 21 avril 2006.
- Reimer, P.J., Baillie, M.G.L., Bard, E., Bayliss, A., Beck, J.W., Bertrand, C.J.H., Blackwell, P.G., Buck, C.E., Burr, G.S., Cutler, K.B., Damon, P.E., Edwards, R.L., Fairbanks, R.G., Friedrich, M., Guilderson, T.P., Hogg, A.G., Hughen, K.A., Kromer, B., McCormac, F.G., Manning, S.W., Ramsey, C.B., Reimer, R.W., Remmele, S., Southon, J.R., Stuiver, M., Talamo, S., Taylor, F.W., van der Plicht, J. et Weyhenmeyer, C.E., 2004. IntCal04 Terrestrial radiocarbon age calibration, 26-0 ka BP. Radiocarbon, 46 : 1029-1058.
- Stork, A., 1963. Plant immigration in front of retreating glaciers, with examples from the Kebnekajse area, Northern Sweden. Geografiska Annaler, 45A : 1-21.
- Stuiver, M., Reimer, P.J., Bard, E., Beck, J.W., Burr, G.S., Hughen, K.A., Kromer, B., McCormac, F.G., van der Plicht, J. et Spurk, M., 1998. INTCAL98 Radiocarbon age calibration 24,000-0 cal BP. Radiocarbon, 40 : 1041-1083.
- Stuiver, M., Reimer, P.J. et Reimer, R., 2005. CALIB 5.0.1. Disponible en ligne à http://www.calib.org, page consultée le 21 avril 2006.
- Telford, R.J., Heegaard, E. et Birks, H.J.B., 2004. The intercept is a poor estimate of a calibrated radiocarbon age. The Holocene, 14 : 296-298.
- Ugolini, F.C., 1986. Pedogenic zonation in the well-drained soils of the arctic regions. Quaternary Research, 26 : 100-120.
- Ugolini, F.C., Stoner, M.G. et Marrett, D.J., 1987. Arctic pedogenesis : 1. Evidence for contemporary podzolization. Soil Science, 144 : 90-100.
- Viereck, L.A., 1966. Plant succession and soil development on gravel outwash of the Muldrow Glacier, Alaska. Ecological Monographs, 36 : 181-199.
- Viereck, L.A., 1970. Forest succession and soil development adjacent to the Chena River in interior Alaska. Arctic and Alpine Research, 2 : 1-26.
- Vreeken, W.J., 1975. Principal kinds of chronosequences and their significance in soil history. Journal of Soil Science, 26 : 378-394.
- Walker, A.L., 1983. The effects of magnetite on oxalate- and ditionite-extractable iron. Soil Science Society of America Journal, 47 : 1022-1026.
- Walkley, A. et Black, I.A., 1934. An examination of Degtjareff method for determining soil organic matter and a proposed modification of the chromic acid titration method. Soil Science, 37 : 29-37.
- Wang, C., 1990. The Canadian approach to identifying the spodic horizon, p. 387-394. In J.M. Kimble et R.D. Yeck, édit., Fifth international soil correlation meeting (ISCOM V) : characterization, classification, and utilisation of Spodosols. United States Department of Agriculture, Soil Conservation Service, Lincoln.
- Wang, Y. et Hsieh, Y.P., 2002. Uncertainties and novel prospects in the study of the soil carbon dynamics. Chemosphere, 49 : 791-804.
- Wells, S.G. et Harvey, A.M., 1987. Sedimentologic and geomorphic variations in storm-generated alluvial fans, Howgill Fells, northwest England. Geological Society of America Bulletin, 98 : 182-198.
- Woodcock, J.R., 1960. Geology of the Richmond Gulf area, New Quebec. Proceedings of the Geological Association of Canada, 12 : 21-39.
List of figures
Localisation de la région du lac Guillaume-Delisle et des sites d’étude. Les noms de lieux mentionnés dans le texte sont également indiqués.
Location of the Lake Guillaume-Delisle area and of the two study sites. Place names mentioned in the text are also indicated.
Blocs diagrammes montrant les diverses étapes de la formation des cônes de déjection dans la vallée de l’Amphithéâtre en relation avec le relèvement isostatique postglaciaire.
3D sketch showing various stages of alluvial fan development in the Amphithéâtre valley in relationship with postglacial land emersion.
Durée de la pédogenèse et propriétés chimiques des sols sableux bien drainés et exempts de l’influence de combes à neige dans les vallées de l’Amphithéâtre et de l’Inverbule. (A) Corg, (B) Fep, (C) Alp, (D) Fep+Alp et (E) pH CaCl2. L’échelle de l’ordonnée diffère en (A) et (E).
Time and soil chemical properties of well-drained soils not located in late-lying snowpatches environments in the Amphithéâtre and Inverbule valleys. (A) Corg, (B) Fep, (C) Alp, (D) Fep+Alp and (E) pH CaCl2. The ordinate scale differs in (A) and (E).