Wood Buffalo Environmental Association Research Articles

Field intercomparison of continuous ambient FRM and FEM NO2 instruments in the Athabasca Oil Sands Region, Alberta, Canada and the potential impact on ambient regulatory compliance

ABSTRACT The Canadian Federal Government promulgated new and lower NO2 Ambient Air Quality Standards (CAAQS) that went into effect in 2020 with additional decreases scheduled for 2025. The new hourly and annual NO2 CAAQS are 60 and 17 ppb, respectively, and the 2025 hourly and annual CAAQS are 42 and 12 ppb, respectively. The province of Alberta has also promulgated Ambient Air Quality Objectives (AAAQO) for NO2 currently set to 159 and 24 ppb on an hourly and annual basis, respectively. The Wood Buffalo Environmental Association (WBEA) in northeastern Alberta, Canada monitors NO2 at 21 community and industrial sites throughout the Athabasca Oil Sands Region (AOSR), for regulatory compliance using Thermo-Environmental (TEI) Model 42i Federal Reference Method (FRM) designated NO-NO2-NOx analyzers. The 42i measures NO directly via NO-O3 chemiluminescence, and NOx following the reduction of oxidized nitrogen to NO by a heated internal molybdenum converter. The difference between the NOx and NO channels is reported as NO2. This study presents the results of a three-year (2018–2021) WBEA comparison of four continuous NO2 analyzers: TEI 42i FRM; the API Model T500U cavity attenuated phase shift (CAPS) Federal Equivalent Method (FEM); a total reactive odd nitrogen analyzer (TEI Model 42i-Y); and a TEI 42i equipped with an external photolytic converter. The study showed that NO2 data from all analyzers were highly correlated and in general agreement, with r2 values (vs. the CAPS) ranging from 0.990–0.997 and slopes ranging from 0.933–0.992. Mean NO2 concentrations over the study period ranged from 7.2–7.5 ppb. Differences between the TEI 42i, TEI 42i-Y, and PhoNO, relative to the CAPS were all positive and highly significant (p < 0.0001), based upon nonparametric tests. The potential impact from the selection of different FRM/FEM measurement methods on current and future Canadian 2025 regulatory compliance in the region is evaluated. Implications: The study objective was to compare/evaluate different regulatory NO2 measurement techniques from a regional monitoring authority in a routine network operational context. Relatively small NO2 differences resulted in significant differences with respect to regulatory compliance triggers, particularly hourly standards based on daily extreme value statistics (e.g., 99th percentiles). For example, mean hourly NO2 △ differences ranged from 0.02–0.26 ppb over the study period but resulted in 2–3 ppb differences in the 3-year hourly CAAQS metrics. These differences could affect regulatory CAAQS and LARP compliance (management level) at monitoring sites observed during 2019 annual and 2020 hourly LARP trigger exceedances.

Suitability Assessment of Weather Networks for Wind Data Measurements in the Athabasca Oil Sands Area

The Athabasca Oil Sands Area (AOSA) in Alberta, Canada, is considered to have a high density of weather stations. Therefore, our objective was to determine an optimal network for the wind data measurement that could sufficiently represent the wind variability in the area. We used available historical data records of the weather stations in the three networks in AOSA, i.e., oil sands monitoring (OSM) water quantity program (WQP) and Wood Buffalo Environmental Association (WBEA) edge sites (ES) and meteorological towers (MT) of the air program. Both graphical and quantitative methods were implemented to find the correlations and similarities in the measurements between weather stations in each network. The graphical method (wind rose diagram) was found as a functional tool to understand the patterns of wind directions, but it was not appropriate to quantify and compare between wind speed data of weather stations. Therefore, we applied the quantitative method of the Pearson correlation coefficient (r) and absolute average error (AAE) in finding a relationship between the wind data of station pairs and the percentage of similarity (PS) method in quantifying the closeness/similarity. In the correlation analyses, we found weak to strong correlations in the wind data of OSM WQP (r = 0.04–0.69) and WBEA ES (r = 0.32–0.77), and a strong correlation (r = 0.33–0.86) in most of the station pairs of the WBEA MT network. In the case of AAE, we did not find any acceptable value within the standard operating procedure (SOP) threshold when logically combining the values of the u and v components together. In the similarity analysis, minor similarities were identified between the stations in the three networks. Hence, we presumed that all weather stations would be required to measure wind data in the AOSA.

Atmospheric sulfur and nitrogen deposition in the Athabasca oil sands region is correlated with foliar nutrient levels and soil chemical properties

