Stationary pattern of vortices or strings in biological systems: Lattice version of the Lotka-Volterra model

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 spatiotemporal 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 semiquantitative description of the (quasi-)stationary state.

Population persistence under two conservation measures: Paradox of habitat protection in a patchy environment

Anthropogenic modification of natural habitats is a growing threat to biodiversity and ecosystem services. The protection of biospecies has become increasingly important. Here, we pay attention to a single species as a conservation target. The species has three processes: reproduction, death and movement. Two different measures of habitat protection are introduced. One is partial protection in a single habitat (patch); the mortality rate of the species is reduced inside a rectangular area. The other is patch protection in a two-patch system, where only the mortality rate in a particular patch is reduced. For the one-patch system, we carry out computer simulations of a stochastic cellular automaton for a "contact process". Individual movements follow random walking. For the two-patch system, we assume an individual migrates into the empty cell in the destination patch. The reaction-diffusion equation (RDE) is derived, whereby the recently developed "swapping migration" is used. It is found that both measures are mostly effective for population persistence. However, comparing the results of the two measures revealed different behaviors. ⅰ) In the case of the one-patch system, the steady-state densities in protected areas are always higher than those in wild areas. However, in the two-patch system, we have found a paradox: the densities in protected areas can be lower than those in wild areas. ⅱ) In the two-patch system, we have found another paradox: the total density in both patches can be lower, even though the proportion of the protected area is larger. Both paradoxes clearly occur for the RDE with swapping migration.

Metapopulation model for a prey-predator system: Nonlinear migration due to the finite capacities of patches

Many species live in spatially separated patches, and individuals can migrate between patches through paths. In real ecosystems, the capacities of patches are finite. If a patch is already occupied by the individuals of some species, then the migration into the patch is impossible. In the present paper, we deal with prey-predator system composed of two patches. Each patch contains a limited number of cells, where the cell is either empty or occupied by an individual of prey or predator. We introduce "swapping migration" defined by the exchange between occupied and empty cells. An individual can migrate, only when there are empty cells in the destination patch. Reaction-migration equations in prey-predator system are presented, where the migration term forms nonlinear function of densities. We numerically solve equilibrium densities, and find that the population dynamics are largely affected by nonlinear migration. Not only extinction points but also the responses to the environmental changes crucially depend on the patch capacities.

Metapopulation model of rock-scissors-paper game with subpopulation-specific victory rates stabilized by heterogeneity

Recently, metapopulation models for rock-paper-scissors games have been presented. Each subpopulation is represented by a node on a graph. An individual is either rock (R), scissors (S) or paper (P); it randomly migrates among subpopulations. In the present paper, we assume victory rates differ in different subpopulations. To investigate the dynamic state of each subpopulation (node), we numerically obtain the solutions of reaction-diffusion equations on the graphs with two and three nodes. In the case of homogeneous victory rates, we find each subpopulation has a periodic solution with neutral stability. However, when victory rates between subpopulations are heterogeneous, the solution approaches stable focuses. The heterogeneity of victory rates promotes the coexistence of species.

Dynamics of Intraguild Predation Systems with Intraspecific Competition

This paper considers intraguild predation (IGP) systems where species in the same community kill and eat each other and there is intraspecific competition in each species. The IGP systems are characterized by a lattice gas model, in which reaction between sites on the lattice occurs in a random and independent way. Global dynamics of the model with two species demonstrate mechanisms by which IGP leads to survival/extinction of species. It is shown that an intermediary level of predation promotes survival of species, while over-predation or under-predation could result in species extinction. An interesting result is that increasing intraspecific competition of one species can lead to extinction of one or both species, while increasing intraspecific competitions of both species would result in coexistence of species in facultative predation. Initial population densities of the species are also shown to play a role in persistence of the system. Then the analysis is extended to IGP systems with one species. Numerical simulations confirm and extend our results.

Metapopulation model for rock-paper-scissors game: Mutation affects paradoxical impacts

