Vincent Cloutier

A Review and Evaluation of the Impacts of Climate Change on Geogenic Arsenic in Groundwater from Fractured Bedrock Aquifers

Climate change is expected to affect the groundwater quality by altering recharge, water table elevation, groundwater flow, and land use. In fractured bedrock aquifers, the quality of groundwater is a sensitive issue, particularly in areas affected by geogenic arsenic contamination. Understanding how climate change will affect the geochemistry of naturally occurring arsenic in groundwater is crucial to ensure sustainable use of this resource, particularly as a source of drinking water. This paper presents a review of the potential impacts of climate change on arsenic concentration in bedrock aquifers and identifies issues that remain unresolved. During intense and prolonged low flow, the decline in the water table is expected to increase the oxidation of arsenic-bearing sulfides in the unsaturated zone. In addition, reduced groundwater flow may increase the occurrence of geochemically evolved arsenic-rich groundwater and enhance arsenic mobilization by reductive dissolution and alkali desorption. In contrast, the occurrence of extreme recharge events is expected to further decrease arsenic concentrations because of the greater dilution by oxygenated, low-pH water. In some cases, arsenic mobilization could be indirectly induced by climate change through changes in land use, particularly those causing increased groundwater withdrawals and pollution. The overall impact of climate change on dissolved arsenic will vary greatly according to the bedrock aquifer properties that influence the sensitivity of the groundwater system to climate change. To date, the scarcity of data related to the temporal variability of arsenic in fractured bedrock groundwater is a major obstacle in evaluating the future evolution of the resource quality.

Mobility and speciation of geogenic arsenic in bedrock groundwater from the Canadian Shield in western Quebec, Canada

High arsenic concentrations occur in groundwater collected from a fractured crystalline bedrock aquifer in western Quebec (Canada). Sampling and analysis of water from 59 private wells reveal that more than half of the bedrock wells exceed the Canadian guideline value of 10μg/l for arsenic, whereas shallow wells in unconsolidated surficial deposits are not affected by the contamination. The weathering of arsenic-bearing sulfides present along the mineralized fault zone is considered to be the primary source of arsenic in groundwater. High-arsenic wells are generally characterized by mildly reducing conditions (Eh<250mV), weak alkaline conditions (pH>7.4), low Ca/Na ratios, elevated dissolved Fe and Mn concentrations and high proportions of As(III). Private bedrock wells are open boreholes that likely receive groundwater from multiple contributing fractures. Hence, it is proposed that dissolved arsenic is mainly derived from the contribution to the well discharge of reducing and alkaline geochemically evolved groundwater that contains arsenic as As(III). Geochemically evolved groundwater provides favorable conditions to release arsenic by reductive dissolution of iron and manganese oxyhydroxides and alkaline desorption from mineral surfaces. Thus, high-arsenic wells would contain a high proportion of geochemically evolved groundwater, while oxidizing low-pH recharge water causes dilution and sequestration of arsenic. In relation with the chemical evolution of groundwater along the flow path, most contaminated wells are located in confined areas whereas most of the wells located in unconfined recharge areas are not contaminated. The occurrence of boreholes with high dissolved arsenic as As(V) and oxidizing conditions is attributed to extensive sulfide oxidation and alkaline desorption. This work shows that the determination of arsenic speciation provides a valuable tool to investigate the behavior of arsenic in bedrock groundwater.

Groundwater Quality, Geochemical Processes and Groundwater Evolution in the Chateauguay River Watershed, Quebec, Canada

A hydrogeochemical study was carried out in the Quebec portion of the Chateauguay River watershed. The objective was to characterize the chemical composition of groundwater in order to evaluate its quality and assess geochemical variations related to the geological and hydrogeological settings. Bulk groundwater samples were collected from 144 wells distributed evenly over the study area. Nine of the wells were sampled with a multi-level packer system, for a total of 22 multi-level samples. Samples were analysed for a comprehensive set of chemical inorganic parameters: dissolved major, minor and trace constituents, bacterial content, and stable (δ2H, δ13C, δ18O) and radioactive (3H, 14C) isotopes. Major dissolved constituents Ca, Mg, Na, K, Cl, SO4 and HCO3– ions represent more than 92% of total dissolved solids and their concentrations seem controlled by both hydrogeological and geological factors. Most water quality problems are related to aesthetic standards for potable water use (hardness, total dissolved solids (TDS), Fe and Mn concentrations) and TDS and Cl for irrigation use. Analyses of tritium (3H) and 14C confirm the inferred recharge zones and indicate the presence of variable water ages. Groundwater shows a wide range of compositions as indicated by 12 water types defined on the basis of major ions with a weak variation of chemical composition with depth. The predominant water type, Ca-HCO3, occurs in most geological and hydrogeological settings. Principal Component Analysis (PCA) and geochemical graphs were used to identify the major processes that exert a control over the chemical composition of the groundwater. Approximately 80% of the geochemical variation can be explained by mixing between fresh recharge water with more saline water associated with the former Champlain Sea, which invaded the aquifer. Secondary processes are related to ion exchange and the potential dissolution of minerals. A cross-section along a major flow path shows that the geochemical evolution of groundwater leads to relations between water type groups, geochemical processes, and groundwater flow conditions. The hydrogeochemical conceptual model infers that carbonate dissolution during recharge leads to one end-member, a Ca-HCO3 type water, which further evolves along its flow path due through ion exchange and mixing with remnant Champlain Sea water (Na-Cl), the other end-member.

