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.     

The assembly of ecological communities inferred from taxonomic and functional composition

Among-site variation in metacommunities (beta diversity) is typically correlated with the distance separating the sites (spatial lag). This distance decay in similarity pattern has been linked to both niche-based and dispersal-based community assembly hypotheses. Here we show that beta diversity patterns in community composition, when supplemented with functional-trait information, can be used to diagnose assembly processes. First, using simulated data, we show how the relationship between distance decay patterns in taxonomic and functional measures of community composition can be used to predict the influence of a given trait on community assembly. We then use the patterns generated by the simulation as a template to show that the sorting of benthic macroinvertebrate metacommunities in headwater streams is likely influenced by different sets of functional traits at regional and local scales. We suggest that functional-trait databases and spatially referenced taxonomic surveys can be used to predict the spatial scales at which different aspects of interspecific functional variation are involved in niche-based community assembly while accounting for the influence of dispersal-based community assembly processes.     

Biodiversity–productivity relationship in ponds: Community and metacommunity patterns along time and environmental gradients

Primary production correlates with diversity in various ways. These patterns may result from the interaction of various mechanisms related to the environmental context and the spatial and temporal scale of analysis. However, empirical evidence on diversity‐productivity patterns typically considers single temporal and spatial scales, and does not include the effect of environmental variables. In a metacommunity of macrophytes in ephemeral ponds, we analysed the diversity‐productivity relationship patterns in the field, the importance of the environmental variables of pond size and heterogeneity on such relationship, and the variation of these patterns at local (community level) and landscape scales (metacommunity level) across 52 ponds on twelve occasions, over five years (2005–2009). Combining all sampling dates, there were 377 ponds and 1954 sample‐unit observations. Vegetation biomass was used as a proxy for productivity, and biodiversity was represented by species richness, evenness, and their interaction. Environmental variables comprised pond area, depth and internal heterogeneity. Productivity and species richness were not directly related at the metacommunity level, and were positively related at the community level. Taking environmental variables into account revealed positive species richness‐productivity relationships at the metacommunity level and positive quadratic relationships at the community level. Productivity showed both positive and negative linear and nonlinear relationships with the size and heterogeneity of ponds. We found a weak relationship between productivity and evenness. The identity of variables associated with productivity changed between spatial scales and through time. The pattern of relationships between productivity and diversity depends on spatial scale and environmental context, and changes idiosyncratically through time within the same ecosystem. Thus, the diversity‐productivity relationship is not only a property of the study system, but also a consequence of environmental variations and the temporal and spatial scale of analysis.     

Using functional diversity patterns to explore metacommunity dynamics: a framework for understanding local and regional influences on community structure

The metacommunity concept, describing how local and regional scale processes interact to structure communities, has been successfully applied to patterns of taxonomic diversity. Functional diversity has proved useful for understanding local scale processes, but has less often been applied to understanding regional scale processes. Here, we explore functional diversity patterns within a metacommunity context to help elucidate how local and regional scale processes influence community assembly. We detail how each of the four metacommunity perspectives (species sorting, mass effects, patch dynamics, neutral) predict different patterns of functional beta‐ and alpha‐diversity and spatial structure along two key gradients: dispersal limitation and environmental conditions. We then apply this conceptual model to a case study from alpine tundra plant communities. We sampled species composition in 17 ‘sky islands’ of alpine tundra in the Colorado Rocky Mountains, USA that differed in geographic isolation and area (key factors related to dispersal limitation) and temperature and elevation (key environmental factors). We quantified functional diversity in each site based on specific leaf area, leaf area, stomatal conductance, plant height and chlorophyll content. We found that colder high elevation sites were functionally more similar to each other (decreased functional beta‐diversity) and had lower functional alpha‐diversity. Geographic isolation and area did not influence functional beta‐ or alpha‐diversity. These results suggest a strong role for environmental conditions structuring alpine plant communities, patterns consistent with the species sorting metacommunity perspective. Incorporating functional diversity into metacommunity theory can help elucidate how local and regional factors structure communities and provide a framework for observationally examining the role of metacommunity dynamics in systems where experimental approaches are less tractable.     

