Improving the assessment of species compositional dissimilarity in a priori ecological classifications: evaluating map scale,sampling intensity and improvement in a hierarchical classification

Question: Can species compositional dissimilarity analyses be used to assess and improve the representation of biodiversity patterns in a priori ecological classifications? Location: The case study examined the northern‐half of the South‐east Queensland Bioregion, eastern Australia. Methods: Site‐based floristic presence–absence data were used to construct species dissimilarity matrices (Kulczynski metric) for three levels of Queensland's bioregional hierarchy – subregions (1:500 000 scale), land zones (1:250 000 scale) and regional ecosystems (1:100 000 scale). Within‐ and between‐class dissimilarities were compiled for each level to elucidate species compositional patterns. Randomized subsampling was used to determine the minimum site sampling intensity for each hierarchy level, and the effects of lumping and splitting illustrated for several classes. Results: Consistent dissimilarity estimates were obtained with five or more sites per regional ecosystem, 10 or more sites per land zone, and more than 15 sites per subregion. On average, subregions represented 4% dissimilarity in floristic composition, land zones approximately 10%, and regional ecosystems over 19%. Splitting classes with a low dissimilarity increased dissimilarity levels closer to average, while merging ecologically similar classes with high dissimilarities reduced dissimilarity levels closer to average levels. Conclusions: This approach demonstrates a robust and repeatable means of analysing species compositional dissimilarity, determining site sampling requirements for classifications and guiding decisions about ‘lumping’ or ‘splitting’ of classes. This will allow more informed decisions on selecting and improving classifications and map scales in an ecologically and statistically robust manner.     

Compositional dissimilarity as a robust measure of ecological distance

The robustness of quantitative measures of compositional dissimilarity between sites was evaluated using extensive computer simulations of species' abundance patterns over one and two dimensional configurations of sample sites in ecological space. Robustness was equated with the strength over a range of models, of the linear and monotonic (rank-order) relationship between the compositional dissimilarities and the corresponding Euclidean distances between sites measured in the ecological space. The range of models reflected different assumptions about species' response curve shape, sampling pattern of sites, noise level of the data, species' interactions, trends in total site abundance, and beta diversity of gradients.The Kulczynski, Bray-Curtis and Relativized Manhattan measures were found to have not only a robust monotonic relationship with ecological distance, but also a robust linear (proportional) relationship until ecological distances became large. Less robust measures included Chord distance, Kendall's coefficient, Chisquared distance, Manhattan distance, and Euclidean distance.A new ordination method, hybrid multidimensional scaling (HMDS), is introduced that combines metric and nonmetric criteria, and so takes advantage of the particular properties of robust dissimilarity measures such as the Kulczynski measure.We thank M. P. Austin for encouraging this study, and I. C. Prentice, E. Van der Maarel, and an anonymous reviewer for helpful comments. E. M. Adomeit provided technical assistance.     

Beta-diversity of ectoparasites at two spatial scales: nested hierarchy,geography and habitat type

Beta-diversity of biological communities can be decomposed into (a) dissimilarity of communities among units of finer scale within units of broader scale and (b) dissimilarity of communities among units of broader scale. We investigated compositional, phylogenetic/taxonomic and functional beta-diversity of compound communities of fleas and gamasid mites parasitic on small Palearctic mammals in a nested hierarchy at two spatial scales: (a) continental scale (across the Palearctic) and (b) regional scale (across sites within Slovakia). At each scale, we analyzed beta-diversity among smaller units within larger units and among larger units with partitioning based on either geography or ecology. We asked (a) whether compositional, phylogenetic/taxonomic and functional dissimilarities of flea and mite assemblages are scale dependent; (b) how geographical (partitioning of sites according to geographic position) or ecological (partitioning of sites according to habitat type) characteristics affect phylogenetic/taxonomic and functional components of dissimilarity of ectoparasite assemblages and (c) whether assemblages of fleas and gamasid mites differ in their degree of dissimilarity, all else being equal. We found that compositional, phylogenetic/taxonomic, or functional beta-diversity was greater on a continental rather than a regional scale. Compositional and phylogenetic/taxonomic components of beta-diversity were greater among larger units than among smaller units within larger units, whereas functional beta-diversity did not exhibit any consistent trend regarding site partitioning. Geographic partitioning resulted in higher values of beta-diversity of ectoparasites than ecological partitioning. Compositional and phylogenetic components of beta-diversity were higher in fleas than mites but the opposite was true for functional beta-diversity in some, but not all, traits.     

Environmental drivers of species composition and functional diversity of dung beetles along the Atlantic Forest–Pampa transition zone

