Characterization and distribution of metal and nonmetal elements in the Alberta oil sands region of Canada

Atmospheric deposition mapping of particulate elements in the Canadian Athabasca oil sands region

This study used a deposition modeling framework to generate gridded dry, wet, and total (dry + wet) deposition fluxes of 27 particulate elements over the Canadian Athabasca oil sands region and its surrounding areas for the years 2016-2017. The framework employed the element concentrations from the CALPUFF dispersion model outputs that were bias-corrected against measured concentrations, modeled dry deposition velocities, precipitation analysis data, and literature values of element-specific fine mode fractions and scavenging ratios by rain and snow. The annual total deposition (mg/m²/year) of all elements (EM) across the domain ranged from 4.49 to 5450 and the mean and median deposition were 60.9 and 31.0, respectively. Total EM deposition decreased rapidly within a short distance from the oil sands mining area. Annual mean total deposition (mg/m²/year) of EM was 717 in Zone 1 (within 30 km from a reference point, representing the center of the oil sands mining area), 115 in Zone 2 (30-100 km from the reference point), and 35.4 in Zone 3 (beyond 100 km from the reference point). The deposition of individual elements was primarily governed by their respective concentrations and among all elements the annual mean total deposition (μg/m²/year) over the domain varied five orders of magnitude ranging from 0.758 (Ag) to 20,000 (Si). Annual mean dry and wet deposition (mg/m²/year) of EM over the domain were 15.7 and 45.2, respectively. Aside from S, which has relatively lower precipitation scavenging efficiencies, wet deposition was the dominant deposition type in the region contributing from 51% (Pb) to 86% (Ca) of the respective total deposition. Total EM deposition over the domain in the warm season (66.2 mg/m²/year) was slightly higher than that in the cold season (55.6 mg/m²/year). Deposition of individual elements in Zone 1 were generally lower than their deposition at other sites across North America.

A numerical modeling framework for simulating the key in-stream fate processes of PAH decay in Muskeg River Watershed, Alberta, Canada

Most previous water quality studies oversimplified in-stream processes for modeling the fate and transport of critical organic contaminants, such as Polycyclic Aromatic Hydrocarbons (PAHs). Taking four selected PAHs as representative organic contaminants, we developed a numerical modeling framework using a Water Quality Analysis Simulation Program 8 (WASP8) and a well-established watershed model, i.e., Soil and Water Assessment Tool (SWAT) to: (1) address the influence of in-stream processes, including direct photolysis, volatilization, partitioning of PAHs to suspended solids, and DOC complexation processes on PAH concentrations; and (2) establish relationships between spatiotemporal distribution of environmental factors (e.g., ice coverage, water temperature, wind, and light attenuation), in-stream processes, and PAH concentrations at a watershed scale. Using calibrated SWAT and WASP8 models, we evaluated the impacts of seasonal changes in environmental factors on in-stream processes in the Muskeg River watershed, which is part of the Athabasca Oil Sands Region (AOSR), the third-largest crude oil reserves of the world in western Canada. Among four selected PAHs, simulation results suggest that Naphthalene primarily decay in the water through volatilization or direct photolysis. For Phenanthrene, Pyrene, and Chrysene, DOC complexation, volatilization, and direct photolysis all contribute to their decay in the water, with a strong dependence on seasonality. Model simulations indicated that direct photolysis and volatilization rates are meager in cold seasons, mainly due to low river temperature and ice coverage. However, these processes gradually resume when entering the warm season. In summary, the model simulation results suggest that critical in-stream processes such as direct photolysis, volatilization, and partitioning and their relationship with environmental factors should be considered when simulating the fate and transport of organic contaminants in the river systems. Our results also reveal that the relationship between environmental factors and fate processes affecting PAH concentrations can vary across a watershed and in different seasons.

A decadal synthesis of atmospheric emissions, ambient air quality, and deposition in the oil sands region

