Bird communities in a degraded forest biodiversity hotspot of East Africa

Tropical forests suffer severe habitat destruction. Thus, tropical forests frequently consist today of only a few small remnants that are often embedded within a matrix of agricultural fields and tree plantations. Forest specialist species have experienced severe population declines under these circumstances. We studied bird communities based on census plots set up in a near-natural forest block, as well as degraded forest patches, tree plantations, and agricultural fields, across the Taita Hills in southern Kenya. We classified each bird species according its ecology and behavior. We quantified the land cover and landscape configuration around each census plot. Typical forest species were mainly observed in the near-natural forest block, and to a lower extent in degraded forest patches. Plantations were almost devoid of birds. Bird communities of small forest fragments were more similar to that of agricultural land than the near-natural forest block. Most frugivorous, insectivorous and nectarivorous birds occurred in forest habitats, while granivorous bird species dominated the bird communities of agricultural land. The surrounding landscape had a marginal impact on bird species composition at local sites. Our study showed that the preservation of near-natural cloud forest, including small forest patches, is essential for the conservation of forest-dependent species, and that plantations do not serve as surrogate habitats.


Introduction
A major driver leading to current biodiversity loss is the destruction of natural habitats and their transformations into agricultural land, settlements and plantations (Jantz et al. 2015;Maxwell et al. 2016). Tropical ecosystems are particularly affected by these threats, which are mainly driven by high human demographic pressure, weak governance structures and poor land management (Williams 2013;Nzau et al. 2020;Teucher et al. 2020). As a consequence, many remaining tropical ecosystems exist today as small and isolated remnants (Haddad et al. 2015) which provide comparatively low habitat quality (Sala et al. 2000;Wilkie et al. 2011). The destruction of habitats and the deterioration of habitat quality mostly result in a decrease of population size of species, and, as a consequence, these local populations often suffer population fluctuations, reduced population viability and low individual fitness (Reed and Frankham 2003;Melbourne and Hastings 2008).
Habitat destruction and the transformation of natural ecosystems into anthropogenic landscapes create surrogate habitats for many species (Greenler and Ebersole 2015). For example, clearing of natural forests and their transformation into plantations of exotic trees may provide surrogate habitats for species, including forest-dependent taxa (Berndt et al. 2008;Calviño-Cancela 2013). Or, the transformation of natural forest into settlements frequently includes the creation of gardens consisting of hedges, trees and flowers, and thus also provide suitable habitats for many species (Batáry et al. 2010;Heath et al. 2017;Wilson et al. 2017). In addition, the fragmentation of a formerly interconnected forest creates transition zones between forest patches and non-forest ecosystems that also provide valuable habitats for many species (Morante-Filho et al. 2018). However, such transformed and disturbed habitats are usually of lower habitat value for species depending on specific habitat conditions, and may even constitute ecological traps for some species (Barlow et al. 2007). Apart from habitat availability and habitat quality, the surrounding landscape may strongly impact the viability of populations, and thus shape species communities. Studies have shown that neighbouring ecosystems can have a greater effect on the local population than the quality of the corresponding habitat (Metzger 1997).
Species respond differently to environmental changes, such as habitat destruction and the deterioration of habitat quality (Bryce et al. 2002). Habitat specialists, such as species adapted to live in a dense and intact cloud forest, mostly occur exclusively in a specific habitat with specific resources. In consequence, specialist species are frequently restricted to still intact natural habitats and rarely use surrogate habitats, and do not live in small and isolated habitat fragments. Thus, species communities found across anthropogenic landscapes are frequently dominated by some few habitat generalist species (Kitahara and Fuji 1994;Abrahamczyk et al. 2008;Börschig et al. 2013).
In this study, we performed census counts on birds in Taita Hills in southern Kenya. The Taita Hills are part of the Eastern Afromontane biodiversity hotspot and host a unique flora and fauna (Mittermeier et al. 2011). The transformation of cloud forest into agricultural land and plantations of exotic trees led to an almost complete destruction of the original forest Lehouck et al. 2009). Today, the landscape of Taita Hills consists of a mosaic of settlements, agricultural fields, tree plantations and cloud forest patches. We analyzed bird communities across a fraction of the Taita Hills, including a near-natural forest, small forest fragments, plantations of exotic trees, and smallholder agricultural fields. We grouped bird species into guilds, according its ecology and behavior. We considered the land cover and landscape configuration around each census plot. Based on the results derived from our census counts we addressed the following questions: 1. Do bird communities differ among habitat types, with respect to species and guild composition? 2. Do small and degraded forest fragments and plantations of exotic trees provide suitable surrogate habitats for forest-dependent birds? 3. Does the landscape around a local census plot affect bird species community structure and species diversity?

