Predatory fish loss affects the structure and functioning of a model marine food web

The rate of species loss is increasing on a global scale and predators are most at risk from human-induced extinction. The effects of losing predators are difficult to predict, even with experimental single species removals, because different combinations of species interact in unpredictable ways. We tested the effects of the loss of groups of common predators on herbivore and algal assemblages in a model benthic marine system. The predator groups were fish, shrimp and crabs. Each group was represented by at least two characteristic species based on data collected at local field sites. We examined the effects of the loss of predators while controlling for the loss of predator biomass. The identity, not the number of predator groups, affected herbivore abundance and assemblage structure. Removing fish led to a large increase in the abundance of dominant herbivores, such as Ampithoids and Caprellids. Predator identity also affected algal assemblage structure. It did not, however, affect total algal mass. Removing fish led to an increase in the final biomass of the least common taxa (red algae) and reduced the mass of the dominant taxa (brown algae). This compensatory shift in the algal assemblage appeared to facilitate the maintenance of a constant total algal biomass. In the absence of fish, shrimp at higher than ambient densities had a similar effect on herbivore abundance, showing that other groups could partially compensate for the loss of dominant predators. Crabs had no effect on herbivore or algal populations, possibly because they were not at carrying capacity in our experimental system. These findings show that contrary to the assumptions of many food web models, predators cannot be classified into a single functional group and their role in food webs depends on their identity and density in 'real' systems and carrying capacities.

Body mass and sex-biased parasitism in wood mice Apodemus sylvaticus

Male sex-biased parasitism (SBP) occurs across a range of mammalian taxa and two contrasting sets of hypotheses have been suggested for its establishment. The first invokes body size per se and suggests that larger individuals are either a larger target for parasites, trade off growth at the expense of immunity or cope better with parasitism than smaller individuals. The second suggests a sex-specific handicap whereby males have reduced immunocompetence compared to females due to the immunodepressive effects of testosterone. The current study investigated whether sex-biased parasitism is driven by host 'body size' or 'sex' using a rodent-tick (Apodemus sylvaticus-. Ixodes ricinus) system. Moreover, the presence or absence of large mammals at study sites were used to control the presence of immature ticks infesting wood mice, allowing the impacts of parasitism on host body mass and female reproduction to be assessed. As expected, male mice had greater tick loads than females and analyses suggested this sex-bias was driven by body mass as opposed to sex. It is therefore likely that larger individuals are a larger target for parasites, trade off growth at the expense of immunity or adapt behavioural responses to parasitism based on their body size. Parasite load had no effect on host body mass or female reproductive output suggesting individuals may alter behaviour or life history strategies to compensate for costs incurred through parasitism. Overall, this study lends support to the 'body size' hypothesis for the formation of sex-biased parasitism.

Body mass-abundance relationships are robust to cascading effects in marine food webs

Body mass has been shown to scale negatively with abundance in a wide range of habitats and ecosystems. It is believed that this relationship has important consequences for the distribution and maintenance of energy in natural communities. Some studies have shown that the relationship between body mass and abundance may be robust to major food web perturbations, fuelling the belief that natural processes may preserve the slope of this relationship and the associated cycling of energy and nutrients. Here, we use data from a long-term experimental food web manipulation to examine this issue in a semi-natural environment. Similar communities were developed in large experimental mesocosms over a six month period. Some of the mesocosms were then subjected to species removals, based on the mean strength of their trophic interactions in the communities. In treatments where the strongest interactors were removed, a community-level trophic cascade occurred. The biomass density of invertebrates increased dramatically in these communities, which led to a suppression of primary production. In spite of these widespread changes in ecosystem functioning, the slope of the relationship between body mass and abundance remained unchanged. This was the case whether average species body mass and abundance or individual organism size spectra were considered. An examination of changes in species composition before and after the experimental manipulations revealed an important mechanism for maintaining the body mass-abundance relationship. The manipulated communities all had a higher species turnover than the intact communities, with the highest turnover in communities that experienced cascading effects. As some species increased in body mass and abundance, new species filled the available size-abundance niches that were created. This maintained the overall body mass-abundance relationship and provided a stabilising structure to these experimental communities.

Detecting the effects of diversity on measures of ecosystem function: experimental design, null models and empirical observations

