Density dependent dispersal decisions and the Allee effect

The decision to move between patches in the environment is among the most important life history choices an organism can make. I derive a new density dependent dispersal rule, and examine how dispersal decisions based on avoiding fitness loss associated with an Allee effect or competitive effects impact upon population dynamics in spatially structured populations with qualitatively different dynamics. I also investigate the effects of the number of patches in the system and a limit to the patch sampling time available to dispersers. Dispersing to avoid competitive pressures can destabilise otherwise stable population dynamics, and stabilise chaotic dynamics. Dispersing to avoid an Allee effect does not qualitatively change local population dynamics until eventually driving unstable populations to global extinction with a sufficiently high fitness threshold. A time limit for sampling can stabilise dynamics if dispersal is based on escaping the Allee effect, and rescue populations from global extinction. The results are sensitive to the number of patches available in the environment and suggest that dispersal to avoid an Allee effect will only arise under biologically plausible conditions, i.e. where there is a limit to the number of dispersal attempts that can be made between generations.     

Habitat selection reduces extinction of populations subject to Allee effects

Theoretical studies indicate that a single population under an Allee effect will decline to extinction if reduced below a particular threshold, but the existence of multiple local populations connected by random dispersal improves persistence of the global population. An additional process that can facilitate persistence is the existence of habitat selection by dispersers. Using analytic and simulation models of population change, I found that when habitat patches exhibiting Allee effects are connected by dispersing individuals, habitat selection by these dispersers increases the likelihood that patches persist at high densities, relative to results expected by random settlement. Populations exhibiting habitat selection also attain equilibrium more quickly than randomly dispersing populations. These effects are particularly important when Allee effects are large and more than two patches exist. Integrating habitat selection into population dynamics may help address why some studies have failed to find extinction thresholds in populations, despite well-known Allee effects in many species.     

Allee effects in stochastic populations

The Allee effect, or inverse density dependence at low population sizes, could seriously impact preservation and management of biological populations. The mounting evidence for widespread Allee effects has lately inspired theoretical studies of how Allee effects alter population dynamics. However, the recent mathematical models of Allee effects have been missing another important force prevalent at low population sizes: stochasticity. In this paper, the combination of Allee effects and stochasticity is studied using diffusion processes, a type of general stochastic population model that accommodates both demographic and environmental stochastic fluctuations. Including an Allee effect in a conventional deterministic population model typically produces an unstable equilibrium at a low population size, a critical population level below which extinction is certain. In a stochastic version of such a model, the probability of reaching a lower size a before reaching an upper size b , when considered as a function of initial population size, has an inflection point at the underlying deterministic unstable equilibrium. The inflection point represents a threshold in the probabilistic prospects for the population and is independent of the type of stochastic fluctuations in the model. In particular, models containing demographic noise alone (absent Allee effects) do not display this threshold behavior, even though demographic noise is considered an "extinction vortex". The results in this paper provide a new understanding of the interplay of stochastic and deterministic forces in ecological populations.     

Interplay between Allee effects and collective movement in metapopulations

Population growth can be positively or negatively dependent on density. Therefore, the distribution pattern of individuals in a patchy environment can greatly affect the growth of each subpopulation and thereby of the metapopulation. When population growth presents positive density‐dependence (Allee effect), the distribution pattern becomes crucial, as small populations have an increased extinction risk. The way in which individuals move between patches largely determines the distribution pattern and thereby the population dynamics. Collective movement, in particular, should be expected to increase the potential number of colonisers and therefore the probability of colonising success. Here, we use mathematical modelling (differential equations and stochastic simulations) to study how collective movement can influence metapopulation dynamics when Allee effects are at stake. The models are inspired by the two‐spotted spider mite, a phytophagous pest of recognised agricultural importance. This sub‐social mite displays trail laying/following behaviour that can provoke collective movement. Moreover, experimental evidence suggests that it is subject to Allee effects. In the first part of this study we present a single‐species population growth model incorporating Allee effects, and study its properties. In the second part, this growth model is integrated into a larger simulation model consisting of a set of interconnected patches, in which the individuals move from one patch to the other either independently or collectively. Our results show that collective movement is more advantageous than independent dispersal only when Allee effects are present and strong enough. Furthermore they provide a theoretical framework that allows the quantification of the interplay between Allee effects and collective movement.     

Interplay between local dynamics and dispersal in discrete-time metapopulation models

