The role of habitat disturbance and recovery in metapopulation persistence

Classical metapopulation theory assumes a static landscape. However, empirical evidence indicates many metapopulations are driven by habitat succession and disturbance. We develop a stochastic metapopulation model, incorporating habitat disturbance and recovery, coupled with patch colonization and extinction, to investigate the effect of habitat dynamics on persistence. We discover that habitat dynamics play a fundamental role in metapopulation dynamics. The mean number of suitable habitat patches is not adequate for characterizing the dynamics of the metapopulation. For a fixed mean number of suitable patches, we discover that the details of how disturbance affects patches and how patches recover influences metapopulation dynamics in a fundamental way. Moreover, metapopulation persistence is dependent not only on the average lifetime of a patch, but also on the variance in patch lifetime and the synchrony in patch dynamics that results from disturbance. Finally, there is an interaction between the habitat and metapopulation dynamics, for instance declining metapopulations react differently to habitat dynamics than expanding metapopulations. We close, emphasizing the importance of using performance measures appropriate to stochastic systems when evaluating their behavior, such as the probability distribution of the state of the metapopulation, conditional on it being extant (i.e., the quasistationary distribution).     

Patch Occupancy and Potential Metapopulation Dynamics of Three Forest Mammals in Fragmented Afromontane Forest in South Africa

Abstract: We investigated the persistence of three medium-sized (2–9 kg), rare forest mammals in the fragmented mist-belt Podocarpus forests of the midlands of KwaZulu-Natal Province, South Africa. We recorded patch occupancy of blue duiker (   Philantomba monticola ), tree hyrax (   Dendrohyrax arboreus ), and samango monkey ( Cercopithecus mitis labiatus ) in 199 forest patches. Their rarity is ascribed to the fragmentation and destruction of their forest habitat. Incidence functions, derived from presence and absence data, were formulated as generalized linear models, and environmental effects were included in the fitted logistic models. The small and mostly solitary hyrax and duiker persisted in smaller patches than the large and social monkey. Although this result follows expectations based on relative home-range sizes of each species, the incidence probability of the samango monkey was invariant with increasing isolation, whereas a gradual decrease with increasing isolation was observed for the hyrax and duiker. Group dynamics may inhibit dispersal and increase the isolation effect in social species such as samango monkeys. A mainland-island metapopulation model adequately describes patterns of patch occupancy by the hyrax and duiker, but the monkeys' poor dispersal ability and obvious area-dependent extirpation suggest that they exist in transient, nonequilibrium (declining) metapopulations. Through identification of large forest patches for careful protection and management, the survival of all three species—especially the monkey—could be prolonged. Because no functional metapopulation may exist for the monkey, however, this is an emergency measure. For the duiker and hyrax, larger patches should form part of a network of smaller and closer patches in a natural matrix.     

A model for behavioral regulation of metapopulation dynamics

Habitat fragmentation can decrease local population persistence by reducing connectivity, which is a function of dispersal of individuals among habitat fragments. Dispersal is often treated as diffusion in population models, even though for many species it is a result of a series of behavioral decisions. We developed a metapopulation model to explore the potential importance of dispersal behaviors in driving metapopulation dynamics. We incorporated types of behavior that affect dispersal—colonization inhibiting, colonization enhancing, extinction inhibiting, extinction enhancing, rescue enhancing, rescue inhibiting—into Levins’ (1969) metapopulation model and projected occupancy rates for a variety of parameter values. Examples from the literature of behaviors associated with each of these parameters are provided. Our model simplifies into previously published metapopulation models that incorporate only a single behavior, and we present a density-dependent rescue function that leads to multiple non-zero equilibria. We found a variety of behavioral effects on metapopulations. Rescue enhancement fills patches faster than does colonization enhancement or extinction inhibition, and declines in patch occupancy are moderate with extinction enhancement, but colonization inhibition causes metapopulation extinction. We also found that with colonization and extinction inhibitions, equilibrium patch occupancy is inversely related to patch turnover rate. With density-dependent rescue, persistence depends not only on the strength of the strong rescue effect, but also on having a sufficient initial fraction of patches occupied; the stronger the rescue effect, the lower this fraction can be. This study suggests that dispersal behavior can have strong influences on metapopulation dynamics. It confirms the importance of understanding the relationship between landscape structure and dispersal behavior in understanding population persistence.     

