Dissimilar effects of human and elephant disturbance on woodland structure and functional bird diversity in the mopane woodlands of Zambia

Humans and elephants are major disturbance agents in the African savanna woodlands. While both species remove trees, humans selectively harvest larger stems, which are less vulnerable to elephants. Increasing human pressures raise the question of how the altered disturbance regime will modify woodland structure, and in turn biodiversity and ecosystem function. Here we investigate this process in the mopane woodlands of Zambia by examining relationships between woodland structure, species and functional bird diversity, and human and elephant disturbance intensity. We conducted a single-season comparison of 178 plots from 45 sites using Bayesian mixed models. The effect of elephants on tree density (− 7.7 ± 1.6%; deviation from intercept) and bird species richness (− 15 ± 6%) was greater than that of humans (density: − 3.5 ± 1.5%; bird richness: − 11.6 ± 4.7%). Despite this, elephants did not significantly affect woody biomass or functional bird diversity, whereas humans had a negative effect on both (woody biomass: − 9.3 ± 2.3%; functional bird diversity: − 5 ± 2%). Elephants were associated with reductions in species and functional turnover (5.3 ± 2.5% and 6 ± 3%, respectively). Replacement of elephants by humans is likely to reduce woody biomass and functional bird diversity affecting the woodland structure, sustainability, and functioning. Concentrated elephant disturbance could lead to spatial homogeneity in bird species and functional compositions, potentially reducing the spatial resilience of bird communities. This is the first study to highlight how the difference between elephant and human disturbances leads to dissimilar effects on biodiversity.


