Anne Lewerentz, Markus Hoffmann, Juliano Sarmento Cabral

• Investigating diversity gradients helps to understand biodiversity drivers and threats. However, one diversity gradient is rarely assessed, namely how plant species distribute along the depth gradient of lakes. Here, we provide the first comprehensive characterization of depth diversity gradient (DDG) of alpha, beta, and gamma species richness of submerged macrophytes across multiple lakes. We characterize the DDG for additive richness components (alpha, beta, gamma), assess environmental drivers, and address temporal change over recent years.Investigating diversity gradients helps to understand biodiversity drivers and threats. However, one diversity gradient is rarely assessed, namely how plant species distribute along the depth gradient of lakes. Here, we provide the first comprehensive characterization of depth diversity gradient (DDG) of alpha, beta, and gamma species richness of submerged macrophytes across multiple lakes. We characterize the DDG for additive richness components (alpha, beta, gamma), assess environmental drivers, and address temporal change over recent years. We take advantage of yet the largest dataset of macrophyte occurrence along lake depth (274 depth transects across 28 deep lakes) as well as of physiochemical measurements (12 deep lakes from 2006 to 2017 across Bavaria), provided publicly online by the Bavarian State Office for the Environment. We found a high variability in DDG shapes across the study lakes. The DDGs for alpha and gamma richness are predominantly hump-shaped, while beta richness shows a decreasing DDG. Generalized additive mixed-effect models indicate that the depth of the maximum richness (Dmax) is influenced by light quality, light quantity, and layering depth, whereas the respective maximum alpha richness within the depth gradient (Rmax) is significantly influenced by lake area only. Most observed DDGs seem generally stable over recent years. However, for single lakes we found significant linear trends for Rmax and Dmax going into different directions. The observed hump-shaped DDGs agree with three competing hypotheses: the mid-domain effect, the mean–disturbance hypothesis, and the mean–productivity hypothesis. The DDG amplitude seems driven by lake area (thus following known species–area relationships), whereas skewness depends on physiochemical factors, mainly water transparency and layering depth. Our results provide insights for conservation strategies and for mechanistic frameworks to disentangle competing explanatory hypotheses for the DDG.…