Adjustment of Mathematical Equations to Determine the Area Leaf of Cordia myxa L.

Michelle Teixeira Araujo¹, Ellen Maurício Faria¹, Gabriel Linhares Cunha¹, Paloma de Moraes Nunes¹, Vinicius de Souza Oliveira², Karina Tiemi Hassuda dos Santos¹, Omar Schmildt¹, Antelmo Ralph Falqueto¹, Mônica Maria Pereira Tognella¹, Edilson Romais Schmildt¹

¹Center North of Espírito Santo, Federal University of Espírito Santo, São Mateus, Espírito Santo, Brazil.
²Center for Agricultural Sciences and Engineering, Federal University of Espírito Santo, Alegre, Espírito Santo, Brazil.

DOI: 10.4236/as.2020.117038                

Abstract

The objective of this study was to obtain mathematical equations through the leaf dimensions to estimate the area leaf of Cordia myxa L. 392 leaves of plants were collected, located in the municipality of São Mateus, North of the State of Espírito Santo, Brazil, 300 of which were used to adjust the equations and 92 for validation. Of all the leaves, the largest length (L) and maximum width (W) were measured, product of multiplying the length with the width (LW) and the real leaf area (RLA). Models equation first degree linear, quadratic and power were fitted, where RLA was the dependent variable as a function of L, W and LW as an independent variable. All equations were validated according to appropriate criterion. Thus, the power model ELA = 0.84(LW)^0.9749, based on the product of multiplication of length with width (LW) is the most suitable for estimating the leaf area of Cordia myxa L.

Keywords

Non-Destructive Method, Validated, Modeling

1. Introduction

Cordia myxa L. is a species native to Asia, popularly known as Assyrian plum and belonging to the Boraginaceae family, widespread especially in the tropics and temperate regions of the world, an ecological characteristic present in Boraginaceae is the ability of many of its species to establish themselves in areas very disturbed, therefore being useful species for the reforestation of areas where the anthropic action is evident [1].

The measurement of the area leaf can be an important parameter in studies related to plant morphology, anatomy and ecophysiology, as it allows to obtain a fundamental indicator for understanding the plant’s responses to environmental factors, since the leaves are the main organs responsible for photosynthetic process, its knowledge being able to estimate the loss of water through transpiration, in addition to the absorption and conversion of light energy into chemical energy [2] [3].

The determination of the area leaf can be carried out by direct or indirect, destructive or non-destructive methods. The direct methods are those that use measurements made directly on the leaves, where they often require the removal of the leaf or other structures, which may not be possible, due to the limited number of plants, in the experimental plot, in its turn the indirect methods are based on the known correlation between a measured variable and the leaf area, obtaining precise and easy mathematical equations, this being a non-destructive method and the measures can be taken on the plant, without the need to remove structures [4].

The use of mathematical equations to estimate leaf area has been reported for several plant species [5] - [11]. However, no literature is lack of studies that seek to model the area leaf for Cordia myxa L.

Thus, the objective of this study was to obtain mathematical equations and indicate the most appropriate for estimating the area leaf of Cordia myxa L. in a non-destructive way from the linear dimensions of the leaves.

2. Material and Methods

The study was carried out at the Federal University of Espírito Santo, Campus São Mateus, located in the North of the State of Espírito Santo, Brazil, with the following geographical coordinates: 18˚40'36'' south latitude and 39˚51'35'' longitude East. The climate of the region is characterized according to Köppen as tropical AW, with dry winter and rain during the summer [12].

The experiment was based on the collection of 392 leaves of Cordia myxa L. healthy and at different stages of development, in the four cardinal points, which did not present damage or attack from diseases or pests, as suggested by Oliveira et al. [6], soon the leaves were packed in plastic bags and quickly transported to the laboratory. In the laboratory, of each leaf, the largest length (L) and the largest width (W) of the leaf blade was determined using a ruler graduated in cm (Figure 1). The length was measured along the main midrib, which is the distance between the base of the leaf at the point of insertion of the petiole to its apex, and the width, considered in the middle part of the leaf. Subsequently, the real leaf area (RLA) of each leaf was determined in cm², with the aid of the Licor LI-3100 table meter. The product of multiplying the length with the width (LW in cm²) was also calculated. Descriptive statistics of the data were determined by obtaining the values of maximum, minimum, mean, amplitude, standard deviation (SD) and coefficient of variation (CV).