The rock-paper-scissors (RPS) game is known as one of the simplest cyclic dominance models. This game is key to understanding biodiversity. Three species, rock (R), paper (P) and scissors (S), can coexist in nature. In the present paper, we first present a metapopulation model for RPS game with mutation. Only mutation from R to S is allowed. The total population consists of spatially separated patches, and the mutation occurs in particular patches. We present reaction-diffusion equations which have two terms: reaction and migration terms. The former represents the RPS game with mutation, while the latter corresponds to random walk. The basic equations are solved analytically and numerically. It is found that the mutation induces one of three phases: the stable coexistence of three species, the stable phase of two species, and a single-species phase. The phase transitions among three phases occur by varying the mutation rate. We find the conditions for coexistence are largely changed depending on metapopulation models. We also find that the mutation induces different paradoxes in different patches.

Transition of interaction outcomes in a facilitation-competition system of two species

A facilitation-competition system of two species is considered, where one species has a facilitation effect on the other but there is spatial competition between them. Our aim is to show mechanism by which the facilitation promotes coexistence of the species. A lattice gas model describing the facilitation-competition system is analyzed, in which nonexistence of periodic solution is shown and previous results are extended. Global dynamics of the model demonstrate essential features of the facilitation-competition system. When a species cannot survive alone, a strong facilitation from the other would lead to its survival. However, if the facilitation is extremely strong, both species go extinct. When a species can survive alone and its mortality rate is not larger than that of the other species, it would drive the other one into extinction. When a species can survive alone and its mortality rate is larger than that of the other species, it would be driven into extinction if the facilitation from the other is weak, while it would coexist with the other if the facilitation is strong. Moreover, an extremely strong facilitation from the other would lead to extinction of species. Bifurcation diagram of the system exhibits that interaction outcome between the species can transition between competition, amensalism, neutralism and parasitism in a smooth fashion. A novel result of this paper is the rigorous and thorough analysis, which displays transparency of dynamics in the system. Numerical simulations validate the results.

Heterogeneous network promotes species coexistence: metapopulation model for rock-paper-scissors game

Understanding mechanisms of biodiversity has been a central question in ecology. The coexistence of three species in rock-paper-scissors (RPS) systems are discussed by many authors; however, the relation between coexistence and network structure is rarely discussed. Here we present a metapopulation model for RPS game. The total population is assumed to consist of three subpopulations (nodes). Each individual migrates by random walk; the destination of migration is randomly determined. From reaction-migration equations, we obtain the population dynamics. It is found that the dynamic highly depends on network structures. When a network is homogeneous, the dynamics are neutrally stable: each node has a periodic solution, and the oscillations synchronize in all nodes. However, when a network is heterogeneous, the dynamics approach stable focus and all nodes reach equilibriums with different densities. Hence, the heterogeneity of the network promotes biodiversity.

Survival behavior in the cyclic Lotka-Volterra model with a randomly switching reaction rate

We study the influence of a randomly switching reproduction-predation rate on the survival behavior of the nonspatial cyclic Lotka-Volterra model, also known as the zero-sum rock-paper-scissors game, used to metaphorically describe the cyclic competition between three species. In large and finite populations, demographic fluctuations (internal noise) drive two species to extinction in a finite time, while the species with the smallest reproduction-predation rate is the most likely to be the surviving one (law of the weakest). Here we model environmental (external) noise by assuming that the reproduction-predation rate of the strongest species (the fastest to reproduce and predate) in a given static environment randomly switches between two values corresponding to more and less favorable external conditions. We study the joint effect of environmental and demographic noise on the species survival probabilities and on the mean extinction time. In particular, we investigate whether the survival probabilities follow the law of the weakest and analyze their dependence on the external noise intensity and switching rate. Remarkably, when, on average, there is a finite number of switches prior to extinction, the survival probability of the predator of the species whose reaction rate switches typically varies nonmonotonically with the external noise intensity (with optimal survival about a critical noise strength). We also outline the relationship with the case where all reaction rates switch on markedly different time scales.

Effect of directional migration on Lotka-Volterra system with desert

