Effects of nonnative species on the stability of riverine fish communities

Despite the increasing ubiquity of biological invasions worldwide, little is known about the scale-dependent effects of nonnative species on real-world ecological dynamics. Here, using an extensive time series dataset of riverine fish communities across different biogeographic regions of the world, we assessed the effects of nonnative species on the temporal variability and synchrony in abundance at different organizational levels (population, metapopulation, community and metacommunity) and spatial scales (stream reach and river basin). At the reach scale, we found that populations of nonnative species were more variable over time than native species, and that this effect scaled up to the community level – significantly destabilizing the dynamics of riverine fish communities. Nonnative species not only contributed to reduced community stability, but also increased variability of native populations. By contrast, we found no effect of nonnative species dominance on local interspecific synchrony among native species. At the basin scale, nonnative metapopulations were again more variable than the native ones. However, neither native metapopulations nor metacommunities showed differences in temporal variability or synchrony as nonnative species dominance increased basin-wide. This suggests a ‘dilution effect’ where the contribution to regional stability of local native populations from sites displaying low levels of invasion reduced the destabilizing effects of nonnative species. Overall, our results indicate that accounting for the destabilizing effect of nonnative species is critical to understanding native species persistence and community stability.     

Nitrogen addition does not reduce the role of spatial asynchrony in stabilising grassland communities

While nitrogen (N) amendment is known to affect the stability of ecological communities, whether this effect is scale‐dependent remains an open question. By conducting a field experiment in a temperate grassland, we found that both plant richness and temporal stability of community biomass increased with spatial scale, but N enrichment reduced richness and stability at the two scales considered. Reduced local‐scale stability under N enrichment arose from N‐induced reduction in population stability, which was partly attributable to the decline in local species richness, as well as reduction in asynchronous local population dynamics across species. Importantly, N enrichment did not alter spatial asynchrony among local communities, which provided similar spatial insurance effects at the larger scale, regardless of N enrichment levels. These results suggest that spatial variability among local communities, in addition to local diversity, may help stabilise ecosystems at larger spatial scales even in the face of anthropogenic environmental changes.     

Metacommunity Structure of Stream Insects across Three Hierarchical Spatial scales

