Current and Recent Research Projects

Risk to Water Quality and Quantity

Mapping the Likelihood of Flowing Artesian Conditions in the Okanagan Basin and Lower Fraser Valley, British Columbia (MSc, Johnson, 2021) Drilling into an aquifer that exists under flowing artesian conditions can have significant environmental and socio-economic consequences. Unfortunately, the risk of drilling a flowing well in different regions of British Columbia is largely unknown. This study used Multi-Criteria Analysis and the Analytical Hierarchy Process to map the likelihood of flowing well occurrence in three hydrogeologic contexts based on aquifer type and location; 1) Okanagan Bedrock, 2) Okanagan Unconsolidated and 3) Fraser Valley Unconsolidated. A novel approach was used to determine the weights of each factor by using the Wilcoxon Rank Sum Test. Each map was determined to reasonably display and predict the extent of flowing artesian conditions by using the Wilcoxon Rank Sum Test to a significance level of 15%. Through site specific evaluations, Well Drilling Advisories are recommended to be developed by the Province for Naramata, Armstrong, and Aquifer 33 in Langley to alert drillers of flowing artesian conditions.

Shale Gas Development Impacts on Water Resources. For the past several years we have been conducting hydrogeological and hydrogeochemical studies in Northeast British Columbia to assess various impacts of shale gas development on water resources. Shale Gas Development in NE BC was used as a case study area for evaluating risk and resilience in the context of the Water-Energy-Food Nexus. This part of the research was funded in part by the Research Institute for Humanity and Nature (Kyoto, Japan) as part of a larger circum-Pacific project on nexus issues, and the Pacific Institute for Climate Solutions (PICS). The research has been carried out in partnership with the BC Ministry of Forests, Lands and Natural Resources, BC Ministry of Environment, Geoscience BC, the BC Oil and Gas Commission, the School for Public Policy at SFU (among others).

Establishing Regional Groundwater Chemistry Baselines for Northeast BC to Support Interpretation of Groundwater Analytical Results (MSc, Martinolich, in progress) This research aims to establish groundwater chemistry baselines for the Peace River Region (PRR), by carrying out a paired hydrogeochemical and multivariate statistical analysis using a regional groundwater chemistry dataset. Establishing groundwater chemistry baselines for the PRR will support ongoing hydrogeological research, provide knowledge for groundwater monitoring, remediation and reclamation programs, and aid recognition of incipient groundwater pollution.

Geochemical Compositions of Subsurface Materials and the Potential Impacts on Groundwater Quality from Fugitive Methane within the Peace Region, British Columbia (MSc, Allen, A., 2021) Methane accumulation in groundwater from gas migration induced by hydrocarbon extraction can lead to degradation of water quality. This research focuses on the potential impact elevated concentrations of methane can have on shallow groundwater systems based on the geochemical properties of the sediments and rocks that are present. Bulk chemical characterization of major oxides, trace elements, and rare earth elements indicate that lithologically similar geologic materials in the Peace Region of British Columbia share some geochemical similarities at a regional scale. Investigation of hydrogeochemical characteristics through sequential extraction experiments suggests there is little difference between materials regarding their potential for mobilization of trace elements into groundwater. PHREEQCi numerical models of methane interacting with these materials indicate there is a greater impact in reducing compared to oxidizing conditions, and the presence of carbonate and iron oxyhydroxide minerals mediate the influence large fluxes of methane have on the hydrochemistry.

Characterizing the Contamination Risk in the Surface/Shallow Subsurface Zone due to Wastewater Spills and Leaks from Shale Gas Activities (MSc, Rosales-Ramirez, 2020). The study area is the Peace River Region in northeastern British Columbia where shale gas activities are rapidly expanding. Numerical models (TOUGH2) are being developed to define the extent of the contamination footprint, as well as to determine the rate of transport and fate of the wastewater in case of accidental release.  The results will be integrated into the shallow groundwater intrinsic vulnerability map (above) to evaluate and quantify the risk from these events in an effort to define areas with low-high risk response to wastewater spills or leakages. The resulting maps can be used in decision making to reduce the risk of potential water contamination, as well to raise awareness and promote better monitoring practices.

