Wave action modifies the effects of consumer diversity and warming on algal assemblages

To understand the consequences of biodiversity loss, it is necessary to test how biodiversity-ecosystem functioning relationships may vary with predicted environmental change. In particular, our understanding will be advanced by studies addressing the interactive effects of multiple stressors on the role of biodiversity across trophic levels. Predicted increases in wave disturbance and ocean warming, together with climate-driven range shifts of key consumer species, are likely to have profound impacts on the dynamics of coastal marine communities. We tested whether wave action and temperature modified the effects of gastropod grazer diversity (Patella vulgata, Littorina littorea and Gibbula umbilicalis) on algal assemblages in experimental rock pools. The presence or absence of L. littorea appeared to drive changes in microalgal and macroalgal biomass and macroalgal assemblage structure. Macroalgal biomass also decreased with increasing grazer species richness, but only when wave action was enhanced. Further, independently of grazer diversity, wave action and temperature had interactive effects on macroalgal assemblage structure. Warming also led to a reversal of grazer-macroalgal interaction strengths from negative to positive, but only when there was no wave action. Our results show that hydrodynamic disturbance can exacerbate the effects of changing consumer diversity, and may also disrupt the influence of other environmental stressors on key consumer-resource interactions. These findings suggest that the combined effects of anticipated abiotic and biotic change on the functioning of coastal marine ecosystems, although difficult to predict, may be substantial.

Choosing and using diversity indices: Insights for ecological applications from the German Biodiversity Exploratories

Biodiversity, a multidimensional property of natural systems, is difficult to quantify partly because of the multitude of indices proposed for this purpose. Indices aim to describe general properties of communities that allow us to compare different regions, taxa, and trophic levels. Therefore, they are of fundamental importance for environmental monitoring and conservation, although there is no consensus about which indices are more appropriate and informative. We tested several common diversity indices in a range of simple to complex statistical analyses in order to determine whether some were better suited for certain analyses than others. We used data collected around the focal plant Plantago lanceolata on 60 temperate grassland plots embedded in an agricultural landscape to explore relationships between the common diversity indices of species richness (S), Shannon's diversity (H'), Simpson's diversity (D1), Simpson's dominance (D2), Simpson's evenness (E), and Berger–Parker dominance (BP). We calculated each of these indices for herbaceous plants, arbuscular mycorrhizal fungi, aboveground arthropods, belowground insect larvae, and P. lanceolata molecular and chemical diversity. Including these trait-based measures of diversity allowed us to test whether or not they behaved similarly to the better studied species diversity. We used path analysis to determine whether compound indices detected more relationships between diversities of different organisms and traits than more basic indices. In the path models, more paths were significant when using H', even though all models except that with E were equally reliable. This demonstrates that while common diversity indices may appear interchangeable in simple analyses, when considering complex interactions, the choice of index can profoundly alter the interpretation of results. Data mining in order to identify the index producing the most significant results should be avoided, but simultaneously considering analyses using multiple indices can provide greater insight into the interactions in a system.

Consistent effects of consumer species loss across different habitats

Our knowledge of the effects of consumer species loss on ecosystem functioning is limited by a paucity of manipulative field studies, particularly those that incorporate inter-trophic effects. Further, given the ongoing transformation of natural habitats by anthropogenic activities, studies should assess the relative importance of biodiversity for ecosystem processes across different environmental contexts by including multiple habitat types. We tested the context-dependency of the effects of consumer species loss by conducting a 15-month field experiment in two habitats (mussel beds and rock pools) on a temperate rocky shore, focussing on the responses of algal assemblages following the single and combined removals of key gastropod grazers (Patella vulgata, P. ulyssiponensis, Littorina littorea and Gibbula umbilicalis). In both habitats, the removal of limpets resulted in a larger increase in macroalgal richness than that of either L. littorea or G. umbilicalis. Further, by the end of the study, macroalgal cover and richness were greater following the removal of multiple grazer species compared to single species removals. Despite substantial differences in physical properties and the structure of benthic assemblages between mussel beds and rock pools, the effects of grazer loss on macroalgal cover, richness, evenness and assemblage structure were remarkably consistent across both habitats. There was, however, a transient habitat-dependent effect of grazer removal on macroalgal assemblage structure that emerged after three months, which was replaced by non-interactive effects of grazer removal and habitat after 15 months. This study shows that the effects of the loss of key consumers may transcend large abiotic and biotic differences between habitats in rocky intertidal systems. While it is clear that consumer diversity is a primary driver of ecosystem functioning, determining its relative importance across multiple contexts is necessary to understand the consequences of consumer species loss against a background of environmental change.

