Species abundance patterns in complex evolutionary dynamics

Replicator equations induced by microscopic processes in nonoverlapping population playing bimatrix games

This paper is concerned with exploring the microscopic basis for the discrete versions of the standard replicator equation and the adjusted replicator equation. To this end, we introduce frequency-dependent selection-as a result of competition fashioned by game-theoretic consideration-into the Wright-Fisher process, a stochastic birth-death process. The process is further considered to be active in a generation-wise nonoverlapping finite population where individuals play a two-strategy bimatrix population game. Subsequently, connections among the corresponding master equation, the Fokker-Planck equation, and the Langevin equation are exploited to arrive at the deterministic discrete replicator maps in the limit of infinite population size.

Statistical physics of community ecology: a cavity solution to MacArthur's consumer resource model

A central question in ecology is to understand the ecological processes that shape community structure. Niche-based theories have emphasized the important role played by competition for maintaining species diversity. Many of these insights have been derived using MacArthur's consumer resource model (MCRM) or its generalizations. Most theoretical work on the MCRM has focused on small ecosystems with a few species and resources. However theoretical insights derived from small ecosystems many not scale up large ecosystems with many resources and species because large systems with many interacting components often display new emergent behaviors that cannot be understood or deduced from analyzing smaller systems. To address these shortcomings, we develop a statistical physics inspired cavity method to analyze MCRM when both the number of species and the number of resources is large. Unlike previous work in this limit, our theory addresses resource dynamics and resource depletion and demonstrates that species generically and consistently perturb their environments and significantly modify available ecological niches. We show how our cavity approach naturally generalizes niche theory to large ecosystems by accounting for the effect of collective phenomena on species invasion and ecological stability. Our theory suggests that such phenomena are a generic feature of large, natural ecosystems and must be taken into account when analyzing and interpreting community structure. It also highlights the important role that statistical-physics inspired approaches can play in furthering our understanding of ecology.

Ecological communities with Lotka-Volterra dynamics

Ecological communities in heterogeneous environments assemble through the combined effect of species interaction and migration. Understanding the effect of these processes on the community properties is central to ecology. Here we study these processes for a single community subject to migration from a pool of species, with population dynamics described by the generalized Lotka-Volterra equations. We derive exact results for the phase diagram describing the dynamical behaviors, and for the diversity and species abundance distributions. A phase transition is found from a phase where a unique globally attractive fixed point exists to a phase where multiple dynamical attractors exist, leading to history-dependent community properties. The model is shown to possess a symmetry that also establishes a connection with other well-known models.

Multiple peaks of species abundance distributions induced by sparse interactions

We investigate the replicator dynamics with "sparse" symmetric interactions which represent specialist-specialist interactions in ecological communities. By considering a large self-interaction u, we conduct a perturbative expansion which manifests that the nature of the interactions has a direct impact on the species abundance distribution. The central results are all species coexistence in a realistic range of the model parameters and that a certain discrete nature of the interactions induces multiple peaks in the species abundance distribution, providing the possibility of theoretically explaining multiple peaks observed in various field studies. To get more quantitative information, we also construct a non-perturbative theory which becomes exact on tree-like networks if all the species coexist, providing exact critical values of u below which extinct species emerge. Numerical simulations in various different situations are conducted and they clarify the robustness of the presented mechanism of all species coexistence and multiple peaks in the species abundance distributions.

The Effect of Intra- and Interspecific Competition on Coexistence in Multispecies Communities

For two competing species, intraspecific competition must exceed interspecific competition for coexistence. To generalize this well-known criterion to multiple competing species, one must take into account both the distribution of interaction strengths and community structure. Here we derive a multispecies generalization of the two-species rule in the context of symmetric Lotka-Volterra competition and obtain explicit stability conditions for random competitive communities. We then explore the influence of community structure on coexistence. Results show that both the most and least stabilized cases have striking global structures, with a nested pattern emerging in both cases. The distribution of intraspecific coefficients leading to the most and least stabilized communities also follows a predictable pattern that can be justified analytically. In addition, we show that the size of the parameter space allowing for feasible communities always increases with the strength of intraspecific effects in a characteristic way that is independent of the interspecific interaction structure. We conclude by discussing possible extensions of our results to nonsymmetric competition.

