Local Population Dynamics in Metapopulation Models: Implications for Conservation

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

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.     

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.     

Geometrical Factors and Metapopulation Dynamics of the Bush Cricket, Metrioptera bicolor Philippi (Orthoptera: Tettigoniidae)

This paper presents a metapopulation study of the bush cricket, Metrioptera bicolor , living in a recently fragmented landscape. The species inhabits grass and heathland patches of varying area and isolation. Analyses are made of how these geometrical factors affect local population size and density, distribution pattern, and the probability of local extinction and colonization. The proportion of available patches occupied varied between 72 and 79% during 1985–1990. Unoccupied patches were smaller and more isolated than those that were occupied. Patches where populations became extinct during this period were smaller than those with persisting populations. Since local population size was well correlated with patch area, it was clear that stochastic extinctions only occurred in small populations. Critical patch size for population extinction was approximately half a hectare. Colonized patches were less isolated than those that had not been colonized. Critical inter-patch distance for colonization was about 100 meters. The turnover was restricted to an identifiable share of the available patches. Only 33% of the patches were so small that extinction due to stochastic causes could be considered highly probable. This metapopulation will therefore most likely persist over a considerable period in its present spatial structure. There are apparent threats of further fragmentation, however, and nothing is known about the likelihood of large-scale extinctions resulting from extremely unfavorable weather conditions. Nevertheless, our results show that it is appropriate to include geometrical factors in metapopulation models.     

Spatial scale of local breeding habitat quality and adjustment of breeding decisions

Experimental studies provide evidence that, in spatially and temporally heterogeneous environments, individuals track variation in breeding habitat quality to adjust breeding decisions to local conditions. However, most experiments consider environmental variation at one spatial scale only, while the ability to detect the influence of a factor depends on the scale of analysis. We show that different breeding decisions by adults are based on information about habitat quality at different spatial scales. We manipulated (increased or decreased) local breeding habitat quality through food availability and parasite prevalence at a small (territory) and a large (patch) scale simultaneously in a wild population of Great Tits (Parus major). Females laid earlier in high-quality large-scale patches, but laying date did not depend on small-scale territory quality. Conversely, offspring sex ratio was higher (i.e., biased toward males) in high-quality, small-scale territories but did not depend on large-scale patch quality. Clutch size and territory occupancy probability did not depend on our experimental manipulation of habitat quality, but territories located at the edge of patches were more likely to be occupied than central territories. These results suggest that integrating different decisions taken by breeders according to environmental variation at different spatial scales is required to understand patterns of breeding strategy adjustment.     

Use of stochastic patch occupancy models in the California red-legged frog for Bayesian inference regarding past events and future persistence

Assessing causes of population decline is critically important to management of threatened species. Stochastic patch occupancy models (SPOMs) are popular tools for examining spatial and temporal dynamics of populations when presence–absence data in multiple habitat patches are available. We developed a Bayesian Markov chain method that extends existing SPOMs by focusing on past environmental changes that may have altered occupancy patterns prior to the beginning of data collection. Using occupancy data from 3 creeks, we applied the method to assess 2 hypothesized causes of population decline—in situ die-off and residual impact of past source population loss—in the California red-legged frog. Despite having no data for the 20–30 years between the hypothetical event leading to population decline and the first data collected, we were able to discriminate among hypotheses, finding evidence that in situ die-off increased in 2 of the creeks. Although the creeks had comparable numbers of occupied segments, owing to different extinction–colonization dynamics, our model predicted an 8-fold difference in persistence probabilities of their populations to 2030. Adding a source population led to a greater predicted persistence probability than did decreasing the in situ die-off, emphasizing that reversing the deleterious impacts of a disturbance may not be the most efficient management strategy. We expect our method will be useful for studying dynamics and evaluating management strategies of many species.     

Extinction debt of forest plants persists for more than a century following habitat fragmentation

Following habitat fragmentation individual habitat patches may lose species over time as they pay off their "extinction debt." Species with relatively low rates of population extinction and colonization ("slow" species) may maintain extinction debts for particularly prolonged periods, but few data are available to test this prediction. We analyzed two unusually detailed data sets on forest plant distributions and land-use history from Lincolnshire, United Kingdom, and Vlaams-Brabant, Belgium, to test for an extinction debt in relation to species-specific extinction and colonization rates. Logistic regression models predicting the presence-absence of 36 plant species were first parameterized using data from Lincolnshire, where forest cover has been relatively low (approximately 5-8%) for the past 1000 years. Consistent with extinction debt theory, for relatively slow species (but not fast species) these models systematically underpredicted levels of patch occupancy in Vlaams-Brabant, where forest cover was reduced from approximately 25% to <10% between 1775 and 1900 (it is presently 6.5%). As a consequence, the ability of the Lincolnshire models to predict patch occupancy in Vlaams-Brabant was worse for slow than for fast species. Thus, more than a century after forest fragmentation reached its current level an extinction debt persists for species with low rates of population turnover.     

