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

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