Next Article in Journal
Reduction and Reuse of Forestry and Agricultural Bio-Waste through Innovative Green Utilization Approaches: A Review
Previous Article in Journal
Effects of Heat Treatment on the Chemical Composition and Microstructure of Cupressus funebris Endl. Wood
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constructing and Validating Estimation Models for Individual-Tree Aboveground Biomass of Salix suchowensis in China

State Key Laboratory of Tree Genetics and Breeding, Co-Innovation Center for Sustainable Forestry in Southern China, Key Laboratory of Tree Genetics and Biotechnology of Educational Department of China, Nanjing Forestry University, Nanjing 210037, China
*
Author to whom correspondence should be addressed.
Submission received: 1 July 2024 / Revised: 30 July 2024 / Accepted: 5 August 2024 / Published: 6 August 2024
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)

Abstract

:
Biomass serves as a crucial indicator of plant productivity, and the development of biomass models has become an efficient way for estimating tree biomass production rapidly and accurately. This study aimed to develop a rapid and accurate model to estimate the individual aboveground biomass of Salix suchowensis. Growth parameters, including plant height (PH), ground diameter (GD), number of first branches (NFB), number of second branches (NSB) and aboveground fresh biomass weight (FW), were measured from 892 destructive sample trees. Correlation analysis indicated that GD had higher positive correlations with FW than PH, NFB and NSB. According to the biological features and field observations of S. suchowensis, the samples were classified into three categories: single-stemmed type, first-branched type and second-branched type. Based on the field measurement data, regression models were constructed separately between FW and each growth trait (PH, GD, NFB and NSB) using linear and nonlinear regression functions (linear, exponential and power). Then, multiple power regression and multiple linear regression were conducted to estimate the fresh biomass of three types of S. suchowensis. For the single-stemmed plant type, model M1 with GD as the single parameter had the highest adj R2, outperforming the other models. Among the 16 constructed biomass-estimating equations for the first-branched plant type, model M32 FW = 0.010371 × PH1.15862 × GD1.250581 × NFB0.190707 was found to have the best fit, with the highest coefficient of determination (adj R2 = 0.6627) and lowest Akaike Information Criterion (AIC = 5997.3081). When it comes to the second-branched plant type, the best-fitting equation was proved to be the multiple power model M43 with PH, GD, NFB and NSB as parameters, which had the highest adj R2 value and best-fitting effect. The stability and reliability of the models were confirmed by the F-test, repeated k-fold cross-validation and paired-sample t-tests. The models developed in this study could provide efficient tools for accurately estimating the total aboveground biomass for S. suchowensis at individual tree levels. The results of this study could also be useful for optimizing the economic productivity of shrub willow plantations.

1. Introduction

Global warming is a worldwide environmental crisis, with carbon dioxide identified as the primary contributor. The extensive use of fossil fuels for heat and electricity production, as well as for transportation, accounts for approximately 80% of global greenhouse gas emissions [1,2]. The escalating issues of global warming and climate change have garnered interest in the use of green biomass, particularly woody plants, due to their ecological and environmental benefits [3]. Biomass energy is green, renewable, abundant and widely distributed, playing an important role in alleviating energy shortages, optimizing energy structures and improving the ecological environment [4]. Consequently, the development and utilization of biomass energy is an inevitable strategy for achieving a low-carbon economy and sustainable development [5,6]. Among various renewable biomasses, forest biomass stands out due to its rich reserves and diverse utilization forms, making it an ideal alternative to fossil fuels [7]. Biomass directly reflects plant growth and serves as a crucial indicator for evaluating plant productivity. It also forms the basis of ecosystem matter and energy dynamics [8]. The timely and accurate estimation of forest biomass is essential for assessing biomass resources and carbon stocks, formulating biomass energy policies and conducting climate change studies [9].
Shrubs play a significant role in ecosystems and are vital components of forest resources, characterized by well-developed root systems, high biological productivity and wide ecological adaptability [10,11]. The well-developed root systems of shrubs and the abundance of inter-root microorganisms enhance the production, decomposition and accumulation of soil humus, which is the major driver of ecosystem productivity and diversity [12]. Although shrubs constitute a smaller proportion of forest ecosystems compared to trees, they are better-adapted to harsh climatic conditions and contribute significantly to environmental and ecological conservation in fragile areas such as deserts, plateaus and regions prone to rocky desertification. They prevent wind erosion, stabilize sand, conserve soil and water, and ameliorate saline and alkaline soils [13]. Despite their importance, research on shrub biomass modeling remains limited due to the smaller size and proportion of shrubs, the vast diversity of species and significant morphological and structural differentiation [14]. Previous studies have often overlooked the biomass and carbon stocks of shrubs and young forests, likely due to the absence of predictive equations for shrub biomass, which hindered the accurate quantification of forest biomass and its carbon stocks [9].
Salix suchowensis is a deciduous shrub of the genus Willow in the Populus family, with a wide distribution range, rapid growth and easy reproduction, which is not only an important urban and landscape greening tree species, but also an important energy and industrial tree species [15]. S. suchowensis grows fast and can reach the height of about 3 m, and efficiently converts solar energy into biomass [16]. Compared with many other annual energy crops, S. suchowensis has a short growing period and can regenerate from stumps after cutting, allowing for multiple harvests from the same plants, thus making it an excellent short-rotation coppice [17]. However, research on the biomass performance and growth model of S. suchowensis is scarce. There is a need to develop biomass estimation models for S. suchowensis across various ecological zones to support decision-making in natural resource management [18].
Traditional plant biomass measurements involve manually weighing the fresh weight of trees, a method that is both destructive and irreversible, and thus suitable only for small sample sizes [19]. In contrast, indirect modeling methods use simple, easily measurable indicators, such as diameter at breast height and plant height to establish biomass prediction models. This approach overcomes the limitations of the harvesting method and enables large-scale biomass assessments [20]. Currently, most studies, both domestic and international, estimate forest tree biomass and develop estimation models using related phenotypic parameters [21]. There are various mathematical models for predicting shrub biomass, including linear models, polynomial models and power exponential allometric growth models. Common predictive indicators include the diameter, height, crown width and number of branches, as well as composite factors like cross-sectional area (diameter squared) and plant volume (product of diameter squared and height). For instance, an allometric growth model for the aboveground biomass of four bamboo species was established using a power function with plant height, diameter at breast height and age as the input parameters [22]. Another study used diameter, tree height, wood density and crown width to develop a biomass estimation model [23]. Additionally, the allometric growth relationship between the dry biomass of the aboveground tree and diameter at breast height during individual tree development has also been reported [24].
In this study, we analyzed the variation in plant height, ground diameter, fresh weight biomass, number of first and second branches, and the ratio of the main stem to branches for S. suchowensis. The fresh weight biomass prediction models were established for S. suchwensis seedlings by performing linear, nonlinear and stepwise regression analyses using plant height, ground diameter, number of primary branches and number of secondary branches as the input parameters, and the F-test, repeated k-fold cross-validation and paired-sample T-test were applied on the model to test the generality of the model, thus providing a reference for the implementation of rapid biomass measurements of S. suchowensis.

2. Materials and Methods

2.1. Field Trial