Metapopulation persistence in a dynamic landscape: more habitat or better stewardship?

Habitat loss and fragmentation has created metapopulations where there were once continuous populations. Ecologists and conservation biologists have become interested in the optimal way to manage and conserve such metapopulations. Several authors have considered the effect of patch disturbance and recovery on metapopulation persistence, but almost all such studies assume that every patch is equally susceptible to disturbance. We investigated the influence of protecting patches from disturbance on metapopulation persistence, and used a stochastic metapopulation model to answer the question: How can we optimally trade off returns from protection of patches vs. creation of patches? We considered the problem of finding, under budgetary constraints, the optimal combination of increasing the number of patches in the metapopulation network vs. increasing the number of protected patches in the network. We discovered that the optimal trade-off is dependent upon all of the properties of the system: the species dynamics, the dynamics of the landscape, and the relative costs of each action. A stochastic model and accompanying methodology are provided allowing a manager to determine the optimal policy for small metapopulations. We also provide two approximations, including a rule of thumb, for determining the optimal policy for larger metapopulations. The method is illustrated with an example inspired by information for the greater bilby, Macrotis lagotis, inhabiting southwestern Queensland, Australia. We found that given realistic costs for each action, protection of patches should be prioritized over patch creation for improving the persistence of the greater bilby during the next 20 years.     

Local Population Dynamics in Metapopulation Models: Implications for Conservation

Abstract: Habitat fragmentation and the division of populations into spatially separated units have led to the increasing use of metapopulation models to characterize these populations. One prominent model that has served as a heuristic tool was introduced by Levins and is based on a collection of simplifying assumptions that exclude information on the dynamics and spatial distribution of local populations. Levins's and similar models predict the proportion of occupied habitat patches at equilibrium and the conditions needed to avoid total extinction. There are many obvious concerns about using such models, including how realistic alterations might change the predictions and whether occupancy has any relationship to population-level processes. Although many of the assumptions of these simple models are known to be unrealistic, we do not know how the assumptions affect model predictions. We simulated a metapopulation, and our results show that assumptions such as homogeneity of habitat patches, random migration among patches, equivalent extinction probabilities in all patches, and a large number of patches can lead to large overestimations of habitat occupancy. But when we explicitly modeled the underlying population dynamics within each patch, we found (1) that there was a strong correlation between proportion of occupied patches and total metapopulation size and (2) that the distribution of individuals among patches was relatively insensitive to model assumptions. Thus, our results show that although realistic modifications will change model predictions for occupancy, occupancy and population trends will be correlated. These correlations between occupancy and population size suggest that occupancy models may have some utility in conservation applications.     

Inferring Metapopulation Dynamics from Patch-Level Incidence of Florida Scrub Plants

Spatial structure and dynamics of multiple populations may explain species distribution patterns in patchy communities with heterogeneous disturbance regimes, especially when species have poor dispersal. The endemic-rich Florida (U.S.A.) rosemary scrub occupies about 4% of the west portion of Archbold Biological Station and occurs scattered within a matrix of less xeric vegetation. Longer fire-return times and higher frequency of open patches in rosemary scrub provide favorable habitat for many plant species. Occupancy of 123 species of vascular plants and ground lichens in 89 patches was determined by repeated site surveys. About two-thirds of the species occurring at more than 14 patches had a significant logistic regression of presence on time-since-fire, patch size, patch isolation, or their interactions. Species with presence related to the interaction between patch isolation and patch size were primarily herbs and small shrubs specializing in rosemary scrub. These results suggest the importance of spatial characteristics of the landscape for population turnover of these species. An incidence-based metapopulation model was used to predict extinction and colonization probabilities of those species with presence in rosemary scrub patches related to the studied spatial variables. This is the first attempt to apply incidence-based metapopulation models to plants. The results showed stronger effects of patch size and patch isolation on extinction probabilities of herbs than on those of woody species. Because of their effect on spatial heterogeneity and habitat availability, fire suppression and habitat destruction may decrease persistence probabilities for these rosemary scrub specialists, many of which are endangered species.     

Quantifying the importance of patch-specific changes in habitat to metapopulation viability of an endangered songbird

