Effects of Habitat Fragmentation on Effective Dispersal of Florida Scrub‐Jays

Abstract: Studies comparing dispersal in fragmented versus unfragmented landscapes show that habitat fragmentation alters the dispersal behavior of many species. We used two complementary approaches to explore Florida Scrub‐Jay (Aphelocoma c?rulescens) dispersal in relation to landscape fragmentation. First, we compared dispersal distances of color‐marked individuals in intensively monitored continuous and fragmented landscapes. Second, we estimated effective dispersal relative to the degree of fragmentation (as inferred from two landscape indexes: proportion of study site covered with Florida Scrub‐Jay habitat and mean distance to nearest habitat patch within each study site) by comparing genetic isolation‐by‐distance regressions among 13 study sites having a range of landscape structures. Among color‐banded individuals, dispersal distances were greater in fragmented versus continuous landscapes, a result consistent with other studies. Nevertheless, genetic analyses revealed that effective dispersal decreases as the proportion of habitat in the landscape decreases. These results suggest that although individual Florida Scrub‐Jays may disperse farther as fragmentation increases, those that do so are less successful as breeders than those that disperse short distances. Our study highlights the importance of combining observational data with genetic inferences when evaluating the complex biological and life‐history implications of dispersal.     

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.     

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.     

Corridor Length and Patch Colonization by a Butterfly, Junonia coenia

Abstract: Corridors have been proposed to reduce isolation and increase population persistence in fragmented landscapes, yet little research has evaluated the types of landscapes in which corridors will be most effective. I tested the hypothesis that corridors increase patch colonization by a butterfly, Junonia coenia , regardless of the butterfly's initial distance from a patch. I chose J. coenia because it has been shown to move between patches preferentially through corridors. Individuals were released 16–192 m away from open experimental patches into adjacent open corridors or forest. Neither corridors nor distance had a significant effect on patch colonization, but there was a significant interaction between the presence or absence of corridors and distance. At small distances (16–64 m), J. coenia was more likely to colonize open patches when released within forest than within open corridors, most likely because J. coenia used corridors as habitat. Nevertheless, patch colonization by butterflies released within forest decreased rapidly as distance from patches increased, as predicted by a null model of random movement. Colonization did not change with distance in the corridor, and at long distances (128–192 m), butterflies released in corridors were twice as likely to colonize open patches as those released in forest. These results suggest that one critical factor, interpatch distance, may determine the relative effectiveness of corridors and other landscape configurations, such as stepping stones, in reducing isolation in fragmented landscapes. When distances between patches are short compared to an animal's movement ability, a stepping stone approach may most effectively promote dispersal. Alternatively, the conservation value of corridors is highest relative to other habitat configurations when longer distances separate patches in fragmented landscapes.     

Habitat fragmentation affects movement behavior of migrating juvenile common toads

Displacement characteristics in animals are crucial drivers of successful movement decisions in resources acquisition, migration, and dispersal. As landscape structure is modified by human activity, mobility patterns are likely to evolve in response to habitat fragmentation. In species with complex life cycles that involve obligatory migrations between different habitats, one can predict that movement propensity would be promoted by fragmentation as long as it allows to reaching a habitat patch. Here, we compare the movement characteristics of naive toadlets sampled in populations distributed over a fragmentation gradient to test the hypothesis of a positive correlation between fragmentation and mobility levels. We studied toadlet movement in experimental arenas providing small patches of suitable conditions. We recorded the use of these patches (patch behavior) or the absence of their use (overtaking behavior). The more fragmented the original landscape, the more prone the toadlets were to combine these two behaviors, thus showing a higher motivation to explore. Moreover, the more fragmented the original landscape, the less the toadlets exhibited patch behavior. As the toadlets were reared in a common environment, the behavioral differences detected, relating to the level of fragmentation, resulted from inheritance. Our results thus illustrate that fragmentation is likely to create cross-generational transmittable variations in movement characteristics.     

Dispersal in Spatially Explicit Population Models