An F1 population consisting of 892 S. suchowensis individuals was established and planted at the Baima Base of Nanjing Forestry University (31°36′46″ N, 119°9′16″ E) in Nanjing, China (Figure 1). This location is situated in the northern subtropical monsoon climate zone, characterized by yellow loam, an average annual temperature of 16.4 °C and an average annual precipitation of 1204.3 mm, with four distinct seasons (Table 1). Seeds obtained from the hybridization of SB48 × SB13 were sown in March 2023. Seedlings were transported to the Baima base for potting on 30 April 2023, and subsequently transplanted to the field on 20 May 2023, with a planting configuration of 0.5 m × 0.5 m. To minimize edge effects, protective rows were set up around the experimental plots.
On 12 March 2024, the plant height and ground diameter of S. suchowensis were measured using a height measurement pole (precision 0.1 cm) and a vernier caliper (precision 0.01 mm), respectively. Additionally, the number of first and second branches was investigated. The aboveground parts of the plants were then harvested, and the samples were immediately transported to the laboratory. The fresh weight of the whole plant, the main stem and branches were measured manually on the same day.

2.2. Data Processing and Analysis

Firstly, the obtained data were screened to remove abnormal data, such as dead or missing plants. Key parameters were measured, including plant height (PH, the distance from the plant shoot tip to the ground, in cm), ground diameter (GD, the diameter of the plant at 5 cm above the ground, in mm), number of first branches (NFB) and number of second branches (NSB).
Before data analysis and model construction, based on the biological characteristics of S. suchowensis and the results of the field investigation, we classified the samples into three categories according to the plant type (Figure 2), namely, single-stemmed type, containing only the main stem; first-branched type, consisting of the main stem and the first branches occurring directly from the main stem; and second-branched type, consisting of the main stem, the first-order branches and the second-order branches. Then, the four parameters, PH, GD, NFB and NSB, were selected as independent variables, and the fresh weight (FW) of S. suchowensis was selected as the dependent variable, and the R software was utilized to conduct the collinearity tests, correlation analyses, linear and nonlinear regression analyses, and cross-validation.
To construct the estimation models of S. suchowensis fresh weight biomass, it was essential to perform collinearity tests on the independent variables. Initially, the correlation coefficients between the independent variables were calculated, followed by the computation of the variance inflation factor (VIF). If the VIF value is over 10, the collinearity between independent variables is considered to exist; if the VIF value is less than 10, the collinearity between the independent variables can be ignored for subsequent analyses. The equation used was as follows:
V I F = 1 1 r 2
where r is the correlation coefficient between independent variables. Furthermore, the correlation between biomass and each factor was tested as well.

2.3. Model Fitting and Evaluation

We used the R software (v.4.4.0) to calculate the basic statistics for the phenotypic data and conduct the regression analysis, including linear regression, multiple nonlinear regression and multiple linear regression. The unitary regression models applied a linear function y = a + bx, the power function y = abx, and the exponential function y = aebx for a total of three types of functions. The multiple power models applied the power function y = ax1bx2cx3dx4e, where: y is the predictor variable (fresh weight biomass); x is the fitted variable (PH, GD, NFB, NSB); a is the regression constant; and b, c, d and e are the regression coefficients. The parameters of the models were estimated by the nonlinear least-squares (NLS) method using the “gslnls” package (v.1.3.2) [25], and the “MASS” package (v.7.3-61) [26] was used to perform multiple linear stepwise regression analyses on the following indicators: plant height, ground diameter, number of first branches and number of second branches, respectively. Box plots, line plots and scatter plots were produced using the “ggplot2” R package (v.3.5.1) [27].
To assess model performance, we selected the adjusted coefficient of determination (adj R2), Akaike Information Criterion (AIC), absolute mean residual (AMR) and root mean square error (RMSE). Higher adj R2 values and lower AIC, AMR and RMSE values generally indicate a better model fit. The specific equations used for evaluation are as follows:
a d j   R 2 = 1 i = 1 n   ( y i y ^ i ) 2 / ( n k 1 ) i = 1 n   ( y i y ¯ i ) 2 / ( n 1 )
AIC = n ln i = 1 n   ( y i y ^ i ) 2 n + n + n ln 2 π + 2 k
AMR = 1 n i = 1 n   y i y ^ i
R M S E = 1 n k i = 1 n   y i y ^ i 2  
where n is the number of samples; yi is the measured value of biomass; ŷi is the predicted value of sample i; ȳi is the mean of the measured value; and k is the number of parameters of the fitted equation.

2.4. Validation and Significance Tests

After fitting the biomass prediction models, we conducted the F-test, the repeat k-fold cross-validation and the paired-sample T-test on the regression models to accurately assess the generalization ability of the model and avoid overfitting and underfitting problems of the prediction models. The F-test helps determine whether the differences between models or data are significant or not, thus providing a statistical basis for model selection and data analysis. Repeated k-fold cross-validation involves splitting the data into multiple subsets (folds), with one subset used for model training and the remaining subsets for validation. This iterative process ensured that each data point is used for both training and validation, thereby enhancing the reliability and generalization of the models. The paired-sample T-test is able to determine whether the means under the two conditions are significantly different or not. Data processing and analysis for cross-validation were performed by the ‘caret’ R package (v.6.0-94) [28].

3. Results

3.1. Variations in the Growth Phenotypes

Initially, the total sample size for biomass calculations included 892 individuals. Following the exclusion of missing and non-viable data points, the refined dataset used for analysis comprised 885 individuals. The frequency distributions of directly measured phenotypic traits of S. suchowensis (FW, PH, GD, NFB) exhibited a well bell-shaped curve, except for NSB, which showed a skewed normal distribution (Figure 3). The statistics for biomass-related phenotypic data are presented in Table 2. Plant height and ground diameter are important parameters in describing plant growth and are integral to establishing biomass prediction models. The data revealed that PH had a mean value of 181.44 ± 17.74 cm, with a coefficient of variation (CV) of 9.78%, while GD averaged 13.24 ± 2.01 mm, with a CV of 15.19%. The fresh weight (FW) ranged from 22.25 g to 375.95 g, with a mean of 174.60 g and a CV of 38.50%, indicating considerable variability compared to PH and GD. Additionally, weights such as the weight of the main stem (WMS), total branches (WTB), first branches (WFB) and second branches (WSB) were also determined, serving as significant indicators of biomass accumulation and growth, all demonstrating large variation. The number of first branches (NFB) and second branches (NSB) exhibited large differences among individuals within the sample. Moreover, the fresh weight of the stem and branches are both important components of the aboveground fresh weight of the stand, accounting for 55% and 45% of the aboveground biomass, respectively, which are essential components of biomass prediction models.

3.2. Correlation and Multicollinearity Analyses

The construction of prediction models was based on the significant correlations between the dependent and independent variables. Four factors were selected to analyze the correlation with biomass. As is vividly shown in Figure 4, FW was significantly correlated with all independent variables. The correlation coefficients of GD and NFB with FW exceeded 0.5, which were 0.730 and 0.531, respectively. All of the factors can be used for biomass prediction model construction. PH exhibited a positive correlation with GD but a negative correlation with NFB and NSB. This suggested that an increase in ground diameter correlates with increased plant height, while an increase in branch number tended to suppress the height growth of S. suchowensis. NFB had a positive correlation with NSB, and GD was positively correlated with NFB and NSB.
To ensure the accuracy of the regression analysis, collinearity tests were conducted on the independent variables for multiple linear regression. The VIF values among the four independent variables are presented in Table 3. All of the VIF values were less than 10, ranging from 1.001 to 1.313, indicating that multicollinearity among the variables is negligible and does not compromise the predictive accuracy of the regression model. Thus, the multiple power regression and the multiple stepwise regression analyses could be carried out.

3.3. Construction of Biomass Models with a Single Variable

