The effects of stochastic population dynamics on food web structure

We develop a stochastic, individual-based model for food web simulation which in the large-population limit reduces to the well-studied Webworld model, which has been used to successfully construct model food webs with several realistic features. We demonstrate that an almost exact match is found between the population dynamics in fixed food webs, and that the demographic fluctuations have systematic effects when the new model is used to construct food webs due to the presence of species with small populations.     

The impact of nonlinear functional responses on the long-term evolution of food web structure

We investigate the long-term web structure emerging in evolutionary food web models when different types of functional responses are used. We find that large and complex webs with several trophic layers arise only if the population dynamics is such that it allows predators to focus on their best prey species. This can be achieved using modified Lotka-Volterra or Holling/Beddington functional responses with effective couplings that depend on the predator's efficiency at exploiting the prey, or a ratio-dependent functional response with adaptive foraging. In contrast, if standard Lotka-Volterra or Holling/Beddington functional responses are used, long-term evolution generates webs with almost all species being basal, and with additionally many links between these species. Interestingly, in all cases studied, a large proportion of weak links result naturally from the evolution of the food webs.     

Rapid contemporary evolution and clonal food web dynamics

Character evolution that affects ecological community interactions often occurs contemporaneously with temporal changes in population size, potentially altering the very nature of those dynamics. Such eco-evolutionary processes may be most readily explored in systems with short generations and simple genetics. Asexual and cyclically parthenogenetic organisms such as microalgae, cladocerans and rotifers, which frequently dominate freshwater plankton communities, meet these requirements. Multiple clonal lines can coexist within each species over extended periods, until either fixation occurs or a sexual phase reshuffles the genetic material. When clones differ in traits affecting interspecific interactions, within-species clonal dynamics can have major effects on the population dynamics. We first consider a simple predator–prey system with two prey genotypes, parametrized with data from a well-studied experimental system, and explore how the extent of differences in defence against predation within the prey population determine dynamic stability versus instability of the system. We then explore how increased potential for evolution affects the community dynamics in a more general community model with multiple predator and multiple prey genotypes. These examples illustrate how microevolutionary ‘details’ that enhance or limit the potential for heritable phenotypic change can have significant effects on contemporaneous community-level dynamics and the persistence and coexistence of species.     

Plankton population dynamics: food web interactions and abiotic constraints

1. In this introduction, I try to follow some developments in plankton ecology, how they have led to current research topics, and how the contributions in this issue of Freshwater Biology are related to these fields of research.
2. Due to several favourable features, such as small size, short generation time and a relatively homogeneous habitat, planktonic organisms remain ideal subjects for theoretical and experimental population ecology.
3. Important current research topics involve: (1) the control of plankton communities by external abiotic factors; (2) bottom-up (limitation by resources) and top-down (control by predators) effects in the food web; (3) the importance of dormant resting stages and benthic–pelagic coupling in plankton dynamics; (4) costs and benefits of the mixotrophic strategy, i.e. the ability to combine a phototrophic and a phagotrophic mode of nutrition.     

Simulating food web dynamics along a gradient: quantifying human influence

Realistically parameterized and dynamically simulated food-webs are useful tool to explore the importance of the functional diversity of ecosystems, and in particular relations between the dynamics of species and the whole community. We present a stochastic dynamical food web simulation for the Kelian River (Borneo). The food web was constructed for six different locations, arrayed along a gradient of increasing human perturbation (mostly resulting from gold mining activities) along the river. Along the river, the relative importance of grazers, filterers and shredders decreases with increasing disturbance downstream, while predators become more dominant in governing eco-dynamics. Human activity led to increased turbidity and sedimentation which adversely impacts primary productivity. Since the main difference between the study sites was not the composition of the food webs (structure is quite similar) but the strengths of interactions and the abundance of the trophic groups, a dynamical simulation approach seemed to be useful to better explain human influence. In the pristine river (study site 1), when comparing a structural version of our model with the dynamical model we found that structurally central groups such as omnivores and carnivores were not the most important ones dynamically. Instead, primary consumers such as invertebrate grazers and shredders generated a greater dynamical response. Based on the dynamically most important groups, bottom-up control is replaced by the predominant top-down control regime as distance downstream and human disturbance increased. An important finding, potentially explaining the poor structure to dynamics relationship, is that indirect effects are at least as important as direct ones during the simulations. We suggest that our approach and this simulation framework could serve systems-based conservation efforts. Quantitative indicators on the relative importance of trophic groups and the mechanistic modeling of eco-dynamics could greatly contribute to understanding various aspects of functional diversity.     

