Diversity patterns from sequentially restored grasslands support the ‘environmental heterogeneity hypothesis’

The ‘environmental heterogeneity hypothesis’ (EHH) has been proposed as a mechanism that enables species coexistence through resource partitioning. In accordance with this hypothesis, plant diversity is predicted to increase with variability in resources, but there has been weak support for this hypothesis from experimental studies. The objectives of this research were to 1) characterize how resource availability and heterogeneity (coefficient of variation) change as plant communities develop using sequentially restored grasslands, 2) determine if resource heterogeneity relates to plant diversity (effective number of species, richness and evenness) and 3) reveal if the strength of resource heterogeneity–diversity relationships is different among levels of resource availability. We quantified means and coefficients of variation in soil nitrate and light availability in grasslands established on former agricultural lands for different times and their relationship to plant diversity using a geostatistically‐informed design. Nitrate availability decreased exponentially with restoration age, but no directional change in nitrate heterogeneity across the chronosequence occurred due to high resource variability in some restorations. Light availability also decreased exponentially across the chronosequence, but there was no directional change in light heterogeneity. Nitrate heterogeneity was positively correlated with both plant richness and plant effective number of species at high levels of nitrate availability. However, no nitrate heterogeneity correlation was detected at low levels of nitrate availability. Light heterogeneity was positively correlated with plant effective number of species at low levels of light availability. However, no light heterogeneity correlation was detected at high levels of light availability. Plant evenness was not correlated with resource heterogeneity at any resource availability level. These results support the positive heterogeneity–diversity relationship predicted by EHH, and uniquely that this relationship develops within a decade of plant community development, but can be obscured by resource availability.     

Metagenomic and small-subunit rRNA analyses reveal the genetic diversity of bacteria, archaea, fungi, and viruses in soil

Recent studies have highlighted the surprising richness of soil bacterial communities; however, bacteria are not the only microorganisms found in soil. To our knowledge, no study has compared the diversities of the four major microbial taxa, i.e., bacteria, archaea, fungi, and viruses, from an individual soil sample. We used metagenomic and small-subunit RNA-based sequence analysis techniques to compare the estimated richness and evenness of these groups in prairie, desert, and rainforest soils. By grouping sequences at the 97% sequence similarity level (an operational taxonomic unit [OTU]), we found that the archaeal and fungal communities were consistently less even than the bacterial communities. Although total richness levels are difficult to estimate with a high degree of certainty, the estimated number of unique archaeal or fungal OTUs appears to rival or exceed the number of unique bacterial OTUs in each of the collected soils. In this first study to comprehensively survey viral communities using a metagenomic approach, we found that soil viruses are taxonomically diverse and distinct from the communities of viruses found in other environments that have been surveyed using a similar approach. Within each of the four microbial groups, we observed minimal taxonomic overlap between sites, suggesting that soil archaea, bacteria, fungi, and viruses are globally as well as locally diverse.     

Metagenomic and Small-Subunit rRNA Analyses Reveal the Genetic Diversity of Bacteria, Archaea, Fungi, and Viruses in Soil

Recent studies have highlighted the surprising richness of soil bacterial communities; however, bacteria are not the only microorganisms found in soil. To our knowledge, no study has compared the diversities of the four major microbial taxa, i.e., bacteria, archaea, fungi, and viruses, from an individual soil sample. We used metagenomic and small-subunit RNA-based sequence analysis techniques to compare the estimated richness and evenness of these groups in prairie, desert, and rainforest soils. By grouping sequences at the 97% sequence similarity level (an operational taxonomic unit [OTU]), we found that the archaeal and fungal communities were consistently less even than the bacterial communities. Although total richness levels are difficult to estimate with a high degree of certainty, the estimated number of unique archaeal or fungal OTUs appears to rival or exceed the number of unique bacterial OTUs in each of the collected soils. In this first study to comprehensively survey viral communities using a metagenomic approach, we found that soil viruses are taxonomically diverse and distinct from the communities of viruses found in other environments that have been surveyed using a similar approach. Within each of the four microbial groups, we observed minimal taxonomic overlap between sites, suggesting that soil archaea, bacteria, fungi, and viruses are globally as well as locally diverse.     

