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

Emissions of organic substances with potential toxicity to humans and the environment are a major concern surrounding the rapid industrial development in the Athabasca oil sands region (AOSR). Although concentrations of polycyclic aromatic hydrocarbons (PAHs) in some environmental samples have been reported, a comprehensive picture of organic contaminant sources, pathways, and sinks within the AOSR has yet to be elucidated. We sought to use a dynamic multimedia environmental fate model to reconcile the emissions and residue levels reported for three representative PAHs in the AOSR. Data describing emissions to air compiled from two official sources result in simulated concentrations in air, soil, water, and foliage that tend to fall close to or below the minimum measured concentrations of phenanthrene, pyrene, and benzo(a)pyrene in the environment. Accounting for evaporative emissions (e.g., from tailings pond disposal) provides a more realistic representation of PAH distribution in the AOSR. Such indirect emissions to air were found to be a greater contributor of PAHs to the AOSR atmosphere relative to reported direct emissions to air. The indirect pathway transporting uncontrolled releases of PAHs to aquatic systems via the atmosphere may be as significant a contributor of PAHs to aquatic systems as other supply pathways. Emission density estimates for the three PAHs that account for tailings pond disposal are much closer to estimated global averages than estimates based on the available emissions datasets, which fall close to the global minima. Our results highlight the need for improved accounting of PAH emissions from oil sands operations, especially in light of continued expansion of these operations.