Foliage Chemistry of Pinus baksiana in the Athabasca Oil Sands Region , Alberta , Canada

Industrial emissions in the Athabasca Oil Sands Region (AOSR), Alberta, Canada, have caused concerns about the effect of oil sands operations on the surrounding terrestrial environments, including jack pine (Pinus banksiana Lamb.) stands. We collected jack pine needles from 19 sites in the AOSR (13–128 km from main operations) for foliar chemical analyses to investigate the environmental impact on jack pine. Pine needles from three age classes, the current annual growth (CAG, 2011), one year and two year old pine needles, were collected. Samples were analyzed for total carbon (TC), nitrogen (TN), and sulfur (TS), inorganic S (SO4-S), base cations (Ca, Mg, Na), and other elements (B, Cu, Fe, Mn, P, Zn); CAG needles were also analyzed for their nitrogen and carbon isotopic compositions. Only TN, TS, Ca, B, Zn, and Fe contents showed weak but significant increases with proximity to the major oil sands operations. C and N isotopic compositions showed no trend with distance or TC and TN contents. Total S contents in CAG of pine foliage increased significantly with proximity to the main industrial operation while foliar inorganic S to organic S ratios (SO4-S/Sorg) ranged consistently between 0.13 and 0.32, indicating low to moderately high S loading. Hence, this study suggests some evidence of uptake of S emissions in close proximity to anthropogenic sources, although the reported values have not reached a level of environmental concern.


Introduction
The Athabasca Oil Sands Region (AOSR) in northeastern Alberta, Canada, is an enormous unconventional oil resource that constitutes the third largest oil reserve in the world.This oil resource has experienced an unprecedented expansion of mining projects (open pit mining) and in-situ operations (e.g., steam-assisted gravity drainage) over the last decades that has resulted in an increase in industrial emissions [1].Mining, bitumen upgrading, vehicle exhaust from heavy haulers operating in the open pit mines and other off-and on-road vehicles, are accompanied by significant industrial nitrogen (N) and sulfur (S) emissions.While earlier air quality monitoring in the AOSR focused on sulfur oxide (SO x ) emissions, current monitoring and modeling also includes nitrogen oxide (NO x ) emissions, base cations, and various other air quality parameters [1].Potential consequences of elevated S and N emissions include acid deposition and soil acidification [2,3], which are of particular concern in the boreal forest ecosystems of the AOSR and the adjacent province Saskatchewan [4][5][6][7][8][9][10].Jack pine forests growing on sandy soils with low base cation weathering rates were identified as most sensitive to acidic deposition in the AOSR [9,11].However, base cation deposition (Ca 2+ , Mg 2+ , Na + ) can mitigate impacts of acid deposition, and high base cation deposition from fugitive mine dust sources have been recently reported for the AOSR [12,13].Hence, the chemical analysis of jack pine foliage including main pollutants indicative of acid deposition as well as base cation deposition may serve as an effective indicator for the extent of environmental effects related to increased industrial emissions and deposition in the AOSR.
In addition to the determination of the contents of total nitrogen (TN) and carbon (TC) in pine needles, the nitrogen and carbon stable isotopic composition (δ 15 N, δ 13 C) may provide information about N and C sources and environmental or climatic stress conditions, and have been previously used in dendro-ecological investigations [14][15][16][17].Globally, δ 15 N values in atmospheric nitrogen deposition (e.g., NO 3 − , NH 4 + ) may range from −15‰ to +15‰ [18], but particulate matter and dry deposition are typically enriched in 15 N compared to δ 15 N values of dissolved N compounds in rainfall [19][20][21][22].Indeed, elevated δ 15 N values of up to 17.9‰ for δ 15 N-NO 3 and up to 20.1‰ for δ 15 N-NH 4 were found in stack fine particulate matter (PM 2.5 ) emissions and in atmospheric deposition in the vicinity of industrial developments in the AOSR [23], whereas nitrogen in bulk deposition at distant background sites revealed low δ 15 N values [22].Proemse et al. (2013) [22] showed that background atmospheric nitrate deposition (i.e., NO 3 -N deposition rates <0.1 kg•N•ha −1 ) were associated with an average δ 15 N-NO 3 value of −3.2‰ ± 1.5‰ during summer and 0.3‰ ± 1.7‰ during winter.δ 15 N-NH 4 in bulk deposition from background sites (NH 4 -N deposition rates <0.5 kg•N•ha −1 ) was characterized by average values of −3.6‰ ± 0.9‰ during summer and of −3.2‰ ± 2.2‰ during winter.δ 15 N values in bio-indicators such as lichens in the AOSR were also proven to be a useful qualitative indicator of N uptake within 3 km distance from industrial emissions [24].
Forest ecosystems are particularly sensitive to sulfur dioxide (SO 2 ), which is absorbed by trees through the stomata [25].Although both atmospheric and pedospheric sources have to be considered, sulfur concentrations in leaves tend to be more influenced by direct uptake of ambient SO 2 as opposed to sulfate (SO 4 ) uptake by roots [26].Total sulfur concentrations in pine needles were found to be a good indicator of SO 2 concentrations in ambient air [27].Conifers absorb gaseous SO 2 or take up SO 4 via roots, and both compounds are subsequently reduced to organic S forms (S org ), with the remainder being stored as sulfate [28,29].In either case, an over-supply of S may cause an accumulation of SO 4 in the plant tissue because the tree is unable to convert all SO 2 or SO 4 into S org .Hence, the ratio of total inorganic sulfur (SO 4 -S) to S org in pine needles has been suggested as an effective indicator for S loading and S stress [30].
Trace elements boron (B), copper (Cu), iron (Fe), manganese (Mn), and zinc (Zn) have been associated with the large number of sources and processes (natural and anthropogenic) in the AOSR and are being used as indicators of air pollutants [31].Copper and Zn are also on the US Environmental Protection Agency (EPA) list of priority pollutants.Phosphorus (P) is an essential nutrient for plant growth and P limitations in the boreal forest may be more widespread than previously thought and has been shown to increase with N additions.Thus, changes in P may affect forest health and in particular tree growth.
In 1998, a forest health monitoring program was established in the AOSR by the Wood Buffalo Environmental Association (WBEA), complementing the Terrestrial Environmental Effects Monitoring (TEEM) program, in order to identify air pollution effects on forest ecosystem health [31].As part of this framework, foliar samples were collected from jack pine trees (Pinus banksiana Lamb.) at 19 different sites in the AOSR in 2011, spanning distances from 13 to 128 km to the major oil sands developments.The objectives were: (1) to identify differences in foliar chemistry for the main pollutants (C, N, S), base cations (Ca, Mg, Na), trace metal pollutants (Cu, Fe, Mn, Zn), and other essential nutrients (B, P), and if differences exist, to determine whether the differences are attributed to oil sands development that potentially modified the ecosystem health of jack pine forests; (2) to investigate differences in the nitrogen and carbon isotopic signatures in pine needles, and if they provide information of N and C sources; and (3) to investigate S stress in jack pine trees using the SO 4 -S to S org ratio.

