Abstract

Large areas of tropical forest now exist as remnants scattered across agricultural landscapes, and so understanding the impacts of forest fragmentation is important for biodiversity conservation. We examined species richness and nestedness among tropical forest remnants in birds (meta-analysis of published studies) and insects (field data for fruit-feeding Lepidoptera (butterflies and moths) and ants). Species–area relationships were evident in all four taxa, and avian and insect assemblages in remnants typically were nested subsets of those in larger areas. Avian carnivores and nectarivores and predatory ants were more nested than other guilds, implying that the sequential loss of species was more predictable in these groups, and that fragmentation alters the trophic organization of communities. For butterflies, the ordering of fragments to achieve maximum nestedness was by fragment area, suggesting that differences among fragments were driven mainly by extinction. In contrast for moths, maximum nestedness was achieved by ordering species by wing length; species with longer wings (implying better dispersal) were more likely to occur at all sites, including low diversity sites, suggesting that differences among fragments were driven more strongly by colonization. Although all four taxa exhibited high levels of nestedness, patterns of species turnover were also idiosyncratic, and thus even species-poor sites contributed to landscape-scale biodiversity, particularly for insects.

1. Introduction

Tropical forests are becoming increasingly fragmented by the widespread and rapid intensification of anthropogenic activities [1–4]. Formerly extensive tracks of continuous forests now exist as patchworks of isolated remnants scattered across inhospitable landscapes of non-forest habitats. This has resulted in the remaining forest patches supporting increasingly isolated populations of forest-dependent species [5]. Tropical rainforests support the majority of global biodiversity and contain large numbers of endemic species, and so understanding the impacts of forest fragmentation in these areas is crucial to the conservation of biodiversity [6–8].

There is a large literature examining the effects of tropical forest fragmentation on animals [1,9], including vertebrates (e.g. [10–14]) and insects [15–20]. Many previous studies have focused on describing changes in species richness following fragmentation, and have typically assessed the conservation value of forest remnants through the use of species–area relationships (SARs), which measure the rate at which species richness is lost as fragment size declines [21]. However, it is not clear how consistent SARs are within or among different taxa. Moreover, SARs may give a distorted view of the broader impacts of fragmentation on forest remnants because they provide no information on how faunal composition is altered by habitat fragmentation. Species losses may be offset by replacements from the surrounding matrix or elsewhere, so that SARs may underestimate changes in species composition following fragmentation, but it is unclear to what extent this might affect SARs in fragmented tropical forests.

An additional commonly observed property of biotas in fragmented habitats is that they are nested, such that species present at species-poor sites are subsets of those present at more species-rich sites [22,23]. Examination of the degree of nestedness has been used to infer the degree to which changes in species richness following fragmentation follow a predictable sequence [10,24–29]. However, few studies have evaluated whether different taxonomic groups follow the same patterns of nestedness or whether these patterns are generated by the same underlying mechanisms in each case [23].

Nestedness of species assemblages in forest fragments may arise as a result of four different processes (passive sampling, selective extinction, selective colonization and habitat nestedness [29]), which predict different patterns of nestedness among species and among fragments. Passive sampling could generate nestedness as a consequence of the fact that rare species are less likely to be sampled in a given area than common species, leading to the prediction that species with low average abundance will be a nested subset of those with high average abundance. Selective extinction predicts that in systems experiencing species loss or ‘relaxation’ post-fragmentation, area will be the main factor explaining species nestedness because species with large minimum-area requirements have greater extinction risk [30]. Selective colonization predicts that fragment isolation will create nested subsets of species through dispersal limitation [31,32]. Habitat nestedness considers the nestedness of species assemblages to arise as a result of their reliance on habitats that have a nested distribution [30]. Differences in the processes generating nestedness have important implications for conservation and can be used to direct management efforts [29], but are only poorly understood.

In addition to the characteristics of fragments, the traits of species may also provide useful information for assessing the importance of different processes in generating patterns of nestedness. For instance, if susceptibility to extinction is a strong determinant of nestedness, then traits linked to extinction proneness (e.g. geographical range size, habitat specificity [33,34]) may order species occurrence patterns. In contrast, if dispersal ability is a major driver, then traits linked to vagility (e.g. wing length [35,36]) may be more important in ordering species. However, this approach has received little attention [29].