A major challenge in community ecology is to understand the underlying factors driving metacommunity (i.e., a set of local communities connected through species dispersal) dynamics. However, little is known about the effects of varying spatial scale on the relative importance of environmental and spatial (i.e., dispersal related) factors in shaping metacommunities and on the relevance of different dispersal pathways. Using a hierarchy of insect metacommunities at three spatial scales (a small, within‐stream scale, intermediate, among‐stream scale, and large, among‐sub‐basin scale), we assessed whether the relative importance of environmental and spatial factors shaping metacommunity structure varies predictably across spatial scales, and tested how the importance of different dispersal routes vary across spatial scales. We also studied if different dispersal ability groups differ in the balance between environmental and spatial control. Variation partitioning showed that environmental factors relative to spatial factors were more important for community composition at the within‐stream scale. In contrast, spatial factors (i.e., eigenvectors from Moran's eigenvector maps) relative to environmental factors were more important at the among‐sub‐basin scale. These results indicate that environmental filtering is likely to be more important at the smallest scale with highest connectivity, while dispersal limitation seems to be more important at the largest scale with lowest connectivity. Community variation at the among‐stream and among‐sub‐basin scales were strongly explained by geographical and topographical distances, indicating that overland pathways might be the main dispersal route at the larger scales among more isolated sites. The relative effect of environmental and spatial factors on insect communities varied between low and high dispersal ability groups; this variation was inconsistent among three hierarchical scales. In sum, our study indicates that spatial scale, connectivity, and dispersal ability jointly shape stream metacommunities.     

Scaling‐up biodiversity‐ecosystem functioning research

A rich body of knowledge links biodiversity to ecosystem functioning (BEF), but it is primarily focused on small scales. We review the current theory and identify six expectations for scale dependence in the BEF relationship: (1) a nonlinear change in the slope of the BEF relationship with spatial scale; (2) a scale‐dependent relationship between ecosystem stability and spatial extent; (3) coexistence within and among sites will result in a positive BEF relationship at larger scales; (4) temporal autocorrelation in environmental variability affects species turnover and thus the change in BEF slope with scale; (5) connectivity in metacommunities generates nonlinear BEF and stability relationships by affecting population  synchrony at local and regional scales; (6) spatial scaling in food web structure and diversity will generate scale dependence in ecosystem functioning. We suggest directions for synthesis that combine approaches in metaecosystem and metacommunity ecology and integrate cross‐scale feedbacks. Tests of this theory may combine remote sensing with a generation of networked experiments that assess effects at multiple scales. We also show how anthropogenic land cover change may alter the scaling of the BEF relationship. New research on the role of scale in BEF will guide policy linking the goals of managing biodiversity and ecosystems.     

Dispersal promotes compensatory dynamics and stability in forced metacommunities

Understanding the factors that govern the stability of populations and communities has gained increasing importance as habitat fragmentation and environmental perturbations continue to escalate due to human activities. Dispersal is commonly viewed as essential to the maintenance of diversity in spatially subdivided communities, but few experiments have explored how dispersal interacts with the spatiotemporal components of environmental perturbations to determine community-level stability. We examined these processes using an experimental planktonic system composed of three competing species of zooplankton. We subjected zooplankton metacommunities to varying levels of dispersal and pH perturbations that varied in their degree of spatial synchrony. We show that dispersal can reverse the destabilizing effects of environmental forcing when perturbations are spatially asynchronous. Asynchrony in pH perturbations generated spatially and temporally varying species refugia that promoted source-sink dynamics and allowed prolonged persistence of zooplankton species that were otherwise extirpated in synchronously varying metacommunities. This, in turn, increased local species diversity, promoted compensatory population dynamics, and enhanced local community-level stability. Our results indicate that patterns of spatial covariation in environmental variability are critical to predicting the effects of dispersal on the dynamics and persistence of communities.     

The metacommunity concept: a framework for multi-scale community ecology

The metacommunity concept is an important way to think about linkages between different spatial scales in ecology. Here we review current understanding about this concept. We first investigate issues related to its definition as a set of local communities that are linked by dispersal of multiple potentially interacting species. We then identify four paradigms for metacommunities: the patch‐dynamic view, the species‐sorting view, the mass effects view and the neutral view, that each emphasizes different processes of potential importance in metacommunities. These have somewhat distinct intellectual histories and we discuss elements related to their potential future synthesis. We then use this framework to discuss why the concept is useful in modifying existing ecological thinking and illustrate this with a number of both theoretical and empirical examples. As ecologists strive to understand increasingly complex mechanisms and strive to work across multiple scales of spatio‐temporal organization, concepts like the metacommunity can provide important insights that frequently contrast with those that would be obtained with more conventional approaches based on local communities alone.     

Latitudinal patterns of forest ecosystem stability across spatial scales as affected by biodiversity and environmental heterogeneity

Our planet is facing a variety of serious threats from climate change that are unfolding unevenly across the globe. Uncovering the spatial patterns of ecosystem stability is important for predicting the responses of ecological processes and biodiversity patterns to climate change. However, the understanding of the latitudinal pattern of ecosystem stability across scales and of the underlying ecological drivers is still very limited. Accordingly, this study examines the latitudinal patterns of ecosystem stability at the local and regional spatial scale using a natural assembly of forest metacommunities that are distributed over a large temperate forest region, considering a range of potential environmental drivers. We found that the stability of regional communities (regional stability) and asynchronous dynamics among local communities (spatial asynchrony) both decreased with increasing latitude, whereas the stability of local communities (local stability) did not. We tested a series of hypotheses that potentially drive the spatial patterns of ecosystem stability, and found that although the ecological drivers of biodiversity, climatic history, resource conditions, climatic stability, and environmental heterogeneity varied with latitude, latitudinal patterns of ecosystem stability at multiple scales were affected by biodiversity and environmental heterogeneity. In particular, α diversity is positively associated with local stability, while β diversity is positively associated with spatial asynchrony, although both relationships are weak. Our study provides the first evidence that latitudinal patterns of the temporal stability of naturally assembled forest metacommunities across scales are driven by biodiversity and environmental heterogeneity. Our findings suggest that the preservation of plant biodiversity within and between forest communities and the maintenance of heterogeneous landscapes can be crucial to buffer forest ecosystems at higher latitudes from the faster and more intense negative impacts of climate change in the future.     

The sugarcane borer Diatraea saccharalis (Fabr.) (Lep., Crambidae) and its parasitoids: a synchrony approach to spatial and temporal dynamics

Abstract:  A data set on Diatraea saccharalis and its parasitoids, Cotesia flavipes and tachinid flies, was analysed at five spatial scales – sugarcane mill, region, intermediary, farm and zone – to determine the role of spatial scale in synchrony patterns, and on temporal population variability. To analyse synchrony patterns, only the three highest spatial scales were considered, but for temporal population variability, all spatial scales were adopted. The synchrony–distance relationship revealed complex spatial structures depending on both species and spatial scale. Temporal population variability [SD log( x  + 1)] levels were highest at the smallest spatial scales although, in the majority of the cases, temporal variability was inversely dependent on sample size. All the species studied, with a few exceptions, presented spatial synchrony independent of spatial scale. The tachinid flies exhibited stronger synchrony dynamics than D. saccharalis and C. flavipes in all spatial scales with the latter displaying the weakest synchrony levels, except when mill spatial scales were compared. In some cases spatial synchrony may at first decay and then increase with distance, but the presence of such patterns can change depending on the spatial scale adopted.     

Testing for temporal coherence across spatial extents: the roles of climate and local factors in regulating stream macroinvertebrate community dynamics

Temporal coherence or spatial synchrony refers to the tendency of population, community or ecosystem dynamics to behave similarly among locations through time as a result of spatially‐correlated environmental stochasticity (Moran effect), dispersal or trophic interactions. While terrestrial studies have treated synchrony mainly as a population‐level concept, the majority of freshwater studies have focused on community‐level patterns, particularly in lake planktonic communities. We used spatially and temporally hierarchical data on benthic stream invertebrates across six years, with three seasonal samples a year, in 11 boreal streams to assess temporal coherence at three spatial extents: 1) among regions (watersheds), 2) among streams within a region, and 3) among riffles within a stream, using the average of correlation coefficients for stream/riffle pairs across years. Our results revealed the primacy of strongly synchronized climatic factors (precipitation, air temperature) in inducing temporal coherence of macroinvertebrate assemblages across geographically distinct sites (i.e. Moran effect). Coherence tended to decrease with increasing spatial extent, but positive coherence was detected for most biological variables even at the largest extent (about 350 km). The generally high level of coherence reflected the strong seasonality of boreal freshwater communities. A hydrologically exceptional year enhanced the synchrony of biological variables, particularly total macroinvertebrate abundance. Regionally low precipitation in that year led to a substantial decrease in benthic densities across a broad spatial extent, followed by a rapid post‐drought recovery. Coherence at the among‐riffle (within‐stream) extent was lower than expected, implying that local‐scale habitat filters determine community dynamics at smaller spatial extents. Thus, temporal coherence of stream benthic communities appears to be controlled by partly different processes at different spatial scales.     