This review is part of a series synthesizing peer-reviewed literature from the past decade on environmental monitoring in the oil sands region (OSR) of northeastern Alberta. It focuses on atmospheric emissions, air quality, and deposition in and downwind of the OSR. Most published monitoring and research activities were concentrated in the surface-mineable region in the Athabasca OSR. Substantial progress has been made in understanding oil sands (OS)-related emission sources using multiple approaches: airborne measurements, satellite measurements, source emission testing, deterministic modeling, and source apportionment modeling. These approaches generally yield consistent results, indicating OS-related sources are regional contributors to nearly all air pollutants. Most pollutants exhibit enhanced air concentrations within ~20 km of surface-mining activities, with some enhanced >100 km downwind. Some pollutants (e.g., sulfur dioxide, nitrogen oxides) undergo transformations as they are transported through the atmosphere. Deposition rates of OS-related substances primarily emitted as fugitive dust are enhanced within ~30 km of surface-mining activities, whereas gaseous and fine particulate emissions have a more diffuse deposition enhancement pattern extending hundreds of kilometers downwind. In general, air quality guidelines are not exceeded, although these single-pollutant thresholds are not comprehensive indicators of air quality. Odor events have occurred in communities near OS industrial activities, although it can be difficult to attribute events to specific pollutants or sources. Nitrogen, sulfur, polycyclic aromatic compounds (PACs), and base cations from OS sources occur in the environment, but explicit and deleterious responses of organisms to these pollutants are not as apparent across all study environments; details of biological monitoring are discussed further in other papers in this special series. However, modeling of critical load exceedances suggests that, at continued emission levels, ecological change may occur in future. Knowledge gaps and recommendations for future work to address these gaps are also presented. Integr Environ Assess Manag 2022;18:333-360. © 2021 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

The Gut Microbial Community Structure of the North American River Otter (Lontra canadensis) in the Alberta Oil Sands Region in Canada: Relationship with Local Environmental Variables and Metal Body Burden

The Alberta Oil Sands Region in Canada is home to one of the largest oil bitumen deposits in the world. The North American river otter (Lontra canadensis) is a top predator with a small home range and is sensitive to disturbances; it has been designated as a sentinel species for the potential impacts of the natural resource exploitation on freshwater ecosystems in the Alberta Oil Sands Region. With an increasing interest in noninvasive biomarkers, recent studies suggest that gut microbiota can be used as a potential biomarker of early biological effects on aquatic wildlife. The goal of the present study was to determine the river otter gut microbial structure related to environmental variables characterizing mining activities and metal body burden. We obtained 18 trapped animals from and surrounding the surface mineable area of the Alberta Oil Sands Region. The gut microbial community structure was characterized using high-throughput sequencing of 16S rRNA gene amplicon analyses. Trace metal concentrations in the liver were measured by inductively coupled plasma-mass spectrometry. Our study revealed that the gut bacteria of river otters in the Alberta Oil Sands Region clustered in 4 groups dominated by Peptostreptococcaceae, Carnobacteriaceae, Enterobacteriaceae, Clostridiaceae, and Nostocaceae. We show that arsenic, barium, rubidium, liver-body weight ratio, and δ¹⁵ N were associated with each cluster. When comparing affected versus less affected sites, we show that river otter gut bacterial community and structure are significantly related to trophic level of the river otter but not to Alberta Oil Sands Region mining activities. Our study reveals that the gut bacterial dynamics can provide insights into the diet and habitat use of river otters but that more work is needed to use it as a pollution biomarker. Environ Toxicol Chem 2020;39:2516-2526. © 2020 SETAC.

Introduction to the virtual special issue monitoring ecological responses to air quality and atmospheric deposition in the Athabasca Oil Sands region the wood Buffalo environmental Association's Forest health monitoring program

