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.     

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.     

A robust-design formulation of the incidence function model of metapopulation dynamics applied to two species of rails

The incidence function model (IFM) uses area and connectivity to predict metapopulation dynamics. However, false absences and missing data can lead to underestimates of the number of sites contributing to connectivity, resulting in overestimates of dispersal ability and turnovers (extinctions plus colonizations). We extend estimation methods for the IFM by using a hierarchical Bayesian model to account both for false absences due to imperfect detection and for missing data due to sites not surveyed in some years. We compare parameter estimates, measures of metapopulation dynamics, and forecasts using stochastic patch occupancy models (SPOMs) among three IFM models: (1) a Bayesian formulation assuming no false absences and omitting site-year combinations with missing data; (2) a hierarchical Bayesian formulation assuming no false absences but incorporating missing data; and (3) a hierarchical Bayesian formulation allowing for imperfect detection and incorporating missing data. We fit the models to multiyear data sets of occupancy for two bird species that differ in body size and presumed dispersal ability but inhabit the same network of sites: the small Black Rail (Laterallus jamaicensis) and the medium-sized Virginia Rail (Rallus limicola). Incorporating missing data affected colonization parameters and led to lower estimates of dispersal ability for the Black Rail. Detection rates were high for the Black Rail in most years but moderate for the Virginia Rail. Incorporating imperfect detection resulted in higher occupancy and lower turnover rates for both species, with largest effects for the Virginia Rail. Forecasts using SPOMs were sensitive to both missing data and false absences; persistence in models assuming no false absences was more optimistic than from robust models. Our results suggest that incorporating false absences and missing data into the IFM can improve (1) estimates of dispersal ability and the effect of connectivity on colonization, (2) the scaling of extinction risk with patch area, and (3) forecasts of occupancy and turnover rates.     

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.     

Metacommunity dynamics over 16 years in a pyrogenic shrubland

Metacommunity theory allows predictions about the dynamics of potentially interacting species' assemblages that are linked by dispersal, but strong empirical tests of the theory are rare. We analyzed the metacommunity dynamics of Florida rosemary scrub, a patchily distributed pyrogenic community, to test predictions about turnover rates, community nestedness, and responses to patch size, arrangement, and quality. We collected occurrence data for 45 plant species from 88 rosemary scrub patches in 1989 and 2005 and used growth form, mechanism of regeneration after fire, and degree of habitat specialization to categorize species by life history. We tested whether patch size, fire history, and structural connectivity (a measure of proximity and size of surrounding patches) could be used to predict apparent extinctions and colonizations. In addition, we tested the accuracy of incidence-function models built with the patch survey data from 1989. After fire local extinction rates were higher for herbs than woody plants, higher for species that regenerated only from seed than species able to resprout, and higher for generalist than specialist species. Fewer rosemary specialists and a higher proportion of habitat generalists were extirpated on recently burned patches than on patches not burned between 1989 and 2005. Nestedness was highest for specialists among all life-history groups. Estimated model parameters from 1989 predicted the observed (1989-2005) extinction rates and the number of patches with persistent populations of individual species. These results indicate that species with different life-history strategies within the same metacommunity can have substantially different responses to patch configuration and quality. Real metacommunities may not conform to certain assumptions of simple models, but incidence-function models that consider only patch size, configuration, and quality can have significant predictive accuracy.     

The role of habitat disturbance and recovery in metapopulation persistence

Classical metapopulation theory assumes a static landscape. However, empirical evidence indicates many metapopulations are driven by habitat succession and disturbance. We develop a stochastic metapopulation model, incorporating habitat disturbance and recovery, coupled with patch colonization and extinction, to investigate the effect of habitat dynamics on persistence. We discover that habitat dynamics play a fundamental role in metapopulation dynamics. The mean number of suitable habitat patches is not adequate for characterizing the dynamics of the metapopulation. For a fixed mean number of suitable patches, we discover that the details of how disturbance affects patches and how patches recover influences metapopulation dynamics in a fundamental way. Moreover, metapopulation persistence is dependent not only on the average lifetime of a patch, but also on the variance in patch lifetime and the synchrony in patch dynamics that results from disturbance. Finally, there is an interaction between the habitat and metapopulation dynamics, for instance declining metapopulations react differently to habitat dynamics than expanding metapopulations. We close, emphasizing the importance of using performance measures appropriate to stochastic systems when evaluating their behavior, such as the probability distribution of the state of the metapopulation, conditional on it being extant (i.e., the quasistationary distribution).     

