Canopy Leaf Traits, Basal Area, and Age Predict Functional Patterns of Regenerating Communities in Secondary Subtropical Forests

Secondary forests originate from natural regeneration after fallow (succession) or restoration. Species assembly in these communities, which can affect ecosystem functions and successional trajectories, is very unpredictable. Trait-based trajectories can shed light on the recovery of ecosystem functions and enable predictions of how the regenerating communities will change with forest age. Regeneration communities are affected by initial conditions and also by canopy structure and functional traits that alter dispersers' attractiveness and coexistence mechanisms. Here we evaluated how community functional traits change over time and tested if functional diversity and composition of the established canopy, as well as the structure of the canopy and forest age, influence the functional structure of regenerating tree communities when compared to their reference forests. For this, we calculated dissimilarity in trait composition (community-weighted means) and in functional diversity of regenerating communities of each succession/restoration stand, using the tree stratum of nearby mature forests as baseline values. Functional trait information comprises leaf, wood density, and reproductive traits from tree species. Our community data contain information from natural successional forests and restoration sites, in the South-Brazilian Atlantic Forest. Predictor variables of functional dissimilarities were forest age, canopy structural variables, canopy functional composition, and functional diversity. Results showed leaf traits (leaf dry matter content, leaf nitrogen content, leaf nitrogen-phosphorus ratio) and seed mass varying with forest age. Canopy functional composition based on leaf traits and total basal area significantly predicted multiple trait functional dissimilarity between the regeneration component of secondary forests and their reference community values. Dissimilarity increased when the canopy was composed of species with more acquisitive traits. Difference in functional diversity was only influenced by forest age. Mid-stage secondary forests showed lower functional diversity than early-stage forests. Our results indicated the importance of canopy traits on the natural regeneration of secondary subtropical forests. If functional similarity with reference forests is a desired objective in order to recover ecosystem functions through natural regeneration, leaf functional traits of canopy trees that establish or are planted in degraded areas must be considered in the successional processes.

Secondary forests originate from natural regeneration after fallow (succession) or restoration. Species assembly in these communities, which can affect ecosystem functions and successional trajectories, is very unpredictable. Trait-based trajectories can shed light on the recovery of ecosystem functions and enable predictions of how the regenerating communities will change with forest age. Regeneration communities are affected by initial conditions and also by canopy structure and functional traits that alter dispersers' attractiveness and coexistence mechanisms. Here we evaluated how community functional traits change over time and tested if functional diversity and composition of the established canopy, as well as the structure of the canopy and forest age, influence the functional structure of regenerating tree communities when compared to their reference forests. For this, we calculated dissimilarity in trait composition (community-weighted means) and in functional diversity of regenerating communities of each succession/restoration stand, using the tree stratum of nearby mature forests as baseline values. Functional trait information comprises leaf, wood density, and reproductive traits from tree species. Our community data contain information from natural successional forests and restoration sites, in the South-Brazilian Atlantic Forest. Predictor variables of functional dissimilarities were forest age, canopy structural variables, canopy functional composition, and functional diversity. Results showed leaf traits (leaf dry matter content, leaf nitrogen content, leaf nitrogen-phosphorus ratio) and seed mass varying with forest age. Canopy functional composition based on leaf traits and total basal area significantly predicted multiple trait functional dissimilarity between the regeneration component of secondary forests and their reference community values. Dissimilarity increased when the canopy was composed of species with more acquisitive traits. Difference in functional diversity was only influenced by forest age. Mid-stage secondary forests showed lower functional diversity than early-stage forests. Our results indicated the importance of canopy traits on the natural regeneration of secondary subtropical forests.

