“Diminishing Returns” in the Scaling Between Leaf Area and Twig Size in Three Forest Communities Along an Elevation Gradient of Wuyi Mountain, China

: Background and aims: The “diminishing returns” hypothesis postulates that the scaling exponent governing the lamina area versus lamina mass scaling relationships has, on average, a numerical value less than one. Theoretically, a similar scaling relationship may exist at the twig level. However, this possibility has not been explored empirically. Methods: We tested both hypotheses by measuring the lamina area and mass, petiole mass of individual leaves, and the total foliage area and stem mass of individual current ‐ year shoots (twigs) of 64 woody species growing in three characteristic forest community types:(1) Evergreen broad ‐ leaved, (2) mixed coniferous and broad ‐ leaved, and (3) deciduous. Key results: The results demonstrate that lamina area vs. mass and lamina area vs. petiole mass differ significantly among the three forest types at both the individual leaf and twig levels. Nevertheless, the scaling exponents of lamina area vs. mass were <1.0 in each of the three community types, as were the corresponding exponents for lamina area vs. petiole mass, both within and across the three community types. Similar trends were observed at the individual twig level. The numerical values of the scaling exponent for lamina area vs. petiole mass and total foliage area vs. stem mass per twig decreased with increased elevation. Conclusions : These data support the “diminishing returns” hypothesis at both the individual leaf level and at the individual twig level, phenomena that can inform future inquiries into the mechanistic basis of biomass allocation patterns to physiological (leaf) and mechanical (stem) plant organs. and stem mass from the current ‐ year shoots (i.e., twigs) of 64 woody plant species collected from three different types of forests (i.e., evergreen, mixed coniferous and broad ‐ leaved, and deciduous) along an elevation gradient in China. Several key environmental conditions (i.e., light intensity, temperature, and nutrient availabilities) were also measured and recorded. The data were used to (1) test whether the “diminishing returns” hypothesis holds at both the individual leaf and twig levels, and (2) determine the scaling relationship between lamina area and total supporting mass (the sun of petiole and stem mass) at the twig level. with simple biomechanical considerations. These findings can serve toinformfuture inquiries into the mechanisticbasis of twig biomass allocation patterns between


Introduction
Leaf area varies by almost six orders of magnitude across species, ranging from less than 1.0 mm 2 tomorethan 1.0 m 2 [1]. As the primary light-harvesting organ for the majority of species, it is importantto understand, both empirically and theoretically, how variation in leaf size (measured as area or mass) affects the capacity to capture solar energy and how this in turn affects ecosystem functions [2]. One approach is to analyze the size-dependent (scaling) relationships between leaf lamina area (A) with respect to leaf mass and the mass (M)of supporting structures, such as the petiole and stem, to determine whether they conform to the general scaling equation A= β M α , where β is the normalization constant and α is the scaling exponent [3][4][5].
In this context, the "diminishing returns" hypothesis predicts that, at the individual leaf level, the interspecific scaling exponent governing the relationship between lamina area (al) and lamina mass (ml), on average, will be less than unity (i.e.,α < 1.0), such that incremental gains in lamina mass fail to produce proportional increases in lamina area [5]. If true, specific leaf area (SLA =al/mL), which is a key index of photosynthetic efficiency [6,7], is predicted to decrease with increasing lamina mass across species. Furthermore, it is reasonable to suppose that a larger lamina area or mass would be associated with a more massive petiole, as reported by some early empirical work [8]. For example, using data collected from 44 species, Niinemets et al. [9] found that larger leaves tend to have disproportionately more massive petioles compared with smaller leaves. Pickup et al. [10]reached a similar conclusion. Taken together, the "diminishing returns "hypothesis and empirical data support the notion that increases in lamina area will fail to keep pace with increasing lamina and petiole mass (mp) at the individual leaf level.
Lamina area and lamina mass are crucial leaf traits underpinning the metabolic performance of plants [11]. The "diminishing returns" economic principle provides one explanation for why the scaling of total plant growth as a function of total mass is governed by a scaling exponent with a value significantly less than one [12].Although the "diminishing returns" hypothesis has been tested at the individual leaf level within and among species [5,13], it remains unclear whether this hypothesis holds for plants growing along major environmental gradients or different community types. For example, the scaling exponent of al vs. ml can shift from<1.0 to >1.0 with increasing elevation, as shown by Pan et al. [14], whereas the corresponding scaling exponent for bamboos is <1.0andmore or less constant along an elevation gradient, as reported by Sun et al. [11]. Evidently, more comprehensive field data are needed to assess the generalizable validity of the "diminishing returns" hypothesis.
Another mostly untested aspect of this hypothesis is whether it holds true at the individual twig level, i.e., does total foliage area per twig scale at <1.0 with respect to total lamina and petiole mass (ML and MP, respectively)? In this regard, Niklas et al. [5]reported that total foliage area (AL) per twig is proportional to both the leaf number per twig (NL) and mean lamina area of individual leaves (i.e., AL  NL × al), and that ML is also proportional to NL and the mean lamina mass of individual leaves (i.e., ML  NL × ml). If individual lamina area scales disproportionately (allometrically) with individual lamina mass (on average) as proposed by the "diminishing returns "hypothesis (i.e., alml α< 1.0 ), it is not unreasonable to speculate that total foliage area at the twig level could be recast as AL NL × ml α< 1.0 ML α< 1.0 . Similar reasoning leads to the conclusion that a negative scaling relationship exists between AL and total petiole mass (MP). Yet, to the best of our knowledge, there has been no systematic, extensive testing of the "diminishing returns "hypothesis conducted simultaneously at the individual leaf level and the individual twig level for different plant community types along an elevation gradient.
In order to fill this gap in our knowledge, we report on the lamina size (i.e., measured as lamina surface area and mass), petiole mass, and stem mass from the current-year shoots (i.e., twigs) of 64 woody plant species collected from three different types of forests (i.e., evergreen, mixed coniferous and broad-leaved, and deciduous) along an elevation gradient in China. Several key environmental conditions (i.e., light intensity, temperature, and nutrient availabilities) were also measured and recorded. The data were used to (1) test whether the "diminishing returns" hypothesis holds at both the individual leaf and twig levels, and (2) determine the scaling relationship between lamina area and total supporting mass (the sun of petiole and stem mass) at the twig level.

