# Why Does Not the Leaf Weight-Area Allometry of Bamboos Follow the 3/2-Power Law?

^{1}Co-Innovation Centre for Sustainable Forestry in Southern China, Bamboo Research Institute, Nanjing Forestry University, Nanjing, China^{2}Department of Mathematical Sciences, Centre for Invasion Biology, African Institute for Mathematical Sciences, Stellenbosch University, Matieland, South Africa^{3}Center for Integrative Conservation, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Mengla, China^{4}Western Triangle Agricultural Research Centre, Montana State University, Conrad, MT, United States^{5}Department of Biosciences Engineering, University of Antwerp, Antwerp, Belgium

The principle of similarity (Thompson, 1917) states that the weight of an organism follows the 3/2-power law of its surface area and is proportional to its volume on the condition that the density is constant. However, the allometric relationship between leaf weight and leaf area has been reported to greatly deviate from the 3/2-power law, with the irregularity of leaf density largely ignored for explaining this deviation. Here, we choose 11 bamboo species to explore the allometric relationships among leaf area (*A*), density (ρ), length (*L*), thickness (*T*), and weight (*W*). Because the edge of a bamboo leaf follows a simplified two-parameter Gielis equation, we could show that *A* ∝ *L*^{2} and that *A* ∝ *T*^{2}. This then allowed us to derive the density-thickness allometry ρ ∝ *T*^{b} and the weight-area allometry *W* ∝ *A*^{(b+3)/2} ≈ *A*^{9/8}, where *b* approximates −3/4. Leaf density is strikingly negatively associated with leaf thickness, and it is this inverse relationship that results in the weight-area allometry to deviate from the 3/2-power law. In conclusion, although plants are prone to invest less dry mass and thus produce thinner leaves when the leaf area is sufficient for photosynthesis, such leaf thinning needs to be accompanied with elevated density to ensure structural stability. The findings provide the insights on the evolutionary clue about the biomass investment and output of photosynthetic organs of plants. Because of the importance of leaves, plants could have enhanced the ratio of dry material per unit area of leaf in order to increase the efficiency of photosynthesis, relative the other parts of plants. Although the conclusion is drawn only based on 11 bamboo species, it should also be applicable to the other plants, especially considering previous works on the exponent of the weight-area relationship being less than 3/2 in plants.

## Introduction

The leaf of a plant is an important organ for transpiration, photosynthesis and heat balance (Parkhurst and Loucks, 1972; Goudriaan and Laar, 1994; Wright et al., 2004a,b). Leaf traits such as leaf area, length, width, weight, thickness, density, shape and leaf mass per area (LMA) are thus important functional indices of the leaf economic spectrum (Wright et al., 2004b; Poorter and Rozendaal, 2008). Due to their simplicity, these features have been extensively used by physiologists and ecologists to explain plant growth, developmental rate, water and nutrient use, and plant distributions (Cornelissen et al., 2003; Wright et al., 2004a,b, 2005; Shipley et al., 2006), such as for fruit development (Montero et al., 2000; Iii and Martinson, 2003; Demirsoy et al., 2004, 2005; Fallovo et al., 2008; Buttaro et al., 2015) and crop yield (Stoppani et al., 2003; Salerno et al., 2005; Peksen, 2007). These leaf functional traits are tightly related with each other. For instance, the leaf area, an important canopy parameter that directly affects light interception, light penetration, leaf energy balance and the distribution of solar radiation (Chiariello, 1984; Niinemets and Kull, 1994), can be accurately estimated by a simple equation of leaf dimensions (leaf width and length) (Gamiely et al., 1991; Whitworth et al., 1992; Uzun and Çelik, 1999), especially in crops such as rabbiteye blueberries, squash, strawberry and grapes (NeSmith, 1991, 1992; Montero et al., 2000; Demirsoy et al., 2005). The LMA, which is strongly correlated with many leaf functional, biochemical and structural traits (Poorter et al., 2009; Villar et al., 2013; Puglielli et al., 2015a,b; Millamoreno et al., 2016; John et al., 2017), tends to be small for fast-growing species (Wright et al., 2004a; Griffith et al., 2016) and is tightly related to gross leaf structure (Cornelissen and Thompson, 1997; Cornwell et al., 2008; Onoda et al., 2011).