A test of the substitution–habitat hypothesis in amphibians

Most examples that support the substitution‐habitat hypothesis (human‐made habitats act as substitutes of original habitat) deal with birds and mammals. We tested this hypothesis in 14 amphibians by using percentage occupancy as a proxy of habitat quality (i.e., higher occupancy percentages indicate higher quality). We classified water body types as original habitat (no or little human influence) depending on anatomical, behavioral, or physiological adaptations of each amphibian species. Ten species had relatively high probabilities (0.16–0.28) of occurrence in original habitat, moderate probability of occurrence in substitution habitats (0.11–0.14), and low probability of occurrence in refuge habitats (0.05–0.08). Thus, the substitution–habitat hypothesis only partially applies to amphibians because the low occupancy of refuges could be due to the negligible human persecution of this group (indicating good conservation status). However, low occupancy of refuges could also be due to low tolerance of refuge conditions, which could have led to selective extinction or colonization problems due to poor dispersal capabilities. That original habitats had the highest probabilities of occupancy suggests amphibians have a good conservation status in the region. They also appeared highly adaptable to anthropogenic substitution habitats.     

Minimum viable metapopulation size, extinction debt, and the conservation of a declining species.

A key question facing conservation biologists is whether declines in species' distributions are keeping pace with landscape change, or whether current distributions overestimate probabilities of future persistence. We use metapopulations of the marsh fritillary butterfly Euphydryas aurinia in the United Kingdom as a model system to test for extinction debt in a declining species. We derive parameters for a metapopulation model (incidence function model, IFM) using information from a 625-km2 landscape where habitat patch occupancy, colonization, and extinction rates for E. aurinia depend on patch connectivity, area, and quality. We then show that habitat networks in six extant metapopulations in 16-km2 squares were larger, had longer modeled persistence times (using IFM), and higher metapopulation capacity (lambdaM) than six extinct metapopulations. However, there was a > 99% chance that one or more of the six extant metapopulations would go extinct in 100 years in the absence of further habitat loss. For 11 out of 12 networks, minimum areas of habitat needed for 95% persistence of metapopulation simulations after 100 years ranged from 80 to 142 ha (approximately 5-9% of land area), depending on the spatial location of habitat. The area of habitat exceeded the estimated minimum viable metapopulation size (MVM) in only two of the six extant metapopulations, and even then by only 20%. The remaining four extant networks were expected to suffer extinction in 15-126 years. MVM was consistently estimated as approximately 5% of land area based on a sensitivity analysis of IFM parameters and was reduced only marginally (to approximately 4%) by modeling the potential impact of long-distance colonization over wider landscapes. The results suggest a widespread extinction debt among extant metapopulations of a declining species, necessitating conservation management or reserve designation even in apparent strongholds. For threatened species, metapopulation modeling is a potential means to identify landscapes near to extinction thresholds, to which conservation measures can be targeted for the best chance of success.     

The importance of incorporating functional habitats into conservation planning for highly mobile species in dynamic systems

The distribution of mobile species in dynamic systems can vary greatly over time and space. Estimating their population size and geographic range can be problematic and affect the accuracy of conservation assessments. Scarce data on mobile species and the resources they need can also limit the type of analytical approaches available to derive such estimates. We quantified change in availability and use of key ecological resources required for breeding for a critically endangered nomadic habitat specialist, the Swift Parrot (Lathamus discolor). We compared estimates of occupied habitat derived from dynamic presence‐background (i.e., presence‐only data) climatic models with estimates derived from dynamic occupancy models that included a direct measure of food availability. We then compared estimates that incorporate fine‐resolution spatial data on the availability of key ecological resources (i.e., functional habitats) with more common approaches that focus on broader climatic suitability or vegetation cover (due to the absence of fine‐resolution data). The occupancy models produced significantly (P < 0.001) smaller (up to an order of magnitude) and more spatially discrete estimates of the total occupied area than climate‐based models. The spatial location and extent of the total area occupied with the occupancy models was highly variable between years (131 and 1498 km²). Estimates accounting for the area of functional habitats were significantly smaller (2–58% [SD 16]) than estimates based only on the total area occupied. An increase or decrease in the area of one functional habitat (foraging or nesting) did not necessarily correspond to an increase or decrease in the other. Thus, an increase in the extent of occupied area may not equate to improved habitat quality or function. We argue these patterns are typical for mobile resource specialists but often go unnoticed because of limited data over relevant spatial and temporal scales and lack of spatial data on the availability of key resources. Understanding changes in the relative availability of functional habitats is crucial to informing conservation planning and accurately assessing extinction risk for mobile resource specialists.     

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.     