Study area
The Taita Hills in southern Kenya (3° 25′ S, 38° 20′ E) consists of a mosaic of ecosystems, with near-natural forest blocks, small and degraded forest fragments, plantations of exotic trees (with representatives of the genus Pinus, Cypressus and Eucalyptus), agricultural land (crop fields with gardens), and settlements (Capitani et al. 2019;Teucher et al. 2020). The cloud forest cover of Taita Hills had been widely transformed into anthropogenic landscapes since the pre-colonial era (about 200 years ago), and continues until today, despite intense effort to protect the natural vegetation Lehouck et al. 2009;Teucher et al. 2020). The remaining natural cloud forest consists of 12 forest fragments, which differ in size (ranging from 1 to 200 ha), geographic isolation, and habitat quality (Lehouck et al. 2009;Callens et al. 2011). In our study, we selected the following four habitat types: 1. Near-natural cloud forest block (here Chawia forest), 2. Small and degraded forest fragments (here Fururu and Ndiwenyi forest patches), 3. Plantations of exotic trees, and 4. Agricultural fields. Chawia represents the third largest natural cloud forest in Taita Hills (Pellikka et al. 2009) with still largely intact plant species composition with Tabernaemontana stapfiana, Albizia gummifera, Phoenix reclinata and Macaranga conglomerata being the dominant tree species (Omoro et al. 2011). However, planting of exotic tree species across this forest patch has diminished its habitat quality Omoro et al. 2011). Fururu and Ndiwenyi represent two small and degraded forest fragments, which are interspersed with exotic trees such as eucalyptus, cypress and pine (Omoro et al. 2010). These forests are much more disturbed compared to the near-natural Chawia forest block. The selected tree plantations are pure stands of exotic tree species, such as woodlot with trunks with Acacia mearnsii, cypress or eucalyptus. Most of the undergrowth in these plantations was cleared, grazed, and/or burnt. The agricultural land are smallholder farms with fields of maize and beans. The geographic location of our study area with the habitat types and census plots set across the study area is shown in Fig. 1.

