Acessibilidade / Reportar erro

Regression models for estimating leaf area of seedlings and adult individuals of Neotropical rainforest tree species

Modelos de regressão para estimativa da área foliar de indivíduos jovens e adultos de espécies arbóreas de floresta neotropical

Abstract

Individual leaf area (LA) is a key variable in studies of tree ecophysiology because it directly influences light interception, photosynthesis and evapotranspiration of adult trees and seedlings. We analyzed the leaf dimensions (length – L and width – W) of seedlings and adults of seven Neotropical rainforest tree species (Brosimum rubescens, Manilkara maxima, Pouteria caimito, Pouteria torta, Psidium cattleyanum, Symphonia globulifera and Tabebuia stenocalyx) with the objective to test the feasibility of single regression models to estimate LA of both adults and seedlings. In southern Bahia, Brazil, a first set of data was collected between March and October 2012. From the seven species analyzed, only two (P. cattleyanum and T. stenocalyx) had very similar relationships between LW and LA in both ontogenetic stages. For these two species, a second set of data was collected in August 2014, in order to validate the single models encompassing adult and seedlings. Our results show the possibility of development of models for predicting individual leaf area encompassing different ontogenetic stages for tropical tree species. The development of these models was more dependent on the species than the differences in leaf size between seedlings and adults.

Keywords:
leaf allometry; leaf size; leaf shape; ontogenetic stage; tree growth

Resumo

Área foliar individual (AF) é uma variável chave em estudos sobre a ecofisiologia de arbóreas, porque influencia diretamente a interceptação de luz, a fotossíntese e a evapotranspiração das árvores adultas e das mudas. Foram analisadas as dimensões foliares (comprimento - C e largura - L) de indivíduos adultos e de mudas de sete espécies arbóreas de florestas neotropicais (Brosimum rubescens, Manilkara maxima, Pouteria caimito, Pouteria torta, Psidium cattleyanum, Symphonia globulifera e Tabebuia stenocalyx), com o objetivo de testar a viabilidade de modelos de regressão linear para estimar a AF de indivíduos adultos e mudas. No sul da Bahia, Brasil, um primeiro conjunto de dados foi coletado entre março e outubro de 2012. A partir das sete espécies analisadas, apenas duas (P. cattleyanum e T. stenocalyx) apresentaram relações muito semelhantes entre e AF e CL, em ambos os estádios ontogenéticos. Para estas duas espécies, um segundo conjunto de dados foi coletado em agosto de 2014, a fim de validar os modelos únicos que englobam folhas de indivíduos adultos e mudas. Nossos resultados mostram a possibilidade de desenvolvimento de modelos para a predição da área foliar, abrangendo diferentes estádios ontogenéticos para espécies arbóreas tropicais. O desenvolvimento destes modelos foi mais dependente das espécies do que das diferenças entre o tamanho das folhas de mudas e de indivíduos adultos.

Palavras-chave:
alometria foliar; tamanho da folha; formato da folha; estádio ontogenético; crescimento arbóreo

1 Introduction