Neighborhood and habitat effects on vital rates: expansion of the Barred Owl in the Oregon coast ranges

In this paper, we modify dynamic occupancy models developed for detection-nondetection data to allow for the dependence of local vital rates on neighborhood occupancy, where neighborhood is defined very flexibly. Such dependence of occupancy dynamics on the status of a relevant neighborhood is pervasive, yet frequently ignored. Our framework permits joint inference about the importance of neighborhood effects and habitat covariates in determining colonization and extinction rates. Our specific motivation is the recent expansion of the Barred Owl (Strix varia) in western Oregon, USA, over the period 1990-2010. Because the focal period was one of dramatic range expansion and local population increase, the use of models that incorporate regional occupancy (sources of colonists) as determinants of dynamic rate parameters is especially appropriate. We began our analysis of 21 years of Barred Owl presence/nondetection data in the Tyee Density Study Area (TDSA) by testing a suite of six models that varied only in the covariates included in the modeling of detection probability. We then tested whether models that used regional occupancy as a covariate for colonization and extinction outperformed models with constant or year-specific colonization or extinction rates. Finally we tested whether habitat covariates improved the AIC of our models, focusing on which habitat covariates performed best, and whether the signs of habitat effects are consistent with a priori hypotheses. We conclude that all covariates used to model detection probability lead to improved AIC, that regional occupancy influences colonization and extinction rates, and that habitat plays an important role in determining extinction and colonization rates. As occupancy increases from low levels toward equilibrium, colonization increases and extinction decreases, presumably because there are more and more dispersing juveniles. While both rates are affected, colonization increases more than extinction decreases. Colonization is higher and extinction is lower in survey polygons with more riparian forest. The effects of riparian forest on extinction rates are greater than on colonization rates. Model results have implications for management of the invading Barred Owl, both through habitat alteration and removal.     

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     

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.     

Estimates of minimum patch size depend on the method of estimation and the condition of the habitat

Minimum patch size for a viable population can be estimated in several ways. The density-area method estimates minimum patch size as the smallest area in which no new individuals are encountered as one extends the arbitrary boundaries of a study area outward. The density-area method eliminates the assumption of no variation in density with size of habitat area that accompanies other methods, but it is untested in situations in which habitat loss has confined populations to small areas. We used a variant of the density area method to study the minimum patch size for the gopher tortoise (Gopherus polyphemus) in Florida, USA, where this keystone species is being confined to ever smaller habitat fragments. The variant was based on the premise that individuals within populations are likely to occur at unusually high densities when confined to small areas, and it estimated minimum patch size as the smallest area beyond which density plateaus. The data for our study came from detailed surveys of 38 populations of the tortoise. For all 38 populations, the areas occupied were determined empirically, and for 19 of them, duplicate surveys were undertaken about a decade apart. We found that a consistent inverse density area relationship was present over smaller areas. The minimum patch size estimated from the density-area relationship was at least 100 ha, which is substantially larger than previous estimates. The relative abundance of juveniles was inversely related to population density for sites with relatively poor habitat quality, indicating that the estimated minimum patch size could represent an extinction threshold. We concluded that a negative density area relationship may be an inevitable consequence of excessive habitat loss. We also concluded that any detrimental effects of an inverse density area relationship may be exacerbated by the deterioration in habitat quality that often accompanies habitat loss. Finally, we concluded that the value of any estimate of minimum patch size as a conservation tool is compromised by excessive habitat loss.     

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     

Role of Patch Size, Disease, and Movement in Rapid Extinction of Bighorn Sheep

Abstract: The controversy (  Berger 1990, 1999 ; Wehausen 1999 ) over rapid extinction in bighorn sheep ( Ovis canadensis ) has focused on population size alone as a correlate to persistence time. We report on the persistence and population performance of 24 translocated populations of bighorn sheep. Persistence in these sheep was strongly correlated with larger patch sizes, greater distance to domestic sheep, higher population growth rates, and migratory movements, as well as to larger population sizes. Persistence was also positively correlated with larger average home-range size ( p = 0.058, n = 10 translocated populations) and home-range size of rams ( p = 0.087, n = 8 translocated populations). Greater home-range size and dispersal rates of bighorn sheep were positively correlated to larger patches. We conclude that patch size and thus habitat carrying capacity, not population size per se, is the primary correlate to both population performance and persistence. Because habitat carrying capacity defines the upper limit to population size, clearly the amount of suitable habitat in a patch is ultimately linked to population size. Larger populations (250+ animals) were more likely to recover rapidly to their pre-epizootic survey number following an epizootic ( p = 0.019), although the proportion of the population dying in the epizootic also influenced the probability of recovery ( p = 0.001). Expensive management efforts to restore or increase bighorn sheep populations should focus on large habitat patches located ≥23 km from domestic sheep, and less effort should be expended on populations in isolated, small patches of habitat.     

