Evaluating officially reported polycyclic aromatic hydrocarbon emissions in the Athabasca oil sands region with a multimedia fate model

Unveiling underutilized public data on Canadian industrial pollutant transfers and disposals

Over 50 nations worldwide have Pollutant Release and Transfer Registers (PRTRs), including Canada's National Pollutant Release Inventory (NPRI), which comprise large public datasets of chemical releases to air, water, and land and also transfers to various on and off-site waste management practices. These inventories aim to support a myriad of audiences in pollution-related decision-making. While the Organisation for Economic Cooperation and Development (OECD) framed a role for PRTRs as indicators for Sustainable Development Goal (SDG) 12 - the sound management of chemicals and wastes, research to date has focused only on air and water releases, omitting vast PRTR data on pollutant transfers to waste management. For Canada's NPRI, 30 years of waste management transfers data on 250+ chemicals has been collected but rarely used in environmental research. Here we show how this overlooked NPRI data may be used to inform snapshots and trends in progress towards SDG 12 using the OECD's framework. Results show that over 28 million (M) tonnes (t) of NPRI pollutants have been transferred from industrial facilities to various waste management practices from 2006 to 2021, of which ~10M t were transferred off-site for waste management operations both within and outside Canada. Time trends show pollutant transfer quantities are increasing, driven by on-site disposals to tailings and waste rock management (of mainly phosphorous, manganese and other metals) and underground injection (of mainly hydrogen sulphide). New route maps reveal that interprovincial and international pollutant transfers are common, and that chain of custody analyses is a burgeoning opportunity but hampered by data limitations. The findings create a state of the knowledge launching point for mainstreaming the use of this overlooked data from both Canada's NPRI and PRTRs around the world, to better track both national and international progress towards sound management of chemicals in waste and SDG 12.Implications: Pollutant transfers and disposals data has not previously been widely used in environmental research. This paper shows how it can be, in the context of SDG 12. Doing so can inspire uptake by researchers and a range of other public users, both strengthening the justification for collecting this data, and bolstering public participation in environmental decision-making from a local to global scale. Doing so also provides the foundation for more in-depth analysis on the domestic and international transboundary movement of Canadian industrial pollutants in waste in the lens of SDG 12 - a topic that was beyond scope here but addressed elsewhere.

Oil sands process-affected water composition effect on Henry's law constants for polycyclic aromatic compounds: Theory and experiment

Oil sands process-affected water (OSPW) is a source of atmospheric emission for polycyclic aromatic compounds (PACs), compounds known to have toxic effects on humans. Estimating emissions and assessing the chemical fate of PACs requires measured or predicted physical-chemical properties such as Henry's law constants (H), that can be used to predict chemical transfer into the atmosphere. OSPW is a complex water-based mixture that is highly variable in composition and nature and contains both organic and inorganic ions. This study uses COSMO-RS solvation theory to estimate and compare Henry's law constants for a set of PACs in both water and theoretically modelled OSPW, to assess the expected deviation that occurs from pure water H values due to the ionic content within OSPW. Experimental measurements of Henry's law constants for PACs in pure water and OSPW using EVA-coated passive dosing and sampler beads were also made in support of our theoretical predictions. For the theory work, OSPW composition data for the Athabasca oil sands in Alberta were used to model a simulated OSPW environment with realistic sodium, chloride, fluoride, sulfate, potassium, bicarbonate, and naphthenic acid concentrations. Theory results indicate that the combined presence of these ions at OSPW concentrations has a negligible effect on H values, causing on average a 3% or 0.014 log unit deviation. By comparison, temperature has a much larger influence on H values, with estimations showing an average 0.20 log unit increase for a 5 °C increase in temperature. The experimental results demonstrate that Henry's law constants can be accurately and precisely measured with this technique in pure water but with less precision in OSPW. Nevertheless, the experimental results support the conclusion that Henry's law constants for OSPW can be accurately estimated assuming a pure water phase.

Measurement constrained emission estimates of alkylated polycyclic aromatic hydrocarbons in the Canadian Athabasca oil sands region

Alkylated polycyclic aromatic hydrocarbons (APAH) are important contaminants of crude oil production and exhibit similar toxicity to their parent compounds. This study developed an emission inventory of APAH in a major oil sands development region of Alberta, Canada, and validated the inventory with ambient concentration measurements through dispersion modeling. The initial estimate of regional total annual emissions of 21 APAH species was 362 tonnes/year in the last decade, of which 309 and 53 tonnes/year were in particle-bound and gas-phase APAH, respectively. Fugitive dust from oil sands mining activities is the primary source of particle-bound APAH, emitting 274 tonnes/year. Other major sources of APAH include point sources (31), tailings ponds (21), anthropogenic fuel consumption from mine fleet (17), and local transportation (13). The group of species with highest emissions was C1-C4 alkylnaphthalenes (53%), followed by C1-C4 alkylphenanthrenes/anthracenes (19%), C1-C4 fluorenes (13%), and C1-C4 fluoranthenes/pyrenes and C1-C4 benz[a]anthracenes/chrysene/triphenylenes (7% each). CALPUFF dispersion modeling was performed using the APAH emissions as model input. The model-predicted annual average ambient APAH concentrations at 17 monitoring sites were 1%-52% (19% on average) lower than the measurements. Inverse dispersion modeling was then applied to adjust APAH emissions higher by 19% for each of the 21 APAH species, which resulted in a revised estimate of APAH emissions to 431 tonnes/year. With the revised emissions as model input, model bias in the predicted ambient concentration was reduced from -19% to -8%. The model results showed the highest concentrations of APAH were near tailings ponds and open mining faces and downwind areas, with total APAH concentrations being higher than 50 ng/m3.

Sediment record of polycyclic aromatic compounds and black carbon over the last ~400 years in Sanjiaolongwan Maar Lake, northeast China

Fuel usage is an important catalyst for socio-economic development and human well-being. Human activities have resulted in significant increases in emissions from biomass burning (BB) and fossil fuel (FF) combustion which have significantly adversely affected human, ecosystem, and planetary health in this era of the Anthropocene. Sanjiaolongwan Maar Lake (SJLW), as a typical crater lake, uniquely receives atmospheric deposition from long-distance transport, and thus, its sediments reflect environmental change and human impacts on a broad scale. In this study, the concentrations and compositions of combustion products, including polycyclic aromatic compounds (PACs, i.e., polycyclic aromatic hydrocarbons (PAHs) and their oxygenated (OPAHs) and nitrogen heterocyclic derivatives (AZAs)) and black carbon (BC and its constituents char and soot), in SJLW over the past 400 years were investigated. The results showed that the PACs and soot concentrations and fluxes in SJLW have rapidly increased since 1950. The concentrations of the total PACs increased ~4 times after the 1950s. Such a fast increase is consistent with the rapid industrialization after the establishment of the People's Republic of China (PRC), which has further accelerated beginning with the implementation of the reform and opening up policy of the PRC in 1978. Moreover, the variations in the compositions of PACs, as well as the decrease in the char/soot ratio, demonstrate a transition in energy usage from BB to FF combustion. The decrease in the benzo[e]pyrene/benzo[a]pyrene ratio indicated an increase in local emissions (because of increasing industrialization in northeast China). The temporal profile of perylene concentrations, fluxes, and perylene/5-ring PAHs ratios strongly suggest that perylene mainly originated from non-pyrogenic sources. The records of PACs and BC in SJLW offer valuable perspectives on human impacts and provide important references for the start of the Anthropocene.

Sources, spatial-distributions and fluxes of PAH-contaminated dusts in the Athabasca oil sands region

Atmospheric deposition of polycyclic aromatic hydrocarbons (PAHs) has increased in northern Alberta, Canada, due to industrial development in the Athabasca oil sands region (AOSR). However, the sources, summertime deposition fluxes and associated spatial patterns are poorly characterized, and the magnitude of contamination has not been directly contrasted with comparable measurements around large Canadian cities. PAHs were measured in Sphagnum moss collected from 30 bogs in the AOSR and compared with reference moss collected from various remote, rural and near-urban sites in Alberta and Ontario. At all 39 locations, strong correlations between depositional fluxes of PAHs and accumulation rates of ash (n = 117, r = 0.877, p < 0.001) implied that the main source of PAHs to moss was atmospheric deposition of particles. Average PAH concentrations at near-field AOSR sites (mean [SD], 62.4 [24.3] ng/g) were significantly higher than at far-field AOSR sites (44.9 [20.8] ng/g; p = 0.038) or the 7 reference sites in Alberta (20.6 [3.5] ng/g; p < 0.001). In fact, average PAH concentrations across the entire AOSR (7,850 km2) were approximately twice as high as in London, Ontario, or near petroleum upgrading and major traffic corridors in Edmonton, Alberta. A chemical mass balance model estimated that both delayed petcoke (33 % of PAHs) and fine tailings (38 % of PAHs) were the major sources of PAHs in the AOSR. Over the 2015 summer growing season, we estimate that 101-110 kg of PAHs (on 14,300-17,300 tonnes of PAH-containing dusts) were deposited to the AOSR within a 50 km radius of surface mining. Given that the highest PAH deposition was to the northern quadrant of the AOSR, which includes the First Nations community of Fort MacKay, further dust control measures should be considered to protect human and environmental health in the region.