Abstract: Ruckelshaus et al. (1997) outlined a simulation model of dispersal between patches in a fragmented landscape. They showed that dispersal success—the proportion of dispersers successfully locating a patch—was particularly sensitive to errors in dispersal mortality and concluded that this limits the utility of spatially explicit population models in conservation biology. I contend that, although they explored error propagation in a simple dispersal model, they did not explore how errors are propagated in spatially explicit population models, as no consideration of population processes was included. I developed a simple simulation model to investigate the effect of varying dispersal success on predictions of patch occupancy and population viability, the conventional outputs of spatially explicit population models. The model simulates births and deaths within habitat patches and dispersal as the transfer of individuals between them. Model predictions were sensitive to changes in dispersal success across a restricted range of within-patch growth rates, which depended on the dispersal initiation mechanism, patch carrying capacities, and number of generations simulated. Predictions of persistence and patch occupancy were generally more sensitive to changes in dispersal success (1) under presaturation rather than saturation dispersal; (2) at lower patch carrying capacities; and (3) over longer time periods. The framework I present provides a means of assessing, quantitatively, the regions of parameter space for which differences in dispersal success are likely to have a large effect on population model outputs. Investigating the effect of the representation of dispersal behavior within the demographic and landscape context provides a more useful assessment of whether our lack of knowledge is likely to cause unacceptable uncertainty in the predictions of spatially explicit population models.     

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.     

Effectiveness of Corridors Relative to Enlargement of Habitat Patches

Abstract:  The establishment of biological corridors between two otherwise isolated habitat patches is a common yet contentious strategy for conserving populations in fragmented landscapes. We compared the effectiveness of corridors with the effectiveness of an alternate conservation strategy, the enlargement of existing habitat patches. We used a spatially explicit population model that simulated population size in two kinds of patches. One patch had a corridor that connected it to a larger "source" patch and the other patch was unconnected and enlarged at the periphery by an area the same size as the corridor. Patch isolation, corridor width, patch size, and the probability that individuals would cross the border from habitat to matrix were varied independently. In general, population size was greater in enlarged patches than in connected patches when patches were relatively large and isolated. Corridor width and the probability of crossing the border from habitat to matrix did not affect the relative benefit of corridors versus patch enlargement. Although biological corridors may mitigate potential effects of inbreeding depression at long time scales, our results suggest that they are not always the best method of conserving fragmented populations.     

Estimating Extinction Risk with Metapopulation Models of Large‐Scale Fragmentation

Habitat loss is the principal threat to species. How much habitat remains—and how quickly it is shrinking—are implicitly included in the way the International Union for Conservation of Nature determines a species’ risk of extinction. Many endangered species have habitats that are also fragmented to different extents. Thus, ideally, fragmentation should be quantified in a standard way in risk assessments. Although mapping fragmentation from satellite imagery is easy, efficient techniques for relating maps of remaining habitat to extinction risk are few. Purely spatial metrics from landscape ecology are hard to interpret and do not address extinction directly. Spatially explicit metapopulation models link fragmentation to extinction risk, but standard models work only at small scales. Counterintuitively, these models predict that a species in a large, contiguous habitat will fare worse than one in 2 tiny patches. This occurs because although the species in the large, contiguous habitat has a low probability of extinction, recolonization cannot occur if there are no other patches to provide colonists for a rescue effect. For 4 ecologically comparable bird species of the North Central American highland forests, we devised metapopulation models with area‐weighted self‐colonization terms; this reflected repopulation of a patch from a remnant of individuals that survived an adverse event. Use of this term gives extra weight to a patch in its own rescue effect. Species assigned least risk status were comparable in long‐term extinction risk with those ranked as threatened. This finding suggests that fragmentation has had a substantial negative effect on them that is not accounted for in their Red List category. Estimación del Riesgo de Extinción Mediante Modelos Metapoblacionales de Fragmentación a Gran Escala     

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.     

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.     

Landscape history, fragmentation, and patch occupancy: models for a forest bird with limited dispersal.

We developed and tested patch occupancy models for an endemic understory bird with limited dispersal ability, the Chucao Tapaculo (Scelorchilus rubecula), in two South American temperate rain forest landscapes that differed in levels and duration of forest loss. We assessed cover changes since 1961 in each landscape and surveyed patches for Chucao Tapaculo occupancy. We then developed incidence-based predictive models independently for each landscape and tested each model reciprocally in the alternative study area. We thereby assessed the domain of model applicability and identified those predictor variables with general effects and those that varied between the two landscapes. The two models were consistent regarding variable selection, and predictive accuracy of each model was high in the landscape where training data were collected. However, the models differed substantially in the magnitudes of effects related to patch size, with larger unoccupied patches observed in the landscape with the more advanced stage of fragmentation. Due to this discrepancy, each model performed poorly when applied to the alternative landscape, potentially reflecting the contrasting stages of habitat loss. Although it was impossible to dissociate effects of level and duration of forest loss, we viewed the landscapes as representing two extremes along a continuum of fragmentation, providing insights into potential trajectories for portions of the biome where deforestation is occurring. Further, our data suggest that static equilibrium models developed from distribution patterns in recently fragmented landscapes may overestimate persistence when used as a forecasting tool, or when extrapolated to alternative landscapes where fragmentation is more advanced. Thus, we suggest that landscapes used as standards for model building should be selected with caution. We recommend that distribution patterns be obtained from landscapes where fragmentation is advanced, preferably with histories of fragmentation long enough that time-delayed extinctions already would have occurred.     

