Mapping Indicators of Toxicity for Polycyclic Aromatic Compounds in the Atmosphere of the Athabasca Oil Sands Region

Influence of oil extraction on concentration distributions, migration, secondary formation and carcinogenic risk of NPAHs and OPAHs in air and soil in an oilfield development area in China

Oil extraction leads to environmental pollution from the oilfields and dweller activities, however, knowledge of the concentration distributions, migration, secondary formation and toxicity of nitrated/oxygenated polycyclic aromatic hydrocarbons (N/OPAHs) in oilfield regions is limited. In this research, atmospheric and soil samples in 7 different location types in an important oil industrial base in China were gathered. The ΣNPAHs and ΣOPAHs in the air ranged from 0.05 to 2.47 ng/m3 and 0.14-22.72 ng/m3, respectively, and in soil ranged from 0.22 to 17.81 ng/g and 9.69-66.86 ng/g, respectively. Both NPAHs and OPAHs in the atmosphere exhibited higher concentrations during winter. The atmospheric NPAH concentrations decreased exponentially with distance from urban area especially in the summer, revealing the impact of vehicles on the air in the Yellow River Delta area. High NPAH and OPAH concentrations were found only in soil near oil extraction facilities, indicating that the impact of oil extraction is limited to the soil near the extraction facilities. The air-soil exchanges of N/OPAHs were assessed through fugacity fraction analysis, and NPAHs were in the equilibrium-deposition state and OPAHs were in the net-deposition state in the winter. Higher incremental lifetime cancer risk (ILCR) occurred at the urban, industrial, and oilfield sites in the atmospheric samples, and the soil samples had the largest ILCR values in the oilfield sites. However, ILCR values for both air and soil did not exceed the threshold of 10-6.

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.

Trace Metals in Global Air: First Results from the GAPS and GAPS Megacities Networks

Trace metals, as constituents of ambient air, can have impacts on human and environmental health. The Global Atmospheric Passive Sampling (GAPS) and GAPS Megacities (GAPS-MC) networks investigated trace metals in the air at 51 global locations by deploying polyurethane foam disk passive air samplers (PUF-PAS) for periods of 3-12 months. Aluminum and iron exhibited the highest concentrations in air (x̅ = 3400 and 4630 ng/m3, respectively), with notably elevated values at a rural site in Argentina thought to be impacted by resuspended soil. Urban sites had the highest levels of toxic Pb and Cd, with enrichment factors suggesting primarily anthropogenic influences. High levels of As at rural sites were also observed. Elevated trace metal concentrations in cities are associated with local emissions and higher PM2.5 and PM10 concentrations. Brake and tire wear-associated metals Sb, Cu, and Zn are significantly correlated and elevated at urban locations relative to those at background sites. These data demonstrate the versatility of PUF-PAS for measuring trace metals and other particle-associated pollutants in ambient air in a cost-effective and simple manner. The data presented here will serve as a global baseline for assessing future changes in ambient air associated with industrialization, urbanization, and population growth.

Characterization and Risk Assessment of PM2.5-Bound Polycyclic Aromatic Hydrocarbons and their Derivatives Emitted from a Typical Pesticide Factory in China

Polycyclic aromatic hydrocarbons (PAHs) and their derivatives have received extensive attention due to their negative effects on the environment and on human health. However, few studies have performed comprehensive assessments of PAHs emitted from pesticide factories. This study assessed the concentration, composition, and health risk of 52 PM_2.5-bound PAHs during the daytime and nighttime in the vicinity of a typical pesticide factory. The total concentration of 52 PAHs (Σ52PAHs) ranged from 53.04 to 663.55 ng/m³. No significant differences were observed between daytime and nighttime PAH concentrations. The average concentrations of twenty-two parent PAHs, seven alkylated PAHs, ten oxygenated PAHs, and twelve nitrated PAHs were 112.55 ± 89.69, 18.05 ± 13.76, 66.13 ± 54.79, and 3.90 ± 2.24 ng/m³, respectively. A higher proportion of high-molecular-weight (4-5 rings) PAHs than low-molecular-weight (2-3 rings) PAHs was observed. This was likely due to the high-temperature combustion of fuels. Analysis of diagnostic ratios indicated that the PAHs were likely derived from coal combustion and mixed sources. The total carcinogenic equivalent toxicity ranged from 15.93 to 181.27 ng/m³. The incremental lifetime cancer risk from inhalation, ingestion, and dermal contact with the PAHs was 2.33 × 10^-3 for men and 2.53 × 10^-3 for women, and the loss of life expectancy due to the PAHs was 11,915 min (about 0.023 year) for men and 12,952 min (about 0.025 year) for women. These results suggest that long-term exposure to PM_2.5 emissions from a pesticide factory has significant adverse effects on health. The study results support implementing the characterization of PAH emissions from pesticide factories and provides a scientific basis for optimizing the living environment around pesticide factories.

