Impact of the Wisconsinian Glaciation on Canadian Continental Groundwater Flow

• Main Author: Lemieux, Jean-Michel
• Format: Doctoral or Postdoctoral Thesis
• Language: English
• Published: University of Waterloo 2006
• Subjects: Wisconsinian glaciation; groundwater; climate change; numerical modelling; Earth Sciences; Ice; Ice Sheet; permafrost
• Online Access: http://hdl.handle.net/10012/2654

During the Quaternary period, cyclic glaciations have occurred over a global scale as the result of a climatic variability that affected the Earth's atmospheric, oceanic and glacial systems. Quaternary glaciations and their associated dramatic climatic conditions, such as kilometers-thick ice sheet formation and permafrost migration, are suspected to have had a large impact on the groundwater flow system over the entire North American continent. Because of the myriad of complex flow-related processes involved during a glaciation period, numerical models have become powerful tools to examine groundwater flow system evolution in this context. In this study, a series of key processes pertaining to coupled groundwater flow and glaciation modelling, such as density-dependent (i.e., brine) flow, hydromechanical loading, subglacial infiltration, isostasy, sea-level change and permafrost development, are included in the numerical model HydroGeoSphere to simulate groundwater flow over the Canadian landscape during the Wisconsinian glaciation (~ -120 kyr to present). The primary objective is to demonstrate the immense impact caused by glacial advances and retreats during the Wisconsinian glaciation on the dynamical evolution of groundwater flow systems over the Canadian landscape, including surface/subsurface water exchanges (i.e., recharge and discharge fluxes) both in the subglacial and the periglacial environments. The major findings of this study are that subglacial meltwater infiltration into the subsurface dominates when the ice sheet is growing and, conversely, groundwater exfiltrates during ice sheet regression. This conclusion, which seems to be opposite to the classical hydromechanical loading theory, is a consequence of the interaction between the subglacial boundary conditions and the elastic properties of the rocks. Subglacial infiltration rates during ice sheet progression can reach up to three orders of magnitude higher than the infiltration rate into the periglacial environment and the current recharge ...