Figure 1. Representation of the length (L) along the midrib and the maximum width (W) of leaves of Cordia myxa L.

For the modeling, 300 leaves were used in which RLA was defined with the dependent variable as a function of L, W or LW as independent variables, being tested the models first degree linear represented by $ELA={\hat{\beta }}_0+{\hat{\beta }}_{1}x$, quadratic represented by $ELA={\hat{\beta }}_0+{\hat{\beta }}_{1}x+{\hat{\beta }}_2{x}^2$ and power represented by $ELA={\hat{\beta }}_0{x}^{{\hat{\beta }}_1}$, totaling nine equations and their respective coefficient of determination (R²).

For validation, the L, W and LW values of 92 leaves that were sampled for this purpose were replaced in the equations proposed in the modeling, thus obtaining the estimated leaf area (AFE), in cm². A simple linear equation ( $ELA={\hat{\beta }}_0+{\hat{\beta }}_{1}x$ ) was thus adjusted, where ELA was used as a dependent variable as a function of RLA. The hypotheses $H_0:{\beta }_0=0$ versus $H_a:{\beta }_0\ne 0$ and $H_0:{\beta }_1=1$ versus $H_a:{\beta }_1\ne 1$ were tested, using the Student t test at 5% probability. The mean absolute error (MAE), the root of the root mean square error (RMSE) and the Willmott d index [13] were also determined for all equations estimated using expressions 1, 2 and 3.

$MAE=\frac{{{\displaystyle \sum }}_{i=1}^n\left|ELA-RLA\right|}{n}$ (1)

$RMSE=\sqrt{\frac{{{\displaystyle \sum }}_{i=1}^n{\left(ELA-RLA\right)}^2}{n}}$ (2)

$d=1-\left[\frac{{{\displaystyle \sum }}_{i=1}^n{\left(ELA-RLA\right)}^2}{{{\displaystyle \sum }}_{i=1}^n{\left(\left|ELA-\overline{RLA}\right|+\left|RLA-\overline{RLA}\right|\right)}^2}\right]$ (3)

where: ELA, are the estimated values of leaf area; RLA, are the real values of leaf area; $\overline{RLA}$ is the average of the real leaf area values; n, is the number of leaves sampled for validation, 92 in the present study.

The criterion for selecting the best models was based on linear coefficient ( ${\beta }_0$ ) is not different from zero, slope ( ${\beta }_1$ ) is not different from one, values of the mean absolute error (MAE) and root mean square error (RMSE) closest to zero and Willmott’s d index [13] closest to one.

Statistical analyzes were performed with the aid of the R software [14], with scripts developed for the data package ExpDes.pt version 1.2 [15].

3. Results and Discussion

High variability of the sampled data was observed, verified by the high amplitude, with values above the average for all characteristics (Table 1), in addition, the values of the coefficient of variation (CV) were higher than 25.30% for both the sample used for modeling as for the sample used for modeling, these values being considered high or very high [16]. This high variation was also found by several authors when studying mathematical models to estimate the leaf area of other species [7] [10] [11]. These data point out that in the sample collected there are leaves of different stages of development and consequently of different sizes, which is a fundamental characteristic in the study of modeling, since the adjusted equations can later be used throughout the crop cycle [5] [8].

The nine equations adjusted to estimate the leaf area of Cordia myxa L. through their linear dimensions can be seen in Table 2. Note that there was a better relationship between the real leaf area (AFR) and the product of multiplying the length with the width (LW) in the three models studied with coefficient of determination (R²) greater than 0.98, showing that more than 98% of the area leaf of Cordia myxa L. leaves can be explained by LW. This better relationship can be explained since the length and width are different measures and when combined they become more appropriate than when used individually [17]. However, we must not pay attention only to the high values of R² when choosing the best equation, because its use, without any other validation criteria, may imply an erroneous estimate of the leaf area [18].

When analyzing the coefficients of the nine adjusted equations (Table 3), five equations showed results not significantly different from zero for the linear coefficient ( ${\beta }_0$ ) and not significantly different from one for the slope ( ${\beta }_1$ ). This finding suggests that for these equations, as 1 cm² in the real leaf area increases, the equation will add 1 cm² in the estimated leaf area and when the real leaf area is 0 cm², the model will estimate a leaf area of 0 cm² [9]. However, despite the