A multivariate analysis of beta diversity across organisms and environments

We examined variability in hierarchical beta diversity across ecosystems, geographical gradients, and organism groups using multivariate spatial mixed modeling analysis of two independent data sets. The larger data set comprised reported ratios of regional species richness (RSR) to local species richness (LSR) and the second data set consisted of RSR: LSR ratios derived from nested species-area relationships. There was a negative, albeit relatively weak, relationship between beta diversity and latitude. We found only relatively subtle differences in beta diversity among the realms, yet beta diversity was lower in marine systems than in terrestrial or freshwater realms. Beta diversity varied significantly among organisms' major characteristics such as body mass, trophic position, and dispersal type in the larger data set. Organisms that disperse via seeds had highest beta diversity, and passively dispersed organisms showed the lowest beta diversity. Furthermore, autotrophs had lower beta diversity than organisms higher up the food web; omnivores and carnivores had consistently higher beta diversity. This is evidence that beta diversity is simultaneously controlled by extrinsic factors related to geography and environment, and by intrinsic factors related to organism characteristics.

The extended phenotype of Scots pine Pinus sylvestris structures the understorey assemblage

Heritable variation in plant secondary compounds in dominant species has been hypothesised to effect ecosystem function and the structure of associated assemblages of plants, microbes and animals. The functioning of this extended phenotype in relation to the understorey vegetation composition was tested within a boreal forest system dominated by Pinus sylvestris which contains a range of monoterpenes, the composition of which is largely under genetic control. A variance partitioning approach was adopted to identify the relative importance of tree chemistry, environment, spatial location and tree architecture in controlling the distribution of species in the ground flora under individual trees. The monoterpene composition of the pine needles appeared to contribute significantly to controlling understorey vegetation composition, but was less important than environmental factors, though similar to spatial factors. Thus there appears to be a link between variation in the chemical composition of the single, dominant tree species within this system and the pattern of occurrence and abundance in other species at the same trophic level.

Young, Old, and Weathered Carbon—Part 2: Using Radiocarbon and Stable Isotopes to Identify Terrestrial Carbon Support of the Food Web in an Alkaline, Humic Lake

Carbon (C) and nitrogen (N) stable isotope analysis (SIA) has been used to identify the terrestrial subsidy of freshwater food webs. However, SIA fails to differentiate between the contributions of old and recently fixed terrestrial C and consequently cannot fully determine the source, age, and biochemical quality of terrestrial carbon. Natural abundance radiocarbon (∆14C) was used to examine the age and origin of carbon in Lower Lough Erne, Northern Ireland. 14C and stable isotope values were obtained from invertebrate, algae, and fish samples, and the results indicate that terrestrial organic C is evident at all trophic levels. High winter δ15N values in calanoid zooplankton (δ15N = 24‰) relative to phytoplankton and particulate organic matter (δ15N = 6‰ and 12‰, respectively) may reflect several microbial trophic levels between terrestrial C and calanoid invertebrates. Winter and summer calanoid ∆14C values show a seasonal switch between autochthonous and terrestrial carbon sources. Fish ∆14C values indicate terrestrial support at the highest trophic levels in littoral and pelagic food webs. 14C therefore is useful in attributing the source of carbon in freshwater in addition to tracing the pathway of terrestrial carbon through the food web.

Constructing random matrices to represent real ecosystems

Models of complex systems with n components typically have order n² parameters because each component can potentially interact with every other. When it is impractical to measure these parameters, one may choose random parameter values and study the emergent statistical properties at the system level. Many influential results in theoretical ecology have been derived from two key assumptions: that species interact with random partners at random intensities and that intraspecific competition is comparable between species. Under these assumptions, community dynamics can be described by a community matrix that is often amenable to mathematical analysis. We combine empirical data with mathematical theory to show that both of these assumptions lead to results that must be interpreted with caution. We examine 21 empirically derived community matrices constructed using three established, independent methods. The empirically derived systems are more stable by orders of magnitude than results from random matrices. This consistent disparity is not explained by existing results on predator-prey interactions. We investigate the key properties of empirical community matrices that distinguish them from random matrices. We show that network topology is less important than the relationship between a species’ trophic position within the food web and its interaction strengths. We identify key features of empirical networks that must be preserved if random matrix models are to capture the features of real ecosystems.

