Next Article in Journal
Project Implementation Management Modalities and Their Implications on Sustainability of Water Services in Rural Areas in Ethiopia: Are Community-Managed Projects More Effective?
Previous Article in Journal
Executive Compensation, Sustainability, Climate, Environmental Concerns, and Company Financial Performance: Evidence from Indonesian Commercial Banks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Above-Ground Biomass Models of Caragana korshinskii and Sophora viciifolia in the Loess Plateau, China

1
State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Chinese Academy of Sciences and Ministry of Water Resources, Yangling 712100, Shaanxi, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
College of Resource and Environment, Northwest Agriculture & Forest University, Yangling 712100, Shaanxi, China
*
Author to whom correspondence should be addressed.
Sustainability 2019, 11(6), 1674; https://0-doi-org.brum.beds.ac.uk/10.3390/su11061674
Submission received: 14 February 2019 / Revised: 15 March 2019 / Accepted: 16 March 2019 / Published: 20 March 2019

Abstract

:
The quantification of above-ground biomass is based on the calculation of carbon storage, which is important for the balance of carbon cycling. However, the allometric models of shrubs for calculating the above-ground biomass of shrubs in the Loess Plateau are scarce. In order to solve this issue, this study analyzed some highly correlated variables, including height (H), branch diameters (D), canopy volume (Cv), canopy area (Ca), and then built a regression model to predict the above-ground biomass in two common shrubs (Caragana korshinskii and Sophora viciifolia) in the Loess Plateau, China. The results show that the above-ground biomass of these two shrubs can be accurately predicted by H and D, and then we can use allometric model (y = axb) to calculate shrub above-ground biomass (including leaf biomass and branch biomass). Furthermore, the correlation between leaf biomass and branch biomass in Caragana korshinskii and Sophora viciifolia indicates that the components of above-ground biomass are closely related to each other. In addition, there is a strong linear relationship (p < 0.01) between the observed and estimated biomass values, which confirms the data accuracy of the above-ground biomass estimation models. In summary, these two biomass estimation models provide an accurate way to estimate the quantification of carbon for shrubs in the Loess Plateau.

1. Introduction

Above-ground biomass, the amount of organic matter in dead and living material, is a key plant property resulting in the balance of carbon cycling [1,2,3,4,5]. The shrub biomass model is an important method for estimating shrub biomass. In the past decades, above-ground biomass has received more attention in terms of carbon storage in ecosystems [4,5,6]. Accordingly, there are a multitude of methods for above-ground biomass estimation by harvesting individuals with morphometric variables through an allometric model (including trees and shrubs) [4,7,8]. Furthermore, advanced techniques and methods are widely used to determine above-ground biomass [9,10]. For example, height (H), branch diameters (D), crown volume (Cv), and crown cover (C) may be employed in biomass estimation models as the independent variables [11,12,13,14,15]. As for trees, stem-related variables, specifically height (H) and diameter (D), have proven to accurately predict individual biomass [4,10,14]. In contrast, shrub biomass is an important component for the total biomass in a terrestrial ecosystem. In multi-stemmed shrub biomass estimations, each stem is usually considered as a separate shrub, thus, by applying the standard relationship between D and H, one can estimate the shrub biomass [16,17,18]. Furthermore, many scientists estimate shrub biomass using the same formula as with tree biomass estimation, which is confirmed in the literature [4,10,14]. However, there is little literature to report the above-ground biomass estimation for shrubs because of the lack of methodology and difficulty in calculation.
In the Loess Plateau, above-ground biomass has greatly increased because of the Grain for Green Project from 1999, and much more attention has been placed on developing generalized allometric models for shrubs [15,19,20]. In particular, there are several common shrub species in this region where allometric models have not yet been developed. Therefore, developing more accurate regression equations for shrub biomass estimation is essential for estimating carbon storage in this region. In this way, the choice of fitting methods and models is the most important issue to be solved [9,10]. Previous studies have reported that the biomass equations can be separated into additive and non-additive equations [21,22,23]. In additive equations, the estimation for the biomass components is calculated to predict the total biomass [4,14]. In contrast, in non-additive equations, the biomass components are neglected, and the biomass components are measured separately [23,24]. Therefore, the calculation of shrub biomass (n = 57) should be based on the same method and variables, the most appropriate regression equation, and regression coefficient test for biomass components (n = 58). This study aimed to find a scientific method for quantifying shrub biomass, which is easy to use in the field, and produces reliable results for forecasting carbon storage in the Loess Plateau.
In this study, two common shrubs (Caragana korshinskii and Sophora viciifolia) were investigated in the Loess Plateau, and some ecological variables were tested to predict shrub biomass. The aim was: (1) to develop equations for biomass components based on easily measured variables; (2) to calculate parameters for estimating biomass components by fitting regression models; and (3) to test the constructed models in the field. Although this study included only two shrub species, it can provide a basis for the standardization of useful variables to predict shrub biomass in the Loess Plateau.

2. Materials and Methods

2.1. Study Site

