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).     

The control of emigration and its consequences for the survival of populations

Dispersal is the key process enhancing the long-term persistence of metapopulations in heterogeneous and dynamic landscapes. However, any individual emigrating from a occupied patch also increases the risk of local population extinction. The consequences of this increase for metapopulation persistence likely depend on the control of emigration. In this paper, we present results of individual-based simulations to compare the consequences of density-independent (DIE) and density-dependent (DDE) emigration on the extinction risk of local populations and a two-patch metapopulation. (1) For completely isolated patches extinction risk increases linearly with realised emigration rates in the DIE scenario. (2) For the DDE scenario extinction risk is nearly insensitive to emigration as longs as emigration probabilities remain below ≈0.2. Survival chances are up to half an order of magnitude larger than for populations with DIE. (3) For low dispersal mortality both modes of emigration increase survival of a metapopulation by ca. one order of magnitude. (4) For high dispersal mortality only DDE can improve the global survival chances of the metapopulation. (5) With DDE individuals are only removed from a population at high population density and the risk of extinction due to demographic stochasticity is thus much smaller compared to the DIE scenario.With density-dependent emigration prospects of metapopulations survival may thus be much higher compared to a system with density-independent emigration. Consequently, the knowledge about the factors driving emigration may significantly affect our conclusions concerning the conservation status of species.     

Effects of Climate Change on Population Persistence of Desert-Dwelling Mountain Sheep in California

Abstract:  Metapopulations may be very sensitive to global climate change, particularly if temperature and precipitation change rapidly. We present an analysis of the role of climate and other factors in determining metapopulation structure based on presence and absence data. We compared existing and historical population distributions of desert bighorn sheep ( Ovis canadensis ) to determine whether regional climate patterns were correlated with local extinction. To examine all mountain ranges known to hold or to have held desert bighorn populations in California and score for variables describing climate, metapopulation dynamics, human impacts, and other environmental factors, we used a geographic information system (GIS) and paper maps. We used logistic regression and hierarchical partitioning to assess the relationship among these variables and the current status of each population (extinct or extant). Parameters related to climate—elevation, precipitation, and presence of dependable springs—were strongly correlated with population persistence in the twentieth century. Populations inhabiting lower, drier mountain ranges were more likely to go extinct. The presence of domestic sheep grazing allotments was negatively correlated with population persistence. We used conditional extinction probabilities generated by the logistic-regression model to rank native, naturally recolonized, and reintroduced populations by vulnerability to extinction under several climate-change scenarios. Thus risk of extinction in metapopulations can be evaluated for global-climate-change scenarios even when few demographic data are available.     

The inflationary effects of environmental fluctuations ensure the persistence of sink metapopulations

Under current rates of environmental change many populations may be found in habitats of low quality and low conservation value, creating population sinks. We test recent theory that suggests, surprisingly, that stochastic environmental variability may enhance the long-term persistence of sink metapopulations. Using experimental populations of Paramecium aurelia we show that it is possible for a metapopulation comprised entirely of sink populations to persist for many generations in a random environment. In accordance with the theory, we show that positive temporal autocorrelation and low spatial correlation in the environment can ensure the long-term persistence and enhance the mean and maximum abundance of sink metapopulations. High levels of spatial correlation in the environment created strong population synchrony and limited the persistence time of the sink metapopulations. These results have important implications for the development of a theory underlying the synergistic effects of habitat fragmentation and environmental change on population persistence.     

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.     

Populations in small, ephemeral habitat patches may drive dynamics in a Daphnia magna metapopulation

Migration is the key process to understand the dynamics and persistence of a metapopulation. Many metapopulation models assume a positive correlation between habitat patch size or stability and the number of emigrants. However, few empirical data exist, and habitat patch size and habitat stability may affect dispersal differently than they affect local persistence. Here, we studied the production of the migration stage (i.e., resting eggs called ephippia) of the cladoceran Daphnia magna in a metapopulation consisting of 530 rock pool habitat patches over 25 years. Earlier, the functioning of this metapopulation was explained with a Levins-type metapopulation model or with a mainland-island metapopulation model, based on local extinction and colonization data or time series data, respectively. We used pool volume, hydroperiod length, and number of desiccation events to calculate per-pool production of ephippia (i.e., migration stages). We estimated that populations in small and ephemeral habitat patches produced more than half of the 250 000 to 1 million ephippia that were produced in the metapopulation as a whole per year between 1982 and 2006. Furthermore, these small populations contributed approximately 90% of the ephippia exposed during desiccation events, while the contribution of the long-lived populations in large pools was minimal. We term this an "inverse mainland-island" type metapopulation and propose that populations in small, ephemeral habitat patches may also be the driving force for metapopulation dynamics in other systems.     