14C as a tool to trace terrestrial carbon in a complex lake system: implications for foodweb structure and carbon cycling

Carbon and nitrogen stable isotope analysis (SIA) has identified the terrestrial subsidy of freshwater food-webs but relies on different 13C fractionation in aquatic and terrestrial primary producers. However dissolved inorganic carbon (DIC) is partly comprised of 13C depleted respiration of terrestrial C and ‘old’ C derived from weathering of catchment geology. SIA thus fails to differentiate between the contribution of old and recently fixed terrestrial C. DIC in alkaline lakes is partially derived from weathering of 14C-free carbonaceous bedrock This
yields an artificial age offset leading samples to appear significantly older than their actual age. As such, 14C can be used as a biomarker to identify the proportion of autochthonous C in the food-web. With terrestrial C inputs likely to increase, the origin and utilisation of ‘old’ or ‘recent’ allochthonous C in the food-web can also be determined. Stable isotopes and 14C were measured for biota, particulate organic matter (POM), DIC and dissolved organic carbon (DOC) from Lough Erne, Northern Ireland, a humic but alkaline lake. High winter δ15N values in calanoid zooplankton (δ15N =24‰) relative to phytoplankton and POM (δ15N =6‰ and 12‰ respectively) may reflect several microbial trophic levels between terrestrial C and calanoids. Furthermore winter calanoid 14C ages are consistent with DOC from inflowing rivers (87 and 75 years BP respectively) but not phytoplankton (355 years BP). Summer calanoid δ13N, δ15N and 14C (312 years BP) indicate greater reliance on phytoplankton. There is also temporal and spatial variation in DIC, DOC and POM C isotopes.

14C as a tool to trace terrestrial carbon in a complex lake system: implications for food-web structure and carbon cycling

Globally lakes bury and remineralise significant quantities of terrestrial C, and the associated flux of terrestrial C strongly influences their functioning. Changing deposition chemistry, land use and climate induced impacts on hydrology will affect soil biogeochemistry and terrestrial C export¹ and hence lake ecology with potential feedbacks for regional and global C cycling. C and nitrogen stable isotope analysis (SIA) has identified the terrestrial subsidy of freshwater food webs. The approach relies on different ¹³C fractionation in aquatic and terrestrial primary producers, but also that inorganic C demands of aquatic primary producers are partly met by ¹³C depleted C from respiration of terrestrial C, and ‘old’ C derived from weathering of catchment geology. SIA thus fails to differentiate between the contributions of old and recently fixed terrestrial C. Natural abundance ¹⁴C can be used as an additional biomarker to untangle riverine food webs² where aquatic and terrestrial δ ¹³C overlap, but may also be valuable for examining the age and origin of C in the lake. Primary production in lakes is based on dissolved inorganic C (DIC). DIC in alkaline lakes is partially derived from weathering of carbonaceous bedrock, a proportion of which is¹⁴C-free. The low ¹⁴C activity yields an artificial age offset leading samples to appear hundreds to thousands of years older than their actual age. As such, ¹⁴C can be used to identify the proportion of autochthonous C in the food-web. With terrestrial C inputs likely to increase, the origin and utilisation of ‘fossil’ or ‘recent’ allochthonous C in the food-web can also be determined. Stable isotopes and 14C were measured for biota, particulate organic matter (POM), DIC and dissolved organic carbon (DOC) from Lough Erne, Northern Ireland, a humic alkaline lake. Temporal and spatial variation was evident in DIC, DOC and POM C isotopes with implications for the fluctuation in terrestrial export processes. Ramped pyrolysis of lake surface sediment indicates the burial of two C components. ¹⁴C activity (507 ± 30 BP) of sediment combusted at 400˚C was consistent with algal values and younger than bulk sediment values (1097 ± 30 BP). The sample was subsequently combusted at 850˚C, yielding 14C values (1471 ± 30 BP) older than the bulk sediment age, suggesting that fossil terrestrial carbon is also buried in the sediment. Stable isotopes in the food web indicate that terrestrial organic C is also utilised by lake organisms. High winter δ ¹⁵N values in calanoid zooplankton (δ ¹⁵N = 24%¸) relative to phytoplankton and POM (δ ¹⁵N = 6h and 12h respectively) may reflect several microbial trophic levels between terrestrial C and calanoids. Furthermore winter calanoid 14C ages are consistent with DOC from an inflowing river (75 ± 24 BP), not phytoplankton (367 ± 70 BP). Summer calanoid δ ^13C, δ ¹⁵N and ¹⁴C (345 ± 80 BP) indicate greater reliance on phytoplankton.