The study area, Yanhe River catchment (36°23′–37°17′ N, 108°45′–110°28′ E), north of Shaanxi Province, is located in the middle of the Yellow River catchment. Hills cover 90% of the region (which is 7687 km2 in total area). The Yanhe River catchment has a periodic flooding and heavy seasonal rainfall with a semi-arid climate. The average annual rainfall from 1970–2000 was about 497 mm, and there are distinct rainy and dry seasons. The average annual temperature is between 7 °C to 9 °C. The rainy season is from July to October, and August accounts for 23% of the total annual rainfall. Most of this area lies at an altitude between 900 to 1500 m. The topography of this area is characterized by cliffs with steep slopes (about 30%–40%). According to the Chinese Soil Taxonomy (Chinese Soil Taxonomy Research Group, 2001), the main soil type is typical Haplic-ustic Cambisol in this area. The dominant vegetation types are forest land, shrubland, natural grassland, and cropland, and the most abundant plants are medium or small shrubs, including Caragana korshinskii (non-native shrub artificial planted shrub) and Sophora viciifolia (native shrub) [25,26].

2.2. Sampling Design

The field surveys were conducted in the middle of August 2016 in the Loess Plateau, and two representative multi-branched shrubs, Caragana korshinskii and Sophora viciifolia, were selected across the whole study area (Figure 1). Keeping the environmental and climatic conditions as constant as possible, six sites were selected in this study (Table 1) and each site was 50 × 50 m2. There were five quadrats in each plot, as shown in the Figure 1, and the distance among the sample plots was about 1 km. During the field sampling, each quadrat was set as 5 × 5 m2. Finally, more than 800 leaves and branches of these two shrubs were harvested. More than 100 individual shrubs were chosen by grid method to cover represent each species in the field. The two common shrubs in the Loess Plateau were selected Each branch was considered as being the main upright portion of the shrubs [27].
Individual shrubs were harvested at ground level, and fresh subsamples of leaves and branches were cut down and stored in sealed plastic bags. Then, the samples were oven-dried in forced-air ovens at 80 °C until constant weight. The water weight was subtracted from the individual fresh mass to get the leaf and branch biomass (Table 1).
According to above-ground biomass, several variables were quantified to be used to estimate leaf and branch biomass: the maximum height (Hmax), defined as the distance between the highest crown point and the ground surface; number of branches emerging from the root; the maximum (D1) and minimum (D2) canopy diameters were measured in 12 transects of 100 m × 100 m. With these easy and measurable variables, canopy volume (Cv) and canopy area (Ca) were calculated as [28]:
C a = π 4 ( D 1 D 2 ) C F = e ( S E E × S E E / 2 )
C v = C a H max

2.3. Biomass Equations

Regression analyses were used to examine the relationship among above-ground biomass components and the measured variables, including diameter (D), height (H), canopy volume (Cv), and canopy area (Ca). A graphical and statistical analysis indicated that D, D2H explained most of the variation in biomass, and the optimum biomass components for the regression equation model was calculated by least-squares procedures using D and D2H. Compared with regression coefficients, the following equations were presented for describing the relationship between above-ground biomass components:
W = a + b x + c x 2
W = a x 2
W = a + b ln x
In addition, the optimum statistical model was selected according to the Akaike Information Criterion (AIC) which penalizes the number of parameters [29,30]. To reflect the accuracy of the statistical model, the predictive mean squared error (MSE) was chosen to be the alternative statistic. Besides, the total variance of equations was conventionally determined by the summation of squared residuals. Thus, the standard error of the estimation (SEE), the coefficient of determination (R2), and F value were used to make the goodness fitting. The equations were originally described in previous reports. The parameter as follows [31,32]:
R 2 = 1 R S S T S S
y 2 = [ i = 1 n y 1 ] / n
R S S = i = 1 n ( y 1 y 2 ) 2
T S S = i = 1 n ( y 1 y 2 ) 2
F = T S S R S S R S S / ( n 2 )
where y1 refers to the estimated biomass (g·m−2); yi refers to the observed biomass (g·m−2); y2 refers to the mean of the observed value of yi (g·m−2); and n refers to the number of observations.
Finally, in order to examine the constructed model in the field, correction factors (CF) were computed so as to be applied to above-ground biomass in the constructed model, and CF were calculated as follows [31,32]:
a b o v e g r o u n d b i o m a s s = C F × e ln ( a b o v e g r o u n d b i o m a s s )
C F = e ( S E E × S E E / 2 )
where SEE refers to the standard error of the estimation; CF refers to correction factor.

2.4. Statistical Analysis

SPSS 20.0 (SPSS Inc., Chicago, IL, USA) was used for all analyses. Before the parametric tests, a homogeneity of variance test was completed. Analysis of variance (ANOVA) was used to analyze the differences between the Caragana korshinskii and Sophora viciifolia; multiple comparisons were used to determine significant analysis by Tukey’s test (p < 0.01 or p < 0.05). Average and standard deviation were determined to fit to a Gaussian distribution using the Kolmogorov–Smirnov test. In Gaussian distribution, the kurtosis of residuals and the standardized values in a N (0,1) distribution must be guaranteed. Finally, all the above-ground biomass components data were randomly divided into two groups with a Chow Test, which was a validation method to determine whether the above-ground biomass components data had the same regression model compared with the random sub-samples.

3. Results

3.1. Correlation Analysis of Biomass Variables

