Simulating groundwater-peatland interactions in depression and slope peatlands in southern Quebec (Canada)

It is crucial to understand hydrogeological interactions between aquifers and peatlands in order to grasp the influence of aquifers in peatland water budgets, to understand the role of groundwater in the evolution or organic matter deposition, and to quantify how a peatland can sustain groundwater levels in a superficial aquifer. These questions have rarely been addressed in literature and there is currently no understanding of which process dominates aquifer-peatland exchanges in different geomorphological settings. The main purpose of the study was to use groundwater flow modeling to answer these questions in two contrasted geological contexts of southern Quebec (Canada). During a three-year study, six peatlands have been instrumented in the Becancour (Centre-du-Quebec) and Amos (Abitibi-Temiscamingue) regions of southern Quebec (Canada). At each site, either one or two transects of six piezometer nests (at 1.20 m depth in the organic deposits and in the mineral deposits below the peat) have been installed, for a total of twelve aquifer-peatland transects of approximately 500 m. The stratigraphy and geometry of the peatland-aquifer system, as well as the hydrodynamic properties of the organic and mineral deposits have been measured at all sites. Groundwater levels have been recorded from autumn 2010 to summer 2012. The Becancour peatlands have developed in depressions while the Amos peatlands have developed through the paludification of esker slopes. The maximum peat thickness measured in the Bécancour peatlands is 6.4 m while it is 4.5 m in the Amos region. In both regions, peatlands are fringed by sandy deposits that extend at least partly under the organic deposits. The thickness of these underlying deposits is not well defined, but available data suggests a metric scale thickness in areas close to the adjacent superficial aquifer. Field data is used to create 2D numerical models in Modflow to simulate flow between the shallow groundwater and the peatland on four peatland transects considered representative of the overall variability observed at the field sites. The models are first calibrated to reproduce measured heads, head gradients and temporal variations. In order to assess typical flow patterns and exchanges, a global sensitivity analysis of the model are performed to identify which parameters and processes control the exchanged fluxes. Results show that for depression peatlands, exchanged aquifer-peatland fluxes occur on short distances near the peatland border. For slope peatlands, exchanged fluxes are distributed further inside the peatland. Local hydrostratigraphy as well as peat and mineral deposits hydraulic properties control aquifer-peatland exchanges. Peat recharge is a challenge to represent, but appears to have a similar effect on the four simulated peatlands. Based on the influence of each parameter on the flow, a graphical tool is proposed to help estimate the exchanges between groundwater and peatlands when limited data is available.