Network isolation and local diversity in neutral metacommunities

Biologists seek an understanding of the biological and environmental factors determining local community diversity. Recent advances in metacommunity ecology, and neutral theory in particular, highlight the importance of dispersal processes interacting with the spatial structure of a landscape for generating spatial patterns and maintaining biodiversity. The relative spatial isolation of a community is traditionally thought to have a large influence on local diversity. However, isolation remains an elusive concept to quantify, particularly in metacommunities with complex spatial structure. We represent the metacommunity as a network of local communities, and use network centrality measures to quantify the isolation of a local community. Using spatially explicit neutral theory, we examine how node position predicts variation in alpha diversity across a metacommunity. We find that diversity increases with node centrality in the network, but only when centrality is measured on a given scale in the network that widens with increasing dispersal rates and narrows with increasing evolutionary rates. More generally, complex biodiversity patterns form only when the underlying geography has structure on this critical scale. This provides a framework for understanding the influence of spatial geographic structure on global biodiversity patterns.     

A conceptual framework for the spatial analysis of functional trait diversity

Functional trait diversity is a popular tool in modern ecology, mainly used to infer assembly processes and ecosystem functioning. Patterns of functional trait diversity are shaped by ecological processes such as environmental filtering, species interactions and dispersal that are inherently spatial, and different processes may operate at different spatial scales. Adding a spatial dimension to the analysis of functional trait diversity may thus increase our ability to infer community assembly processes and to predict change in assembly processes following disturbance or land‐use change. Richness, evenness and divergence of functional traits are commonly used indices of functional trait diversity that are known to respond differently to large‐scale filters related to environmental heterogeneity and dispersal and fine‐scale filters related to species interactions (competition). Recent developments in spatial statistics make it possible to separately quantify large‐scale patterns (variation in local means) and fine‐scale patterns (variation around local means) by decomposing overall spatial autocorrelation quantified by Moran's coefficient into its positive and negative components using Moran eigenvector maps (MEM). We thus propose to identify the spatial signature of multiple ecological processes that are potentially acting at different spatial scales by contrasting positive and negative components of spatial autocorrelation for each of the three indices of functional trait diversity. We illustrate this approach with a case study from riparian plant communities, where we test the effects of disturbance on spatial patterns of functional trait diversity. The fine‐scale pattern of all three indices was increased in the disturbed versus control habitat, suggesting an increase in local scale competition and an overall increase in unexplained variance in the post‐disturbance versus control community. Further research using simulation modeling should focus on establishing the proposed link between community assembly rules and spatial patterns of functional trait diversity to maximize our ability to infer multiple processes from spatial community structure.     

Integrating trait and phylogenetic distances to assess scale‐dependent community assembly processes

Biodiversity is structured by multiple mechanisms that are dependent, at least in part, on ecological similarities and differences among species. Integrating traits and phylogenies in diversity metrics may provide deeper insight into community assembly processes across spatial scales. However, different traits are influenced by processes at different spatial scales, and it is not clear how trait‐spatial scale mismatches skew our ability to detect assembly patterns. An additional complexity is how phylogenetic distances, which might capture unmeasured traits, reflect spatially dependent processes. Here we analyze a freshwater zooplankton dataset from 91 ponds and show that different traits are associated with processes at different spatial scales. We first assessed the response of individual traits to processes at both α‐ and β‐scales, and then quantified the power of different combinations of traits and phylogenetic distances to reveal environmental and spatial drivers of α‐ and β‐diversity. We found that explanatory power was maximised when we accounted for environmental and spatial drivers with single, but different traits for α‐ and β‐diversity. Using the most appropriate trait for each spatial scale outperformed phylogenetic information, but phylogenetic information outperformed the same traits when these were used at the wrong spatial scale, and all outperformed taxonomic analyses that ignore trait and phylogenetic information. We demonstrate that accounting for species’ similarities and differences provides important information about dominant assembly mechanisms at different spatial scales, and that phylogeny is especially useful when measured traits are uninformative at a given spatial scale or when there is lack of trait data. Our study also indicates, however, that trait‐scale mismatches among phylogenetically conserved traits may affect the performance of phylogenetic indices compared to indices that account only for the best single trait at each spatial scale.     