Atmospheric emissions estimation of polycyclic aromatic compounds from an oil sands tailings pond using passive air samplers

A mapping study targeting emissions of polycyclic aromatic compounds (PACs) from an oil sands tailings pond was undertaken in the Athabasca Oil sands Region (AOSR). Ten passive air samplers comprising polyurethane foam (PUF) disks were deployed around the perimeter of Suncor Tailings Pond 2/3 for a five-week period to generate time-integrated concentrations in air for PACs, which included ∑unsubstituted polycyclic aromatic hydrocarbons (PAHs), ∑alkylated PAHs (alk-PAHs), and ∑dibenzothiophenes (DBTs) (both unsubstituted and alkylated). Concentrations in air ranged from 13 to 70, 220-970, and 30-210 ng/m3, respectively, and were elevated in samplers downwind of the tailings pond. PAC emissions to air from the pond were estimated using only the air-side concentration information by applying a simplified Gaussian dispersion model and found to be 896 μg/m2/day. ∑alk-PAHs and ∑DBTs had the highest contribution to the total PAC fluxes (79% and 16%, respectively). This flux estimate for PACs is equivalent to 460 kg on an annual basis and 35 000 kg/year when scaled to represent all tailings ponds in the region. The results generally agree with fluxes estimated from coupled high volume air sampling data and tailings pond water concentrations from the same field study but which are complicated due to uncertainties associated with the use of pure water Henry's Law values for tailings pond water as well as the potential for surface oily films on the tailings ponds to impact water-air exchange of PACs. Overall, these findings support the use of relatively simple and electricity-free PUF disk samplers for mapping and estimating emissions from area sources such as tailings ponds, using only air-side concentration information.

Temporal Trends of Legacy and Current-Use Halogenated Flame Retardants in Lake Ontario in Relation to Atmospheric Loadings, Sources, and Environmental Fate

Since the phase-out of polybrominated diphenyl ethers (PBDEs), large amounts of alternative halogenated flame retardants (AHFRs) have been introduced to the market. Due to their persistence and toxicity, halogenated flame retardants (HFRs) have become a concern for the ecosystem and human health. However, there remains limited assessment of the atmospheric loadings, sources, and environmental fate of HFRs in Lake Ontario, which receives urban-related inputs and cumulative chemical inputs from the upstream Great Lakes from Canada and the United States. We combined long-term measurements with a modified multimedia model based on site-specific environmental parameters from Lake Ontario to understand the trends and fate of HFRs. All HFRs were detected in the air, precipitation, lake trout, and herring gull egg samples throughout the sampling periods. General decreasing trends were found for PBDEs, while the temporal trends for AHFRs were not clear. Physical-chemical properties and emissions significantly influence the levels, profiles, and trends. Using the probabilistic modeling, HFR concentrations in lake water and sediment were predicted to be close to the measurement, suggesting a good performance for the modified model. The loadings from tributaries and wastewater effluent were the primary input pathways. Transformations in the water and sedimentation were estimated to be the dominant output pathway for the three HFRs.

A multi-pathway exposure assessment for polycyclic aromatic hydrocarbons among residents in the Athabasca oil sands region, Canada

Due to increasing emissions from ongoing development of the oil sands in Northern Alberta, Canada, there is concern that local residents and organisms are experiencing elevated exposures to hazardous contaminants. We modified an existing human bioaccumulation model (ACC-Human) to represent the local food chain in the Athabasca oil sands region (AOSR), the focus of oil sands development in Alberta. We used the model to assess the potential exposure to three polycyclic aromatic hydrocarbons (PAHs) among local residents that have a high intake of locally sourced traditional foods. To place these estimates into context, we complemented them with estimated PAH intake through market foods and smoking. Our approach was able to produce realistic body burdens of the PAHs in aquatic and terrestrial wildlife and in humans, both in magnitude and with respect to the relative difference between smokers and non-smokers. Over the model simulation period (1967-2009), market food was the dominant dietary exposure route for phenanthrene and pyrene, while local food, and in particular local fish, dominated the intake of benzo[a]pyrene. Exposure to benzo[a]pyrene therefore was also predicted to increase over time in concert with expanding oil sands operations. Those smoking at the average rate of Northern Albertans take in an additional amount of all three PAHs that is at least as large as dietary intake. Estimated daily intake rates are below toxicological reference thresholds for all three PAHs. However, daily intake of BaP in adults is only ∼20 fold below those thresholds and is predicted to increase. Key uncertainties in the assessment included the effect of food preparation on the PAH content in food (e.g., smoking of fish), the limited availability of market food contamination data specific to Canada, and the PAH content of the vapor phase of first-hand cigarette smoke. Considering the satisfactory model evaluation, ACC-Human AOSR should be suited to making predictions of future contaminant exposure based on development scenarios in the AOSR or in response to potential emission reduction efforts. It should also be applicable to other organic contaminants of concern released by oil sands operations.

The last 25 years of research on bioflocculants for kaolin flocculation with recent trends and technical challenges for the future

The generation of kaolin-containing wastewater is an inevitable consequence in a number of industries including mining, wastewater treatment, and bitumen processing. In some cases, the production of kaolin tailings waste during the production of bitumen or phosphate is as high as 3 times greater than the actual produced product. The existing inventory of nearly five billion barrels of oil sands tailings alone represents a massive storage and reclamation challenge, as well as a significant economic and environmental liability. Current reclamation options like inorganic coagulants and organic synthetic polymers may settle kaolin effectively, but may themselves pose an additional environmental hazard. Bioflocculants are an emerging alternative, given the inherent safety and biodegradability of their bio-based compositions. This review summarizes the different research attempts towards a better bioflocculant of kaolin, with a focus on the bioflocculant source, composition, and effective flocculating conditions. Bacillus bacteria were the most prevalent single species for bioflocculant production, with wastewater also hosting a large number of bioflocculant-producing microorganisms while serving as an inexpensive nutrient. Effective kaolin flocculation could be obtained over a broad range of pH values (1-12) and temperatures (5-95°C). Uronic acid and glutamic acid were predominant sugars and amino acids, respectively, in a number of effective bioflocculants, potentially due to their structural and charge similarities to effective synthetic polymers like polyacrylamide. Overall, these results demonstrate that bioflocculants can be produced from a wide range of microorganisms, can be composed of polysaccharides, protein or glycoproteins and can serve as effective treatment options for kaolin. In some cases, the next obstacle to their wide-spread application is scaling to industrially relevant volumes and their deployment strategies.

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 re-analysis and review of elemental and polycyclic aromatic compound deposition in snow and lake sediments from Canada's Oil Sands Region integrating industrial performance and climatic variables

Much of the research from Canada's oil sands region (OSR) shows contaminants of concern (CoCs) throughout the ambient environment surrounding the industrial facilities. While there are some well-established sources of the CoCs, there is also spatial and temporal variability suggesting activity intensity, changes in technology, types and amounts of fuels combusted at the facilities, and climate may affect the results of deposition studies. This study re-analysed published data on the deposition of elements and polycyclic aromatic compounds (PACs) in snow and the sediments of some lakes by incorporating production data from facilities and climate. Using the Elastic Net (EN) regularized regression, variables describing potential associations between facility-specific activity and climate on the deposition of CoCs were identified. Among the selected variables, the combustion of delayed petroleum coke at the Suncor Basemine was associated with the deposition of CoCs, including elements in snow and in some lakes. Similarly, combustion of petroleum coke at Syncrude Mildred Lake was also identified in some models. In both cases, the effects of petroluem coke combustion are likely associated with the emission and deposition of fly ash. The mass of stored petroleum coke was not selected in snow CoC models, but the speed of the wind was a common driver for PACs. However, the mass of stockpiled petcoke was more closely associated with both elements and PACs in lake sediments. While the potential influence of other variables on the occurrence of CoCs in the OSR was also identified, including the production of crude bitumen and synthetic crude, the use of process and natural gases, temperature, and precipitation, these analyses support much of the earlier work and provides additional nuance. While more work is required, these results suggest facility-specific production and climatic data can be coupled with existing approaches to improve the identification of sources of CoCs in Canada's OSR and practices associated with their release.