Using water stable isotopes for tracing surface and groundwater flow systems in the Barlow-Ojibway Clay Belt, Quebec, Canada

This study aims to improve the understanding of surface and groundwater flow systems based on water stable isotope data in a 19,549 km2 region of the Barlow-Ojibway Clay Belt, in western Quebec, Canada. The available geochemical database contains 645 samples including precipitation, snow cores, surface waters, groundwater and springs. All samples were analyzed for water stable isotopes (δ2H-δ18O) and complementary tritium analyses were conducted on 98 groundwater and spring samples. Precipitations depict a clear temperature-dependent seasonal pattern and define a local meteoric water line (LMWL) without a latitudinal trend in δ2H-δ18O. Samples collected from the snowpack plot on the LMWL, suggesting that the bulk snowpack preserves the isotopic composition of precipitation throughout the frozen period, prior to the spring snowmelt. Surface water samples define a local evaporation line (LEL), and evaporation over inflow (E/I) ratios range between 0 and 36%. Groundwater and spring samples are evenly distributed around the LMWL, suggesting that evaporation processes are limited prior to infiltration and that surface waters do not significantly contribute to groundwater recharge. Shallow unconfined aquifers present a greater variability in δ2H-δ18O compared to confined aquifers located farther down gradient, suggesting the mixing of varied recharge waters along the regional groundwater flow system. A three-component mixing model based on isotopic and specific electrical conductivity data allows the quantification of such mixing processes. The interpretation of isotopic data constrains a regional-scale conceptual model of groundwater flow systems and describes processes related to the timing of recharge, evaporation, mixing and discharge.

A graphical approach for documenting peatland hydrodiversity and orienting land management strategies

This study focuses on the development of an approach to document the hydrological characteristics of peatlands and understand their potential influence on runoff processes and groundwater flow dynamics. Spatial calculations were performed using geographic information systems data in order to evaluate the distribution of peatlands according to (a) neighbouring hydrogeological units and (b) their position within the hydrographic network. The data obtained from these calculations were plotted in a multiple trilinear diagram (two ternary plots projected into a diamond‐shaped diagram) that illustrates the position of a given peatland within the hydrogeological environment. The data allow for the segregation of peatlands according to groups sharing similarities as well as the identification of peatlands that are most likely to have similar hydrological functions. The approach was tested in a 19,549 km2 region of the southern portion of the Barlow‐Ojibway Clay Belt (in Abitibi‐Temiscamingue, Canada) and lead to a conceptual model representing the hydrological interactions between peatlands, aquifers, and surface waters. This approach allows for a geographic information systems‐based differentiation of headwater peatland complexes that are likely to interact with aquifers and to supply continuous baseflow to small streams from lowland peatland complexes of the clay plain that are isolated from surrounding aquifers but that can act as storage reservoirs within the hydrographic network. The typology is further used to discuss land management strategies aimed at preserving peatland hydrodiversity within the study region. The proposed approach relies on widely applicable hydrogeological and hydrographic criteria and provides a tool that could be used for assessing peatland hydrodiversity in other regions of the planet.

A review of GIS-integrated statistical techniques for groundwater quality evaluation and protection

Water quality evaluation is critically important for the protection and sustainable management of groundwater resources, which are variably vulnerable to ever-increasing human-induced physical and chemical pressures (e.g., overexploitation and pollution of aquifers) and to climate change/variability. Preceding studies have applied a variety of tools and techniques, ranging from conventional to modern, for characterization of the groundwater quality worldwide. Recently, geographic information system (GIS) technology has been successfully integrated with the advanced statistical/geostatistical methods, providing improved interpretation capabilities for the assessment of the water quality over different spatial scales. This review intends to examine the current standing of the GIS-integrated statistical/geostatistical methods applied in hydrogeochemical studies. In this paper, we focus on applications of the time series modeling, multivariate statistical/geostatistical analyses, and artificial intelligence techniques used for groundwater quality evaluation and aquifer vulnerability assessment. In addition, we provide an overview of salient groundwater quality indices developed over the years and employed for the assessment of groundwater quality across the globe. Then, limitations and research gaps of the past studies are outlined and perspectives of the future research needs are discussed. It is revealed that comprehensive applications of the GIS-integrated advanced statistical methods are generally rare in groundwater quality evaluations. One of the major challenges in future research will be implementing procedures of statistical methods in GIS software to enhance analysis capabilities for both spatial and temporal data (multiple sites/stations and time frames) in a simultaneous manner.