Hierarchical Dynamics of Ecological Communities: Do Scales of Space and Time Match?

Theory posits that community dynamics organize at distinct hierarchical scales of space and time, and that the spatial and temporal patterns at each scale are commensurate. Here we use time series modeling to investigate fluctuation frequencies of species groups within invertebrate metacommunities in 26 boreal lakes over a 20-year period, and variance partitioning analysis to study whether species groups with different fluctuation patterns show spatial signals that are commensurate with the scale-specific fluctuation patterns identified. We identified two groups of invertebrates representing hierarchically organized temporal dynamics: one species group showed temporal variability at decadal scales (slow patterns of change), whilst another group showed fluctuations at 3 to 5-year intervals (faster change). This pattern was consistently found across all lakes studied. A spatial signal was evident in the slow but not faster-changing species groups. As expected, the spatial signal for the slow-changing group coincided with broad-scale spatial patterns that could be explained with historical biogeography (ecoregion delineation, and dispersal limitation assessed through a dispersal trait analysis). In addition to spatial factors, the slow-changing groups correlated with environmental variables, supporting the conjecture that boreal lakes are undergoing environmental change. Taken together our results suggest that regionally distinct sets of taxa, separated by biogeographical boundaries, responded similarly to broad-scale environmental change. Not only does our approach allow testing theory about hierarchically structured space-time patterns; more generally, it allows assessing the relative role of the ability of communities to track environmental change and dispersal constraints limiting community structure and biodiversity at macroecological scales.     

