Do invasive species undergo metapopulation dynamics? A case study of the invasive Caspian gull,Larus cachinnans,in Poland

Aim The mechanisms of initial dispersal and habitat occupancy by invasive alien species are fundamental ecological problems. Most tests of metapopulation theory are performed on local population systems that are stable or in decline. In the current study we were interested in the usefulness of metapopulation theory to study patch occupancy, local colonization, extinction and the abundance of the invasive Caspian gull (Larus cachinnans) in its initial invasion stages. Location Waterbodies in Poland. Methods Characteristics of the habitat patches (waterbodies, 35 in total) occupied by breeding pairs of Caspian gulls and an equal sample of randomly selected unoccupied patches were compared with t‐tests. Based on presence–absence data from 1989 to 2006 we analysed factors affecting the probability of local colonization, extinction and the size of local populations using generalized linear models. Results Occupied habitat patches were significantly larger and less isolated (from other habitat patches and other local populations) and were located closer to rivers than empty patches. The proximity of local food resources (fish ponds, refuse dumps) positively affected the occurrence of breeding pairs. The probability of colonization was positively affected by patch area, and negatively by distances to fish ponds, nearest habitat patch, nearest breeding colony and to a river, and by higher forest cover around the patch boundaries. The probability of extinction was lower in patches with a higher number of breeding pairs and with a greater area of islets. The extinction probability increased with distances to other local populations, other habitat patches, fish ponds and to refuse dumps and with a higher cover of forest around the patch boundaries. The size of the local population decreased with distances to the nearest habitat patch, local population, river, fish pond and refuse dump. Local abundance was also positively affected by the area of islets in the patch. Main conclusions During the initial stages of the invasion of Caspian gulls in Poland the species underwent metapopulation‐like dynamics with frequent extinctions from colonized habitat patches. The results prove that metapopulation theory may be a useful conceptual framework for predicting which habitats are more vulnerable to invasion.     

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.     

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     

Single-species metapopulation dynamics: concepts, models and observations

This paper outlines a conceptual and theoretical framework for single-species metapopulation dynamics based on the Levins model and its variants. The significance of the following factors to metapopulation dynamics are explored: evolutionary changes in colonization ability; habitat patch size and isolation; compensatory effects between colonization and extinction rates; the effect of immigration on local dynamics (the rescue effect); and heterogeneity among habitat patches. The rescue effect may lead to alternative stable equilibria in metapopulation dynamics. Heterogeneity among habitat patches may give rise to a bimodal equilibrium distribution of the fraction of patches occupied in an assemblage of species (the core-satellite distribution). A new model of incidence functions is described, which allows one to estimate species' colonization and extinction rates on islands colonized from mainland. Four distinct kinds of stochasticity affecting metapopulation dynamics are discussed with examples. The concluding section describes four possible scenarios of metapopulation extinction.     

Testing metapopulation concepts: effects of patch characteristics and neighborhood occupancy on the dynamics of an endangered lagomorph

Metapopulation ecology is a field that is richer in theory than in empirical results. Many existing empirical studies use an incidence function approach based on spatial patterns and key assumptions about extinction and colonization rates. Here we recast these assumptions as hypotheses to be tested using 18 years of historic detection survey data combined with four years of data from a new monitoring program for the Lower Keys marsh rabbit. We developed a new model to estimate probabilities of local extinction and colonization in the presence of nondetection, while accounting for estimated occupancy levels of neighboring patches. We used model selection to identify important drivers of population turnover and estimate the effective neighborhood size for this system. Several key relationships related to patch size and isolation that are often assumed in metapopulation models were supported: patch size was negatively related to the probability of extinction and positively related to colonization, and estimated occupancy of neighboring patches was positively related to colonization and negatively related to extinction probabilities. This latter relationship suggested the existence of rescue effects. In our study system, we inferred that coastal patches experienced higher probabilities of extinction and colonization than interior patches. Interior patches exhibited higher occupancy probabilities and may serve as refugia, permitting colonization of coastal patches following disturbances such as hurricanes and storm surges. Our modeling approach should be useful for incorporating neighbor occupancy into future metapopulation analyses and in dealing with other historic occupancy surveys that may not include the recommended levels of sampling replication.     

Spatial and temporal variation in patch occupancy and population density in a model system of an arctic Collembola species assemblage

