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

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

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

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

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.     

Meso-scale hydrodynamic and reproductive asynchrony affects the source–sink metapopulation structure of the coastal polychaete Pectinaria koreni

The polychaete Pectinaria koreni exhibits a complex life cycle characterized by non-overlapping generations and widespread larval dispersal. To explore how “local” metapopulation genetic structure varies spatially and temporally during population turnover, we combined observations on demography, larval dispersal through hydrodynamic modelling and population genetics of successive age cohorts in the Baie de Seine (eastern English Channel, France). Mature adults (March), newly settled (July) and later-stage juveniles (September) were sampled in 2003 on the edge and in the main demes of the metapopulation. Demes displayed an asynchronous dynamics due to variations in habitat quality affecting reproductive timing (e.g. three distinct spawning events observed) and in local larval supply linked to temporal fluctuations of hydrodynamism. Two-source populations were identified among dense areas with the greatest larval retention and self-recruitment rates: one with a single recruitment event, stable temporal genetic variation and a strong spatial genetic re-homogenization during turnover, and the other with two recruitment events and significant allele frequency changes over time. Sink habitats displayed single recruitment event and experienced strong inter-generational (juveniles vs adults) genetic changes due to genetic drift associated with strong winter mortalities. Altogether, results suggested that adult spawning asynchrony and sweepstakes reproductive success, together with genetic drift, played a greater role than hydrodynamics itself in determining effective recruitment rates at some sites and generating genetic patchiness.     

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.     

Cougar Dispersal Patterns, Metapopulation Dynamics, and Conservation

Abstract: We examined cougar (   Puma concolor ) dispersal, emigration, and immigration in the San Andres Mountains, New Mexico, from 1985 to 1995 to quantify the effects of dispersal on the local population and surrounding subpopulations. We captured, tagged, and radio-collared animals to detect the arrival of new immigrants and dispersal characteristics of progeny. We found that cougars in southern New Mexico exhibited a metapopulation structure in which cougar subpopulations were separated by expanses of noncougar habitat and linked by dispersers. Of 43 progeny (n = 20 males , 23 females ) studied after independence, only 13 females exhibited philopatric behavior. Males dispersed significantly farther than females, were more likely to traverse large expanses of noncougar habitat, and were probably most responsible for nuclear gene flow between habitat patches. We estimated that an average of 8.5 progeny (i.e., cougars born in the study area) successfully emigrated from and 4.3 cougars successfully immigrated to the San Andres Mountains each year. Concurrently, an average of 4.1 progeny were recruited into the San Andres cougar population. Protected cougar subpopulations can contribute to metapopulation persistence by supplying immigrants to surrounding subpopulations that are affected by fragmentation or offtake by humans. Cougar population dynamics and dispersal behavior dictate that cougar management and conservation should be considered on a regional scale.     

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.     

Why are metapopulations so rare?

Roughly 40 years after its introduction, the metapopulation concept is central to population ecology. The notion that local populations and their dynamics may be coupled by dispersal is without any doubt of great importance for our understanding of population-level processes. A metapopulation describes a set of subpopulations linked by (rare) dispersal events in a dynamic equilibrium of extinctions and recolonizations. In the large body of literature that has accumulated, the term "metapopulation" is often used in a very broad sense; most of the time it simply implies spatial heterogeneity. A number of reviews have recently addressed this problem and have pointed out that, despite the large and still growing popularity of the metapopulation concept, there are only very few empirical examples that conform with the strict classical metapopulation (CM) definition. In order to understand this discrepancy between theory and observation, we use an individual-based modeling approach that allows us to pinpoint the environmental conditions and the life-history attributes required for the emergence of a CM structure. We find that CM dynamics are restricted to a specific parameter range at the border between spatially structured but completely occupied and globally extinct populations. Considering general life-history attributes, our simulations suggest that CMs are more likely to occur in arthropod species than in (large) vertebrates. Since the specific type of spatial population structure determines conservation concepts, our findings have important implications for conservation biology. Our model suggests that most spatially structured populations are panmictic, patchy, or of mainland-island type, which makes efforts spent on increasing connectivity (e.g., corridors) questionable. If one does observe a true CM structure, this means that the focal metapopulation is on the brink of extinction and that drastic conservation measures are needed.     

Information processing in models for density-dependent emigration: A comparison