INTRODUCTION
Secondary forests are important for biodiversity conservation and ecosystem services, such as biomass stocks and water cycling, given the widespread deforestation and forest degradation processes in the tropics and subtropics . These forests are originated from natural regeneration processes after the abandonment of anthropogenic land-use types (succession) or after intentional actions to assist the recovery of ecosystems (restoration). Regeneration trajectories from both succession or restoration processes and the recovery of functioning and composition of these forests depend on ecosystem resilience, land-use history and landscape context, including many local and regional features (Holl and Aide, 2011;Crouzeilles et al., 2017;Holl et al., 2018). Even though many abiotic and biotic variables, as well as socio-economic context, are key drivers associated with the success of restoration strategies and regeneration trajectories (Pichancourt et al., 2014;Meli et al., 2017), biodiversity-based ecosystem services should be considered in restoration planning and monitoring to optimize forest restoration targets (Kollmann et al., 2016;Reij and Garrity, 2016;Rosenfield and Müller, 2020a). Ecosystem functioning and services are key goals in restoration projects, therefore biodiversity and trait-based researches have a high potential in informing the success of restoration efforts, embracing species parameters, interactions, and ecosystem functioning (Montoya et al., 2012;Rosenfield and Müller, 2019).
Successful regeneration of secondary forests should consider multifaceted parameters, including biodiversity, species composition, ecosystem functioning and services, to achieve a more holistic restoration. However, predicting secondary forest regeneration under natural conditions or even after active restoration actions is a big challenge. Even when considering the effect of distinct factors, studies from successional and restoration forests have shown that time is mandatory for tropical forest recovery (Holl et al., 2018), indicating that forest structure, species richness, and biomass are highly resilient in forests under natural regeneration processes (Martin et al., 2013;Zanini et al., 2014;Poorter et al., 2016;Meli et al., 2017;Rozendaal et al., 2019). Species composition, however, when compared to structural parameters or species richness, is the most unpredictable factor along succession, taking centuries to recover or even not recovering to similar pre-disturbances patterns at all (Chazdon, 2008;Liebsch et al., 2008;Rozendaal et al., 2019). Considering this and the importance of species characteristics and assemblies to ecosystem functions (Diaz et al., 2007), understanding and predicting the recovery of plant functional traits along the successional process of secondary forests is of most importance. A central question in recovering these forests is if secondary forests can achieve similar ecosystem functionality regardless of reaching similar values of former species composition. This is especially important in tropical and subtropical forests, where the number of rare species (with potential distinct and unique traits) is high when compared to the amount of common, abundant species.
Species establishment is a well-known barrier in early regeneration stages and, in forest restoration projects, planting seedlings is widely applied as a restoration tool (Holl and Brancalion, 2020). The choice of species to be used in these plantings can be a challenge, as they have distinct chances to survive and might further influence regeneration trajectories by altering the conditions for secondary species establishment or by determining coexistence mechanisms (Muscarella et al., 2017;Fernandes Neto et al., 2019;Charles, 2020). In areas of natural succession, earlier stages of forest under moist and wet conditions often have species that present a higher investment in leaf traits with high potential for acquiring resources, such as light and soil nutrients (e.g., high specific leaf area and leaf nitrogen content), and consequently grow faster (e.g., low wood density), quickly occupying the available space (Lohbeck et al., 2013;Boukili and Chazdon, 2017). However, under dry conditions, forests in earlier stages of succession, which are under more severe environmental conditions, are often dominated by species with conservative characteristics, showing denser wood stems and thicker leaves (Lohbeck et al., 2013;Poorter et al., 2019). Even though these general patterns are predictable to one extent, successional trajectories in the long-term might follow stochastic patterns or converge toward the functional composition of reference forests (Norden et al., 2015), which depends on the trait or function being analyzed (Boukili and Chazdon, 2017). Convergent trajectories regarding functional trait composition during forest development seem to prevail, but there are only few studies and less information about functional diversity trajectories along the successional time (Teixeira et al., 2020).
In addition to time and site conditions that influence community recovery trajectories, plant traits associated with the initial colonizers might influence the ongoing process (Mesquita et al., 2001), potentially causing priority effects (Fukami, 2015). Under natural regeneration, initial colonizers are influenced by land-use legacy and drive successional trajectories: colonizers of intensively used sites (e.g., pasture) arrest succession and trajectories might be unpredictable when compared to regeneration trajectories on sites that were clearcut and abandoned without any type of soil-use (Mesquita et al., 2015;Fernandes Neto et al., 2019). Besides purely structural features of the canopy (e.g., height and cover), which affect microclimatic conditions associated with light and soil (Chazdon, 2008), morphological, physiological, and phenological traits of initial colonizers might influence the ongoing regeneration community, altering dispersal attractiveness (Viani et al., 2015) and coexistence mechanisms associated to growth and survival. Under the lens of priority effects, these mechanisms are associated with niche preemption and niche modification (Fukami, 2015). Some initial colonizers have traits that determine their own success in establishing and becoming dominant across time, also slowing down the colonization by other species with different traits (Jakovac et al., 2014;Fernandes Neto et al., 2019). This may result in stable states that retard the successional process toward reference conditions, whereas other traits can lead to community attributes that improve colonization by later species or ecosystem functions, thus speeding the successional process (Weidlich et al., 2018;Fernandes Neto et al., 2019). Forest regeneration success (trajectory and speed) might thus be associated with characteristics of first colonizers along the successional process (e.g., Werden et al., 2018), besides the already known influence of structural features, landscape, and abiotic conditions. Beyond the potential priority effects from first colonizers of community reassembly, the established forest canopy and their characteristics might continuously affect the regeneration mechanisms within forest dynamics (Muscarella et al., 2017).
By understanding the influence of some key traits of initial colonizers on forest regeneration patterns, relevant restoration actions can be employed. For example, by choosing tree species for planting through traits that may best influence subsequent regeneration trajectories, one could enhance the chance of restoring ecosystem functioning in a similar way as found in pre-disturbance communities, i.e., the target forest communities (Laughlin, 2014;Laughlin et al., 2018). Here we aim to assess if functional trait composition and diversity of regenerating tree communities in successional forests increase their similarity in relation to reference forests over time, also considering the characteristics of the established canopy as potential drivers of observed trajectories. More specifically, we are interested in evaluating (i) how community functional traits change with the age of secondary forests; and whether (ii) the functional diversity and composition of the established canopy, the structure of the canopy, or the time since forest abandonment influence the similarity of regenerating tree communities in relation to the functional structure of reference forests. This similarity was evaluated considering trait composition and functional diversity differences between each successional forest and their nearby reference forest. We are interested in possibly making predictions of trait-based trajectories of regenerating tree communities based on the leaf, reproductive, and wood characteristics of initial colonizers or planted trees (which are the main components of the established canopy), in order to recover reference baseline values of ecosystem functioning.