Microsatellite mutation frequencies in river otters (Lontra Canadensis) from the Athabasca Oil Sands region are correlated to polycyclic aromatic compound tissue burden

Mining activities in the Athabasca oil sands region (AOSR) have contributed to an increase of polycyclic aromatic compounds (PACs) locally. However, many PACs found in the AOSR, and the combined effects of PAC mixtures have not been evaluated for genotoxicity in wildlife. Here, we examine whether mutation frequencies in AOSR river otters are correlated to PAC tissue burdens. We used single-molecule polymerase chain reaction (SM-PCR) to measure the mutant frequency of unstable DNA microsatellite loci in the bone marrow of wild river otters (n = 11) from the AOSR. Microsatellite mutation frequencies were regressed against liver PAC burden (total, low/high molecular weight [LMW/HMW], and parent/alkylated PACs), and to the distances from where the samples were collected to nearby bitumen upgraders. We found that microsatellite mutation frequency was positively correlated with total liver PAC burden. LMW and alkylated PACs were detected at higher levels and had a stronger positive relationship with mutation frequency than HMW (alkylated and parent) PACs. There were no significant relationships detected between mutation frequency and LMW parent PACs or the distance from bitumen upgraders. Furthermore, pyrogenic and petrogenic signatures suggest PACs in animals with high mutation frequencies were associated with combustion processes; although further investigation is warranted, due to limitations of diagnostic ratio determination with biotic models. Our findings support the hypothesis that PACs found in the AOSR increase mutation frequency in wildlife. Further investigation is required to determine if the elevated PAC levels associated with higher mutation frequency are due to natural exposure or elevated human activity.

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).

Polycyclic aromatic compounds in the Canadian Environment: Aquatic and terrestrial environments

Polycyclic aromatic compounds (PACs) are ubiquitous across environmental media in Canada, including surface water, soil, sediment and snowpack. Information is presented according to pan-Canadian sources, and key geographical areas including the Great Lakes, the Alberta Oil Sands Region (AOSR) and the Canadian Arctic. Significant PAC releases result from exploitation of fossil fuels containing naturally-derived PACs, with anthropogenic sources related to production, upgrading and transport which also release alkylated PACs. Continued expansion of the oil and gas industry indicates contamination by PACs may increase. Monitoring networks should be expanded, and include petrogenic PACs in their analytical schema, particularly near fuel transportation routes. National-scale roll-ups of emission budgets may not expose important details for localized areas, and on local scales emissions can be substantial without significantly contributing to total Canadian emissions. Burning organic matter produces mainly parent or pyrogenic PACs, with forest fires and coal combustion to produce iron and steel being major sources of pyrogenic PACs in Canada. Another major source is the use of carbon electrodes at aluminum smelters in British Columbia and Quebec. Temporal trends in PAC levels across the Great Lakes basin have remained relatively consistent over the past four decades. Management actions to reduce PAC loadings have been countered by increased urbanization, vehicular emissions and areas of impervious surfaces. Major cities within the Great Lakes watershed act as diffuse sources of PACs, and result in coronas of contamination emanating from urban centres, highlighting the need for non-point source controls to reduce loadings.

Alkylated polycyclic aromatic hydrocarbons are the largest contributor to polycyclic aromatic compound concentrations in traditional foods of the Bigstone Cree Nation in Alberta, Canada

Rising global demand for energy promotes extensive mining of natural resources, such as oil sands extractions in Alberta, Canada. These extractive activities release hazardous chemicals into the environment, such as polycyclic aromatic compounds (PACs), which include the parent polycyclic aromatic hydrocarbons (PAHs), alkylated PAHs, and sulfur-containing heterocyclic dibenzothiophenes (DBTs). In areas adjacent to industrial installations, Indigenous communities may be exposed to these PACs through the consumption of traditional foods. Our objective was to evaluate and compare the concentrations of total PACs (∑PAC), expressed as the sum of the 16 U.S. EPA priority PAHs (∑PAH), 49 alkylated PAHs (∑alkyl-PAH), and 7 DBTs (∑DBT) in plant and animal foods collected in 2015 by the Bigstone Cree Nation in Alberta, Canada. We analyzed 42 plant tissues, 40 animal muscles, 5 ribs, and 4 pooled liver samples. Concentrations of ∑PAC were higher in the lichen, old man's beard (Usnea spp.) (808 ± 116 ng g^-1 w.w.), than in vascular plants, and were also higher in smoked moose (Alces alces) rib (461 ± 120 ng g^-1 w.w.) than in all other non-smoked animal samples. Alkylated-PAHs accounted for between 63% and 95% of ∑PAC, while the concentrations of ∑PAH represented 4%-36% of ∑PAC. Contributions of ∑DBT to ∑PAC were generally lowest, ranging from <1% to 14%. While the concentrations of benzo(a)pyrene (B[a]P) and ∑PAH4 (∑benzo[a]anthracene, chrysene, benzo[b]fluoranthene, and B[a]P) in all samples were below guideline levels for human consumption as determined by the European Commission, guideline levels for the more prevalent alkylated PAHs are not available. Given the predominance of alkylated PAHs in all food samples and the potentially elevated toxicity relative to parent PAHs of this class of PACs, it is critical to consider a broader range of PACs other than just parent PAHs in research conducted close to oil sands mining activities.

Polycyclic aromatic compounds (PACs) in the Canadian environment: Ambient air and deposition

Polycyclic aromatic compounds (PACs) in Canadian air and deposition were examined at the national scale for the first time in over twenty-five years. Air concentrations spanned four orders of magnitude, and were highest near industrial emitters and lowest in the Arctic. Declines in unsubstituted PAHs were observed at locations close to industrial facilities that had reduced emissions, but trends elsewhere were modest or negligible. Retene concentrations are increasing at several locations. Ambient concentrations of benzo[a]pyrene exceeded Ontario's health-based guideline in many urban/industrial areas. The estimated toxicity of the ambient PAC mixture increased by up to a factor of six when including compounds beyond the US EPA PAHs. Knowledge of PAC deposition is limited to the Laurentian Great Lakes and Athabasca Oil Sands regions. The atmosphere remained a net source of PAHs to the Great Lakes, though atmospheric inputs were decreasing with halving times of 26-30 years. Chemical transport modelling substantially overestimated wet deposition, but model performance is unknown for dry deposition. Sources from Asia, Europe and North America contributed to Arctic and Sub-Arctic concentrations, whereas transboundary or long-range transport have not been assessed outside Canada's north. Climate-related impacts from re-emission and forest fires were implicated in maintaining air concentrations in the high Arctic that were not consistent with global emissions reductions. Industrial emission decreases were substantial at the national scale, but their influence on the environment was limited to areas near relevant facilities. When examined through the lens of ambient levels at the local scale, evidence suggested that contributions from residential wood combustion and motor vehicles were smaller and larger, respectively, than those reported in national inventories. Future work aimed at characterizing PACs beyond the EPA PAHs, improving measurement coverage, elucidating deposition phenomena, and refining estimates of source contributions would assist in reducing remaining knowledge gaps about PACs in Canada.

Fugitive emissions of polycyclic aromatic compounds from an oil sands tailings pond based on fugacity and inverse dispersion flux calculations

Alberta's oil sands tailings ponds are suspected to be a source of fugitive emissions of polycyclic aromatic compounds (PACs) to the atmosphere. Here we report, for the first time, fluxes of 6 parent and 21 alkylated PACs based on the measured co-located air and water concentrations using a two-film fugacity-based model (FUG), an inverse dispersion model (DISP) and a simple box model (BOX). Air samples were collected at the Suncor Tailings Pond 2/3 using a high volume air sampler from the "pond" and towards the pond ("non-pond") directions separately. Mean ∑₂₇PACs in air from the "pond" direction was greater than the "non-pond" direction by a factor of 17. Water-air fugacity ratio of 20 PACs quantifiable in water indicated net volatilization from water. Dispersion and box model results also indicated upward fluxes of 22 PACs. Correlation between the estimated flux results of BOX and DISP model was statistically significant (r = 0.99 and p < 0.05), and correlation between FUG and DISP results ranged from 0.54 to 0.85. In this first-ever assessment of PAC fluxes from tailings pond, the three models confirmed volatilization fluxes of PACs indicating Suncor Tailings Pond 2/3 is a source of PAC emissions to the atmosphere. This study addressed a key data gap identified in the Joint Oil Sands Monitoring Emissions Inventory Compilation Report (Government of Alberta and Canada, 2016) which is the lack of consistent real-world tailings pond fugitive emission monitoring of organic chemicals. Our findings highlight the need for measurements from other tailings ponds to determine their overall contribution in releasing PACs to the atmosphere. This paper presents a practical method for estimating PAC emissions from other tailings ponds, which can provide a better understanding of these fugitive emissions, and thereby help to improve the overall characterization of emissions in the oil sands region.

