Simulating Active Layer Thaw in a Boreal Environment

Abstract

A large part of the boreal zone of the western Canadian Arctic is underlain by ice-rich discontinuous permafrost which when thawed, can lead to settlement of the ground surface that has implications for the integrity of northern infrastructure, including oil and gas pipelines. A simple yet physically-based model is desired to simulate thawing of the active layer in different materials commonly found along the Mackenzie Valley pipeline corridor. Stefan’s algorithm determines the phase change of soil moisture using ground surface temperature as the upper boundary condition and conduction to transfer heat to the freeze-thaw front. It is tested on a permafrost site near Wrigley, Northwest Territories, where the computed thaw penetration compares satisfactorily with field data. To further explore the effects of climate and soil types on active layer depth, three representative sites in the Mackenzie valley where ground surface temperatures are available were selected for simulation of ground thaw, under two summer conditions. Results of the simulation demonstrate the sensitivity of active layer thaw to (1) soil materials due to differential thermal properties, (2) moisture content, which largely controls the latent heat requirement for phase change, and (3) inter-annual variations in ground surface temperature. Given the strong potential for environmental changes in the vast boreal region, the model allows the active layer thaw responses to be easily assessed.

Résumé

Simulation de la couche active de dégel d’un environnement boréal. Une vaste partie de la zone boréale de l’Arctique canadien occidental repose sur un pergélisol discontinu riche en glace qui, en dégelant, peut causer des problèmes de surface du sol pouvant compromettre l’intégrité des structures nordiques, dont les pipelines de pétrole et de gaz. Un modèle physique simple est nécessaire pour la simulation du dégel de la couche active des divers dépôts meubles trouvés le long du corridor du pipeline de la vallée du Mackenzie. L’algorithme de Stefan détermine la phase du changement d’humidité du sol en utilisant la température de surface comme la condition limite supérieure, et la conduction au transfert de la chaleur jusqu’au front de gel-dégel. Cet algorithme est appliqué à un site de pergélisol près de Wrigley, Territoires-du-Nord-Ouest, où la simulation de la profondeur de dégel se compare bien aux données de terrain. Afin d’explorer plus encore les effets du climat et du type de sol sur la profondeur de la couche active, trois sites représentatifs de la vallée du Mackenzie, où les températures de la surface du sol sont disponibles, ont été sélectionnés pour une simulation de dégel du sol, selon deux conditions estivales. Les résultats de la simulation montrent la sensibilité de la couche active de dégel (1) aux propriétés thermiques différentielles des matériaux du sol, (2) au contenu en humidité, qui contrôle la chaleur latente provoquant le changement de phase, et (3) aux variations interannuelles des températures de la surface du sol. En raison du fort potentiel de changements environnementaux dans la vaste région boréale, le modèle permet d’évaluer aisément le dégel de la couche active qui en constitue une réaction.