The observations of this study indicate that the biomass variables of Caragana korshinskii were significantly higher than those of Sophora viciifolia (p < 0.05) except for height (Table 2), which presented no significant difference between Caragana korshinskii and Sophora viciifolia (p > 0.05). In addition, the correlation analysis indicates that there was a significant correlation between biomass variables of Caragana korshinskii and Sophora viciifolia shrubs (Table 3).
Positive correlations were observed between H, D, leaf biomass, and branch biomass of Caragana korshinskii and Sophora viciifolia shrubs (p < 0.05). There were significantly positive correlations between canopy volume (Cv), leaf biomass, and branch biomass (p < 0.05), while canopy area (Ca) had no significant positive correlations with leaf biomass or branch biomass (p > 0.05). Furthemore, H and D had significantly positive correlations with leaf biomass and branch biomass of the Caragana korshinskii and Sophora viciifolia shrubs (p < 0.05). Thus, H and D were selected to estimate the leaf biomass and branch biomass for the biomass equations.
The response curves of leaf biomass to branch biomass were linear relationship (R2 = 0.7753, p < 0.001; R2 = 0.8331, p < 0.001) (Figure 2), which might give an explanation for the correlations between leaf biomass and branch biomass of the Caragana korshinskii and Sophora viciifolia shrubs.
The confidence interval (CI) and average for above-ground biomass components was 77.53% and 23.8 g for Caragana korshinskii, and 83.31% and 13.2 g for Sophora viciifolia. In other words, leaf and branch biomass components accounted for 77.53% and 83.31%, respectively. However, the line of slope of Caragana korshinskii was above 1.0 and the line of slope of Sophora viciifolia was below 1.0. Therefore, the major part of the biomass of the Caragana korshinskii shrubs was concentrated in the branches and for the Sophora viciifolia shrubs it was concentrated in the leaves.
In addition, the relationship between D, D2H and above-ground biomass was calculated by unit regression (Figure 3 and Figure 4). The results show the positive exponential model between D and above-ground biomass, and the positive linear relationship between D2H and above-ground biomass. Furthermore, the SEE of the above-ground biomass is randomly distributed on the zero level with the increasing D and D2H, which means that the heteroscedasticity of the data was weaker, and it was suitable for estimation model analysis, which indicated that D or D2H with the biomass equations could accurately estimate the above-ground biomass.

3.2. Allometric Model for Biomass

On the basis of the independent variables, to build as many models as possible, the same methodology and test were used to find the optimum model (Table 4). The correlations between above-ground biomass and various independent variables were generally greater. However, the ability to predict the biomass of Caragana korshinskii was better than that of the Sophora viciifolia shrubs (with the higher correlation coefficient), and so this statistical test was applied to these equations.
In all the equations, the Caragana korshinskii and Sophora viciifolia shrubs have the same regression model across all samples. R2 values of the exponential model were higher than the other models with the lower SEE, and the predictability (ANOVA test) was significant (p < 0.05). In the multiple-variables model, H and D were selected, as expected, the optimum model was also model (2), explaining about 90% of the above-ground biomass of the Caragana korshinskii and Sophora viciifolia shrubs (p < 0.001). In addition, the exponential model (2) had the highest R2 values, although R2 values were generally high in models (1) and (3). In most cases, R2 values of D2H were higher than D when models (2) were used to investigate the Caragana korshinskii and Sophora viciifolia shrubs. Thus, model (2) was the best equation to estimate the above-ground biomass using D2H with the higher R2, Adjust R2, F, and the lower SEE. Furthermore, the allometric model of the Caragana korshinskii shrubs was W = 0.076x0.703, the allometric model of the Sophora viciifolia shrubs was W = 0.745x0.080.

3.3. Biomass Model Verification

In this part, these independent variables were firstly checked to see whether they could improve the fitness of the regression model. Secondly, the “intercept” was tested to see whether it had a statistically significantly difference from zero using the t-test. If this was not the case, we removed the “intercept” from our models using the statistical points. Finally, some of the “intercepts” were recommended to remain even though the value was not significantly different from zero, because they help to make a precise prediction. Thus, the mean error (ME), mean absolute error (MAE), total relative error (TRE), mean sybranchatic error (MSE), and mean absolute percentage error (MPSE) were calculated by three models, and all these parameters were low, which indicated that the allometric model had an ecological significance according to one-way ANOVA (Table 5).
In order to test the biomass estimation model in the field, the observed and estimated biomass values of the Caragana korshinskii (n = 125) and Sophora viciifolia (n = 129) shrubs were compared (Figure 5). The results show that the response curves of the observed and estimated values of the Caragana korshinskii and Sophora viciifolia shrubs had a linear relationship (R2 = 0.8499, p < 0.001, SEE = 26.3856; R2 = 0.9007, p < 0.001, SEE = 24.1357), which indicated that the estimated values were quite close to the observed values, suggesting that this allometric model could be used to estimate the above-ground biomass of the Caragana korshinskii and Sophora viciifolia shrubs in the Loess Plateau.

4. Discussion

