We investigate the competing effects and relative importance of intrinsic demographic and environmental variability on the evolutionary dynamics of a stochastic two-species Lotka-Volterra model by means of Monte Carlo simulations on a two-dimensional lattice. Individuals are assigned inheritable predation efficiencies; quenched randomness in the spatially varying reaction rates serves as environmental noise. We find that environmental variability enhances the population densities of both predators and prey while demographic variability leads to essentially neutral optimization.
We perform individual-based Monte Carlo simulations in a community consisting of two predator species competing for a single prey species, with the purpose of studying biodiversity stabilization in this simple model system. Predators are characterized with predation efficiency and death rates, to which Darwinian evolutionary adaptation is introduced. Competition for limited prey abundance drives the populations optimization with respect to predation efficiency and death rates. We study the influence of various ecological elements on the final state, finding that both indirect competition and evolutionary adaptation are insufficient to yield a stable ecosystem. However, stable three-species coexistence is observed when direct interaction between the two predator species is implemented.
It is well-established that including spatial structure and stochastic noise in models for predator-prey interactions invalidates the classical deterministic Lotka-Volterra picture of neutral population cycles. In contrast, stochastic models yield long-lived, but ultimately decaying erratic population oscillations, which can be understood through a resonant amplification mechanism for density fluctuations. In Monte Carlo simulations of spatial stochastic predator-prey systems, one observes striking complex spatio-temporal structures. These spreading activity fronts induce persistent correlations between predators and prey. In the presence of local particle density restrictions (finite prey carrying capacity), there exists an extinction threshold for the predator population. The accompanying continuous non-equilibrium phase transition is governed by the directed-percolation universality class. We employ field-theoretic methods based on the Doi-Peliti representation of the master equation for stochastic particle interaction models to (i) map the ensuing action in the vicinity of the absorbing state phase transition to Reggeon field theory, and (ii) to quantitatively address fluctuation-induced renormalizations of the population oscillation frequency, damping, and diffusion coefficients in the species coexistence phase.
We simulate an individual-based model that represents both the phenotype and genome of digital organisms with predator-prey interactions. We show how open-ended growth of complexity arises from the invariance of genetic evolution operators with respect to changes in the complexity, and that the dynamics which emerges is controlled by a non-equilibrium critical point. The mechanism is analogous to the development of the cascade in fluid turbulence.
Trait-mediated indirect effects are increasingly acknowledged as important components in the dynamics of ecological systems. The hamiltonian form of the LV equations is traditionally modified by adding density dependence to the prey variable and functional response to the predator variable. Enriching these non-linear elements with a trait-mediation added to the carrying capacity of the prey creates the dynamics of critical transitions and hysteretic zones.
Cooperative interactions pervade the dynamics of a broad rage of many-body systems, such as ecological communities, the organization of social structures, and economic webs. In this work, we investigate the dynamics of a simple population model that is driven by cooperative and symmetric interactions between two species. We develop a mean-field and a stochastic description for this cooperative two-species reaction scheme. For an isolated population, we determine the probability to reach a state of fixation, where only one species survives, as a function of the initial concentrations of the two species. We also determine the time to reach the fixation state. When each species can migrate into the population and replace a randomly selected individual, the population reaches a steady state. We show that this steady-state distribution undergoes a unimodal to trimodal transition as the migration rate is decreased beyond a critical value. In this low-migration regime, the steady state is not truly steady, but instead fluctuates strongly between near-fixation states of the two species. The characteristic time scale of these fluctuations diverges as $lambda^{-1}$.