Migration is observed across many species. Several authors have studied ecological migration by applying cellular automaton (CA). In this paper, we present a directional migration model with desert on a one-dimensional lattice where a traffic CA model and a lattice Lotka-Volterra system are connected. Here predators correspond to locomotive animals while prey is immobile plants. Predators migrate between deserts and fertile lands repeatedly. Computer simulations reveal the two types of phase transition: coexistence of both species and prey dominance, which is caused by both benefit and cost of migration. In the coexistence phase, the steady-state density of predators usually increases by migration as long as the desert size is small and their mortality rate is low. In contrast, the prey density increases, even if the desert size becomes large. Such a paradox comes from the indirect effect: predators go extinct by the increase of desert size, so that the plant density can increase. Moreover, we find several self-organized spatial patterns: 1) predators form a stripe pattern; namely swarms. 2) The velocity of predators is high on deserts, but very low on fertile land. 3) Predators give birth only on fertile lands.

Microhabitat locality allows multi-species coexistence in terrestrial plant communities

Most terrestrial plant communities exhibit relatively high species diversity and many competitive species are ubiquitous. Many theoretical studies have been carried out to investigate the coexistence of a few competitive species and in most cases they suggest competitive exclusion. Theoretical studies have revealed that coexistence of even three or four species can be extremely difficult. It has been suggested that the coexistence of many species has been achieved by the fine differences in suitable microhabitats for each species, attributing to niche-separation. So far there is no explicit demonstration of such a coexistence in mathematical and simulation studies. Here we built a simple lattice Lotka-Volterra model of competition by incorporating the minute differences of suitable microhabitats for many species. By applying the site variations in species-specific settlement rates of a seedling, we achieved the coexistence of more than 10 species. This result indicates that competition between many species is avoided by the spatial variations in species-specific microhabitats. Our results demonstrate that coexistence of many species becomes possible by the minute differences in microhabitats. This mechanism should be applicable to many vegetation types, such as temperate forests and grasslands.

How turbulence regulates biodiversity in systems with cyclic competition

Cyclic, nonhierarchical interactions among biological species represent a general mechanism by which ecosystems are able to maintain high levels of biodiversity. However, species coexistence is often possible only in spatially extended systems with a limited range of dispersal, whereas in well-mixed environments models for cyclic competition often lead to a loss of biodiversity. Here we consider the dispersal of biological species in a fluid environment, where mixing is achieved by a combination of advection and diffusion. In particular, we perform a detailed numerical analysis of a model composed of turbulent advection, diffusive transport, and cyclic interactions among biological species in two spatial dimensions and discuss the circumstances under which biodiversity is maintained when external environmental conditions, such as resource supply, are uniform in space. Cyclic interactions are represented by a model with three competitors, resembling the children's game of rock-paper-scissors, whereas the flow field is obtained from a direct numerical simulation of two-dimensional turbulence with hyperviscosity. It is shown that the space-averaged dynamics undergoes bifurcations as the relative strengths of advection and diffusion compared to biological interactions are varied.

Population dynamics of intraguild predation in a lattice gas system

In the system of intraguild predation (IGP) we are concerned with, species that are in a predator-prey relationship, also compete for shared resources (space or food). While several models have been established to characterize IGP, mechanisms by which IG prey and IG predator can coexist in IGP systems with spatial competition, have not been shown. This paper considers an IGP model, which is derived from reactions on lattice and has a form similar to that of Lotka-Volterra equations. Dynamics of the model demonstrate properties of IGP and mechanisms by which the IGP leads to coexistence of species and occurrence of alternative states. Intermediate predation is shown to lead to persistence of the predator, while extremely big predation can lead to extinction of one/both species and extremely small predation can lead to extinction of the predator. Numerical computations confirm and extend our results. While empirical observations typically exhibit coexistence of IG predator and IG prey, theoretical analysis in this work demonstrates exact conditions under which this coexistence can occur.

Cyclic dominance in evolutionary games: a review

