Slopes of Avian Species-Area Relationships, Human Population Density, and Environmental Factors

There is increasing interest in how humans influence spatial patterns in biodiversity. One of the most frequently noted and marked of these patterns is the increase in species richness with area, the species-area relationship (SAR). SARs are used for a number of conservation purposes, including predicting extinction rates, setting conservation targets, and identifying biodiversity hotspots. Such applications can be improved by a detailed understanding of the factors promoting spatial variation in the slope of SARs, which is currently the subject of a vigorous debate. Moreover, very few studies have considered the anthropogenic influences on the slopes of SARs; this is particularly surprising given that in much of the world areas with high human population density are typically those with a high number of species, which generates conservation conflicts. Here we determine correlates of spatial variation in the slopes of species-area relationships, using the British avifauna as a case study. Whilst we focus on human population density, a widely used index of human activities, we also take into account (1) the rate of increase in habitat heterogeneity with increasing area, which is frequently proposed to drive SARs, (2) environmental energy availability, which may influence SARs by affecting species occupancy patterns, and (3) species richness. We consider environmental variables measured at both local (10 km x 10 km) and regional (290 km x 290 km) spatial grains, but find that the former consistently provides a better fit to the data. In our case study, the effect of species richness on the slope SARs appears to be scale dependent, being negative at local scales but positive at regional scales. In univariate tests, the slope of the SAR correlates negatively with human population density and environmental energy availability, and positively with the rate of increase in habitat heterogeneity. We conducted two sets of multiple regression analyses, with and without species richness as a predictor. When species richness is included it exerts a dominant effect, but when it is excluded temperature has the dominant effect on the slope of the SAR, and the effects of other predictors are marginal.

Fractal species distributions do not produce power-law species-area relationships

We derive the species-area relationship (SAR) expected from an assemblage of fractally distributed species. If species have truly fractal spatial distributions with different fractal dimensions, we show that the expected SAR is not the classical power-law function, as suggested recently in the literature. This analytically derived SAR has a distinctive shape that is not commonly observed in nature: upward-accelerating richness with increasing area (when plotted on log-log axes). This suggests that, in reality, most species depart from true fractal spatial structure. We demonstrate the fitting of a fractal SAR using two plant assemblages (Alaskan trees and British grasses). We show that in both cases, when modelled as fractal patterns, the modelled SAR departs from the observed SAR in the same way, in accord with the theory developed here. The challenge is to identify how species depart from fractality, either individually or within assemblages, and more importantly to suggest reasons why species distributions are not self-similar and what, if anything, this can tell us about the spatial processes involved in their generation.

Cichlid species-area relationships are shaped by adaptive radiations that scale with area

A positive relationship between species richness and island size is thought to emerge from an equilibrium between immigration and extinction rates, but the influence of species diversification on the form of this relationship is poorly understood. Here, we show that within-lake adaptive radiation strongly modifies the species-area relationship for African cichlid fishes. The total number of species derived from in situ speciation increases with lake size, resulting in faunas orders of magnitude higher in species richness than faunas assembled by immigration alone. Multivariate models provide evidence for added influence of lake depth on the species-area relationship. Diversity of clades representing within-lake radiations show responses to lake area, depth and energy consistent with limitation by these factors, suggesting that ecological factors influence the species richness of radiating clades within these ecosystems. Together, these processes produce lake fish faunas with highly variable composition, but with diversities that are well predicted by environmental variables.

Species-area and species-individual relationships for tropical trees: A comparison of three 50-ha plots