Novel eco-friendly methodology to determine polycyclic aromatic hydrocarbons in polyurethane foam for air monitoring: Application to spatial and temporal distribution survey

This study aimed at developing a new method for the extraction of polycyclic aromatic hydrocarbons (PAHs) in polyurethane foam (PUF). In the field of PAH monitoring, passive samplers using PUF disks are widely used. However, current extraction methods are time and solvent consuming. This new method employs 3 times a sixteenth of the PUF disk, with method detection limits (MDL) values below 5 and 13 ng/sampler for 3- and 4-rings PAHs, respectively. The use of only parts of the disk allows extraction by ultrasounds using exclusively 120 mL of ethanol, making it environmentally friendly. Ethanolic extracts are then purified and concentrated using microextraction by packed sorbent (MEPS) before GC-MS analyses. This method was applied for an environmental survey in a French urban area with an oceanic climate. Variations in PAH concentrations were observed depending on the site studied (urban, traffic, periurban and rural), as well as temporal variations.

Emissions and light absorption of PM2.5-bound nitrated aromatic compounds from on-road vehicle fleets

Vehicle emissions are an important source of nitrated aromatic compounds (NACs) in particulate size smaller 2.5 μm (PM_2.5), which adversely affect human health and biodiversity, especially in urban areas. In this study, filter-based PM_2.5 samples were collected during October 14-19, 2019, in a busy urban tunnel (approximately 35,000 vehicles per day) in south China to identify PM_2.5-bound NACs. Among them, 2,8-dinitrodibenzothiophene, 3-nitrodibenzofuran and 2-nitrodibenzothiophene were the most abundant nitrated polycyclic aromatic hydrocarbons (NPAHs), while 2-methyl-4-nitrophenol, 2,4-dinitrophenol, 3-methyl-4-nitrophenol and 4-nitrophenol were the most abundant nitrophenols (NPs). The observed mean fleet emission factors (EFs) of NPAHs and NPs were 2.2 ± 2.1 and 7.7 ± 4.1 μg km^-1, and were 2.9 ± 2.7 and 10.2 ± 5.4 μg km^-1 if excluding electric and liquefied petroleum gas vehicles, respectively. Regression analysis revealed that diesel vehicles (DVs) had NPAH-EFs (55.3 ± 5.3 μg km^-1) approximately 180 times higher than gasoline vehicles (GVs) (0.3 ± 0.2 μg km^-1), and NP-EFs (120.6 ± 25.8 μg km^-1) approximately 30 times higher than GVs (4.1 ± 0.2 μg km^-1), and thus 89% NPAH emissions and 56% NP emissions from the onroad fleets were contributed by DVs although DVs only accounted for 3.3% in the fleets. Methanol solution-based light absorption measurements demonstrated that the mean incremental light absorption for methanol-soluble brown carbon at 365 nm was 6.8 ± 2.2 Mm^-1, of which the 44 detected NACs only contributed about 1%. The mean EF of the 7 toxic NACs was approximately 3% that of the 16 priority PAHs; However, their benzo(a)pyrene toxic equivalence quotients (TEQ_BaP) could reach over 25% that of the PAHs. Moreover, 6-nitrochrysene mainly from DVs contributed 93% of the total TEQ_BaP of the NACs. This study demonstrated that enhancing DV emission control in urban areas could benefit the reduction of exposure to air toxins such as 6-nitrochrysene.

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

Toxicity screening of air extracts representing different source sectors in the Greater Toronto and Hamilton areas: In vitro oxidative stress, pro-inflammatory response, and toxicogenomic analysis

In the present study, the suitability and sensitivity of different in vitro toxicity endpoints were determined to evaluate and distinguish the specific contributions of polycyclic aromatic carbon (PAC) mixtures from various sites in Toronto (Canada), to pulmonary toxicity. Air samples were collected for two-month periods from April 2014 to March 2015 from one location, and from August 2016 to August 2017 from multiple locations reflecting different geographical areas in Toronto, and the Greater Toronto Area, with varying source emissions including background, traffic, urban, industrial and residential sites. Relative concentrations of PACs and their derivatives in these air samples were characterised. In vitro cytotoxicity, pro-inflammatory, and oxidative stress assays were employed to assess the acute pulmonary effects of urban-air-derived air pollutants. In addition, global transcriptional profiling was utilized to understand how these chemical mixtures exert their harmful effects. Lastly, the transcriptomic data and the chemical profiles for each site and season were used to relate the biological response back to individual constituents. Site-specific responses could not be derived; however, the Spring season was identified as the most responsive through benchmark concentration analysis. A combination of correlational analysis and principal component analysis revealed that nitrated and oxygenated polycyclic aromatic hydrocarbons (PAHs) drive the response at lower concentrations while specific PAHs drive the response at the highest concentration tested. Unsubstituted PAHs are the current targets for analysis as priority pollutants. The present study highlights the importance of by-products of PAH degradation in the assessment of risk. The study also demonstrates the usefulness of in vitro toxicity assays to derive meaningful data in support of risk assessment.