Integrating over uncertainty in spatial scale of response within multispecies occupancy models yields more accurate assessments of community composition

Species abundance and community composition are affected not only by the local environment, but also by broader landscape and regional context. Yet, determining the spatial scales at which landscapes affect species remains a persistent challenge, hindering our ability to understand how environmental gradients shape communities. This problem is amplified by rare species and imperfect species detection. Here, we present a Bayesian framework that allows uncertainty surrounding the ‘true’ spatial scale of species’ responses (i.e. changes in presence/absence) to be integrated directly into a community hierarchical model. This scale‐selecting multispecies occupancy model (ssMSOM) estimates the scale of response, and shows high accuracy and correct levels of uncertainty in parameter estimates across a broad range of simulation conditions. An ssMSOM can be run in a matter of minutes, as opposed to the many hours required to run normal multispecies occupancy models at all queried spatial scales, and then conduct model selection – a problem that up to now has prohibited scale of response from being rigorously evaluated in an occupancy framework. Alternatives to the ssMSOM, such as GLM‐based approaches frequently fail to detect the correct spatial scale and magnitude of response, and are often falsely confident by favoring the incorrect parameter estimates, especially as species’ detection probabilities deviate from perfect. We further show how trait information can be leveraged to understand how individual species’ scales of response vary within communities. Integrating spatial scale selection directly into hierarchical community models provides a means of formally testing hypotheses regarding spatial scales of response, and more accurately determining the environmental drivers that shape communities.     

Population synchrony decreases with richness and increases with environmental fluctuations in an experimental metacommunity

Fluctuations of local but connected populations may show correlation or synchrony whenever they experience significant dispersal or correlated environmental biotic and abiotic variability. Synchrony may be an important variable in multispecies systems, but its nature and implications have not been explicitly examined. Because the number of locally coexisting species (richness) affects the population variability of community members, we manipulated richness under different regimes of environmental fluctuation (EF). We predicted that the temporal synchrony of populations in a species should decline with increasing richness of the metacommunity they live in. Additionally, we predicted that specialist species that are sensitive to a specific environmental factor would show higher synchronization when EF increases. We thus created experimental communities with varied richness, EF, and species specialization to examine the synchronizing effects of these factors on three aquatic invertebrate species. We created four levels of richness and three levels of EF by manipulating the salinity of the culture media. Monocultures exhibited higher population synchrony than metacommunities of 2–4 species. Furthermore, we found that species responded differently to EF treatments: high EF enhanced population synchrony for the specialist and intermediate species, but not for the generalist species. Our findings emphasize that the magnitude of EF and species richness both contribute to determine population synchrony, and importantly, our results suggest that biotic diversity may actually stabilize metacommunities by disrupting synchrony.     