1 Species-accumulation curves for woody plants were calculated in three tropical forests, based on fully mapped 50-ha plots in wet, old-growth forest in Peninsular Malaysia, in moist, old-growth forest in central Panama, and in dry, previously logged forest in southern India. A total of 610 000 stems were identified to species and mapped to < Im accuracy. Mean species number and stem number were calculated in quadrats as small as 5 m x 5 m to as large as 1000 m x 500 m, for a variety of stem sizes above 10 mm in diameter. Species-area curves were generated by plotting species number as a function of quadrat size; species-individual curves were generated from the same data, but using stem number as the independent variable rather than area. 2 Species-area curves had different forms for stems of different diameters, but species-individual curves were nearly independent of diameter class. With < 10(4) stems, species-individual curves were concave downward on log-log plots, with curves from different forests diverging, but beyond about 104 stems, the log-log curves became nearly linear, with all three sites having a similar slope. This indicates an asymptotic difference in richness between forests: the Malaysian site had 2.7 times as many species as Panama, which in turn was 3.3 times as rich as India. 3 Other details of the species-accumulation relationship were remarkably similar between the three sites. Rectangular quadrats had 5-27% more species than square quadrats of the same area, with longer and narrower quadrats increasingly diverse. Random samples of stems drawn from the entire 50 ha had 10-30% more species than square quadrats with the same number of stems. At both Pasoh and BCI, but not Mudumalai. species richness was slightly higher among intermediate-sized stems (50-100mm in diameter) than in either smaller or larger sizes, These patterns reflect aggregated distributions of individual species, plus weak density-dependent forces that tend to smooth the species abundance distribution and 'loosen' aggregations as stems grow. 4 The results provide support for the view that within each tree community, many species have their abundance and distribution guided more by random drift than deterministic interactions. The drift model predicts that the species-accumulation curve will have a declining slope on a log-log plot, reaching a slope of O.1 in about 50 ha. No other model of community structure can make such a precise prediction. 5 The results demonstrate that diversity studies based on different stem diameters can be compared by sampling identical numbers of stems. Moreover, they indicate that stem counts < 1000 in tropical forests will underestimate the percentage difference in species richness between two diverse sites. Fortunately, standard diversity indices (Fisher's sc, Shannon-Wiener) captured diversity differences in small stem samples more effectively than raw species richness, but both were sample size dependent. Two nonparametric richness estimators (Chao. jackknife) performed poorly, greatly underestimating true species richness.

Slopes of Avian Species-Area Relationships, Human Population Density, and Environmental Factors

There is increasing interest in how humans influence spatial patterns in biodiversity. One of the most frequently noted and marked of these patterns is the increase in species richness with area, the species–area relationship (SAR). SARs are used for a number of conservation purposes, including predicting extinction rates, setting conservation targets, and identifying biodiversity hotspots. Such applications can be improved by a detailed understanding of the factors promoting spatial variation in the slope of SARs, which is currently the subject of a vigorous debate. Moreover, very few studies have considered the anthropogenic influences on the slopes of SARs; this is particularly surprising given that in much of the world areas with high human population density are typically those with a high number of species, which generates conservation conflicts. Here we determine correlates of spatial variation in the slopes of species–area relationships, using the British avifauna as a case study. Whilst we focus on human population density, a widely used index of human activities, we also take into account (1) the rate of increase in habitat heterogeneity with increasing area, which is frequently proposed to drive SARs, (2) environmental energy availability, which may influence SARs by affecting species occupancy patterns, and (3) species richness. We consider environmental variables measured at both local (10 km × 10 km) and regional (290 km × 290 km) spatial grains, but find that the former consistently provides a better fit to the data. In our case study, the effect of species richness on the slope SARs appears to be scale dependent, being negative at local scales but positive at regional scales. In univariate tests, the slope of the SAR correlates negatively with human population density and environmental energy availability, and positively with the rate of increase in habitat heterogeneity. We conducted two sets of multiple regression analyses, with and without species richness as a predictor. When species richness is included it exerts a dominant effect, but when it is excluded temperature has the dominant effect on the slope of the SAR, and the effects of other predictors are marginal.

Tree species-area and species-diameter relationships at three lowland rain forest sites in Sumatra

Data from three forest sites in Sumatra (Batang Ule, Pasirmayang and Tebopandak) have been analysed and compared for the effects of sample area cut-off, and tree diameter cut-off. An 'extended inverted exponential model' is shown to be well suited to fitting tree-species-area curves. The model yields species carrying capacities of 680 for Batang Ule, 380 species for Pasirmayang, and 35 for Tebopandak (tree diameter >10cm). It would seem that in terms of species carrying capacity, Tebopandak and Pasirmayang are rather similar, and both less diverse than the hilly Batang Ule site. In terms of conservation policy, this would mean that rather more emphasis should be put on conserving hilly sites on a granite substratum. For Pasirmayang with tree diameter >3cm, the asymptotic species number estimate is 567, considerably higher than the estimate of 387 species for trees with diameter >10cm. It is clear that the diameter cut-off has a major impact on the estimate of the species carrying capacity. A conservative estimate of the total number of tree species in the Pasirmayang region is 632 species! In sampling exercises, the diameter cut-off should not be chosen lightly, and it may be worth adopting field sampling procedures which involve some subsampling of the primary sample area, where the diameter cut-off is set much lower than in the primary plots.

The imprint of the geographical, evolutionary and ecological context on species-area relationships

Species-area relationships (SAR) are fundamental in the understanding of biodiversity patterns and of critical importance for predicting species extinction risk worldwide. Despite the enormous attention given to SAR in the form of many individual analyses, little attempt has been made to synthesize these studies. We conducted a quantitative meta-analysis of 794 SAR, comprising a wide span of organisms, habitats and locations. We identified factors reflecting both pattern-based and dynamic approaches to SAR and tested whether these factors leave significant imprints on the slope and strength of SAR. Our analysis revealed that SAR are significantly affected by variables characterizing the sampling scheme, the spatial scale, and the types of organisms or habitats involved. We found that steeper SAR are generated at lower latitudes and by larger organisms. SAR varied significantly between nested and independent sampling schemes and between major ecosystem types, but not generally between the terrestrial and the aquatic realm. Both the fit and the slope of the SAR were scale-dependent. We conclude that factors dynamically regulating species richness at different spatial scales strongly affect the shape of SAR. We highlight important consequences of this systematic variation in SAR for ecological theory, conservation management and extinction risk predictions.

Tree species composition and environmental relationships in a Neotropical swamp forest in Southeastern Brazil

We examined the relationships between topography, soil properties and tree species composition in a Neotropical swamp forest in southeastern Brazil. Plots were sampled in the forest, encompassing three different soil ground water regimes along the topographical declivity. All non-climbing plant individuals with trunk height >1.3 m were sampled. A canonical correspondence analysis-CCA-of the species-environmental relationships grouped tree species according to drainage and chemical soil conditions. A total of 86 species were found, being 77 species in the inferior, 40 species in the intermediate and 35 species in the superior topographic section. Some species were among the 10 most abundant ones, both in the overall sampled area and in each topographical section, with alternation events occurring only with their abundance position. However, substantial differences in floristic composition between sections were detected in a fine spatial scale, due to higher number of species, diversity index (H′) and species unique (exclusives) in the inferior topographic section. These higher values can be attributed to its higher spatial heterogeneity that included better drained and seasonally waterlogged soils, higher soil fertility and lower acidity. The increase of the soil water saturation and the uniform conditions derived from the superficial water layer has led to a lower number of species and an increase on the palm trees abundance in the intermediate and superior sections. Our results showed that at a small spatial scale niche differentiation must be an important factor related to the increase of the local diversity. The wide distribution of the most abundant species in the studied area and the increase of local diversity corroborate the pattern of distribution of species in larger scales of swamp forests, in which the most abundant species repeat themselves in high densities in different remnants. However, the floristic composition of each remnant is strongly variable, contributing to the increase of regional diversity. © 2008 Springer Science+Business Media B.V.

Unraveling the drivers of community dissimilarity and species extinction in fragmented landscapes

Communities in fragmented landscapes are often assumed to be structured by species extinction due to habitat loss, which has led to extensive use of the species-area relationship (SAR) in fragmentation studies. However, the use of the SAR presupposes that habitat loss leads species to extinction but does not allow for extinction to be offset by colonization of disturbed-habitat specialists. Moreover, the use of SAR assumes that species richness is a good proxy of community changes in fragmented landscapes. Here, we assessed how communities dwelling in fragmented landscapes are influenced by habitat loss at multiple scales; then we estimated the ability of models ruled by SAR and by species turnover in successfully predicting changes in community composition, and asked whether species richness is indeed an informative community metric. To address these issues, we used a data set consisting of 140 bird species sampled in 65 patches, from six landscapes with different proportions of forest cover in the Atlantic Forest of Brazil. We compared empirical patterns against simulations of over 8 million communities structured by different magnitudes of the power-law SAR and with species-specific rules to assign species to sites. Empirical results showed that, while bird community composition was strongly influenced by habitat loss at the patch and landscape scale, species richness remained largely unaffected. Modeling results revealed that the compositional changes observed in the Atlantic Forest bird metacommunity were only matched by models with either unrealistic magnitudes of the SAR or by models ruled by species turnover, akin to what would be observed along natural gradients. We show that, in the presence of such compositional turnover, species richness is poorly correlated with species extinction, and z values of the SAR strongly underestimate the effects of habitat loss. We suggest that the observed compositional changes are driven by each species reaching its individual extinction threshold: either a threshold of forest cover for species that disappear with habitat loss, or of matrix cover for species that benefit from habitat loss.

Sampling-Design Effects on Properties of Species-Area Relationships - A Case Study from Estonian Dry Grassland Communities

Despite widespread use of species-area relationships (SARs), dispute remains over the most representative SAR model. Using data of small-scale SARs of Estonian dry grassland communities, we address three questions: (1) Which model describes these SARs best when known artifacts are excluded? (2) How do deviating sampling procedures (marginal instead of central position of the smaller plots in relation to the largest plot; single values instead of average values; randomly located subplots instead of nested subplots) influence the properties of the SARs? (3) Are those effects likely to bias the selection of the best model? Our general dataset consisted of 16 series of nested-plots (1 cm(2)-100 m(2), any-part system), each of which comprised five series of subplots located in the four corners and the centre of the 100-m(2) plot. Data for the three pairs of compared sampling designs were generated from this dataset by subsampling. Five function types (power, quadratic power, logarithmic, Michaelis-Menten, Lomolino) were fitted with non-linear regression. In some of the communities, we found extremely high species densities (including bryophytes and lichens), namely up to eight species in 1 cm(2) and up to 140 species in 100 m(2), which appear to be the highest documented values on these scales. For SARs constructed from nested-plot average-value data, the regular power function generally was the best model, closely followed by the quadratic power function, while the logarithmic and Michaelis-Menten functions performed poorly throughout. However, the relative fit of the latter two models increased significantly relative to the respective best model when the single-value or random-sampling method was applied, however, the power function normally remained far superior. These results confirm the hypothesis that both single-value and random-sampling approaches cause artifacts by increasing stochasticity in the data, which can lead to the selection of inappropriate models.

Predicting species diversity in tropical forests

A fundamental question in ecology is how many species occur within a given area. Despite the complexity and diversity of different ecosystems, there exists a surprisingly simple, approximate answer: the number of species is proportional to the size of the area raised to some exponent. The exponent often turns out to be roughly 1/4. This power law can be derived from assumptions about the relative abundances of species or from notions of self-similarity. Here we analyze the largest existing data set of location-mapped species: over one million, individually identified trees from five tropical forests on three continents. Although the power law is a reasonable, zeroth-order approximation of our data, we find consistent deviations from it on all spatial scales. Furthermore, tropical forests are not self-similar at areas ≤50 hectares. We develop an extended model of the species-area relationship, which enables us to predict large-scale species diversity from small-scale data samples more accurately than any other available method.

The geographical structure of British bird distributions:diversity, spatial turnover and scale

1. Using data on the spatial distribution of the British avifauna, we address three basic questions about the spatial structure of assemblages: (i) Is there a relationship between species richness (alpha diversity) and spatial turnover of species (beta diversity)? (ii) Do high richness locations have fewer species in common with neighbouring areas than low richness locations?, and (iii) Are any such relationships contingent on spatial scale (resolution or quadrat area), and do they reflect the operation of a particular kind of species-area relationship (SAR)?

2. For all measures of spatial turnover, we found a negative relationship with species richness. This held across all scales, with the exception of turnover measured as beta (sim).

3. Higher richness areas were found to have more species in common with neighbouring areas.

4. The logarithmic SAR fitted better than the power SAR overall, and fitted significantly better in areas with low richness and high turnover.

5. Spatial patterns of both turnover and richness vary with scale. The finest scale richness pattern (10 km) and the coarse scale richness pattern (90 km) are statistically unrelated. The same is true of the turnover patterns.

6. With coarsening scale, locations of the most species-rich quadrats move north. This observed sensitivity of richness 'hotspot' location to spatial scale has implications for conservation biology, e.g. the location of a reserve selected on the basis of maximum richness may change considerably with reserve size or scale of analysis.

7. Average turnover measured using indices declined with coarsening scale, but the average number of species gained or lost between neighbouring quadrats was essentially scale invariant at 10-13 species, despite mean richness rising from 80 to 146 species (across an 81-fold area increase). We show that this kind of scale invariance is consistent with the logarithmic SAR.

Influence of multi-scale landscape structure on the occurrence of carnivorous mammals in a human-modified savanna, Brazil

So Paulo is the most developed state in Brazil and contains few fragments of native ecosystems, generally surrounded by intensive agriculture lands. Despite this, some areas still shelter large native animals. We aimed at understanding how medium and large carnivores use a mosaic landscape of forest/savanna and agroecosystems, and how the species respond to different landscape parameters (percentage of landcover and edge density), in a multi-scale perspective. The response variables were: species richness, carnivore frequency and frequency for the three most recorded species (Puma concolor, Chrysocyon brachyurus and Leopardus pardalis). We compared 11 competing models using Akaike`s information criterion (AIC) and assessed model support using weight of AIC. Concurrent models were combinations of landcover types (native vegetation, ""cerrado"" formations, ""cerrado"" and eucalypt plantation), landscape feature (percentage of landcover and edge density) and spatial scale. Herein, spatial scale refers to the radius around a sampling point defining a circular landscape. The scales analyzed were 250 (fine), 1,000 (medium) and 2,000 m (coarse). The shape of curves for response variables (linear, exponential and power) was also assessed. Our results indicate that species with high mobility, P. concolor and C. brachyurus, were best explained by edge density of the native vegetation at a coarse scale (2,000 m). The relationship between P. concolor and C. brachyurus frequency had a negative power-shaped response to explanatory variables. This general trend was also observed for species richness and carnivore frequency. Species richness and P. concolor frequency were also well explained by a second concurrent model: edge density of cerrado at the fine (250 m) scale. A different response was recorded for L. pardalis, as the frequency was best explained for the amount of cerrado at the fine (250 m) scale. The curve of response was linearly positive. The contrasting results (P. concolor and C. brachyurus vs L. pardalis) may be due to the much higher mobility of the two first species, in comparison with the third. Still, L. pardalis requires habitat with higher quality when compared with other two species. This study highlights the importance of considering multiple spatial scales when evaluating species responses to different habitats. An important and new finding was the prevalence of edge density over the habitat extension to explain overall carnivore distribution, a key information for planning and management of protected areas.