The effective size of a metapopulation living in a heterogeneous patch network

I analyze stochastic patch occupancy models (SPOMs), which record habitat patches as empty or occupied. A problem with SPOMs has been that if the spatial structure of a heterogeneous habitat patch network is taken into account, the computational effort needed to analyze a SPOM grows as a power of 2n, where n is the number of habitat patches. I propose a computationally feasible approximation method, which approximates the behavior of a heterogeneous SPOM by an "ideal" metapopulation inhabiting a network of identical and equally connected habitat patches. The transformation to the ideal metapopulation is based on weighting the individual patch occupancies by the dynamic values of the habitat patches, which may be calculated from the deterministic mean-field approximation of the original SPOM. Conceptually, the method resembles the calculation of the effective size of a population in the context of population genetics. I demonstrate how the method may be applied to SPOMs with flexible structural assumptions and with spatially correlated and temporally varying parameter values. I apply the method to a real habitat patch network inhabited by the Glanville fritillary butterfly, illustrating that the metapopulation dynamics of this species are essentially driven by temporal variability in the environmental conditions.     

Coherence,conservation and patch‐occupancy analysis

Spatial coherence (synchrony) among subpopulations poses a danger to the metacommunity, as it increases the risk of regional extinction. When this effect is significant, the use of inference techniques based on the stochastic patch occupancy model (SPOM) may be inadequate, since SPOMs assume that each habitat patch is either occupied or empty, thereby neglecting the intra‐patch dynamics. Here we suggest a general classification of the dynamics that allows the identification, in a model‐independent manner, of the regimes where coherence effects are strong. We also present a new technique, based on patch occupancy (presence/absence) data, for identifying the role of spatial coherence in the stabilization of a metapopulation. If the chance of a local extinction grows with the connectivity, this implies that spatial synchronization is too strong and that regional‐scale extinction becomes possible. When this scenario occurs, a decrease in the movement of individuals (habitat fragmentation, reduced dispersal rates) has a positive effect on the sustainability of the spatially distributed population. The results of individual based simulations of a spatially structured population are analyzed with SPOM and the regime where the two‐state approximation fails is identified.     

The limiting behaviour of a mainland-island metapopulation

Stochastic patch occupancy models (SPOMs) are a class of discrete time Markov chains used to model the presence/absence of a population in a collection of habitat patches. This class of model is popular with ecologists due to its ability to incorporate important factors of the habitat patch network such as connectivity and distance between patches as well as heterogeneity in patch characteristics. We present an asymptotic examination of a simple type of SPOM called the mainland-island model. In this model a single patch called the mainland is connected to a large number of smaller patches called islands and each island is only connected to the mainland. We discuss the limiting behaviour of the SPOM as the number of islands increases and the size of the islands decrease relative to the mainland. We demonstrate that a variety of limiting behaviours is possible depending on the scaling of the island size and on the heterogeneity of habitat quality.     

Temporal patch occupancy dynamics of the Siberian flying squirrel in a boreal forest landscape

Ability to predict species distribution in a landscape is of crucial importance for natural resource management and species conservation. Therefore, the understanding of species habitat requirements and spatio-temporal dynamics in occurrence is needed. We examined patch occupancy patterns of the Siberian flying squirrel Pteromys volans in northern Finland across a seven year study period. Forest patches dominated by mature spruce ( Picea abies ) in a study area (375 km²) were surveyed to monitor the presence or absence of the flying squirrel. The patch occupancy pattern was dynamic: about half of the habitat patches were occupied at least once during the study period and more patches were colonised than were abandoned. Patches that were continuously occupied (i.e. occupied during all sample periods) were typically of high quality (based on habitat and landscape characteristics), continuously unoccupied patches were usually of low quality, and intermediate quality patches were occupied intermittently. The variables explaining patch occupancy were similar each year, and a statistical model based on data from the year 2000 also predicted occupancy in 2004 with similar accuracy. However, data from a single survey were inadequate for identifying patches used intermittently by flying squirrels. Despite inconsistent occupancy, these patches may be important for the local persistence of flying squirrels. The dynamic occupancy pattern may thus affect estimates of suitable habitat area and identification of functional patch networks for landscape planning. These results emphasise the need for follow-up studies to better understand population patterns and processes in time.     

From individual behavior to metapopulation dynamics: unifying the patchy population and classic metapopulation models