Human activities are causing a rapid loss of biodiversity, which impairs ecosystem functions and services. Therefore, understanding which processes shape how biodiversity is distributed along spatial and environmental gradients is a first step to guide conservation and management efforts. We aimed to determine the relative explanatory importance of biogeographic, environmental, landscape and spatial variables on assemblage dissimilarities and functional diversity of dung beetles along the Atlantic Forest–Pampa (i.e. forest–grassland) transition zone located in Southeast South America. We described each site according to their biogeographic position, environmental conditions, landscape features and spatial patterns. The compositional dissimilarity was partitioned into turnover and nestedness components of β‐diversity. Mantel tests and generalised dissimilarity models were used to relate β‐diversity and its components to biogeographic, environmental, landscape and spatial variables. Variation partitioning analysis was used to estimate the pure and shared variation in species composition and functional diversity explained by the four categories of predictors. Biome domain was the main factor causing dung beetle compositional dissimilarity, with a high species replacement between Atlantic Forest and Pampa. Biogeographic, environmental, landscape and spatial distances also affected the patterns of dung beetle dissimilarity and β‐diversity components. The shared effects of the four sets of predictors explained most of the variation in dung beetle composition. A similar response pattern was found for dung beetle functional diversity, which excluded biogeographic effects. Only the pure effects of environmental and spatial predictors were significant for species composition and functional diversity. Our results indicate that dung beetle species composition and functional diversity are jointly driven by environmental, landscape and spatial predictors with higher pure environmental and spatial effects. The forest–grassland transition zone promotes a strong species and trait replacement highly influenced both by environmental filtering and dispersal limitation.     

Forest complexity drives dung beetle assemblages along an edge-interior gradient in the southwest Amazon rainforest

1. The habitat heterogeneity hypothesis predicts that heterogeneous habitats may provide more niches and diverse ways of exploiting environmental resources, thereby allowing more species to coexist, persist and diversify. 2. We aimed to investigate how an edge-interior gradient related to forest complexity influences species composition, abundance and richness of dung beetles in the western Amazon rainforest. We expected dung beetle abundance and richness to increase along the forest edge-interior gradient, in accordance with the habitat heterogeneity hypothesis. We also expected strong changes in species composition driven by species turnover in the forest interior and nestedness along the forest edges. We sampled dung beetles using baited pitfall traps across an edge-interior gradient. We also assessed the variation in forest features along the edge-interior gradient to identify changes in forest complexity. 3. Both species richness and abundance of dung beetles increased along the forest edge-interior, following the gradient of forest complexity. The Sorensen dissimilarity of dung beetle assemblages was higher among sampling units placed near the forest edge, although neither turnover, nor nestedness was different between the extremes of the forest edge-interior gradient. There was a clear compositional change along the edge-interior gradient mostly driven by species turnover. Individual indicator value analysis revealed that species were strongly associated with the forest interior conditions. 4. The simplification of the Amazon rainforest near clearings causes compositional changes in dung beetle assemblages. These changes are characterised by species-poor and low-abundance assemblages and may impair dung beetle ecological functions and therefore forest recovery.     

A guide through a family of phylogenetic dissimilarity measures among sites

Ecological studies have now gone beyond measures of species turnover towards measures of phylogenetic and functional dissimilarity. This change of perspective has a main objective: disentangling the processes that drive species distributions from local to broad scales. A fundamental difference between phylogenetic and functional analyses is that phylogeny is intrinsically dependent on a tree‐like structure whereas functional data can, most of time, only be forced to adhere a tree structure, not without some loss of information. When the branches of a phylogenetic tree have lengths, then each evolutionary unit on these branches can be considered as a basic entity on which dissimilarities among sites should be measured. Several of the recent measures of phylogenetic dissimilarities among sites thus are traditional dissimilarity indices where species are replaced by evolutionary units. The resulting indices were named PD‐dissimilarity indices, in reference to early work on the phylogenetic diversity (PD) measure. Here I review and compare indices and ordination approaches that, although first developed to analyse the differences in the species compositions of sites, can be adapted to describe PD‐dissimilarities among sites. Using simulations of species distributions along environmental gradients, I compare indices, associated with permutation tests and null models, in their ability to reveal existing phylogenetic patterns along the gradients. As an illustration, I show that the amount of bat PD‐dissimilarities along a disturbance gradient in Selva Lacandona of Chiapas, Mexico is dependent on whether species' abundance is considered, and on the PD‐dissimilarity index used. Overall, the family of PD‐dissimilarity indices has a critical potential for future analyses of phylogenetic diversity as it benefits from decades of research on the measure of species dissimilarity. I provide clues to help to choose among many potential indices, identifying which indices satisfy minimal basic properties, and analysing their sensitivity to abundance, size, diversity and joint absences.     

A pathway for multivariate analysis of ecological communities using copulas

