The dynamics of two competing species in a finite size community is one of the most studied problems in population genetics and community ecology. Stochastic fluctuations lead, inevitably, to the extinction of one of the species, but the relevant timescale depends on the underlying dynamics. The persistence time of the community has been calculated for neutral models, where the only drive of the system is drift (demographic stochasticity) and for models with strong selection. Following recent analyses that stress the importance of environmental stochasticity in empirical systems, we present here a general theory of persistence time of two-species community where drift, environmental variations and time independent selective advantage are all taken into account.

The dynamics of two competing species in a finite size community is one of the most studied problems in population genetics and community ecology. Stochastic fluctuations lead, inevitably, to the extinction of one of the species, but the relevant timescale depends on the underlying dynamics. The persistence time of the community has been calculated for neutral models, where the only drive of the system is drift (demographic stochasticity) and for models with strong selection. Following recent analyses that stress the importance of environmental stochasticity in empirical systems, we present here a general theory of persistence time of two-species community where drift, environmental variations and time independent selective advantage are all taken into account.

A fundamental problem in the fields of population genetics, evolution, and community ecology, is the fate of a single mutant, or invader, introduced in a finite population of wild types. For a fixed-size community of N individuals, with Markovian, zero-sum dynamics driven by stochastic birth-death events, the mutant population eventually reaches either fixation or extinction. The classical analysis, provided by Kimura and his coworkers, is focused on the neutral case, [where the dynamics is only due to demographic stochasticity (drift)], and on time-independent selective forces (deleterious/beneficial mutation). However, both theoretical arguments and empirical analyses suggest that in many cases the selective forces fluctuate in time (temporal environmental stochasticity). Here we consider a generic model for a system with demographic noise and fluctuating selection. Our system is characterized by the time-averaged (log)-fitness s and zero-mean fitness fluctuations. These fluctuations, in turn, are parameterized by their amplitude γ and their correlation time δ. We provide asymptotic (large N) formulas for the chance of fixation, the mean time to fixation and the mean time to absorption. Our expressions interpolate correctly between the constant selection limit γ → 0 and the time-averaged neutral case s=0.

Competition is the main driver of population dynamics, which shapes the genetic composition of populations and the assembly of ecological communities. Neutral models assume that all the individuals are equivalent and that the dynamics is governed by demographic (shot) noise, with a steady state species abundance distribution (SAD) that reflects a mutation-extinction equilibrium. Recently, many empirical and theoretical studies emphasized the importance of environmental variations that affect coherently the relative fitness of entire populations. Here we consider two generic time-averaged neutral models; in both the relative fitness of each species fluctuates independently in time but its mean is zero. The first (model A) describes a system with local competition and linear fitness dependence of the birth-death rates, while in the second (model B) the competition is global and the fitness dependence is nonlinear. Due to this nonlinearity, model B admits a noise-induced stabilization mechanism that facilitates the invasion of new mutants. A self-consistent mean-field approach is used to reduce the multispecies problem to two-species dynamics, and the large-N asymptotics of the emerging set of Fokker-Planck equations is presented and solved. Our analytic expressions are shown to fit the SADs obtained from extensive Monte Carlo simulations and from numerical solutions of the corresponding master equations.

Neutral models, in which individual agents with equal fitness undergo a birth-death-mutation process, are very popular in population genetics and community ecology. Usually these models are applied to populations and communities with spatial structure, but the analytic results presented so far are limited to well-mixed or mainland-island scenarios. Here we combine analytic results and numerics to obtain an approximate solution for the species abundance distribution and the species richness for the neutral model on continuous landscape. We show how the regional diversity increases when the recruitment length decreases and the spatial segregation of species grows. Our results are supported by extensive numerical simulations and allow one to probe the numerically inaccessible regime of large-scale systems with extremely small mutation/speciation rates. Model predictions are compared with the findings of recent large-scale surveys of tropical trees across the Amazon basin, yielding new bounds for the species richness (between 13100 and 15000) and the number of singleton species (between 455 and 690).