Evolution of complex dynamics in spatially structured populations

Dynamics of populations depend on demographic parameters which may change during evolution. In simple ecological models given by one-dimensional difference equations, the evolution of demographic parameters generally leads to equilibrium population dynamics. Here we show that this is not true in spatially structured ecological models. Using a multi-patch metapopulation model, we study the evolutionary dynamics of phenotypes that differ both in their response to local crowding, i.e. in their competitive behaviour within a habitat, and in their rate of dispersal between habitats. Our simulation results show that evolution can favour phenotypes that have the intrinsic potential for very complex dynamics provided that the environment is spatially structured and temporally variable. These phenotypes owe their evolutionary persistence to their large dispersal rates. They typically coexist with phenotypes that have low dispersal rates and that exhibit equilibrium dynamics when alone. This coexistence is brought about through the phenomenon of evolutionary branching, during which an initially uniform population splits into the two phenotypic classes.     

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     

Evolution of Dispersal in a Structured Metapopulation Model in Discrete Time

In this article, a structured metapopulation model in discrete time with catastrophes and density-dependent local growth is introduced. The fitness of a rare mutant in an environment set by the resident is defined, and an efficient method to calculate fitness is presented. With this fitness measure evolutionary analysis of this model becomes feasible. This article concentrates on the evolution of dispersal. The effect of catastrophes, dispersal cost, and local dynamics on the evolution of dispersal is investigated. It is proved that without catastrophes, if all population–dynamical attractors are fixed points, there will be selection for no dispersal. A new mechanism for evolutionary branching is also found: Even though local population sizes approach fixed points, catastrophes can cause enough temporal variability, so that evolutionary branching becomes possible.     

Evolution of migration in a metapopulation

In this paper a general deterministic discrete-time metapopulation model with a finite number of habitat patches is analysed within the framework of adaptive dynamics. We study a general model and prove analytically that (i) if the resident populations state is a fixed point, then the resident strategy with no migration is an evolutionarily stable strategy, (ii) a mutant population with no migration can invade any resident population in a fixed point state, (iii) in the uniform migration case the strategy not to migrate is attractive under small mutational steps so that selection favours low migration. Some of these results have been previously observed in simulations, but here they are proved analytically in a general case. If the resident population is in a two-cyclic orbit, then the situation is different. In the uniform migration case the invasion behaviour depends both on the type of the residents attractor and the survival probability during migration. If the survival probability during migration is low, then the system evolves towards low migration. If the survival probability is high enough, then evolutionary branching can happen and the system evolves to a situation with several coexisting types. In the case of out-of-phase attractor, evolutionary branching can happen with significantly lower survival probabilities than in the in-phase attractor case. Most results in the two-cyclic case are obtained by numerical simulations. Also, when migration is not uniform we observe in numerical simulations in the two-cyclic orbit case selection for low migration or evolutionary branching depending on the survival probability during migration.     

Joint evolution of specialization and dispersal in structured metapopulations

We study the joint evolution of dispersal and specialization concerning resource usage in a mechanistically underpinned structured discrete-time metapopulation model. We show that dispersal significantly affects the evolution of specialization and that specialization is a key factor that determines the possibility of evolutionary branching in dispersal propensity. Allowing both dispersal propensity and specialization to evolve as a consequence of natural selection is necessary in order to understand the evolutionary dynamics. The joint evolution of dispersal and specialization forms a natural evolutionary path leading to the coexistence of generalists and specialists. We show that in this process, the number of different patch types and the resource distribution are essential.     

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.     

Area-dependent migration by ringlet butterflies generates a mixture of patchy population and metapopulation attributes