Polycyclic aromatic compounds (PACs) in the Canadian environment: Sources and emissions

Twenty-five years after the first look at polycyclic aromatic compounds (PACs) in Canada, this article presents current knowledge on Canadian PAC emission sources. The analysis is based on national inventories (the National Pollutant Release Inventory (NPRI) and the Air Pollutant Emissions Inventory (APEI)), an analysis of Canadian forest fires, and several air quality model-ready emissions inventories. Nationally, forest fires continue to dominate PAC emissions in Canada, however there is uncertainty in these estimates. Though forest fire data show a steady average in the total annual area burned historically, an upward trend has developed recently. Non-industrial sources (home firewood burning, mobile sources) are estimated to be the second largest contributor (∼6-8 times lower than forest fires) and show moderate decreases (25%-65%) in the last decades. Industrial point sources (aluminum production, iron/steel manufacturing) are yet a smaller contributor and have seen considerable reductions (90% +) in recent decades. Fugitive emissions from other industrial sources (e.g. disposals by the non-conventional oil extraction and wastewater sectors, respectively) remain a gap in our understanding of total PAC emissions in Canada. Emerging concerns about previously unrecognized sources such as coal tar-sealed pavement run-off, climate change are discussed elsewhere in this special issue. Results affirm that observations at the annual/national scale are not always reflective of regional/local or finer temporal scales. When determining which sources contribute most to human and ecosystem exposure in various contexts, examination at regional and local scales is needed. There is uncertainty overall in emissions data stemming in part from various accuracy issues, limitations in the scope of the various inventories, and inventory gaps, among others.

Improving Insights on Air Pollutant Mixtures and Their Origins by Enhancing Local Monitoring in an Area of Intensive Resource Development

An "event-based" approach to characterize complex air pollutant mixtures was applied in the Oil Sands region of northern Alberta, Canada. This approach was developed to better-inform source characterization and attribution of the air pollution in the Indigenous community of Fort McKay, within the context of the lived experience of residents. Principal component analysis was used to identify the characteristics of primary pollutant mixtures, which were related to hydrocarbon emissions, fossil fuel combustion, dust, and oxidized and reduced sulfur compounds. Concentration distributions of indicator compounds were used to isolate sustained air pollution "events". Diesel-powered vehicles operating in the mines were found to be an important source during NO_x events. Industry-specific volatile organic compound (VOC) profiles were used in a chemical mass balance model for source apportionment, which revealed that nearby oil sands operations contribute to 86% of the total mass of nine VOC species (2-methylpentane, hexane, heptane, octane, benzene, toluene, m,p-xylene, o-xylene, and ethylbenzene) during VOC events. Analyses of the frequency distribution of air pollution events indicate that Fort McKay is regularly impacted by multiple mixtures simultaneously, underscoring the limitations of an exceedance-based approach relying on a small number of air quality standards as the only tool to assess risk.

How does climate variability affect water quality dynamics in Canada's oil sands region?

In Canada's oil sands region, classic boreal hydrology (i.e., winter low flow followed by peaks during spring freshet and then summer flow recession) combined with erosion of both natural and anthropogenically-exposed bitumen results in seasonal and inter-annual variability in stream water chemistry. Using data collected from all seasons over three years (2012-2015), we investigated the mechanisms driving spatial and temporal change in the concentration of 26 water quality parameters for six rivers draining Canada's oil sands region. Mantel tests showed a strong spatial aggregation of climatic drivers (average daily precipitation, accumulated precipitation, snow water equivalent) associated with west versus east discharge patterns. Wavelet analysis highlighted unique watershed attributes, in particular the importance of developed area in lowering responsiveness to seasonal precipitation. Concentrations of most chemical parameters (20 of 23) showed distinct temporal patterns that were correlated with seasonal changes in hydrology which, in turn, were related to changes in weather. Comparison of concentrations observed in this study with those reported in the scientific literature for the same watersheds showed 81% of comparisons differed significantly. This was likely due to the short duration of previous field campaigns and thus the sampling of a very narrow window of the annual streamflow regime.

Waste Containment Ponds Are a Major Source of Secondary Organic Aerosol Precursors from Oil Sands Operations

The surface mining and bitumen extraction of oil sands (OS) generates over one million barrels of heavy oil each day in the Alberta Oil Sands Region of Canada. Recent observations suggest that emissions from OS development contribute to secondary organic aerosol (SOA) formation, but the chemical composition, mass fluxes, and sources of those emissions are poorly delineated. Here, we simulated OS extraction and used comprehensive two-dimensional gas chromatography to quantify and characterize direct air emissions, bitumen froth, residual wastewater, and tailings components, ultimately enabling fate modeling of over 1500 chromatographic features simultaneously. During the non-ice cover season, tailings ponds emissions contributed 15 000-72 000 metric tonnes of hydrocarbon SOA precursors, translating to 3000-13 000 tonnes of SOA, whereas direct emissions during the extraction process itself were notably smaller (960 ± 500 tonnes SOA yr^-1). These results suggest that tailings pond waste management practices should be targeted to reduce environmental emissions.

Forest health effects due to atmospheric deposition: Findings from long-term forest health monitoring in the Athabasca Oil Sands Region

Oil sands developments release acidifying compounds (SO₂ and NO₂) with the potential for acidifying deposition and impacts to forest health. This article integrates the findings presented in the Oil Sands Forest Health Special Issue, which reports on the results of 20 years of forest health monitoring, and addresses the key questions asked by WBEA's Forest Health Monitoring (FHM) Program: 1) is there evidence of deposition affecting the environment?, 2) have there been changes in deposition or effects over time?, 3) do acid deposition levels require management intervention?, 4) what are major sources of deposited substances? and 5) how can the program be improved? Deposition of sulphur, nitrogen, base cations (BC), polycyclic aromatic compounds and trace elements decline exponentially with distance from sources. There is little evidence for acidification effects on forest soils or on understory plant communities or tree growth, but there is evidence of nitrogen accumulation in jack pine needles and fertilization effects on understory plant communities. Sulphur, BC and trace metal concentrations in lichens increased between 2008 and 2014. Source apportionment studies suggest fugitive dust in proximity to mining is a primary source of BC, trace element and organic compound deposition, and BC deposition may be neutralizing acidifying deposition. Sulphur accumulation in soils and nitrogen effects on vegetation may indicate early stages of acidification. Deposition estimates for sites close to emissions sources exceed proposed regulatory trigger levels, suggesting a detailed assessment of acidification risk close to the emission sources is warranted. However, there is no evidence of widespread acidification as suggested by recent modeling studies, likely due to high BC deposition. FHM Program evolution should include continued integration with modeling approaches, ongoing collection and assessment of monitoring data and testing for change over time, and addition of monitoring sites to fill gaps in regional coverage.

Application of thermal desorption methods for airborne polycyclic aromatic hydrocarbon measurement: A critical review

Thermal desorption (TD) is a universal solvent-free pre-concentration technique. It is often used to pre-concentrate semi-volatile and volatile organic compounds in various sample types. Polycyclic aromatic hydrocarbons (PAHs) are widespread contaminants from incomplete combustion of organic matter and fossil fuel, which have carcinogenic effects on human health. Conventional methods for determining PAHs, represented by solvent extraction, are gradually being replaced by solvent-free methods, typically the TD technique, because of TD's many advantages, including time savings and environmentally friendly treatment. This work presents an extensive review of the universal methods used to determine PAHs in the atmosphere based on the TD technique. The methods currently used for collection and detection of both gas- and particle-phase PAHs in the air are critically reviewed. In addition, the operating parameters of the TD unit are summarized and discussed. The design shortcomings of existing studies and the problems that researchers should address are presented, and promising alternatives are suggested. This paper also discusses important parameters, such as reproducibility and limit of detection, that form a crucial part of quality assurance. Finally, the limitations and the future prospects of the TD technique for use in airborne PAH analyses are addressed. This is the first review of the latest developments of the TD technique for analysis of PAHs and their derivatives in the atmosphere.

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.

Athabasca oil sands region snow contains efficient micron and nano-sized ice nucleating particles