Investigating the Role of Buried Valley Aquifer Systems in the Regional Hydrogeology of the Peace River Region of Northeast British Columbia (MSc, Morgan, 2018). Geological and numerical flow models were developed to explore the hydraulic role of buried valley aquifers in regional groundwater flow and assess the potential groundwater resource. The study area was the central Peace Region in Northeast British Columbia. The reservoir software Petrel was used to construct the geological model of a buried valley network by integrating interpretations from an airborne electromagnetic survey (SkyTEM) and borehole gamma-ray and lithology logs. This detailed geological model and a simplified geological model were used to develop two numerical flow models in MODFLOW. The modelling results suggest that permeable deposits exist within the buried valleys, but are not regionally connected throughout the whole network, and thus do not play a significant role in the regional groundwater flow regime. However, extensive permeable deposits within the buried valleys appear to exist at smaller scales, and may offer a viable water source in the area.

Numerical Modeling of Deep-Aquifer Wastewater Disposal in Northeast British Columbia (MSc, Simons, 2018). Dense, saline wastewater generated during oil and gas activities (hydraulic fracturing and production) is commonly disposed of in deep formations, but the migration of this wastewater after entering the subsurface is poorly understood. This study uses numerical models to simulate wastewater disposal in the Paddy-Cadotte of Northeast British Columbia using both single-well axisymmetric box models and a regional model of the formation in which multiple disposal and water source wells operate. A sensitivity analysis performed on the box models reveals that dispersivity and permeability exert the strongest control on overall wastewater distribution. Models show that wastewater migrates further than predicted using a simple volumetric calculation, and extends further along the base of the formation than the top due to variations in fluid density. Interference between disposal and source wells is observed to influence wastewater migration, while formation dip and regional groundwater flow have no discernible impact.

Springs in Northeast British Columbia (Bystron, MSc, 2018) Springs in the Peace River Regional District are important fresh water sources for domestic and industrial use, however, little is known about their hydrogeology. This study investigated the controls on spring occurrence and delineated the spring source areas (SSA). Samples from 36 springs were collected and analyzed for chemical and isotopic composition. The chemistry suggested 20 Quaternary and 16 Bedrock aquifer sourced springs. A GIS-based multi-criteria decision analysis integrated spring-related factors (slope, curvature, hydrogeological features, topographic wetness index, surficial geology, and drift thickness) to produce maps showing the likelihood for occurrence of Quaternary or Bedrock-sourced springs. There was a ~70% success rate for correctly identifying sources of the sampled springs. A GIS-based tool was applied to determine topographic SSA’s for 33 springs with varying degrees of success. Although both the maps and SSAs contain uncertainty, there is sufficient information to support decision-making in water resource protection and management.

Investigating Regional Groundwater Flow Influences on Slope Stability in Unlithified Materials (MSc, Dandurand, 2018). Elevated pore pressures play an integral part in slope failure; however, the nature of the regional groundwater flow regime is often not incorporated in slope stability analysis. This research aims to improve our understanding of how the dynamic nature of the groundwater flow regime exerts control on the regional occurrence of landslides in unlithified materials. An integrated hydrogeological – geotechnical methodology is used to investigate the hydrogeological controls on slope stability in Northeast British Columbia. A two-part approach utilized both a series of steady-state groundwater models developed using SVFLUX to investigate the role of geologic contact geometry and hydrogeological characteristics, and a transient groundwater model developed using MIKE SHE to investigate the role of climate. The analysis of these groundwater results highlights implications of regional and local scale hydrogeological processes on pore water pressures. Several challenges were encountered in regard to investigating processes at fine scale resolution within in a large scale groundwater flow model, and an inability to export unsaturated zone results from MIKE SHE limited the slope stability modeling in SV SLOPE.