Ranking individual habitat patches as connectivity providers: Integrating network analysis and patch removal experiments

Here we propose an integrated framework for modeling connectivity that can help ecologists, conservation planners and managers to identify patches that, more than others, contribute to uphold species dispersal and other ecological flows in a landscape context. We elaborate, extend and partly integrate recent network-based approaches for modeling and supporting the management of fragmented landscapes. In doing so, experimental patch removal techniques and network analytical approaches are merged into one integrated modeling framework for assessing the role of individual patches as connectivity providers. In particular, we focus the analyses on the habitat availability metrics PC and IIC and on the network metric Betweenness Centrality. The combination and extension of these metrics jointly assess both the immediate connectivity impacts of the loss of a particular patch and the resulting increased vulnerability of the network to subsequent disruptions. In using the framework to analyze the connectivity of two real landscapes in Madagascar and Catalonia (NE Spain), we suggest a procedure that can be used to rank individual habitat patches and show that the combined metrics reveal relevant and non-redundant information valuable to assert and quantify distinctive connectivity aspects of any given patch in the landscape. Hence, we argue that the proposed framework could facilitate more ecologically informed decision-making in managing fragmented landscapes. Finally, we discuss and highlight some of the advantages, limitations and key differences between the considered metrics.     

Genes and song: genetic and social connections in fragmented habitat in a woodland bird with limited dispersal

Understanding the processes leading to population declines in fragmented landscapes is essential for successful conservation management. However, isolating the influence of disparate processes, and dispersal in particular, is challenging. The Grey Shrike-thrush, Colluricincla harmonica, is a sedentary woodland-dependent songbird, with learned vocalizations whose incidence in suitable habitat patches falls disproportionally with decline in tree cover in the landscape. Although it has been suggested that gaps in tree cover might act as barriers to its dispersal, the species remains in many remnants of native vegetation in agricultural landscapes, suggesting that it may have responded to habitat removal and fragmentation by maintaining or even increasing dispersal distances. We quantified population connectivity of the Grey Shrike-thrush in a system fragmented over more than 120 years using genetic (microsatellites) and acoustic (song types) data. First, we tested for population genetic and acoustic structure at regional and local scales in search of barriers to dispersal or gene flow and signals of local spatial structuring indicative of restricted dispersal or localized acoustic similarity. Then we tested for effects of habitat loss and fragmentation on genetic and acoustic connectivity by fitting alternative models of mobility (isolation-by-distance [the null model] and reduced and increased movement models) across treeless vs. treed areas. Birds within -5 km of each other had more similar genotypes and song types than those farther away, suggesting that dispersal and song matching are limited in the region. Despite restricted dispersal detected for females (but not males), populations appeared to be connected by gene flow and displayed some cultural (acoustic) connectivity across the region. Fragmentation did not appear to impact greatly the dispersal of the Grey Shrike-thrush: none of the mobility models fit the genetic distances of males, whereas for females, an isolation-by-distance model could not be rejected in favor of the models of reduced or increased movement through treeless gaps. However, dissimilarities of the song types were more consistent with the model of reduced cultural connectivity through treeless areas, suggesting that fragmentation impedes song type sharing in the Grey Shrike-thrush. Our paper demonstrates that habitat fragmentation hinders important population processes in an Australian woodland bird even though its dispersal is not detectably impacted.     

The value of small size: loss of forest patches and ecological thresholds in southern Madagascar.