Leaf dry weight (*W*) and leaf surface area (*A*) are two important traits for the vast majority of vascular plants (Roderick et al., 1999; Sack et al., 2003; Calvo-Alvarado et al., 2008; Pan et al., 2013). The relationship between *A* and *W* follows a power law, *W* = α *A*^{β}, where α is the normalization constant and β the scaling exponent (Niklas et al., 2007; Niklas and Christianson, 2011), suggesting that the light-capturing surface per unit investment of dry mass is an important factor to leaf size.

Bamboos, the subfamily Bambusoideae of the grass family Poaceae (or Graminaceae), have been widely used for food, handicraft and construction materials. It includes 75 genera and 1,300 species, covering 25 million hectares worldwide (Liese and Köhl, 2015), with more than 500 species distributed in China. Besides the apparent differences in leaf size among bamboo species, the bamboo leaf appears to be very similar in shape across species, and its bilateral symmetry can be accurately described by a simplified Gielis equation (Lin et al., 2016). To date, the relationship between leaf area and leaf weight of bamboos have not been systematically examined. The principle of similarity (Thompson, 1917) states that the weight of an organism should follow a 3/2-power law with its surface area or proportional to its volume if the density does not change. However, the leaf weight-area allometry has been reported to greatly deviate from the 3/2-power law, with the reason yet to be clarified. We argue that this could be due to the currently unspecified scaling relationship between leaf thickness and density. Here, using data measured for 11 bamboo species, we first explore the relationship between the leaf thickness and leaf density and then derive the allometric relationship between leaf area and leaf weight, through which providing an overall explanation for the allometric relationships among leaf area, density, length, thickness and weight.

## Materials and Methods

### Sample Collection

We measured leaf surface area and leaf biomass values of at least 100 leaves for each of the 11 bamboo species located in the Nanjing Forestry University campus used in this study (Table 1). Leaf biomass was measured as fresh weight using ME204 (METTLER TOLEDO Equipment Limited Company, Shanghai, China; *d* = 0.0001 g) for species 1–7 in 2016, and an electronic scale with precision 0.01 g (JM-A3002; Chaozeheng Equipment Limited Company, Zhuji, Zhejiang, China) for species 8–11 in 2014. The leaves were weighed as soon as they were picked from the plants, so the loss of water in leaves can be neglected. Species 8–11 belong to *Indocalamus*, and their leaf sizes are all large. Thus, we used a low-precision electronic scale that can directly weigh large leaves.

### Calculations and Statistics

To connect these leaf functional traits, let *W* denote leaf weight, *V* leaf volume, ρ leaf density, *A* leaf surface area, *L* leaf length and *T* leaf (mean) thickness. Leaf shapes of bamboos have been demonstrated to follow a simplified Gielis equation (Gielis, 2003; Shi et al., 2015; Lin et al., 2016):

Here *r* and φ represent the polar radius and angle from the horizontal axis (i.e., polar coordinates), *l* and *n* are constants. Figure 1 illustrates the scanned edge of a bamboo leaf and the prediction from the fitted Gielis equation. Leaf length (*L*) can be described by the following equation (Shi et al., 2015):

Assuming that parameter *n* is a constant for different individuals of a same species, we could estimate the leaf surface area (*A*) as:

(see Appendix for mathematical proofs).

**Figure 1**. Comparison between the scanned and predicted leaf edges of *P. incarnata*. The gray curve represents the scanned leaf edge, and the red solid curve represents the predicted leaf edge by the simplified Gielis equation.

By definition, we have

According to the principle of similarity (Thompson, 1917), we have

The above equation will be tested using the experimental data (see below for details). Let us assume that leaf density and thickness follow the following relationship:

Alternatively, this can be described as:

Here, α, β, *a*, and *b* are constants. For convenience, we can rewrite Equation (7) to

where *p* represents the proportionality. It is apparent that the slope *b* can determine the scaling relationship between leaf weight and leaf area.

If leaf density was not related to leaf thickness, *b* = 0, we would expect to have *W* ∝ *A*^{3/2}, directly following the principle of similarity (Thompson, 1917).

In the aforementioned derivation, the key assumption is that leaf area is proportional to the 2-power of leaf thickness (i.e., Equation 5). In order to examine this hypothesis, we measured the leaf thickness of 100 leaves of *P. incarnata* randomly chosen from 210 leaves that had been used. Because the thickness at different locations on a leaf is different, the mean thicknesses was measured at sampled points that are not located on the main vein (Figure 2). Fresh leaves were punched early in the morning, and actual leaf thickness (LT) measured using hand sectioning. Leaf cross sections were photographed and measured under a light microscopy (DM 2500, a Leica Microsysterms CMS GmbH, Wetziar, Germany) at 20× magnification. We measured the distance between the adaxial and abaxial surfaces of the leaf at the base, middle and tip. Figure 2 illustrates the locations of sampled points in a leaf for measuring the mean thickness. Leaf areas were measured based on scanned images (Shi et al., 2015).

