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     

Spatially structured metapopulation models: global and local assessment of metapopulation capacity

We model metapopulation dynamics in finite networks of discrete habitat patches with given areas and spatial locations. We define and analyze two simple and ecologically intuitive measures of the capacity of the habitat patch network to support a viable metapopulation. Metapopulation persistence capacity lambda(M) defines the threshold condition for long-term metapopulation persistence as lambda(M)>delta, where delta is defined by the extinction and colonization rate parameters of the focal species. Metapopulation invasion capacity lambda(I) sets the condition for successful invasion of an empty network from one small local population as lambda(I)>delta. The metapopulation capacities lambda(M) and lambda(I) are defined as the leading eigenvalue or a comparable quantity of an appropriate "landscape" matrix. Based on these definitions, we present a classification of a very general class of deterministic, continuous-time and discrete-time metapopulation models. Two specific models are analyzed in greater detail: a spatially realistic version of the continuous-time Levins model and the discrete-time incidence function model with propagule size-dependent colonization rate and a rescue effect. In both models we assume that the extinction rate increases with decreasing patch area and that the colonization rate increases with patch connectivity. In the spatially realistic Levins model, the two types of metapopulation capacities coincide, whereas the incidence function model possesses a strong Allee effect characterized by lambda(I)=0. For these two models, we show that the metapopulation capacities can be considered as simple sums of contributions from individual habitat patches, given by the elements of the leading eigenvector or comparable quantities. We may therefore assess the significance of particular habitat patches, including new patches that might be added to the network, for the metapopulation capacities of the network as a whole. We derive useful approximations for both the threshold conditions and the equilibrium states in the two models. The metapopulation capacities and the measures of the dynamic significance of particular patches can be calculated for real patch networks for applications in metapopulation ecology, landscape ecology, and conservation biology.     

On fitness in structured metapopulations

In this paper we derive a general expression measuring fitness in general structured metapopulation models. We apply the theory to a model structured by local population size and in which local dynamics is explicitly modelled. In particular, we calculate the evolutionarily stable dispersal strategy for individuals that can assess the local population density in the case where only dispersal is subject to evolutionary control but all other model ingredients are assumed fixed. We show that there exists a threshold size such that at ESS everyone should stay as long as the population size is below the threshold and everyone should disperse immediately as the population size reaches the threshold. Received: 13 August 1999 / Revised version: 2 May 2001 / Published online: 12 October 2001     

Viability analyses with habitat-based metapopulation models

Population viability analysis (PVA) models incorporate spatial dynamics in different ways. At one extreme are the occupancy models that are based on the number of occupied populations. The simplest occupancy models ignore the location of populations. At the other extreme are individual-based models, which describe the spatial structure with the location of each individual in the population, or the location of territories or home ranges. In between these are spatially structured metapopulation models that describe the dynamics of each population with structured demographic models and incorporate spatial dynamics by modeling dispersal and temporal correlation among populations. Both dispersal and correlation between each pair of populations depend on the location of the populations, making these models spatially structured. In this article, I describe a method that expands spatially structured metapopulation models by incorporating information about habitat relationships of the species and the characteristics of the landscape in which the metapopulation exists. This method uses a habitat suitability map to determine the spatial structure of the metapopulation, including the number, size, and location of habitat patches in which subpopulations of the metapopulation live. The habitat suitability map can be calculated in a number of different ways, including statistical analyses (such as logistic regression) that find the relationship between the occurrence (or, density) of the species and independent variables which describe its habitat requirements. The habitat suitability map is then used to calculate the spatial structure of the metapopulation, based on species-specific characteristics such as the home range size, dispersal distance, and minimum habitat suitability for reproduction. Received: April 1, 1999 / Accepted: October 29, 1999     

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.     

Stochastic models for mainland-island metapopulations in static and dynamic landscapes

This paper has three primary aims: to establish an effective means for modelling mainland-island metapopulations inhabiting a dynamic landscape; to investigate the effect of immigration and dynamic changes in habitat on metapopulation patch occupancy dynamics; and to illustrate the implications of our results for decision-making and population management. We first extend the mainland-island metapopulation model of Alonso and McKane [Bull. Math. Biol. 64:913–958, 2002] to incorporate a dynamic landscape. It is shown, for both the static and the dynamic landscape models, that a suitably scaled version of the process converges to a unique deterministic model as the size of the system becomes large. We also establish that, under quite general conditions, the density of occupied patches, and the densities of suitable and occupied patches, for the respective models, have approximate normal distributions. Our results not only provide us with estimates for the means and variances that are valid at all stages in the evolution of the population, but also provide a tool for fitting the models to real metapopulations. We discuss the effect of immigration and habitat dynamics on metapopulations, showing thatmainland-like patches heavily influence metapopulation persistence, and we argue for adopting measures to increase connectivity between this large patch and the other island-like patches. We illustrate our results with specific reference to examples of populations of butterfly and the grasshopperBryodema tuberculata.     