The Role of Datasets on Scientific Influence within Conflict Research

We inductively tested if a coherent field of inquiry in human conflict research emerged in an analysis of published research involving “conflict” in the Web of Science (WoS) over a 66-year period (1945–2011). We created a citation network that linked the 62,504 WoS records and their cited literature. We performed a critical path analysis (CPA), a specialized social network analysis on this citation network (~1.5 million works), to highlight the main contributions in conflict research and to test if research on conflict has in fact evolved to represent a coherent field of inquiry. Out of this vast dataset, 49 academic works were highlighted by the CPA suggesting a coherent field of inquiry; which means that researchers in the field acknowledge seminal contributions and share a common knowledge base. Other conflict concepts that were also analyzed—such as interpersonal conflict or conflict among pharmaceuticals, for example, did not form their own CP. A single path formed, meaning that there was a cohesive set of ideas that built upon previous research. This is in contrast to a main path analysis of conflict from 1957–1971 where ideas didn’t persist in that multiple paths existed and died or emerged reflecting lack of scientific coherence (Carley, Hummon, and Harty, 1993). The critical path consisted of a number of key features: 1) Concepts that built throughout include the notion that resource availability drives conflict, which emerged in the 1960s-1990s and continued on until 2011. More recent intrastate studies that focused on inequalities emerged from interstate studies on the democracy of peace earlier on the path. 2) Recent research on the path focused on forecasting conflict, which depends on well-developed metrics and theories to model. 3) We used keyword analysis to independently show how the CP was topically linked (i.e., through democracy, modeling, resources, and geography). Publically available conflict datasets developed early on helped shape the operationalization of conflict. In fact, 94% of the works on the CP that analyzed data either relied on publically available datasets, or they generated a dataset and made it public. These datasets appear to be important in the development of conflict research, allowing for cross-case comparisons, and comparisons to previous works.     

Enhancing the performance and interpretation of freshwater biological indices: An application in arid zone streams

We assessed the performance of biological indices developed for invertebrate assemblages occurring in arid zone streams: a multimetric index (MMI) and an O/E index of taxonomic completeness. Our overall goal was to advance our understanding of the factors that affect performance and interpretation of biological indices. Our specific objectives were to (1) develop biological indices that are insensitive to natural environmental gradients, (2) develop a general method to determine if the biological potential of an assessed site is adequately represented by the population of reference sites, (3) develop a robust method to select metrics for inclusion in MMIs that ensures maximum independence of metrics, and (4) determine if a fundamental sample property (the evenness of taxa counts within a sample) affects index performance. Random Forest modeling revealed that both individual metrics and taxa composition were strongly associated with natural environmental heterogeneity, which meant both the MMI and O/E index needed to be based on site-specific expectations. We produced a precise, responsive, and ecologically robust MMI by using principal components analysis to identify 7 statistically independent metrics from a list of 31 candidate assemblage-level metrics. However, the O/E index we developed was relatively imprecise compared with O/E indices developed for other regions. This imprecision may be the consequence of low predictability in local taxa composition associated with the relatively high spatial isolation of aquatic habitats within arid regions. We were also able to assess the likelihood that the biological potential of assessed sites were adequately characterized by the population of reference sites by developing and applying a multivariate, nearest-neighbor test that determined if an assessed site occurred within the environmental space of the reference site network. This approach is robust and applicable to all biological indices. We also demonstrate that the evenness of taxa counts within a sample is positively related to estimates of sample taxa richness and thus the scores of both indices. The relationship between richness and sample evenness can potentially compromise inferences regarding biological condition, and post hoc adjustments for the effects of evenness on index scores might be desirable. Further improvements in the performance and interpretation of biological indices will require simultaneous consideration of the effects of incomplete sampling on characterization of biological assemblages and the physical and biological factors that influence community assembly.     