Spatially structured populations in patchy habitats show much variation in migration rate, from patchy populations in which individuals move repeatedly among habitat patches to classic metapopulations with infrequent migration among discrete populations. To establish a common framework for population dynamics in patchy habitats, we describe an individual-based model (IBM) involving a diffusion approximation of correlated random walk of individual movements. As an example, we apply the model to the Glanville fritillary butterfly (Melitaea cinxia) inhabiting a highly fragmented landscape. We derive stochastic patch occupancy model (SPOM) approximations for the IBMs assuming pure demographic stochasticity, uncorrelated environmental stochasticity, or completely correlated environmental stochasticity in local dynamics. Using realistic parameter values for the Glanville fritillary, we show that the SPOMs mimic the behavior of the IBMs well. The SPOMs derived from IBMs have parameters that relate directly to the life history and behavior of individuals, which is an advantage for model interpretation and parameter estimation. The modeling approach that we describe here provides a unified framework for patchy populations with much movements among habitat patches and classic metapopulations with infrequent movements.     

Influence of Habitat Quality and Patch Size on Occupancy and Persistence in two Populations of the Apollo Butterfly (Parnassius apollo)

Recent studies on butterflies emphasize habitat characteristics together with metapopulation parameters (patch area and isolation) giving a more thorough understanding of processes influencing population persistence and patch occupancy, than either of them alone. We studied a coastal and an archipelago population of the Apollo butterfly (Parnassius apollo) in SW Finland. Larvae were surveyed for four years in both populations. Counting larvae on three consecutive days and temporarily removing them tested the survey accuracy. The removals showed four times higher larval abundance in the archipelago than on the coast. Survey methods were reliable, provided that empty patch status was not based on single visits only, if larval abundance was low. On the coast, large patches, and patches with high host-plant abundance were often occupied. In the archipelago, patches rich in host-plant were often occupied whereas patch area did not affect patch occupancy. In both populations, the probability of patches being occupied for three consecutive years increased with increasing host-plant abundance and patch area. Conservation of P. apollo depends on securing host-plant abundance on large enough patches in both study systems. In these systems, even crude habitat measures prove useful for understanding ecological processes behind observed patterns.     

Quantitative predictions for patch occupancy of capercaillie in fragmented habitats

A major conclusion of studying metapopulation biology is that species conservation should favor regional rather than local population persistence. Regional persistence is tightly linked to size, spatial configuration and quality of habitat patches. Hence it is important for the management of endangered species that priority patches can be identified. We developed a predictive model of patch occupancy by capercaillie, a threatened grouse species, based on a single snapshot of data. We used logistic regression to predict patch occupancy as a function of patch size, isolation, connectivity, relative altitude, and biogeographical area. The probability of a patch being occupied increased with patch size and increasing altitude, and decreased with increasing distance to the next occupied patch. Patch size was the most important predictor although occupied patches varied considerably in size. Our model only uses data on the number, size and spatial configuration of habitat patches. It is a useful tool to designate priority areas for conservation, i.e. large core patches with high resilience in habitat quality, smaller island‐patches that still have high probability of being inhabited or becoming recolonised, and patches functioning as “stepping stones”. If capercaillie is to be preserved, habitat suitability needs to be maintained in a functional network of patches that account for size and inter‐patch distance thresholds as found in this study. We suggest that similar area‐isolation relationships are valid for almost any region within the distribution range of capercaillie. The thresholds for occupancy are however likely to depend on characteristics of the respective landscape. The outcome of our study emphasises the need for future investigations that explore the relationship between patch occupancy, matrix quality and its resistance to dispersing individuals.     

The larval ecology of the butterfly Euphydryas desfontainii (Lepidoptera: Nymphalidae) in SW-Portugal: food plant quantity and quality as main predictors of habitat quality

