Refine
Year of publication
Document Type
- Article (106)
- Postprint (10)
- Monograph/Edited Volume (2)
- Conference Proceeding (2)
- Part of a Book (1)
- Part of Periodical (1)
- Review (1)
Keywords
- functional traits (7)
- coexistence (6)
- community ecology (6)
- population dynamics (6)
- lake (5)
- plankton (5)
- predator-prey dynamics (5)
- top-down control (5)
- biodiversity (4)
- functional diversity (4)
Institute
- Institut für Biochemie und Biologie (103)
- Institut für Geowissenschaften (10)
- Institut für Physik und Astronomie (3)
- Mathematisch-Naturwissenschaftliche Fakultät (3)
- Extern (2)
- Zentrum für Umweltwissenschaften (2)
- Dezernat 2: Studienangelegenheiten (1)
- Institut für Umweltwissenschaften und Geographie (1)
Phenotypic plasticity can increase individual fitness when environmental conditions change over time. Inducible defences are a striking example, allowing species to react to fluctuating predation pressure by only expressing their costly defended phenotype under high predation risk. Previous theoretical investigations have focused on how this affects predator–prey dynamics, but the impact on competitive outcomes and broader community dynamics has received less attention. Here we use a small food web model, consisting of two competing plastic autotrophic species exploited by a shared consumer, to study how the speed of inducible defences across three trade-off constellations affects autotroph coexistence, biomasses across trophic levels, and temporal variability. Contrary to the intuitive idea that faster adaptation increases autotroph fitness, we found that higher switching rates reduced individual fitness as it consistently provoked more maladaptive switching towards undefended phenotypes under high predation pressure. This had an unexpected positive impact on the consumer, increasing consumer biomass and lowering total autotroph biomass. Additionally, maladaptive switching strongly reduced autotroph coexistence through an emerging source-sink dynamic between defended and undefended phenotypes. The striking impact of maladaptive switching on species and food web dynamics indicates that this mechanism may be of more critical importance than previously recognized.
Phenotypic plasticity can increase individual fitness when environmental conditions change over time. Inducible defences are a striking example, allowing species to react to fluctuating predation pressure by only expressing their costly defended phenotype under high predation risk. Previous theoretical investigations have focused on how this affects predator–prey dynamics, but the impact on competitive outcomes and broader community dynamics has received less attention. Here we use a small food web model, consisting of two competing plastic autotrophic species exploited by a shared consumer, to study how the speed of inducible defences across three trade-off constellations affects autotroph coexistence, biomasses across trophic levels, and temporal variability. Contrary to the intuitive idea that faster adaptation increases autotroph fitness, we found that higher switching rates reduced individual fitness as it consistently provoked more maladaptive switching towards undefended phenotypes under high predation pressure. This had an unexpected positive impact on the consumer, increasing consumer biomass and lowering total autotroph biomass. Additionally, maladaptive switching strongly reduced autotroph coexistence through an emerging source-sink dynamic between defended and undefended phenotypes. The striking impact of maladaptive switching on species and food web dynamics indicates that this mechanism may be of more critical importance than previously recognized.
Quantifying the capacity for contemporary trait changes to drive intermittent predator-prey cycles
(2022)
A large and growing body of theory has demonstrated how the presence of trait variation in prey or predator populations may affect the amplitude and phase of predator-prey cycles. Less attention has been given to so-called intermittent cycles, in which predator-prey oscillations recurrently disappear and re-appear, despite such dynamics being observed in empirical systems and modeling studies. A comprehensive understanding of the conditions under which trait changes may drive intermittent predator-prey dynamics, as well as their potential ecological implications, is therefore missing. Here we provide a first systematic analysis of the eco-evolutionary conditions that may give rise to intermittent predator-prey cycles, investigating 16 models that incorporate different types of trait variation within prey, predators, or both. Our results show that intermittent dynamics often arise through predator-prey coevolution, but only very rarely when only one trophic level can adapt. Additionally, the frequency of intermittent cycles depends on the source of trait variation (genetic variation or phenotypic plasticity) and the genetic architecture (Mendelian or quantitative traits), with intermittency occurring most commonly through Mendelian evolution, and very rarely through phenotypic plasticity. Further analysis identified three necessary conditions for when trait variation can drive intermittent cycles. First, the intrinsic stability of the predator-prey system must depend on the traits of prey, predators, or both. Second, there must be a mechanism causing the recurrent alternation between stable and unstable states, leading to a "trait" cycle superimposed on the population dynamics. Finally, these trait dynamics must be significantly slower than the predator-prey cycles. We show how these conditions explain all the abovementioned patterns. We further show an important unexpected consequence of these necessary conditions: they are most easily met when intraspecific trait variation is at high risk of being lost. As trait diversity is positively associated with ecosystem functioning, this can have potentially severe negative consequences. This novel result highlights the importance of identifying and understanding intermittent cycles in theoretical studies and natural systems. The new approach for detecting and quantifying intermittency we develop here will be instrumental in enabling future study.
Trophic transfer efficiency (TTE) is usually calculated as the ratio of production rates between two consecutive trophic levels. Although seemingly simple, TTE estimates from lakes are rare. In our review, we explore the processes and structures that must be understood for a proper lake TTE estimate.
We briefly discuss measurements of production rates and trophic positions and mention how ecological efficiencies, nutrients (N, P) and other compounds (fatty acids) affect energy transfer between trophic levels and hence TTE.
Furthermore, we elucidate how TTE estimates are linked with size-based approaches according to the Metabolic Theory of Ecology, and how food-web models can be applied to study TTE in lakes.
Subsequently, we explore temporal and spatial heterogeneity of production and TTE in lakes, with a particular focus on the links between benthic and pelagic habitats and between the lake and the terrestrial environment.
We provide an overview of TTE estimates from lakes found in the published literature. Finally, we present two alternative approaches to estimating TTE. First, TTE can be seen as a mechanistic quantity informing about the energy and matter flow between producer and consumer groups.
This approach is informative with respect to food-web structure, but requires enormous amounts of data. The greatest uncertainty comes from the proper consideration of basal production to estimate TTE of omnivorous organisms.
An alternative approach is estimating food-chain and food-web efficiencies, by comparing the heterotrophic production of single consumer levels or the total sum of all heterotrophic production including that of heterotrophic bacteria to the total sum of primary production.
We close the review by pointing to a few research questions that would benefit from more frequent and standardized estimates of TTE in lakes.
Ciliates represent a crucial link between phytoplankton and bacteria and mesozooplankton in pelagic food webs, but little is known about the processes influencing the dynamics of individual species.
Using long-term, high-frequency observations, we compared the diversity and the temporal variability in biomass and species composition of the ciliate community in large, deep, mesotrophic Lake Constance to that of the phytoplankton and rotifer communities in the same lake.
Furthermore, we used boosted regression trees to evaluate possible environmental predictors (temperature, three prey groups, four predator/competitor groups) influencing ciliate net growth.
The biomass of all ciliate species showed a common, recurrent seasonal pattern, often with peaks in spring and summer.
The ciliate community was more diverse than the rotifer community, exhibited highly synchronous dynamics and its species were regularly encountered during the season. The top-down control by copepods likely contributes to the ciliates' synchronized decline prior to the clear-water phase when food concentration is still high.
The high temporal autocorrelation of the ciliate biomasses together with the inter-annual recurrent seasonal patterns and the low explanatory power of the environmental predictors suggest that the dynamics of individual ciliate species are strictly controlled, yet it remains difficult to determine the responsible factors.
It is well known that functional diversity strongly affects ecosystem functioning. However, even in rather simple model communities consisting of only two or, at best, three trophic levels, the relationship between multitrophic functional diversity and ecosystem functioning appears difficult to generalize, because of its high contextuality. In this study, we considered several differently structured tritrophic food webs, in which the amount of functional diversity was varied independently on each trophic level. To achieve generalizable results, largely independent of parametrization, we examined the outcomes of 128,000 parameter combinations sampled from ecologically plausible intervals, with each tested for 200 randomly sampled initial conditions. Analysis of our data was done by training a random forest model. This method enables the identification of complex patterns in the data through partial dependence graphs, and the comparison of the relative influence of model parameters, including the degree of diversity, on food-web properties. We found that bottom-up and top-down effects cascade simultaneously throughout the food web, intimately linking the effects of functional diversity of any trophic level to the amount of diversity of other trophic levels, which may explain the difficulty in unifying results from previous studies. Strikingly, only with high diversity throughout the whole food web, different interactions synergize to ensure efficient exploitation of the available nutrients and efficient biomass transfer to higher trophic levels, ultimately leading to a high biomass and production on the top level. The temporal variation of biomass showed a more complex pattern with increasing multitrophic diversity: while the system initially became less variable, eventually the temporal variation rose again because of the increasingly complex dynamical patterns. Importantly, top predator diversity and food-web parameters affecting the top trophic level were of highest importance to determine the biomass and temporal variability of any trophic level. Overall, our study reveals that the mechanisms by which diversity influences ecosystem functioning are affected by every part of the food web, hampering the extrapolation of insights from simple monotrophic or bitrophic systems to complex natural food webs.
Trait variation among heterospecific and conspecific organisms may substantially affect community and food web dynamics. While the relevance of competition and feeding traits have been widely studied for different consumer species, studies on intraspecific differences are more scarce, partly owing to difficulties in distinguishing different clones of the same species. Here, we investigate how intraspecific trait variation affects the competition between the freshwater ciliates Euplotes octocarinatus and Coleps hirtus in a nitrogen-limited chemostat system. The ciliates competed for the microalgae Cryptomonas sp. (Cry) and Navicula pelliculosa (Nav), and the bacteria present in the cultures over a period of 33 days. We used monoclonal Euplotes and three different Coleps clones (Col 1, Col 2, and Col 3) in the experiment that could be distinguished by a newly developed rDNA-based molecular assay based on the internal transcribed spacer (ITS) regions. While Euplotes feeds on Cry and on bacteria, the Coleps clones cannot survive on bacteria alone but feed on both Cry and Nav with clone-specific rates. Experimental treatments comprised two-species mixtures of Euplotes and one or all of the three different Coleps clones, respectively. We found intraspecific variation in the traits "selectivity" and "maximum ingestion rate" for the different algae to significantly affect the competitive outcome between the two ciliate species. As Nav quickly escaped top-down control and likely reached a state of low food quality, ciliate competition was strongly determined by the preference of different Coleps clones for Cry as opposed to feeding on Nav. In addition, the ability of Euplotes to use bacteria as an alternative food source strengthened its persistence once Cry was depleted. Hence, trait variation at both trophic levels codetermined the population dynamics and the outcome of species competition.
Seagrass beds are important habitats in coastal areas but increasingly decline in area and quality, thus conservation measures are urgently needed. Quantitative food webs, describing the biomass distribution and energy fluxes among trophic groups, reveal structural and functional aspects of ecosystems. Their knowledge can improve ecological conservation. For the recently discovered large warm-temperate seagrass (Zostera japonica) habitat in China's Yellow River Delta wetland, we used delta C-13 and delta N-15 measurements and a Bayesian isotope mixing model to construct its food web diagram with quantitative estimations of consumer diet compositions, comprising detritus and 14 living trophic groups from primary producers to fish. We then estimated the quantitative food web fluxes based on biomass measurements and calculated corresponding ecosystem functions. Pelagic producers were significantly C-13-depleted compared to benthic sources. Consumers (except zooplankton) were increasingly C-13-depleted with increasing trophic positions even though the consumed benthic production surpassed the pelagic one. Bivalves dominated consumer biomasses and fluxes and were the first to connect the pelagic and benthic pathways, whereas zooplankton and gastropods were specialized on the two pathways, respectively. We found flat biomass and production pyramids indicating low trophic transfer efficiencies. Generally, the energetic structure of the quantitative food web was consistent with the stable isotope analysis, and the estimated net primary production and most estimated production to biomass ratios of the trophic groups fell within literature ranges. This study provides a systematical understanding of the quantitative trophic ecology of a seagrass bed and facilitates synergistic knowledge on management, conservation, and restoration.
Biodiversity decline causes a loss of functional diversity, which threatens ecosystems through a dangerous feedback loop: This loss may hamper ecosystems’ ability to buffer environmental changes, leading to further biodiversity losses. In this context, the increasing frequency of human-induced excessive loading of nutrients causes major problems in aquatic systems. Previous studies investigating how functional diversity influences the response of food webs to disturbances have mainly considered systems with at most two functionally diverse trophic levels. We investigated the effects of functional diversity on the robustness, that is, resistance, resilience, and elasticity, using a tritrophic—and thus more realistic—plankton food web model. We compared a non-adaptive food chain with no diversity within the individual trophic levels to a more diverse food web with three adaptive trophic levels. The species fitness differences were balanced through trade-offs between defense/growth rate for prey and selectivity/half-saturation constant for predators. We showed that the resistance, resilience, and elasticity of tritrophic food webs decreased with larger perturbation sizes and depended on the state of the system when the perturbation occurred. Importantly, we found that a more diverse food web was generally more resistant and resilient but its elasticity was context-dependent. Particularly, functional diversity reduced the probability of a regime shift toward a non-desirable alternative state. The basal-intermediate interaction consistently determined the robustness against a nutrient pulse despite the complex influence of the shape and type of the dynamical attractors. This relationship was strongly influenced by the diversity present and the third trophic level. Overall, using a food web model of realistic complexity, this study confirms the destructive potential of the positive feedback loop between biodiversity loss and robustness, by uncovering mechanisms leading to a decrease in resistance, resilience, and potentially elasticity as functional diversity declines.
Biodiversity decline causes a loss of functional diversity, which threatens ecosystems through a dangerous feedback loop: This loss may hamper ecosystems’ ability to buffer environmental changes, leading to further biodiversity losses. In this context, the increasing frequency of human-induced excessive loading of nutrients causes major problems in aquatic systems. Previous studies investigating how functional diversity influences the response of food webs to disturbances have mainly considered systems with at most two functionally diverse trophic levels. We investigated the effects of functional diversity on the robustness, that is, resistance, resilience, and elasticity, using a tritrophic—and thus more realistic—plankton food web model. We compared a non-adaptive food chain with no diversity within the individual trophic levels to a more diverse food web with three adaptive trophic levels. The species fitness differences were balanced through trade-offs between defense/growth rate for prey and selectivity/half-saturation constant for predators. We showed that the resistance, resilience, and elasticity of tritrophic food webs decreased with larger perturbation sizes and depended on the state of the system when the perturbation occurred. Importantly, we found that a more diverse food web was generally more resistant and resilient but its elasticity was context-dependent. Particularly, functional diversity reduced the probability of a regime shift toward a non-desirable alternative state. The basal-intermediate interaction consistently determined the robustness against a nutrient pulse despite the complex influence of the shape and type of the dynamical attractors. This relationship was strongly influenced by the diversity present and the third trophic level. Overall, using a food web model of realistic complexity, this study confirms the destructive potential of the positive feedback loop between biodiversity loss and robustness, by uncovering mechanisms leading to a decrease in resistance, resilience, and potentially elasticity as functional diversity declines.
Trait-based approaches have broadened our understanding of how the composition of ecological communities responds to environmental drivers. This research has mainly focussed on abiotic factors and competition determining the community trait distribution, while effects of trophic interactions on trait dynamics, if considered at all, have been studied for two trophic levels at maximum. However, natural food webs are typically at least tritrophic. This enables indirect interactions of traits and biomasses among multiple trophic levels leading to underexplored effects on food web dynamics. Here, we demonstrate the occurrence of mutual trait adjustment among three trophic levels in a natural plankton food web (Lake Constance) and in a corresponding mathematical model. We found highly recurrent seasonal biomass and trait dynamics, where herbivorous zooplankton increased its size, and thus its ability to counter phytoplankton defense, before phytoplankton defense actually increased. This is contrary to predictions from bitrophic systems where counter-defense of the consumer is a reaction to prey defense. In contrast, counter-defense of carnivores by size adjustment followed the defense of herbivores as expected. By combining observations and model simulations, we show how the reversed trait dynamics at the two lower trophic levels result from a "trophic biomass-trait cascade" driven by the carnivores. Trait adjustment between two trophic levels can therefore be altered by biomass or trait changes of adjacent trophic levels. Hence, analyses of only pairwise trait adjustment can be misleading in natural food webs, while multitrophic trait-based approaches capture indirect biomass-trait interactions among multiple trophic levels.
The shape of a defense-growth trade-off governs seasonal trait dynamics in natural phytoplankton
(2020)
Theory predicts that trade-offs, quantifying costs of functional trait adjustments, crucially affect community trait adaptation to altered environmental conditions, but empirical verification is scarce. We evaluated trait dynamics (antipredator defense, maximum growth rate, and phosphate affinity) of a lake phytoplankton community in a seasonally changing environment, using literature trait data and 21 years of species-resolved high-frequency biomass measurements. The trait data indicated a concave defense-growth trade-off, promoting fast-growing species with intermediate defense. With seasonally increasing grazing pressure, the community shifted toward higher defense levels at the cost of lower growth rates along the trade-off curve, while phosphate affinity explained some deviations from it. We discuss how low fitness differences of species, inferred from model simulations, in concert with stabilizing mechanisms, e.g., arising from further trait dimensions, may lead to the observed phytoplankton diversity. In conclusion, quantifying trade-offs is key for predictions of community trait adaptation and biodiversity under environmental change.
The shape of a defense-growth trade-off governs seasonal trait dynamics in natural phytoplankton
(2020)
Theory predicts that trade-offs, quantifying costs of functional trait adjustments, crucially affect community trait adaptation to altered environmental conditions, but empirical verification is scarce. We evaluated trait dynamics (antipredator defense, maximum growth rate, and phosphate affinity) of a lake phytoplankton community in a seasonally changing environment, using literature trait data and 21 years of species-resolved high-frequency biomass measurements. The trait data indicated a concave defense-growth trade-off, promoting fast-growing species with intermediate defense. With seasonally increasing grazing pressure, the community shifted toward higher defense levels at the cost of lower growth rates along the trade-off curve, while phosphate affinity explained some deviations from it. We discuss how low fitness differences of species, inferred from model simulations, in concert with stabilizing mechanisms, e.g., arising from further trait dimensions, may lead to the observed phytoplankton diversity. In conclusion, quantifying trade-offs is key for predictions of community trait adaptation and biodiversity under environmental change.
Parasites, such as bacterial viruses (phages), can have large effects on host populations both at the ecological and evolutionary levels. In the case of cyanobacteria, phages can reduce primary production and infected hosts release intracellular nutrients influencing planktonic food web structure, community dynamics, and biogeochemical cycles. Cyanophages may be of great importance in aquatic food webs during large cyanobacterial blooms unless the host population becomes resistant to phage infection. The consequences on plankton community dynamics of the evolution of phage resistance in bloom forming cyanobacterial populations are still poorly studied. Here, we examined the effect of different frequencies of a phage-resistant genotype within a filamentous nitrogen-fixing Nodularia spumigena population on an experimental plankton community. Three Nodularia populations with different initial frequencies (0%, 5%, and 50%) of phage-resistant genotypes were inoculated in separate treatments with the phage 2AV2, the green alga Chlorella vulgaris, and the rotifer Brachionus plicatilis, which formed the experimental plankton community subjected to either nitrogen-limited or nitrogen-rich conditions. We found that the frequency of the phage-resistant Nodularia genotype determined experimental community dynamics. Cyanobacterial populations with a high frequency (50%) of the phage-resistant genotype dominated the cultures despite the presence of phages, retaining most of the intracellular nitrogen in the plankton community. In contrast, populations with low frequencies (0% and 5%) of the phage-resistant genotype were lysed and reduced to extinction by the phage, transferring the intracellular nitrogen held by Nodularia to Chlorella and rotifers, and allowing Chlorella to dominate the communities and rotifers to survive. This study shows that even though phages represent minuscule biomass, they can have key effects on community composition and eco-evolutionary feedbacks in plankton communities.
Estimating parameters from multiple time series of population dynamics using bayesian inference
(2019)
Empirical time series of interacting entities, e.g., species abundances, are highly useful to study ecological mechanisms. Mathematical models are valuable tools to further elucidate those mechanisms and underlying processes. However, obtaining an agreement between model predictions and experimental observations remains a demanding task. As models always abstract from reality one parameter often summarizes several properties. Parameter measurements are performed in additional experiments independent of the ones delivering the time series. Transferring these parameter values to different settings may result in incorrect parametrizations. On top of that, the properties of organisms and thus the respective parameter values may vary considerably. These issues limit the use of a priori model parametrizations. In this study, we present a method suited for a direct estimation of model parameters and their variability from experimental time series data. We combine numerical simulations of a continuous-time dynamical population model with Bayesian inference, using a hierarchical framework that allows for variability of individual parameters. The method is applied to a comprehensive set of time series from a laboratory predator-prey system that features both steady states and cyclic population dynamics. Our model predictions are able to reproduce both steady states and cyclic dynamics of the data. Additionally to the direct estimates of the parameter values, the Bayesian approach also provides their uncertainties. We found that fitting cyclic population dynamics, which contain more information on the process rates than steady states, yields more precise parameter estimates. We detected significant variability among parameters of different time series and identified the variation in the maximum growth rate of the prey as a source for the transition from steady states to cyclic dynamics. By lending more flexibility to the model, our approach facilitates parametrizations and shows more easily which patterns in time series can be explained also by simple models. Applying Bayesian inference and dynamical population models in conjunction may help to quantify the profound variability in organismal properties in nature.
The size structure of autotroph communities - the relative abundance of small vs. large individuals - shapes the functioning of ecosystems. Whether common mechanisms underpin the size structure of unicellular and multicellular autotrophs is, however, unknown. Using a global data compilation, we show that individual body masses in tree and phytoplankton communities follow power-law distributions and that the average exponents of these individual size distributions (ISD) differ. Phytoplankton communities are characterized by an average ISD exponent consistent with three-quarter-power scaling of metabolism with body mass and equivalence in energy use among mass classes. Tree communities deviate from this pattern in a manner consistent with equivalence in energy use among diameter size classes. Our findings suggest that whilst universal metabolic constraints ultimately underlie the emergent size structure of autotroph communities, divergent aspects of body size (volumetric vs. linear dimensions) shape the ecological outcome of metabolic scaling in forest vs. pelagic ecosystems.
In ecological communities, especially the pelagic zones of aquatic ecosystems, certain bodysize ranges are often over-represented compared to others. Community size spectra, the distributions of community biomass over the logarithmic body-mass axis, tend to exhibit regularly spaced local maxima, called "domes", separated by steep troughs. Contrasting established theory, we explain these dome patterns as manifestations of top-down trophic cascades along aquatic food chains. Compiling high quality size-spectrum data and comparing these with a size-spectrum model introduced in this study, we test this theory and develop a detailed picture of the mechanisms by which bottom-up and top-down effects interact to generate dome patterns. Results imply that strong top-down trophic cascades are common in freshwater communities, much more than hitherto demonstrated, and may arise in nutrient rich marine systems as well. Transferring insights from the general theory of nonlinear pattern formation to domes patterns, we provide new interpretations of past lake-manipulation experiments.
Diverse communities can adjust their trait composition to altered environmental conditions, which may strongly influence their dynamics. Previous studies of trait-based models mainly considered only one or two trophic levels, whereas most natural system are at least tritrophic. Therefore, we investigated how the addition of trait variation to each trophic level influences population and community dynamics in a tritrophic model. Examining the phase relationships between species of adjacent trophic levels informs about the strength of top-down or bottom-up control in non-steadystate situations. Phase relationships within a trophic level highlight compensatory dynamical patterns between functionally different species, which are responsible for dampening the community temporal variability. Furthermore, even without trait variation, our tritrophic model always exhibits regions with two alternative states with either weak or strong nutrient exploitation, and correspondingly low or high biomass production at the top level. However, adding trait variation increased the basin of attraction of the high-production state, and decreased the likelihood of a critical transition from the high- to the lowproduction state with no apparent early warning signals. Hence, our study shows that trait variation enhances resource use efficiency, production, stability, and resilience of entire food webs.
The intrinsic predictability of ecological time series and its potential to guide forecasting
(2019)
Predator-prey cycles rank among the most fundamental concepts in ecology, are predicted by the simplest ecological models and enable, theoretically, the indefinite persistence of predator and prey(1-4). However, it remains an open question for how long cyclic dynamics can be self-sustained in real communities. Field observations have been restricted to a few cycle periods(5-8) and experimental studies indicate that oscillations may be short-lived without external stabilizing factors(9-19). Here we performed microcosm experiments with a planktonic predator-prey system and repeatedly observed oscillatory time series of unprecedented length that persisted for up to around 50 cycles or approximately 300 predator generations. The dominant type of dynamics was characterized by regular, coherent oscillations with a nearly constant predator-prey phase difference. Despite constant experimental conditions, we also observed shorter episodes of irregular, non-coherent oscillations without any significant phase relationship. However, the predator-prey system showed a strong tendency to return to the dominant dynamical regime with a defined phase relationship. A mathematical model suggests that stochasticity is probably responsible for the reversible shift from coherent to non-coherent oscillations, a notion that was supported by experiments with external forcing by pulsed nutrient supply. Our findings empirically demonstrate the potential for infinite persistence of predator and prey populations in a cyclic dynamic regime that shows resilience in the presence of stochastic events.