Trait evolution, community assembly, and the phylogenetic structure of ecological communities

Taxa co-occurring in communities often represent a nonrandom sample, in phenotypic or phylogenetic terms, of the regional species pool. While heuristic arguments have identified processes that create community phylogenetic patterns, further progress hinges on a more comprehensive understanding of the interactions between underlying ecological and evolutionary processes. We created a simulation framework to model trait evolution, assemble communities (via competition, habitat filtering, or neutral assembly), and test the phylogenetic pattern of the resulting communities. We found that phylogenetic community structure is greatest when traits are highly conserved and when multiple traits influence species membership in communities. Habitat filtering produces stronger phylogenetic structure when taxa with derived (as opposed to ancestral) traits are favored in the community. Nearest-relative tests have greater power to detect patterns due to competition, while total community relatedness tests perform better with habitat filtering. The size of the local community relative to the regional pool strongly influences statistical power; in general, power increases with larger pool sizes for communities created by filtering but decreases for communities created by competition. Our results deepen our understanding of processes that contribute to phylogenetic community structure and provide guidance for the design and interpretation of empirical research.     

A comparison of statistical approaches to analyzing community convergence between natural and constructed oyster reefs

Evaluating the success of habitat creation or restoration depends primarily on the selection of appropriate goals, relevant metrics and robust analytic approaches. For intertidal oyster reefs, the goal of restoring ecological function often is as important as the production of harvestable oysters, especially since oysters are the habitat. Assessing differences in resident faunal composition between created and natural reefs is one possible metric for evaluating ecological success. Yearly changes in the resident faunal composition on constructed and natural intertidal oyster reefs at one South Carolina restoration site were analyzed with a variety of statistical approaches to determine the most effective method(s) for documenting possible convergence in the similarity of reef assemblages over time. Two datasets were defined by the level of taxonomic identification, all taxa or a subset of common taxa, and the level of taxa reduction; all taxa, taxa > 1% of total abundance, and taxa significantly contributing to variation. Data were analyzed using “classic” multivariate analysis of variance (MANOVA), null model analysis of co-occurrence (ECOSIM), nonparametric analysis of similarity (ANOSIM), and permutation tests for multivariate analysis of similarity (PERMANOVA). Taxa abundance was used to weight MANOVA and ECOSIM analyses, while the Bray-Curtis dissimilarity index was used to weight ANOSIM and PERMANOVA analyses. Initial constraints on the analytic design and data manipulations resulted in only one test where convergence of the constructed and natural reef assemblages was indicated. Prescribed reductions in the suite of taxa considered did not alter appreciably the results. The analytic approaches varied in suitability and effectiveness at discriminating among changes in compositional similarity, even when initial constraints were relaxed. MANOVA results indicated either no difference or a significant difference in resident faunal composition between reefs, but were compromised by the inability to transform the data sufficiently to test for multivariate homogeneity violations even in analyses with reduced taxa numbers. Interpretation of ECOSIM results suggested fewer taxa in common even on natural reefs and were affected by a lack of design alternatives and the possible inflation of Type I error that weighting by abundance may cause. ANOSIM results identified no significant reef treatment effects and also suffered from design constraints and an inability to generate enough permutations to test for significant differences in datasets with relatively small sample sizes. All test results from PERMANOVA analyses except one indicated unambiguously that resident faunal assemblages on constructed reefs generally were not yet similar to natural reefs even after 7 years. Convergence of constructed and natural reef resident assemblages was suggested by PERMANOVA tests only for the dataset with the fewest taxa. The negligible limitations of PERMANOVA, flexible design options, and ability to generate significance tests for small sample sizes make the approach powerful. The ongoing development of effective statistical approaches for testing the significance of taxonomic compositional changes among habitats makes the determination of whether restoration projects are successful less dependent on the choice of analytic technique. More critical, biological questions including whether convergence of taxa abundance and composition is a valid indicator of similar ecological function remain to be answered.     