Study Region
The study was conducted in the southern Brazilian Atlantic Forest domain. The most dominant climate in the region is classified as Cfa temperate humid climate, lacking a marked dry season, with hot summers above 22 • C (Peel et al., 2007). In the eastern highlands, the climate type is Cfb, with the average temperature of the warmest month not exceeding 22 • C, and the annual isotherm below 18 • C. The southern Atlantic Forest encompasses different forest types such as the Atlantic moist forest, Araucaria forest (a mixed subtropical forest type located in the highlands), and seasonal forests, with high levels of beta diversity (Bergamin et al., 2017).
We included seven study sites from which we collected data on the regeneration tree communities in secondary forests and on reference tree communities (one for each site) in nearby mature forests. Secondary forests of four of our study sites consist of stands under natural succession after land abandonment (Zanini et al., 2014;Vicente-Silva et al., 2016), whereas the remaining three are sites undergoing ecological restoration after planting interventions of tree seedlings (Rosenfield and Müller, 2019) ( Table 1). Planting of seedlings was conducted under legal requirement, using native species with distinct number of species and density (see Rosenfield and Müller, 2019). Forest stand age of restoration sites was determined by considering the regeneration time after the end of interventions, whereas for the stands under natural succession, the information was gathered through interviews and crosschecked with satellite images (Zanini et al., 2014;Vicente-Silva et al., 2016). Landuse history of the studied secondary forests was not intensive in terms of mechanization: Maquiné areas were used for slash-andburn small subsistence agriculture; Pró-Mata areas were used for subsistence agriculture and pasture, and the restoration areas were formerly used for cattle grazing (Cachoeirinha), Eucalyptus plantation (Canela), and grape production (Santa Tereza). The age range of secondary forests, with the respective forest type and other relevant information of the study areas, are shown in Table 1. For each study site, we identified and surveyed a nearby reference site (old-growth or mature forest) to be used as baseline values of a desired tree functional structure (trait composition and functional diversity). For the reference forests, we set an approximated age since the interruption of human interventions.

Data Sampling
In each study site, we surveyed tree and treelet species in both successional/restoration and reference forests. Reference forests were described by the tree stratum, including information of trees with diameter at breast height (DBH) ≥ 5 cm, in order to represent adult tree communities of mature forests. Secondary forests (successional/restoration stands) were described by the canopy (trees of the upper stratum, including all individuals with DBH ≥ 10 cm) and the regeneration component. The upper stratum of trees in secondary forests was defined as canopy, although we recognize that in very initial successional forests the crown of such individuals might not yet form a closed canopy. Structural characteristics of this stratum were also evaluated, including tree height, abundance and total basal area. By describing these features we intend to explicitly analyze the potential effect of these canopy trees (overstory) on the understory regeneration. The regeneration component was considered to be saplings from tree and treelet species with height ≥ 1.50 m and DBH < 5 cm. Sampling design differed between study sites, so there is a variation in the number of stands and plots depending on the study site (Table 1). For the tree stratum in reference areas and the canopy of secondary forests, plots were 100 m 2 in size, but the number of plots per stand varied depending on site characteristics. For the regeneration component, the sampling within each 100 m 2 plot varied from 12 to 40 m 2 : in areas where restoration interventions were conducted (Cachoeirinha, Canela, and Santa Tereza), sampling consisted of three subplots of 4 m 2 (i.e., 12 m 2 per plot), while in the successional areas (Maquiné and São Francisco de Paula), sampling consisted of four subplots of 10 m 2 in size (i.e., 40 m 2 per plot). For our analysis, sampling units consisted of each secondary forest stand with a distinct age in each study site (i.e., six sampling units in restoration sites and 25 in successional sites; Table 1). Data from reference forests were only considered as baseline values, being used to calculate the dissimilarity in functional composition and diversity with the secondary forest stands within each study site (see Statistical analysis for more details). Data on species composition was averaged for number of individuals per 100 m 2 per forest stand.