Rock is wrapped by paper, paper is cut by scissors and scissors are crushed by rock. This simple game is popular among children and adults to decide on trivial disputes that have no obvious winner, but cyclic dominance is also at the heart of predator-prey interactions, the mating strategy of side-blotched lizards, the overgrowth of marine sessile organisms and competition in microbial populations. Cyclical interactions also emerge spontaneously in evolutionary games entailing volunteering, reward, punishment, and in fact are common when the competing strategies are three or more, regardless of the particularities of the game. Here, we review recent advances on the rock-paper-scissors (RPS) and related evolutionary games, focusing, in particular, on pattern formation, the impact of mobility and the spontaneous emergence of cyclic dominance. We also review mean-field and zero-dimensional RPS models and the application of the complex Ginzburg-Landau equation, and we highlight the importance and usefulness of statistical physics for the successful study of large-scale ecological systems. Directions for future research, related, for example, to dynamical effects of coevolutionary rules and invasion reversals owing to multi-point interactions, are also outlined.

Self-organizing patterns in an evolutionary rock-paper-scissors game for stochastic synchronized strategy updates

We study a spatial evolutionary rock-paper-scissors game with synchronized strategy updating. Players gain their payoff from games with their four neighbors on a square lattice and can update their strategies simultaneously according to the logit rule, which is the noisy version of the best-response dynamics. For the synchronized strategy update two types of global oscillations (with an ordered strategy arrangement and periods of three and six generations) can occur in this system in the zero noise limit. At low noise values, all nine oscillating phases are present in the system by forming a self-organizing spatial pattern due to the comprising invasion and speciation processes along the interfaces separating the different domains.

Characterization of spiraling patterns in spatial rock-paper-scissors games

The spatiotemporal arrangement of interacting populations often influences the maintenance of species diversity and is a subject of intense research. Here, we study the spatiotemporal patterns arising from the cyclic competition between three species in two dimensions. Inspired by recent experiments, we consider a generic metapopulation model comprising "rock-paper-scissors" interactions via dominance removal and replacement, reproduction, mutations, pair exchange, and hopping of individuals. By combining analytical and numerical methods, we obtain the model's phase diagram near its Hopf bifurcation and quantitatively characterize the properties of the spiraling patterns arising in each phase. The phases characterizing the cyclic competition away from the Hopf bifurcation (at low mutation rate) are also investigated. Our analytical approach relies on the careful analysis of the properties of the complex Ginzburg-Landau equation derived through a controlled (perturbative) multiscale expansion around the model's Hopf bifurcation. Our results allow us to clarify when spatial "rock-paper-scissors" competition leads to stable spiral waves and under which circumstances they are influenced by nonlinear mobility.

The paradox of enrichment in phytoplankton by induced competitive interactions

The biodiversity loss of phytoplankton with eutrophication has been reported in many aquatic ecosystems, e.g., water pollution and red tides. This phenomenon seems similar, but different from the paradox of enrichment via trophic interactions, e.g., predator-prey systems. We here propose the paradox of enrichment by induced competitive interactions using multiple contact process (a lattice Lotka-Volterra competition model). Simulation results demonstrate how eutrophication invokes more competitions in a competitive ecosystem resulting in the loss of phytoplankton diversity in ecological time. The paradox is enhanced under local interactions, indicating that the limited dispersal of phytoplankton reduces interspecific competition greatly. Thus, the paradox of enrichment appears when eutrophication destroys an ecosystem either by elevated interspecific competition within a trophic level and/or destabilization by trophic interactions. Unless eutrophication due to human activities is ceased, the world's aquatic ecosystems will be at risk.

Diverging fluctuations in a spatial five-species cyclic dominance game

A five-species predator-prey model is studied on a square lattice where each species has two prey and two predators on the analogy to the rock-paper-scissors-lizard-Spock game. The evolution of the spatial distribution of species is governed by site exchange and invasion between the neighboring predator-prey pairs, where the cyclic symmetry can be characterized by two different invasion rates. The mean-field analysis has indicated periodic oscillations in the species densities with a frequency becoming zero for a specific ratio of invasion rates. When varying the ratio of invasion rates, the appearance of this zero-eigenvalue mode is accompanied by neutrality between the species associations. Monte Carlo simulations of the spatial system reveal diverging fluctuations at a specific invasion rate, which can be related to the vanishing dominance between all pairs of species associations.

