No Arabic abstract
We use dynamical generating functionals to study the stability and size of communities evolving in Lotka-Volterra systems with random interaction coefficients. The size of the eco-system is not set from the beginning. Instead, we start from a set of possible species, which may undergo extinction. How many species survive depends on the properties of the interaction matrix; the size of the resulting food web at stationarity is a property of the system itself in our model, and not a control parameter as in most studies based on random matrix theory. We find that prey-predator relations enhance stability, and that variability of species interactions promotes instability. Complexity of inter-species couplings leads to reduced sizes of ecological communities. Dynamically evolved community size and stability are hence positively correlated.
In his seminal work in the 1970s Robert May suggested that there was an upper limit to the number of species that could be sustained in stable equilibrium by an ecosystem. This deduction was at odds with both intuition and the observed complexity of many natural ecosystems. The so-called stability-diversity debate ensued, and the discussion about the factors making an ecosystem stable or unstable continues to this day. We show in this work that dispersal can be a destabilising influence. To do this, we combine ideas from Alan Turings work on pattern formation with Mays random-matrix approach. We demonstrate how a stable equilibrium in a complex ecosystem with two trophic levels can become unstable with the introduction of dispersal in space. Conversely, we show that Turing instabilities can occur more easily in complex ecosystems with many species than in the case of only a few species. Our work shows that adding more details to the model of May gives rise to more ways in which an equilibrium can become unstable. Making Mays simple model more realistic is therefore unlikely to remove the upper bound on complexity.
Global dynamical behaviors of the competitive Lotka-Volterra system even in 3-dimension are not fully understood. The Lyapunov function can provide us such knowledge once it is constructed. In this paper, we construct explicitly the Lyapunov function in three examples of the competitive Lotka-Volterra system for the whole state space: (1) the general 2-dimensional case; (2) a 3-dimensional model; (3) the model of May-Leonard. The dynamics of these examples include bistable case and cyclical behavior. The first two examples are the generalized gradient system defined in the Appendixes, while the model of May-Leonard is not. Our method is helpful to understand the limit cycle problems in general 3-dimensional case.
Forty years ago, Robert May questioned a central belief in ecology by proving that sufficiently large or complex ecological networks have probability of persisting close to zero. To prove this point, he analyzed large networks in which species interact at random. However, in natural systems pairs of species have well-defined interactions (e.g., predator-prey, mutualistic or competitive). Here we extend Mays results to these relationships and find remarkable differences between predator-prey interactions, which increase stability, and mutualistic and competitive, which are destabilizing. We provide analytic stability criteria for all cases. These results have broad applicability in ecology. For example, we show that, surprisingly, the probability of stability for predator-prey networks is decreased when we impose realistic food web structure or we introduce a large preponderance of weak interactions. Similarly, stability is negatively impacted by nestedness in bipartite mutualistic networks.
The dynamics of populations is frequently subject to intrinsic noise. At the same time unknown interaction networks or rate constants can present quenched uncertainty. Existing approaches often involve repeated sampling of the quenched disorder and then running the stochastic birth-death dynamics on these samples. In this paper we take a different view, and formulate an effective jump process, representative of the ensemble of quenched interactions as a whole. Using evolutionary games with random payoff matrices as an example, we develop an algorithm to simulate this process, and we discuss diffusion approximations in the limit of weak intrinsic noise.
We study a minimal model for the growth of a phenotypically heterogeneous population of cells subject to a fluctuating environment in which they can replicate (by exploiting available resources) and modify their phenotype within a given landscape (thereby exploring novel configurations). The model displays an exploration-exploitation trade-off whose specifics depend on the statistics of the environment. Most notably, the phenotypic distribution corresponding to maximum population fitness (i.e. growth rate) requires a non-zero exploration rate when the magnitude of environmental fluctuations changes randomly over time, while a purely exploitative strategy turns out to be optimal in two-state environments, independently of the statistics of switching times. We obtain analytical insight into the limiting cases of very fast and very slow exploration rates by directly linking population growth to the features of the environment.