We employed an experimental model system to investigate the mechanisms underlying patterns of patch occupancy and population density in a high arctic assemblage of Collembola species inhabiting a sedge tussock landscape on Svalbard. The replicate model systems consisted of 5 cores of the tussocks (habitat patches) imbedded in a barren matrix. Four of the patches were open so that animals could migrate between them, while there was one closed patch per system to test the effect of migration on extinction rate. There were model systems of two types: one with long and one with short inter‐patch distances to test the effect of patch isolation on colonisation and extinction rates. Each of the four most common collembolan species at the field site were introduced to two open patches per system (source patches), with the other two functioning as colonisation patches for the species. The experiment was run in an ecotrone over three identical, simulated arctic summers separated by winters of 3 weeks. Six replicates of systems with short and long inter‐patch distances were sampled at the end of each summer. The species varied markedly in their performance in both open arenas and closed patches, indicating differential responses to patch humidity, consistent with their differential distribution along the moisture gradient in the field site. The extinction – colonisation dynamics differed markedly between species as predicted from our field studies. This could partly be ascribed to differential dispersal and colonisation ability, but also to different tolerance to spatially variable patch quality and/or tendency for aggregative behaviour. Three of the species exhibited dynamics that superficially resemble what could be expected from classical metapopulation dynamics. However, there was a striking discrepancy between what would be expected from the effect of migration on the extinction rate of isolated patches (in particular closed patches) and the observed rates. Thus, metapopulation processes, such as stochastic colonisation and extinction events due to demographic stochasticity, were relatively unimportant compared to other sources of spatial variability among which subtle differences in patch quality are probably most important. We discuss the value of combining field studies with model system experiments, in particular when habitat quality cannot easily be measured in the field. However, our field and laboratory studies also emphasise the need for a thorough knowledge of species‐specific life history traits for making biologically sound interpretations based on both observational and experimental data.     

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.     

Rules of thumb for conservation of metapopulations based on a stochastic winking-patch model

From a theoretical viewpoint, nature management basically has two options to prolong metapopulation persistence: decreasing local extinction probabilities and increasing colonization probabilities. This article focuses on those options with a stochastic, single-species metapopulation model. We found that for most combinations of local extinction probabilities and colonization probabilities, decreasing the former increases metapopulation extinction time more than does increasing the latter by the same amount. Only for relatively low colonization probabilities is an effort to increase these probabilities more beneficial, but even then, decreasing extinction probabilities does not seem much less effective. Furthermore, we found the following rules of thumb. First, if one focuses on extinction, one should preferably decrease the lowest local extinction probability. Only if the extinction probabilities are (almost) equal should one prioritize decreases in the local extinction probability of the patch with the best direct connections to and from other patches. Second, if one focuses on colonization, one should preferably increase the colonization probability between the patches with the lowest local extinction probability. Only if the local extinction probabilities are (almost) equal should one instead prioritize increases in the highest colonization probability (unless extinction probabilities and colonization probabilities are very low). The rules of thumb have an important common denominator: the local extinction process has a greater bearing on metapopulation extinction time than colonization.     

Metapopulation models for seasonally migratory animals

Metapopulation models are widely used to study species that occupy patchily distributed habitat, but are rarely applied to migratory species, because of the difficulty of identifying demographically independent subpopulations. Here, we extend metapopulation theory to describe the directed seasonal movement of migratory populations between two sets of habitat patches, breeding and non-breeding, with potentially different colonization and extinction rates between patch types. By extending the classic metapopulation model, we show that migratory metapopulations will persist if the product of the two colonization rates exceeds the product of extinction rates. Further, we develop a spatially realistic migratory metapopulation model and derive a landscape metric-the migratory metapopulation capacity-that determines persistence. This new extension to metapopulation theory introduces an important tool for the management and conservation of migratory species and may also be applicable to model the dynamics of two host-parasite systems.     

On optimal size and number of reserves for metapopulation persistence