Global attractors and extinction dynamics of cyclically competing species

Transitions to absorbing states are of fundamental importance in nonequilibrium physics as well as ecology. In ecology, absorbing states correspond to the extinction of species. We here study the spatial population dynamics of three cyclically interacting species. The interaction scheme comprises both direct competition between species as in the cyclic Lotka-Volterra model, and separated selection and reproduction processes as in the May-Leonard model. We show that the dynamic processes leading to the transient maintenance of biodiversity are closely linked to attractors of the nonlinear dynamics for the overall species' concentrations. The characteristics of these global attractors change qualitatively at certain threshold values of the mobility and depend on the relative strength of the different types of competition between species. They give information about the scaling of extinction times with the system size and thereby the stability of biodiversity. We define an effective free energy as the negative logarithm of the probability to find the system in a specific global state before reaching one of the absorbing states. The global attractors then correspond to minima of this effective energy landscape and determine the most probable values for the species' global concentrations. As in equilibrium thermodynamics, qualitative changes in the effective free energy landscape indicate and characterize the underlying nonequilibrium phase transitions. We provide the complete phase diagrams for the population dynamics and give a comprehensive analysis of the spatio-temporal dynamics and routes to extinction in the respective phases.

Coexistence and survival in conservative Lotka-Volterra networks

Analyzing coexistence and survival scenarios of Lotka-Volterra (LV) networks in which the total biomass is conserved is of vital importance for the characterization of long-term dynamics of ecological communities. Here, we introduce a classification scheme for coexistence scenarios in these conservative LV models and quantify the extinction process by employing the Pfaffian of the network's interaction matrix. We illustrate our findings on global stability properties for general systems of four and five species and find a generalized scaling law for the extinction time.

Effects of competition on pattern formation in the rock-paper-scissors game

We investigate the impact of cyclic competition on pattern formation in the rock-paper-scissors game. By separately considering random and prepared initial conditions, we observe a critical influence of the competition rate p on the stability of spiral waves and on the emergence of biodiversity. In particular, while increasing values of p promote biodiversity, they may act detrimentally on spatial pattern formation. For random initial conditions, we observe a phase transition from biodiversity to an absorbing phase, whereby the critical value of mobility grows linearly with increasing values of p on a log-log scale but then saturates as p becomes large. For prepared initial conditions, we observe the formation of single-armed spirals, but only for values of p that are below a critical value. Once above that value, the spirals break up and form disordered spatial structures, mainly because of the percolation of vacant sites. Thus there exists a critical value of the competition rates p(c) for stable single-armed spirals in finite populations. Importantly though, p(c) increases with increasing system size because noise reinforces the disintegration of ordered patterns. In addition, we also find that p(c) increases with the mobility. These phenomena are reproduced by a deterministic model that is based on nonlinear partial differential equations. Our findings indicate that competition is vital for the sustenance of biodiversity and the emergence of pattern formation in ecosystems governed by cyclical interactions.

Three- and four-state rock-paper-scissors games with diffusion

Cyclic dominance of three species is a commonly occurring interaction dynamics, often denoted the rock-paper-scissors (RPS) game. Such a type of interactions is known to promote species coexistence. Here, we generalize recent results of Reichenbach [Nature (London) 448, 1046 (2007)] of a four-state variant of the RPS game. We show that spiral formation takes place only without a conservation law for the total density. Nevertheless, in general, fast diffusion can destroy species coexistence. We also generalize the four-state model to slightly varying reaction rates. This is shown both analytically and numerically not to change pattern formation, or the effective wavelength of the spirals, and therefore not to alter the qualitative properties of the crossover to extinction.

Time correlation function in systems with two coexisting biological species