Introduction
Humans and elephants are integral and ancient parts of the African savanna woodlands, structuring the ecosystem through disturbance (Ellis 2011;Charles-Dominique et al. 2016). However, in recent times their relative influence as disturbance agents has radically changed (Boer et al. 2013); humans have replaced elephants as the dominant disturbance agents with unknown consequences for biodiversity and ecosystem function. As elephants and humans both remove trees (Mograbi et al. 2017), they may appear functionally interchangeable (Hempson et al. 2017), but there are important differences in the way they use savannas and the size of the trees they preferentially remove (Staver and Bond 2014;Woollen et al. 2016;Mograbi et al. 2017). As such, elephants and humans may be modifying woodland structure and thus its biodiversity in dissimilar ways, but the identification of such discrepancies-and the assessment of their effects-has received surprisingly little attention.
Here we present a spatial comparison of the effects of disturbance by humans and elephants on (i) structural attributes of the woody vegetation (ii) species and functional diversity of birds in the mopane woodlands of Zambia.
Being dominant land herbivores and bulk feeders, elephants alter the woodland structure by removing trees for forage through pollarding, uprooting, and debarking (O'Connor et al. 2007). They target trees predominantly in the 5-9 m height range-the elephant browse-trap (Asner and Levick 2012;Staver and Bond 2014)-and, through pushing over and pollarding, they keep patches of trees in suppressed, stunted forms (Styles and Skinner 2000). Elephants repeatedly visit the stunted trees for the nutritious leaves that resprout in response to the extensive browsing (Kohi et al. 2011). However, certain treesdepending on their growth rates and the densities of elephants-episodically escape the elephants' 'browse trap' and grow large enough to become less vulnerable to elephant disturbance (Shannon et al. 2011;Staver and Bond 2014). These large trees possibly counter-balance the woody biomass lost due to elephant-mediated tree removal and suppression at the landscape-scale (Ben-Shahar 1996). These processes maintain the open-canopied structure of savanna woodlands and prevent woody encroachment . Furthermore, the elephantmaintained open woodlands containing large trees are associated with woodland-specific co-occurring species like the cavity-nesting and ground-foraging Lilian's Lovebird Agapornis lilianae (Mzumara et al. 2014), hence, elephants also facilitate the distinct species and functional diversity of the savanna woodlands (Botes et al. 2006;Pringle 2008).
In contrast to elephants, humans selectively harvest-consume (Archibald 2016)-large trees for fuelwood, construction and timber, especially those that produce good quality, slow-burning commercial charcoal (Woollen et al. 2016). The selective removal of large trees means that woody biomass in the humanutilized woodlands declines at a much faster rate than tree density (Woollen et al. 2016). As a result, the woodlands may transform into novel ecosystems composed of stems of smaller sizes, with increased shrub vegetation in the areas where large trees formerly persisted (Hosier 1993;Kalema and Witkowski 2012). Furthermore, once the large trees are depleted, the smaller, less desirable trees are also harvested to meet the increasing demands for fuelwood, further reducing the remaining woody biomass (Hosier 1993). There is therefore potentially no escape for trees from the human trap. The human-utilized woodlands devoid of large trees and increased ground 123 vegetation are known to cause a switch to nonwoodland and shrubby habitat-preferring faunal communities (Botes et al. 2006;Seymour and Dean 2010).
Increased hunting and anthropogenic land-use change have reduced elephant population sizes and compressed their distributions into a small network of protected areas-a fraction of their historical ranges (Robson et al. 2017). Humans have therefore replaced elephants as one of the dominant disturbance agents in African savanna woodlands (Hempson et al. 2017). To understand the consequences of this replacement, a comparison of the effects of human and elephant disturbance on woodland structure and biodiversity is required. However, to our knowledge, there are only two human-elephant effect comparisons. Botes et al. (2006) showed that dung beetle assemblages in elephant disturbed forests were more similar to undisturbed forests (73% similar) than the assemblages in human-disturbed forests were (50% similar). Mograbi et al. (2017) compared biennial tree fall rates and demonstrated that elephant and human disturbances affect all woody height classes with both inducing shrubland conversion. While Botes et al. (2006) did not investigate functional diversity patterns, Mograbi et al. (2017) is restricted to a single attribute of the woodland structure-the tree density. Comparisons incorporating multiple attributes of the woodland structure-woody biomass, tree size, and ground cover, and biodiversity-species and functional a-and b-diversity are required in order to fully assess potential differences between human-and elephant-mediated change. Separate studies of human and elephant disturbance suggest that while anthropogenic impacts on woody biomass and biodiversity are often negative (but see Luck 2007;Winfree et al. 2007;Andela et al. 2017 for exceptions), elephant impacts are variable, and dependent upon elephant density as well as factors such as fencing, proximity to artificial water bodies, rainfall and soil character (Guldemond and Van Aarde 2008;Guldemond et al. 2017). Furthermore, studies have demonstrated that increased elephant concentrations may reduce landscape heterogeneity through size-structured tree removal (Asner and Levick 2012), with possible effects on the landscape-wide species and functional b-diversity of faunal communities (Herremans 1995;Cumming et al. 1997). Thus, human and elephant disturbance may be associated with divergent consequences, an understanding of which is paramount for the effective management of the savanna woodlands.
We hypothesized that elephant and human disturbance would reduce tree density but have different effects on other attributes-woody biomass, tree size, and ground cover-of the woodland habitat, and consequently bring about dissimilar impacts on species and functional diversity of birds. We tested this hypothesis in the mesic mopane woodlands of Zambia where elephant-dominated areas (unfenced national parks and game reserves), human-utilized lands (non-protected areas near urban centres), and low-disturbance regions (remote non-protected areas) are spatially exclusive. We chose birds as the taxonomic group for this study because they are relatively visible and easy to identify (Bibby 1999), are sensitive to habitat structural modifications (Owens and Bennett 2000) and are good indicators of effects of habitat filtering process (Benítez-López et al. 2017).

Study area
The study area spans the distributional range of mopane woodlands in Zambia -identified using White's African vegetation map (White 1983)which includes the eastern part of the country from north to south, covering the Luangwa, Luano and Zambezi valleys (Fig. 1). The area includes three national parks (North Luangwa, South Luangwa, and Lower Zambezi), two major charcoal and timber production regions (Chirundu and Chipata; (Gumbo et al. 2013), and remote areas of low human-utilization (Luano). The Zambian mopane region is mainly occupied by alkaline alluvial soils, and receives mean annual rainfall of 750-1000 mm in central-north and 500-750 mm in the south (Nkhuwa et al. 2016). Monospecific stands of Colophospermum mopane dominate the region, along with occasionally cooccurring species such as Combretum apiculatum, C. zeyheri, Terminalia prunioides and T. sericea (Martini et al. 2016).
At the beginning of the 2015 dry season (May-July) we measured woodland vegetation structure and disturbance in 178 plots clustered within 45 sites (grids of 1 km 2 ). At the centre of 120 plots from 30 sites, we collected bird diversity data using pointcounts. Each site (except one for the woodland structure, where we could measure only 2 plots) contained 4 circular plots 500 m apart. The sites were distributed across the mopane landscape from south (Sinazeze; -17.46732, 27.29626) to north (North Luangwa;-11.91523, 32.40934), were at least 250 m away from any tar road, large river or waterbody, were at least 1 km away from each other, and had a mean altitude below 750 m.

Woodland structure
In each circular plot, for all tree stems C 10 cm DBH (diameter at breast height, 1.3 m) we identified species, measured DBH and counted the total number of stems where 2 cm \ DBH \ 10 cm. To estimate ground vegetation cover (hereafter ground cover), we used four smaller quadrats of 1 m 2 placed within the circular plot, one in each cardinal direction, and at 10 m from the centre. Each ground cover quadrat was subdivided into 4 blocks such that in each circular plot, the quadrats together constituted a total of 16 blocks. The number of blocks occupied by ground vegetation was used as an index of ground cover. For unidentified trees, photographs and sample specimens were collected for identifications with help of botanists in Lusaka and online floras (Bingham et al. 2017;Hyde et al. 2017).
We computed five variables to represent woodland structure 1. Density (N, total stems ha -1 ), the number of tree stems C 10 cm DBH. 2. Small stem density (N s , stems ha -1 ), the number of tree stems where 2 \ DBH \ 10 cm. 3. Woody biomass (B, tC ha -1 ), calculated following the allometric Eq. (0.0267 9 DBH 2.59 ) in Ryan et al. (2011). 4. Tree size (D, cm), mean DBH of all tree stems. 5. Ground cover (C), the proportion of the total number of blocks covered by ground vegetation. top right) contained 178 circular sampling plots 500 m apart. Point counts for birds were undertaken in the centre of plots (bottom right). Within each plot, four 1 m 2 quadrats divided into 16 blocks (bottom right) were used for ground vegetation surveys. The inset shows regional context-mopane woodlands (light grey) in Zambia (black) and southern Africa

