A landscape approach to examining spatial patterns of limnological variables and long-term environmental change in a southern Canadian lake district

1. We explored patterns of limnological variables (physical, chemical and biological) with relation to landscape position (expressed as lake order) in 86 study lakes located on shield bedrock in south‐central Ontario, Canada. 2. Using anova s with lake order as the categorical variable, landscape position explained significant amounts of variation in major ion chemistry, physical and catchment characteristics, hypolimnetic oxygen, and community composition in algal (diatom, chrysophyte) and invertebrate (chironomid) assemblages preserved in surficial sediments. Several nutrient variables (TP, total phosphorus and TN, total nitrogen) and dissolved organic carbon did not have significant relationships with lake order. 3. The strongest relationships with lake order (as a fraction of variation explained in anova s) included silica concentrations (r² = 0.40) and SO₄ (r² = 0.29) concentrations, surface area (r² = 0.50) and hypolimnetic oxygen (r² = 0.29). 4. Bedrock geology (carbonate metasedimentary versus non‐carbonate bedrock) had strong influences on spatial gradients of pH and major ion chemistry. It was difficult to separate geological influences from spatial influences on limnological variables in this study, as drainage patterns in the region are highly influenced by surface features of underlying geological formations because of the very thin glacial till or exposed bedrock that exists in most catchments. 5. Patterns of limnological variables indicated that low‐order, headwater lakes had the lowest concentrations of major ions, and, from algal inferences of pH change, had been most susceptible to acidic deposition. High‐order, downstream lakes were larger and deeper, and had higher concentrations of hypolimnetic oxygen, indicating that optimal lake trout habitat was primarily located in high‐order lakes. 6. Variance partitioning analyses indicated that lake order as a metric of landscape position explained comparable portions of community variation in algal and invertebrate assemblages compared with geographic position (latitude, longitude) and Cartesian coordinate position (e.g. x, y, x², y², etc.) metrics. Lake order explained more community variation in chironomid assemblages compared with other landscape metrics, possibly because of the strong relationships between lake order and lake morphometry variables.