Forest and Grassland Ecosystem

Predictability of Functional Diversity Depends on the Number of Traits

  • ZHANG Zihao 1 ,
  • HOU Jihua , 1, * ,
  • HE Nianpeng , 2, 3, 4, *
Expand
  • 1.Key Laboratory of Forest Resources and Ecosystem Process, Beijing Forestry University, Beijing 100083, China
  • 2.Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
  • 3.College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
  • 4.Key Laboratory of Vegetation Ecology of Ministry of Education, Institute of Grassland Science, Northeast Normal University, Changchun 130024, China
*HOU Jihua, E-mail:;
HE Nianpeng, E-mail:

Received date: 2020-10-27

  Accepted date: 2021-01-04

  Online published: 2021-07-30

Supported by

The National Natural Science Foundation of China(31872683)

The National Natural Science Foundation of China(31800368)

The National Key Research and Development Program of China(2017YFA0604803)

Abstract

Analysis of functional diversity, based on plant traits and community structure, provides a promising approach for exploration of the adaptive strategies of plants and the relationship between plant traits and ecosystem functioning. However, it is unclear how the number of plant traits included influences functional diversity, and whether or not there are quantitatively dependent traits. This information is fundamental to the correct use of functional diversity metrics. Here, we measured 34 traits of 366 plant species in nine forests from the tropical to boreal zones in China. These traits were used to calculate seven functional diversity metrics: functional richness (functional attribute diversity (FAD), modified FAD (MFAD), convex hull hypervolume (FRic)), functional evenness (FEve), and functional divergence (functional divergence (FDiv), functional dispersion (FDis), quadratic entropy (RaoQ)). Functional richness metrics increased with an increase in trait number, whereas the relationships between the trait divergence indexes (FDiv and FDis) and trait number were inconsistent. Four of the seven functional diversity indexes (FAD, MFAD, FRic, and RaoQ) were comparable with those in previous studies, showing predictable trends with a change in trait number. We verified our hypothesis that the number of traits strongly influences functional diversity. The relationships between these predictable functional diversity metrics and the number of traits facilitated the development of a standard protocol to enhance comparability across different studies. These findings can support integration of functional diversity index data from different studies at the site to the regional scale, and they focus attention on the influence of quantitative selection of traits on functional diversity analysis.

Cite this article

ZHANG Zihao , HOU Jihua , HE Nianpeng . Predictability of Functional Diversity Depends on the Number of Traits[J]. Journal of Resources and Ecology, 2021 , 12(3) : 332 -345 . DOI: 10.5814/j.issn.1674-764x.2021.03.003

1 Introduction

Functional traits are core plant properties, and are closely related to the colonization, survival, growth, and mortality of plants (Cornelissen et al., 2003). In particular, these traits allow plants to optimize productivity (Craven et al., 2018; Xu et al., 2018b; Kordbacheh et al., 2019; Mahaut et al., 2020). Recently, scientists have demonstrated that community traits on per land area using community-weighted means (He et al., 2019a), scaled up from the organ or species level, could be used to interpret certain ecosystem processes and functions at larger spatial scales (He et al., 2019a). To date, trait-based approaches have been used to explore various ecological issues from the species or community level to the ecosystem level (Faucon et al., 2017; Greenwood et al., 2017; Ma et al., 2018).
There are trade-offs or synergies that exist among traits (Cornelissen et al., 2003; de Deyn et al., 2008; Reich, 2014; Li et al., 2015; He et al., 2019b), and studing multiple traits has hence become important. The classic leaf economic spectrum demonstrates the trade-off strategy of plant resources (Wright et al., 2004). The adaptation of a single trait to an environmental change may not be important as the adaptive evolution of the specific trait may lag behind environmental change, and the adaptation rather manifested as a coordinated adaptation of multiple traits. Plants adapt to the environment by adjusting multiple traits, with the complex relationships between these traits crucial for maintaining the ecosystem structure and function (Fontana et al., 2016). Furthermore, there is consensus that species diversity influences ecosystem functioning, as well as the diversity of functional traits within communities (Leps et al., 2006). The combination of functional diversity, community structure, and plant traits could help us explore how multiple traits determine ecosystem processes and functioning (Fontana et al., 2016).
Various indices have been proposed to measure plant functional diversity (Pérez-Harguindeguy et al., 2013). In general, functional diversity can be subdivided into functional richness, evenness, and divergence (Table 1). Functional richness indices can quantify the size of the niche space utilized by these measured traits per volume of multidimensional space. Functional evenness can describe the distribution of the trait space occupied by species traits. It refers to whether the distance between the traits of any two nearest species is equal in the trait space. The functional divergence describes the formation of species clusters in the character space, where higher dispersion indicates that the species clusters are distributed at the edge of the character space (Mason and Mouillot, 2005; Laliberte and Legendre, 2010; Schleuter et al., 2010) (The detailed formula is provided in the supplementary file). Through the research of functional diversity, certain diversity indexes have been developed, including functional attribute diversity (FAD) (Walker et al., 2005), convex hull hypervolume (FRic) (Laliberte and Legendre, 2010), functional evenness (FEve) (Laliberte and Legendre, 2010), functional divergence (FDiv) (Mason and Mouillot, 2005; Laliberte and Legendre, 2010; Schleuter et al., 2010), and quadratic entropy (RaoQ) (Botta-Dukát, 2005).
Functional diversity can be used to accurately assess and predict ecosystem functioning, because multiple traits collectively accomplish this (Kearney and Porter, 2006). Similarly, functional evenness and functional divergence are used to elucidate the variation in functional and adaptive mechanisms (Komac et al., 2015). Functional diversity has received extensive attention recently; and research on the relationships between functional diversity and ecosystem function has rapidly increased (Ross et al., 2017). However, dispute remains over how to scientifically calculate and use these metrics. In practice, researchers tend to select certain traits based on their own objectives or the availability of data, leading to inconsistent results across studies that make comparisons among studies difficult (Petchey and Gaston, 2002; Wright et al., 2006; Pasari et al., 2013). Knowledge about how to use a wide variety of plant traits is crucial because different traits will generate different results (Petchey and Gaston, 2006). Therefore, development of a uniform protocol is needed that could be used consistently, or that could enhance comparability across different studies. To accomplish this, it is necessary to identify how the number of selected traits influences functional diversity indexes from the site to the regional scale, and whether these influences are consistent across the different functional diversity metrics.
To explore these issues, we measured 34 traits of 366 plant species from nine forest communities along the North-South Transect of Eastern China (NSTEC), extending from the tropical to the boreal zones. Morphological, chlorophyll, stomatal, and anatomical leaf traits were measured, in addition to multiple element content, and other traits (Fig. 1). The main objectives of this study were to: 1) demonstrate the predictability of the functional diversity metrics by altering the numbers of selected traits, and 2) develop a uniform protocol for metrics that predict functional diversity consistently from the site to the regional scale for use in future studies.
Fig. 1 The opportunities and challenges associated with functional diversity
The conceptual framework of functional diversity holds promise as a way to explore relationships in natural communities. Specifically, these functional diversity are predictable (invariably, regularly, or irregularly) under different scenarios of trait number, diversity in importance of different traits, and differences in species diversity. Development of a functional diversity protocol is needed to promote the applicability of these metrics.

