Modeling Develops to Estimate Leaf Area and Leaf Biomass of Lagerstroemia speciosa in West Vanugach Reserve Forest of Bangladesh

Leaf area and leaf biomass have an important influence on the exchange of energy, light interception, carbon cycling, plant growth, and forest productivity. This study showed development and comparison of models for predicting leaf area and leaf biomass of Lagerstroemia speciosa on the basis of diameter at breast height and tree height as predictors. Data on tree parameters were collected randomly from 312 healthy, well-formed tree species that were considered specifically for full tree crowns. Twenty-four different forms of linear and powermodelswere compared in this study to select the bestmodel. Twomodels (M 10 andM 22 ) for the estimation of leaf area and leaf biomasswere selected based onR, adjustedR, rootmean squared error, corrected akaike information criterion, Bayesian information criterion and Furnival’s index, and the three assumptions of linear regression. The models were validated with a test data set having the same range of DBH and tree height of the sampled data set on the basis of linear regressionMorisita’s similarity index. So, the robustness of themodels suggests their further application for leaf area and biomass estimation ofL. speciosa in West Vanugach reserve forest of Bangladesh.


Introduction
Leaf area (A 1 ) and leaf biomass (B 1 ) estimation are significant basics of studying gas-exchange processes and modeling ecosystems.It is key traits in ecophysiological studies that determine assessing photosynthetic efficiency, evapotranspiration, atmospheric deposition, biogenic volatile organic emissions, light interception, and other ecosystem processes.Leaf area is valuable for the evaluation and understanding of individual tree growth models [1], biogeochemical models [2], and gap models [3].It is defined as the one-sided projected surface area which is an important consequence for the interception of radiant energy, the absorption of carbon dioxide, and the circulation of water between the foliage and the atmosphere [4].Leaf biomass constitutes one of the most important pools of essential nutrients, which is vital for forest nutrient cycling [5], including carbon cycling in a forest ecosystem.Leaf biomass estimates were considerably improved when additional biometric information relating to crown structure was added, whereas canopy B 1 is the product of leaf dry matter content and leaf area index [6].Interspecific variation in A 1 and B 1 has been connected with climatic variation, geology, altitude, or latitude, where heat stress, cold stress, drought stress, and high-radiation stress all tend to select relatively small leaves.Within climatic zones, leaf-size variation can also be linked to allometric factors (plant size, twig size, anatomy, and architecture) and ecological strategy [7].
Measurement of the destructive (direct) method of A 1 and B 1 is very time consuming, labour intensive, and ecounfriendly and depends on very small samples.But, nondestructive (indirect) methods were found user friendly and less expensive and can give accurate A 1 and B 1 estimates [8,9].To quantify the  1 and  1 of individual trees, the relationship between  1 and  1 with DBH (diameter at breast height) and height are widely used [10].This relationship differs between tree species and to some degree also within a tree species [11,12] and needs to be established for each species in question.For the development of models, the power function ( =   ) form is widely used in biological sciences particularly  for these relationships [13].But, in practice, the default nonlinear technique assumes homogeneity of errors that cannot be safely assumed with most model data.At present, models are developed by fitting a linear relationship between log-transformed diameter and leaf data.Although the logtransformed linear equation is mathematically equivalent to the power equation, they are not identical in the statistical sense [14].However, many model characteristics of organisms are multiplicative by nature and thus fitting models to logtransformed data is perfectly acceptable because accounting for proportional rather than absolute variation is most important [15].Lagerstroemia speciosa (L.) is the dominant tree both spatially and vertically in the forests of Bangladesh.It is well distributed in natural forests as well as in plantation forests [16].It is also found at low to medium altitudes in comparatively open habitats, in disturbed or secondary forest and grassland.The habitat may vary from well drained to occasionally flooded but not peat soil and it is resistant to fire.L. speciosa is the most important wood species in parts of South and Southeast Asia [16], but to estimate A 1 and B 1 relationships of this species is not available.So, the main goal of my study was to develop and compare models of L. speciosa that will help to estimate A 1 and B 1 of the species in reserve forest of Bangladesh.1).The soil of this forest is alluvial brown sandy clay loam to clay loam.Average annual rainfall is 4,332 mm and January is the coldest month (minimum temperature around 11 ∘ C), while May and October are the hottest months (average maximum temperature around 33 ∘ C).The relative humidity is about 75% during December and over 92% during July-August [17].This forest stand is a mixed forest stand and its vegetation is recognized as old growth plantation forest.It includes 78 species of trees, 14 species of shrubs, 42 species of herbs, and 25 species of climbers [18].

