Metapopulation persistence in dynamic landscapes: the role of dispersal distance

Species associated with transient habitats need efficient dispersal strategies to ensure their regional survival. Using a spatially explicit metapopulation model, we studied the effect of the dispersal range on the persistence of a metapopulation as a function of the local population and landscape dynamics (including habitat patch destruction and subsequent regeneration). Our results show that the impact of the dispersal range depends on both the local population and patch growth. This is due to interactions between dispersal and the dynamics of patches and populations via the number of potential dispersers. In general, long-range dispersal had a positive effect on persistence in a dynamic landscape compared to short-range dispersal. Long-range dispersal increases the number of couplings between the patches and thus the colonisation of regenerated patches. However, long-range dispersal lost its advantage for long-term persistence when the number of potential dispersers was low due to small population growth rates and/or small patch growth rates. Its advantage also disappeared with complex local population dynamics and in a landscape with clumped patch distribution.     

Metapopulation models for extinction threshold in spatially correlated landscapes

Simple analytical models assuming homogeneous space have been used to examine the effects of habitat loss and fragmentation on metapopulation size. The models predict an extinction threshold, a critical amount of suitable habitat below which the metapopulation goes deterministically extinct. The consequences of non-random loss of habitat for species with localized dispersal have been studied mainly numerically. In this paper, we present two analytical approaches to the study of habitat loss and its metapopulation dynamic consequences incorporating spatial correlation in both metapopulation dynamics as well as in the pattern of habitat destruction. One approach is based on a measure called metapopulation capacity, given by the dominant eigenvalue of a "landscape" matrix, which encapsulates the effects of landscape structure on population extinctions and colonizations. The other approach is based on pair approximation. These models allow us to examine analytically the effects of spatial structure in habitat loss on the equilibrium metapopulation size and the threshold condition for persistence. In contrast to the pair approximation based approaches, the metapopulation capacity based approach allows us to consider species with long as well as short dispersal range and landscapes with spatial correlation at different scales. The two methods make dissimilar assumptions, but the broad conclusions concerning the consequences of spatial correlation in the landscape structure are the same. Our results show that increasing correlation in the spatial arrangement of the remaining habitat increases patch occupancy, that this increase is more evident for species with short-range than long-range dispersal, and that to be most beneficial for metapopulation size, the range of spatial correlation in landscape structure should be at least a few times greater than the dispersal range of the species.     

The role of landscape-dependent disturbance and dispersal in metapopulation persistence

The fundamental processes that influence metapopulation dynamics (extinction and recolonization) will often depend on landscape structure. Disturbances that increase patch extinction rates will frequently be landscape dependent such that they are spatially aggregated and have an increased likelihood of occurring in some areas. Similarly, landscape structure can influence organism movement, producing asymmetric dispersal between patches. Using a stochastic, spatially explicit model, we examine how landscape-dependent correlations between dispersal and disturbance rates influence metapopulation dynamics. Habitat patches that are situated in areas where the likelihood of disturbance is low will experience lower extinction rates and will function as partial refuges. We discovered that the presence of partial refuges increases metapopulation viability and that the value of partial refuges was contingent on whether dispersal was also landscape dependent. Somewhat counterintuitively, metapopulation viability was reduced when individuals had a preponderance to disperse away from refuges and was highest when there was biased dispersal toward refuges. Our work demonstrates that landscape structure needs to be incorporated into metapopulation models when there is either empirical data or ecological rationale for extinction and/or dispersal rates being landscape dependent.     

Asymptotically exact analysis of stochastic metapopulation dynamics with explicit spatial structure

We describe a mathematically exact method for the analysis of spatially structured Markov processes. The method is based on a systematic perturbation expansion around the deterministic, non-spatial mean-field theory, using the theory of distributions to account for space and the underlying stochastic differential equations to account for stochasticity. As an example, we consider a spatial version of the Levins metapopulation model, in which the habitat patches are distributed in the d-dimensional landscape Rd in a random (but possibly correlated) manner. Assuming that the dispersal kernel is characterized by a length scale L, we examine how the behavior of the metapopulation deviates from the mean-field model for a finite but large L. For example, we show that the equilibrium fraction of occupied patches is given by p(0)+c/L(d)+O(L(-3d/2)), where p(0) is the equilibrium state of the Levins model and the constant c depends on p(0), the dispersal kernel, and the structure of the landscape. We show that patch occupancy can be increased or decreased by spatial structure, but is always decreased by stochasticity. Comparison with simulations show that the analytical results are not only asymptotically exact (as L-->infinity), but a good approximation also when L is relatively small.     

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.     