The effects of synchronous dispersal on discrete-time metapopulation dynamics with local (patch) dynamics of the same (compensatory or overcompensatory) or mixed (compensatory and overcompensatory) types are explored. Single-species metapopulation models behave as single-species single-patch models, whenever all local patches are governed by compensatory dynamics. Dispersal gives rise to multiple attractors with complex basin structures, whenever some local patches are under overcompensatory dynamics. In mixed systems, dispersal is capable of altering the local dynamics from compensatory to overcompensatory dynamics and vice versa. Examples are provided of metapopulation models supporting multiple attractors with intermingled basins of attraction.     

Single-species models of the Allee effect: extinction boundaries,sex ratios and mate encounters

We critically review and classify models of single-species population dynamics subject to the demographic Allee effect with emphasis on non-spatial, deterministic approach. Inclusion of spatial movement and stochastic phenomena does not substantially change the behaviour; stochasticity only "blurs" step-like character of the Allee effect into a sigmoidal form. The outcome of all non-spatial, deterministic models is either unconditional extinction, extinction-survival scenario (ES), or unconditional survival. Three major model classes are recognized: (1) one-dimensional heuristic models, (2) one-dimensional models with mating probability and fixed sex ratio, and (3) two-sex models with variable adult sex ratio. Each class is characterized by the shape of extinction boundary which separates extinction from survival in the ES scenario. The latter two classes may give better predictions of extinction thresholds than heuristic models but require specific information and are data intensive. In one-dimensional models with fixed sex ratio, population cannot survive if density/number of males decreases below some threshold while there is no such restriction on females. Individual-based models seem to be most capable of explaining mechanisms leading to the Allee effect.     

Stochastic Sensitivity Analysis of Noise-Induced Extinction in the Ricker Model with Delay and Allee Effect

A susceptibility of population systems to the random noise is studied on the base of the conceptual Ricker-type model taking into account the delay and Allee effect. This two-dimensional discrete model exhibits the persistence in the form of equilibria, discrete cycles, closed invariant curves, and chaotic attractors. It is shown how the Allee effect constrains the persistence zones with borders defined by crisis bifurcations. We study the role of random noise on the contraction and destruction of these zones. This phenomenon of the noise-induced extinction is investigated with the help of direct numerical simulations and semi-analytical approach based on the stochastic sensitivity functions. Stochastic transitions from the persistence regimes to the extinction are studied by the analysis of the mutual arrangement of the basins of attraction and confidence domains.     

Allee effects limit population viability of an annual plant

ABSTRACT Allee effects may be experienced by plants when populations are too small or isolated to receive sufficient pollinator services to replace themselves. This article reports experimental data from an annual herb, Clarkia concinna, documenting that small patches suffered reproductive failure due to lack of effective pollination when critical thresholds of isolation were exceeded. In contrast, sufficiently large patches attracted pollinators regardless of their degree of isolation. These data accord with data on patch extinctions showing that small and isolated patches have a higher extinction rate than do large patches and with observations showing chronically low reproductive success in such patches prior to extinction. While not conclusively demonstrating that Allee effects cause extinction in small and isolated patches, the data are suggestive. Although threshold effects have been postulated in several mathematical models of population viability, this is the first report of data from natural populations that display the occurrence of such thresholds. These results have implications for the management of endangered plants, which often are restricted to isolated, small populations, as well as suggesting a potential limit to spatial spread in plant populations dependent on animal vectors for reproduction.     

Generalizing Levins metapopulation model in explicit space: models of intermediate complexity

A recent study [Harding and McNamara, 2002. A unifying framework for metapopulation dynamics. Am. Nat. 160, 173-185] presented a unifying framework for the classic Levins metapopulation model by incorporating several realistic biological processes, such as the Allee effect, the Rescue effect and the Anti-rescue effect, via appropriate modifications of the two basic functions of colonization and extinction rates. Here we embed these model extensions on a spatially explicit framework. We consider population dynamics on a regular grid, each site of which represents a patch that is either occupied or empty, and with spatial coupling by neighborhood dispersal. While broad qualitative similarities exist between the spatially explicit models and their spatially implicit (mean-field) counterparts, there are also important differences that result from the details of local processes. Because of localized dispersal, spatial correlation develops among the dynamics of neighboring populations that decays with distance between patches. The extent of this correlation at equilibrium differs among the metapopulation types, depending on which processes prevail in the colonization and extinction dynamics. These differences among dynamical processes become manifest in the spatial pattern and distribution of “clusters” of occupied patches. Moreover, metapopulation dynamics along a smooth gradient of habitat availability show significant differences in the spatial pattern at the range limit. The relevance of these results to the dynamics of disease spread in metapopulations is discussed.     

