The Leaf Economics Spectrum Constrains Phenotypic Plasticity Across a Light Gradient

The leaf economics spectrum (LES) characterizes multivariate correlations that confine the global diversity of leaf functional traits onto a single axis of variation. Although LES is well established for traits of sun leaves, it is unclear how well LES characterizes the diversity of traits for shade leaves. Here, we evaluate LES using the sun and shade leaves of 75 woody species sampled at the extremes of a within-canopy light gradient in a subtropical forest. Shading significantly decreased the mean values of LMA and the rates of photosynthesis and dark respiration, but had no discernable effect on nitrogen and phosphorus content. Sun and shade leaves manifested the same relationships among N_mass, P_mass, A_mass, and R_mass (i.e., the slopes of log–log scaling relations of LES traits did not differ between sun and shade leaves). However, the difference between the normalization constants of shade and sun leaves was correlated with functional trait plasticity. Although the generality of this finding should be evaluated further using larger datasets comprising more phylogenetically diverse taxa and biomes, these findings support a unified LES across shade as well as sun leaves.

Introduction

A central premise of trait-based ecology is that performance and functioning across levels of organization, from organs to ecosystems, are predictable from organ-level functional traits (Niklas, 1994; Westoby et al., 2002; Hallik et al., 2009; Adler et al., 2014; Wullschleger et al., 2014; Fisher et al., 2015; Christoffersen et al., 2016; Díaz et al., 2016; McDowell et al., 2018). The leaf economics spectrum (LES) offers substantial potential for quantifying the broad diversity of leaf functional traits to predict changes in community composition in response to global change, as well as how these changes in composition will feedback to ecosystem function (Reich et al., 1998; Wright et al., 2004). However, the data used to construct the LES is primarily based on the functional traits of sun leaves (i.e., leaves exposed to direct sunlight) and it is unclear whether the trends established by the LES extend to shade leaves.

The LES characterizes evolutionary convergence on robust and consistent relationships among the form, function, chemistry, and longevity of leaves from diverse ecosystems and climates, confining this diversity onto a single axis of variation. Specifically, the LES describes trade-offs among leaf mass per area (LMA), nitrogen and phosphorus contents per unit mass (N_mass and P_mass, respectively), assimilation and respiration rates per unit mass (A_mass and R_mass, respectively), and a number of other critically important functional traits that characterize plant physiology and ecology. The LES is shaped by the joint effects of environmental filtering and biophysical constraints on leaf carbon economics, yielding emergent functional tradeoffs between assimilation rate, leaf longevity, and leaf construction costs (Chabot and Hicks, 1982; Westoby et al., 2002; Kikuzawa and Lechowicz, 2006; Michaletz et al., 2015, 2016; Anderegg et al., 2018). These trade-offs have been demonstrated primarily using global-scale data, and thereby provide a general mechanistic framework for understanding and predicting variations of plant community structure and function across resource availability gradients. For example, under low-light conditions, the LES predicts that leaves with higher LMA will have lower N_mass, P_mass, and A_mass, but longer life spans (Reich et al., 1998; Wright et al., 2004). At the other end of the spectrum, leaves exposed to high light conditions are predicted to have lower LMA and higher N_mass, P_mass, R_mass, and A_mass, but shorter lifespans (Chabot and Hicks, 1982; Mooney and Gulmon, 1982; Merino et al., 1984; Wright et al., 2004; Kikuzawa et al., 2013; Reich, 2014).

Despite the apparent broad applicability of the LES, a major concern is that its predictions are primarily based on data derived from sun leaves (e.g., Wright et al., 2004; Keenan and Niinemets, 2016). As noted previously, it remains uncertain as to whether the trends manifest in the LES hold equally well for shade leaves. Light availability varies more than 50-fold within dense canopies and more than 10-fold within open canopies (Hirose et al., 1997; Koike et al., 2001; Rambal, 2001; Valladares, 2003; Niinemets and Anten, 2009; Niinemets et al., 2015). Light plays a key role in leaf growth and physiology. One of the most profound effects of light on plant development is the suppression of the expression of a series of light-responsive genes (Li et al., 1995). Therefore, sun leaves might exhibit better trait full genotypic expression than shade leaves. This might be one of the reasons for why large-scale LES analyses avoided shade leaves. It is crucial to evaluate whether LES expands to shade leaves, including those from low light environments commonly found in within-canopy light gradients.

Furthermore, as formalized by the carbon economics theory (Chabot and Hicks, 1982; Blonder et al., 2011; Michaletz et al., 2015, 2016), leaf functional traits are expected to vary in response to within-canopy light gradients in order to maximize leaf carbon return on investment. Indeed, it has long been known that many leaf functional traits exhibit plasticity and/or acclimation across within-canopy light gradients (Iio et al., 2005; Lloyd et al., 2010; Kattge et al., 2011; Kenzo et al., 2015). For example, LMA and nitrogen content per area generally decrease with decreasing irradiance (Hirose and Werger, 1987; Givnish, 1988; Hollinger, 1996; Pons and Anten, 2004; Posada et al., 2009; Fajardo and Siefert, 2016), which helps to maximize leaf net carbon gain because of an accompanied lower carbon cost (Anten, 2016). In general, because of the reallocation of nitrogen, leaf N content (the “investment”) along light gradients within plants does not agree with theoretical expectations that N declines with decreasing light availability from the top toward the bottom of the canopy (Hirose and Werger, 1987; Givnish, 1988; Anten, 2016). Hence, large light-dependent sensitivity in LMA appears to be responsible for the observed differences in the photosynthetic production and nutrient content of woody plants (Niinemets et al., 2015; Messier et al., 2016). However, relatively little attention has been paid to LES of shade leaves. Using a worldwide database (spanning 167 species with LMA measurements, 75 species with N_mass measurements, and 44 species with A_mass measurements), Keenan and Niinemets (2016) concluded that LES, thought to comprise data for sun leaves only, might instead include a large proportion of data for leaves growing under shaded conditions, implicitly suggesting that LES may apply generally to all leaves.

Each trait underlying the LES may be affected by canopy conditions and exhibit plasticity (Messier et al., 2010; Sandel et al., 2010). This suggests that values of a certain LES trait in sun and shade leaves might be consistently correlated across all species but that the scaling exponents and normalization constants of the relationship might differ numerically (Wright et al., 2001; Wright and Sutton-Grier, 2012; Niklas and Hammond, 2019). For example, the numerical values of the scaling exponents of some LES relationships have been shown to differ for plants growing in environments with contrasting water and nutrient availabilities (Wright et al., 2001). In addition, Wright and Sutton-Grier (2012) have suggested that trends in the traits of individual species do not consistently predict those of other species. Therefore, it is unclear whether the scaling exponents and normalization constants of LES trait-trait relationships are consistent between sun and shade leaves.