Individual leaf area (LA) is a key variable in forest ecophysiological studies, because it directly influences light interception, photosynthesis and evapotranspiration of adult trees, saplings and seedlings (Niinemets, 2010Niinemets, Ü., 2010. A review of light interception in plant stands from leaf to canopy in different plant functional types and in species with varying shade tolerance. Ecological Research, vol. 25, no. 4, pp. 693-714. http://dx.doi.org/10.1007/s11284-010-0712-4.
http://dx.doi.org/10.1007/s11284-010-071...
; Zhang and Pan, 2011Zhang, L. and Pan, L., 2011. Allometric models for leaf area estimation across different leaf-age groups of evergreen broadleaved trees in a subtropical forest. Photosynthetica, vol. 49, no. 2, pp. 219-226. http://dx.doi.org/10.1007/s11099-011-0027-x.
http://dx.doi.org/10.1007/s11099-011-002...
). Furthermore, LA is an important variable for seedling establishment and for sun/shade acclimation of tropical rainforest tree species (Varone and Gratani, 2009Varone, L. and Gratani, L., 2009. Leaf expansion in Rhamnus alaternus L. by leaf morphological, anatomical and physiological analysis. Trees (Berlin), vol. 23, no. 6, pp. 1255-1262. http://dx.doi.org/10.1007/s00468-009-0365-5.
http://dx.doi.org/10.1007/s00468-009-036...
; Coste et al., 2011Coste, S., Roggy, J.C., Schimann, H., Epron, D. and Dreyer, E., 2011. A cost–benefit analysis of acclimation to low irradiance in tropical rainforest tree seedlings: leaf life span and payback time for leaf deployment. Journal of Experimental Botany, vol. 62, no. 11, pp. 3941-3955. http://dx.doi.org/10.1093/jxb/err092. PMid:21511904.
http://dx.doi.org/10.1093/jxb/err092...
; Houter and Pons, 2012Houter, N. and Pons, T.L., 2012. Ontogenetic changes in leaf traits of tropical rainforest trees differing in juvenile light requirement. Oecologia, vol. 169, no. 1, pp. 33-45. http://dx.doi.org/10.1007/s00442-011-2175-x. PMid:22038060.
http://dx.doi.org/10.1007/s00442-011-217...
; Laurans et al., 2012Laurans, M., Martin, O., Nicolini, E. and Vincent, G., 2012. Functional traits and their plasticity predict tropical trees regeneration niche even among species with intermediate light requirements. Journal of Ecology, vol. 100, no. 6, pp. 1440-1452. http://dx.doi.org/10.1111/j.1365-2745.2012.02007.x.
http://dx.doi.org/10.1111/j.1365-2745.20...
). The direct determination or estimation of LA is of great importance to evaluate the performance of seedlings in tree plantations and in natural forests (Coste et al., 2011Coste, S., Roggy, J.C., Schimann, H., Epron, D. and Dreyer, E., 2011. A cost–benefit analysis of acclimation to low irradiance in tropical rainforest tree seedlings: leaf life span and payback time for leaf deployment. Journal of Experimental Botany, vol. 62, no. 11, pp. 3941-3955. http://dx.doi.org/10.1093/jxb/err092. PMid:21511904.
http://dx.doi.org/10.1093/jxb/err092...
; Houter and Pons, 2012Houter, N. and Pons, T.L., 2012. Ontogenetic changes in leaf traits of tropical rainforest trees differing in juvenile light requirement. Oecologia, vol. 169, no. 1, pp. 33-45. http://dx.doi.org/10.1007/s00442-011-2175-x. PMid:22038060.
http://dx.doi.org/10.1007/s00442-011-217...
; Laurans et al., 2012Laurans, M., Martin, O., Nicolini, E. and Vincent, G., 2012. Functional traits and their plasticity predict tropical trees regeneration niche even among species with intermediate light requirements. Journal of Ecology, vol. 100, no. 6, pp. 1440-1452. http://dx.doi.org/10.1111/j.1365-2745.2012.02007.x.
http://dx.doi.org/10.1111/j.1365-2745.20...
). The methods used for LA estimation can be either destructive or nondestructive. Destructive methods involve leaf harvesting, which can be done using automated leaf area meters, photocopies, planimeters (Tsuda, 1999Tsuda, M., 1999. Errors in leaf area measurement with an automatic area meter due to leaf chlorophyll in crop plants. Annals of Botany, vol. 84, no. 6, pp. 799-801. http://dx.doi.org/10.1006/anbo.1999.0984.
http://dx.doi.org/10.1006/anbo.1999.0984...
) and softwares that calculate LA based on scanned images and digital photos (O’Neal et al., 2002O’Neal, M.E., Landis, D.A. and Isaacs, R., 2002. An inexpensive, accurate method for measuring leaf area and defoliation through digital image analysis. Journal of Economic Entomology, vol. 95, no. 6, pp. 1190-1194. http://dx.doi.org/10.1603/0022-0493-95.6.1190. PMid:12539831.
http://dx.doi.org/10.1603/0022-0493-95.6...
; Godoy et al., 2007Godoy, L.J.G., Yanagiwara, R.S., Bôas, R.L.V., Backes, C. and Lima, C.P., 2007. Análise da imagem digital para estimativa da área foliar em plantas de laranja “Pera”. Revista Brasileira de Fruticultura, vol. 29, no. 3, pp. 420-424. http://dx.doi.org/10.1590/S0100-29452007000300004.
http://dx.doi.org/10.1590/S0100-29452007...
; Toebe et al., 2010Toebe, M., Brum, B., Lopes, S.J., Cargnelutti-Filho, A. and Silveira, T.R., 2010. Estimativa da área foliar de Crambe abyssinica por discos foliares e por fotos digitais. Ciência Rural, vol. 40, no. 2, pp. 475-478. http://dx.doi.org/10.1590/S0103-84782010000200036.
http://dx.doi.org/10.1590/S0103-84782010...
). Nondestructive methods do not require leaf destruction and are based on regression models to estimate LA from leaf linear dimensions, as length (L) and width (W) (Kandiannan et al., 2002Kandiannan, K., Kailasam, C., Chandaragiri, K.K. and Sankaran, N., 2002. Allometric model for leaf area estimation in black pepper (Piper nigrum L.). Journal Agronomy & Crop Science, vol. 188, no. 2, pp. 138-140. http://dx.doi.org/10.1046/j.1439-037X.2002.00540.x.
http://dx.doi.org/10.1046/j.1439-037X.20...
; Antunes et al., 2008Antunes, W.C., Pompelli, M.F., Carretero, D.M. and Damatta, F.M., 2008. Allometric models for non-destructive leaf area estimation in coffe (Coffea arabica and Coffea canephora). Annals of Applied Biology, vol. 153, no. 1, pp. 33-40. http://dx.doi.org/10.1111/j.1744-7348.2008.00235.x.
http://dx.doi.org/10.1111/j.1744-7348.20...
; Rouphael et al., 2010Rouphael, Y., Mouneimne, A.H., Ismail, A., Mendoza-Degyves, E., Rivera, C.M. and Colla, G., 2010. Modeling individual leaf area of rose (Rosa hybrida L.) based on leaf length and width measurement. Photosynthetica, vol. 48, no. 1, pp. 9-15. http://dx.doi.org/10.1007/s11099-010-0003-x.
http://dx.doi.org/10.1007/s11099-010-000...
; Zhang and Pan, 2011Zhang, L. and Pan, L., 2011. Allometric models for leaf area estimation across different leaf-age groups of evergreen broadleaved trees in a subtropical forest. Photosynthetica, vol. 49, no. 2, pp. 219-226. http://dx.doi.org/10.1007/s11099-011-0027-x.
http://dx.doi.org/10.1007/s11099-011-002...
; Pompelli et al., 2012Pompelli, M.F., Antunes, W.C., Ferreira, D.T.R.G., Cavalcante, P.G.S., Wanderley-Filho, H.C.L. and Endres, L., 2012. Allometric models for non-destructive leaf area estimation of Jatropha curcas. Biomass and Bioenergy, vol. 36, pp. 77-85. http://dx.doi.org/10.1016/j.biombioe.2011.10.010.
http://dx.doi.org/10.1016/j.biombioe.201...
; Souza and Habermann, 2014Souza, M.C. and Habermann, G., 2014. Non-destructive equations to estimate the leaf area of Styrax pohliiand and Styrax ferrugineus. Brazilian Journal of Biology = Revista Brasileira de Biologia, vol. 74, no. 1, pp. 222-225. http://dx.doi.org/10.1590/1519-6984.17012. PMid:25055106.
http://dx.doi.org/10.1590/1519-6984.1701...
).

Depending on the objective of the study, destructive methods can be more advantageous than nondestructive ones. For example, from a practical point of view, destructive methods that use leaf x area meters are the fastest among all, enabling assessment of a large number of leaves within a short time (Tsuda, 1999Tsuda, M., 1999. Errors in leaf area measurement with an automatic area meter due to leaf chlorophyll in crop plants. Annals of Botany, vol. 84, no. 6, pp. 799-801. http://dx.doi.org/10.1006/anbo.1999.0984.
http://dx.doi.org/10.1006/anbo.1999.0984...
). In addition, scanned images and digital photographs allow the construction of image banks, which can then be used, for example, for studies aiming to analyze the color and shape of the leaves. On the other hand, the development of regression models to estimate LA based on easily measured variables is a practical and fast methodology for using at the field. The regression models make possible to obtain reliable estimates of individual leaf areas using simple tools, as ruler and calculator. This implies the possibility to analyze the development of a single leaf and prevents loss of photosynthetic leaf surface and the resulting energy cost used in the formation of plant tissues.

