Eur. Phys. J. B (2012) 85: 141
https://doi.org/10.1140/epjb/e2012-20918-4

Regular Article

On the relationship between cyclic and hierarchical three-species predator-prey systems and the two-species Lotka-Volterra model

Department of Physics, Virginia Tech, Blacksburg, Virginia 24061-0435, USA

^a e-mail: tauber@vt.edu

Received: 10 November 2011
Published online: 26 April 2012

Abstract

Stochastic spatial predator-prey competition models represent paradigmatic systems to understand the emergence of biodiversity and the stability of ecosystems. We aim to clarify the relationship and connections between interacting three-species models and the classic two-species Lotka-Volterra (LV) model that entails predator-prey coexistence with long-lived population oscillations. To this end, we utilize mean-field theory and Monte Carlo simulations on two-dimensional square lattices to explore the temporal evolution characteristics of two different interacting three-species predator-prey systems, namely: (1) a cyclic rock-paper-scissors (RPS) model with conserved total particle number but strongly asymmetric reaction rates that lets the system evolve towards one “corner” of configuration space; (2) a hierarchical “food chain” where an additional intermediate species is inserted between the predator and prey in the LV model. For the asymmetric cyclic model variant (1), we demonstrate that the evolutionary properties of both minority species in the (quasi-)steady state of this stochastic spatial three-species “corner” RPS model are well approximated by the two-species LV system, with its emerging characteristic features of localized population clustering, persistent oscillatory dynamics, correlated spatio-temporal patterns, and fitness enhancement through quenched spatial disorder in the predation rates. In contrast, we could not identify any regime where the hierarchical three-species model (2) would reduce to the two-species LV system. In the presence of pair exchange processes, the system remains essentially well-mixed, and we generally find the Monte Carlo simulation results for the spatially extended hierarchical model (2) to be consistent with the predictions from the corresponding mean-field rate equations. If spreading occurs only through nearest-neighbor hopping, small population clusters emerge; yet the requirement of an intermediate species cluster obviously disrupts spatio-temporal correlations between predator and prey, and correspondingly eliminates many of the intriguing fluctuation phenomena that characterize the stochastic spatial LV system.