The Athabasca Oil Sands Region (AOSR) in Alberta, Canada, is an important source of atmospheric pollutants, such as aerosols, that have repercussions on both the climate and human health. We show that the mean freezing temperature of snow-borne particles from AOSR was elevated (-7.1 ± 1.8 °C), higher than mineral dust which freezes at ∼ -15 °C and is recognized as one of the most relevant ice nuclei globally. Ice nucleation of nanosized snow samples indicated an elevated freezing ability (-11.6 ± 2.0 °C), which was statistically much higher than snow-borne particles from downtown Montreal. AOSR snow had a higher concentration (∼2 orders of magnitude) of >100 nm particles than Montreal. Triple quadrupole ICP-(QQQ)-MS/MS analysis of AOSR and Montreal snow demonstrated that most concentrations of metals, including those identified as emerging nanoparticulate contaminants, were much more elevated in AOSR in contrast to Montreal: 34.1, 34.1, 16.6, 5.8, 0.3, 0.1, and 9.4 mg/m³ for Cr, Ni, Cu, As, Se, Cd, and Pb respectively, in AOSR and 1.3, 0.3, 2.0, <0.03, 0.1, 0.03, and 1.2 mg/m³ in Montreal snow. High-resolution Scanning Transmission Electron Microscopy/Energy-dispersive X-ray Spectroscopy (STEM-EDS) imaging provided evidence for various anthropogenic nano-materials, including carbon nanotubes resembling structures, in AOSR snow up to 7-25 km away from major oil sands upgrading facilities. In summary, particles characterized as coming from oil sands are more efficient at ice nucleation. We discuss the potential impacts of AOSR emissions on atmospheric and microphysical processes (ice nucleation and precipitation) both locally and regionally.

Extracts of Passive Samplers Deployed in Variably Contaminated Wetlands in the Athabasca Oil Sands Region Elicit Biochemical and Transcriptomic Effects in Avian Hepatocytes

Recent contaminant monitoring in boreal wetlands situated in Alberta's Athabasca oil sands region revealed increased concentrations of polycyclic aromatic compounds (PACs) in passive sampling devices deployed in wetlands close to bitumen surface mining operations. In this study, graded concentrations of semipermeable membrane device (SPMD) extracts, collected from 4 wetlands with variable burdens of PACs, were administered to chicken and double-crested cormorant (DCCO) embryonic hepatocytes to determine effects on 7-ethoxyresorufin-O-deethylase (EROD) activity and mRNA expression. Concentrations and composition of PACs detected in SPMDs varied among sites, and the proportion of alkyl PACs was greater than parent compounds at all sites. ΣPACs was the highest in SPMDs deployed within 10 km of mining activity (near-site wetland [5930 ng SPMD^-1]) compared to those ∼50 km south (far-site wetland [689 ng SPMD^-1]). Measures of EROD activity and Cyp1a4 mRNA expression allowed the ranking of wetland sites based on aryl hydrocarbon receptor-mediated end points; EROD activity and Cyp1a4 mRNA induction were the highest at the near-site wetland. ToxChip PCR arrays (one chicken and one DCCO) provided a more exhaustive transcriptomic evaluation across multiple toxicological pathways following exposure to the SPMD extracts. Study sites with the greatest PAC concentrations had the most genes altered on the chicken ToxChip (12-15/43 genes). Exposure of avian hepatocytes to SPMD extracts from variably contaminated wetlands highlighted traditional PAC-related toxicity pathways as well as other novel mechanisms of action. A novel combination of passive sampling techniques and high-throughput toxicity evaluation techniques shows promise in terms of identifying hotspots of chemical concern in the natural environment.

Deposition and Source Identification of Nitrogen Heterocyclic Polycyclic Aromatic Compounds in Snow, Sediment, and Air Samples from the Athabasca Oil Sands Region

Polycyclic aromatic compounds (PACs) can have multiple sources in the Athabasca Oil Sands Region (AOSR). The current study was designed to identify and explore the potential of nitrogen heterocyclic PACs (NPACs) as source indicators in snowpack, lake sediment and passive air samples from the AOSR during 2014-2015. Source samples including petroleum coke (petcoke), haul road dust, and unprocessed oil sands were also analyzed. Samples were analyzed using comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry, and liquid chromatography-high resolution Orbitrap mass spectrometry. Over 200 NPACs were identified and classified into at least 24 isomer groups, including alkylated carbazoles, benzocarbazoles, and indenoquinolines. Levels of NPACs in environmental samples decreased with distance from the main developments and with increasing depth in lake sediments but were detected within 50 km from the major developments. The composition profiles of several NPAC isomer classes, such as dimethylcarbazoles, showed that petcoke had a distinct distribution of NPACs compared to the haul road dust and unprocessed oil sands ores and was the most similar source material to near-field environmental samples. These results suggest that petcoke is a major contributing source for the identified NPACs and that these compounds have the potential to be used as source indicators for future research in the AOSR.

Tailings ponds of the Athabasca Oil Sands Region, Alberta, Canada, are likely not significant sources of total mercury and methylmercury to nearby ground and surface waters

Tailings ponds created during industrial bitumen extraction from the Athabasca Oil Sands Region (AOSR), Alberta, Canada, have been shown to contain numerous contaminants, such as polycyclic aromatic compounds and naphthenic acids, and to slowly leak into adjacent ground and surface waters. Despite elevated concentrations of total mercury (THg) in nearby Athabasca River waters downstream of the AOSR developments, to date there are no published studies of THg or methylmercury (MeHg; a potent neurotoxin) in the AOSR tailings ponds. Here we present concentrations of THg and MeHg, as well as various water chemistry parameters, within four AOSR tailings ponds. Concentrations of SO₄^2-, NH₃, Na, and Cl were elevated in tailings ponds relative to nearby freshwaters. Surface water concentrations of THg (filtered: 0.15-0.57 ng/L) and MeHg (unfiltered: <0.02-0.53 ng/L; filtered: <0.02-0.32 ng/L), though, were generally low in the tailings ponds, with the highest concentrations observed in the oldest pond. In the mature fine tailings that settle out in the ponds, concentrations of THg (37.0-197 ng/g) and MeHg (0.10-0.52 ng/g) were also low, with the highest concentrations again in the oldest pond. We calculated that if all the dissolved THg and MeHg potentially leaking annually from the tailings ponds entered the nearby Athabasca River, river THg and MeHg concentrations would increase by only 0.01% and 0.03%, respectively. Overall, these ponds are likely not significant sources of THg or MeHg to nearby ground and surface waters, although due to the potential for Hg methylation to occur in the ponds themselves, other tailings ponds in the AOSR should be monitored to ensure that concentrations of MeHg in them are also low.

Multimedia fate modeling of perfluorooctanoic acid (PFOA) and perfluorooctane sulphonate (PFOS) in the shallow lake Chaohu, China

Freshwater shallow lake ecosystems provide valuable ecological services to human beings. However, these systems are subject to severe contamination from anthropogenic sources. Per- and polyfluoroalkyl substances (PFASs), including perfluorooctanoic acid (PFOA) and perfluorooctane sulphonate (PFOS), are among the contaminants that have received substantial attention, primarily due to abundant applications, environment persistence, and potential threats to ecological and human health. Understanding the environmental behavior of these contaminants in shallow freshwater lake environments using a modeling approach is therefore critical. Here, we characterize the fate, transport and transformation of both PFOA and PFOS in the fifth largest freshwater lake in China (Chaohu) during a two-year period (2013-2015) using a fugacity-based multimedia fate model. A reasonable agreement between the measured and modeled concentrations in various compartments confirms the model's reliability. The model successfully quantifies the environmental processes and identifies the major sources and input pathways of PFOA and PFOS to the Chaohu water body. Sensitivity analysis reveals the critical role of nonlinear Freundlich sorption, which contributes to a variable fraction of the model true uncertainty in different compartments (8.1%-93.6%). Through additional model scenario analyses, we further elucidate the importance of nonlinear Freundlich sorption that is essential for the reliable model performance. We also reveal the distinct composition of emission sources for the two contaminants, as the major sources are indirect soil volatilization and direct release from human activities for PFOA and PFOS, respectively. The present study is expected to provide implications for local management of PFASs pollution in Lake Chaohu and to contribute to developing a general model framework for the evaluation of PFASs in shallow lakes.

Fluorescence-linked immunosorbent assay for detection of phenanthrene and its homolog

A competitive fluorescence-linked immunosorbent assay (FLISA) was developed using rhodamine B isothiocyanate (RBITC) as the model fluorescent dye conjugate monoclonal antibody (McAb) for detection of Phe and its homolog (acenaphthene, fluorene, fluoranthene, pyrene and indeno [1,2,3-cd] pyrene) in water samples. The detection range of the assay for Phe was from 2.10 to 91.95 ng/mL. The limit of detection was 1.05 ng/mL, which was approximately 2-fold lower than that of traditional ic-ELISA. Compared with traditional ic-ELISA, more than 70 min was saved because of only one immunoreaction step was needed to accomplish the assay. The average recoveries of Phe and its homolog from domestic water, contaminated water and natural water were 100.7%, 100.8% and 101.2% respectively. The accuracy and precision of the developed FLISA were validated with GC-MS/MS. There were good correlation between the two methods from tap water, contaminated water and river water samples were 0.9994, 0.9935 and 0.9967, respectively. The results suggested that the proposed FLISA could be a potential alternative format for rapid, sensitive, and quantitative detection of Phe and its homolog in environmental water.

