Coupled cryo-hydrogeological modelling of permafrost dynamics near Umiujaq (Nunavik, Canada)

A two-dimensional (2D) cryo-hydrogeological numerical model of groundwater flow, coupled with advective-conductive heat transport with phase change, has been developed to study permafrost dynamics around an ice-rich permafrost mound in the Tasiapik Valley near Umiujaq, Nunavik (Quebec), Canada. Permafrost is degrading in this valley due to climate warming observed in Nunavik over the last two decades. Ground temperatures measured along thermistor cables in the permafrost mound show that permafrost thaw is occurring both at the permafrost table and base, and that heat fluxes at the permafrost base are up to ten times higher than the expected geothermal heat flux. Based on a vertical cross-section extracted from a 3D geological model of the valley, the numerical model was first calibrated using observed temperatures and heat fluxes. Comparing simulations with and without groundwater flow, advective heat transport due to groundwater flow in the subpermafrost aquifer is shown to play a critical role in permafrost dynamics and can explain the high apparent heat flux at the permafrost base. Advective heat transport leads to warmer subsurface temperatures in the recharge area, while the cooled groundwater arriving in the downgradient discharge zone maintains cooler temperatures than those resulting from thermal conduction alone. Predictive simulations incorporating a regional climate-change scenario suggest the active layer thickness will increase over the coming decades by about 12 cm/year, while the depth to the permafrost base will decrease by about 80 cm/year. Permafrost within the valley is predicted to completely thaw by around 2040.