Here, we evaluate the effects of within-canopy shading on LES. Specifically, we examine relationships among five key leaf functional traits (i.e., LMA, N_mass, P_mass, A_mass, and R_mass) to address two key questions. First, to what degree are leaf functional traits plastic between sun- and shade- leaves? Second, do LES relationships exist for both sun- and shade-leaves? And, if they do, are the scaling relationships between sun and shade leaves consistent? We hypothesize that (i) shade induced morphological plasticity is a stronger determinant of leaf acclimation to light availability compared to chemical (N_mass and P_mass) or physiological plasticity (A_mass and R_mass) (see Vaz et al., 2010; Messier et al., 2016); and (ii) LES relationships should hold for shade leaves, but the numerical values of the scaling exponents and normalization constants of some trait scaling relationships may differ between sun and shade leaves.

Materials and Methods

Study Site

This study was conducted in the Yangjifeng National Nature Reserve of Jiangxi, southeastern China (117°11′30″∼117°28′40″ E, 27°51′10″∼ 28°02′20″ N, 1540.9 m a.s.l.). The reserve has a typical humid mid-subtropical monsoon climate. The range of annual precipitation average is 1870–2191 mm, 72–75% of which falls from April to June. The range of average annual temperature is 11.4–18.5°C. The soil is a haplic lixisol. In August 2015, a 500 m × 500 m plot was established. The tree density was 703.5 trees ha^–1, mean diameter at breast height (DBH) was 11.77 ± 0.10 cm (mean ± SE), mean height was 8.54 ± 0.06 m, mean crown base height was 2.85 ± 0.02 m, and mean crown diameter was 3.18 ± 0.02 m. The dominant species were Castanopsis fargesii, Alniphyllum fortune, Litsea cubeba, Castanopsis carlesii, Elaeocarpus sylvestris, and Schima superba.

Gas-Exchange Measurements

Branches for gas exchange measurements were collected in August 2017. Seventy-five species, from 61 genera and 33 families, were selected based on abundance (representing 88.12% of individuals with DBH ≥ 1 cm in the plot; Supplementary Table S1). Three individuals were sampled from each species. To minimize the confounding effects of leaf age (Wright et al., 2006) and branch length (Niklas and Cobb, 2010), three current-year branches (twigs) were selected from each of the canopy top (sun leaves) and canopy bottom (shade leaves) for each individual following the protocol of Sack et al. (2006). Branches were collected from canopies and immediately placed in water to reduce water loss. Embolisms were removed by re-cutting branch ends under water (Yoder et al., 1994; Mori et al., 2010; Michaletz et al., 2016; Kelly et al., 2017; Malhi et al., 2017). There were no significant differences between branch traits from sun and shade canopy positions (P > 0.05, t-tests) (Supplementary Table S2). For each sampled tree, the photosynthetic photon flux densities (PPFDs) of leaves at the sampled canopy locations were measured on a sunny day with a LI-250A light meter (LICOR, Lincoln, NE, United States). All the PPFDs of the sun leaves (1172.25 ± 25.73 μmol m^–2 s^–1) were significantly higher than that in the shade leaves (82.88 ± 3.23 μmol m^–2 s^–1) (Supplementary Figure S1A). The PPFDs at the shade leaves consisted of about 7.15% ± 0.36 (Supplementary Figure S1A) of full sunlight in the sun leaves.

Gas-exchange measurements were conducted on six to nine leaves from each branch. Leaves showing no evidence of disease or herbivory were used for all measurements and kept hydrated throughout the time measurements were taken and recorded. Maximum net assimilation rates (A_area; μmol m^–2 s^–1) and respiration rates (R_area) were estimated from light response curves obtained using a LI-6800 portable photosynthesis system and calculated by Photosynthesis (LICOR, Lincoln, NE, United States). During measurements, leaf temperatures were maintained at 25°C, chamber CO₂ was maintained at 400 p.p.m, and flow rates were maintained at 500 mmol s^–1. To obtain light response curves, net assimilation rates were measured at PPFDs of 2000, 1600, 1200, 900, 600, 300, 200, 100, 80, 50, 20, and 0 μmol m^–2 s^–1 (Supplementary Figure S1B).

Leaf Mass per Area, Leaf Nitrogen, and Leaf Phosphorus Content Measurements

All of the leaves from each branch were detached and scanned using a flatbed scanner (Epson V39, Epson, Japan), yielding a total of data from 7417 leaves. Leaf area was measured from scanned images using Image J (National Institute of Health, Bethesda, ML, United States). Leaf samples were dried subsequently at 75°C for 48 h before measuring dry weight. Samples were then ground and passed through a 100-mesh sieve (0.15 mm). The leaf nitrogen content per unit mass (N_mass,%) was determined using an Element Analyzer (VARIO EL III Element Analyzer, Elementar, Germany). The leaf phosphorus content per unit mass (P_mass,%) was measured using the molydate/ascorbic acid method and a continuous flow analyzer (SKALAR SAN++, Netherlands) after H₂SO₄-HClO₄ (4:1, v:v) digestion.

Data Analysis

Leaf Trait Plasticity

To evaluate whether the interspecies traits of sun leaves are significantly different than those of shade leaves, traits were compared using t-tests and linear regression analysis. For each trait, an index of maximum within-canopy plasticity was calculated as the quotient for sun leaves divided by that of shade leaves following the protocol of Sack et al. (2006) after calculating the average of a trait in sun and shade leaves at three taxonomic levels (i.e., species, genus, and family). Linear mixed effect models were fit using the lme4 packages in R (v.3.5.1; R Foundation for Statistical Computing, Vienna, Austria) to determine the extent to which the three different taxonomic levels contributed to trait variation. Analyses showed that there were no significant differences of leaf trait plasticity among the three taxonomic levels (Table 1 and Supplementary Table S3) and that among-individual and among-species differences together consistently explained over 50% of the total trait variance with the exception LMA (i.e., 48.57%) (Supplementary Figure S2). Thus, leaf trait plasticity was expressed as the average plasticity of all species.

TABLE 1

Table 1. Summary of traits for sun and shade leaves from 75 tree species in a subtropical forest.

PCA Analysis

To check how sun and shade leaf traits variability is organized, a principle component analysis (PCA) was conducted using species as data points. All variables were log-transformed before analysis to improve normality. PCAs were performed for the five traits using the ape (Paradis et al., 2004) and vegan (Dixon, 2003) packages in R (v.3.5.1; R Foundation for Statistical Computing, Vienna, Austria). Bivariate correlations between each trait and PCA axis 1 (PC1) and axis 2 (PC2) were evaluated using Pearson’s correlation coefficient (r).