Corresponding to theory, the persistence of metapopulations in fragmented landscapes depends on the area of suitable habitat patches and their degree of isolation, mediating the individual exchange between habitats. More recently, habitat quality has been highlighted as being equally important. We therefore assess the role of habitat area, isolation and quality for the occupancy of larval stages of the regionally threatened butterfly Euphydryas desfontainii occurring in grassland habitats comprising the host plant Dipsascus comosus. We put a special focus on habitat quality which was determined on two spatial scales: the landscape (among patches) and the within-patch level. On the landscape level, occupancy of caterpillars was determined by a presence-absence analysis at 28 host plant patches. On the within-patch level, oviposition site selection was studied by comparing 159 host plants with egg clutches to a random sample of 253 unoccupied host plants within six habitat patches. The occupancy of caterpillars and presence of egg clutches on host plants was then related to several predictors such as patch size and isolation on the landscape level and host plant characteristics and immediate surroundings on the within patch level. On the landscape level, only host plant abundance was related to the presence of caterpillars, while size and isolation did not differ between occupied and unoccupied patches. However, the weak discrimination of larval stages among patches changed on the within-patch level: here, several microclimatic predictors such as sunshine hours and topography, host plant morphology and phenology as well as further potential host plants in the immediate surroundings of the plant chosen for oviposition strongly determined the presence of egg clutches. We strongly suggest promoting the presence of the host plant in topographically and structurally rich habitat patches to offer potential for microclimatic compensation for a species considered threatened by climate change.     

Colonization and extinction dynamics of an annual plant metapopulation in an urban environment

The metapopulation framework considers that the spatiotemporal distribution of organisms results from a balance between the colonization and extinction of populations in a suitable and discrete habitat network. Recent spatially realistic metapopulation models have allowed patch dynamics to be investigated in natural populations but such models have rarely been applied to plants. Using a simple urban fragmented population system in which favourable habitat can be easily mapped, we studied patch dynamics in the annual plant Crepis sancta (Asteraceae). Using stochastic patch occupancy models (SPOMs) and multi‐year occupancy data we dissected extinction and colonization patterns in our system. Overall, our data were consistent with two distinct metapopulation scenarios. A metapopulation (sensu stricto) dynamic in which colonization occurs over a short distance and extinction is lowered by nearby occupied patches (rescue effect) was found in a set of patches close to the city centre, while a propagule rain model in which colonization occurs from a large external population was most consistent with data from other networks. Overall, the study highlights the importance of external seed sources in urban patch dynamics. Our analysis emphasizes the fact that plant distributions are governed not only by habitat properties but also by the intrinsic properties of colonization and dispersal of species. The metapopulation approach provides a valuable tool for understanding how colonization and extinction shape occupancy patterns in highly fragmented plant populations. Finally, this study points to the potential utility of more complex plant metapopulation models than traditionally used for analysing ecological and evolutionary processes in natural metapopulations.     

Patch occupancy in the endangered butterfly Lycaena helle in a fragmented landscape: effects of habitat quality, patch size and isolation

While there is agreement that both habitat quality and habitat network characteristics (such as patch size and isolation) contribute to the occupancy of patches by any given species, the relative importance of these factors is under debate. This issue is of fundamental ecological importance, and moreover of special concern for conservation biologists aiming at preserving endangered species. Against this background we investigated patch occupancy in the violet copper Lycaena helle, one of the rarest butterfly species in Central Europe, in the Westerwald area (Rhineland-Palatinate, Western Germany). Occupied (n = 102) differed from vacant (n = 128) patches in altitude, size, connectivity, availability of wind shelter, in the abundance of the larval host-plant, in the abundance of a grass species indicating favorable habitat conditions and in the abundance of nitrophilous plants. Overall, patch occupancy was primarily determined by patch size, connectivity and the abundance of the larval host plant, while all other parameters of habitat quality were of subordinate importance. Therefore, our findings suggest that even for extremely sedentary species such as L. helle habitat networks are decisive and—next to the preservation of habitat quality—need to be an integral part of any conservation management for this species.     

The influence of patch demographics on metapopulations, with particular reference to successional landscapes

The classical view of metapopulations relates the regional abundance of a species to the balance between the extinction and colonization dynamics of identical local populations. Species in successional landscapes may represent the most appropriate examples of classical metapopulations. However, Levins‐type metapopulation models do not explicitly separate population loss due to successional habitat change from other causes of extinction. A further complication is that the chance of population loss due to successional habitat change may be related to the age of a patch. I developed simple patch occupancy models to include succession and included consideration of patch age structure to address two related questions: what are the implications of changes in patch demographic rates and when is a move to a structured patch occupancy model justified? Age‐related variation in patch demography could increase or decrease the equilibrium fraction of the available habitat occupied by a species when compared to the predictions of an unstructured model. Metapopulation persistence was enhanced when the age class of patches with the highest species occupancy suffered relatively low losses to habitat succession. Conversely, when the age class of patches with the highest species occupancy also had relatively high successional loss rates, extinction thresholds were higher that would be predicted by a simple unstructured model. Hence age‐related variation in patch successional rate introduces biases into the predictions of simple unstructured models. Such biases can be detected from field surveys of the fraction of occupied and unoccupied patches in each age class. Where a bias is demonstrated, unstructured models will not be adequate for making predictions about the effects of changing parameters on metapopulation size. Thinking in successional terms emphasizes how landscapes might be managed to enhance or reduce the patch occupancy by any particular metapopulation     

