Emergence of diversity in a model ecosystem

Diversity and ecotones in a model ecosystems of sessile species

Sessile species compete for space and accessible light, with directed interactions evident in one species overgrowing another and with multispecies systems characterized by nontransitive relationships. Such patterns are observed in coral reefs or lichens on rock surfaces. Open systems with episodic invasions of such species have been predicted to exhibit a stable high-diversity state when the interaction probability is below a certain critical threshold. Here, we explore this metastable high-diversity state and find that the diversity in the high-diversity state scales with the square root of the system area. When introducing two different environments, we predict a hugely increased diversity along mutual environment border. Further, the presence of spatially segregated environments is predicted to allow for increased robustness of the high-diversity state.

Emerging diversity in a population of evolving intransitive dice

Exploiting the mathematical curiosity of intransitive dice, we present a simple theoretical model for coevolution that captures scales ranging from the genome of the individual to the system-wide emergence of species diversity. We study a set of evolving agents that interact competitively in a closed system, in which both the dynamics of mutations and competitive advantage emerge directly from interpreting a genome as the sides of a die. The model demonstrates sympatric speciation where new species evolve from existing ones while in contact with the entire ecosystem. Allowing free mutations both in the genomes and the mutation rates, we find, in contrast to hierarchical models of fitness, the emergence of a metastable state of finite mutation rate and diversity.

Why is cyclic dominance so rare?

Natural populations can contain multiple types of coexisting individuals. How does natural selection maintain such diversity within and across populations? A popular theoretical basis for the maintenance of diversity is cyclic dominance, illustrated by the rock-paper-scissor game. However, it appears difficult to find cyclic dominance in nature. Why is this the case? Focusing on continuously produced novel mutations, we theoretically addressed the rareness of cyclic dominance. We developed a model of an evolving population and studied the formation of cyclic dominance. Our results showed that the chance for cyclic dominance to emerge is lower when the newly introduced type is similar to existing types compared to the introduction of an unrelated type. This suggests that cyclic dominance is more likely to evolve through the assembly of unrelated types whereas it rarely evolves within a community of similar types.

Spatial ecology of territorial populations

Many ecosystems, from vegetation to biofilms, are composed of territorial populations that compete for both nutrients and physical space. What are the implications of such spatial organization for biodiversity? To address this question, we developed and analyzed a model of territorial resource competition. In the model, all species obey trade-offs inspired by biophysical constraints on metabolism; the species occupy nonoverlapping territories, while nutrients diffuse in space. We find that the nutrient diffusion time is an important control parameter for both biodiversity and the timescale of population dynamics. Interestingly, fast nutrient diffusion allows the populations of some species to fluctuate to zero, leading to extinctions. Moreover, territorial competition spontaneously gives rise to both multistability and the Allee effect (in which a minimum population is required for survival), so that small perturbations can have major ecological effects. While the assumption of trade-offs allows for the coexistence of more species than the number of nutrients-thus violating the principle of competitive exclusion-overall biodiversity is curbed by the domination of "oligotroph" species. Importantly, in contrast to well-mixed models, spatial structure renders diversity robust to inequalities in metabolic trade-offs. Our results suggest that territorial ecosystems can display high biodiversity and rich dynamics simply due to competition for resources in a spatial community.

Trade-off shapes diversity in eco-evolutionary dynamics

We introduce an Interaction- and Trade-off-based Eco-Evolutionary Model (ITEEM), in which species are competing in a well-mixed system, and their evolution in interaction trait space is subject to a life-history trade-off between replication rate and competitive ability. We demonstrate that the shape of the trade-off has a fundamental impact on eco-evolutionary dynamics, as it imposes four phases of diversity, including a sharp phase transition. Despite its minimalism, ITEEM produces a remarkable range of patterns of eco-evolutionary dynamics that are observed in experimental and natural systems. Most notably we find self-organization towards structured communities with high and sustained diversity, in which competing species form interaction cycles similar to rock-paper-scissors games.

Models of life: epigenetics, diversity and cycles

This review emphasizes aspects of biology that can be understood through repeated applications of simple causal rules. The selected topics include perspectives on gene regulation, phage lambda development, epigenetics, microbial ecology, as well as model approaches to diversity and to punctuated equilibrium in evolution. Two outstanding features are repeatedly described. One is the minimal number of rules to sustain specific states of complex systems for a long time. The other is the collapse of such states and the subsequent dynamical cycle of situations that restitute the system to a potentially new metastable state.

Characterization of phase transitions in a model ecosystem of sessile species