**Figure 2**. An illustration of how the mean thickness of a bamboo leaf is measured. The yellow points represent the locations for measuring the thickness. The data on the main vein were neglected because the values are extremely higher than those apart from the main vein.

The scaling relationship was demonstrated using the reduced major axis (RMA; Milla and Reich, 2007). The reduced major axis is a linear-function fitting method. Different from the linear regression based on the ordinary least squares, the estimate of slope in the RMA is equal to the square root of the variance of response variables divided by the variance of independent variables (Hui et al., 2010). For the RMA and the ordinary least-squares (OLS), the applied target functions for minimizing the residuals to obtain the estimates of parameters are both residual sum of squares (RSS). The bootstrap percentile method (Efron and Tibshirani, 1993; Sandhu et al., 2011) was used to calculate the 95% confidence interval of the slope of the linear model. If a hypothesized power is not included within the 95% confidence interval of the slope (exponent) estimate, this will indicate an allometric relationship with the slope other than the tested value (in this case, 3/2). All analyses were done using R (version 3.2.2) (R Core Team, 2015).

The raw data used in the present study of leaf area, weight and thickness can be found in the Supplementary Material.

## Results

### Proportionality to Leaf Density and Proportionality to Leaf Thickness

Using 1864 data points, we obtained ln(ρ) = −5.00 −0.736 ln(*T*), RSS = 26.7, and *R*^{2} = 0.94 (Figure 3). That is, there is an inverse relationship between leaf thickness and leaf density. The 95% confidence interval of the slope is (−0.744, −0.727) based on 3000 bootstrap replications. Interestingly, the estimate of slope (i.e., *b* in Equation 9) is very approximate −3/4, which frequently appears in the investigation of the self-thinning law and metabolic theoretical ecology.

**Figure 3**. The linear fitting between the natural logarithm of the proportionality to leaf thickness and that of leaf density. Different colors represent different bamboo species, and there are totally 11 species. The proportionality of leaf thickness was obtained by *A*^{1/2}, and the proportionality of leaf density was obtained by *W*·*A*^{−3/2}.

### Leaf Mean Thickness (T) and Leaf Surface Area (A)

Using the leaf area and thickness data of 100 leaves of *P. incarnata*, we obtained ln(*A*) = 12.38 + 2.13 ln(*T*), RSS = 1.07, and *R*^{2} = 0.85 (Figure 4). The 95% confidence interval of slope is (1.98, 2.30) including the theoretical assumption of 2.

**Figure 4**. The linear fitting between the natural logarithm of leaf mean thickness and that of leaf area.

### Leaf Surface Area (A) and Leaf Fresh Weight (W)

Using 1864 data points, we obtained ln(*W*) = −5.03 + 1.147 ln(*A*), RSS = 26.4, and *R*^{2} = 0.99 (Figure 5). The 95% confidence interval of slope is (1.143, 1.152) based on 3000 bootstrap replications. The upper bound of the 95% confidence interval 1.152 is lower than 3/2, so the leaf weight-area allometry does not follow the 3/2-power law.

**Figure 5**. The linear fitting between the natural logarithm of leaf area and that of leaf fresh weight. Different colors represent different bamboo species, and there are totally 11 species.

## Discussion

In this work, a negative relationship between leaf thickness and density was found in 11 bamboo species, consistent with the results for sclerophytes, mesophytes and succulents where the slope ranges from −0.46 to −0.82 (Vendramini et al., 2002). This negative correlation indicates that thinner leaves are also denser (Yano and Terashima, 2004), with more densely packed leaf cells (Pyankov et al., 1999; Vasfilov, 2012), allowing the thinner leaves to catch more light. Of course, leaf thickness and density can be affected by environmental factors, especially light. Leaf thickness tends to increase with decreasing rainfall, humidity, soil fertility and light (Beadle, 1966; Sobrado, 2008). Evergreen leaves are mostly thicker than deciduous leaves (Mooney and Rundel, 1979). Witkowski and Lamont (1991) also demonstrated that leaf area and thickness often vary independently with leaf position within a plant and between species. Here we showed that the correlation between leaf thickness and density was approximate −3/4, similar to that of the herbaceous plants (Vasfilov, 2012).