Now, constructing an empirical model to estimate forest biomass from easily measurable indicators (height, DBH, etc.) is the most effective and environmentally friendly way to evaluate forest biomass and carbon storage [33,34,35]. This study also used this method and made a contribution to above-ground biomass prediction by developing the equation models. Variables for the above-ground biomass components (leaf and branch) with three procedures were estimated: (1) estimating leaf and branch biomass by using the same method; (2) adding the optimum regression models of each leaf and branch biomass; and (3) the joint generalized least squares was used to iteratively force the sum of the coefficients to satisfy the equation models [21,22,23,24,36]. Finally, the estimated and observed values were used to test the biomass estimation model.
The results of this study show that that the biomass variables of Caragana korshinskii were significantly higher than those of Sophora viciifolia (p < 0.05) except for height (Table 2), which may be because Caragana korshinskii and Sophora viciifolia have different growth characteristics, namely, they have different transpiration rate and water use efficiency [37,38]. A previous study found that it was not ideal to use a single variable to simulate the biomass regression equation. Only the proper combination of each variable can estimate the biomass of shrubs accurately [39]. In the present study, H and D were selected to estimate the leaf biomass and branch biomass by the best fitting biomass equations for Caragana korshinskii and Sophora viciifolia, which confirmed that a single variable cannot reflect the above-ground biomass.
In this study, the biomass equation models presented a cross-validated R2, which indicated the suitability to calculate the above-ground biomass of shrubs. Therefore, we can conclude that the allometric model of the Caragana korshinskii shrubs was W = 0.076 x0.703, and the allometric model of the Sophora viciifolia shrubs was W = 0.745x0.080 (Table 4). These results are in agreement with the other studies [4,10,15], which also certify the allometric model for above-ground biomass estimation.
Furthermore, other studies have shown that the optimal biomass estimation models of different shrub species were different, while the biomass estimation models for the same species were mostly the same; this is due to the differences between species [40]. However, in this study, the biomass models of Caragana korshinskii and Sophora viciifolia were both based on the form y = axb, which may be due to Caragana korshinskii and Sophora viciifolia both belonging to the leguminous family.
The fitting of the allometric model is based on logarithmic transformation of variables, which may produce a higher homogeneity of variance [4,10]. Furthermore, the present study found that the harvested biomass of the percentage of outliers in all samples (measured) was highly significant (p < 0.001) (Table 5). The percentage of outliers was more than 10% in a few cases, so in future studies, increasing the sample scale might act to reduce the number of outliers, though there was a normal level of outliers. In addition, as regards the progress of these two shrubs’ growth, the exponential biomass–size relationships may be transformed into linear models [12,14,32]. In view of these considerations, other important variables related to the growth of shrubs should be considered in above-ground biomass estimation [12,14].

5. Conclusions

This study explored the above-ground biomass models of two common shrubs (Caragana korshinskii and Sophora viciifolia) in the Loess Plateau, China. For both shrubs, models based on y = axb were suitable to estimate above-ground biomass using the height and crown diameter. In addition, the response curves of the observed and estimated values of the Caragana korshinskii and Sophora viciifolia shrubs also confirmed that our above-ground biomass estimation models were suitable. Because of the impact of climate, soil, terrain, and other factors on the growth of shrubs, especially in arid and semi-arid areas, the growth of shrubs can be greatly affected by the soil moisture status because of the great variation in precipitation resources in interannual or annual distribution; therefore, the same species may have different estimation models under different site conditions. Hence, the above-ground biomass models for the Caragana korshinskii and Sophora viciifolia shrubs may be unsuitable in other regions. In the future, studies in this field should attempt to use environmental factors as the explanatory variables.

Author Contributions

Project Administration, Supervision, Writing-Review & Editing, S.A.; Methodology, Software, Investigation, Formal Analysis, Resources, Data Curation, Writing-Original Draft Preparation, Y.D.; Investigation, Software, Data Curation, Y.Y.

Funding

