The ideal free distribution as an evolutionarily stable strategy

We examine the evolutionary stability of strategies for dispersal in heterogeneous patchy environments or for switching between discrete states (e.g. defended and undefended) in the context of models for population dynamics or species interactions in either continuous or discrete time. There have been a number of theoretical studies that support the view that in spatially heterogeneous but temporally constant environments there will be selection against unconditional, i.e. random, dispersal, but there may be selection for certain types of dispersal that are conditional in the sense that dispersal rates depend on environmental factors. A particular type of dispersal strategy that has been shown to be evolutionarily stable in some settings is balanced dispersal, in which the equilibrium densities of organisms on each patch are the same whether there is dispersal or not. Balanced dispersal leads to a population distribution that is ideal free in the sense that at equilibrium all individuals have the same fitness and there is no net movement of individuals between patches or states. We find that under rather general assumptions about the underlying population dynamics or species interactions, only such ideal free strategies can be evolutionarily stable. Under somewhat more restrictive assumptions (but still in considerable generality), we show that ideal free strategies are indeed evolutionarily stable. Our main mathematical approach is invasibility analysis using methods from the theory of ordinary differential equations and nonnegative matrices. Our analysis unifies and extends previous results on the evolutionary stability of dispersal or state-switching strategies.     

Evolutionary stability of ideal free dispersal strategies in patchy environments

A central question in the study of the evolution of dispersal is what kind of dispersal strategies are evolutionarily stable. Hastings (Theor Pop Biol 24:244-251, 1983) showed that among unconditional dispersal strategies in a spatially heterogeneous but temporally constant environment, the dispersal strategy with no movement is convergent stable. McPeek and Holt's (Am Nat 140:1010-1027, 1992) work suggested that among conditional dispersal strategies in a spatially heterogeneous but temporally constant environment, an ideal free dispersal strategy, which results in the ideal free distribution for a single species at equilibrium, is evolutionarily stable. We use continuous-time and discrete-space models to determine when the dispersal strategy with no movement is evolutionarily stable and when an ideal free dispersal strategy is evolutionarily stable, both in a spatially heterogeneous but temporally constant environment.     

Evolutionary Convergence to Ideal Free Dispersal Strategies and Coexistence

We study a two species competition model in which the species have the same population dynamics but different dispersal strategies and show how these dispersal strategies evolve. We introduce a general dispersal strategy which can result in the ideal free distributions of both competing species at equilibrium and generalize the result of Averill et al. (2011). We further investigate the convergent stability of this ideal free dispersal strategy by varying random dispersal rates, advection rates, or both of these two parameters simultaneously. For monotone resource functions, our analysis reveals that among two similar dispersal strategies, selection generally prefers the strategy which is closer to the ideal free dispersal strategy. For nonmonotone resource functions, our findings suggest that there may exist some dispersal strategies which are not ideal free, but could be locally evolutionarily stable and/or convergent stable, and allow for the coexistence of more than one species.     

Evolution of unconditional dispersal in periodic environments

Organisms modulate their fitness in heterogeneous environments by dispersing. Prior work shows that there is selection against 'unconditional' dispersal in spatially heterogeneous environments. 'Unconditional' means individuals disperse at a rate independent of their location. We prove that if within-patch fitness varies spatially and between two values temporally, then there is selection for unconditional dispersal: any evolutionarily stable strategy (ESS) or evolutionarily stable coalition (ESC) includes a dispersive phenotype. Moreover, at this ESS or ESC, there is at least one sink patch (i.e. geometric mean of fitness less than one) and no sources patches (i.e. geometric mean of fitness greater than one). These results coupled with simulations suggest that spatial-temporal heterogeneity is due to abiotic forcing result in either an ESS with a dispersive phenotype or an ESC with sedentary and dispersive phenotypes. In contrast, the spatial-temporal heterogeneity due to biotic interactions can select for higher dispersal rates that ultimately spatially synchronize population dynamics.     

Evolutionary stability of ideal free nonlocal dispersal