Species Description. Lagerstroemia speciosa (L.
) is a medium-sized tree growing to 20 meters tall, much branched deciduous tree, with smooth and flaky bark [16].Its leaves are opposite, elliptic, or oblong-lanceolate, 8.5-20.0cm long and 3.0-7.5 cm broad, acuminate at the apex, acute to rounded at the base, and glabrous and finely reticulate on both surfaces, have lateral nerves 8-12 on either half, prominent beneath, and petioles 5-9 mm long, and are stout.Its wood is strong and fairly durable and is used for carts, furniture, and house posts.The root is considered as an astringent, stimulant, and febrifuge and due to the tree's dense and wide spreading root system, it is used in erosion control.Its extracts have been used as traditional medicines and are effective in controlling diabetes and obesity.A leaf poultice is used to relief malarial fever and is applied on cracked feet.It is planted in forests, edges of forest streams, swamps, ditches, and river banks.It is growing in South East Asia, Bangladesh, India, and Philippines and it is also widely cultivated as an ornamental plant in tropical and subtropical areas [19].

Field Measurements.
Leaf, DBH, and tree height data were collected for three hundred twelve sampled healthy wellformed L. speciosa tree species with full crown by maintaining the protocols for standardized and easy measurement of plant functional traits worldwide developed by Cornelissen et al. 2003.Field sampling occurred in the rainy season during which A 1 and B 1 were maximum.Diameter at breast height (DBH; at 1.3 m height from ground level) and tree height () was measured with Suunto clinometers (Suunto PM5-360, Vantaa, Finland) and tree calipers (Haglöf Graduated Aluminum Tree Caliper, Ben Meadows, Janesville, Wisconsin, USA).Data on tree parameters are usually measured and collected after harvesting for improving data accuracy [6,20].Though using nondestructive methods is a highly reliable option [21] where whole tree removal is not possible, as is the case in Bangladeshi reserve forests [22], leaf data were collected manually (by climbing the tree).To retain high data accuracy, stratified random sampling was used for calculation of the leaf number per individual tree.Thereafter, the following steps were performed to count the leaves of each sampled tree: (I) the diameter of all of the main branches was measured; (II) the diameter of the main branch closest to the mean diameter was selected as the model main branch (MMB); (III) subbranches of all main branches were counted; (IV) three subbranches were selected randomly from the MMB; (V) twigs of each subbranch were counted and the mean twig number per subbranch was calculated; (VI) three twigs were selected randomly from each selected subbranch; (VII) the average number of leaves per twig was calculated; (VIII) the leaf number per subbranch was calculated by multiplying the mean twig number per subbranch and mean leaf number per twig; and (IX) the total leaf number per tree was calculated by multiplying the total number of subbranches and estimated leaf number per subbranch [9].The leaves were categorized into small, medium, and large size and fifteen leaves of each category were collected according to different directions and heights of a tree.A total of fortyfive samples of leaves were taken from each sampled tree and packed in a plastic bag.The leaf area of each leaf sample was measured with a leaf area meter (CI-202, CID, Inc., Vancouver, Washington, USA) for small and medium leaves and Adobe Photoshop CS5 software for large leaves.Data of DBH, height, leaf samples, and total leaf number per tree from additional 40 trees were collected and used as a test data set for validating the models.The projected leaf area for each sample tree (A 1 , m 2 ) was calculated by multiplying the average leaf area and total leaf number.Thickness of each sample leaf was estimated with digital caliper (Absolute Digimatic CD-6  CS, Mitutoyo Corporation, Kanagawa, Japan) and the fresh mass of each leaf was measured with a balance meter.Measured leaves were oven-dried at 65 ∘ C for 72 h and weighed to determine fresh mass: dry mass ratios.Total leaf volume per tree was estimated by multiplying total leaf area and leaf thickness.Leaf density was calculated by dividing by the average leaf fresh mass and the leaf volume of forty-five sample leaf in a tree.Total leaf biomass ( 1 , kg) per tree was estimated by multiplying the total leaf volume and leaf density.