Here we present a novel experimental approach to examine the relationship between diversity and ecosystem Function. We develop four null predictive models, with which to differentiate between the 'sampling effect' - the chance inclusion of a highly productive species, and 'species complementarity' - the complementary use of resources by species that differ in their niche or resource use. We investigate the effects of manipulating species and functional richness on ecosystem function in marine benthic system and using empirical data from our own experiments we illustrate the application of these models.

Trait-mediated indirect interactions in a marine intertidal system as quantified by functional responses

Studies of trait-mediated indirect interactions (TMIIs) typically focus on effects higher predators have on per capita consumption by intermediate consumers of a third, basal prey resource. TMIIs are usually evidenced by changes in feeding rates of intermediate consumers and/or differences in densities of this third species. However, understanding and predicting effects of TMIIs on population stability of such basal species requires examination of the type and magnitude of the functional responses exhibited towards them. Here, in a marine intertidal system consisting of a higher-order fish predator, the shanny Lipophrys pholis, an intermediate predator, the amphipod Echinogammarus marinus, and a basal prey resource, the isopod Jaera nordmanni, we detected TMIIs, demonstrating the importance of habitat complexity in such interactions, by deriving functional responses and exploring consequences for prey population stability. Echinogammarus marinus reacted to fish predator diet cues by reducing activity, a typical anti-predator response, but did not alter habitat use. Basal prey, Jaera nordmanni, did not respond to fish diet cues with respect to activity, distribution or aggregation behaviour. Echinogammarus marinus exhibited type II functional responses towards J. nordmanni in simple habitat, but type III functional responses in complex habitat. However, while predator cue decreased the magnitude of the type II functional response in simple habitat, it increased the magnitude of the type III functional response in complex habitat. These findings indicate that, in simple habitats, TMIIs may drive down consumption rates within type II responses, however, this interaction may remain de-stabilising for prey populations. Conversely, in complex habitats, TMIIs may strengthen regulatory influences of intermediate consumers on prey populations, whilst potentially maintaining prey population stability. We thus highlight that TMIIs can have unexpected and complex ramifications throughout communities, but can be unravelled by considering effects on intermediate predator functional response types and magnitudes.

Determining trophic niche width: An experimental test of the stable isotope approach

Determining the trophic niche width of an animal population and the relative degree to which a generalist population consists of dietary specialists are long-standing problems of ecology. It has been proposed that the variance of stable isotope values in consumer tissues could be used to quantify trophic niche width of consumer populations. However, this promising idea has not yet been rigorously tested. By conducting controlled laboratory experiments using model consumer populations (Daphnia sp., Crustacea) with controlled diets, we investigated the effect of individual- and population-level specialisation and generalism on consumer d C mean and variance values. While our experimental data follow general expectations, we extend current qualitative models to quantitative predictions of the dependence of isotopic variance on dietary correlation time, a measure for the typical time over which a consumer changes its diet. This quantitative approach allows us to pinpoint possible procedural pitfalls and critical sources of measurement uncertainty. Our results show that the stable isotope approach represents a powerful method for estimating trophic niche widths, especially when taking the quantitative concept of dietary correlation time into account. © 2012 The Authors.

Testing metabolic scaling theory using intraspecific allometries in Antarctic microarthropods

Quantitative scaling relationships among body mass, temperature and metabolic rate of organisms are still controversial, while resolution may be further complicated through the use of different and possibly inappropriate approaches to statistical analysis. We propose the application of a modelling strategy based on the theoretical approach of Akaike's information criteria and non-linear model fitting (nlm). Accordingly, we collated and modelled available data at intraspecific level on the individual standard metabolic rate of Antarctic microarthropods as a function of body mass (M), temperature (T), species identity (S) and high rank taxa to which species belong (G) and tested predictions from metabolic scaling theory (mass-metabolism allometric exponent b = 0.75, activation energy range 0.2-1.2 eV). We also performed allometric analysis based on logarithmic transformations (lm). Conclusions from lm and nlm approaches were different. Best-supported models from lm incorporated T, M and S. The estimates of the allometric scaling exponent linking body mass and metabolic rate resulted in a value of 0.696 +/- 0.105 (mean +/- 95% CI). In contrast, the four best-supported nlm models suggested that both the scaling exponent and activation energy significantly vary across the high rank taxa (Collembola, Cryptostigmata, Mesostigmata and Prostigmata) to which species belong, with mean values of b ranging from about 0.6 to 0.8. We therefore reached two conclusions: 1, published analyses of arthropod metabolism based on logarithmic data may be biased by data transformation; 2, non-linear models applied to Antarctic microarthropod metabolic rate suggest that intraspecific scaling of standard metabolic rate in Antarctic microarthropods is highly variable and can be characterised by scaling exponents that greatly vary within taxa, which may have biased previous interspecific comparisons that neglected intraspecific variability.