We study the evolutionary stability of nonlocal dispersal strategies that can produce ideal free population distributions, that is, distributions where all individuals have equal fitness and there is no net movement of individuals at equilibrium. We find that the property of producing ideal free distributions is necessary and often sufficient for evolutionary stability. Our results extend those already developed for discrete diffusion models on finite patch networks to the case of nonlocal dispersal models based on integrodifferential equations. The analysis is based on the use of comparison methods and the construction of sub- and supersolutions.     

Evolutionary stability of ideal free nonlocal dispersal

We study the evolutionary stability of nonlocal dispersal strategies that can produce ideal free population distributions, that is, distributions where all individuals have equal fitness and there is no net movement of individuals at equilibrium. We find that the property of producing ideal free distributions is necessary and often sufficient for evolutionary stability. Our results extend those already developed for discrete diffusion models on finite patch networks to the case of nonlocal dispersal models based on integrodifferential equations. The analysis is based on the use of comparison methods and the construction of sub- and supersolutions.     

Equilibrium structure and stability in a frequency-dependent,two-population diploid model

We investigate the equilibrium structure for an evolutionary genetic model in discrete time involving two monoecious populations subject to intraspecific and interspecific random pairwise interactions. A characterization for local stability of an equilibrium is found, related to the proximity of this equilibrium with evolutionarily stable strategies (ESS). This extends to a multi-population framework a principle initially proposed for single populations, which states that the mean population strategy at a locally stable equilibrium is as close as possible to an ESS.     

Adaptive diversification of germination strategies.

Evolution of the germination rate (the proportion of newly produced and dormant seeds that germinates every year) of annual plants is investigated, when the environment is temporally stochastic and spatially heterogeneous. The environment consists of two habitats with synchronous stochastic variation in the annual yield and permanent difference in constant seed survival rates. Density dependence operates within the habitats, which are connected via restricted seed dispersal. We find that instead of a single common evolutionarily stable strategy the coexistence of several germination strategies is possible and that in an initially monomorphic population evolutionary branching may occur. During evolutionary branching the population undergoes disruptive selection and splits into two branches of different lineages that converge to the evolutionarily stable coalition of different germination strategies. It is shown that spatial heterogeneity and restricted dispersal are essential for evolutionary branching. Disruptive selection on the germination rate presents yet another possibility for parapatric speciation.     

On several conjectures from evolution of dispersal

We address several conjectures raised in Cantrell et al. [Evolution of dispersal and ideal free distribution, Math. Biosci. Eng. 7 (2010), pp. 17-36 [ 9 ]] concerning the dynamics of a diffusion-advection-competition model for two competing species. A conditional dispersal strategy, which results in the ideal free distribution of a single population at equilibrium, was found in Cantrell et al. [ 9 ]. It was shown in [ 9 ] that this special dispersal strategy is a local evolutionarily stable strategy (ESS) when the random diffusion rates of the two species are equal, and here we show that it is a global ESS for arbitrary random diffusion rates. The conditions in [ 9 ] for the coexistence of two species are substantially improved. Finally, we show that this special dispersal strategy is not globally convergent stable for certain resource functions, in contrast with the result from [ 9 ], which roughly says that this dispersal strategy is globally convergent stable for any monotone resource function.     

The ideal free distribution as an evolutionarily stable state in density‐dependent population games

In classical games that have been applied to ecology, individual fitness is either density independent or population density is fixed. This article focuses on the habitat selection game where fitness depends on the population density that evolves over time. This model assumes that changes in animal distribution operate on a fast time scale when compared to demographic processes. Of particular interest is whether it is true, as one might expect, that resident phenotypes who use density‐dependent optimal foraging strategies are evolutionarily stable with respect to invasions by mutant strategies. In fact, we show that evolutionary stability does not require that residents use the evolutionarily stable strategy (ESS) at every population density; rather it is the combined resident–mutant system that must be at an evolutionary stable state. That is, the separation of time scales assumption between behavioral and ecological processes does not imply that these processes are independent. When only consumer population dynamics in several habitats are considered (i. e. when resources do not undergo population dynamics), we show that the existence of optimal foragers forces the resident‐mutant system to approach carrying capacity in each habitat even though the mutants do not die out. Thus, the ideal free distribution (IFD) for the single‐species habitat selection game becomes an evolutionarily stable state that describes a mixture of resident and mutant phenotypes rather than a strategy adopted by all individuals in the system. Also discussed is how these results are affected when animal distribution and demographic processes act on the same time scale.     