Site Selection and Foliage Sampling
The TEEM program established a rigorous sampling and site selection criteria protocol that identified ecological analogues with same soil, vegetation community, and orientation of sampling plot [32].In 2011, the sampling design was updated resulting in a total of 25 sites suitable for foliage sampling [31,33].Unfortunately, six sites burnt down during forest fires in summer 2011, leaving 19 sites for foliage collection for this study [34].
Foliar (needle) samples were collected in August/September 2011 before physiological changes due to autumn conditions occurred to allow for year-to-year comparison.Samples were collected at 19 sites (Figure 1) from five jack pine trees per plot resulting in five replicate samples for each of the three age classes and a total of 15 needle samples per sites.In addition, a duplicate set for each age class was obtained per site.Branches were clipped from the upper third of the canopy facing the regional emission sources and needles were identified for the current annual growth (CAG), year 1 (2010), and year 2 (2009) age classes.Distances between sampling sites and one of the major emission stacks nearest to site JP212 [35] ranged from 13 to 128 km.Prevailing winds in this region are from the West and Northwest [36].The most distant site JP201 (128 km) and site JP106 (74 km) were therefore located upwind of the operations, whereas JP316 (119 km) and JP213 (113 km) represent the most distant sampling sites downwind of the operations (Figure 1).The TEEM program established a rigorous sampling and site selection criteria protocol that identified ecological analogues with same soil, vegetation community, and orientation of sampling plot [32].In 2011, the sampling design was updated resulting in a total of 25 sites suitable for foliage sampling [31,33].Unfortunately, six sites burnt down during forest fires in summer 2011, leaving 19 sites for foliage collection for this study [34].
Foliar (needle) samples were collected in August/September 2011 before physiological changes due to autumn conditions occurred to allow for year-to-year comparison.Samples were collected at 19 sites (Figure 1) from five jack pine trees per plot resulting in five replicate samples for each of the three age classes and a total of 15 needle samples per sites.In addition, a duplicate set for each age class was obtained per site.Branches were clipped from the upper third of the canopy facing the regional emission sources and needles were identified for the current annual growth (CAG), year 1 (2010), and year 2 (2009) age classes.Distances between sampling sites and one of the major emission stacks nearest to site JP212 [35] ranged from 13 to 128 km.Prevailing winds in this region are from the West and Northwest [36].The most distant site JP201 (128 km) and site JP106 (74 km) were therefore located upwind of the operations, whereas JP316 (119 km) and JP213 (113 km) represent the most distant sampling sites downwind of the operations (Figure 1).