The possibility of the negative correlation between leaf thickness and density could stabilize LMA. Identifying the contributions of leaf thickness (*T*) and leaf density (ρ) to LMA is important in plant science. For example, Choong et al. (1992) found that the variation in LMA can be primarily driven by the variation of *T*. In contrast, Castro-Díez et al. (2000) confirmed that the LMA was correlated with ρ but not with *T* based on the observations of leaf inner structures of 52 European woody species. Xiong et al. (2016) also demonstrated that leaf density is the major cause of the variation in LMA across rice varieties. Villar et al. (2013) reported that ρ and *T* contributed equally to explaining the variation in LMA.

The correlation between leaf thickness and leaf area has been poorly studied. Niklas et al. (2007, 2009) found the leaf thickness increases with leaf area, but decreases for some species (e.g., *Auranticarpa rhombifolia* (A. Cunn. ex Hook.) L. W. Cayzer et al., *Ginkgo biloba* L., and *Populus nigra* L.). The estimate of slope between the log-transformed data of leaf area and mean leaf thickness in the present study is 2.13 based on 100 leaves of *P. incarnata*, supporting our assumption of a slope of 2 (95% confidence interval from 1.98 to 2.30).

Because the scaling exponent governing the relationship between leaf dry and fresh weight across species is statistically indistinguishable from unity (Figure 6; Niklas et al., 2007; Shi et al., 2015), we directly used fresh weight (*W*_{F}) rather than dry weight (*W*_{D}):

**Figure 6**. Comparison between the observed (points) and predicted fresh weights by using dry weights (the straight line) of *I. barbatus*. 30 leaves were randomly picked from these plants in the Nanjing Forestry University campus on 11 August, 2015. We measured their fresh weights as soon as we picked them from the plants (17:00 or so). Then these leaves were transferred to the oven at 60°C for 48 h. We measured their corresponding dry weights at 18:00 on 13 August, 2015. Then a linear regression was conducted between fresh weight (*W*_{F}) and dry weight (*W*_{D}). The estimate of slope is 2.0830 ± 0.0178, and *R*^{2} = 0.9979. That is, there is a strong proportional relationship between fresh and dry weights.

Price and Enquist (2007) have proposed a model to account for the scaling of leaf area and weight based on the WBE model (West et al., 1999). The overall scaling exponent of leaf area to leaf weight was allometric (β = 1.08 in Milla et al., 2008; β = 1.10 in Milla and Reich, 2007; β = 1.02 in Niklas et al., 2007; β is from 0.859 to 1.299 in Pan et al., 2013). Here, we obtained ln (*W*) = −5.03 + 1.147 ln(*A*), RSS = 26.1, and *R*^{2} = 0.99 using 1864 data points. The 95% confidence interval of slope is (1.143, 1.152) based on 3000 bootstrap replications. Our scaling exponent 1.147 is in accordance with those reported in other studies (Niklas et al., 2007). Milla and Reich (2007) discovered the average exponent (1.10 with 95% CI 1.08–1.13) was significantly greater than 1 but lower than 3/2, and found that within each species large leaves deploy less surface area per unit dry mass than small ones. Niklas et al. (2007) also found that large leaves deploy less light-absorbing area per unit dry mass investment than small leaves in 943 species. Leaves are subject to strong selection pressure along gradients of aridity, solar radiation and nutrient availability that affect their size and shape (Price and Enquist, 2007). This leads to a wide range of variation in leaf area to leaf weight scaling exponent among species (Milla and Reich, 2007). The scaling relationship reflects the outcome of an evolutionary trade-off among ancestral metabolic, morphological and anatomical traits shared by all vascular plants (Niklas et al., 2007).

For bamboos, leaf area was demonstrated to be proportional to leaf mean thickness squared and also to be proportional to leaf length squared. Consequently, leaf volume can be expressed by leaf area to the power 3/2, which means that leaf volume, leaf area, leaf mean thickness (or leaf length) follow the similarity of principle of Thompson. The allometric relationship between leaf weight and leaf area based on the experimental data of 11 bamboo species was demonstrated to significantly deviate from the 3/2–power law. In other words, leaf weight is not proportional to leaf volume because leaf volume is the proportionality of leaf area to the power 3/2. To find the reason of the deviation to the 3/2–power law, we analyzed the relationship between the proportionality of leaf density and that of leaf mean thickness, and found a significant negative linear relationship between the log-transformed data of these two variables. That is, leaf density decreases with leaf size (represented by leaf mean thickness or leaf length or leaf area) increasing, which answers why leaf weight does not have a proportional relationship with leaf area to the power 3/2. The main conclusion in the present study that the negative scaling exponent of leaf thickness to density results in the deviation of the scaling exponent of leaf area to weight from the 3/2–power law could be potentially extended to other plants. However, the negative scaling exponent of leaf thickness to density (as demonstrated to be approximate to −3/4 for bamboos) might vary with different species of plants. In addition, whether the leaf area-thickness allometry follows the 2–power law for other plants also merits further investigation.

