Improvement of vegetation structure enhances bird functional traits and habitat resilience in an area of ongoing restoration in the Atlantic Forest

: Ecological restoration is a traditional option for recovering biodiversity and ecosystem functions. Birds perform pollination, seed dispersal, and pest-control services, which catalyze increases in habitat structure. Habitat complexity changes bird composition, but there is little evidence of its effects on bird functional diversity in Neotropical restorations. We tested whether bird functional diversity and composition respond to increased habitat complexity. Point-counts were performed (January-December 2015) in an area undergoing restoration (536 ha) in the Atlantic Forest of southeastern Brazil, in restorations with less and more structured vegetation and pastures and forest-fragments. The functional bird traits considered were diet, habitat, biomass, environmental sensitivity, and foraging strata. Increased habitat complexity was evaluated using plant characteristics (exotic grass, canopy, herbaceous cover, and diameter at breast height). A total of 172 bird species (5% endemic; 12% migratory) were recorded. Increased vegetation structure in both restored sites and forest-fragments drove a reorganization and addition of functional bird traits, which positively influenced functional richness, dispersion, and evenness. Shifts in plant-characteristics rearranged bird functional traits (diet-forest-dependence and diet-strata-foraging). The rapid development of vegetation structure is a key factor for restoration because it provides additional habitat for semi-dependent forest birds and enhances resilience and sustainability in new man-made