Insights into sources and occurrence of oxy- and nitro-PAHs in the alberta oil sands region using a network of passive air samplers

Mining-related activities in the Alberta Oil Sands Region (AOSR) are known to emit polycyclic aromatic hydrocarbons (PAHs) and related compounds to ambient air. This is a concern due to the toxicity of PAHs, including their transformation products such as nitrated (NPAHs) and oxygenated (OPAHs) PAHs. This is the first study that provided a more extensive outlook into the sources, occurrence in air, and spatial and seasonal patterns of NPAHs and OPAHs in the AOSR by using passive air sampling. A sampling campaign from 2013 to 2016 revealed concentrations of NPAHs that were much lower than those of OPAHs. The highest concentrations of NPAHs were concentrated in the region associated with extensive mining activities, with ∑NPAH concentrations ranging from 20 to 250 pg/m³. Within the oil sands (OS) mineable area, NPAHs associated with primary release appear more commonly, while NPAHs produced via oxidative transformation are predominant outside of this area. The concentrations of ∑OPAH ranged from 400 to 2400 pg/m³, with the highest air concentrations in the region located south of the main OS activity zone, with peak concentrations attributed to a 2016 forest fire event. Uptake of PAHs from ambient air and their subsequent conversion to generate OPAHs is believed to play an important role in wildfire emissions of OPAHs. The seasonal trend investigation was inconclusive, with NPAHs slightly higher during the winter, while OPAHs were slightly elevated during summer. A preliminary comparison of ambient concentrations of OPAHs and NPAHs in the AOSR to measurements in the Greater Toronto Area revealed a similar range of concentrations, but also a unique presence of certain NPAHs such as 4-nitrobiphenyl, 2-nitrodibenzothiophene, 2,8-dinitrodibenzothiophene and 6-nitrobenzo-(a)-pyrene. This indicates that AOSR might have its own NPAH profile - creating the need to better understand associated NPAH toxicity and propensity for long range transport.

The effects of plume episodes on PAC profiles in the athabasca oil sands region

Summer intensive air measurements of alkylated polycyclic aromatic compounds (Alk-PACs), nitrated polycyclic aromatic hydrocarbons (NPAHs), and oxygenated polycyclic aromatic hydrocarbons (OPAHs) was conducted during the summer of 2013 at an air monitoring site near the community of Fort McKay in the Athabasca oil sands region (AOSR). This study uses the ambient air measurements in conjunction with supplementary meteorological and air quality data from coordinated ground- and aircraft-based sampling over the same period to characterize diurnal variations and changes in the organic air pollutant profiles associated with the plume episodes. Principal component analysis showed a distinct PAC profile during plume episodes, driven mainly by higher fluorenone (FLO) and 9,10-anthraquinone (ANQ) concentrations. During the plume episodes (August 23-24), means of NPAHs and OPAHs concentrations were 120 and 2020 pg/m³, respectively, which were 2.7 and 2.5 times higher than those measured on the other days, while Alk-PACs did not reach maxima. The relative constancy of Alk-PACs during the plume episodes and baseline air quality periods likely reflects a continuous and broad emission of Alk-PACs from the oil sands mining activities. Only four OPAHs, including FLO, ANQ, benzo(a)fluorenone, and benzanthrone, exhibited higher average daytime than nighttime concentrations (p-value < 0.05). Categorizing air samples into clean and polluted conditions demonstrated that the polluted condition air samples were characterized by higher percent composition of alkylated fluorenes, FLO, MANQ, and photochemically-derived 1M4NN. A comparison of PAC profiles in air samples and oil sand ore samples suggests that the NPAHs were likely influenced by atmospheric formation while the OPAHs were impacted by a combination of primary sources and atmospheric formation. The strong correlations found between a number of NPAHs and OPAHs, and PM_2.5 and NO_x in this study could support the modelling of ambient air burdens of these compounds across the region.

Polycyclic aromatic compounds (PACs) in the Canadian environment: A review of sampling techniques, strategies and instrumentation