The structure of a plant-pollinator food web

The pollination biology literature is dominated by examples of specialization between plants and their pollinators. However, a recent review shows that it is generalization that prevails in the field, with most plants having a number of pollinators and most pollinators visiting a number of plants. Consequently, the vast majority of plant–pollinator interactions are embedded in a complex web of plant–pollinator interactions. These plant-pollinator webs can be studied in the manner of conventional food webs and the aim of this paper is to illustrate how contemporary methods of web construction and analysis can be applied to plant-pollinator communities.     

Prey preference,intraguild predation and population dynamics of an arthropod food web on plants

The theory of intraguild predation (IGP) largely studies effects on equilibrium densities of predators and prey, while experiments mostly concern transient dynamics. We studied the effects of an intraguild (IG) predator, the bug Orius laevigatus, on the population dynamics of IG-prey, the predatory mite Phytoseiulus persimilis, and a shared prey, the phytophagous two-spotted spider mite Tetranychus urticae, as well as on the performance of cucumber plants in a greenhouse. The interaction of the predatory mite and the spider mite is highly unstable, and ends either by herbivores overexploiting the plant or predators exterminating the herbivores. We studied the effect of IGP on the transient dynamics of this system, and compared the dynamics with that predicted by a simple population-dynamical model with IGP added. Behavioural studies showed that the predatory bug and the predatory mite were both attracted to plants infested by spider mites and that the two predators did not avoid plants occupied by the other predator. Observations on foraging behaviour of the predatory bug showed that it attacks and kills large numbers of predatory mites and spider mites. The model predicts strong effects of predation and prey preference by the predatory bugs on the dynamics of predatory mites and spider mites. However, experiments in which the predatory bug was added to populations of predatory mites and spider mites had little or no effect on numbers of both mite species, and cucumber plant and fruit weight.     

Consequences of symbiosis for food web dynamics

Basic Lotka-Volterra type models in which mutualism (a type of symbiosis where the two populations benefit both) is taken into account, may give unbounded solutions. We exclude such behaviour using explicit mass balances and study the consequences of symbiosis for the long-term dynamic behaviour of a three species system, two prey and one predator species in the chemostat. We compose a theoretical food web where a predator feeds on two prey species that have a symbiotic relationships. In addition to a species-specific resource, the two prey populations consume the products of the partner population as well. In turn, a common predator forages on these prey populations. The temporal change in the biomass and the nutrient densities in the reactor is described by ordinary differential equations (ODE). Since products are recycled, the dynamics of these abiotic materials must be taken into account as well, and they are described by odes in a similar way as the abiotic nutrients. We use numerical bifurcation analysis to assess the long-term dynamic behaviour for varying degrees of symbiosis. Attractors can be equilibria, limit cycles and chaotic attractors depending on the control parameters of the chemostat reactor. These control parameters that can be experimentally manipulated are the nutrient density of the inflow medium and the dilution rate. Bifurcation diagrams for the three species web with a facultative symbiotic association between the two prey populations, are similar to that of a bi-trophic food chain; nutrient enrichment leads to oscillatory behaviour. Predation combined with obligatory symbiotic prey-interactions has a stabilizing effect, that is, there is stable coexistence in a larger part of the parameter space than for a bi-trophic food chain. However, combined with a large growth rate of the predator, the food web can persist only in a relatively small region of the parameter space. Then, two zero-pair bifurcation points are the organizing centers. In each of these points, in addition to a tangent, transcritical and Hopf bifurcation a global heteroclinic bifurcation is emanating. This heteroclinic cycle connects two saddle equilibria where the predator is absent. Under parameter variation the period of the stable limit cycle goes to infinity and the cycle tends to the heteroclinic cycle. At this global bifurcation point this cycle breaks and the boundary of the basin of attraction disappears abruptly because the separatrix disappears together with the cycle. As a result, it becomes possible that a stable two-nutrient–two-prey population system becomes unstable by invasion of a predator and eventually the predator goes extinct together with the two prey populations, that is, the complete food web is destroyed. This is a form of over-exploitation by the predator population of the two symbiotic prey populations. When obligatory symbiotic prey-interactions are modelled with Liebigs minimum law, where growth is limited by the most limiting resource, more complicated types of bifurcations are found. This results from the fact that the Jacobian matrix changes discontinuously with respect to a varying parameter when another resource becomes most limiting.Revised version: 21 July 2003     