Environmental fate and exposure models: advances and challenges in 21st century chemical risk assessment

Environmental fate and exposure models are a powerful means to integrate information on chemicals, their partitioning and degradation behaviour, the environmental scenario and the emissions in order to compile a picture of chemical distribution and fluxes in the multimedia environment. A 1995 pioneering book, resulting from a series of workshops among model developers and users, reported the main advantages and identified needs for research in the field of multimedia fate models. Considerable efforts were devoted to their improvement in the past 25 years and many aspects were refined; notably the inclusion of nanomaterials among the modelled substances, the development of models at different spatial and temporal scales, the estimation of chemical properties and emission data, the incorporation of additional environmental media and processes, the integration of sensitivity and uncertainty analysis in the simulations. However, some challenging issues remain and require research efforts and attention: the need of methods to estimate partition coefficients for polar and ionizable chemical in the environment, a better description of bioavailability in different environments as well as the requirement of injecting more ecological realism in exposure predictions to account for the diversity of ecosystem structures and functions in risk assessment. Finally, to transfer new scientific developments into the realm of regulatory risk assessment, we propose the formation of expert groups that compare, discuss and recommend model modifications and updates and help develop practical tools for risk assessment.

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.

Polycyclic aromatic hydrocarbons contamination in surface soil of China: A review

This paper reviews the concentration, distribution, source, and potential risk of polycyclic aromatic hydrocarbons (PAHs) in surface soils of China through analysis of data from >6000 surface soil samples in nearly 100 references. The mean value of total 16 PAHs was 730ng·g^-1 in surface soil in China, a relatively lower or moderate level than other countries. Based on the Maliszewska-Kordybach classification criteria, the proportions of heavily contaminated, contaminated, weakly contaminated, and non-contaminated soil samples were 21.4%, 11.9%, 49.5%, and 17.2%, respectively. There was a clear geographical distribution, with concentrations of the total 16 PAHs descending in the following order: Northeast China (1467ng·g^-1)>North China (911ng·g^-1)>East China (737ng·g^-1)>South China (349ng·g^-1)>West China (209ng·g^-1). Moreover, it was found that the PAH concentrations in surface soil in China descended along the urban-suburban-rural gradient. The concentration and distribution of PAHs were mainly related to the degree of economic development, population density, climatic conditions, and soil organic matter, and the divergence of regional economic patterns and climatic conditions was the main reason for the observed PAH distribution in the soils. Traffic emissions, coal and biomass combustion mainly contributed to the PAH contamination of surface soil in China during the process of urbanization and industrialization, and the average Benzo(a)pyrene equivalent concentration of ∑PAH_7c (seven carcinogenic PAHs) was 99ng·g^-1, which indicated the soil samples had a small potential carcinogenic risk. Despite soil pollution being generally low, PAH concentrations in some areas were relatively high, therefore it is necessary to produce strategies, such as establishing effective guidelines and developing environmental-friendly technology to reduce PAH emissions, and prevent further contamination.

Characterizing the PAHs in surface waters and snow in the Athabasca region: Implications for identifying hydrological pathways of atmospheric deposition

The composition of polycyclic aromatic hydrocarbons present in snow and surface waters in the Athabasca Oil Sands Region (AOSR) was characterized in order to identify major contributors to the organics detected in rivers and lakes in the region. PAH concentrations, measured by three monitoring programs in 2011, were used to compare the PAH compositions of snow and surface waters across the AOSR. The 2011 dataset includes total (dissolved+particulate) concentrations of thirty-four parent and alkylated PAH compounds in 105 snow, 272 river, and 3 lake samples. The concentration of PAHs in rivers varies seasonally, with the highest values observed in July. The timing of increases in PAH concentrations in rivers coincides with the high river discharge during the spring freshet, indicating that this major hydrological event may play an important role in delivering PAHs to rivers. However, the composition of PAHs present in rivers during this period differs from the composition of PAHs present in snow, suggesting that direct runoff and release of PAHs accumulated on snow may not be the major source of PAHs to the Athabasca River and its tributaries. Instead, snowmelt may contribute indirectly to increases in PAHs due to hydrological processes such as erosion of stream channels, remobilization of PAH-containing sediments, increased catchment runoff, and snowmelt-induced groundwater inputs during this dynamic hydrologic period. Better understanding of transformations of PAH profiles during transport along surface and subsurface flow paths in wetland-dominated boreal catchments would improve identification of potential sources and pathways in the region. The compositional differences highlight the challenges in identifying the origins of PAHs in a region with multiple potential natural and anthropogenic sources particularly when the potential transport pathways include air, soil and water.

Episodic acidification of 5 rivers in Canada's oil sands during snowmelt: A 25-year record

Episodic acidification during snowmelt is a natural phenomenon that can be intensified by acidic deposition from heavy industry. In Canada's oil sands region, acid deposition is estimated to be as much as 5% of the Canadian total and large tracks of northeastern Alberta are considered acid-sensitive because of extensive peatland habitats with poorly weathered soils. To identify the frequency, duration and severity of acidification episodes during snowmelt (the predominant hydrological period for delivery of priority pollutants from atmospheric oil sands emissions to surface waters), a 25-year record (1989 to 2014) of automated water quality data (pH, temperature, conductivity) was assembled for 3 rivers along with a shorter record (2012-2014) for another 2 rivers. Acidic episodes (pH<7, ANC<0) were recorded during 39% of all 83 snowmelt events. The severity (duration x magnitude) of episodic acidification increased exponentially over the study period (r²=0.56, P<0.01) and was strongly correlated (P<0.01) with increasing maximum air temperature and weakly correlated with regional land development (P=0.06). Concentrations of aluminum and 11 priority pollutants (Sb, As, Be, Cd, Cr, Cu, Pb, Se, Ag, Tl and Zn) were greatest (P<0.01) during low (<6.5) pH episodes, particularly when coincident with high discharge, such that aluminum and copper concentrations were at times high enough to pose a risk to juvenile rainbow trout (Oncorhynchus mykiss). Although low pH (pH<6.5) was observed during only 8% of 32 acidification episodes, when present, low pH typically lasted 10days. Episodic surface water acidification during snowmelt, and its potential effects on aquatic biota, is therefore an important consideration in the design of long-term monitoring of these typically alkaline (pH=7.72±0.05) rivers.

The toxicity of oil sands process-affected water (OSPW): A critical review

Large volumes of oil sands process-affected water (OSPW) are produced by the surface-mining oil sands industry in Alberta. Both laboratory and field studies have demonstrated that the exposure to OSPW leads to many physiological changes in a variety of organisms. Adverse effects include compromised immunological function, developmental delays, impaired reproduction, disrupted endocrine system, and higher prevalence of tissue-specific pathological manifestations. The composition of OSPW varies with several factors such as ore sources, mining process, and tailings management practices. Differences in water characteristics have confounded interpretation or comparison of OSPW toxicity across studies. Research on individual fractions extracted from OSPW has helped identify some target pollutants. Naphthenic acids (NAs) are considered as the major toxic components in OSPW, exhibiting toxic effects through multiple modes of action including narcosis and endocrine disruption. Other pollutants, like polycyclic aromatic hydrocarbons (PAHs), metals, and ions may also contribute to the overall OSPW toxicity. Studies have been conducted on OSPW as a whole complex effluent mixture, with consideration of the presence of unidentified components, and the interactions (potential synergistic or antagonistic reactions) among chemicals. This review summarizes the toxicological data derived from in vitro and in vivo exposure studies using different OSPW types, and different taxa of organisms. In general, toxicity of OSPW was found to be dependent on the OSPW type and concentration, duration of exposures (acute versus sub chronic), and organism studied.

Effects of freeze-thawing cycles on desorption behaviors of PAH-contaminated soil in the presence of a biosurfactant: a case study in western Canada