2 Materials and methods

2.1 Study sites

Nine typical forest communities were selected for field sampling along the NSTEC (Fig. 2), which is the 15th standard transect of the International Geosphere-Biosphere Program (Zhang and Yang, 1995; Wang et al., 2016). To minimize the effects of human disturbance, we set up these sampling plots in national nature reserves. From south to north, the forests sampled were: tropical, south-subtropical evergreen broad-leaved, subtropical evergreen broad-leaved, and north-subtropical evergreen deciduous broad-leaved mixed forests. The study area covered almost all the main types of forest vegetation in the Northern Hemisphere (He et al., 2018). Across the study sites, the annual average temperature and precipitation ranged from -4.4 to 20.9 °C and 481.6 to 2449.0 mm, respectively, with the precipitation mainly concentrating in summer, from June to August (Xu et al., 2018a).
Fig. 2 The spatial distribution of nine forest communities from the boreal to the tropical zone in China

Note: A total of 34 functional traits and 366 species were sampled across the nine forest sites (orange circles).

2.2 Field sampling and measurement of traits

Field sampling was carried out from the end of July to the beginning of August 2013. We set up four representative plots of 30 m × 40 m in each forest type to collect mature leaves of common tree species within and around the plots. The specific sampling method involved selection of healthy, disease-free, mature trees of different species. Thereafter, four branches at different levels and orientations, specifically from four different directions in the upper, southeast, and northwest of the canopy, were collected from each tree using high-pruning shears or artificial tree-climbing methods. A mixture of leaf samples was obtained from the cut branches, placed in plastic bags, and transported to a laboratory for processing. One tree in each plot represented one replicate; consequently, four trees were considered as replicates per site (Li et al., 2018).
2.2.1 Selection of metrics for plant traits
We measured 34 plant traits of 366 species in the nine forest communities, using a consistent method for recording specific traits across all forest types (cold-temperate to tropical) (He et al., 2020; Table S1 and Fig. S2).
2.2.2 Measurement of leaf morphological traits
Six to ten healthy, intact leaves were flattened out and scanned using a Cano Scan LiDE 110 (Japan), before determining their area using ImageJ software (Nanes, 2015). The leaves were dried to a constant weight in an oven and were weighed to calculate the specific leaf area based on the leaf area and dry weight (Wang et al., 2016).
2.2.3 Measurement of C and N content
An elemental analyzer method (vario MAX CN Elemental Analyzer, Elementar, Germany) was used to determine foliar N content. The leaves were first cleaned to remove soil and other impurities and then dried. The N content in 0.3-0.5 mg of the dried samples was determined (Zhao et al., 2014).
2.2.4 Measurement of chlorophyll traits
Fresh leaf samples (0.1 g) were cut into pieces and, together, with 10 mL of 95% ethanol, were homogenized in a mortar. An additional 5 mL of 95% ethanol was added to the homogenate, which was then filtered, and the filtrate diluted to 50 mL with ethanol. The chlorophyll-ethanol solution was injected into a cuvette with a light path of 1 cm. Pure ethanol was injected into another cuvette of the same specification to serve as a control. The chlorophyll content was calculated according to Lambert Beer's law (Li et al., 2018).
2.2.5 Measurement of stomatal traits
From the samples, five to ten leaves were cut into 1 cm x 0.5 cm sections along the main vein and fixed in FAA solution (75% alcohol : formalin : glacial acetic acid : glycerol = 90 : 5 : 5 : 5) (Field sampling was completed). These sections were dried, pasted to the sample stage, and observed using a scanning electron microscope (Hitachi s-3400 II, Hitachi, Japan). Three small blocks were randomly selected from the sample stage for observation. Two photographs were taken in each of the blocks, resulting in six replicate images for each species. The number of stomata in each image was counted and the area of each image was determined to calculate stomatal density. Five pores in each image were randomly selected to determine the average stomatal pore width, and length, and stomatal area (Liu et al., 2018).
2.2.6 Measurement of anatomical traits
Samples were removed from the FAA fixative and permanent cross-section slides were made using paraffin sectioning, with a total of three replicate sections per species. Anatomical data was acquired using scanning electron microscopy (Liu et al., 2019).