## Author Contributions

SL and PS designed the experiment; LS, FL, YS, QW, and YD carried out this experiment; PS analyzed the data; PS, SL, CH, JG, and GR wrote the manuscript. All authors read and commented on this manuscript.

## Conflict of Interest Statement

The reviewer Z-XF declared a shared affiliation, with no collaboration, with one of the authors YS, to the handling Editor.

The other 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.

## Acknowledgments

We are thankful to two reviewers and the handling editor for their constructive and valuable comments on the earlier version of this paper. This work was financially supported by the National Key Research & Development Program of China (2016YFD0600901); Jiangsu Province Support Project (LYSX [2016]04; BE2016304); the National Natural Science Foundation for Young Scholars of China (31000294) and the Priority Academic Program Development of Jiangsu Higher Education Institutions.

## Supplementary Material

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

## References

Beadle, N. C. W. (1966). Soil phosphate and its role in molding segments of the Australian flora and vegetation, with special reference to xeromorphy and sclerophylly. *Ecology* 47, 992–1007. doi: 10.2307/1935647

Buttaro, D., Rouphael, Y., Rivera, C. M., Colla, G., and Gonnella, M. (2015). Simple and accurate allometric model for leaf area estimation in *Vitis vinifera* L. *genotypes. Photosynthetica* 53, 1–8. doi: 10.1007/s11099-015-0117-2

Calvo-Alvarado, J. C., McDowell, N. G., and Waring, R. H. (2008). Allometric relationships predicting foliar biomass and leaf area: sapwood area ratio from tree height in five Costa Rican rain forest species. *Tree Physiol.* 28, 1601–1608. doi: 10.1093/treephys/28.11.1601

Castro-Díez, P., Puyravaud, J. P., and Cornelissen, J. H. C. (2000). Leaf structure and anatomy as related to leaf mass per area variation in seedlings of a wide range of woody plant species and types. *Oecologia* 124, 476–486. doi: 10.1007/PL00008873

Chiariello, N. (1984). “Leaf energy balance in the wet lowland tropics,” in *Physiological Ecology of Plants of the Wet Tropics, Tasks for Vegetation Science, 12*, eds E. Medina, H. A. Mooney, C. VáZquez-Yánes (Dordrecht: Springer), 85–98.

Choong, M. F., Lucas, P. W., Ong, J. S. Y., Pereira, B., Tan, H. T. W., and Turner, I. M. (1992). Leaf fracture toughness and sclerophylly: their correlations and ecological implications. *New Phytol.* 121, 597–610. doi: 10.1111/j.1469-8137.1992.tb01131.x

Cornelissen, J. H. C., Lavorel, S., Garnier, E., Díaz, S., Buchmann, N., Gurvich, D. E., et al. (2003). A handbook of protocols for standardised and easy measurement of plant functional traits worldwide. *Aust. J. Bot.* 51, 335–380. doi: 10.1071/BT02124

Cornelissen, J. H. C., and Thompson, K. (1997). Functional leaf attributes predict litter decomposition rate in herbaceous plants. *New Phytol.* 135, 109–114. doi: 10.1046/j.1469-8137.1997.00628.x

Cornwell, W. K., Cornelissen, J. H. C., Amatangelo, K., Dorrepaal, E., Eviner, V. T., Godoy, O., et al. (2008). Plant species traits are the predominant control on litter decomposition rates within biomes worldwide. *Ecol. Lett.* 11, 1065–1071. doi: 10.1111/j.1461-0248.2008.01219.x

Demirsoy, H., Demirsoy, L., and Öztürk, A. (2005). Improved model for the non-destructive estimation of strawberry leaf area. *Fruits* 60, 69–73. doi: 10.1051/fruits:2005014

Demirsoy, H., Demirsoy, L., Uzun, S., and Ersoy, B. (2004). Non-destructive leaf area estimation in peach. *Europ. J. Hort. Sci.* 69, 144–146.