Community assembly along a species pool gradient: implications for multiple-scale patterns of species diversity

Various ecological processes influence patterns of species diversity at multiple spatial scales. One process that is potentially important but rarely considered is community assembly. I assembled model communities using species pools of differing size to examine how the history of community assembly may affect multi-scale diversity patterns. The model contained three scales at which diversity could be measured: local community, metacommunity, and species pool. Local species saturation occurred, as expected from the competition and predation built in the model. However, local communities did not become resistant to invasions except when the species pool was very small. Depending on dispersal rate and trophic level, the larger the species pool, the harder it was to predict which species invades which local community at a given time. Consequently, local-community dissimilarity maintained by assembly history increased linearly with pool size, even though local diversity was decoupled from pool size. These results have two implications for multi-scale diversity patterns. First, assembly history may provide an explanation for scale-dependent relationships between local and regional diversity: assembly causes the relationship to be curvilinear at one scale (local community), while linear at another (metacommunity). Second, assembly history influences how -diversity is partitioned into - and -diversity: assembly causes the relative contribution of to increase with pool size. Overall, this study suggests that community assembly history interacts with species pool size to regulate multi-scale patterns of species diversity.     

Host patch traits have scale-dependent effects on diversity in a stickleback parasite metacommunity

Many metacommunities are distributed across habitat patches that are themselves aggregated into groups. Perhaps the clearest example of this nested metacommunity structure comes from multi-species parasite assemblages, which occupy individual hosts that are aggregated into host populations. At both spatial scales, we expect parasite community diversity in a given patch (either individual host or population) to depend on patch characteristics that affect colonization rates and species sorting. But, are these patch effects consistent across spatial scales? Or, do different processes govern the distribution of parasite community diversity among individual hosts, versus among host patches? To answer these questions, we document the distribution of parasite richness among host individuals and among populations in a metapopulation of threespine stickleback Gasterosteus aculeatus. We find some host traits (host size, gape width) are associated with increased parasite richness at both spatial scales. Other patch characteristics affect parasite richness only among individuals (sex), or among populations (lake size, lake area, elevation and population mean heterozygosity). These results demonstrate that some rules governing parasite richness in this metacommunity are shared across scales, while others are scale-specific.     

A metacommmunity approach to co-occurrence patterns and the core-satellite hypothesis in a community of tropical arboreal ants

Inferring mechanisms of community assembly from co-occurrence patterns is difficult in systems where assembly processes occur at multiple spatial scales and among species with heterogeneous dispersal abilities. Here, we demonstrate that local scale analysis of co-occurrence patterns is inadequate to fully describe assembly mechanisms and instead utilize a metacommunity and core-satellite approach. We generated a co-occurrence and life-history data set for a community of twig-nesting ants on coffee plants across 36 sites within a tropical agroecosystem to test the following three hypotheses: (1) twig-nesting ant species compete for nest-sites, (2) they are structured as a metacommunity, and (3) core species show segregated patterns, while satellite species show random patterns of co-occurrence. Species were divided into four groups: core species that are well distributed regionally and dominant locally, regional dominants that are well distributed regionally but do not dominate locally, local dominants that are dominant locally but are not widely distributed, and satellites that are neither widely distributed nor dominant locally. Only the most abundant species in the community, Pseudomyrmex simplex, was classified as a core species. Regional dominants, local dominants, and satellite species show random patterns of co-occurrence. However, when P. simplex is included in the co-occurrence matrix, patterns become aggregated for all three species groupings. This suggests that P. simplex “assembles” the community by providing a core metapopulation that other species track. Analyzing co-occurrence patterns among candidate subsets of species, at multiple spatial scales, and linking them to species traits substantially improves the explanatory power of co-occurrence analyses in complex metacommunities.     

Ontogenetic trait variation and metacommunity effects influence species relative abundances during tree community assembly