Influence of Landscape Elements in Riparian Buffers on the Conservation of Semiaquatic Amphibians

Abstract:  Studies on riparian buffers have usually focused on the amount of land needed as habitat for the terrestrial life stages of semiaquatic species. Nevertheless, the landscape surrounding wetlands is also important for other key processes, such as dispersal and the dynamics of metapopulations. Multiple elements that influence these processes should therefore be considered in the delineation of buffers. We analyzed landscape elements (forest cover, density of roads, and hydrographic network) in concentric buffers to evaluate the scale at which they influence stream amphibians in 77 distinct landscapes. To evaluate whether our results could be generalized to other contexts, we determined whether they were consistent across the study areas. Amphibians required buffers of 100–400 m of suitable terrestrial habitat, but interspecific differences in the amount of habitat were large. The presence of amphibians was related to roads and the hydrographic network at larger spatial scales (300–1500 m), which suggests that wider buffers are needed with these elements. This pattern probably arose because these elements influence dispersal and metapopulation persistence, processes that occur at large spatial scales. Furthermore, in some cases, analyses performed on different sets of landscapes provided different results, which suggests caution should be used when conservation recommendations are applied to disparate areas. Establishment of riparian buffers should not be focused only on riparian habitat, but should take a landscape perspective because semiaquatic species use multiple elements for different functions. This approach can be complex because different landscape elements require different spatial extents. Nevertheless, a shift of attention toward the management of different elements at multiple spatial scales is necessary for the long-term persistence of populations.     

Dispersal connectivity and reserve selection for marine conservation

Although larval dispersal is crucial for the persistence of most marine populations, dispersal connectivity between sites is rarely considered in designing marine protected area networks. In particular the role of structural characteristics (known as topology) for the network of larval dispersal routes in the conservation of metapopulations has not been addressed. To determine reserve site configurations that provide highest persistence values with respect to their connectivity characteristics, we model nine connectivity topological models derived from graph theory in a demographic metapopulation model. We identify reserve site configurations that provide the highest persistence values for each of the metapopulation connectivity models. Except for the minimally connected and fully connected populations, we observed two general ‘rules of thumb’ for optimising the mean life time for all topological models: firstly place the majority of reserves, so that they are neighbours of each other, on the sites where the number of connections between the populations is highest (hub), secondly when the reserves have occupied the majority of the vertices in the hub, then select another area of high connectivity and repeat. If there are no suitable hubs remaining then distribute the remaining reserves to isolated locations optimising contact with non-reserved sites.     

Colonization and extinction in dynamic habitats: an occupancy approach for a Great Plains stream fish assemblage

Despite the importance of habitat in determining species distribution and persistence, habitat dynamics are rarely modeled in studies of metapopulations. We used an integrated habitat-occupancy model to simultaneously quantify habitat change, site fidelity, and local colonization and extinction rates for larvae of a suite of Great Plains stream fishes in the Arikaree River, eastern Colorado, USA, across three years. Sites were located along a gradient of flow intermittency and groundwater connectivity. Hydrology varied across years: the first and third being relatively wet and the second dry. Despite hydrologic variation, our results indicated that site suitability was random from one year to the next. Occupancy probabilities were also independent of previous habitat and occupancy state for most species, indicating little site fidelity. Climate and groundwater connectivity were important drivers of local extinction and colonization, but the importance of groundwater differed between periods. Across species, site extinction probabilities were highest during the transition from wet to dry conditions (range: 0.52-0.98), and the effect of groundwater was apparent with higher extinction probabilities for sites not fed by groundwater. Colonization probabilities during this period were relatively low for both previously dry sites (range: 0.02-0.38) and previously wet sites (range: 0.02-0.43). In contrast, no sites dried or remained dry during the transition from dry to wet conditions, yielding lower but still substantial extinction probabilities (range: 0.16-0.63) and higher colonization probabilities (range: 0.06-0.86), with little difference among sites with and without groundwater. This approach of jointly modeling both habitat change and species occupancy will likely be useful to incorporate effects of dynamic habitat on metapopulation processes and to better inform appropriate conservation actions.     