The downward spiral: eco‐evolutionary feedback loops lead to the emergence of ‘elastic’ ranges

In times of severe environmental changes and resulting shifts in the geographical distribution of animal and plant species it is crucial to unravel the mechanisms responsible for the dynamics of species’ ranges. Without such a mechanistic understanding, reliable projections of future species distributions are difficult to derive. Species’ ranges may be highly dynamic. One particularly interesting phenomenon is range contraction following a period of expansion, referred to as ‘elastic’ behaviour. It has been proposed that this phenomenon occurs in habitat gradients, which are characterized by a negative cline in selection for dispersal from the range core towards the margin, as one may find, for example, with increasing patch isolation. Using individual‐based simulations and numerical analyses we show that Allee effects are an important determinant of range border elasticity. If only intra‐specific processes are considered, Allee effects are even a necessary condition for ranges to exhibit elastic behavior. The eco‐evolutionary interplay between dispersal evolution, Allee effects and habitat isolation leads to lower colonization probability and higher local extinction risk after range expansions, which result in an increasing amount of marginal sink patches and consequently, range contraction. We also demonstrate that the nature of the gradient is crucial for range elasticity. Gradients which do not select for lower dispersal at the margin than in the core (especially gradients in patch size, demographic stochasticity and extinction rate) do not lead to elastic range behavior. Thus, we predict that range contractions are likely to occur after periods of expansion for species living in gradients of increasing patch isolation, which suffer from Allee effects.     

Allee effects and resilience in stochastic populations

Allee effects, or positive functional relationships between a population’s density (or size) and its per unit abundance growth rate, are now considered to be a widespread if not common influence on the growth of ecological populations. Here we analyze how stochasticity and Allee effects combine to impact population persistence. We compare the deterministic and stochastic properties of four models: a logistic model (without Allee effects), and three versions of the original model of Allee effects proposed by Vito Volterra representing a weak Allee effect, a strong Allee effect, and a strong Allee effect with immigration. We employ the diffusion process approach for modeling single-species populations, and we focus on the properties of stationary distributions and of the mean first passage times. We show that stochasticity amplifies the risks arising from Allee effects, mainly by prolonging the amount of time a population spends at low abundance levels. Even weak Allee effects become consequential when the ubiquitous stochastic forces affecting natural populations are accounted for in population models. Although current concepts of ecological resilience are bound up in the properties of deterministic basins of attraction, a complete understanding of alternative stable states in ecological systems must include stochasticity.     

Allee threshold and extinction threshold for spatially explicit metapopulation dynamics with Allee effects

In this paper, we investigate a spatially explicit metapopulation model with Allee effects. We refer to the patch occupancy model introduced by Levins (Bull Entomol Soc Am 15:237–240, 1969) as a spatially implicit metapopulation model, i.e., each local patch is either occupied or vacant and a vacant patch can be recolonized by a randomly chosen occupied patch from anywhere in the metapopulation. When we transform the model into a spatially explicit one by using a lattice model, the obtained model becomes theoretically equivalent to a “lattice logistic model” or a “basic contact process”. One of the most popular or standard metapopulation models with Allee effects, developed by Amarasekare (Am Nat 152:298–302, 1998), supposes that those effects are introduced formally by means of a logistic equation. However, it is easier to understand the ecological meaning of associating Allee effects with this model if we suppose that only the logistic colonization term directly suffers from Allee effects. The resulting model is also well defined, and therefore we can naturally examine it by Monte Carlo simulation and by doublet and triplet decoupling approximation. We then obtain the following specific features of one-dimensional lattice space: (1) the metapopulation as a whole does not have an Allee threshold for initial population size even when each local population follows the Allee effects; and (2) a metapopulation goes extinct when the extinction rate of a local population is lower than that in the spatially implicit model. The real ecological metapopulation lies between two extremes: completely mixing interactions between patches on the one hand and, on the other, nearest neighboring interactions with only two nearest neighbors. Thus, it is important to identify the metapopulation structure when we consider the problems of invasion species such as establishment or the speed of expansion.     

Determinants of extinction in fragmented plant populations: Crepis sancta (asteraceae) in urban environments