We study a stochastic lattice model describing the dynamics of coexistence of two interacting biological species. The model comprehends the local processes of birth, death, and diffusion of individuals of each species and is grounded on interaction of the predator-prey type. The species coexistence can be of two types: With self-sustained coupled time oscillations of population densities and without oscillations. We perform numerical simulations of the model on a square lattice and analyze the temporal behavior of each species by computing the time correlation functions as well as the spectral densities. This analysis provides an appropriate characterization of the different types of coexistence. It is also used to examine linked population cycles in nature and in experiment.

Self-organizing patterns maintained by competing associations in a six-species predator-prey model

Formation and competition of associations are studied in a six-species ecological model where each species has two predators and two prey. Each site of a square lattice is occupied by an individual belonging to one of the six species. The evolution of the spatial distribution of species is governed by iterated invasions between the neighboring predator-prey pairs with species specific rates and by site exchange between the neutral pairs with a probability X . This dynamical rule yields the formation of five associations composed of two or three species with proper spatiotemporal patterns. For large X a cyclic dominance can occur between the three two-species associations whereas one of the two three-species associations prevails in the whole system for low values of X in the final state. Within an intermediate range of X all the five associations coexist due to the fact that cyclic invasions between the two-species associations reduce their resistance temporarily against the invasion of three-species associations.

Noise and correlations in a spatial population model with cyclic competition

Noise and spatial degrees of freedom characterize most ecosystems. Some aspects of their influence on the coevolution of populations with cyclic interspecies competition have been demonstrated in recent experiments [e.g., B. Kerr, Nature (London) 418, 171 (2002)10.1038/nature00823]. To reach a better theoretical understanding of these phenomena, we consider a paradigmatic spatial model where three species exhibit cyclic dominance. Using an individual-based description, as well as stochastic partial differential and deterministic reaction-diffusion equations, we account for stochastic fluctuations and spatial diffusion at different levels and show how fascinating patterns of entangled spirals emerge. We rationalize our analysis by computing the spatiotemporal correlation functions and provide analytical expressions for the front velocity and the wavelength of the propagating spiral waves.

Segregation process and phase transition in cyclic predator-prey models with an even number of species

We study a spatial cyclic predator-prey model with an even number of species (for n=4, 6, and 8) that allows the formation of two defensive alliances consisting of the even and odd label species. The species are distributed on the sites of a square lattice. The evolution of spatial distribution is governed by iteration of two elementary processes on neighboring sites chosen randomly: if the sites are occupied by a predator-prey pair then the predator invades the prey's site; otherwise the species exchange their sites with a probability X . For low X values, a self-organizing pattern is maintained by cyclic invasions. If X exceeds a threshold value, then two types of domain grow up that are formed by the odd and even label species, respectively. Monte Carlo simulations indicate the blocking of this segregation process within a range of X for n=8.

Cyclical interactions with alliance-specific heterogeneous invasion rates

We study a six-species Lotka-Volterra-type system on different two-dimensional lattices when each species has two superior and two inferior partners. The invasion rates from predator sites to a randomly chosen neighboring prey site depend on the predator-prey pair, whereby cyclic symmetries within the two three-species defensive alliances are conserved. Monte Carlo simulations reveal an unexpected nonmonotonous dependence of alliance survival on the difference of alliance-specific invasion rates. This behavior is qualitatively reproduced by a four-point mean-field approximation. The study addresses fundamental problems of stability for the competition of two defensive alliances and thus has important implications in natural and social sciences.

Evolutionary ecology in silico: Does mathematical modelling help in understanding 'generic' trends?

Motivated by the results of recent laboratory experiments, as well as many earlier field observations, that evolutionary changes can take place in ecosystems over relatively short ecological time scales, several 'unified' mathematical models of evolutionary ecology have been developed over the last few years with the aim of describing the statistical properties of data related to the evolution of ecosystems. Moreover, because of the availability of sufficiently fast computers, it has become possible to carry out detailed computer simulations of these models. For the sake of completeness and to put these recent developments in perspective, we begin with a brief summary of some older models of ecological phenomena and evolutionary processes. However, the main aim of this article is to review critically these 'unified' models, particularly those published in the physics literature, in simple language that makes the new theories accessible to a wider audience.