Stratigraphic sequence map for groundwater assessment and protection of unconsolidated aquifers: A case example in the Abitibi-Témiscamingue region, Québec, Canada

Quantifying water resources at various scales is crucial for ensuring safe access to potable water for present and future generations. Fitting into this framework, this study presents a GIS-based approach aimed at allowing the evaluation of available groundwater resources within unconsolidated aquifers set in vast and heterogeneous shield regions. The approach was developed in a 19,397 km2 region located in Abitibi-Témiscamingue (Québec, Canada) where unconsolidated aquifers are set in an irregular geologic framework owing to the rugged Canadian Shield topography and to the diversity of glacial and post-glacial events that shaped the landscape. Sparse and unevenly distributed stratigraphic borehole data were used for constructing a GIS-based model representing the total overburden thickness (unconsolidated sediments covering the bedrock), as well as the thickness and extent of the fine-grained deep-water sediments (silt-clay) of the Barlow-Ojibway proglacial Lake on a 100 m × 100 m mesh at the regional scale. The study suggests a new set of functions to map drift thickness in areas of highly discontinuous drift cover, using GIS tools. These data were used jointly with surficial deposits maps in order to develop a regional-scale two-dimensional model where 15 distinct stratigraphic sequences allow representation of the architecture of unconsolidated geological units. The approach allows simplifying data representation and ensuring the consistency between various regional-scale hydrogeological maps. The data is used to produce a stratigraphic sequence map for representing the architecture of aquifer– aquitard systems at the regional scale in a manner that is intelligible to non-specialists. The map and related data are first discussed for documenting the extent and volume of regional aquifers. The stratigraphic sequence map is subsequently discussed as a tool for supporting political decision makers dealing with issues related to groundwater resource assessment, protection and sustainable development.

Conceptual model of regional groundwater flow based on hydrogeochemistry (Montérégie Est, Québec, Canada)

The groundwater geochemistry of the fractured rock aquifer system in the Montérégie Est region, southern Quebec, Canada, was studied as part of a regional groundwater resources assessment. The 9218 km² study area included three major watersheds that were divided into five hydrogeological contexts: Northern St. Lawrence Lowlands, Southern St. Lawrence Lowlands, Appalachian Uplands, Appalachian Piedmont and Monteregian Hills. A large part of this study area was invaded by the Champlain Sea from 13,000 to 11,000 years ago. Study objectives were to identify the mechanisms controlling groundwater composition and to support the understanding of the aquifer hydrodynamics. Groundwater from 206 wells drilled into the rock aquifer was sampled and analyzed for conventional parameters and isotopic analyses were also done on selected samples (δ2H, δ18O and 3H of water; δ13C and 14C of dissolved inorganic carbon). The interpretation of geochemical results was based on a multivariate statistical analysis, which led to the definition of eight water groups. The study allowed the delineation of a 2200-km² zone containing brackish groundwater of marine origin in the northwestern part of the study area. This zone is surrounded by sodic and alkaline groundwater originating from Na+-Ca2+ ionic exchange. Young groundwater and therefore recharge zones were only encountered in the southern part of the Lowlands, in the northern part of the Piedmont and in the Appalachian Uplands. In the southern part of Lowlands, recharge is presumed to be slow and water composition shows the influence of the former presence of the Champlain Sea. Relatively deep groundwater circulation was also inferred to occur from the Appalachian Uplands toward mixing zones mainly located to the west at the Appalachian frontal thrust faults and around the Monteregian Hills. The geochemical interpretation provided indications on regional recharge and discharge zones as well as groundwater flow, which could not have been determined otherwise.

The role of hydrogeological setting in two Canadian peatlands investigated through 2D steady-state groundwater flow modelling

ABSTRACT This study investigated how hydrogeological setting influences aquifer–peatland connections in slope and basin peatlands. Steady-state groundwater flow was simulated using Modflow on 2D transects for an esker slope peatland and for a basin peatland in southern Quebec (Canada). Simulations investigated how hydraulic heads and groundwater flow exported toward runoff from the peatland can be influenced by recharge, hydraulic properties, and heterogeneity. The slope peatland model was strongly dominated by horizontal flow from the esker. This suggests that slope peatlands are dependent on the hydrogeological conditions of the adjacent aquifer reservoir, but are resilient to hydrological changes. The basin peatland produced groundwater outflow to the surface aquifer. Lateral and vertical peat heterogeneity due to peat decomposition or compaction were identified as having a significant influence on fluxes. These results suggest that basin peatlands are more dependent on recharge conditions, and could be more susceptible to land use and climate changes.