The role of mixotrophic protists in the population dynamics of the microbial food web in a small artificial pond

• 1 The seasonal variations of protist and rotifer populations were monitored over 1 year in a small artificial pond. Grazing rates on fluorescently labelled bacteria were also determined.
• 2 The data showed population dynamics similar to other small freshwater bodies; diatoms were numerous during the spring, chlorophytes dominated during the summer months, and mixotrophs, in particular Gymnodinium, dominated during the autumn and winter.
• 3 The mixotrophic dinoflagellates were responsible for a high chlorophyll concentration during the autumn and winter. Mixotrophs were important consumers of bacteria, particularly during the autumn when population densities of pure heterotrophs were low. Mixotrophs were an important component of the microbial food web in this pond.

     

The phylogenetic component of food web structure and intervality

Despite the exceptional complexity formed by species and their interactions in ecological networks, such as food webs, regularities in the network structures are repeatedly demonstrated. The interactions are determined by the characteristics of a species. The characteristics are in turn determined by the species’ phylogenetic relationships, but also by factors not related to evolutionary history. Here, we test whether species’ phylogenetic relationships provides a significant proxy for food web intervality. We thereafter quantify the degree to which different species traits remain valuable predictors of food web structure after the baseline effect of species’ relatedness has been removed. We find that the phylogenetic relationships provide a significant background from which to estimate food web intervality and thereby structure. However, we also find that there is an important, non-negligible part of some traits, e.g., body size, in food webs that is not accounted for by the phylogenetic relationships. Additionally, both these relationships differ depending if a predator or a prey perspective is adopted. Clearly, species’ evolutionary history as well as traits not determined by phylogenetic relationships shapes predator-prey interactions in food webs, and the underlying evolutionary processes take place on slightly different time scales depending on the direction of predator-prey adaptations.     

The influence of recombination on the population structure and evolution of the human pathogen Neisseria meningitidis.

The extent to which recombination disrupts the bifurcating treelike phylogeny and clonal structure imposed by binary fission on bacterial populations remains contentious. Here, we address this question with a study of nucleotide sequence data from 107 isolates of the human pathogen Neisseria meningitidis. Gene fragments from 12 house-keeping loci distributed around the meningococcal chromosome were analyzed, showing that (1) identical alleles are disseminated among genetically diverse isolates, with no evidence for linkage disequilibrium; (2) different loci give distinct and incongruent phylogenetic trees; and (3) allele sequences are incompatible with a bifurcating treelike phylogeny at all loci. These observations are consistent with the hypothesis that meningococcal populations comprise organisms assembled from a common gene pool, with alleles and allele fragments spreading independently, together with the occasional importation of genetic material from other species. Further, they support the view that recombination is an important genetic mechanism in the generation new meningococcal clones and alleles. Consequently, for anything other than the short-term evolution of this species, a bifurcating treelike phylogeny is not an appropriate model.     

The effect of population structure on the rate of evolution

Ecological factors exert a range of effects on the dynamics of the evolutionary process. A particularly marked effect comes from population structure, which can affect the probability that new mutations reach fixation. Our interest is in population structures, such as those depicted by ‘star graphs’, that amplify the effects of selection by further increasing the fixation probability of advantageous mutants and decreasing the fixation probability of disadvantageous mutants. The fact that star graphs increase the fixation probability of beneficial mutations has lead to the conclusion that evolution proceeds more rapidly in star-structured populations, compared with mixed (unstructured) populations. Here, we show that the effects of population structure on the rate of evolution are more complex and subtle than previously recognized and draw attention to the importance of fixation time. By comparing population structures that amplify selection with other population structures, both analytically and numerically, we show that evolution can slow down substantially even in populations where selection is amplified.     

The influence of habitat heterogeneity on host-pathogen population dynamics