Dispersal depression with habitat fragmentation in the bog fritillary butterfly

Habitat fragmentation is expected to impose strong selective pressures on dispersal rates. However, evolutionary responses of dispersal are not self-evident, since various selection pressures act in opposite directions. Here we disentangled the components of dispersal behavior in a metapopulation context using the Virtual Migration model, and we linked their variation to habitat fragmentation in the specialist butterfly Proclossiana eunomia. Our study provided a nearly unique opportunity to study how habitat fragmentation modifies dispersal at the landscape scale, as opposed to microlandscapes or simulation studies. Indeed, we studied the same species in four landscapes with various habitat fragmentation levels, in which large amounts of field data were collected and analyzed using similar methodologies. We showed the existence of quantitative variations in dispersal behavior correlated with increased fragmentation. Dispersal propensity from habitat patches (for a given patch size), and mortality during dispersal (for a given patch connectivity) were lower in more fragmented landscapes. We suggest that these were the consequences of two different evolutionary responses of dispersal behavior at the individual level: (1) when fragmentation increased, the reluctance of individuals to cross habitat patch boundaries also increased; (2) when individuals dispersed, they flew straighter in the matrix, which is the best strategy to improve dispersal success. Such evolutionary responses could generate complex nonlinear patterns of dispersal changes at the metapopulation level according to habitat fragmentation. Due to the small size and increased isolation of habitat patches in fragmented landscapes, overall emigration rate and mortality during dispersal remained high. As a consequence, successful dispersal at the metapopulation scale remained limited. Therefore, to what extent the selection of individuals with a lower dispersal propensity and a higher survival during dispersal is able to limit detrimental effects of habitat fragmentation on dispersal success is unknown, and any conclusion that metapopulations would compensate for them is flawed.     

Effect of Rotational Shepherding on Demographic and Genetic Connectivity of Calcareous Grassland Plants

Response to habitat fragmentation may not be generalized among species, in particular for plant communities with a variety of dispersal traits. Calcareous grasslands are one of the most species‐rich habitats in Central Europe, but abandonment of traditional management has caused a dramatic decline of calcareous grassland species. In the Southern Franconian Alb in Germany, reintroduction of rotational shepherding in previously abandoned grasslands has restored species diversity, and it has been suggested that sheep support seed dispersal among grasslands. We tested the effect of rotational shepherding on demographic and genetic connectivity of calcareous grassland specialist plants and whether the response of plant populations to shepherding was limited to species dispersed by animals (zoochory). Specifically, we tested competing dispersal models and source and focal patch properties to explain landscape connectivity with patch‐occupancy data of 31 species. We fitted the same connectivity models to patch occupancy and nuclear microsatellite data for the herb Dianthus carthusianorum (Carthusian pink). For 27 species, patch connectivity was explained by dispersal by rotational shepherding regardless of adaptations to zoochory, whereas population size (16% species) and patch area (0% species) of source patches were not important predictors of patch occupancy in most species. [Correction made after online publication, February 25, 2014: Population size and patch area percentages were mistakenly inverted, and have now been fixed.] Microsite diversity of focal patches significantly increased the model variance explained by patch occupancy in 90% of the species. For D. carthusianorum, patch connectivity through rotational shepherding explained both patch occupancy and population genetic diversity. Our results suggest shepherding provides dispersal for multiple plant species regardless of their dispersal adaptations and thus offers a useful approach to restore plant diversity in fragmented calcareous grasslands. Efectos del Pastoreo Rotacional sobre la Conectividad Genética y Demográfica de Plantas de Pastizales Calcáreos     

Connectivity or demography: Defining sources and sinks in coral reef fish metapopulations

The identity of an individual patch as a source or a sink within a metapopulation is a function of its ability to produce individuals and to disperse them to other patches. In marine systems patch identity is very often defined by dispersal ability alone—upstream patches are sources—while issues of variable habitat quality (which affects local production) are ignored. This can have important ramifications for the science of marine reserve siting. This study develops a spatially explicit source–sink metapopulation model for reef fish and uses it to evaluate the relative importance of connectivity versus demography and how this depends upon the level of local larval retention and the strength of density-dependent recruitment. Elasticity analyses indicated that patch contribution (source or sink) was more sensitive to demographic parameters (particularly survival) than connectivity and this effect was conserved even under strong levels of density-dependence and was generally strengthened as local retention increased. Variability in the relationship between parameter elasticity and local retention was shown to be dependent upon the magnitude of connectivity for an individual patch relative to a critical connectivity value. The proportion of larvae lost due to transport processes was an important parameter which directly affected the magnitude of this critical connectivity value. Patches with connectivity values less than the critical value contributed to the metapopulation largely via production (i.e., local demographics most important). As local retention increased, so did the importance of demographic parameters in these patches. Patches with connectivity values greater than the critical value contributed largely via dispersal of larvae and thus the importance of local demographics decreased as local retention increased.     

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.     