Pairwise Scaling Relationships

We examined bivariate scaling relationships between LMA, N_mass, P_mass, A_mass, and R_mass, and compared our results with those of a published dataset of global trait variation (Wright et al., 2004). The formula y = βx^α (Equation 1) was used to examine bivariate trait relationships. To linearize this relationship, Equation (1) was log-transformed to give the linear form log (y) = log (β) + α log (x) (Equation 2), where y and x are interdependent leaf functional traits of interest, β is the normalization constant, and α is the scaling exponent (the slope of the log-log linear bivariate plot). The linear formula was then fit to data using model II standardized major axis (SMA) regression (Warton et al., 2012) using the smatr package 3.4-3 in R (v.3.5.1; R Foundation for Statistical Computing, Vienna, Austria). When no significant numerical differences between two or more scaling exponents (P > 0.05) for comparable bivariate comparisons were observed (e.g., sun vs. shade leaves), a common scaling exponent was then estimated.

Results

Variation of Traits for Sun and Shade Leaves

Within a given canopy location (sun or shade), we observed more than one order of magnitude variation in each of the five leaf traits (Table 1). The most variable trait was A_mass, with a 14-fold range in sun leaves and a 7-fold range in shade leaves (Table 1 and Figure 1). In contrast, the least variable trait was N_mass, with a 4-fold range in sun leaves and a 1-fold range in shade leaves (Table 1 and Figure 1).

FIGURE 1

Figure 1. Histogram with normal distribution and shade-sun trait relationships for functional traits of sun and shade leaves from the same canopy. The dashed lines indicate 1:1. (A) Leaf mass per area (LMA); (B) mass-based net assimilation rate (A_mass); (C) mass-based respiratory rate (R_mass); (D) mass-based nitrogen content (N_mass); (E) mass-based phosphorus content (P_mass). *indicates a significant difference at P < 0.05 and ***indicates a significant difference at P < 0.001.

Shading had different effects on each of the leaf functional traits (Table 1 and Figure 1). The LMA, A_mass, and R_mass of shade leaves were significantly lower compared with those of sun leaves (Figure 1). Specifically, the mean values of A_mass and R_mass respectively were 113.79 nmol g^–1 s^–1 and 9.50 nmol g^–1 s^–1 in sun leaves, 66.51 nmol g^–1 s^–1 and 6.25 in nmol g^–1 s^–1 in shade leaves (Table 1). The plasticity of A_mass and R_mass was 1.80 ± 0.07 and 1.63 ± 0.09, respectively. The mean value of LMA decreased from 83.68 g m^–2 in sun leaves to 75.25 g m^–2 in shade. The plasticity for LMA was 1.13 ± 0.03 (Table 1 and Figure 1). In contrast, no significant differences between sun and shade leaves were detected for N_mass and P_mass, although the plasticity of N_mass (i.e., 1.02 ± 0.01) was somewhat greater than that of P_mass (i.e., 0.98 ± 0.02) (Table 1).

PCAs of Sun and Shade Leaf Traits

The PC1 axis captured 73.32 and 69.24% of the total variation in the five traits of sun and shade leaves, respectively, with all traits contributing substantially to this axis (Table 2). The PC1 and PC2 axes captured 81.04% of the total variability across all species (Table 2 and Figure 2). When sun and shade leaves were combined, PC1 captured 67.90% of the total variation in the five leaf traits, with all traits contributing substantially to this axis (Table 2 and Figure 2).

TABLE 2

Table 2. Bivariate relationships between individual traits and the scores of the first and second PC for sun leaves, shade leaves and whole-canopy economic spectra.

FIGURE 2

Figure 2. Principal component analysis (PCA) of morphological, chemical, and physiological traits of sun and shade leaves for 75 woody species. Sun and shade group was separated by ellipsoid at 68% normal probability. The PC1 scores between sun and shade leaves were not significantly different (P > 0.05).

All traits contributed substantially to the PC1 axis (Table 2). Leaf N_mass and P_mass and physiological rates were negatively correlated with PC1, whereas LMA was positively correlated with PC1 (Table 2).

Scaling Relationships of Sun and Shade Leaf Traits

Consistent with sun leaves, robust correlations among leaf functional traits were also observed in shade leaves (Table 3, 4). With the exception of LMA vs. N_mass and LMA vs. P_mass, the numerical values of scaling exponents were significantly different from ± 1.0 (Table 4), indicating that traits manifest allometric rather than isometric scaling relationships. The scaling exponents of all the bivariate trait relationships did not differ significantly between the sun and shade leaves (Table 4 and Figure 3). Indeed, generalizable scaling relationships were evident within each of the bivariate trait relationships. In the particular case of photosynthetic rates, the scaling exponents governing the N content for sun and shade leaves were 1.38 and 2.03 (common slope = 1.67), and 1.41 and 2.07 (common slope = 1.71), respectively (Table 3 and Figure 3). However, shading decreased the numerical values of the normalization constants for physiological metabolic rates vs. chemical traits and LMA scaling relationships, as well as decreased the numerical values of the normalization constants of A_mass and R_mass (Table 4 and Figure 3). In contrast, the normalization constants between N_mass (and P_mass) vs. LMA as well as N_mass vs. P_mass were unaffected for leaves collected in the shade (Table 4 and Figure 3).

TABLE 3

Table 3. Scaling exponents for bivariate relationships among mass-based traits of sun or shade leaves.

TABLE 4

Table 4. Scaling exponents for bivariate relationships among mass-based traits of sun and shade leaves.

FIGURE 3

Figure 3. Leaf Economics Spectrum (LES) scaling relationships for sun and shade leaves from this study compared with those from the global dataset of Wright et al. (2004). The numbers in parenthesis indicate 95% confidence intervals of scaling exponents. ns indicates no significant difference between intercepts of sun and shade leaves, ** indicates a significant difference at P < 0.01, and *** indicates a significant difference at P < 0.001.

The difference of the normalization constants between sun and shade leaves was significantly positively correlated with the differences of plasticity of the related functional traits (Figure 4).

FIGURE 4

Figure 4. Bivariate plot of the difference of the normalization constants between sun and shade leaves against the difference of plasticity of two related functional traits. D-value: the difference value.

Discussion

The goal of this paper was to evaluate the universality of LES relationships in the context of within-canopy light gradients. This goal was approached in two ways. First, we examined the leaf functional trait plasticity for sun and shade leaves. Second, we evaluated whether species-level differences in trait plasticity affected the canonical correlations predicted by the LES. In general, our results demonstrate only minimal to moderate variation in traits across leaves collected from different canopy light gradients (i.e., the relationships between sun and shade leaves scaling relationships do not statistically significantly differ), indicating that phenotypic trait variation across species is constrained by convergent leaf economic spectra.