2.3 Selection of metrics for functional diversity

Functional diversity is used to reflect the overall difference or diversity of traits in plant communities. At the beginning of the formation of the functional diversity, the functional diversity index should generally meet the following criteria: 1) multiple traits must be processed simultaneously, 2) species traits and species richness influence the functional diversity of the community, and 3) increasing or decreasing the number of new species or traits affects functional diversity (Villéger et al., 2008). These criteria reinforce the fact that functional diversity requires consideration of both the number of species and traits. Although many types of functional diversity exist, we selected three kinds of diversity corresponding to the three components of functional diversity: functional richness, functional evenness, and functional divergence (Table 1).
Table 1 The protocol of the functional diversity index
Functional diversity Specific indexes Description Predictable Fitted equation
Functional richness Functional Attribute Diversity (FAD) The sum of the distances of species in trait space Yes Y=513.356x0.5224†
Modified FAD (MFAD) Modified FAD, includes species diversity Yes Y=12.4681x0.5224
Convex hull hypervolume (FRic) Convex hull hypervolume Yes Y=0.6849x+5.3610
Functional evenness Functional Evenness (FEve) Distribution rule of trait space occupied by traits No
Functional divergence Functional Divergence (FDiv) Dispersion of functional traits No
Functional Dispersion (FDis) Dispersion of functional traits No
Rao's Quadratic entropy (RaoQ) Both trait richness and trait dispersion Yes Y=0.7969x-0.2641

Note: These prediction equations were reduced on all data of the nine forest communities from the tropical to the boreal zone. Detail information for each forest is presented in the supplementary files.

Functional richness is mainly used to quantify how much niche space is occupied by all species in the community (Laliberte and Legendre, 2010). Low trait richness indicates that the available niche in the community is occupied, meaning that some of the resources in the community are not utilized, resulting in lower productivity. In this study, we used the following functional richness indexes: functional attribute diversity (FAD), modified functional attribute diversity (MFAD) (Schmera et al., 2009), and convex hull hypervolume (FRic) (Villéger et al., 2008).
The evenness of traits (FEve) is the distribution law of the mean value of species traits in the trait space. This index quantifies the degree of trait uniformity. Low trait evenness indicates that some resources are over-utilized, whereas other resources are rarely utilized or have not yet been utilized.
Trait dispersion can be used to describe the degree of convergence of traits in the functional space. High functional dispersion indicates that the species group is distributed at the edge of the trait space. This metrics can be used to indicate the degree to which community resources differ as well as the degree of competition. We selected the following three trait dispersion indexes for used in our study: functional divergence (FDiv), functional dispersion (FDis), and Rao's Quadratic entropy (RaoQ).

2.4 Statistical analysis

The seven functional diversity metrics mentioned above (specific index in Table 1) were calculated using R language and FDiversiry software (Casanoves et al., 2011). Based on the well-matched trait and species data collected, we developed a calculation scheme to explore our objectives. The importance of different traits was attenuated by random selection and repeated simulation. At the same time, the multi-site simulated sampling ensured consistency of species diversity. In this way, we could focus more on the effect of trait number on functional diversity. In practice, the procedure was as follows: 1) according to the calculated metric of functional diversity, corresponding data were consistently compiled; 2) different sampling scenarios were simulated and specific trait indicators were extracted from the data table; 3) R language was used to randomly select trait indicators according to a specific number, and these metrics were calculated; and 4) the selected number of plant traits was increased one at a time from 3 to 34, with these traits being randomly selected 100 times, and the average was calculated as the result of each cycle.

3 Results and discussion

