No Arabic abstract
We study the extinction risk of a fragmented population residing on a network of patches coupled by migration, where the local patch dynamics include the Allee effect. We show that mixing between patches dramatically influences the populations viability. Slow migration is shown to always increase the populations global extinction risk compared to the isolated case. At fast migration, we demonstrate that synchrony between patches minimizes the populations extinction risk. Moreover, we discover a critical migration rate that maximizes the extinction risk of the population, and identify an early-warning signal when approaching this state. Our theoretical results are confirmed via the highly-efficient weighted ensemble method. Notably, our analysis can also be applied to studying switching in gene regulatory networks with multiple transcriptional states.
We study the dynamics of colonization of a territory by a stochastic population at low immigration pressure. We assume a sufficiently strong Allee effect that introduces, in deterministic theory, a large critical population size for colonization. At low immigration rates, the average pre-colonization population size is small thus invalidating the WKB approximation to the master equation. We circumvent this difficulty by deriving an exact zero-flux solution of the master equation and matching it with an approximate non-zero-flux solution of the pertinent Fokker-Planck equation in a small region around the critical population size. This procedure provides an accurate evaluation of the quasi-stationary probability distribution of population sizes in the pre-colonization state, and of the mean time to colonization, for a wide range of immigration rates. At sufficiently high immigration rates our results agree with WKB results obtained previously. At low immigration rates the results can be very different.
We investigate in detail the model of a trophic web proposed by Amaral and Meyer [Phys. Rev. Lett. 82, 652 (1999)]. We focused on small-size systems that are relevant for real biological food webs and for which the fluctuations are playing an important role. We show, using Monte Carlo simulations, that such webs can be non-viable, leading to extinction of all species in small and/or weakly coupled systems. Estimations of the extinction times and survival chances are also given. We show that before the extinction the fraction of highly-connected species (omnivores) is increasing. Viable food webs exhibit a pyramidal structure, where the density of occupied niches is higher at lower trophic levels, and moreover the occupations of adjacent levels are closely correlated. We also demonstrate that the distribution of the lengths of food chains has an exponential character and changes weakly with the parameters of the model. On the contrary, the distribution of avalanche sizes of the extinct species depends strongly on the connectedness of the web. For rather loosely connected systems we recover the power-law type of behavior with the same exponent as found in earlier studies, while for densely-connected webs the distribution is not of a power-law type.
Many populations in nature are fragmented: they consist of local populations occupying separate patches. A local population is prone to extinction due to the shot noise of birth and death processes. A migrating population from another patch can dramatically delay the extinction. What is the optimal migration rate that minimizes the extinction risk of the whole population? Here we answer this question for a connected network of model habitat patches with different carrying capacities.
We consider non-demographic noise in the form of uncertainty in the reaction step size, and reveal a dramatic effect this noise may have on the stability of self-regulating populations. Employing the reaction scheme mA->kA, but allowing, e.g., the product number k to be a-priori unknown and sampled from a given distribution, we show that such non-demographic noise can greatly reduce the populations extinction risk compared to the fixed k case. Our analysis is tested against numerical simulations, and by using empirical data of different species, we argue that certain distributions may be more evolutionary beneficial than others.
In recent years non-demographic variability has been shown to greatly affect dynamics of stochastic populations. For example, non-demographic noise in the form of a bursty reproduction process with an a-priori unknown burst size, or environmental variability in the form of time-varying reaction rates, have been separately found to dramatically impact the extinction risk of isolated populations. In this work we investigate the extinction risk of an isolated population under the combined influence of these two types of non-demographic variation. Using the so-called momentum-space WKB approach we arrive at a set of time-dependent Hamilton equations. In order to account for the explicit time dependence, we find the instanton of the time-perturbed Hamiltonian numerically, where analytical expressions are presented in particular limits using various perturbation techniques. We focus on two classes of time-varying environments: periodically-varying rates corresponding to seasonal effects, and a sudden decrease in the birth rate corresponding to a catastrophe. All our theoretical results are tested against numerical Monte Carlo simulations with time-dependent rates and also against a numerical solution of the corresponding time-dependent Hamilton equations.