A growing number of programs seek to facilitate species conservation using incentive-based mechanisms. Recently, a market-based incentive program for the federally endangered Golden-cheeked Warbler (Dendroica chrysoparia) was implemented on a trial basis at Fort Hood, an Army training post in Texas, USA. Under this program, recovery credits accumulated by Fort Hood through contracts with private landowners are used to offset unintentional loss of breeding habitat of Golden-cheeked Warblers within the installation. Critical to successful implementation of such programs is the ability to value, in terms of changes to overall species viability, both habitat loss and habitat restoration or protection. In this study, we sought to answer two fundamental questions: Given the same amount of change in breeding habitat, does the change in some patches have a greater effect on metapopulation persistence than others? And if so, can characteristics of a patch (e.g., size or spatial location) be used to predict how the metapopulation will respond to these changes? To answer these questions, we describe an approach for using sensitivity analysis of a metapopulation projection model to predict how changes to specific habitat patches would affect species viability. We used a stochastic, discrete-time projection model based on stage-specific estimates of survival and fecundity, as well as various assumptions about dispersal among populations. To assess a particular patch's leverage, we quantified how much metapopulation viability was expected to change in response to changing the size of that patch. We then related original patch size and distance from the largest patch to each patch's leverage to determine if general patch characteristics could be used to develop guidelines for valuing changes to patches within a metapopulation. We found that both the characteristic that best predicted patch leverage and the magnitude of the relationship changed under different model scenarios. Thus, we were unable to find a consistent set of relationships, and therefore we emphasize the dangers in relying on general guidelines to assess patch value. Instead, we provide an approach that can be used to quantitatively evaluate patch value and identify critical needs for future research.     

Comparison of Two Types of Metapopulation Models in Real and Artificial Landscapes

Abstract: Application of metapopulation models is becoming increasingly widespread in the conservation of species in fragmented landscapes. We provide one of the first detailed comparisons of two of the most common modeling techniques, incidence function models and stage-based matrix models, and test their accuracy in predicting patch occupancy for a real metapopulation. We measured patch occupancies and demographic rates for regional populations of the Florida scrub lizard (   Sceloporus woodi ) and compared the observed occupancies with those predicted by each model. Both modeling strategies predicted patch occupancies with good accuracy ( 77–80%) and gave similar results when we compared hypothetical management scenarios involving removal of key habitat patches and degradation of habitat quality. To compare the two modeling approaches over a broader set of conditions, we simulated metapopulation dynamics for 150 artificial landscapes composed of equal-sized patches (2–1024 ha) spaced at equal distances (50–750 m). Differences in predicted patch occupancy were small to moderate (<20%) for about 74% of all simulations, but 22% of the landscapes had differences openface> 50%. Incidence function models and stage-based matrix models differ in their approaches, assumptions, and requirements for empirical data, and our findings provide evidence that the two models can produce different results. We encourage researchers to use both techniques and further examine potential differences in model output. The feasibility of obtaining data for population modeling varies widely among species and limits the modeling approaches appropriate for each species. Understanding different modeling approaches will become increasingly important as conservation programs undertake the challenge of managing for multiple species in a landscape context.     

Consistent scaling of persistence time in metapopulations

Recent theory and experimental work in metapopulations and metacommunities demonstrates that long-term persistence is maximized when the rate at which individuals disperse among patches within the system is intermediate; if too low, local extinctions are more frequent than recolonizations, increasing the chance of regional-scale extinctions, and if too high, dynamics exhibit region-wide synchrony, and local extinctions occur in near unison across the region. Although common, little is known about how the size and topology of the metapopulation (metacommunity) affect this bell-shaped relationship between dispersal rate and regional persistence time. Using a suite of mathematical models, we examined the effects of dispersal, patch number, and topology on the regional persistence time when local populations are subject to demographic stochasticity. We found that the form of the relationship between regional persistence time and the number of patches is consistent across all models studied; however, the form of the relationship is distinctly different among low, intermediate, and high dispersal rates. Under low and intermediate dispersal rates, regional persistence times increase logarithmically and exponentially (respectively) with increasing numbers of patches, whereas under high dispersal, the form of the relationship depends on local dynamics. Furthermore, we demonstrate that the forms of these relationships, which give rise to the bell-shaped relationship between dispersal rate and persistence time, are a product of recolonization and the region-wide synchronization (or lack thereof) of population dynamics. Identifying such metapopulation attributes that impact extinction risk is of utmost importance for managing and conserving the earth's evermore fragmented populations.     