Insular shifts in body size of rice frogs in the Zhoushan Archipelago, China

1. Differences in body size between mainland and island populations have been reported for reptiles, birds and mammals. Despite widespread recognition of insular shifts in body size in these taxa, there have been no reports of such body size shifts in amphibians. 2. We provide the first evidence of an insular shift in body size for an amphibian species, the rice frog Rana limnocharis. We found significant increases in body size of rice frogs on most sampled islands in the Zhoushan archipelago when compared with neighbouring mainland China. 3. Large body size in rice frogs on islands was significantly related to increased population density, in both breeding and non-breeding seasons. Increases in rice frog density were significantly related to higher resource availability on islands. Increased resource availability on islands has led to higher carrying capacities, which has subsequently facilitated higher densities and individual growth rates, resulting in larger body size in rice frogs. We also suggest that large body size has evolved on islands, as larger individuals are competitively superior under conditions of harsh intraspecific competition at high densities. 4. Increases in body size in rice frogs were not related to several factors that have been implicated previously in insular shifts in body size in other taxa. We found no significant relationships between body size of rice frogs and prey size, number of larger or smaller frog species, island area or distance of islands from the mainland. 5. Our findings contribute to the formation of a broad, repeatable ecological generality for insular shifts in body size across a range of terrestrial vertebrate taxa, and provide support for recent theoretical work concerning the importance of resource availability for insular shifts in body size.     

A METHOD FOR ASSESSING PHYLOGENETIC LEAST SQUARES MODELS FOR SHAPE AND OTHER HIGH‐DIMENSIONAL MULTIVARIATE DATA

Studies of evolutionary correlations commonly use phylogenetic regression (i.e., independent contrasts and phylogenetic generalized least squares) to assess trait covariation in a phylogenetic context. However, while this approach is appropriate for evaluating trends in one or a few traits, it is incapable of assessing patterns in highly multivariate data, as the large number of variables relative to sample size prohibits parametric test statistics from being computed. This poses serious limitations for comparative biologists, who must either simplify how they quantify phenotypic traits, or alter the biological hypotheses they wish to examine. In this article, I propose a new statistical procedure for performing ANOVA and regression models in a phylogenetic context that can accommodate high‐dimensional datasets. The approach is derived from the statistical equivalency between parametric methods using covariance matrices and methods based on distance matrices. Using simulations under Brownian motion, I show that the method displays appropriate Type I error rates and statistical power, whereas standard parametric procedures have decreasing power as data dimensionality increases. As such, the new procedure provides a useful means of assessing trait covariation across a set of taxa related by a phylogeny, enabling macroevolutionary biologists to test hypotheses of adaptation, and phenotypic change in high‐dimensional datasets.     

Hierarchical determinants of the functional richness, evenness and divergence of a vertebrate predator assemblage

Environmental factors are hypothesized to affect the functional diversity of assemblages hierarchically. First, ecological filters constrain the range of functional traits potentially displayed by an assemblage by determining its taxonomic composition. Second, some factors determine actual functional traits through the phenotypic plasticity of species. Little is known about the relative importance of each set of factors and, by using a priori functional information, most empirical studies report only the variability expected from species turnover and neglect that associated to phenotypic plasticity. Herein, we use structural equation models to assess the factors driving the functional richness, evenness, and divergence of a top‐predator assemblage faced with high variability in resource availability and assemblage structure (i.e. species richness and abundance). We measured actual functional traits (i.e. diet composition and predation pressure) in the field and contrasted the effects of environmental filters and phenotypic plasticity (i.e. behavioral responses) by controlling for species turnover and sample size with a null model. We found that a priori estimations (i.e. null‐model expectations) tended to significantly underestimate (richness and divergence) or overestimate (evenness) functional‐diversity components, explaining just a fraction (13–63%) of the variability in observed values. Furthermore, while species richness strongly affected functional richness (positively) and evenness (negatively), and resource availability slightly affected functional divergence, via compositional effects, changes in functional‐diversity components attributable to behavioral responses of predators showed little or no association with environmental variables. As a consequence, results indicated that in productive years, functionally‐distinctive species exerting relatively low predation pressure entered the assemblage, increasing functional richness and decreasing functional evenness. However, the strong behavioral responses of functionally dominant species buffered such compositional changes, affecting to different extents the three functional‐diversity components. Thus, we argue in favor of considering phenotypic plasticity in future studies of functional diversity.     

