Emergence of Holling type III zooplankton functional response: Bringing together field evidence and mathematical modelling

Food-web population models are rather sensitive to parameterization of functional response in predation terms. Theoretical studies predict enhancing of ecosystems’ stability for a functional response of sigmoid type (Holling type III). The choice of a correct type of response is especially important for modelling outcome of grazing control of algal blooms by zooplankton in nutrient-rich ecosystems. Extensive experiments on zooplankton feeding in laboratories show non-sigmoid nature of response for most herbivorous zooplankton species. As a consequence, there is a strong opinion in literature that the implementation of Holling III type grazing in plankton models is biologically meaningless. I argue, however, that such an ‘evident’ claim might be wrong and sigmoid functional responses in real plankton communities would emerge more often than was suggested earlier. Especially, this concerns plankton models without vertical resolution, which ignore heterogeneity in vertical distribution of species. Having conducted extensive literature search of data on zooplankton feeding in situ, I show that vertical heterogeneity in food distribution as well as active food searching behaviour of zooplankton can modify the type of functional response. In particular, the rate of food intake by the whole zooplankton population in the column, as a function of total amount of food, often exhibits a sigmoid behaviour, instead of a non-sigmoid one postulated previously based on laboratory experiments. This conceptual discrepancy is due to the ability of zooplankton to feed mostly in layers with high algal density. I propose a generic model explaining the observed alteration of type between overall and local functional responses. I show that emergence of Holling type III in plankton systems is due to mechanisms different from those well known in the ecological literature (e.g. food search learning, existence of alternative food, refuge for prey).     

Linking saturation,stability and sustainability in food webs with observed equilibrium structure

Stability of a dynamic equilibrium in a predator-prey system depends both on the type of functional response and on the point of equilibrium on the response curve. Saturation effects from Holling type II responses are known to destabilise prey populations, while a type III (sigmoid) response curve has been shown to provide stability at lower levels of saturation. These effects have also been shown in multi-trophic model systems. However, stability analyses of observed equilibria in real complex ecosystems have as yet not assumed non-linear functional responses. Here, we evaluate the implications of saturation in observed balanced material-flow structures, for system stability and sustainability. We first make the effects of the non-linear functional responses on the interaction strengths in a food web transparent by expressing the elements of Jacobian ‘community’ matrices for type II and III systems as simple functions of their linear (type I) counterparts. We then determine the stability of the systems and distinguish two critical saturation levels: (1) a level where the system is just as stable as a type I system and (2) a level above which the system cannot be stable unless it is subsidised, separating a stable materially sustainable regime from an unsustainable one. We explain the stabilising and destabilising effects in terms of the feedbacks in the systems. The results shed light on the robustness of observed patterns of interaction strengths in complex food webs and suggest the implausibility of saturation playing a significant role in the equilibrium dynamics of sustainable ecosystems.     

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.     

Body masses,functional responses and predator–prey stability

The stability of ecological communities depends strongly on quantitative characteristics of population interactions (type‐II vs. type‐III functional responses) and the distribution of body masses across species. Until now, these two aspects have almost exclusively been treated separately leaving a substantial gap in our general understanding of food webs. We analysed a large data set of arthropod feeding rates and found that all functional‐response parameters depend on the body masses of predator and prey. Thus, we propose generalised functional responses which predict gradual shifts from type‐II predation of small predators on equally sized prey to type‐III functional‐responses of large predators on small prey. Models including these generalised functional responses predict population dynamics and persistence only depending on predator and prey body masses, and we show that these predictions are strongly supported by empirical data on forest soil food webs. These results help unravelling systematic relationships between quantitative population interactions and large‐scale community patterns.     

Enrichment can damp population cycles: a balance of inflexible and flexible interactions

Destabilization of one predator–one prey systems with an increase in nutrient input has been viewed as a paradox. We report that enrichment can damp population cycles by a food‐web structure that balances inflexible and flexible interaction links (i.e. specialist and generalist predators). We modeled six predator–prey systems involving three or four species in which the predators practice optimal foraging based on prey profitability determined by handling time. In all models, the balance of interaction links simultaneously decreased the amplitude of population oscillations and increased the minimum density with increasing enrichment, leading to a potential theoretical resolution of the paradox of enrichment in non‐equilibrium dynamics. The stabilization mechanism was common to all of the models. Important previous studies on the stability of food webs have also demonstrated that a balance of interaction strengths stabilizes systems, suggesting a general rule of ecosystem stability.     

Consumer-resource dynamics: quantity, quality, and allocation

Background