A global synthesis of the effects of diversified farming systems on arthropod diversity within fields and across agricultural landscapes

Agricultural intensification is a leading cause of global biodiversity loss, which can reduce the provisioning of ecosystem services in managed ecosystems. Organic farming and plant diversification are farm management schemes that may mitigate potential ecological harm by increasing species richness and boosting related ecosystem services to agroecosystems. What remains unclear is the extent to which farm management schemes affect biodiversity components other than species richness, and whether impacts differ across spatial scales and landscape contexts. Using a global metadataset, we quantified the effects of organic farming and plant diversification on abundance, local diversity (communities within fields), and regional diversity (communities across fields) of arthropod pollinators, predators, herbivores, and detritivores. Both organic farming and higher in‐field plant diversity enhanced arthropod abundance, particularly for rare taxa. This resulted in increased richness but decreased evenness. While these responses were stronger at local relative to regional scales, richness and abundance increased at both scales, and richness on farms embedded in complex relative to simple landscapes. Overall, both organic farming and in‐field plant diversification exerted the strongest effects on pollinators and predators, suggesting these management schemes can facilitate ecosystem service providers without augmenting herbivore (pest) populations. Our results suggest that organic farming and plant diversification promote diverse arthropod metacommunities that may provide temporal and spatial stability of ecosystem service provisioning. Conserving diverse plant and arthropod communities in farming systems therefore requires sustainable practices that operate both within fields and across landscapes.     

Asynchrony among local communities stabilises ecosystem function of metacommunities

Temporal stability of ecosystem functioning increases the predictability and reliability of ecosystem services, and understanding the drivers of stability across spatial scales is important for land management and policy decisions. We used species‐level abundance data from 62 plant communities across five continents to assess mechanisms of temporal stability across spatial scales. We assessed how asynchrony (i.e. different units responding dissimilarly through time) of species and local communities stabilised metacommunity ecosystem function. Asynchrony of species increased stability of local communities, and asynchrony among local communities enhanced metacommunity stability by a wide range of magnitudes (1–315%); this range was positively correlated with the size of the metacommunity. Additionally, asynchronous responses among local communities were linked with species’ populations fluctuating asynchronously across space, perhaps stemming from physical and/or competitive differences among local communities. Accordingly, we suggest spatial heterogeneity should be a major focus for maintaining the stability of ecosystem services at larger spatial scales.     

Scale-dependent relationships between native richness, resource stability and exotic cover in dock fouling communities of Washington, USA

Aim In terrestrial plant communities, the relationship between native species diversity and exotic success is typically scale‐dependent. It is often proposed that within local neighbourhoods, high native diversity limits resources, thereby inhibiting exotic success. However, environmental variation that manifests over space or time can create positive correlations between native diversity and exotic success at larger scales. In marine habitats, there have been few multi‐scale surveys of this pattern, so it is unclear how diversity, resource limitation and the environment influence the success of exotic species in these systems. Location Washington, USA. Methods I analysed nested spatial and temporal surveys of fouling communities, which are assemblages of sessile marine invertebrates, to test whether the relationships between native richness, resource availability and exotic cover supported the diversity‐stability and diversity‐resistance theories, to test whether these relationships changed with spatio‐temporal scale, and to explore the temperature preferences of native and exotic fouling species. Results Survey data failed to support diversity‐stability theory: space availability actually increased with native richness at the local neighbourhood scale, and neither space availability nor variability decreased with native richness across larger spatio‐temporal scales. I did find support for diversity‐resistance theory, as richness negatively correlated with exotic cover in local neighbourhoods. Unexpectedly, this negative correlation disappeared at intermediate scales, but emerged again at the regional scale. This scale‐dependent pattern could be partially explained by contrasting water temperature preferences of native and exotic species. Main conclusions Within local neighbourhoods, native diversity may inhibit exotic abundance, but the mechanism is unlikely related to resource limitation. At the largest scale, correlations suggest that native richness is higher in cooler environments, whereas exotic richness is higher in warmer environments. This large‐scale pattern contrasts with the typical plant community pattern, and has important implications for coastal management in the face of global climate change.     