Thirty‐eight years of change in a tropical forest: plot data from Mpanga Forest Reserve,Uganda

Variations in the composition and structure of mid‐altitude, semi‐deciduous topical forest in a 0.64‐ha plot in Mpanga Forest Reserve, Uganda, are described for a 38‐year period to 2006. Repeat surveys of trees in the plot with a girth ≥30 cm at reference height in 1982, 1993 and 2006, following a baseline survey in 1968, indicate only relatively minor fluctuations in density, Shannon diversity, evenness, basal area (BA) and estimated above ground biomass. The largest trees [diameter at breast height (dbh) > 40 cm] and main canopy taxa (e.g. Celtis mildbraedii) accounted for the largest fraction of BA. Mortality was the highest amongst taxa classed as early seral, understorey and/or in the smallest size category (dbh = 9.5–20 cm), while new recruits were predominantly understorey taxa. Only one tree was recorded as felled for human use between the surveys of 1968 and 1993. In contrast, a considerable increase in anthropogenic disturbance was evident at the time of the 2006 survey, and illegal logging now poses a substantial threat to future resource availability and carbon storage in what was for a time one of the most protected areas of forest in Uganda.     

Abundance-based similarity indices and their estimation when there are unseen species in samples

A wide variety of similarity indices for comparing two assemblages based on species incidence (i.e., presence/absence) data have been proposed in the literature. These indices are generally based on three simple incidence counts: the number of species shared by two assemblages and the number of species unique to each of them. We provide a new probabilistic derivation for any incidence-based index that is symmetric (i.e., the index is not affected by the identity ordering of the two assemblages) and homogeneous (i.e., the index is unchanged if all counts are multiplied by a constant). The probabilistic approach is further extended to formulate abundance-based indices. Thus any symmetric and homogeneous incidence index can be easily modified to an abundance-type version. Applying the Laplace approximation formulas, we propose estimators that adjust for the effect of unseen shared species on our abundance-based indices. Simulation results show that the adjusted estimators significantly reduce the biases of the corresponding unadjusted ones when a substantial fraction of species is missing from samples. Data on successional vegetation in six tropical forests are used for illustration. Advantages and disadvantages of some commonly applied indices are briefly discussed.     

Assessing the affinities of fossils using canonical variates and generalized distances

Canonical variate analysis and generalized distances are commonly used multivariate statistical techniques for assessing the comparative morphology of living and fossil primates. Some common pitfalls of these methods when used to analyze fossil specimens are: (1) ignoring the possibility that a fossil belongs to a group other than one of the predefined reference samples (i.e., restricted versus unrestricted approaches to classification), (2) misinterpreting probabilities of group membership (i.e., posterior versus typicality probabilities), and (3) failing to understand how sample sizes influence multivariate ordinations in trying to effectively illustrate the morphometric affinities of a fossil (i.e., weighted versus unweighted analyses with fossils entered indirectly or directly into the analysis). To better understand canonincal variate analysis and generalized distances, the workings of these methods are portrayed graphically as a series of rotations and rescaling of data plotted in a multivariate data space. The influence of sample sizes and the importance of higher-numberer canonical variates are emphasized. Simple examples taken from the literature illustrate how the different approaches to including a fosil in canonical variate analysis may affect the multivariate results upon which physical anthropologists base their interpretations of the fossil's morphology.     

