Clustering in non-parametric multivariate analyses.

Non-parametric multivariate analyses of complex ecological datasets are widely used. Following appropriate pre-treatment of the data inter-sample resemblances are calculated using appropriate measures. Ordination and clustering derived from these resemblances are used to visualise relationships among samples (or variables). Hierarchical agglomerative clustering with group-average (UPGMA) linkage is often the clustering method chosen. Using an example dataset of zooplankton densities from the Bristol Channel and Severn Estuary, UK, a range of existing and new clustering methods are applied and the results compared. Although the examples focus on analysis of samples, the methods may also be applied to species analysis. Dendrograms derived by hierarchical clustering are compared using cophenetic correlations, which are also used to determine optimum  in flexible beta clustering. A plot of cophenetic correlation against original dissimilarities reveals that a tree may be a poor representation of the full multivariate information. UNCTREE is an unconstrained binary divisive clustering algorithm in which values of the ANOSIM R statistic are used to determine (binary) splits in the data, to form a dendrogram. A form of flat clustering, k-R clustering, uses a combination of ANOSIM R and Similarity Profiles (SIMPROF) analyses to determine the optimum value of k, the number of groups into which samples should be clustered, and the sample membership of the groups. Robust outcomes from the application of such a range of differing techniques to the same resemblance matrix, as here, result in greater confidence in the validity of a clustering approach.

Clustering in non-parametric multivariate analyses.

Non-parametric multivariate analyses of complex ecological datasets are widely used. Following appropriate pre-treatment of the data inter-sample resemblances are calculated using appropriate measures. Ordination and clustering derived from these resemblances are used to visualise relationships among samples (or variables). Hierarchical agglomerative clustering with group-average (UPGMA) linkage is often the clustering method chosen. Using an example dataset of zooplankton densities from the Bristol Channel and Severn Estuary, UK, a range of existing and new clustering methods are applied and the results compared. Although the examples focus on analysis of samples, the methods may also be applied to species analysis. Dendrograms derived by hierarchical clustering are compared using cophenetic correlations, which are also used to determine optimum  in flexible beta clustering. A plot of cophenetic correlation against original dissimilarities reveals that a tree may be a poor representation of the full multivariate information. UNCTREE is an unconstrained binary divisive clustering algorithm in which values of the ANOSIM R statistic are used to determine (binary) splits in the data, to form a dendrogram. A form of flat clustering, k-R clustering, uses a combination of ANOSIM R and Similarity Profiles (SIMPROF) analyses to determine the optimum value of k, the number of groups into which samples should be clustered, and the sample membership of the groups. Robust outcomes from the application of such a range of differing techniques to the same resemblance matrix, as here, result in greater confidence in the validity of a clustering approach.

Simple shade plots aid better long-term choices of data pre-treatment in multivariate assemblage studies

Shade plots, simple visual representations of abundance matrices from multivariate species assemblage studies, are shown to be an effective aid in choosing an overall transformation (or other pre-treatment) of quantitative data for long-term use, striking an appropriate balance between dominant and less abundant taxa in ensuing resemblance-based multivariate analyses. Though the exposition is entirely general and applicable to all community studies, detailed illustrations of the comparative power and interpretative possibilities of shade plots are given in the case of two estuarine assemblage studies in south-western Australia: (a) macrobenthos in the upper Swan Estuary over a two-year period covering a highly significant precipitation event for the Perth area; and (b) a wide-scale spatial study of the nearshore fish fauna from five divergent estuaries. The utility of transformations of intermediate severity is again demonstrated and, with greater novelty, the potential importance seen of further mild transformation of all data after differential down-weighting (dispersion weighting) of spatially clumped' or schooled' species. Among the new techniques utilized is a two-way form of the RELATE test, which demonstrates linking of assemblage structure (fish) to continuous environmental variables (water quality), having removed a categorical factor (estuary differences). Re-orderings of sample and species axes in the associated shade plots are seen to provide transparent explanations at the species level for such continuous multivariate patterns.

Structure and diversity of intertidal benthic diatom assemblages in contrasting shores: a case study from the Tagus estuary1