Site Description
The study site is located in the National Natural Reserve of the Wuyi Mountains, southeastern China (27°48.11′-28°00.35′N, 117°39.30′-117°55.47′E). The reserve has a typical subtropical monsoon climate, featuring abundant precipitation and abundant solar radiation. Mean annual precipitation is 2583 mm, with a mean annual temperature of 14.2°C and evaporation of 778mm, and 72-92% relative humidity (yearly average). The Wuyi Mountains are characterized by vertical vegetation zonation, going from evergreen forest at low elevation to mountainous steppe at high elevation. The area belongs to the mid-subtropical monsoon climate, with an annual average temperature of 13.2~14.8 °C, an extreme minimum temperature of −14.2 °C, an extreme maximum temperature of 36.3 °C, and an increase in altitude. The annual, monthly, and seasonal average temperatures decrease with the annual average. The vertical decline rate of temperature is 0.44 °C/100m, that is, for every 100m rise, the annual average temperature drops by 0.44 °C. The average annual precipitation is 2583 mm, and the period from April to June is the precipitation concentration period in the area. The annual evaporation is 778 mm, and the annual average relative humidity is 72-92%. The annual sunshine hours are 774-1144 h, and the annual average frost-free period is 231 days [15].The major soil types are classified as mountain yellow-red soil (400-600 m a.s.l.), mountain yellow soil (600-1300 m a.s.l.), mountain dark yellow-brown soil (1300-1900 m a.s.l.), and mountain meadow soil (above 1900 m a.s.l.) [16].

Leaf and Twig Sampling
Based on the distribution pattern of forest ecosystems at Wuyi Mountain, three types of forest discernable along the elevation gradient were selected: Evergreen broad-leaved forest (EF; 1319m a.s.l.), a mixed forest (MF;1697 m a.s.l.), and a deciduous forest (DF; 1818 m a.s.l.).A total of 32, 20, and 23 species (including overlapping species) were respectively sampled in EF, MF, and DF ( Table  1). The dominant species in the EF were Rhododendron simiarumHance, Schimasuperba, Cyclobalanopsisglauca, Rhododendron.ovatum (Lindl.) Planch. ex Maxim., and Symplocossumuntia; the MF was dominated by Symplocossumuntia, Cyclobalanopsismultiervis, Tsugachinensis, Taxuschinensis, Acerelegantulum, and Illiciumminwanense; and the DF was dominated by Clethrabarbinervis, Photiniabeauverdiana, Acernikoense, and Fraxinuschinensis [17].The stand densities were 3033 trees/hm 2 ,1133 trees/hm 2 , and 2725trees/hm 2 for EF, MF, and DF, respectively, with corresponding mean plant heights of 7.87 m, 10.56 m, and 6.94 m. Three20 m ×20 m plots were established in each forest type, with at least a 20m spacing between plots. In each plot, all healthy woody plants with a diameter at breast height (DBH)>5 cm were marked and their species, DBH, and height recorded. Three to five representative individuals of each plant species were selected, from which five current-year twigs of medium size with mature leaves (but lacking reproductive organs) were randomly selected at the periphery and two sides of each plant's crown. Twigs were removed and placed into sealable plastic bags with their ID-code and brought to the laboratory to separate laminas, petioles, and stems. The length and diameter of each stem were measured with a ruler and a vernier caliper, respectively, and the total number of leaves counted. All individual leaves on twigs were digitally scanned using an Epson scanner (Epson V39, Epson, Japan). After separating the organ components, the fresh lamina mass and area, petiole mass, stem mass, and quantitative characteristics of the twig parts and accessories were measured. All leaves on a given twig were removed to measure total foliage area and mass at the twig level. Immediately afterwards, the laminas, petioles, and stems were oven-dried at 75°C to constant weight to determine their respective dry mass. Each scanned lamina was analyzed usingImageJv1.2 software (National Institute of Health, Bethesda, ML), and its lamina area (al) was recorded.