Divergent effects of land-use,propagule pressure,and climate on woody riparian invasion

Landscape-scale analyses of biological invasion are needed to understand the relative importance of environmental drivers that vary at larger scales, such as climate, propagule pressure, resource availability, and human disturbance. One poorly understood landscape-scale question is, how does human land-use influence riparian plant invasion? To evaluate the relative importance of land-use, climate, propagule pressure, and water availability in riparian invasion, we examined tamarisk (Tamarix ramosissima, T. chinensis, hybrids), Russian olive (Elaeagnus angustifolia), and Siberian elm (Ulmus pumila) occurrence, abundance, and dominance in 238 riparian sites in developed, cultivated, and undeveloped areas of four western USA river basins (281,946 km²). Temperature and propagule pressure from individuals planted nearby largely drove invasive species occurrence, whereas factors likely to affect resource availability (e.g., land-use, precipitation, streamflow intermittency) were more important to abundance and dominance, supporting the argument that species distribution models based on occurrence alone may fail to identify conditions where invasive species have the greatest impact. The role of land-use varied among taxa: urban and suburban land-use increased Siberian elm occurrence, abundance, and dominance, and urban land-use increased Russian olive occurrence, whereas suburban land-use reduced tamarisk dominance. Surprisingly, Siberian elm, which has received scant prior scientific and management attention, occurred as or more frequently than tamarisk and Russian olive (except in undeveloped areas of the Colorado River headwaters) and had higher density and dominance than tamarisk and Russian olive in developed areas. More research is needed to understand the impacts of this largely unrecognized invader on riparian ecosystem services, particularly in urban and suburban areas.     

On the comparison of the strength of morphological integration across morphometric datasets

Evolutionary morphologists frequently wish to understand the extent to which organisms are integrated, and whether the strength of morphological integration among subsets of phenotypic variables differ among taxa or other groups. However, comparisons of the strength of integration across datasets are difficult, in part because the summary measures that characterize these patterns (RV coefficient and r_PLS) are dependent both on sample size and on the number of variables. As a solution to this issue, we propose a standardized test statistic (a z‐score) for measuring the degree of morphological integration between sets of variables. The approach is based on a partial least squares analysis of trait covariation, and its permutation‐based sampling distribution. Under the null hypothesis of a random association of variables, the method displays a constant expected value and confidence intervals for datasets of differing sample sizes and variable number, thereby providing a consistent measure of integration suitable for comparisons across datasets. A two‐sample test is also proposed to statistically determine whether levels of integration differ between datasets, and an empirical example examining cranial shape integration in Mediterranean wall lizards illustrates its use. Some extensions of the procedure are also discussed.     

Resource availability influences the diversity of a functional group of heterotrophic soil bacteria

Resource availability is a key factor regulating biodiversity and ecosystem functioning, but the relationship between resource availability and diversity has only been rarely investigated in microbial communities. The aim of this study was to determine how diversity and community structure of a functional group of soil bacteria are influenced by resource concentration. To achieve this, we used soil microcosms to investigate degradation of benzoate, which served as a model compound, by soil bacterial communities. Microcosms were supplied with ¹³C-labelled benzoate at four concentrations and RNA-stable isotope probing followed by molecular fingerprinting analysis of 16S rRNA genes was employed to identify bacteria able to assimilate benzoate at different concentrations. The composition of the benzoate degrader community differed at different concentrations and there was a significant decrease in taxa evenness at the highest substrate concentration. Active organisms could be grouped into generalists, occurring at all substrate concentrations, specialists, active at one particular benzoate concentration only, and taxa that were active at either the two lowest or two highest concentrations. The study comprises the first explicit demonstration that resource availability has an effect on the diversity of a functional group of heterotrophic soil bacteria.     

Climate change and shrinking salamanders: alternative mechanisms for changes in plethodontid salamander body size