Genetic influences on social dominance: cow wars

Analyzing population dynamics in changing habitats is a prerequisite for population dynamics forecasting. The recent development of metapopulation modeling allows the estimation of dispersal kernels based on the colonization pattern but the accuracy of these estimates compared with direct estimates of the seed dispersal kernel has rarely been assessed. In this study, we used recent genetic methods based on parentage analysis (spatially explicit mating models) to estimate seed and pollen dispersal kernels as well as seed and pollen immigration in fragmented urban populations of the plant species Crepis sancta with contrasting patch dynamics. Using two independent networks, we documented substantial seed immigration and a highly restricted dispersal kernel. Moreover, immigration heterogeneity among networks was consistent with previously reported metapopulation dynamics, showing that colonization was mainly due to external colonization in the first network (propagule rain) and local colonization in the second network. We concluded that the differences in urban patch dynamics are mainly due to seed immigration heterogeneity, highlighting the importance of external population source in the spatio-temporal dynamics of plants in a fragmented landscape. The results show that indirect and direct methods were qualitatively consistent, providing a proper interpretation of indirect estimates. This study provides attempts to link genetic and demographic methods and show that patch occupancy models may provide simple methods for analyzing population dynamics in heterogeneous landscapes in the context of global change.     

Effects of aerial dispersal, habitat specialisation, and landscape structure on spider distribution across fragmented grey dunes

Species distribution patterns have been explained by Hutchinson's niche theory, metapopulation theory and source-sink theory. Empirical verification of this framework, however, remains surprisingly scant. In this paper, we test the hypothesis that landscape characteristics (patch size and connectivity), aerial dispersal ability and niche breadth interact in explaining distribution patterns of 29 spider species inhabiting fragmented grey dunes. Distribution patterns only depended on aerial dispersal potential, and the interaction between patch connectivity and area. Niche breadth, measured as the degree of habitat specialisation in the total coastal dune system, did not contribute to the observed distribution patterns. Additional variation in patch occupancy frequency was strongly species-dependent and was determined by different responses to the degree of patch connectivity for ballooning dispersal. Results from this study suggest that dispersal ability largely affects our perception of a species "fundamental niche", and that source-sink and metapopulation dynamics may have a major impact on the distribution of species. From a conservation point of view, specialised (and hence intrinsically rare) species can be predicted to become rarer if fragmentation increases and connectivity decreases. This study is, to our knowledge, one of the few linking species distribution (and not patch occupancy, species diversity or richness) to landscape ecological (patch connectivity and area) and auto-ecological (niche breadth, dispersal potential) features.     

Rapid assessment of metapopulation viability under climate and land-use change

In order to predict species response to climate and land-use change, numerically fast and easily applicable assessment tools for species survival are required. We present a set of formulae to calculate the mean lifetime of a metapopulation in a spatially heterogeneous and dynamic landscape subject to habitat patch diminution, loss and/or spatial shift of the habitat network. The formulae require as inputs (i) information about the number, location and size of the habitat patches for several time steps to quantify landscape dynamics in terms of patch destruction, diminution or shifting rates and (ii) data on species traits such as their vulnerability to environmental variation and their dispersal ability to quantify local colonisation and extinction rates. We validate the formulae with a spatially explicit simulation. The analysis is complemented by a protocol for the easy use of the approach and practical application examples. A software implementation is available on request from the authors.     

Spatial occurrence and colonisations in patch-tracking metapopulations: local conditions versus dispersal