Evolution of conditional dispersal: evolutionarily stable strategies in spatial models

We consider a two-species competition model in which the species have the same population dynamics but different dispersal strategies. Both species disperse by a combination of random diffusion and advection along environmental gradients, with the same random dispersal rates but different advection coefficients. Regarding these advection coefficients as movement strategies of the species, we investigate their course of evolution. By applying invasion analysis we find that if the spatial environmental variation is less than a critical value, there is a unique evolutionarily singular strategy, which is also evolutionarily stable. If the spatial environmental variation exceeds the critical value, there can be three or more evolutionarily singular strategies, one of which is not evolutionarily stable. Our results suggest that the evolution of conditional dispersal of organisms depends upon the spatial heterogeneity of the environment in a subtle way.     

An explicit approach to evolutionarily stable dispersal strategies: no cost of dispersal

The evolution of dispersal is examined by looking at evolutionarily stable strategies (ESS) for dispersal parameters in discrete time multisite models without any cost of dispersal. ESS are investigated analytically, based on explicit results on sensitivity analysis of matrix models. The basic model considers an arbitrary number of sites and a single age class. An ESS for dispersal parameters is obtained when the spatial reproductive values, calculated at the density-dependent population equilibrium, are equal across sites. From this basic formulation, one derives equivalently that all local populations should be at equilibrium in the absence of migration, and that dispersal between sites should be balanced, i.e., the numbers of individuals arriving to and leaving a site are equal. These results are then generalized to a model with several age classes. Equal age-specific reproductive values do not however imply balanced dispersal in this case. Our results generalize to any number of sites and age classes those available ?M. Doebeli, Dispersal and dynamics, Theoret. Popul. 47 (1995) 82 for two sites and one age class.     

The evolution of dispersal in a two-patch system: some consequences of differences between migrants and residents

We investigate how age-structure and differences in certain demographic traits between residents and immigrants of a single species act to determine the evolutionarily stable dispersal strategy in a two-patch environment that is heterogeneous in space but constant in time. These two factors have been neglected in previous models of the evolution of dispersal, which generally consider organisms with very simple life-cycles and assume that, whatever their origin, individuals in a given habitat have the same bio-demographic characteristics. However, there is increasing empirical evidence that dispersing individuals have different demographic properties from phylopatric ones. We develop a matrix model in which recruitment depends on local population densities. We assume that dispersal entails a proportional cost to immigrant fecundity, which can be compensated by differences in survival rates between immigrants and residents. The evolutionarily stable strategies (ESS) for dispersal are identified using a combination of analytical expressions and numerical simulations. Our results show that philopatry is selected (1) when dispersal rates do not vary in space, (2) when the metapopulation is a source-sink system and (3) when dispersal rates vary in space (asymmetric dispersal) and immigrants do not compensate for their reduced fecundity. We observe that non-zero asymmetric dispersal rates may be evolutionarily stable when (1) immigrants and residents are demographically alike and (2) immigrants compensate totally for their reduced fecundity through an increase in adult survival. Under these conditions, we find that the ESS occurs when the fitnesses at equilibrium in the two habitats, measured in our model by the realized reproductive rates, are each equal to unity. A comparison with previous studies suggests a unifying rule for the evolution of dispersal: the dispersal rates which permit the spatial homogenization of fitnesses are ESSs. This condition provides new insight into the evolutionary stability of source-sink systems. It also supports the hypothesis that immigrants have adapted demographic strategies, rather than the hypothesis that dispersal is costly and immigrants are at a disavantage compared with residents.     

Stabilization through spatial pattern formation in metapopulations with long-range dispersal