The expansion of oil sands resource development in the Athabasca Oil Sands Region in the early 1990's led to concerns regarding the potential ecological and health effects of increased emissions and deposition of acidic substances. Conditions attached to a 1994 approval for an oil sands facility expansion led to the creation of the Wood Buffalo Environmental Association, and its Terrestrial Environmental Effects Monitoring committee. This multi-stakeholder body was tasked with development and operation of an environmental (forest health) monitoring program for the detection of ecological responses to atmospheric emissions and deposition. Initially focused on acid deposition monitoring, jack pine forest, growing on sandy soils with limited acid buffering capacity, was selected as the receptor system. An initial set of 10 monitoring locations was established using the Canadian Acid Rain Network Early Warning System methodology (since increased to 27, with three lost to development). Ecological monitoring is on a 6-year cycle, with concurrent measures of soil, needle and lichen chemistry, and tree and understory condition, together with ongoing measurements of air quality and atmospheric deposition. Because jack pine forest edges facing the emissions sources were expected to be more exposed to acidic emissions, evaluation of stand edge monitoring locations began in 2008. Monitoring of a targeted suite of indicators began in 2012 at 25 jack pine stand edge monitoring sites. This special issue presents the results derived from biophysical sampling campaigns (1998 to 2013), coupled with ongoing ambient atmospheric, deposition and epiphytic lichen monitoring (data through 2017) and source apportionment studies, as well as papers contributed by others engaged in regional research and monitoring programs. The Forest Health Monitoring Program provides data supportive of regulatory and stakeholder evaluations of environmental quality, and is adaptive to new needs, extreme environmental events and technological development while providing continuity of monitoring.

Incorporating a non-reactive heavy metal simulation module into SWAT model and its application in the Athabasca oil sands region

Heavy metal contaminations in an aquatic environment is a serious issue since the exposure to toxic metals can cause a variety of public health problems. A watershed-scale model is a useful tool for predicting and assessing heavy metal fate and transport in both terrestrial and aquatic environments. In this study, we developed a simulation module for non-reactive heavy metals and incorporated it into the widely used Soil Water Assessment Tool (SWAT) model. The simulated processes in the developed model include heavy metal deposition, partitioning in soil and water, and transport by different pathways in both terrestrial and aquatic environments. Three-phase partitioning processes were considered in the module by simulating heavy metals portioning to dissolved organic carbon in the soil and stream. This developed module was used for watershed-scale simulation of heavy metal processes in the Muskeg River watershed (MRW) of the Athabasca oil sands region in western Canada for the first time. The daily streamflow and sediment load from 2015 to 2017 were first calibrated and validated. Subsequently, the daily Lead and Copper loads at the outlet station were used for heavy metal calibration and validation. The performances for the daily heavy metal loads simulation during the whole simulation period can be considered as "satisfactory" based on the recommended model performance criteria with the Nash-Sutcliffe efficiency as 0.41 and 0.71 for Pb and Cu loads, respectively. The simulation results indicate that the spring and summer are hot moments for heavy metal transport and the snowmelt in spring and rainfall runoff events in summer are the main driving forces for the metal transport in the MRW. We believe the developed model can be a useful tool for simulating the fate and transport of non-reative heavy metals at watershed scale and further used to assess management scenarios for mitigating heavy metal pollution in the Athabasca oil sands region.

In vivo solid-phase microextraction sampling combined with metabolomics and toxicological studies for the non-lethal monitoring of the exposome in fish tissue

Various environmental studies have employed the biomonitoring of fish in their aquatic ecosystems in order to identify potential metabolic responses to the exposome. In this study, we applied in vivo solid-phase microextraction (SPME) to perform non-lethal sampling on the muscle tissue of living fish to extract toxicants and various endogenous metabolites. Sixty white suckers (Catastomus commersonii) were sampled from sites upstream, adjacent, and downstream from the oil sands development region of the Athabasca River (Alberta, Canada) in order to track their biochemical responses to potential contaminants. In vivo SPME sampling facilitated the extraction of a wide range of endogenous metabolites, mainly related to lipid metabolism. The obtained results revealed significant changes in the levels of numerous metabolites, including eicosanoids, linoleic acids, and fat-soluble vitamins, in fish sampled in different areas of the river, thus demonstrating SPME's applicability for the direct monitoring of exposure to different environmental toxicants. In addition, several classes of toxins, including petroleum-related compounds, that can cause serious physiological impairment were tentatively identified in the extracts. In vivo SPME, combined with the analysis of contaminants and endogenous metabolites, provided important information about the exposome; as such, this approach represents a potentially powerful and non-lethal tool for identifying the mechanisms that produce altered metabolic pathways in response to the mixtures of different environmental pollutants.