Predicting species abundance is one of the most fundamental pursuits of ecology. Combining the information encoded in functional traits and metacommunities provides a new perspective to predict the abundance of species in communities. We applied a community assembly via trait selection model to predict quadrat-scale species abundances using functional trait variation on ontogenetic stages and metacommunity information for over 490 plant species in a subtropical forest and a lowland tropical forest in Yunnan, China. The relative importance of trait-based selection, mass effects, and stochasticity in shaping local species abundances is evaluated using different null models. We found both mass effects and trait selection contribute to local abundance patterns. Trait selection was detectable at all studied spatial scales (0.04–1 ha), with its strength stronger at larger scales and in the subtropical forest. In contrast, the importance of stochasticity decreased with spatial scale. A significant mass effect of the metacommunity was observed at small spatial scales. Our results indicate that tree community assembly is primarily driven by ontogenetic traits and metacommunity effects. Our findings also demonstrate that including ontogenetic trait variation into predictive frameworks allows ecologists to infer ecological mechanisms operating in community assembly at the individual level.     

A simulation‐based approach to understand how metacommunity characteristics influence emergent biodiversity patterns

To understand controls over biodiversity, it is necessary to take a multi‐scale approach to understand how local and regional factors affect the community assembly processes that drive emergent patterns. This need is reflected in the growing use of the metacommunity concept to interpret multi‐scale measures of biodiversity, including metrics derived from diversity partitioning (e.g. α, β and γ diversity) and variation partitioning (e.g. spatial and environmental components of compositional turnover) techniques. However, studies have shown limited success using these metrics to characterize underlying community assembly dynamics. Here we demonstrate how a metacommunity simulation package (MCSim) can be used to evaluate when and how biodiversity metrics can be used to make inferences about metacommunity characteristics. We examined a wide range of parameter settings representing ecologically relevant scenarios. We used artificial neural networks (ANNs) to assess the sensitivity of diversity and variation partitioning metrics (calculated from simulation outcomes) to metacommunity parameter settings. In the scenarios examined in this study, the niche‐neutral gradient strongly influenced most biodiversity metrics, metacommunity size exhibited a marginal influence over some metrics, and dispersal dynamics only affected a subset of variation partitioning outcomes. Variation partitioning response curves along the niche‐neutral gradient were not monotonic; however, simulation outcomes suggest other biodiversity metrics (e.g. dissimilarity saturation) can be used in combination with variation partitioning metrics to make inferences about metacommunity properties. With the growing availability of archived ecological data, we expect future work will apply simulation‐based techniques to better understand links between biodiversity and the metacommunity characteristics that are presumed to control the underlying community assembly processes.     

Habitat selection behaviour links local and regional scales in aquatic systems

The role of habitat selection behaviour in the assembly of natural communities is an increasingly important theme in ecology. At the same time, ecologists and conservation biologists are keenly interested in scale and how processes at scales from local to regional interact to determine species distributions and patterns of biodiversity. How important is habitat selection in generating observed patterns of distribution and diversity at multiple spatial scales? In theory, habitat selection in response to interacting species can generate both positive and negative covariances among species distributions and create the potential to link processes of community assembly across multiple scales. Here I demonstrate that habitat selection by treefrogs in response to the distribution of fish predators functions at both the regional scale among localities and the local scale among patches within localities, implicating habitat selection as a critical link between local communities and the regional dynamics of metacommunities in complex landscapes.     

A phylogenetic and trait‐based analysis of community assembly in a subtropical forest in central China

Despite several decades of study in community ecology, the relative importance of the ecological processes that determine species co‐occurrence across spatial scales remains uncertain. Some of this uncertainty may be reduced by studying the scale dependency of community assembly in the light of environmental variation. Phylogenetic information and functional trait information are often used to provide potentially valuable insights into the drivers of community assembly. Here, we combined phylogenetic and trait‐based tests to gain insights into community processes at four spatial scales in a large stem‐mapped subtropical forest dynamics plot in central China. We found that all of the six leaf economic traits measured in this study had weak, but significant, phylogenetic signal. Nonrandom phylogenetic and trait‐based patterns associated with topographic variables indicate that deterministic processes tend to dominate community assembly in this plot. Specifically, we found that, on average, co‐occurring species were more phylogenetically and functionally similar than expected throughout the plot at most spatial scales and assemblages of less similar than expected species could only be found on finer spatial scales. In sum, our results suggest that the trait‐based effects on community assembly change with spatial scale in a predictable manner and the association of these patterns with topographic variables, indicates the importance of deterministic processes in community assembly relatively to random processes.     