Data Analysis
The scaling relationships among lamina area and mass, petiole mass, stem mass, and total supporting mass (MSTM; the sum of petiole and stem mass) at the individual leaf and twig level were analyzed using the equation Y1= βY2 α , whereY1and Y2are the interdependent variables of interest, α is the scaling exponent, and β is the normalization constant. All equations were log10-transformed to yield logY1= logβ+ αlogY2. Scaling relationships were quantified using the "smatr" package v 3.4-3 [18] in the R v3.2.3 software platform (R Development Core Team2012; R Foundation for Statistical Computing, Vienna, Austria). The significance level used for testing slope heterogeneity was set a priori at p≤ 0.05 (i.e., slope heterogeneity was rejected for p-values>0.05).

Variation in Lamina Size, Petiole and Stem Mass at theLeaf and Twig Levels
Across the three forest types, at the individual leaf level, mean values of petiole mass (mp), lamina mass (ml),and lamina area (al) individual leaf were 0.010g, 0.168g, and 2228mm 2 , respectively. The mp differed significantly among the three forest communities, but not al, and individual leaf mass (mL)decreased significantly as elevation increased (Table 2).
At the twig level, the mean values of ML, petiole mass (MP), MS, AL, and MSTM were, respectively, 0.824g, 0.052g, 0.107 g, 10994 mm 2 , and 0.154g across the three forest types ( Table 3).The ML of the three forest types declined with elevation, being significantly lowest in DF, while the AL in the EF and DF significantly exceeded that of the MF (Table 2). Table 2.Leaf-level and twig-level biomass allocation (mean ± standard error (SE)) in different forest typesonWuyi Mountain (China).

Type
Leaf-level Twig-level

Thelamina Area vs.Lamina Mass Scaling Relationship
Our results indicate that the numerical value of the scaling exponent governing the relationships between lamina area (al) and lamina mass (ml) is less than 1.0 along an elevational gradient in the evergreen and deciduous forestsof Wuyi Mountain and across all three forest types at the individual leaf level. Therefore,the data supporta key prediction of the "diminishing returns" hypothesis [5,19], i.e., increases in lamina area fail to keep pace with increasesin lamina mass at the individual leaf level.We also found that the scaling exponent foralversus mlis less than 1.0 in the highest elevation forest (the DF), a result inconsistent with Pan et al. [14], who report a scaling exponent>1.0 at a comparableelevation.One plausible explanation for this discrepancyis differences betweenthe species composition and the precipitation and temperatureassociated with the two differentenvironments,both of which can influence the scaling relationship betweenlamina area and lamina mass [5,20,21].
Perhaps more importantly, our analyses demonstrate that the scaling exponent for whole lamina area versus total foliage massis less than 1.0 at the level ofcurrent-year twigs, both within and across the three forest types (Table4; Figure 2A). With increasingleaf vascular lignifications, the dry-weight investmenton a per-leaf area basis should increase, yielding alower SLA for large-sized leaves [6,7]. If this is a general phenomenon, it would implythat leaf area, as an investment trait for capturing more light, would gradually fail to keep pace withincreasing leaf biomass. Indeed, some empirical results indicatethat whole-plant leaf area scales as a 0.90 power law with respect towhole-plant leaf mass in a young Hinoki cypress forest [22].Likewise, it is known that net primary production (NPP) in older forest stands declines [23]. It is possible that the diminishing returns hypothesis helps to explain whythe efficiency of light harvesting decreases as total leaf biomass increases,such that NPPdeclines as forests age. Clearly, more data areneeded to verify the link between diminishing returnsand declining NPP in forest ecosystems.

