Selective logging is known to alter the structural and community composition of tropical forests and may disrupt plant-frugivore interactions. We hypothesized that even after a sufficient period of recovery, logged areas will not possess as complete a suite of species as an unlogged forest, the differences being more marked for biotically-dispersed species. Species of this functional group are expected to occur at lower densities, have lower species richness and diversity, and be smaller in logged forests. To quantify structural and functional differences in tree communities, we sampled 120 randomly placed plots, 60 each in logged and unlogged forest sites. We found significant differences in species richness and diversity between logged and intact forest. Within biotically-dispersed species, bird dispersed species showed a significant reduction in species richness. Consistent with previous studies, trees in logged forests were smaller, although overall density was not different between the two treatments. We posit that selective logging might have pervasive effects on functional aspects of tropical tree communities, which appear to persist even after two decades of logging cessation.
Introduction
Tropical moist forests display large degrees of plant-animal interactions. Tree communities in many tropical rainforests are shaped by the relative roles of dispersers, seed predators, and herbivores that limit seedling growth and recruitment [1,2]. Frugivorous dispersers make up a bulk of the vertebrate population, and up to 80% of tropical rainforest trees may be dispersed by these frugivores [3,4]. Human-induced disturbances such as complete or selective logging and hunting may disrupt these interactions and affect tree community composition and diversity of tropical forests [5].
Selective logging, although considered less harmful to biodiversity [6], is known to alter the physical and community structure of tropical forests [7,8]. In addition to well-established structural changes to tree size, basal area, and biomass persisting even 12 years after logging cessation [7,9], logging may also cause pervasive changes in species composition of tropical moist forests [9]. Long-term studies (>15 years) and predictions based on the trajectory of forest recovery in these sites indicate that regeneration of the full suite of trees may take up to a century following complete termination of logging [10].
Selective logging affects the community dynamics of tropical rainforest trees by creating gaps with altered light and micro-climate conditions, facilitating the growth of small-seeded, light-loving pioneer species [11,12]. Additionally, logging is known to reduce populations of old-growth species [12]. As these species are slow-growing, they may take many years to recover [8,13]. Many old-growth trees are also important fruiting resources for animals, which in turn disperse the seeds of these species. Logging may thus reduce overall densities and alter relative abundances of frugivorous dispersers either directly or indirectly through loss of food trees [1415–16]. Loss of dispersers may eventually result in reduction or loss of biotically-dispersed trees [1718–19], and change the functional composition of the forest towards abiotically-dispersed species.
We examined the implications of logging for the tree-community as a whole. Specifically, our study aimed to discern the differences in tree density, species richness and diversity between an unlogged and a selectively logged forest in north-east India, after a 16-year recovery period. Further, we sought to determine if there are finer-scale differences between biotically and abiotically-dispersed species between these two treatments, and if biotically-dispersed species had lower densities, species richness, and species diversity following logging. Within the biotically-dispersed group, we examined whether richness and diversity were differentially affected within guilds of bird and mammal-dispersed tree species.
Methods
Study area
This study was done from November 2009 to May 2010 in and around Pakke Wildlife Sanctuary and Tiger Reserve (Fig. 1, 861.95 km2; 92º 36' - 93º 09' E; 26º 54' - 27º 16' N) in western Arunachal Pradesh. The main vegetation type is classified as Assam-valley tropical semi-evergreen forest [20]. The dominant tree species are Polyalthia simiarum, Pterospermum acerifolium, Sterculia alata, Stereospermum chelonioides, Ailanthus grandis and Duabanga grandiflora [19]. A high percentage of tree species (64%) are animal-dispersed, and 12% of tree species are wind-dispersed [21].
There are at least 35 predominantly frugivorous bird species [21]. The Great Hornbill (Buceros bicornis), Wreathed Hornbill (Rhyticeros undulatus) and Oriental Pied Hornbill (Anthracoceros albirostris) are the largest avian frugivores in the area. The Mountain Imperial Pigeon (Ducula badia) and the Green Imperial Pigeon (Ducula aenea) are also important large-bodied avian frugivores. There is also a great diversity of smaller frugivores, with more than 30 species known to occur here [22]. The other species that predominantly depend on fruits and seeds are four species of viverrids (Arctictis binturong, Viverra zibetha, Paguma larvata and Paradoxurus hermaphroditus) and five species of rodents (Hystrix brachyura, Atherurus macrourus, Niviventer sp., Rattus sp., Berylmus sp.), excluding four species of diurnal tree squirrels [23].
We selected two blocks (c. 10 km2 each) each in logged (27º 01' N; 92º 91' E and 27º 01' ‘N; 92º57’ E) and unlogged (27º 0' N, 92º 35' E and 26º 56' N; 92º 58' E) areas based on information from previous systematic records (maintained in the form of working plans of the Forest Department), Departmental Timber Operations (DTO) and published literature [13,15]. Both logged sites (27º 01' N, 92º 91' E and 27º 01' ‘N, 92º57’ E) were located in the adjacent peripheries of the protected area and were commercially logged. Up to the late 1970s. The Supreme Court imposed a complete ban on logging in 1996 [17]. Our study sites were about 30 km apart and broadly comparable in terms of geology, rainfall, climate, vegetation, and topography [15]. The commonly logged species were Terminalia myriocarpa, Duabanga grandiflora, Artocarpus chaplasha, Amoora wallichii, Michelia spp., Dysoxylum procerum and Toona ciliata. Trees would be marked for felling depending on the volume sanctioned to each person by the forest department. Sawmills operating in the adjacent parts of the park were either directly supplied with timber, or residents and forest staff would sell their permits to sawmills, which led to timber being commercially exploited. Hence, although areas were largely selectively logged by the government, the irregularities in issuing permits for selective logging, also led to some amount of illegal logging.
Data collection
We set up 120 (10 lgr; 10 m plot) vegetation plots, 60 each in logged and unlogged forest. Within each treatment (logged and unlogged), we had two blocks of 10 km2, within which there were 30 plots randomly placed and separated by a distance of at least 500 m. Random selection of plots was performed using GIS software prior to undertaking fieldwork (Quantum GIS - Open Source Geographic Information System). We used the selected random points as the right lower vertex of the plot. In the protected area there was a 53 km along the east-west axis, from which plots were accessed laterally. In the adjacent reserve forests, the road segments were 18 km and 4 km in length, which we used to access our sites. There were no roads passing through our study plots, except occasional foot trails. Within these plots we enumerated all adult tree species > 25 cm diameter at breast height (DBH). We measured DBH and tree height for all species.
Statistical analyses
We used the statistical software R (version 2.11.0) [24]. We used package Vegan in Program R for calculating species richness, species diversity indices, and species accumulation curves.
Characterising forest structure and composition in logged and unlogged forests
We examined the changes due to logging on forest structure (DBH, tree density) and more importantly, on composition (species richness, diversity, and similarity). As tree height was positively correlated with DBH in both logged (Pearson's R = 0.66, p < 0.001) and unlogged areas (Pearson's R = 0.48, p < 0.001), we retained DBH as the response for further analyses. We used Student's t tests and analysis of variance (ANOVA) to compare the DBH of trees in logged and intact forest.
To examine the effects of logging on forest composition, we used ANOVA to examine differences in overall species richness between logged and intact forest. We used the respective values of the indices of diversity and evenness to compare the two treatments. We assessed this difference between two functional guilds - species dispersed by biotic and abiotic means. Within animal-dispersed species, we used one-way ANOVA to model the effects of logging on the species richness and diversity of trees dispersed by mammals alone, birds alone, and both by mammals and birds. We determined dispersal modes for trees from existing literature [21,25] and from focal tree watches [Velho et al, unpublished data].
We used Fisher's α, S = ln(1 + N// α), where S is the species richness, N is the number of individuals, ln is the natural logarithm, and a is an index of species diversity. We used Morisita-Horn index of dissimilarity to quantify extent of species overlap in logged versus unlogged plots, and between logged and intact forest. We also used the Shannon-Weiner index of diversity (H') to compare logged and intact forest for overall species composition and for differences between the entropy of species within the two functional guilds. The formula for this index is given as
We used Pielou's index J= H' / ln(S), to compare evenness within the two treatments, where H' is Shannon-Weiner index and S is the total number of species in the sample.Results
Structural differences between logged and unlogged habitats
Overall tree density did not differ significantly between the two habitats (mean logged = 601.5 ±83.5 S.E.; mean unlogged = 598.2 ±76.8 S.E.). Density of biotically-dispersed species, however, was higher in intact forest (meanlogged =448.4 ±78.4 S.E.; mean unlogged = 506.9 ±66.1 S.E) and abiotically-dispersed species displayed the opposite trend (mean logged = 141.6 ±62.2 S.E.; meanunlogged = 86.8 ±35.8 S.E), although the results were not significant (p>0.05). As expected, we found a significantly lower DBH in logged forests compared to intact forests (Student's t test, t = −6.28, df = 681.42, p< 0.001) (Table 1). The DBH of biotically and abiotically-dispersed tree species, separately, was also lower in logged forests (Table 1). Diameter at breast height of species that were bird-dispersed (One way ANOVA, F = 38.54, df = 396, p<0.001), and dispersed by both (One way ANOVA, F = 5.05, df = 124, p = 0.03), was significantly lower in logged forests, and DBH of mammal-dispersed species was nearly so (One way ANOVA, F = 3.18, df = 47, p = 0.08).
Table 1:
Summary of changes in structural and compositional aspects of the tree communities in logged versus intact forest. The numbers in subscript represent the standard error (SE) of the mean values presented. The mean values for species richness are estimates. We compared species richness and tree density using ANOVA, and DBH using Student's t test. Arrows represent increase (↑) or decrease (↓) in the attributes of the two functional guilds and statistically significant differences at p < 0.05 are indicated as *. For diversity indices we used the magnitude of effect and not p-values for comparison.
Overall compositional differences between logged and unlogged habitats
Species diversity was higher in intact forests (Fisher's α = 42.17 ± 5.57S.E., df = 94) compared to logged forests (Fisher's α = 27.6 ±.4.01 S.E., df = 72). Morisita-Horn index of dissimilarity was 0.35, indicating a 35% difference in species composition between logged and unlogged sites. Shannon-Wiener index of diversity (H') was 3.96 for unlogged forest and 3.94 for logged forest. Pielou's index of evenness in composition was higher for logged forest (0.93) compared to unlogged forest (0.87).
Although the species accumulation curves were not asymptotic in either case, the curves for logged and unlogged forests indicate that the estimated species richness may be lower for logged forests (Fig. 2, mean logged = 94.7 ±5.8 S.E.; mean unlogged = 139.3 ±8.6 S.E.). We formally tested the relationship between felling regime and estimates of species richness (using the Chao1 for individual plots) and found significant differences between logged and unlogged forests (One-way ANOVA, F = 9.79, df = 118, p = 0.002).
Dispersal mode: How does logging affect biotically versus abiotically dispersed plant species?
We found a lower number of biotically-dispersed species per plot in logged forests, although the results were not statistically significant (mean logged = 4.48 CI 4.22-4.77; mean unlogged = 5.17 CI 4.66-5.51; p = 0.14). The number of abiotically-dispersed individuals per plot were significantly higher in logged habitats compared to intact forests (mean logged= 1.42CI. 1.27-1.58; mean unlogged = 0.87 CI 0.73-1.03; p = 0.005). Diversity of biotically-dispersed species was higher for unlogged plots (Fisher's α = 25.1 ± 4.3 S.E.) compared to logged plots (Fisher's α = 11.1 ± 2.4 S.E.). For abiotically-dispersed species too, Fisher's α was higher for intact forest (Fisher's α = 15.8 ± 3.2 S.E) compared to logged forest (Fisher's α = 19.4 ± 3.9 S.E.). Morisita-horn index of dissimilarity for biotic species was 0.67 indicating 67% difference in tree species composition, and for abiotic species was 0.24, or a 24% difference in tree species composition.
Shannon-Wiener index (H') for biotically-dispersed guild was 3.1 for logged forest and 3.5 for unlogged forest. For abiotic species, the index was lower in unlogged forest (H'= 3.2) compared to logged forest (H'= 3.4). For biotically-dispersed species, Pielou's evenness index was 0.89 in logged areas compared to 0.85 in unlogged forest. Pielou's index was 0.93 in logged forest and 0.85 in unlogged forest for abiotic species.
The species accumulation curves reached an asymptote only for abiotically-dispersed species within the treatments. However, the accumulation curves for biotic species (Fig. 2, mean logged = 72.7 ±5.6 S.E.; mean unlogged = 112.4 ±7.6 S.E.) indicate lower expected species richness in logged areas. For abiotic species (Fig. 2, mean logged = 19.0 ±1.4 S.E.; mean unlogged = 24.9 ±2.6 S.E.), the estimated species richness was higher in logged forest. There were significant differences between the species richness estimates (Chao1) of both abiotically (One-way ANOVA, F = 4.98, df = 118, p = 0.03) and biotically-dispersed tree species (Fig. 3, One-way ANOVA, F = 4.44, df = 118, p = 0.04). Within biotic communities, species richness of bird dispersed trees was significantly lower in logged forest (One-way ANOVA, F = 5.84, df = 118, p = 0.02). There were no significant differences in the species richness of trees dispersed by mammals (One-way ANOVA, F = 0.24, df = 118, p = 0.63) and species dispersed by both birds and mammals (One-way ANOVA, F = 0.01, df = 118, p = 0.94).
Discussion
Our study is one of the few that has looked at the effect of logging differentially on abiotically and biotically-dispersed tree species at a community level. We found marked changes in the functional composition of logged forest compared to adjacent and geo-physically similar intact forest. Logged forests possessed only a subset of the species found in intact forest, the difference being more marked for the biotically-dispersed guild. Canopy openness following logging is expected to favor the growth of small-seeded pioneer species, many of which are wind-dispersed, a change in community composition that persists decades after logging [2627–28]. We did find a greater number of abiotically-dispersed trees per plot and greater density of this guild in logged forest. Diversity of abiotically-dispersed species was also higher in logged forest. However, our results indicate that logged forest plots were more similar in species composition with high evenness, indicating that a few species dominated the community. For the biotically-dispersed guild, logged forest was much less diverse, with a low species overlap with unlogged forest. In addition, logged forest plots were again more similar in species composition, and were dominated by a few species.
Species richness was also markedly lower in logged forest, the difference being so for both functional guilds. Within biotically-dispersed trees, species that were bird-dispersed showed a significant decrease in species richness in logged forests. However, species-area curves were not asymptotic overall or for biotically dispersed species, although estimated species richness was found to be higher in intact forest. Previous studies in the area have shown that populations of large-bodied avian dispersers are lower in logged areas [15] and this may be affecting recruitment [17]. Larger vegetation plots need to be set up on a priority basis to capture the full complement of adult trees and to understand finer scale differences for rare species (especially those dispersed by mammals and by both birds and mammals).
Structurally, our community-wide analysis on tree species showed that biotically-dispersed species occurred at lower densities. Other studies have found that stem density increased > 5 years post-logging, probably because greater canopy openness promoted the establishment and recruitment of pioneers and other competitive colonizers [26]. There was no significant difference in density of biotically-dispersed tree species between the two regimes, but there was a notable increase in abiotically-dispersed tree species. Also consistent with previous studies [9,12,19], size of trees (DBH) was significantly smaller in logged areas, unsurprising since the larger size classes are usually preferentially logged and it may take many decades for slow-growing tropical trees to reach the pre-logging size-classes.
The decrease in densities of biotically-dispersed tree species could alter the spatio-temporal availability of fruits for frugivores. In addition, many frugivores either decline in abundance or disappear from logged patches due loss of habitat or changes in the abundance and quality of food sources [12,13,18,27]. Frugivores are important seed dispersers and affect the spatial patterns of plant recruitment [2,4]. The loss of frugivorous dispersers can alter tree community composition [29], as is well documented in hunted forests [30]. However, few studies have tested for finer scale responses to selective logging in dispersal and recruitment among trees with different dispersal modes, especially comparing bird and mammal dispersed species, although studies of individual species indicate that such differences are expected [31].
Implications for conservation
The moist forests of north-east India continue to suffer high pressures from logging, hunting, and large-scale felling for agriculture. Remaining forests, especially protected areas, serve as vital population sources of biodiversity, especially plant species richness [32]. In addition to the stark loss of forests and wildlife, changes in forest composition and subtle losses of ecosystem function, as seen in our study, are likely to impact the richness of these forests in the long-run. Our study site is also representative of the conservation challenges faced by the global hotspot of the Indo-Myanmar area in dealing with timber over-harvesting [33,34]. In other parts of south-east Asia, selectively-logged and human-modified forests provide important ecosystem services and allow persistence of a wide variety of faunal species [35,36,37]. However, we find that although logged forests had marginally lower tree density overall, they had significantly lower tree species richness and diversity. Our study clearly indicates that selective logging significantly alters the functional composition of forests. Even though 16 years is a relatively short recovery period for slow-growing tropical trees, the differences in density, richness, and diversity of biotically-dispersed tree species between selectively logged and intact forests are of importance for conservation. Even amongst abiotically-dispersed tree species the community composition was more homogenous in logged areas.
The effects of dramatic habitat loss are usually clear, and populations and communities of plants and animals continue to be affected by human modification of landscapes [38]. The more insidious effects of habitat degradation, however, might be just as severe in terms of loss of ecosystem function [38,39]. Lack of post-harvest assessment presents a major impediment to monitoring species recovery in logged forests, especially with regard to functional guilds. Together with resident communities, the Arunachal forest department has instituted a variety of co-management schemes which include creation of village forest reserves and community conservation areas in many parts of the state. Given that these areas are still open to timber and firewood extraction, which directly impact local livelihoods, our study is of particular relevance and indicates that the present timber extraction patterns may be unsustainable. Additional information is required to monitor seedling recruitment and species regeneration. Further, the effects of logging on abiotic factors such as light, moisture availability, and soil characteristics and its effect on community dynamics of regenerating seedlings may also provide important insights on regenerative ability of species. Detailed information about such processes will allow us to better predict the impacts of specific management practices on tree species and communities, and scope for recovery following logging.
Acknowledgements
We obtained funding for this project from the Department of Biotechnology grant (number 6324) provided to Mahesh Sankaran at the National Centre for Biological Sciences. The Arunachal Forest Department granted research permits. We thank Putul Munda, Arun Saikia, Kishore Dorje for their help in the field. We thank Mahesh Sankaran, Jayashree Ratnam, Aparajita Datta, David Edwards, Umesh Srinivasan, Dayani Chakravarthy, Shashank Dalvi and Sachin Sridhara for ideas, discussions and help with analysis. Dr. Haridasan K identified plant specimens. George Gale and Alejandro Estrada helped improve our manuscript. For honouring our collective obligation to save Pakke Tiger Reserve, this paper is dedicated to Tana Tapi and his team of forest watchers.