Many services generated by forest ecosystems provide essential support for human well-being. However, the vulnerability of these services to environmental change such as forest fragmentation are still poorly understood. We present spatial modeling of the generation of ecosystem services in a human-dominated landscape where forest habitat patches, protected by local taboos, are located in a matrix of cultivated land in southern Madagascar. Two ecosystem services dependent on the forest habitats were addressed: (1) crop pollination services by wild and semidomesticated bees (Apoidea), essential for local crop production of, for example, beans, and (2) seed dispersal services based on the presence of ring-tailed lemurs (Lemur catta). We studied the vulnerability of these ecosystem services to a plausible scenario of successive destruction of the smallest habitat patches. Our results indicate that, in spite of the fragmented nature of the landscape, the fraction of the landscape presently covered by both crop pollination and seed dispersal services is surprisingly high. It seems that the taboo system, though indirectly and unintentionally, contributes to upholding the generation of these services by protecting the forest patches. Both services are, however, predicted to be very vulnerable to the successive removal of small patches. For crop pollination, the rate of decrease in cover was significant even when only the smallest habitat patches were removed. The capacity for seed dispersal across the landscape displayed several thresholds with habitat patch removal. Our results suggest that, in order to maintain capacity for seed dispersal across the landscape and crop pollination cover in southern Androy, the geographical location of the remaining forest patches is more crucial than their size. We argue that in heavily fragmented production landscapes, small forest patches should increasingly be viewed as essential for maintaining ecosystem services, such as agricultural production, and also should be considered in the ongoing process of tripling the area of protected habitats in Madagascar.     

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).     

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.     

The Quantitative Incidence Function Model and Persistence of an Endangered Butterfly Metapopulation

The incidence function model is derived from a linear first-order Markov chain of the presence or absence of a species in a habitat patch. The model can be parameterized with "snapshot" presence/absence data from a patch network. Using the estimated parameter values the Markov chain can be iterated in the same or in some other patch network to generate quantitative predictions about transient metapopulation dynamics and the stochastic steady state. We tested the ability of the incidence function model to predict patch occupancy using extensive data on an endangered butterfly, the Glanville fritillary ( Melitaea cinxia ) Parameter values were estimated with data collected from a 50-patch network in 1991. In 1993 we surveyed the entire geographic range of the species in Finland, within an area of 50 × 70 km², with 1502 habitat patches (dry meadows) of which 536 were occupied. Model predictions were generated for the 1502 patches and were compared with the observed pattern of occupancy in 1993. The model predicted patch occupancy well in more than half of the study area, but prediction was poor for one quarter of the area, probably because of regional variation in habitat quality and because metapopulations may have been perturbed away from the steady state. The incidence function model provides a practical tool for making quantitative predictions about metapopulation dynamics of species living in fragmented landscapes.     

Effects of Patch Size and Type of Coffee Matrix on Ithomiine Butterfly Diversity and Dispersal in Cloud-Forest Fragments

Abstract: Determining the permeability of different types of landscape matrices to animal movement is essential for conserving populations in fragmented landscapes. We evaluated the effects of habitat patch size and matrix type on diversity, isolation, and dispersal of ithomiine butterflies in forest fragments surrounded by coffee agroecosystems in the Colombian Andes. Because ithomiines prefer a shaded understory, we expected the highest diversity and abundance in large fragments surrounded by shade coffee and the lowest in small fragments surrounded by sun coffee. We also thought shade coffee would favor butterfly dispersal and immigration into forest patches. We marked 9675 butterflies of 39 species in 12 forest patches over a year. Microclimate conditions were more similar to the forest interior in the shade‐coffee matrix than in the sun‐coffee matrix, but patch size and matrix type did not affect species richness and abundance in forest fragments. Furthermore, age structure and temporal recruitment patterns of the butterfly community were similar in all fragments, independent of patch size or matrix type. There were no differences in the numbers of butterflies flying in the matrices at two distances from the forest patch, but their behavior differed. Flight in the sun‐coffee matrix was rapid and directional, whereas butterflies in shade‐coffee matrix flew slowly. Seven out of 130 recaptured butterflies immigrated into patches in the shade‐coffee matrix, and one immigrated into a patch surrounded by sun coffee. Although the shade‐coffee matrix facilitated movement in the landscape, sun‐coffee matrix was not impermeable to butterflies. Ithomiines exhibited behavioral plasticity in habitat use and high mobility. These traits favor their persistence in heterogeneous landscapes, opening opportunities for their conservation. Understanding the dynamics and resource requirements of different organisms in rural landscapes is critical for identifying management options that address both animals’ and farmers’ needs.