Fragmentation of habitats has been shown to affect biotic interactions between species [37], such as herbivory [38–42] and seed predation [43]. There is also evidence that fragmentation can lead to trophic cascades [44–46], raising concerns that remnants of forest may not be viable in the longer term [47,48]. Nonetheless, small remnants may be important for promoting diversity in highly modified landscapes [49], and may also boost pollination success in surrounding agricultural crops (e.g. coffee [50,51]). Hence, patterns of nestedness among different trophic guilds may have important consequences for ecosystem functioning, but these have seldom been considered.

In this paper, we use information from published studies together with our own field data to examine how different ecological and functional groups of birds and insects respond to tropical forest fragmentation. There is an extensive literature on avian responses to fragmentation, as well as data on ecological traits (e.g. [14,52,53]), which we use to examine whether there is a consensus in SARs and patterns of nestedness among a wide range of tropical bird studies and among different feeding guilds. Most studies to date have focused on single taxa, but it is unclear the degree to which species and taxon responses are idiosyncratic or predictable [54,55]. Hence, we also analyse field data for butterflies, moths and ants from forest fragments in Sabah (Malaysian Borneo). The insect data were collected from many of the same forest fragments, allowing us a rare opportunity to evaluate the importance of different processes in generating SARs and patterns of nestedness among three insect taxa at the same locations.

2. Methods

(a) Avian responses to fragmentation

(i) Source data

We searched for suitable peer-reviewed studies examining avian responses to tropical forest fragmentation in online bibliographic databases (ISI Web of Science, EDINA Biosis); references within these were a source of further studies for consideration. Studies were included if they enumerated the species richness of each forest fragment studied, or provided data from which this could be obtained. Studies were included only if they provided data from three or more fragments. If fragments were compared at more than one distinct geographical location within a single published study, then data for each location were included separately such that each study provided more than one dataset for analysis. The area and isolation (distance from continuous forest) of fragments were taken from the authors' accounts or estimated from figures and maps using SigmaScan software (www.sigmaplot.com). We also recorded the surrounding landscape matrix as either cultivated, grassland, water or other (urban and mixed matrices).

(ii) Species richness within fragments