Chemical and Isotopic Analyses
Samples were dried at 70 °C for 24 h and ground in a zirconium ball mill.Total carbon (TC), and nitrogen (TN), and sulfur (TS) were analyzed by dry combustion with a precision of 2%; (TN and TC by Costech 4010 elemental analyzer (Costech Analytical Technologies Inc., Valencia, CA 91355 -USA), and TS by Eltra Helios S analyzer (ELTRA GmbH, Haan, Germany).Other elements were analyzed by inductively coupled plasma optical emission spectrometry (ICP-OES) following
Isotope abundance ratios are reported in the internationally accepted δ notation defined as where R is the 15 N/ 14 N or 13 C/ 12 C ratio of the sample and a standard respectively.The N isotopic composition (δ 15 N) of total nitrogen (TN) and the C isotopic composition (δ 13 C) of total carbon (TC) of the most recent (CAG, n = 109) foliage samples were determined using an elemental analyzer (EA, Costech 4010) interfaced with an isotope ratio mass spectrometer (IRMS, Finnigan MAT Delta + XL) at the Isotope Science Laboratory, University of Calgary, Canada.Samples loaded in tin cups were flash-combusted with an oxygen pulse in a column held at 1020 • C followed by purification on the generated measurement gases N 2 and CO 2 via gas chromatography (GC) and subsequent IRMS.The δ 15 N values are reported in per mille (‰) relative to atmospheric N 2 (AIR), and δ 13 C values are reported in per mille (‰) relative to Vienna Pee Dee Belemnite (V-PDB).
Accuracies and precisions of δ 15 N and δ 13 C values were ensured using the reference material USGS-40, USGS-41, NBS-22, USGS-32, and USGS-43.Precision and accuracy as 1 sigma of laboratory standards are ±0.2‰ for both δ 15 N and δ 13 C values.

Statistical Analysis and Predicted N and S Deposition Rates
Statistical analysis was performed using the free software R, and the presence of relationships was determined according to the Pearson correlation coefficient.Site averages and standard deviations were calculated for each element and needle year at each site, except for δ 15 N and δ 13 C values that were only determined for CAG needles.Correlation analysis was performed between individual elemental contents and needle age, but site averaged data was compared to distance and modeled N and S deposition rates to investigate the potential effect of industrial operations on a given site.Total N, S, and base cation (Ca 2+ , Mg 2+ , Na + ) deposition data for each site were available from predictions of the CALMET/CALPUFF dispersion model for the Lower Athabasca Region.This model is an ambient air quality model that has been applied to the AOSR in order to assess the magnitude and spatial distribution of atmospheric deposition [1,34].Measured deposition data in bulk deposition and throughfall (bulk deposition beneath the tree canopy) deposition has recently been published [13] but were not available for each sampling site, hence model predictions were used for comparison to foliar data.
Sodium concentrations did not correlate with needle age, averaging 9.6 ± 7.7 mg•kg −1 , 12.3 ± 11.0 mg•kg −1 , and 11.1 ± 8.4 mg•kg −1 for CAG, year 1, and year 2 needles, respectively.Site averaged Na contents did not correlate with distance for any of the needle ages.

N and C Isotopic Compositions of Pine Needles
The nitrogen isotopic composition (δ 15 N) of CAG pine needles ranged from −5.9‰ to 1.3‰ (Figure 4a, b).Plotting isotopic compositions versus the inverse concentration data reveals whether isotopic compositions can be explained by mixing of two or more sources.However, there was no trend of δ 15 N values with the inverse TN content (Figure 4a), nor with distance to the main industrial emission sources (Figure 4b) The carbon isotopic composition (δ 13 C) of CAG needles ranged from Magnesium concentrations decreased slightly with needle age (r = −0.27,p < 0.05), from 0.89 ± 0.16 g•kg −1 for CAG needles, 0.84 ± 0.19 g•kg −1 for year 1 needles, to 0.77 ± 0.19 g•kg −1 in year 2 needles.The highest site averaged CAG value was measured at JP212 (1.05 ± 0.17 g•kg −1 , 13 km), and the lowest at JP308 (0.71 ± 0.12 g•kg −1 , 75 km), but there was no significant correlation with distance.
Sodium concentrations did not correlate with needle age, averaging 9.6 ± 7.7 mg•kg −1 , 12.3 ± 11.0 mg•kg −1 , and 11.1 ± 8.4 mg•kg −1 for CAG, year 1, and year 2 needles, respectively.Site averaged Na contents did not correlate with distance for any of the needle ages.

N and C Isotopic Compositions of Pine Needles
The nitrogen isotopic composition (δ 15 N) of CAG pine needles ranged from −5.9‰ to 1.3‰ (Figure 4a,b).Plotting isotopic compositions versus the inverse concentration data reveals whether isotopic compositions can be explained by mixing of two or more sources.However, there was no trend of δ 15 N values with the inverse TN content (Figure 4a), nor with distance to the main industrial emission sources (Figure 4b) The carbon isotopic composition (δ 13 C) of CAG needles ranged from −29.8‰ to −26.7‰ (Figure 4c,d).Again, no significant trend was observed between δ 13 C values and the inverse TC content or distance (Figure 4c,d).Neither δ 15 N nor δ 13 C values of CAG needles correlated with CALPUFF predicted atmospheric deposition of N and S for the Lower Athabasca Region.

SO4-S/Sorg Ratios in Pine Needles
Site averaged ratios of SO4-S over Sorg in jack pine samples ranged from 0.13 ± 0.03 to 0.32 ± 0.17, 0.09 ± 0.06 to 0.30 ± 0.22, and 0.06 ± 0.03 to 0.18 ± 0.05 for CAG, year 1, and year 2 needles, respectively.The highest site averaged SO4-S/Sorg ratios were observed in CAG needles at JP303 (0.32 ± 0.17, 39 km).Site averaged SO4-S/Sorg ratios for each needle year did not correlate significantly with CALPUFF predicted S deposition rates for the Lower Athabasca Region.Site averaged ratios for all three needle ages did also not correlate significantly with distance (Figure 5).

SO 4 -S/S org Ratios in Pine Needles
Site averaged ratios of SO 4 -S over S org in jack pine samples ranged from 0.13 ± 0.03 to 0.32 ± 0.17, 0.09 ± 0.06 to 0.30 ± 0.22, and 0.06 ± 0.03 to 0.18 ± 0.05 for CAG, year 1, and year 2 needles, respectively.The highest site averaged SO 4 -S/S org ratios were observed in CAG needles at JP303 (0.32 ± 0.17, 39 km).Site averaged SO 4 -S/S org ratios for each needle year did not correlate significantly with CALPUFF predicted S deposition rates for the Lower Athabasca Region.Site averaged ratios for all three needle ages did also not correlate significantly with distance (Figure 5).CAG ratios of jack pine needles (* data from [24]) versus distance to the main oil sands operations.

Discussion
Carbon contents in jack pine needles in the vicinity of the Athabasca oil sands operations ranged from 492 to 574 g•kg −1 , and are similar to C contents reported for conifers [38].The lowest site averaged CAG C content was measured at JP310 (508 g•kg −1 ), which is 41 km south-west of the oil sands operations and Fort McMurray.The highest average CAG C content (548 g•kg −1 ) was measured at JP102 in direct proximity (16 km) of the main operators.The carbon isotopic composition in tree rings and foliage samples has been suggested as an indicator of the regional extent of pollution [16,17].However, neither C nor δ 13 C values correlated strongly with distance to the major industrial developments or with CALPUFF predicted S and N deposition rates.The carbon cycling in jack pine needles therefore appears to be unaffected by any potential environmental impacts of industrial emissions.
Total N contents increased from 10.2 g•kg −1 (CAG) at the most distant site (JP201, 128 km) towards higher values with decreasing distance to the oil sand operations (to 14.0 g•kg −1 at JP104, 14 km) and correlated with modeled N deposition data, suggesting that N contents in pine foliage is affected by oil sands emissions.δ 15 N values in CAG needles are similar to those previously reported for the AOSR (Table 1) and did not show a trend with distance or total N content (Figure 4a,b).There was no evidence for high δ 15 N values of foliar TN as observed for nitrate or ammonium of PM2.5 and atmospheric nitrate and ammonium deposition in close proximity to the oil sands operations as an indicator of industrial emissions [22,23], or as previously demonstrated for the lichen Evernia mesomorpha in close proximity (3 km) to one of the emission stacks [24].δ 15 N values of summer atmospheric nitrate and ammonium deposition in the AOSR at background deposition levels were −3.2‰ ± 1.5‰ and -3.6‰ ± 0.9‰, respectively [22], similar to δ 15 N values observed in total nitrogen of CAG needles.This suggests that these jack pine stands are N-limited and all available N is taken up by vegetation, allowing little or no room for N isotope fractionation [39].However, no δ 15 N values of gaseous NOx emissions (including vehicle emissions) are known for the AOSR, but may constitute a major source of nitrogen for conifer needles [40].[24]) versus distance to the main oil sands operations.