Many studies of metapopulation models assume that spatially extended populations occupy a network of identical habitat patches, each coupled to its nearest neighbouring patches by density-independent dispersal. Much previous work has focused on the temporal stability of spatially homogeneous equilibrium states of the metapopulation, and one of the main predictions of such models is that the stability of equilibrium states in the local patches in the absence of migration determines the stability of spatially homogeneous equilibrium states of the whole metapopulation when migration is added. Here, we present classes of examples in which deviations from the usual assumptions lead to different predictions. In particular, heterogeneity in local habitat quality in combination with long-range dispersal can induce a stable equilibrium for the metapopulation dynamics, even when within-patch processes would produce very complex behaviour in each patch in the absence of migration. Thus, when spatially homogeneous equilibria become unstable, the system can often shift to a different, spatially inhomogeneous steady state. This new global equilibrium is characterized by a standing spatial wave of population abundances. Such standing spatial waves can also be observed in metapopulations consisting of identical habitat patches, i.e. without heterogeneity in patch quality, provided that dispersal is density dependent. Spatial pattern formation after destabilization of spatially homogeneous equilibrium states is well known in reaction–diffusion systems and has been observed in various ecological models. However, these models typically require the presence of at least two species, e.g. a predator and a prey. Our results imply that stabilization through spatial pattern formation can also occur in single-species models. However, the opposite effect of destabilization can also occur: if dispersal is short range, and if there is heterogeneity in patch quality, then the metapopulation dynamics can be chaotic despite the patches having stable equilibrium dynamics when isolated. We conclude that more general metapopulation models than those commonly studied are necessary to fully understand how spatial structure can affect spatial and temporal variation in population abundance.     

The Allee-type ideal free distribution

The ideal free distribution (IFD) in a two-patch environment where individual fitness is positively density dependent at low population densities is studied. The IFD is defined as an evolutionarily stable strategy of the habitat selection game. It is shown that for low and high population densities only one IFD exists, but for intermediate population densities there are up to three IFDs. Population and distributional dynamics described by the replicator dynamics are studied. It is shown that distributional stability (i.e., IFD) does not imply local stability of a population equilibrium. Thus distributional stability is not sufficient for population stability. Results of this article demonstrate that the Allee effect can strongly influence not only population dynamics, but also population distribution in space.     

Continuous-discrete models of the dynamics of an isolated population and of 2 competing species

