Skip to main content

Plant economics spectrum governs leaf nitrogen and phosphorus resorption in subtropical transitional forests

Abstract

Background

Leaf nitrogen (N) and phosphorus (P) resorption is a fundamental adaptation strategy for plant nutrient conservation. However, the relative roles that environmental factors and plant functional traits play in regulating N and P resorption remain largely unclear, and little is known about the underlying mechanism of plant functional traits affecting nutrient resorption. Here, we measured leaf N and P resorption and 13 plant functional traits of leaf, petiole, and twig for 101 representative broad-leaved tree species in our target subtropical transitional forests. We integrated these multiple functional traits into the plant economics spectrum (PES). We further explored whether and how elevation-related environmental factors and these functional traits collectively control leaf N and P resorption.

Results

We found that deciduous and evergreen trees exhibited highly diversified PES strategies, tending to be acquisitive and conservative, respectively. The effects of PES, rather than of environmental factors, dominated leaf N and P resorption patterns along the elevational gradient. Specifically, the photosynthesis and nutrient recourse utilization axis positively affected N and P resorption for both deciduous and evergreen trees, whereas the structural and functional investment axis positively affected leaf N and P resorption for evergreen species only. Specific leaf area and green leaf nutrient concentrations were the most influential traits driving leaf N and P resorption.

Conclusions

Our study simultaneously elucidated the relative contributions of environmental factors and plant functional traits to leaf N and P resorption by including more representative tree species than previous studies, expanding our understanding beyond the relatively well-studied tropical and temperate forests. We highlight that prioritizing the fundamental role of traits related to leaf resource capture and defense contributes to the monitoring and modeling of leaf nutrient resorption. Therefore, we need to integrate PES effects on leaf nutrient resorption into the current nutrient cycling model framework to better advance our general understanding of the consequences of shifting tree species composition for nutrient cycles across diverse forests.

Peer Review reports

Introduction

Nutrient resorption is the process by which plant withdraws nutrients from senescing tissues before litterfall and thus is recognized as an adaptation strategy for plant nutrient conservation [1,2,3,4]. This process prolongs nutrient mean residence time, which enables plants to depend less on current nutrient uptake capacity to mitigate the nutrient limitation of plant production, thereby enhancing overall nutrient use efficiency and improving plant fitness in natural ecosystems [5,6,7]. Nutrient resorption represents a fundamental pathway of nutrient recycling in forest ecosystems, especially for nitrogen (N) and phosphorus (P), which are heavily involved in plant physiological processes and are integral components in many metabolic compounds, such as proteins and Ribonucleic Acid (RNA) [8, 9]. Plant N and P resorption is estimated to provide an average of 31% and 40% of the annual plant N and P requirements for forests globally, respectively [10]. Nutrient resorption efficiency (RE) and nutrient resorption proficiency (RP) have been well-recognized as complementary metrics of leaf nutrient resorption [3, 5, 11]. For leaves, RE indicates the proportion of nutrients that were resorbed before the abscission of green leaves, and RP represents the lowest nutrient concentration in the senescent leaves, which is considered to be the biochemical limit of resorption [5]. Therefore, conducting a quantitative analysis of leaf N and P resorption is necessary for gaining a better understanding of nutrient utilization and adaptation strategies of tree species and of how to simulate the nutrient cycling process of forest ecosystems.

Decades of research have illustrated environmental conditions (including climate and edaphic properties) and plant functional traits as predominant controls over leaf N and P resorption [2, 4, 12]. However, their relative roles in the leaf N and P resorption process are still elusive. Previous studies have revealed that N resorption efficiency (NRE) decreased as mean annual temperature (MAT) and mean annual precipitation (MAP) increased, while P resorption efficiency (PRE) followed opposite trends [12]. However, other studies have suggested that both NRE and PRE decreased as MAT and MAP increased [2, 3]. Because our current understanding of such associations between leaf N and P resorption and climate is primarily derived from data-integration studies at the global and regional scales, these inconsistent conclusions may be attributed to the diverse species compositions in the respective study areas.

Soil nutrient availability is another major environmental driver of leaf N and P resorption [2, 6]. It is recognized as a key factor that mediates the impact of climate on leaf N and P resorption [7, 13]. However, there is still an ongoing debate as to whether and how N and P resorption relates directly to soil nutrient availability. One possible reason is that the relative costs of resorbing N and P back into live tissues from senescing ones, instead of uptaking new N and P from the soil, remain largely uncertain. Particularly, nutrient constraints on plant growth and productivity transition from P limitation at low latitudes to N limitation at high latitudes [4, 14, 15], leading plants to trade off different nutrient acquisition strategies depending on soil nutrient availability in tropical versus temperate regions. Furthermore, our nuanced understanding of nutrient availability and utilization in subtropical transitional forests remains unclear, seriously hindering a more general understanding of N and P resorption.

In addition to the above-mentioned environmental controls, plant functional traits are also the intrinsic biotic factors affecting N and P resorption. Current studies on leaf N and P resorption of different tree species mainly focus on leaf habits (i.e., deciduous and evergreen). Previous studies have demonstrated that evergreen tree species in tropical forests had higher NRE and PRE than deciduous counterparts [2, 16], but the NRE of evergreen trees was seemingly lower, and the PRE showed no significant difference in temperate forests [17, 18]. However, patterns and controls of nutrient resorption in subtropical transitional forests between tropical and temperate forests remain relatively unclear. These unexplained mechanisms call for the urgent need to link functional traits closely related to leaf habit to leaf N and P resorption, such as specific leaf area (SLA), leaf dry matter content (LDMC), and green leaf N and P concentrations [19, 20].

The analysis of individual functional traits can contribute to identifying specific traits that influence N and P resorption. In some cases, however, a certain plant trait may have reached its physiological threshold, rendering it insufficient for mitigating increasing environmental stress [21,22,23]. Under such circumstances, these different combinations of intercorrelating physiological and morphological traits led to the construction of trait coordination networks [24, 25]. Quantifying trait network patterns across leaf habits is essential for understanding the diversity of plant form and function, as well as for identifying the main drivers of ecological processes such as leaf nutrient resorption [25, 26]. The key traits related to leaf N and P economics could reflect trait-based adaptive strategies [27]. The differences between trait networks of deciduous and evergreen trees can characterize their respective manners of acquiring N and P resources [26, 28].

In recent decades, the plant economics spectrum (PES) has been proposed as a method to reflect the pervasive trade-offs and coordination in plant nutrient utilization by using the components and coordination relationships of trait networks to integrate a unified concept across plant taxonomy and growth forms [19, 29]. By integrating multivariate functional traits into PES, ecological trade-offs may become apparent, thus facilitating predictions of some important ecosystem processes such as nutrient resorption [21, 30]. Nonetheless, current studies examining the association between N and P resorption and PES have produced inconsistent patterns, or even uncoupled relationships [6, 27, 31]. This phenomenon may be mainly attributed to particular adaptation strategies for each PES trait associated with N and P economics under environmental selection, forming multiple counteracting pathways affecting resorption processes [31]. Species with similar or opposite strategies converge or diverge along the axis of trait variation according to PES. Deciduous trees tend to adopt an acquisitive strategy, whereas evergreen trees tend to exhibit a conservative strategy [19, 28]. Furthermore, there have been significant advances in our understanding of PES for plant functional traits related to N and P economics at different spatial and biological scales [20, 31], but quantitative studies on how PES drives N and P resorption in leaves are still lacking [20, 27]. Therefore, we expected that PES well describe a trait syndrome of leaf habits that is closely associated with N and P resorption.