Efron, B., and Tibshirani, R. J. (1993). *An Introduction to the Bootstrap, Monogragh on Statistics and Applied Probability*. London: Chapman and Hall.

Fallovo, C., Cristofori, V., De-Gyves, E. M., Rivera, C. M., Rea, R., and Fanasca, S. (2008). Leaf area estimation model for small fruits from linear measurements. *HortScience* 43, 2263–2267.

Gamiely, S., Randle, W. M., Mills, H. A., and Smittle, D. A. (1991). A rapid and non-destructive method for estimating leaf area of onions. *HortScience* 26:206.

Gielis, J. (2003). A generic geometric transformation that unifies a wide range of natural and abstract shapes. *Am. J. Bot.* 90, 333–338. doi: 10.3732/ajb.90.3.333

Goudriaan, J., and Laar, H. H. V. (1994). *Modelling Potential Crop Growth Processes*. Wageningen: Springer Current Issues in Production Ecology.

Griffith, D. M., Quigley, K. M., and Anderson, T. M. (2016). Leaf thickness controls variation in leaf mass per area (LMA) among grazing-adapted grasses in Serengeti. *Physiol. Ecol.* 181, 1035–1040. doi: 10.1007/s00442-016-3632-3

Hui, C., Terblanche, J. S., Chown, S. L., and McGeoch, M. A. (2010). Parameter landscapes unveil the bias in allometric prediction. *Methods Ecol. Evol.* 1, 69–74. doi: 10.1111/j.2041-210X.2009.00005.x

Iii, L. W., and Martinson, T. E. (2003). Nondestructive leaf area estimation of ‘Niagara' and ‘DeChaunac' grapevines. *Sci. Hort.* 98, 493–498. doi: 10.1016/S0304-4238(03)00020-7

John, G. P., Scoffoni, C., Buckley, T. N., Villar, R., Poorter, H., and Sack, L. (2017). The anatomical and compositional basis of leaf mass per area. *Ecol. Lett.* 20, 412–425. doi: 10.1111/ele.12739

Lin, S. Y., Zhang, L., Reddy, G. V. P., Hui, C., Gielis, J., Ding, Y., et al. (2016). A geometrical model for testing bilateral symmetry of bamboo leaf with a simplified Gielis equation. *Ecol. Evol.* 6, 6798–6806. doi: 10.1002/ece3.2407

Milla, R., and Reich, P. B. (2007). The Scaling of Leaf Area and Mass: the cost of light interception increases with leaf size. *Proc. R. Soc. B* 274, 2109–2114. doi: 10.1098/rspb.2007.0417

Milla, R., Reich, P. B., Niinemets, Ü., and Castro-Díez, P. (2008). Environmental and developmental controls on specific leaf area are little modified by leaf allometry. *Funct. Ecol.* 22, 565–576. doi: 10.1111/j.1365-2435.2008.01406.x

Millamoreno, E. A., McKown, A. D., Guy, R. D., and Soolanayakanahally, R. Y. (2016). Leaf mass per area predicts palisade structural properties linked to mesophyll conductance in balsam poplar (*Populus balsamifera* L.). *Botany* 94, 225-239. doi: 10.1139/cjb-2015-0219

Montero, F. J., Juan, J. A. D., Cuesta, A., and Brasa, A. (2000). Nondestructive methods to estimate leaf area in *Vitis vinifera* L. *HortScience* 35, 696–698.

Mooney, H. A., and Rundel, P. W. (1979). Nutrient relations of the evergreen shrub, *Adenostoma fasciculatum*, in the Califonia chaparral. *Bot. Gaz.* 140, 109–113. doi: 10.1086/337064

NeSmith, D. S. (1991). Non-destructive leaf area estimation of rabbit eye blueberries. *HortScience* 26:1332.

Niinemets, Ü., and Kull, K. (1994). Leaf weight per area and leaf size of 85 Estonian woody species in relation to shade tolerance and light availability. *Forest Ecol. Manag.* 70, 1–10. doi: 10.1016/0378-1127(94)90070-1

Niklas, K. J., and Christianson, M. L. (2011). Differences in the scaling of area and mass of *Ginkgo biloba* (Ginkgoaceae) leaves and their relevance to the study of specific leaf area. *Am. J. Bot.* 98, 1381–1386. doi: 10.3732/ajb.1100106

Niklas, K. J., Cobb, E. D., Niinemets, Ü., Reich, P. B., Sellin, A., Shipley, B., et al. (2007). “Diminishing returns” in the scaling of functional leaf traits across and within species groups. *Proc. Natl Acad. Sci. U.S.A.* 104, 8891–8896. doi: 10.1073/pnas.0701135104