INTRODUCTION
The expansion of human activities and the development of agriculture, livestock, and urbanization are responsible for the collapse of pristine tropical forests worldwide (Houghton 1994, Myers et al. 2000, Laurance et al. 2014. Therefore, the restoration of more than two billion hectares of former native vegetation is now a global priority (Minnermayer et al. 2011, Crouzeilles et al. 2016. Ecological restoration aims to recover degraded lands (SER 2004, Brancalion et al. 2013 and rescue a set of their biodiversity, ecological interactions, and ecosystem services (Chazdon et al. 2009, Hobbs et al. 2009, Rey-Benayas et al. 2009), and has been a widely-adopted strategy by ecologists and decision-makers (Brancalion et al. 2013).
Restoration success can be evaluated by the ability to improve biodiversity and ecosystem functions (Sullivan et al. 2018, Batisteli et al. 2018. The restoration of animal communities has long remained in the background, following the idea that fauna would be passively rescued by only recreating suitable conditions (Palmer et al. 1997). Nonetheless, the restoration of faunal assemblages is highly complex in the tropics due to the awesome amount of biodiversity and ecological interactions involved (da Silva et al. 2015). Birds possess characteristics that make them good models for evaluating whether restorations are reaching their goals. These vertebrates are highly diverse (Del Hoyo et al. 2019), extremely vagile (usually the first recolonizers of restored sites) (Munro et al. 2011), and perform vital ecosystem functions (e.g., pollination, seed dispersal, and pest control) (Medellin & Gaona 1999, Morrison & Lindell 2012. The ecosystem functions provided by birds are essential for habitat maintenance and are strictly connected to their multiple ecological traits . Functional traits are individual morphological, physiological and/or phenological attributes (Laliberté & Legendre 2010) that have been used to calculate functional diversity (FD) (Laliberté et al. 2014) and to understand how birds interact in an ecosystem (Violle et al. 2007). Positive responses of bird FD to vegetation diversity can be modulated by multiple associations between animal traits and vegetation structure (Sitters et al. 2016). Likewise, FD metrics are sensitive to vegetation structure of restored sites (Batisteli et al. 2018) and have been considered better metrics for evaluating restoration success than taxonomic composition (Brancalion & Holl 2016). The relationship between FD and habitat resilience is linked to the maintenance of a high number of species with similar functions yet distinct responses to disturbance (Arruda Almeida et al. 2018).
The Atlantic Forest is considered a global biodiversity hotspot (Myers et al. 2000). Currently highly fragmented, the Atlantic Forest persists as small-isolated remnants surrounded by a matrix of pastures and croplands (Ribeiro et al. 2009, Calaboni et al. 2018. This scenario makes the coexistence of bird conservation and agriculture production in this biome a great challenge (Uezu et al. 2005, Uezu & Metzger 2016, Piratelli et al. 2019). On the other hand, the Atlantic Forest offers great opportunities for habitat restoration (Minnermayer et al. 2011), and active restoration has become a central strategy for reducing the effects of forest fragmentation and to avoid local extinctions of forest birds (Uezu & Metzger 2016).
The effects that vegetation structure has on bird FD has been studied in ecoregions across the world. In general, bird FD metrics increase with increasing vegetation structure in temperate forests in North America, Europe and Asia (Bae et al. 2018). However, distinct responses have been found in tropical zones. In Australia, FD indices were found to be positively related to vertical vegetation diversity in humid forests, yet inversely related to vertical vegetation in dry forests (Sitters et al. 2016). A similar pattern was also found in savannahs in Namibia (Seymour et al. 2015). These results illustrate that tropical ecosystems may have different responses than temperate habitats, and extrapolations based on FD indices found in distinct ecoregions can lead to unsuccessful conservation strategies for the megadiversity of the Neotropics (Freitas & Mantovani 2018). Nonetheless, there is a large knowledge gap regarding how bird assemblages and FD respond to the increased vegetation structure from active restoration process in Neotropical ecosystems, which limits the assessment of the efficacy of restoration actions (Ortega-Alvárez & Lindig-Cisneros 2012). Several local-scale active restorations have been undertaken in the Brazilian Atlantic Forest (Rodrigues et al. 2011), but there is little evidence regarding the effects of the development of vegetation structure on bird taxonomic composition and the conservation of their functional traits in these systems.
Here we address how a recent (i.e., less than 10 years old; Twedt et al. 2002) active restoration program with high-diversity plantation may drive the taxonomic and functional diversity of birds, compared to pasture and small native forest fragments. Specifically, we tested whether increased vegetation structure of restored areas changes the arrangement of bird assemblages. Thus, we evaluated whether bird species richness (SR) and FD indexes (based on bird traits regarding diet, biomass, forest dependence, foraging strata, and environmental sensitivity) differ between restored and non-restored habitats due to non-random environmental changes linked to increased vegetation structure. These non-restored habitats were represented by pasture (as a degraded reference) and small native forest fragments (as meta reference). Although from a restoration ecology perspective the restored habitats are young (< 10 years old post-planting) (Twedt et al. 2002), we still expected to find differences in vegetation structure because they were first visually checked and then confirmed as different by comparing the proportion of grass and canopy cover, diameter at breast height and tree morphorichness). Increased vegetation structure in new man-made forests might be a central tool for recovering native bird species and their essential functional traits to increase ecosystem resilience and sustainability. We predicted the following: Less structured restorations, dominated by exotic grass and with little canopy cover and low tree morphorichness (Melo et al. 2007), should maintain bird SR and FD similar to that of pasture, but different from more-structured restorations and small native forest fragments. Degraded habitats dominated by exotic grass cover sustain simplified bird assemblages composed of omnivorous, granivorous, forestindependent and low sensitivity species (Becker et al. 2013, Casas et al. 2016); More structured restorations with higher tree morphorichness, more canopy cover, and less exotic grass cover, should provide more microhabitats, niches, and resources for birds. Thus, they are likely to maintain levels of bird SR and FD above that of pasture and less-structured restored habitats -both dominated by exotic grass -yet below that of small native forest fragments. Active restorations add vegetation structure (e.g., herbaceous and tree traits), which provides habitats for several functional groups of birds (e.g., frugivores, nectarivores and insectivores) (Becker et al. 2013, Batisteli et al. 2018).

Study area
This study was carried out at "Centro de Experimentos Florestais", an area of 526 ha in eastern São Paulo State, southeastern Brazil (23°14'15.18" S, 47°24'3.29" W; 580 m a.s.l). The climate is characterized by dry winters and hot summers, with average monthly precipitation of 56 mm and 160 mm, respectively (Alvares et al. 2013). The vegetation is semideciduous seasonal forest, one of the phytophysiognomies of the Atlantic Forest (Veloso et al. 1991).
The native forest fragments in the area are of different ages and were historically deforested and converted to coffee (Coffea arabica L.) plantations in the early 20 th century (César et al. 2013). Prior to restoration actions, the coffee plantations were abandoned and converted to pastures dominated by invasive exotic grass (Urochloa spp.) for raising cattle by intensive system (Amazonas et al. 2018). These land uses drastically affected the native vegetation arising in the highly fragmented landscape, with few small and isolated secondary forest remnants (Amazonas et al. 2018, Andrade et al. 2018. In this context, a forest restoration program was initiated in 2007 with the aim of restoring a set of the native vegetation, fauna and ecosystem functions and services (Gagetti et al. 2016). A high diversity planting, using ~720,000 seedlings from more than 100 native tree species, was performed randomly with 3 x 2 m spacing, composed of pioneer and secondary species (see Gagetti et al. 2016). The restored area encompasses nearly 400 ha, with planting age ranging from four to 11 years. This area is currently a restored island in a highly fragmented landscape of Atlantic Forest, surrounded by pastures with exotic grass (Urochloa spp.) for cattle ranching, and croplands. Bodies of water (e.g., artificial lakes and swamps) are also found throughout the area.

Vegetation structure
Tree traits, such as canopy cover, diameter at breast height, tree height and increase in biomass following the age of plantings, have been used as indicators of the structural development of vegetation during forest restoration (Melo et al. 2007, Crouzeilles et al. 2016. Herein, we assume that habitat types with greater exotic grass and lower values for tree traits (e.g., canopy cover, diameter at breast height, and tree morphorichness) have less vegetation structural development (hereafter, VSD), while habitat types with higher values for these tree traits and less exotic grass cover have higher VSD. We measured plant characteristics (below) for each of 39 fixed-point bird-count sites to acquire VSD. Fixed-point counting is a commonly used method for evaluating bird assemblages by which researchers consider only birds recorded during a limited sampling time and/or in a defined area (Bibby et al. 2000; more details below). The measured plant characteristics were: (i) herbaceous plants (Herb) -visually estimated percentage of ground covered by herbaceous plants in 20-m radius plots; (ii) exotic grasses (Grass) -visually estimated percentage of ground covered by exotic grasses in 20-m radius plots; (iii) canopy cover (Canopy) -mean of four estimates, one in each cardinal direction (N, S, E, and W) at 10-m from the center of each point, made using percentage of spherodensiometer squares occupied by light passing through the foliage; (iv) diameter at breast height (cm) (DBH) -estimated from tree basal area; (v) tree morphorichness (Ric.tree) -morphologically distinct arboreal tree specimens were classified as tree morphospecies and counted at each point locations; and (vi) tree height (cm) (Tree. height) -measured from the ground to the top of tree foliage (Table I).
Even though the two restored habitats were of similar age, we visually verified in loco that they had distinct characteristics of VSD. This was confirmed by subjecting plant traits (grass cover, canopy cover, herbaceous cover, DBH and tree morphorichness) to a one-way ANOVA followed by Tukey's pos-hoc test p <0.05). The two restoration habitats were subsequently considered as distinct areas, and form here on referred to as the (1) less structured restored habitat and the (2) more structured restored habitat. Forest dependence (Parker III et al. 1996) asindependent (I) for species recorded mainly in grasslands, pastures, and marshes; semidependent (S) for species observed mainly at forest edges and also using open habitats; and dependent (D) for species recorded mainly in forest habitats; (e) Environmental sensitivity (Parker III et al. 1996) S: low (L), medium (M) and high (H) sensitivity to disturbed and degraded lands. Generalists are defined here as those species having high ecological plasticity (diet and/or environmental) and whose populations tend to grow in simplified habitats such as pastures and degraded/open-habitats (e.g., granivores, omnivores, and insectivores, here also considered as forest independent and having low environmental sensitivity). Similarly, specialists are defined as those species with some restrictions in their requirements, either for forest habitats or feeding strategies (forest dependent, and/or high sensitivity or insectivorous and frugivorous forest birds).

Bird samplings, functional traits, and other arrangements
The conservation status of bird species was assessed using official regional (São Paulo 2018) and global (IUCN 2019) red lists. Endemic birds of the Atlantic Forest followed Vale et al. (2018). Migratory birds followed Somenzari et al. (2018) as: (a) migratory (MGT) -species with populations that regularly and seasonally move away from their breeding sites and return every breeding season; and (b) partially migratory (MPR) -species with populations that are partially migratory; and (c) resident (RE), species that occupy the same area throughout the year, including nomadic birds, with minor variation in population structure. Nomenclature followed the Brazilian Ornithological Records Committee (Piacentini et al. 2015).

Statistical analysis
Bird assemblages of the four habitats were first evaluated using the iNEXT function of the homonym package (Hsieh et al. 2016). Individualbased rarefaction curves were plotted against a given number of individuals chosen randomly from observed samples until all individuals had been accumulated, whereas extrapolation curves were plotted to double the degraded reference habitat (Colwell 2012), which in the present case was pasture. Pasture was chosen as a reference because it is a regional habitat that existed before the restoration program. A total of 999 bootstrap replicates were used to estimate 95% confidence intervals, with non-overlapping confidence intervals among the habitats being assumed to reflect significant differences.
Functional distances between pairs of species were then calculated using Gower's distance metric (Gower 1966), which is suitable for calculating both categorical and continuous traits with missing trait values (Legendre & Legendre 1998, Podani 1999. Species richness (SR) and functional diversity (FD) parameters were computed, the latter as functional richness (FRic), functional divergence (FDiv), functional evenness (FEve), functional dispersion (FDis), and Rao's quadratic entropy (Rao's Q) for each point-count site using the dbFD function of the The difference between these two parameters is that FDis represents the mean distance of individual species to the centroid of all species, while Rao's Q is the mean distance between each pair of these species. Higher values of FDis indicate a greater potential for functional complementarity among species. Different from FRic, which is monotonically related to SR, FDis and Rao's Q are influenced only by species abundances (Laliberté & Legendre 2010). Moreover, the vif function in the usdm package (Naimi et al. 2014) was used to calculate the variance inflation factor (VIF) for the vegetation variables and excluded Tree.height because of a problem with multi-collinearity (VIF > 0.4). Moran's I test was then used to test for spatial autocorrelation, which was not found for SR (Moran's I: -0.02, p-value: 0.85) or FD parameters (all p-value: >0.1).
Whether vegetation structural development (VSD) shaped bird predictors as SR and FD measures in restored and non-restored (pasture and forest fragments) habitats was then tested. The dimensionality of the vegetation characteristics was reduced using principal component analysis (PCA), with each of the first two axes (1 st and 2 nd PCs) being used to frame a vegetation structure gradient (VSD 1 and 2) (Batisteli et al. 2018). Multiple regressions models were then applied to test whether VSD values (predictor variables) significantly influenced bird predictors (p < 0.05).
Finally, since the habitats were spatially close (10 to 1200 m), we assumed that birds have the same chance of occurring in all areas, except for their specific habitat requirements, here translated as VSD. Thus, the influence that VSD had on the composition of bird guilds in the four studied habitats was tested. To do this, bird traits were combined to transform them into categorical functional bird groups: (a) diet-forest dependence, which was composed of species with a combination of diet and forest dependence, such as fruit-nectar forest independent (Fn.I), semi-dependent (Fn.S) and dependent (Fn.D), insectivores (In.I, In.S and In.D), omnivores (Om.I, Om.S, and Om.D), granivores (Gf.I, Gf.S, Gf.D), and vertebrate-fish-scavengers (Vfs.I, Vfs.S); and (b) diet-foraging stratum, which comprised frugivores-nectarivores of the canopy (Fc.C), mid-height (Fn.M), generalist (Fn.Mix), understory (Fn.U), insectivores (In.C, In.M, In.Mix, In.U), omnivores (Om.C, Om.M, Om.Mix, Om.U, including Om.G as ground omnivores), granivores (Gf.C, Gf.G, Gf.Mix, Gf.U), and vertebrate-fishscavengers (Vfs.C, Vsf.G, Vsf.Mix). This was done by evaluating the correlation among species abundance and vegetation characteristics in restored and non-restored habitats, assuming a null hypothesis of the absence of correlation.
A non-metric multidimensional scaling (NMDS) analysis with Bray-Curtis distance index was conducted to evaluate the composition of each bird functional group weighted by species abundances. The stress values of the best ordination solution arrived at after 20 tries was examined and the number of dimensions increased to achieve solutions with stress values of < 20.0, due to the acceptable representation of multivariate relationships (McCune & Grace 2002). The same number of dimensions was adopted across all ordinations to allow comparisons. Thus, the vectorfit function in the Vegan package (Oksanen et al. 2019) was implemented to regress each variable on scores for the NMDS axes, calculate the strength of the association (R²), and determine the significance of this R² using a permutation test with 1000 simulations.
For graphical representation, only significant variables (p-value < 0.05) were added to NMDS ordination plots as vectors with directions weighted by the regression coefficient with each axis, and the length of the vector weighted by the R² value. This procedure was able to determine which bird functional groups were associated with the vegetation characteristics presented in the restored and non-restored habitats. All analyses and graphics were performed in R software v.3.5.2 (R Core Team 2019).

Vegetation structure
Vegetation structure differed between restored and non-restored habitats (Table I; Figure 1). The first axis of the PCA (Vegetation axis 1) explained 56.7% of the variation and was positively related to tree morphorichness (Ric.Tree) and canopy cover (Canopy), and negatively associated with grass cover (Grass). The second axis (Vegetation axis 2) explained 18.4% of the variation, with diameter at breast height (DBH) and herbaceous cover (Herb) being positively and negatively related, respectively. Our prediction that the vegetation (VSD) of the two restored habitats, even having the same post-planting age, would be structurally distinct was confirmed. Vegetation structure also differed among these and the reference habitats (pasture and forest fragments). Thus, the sampling habitats will hereinafter be referred to as restoration with more structured and developed vegetation (MS); restoration with less structured and developed vegetation (LS); pasture (PA); and fragments of native vegetation (FF).
Thus, the axes of the PCA created a vegetation structure gradient with forest-fragments (FF) and more structured restoration (MS) being dominated by a higher prevalence of Ric.Tree, Canopy and DBH (with the exception of MS1, MS9, MS11, MS12, and MS13 point-count sites). On the other hand, less structured restoration (LS) and pasture (PA) habitats were dominated by high grass cover (Figure 1; with the exception of LS3, LS 4, LS7 and LS10 point-count sites).

Birds
A total of 9,163 records representing 172 bird species (20 Orders and 44 Families) were sampled (Table SI -Supplementary Material). The observed species richness (SR) corresponded to 89% and 92% of the richness estimated by Jacknnife1 (192.9) and Chao1 (187.7) respectively. Rarefaction curves suggest stabilization only for forest-fragments (Figure 2). Less structured restoration (LS) and pasture (PA) habitats had higher SR than more structured restoration (MS) and forest-fragments (FF), while restorations (both LS and MS) had the greatest bird abundance (Figure 2), although significantly different only between LS and FF ( Figure 2).
Nine of the sampled species (5%) (São Paulo 2018); they were recorded in both PA and LS, and A. aestiva also in MS. The Common Waxbill (Estrilda astrild) was the only exotic introduced species (Sick 1997) (Table SI). Most species (151 species; 88%) were classified as residents, while

Bird functional traits and composition in relation to VSD
Only LS and FF differed significantly in accumulated SR, due to non-overlapping confidence interval curves ( Figure 2). Vegetation structural development (PCA axis 1) did not affect SR (p > 0.05) or functional divergence (FDiv, p > 0.05), but positively affected functional richness (FRic, p = 0.008), functional dispersion (FDis; p < 0.001), functional evenness (FEve; p = 0.01), and Rao's quadratic entropy (Rao's Q, p < 0.001) (Figure 3a-d). Point-counts in MS and FF were more related (Figure 3a-d). However, contrary to our expectation, restoration point-count sites showed two-way results for FRic: (1) low FRic with higher VSD, and (2) high FRic with low VSD ( Figure  3a). Plant traits induced non-aleatory effects on bird assemblages, while bird functional groups were significantly related to VSD (Tables II and  III). Ordination by NMDS revealed a trend for segregating bird composition between restored and non-restored habitats. Point-count sites of both pasture (PA) and forest-fragments (FF) were  (FDis); c) functional evenness; and d) Rao's quadratic entropy and its relationships at each habitat-types. Legend: triangles, less structured restoration; squares, more structured restoration; white circles, pasture; and black circles, forest fragment. Red dashed lines indicate significative multiple linear regression models (p ≤ 0.01). Species richness (SR) and functional divergence (FDiv) were non-significantly related to VSD and not showed. located at the extremes of the NMDS axes, while those of both LS and MS were located in the middle.
Canopy, Ric.tree, and DBH were strongly correlated with NMDS axis 1, in the same direction as diet-strata foraging, while Grass was correlated in the opposite direction; Herb was correlated with NMDS axis 2. Thus, the high predominance of exotic grass in PA, LS, and some point-count sites in MS retained bird assemblages typical of degraded sites, such as non-forest birds represented by groundforaging granivores, omnivores, and insectivores (Table II). The relative increase in the values for plant characteristics (Ric.tree, DBH, and Canopy) in restored habitats was translated into better-structured diet-strata bird groups (insectivores, granivores, omnivores and frugivores-nectarivores), which were mostly generalists foraging from the understory to the canopy. Only the increase in vegetation structure observed in FF and some point-count sites in MS could allow mid-story insectivores and frugivores-nectarivores (Table II, Figure 4). For diet-forest dependents, however, Ric.tree, DBH, and Herb were associated with NMDS axis 1 in the same direction, and Grass in the opposite direction; Canopy was related to NMDS axis 2 (Table III). Therefore, non-forest insectivores, omnivores, granivores, and carnivores were more correlated with the massive presence of Grass cover in PA, LS, and in some point-count sites in MS (Table III). Restored habitats added important plant structure (Ric.tree, DBH, Herb, and Canopy), supporting insectivores, frugivoresnectarivores, and forest semi-dependent granivores. Nevertheless, only forest fragments supported forest-dependent (insectivores, frugivores-nectarivores and omnivores), and forest-semi dependent (omnivores) species (Table III, Figure 5).   The present results add evidence that active restored habitats in the Atlantic Forest need more attention with planning to shade exotic grasses (Parrotta et al. 1997, Melo et al. 2007, since this may increase vegetation structural complexity and thus favor birds (Munro et al. 2011, Becker et al. 2013, Batisteli et al. 2018. Overall, the present study demonstrated that habitats with more structured vegetation had higher values for bird FRic, FDis, FEve, and Rao's Q compared to lesserstructured habitats (Figure 3a-d); however, SR and FDiv were not affected. This study also contributes evidence that FD may be a better metric than taxonomic diversity for evaluating how birds respond to environmental changes in restored ( Here, increases in vegetation structure represented replacement of grass cover by tree morphorichness and canopy cover. Indeed, increased habitat structure allows niche diversification because it offers greater resource partitioning and species coexistence by reducing niche overlap (Sitters et al. 2016), which leads to more diversified bird functional traits , Batisteli et al. 2018). The present study found components of FD (FDiv excepted) to be related to the higher vegetation structural development (VSD) of the forest fragments and the more structured restored sites, compared to pasture and the less structured restoration sites. As SR was unaffected by the VSD gradient, an increase in FD indexes suggests a more parsimonious distribution of functional traits at most point-count sites in restoration habitats and forest fragments, both having more structured vegetation (e.g., canopy cover and tree morphorichness). This was reinforced by significant results for FRic as a function of positive VSD, with these plant characteristics being important for creating a greater amount of niche space in young restorations. This, in turn, catalyzes the ability of restored habitats to become occupied by more diversified traits of birds, independent of any SR effect. For instance, as FRic is monotonically related to SR (Laliberté & Legendre 2010), shifts in FRic without significative alterations in SR suggests that habitats with higher VSD increased functional space for the regional bird species pool, reflecting the amount of niche space efficiently used by birds . Previous studies have also found a similar correlation (Batisteli et al. 2018, Sitters et al. 2016, but see Oliveira et al. 2019.
The present study demonstrated that the effects that canopy cover and tree morphorichness have on bird FRic were strong only in the forest fragments (Figure 3a). Although tree traits, such as DBH and canopy cover, have been considered central characteristics for augmenting the amount of niche space for birds (Batisteli et al. 2018), increased VSD in the present study area was not able to increase the occurrence of birds with unique functional traits. Instead, restorations presented contradictory FRic values in relation to VSD (Figure 3a). For example, low FRic was seen with high VSD and high FRic with low VSD. The former reflects the presence of resources (alpha niches) not completely used by the bird community while the latter reflects better use of niche space ). The first case is expected to occur due to asynchronous timing between habitat provision and bird species occupation (Catterall et al. 2012, Santos-Junior et al. 2016. Forest-dependent bird species slowly colonize restoration areas (Twedt et al. 2002, Catterall et al. 2012. This delay in niche occupation is likely to result in a low occurrence of forest specialist birds, leading to vague or underutilized niches. On the other hand, it suggests potential for the occupation of additional birds (generalists first - Gould & Mackey 2015, Becker et al. 2013, Santos-Junior et al. 2016). On the other hand, to explain high FRic with low VSD, we suspect that the presence of water-bodies at some of the point-count sites may have added some environmental heterogeneity and increased habitat diversification for varied functional traits. A similar effect, but in an inverse habitat order, was found to be related to the presence of vegetation surrounding bodies of water, which increased FD metrics of waterbirds in the Atlantic Forest (Arruda Almeida et al. 2018). Habitats in which vegetation structure was spatially variable was found to support the largest number of bird functional traits in Australia (Sitters et al. 2016). However, in the Amazon Forest, FRic for understory birds was found to be weakly affected by the distance of streams (Oliveira Increases in FDis (which is independent of SR) may be related to a decay in the abundance and/or local extinction of species having more central trait values (Laliberté & Legendre 2010). In fragmented landscapes of the Atlantic Forest, FDis of specialist birds has been strongly linked to continuous forest, whereas that of generalists has been linked to small forest fragments (Anjos et al. 2019). The data of the present study showed that even generalist birds were favored by increased VSD in both the restored sites and forest fragments, giving rise to increased niche complementarity (Mason et al. 2013). Thus, canopy cover and tree morphorichness were responsible for gains in FDis and Rao's Q. Canopy cover was also an important tree trait for increases in FDis and Rao's Q for birds in a restored riparian forest (Batisteli et al. 2018). Shifts in FDis are likely linked to turnover in abundance between open-area and forest semi-dependent species, related to shifts in the vegetation structure. This turnover of generalist versus specialist functional traits was also observed in active restored riparian forests (Batisteli et al. 2018) and in passively regenerated forests (Dias et al. 2015) in the Atlantic Forest, confirming what had previously been found for taxonomic diversity in this biome (Becker et al. 2013, Vogel et al. 2015, Santos-Junior et al. 2016. Moreover, the present results give more evidence that VSD of young restorations affects the abundance and distribution of bird traits. High FEve indicates a better distribution of biomass in niche space, which represents the effective use of resources and increasing resilience against biological invasions due to the occupancy of niches . Thus, habitats with more structured vegetation (e.g., MS and FF) assemble birds with diversified functional traits related to sensitivity, biomass, diet feeding, and diet-foraging strata. Thus, bird species composition in restored habitats may gradually become similar to that of older native forest remnants (Gardali et al. 2006, Munro et al. 2011, Catterall et al. 2012, despite the absence of forest specialized species in early restored systems (Casas et al. 2016, Santos-Junior et al. 2016.
We The shifts in avian composition found in the preset study reflect better grouping of more diverse functional traits from the regional species pool. Indeed, restored and reference habitats were unable to attract rarer/endemic and forest-dependent species. The absence of suitable microhabitats can act as a local biological filter for specialist forest birds (Santos-Junior et al. 2016). The establishment of additional forest bird species may be limited because the small native forest fragments are highly isolated, surrounded by pasturelands and croplands, with more than 5 km from the nearest large forest fragment. This scenario is highly unfavorable for bird conservation (Baum et al. 2004, Martensen et al. 2012, Barbosa et al. 2017), but it is the reality of the greater part of this biome (Ribeiro et al. 2009). Endemic Atlantic Forest birds are unable to disperse long distances across a deforested matrix, which restricts their ability to colonize isolated remnants (Uezu et al. 2005, Martensen et al. 2012). There is evidence that a minimum of 30% native vegetation may increase the occurrence of forest specialist birds (Banks-Leite et al. 2014, Boscolo & Metzger 2011.
The present study provided additional evidence that specialized functional groups are dispersed according to vegetation structure in restored sites (see also Batisteli et al. 2018). Overall, independent forest birds represented by either diet-generalists (omnivores, granivores and insectivores) that forage on the ground of open areas, or foraging strata generalists, were related to habitats with a dominance of exotic grass cover. On the other hand, semidependent and dependent forest birds with some diet-foraging strata specialization (mainly insectivores and frugivores-nectarivores that forage in the canopy, mid-height or understory), were positively influenced by tree traits (e.g., DBH, Rich.tree, Canopy and Herb), which increased linearly from restorations to forest fragments, thus representing a gradient of vegetation structure. Previous studies on both taxonomic (Twedt et al. 2002, Becker et al. 2013, Casas et al. 2016) and functional diversity (Batisteli et al. 2018) have also reported that generalist functional groups were positively related to the dominance of exotic grass cover or weeds, and birds with some habitat or food specialization were immediately recruited to restored sites (Santos-Junior et al. 2016). This recruitment occurs when the dominance of exotic grass is broken by shade provided by canopy cover and other tree traits (Melo et al. 2007), which offer additional microhabitats, food resources, and sites for perching and nesting (Catterall et al. 2012, Gould & Mackey 2015, Santos-Junior et al. 2016, Batisteli et al. 2018, and leads to a rapid increase in the abundance of birds with some diet, habitat or foraging specialization. Habitats with more structured tree traits than pasture and/or open agroecosystems provide multiple and varied microhabitats and resources for insectivorous and frugivorous species (Becker et al. 2013, Godoi et al. 2016, Sekercioglu 2012. Accordingly, pasturelands with more tree and shrub cover increase the occurrence of several forest bird species, including insectivores and frugivores (Godoi et al. 2017). Indeed, the return of insectivorous and frugivorous-nectarivorous birds is fundamental to increasing important ecosystem functions provided by birds. Insectivorous birds exert intensive biological control of herbivorous insect populations (Nyffeler et al. 2018), and thus reduce leaf damage and favor plant growth (Marquis & Whelan 1994). Frugivorousnectarivorous species act as seed dispersers and pollinators (Sekercioglu 2012), both of which are fundamental mechanisms to ensure plant reproduction and to conduct shifts in patterns of diversity and density of plant communities in restored sites (Viani et al. 2015, Casas et al. 2016, Carlo & Morales 2016, mainly those plants with fleshy fruits (Gagetti et al. 2016). Increased bird biomass augments excrement production, which in turn benefits restoration sites by inputting nitrogen into the system (Slavin & Shisler 1983).
Our findings validate our predictions. Habitats dominated by exotic grass sustain simplified functional bird assemblages, similar to that of pasture habitats. However, the addition of canopy cover and other tree attributes in restoration sites produced: (a) vertical structuration of habitats and environmental heterogeneity; (b) additional habitats and resources responsible for increases in bird functional groups with some forest specialization, in accordance with the habitat complexity hypothesis (MacArthur & MacArthur 1961); (c) better arrangement of bird functional traits (FEve and FDis, same without changes in SR) from the existent regional pool; and (d) the attraction of frugivorous-nectarivorous and insectivorous birds that might catalyze ecological succession (Viani et al. 2015) and thus assuring ecosystem sustainability and resilience (Sekercioglu 2006, Mouchet et al. 2010, Maure et al. 2018. The increase of tree traits in restored sites is also an important biological tool for reducing biological invasions of exotic grass cover (mainly canopy cover; Melo et al. 2007) and recovering central ecosystem services provide by native bird groups.
Therefore, we advocate that the rapid development of vegetation structural complexity in restored sites a key factor for provisioning planned additional habitat for birds with some specialization in diet, foraging and/or habitat resource (here called semi-specialist forest birds). To accelerating this process we recommend that restoration procedures should initially focus in the suppression of exotic grass cover. Rapid suppression of exotic grass can be achieved by introducing native tree species with larger and fast-growing canopies to promote more shaded habitats, with the addition of integrating shrub (Roels et al. 2019) and herbaceous plant species to build structurally diverse vertical vegetation. To more rapidly supply habitat and food resources