We consider a model ecosystem of sessile species competing for space. In particular, we consider the system introduced by Mathiesen et al. [J. Mathiesen, N. Mitarai, K. Sneppen, and A. Trusina, Phys. Rev. Lett. 107, 188101 (2011)PRLTAO0031-900710.1103/PhysRevLett.107.188101] where species compete according to a fixed interaction network with links determined by a Bernoulli process. In the limit of a small introduction rate of new species, the model exhibits a discontinuous transition from a high-diversity state to a low-diversity state as the interaction probability between species, γ, is increased from zero. Here we explore the effects of finite introduction rates and system size on the phase transition by utilizing efficient parallel computing. We find that the low state appears for γ>γ_{c}. As γ is increased further, the high state approaches the low state, suggesting the possibility that the two states merge at a high γ. We find that the fraction of time spent in the high state becomes longer with higher introduction rates, but the availability of the two states is rather insensitive to the value of the introduction rate. Furthermore, we establish a relation between the introduction rate and the system size, which preserves the probability for the system to remain in the high-diversity state.

Three is much more than two in coarsening dynamics of cyclic competitions

The classical game of rock-paper-scissors has inspired experiments and spatial model systems that address the robustness of biological diversity. In particular, the game nicely illustrates that cyclic interactions allow multiple strategies to coexist for long-time intervals. When formulated in terms of a one-dimensional cellular automata, the spatial distribution of strategies exhibits coarsening with algebraically growing domain size over time, while the two-dimensional version allows domains to break and thereby opens the possibility for long-time coexistence. We consider a quasi-one-dimensional implementation of the cyclic competition, and study the long-term dynamics as a function of rare invasions between parallel linear ecosystems. We find that increasing the complexity from two to three parallel subsystems allows a transition from complete coarsening to an active steady state where the domain size stays finite. We further find that this transition happens irrespective of whether the update is done in parallel for all sites simultaneously or done randomly in sequential order. In both cases, the active state is characterized by localized bursts of dislocations, followed by longer periods of coarsening. In the case of the parallel dynamics, we find that there is another phase transition between the active steady state and the coarsening state within the three-line system when the invasion rate between the subsystems is varied. We identify the critical parameter for this transition and show that the density of active boundaries has critical exponents that are consistent with the directed percolation universality class. On the other hand, numerical simulations with the random sequential dynamics suggest that the system may exhibit an active steady state as long as the invasion rate is finite.

Speciation, diversification, and coexistence of sessile species that compete for space

Speciation, diversification, and competition between species challenge the stability of complex ecosystems. Laboratory experiments often focus on one or two species competing under conditions where they may grow exponentially. Field studies, in contrast, emphasize multi-species communities characterized by many types of ecological interactions. A general problem is to understand conditions that support a dynamically maintained coexistence of many species in an ecosystem over a long time span. In the present paper we propose a lattice model of multiple competing and evolving sessile species. When allowing the interspecies interactions to mutate, we obtain coexistence of many species in a complex ecosystem, provided that there is a cost for each interaction. The diversity reached by the model incorporating speciation is found to be substantially higher than in the case when entirely new species appear due to immigration from outside of the considered ecosystem. The species self-organize their spatial distribution through competitive interactions to create many patches, implicitly protecting each other from competitively superior species, and speciation in each patch leads the system to high diversity. We also show that species that exist a long time tend to have a relatively small population, as this allows them to avoid encounter with competitive invaders.

Geometry for a penguin-albatross rookery

We introduce a simple ecological model describing the spatial organization of two interacting populations whose individuals are indifferent to conspecifics and avoid the proximity to heterospecifics. At small population densities Φ a nontrivial structure is observed where clusters of individuals arrange into a rhomboidal bipartite network with an average degree of 4. For Φ → 0 the length scale, order parameter, and susceptibility of the network exhibit power-law divergences compatible with hyperscaling, suggesting the existence of a zero-density nontrivial critical point. At larger densities a critical threshold Φ(c) is identified above which the evolution toward a partially ordered configuration is prevented and the system becomes jammed in a fully mixed state.

Disturbance accelerates the transition from low- to high-diversity state in a model ecosystem

The effect of disturbance on a model ecosystem of sessile and mutually competitive species [Mathiesen et al., Phys. Rev. Lett. 107, 188101 (2011); Mitarai et al., Phys. Rev. E 86, 011929 (2012)] is studied. The disturbance stochastically removes individuals from the system, and the created empty sites are recolonized by neighboring species. We show that the stable high-diversity state, maintained by occasional cyclic species interactions that create isolated patches of metapopulations, is robust against small disturbance. We further demonstrate that finite disturbance can accelerate the transition from the low- to high-diversity state by helping the creation of small patches through diffusion of boundaries between species with standoff relations.

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.