Habitat-quality effects on metapopulation dynamics in greater white-toothed shrews, Crocidura russula

The effects of patch size and isolation on metapopulation dynamics have received wide empirical support and theoretical formalization. By contrast, the effects of patch quality seem largely underinvestigated, partly due to technical difficulties in properly assessing quality. Here we combine habitat-quality modeling with four years of demographic monitoring in a metapopulation of greater white-toothed shrews (Crocidura russula) to investigate the role of patch quality on metapopulation processes. Together, local patch quality and connectivity significantly enhanced local population sizes and occupancy rates (R2 = 14% and 19%, respectively). Accounting for the quality of patches connected to the focal one and acting as potential sources improved slightly the model explanatory power for local population sizes, pointing to significant source-sink dynamics. Local habitat quality, in interaction with connectivity, also increased colonization rate (R2 = 28%), suggesting the ability of immigrants to target high-quality patches. Overall, patterns were best explained when assuming a mean dispersal distance of 800 m, a realistic value for the species under study. Our results thus provide evidence that patch quality, in interaction with connectivity, may affect major demographic processes.     

Modelling dispersal with diffusion and habitat selection: Analytical results for highly fragmented landscapes

Quantifying dispersal is crucial both for understanding ecological population dynamics, and for gaining insight into factors that affect the genetic structure of populations. The role of dispersal becomes pronounced in highly fragmented landscapes inhabited by spatially structured populations. We consider a landscape consisting of a set of habitat patches surrounded by unsuitable matrix, and model dispersal by assuming that the individuals follow a random walk with parameters that may be specific to the habitat type. We allow for spatial variation in patch quality, and account for edge-mediated behavior, the latter meaning that the individuals bias their movement towards the patches when close to an edge between a patch and the matrix. We employ a diffusion approximation of the random walk model to derive analytical expressions for various characteristics of the dispersal process. For example, we derive formulae for the time that an individual is expected to spend in its current patch i, and for the time that it will spend in the matrix, both conditional on the individual hitting next a given patch j before hitting any of the other patches or dying. The analytical formulae are based on the assumptions that the landscape is infinitely large, that the patches are circularly shaped, and that the patches are small compared to interpatch distances. We evaluate the effect of these assumptions by comparing the analytical results to numerical results in a real patch network that violates all of the three assumptions. We then consider a landscape that fulfills the assumptions, and show that in this case the analytical results are in a very good agreement with the numerical results. The results obtained here allow the construction of computationally efficient dispersal models that can be used as components of metapopulation models.     

Role of Patch Size, Disease, and Movement in Rapid Extinction of Bighorn Sheep

Abstract: The controversy (  Berger 1990, 1999 ; Wehausen 1999 ) over rapid extinction in bighorn sheep ( Ovis canadensis ) has focused on population size alone as a correlate to persistence time. We report on the persistence and population performance of 24 translocated populations of bighorn sheep. Persistence in these sheep was strongly correlated with larger patch sizes, greater distance to domestic sheep, higher population growth rates, and migratory movements, as well as to larger population sizes. Persistence was also positively correlated with larger average home-range size ( p = 0.058, n = 10 translocated populations) and home-range size of rams ( p = 0.087, n = 8 translocated populations). Greater home-range size and dispersal rates of bighorn sheep were positively correlated to larger patches. We conclude that patch size and thus habitat carrying capacity, not population size per se, is the primary correlate to both population performance and persistence. Because habitat carrying capacity defines the upper limit to population size, clearly the amount of suitable habitat in a patch is ultimately linked to population size. Larger populations (250+ animals) were more likely to recover rapidly to their pre-epizootic survey number following an epizootic ( p = 0.019), although the proportion of the population dying in the epizootic also influenced the probability of recovery ( p = 0.001). Expensive management efforts to restore or increase bighorn sheep populations should focus on large habitat patches located ≥23 km from domestic sheep, and less effort should be expended on populations in isolated, small patches of habitat.     

The Quantitative Incidence Function Model and Persistence of an Endangered Butterfly Metapopulation