We describe a new pathway for multivariate analysis of data consisting of counts of species abundances that includes two key components: copulas, to provide a flexible joint model of individual species, and dissimilarity‐based methods, to integrate information across species and provide a holistic view of the community. Individual species are characterized using suitable (marginal) statistical distributions, with the mean, the degree of over‐dispersion, and/or zero‐inflation being allowed to vary among a priori groups of sampling units. Associations among species are then modeled using copulas, which allow any pair of disparate types of variables to be coupled through their cumulative distribution function, while maintaining entirely the separate individual marginal distributions appropriate for each species. A Gaussian copula smoothly captures changes in an index of association that excludes joint absences in the space of the original species variables. A permutation‐based filter with exact family‐wise error can optionally be used a priori to reduce the dimensionality of the copula estimation problem. We describe in detail a Monte Carlo expectation maximization algorithm for efficient estimation of the copula correlation matrix with discrete marginal distributions (counts). The resulting fully parameterized copula models can be used to simulate realistic ecological community data under fully specified null or alternative hypotheses. Distributions of community centroids derived from simulated data can then be visualized in ordinations of ecologically meaningful dissimilarity spaces. Multinomial mixtures of data drawn from copula models also yield smooth power curves in dissimilarity‐based settings. Our proposed analysis pathway provides new opportunities to combine model‐based approaches with dissimilarity‐based methods to enhance understanding of ecological systems. We demonstrate implementation of the pathway through an ecological example, where associations among fish species were found to increase after the establishment of a marine reserve.     

Effects of freshwater eutrophication on species and functional beta diversity of periphytic algae

Studies about beta diversity and environmental heterogeneity have shown that the strength of the environmental filtering effect may decrease with the increasing scale. These empirical results have related eutrophic aquatic environments to higher values of beta diversity, but never to dissimilarity of species and functional traits of periphytic algae. We tested the hypotheses that periphytic algae have higher dissimilarity of both species and functional traits in eutrophic environments, and that these dissimilarities are related to environmental dissimilarity. To this end, we used richness, density, and four functional traits of periphytic algae and local limnological data from wetlands in the Brazilian savanna (Cerrado). We analyzed the beta diversity and the relationship of species and functional dissimilarities with the environmental dissimilarity and geographic distances. Our hypothesis was confirmed for functional traits dissimilarity and for the importance of the environmental dissimilarity for both species and functional beta diversity. The cultural eutrophication led to a functional homogenization in urban wetlands, which indicates the establishment of species with similar ecological requirements, and consequently, similar ‘roles’ in the ecosystem, and also that sensitive species may have been replaced by tolerant species, leading to declining biodiversity.

     

Extended dissimilarity: a method of robust estimation of ecological distances from high beta diversity data

It is widely accepted that reliable ordination of ecological data requires a strong linear or ordinal relationship between the dissimilarity of sites, based on species composition, and the ecological distance between them. Certain dissimilarity measures, having the property that they take a fixed maximum value when sites have no species in common, have been shown to be strongly correlated with ecological distance. For ecological gradients of moderate length (moderate beta diversity), such measures, in conjunction with non-metric multidimensional scaling, will reliably yield successful ordinations. However, as beta diversity increases, more sites have no species in common, and such measures invariably under-estimate ecological distance for such sites. Thus ordinations of data with high species turnover (high beta diversity) may fail.Extended dissimilarities are defined using an iterative adaptation of flexible shortest path adjustment applied to the matrix of dissimilarities with fixed maximum values. By means of theoretical argument and simulations, this is shown to lead to far stronger correlations between the adjusted site dissimilarity and ecological distance for ecological gradients of greater length than previously considered. Hence ordinations of extended dissimilarities, by means of either metric or non-metric scaling techniques, are shown to outperform corresponding ordinations of unadjusted dissimilarities, with the difference increasing with increasing beta diversity.     

A family of functional dissimilarity measures for presence and absence data

Plot‐to‐plot dissimilarity measures are considered a valuable tool for understanding the complex ecological mechanisms that drive community composition. Traditional presence/absence coefficients are usually based on different combinations of the matching/mismatching components of the 2 × 2 contingency table. However, more recently, dissimilarity measures that incorporate information about the degree of functional differences between the species in both plots have received increasing attention. This is because such “functional dissimilarity measures” capture information on the species' functional traits, which is ignored by traditional coefficients. Therefore, functional dissimilarity measures tend to correlate more strongly with ecosystem‐level processes, as species influence these processes via their traits. In this study, we introduce a new family of dissimilarity measures for presence and absence data, which consider functional dissimilarities among species in the calculation of the matching/mismatching components of the 2 × 2 contingency table. Within this family, the behavior of the Jaccard coefficient, together with its additive components, species replacement, and richness difference, is examined by graphical comparisons and ordinations based on simulated data.