The contribution of patch topology and demographic parameters to population viability analysis predictions: the case of the European tree frog

Population viability analyses (PVA) are increasingly used in metapopulation conservation plans. Two major types of models are commonly used to assess vulnerability and to rank management options: population-based stochastic simulation models (PSM such as RAMAS or VORTEX) and stochastic patch occupancy models (SPOM). While the first set of models relies on explicit intrapatch dynamics and interpatch dispersal to predict population levels in space and time, the latter is based on spatially explicit metapopulation theory where the probability of patch occupation is predicted given the patch area and isolation (patch topology). We applied both approaches to a European tree frog (Hyla arborea) metapopulation in western Switzerland in order to evaluate the concordances of both models and their applications to conservation. Although some quantitative discrepancies appeared in terms of network occupancy and equilibrium population size, the two approaches were largely concordant regarding the ranking of patch values and sensitivities to parameters, which is encouraging given the differences in the underlying paradigms and input data.     

A simple landscape-scale test of a spatially explicit population model: patch occupancy in fragmented south-eastern Australian forests

The results of a landscape‐scale test of ALEX, a widely used metapopulation model for Population Viability Analysis (PVA), are described. ALEX was used to predict patch occupancy by the laughing kookkaburra and the sacred kingfisher in patches of eucalypt forest in south‐eastern Australia. These predictions were compared to field surveys to determine the accuracy of the model. Predictions also were compared to a “naïve” null model assuming no fragmentation effects.
The naïve null model significantly over‐predicted the number of eucalypt patches occupied by the sacred kingfisher, but the observed patch occupancy was not significantly different from that predicted using ALEX. ALEX produced a better fit to the field data than the naïve null model for the number of patches occupied by the laughing kookaburra. Nevertheless, ALEX still significantly over‐predicted the number of occupied patches, particularly remnants dominated by certain forest types – ribbon gum and narrow‐leaved peppermint. The predictions remained significantly different from observations, even when the habitat quality of these patches was reduced to zero. Changing the rate of dispersal improved overall predicted patch occupancy, but occupancy rates for the different forest types remained significantly different from the field observations. The lack of congruence between field data and model predictions could have arisen because the laughing kookaburra may move between an array of patches to access spatially separated food and nesting resources in response to fragmentation. Alternatively, inter‐specific competition may be heightened in a fragmented habitat. These types of responses to fragmentation are not incorporated as part of traditionally applied metapopulation models. Assessments of predictions from PVA models are rare but important because they can reveal the types of species for which forecasts are accurate and those for which they are not. This can assist the collection of additional empirical data to identify important factors affecting population dynamics.     

Monitoring Golden-Cheeked Warblers on Private Lands in Texas

ABSTRACT A majority of North American breeding habitat for neotropical migrants exists on private lands, requiring monitoring strategies focused on habitat in these private holdings. We outline study designs and protocols using repeated Presence-Absence surveys across a gradient of patch sizes to develop a range-wide monitoring program for the endangered golden-cheeked warbler (Dendroica chrysoparia) in Texas, USA. We surveyed 200–400 point-count locations across approximately 30 private properties annually from 2005 to 2008. We used data from our surveyed patches (n = 147) and the Ψ (occupancy), p (detection), and γ = 1 - ɛ parameterization to estimate patch dynamics and associated detection probabilities for golden-cheeked warblers. Patch size had a strong association with patch occupancy, and all patches >160 ha were predicted to be occupied. We found no evidence that large golden-cheeked warbler populations located on public lands in the vicinity of our study area influenced occupancy dynamics. We conducted simulations across a range of detection probabilities to evaluate potential sample sizes for both standard- and removal-based occupancy modeling. Simulations using parameter estimates from our analysis indicated that removal-based sampling is superior to standard sampling. Based on our results, surveying golden-cheeked warbler presence in oak-juniper (Quercus-Juniperus) patches under a removal modeling framework should be considered as one alternative for range-wide monitoring programs because patch-level monitoring would be necessary to estimate proportion of range occupied. Large contiguous patches are rare across the species’ range; hence, conservation and management of the mosaic of smaller patches within a landscape context would be required for maintaining species viability. Thus, we recommend the identification of areas where smaller, contiguous patches represent a significant portion of the available habitat within the local landscape and targeting these areas for habitat maintenance and improvement.     

