Numerical simulation of groundwater flow in regional rock aquifers, southwestern Quebec, Canada

Abstract

The St. Lawrence Lowlands platform, Quebec, Canada, is a densely-populated area, heavily dependent on groundwater resources. In 1999, the Geological Survey of Canada initiated a large-scale hydrogeological assessment study over a 1,500 km2 region northwest of Montreal. The objectives were to define the regional groundwater flow, and to give quantitative estimates of the groundwater dynamic parameters and of the available groundwater resources. The applied approach consisted of defining the hydrogeologic framework, hydraulic properties of the aquifer units, and groundwater dynamic components. Lower Paleozoic sedimentary rocks represent regional aquifer units. Coarse Quaternary fluvio-glacial sediments locally overlay the rock sequence and constitute an interface aquifer unit. Fine marine sediments confine most of the regional aquifers. Collected GIS based information was synthesized in a finite element numerical model. The regional saturated steady-state flow was calibrated under current stress conditions assuming an equivalent porous medium approach. Water budget calculations show that the total groundwater flow in regional aquifers amounts to 97.7 Mm3/y. Infiltration from precipitation provides 86.6% of the groundwater supply, while 9.6% comes from subsurface inflow and the remaining 3.8% is induced recharge from surface waters. Discharge from regional aquifers occurs through flow to streams (76.9%), groundwater withdrawal (18.4%), and underground outflow (4.7%).