Simple rules for ranking and optimally managing metapopulations

We present two ‘rules of thumb’ for metapopulation management. The first identifies an explicit formula for the persistence time of the population, and thus enables the population manager to form a priority species ranking by identifying those species most at risk of extinction. The second identifies an optimal management strategy that gives direction on how to alter the colonisation rate (creation or improvement of habitat corridors) and local extinction rate (restoring habitat quality or expanding habitat) in order to maximise the persistence time under a budgetary constraint. We employ a simple stochastic version of Levins (1969) metapopulation model, which is first calibrated to a more realistic spatial model. Our rules are tested on computer-generated patch networks and a model for malleefowl (Leipoa ocellata) in the Bakara region of South Australia.     

Linking life history to landscape for threatened species conservation in a multiuse region

Landscape-scale conservation that considers metapopulation dynamics will be essential for preventing declines of species facing multiple threats to their survival. Toward this end, we developed a novel approach that combines occurrence records, spatial–environmental data, and genetic information to model habitat, connectivity, and patterns of genetic structure and link spatial attributes to underlying ecological mechanisms. Using the threatened northern quoll (Dasyurus hallucatus) as a case study, we applied this approach to address the need for conservation decision-making tools that promote resilient metapopulations of this threatened species in the Pilbara, Western Australia, a multiuse landscape that is a hotspot for biodiversity and mining. Habitat and connectivity were predicted by different landscape characteristics. Whereas habitat suitability was overwhelmingly driven by terrain ruggedness, dispersal was facilitated by proximity to watercourses. Although there is limited evidence for major physical barriers in the Pilbara, areas with high silt and clay content (i.e., alluvial and hardpan plains) showed high resistance to dispersal. Climate subtlety shaped distributions and patterns of genetic turnover, suggesting the potential for local adaptation. By understanding these spatial–environmental associations and linking them to life-history and metapopulation dynamics, we highlight opportunities to provide targeted species management. To support this, we have created habitat, connectivity, and genetic uniqueness maps for conservation decision-making in the region. These tools have the potential to provide a more holistic approach to conservation in multiuse landscapes globally.     

Extinction, Colonization, and Metapopulations: Environmental Tracking by Rare Species

Extinction and metapopulation theories emphasize that stochastic fluctuations in local populations cause extinction and that local extinctions generate empty habitat patches that are then available for recolonization. Metapopulation persistence depends on the balance of extinction and colonization in a static environment. For many rare and declining species, I argue (1) that extinction is usually the deterministic consequence of the local environment becoming unsuitable (through habitat loss or modification, introduction of a predator, etc.); (2) that the local environment usually remains unsuitable following local extinction, so extinctions only rarely generate empty patches of suitable habitat; and (3) that colonization usually follows improvement of the local environment for a particular species (or long-distance transfer by humans). Thus, persistence depends predominantly on whether organisms are able to track the shifting spatial mosaic of suitable environmental conditions or on maintainance of good conditions locally.     

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.     

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.     

Continuous versus binary representations of landscape heterogeneity in spatially-explicit models of mobile populations

How a landscape is represented is an important structural assumption in spatially-explicit simulation models. Simple models tend to specify just habitat and non-habitat (binary), while more complex models may use multiple levels or a continuum of habitat quality (continuous). How these different representations influence model projections is unclear. To assess the influence of landscape representation on population models, I developed a general, individual-based model with local dispersal and examined population persistence across binary and continuous landscapes varying in the amount and fragmentation of habitat. In binary and continuous landscapes habitat and non-habitat were assigned a unique mean suitability. In continuous landscapes, suitability of each individual site was then drawn from a normal distribution with fixed variance. Populations went extinct less often and abundances were higher in continuous landscapes. Production in habitat and non-habitat was higher in continuous landscapes, because the range of habitat suitability sampled by randomly dispersing individuals was higher than the overall mean habitat suitability. Increasing mortality, dispersal distance, and spatial heterogeneity all increased the discrepancy between continuous and binary landscapes. The effect of spatial structure on the probability of extinction was greater in binary landscapes. These results show that, under certain circumstances, model projections are influenced by how variation in suitability within a landscape is represented. Care should be taken to assess how a given species actually perceives the landscape when conducting population viability analyses or empirical validation of theory.     

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.     