Thyme and isolation for the Sinai baton blue butterfly (<Emphasis Type="Italic">Pseudophilotes sinaicus</Emphasis>)

The distribution of the narrowly endemic butterfly Pseudophilotes sinaicus (Lycaenidae) was studied. Potential habitat within its range was first located and then the quality of that habitat assessed. Degree of shelter, diversity of plant species, and resource area of an individual food plant (Thymus decussatus) all affected habitat quality and together were used to develop an index of habitat suitability applicable to each site. The butterfly's distribution was then studied within the identified network of suitable habitat patches: isolated patches with a small resource area were least likely to contain butterflies. Population size in a patch (as opposed merely to patch occupancy) was affected by resource area and the quality of habitat within that patch. Metapopulation processes and variation in habitat quality therefore appear to combine to describe the distribution of patches occupied by P. sinaicus and their population sizes. This finding provides insights into some of the processes operating on an endemic species throughout its geographical range and has important implications for the conservation of this rare butterfly.     

The importance of forest patch networks for the conservation of the Thorn-tailed Rayaditos in central Chile

Conservation of forest birds in fragmented landscapes requires not only determining the critical patch characteristics influencing local population persistence but also identifying patch networks providing connectivity and suitable habitat conditions necessary to ensure regional persistence. In this study, we assessed the importance of patch attributes, patch connectivity, and network components (i.e., groups of interconnected patches) in explaining the occupancy pattern of the Thorn-tailed Rayadito (Aphrastura spinicauda), a forest bird species of central Chile. Using a daily movement threshold distance, we identified a total of 16 network components of sclerophyllous forest within the study area. Among those components, patch area and vegetation structure-composition were important predictors of patch occupancy. However, the inclusion of patch connectivity and component size (i.e., the area of a network component) into the models greatly increases the models’ accuracy and parsimony. Using the best-fitted model, a total of 33 patches were predicted to be occupied by rayaditos within the study area, but such occupied patches were distributed in only six network components. These results suggest that persistence of rayaditos in central Chile requires the maintenance of large single patches and patch networks providing habitat and connectivity.     

Extinction, colonization, and species occupancy in tidepool fishes

Despite the increasing sophistication of ecological models with respect to the size and spatial arrangement of habitat, there is relatively little empirical documentation of how species dynamics change as a function of habitat size and the fraction of habitat occupied. In an assemblage of tidepool fishes, I used maximum-likelihood estimation to test whether models which included habitat size provided a better fit to empirical data on extinction and colonization probabilities than models that assumed constant probabilities over all habitats. I found species differences in how extinction and colonization probabilities scaled with habitat size (and hence local population size). However, there was little evidence for a relationship between extinction and colonization probabilities and the fraction of occupied tidepools, as assumed in simple metapopulation models. Instead, colonization and extinction were independent of the fraction of occupied tidepools, favoring a MacArthur-Wilson island-mainland model. When I incorporated declines in extinction probability with tidepool volume in a simple simulation model, I found that predicted occupancy could change greatly, especially when colonization was low. However, the predicted fraction of occupied patches in the simulation model changed little when I incorporated the range of values reported here for extinction and colonization and the rate at which they scale with habitat size. Quantifying extinction and colonization patterns of natural populations is fundamental to understanding how species are distributed spatially and whether metapopulation models of species occupancy provide explanatory power for field populations. Received: 14 March 1997 / Accepted: 21 September 1997     

Monitoring Svalbard rock ptarmigan: Distance sampling and occupancy modeling