The oil extraction industry and human activity in north eastern Alberta has been growing steadily since the 1960’s and is a source of air pollution. In the late 1990’s the Wood Buffalo Environmental Association was established to monitor air quality for both public and environmental health. A primary environmental concern was soil acidification caused by sulfur (S) and nitrogen (N) deposition. A network of forest health monitoring (FHM) sites was established in dry jack pine ecosystems to serve as an early indicator of negative impacts. A sampling campaign was executed in 2011 and this study examines soil properties and foliar nutrients in the context of measured and modeled acid deposition. Total N (TN), SO42−, pH, base cation to aluminum ratio (BC:Al), and base saturation (% BS) are reported for the organic layer (LFH) and 3 depths in the mineral soil, while foliar nutrients were analysed from current annual growth in jack pine needles. Atmospheric deposition of S, N, BC, and potential acid input (PAI) in the study area was recently provided by Edgerton et al. (2020) and soil and foliar chemistry was evaluated based on deposition estimates and measurements. Inverse distance weighting was used to examine spatial patterns and regression analysis was used to quantify relationships between variables. The results indicated that S deposition is spatially correlated with foliar SO42− concentration, and LFH SO42−, but not mineral topsoil (0–5 cm) SO42−. Nitrogen deposition was spatially correlated with foliar N concentration, but not LFH or topsoil TN indicating potential uptake by the foliage or rapid uptake by roots in the LFH layer. High BC deposition in the same areas with the highest potential acid inputs (PAI) did not correlate significantly with changes in soil pH. However, LFH pH was significantly related to dry NH3 deposition, which has not been reported previously and requires further investigation.

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.

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 PM2.5 and PM10–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 PM2.5 averaged 8.6 ± 11.8 μg m−3 (mean ± standard deviation), and PM10–2.5 averaged 8.5 ± 9.5 μg m−3. Wind regression analysis indicated that the oil sands production facilities were significant sources of PM2.5 mass and black carbon (BC), and that wildland fires were a significant source of the highest PM2.5 (>10 μg m−3) and BC events. A six-factor positive matrix factorization (PMF) model solution explained 95% of the measured PM2.5 and 78% of the measured ΣPAH. Five sources significantly contributed to PM2.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 PM10–2.5 mass and 86% of the ΣPAH. The predominant sources of PM10–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.

Using Pb isotope ratios of particulate matter and epiphytic lichens from the Athabasca Oil Sands Region in Alberta, Canada to quantify local, regional, and global Pb source contributions

Ambient air particulate matter (PM) was collected at the Wood Buffalo Environmental Association Bertha Ganter Fort McKay monitoring station in the Athabasca Oil Sand Region (AOSR) in Alberta, Canada from February 2010 to July 2011 as part of an air quality source assessment study. Daily 24-hour duration fine (PM2.5) and coarse (PM10–2.5) PM was collected using a sequential dichotomous sampler. 100 pairs of PM2.5 and PM10–2.5 were selected for lead (Pb) concentration and isotope analysis. Pb isotope and concentration results from 250 epiphytic lichen samples collected as far as 160 km from surface mining operations in 2008, 2011, and 2014 were analyzed to examine longer term spatial variations in Pb source contributions. A key finding was recognition of thorogenic 208Pb from eastern Asia in the springtime in the PM2.5 in 2010 and 2011. 206Pb/207Pb and 208Pb/207Pb isotope ratios were used in a three-component mixing model to quantify local, regional, and global Pb sources in the PM and lichen data sets. 47 ± 3% of the Pb in the PM2.5 at AMS-1 was attributed to sources from eastern Asia. Combined results from PM10–2.5 and PM2.5 indicate PM2.5 Pb contributions from eastern Asia (34%) exceed local AOSR sources of PM2.5 Pb (20%), western Canada sources of PM2.5 Pb (19%), and PM10–2.5 Pb from fugitive dust including oil sands (14%), tailings (10%), and haul roads (3%). The lichen analysis indicates regional sources contribute 46% of the Pb, local sources 32%, and global sources 22% over the 2008–2014 timeframe. Local sources dominate atmospheric Pb deposition to lichens at near field sites (0–30 km from mining operations) whereas regional Pb sources are prevalent at distal sites (30–160 km). The Pb isotope methodology successfully quantified trans-Pacific transport of Pb to the AOSR superimposed over the aerosol footprint of the world's largest concentration of bitumen mining and upgrading facilities.

The impact of the 2016 Fort McMurray Horse River Wildfire on ambient air pollution levels in the Athabasca Oil Sands Region, Alberta, Canada