Spreading of families in cyclic predator-prey models

We study the spreading of families in two-dimensional multispecies predator-prey systems, in which species cyclically dominate each other. In each time step randomly chosen individuals invade one of the nearest sites of the square lattice eliminating their prey. Initially all individuals get a family name which will be carried on by their descendants. Monte Carlo simulations show that the systems with several species (N=3,4,5) are asymptotically approaching the behavior of the voter model, i.e., the survival probability of families, the mean size of families, and the mean-square distance of descendants from their ancestor exhibits the same scaling behavior. The scaling behavior of the survival probability of families has a logarithmic correction. In case of the voter model this correction depends on the number of species, while cyclic predator-prey models behave like the voter model with infinite species. It is found that changing the rates of invasions does not change this asymptotic behavior. As an application a three-species system with a fourth-species intruder is also discussed.

Phase transition and selection in a four-species cyclic predator-prey model

We study a four-species ecological system with cyclic dominance whose individuals are distributed on a square lattice. Randomly chosen individuals migrate to one of the neighboring sites if it is empty or invade this site if occupied by their prey. The cyclic dominance maintains the coexistence of all four species if the concentration of vacant sites is lower than a threshold value. Above the threshold, a symmetry breaking ordering occurs via growing domains containing only two neutral species inside. These two neutral species can protect each other from the external invaders (predators) and extend their common territory. According to our Monte Carlo simulations the observed phase transition seems to be equivalent to those found in spreading models with two equivalent absorbing states although the present model has continuous sets of absorbing states with different portions of the two neutral species. The selection mechanism yielding symmetric phases is related to the domain growth process with wide boundaries where the four species coexist.

Nonextensivity of the cyclic lattice Lotka-Volterra model

We numerically show that the lattice Lotka-Volterra model, when realized on a square lattice support, gives rise to a finite production, per unit time, of the nonextensive entropy S(q)=(1- summation operator (i)p(q)(i))/(q-1) (S(1)=- summation operator (i)p(i) ln p(i)). This finiteness only occurs for q=0.5 for the d=2 growth mode (growing droplet), and for q=0 for the d=1 one (growing stripe). This strong evidence of nonextensivity is consistent with the spontaneous emergence of local domains of identical particles with fractal boundaries and competing interactions. Such direct evidence is, to our knowledge, exhibited for the first time for a many-body system which, at the mean field level, is conservative.

Mussel disturbance dynamics: signatures of oceanographic forcing from local interactions

Local interactions, biotic and abiotic, can have a strong influence on the large-scale properties of ecosystems. However, ecological models often explore the influence of local biotic interactions where physical disturbance is included as a large-scale and imposed source of variability but is not allowed to interact with biotic processes at the local scale. In marine intertidal communities dominated by mussels, wave disturbances create gaps in the mussel bed that recover through a successional sequence. We present a lattice model of mussel disturbance dynamics that allows local interactions between wave disturbance and mussel recolonization, in which each cell of the lattice can be empty, occupied by a mussel bed element, or disturbed (which corresponds to a newly disturbed cell that has unstable edges). As in natural ecosystems, wave disturbance can also spread from disturbed to adjacent occupied cells, and recolonization can also spread from occupied to adjacent empty cells. We first validate the local rules from artificial gap experiments and from natural gap monitoring along the Oregon coast. We analyze the properties of the model system as a function of different oceanographic forcings of productivity and disturbance. We show that the mussel bed can go through phase transitions characterized by a large sensitivity of mussel cover and patterns to oceanographic forcings but also that criticality (scale invariance) is observed over wide ranges of parameters, which suggests self-organization. We also show that spatial patterns in the intertidal can provide a robust signature of local processes and can inform about oceanographic regimes. We do so by comparing the large-scale patterns of the simulation (scaling exponents) with field data, which suggest that some experimental sites are close to criticality. Our results suggest that regional patterns in disturbed populations can be explained by local biotic and abiotic processes submitted to oceanographic forcing.