To evaluate the consistency in avian responses to tropical forest fragmentation, we first used linear regression to parameterize and determine the statistical significance of the SAR in each study. Methods for measuring species richness of fragments varied among the published studies and so, to ensure as far as possible that sampling effort was equal among fragments, we used individual-based rarefaction [56] to standardize sample sizes across fragments within each study. In some studies, species richness estimates for fragments were the result of historical inventories rather than discrete sampling regimes and these studies were excluded from the analysis. We then combined data across studies (N = 26) and used general linear modelling to relate the slopes (z) of the SARs to factors associated with extinction (range of fragment areas encompassed by study (i.e. max–min), and colonization (average isolation of fragments from continuous forest). We also controlled for possible confounding effects by including: habitat of the surrounding landscape matrix, time since fragmentation, duration and method of sampling (mist-net or observation), latitude, elevation and region (Neotropic, Africa or Asia) of study in the analysis. All predictor variables were included in the initial model, and the final model was obtained following serial deletion of least significant factors [57].

(iii) Nestedness and avian ecological traits

To examine the importance of different processes in generating patterns of nestedness, we first generated a species by site matrix for each study, indicating the presence or absence of each species within each fragment. We then used Binmatnest software [22] to reorder the matrix to achieve maximum nestedness among fragments and to compute a temperature index (T) scaled to vary between 0 (perfect nestedness, ‘cold’) and 100 (random order, ‘hot’). We chose this method over commonly used alternatives (e.g. the nestedness temperature calculator; [58]) because Binmatnest is a more robust test [29]. As recommended by the authors [22], we used null model 3 to examine whether the assemblages of birds in the different studies were significantly nested. Within each study, we then used Spearman correlations to examine whether or not the ranked order of fragments within the maximally packed matrix was related to fragment area. We next combined data across studies to examine how the degree of nestedness (temperature index, T) was related to variables linked to extinction, colonization or potential confounding effects (using general linear modelling with the same predictor variables as for the analysis of the slopes of SARs above). We based our analysis of nestedness on presence–absence rather than abundance data [59] because we were primarily interested in extinction and colonization processes rather than those generating variation in abundance within occupied sites.

To examine how the life-history traits of species influenced their degree of nestedness among fragments, we conducted separate analyses for subsets of species differing in body mass, geographical range size and feeding guild. Data for body mass were from [60] and data for species' geographical range sizes were from Birdlife International (www.birdlife.org). Species were classified into one of six feeding guilds (carnivore (including raptor, piscivore, faunivore, predator), frugivore, granivore, insectivore (including ant-follower, bark/foliage-gleaner, woodpecker), nectarivore, and omnivore (including scavenger, ‘miscellaneous’, and any species spanning two or more guilds)) according to their diet as reported by authors in the studies. For each study with suitable data, the nestedness of each feeding guild was quantified and compared with the nestedness of the total assemblage. For analyses of body mass and range size, assemblages were divided into two according to the median value, and nestedness was compared between the two groups.

(b) Insect responses to forest fragmentation

(i) Study sites

Leaf litter ants and fruit-feeding moths and butterflies were sampled at 13 sites in Sabah (Northern Borneo), comprising two widely separated ‘control’ sites within a single large area (approx. 1 million ha) of continuous forest (Danum Valley and Maliau Basin protected areas) and eight (Lepidoptera) or nine (ants) sites within forest remnants varying in size from 45 ha to approximately 120 000 ha (figure 1). Sample locations were all in lowland eastern Sabah, to reduce any residual effects of β-diversity prior to fragmentation. At the time of sampling, forest fragments were surrounded by an agricultural landscape containing mainly oil palm (Elaeis guineensis Jacq.) plantations, and had probably been isolated for approximately 50 years [61]. Isolation of each fragment, measured in terms of the minimum distance to continuous forest using a regional map of Sabah, varied from 6 to 69 km (figure 1).

(ii) Insect sampling

Obtaining complete inventories of species at each site was not possible and so, following previous studies [62], we used a standardized sampling technique to determine species richness in terms of the number of species in a uniform sampling area. We used standardized techniques to sample species along 2 km transects established at each study site. All samples were from locations more than 200 m within forest, to reduce edge effects [63]. Fruit-feeding moths (Families Noctuidae and Geometridae) and butterflies (Family Nymphalidae) were sampled with traps baited with rotting banana, hung at stations at 100 m intervals along each transect. This guild comprises about 75 per cent of all nymphalid butterflies recorded on Borneo [64,65] but a much smaller proportion of noctuid and geometrid moths [66–73]. Traps were suspended 2 m above ground level at each station and fresh banana was added to the trap each day to ensure a mixture of fresh and well-rotted bait [74]. Traps were emptied daily for 12 consecutive days, and sampling was carried out on two occasions covering both wet and dry seasons, to account for potential seasonal variation in species abundance [75]; 20 traps × 24 days × 10 sites = 4800 trap-days in total).

Leaf litter ants were sampled at stations every 500 m along each transect. Leaf litter and loose topsoil were collected from five 1 m2 quadrats positioned 10 m apart (forming a cross shape) at each station (five quadrats × 5 stations × 11 sites = 275 quadrats in total). Quadrats were placed at least 10 m apart to ensure that sampled individuals were likely to have come from different nests [16,76,77]. Soil and leaf litter was sieved (mesh = 1 cm2) to remove debris and combined into a single sample per station. These samples were then emptied into mesh bags and hung inside modified Winkler bags in the shade for 3 days to extract ants [78], which were stored in 95 per cent alcohol.

(iii) Species identification

All trapped butterflies were identified where possible in the field (by S.B.; following Otsuka [79]), marked with a felt-tipped pen and released. Recaptured individuals were excluded from subsequent analyses. Individuals of the genera Tanaecia and Euthalia cannot be identified reliably in the field and were collected and identified in the laboratory using published keys and figures [79,80]; where necessary, this included dissection of male genitalia. All trapped moths were collected and identified in the laboratory (by S.B. and C.V.K.) using reference collections at the Forest Research Centre, Sabah, and published keys and figures [66–71].