Discussion
Carbon contents in jack pine needles in the vicinity of the Athabasca oil sands operations ranged from 492 to 574 g•kg −1 , and are similar to C contents reported for conifers [38].The lowest site averaged CAG C content was measured at JP310 (508 g•kg −1 ), which is 41 km south-west of the oil sands operations and Fort McMurray.The highest average CAG C content (548 g•kg −1 ) was measured at JP102 in direct proximity (16 km) of the main operators.The carbon isotopic composition in tree rings and foliage samples has been suggested as an indicator of the regional extent of pollution [16,17].However, neither C nor δ 13 C values correlated strongly with distance to the major industrial developments or with CALPUFF predicted S and N deposition rates.The carbon cycling in jack pine needles therefore appears to be unaffected by any potential environmental impacts of industrial emissions.
Total N contents increased from 10.2 g•kg −1 (CAG) at the most distant site (JP201, 128 km) towards higher values with decreasing distance to the oil sand operations (to 14.0 g•kg −1 at JP104, 14 km) and correlated with modeled N deposition data, suggesting that N contents in pine foliage is affected by oil sands emissions.δ 15 N values in CAG needles are similar to those previously reported for the AOSR (Table 1) and did not show a trend with distance or total N content (Figure 4a,b).There was no evidence for high δ 15 N values of foliar TN as observed for nitrate or ammonium of PM 2.5 and atmospheric nitrate and ammonium deposition in close proximity to the oil sands operations as an indicator of industrial emissions [22,23], or as previously demonstrated for the lichen Evernia mesomorpha in close proximity (3 km) to one of the emission stacks [24].δ 15 N values of summer atmospheric nitrate and ammonium deposition in the AOSR at background deposition levels were −3.2‰ ± 1.5‰ and −3.6‰ ± 0.9‰, respectively [22], similar to δ 15 N values observed in total nitrogen of CAG needles.This suggests that these jack pine stands are N-limited and all available N is taken up by vegetation, allowing little or no room for N isotope fractionation [39].However, no δ 15 N values of gaseous NO x emissions (including vehicle emissions) are known for the AOSR, but may constitute a major source of nitrogen for conifer needles [40].The total sulfur contents in jack pine needles were similar to total S contents in spruce and pine needles observed in northern conifer stands in Finland [40].Younger leaves are more sensitive to atmospheric SO 2 concentrations compared to older ones [25], as reflected by slightly higher (SO 4 -S)/(S org ) ratios in younger pine needles (Table 2).Jaques and Legge [33] determined the following (SO 4 -S)/(S org ) ratios indicative for certain sulfur loading levels in the AOSR: <0.12: background levels; 0.12-0.19:low levels; 0.19-0.25:moderate levels; 0.25-0.35:moderately high levels; 0.35-0.45:high levels; and >0.45: very high levels of S loading.According to this classification, none of the sites are indicative of high levels or very high levels, but also none of the CAG needles were representative of what is considered background levels for S loading.Seven of the sites show moderately high levels, six sites moderate levels and the remainder are indicative of low levels of S loading (Table 2).Although none of the sites show (SO 4 -S)/(S org ) ratios of environmental concern, the increase of TS towards the emission sources demonstrates that industrial S is taken up by jack pine stands.(SO 4 -S)/(S org ) ratios reported for 2008 CAG jack pine needles [24] were in the same range as the 2011 CAG ratios found in this study (Figure 5), demonstrating no recent changes.We conclude that the investigated jack pine stands in the AOSR show currently no significant signs of sulfur stress resulting from elevated sulfur loadings and deposition rates as demonstrated by generally low to moderate (SO 4 -S)/(S org ) ratios, but that further monitoring is warranted.Calcium concentrations in the CAG jack pine foliage increased with decreasing distance to the oil sands operations for all needle ages (r = −0.59,r = −0.62,r = −0.64 for CAG, year 1, and year 2 needles; p < 0.05).There was a linear correlation to the CALPUFF predicted base cation throughfall deposition [34] (CAG: r = 0.66, p < 0.05); however, concentrations were still within the natural range of variation for Ca CAG jack pine foliage in western Canada (Table 3).Magnesium and Na concentrations did not vary with distance from the oil sands operation and were within the range of natural variation found in jack pine foliage in western Canada (Table 3).Watmough et al. [12] found that base cation deposition was dominated by Ca, likely from fugitive dust sources from mining sites, and suggested that effects of acid inputs (N and S) may have been offset by relatively high base cation deposition at the sites within 20 km of the oil sands operation.
Other elemental concentrations in jack pine foliage decreased following the trend P > Mn > Zn > Fe > B > Cu.Lynam et al. (2015) reported elevated crustal element concentrations of Mn, Fe, and Mg in rainfall for Air Monitoring Station (AMS) "6" at 30 km to the nearest upgrading facility emissions by oil sands operations in the AOSR [4,22,47].However, these studies typically lack sites within 10 km distance to the main oil sands operators, and we therefore suggest that future studies on foliage chemistry in the AOSR include sites in close proximity to the emitters.

