Is metapopulation patch occupancy in nature well predicted by the Levins model?

Although the Levins model has made important theoretical contributions to ecology, its empirical support has not been conclusively established yet. We used published colonization and extinction data from 55 metapopulations to calculate their Levins equilibrium patch occupancy. Over all species, there were not significant differences between the observed patch occupancies and the Levins model's estimates. However, invertebrates and vertebrate species with some degree of threat had patch occupancies larger than the model's expectancies. A temporal sampling effect was found for invertebrate species, with departure from the Levins model decreasing as the length of the study period increased. There was a negative relationship between patch occupancy and extinction probability, as expected under the “rescue effect”. The high rates at which invertebrates produce propagules could lead the Levins model to underestimate patch occupancy, whereas the observed patch occupancy of threatened species may be a transient phenomenon that results from extinction probabilities that increase over time. Therefore, the Levins model captures the metapopulation dynamics of a wide range of species in a simple formula whereas its equilibrium point can be used as evidence of metapopulation stability. Although mechanistic models provide more precise and accurate metapopulation predictions, they also can sacrifice the generality and simplicity of the Levins model.     

Evaluating the metapopulation consequences of ecological traps

Ecological traps occur when environmental changes cause maladaptive habitat selection. Despite their relevance to metapopulations, ecological traps have been studied predominantly at local scales. How these local impacts scale up to affect the dynamics of spatially structured metapopulations in heterogeneous landscapes remains unexplored. We propose that assessing the metapopulation consequences of traps depends on a variety of factors that can be grouped into four categories: the probability of encounter, the likelihood of selection, the fitness costs of selection and species-specific vulnerability to these costs. We evaluate six hypotheses using a network-based metapopulation model to explore the relative importance of factors across these categories within a spatial context. Our model suggests (i) traps are most severe when they represent a large proportion of habitats, severely reduce fitness and are highly attractive, and (ii) species with high intrinsic fitness will be most susceptible. We provide the first evidence that (iii) traps may be beneficial for metapopulations in rare instances, and (iv) preferences for natal-like habitats can magnify the effects of traps. Our study provides important insight into the effects of traps at landscape scales, and highlights the need to explicitly consider spatial context to better understand and manage traps within metapopulations.     

Theoretical models of species' borders: single species approaches

The range of potential mechanisms limiting species' distributions in space is nearly as varied and complex as the diversity of life itself. Yet viewed abstractly, a species' border is a geographic manifestation of a species' demographic responses to a spatially and temporally varying world. Population dynamic models provide insight into the different routes by which range limits can arise owing to gradients in demographic rates. In a metapopulation context, for example, range limits may be caused by gradients in extinction rates, colonization rates or habitat availability. We have consider invasion models in uniform and heterogeneous environments as a framework for understanding non-equilibrium range limits, and explore conditions under which invasions may cease to spread leaving behind a stationary range limit. We conclude that non-equilibrial range dynamics need further theoretical and empirical attention.     

Spatial Heterogeneity in Habitat Quality and Cross-Scale Interactions in Metapopulations

Abstract Integration of habitat heterogeneity into spatially realistic metapopulation approaches reveals the potential for key cross-scale interactions. Broad-scale environmental gradients and land-use practices can create autocorrelation of habitat quality of suitable patches at intermediate spatial scales. Patch occupancy then depends not only on habitat quality at the patch scale but also on feedbacks from surrounding neighborhoods of autocorrelated patches. Metapopulation dynamics emerge from how demographic and dispersal processes interact with relevant habitat heterogeneity. We provide an empirical example from a metapopulation of round-tailed muskrats (Neofiber alleni) in which habitat quality of suitable patches was spatially autocorrelated most strongly within 1,000 m, which was within the expected dispersal range of the species. After controlling for factors typically considered in metapopulation studies—patch size, local patch quality, patch connectivity—we use a cross-variogram analysis to demonstrate that patch occupancy by muskrats was correlated with habitat quality across scales ≤1,171 m. We also discuss general consequences of spatial heterogeneity of habitat quality for metapopulations related to potential cross-scale interactions. We focus on spatially correlated extinctions and metapopulation persistence, hierarchical scaling of source–sink dynamics, and dispersal decisions by individuals in relation to information constraints.     

The community context of species' borders: ecological and evolutionary perspectives