The incidence function model is derived from a linear first-order Markov chain of the presence or absence of a species in a habitat patch. The model can be parameterized with "snapshot" presence/absence data from a patch network. Using the estimated parameter values the Markov chain can be iterated in the same or in some other patch network to generate quantitative predictions about transient metapopulation dynamics and the stochastic steady state. We tested the ability of the incidence function model to predict patch occupancy using extensive data on an endangered butterfly, the Glanville fritillary ( Melitaea cinxia ) Parameter values were estimated with data collected from a 50-patch network in 1991. In 1993 we surveyed the entire geographic range of the species in Finland, within an area of 50 × 70 km², with 1502 habitat patches (dry meadows) of which 536 were occupied. Model predictions were generated for the 1502 patches and were compared with the observed pattern of occupancy in 1993. The model predicted patch occupancy well in more than half of the study area, but prediction was poor for one quarter of the area, probably because of regional variation in habitat quality and because metapopulations may have been perturbed away from the steady state. The incidence function model provides a practical tool for making quantitative predictions about metapopulation dynamics of species living in fragmented landscapes.     

Geometrical Factors and Metapopulation Dynamics of the Bush Cricket, Metrioptera bicolor Philippi (Orthoptera: Tettigoniidae)

This paper presents a metapopulation study of the bush cricket, Metrioptera bicolor , living in a recently fragmented landscape. The species inhabits grass and heathland patches of varying area and isolation. Analyses are made of how these geometrical factors affect local population size and density, distribution pattern, and the probability of local extinction and colonization. The proportion of available patches occupied varied between 72 and 79% during 1985–1990. Unoccupied patches were smaller and more isolated than those that were occupied. Patches where populations became extinct during this period were smaller than those with persisting populations. Since local population size was well correlated with patch area, it was clear that stochastic extinctions only occurred in small populations. Critical patch size for population extinction was approximately half a hectare. Colonized patches were less isolated than those that had not been colonized. Critical inter-patch distance for colonization was about 100 meters. The turnover was restricted to an identifiable share of the available patches. Only 33% of the patches were so small that extinction due to stochastic causes could be considered highly probable. This metapopulation will therefore most likely persist over a considerable period in its present spatial structure. There are apparent threats of further fragmentation, however, and nothing is known about the likelihood of large-scale extinctions resulting from extremely unfavorable weather conditions. Nevertheless, our results show that it is appropriate to include geometrical factors in metapopulation models.     

Dispersal depression with habitat fragmentation in the bog fritillary butterfly

Habitat fragmentation is expected to impose strong selective pressures on dispersal rates. However, evolutionary responses of dispersal are not self-evident, since various selection pressures act in opposite directions. Here we disentangled the components of dispersal behavior in a metapopulation context using the Virtual Migration model, and we linked their variation to habitat fragmentation in the specialist butterfly Proclossiana eunomia. Our study provided a nearly unique opportunity to study how habitat fragmentation modifies dispersal at the landscape scale, as opposed to microlandscapes or simulation studies. Indeed, we studied the same species in four landscapes with various habitat fragmentation levels, in which large amounts of field data were collected and analyzed using similar methodologies. We showed the existence of quantitative variations in dispersal behavior correlated with increased fragmentation. Dispersal propensity from habitat patches (for a given patch size), and mortality during dispersal (for a given patch connectivity) were lower in more fragmented landscapes. We suggest that these were the consequences of two different evolutionary responses of dispersal behavior at the individual level: (1) when fragmentation increased, the reluctance of individuals to cross habitat patch boundaries also increased; (2) when individuals dispersed, they flew straighter in the matrix, which is the best strategy to improve dispersal success. Such evolutionary responses could generate complex nonlinear patterns of dispersal changes at the metapopulation level according to habitat fragmentation. Due to the small size and increased isolation of habitat patches in fragmented landscapes, overall emigration rate and mortality during dispersal remained high. As a consequence, successful dispersal at the metapopulation scale remained limited. Therefore, to what extent the selection of individuals with a lower dispersal propensity and a higher survival during dispersal is able to limit detrimental effects of habitat fragmentation on dispersal success is unknown, and any conclusion that metapopulations would compensate for them is flawed.     

Epiphyte metapopulation dynamics are explained by species traits, connectivity, and patch dynamics