Using FW as the dependent variable and PH, GD, NFB and NSB as single independent variables, respectively, we chose three types of regression equations (linear, exponential and power functions) to establish the fresh weight biomass prediction models for S. suchowensis. After fitting the equations, we evaluated the goodness of fit using the adjusted R2, absolute mean residual (AMR), root mean square error (RMSE) and Akaike Information Criterion (AIC) of the fitted model. The model parameters and evaluation indicators are presented in Table 4 and Figure 5, Figure 6 and Figure 7.
For the single-stemmed plants, models using PH as the independent variable, the three functions yielded adj R2 values ranging from 0.4799 to 0.5034, with corresponding AIC values ranging from 151.5592 to 152.3449. For the GD biomass prediction models, the best fit was observed with the linear regression, achieving the highest adj R2 value of 0.7547 and the smallest AIC value of 139.5674, with an RMSE value of 11.5529. Moreover, the models constructed with GD were better than models constructed with PH, mainly due to higher correlation with FW. For the first-branched plants, it is always the regression model constructed by GD that had a better fit when compared to models established by PH and NFB. Of the nine functions for estimating first-branched plant biomass, the GD power model acts better than any other models, achieving an adj R2 value of 0.5423, slightly better than the 0.5373 and 0.5392 of linear and exponential functions. When it comes to the second-branched models with four single indicators, models constructed with GD as a single variable still outperformed other indexes, showing a better fit. Furthermore, among the three regression functions, the power function was more appropriate for the fresh weight biomass estimation for S. suchowensis.

3.4. Construction of Biomass Models with the Multiple Power Regression Model

Combining the biological growth pattern of woody plants and the results of the biomass equations constructed with different growth indicators as single parameters, we chose four growth indicators, namely, PH, GD, NFB and NSB, and used the power function equation as the basic regression function for constructing the multiple regression models of the aboveground fresh biomass of S. suchowensis. The regression results were presented in Table 5. The fitted model for single-stemmed plants comprised PH and GD, the adj R2 value of which was 0.7330 and the RMSE value was 11.2208. For the first-branched type, model M32 obtained an adj R2 value of 0.6627 and the RMSE value was 32.3746, which consisted of three variables, showing a better-fitting effect than models constructed by two variables. Furthermore, the established models with GD as parameters, including model M29 and M31, indicated a better-fitting effect than model M30. For second-branched type, the best-fitting effect was achieved by M43, which consisted of four variables and had the highest adj R2 and lowest RMSE value, at 0.6695 and 35.5435, respectively. The second fit model was M3,9 constructed by PH, GD and NFB. Compared to other two-variable models, M36 with GD and NFB as independent variables performed better.

3.5. Construction of the Biomass Model with the Multiple Linear Regression Model

The multiple linear regression models for S. suchowensis aboveground biomass are presented in Table 6. The constructed multiple linear model M44 obtained an adj R2 value of 0.7366, which was lower than model M1 with GD as the single variable. Except for the single-stemmed plant type, with the introduction of each independent variable, the adj R2 of the regression model gradually increased and the RMSE and AIC values gradually decreased, thus the model fit improved. The best model fitting effect for first-branched type was M46, which was established by GD, NFB and PH, and had the adj R2 value of 0.6329. For the second-branched plants, the optimal multiple linear regression model incorporated four variables GD, NFB, PH and NSB to construct an optimal stepwise regression model with an adj R2 value of 0.6592, an RMSE value of 36.0939 and an AIC value of 2589.6514.

3.6. Validation of the Optimal Biomass Prediction Model

As shown in Table 7 and Figure 8, the adj R2 values between the predicted and measured value of the regression models were all over 0.7. The slopes of the linear equations of the regression model for the single-stemmed type did not have a significant difference from the slopes of y = x, while the intercepts and slopes of the other equations were all significantly different, which indicated that there was a significant difference between the predicted and the actual values of the model. This indicated that the constructed models did not fit the data perfectly. Consequently, to better evaluate the predictive performance of the models, we conducted cross-validation of the selected models and paired-sample T-tests of the observed values against the predicted values of the models.
To further explore the predictive and practical performance of the models, repeated 10-fold cross-validation was conducted to assess the stability and reliability of the models. As can be seen from the cross-validation results in Table 8, the adj R2, AMR values and RMSE values obtained from cross-validation had little difference from the original model. Compared with the original model, the difference in adj R2 of model M1 was less than 0.0025 and the difference in RMSE was less than 1.4091; the difference in adj R2 of multiple power regression model M32 and model M43 was less than 0.0002 and 0.0008, respectively; the difference in the adj R2 of model M46 was less than 0.0002 and the difference in the RMSE was less than 0.2191; the difference in the adj R2 of model M49 was less than 0.0006 and the difference in the RMSE was less than 0.7036. These results proved the stability and predictive capability of the constructed models. The differences observed between the original and cross-validated metrics indicated that our models are reliable and well-suited for predicting the biomass of S. suchowensis. Furthermore, paired-sample T-tests were performed to determine whether there is a significant difference between the predicted values of the constructed models and the sample means of the observations. As can be seen in Figure 9, the results indicated no significant difference between the mean predicted and observed values of any of the nine regression models selected.

4. Discussion