The structure of intertidal benthic diatoms assemblages in the Tagus estuary was investigated during a 2-year survey, carried out in six stations with different sediment texture. Nonparametric multivariate analyses were used to characterize spatial and temporal patterns of the assemblages and to link them to the measured environmental variables. In addition, diversity and other features related to community physiognomy, such as size-class or life-form distributions, were used to describe the diatom assemblages. A total of 183 diatom taxa were identified during cell counts and their biovolume was determined. Differences between stations (analysis of similarity (ANOSIM), R=0.932) were more evident than temporal patterns (R=0.308) and mud content alone was the environmental variable most correlated to the biotic data (BEST, rho=0.863). Mudflat stations were typically colonized by low diversity diatom assemblages (H' similar to 1.9), mainly composed of medium-sized motile epipelic species (250-1,000 mu m(3)), that showed species-specific seasonal blooms (e.g., Navicula gregaria Donkin). Sandy stations had more complex and diverse diatom assemblages (H' similar to 3.2). They were mostly composed by a large set of minute epipsammic species (<250 mu m(3)) that, generally, did not show temporal patterns. The structure of intertidal diatom assemblages was largely defined by the interplay between epipelon and epipsammon, and its diversity was explained within the framework of the Intermediate Disturbance Hypothesis. However, the spatial distribution of epipelic and epipsammic life-forms showed that the definition of both functional groups should not be over-simplified.

Depth-related gradients in community structure and relatedness of bivalves and isopods in the Southern Ocean.

Despite increased research over the last decade, diversity patterns in Antarctic deep-sea benthic taxa and their driving forces are only marginally known. Depth-related patterns of diversity and distribution of isopods and bivalves collected in the Atlantic sector of the Southern Ocean are analysed. The data, sampled by epibenthic sledge at 40 deep-sea stations from the upper continental slope to the hadal zone (774 – 6348 m) over a wide area of the Southern Ocean, comprises 619 species of isopods and 81 species of bivalves,. There were more species of isopods than bivalves in all samples, and species per station varied from 2 to 85 for isopods and from 0 to 18 for bivalves. Most species were rare, with 72% of isopod species restricted to one or two stations, and 45% of bivalves. Among less-rare species bivalves tended to have wider distributions than isopods. The species richness of isopods varied with depth, showing a weak unimodal curve with a peak at 2000 – 4000 m, while the richness of bivalves did not. Multivariate analyses indicate that there are two main assemblages in the Southern Ocean, one shallow and one deep. These overlap over a large depth-range (2000 – 4000 m). Comparing analyses based on the Sørensen resemblance measure (presence/absence) and Γ+ (presence/absence incorporating relatedness among species) indicates that rare species tend to have other closely related species within the same depth band. Analysis of relatedness among species indicates that the taxonomic variety of bivalves tends to decline at depth, whereas that of isopods is maintained. This, it is speculated, may indicate that the available energy at depth is insufficient to maintain a range of bivalve life-history strategies

Relationships between biodiversity and the stability of marine ecosystems: comparisons at a European scale using meta-analysis.

The relationship between biodiversity and stability of marine benthic assemblages was investigated using existing data sets (n = 28) covering various spatial (m-km) and temporal (1973-2006) scales in different benthic habitats (emergent rock, rock pools and sedimentary habitats) through meta-analyses. Assemblage stability was estimated by measuring temporal variances of species richness, total abundance (density or % cover) and community species composition and abundance structure (using multivariate analyses). Positive relationships between temporal variability in species number and richness were generally observed at both quadrat (<1 m2) and site (100 m2) scales, while no relationships were observed by multivariate analyses. Positive relationships were also observed at the scale of site between temporal variability in species number and variability in community structure with evenness estimates. This implies that the relationship between species richness or evenness and species richness variability is slightly positive and depends on the scale of observation, suggesting that biodiversity per se is important for the stability of ecosystems. Changes within community assemblages in terms of structure are, however, generally independent of biodiversity, suggesting no effect of diversity, but the potential impact of individual species, and/or environmental factors. Except for sedimentary and rock pool habitats, no relationship was observed between temporal variation of the aggregated variable of total abundances and diversity at either scale. Overall our results emphasise that relationships depend on scale of measurements, type of habitats and the marine systems (North Atlantic and Mediterranean) considered.

