Time-resolved extinction rates of stochastic populations

Extinction dynamics from metastable coexistences in an evolutionary game

Deterministic evolutionary game dynamics can lead to stable coexistences of different types. Stochasticity, however, drives the loss of such coexistences. This extinction is usually accompanied by population size fluctuations. We investigate the most probable extinction trajectory under such fluctuations by mapping a stochastic evolutionary model to a problem of classical mechanics using the Wentzel-Kramers-Brillouin (WKB) approximation. Our results show that more abundant types in a coexistence may be more likely to go extinct first, in good agreement with previous results. The distance between the coexistence and extinction points is not a good predictor of extinction either. Instead, the WKB method correctly predicts the type going extinct first.

Extinction of oscillating populations

Established populations often exhibit oscillations in their sizes that, in the deterministic theory, correspond to a limit cycle in the space of population sizes. If a population is isolated, the intrinsic stochasticity of elemental processes can ultimately bring it to extinction. Here we study extinction of oscillating populations in a stochastic version of the Rosenzweig-MacArthur predator-prey model. To this end we develop a WKB (Wentzel, Kramers and Brillouin) approximation to the master equation, employing the characteristic population size as the large parameter. Similar WKB theories have been developed previously in the context of population extinction from an attracting multipopulation fixed point. We evaluate the extinction rates and find the most probable paths to extinction from the limit cycle by applying Floquet theory to the dynamics of an effective four-dimensional WKB Hamiltonian. We show that the entropic barriers to extinction change in a nonanalytic way as the system passes through the Hopf bifurcation. We also study the subleading pre-exponential factors of the WKB approximation.

Two-strain competition in quasineutral stochastic disease dynamics

We develop a perturbation method for studying quasineutral competition in a broad class of stochastic competition models and apply it to the analysis of fixation of competing strains in two epidemic models. The first model is a two-strain generalization of the stochastic susceptible-infected-susceptible (SIS) model. Here we extend previous results due to Parsons and Quince [Theor. Popul. Biol. 72, 468 (2007)], Parsons et al. [Theor. Popul. Biol. 74, 302 (2008)], and Lin, Kim, and Doering [J. Stat. Phys. 148, 646 (2012)]. The second model, a two-strain generalization of the stochastic susceptible-infected-recovered (SIR) model with population turnover, has not been studied previously. In each of the two models, when the basic reproduction numbers of the two strains are identical, a system with an infinite population size approaches a point on the deterministic coexistence line (CL): a straight line of fixed points in the phase space of subpopulation sizes. Shot noise drives one of the strain populations to fixation, and the other to extinction, on a time scale proportional to the total population size. Our perturbation method explicitly tracks the dynamics of the probability distribution of the subpopulations in the vicinity of the CL. We argue that, whereas the slow strain has a competitive advantage for mathematically "typical" initial conditions, it is the fast strain that is more likely to win in the important situation when a few infectives of both strains are introduced into a susceptible population.

Non-equilibrium physics and evolution--adaptation, extinction, and ecology: a key issues review

Evolutionary dynamics in nature constitute an immensely complex non-equilibrium process. We review the application of physical models of evolution, by focusing on adaptation, extinction, and ecology. In each case, we examine key concepts by working through examples. Adaptation is discussed in the context of bacterial evolution, with a view toward the relationship between growth rates, mutation rates, selection strength, and environmental changes. Extinction dynamics for an isolated population are reviewed, with emphasis on the relation between timescales of extinction, population size, and temporally correlated noise. Ecological models are discussed by focusing on the effect of spatial interspecies interactions on diversity. Connections between physical processes--such as diffusion, turbulence, and localization--and evolutionary phenomena are highlighted.

Multiple extinction routes in stochastic population models

Isolated populations ultimately go extinct because of the intrinsic noise of elementary processes. In multipopulation systems extinction of a population may occur via more than one route. We investigate this generic situation in a simple predator-prey (or infected-susceptible) model. The predator and prey populations may coexist for a long time, but ultimately both go extinct. In the first extinction route the predators go extinct first, whereas the prey thrive for a long time and then also go extinct. In the second route the prey go extinct first, causing a rapid extinction of the predators. Assuming large subpopulation sizes in the coexistence state, we compare the probabilities of each of the two extinction routes and predict the most likely path of the subpopulations to extinction. We also suggest an effective three-state master equation for the probabilities to observe the coexistence state, the predator-free state, and the empty state.

Switching between phenotypes and population extinction

Many types of bacteria can survive under stress by switching stochastically between two different phenotypes: the "normals" who multiply fast, but are vulnerable to stress, and the "persisters" who hardly multiply, but are resilient to stress. Previous theoretical studies of such bacterial populations have focused on the fitness: the asymptotic rate of unbounded growth of the population. Yet for an isolated population of established (and not very large) size, a more relevant measure may be the population extinction risk due to the interplay of adverse extrinsic variations and intrinsic noise of birth, death and switching processes. Applying a WKB approximation to the pertinent master equation of such a two-population system, we quantify the extinction risk, and find the most likely path to extinction under both favorable and adverse conditions. Analytical results are obtained both in the biologically relevant regime when the switching is rare compared with the birth and death processes, and in the opposite regime of frequent switching. We show that rare switches are most beneficial in reducing the extinction risk.

Extinction rates of established spatial populations

This paper deals with extinction of an isolated population caused by intrinsic noise. We model the population dynamics in a "refuge" as a Markov process which involves births and deaths on discrete lattice sites and random migrations between neighboring sites. In extinction scenario I, the zero population size is a repelling fixed point of the on-site deterministic dynamics. In extinction scenario II, the zero population size is an attracting fixed point, corresponding to what is known in ecology as the Allee effect. Assuming a large population size, we develop a WKB (Wentzel-Kramers-Brillouin) approximation to the master equation. The resulting Hamilton's equations encode the most probable path of the population toward extinction and the mean time to extinction. In the fast-migration limit these equations coincide, up to a canonical transformation, with those obtained, in a different way, by Elgart and Kamenev [Phys. Rev. E 70, 041106 (2004)]. We classify possible regimes of population extinction with and without an Allee effect and for different types of refuge, and solve several examples analytically and numerically. For a very strong Allee effect, the extinction problem can be mapped into the overdamped limit of the theory of homogeneous nucleation due to Langer [Ann. Phys. (NY) 54, 258 (1969)]. In this regime, and for very long systems, we predict an optimal refuge size that maximizes the mean time to extinction.