Transient dynamics of invasive competition: barred owls, spotted owls, habitat, and the demons of competition present

The recent range expansion of Barred Owls (Strix varia) into the Pacific Northwest, where the species now co-occurs with the endemic Northern Spotted Owl (Strix occidentalis caurina), resulted in a unique opportunity to investigate potential competition between two congeneric, previously allopatric species. The primary criticism of early competition research was the use of current species' distribution patterns to infer past processes; however, the recent expansion of the Barred Owl and the ability to model the processes that result in site occupancy (i.e., colonization and extinction) allowed us to address the competitive process directly rather than inferring past processes through current patterns. The purpose of our study was to determine whether Barred Owls had any negative effects on occupancy dynamics of nesting territories by Northern Spotted Owls and how these effects were influenced by habitat characteristics of Spotted Owl territories. We used single-species, multi-season occupancy models and covariates quantifying Barred Owl detections and habitat characteristics to model extinction and colonization rates of Spotted Owl pairs in southern Oregon, USA. We observed a strong, negative association between Barred Owl detections and colonization rates and a strong positive effect of Barred Owl detections on extinction rates of Spotted Owls. We observed increased extinction rates in response to decreased amounts of old forest at the territory core and higher colonization rates when old-forest habitat was less fragmented. Annual site occupancy for pairs reflected the strong effects of Barred Owls on occupancy dynamics with much lower occupancy rates predicted for territories where Barred Owls were detected. The strong Barred Owl and habitat effects on occupancy dynamics of Spotted Owls provided evidence of interference competition between the species. These effects increase the importance of conserving large amounts of contiguous, old-forest habitat to maintain Northern Spotted Owls in the landscape.     

Effects of landscape transformation on bird colonization and extinction patterns in a large‐scale,long‐term natural experiment

Conversion of agricultural land to forest plantations is a major driver of global change. Studies on the impact of forest plantations on biodiversity in plantations and in the surrounding native vegetation have been inconclusive. Consequently, it is not known how to best manage the extensive areas of the planet currently covered by plantations. We used a novel, long‐term (16 years) and large‐scale (30,000 ha) landscape transformation natural experiment (the Nanangroe experiment, Australia) to test the effects of land conversion on population dynamics of 64 bird species associated with woodland and forest. A unique aspect of our study is that we focused on the effects of plantations on birds in habitat patches within plantations. Our study design included 56 treatment sites (Eucalyptus patches where the surrounding matrix was converted from grazed land to pine plantations), 55 control sites (Eucalyptus patches surrounded by grazed land), and 20 matrix sites (sites within the pine plantations and grazed land). Bird populations were studied through point counts, and colonization and extinction patterns were inferred through multiple season occupancy models. Large‐scale pine plantation establishment affected the colonization or extinction patterns of 89% of studied species and thus led to a comprehensive turnover in bird communities inhabiting Eucalyptus patches embedded within the maturing plantations. Smaller bodied species appeared to respond positively to plantations (i.e., colonization increased and extirpation of these species decreased in patches surrounded by plantations) because they were able to use the newly created surrounding matrix. We found that the effects of forest plantations affected the majority of the bird community, and we believe these effects could lead to the artificial selection of one group of species at the expense of another.     

Obstruction of biodiversity conservation by minimum patch size criteria

Minimum patch size criteria for habitat protection reflect the conservation principle that a single large (SL) patch of habitat has higher biodiversity than several small (SS) patches of the same total area (SL > SS). Nonetheless, this principle is often incorrect, and biodiversity conservation requires placing more emphasis on protection of large numbers of small patches (SS > SL). We used a global database reporting the abundances of species across hundreds of patches to assess the SL > SS principle in systems where small patches are much smaller than the typical minimum patch size criteria applied for biodiversity conservation (i.e., ∼85% of patches <100 ha). The 76 metacommunities we examined included 4401 species in 1190 patches. From each metacommunity, we resampled species–area accumulation curves to evaluate how biodiversity responded to habitat existing as a few large patches or as many small patches. Counter to the SL > SS principle and consistent with previous syntheses, species richness accumulated more rapidly when adding several small patches (45.2% SS > SL vs. 19.9% SL > SS) to reach the same cumulative area, even for the very small patches in our data set. Responses of taxa to habitat fragmentation differed, which suggests that when a given total area of habitat is to be protected, overall biodiversity conservation will be most effective if that habitat is composed of as many small patches as possible, plus a few large ones. Because minimum patch size criteria often require larger patches than the small patches we examined, our results suggest that such criteria hinder efforts to protect biodiversity.