A wide variety of sampling techniques and strategies are needed to analyze polycyclic aromatic compounds (PACs) and interpret their distributions in various environmental media (i.e., air, water, snow, soils, sediments, peat and biological material). In this review, we provide a summary of commonly employed sampling methods and strategies, as well as a discussion of routine and innovative approaches used to quantify and characterize PACs in frequently targeted environmental samples, with specific examples and applications in Canadian investigations. The pros and cons of different analytical techniques, including gas chromatography - flame ionization detection (GC-FID), GC low-resolution mass spectrometry (GC-LRMS), high performance liquid chromatography (HPLC) with ultraviolet, fluorescence or MS detection, GC high-resolution MS (GC-HRMS) and compound-specific stable (δ¹³C, δ²H) and radiocarbon (Δ¹⁴C) isotope analysis are considered. Using as an example research carried out in Canada's Athabasca oil sands region (AOSR), where alkylated polycyclic aromatic hydrocarbons and sulfur-containing dibenzothiophenes are frequently targeted, the need to move beyond the standard list of sixteen EPA priority PAHs and for adoption of an AOSR bitumen PAC reference standard are highlighted.

Using tree cores to evaluate historic atmospheric concentrations and trends of polycyclic aromatic compounds in the Oil Sands region of Alberta, Canada

Tree cores and bark were sampled from jack pine trees at 18 sites in the Athabasca Oil Sands Region (AOSR) of Alberta, Canada, to investigate spatial and temporal trends of polycyclic aromatic compounds (PACs). Spatial trends were investigated in the bark samples, where ΣPAC concentrations ranged from 75 to 3615 ng/g. Highest concentrations were observed from trees within 40 km of the nearest mining or upgrading facility perimeter fence, in line with previous deposition studies in the AOSR. The sampled tree cores were separated into segments representing 5 years of growth/atmospheric collection by counting tree rings. A significant increase in PAC concentrations over the lifetime of the tree was observed at sites with the highest PAC concentrations, and the average % increase in concentration from 1970 to 2015 was in line with average % growth in bitumen extraction in the AOSR. Finally, the concentrations in the tree core segments representing collection from 2010 to 2015 were converted into an atmospheric PAC concentration using previously published wood-air partition coefficients. The calculated atmospheric concentrations were within the same range as concentrations reported from the passive atmospheric sampling network in this region. The importance of site location is highlighted, with forest edge sites providing an improved comparison for atmospheric exposure and deposition. This is the first study to use tree cores to calculate an atmospheric concentration of PACs, demonstrating the applicability of this methodology for providing historic atmospheric data.

Passive air sampling for semi-volatile organic chemicals

During passive air sampling, the amount of a chemical taken up in a sorbent from the air without the help of a pump is quantified and converted into an air concentration. In an equilibrium sampler, this conversion requires a thermodynamic parameter, the equilibrium sorption coefficient between gas-phase and sorbent. In a kinetic sampler, a time-averaged air concentration is obtained using a sampling rate, which is a kinetic parameter. Design requirements for kinetic and equilibrium sampling conflict with each other. The volatility of semi-volatile organic compounds (SVOCs) varies over five orders of magnitude, which implies that passive air samplers are inevitably kinetic samplers for less volatile SVOCs and equilibrium samplers for more volatile SVOCs. Therefore, most currently used passive sampler designs for SVOCs are a compromise that requires the consideration of both a thermodynamic and a kinetic parameter. Their quantitative interpretation depends on assumptions that are rarely fulfilled, and on input parameters, that are often only known with high uncertainty. Kinetic passive air sampling for SVOCs is also challenging because their typically very low atmospheric concentrations necessitate relatively high sampling rates that can only be achieved without the use of diffusive barriers. This in turn renders sampling rates dependent on wind conditions and therefore highly variable. Despite the overall high uncertainty arising from these challenges, passive air samplers for SVOCs have valuable roles to play in recording (i) spatial concentration variability at scales ranging from a few centimeters to tens of thousands of kilometers, (ii) long-term trends, (iii) air contamination in remote and inaccessible locations and (iv) indoor inhalation exposure. Going forward, thermal desorption of sorbents may lower the detection limits for some SVOCs to an extent that the use of diffusive barriers in the kinetic sampling of SVOCs becomes feasible, which is a prerequisite to decreasing the uncertainty of sampling rates. If the thermally stable sorbent additionally has a high sorptive capacity, it may be possible to design true kinetic samplers for most SVOCs. In the meantime, the passive air sampling community would benefit from being more transparent by rigorously quantifying and explicitly reporting uncertainty.

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.

Experimental Study of OH-Initiated Heterogeneous Oxidation of Organophosphate Flame Retardants: Kinetics, Mechanism, and Toxicity

