Importance of microhabitat selection by birds for the early recruitment of endangered trees in a fragmented forest

Remnant microhabitats are important for bird habitat selection and plant regeneration in the fragmented habitat. However, empirical information on the consequences of how microhabitat use by birds affects the early recruitment of plants is lacking. In this study, we evaluated whether microhabitat selection by the Black Bulbul (Hypsipetes leucocephalus) (J. F. Gmelin, 1789) impacts the early recruitment of the endangered tree species, the Chinese Yew (Taxus chinensis (Pilger) Rehd), in a fragmented forest over a 4-year period (2011–2012, 2018–2019). Our results showed the main factors affecting H. leucocephalus microhabitat selection were distance to the nearest T. chinensis mature tree, herb cover and density, leaf litter cover, and vegetation type. Moreover, the results of logistic regression also highlighted the importance of elevation, distance to light gap and roads, tree cover in bird microhabitat selection. Furthermore, the seed emergence rate in microhabitats used by birds did not differ from the natural forest, which was related to five factors of bird microhabitat. The Random Forest model showed that seedling emergence rate was increased with leaf litter cover and distance to fallen dead trees, but decreased in relation to herb cover, slope, and elevation. Our results highlight the importance of remnant microhabitats in fragmented forests for sustaining forest ecology and optimal management. The contribution of microhabitats used by birds to plant recruitment provides insights into how frugivore species contribute to plant regeneration, which should be incorporated in future conservation and management practices of fragmented forests.