Regression models encompassing seedlings and adults for estimate LA may be useful in field studies of ecology and ecophysiology of tropical tree species, such as studies focused on differences in the herbivory rate between seedlings and adults (Boege and Marquis, 2005Boege, K. and Marquis, R.J., 2005. Facing herbivory as you grow up: the ontogeny of resistance in plants. Trends in Ecology & Evolution, vol. 20, no. 8, pp. 441-448. http://dx.doi.org/10.1016/j.tree.2005.05.001. PMid:16701415.
http://dx.doi.org/10.1016/j.tree.2005.05...
) and differences in the leaf development and photosynthetic behavior between seedlings and adults (Ishida et al., 2005Ishida, A., Yazaki, K. and Hoe, A.L., 2005. Ontogenetic transition of leaf physiology and anatomy from seedlings to mature trees of a rain forest pioneer tree, Macaranga gigantea. Tree Physiology, vol. 25, no. 5, pp. 513-522. http://dx.doi.org/10.1093/treephys/25.5.513. PMid:15741149.
http://dx.doi.org/10.1093/treephys/25.5....
). It is especially important in highly biodiverse forests, as the Neotropical rainforests (Thomas, 1999Thomas, W.W., 1999. Conservation and monographic research on the flora of tropical America. Biodiversity and Conservation, vol. 8, no. 8, pp. 1007-1015. http://dx.doi.org/10.1023/A:1008857429787.
http://dx.doi.org/10.1023/A:100885742978...
), where the large number of species requires the development of numerous models for the prediction of species-specific individual leaf areas. However, the development of a single regression model to estimate LA of seedlings and adults may not be appropriate because of the biological characteristics of tree species. First, seedlings of the same tree species have inherently smaller individual leaf areas than adults, because they have smaller size (Ishida et al., 2005Ishida, A., Yazaki, K. and Hoe, A.L., 2005. Ontogenetic transition of leaf physiology and anatomy from seedlings to mature trees of a rain forest pioneer tree, Macaranga gigantea. Tree Physiology, vol. 25, no. 5, pp. 513-522. http://dx.doi.org/10.1093/treephys/25.5.513. PMid:15741149.
http://dx.doi.org/10.1093/treephys/25.5....
). Second, many species have differences in leaf morphological traits between leaf seedlings, saplings and adults, which is related to the regeneration niche (Thomas and Ickes, 1995Thomas, S.C. and Ickes, K., 1995. Ontogenetic changes in leaf size in Malaysian Rain Forest Trees. Biotropica, vol. 27, no. 4, pp. 427-434. http://dx.doi.org/10.2307/2388954.
http://dx.doi.org/10.2307/2388954...
; Ishida et al., 2005Ishida, A., Yazaki, K. and Hoe, A.L., 2005. Ontogenetic transition of leaf physiology and anatomy from seedlings to mature trees of a rain forest pioneer tree, Macaranga gigantea. Tree Physiology, vol. 25, no. 5, pp. 513-522. http://dx.doi.org/10.1093/treephys/25.5.513. PMid:15741149.
http://dx.doi.org/10.1093/treephys/25.5....
; Houter and Pons, 2012Houter, N. and Pons, T.L., 2012. Ontogenetic changes in leaf traits of tropical rainforest trees differing in juvenile light requirement. Oecologia, vol. 169, no. 1, pp. 33-45. http://dx.doi.org/10.1007/s00442-011-2175-x. PMid:22038060.
http://dx.doi.org/10.1007/s00442-011-217...
). Even though there are many references of regression models in the literature for predicting of LA in tropical tree species (Antunes et al., 2008Antunes, W.C., Pompelli, M.F., Carretero, D.M. and Damatta, F.M., 2008. Allometric models for non-destructive leaf area estimation in coffe (Coffea arabica and Coffea canephora). Annals of Applied Biology, vol. 153, no. 1, pp. 33-40. http://dx.doi.org/10.1111/j.1744-7348.2008.00235.x.
http://dx.doi.org/10.1111/j.1744-7348.20...
; Zhang and Pan, 2011Zhang, L. and Pan, L., 2011. Allometric models for leaf area estimation across different leaf-age groups of evergreen broadleaved trees in a subtropical forest. Photosynthetica, vol. 49, no. 2, pp. 219-226. http://dx.doi.org/10.1007/s11099-011-0027-x.
http://dx.doi.org/10.1007/s11099-011-002...
; Pompelli et al., 2012Pompelli, M.F., Antunes, W.C., Ferreira, D.T.R.G., Cavalcante, P.G.S., Wanderley-Filho, H.C.L. and Endres, L., 2012. Allometric models for non-destructive leaf area estimation of Jatropha curcas. Biomass and Bioenergy, vol. 36, pp. 77-85. http://dx.doi.org/10.1016/j.biombioe.2011.10.010.
http://dx.doi.org/10.1016/j.biombioe.201...
; Souza and Habermann, 2014Souza, M.C. and Habermann, G., 2014. Non-destructive equations to estimate the leaf area of Styrax pohliiand and Styrax ferrugineus. Brazilian Journal of Biology = Revista Brasileira de Biologia, vol. 74, no. 1, pp. 222-225. http://dx.doi.org/10.1590/1519-6984.17012. PMid:25055106.
http://dx.doi.org/10.1590/1519-6984.1701...
), we are aware of no published studies aiming to compare regression models developed for both seedlings and adult individuals of tropical trees.

In this study, we analyzed leaf dimensions of seedlings and adults of seven Neotropical rainforest tree species with the objective to test the feasibility of single regression models to estimate LA of both adults and seedlings. We hypothesize that integrative predictive models for both ontogenetic stages provide similar adjustments compared to single models for each ontogenetic stage despite differences in leaf size between seedlings and adults.

2 Material and Methods

We analyzed the leaf dimensions of seedlings and adults from seven native species of southern Bahia, Brazil. Species were chosen based on their morphological characteristics, natural distribution, and phylogeny. The selected species have simple leaves, broad distribution in the Neotropics, and belong to five botanical families (Table 1).

Table 1
List of the seven studied species with their scientific names, family and age of seedlings grown in the Floresta Viva Institute nursery, Serra Grande, Bahia, Brazil.

Leaves were collected in the Floresta Viva Institute nursery, Serra Grande, Uruçuca, Bahia, Brazil, from March to October 2012. The leaves collected were all fully expanded, had their edges preserved and did not show herbivory or disease symptoms. For seedlings, 50 leaves were randomly collected from the same number of individuals (one leaf fully expanded and a mature leaf per plant). Seedlings were growing at the nursery of Floresta Viva Institute under 50% artificial shading. They were chosen for the study based on the presence of at least six fully expanded leaves and height between 0.15 and 0.25 m. For adults, 50 leaves, for species, were randomly collected from three to six trees of different sizes. The criterion used to select the adult individuals sampled was sexual maturity, ensuring that all plants were adults. Leaves were collected from different parts of the interior and top of the tree, in order to obtain greater morphological variability for sampling.