Patch Occupancy and Potential Metapopulation Dynamics of Three Forest Mammals in Fragmented Afromontane Forest in South Africa

Abstract: We investigated the persistence of three medium-sized (2–9 kg), rare forest mammals in the fragmented mist-belt Podocarpus forests of the midlands of KwaZulu-Natal Province, South Africa. We recorded patch occupancy of blue duiker (   Philantomba monticola ), tree hyrax (   Dendrohyrax arboreus ), and samango monkey ( Cercopithecus mitis labiatus ) in 199 forest patches. Their rarity is ascribed to the fragmentation and destruction of their forest habitat. Incidence functions, derived from presence and absence data, were formulated as generalized linear models, and environmental effects were included in the fitted logistic models. The small and mostly solitary hyrax and duiker persisted in smaller patches than the large and social monkey. Although this result follows expectations based on relative home-range sizes of each species, the incidence probability of the samango monkey was invariant with increasing isolation, whereas a gradual decrease with increasing isolation was observed for the hyrax and duiker. Group dynamics may inhibit dispersal and increase the isolation effect in social species such as samango monkeys. A mainland-island metapopulation model adequately describes patterns of patch occupancy by the hyrax and duiker, but the monkeys' poor dispersal ability and obvious area-dependent extirpation suggest that they exist in transient, nonequilibrium (declining) metapopulations. Through identification of large forest patches for careful protection and management, the survival of all three species—especially the monkey—could be prolonged. Because no functional metapopulation may exist for the monkey, however, this is an emergency measure. For the duiker and hyrax, larger patches should form part of a network of smaller and closer patches in a natural matrix.     

Metapopulation persistence in a dynamic landscape: more habitat or better stewardship?

Habitat loss and fragmentation has created metapopulations where there were once continuous populations. Ecologists and conservation biologists have become interested in the optimal way to manage and conserve such metapopulations. Several authors have considered the effect of patch disturbance and recovery on metapopulation persistence, but almost all such studies assume that every patch is equally susceptible to disturbance. We investigated the influence of protecting patches from disturbance on metapopulation persistence, and used a stochastic metapopulation model to answer the question: How can we optimally trade off returns from protection of patches vs. creation of patches? We considered the problem of finding, under budgetary constraints, the optimal combination of increasing the number of patches in the metapopulation network vs. increasing the number of protected patches in the network. We discovered that the optimal trade-off is dependent upon all of the properties of the system: the species dynamics, the dynamics of the landscape, and the relative costs of each action. A stochastic model and accompanying methodology are provided allowing a manager to determine the optimal policy for small metapopulations. We also provide two approximations, including a rule of thumb, for determining the optimal policy for larger metapopulations. The method is illustrated with an example inspired by information for the greater bilby, Macrotis lagotis, inhabiting southwestern Queensland, Australia. We found that given realistic costs for each action, protection of patches should be prioritized over patch creation for improving the persistence of the greater bilby during the next 20 years.     

Developing Metapopulation Connectivity Criteria from Genetic and Habitat Data to Recover the Endangered Mexican Wolf

Restoring connectivity between fragmented populations is an important tool for alleviating genetic threats to endangered species. Yet recovery plans typically lack quantitative criteria for ensuring such population connectivity. We demonstrate how models that integrate habitat, genetic, and demographic data can be used to develop connectivity criteria for the endangered Mexican wolf (Canis lupus baileyi), which is currently being restored to the wild from a captive population descended from 7 founders. We used population viability analysis that incorporated pedigree data to evaluate the relation between connectivity and persistence for a restored Mexican wolf metapopulation of 3 populations of equal size. Decreasing dispersal rates greatly increased extinction risk for small populations (<150–200), especially as dispersal rates dropped below 0.5 genetically effective migrants per generation. We compared observed migration rates in the Northern Rocky Mountains (NRM) wolf metapopulation to 2 habitat‐based effective distance metrics, least‐cost and resistance distance. We then used effective distance between potential primary core populations in a restored Mexican wolf metapopulation to evaluate potential dispersal rates. Although potential connectivity was lower in the Mexican wolf versus the NRM wolf metapopulation, a connectivity rate of >0.5 genetically effective migrants per generation may be achievable via natural dispersal under current landscape conditions. When sufficient data are available, these methods allow planners to move beyond general aspirational connectivity goals or rules of thumb to develop objective and measurable connectivity criteria that more effectively support species recovery. The shift from simple connectivity rules of thumb to species‐specific analyses parallels the previous shift from general minimum‐viable‐population thresholds to detailed viability modeling in endangered species recovery planning. Desarrollo de Criterios de Conectividad Metapoblacional a Partir de Datos Genéticos y de Hábitat para Recuperar al Lobo Mexicano en Peligro de Extinción     