Functional Traits
We selected traits related to resource acquisition, structure, establishment, and reproduction. Resource acquisition traits are represented by leaf traits-leaf area (LA, cm 2 ), specific leaf area (SLA, cm 2 /g), leaf dry-matter content (LDMC, mg/g), leaf nitrogen and phosphorus content (LNC and LPC, respectively, %), and the ratio of leaf nitrogen and leaf phosphorus content (LNC: LPC). These traits are related to photosynthetic assimilation, leaf stem structures, and nutrient concentration (Wright et al., 2004;Poorter et al., 2008). Wood density (WD, g/cm 3 ) is related to stem structures, hydraulic conductance, and competitive ability (Poorter et al., 2008;Chave et al., 2014). Plant height is another key trait used to explore species competitive ability, however, in the present study, we did not include it as a trait because we used field-measured canopy height directly as a predictor variable in the models. Finally, reproductive traits are represented by seed mass (SM, g) and fruit size (FS, small: up to 5 mm of a mean size, medium: between 5.1 and 14 mm, large: larger than 14.1 mm (Githiru et al., 2002), which are related to seed dispersal and establishment (Pérez-Harguindeguy et al., 2013).
Leaf trait information was obtained from our Plant Ecology Lab at the Federal University of Rio Grande do Sul (LEVEG-UFRGS), which has a functional trait database based on field samples collected from individuals of the regional tree species pool that had been collected since 2009. Thus, for the present study, we consider trees measured within the region to calculate trait means per species. Collection and measurement protocols follow Pérez-Harguindeguy et al. (2013). Currently, all data is available under request from the TRY database (Kattge et al., 2020). Wood density was mostly obtained from the literature (Chave et al., 2014;Zimermann Oliveira et al., 2019), as well as tree reproductive traits (e.g., Sobral et al., 2006;Galetti et al., 2011;Seger et al., 2013). The list of species per study site and their mean trait values can be seen in Supplementary Table 6.

Statistical Analysis
The community-weighted mean of traits (CWM) was used as a functional composition measure (Garnier et al., 2004), while community functional diversity was represented by Rao's entropy index (Botta-Dukát, 2005). Before CWM and RaoQ calculation, we tested the correlation between functional traits. We excluded LPC due to its high correlation with LNC (r = 0.60) and retained the remaining traits in the analyses: LA, SLA, LDMC, LNC, LNC:LPC, WD, SM, and FS.
For objective (i), we calculated CWM values for each trait, based on a matrix of species abundances, and RaoQ for each sampling unit of secondary forests (regeneration component) and reference sites (tree stratum). We then tested the association between forest age and CWM values for each trait using simple linear regressions, standardizing trait values within each study site before fitting the models. The age of forest was log-scaled prior to the analysis.
For objective (ii), we first log-transformed values of LA, SLA, LNC, LNC:LPC, and SM, and standardized all traits to zero mean and unit variance [using the "decostand" function in the vegan package (Oksanen et al., 2019)]. Following this procedure, with the matrices of species abundance per stand and species by traits, we calculated CWM and RaoQ values for (1) the tree stratum of the reference forests, (2) the regeneration component of the successional/restoration stands, and (3) the canopy of the successional/restoration stands, using the Syncsa package in R (Debastiani, 2020).
Functional composition and diversity of the tree stratum from reference forests were used to set baseline values (target) for each site (Table 1), whereas the regeneration component of the successional/restoration stands was used to set the status of secondary forest regeneration. Based on this, we calculated the dissimilarity between the regeneration component and the respective reference baseline (i.e., within each site) to evaluate how distant (or similar) from the target tree regeneration of secondary forest stands were. For CWM values, we used the pairwise Euclidean distance between each regeneration community and their reference ("CWM.dist") to reflect the functional composition dissimilarity of tree regeneration communities in secondary forest stands with its target forest. For RaoQ values, we directly calculated the difference (RaoQ.diff) between the values found in each reference and regeneration stand: [RaoQ.diff = RaoQ reference -RaoQ regeneration ]. This parameter thus shows if regenerating tree communities are reaching similar functional diversity values (near zero) or if they have lower (positive values) or higher (negative values) diversity than their reference forests. As the values of CWM.dist and RaoQ.diff consider comparisons within each site, we were able to control for differences associated to forest type, species pool, and regeneration process (succession or restoration), and only the differences in terms of functionality (composition and diversity) were considered as response variables in the statistical analyses. We repeated this procedure three times: with all traits analyzed together, with leaf traits only and with reproductive traits plus wood density (named "reproductive traits, " hereafter). Due to high correlations (r ≥ 0.70) between some of the variables generated by this procedure (see Supplementary Table 1), we selected the following response variables to be used in the models: CWM distance for leaf and reproductive traits (CWM.dist leaf and CWM.dist rep ), and RaoQ difference for all traits (RaoQ.diff all ).
Community data of the established canopy in successional/restoration stands was used to determine the predictor variables in our analysis, which were related to functional composition (CWM), diversity (RaoQ), and structure (total basal area, abundance, canopy mean height and variance). We conducted Principal Component Analyses (PCA) with canopy CWM traits in order to synthesize the main functional composition axes associated with this component. We performed three different PCAs: one including all traits, a second using only leaf traits, and a third with only reproductive traits. We then evaluated correlations between the first two axes of each PCA to select variables (functional composition axes) to be used in the models (see below). In order to meet model analysis requirements, we conducted log-transformation on stand age and square-root transformation on basal area and then standardized all predictors to mean zero and unit variance.
For all three response variables described above (CWM.dist leaf , CWM.dist rep , RaoQ.diff all ), we fitted linear mixed models in order to evaluate the effect of predictor variables from the secondary forest canopy on regeneration trait compositional and diversity differences. Site (Table 1) was used as a random effect to control for local conditions. We first determined four main groups of predictor variables, each with a set of operational variables (Supplementary Table 2): Stand age (years after fallow or from the end of restoration activities), Canopy structure (structural variables), Functional composition (PCA axes originated from CWM values), and Functional diversity (RaoQ's values). Within each group of variables, we tested for correlations between variables (Supplementary Tables 3, 4) and only kept variables with r < 0.70. Predictor variables used in our full models were: stand age; basal area (canopy structure variables); PC1 leaf , PC2 leaf , and PC2 rep (canopy functional composition variables); and RaoQ all (canopy functional diversity variable). We additionally used the variance inflation factor (VIF) to test for multicollinearity among these variables before running the regression models of each response variable. We found no inflation among the chosen variables. We also accounted for the presence of spatial autocorrelation between our plots by fitting each model with and without an autocorrelation structure (i.e., linear, spherical, exponential, Gaussian, and rational quadratics relations based on local X and Y coordinates), and then selecting the model with the lowest AICc. This diagnosis showed no significant spatial structure regarding our response variables. Finally, regressions were conducted using the full model for each response variable and then running a backward selection (using the "stepAIC" function). The optimal model was considered the one with the lowest AICc (Burnham and Anderson, 2002). Variables considered significant within each optimal model were analyzed based on their confidence intervals. We performed all statistical analyses in R v.3.6.1 (R Core Team, 2019).

RESULTS
The associations between forest age and functional composition or diversity showed a high range of variation in tree regeneration in secondary forests, even between forests with similar ages (Figure 1). Significant relationships between forest age and functional composition values were observed for leaf dry-matter content and seed mass (positive association) and leaf nitrogen content (negative association) (Figure 1). The ratio between nitrogen and phosphorus (LNP:LPC) also increased with stand age (with a marginally significant p-value of p = 0.057).
Principal Component Analyses indicated functional composition patterns from the canopy of secondary forests and synthesize the differences according to leaf and reproductive traits (Figure 2). For leaf traits, the first PCA axis (PC1 leaf ) showed a variation from communities where trees have mostly acquisitive traits-higher LA, LNC and SLA (positive scores)-to communities with higher dominance of trees with conservative traits and small leaves, especially in restoration sites (warm colors), but also in some successional stands (Figure 2A). The second leaf trait axis (PC2 leaf ) was associated with tough leaves (higher LDMC) and higher limitation of nitrogen (LNC:LPC). For reproductive traits and wood density, large fruit size, and denser wood trees separated communities along the first PCA axis (PC1 rep ), whereas the second axis (PC2 rep ) represented a gradient from larger seed mass values (positive scores) to smaller fruits ( Figure 2B). Axes PC1 leaf , PC2 leaf , and PC2 rep constitute predictor variables used in the models that evaluated drivers of community regeneration dissimilarities with reference targets (see below).
Results from linear mixed models indicated three predictors for CWM.dist leaf , two for CWM.dist rep , and four for RaoQ.diff all in their respective optimal models (Figure 3, Supplementary Table 5). The optimal model for CWM.dist leaf showed a positive influence of PC1 leaf (which represents an acquisitive strategy gradient), and a negative effect of basal area ( Figure 4B) on the functional distance between secondary forest stands and their reference forests ( Figure 4A). This indicates that the regeneration component of secondary forests is functionally more similar to their reference baseline when the canopy trees are composed by species with more conservative traits (lower LA, LNC, and SLA) and under forest stands with higher total basal area of canopy trees. The selected optimal model for CWM.dist rep showed a negative influence of basal area on the functional composition distance (Figure 4C), evidencing a regeneration component functionally more similar to reference forests in communities characterized by larger trees on the upper stratum. The optimal model for functional diversity (RaoQ.diff all ) indicated a significant positive association with stand age (Figures 3, 4D). Values of RaoQ.diff all close to zero indicate similar functional diversity between the regeneration component and their reference forest, whereas positive values indicate lower diversity in the regeneration component. Thus, regenerating tree communities of older secondary forests showed lower functional diversity than the observed in their reference forests and early secondary forests.

DISCUSSION
In this study, we show that functional leaf trait composition and total basal area of the established canopy, as well as forest age, determine the functional patterns of natural tree regeneration in the understory of secondary subtropical forests (ranging from 6 to 45 years old). Leaf attributes of the upper stratum trees in early to intermediate secondary forests, mainly composed by early colonizers in successional forests and/or planted trees in restoration sites, determine how similar the regenerating trees communities are when compared to the observed functional composition patterns in reference mature forests. Community traits similarity was higher in sites where there is a higher amount of canopy tree basal area. We additionally show that the functional diversity of regenerating communities was lower in mid-successional forests (20-45 years) when compared with younger stands and with reference forests. Considering that functional similarity between regenerating communities and their reference forests is the desired target in order to recover ecosystem functions of the formerly existing forests, leaf functional traits of canopy trees that early colonize degraded areas and/or are planted for restoration must be considered as important features in the successional process, especially for restoration purposes.

Functional Composition and Diversity of Regenerating Trees Across Forest Age
Our results indicate an overall dispersion of trait communityweighted means in the successional forests, without clear trajectories along the age of successional forests for many individual traits, and also for functional diversity. However, seed mass, leaf dry-matter content, and the proportion of nitrogen and phosphorus on leaves (both LNC and LNC:LPC) showed significant relationships with forest age. Seed mass had the strongest relationship, showing a continuous increase in community mean values with forest age. This is in accordance with other papers that found a decrease in the proportion of small-seeded species with increased forest age (Warring et al., 2016). Large or heavier seeds are usually found in intermediate to advanced stages of succession, as they are usually associated with animal dispersal (Tabarelli et al., 2008), which takes longer to recover after disturbance. Leaf traits LDMC, LNC, and LNC:LPC also evidenced an expected change toward mature forest communities, showing a prevalence of more conservative attributes across successional time in moist tropical forests (Lohbeck et al., 2013;Fernandes Neto et al., 2019). Leaf nitrogen-phosphorus ratio increases with forest age, evidencing a potential nitrogen limitation in initial successional forests and a phosphorus shortage mainly observed in old-growth forests (Zeng et al., 2017). This ratio can be a useful indicator of nutrient limitation for productivity (Townsend et al., 2007). However, the transition state between nitrogen and phosphorus limitation is a ratio of 14-16, and our community values ranged from 11.7 to 15.8 (mean = 14.1), indicating almost no limitation for phosphorus, but maybe for nitrogen in some early-stage communities. For the other traits studied here (leaf FIGURE 1 | Scatter diagrams of functional trait composition (CWM: community-weighted means; A-J) and diversity (RaoQ: based on all traits; K) from the tree regeneration component of secondary forests of southern Brazilian Atlantic Forest against the age since regeneration started. Successional secondary forest sites are Encantada, Forqueta1, Forqueta2, and ProMata (green-blue colors), while restoration sites are Cachoeirinha, Canela, and StTereza (warm colors). Values for the tree stratum of reference forests are also given and correspond to the oldest forests in the right-hand side of axis x (indicated in triangles). Regression lines indicate significant associations (p < 0.05 or p = 0.057 for (E); see text for further information). LA, leaf area; SLA, specific leaf area; LDMC, leaf dry matter content; LNC, leaf nitrogen content; LNC:LPC, leaf nitrogen-phosphorus ratio; SM, seed mass; FS, fruit size (sm: small; md: medium; lg: large); WD, wood density.
Frontiers in Forests and Global Change | www.frontiersin.org FIGURE 2 | Ordination diagrams of Principal Component Analysis (PCA) for community-weighted means of leaf (A) and reproductive (B) traits of canopy trees in secondary forests of southern Brazilian Atlantic Forest. Successional secondary forest sites are Encantada, Forqueta1, Forqueta2, and ProMata (green-blue colors), while restoration sites are Cachoeirinha, Canela, and StTereza (warm colors). Circle sizes indicate an increase in the regeneration age of each forest stand. LA, leaf area; SLA, specific leaf area; LDMC, leaf dry matter content; LNC, leaf nitrogen content; LNC:LPC, leaf nitrogen-phosphorus ratio; SM, seed mass; FS, fruit size (sm: small; md: medium; lg: large); WD, wood density. The axes PC1 and PC2 of leaf traits, and PC2 of reproductive traits were considered as predictor variables of functional composition distance and diversity difference between the regeneration component of these secondary forests and their reference baseline (response variables: CWM.dist and RaoQ.diff in Figures 3, 4). FIGURE 3 | Standardized coefficients of optimal linear mixed models selected for each response variable considered for secondary forests in southern Brazilian Atlantic Forest. Response variables are functional traits composition distance (CWM.dist in terms of leaf or reproductive traits) and diversity difference (RaoQ.diff) between the regeneration component of secondary forest stands and their reference baseline. Selected predictor variables are the age of forest stand since regeneration started (Age) and parameters derived from the canopy of studied secondary forests: total basal area of canopy trees (BA), distinct spectra of strategies (PC1 leaf and PC2 leaf shown in Figure 2A, with a clear acquisitive strategy gradient in PC1), and the functional diversity based on all traits (RaoQ.all). Points are the average estimates of the models and bars correspond to 95% confident intervals (CI). Variables that do not overlap zero (indicated in black) have significant effects on the response variable.
area, SLA, fruit size, and WD), we found a large variation across the studied forests, which prevented us to evidence any trend with forest age. One could improve community assembly predictions and understanding of natural regeneration process along with forest age development, by considering intraspecific trait variability. Plant functional traits vary among individuals of the same species, especially when environmental conditions change (Albert et al., 2010;Siefert et al., 2015), but here we were not able to account for intraspecific variability within the studied species.
FIGURE 4 | Scatterplots of significant predictor variables of functional trait composition distances (CWM.dist) and functional diversity differences (RaoQ.diff) between the regeneration component of secondary forest stands and their reference baseline in southern Brazilian Atlantic Forest. Successional secondary forest sites are Encantada, Forqueta1, Forqueta2, and ProMata (green-blue colors), while restoration sites are Cachoeirinha, Canela, and StTereza (warm colors). Optimal models are presented in Figure 3 and Supplementary Table 5. In plots (A-C), values in axis y indicate the Euclidean distance between community-weighted means (CWM) of leaf traits (CWM.dist leaf ) and reproductive plus wood density traits (CWM.dist rep ) of the regeneration component of secondary forests and their respective reference baseline (i.e., comparisons were restricted to each site, see Table 1). In plot (D), values in axis y indicate functional diversity differences: RaoQ.diff all = RaoQ reference -RaoQ regeneration ; values close to zero suggest similar functional diversity between the regeneration component of secondary forest and their reference baseline. Regression lines indicate significant associations (p < 0.05; see Supplementary Table 5). Circle sizes indicate an increase in the regeneration age of each forest stand.
Time since abandonment in secondary succession and/or after restoration interventions is crucial for natural regeneration processes in tropical and subtropical forests. Tree stem density, biomass, height, and species richness can achieve similar values to old-growth forests in 30 to 50 years of forest regeneration (Guariguata and Osterag, 2001;Letcher and Chazdon, 2009;Poorter et al., 2016;Rozendaal et al., 2019), however successional trajectories can be either deterministic or stochastic. Temporal dynamic of successional trajectories varies widely within and between sites under similar land-use history, challenging the predictability even for structural features (e.g., stem density and basal area) in highly diverse tropical forests (Norden et al., 2015). This unpredictability might be associated with plot identity (i.e., involving processes that vary under local environmental conditions and interactions that occur inside each plot), which can overlap the contribution of age on successional dynamics (Norden et al., 2015). Random assembly patterns might be expected due to the high diversity of tropical and subtropical forests per se and to the often-patchy recruitment patterns of species (Fridley, 2013), especially when we think of low-abundant or rare species that predominate in these forest communities.

Determinants of Similarities in Functional Composition and Diversity Between Regenerating Trees and Reference Forests
Tree regenerating communities were more similar to reference forest values in terms of leaf traits as the canopy of established communities was composed mostly by species with conservative plant strategies: lower values of SLA, LA, and LNC. A study that experimentally tested distinct sets of planting trees (animaldispersed vs. wind-dispersed tree species), aiming to restore tropical forests in Mexico, did not find differences in phylogenetic and functional diversity of early communities of recruits (eight years old), assuming that the initial composition did not influence community assembly (Li et al., 2018). Animal-dispersed species are often pointed out as an important component of forest recovery, due to their potential in attracting dispersers from nearby forest remnants and thus increasing the proportion of propagules' arrival (Boukili and Chazdon, 2017;Meli et al., 2017); leaf trait composition, however, is less frequently considered. Here, we found consistent influence of leaf attributes of the canopy on the functional composition of tree regenerating communities. Leaf functional traits of initial colonizers and/or of planted trees affect abiotic conditions-by altering shading patterns and nutrient content in litter and on the soil surface (Poorter, 2009;Rosenfield and Müller, 2020b;Teixeira et al., 2020)-and biotic interactions associated to competition for light and the use of soil resources (Kunstler et al., 2016;Muscarella et al., 2017). These conditions are relevant for mechanisms of tree establishment and growth and might be more relevant than the dispersal/arrival of individuals for the regeneration of secondary forests toward reference baseline values. Dispersal seems not to be a barrier to forest regeneration where the surrounding landscape still have forest remnants to serve as font of propagules and habitat for dispersers (Meiners et al., 2015), and our studied areas do not seem to be limited in terms of propagules arriving from surrounding forests (Zanini et al., 2014;Vicente-Silva et al., 2016;Rosenfield and Müller, 2019).
Functional composition distances between the regeneration component and the reference forests were also associated with the total basal area of the canopy stratum, a parameter of forest structure. Secondary forest stands with a greater sum of basal area of upper trees (canopy) had a regeneration component more similar to the expected target than less structured forests (i.e., with lower values of total basal area). The similarity was greater both in terms of leaf traits and reproductive plus wood density traits, indicating a secondary forest trajectory toward the expected target in terms of the functional composition and, consequently, ecosystem functionality. The functional composition of forest communities is often associated with ecosystem processes and temporal dynamics, representing the effect of the attributes of the most abundant species (mass-ratio theory), instead of the variability of traits (niche complementarity hypothesis), on ecosystem functions (Garnier et al., 2004;Van der Sande et al., 2018;Bordin and Müller, 2019;Rosenfield and Müller, 2020b). Greater similarity with respect to functional composition in these secondary forests may indicate that important ecosystem functions associated with these traits, and encompassing different compartments of the system (e.g., nutrient cycling, trophic interactions), may be recovering as well (Rosenfield and Müller, 2020a). Less structured secondary forests often fail to offer adequate conditions, such as more shadowed environments, for the establishment of secondary tree species (Holl et al., 2018). Regenerating trees in such communities might be composed by species with distinct strategies, e.g., pioneer species, thus increasing the distance from reference forests in terms of leaf and reproductive traits.
Differences in functional diversity values of tree regeneration communities in secondary forests when compared to reference values were predicted by forest age. The best model included functional composition and diversity variables, but only age was significant and with a low predictability power. Forest age was positively associated with the differences in RaoQ, indicating that early-stages of succession (6-12 years) have more similar values to reference forests (points close to zero located in the left side of Figure 4D) than mid-stage secondary forests (15-45 years), that still showed lower values of functional diversity in the regeneration component (positive values of RaoQ.diff; right side of Figure 4D). This is partially in accordance with a study that found decreasing values of functional diversity along succession (from 10 to 40 years) in tropical forests, but higher functional diversity in old-growth forests (Lasky et al., 2014). In this study, the pattern of decreasing functional diversity along succession was found when analyzing all traits together in a single index (similar to what was performed in our study), however diversity values based on single traits showed slight increases with stand-age (Lasky et al., 2014). In our study, RaoQ for all traits was highly correlated with RaoQ calculated for leaf traits and reproductive traits separately, so the observed pattern was consistent among the sets of chosen traits and not skewed by one or another selected set (Zhu et al., 2017).
The results found for functional diversity patterns can be related to transient community stages, where pioneer and secondary tree species co-occur in the regeneration component of early-stage successional forests without the dominance of only a few species. This leads to functional diversity values to be similar to the values found in the tree stratum of reference forests. High rates of species arrival and the absence of severe abiotic or biotic filters, arresting establishment, might be contributing to high functional diversity values in the early-stage secondary forests studied here. Along successional time, some shade-tolerant tree and treelets (e.g., Cupania vernalis, Casearia silvestris, Mollinedia shottiana, Piper aduncum, Psychotria suterela) eventually become more dominant in mid-stage forests, due to high colonization rate and low mortality in the understory. This dominance affects RaoQ values, which considers pair-wise traits' distance weighted by species abundance (Botta-Dukát, 2005). As self-thinning processes occur (Stephenson et al., 2011), decreasing the dominance of these species, there should be an increase in functional diversity, increasing the similarity with reference forests. However, functional diversity does not often vary across succession (e.g., Lohbeck et al., 2012;Böhnke et al., 2014), as the turnover of species along the regeneration time can maintain similar patterns of trait dispersion within communities.

CONCLUSION
A great portion of current subtropical and tropical forest coverage comprises secondary forests, which increases the role of these forests for biodiversity conservation and the provision of ecosystem services Martin et al., 2013;Pichancourt et al., 2014). Even though many studies have shown that the recovery of forest structure parameters and even species diversity is possible in a relatively short time of succession (e.g., Guariguata and Osterag, 2001;Letcher and Chazdon, 2009;Poorter et al., 2016), species composition is hard to recover or restore (e.g., Dent et al., 2013;Suganuma et al., 2014;Rozendaal et al., 2019). Understanding whether functional composition and diversity of regenerating tree communities increase in similarity when compared to reference forests is essential in order to evaluate the potential return of ecosystem functions and services from these secondary forests. Our study offers an opportunity to understand the effect of first colonizers and/or planted species in the functionality patterns of regenerating tree communities in the understory of secondary forests originated from both passive and active restoration processes.
We provided some perspectives about the importance of forest structure and leaf trait composition of canopy species in the functional patterns of the future tree component of secondary forests (i.e., tree regeneration component). These perspectives go beyond any potential influence of climatic and soil conditions, as they concern natural regeneration of forests under sites that differ in terms of environment and species pool. Leaf trait composition of trees from the upper stratum influences the similarity of ongoing successional communities, indicating potential priority effects of initial assemblies on natural regeneration processes. Thus, despite the great potential of natural regeneration in tropical and subtropical moist forests, under the scope of ecological restoration, practitioners should consider leaf attributes that characterize plants with acquisitive strategies that grow fast to compose the set of species for early-stage communities, by planting or by thinning in the first years of recovery. This might turn initial canopy structure more similar to the observed in natural successional forests and might accelerate the increment of tree stem basal area. As we saw, forest stands with higher basal area of canopy trees can improve the functional similarity of regenerating trees with the expected functional target. In addition, the presence of species with conservative strategies in the canopy of early and intermediate successional forests seems to improve similarities of tree regeneration communities with target functional patterns; thus, species with such strategies (low SLA, LA, LNC) should also be considered in restoration plantings, especially after the planting of fast-growing species. Finally, the influence of recovery time on the regeneration of functional patterns along the successional and restoration processes is uncertain, making the monitoring of the tree regeneration component necessary to inform if, when, and where secondary forests are becoming functionally similar to their reference forest ecosystems.

DATA AVAILABILITY STATEMENT
The data analyzed in this study is subject to the following licenses/restrictions: Most data are available as Supplementary Material; additional information might be available under request to the authors. Requests to access these datasets should be directed to Sandra Cristina Müller, sandra.muller@ufrgs.br.

AUTHOR CONTRIBUTIONS
SM conceived the idea of the paper. RB, KB, JK, and MR participated in the collection of field data. SM, MR, KB, and JK performed statistical analyses. MR and KB prepared the figures. SM and MR led the writing of the manuscript. All authors contributed to the collection and compilation of functional trait data, discussed the results and thoroughly revised the manuscript, and contributed critically to the drafts and gave final approval for publication.

FUNDING
This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior-Brazil (CAPES) through the Programa de Excelência Acadêmica (PROEX) and the Portal de Periódicos, by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq grant number 309659/2019-1 to SCM) and by the Fundação de Amparo à Pesquisa do Estado do Rio Grande do Sul (FAPERGS process 2218-2551/12-2).