Seasonality, species richness and poor dispersion mediate intraspecific trait variability in stonefly community responses along an elevational gradient

1. Across environmental gradients, some functional strategies are favoured over others resulting in differences in local species composition and distribution of any given functional trait. This generates among-Species Trait Variability (STV) across the gradient (reflecting species turnover), as well as spatial Intraspecific Trait Variability (ITV), which together contribute to Community Trait Variability (CTV). Understanding functional trait responses of freshwater biota is critical for elucidating the mechanisms that operate during community assembly and addressing community level responses to climate warming. 2. In freshwater insects, temperature and ecological conditions regulate growth rates. Here, we examined size-related trait variability components of short-dispersing stoneflies (Insecta: Plecoptera) within and among multispecies stonefly assemblages (henceforth referred to as communities), located along an elevational gradient in Sierra Nevada (Spain) to assess the components of trait-elevation matching (STV, ITV, CTV). 3. Relative to downstream reaches, headwaters hypothetically show greater internal physical habitat heterogeneity (large exposed rocks, tree roots, woody debris) and greater local species richness of cool-water adapted stoneflies. As a result, we expected CTV to increase with elevation, with a higher contribution of STV than ITV in upper reaches due to higher species richness and trait divergence among them. In contrast, the pattern of ITV should vary depending on whether the main driver across elevation is species richness (decrease due to niche partitioning) or ecological conditions (increase associated to physical habitat heterogeneity). 4. Eleven streams were sampled every 2 weeks over 1 year and five size-related traits were measured on 6,893 individuals belonging to 16 species. Overall, STV made a much larger contribution to CTV than ITV, which only became an important contributor in species-poor assemblages, especially in winter and autumn. Within each season, the trait-elevation matching was weak. Functional species-specific responses to the elevational gradient most often did not mirror the response of the entire communities, since highly variable (and even opposed or discontinuous) contributions to ITV were found across species. 5. Overall, ITV contributed little to trait-elevation matching, which suggests high vulnerability of short-dispersing stoneflies to habitat shifting associated to future climate warming due to limited phenotypic plasticity. Moreover, varying functional trait patterns across species emphasise species-specific responses to climate change as opposed to a unique whole community response. However, further trait-based macroecological studies across freshwater macroinvertebrates lineages are needed to compare patterns and establish generalisations.