Interpretation of spatially structured population systems is critically dependent on levels of migration between habitat patches. If there is considerable movement, with each individual visiting several patches, there is one ”patchy population”; if there is intermediate movement, with most individuals staying within their natal patch, there is a metapopulation; and if (virtually) no movement occurs, then the populations are separate (Harrison 1991, 1994). These population types actually represent points along a continuum of much to no mobility in relation to patch structure. Therefore, interpretation of the effects of spatial structure on the dynamics of a population system must be accompanied by information on mobility. We use empirical data on movements by ringlet butterflies, Aphantopus hyperantus, to investigate two key issues that need to be resolved in spatially-structured population systems. First, do local habitat patches contain largely independent local populations (the unit of a metapopulation), or merely aggregations of adult butterflies (as in patchy populations)? Second, what are the effects of patch area on migration in and out of the patches, since patch area varies considerably within most real population systems, and because human landscape modification usually results in changes in habitat patch sizes? Mark-release-recapture (MRR) data from two spatially structured study systems showed that 63% and 79% of recaptures remained in the same patch, and thus it seems reasonable to call both systems metapopulations, with some capacity for separate local dynamics to take place in different local patches. Per capita immigration and emigration rates declined with increasing patch area, while the resident fraction increased. Actual numbers of emigrants either stayed the same or increased with area. The effect of patch area on movement of individuals in the system are exactly what we would have expected if A. hyperantus were responding to habitat geometry. Large patches acted as local populations (metapopulation units) and small patches simply as locations with aggregations (units of patchy populations), all within 0.5 km². Perhaps not unusually, our study system appears to contain a mixture of metapopulation and patchy-population attributes.     

An epidemic model in a patchy environment

An epidemic model is proposed to describe the dynamics of disease spread among patches due to population dispersal. We establish a threshold above which the disease is uniformly persistent and below which disease-free equilibrium is locally attractive, and globally attractive when both susceptible and infective individuals in each patch have the same dispersal rate. Two examples are given to illustrate that the population dispersal plays an important role for the disease spread. The first one shows that the population dispersal can intensify the disease spread if the reproduction number for one patch is large, and can reduce the disease spread if the reproduction numbers for all patches are suitable and the population dispersal rate is strong. The second example indicates that a population dispersal results in the spread of the disease in all patches, even though the disease can not spread in each isolated patch.     

Evolutionary branching of dispersal strategies in structured metapopulations

 Dispersal polymorphism and evolutionary branching of dispersal strategies has been found in several metapopulation models. The mechanism behind those findings has been temporal variation caused by cyclic or chaotic local dynamics, or temporally and spatially varying carrying capacities. We present a new mechanism: spatial heterogeneity in the sense of different patch types with sufficient proportions, and temporal variation caused by catastrophes. The model where this occurs is a generalization of the model by Gyllenberg and Metz (2001). Their model is a size-structured metapopulation model with infinitely many identical patches. We present a generalized version of their metapopulation model allowing for different types of patches. In structured population models, defining and computing fitness in polymorphic situations is, in general, difficult. We present an efficient method, which can be applied also to other structured population or metapopulation models. Received: 6 March 2001 / Revised version: 12 February 2002 / Published online: 17 July 2002     

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.     

Metapopulation persistence in dynamic landscapes: the role of dispersal distance

Species associated with transient habitats need efficient dispersal strategies to ensure their regional survival. Using a spatially explicit metapopulation model, we studied the effect of the dispersal range on the persistence of a metapopulation as a function of the local population and landscape dynamics (including habitat patch destruction and subsequent regeneration). Our results show that the impact of the dispersal range depends on both the local population and patch growth. This is due to interactions between dispersal and the dynamics of patches and populations via the number of potential dispersers. In general, long-range dispersal had a positive effect on persistence in a dynamic landscape compared to short-range dispersal. Long-range dispersal increases the number of couplings between the patches and thus the colonisation of regenerated patches. However, long-range dispersal lost its advantage for long-term persistence when the number of potential dispersers was low due to small population growth rates and/or small patch growth rates. Its advantage also disappeared with complex local population dynamics and in a landscape with clumped patch distribution.     

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.     

The evolution of phenotypic polymorphism: randomized strategies versus evolutionary branching

A population is polymorphic when its members fall into two or more categories, referred to as alternative phenotypes. There are many kinds of phenotypic polymorphisms, with specialization in reproduction, feeding, dispersal, or protection from predators. An individual's phenotype might be randomly assigned during development, genetically determined, or set by environmental cues. These three possibilities correspond to a mixed strategy of development, a genetic polymorphism, and a conditional strategy. Using the perspective of adaptive dynamics, I develop a unifying evolutionary theory of systems of determination of alternative phenotypes, focusing on the relative possibilities for random versus genetic determination. The approach is an extension of the analysis of evolutionary branching in adaptive dynamics. It compares the possibility that there will be evolutionary branching, leading to genetic polymorphism, with the possibility that a mixed strategy evolves. The comparison is based on the strength of selection for the different outcomes. An interpretation of the resulting criterion is that genetic polymorphism is favored over random determination of the phenotype if an individual's heritable genotype is an adaptively advantageous cue for development. I argue that it can be helpful to regard genetic polymorphism as a special case of phenotypic plasticity.     