The environmental risks and health impacts associated with particulate organophosphate flame retardants (OPFRs), which are ubiquitous in the global atmosphere, have not been adequately assessed due to the lack of data on the reaction kinetics, products, and toxicity associated with their atmospheric transformations. Here, the importance of such transformations for OPFRs are explored by investigating the reaction kinetics, degradation chemical mechanisms, and toxicological evolution of two OPFRs (2-ethylhexyl diphenyl phosphate (EHDP) and diphenyl phosphate (DPhP)) coated on (NH₄)₂SO₄ particles upon heterogeneous OH oxidation. The derived reaction rate constants for the heterogeneous loss of EHDP and DPhP are (1.12 ± 0.22) × 10^-12 and (2.33 ± 0.14) × 10^-12 cm³ molecules^-1 s^-1, respectively. Using recently developed real-time particle chemical composition measurements, particulate products from heterogeneous photooxidation and the associated degradation mechanisms for particulate OPFRs are reported for the first time. Subsequent cytotoxicity analysis of the unreacted and oxidized OPFR particles indicated that the overall particle cytotoxicity was reduced by up to 94% with heterogeneous photooxidation, likely due to a significantly lower cytotoxicity associated with the oxidized OPFR products relative to the parent OPFRs. The present work not only provides guidance for future field sampling for the detection of transformation products of OPFRs, but also strongly supports the ongoing risk assessment of these emerging chemicals and most critically, their products.

The Use of Geographic Information Systems for Spatial Ecological Risk Assessments: An Example from the Athabasca Oil Sands Area in Canada

There is an acknowledged need in ecotoxicology for methods that integrate spatial analyses in risk assessment. This has resulted in the emergence of landscape ecotoxicology, a subdiscipline of ecotoxicology. However, landscape ecotoxicology has yet to become common practice in risk assessment due to the underdevelopment of techniques and a lack of standardized methods. In the present study, we demonstrate how geographic information systems (GISs) can serve as a standardized platform to integrate data, assess spatial patterns of ecotoxicological data for multiple species, and assess relationships between chemical mixture exposures and effects on biota for landscape ecotoxicological risks assessment. We use data collected under the Joint Oil Sands Monitoring Program in the Athabasca Oil Sands Region in Alberta, Canada. This dataset is composed of concentrations of contaminants including metals and polycyclic aromatic compounds, and health endpoints measured in 1100 biological samples, including tree swallows, amphibians, gull and tern eggs, plants, and mammals. We present 3 examples using a GIS as a platform and geospatial analysis to: 1) integrate data and assess spatial patterns of contaminant exposure in the region, 2) assess spatial patterns of exposures to complex mixtures, and 3) examine patterns of exposures and responses across the landscape. We summarize the methods used in the present study into a workflow for ease of use. The GIS methods allow researchers to identify hot spots of contamination, use georeferenced monitoring data to derive quantitative exposure-response relationships, and assess complex exposures with more realism. Environ Toxicol Chem 2019;38:2797-2810. © 2019 SETAC.

Reanalysis of aerial deposition of metals and polycyclic aromatic compounds to snow in the Athabasca Oil Sands Region of Alberta Canada

Oil sands mining and bitumen upgrading activities in the Athabasca Oil Sands Region (AOSR) have been identified as sources of metals and polycyclic aromatic compounds (PAC) being deposited to the regional snowpack. We performed an independent reanalysis of publicly available AOSR snow pack data to: replicate previous results; to provide new insights into the spatial and temporal patterns of metal and PAC deposition; and, to determine whether certain metals or PACs were associated with specific oil sands mining or upgrading activities. Using PAC ratios, we use a K-means clustering approach to classify snowpack data into two combustion-dominated classes, and three classes associated with oil sands mining and bitumen upgrading. Snow samples dominated by "oil sands mine" emissions are consistent with a petrogenic source and exhibited low UNS ratios and high DBT ratios. Snow samples dominated by "petroleum coke" emissions had the highest BaP ratios, high DBT ratios, and were collected nearest the upgrader complexes. Metals data indicate snow samples dominated by oil sands mine emissions are consistent with an Athabasca Sands type composition. Those dominated by emissions from petroleum coke show enrichment of biophile metals V, Ni, and M. We conclude that previous studies have over-estimated environmental loadings of PACs, their spatial extent, and direction of their trend over time. These differences are attributed to the use of arithmetic rather than geometric spatial averaging, use of an arbitrary location (AR6) to determine the extent of metals and PAC deposition, and because previous studies neglected to account for metals and PACs being deposited from non-oil sands sources. Oil sands operators continue to reduce their emissions intensity, however there is an emerging consensus that mitigating fugitive emissions from petroleum coke stockpiles may represent the greatest opportunity to reduce environmental loadings of PACs in the AOSR.

Polycyclic aromatic compounds in urban air and associated inhalation cancer risks: A case study targeting distinct source sectors

