Marie Larocque

Methane Baseline Concentrations and Sources in Shallow Aquifers from the Shale Gas-Prone Region of the St. Lawrence Lowlands (Quebec, Canada)

Hydraulic fracturing is becoming an important technique worldwide to recover hydrocarbons from unconventional sources such as shale gas. In Quebec (Canada), the Utica Shale has been identified as having unconventional gas production potential. However, there has been a moratorium on shale gas exploration since 2010. The work reported here was aimed at defining baseline concentrations of methane in shallow aquifers of the St. Lawrence Lowlands and its sources using δ(13)C methane signatures. Since this study was performed prior to large-scale fracturing activities, it provides background data prior to the eventual exploitation of shale gas through hydraulic fracturing. Groundwater was sampled from private (n = 81), municipal (n = 34), and observation (n = 15) wells between August 2012 and May 2013. Methane was detected in 80% of the wells with an average concentration of 3.8 ± 8.8 mg/L, and a range of <0.0006 to 45.9 mg/L. Methane concentrations were linked to groundwater chemistry and distance to the major faults in the studied area. The methane δ(1)(3)C signature of 19 samples was > -50‰, indicating a potential thermogenic source. Localized areas of high methane concentrations from predominantly biogenic sources were found throughout the study area. In several samples, mixing, migration, and oxidation processes likely affected the chemical and isotopic composition of the gases, making it difficult to pinpoint their origin. Energy companies should respect a safe distance from major natural faults in the bedrock when planning the localization of hydraulic fracturation activities to minimize the risk of contaminating the surrounding groundwater since natural faults are likely to be a preferential migration pathway for methane.

Estimating Flow Using Tracers and Hydraulics in Synthetic Heterogeneous Aquifers

Regional ground water flow is most usually estimated using Darcys law, with hydraulic conductivities estimated from pumping tests, but can also be estimated using ground water residence times derived from radioactive tracers. The two methods agree reasonably well in relatively homogeneous aquifers but it is not clear which is likely to produce more reliable estimates of ground water flow rates in heterogeneous systems. The aim of this paper is to compare bias and uncertainty of tracer and hydraulic approaches to assess ground water flow in heterogeneous aquifers. Synthetic two-dimensional aquifers with different levels of heterogeneity (correlation lengths, variances) are used to simulate ground water flow, pumping tests, and transport of radioactive tracers. Results show that bias and uncertainty of flow rates increase with the variance of the hydraulic conductivity for both methods. The bias resulting from the nonlinearity of the concentration-time relationship can be reduced by choosing a tracer with a decay rate similar to the mean ground water residence time. The bias on flow rates estimated from pumping tests is reduced when performing long duration tests. The uncertainty on ground water flow is minimized when the sampling volume is large compared to the correlation length. For tracers, the uncertainty is related to the ratio of correlation length to the distance between sampling wells. For pumping tests, it is related to the ratio of correlation length to the pumping tests radius of influence. In regional systems, it may be easier to minimize this ratio for tracers than for pumping tests.

Groundwater recharge trends in CanadaGeological Survey of Canada Contribution 20090009.

Groundwater plays a major role in social and economic development and in human and ecosystem health. However, little is known about the potential impacts of climate change on this resource in Canada, namely if groundwater recharge is increasing or decreasing over time. This paper focuses on trend statistical analysis of historical series of baseflow and groundwater levels and their field significance as indicators of recharge. Monitoring wells are mainly located in the southern half of western Canada, where few gauging stations either are available or provide significant trends. Both data sets are thus complementary. Results show that most available groundwater level series have significant trends (80%), whereas most available baseflow series have not (3%–33%). However, groundwater level series usually show smaller slope magnitudes than baseflow series. Mixed trends are often observed across Canada for a given variable, period, or series length, although some regions can have marked trends. For instance, ...

Environmental factors explaining the vegetation patterns in a temperate peatland