Data collection
In total, we set up 117 bird census plots across the four habitat types, with the following number per habitat type: Near-natural forest block Chawia (30 plots), degraded forest patches Fururu (18 plots) and Ndiwenyi (12 plots), plantations of exotic trees (26 plots), and agricultural land (21 plots). Plots were placed across the landscape, with minimum distances of at least 200 m between plots, to minimize potential effects from autocorrelation. In addition, each plot was set at least 50 m inside of each habitat type to reduce potential edge effects. Furthermore, plots were set at sites which were easily accessible (e.g. along 1 3 existing roads, tracks and paths, especially inside of the cloud forest block Chawia and in the two forest fragments). The distribution of all plots is displayed in Fig. 1. We recorded birds at these plots during August 2018. Each bird count session lasted 10 min. Bird counts were performed during morning (between 7 and 10am) and afternoon (between 3:30 and 6 pm). We visited each plot at least three times, including morning and afternoon sessions for each plot (cf. Bibby et al. 1998Bibby et al. , 2000. Bird recordings were performed during similar weather conditions (sunny to partly cloudy weather) and by the same team consisting of three ornithologists. The observers considered all birds within a radius of about 50 m, but excluded birds which were flying over the plot. Most bird individuals were identified to species level based on calls and visual observation (cf. Bibby et al. 1998Bibby et al. , 2000. The consideration of bird calls reduces potential biases from higher visual detectability in open landscapes than in dense cloud forests (Kellner and Swihart, 2014). The complete list of all bird species observed including the number of individuals assessed is given in Appendix S1

Guilds
We grouped all birds observed into guilds, according to their diet and habitat preference (i.e. forest-dependency). We considered the following diet classes according to the major food source: raptors, insectivores, nectarivores, granivores, frugivores, and omnivores Fig. 1 The study area in Kenya (star in small inlet map), locations of census plots across the study area (left map) and the plot locations for study region Fururu and Ndiwenyi in the upper map and for Chawia in the bottom map on the right. Plots are indicated according to referencing habitat type for near-natural forest (square), degraded forest fragments (cross), agricultural land (circle), and plantations of exotic trees (triangle) (Kissling et al. 2007). Dietary items for all recorded birds were referenced from Kissling et al. (2007), by considering the major food item for each bird species. Dietary preferences may strongly shape bird assemblages (see Githiru et al. 2005). We grouped species according habitat preferences as follows: forest specialists (i.e., forest interior species), forest generalists (i.e., also found in gardens and shrubs beyond the forest), savannah species/forest visitors, and non-forest species (i.e., occurring mainly in open landscapes) (Bennun et al. 1996). Guild classification is compiled in Appendix S1.

Land cover analyses
Land cover assessment was conducted based on multispectral satellite imagery of the SPOT series (Centre National d´Etudes Spatiales 2015). We used SPOT 7 images, with a 6-m geometric resolution and 4-band spectral resolution (blue, green, red and NIR). These satellite images were taken during October 2018. The images were radiometrically corrected using the Geosud ToA Reflectance Plugin within QGIS (Ose 2015). Initial land cover classification was done using the semi-automatic classification plugin for QGIS (Congedo 2016). Further details on the procedure are given in Teucher et al. (2020). We considered the following land cover categories (according to the FAO Land Cover Classification System, Di Gregorio (2016), not identical with the naming of the habitat types used in this study so far): Indigenous forest (hereafter referring to natural forest), mixed forest (i.e. natural forest mixed with exotic trees), exotic tree plantation, bushland/shrubland/ woodland (hereafter referred as bushland), cropland, built-up area and bare soil/rock and other unconsolidated material, and grassland. The classes 'cloud' and 'cloud shadow' were assigned to the class 'unclassified' (1.89 km 2 , 1.73% only). Classification accuracy was assessed with an error matrix and Kappa statistics (Monserud and Leemans 1992; Pontius 2000). We achieved an overall accuracy of 80.7% according ground truthing data, expert knowledge and inspection of Google Earth images. Due to the fact that the neighbouring habitats and the surrounding landscape structure (e.g. heterogeneity and connectedness) strongly impact and shape species community structures in a habitat and can have even greater effects on a local population than the quality of the habitat (Metzger 1997;Mitchell et al. 2006), we considered data on adjoining land cover (within a 50-m circle) and the surrounding landscape structure (within a 1000-m circle). Thus, we calculated the proportional share of each land cover category using the zonal statistics function in QGIS for a 50-m circular buffer around each plot. For a 1000-m circular buffer we calculated the landscape configuration with the program Fragstats 4.2 (McGarigal 2015). Here we considered the following metrics: 1. Clumpiness index (landscape texture/aggregation in terms of the spatial intermixing of different land cover categories); 2. Standard deviation of the clumpiness index (a measure of variability); and, 3. Patch cohesion index (measure of the physical connectedness of the land cover categories). Further details on landscape configuration and composition can be found in Teucher et al. (2020). Details on land cover and landscape configuration for each plot are provided in Appendix S2.

Statistics
To infer habitat specific differences in species richness and abundances of bird feeding and habitat guilds we used general linear modelling (GLM) analysis with habitat types as predictors and abundance and numbers of species as response variables. Additional response variables were the first and second dominant eigenvectors of a principle coordinates analysis (Bray-Curtis dissimilarities) applied to the 65 species × 117 plots abundance matrix. These eigenvectors captured the variability in community structure among the habitat types. In a second GLM analysis, landscape features (clumpiness, SD of clumpiness, and connectivity) served as predictors. Predictors and response variables were only moderately intercorrelated (Appendix S3). Canonical correspondence analysis in combination with two-way Permanova (using Bray-Curtis dissimilarity) served to infer differences in bird community composition among habitat categories and study sites. We used t-tests for simple pairwise comparisons. We identified habitat specific bird guilds using correspondence analysis (seriation) in combination with k-means clustering. Calculations were done with Statistica 12.0.

Results
In total, we recorded 2,221 bird individuals, which belong to 65 bird species (Table 1). Plots in disturbed forest fragments and agricultural land provided highest species richness (Table 1). Agricultural land and near-natural forest (Chawia) yielded the highest number of individuals (Table 1). However, the average numbers of species per plot did not significantly differ between near-natural forest (Chawia), the forest fragments (Fururu and Ndiwenyi) and agricultural land (Table 1). In turn, bird abundances of plots in the agricultural land were significantly higher (t-test P < 0.01) than in the near-natural forest and the forest fragments (Table 1). In the plantations, only two individuals of the Chin-spot Batis (Batis molitor), a cooperative breeding insectivore bird species, were recorded. All other birds found there were flying over and thus were not considered in our analyses.
When classifying birds according to their feeding guilds, frugivores and insectivores achieved the highest average abundances in near-natural forest (Fig. 2a). Granivores and nectarivores dominated in the agricultural landscape (Fig. 2a), and insectivores and frugivores reached highest abundances in the forest fragments (Fig. 2a). The agricultural landscape was highly dominated by forest visitors, while forest specialists and generalists dominated in the near-natural forest (Fig. 2b). The forest fragments were evenly inhabited by all four habitat guilds (ANOVA: P > 0.05; Fig. 2b).
Two-way Permanova revealed significant differences in bird community composition between the three habitat types (excluding plantations) and the three study sites (Table 2). Of the 65 species, 28 (43.1%) were habitat specific (Fig. 3). Nevertheless, as much as 18 species occurred in all three habitat types (near-natural forest, degraded forest fragments, and agricultural land, Fig. 3). A high compositional similarity also occurred between the Table 1 Basic richness and abundance data. Given are number of sites, species richness (total and mean (± SD)), and abundance (total and mean (± SD)) according to the three habitat types studied 1 3 agricultural sites and the forest fragments (Fig. 3). A substantial number of forest generalists and forest visitors were present in all three habitat types (Fig. 3). Forest specialists were not seen in open agricultural landscapes (Fig. 3). Interestingly, 5 out of the 19 open landscape species were only recorded in the forest (Fig. 3). Correspondence analysis identified four major groups of habitat preferences (Appendix S3, Figs. A1, A2). Two of these groups with 17 and 15 species, respectively, consist of pure and partial open landscape dwellers, the latter with occasional occurrences in the forest fragments and the near-natural forest (Appendix S3). A third group of 26 species are true habitat generalists (Appendix S3, Fig. A1). Only nine species are true forest specialists (Appendix S3, Fig. A1). Canonical correspondence analysis ordered the bird feeding and habitat guilds along the dominant axis that correlated with the gradient from near-natural forest to agricultural open landscape (Appendix S3, Fig. A2). Finally, we analyzed potential effects of land cover type and land configuration to bird species richness, abundance and community composition (Table 3). We found a moderate negative correlation of the percentage of mixed forest with overall bird abundances (Table 3). The percentage of forest was positively linked to bird species richness (Table 3). Further, the pairs of forest and percentage of building, and of forest and agricultural land had opposite effects on community structure (Table 3). The variability in habitat    clumpiness was negatively correlated with species richness and also affected the composition of the bird species community (Table 3).

Discussion
We found differences in bird species richness and community structure among the four habitat types. Disturbed forest fragments and the agricultural land provided higher species richness than the near-natural forest block and plantations of exotic trees. The composition of the bird species communities was relatively similar in the degraded forest fragments and the surrounding agricultural land, while the bird community composition found in the near-natural forest block diverged considerably, where forest-dependent bird species occur in comparatively high densities. In the following we will discuss 1. The value of large nearnatural forest blocks; 2. The relevance of forest fragments and plantations as potential habitat surrogates for forest species; and 3. The role of the environment in which habitats are embedded.

The value of large and near-natural forest
Nine out of the 65 bird species observed were exclusively found in the near-natural forest block, Chawia. A closer examination of species composition across the four studied habitat types indicated that forest specialists are mainly found in the near-natural forest block, and did not occur in open agricultural land and only rarely colonize degraded forest fragments (although possibly as temporary visitors rather than as permanent residents). For example, the Sharpe's Starling Cinnyricinclus sharpii is a scarce forest specialist of the Afrotorpical highland biome (Bennun and Njoroge 1999) that relies on fruiting forest trees (Lewis and Pomerory 1989). This regionally near-threatened species ) was exclusively observed in the interior of the large near-natural forest block, Chawia, but was not found in the small and degraded forest fragments. Our data revealed further differences in bird species community structures. Frugivores and insectivores were found in the near-natural forest block Chawia, while granivorous bird species dominated the bird species community in the agricultural land. The degraded forest fragments supported frugivores and insectivores, with less of forest-dependent guilds. Previous studies on bird species community structures in Taita Hills forest fragments underline that the abundance of bird species with a broad variety of food demands (insectivores, frugivores, nectarivores) are frequently rare in fragments with low habitat quality, while omnivore bird species remained stable across all forest fragments, irrespective of the degree of habitat quality (Ulrich et al. 2016). In conclusion, large near-natural forest blocks are of relevance to maintain forest-dependent bird species, requiring specific forest resources and habitat pre-conditions.

Diversity in small and disturbed forest patches
Our study showed higher bird species richness for small and disturbed forest fragments and agricultural land if compared with the near-natural forest block. This goes in line with other studies indicating higher bird species diversity in disturbed habitats than in rather intact ones (Fischer and Lindenmayer 2002;Woltmann 2003;Bett et al. 2016). These differences in species richness between disturbed and intact habitats are 1 3 primarily attributed to the additional occurrence of widespread species in disturbed habitats (Kitahara and Fujii 1994;Börschig et al. 2013; but also see Lenda et al. 2012;Tropek et al. 2013 for opposite relationships). Similarly, transition zones between two habitat types, such as forest edges, are often richer in species than the interior of an intact forest (Habel et al. 2018). This is due to the fact that transition ecotones as well as disturbed (forest) habitat fragments mostly provide a larger variety of ecological niches and resources, and might become colonized from both forest and open-land dwelling species.
Our study showed that forest-dependent birds also occurred in small and degraded forest fragments, but in lower abundances. This could mean either that at least some habitat specialists are still able to survive in disturbed habitats, or that forest specialist species are visitors of small forest patches. However, small and isolated forest fragments are generally considered as having reduced habitat quality, and might not sustain rich communities for longer time (Ulrich et al. 2016). Despite the fact that the forest fragments do not harbour some of the forest-dependent bird species and show divergent forest community structure if compared with the bird species community of the near natural Chawia forest block (see also Simberloff and Abele 1976;Muwal et al. 2012), these smaller fragments may serve as important stepping stones. This allows species to migrate among three larger forest blocks across Taita Hills (see Pellikka et al. 2009) and thus enhance the opportunity for biotic exchange and the persistence of species at the landscape level (Blake and Karr 1984).

Anthropogenic surrogate habitats for forest species
We found very low bird species richness and abundance in the plantations of exotic trees. All plantations assessed in this study were stands of adult Acacia mearnsii, eucalyptus and cypress species. The undergrowth was mainly used for cattle grazing, and the harvested wood as an energy source (i.e., firewood) and for housing construction. Thus, the plantations were mostly very light and without a shrub layer. The lack of plant undercover and flowering plants might lead to a rather low number of insects. As a consequence, such habitat structures exclude frugivorous, granivorous and nectarivorous bird species. Previous studies showed that bird species richness and abundances in plantations (versus natural forest) is strongly driven by local and temporal food availability (Mitra and Sheldon 1993). In contrast to our results, other studies have shown that tree plantations, even when consisting of a single exotic tree species, can provide a suitable habitat for forest species (Dent and Wright2009;Castano-Villa et al. 2019) and thus act as important surrogate habitat, which become relevant in regions suffering under high rates of deforestation and forest degradation (e.g. the Taita Hills). Thus, there is a controversial debate on the ecological value of tree plantations for (forest) species (Marsden et al. 2001;Berndt et al. 2008;Farwig et al. 2008;Calviño-Cancela 2013). The ecological value of plantations can be significantly improved by an increase of habitat heterogeneity through planting of native plant species undergrowth (Hartley 2002;Nájera and Simonetti 2010). However, this is completely lacking in our study sites. Given that the plantations and small forest patches only hold a small proportion of species richness and abundance of forest-dependent species, it remains questionable to what extent the unique biodiversity of Taita Hills can be preserved in the long term, as only three larger forest blocks are left (Lens et al. 2002).

Impacts from the surrounding environment
We found that the effects of the land cover around a census plot on bird species richness and diversity were of only little relevance. Less habitat clumpiness was positively correlated with species richness, and also affected the composition of the bird species community significantly. Clumpiness quantifies landscape texture in terms of the spatial intermixing of different land cover categories. More habitat clumpiness means more like adjacencies in the landscape, higher connectivity within habitats and less habitat heterogeneity. In turn, less clumpiness (i.e. increasing habitat composition and/ or configuration in the landscape, and more disturbance) is associated with higher species richness and population persistence. This goes in line with other studies analysing landscape configuration (e.g. clumpiness) and biodiversity (see Olsoy et al. 2016). The effect caused by the surrounding environment could be due to the fine-grained mosaic landscape structure in combination with the comparatively high mobility of birds.

Caveats
As any comparative study, ours might be affected by various shortcomings that could limit the significance of our results and conclusions: 1. Study duration: Our surveys were conducted within a short period of time, and seasonal fluctuations might significantly change the spatial distribution of bird species (see Mulwa et al. 2013). 2. Limited sample size and restricted sample coverage: Our study represents the situation only for a very restricted area within the Taita Hills. Therefore, the generality of these results and correlations obtained should be considered with caution. 3. Autocorrelation: At the regional scale, habitat types are not randomly distributed, but occur clustered depending on land use and soil features. Consequently, any comparison of habitat types might be affected from spatial non-independence of data points (Hawkins 2012). 4. Diverging detectability in different habitat types: Lower species richness in near-natural forest could have been the result of under-detection of birds due to high vegetation density or tree heights in near-natural forest compared to disturbed forest fragments and agricultural land, and also because of forest specialist species tend to occur in lower densities than generalist species. We partially countered this bias by including bird calls in addition to visual observations (cf. Bibby et al. 1998Bibby et al. , 2000.