The colonization-extinction dynamics of many species are affected by the dynamics of their patches. For increasing our understanding of the metapopulation dynamics of sessile species confined to dynamic patches, we fitted a Bayesian incidence function model extended for dynamic landscapes to snapshot data on five epiphytic lichens among 2083 mapped oaks (dynamic patches). We estimate the age at which trees become suitable patches for different species, which defines their niche breadth (number of suitable trees). We show that the colonization rates were generally low, but increased with increasing connectivity in accordance with metapopulation theory. The rates were related to species traits, and we show, for the first time, that they are higher for species with wide niches and small dispersal propagules than for species with narrow niches or large propagules. We also show frequent long-distance dispersal in epiphytes by quantifying the relative importance of local dispersal and background deposition of dispersal propagules. Local stochastic extinctions from intact trees were negligible in all study species, and thus, the extinction rate is set by the rate of patch destruction (tree fall). These findings mean that epiphyte metapopulations may have slow colonization-extinction dynamics that are explained by connectivity, species traits, and patch dynamics.     

A Resistant-Kernel Model of Connectivity for Amphibians that Breed in Vernal Pools

Abstract:  Pool-breeding amphibian populations operate at multiple scales, from the individual pool to surrounding upland habitat to clusters of pools. When metapopulation dynamics play a role in long-term viability, conservation efforts limited to the protection of individual pools or even pools with associated upland habitat may be ineffective over the long term if connectivity among pools is not maintained. Connectivity becomes especially important and difficult to assess in regions where suburban sprawl is rapidly increasing land development, road density, and traffic rates. We developed a model of connectivity among vernal pools for the four ambystomatid salamanders that occur in Massachusetts and applied it to the nearly 30,000 potential ephemeral wetlands across the state. The model was based on a modification of the kernel estimator (a density estimator commonly used in home range studies) that takes landscape resistance into account. The model was parameterized with empirical migration distances for spotted salamanders ( Ambystoma maculatum ), dispersal distances for marbled salamanders ( A. opacum ), and expert-derived estimates of landscape resistance. The model ranked vernal pools in Massachusetts by local, neighborhood, and regional connectivity and by an integrated measure of connectivity, both statewide and within ecoregions. The most functionally connected pool complexes occurred in southeastern and northeastern Massachusetts, areas with rapidly increasing suburban development. In a sensitivity analysis estimates of pool connectivity were relatively insensitive to uncertainty in parameter estimates, especially at the local and neighborhood scales. Our connectivity model could be used to prioritize conservation efforts for vernal-pool amphibian populations at broader scales than traditional pool-based approaches.     

Spatial Structure and Population Extinction: A Study with Drosophila Flies

Abstract: The total amount of habitat and also its distribution and subdivision affect the extinction probability of a resident population Two species of Drosophila are studied in spatial configurations of a single large habitat patch, single small habitat patches, and two small but connected habitat patches in which a low rate of migration, roughly one fly per generation, is possible. The single large habitat patch shows the lowest extinction rate lower than the combined rate of two small patches of the same total size. For one of the species, the "corridor" between the pair of small patches seems to produce a "rescue effect" that lowers extinction rates, probably due to a decrease in the coefficient of variation in fluctuations of the population sire in this coupled system. The systems seem to have been influenced by demographic stochasticity, based on the relationship of population size to extinction probability.     

The effect of migration on the viability,dynamics and structure of two coexisting metapopulations

Two species of butterflies, Euphydryas aurinia and Melitaea phoebe, coexist as two metapopulations in a 38-patch network in Hebei Province, China. A Markovian model, whose transition matrix is the product of two matrices which represent the local extinction and recolonization process respectively, is used to describe the metapopulation dynamics. The application of this model to the metapopulation, consisting of 12 local populations in the northern subregion, shows that the expected life times of E. aurinia and M. phoebe are 160 and 121 years respectively and usually nearly half of the patches are occupied by E. aurinia, while only 1–3 patches are occupied by M. phoebe. We claim that E. aurinia can persist for a long time while M. phoebe faces relatively big extinction risk. By comparing the population dynamics with and without migration, we find M. phoebe benefits much more from migration than E. aurinia. Most patches are occupied mainly by local populations for E. aurinia, while by immigrants from the 8th patch for M. phoebe, meaning that E. aurinia has a classical metapopulation structure while M. phoebe has a source–sink metapopulation structure.