Does colonization asymmetry matter in metapopulations?

Despite the considerable evidence showing that dispersal between habitat patches is often asymmetric, most of the metapopulation models assume symmetric dispersal. In this paper, we develop a Monte Carlo simulation model to quantify the effect of asymmetric dispersal on metapopulation persistence. Our results suggest that metapopulation extinctions are more likely when dispersal is asymmetric. Metapopulation viability in systems with symmetric dispersal mirrors results from a mean field approximation, where the system persists if the expected per patch colonization probability exceeds the expected per patch local extinction rate. For asymmetric cases, the mean field approximation underestimates the number of patches necessary for maintaining population persistence. If we use a model assuming symmetric dispersal when dispersal is actually asymmetric, the estimation of metapopulation persistence is wrong in more than 50% of the cases. Metapopulation viability depends on patch connectivity in symmetric systems, whereas in the asymmetric case the number of patches is more important. These results have important implications for managing spatially structured populations, when asymmetric dispersal may occur. Future metapopulation models should account for asymmetric dispersal, while empirical work is needed to quantify the patterns and the consequences of asymmetric dispersal in natural metapopulations.     

Effect of habitat fragmentation on dispersal in the butterfly Proclossiana eunomia

Comparison of dispersal rates of the bog fritillary butterfly between continuous and fragmented landscapes indicates that between patch dispersal is significantly lower in the fragmented landscape, while population densities are of the same order of magnitude. Analyses of the dynamics of the suitable habitat for the butterfly in the fragmented landscape reveal a severe, non linear increase in spatial isolation of patches over a time period of 30 years (i.e. 30 butterfly generations), but simulations of the butterfly metapopulation dynamics using a structured population model show that the lower dispersal rates in the fragmented landscape are far above the critical threshold leading to metapopulation extinction. These results indicate that changes in individual behaviour leading to the decrease of dispersal rates in the fragmented landscape were rapidly selected for when patch spatial isolation increased. The evidence of such an adaptive answer to habitat fragmentation suggests that dispersal mortality is a key factor for metapopulation persistence in fragmented landscapes. We emphasise that landscape spatial configuration and patch isolation have to be taken into account in the debate about large-scale conservation strategies.     

Dynamics of Adaptation in Spatially Heterogeneous Metapopulations

The selection pressure experienced by organisms often varies across the species range. It is hence crucial to characterise the link between environmental spatial heterogeneity and the adaptive dynamics of species or populations. We address this issue by studying the phenotypic evolution of a spatial metapopulation using an adaptive dynamics approach. The singular strategy is found to be the mean of the optimal phenotypes in each habitat with larger weights for habitats present in large and well connected patches. The presence of spatial clusters of habitats in the metapopulation is found to facilitate specialisation and to increase both the level of adaptation and the evolutionary speed of the population when dispersal is limited. By showing that spatial structures are crucial in determining the specialisation level and the evolutionary speed of a population, our results give insight into the influence of spatial heterogeneity on the niche breadth of species.     

CHAOS AND UNPREDICTABILITY IN EVOLUTION

The possibility of complicated dynamic behavior driven by nonlinear feedbacks in dynamical systems has revolutionized science in the latter part of the last century. Yet despite examples of complicated frequency dynamics, the possibility of long‐term evolutionary chaos is rarely considered. The concept of “survival of the fittest” is central to much evolutionary thinking and embodies a perspective of evolution as a directional optimization process exhibiting simple, predictable dynamics. This perspective is adequate for simple scenarios, when frequency‐independent selection acts on scalar phenotypes. However, in most organisms many phenotypic properties combine in complicated ways to determine ecological interactions, and hence frequency‐dependent selection. Therefore, it is natural to consider models for evolutionary dynamics generated by frequency‐dependent selection acting simultaneously on many different phenotypes. Here we show that complicated, chaotic dynamics of long‐term evolutionary trajectories in phenotype space is very common in a large class of such models when the dimension of phenotype space is large, and when there are selective interactions between the phenotypic components. Our results suggest that the perspective of evolution as a process with simple, predictable dynamics covers only a small fragment of long‐term evolution.     