123
Bird species and functional diversity variables Bird data were collected at the centre of each circular plot between 6 and 9 am by the point-sampling presence-absence survey method (Gregory et al. 2004), where the bird species seen or heard during the sampling period of 15 min were listed. For each recorded species, we collected trait information from Elton Traits 1.0 (Wilman et al. 2014). The trait values were the relative percentages of 10 diet and seven foraging strata categories, activity time (categorical, 1 or 0) and body mass (g)-table of traits provided in supplementary information (SI; Table S1). For four unidentified species, mean values were used (the analysis was also conducted after removing the unidentified species; this did not make any significant difference, so only results using the mean values are presented).
We computed aand b-diversities for species and functional type from the survey data. Since the main objective of this study was to examine changes in composition due to species or functional replacement, we used only the turnover component of the multicomponent (nestedness and turnover) b-diversity (Baselga 2010). Species a-diversity (species richness, S) is the total number of species recorded in each plot and is conditional upon the assumption that the bias due to imperfect detection would be uniform (since sampling effort was constant). Species b-diversity was derived by computing the incidence-based mean pairwise Simpson dissimilarity in species-presence between communities (S bsim ), following Baselga (2010).
For functional diversity, we first made a Bray-Curtis pairwise functional distance matrix of species. We then computed the mean functional pairwisedistance (FD) between species in a community as a measure of functional a-diversity, following Swenson (2014). Further, using a functional dendrogram (Figure S2 in SI), derived from the functional distance matrix, and the presence-absence matrix of bird species in each community, we calculated community-level pairwise Simpson dissimilarity in functional space as a measure of functional b-diversity (FD bsim ), following Leprieur et al. (2012) and Melo et al. (2014).

Human and elephant disturbance
We counted the number of stems impacted by elephants, evidenced by knocked down, pollarded, broken, or bark-stripped stems, and affected by humans, evidenced by cut stems. As an indicator of disturbance, we computed the % of the total number of stems in a plot that were affected by elephants (I Ele ) and humans (I Human ). We also derived mean pairwise geographical distance (g) between plots to account for spatial autocorrelation.