The only resident terrestrial herbivorous bird species in high-Arctic Svalbard, Norway is the endemic Svalbard rock ptarmigan (Lagopus muta hyperborea) of which little is known of its population dynamics. We assessed temporal and spatial variability of the pre-breeding population of Svalbard rock ptarmigan males using: 1) distance sampling to estimate density (2000–2009) and 2) occupancy modeling to determine the proportion of survey points being occupied in relation to a habitat index for ptarmigan habitat suitability (2005–2009). Data were collected using a point-transect sampling design. We split the analysis according to type of survey point (non-random, random, and survey points combined). Our estimated spring densities were low (1.3–3.1 territorial male/km², non-random survey points, 2000–2009) with limited annual variability. The best models describing occupancy rates of territorial males at 2 different spatial scales (ptarmigan males observed ≤250 m and ≤450 m from the sampling point) were independent of spatial scales and the type of survey points. Occupancy dynamics were related to the habitat index whereas detection probability was year dependent. Extinction probability was negatively related to habitat quality (good habitats had lower extinction probability). We could not estimate the habitat effect on colonization precisely because initial occupancy rates were high at both spatial scales (estimated average initial occupancy at scale ≤250 m = 0.96; scale ≤450 m = 0.97). Colonization appeared to be positively related to the habitat index for the random survey points (including mainly marginal habitats), but the small sample size led to large uncertainty in the parameter estimate. Detection probabilities varied greatly between study years, thus demonstrating the importance of estimating detection probability annually. We recommend that future surveys are stratified with respect to habitat quality and to integrate the 2 methodologies in population monitoring of Svalbard rock ptarmigan. © 2011 The Wildlife Society.     

Patch occupancy of North American mammals: is patchiness in the eye of the beholder?

Aim Intraspecific variation in patch occupancy often is related to physical features of a landscape, such as the amount and distribution of habitat. However, communities occupying patchy environments typically exhibit non‐random distributions in which local assemblages of species‐poor patches are nested subsets of assemblages occupying more species‐rich patches. Nestedness of local communities implies interspecific differences in sensitivity to patchiness. Several hypotheses have been proposed to explain interspecific variation in responses to patchiness within a community, including differences in (1) colonization ability, (2) extinction proneness, (3) tolerance to disturbance, (4) sociality and (5) level of adaptation to prevailing environmental conditions. We used data on North American mammals to compare the performance of these ‘ecological’ hypotheses and the ‘physical landscape’ hypothesis. We then compared the best of these models against models that scaled landscape structure to ecologically relevant attributes of individual species. Location North America. Methods We analysed data on prevalence (i.e. proportion of patches occupied in a network of patches) and occupancy for 137 species of non‐volant mammals and twenty networks consisting of four to seventy‐five patches. Insular and terrestrial networks exhibited significantly different mean levels of prevalence and occupancy and thus were analysed separately. Indicator variables at ordinal and family levels were included in models to correct for effects caused by phylogeny. Akaike's information criterion was used in conjunction with ordinary least squares and logistic regression to compare hypotheses. Results A patch network's physical structure, indexed using patch area and isolation, received the greatest support among models predicting the prevalence of species on insular networks. Niche breadth (diet and habitat) received the greatest support for predicting prevalence of species occupying terrestrial networks. For both insular and terrestrial systems, physical features (patch area and isolation) received greater support than any of the ecological hypotheses for predicting species occupancy of individual patches. For terrestrial systems, scaling patch area by its suitability to a focal species and by individual area requirements of the species, and scaling patch isolation by species‐specific dispersal ability and niche breadth, resulted in models of patch occupancy that were superior to models relying solely on physical landscape features. For all selected models, unexplained levels of variation were high. Main conclusions Stochasticity dominated the systems we studied, indicating that random events are probably quite important in shaping local communities. With respect to deterministic factors, our results suggest that forces affecting species prevalence and occupancy may differ between insular and terrestrial systems. Physical features of insular systems appeared to swamp ecological differences among species in determining prevalence and occupancy, whereas species with broad niches were disproportionately represented in terrestrial networks. We hypothesize that differential extinction over long time periods in highly variable networks has driven nestedness of mammalian communities on islands, whereas differential colonization over shorter time‐scales in more homogeneous networks probably governed the local structure of terrestrial communities. Our results also demonstrate that integration of a species' ecological traits with physical features of a patch network is superior to reliance on either factor separately when attempting to predict the species' probability of patch occupancy in terrestrial systems.     

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.