Many regions in Canada are facing increasing environmental threats posed by oil and gas exploitation and transportation. These contaminated lands are inevitably subjected to seasonal and diurnal freeze-thawing cycles (FTCs). However, knowledge about the effect of FTCs on the behaviours of hydrophobic contaminants during the aging process of soil is limited. This study investigated the desorption characteristics of phenanthrene in aging soils in the presence of the biosurfactant rhamnolipid under diurnal and seasonal FTC treatments. It was found that the presence of rhamnolipid in soil during the aging process was able to increase the desorption efficiency of phenanthrene. In the presence of rhamnolipid above 100 mg L^-1, FTCs could inhibit the sequestration of phenanthrene. Soil moisture and rhamnolipid concentration are two major factors affecting this effect. High moisture content and FTC frequency could lead to lower desorption in the early stage of FTCs due to the increased specific surface area. The sequestration of phenanthrene was less effectively hindered under seasonal FTCs than diurnal FTCs. The results from this study have important implications for understanding the role of surfactants in cold-region soil aging, and for the improvement of site remediation strategies of PAH contaminated soil in cold regions.

Heterocyclic Aromatics in Petroleum Coke, Snow, Lake Sediments, and Air Samples from the Athabasca Oil Sands Region

The aromatic fractions of snow, lake sediment, and air samples collected during 2011-2014 in the Athabasca oil sands region were analyzed using two-dimensional gas chromatography following a nontargeted approach. Commonly monitored aromatics (parent and alkylated-polycyclic aromatic hydrocarbons and dibenzothiophenes) were excluded from the analysis, focusing mainly on other heterocyclic aromatics. The unknowns detected were classified into isomeric groups and tentatively identified using mass spectral libraries. Relative concentrations of heterocyclic aromatics were estimated and were found to decrease with distance from a reference site near the center of the developments and with increasing depth of sediments. The same heterocyclic aromatics identified in snow, lake sediments, and air were observed in extracts of delayed petroleum coke, with similar distributions. This suggests that petroleum coke particles are a potential source of heterocyclic aromatics to the local environment, but other oil sands sources must also be considered. Although the signals of these heterocyclic aromatics diminished with distance, some were detected at large distances (>100 km) in snow and surface lake sediments, suggesting that the impact of industry can extend >50 km. The list of heterocyclic aromatics and the mass spectral library generated in this study can be used for future source apportionment studies.

Differences between measured and reported volatile organic compound emissions from oil sands facilities in Alberta, Canada

Large-scale oil production from oil sands deposits in Alberta, Canada has raised concerns about environmental impacts, such as the magnitude of air pollution emissions. This paper reports compound emission rates (E) for 69-89 nonbiogenic volatile organic compounds (VOCs) for each of four surface mining facilities, determined with a top-down approach using aircraft measurements in the summer of 2013. The aggregate emission rate (aE) of the nonbiogenic VOCs ranged from 50 ± 14 to 70 ± 22 t/d depending on the facility. In comparison, equivalent VOC emission rates reported to the Canadian National Pollutant Release Inventory (NPRI) using accepted estimation methods were lower than the aE values by factors of 2.0 ± 0.6, 3.1 ± 1.1, 4.5 ± 1.5, and 4.1 ± 1.6 for the four facilities, indicating underestimation in the reported VOC emissions. For 11 of the combined 93 VOC species reported by all four facilities, the reported emission rate and E were similar; but for the other 82 species, the reported emission rate was lower than E The median ratio of E to that reported for all species by a facility ranged from 4.5 to 375 depending on the facility. Moreover, between 9 and 53 VOCs, for which there are existing reporting requirements to the NPRI, were not included in the facility emission reports. The comparisons between the emission reports and measurement-based emission rates indicate that improvements to VOC emission estimation methods would enhance the accuracy and completeness of emission estimates and their applicability to environmental impact assessments of oil sands developments.

Potential health risks posed by polycyclic aromatic hydrocarbons in muscle tissues of fishes from the Athabasca and Slave Rivers, Canada

Polycyclic aromatic hydrocarbons (PAHs) are released to the environment from oil sands operations and from natural sources in Alberta, Canada. Concentrations of 16 USEPA priority PAHs were measured in tissues of fishes collected from three locations on the Athabasca River in Alberta and two downstream locations on the Slave River in the Northwest Territories, Canada. A total of 425 individual fish were collected including 89 goldeye (Hiodon alosoides), 93 whitefish (Coregonus clupeaformis), 104 northern pike/jackfish (Esox lucius), 96 walleye (Sander vitreus) and 43 burbot/loche mariah/mariah (Lota lota). Fish were sampled during the summer and fall of 2011 and spring of 2012. Dorsal muscle of fishes from upstream reaches of the Athabasca River, close to oil sands extraction and upgrading activities, contained greater concentrations of individual PAHs than concentrations in muscle of fishes from further downstream in the Slave River. Concentrations of the sum of USEPA indicator PAHs (∑PAHs) in fishes collected in the vicinity of Fort McKay, closest to oil sands activities, varied among seasons with average concentrations ranging from 11 (burbot, summer) to 1.2 × 10² ng/g, wm (burbot, spring) with a mean of 48 ng/g, wm. Concentrations of ∑PAHs in fishes collected in the vicinity of Fort Resolution, the location most distant from oil sands activities, also varied among species and seasons, with average concentrations ranging from 4.3 (whitefish, summer) to 33 ng/g, wm (goldeye, summer) with a mean of 13 ng/g, wm. Significant differences in concentrations of ∑PAHs in muscle were observed within goldeye, jackfish, walleye and whitefish among sites. Health risks posed by PAHs to humans were assessed probabilistically using a B[a]P equivalents approach (B[a]P_eq). The average lifetime risk of additional cancers for humans who consumed fish was deemed to be within an 'acceptable' range of risk (i.e., less than 10^-6).

Stream invertebrate community structure at Canadian oil sands development is linked to concentration of bitumen-derived contaminants

In Canada, the Athabasca oil sands deposits are a source of bitumen-derived contaminants, reaching the aquatic environment via various natural and anthropogenic pathways. The ecological effects of these contaminants are under debate. To quantify the effects of bitumen-derived contaminants we monitored the aquatic exposure of polycyclic aromatic hydrocarbons (PAHs), metals, and naphthenic acids as well as the invertebrate community in the Athabasca River and its tributaries. PAH concentrations over 3 consecutive years were related to discharge and were highest in the year with high autumn rainfall. In the year with the highest PAH concentrations, these were linked with adverse effects on the aquatic invertebrate communities. We observed relative effects of the composition and concentration of contaminants on the invertebrate fauna. This is reflected by the composition and abundance of invertebrate species via the use of the species' traits "physiological sensitivity" and "generation time". Applying the SPEAR approach we observed alterations of community structure in terms of an increased physiological sensitivity and a decrease of generation time for the average species. These effects were apparent at concentrations 100 times below the acute sensitivity of the standard test organism Daphnia magna. To rapidly identify oil sands related effects in the field we designed a biological indicator system, SPEAR_oil, applicable for future routine monitoring.

Air-water partition coefficients for a suite of polycyclic aromatic and other C10 through C20 unsaturated hydrocarbons

The air-water partition coefficients (Kaw) for 86 large polycyclic aromatic hydrocarbons and their unsaturated relatives were estimated using high-level G4(MP2) gas and aqueous phase calculations with the SMD, IEFPCM-UFF, and CPCM solvation models. An extensive method validation effort was undertaken which involved confirming that, via comparisons to experimental enthalpies of formation, gas-phase energies at the G4(MP2) level for the compounds of interest were at or near thermochemical accuracy. Investigations of the three solvation models using a range of neutral and ionic compounds suggested that while no clear preferential solvation model could be chosen in advance for accurate Kaw estimates of the target compounds, the employment of increasingly higher levels of theory would result in lower Kaw errors. Subsequent calculations on the polycyclic aromatic and unsaturated hydrocarbons at the G4(MP2) level revealed excellent agreement for the IEFPCM-UFF and CPCM models against limited available experimental data. The IEFPCM-UFF-G4(MP2) and CPCM-G4(MP2) solvation energy calculation approaches are anticipated to give Kaw estimates within typical experimental ranges, each having general Kaw errors of less than 0.5 log10 units. When applied to other large organic compounds, the method should allow development of a broad and reliable Kaw database for multimedia environmental modeling efforts on various contaminants.

Products of biotransformation of polycyclic aromatic hydrocarbons in fishes of the Athabasca/Slave river system, Canada

Concentrations of products of biotransformation of polycyclic aromatic hydrocarbons (PBPAH) were measured in bile of five fishes of nutritional, cultural and ecological relevance from the Athabasca/Slave river system. Samples were collected in Alberta and the Northwest Territories, Canada, during three seasons. As a measure of concentrations of PBPAHs to which fishes are exposed and to gain information on the nature and extent of potential exposures of people or piscivorous wildlife, concentrations of biotransformation products of two- and three-ringed, four-ringed and five-ringed PAHs were measured using synchronous fluorescence spectroscopy. Spatial and seasonal differences were observed with greater concentrations of PBPAHs in samples of bile of fish collected from Fort McKay as well as greater concentrations of PBPAHs in bile of fish collected during summer compared to those collected in other seasons. Overall, PBPAHs were greater in fishes of lower trophic levels and fishes more closely associated with sediments. In particular, goldeye (Hiodon alosoides), consistently contained greater concentrations of all the PBPAHs studied.