Harvested leaves were taken to the Plant Physiology Laboratory of the Universidade Estadual de Santa Cruz (UESC), Ilhéus, Bahia, Brazil, where they were washed and dried with paper towel, and had the image of the adaxial leaf surface digitalized. For leaves with length less than 35 cm, images of the adaxial leaf surface were scanned using a desktop scanner (HP Scanjet 2400 scanner), and digitalized in TIFF format with 100 dpi resolution. For leaves with length greater than 35 cm, adaxial leaf surface was photographed with a digital camera (Sony Alpha 200), and saved in JPEG format with 72 dpi resolution. A wooden white plate covered with an anti-reflective glass was used as the background of the photos. In the corner of the plate, a figure of a black cardboard, with an area of 36 cm2 (6 cm × 6 cm), was fixed to calibrate the image scale. Images were then analyzed using the free software ImageJ (Abramoff et al., 2004ABRAMOFF, M.D., MAGALHAES, P.J. and RAM, S.J., 2004. Image Processing with ImageJ. Biophotonics International, vol. 11, no. 7, pp. 36-42.). For each leaf, the leaf length (L), maximum leaf width (W) and LA were measured.

To determine the feasibility of using a single model for seedlings and adults of the same species, we performed an analysis of covariance (ANCOVA), comparing values in terms of a categorical variable with two levels. The use of a single model was chosen when there was no significant difference (p > 0.10) between the models for each of the ontogenetic stages (Logan, 2010Logan, M., 2010. Biostatistical design and analysis using R: a practical guide. Oxford: Wiley-Blackwell.).

The models for predicting LA were obtained by linear regression analysis, considering as dependent variable the observed LA and as independent variables the leaf linear dimensions (L2 and W2) and the product LW. Individual models were developed for each of the species and ontogenetic stages (seedlings and adults). The coefficient of determination (R2) was used as a measure of goodness-of-fit. Also, to assess the goodness-of-fit of the equations, the relative root mean square errors (RRMSE, %) were calculated dividing the root of mean square errors (RMSE) by the mean values of measured leaf areas, and multiplying the resulting value by 100 (Richardson et al., 2002Richardson, A.D., Duigan, S.P. and Berlyn, G.P., 2002. An evaluation of noninvasive methods to estimate foliar chlorophyll content. The New Phytologist, vol. 153, no. 1, pp. 185-194. http://dx.doi.org/10.1046/j.0028-646X.2001.00289.x.
http://dx.doi.org/10.1046/j.0028-646X.20...
). The predictive performance of the regression models were validated using the cross-validation method (Efron and Tibshirani, 1993Efron, B. and Tibshirani, R.J., 1993. An introduction to the bootstrap. London: Chapman & Hall.). For this, the data was randomly divided into k groups (in this work we used k = 10) and the model was adjusted k times, so that each time one of k groups is used as a test set for the model that was adjusted with corresponding training sets (the remaining n - k data division of each data set). It was expected that the estimate of the prediction error obtained by cross-validation based on model adjustment with only a portion of the data is greater than the estimated error in the model, when all observations are used in the fit. The relative difference between the error of the full model and the estimated error of the cross-validation (prediction error PE) represents a measure of the predictive accuracy of the fitted model (Maindonald and Braun, 2003Maindonald, J. and Braun, J., 2003. Data analysis and graphics using R: an example-based approach. Cambridge: Cambridge University Press.).

For the species in which it was possible to obtain single models for the two ontogenetic stages, the models were validated using a second sample of leaves collected in the same place, but from different plants and seedling lots, in August 2014. The validation of the models was done plotting the values of predicted (PLA) and observed (OLA) values of LA. In addition, the relative root mean square calibration errors (RRMSEc, %) were calculated for seedlings and adults. The values of RRMSEc were calculated dividing the root of mean square calibration errors (RMSEc) of the OLA and PLA by the mean values of OLA, and multiplying the resulting value by 100 (Richardson et al., 2002Richardson, A.D., Duigan, S.P. and Berlyn, G.P., 2002. An evaluation of noninvasive methods to estimate foliar chlorophyll content. The New Phytologist, vol. 153, no. 1, pp. 185-194. http://dx.doi.org/10.1046/j.0028-646X.2001.00289.x.
http://dx.doi.org/10.1046/j.0028-646X.20...
).

3 Results

For all species, the average values of L, W and LA were significantly higher in adults than in seedlings (p < 0.01, t-test, Table 2). Leaves of T. stenocalyx and P. torta showed the highest values of L, W and LA for adults. For seedlings, the highest values of L, W and LA were exhibited by P. torta. The smallest values of L, W and LA were found in the leaves of B. rubescens, P. caimito and S. globulifera, respectively. The highest and smallest differences of the measured values of L, W and LA between adults and seedlings were observed for T. stenocalyx and S. globulifera, respectively. For T. stenocalyx, the mean values of L, W and LA were respectively 3.8, 4.9 and 20.6-fold higher in adults than seedlings. In contrast, for S. globulifera the mean values of L, W and LA were respectively, 1.1-, 1.3- and 1.5-fold higher in adults than seedlings. Leaf shapes (L:W) varied from 1.95 (P. cattleyanum, adults) to 3.73 (P. caimito, seedlings). The L:W ratio was higher in seedlings than in adults for most species, with exception of M. maxima and P. torta. The highest difference between adults and seedlings for L:W was obtained for P. torta (adult/seedling = 1.57) and the smallest difference for M. maxima and P. cattleyanum (adult/seedling = 1.12).

Table 2
Mean values (minimum - maximum) of maximum leaf length (L), width (W), leaf area (LA) and leaf shape (L:W) for seedlings and adults of seven rainforest tree species. Ontogenetic stage (OS): A – adult, S – seedling (n = 50).