The dominant paradigm for modeling the complexities of interacting populations and food webs is a system of coupled ordinary differential equations in which the state of each species, population, or functional trophic group is represented by an aggregated numbers-density or biomass-density variable. Here, using the metaphysiological approach to model consumer-resource interactions, we formulate a two-state paradigm that represents each population or group in a food web in terms of both its quantity and quality.

Methodology and Principal Findings

The formulation includes an allocation function controlling the relative proportion of extracted resources to increasing quantity versus elevating quality. Since lower quality individuals senesce more rapidly than higher quality individuals, an optimal allocation proportion exists and we derive an expression for how this proportion depends on population parameters that determine the senescence rate, the per-capita mortality rate, and the effects of these rates on the dynamics of the quality variable. We demonstrate that oscillations do not arise in our model from quantity-quality interactions alone, but require consumer-resource interactions across trophic levels that can be stabilized through judicious resource allocation strategies. Analysis and simulations provide compelling arguments for the necessity of populations to evolve quality-related dynamics in the form of maternal effects, storage or other appropriate structures. They also indicate that resource allocation switching between investments in abundance versus quality provide a powerful mechanism for promoting the stability of consumer-resource interactions in seasonally forcing environments.

Conclusions/Significance

Our simulations show that physiological inefficiencies associated with this switching can be favored by selection due to the diminished exposure of inefficient consumers to strong oscillations associated with the well-known paradox of enrichment. Also our results demonstrate how allocation switching can explain observed growth patterns in experimental microbial cultures and discuss how our formulation can address questions that cannot be answered using the quantity-only paradigms that currently predominate.     

Topological properties of food webs: from real data to community assembly models

We explore patterns of trophic connections between species in the largest and highest-quality empirical food webs to date, introducing a new topological property called the link distribution frequency (i.e. degree distribution), defined as the frequency of species S _L with L links. Non-trivial differences are shown in link distribution frequencies between species-rich and species-poor communities, which might have important consequences for the responses of ecosystems to disturbances. Coarse-grained topological properties observed, as species richness-connectance and number of links-species richness relationships, provide no support for the theory of links-species scaling law or constant connectance across empirical food webs investigated. We further explore these observations by means of simulated food webs resulting from multitrophic assembly models using different functional responses between species. Species richness-connectance and links-species richness relationships of empirical food webs are reproduced by our models, but degree distributions are not properly predicted, suggesting the need of new theoretical approximations to food web assembly. The best agreement between empirical and simulated webs occurs for low values of interaction strength between species, corroborating previous empirical and theoretical findings where weak interactions govern food web dynamics.     

Metabolic adjustment enhances food web stability

Understanding ecosystem stability is one of the greatest challenges of ecology. Over several decades, it has been shown that allometric scaling of biological rates and feeding interactions provide stability to complex food web models. Moreover, introducing adaptive responses of organisms to environmental changes (e.g. like adaptive foraging that enables organisms to adapt their diets depending on resources abundance) improved species persistence in food webs. Here, we introduce the concept of metabolic adjustment, i.e. the ability of species to slow down their metabolic rates when facing starvation and to increase it in time of plenty. We study the reactions of such a model to nutrient enrichment and the adjustment speed of metabolic rates. We found that increasing nutrient enrichment leads to a paradox of enrichment (increase in biomasses and oscillation amplitudes and ultimately extinction of species) but metabolic adjustment stabilises the system by dampening the oscillations. Metabolic adjustment also increases the average biomass of the top predator in a tri‐trophic food chain. In complex food webs, metabolic adjustment has a stabilising effect as it promotes species survival by creating a large diversity of metabolic rates. However, this stabilising effect is mitigated in enriched ecosystems. Phenotypic plasticity of organisms must be considered in food web models to better understand the response of organisms to their environment. As metabolic rate is central in describing biological rates, we must pay attention to its variations to fully understand the population dynamics of natural communities.     

Revisiting the influence of learning in predator functional response,how it can lead to shapes different from type III

Predator/parasitoid functional response is one of the main tools used to study predation behavior, and in assessing the potential of biological control candidates. It is generally accepted that predator learning in prey searching and manipulation can produce the appearance of a type III functional response. Holling proposed that in the presence of alternative prey, at some point the predator would shift the preferred prey, leading to the appearance of a sigmoid function that characterized that functional response. This is supported by the analogy between enzyme kinetics and functional response that Holling used as the basis for developing this theory. However, after several decades, sigmoidal functional responses appear in the absence of alternative prey in most of the biological taxa studied. Here, we propose modeling the effect of learning on the functional response by using the explicit incorporation of learning curves in the parameters of the Holling functional response, the attack rate (a), and the manipulation time (h). We then study how the variation in the parameters of the learning curves causes variations in the shape of the functional response curve. We found that the functional response product of learning can be either type I, II, or III, depending on what parameters act on the organism, and how much it can learn throughout the length of the study. Therefore, the presence of other types of curves should not be automatically associated with the absence of learning. These results are important from an ecological point of view because when type III functional response is associated with learning, it is generally accepted that it can operate as a stabilizing factor in population dynamics. Our results, to the contrary, suggest that depending on how it acts, it may even be destabilizing by generating the appearance of functional responses close to type I.     