The influence of spatial heterogeneity on the population dynamics of a naturally occurring invertebrate host-pathogen system was experimentally investigated. At ten week intervals over a two year period, I quantified the spatial distribution of natural populations of the terrestrial isopod crustacean Porcellio scaber infected with the isopod iridescent virus (IIV). During the seasonally dry periods of summer and early fall in central California, isopod populations were highly aggregated and the degree of patchiness and distance between inhabited patches was greatest. Coincident with increased patchiness and patch spacing was an increase in isopod density within patches. During the wet seasons of winter and spring, isopod population patchiness, inter-patch spacing, and within-patch density was low. Seasonal changes in virus prevalence were negatively correlated with within-patch density, patchiness, and inter-patch spacing. The influence of the spatial distribution of isopods on virus prevalence was also tested in field experiments. The virus was introduced into arrays of artificial habitat patches colonized by isopods in which interpatch distance was varied. The prevalence of resulting infections was monitored at weekly intervals. In addition, dispersal rates between artificial patches and natural patches were quantified and compared. The results showed that isopods in treatments with the smallest inter-patch spacing had the highest virus prevalence, with generally lower prevalence among isopods in more widely spaced patches. The spacing of experimental patches significantly affected virus prevalence, although the experiments did not resolve a clear relationship between patch spacing and virus prevalence. Rates of dispersal between patches decreased with increased patch spacing, and these rates did not differ significantly from dispersal between natural patches. The results suggest that rates of dispersal between isopod subpopulations may be an important component of the infection dynamics in this system. I discuss the consequences of these findings for host-pathogen dynamics in fragmented habitats, and for other ecological interactions in spatially heterogeneous habitats.     

Successional dynamics in the seasonally forced diamond food web

Plankton seasonal succession is a classic example of nonequilibrium community dynamics. Despite the fact that it has been well studied empirically, it lacks a general quantitative theory. Here we investigate a food web model that includes a resource, two phytoplankton, and a shared grazer-the diamond food web-in a seasonal environment. The model produces a number of successional trajectories that have been widely discussed in the context of the verbal Plankton Ecology Group model of succession, such as a spring bloom of a good competitor followed by a grazer-induced clear-water phase, setting the stage for the late-season dominance of a grazer-resistant species. It also predicts a novel, counterintuitive trajectory where the grazer-resistant species has both early- and late-season blooms. The model often generates regular annual cycles but sometimes produces multiyear cycles or chaos, even with identical forcing each year. Parameterizing the model, we show how the successional trajectory depends on nutrient supply and the length of the growing season, two key parameters that vary among water bodies. This model extends nonequilibrium theory to food webs and is a first step toward a quantitative theory of plankton seasonal succession.     

The influence of food availability on population dynamics of a supralittoral isopod, Ligia dilatata brandt

Size frequency analysis of monthly samples of the isopod Ligia dilatata Brandt in the rocky supralittoral region on the west coast of the Cape Peninsula, South Africa, over a 20-month period shows that L. dilatata lives for 2 yr and that annual recruitment occurs from spring until autumn. Females start reproducing at 12 months but, unlike males, probably do not survive to breed twice. The brood period is 5–6 wk. Fecundity is described by the regression equation N = 1.9L – 3.7 (N =no. of eggs in brood pouch, L = length of females in mm). Growth is slow during summer but faster in winter when food is more plentiful. Two distinct phases of mortality are evident; a constant slow mortality for the first 9 months after recruitment, and a faster second phase culminating in the elimination of the cohort after a further 11 months. The in situ temperature range is only 6.2°C (mean: 14.6°C) due to the insulating effect of kelp debris and appears to be less important than food availability in influencing the regular annual cycles observed.     

Chaotic dynamics of a food web in a chemostat

We analyze a mathematical model of a simple food web consisting of one predator and two prey populations in a chemostat. Monod's model is employed for the dependence of the specific growth rates of the two prey populations on the concentration of the rate-limiting substrate and a generalization of Monod's model for the dependence of the specific growth rate of the predator on the concentrations of the prey populations. We use numerical bifurcation techniques to determine the effect of the operating conditions of the chemostat on the dynamics of the system and construct its operating diagram. Chaotic behavior resulting from successive period doublings is observed. Multistability phenomena of coexistence of steady and periodic states at the same operating conditions are also found.     

The effects of mixotrophy on the stability and dynamics of a simple planktonic food web model

Recognition of the microbial loop as an important part of aquatic ecosystems disrupted the notion of simple linear food chains. However, current research suggests that even the microbial loop paradigm is a gross simplification of microbial interactions due to the presence of mixotrophs-organisms that both photosynthesize and graze. We present a simple food web model with four trophic species, three of them arranged in a food chain (nutrients-autotrophs-herbivores) and the fourth as a mixotroph with links to both the nutrients and the autotrophs. This model is used to study the general implications of inclusion of the mixotrophic link in microbial food webs and the specific predictions for a parameterization that describes open ocean mixed layer plankton dynamics. The analysis indicates that the system parameters reside in a region of the parameter space where the dynamics converge to a stable equilibrium rather than displaying periodic or chaotic solutions. However, convergence requires weeks to months, suggesting that the system would never reach equilibrium in the ocean due to alteration of the physical forcing regime. Most importantly, the mixotrophic grazing link seems to stabilize the system in this region of the parameter space, particularly when nutrient recycling feedback loops are included.     