This study was supported by the National Natural Sciences Foundation of China (41671280) and Special-Funds of Scientific Research Programs of State Key Laboratory of Soil Erosion and Dryland Farming in the Loess Plateau (A314021403-C6).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Johnson, D.A.; Tieszen, L.L. Aboveground biomass allocation, leaf growth, and photosynthesis patterns in tundra plant forms in arctic Alaska. Oecologia 1976, 24, 159–173. [Google Scholar] [CrossRef] [PubMed]
  2. Chapin, F.S.; Matson, P.A.; Vitousek, P. Principles of Terrestrial Ecosystem Ecology; Springer Science Business Media: Berlin, Germany, 2011; pp. 157–181. [Google Scholar]
  3. Conti, G.; Enrico, L.; Casanoves, F.; Díaz, S. Shrub biomass estimation in the semiarid Chaco forest: A contribution to the quantification of an underrated carbon stock. Ann. For. Sci. 2013, 70, 515–524. [Google Scholar] [CrossRef]
  4. Chave, J.; Réjou-Méchain, M.; Búrquez, A.; Chidumyo, E.; Colgan, M.S.; Delitti, W.B.C. Improved allometric models to estimate the aboveground biomass of tropical trees. Glob. Chang. Biol. 2014, 20, 3177–3190. [Google Scholar] [CrossRef]
  5. Thomas, S.M.; Abbott, K.C.; Moloney, K.A. Effects of aboveground herbivory on plants with long-term belowground biomass storage. Theor. Ecol. 2017, 10, 35–50. [Google Scholar] [CrossRef]
  6. Slik, J.W.; Paoli, G.; McGuire, K.; Amaral, L.; Barroso, J.; Bastian, M. Large trees drive forest aboveground biomass variation in moist lowland forests across the tropics. Glob. Ecol. Biogeogr. 2013, 22, 1261–1271. [Google Scholar] [CrossRef]
  7. Zheng, D.; Rademacher, J.; Chen, J.; Crow, T.; Bresee, M.; Moine, J.L.; Ryu, S.R. Estimating aboveground biomass using Landsat 7 ETM+ data across a managed landscape in northern Wisconsin, USA. Remote Sens. Environ. 2004, 93, 402–411. [Google Scholar] [CrossRef]
  8. Arturo, A.S.C.; Javier, T.G.E.; Alberto, A.C.O.; Javier, J.P.; Aurelio, G.T.M.; René, V.L. Construction of aboveground biomass models with remote sensing technology in the intertropical zone in Mexico. J. Geogr. Sci. 2012, 22, 669–680. [Google Scholar] [CrossRef]
  9. Segura, M.; Kanninen, M.; Suárez, D. Allometric models for estimating aboveground biomass of shade trees and coffee bushes grown together. Agrofor. Syst. 2006, 68, 143–150. [Google Scholar] [CrossRef]
  10. Zolkos, S.G.; Goetz, S.J.; Dubayah, R. A meta-analysis of terrestrial aboveground biomass estimation using lidar remote sensing. Remote Sens. Environ. 2013, 128, 289–298. [Google Scholar] [CrossRef]
  11. Baskerville, G.L. Use of logarithmic regression in the estimation of plant biomass. Can. J. Forest Res. 1972, 2, 49–53. [Google Scholar] [CrossRef]
  12. Parton, W.J.; Scurlock, J.M.O.; Ojima, D.S.; Gillmanov, T.G.; Scholes, R.J.; Schimel, D.S.; Kirchner, T.; Meanaut, J.C.; Seastedt, T.; Garcia Moya, E. Observations and modeling of biomass and soil organic matter dynamics for the grassland biome worldwide. Glob. Biogeochem. Cycles 1993, 7, 785–809. [Google Scholar] [CrossRef]
  13. Jia, G.M.; Liu, B.R.; Wang, G.; Zhang, B. The microbial biomass and activity in soil with shrub (Caragana korshinskii K.) plantation in the semi-arid Loess Plateau in China. Eur. J. Soil Biol. 2010, 46, 6–10. [Google Scholar] [CrossRef]
  14. Englhart, S.; Keuck, V.; Siegert, F. Aboveground biomass retrieval in tropical forests-The potential of combined X-and L-band SAR data use. Remote Sens. Environ. 2011, 115, 1260–1271. [Google Scholar] [CrossRef]
  15. Cheng, J.; Wu, G.L.; Zhao, L.P.; Li, Y.; Li, M.; Cheng, J.M. Cumulative effects of 20-year exclusion of livestock grazing on above-and belowground biomass of typical steppe communities in arid areas of the Loess Plateau, China. Plant Soil Environ. 2011, 57, 96–102. [Google Scholar] [CrossRef]
  16. Tang, S.; Zhang, H.; Xu, H. Study on establish and estimate method of compatible biomass model. Sci. Silvae Sin. 2000, 36, 19–27. [Google Scholar]
  17. Yoshida, T.; Oshima, Y.; Matsumura, Y. Gasification of biomass model compounds and real biomass in supercritical water. Biomass Bioenergy 2004, 26, 71–78. [Google Scholar]
  18. Roos, A.M.; Schellekens, T.; Kooten, T.V.; Wolfshaar, K.V.D.; Claessen, D.; Persson, L. Simplifying a physiologically structured population model to a stage-structured biomass model. Theor. Popul. Biol. 2008, 73, 47–62. [Google Scholar] [CrossRef]
  19. Wang, Q.; Bai, Y.; Gao, H.; He, J.; Chen, H.; Chesney, R.C.; Kuhn, N.J.; Li, H.W. Soil chemical properties and microbial biomass after 16 years of no-tillage farming on the Loess Plateau, China. Geoderma 2008, 144, 502–508. [Google Scholar] [CrossRef]
  20. Wang, X.C.; Li, J. Evaluation of crop yield and soil water estimates using the EPIC model for the Loess Plateau of China. Math. Comput. Model. 2010, 51, 1390–1397. [Google Scholar]
  21. Lu, D. Aboveground biomass estimation using Landsat TM data in the Brazilian Amazon. Int. J. Remote Sens. 2005, 26, 2509–2525. [Google Scholar] [CrossRef]
  22. Litton, C.M.; Boone, K.J. Allometric models for predicting aboveground biomass in two widespread woody plants in Hawaii. Biotropica 2008, 40, 313–320. [Google Scholar] [CrossRef]
  23. Anaya, J.A.; Chuvieco, E.; Palacios-Orueta, A. Aboveground biomass assessment in Colombia: A remote sensing approach. Forest Ecol. Manag. 2009, 257, 1237–1246. [Google Scholar] [CrossRef]
  24. Lambert, M.C.; Ung, C.H.; Raulier, F. Canadian national tree aboveground biomass equations. Can. J. For. Res. 2005, 35, 1996–2018. [Google Scholar] [CrossRef]
  25. Jiang, W.; Cheng, Y.; Yang, X.; Yang, S.L. Chinese Loess Plateau vegetation since the Last Glacial Maximum and its implications for vegetation restoration. J. Appl. Ecol. 2013, 50, 440–448. [Google Scholar] [CrossRef]
  26. Lang, Y.; Zhiyuan, R.; Chong, Z.; Liu, W. Vegetation cover, climate and human activities on the Loess Plateau. Resour. Sci. 2014, 36, 166–174. (In Chinese) [Google Scholar]
  27. Ma, K.P.; Huang, J.H.; Yu, S.L.; Chen, L.Z. Plant community diversity in dongliang mountain, Beijing, China: II. Sepcies richness, evennes and species diersityes. Acta Ecol. Sin. 1995, 15, 268–277. (In Chinese) [Google Scholar]
  28. Peek, J.M. Relation of canopy area and volume to production of three woody species. Ecology 1970, 51, 1098–1101. [Google Scholar] [CrossRef]
  29. Yamaoka, K.; Nakagawa, T.; Uno, T. Application of Akaike’s information criterion (AIC) in the evaluation of linear pharmacokinetic equations. J. Pharmacokinet. Biopharm. 1978, 6, 165–175. [Google Scholar] [CrossRef]
  30. Posada, D.; Buckley, T.R. Model selection and model averaging in phylogenetics: Advantages of Akaike information criterion and Bayesian approaches over likelihood ratio tests. Syst. Biol. 2004, 53, 793–808. [Google Scholar] [CrossRef]
  31. Muukkonen, P.; Makipaa, R. Empirical biomass models of understorey vegetation in boreal forests according to stand and site attributes. Boreal Environ. Res. 2006, 11, 355–369. [Google Scholar]
  32. Gertrudix, R.R.P.; Montero, G.; Delrio, M. Biomass models to estimate carbon stocks for hardwood tree species. For. Syst. 2012, 21, 42–52. [Google Scholar] [CrossRef]
  33. Luo, Y.J.; Zhang, X.Q.; Wang, X.K.; Zhu, J.H.; Hou, Z.H.; Zhang, Z.J. Forest biomass estimation methods and their prospects. Sci. Silvae Sin. 2009, 8, 129–134. [Google Scholar]
  34. Zeng, W.S. A review of studies of shrub bio-mass modeling. World For. Res. 2015, 28, 31–36. [Google Scholar]
  35. Komiyama, A.; Ong, J.E.; Poungparn, S. Allometry, bio-mass and productivity of mangrove forests: A review. Aquat. Bot. 2008, 89, 128–137. [Google Scholar] [CrossRef]
  36. Ter-Mikaelian, M.T.; Korzukhin, M.D. Biomass equations for sixty-five North American tree species. For. Ecol. Manag. 1997, 97, 1–24. [Google Scholar] [CrossRef] [Green Version]
  37. Bu, C.F.; Liu, G.B.; Zhang, W.H. Growth of characteristics of Sophora viciifolia population in the hilly and gully region of Loess Plateau. Acta Bot. Boreali-Occident. Sin. 2004, 10, 1792–1797. [Google Scholar]
  38. Bu, C.F.; Liu, G.B. Study on the transpiration of Sophara viciifolia in the hilly and gully Loess Plateau. Bull. Bot. Res. 2005, 1, 64–68. [Google Scholar]
  39. Zeng, Z.Y.; Liu, Q.J.; Zeng, H.Q. Study of biomass models of some kinds of shrubs in Qianyanzhou, Jiangxi Province. J. Fujian For. Sci. Technol. 2005, 4, 68–72. [Google Scholar]
  40. Luo, Y.K.; Fang, J.Y.; Hu, H.F. Biomass estimation models and allocation patterns of 14 shrub species in Mountain Luya, Shanxi, China. Chin. J. Plant Ecol. 2017, 1, 115–125. [Google Scholar]