We observed moderate levels of leaf functional trait plasticity across species within-canopy light gradients, indicating that leaf functional traits across a wide range of light conditions appear to be in equilibrium with the environment. These results are consistent with those of previous studies. For example, Sack et al. (2003) reported a limited LMA plasticity for mixed Quercus forests in southern Spain (i.e., average within-canopy plasticity was 1.5 or less). Wright et al. (2006) quantified the leaf N and P contents for 28 woody species along vertical canopy profiles and found that P_mass varied in a manner consistent with optimal allocation, whereas LMA and N_mass decreased with canopy depth in only two and five species, respectively. Iio et al. (2005) found that there was no significant difference in N_mass between the best-lit and the most shaded leaves within the crown of a Fagus crenata tree. We also found that shading had the most influence on leaf physiological traits (A_mass and R_mass) followed by morphological traits (i.e., 1.13 ± 0.04 times), but shading had no disenable effect on chemical (N and P) traits (Table 1). We suggest that these modest phenotypic patterns may reflect adaptations to maximize photosynthetic carbon gain.

As predicted, regardless of whether the data were drawn from sun or shade leaves, species at the fast-return end of the LES had high leaf nutrient contents, high A_mass and R_mass, and low LMA. At the slow-return end of the spectrum, species had high LMA, low nutrient contents, and low A_mass and R_mass (Figure 2). Many other studies have reported that trends in the LES apply equally well to species growing in resource-poor environments. For example, Liu et al. (2010) found that correlated pairs of leaf and root traits held equally well for species adapted to semi-arid to arid environments. Correlations among leaves traits in forest understory ferns have also been shown to hold true for patterns reported for seed plants (Karst and Lechowicz, 2006). Therefore, resource availability appears to have only modest effects on LES trends (Wright et al., 2004).

It is noteworthy that sun and shade leaves do not depart from each other on the PC1 axis (Figure 2), indicating that trait plasticity does not change LES relationships between the two leaf-types, thereby supporting the expectations of the leaf-level carbon economics theory (Chabot and Hicks, 1982; Blonder et al., 2011; Michaletz et al., 2015, 2016). In contrast, Zhao et al. (2016) found that sun-acclimated species had a greater capacity to transport nutrients and higher rates of resource acquisition, whereas the reverse was true for shade- acclimated species manifesting slow-growth strategies. It is understandable that the economics strategies of sun and shade leaves will differ among different plant functional groups, such as tree and shrub species. Ordoñez et al. (2010) have demonstrated that qualitative attributes, such as growth form, woodiness, and leaf habit, can differ in accord with different environmental factors. In this context, our data show that deciduous and evergreen species cluster on the opposing ends of the LES likely because of their different acquisition strategies (Supplementary Table S4).

The LES does not consider the effects of within-canopy light availability on whole-plant performance. It implicitly assumes that within-canopy light gradients do not alter general scaling relationships. Our analyses support this supposition. The data show that the numerical values of the scaling exponents of key LES relationships are not statistically different between sun and shade leaves (Table 4 and Figure 3). For example, the scaling exponents for the A_mass vs. N_mass are numerically indistinguishable in sun and shade leaves (with a common slope of 1.69, 95% CIs = 1.48 and 1.93, P > 0.05). These similarities indicate that the LES predicts coordinated shifts among traits due to the fundamental biological trade-offs that underpin the LES and constrains leaf phenotypic plasticity across the subtropical forest canopy light gradients examined in this study (Reich et al., 2003; Shipley et al., 2006). Specifically, differences in species-level responses to changes in light conditions may not be sufficiently different to alter the relative positions of species along trait axes. This finding is consistent with the study of Keenan and Niinemets (2016) who report that the scaling relationships between LMA and both A_mass and Nmass are similar in sun and shade leaves (see Figure 3 in Keenan and Niinemets, 2016). Furthermore, the scaling exponents between key leaf functional traits of shade and sun leaves do not differ numerically, even when evergreen and deciduous woody species are examined separately (Supplementary Table S4).

