Non-neutral theory of biodiversity

We present a non-neutral stochastic model for the dynamics taking place in a meta-community ecosystems in presence of migration. The model provides a framework for describing the emergence of multiple ecological scenarios and behaves in two extreme limits either as the unified neutral theory of biodiversity or as the Bak-Sneppen model. Interestingly, the model shows a condensation phase transition where one species becomes the dominant one, the diversity in the ecosystems is strongly reduced and the ecosystem is non-stationary. This phase transition extend the principle of competitive exclusion to open ecosystems and might be relevant for the study of the impact of invasive species in native ecologies.

💡 Research Summary

The paper introduces a stochastic, non‑neutral model for meta‑community dynamics that explicitly incorporates species‑specific fitness and migration. The authors begin by reviewing the limitations of the Unified Neutral Theory of Biodiversity (UNTB), which assumes ecological equivalence among individuals, and the Bak‑Sneppen model, which captures evolutionary avalanches but lacks explicit migration. Their new framework bridges these two extremes by defining a meta‑community composed of L local patches, each containing N individuals. At each time step two processes occur: (1) internal reproduction, where an individual of species i is chosen with probability proportional to its fitness w_i, and (2) migration, where an individual drawn from a global pool replaces a local individual with rate ε. By varying ε, the model interpolates between a pure neutral regime (ε → ∞) and a pure selection‑driven regime (ε → 0) that reproduces Bak‑Sneppen‑type criticality.

A central result is the discovery of a condensation phase transition. When the product of the migration rate ε and the tail of the fitness distribution exceeds a critical threshold θ, a single species with sufficiently high fitness dominates the entire meta‑community. In this “condensed” phase, species richness collapses, the Shannon diversity index drops sharply, and the system becomes non‑stationary, exhibiting long‑range temporal correlations rather than the exponential relaxation typical of neutral dynamics. Analytically, the authors derive a master equation for the probability distribution of species abundances and identify two classes of solutions: a stationary, neutral‑like distribution and a non‑stationary, condensed solution where one species occupies a macroscopic fraction of all individuals. The critical point depends on the shape of the fitness distribution; for power‑law tails P(w) ∝ w^‑α, smaller α (heavier tails) lower the critical ε, making condensation more likely even with modest migration.

Through extensive Monte‑Carlo simulations (N = 10^3, L = 100, ε ranging from 10^‑4 to 1, α from 1.5 to 3), the authors verify the analytical predictions. Prior to condensation, species abundance follows log‑normal or exponential tails, and the occupancy‑abundance relationship is weakly positive. After condensation, a single spike appears in the abundance distribution, representing the dominant species, while all others are relegated to negligible frequencies. Time‑correlation functions shift from exponential decay with characteristic time ~1/ε to power‑law decay (C(t) ∝ t^‑β), indicating the emergence of memory effects.

The ecological implications are highlighted by linking the condensation transition to invasive species dynamics. An introduced organism with high fitness can, through a few migration events, trigger a rapid takeover of an island or fragmented habitat, mirroring the model’s condensed phase. Conversely, managing migration pathways (e.g., physical barriers, quarantine measures) can keep ε below the critical value, preserving biodiversity. The authors suggest that the model provides a quantitative framework for assessing invasion risk and designing conservation strategies in open ecosystems.

Finally, the paper acknowledges limitations: fitness is treated as static, interspecific interactions such as predation or mutualism are omitted, and spatial structure is reduced to a mean‑field coupling via ε. Future work is proposed to incorporate evolving fitness landscapes, multi‑trophic networks, and explicit spatial networks of patches. Despite these simplifications, the non‑neutral theory presented successfully unifies neutral and selection‑driven perspectives, offering a versatile tool for understanding biodiversity patterns, phase‑like transitions in ecological communities, and the impact of migration on ecosystem stability.