Niklas, K. J., Cobb, E. D., and Spatz, H. C. (2009). Predicting the allometry of surface area and dry mass. *Am. J. Bot.* 96, 531–536. doi: 10.3732/ajb.0800250

Onoda, Y., Westoby, M., Adler, P. B., Choong, A. M., Clissold, F. J., Cornelissen, J. H., et al. (2011). Global patterns of leaf mechanical properties. *Ecol. Lett.* 14, 301–312. doi: 10.1111/j.1461-0248.2010.01582.x

Pan, S., Liu, C., Zhang, W., Xu, S., Wang, N., Li, Y., et al. (2013). The scaling relationships between leaf mass and leaf area of vascular plant species change with altitude. *PLoS ONE* 8:e76872. doi: 10.1371/journal.pone.0076872

Parkhurst, D. F., and Loucks, O. L. (1972). Optimal leaf size in relation to environment. *J. Ecol.* 60, 505–537. doi: 10.2307/2258359

Peksen, E. (2007). Non-destructive leaf area estimation model for faba bean (*Vicia faba* L.). *Sci. Hort.* 113, 322–328. doi: 10.1016/j.scienta.2007.04.003

Poorter, H., Niinemets, Ü., Poorter, L., Wright, I. J., and Villar, R. (2009). 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

Poorter, L., and Rozendaal, D. M. A. (2008). Leaf size and leaf display of thirty-eight tropical tree species. *Oecologia* 158, 35–46. doi: 10.1007/s00442-008-1131-x

Price, C. A., and Enquist, B. J. (2007). Scaling mass and morphology in leaves: an extension of the WBE model. *Ecology* 88, 1132–1141. doi: 10.1890/06-1158

Puglielli, G., Crescente, M. F., Frattaroli, A. R., and Gratani, L. (2015a). Leaf mass per area (LMA) as a possible predictor of adaptive strategies in two species of *Sesleria* (Poaceae): analysis of morphological, anatomical and physiological leaf traits. *Ann. Bot. Fenni.* 52, 135–143. doi: 10.5735/085.052.0201

Puglielli, G., Crescente, M. F., Frattaroli, A. R., and Gratani, L. (2015b). Morphological, anatomical and physiological leaf trait plasticity of *Sesleria nitida* (poaceae) in open vs. shaded conditions. *Pol. J. Ecol.* 63, 10–22. doi: 10.3161/15052249PJE2015.63.1.002

Pyankov, V. I., Kondratchuk, A. V., and Shipley, B. (1999). Leaf structure and specific leaf mass: the alpine desert plants of the Eastern Pamirs, Tadjikistan, *New Phytol.* 143, 131–142. doi: 10.1046/j.1469-8137.1999.00435.x

R Core Team (2015). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing. Available Online at: https://www.R-project.org/.

Roderick, M. L., Berry, S. L., Noble, I. R., and Farquhar, G. D. (1999). A theoretical approach to linking the composition and morphology with the function of leaves. *Funct. Ecol.* 13, 683–695. doi: 10.1046/j.1365-2435.1999.00368.x

Sack, L., Cowan, P. D., Jaikumar, N., and Holbrook, N. M. (2003). The “Hydrology” of leaves: coordination of structure and function in temperate woody species. *Plant Cell Environ.* 26, 1343–1356. doi: 10.1046/j.0016-8025.2003.01058.x

Salerno, A., Rivera, C. M., Rouphael, Y., Colla, G., Cardarelli, M., Pierandrei, F., et al. (2005). Leaf area estimation of radish from linear measurements. *Adv. Hort. Sci.* 19, 213–215.

Sandhu, H. S., Shi, P., Kuang, X., Xue, F., and Ge, F. (2011). Applications of the bootstrap to insect physiology. *Fla. Entomol.* 94, 1036–1041. doi: 10.1653/024.094.0442

Shi, P. J., Xu, Q., Sandhu, H. S., Gielis, J., Ding, Y. L., Li, H. R., et al. (2015). Comparison of dwarf bamboos (*Indocalamus sp*.) leaf parameters to determine relationship between spatial density of plants and total leaf area per plant. *Ecol. Evol.* 5, 4578–4589. doi: 10.1002/ece3.1728

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

