Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Effect of habitat area and isolation on fragmented animal populations

Laura R. Prugh, Karen E. Hodges, Anthony R. E. Sinclair, and Justin S. Brashares
PNAS December 30, 2008 105 (52) 20770-20775; https://doi.org/10.1073/pnas.0806080105
Laura R. Prugh
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: prugh@nature.berkeley.edu
Karen E. Hodges
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anthony R. E. Sinclair
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Justin S. Brashares
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by Monica G. Turner, University of Wisconsin, Madison, WI, and approved October 29, 2008 (received for review June 25, 2008)

Related Articles

  • Importance of matrix habitats in maintaining biological diversity
    - Jan 07, 2009
  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

Habitat destruction has driven many once-contiguous animal populations into remnant patches of varying size and isolation. The underlying framework for the conservation of fragmented populations is founded on the principles of island biogeography, wherein the probability of species occurrence in habitat patches varies as a function of patch size and isolation. Despite decades of research, the general importance of patch area and isolation as predictors of species occupancy in fragmented terrestrial systems remains unknown because of a lack of quantitative synthesis. Here, we compile occupancy data from 1,015 bird, mammal, reptile, amphibian, and invertebrate population networks on 6 continents and show that patch area and isolation are surprisingly poor predictors of occupancy for most species. We examine factors such as improper scaling and biases in species representation as explanations and find that the type of land cover separating patches most strongly affects the sensitivity of species to patch area and isolation. Our results indicate that patch area and isolation are indeed important factors affecting the occupancy of many species, but properties of the intervening matrix should not be ignored. Improving matrix quality may lead to higher conservation returns than manipulating the size and configuration of remnant patches for many of the species that persist in the aftermath of habitat destruction.

  • incidence function
  • island biogeography
  • logistic regression
  • metaanalysis
  • occupancy

Habitat loss and fragmentation are major threats to terrestrial biodiversity (1). Globally, ≈40% of land has been converted for agricultural use (2), and regions as diverse as the eastern United States, the Philippines, and Ghana have lost >90% of their natural habitat (3, 4). Conservation theory and practice are founded on the principle that large habitat patches have more species than small ones and connected patches have more species than isolated ones (5). Although few would dispute this basic premise, we still do not know the general value of patch area and isolation as predictors of species occupancy in fragmented terrestrial systems. Despite hundreds of patch occupancy studies over >4 decades, there has been no quantitative synthesis of these findings. Several syntheses have examined species-area and diversity relationships (6, 7), but the species occupancy patterns that underlie diversity patterns in fragmented landscapes have been overlooked (8). How important is patch isolation relative to patch size in determining where species occur, and how consistent are these effects across diverse taxonomic groups? These are foundational, yet unanswered, questions for ecology and conservation biology.

We synthesized patch occupancy data from 89 studies of terrestrial fauna on 6 continents (Table S1) to determine how patch area and isolation affect species' occurrence patterns. Collectively, these studies recorded the occurrence of 785 animal species (Table 1) in 1,015 population networks surveyed in 12,370 discrete habitat patches. We use the term “population network” to refer to a spatially-structured population that occupies habitat patches embedded in a matrix of land cover deemed unsuitable for the species. Most of the population networks in our analysis were probably metapopulations linked by dispersal events, but few studies provided the necessary evidence for these linkages. Animals ranged in size from the 0.2-mm false spider mite (Pentamerismus sp.) to the 500-kg American bison (Bos bison), and habitat patches varied in isolation and area by 8 and 12 orders of magnitude, respectively (within-study order of magnitudes ranged from 0.2 to 4.4 for isolation and 0.8 to 5.4 for area). The dataset presents a global representation of faunal occurrence patterns in ecosystems made patchy both by humans (via agriculture, forestry, and urbanization; 72% of studies) and natural processes (28% of studies). We begin with the broad question: are patch area and isolation good predictors of occupancy for animals in terrestrial habitat fragments?

View this table:
  • View inline
  • View popup
Table 1.

Numbers of species included in the metaanalysis

Results and Discussion