Density-dependent emigration has been recognized as a fitness enhancing strategy. Yet, especially in the modelling literature there is no consensus about how density-dependent emigration should quantitatively be incorporated into metapopulation models. In this paper we compare the performance of five different dispersal strategies (defined by the functional link between density and emigration probability). Four of these strategies are based on published functional relationships between local population density and emigration probability, one assumes density-independent dispersal. We use individual-based simulations of time-discrete metapopulation dynamics and conduct evolution experiments for a broad range of values for dispersal mortality and environmental stochasticity. For each set of these conditions we analyze the evolution of emigration rates in ‘monoculture experiments’ (with only one type of dispersal strategy used by all individuals in the metapopulation) as well as in selection experiments that allow a pair-wise comparison of the performance of each functional type. We find that a single-parameter ‘asymptotic threshold’ strategy - derived from the marginal value theorem - with a decelerating increase of emigration rate with increasing population density, out-competes any other strategy, i.e. density-independent emigration, a ‘linear threshold’ strategy and a flexible three-parameter strategy. Only when environmental conditions select for extremely high emigration probabilities (close to one), strategies may perform approximately equally. A simple threshold strategy derived for the case of continuous population growth performs even worse than the density-independent strategy. As the functional type of the dispersal function implemented in metapopulation models may severely affect predictions concerning the survival of populations, range expansion, or community changes we clearly recommend to carefully select adequate functions to model density-dependent dispersal.     

Spatiotemporal dynamics in variable population interactions with density-dependent interaction coefficients

In this article, I present a two-patch metapopulation model with locally explicit dynamics to study the effect of spatial heterogeneity and dispersal upon population interactions with variable or conditional outcomes. These are interactions that may be either detrimental or beneficial for each species depending on the balance of the density-dependent costs and benefits involved. The local dynamics respond to density-dependent α-interaction functions that may change sign, thus yielding a diversity of possible local outcomes for the association in terms of type of interaction and in the number of stable solutions. The spatiotemporal model predicts that the fragmentation of space and dispersal between patches may cause further variation in these outcomes. First, the demographic performance of a species in the association is enhanced if migrations cause a proportional increase of individuals of its own species; being so, a victim may become a mutualist or an exploiter, an excluded species may invade, and a good competitor may overcome its own carrying capacity: the ‘enhancement effect of dispersal’; a sort of rescue effect in source-sink dynamics. The underlying mechanisms involve an interplay between density-dependent effects of dispersal per se and the relative local and global average α-interaction functions, which involve costs and benefits at both the local and regional level that may either counteract or reinforce each other; thus, localities and/or populations may change dynamically their sink or source role in the spatial dynamics. A significant insight arises herewith: in the context of variable or conditional interactions the concept of the role of a species does not make strict sense; it becomes a spatiotemporal dynamic quality. Second, regardless of which species disperses, bifurcation of equilibria may occur in those patches that receive the migrating individuals, and annihilation of equilibria in those from where migration leaves; thus, the number of equilibria increases or decreases accordingly.     

Risk spreading, connectivity, and optimal reserve spacing

Two important processes determining the dynamics of spatially structured populations are dispersal and the spatial covariance of demographic fluctuations. Spatially explicit approaches to conservation, such as reserve networks, must consider the tension between these two processes and reach a balance between distances near enough to maintain connectivity, but far enough to benefit from risk spreading. Here, we model this trade-off. We show how two measures of metapopulation persistence depend on the shape of the dispersal kernel and the shape of the distance decay in demographic covariance, and we consider the implications of this trade-off for reserve spacing. The relative rates of distance decay in dispersal and demographic covariance determine whether the long-run metapopulation growth rate, and quasi-extinction risk, peak for adjacent patches or intermediately spaced patches; two local maxima in metapopulation persistence are also possible. When dispersal itself fluctuates over time, the trade-off changes. Temporal variation in mean distance that propagules are dispersed (i.e., propagule advection) decreases metapopulation persistence and decreases the likelihood that persistence will peak for adjacent patches. Conversely, variation in diffusion (the extent of random spread around mean dispersal) increases metapopulation persistence overall and causes it to peak at shorter inter-patch distances. Thus, failure to consider temporal variation in dispersal processes increases the risk that reserve spacings will fail to meet the objective of ensuring metapopulation persistence. This study identifies two phenomena that receive relatively little attention in empirical work on reserve spacing, but that can qualitatively change the effectiveness of reserve spacing strategies: (1) the functional form of the distance decay in covariance among patch-specific demographic rates and (2) temporal variation in the shape of the dispersal kernel. The sensitivity of metapopulation recovery and persistence to how covariance of vital rates decreases with distance suggests that estimating the shape of this function is likely to be as important for effective reserve design as estimating connectivity. Similarly, because temporal variation in dispersal dynamics influences the effect of reserve spacing, approaches to reserve design that ignore such variation, and rely instead on long-term average dispersal patterns, are likely to lead to lower metapopulation viability than is actually achievable.     

Local Population Dynamics in Metapopulation Models: Implications for Conservation

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

Consistent scaling of persistence time in metapopulations

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

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.     

Regional Consequences of Local Population Demography and Genetics in Relation to Habitat Management in Gentiana pneumonanthe