Shallow Groundwater Intrinsic Vulnerability Mapping for North East British Columbia (Post-Doctoral Fellow, Holding, 2015) Simon Fraser University (SFU), with financial support from the BC Ministry of Forests, Lands and Natural Resource Operations (FLNRO) and the Pacific Institute for Climate Solutions (PICS), developed a Shallow Groundwater Intrinsic Vulnerability Map of Northeast British Columbia (BC).  The assessment was conducted in response to mounting concerns surrounding water management and protection in Northeast BC in relation to shale gas development. The intent of the mapping is to characterise the intrinsic vulnerability of near surface geological materials to contamination originating at land surface. The resulting map is intended to support agencies in the development of policies and regulations that protect groundwater quality.


Aquifer-Stream Hydraulic Connectivity

Groundwater - Surface Water Exchange Dynamics in Low- Gradient and Tidally Influenced Streams in the Lower Fraser Valley (PhD, Mitton, in progress) The study aims to build an understanding of groundwater-surface water exchanges in tidally influenced and low-gradient stream systems, and how hydrogeological and hydrological processes interact to drive hydrological variability. A main goal is to link aquifer and stream dynamics mechanistically to environmental response metrics (e.g., hydraulic-habitat, geochemical, and biological variables) to support management decisions related to the evaluation of potential impacts of hydrological alteration (e.g., groundwater pumping) and the determination of environmental flow needs (EFNs) in groundwater dependent, low- gradient, and tidally influenced streams.

An Ecological Assessment of Aquifer-Stream Connectivity in a Groundwater Dependent Stream (MSc, Mitton, 2021) The ecological roles of groundwater are rarely properly addressed in water management due to a lack of understanding of the mechanisms influencing river ecosystems and inadequate methods to assess potential impacts of activities like groundwater pumping. Here, a hydroecological field study was undertaken in Bertrand Creek, South Coast region of British Columbia, to characterize aquifer-stream exchanges, their drivers, how they influence aquatic habitat and macroinvertebrates, and to develop a method to assess potential groundwater pumping impacts. Data collected from May to October from 2018 to 2020 showed declining groundwater levels in response to hotter and drier conditions, but groundwater discharge was sufficient for streamflow to persist downstream. In 2020, groundwater cooled the stream and maintained habitat conditions that allowed macroinvertebrates to survive the summer. Habitat models built with these data were integrated with an analytical streamflow depletion model to produce a tool for assessing the ecological impacts of groundwater pumping.

Quantifying Aquifer-Stream Exchanges Along Bertrand Creek, British Columbia and Washington State, Using a Numerical Groundwater Flow Model (BSc USRA project, Nott, 2020) The purpose of this study was to quantify aquifer-stream exchanges along Bertrand Creek, in order to identify reaches where streamflow drought may occur. A steady-state three-dimensional numerical groundwater flow model of the watershed was developed. See Research Highlights for a link to the report.

Aquifer-Stream Connectivity at Otter Park, Langley, BC (BSc Honours, Johnson, 2018) Streams that are hydraulically connected to an aquifer can be adversely impacted by groundwater abstraction, which may result in streamflow depletion.  Characterisation of the potential for streamflow depletion is thus an important aspect of water management. The purpose of this study is to characterise the aquifer and streambed sediment properties at Otter Park in Langley, British Columbia, in order to evaluate streamflow depletion due to pumping. The aquifer and streambed sediment properties were characterised through a series of methods, including grain size analysis, slug and bail tests, and a 30 hour pumping test. Streambed sediments have a hydraulic conductivity approximately one order of magnitude lower than the aquifer sediments. Hydraulic connectivity between the aquifer and the stream was assessed using stream discharge, vertical hydraulic gradients in nested instream piezometers, and streambed temperature. Several lines of evidence indicate that streamflow depletion occurred; however, the variation in the streamflow data limits the conclusions that can made without a longer pumping test. 