Patch area and isolation were surprisingly poor predictors of occupancy across species. All together, 38% of patches were occupied (n = 55,855 occupancy records). We ran 4 logistic regression models for each population network to determine how well patch area and isolation predicted occupancy: (i) area only, (ii) isolation only, (iii) area + isolation, and (iv) area × isolation (see Methods). The amount of deviance (i.e., variation) in occupancy explained by each model (pseudoR2, or pR2) represents the ability of patch area and/or isolation to predict species occurrence patterns; this statistic is analogous to the R2 of linear regression (9). The pR2 values from each of the 1,015 population networks created distributions that were skewed toward 0, and few species were strongly influenced by patch area or isolation (Fig. 1A). The full area × isolation model accounted for a median of 25% of the deviance in occupancy, indicating that at least 75% of occupancy deviance was caused by other factors for most population networks. Area was a better predictor of occupancy than isolation: area alone accounted for a median of 13% of the deviance in occupancy versus 3% accounted for by isolation alone. These results raise 2 questions: first, why does area explain more deviance in occupancy than isolation, and second, why are both metrics such poor predictors of occupancy, given that habitat destruction is the dominant threat to species persistence?

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Strength of patch area and isolation effects on fragmented animal populations. Area explained more deviance in occupancy than isolation when all measures of isolation were included (A), whereas effects were equal when analyses were restricted to studies that used distance to nearest source population (NS) as the isolation measure (B). Isolation was measured as the distance to nearest patch, mainland, or source population in A. Four logistic regression models were run for each population network (when possible), with the following predictors of occupancy: (i) patch area, (ii) isolation, (iii) area + isolation, or (iv) area × isolation. The pR2 value (a goodness-of-fit measure analogous to R2 of linear regression) was recorded for each model. Box plots show medians (horizontal lines), interquartile ranges (boxes), the extent of nonoutlier datapoints (whiskers), and outliers (points).

The stronger effect of patch area can be explained by examining the type of isolation measure used in the analysis. Patch isolation is measured in many ways (10), and we used 3 categories to define the isolation measure used in each study: distance to nearest habitat patch of any size (NH); distance to nearest large patch, or “mainland” (NM); or distance to nearest occupied source patch (NS). NH and NM, which are landscape isolation measures, were used for 94% of the population networks, whereas NS (a demographic isolation measure) was used for only 6%. When we limited analyses to studies that used NS, the importance of isolation in predicting occupancy was equal to that of area (Fig. S1). Likewise, for the 82 population networks for which we could calculate both NS and NH, NS had a much stronger effect on occupancy (paired t test on slopes of logistic regression, P < 0.0001). Previous studies that found stronger effects of area on occupancy often used landscape isolation measures (7, 11) or “isolated/not isolated” categories for each patch (12). NS is more relevant to metapopulation dynamics because the proximity of source populations, not habitat patches per se, affects colonization probabilities.

Although demographic isolation was a better predictor of occupancy than landscape isolation, area and demographic isolation were still poor predictors: the full area × isolation model explained a median of only 24% of the deviance in occupancy across the subset of studies (Fig. 1B). NS may be superior to NND as an isolation measure, but a previous study (10) indicated that NS was inferior to more complex measures that take into account the areas and distances of all potential source patches, known as model Si. Data were available to calculate Si for 24 population networks (including the 2 metapopulations used in the previous study), and our larger analysis showed that Si did not explain any more variance in occupancy than NS did (see SI Text). Thus, NS appears to be as good a measure of isolation as the more complex measures that have been developed.

The overall poor performance of area and isolation as predictors of species occupancy could be caused by limitations of the data, ecological factors, or a combination of both. Snapshot occupancy data contain less information than population density or turnover data and may fail to capture the demographic advantages of large, highly connected habitat patches. Indeed, population networks that were monitored over multiple years had stronger relationships with area than those monitored for only 1 year (mean pR2 = 0.27 for multiyear records vs. 0.17 for single-year, F1,952 = 44.5, P < 0.0001). However, the relationship with isolation did not improve with additional years of monitoring (mean pR2 = 0.09 for both types of studies, F1,699 = 0.02, P = 0.90). Although this limitation likely explains some of the poor model fit, detailed examinations have suggested that occupancy data captures much of the information provided by population density data (13, 14). We propose 4 additional explanations for the low predictive power of patch area and isolation: (i) patches studied were of an inappropriate scale, such that areas and distances were not matched to focal species' body sizes and dispersal abilities; (ii) particular taxonomic groups or species with certain life history traits were less sensitive than others; (iii) most of the species were “survivors,” able to tolerate disturbance and not threatened with extinction; or (iv) the habitat island paradigm is not adequate in fragmented terrestrial systems because of strong effects of the matrix surrounding patches. We consider each of these.

Scale.

Assuming that researchers chose patch sizes that were appropriate to the body size of the study species, we would expect population networks from multispecies studies to be less sensitive to patch area, on average, than population networks from single-species studies. Most studies in this metaanalysis recorded the occupancy of multiple species in a single landscape, with body size ranges of up to 5 orders of magnitude, thus reducing the likelihood that patch sizes could be scaled appropriately for all species studied. Contrary to expectation, there was no relationship between the number of species in a study and the sensitivity of those species to patch area (n = 86 studies, R2 = 0.002, P = 0.68). In addition, mean patch area was roughly scaled to body size across population networks (r = 0.52, n = 971). The small effect of area on occupancy therefore does not appear to be caused by inappropriate scaling.