Effects of spatial frequency distributions on amplitude death in an array of coupled Landau-Stuart oscillators

The influences of spatial frequency distributions on complete amplitude death are explored by studying an array of diffusively coupled oscillators. We found that with all possible sets of spatial frequency distributions, the two critical coupling strengths ε(c1) (lower-bounded value) and ε(c2) (upper-bounded value) needed to get complete amplitude death exhibit a universal power law and a log-normal distribution respectively, which has long tails in both cases. This is significant for dynamics control, since large variations of ε(c1) and ε(c2) are possible for some spatial arrangements. Moreover, we explore optimal spatial distributions with the smallest (largest) ε(c1) or ε(c2).

Multivariate Moran process with Lotka-Volterra phenomenology

For a population with any given number of types, we construct a new multivariate Moran process with frequency-dependent selection and establish, analytically, a correspondence to equilibrium Lotka-Volterra phenomenology. This correspondence, on the one hand, allows us to infer the phenomenology of our Moran process based on much simpler Lokta-Volterra phenomenology and, on the other, allows us to study Lotka-Volterra dynamics within the finite populations of a Moran process. Applications to community ecology, population genetics, and evolutionary game theory are discussed.

The statistical mechanics of community assembly and species distribution

• Theoretically, communities at or near their equilibrium species number resist entry of new species. Such 'biotic resistance' recently has been questioned because of successful entry of alien species into diverse natural communities. • Data on 10,409 naturalizations of 5350 plant species over 16 sites dispersed globally show exponential distributions both for species over sites and for sites over number of species shared. These exponentials signal a statistical mechanics of species distribution, assuming two conditions. First, species and sites are equivalent, either identical ('neutral') or so complex that the chance a species is in the right place at the right time is vanishingly small ('idiosyncratic'); the range of species and sites in our data disallows a neutral explanation. Secondly, the total number of naturalizations is fixed in any era by a 'regulator'. • Previous correlation of species naturalization rates with net primary productivity over time suggests that the regulator is related to productivity. • We conclude that biotic resistance is a moving ceiling, with resistance controlled by productivity. The general observation that the majority of species occur naturally at only a few sites, and only a few species occur at many sites, now has a quantitative (exponential) character, offering the study of species' distributions a previously unavailable rigor.

Rank abundance relations in evolutionary dynamics of random replicators

We present a nonequilibrium statistical mechanics description of rank abundance relations (RAR) in random community models of ecology. Specifically, we study a multispecies replicator system with quenched random interaction matrices. We here consider symmetric interactions as well as asymmetric and antisymmetric cases. RARs are obtained analytically via a generating functional analysis, describing fixed-point states of the system in terms of a small set of order parameters, and in dependence on the symmetry or otherwise of interactions and on the productivity of the community. Our work is an extension of Tokita [Phys. Rev. Lett. 93, 178102 (2004)], where the case of symmetric interactions was considered within an equilibrium setup. The species abundance distribution in our model come out as truncated normal distributions or transformations thereof and, in some case, are similar to left-skewed distributions observed in ecology. We also discuss the interaction structure of the resulting food-web of stable species at stationarity, cases of heterogeneous cooperation pressures as well as effects of finite system size and of higher-order interactions.

Species clustering in competitive Lotka-Volterra models

We study the properties of general Lotka-Volterra models with competitive interactions. The intensity of the competition depends on the position of species in an abstract niche space through an interaction kernel. We show analytically and numerically that the properties of these models change dramatically when the Fourier transform of this kernel is not positive definite, due to a pattern-forming instability. We estimate properties of the species distributions, such as the steady number of species and their spacings, for different types of interactions, including stretched exponential and constant kernels.