Passive air sampling was conducted in Toronto and the Greater Toronto Area from 2016 to 2017 for 6 periods, in order to investigate ambient levels of polycyclic aromatic compounds (PACs) associated with different source types. The selected sampling sites (n = 8) cover geographical areas with varying source emissions including background, traffic, urban, industrial and residential sites. Passive air samples were analyzed for PACs which include PAHs, alkylated PAHs (alk-PAHs), dibenzothiophene and alkylated dibenzothiophenes (DBTs) and results for PAHs were used to calculate inhalation cancer risks using different approaches. The samples were also characterized for PAH derivatives including nitrated PAHs (NPAHs) and oxygenated PAHs (OPAHs). Concentrations of Σalk-PAHs and DBTs, which are known to be enriched in fossil fuels, as well as ΣNPAHs, were highest at a traffic site (MECP) located adjacent to the 18-lane Highway 401 that runs across Toronto. Except for an industrial site (HH/BU), PAC compositions were similar across the sampling sites with Σalk-PAHs being the most abundant class of PACs suggesting traffic emission was a major contributor to PACs in the atmosphere of Toronto. The industrial site exhibited a distinct chemical composition with ΣPAHs dominating over Σalk-PAHs and with elevated levels of fluoranthene, 9-nitroanthracene, and 9,10-anthraquinone, which likely reflects emissions from nearby industrial sources. MECP and HH/BU exhibited higher lifetime excess inhalation cancer risks indicating an association with traffic and industrial sources. The importance of the traffic sector as a source of PACs to ambient air is further supported by strong correlations of the ΣPAHs, Σalk-PAHs, DBTs, and ΣOPAHs with NO_x. This study highlights the importance of traffic as an emission source of PACs to urban air and the relevance of PAC classes other than just unsubstituted PAHs that are important but currently not included in air quality guidelines or for assessing inhalation cancer risks.

Seasonal variations of NPAHs and OPAHs in PM2.5 at heavily polluted urban and suburban sites in North China: Concentrations, molecular compositions, cancer risk assessments and sources

16 nitrated polycyclic aromatic hydrocarbons (NPAHs) and 5 oxygenated polycyclic aromatic hydrocarbons (OPAHs) in PM_2.5 at two locations in Northern China were analyzed by Gas Chromatography-Mass Spectrometry (GC-MS). Sampling was conducted at an urban site in Shandong University in Jinan (SDU) and a suburban site in Qixingtai in Jinan (QXT) in March, June, September and December in 2016. Overall, the concentrations of NPAHs and OPAHs were higher at SDU (1.88 and 9.49 ng/m³, respectively) than QXT (1.57 and 6.90 ng/m³, respectively), and the NPAHs and OPAHs concentrations were significantly higher during the winter than the other seasons at both sites. The incremental lifetime cancer risk (ILCR) values were lower than 10^-6 for all sites, seasons and age groups (ranging between 1.85E-08 and 2.56E-07), so there was no risk of carcinogenesis due to exposure to these pollutants. Total cancer risk at SDU was higher than QXT and NPAHs have the highest carcinogenic risk for adults aged from 30 to 70 years. The positive matrix factorization (PMF) results revealed that coal/biomass combustion, diesel vehicle emissions, gasoline vehicle emissions and secondary formation were the main sources of NPAHs and OPAHs at SDU and QXT. Coal/biomass combustion contributed more in spring, autumn and winter; diesel vehicle emission contributed the most in summer; secondary formation made greatest contributions in winter; the contributions of gasoline vehicle emission were similar in summer, autumn and winter. Diagnostic ratios clearly demonstrated that secondary formation is more active in winter than in other seasons, and the reactions of PAHs and OH radical were the dominant secondary formation pathway at both SDU and QXT. In addition, the potential source contribution function (PSCF) identified that the Beijing-Tianjin-Hebei region, Shandong province, Bohai Sea, Yellow Sea, Anhui province and Henan province were the main source regions of NPAHs and OPAHs in Jinan.

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.

Source apportionment of ambient fine and coarse particulate matter polycyclic aromatic hydrocarbons at the Bertha Ganter-Fort McKay community site in the Oil Sands Region of Alberta, Canada