Likewise, isolation sensitivity was not strongly influenced by dispersal limitation. A dispersal limitation index was calculated as: Embedded Image where distmax is the maximum distance between patches in the landscape and dispmax is the maximum recorded dispersal distance for the species (n = 189). A species with a negative limitation index should not be dispersal limited. For example, a limitation index of −2 indicates that dispersal ability is 2 orders of magnitude greater than the maximum distance between patches. Limitation indices ranged from −3.2 for the red-tailed hawk (Buteo jamaicensis) in Mexico (41) to 2.9 for the northern pocket gopher (Thomomys talpoides) in Washington (42), and nearly half (49%) of the species were not dispersal-limited according to this index. Dispersal records are notoriously poor and often rely on chance observations, so maximum dispersal distances are likely underestimated for many species. Thus, the number of dispersal-limited species may be fewer than indicated by this index. If a lack of dispersal limitation explains the low influence of isolation on occupancy, then sensitivity to isolation should increase as dispersal limitation increases. There was a positive relationship between isolation sensitivity and dispersal limitation across species, indicating that the poor performance of isolation as a predictor of occupancy may be caused due by inappropriate scaling. However, the relationship between isolation sensitivity and dispersal limitation was weak, particularly for mammals and invertebrates (Fig. 2). Use of mean dispersal distances and mean or minimum patch distances to create the dispersal limitation index produced similar results. Although the scale of patch area and isolation may have been inappropriate for some species, scale issues are insufficient to explain the weak overall effect of area and isolation on species occurrence patterns.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Relationship between the dispersal limitation of a species and its sensitivity to patch isolation. The ability of patch isolation to predict occupancy was weakly related to dispersal limitation for all species combined (F1,192 = 4.04, R2 = 0.02, P = 0.05). See Scale in Results and Discussion for calculation of dispersal limitation. Relationships were stronger for amphibians and birds (A) than for invertebrates and mammals (B) but were not significant for any individual taxonomic group (amphibians: n = 7, R2 = 0.45, P = 0.1; birds: n = 31, R2 = 0.06, P = 0.2; mammals: n = 142, R2 = 0.02, P = 0.06; invertebrates: n = 13, R2 < 0.001, P = 0.98). Patch isolation included all 3 measures (nearest patch, mainland, or source). When restricted to nearest source, the relationship was similarly weak (F1,30 = 3.31, R2 = 0.10, P = 0.08). Alternative analyses using the slope parameter as a measure of effect size rather than pR2 showed no relationship between dispersal limitation and isolation sensitivity (F1,175 = 1.82, R2 = 0.01, P = 0.18 with all measures; F1,30 = 1.25, R2 = 0.04, P = 0.27 with nearest source only).

Species Traits.

Sensitivity to area and isolation may have been low if particularly sensitive taxonomic groups, or species with certain life history traits, were poorly represented in the database. We examined the effects of taxonomic group, diet, specialization, habit (terrestrial or arboreal), and fecundity on area and isolation sensitivity by using general linear models. None of these factors significantly affected isolation sensitivity, but area sensitivity was affected by taxonomic group, diet, and habit (Fig. 3). Birds and mammals were most sensitive to area whereas amphibians were least sensitive, carnivores (including insectivores) were more sensitive than omnivores, and arboreal species were more sensitive than terrestrial species. It is notable that specialists, which should be more restricted to habitat patches than generalists, did not have increased sensitivity to patch area or isolation, despite inclusion of several individual studies that showed higher area sensitivity of specialists (e.g., ref. 43). This finding contrasts the idea that the habitat island paradigm fits terrestrial systems best when species that use matrix habitats are excluded from analyses (15, 16). The groups with relatively high area sensitivity (birds, mammals, and carnivores) were those with relatively high metabolic requirements rather than those with habitat restrictions; arboreal species were the only area-sensitive group with obvious ties to habitat requirements. Difficulty in defining and delineating habitat may partially explain the low area sensitivity of amphibians, because different habitats are often used during different life stages. Thus, delineated patches may underestimate the habitat actually used. These traits explain some of the variation in area sensitivity across species, but they do not explain the overall low effect of area on occupancy: area-sensitive groups were not poorly represented in the database and in fact tended to be better represented (Fig. 3).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Influence of species traits on the strength of patch area effects. The ability of patch area to predict occupancy was affected by the taxonomic group, diet, and habit of the species (full model F7,930 = 14.2, R2 = 0.10, P < 0.0001; taxon F4 = 12.5, P < 0.0001, diet F2 = 5.2, P = 0.005, habit F1 = 7.7, P = 0.007). Diets were grouped such that “carnivore” included insectivores and parasitoids and “herbivore” included frugivores, nectivores, granivores, and detritivores. Habit was grouped such that “terrestrial” included fossorial and semiaquatic species. Least-squared means and SE bars are shown.