We studied the relative importance of local variables and dispersal for the occurrence and colonisation of the epiphytic bryophytes Orthotricum speciosum (spore dispersed), and O. obtusifolium (spore and asexual gemmae) on aspen trees ('patches') in two forest landscapes (one old-growth and one fragmented) using multiple logistic regression. The relative importance of dispersal was quantified as the reduction of residual deviance for a connectivity variable. In modelling dispersal, we assumed that trees with low local abundance were recent colonisations, and that trees with high local abundance were diaspore sources for colonisation. The occurrence of O. speciosum in the fragmented landscape was most affected by shading, but also by connectivity, aspen diameter and vitality. In the old-growth landscape, connectivity was the single most important variable for recent colonisations, but its effect was lower than the sum of the effects of all local environmental variables. The occurrence of O. obtusifolium in the fragmented landscape was related to similar variables but the relative importance of these variables was different, and connectivity did not affect the probability of a recent colonisation in this species. We describe the epiphyte-tree system in the patch-tracking metapopulation model. In this model colonisations are distance dependent, but in contrast to the classical metapopulation model local extinctions are caused by deterministic patch destruction – once the epiphyte has colonised the tree it remains until the tree dies.     

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.     

Heterogeneity in patch quality buffers metapopulations from pathogen impacts

Many wildlife species persist on a network of ephemerally occupied habitat patches connected by dispersal. Provisioning of food and other resources for conservation management or recreation is frequently used to improve local habitat quality and attract wildlife. Resource improvement can also facilitate local pathogen transmission, but the landscape-level consequences of provisioning for pathogen spread and habitat occupancy are poorly understood. Here, we develop a simple metapopulation model to investigate how heterogeneity in patch quality resulting from resource improvement influences long-term metapopulation occupancy in the presence of a virulent pathogen. We derive expressions for equilibrium host–pathogen outcomes in terms of provisioning effects on individual patches (through decreased patch extinction rates) and at the landscape level (the fraction of high-quality, provisioned patches), and highlight two cases of practical concern. First, if occupancy in the unprovisioned metapopulation is sufficiently low, a local maximum in occupancy occurs for mixtures of high- and low-quality patches, such that further increasing the number of high-quality patches both lowers occupancy and allows pathogen invasion. Second, if the pathogen persists in the unprovisioned metapopulation, further provisioning can result in all patches becoming infected and in a global minimum in occupancy. This work highlights the need for more empirical research on landscape-level impacts of local resource provisioning on pathogen dynamics.     

Matrix Matters:?Differences of Grand Skink Metapopulation Parameters in Native Tussock Grasslands and Exotic Pasture Grasslands

Modelling metapopulation dynamics is a potentially very powerful tool for conservation biologists. In recent years, scientists have broadened the range of variables incorporated into metapopulation modelling from using almost exclusively habitat patch size and isolation, to the inclusion of attributes of the matrix and habitat patch quality. We investigated the influence of habitat patch and matrix characteristics on the metapopulation parameters of a highly endangered lizard species, the New Zealand endemic grand skink (Oligosoma grande) taking into account incomplete detectability. The predictive ability of the developed zxmetapopulation model was assessed through cross-validation of the data and with an independent data-set. Grand skinks occur on scattered rock-outcrops surrounded by indigenous tussock (bunch) and pasture grasslands therefore implying a metapopulation structure. We found that the type of matrix surrounding the habitat patch was equally as important as the size of habitat patch for estimating occupancy, colonisation and extinction probabilities. Additionally, the type of matrix was more important than the physical distance between habitat patches for colonisation probabilities. Detection probability differed between habitat patches in the two matrix types and between habitat patches with different attributes such as habitat patch composition and abundance of vegetation on the outcrop. The developed metapopulation models can now be used for management decisions on area protection, monitoring, and the selection of translocation sites for the grand skink. Our study showed that it is important to incorporate not only habitat patch size and distance between habitat patches, but also those matrix type and habitat patch attributes which are vital in the ecology of the target species.     

Effect of habitat fragmentation on dispersal in the butterfly Proclossiana eunomia

Comparison of dispersal rates of the bog fritillary butterfly between continuous and fragmented landscapes indicates that between patch dispersal is significantly lower in the fragmented landscape, while population densities are of the same order of magnitude. Analyses of the dynamics of the suitable habitat for the butterfly in the fragmented landscape reveal a severe, non linear increase in spatial isolation of patches over a time period of 30 years (i.e. 30 butterfly generations), but simulations of the butterfly metapopulation dynamics using a structured population model show that the lower dispersal rates in the fragmented landscape are far above the critical threshold leading to metapopulation extinction. These results indicate that changes in individual behaviour leading to the decrease of dispersal rates in the fragmented landscape were rapidly selected for when patch spatial isolation increased. The evidence of such an adaptive answer to habitat fragmentation suggests that dispersal mortality is a key factor for metapopulation persistence in fragmented landscapes. We emphasise that landscape spatial configuration and patch isolation have to be taken into account in the debate about large-scale conservation strategies.     