Sobrado, M. A. (2008). Leaf characteristics and diurnal variation of chlorophyll fluorescence in leaves of the ‘bana' vegetation of the amazon region. *Photosynthetica* 46:202. doi: 10.1007/s11099-008-0033-9

Stoppani, M. I., Wolf, R., Francescangeli, N., and Marti, H. R. (2003). A nondestructive and rapid method for estimating leaf area of broccoli. *Adv. Hort. Sci.* 17, 173–175.

Uzun, S., and Çelik, H. (1999). Leaf area prediction models for different horticultural plants, *Turk. J. Agric. For.* 23, 645–650.

Vasfilov, S. P. (2012). Analysis of the causes of variability of the dry leaf mass-per-area ratio. *Biol. Bull. Rev.* 3, 238–253. doi: 10.1134/S2079086412030061

Vendramini, F., Diaz, S., Gurvich, D. E., Wilson, P. J., Thompson, K., and Hodgson, J. G. (2002). Leaf traits as indicators of resource-use strategy in floras with succulent species. *New Phytol.* 154, 147–157. doi: 10.1046/j.1469-8137.2002.00357.x

Villar, R., Ruiz-Robleto, J., Ubera, J. L., and Poorter, H. (2013). Exploring variation in leaf mass per area (LMA) from leaf to cell: an anatomical analysis of 26 woody species. *Am. J. Bot.* 100, 1969–1980. doi: 10.3732/ajb.1200562

West, G. B., Brown, J. H., and Enquist, B. J. (1999). A general model for the structure and allometry of plant vascular systems. *Nature* 400, 664–667. doi: 10.1038/23251

Whitworth, J. L., Mauromoustakos, A., and Smith, M. W. (1992). A non-destructive method for estimation of leaf area in pecan, *HortScience* 27:851.

Witkowski, E. T. F., and Lamont, B. B. (1991). Leaf specific mass confounds leaf density and thickness. *Oecologia* 88, 486–493. doi: 10.1007/BF00317710

Wright, I. J., Groom, P. K., Lamont, B. B., Poot, P., Prior, L. D., Reich, P. B., et al. (2004a). Short communication: leaf trait relationships in Australian plant species. *Funct. Plant Biol.* 31, 551–558. doi: 10.1071/FP03212

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

Wright, I. J., Reich, P. B., Cornelissen, J. H. C., Falster, D. S., Groom, P. K., Hikosaka, K., et al. (2005). Modulation of leaf economic traits and trait relationships by climate. *Glob. Ecol. Biogeogr.* 14, 411–421. doi: 10.1111/j.1466-822x.2005.00172.x

Xiong, D., Wang, D., Liu, X., Peng, S., Huang, J., and Li, Y. (2016). Leaf density explains variation in leaf mass per area in rice between cultivars and nitrogen treatments. *Ann. Bot.* 117, 963–971. doi: 10.1093/aob/mcw022

Yano, S., and Terashima, I. (2004). Developmental process of sun and shade leaves in *Chenopodium album* L. *Plant Cell Environ.* 27, 781–793. doi: 10.1111/j.1365-3040.2004.01182.x

## Appendix

### The Relationship Between Leaf Surface Area (*A*) and Leaf Length (*L*)

where “Gamma” represents the gamma function, and “HypergeometricPFQ” represents the generalized hypergeometric function. The detailed definitions for these two functions can be found in Wolfram MathWorld (http://mathworld.wolfram.com/). Because *n* is a constant, we have:

Leaf length has demonstrated to be proportional to parameter *l* (Shi et al., 2015)

Then we have:

Keywords: allometry, Gielis equation, natural logarithm, proportionality, reduced major axis

Citation: Lin S, Shao L, Hui C, Song Y, Reddy GVP, Gielis J, Li F, Ding Y, Wei Q and Shi P (2018) Why Does Not the Leaf Weight-Area Allometry of Bamboos Follow the 3/2-Power Law? *Front. Plant Sci*. 9:583. doi: 10.3389/fpls.2018.00583

Received: 02 November 2017; Accepted: 13 April 2018;

Published: 04 May 2018.

Edited by:

Jian-Guo Huang, South China Botanical Garden (CAS), ChinaReviewed by:

Ze-Xin Fan, Xishuangbanna Tropical Botanical Garden (CAS), ChinaChiara Cirillo, University of Naples Federico II, Italy

Roland Valcke, University of Hasselt, Belgium

Copyright © 2018 Lin, Shao, Hui, Song, Reddy, Gielis, Li, Ding, Wei and Shi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Peijian Shi, peijianshi@gmail.com