Species Are Survivors.

The landscapes in this metaanalysis had been fragmented for at least 30 years, so it is possible that area and isolation effects were weak because the most sensitive species had disappeared before surveys. Indeed, only 5% of the mammals, birds, and amphibians in the metaanalysis were classified as threatened by the International Union for Conservation of Nature (IUCN). Worldwide, 20% of species in these taxonomic groups are threatened (17). However, IUCN threat status was unrelated to area and isolation sensitivity (n = 641, R2 = 0.001, P = 0.39), indicating that low pR2 values were not caused by a lack of threatened species in the database. For example, the critically endangered mountain yellow-legged frog (Rana muscosa) and the endangered Amani sunbird (Anthreptes pallidigaster) were both insensitive to patch area (pR2 = 0.01 and 0.0003, respectively). Conversely, the 2 most area-sensitive species (the wood thrush, Hylocichla mustelina, and the veery, Catharus fuscescens) are ranked as “least concern” by the IUCN but are experiencing declines throughout their large ranges (ref. 18 and http://bna.birds.cornell.edu/BNA). Species in the database are “survivors” in that they have persisted through initial habitat destruction, but many of their populations are far from secure. Doomed populations can take far >30 years to disappear after fragmentation (19). This extinction debt could also lessen area effects, but we found that the reverse was true: population networks in systems with a natural matrix, which have generally been fragmented for thousands of years, were actually the least sensitive to area.

Patches Are Not Islands.

Unlike true oceanic islands, terrestrial habitat islands are not surrounded by a uniformly hostile matrix. Animals venturing outside of patches may find sufficiently benign conditions to live and reproduce, rendering the notion of the patch less relevant. We found strong effects of the matrix type on area and isolation sensitivity (Fig. 4). These broad patterns are particularly striking because the hostility of a matrix can vary widely among species within a landscape (20), and we therefore did not expect to find consistent effects of matrix type across species. Area sensitivity was higher in landscapes with human-dominated matrix types than in areas with natural matrix types. In contrast, isolation sensitivity was highest for population networks in forest patches surrounded by clear cuts. In these forestry-dominated landscapes, isolation effects were actually stronger than area effects even when landscape isolation measures were included.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Effect of predominant land cover in the matrix surrounding habitat patches on the sensitivity of species to patch area (●) and isolation (○). Sensitivity was measured as the proportion of deviance in occupancy accounted for by patch area or isolation in logistic regression analyses (pR2). Analyses were weighted such that each study (i.e., landscape) contributed equally to the results (n = 52 agricultural landscapes, 7 forestry clearcut, 5 urban, and 25 natural). Patch area was a better predictor of occupancy in landscapes made patchy by human activities (agriculture, forestry, urbanization) than by natural processes (F3, 950 = 5.2, P = 0.002). Patch isolation (all measures included) best predicted occupancy in forest patches surrounded by clear cuts (F3, 697 = 40.3, P < 0.0001). Results were the same when analyses were restricted to studies using a demographic isolation measure (F3, 77 = 11.7, P < 0.0001). Error bars show 95% confidence intervals.

The patch/nonpatch dichotomy appears to be a gross oversimplification for many species in fragmented landscapes. Irrespective of matrix effects, habitat patches often vary substantially in quality (21, 22). Many studies in this metaanalysis reported a strong influence of habitat quality on occupancy patterns. For example, patch area and isolation had little effect on the occupancy of most amphibians, whereas local habitat factors had stronger effects (45, 46). If habitat quality varies among patches, patch area may be a poor proxy of population size and thus unable to predict extinction probabilities accurately (23).

Habitat heterogeneity not only affects occupancy probabilities, but it also increases the difficulty of defining and delineating habitat patches in terrestrial landscapes (24). For example, extremely small patches may be ignored, and omission of small patches could lead to underestimates of area effects if those patches are empty. Furthermore, if some “nonhabitat” matrix areas are used by the animals, used patches may be both larger and closer together than those identified a priori by the researchers, thus contributing to the poor predictive power of area and isolation we detected.

Because habitat loss and fragmentation frequently co-occur, patch area and isolation may combine to reflect the influence of total habitat amount on a species (25). Previous work indicates that the total amount of habitat in a landscape may be the main determinant of population size and viability (26), and the spatial configuration of habitat does not influence occupancy until the proportion of habitat falls below a threshold (27–29). We were unable to test this hypothesis quantitatively because the total amount of habitat and nonhabitat in each landscape was unknown for most studies. However, maps of study areas in our metaanalysis generally showed that delineated habitat patches represented a small proportion of each landscape's total area. It is therefore unlikely that the poor predictive power of patch metrics across studies was caused by large amounts of identified habitat in each landscape.

Our empirical analysis of hundreds of animal species occurring in fragmented terrestrial landscapes throughout the world shows that patch area and isolation are poor predictors of occupancy in many cases. This finding does not mean patch metrics should be ignored altogether. Most species were more likely to be found in larger and less isolated patches, and some species were highly sensitive to patch area and isolation. However, the characteristics of intervening areas between patches, and of patches themselves, should not be ignored. The potential value of matrix areas for conservation has long been recognized (30–32), and our synthesis shows that landscape context has a strong effect on occupancy patterns across many taxonomic groups and ecosystems worldwide. The burgeoning body of research addressing matrix attributes should greatly increase our ability to predict and enhance species' persistence probabilities in fragmented systems (33–36). Identification of key features that determine whether a species will occupy an area should also improve both occupancy modeling and habitat restoration efforts. Minimizing the impact of threats that accompany habitat destruction, such as the spread of exotic species, overexploitation, and degradation of habitat within patches may have larger conservation returns than focusing on the amount and configuration of remaining habitat (37, 38), at least for those populations that persist in the aftermath of habitat destruction.

Methods

Data Acquisition.

In contrast to traditional metaanalyses, which are constrained by the often inconsistent statistical summaries reported in original publications, we maximized the size, accuracy, and standardization of our dataset by obtaining raw data directly from authors and conducting statistical analyses ourselves. Studies were found by a comprehensive search using the Web of Science in March 2005 with the terms “patch occupancy,” “habitat occupancy,” “metapopulation,” “island biogeography,” and “incidence function.” Forward and backward citations of articles found in these searches were also used to locate studies. Studies were included if 10 or more discrete terrestrial habitat patches (of the same general habitat type) were surveyed and raw occupancy data for birds, mammals, invertebrates, reptiles, or amphibians could be obtained. True island archipelagos were excluded. Hundreds of articles were screened, and 280 were examined in detail, 109 of which were found to be suitable for inclusion. Twenty of these studies were not included because raw data could not be obtained from the authors or their papers. The landscapes of the 89 studies included in the metaanalysis had been fragmented for a minimum of 30 years before surveys.

Logistic Regression Models.

Patch occupancy data were analyzed in the program R by using logistic regression models. In each model, presence/absence was used as the dependent variable for single-year studies (n = 67) and years present versus years surveyed was used for multiyear studies (n = 22). A binomial distribution was assumed, and data were generally not overdispersed or underdispersed (dispersion = 1.08 ± 0.51 SD). When patch area and isolation were both available, 4 models were run for each population network with the following predictors: (i) patch area only, (ii) patch isolation only, (iii) area plus isolation (additive), and (iv) area × isolation (interactive). Univariate models were run when only area or isolation was available. Patch area (hectares) and isolation (meters) were log-transformed before analyses to allow for direct comparison of area and isolation regression coefficients; these coefficients (i.e., slopes of the logistic regressions) indicate the change in occupancy probability given an order of magnitude change in area or isolation. Analyses were also conducted on untransformed data, and results were similar. Population networks were treated independently; adding “study” as a random effect did not change results.

In addition to the slopes, intercepts, and associated errors, we also calculated the pR2 of each model as: Embedded Image Null and residual deviances are analogous to the sum-of-squares estimates of linear regression (39). Although pR2 does not have all of the statistical properties of the ordinary least-squares R2 (such as a straightforward connection to the F statistic), it provides analogous goodness-of-fit estimates (9).

Paired t Tests on Slopes.

To determine the relative strength of area and isolation effects, we used paired t tests of the coefficients (i.e., slopes) of univariate logistic regression models by using area and isolation as predictors (n = 590 population networks with both variables). A paired t test of area and isolation slopes should not differ from 0 if the strengths of area and isolation effects on occupancy are equal, because area slopes should be positive (occupancy increases with patch area) and isolation slopes should be negative (occupancy decreases with patch isolation). Positive differences in slopes indicate stronger area effects.

Weighting.

In traditional metaanalyses, effect sizes are weighted by the inverse of the associated error estimates, so that studies with more precise estimates are given more weight in analyses (40). We did not weight analyses presented here, but see Figs. S2–S4 and SI Text for alternative analyses with weighted effect sizes. Weighting was not possible when comparing area and isolation slopes because there were 2 estimates of error (1 for each slope) and 1 estimate of effect size (the difference between the slopes). We did not weight pR2 values shown in Figs. 1–3 to facilitate ease of interpretation and display the distribution of pR2 values. Results from weighted and unweighted analyses were very similar.

Species and Landscape Traits.

The effect of species traits on sensitivity to patch area and isolation was examined by using general linear models. Because body size, maximum lifespan, and fecundity were highly correlated (r = −0.65 for body size and fecundity, r = 0.76 for body size and lifespan), only fecundity was included in models because it had the strongest relationship with area and isolation sensitivity. The starting model included the following predictors: taxonomic group, diet, fecundity, specialization, and habit (see Table S2 for a full list of species and traits). Habit was categorized as arboreal (species primarily living in trees) and terrestrial (species primarily living on the ground, below ground, or in low shrubs). Diet was categorized as carnivorous (including insectivores and parasitoids), herbivorous (including frugivores, nectivores, granivores, and detritivores), or omnivorous (consuming plant and animal material). Reclassifying detritivores as omnivores did not affect results. Species were classified as specialists (primarily use 1 habitat or food item, e.g., monophagous insects) or generalists (use several habitats or foods). Sources for all species traits are available on request. The stepAIC modeling function in program R was used for analyses, in which forward and backward stepwise regression and Akaike Information Criterion (AIC) values were used to rank models. A Gaussian link was used. The model with the lowest AIC value was retained as the best model (Fig. 2 and Fig. S4).

The effect of matrix quality on area and isolation sensitivity was examined by an ANOVA comparing the pR2 values among population networks in landscapes with different matrix types. Matrix types were classified as natural (e.g., meadows, forests), urban, agricultural, or clearcut. Landscapes containing multiple matrix types were classified according to the most prevalent type. ANOVAs were weighted such that each landscape contributed equally to the analysis to avoid overrepresentation of landscapes with occupancy records for large numbers of species (16).

Acknowledgments

We thank all of the researchers cited in SI Text for raw data for this metaanalysis; A. Jacob and A. Pelletier for help building the database; L. Harmon for assistance with programming in R; and D. Srivastava, D. Schluter, J. Goheen, G. Stewart, and anonymous reviewers for helpful comments and advice. This work was supported by a grant from Environment Canada.

Footnotes

  • 1To whom correspondence should be addressed. E-mail: prugh{at}nature.berkeley.edu
  • Author contributions: L.R.P. and K.E.H. designed research; L.R.P. and K.E.H. performed research; L.R.P. analyzed data; and L.R.P., K.E.H., A.R.E.S., and J.S.B. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/cgi/content/full/0806080105/DCSupplemental.

  • © 2008 by The National Academy of Sciences of the USA
View Abstract

References

  1. 1.↵
    1. Baille J,
    2. Hilton-Taylor C,
    3. Stuart SN
    (2004) 2004 IUCN Red List of Threatened Species: A Global Assessment (International Union for Conservation of Nature, Cambridge, UK).
  2. 2.↵
    1. Foley JA,
    2. et al.
    (2005) Global consequences of land use. Science 309:570–574.
    OpenUrlAbstract/FREE Full Text
  3. 3.↵
    1. Noss RF,
    2. Roe ETL,
    3. Scott JM
    (1995) Endangered Ecosystems in the United States: A Preliminary Assessment of Loss and Degradation (U.S. Department of the Interior, National Biological Service, Washington, DC).
  4. 4.↵
    1. World Resources Institute
    (1998) A Guide to World resources 1998–1999: Environmental Change and Human Health (World Resources Institute, Washington DC).
  5. 5.↵
    1. MacArthur RH,
    2. Wilson EO
    (1967) The Theory of Island Biogeography (Princeton Univ Press, Princeton).
  6. 6.↵
    1. Drakare S,
    2. Lennon JJ,
    3. Hillebrand H
    (2006) The imprint of the geographical, evolutionary and ecological context on species-area relationships. Ecol Lett 9:215–227.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Watling JI,
    2. Donnelly MA
    (2006) Fragments as islands: A synthesis of faunal responses to habitat patchiness. Conserv Biol 20:1016–1025.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Ovaskainen O,
    2. Hanski I
    (2003) The species-area relationship derived from species-specific incidence functions. Ecol Lett 6:903–909.
    OpenUrlCrossRef
  9. 9.↵
    1. Hagle TM,
    2. Mitchell GE
    (1992) Goodness-of-fit measures for probit and logit. Am J Pol Sci 36:762–784.
    OpenUrlCrossRef
  10. 10.↵
    1. Moilanen A,
    2. Nieminen M
    (2002) Simple connectivity measures in spatial ecology. Ecology 83:1131–1145.
    OpenUrlCrossRef
  11. 11.↵
    1. Pavlacky DC,
    2. Anderson SH
    (2007) Does avian species richness in natural patch mosaics follow the forest fragmentation paradigm? Anim Conserv 10:57–68.
    OpenUrl
  12. 12.↵
    1. Ferraz G,
    2. et al.
    (2007) A large-scale deforestation experiment: Effects of patch area and isolation on Amazon birds. Science 315:238–241.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Zuckerberg B,
    2. Porter WF,
    3. Corwin K
    (2008) The consistency and stability of abundance-occupancy relationships in large-scale population dynamics. J Anim Ecol, in press.
  14. 14.↵
    1. He FL,
    2. Gaston KJ
    (2000) Estimating species abundance from occurrence. Am Nat 156:553–559.
    OpenUrlCrossRef
  15. 15.↵
    1. Cook WM,
    2. Lane KT,
    3. Foster BL,
    4. Holt RD
    (2002) Island theory, matrix effects, and species richness patterns in habitat fragments. Ecol Lett 5:619–623.
    OpenUrlCrossRef
  16. 16.↵
    1. Bender DJ,
    2. Contreras TA,
    3. Fahrig L
    (1998) Habitat loss and population decline: A meta-analysis of the patch size effect. Ecology 79:517–533.
    OpenUrlCrossRef
  17. 17.↵
    1. IUCN Species Survival Commission
    (2007) 2007 IUCN Red List of Threatened Species (International Union for Conservation of Nature, Cambridge, UK).
  18. 18.↵
    1. Poole A
    (2005) Birds of North America (Cornell Laboratory of Ornithology, Ithaca, NY).
  19. 19.↵
    1. Vellend M,
    2. et al.
    (2006) Extinction debt of forest plants persists for more than a century following habitat fragmentation. Ecology 87:542–548.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Ricketts TH
    (2001) The matrix matters: Effective isolation in fragmented landscapes. Am Nat 158:87–99.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Fleishman E,
    2. Ray C,
    3. Sjogren-Gulve P,
    4. Boggs CL,
    5. Murphy DD
    (2002) Assessing the roles of patch quality, area, and isolation in predicting metapopulation dynamics. Conserv Biol 16:706–716.
    OpenUrlCrossRef
  22. 22.↵
    1. Thomas JA,
    2. et al.
    (2001) The quality and isolation of habitat patches both determine where butterflies persist in fragmented landscapes. Proc R Soc London Ser B 268:1791–1796.
    OpenUrlPubMed
  23. 23.↵
    1. Pellet J,
    2. Fleishman E,
    3. Dobkin DS,
    4. Gander A,
    5. Murphy DD
    (2007) An empirical evaluation of the area and isolation paradigm of metapopulation dynamics. Biol Conserv 136:483–495.
    OpenUrlCrossRef
  24. 24.↵
    1. Turner MG
    (2005) Landscape ecology: What is the state of the science? Annu Rev Ecol Evol System 36:319–344.
    OpenUrl
  25. 25.↵
    1. With KA,
    2. King AW
    (1997) The use and misuse of neutral landscape models in ecology. Oikos 79:219–229.
    OpenUrlCrossRef
  26. 26.↵
    1. Fahrig L
    (2003) Effects of habitat fragmentation on biodiversity. Annu Rev Ecol Evol System 34:487–515.
    OpenUrl
  27. 27.↵
    1. Flather CH,
    2. Bevers M
    (2002) Patchy reaction-diffusion and population abundance: The relative importance of habitat amount and arrangement. Am Nat 159:40–56.
    OpenUrlPubMed
  28. 28.↵
    1. Rhodes JR,
    2. et al.
    (2008) Regional variation in habitat-occupancy thresholds: A warning for conservation planning. J Appl Ecol 45:549–557.
    OpenUrl
  29. 29.↵
    1. Andren H
    (1994) Effects of habitat fragmentation on birds and mammals in landscapes with different proportions of suitable habitat: A review. Oikos 71:355–366.
    OpenUrlCrossRef
  30. 30.↵
    1. Franklin JF
    (1993) Preserving biodiversity: Species, ecosystems, or landscapes. Ecol Appl 3:202–205.
    OpenUrlCrossRef
  31. 31.↵
    1. Fahrig L,
    2. Merriam G
    (1994) Conservation of fragmented populations. Conserv Biol 8:50–59.
    OpenUrlPubMed
  32. 32.↵
    1. Fahrig L
    (2001) How much habitat is enough? Biol Conserv 100:65–74.
    OpenUrlCrossRef
  33. 33.↵
    1. Rosenzweig ML
    (2001) Loss of speciation rate will impoverish future diversity. Proc Natl Acad Sci USA 98:5404–5410.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Daily GC,
    2. Ehrlich PR,
    3. Sanchez-Azofeifa GA
    (2001) Countryside biogeography: Use of human-dominated habitats by the avifauna of southern Costa Rica. Ecol Appl 11:1–13.
    OpenUrlCrossRef
  35. 35.↵
    1. Baum KA,
    2. Haynes KJ,
    3. Dillemuth FP,
    4. Cronin JT
    (2004) The matrix enhances the effectiveness of corridors and stepping stones. Ecology 85:2671–2676.
    OpenUrlCrossRef
  36. 36.↵
    1. Sekercioglu CH,
    2. Loarie SR,
    3. Oviedo Brenes F,
    4. Ehrlich PR,
    5. Daily GC
    (2007) Persistence of forest birds in the Costa Rican agricultural countryside. Conserv Biol 21:482–494.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Brashares JS,
    2. Arcese P,
    3. Sam MK
    (2001) Human demography and reserve size predict wildlife extinction in West Africa. Proc R Soc London Ser B 268:2473–2478.
    OpenUrlPubMed
  38. 38.↵
    1. Li YM,
    2. Wilcove DS
    (2005) Threats to vertebrate species in China and the United States. Bioscience 55:147–153.
    OpenUrl
  39. 39.↵
    1. Chambers JM,
    2. Hastie TJ
    (1992) Statistical Models in S (Wadsworth & Brooks, Pacific Grove, CA).
  40. 40.↵
    1. Scheiner S,
    2. Gurevitch J
    1. Gurevitch J,
    2. Hedges LV
    (2001) in Design and Analysis of Ecological Experiments, eds Scheiner S, Gurevitch J (Oxford Univ Press, New York), pp 347–369.
  41. 41.↵
    1. Martinez-Morales MA
    (2005) Nested species assemblages as a tool to detect sensitivity to forest fragmentation: the case of cloud forest birds. Oikos 108:634–642.
    OpenUrl
  42. 42.↵
    1. Lomolino MV,
    2. Perault DR
    (2001) Island biogeography and landscape ecology of mammals inhabiting fragmented, temperate rain forests. Glob Ecol Biogeogr 10:113–132.
    OpenUrl
  43. 43.↵
    1. Steffan-Dewenter I,
    2. Tscharntke T
    (2000) Butterfly community structure in fragmented habitats. Ecol Lett 3:449–456.
    OpenUrlCrossRef
  44. 45.↵
    1. Knapp RA,
    2. Matthews KR,
    3. Preisler HK,
    4. Jellison R
    (2003) Developing probabilistic models to predict amphibian site occupancy in a patchy landscape. Ecol Appl 13:1069–1082.
    OpenUrl
  45. 46.↵
    1. Van Buskirk J
    (2005) Local and landscape influence on amphibian occurrence and abundance. Ecology 86:1936–1947.
    OpenUrl
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Effect of habitat area and isolation on fragmented animal populations
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Effect of habitat area and isolation on fragmented animal populations
Laura R. Prugh, Karen E. Hodges, Anthony R. E. Sinclair, Justin S. Brashares
Proceedings of the National Academy of Sciences Dec 2008, 105 (52) 20770-20775; DOI: 10.1073/pnas.0806080105

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Effect of habitat area and isolation on fragmented animal populations
Laura R. Prugh, Karen E. Hodges, Anthony R. E. Sinclair, Justin S. Brashares
Proceedings of the National Academy of Sciences Dec 2008, 105 (52) 20770-20775; DOI: 10.1073/pnas.0806080105
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 105 (52)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results and Discussion
    • Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Abstract depiction of a guitar and musical note
Science & Culture: At the nexus of music and medicine, some see disease treatments
Although the evidence is still limited, a growing body of research suggests music may have beneficial effects for diseases such as Parkinson’s.
Image credit: Shutterstock/agsandrew.
Scientist looking at an electronic tablet
Opinion: Standardizing gene product nomenclature—a call to action
Biomedical communities and journals need to standardize nomenclature of gene products to enhance accuracy in scientific and public communication.
Image credit: Shutterstock/greenbutterfly.
One red and one yellow modeled protein structures
Journal Club: Study reveals evolutionary origins of fold-switching protein
Shapeshifting designs could have wide-ranging pharmaceutical and biomedical applications in coming years.
Image credit: Acacia Dishman/Medical College of Wisconsin.
White and blue bird
Hazards of ozone pollution to birds
Amanda Rodewald, Ivan Rudik, and Catherine Kling talk about the hazards of ozone pollution to birds.
Listen
Past PodcastsSubscribe
Goats standing in a pin
Transplantation of sperm-producing stem cells
CRISPR-Cas9 gene editing can improve the effectiveness of spermatogonial stem cell transplantation in mice and livestock, a study finds.
Image credit: Jon M. Oatley.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Latest Articles
  • Archive

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490