Anthropogenic landscape change promotes asymmetric dispersal and limits regional patch occupancy in a spatially structured bird population

1. Local extinctions in habitat patches and asymmetric dispersal between patches are key processes structuring animal populations in heterogeneous environments. Effective landscape conservation requires an understanding of how habitat loss and fragmentation influence demographic processes within populations and movement between populations. 2. We used patch occupancy surveys and molecular data for a rainforest bird, the logrunner (Orthonyx temminckii), to determine (i) the effects of landscape change and patch structure on local extinction; (ii) the asymmetry of emigration and immigration rates; (iii) the relative influence of local and between-population landscapes on asymmetric emigration and immigration; and (iv) the relative contributions of habitat loss and habitat fragmentation to asymmetric emigration and immigration. 3. Whether or not a patch was occupied by logrunners was primarily determined by the isolation of that patch. After controlling for patch isolation, patch occupancy declined in landscapes experiencing high levels of rainforest loss over the last 100 years. Habitat loss and fragmentation over the last century was more important than the current pattern of patch isolation alone, which suggested that immigration from neighbouring patches was unable to prevent local extinction in highly modified landscapes. 4. We discovered that dispersal between logrunner populations is highly asymmetric. Emigration rates were 39% lower when local landscapes were fragmented, but emigration was not limited by the structure of the between-population landscapes. In contrast, immigration was 37% greater when local landscapes were fragmented and was lower when the between-population landscapes were fragmented. Rainforest fragmentation influenced asymmetric dispersal to a greater extent than did rainforest loss, and a 60% reduction in mean patch area was capable of switching a population from being a net exporter to a net importer of dispersing logrunners. 5. The synergistic effects of landscape change on species occurrence and asymmetric dispersal have important implications for conservation. Conservation measures that maintain large patch sizes in the landscape may promote asymmetric dispersal from intact to fragmented landscapes and allow rainforest bird populations to persist in fragmented and degraded landscapes. These sink populations could form the kernel of source populations given sufficient habitat restoration. However, the success of this rescue effect will depend on the quality of the between-population landscapes.     

Modelling epiphyte metapopulation dynamics in a dynamic forest landscape

We combine simulations with spatial statistics to estimate the parameters of a metapopulation model for the epiphytic lichen Lobaria pulmonaria specializing on aspen ( Populus tremula ) and goat willow ( Salix caprea ) in Fennoscandian boreal old-growth forests. We estimated the parameters of a forest landscape model (FIN-LANDIS) by repeatedly running simulations and selecting the set of parameters for tree ecology and fire regime that reproduced empirical host tree density and spatial patterns. Second, we tested which variables were important in epiphyte colonization and estimated the dispersal kernel. Third, we run a metapopulation model for the lichen across the estimated landscape scenarios and selected values for the remaining parameters that reproduced the empirical patterns of epiphyte occurrence. There was little variation in predicted dynamics, occupancy and spatial patterns between replicate metapopulation simulations. However, more data would be required for accurately estimating the parameters of FIN-LANDIS primarily because of the inherent stochasticity in large scale forest fires. Following the beginning of fire suppression in the study area 150 years ago, the model predicts that lichen occupancy first increases but subsequently declines. The lower occupancy in the past than at present is explained by high rate of tree destruction by fires, which increases local extinction rate in patch-tracking metapopulations. In the absence of fires, the occupancy increases because of lower extinction rate, but without forest fires or alternative means of host tree regeneration, the lichen is predicted to go ultimately extinct because of severely reduced density of aspen and goat willow.     

Dispersal network heterogeneity promotes species coexistence in hierarchical competitive communities