Investigating the Thermo-Hydraulic Regime of Union Creek, Langley BC (BSc Honours, Whistler, 2018) Groundwater abstraction near a stream can impact the surface water by changing the direction and magnitude of the groundwater flux and the thermal regime of the stream. In this study, a thermo-hydraulic survey was conducted at Steele Park, Langley B.C, to assess the thermal and hydraulic regime of Union Creek through the summer season and under non-pumping and pumping conditions. Groundwater flux measurements across the streambed and cross-correlation analyses between streambed temperature and air temperature, show Union Creek to be primarily a losing stream. The temperature signals recorded at the streambed at upstream locations lagged behind  those downstream due to a more immediate response to radiative heating in shallower waters downstream. From non-pumping to pumping, streambed temperatures and stream stage did not change, and unfortunately, the vertical flux results were uncertain. Thermal infrared images did not show groundwater discharge into the stream, but suggest seepage of warm water from the bank storage during summer. Analysis and interpretation show the significance of properly characterizing hydraulic properties to investigate the hydraulic connectivity with confidence.

Investigation of Aquifer-Stream Connectivity at Steele Park, Langley BC (BSc Honours, Hall, 2017) Abstraction of groundwater from a pumping well located beside a stream can result in sourcing of the pumped water directly from the stream and consequent depletion of stream discharge. The ability to accurately estimate stream depletion due to pumping is necessary for water rights management. A 48-hour pumping test was conducted at Steele Park, Langley BC to assess aquifer-stream connectivity and the effects of pumping on streamflow and groundwater seepage. The site was also modeled using a suite of available analytical solutions for streamflow depletion. Analysis of the field data highlighted the importance of accurately characterizing the hydraulic properties of the streambed sediments and properly designing and implementing the field test. Techniques for measuring field parameters crucial for direct comparison between field data and the analytical streamflow depletion models are discussed.


Climate Change Impacts on Water

A Morphological Approach to Evaluating Risk to Freshwater Resources in Coastal Deltas (PhD, Anderson, in progress). Today, climate change and rising sea levels threaten to sink coastal deltas, potentially leaving millions of people to find new homes and causing vast ecosystem loss. Population growth and increased food production have led to increased reliance on groundwater resources to meet freshwater demands. However, groundwater overdrafts have resulted in freshwater sources becoming contaminated with salinity, making them unsuitable for use. This research tackles the freshwater availability challenges facing coastal deltas by using geologic history to understand the processes controlling freshwater contamination in deltas. Population and socioeconomic factors are used to determine where freshwater contamination will most impact coastal populations. This research will bring attention to the challenges in freshwater protection and will help determine where science and engineering efforts can be most beneficial to society. 

Reconstructing Groundwater Levels from Tree Ring Widths in the North American Cordillera (PhD, Hunter, 2021) Water resources around the world are periodically impacted by droughts, which can reduce water available for human use and put stress on the environment. The effects of long-term changes in climate on groundwater resources are difficult to determine, as most observed groundwater level records are short (< 60 years). However, observational records can be extended using proxy records such as tree ring widths. This research uses two methods to study past groundwater changes (1850-2002) in the North American Cordillera: a nested principal component regression technique to reconstruct past groundwater levels from tree ring widths, and a nonparametric approach which reconstructs daily temperature and precipitation based on annual tree ring widths, to be used as input into a hydrological model for simulating groundwater levels. The past spatial and temporal occurrence of droughts in this region is examined, as are the associations between past groundwater levels and climate teleconnections. The groundwater level reconstructions demonstrate that the relationships between groundwater levels and snowpack, and tree ring widths and snowpack, are important for creating groundwater level reconstructions in mountain regions. In addition, separating groundwater observation wells by dominant aquifer-stream system strengthens the reconstructions. Analysis of past groundwater drought (between 1850 and 2004) shows that drought events were not synchronous across the North American Cordillera, with variations in timing occurring between the north and south cordillera, as well as between different aquifer-stream systems. Teleconnections patterns that affect observed groundwater levels in this region were also found to influence past groundwater levels. Hydrologic modeling results show that the nonparametric approach was useful for creating simulations of past groundwater levels, and that in the past time periods analyzed (1920-1950 and 1950-1980), peak recharge occurred earlier in the year than in the observed climate; this earlier peak in recharge is likely due to warmer spring temperatures and earlier snowmelt during the paleo-periods. These results show that the relationship between groundwater levels and tree rings can be used to successfully reconstruct and simulate past groundwater levels in mountain regions, providing a long-term perspective into groundwater variability throughout the past. 