Species distributional limits may coincide with hard dispersal barriers or physiological thresholds along environmental gradients, but they may also be influenced by species interactions. We explore a number of models of interspecific interactions that lead to (sometimes abrupt) distribution limits in the presence and absence of environmental gradients. We find that gradients in competitive ability can lead to spatial segregation of competitors into distinct ranges, but that spatial movement tends to broaden the region of sympatry between the two species, and that Allee effects tend to sharpen these boundaries. We generalize these simple models to include metapopulation dynamics and other types of interactions including predator–prey and host–parasite interactions. We derive conditions for range limits in each case. We also consider models that include coevolution and gene flow and find that character displacement along environmental gradients can lead to stable parapatric distributions. We conclude that it is essential to consider coevolved species interactions as a potential mechanism limiting species distributions, particularly when barriers to dispersal are weak and environmental gradients are gradual.     

The effect of conspecific attraction on metapopulation dynamics

Random dispersal direction is assumed in all current metapopulation models. This assumption is called into question by recent experiments demonstrating that some species disperse preferentially to sites occupied by conspecifies. We incorporate conspecific attraction into two metapopulation models which differ in type of dispersal, the Levins model and a two-dimensional stepping-stone model. In both models, conspecific attraction lowers the proportion of occupied habitat patches within a metapopulation at equilibrium.     

How populations persist when asexuality requires sex: the spatial dynamics of coping with sperm parasites

The twofold cost of sex implies that sexual and asexual reproduction do not coexist easily. Asexual forms tend to outcompete sexuals but may eventually suffer higher extinction rates, creating tension between short- and long-term advantages of different reproductive modes. The 'short-sightedness' of asexual reproduction takes a particularly intriguing form in gynogenetic species complexes, in which an asexual species requires sperm from a related sexual host species to trigger embryogenesis. Asexuals are then predicted to outcompete their host, after which neither species can persist. We examine whether spatial structure can explain continued coexistence of the species complex, and assess the evidence based on data on the Amazon molly (Poecilia formosa). A modification of the Levins metapopulation model creates two regions of good prospects for coexistence, connected by a region of poorer patch occupancy levels. In the first case, mate discrimination and/or niche differentiation keep local extinction rates low, and most patches contain both species; the other possibility resembles host-parasite dynamics where parasites frequently drive the host locally extinct. Several dynamical features are counterintuitive and relate to the parasitic nature of interactions in the species complex: for example, high local extinction rates of the asexual species can be beneficial for its own persistence. This creates a link from the evolution of sexual reproduction to that of prudent predation.     

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.     

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.     

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

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

Metapopulation dynamics: brief history and conceptual domain

We review the early development of metapopulation ideas, which culminated in the well-known model by Levins in 1969. We present a survey of metapopulation terminology and outline the kinds of studies that have been conducted on single-species and multispecies metapopulations. Metapopulation studies have important conceptual links with the equilibrium theory of island biogeography and with studies on the dynamics of species living in patchy environments. Metapopulation ideas play an increasingly important role in landscape ecology and conservation biology.     

异质种群动态模型:破碎化景观动态模拟的新途径

景观破碎化导致物种以异质种群方式存活,使得基于异质种群动态模拟破碎化景观动态成为可能。异质种群动态模型的发展为景观动态模拟奠定了良好基础。根据空间处理方式的不同,异质种群模型可分为三大类,可不同程度地用于描述破碎化景观动态。(1)空间不确定异质种群模型,假定所有局域种群间均等互联,模型中不包含空间信息,仅能用于景观斑块动态描述;(2)空间确定异质种群模型,假设局域种群在二维空间上以规则格子形式排列,是一种准现实的空间处理方式,可用于景观动态的简单描述;(3)空间现实异质种群模型,包含了破碎化景观中局域种群的几何特征,可直接用于真实景观动态的模拟研究。空间现实的和基于个体的异质种群模型不但是未来异质种群模型发展的主流,也将成为未来破碎化景观动态研究的重要工具。为了更加准确完整地描述破碎化景观动态,不但应该综合运用已有的各种异质种群模型方法,更要引进新模型来刎画多物种、多变量、高维度、复杂连接的破碎化景观格局与过程。     

Evaluating the potential causes of range limits of birds of the Colombian Andes