The development and application of biomass energy can not only achieve the harmless, reduced and resourceful utilization of various organic wastes, thus improving the ecological environment, but also contribute to optimizing the energy structure, and realizing the synergy between carbon peaking, carbon neutrality and environmental pollution control [29]. Bioenergy derived from S. suchowensis biomass is more economically competitive than many other shrub species due to its sustainable high production, multiple harvest cycles and the high input–output ratio of energy exchange, making it an excellent short-rotation coppice (SRC) species [3,15,30,31]. However, the lack of reliable biomass estimation models hinders the utilization of biomass energy of S. suchowensis. The estimation of forest biomass basically relies on two methods: the direct measurement method and the indirect estimation method. The direct measurement method, known as the harvesting method, has the advantage of high measurement accuracy, but with the disadvantage of large effort, unrecoverable damage to the environment and high cost [32,33]. Consequently, indirect modeling methods are usually used to estimate biomass. Consequently, indirect modeling methods, particularly allometric growth models, are preferred for biomass estimation, having been applied across various tree species [21,34,35]. Sangsuree et al. employed allometric models to estimate the seedling biomass in the tropical dry evergreen forests of Thailand [36]. The biomass of Pinus radiata, Eucalyptus globulus and Eucalyptus nitens were estimated by the measured phenotypic data. Therefore, to maximize the S. suchowensis biomass, it is necessary to excavate its biomass model. The main objectives of this study were to accurately determine the fresh aboveground biomass of S. suchowensis in Nanjing, China, by traditional harvesting and weighing methods. Firstly, the samples were classified into three types by plant type and the branching conditions. Then, the phenotypic data were analyzed to assess the correlation and collinearity between biomass and variables to test the accuracy of biomass prediction by dynamic variables and equations, and to screen out and cross-validate species-specific prediction equations suitable for estimating the fresh weight biomass of S. suchowensis in subtropical regions of Nanjing.
GD and PH are important indicators for assessing tree biomass, reflecting radial growth and longitudinal growth, and demonstrating strong correlations with biomass accumulation. Many studies have selected PH, GD or the combined parameter (GD2 × PH) as independent variables to predict the biomass of different organs or the whole plant [37]. In this study, all parameters exhibited highly significant correlations with fresh weight biomass, which were suitable for fitting the biomass prediction model. Consistent with previous studies, the results of the single-variable model showed that GD is the most suitable measure for biomass prediction equations for S. suchowensis in Nanjing, demonstrating that S. suchowensis biomass production was mainly related to radial growth. Several studies have reported that the aboveground biomass estimation models constructed with GD alone or composite variables based on GD have good predictive abilities [30]. For instance, ground diameter, plant height and wood density were used to model the optimal biomass of 205 seedlings of 20 species in Thailand [36]. Zeng et al. found that the bivariate (D and H) model was slightly better than the univariate model in estimating single wood biomass of five conifers [21]. Similarly, the function y = aDbHc was found to be the most significant model for estimating the aboveground biomass of Larix sibirica in the forest steppe of Mongolia [38]. In line with previous studies, power function models constructed based on GD or parameters with GD included fitted well in this study. These findings collectively highlight the good predictive ability of GD and its derivative indexes in aboveground biomass prediction. There are several key factors contributing to this situation. (1) Morphologically, trees usually resemble cylinders or cones, whose volume is related to height and diameter, and the power function can well describe this heavy geometry. In shrubs, when the shrub growth reaches a certain height, its growth transforms into radial growth [39]. (2) The growth of woody plants follows a certain biological pattern that often manifests as a power relationship between various growth indexes [40]. (3) The growth process of an organism is often complex and multi-scaled, and the power function is sufficiently scale-invariant to capture this complexity succinctly [41].
However, compared to trees, shrubs have more complex morphological characteristics, with more clumped branches and lateral branches, which must be considered when establishing prediction models for the aboveground biomass of shrubs [42,43]. Accurate predictions of shrub biomass are challenging due to the abundant sprouts, branches and the dispersed characteristics [33,44,45]. In estimating the aboveground biomass of Rhododendron mariesii and Eurya muricata, incorporating crown area or the crown volume into the models increased prediction accuracy, which could reflect the morphological and structural characteristics of shrubs better than using D, H and D2H alone [46]. Therefore, selecting appropriate independent variables specific to particular species and regions is crucial for optimizing the shrub biomass prediction accuracy. During the investigation of field phenotypes, we found that S. suchowensis exhibits a higher number of lateral branches compared to other shrubs with pronounced main stem effects. Therefore, based on biological features and field observations, the samples were classified into three categories: single-stemmed type, first-branched type and second-branched type. Subsequently, the single-variable linear and nonlinear regression, as well as multiple power and multiple linear regression was conducted to figure out the effects of branches on aboveground biomass and refine the prediction model. Table 4 illustrates that the fitted models differ across various strain or branching types, even when the independent variables are the same. Furthermore, in the multiple regression analysis of first-branched and second-branched types, the inclusion of more variables generally increased the adj R2 value, indicating an improvement in the predictive ability of the models, which acted a similar function with the crown area [32]. These findings underscored the importance of considering lateral branches alongside traditional metrics when modeling shrub biomass, thereby improving the accuracy of biomass prediction models for S. suchowensis.
Indirect measurement methods for estimating biomass through constructing biomass equations based on easily measurable traits have been widely used to predict forest productivity and carbon sinks, with little environmental damage and great economic accuracy [37,38]. Forest biomass and carbon stocks could be effectively quantified using dynamic equations [35]. Many mathematical functions have been utilized in biomass estimation, such as the logarithmic, exponential, power and quadratic functions, and the most commonly used model is the power function (allometric growth function) [24]. Previous studies have demonstrated that the power function with GD2 × PH as the predictive variable was one of the best models to estimate the forest biomass [23,24]. In this study, three types of linear and nonlinear functions, including linear, exponential and power functions, as well as the multiple power and multiple linear regression were utilized to develop the predictive model for the aboveground biomass of S. suchowensis within a single growing period. The effectiveness of these mathematical equations varied across different parameters. Our findings indicated that the most suitable equations for fitting univariate parameters varied, where the linear equation with GD as the predictive indicator was most suitable for estimating the aboveground biomass of the single-stemmed plant type, while the power function was optimal for aboveground biomass estimation using GD, exhibiting a higher adj R2 value than other parameters, consistent with previous studies. Although single-variable function models have been suggested for estimating tree aboveground biomass [8,18], research has shown that the inclusion of more variables could enhance the accuracy of modeling shrub biomass. The predictive accuracy of the model improves as its complexity increases [46]. Zeng et al. estimated the biomass of organs and the total plants of four shrubs [11]. They found that the quadratic model with D2 (breast diameter) as the independent variable was the best predictive model for the total biomass of Loropetalum chinense and Gardenia jasminoides, while the best predictive quadratic model for white oak biomass involved V (D2H) as the independent variable. In our study, to improve the predictive ability of biomass models, we conducted multiple power regression and multiple linear regression on the aboveground biomass estimation of three plant types. The results showed that, except for the single-stemmed multi-parameter model, which performed similarly to the unitary model, the fitting effect of the other two plant types was superior to the other models with fewer parameters. This may be attributed to the complex biological characteristics of S. suchowensis.
To assess the applicability of the constructed biomass prediction models, we employed the F-test, repeated k-fold cross-validation and the paired-sample T-tests on the selected models. As obviously seen from Table 7 and Table 8, Figure 8 and Figure 9, although the significance test showed that there were significant differences between the regression function of predicted value and measured value and y = x function, which meant the predictive models could not perfectly predict the aboveground biomass, the paired-sample T-tests showed that there was little difference between the mean value of the predicted value and the measured value. Additionally, the cross-validation indicated that the selected models for three kinds of plant types of S. suchowensis exhibited good stability and predictive ability to estimate the aboveground biomass of S. suchowensis. Consequently, these models can be applied to make predictions for fresh biomass of the fresh aboveground portion of S. suchowensis. Furthermore, factors influencing plant biomass production are divided into environmental and genetic factors. Environmental factors include interregional climate and soil differences, while genetic factors contain photosynthetic efficiency, energy conversion efficiency and other endogenous factors [47]. The photosynthesis and growth of trees are strongly influenced by these factors, which in turn affects tree biomass, resulting in biomass modeling that is highly site-specific [48,49]. The models developed in this study are based on measurements and calculations of aboveground biomass of S. suchowensis within the study area; thus, validation is necessary before applying these models to regions with different terrains and climatic conditions [14].

5. Conclusions

Forests are the main carbon sequestration system around the world, and shrubs play a crucial role in ecological restoration due to their good fuel properties and carbon storage capacity. Numerous studies have proven the high biomass and carbon storage potential of shrubs. This study aimed to develop a non-destructive estimation model for predicting the aboveground biomass of S. suchowensis in subtropical China. The study identified correlations between the aboveground biomass and GD, PH, NFB and NSB. Additionally, S. suchowensis was classified into three plant types based on the morphological characteristics and the number of branches. By comparing the fitting results of linear regression, nonlinear regression and multiple power and multiple linear regression, the optimal model for estimating the aboveground biomass of S. suchowensis in three plant types was identified. For the single-stemmed plant type, the linear model M1 with GD as a single variable was the most appropriate. On the other hand, for the first-branched type, the multiple power model M32 was better for the biomass estimation, which included more indexes. Furthermore, the best biomass prediction model M49 for the second-branched type was established using multiple power regression analysis with PH, GD, NFB and NSB as variables. The stability and predictive accuracy of these models were verified through the F-test, cross-validation and paired-sample T-test. It is hoped that this study offers a reference for exploring shrub aboveground biomass in regions with similar climate and soil conditions to the study area. Future research should incorporate climate and environmental factors into the model construction to explore the variations in the prediction models across different regions, and expand the applicability and predictive capacity of the developed models.

Author Contributions

