A chemical and thermal modelling study of the active layer on Herschel Island, Yukon Territory

An essential requirement in evaluating the potential impacts of climate warming in Polar Regions, is a greater understanding of the nature and behavior of the active layer and the permafrost system. Numerous studies have documented the inherent capacity of permafrost to limit significant degradation from seasonal thawing because of the ability of overlying soil layers to dissipate influxes of heat. However, information on the degree to which these layers are able to act as buffers to permafrost degradation is limited, and thus further research on this system is needed.The research presented in this thesis in the form of two manuscripts. The first manuscript (Chapter 3) discusses the role of chemical processes in the soil weathering regime on Herschel Island. Three distinct cryostratigraphic layers were observed in the soil profiles, including; (1) the modern active-layer located immediately below the ground surface (0 – 39 cm); (2) a paleo-active (transition layer), extending from the bottom of the modern active layer to the top of an undisturbed permafrost layer (39 – 192 cm); and (3) the underlying permafrost. These layers are marked by two thaw unconformities that separate the modern active layer from the transition layer and the transition layer from the permafrost layer, respectively. The two thaw unconformities were distinguished by abrupt breaks in the cryogenic texture and peaks in the concentrations of several chemical components. The upper thaw unconformity, located at the base of the modern active layer, coincided with a doubling in the concentration of major soluble cations and a tripling in the organic matter content compared to those of the overlying soil horizons. The lower thaw unconformity is marked by the highest concentration of Al2O3 and K2O in the entire profile and the second highest concentration of Na2O. Between these two unconformities, (i.e. in the transition zone), soluble cation concentrations varied considerably, the bulk soil SiO2 concentration increased significantly and there was a noticeable decrease in the concentrations of Al2O3, Fe2O3, CaO, MgO and K2O in the depth interval of 82-107 cm. Chemical weathering indices recorded an overall decrease in the intensity of chemical weathering with depth; the soil layers immediately above the first thaw unconformity were the most weathered and the deeper permafrost layer was the least weathered.The second manuscript (Chapter 4) employs a heat conduction algorithm (Stefan equation) based on the composition and thermal properties of the soils discussed in Chapter 3, and measured active layer depth values, to evaluate heat transfer from the ground surface to the underlying permafrost. Thermal constants (the square root of the ratio of twice the thermal conductivity of the unfrozen soil to the latent heat of fusion of ice) were calculated using the De Vries method and were compared to empirical thermal constants derived from measurements of active layer depths and the corresponding thawing degree days using the Stefan equation. The theoretical constants are approximately 38% lower than those derived empirically, a difference that is likely due to the fact that the Stefan equation does not consider advective heat transfer and that the active layer depths were probably overestimated because of the method of measurement (resistance to penetration of a metal probe). Third generation Canadian Global Coupled Models (CGCM 3.1/T63) from the Canadian Centre for Climate Modelling and Analysis (CCCma), based on the IPCC SRES A2 and B1 scenarios, were used in conjunction with the theoretically and empirically derived thermal constants to predict active layer depths for 2050 and 2100. These calculations suggest that the active layer will thicken by up to 41.76 % before the end of this century. This thesis provides new insights into the important role that the active and transition layers play in controlling the response of permafrost systems to climate change.Une condition essentielle pour évaluer les impacts potentiels du réchauffement climatique dans les régions polaires est l'acquisition d'une meilleure compréhension de la nature et du comportement de la couche active du pergélisol. La recherche présentée dans cette thèse est présentés sous la forme de deux manuscrits. Le premier manuscrit examine le rôle des processus chimiques dans le régime d'altération du sol sur l'île Herschel. Trois couches distinctes cryostratigraphiques ont été observées dans les profils de sol: (1) la couche active moderne située immédiatement sous la surface du sol (0 - 39 cm), (2) la couche paléo-active (couche de transition) s'étendant de la base de la couche moderne active à la partie supérieure de la couche de pergélisol non perturbée (39 - 192 cm) et (3) le pergélisol sous-jacent. Ces couches sont marquées par deux discordances de dégel séparant la couche active moderne de la couche de transition et la couche de transition de la couche de pergélisol, respectivement. Ces deux discordances de dégel ont été identifiées grâce à la présence de ruptures brutales dans la texture cryogénique et par des pics dans la concentration de plusieurs composants chimiques. La discordance de dégel supérieure, située à la base de la couche moderne active, coïncide avec un doublement de la concentration des principaux cations solubles et un triplement de la teneur en matière organique par rapport aux horizons sus-jacents. La discordance de dégel inférieure est marquée par la plus forte concentration en Al2O3 et en K2O de l'ensemble du profil et par la deuxième plus forte concentration en Na2O. Entre ces deux discordances, les concentrations en cations solubles varient considérablement, la concentration en SiO2 du sol brut s'accroit de façon significative et on observe une diminution notable des concentrations en Al2O3, Fe2O3, CaO, MgO et K2O dans l'intervalle de profondeur de 82-107 cm. Les indices d'altération chimique ont enregistré une diminution globale de l'intensité de l'altération chimique avec la profondeur ; les couches de sol immédiatement au-dessus de la première discordance de dégel sont les plus altérées et la couche profonde du pergélisol est la moins altérée.Le second manuscrit utilise un algorithme de conduction de la chaleur basé sur la composition et les propriétés thermiques des sols examinés au chapitre 3 et sur les mesures des valeurs de profondeur de la couche active afin d'évaluer le transfert de chaleur de la surface du sol au pergélisol sous-jacent. Des constantes thermiques (la racine carrée du rapport entre deux fois la conductivité thermique du sol non gelé et la chaleur latente de fusion de la glace) ont été calculées en utilisant la méthode de De Vries et comparées aux constantes thermiques empiriques issues des mesures de profondeur de la couche active et des degrés-jours de dégel correspondant en utilisant l'équation de Stefan. Les constantes théoriques sont environ 38% inférieures à celles obtenues empiriquement, une différence qui s'explique probablement par le fait que l'équation de Stefan ne considère pas le transfert de chaleur par advection et que les profondeurs de la couche active ont probablement été surestimées en raison de la méthode de mesure (résistance à la pénétration d'une sonde métallique). La troisième génération du modèle couplé climatique global canadien (CGCM 3.1/T63) du Centre canadien de la modélisation et de l'analyse climatique (CCmaC), basée sur les scénarios A2 et B1 du IPCC SRES, a été utilisée en conjonction avec les constantes thermiques théoriques et empiriques pour prédire les profondeurs de la couche active pour 2050 et 2100. Ces calculs suggèrent que la couche active s'épaissira de jusqu'à 41,76% avant la fin du siècle.Ce mémoire apporte un nouvel éclairage sur le rôle important que les couches actives et de transition jouent dans le contrôle de la réponse des systèmes de pergélisol au changement climatique