Conservation Corridors and Contagious Disease: A Cautionary Note

Recent conservation proposals frequently include the establishment of corridors to connect isolated patches of wildlife habitat. Much attention has been focused on the potential benefits of corridors with little note given to potentially adverse consequences. A simulation model is developed here to study the effect of corridors on the survival of a metapopulation in the presence of a fatal disease that is spread by direct contact between susceptible and infected individuals. For the disease modeled here, a landscape of patches connected by corridors generally suffers fewer metapopulation extinctions than a landscape of isolated patches. However, under a narrow range of conditions, results suggest that corridors may dramatically increase the probability of metapopulation extinction. This occurs when disease-induced mortality is low enough to allow infected individuals to spread the disease, but high enough to reduce population levels to the point that random demographic and environmental events cause frequent metapopulation extinctions. This has important implications for the design and management of conservation reserve networks. Although discussion focuses primarily on conservation corridors, the model results apply to any management techniques that increase the movement of individuals among populations.     

The Application of Neutral Landscape Models in Conservation Biology

Neutral landscape models, derived from percolation theory in the field of landscape ecology, are grid-based maps in which complex habitat distributions are generated by random or fractal algorithms. This grid-based representation of landscape structure is compatible with the raster-based format of geographical information systems (GIS), which facilitates comparisons between theoretical and real landscapes. Neutral landscape models permit the identification of critical thresholds in connectivity, which can be used to predict when landscapes will become fragmented. The coupling of neutral landscape models with generalized population models, such as metapopulation theory, provides a null model for generating predictions about population dynamics in fragmented landscapes. Neutral landscape models can contribute to the following applications in conservation: (1) incorporation of complex spatial patterns in (meta)population models; (2) identification of species' perceptions of landscape structure; (3) determination of landscape connectivity; (4) evaluation of the consequences of habitat fragmentation for population subdivision; (5) identification of the domain of metapopulation dynamics; (6) prediction of the occurrence of extinction thresholds; ( 7) determination of the genetic consequences of habitat fragmentation; and (8) reserve design and ecosystem management. This generalized, spatially explicit framework bridges the gap between spatially implicit, patch-based models and spatially realistic GIS applications which are usually parameterized for a single species in a specific landscape. Development of a generalized, spatially explicit framework is essential in conservation biology because we will not be able to develop individual models for every species of management concern.     

Population viability of the Siberian Tiger in a changing landscape: Going, going and gone?

The Amur tiger (Panthera tigris altaica) is a flagship species of the boreal forest ecosystem in northeastern China and Russia Far East. During the past century, the tiger population has declined sharply from more than 3000 to fewer than 600 individuals, and its habitat has become much smaller and greatly fragmented. Poaching, habitat degradation, habitat loss, and habitat fragmentation have been widely recognized as the primary causes for the observed population decline. Using a population viability analysis tool (RAMAS/GIS), we simulated the effects of poaching, habitat degradation, habitat loss, and habitat fragmentation on the population dynamics and extinction risk of the Amur tiger, and then explored the relative effectiveness of three conservation strategies involving improving habitat quality and establishing movement corridors in China and Russia. A series of controlled simulation experiments were performed based on the current spatial distribution of habitat and field-observed vital rates. Our results showed that the Amur tiger population could be viable for the next 100 years if the current habitat area and quality were well-maintained, with poaching strictly prohibited of the tigers and their main prey species. Poaching and habitat degradation (mainly prey scarcity) had the largest negative impacts on the tiger population persistence. While the effect of habitat loss was also substantial, habitat fragmentation per se had less influence on the long-term fate of the tiger population. However, to sustain the subpopulations in both Russia and China would take much greater conservation efforts. The viability of the Chinese population of tigers would rely heavily on its connectivity with the largest patch on the other side of the border. Improving the habitat quality of small patches only or increasing habitat connectivity through movement corridors alone would not be enough to guarantee the long-term population persistence of the Amur tiger in both Russia and China. The only conservation strategy that allowed for long-term persistence of tigers in both countries required both the improvement of habitat quality and the establishment of a transnational reserve network. Our study provides new insights into the metapopulation dynamics and persistence of the Amur tiger, which should be useful in landscape and conservation planning for protecting the biggest cat species in the world.