¹ Monteith, D.T et al., (2007) Dissolved organic carbon trends resulting from changes in atmospheric deposition chemistry. Nature, 450:537-535

² Caraco, N., et al.,(2010) Millennial-aged organic carbon subsidies to a modern river food web. Ecology,91: 2385-2393.

Differential ecological impacts of an invader and a native under environmental change are revealed by comparative functional responses.

Predicting the ecological impacts of damaging invasive species under relevant environmental contexts is a major challenge, for which comparative functional responses (the relationship between resource availability and consumer uptake rate) have great potential. Here, the functional responses of Gammarus pulex, an ecologically damaging invader in freshwaters in Ireland and other islands, were compared with those of a native trophic equivalent Gammarus duebeni celticus. Experiments were conducted at two dissolved oxygen concentrations (80 and 50 % saturation), representative of anthropogenic water quality changes, using two larval prey, blackfly (Simuliidae spp.) and mayfly (Baetis rhodani). Overall, G. pulex had higher Type II functional responses and hence predatory impacts than G. d. celticus and the functional responses of both predators were reduced by lowered oxygen concentration. However, this reduction was of lower magnitude for the invader as compared to the native. Further, the invader functional response at low oxygen was comparable to that of the native at high oxygen. Attack rates of the two predators were similar, with low oxygen reducing these attack rates, but this effect occurred more strongly for blackfly than mayfly prey. Handling times were significantly lower for the invader compared with the native, and significantly higher at low oxygen, however, the effect of lowered oxygen on handling times was minimal for the invader and pronounced for the native. Maximum feeding rates were significantly greater for the invader compared with the native, and significantly reduced at low oxygen, with this effect again lesser for the invader as compared to the native. The greater functional responses of the invader corroborate with its impacts on recipient macroinvertebrate communities when it replaces the native. Further, our experiments predict that the impact of the invader will be less affected than the native under altered oxygen regimes driven by anthropogenic influences.

Ecological communities are vulnerable to realistic extinction sequences

Loss of species will directly change the structure and potentially the dynamics of ecological communities, which in turn may lead to additional species loss (secondary extinctions) due to direct and/or indirect effects (e.g. loss of resources or altered population dynamics). Furthermore, the vulnerability of food webs to repeated species loss is expected to be affected by food web topology, species interactions, as well as the order in which species go extinct. Species traits such as body size, abundance and connectivity might determine a species' vulnerability to extinction and, thus, the order in which species go primarily extinct. Yet, the sequence of primary extinctions, and their effects on the vulnerability of food webs to secondary extinctions, when species abundances are allowed to respond dynamically, has only recently become the focus of attention. Here, we analyse and compare topological and dynamical robustness to secondary extinctions of model food webs, in the face of 34 extinction sequences based on species traits. Although secondary extinctions are frequent in the dynamical approach and rare in the topological approach, topological and dynamical robustness tends to be correlated for many bottom-up directed, but not for top-down directed deletion sequences. Furthermore, removing species based on traits that are strongly positively correlated to the trophic position of species (such as large body size, low abundance, high net effect) is, under the dynamical approach, found to be as destructive as removing primary producers. Such top-down oriented removal of species are often considered to correspond to realistic extinction scenarios, but earlier studies, based on topological approaches, have found such extinction sequences to have only moderate effects on the remaining community. Thus, our result suggests that the structure of ecological communities, and therefore the integrity of important ecosystem processes could be more vulnerable to realistic extinction sequences than previously believed.

Size-balanced community reorganization in response to nutrients and warming