Patterns of synchrony in natterjack toad breeding activity and reproductive success at local and regional scales

Empirical studies of the spatiotemporal dynamics of populations are required to better understand natural fluctuations in abundance and reproductive success, and to better target conservation and monitoring programmes. In particular, spatial synchrony in amphibian populations remains little studied. We used data from a comprehensive three year study of natterjack toad Bufo calamita populations breeding at 36 ponds to assess whether there was spatial synchrony in the toad breeding activity (start and length of breeding season, total number of egg strings) and reproductive success (premetamorphic survival and production of metamorphs). We defined a novel approach to assess the importance of short-term synchrony at both local and regional scales. The approach employs similarity indices and quantifies the interaction between the temporal and spatial components of populations using mixed effects models. There was no synchrony in the toad breeding activity and reproductive success at the local scale, suggesting that populations function as individual clusters independent of each other. Regional synchrony was apparent in the commencement and duration of the breeding season and in the number of egg strings laid (indicative of female population size). Regional synchrony in both rainfall and temperature are likely to explain the patterns observed (e.g. Moran effect). There was no evidence supporting regional synchrony in reproductive success, most likely due to spatial variability in the environmental conditions at the breeding ponds, and to differences in local population fitness (e.g. fecundity). The small scale asynchronous dynamics and regional synchronous dynamics in the number of breeding females indicate that it is best to monitor several populations within a subset of regions. Importantly, variations in the toad breeding activity and reproductive success are not synchronous, and it is thus important to consider them both when assessing the conservation status of pond-breeding amphibians. © 2012 The Authors. Ecography © 2012 Nordic Society Oikos.

Stressor intensity determines antagonistic interactions between species invasion and multiple stressor effects on ecosystem functioning

Biological invasions, nutrient enrichment and ocean warming are known to threaten biodiversity and ecosystem functioning. The independent effects of these ecological stressors are well studied, however, we lack understanding of their cumulative effects, which may be additive, antagonistic or synergistic. For example, the impacts of biological invasions are often determined by environmental context, which suggests that the effects of invasive species may vary with other stressors such as pollution or climate change. This study examined the effects of an invasive seaweed (Sargassum muticum) on the structure and functioning of a benthic marine assemblage and tested explicitly whether these effects varied with nutrient enrichment and ocean warming. Overall, the presence of Sargassum muticum increased assemblage productivity rates and warming altered algal assemblage structure, which was characterised by a decrease in kelp and an increase in ephemeral green algae. The effects of Sargassum muticum on total algal biomass accumulation, however, varied with nutrient enrichment and warming producing antagonistic cumulative effects on total algal biomass accumulation. These findings show that the nature of stressor interactions may vary with stressor intensity and among response variables, which leads to less predictable consequences for the structure and functioning of communities.

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.

Ecological dynamics of extinct species in empty habitat networks. 1. The role of habitat pattern and quantity, stochasticity and dispersal

We examined a remnant host plant (Primula veris L.) habitat network that was last inhabited by the rare butterfly Hamearis lucina L. in north Wales in 1943, to assess the relative contribution of several spatial parameters to its regional extinction. We first examined relationships between P. veris characteristics and H. lucina eggs in surviving H. lucina populations, and used these to predict the suitability and potential carrying capacity of the habitat network in north Wales. This resulted in an estimate of roughly 4500 eggs (ca 227 adults). We developed a discrete space, discrete time metapopulation model to evaluate the relative contribution of dispersal distance, habitat and environmental stochasticity as possible causes of extinction. We simulated the potential persistence of the butterfly in the current network as well as in three artificial (historical and present) habitat networks that differed in the quantity (current and X3) and fragmentation of the habitat (current and aggregated). We identified that reduced habitat quantity and increased isolation would have increased the probability of regional extinction, in conjunction with environmental stochasticity and H. lucina's dispersal distance. This general trend did not change in a qualitative manner when we modified the ability of dispersing females to stay in, and find suitable habitats (by changing the size of the grid cells used in the model). Contrary to most metapopulation model predictions, system persistence declined with increasing migration rate, suggesting that the mortality of migrating individuals in fragmented landscapes may pose significant risks to system-wide persistence. Based on model predictions for the present landscape we argue that a major programme of habitat restoration would be required for a re-established metapopulation to persist for > 100 years.