To date, although the individual effects of environmental factors and plant functional traits on leaf N and P resorption have been well documented, their relative roles at different scales remain largely unclear [2, 16, 32]. Global-scale studies have revealed that MAT and MAP significantly affect leaf N and P resorption, regional-scale studies have demonstrated that soil nutrient availability and leaf life span are dominant factors controlling leaf N and P resorption, but local-scale studies have indicated that leaf N and P resorption are tightly linked with tree species composition and plant functional traits [6, 12, 16, 32]. More critically, different spatial scales are usually accompanied by multiple differences in environmental conditions and clear distinct tree species compositions, which complicates the issue of dominant factors for nutrient resorption, and further hinders our robust understanding of regulatory mechanisms of plant nutrient resorption strategies in these diverse environments [23, 32, 33]. These aforementioned substantial complexity and uncertainty call for simultaneous studies on plant functional traits under continuous, wide-ranging environmental conditions. Elevational gradient precisely provides a continuously changing environment and bridges the extensive gaps in tree species composition and environmental isolation with its gradual change, is thus an ideal platform for exploring how plant functional traits and environmental controls over nutrient resorption simultaneously [34, 35].

In this study, we have investigated and quantified relationships between plant functional traits and leaf N and P resorption of broad-leaved tree species along an elevational gradient to get a more nuanced understanding of the drivers of nutrient resorption in subtropical China. To this end, we attempted to test the following hypotheses: (a) leaf habits directly explain more variations in leaf N and P resorption of broad-leaved tree species than elevation, meaning that N and P resorption should be more strongly controlled by morphological and physiological trait combinations than by climate and soil conditions associated with elevational gradients; (b) deciduous and evergreen trees differentiate along an integrated PES, arraying on the acquisitive side and the conservative side, respectively; (c) PES is strongly coupled with N and P resorption, and these mediating associations can be largely attributed to one or more crucial traits.

Materials and methods

Study area

This study was carried out on the southern slope of Shennongjia Mountain (31°19′4″ N, 110°29′44″ E), northwest of Hubei Province, in central China. The region belongs to the transition zone between the (sub-) tropical and temperate climates and is an important biodiversity hotspot in China and globally. The mean annual temperature of this area is 10.6 °C, and annual precipitation ranges from 1306 to 1722 mm. The dominant soil classes are mountain yellow–brown soil and mountain brown soil, with a typical subtropical forest vegetation of deciduous and evergreen broad-leaved mixed forest [36].

Experimental design and field sampling

We conducted a field forest inventory on the southern slope of Shennongjia Mountain in July 2019 to identify common broad-leaved tree species that occurred with high frequency. We set up six sampling sites at intervals of 200 m from 800 to 1800 m above sea level (Fig. 1). We sequentially sampled all of the common representative broad-leaved tree species identified through pre-inventory on six sites (± 20 m elevation) and previous studies [37, 38]. For each site, we followed similar previous studies to sample [26, 31]. Specifically, we treated an individual mature tree as a sampling unit, with three repeat units for each tree species, at least 15 m apart between each two sampling units. Overall, we sampled 101 broad-leaved species, belonging to 33 families and 60 genera and formed 205 elevation-species combinations. The identification of tree species was based on taxonomic expertise, and the nomenclature followed the Flora of China [39], and all selected tree species were classified as evergreen (n = 43) or deciduous (n = 58) according to their leaf habit. Detailed sampling information is shown in Table S1.

Fig. 1
figure 1