Local populations are subject to recurrent extinctions, and small populations are particularly prone to extinction. Both demographic (stochasticity and the Allee effect) and genetic factors (drift load and inbreeding depression) potentially affect extinction. In fragmented populations, regular dispersal may boost population sizes (demographic rescue effect) or/and reduce the local inbreeding level and genetic drift (genetic rescue effect), which can affect extinction risks. We studied extinction processes in highly fragmented populations of the common species Crepis sancta (Asteraceae) in urban habitats exhibiting a rapid turnover of patches. A four-year demographic monitoring survey and microsatellite genotyping of individuals allowed us to study the determinants of extinction. We documented a low genetic structure and an absence of inbreeding (estimated by multilocus heterozygosity), which suggest that genetic factors were not a major cause of patch extinction. On the contrary, local population size was the main factor in extinction, whereas connectivity was shown to decrease patch extinction, which we interpreted as a demographic rescue effect that was likely due to better pollination services for reproduction. This coupling of demographic and genetic tools highlighted the importance of dispersal in local patch extinctions of small fragmented populations connected by gene flow.     

Establishing a beachhead: a stochastic population model with an Allee effect applied to species invasion

We formulated a spatially explicit stochastic population model with an Allee effect in order to explore how invasive species may become established. In our model, we varied the degree of migration between local populations and used an Allee effect with variable birth and death rates. Because of the stochastic component, population sizes below the Allee effect threshold may still have a positive probability for successful invasion. The larger the network of populations, the greater the probability of an invasion occurring when initial population sizes are close to or above the Allee threshold. Furthermore, if migration rates are low, one or more than one patch may be successfully invaded, while if migration rates are high all patches are invaded.     

Evolutionary suicide and evolution of dispersal in structured metapopulations

 We study the evolution of dispersal in a structured metapopulation model. The metapopulation consists of a large (infinite) number of local populations living in patches of habitable environment. Dispersal between patches is modelled by a disperser pool and individuals in transit between patches are exposed to a risk of mortality. Occasionally, local catastrophes eradicate a local population: all individuals in the affected patch die, yet the patch remains habitable. We prove that, in the absence of catastrophes, the strategy not to migrate is evolutionarily stable. Under a given set of environmental conditions, a metapopulation may be viable and yet selection may favor dispersal rates that drive the metapopulation to extinction. This phenomenon is known as evolutionary suicide. We show that in our model evolutionary suicide can occur for catastrophe rates that increase with decreasing local population size. Evolutionary suicide can also happen for constant catastrophe rates, if local growth within patches shows an Allee effect. We study the evolutionary bifurcation towards evolutionary suicide and show that a discontinuous transition to extinction is a necessary condition for evolutionary suicide to occur. In other words, if population size smoothly approaches zero at a boundary of viability in parameter space, this boundary is evolutionarily repelling and no suicide can occur. Received: 10 November 2000 / Revised version: 13 February 2002 / Published online: 17 July 2002     

Non-equilibria in small metapopulations: comparing the deterministic Levins model with its stochastic counterpart

In this paper, we examine, for small metapopulations, the stochastic analog of the classical Levins metapopulation model. We study its basic model output, the expected time to metapopulation extinction, for systems which are brought out of equilibrium by imposing sudden changes in patch number and the colonization and extinction parameters. We find that the expected metapopulation extinction time shows different behavior from the relaxation time of the original, deterministic, Levins model. This relaxation time is therefore limited in value for predicting the behavior of the stochastic model. However, predictions about the extinction time for deterministically unviable cases remain qualitatively the same. Our results further suggest that, if we want to counteract the effects of habitat loss or increased dispersal resistance, the optimal conservation strategy is not to restore the original situation, that is, to create habitat or decrease resistance against dispersal. As long as the costs for different management options are not too dissimilar, it is better to improve the quality of the remaining habitat in order to decrease the local extinction rate.     

Predator-prey interactions in a nonequilibrium context: the metapopulation approach to modeling "hide-and-seek" dynamics in a spatially explicit tri-trophic system

Predators and prey are usually heterogeneously distributed in space so that the ability of the predators to respond to the distribution of their prey may have a profound influence on the stability and persistence of a predator‐prey system. A special type of dynamics is “hide‐and‐seek” characterized by a high turnover rate of local populations of prey and predators, because once the predators have found a patch of prey they quickly overexploit it, whereupon the starving predators either should move to better places or die. Continued persistence of prey and predators thus hinges on a long‐term balance between local extinctions and founding of new subpopulations. The colonization rate depends on the rate of emigration from occupied patches and the likelihood of successfully arriving at a suitable new patch, while extinction rate depends on the local population dynamics. Since extinctions and colonizations are both discrete probabilistic events, these phenomena are most adequately modeled by means of a stochastic model. In order to demonstrate the qualitative differences between a deterministic and stochastic approach to population dynamics, a spatially explicit tritrophic predator‐prey model is developed in a deterministic and a stochastic version. The model is parameterized using data for the two‐spotted spider mite (Tetranychus urticae) and the phytoseiid mite predator Phytoseiulus persimilis inhabiting greenhouse cucumbers.
Simulations show that the deterministic and stochastic approaches yield different results. The deterministic version predicts that the populations will exhibit violent fluctuations, implying that the system is fundamentally unstable. In contrast, the stochastic version predicts that the two species will be able to coexist in spite of frequent local extinctions of both species, provided the system consists of a sufficiently large number of local populations. This finding is in agreement with experimental results. It is therefore concluded that demographic stochasticity in combination with dispersal is capable of producing and maintaining sufficient asynchrony between local populations to ensure long‐term regional (metapopulation) persistence.     