Y.C. conceived and designed the research; W.F., C.N., C.H. and P.Z. performed the experiments; W.F. analyzed the data and drafted the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (2021YFD2201100) and the National Natural Science Foundation of China (32071795).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Perera, F. Pollution from Fossil-Fuel Combustion Is the Leading Environmental Threat to Global Pediatric Health and Equity: Solutions Exist. Int. J. Environ. Res. Public Health 2018, 15, 16. [Google Scholar] [CrossRef] [PubMed]
  2. Shuit, S.H.; Tan, K.T.; Lee, K.T.; Kamaruddin, A.H. Oil Palm Biomass as a Sustainable Energy Source: A Malaysian Case Study. Energy 2009, 34, 1225–1235. [Google Scholar] [CrossRef]
  3. Kopp, R.F.; Abrahamson, L.P.; White, E.H.; Volk, T.A.; Nowak, C.A.; Fillhart, R.C. Willow Biomass Production during Ten Successive Annual Harvests. Biomass Bioenergy 2001, 20, 1–7. [Google Scholar] [CrossRef]
  4. Widagdo, F.R.A.; Li, F.; Zhang, L.; Dong, L. Aggregated Biomass Model Systems and Carbon Concentration Variations for Tree Carbon Quantification of Natural Mongolian Oak in Northeast China. Forests 2020, 11, 397. [Google Scholar] [CrossRef]
  5. Maud, H.; William, R.; Lauren, M.; Zhang, C.S.; Chang, S.J.; Michael, C.H.; Leslie, P.; Narender, N. Short-Rotation Woody Crops for Bioenergy and Biofuels Applications. In Vitro Cell. Dev. Biol. 2009, 45, 619–629. [Google Scholar] [CrossRef]
  6. Mleczek, M.; Rutkowski, P.; Rissmann, I.; Kaczmarek, Z.; Golinski, P.; Szentner, K.; Strażyńska, K.; Stachowiak, A. Biomass Productivity and Phytoremediation Potential of Salix alba and Salix viminalis. Biomass Bioenergy 2010, 34, 1410–1418. [Google Scholar] [CrossRef]
  7. Antar, M.; Lyu, D.; Nazari, M.; Shah, A.; Zhou, X.; Smith, D.L. Biomass for a Sustainable Bioeconomy: An Overview of World Biomass Production and Utilization. Renew. Sustain. Energy Rev. 2021, 139, 110691. [Google Scholar] [CrossRef]
  8. Dai, E.; Wu, Z.; Ge, Q.; Xi, W.; Wang, X. Predicting the Responses of Forest Distribution and Aboveground Biomass to Climate Change under RCP Scenarios in Southern China. Glob. Chang. Biol. 2016, 22, 3642–3661. [Google Scholar] [CrossRef] [PubMed]
  9. Cavanaugh, K.C.; Gosnell, J.S.; Davis, S.L.; Ahumada, J.; Boundja, P.; Clark, D.B.; Mugerwa, B.; Jansen, P.A.; O’Brien, T.G.; Rovero, F.; et al. Carbon Storage in Tropical Forests Correlates with Taxonomic Diversity and Functional Dominance on a Global Scale. Glob. Ecol. Biogeogr. 2014, 23, 563–573. [Google Scholar] [CrossRef]
  10. Treml, V.; Hejda, T.; Kašpar, J. Differences in Growth between Shrubs and Trees: How Does the Stature of Woody Plants Influence Their Ability to Thrive in Cold Regions? Agric. For. Meteorol. 2019, 271, 54–63. [Google Scholar] [CrossRef]
  11. Zeng, H.-Q.; Liu, Q.-J.; Feng, Z.-W.; Ma, Z.-Q. Biomass Equations for Four Shrub Species in Subtropical China. J. For. Res. 2010, 15, 83–90. [Google Scholar] [CrossRef]
  12. Nilsson, M.-C.; Wardle, D.A. Understory Vegetation as a Forest Ecosystem Driver: Evidence from the Northern Swedish Boreal Forest. Front. Ecol. Environ. 2005, 3, 421–428. [Google Scholar] [CrossRef]
  13. Xu, H.; Wang, Z.; Li, Y.; He, J.; Wu, X. Dynamic Growth Models for Caragana korshinskii Shrub Biomass in China. J. Environ. Manag. 2020, 269, 110675. [Google Scholar] [CrossRef] [PubMed]
  14. Huff, S.; Ritchie, M.; Temesgen, H. Allometric Equations for Estimating Aboveground Biomass for Common Shrubs in Northeastern California. For. Ecol. Manag. 2017, 398, 48–63. [Google Scholar] [CrossRef]
  15. Yang, G.; Li, X.; Xu, Q.; Ling, J.; Yin, T. Effects of Cutting Size on the Growth and Wood Property Traits of Short-Rotation Coppice Willows. Can. J. For. Res. 2020, 50, 1365–1372. [Google Scholar] [CrossRef]
  16. Jia, H.; Wang, L.; Li, J.; Sun, P.; Lu, M.; Hu, J. Physiological and Metabolic Responses of Salix sinopurpurea and Salix suchowensis to Drought Stress. Trees 2020, 34, 563–577. [Google Scholar] [CrossRef]
  17. Karp, A.; Shield, I. Bioenergy from Plants and the Sustainable Yield Challenge. New Phytol. 2008, 179, 15–32. [Google Scholar] [CrossRef] [PubMed]
  18. Mandal, G.; Joshi, S.P. Estimation of Above-Ground Biomass and Carbon Stock of an Invasive Woody Shrub in the Subtropical Deciduous Forests of Doon Valley, Western Himalaya, India. J. For. Res. 2015, 26, 291–305. [Google Scholar] [CrossRef]
  19. Salis, S.M.; Assis, M.A.; Mattos, P.P.; Pião, A.C.S. Estimating the Aboveground Biomass and Wood Volume of Savanna Woodlands in Brazil’s Pantanal Wetlands Based on Allometric Correlations. For. Ecol. Manag. 2006, 228, 61–68. [Google Scholar] [CrossRef]
  20. Güner, Ş.T.; Diamantopoulou, M.J.; Poudel, K.P.; Çömez, A.; Özçelik, R. Employing Artificial Neural Network for Effective Biomass Prediction: An Alternative Approach. Comput. Electron. Agric. 2022, 192, 106596. [Google Scholar] [CrossRef]
  21. Zeng, W.; Fu, L.; Xu, M.; Wang, X.; Chen, Z.; Yao, S. Developing Individual Tree-Based Models for Estimating Aboveground Biomass of Five Key Coniferous Species in China. J. For. Res. 2018, 29, 1251–1261. [Google Scholar] [CrossRef]
  22. Li, L.-E.; Yen, T.-M.; Lin, Y.-J. A Generalized Allometric Model for Predicting Aboveground Biomass across Various Bamboo Species. Biomass Bioenergy 2024, 184, 107215. [Google Scholar] [CrossRef]
  23. Djomo, A.N.; Chimi, C.D. Tree Allometric Equations for Estimation of above, below and Total Biomass in a Tropical Moist Forest: Case Study with Application to Remote Sensing. For. Ecol. Manag. 2017, 391, 184–193. [Google Scholar] [CrossRef]
  24. Zianis, D.; Mencuccini, M. On Simplifying Allometric Analyses of Forest Biomass. For. Ecol. Manag. 2004, 187, 311–332. [Google Scholar] [CrossRef]
  25. Chau, J.; gslnls: GSL Multi-Start Nonlinear Least-Squares Fitting. R Package Version 1.3.2. 2024. Available online: https://CRAN.R-project.org/package=gslnls (accessed on 20 June 2024).
  26. Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S, 4th ed.; Statistics and Computing; Springer: New York, NY, USA, 2002; pp. 107–250. [Google Scholar]
  27. Wickham, H. ggplot2: Elegant Graphics for Data Analysis, 2nd ed.; Springer: New York, NY, USA, 2016; pp. 75–186. [Google Scholar]
  28. Kuhn, M. Building Predictive Models in R Using the Caret Package. J. Stat. Softw. 2008, 28, 1–26. [Google Scholar] [CrossRef]
  29. Rozzi, E.; Minuto, F.D.; Lanzini, A.; Leone, P. Green Synthetic Fuels: Renewable Routes for the Conversion of Non-Fossil Feedstocks into Gaseous Fuels and Their End Uses. Energies 2020, 13, 420. [Google Scholar] [CrossRef]
  30. Börjesson, P.; Tufvesson, L.M. Agricultural Crop-Based Biofuels—Resource Efficiency and Environmental Performance Including Direct Land Use Changes. J. Clean. Prod. 2011, 19, 108–120. [Google Scholar] [CrossRef]
  31. Sassner, P.; Galbe, M.; Zacchi, G. Bioethanol Production Based on Simultaneous Saccharification and Fermentation of Steam-Pretreated Salix at High Dry-Matter Content. Enzym. Microb. Technol. 2006, 39, 756–762. [Google Scholar] [CrossRef]
  32. Gonzalez-Benecke, C.A.; Fernández, M.P.; Gayoso, J.; Pincheira, M.; Wightman, M.G. Using Tree Height, Crown Area and Stand-Level Parameters to Estimate Tree Diameter, Volume, and Biomass of Pinus Radiata, Eucalyptus Globulus and Eucalyptus Nitens. Forests 2022, 13, 2043. [Google Scholar] [CrossRef]
  33. McGinnis, T.W.; Shook, C.D.; Keeley, J.E. Estimating Aboveground Biomass for Broadleaf Woody Plants and Young Conifers in Sierra Nevada, California, Forests. West. J. Appl. For. 2010, 25, 203–209. [Google Scholar] [CrossRef]
  34. Yen, T.-M. Predicting Aboveground Biomass Yield for Moso Bamboo (Phyllostachys pubescens) Plantations Based on the Diameter Distribution Model. Eur. J. For. Res. 2023, 142, 1341–1351. [Google Scholar] [CrossRef]
  35. Chave, J.; Réjou-Méchain, M.; Búrquez, A.; Chidumayo, E.; Colgan, M.S.; Delitti, W.B.C.; Duque, A.; Eid, T.; Fearnside, P.M.; Goodman, R.C.; et al. Improved Allometric Models to Estimate the Aboveground Biomass of Tropical Trees. Glob. Chang. Biol. 2014, 20, 3177–3190. [Google Scholar] [CrossRef] [PubMed]
  36. Thippawan, S.; Chowtiwuttakorn, K.; Pongpattananurak, N.; Kraichak, E. Allometric Models to Estimate the Biomass of Tree Seedlings from Dry Evergreen Forest in Thailand. Forests 2023, 14, 725. [Google Scholar] [CrossRef]
  37. Roxburgh, S.H.; Paul, K.I.; Clifford, D.; England, J.R.; Raison, R.J. Guidelines for Constructing Allometric Models for the Prediction of Woody Biomass: How Many Individuals to Harvest? Ecosphere 2015, 6, art38. [Google Scholar] [CrossRef]
  38. Battulga, P.; Tsogtbaatar, J.; Dulamsuren, C.; Hauck, M. Equations for Estimating the Above-Ground Biomass of Larix Sibirica in the Forest-Steppe of Mongolia. J. For. Res. 2013, 24, 431–437. [Google Scholar] [CrossRef]
  39. Grime, J.P. Plant Strategies, Vegetation Processes, and Ecosystem Properties; Environmental Science, Biology; John Wiley & Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  40. West, G.B.; Brown, J.H.; Enquist, B.J. A General Model for the Structure and Allometry of Plant Vascular Systems. Nature 1999, 400, 664–667. [Google Scholar] [CrossRef]
  41. West, G.B.; Brown, J.H.; Enquist, B.J. A General Model for the Origin of Allometric Scaling Laws in Biology. Science 1997, 276, 122–126. [Google Scholar] [CrossRef]
  42. Paton, D.; Nuñez, J.; Bao, D.; Muñoz, A. Forage Biomass of 22 Shrub Species from Monfragüe Natural Park (SW Spain) Assessed by Log–Log Regression Models. J. Arid Environ. 2002, 52, 223–231. [Google Scholar] [CrossRef]
  43. Yuen, J.Q.; Fung, T.; Ziegler, A.D. Review of Allometric Equations for Major Land Covers in SE Asia: Uncertainty and Implications for above- and below-Ground Carbon Estimates. For. Ecol. Manag. 2016, 360, 323–340. [Google Scholar] [CrossRef]
  44. Chen, J.; Fang, X.; Wu, A.; Xiang, W.; Lei, P.; Ouyang, S. Allometric Equations for Estimating Biomass of Natural Shrubs and Young Trees of Subtropical Forests. New For. 2024, 55, 15–46. [Google Scholar] [CrossRef]
  45. Abich, A.; Mucheye, T.; Tebikew, M.; Gebremariam, Y.; Alemu, A. Species-Specific Allometric Equations for Improving Aboveground Biomass Estimates of Dry Deciduous Woodland Ecosystems. J. For. Res. 2019, 30, 1619–1632. [Google Scholar] [CrossRef]
  46. Qiu, R.C.; Miao, Y.L.; Zhang, M.; Li, H.; Sun, H. Modeling and verification of maize biomass based on linear regression anaysis. Trans. CSAE 2018, 34, 131–137. [Google Scholar] [CrossRef]
  47. Raich, J.W.; Russell, A.E.; Kitayama, K.; Parton, W.J.; Vitousek, P.M. Temperature Influences Carbon Accumulation in Moist Tropical Forests. Ecology 2006, 87, 76–87. [Google Scholar] [CrossRef] [PubMed]
  48. Riofrío, J.; Herrero, C.; Grijalva, J.; Bravo, F. Aboveground Tree Additive Biomass Models in Ecuadorian Highland Agroforestry Systems. Biomass Bioenergy 2015, 80, 252–259. [Google Scholar] [CrossRef]
  49. Dong, L.; Zhang, Y.; Zhang, Z.; Xie, L.; Li, F. Comparison of Tree Biomass Modeling Approaches for Larch (Larix Olgensis Henry) Trees in Northeast China. Forests 2020, 11, 202. [Google Scholar] [CrossRef]