Table 1. Descriptive statistics with value minimum, maximum, mean, amplitude, standard deviation (SD) and coefficient of variation (CV) of the variables: length (L); width (W); product of the length and width (LW) and real leaf area (RLA) of leaves of Cordia myxa L.

Table 2. Equation with linear adjustment of first degree, quadratic and power and its respective coefficient of determination (R²) using the real leaf area (RLA) as dependent variable, in function of length (L), width (W), product of length with width (LW) of leaves of Cordia myxa L.

Table 3. Linear coefficient (β₀), slope (β₁) and determination coefficient (R²), obtained from equations with first line linear adjustment between real leaf area (RLA) and estimated leaf area (ELA) of linear equations of first degree, quadratic and power for the independent variables length (L), width (W) and product of length and width (LW), besides, mean absolute error (MAE), root mean square error (RMSE) and Willmott d index of leaves of Cordia myxa L. used for validation.

(1) ^nsValue of β₀ does not differ from zero by Student’s t-test, at a level of 5%; (1) *Value of β₀ differs from zero by Student’s t-test, at a level of 5%; (2) ^nsValue of β₁ does not differ from one, by Student’s t-test, at a level of 5%; (2) *Value of β₁ differs from one, by Student’s t-test, at a level of 5%.

fact that five equations have no significant effect for ${\beta }_0$ and ${\beta }_1$, the power model equations adjusted from the product of the length with the width proved to be more effective in determining the leaf area of Cordia myxa L. because it indicates less error in the estimation verified by the values of the mean absolute error (MAE) and root mean square error (RMSE) closer to zero, besides, of index values d closer to one.

Figure 2. Adjusted equation of power model and its coefficient of determination (R²), using the real leaf area (RLA) as dependent variable, as a function of the product of length and width (LW) of leaves of Cordia myxa L.

The equations adjusted with only one measurement of the leaf surface (L or W) are easier to use in practice [8]. However, this type of equation is less accurate and may require different adjustments throughout the vegetative cycle of the plants, thus requiring different equations for the different stages of development [19]. Thus, the results presented here show us that the equation $ELA=0.84{\left(LW\right)}^{0.9749}$ whose behavior can be seen in Figure 2 is the most accurate in estimating the leaf area of Cordia myxa L. This equation was chosen because it presents the best results when submitted to appropriate validation methods [20].

4. Conclusion

The power model equation represented by $ELA=0.84{\left(LW\right)}^{0.9749}$ adjusted from the product of the length with the width is the most accurate in estimating the leaf area of Cordia myxa L. and can be used without the easy way without the need to destroy the leaves.

Acknowledgements