Figure 1. Location of the study area and layout of the shrub plot for the measurement of vegetation characteristics in the Loess Plateau.
Figure 1. Location of the study area and layout of the shrub plot for the measurement of vegetation characteristics in the Loess Plateau.
Sustainability 11 01674 g001
Figure 2. Relationships between the branch biomass and leaf biomass of Caragana korshinskii (A, n = 861) and Sophora viciifolia (B, n = 828). In the Figure: BB—branch biomass, LB—leaf biomass.
Figure 2. Relationships between the branch biomass and leaf biomass of Caragana korshinskii (A, n = 861) and Sophora viciifolia (B, n = 828). In the Figure: BB—branch biomass, LB—leaf biomass.
Sustainability 11 01674 g002
Figure 3. Relationships between the above-ground biomass and its morphological variables for Caragana korshinskii. (A,B) AB—above-ground biomass, (C,D) SEE—standard error of the estimation; sample content = 861.
Figure 3. Relationships between the above-ground biomass and its morphological variables for Caragana korshinskii. (A,B) AB—above-ground biomass, (C,D) SEE—standard error of the estimation; sample content = 861.
Sustainability 11 01674 g003
Figure 4. Relationships between the above-ground biomass and its morphological variables of Sophora viciifolia. (A,B) AB—above-ground biomass, (C,D) SEE—standard error of the estimation; sample content = 828.
Figure 4. Relationships between the above-ground biomass and its morphological variables of Sophora viciifolia. (A,B) AB—above-ground biomass, (C,D) SEE—standard error of the estimation; sample content = 828.
Sustainability 11 01674 g004
Figure 5. Relationships between the observed and estimated values of above-ground biomass of Caragana korshinskii and Sophora viciifolia. Capacity of model (2) for predicting shrubs above-ground biomass, considering all species together. Each dot corresponds to an individual weighted shrub in W = 0.076 × 0.703, W = 0.745x0.080. (A) Caragana korshinskii, sample content = 125, (B) Sophora viciifolia, sample content = 129, Number of plants sampled (i.e., measured), harvested (for biomass estimation), R2, F-value of ANOVA test, and average and standard deviation of residuals are shown. All models were highly significant (p < 0.001).
Figure 5. Relationships between the observed and estimated values of above-ground biomass of Caragana korshinskii and Sophora viciifolia. Capacity of model (2) for predicting shrubs above-ground biomass, considering all species together. Each dot corresponds to an individual weighted shrub in W = 0.076 × 0.703, W = 0.745x0.080. (A) Caragana korshinskii, sample content = 125, (B) Sophora viciifolia, sample content = 129, Number of plants sampled (i.e., measured), harvested (for biomass estimation), R2, F-value of ANOVA test, and average and standard deviation of residuals are shown. All models were highly significant (p < 0.001).
Sustainability 11 01674 g005
Table 1. The basic characteristics of sample plots.
Table 1. The basic characteristics of sample plots.
ShrubsSiteSample ContentLongitude and LatitudeSlope/(%)Altitude/(m)Soil Temperature/(℃)
0–10 cm
10–20 cm
Accompanying Species
Caragana korshinskii19336°47′30″ NSE115121.8 ± 1.7 a20.1 ± 0.6 bLespedeza daurica, Artemisia giraldii, Patriniascaniosaefolia
109°16′10″ E15~23
216136°49′01″ NSW132922.6 ± 0.8 a20.7 ± 1.3 bLespedeza daurica, Artemisia giraldii, Patriniascaniosaefolia
109°14′46″ E10~20
316136°49′04″ NNE134727.4 ± 1.9 a25.6 ± 0.8 bLespedeza daurica, Artemisia giraldii, Melilotussuavcolen
109°14′43″ E23~35
416036°43′54″ NNE121920.8 ± 1.2 a20.7 ± 1.1 aArtemisia vestita, Lespedeza daurica, Artemisia giraldii
109°15′30″ E9~18
515836°45′36″ NSW130121.2 ± 0.5 a19.9 ± 0.3 bArtemisia vestita, Lespedeza daurica, Artemisiascoparia
109°16′01″ E10~20
612836°44′49″ NSW123622.0 ± 0.8 a20.5 ± 0.6 bArtemisia vestita, Artemisiascoparia, Artemisia frigida
109°16′13″ E9~13
Sophora
viciifolia
111636°47′20″ NS125624.3 ± 5.1 a22.1 ± 2.3 bLespedeza daurica, Bothriochloaischaemum, Stipabungeana
109°16′30″ E19~35
218036°47′35″ NW124823.4 ± 2.1 a21.2 ± 1.9 bLespedeza daurica, Bothriochloaischaemum, Stipabungeana
109°16′41″ E14~20
313936°47′58″ NSE127125.5 ± 3.2 a23.3 ± 1.7 bLespedeza daurica, Bothriochloaischaemum, Lespedeza daurica
109°16′47″ E8~12
416536°44′42″ NNE128325.2 ± 1.8 a24.4 ± 1.3 bLespedeza daurica, Bothriochloaischaemum, Thyme
109°14′35″ E11~19
59636°44′56″ NNW126927.9 ± 0.6 a25.3 ± 0.4 bLespedeza daurica, Artemisia giraldii, Bothriochloaischaemum
109°14′01″ E8~23
613236°45′03″ NNE115626.9 ± 0.9 a26.1 ± 1.1 aArtemisia giraldii, Patriniascaniosaefolia, Lespedeza daurica
109°15′11″ E18~30
Note: Values followed by lowercase letters within columns are significantly different for 0.05 using Tukey’s method. The same applies below.
Table 2. Ecological parameters for the Caragana korshinskii and Sophora viciifolia.
Table 2. Ecological parameters for the Caragana korshinskii and Sophora viciifolia.
ShrubsSiteHeight/(cm)Diameter/(mm)Canopy Area/(cm2 104)Canopy Volume/(cm3 106)Coverage/(%)Leaf Biomass/(g)Branch Biomass/(g)
Caragana korshinskii1156 ± 325.2 ± 0.62.91 ± 0.214.54 ± 0.5658 ± 6358 ± 35589 ± 39
2137 ± 214.8 ± 0.52.32 ± 0.583.18 ± 0.8967 ± 9321 ± 29512 ± 35
3125 ± 354.6 ± 0.91.66 ± 0.352.08 ± 0.5283 ± 8305 ± 16456 ± 29
4113 ± 193.2 ± 0.31.34 ± 0.261.51 ± 0.4192 ± 12294 ± 32478 ± 34
5109 ± 263.4 ± 0.41.22 ± 0.101.33 ± 0.3285 ± 5287 ± 18423 ± 25
6151 ± 145.1 ± 0.82.77 ± 0.434.18 ± 0.7862 ± 8389 ± 42563 ± 37
CV/(%)14.7919.8135.9148.9618.7312.2912.76
Mean132 ± 19 a 4.4 ± 0.9 a2.04 ± 0.73 a2.80 ± 1.37 a75 ± 14 a326 ± 40 a503 ± 64 a
F89.69156.98178.32135.69103.56158.96147.89
Sophora viciifolia1123 ± 114.2 ± 0.21.15 ± 0.091.41 ± 0.4523 ± 3324 ± 36456 ± 52
296 ± 232.9 ± 0.60.83 ± 0.120.79 ± 0.2135 ± 5257 ± 23356 ± 43
3108 ± 343.1 ± 0.50.90 ± 0.250.97 ± 0.1938 ± 9269 ± 18324 ± 51
4145 ± 254.8 ± 0.81.92 ± 0.082.78 ± 0.2320 ± 6395 ± 25423 ± 41
5112 ± 183.5 ± 0.41.09 ± 0.151.22 ± 0.3543 ± 8298 ± 31389 ± 38
697 ± 163.1 ± 0.30.86 ± 0.140.83 ± 0.1439 ± 8278 ± 17362 ± 31
CV/(%)16.8120.7936.4756.0528.2316.6912.51
Mean114 ± 18 a3.6 ± 0.7 b1.13 ± 0.41 b1.33 ± 0.75 b33 ± 9 b304 ± 51 a385 ± 48 b
F96.38156.41147.85125.87214.79134.57126.84
Note: Values followed by lowercase letters within columns are significantly different for 0.05 using Tukey’s method. The same applies below.
Table 3. Correlation analysis of each variable of Caragana korshinskii and Sophora viciifolia.
Table 3. Correlation analysis of each variable of Caragana korshinskii and Sophora viciifolia.
ShrubsItemHeightDiameterCanopy AreaCanopy VolumeLeaf BiomassBranch Biomass
Caragana korshinskiiHeight1.000
Diameter0.724 *1.000
Canopy area0.523 *0.3211.000
Canopy volume0.4690.4560.759 **1.000
Leaf biomass0.826 **0.856 **0.2230.566 *1.000
Branch biomass0.879 **0.891 **03960.547 *0.899 **1.000
Sophora viciifoliaHeight1.000
Diameter0.789 **1.000
Canopy area0.587 *0.563 *1.000
Canopy volume0.4360.525 *0.689**1.000
Leaf biomass0.875 **0.845 **0.3680.698 **1.000
Branch biomass0.786 **0.756 **0.4230.513 *0.903 **1.000
*, p < 0.05; **, p < 0.01.
Table 4. Biomass regression models of Caragana korshinskii and Sophora viciifolia.
Table 4. Biomass regression models of Caragana korshinskii and Sophora viciifolia.
ShrubsOutliers/%Equation Parameters
VariateFormulaEquationsa db dcR2Adjust R2SEEF
Caragana korshinskiia13.6xc = D1W = a + bx + cx213.424−6.2031.3080.729 ***0.70135.154560.709
9.2x = D2HW = a + bx + cx24.2820.0050.0070.770 ***0.62533.368968.813
6.8x = D2W = ax b0.8641.739 0.765 ***0.63231.6821117.664
10.7x = D2HW = ax b0.0760.703 0.898 ***0.7718.3691748.984
8.5x = D3W = a + blnx−18.35422.196 0.695 ***0.52435.389483.002
9.3x = D2HW = a + blnx−49.7529.039 0.775 ***0.60532.347777.264
Sophora viciifoliab6.2x = D1W = a + bx + cx20.0010.1252.7270.683 ***0.59831.157130.596
3.8x = D2HW = a + bx + cx20.1500.1531.3710.829 ***0.7727.4921453.608
13.6x = D2W = ax b0.1994.607 0.727 ***0.65829.568152.494
5.3x = D2HW = ax b0.7450.080 0.868 ***0.7536.9872614.672
4.9x = D3W = a + blnx2.2264.653 0.681 ***0.60234.714145.606
8.2x = D2HW = a + blnx6.025−27.358 0.805 ***0.64711.4931575.173
Note: a Caragana korshinskii sample content = 861; b Sophora viciifolia sample content = 828. Models for all dependent variables are of these three forms, where y is the dependent variable (above-ground biomass), c x is the predictor variable (D or D2H), and d a and b are constants in the equation. Percentage of outliers, R2, F-value of ANOVA test, and average and standard deviation of residuals are shown. *, p < 0.05; **, p < 0.01; ***, p < 0.001.
Table 5. Testing of accuracy of biomass regression models of Caragana korshinskii and Sophora viciifolia.
Table 5. Testing of accuracy of biomass regression models of Caragana korshinskii and Sophora viciifolia.
ShrubsEquationsGoodness Indicator of Fitting
MEa/gMAE b/gTREc/%MSEd/%MPSEe/%
Caragana korshinskiiW = a + bx + cx213.25643.2572.36915.03632.568
W = ax b0.80315.369−0.063−1.05815.036
W = a + blnx9.12532.0258.2638.21423.256
Sophora viciifoliaW = a + bx + cx210.03934.177−0.1593.69212.035
W = axb0.23619.256−1.2562.3699.285
W = a + blnx5.12323.0545.0328.03829.877
Note: Diameter at breast height versus total tree height models for Caragana korshinskii and Sophora viciifolia used to develop the allometric models. a Mean while minimum values of mean error (ME), b mean absolute error (MAE), c total relative error (TRE), d mean sybranchatic error (MSE), and e mean absolute percentage error (MPSE).

Share and Cite

MDPI and ACS Style

Dou, Y.; Yang, Y.; An, S. Above-Ground Biomass Models of Caragana korshinskii and Sophora viciifolia in the Loess Plateau, China. Sustainability 2019, 11, 1674. https://0-doi-org.brum.beds.ac.uk/10.3390/su11061674

AMA Style

Dou Y, Yang Y, An S. Above-Ground Biomass Models of Caragana korshinskii and Sophora viciifolia in the Loess Plateau, China. Sustainability. 2019; 11(6):1674. https://0-doi-org.brum.beds.ac.uk/10.3390/su11061674

Chicago/Turabian Style

Dou, Yanxing, Yang Yang, and Shaoshan An. 2019. "Above-Ground Biomass Models of Caragana korshinskii and Sophora viciifolia in the Loess Plateau, China" Sustainability 11, no. 6: 1674. https://0-doi-org.brum.beds.ac.uk/10.3390/su11061674

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