Understanding the mechanisms of biodiversity maintenance is a fundamental issue in ecology. The possibility that species disperse within the landscape along differing paths presents a relatively unexplored mechanism by which diversity could emerge. By embedding a classical metapopulation model within a network framework, we explore how access to different dispersal networks can promote species coexistence. While it is clear that species with the same demography cannot coexist stably on shared dispersal networks, we find that coexistence is possible on unshared networks, as species can surprisingly form self‐organised clusters of occupied patches with the most connected patches at the core. Furthermore, a unimodal biodiversity response to an increase in species colonisation rates or average patch connectivity emerges in unshared networks. Increasing network size also increases species richness monotonically, producing characteristic species–area curves. This suggests that, in contrast to previous predictions, many more species can co‐occur than the number of limiting resources.     

Pattern does not equal process: what does patch occupancy really tell us about metapopulation dynamics?

Patch occupancy surveys are commonly used to parameterize metapopulation models. If isolation predicts patch occupancy, this is generally attributed to a balance between distance-dependent recolonization and spatially independent extinctions. We investigated whether similar patterns could also be generated by a process of spatially correlated extinctions following a unique colonization event (analogous to nonequilibrium processes in island biogeography). We simulated effects of spatially correlated extinctions on patterns of patch occupancy among pikas (Ochotona princeps) at Bodie, California, using randomly located extinction disks to represent the likely effects of predation. Our simulations produced similar patterns to those cited as evidence of balanced metapopulation dynamics. Simulations using a variety of disk sizes and patch configurations confirmed that our results are potentially applicable to a broad range of species and sites. Analyses of the observed patterns of patch occupancy at Bodie revealed little evidence of rescue effects and strong evidence that most recolonizations are ephemeral in nature. Persistence will be overestimated if static or declining patterns of patch occupancy are mistakenly attributed to dynamically stable metapopulation processes. Consequently, simple patch occupancy surveys should not be considered as substitutes for detailed experimental tests of hypothesized population processes, particularly when conservation concerns are involved.     

Metapopulation theory for fragmented landscapes

We review recent developments in spatially realistic metapopulation theory, which leads to quantitative models of the dynamics of species inhabiting highly fragmented landscapes. Our emphasis is in stochastic patch occupancy models, which describe the presence or absence of the focal species in habitat patches. We discuss a number of ecologically important quantities that can be derived from the full stochastic models and their deterministic approximations, with a particular aim of characterizing the respective roles of the structure of the landscape and the properties of the species. These quantities include the threshold condition for persistence, the contributions that individual habitat patches make to metapopulation dynamics and persistence, the time to metapopulation extinction, and the effective size of a metapopulation living in a heterogeneous patch network.     

Effects of matrix heterogeneity on animal dispersal: from individual behavior to metapopulation-level parameters

Mounting theoretical and empirical evidence shows that matrix heterogeneity may have contrasting effects on metapopulation dynamics by contributing to patch isolation in nontrivial ways. We analyze the movement properties during interpatch dispersal in a metapopulation of Iberian lynx (Lynx pardinus). On a daily temporal scale, lynx habitat selection defines two types of matrix habitats where individuals may move: open and dispersal habitats (avoided and used as available, respectively). There was a strong and complex impact of matrix heterogeneity on movement properties at several temporal scales (hourly and daily radiolocations and the entire dispersal event). We use the movement properties on the hourly temporal scale to build a simulation model to reconstruct individual dispersal events. The two most important parameters affecting model predictions at both the individual (daily) and metapopulation scales were related to the movement capacity (number of movement steps per day and autocorrelation in dispersal habitat) followed by the parameters representing the habitat selection in the matrix. The model adequately reproduced field estimates of population-level parameters (e.g., interpatch connectivity, maximum and final dispersal distances), and its performance was clearly improved when including the effect of matrix heterogeneity on movement properties. To assume there is a homogeneous matrix results in large errors in the estimate of interpatch connectivity, especially for close patches separated by open habitat or corridors of dispersal habitat, showing how important it is to consider matrix heterogeneity when it is present. Movement properties affect the interaction of dispersing individuals with the landscape and can be used as a mechanistic representation of dispersal at the metapopulation level. This is so when the effect of matrix heterogeneity on movement properties is evaluated under biologically meaningful spatial and temporal scales.     

Viability analyses with habitat-based metapopulation models

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