A comprehensive filter-based particulate matter polycyclic aromatic hydrocarbon (PAH) source apportionment study was conducted at the Wood Buffalo Environmental Association Bertha Ganter-Fort McKay (BGFM) community monitoring station from 2014 to 2015 to quantify ambient concentrations and identify major sources. The BGFM station is located in close proximity to several surface oil sands production facilities and was previously found to be impacted by their air emissions. 24-hour integrated PM_2.5 and PM_10-2.5 samples were collected on a 1-in-3-day schedule yielding 108 complete organic/inorganic filter sets for source apportionment modeling. During the study period PM_2.5 averaged 8.6 ± 11.8 μg m^-3 (mean ± standard deviation), and PM_10-2.5 averaged 8.5 ± 9.5 μg m^-3. Wind regression analysis indicated that the oil sands production facilities were significant sources of PM_2.5 mass and black carbon (BC), and that wildland fires were a significant source of the highest PM_2.5 (>10 μg m^-3) and BC events. A six-factor positive matrix factorization (PMF) model solution explained 95% of the measured PM_2.5 and 78% of the measured ΣPAH. Five sources significantly contributed to PM_2.5 including: Biomass Combustion (3.57 μg m^-3; 40%); Fugitive Dust (1.86 μg m^-3; 28%); Upgrader Stack Emissions (1.44 μg m^-3; 21%); Petrogenic PAH (1.20 μg m^-3; 18%); and Transported Aerosol (0.43 μg m^-3 and 6%). However, the analysis indicated that only the pyrogenic PAH source factor significantly contributed (78%) to the measured ΣPAH. A five-factor PMF model dominated by fugitive dust sources explained 98% of PM_10-2.5 mass and 86% of the ΣPAH. The predominant sources of PM_10-2.5 mass were (i) Haul Road Dust (4.82 μg m^-3; 53%), (ii) Mixed Fugitive Dust (2.89 μg m^-3; 32%), (iii) Fugitive Oil Sand (0.88 μg m^-3; 10%), Mobile Sources (0.23 μg m^-3; 2%), and Organic Aerosol (0.06 μg m^-3; 1%). Only the Organic Aerosol source significantly contributed (86%) to the measured ΣPAH.

Flame retardants in urban air: A case study in Toronto targeting distinct source sectors

Based on distinct land-use categories, a sampling campaign was carried out at eight locations across Toronto and the Greater Toronto Area in 2016-2017. Source sectors' dependent patterns of atmospheric concentrations of 9 organophosphate esters (OPEs), 9 polybrominated diphenyl ethers (PBDEs) and 5 novel flame retardants (NFRs) showed dominance of OPEs and PBDEs at highly commercialised urban and traffic sites, while NFRs, were dominant at residential sites. Overall, average concentrations of Σ₉OPEs (1790 pg/m³) were two orders of magnitude higher than Σ₉PBDEs (9.17 pg/m³) and Σ₅NFRs (8.14 pg/m³). The atmospheric concentrations of given chemical classes also showed a general trend of lower levels in winter as compared to summer months. Statistically significant negative correlations between the natural logarithm of concentrations and inverse of temperature for some OPEs and PBDEs highlighted the role of volatilization from local sources at given sites as primarily influencing their atmospheric concentrations. Overall, this study adds to the current knowledge of urban settings as a major emitter of the chemicals of emerging concern and their replacements, as well as the ongoing problem of phased out PBDEs due to their presence in existing inventories of commercial/recycled products. It is recommended that long-term monitoring programs targeting flame retardants (FRs) include urban sites, which provide an early indicator of effectiveness of control measures of targeted FRs, while at the same time providing information on emission sources and trends of replacement FR chemicals.

Source apportionment of an epiphytic lichen biomonitor to elucidate the sources and spatial distribution of polycyclic aromatic hydrocarbons in the Athabasca Oil Sands Region, Alberta, Canada

The sources and spatial distribution of polycyclic aromatic hydrocarbons (PAHs) atmospheric deposition in the boreal forests surrounding bitumen production operations in the Athabasca Oil Sands Region (AOSR), Alberta, Canada were investigated as part of a 2014 passive in-situ bioindicator source apportionment study. Epiphytic lichen species Hypogymnia physodes samples (n = 127) were collected within a 150 km radius of the main surface oil sand production operations and analyzed for total sulfur, total nitrogen, forty-three elements, twenty-two PAHs, ten groups of C1-C2-alkyl PAHs and dibenzothiophenes (polycyclic aromatic compounds; PACs), five C1- and C2-alkyldibenzothiophenes, and retene. The ΣPAH + PAC in H. physodes ranged from 54 to 2778 ng g^-1 with a median concentration of 317 ng g^-1. Source apportionment modeling found an eight-factor solution that explained 99% of the measured ΣPAH + PAC lichen concentrations from four anthropogenic oil sands production sources (Petroleum Coke, Haul Road Dust, Stack Emissions, Raw Oil Sand), two local/regional sources (Biomass Combustion, Mobile Source), and two lichen biogeochemical factors. Petroleum Coke and Raw Oil Sand dust were identified as the major contributing sources of ΣPAH + PAC in the AOSR. These two sources accounted for 63% (43.2 μg g^-1) of ΣPAH + PAC deposition to the entire study domain. Of this overall 43.2 μg g^-1 contribution, approximately 90% (39.9 μg g^-1) ΣPAH + PAC was deposited within 25 km of the closest oil sand production facility. Regional sources (Biomass Combustion and Mobile Sources) accounted for 19% of ΣPAH + PAC deposition to the entire study domain, of which 46% was deposited near-field to oil sand production operations. Source identification was improved over a prior lichen-based study in the AOSR through incorporation of PAH and PAC analytes in addition to inorganic analytes.