Scaling Relationships ofLamina Area andPetiole Mass
Petiole investment hasan important influence on biomass allocation atthe twig level [9,[24][25][26].For example, previous studies have shown that petiole investment is size-dependent [9], varying from ca. 2 to 30%, and by even more than 60% in the largest leaves [8,9,25]. However, the scaling relationships between lamina area and petiole mass haveseldom been determined for woody plants at both the within-leaf and within-twig level.
Once again, our data are consistent with asize-dependent (scaling) relationship between lamina area andpetiole mass predicted by the diminishing returns hypothesis, i.e., increases in lamina area (al) fail to keep pace with increases in petiole biomass(mp) atboth the individual leaf and twig levels,both within and acrossdifferent forest types (Tables 3 and 4; Figures1 and 2). These results are interesting from a biomechanical as well as a metabolic perspective, becausea petiole mechanically behaves as an elastic cantilevered beam that must support the static load of the leaf lamina [24]. The capacity of a cantilevered beam to static loads becomes less efficient as the beam becomes longer [27]. Therefore, with increasing lamina area, a disproportionately larger investment must be provided to provide mechanical support [20]. Similarly, fluid mechanics shows that the dynamic drag force induced by wind pressure is proportional to lamina area [24,27] and that this force is transmitted to the petiole [20,24,27,28]. Taken together, these two facts help to explain why petiole mass increases disproportionately with respect to lamina area and mass.Indeed, consistent with our results, prior studies have shown that the scaling exponent relating lamina area and petiole mass is significantly <1.0 in different subtropical forest types [26,29]. Our data similarly show that lamina area fails to keep pace with increases in total supporting biomass (i.e., the sum of stem and petiole mass) within and across the three forest communities examined in this study.

Scaling Relationshipsof Petiole and Stem Biomass Along anElevation Gradient
Previous studieshave generallyneglected the effects of abiotic factors on structural biomass allocation patterns. In the present study, we examined the scaling relationship between total petiole mass and stem mass (i.e., MP and MS) at the whole twig level across three forest communities. We have shown that this relationshipshiftssystematically, in that the numerical value of the scaling exponent decreases with increasing elevation (Table4; Figure 2E), i.e., for any given petiole biomass, woody plants in the evergreen and mixed forests allocate more stem biomass than those in the deciduous forest.
A likely explanation for this result is theadaptationsof deciduous versus non-deciduous species to changes in precipitation and temperature along anelevationalgradient. Arguably, it is reasonable to speculate that the risk of hydraulic failurewill increase (and access of soil nutrients will decrease) with increasing elevation [30,31]. A case in point is the deciduous forest, which is distributed at high elevationsin our study region. Within this region, water, temperature, and other stress factors measurablyincrease. Species adapted to this ecosystem tend to have reduced leaf biomass and increased investments in supporting biomass, presumably to maintain a functioning nutrient and hydraulic transport system [32]. In addition, twigs also require additional investment in stem biomassto provide adequate support tocope with wind and snow mechanical loadings [24]. Otherwise, the high-altitude site is often windy inwinters and suffers from winter snow cover. Windy environmentsnecessitate a high investment in stems or petioles, as suggested by Niklas [24,28].Yet, interpreting the relationship between leaf biomass and stem biomass at the twig level along anelevation gradient must be done cautiously. For example, Yang et al. [33] report that leaf biomass scales is ometrically with respect to stem biomass along anelevationalgradient.In this particular case, the scalingrelationship between lamina and petiole mass in tandem with the systematic shift in petiole versusstem masscan conceivably lead to variationsin theinvestment of lamina area in different forest types, even with a constant biomass allocation between leaf and stem biomass.

Conclusions
We conclude that there is a diminishing return between increase biomass and leaf area as elevation increases. The predictions of the diminishing returns hypothesisappear tohold for twigs as well as leaves across diversewoody plants species. Furthermore, biomass investments in petiole and stem tissues varywith elevation in a manner that is consistent with the predictions of the diminishing returns hypothesis when juxtaposed with simple biomechanical considerations. These findings can serve toinformfuture inquiries into the mechanisticbasis of twig biomass allocation patterns between metabolically active tissuesand structurallysupportive tissues, and future attemptsto model twig ecophysiological traits.
We anticipatethat biomassinvestment in supporting (petiole and stem) structures will be shown to be responsive to environmental gradientsand asleaf area increases within and across different species.More comprehensive data sets,for many co-occurring species under a variety of environmental conditions, are still needed to assessthe generality of the "diminishing returns" hypothesis in plant scaling theory.