Snow Drought and Streamflow Drought in Western North America in the Context of a Warming Climate (PhD, Dierauer, 2018) In western North America, snowpack supplies much of the water used for irrigation and for municipal and industrial uses, and snowmelt recharges groundwater and provides ecosystem-sustaining baseflow during low flow periods. Continued climate warming is expected to have large impacts on snowmelt hydrology, with subsequent impacts on low flows and snow and streamflow drought regimes. This research combined two separate methodologies, a data-driven (downward) approach and a process-based (upward) approach, to improve our understanding of snow drought and streamflow drought in the context of a warming climate. The data-driven approach combined observed hydroclimatic time series with multiple statistical methods, including bivariate and partial correlation and temporal and spatial analogs. The process-based approach combined climate change projections and hydrological modelling. The two approaches yielded consistent results that, together, illustrate that snow drought, low flows, and streamflow drought are sensitive to winter climate conditions, particularly precipitation and thawing degrees. In the context of climate warming, increased winter season thawing degrees leads to increased warm (temperature-driven) snow drought, shorter and less severe winter low flows, longer and more severe summer low flows, and increased summer streamflow drought risk. Further, both approaches showed that the response of snowmelt hydrology to climate warming is non-linear, and regions with winter temperatures near 0°C exhibit substantially larger impacts from +2°C of warming compared to regions with winter temperatures far below 0°C. Temperature-driven shifts in snow drought, low flows, and streamflow drought regimes will have widespread implications for surface water supply security. Increased frequency of warm snow droughts will likely lead to an increased frequency of mid-winter melt events, which will create challenges for water management. As summer low flow periods become more severe and snow-drought related summer streamflow droughts become more frequent, the potential for more severe summer water shortages increases. The most severe shortages will likely occur due to the co-occurrence of warm and dry conditions.


Bonding environment and distribution of arsenic and selenium in pyrite and controls on element mobilization. (PhD, J. Ross-Lindeman, in progress) is investigating the reaction mechanisms of arsenic and selenium oxidation and mobilization during arsenic- and selenium-rich pyrite oxidation, a process that can potentially introduce contaminants to surface water and groundwater. This research combines chemical and mineralogical characterization of arsenic and selenium in pyrites, including trace element concentration, speciation, and distribution, with experiments on pyrites to understand the effects of arsenic or selenium substitution on oxidation reaction mechanisms.

Investigating the Geochemistry of Selenium in the Residual from Biologically Treated Mine-Impacted Waters (MSc, L. Devolder, 2018) The goal of this study is to determine which selenium species are present in the residual and under what conditions it may become mobile in aqueous phase. Detailed selenium speciation of the residual is being determined using sequential extraction procedures in the laboratory and X-ray Absorption Near Edge Structure analysis at two separate synchrotron facilities.  A series of aqueous geochemical batch experiments that simulate different redox and pH conditions are being used to determine under what conditions the selenium may mobilize and if certain species are more vulnerable than others. This study aims to help quantify the risk posed by the residual and eventually develop safe disposal methods. 

Sources and Distribution of Arsenic in Groundwater on the Gulf Islands of British Columbia (BSc Honours, A. Allen, 2018). Arsenic mobilization to groundwater results in significant impacts on groundwater quality which can lead to adverse health effects. Understanding in what form the arsenic occurs is critical to determining the potential for arsenic contamination. Mineralogical and bulk chemical analyses, a cation exchange capacity experiment, and a chemical sequential extraction procedure for arsenic speciation were performed to determine the abundance and occurrence of arsenic in rocks from several formations of the Nanaimo Group in the Gulf Islands of British Columbia. Through these methods it was found that arsenic is strongly associated with iron oxides, oxyhydroxides, and sulfides. Results indicate finer grained formations have higher arsenic content, and a minor correlation exists between formations which contain elevated arsenic in groundwater and higher abundances of more easily mobilized arsenic.