Recent Warming, Rather than Industrial Emissions of Bioavailable Nutrients, Is the Dominant Driver of Lake Primary Production Shifts across the Athabasca Oil Sands Region

Freshwaters in the Athabasca Oil Sands Region (AOSR) are vulnerable to the atmospheric emissions and land disturbances caused by the local oil sands industry; however, they are also affected by climate change. Recent observations of increases in aquatic primary production near the main development area have prompted questions about the principal drivers of these limnological changes. Is the enhanced primary production due to deposition of nutrients (nitrogen and phosphorus) from local industry or from recent climatic changes? Here, we use downcore, spectrally-inferred chlorophyll-a (VRS-chla) profiles (including diagenetic products) from 23 limnologically-diverse lakes with undisturbed catchments to characterize the pattern of primary production increases in the AOSR. Our aim is to better understand the relative roles of the local oil sands industry versus climate change in driving aquatic primary production trends. Nutrient deposition maps, generated using geostatistical interpolations of spring-time snowpack measurements from a grid pattern across the AOSR, demonstrate patterns of elevated total phosphorus, total nitrogen, and bioavailable nitrogen deposition around the main area of industrial activity. However, this pattern is not observed for bioavailable phosphorus. Our paleolimnological findings demonstrate consistently greater VRS-chla concentrations compared to pre-oil sands development levels, regardless of morphological and limnological characteristics, landscape position, bioavailable nutrient deposition, and dibenzothiophene (DBT)-inferred industrial impacts. Furthermore, breakpoint analyses on VRS-chla concentrations across a gradient of DBT-inferred industrial impact show limited evidence of a contemporaneous change among lakes. Despite the contribution of bioavailable nitrogen to the landscape from industrial activities, we find no consistency in the spatial pattern and timing of VRS-chla shifts with an industrial fertilizing signal. Instead, significant positive correlations were observed between VRS-chla and annual and seasonal temperatures. Our findings suggest warmer air temperatures and likely decreased ice covers are important drivers of enhanced aquatic primary production across the AOSR.