Although ombrotrophic temperate peatlands are important ecosystems for maintaining biodiversity in eastern North America, the environmental factors influencing their flora are only partly understood. The relationships between plant species distribution and environmental factors were thus studied within the oldest temperate peatland of Québec. Plant assemblages were identified by cluster analysis while CCA was used to related vegetation gradients to environmental factors. Five assemblages were identified; three typical of open bog and two characterized by more minerotrophic vegetation. Thicker peat deposit was encounter underlying the bog assemblages while higher water table level and percentage of free surface water distinguished the minerotrophic assemblages. Overall, the floristic patterns observed were spatially structured along the margins and the expanse. The most important environmental factors explaining this spatial gradient were the percentage of free surface water and the highest water-table level.

On the Variable-Lag and Variable-Velocity Cell-to-Cell Routing Schemes for Climate Models

Abstract The cell‐to‐cell channel routing schemes used in General Circulation Models (GCMs) and Regional Climate Models (RCMs) are revisited. A simpler parsimonious routing scheme based on Askews formula (1970) for computing time‐evolving channel lags is implemented and tested against observations and compared with the variable‐velocity scheme of Arora and Boer (1999). The variable‐lag routing scheme agrees very well with the variable‐velocity scheme. The variable‐lag scheme has the advantage of using fewer parameters, which is a major advantage at fine resolution over a large domain, where the uncertainty associated with parameters can be quite large. The spatial resolutions of RCMs are much finer than those of GCMs and hence there is a need for channel routing at fine spatial resolutions. The task of extending the cell‐to‐cell routing schemes developed for large‐scale routing, as in GCMs, to finer spatial scales, as in RCMs, is addressed. The sensitivity of the variable‐lag scheme to the routing time interval is studied. The choice of the routing time interval is very critical and varies with the spatial resolution as in any hydrological model. A simple method for determining the appropriate range of routing intervals at different spatial resolutions for the variable‐lag scheme is presented.

The potential of major ion chemistry to assess groundwater vulnerability of a regional aquifer in southern Quebec (Canada)

Groundwater vulnerability mapping provides useful but limited information for developing protection plans of the resource. Classical vulnerability ranking methods often do not take into account complex hydrostratigraphy and never consider groundwater flow dynamics. The objective of this work was to test the potential of major ion chemistry to assess regional-scale intrinsic groundwater vulnerability. Because it reflects water–sediment and water–rock interactions, the new vulnerability index reflects both infiltration processes and groundwater hydrodynamics. The method was applied on a regional fractured bedrock aquifer located in the Becancour region of southern Quebec (Canada). In this region, hydrogeochemistry shows that freshly recharged groundwater evolves from (Ca, Mg)–HCO3 and Ca–SO4 to Na–HCO3 type with gradually increasing confinement conditions in the fractured aquifer and tends to Na–Cl type locally by mixing with trapped marine pore-water. The new method identified recharge areas as those of the highest vulnerability and gradually decreasing vulnerability as confinement of the aquifer increased. It also highlights local discontinuities in confinement that differ from the regional pattern. Results showed a good correlation between groundwater vulnerability estimated with the new method and nitrate occurrence in groundwater. Eighty-two per cent of all samples presenting detectable nitrate concentrations were characterized by a Hydrogeochemical Vulnerability Index greater than 9 (maximum is 10). The ability of the new vulnerability method to identify areas vulnerable to detectable nitrate concentrations was much higher than that deriving from the DRASTIC method. This work confirms that major ions chemistry contains significant information about groundwater vulnerability and could be used to improve groundwater resource management.

222Rn activity in groundwater of the St. Lawrence Lowlands, Quebec, eastern Canada: relation with local geology and health hazard