CNPq, CAPES and FAPES for financial support.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Barroso, I.C.E., Oliveira, F. and Ciarelli, D.M. (2009) Morfologia da unidade de dispersão e germinação de Cordia sellowiana Cham e Cordia myxa L. Bragantia, 68, 241-249. https://doi.org/10.1590/S0006-87052009000100026
[2]	Favarin, J.L., Dourado Neto, D., García, A.G., Nova, N.A.V. and Favirin, M.G.G.V. (2002) Equações para estimativa do índice de área foliar do cafeeiro. Pesquisa agropecuária brasileira, 37, 769-773. https://doi.org/10.1590/S0100-204X2002000600005
[3]	Lopes, M.C., Andrade, I., Pedroso, V. and Martins, S. (2004) Modelos empírico para a estimativa de área foliar da videira na casta Jaen. Ciência e Tecnologia Vitivinicultura, 19, 61-75.
[4]	Benincasa, M.M.P. (1988) Análise do crescimento de plantas: Noções básicas. FUNEP, Jaboticabal, 42 p.
[5]	Schmildt, E.R., Schmildt, O., Alexandre, R.S., Fernandes, A.A. and Czepak, M.P. (2016) Modelos alométricos na determinação da área foliar de Bauhinia monandra Kurz. Comunicata Scientiae, 7, 415-419. https://doi.org/10.14295/cs.v7i3.1095
[6]	Oliveira, P.S., Silva, W., Costa, A.A.M., Schmildt, E.R. and Vitória, E.L. (2017) Leaf Area Estimation in Litchi by Means of Allometric Relationships. Revista Brasileira de Fruticultura, 39, 1-6. https://doi.org/10.1590/0100-29452017403
[7]	Cargnelutti Filho, A., Toebe, M., Burin, C., Neu, I.M.M. and Alves, B.M. (2018) Número de folhas para modelar a área foliar de mucuna cinza por dimensões foliares. Revista de Ciências Agroveterinárias, 17, 571-578. https://doi.org/10.5965/223811711732018571
[8]	Oliveira, V.S., Carvalho, C.F.M., França, J.M., Pinto, F.B., Santos, K.T.H., Santos, J.S.H., Santos, G.P., Pinheiro, A.P.B., Schmildt, O., Czepak, M.P., Arantes, S.D., Alexandre, R.S., Amaral, J.A.T., Vitória, E.L. and Schmildt, E.R. (2019a). Alometric Model for Estimation of Leaf Area of Garcinia brasiliensis Mart. through Non-Destructive Method. Journal of Agricultural Science, 11, 154-161. https://doi.org/10.5539/jas.v11n10p154
[9]	Oliveira, V.S., Pinheiro, A.P.B., Silva, J.V.G., Jardim, A.S., Fernandes, A.A., Schmildt, O., Arantes, S.D., Posse, R.P. and Schmildt, E.R. (2019) Uso de imagens digitalizadas para elaboração de equações que estime a área foliar de mudas de Eugenia uniflora L. International Journal of Development Research, 9, 31062-31064.
[10]	Pinheiro, A.P.B., Oliveira, V.S., Santos, K.T.H., Santos, J.S.H., Santos, G.P., Silva, J.V.G., Jardim, A.S., Longue, L.L., Nunes, S.F., Azeredo, A.L.R., Pires, F.R., Fernandes, A.A., Schmildt, O., Czepak, M.P. and Schmildt, E.R. (2019) Estimation Leaf Area by Composite Leaves of Canavalia rosea Seedlings through Linear Dimensions from Last Leaflet. Journal of Agricultural Science, 11, 299-308. https://doi.org/10.5539/jas.v11n9p299
[11]	Santos, J.S.H., Santos, K.T.H., Oliveira, V.S., Santos, G.P., Menezes, L.F.T., Czepak, M.P., Falqueto, A.R., Aoyama, E.M., Schmildt, O. and Schmildt, E.R. (2020) Regression Models for Prediction of Leaf Area in Purple Ipe [Tabebuia impetiginosa (Mart.)]. Australian Journal of Crop Science, 12, 654-659. https://doi.org/10.21475/ajcs.20.14.04.p2291
[12]	Alvares, C.A., Stape, J.L., Sentelhas, P.C., Gonçalves, J.L.M. and Sparovek, G. (2014) Köppen’s Climate Classification Map for Brazil. Meteorologische Zeitschrift, 22, 711-728. https://doi.org/10.1127/0941-2948/2013/0507
[13]	Willmott, C.J. (1981) On the Validation of Models. Physical Geography, 2, 184-194. https://doi.org/10.1080/02723646.1981.10642213
[14]	R Core Team (2020) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna.
[15]	Ferreira, E.B., Cavalcanti, P.P. and Nogueira, D.A. (2018) Package “ExpDes.pt”.
[16]	Pimentel-Gomes, F. (2009) Curso de estatística experimental. 15th Edition, Fealq, Piracicaba, 451 p.
[17]	Lavanhole, D.F., Oliveira, P.S., Vitória, E.L. and Aoyama, E.M. (2018) Estimativa de área foliar por meio de relações alométricas em Aechmea blanchetiana (Baker) L. B. SM sob distintas condições de luminosidade. Iheringia: Série Botanica, 73, 363-373. https://doi.org/10.21826/2446-8231201873313
[18]	Antunes, W.C., Pompelli, M.F., Carretero, D.M. and Damatta, F.M. (2008) Allometric Models for Non-Destructive Leaf Area Estimation in Coffee (Coffea arabica and Coffea canephora). Annals of Applied Biology, 153, 33-40. https://doi.org/10.1111/j.1744-7348.2008.00235.x
[19]	Montero, F.J., Juan, J.A., Cuesta, A. and Brasa, A. (2000) Non Destructive Methods to Estimate Leaf Area in Vitis vinifera L. HortScience, 35, 696-698. https://doi.org/10.21273/HORTSCI.35.4.696
[20]	Fascella, G., Darwich, S. and Rouphael, Y. (2013) Validation of a Leaf Area Prediction Model Proposed for Rose. Chilean Journal of Agricultural Research, 73, 73-76. https://doi.org/10.4067/S0718-58392013000100011