Evaluating Computational and Structural Approaches to Predict Transformation Products of Polycyclic Aromatic Hydrocarbons

Polycyclic aromatic hydrocarbons (PAHs) undergo transformation reactions with atmospheric photochemical oxidants, such as hydroxyl radicals (OH•), nitrogen oxides (NOx), and ozone (O3). The most common PAH-transformation products (PAH-TPs) are nitrated, oxygenated, and hydroxylated PAHs (NPAHs, OPAHs, and OHPAHs, respectively), some of which are known to pose potential human health concerns. We sampled four theoretical approaches for predicting the location of reactive sites on PAHs (i.e., the carbon where atmospheric oxidants attack), and hence the chemoselectivity of the PAHs. All computed results are based on density functional theory (B3LYP/6-31G(d) optimized structures and energies). The four approaches are (1) Clar's prediction of aromatic resonance structures, (2) thermodynamic stability of all OHPAH adduct intermediates, (3) computed atomic charges (Natural Bond order, ChelpG, and Mulliken) at each carbon on the PAH, and (4) average local ionization energy (ALIE) at atom or bond sites. To evaluate the accuracy of these approaches, the predicted PAH-TPs were compared to published laboratory observations of major NPAH, OPAH, and OHPAH products in both gas and particle phases. We found that the Clar's resonance structures were able to predict the least stable rings on the PAHs but did not offer insights in terms of which individual carbon is most reactive. The OHPAH adduct thermodynamics and the ALIE approaches were the most accurate when compared to laboratory data, showing great potential for predicting the formation of previously unstudied PAH-TPs that are likely to form in the atmosphere.

Meltwater from snow contaminated by oil sands emissions is toxic to larval fish, but not spring river water

To assess the toxicity of winter-time atmospheric deposition in the oil sands mining area of Northern Alberta, embryo-larval fathead minnow (Pimephales promelas) were exposed to snowmelt samples. Snow was collected in 2011-2014 near (<7km) oil sands open pit mining operations in the Athabasca River watershed and at sites far from (>25km) oil sands mining. Snow was shipped frozen back to the laboratory, melted, and amended with essential ions prior to testing. Fertilized fathead minnow eggs were exposed (<24h post-fertilization to 7-16days post-hatch) to a range of 25%-100% snowmelt. Snow samples far from (25-277km away) surface mining operations and upgrading facilities did not affect larval fathead minnow survival at 100%. Snow samples from sites near surface mining and refining activities (<7km) showed reduced larval minnow survival. There was some variability in the potencies of snow year-to-year from 2011 to 2014, and there were increases in deformities in minnows exposed to snow from 1 site on the Steepbank River. Although exposure to snowmelt from sites near oil sands surface mining operations caused effects in larval fish, spring melt water from these same sites in late March-May of 2010, 2013 and 2014 showed no effects on larval survival when tested at 100%. Snow was analyzed for metals, total naphthenic acid concentrations, parent PAHs and alkylated PAHs. Naphthenic acid concentrations in snow were below those known to affect fish larvae. Concentrations of metals in ion-amended snow were below published water quality guideline concentrations. Compared to other sites, the snowmelt samples collected close to mining and upgrading activities had higher concentrations of PAHs and alkylated PAHs associated with airborne deposition of fugitive dusts from mining and coke piles, and in aerosols and particles from stack emissions.

CAPSULE

Snow collected close to oil sands surface mining sites is toxic to larval fathead minnows in the lab; however spring melt water samples from the same sites do not reduce larval fish survival.

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

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.

Measurement of polyurethane foam - air partition coefficients for semivolatile organic compounds as a function of temperature: Application to passive air sampler monitoring

Polyurethane foam - air partition coefficients (K_PUF-air) for 9 polycyclic aromatic hydrocarbons (PAHs), 10 alkyl-substituted PAHs, 4 organochlorine pesticides (OCPs) and dibenzothiophene were measured as a function of temperature over the range 5 °C-35 °C, using a generator column approach. Enthalpies of PUF-to-air transfer (ΔH_PUF-air, kJ/mol) were determined from the slopes of log K_PUF-air versus 1000/T (K), and have an average value of 81.2 ± 7.03 kJ/mol. The log K_PUF-air values at 22 °C ranged from 4.99 to 7.25. A relationship for log K_PUF-air versus log K_OA was shown to agree with a previous relationship based on only polychlorinated biphenyls (PCBs) and derived from long-term indoor uptake study experiments. The results also confirm that the existing K_OA-based model for predicting log K_PUF-air values is accurate. This new information is important in the derivation of uptake profiles and effective air sampling volumes for PUF disk samplers so that results can be reported in units of concentration in air.