The use of a single model for both ontogenetic stages (seedling and adult) was possible only for two species (P. cattleyanum and T. stenocalyx) (Table 3). The single model based on the product LW showed the high values of R2 (> 0.99), followed by the RRMSE below 6.0% and PE below 7.0%, are indicatives of the high accuracy of this models (Table 3). For T. stenocalyx, a single model also proved feasible based on L2, but the RRMSE and PE values are higher than 35% and 16%, respectively (Table 4). For the two species, only the models based on LW were tested against a new set of leaves (Table 5). When the predicted (PLA) and observed (OLA) values of LA from a new set of leaves were plotted, high values of R2 and low RRMSEc were observed for both species (Figure 1, Table 5).

Table 3
Species, values of p-interaction and p-level (ANCOVA) for leaves of seedling and adult individuals (n = 100).
Table 4
Integrative models for leaves of seedling and adult individuals of P. cattleyanum and T. stenocalyx, with their R2, p-values, RRMSE and PE (n = 100).
Table 5
Relationships between predicted (PLA) and observed (OLA) leaf areas for seedlings and adults trees of P. cattleyanum and T. stenocalyx. The models used to calculate PLA are in Table 4.
Figure 1
Relationships between predicted (PLA) and observed (OLA) leaf areas for seedlings and adults trees of P. cattleyanum (a) and T. stenocalyx (b). Solid line represent the 1:1 relationship between PLA and OLA. The models used to calculate PLA are in Table 4.

4 Discussion

Seedlings of the same tree species have smaller LA than adults because they have smaller size. As expected, the values of leaf dimensions (L, W and LA) were much higher in adults than in seedlings for all the seven species analyzed in our study. Moreover, the differences between average values of LA for adults and seedlings of T. stenocalyx and M. maxima indicate a large variation in the leaf dimensions among tree species according to the ontogenic stage. Ishida et al. (2005)Ishida, A., Yazaki, K. and Hoe, A.L., 2005. Ontogenetic transition of leaf physiology and anatomy from seedlings to mature trees of a rain forest pioneer tree, Macaranga gigantea. Tree Physiology, vol. 25, no. 5, pp. 513-522. http://dx.doi.org/10.1093/treephys/25.5.513. PMid:15741149.
http://dx.doi.org/10.1093/treephys/25.5....
, for instance, found that LA of adult individuals of Macaranga gigantea was around 10 times higher than seedlings; an intermediate value between the extremes found in our study.

Compared with L, W and LA, very small differences were observed between seedlings and adults for the average values of L:W. We are unaware of published studies reporting differences of L:W between seedlings and adults for tropical trees. Comparing new and old leaves of eight evergreen subtropical trees in China, Zhang and Pan (2011)Zhang, L. and Pan, L., 2011. Allometric models for leaf area estimation across different leaf-age groups of evergreen broadleaved trees in a subtropical forest. Photosynthetica, vol. 49, no. 2, pp. 219-226. http://dx.doi.org/10.1007/s11099-011-0027-x.
http://dx.doi.org/10.1007/s11099-011-002...
reported values of L:W very close of those we obtained in our study. It is noteworthy that in our study, only two of the seven species (M. maxima and P. torta) had values of L:W greater in adults than in seedlings. In addition, between the two species of the genus Pouteria, the values of L:W ratio are opposed. The average value of L:W was higher in seedlings than in adults for P. torta, and higher in seedlings than in adults for P. caimito. As an indicative of the leaf shapes, the average values of L:W in our study indicate that all leaves tends to have elliptical shapes, as the average values for adults and seedlings were, respectively, 2.5 and 2.8. These values are very close of those reported by Zhang and Pan (2011)Zhang, L. and Pan, L., 2011. Allometric models for leaf area estimation across different leaf-age groups of evergreen broadleaved trees in a subtropical forest. Photosynthetica, vol. 49, no. 2, pp. 219-226. http://dx.doi.org/10.1007/s11099-011-0027-x.
http://dx.doi.org/10.1007/s11099-011-002...
and higher than the values obtained by Rouphael et al. (2010)Rouphael, Y., Mouneimne, A.H., Ismail, A., Mendoza-Degyves, E., Rivera, C.M. and Colla, G., 2010. Modeling individual leaf area of rose (Rosa hybrida L.) based on leaf length and width measurement. Photosynthetica, vol. 48, no. 1, pp. 9-15. http://dx.doi.org/10.1007/s11099-010-0003-x.
http://dx.doi.org/10.1007/s11099-010-000...
for Rosa hybrida. Moreover, the ratios between adults and seedlings for L:W in our study were low and much smaller than the differences obtained for LA. In this study, the differences between adults and seedlings for L:W ranged from 1.1, for M. maxima and P. cattleyanum, to 1.6, for P. torta.