Sampling diagram of broad-leaved tree species on the southern slope of Shennongjia Mountain. Asterisks represent the six sampling transects along the elevational gradient. The size of the pie chart represents the total number of tree species as shown in the left-most concentric circle, and the labels in the pie charts represent the numbers of deciduous and evergreen tree species sampled at each elevation. The background of the mountain was designed by pch.vector, Freepik (http://www.freepik.com/). Tree illustrations are from Tracey Saxby, Dylan Taillie, Kim Kraeer, and Lucy Van Essen-Fishman, Integration and Application Network, University of Maryland Center for Environmental Science (http://ian.umces.edu/imagelibrary/)

We took the tree height as the selection criteria to identify mature and healthy individual trees as sampling units, referring to the Flora of China [39]. All sampling units are labeled to ensure the source of the mature leaf samples corresponds to that of the senescent leaf samples. In July and August of 2019 and 2020, we collected 30–50 undamaged, relatively large, and fully expanded fresh mature leaf samples from the different canopy positions of each sampling unit, and we cut off representative twigs from each sample. We wore gloves during sampling to minimize sample contamination. Based on the previous field investigation experience in our target subtropical forests and literature records [40, 41], the peak period of leaf abscission for deciduous tree species is October and November in autumn, and that of evergreen tree species is October and November in autumn or April and May in spring. Therefore, in October and November of 2019 and 2020, we collected senescent leaves of all deciduous trees and partial evergreen trees, and the remaining senescent leaves of evergreen trees were collected in April and May 2020 to ensure that all senescent leaf samples were from the peak period of leaf abscission. We collected yellowed senescent leaves by shaking trees or picking up newly fallen leaves on the ground to avoid nutrient leaching from partly decomposed leaf litter on the forest floor, which could lead to an underestimation of leaf N and P concentrations [22, 42].

Measurement of plant functional traits

Here, we selected and determined 13 fundamental traits of fresh leaves according to previous standard protocols [43,44,45] (Table S2). Mature fresh leaf samples were stored in sealed plastic bags immediately after collection in the field, and they were kept cool until being brought to the laboratory. On the day of sampling, we selected 3 intact fresh leaves from each sampling unit, separating the petioles for the following physical trait measurements [42, 46]. We then dried the remaining leaves in the oven at 65 °C for 48 h for subsequent chemical analysis [36, 46]. For each leaf subsample, we measured the relative content of Chlorophyll (SPAD) using a hand-held portable chlorophyll meter (SPAD-502), which was repeated at least three times for each leaf to avoid the possible impacts of leaf thickness and hair coats on SPAD. Leaf thickness (LT, mm) was quantified with a vernier caliper, and care was taken to avoid the main veins. Leaf samples were scanned by a leaf area meter (CI-203, USA) to determine leaf area (LA, cm2). The leaf subsamples were then weighed to determine the leaf fresh mass, and then they were dried in the oven at 65 °C for 48 h to constant mass to determine the leaf dry mass [36, 46]. This data was then used to calculate specific leaf area (SLA, leaf area/leaf dry mass, cm2 g−1) and leaf dry matter content (LDMC, leaf dry mass/leaf fresh mass). The length (PL, mm) and width (PW, mm) of each leaf sample petiole were measured by a vernier caliper. The twig dry matter content (TDMC) of each tree was calculated using the same method as the leaf. Then twig wood density (TWD, g cm−3) was calculated using the drainage method [47].

On the day of sampling, fresh and senescent leaves were dried in the oven at 65 °C for 48 h to constant mass [36, 46], and then they were ground with a laboratory mill for analysis of chemical elements. The leaf carbon (C) and N concentrations per unit mass for fresh green leaves (Cgr, Ngr, mg g−1) and senescent leaves (Csen, Nsen, mg g−1) were determined by an element analyzer (vario EL cube CHNOS Elemental Analyzer, Elementar Analysensysteme GmbH, Hanau, Germany). We also calculated the ratio of leaf C to N concentrations (CNgr) to serve as an indicator trait. The leaf P concentration per unit mass for fresh leaves (Pgr, mg g−1) and senescent leaves (Psen, mg g−1) were determined using inductively-coupled plasma spectrometry after digestion of samples in HNO3 (iCAP 6300 ICP-OES Spectrometer, Thermo Fisher, USA). More details are provided in Table S2.

Measurement of environmental factors

We measured a total of 11 environmental factors in this study (Table S3). At each elevation site, we used HOBO Onset microclimatic recorders (Onset Computer Corporation, USA) to measure actual microclimatic variables. These included air temperature (AT, °C), soil temperature (ST, °C), and soil moisture (SM, %). At each level, we selected five 2 m × 2 m plots with high frequencies of sampled trees and collected soil samples at a depth of 0–10 cm [18]. Visible plant materials and stones were removed from all soil samples. Next, the soil samples were air-dried at room temperature and either passed through a 100-mesh sieve for subsequent chemical analysis or frozen at − 20 °C prior to microbiological assay. Soil total N concentration (SN, mg g−1) was measured using the Kjeldahl method after digesting samples with H2SO4, and available N concentration (SNA, mg g−1) was determined using the diffusion method after alkaline hydrolysis. Soil total P concentration (SP, mg g−1) was quantified after wet digestion with concentrated HF and HClO4, and soil available P concentration (SPA, mg g−1) was determined in NaHCO3 extraction, using the molybdenum antimony resistance colorimetry method. Soil organic C concentration (SOC, mg g−1) was measured by titrimetry after oxidation with a mixture of potassium dichromate and sulfuric acid. Soil pH was measured using a pH electrode (PB-10, Sartorius, Germany) in a soil–water suspension (soil: water = 1:2.5 [v/w]). Soil microbial biomass C and N concentration (MBC, MBN, mg g−1) were quantified using chloroform fumigation (multi N/C 3100, TOC/TNb analyzer, Analytik-Jena AG, Germany). The specific and detailed experimental procedures can be found in Protocols for Standard Biological Observation and Measurement in Terrestrial Ecosystems [46]. More data details are provided in Table S3.

Nutrient resorption calculation

We used two fundamental and complementary metrics to quantify nutrient resorption: resorption efficiency (RE, %) and resorption proficiency (RP, mg g−1) [2, 6]. RE was calculated as the ratio of the difference in nutrient concentrations between green and senescent leaves to green leaf nutrient concentrations. Considering the mass loss that occurs during leaf senescence, we used a mass loss correction factor (MLCF, 0.78 for evergreen broad-leaved tree species and 0.784 for deciduous broad-leaved tree species) to correct nutrient concentrations in senescent leaves [3]. RP was measured as the nutrient concentration in senescent leaves, with low litter nutrient concentrations corresponding to high RP and vice versa [5]. At each elevational site, we calculated the average value of each trait and nutrient resorption for each species.

Data analysis

We checked the data for normality and log-transformed it to improve the linearity of the relationships if necessary before conducting the analysis. To explore the effects of leaf habit and elevation on RE and RP, we used a two-way analysis of variance (ANOVA) followed by t-tests and Tukey’s post-hoc test. To explore the relative importance of environmental factors and plant functional traits on leaf N and P resorption, we have added variance partition analysis. Based on the Pearson correlation, we further examined the intercorrelation of environmental factors and the relationships between individual environmental factors and nutrient resorption using Mantel’s test in the VEGAN package [48]. In the VEGAN package we also then performed redundancy analysis (RDA) to test overall associations between environmental factors and nutrient resorption [48].

Since tree species phylogeny may influence the process of nutrient resorption [5], we used the V.PHYLOMAKER package [49], based on the published phylogenetic trees as background trees, to extract the common target nodes through tree species names and generated our target phylogenetic trees. The PLANTLIST package [50] was used to match the same canonical family, genus, and species information as the V.PHYLOMAKER package. We used Blomberg’s K value and Pagel’s λ value to evaluate the phylogenetic signal on each plant trait using the PICANTE and PHYTOOLS packages [51, 52]. A large K value and λ value with p < 0.05 indicated conservatism for a trait. We described plant functional trait coordination networks according to statistically significant Pearson correlations among traits and illustrated them using the IGRAPH package [53]. We adopted two indicators to characterize the centrality of networks: D (degree, the number of adjacent edges for each vertex) and D W (the sum of the edge weights of the adjacent edges for each vertex) [54].

To integrate these interrelated traits into PES, we performed a principal component analysis (PCA) using the basic function ‘prcomp’. Corresponding illustrations were made in the GGBIPLOT package [55]. To evaluate the difference in functional traits across leaf habits, we compared the trait scores of the first two axes (PC1 and PC2) between deciduous and evergreen trees using t-tests. Furthermore, to assess and compare the directions and strengths of principal components on N and P resorption across leaf habits, we conducted linear mixed effect models (LMMs) with tree species as random factors and the PC1 and PC2 scores as fixed effects using the LME4 package [56]. LMMs were used to deal with elevation-species pseudo-replication. To better identify and quantify which traits or trait combinations varied in relation to N and P resorption, we re-applied the similar LMMs with tree species as random factors and the 13 functional traits mentioned above as fixed effects to evaluate the role of each single trait. We removed the variables with the variance inflation factor (VIF) value > 10 to avoid the potential multicollinearity of the models. These were visualized using the GGPLOT2 package [57]. All analyses and illustrations were performed in R 4.0.3 [58].

Results

Patterns of leaf N and P resorption across elevations and leaf habits

The overall NRE and PRE of broad-leaved trees were 39.65% and 50.76%, respectively. NRE was significantly higher in deciduous trees than in evergreen trees (p < 0.05, Fig. 2), but PRE displayed no significant difference across leaf habits (p > 0.05, Fig. 2). Mean leaf litter N and P concentrations of broad-leaved trees were 13.08 mg g−1 and 0.74 mg g−1, respectively, generally higher than the threshold of incomplete resorption proposed by Killingbeck [5], indicating low N and P resorption proficiency (NRP and PRP). Moreover, deciduous trees exhibited higher N and P concentrations in senescent leaves than did evergreen trees, which suggests that deciduous trees have lower resorption proficiency (p < 0.05, Fig. 2). By contrast, NRE, NRP, PRE, and PRP were all insignificantly associated with elevation (p > 0.05, Table 1, Fig. S1). Furthermore, when considering the interaction of elevation and leaf habit, we found no significant relationship between NRE, NRP, and the interaction, but PRE and PRP were significantly affected by the interaction (Table 1).

Fig. 2
figure 2

N and P resorption efficiency (RE) and proficiency (RP) across leaf habits. Asterisks indicate significant differences between evergreen and deciduous tree species (**p < 0.01, ***p < 0.001), and NS indicates no significant difference. The dotted lines represent the thresholds for incomplete (N > 10 mg g−1, P > 0.5 mg g−1 for evergreen tree species, P > 0.8 mg g−1 for deciduous tree species) and complete (N < 7 mg g−1, P < 0.4 mg g−1 for evergreen tree species, P < 0.5 mg g.−1 for deciduous tree species) resorption, as defined by Killingbeck (1996)

Table 1 Results of two-way ANOVA for elevation and leaf habit effects on N and P resorption efficiency and proficiency

Environmental effects on leaf N and P resorption

We found that environmental factors only explained 0.9% of the variation of leaf N and P resorption, while plant functional traits explained 13.8% of that, indicating that the relative effects of environmental factors are much weaker than plant functional traits (Fig. S2). In addition, the RDA results revealed that environmental factors played a minor role in N and P resorption along the elevational gradient. Environmental factors merely contributed 4.19% (R 2 = 4.19%, Fig. S3) of variation in N and P resorption. Furthermore, Mantel’s test suggested that individual environmental factors do not significantly control N and P resorption (p > 0.05, r < 0.2, Fig. 3). Given the minor effects of the environmental factors, we mainly focused on the roles of plant functional traits in N and P resorption below analysis.

Fig. 3
figure 3

Environmental factors play a minor role in nutrient resorption. Correlation between environmental factors and their effects on N and P resorption. Asterisks indicate the significance of effects among environmental factors (*p < 0.05, **p < 0.01, ***p < 0.001). The edge color (r-size) corresponds to the coefficients between resorption characteristics and environmental variables. Edge width denotes the statistical significance (p-value) between each trait and environmental variable. AT, Air temperature; ST, Soil temperature; SM, Soil moisture; MBC, Soil microbial biomass carbon content; MBN, Soil microbial biomass nitrogen content; pH, Soil pH; SOC, Soil organic carbon content; SNA, Soil available nitrogen content; SP, Soil total phosphorus content; SPA, Soil available phosphorus content

Associations between traits and leaf N and P resorption

We first constructed a phylogeny tree of our targeted tree species to assess phylogenetic affiliation (Fig. S4) and found that plant functional traits did not exhibit relatively high degrees of phylogenetic signal, with Blomberg’s K value ranging from 0.142 (SPAD) to 0.374 (TDMC) and Pagel’s λ value ranging from 0.082 (LA) to 0.787 (Cgr) (Table S4). This implied preference towards random trait evolution, with evolutionary patterns of traits likely to be influenced by non-phylogeny. Most of the plant functional traits correlated significantly with one another (p < 0.05, Fig. 4). For all trait networks, SLA exhibited the highest connectedness, followed by PL, TDMC, and green leaf nutrient concentrations. The trait networks of deciduous trees displayed higher transitivity than those of evergreen trees (Table S5).

Fig. 4
figure 4

Trait correlation networks of broad-leaved trees for (a) all, (b) deciduous, and (c) evergreen tree species. Traits in the network were represented as vertices and their correlations as the edges linking them. Different node colors represent different types of plant organs or tissues. Node size denotes the degree of connectedness. Red and blue edges represent positive and negative correlations (p < 0.05), respectively. SPAD, Chlorophyll relative content; LDMC, Leaf dry matter content; LA, Leaf area; SLA, Specific leaf area; LT, Leaf thickness; PL, Petiole length; PW, Petiole width; TDMC, Twig dry matter content; TWD, Twig wood density; Cgr, Carbon concentration in green leaves; Ngr, Nitrogen concentration in green leaves; Pgr, Phosphorus concentration in green leaves; CNgr, The ratio of carbon to nitrogen in green leaves

Given the close associations among traits, we performed PCA and found that PES could predict variation in N and P resorption well (Fig. 5 and Fig. 6). PC1 and PC2 explained 34.2% and 17.2% of trait variation, respectively. The 10 plant functional traits besides PW, TDMC, and Cgr all contributed significantly to PC1 (p < 0.05, Table S6). Along the PC1 axis, traits representing a conservative resource utilization strategy (SPAD, LDMC, LT, TWD, and CNgr) occupied the right end, while traits representing an acquisitive strategy (LA, SLA, PL, Ngr, and Pgr) were on the left. PC2 correlated significantly and positively with LA, LT, and PW, which represented a conservative strategy of tree structure and defense, and negatively with LDMC, TDMC, and TWD, which related to an acquisitive strategy of photosynthesis and nutrient supply. Most evergreen trees were distributed on the conservative side, while most deciduous species clustered on the acquisitive side (p < 0.05, Fig. 5).

Fig. 5
figure 5

Principal component analysis (PCA) of 13 plant functional traits. The confidence ellipses are drawn at the 90% level. The boxplots at the upper and right margins illustrate the differences in the distributions of tree species across leaf habits along the first and second PCA axes. Asterisks indicate significant differences between evergreen and deciduous trees (***p < 0.001). Abbreviations for plant functional traits are provided in Fig. 4

Fig. 6
figure 6

The effect size of the first (PC1: photosynthesis and nutrient resource utilization) and second (PC2: structural and functional investment) principal components of 13 plant functional traits of evergreen and deciduous trees on leaf N and P (a) resorption efficiency (RE) and (b) resorption proficiency (RP) in linear mixed effect models. The horizontal error bars show the 95% confidence intervals of the fixed effect size from linear mixed effect models. Asterisks indicate the significance of the effects of PC1 and PC2 on RE and RP (*p < 0.05, **p < 0.01, ***p < 0.001). Yellow and green represent deciduous and evergreen trees, respectively

The PC1 axis exerted significant negative effects on NRE, NRP, PRE, and PRP for both deciduous and evergreen trees (p < 0.05, Fig. 6), except that NRE of evergreen trees was not significantly affected. The PC2 axis showed no obvious impact on deciduous trees, while it positively affected NRE and negatively affected NRP and PRP of evergreen species (p < 0.05, Fig. 6). For both deciduous and evergreen trees, SLA, N, and P concentrations in green leaves explained N and P resorption the most (Fig. 7).

Fig. 7
figure 7

Effects of plant functional traits of deciduous (D) and evergreen (E) trees on leaf N and P (a) resorption efficiency (RE) and (b) resorption proficiency (RP) in linear mixed effect models. The color and size shown indicate the direction and strength of the coefficients from linear mixed effect models. Asterisks indicate the significance of the effects of PC1 and PC2 on RE and RP (·p < 0.01, *p < 0.05, **p < 0.01, ***p < 0.001). Abbreviations for plant functional traits are provided in Fig. 4

Discussion

Here, we have mainly focused on leaf N and P resorption for more representative species along an elevational gradient than in past studies [22, 23, 35], considering environmental factors and plant functional traits simultaneously to explore their roles. Our field observations confirmed previous modeling predictions of N and P limitation that our target subtropical forest is potentially less nutrient-limited [14]. We found large differences in leaf N and P resorption as a result of various functional trait combinations derived from leaf habit of tree species, rather than diverse environmental conditions along the elevational gradient. The highly diversified strategies of nutrient use and conservation in deciduous and evergreen trees (i.e., acquisitive and conservative strategies, respectively) led to pronounced divergence in leaf N and P resorption. These findings clearly highlight the urgent need to recognize that distinct controls of nutrient cycling dominate at different scales and suggest that biogeochemical cycle models should incorporate multiple traits and their coordination rather than just including individual traits.

Deciduous tree species tend to resorb N and P more effectively but less proficiently than evergreen counterparts

We found that deciduous trees resorbed N more effectively than evergreen counterparts, but there was only a slight difference in resorbing P of trees (Fig. 2), which was similar to findings in previous studies from Aert [11] and Tang et al. [59], but contradictory to global-scale findings by Yuan and Chen [12] and Vergutz et al. [3]. Different from insignificant differences in NRP and PRP between deciduous and evergreen trees [59, 60], our study demonstrated that deciduous trees were less proficient in resorbing N and P than evergreen trees (Fig. 2), concurred with Aerts [11]. On one hand, this may be derived from the observation that mature leaf N and P exhibited much higher concentrations in deciduous trees (Ngr = 18.67 mg g−1, Pgr = 1.39 mg g−1) compared to evergreen trees (Ngr = 14.16 mg g−1, Pgr = 1.01 mg g−1) in our study. Consequently, leaf litter N and P concentrations, as the final products reflecting the nutrient status of mature fresh leaves after resorption, are likely to display significant differences. Moreover, the leaves of deciduous tree species are typically composed of more soluble compounds, such as proteins and carbohydrates [61, 62]. These compounds break down more easily thus achieving more efficient nutrient resorption for deciduous tree species. In contrast, the leaves of evergreen tree species typically contain more insoluble structural compounds such as cellulose and lignin, hindering leaf nutrient resorption for evergreen species [61, 63].

On the other hand, this could be attributed to the fact that leaf N and P resorption between deciduous and evergreen trees varied with leaf longevity and duration of leaf abscission [2, 6, 64]. Deciduous trees, characterized by frequent leaf abscission and annual new leaf regeneration, tend to allocate more nutrients to facilitate faster growth and exhibit higher leaf production within their comparatively shorter growing seasons [65, 66]. Additionally, owing to their shorter abscission period and lower temperatures during the abscission season, the reaction time and activity of the decomposing enzymes for leaves may be insufficient, resulting in diminished soil nutrient availability and necessitating enhanced nutrient resorption to sustain growth and productivity [22]. In contrast, evergreen trees grow slowly and have long leaf lifespans, they develop adaptation strategies with low nutrient requirements to tolerate disadvantageous climatic and nutritional conditions, by maximizing the growth period and minimizing nutrient loss [6, 17]. Overall, due to the variable functional traits and nutrient strategies, N and P resorption varied across leaf habits [28, 67].

More notably, the strict classification of tree species into deciduous and evergreen could obscure important ecological variations in traits within each leaf habit, and the prominent distinction across leaf habits is essentially the result of a combination of many plant functional traits [64, 67, 68]. Therefore, considering the effect of leaf habit, a classification indicator only may be not enough, further analyses on the relationship between plant functional traits and N and P resorption are necessary.

Dominant effects of PES on leaf N and P resorption

We found no clear phylogenic signal of plant functional traits, which aligns with the point that the role of phylogeny was negligible when specific traits were considered [31, 69]. Other studies have also revealed that N and P resorption are not affected by phylogeny [70]. Therefore, we posit that the associations between plant functional traits and nutrient resorption in this study were phylogenetically independent. The strong associations of plant functional traits suggest that the evolution of traits is intimately linked [19, 29]. Trade-off strategies of resource acquisition and investment are defined as the trait combination form, that is, PES. One major dimension of our PES (PC1) represented photosynthesis and nutrient resource utilization, with a gradient of acquisitive strategy to conservative strategy running from left to right. The other dimension (PC2) reflected the structural and functional investment, with a gradient of investment preferences for photosynthesis and nutrient acquisition and for tree growth and structural conservation from bottom to top (Fig. 5). Evergreen and deciduous trees here were mainly distributed in the upper right and lower left parts of the PCA, respectively, which suggested that each leaf habit had obvious preferences for plant functional trait combinations along with nutrient strategies [67]. The findings support our hypothesis that deciduous and evergreen trees differentiate along an integrated PES, arraying on the acquisitive side and the conservative side, respectively.

Plant functional trait networks and PES approaches could be particularly useful in considering complex and dynamic processes such as nutrient resorption [20, 31, 45]. We did find significant correlations between PES and resorption for both N and P in our subtropical forest, thus supporting our hypothesis that trait coordination taken along with multi-trait variation on the PCA axes collectively determine trade-offs in nutrient resource utilization patterns. Remarkably, photosynthetic and nutrient resource utilization capacity positively affected N and P resorption for both deciduous and evergreen trees, whereas structure and defense investments only positively affected N and P resorption for evergreen trees (p < 0.05, Fig. 6). The differences in strategies between deciduous and evergreen trees is understandable, as they differ in phenology, leaf lifespan, and nutrient use efficiency [17, 68]. Deciduous trees are fast-growing species with cheaper tissue investments and rapid resource returns, but evergreen trees are slow-growing species with expensive tissue investments and slower resource returns on that investment [6, 66]. Deciduous species have an acquisitive strategy, and they accordingly have stronger nutrient resorption efficiency and weaker nutrient resorption proficiency. Deciduous trees use N and P exploitatively and have higher carbon assimilation and transpiration rates, which can improve photosynthetic efficiency and fix nutrients in a short period to enable rapid growth. Thus, deciduous trees are characterized by high SLA, high green leaf N, and P concentration, as well as by cheap but plentiful tissue expenditures (i.e., higher PL and LA) [26, 66]. Conversely, evergreen species adopt a conservative strategy, and so their N and P resorption process is relatively weak. Evergreen trees consume N and P economically, but they allocate expensive resources to build tissue with thick laminas and high tissue densities (i.e., high LT, LDMC, and TWD) to withstand physical and mechanical damage. Meanwhile, they slow growth to preserve nutrients, and they have higher SPAD and CNgr to reduce metabolic costs [26, 71, 72].

While the fundamental role of SLA and green leaf nutrients in mediating N or P resorption has been increasingly recognized, their explicit roles remain to be explored [21, 45, 73]. Here, we found that SLA and green leaf nutrient concentrations are the traits most affecting leaf N and P resorption, matching previous studies that SLA underlined the role of functional traits on leaf N and P utilization [21, 45, 74]. On one hand, while trait correlation networks varied with leaf habit, these traits had markedly higher values of centrality for both deciduous and evergreen trees. On the other hand, they were the most important components of PC1 and PC2 axes. Additionally, in the multivariate linear mixed effects models, SLA, Ngr, and Pgr had significant effects on N and P resorption (Fig. 6). These findings collectively suggested that plant functional traits related to leaf resource capture and defense play dominant roles in mediating the trait coordination of N and P economics. Specifically, we found that SLA, Ngr, and Pgr correlated positively with N and P resorption efficiency and senescent leaf N and P concentrations, which agreed with previous studies [3, 45, 74]. Here, we assumed that in our targeted subtropical forests, evergreen trees showed an adaptive response to low resource availability, and deciduous trees adapted to reduce water loss during winter periods. Evergreen trees produce leaves with a long lifespan that exhibit a slow metabolism, low SLA, and low nutrient concentrations. Relative to evergreen trees, deciduous trees produce short-lived leaves that exhibit traits related to rapid metabolism, i.e., high SLA and high nutrient concentrations. This further accounts for the higher nutrient resorption efficiency of deciduous trees. Finally, we found the notable phenomenon that individual traits displayed a weaker relationship with nutrient resorption than did PES, which suggests that plant functional traits operate in a coordinated manner to produce divergent mechanisms for nutrient utilization under different nutrient conditions. This finding underscores the importance of considering multiple trait dimensions, rather than individual traits one at a time.

Subordinate role of environmental factors in leaf N and P resorption

Both efficiency and proficiency of N and P resorption showed no significant association with elevation, and only PRE and PRP exhibited marginally significant differences when considering the interaction of elevation and leaf habit (Table 1, Fig. S1). Variance partition analysis also demonstrated that environmental factors were less important than plant functional traits in driving variation of leaf N and P resorption (Fig. S2). These findings aligned with our hypothesis that in our target subtropical forests, elevation-related environmental factors directly explain less variations in leaf N and P resorption of broad-leaved tree species than leaf habits-related plant functional traits. Mantel’s test and RDA further verified the negligible effects of environmental factors, including microclimates and soil nutrients, on leaf N and P resorption (Fig. 3, Fig. S3). Broad-scale studies have indicated that N and P resorption is closely related to soil nutrient availability [6, 18, 75]. At first sight, these seem to imply that tree species from nutrient-poor environments resorb nutrients more efficiently, with the important caveat that study species came from habitats differing widely in climate and soil. These tree species have undergone different environmental selection and genetic adaption processes, thus tree species adaptation and acclimation are accompanying [35, 76]. There has been compelling evidence that plant functional traits are more likely to reflect plant nutrient utilization strategies and thus directly explain more to the variation in leaf N and P resorption than climate and soil along elevational gradients [22, 59, 77].

Therefore, we confirm our first hypothesis and argue that environmental factors play a minor direct role in N and P resorption along elevational gradients. In other words, trees alter their N and P resorption along elevational gradients due to changes in tree species composition, essentially by changing plant trait combinations. Notably, similar to previous studies, we sampled only 0–10 cm of topsoil, which, although sufficient in the majority of cases [18], may weaken the importance of soil properties in leaf nutrient resorption. Therefore, it could be useful to sample deeper soil to improve our understanding of soil properties on nutrient resorption in further studies. Moreover, previous studies indicated that the abundance of trees with conservative strategies, which have lower nutrient resorption efficiency, increases at high altitudes in response to the low temperatures that hamper nutrient mineralization by slowing soil microbial activities [20, 62, 73]. This finding highlights that nutrient resorption is not merely a process but a plant nutrient utilization for adapting to diverse environments. Consequently, future research should pay more attention to the trade-offs of above- and below-ground nutrient acquisition pathways when exploring plant trait-resorption linkage to advance our understanding of nutrient cycling in forest ecosystems.

Conclusions

By collecting a greater number of representative tree species than in previous studies, we have simultaneously elucidated the relative contributions of environmental factors and plant functional traits to N and P resorption along an elevational gradient in the subtropical region. We found lower leaf N and P resorption in the subtropical transitional forest than previously assumed, expanding our current knowledge of leaf N and P resorption beyond the relatively well-studied tropical and temperate forests. Deciduous species tend to resorb nutrients more effectively but less proficiently than evergreen counterparts for both N and P. These imply highly diversified nutrient conservation strategies across leaf habits, with evergreen trees preferring conservative strategies that invest in growth and structure and deciduous trees preferring acquisitive strategies for photosynthetic and nutrient resources. Our work further highlighted the overriding role of plant functional traits in the transitional subtropical forest, and that incorporating PES, especially SLA and green leaf nutrient concentrations, into the monitoring and modeling of leaf nutrient resorption could notably improve our understanding of species-driven nutrient cycling. These findings underscored the potential importance of plant functional trait coordination and thus call for these PES effects on leaf nutrient resorption to be more explicitly incorporated into a current nutrient cycling model framework to better forecast the consequences of shifting tree species composition for biogeochemical cycles in forest ecosystems since nutrient resorption has been commonly involved in many biogeochemical models.

Availability of data and materials

The data that support the findings of this study are available from the corresponding author upon reasonable request.

References

  1. Liu B, Fan X, Meng D, Liu Z, Gao D, Chang Q, Bai E. Ectomycorrhizal trees rely on nitrogen resorption less than arbuscular mycorrhizal trees globally. Ecol Lett. 2023;00:1–11.

    Google Scholar 

  2. Brant AN, Chen HYH. Patterns and mechanisms of nutrient resorption in plants. Crit Rev Plant Sci. 2015;34(5):471–86.

    Article  CAS  Google Scholar 

  3. Vergutz L, Manzoni S, Porporato A, Novais RF, Jackson RB. Global resorption efficiencies and concentrations of carbon and nutrients in leaves of terrestrial plants. Ecol Monogr. 2012;82(2):205–20.

    Article  Google Scholar 

  4. Augusto L, Achat DL, Jonard M, Vidal D, Ringeval B. Soil parent material-A major driver of plant nutrient limitations in terrestrial ecosystems. Glob Change Biol. 2017;23(9):3808–24.

    Article  Google Scholar 

  5. Killingbeck KT. Nutrients in senesced leaves: Keys to the search for potential resorption and resorption proficiency. Ecology. 1996;77(6):1716–27.

    Article  Google Scholar 

  6. Achat DL, Pousse N, Nicolas M, Augusto L. Nutrient remobilization in tree foliage as affected by soil nutrients and leaf life span. Ecol Monogr. 2018;88(3):408–28.

    Article  Google Scholar 

  7. Sun X, Li D, Lü X, Fang Y, Ma Z, Wang Z, Chu C, Li M, Chen H. Widespread controls of leaf nutrient resorption by nutrient limitation and stoichiometry. Funct Ecol. 2023;37(6):1653–62.

    Article  CAS  Google Scholar 

  8. Prieto I, Querejeta JI. Simulated climate change decreases nutrient resorption from senescing leaves. Glob Change Biol. 2019;26:1795–807.

    Article  Google Scholar 

  9. Suriyagoda LDB, Ryan MH, Gille CE, Dayrell RLC, Finnegan PM, Ranathunge K, Nicol D, Lambers H. Phosphorus fractions in leaves. New Phytol. 2023;237(4):1122–35.

    Article  CAS  PubMed  Google Scholar 

  10. Cleveland CC, Houlton BZ, Smith WK, Marklein AR, Reed SC, Parton W, Del Grosso SJ, Running SW. Patterns of new versus recycled primary production in the terrestrial biosphere. Proc Natl Acad Sci. 2013;110(31):12733–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Aerts R. Nutrient resorption from senescing leaves of perennials: are there general patterns? J Ecol. 1996;84(4):597–608.

    Article  Google Scholar 

  12. Yuan ZY, Chen HYH. Global-scale patterns of nutrient resorption associated with latitude, temperature and precipitation. Glob Ecol Biogeogr. 2009;18(1):11–8.

    Article  Google Scholar 

  13. Suseela V, Tharayil N, Xing B, Dukes JS. Warming and drought differentially influence the production and resorption of elemental and metabolic nitrogen pools in Quercus rubra. Glob Change Biol. 2015;21(11):4177–95.

    Article  Google Scholar 

  14. Du E, Terrer C, Pellegrini AF, Ahlström A, van Lissa CJ, Zhao X, Xia N, Wu X, Jackson RB. Global patterns of terrestrial nitrogen and phosphorus limitation. Nat Geosci. 2020;13(3):221–6.

    Article  CAS  Google Scholar 

  15. Deng M, Liu L, Jiang L, Liu W, Wang X, Li S, Yang S, Wang B. Ecosystem scale trade-off in nitrogen acquisition pathways. Nature Ecology & Evolution. 2018;2(11):1724–34.

    Article  Google Scholar 

  16. Xu M, Zhu Y, Zhang S, Feng Y, Zhang W, Han X. Global scaling the leaf nitrogen and phosphorus resorption of woody species: Revisiting some commonly held views. Sci Total Environ. 2021;788: 147807.

    Article  CAS  PubMed  Google Scholar 

  17. Aerts R, Chapin FS. The mineral nutrition of wild plants revisited: A re-evaluation of processes and patterns. Adv Ecol Res. 1999;30:1–67.

    Article  Google Scholar 

  18. Yan T, Zhu J, Yang K. Leaf nitrogen and phosphorus resorption of woody species in response to climatic conditions and soil nutrients: a meta-analysis. Journal of Forestry Research. 2018;29(4):905–13.

    Article  CAS  Google Scholar 

  19. Díaz S, Kattge J, Cornelissen JH, Wright IJ, Lavorel S, Dray S, Reu B, Kleyer M, Wirth C, Prentice IC. The global spectrum of plant form and function. Nature. 2016;529(7585):167–71.

    Article  PubMed  Google Scholar 

  20. Sartori K, Violle C, Vile D, Vasseur F, Villemereuil P, Bresson J, Gillespie L, Fletcher LR, Sack L, Kazakou E. Do leaf nitrogen resorption dynamics align with the slow-fast continuum? A test at the intraspecific level. Funct Ecol. 2022;36(5):1315–28.

    Article  CAS  Google Scholar 

  21. Wood TE, Lawrence D, Wells JA. Inter-specific variation in foliar nutrients and resorption of nine canopy-tree species in a secondary neotropical rain forest. Biotropica. 2011;43(5):544–51.

    Article  Google Scholar 

  22. Liu B, Gao D, Chang Q, Liu Z, Fan X, Meng D, Bai E. Leaf enzyme plays a more important role in leaf nitrogen resorption efficiency than soil properties along an elevation gradient. J Ecol. 2022;110(11):2603–14.

    Article  CAS  Google Scholar 

  23. Hättenschwiler S, Aeschlimann B, Coûteaux MM, Roy J, Bonal D. High variation in foliage and leaf litter chemistry among 45 tree species of a neotropical rainforest community. New Phytol. 2008;179(1):165–75.

    Article  PubMed  Google Scholar 

  24. Poorter H, Lambers H, Evans JR. Trait correlation networks: a whole-plant perspective on the recently criticized leaf economic spectrum. New Phytol. 2014;201(2):378–82.

    Article  PubMed  Google Scholar 

  25. Messier J, Lechowicz MJ, McGill BJ, Violle C, Enquist BJ. Interspecific integration of trait dimensions at local scales: the plant phenotype as an integrated network. J Ecol. 2017;105(6):1775–90.

    Article  Google Scholar 

  26. Li J, Chen X, Niklas KJ, Sun J, Wang Z, Zhong Q, Hu D, Cheng D. A whole-plant economics spectrum including bark functional traits for 59 subtropical woody plant species. J Ecol. 2022;110(1):248–61.

    Article  CAS  Google Scholar 

  27. Yu L, Huang Z, Li Z, Korpelainen H, Li C. Sex-specific strategies of nutrient resorption associated with leaf economics in Populus euphratica. J Ecol. 2022;110(9):2062–73.

    Article  CAS  Google Scholar 

  28. Zhao YT, Ali A, Yan ER. The plant economics spectrum is structured by leaf habits and growth forms across subtropical species. Tree Physiol. 2017;37(2):173–85.

    PubMed  Google Scholar 

  29. Weigelt A, Mommer L, Andraczek K, Iversen CM, Bergmann J, Bruelheide H, Fan Y, Freschet GT, Guerrero-Ramírez NR, Kattge J, et al. An integrated framework of plant form and function: the belowground perspective. New Phytol. 2021;232(1):42–59.

    Article  PubMed  Google Scholar 

  30. Joswig JS, Wirth C, Schuman MC, Kattge J, Reu B, Wright IJ, Sippel SD, Rüger N, Richter R, Schaepman ME. Climatic and soil factors explain the two-dimensional spectrum of global plant trait variation. Nature Ecology & Evolution. 2022;6(1):36–50.

    Article  Google Scholar 

  31. Freschet GT, Cornelissen JHC, van Logtestijn RSP, Aerts R. Substantial nutrient resorption from leaves, stems and roots in a subarctic flora: what is the link with other resource economics traits? New Phytol. 2010;186(4):879–89.

    Article  CAS  PubMed  Google Scholar 

  32. Reed SC, Townsend AR, Davidson EA, Cleveland CC. Stoichiometric patterns in foliar nutrient resorption across multiple scales. New Phytol. 2012;196(1):173–80.

    Article  CAS  PubMed  Google Scholar 

  33. Drenovsky RE, Pietrasiak N, Short TH. Global temporal patterns in plant nutrient resorption plasticity. Glob Ecol Biogeogr. 2019;28(6):728–43.

    Article  Google Scholar 

  34. Pepin NC, Arnone E, Gobiet A, Haslinger K, Kotlarski S, Notarnicola C, Palazzi E, Seibert P, Serafin S, Schöner W, et al. Climate changes and their elevational patterns in the mountains of the world. Reviews of Geophysics. 2022;60:e2020RG000730.

    Article  Google Scholar 

  35. Gerdol R, Iacumin P, Brancaleoni L, Wang F. Differential effects of soil chemistry on the foliar resorption of nitrogen and phosphorus across altitudinal gradients. Funct Ecol. 2019;33(7):1351–61.

    Article  Google Scholar 

  36. Ge J, Ma B, Xu W, Zhao C, Xie Z. Temporal shifts in the relative importance of climate and leaf litter traits in driving litter decomposition dynamics in a Chinese transitional mixed forest. Plant Soil. 2022;477(1–2):679–92.

    Article  CAS  Google Scholar 

  37. Xie Z, Shen G. The outstanding universal value and conservation of Hubei Shennongjia. Beijing: Science Press; 2021.

  38. Luo L, Shen G, Xie Z, Yu J. Leaf functional traits of four typical forests along the altitudinal gradients in Mt. Shennongjia Acta Ecologica Sinica. 2011;31(21):6420–8.

    Google Scholar 

  39. The Editorial Board of Flora of China: Flora of China. Beijing: Science Press; 2004.

  40. Zhang L, Luo T, Zhu H, Daly C, Deng K. Leaf life span as a simple predictor of evergreen forest zonation in China. J Biogeogr. 2010;37(1):27–36.

    Article  CAS  Google Scholar 

  41. Ma B, Zhao C, Ge J, Xu W, Xiong G, Shen G, Xie Z. A dataset of 17 dominant plants phenological observation in Shennongjia (2009–2018). China Scientific Data. 2020;5(1):20–30.

    Google Scholar 

  42. Halbritter AH, De Boeck HJ, Eycott AE, Reinsch S, Robinson DA, Vicca S, Berauer B, Christiansen CT, Estiarte M, Grunzweig JM, et al. The handbook for standardized field and laboratory measurements in terrestrial climate change experiments and observational studies (ClimEx). Methods Ecol Evol. 2020;11(1):22–37.

    Article  Google Scholar 

  43. Funk JL, Larson JE, Ames GM, Butterfield BJ, Cavender-Bares J, Firn J, Laughlin DC, Sutton-Grier AE, Williams L, Wright J. Revisiting the Holy Grail: using plant functional traits to understand ecological processes. Biol Rev. 2017;92(2):1156–73.

    Article  PubMed  Google Scholar 

  44. Pérez-Harguindeguy N, Diaz S, Gamier E, Lavorel S, Poorter H, Jaureguiberry P, Bret-Harte M, Comwell W, Craine J, Gurvich D. New handbook for stand-ardised measurement of plant functional traits worldwide. Aust J Bot. 2013;61(3):167–234.

    Article  Google Scholar 

  45. Zhang JL, Zhang SB, Chen YJ, Zhang YP, Poorter L, Bonser S. Nutrient resorption is associated with leaf vein density and growth performance of dipterocarp tree species. J Ecol. 2015;103(3):541–9.

    Article  Google Scholar 

  46. China Ecosystem Research Network Science Committee. Protocols for standard biological observation and measurement in terrestrial ecosystems. Beijing: China Environmental Science Press; 2007.

    Google Scholar 

  47. Pérez-Harguindeguy N, Díaz S, Garnier E, Lavorel S, Poorter H, Jaureguiberry P, Bret-Harte MS, Cornwell WK, Craine JM, Gurvich DE, et al. New handbook for standardised measurement of plant functional traits worldwide. Aust J Bot. 2013;61(3):167–234.

    Article  Google Scholar 

  48. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin P, O’Hara R, Simpson G, Solymos P. vegan: Community ecology package. R package version 2.6–4. 2022.

  49. Jin Y, Qian H. V.PhyloMaker: an R package that can generate very large phylogenies for vascular plants. Ecography. 2019;42(8):1353–9.

    Article  Google Scholar 

  50. Zhang J LB, Liu S, Feng Z, Jiang K. plantlist: Looking up the status of plant scientific names based on the Plant List Database, searching the Chinese names and making checklists of plants. R package version 0.8.0. 2022

  51. Kembel SW, Cowan PD, Helmus MR, Cornwell WK, Morlon H, Ackerly DD, Blomberg SP, Webb CO. Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010;26(11):1463–4.

    Article  CAS  PubMed  Google Scholar 

  52. Revell LJ. Phytooln updated R ecosystem for phylogenetic comparative methods (and other things). PeerJ. 2024;12:e16505.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Csardi G, Nepusz T. The igraph software package for complex network research. InterJournal, Complex Systems. 2006;1695(5):1–9.

    Google Scholar 

  54. Xie J, Wang Z, Li Y. Stomatal opening ratio mediates trait coordinating network adaptation to environmental gradients. New Phytol. 2022;235(3):907–22.

    Article  PubMed  Google Scholar 

  55. Vu VQ, Friendly M. ggbiplot: A grammar of graphics implementation of biplots. R package version 0.6.2. 2024.

  56. Bates D, Mächler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models Using lme4. J Stat Softw. 2015;67(1):1–48.

    Article  Google Scholar 

  57. Wickham H. Package ‘ggplot2’: elegant graphics for data analysis. New York: Springer-Verlag; 2016.

    Google Scholar 

  58. R Core Team: R: A language and environment for statistical computing. Vienna:  R Foundation for Statistical Computing; 2023.

  59. Tang L, Han W, Chen Y, Fang J. Resorption proficiency and efficiency of leaf nutrients in woody plants in eastern China. Journal of Plant Ecology. 2013;6(5):408–17.

    Article  Google Scholar 

  60. Kang H, Xin Z, Berg B, Burgess PJ, Liu Q, Liu Z, Li Z, Liu C. Global pattern of leaf litter nitrogen and phosphorus in woody plants. Ann For Sci. 2010;67(8):811.

    Article  Google Scholar 

  61. Estiarte M, Campiolli M, Mayol M, Penuelas J. Variability and limits in resorption of nitrogen and phosphorus during foliar senescence. Plant Commun. 2023;4(2):100503.

  62. Tsujii Y, Onoda Y, Kitayama K. Phosphorus and nitrogen resorption from different chemical fractions in senescing leaves of tropical tree species on Mount Kinabalu. Borneo Oecologia. 2017;185(2):171–80.

    Article  PubMed  Google Scholar 

  63. McGroddy ME, Daufresne T, Hedin LO. Scaling of C: N: P stoichiometry in forests worldwide: Implications of terrestrial redfield-type ratios. Ecology. 2004;85(9):2390–401.

    Article  Google Scholar 

  64. Craine JM, Mack MC. Nutrients in senesced leaves: comment. Ecology. 1998;79(5):1818–20.

    Google Scholar 

  65. Estiarte M, Peñuelas J. Alteration of the phenology of leaf senescence and fall in winter deciduous species by climate change: effects on nutrient proficiency. Glob Change Biol. 2015;21(3):1005–17.

    Article  Google Scholar 

  66. Reich PB. The world-wide “fast-slow” plant economics spectrum: a traits manifesto. J Ecol. 2014;102(2):275–301.

    Article  Google Scholar 

  67. Powers JS, Tiffin P. Plant functional type classifications in tropical dry forests in Costa Rica: leaf habit versus taxonomic approaches. Funct Ecol. 2010;24(4):927–36.

    Article  Google Scholar 

  68. van Ommen KA, Douma J, Ordonez JC, Reich PB, Van Bodegom P. Global quantification of contrasting leaf life span strategies for deciduous and evergreen species in response to environmental conditions. Glob Ecol Biogeogr. 2012;21(2):224–35.

    Article  Google Scholar 

  69. Keller AB, Phillips RP. Leaf litter decay rates differ between mycorrhizal groups in temperate, but not tropical, forests. New Phytol. 2019;222(1):556–64.

    Article  CAS  PubMed  Google Scholar 

  70. Huang X, Lu Z, Xu X, Wan F, Liao J, Wang J. Global distributions of foliar nitrogen and phosphorus resorption in forest ecosystems. Sci Total Environ. 2023;871: 162075.

    Article  CAS  PubMed  Google Scholar 

  71. Ouédraogo D-Y, Fayolle A, Gourlet-Fleury S, Mortier F, Freycon V, Fauvet N, Rabaud S, Cornu G, Bénédet F, Gillet J-F, et al. The determinants of tropical forest deciduousness: disentangling the effects of rainfall and geology in central Africa. J Ecol. 2016;104(4):924–35.

    Article  Google Scholar 

  72. Kaproth MA, Fredericksen BW, González-Rodríguez A, Hipp AL, Cavender-Bares J. Drought response strategies are coupled with leaf habit in 35 evergreen and deciduous oak (Quercus) species across a climatic gradient in the Americas. New Phytol. 2023;239(3):888–904.

    Article  PubMed  Google Scholar 

  73. González-Zurdo P, Escudero A, Mediavilla S. N resorption efficiency and proficiency in response to winter cold in three evergreen species. Plant Soil. 2015;394(1):87–98.

    Article  Google Scholar 

  74. Wright IJ, Westoby M. Nutrient concentration, resorption and lifespan: leaf traits of Australian sclerophyll species. Funct Ecol. 2003;17(1):10–9.

    Article  Google Scholar 

  75. Yuan Z, Chen HY. Negative effects of fertilization on plant nutrient resorption. Ecology. 2015;96(2):373–80.

    Article  CAS  PubMed  Google Scholar 

  76. Kobe RK, Lepczyk CA, Iyer M. Resorption efficiency decreases with increasing green leaf nutrients in a global data set. Ecology. 2005;86(10):2780–92.

    Article  Google Scholar 

  77. Sigdel SR, Liang E, Rokaya MB, Rai S, Dyola N, Sun J, Zhang L, Zhu H, Chettri N, Chaudhary RP, et al. Functional traits of a plant species fingerprint ecosystem productivity along broad elevational gradients in the Himalayas. Funct Ecol. 2022;37(2):383–94.

    Article  Google Scholar 

Download references

Acknowledgements

We would like to thank Professor Björn Berg at the University of Helsinki for his insightful comments on this manuscript. We are grateful to Dr. Savannah Grace at the University of Florida for her assistance with the English language and grammatical editing of the manuscript.

Funding

The National Natural Science Foundation of China (Grant No.32271641 and No.31600360) financed this study.

Author information

Authors and Affiliations

Authors

Contributions

BM, JG, and ZX conceived and designed the study. BM and JG conducted the fieldwork and lab experiments with contributions from CZ, WX, KX, and ZX. BM conducted the data analyses. BM and JG led the writing of the manuscript. All authors contributed substantially to the final writing.

Corresponding author

Correspondence to Jielin Ge.

Ethics declarations

Ethics approval and consent to participate

The collection of leaf samples was permitted by Xingshan County Forestry Bureau, Yichang, Hubei Province, China. We declare that this study complies with the guidelines and legislation of the People’s Republic of China, the IUCN Policy Statement on Research Involving Species at Risk of Extinction, and the Convention on the Trade in Endangered Species of Wild Fauna and Flora.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ma, B., Ge, J., Zhao, C. et al. Plant economics spectrum governs leaf nitrogen and phosphorus resorption in subtropical transitional forests. BMC Plant Biol 24, 764 (2024). https://doi.org/10.1186/s12870-024-05484-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05484-9

Keywords