Biodiversity-Based Empirical Critical Loads of Nitrogen Deposition in the Athabasca Oil Sands Region

Anthropogenic nitrogen (N) emissions can have considerable effects on terrestrial ecosystems, with chronic N deposition leading to changes in plant species composition. The Athabasca Oil Sands Region (AOSR) represents a large point source of N emissions, which has prompted concern for surrounding habitats. The objective of this study was to determine the relative importance of N deposition as a driver of plant species community composition against bioclimatic and soil chemical variables. Further, we sought to identify community thresholds in plant species composition across a N deposition gradient. This assessment was performed for 46 Jack pine (Pinus banksiana Lamb.)-dominant forest sites surrounding the AOSR spanning Alberta and Saskatchewan. In total, 35 environmental variables were evaluated using redundancy analysis (RDA), followed by gradient forest analysis applied to plant species abundance data. Soil chemical variables accounted for just over 26% of the total explainable variation in the dataset, followed by bioclimatic variables (19%) and deposition variables (5%), but joint effects between variables also explained a significant portion of the total variation (p < 0.001). Total deposited nitrogen (TDN), and sulphur (TDS) along with bioclimatic and soil chemical variables, were identified as important variables in gradient forest analysis. A single, definitive threshold across TDN was identified at approximately 5.6 kg N ha−1 yr−1 (while a TDS threshold was found at 14.4 kg S ha−1 yr−1). The TDN threshold range was associated primarily with changepoints for several vascular species (Pyrola asarifolia, Pyrola chlorantha, Cornus canadensis, and Arctostaphylos uva-ursi) and bryophyte and lichen species (Pleurozium schreberi, Vulpicida pinastri, and Dicranum polysetum). These results suggest that across Jack pine-dominant forests surrounding the AOSR, the biodiversity-based empirical critical load of nutrient N is 5.6 kg N ha−1 yr−1.