The influence of patch demographics on metapopulations, with particular reference to successional landscapes

The classical view of metapopulations relates the regional abundance of a species to the balance between the extinction and colonization dynamics of identical local populations. Species in successional landscapes may represent the most appropriate examples of classical metapopulations. However, Levins‐type metapopulation models do not explicitly separate population loss due to successional habitat change from other causes of extinction. A further complication is that the chance of population loss due to successional habitat change may be related to the age of a patch. I developed simple patch occupancy models to include succession and included consideration of patch age structure to address two related questions: what are the implications of changes in patch demographic rates and when is a move to a structured patch occupancy model justified? Age‐related variation in patch demography could increase or decrease the equilibrium fraction of the available habitat occupied by a species when compared to the predictions of an unstructured model. Metapopulation persistence was enhanced when the age class of patches with the highest species occupancy suffered relatively low losses to habitat succession. Conversely, when the age class of patches with the highest species occupancy also had relatively high successional loss rates, extinction thresholds were higher that would be predicted by a simple unstructured model. Hence age‐related variation in patch successional rate introduces biases into the predictions of simple unstructured models. Such biases can be detected from field surveys of the fraction of occupied and unoccupied patches in each age class. Where a bias is demonstrated, unstructured models will not be adequate for making predictions about the effects of changing parameters on metapopulation size. Thinking in successional terms emphasizes how landscapes might be managed to enhance or reduce the patch occupancy by any particular metapopulation     

Metapopulation viability analysis of the bog fritillary butterfly using RAMAS/GIS

Many species living in man-shaped landscapes are restricted to small natural habitat patches and form metapopulations; predicting their future is a central issue in applied ecology. We examined the viability of the bog fritillary butterfly Proclossiana eunomia Esper, a specialist glacial relict species, in a highly fragmented landscape (<1% of suitable habitat in 10 km²), by way of population viability analysis. We used comprehensive data from a long-term study in which a patchy population was monitored during ten consecutive years to parameterise a spatially structured metapopulation model using commercially available platform RAMAS/GIS 3.0. Population growth rate was density-dependent and modulated by various climatic variables acting on different developmental stages of the butterfly. Density dependence was probably related to larval parasitism by a specific parasitoid. Population size was negatively affected by an increase in the mean temperature. Dispersal was modelled as the observed proportion of movements between patches, taking into account the probability of emigration out of a given patch. Our model provided results close to the picture of the system drawn from the field data and was considered as useful in making predictions about the metapopulation. Demographic parameters proved to have a far higher impact on metapopulation persistence than dispersal or correlation of local dynamics. Scenarios simulating both global warming and management of habitat patches by rustic herbivore grazing indicated a decrease in the viability of the metapopulation. Our results prompted the regional nature conservation agency to modify the planned management regime. We urge conservation biologists to use structured population models including local population dynamics for viability analysis targeted to such threatened metapopulations in highly fragmented landscapes.     

Habitat Deterioration, Habitat Destruction, and Metapopulation Persistence in a Heterogenous Landscape

Levins's unstructured metapopulation model predicts that the equilibrium fraction of empty habitat patches is a constant function of the fractionhof suitable patches in the landscape and that this constant equals the threshold value for metapopulation persistence. Levins's model thus suggests that the minimum amount of suitable habitat necessary for metapopulation persistence can be estimated from the fraction of empty patches at steady state. In this paper we construct several more realistic structured metapopulation models that include variation in patch quality and the rescue effect. These models predict both positive and negative correlations between the fractions of suitable patches and empty patches. The type of correlation depends in an intricate manner on the strength of the rescue effect and on the quality distribution of the patches to be destroyed. Empty patches can be considered as the resource limiting metapopulation growth. Our results demonstrate that the correlation between the fractions of suitable patches and empty patches is positive if and only if the average value of the resource decreases as the number of patches increases.     

Evolution of specialization in resource utilization in structured metapopulations