Dissecting spatiotemporal patterns of functional diversity through the lens of Darwin's naturalization conundrum

Darwin's naturalization conundrum describes the paradigm that community assembly is regulated by two opposing processes, environmental filtering and competitive interactions, which predict both similarity and distinctiveness of species to be important for establishment. Our goal is to use long‐term, large‐scale, and high‐resolution temporal data to examine diversity patterns over time and assess whether environmental filtering or competition plays a larger role in regulating community assembly processes. We evaluated Darwin's naturalization conundrum and how functional diversity has changed in the Laurentian Great Lakes fish community from 1870 to 2010, which has experienced frequent introductions of non‐native species and extirpations of native species. We analyzed how functional diversity has changed over time by decade from 1870 to 2010 at three spatial scales (regional, lake, and habitat) to account for potential noninteractions between species at the regional and lake level. We also determined which process, environmental filtering or competitive interactions, is more important in regulating community assembly and maintenance by comparing observed patterns to what we would expect in the absence of an ecological mechanism. With the exception of one community, all analyses show that functional diversity and species richness has increased over time and that environmental filtering regulates community assembly at the regional level. When examining functional diversity at the lake and habitat level, the regulating processes become more context dependent. This study is the first to examine diversity patterns and Darwin's conundrum by integrating long‐term, large‐scale, and high‐resolution temporal data at multiple spatial scales. Our results confirm that Darwin's conundrum is highly context dependent.     

Assembly mechanisms determining high species turnover in aquatic communities over regional and continental scales

Niche and neutral processes drive community assembly and metacommunity dynamics, but their relative importance might vary with the spatial scale. The contribution of niche processes is generally expected to increase with increasing spatial extent at a higher rate than that of neutral processes. However, the extent to what community composition is limited by dispersal (usually considered a neutral process) over increasing spatial scales might depend on the dispersal capacity of composing species. To investigate the mechanisms underlying the distribution and diversity of species known to have great powers of dispersal (hundreds of kilometres), we analysed the relative importance of niche processes and dispersal limitation in determining beta‐diversity patterns of aquatic plants and cladocerans over regional (up to 300 km) and continental (up to 3300 km) scales. Both taxonomic groups were surveyed in five different European regions and presented extremely high levels of beta‐diversity, both within and among regions. High beta‐diversity was primarily explained by species replacement (turnover) rather than differences in species richness (i.e. nestedness). Abiotic and biotic variables were the main drivers of community composition. Within some regions, small‐scale connectivity and the spatial configuration of sampled communities explained a significant, though smaller, fraction of compositional variation, particularly for aquatic plants. At continental scale (among regions), a significant fraction of compositional variation was explained by a combination of spatial effects (exclusive contribution of regions) and regionally‐structured environmental variables. Our results suggest that, although dispersal limitation might affect species composition in some regions, aquatic plant and cladoceran communities are not generally limited by dispersal at the regional scale (up to 300 km). Species sorting mediated by environmental variation might explain the high species turnover of aquatic plants and cladocerans at regional scale, while biogeographic processes enhanced by dispersal limitation among regions might determine the composition of regional biotas.     

From species distributions to meta‐communities

The extent that biotic interactions and dispersal influence species ranges and diversity patterns across scales remains an open question. Answering this question requires framing an analysis on the frontier between species distribution modelling (SDM), which ignores biotic interactions and dispersal limitation, and community ecology, which provides specific predictions on community and meta‐community structure and resulting diversity patterns such as species richness and functional diversity. Using both empirical and simulated datasets, we tested whether predicted occurrences from fine‐resolution SDMs provide good estimates of community structure and diversity patterns at resolutions ranging from a resolution typical of studies within reserves (250 m) to that typical of a regional biodiversity study (5 km). For both datasets, we show that the imprint of biotic interactions and dispersal limitation quickly vanishes when spatial resolution is reduced, which demonstrates the value of SDMs for tracking the imprint of community assembly processes across scales.     

