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     

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.     

Do invasive species undergo metapopulation dynamics? A case study of the invasive Caspian gull,Larus cachinnans,in Poland

Aim The mechanisms of initial dispersal and habitat occupancy by invasive alien species are fundamental ecological problems. Most tests of metapopulation theory are performed on local population systems that are stable or in decline. In the current study we were interested in the usefulness of metapopulation theory to study patch occupancy, local colonization, extinction and the abundance of the invasive Caspian gull (Larus cachinnans) in its initial invasion stages. Location Waterbodies in Poland. Methods Characteristics of the habitat patches (waterbodies, 35 in total) occupied by breeding pairs of Caspian gulls and an equal sample of randomly selected unoccupied patches were compared with t‐tests. Based on presence–absence data from 1989 to 2006 we analysed factors affecting the probability of local colonization, extinction and the size of local populations using generalized linear models. Results Occupied habitat patches were significantly larger and less isolated (from other habitat patches and other local populations) and were located closer to rivers than empty patches. The proximity of local food resources (fish ponds, refuse dumps) positively affected the occurrence of breeding pairs. The probability of colonization was positively affected by patch area, and negatively by distances to fish ponds, nearest habitat patch, nearest breeding colony and to a river, and by higher forest cover around the patch boundaries. The probability of extinction was lower in patches with a higher number of breeding pairs and with a greater area of islets. The extinction probability increased with distances to other local populations, other habitat patches, fish ponds and to refuse dumps and with a higher cover of forest around the patch boundaries. The size of the local population decreased with distances to the nearest habitat patch, local population, river, fish pond and refuse dump. Local abundance was also positively affected by the area of islets in the patch. Main conclusions During the initial stages of the invasion of Caspian gulls in Poland the species underwent metapopulation‐like dynamics with frequent extinctions from colonized habitat patches. The results prove that metapopulation theory may be a useful conceptual framework for predicting which habitats are more vulnerable to invasion.     

Asynchronization of local population dynamics and persistence of a metapopulation: a lesson from an endangered composite plant, Aster kantoensis