It is widely accepted that global warming will adversely affect ecological communities. As ecosystems are simultaneously exposed to other anthropogenic influences, it is important to address the effects of climate change in the context of many stressors. Nutrient enrichment might offset some of the energy demands that warming can exert on organisms by stimulating growth at the base of the food web. It is important to know whether indirect effects of warming will be as ecologically significant as direct physiological effects. Declining body size is increasingly viewed as a universal response to warming, with the potential to alter trophic interactions. To address these issues, we used an outdoor array of marine mesocosms to examine the impacts of warming, nutrient enrichment and altered top-predator body size on a community comprised of the predator (shore crab Carcinus maenas), various grazing detritivores (amphipods) and algal resources. Warming increased mortality rates of crabs, but had no effect on their moulting rates. Nutrient enrichment and warming had near diametrically opposed effects on the assemblage, confirming that the ecological effects of these two stressors can cancel each other out. This suggests that nutrient-enriched systems might act as an energy refuge to populations of species under metabolic constraints due to warming. While there was a strong difference in assemblages between mesocosms containing crabs compared to mesocosms without crabs, decreasing crab size had no detectable effect on the amphipod or algal assemblages. This suggests that in allometrically balanced communities, the expected long-term effect of warming (declining body size) is not of similar ecological consequence to the direct physiological effects of warming, at least not over the six week duration of the experiment described here. More research is needed to determine the long-term effects of declining body size on the bioenergetic balance of natural communities.

Impact of biodiversity loss on production in complex marine food webs mitigated by prey-release

Public concern over biodiversity loss is often rationalized as a threat to ecosystem functioning, but biodiversity-ecosystem functioning (BEF) relations are hard to empirically quantify at large scales. We use a realistic marine food-web model, resolving species over five trophic levels, to study how total fish production changes with species richness. This complex model predicts that BEF relations, on average, follow simple Michaelis-Menten curves when species are randomly deleted. These are shaped mainly by release of fish from predation, rather than the release from competition expected from simpler communities. Ordering species deletions by decreasing body mass or trophic level, representing 'fishing down the food web', accentuates prey-release effects and results in unimodal relationships. In contrast, simultaneous unselective harvesting diminishes these effects and produces an almost linear BEF relation, with maximum multispecies fisheries yield at approximate to 40% of initial species richness. These findings have important implications for the valuation of marine biodiversity.

Comparative functional responses predict the invasiveness and ecological impacts of alien herbivorous snails.

Understanding determinants of the invasiveness and ecological impacts of alien species is amongst the most sought-after and urgent research questions in ecology. Several studies have shown the value of comparing the functional responses (FRs) of alien and native predators towards native prey, however, the technique is under-explored with herbivorous alien species and as a predictor of invasiveness as distinct from ecological impact. Here, in China, we conducted a mesocosm experiment to compare the FRs among three herbivorous snail species: the golden apple snail, Pomacea canaliculata, a highly invasive and high impact alien listed in “100 of the World's Worst Invasive Alien Species”; Planorbarius corneus, a non-invasive, low impact alien; and the Chinese native snail, Bellamya aeruginosa, when feeding on four locally occurring plant species. Further, by using a numerical response equation, we modelled the population dynamics of the snail consumers. For standard FR parameters, we found that the invasive and damaging alien snail had the highest “attack rates” a, shortest “handling times” h and also the highest estimated maximum feeding rates, 1/hT, whereas the native species had the lowest attack rates, longest handling times and lowest maximum feeding rates. The non-invasive, low impact alien species had consistently intermediate FR parameters. The invasive alien species had higher population growth potential than the native snail species, whilst that of the non-invasive alien species was intermediate. Thus, while the comparative FR approach has been proposed as a reliable method for predicting the ecological impacts of invasive predators, our results further suggest that comparative FRs could extend to predict the invasiveness and ecological impacts of alien herbivores and should be explored in other taxa and trophic groups to determine the general utility of the approach.

Deriving ecological quality standards for lakes using hypolimnetic oxygen depletion

The ecological quality of lakes and other surface water bodies in the European Union is determined by the quality of the structure and functioning of the aquatic ecosystem. The depletion rate of oxygen in the hypolimnion is an important process in thermally stratified lakes and the distribution of consumption between water and sediment an important structural characteristic. It is shown that the variation of volumetric oxygen consumption rate with trophic state can be used to select lake water total phosphorus and chlorophyll concentrations that correspond to changes in the functioning of the lake. Lake morphometry has little effect on this aspect of lake function and the relative amount of oxygen consumption in the water and sediment changes only a little with trophic state, most of the consumption being in the water. Suggestions for the reference condition, good and moderate ecological quality are made using the changes in this aspect of lake function and they are presented as lake water total phosphorus and chlorophyll concentration.