One hundred ninety-eight groundwater wells were sampled to measure the (222)Rn activity in the region between Montreal and Quebec City, eastern Canada. The aim of this study was to relate the spatial distribution of (222)Rn activity to the geology and the hydrogeology of the study area and to estimate the potential health risks associated with (222)Rn in the most populated area of the Province of Quebec. Most of the groundwater samples show low (222)Rn activities with a median value of 8.6 Bq/L. Ninety percent of samples show (222)Rn activity lower than 100 Bq/L, the exposure limit in groundwater recommended by the World Health Organization. A few higher (222)Rn activities (up to 310 Bq/L) have been measured in wells from the Appalachian Mountains and from the magmatic intrusion of Mont-Saint-Hilaire, known for its high level of indoor radon. The spatial distribution of (222)Rn activity seems to be related mainly to lithology differences between U-richer metasediments of the Appalachian Mountains and magmatic intrusions and the carbonaceous silty shales of the St. Lawrence Platform. Radon is slightly enriched in sodium-chlorine waters that evolved at contact with clay-rich formations. (226)Ra, the parent element of (222)Rn could be easily adsorbed on clays, creating a favorable environment for the production and release of (222)Rn into groundwater. The contribution of groundwater radon to indoor radon or by ingestion is minimal except for specific areas near Mont-Saint-Hilaire or in the Appalachian Mountains where this contribution could reach 45% of the total radioactive annual dose.

Full range determination of 222Rn at the watershed scale by liquid scintillation counting

(222)Rn has been increasingly used to identify groundwater contribution to surface water. Particular attention has been paid to analytical protocols and counting parameters used for liquid alpha scintillation measurements over a range of activities covering river and groundwater domains. Direct measurements and Rn-extraction protocols are optimized, and scintillometer efficiency is calibrated using international standards over the 0.5-35 Bq/L range. The interval of activities was performed in surface water and groundwater from a small Canadian watershed.

Simulation of Distributed Base Flow Contributions to Streamflow Using a Hillslope-Based Catchment Model Coupled to a Regional-Scale Groundwater Model

AbstractThe estimation of base flow contributions to streamflow, and in particular the partitioning of base flow into shallow and deep aquifer components, is a challenging problem in water resources management. In this study, the newly developed hillslope-storage Boussinesq analytic-element (hsB/AE) model is applied to the Allen River headwater catchment in southern Quebec, Canada, which is a major recharge area for the transboundary Chateauguay River watershed. This first application to a real catchment serves as an illustration of the basic principles behind the coupled model and an investigation into interaquifer interactions and the origins of base flow contributions to streamflow. The hsB component of the model represents local groundwater flow in sloping, shallow Quaternary sediments, whereas the AE component represents the underlying regional bedrock aquifer. The catchment is partitioned into hillslopes of regular convergent, divergent, or combined planform shapes. Simulations are run for the summe...

Dynamics of groundwater floodwaves and groundwater flood events in an alluvial aquifer

In gravelly floodplains, streamflood events induce groundwater floodwaves that propagate through the alluvial aquifer. Understanding groundwater floodwave dynamics can contribute to groundwater flood risk management. This study documents groundwater floodwaves on a flood event basis to fully assess environmental factors that control their propagation velocity, their amplitude and their extension in the floodplain, and examines the expression of groundwater flooding in the Matane River floodplain (Québec, Canada). An array of 15 piezometers equipped with automated level sensors and a river stage gauge monitoring at 15-minute intervals from September 2011 to September 2014 were installed within a 0.04-km2 area of the floodplain. Cross-correlation analyses were performed between piezometric and river-level time series for 54 flood events. The results reveal that groundwater floodwave propagation occurs at all flood magnitudes. The smaller floods produced a clear groundwater floodwave through the floodplain, while the largest floods affected local groundwater flow orientation by generating an inversion of the hydraulic gradient. Propagation velocities ranging from 8 to 13 m/h, which are two to three orders of magnitude higher than groundwater velocity, were documented while the induced pulse propagated across the floodplain to more than 230 m from the channel. Propagation velocity and amplitude attenuation of the groundwater floodwaves depend both on flood event characteristics and the aquifer characteristics. Groundwater flooding events are documented at discharge below bankfull (< 0.5 Qbf). This study highlights the role of flood event hydrographs and environmental variables on groundwater floodwave properties and the complex relationship between flood event discharge and groundwater flooding. The role that groundwater floodwaves play in flood mapping and the ability of analytical solutions to reproduce them are also discussed.