Drivers of the beta diversity of spider assemblages in southern Brazilian temporary wetlands

1. The composition of local assemblages is assembled by an interplay of species sorting, mass effects and dispersal limitation processes. The contributions of assembly processes to metacommunity structure can change with ecosystem type and specificities of the study area. Spider composition is influenced by environmental features such as habitat structure and climate, and also by spatial distances between patches. However, little is known about the roles of assembly processes in spider metacommunity structure in wetlands. 2. The beta diversity patterns of spider assemblages were assessed in 24 temporary wetlands distributed along a latitudinal gradient in southern Brazil. The study also assessed the individual correspondence of beta diversity (and its turnover and richness components) with dissimilarities in habitat structure and climate, as well as with geographic distances, using Mantel and partial Mantel correlation tests and multivariate correlograms. 3. Turnover was the most important component of spider beta diversity. Mantel tests detected significant correlations of spider beta diversity with habitat structure. Partial Mantel tests detected significant relationships only between spider beta diversity (and the richness component) and geographic distances. Additionally, spider composition was more similar than chance on smaller scales. 4. These results evidenced a complex interplay of assembly processes explaining spider metacommunity structure in temporary wetlands. Although species-sorting processes associated with habitat structure were important in structuring local spider composition, mass effects and dispersal limitation across climatic zones played an important role on a broader scale.     

Disentangling the effects of local and regional processes on biodiversity patterns through taxon‐contingent metacommunity network analysis

Metacommunity theory, which has gained a central position in ecology, accounts for the role of migration in patterns of diversity among communities at different scales. Community isolation has a main role in this theory, but is difficult to estimate empirically, partly due to the taxon‐dependent nature of dispersal. Landscapes could be perceived as either fragmented or connected for organisms with contrasting dispersal abilities. Indeed, the dispersal ability of a taxon, and the spatial scale at which eco‐evolutionary processes shape local diversity, determine a taxon‐dependent metacommunity network. In this paper, we introduce a methodology using graph theory to define this taxon‐dependent metacommunity network and then to estimate the isolation of local communities. We analyzed the relative importance of local conditions versus community isolation as determinants of community richness for 25 taxa inhabiting 18 temporary ponds. Although local factors have been the foci of most previous empirical and theoretical considerations, we demonstrate that the metacommunity network is an equally important contributor to local diversity. We also found that the relative effect of local conditions and the metacommunity network depend on body size and taxon abundance. Local diversity of larger species was more affected by patch isolation, while taxon abundances were associated with positive or negative effects of isolation. Our results provide empirical support for the proposed role of metacommunity networks as determinants of community diversity and show the taxon‐dependent nature of these networks.     

Shared or distinct responses between intermediate and satellite stream fish species in an altered Amazonian River?

Environmental and spatial variables can distinctly influence the occupancy frequency distributions in stream fish. From a metacommunity context, we tested the following hypothesis, intermediate species are governed by dispersal and niche-based processes; in contrast, satellite species are governed by niche-based processes. To test this, we separately analyzed three data sets, the entire metacommunity, the intermediate species and the satellite species, using a forward selection of explanatory variables, and a partial Redundancy Distance Analysis. The fish and 31 variables of 52 stream reaches of a Brazilian river basin in the Western Amazon were collected during the dry period of 2012. The results for all of the data set revealed two different patterns: on one side, satellite species revealed that niche and dispersal-based processes were the most important; on the other side, for intermediate species and for all of the species set, only dispersal-based processes were the most important. For the data set including all of the species and the intermediate species, the variance was explained mainly by landscape scale variables. By contrast, the variance within the satellite species set was explained by local scale variables. Management efforts for intermediate species should be taking at larger scale, but they are usually less critical for the maintenance of aquatic biodiversity; on the other hand, management efforts for satellite species should be taken at smaller scale and based on specific biological and ecological information for the focal species.