Statistical analyses
For descriptive analysis of differential impacts and size-class specific variances, we classified the plots into three categories: Low disturbance (n = 85, (I Ele and I Human \ 20%), Human-dominated (n = 51, 100 C I Human C 20% and 0 B I Ele \ 20%), and Elephant-dominated (n = 40, 100 C I Ele C 20% and 0 B I Human \ 20%). The number of plots where both disturbances occurred were few (n = 2, I Ele and I Human [ 20%) and hence were excluded.
The human population density of the study region has doubled since the 1970s (Ohadike and Tesfaghiorghis 1975;Hoare and Toit 1999) whilst the population of elephants has reduced 10-fold (Kampamba et al. 2003;Chomba et al. 2012). As human population increased, elephant densities reduced and their distribution became restricted to the few protected areas (Boer et al. 2013). In this study, therefore, the elephant-dominated areas-mainly national parks, represent the early stage in this transition to human domination, whilst non-protected areas near towns represent the final stage of transition. Areas with low activity of both elephants and humans-principally game management areas-are considered ''transitional'' (Table 1).
To examine the effects of disturbance on woodland structure and bird diversity, we constructed mixed models. We used each of the woodland structure (Eq. 1) and bird diversity (Eq. 2) variables, respectively, as responses, disturbance variables and geographical distance as main effects, and sites within the plots as nested random effects.
where, W is the woodland structure variable {N, N s , B, D, and C}, D is the diversity variable {S, S bsim , FD, FD bsim } for plot i {1..0.178}, a and b are parameters to be estimated from the model, a j is random effect of site j {1..0.45} within which the plots are clustered, and a j $ Normal l a r 2 a À Á . To explain the community-level patterns of bird diversity better, we also investigated the effects of disturbance on species presence (Eq. 3).
where Sp is the species presence {species 1…132}. For all woodland structure and bird diversity variables (except species richness), we fitted a linear mixed model with normal distribution. We assessed the habitat and diversity variables for normality, and subsequently log 10 -(hereafter log) transformed the non-normal variables, N, N s , B, and D. For bird species richness and incidence, we fitted Poisson and Bernoulli random-effect models, respectively.
We also ran the above models including environmental variables (Relative elevation, Mean Annual Precipitation-MAP, Mean Annual Precipitation-MAT, and Soil-soil type cluster) as linear predictors.
We analysed models using Bayesian methods (Kéry and Royle 2016) implemented in WinBUGS. We provided non-informative flat priors and used three Markov Chain Monte Carlo (MCMC) chains, with 75,000 total iterations each. The first 25,000 iterations were discarded as burn-into remove transient states and attain equilibrium distribution. From the remaining 50,000 total iterations, we selected samples after every 50 iterations (thinning), thus yielding 3000 samples from all three chains for posterior inference. We checked model convergence using the Gelman-Rubin convergence diagnostic (Gelman and Rubin 1992), with potential scale reduction factor values approaching 1 considered acceptable (Kéry and Royle 2016). We used the 95% Bayesian Credible Interval (CI) to indicate significant effects and % deviation from the intercept as the standard effect size (standardised coefficients ± standard error).
Analyses were conducted using R statistical software version 3.4.2 (R Core Team 2017). We used the vegan package (Oksanen et al. 2016

123
The stem-size distribution in this study followed the typical reverse-j shaped, characterized by the dominance of stems in smaller size classes followed by decreasing number of stems of larger size classes. Stemsize distribution differed between human-and elephantdominated areas. In comparison to the low disturbance plots, the human-dominated plots contained fewer number of large-size stems (Fig. 2a), whereas, the Elephant-dominated plots comprised a relatively longer tail of stems in large-size classes (Fig. 2b).
As expected, both I Human and I Ele were negatively associated with stem density, but they affected woody biomass, mean DBH, and small-stem density differently (Fig. 3, coefficient plots in SI- Figure S3), resulting in different stand structures (Fig. 4a-c). Although the effect of I Ele (-7.7 ± 1.6%) on stem density was twice that of I Human (-3.5 ± 1.4%), I Ele did not significantly affect the woody biomass, whereas I Human had a negative effect (-9.3 ± 2.7%). Further, mean DBH was associated positively with I Ele (5.4 ± 1.2%) and negatively with I Human (-1.8 ± 1.1%; only significant at 90% CI). The small-stem density significantly increased with I Human (15.7 ± 7.5%) but did not change with I Ele . Ground cover increased with mean geographic distance between plots and did not show any significant relationship with I Human or I Ele .
I Human and I Ele were both associated with declining bird species richness but affected the other bird diversity variables differently ( Fig. 5; coefficient plots in SI- Figure S4). I Ele had a stronger effect on species richness (-15.5 ± 6.2%) than I Human (-11.7 ± 4.7%) and was associated with a reduction in S bsim (-5.4 ± 2.5%), while I Human showed no effect.
Despite the negative effects on aand b-diversities, I Ele did not significantly affect FD, although FD bsim did reduce with increasing I Ele (-4.9 ± 2.8%). In contrast, I Human was associated with reductions in FD (-4.3 ± 1.8%) but had no effect on FD bsim .
The bird species incidence models suggested that more species showed a significantly negative association with I Human (n = 31) than I Ele (n = 18). The species negatively affected by I Ele were primarily insectivorous habitat generalists (Terpsiphone viridis and Tchagra australis), non-mopane species (Prionops retzii, Merops pusillus) and species associated with regrowing woodlands (Camaroptera brachyura and Prinia subflava). Those that increased with I Ele were mainly woodland-preferring seed eaters and ground foragers (mopane specialists Agapornis lilianae and Emberiza flaviventris), woodland-specific insectivores (Lamprotornis mevesii and Nectarinia amethystina), and plant-eating habitat generalists (Estrilda astrild).
Inclusion of environment variables did not make any difference to our results; in other words, environmental variables are not contributing to the differences we detected.

Discussion
Our study is the first to demonstrate that the replacement of elephants with humans as a major agent of Fig. 3 The effect of disturbance regime on woodland structure in 20 m radius plots. Lines denote predicted estimates, with solid lines indicating significant, and dashed lines non-significant, effects (within 95% CI). The shaded areas are the 95% CIs. Colours indicate type of the impact-elephant (blue) and human (orange). Both human and elephant disturbance were associated with reducing stem density but had dissimilar effects on other attributes of the woodland structure Fig. 4 Examples of mopane woodland under different disturbance regimes. The mopane woodlands with relatively low disturbance (a) have relatively high stem density and a mixture of stems of different size classes. In the more human-affected woodlands (b), large trees are selectively harvested, leading to reduced density, loss of woody biomass, and an increase in small-stem density. Elephant-dominated woodlands (c) contain smaller stems which are frequently browsed by elephants, and occasional stems of larger size classes which have escaped elephant impact. These large stems maintain the plot-level woody biomass 123 disturbance can alter the structure of savannas and have negative effects on bird biodiversity and functionality.

Effects on woodland structure
As expected, both humans and elephants were associated with reduced tree densities. The treefall rates reported in other studies (Asner and Levick 2012;Mograbi et al. 2017) are consistent with the patterns observed here, and are 7.7 ± 1.6% ha -1 by elephants and 3.5 ± 1.4% ha -1 by humans. This study provides mopane woodland specific support to the findings of Mograbi et al., (2017) in Kruger National Park that the magnitude of tree removal by elephants is generally higher than that of humans.
Despite a greater negative effect on tree density, elephant disturbance did not significantly reduce woody biomass, since elephant disturbance was associated with an increase in mean tree size. The relatively more frequent presence of larger trees ([ 35 cm DBH or 11 m height) in the elephantdisturbed areas counterbalanced the lost woody biomass as elephants damaged smaller trees. Elephants are known to predominantly utilize the smaller stems for browse and selectively revisit and reuse the previously utilized woodland patches which produce more nutritious and easily accessible browse. Certain trees, however, episodically escape from the elephant browse-trap (Asner and Levick 2012; Staver and Bond 2014) becoming larger and less vulnerable to future elephant impacts. This phenomenon of the escape from the elephant impacts may explain the higher number of large trees in elephant-impacted areas compared to those of humans, evident from higher mean tree size and long tail of large-size stem classes. However, as elephant density increases, the chance of trees escaping the browse-trap will diminish, with a concomitant loss of woody biomass (Ben-Shahar 1996). This density-dependent impact of elephants is influenced by rainfall, which mitigates the damage by promoting higher woody cover (Hempson et al. 2015), and by fences, which aggravate it by confining elephants and thereby concentrating their activities (Guldemond and Van Aarde 2008). Our study areas were not fenced, and since we sampled north-south regions of the mopane distribution in Zambia, we covered a representative range of rainfall. The observations made during this study are therefore considered to be generalizable, at least across Zambia and Fig. 5 The effect of disturbance agent on species and functional bird diversity. Lines denote predicted estimates with solid lines indicating significant and dashed lines non-significant effects (within 95% CI). Shaded area shows the 95% CI. Both human (red line) and elephant (blue) disturbances were associated reductions in species richness, but only human disturbances were associated with significant declines in functional bird diversity regions with similar climatic and edaphic conditions. In addition, our finding that elephant disturbance has no significant impact on woody biomass supports the consensus in the literature that, at moderate densities at least, elephants have little impact on woody biomass (Ben-Shahar 1996; Guldemond and Van Aarde 2008; Fig. 6 The effect of human and elephant disturbance on species presence. The dots represent scaled model coefficients (proportion of deviance from the intercept), horizontal lines indicate 95% CI, and colours denote direction of the significant (95% CI) relationship. For both disturbance there are more losses than gains in species incidence probability. Only species with at least one significant effect (n = 60) are shown 123 Kalwu et al. 2010). We acknowledge that what constitutes 'moderate' density is debatable, and that the definition varies between studies (Ben-Shahar 1996;Holdo et al. 2009;Kalwu et al. 2010); further research would be needed to draw objective thresholds.
While elephants mainly affect smaller trees, human use of the woodland is primarily driven by the need for timber and fuelwood, especially for charcoal production, leading to upsurge in selective logging of large stems (Hosier 1993;Ndegwa et al. 2016;Woollen et al. 2016). Therefore, in human-dominated landscapes, large trees face high selective pressure and rarely escape the human impact (Woollen et al. 2016;Mograbi et al. 2017). Hence, at the landscape level, human-affected areas are likely to undergo a decline in number of large stems and woody biomass, leaving behind small undesired or coppicing stems (Backéus et al. 2006;Fontodji et al. 2011;Ndegwa et al. 2016;Kiruki et al. 2017). As humans replace elephants as dominant agents of disturbance, the mopane woodlands are, therefore, more likely to become lowbiomass degraded systems, degraded from the perspective of carbon storage (McNicol et al. 2018).

Effects on bird diversity
Human and elephant disturbance both caused reductions in bird species richness-likely by removing trees that birds exploit for food and nesting (Skowno and Bond 2003)-but had dissimilar effects on other bird diversity attributes; functional diversity only reduced in response to human disturbance, whereas the species and functional b-diversities decreased with elephant disturbance.
The open-canopied mopane woodlands provide a characteristic habitat for cavity nesting, ground foraging, and canopy-specialist birds-the woodland species (Herremans 1995). A less disturbed woodland containing trees of diverse size classes has high habitat heterogeneity and supports bird communities composed of the dominant woodland species as well as the locally infrequent non-woodland species (Herremans 1995). While, the dominant woodland species form the core of the bird community, the infrequent nonwoodland species largely determine the compositional variability-the spatial b-diversity (Socolar et al. 2016). Although both human and elephant disturbances were associated with reductions in bird species richness, the type of species that were affected differed between human and elephant disturbances. Elephant disturbance was positively associated with the woodland bird species-functionally unique to the woodland bird communities, whereas human disturbance was positively associated with the infrequent nonwoodland species.
As the elephant-dominated areas maintained the dominant-and functionally unique-species of the community and lost only the infrequent bird species, the bird communities underwent spatial taxonomic and functional homogenization in response to elephant disturbance, i.e. a reduction in the turnover component of species and functional b-diversity (Socolar et al. 2016). Despite a reduction in functional turnover, there was no significant difference in functional alpha diversity, demonstrating high functional redundancy of bird communities in elephant-disturbed areas. In contrast, in the human-dominated areas, where the woodland structure is transformed into shrubby system, the dominant woodland bird species are reduced, and the infrequent non-woodland species are retained. This may explain the reduction in bird species richness and functional diversity without any significant effects on b-diversity. In circumstances where the woodlands are completely transformed into novel systems (e.g. croplands), the majority of the dominant woodland bird species may be replaced by the infrequent nonwoodland bird species causing significant increase in b-diversity and compositional drift of bird communities (Socolar et al. 2016).
Our finding that the human-led habitat modification reduces species and functional diversity is consistent with other studies (Flynn et al. 2009;Luck et al. 2013). However, the finding of dissimilar effects of human and elephant mediated habitat modification on species and functional diversity is novel. The results here indicate that despite losing some species, elephant disturbed woodlands maintain functional bird diversity. However, given that spatial functional b-diversity is reduced, the functional resilience of bird communities in elephant disturbed areas may be compromised. Importantly, elephant-dominated landscapes contain species of high conservation value like Agapornis lilianae (endemic to the mopane woodlands of the Zambezi basin and near threatened) and woodland specialist bird species such as Emberiza flaviventris, Parus niger and Thripias namaquus.
These findings imply that to manage the negative effects of human disturbance on woodland birds, large trees (C 0.35 m DBH or C 11 m height-the mean values in elephant-dominated areas) should be retained in the landscape. In our study area, this might be achieved by altering the charcoal and timber license conditions issued by the Forest Department in Zambia. Finally, this study is a snapshot of the complex dynamics of human and elephant interactions in the African woodlands. Although it only examined spatial variability, it provides an indication of how woodland structure and bird diversity may respond to changing disturbance regimes over time. More studies with spatial and temporal replications are required to better understand the dynamics of human and elephant mediated changes in woodland structure and biodiversity.
In summary, this study concludes the following: (i) Although both modes of disturbance are associated with reduction in tree density, elephant disturbance does not necessarily reduce woody biomass, due to an increase in mean tree diameter. This increase is driven by the size-structured episodic escape of certain trees from the elephant browse trap (Asner and Levick 2012;Staver and Bond 2014). The escaped trees grow larger, and consequently become less vulnerable to elephant disturbance. In the densities at which they occur in the Zambian mopane region (0.44 km -2 ; DNPW 2016), elephants currently do not affect large tree populations substantially. (ii) In the human-utilized areas, on the other hand, humans selectively harvest larger trees for fuelwood (especially charcoal production) and timber, thereby reducing mean tree size and woody biomass at the same time as tree density and promoting an increase in the density of small stems. (iii) The dissimilar effects of human and elephant disturbance on woodland structure are apparent in patterns of bird diversity. Although both human and elephant disturbance were associated with reductions in bird species richness, only human disturbance was associated with reductions in functional bird diversity. The woodland-specific dominant bird functional groups (ground foraging plant and insectivores) which were not affected by elephant disturbance showed a negative response to human disturbance. In other words, elephant disturbance is associated with declines in non-woodland non-ubiquitous species which contribute very little to the functional diversity of the woodland bird communities, whereas human disturbance is associated with declines in ubiquitous species of the woodland resulting in functional decay. (iv) Although elephant disturbance does not reduce functional bird diversity, it does decrease species as well as functional bdiversity, possibly via a reduction in landscape-scale habitat heterogeneity.
Acknowledgements The authors thank the two anonymous reviewers for their critical and constructive comments, which helped to improve the manuscript. Fieldwork of this study was funded by World Parrot Trust. Fieldwork permissions in national parks was facilitated by the BirdWatch Zambia and assistance was provided by the Zambia Wildlife Authority. We thank Mike Bingham, Lusaka, for help in plant identifications. We acknowledge support from the Abrupt Changes in Ecosystem Services project (NE/K010395/1) which was funded with support from the Ecosystem Services for Poverty Alleviation (ESPA) programme. The ESPA programme is funded by the Department for International Development (DFID), the Economic and Social Research Council (ESRC) and the Natural Environment Research Council (NERC). We also thank Dr Andrew Dobson and Dr Rose Pritchard for proof reading this manuscript and providing comments.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Author contributions HT, CR, TG, RM and CLP conceived the ideas and designed methodology; HT, TG and CP collected the data; HT analysed the data; HT led the writing of the manuscript. All authors contributed critically to the drafts and gave final approval for publication. The authors declare no conflict of interest.
Data accessibility Data and R codes used in this study are available from the Git Repository. https://bitbucket.org/ed_ trop_land_use/2018_birdhab_zambia/src.