Effects of a harmful algal bloom on the community ecology, movements and spatial distributions of fishes in a microtidal estuary

Harmful algal blooms can adversely affect fish communities, though their impacts are highly context-dependent and typically differ between fish species. Various approaches, comprising univariate and multivariate analyses and multimetric Fish Community Indices (FCI), were employed to characterise the perceived impacts of a Karlodinium veneficum bloom on the fish communities and ecological condition of the Swan Canning Estuary, Western Australia. The combined evidence suggests that a large proportion of the more mobile fish species in the offshore waters of the bloom-affected area relocated to other regions during the bloom. This was indicated by marked declines in mean species richness, catch rates and FCI scores in the bloom region but concomitant increases in these characteristics in more distal regions, and by pronounced and atypical shifts in the pattern of inter-regional similarities in fish community composition during the bloom. The lack of any significant changes among the nearshore fish communities revealed that bloom impacts were less severe there than in deeper, offshore waters. Nearshore habitats, which generally are in better ecological condition than adjacent offshore waters in this system, may provide refuges for fish during algal blooms and other perturbations, mirroring similar observations of fish avoidance responses to such stressors in estuaries worldwide.

Effects of a harmful algal bloom on the community ecology, movements and spatial distributions of fishes in a microtidal estuary

Harmful algal blooms can adversely affect fish communities, though their impacts are highly context-dependent and typically differ between fish species. Various approaches, comprising univariate and multivariate analyses and multimetric Fish Community Indices (FCI), were employed to characterise the perceived impacts of a Karlodinium veneficum bloom on the fish communities and ecological condition of the Swan Canning Estuary, Western Australia. The combined evidence suggests that a large proportion of the more mobile fish species in the offshore waters of the bloom-affected area relocated to other regions during the bloom. This was indicated by marked declines in mean species richness, catch rates and FCI scores in the bloom region but concomitant increases in these characteristics in more distal regions, and by pronounced and atypical shifts in the pattern of inter-regional similarities in fish community composition during the bloom. The lack of any significant changes among the nearshore fish communities revealed that bloom impacts were less severe there than in deeper, offshore waters. Nearshore habitats, which generally are in better ecological condition than adjacent offshore waters in this system, may provide refuges for fish during algal blooms and other perturbations, mirroring similar observations of fish avoidance responses to such stressors in estuaries worldwide.

Analyses of sublittoral macrobenthic community change in a marine nature reserve using similarity profiles (SIMPROF).

Sublittoral macrobenthic communities in the Skomer Marine Nature Reserve (SMNR), Pembrokeshire, Wales, were sampled at 10 stations in 1993, 1996, 1998, 2003, 2007 and 2009 using a Day grab and a 0.5 mm mesh. The time series is analysed using Similarities Profiles (SIMPROF) tests and associated methods. Q-mode analysis using clustering with Type 1 SIMPROF addresses multivariate structure among samples, showing that there is clear structure associated with differences among years. Inverse (r-mode) analysis using Type 2 SIMPROF decisively rejects a hypothesis that species are not associated with each other. Clustering of the variables (species) with Type 3 SIMPROF identifies groups of species which covary coherently through the time-series. The time-series is characterised by a dramatic decline in abundances and diversity between the 1993 and 1996 surveys. By 1998 there had been a shift in community composition from the 1993 situation, with different species dominating. Communities had recovered in terms of abundance and species richness, but different species dominated the community. No single factor could be identified which unequivocally explained the dramatic changes observed in the SMNR. Possible causes were the effects of dispersed oil and dispersants from the Sea Empress oil spill in February 1996 and the cessation of dredge-spoil disposal off St Anne’s Head in 1995, but the most likely cause was severe weather. With many species, and a demonstrable recovery from an impact, communities within the SMNR appear to be diverse and resilient. If attributable to natural storms, the changes observed here indicate that natural variability may be much more important than is generally taken into account in the design of monitoring programmes.