Statistical Analysis.
Regression analysis was used to evaluate the appropriateness of each of the two independent variables (): height (m) and DBH (cm), as estimators of total  1 and  1 (dependent variables = y) for L. speciosa by using  statistical software version 2.15.2 and SPSS 17 software (SPSS, ver.17, Chicago, IL, USA).The data with two regression models (linear and power) was analyzed.At first, scatter plots were used to see whether the relationship between independent and dependent variables was linear.Secondly, given the presence of heteroscedasticity, power models transformed the data for linear regression using natural logarithmic.A total of twenty-four models were developed using DBH and height for the best model selection (Table 2).The performance of the developed models was evaluated by examining the goodness of fit ( 2 and adjusted  2 ), root mean squared error (RMSE), corrected akaike information criterion (AICc), Bayesian information criterion (BIC), and Furnival's index (FI).Thus, to select between models, the corrected akaike information criterion (AICc) for each model was then computed as where  is the sample size,  is the number of parameters in the statistical model, and  is the corresponding likelihood [23].The formula for the BIC (Bayesian information criterion) is where σ2  is the error variance [23].For models with different dependent variables using  2 to compare them could sometimes produce false results [24].So, Furnival's index [25] which is able to compare models of different dependent variables or weights was used to compare the logarithmic and nonlogarithmic models.The index is calculated as follows: where   () is the derivative of the dependent variable with respect to  1 and  1 , MSE is the mean square error of the fitted model, and the square bracket ([]) is the geometric mean.After that, the three regression assumptions were used to judge the consistency of the models for the residues: (I) no outlier, (II) the distribution of these residuals should be normal with mean = 0 and a constant variance and (III) the data points must be independent [26].The Durbin-Watson test was used to check for autocorrelation.Though logarithmic transformation is reported to increase the statistical validity of regression analysis by homogenising variance, it introduces a slight downward bias when data are back-transformed to arithmetic units [27].To account for the bias, the back-transformed results from logarithmic unit are usually multiplied by a correction factor.Consequently, a correction factor (CF) was calculated for all logarithmic models.The CF is given by the following [28]: where SEE is the standard error of the estimate.
For validation of the models, the observed  1 and  1 data regressed from a test data set ( = 40) against the predicted  1 and  1 using linear regression.As a measure of distance between observed and predicted values, Morisita's [29] dissimilarity index (  ) was calculated as follows: where   = observed  1 and  1 of the sampled data and    = predicted  1 and  1 .Note that 1-  is essentially Morisita's measure of niche overlap or similarity [30].