Plasma proteome profiles of White Sucker (Catostomus commersonii) from the Athabasca River within the oil sands deposit

There are questions about the potential for oil sands related chemicals to enter the Athabasca River, whether from tailing ponds, atmospheric deposition, precipitation, or transport of mining dust, at concentrations sufficient to negatively impact the health of biota. We applied shotgun proteomics to generate protein profiles of mature male and female White Sucker (Catostomus commersonii) that were collected from various sites along the main stem of the Athabasca River in 2011 and 2012. On average, 399±131 (standard deviation) proteins were identified in fish plasma from each location in both years. Ingenuity Pathway Analysis software was used to determine the proteins' core functions and to compare the datasets by location, year, and sex. Principal component analysis (PCA) was used to determine if variation in the number of proteins related to a core function among all male and female individuals from both sampling years was affected by location. The core biological functions of plasma proteins that were common to both sampling years for males and females from each location were also estimated separately (based on Ingenuity's Knowledge Base). PCA revealed site-specific differences in the functional characteristics of the plasma proteome from white sucker sampled from downstream of oil sands extraction facilities compared with fish from upstream. Plasma proteins that were unique to fish downstream of oil sands extraction were related to lipid metabolism, small molecule biochemistry, vitamin and mineral metabolism, endocrine system disorders, skeletal and muscular development and function, neoplasia, carcinomas, and gastrointestinal disease.

Airborne Petcoke Dust is a Major Source of Polycyclic Aromatic Hydrocarbons in the Athabasca Oil Sands Region

Oil sands mining has been linked to increasing atmospheric deposition of polycyclic aromatic hydrocarbons (PAHs) in the Athabasca oil sands region (AOSR), but known sources cannot explain the quantity of PAHs in environmental samples. PAHs were measured in living Sphagnum moss (24 sites, n = 68), in sectioned peat cores (4 sites, n = 161), and snow (7 sites, n = 19) from ombrotrophic bogs in the AOSR. Prospective source samples were also analyzed, including petroleum coke (petcoke, from both delayed and fluid coking), fine tailings, oil sands ore, and naturally exposed bitumen. Average PAH concentrations in near-field moss (199 ng/g, n = 11) were significantly higher (p = 0.035) than in far-field moss (118 ng/g, n = 13), and increasing temporal trends were detected in three peat cores collected closest to industrial activity. A chemical mass-balance model estimated that delayed petcoke was the major source of PAHs to living moss, and among three peat core the contribution to PAHs from delayed petcoke increased over time, accounting for 45-95% of PAHs in contemporary layers. Petcoke was also estimated to be a major source of vanadium, nickel, and molybdenum. Scanning electron microscopy with energy-dispersive X-ray spectroscopy confirmed large petcoke particles (>10 μm) in snow at near-field sites. Petcoke dust has not previously been considered in environmental impact assessments of oil sands upgrading, and improved dust control from growing stockpiles may mitigate future risks.

Positive relationship detected between soil bioaccessible organic pollutants and antibiotic resistance genes at dairy farms in Nanjing, Eastern China

Co-contaminated soils by organic pollutants (OPs), antibiotics and antibiotic resistance genes (ARGs) have been becoming an emerging problem. However, it is unclear if an interaction exists between mixed pollutants and ARG abundance. Therefore, the potential relationship between OP contents and ARG and class 1 integron-integrase gene (intI1) abundance was investigated from seven dairy farms in Nanjing, Eastern China. Phenanthrene, pentachlorophenol, sulfadiazine, roxithromycin, associated ARG genes, and intI1 had the highest detection frequencies. Correlation analysis suggested a stronger positive relationship between the ARG abundance and the bioaccessible OP content than the total OP content. Additionally, the significant correlation between the bioaccessible mixed pollutant contents and ARG/intI1 abundance suggested a direct/indirect impact of the bioaccessible mixed pollutants on soil ARG dissemination. This study provided a preliminary understanding of the interaction between mixed pollutants and ARGs in co-contaminated soils.

Isotopic Evidence for Oil Sands Petroleum Coke in the Peace-Athabasca Delta

The continued growth of mining and upgrading activities in Canada's Athabasca oil sands (AOS) region has led to concerns about emissions of contaminants such as polycyclic aromatic hydrocarbons (PAHs). Whereas a recent increase in PAH emissions has been demonstrated within around 50 km of the main center of surface mining and upgrading operations, the exact nature of the predominant source(s) and the geographical extent of the deposition are still under debate. Here, we report a century-long source apportionment of PAHs using dual (δ(2)H, δ(13)C) compound-specific isotope analysis on phenanthrene deposited in a lake from the Athabasca sector of the Peace-Athabasca Delta situated ∼150 km downstream (north) of the main center of mining operations. The isotopic signatures in the core were compared to those of the main potential sources in this region (i.e., unprocessed AOS bitumen, upgrader residual coke, forest fires, coal, gasoline and diesel soot). A significant concurrent increase (∼55.0‰) in δ(2)H and decrease (∼1.5‰) in δ(13)C of phenanthrene over the last three decades pointed to an increasingly greater component of petcoke-derived PAHs. This study is the first to quantify long-range (i.e., >100 km) transport of a previously under-considered anthropogenic PAH source in the AOS region.

PAH Measurements in Air in the Athabasca Oil Sands Region

Polycyclic aromatic hydrocarbon (PAH) measurements were conducted by Wood Buffalo Environmental Association (WBEA) at four community ambient Air quality Monitoring Stations (AMS) in the Athabasca Oil Sands Region (AOSR) in Northeastern Alberta, Canada. The 2012 and 2013 mean concentrations of a subset of the 22 PAH species were 9.5, 8.4, 8.8, and 32 ng m(-3) at AMS 1 (Fort McKay), AMS 6 (residential Fort McMurray), AMS 7 (downtown Fort McMurray), and AMS 14 (Anzac), respectively. The average PAH concentrations in Fort McKay and Fort McMurray were in the range of rural and semirural areas, but peak values reflect an industrial emission influence. At these stations, PAHs were generally associated with NO, NO2, PM2.5, and SO2, indicating the emissions were from the combustion sources such as industrial stacks, vehicles, residential heating, and forest fires, whereas the PAH concentrations at AMS 14 (∼35 km south of Fort McMurray) were more characteristic of urban areas with a unique pattern: eight of the lower molecular weight PAHs exhibited strong seasonality with higher levels during the warmer months. Enthalpies calculated from Clausius-Clapeyron plots for these eight PAHs suggest that atmospheric emissions were dominated by temperature-dependent processes such as volatilization at warm temperatures. These findings point to the potential importance of localized water-air and/or surface-air transfer on observed PAH concentrations in air.

First results from the oil sands passive air monitoring network for polycyclic aromatic compounds

Results are reported from an ongoing passive air monitoring study for polycyclic aromatic compounds (PACs) in the Athabasca oil sands region in Alberta, Canada. Polyurethane foam (PUF) disk passive air samplers were deployed for consecutive 2-month periods from November 2010 to June 2012 at 17 sites. Samples were analyzed for polycyclic aromatic hydrocarbons (PAHs), alkylated PAHs, dibenzothiophene and its alkylated derivatives (DBTs). Relative to parent PAHs, alkylated PAHs and DBTs are enriched in bitumen and therefore considered to be petrogenic markers. Concentrations in air were in the range 0.03-210 ng/m(3), 0.15-230 ng/m(3) and 0.01-61 ng/m(3) for ∑PAHs, ∑alkylated PAHs and ΣDBTs, respectively. An exponential decline of the PAC concentrations in air with distance from mining areas and related petrogenic sources was observed. The most significant exponential declines were for the alkylated PAHs and DBTs and attributed to their association with mining-related emissions and near-source deposition, due to their lower volatility and greater association with depositing particles. Seasonal trends in concentrations in air for PACs were not observed for any of the compound classes. However, a forest fire episode during April to July 2011 resulted in greatly elevated PAH levels at all passive sampling locations. Alkylated PAHs and DBTs were not elevated during the forest fire period, supporting their association with petrogenic sources. Based on the results of this study, an "Athabasca PAC profile" is proposed as a potential source marker for the oil sands region. The profile is characterized by ∑PAHs/∑Alkylated PAHs = ∼0.2 and ∑PAHs/∑DBTs = ∼5.