Backgrounds
Understanding how habitat fragmentation affects the persistence of plant populations is a central component of forest ecology and management (Taubert et al. 2018;Liu et al. 2019;Peters et al. 2019). In the case of fleshy-fruited trees, plant regeneration in remnant forest patches depends on the coupling of seed dispersal and seedling recruitment processes (Cordeiro and Howe 2003;Bregman et al. 2016;Farwig et al. 2017). Empirical evidence suggests that habitat fragmentation could disrupt plant regeneration by altering seed dispersal processes and reducing the availability of suitable microhabitat for regeneration (Bomfim et al. 2018;Emer et al. 2018;Simmons et al. 2018;Marjakangas et al. 2019). Nevertheless, the fundamental question remains of whether recruitment failure in fragmented forests is caused by a greater limitation of seeds available for dispersal or by post-dispersal processes (i.e., seedling establishment and
For bird-dispersed plants, birds transporting seeds away from the mother plants affect seed deposition in fragmented forests (Schupp et al. 2017;Li et al. 2020) and, hence, patterns of plant regeneration (Donoso et al. 2016;García-Cervigón et al. 2018). Often, birds exhibit complex behavioral pattern in response to forest fragmentation, depending on the distribution of food and other resources (e.g. shelter, nesting, vigilance or resting sites ;Cody 1985;Côrtes and Uriarte 2013). For example, the loss of food resources by fragmentation might cause the number of frugivorous bird species to decline, thus reducing the amount of seeds removed and disrupting seed dispersal process (Pérez-Méndez et al. 2015;Farwig et al. 2017;Zwolak 2018). After foraging, birds often exhibit highly specific microhabitat selection. Consequently, a disproportionate number of seeds are deposited at sites selected by bird dispersers, negatively impacting future plant regeneration (Sasal and Morales 2013;Li et al. 2019).
Clearly, seed deposition patterns by birds depend on their behavioral decisions (Schupp 1993;Jordano and Schupp 2000;Cousens et al. 2010). Key behavioral processes include seed removal and post-foraging microhabitat use (Schupp et al. 2017). However, while the regeneration of plant populations in fragmented forest depends on where seeds are deposited by birds (Spiegel and Nathan 2007;Lehouck et al. 2009;Carlo et al. 2013), sites must also be suitable for the early regeneration of plants (Puerta-Piñero et al. 2012;Schupp et al. 2017). Although many studies have highlighted the role of birds in seed removal and seed deposition in fragmented forests (Farwig et al. 2017;Bomfim et al. 2018;Marjakangas et al. 2019), empirical evidence of the consequence of bird microhabitat use on seed germination in fragmented forests is lacking.
In Southeast China, the endangered Chinese Yew, Taxus chinensis (Pilger) Rehd, is a dominant tree species in fragmented forests, and it mainly depends on bird dispersal for regeneration (Li et al. 2019). Over 80% seeds of T. chinensis were removed by the Black Bulbul, Hypsipetes leucocephalus (J. F. Gmelin, 1789) (resident species in Fujian, weight: 41-62 g), indicating that it is the most important species in T. chinensis-bird mutualism (Li et al. 2015). Here, we evaluated whether post-foraging microhabitat selection by H. leucocephalus impacts the early recruitment T. chinensis in a fragmented forest over a 4-year period (2011)(2012)(2018)(2019). Specifically, we examined: (1) which factors determine bird microhabitat use; and (2) how the microhabitat selection of H. leucocephalus impacts the early recruitment of T. chinensis in a fragmented forest. The results of this study are expected to demonstrate the importance of frugivores in facilitating the regeneration of tree species, with implications on conservation and management practices in remnant fragmented forests.

Species and study site
T. chinensis is a dioecious and wind-pollinated species that is distributed in evergreen broadleaf forests. Every year, female plants bear axillary cones which, in autumn, develop into fleshy arils (commonly, although incorrectly, referred to as "fruits") that contain a single seed. An average tree bears more than 4000 of these "fruits" (Li et al. 2015(Li et al. , 2019. This study was conducted in a yew ecological garden (elevation 895-1218 m above sea level [a.s.l.], slope gradient 27°), located in the southern experimental area of the Meihua Mountain National Nature Reserve (25°15′-25°35′N, 116°45′-116°57′E) in the west part of Fujian Province, China. This site contains the largest natural population of T. chinensis in China (approximately 490 adults, distributed in the evergreen broadleaf forest), including 200 trees that are > 500 years old. A national forest garden of 15 ha was established by the government in 2003 to protect these endangered trees. Due to long-term of human use, the vegetation around the forest garden is highly fragmented. The remnant evergreen broadleaf forest patch is interlaced with bamboo patches and mixed bamboo and broadleaf patches to form a fragmented forest. The dominant tree species of the remnant evergreen broad-leaved forest is T. chinensis (Additional file 1: Fig. S1).

Microhabitat selection by Hypsipetes leucocephalus
To study post-foraging microhabitat selection of H. leucocephalus, field observations were made after the birds departed mature T. chinensis plants. We observed the post-foraging perching position of H. leucocephalus using a field scope (Leica 70, Germany) at distances of 50-100 m from the opposite mountain slopes. When the position of birds was recorded, we collected regurgitated seeds in the canopy of bird preferred microhabitats (regurgitated seeds refer to cleaned seeds without any aril/coat that had been totally digested by birds). Because we tried to show the relationship between bird habitat use and plant recruitment, we chose the site with regurgitated seeds as bird-preferred microhabitat and set 1 m × 1 m quadrats. To test whether bird selected habitat, we also set the quadrats in other available areas, and the position were confirmed by random number table. Totally, 30 used and 30 available quadrats were set to collect information on microhabitat factors in 2011. To exclude year-to-year variation in bird habitat selection, we recorded perching frequency at these 60 sites in the other study years (2012,(2018)(2019).
In both bird-use and available quadrats, we measured three qualitative factors: aspect (shade slope; sunny slope), vegetation type (bamboo forest; Chinese Yew forest; farmland; mixed bamboo and broadleaf forest) and heterogeneous tree species (other tree species, except T. chinensis trees). Moreover, we also measured 15 quantitative factors: elevation, slope, distance to water, distance to light gap, distance to roads, distance to nearest T. chinensis tree, distance to nearest heterogeneous tree, distance to nearest T. chinensis mature tree, distance to fallen dead tree, herb cover, herb density, shrub cover, shrub density, tree cover, and leaf litter cover.
For analyzing microhabitat selection by birds, we first compared three qualitative factors by Chi square test. The other 15 quantitative factors evaluated between bird used, and available quadrats were first analyzed by a t-test. All quantitative variables were evaluated with Principal Component Analysis (PCA) based on their correlation matrix with a varimax rotation to screen out the key factors in microhabitat selection of H. leucocephalus. PCA is a multivariate technique that produces a simplified, reduced expression of the original data with complex relationships, and has been widely applied in studies of wildlife habitats (Fowler et al. 1998). We also used logistic regression to explore the role of microhabitat factors for bird habitat selection.

Effects of microhabitat selection by birds on the early recruitment of Taxus chinensis
Independent of the seed dispersal study, seedling emergence was assessed experimentally in 2017 and 2018 beneath the 30 sites used by H. leucocephalus. At each site, 200 seeds were sown at a depth of 1 cm to avoid predation by rodents. The germinated seedlings were monitored weekly from spring to fall of the following year. Plant early recruitment was computed as the fraction of germinated seedlings that survived to the end of the first fall.
To study the effects of bird microhabitat selection on the early recruitment of plants, we used the t-test to compare the seedling emergence rate in the fragmented forest with natural conditions (seedling emergence rate: 10.86%, with 1000 seeds sowing under the canopy of 10 microhabitats in the natural habitat; Gao 2006). Random Forest model is an ensemble machine-learning method for classification and regression that operates by constructing a multitude of decision trees. It is appropriate for illustrating the nonlinear effect of variables, can handle complex interactions among variables and is not affected by multicollinearity. Random Forest can assess the effects of all explanatory variables simultaneously and automatically ranks the importance of variables in descending order. Then, we used the Random Forest (RF) algorithm to evaluate the quadrat habitat factors selected by birds in relation to the number of germinated seedlings (R package Random Forest) (Breiman 2001).

Results
After foraging, H. leucocephalus exhibited strong microhabitat selection. The Chi square test showed that vegetation type (Chi square test: χ 2 = 6.300, p = 0.043) and slope aspect (Chi square test: χ 2 = 9.600, p = 0.002) varied between microhabitats used by birds and those that were available. H. leucocephalus preferred shade slope and bamboo patches. Highly significant differences were detected between microhabitats used by birds and those that were available when considering distance to the nearest T. chinensis tree, distance to the nearest heterogeneous tree, distance to the nearest T. chinensis mature tree, shrub density, and leaf litter cover (Table 1).
Importantly, the PCA results highlighted that the main factors affecting H. leucocephalus microhabitat selection were distance to the nearest T. chinensis mature tree, herb cover, herb density, leaf litter cover, and vegetation type. Distance to light gaps and nearest heterogeneous trees were also important for bird microhabitat selection (Table 2). Moreover, the results by logistic regression showed habitat selection of birds was only affected by elevation, distance to light gap and roads, tree cover (Table 3; the other eleven factors did not affect habitat selection, All p >0.05). As a consequence of bird habitat selection, the sowing experiment first showed a 5%-15.5% seedling emergence rate in the H. leucocephalus microhabitat, which was not significantly different from natural conditions (t = 1.679, p = 0.104). Considering the seed germination related to habitat selection, the Random Forest model showed that seedling emergence rate increased with leaf litter cover and distance to fallen dead trees, but decreased in relation to herb cover, slope, and elevation (Random Forest: 62.43% of germinated seedlings could be explained by five variables) (Fig. 1).

Discussion
After foraging, H. leucocephalus exhibited strong microhabitat selection. Consequently, H. leucocephalus microhabitats influence the early recruitment of T. chinensis.
Because microhabitat characteristics vary in fragmented forests, they influence the microhabitat selection of birds. An optimum suitable microhabitat provides safe shelter for bird to avoid predation and an opportunity to access reliable food resources (Cody 1985). With forest fragmentation, the remnant microhabitat was important for the ecology and management of the fragmented forest. In our fragmented forest, bamboo patches and shrub density potentially supply safe shelter for H. leucocephalus. Microhabitats close to T. chinensis trees, heterogeneous trees, and T. chinensis mature trees could meet the two requirements of safety and food for H. leucocephalus. The requirement for optimum temperature is also an important component in the strategy of bird habitat selection (Moore 1945). Suitable microhabitat temperature is maintained by shrubs and herbs, providing shelter from wind and rain (Kelty and Lustick 1977;Cody 1985). In the current study, shrub density, distance to light gaps, herb cover, herb density, and leaf litter cover were preferred by H. leucocephalus, possibly because they supply suitable microhabitat temperature.
The remnant microhabitats selected by birds were important for plant recruitment (García-Cervigón et al. 2018). The seed emergence rate of T. chinensis beneath the microhabitats used by birds showed no significant difference to natural conditions; thus, H. leucocephalus is likely important for the early recruitment of T. chinensis. Furthermore, the emergence of T. chinensis seeds was influenced by bird microhabitat factors. Sites with low herb cover facilitated seed emergence, due to the low competitive ability of trees. T. chinensis (Li et al. 2015). Living with dense leaf litter cover could also provide an important supply of nutrients for seed germination. Because of the requirements for T. chinensis seedlings to germinate (Li et al. 2015), shaded slopes might supply a shaded microenvironment for seed germination.

Conclusions
This study demonstrated the microhabitats used by H. leucocephalus affected the early recruitment of T. chinensis. Our results also highlight the importance of remnant microhabitats in fragmented forest for  sustaining forest ecology and enhancing management practices. The contribution of habitats used by birds to sites of plant recruitment could be used to determine how frugivore species facilitate plant regeneration. Such information could be incorporated in future conservation and management practices to facilitate the regeneration of fragmented forests. However, our study partially reflected the consequence of bird habitat selection on plant recruitment, owing to sowing the seeds in bird preferred habitat. Researchers need to consider in future studies including the sowing experiment in both bird preferred and available habitat with bird regurgitated seeds and natural fallen seeds. Furthermore, future studies also need to compare how bird habitat selection affect plant recruitment in both fragmented and continuous habitat, which could explore the effects of habitat fragmentation on bird-plant mutualism. Fig. 1 Microhabitat selection by Hypsipetes leucocephalus and how these parameters (distance to fallen dead tree, distance to leaf litter cover, herb cover, slope, and elevation) affect the seedling emergence rate of Taxus chinensis in a yew ecological garden representing a fragmented forest in Southeast China. Results were determined using the Random Forest algorithm, and show the partial effects of the five independent variables on seedling emergence rate