Figure 1 .
Figure 1.Map of study site location in northeastern Alberta, Canada, showing the 19 sampling sites where foliage samples were collected in 2011.The grey shaded areas represent active mining operations as of 2013.Main operations are centered near JP212 and JP102.

Figure 1 .
Figure 1.Map of study site location in northeastern Alberta, Canada, showing the 19 sampling sites where foliage samples were collected in 2011.The grey shaded areas represent active mining operations as of 2013.Main operations are centered near JP212 and JP102.

Figure 2 .
Figure 2. Site averaged total S contents in CAG, year 1, and year 2 needles versus distance to the main oil sands operations.3.1.2.Base Cations (Ca, Mg, Na)

Figure 2 .
Figure 2. Site averaged total S contents in CAG, year 1, and year 2 needles versus distance to the main oil sands operations.

Figure 2 .
Figure 2. Site averaged total S contents in CAG, year 1, and year 2 needles versus distance to the main oil sands operations.3.1.2.Base Cations (Ca, Mg, Na)

Forests 2016, 7 , 312 7 of 14 − 29 .
8‰ to −26.7‰ (Figure4c, d).Again, no significant trend was observed between δ 13 C values and the inverse TC content or distance (Figure4c, d).Neither δ 15 N nor δ 13 C values of CAG needles correlated with CALPUFF predicted atmospheric deposition of N and S for the Lower Athabasca Region.