Ecological dynamics of extinct species in empty habitat networks. 2. The role of host plant dynamics

This paper explores the relative effects of host plant dynamics and butterfly-related parameters on butterfly persistence. It considers an empty habitat network where a rare butterfly (Cupido minimus) became extinct in 1939 in part of its historical range in north Wales, UK. Surviving populations of the butterfly in southern Britain were visited to assess use of its host plant (Anthyllis vulneraria) in order to calibrate habitat suitability and carrying capacity in the empty network in north Wales. These data were used to deduce that only a portion ( similar to 19%) of the host plant network from north Wales was likely to be highly suitable for oviposition. Nonetheless, roughly 65,460 eggs (3273 adult equivalents) could be expected to be laid in north Wales, were the empty network to be populated at the same levels as observed on comparable plants in surviving populations elsewhere. Simulated metapopulations of C. minimus in the empty network revealed that time to extinction and patch occupancy were significantly influenced by carrying capacity, butterfly mean dispersal distance and environmental stochasticity, although for most reasonable parameter values, the model system persisted. Simulation outputs differed greatly when host plant dynamics was incorporated into the modelled butterfly dynamics. Cupido minimus usually went extinct when host plant were at low densities. In these simulations host plant dynamics appeared to be the most important determinant of the butterfly's regional extirpation. Modelling the outcome of a reintroduction programme to C. minimus variation at high quality locations, revealed that 65% of systems survived at least 100 years. Given the current amount of resources of the north Wales landscape, the persistence of C. minimus under a realistic reintroduction programme has a good chance of being successful, if carried out in conjunction with a host plant management programme.

Fractal species distributions do not produce power-law species-area relationships

We derive the species-area relationship (SAR) expected from an assemblage of fractally distributed species. If species have truly fractal spatial distributions with different fractal dimensions, we show that the expected SAR is not the classical power-law function, as suggested recently in the literature. This analytically derived SAR has a distinctive shape that is not commonly observed in nature: upward-accelerating richness with increasing area (when plotted on log-log axes). This suggests that, in reality, most species depart from true fractal spatial structure. We demonstrate the fitting of a fractal SAR using two plant assemblages (Alaskan trees and British grasses). We show that in both cases, when modelled as fractal patterns, the modelled SAR departs from the observed SAR in the same way, in accord with the theory developed here. The challenge is to identify how species depart from fractality, either individually or within assemblages, and more importantly to suggest reasons why species distributions are not self-similar and what, if anything, this can tell us about the spatial processes involved in their generation.

A metapopulation model of species boundaries

The majority, if not all, species have a limited geographic range bounded by a distribution edge. Violent ecotones such as sea coasts clearly produce edges for many species; however such ecotones, while sufficient for the formation of an edge, are not always necessary. We demonstrate this by simulation in discrete time of a spatially structured finite size metapopulation subjected to a spatial gradient in per-unit-time population extinction probability together with spatially structured dispersal and recolonisation. We find that relatively sharp edges separating a homeland or main geographical range from an outland or zone of relatively sparse and ephemeral colonisation can form in gradual environmental gradients. The form and placing of the edge is an emergent property of the metapopulation dynamics. The sharpness of the edge declines with increasing dispersal distance, and is dependent on the relative scales of dispersal distance and gradient length. The space over which the edge develops is short relative to the potential species range. The edge is robust against changes in both the shape of the environmental gradient and to a lesser extent to alterations in the kind of dispersal operating. Persistence times in the absence of environmental gradients are virtually independent of the shape of the dispersal function describing migration. The common finding of bell shaped population density distributions across geographic ranges may occur without the strict necessity of a niche mediated response to a spatially autocorrelated environment.

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.