Altering the number of traits had different effects on the metrics of functional diversity. The specific indexes FAD, MFAD, FRic, and RaoQ were predictable using the equations indicated in Fig. 3 and Table 1. Theoretically, the niche space of functional richness quantification depends on the traits of species (Fig. 3). The more functional traits included, the larger the niche space that can be quantified (Mason et al., 2005). Therefore, the trait richness metrics could be measured and used to predict the total resource or niche space of different communities (de la Riva et al., 2018; Solefack et al., 2018)( Figs. S1-S4). Using a different number of traits to calculate functional diversity should yield Four of the seven functional diversity metrics included in this study were predictable, namely functional attribute diversity (FAD), modified FAD (MFAD), convex hull hypervolume (FRic), and Rao's quadratic entropy (RaoQ). When the number of traits increased from 3 to 34, 100 random combinations were calculated on average (The supplementary files provide these equations for the different regions).
Fig. 3 Relationship between functional diversity metrics and trait number different results.
The three metrics of trait richness and RaoQ increased with the number of selected traits. Specifically, FRic showed a linear increasing trend, and FAD and MFAD were well fitted as power functions (Table 1 and Fig. 3). Furthermore, RaoQ was linearly correlated with the number of selected traits, similar to FRic. A plausible explanation for these results is that the calculated equation (Eqs.1-16) contained information on both functional richness and trait dispersion; however functional richness was dominant. Furthermore, the general equations used to determine FAD, MFAD, FRic, and RaoQ were similar across the forest communities sampled (from the tropical to the boreal zone); however, some data sets differed among these functional richness metrics (Tables S2-S5).
The number of traits is a non-negligible factor for the study functional diversity; however, there is no implementation standard that can be referenced at present. Consequently, small numbers of traits (three or four) have been used to calculate functional diversity in previous studies (Mensah et al., 2016; Staples et al., 2019). Researchers should design experimental protocols for their own experimental purposes more rigorously and scientifically. Some researchers may consider functional redundancy and abandon the use of larger number of functional traits in their analyses (Correia et al., 2018; Lozanovska et al., 2018; Kearsley et al., 2019). Thus, it is important to determine which functional traits perform the more important roles in the community, along with identifying the importance of the traits (Garcia et al., 2015). In natural communities, functional redundancy must exist, and could facilitate the stability and anti-interference ability of the community (Pillar et al., 2013). In this study, the predictive equations of the three trait richness metrics (FAD, MFAD, and FRic) and RaoQ increased non-linearly with the number of traits. This may be related to functional redundancy between traits, and requires further study. Furthermore, researchers are encouraged to select more traits to explore ecosystem processes and functions, to reduce potential errors, and facilitate comparisons across studies in the future.
Using a database of consistently measured traits of plant species found in a range of forest types, we explored the predictability of functional diversity with respect to trait richness, trait evenness, and trait divergence. Our results demonstrated that only the three trait richness metrics and RaoQ were predictable, in relation to the influence of the number of traits on functional diversity. The relationship between each of the three functional richness metrics and number of traits, as well as the relationship between entropy and trait number, are expected to provide standard protocols to enable incorporation of different studies from the site to the regional scale in future. However, at present, uniform standards on what number of traits influences other functional diversity indexes do not exist, nor is there information about how this is achieved. More research is required to improve both the theory of, and methodology used to calculate functional diversity, even though the conceptual framework of functional diversity is promising.

Acknowledgements

The data of this research have been shared from Functional Trait database of terrestrial ecosystems in China (China_Traits). For further data sharing contact HE Nianpeng (henp@igsnrr.ac.cn).

Appendix Ⅰ

Fig. S1 The relationship between convex hull hypervolume (FRic) and the selected number of traits in different typical forests

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DongLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing.

Fig. S2 The relationship between quadratic entropy (RaoQ) and the number of traits in different typical forests

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DongLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing.

Fig. S3 The relationships between functional Attribute Diversity (FAD)and the number of traits in different typical forests

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DongLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing.

Fig. S4 The relationships between modified Functional Attribute Diversity (MFAD) and the number of traits in different typical forests

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DongLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing.

Fig. S5 The relationships between functional evenness (FEve) and the number of traits in different typical forests

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DongLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing.

Fig. S6 The relationship between functional divergence (FDiv) and the number of traits in different typical forests

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DongLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing.

Fig. S7 The relationship between functional dispersion (FDis) and the number of traits in different typical forests

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DongLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing.

Appendix Ⅱ

Table S1 The list of these selected traits and their abbreviations
Category Traits Unit Abbreviation Category Traits Unit Abbreviation



Morphological traits
Leaf area cm2 LA


Stomatal traits
Stomatal pore length μm PL
Leaf dry weight g LDW Stomatal pore width μm PW
Leaf water content % LWC Stomatal area mm2 SA
Specific leaf area mm2 mg-1 SLA Stomatal length μm SL
Leaf thickness mm LT Stomatal density pores per mm2 SD





Stoichiometric
Leaf carbon content % LCC Stomatal area fraction % P
Leaf nitrogen content % N



Anatomical traits
upper epidermal cells width μm UEW
Ratio of C:N NA C/N Palisade tissue thickness (PT) μm PT
Leaf Ferrum content mg g-1 Fe Sponge tissue thickness (ST) μm ST
Leaf Kalium content mg g-1 K Ratio of PT:ST NA PT/ST
Leaf Magnesium content mg g-1 Mg lower epidermal cells width μm LEW
Leaf Phosphorus content mg g-1 P vessels width μm VW
Ratio of C:P NA C/P