We study the evolution of resource utilization in a structured discrete-time metapopulation model with an infinite number of patches, prone to local catastrophes. The consumer faces a trade-off in the abilities to consume two resources available in different amounts in each patch. We analyse how the evolution of specialization in the utilization of the resources is affected by different ecological factors: migration, local growth, local catastrophes, forms of the trade-off and distribution of the resources in the patches. Our modelling approach offers a natural way to include more than two patch types into the models. This has not been usually possible in the previous spatially heterogeneous models focusing on the evolution of specialization.     

Evolution of specialization in resource utilization in structured metapopulations

We study the evolution of resource utilization in a structured discrete-time metapopulation model with an infinite number of patches, prone to local catastrophes. The consumer faces a trade-off in the abilities to consume two resources available in different amounts in each patch. We analyse how the evolution of specialization in the utilization of the resources is affected by different ecological factors: migration, local growth, local catastrophes, forms of the trade-off and distribution of the resources in the patches. Our modelling approach offers a natural way to include more than two patch types into the models. This has not been usually possible in the previous spatially heterogeneous models focusing on the evolution of specialization.     

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.     

Hopf bifurcation in a structured population model for the sexual phase of monogonont rotifers

 We are studying a population of monogonont rotifers in the context of non-linear age-dependent models. In the sexual phase of their reproductive cycle we consider the population structured by age, and composed of three subclasses: virgin mictic females, mated mictic females, and haploid males. The model system has a unique stationary population density which is stable as long as a parameter, related to male-female encounter rate, remains below a critical value. When the parameter increases beyond this critical value, the stationary solution becomes unstable and a stable limit cycle (isolated periodic orbit) appears. The occurrence of this supercritical Hopf bifurcation is shown analytically. Received: 2 August 2001 / Revised version: 3 January 2002 / Published online: 26 June 2002     

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.     

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.     

How does temporal variation in habitat connectivity influence metapopulation dynamics?

Metapopulation persistence depends on connectivity between habitat patches. While emphasis has been placed on the spatial dynamics of connectivity, much less has been placed on its short‐term temporal dynamics. In many terrestrial and aquatic ecosystems, however, transient (short‐term) changes in connectivity occur as habitat patches are connected and disconnected due, for example, to climatic or hydrological variability. We evaluated the implications of transient connectivity using a network‐based metapopulation model and a series of scenarios representing temporal changes in connectivity. The transient loss of connectivity can influence metapopulation persistence, and more strongly autocorrelated temporal dynamics affect metapopulation persistence more severely. Given that many ecosystems experience short‐term and temporary loss of habitat connectivity, it is important that these dynamics are adequately represented in metapopulation models; failing to do so may yield overly optimistic‐estimates of metapopulation persistence in fragmented landscapes.     

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.     

A Stochastic Metapopulation Model with Variability in Patch Size and Position

Analytically tractable metapopulation models usually assume that every patch is identical, which limits their application to real metapopulations. We describe a new single species model of metapopulation dynamics that allows variation in patch size and position. The state of the metapopulation is defined by the presence or absence of the species in each patch. For a system of n patches, this gives 2ⁿ possible states. We show how to construct and analyse a matrix describing transitions between all possible states by first constructing separate extinction and colonisation matrices. We illustrate the model′s application to metapopulations by considering an example of malleefowl, Leipoa ocellata, in southern Australia, and calculate extinction probabilities and quasi-stationary distributions. We investigate the relative importance of modelling the particular arrangement of patches and the variation in patch sizes for this metapopulation and we use the model to examine the effects of further habitat loss on extinction probabilities.     

Metapopulation dynamics of the bog fritillary butterfly: modelling the effect of habitat fragmentation

Population viability analysis (PVA) and metapopulation theory are valuable tools to model the dynamics of spatially structured populations. In this article we used a spatially realistic population dynamic model to simulate the trajectory of a Proclossiana eunomia metapopulation in a network of habitat patches located in the Belgian Ardenne. Sensitivity analysis was used to evaluate the relative influence of the different parameters on the model output. We simulated habitat loss by removing a percentage of the original habitat, proportionally in each habitat patch. Additionally, we evaluated isolation and fragmentation effects by removing and dividing habitat patches from the network, respectively. The model predicted a slow decline of the metapopulation size and occupancy. Extinction risks predicted by the model were highly sensitive to environmental stochasticity and carrying capacity. For a determined level of habitat destruction, the expected lifetime of the metapopulation was highly dependent on the spatial configuration of the landscape. Moreover, when the proportion of removed habitat is above 40% of the original habitat, the loss of whole patches invariably leads to the strongest reduction in metapopulation viability.     

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.