Aim To evaluate how factors acting at different spatial scales influence range limits in bird species of the Colombian Andes. Location Andes Mountains of Colombia. Methods We used Maxent , a climate envelope model (CEM), and environmental and geographic information to study range‐filling (i.e. the extent to which a species occurs in all the areas in which it is predicted to occur) in 70 range‐restricted bird species of the Colombian Andes. Environmental data were taken from the WorldClim database, and species occurrence data were taken from museum data collated by the BioMap project, an observational database, and the literature. We evaluated how climate and geographic barriers may shape range limits at two scales. Results At a broad extent (i.e. across the three main cordilleras within the Colombian Andes), we find that CEMs predict there to be suitable environmental conditions for particular species in regions where the species is absent, possibly as a result of dispersal limitation or biotic interactions. In contrast, at a finer scale (within a given cordillera), species generally occur across the entire area predicted to be suitable by a given CEM. Geographic discontinuities within cordilleras do not generally correspond to range limits; instead, range limits correspond to changes in environmental conditions. Main conclusions Our results suggest that different mechanisms influence the presence of species at different scales. Dispersal limitation, potentially combined with species interactions, may influence range limits at a broad extent (the entire Colombian Andes), while strong environmental gradients correspond to range limits at a finer scale (within a cordillera).     

Transient dynamics in metapopulation response to perturbation

Transient time in population dynamics refers to the time it takes for a population to return to population-dynamic equilibrium (or close to it) following a perturbation in the environment or in population size. Depending on the direction of the perturbation, transient time may either denote the time until extinction (or until the population has decreased to a lower equilibrium level), or the recovery time needed to reach a higher equilibrium level. In the metapopulation context, the length of the transient time is set by the interplay between population dynamics and landscape structure. Assuming a spatially realistic metapopulation model, we show that transient time is a product of four factors: the strength of the perturbation, the ratio between the metapopulation capacity of the landscape and a threshold value determined by the properties of the species, and the characteristic turnover rate of the species, adjusted by a factor depending on the structure of the habitat patch network. Transient time is longest following a large perturbation, for a species which is close to the threshold for persistence, for a species with slow turnover, and in a habitat patch network consisting of only a few dynamically important patches. We demonstrate that the essential behaviour of the n-dimensional spatially realistic Levins model is captured by the one-dimensional Levins model with appropriate parameter transformations.     

Local populations of different sizes, mechanistic rescue effect and patch preference in the Levins metapopulation model

In this paper three extensions of the Levins metapopulation model are discussed: (1) It is shown that the Levins model is still valid if patches contain local populations of different sizes with different colonization and extinction rates. (2) A more mechanistic formulation of the rescue effect is presented. (3) The addition of preference of dispersers for occupied or empty patches and its consequences for conservation strategies are studied.     

Colonization and persistence ability explain the extent to which plant species fill their potential range

Aim How species traits and environmental conditions affect biogeographical dynamics is poorly understood. Here we test whether estimates of a species’ evolutionary age, colonization and persistence ability can explain its current ‘range filling’ (the ratio between realized and potential range size). Location Fynbos biome (Cape Floristic Region, South Africa). Methods For 37 species of woody plants (Proteaceae), we estimate range filling using atlas data and distribution models, evolutionary age using molecular phylogenies, and persistence ability using estimates of individual longevity (which determines the probability of extinction of local populations). Colonization ability is estimated from validated process‐based seed dispersal models, the arrangement of potential habitat, and data on local abundance. To relate interspecific variation in range filling to evolutionary age, colonization and persistence ability, we use two complementary model types: phenomenological linear models and the process‐based metapopulation model of Levins. Results Linear model analyses show that range filling increases with a species’ colonization and persistence ability but is not affected by species age. Moreover, colonization ability is a better predictor of range filling than its component variables (local abundance and dispersal ability). The phylogenetically independent interaction between colonization and persistence ability is significant (P < 0.05) for 97% of 180 alternative phylogenies. While the selected linear model explains 42% of the variance in arcsine transformed range filling, the Levins model performs more poorly. It overestimates range filling for realistic parameter values and produces unrealistic parameter estimates when fitted statistically. Main conclusions Colonization and local extinction seem to shape Proteaceae range dynamics on ecological rather than macroevolutionary time‐scales. Our results suggest that the positive abundance–range size relationship in this group is due primarily to the effect of abundance on colonization. In summary, this study contributes to a process‐based understanding of range dynamics and highlights the importance of colonization for the future survival of Fynbos Proteaceae.     

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

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