There are several published studies of nondestructive estimations of LA involving measurements of leaf dimensions as independent variables, but using different regression models, including power (Kandiannan et al., 2002Kandiannan, K., Kailasam, C., Chandaragiri, K.K. and Sankaran, N., 2002. Allometric model for leaf area estimation in black pepper (Piper nigrum L.). Journal Agronomy & Crop Science, vol. 188, no. 2, pp. 138-140. http://dx.doi.org/10.1046/j.1439-037X.2002.00540.x.
http://dx.doi.org/10.1046/j.1439-037X.20...
; Antunes et al., 2008Antunes, W.C., Pompelli, M.F., Carretero, D.M. and Damatta, F.M., 2008. Allometric models for non-destructive leaf area estimation in coffe (Coffea arabica and Coffea canephora). Annals of Applied Biology, vol. 153, no. 1, pp. 33-40. http://dx.doi.org/10.1111/j.1744-7348.2008.00235.x.
http://dx.doi.org/10.1111/j.1744-7348.20...
; Pompelli et al., 2012Pompelli, M.F., Antunes, W.C., Ferreira, D.T.R.G., Cavalcante, P.G.S., Wanderley-Filho, H.C.L. and Endres, L., 2012. Allometric models for non-destructive leaf area estimation of Jatropha curcas. Biomass and Bioenergy, vol. 36, pp. 77-85. http://dx.doi.org/10.1016/j.biombioe.2011.10.010.
http://dx.doi.org/10.1016/j.biombioe.201...
; Wang and Zhang, 2012Wang, Z. and Zhang, L., 2012. Leaf shape alters the coefficients of leaf area estimation models for Saussurea stoliczkai in central Tibet. Photosynthetica, vol. 50, no. 3, pp. 337-342. http://dx.doi.org/10.1007/s11099-012-0039-1.
http://dx.doi.org/10.1007/s11099-012-003...
; Misle et al., 2013Misle, E., Kahlaoui, B., Hachicha, M. and Alvarado, P., 2013. Leaf area estimation in muskmelon by allometry. Photosynthetica, vol. 51, no. 4, pp. 613-620. http://dx.doi.org/10.1007/s11099-013-0062-x.
http://dx.doi.org/10.1007/s11099-013-006...
) and second-order polynomial models (Mazzini et al., 2010Mazzini, R.B., Ribeiro, R.V. and Pio, R.M., 2010. A simple and non-destructive model for individual leaf area estimation in citrus. Fruits, vol. 65, no. 5, pp. 269-275. http://dx.doi.org/10.1051/fruits/2010022.
http://dx.doi.org/10.1051/fruits/2010022...
). In our study, we adopted the use of linear regression models because they are simpler and provide good estimates for a great number of plant species with different leaf areas and leaf shapes when compared with non-linear regression models (Cristofori et al., 2007Cristofori, V., Rouphael, Y., Mendoza-gyves, E. and Bignami, C., 2007. A simple model for estimating leaf area of hazelnut from linear measurements. Scientia Horticulturae, vol. 113, no. 2, pp. 221-225. http://dx.doi.org/10.1016/j.scienta.2007.02.006.
http://dx.doi.org/10.1016/j.scienta.2007...
; Mazzini et al., 2010Mazzini, R.B., Ribeiro, R.V. and Pio, R.M., 2010. A simple and non-destructive model for individual leaf area estimation in citrus. Fruits, vol. 65, no. 5, pp. 269-275. http://dx.doi.org/10.1051/fruits/2010022.
http://dx.doi.org/10.1051/fruits/2010022...
; Rouphael et al., 2010Rouphael, Y., Mouneimne, A.H., Ismail, A., Mendoza-Degyves, E., Rivera, C.M. and Colla, G., 2010. Modeling individual leaf area of rose (Rosa hybrida L.) based on leaf length and width measurement. Photosynthetica, vol. 48, no. 1, pp. 9-15. http://dx.doi.org/10.1007/s11099-010-0003-x.
http://dx.doi.org/10.1007/s11099-010-000...
; Wang and Zhang, 2012Wang, Z. and Zhang, L., 2012. Leaf shape alters the coefficients of leaf area estimation models for Saussurea stoliczkai in central Tibet. Photosynthetica, vol. 50, no. 3, pp. 337-342. http://dx.doi.org/10.1007/s11099-012-0039-1.
http://dx.doi.org/10.1007/s11099-012-003...
; Souza and Habermann, 2014Souza, M.C. and Habermann, G., 2014. Non-destructive equations to estimate the leaf area of Styrax pohliiand and Styrax ferrugineus. Brazilian Journal of Biology = Revista Brasileira de Biologia, vol. 74, no. 1, pp. 222-225. http://dx.doi.org/10.1590/1519-6984.17012. PMid:25055106.
http://dx.doi.org/10.1590/1519-6984.1701...
). In addition, we tested two geometric formulas for estimating LA from simple linear regression models. Initially, we attempted to use simple linear measurements (L and W), based on the square area formula (L2 and W2). This approach was tested with the purpose of reducing the effort for field measurements. Using only a single linear measurement (L or W) the time needed for field measurements is half of what is needed considering two linear measurements (L and W). The second approach was the use of the rectangle area (L × W). Although more time and effort is needed for collecting data, the use of the rectangle area (LW) showed better adjustments for predicting LA. These results were closely related with the leaf shapes, as the values of the inclinations were about 0.7 for L × W and 0.3 for L2 and W2.

The use of a single model, regardless of the plant ontogenetic stage, may represent less work and high reliability. For instance, leaves of T. stenocalyx showed the highest differences in L, W and LA between adults and seedlings, compared to the other species. The average LA of adults was almost 21-fold higher than seedlings. The highest and smallest measured values of LA in seedlings and adults of this species were 35.8 cm2 and 88.7 cm2, respectively. Moreover, the difference between adults and seedlings for L:W was exactly the average found for the set of all seven species (1.3). Therefore, the development of simple predictive models for LA for both ontogenetic stages seems to be more dependent on the species than the differences in leaf size between seedlings and adults.

The development of regression models encompassing different ontogenetic stages for tropical tree species may have numerous practical advantages in field studies of ecology and ecophysiology of tropical trees. However, in tropical rainforests, seedlings, saplings and adults can vary widely in leaf shape and leaf areas (Thomas and Ickes, 1995Thomas, S.C. and Ickes, K., 1995. Ontogenetic changes in leaf size in Malaysian Rain Forest Trees. Biotropica, vol. 27, no. 4, pp. 427-434. http://dx.doi.org/10.2307/2388954.
http://dx.doi.org/10.2307/2388954...
; Ishida et al., 2005Ishida, A., Yazaki, K. and Hoe, A.L., 2005. Ontogenetic transition of leaf physiology and anatomy from seedlings to mature trees of a rain forest pioneer tree, Macaranga gigantea. Tree Physiology, vol. 25, no. 5, pp. 513-522. http://dx.doi.org/10.1093/treephys/25.5.513. PMid:15741149.
http://dx.doi.org/10.1093/treephys/25.5....
; Houter and Pons, 2012Houter, N. and Pons, T.L., 2012. Ontogenetic changes in leaf traits of tropical rainforest trees differing in juvenile light requirement. Oecologia, vol. 169, no. 1, pp. 33-45. http://dx.doi.org/10.1007/s00442-011-2175-x. PMid:22038060.
http://dx.doi.org/10.1007/s00442-011-217...
). Our results demonstrate that it is not always possible to use models developed for adults in studies with seedlings, and vice-versa. This fact is very critical, since most of the models available in the literature have been developed based upon leaves of adult plants (Cristofori et al., 2007Cristofori, V., Rouphael, Y., Mendoza-gyves, E. and Bignami, C., 2007. A simple model for estimating leaf area of hazelnut from linear measurements. Scientia Horticulturae, vol. 113, no. 2, pp. 221-225. http://dx.doi.org/10.1016/j.scienta.2007.02.006.
http://dx.doi.org/10.1016/j.scienta.2007...
; Mazzini et al., 2010Mazzini, R.B., Ribeiro, R.V. and Pio, R.M., 2010. A simple and non-destructive model for individual leaf area estimation in citrus. Fruits, vol. 65, no. 5, pp. 269-275. http://dx.doi.org/10.1051/fruits/2010022.
http://dx.doi.org/10.1051/fruits/2010022...
; Rouphael et al., 2010Rouphael, Y., Mouneimne, A.H., Ismail, A., Mendoza-Degyves, E., Rivera, C.M. and Colla, G., 2010. Modeling individual leaf area of rose (Rosa hybrida L.) based on leaf length and width measurement. Photosynthetica, vol. 48, no. 1, pp. 9-15. http://dx.doi.org/10.1007/s11099-010-0003-x.
http://dx.doi.org/10.1007/s11099-010-000...
; Wang and Zhang, 2012Wang, Z. and Zhang, L., 2012. Leaf shape alters the coefficients of leaf area estimation models for Saussurea stoliczkai in central Tibet. Photosynthetica, vol. 50, no. 3, pp. 337-342. http://dx.doi.org/10.1007/s11099-012-0039-1.
http://dx.doi.org/10.1007/s11099-012-003...
; Souza and Habermann, 2014Souza, M.C. and Habermann, G., 2014. Non-destructive equations to estimate the leaf area of Styrax pohliiand and Styrax ferrugineus. Brazilian Journal of Biology = Revista Brasileira de Biologia, vol. 74, no. 1, pp. 222-225. http://dx.doi.org/10.1590/1519-6984.17012. PMid:25055106.
http://dx.doi.org/10.1590/1519-6984.1701...
). The use of models developed with leaves of adults in studies with seedlings can cause serious errors in the data analysis and interpretation. Furthermore, in this study, we only test square-based (L2, W2) and rectangular (L × W) formulas and adopted adjustments based on simple linear regression models. Considering the changes in the leaf shapes due to ontogenetic changes, we also encourage the development of further studies, involving a larger number of species with different geometric equations and regression formulas.

