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

Spatial variation and chemical reactivity of dusts from open-pit bitumen mining using trace elements in snow

The chemical reactivity of trace elements (TEs) in dusts from bitumen mining, upgrading and related industrial activities in the Athabasca Bituminous Sands region (ABS), Alberta, Canada, was evaluated using the acid-soluble fraction of snow. Samples were collected at 14 sites along the Athabasca River (AR) and its tributaries, and at 3 remote locations. Following metal-free, ultra-clean procedures for processing and analysis, samples were leached with nitric acid (pH < 1), filtered (<0.45 μm), and analyzed using ICP-MS. Insoluble particles (>0.45 μm) were examined using SEM-EDS. Along the river, acid-soluble concentrations of TEs varied by 6 orders of magnitude, from 1 mg/L (Al) to less than 1 ng/L (Tl). Conservative (Al, Y, La, Th) and mobile (Li, Be, Cs, Sr) lithophile elements, those enriched in bitumen (V, Ni, Mo), and potentially toxic chalcophile elements (As, Cd, Pb, Sb, Tl) showed considerable spatial variation. Normalizing the concentrations of TEs in samples collected near industry to the corresponding concentrations in snow from the reference site (UTK), resulted in enrichments of V and most of the lithophile elements. Dust reactivity, quantified as the ratio of acid-soluble to total concentrations, was less than 50% suggesting limited bioaccessibility. The large differences in behaviour between Cd and Pb versus Ni and V could be due to the occurrence of the former pair in carbonate or sulfide minerals, versus acid-insoluble petcoke particles for the latter couple. Spatial variations in the reactivity of TEs most likely reflect the range in diversity and chemical stability of dust particles, and variations in their abundance in primary source areas. The leaching conditions employed here are extreme (pH < 1) and intended to identify an upper limit of chemical reactivity, with far less dust dissolution expected when these dusts encounter natural waters of the area which range in pH from 4 to 8.

Co-PBK: a computational biomonitoring tool for assessing chronic internal exposure to chemicals and metabolites

Toxic chemicals are released into the environment through diverse human activities. An increasing number of chronic diseases are associated with ambient pollution, thus posing a threat to people. Given the high consumption of resources for human biomonitoring, this study proposed coupled physiologically-based kinetic (co-PBK) modeling matrices as a biomonitoring tool for simplifying chronic internal exposure estimates of environmental chemicals and their metabolites using naphthalene (NAP) and its metabolites (i.e., 1-OHN and 2-OHN) as simulation examples. According to the simulation of the steady-state mass among various organs/tissues via the co-PBK modeling matrices, fat had the highest potential bioaccumulation of NAP and its metabolites. With respect to body fluids, 1-OHN and 2-OHN tended to bioaccumulate more in the bile than in the urine. According to the sensitivity analysis, the calculated sensitivity factors for the first-order kinetics-based rate constants imply that due to the biotransformation process, target organs/tissues (e.g., liver and kidneys) would be continuously exposed to more NAP metabolites under chronic exposure. Meanwhile, 1-OHN may be more stably transported to the urine than 2-OHN for further human biomonitoring during long-term internal exposure. According to the case study of simulating population chronic exposure to NAP in Shenzhen, the co-PBK modeling estimated the population exposure to NAP with an intake rate of 8.77 × 10-2 mg d-1 and the aggregated urinary concentration of NAP metabolites of 2.60 μg L-1. Furthermore, the accuracy of the urinary levels between the real-world data and the values simulated by the co-PBK modeling was assessed and the root-mean-square error of c1-OHN,urine was found to be lower than that of c2-OHN,urine.

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.

Incorporating Industrial and Climatic Covariates into Analyses of Fish Health Indicators Measured in a Stream in Canada’s Oil Sands Region

Industrial and other human activities in Canada’s oil sands region (OSR) influence the environment. However, these impacts can be challenging to separate from natural stresses in flowing waters by comparing upstream reference sites to downstream exposure locations. For example, health indicators of lake chub (Couesius plumbeus) compared between locations in the Ells River (Upper and Lower) in 2013 to 2015 and 2018 demonstrated statistical differences. To further examine the potential sources of variation in fish, we also analyzed data at sites over time. When fish captured in 2018 were compared to pooled reference years (2013–2015), results indicated multiple differences in fish, but most of the differences disappeared when environmental covariates were included in the Elastic Net (EN) regularized regression models. However, when industrial covariates were included separately in the EN, the large differences in 2018 also disappeared, also suggesting the potential influence of these covariables on the health of fish. Further ENs incorporating both environmental and industrial covariates along with other variables which may describe industrial and natural influences, such as spring or summer precipitation and summer wind speeds and distance-based penalty factors, also support some of the suspected and potential mechanisms of impact. Further exploratory analyses simulating changes from zero and the mean (industrial) activity levels using the regression equations respectively suggest effects exceeding established critical effect sizes (CES) for fish measurements may already be present or effects may occur with small future changes in some industrial activities. Additional simulations also suggest that changing regional hydrological and thermal regimes in the future may also cause changes in fish measurements exceeding the CESs. The results of this study suggest the wide applicability of the approach for monitoring the health of fish in the OSR and beyond. The results also suggest follow-up work required to further evaluate the veracity of the suggested relationships identified in this analysis.