Abstract:  A joint demographic and population genetics stage-based model for a subdivided population was applied to Gentiana pneumonanthe , an early successional perennial herb, at a regional (metapopulation) scale. We used numerical simulations to determine the optimal frequency of habitat disturbance (sod cutting) and the intensity of gene flow among populations of G. pneumonanthe to manage both population viability and genetic diversity in this species. The simulations showed that even small populations that initially had near-equal allele frequencies could, if managed properly through sod cutting every 6 to 7 years, sustain their high genetic variation over the long run without gene flow. The more the allele frequencies in the small populations are skewed, however, the higher the probability that in the absence of gene flow, some alleles will be lost and within-population genetic variation will decrease even under proper management. This implies that although local population dynamics should be the major target for management, regional dynamics become important when habitat fragmentation and decreased population size lead to the loss of local genetic diversity. The recommended strategy to improve genetic composition of small populations is the introduction of seeds or seedlings of nonlocal origin.     

Integrating Landscape and Metapopulation Modeling Approaches: Viability of the Sharp-Tailed Grouse in a Dynamic Landscape

Abstract:  The lack of management experience at the landscape scale and the limited feasibility of experiments at this scale have increased the use of scenario modeling to analyze the effects of different management actions on focal species. However, current modeling approaches are poorly suited for the analysis of viability in dynamic landscapes. Demographic (e.g., metapopulation) models of species living in these landscapes do not incorporate the variability in spatial patterns of early successional habitats, and landscape models have not been linked to population viability models. We link a landscape model to a metapopulation model and demonstrate the use of this model by analyzing the effect of forest management options on the viability of the Sharp-tailed Grouse (  Tympanuchus phasianellus ) in the Pine Barrens region of northwestern Wisconsin (U.S.A.). This approach allows viability analysis based on landscape dynamics brought about by processes such as succession, disturbances, and silviculture. The landscape component of the model (LANDIS) predicts forest landscape dynamics in the form of a time series of raster maps. We combined these maps into a time series of patch structures, which formed the dynamic spatial structure of the metapopulation component (RAMAS). Our results showed that the viability of Sharp-tailed Grouse was sensitive to landscape dynamics and demographic variables such as fecundity and mortality. Ignoring the landscape dynamics gave overly optimistic results, and results based only on landscape dynamics (ignoring demography) lead to a different ranking of the management options than the ranking based on the more realistic model incorporating both landscape and demographic dynamics. Thus, models of species in dynamic landscapes must consider habitat and population dynamics simultaneously.     

The effect of summer drought on the predictability of local extinctions in a butterfly metapopulation

The ecological impacts of extreme climatic events on population dynamics and community composition are profound and predominantly negative. Using extensive data of an ecological model system, we tested whether predictions from ecological models remain robust when environmental conditions are outside the bounds of observation. We observed a 10-fold demographic decline of the Glanville fritillary butterfly (Melitaea cinxia) metapopulation on the Åland islands, Finland in the summer of 2018 and used climatic and satellite data to demonstrate that this year was an anomaly with low climatic water balance values and low vegetation productivity indices across Åland. Population growth rates were strongly associated with spatiotemporal variation in climatic water balance. Covariates shown previously to affect the extinction probability of local populations in this metapopulation were less informative when populations were exposed to severe drought during the summer months. Our results highlight the unpredictable responses of natural populations to extreme climatic events.     

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

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

Applying Metapopulation Theory to Conservation of Migratory Birds

Abstract: Metapopulation theory has proven useful for understanding the population structure and dynamics of many species of conservation concern. The metapopulation concept has been applied almost exclusively to nonmigratory species, however, for which subpopulation demographic independence—a requirement for a classically defined metapopulation—is explicitly related to geographic distribution and dispersal probabilities. Defining the degree of demographic independence among subpopulations of migratory animals, and thus the applicability of metapopulation theory as a conceptual framework for understanding population dynamics, is much more difficult. Unlike nonmigratory species, subpopulations of migratory animals cannot be defined as synonymous with geographic areas. Groups of migratory birds that are geographically separate at one part of the annual cycle may occur together at others, but co-occurrence in time and space does not preclude the demographic independence of subpopulations. I suggest that metapopulation theory can be applied to migratory species but that understanding the degree of subpopulation independence may require information about both spatial distribution throughout the annual cycle and behavioral mechanisms that may lead to subpopulation demographic independence. The key for applying metapopulation theory to migratory animals lies in identifying demographically independent subpopulations, even as they move during the annual cycle and potentially co-occur with other subpopulations. Using examples of migratory bird species, I demonstrate that spatial and temporal modes of subpopulation independence can interact with behavioral mechanisms to create demographically independent subpopulations, including cases in which subpopulations are not spatially distinct in some parts of the annual cycle.     

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.