An increasing number of studies have demonstrated relationships between climate trends and body size change of organisms. In many cases, climate might be expected to influence body size by altering thermoregulation, energetics or food availability. However, observed body size change can result from a variety of ecological processes (e.g. growth, selection, population dynamics) or imperfect observation of biological systems. We used two extensive datasets to evaluate alternative mechanisms for recently reported changes in the observed body size of plethodontid salamanders. We found that mean adult body size of salamanders can be highly sensitive to survey conditions, particularly rainfall. This systematic bias in the detection of larger or smaller individuals could result in a signature of body size change in relation to reported climate trends when it is simply observation error. We also identify considerable variability in body size distributions among years and find that individual growth rates can be strongly influenced by weather. Finally, our study demonstrates that measures of mean adult body size can be highly variable among surveys and that large sample sizes may be required to make reliable inferences. Identifying the effects of climate change is a critical area of research in ecology and conservation. Researchers should be aware that observed changes in certain organisms can result from multiple ecological processes or systematic bias due to nonrandom sampling of populations.     

Winter severity determines functional trait composition of phytoplankton in seasonally ice‐covered lakes

How climate change will affect the community dynamics and functionality of lake ecosystems during winter is still little understood. This is also true for phytoplankton in seasonally ice‐covered temperate lakes which are particularly vulnerable to the presence or absence of ice. We examined changes in pelagic phytoplankton winter community structure in a north temperate lake (Müggelsee, Germany), covering 18 winters between 1995 and 2013. We tested how phytoplankton taxa composition varied along a winter‐severity gradient and to what extent winter severity shaped the functional trait composition of overwintering phytoplankton communities using multivariate statistical analyses and a functional trait‐based approach. We hypothesized that overwintering phytoplankton communities are dominated by taxa with trait combinations corresponding to the prevailing winter water column conditions, using ice thickness measurements as a winter‐severity indicator. Winter severity had little effect on univariate diversity indicators (taxon richness and evenness), but a strong relationship was found between the phytoplankton community structure and winter severity when taxon trait identity was taken into account. Species responses to winter severity were mediated by the key functional traits: motility, nutritional mode, and the ability to form resting stages. Accordingly, one or the other of two functional groups dominated the phytoplankton biomass during mild winters (i.e., thin or no ice cover; phototrophic taxa) or severe winters (i.e., thick ice cover; exclusively motile taxa). Based on predicted milder winters for temperate regions and a reduction in ice‐cover durations, phytoplankton communities during winter can be expected to comprise taxa that have a relative advantage when the water column is well mixed (i.e., need not be motile) and light is less limiting (i.e., need not be mixotrophic). A potential implication of this result is that winter severity promotes different communities at the vernal equinox, which may have different nutritional quality for the next trophic level and ecosystem‐scale effects.     

The relationship between species richness and evenness: a meta-analysis of studies across aquatic ecosystems

Biological diversity comprises both species richness, i.e., the number of species in a community, and evenness, measuring how similar species are in their abundances. The relationship between species richness and evenness (RRE) across communities remains, however, a controversial issue in ecology because no consistent pattern has been reported. We conducted a systematic meta-review of RRE in aquatic ecosystems along regional to continental gradients and across trophic groups, differing in body size by 13 orders of magnitude. Hypotheses that RRE responded to latitudinal and scale variability across trophic groups were tested by regression analyses. Significant correlations of species richness and evenness only existed in 71 out of 229 datasets. Among the RRE, 89 were negative and 140 were positive. RRE did not vary with latitude but showed a positive response to scale. In a meta-analysis with ecosystem type as a single explaining variable, RRE did not vary among ecosystem types, i.e. between marine and freshwater. Finally, autotrophs had more positive RRE than heterotrophs. The weak RRE in many aquatic datasets suggests that richness and evenness often reflect independent components of biodiversity, highlighting that richness alone may be an incomplete surrogate for biodiversity. Our results further elucidate that RRE is driven by organismal and environmental properties, both of which must be considered to gain a deeper understanding of large-scale patterns of biodiversity.