Model Development.
DBH and tree height of sampled individuals of L. speciosa ranged from 7.52 to 65.16 cm and 3.27 to 19.01 m, respectively, and the  1 and  1 of sampled trees ranged from 36.20 to 141.31 cm 2 and 0.35 to 1.26 g, respectively.In addition, the range of DBH and H data in the model test data set fell within the boundary of the sampled data set (Table 1).Table 3 summarizes the different models derived to predict  1 and  1 with DBH and H as the independent variable.For the development of the best model, 24 different model forms were determined (Table 2).Model M 10 and M 22 were found to be the best predictor for  1 and  1 estimation, respectively.Whereas the RMSE value was the lowest and the goodness of fit ( 2 and adjusted  2 ) was highly significant and explained more than 96% variation for the selected models (M 10 and M 22 ).To compare the models, the AICc and BIC values were calculated and the results are presented in Table 3.The AICc and BIC values for models M 10 and M 22 were lower than that of other tested models indicating statistical robustness of the selected models.Based on Furnival's index (FI) of fit, models M 10 and M 22 performed better in both  1 and  1 estimation, whereas the FI values for both models were lower (0.36 and 0.47, respectively) than that of other tested models.
To avoid any confusion, the three regression assumptions were checked to select the best model.A good model should conform to all valid statistical assumptions [26].First assumption, residual analysis of models M 10 and M 2 , showed that the values of standard residual were 2.93 and 2.48, respectively.A model to be outlier free requires that the value of standard residual must be below 3 [26].In both models this requirement was fulfilled.The normality of residuals was measured for selected models by histogram and normal  probability plot.The distribution of the residual satisfies the normality assumptions for both models M 10 (Figures 2(a Second assumption, diagnosis for constant variance, was done by plotting the regression standardized residuals and predicted values in a scatter plot.It was shown that the scatter of the residual points was randomly distributed (Figures 2(c) and 3(c)).Thus it can be concluded that the variance was constant for both models.Third assumption, the Durbin-Watson statistic, lies in the range 0.0-4.0.A value of 2 or nearly 2 indicates that there is independence of the data points [31].This range of the selected models was 2.003 for  1 and 1.971 for  1 .Thus, the independence assumption was satisfied for both models.However, due to the systematic bias introduced by logarithmic transformation of data, correction factor is normally calculated for logarithmic transformed models so as to account for the bias.Correction factors showed a rather narrow variation for  1 and  1 (Table 3).

Model Validation.
In the model validation, the goodness of fit ( 2 ) showed that there was a highly reliable relationship between predicted and observed data for  1 and  1 .The  2 values between them for models using  1 or  1 were 0.9716 and 0.9736, respectively (Figure 4).Morisita's [29] index of similarity (1 −   ) was 0.994836 for  1 and 0.993567 for  1 (almost 1), which indicated that the models can predict  1 and  1 , outside the sampled data set, effectively.

Discussion
Different combinations and forms of model of diameter and tree height were used to select the best models (M 10 and M 22 ) that are more than 96% variation and lower RMSE, AICc, BIC, FI value in measuring  1 and  1 of L. speciosa.After that, logarithmic transformation induces a systematic bias in the estimation, which was corrected using a CF in the final model [27].Given the paucity of models for this forest, I expect these models to prove valuable in future research that requires estimates of  1 and  1 of this species.Burton et al. [32] found that DBH is the best predictor for estimating  1 and  1 ( 2 > 0.90).Using a nondestructive sampling technique for  1 estimation, Grace [33] found that DBH could explain 91% variation in A 1 of Acacia Koa.Although sapwood area (AS) or sapwood volume has proved a good predictor of tree canopy properties in different temperate in the past [34,35], it has not been measured in this study or any national forest inventories in Bangladesh.However, Turner et al. [36] found that, in Pseudotsuga menziesii, the A 1 estimates based on DBH were more accurate than those on AS.So, I excluded AS from this study.Other studies found higher [37][38][39][40][41]  also show that allometric models are explained for more than 95% of the variation based on DBH with A 1 and B 1 .L. speciosa is the dominant and medium sized deciduous tree species in tropical semievergreen forest of Bangladesh.Thus, this study is important to tree growth models of this species and ecological purposes, including light interception, transpiration rate, biomass estimation, and carbon storage.As a result, I suggest estimating  1 and  1 of this species based on available standard inventory measurements such as DBH and tree height in Bangladesh.

Conclusion
The consequences of this result showed a strong statistical dependence and best fitted models (M 10 and M 22 ) between A 1 and B 1 with DBH and H.These models can be used effectively at the stand level for estimation of  1 and  1 of L. speciosa quickly, accurately, and nondestructively.Therefore, I propose that the models obtained and validated in this study could be employed with confidence by foresters, forest ecologists, and other scientists for estimating leaf area and leaf biomass in Bangladeshi reserve forests.

Figure 1 :
Figure 1: Study site in West Vanugach reserve forest of Bangladesh.

Figure 2 :
Figure 2: Residual distribution (a), normal P-P plot of regression standardized residuals (b), and scatter plot of residuals (c) for model M 10 .

Figure 3 :
Figure 3: Residual distribution (a), normal P-P plot of regression standardized residuals (b), and scatter plot of residuals (c) for model M 22 .

Figure 4 :
Figure 4: Linear regression between observed and predicted leaf area (a); between observed and predicted leaf biomass (b).

Table 1 :
General characteristics of L. speciosa species.

Table 3 :
Different models tested for the estimation of leaf area ( 1 , m 2 ) and leaf biomass ( 1 , kg) with diameter at breast height (DBH, cm) and tree height (, m) of L. speciosa species.: root mean squared error, FI: Furnival's index, AICc: corrected akaike information criterion, BIC: Bayesian information criterion, and CF: correction factor. RMSE