Fragmentation alters beta‐diversity patterns of habitat specialists within forest metacommunities

The relative importance of local, regional and historical factors in controlling the spatial patterns of plant species distribution is still poorly known and challenging for conservation ecology. We conducted an empirical study to link the spatial variation of species and environments among forest patches embedded in contrasted agricultural matrices. We compared how forest herb communities responded to spatial environmental gradients and past forest cover. We found low values of β‐diversity in both unfragmented and highly fragmented systems, independently from local and regional diversities. As fragmentation increased, the spatial structure of local plant communities was more complex and spatial effects explained an increasing proportion of β‐diversity, suggesting that the importance of dispersal limitations increased and played out at broad spatial scales. However, where spatio‐temporal isolation of forest patches was the highest, local species assemblages could not be explained, suggesting that the metacommunity functioning was disrupted. Where the historical continuity was high, local environmental characteristics explained a significant amount of species assemblages within metacommunities, suggesting habitat‐selection processes. Beta‐diversity and variations in presence–absence of species were also influenced by the intensity of landscape management, via the permeability of both forest edges and the matrix. This spatially‐explicit analysis of metacommunities revealed that forest fragmentation impacts beta‐diversity by altering not only the relative importance of deterministic and stochastic processes, but also the spatial scales at which they act. These results provide empirical support for the conservation of ancient forests and the maintenance of a high connectedness between fragments within agricultural landscapes.     

Frog community composition-environment relationships vary over time: Are snapshot studies of metacommunity dynamics useful?

There have been important advances in understanding the relative importance of environmental and spatial processes for the variation in species composition across a set of local communities linked by dispersal (i.e. metacommunities). However, community composition-environment relationships change over time, and the mechanisms shaping such temporal variation in metacommunities encompassing large environmental gradients remain poorly understood. If the ability of statistical models to predict community composition-environment relationships depends on the sampling year, snapshot metacommunity studies would have limited implications, both theoretical and applied. Here, we partitioned the variation in compositional data of frog communities and asked if the relative importance of environmental and spatial components change over years at broad spatial scales (hereafter, protected areas in coastal and inland regions) of southeastern Brazil. These regions have marked differences in environmental characteristics as well as the size and composition of their regional species pool. Our results showed that the factors explaining the temporal variability in community composition-environment relationships were congruent for the inland region, which is less productive and characterized by harsh environmental conditions. In contrast, the relative importance of environmental and spatial components changed over years in the coastal region, which has more productive environments and benign conditions. Although snapshot studies will continue to provide important information about metacommunity dynamics, researchers have to be better able to incorporate the temporal variation inherent in community composition-environment relationships, which may be especially important in productive environments.     

Effect of (a)synchronous light fluctuation on diversity,functional and structural stability of a marine phytoplankton metacommunity

Disentangling the mechanisms that maintain the stability of communities and ecosystem properties has become a major research focus in ecology in the face of anthropogenic environmental change. Dispersal plays a pivotal role in maintaining diversity in spatially subdivided communities, but only a few experiments have simultaneously investigated how dispersal and environmental fluctuation affect community dynamics and ecosystem stability. We performed an experimental study using marine phytoplankton species as model organisms to test these mechanisms in a metacommunity context. We established three levels of dispersal and exposed the phytoplankton to fluctuating light levels, where fluctuations were either spatially asynchronous or synchronous across patches of the metacommunity. Dispersal had no effect on diversity and ecosystem function (biomass), while light fluctuations affected both evenness and community biomass. The temporal variability of community biomass was reduced by fluctuating light and temporal beta diversity was influenced interactively by dispersal and fluctuation, whereas spatial variability in community biomass and beta diversity were barely affected by treatments. Along the establishing gradient of species richness and dominance, community biomass increased but temporal variability of biomass decreased, thus highest stability was associated with species-rich but highly uneven communities and less influenced by compensatory dynamics. In conclusion, both specific traits (dominance) and diversity (richness) affected the stability of metacommunities under fluctuating conditions.     