The paper presents the analysis of various mathematical models for dynamics of isolated population and for competition between two species. It is assumed that mortality is continuous and birth of individuals of new generations takes place in certain fixed moments. Influence of winter upon the population dynamics and conditions of classic discrete model "deduction" of population dynamics (in particular, Moran-Ricker and Hassel's models) are investigated. Dynamic regimes of models under various assumptions about the birth and death rates upon the population states are also examined. Analysis of models of isolated population dynamics with nonoverlapping generations showed the density changes regularly if the birth rate is constant. Moreover, there exists a unique global stable level and population size stabilizes asymptotically at this equilibrium, i.e. cycle and chaotic regimes in various discrete models depend on correlation between individual productivity and population state in previous time. When the correlation is exponential upon mean population size the discrete Hassel model is realized. Modification of basis model, based on the assumption that during winter survival/death changes are constant, showed that population size at global level is stable. Generally, the dependence of population rate upon "winter parameters" has nonlinear character. Nonparametric models of competition between two species does not vary if the individual productivity is constant. In a phase space there are several stable stationary states and population stabilizes at one or other level asymptotically. So, in discrete models of competition between two species oscillation can be explained by dependence of population growth rate on the population size at previous times.     

A novel fitness proxy in structured locally finite metapopulations with diploid genetics, with an application to dispersal evolution

Many studies of evolutionarily stable strategies (ESS) for technical reasons make the simplification that reproduction is clonal. A post-hoc justification is that in the simplest eco-evolutionary models more realistic genetic assumptions, such as haploid sexual or diploid sexual cases, yield results compatible with the clonal ones. For metapopulations the technical reasons were even more poignant thanks to the lack of accessible fitness proxies for the diploid case. However, metapopulations are also precisely the sort of ecological backdrop for which one expect discrepancies between the evolutionary outcomes derived from clonal reproduction and diploid genetics, because substantially many mutant homozygotes appear locally even though the mutant is rare globally. In this paper we devise a fitness proxy applicable to the haploid sexual and diploid sexual case, in the style of Metz and Gyllenberg [Metz, J.A.J., Gyllenberg, M., 2001. How should we define fitness in structured metapopulation models? Including an application to the calculation of ES dispersal strategies. Proc. R. Soc. Lond. B 268, 499-508], that can cope with local population fluctuations due to environmental and demographic stochasticity. With the use of this fitness proxy we find that in dispersal evolution the studied clonal model is equivalent with the haploid sexual model, and that there are indeed many differences between clonal and diploid ESS dispersal rates. In a homogenous landscape the discrepancy is but minor (less than 2%), but the situation is different in a heterogeneous landscape: Not only is the quantitative discrepancy between the two types of ESSs appreciable (around 10%-20%), but more importantly, at the same parameter values, evolutionarily stability properties may differ. It is possible, that the singular strategy is evolutionarily stable in the clonal case but not in the diploid case, and vice versa.     

The evolution of patch selection in stochastic environments

A null model for habitat patch selection in spatially heterogeneous environments is the ideal free distribution (IFD), which assumes individuals have complete knowledge about the environment and can freely disperse. Under equilibrium conditions, the IFD predicts that local population growth rates are zero in all occupied patches, sink patches are unoccupied, and the fraction of the population selecting a patch is proportional to the patch's carrying capacity. Individuals, however, often experience stochastic fluctuations in environmental conditions and cannot respond to these fluctuations instantaneously. An evolutionary stability analysis for fixed patch-selection strategies reveals that environmental uncertainty disrupts the classical IFD predictions: individuals playing the evolutionarily stable strategy may occupy sink patches, local growth rates are negative and typically unequal in all patches, and individuals prefer higher-quality patches less than predicted by their carrying capacities. Spatial correlations in environmental fluctuations can enhance or marginalize these trends. The analysis predicts that continually increasing environmental variation first selects for range expansion, then selects for persisting coupled sink populations, and ultimately leads to regional extinction. In contrast, continually increasing habitat degradation first selects for range contraction and may select for persisting coupled sink populations before regional extinction. These results highlight the combined roles of spatial and temporal heterogeneity on the evolution of habitat selection.     

Evolutionary games for spatially dispersed populations

An evolutionary game model is developed that incorporates both spatial dispersion and density effects in the evolutionary dynamic. It is shown that a stable equilibrium (e.g. an evolutionarily stable strategy) of the non-dispersed frequency dynamic becomes a stable equilibrium of the larger system if population density stabilizes at these fixed frequencies. It is also shown, by example, that other equilibria, whose frequencies change from one location to another, may appear when dispersal rates are relatively small.Research supported by Natural Sciences and Engineering Research Council of Canada Operating Grant A6187Research supported by Natural Sciences and Engineering Research Council of Canada Operating Grant A7822     

Classical metapopulation dynamics and eco‐evolutionary feedbacks in dendritic networks

Eco‐evolutionary dynamics are now recognized to be highly relevant for population and community dynamics. However, the impact of evolutionary dynamics on spatial patterns, such as the occurrence of classical metapopulation dynamics, is less well appreciated. Here, we analyse the evolutionary consequences of spatial network connectivity and topology for dispersal strategies and quantify the eco‐evolutionary feedback in terms of altered classical metapopulation dynamics. We find that network properties, such as topology and connectivity, lead to predictable spatio‐temporal correlations in fitness expectations. These spatio‐temporally stable fitness patterns heavily impact evolutionarily stable dispersal strategies and lead to eco‐evolutionary feedbacks on landscape level metrics, such as the number of occupied patches, the number of extinctions and recolonizations as well as metapopulation extinction risk and genetic structure. Our model predicts that classical metapopulation dynamics are more likely to occur in dendritic networks, and especially in riverine systems, compared to other types of landscape configurations. As it remains debated whether classical metapopulation dynamics are likely to occur in nature at all, our work provides an important conceptual advance for understanding the occurrence of classical metapopulation dynamics which has implications for conservation and management of spatially structured populations.