Modelling dispersal with diffusion and habitat selection: Analytical results for highly fragmented landscapes

Quantifying dispersal is crucial both for understanding ecological population dynamics, and for gaining insight into factors that affect the genetic structure of populations. The role of dispersal becomes pronounced in highly fragmented landscapes inhabited by spatially structured populations. We consider a landscape consisting of a set of habitat patches surrounded by unsuitable matrix, and model dispersal by assuming that the individuals follow a random walk with parameters that may be specific to the habitat type. We allow for spatial variation in patch quality, and account for edge-mediated behavior, the latter meaning that the individuals bias their movement towards the patches when close to an edge between a patch and the matrix. We employ a diffusion approximation of the random walk model to derive analytical expressions for various characteristics of the dispersal process. For example, we derive formulae for the time that an individual is expected to spend in its current patch i, and for the time that it will spend in the matrix, both conditional on the individual hitting next a given patch j before hitting any of the other patches or dying. The analytical formulae are based on the assumptions that the landscape is infinitely large, that the patches are circularly shaped, and that the patches are small compared to interpatch distances. We evaluate the effect of these assumptions by comparing the analytical results to numerical results in a real patch network that violates all of the three assumptions. We then consider a landscape that fulfills the assumptions, and show that in this case the analytical results are in a very good agreement with the numerical results. The results obtained here allow the construction of computationally efficient dispersal models that can be used as components of metapopulation models.     

Inferring Metapopulation Dynamics from Patch-Level Incidence of Florida Scrub Plants

Spatial structure and dynamics of multiple populations may explain species distribution patterns in patchy communities with heterogeneous disturbance regimes, especially when species have poor dispersal. The endemic-rich Florida (U.S.A.) rosemary scrub occupies about 4% of the west portion of Archbold Biological Station and occurs scattered within a matrix of less xeric vegetation. Longer fire-return times and higher frequency of open patches in rosemary scrub provide favorable habitat for many plant species. Occupancy of 123 species of vascular plants and ground lichens in 89 patches was determined by repeated site surveys. About two-thirds of the species occurring at more than 14 patches had a significant logistic regression of presence on time-since-fire, patch size, patch isolation, or their interactions. Species with presence related to the interaction between patch isolation and patch size were primarily herbs and small shrubs specializing in rosemary scrub. These results suggest the importance of spatial characteristics of the landscape for population turnover of these species. An incidence-based metapopulation model was used to predict extinction and colonization probabilities of those species with presence in rosemary scrub patches related to the studied spatial variables. This is the first attempt to apply incidence-based metapopulation models to plants. The results showed stronger effects of patch size and patch isolation on extinction probabilities of herbs than on those of woody species. Because of their effect on spatial heterogeneity and habitat availability, fire suppression and habitat destruction may decrease persistence probabilities for these rosemary scrub specialists, many of which are endangered species.     

Design and validation of a spatially explicit simulation model for bottomland hardwood forests

Species composition in forests depends on the interaction of species traits and species availability. Yet many forest simulation models focus only on interactions of adult trees and saplings, ignoring how species become members of the community. We modify a published forest model for bottomland hardwood forests (program SWAMP [Phipps, R.L., 1979. Simulation of wetlands forest vegetation dynamics. Ecol. Modell. 7, 257–288]) to make it spatially explicit and incorporate explicit seed production and dispersal algorithms. The resulting individual-based, spatially explicit forest simulator (YAFSIM) combines mechanistic seed dispersal with growth and mortality of trees to track forest dynamics over time. We describe the structure of the model and test its validity for dynamics in small bottomland hardwood patches in the Mississippi Alluvial Valley. Dynamics of species composition and basal areas of trees predicted by Yazoo Forest Simulator (YAFSIM) were similar to those of natural second- and old-growth bottomland forests. However, diversity of simulated forest patches declined over time largely because of random dynamics acting on small, isolated populations.     