Community assembly is a race between immigration and adaptation: eco‐evolutionary interactions across spatial scales

Both ecological and evolutionary mechanisms have been proposed to describe how natural communities become assembled at both regional and biogeographical scales. Yet, these theories have largely been developed in isolation. Here, we unite these separate views and develop an integrated eco‐evolutionary framework of community assembly. We use a simulation approach to explore the factors determining the interplay between ecological and evolutionary mechanisms systematically across spatial scales. Our results suggest that the same set of ecological and evolutionary processes can determine community assembly at both regional and biogeographical scales. We find that the importance of evolution and community monopolization effects, defined as the eco‐evolutionary dynamics that occur when local adaptation of early established immigrants is fast enough to prevent the later immigration of better pre‐adapted species, are not restricted to adaptive radiations on remote islands. They occur at dispersal rates of up to ten individuals per generation, typical for many species at the scale of regional metacommunities. Dispersal capacity largely determines whether ecological species sorting or evolutionary monopolization structure metacommunity diversity and distribution patterns. However, other factors related to the spatial scale at which community assembly processes are acting, such as metacommunity size and the proportion of empty patches, also affect the relative importance of ecology versus evolution. We show that evolution often determines community assembly, and this conclusion is robust to a wide range of assumptions about spatial scale, mode of reproduction, and environmental structure. Moreover, we found that community monopolization effects occur even though species fully pre‐adapted to each habitat are abundant in the metacommunity, a scenario expected a priori to prevent any meaningful effect of evolution. Our results strongly support the idea that the same eco‐evolutionary processes underlie community assembly at regional and biogeographical scales.     

Species turnover in lake littorals: spatial and temporal variation of benthic macroinvertebrate diversity and community composition

Aims Despite wide consensus that ecological patterns and processes should be studied at multiple spatial scales, the temporal component of diversity variation has remained poorly examined. Specifically, rare species may exhibit patterns of diversity variation profoundly different from those of dominant taxa. Location Southern Finland. Methods We used multiplicative partitioning of true diversities (species richness, Shannon diversity) to identify the most important scale(s) of variation of benthic macroinvertebrate communities across several hierarchical scales, from individual samples to multiple littorals, lakes and years. We also assessed the among‐scale variability of benthic macroinvertebrate community composition by using measures of between‐ and within‐group distances at hierarchical scales. Results On average, a single benthic sample contained 23% of the total regional macroinvertebrate species pool. For both species richness and Shannon diversity, beta‐diversity was clearly the major component of regional diversity, with within‐littoral beta‐diversity (β₁) being the largest component of gamma‐diversity. The interannual component of total diversity was small, being almost negligible for Shannon index. Among‐sample (within‐littoral) diversity was related to variation of substratum heterogeneity at the same scale. By contrast, only a small proportion of rare taxa was found in an average benthic sample. Thus, dominant species among lakes and years were about the same, whereas rare species were mostly detected in a few benthic samples in one lake (or year). For rare species, the temporal component of diversity was more important than spatial turnover at most scales. Main conclusions While individual species occurrences and abundances, particularly those of rare taxa, may vary strongly through space and time, patterns of dominance in lake littoral benthic communities are highly predictable. Consequently, many rare species will be missed in temporally restricted samples of lake littorals. In comprehensive biodiversity surveys, interannual sampling of littoral macroinvertebrate communities is therefore needed.