Individual behavioral variation in predator–prey models

The role of individual behavioral variation in community dynamics was studied. Behavioral variation in this study does not refer to differences in average responses (e.g., average response between presence and absence of antipredator behavior). Rather it refers to the variation around the average response that is not explained by trivial experimental treatments. First, the effect of behavioral variation was examined based on Jensen’s inequality. In cases of commonly used modeling framework with type II functional response, neglecting behavioral variation (a component of encounter rate) causes overestimation of predation effects. The effect of this bias on community processes was examined by incorporating the behavioral variation in a commonly used consumer-resource model (Rosenzweig–MacArthur model). How such a consideration affects a model prediction (paradox of enrichment) was examined. The inclusion of behavioral variation can both quantitatively and qualitatively alter the model characteristics. Behavioral variation can substantially increase the stability of the community with respect to enrichment.     

Predicting the effects of temperature on food web connectance

Few models concern how environmental variables such as temperature affect community structure. Here, we develop a model of how temperature affects food web connectance, a powerful driver of population dynamics and community structure. We use the Arrhenius equation to add temperature dependence of foraging traits to an existing model of food web structure. The model predicts potentially large temperature effects on connectance. Temperature-sensitive food webs exhibit slopes of up to 0.01 units of connectance per 1°C change in temperature. This corresponds to changes in diet breadth of one resource item per 2°C (assuming a food web containing 50 species). Less sensitive food webs exhibit slopes down to 0.0005, which corresponds to about one resource item per 40°C. Relative sizes of the activation energies of attack rate and handling time determine whether warming increases or decreases connectance. Differences in temperature sensitivity are explained by differences between empirical food webs in the body size distributions of organisms. We conclude that models of temperature effects on community structure and dynamics urgently require considerable development, and also more and better empirical data to parameterize and test them.     

Estimating morphologically determined connectance and structure for food webs of freshwater invertebrates

1. Connectance is a parameter of central importance in determining food-web structure, but the processes determining its value remain unclear. In evaluating possible explanations it is useful to know what patterns, and values, of connectance occur in food webs assembled at random from a set of species in a regional species pool; i.e. where the number of links is determined by the morphological features of the species present, not by the immediate effects of energetics or stability on the particular web. 2. This study examines, by means of laboratory experiments, the occurrence of potential feeding interactions among a set of freshwater invertebrate species randomly selected from different freshwater sites in a geographical region. The results from pairwise feeding trials are used to construct two ‘theoretical’ food webs, in which the patterns and values of connectance are examined. 3. Analyses of these webs indicate that their structure is consistent with the observed values in previously documented ‘real’ webs. Directed connectance values of 0.12–0.16 (or less) suggest that the assembled webs are no more connected than many freshwater webs from natural systems. The number of links per species increases curvilinearly with the number of species, during web assembly, consistent with recent hypotheses. 4. These results also indicate that quantifying, and understanding the determinants of, trophic generalism or specialism does have implications for understanding how connectance is constrained in real webs. Freshwater invertebrates seem to be relatively generalist, and freshwater food webs perhaps correspondingly highly connected. Such arguments have implications for interpreting other aspects of food-web structure in these systems, and for parameterizing models that are based on connectance.     

Relation between complexity and stability in food webs with adaptive behavior

We investigate the influence of functional responses (Lotka-Volterra or Holling type), initial topological web structure (randomly connected or niche model), adaptive behavior (adaptive foraging and predator avoidance) and the type of constraints on the adaptive behavior (linear or nonlinear) on the stability and structure of food webs. Two kinds of stability are considered: one is the network robustness (i.e., the proportion of species surviving after population dynamics) and the other is the species deletion stability. When evaluating the network structure, we consider link density as well as the trophic level structure. We show that the types of functional responses and initial web structure do not have a large effect on the stability of food webs, but foraging behavior has a large stabilizing effect. It leads to a positive complexity-stability relationship whenever higher "complexity" implies more potential prey per species. The other type of adaptive behavior, predator avoidance behavior, makes food webs only slightly more stable. The observed link density after population dynamics depends strongly on the presence or absence of adaptive foraging, and on the type of constraints used. We also show that the trophic level structure is preserved under population dynamics with adaptive foraging.     

Transient dynamics and food-web complexity in the Lotka-Volterra cascade model