Habitat fragmentation is generally considered to be detrimental to the persistence of natural populations. In nature management, one therefore tends to prefer few large nature reserves over many small nature reserves having equal total area. This paper examines whether this preference is warranted in a metapopulation framework with circular reserves (patches) by formulating the dependence of metapopulation persistence on the size and number of reserves, both of which depend on reserve radius if the total area is kept constant. Two measures of metapopulation persistence are used: R(0), the number of patches colonized by an occupied patch during its lifetime as an occupied patch, and T(e), the expected time to extinction. These two measures are functions of the extinction and colonization rates of the metapopulation. Several mechanisms for the extinction and colonization processes are formulated from which the dependence of these rates on reserve radius is calculated. It turns out that T(e)generally increases with reserve radius for all mechanisms, which supports the preference of few large reserves. However, R(0)supports this preference only in the case of some special, rather unrealistic, mechanisms. In many other, more realistic, cases an intermediate reserve size exists for which metapopulation persistence measured by R(0)is optimal.     

用空间直观模型是否足以从斑块占据性资料中推断集合种群的动态过程?

在集合种群的研究中,经常要根据空间占据性数据应用斑块模型来推断种群的动态过程,在保护生物学应用中,斑块占据性模型的参数估测对于阐释集合种群动态和预测种群对生境破坏的反应极为重要。我们探讨了一种广泛应用的空间直观模型——率函数模型(Incidence function model)中参数估测的不确定性问题,通过构建由50个斑块组成的网络和两个假想的已知参数的集合种群,应用模拟模型产生集合种群随时间变化的斑块占据性数据系列：即快照(snapshot)。然后,根据这些快照,应用率函数模型和最大似然法估测种群动态参数。此外,我们还给出了传统的率函数模型的一个变形,这个变形包含了目标区效应(Target area effect)：即一个斑块的占据概率不但取决于空间隔离度,也取决于斑块本身面积的大小。结果表明：根据同一个集合种群不同的快照所估测的参数可以有很大差异,一个快照得出的参数提示的是占据性强但存活率低的集合种群,而另一个快照可能反映的是一个占据性弱但存活率高的集合种群。应用传统的率函数模型于一个包含了目标区效应的集合种群,导致斑块大小相关的灭绝率参数估测的正偏差。因此,仅根据一个快照的空间占据性数据来推测集合种群的过程有很大的不确定性[动物学报49(6)：787～794,2003]。     

Spatial patterns of coexistence of competing species in patchy habitat

The single-species spatially realistic patch occupancy metapopulation model is, in this study, extended to a metacommunity of many competing species. Competition is assumed to reduce the local carrying capacity (effective patch area), which in turn increases local extinction rates and reduces colonization rates because of smaller population sizes. Each species is described by three parameters: pre-competitive abundance (equilibrium incidence of patch occupancy, which reflects the rate of colonization in relation to extinction rate), the spatial range of migration, and competitive ability. The model ignores spatio–temporal correlations caused by interspecific interactions, because in metacommunities of unequal competitors inhabiting heterogeneous landscapes, correlations in the occurrence of species are driven more by patch heterogeneity than by competition. The model allows the calculation of multispecies equilibria in patchy habitats without simulations. In general, the number of coexisting species in the metacommunity increases with decreasing strength of competition, increasing rate of colonization, and decreasing range of migration. Habitat heterogeneity in the form of spatial variation in patch areas tends to facilitate coexistence. Poor competitors may coexist with superior competitors in the patch network if the former have higher colonization rates (competition–colonization trade-off). When migration distances are short, competition leads to spatial pattern formation: Species tend to have restricted spatial distributions in the network, but contrary to intuitive expectations, often the distributions of many species are nested. Having more dispersive species enhances both local and global diversity, whereas more local migration decreases local but increases global diversity.     

Long-term dynamics in a metapopulation of the American pika

A 20-yr study of a metapopulation of the American pika revealed a regional decline in occupancy in one part of a large network of habitat patches. We analyze the possible causes of this decline using a spatially realistic metapopulation model, the incidence function model. The pika metapopulation is the best-known mammalian example of a classical metapopulation with significant population turnover, and it satisfies closely the assumptions of the incidence function model, which was parameterized with data on patch occupancy. The model-predicted incidences of patch occupancy are consistent with observed incidences, and the model predicts well the observed turnover rate between four metapopulation censuses. According to model predictions, the part of the metapopulation where the decline has been observed is relatively unstable and prone to large oscillations in patch occupancy, whereas the other part of the metapopulation is predicted to be persistent. These results demonstrate how extinction-colonization dynamics may produce spatially correlated patterns of patch occupancy without any spatially correlated processes in local dynamics or extinction rate. The unstable part of the metapopulation gives an empirical example of multiple quasi equilibria in metapopulation dynamics. Phenomena similar to those observed here may cause fluctuations in species' range limits.     