Allee效应对具有生境恢复的集合种群的影响

Allee效应与种群的灭绝密切相关,其研究对生态保护和管理至关重要。Allee效应对物种续存是潜在的干扰因素,濒危物种更容易受其影响,可能会增加生存于生境破碎化斑块的濒危物种的死亡风险,因此研究Allee效应对种群的动态和续存的影响是必要的。从包含由生物有机体对环境的修复产生的Allee效应的集合种群模型出发,引入由其他机制形成的Allee效应,建立了常微分动力系统模型和基于网格模型的元胞自动机模型。通过理论分析和计算机模拟表明:(1)强Allee效应不利于具有生境恢复的集合种群的续存;(2)生境恢复有利于种群续存;(3)局部扩散影响了集合种群的空间结构、动态行为和稳定性,生境斑块之间的局部作用将会减缓或消除集合种群的Allee效应,有利于集合种群的续存。     

Colonization–extinction dynamics of epixylic lichens along a decay gradient in a dynamic landscape

Metapopulation models are often used for understanding and predicting species dynamics in fragmented landscapes. Several models have been proposed depending on e.g. the relative importance of patch dynamics on the metapopulation dynamics. Dead wood is a dynamic substrate patch, and species that are confined to such patches have experienced a high degree of habitat loss in managed forests. Little is, however, known about how the population dynamics of epixylic species are affected by the fast dynamics of their substrate patches. We quantified the effect of local patch conditions and metapopulation processes on colonizations and extinctions of epixylic lichen species in a managed boreal forest landscape. This was done by twice surveying seven lichen metapopulations on 293 stumps in 30 stands of ages covering the duration of the dynamic patches (stumps). We also investigated the relative importance of local stochastic extinctions from stumps that remained available, and deterministic extinctions due to stump surface disappearance. We found importance of a decay gradient, surrounding metapopulation size, and local population sizes, in driving the colonization–extinction dynamics of epixylic lichens. The species were sorted along the stump decay gradient. Increasing surrounding metapopulation size was associated with increased colonization rates, and increasing local population size decreased lichen extinction rates. Finally, both local stochastic extinctions and deterministic extinctions due to patch disappearance occur, confirming that the long‐term persistence of epixylic lichens depends on colonization rates that compensate for stochastic population extinctions as well as deterministic extinctions.     

Contributions of high- and low-quality patches to a metapopulation with stochastic disturbance

Studies of time-invariant matrix metapopulation models indicate that metapopulation growth rate is usually more sensitive to the vital rates of individuals in high-quality (i.e., good) patches than in low-quality (i.e., bad) patches. This suggests that, given a choice, management efforts should focus on good rather than bad patches. Here, we examine the sensitivity of metapopulation growth rate for a two-patch matrix metapopulation model with and without stochastic disturbance and found cases where managers can more efficiently increase metapopulation growth rate by focusing efforts on the bad patch. In our model, net reproductive rate differs between the two patches so that in the absence of dispersal, one patch is high quality and the other low quality. Disturbance, when present, reduces net reproductive rate with equal frequency and intensity in both patches. The stochastic disturbance model gives qualitatively similar results to the deterministic model. In most cases, metapopulation growth rate was elastic to changes in net reproductive rate of individuals in the good patch than the bad patch. However, when the majority of individuals are located in the bad patch, metapopulation growth rate can be most elastic to net reproductive rate in the bad patch. We expand the model to include two stages and parameterize the patches using data for the softshell clam, Mya arenaria. With a two-stage demographic model, the elasticities of metapopulation growth rate to parameters in the bad patch increase, while elasticities to the same parameters in the good patch decrease. Metapopulation growth rate is most elastic to adult survival in the population of the good patch for all scenarios we examine. If the majority of the metapopulation is located in the bad patch, the elasticity to parameters of that population increase but do not surpass elasticity to parameters in the good patch. This model can be expanded to include additional patches, multiple stages, stochastic dispersal, and complex demography.