Measuring the distribution of trace elements amongst dissolved colloidal species as a fingerprint for the contribution of tributaries to large boreal rivers

Organic and inorganic colloids play important roles governing the speciation, transport, and bioaccessibility of trace elements in aquatic systems. These carriers are especially important in the boreal zone, where rivers that contain high concentrations of iron and organic matter are prevalent. The distribution of trace elements amongst different colloidal species (or "speciation profile") can therefore be useful as a fingerprint to detect different trace element sources and for tracking colloid transformations, with implications for bioaccessibility. Asymmetrical flow field-flow fractionation coupled to an inductively coupled plasma mass spectrometer was applied to detect the source of trace elements based on their speciation profile along a 125-km stretch of a large river in the Canadian boreal forest. Both the concentration and proportion of bound trace elements were increased by tributary inputs: bound As, Co, Fe, Mn, Pb, U, and Zn increased monotonically from upstream to downstream, increasingly resembling the speciation profile of tributaries. Principal component (PC) analysis also revealed tributary contributions of bound Cu, Ni, Th, V, and Y reflecting their higher concentrations in tributaries, and PC scores also increased monotonically from upstream-downstream. Monotonically decreasing concentrations of mainly ionic and small (i.e. <ca. 300 Da) As, Ba, Mo, and U species were also observed from upstream-downstream.

Degradation of naphthenic acid model compounds in aqueous solution by UV activated persulfate: Influencing factors, kinetics and reaction mechanisms

Naphthenic acids (NAs) are one of the constituents of concerns in oil sands process water (OSPW) because of their persistence and recalcitrance. Herein, we investigated the degradation of five model NA compounds by UV-activated persulfate (UV/persulfate) process under medium-pressure UV lamp irradiation at pH 8.0. UV/persulfate process showed higher degradation efficiency towards cyclohexanoic acid (CHA) compared to UV/H₂O₂ process under the same experimental conditions. CHA (0.39 mM) was completely removed within 30 min when 2 mM persulfate was used as oxidant, while more than 60 min were needed for the UV/H₂O₂ process. The removal of CHA decreased from 100% to 10% when 300 mM tert-butyl alcohol (TBA) was used as the scavenger, indicating that hydroxyl radical (OH) was responsible for the CHA degradation in the UV/persulfate process. Sulfate (SO₄^-) radicals reacted slowly with CHA in the UV/persulfate process with a second-order rate constant of k = 5.3 × 10⁷ M^-1s^-1. Relative kinetics studies using binary mixtures of model NA compounds showed similar structure-reactivity to that under UV/H₂O₂ process. NAs with long carbon chain, cyclic ring, and aromatic ring were more reactive in the UV/persulfate process. The presence of high concentration of chloride ions dramatically inhibited the reaction. The OH radicals in the UV/persulfate process were generated by capturing OH^- in solutions, as evidenced by the decrease of the pH value from 8.0 to 2.8 before and after treatment, respectively, in a pure water matrix. Primary intermediate products (oxy-CHA, hydroxyl-CHA, and dihydroxyl-CHA) of UV/persulfate process were confirmed by UPLC-MS.

Total Mercury and Methylmercury in Lake Water of Canada's Oil Sands Region

Increased delivery of mercury to ecosystems is a common consequence of industrialization, including in the Athabasca Oil Sands Region (AOSR) of Canada. Atmospheric mercury deposition has been studied previously in the AOSR; however, less is known about the impact of regional industry on toxic methylmercury (MeHg) concentrations in lake ecosystems. We measured total mercury (THg) and MeHg concentrations for five years from 50 lakes throughout the AOSR. Mean lake water concentrations of THg (0.4-5.3 ng L^-1) and MeHg (0.01-0.34 ng L^-1) were similar to those of other boreal lakes and <5% of all samples exceeded Provincial water quality guidelines. Lakes with the highest THg concentrations were found >100 km northwest of oil sands mines and received runoff from geological formations high in metals concentrations. MeHg concentrations were highest in those lakes, and in smaller productive lakes closer to oil sands mines. Simulated annual average direct deposition of THg to sampled lakes using an atmospheric chemical transport model showed <2% of all mercury deposited to sampled lakes was emitted from oil sands activities. Consequently, spatial patterns of mercury in AOSR lakes were likely most influenced by watershed and lake conditions, though mercury concentrations in these lakes may be perturbed with future development and climatic change.