The numerical values of the scaling exponents for each of the leaf morphology vs. physiological scaling relationship in our study are all numerically larger than those in the world-wide data set, and fall within the lower region of the world-wide pattern of LMA (Table 4 and Figure 3), which consists of data drawn from plants growing at comparatively high temperatures in lower latitudes (Poorter et al., 2010). For the world-wide sun leaf dataset, the scaling relationship for LMA vs. N_mass (Table 4) is negatively allometric, i.e., a = -1.28 (95% CIs = -1.32 and 1.24 (Wright et al., 2004). Extensive N deposition tends to increase leaf N content (Peñuelas et al., 2015), which could explain the numerical difference observed between our data and those using the world-wide data set. The scaling exponents of the N_mass vs. P_mass scaling relationship for sun and shade leaves (i.e., 0.77 and 0.78, respectively) are similar to those reported by Niklas and Cobb (2005) who report a leaf N_mass vs. P_mass scaling exponent approximately equal to 0.75.

However, the scaling exponent of N_mass vs. P_mass in this study is numerically greater than that in the world-wide data set (with a slope 0.66) (Table 4 and Figure 3). A number of factors may help to explain this difference, particularly since our data are drawn from species living in a subtropical community. For example, Tian et al. (2017) report large variations in leaf N_mass vs. P_mass scaling relationships across latitudinal zones, wherein scaling exponents tend to numerically increase from boreal to tropical regions. In addition, the growing season and leaf longevity are longer in subtropical than in boreal regions, and leaf growth rates and P-demands are greater during the growing season with increasing latitude (Yan et al., 2017). Furthermore, soil P availability relative to N availability tends to increase with increasing latitude or from humid to arid regions, resulting in decreasing leaf N:P ratios (Han et al., 2005). These relationships help to explain at least in part why the leaf N_mass vs. P_mass scaling exponent reported here is numerically greater than that observed using a world-wide spectrum of data.

The N:P ratio in our study was also higher (18.37 ± 0.36) than that reported by Han et al. (2005) (i.e., N:P = 14.4) for plants across China, or He et al. (2008) for temperate and cold grasses in North China (i.e., N:P = 15.3). Most studies have concluded that P rather than N is the main factor limiting plant growth in the subtropical region (e.g., Reich et al., 1997; Han et al., 2005), which is consistent with the higher N:P ratio observed for the leaves examined in our study. This is also consistent with the larger scaling exponents observed for the LMA and A_mass vs. P_mass relationships (Table 4 and Figure 3), and the lower R_mass reported here compared to those for the world-wide data set (Table 4 and Figure 3).

Of particular interest in this regard is the observation that, although vastly different canopy irradiance gradient exhibit similar scaling exponents for any given trait relationship, the numerical values of the normalization constants of LES relationships differ significantly (Figure 3). These results may reflect quantifiable differences in species evolutionary and individual life histories (e.g., light) (Niklas and Hammond, 2019). Although the normalization constant is as biologically meaningful as the scaling exponent, little attention has been paid to how or why the normalization constant differs across data sets, or how it changes during the course of evolution by natural selection (Niklas and Hammond, 2019). In our study, the differences in the normalization constants of traits between sun and shade leaves were positively correlated with the differences in their related plastic index (Figure 4). This suggests that the variations in normalization constants might be attributable to differences in trait plasticity. N was allocated to metabolic functions in sun leaves and then a leaf with very high net gas exchange characteristics was constructed. Whereas shade leaves appeared to show excess uptake of N, which could be used for future use or allocated for non-photosynthesis (Chapin, 1980; Ågren, 2008). Therefore, at any given N_mass, A_mass, and R_mass were higher in sun leaves than in shade leaves. Furthermore, the constants for A_mass vs. R_mass and leaf morphology vs. N_mass are not significantly different for sun and shade leaves for either evergreen or deciduous woody species (Supplementary Table S5), indicating that differences in normalization constants for A_mass vs. R_mass and leaf morphology vs. N_mass across species are likely caused by differences of leaf habits among species in our data set. It is important to note that this inferences is drawn based on data drawn from subtropical communities. Because our study deals with only a relatively small community, more studies in different biomes are needed to confirm the generality of our findings.

The effects of within-canopy light gradients on leaf traits have been incorporated into canopy process (and land surface) models to predict canopy carbon flux and forest production (e.g., Ryu et al., 2011; Bonan et al., 2012). However, only some dynamic global vegetation models concerned with photosynthesis at the leaf level explicitly or implicitly consider the effects of shading within-canopy (Prentice et al., 2007). Our data contribute to a general LES for both sun and shade leaves and can inform future attempts at modeling whole canopy processes, particularly attempts to couple leaf photosynthesis and position within-canopy.

Conclusion

Although the light availability gradient throughout canopies significantly affects physiological and morphological leaf traits, the scaling exponents of pairwise functional trait relationships for five mass-based leaf traits in sun and shade leaves do not numerically differ. These results support the proposition that a “canonical” LES exists and can be applied to shade leaves. Future progress toward understanding the scaling relationships elucidated by the LES nevertheless requires additional data with which to assess the effects of environmental factors on leaf traits, because it is becoming increasingly clear that the LES relationships are dependent on the species composition of communities.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Author Contributions

XC, MW, QZ, and DC designed the research. XC, JS, and ML performed the research and analyzed the data. XC, MW, KN, SM, and DC wrote the manuscript. All authors read and approved the final manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (31722007), the National Key Research and Development Program of China (2017YFC0505400), and the Fujian Natural Science Funds for Distinguished Young Scholars (2018J07003).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00735/full#supplementary-material

References

Adler, P. B., Fajardo, A., Kleinhesselink, A. R., and Kraft, N. J. B. (2014). Trait-based tests of coexistence mechanisms. Ecol. Lett. 16, 1294–1306. doi: 10.1111/ele.12157

PubMed Abstract | CrossRef Full Text | Google Scholar

Ågren, G. I. (2008). Stoichiometry and nutrition of plant growth in natural communities. Annu. Rev. Ecol. Evol. Syst. 39, 153–170. doi: 10.1146/annurev.ecolsys.39.110707.173515

CrossRef Full Text | Google Scholar

Anderegg, L. D. L., Berner, L. T., Badgley, G., Sethi, M. L., Law, B. E., HilleRisLambers, J., et al. (2018). Within-species patterns challenge our understanding of the leaf economics spectrum. Ecol. Lett. 215, 734–744. doi: 10.1111/ele.12945

PubMed Abstract | CrossRef Full Text | Google Scholar

Anten, N. P. R. (2016). “Optimization and game theory in canopy models,” in Canopy Photosynthesis: From Basics to Applications, eds K. Hikosaka, Ü Niinemets, and N. P. R. Anten (Dordrecht: Springer Netherlands), 355–377. doi: 10.1007/978-94-017-7291-4_13

CrossRef Full Text | Google Scholar

Blonder, B., Violle, C., Bentley, L. P., and Enquist, B. J. (2011). Venation networks and the origin of the leaf economics spectrum. Ecol. Lett. 14, 91–100. doi: 10.1111/j.1461-0248.2010.01554.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonan, G. B., Oleson, K. W., Fisher, R. A., Lasslop, G., and Reichstein, M. (2012). Reconciling leaf physiological traits and canopy flux data: use of the TRY and FLUXNET databases in the Community Land Model version 4. J. Geophys Res. Biogeol. 117:G02026. doi: 10.1029/2011JG001913

CrossRef Full Text | Google Scholar

Chabot, B. F., and Hicks, D. J. (1982). The ecology of leaf life spans. Annu. Rev. Ecol. Syst. 13, 229–259. doi: 10.1146/annurev.es.13.110182.001305

CrossRef Full Text | Google Scholar

Chapin, F. S. I. (1980). The mineral nutrition of wild plants. Annu. Rev. Ecol. Evol. Syst. 11, 233–260. doi: 10.1146/annurev.es.11.110180.001313

CrossRef Full Text | Google Scholar

Christoffersen, B. O., Gloor, M., Fauset, S., Fyllas, N. M., Galbraith, D. R., Baker, T. R., et al. (2016). Linking hydraulic traits to tropical forest function in a size-structured and trait-driven model (TFS v. 1-Hydro). Geosci. Mod. Dev. 9, 4227–4255. doi: 10.5194/gmd-9-4227-2016

CrossRef Full Text | Google Scholar

Díaz, S., Kattge, J., Cornelissen, J. H., Wright, I. J., Lavorel, S., Dray, S., et al. (2016). The global spectrum of plant form and function. Nature 529, 167–171. doi: 10.1038/nature16489

PubMed Abstract | CrossRef Full Text | Google Scholar

Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x

CrossRef Full Text | Google Scholar

Fajardo, A., and Siefert, A. (2016). Temperate rain forest species partition fine-scale gradients in light availability based on their leaf mass per area (LMA). Ann. Bot. 118, 1307–1315. doi: 10.1093/aob/mcw184

PubMed Abstract | CrossRef Full Text | Google Scholar

Fisher, R. A., Muszala, S., Verteinstein, M., Lawrence, P., Xu, C., McDowell, N. G., et al. (2015). Taking off the training wheels: the properties of a dynamic vegetation model without climate envelopes. Geosci. Model Dev. 8, 3293–3357. doi: 10.5194/gmdd-8-3293-2015

CrossRef Full Text | Google Scholar

Givnish, T. J. (1988). Adaptation to sun and shade: a whole-plant perspective. Funct. Plant Biol. 15, 63–92. doi: 10.1071/PP9880063

CrossRef Full Text | Google Scholar

Hallik, L., Kull, O., Nilson, T., and Peñuelas, J. (2009). Spectral reflectance of multispecies herbaceous and moss canopies in the boreal forest understory and open field. Can. J. Remote. Sens. 35, 474–485. doi: 10.5589/m09-040

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, W., Fang, J., Guo, D., and Zhang, Y. (2005). Leaf Nitrogen and Phosphorus Stoichiometry across 753 Terrestrial Plant Species in China. New Phytol. 168, 377–385. doi: 10.1111/j.1469-8137.2005.01530.x

PubMed Abstract | CrossRef Full Text | Google Scholar

He, J., Wang, L., Flynn, D., Wang, X., Ma, W., and Fang, J. (2008). Leaf nitrogen: phosphorus stoichiometry across Chinese grassland biomes. Oecologia 155, 301–310. doi: 10.1007/s00442-007-0912-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirose, T., Ackerly, D. D., Traw, M. B., Ramseier, D., Bazzaz, F. A., et al. (1997). CO2 elevation, canopy photosynthesis, and optimal leaf area index. Ecology 78, 2339–2350. doi: 10.2307/2265896

CrossRef Full Text | Google Scholar

Hirose, T., and Werger, M. J. A. (1987). Maximizing daily canopy photosynthesis with respect to the leaf nitrogen allocation pattern in the canopy. Oecologia 72, 520–526. doi: 10.1007/BF00378977

PubMed Abstract | CrossRef Full Text | Google Scholar

Hollinger, D. Y. (1996). Optimality and nitrogen allocation in a tree canopy. Tree Physiol. 16, 627–634. doi: 10.1093/treephys/16.7.627

PubMed Abstract | CrossRef Full Text | Google Scholar

Iio, A., Fukasawa, H., Nose, Y., Kato, S., and Kakubari, Y. (2005). Vertical, horizontal and azimuthal variations in leaf photosynthetic characteristics within a Fagus crenata crown in relation to light acclimation. Tree Physiol. 25, 533–544. doi: 10.1093/treephys/25.5.533

PubMed Abstract | CrossRef Full Text | Google Scholar

Karst, A. L., and Lechowicz, M. J. (2006). Are correlations among foliar traits in ferns consistent with those in the seed plants? New Phytol. 173, 306–312. doi: 10.1111/j.1469-8137.2006.01914.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kattge, J., Díaz, S., Lavoreli, S., Prentice, C., Leadley, P., Boenisch, G., et al. (2011). TRY - a global database on plant traits. Glob. Change Biol. 17, 2905–2935.

Google Scholar

Keenan, T. F., and Niinemets, Ü. (2016). Global leaf trait estimates biased due to plasticity in the shade. Nat. Plants 3:16201. doi: 10.1038/nplants.2016.201

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelly, G., Lugassi, N., Belausov, E., Wolf, D., Khamaisi, B., Brandsma, D., et al. (2017). The Solanum tuberosum KST1 partial promoter as a tool for guard cell expression in multiple plant species. J. Exp. Bot. 68, 2885–2897. doi: 10.1093/jxb/erx159

PubMed Abstract | CrossRef Full Text | Google Scholar

Kenzo, T., Inoue, Y., Yoshimura, M., Yamashita, M., Tanaka-Oda, A., Ichie, T., et al. (2015). Height-related changes in leaf photosynthetic traits in diverse Bornean tropical rain forest trees. Oecologia 1771, 191–202. doi: 10.1007/s00442-014-3126-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kikuzawa, K., and Lechowicz, M. J. (2006). Toward synthesis of relationships among leaf longevity, instantaneous photosynthetic rate, lifetime leaf carbon gain, and the gross primary production of forests. Am. Nat. 168, 373–383. doi: 10.1086/506954

PubMed Abstract | CrossRef Full Text | Google Scholar

Kikuzawa, K., Onoda, Y., Wright, I. J., and Reich, P. B. (2013). Mechanisms underlying global temperature-related patterns in leaf longevity. Glob. Ecol. Biogeogr. 22, 982–993. doi: 10.1111/geb.12042

CrossRef Full Text | Google Scholar

Koike, T., Kitao, M., Maruyama, Y., and Lei, T. T. (2001). Leaf morphology and photosynthetic adjustments among deciduous broad-leaved trees within the vertical canopy profile. Tree Physiol. 21, 951–958. doi: 10.1093/treephys/21.12-13.951

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H. M., Culligan, K., and Chory, D. J. (1995). CUE1: a mesophyll cell-specific positive regulator of light-controlled gene expression in arabidopsis. Plant Cell. 7, 1599–1610. doi: 10.2307/3870022

CrossRef Full Text | Google Scholar

Liu, G., Freschet, G. T., Pan, X., Cornelissen, J. H., Li, Y., Dong, M., et al. (2010). Coordinated variation in leaf and root traits across multiple spatial scales in Chinese semi-arid and arid ecosystems. New Phytol. 188, 543–553. doi: 10.1111/j.1469-8137.2010.03388.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lloyd, J., Patiño, S., Paiva, R. Q., Nardoto, G. B., Quesada, C. A., Santos, A. J. B., et al. (2010). Optimisation of photosynthetic carbon gain and within-canopy gradients of associated foliar traits for Amazon forest trees. Biogeosciences 7, 1833–1859. doi: 10.5194/bg-7-1833-2010

CrossRef Full Text | Google Scholar

Malhi, Y., Girardin, C. A. J., Goldsmith, G. R., Doughty, C. E., Salinas, N., Metcalfe, D. B., et al. (2017). The variation of productivity and its allocation along a tropical elevation gradient: a whole carbon budget perspective. New Phytol. 21, 1019–1032. doi: 10.1111/nph.14189

PubMed Abstract | CrossRef Full Text | Google Scholar

McDowell, N. G., Michaletz, S. T., Bennett, K. E., Solander, K. C., Xu, C., Maxwell, R. M., et al. (2018). Predicting chronic climate-driven disturbances and their mitigation. Trends Ecol. Evol. 33, 15–27. doi: 10.1016/j.tree.2017.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Merino, J., Field, C., and Mooney, H. A. (1984). Construction and maintenance costs of mediterranean-climate evergreen and deciduous leaves I. Growth and CO2 exchange analysis. Oecologia 53, 208–213. doi: 10.1007/BF00545665

PubMed Abstract | CrossRef Full Text | Google Scholar

Messier, J., Mcgill, B. J., Enquist, B. J., and Lechowicz, M. J. (2016). Trait variation and integration across scales: is the leaf economic spectrum present at local scales? Ecography 13, 685–697. doi: 10.1111/ecog.02006

CrossRef Full Text | Google Scholar

Messier, J., McGill, B. J., and Lechowicz, M. J. (2010). How do traits vary across ecological scales? A case for trait-based ecology. Ecol. Lett. 13, 838–848. doi: 10.1111/j.1461-0248.2010.01476.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Michaletz, S. T., Weiser, M. D., Mcdowell, N. G., Zhou, J., Kaspari, M., Helliker, B. R., et al. (2016). The energetic and carbon economic origins of leaf thermoregulation. Nat. Plants 2:16129. doi: 10.1038/nplants.2016.129

PubMed Abstract | CrossRef Full Text | Google Scholar

Michaletz, S. T., Weiser, M. D., Zhou, J., Kaspari, M., Helliker, B. R., Enquist, B. J., et al. (2015). Plant thermoregulation: energetics, trait-environment interactions, and carbon economics. Trends Ecol. Evol. 30, 714–724. doi: 10.1016/j.tree.2015.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Mooney, H. A., and Gulmon, S. L. (1982). Constraints on leaf structure and function in reference to herbivory. Bioscience 32, 198–206. doi: 10.2307/1308943

CrossRef Full Text | Google Scholar

Mori, S., Yamaji, K., Ishida, A., Prokushkin, S. G., Masyagina, O. V., Hagihara, A., et al. (2010). Mixed-power scaling of whole-plant respiration from seedlings to giant trees. Proc. Natl. Acad. Sci. U.S.A. 107, 1447–1451. doi: 10.1073/pnas.0902554107

PubMed Abstract | CrossRef Full Text | Google Scholar

Niinemets, Ü, and Anten, N. P. R. (2009). Photosynthesis In Silico: Understanding Complexity from Molecules to Ecosystems, eds A. Laisk, L. Nedbal, and J. Govindjee (Berlin: Springer Verlag), 363–399.

Google Scholar

Niinemets, Ü, Keenan, T. F., and Hallik, L. A. (2015). A worldwide analysis of within-canopy variations in leaf structural, chemical and physiological traits across plant functional types. New Phytol. 205, 973–993. doi: 10.1111/nph.13096

PubMed Abstract | CrossRef Full Text | Google Scholar

Niklas, K. J. (1994). Plant Allometry: the Scaling of Form and Process. Chicago, IL: University of Chigo Press, doi: 10.1016/S0169-5347(00)89015-1

CrossRef Full Text | Google Scholar

Niklas, K. J., and Cobb, E. D. (2005). N, P, and C stoichiometry of Eranthis hyemalis (Ranunculaceae) and the allometry of plant growth. Am. J. Bot. 92, 1256–1263. doi: 10.3732/ajb.92.8.1256

PubMed Abstract | CrossRef Full Text | Google Scholar

Niklas, K. J., and Cobb, E. D. (2010). Ontogenetic changes in the numbers of short- vs. long-shoots account for decreasing specific leaf area in acer rubrum (aceraceae) as trees increase in size. Am. J. Bot. 97, 27–37. doi: 10.3732/ajb.0900249

PubMed Abstract | CrossRef Full Text | Google Scholar

Niklas, K. J., and Hammond, S. T. (2019). On the interpretation of the normalization constant in the scaling equation. Front. Ecol. Evol. 6:212. doi: 10.3389/fevo.2018.00212

CrossRef Full Text | Google Scholar

Ordoñez, J. C., Van Bodegom, P. M., Witte, J. P. M., Bartholomeus, R. P., van Dobben, H. F., et al. (2010). Leaf habit and woodiness regulate different leaf economy traits at a given nutrient supply. Ecology 91, 3218–3228. doi: 10.1890/09-1509.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Paradis, E., Claude, J., and Strimmer, K. (2004). APE: analyses of Phylogenetics and Evolution in R language. Bioinformatics 20, 289–290. doi: 10.1093/bioinformatics/btg412

PubMed Abstract | CrossRef Full Text | Google Scholar

Peñuelas, J., Sardans, J., Rivas-Ubach, A., and Janssens, I. A. (2015). The human-induced imbalance between C, N and P in Earth’s life system. Glob. Change Biol. 18, 3–6. doi: 10.1111/j.1365-2486.2011.02568.x

CrossRef Full Text | Google Scholar

Pons, T. L., and Anten, N. P. R. (2004). Is plasticity in partitioning of photosynthetic resources between and within leaves important for whole-plant carbon gain in canopies? Funct. Ecol. 18, 802–811. doi: 10.1111/j.0269-8463.2004.00895.x

CrossRef Full Text | Google Scholar

Poorter, H., Niinemets, Ü, Poorter, L., Wright, I. J., and Villar, R. (2010). Causes and consequences of variation in leaf mass per area (lma): a meta-analysis. New Phytol. 182, 565–588. doi: 10.1111/j.1469-8137.2009.02830.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Posada, J. M., Lechowicz, M. J., and Kitajima, K. (2009). Optimal photosynthetic use of light by tropical tree crowns achieved by adjustment of individual leaf angles and nitrogen content. Ann. Bot. 103, 795–805. doi: 10.1093/aob/mcn265

PubMed Abstract | CrossRef Full Text | Google Scholar

Prentice, I. C., Bondeau, A., Cramer, W., Harrison, S. P., Hickler, T., Lucht, W., et al. (2007). Dynamic Global Vegetation Modeling: Quantifying Terrestrial Ecosystem Responses to Large-Scale Environmental Change in Terrestrial Ecosystems in a Changing World. Berlin: Springer Berlin Heidelberg, 175–192. doi: 10.1007/978-3-540-32730-1_15

CrossRef Full Text | Google Scholar

Rambal, S. (2001). “Productivity of mediterranean-type ecosystems,” in Terrestrial Global Productivity: Past, Present, and Future, eds H. A. Mooney, B. Saugier, and J. Roy (San Diego, CA: Academic Press, Inc), 315–344. doi: 10.1046/j1442-9993.2002.12041.x

CrossRef Full Text | Google Scholar

Reich, P. B. (2014). The world-wide ‘fast–slow’ plant economics spectrum: a traits manifesto. J. Ecol. 102, 275–301. doi: 10.1111/1365-2745.12211

CrossRef Full Text | Google Scholar

Reich, P. B., Walters, M. B., and Ellsworth, D. S. (1997). From tropics to tundra: global convergence in plant functioning. Proc. Natl. Acad. Sci. U.S.A. 94, 13730–13734. doi: 10.1073/pnas.94.25.13730

PubMed Abstract | CrossRef Full Text | Google Scholar

Reich, P. B., Walters, M. B., Ellsworth, D. S., Vose, J. M., Volin, J. C., Gresham, C., et al. (1998). Relationships of leaf dark respiration to leaf nitrogen, specific leaf area and leaf life-span: a test across biomes and functional groups. Oecologia 114, 471–482. doi: 10.1007/s004420050471

PubMed Abstract | CrossRef Full Text | Google Scholar

Reich, P. B., Wright, I. J., Cavender-Bares, J., Craine, J. M., Oleksyn, J., Westoby, M., et al. (2003). The evolution of plant functional variation: traits, spectra, and strategies. Int. J. Plant. Sci. 164, 143–164. doi: 10.1086/374368

CrossRef Full Text | Google Scholar

Ryu, Y., Baldocchi, D. D., Kobayashi, H., van Ingen, C., Li, J., Black, T. A., et al. (2011). Integration of MODIS land and atmosphere products with a coupled-process model to estimate gross primary productivity and evapotranspiration from 1 km to global scales. G. Bio. C 25:Gb4017. doi: 10.1029/2011GB004053

CrossRef Full Text | Google Scholar

Sack, L., Grubb, P. J., and Marañón, T. (2003). The functional morphology of juvenile plants tolerant of strong summer drought in shaded forest understories in southern Spain. Plant Ecol. 168, 139–163. doi: 10.2307/20146470

CrossRef Full Text | Google Scholar

Sack, L., Melcher, P. J., Liu, W. H., Middleton, E., and Pardee, T. (2006). How strong is intracanopy leaf plasticity in temperate deciduous trees? Am. J. Bot. 93, 829–839. doi: 10.3732/ajb.93.6.829

PubMed Abstract | CrossRef Full Text | Google Scholar

Sandel, B., Goldstein, L. J., Kraft, N. J., Okie, J. G., Shuldman, M. I, Ackerly, D. D., et al. (2010). Contrasting trait responses in plant communities to experimental and geographic variation in precipitation. New Phytol. 188, 565–575. doi: 10.1111/j.1469-8137.2010.03382.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shipley, B., Lechowicz, M. J., Wright, I., and Reich, P. B. (2006). Fundamental trade-offs generating the worldwide leaf economics spectrum. Ecology 87, 535–541. doi: 10.1890/05-1051

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, D., Yan, Z., Niklas, K. J., Han, W., Kattge, J., Reich, P. B., et al. (2017). Global leaf nitrogen and phosphorus stoichiometry and their scaling exponent. Natl. Sci. Rev. 5, 728–739. doi: 10.1093/nsr/nwx142

CrossRef Full Text | Google Scholar

Valladares, F. (2003). “Light heterogeneity and plants: from ecophysiology to species coexistence and biodiversity,” in Progress in Botany, eds K. Esser, U. Luttge, W. Beyschlag, and F. Hellwig (Berlin, Germany: Springer Verlag), 439–471.

Google Scholar

Vaz, M., Pereira, J. S., Gazarini, L. C., David, T. S., David, J. S., Rodrigues, A., et al. (2010). Drought-induced photosynthetic inhibition and autumn recovery in two Mediterranean oak species (Quercus ilex and Quercus suber). Tree Physiol. 30, 946–956. doi: 10.1093/treephys/tpq044

PubMed Abstract | CrossRef Full Text | Google Scholar

Warton, D. I, Duursma, R. A., Falster, D. S., and Taskinen, S. (2012). Smatr 3– an R package for estimation and inference about allometric lines. Methods Ecol. Evol. 3, 257–259. doi: 10.1111/j.2041-210X.2011.00153.x

CrossRef Full Text | Google Scholar

Westoby, M., Falster, D. S., Moles, A. T., Vesk, P. A., and Wright, I. J. (2002). Plant ecological strategies: some leading dimensions of variation between species. Annu. Rev. Ecol. Syst. 33, 125–159. doi: 10.1146/annurev.ecolsys.33.010802.150452

CrossRef Full Text | Google Scholar

Wright, I. J., Leishman, M. R., Read, C., Read, C., and Westoby, M. (2006). Gradients of light availability and leaf traits with leaf age and canopy position in 28 Australian shrubs and trees. Funct. Plant Biol. 33, 407–419. doi: 10.1071/FP05319

CrossRef Full Text | Google Scholar

Wright, I. J., Reich, P. B., and Westoby, M. (2001). Strategy shifts in leaf physiology, structure and nutrient content between species of high- and low-rainfall and high- and low-nutrient habitats. Funct. Ecol. 15, 423–434. doi: 10.2307/826662

CrossRef Full Text | Google Scholar

Wright, I. J., Reich, P. B., Westoby, M., Ackerly, D. D., Baruch, Z., Bongers, F., et al. (2004). The worldwide leaf economics spectrum. Nature 428, 821–827. doi: 10.1038/nature02403

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, J. P., and Sutton-Grier, A. (2012). Does the leaf economic spectrum hold within local species pools across varying environmental conditions? Funct. Ecol. 26, 1390–1398. doi: 10.1111/1365-2435.12001

CrossRef Full Text | Google Scholar

Wullschleger, S. D., Epstein, H. E., Box, E. O., Euskirchen, E. S., Goswami, S., Iversen, C. M., et al. (2014). Plant functional types in Earth system models: past experiences and future directions for application of dynamic vegetation models in high-latitude ecosystems. Ann. Bot. 114, 1–16. doi: 10.1093/aob/mcu077

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, Z., Tian, D., Han, W., Tang, Z., Fang, J., et al. (2017). An assessment on the uncertainty of the nitrogen to phosphorus ratio as a threshold for nutrient limitation in plants. Ann. Bot. 120, 937–942. doi: 10.1093/aob/mcx106

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoder, B. J., Ryan, M. G., Waring, R. H., Schoettle, A. W., Kaufmann, M. R., et al. (1994). Evidence of reduced photosynthetic rates in old trees. For. Sci. 40, 513–527. doi: 10.1093/forestscience/40.3.513

CrossRef Full Text | Google Scholar

Zhao, Y. T., Ali, A., and Yan, E. R. (2016). The plant economics spectrum is structured by leaf habits and growth forms across subtropical species. Tree Physiol. 37, 173–185. doi: 10.1093/treephys/tpw098

PubMed Abstract | CrossRef Full Text | Google Scholar