In summary, our results indicate that from the seven species analyzed, only two (P. cattleyanum and T. stenocalyx) had very similar relationships between LW and LA in both ontogenetic stages. Our study shows that it is possible to develop predictive models for estimating leaf area encompassing different ontogenetic stages for tropical trees, but the reliability of these models depends more on the species itself than the differences in leaf size between seedlings and adults.

Acknowledgements

We thank Gerson J. Sales Neto, Nilson A. dos Santos and Rones F. Souza of Floresta Viva Institute for assistance with data collection. We also thank Dr. Maria C. Sanches of Universidade Federal de Uberlândia for her helpful suggestions on the manuscript. Funding for Ediófila B. Rocha during this study was provided by a scholarship from Capes (Brazilian Higher Education Council). The corresponding author gratefully acknowledges CNPq (Brazilian National Council for Scientific and Technological Development), Brazil, for the award of a fellowship of scientific productivity. This study was supported by CNPq (Proc. 561933/2010-3).

  • (With 1 figure)

References

  • ABRAMOFF, M.D., MAGALHAES, P.J. and RAM, S.J., 2004. Image Processing with ImageJ. Biophotonics International, vol. 11, no. 7, pp. 36-42.
  • Antunes, W.C., Pompelli, M.F., Carretero, D.M. and Damatta, F.M., 2008. Allometric models for non-destructive leaf area estimation in coffe (Coffea arabica and Coffea canephora). Annals of Applied Biology, vol. 153, no. 1, pp. 33-40. http://dx.doi.org/10.1111/j.1744-7348.2008.00235.x
    » http://dx.doi.org/10.1111/j.1744-7348.2008.00235.x
  • Boege, K. and Marquis, R.J., 2005. Facing herbivory as you grow up: the ontogeny of resistance in plants. Trends in Ecology & Evolution, vol. 20, no. 8, pp. 441-448. http://dx.doi.org/10.1016/j.tree.2005.05.001 PMid:16701415.
    » http://dx.doi.org/10.1016/j.tree.2005.05.001
  • Coste, S., Roggy, J.C., Schimann, H., Epron, D. and Dreyer, E., 2011. A cost–benefit analysis of acclimation to low irradiance in tropical rainforest tree seedlings: leaf life span and payback time for leaf deployment. Journal of Experimental Botany, vol. 62, no. 11, pp. 3941-3955. http://dx.doi.org/10.1093/jxb/err092 PMid:21511904.
    » http://dx.doi.org/10.1093/jxb/err092
  • Cristofori, V., Rouphael, Y., Mendoza-gyves, E. and Bignami, C., 2007. A simple model for estimating leaf area of hazelnut from linear measurements. Scientia Horticulturae, vol. 113, no. 2, pp. 221-225. http://dx.doi.org/10.1016/j.scienta.2007.02.006
    » http://dx.doi.org/10.1016/j.scienta.2007.02.006
  • Efron, B. and Tibshirani, R.J., 1993. An introduction to the bootstrap. London: Chapman & Hall.
  • Godoy, L.J.G., Yanagiwara, R.S., Bôas, R.L.V., Backes, C. and Lima, C.P., 2007. Análise da imagem digital para estimativa da área foliar em plantas de laranja “Pera”. Revista Brasileira de Fruticultura, vol. 29, no. 3, pp. 420-424. http://dx.doi.org/10.1590/S0100-29452007000300004
    » http://dx.doi.org/10.1590/S0100-29452007000300004
  • Houter, N. and Pons, T.L., 2012. Ontogenetic changes in leaf traits of tropical rainforest trees differing in juvenile light requirement. Oecologia, vol. 169, no. 1, pp. 33-45. http://dx.doi.org/10.1007/s00442-011-2175-x PMid:22038060.
    » http://dx.doi.org/10.1007/s00442-011-2175-x
  • Ishida, A., Yazaki, K. and Hoe, A.L., 2005. Ontogenetic transition of leaf physiology and anatomy from seedlings to mature trees of a rain forest pioneer tree, Macaranga gigantea. Tree Physiology, vol. 25, no. 5, pp. 513-522. http://dx.doi.org/10.1093/treephys/25.5.513 PMid:15741149.
    » http://dx.doi.org/10.1093/treephys/25.5.513
  • Kandiannan, K., Kailasam, C., Chandaragiri, K.K. and Sankaran, N., 2002. Allometric model for leaf area estimation in black pepper (Piper nigrum L.). Journal Agronomy & Crop Science, vol. 188, no. 2, pp. 138-140. http://dx.doi.org/10.1046/j.1439-037X.2002.00540.x
    » http://dx.doi.org/10.1046/j.1439-037X.2002.00540.x
  • Laurans, M., Martin, O., Nicolini, E. and Vincent, G., 2012. Functional traits and their plasticity predict tropical trees regeneration niche even among species with intermediate light requirements. Journal of Ecology, vol. 100, no. 6, pp. 1440-1452. http://dx.doi.org/10.1111/j.1365-2745.2012.02007.x
    » http://dx.doi.org/10.1111/j.1365-2745.2012.02007.x
  • Logan, M., 2010. Biostatistical design and analysis using R: a practical guide. Oxford: Wiley-Blackwell.
  • Maindonald, J. and Braun, J., 2003. Data analysis and graphics using R: an example-based approach. Cambridge: Cambridge University Press.
  • Mazzini, R.B., Ribeiro, R.V. and Pio, R.M., 2010. A simple and non-destructive model for individual leaf area estimation in citrus. Fruits, vol. 65, no. 5, pp. 269-275. http://dx.doi.org/10.1051/fruits/2010022
    » http://dx.doi.org/10.1051/fruits/2010022
  • Misle, E., Kahlaoui, B., Hachicha, M. and Alvarado, P., 2013. Leaf area estimation in muskmelon by allometry. Photosynthetica, vol. 51, no. 4, pp. 613-620. http://dx.doi.org/10.1007/s11099-013-0062-x
    » http://dx.doi.org/10.1007/s11099-013-0062-x
  • Niinemets, Ü., 2010. A review of light interception in plant stands from leaf to canopy in different plant functional types and in species with varying shade tolerance. Ecological Research, vol. 25, no. 4, pp. 693-714. http://dx.doi.org/10.1007/s11284-010-0712-4
    » http://dx.doi.org/10.1007/s11284-010-0712-4
  • O’Neal, M.E., Landis, D.A. and Isaacs, R., 2002. An inexpensive, accurate method for measuring leaf area and defoliation through digital image analysis. Journal of Economic Entomology, vol. 95, no. 6, pp. 1190-1194. http://dx.doi.org/10.1603/0022-0493-95.6.1190 PMid:12539831.
    » http://dx.doi.org/10.1603/0022-0493-95.6.1190
  • Pompelli, M.F., Antunes, W.C., Ferreira, D.T.R.G., Cavalcante, P.G.S., Wanderley-Filho, H.C.L. and Endres, L., 2012. Allometric models for non-destructive leaf area estimation of Jatropha curcas. Biomass and Bioenergy, vol. 36, pp. 77-85. http://dx.doi.org/10.1016/j.biombioe.2011.10.010
    » http://dx.doi.org/10.1016/j.biombioe.2011.10.010
  • Richardson, A.D., Duigan, S.P. and Berlyn, G.P., 2002. An evaluation of noninvasive methods to estimate foliar chlorophyll content. The New Phytologist, vol. 153, no. 1, pp. 185-194. http://dx.doi.org/10.1046/j.0028-646X.2001.00289.x
    » http://dx.doi.org/10.1046/j.0028-646X.2001.00289.x
  • Rouphael, Y., Mouneimne, A.H., Ismail, A., Mendoza-Degyves, E., Rivera, C.M. and Colla, G., 2010. Modeling individual leaf area of rose (Rosa hybrida L.) based on leaf length and width measurement. Photosynthetica, vol. 48, no. 1, pp. 9-15. http://dx.doi.org/10.1007/s11099-010-0003-x
    » http://dx.doi.org/10.1007/s11099-010-0003-x
  • Souza, M.C. and Habermann, G., 2014. Non-destructive equations to estimate the leaf area of Styrax pohliiand and Styrax ferrugineus. Brazilian Journal of Biology = Revista Brasileira de Biologia, vol. 74, no. 1, pp. 222-225. http://dx.doi.org/10.1590/1519-6984.17012 PMid:25055106.
    » http://dx.doi.org/10.1590/1519-6984.17012
  • Thomas, S.C. and Ickes, K., 1995. Ontogenetic changes in leaf size in Malaysian Rain Forest Trees. Biotropica, vol. 27, no. 4, pp. 427-434. http://dx.doi.org/10.2307/2388954
    » http://dx.doi.org/10.2307/2388954
  • Thomas, W.W., 1999. Conservation and monographic research on the flora of tropical America. Biodiversity and Conservation, vol. 8, no. 8, pp. 1007-1015. http://dx.doi.org/10.1023/A:1008857429787
    » http://dx.doi.org/10.1023/A:1008857429787
  • Toebe, M., Brum, B., Lopes, S.J., Cargnelutti-Filho, A. and Silveira, T.R., 2010. Estimativa da área foliar de Crambe abyssinica por discos foliares e por fotos digitais. Ciência Rural, vol. 40, no. 2, pp. 475-478. http://dx.doi.org/10.1590/S0103-84782010000200036
    » http://dx.doi.org/10.1590/S0103-84782010000200036
  • Tsuda, M., 1999. Errors in leaf area measurement with an automatic area meter due to leaf chlorophyll in crop plants. Annals of Botany, vol. 84, no. 6, pp. 799-801. http://dx.doi.org/10.1006/anbo.1999.0984
    » http://dx.doi.org/10.1006/anbo.1999.0984
  • Varone, L. and Gratani, L., 2009. Leaf expansion in Rhamnus alaternus L. by leaf morphological, anatomical and physiological analysis. Trees (Berlin), vol. 23, no. 6, pp. 1255-1262. http://dx.doi.org/10.1007/s00468-009-0365-5
    » http://dx.doi.org/10.1007/s00468-009-0365-5
  • Wang, Z. and Zhang, L., 2012. Leaf shape alters the coefficients of leaf area estimation models for Saussurea stoliczkai in central Tibet. Photosynthetica, vol. 50, no. 3, pp. 337-342. http://dx.doi.org/10.1007/s11099-012-0039-1
    » http://dx.doi.org/10.1007/s11099-012-0039-1
  • Zhang, L. and Pan, L., 2011. Allometric models for leaf area estimation across different leaf-age groups of evergreen broadleaved trees in a subtropical forest. Photosynthetica, vol. 49, no. 2, pp. 219-226. http://dx.doi.org/10.1007/s11099-011-0027-x
    » http://dx.doi.org/10.1007/s11099-011-0027-x

Publication Dates

  • Publication in this collection
    17 May 2016
  • Date of issue
    Oct-Dec 2016

History

  • Received
    28 Mar 2015
  • Accepted
    15 Sept 2015
Instituto Internacional de Ecologia R. Bento Carlos, 750, 13560-660 São Carlos SP - Brasil, Tel. e Fax: (55 16) 3362-5400 - São Carlos - SP - Brazil
E-mail: bjb@bjb.com.br