Comparison of Two Types of Metapopulation Models in Real and Artificial Landscapes

Abstract: Application of metapopulation models is becoming increasingly widespread in the conservation of species in fragmented landscapes. We provide one of the first detailed comparisons of two of the most common modeling techniques, incidence function models and stage-based matrix models, and test their accuracy in predicting patch occupancy for a real metapopulation. We measured patch occupancies and demographic rates for regional populations of the Florida scrub lizard (   Sceloporus woodi ) and compared the observed occupancies with those predicted by each model. Both modeling strategies predicted patch occupancies with good accuracy ( 77–80%) and gave similar results when we compared hypothetical management scenarios involving removal of key habitat patches and degradation of habitat quality. To compare the two modeling approaches over a broader set of conditions, we simulated metapopulation dynamics for 150 artificial landscapes composed of equal-sized patches (2–1024 ha) spaced at equal distances (50–750 m). Differences in predicted patch occupancy were small to moderate (<20%) for about 74% of all simulations, but 22% of the landscapes had differences openface> 50%. Incidence function models and stage-based matrix models differ in their approaches, assumptions, and requirements for empirical data, and our findings provide evidence that the two models can produce different results. We encourage researchers to use both techniques and further examine potential differences in model output. The feasibility of obtaining data for population modeling varies widely among species and limits the modeling approaches appropriate for each species. Understanding different modeling approaches will become increasingly important as conservation programs undertake the challenge of managing for multiple species in a landscape context.     

Effects of matrix characteristics and interpatch distance on functional connectivity in fragmented temperate rainforests

The connectivity of remnant patches of habitat may affect the persistence of species in fragmented landscapes. We evaluated the effects of the structural connectivity of forest patches (i.e., distance between patches) and matrix class (land-cover type) on the functional connectivity of 3 bird species (the White-crested Elaenia [Elaenia albiceps], the Green-backed Firecrown Hummingbird [Sephanoides sephaniodes], and the Austral Thrush [Turdus falklandii]). We measured functional connectivity as the rate at which each species crossed from one patch to another. We also evaluated whether greater functional connectivity translated into greater ecological connectivity (dispersal of fruit and pollen) by comparing among forest patches fruit set of a plant pollinated by hummingbirds and abundance of seedlings and adults of 2 plants with bird- and wind-dispersed seeds. Interpatch distance was strongly associated with functional connectivity, but its effect was not independent of matrix class. For one of the bird-dispersed plants, greater functional connectivity for White-crested Elaenias and Austral Thrushes (both frugivores) was associated with higher densities of this plant. The lack of a similar association for the wind-dispersed species suggests this effect is linked to the dispersal vector. The abundance of the hummingbird-pollinated species was not related to the presence of hummingbirds. Interpatch distance and matrix class affect animal movement in fragmented landscapes and may have a cascading effect on the distribution of some animal-dispersed species. On the basis of our results, we believe effort should be invested in optimizing patch configuration and modifying the matrix so as to mitigate the effects of patch isolation in fragmented landscapes.     

Influence of habitat quality, population size, patch size, and connectivity on patch-occupancy dynamics of the middle spotted woodpecker

Despite extensive research on the effects of habitat fragmentation, the ecological mechanisms underlying colonization and extinction processes are poorly known, but knowledge of these mechanisms is essential to understanding the distribution and persistence of populations in fragmented habitats. We examined these mechanisms through multiseason occupancy models that elucidated patch-occupancy dynamics of Middle Spotted Woodpeckers (Dendrocopos medius) in northwestern Spain. The number of occupied patches was relatively stable from 2000 to 2010 (15-24% of 101 patches occupied every year) because extinction was balanced by recolonization. Larger and higher quality patches (i.e., higher density of oaks >37 cm dbh [diameter at breast height]) were more likely to be occupied. Habitat quality (i.e., density of large oaks) explained more variation in patch colonization and extinction than did patch size and connectivity, which were both weakly associated with probabilities of turnover. Patches of higher quality were more likely to be colonized than patches of lower quality. Populations in high-quality patches were less likely to become extinct. In addition, extinction in a patch was strongly associated with local population size but not with patch size, which means the latter may not be a good surrogate of population size in assessments of extinction probability. Our results suggest that habitat quality may be a primary driver of patch-occupancy dynamics and may increase the accuracy of models of population survival. We encourage comparisons of competing models that assess occupancy, colonization, and extinction probabilities in a single analytical framework (e.g., dynamic occupancy models) so as to shed light on the association of habitat quality and patch geometry with colonization and extinction processes in different settings and species.