Figure 4 .
Figure 4. Nitrogen and carbon isotope results from the CAG pine needles.(a) δ 15 N values of total nitrogen (TN) versus the inverse TN (wt %) content; (b) site averaged δ 15 N values (black diamonds) and TN values (red diamonds) versus distance to the major emission sources in the AOSR; (c) δ 13 C values of total carbon (TC) versus the inverse TC (wt %) content; (d) site averaged δ 13 C values (black diamonds) and TC values (red diamonds) versus distance to the major emission sources in the AOSR.

Figure 4 .
Figure 4. Nitrogen and carbon isotope results from the CAG pine needles.(a) δ 15 N values of total nitrogen (TN) versus the inverse TN (wt %) content; (b) site averaged δ 15 N values (black diamonds) and TN values (red diamonds) versus distance to the major emission sources in the AOSR; (c) δ 13 C values of total carbon (TC) versus the inverse TC (wt %) content; (d) site averaged δ 13 C values (black diamonds) and TC values (red diamonds) versus distance to the major emission sources in the AOSR.

Figure 5 .
Figure 5. Site averaged SO4-S to Sorg ratios for CAG needles (2011, this study) in comparison to 2008CAG ratios of jack pine needles (* data from[24]) versus distance to the main oil sands operations.

Table 1 .
Average nitrogen contents and isotopic composition of jack pine needles in the AOSR.Site averaged SO 4 -S to S org ratios for CAG needles (2011, this study) in comparison to 2008 CAG ratios of jack pine needles (* data from

Table 1 .
Average nitrogen contents and isotopic composition of jack pine needles in the AOSR.

Table 2 .
(SO 4 -S)/(S org ) ratios (site averages and standard deviations) of jack pine needles in the AOSR.