How does the long-term behaviour near equilibrium of model food webs correlate with their short-term transient dynamics? Here, simulations of the Lotka -Volterra cascade model of food webs provide the first evidence to answer this question. Transient behaviour is measured by resilience, reactivity, the maximum amplification of a perturbation and the time at which the maximum amplification occurs. Model food webs with a higher probability of local asymptotic stability may be less resilient and may have a larger transient growth of perturbations. Given a fixed connectance, the sizes and durations of transient responses to perturbations increase with the number of species. Given a fixed number of species, as connectance increases, the sizes and durations of transient responses to perturbations may increase or decrease depending on the type of link that is varied. Reactivity is more sensitive to changes in the number of donor-controlled links than to changes in the number of recipient-controlled links, while resilience is more sensitive to changes in the number of recipient-controlled links than to changes in the number of donor-controlled links. Transient behaviour is likely to be one of the important factors affecting the persistence of ecological communities.     

Cannibalism in an age-structured predator-prey system

Recently, Kohlmeier and Ebenhöh showed that cannibalism can stabilize population cycles in a Lotka-Volterra type predator-prey model. Population cycles in their model are due to the interaction between logistic population growth of the prey and a hyperbolic functional response. In this paper, we consider a predator-prey system where cyclic population fluctuations are due to the age structure in the predator species. It is shown that cannibalism is also a stabilizing mechanism when population oscillations are due to this age structure. We conclude that in predator-prey systems, cannibalism by predators can stabilize both externally generated (consumer-resource) as well as internally generated (agestructure) fluctuations.     

Current food web models cannot explain the overall topological structure of observed food webs

Topological food webs illustrating “who eats whom” in different systems exhibit similar, non‐random, structures suggesting that general rules govern food web structure. Current food web models correctly predict many measures of food web topology from knowledge of species richness and connectance (fraction of possible predator–prey links that actually occur), together with assumptions about the ecological rules governing “who eats whom”. However, current measures are relatively insensitive to small changes in topology. Improvement of, and discrimination among, current models requires development of new measures of food web structure. Here I examine whether current food web models (cascade, niche, and nested hierarchy models, plus a random null model) can predict a new measure of food web structure, structural stability. Structural stability complements other measures of food web topology because it is sensitive to changes in topology that other measures often miss. The cascade and null models respectively over‐ and underpredict structural stability for a set of 17 high‐quality food webs. While the niche and nested hierarchy models provide unbiased predictions on average, their 95% confidence intervals frequently fail to include the observed data. Observed structural stabilities for all models are overdispersed compared to model predictions, and predicted and observed structural stabilities are uncorrelated, indicating that important sources of variation in structural stability are not captured by the models. Crucially, poor model performance arises because observed variation in structural stability is unrelated to variation in species richness and connectance. In contrast, almost all other measures of food web topology vary with species richness and connectance in natural webs. No model that takes species richness and connectance as the only input parameters can reproduce observed variation in structural stability. Further progress in predicting and explaining food web topology will require fundamentally new models based on different input parameters.     

Does structural sensitivity alter complexity–stability relationships?

Structural sensitivity, namely the sensitivity of a model dynamics to slight changes in its mathematical formulation, has already been studied in some models with a small number of state variables. The aim of this study is to investigate the impact of structural sensitivity in a food web model. Especially, the importance of structural sensitivity is compared to that of trophic complexity (number of species, connectance), which is known to strongly influence food web dynamics. Food web structures are built using the niche model. Then food web dynamics are modeled using several type II functional responses parameterized to fit the same predation fluxes. Food web persistence was found to be mostly determined by trophic complexity. At the opposite, even if food web connectance promotes equilibrium dynamics, their occurrence is mainly driven by the choice of the functional response. These conclusions are robust to changes in some parameter values, the fitting method and some model assumptions. In a one-prey/one-predator system, it was shown that the possibility that multiple stable states coexist can be highly structural sensitive. Quantifying this type of uncertainty at the scale of ecosystem models will be both a natural extension to this work and a challenging issue.     

Remarks on antipredator behavior and food chain dynamics

When consumers feeding on a resource spend time in avoiding high risks of predation, the predator functional response declines with predator density. While this is well established, less attention has been paid to the dependence of the consumer functional response on predator density. Here we show how the separation of behavioral and ecological timescales allows one to determine both responses starting from an explicit behavioral model. Within the general set-up considered in this paper, the two functional responses can tend toward Holling type II responses when consumers react only weakly to predation. Thus, the main characteristics of the standard Rosenzweig-MacArthur tritrophic food chain (logistic resource and Holling type II consumer and predator) remain valid also when consumers have weak antipredator behavior. Moreover, through numerical analysis, we show that in a particular but interesting case pronounced antipredator behaviors stabilize the system.