Ants were identified (by N.T.) first to genus following existing keys [81,82], then to species using reference collections at the Natural History Museum (London), and Borneensis Collection (Universiti Malaysia Sabah), and online resources (http://www.antweb.org; http://www.antbase.net). Undescribed species that occur in www.antweb.org were given the same number as the online collection, and morphospecies not featured were given new reference numbers and submitted to the collection. Voucher specimens of ants and Lepidoptera have been deposited at the Forest Research Centre, Sandakan, Sabah.

(iv) Species richness and nestedness

Species richness is highly sensitive to sample size and even though sampling effort was standardized, sampling efficiency or total abundance may have varied across transects. The numbers of species recorded on different transects could therefore have reflected the numbers of individuals sampled rather than the numbers of species inhabiting different patches. Hence, we used a jackknife estimate to assess the probable species pool at each site (Smax ± 1 s.e.) from the number of species observed (Sobs), using the formula Smax = Sobs + a(n−1/n) where n is the number of samples and a the number of species in only one sample [83]. In addition to species richness at each study site, we also analysed the nestedness of assemblages across sites, using Binmatnest software (see above). Data from the two ‘control’ sites in continuous forest (Danum Valley and Maliau Basin) were combined prior to running Binmatnest. In keeping with previous studies (e.g. [27]), we restricted this analysis to species sampled five or more times in total, to increase the reliability of the presence/absence data in our species by site matrices. In addition, to check how robust our results were, we also used Nodf software [59] to analyse nestedness for all species sampled, using quantitative abundance-based matrices.

(v) Vegetation structure in fragments

To investigate relationships between nestedness and vegetation structure within each forest fragment, we measured the following variables at each sampling station (slightly different methods were used for ants and Lepidoptera). For Lepidoptera (n = 200 stations), we measured the height of the two nearest large trees (girth at breast height (g.b.h.) more than 60 cm) within 30 m in each of four quadrants centred at the station (maximum of eight trees per station), the number of quadrants with one or more trees within 30 m (maximum of four) and we estimated vegetation cover (%) at ground and canopy levels (see [27] for further details). For ants (n = 55 stations), we measured the girth (in cm) and distance (in m) of the eight nearest large trees (g.b.h. more than 60 cm) within 20 m of the centre of each station, and the number of these trees that were dipterocarps (which is a good indicator of the degree to which forest habitats within fragments are similar to undisturbed habitats within continuous forest [84]). We also measured leaf litter depth (cm) at each station (mean of five measurements, taken within each 1 m2 quadrat with a wooden rule prior to sampling of ants), and we estimated vegetation cover (%) at canopy and ground levels. Vegetation measurements were normalized where necessary (i.e. arcsine transformation of percentages) and analysed by principal components analysis (PCA [84]; Lepidoptera and ants were sampled at slightly different locations along transects and so the two datasets were analysed separately).

(vi) Ecological and life-history traits

To gauge the importance of stochastic effects in driving patterns of nestedness, we calculated the average abundance of each species of butterfly and moth, across those sites where each species occurred. We did not do this for ants, because individuals from each nest are highly aggregated and non-independent; hence, numbers of individuals sampled at a site are unlikely to indicate the probability of recording a species as present within other occupied sites. As an index of vagility, we measured butterfly forewing length in the field (±1 mm using vernier calipers), and obtained data on moth forewing length from published sources [66–71]. To assess the vulnerability of species to local extinction, butterflies and moths were assigned to four geographical range size groups ([27]; endemic to Borneo, restricted to Sundaland, restricted to the Oriental region, occurring beyond the Oriental region (using data from [66–71,85])). To examine host–plant diversity, we assigned species to two groups in terms of larval host plant specificity (narrow, all host plants within a single plant family; broad, exploits two or more families of host plants; data from [86]).

To examine how the nestedness of ant species was related to different ecological traits, sampled ants were assigned to one of four feeding guilds (predator, granivore, omnivore, scavenger) and four nesting site guilds (soil, leaf litter, decaying wood, generalist) using published information ([87] and the Australian Ant Online database (http://anic.ento.csiro.au/ants/). Ants were assigned to guilds at genus level because detailed information at species level is not currently available. We then analysed nestedness separately for each guild.

3. Results

(a) Avian responses to fragmentation

(i) Species–area relationships

Data were extracted from 26 published studies for subsequent re-analysis (see the electronic supplementary material, appendix S1). Studies examined between 3 and 35 fragments (median = 6 fragments), and all three major tropical regions were represented. Significant SARs were present in 13 out of 26 studies (50%; double-log power function, log S/log A), of which the range of slope values (z) was 0.048–0.81 and the range of intercept values (c) was 0.81–1.76 (untransformed: 6.5–57.5 ha). The number of fragments within a study was not significantly associated with any regression parameters (p > 0.05 in each case), and isolation did not explain a significant amount of variation in S beyond that accounted for by area (p > 0.05 in all cases).

However, studies with a wider range of isolation values among fragments had steeper slopes (rs = 0.36, p < 0.05) as did those at lower latitudes (rs = −0.57, p < 0.01). Studies from Africa had significantly shallower slopes than did studies from the Neotropics or Asia (Kruskal–Wallis: χ2 = 9.2, p < 0.01; Steel-type multiple comparison: p < 0.05). No other variables showed significant relationships with the slopes of SARs.

(ii) Nestedness

Data suitable for inclusion into nestedness analyses were obtained for 1570 records (1057 species) from 15 studies. Data on geographical range sizes were available for all species, and data on body masses were available for 876 species. All but two of the 15 studies were significantly nested, and nestedness index values (T) across studies were low (range = 5.0–35.3, mean ± 95% CI = 19.3 ± 6.3) indicating a high degree of nestedness in most cases (figure 2). For those studies showing nestedness, the rank order in which fragments were arranged to achieve maximum nestedness showed a significant negative relationship with fragment area in seven of 14 cases (figure 2c); a further five studies, mostly with few fragments, showed non-significant negative associations with area.

The types of analyses we carried out on published datasets for birds. Data [88] showed significant patterns of species richness and composition among fragments. (a) Fragments arranged to achieve maximum nestedness of species presences (grey cells) and absences (white cells) in the fragment (rows) × species (columns) matrix. Values to the left of the matrix show fragment areas (ha), all other values show species richness. The matrix was significantly nested (T = 4.48, p < 0.001). (b) SAR fitted by the power function: log10S = log10 0.81 (s.e. ± 0.05) + 0.254 (s.e. ± 0.05) × log10 area. R2 = 0.76, p = 0.001. (c) The order in which fragments were arranged to achieve maximum nestedness was negatively correlated with their area (rs = −0.95, p < 0.001).

Studies that were more highly nested (i.e. lower T values) had steeper SAR slopes (rs = −0.65, p < 0.05) indicating that greater net rates of species loss following fragmentation were associated with more predictable sequences of species turnover. However, the degree of nestedness was not significantly related to the average isolation of fragments or to variation in sampling methodology or location of study (general linear model; p > 0.1 in all cases). Separate analysis for different subsets of species indicated that the degree of nestedness was significantly different among the six avian feeding guilds (χ2 = 11.1, p < 0.05) with carnivores and nectarivores significantly more nested than either granivores or insectivores (Steel-type multiple comparison: p < 0.01; figure 3). However, there was no significant difference in the nestedness of large and small bodied species (Mann–Whitney; z = 1.48, p = 0.2) or between species with large and small geographical ranges (z = 0.60, p = 0.6).

(b) Responses of insects to forest fragmentation

(i) Differences in vegetation structure among sites

The PCA analyses of vegetation structure for sites where Lepidoptera and ants were sampled (n = 200 stations and 55 stations, respectively) produced one factor score for Lepidoptera sites (LEP-PRIN1) that accounted for 71 per cent of the variation in the data, and three independent factor scores for ant sites (ANT-PRIN1, ANT-PRIN2, ANT-PRIN3) that together explained 76 per cent of variation. For Lepidoptera sites, a high LEP-PRIN1 score represented tall dense forest with a closed canopy and little ground vegetation. For ant sites, a high ANT-PRIN1 score represented forest with many large trees, a high ANT-PRIN2 score represented forest with high canopy cover and many dipterocarps, and a high ANT-PRIN3 score represented forest with deep leaf litter and low ground cover. LEP-PRIN1 and ANT-PRIN2 were both positively related to fragment size (Spearman correlation, Lepidoptera, rs = 0.69, n = 10, p < 0.05; ants, rs = 0.61, n = 11, p < 0.05) showing that larger fragments were more likely to support closed-canopy forest with a high proportion of dipterocarp trees.

(ii) Species richness in forest fragments

Excluding recaptures, we sampled 6833 individuals of 84 species of butterfly, 2112 individuals of 107 species of moth and 19 998 individuals of 200 species/morphospecies of ant (figure 4; see the electronic supplementary material, appendices S3 and S4). Estimation of Smax indicated that we sampled 77 per cent of the species present at each site for butterflies (range = 69–84%), 70 per cent for moths (range = 64–79%) and 71 per cent for ants (range = 67–81%). Excluding continuous forest locations, Smax for butterflies was significantly and positively related to log10 fragment area (figure 5) and negatively related to isolation (general linear model; F2,5 = 35.9, p = 0.001) but not to vegetation structure. For both moths and ants, there were significant positive relationships between Smax and log10 fragment area (figure 5; moths, F1,6 = 8.1 p < 0.05; ants, F1,6 = 13.4, p = 0.01) but no relationships with fragment isolation or vegetation structure (p > 0.3 in all cases). Comparison of SARs among the three insect taxa revealed very similar z values, although intercept values varied among taxa (species richness of butterflies < moths < ants; figure 5).

Final estimated size of the species pool Smax ± 1 s.d. for leaf-litter ants (triangles) and fruit-feeding moths (open circles) and butterflies (solid circles) in rainforest fragments in Sabah, Borneo, in relation to log10 fragment area.

(iii) Patterns of nestedness

Excluding species with fewer than five individuals, we analysed patterns of nestedness for 48 of 84 recorded species of butterflies, 40 of 108 species of moths and 99 of 200 species of ants. The species occupancy data filled 49.4, 45.7 and 47.4 per cent of the overall species by site matrix for butterflies, moths and ants, respectively, and the observed temperatures of the maximally packed matrices (38.1°C, 31.7°C and 29.2°C, respectively) were significantly cooler than random (p < 0.01). Thus, all three insect taxa were significantly nested. The ant species data were then split according to feeding guild (four groups) and nesting locations (four groups), and the nestedness analysis was re-run separately for each group. Only predators were significantly nested (T = 15.0°C, n = 23 species, p = 0.01).

Examination of the species by site matrix showed that for butterflies, sites were nested by area, such that smaller sites contained species assemblages that were subsets of those present at larger sites (rs = −0.79, n = 9, p = 0.01), but this relationship was not evident in either moths or ants (p > 0.1). There was also no evidence for any insect taxon that sites were nested according to either site isolation or vegetation structure (p > 0.1). We also investigated whether the ordering of Lepidoptera in the species by site matrix was related to species traits. Moth species with large wings occurred at more sites, including those with few other species (rs = −0.58, n = 34, p < 0.001), but this effect was not evident in butterflies (rs = −0.25, n = 48, p = 0.1). In addition, there was some evidence that butterfly species with more restricted geographical ranges were confined mainly to those few sites with high species richness (Kruskal–Wallis χ2 = 5.30, n = 48, p = 0.07), but this effect was not evident in moths (χ2 = 3.5, n = 35, p = 0.3).

4. Discussion

(a) Comparisons of SARs among taxa

Our analysis of 26 bird studies, and data for more than 390 insect species from 13 sites in Sabah revealed a general consensus in terms of species responses to habitat fragmentation. Significant SARs were common; approximately half of the bird studies and all three insect studies produced significant positive relationships between log10S and log10 area. The three insect taxa were studied in fragments from the same region and so it was perhaps not surprising that the SARs were very similar. Nonetheless, figure 5 shows that the slopes were almost identical across the three taxa, and the intercept values reflected differences in overall species richness between the groups (ants > moths > butterflies). Our studies on birds came from a heterogeneous dataset encompassing a wide range of study regions, matrix types and sizes and numbers of fragments. However, with the exception of region and latitude, variation in SARs in our bird studies was not related to these potentially confounding variables. Thus, we detected more variation in SARs within a single taxon sampled in different locations than we found among taxa studied at the same site. In our bird studies, we found no relationship between time-since-fragmentation and slope of SARs, implying that any differences in relaxation times of assemblages following fragmentation (i.e. extinction debt; [89]) did not affect the slopes of SARs. Previous studies have estimated that it takes approximately 50 years for half the species to be lost from a 1000 ha fragment [90]. Given the large sizes of many of our study fragments, it is likely that these fragments were not at equilibrium at the time of sampling, implying more extinctions from fragmented landscapes in future than were recorded by these studies.

In all taxa, area explained most of the variation in species richness (S) among fragments. However, in birds, the slope of the SARs (z) was steeper in studies that exhibited a great range of isolation values, and in butterflies, isolation explained additional variation in species richness after accounting for area. In addition, there was no effect of vegetation structure on SARs in any insect taxon. Thus, we conclude that extinction process were the most important drivers of species turnover following fragmentation, but that colonization processes were also evident.

(b) Nestedness patterns in different taxa

A high degree of nestedness was evident in all study taxa, and birds (mean = 19°C) were more nested than insects (mean = 33°C), implying more predictable patterns of species loss following fragmentation in birds. In insects there was no relationship between nestedness and abundance, indicating that observed patterns of nestedness were not an artefact of incomplete species inventories. In 50 per cent of the bird studies, and in butterflies, sites were nested by area. For all four study taxa, there was no suggestion that fragment isolation or vegetation structure of fragments affected the ranking of sites in the maximally ordered matrix. Thus, the findings from analyses of nestedness support those from analyses of SARs, and we conclude that extinction processes are most important in driving changes in species assemblages of birds and butterflies following fragmentation.

We found that in birds, more nested systems were associated with steeper SARs. This indicates that greater reduction in the number of species following fragmentation was associated with more predictable losses of species. Compared with insects, there was little evidence for new species of birds turning up in small fragments that were not present in larger areas. Compared with insects, birds may be more vulnerable to fragmentation effects, if they do not disperse across landscape matrices and so fail to colonize small fragments, or if fragments are too small to support viable populations of many species. Thus responses of birds to fragmentation may be less stochastic than those of insects.

Nestedness of birds and ants differed among feeding guilds, with carnivorous and nectivorous birds and predatory ants being the most nested whereas granivorous birds were the least nested, indicating that they may be less sensitive to reductions in forest area or more likely to colonize from the matrix [24]. These results are consistent with previous studies of birds on impacts of commercial selective logging [91,92] and fragmentation [53] showing that avian feeding guilds differed in their sensitivity to habitat disturbance, and that granivores were least affected. The sensitivity of avian carnivores and predatory ants to fragmentation may reflect their relatively high trophic status, for example, if their feeding requirements mean they have large minimum-area requirements for persistence. The sensitivity of avian nectarivores to habitat fragmentation that we detected may reflect that they are generally small specialist species with restricted ranges, which might be expected to be more adversely affected by fragmentation ([24,93]; but see [53]). In contrast, avian granivores and some guilds of ants may be able to exploit resources in matrix habitats [91].

(c) Relative importance of extinction processes

Our analyses of SARs and nestedness patterns in birds and insects following fragmentation allow us to determine the relative importance of different processes in driving changes in species richness and composition in fragments. Significant SARs with relatively little impact of fragment isolation, and no effect of vegetation structure on these relationships implies a primary role of extinction (i.e. area) on species richness. However, area was only an important driver in nestedness patterns in birds and butterflies, implying other factors were important for moths and ants.

In insects, there was no evidence that the nestedness ranking of sites was associated with fragment isolation or vegetation structure, showing that colonization and habitat quality were not important drivers of nestedness. It could be that our data did not adequately capture the effects of forest fragmentation on habitat quality, although this seems unlikely given that there was a significant relationship between differences in vegetation structure and β-diversity of fruit-feeding butterflies among sites [27] (figure 6). The ranking of species of Lepidoptera in the maximally ordered site by species matrix showed that different traits were important for butterflies versus moths. For moths, species were nested by wing length showing that species with longer wings (which may be better dispersers [35,36]) occurred at a wide range of sites, including those small sites with few species present. This indicates that some moths are able to disperse across the agricultural matrix [94]. By contrast, butterfly species were not ranked by dispersal ability, but there was a trend for them to be nested by geographical range size, reflecting the fact that endemic species occurred only in the large, species-rich sites (area more than 4000 ha) and were not recorded in smaller fragments [27]. We interpret these results as evidence for colonization processes being important for moths, but extinction processes being more important for butterflies. Thus, even in two very closely related taxa from within a single feeding guild (fruit-feeding butterflies and moths), factors driving responses to fragmentation were different.

Maximally ordered species presence–absence matrix for (a) fruit-feeding moths and (b) fruit-feeding butterflies in rainforest fragments in Sabah, Borneo. Filled squares indicate presence, white squares indicate absence; two species of moth and five species of butterfly present at all nine sites have been compressed into a single column.

(d) Implications for conservation and the value of small fragments

Our analyses showed that bird and insect assemblages in tropical forest fragments were highly nested, particularly in birds, such that species in small fragments were subsets of those in larger sites. However, responses of species were also idiosyncratic, particularly in insects where new species occurred in small fragments that had not been present in larger fragments. Thus, conservation of small fragments may be less important for birds than for insects in terms of conserving regional diversity. This supports findings from other studies for birds which concluded that small fragments were of little conservation value [95,96] whereas studies on insects and plants have highlighted the conservation value of small fragments [27,49,97,98].

While small fragments may support species not present in large fragments, differences among feeding guilds in patterns of nestedness following habitat fragmentation indicate that ecosystem functioning is likely to be altered in small fragments. For example, herbivory is reduced in small fragments [38–40,42], and the long-term viability of fragments may be reduced if seedling recruitment is affected [48].

Resources for conserving species are usually focused on large sites [99], but the idiosyncratic responses of some taxa to fragmentation indicate that the conservation of small fragments may promote regional diversity [27], and also increase productivity in surrounding agricultural areas [100,101]. Small fragments often experience increased habitat disturbance (this study; [27]) and altered micro-climates, and whether or not changes in trophic organization affect the long-term viability of these remnants requires further study.

Acknowledgements

We thank Glen Reynolds, the Royal Society South East Asia Rainforest Research Programme, Nazirah Mustaffa, and staff at Danum Valley Field Centre for logistical support. The project received funding from the UK Government Darwin Initiative (Defra) and the UK Biotechnology and Biological Sciences Research Council. This paper constitutes publication number A/578 of the Royal Society South East Asia Rainforest Research Programme.

2006Impacts of rain forest fragmentation on butterflies in Northern Borneo: species richness, turnover and the value of small fragments. J. Appl. Ecol.43, 967–977.doi:10.1111/j.1365-2664.2006.01209.x (doi:10.1111/j.1365-2664.2006.01209.x)

2010Nestedness for different reasons: the distributions of birds, lizards and small mammals on islands of an inundated lake. Div. Distrib.16, 862–873.doi:10.1111/j.1472-4642.2010.00682.x (doi:10.1111/j.1472-4642.2010.00682.x)

2005Temporal variation in abundance and diversity of butterflies in Bornean rain forests: opposite impacts of logging recorded in different seasons. J. Trop. Ecol.21, 417–425.doi:10.1017/S0266467405002361 (doi:10.1017/S0266467405002361)

1997Species-specific predictions of the impact of habitat fragmentation: local extinction of birds in the box-ironbark forests of central Victoria, Australia. Biol. Conserv.82, 147–155.doi:10.1016/S0006-3207(97)00028-1 (doi:10.1016/S0006-3207(97)00028-1)