Fragmentation of a large habitat makes local populations less linked to others, and a whole population structure changes to a metapopulation. The smaller a local population is, the more strengthened extinction factors become. Then, frequent extinctions of local populations threaten persistence of the metapopulation unless recolonizations occur rapidly enough after local extinctions. Spatially structured models have been more widely used for predicting future population dynamics and for assessing the extinction risk of a metapopulation. In this article, we first review such spatially structured models that have been applied to conservation biology, focusing on effects of asynchronization among local population dynamics on persistence of the whole metapopulation. Second, we introduce our ongoing project on extinction risk assessment of an endangered composite biennial plant, Aster kantoensis, in the riverside habitat, based on a lattice model for describing its spatiotemporal population dynamics. The model predicted that the extinction risk of A. kantoensis depends on both the frequency of flood occurrence and the time to coverage of a local habitat by other competitively stronger perennials. Finally, we present a measure (Hassell and Pacala's CV ²) for quantifying the effect of asynchronization among local population dynamics on the persistence of a whole metapopulation in conservation ecology. Received: January 12, 2000 / Accepted: February 8, 2000     

Population turnover and habitat dynamics in Notonecta (Hemiptera: Notonectidae) metapopulations

Simple metapopulation models assume that local populations occur in patches of uniform quality habitat separated by non-habitat. However field metapopulations tend to show considerable spatial and temporal variation in patch quality, and hence probability of occupancy. This may have implications for the adequacy of simple metapopulation models in describing and predicting regional population dynamics of natural systems. This study investigated the effects of habitat characteristics on landscape-scale occupancy dynamics of two species of backswimmer (Notonecta, Hemiptera: Notonectidae) in small freshwater ponds. The results demonstrated clear links between habitat, pond occupancy and population turnover, particularly local extinction. There were considerable changes in the habitat of individual ponds between years, but local changes were not spatially correlated and the frequency distribution of habitat conditions at the landscape level remained similar in different years. Stable occupancy levels of Notonecta species appears to result from a balance of the rates of creation and loss of suitable habitat due to spatially uncorrelated habitat change. Systems such as this, where turnover is driven by habitat dynamics, demonstrate the potential value of incorporating the dynamics of habitat change into metapopulation models. Such developments are likely to improve predictions of landscape-scale occupancy dynamics, whilst also allowing patch-level predictions of occupancy, based on local habitat conditions. Received: 18 August 1999 / Accepted: 3 December 1999     

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     

Patterns of commonness and rarity in central European birds: reliability of the core-satellite hypothesis within a large scale

The frequency distribution of species’ area of occupancy is often bimodal, most species being either very rare or very common in terms of number of occupied sites. This pattern has been attributed to the nonlinearity associated with metapopulation dynamics of the species, but there are also other explanations comprising sampling artifact and frequency distribution of suitable habitats. We tested whether the bimodal frequency distribution of occupied squares in central European birds could be derived solely from the frequency distribution of species population sizes (i.e. the sampling artifact hypothesis) or from the spatial distribution of their preferred habitats. Both models predict high proportion of very common species, i.e. the right side of frequency distribution. Bimodality itself is well predicted by models based on random placement of individuals according to their abundances but neither model predicts the observed prevalence of rare species. Even the combined models that assume random placement of individuals within the squares with suitable habitat do not predict such a high proportion of rare species. The observed distribution is more aggregated, rare species occupying a smaller portion of suitable habitat than predicted on the basis of their abundance. The pattern is consistent with metapopulation processes involving local population extinctions. The involvement of these processes is supported by two further observations. First, species rarity is associated with significant population trend and/or location on the edge of their ranges within central Europe, both situations presumably associated with metapopulation processes. Second, suitable habitats seem to be either saturated or almost unoccupied, which is consistent with the predictions of the metapopulation model based on nonlinear dynamics of extinction and colonization. Although the habitat suitability is an important determinant of species distribution, the rarity of many species of birds within this scale of observation seems to be affected by other factors, including local population extinctions associated with fragmentation of species’ habitats.     

Local and Regional Determinants of Colonisation-Extinction Dynamics of a Riparian Mainland-Island Root Vole Metapopulation

The role of local habitat geometry (habitat area and isolation) in predicting species distribution has become an increasingly more important issue, because habitat loss and fragmentation cause species range contraction and extinction. However, it has also become clear that other factors, in particular regional factors (environmental stochasticity and regional population dynamics), should be taken into account when predicting colonisation and extinction. In a live trapping study of a mainland-island metapopulation of the root vole (Microtus oeconomus) we found extensive occupancy dynamics across 15 riparian islands, but yet an overall balance between colonisation and extinction over 4 years. The 54 live trapping surveys conducted over 13 seasons revealed imperfect detection and proxies of population density had to be included in robust design, multi-season occupancy models to achieve unbiased rate estimates. Island colonisation probability was parsimoniously predicted by the multi-annual density fluctuations of the regional mainland population and local island habitat quality, while extinction probability was predicted by island population density and the level of the recent flooding events (the latter being the main regionalized disturbance regime in the study system). Island size and isolation had no additional predictive power and thus such local geometric habitat characteristics may be overrated as predictors of vole habitat occupancy relative to measures of local habitat quality. Our results suggest also that dynamic features of the larger region and/or the metapopulation as a whole, owing to spatially correlated environmental stochasticity and/or biotic interactions, may rule the colonisation – extinction dynamics of boreal vole metapopulations. Due to high capacities for dispersal and habitat tracking voles originating from large source populations can rapidly colonise remote and small high quality habitat patches and re-establish populations that have gone extinct due to demographic (small population size) and environmental stochasticity (e.g. extreme climate events).     

Long-term dynamics in a metapopulation of the American pika

A 20-yr study of a metapopulation of the American pika revealed a regional decline in occupancy in one part of a large network of habitat patches. We analyze the possible causes of this decline using a spatially realistic metapopulation model, the incidence function model. The pika metapopulation is the best-known mammalian example of a classical metapopulation with significant population turnover, and it satisfies closely the assumptions of the incidence function model, which was parameterized with data on patch occupancy. The model-predicted incidences of patch occupancy are consistent with observed incidences, and the model predicts well the observed turnover rate between four metapopulation censuses. According to model predictions, the part of the metapopulation where the decline has been observed is relatively unstable and prone to large oscillations in patch occupancy, whereas the other part of the metapopulation is predicted to be persistent. These results demonstrate how extinction-colonization dynamics may produce spatially correlated patterns of patch occupancy without any spatially correlated processes in local dynamics or extinction rate. The unstable part of the metapopulation gives an empirical example of multiple quasi equilibria in metapopulation dynamics. Phenomena similar to those observed here may cause fluctuations in species' range limits.     

Anthropogenic, ecological and genetic factors in extinction and conservation

Anthropogenic factors constitute the primary deterministic causes of species declines, endangerment and extinction: land development, overexploitation, species translocations and introductions, and pollution. The primary anthropogenic factors produce ecological and genetic effects contributing to extinction risk. Ecological factors include environmental stochasticity, random catastrophes, and metapopulation dynamics (local extinction and colonization) that are intensified by habitat destruction and fragmentation. Genetic factors include hybridization with nonadapted gene pools, and selective breeding and harvesting. In small populations stochastic factors are especially important, including the ecological factors of Allee effect, edge effects, and demographic stochasticity, and the genetic factors of inbreeding depression, loss of genetic variability, and fixation of new deleterious mutations. All factors affecting extinction risk are expressed, and can be evaluated, through their operation on population dynamics.     

Finite metapopulation models with density-dependent migration and stochastic local dynamics

The effects of small density-dependent migration on the dynamics of a metapopulation are studied in a model with stochastic local dynamics. We use a diffusion approximation to study how changes in the migration rate and habitat occupancy affect the rates of local colonization and extinction. If the emigration rate increases or if the immigration rate decreases with local population size, a positive expected rate of change in habitat occupancy is found for a greater range of habitat occupancies than when the migration is density-independent. In contrast, the reverse patterns of density dependence in respective emigration and immigration reduce the range of habitat occupancies where the metapopulation will be viable. This occurs because density-dependent migration strongly influences both the establishment and rescue effects in the local dynamics of metapopulations.     

Population dynamics can be more important than physiological limits for determining range shifts under climate change

Evidence is accumulating that species' responses to climate changes are best predicted by modelling the interaction of physiological limits, biotic processes and the effects of dispersal‐limitation. Using commercially harvested blacklip (Haliotis rubra) and greenlip abalone (Haliotis laevigata) as case studies, we determine the relative importance of accounting for interactions among physiology, metapopulation dynamics and exploitation in predictions of range (geographical occupancy) and abundance (spatially explicit density) under various climate change scenarios. Traditional correlative ecological niche models (ENM) predict that climate change will benefit the commercial exploitation of abalone by promoting increased abundances without any reduction in range size. However, models that account simultaneously for demographic processes and physiological responses to climate‐related factors result in future (and present) estimates of area of occupancy (AOO) and abundance that differ from those generated by ENMs alone. Range expansion and population growth are unlikely for blacklip abalone because of important interactions between climate‐dependent mortality and metapopulation processes; in contrast, greenlip abalone should increase in abundance despite a contraction in AOO. The strongly non‐linear relationship between abalone population size and AOO has important ramifications for the use of ENM predictions that rely on metrics describing change in habitat area as proxies for extinction risk. These results show that predicting species' responses to climate change often require physiological information to understand climatic range determinants, and a metapopulation model that can make full use of this data to more realistically account for processes such as local extirpation, demographic rescue, source‐sink dynamics and dispersal‐limitation.     

Metapopulation dynamics: Does it help to have more of the same?

With the interest in conservation biology shifting towards processes from patterns, and to populations from communities, the theory of metapopulation dynamics is replacing the equilibrium theory of island biogeography as the population ecology paradigm in conservation biology. The simplest models of metapopulation dynamics make predictions about the effects of habitat fragmentation - size and isolation of habitat patches - on metapopulation persistence. The simple models may be enriched by considerations of the effects of demographic and environmental stochasticity on the size and extinction probability of local populations. Environmental stochasticity affects populations at two levels: it makes local extinctions more probable, and it also decreases metapopulation persistence time by increasing the correlation of extinction events across populations. Some controversy has arisen over the significance of correlated extinctions, and how they may affect the optimal subdivision of metapopulations to maximize their persistence time.     

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     

用空间直观模型是否足以从斑块占据性资料中推断集合种群的动态过程?

在集合种群的研究中,经常要根据空间占据性数据应用斑块模型来推断种群的动态过程,在保护生物学应用中,斑块占据性模型的参数估测对于阐释集合种群动态和预测种群对生境破坏的反应极为重要。我们探讨了一种广泛应用的空间直观模型——率函数模型(Incidence function model)中参数估测的不确定性问题,通过构建由50个斑块组成的网络和两个假想的已知参数的集合种群,应用模拟模型产生集合种群随时间变化的斑块占据性数据系列：即快照(snapshot)。然后,根据这些快照,应用率函数模型和最大似然法估测种群动态参数。此外,我们还给出了传统的率函数模型的一个变形,这个变形包含了目标区效应(Target area effect)：即一个斑块的占据概率不但取决于空间隔离度,也取决于斑块本身面积的大小。结果表明：根据同一个集合种群不同的快照所估测的参数可以有很大差异,一个快照得出的参数提示的是占据性强但存活率低的集合种群,而另一个快照可能反映的是一个占据性弱但存活率高的集合种群。应用传统的率函数模型于一个包含了目标区效应的集合种群,导致斑块大小相关的灭绝率参数估测的正偏差。因此,仅根据一个快照的空间占据性数据来推测集合种群的过程有很大的不确定性[动物学报49(6)：787～794,2003]。     

Dynamic impacts of feral mink predation on vole metapopulations in the outer archipelago of the Baltic Sea

Predation impacts by introduced predators are predicted to be most intense in island ecosystems, and also variable depending on environmental conditions, but large-scale experimental field testing is rare. In this study we examine the factors that determine the distribution and abundance of vole metapopulations preyed upon by feral American mink Mustela vison in the outer Finnish archipelago of the Baltic Sea. Specifically, we follow the dynamics of field voles Microtus agrestis and bank voles Clethrionomys glareolus on 40 small islands under variable rainfall as part of a large-scale mink removal experiment. For both vole species occupancy rates were negatively influenced by island isolation, as were extinction events for field voles. High summer rainfall in 1998 corresponded to large vole populations where mink were absent, populations that then crashed in 1999 and 2000 when below average rains fell during the summer breeding season. Where mink were present however, vole abundance remained more constant between years with no boom-bust apparent. We conclude that weather and predation may drive vole abundance whereas habitat patchiness and metapopulation processes more strongly drive vole distributions. There may also be potential for interaction between these factors: because feral mink prevent rapid vole population growth after good summer rains, and vole dispersal is influenced by population size, feral mink may be changing vole dispersal patterns to disrupt the natural metapopulation dynamic. Hence this indirect impact of mink could lead to gradual erosion of vole populations in the outer archipelago by reducing recolonisation processes.     

Metapopulation dynamics with quasi-local competition

Stepping-stone models for the ecological dynamics of metapopulations are often used to address general questions about the effects of spatial structure on the nature and complexity of population fluctuations. Such models describe an ensemble of local and spatially isolated habitat patches that are connected through dispersal. Reproduction and hence the dynamics in a given local population depend on the density of that local population, and a fraction of every local population disperses to neighboring patches. In such models, interesting dynamic phenomena, e.g. the persistence of locally unstable predator-prey interactions, are only observed if the local dynamics in an isolated patch exhibit non-equilibrium behavior. Therefore, the scope of these models is limited. Here we extend these models by making the biologically plausible assumption that reproductive success in a given local habitat not only depends on the density of the local population living in that habitat, but also on the densities of neighboring local populations. This would occur if competition for resources occurs between neighboring populations, e.g. due to foraging in neighboring habitats. With this assumption of quasi-local competition the dynamics of the model change completely. The main difference is that even if the dynamics of the local populations have a stable equilibrium in isolation, the spatially uniform equilibrium in which all local populations are at their carrying capacity becomes unstable if the strength of quasi-local competition reaches a critical level, which can be calculated analytically. In this case the metapopulation reaches a new stable state, which is, however, not spatially uniform anymore and instead results in an irregular spatial pattern of local population abundance. For large metapopulations, a huge number of different, spatially non-uniform equilibrium states coexist as attractors of the metapopulation dynamics, so that the final state of the system depends critically on the initial conditions. The existence of a large number of attractors has important consequences when environmental noise is introduced into the model. Then the metapopulation performs a random walk in the space of all attractors. This leads to large and complicated population fluctuations whose power spectrum obeys a red-shifted power law. Our theory reiterates the potential importance of spatial structure for ecological processes and proposes new mechanisms for the emergence of non-uniform spatial patterns of abundance and for the persistence of complicated temporal population fluctuations.     

Eco-evolutionary dynamics of dispersal in spatially heterogeneous environments

Ecology Letters (2011) 14: 1025-1034 ABSTRACT: Evolutionary changes in natural populations are often so fast that the evolutionary dynamics may influence ecological population dynamics and vice versa. Here we construct an eco-evolutionary model for dispersal by combining a stochastic patch occupancy metapopulation model with a model for changes in the frequency of fast-dispersing individuals in local populations. We test the model using data on allelic variation in the gene phosphoglucose isomerase (Pgi), which is strongly associated with dispersal rate in the Glanville fritillary butterfly. Population-specific measures of immigration and extinction rates and the frequency of fast-dispersing individuals among the immigrants explained 40% of spatial variation in Pgi allele frequency among 97 local populations. The model clarifies the roles of founder events and gene flow in dispersal evolution and resolves a controversy in the literature about the consequences of habitat loss and fragmentation on the evolution of dispersal.     

Simple discrete-time metapopulation models of patch occupancy

Simple models in theoretical ecology have a long-standing history of being used to understand how specific processes influence population dynamics as well as providing a foundation for future endeavors. The Levins model is the seminal example of this for continuous-time metapopulation dynamics. However, many natural populations have a distinct separation between processes and data is not collected continuously leading to the need for using a discrete-time model. Our goal is to develop a simple discrete-time metapopulation model of patch occupancy using difference equations. In our formulation, we consider the two fundamental processes of colonization and extinction that will be treated as sequential events and will only consider patch occupancy. To achieve this, we use a composition of two functions where one will reflect the extinction process and the other for the colonization process. Under some mild assumptions, we are able determine the dynamic behavior of the metapopulation. In addition, we provide numerous examples for the functions used to emulate the colonization and extinction processes. Our results illustrate that the dynamics of the model are tied to properties such as convexity and monotonicity of the colonization and extinction functions. In particular, if the model is non-monotone, then complex dynamics can arise such as cyclic and even chaotic behavior. Overall, our approach shows how certain properties of the colonization and extinction functions can influence metapopulation dynamics.     

Occupancy versus colonization–extinction models for projecting population trends at different spatial scales

Understanding spatiotemporal population trends and their drivers is a key aim in population ecology. We further need to be able to predict how the dynamics and sizes of populations are affected in the long term by changing landscapes and climate. However, predictions of future population trends are sensitive to a range of modeling assumptions. Deadwood‐dependent fungi are an excellent system for testing the performance of different predictive models of sessile species as these species have different rarity and spatial population dynamics, the populations are structured at different spatial scales, and they utilize distinct substrates. We tested how the projected large‐scale occupancies of species with differing landscape‐scale occupancies are affected over the coming century by different modeling assumptions. We compared projections based on occupancy models against colonization–extinction models, conducting the modeling at alternative spatial scales and using fine‐ or coarse‐resolution deadwood data. We also tested effects of key explanatory variables on species occurrence and colonization–extinction dynamics. The hierarchical Bayesian models applied were fitted to an extensive repeated survey of deadwood and fungi at 174 patches. We projected higher occurrence probabilities and more positive trends using the occupancy models compared to the colonization–extinction models, with greater difference for the species with lower occupancy, colonization rate, and colonization:extinction ratio than for the species with higher estimates of these statistics. The magnitude of future increase in occupancy depended strongly on the spatial modeling scale and resource resolution. We encourage using colonization–extinction models over occupancy models, modeling the process at the finest resource‐unit resolution that is utilizable by the species, and conducting projections for the same spatial scale and resource resolution at which the model fitting is conducted. Further, the models applied should include key variables driving the metapopulation dynamics, such as the availability of suitable resource units, habitat quality, and spatial connectivity.