Modelling the spatial dynamics and persistence of the leaf beetle Gonioctena olivacea in dynamic habitats

In dynamic landscapes natural and anthropogenic disturbance as well as succession are responsible for the emergence and subsequent disappearance of suitable habitat patches. Species inhabiting such landscapes are faced with varying number and spatial configuration of patches. A stochastic, spatially explicit simulation model was developed in order to analyse the persistence of the leaf beetle Gonioctena olivacea in a system of dynamic patches of its host plant Cytisus scoparius . The model was parameterized with data from a three-year field study on the spatial configuration, distribution, and turnover of the host plant patches as well as the patch occupancy, extinction, and colonization rates of the beetle. The simulations showed large fluctuations in the occurrence of the beetle in the patches. High levels of occupancy were related to high aggregation of the patches within the landscape. The velocity of patch turnover was found to have a severe effect on the persistence of the beetle metapopulation. Enhancing the turnover rate by only a few patches, the mean time to extinction decreases rapidly. Moreover, the results revealed that not necessarily an effect of connectivity can be detected in the analysis of occupancy patterns in dynamic landscapes, although the colonization of patches is clearly connectivity-dependent. In general, this modelling study demonstrates the importance of detailed information on patch turnover. The amount and spatial distribution of suitable habitat is a major driver of metapopulation dynamics of species in dynamic landscapes.     

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     

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.     

Dispersal-induced desynchronization: from metapopulations to metacommunities

The conceptualization of fragmented populations in terms of metapopulation theory has become standard over the last three decades. It is well known that increases in between‐patch migration rates cause more synchronous population fluctuations and that this coherence increases the risk of global metapopulation extinction. Because species’ migration rates and the probability of individuals surviving migration events depend on the effective distance between patches, the benefit of improving conservation corridors or the matrix between habitat patches has been questioned. As populations occur in the context of larger communities, moving from a metapopulation to a metacommunity model framework is a natural extension to address the generality of these conclusions. We show how considering a metacommunity can modify the conclusion that decreasing the effective distance between habitat patches (via improving matrix quality or other measures) necessarily increases the degree of metapopulation synchrony. We show that decreases in effective between‐patch distance may deter population synchrony because of the simultaneous effect this change has on the migration patterns of other species. These results indicate that species interactions need to be considered when the effect of conservation measures on population synchrony, and ultimately persistence, is addressed.     

Simple discrete-time metapopulation models of patch occupancy

Simple models in theoretical ecology have a long-standing history of being used to understand how specific processes influence population dynamics as well as providing a foundation for future endeavors. The Levins model is the seminal example of this for continuous-time metapopulation dynamics. However, many natural populations have a distinct separation between processes and data is not collected continuously leading to the need for using a discrete-time model. Our goal is to develop a simple discrete-time metapopulation model of patch occupancy using difference equations. In our formulation, we consider the two fundamental processes of colonization and extinction that will be treated as sequential events and will only consider patch occupancy. To achieve this, we use a composition of two functions where one will reflect the extinction process and the other for the colonization process. Under some mild assumptions, we are able determine the dynamic behavior of the metapopulation. In addition, we provide numerous examples for the functions used to emulate the colonization and extinction processes. Our results illustrate that the dynamics of the model are tied to properties such as convexity and monotonicity of the colonization and extinction functions. In particular, if the model is non-monotone, then complex dynamics can arise such as cyclic and even chaotic behavior. Overall, our approach shows how certain properties of the colonization and extinction functions can influence metapopulation dynamics.     

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.     

Colonization in metapopulations: a review of theory and observations

In metapopulation dynamics turnover of populations in isolated patches may be frequent. Regional survival of a species in such a system with frequent extinctions hinges on its colonization ability. Colonization is more than just dispersal; when a propagule reaches a new patch it faces higher extinction probabilities than does an established population. Extinction models as well as empirical data suggest that a large propagule with a potential for rapid increase in a varying environment, or with a low mortality rate in an environment perceived as constant, has a higher probability of successful colonization. Large variation in population size when it is still small increases the risk of failure. Factors introducing such variation are demographic stochasticity and environmental variation. It is hard to single out demographic traits that ensure good colonizing ability, since colonization can be achieved in many different ways, but generalists and species with self-fertilization seem to be superior.