Estimating bioaccessibility of trace elements in particles suspended in the Athabasca River using sequential extraction

Employing protocols developed for polar snow and ice, water samples were collected upstream, midstream and downstream of open pit bitumen mines and upgraders along the Lower Athabasca River (AR). The purpose was to: i) estimate the bioaccessibility of trace elements associated with particulate matter in the AR using sequential extraction, and ii) determine whether their forms have been measurably impacted by industrial activities. Of the trace metals known to be enriched in bitumen (V, Ni, Mo and Re), a substantial proportion of V (78-93%) and Ni (35-81%) was found in the residual fraction representing stable minerals. In contrast, Mo and Re were partitioned mainly into more reactive forms (water soluble, acid extractable, reducible and oxidisable). Comparing the non-residual fractions in upstream versus downstream sites, only water soluble Re was significantly (P = 0.005) greater downstream of industry. In respect to the potentially toxic chalcophile elements (Cu, Pb and Tl), no measurable change was observed in Cu and Pb distribution in upstream versus downstream sites. Only residual Tl was found at upstream and midstream sites, whereas a significant proportion of Tl was also present in the reducible fraction in downstream sites. Overall, a greater proportion of trace metals in the residual fraction at midstream sites appears to be due to inputs of atmospheric dust, clearly evident in microscopic images: energy dispersive spectroscopy and x-ray diffraction analyses showed that these particles were predominantly silicates, which are assumed to have limited bioaccessibility.

Elevated exposure, uptake and accumulation of polycyclic aromatic hydrocarbons by nestling tree swallows (Tachycineta bicolor) through multiple exposure routes in active mining-related areas of the Athabasca oil sands region

In the Athabasca Oil Sands (OS) Region, the exposure (by air, water, diet), uptake and deposition of polycyclic aromatic compounds (PACs), including parent and alkylated hydrocarbons (PAHs) and dibenzothiophenes (DBTs), was assessed in nestling tree swallows (Tachycineta bicolor) at mining-related (OS1, OS2) and reference (REF) sites. The OS sites did not receive oil-sands processed waters (OSPW) and were ≥60km from the reference sites. Most of the 42 PACs (≤98%) were detected in all matrices. Swallows at the OS sites were exposed to higher air and water concentrations of individual PAC congeners, ΣPACs, Σparent-PAHs, Σalkyl-PAHs and ΣDBTs. Compared to reference nestlings (ΣPACs: 13-27ng/g wet weight (ww)), PACs were significantly higher in OS nestlings (31-106ng/gww) that also accumulated higher concentrations of major PAHs (i.e., naphthalene, C1-naphthalene, C2-naphthalene, C1-fluorenes, C2-fluorenes, C1-phenanthrenes) measured in 60% of nestlings. Uptake and deposition of PAHs in the birds' muscle was related to diet (δ¹⁵N: C1-naphthalenes, C2-naphthalenes, C1-fluorenes), water (C1-phenanthrenes), and air through inhalation and feather preening (C1-fluorenes), but fecal concentrations were not well explained by diet or environmental concentrations. While PAH concentrations were much higher in muscle than feces, they were highly correlated (p≤0.001 for all). Thus feces may represent a non-lethal method for characterizing PAH exposure of birds, with muscle characterizing accumulation and sources of PAH exposure. Tree swallows in the Athabasca OS Region are exposed to many PACs, accumulating higher concentrations when developing in close proximity to mining activity through diet, aerial deposition and mining-impacted freshwater sources (e.g., wetlands).

Sources of Methylmercury to Snowpacks of the Alberta Oil Sands Region: A Study of In Situ Methylation and Particulates

