Diatom–environmental relationships and nutrient transfer functions from contrasting shallow and deep limestone lakes in Ontario, Canada

Abstract

Diatom assemblages were analysed in the surface sediments of 44 alkaline lakes in south-western Ontario, Canada, and combined with a pre-existing 58 south-eastern Ontario lake set: (1) to determine if shallow, polymictic Ontario lakes contain different diatom assemblages from deeper, dimictic lakes, and if so, which environmental variables most influence assemblages; (2) to improve the existing transfer functions; (3) to construct and compare transfer functions separately for dimictic, deep lakes and for polymictic, shallow lakes. Polymictic and dimictic lakes covered a similar nutrient range (spring total phosphorus (TP)=4–54 μg/l, spring total nitrogen (TN)=200–927 μg/l; n=101) and spring pH levels (7.6–9; n=101). However, polymictic lakes were shallower (median mean depth = 2.9 m vs. 7.3 m in dimictic lakes). Benthic diatoms (average 60% relative abundance) dominated the polymictic lakes, whereas planktonic diatoms (average 60%) dominated dimictic lakes. A Canonical Correspondence Analyses with forward selection (p < 0.05, 999 Monte Carlo permutations) identified TP, alkalinity, watershed to volume ratios and lake depth as the most important measured environmental variables influencing diatom distribution in both polymictic and dimictic lakes. Additionally, pH was identified as an important variable in polymictic lakes, whereas TN was also forward selected in the dimictic lakes. Adding more lakes to the original southern Ontario calibration set improved the TN transfer function (r2 jack=0.42, root mean squared error of prediction (RMSEP)jack=0.11 [log μg TN/l]), although there was a high systematic error in the revised model (r2 residual = 0.48). However, the strongest TP model was derived from the polymictic lakes (r2 boot =0.44, RMSEPboot=0.20 [log μg TP/l]), which was the smallest lake set (n=30) with the lowest number of diatom species. The stronger TP model from the polymictic lakes may be partly due to the relatively low macrophyte cover in our polymictic lakes, which may lead to stronger benthic–pelagic coupling than in lakes with large macrophyte populations. Additionally, our study suggests that the Chrysophyceae cyst:diatom frustule ratio may be useful for indicating trends in TP levels of ≤35 μg/l in alkaline lakes that are dimictic, but is not necessarily indicative of trophic state changes in shallow, polymictic lakes. Our study demonstrates that it may be important to construct separate diatom-based nutrient transfer functions for polymictic and dimictic lakes.