Colonization and extinction dynamics of an annual plant metapopulation in an urban environment

The metapopulation framework considers that the spatiotemporal distribution of organisms results from a balance between the colonization and extinction of populations in a suitable and discrete habitat network. Recent spatially realistic metapopulation models have allowed patch dynamics to be investigated in natural populations but such models have rarely been applied to plants. Using a simple urban fragmented population system in which favourable habitat can be easily mapped, we studied patch dynamics in the annual plant Crepis sancta (Asteraceae). Using stochastic patch occupancy models (SPOMs) and multi‐year occupancy data we dissected extinction and colonization patterns in our system. Overall, our data were consistent with two distinct metapopulation scenarios. A metapopulation (sensu stricto) dynamic in which colonization occurs over a short distance and extinction is lowered by nearby occupied patches (rescue effect) was found in a set of patches close to the city centre, while a propagule rain model in which colonization occurs from a large external population was most consistent with data from other networks. Overall, the study highlights the importance of external seed sources in urban patch dynamics. Our analysis emphasizes the fact that plant distributions are governed not only by habitat properties but also by the intrinsic properties of colonization and dispersal of species. The metapopulation approach provides a valuable tool for understanding how colonization and extinction shape occupancy patterns in highly fragmented plant populations. Finally, this study points to the potential utility of more complex plant metapopulation models than traditionally used for analysing ecological and evolutionary processes in natural metapopulations.     

Regional and local scale metapopulation dynamics in the interaction between Callosobruchus maculatus and Anisopteromalus calandrae

Habitat structure increases the persistence of many extinction‐prone resource–consumer interactions. Metapopulation theory is one of the leading approaches currently used to explain why local, ephemeral populations persist at a regional scale. Central to the metapopulation concept is the amount of dispersal occurring between patches, too much or too little can result in regional extinction. In this study, the role of dispersal on the metapopulation dynamics of an over‐exploitative host–parasitoid interaction is assessed. In the absence of the parasitoid the highly vagile bruchid, Callosobruchus maculatus, can maintain a similar population size regardless of the permeability of the inter‐patch matrix and exhibits strong negative density‐dependence. After the introduction of the parasitoid the size of the bruchid population decreases with a corresponding increase in the occurrence of empty patches. In this case, limiting the dispersal of both species decouples the interaction to a greater extent and results in larger regional bruchid populations. Given the disparity between the dispersal rates of the two species, it is proposed that the more dispersive host benefits from the reduction in landscape permeability by increasing the opportunity to colonise empty patches and rescue extinction prone populations. Associated with the introduction of the parasitoid is a shift in the strength of density‐dependence as the population moves from bottom–up towards top–down regulation. The importance of local and regional scale measurements is apparent when the role of individual patches on regional dynamics is considered. By only taking regional dynamics into account the importance of dispersal regime on local dynamics is overlooked. Similarly, when local dynamics were examined, patches were found to have different influences on regional dynamics depending on dispersal regime and patch location.     

Metapopulation persistence in fragmented landscapes: significant interactions between genetic and demographic processes

We formulated a mathematical model in order to study the joint influence of demographic and genetic processes on metapopulation viability. Moreover, we explored the influence of habitat structure, matrix quality and disturbance on the interplay of these processes. We showed that the conditions that allow metapopulation persistence under the synergistic action of genetic and demographic processes depart significantly from predictions based on a mere superposition of the effects of each process separately. Moreover, an optimal dispersal rate exists that maximizes the range of survival rates of dispersers under which metapopulation persists and at the same time allows the largest sustainable patch removal and patch‐size reduction. The relative impact of patch removal and patch‐size reduction depends both on matrix quality and the dispersal strategy of the species: metapopulation persistence is more affected by patch‐size reduction (patch removal) for low (high)‐dispersing species, in presence of a low (high) quality matrix. Avoidance of inbreeding, through increased dispersal when the rate of inbreeding in a population is large, has positive effects on low‐dispersing species, but impairs the persistence of high‐dispersing species. Finally, size heterogeneity between patches largely influences metapopulation dynamics; the presence of large patches, even at the expense of other patches being smaller, can have positive effects on persistence in particular for species of low dispersing ability.