An unprecedented wildfire impacted the northern Alberta city of Fort McMurray in May 2016 causing a mandatory city wide evacuation and the loss of 2,400 homes and commercial structures. A two-hectare wildfire was discovered on May 1, grew to ~157,000ha by May 5, and continued to burn an estimated ~590,000ha by June 13. A comprehensive air monitoring network operated by the Wood Buffalo Environmental Association (WBEA) in and around Fort McMurray provided essential health-related real-time air quality data to firefighters during the emergency, and provided a rare opportunity to elucidate the impact of gaseous and particulate matter emissions on near-field communities and regional air pollution concentrations. The WBEA network recorded 188 fire-related exceedances of 1-hr and 24-hr Alberta Ambient Air Quality Objectives. Two air monitoring sites within Fort McMurray recorded mean/maximum 1-hr PM2.5 concentrations of 291/5229μgm−3 (AMS-6) and 293/3259μgm−3 (AMS-7) during fire impact periods. High correlations (r2=0.83–0.97) between biomass combustion related gases (carbon monoxide (CO), non-methane hydrocarbons (NMHC), total hydrocarbons (THC), total reduced sulfur (TRS), ammonia) and PM2.5 were observed at the sites. Filter-based 24-hr integrated PM2.5 samples collected every 6 days showed maximum concentrations of 267μgm−3 (AMS-6) and 394μgm−3 (AMS-7). Normalized excess emission ratios relative to CO were 149.87±3.37μgm−3ppm−1 (PM2.5), 0.274±0.002ppmppm−1 (THC), 0.169±0.001ppmppm−1 (NMHC), 0.104±0.001ppmppm−1 (CH4), 0.694±0.007ppbppm−1 (TRS), 0.519±0.040ppbppm−1 (SO2), 0.412±0.045ppbppm−1 (NO), 1.968±0.053ppbppm−1 (NO2), and 2.337±0.077ppbppm−1 (NOX). A subset of PM2.5 filter samples was analyzed for trace elements, major ions, organic carbon, elemental carbon, and carbohydrates. Sample mass reconstruction and fire specific emission profiles are presented and discussed. Potential fire-related photometric ozone instrument positive interferences were observed and were positively correlated with NO and NMHC.

Source apportionment of ambient fine and coarse particulate matter at the Fort McKay community site, in the Athabasca Oil Sands Region, Alberta, Canada.

An ambient air particulate matter sampling study was conducted at the Wood Buffalo Environmental Association (WBEA) AMS-1 Fort McKay monitoring station in the Athabasca Oil Sand Region (AOSR) in Alberta, Canada from February 2010 to July 2011. Daily 24h integrated fine (PM2.5) and coarse (PM10-2.5) particulate matter was collected using a sequential dichotomous sampler. Over the duration of the study, 392 valid daily dichotomous PM2.5 and PM10-2.5 sample pairs were collected with concentrations of 6.8±12.9μgm-3 (mean±standard deviation) and 6.9±5.9μgm-3, respectively. A subset of 100 filter pairs was selected for element analysis by energy dispersive X-ray fluorescence and dynamic reaction cell inductively coupled plasma mass spectrometry. Application of the U.S. EPA positive matrix factorization (PMF) receptor model to the study data matrix resolved five PM2.5 sources explaining 96% of the mass including oil sands upgrading (32%), fugitive dust (26%), biomass combustion (25%), long-range Asian transport lead source (9%), and winter road salt (4%). An analysis of historical PM2.5 data at this site shows that the impact of smoke from wildland fires was particularly high during the summer of 2011. PMF resolved six PM10-2.5 sources explaining 99% of the mass including fugitive haul road dust (40%), fugitive oil sand (27%), a mixed source fugitive dust (16%), biomass combustion (12%), mobile source (3%), and a local copper factor (1%). Results support the conclusion of a previous epiphytic lichen biomonitor study that near-field atmospheric deposition in the AOSR is dominated by coarse fraction fugitive dust from bitumen mining and upgrading operations, and suggest that fugitive dust abatement strategies targeting the three major sources of PM10-2.5 (e.g., oil sand mining, haul roads, bulk material stockpiles) would significantly reduce near-field atmospheric deposition gradients in the AOSR and reduce ambient PM concentrations in the Fort McKay community.

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.

Total Gaseous Mercury Concentration Measurements at Fort McMurray, Alberta, Canada

Observations are described from total gaseous mercury (TGM) concentrations measured at the Wood Buffalo Environmental Association (WBEA) Fort McMurray—Patricia McInnes air quality monitoring station—from 21 October 2010 through 31 May 2013, inclusively. Fort McMurray is approximately 380 km north-northeast of Edmonton, Alberta, and approximately 30 km south of major Canadian oil sands developments. The average TGM concentration over the period of this study was 1.45 ± 0.18 ng∙m−3. Principal component analysis suggests that observed TGM concentrations are correlated with meteorological conditions including temperature, relative humidity, and solar radiation, and also ozone concentration. There is no significant correlation between ambient concentrations of TGM and anthropogenic pollutants, such as nitrogen oxides (NOX) and sulphur dioxide (SO2). Principal component analysis also shows that the highest TGM concentrations observed are a result of forest fire smoke near the monitoring station. Back trajectory analysis highlights the importance of long-range transport, indicating that unseasonably high TGM concentrations are generally associated with air from the southeast and west, while unseasonably low TGM concentrations are a result of arctic air moving over the monitoring station. In general, TGM concentration appears to be driven by diel and seasonal trends superimposed over a combination of long-range transport and regional surface-air flux of gaseous mercury.