Stability of pitcher-plant microfaunal populations depends on food web structure

Enrichment (increasing K) destabilizes simple consumer–resource interactions, but increasing food web complexity in various ways can remove this paradox of enrichment. We varied resources and number of omnivorous predators (mosquitoes) and tested for effects on the stability (persistence and temporal variability) of microfaunal populations living in pitcher plants. Top-down (omnivorous) effects were destabilizing, decreasing the persistence time of a rotifer, Habrotrocha rosa, and perhaps a microflagellate, Bodo sp. Enrichment effects were more complex, in part due to effects of shredding midges on resource availability, and in part due to interactions with predation. The persistence of Bodo increased with resource availability (more bacteria due to shredding by midges; no paradox of enrichment). Increasing resources by adding ants decreased persistence of H. rosa when mosquitoes were rare (paradox of enrichment), but the effect was reversed in leaves with significant colonization by mosquitoes. Thus, in the microfaunal community of pitcher plants, omnivorous predation tends to be destabilizing, and also tends to remove the paradox of enrichment.     

Relevance of evolutionary history for food web structure

Explaining the structure of ecosystems is one of the great challenges of ecology. Simple models for food web structure aim at disentangling the complexity of ecological interaction networks and detect the main forces that are responsible for their shape. Trophic interactions are influenced by species traits, which in turn are largely determined by evolutionary history. Closely related species are more likely to share similar traits, such as body size, feeding mode and habitat preference than distant ones. Here, we present a theoretical framework for analysing whether evolutionary history--represented by taxonomic classification--provides valuable information on food web structure. In doing so, we measure which taxonomic ranks better explain species interactions. Our analysis is based on partitioning of the species into taxonomic units. For each partition, we compute the likelihood that a probabilistic model for food web structure reproduces the data using this information. We find that taxonomic partitions produce significantly higher likelihoods than expected at random. Marginal likelihoods (Bayes factors) are used to perform model selection among taxonomic ranks. We show that food webs are best explained by the coarser taxonomic ranks (kingdom to class). Our methods provide a way to explicitly include evolutionary history in models for food web structure.     

Biological vs. physical mixing effects on benthic food web dynamics

Biological particle mixing (bioturbation) and solute transfer (bio-irrigation) contribute extensively to ecosystem functioning in sediments where physical mixing is low. Macrobenthos transports oxygen and organic matter deeper into the sediment, thereby likely providing favourable niches to lower trophic levels (i.e., smaller benthic animals such as meiofauna and bacteria) and thus stimulating mineralisation. Whether this biological transport facilitates fresh organic matter assimilation by the metazoan lower part of the food web through niche establishment (i.e., ecosystem engineering) or rather deprives them from food sources, is so far unclear. We investigated the effects of the ecosystem engineers Lanice conchilega (bio-irrigator) and Abra alba (bioturbator) compared to abiotic physical mixing events on survival and food uptake of nematodes after a simulated phytoplankton bloom. The (13)C labelled diatom Skeletonema costatum was added to 4 treatments: (1) microcosms containing the bioturbator, (2) microcosms containing the bio-irrigator, (3) control microcosms and (4) microcosms with abiotic manual surface mixing. Nematode survival and subsurface peaks in nematode density profiles were most pronounced in the bio-irrigator treatment. However, nematode specific uptake (Δδ(13)C) of the added diatoms was highest in the physical mixing treatment, where macrobenthos was absent and the diatom (13)C was homogenised. Overall, nematodes fed preferentially on bulk sedimentary organic material rather than the added diatoms. The total C budget (μg C m(-2)), which included TO(13)C remaining in the sediment, respiration, nematode and macrobenthic uptake, highlighted the limited assimilation by the metazoan benthos and the major role of bacterial respiration. In summary, bioturbation and especially bio-irrigation facilitated the lower trophic levels mainly over the long-term through niche establishment. Since the freshly added diatoms represented only a limited food source for nematodes, the macrobenthic effect was more pronounced in niche establishment than the negative structuring effects such as competition.