Snowpacks in the Alberta Oil Sands Region (AOSR) of Canada contain elevated loadings of methylmercury (MeHg; a neurotoxin that biomagnifies through foodwebs) due to oil sands related activities. At sites ranging from 0 to 134 km from the major AOSR upgrading facilities, we examined sources of MeHg by quantifying potential rates of MeHg production in snowpacks and melted snow using mercury stable isotope tracer experiments, as well as quantifying concentrations of MeHg on particles in snowpacks (pMeHg). At four sites, methylation rate constants were low in snowpacks (k_m = 0.001-0.004 d^-1) and nondetectable in melted snow, except at one site (k_m = 0.0007 d^-1). The ratio of methylation to demethylation varied between 0.3 and 1.5, suggesting that the two processes are in balance and that in situ production is unlikely an important net source of MeHg to AOSR snowpacks. pMeHg concentrations increased linearly with distance from the upgraders (R² = 0.71, p < 0.0001); however, snowpack total particle and pMeHg loadings decreased exponentially over this same distance (R² = 0.49, p = 0.0002; R² = 0.56, p < 0.0001). Thus, at near-field sites, total MeHg loadings in snowpacks were high due to high particle loadings, even though particles originating from industrial activities were not MeHg rich compared to those at remote sites. More research is required to identify the industrial sources of snowpack particles in the AOSR.

Peat Bogs Document Decades of Declining Atmospheric Contamination by Trace Metals in the Athabasca Bituminous Sands Region

Peat cores were collected from five bogs in the vicinity of open pit mines and upgraders of the Athabasca Bituminous Sands, the largest reservoir of bitumen in the world. Frozen cores were sectioned into 1 cm slices, and trace metals determined in the ultraclean SWAMP lab using ICP-QMS. The uppermost sections of the cores were age-dated with ²¹⁰Pb using ultralow background gamma spectrometry, and selected plant macrofossils dated using ¹⁴C. At each site, trace metal concentrations as well as enrichment factors (calculated relative to the corresponding element/Th ratio of the Upper Continental Crust) reveal maximum values 10 to 40 cm below the surface which shows that the zenith of atmospheric contamination occurred in the past. The age-depth relationships show that atmospheric contamination by trace metals (Ag, Cd, Sb, Tl, but also V, Ni, and Mo which are enriched in bitumen) has been declining in northern Alberta for decades. In fact, the greatest contemporary enrichments of Ag, Cd, Sb, and Tl (in the top layers of the peat cores) are found at the control site (Utikuma) which is 264 km SW, suggesting that long-range atmospheric transport from other sources must be duly considered in any source assessment.

Silver-Ion Solid Phase Extraction Separation of Classical, Aromatic, Oxidized, and Heteroatomic Naphthenic Acids from Oil Sands Process-Affected Water

The separation of classical, aromatic, oxidized, and heteroatomic (sulfur-containing) naphthenic acid (NA) species from unprocessed and ozone-treated oil sands process-affected water (OSPW) was performed using silver-ion (Ag-ion) solid phase extraction (SPE) without the requirement of pre-methylation for NAs. OSPW samples before SPE and SPE fractions were characterized using ultra performance liquid chromatography ion mobility time-of-flight mass spectrometry (UPLC-IM-TOFMS) to corroborate the separation of distinct NA species. The mass spectrum identification applied a mass tolerance of ±1.5 mDa due to the mass errors of NAs were measured within this range, allowing the identification of O2S-NAs from O2-NAs. Moreover, separated NA species facilitated the tandem mass spectrometry (MS/MS) characterization of NA compounds due to the removal of matrix and a simplified composition. MS/MS results showed that classical, aromatic, oxidized, and sulfur-containing NA compounds were eluted into individual SPE fractions. Overall results indicated that the separation of NA species using Ag-ion SPE is a valuable method for extracting individual NA species that are of great interest for environmental toxicology and wastewater treatment research, to conduct species-specific studies. Furthermore, the separated NA species on the milligram level could be widely used as the standard materials for environmental monitoring of NAs from various contamination sites.