Figure 1. Geographic location of the experimental and study area.
Figure 1. Geographic location of the experimental and study area.
Forests 15 01371 g001
Figure 2. Schematic diagrams depicting the classification of three distinct plant types among the sampled individuals.
Figure 2. Schematic diagrams depicting the classification of three distinct plant types among the sampled individuals.
Forests 15 01371 g002
Figure 3. Box plots and distribution curves for the five phenotypic traits. The colors of the y-axis correspond to those used in the images and legends.
Figure 3. Box plots and distribution curves for the five phenotypic traits. The colors of the y-axis correspond to those used in the images and legends.
Forests 15 01371 g003
Figure 4. Heat map of the correlation matrix for the five growth traits.
Figure 4. Heat map of the correlation matrix for the five growth traits.
Forests 15 01371 g004
Figure 5. Aboveground biomass regression plots for the type of single-stemmed Salix suchowensis. Blue lines are fitted lines; the blue and grey ribbons indicate the prediction of the interval and the confidence interval, respectively. (ac) Aboveground biomass predicted by GD; (df) aboveground biomass predicted by PH.
Figure 5. Aboveground biomass regression plots for the type of single-stemmed Salix suchowensis. Blue lines are fitted lines; the blue and grey ribbons indicate the prediction of the interval and the confidence interval, respectively. (ac) Aboveground biomass predicted by GD; (df) aboveground biomass predicted by PH.
Forests 15 01371 g005
Figure 6. Aboveground biomass regression plots for the first-branched Salix suchowensis. Blue lines are fitted lines; the blue and grey ribbons indicate the prediction interval and the confidence interval, respectively. (ac) Aboveground biomass predicted by GD; (df), aboveground biomass predicted by PH; (gi), aboveground biomass predicted by NFB.
Figure 6. Aboveground biomass regression plots for the first-branched Salix suchowensis. Blue lines are fitted lines; the blue and grey ribbons indicate the prediction interval and the confidence interval, respectively. (ac) Aboveground biomass predicted by GD; (df), aboveground biomass predicted by PH; (gi), aboveground biomass predicted by NFB.
Forests 15 01371 g006
Figure 7. Aboveground biomass regression plots for the second-branched Salix suchowensis. Blue lines are fitted lines; the blue and grey ribbons indicate the prediction interval and the confidence interval, respectively. (ac) Aboveground biomass predicted by GD; (df) aboveground biomass predicted by PH; (gi) aboveground biomass predicted by NFB; (jl) aboveground biomass predicted by NSB.
Figure 7. Aboveground biomass regression plots for the second-branched Salix suchowensis. Blue lines are fitted lines; the blue and grey ribbons indicate the prediction interval and the confidence interval, respectively. (ac) Aboveground biomass predicted by GD; (df) aboveground biomass predicted by PH; (gi) aboveground biomass predicted by NFB; (jl) aboveground biomass predicted by NSB.
Forests 15 01371 g007
Figure 8. Comparison between observed and predicted aboveground biomass of Salix suchowensis. Green lines are fitted lines, while dashed lines are reference lines 1:1. (ac) Single-stemmed type (models M1, M28 and M44); (df) first-branched type (models M10, M32 and M46); (gi) second-branched type (models M20, M43 and M49).
Figure 8. Comparison between observed and predicted aboveground biomass of Salix suchowensis. Green lines are fitted lines, while dashed lines are reference lines 1:1. (ac) Single-stemmed type (models M1, M28 and M44); (df) first-branched type (models M10, M32 and M46); (gi) second-branched type (models M20, M43 and M49).
Forests 15 01371 g008
Figure 9. Paired-sample T-test on the measured and predicted value. (ac) Single-stemmed type (models M1, M28 and M44); (df) first-branched type (models M10, M32 and M46); (gi) second-branched type (models M20, M43 and M49); ns, not significant.
Figure 9. Paired-sample T-test on the measured and predicted value. (ac) Single-stemmed type (models M1, M28 and M44); (df) first-branched type (models M10, M32 and M46); (gi) second-branched type (models M20, M43 and M49); ns, not significant.
Forests 15 01371 g009
Table 1. Environmental conditions of the experimental field.
Table 1. Environmental conditions of the experimental field.
SiteLongitude (E)Latitude (N)Annual
Temperature/°C
Min
Temperature/°C
Max
Temperature/°C
Annual
Rainfall/mm
Soil
Type
Plant Space/mTrial
Size
Nanjing119°9′16″31°36′46″16.4−10371204.3yellow loam0.5 × 0.5892
Table 2. Growth statistics of the Salix suchowensis F1 full-sib family.
Table 2. Growth statistics of the Salix suchowensis F1 full-sib family.
ParameterSample SizeMin.Max.RangeMeanSDCV/%
PH (cm)885103.00235.00132.00181.4417.749.78
GD (mm)5.9819.6213.6413.242.0115.19
FW (g)22.25375.95353.70157.8260.7638.50
WMS (g)11.55170.75159.2080.7723.6729.31
WTB (g)0.00301.65301.6577.0449.5164.27
WFB (g)0.00299.20299.2075.1748.3064.25
WSB (g)0.0049.6549.651.864.41236.82
NFB0.0023.0023.006.443.8760.16
NSB0.0011.0011.000.681.41206.77
WMS/FW0.171.000.830.550.1730.63
WTB/FW0.000.830.830.450.1737.85
Note: PH, plant height, in cm; GD, ground diameter, in mm; FW, fresh weight of the aboveground biomass, in g; WMS, weight of main stem, in g; WTB, weight of total branches, in g; WFB, weight of first branches, in g; WSB, weight of second branches, in g; NFB, number of first branches, in count; NSB, number of second branches, in count; WMS/FW, the ratio of WMS to FW; WTB/FW, the ratio of WTB to FW; SD, the standard deviation, in g; CV, the coefficient of variation.
Table 3. VIF values for the PH, GD, NFB and NSB of Salix suchowensis.
Table 3. VIF values for the PH, GD, NFB and NSB of Salix suchowensis.
ValueParameterPHGDNFBNSB
VIF valuePH-1.3131.0011.001
GD1.313-1.1341.033
NFB1.0011.134-1.055
NSB1.0011.0331.055-
Table 4. Construction and evaluation of the aboveground biomass model based on linear and nonlinear regression method with a single variable.
Table 4. Construction and evaluation of the aboveground biomass model based on linear and nonlinear regression method with a single variable.
Plant TypeNumber of ModelsVariableFunctionModeladj R2AMRRMSEAIC
Single-stemmedM1GD (mm)linear−89.254 + 14.994 × GD0.75478.275311.5529139.5674
M2PH (cm)linear−89.2289 + 0.9011 × PH0.479912.522816.8229152.3449
M3GD (mm)power0.5717 × GD2.03130.75207.967511.6163139.7537
M4PH (cm)power0.001256 × PH0.4747030.494011.968616.5939151.8788
M5GD (mm)exponential11.66273 × e0.16778×GD0.74158.188611.8608140.4619
M6PH (cm)exponential10.458999 × e0.010688×PH0.503411.790816.4386151.5592
First-branchedM7GD (mm)linear−118.0284 + 20.4825 × GD0.537330.338438.04256188.4420
M8PH (cm)linear−138.2263 + 1.5877 × PH0.246037.984948.56326486.8093
M9NFB (count)linear105.2827 + 7.3416 × NFB0.220339.361049.38426507.2940
M10GD (mm)power1.38538 × GD1.814980.542330.070937.83776181.8457
M11PH (cm)power0.00597 × PH0.156430.242138.113748.68856489.9570
M12NFB (count)power93.07986 × NFB0.285110.225439.487849.22326503.3047
M13GD (mm)exponential26.392691 × e0.130215×GD0.539230.331837.96646185.9965
M14PH (cm)exponential20.2828 × e0.01067×PH0.235638.277148.89706495.1793
M15NFB (count)exponential115.4 × e0.04135×NFB0.203139.673549.92566520.6188
Second-branchedM16GD (mm)linear−131.94 + 22.79 × GD0.464136.992245.79922703.0155
M17PH (cm)linear−126.9487 + 1.7016 × PH0.213742.803055.47642801.5444
M18NFB (count)linear116.6239 + 8.3883 × NFB0.273142.388253.34122781.3707
M19NSB (count)linear164.133 + 7.529 × NSB0.034549.545961.47622854.3284
M20GD (mm)power1.822 × GD1.75050.467236.748245.66762701.5364
M21PH (cm)power0.02423 × PH1.714610.212342.813055.52732802.0162
M22NFB (count)power91.59733 × NFB0.349080.271742.236553.39122781.8525
M23NSB (count)power168.49783 × NSB0.113720.034449.355661.48012854.3606
M24GD (mm)exponential33.425143 × e0.121141×GD0.462836.844745.85642703.6568
M25PH (cm)exponential32.738980 × e0.00938×PH0.208242.906755.67292803.3622
M26NFB (count)exponential131.2 × e0.0404×NFB0.260142.911653.81762785.9410
M27NSB (count)exponential166.05339 × e0.03769×NSB0.033549.585361.50892854.6016
Table 5. Construction and evaluation of the aboveground fresh biomass model based on the multiple power function.
Table 5. Construction and evaluation of the aboveground fresh biomass model based on the multiple power function.
Plant TypeNumber of ModelsVariableModeladj R2AMRRMSEAIC
Single-stemmedM28PH, GD0.65878 × PH−0.04224 × GD2.064030.73308.054111.2208141.7487
First-branchedM29PH, GD0.04184 × PH0.76849 × GD1.620380.566929.246136.74796149.1357
M30PH, NFB0.0013373 × PH2.143829 × NFB0.29075750.498430.520639.54536238.7903
M31GD, NFB1.85901 × GD1.59945 × NFB0.157140.609927.450334.87246085.1222
M32PH, GD, NFB0.010371 × PH1.15862 × GD1.250581 × NFB0.1907070.662725.540532.37465997.3081
Second-branchedM33PH, GD0.05619 × PH0.77712 × GD1.53530.502134.767843.97282685.1083
M34PH, NFB0.006992 × PH1.822947 × NFB0.3465520.508933.463243.67382681.6003
M35PH, NSB0.01443 × PH1.79633 × NFB0.139010.271640.513253.18802782.9024
M36GD, NFB2.20095 × GD1.49342 × NFB0.248350.603331.294739.25172626.7298
M37GD, NSB2.01224 × GD1.70305 × NSB0.039170.470436.464645.35332700.9963
M38NFB, NSB90.48203 × NFB0.33657 × NSB0.055570.279241.925652.90912780.1999
M39PH, GD, NFB0.01798 × PH1.07279 × GD1.18201 × NFB0.27010.668027.408135.76642581.9485
M40PH, GD, NSB0.04731 × PH0.84888 × GD1.44434 × NSB^0.057740.511034.200543.40582681.4509
M41PH, NFB, NSB0.005093 × PH1.879635 × NFB0.329125 × NSB0.0834160.529532.336642.57572671.5259
M42GD, NFB, NSB2.240458 × GD1.485723 × NFB0.246951 × NSB0.0079620.601931.282839.16352628.5872
M43PH, GD, NFB, NSB0.01635 × PH1.10962 × GD1.14101 × NFB0.26604 × NSB0.029360.669527.379535.54352581.7529
Note: PH, plant height, in cm; GD, ground diameter, in mm; NFB, number of first branches, in count; NSB, number of second branches, in count.
Table 6. The multiple linear regression models for aboveground biomass prediction.
Table 6. The multiple linear regression models for aboveground biomass prediction.
Plant TypeNumber of ModelsVariableModeladj R2AMRRMSEAIC
Single-stemmedM1GD−89.254 + 14.994 × GD0.75478.275311.5529139.5674
M44PH, GD−92.31803 + 14.42325 × GD + 0.05036 × PH0.73668.045311.1464141.5227
First-branchedM7GD−118.0284 + 20.4825 × GD0.537330.338438.04256188.4420
M45GD, NFB−111.1329 + 18.1047 × GD + 3.992 × NFB0.595528.124435.51266107.3536
M46GD, NFB, PH−200.504 + 14.22424 × GD + 4.74834 × NFB + 0.7471 × PH0.632926.902933.77586049.0851
Second-branchedM16GD−131.94 + 22.79 × GD0.464136.992245.79922703.0155
M47GD, NFB−135.1219 + 19.6094 × GD + 6.0419 × NFB0.599131.687339.45572629.3932
M48GD, NFB, PH−255.6368 + 15.2791 × GD + 6.62 × NFB + 0.9682 × PH0.653828.458636.52502592.7359
M49GD, NFB, PH, NSB−262.7903 + 14.6867 × GD + 6.4681 × NFB + 1.0191 × PH + 3.1043 × NSB0.659228.134236.09392589.6514
Note: PH, plant height, in cm; GD, ground diameter, in mm; NFB, number of first branches, in count; NSB, number of second branches, in count.
Table 7. Significance test for the intercept and coefficient of the linear regression equation against the y = x line.
Table 7. Significance test for the intercept and coefficient of the linear regression equation against the y = x line.
Plant TypeNumber of ModelsRegression
Equation
InterceptCoefficientadj R2
Single-stemmedM10.7854x + 17.7427p < 0.0010.10.8074
Single-stemmedM280.7777x + 18.4189p < 0.0010.10.8086
Single-stemmedM440.7860x + 17.6952p < 0.0010.10.8077
First-branchedM100.5443x + 68.2393p < 0.001p < 0.0010.7506
First-branchedM320.67136x + 49.10486p < 0.001p < 0.0010.7722
First-branchedM460.63528x + 54.62628p < 0.001p < 0.0010.7675
Second-branchedM200.47162x + 96.03041p < 0.001p < 0.0010.7486
Second-branchedM430.67697x + 58.68191p < 0.001p < 0.0010.7786
Second-branchedM490.66585x + 60.72916p < 0.001p < 0.0010.7758
Table 8. Cross-validation of the aboveground biomass prediction models for Salix suchowensis.
Table 8. Cross-validation of the aboveground biomass prediction models for Salix suchowensis.
Plant TypeNumber of ModelsFunctionVariableadj R2AMRRMSEAIC
Single-stemmedM1linearGD0.75828.289912.9620125.2333
Single-stemmedM28multiple powerPH, GD0.74098.098413.3771126.6291
Single-stemmedM44multiple linearPH, GD0.74338.106113.3424126.6930
First-branchedM10powerGD0.542330.071137.96195564.0003
First-branchedM32multiple powerPH, GD, NFB0.662925.526132.58335397.9424
First-branchedM46multiple linearPH, GD, NFB0.633126.893433.99495444.5853
Second-branchedM20powerGD0.467336.738446.02722431.7251
Second-branchedM43multiple powerPH, GD, NFB, NSB0.670327.338536.21972323.7875
Second-branchedM49multiple linearPH, GD, NFB, NSB0.659828.096536.79752331.1340
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fu, W.; Niu, C.; Hu, C.; Zhang, P.; Chen, Y. Constructing and Validating Estimation Models for Individual-Tree Aboveground Biomass of Salix suchowensis in China. Forests 2024, 15, 1371. https://doi.org/10.3390/f15081371

AMA Style

Fu W, Niu C, Hu C, Zhang P, Chen Y. Constructing and Validating Estimation Models for Individual-Tree Aboveground Biomass of Salix suchowensis in China. Forests. 2024; 15(8):1371. https://doi.org/10.3390/f15081371

Chicago/Turabian Style

Fu, Wei, Chaoyue Niu, Chuanjing Hu, Peiling Zhang, and Yingnan Chen. 2024. "Constructing and Validating Estimation Models for Individual-Tree Aboveground Biomass of Salix suchowensis in China" Forests 15, no. 8: 1371. https://doi.org/10.3390/f15081371

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop