An official website of the United States government
We apply a perturbative Doi-Peliti field-theoretical analysis to the stochastic spatially extended symmetric Rock-paper-Scissors (RPS) and May-Leonard (ML) models, in which three species compete cyclically. Compared to the two-species Lotka-Volterra predator-prey (LV) model, according to numerical simulations, these cyclical models appear to be less affected by intrinsic stochastic fluctuations. Indeed, we demonstrate that the qualitative features of the ML model are insensitive to intrinsic reaction noise. In contrast, and although not yet observed in numerical simulations, we find that the RPS model acquires significant fluctuation-induced renormalizations in the perturbative regime, similar to the LV model. We also study the formation of spatio-temporal structures in the framework of stability analysis and provide a clearcut explanation for the absence of spatial patterns in the RPS system, whereas the spontaneous emergence of spatio-temporal structures features prominently in the LV and the ML models.
more »
« less
Lotka-Volterra predator-prey model with periodically varying carrying capacity
We study the stochastic spatial Lotka-Volterra model for predator-prey interaction subject to a periodically varying carrying capacity. The Lotka-Volterra model with on-site lattice occupation restrictions (i.e., finite local carrying capacity) that represent finite food resources for the prey population exhibits a continuous active-to-absorbing phase transition. The active phase is sustained by the existence of spatio-temporal patterns in the form of pursuit and evasion waves. Monte Carlo simulations on a two-dimensional lattice are utilized to investigate the effect of seasonal variations of the environment on species coexistence. The results of our simulations are also compared to a mean-field analysis in order to specifically delineate the impact of stochastic fluctuations and spatial correlations. We find that the parameter region of predator and prey coexistence is enlarged relative to the stationary situation when the carrying capacity varies periodically. The (quasi-)stationary regime of our periodically varying Lotka-Volterra predator-prey system shows qualitative agreement between the stochastic model and the mean-field approximation. However, under periodic carrying capacity-switching environments, the mean-field rate equations predict period-doubling scenarios that are washed out by internal reaction noise in the stochastic lattice model. Utilizing visual representations of the lattice simulations and dynamical correlation functions, we study how the pursuit and evasion waves are affected by ensuing resonance effects. Correlation function measurements indicate a time delay in the response of the system to sudden changes in the environment. Resonance features are observed in our simulations that cause prolonged persistent spatial correlations. Different effective static environments are explored in the extreme limits of fast and slow periodic switching. The analysis of the mean-field equations in the fast-switching regime enables a semi-quantitative description of the (quasi-)stationary state.
more »
« less
- Award ID(s):
- 2128587
- PAR ID:
- 10505692
- Publisher / Repository:
- The American Physical Society
- Date Published:
- Journal Name:
- Physical Review E
- Volume:
- 107
- Issue:
- 6
- ISSN:
- 2470-0045
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Stochastic reaction-diffusion models are employed to represent many complex physical, biological, societal, and ecological systems. The macroscopic reaction rates describing the large-scale, long-time kinetics in such systems are effective, scale-dependent renormalized parameters that need to be either measured experimentally or computed by means of a microscopic model. In a Monte Carlo simulation of stochastic reaction-diffusion systems, microscopic probabilities for specific events to happen serve as the input control parameters. To match the results of any computer simulation to observations or experiments carried out on the macroscale, a mapping is required between the microscopic probabilities that define the Monte Carlo algorithm and the macroscopic reaction rates that are experimentally measured. Finding the functional dependence of emergent macroscopic rates on the microscopic probabilities (subject to specific rules of interaction) is a very difficult problem, and there is currently no systematic, accurate analytical way to achieve this goal. Therefore, we introduce a straightforward numerical method of using lattice Monte Carlo simulations to evaluate the macroscopic reaction rates by directly obtaining the count statistics of how many events occur per simulation time step. Our technique is first tested on well-understood fundamental examples, namely, restricted birth processes, diffusion-limited two-particle coagulation, and two-species pair annihilation kinetics. Next we utilize the thus gained experience to investigate how the microscopic algorithmic probabilities become coarse-grained into effective macroscopic rates in more complex model systems such as the Lotka-Volterra model for predator-prey competition and coexistence, as well as the rock-paper-scissors or cyclic Lotka-Volterra model and its May-Leonard variant that capture population dynamics with cyclic dominance motifs. Thereby we achieve a more thorough and deeper understanding of coarse graining in spatially extended stochastic reaction-diffusion systems and the nontrivial relationships between the associated microscopic and macroscopic model parameters, with a focus on ecological systems. The proposed technique should generally provide a useful means to better fit Monte Carlo simulation results to experimental or observational data.more » « less
-
For species primarily regulated by a common predator, the P * rule of Holt & Lawton (Holt & Lawton, 1993. Am. Nat. 142 , 623–645. ( doi:10.1086/285561 )) predicts that the prey species that supports the highest mean predator density ( P *) excludes the other prey species. This prediction is re-examined in the presence of temporal fluctuations in the vital rates of the interacting species including predator attack rates. When the fluctuations in predator attack rates are temporally uncorrelated, the P * rule still holds even when the other vital rates are temporally auto-correlated. However, when temporal auto-correlations in attack rates are positive but not too strong, the prey species can coexist due to the emergence of a positive covariance between predator density and prey vulnerability. This coexistence mechanism is similar to the storage effect for species regulated by a common resource. Negative or strongly positive auto-correlations in attack rates generate a negative covariance between predator density and prey vulnerability and a stochastic priority effect can emerge: with non-zero probability either prey species is excluded. These results highlight how temporally auto-correlated species’ interaction rates impact the structure and dynamics of ecological communities.more » « less
-
Abstract Despite the well known scale‐dependency of ecological interactions, relatively little attention has been paid to understanding the dynamic interplay between various spatial scales. This is especially notable in metacommunity theory, where births and deaths dominate dynamics within patches (the local scale), and dispersal and environmental stochasticity dominate dynamics between patches (the regional scale). By considering the interplay of local and regional scales in metacommunities, the fundamental processes of community ecology—selection, drift, and dispersal—can be unified into a single theoretical framework. Here, we analyze three related spatial models that build on the classic two‐species Lotka–Volterra competition model. Two open‐system models focus on a single patch coupled to a larger fixed landscape by dispersal. The first is deterministic, while the second adds demographic stochasticity to allow ecological drift. Finally, the third model is a true metacommunity model with dispersal between a large number of local patches, which allows feedback between local and regional scales and captures the well studied metacommunity paradigms as special cases. Unlike previous simulation models, our metacommunity model allows the numerical calculation of equilibria and invasion criteria to precisely determine the outcome of competition at the regional scale. We show that both dispersal and stochasticity can lead to regional outcomes that are different than predicted by the classic Lotka–Volterra competition model. Regional exclusion can occur when the nonspatial model predicts coexistence or founder control, due to ecological drift or asymmetric stochastic switching between basins of attraction, respectively. Regional coexistence can result from local coexistence mechanisms or through competition‐colonization or successional‐niche trade‐offs. Larger dispersal rates are typically competitively advantageous, except in the case of local founder control, which can favor intermediate dispersal rates. Broadly, our models demonstrate the importance of feedback between local and regional scales in competitive metacommunities and provide a unifying framework for understanding how selection, drift, and dispersal jointly shape ecological communities.more » « less
-
Zhu, Shanfeng (Ed.)Understanding microbial interactions is fundamental for exploring population dynamics, particularly in microbial communities where interactions affect stability and host health. Generalized Lotka-Volterra (gLV) models have been widely used to investigate system dynamics but depend on absolute abundance data, which are often unavailable in microbiome studies. To address this limitation, we introduce an iterative Lotka-Volterra (iLV) model, a novel framework tailored for compositional data that leverages relative abundances and iterative refinements for parameter estimation. The iLV model features two key innovations: an adaptation of the gLV framework to compositional constraints and an iterative optimization strategy combining linear approximations with nonlinear refinements to enhance parameter estimation accuracy. Using simulations and real-world datasets, we demonstrate that iLV surpasses existing methodologies, such as the compositional LV (cLV) and the generalized LV (gLV) model, in recovering interaction coefficients and predicting species trajectories under varying noise levels and temporal resolutions. Applications to the lynx-hare predator-prey,Stylonychia pustula-P. caudatummixed culture, and cheese microbial systems revealed consistency between predicted and observed relative abundances showcasing its accuracy and robustness. In summary, the iLV model bridges theoretical gLV models and practical compositional data analysis, offering a robust framework to infer microbial interactions and predict community dynamics using relative abundance data, with significant potential for advancing microbial research.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1103/PhysRevE.107.064144