Other traits
Soluble sugar content (SSC) mg g-1 SSC


Chlorophyll traits
Leaf Chlorophyll a content mg g-1 Chl a Starch content (SC) mg g-1 SC
Leaf Chlorophyll b content mg g-1 Chl b Non-structure carbohydrate mg g-1 NSC
Total chlorophyll content mg g-1 Chl Ratio of SSC:SC NA SSC/SC
Ratio of Chla: Chlb NA Chl a/b Leaf calorific value kJ cm-2 LCV
Table S2 The predictive equations of trait richness index (FRic) prediction equations at different sites (mean ±SE)
Area Site Fitted equation (y = kx + b) k R2
Temperate zone HZ y = 0.038x + 0.8043 0.038±0 a‡ 0.954
LS y = 0.8472x + 4.5882 0.8472±0.0587 b 0.996
CB y = 0.7285x + 4.8064 0.7285±0.2366 b 0.9959
Subtropics DL y = 0.2274x + 2.8132 0.2274±0.0925 a 0.9986
TY y = 0.253x + 2.0453 0.253±0.0956 a 0.9918
SN y = 1.5507x + 7.3856 1.5507±0.1718 c 0.9956
Tropic JL y = 1.2865x + 8.2721 1.2865±0.2096 c 0.9989
DH y = 0.8329x + 8.2909 0.8329±0.5388 b 0.9902
JF y = 0.7927x + 12.1820 0.7927±0.1592 b 0.9849

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DangLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing. There is no significant difference between data with the same alphabet in each column (P < 0.01).

Table S3 The predictive equations of quadratic entropy (RaoQ) at different sites (mean ±SE)
Area Site Fitted equation (y = kx + b) k R2
Temperate zone HZ y = 0.5036x + 0.2167 0.5036±0 ab 0.9996
LS y = 1.1944x - 0.8658 1.1944±0.1098 f 0.9997
CB y = 1.141x - 0.8301 1.141±0.2149 f 0.9997
Subtropics DL y = 0.4314x - 0.0682 0.4314±0.0366 a 0.9998
TY y = 0.9301x - 0.7576 0.9301±0.0862 e 0.9997
SN y = 0.8623x + 0.2462 0.8623±0.1035 de 0.9998
Tropic JL y = 0.6646x - 0.1202 0.6646±0.1134 bc 0.9999
DH y = 0.699x - 0.0952 0.699±0.0609 cd 0.9999
JF y = 0.7051x + 0.0154 0.7051±0.0433 cd 0.9999

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DangLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing. There is no significant difference between data with the same alphabet in each column (P < 0.01).

Table S4 The predictive equations of functional attribute diversity (FAD) at different sites (mean ± SE).
Area Site Fitted equation (y = Axk) k R2
Temperate zone HZ y = 4.1931x0.5024 0.5024±0 a 0.9999
LS y = 62.034x0.5281 0.5281±0.0075 bc 0.9995
CB y = 94.622x0.5366 0.5366±0.0063 c 0.9997
Subtropics DL y = 43.469x0.5172 0.5172±0.0023 b 0.9997
TY y = 28.148x0.5326 0.5326±0.0163 b 0.9992
SN y = 535.82x0.5195 0.5195±0.0050 b 0.9998
Tropic JL y = 1199x0.5251 0.5251±0.0020 bc 0.9997
DH y = 345.39x0.5232 0.5232±0.0018 bc 0.9999
JF y = 3141.7x0.5339 0.5339±0.0015 c 0.9996

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DangLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing.

There is no significant difference between data with the same alphabet in each column (P < 0.01).

Table S5 The predictive equations of modified functional attribute diversity (MFAD) at different sites (mean ±SE)
Area Site Fitted equation (y = Axk) k R2
Temperate zone HZ y = 1.3977x0.5024 0.5024±0 a 0.9999
LS y = 5.894x0.5281 0.5281±0.0075 bc 0.9995
CB y = 7.8935x0.5315 0.5315±0.0063 c 0.9998
Subtropics DL y = 4.9451x0.5168 0.5168±0.0023 b 0.9997
TY y = 3.426x0.5284 0.5284±0.0163 b 0.9993
SN y = 17.892x0.5199 0.5199±0.0050 b 0.9998
Tropic JL y = 26.969x0.5251 0.5251±0.0020 bc 0.9997
DH y = 14.092x0.5233 0.5233±0.0018 bc 0.9999
JF y = 43.226x0.5338 0.5338±0.0015 c 0.9996

Note: HZ, HuZhong; LS, LiangShui; CB, ChangBai; DL, DangLing; TY, TaiYue; SN, ShenNong; JL, JiuLian; DH, DingHu; JF, JianFengLing.

There is no significant difference between data with the same alphabet in each column (P < 0.01).

Appendix Ⅲ

The formula of trait diversity index
1. FAD is commonly calculated as Eq. 1 and Eq. 2:
$E{{D}_{ij}}=\sqrt{\underset{t=1}{\overset{T}{\mathop \sum }}\,{{\left( {{x}_{tj}}-{{x}_{ti}} \right)}^{2}}}$
$FAD=\sum{i=1}{S} \sum{j>1}{S} ED_{ij}$
where T is the number of traits, and xti and xtj are the values of species i and species j for trait t, respectively. EDij is the Euclidean distance between species i and species j, and S is the number of species.
2. MFAD is a multidimensional trait richness index derived from the development of the FAD index, which considers the impact of species diversity (Eq. 3).
$MFAD=\frac{FAD}{S}$
3. FRic is calculated by the area or volume of the smallest polygon generated within the multidimensional convex hull volume formed by multiple traits (Villeger et al. 2008). This is a complex calculation that can be done in both Qhull software (http://www.Pricklysoft.org/software/raithull.html) and R languages (http://www.ecolag.univmontp2.fr/index.php?option=com_content&task=view&id=219&Itemid=125).
4. The evenness of traits (FEve) is the distribution law of the mean value of species traits in the trait space (Eqs.4-6).
$FEve=\frac{\sum\limits_{b=1}^{S-1}{\min }\left( PE{{W}_{b}}\ \times \ \frac{1}{S-1} \right)-\frac{1}{S-1}}{1-\frac{1}{S-1}}$
$PE{{W}_{b}}=\frac{E{{W}_{b}}}{\sum\limits_{b=1}^{S-1}{E{{W}_{b}}}}$
$E{{W}_{b}}=\frac{{{d}_{ij}}}{{{w}_{i}}+{{w}_{j}}}$
where S is the number of species, PEWb is the local weighted average uniformity, EWb is the weighted average uniformity, wi is the relative abundance of species i, and dij is the Euclidean distance between species i and species j.
5. FDiv was calculated as follows:
${{g}_{k}}=\frac{1}{S}\underset{i=1}{\overset{S}{\mathop \sum }}\,{{x}_{ik}}$
$d{{G}_{i}}=\sqrt{\underset{k=1}{\overset{T}{\mathop \sum }}\,{{\left( {{x}_{ik}}-{{g}_{k}} \right)}^{2}}}$
$\overline{dG}=\frac{1}{S}\underset{i=1}{\overset{S}{\mathop \sum }}\,d{{G}_{i}}$
$\Delta d=\underset{i=1}{\overset{S}{\mathop \sum }}\,{{w}_{i}}\text{ }\!\!\times\!\!\text{ }\left( d{{G}_{i}}-\overline{dG} \right)$
$\Delta \left| d \right|=\underset{i=1}{\overset{S}{\mathop \sum }}\,{{w}_{i}}\text{ }\!\!\times\!\!\text{ }\left| d{{G}_{i}}-\overline{dG} \right|$
$FDiv=\frac{\Delta d+\overline{dG}}{\Delta \left| d \right|+\overline{dG}}$
where S is the number of species, xik is the value of species i trait k, gk is the center of trait k, T is the number of traits, dGi is the average distance between species i and the center of gravity, and d is the dispersion degree with multiple degrees as the weight, wi is the relative abundance of species i; $\bar{dG}$is the mean of dGi.
6. FDis is a new index from FDiv's further improvement.
$c=[c_{i}]=\frac{\sum a_{j}}{\sum a_{j} x_{ij}}$
$FDis=\frac{\mathop{\sum }^{}{{a}_{j}}{{z}_{j}}}{\mathop{\sum }^{}{{a}_{j}}}$
where aj is the relative abundance of species j and xij is the ith trait of species j. zj is the weighted distance from species j to center c. Where aj is the relative abundance of species j and xij is the ith trait of species j. zj is the weighted distance from species j to center c.
7. Rao defines a quadratic entropy equation that can be used to quantify the diversity and variability within a system.
$RaoQ=\underset{i=1}{\overset{S-1}{\mathop \sum }}\,\underset{j=i+1}{\overset{S}{\mathop \sum }}\,{{d}_{ij}}{{p}_{i}}{{p}_{j}}$
where S is the number of species, xij is the t-characteristic of species i, pi is the relative abundance of species i, and dij is the distance of species i and j in the trait space.
1
Biswas S R, Mallik A U, Braithwaite N T , et al. 2019. Effects of disturbance type and microhabitat on species and functional diversity relationship in stream-bank plant communities. Forest Ecology and Management, 432:812-822.

DOI

2
Botta-Dukát Z . 2005. Rao’s quadratic entropy as a measure of functional diversity based on multiple traits. Journal of Vegetation Science, 16(5):533-540.

DOI

3
Casanoves F, Pla L, di Rienzo J A , et al. 2011. FDiversity: A software package for the integrated analysis of functional diversity. Methods in Ecology and Evolution, 2(3):233-237.

DOI

4
Cornelissen J H C, Lavorel S, Garnier E , et al. 2003. A handbook of protocols for standardised and easy measurement of plant functional traits worldwide. Australian Journal of Botany, 51(4):335-380.

DOI

5
Cornwell W K, Schwilk D W, Ackerly D D . 2006. A trait-based test for habitat filtering: Convex hull volume. Ecology, 87(6):1465-1471.

PMID

6
Correia D L P, Raulier F, Bouchard M , et al. 2018. Response diversity, functional redundancy, and post-logging productivity in northern temperate and boreal forests. Ecological Applications, 28(5):1282-1291.

DOI

7
Craven D, Eisenhauer N, Pearse W D , et al. 2018. Multiple facets of biodiversity drive the diversity-stability relationship. Nature Ecology & Evolution, 2(10):1579-1587.

8
de Deyn G B, Cornelissen J H C, Bardgett R D . 2008. Plant functional traits and soil carbon sequestration in contrasting biomes. Ecology Letters, 11(5):516-531.

DOI

9
de la Riva E G, Violle C, Pérez-Ramos I M , et al. 2018. A multidimensional functional trait approach reveals the imprint of environmental stress in Mediterranean woody communities. Ecosystems, 21(2):248-262.

DOI

10
Faucon M P, Houben D, Lambers H . 2017. Plant functional traits: Soil and ecosystem services. Trends in Plant Science, 22(5):385-394.

DOI

11
Fontana S, Petchey O L, Pomati F . 2016. Individual-level trait diversity concepts and indices to comprehensively describe community change in multidimensional trait space. Functional Ecology, 30(5):808-818.

DOI

12
Garcia L C, Cianciaruso M V, Ribeiro D B , et al. 2015. Flower functional trait responses to restoration time. Applied Vegetation Science, 18(3):402-412.

DOI

13
Greenwood S, Ruiz-Benito P, Martinez-Vilalta J , et al. 2017. Tree mortality across biomes is promoted by drought intensity, lower wood density and higher specific leaf area. Ecology Letters, 20(4):539-553.

DOI PMID

14
He N P, Li Y, Liu C C , et al. 2020. Plant trait networks: Improved resolution of the dimensionality of adaptation. Trends in Ecology & Evolution, 35(10):908-918.

DOI

15
He N P, Liu C C, Piao S L , et al. 2019a. Ecosystem traits linking functional traits to macroecology. Trends in Ecology & Evolution, 34(3):200-210.

DOI

16
He N P, Liu C C, Tian M , et al. 2018. Variation in leaf anatomical traits from tropical to cold-temperate forests and linkage to ecosystem functions. Functional Ecology, 32(1):10-19.

DOI

17
He P C, Wright I J, Zhu S D , et al. 2019b. Leaf mechanical strength and photosynthetic capacity vary independently across 57 subtropical forest species with contrasting light requirements. New Phytologist, 223(2):607-618.

DOI

18
Kearney M, Porter W P . 2006. Ecologists have already started rebuilding community ecology from functional traits. Trends in Ecology & Evolution, 21(9):481-482.

DOI

19
Kearsley E, Hufkens K, Verbeeck H , et al. 2019. Large-sized rare tree species contribute disproportionately to functional diversity in resource acquisition in African tropical forest. Ecology and Evolution, 9(8):4349-4361.

DOI

20
Komac B, Pladevall C, Domènech M , et al. 2015. Functional diversity and grazing intensity in sub-alpine and alpine grasslands in Andorra. Applied Vegetation Science, 18(1):75-85.

DOI

21
Kordbacheh F, Jarchow M, English L , et al. 2019. Productivity and diversity of annually harvested reconstructed prairie communities. Journal of Applied Ecology, 56(2):330-342.

DOI PMID

22
Kunstler G, Falster D, Coomes D A , et al. 2016. Plant functional traits have globally consistent effects on competition. Nature, 529(7585):204-207.

DOI PMID

23
Laliberté E, Legendre P . 2010. A distance-based framework for measuring functional diversity from multiple traits. Ecology, 91(1):299-305.

PMID

24
Leps J, de Bello F, Lavorel S , et al. 2006. Quantifying and interpreting functional diversity of natural communities: Practical considerations matter. Preslia, 78:481-501.

25
Li L, McCormack M L, Ma C G , et al. 2015. Leaf economics and hydraulic traits are decoupled in five species-rich tropical-subtropical forests. Ecology Letters, 18(9):899-906.

DOI PMID

26
Li Y, He N P, Hou J H , et al. 2018. Factors influencing leaf chlorophyll content in natural forests at the biome scale. Frontiers in Ecology and Evolution, 6:64. DOI: 10.3389/fevo.2018.00064.

DOI

27
Liu C C, He N P, Zhang J H , et al. 2018. Variation of stomatal traits from cold temperate to tropical forests and association with water use efficiency. Functional Ecology, 32(1):20-28.

DOI

28
Liu C C, Li Y, Xu L , et al. 2019. Variation in leaf morphological, stomatal, and anatomical traits and their relationships in temperate and subtropical forests. Scientific Reports, 9:5803. DOI: 10.1038/s41598-019-42335-2.

DOI

29
Lozanovska I, Ferreira M T, Segurado P , et al. 2018. Limited resilience in hotspots of functional richness: The Mediterranean riparian shrublands. Aquatic Sciences, 80(3):1-12.

DOI

30
Ma Z Q, Guo D L, Xu X L , et al. 2018. Evolutionary history resolves global organization of root functional traits. Nature, 555(7694):94-97.

DOI

31
Mahaut L, Fort F, Violle C , et al. 2020. Multiple facets of diversity effects on plant productivity: Species richness, functional diversity, species identity and intraspecific competition. Functional Ecology, 34(1):287-298.

DOI PMID

32
Mason N W H, Mouillot D, Lee W G , et al. 2005. Functional richness, functional evenness and functional divergence: The primary components of functional diversity. Oikos, 111(1):112-118.

DOI

33
Mensah S, Veldtman R, Assogbadjo A E , et al. 2016. Tree species diversity promotes aboveground carbon storage through functional diversity and functional dominance. Ecology and Evolution, 6(20):7546-7557.

DOI

34
Nanes B A . 2015. Slide set: Reproducible image analysis and batch processing with Image[J]. BioTechniques, 59:269-278.

35
Pasari J R, Levi T, Zavaleta E S , et al. 2013. Several scales of biodiversity affect ecosystem multifunctionality. Proceedings of the National Academy of Sciences of the USA, 110(25):10219-10222.

DOI

36
Pérez-Harguindeguy N, Díaz S, Garnier E , et al. 2013. New handbook for standardised measurement of plant functional traits worldwide. Australian Journal of Botany, 61(3):167-234.

DOI

37
Petchey O L, Gaston K J . 2002. Functional diversity (FD), species richness and community composition. Ecology Letters, 5(3):402-411.

DOI

38
Petchey O L, Gaston K J . 2006. Functional diversity: Back to basics and looking forward. Ecology Letters, 9(6):741-758.

PMID

39
Petchey O L, Hector A, Gaston K J . 2004. How do different measures of functional diversity perform? Ecology, 85(3):847-857.

DOI

40
Pillar V D, Blanco C C, Müller S C , et al. 2013. Functional redundancy and stability in plant communities. Journal of Vegetation Science, 24(5):963-974.

DOI

41
Poos M S, Walker S C, Jackson D A . 2009. Functional-diversity indices can be driven by methodological choices and species richness. Ecology, 90(2):341-347.

DOI

42
Reich P B . 2014. The world-wide ‘fast-slow’ plant economics spectrum: A traits manifesto. Journal of Ecology, 102(2):275-301.

DOI

43
Ross S R P J, Hassall C, Hoppitt W J E , et al. 2017. Incorporating intraspecific trait variation into functional diversity: Impacts of selective logging on birds in Borneo. Methods in Ecology and Evolution, 8(11):1499-1505.

DOI

44
Schleuter D, Daufresne M, Massol F , et al. 2010. A user’s guide to functional diversity indices. Ecological Monographs, 80(3):469-484.

DOI

45
Schmera D, Erős T, Podani J . 2009. A measure for assessing functional diversity in ecological communities. Aquatic Ecology, 43(1):157-167.

DOI

46
Solefack M C M, Fedoung E F, Temgoua L F . 2018. Factors determining floristic composition and functional diversity of plant communities of Mount Oku forests, Cameroon. Journal of Asia-Pacific Biodiversity, 11(2):284-293.

DOI

47
Sonkoly J, Kelemen A, Valkó O , et al. 2019. Both mass ratio effects and community diversity drive biomass production in a grassland experiment. Scientific Reports, 9:1848. DOI: 10.1038/s41598-018-37190-6.

DOI PMID

48
Staples T L, Dwyer J M, England J R , et al. 2019. Productivity does not correlate with species and functional diversity in Australian reforestation plantings across a wide climate gradient. Global Ecology and Biogeography, 28(10):1417-1429.

49
Villéger S, Mason N W H, Mouillot D . 2008. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecology, 89(8):2290-2301.

DOI

50
Walker B, Kinzig A, Langridge J . 2005. Plant attribute diversity, resilience, and ecosystem function: The nature and significance of dominant and minor species. Ecosystems, 2(2):95-113.

DOI

51
Wang R L, Yu G R, He N P , et al. 2016. Latitudinal variation of leaf morphological traits from species to communities along a forest transect in Eastern China. Journal of Geographical Sciences, 26(1):15-26.

DOI

52
Wright I J, Reich P B, Westoby M , et al. 2004. The worldwide leaf economics spectrum. Nature, 428(6985):821-827.

PMID

53
Wright J P, Naeem S, Hector A , et al. 2006. Conventional functional classification schemes underestimate the relationship with ecosystem functioning. Ecology Letters, 9(2):111-120.

DOI

54
Xu L, Wang C Y, Zhu J X , et al. 2018a. Latitudinal patterns and influencing factors of soil humic carbon fractions from tropical to temperate forests. Journal of Geographical Sciences, 28(1):15-30.

DOI

55
Xu Z W, Li M H, Zimmermann N E , et al. 2018b. Plant functional diversity modulates global environmental change effects on grassland productivity. Journal of Ecology, 106(5):1941-1951.

DOI

56
Zhang X S, Yang D A . 1995. Allocation and study on global change transects in China. Quaternary Sciences, 15(1):43-52.

57
Zhao N, He N P, Wang Q F , et al. 2014. The altitudinal patterns of leaf C : N : P stoichiometry are regulated by plant growth form, climate and soil on Changbai Mountain, China. Plos One, 9(4):e95196. DOI: 10.1371/journal.pone.0095196.

DOI

Outlines

/