Next Article in Journal
Contrasting Response to Drought and Climate of Planted and Natural Pinus pinaster Aiton Forests in Southern Spain
Previous Article in Journal
The Characteristics of Soil C, N, and P Stoichiometric Ratios as Affected by Geological Background in a Karst Graben Area, Southwest China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Removing the Scaling Error Caused by Allometric Modelling in Forest Biomass Estimation at Large Scales

1
Faculty of Health Sciences, University of Ottawa, Ottawa, ON K1N 6N5, Canada
2
Research Center for Ecological Forecasting and Global Change, Northwest A&F University, Yangling 712100, China
3
Ecological Modeling and Carbon Science, Department of Biology Science, University of Quebec at Montreal, Montreal, QC H3C 3P8, Canada
*
Author to whom correspondence should be addressed.
Forests 2019, 10(7), 602; https://doi.org/10.3390/f10070602
Submission received: 14 June 2019 / Accepted: 20 July 2019 / Published: 22 July 2019
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)

Abstract

:
To estimate the responses of forest ecosystems, most relationships in biological systems are described by allometric relationships, the parameters of which are determined based on field measurements. The use of existing observed data errors may occur during the scaling of fine-scale relationships to describe ecosystem properties at a larger ecosystem scale. Here, we analyzed the scaling error in the estimation of forest ecosystem biomass based on the measurement of plots (biomass or volume per hectare) using an improved allometric equation with a scaling error compensator. The efficiency of the compensator on reducing the scaling error was tested by simulating the forest stand populations using pseudo-observation. Our experiments indicate that, on average, approximately 94.8% of the scaling error can be reduced, and for a case study, an overestimation of 3.6% can be removed in practice from a large-scale estimation for the biomass of Pinus yunnanensis Franch.

1. Introduction

Most relationships in biological systems are nonlinear between two organs, subsystems, and groups [1]. These relationships are usually described as allometric [2,3,4] and their parameters can be determined by regression analysis based on observations [5,6]. That is, our ecological knowledge mostly comes from small, easily measured sample plots [7,8]. These plot observations and estimates can be extrapolated to large-scale estimates [9] and global changes [10,11]. The technology of upscaling using allometric relationships was developed gradually. Originally, the allometry theory was used to study individual living bodies, e.g., single tree, and the relative growth of their different organs [2]. This theory has been widely applied in various fields, from biology to sociology [12,13]. In the research field of forest ecosystems, allometric equations were further extended for depicting the relationship between different tree parts in order to address various issues [5,14,15]. Since the 1980s researchers have compiled several large databases across the world [16,17,18] according to field measurements of forest biomass which had been converted from tree level to plot level by measurers. These hectare-based data provide the opportunity to test allometric relationships at both the plot-scale and larger scales. Using this methodology of model testing, the analysis in our study was carried out based on plot measurements (biomass or volume per hectare) rather than individual trees—that is, to extrapolate biomass from plot-scale to large-scale.
When describing the properties of a large-scale ecosystem using existing fine-scale knowledge, the scaling error occurs during the transforming processing [19,20,21] (also see Appendix A). This causes prediction bias on large-scales [20] because of the method of weighted mean [19,22]. Specifically, during the scaling of ecological properties from field plots up to a region, the variability among the plots can be aggregated [23] by using a nonlinear equation that only reflects the behavior of the plots [23,24,25]. The scaling errors usually include two types that respectively correspond to two stages in forest biomass estimation. The first stage is volume prediction, which is carried out based on large sampling and hierarchical data structure in national forest inventories (NFI) [26], in which the sampling error would exist [27]. Following it, the second stage converts forest biomass from volume information (m3 ha−1) released in NFI reports. As a simple conversion, the biomass expansion factor (BEF) and biomass conversion expansion factor (BCEF) are generally used to estimate forest biomass [28]. For more accurate estimations, the volume-based biomass equations were also frequently employed in past decades [29,30,31,32,33,34]. This requires removing the scaling error at the regional level. Our study focuses on the second stage and second scaling error. This stage is implemented especially for those regions or countries that have NFI without biomass data. In this stage, the scaling error is not involved with trees because upscaling is only based on the means per unit area.
According to the computational principle, the form of the allometric equation affects the magnitude of scaling error [23]. Zhou et al. [25] have reported a quantitative analysis to reduce scaling error caused by power-law equations. However, besides the general power-law equation, an improved allometric equation [35] has not been addressed. The improved equation incorporates a variable allometric ratio (VAR) instead of a constant ratio. It is more complex in form than conventional power-law equations. Because it can explain statistically significant variation and give well-fitting curves according to field measurements [35,36] for some special species [37], this equation increases the accuracy of biomass estimates, particularly for the species (Pinus yunnanensis Franch.) addressed in our study. Our analysis adopted the improved allometric equation and presented a quantitative approach on the issues of scaling error and error compensation. The presented method provides guidance for statistical and mechanistic analyses of errors for scaling forest biomass from fine-scale to large-scales.
Our objectives are to test the margin of scaling error, to derive a user-friendly compensator, and to examine its ability in correcting scaling error in the application of the improved allometric equation for large-scale estimation of forest ecosystem biomass. We only focus on correcting the scaling error itself and leave out other types of errors in the biomass estimation.

2. Methods

2.1. Derivation

Ruark et al. [35] presented an improved allometric VAR equation y = rxkeux where y denotes total biomass density (t ha−1), x is stand stocking volume (m3 ha−1), and r, k, and u are parameters greater than zero. The error compensator for this equation is derived by expanding the scaling error expression (error = Y1 − Y2, see Appendix A) as Taylor series, and utilizing the central moment method to solve the items in Taylor series. Eventually the difference of scaling error before and after compensation is proved by an analytical expression. This expression includes two parameters that are used in the allometric equation. The variables in the expression reflect the feature of the estimated forest ecosystem, i.e., nonlinearity of volume-to-biomass relationship.

2.2. Simulation

To test the algorithm and efficiency of the error compensator, we construct a pseudo population, and compare our biomass estimate with the true value of the population. The error compensation is simulated by scaling tree biomass from stand level up to a large region consisting of 100 stands with different stages of development (see Supporting Information). The area and volume of the 100 stands are produced by pseudo-random number generator under uniform probability distribution. The area ranges from 1.0 to 10.0 ha and stand stocking volume ranges from 15 to 300 (m3 ha−1). Thus, the true value of the total amount of regional forest biomass is recognized by accumulating all stand biomasses. Based on the pseudo population, numerous trials provide a uniform distribution that gives how much error can be reduced and information on how many times the convergence of compensated percentage will be reached.

2.3. Data for Case Study

A case study of forest biomass estimation is carried out to test the difference that could come from scaling error at the provincial scale. We estimated the biomass of Pinus yunnanensis in Yunnan province of China. This province is located in the southwest border of China, between latitude 21°8′32″–29°15′8″ N and longitude 97°31′39″–106°11′47″ E. The total forested area and stem volume for the species is 2.92 million ha and 222.43 million m3 in the province [38]. The data of area and volume were classified as different age groups: young, middle, near-mature, mature, and over-mature. The range of mean stocking volume was estimated according to a literature review [38,39,40]. The variance of stocking volume is computed as the between-group variance (refer to 3.3. Comparison), yet the within-group variance could not be computed because the detailed data are not released to public in China. The field measurements used in the case study were collected from a published database [18]. This database compiled 34 plot measurements reported by 24 studies that carried out field studies in different years. Due to the costly destructive measurement of forest biomass, each of the studies only addressed one or a few plot works. All measurements have been scaled from tree-level to plot-level and hectare-based values by data measurers. All biomass values are measured by oven-drying different parts of tree samples [18].

3. Results and Discussions

3.1. Error Compensator

All symbols and their descriptions refer to Table 1. The improved allometric equation g(x) = xkeux can be approximated at x = μ as
g ( x ) q = 0 2 [ ( 1 / q ! )   g ( q ) ( μ ) ( x μ ) q ] ,
by using a 2nd-order Taylor series expansion and omitting the remainder term. Taking expectations on both sides, then we see that
μ g E { q = 0 2   [ ( 1 / q ! )   g ( q ) ( μ ) ( x μ ) q ] } ,
After expanding the right side of Equation (2), and then substituting g(x) = xk eux into the expansion, it gives the following expression:
μgμke + 0.5μk – 2 e(k2 − 2kuμk + u2μ2)σ2,
By substituting μg into error = Ar (μgμk e), it leads to a simple form (also refer to Appendix A)
error ≈ 0.5Ak − 2e(k2k − 2kuμ + u2μ2)σ2,
for any probability distribution. Then we have the simple form of the error per unit area as:
Φ(μ, σ2) = 0.5k − 2e (k2k − 2kuμ + u2μ2)σ2,
This is the compensator for correcting the scaling error. Thus, the ratio |Φ(μ, σ2)|/(ke) of regional biomass density becomes:
η(μ, σ2) = 0.5 μ−2(k2k − 2kuμ + u2μ2)σ2,
Rewrite it, we have:
η(μ, σ2) = 0.5μ−2k(k − 1)σ2 + u(0.5u−1)σ2,
where the domains are xminμ ≤ xmax, 0 ≤ σ2 ≤ Dm, and 0.8 ≤ k ≤ 1 (in this study). Notes: (1) Dm is derived by assuming that xi has the largest variance for a given expectation if and only if xi equals either xmin or xmax; (2) μ (xs) and νs are used for calculating μg.
For a distribution between xmin and xmax, the definition of variance implies that a variance will be largest for a given expectation while all samples (xi) are either xmin or xmax. The maximum variance can be expressed as
Dm = −(μ − xmin)2 + (μ − xmin) (xmax − xmin),
where μ represents regional stocking volume. This is for any probability distribution where xmin and xmax are known. Here μ and Dm define a parabola going downwards in the μσ2 plane. Dm becomes a boundary of the curved surface η(μ, σ2) (Figure 1). The derivation of Dm is detailed in Appendix A.

3.2. Efficiency of Reducing the Error

The efficiency and ability of the compensator to reduce scaling error is tested by 5000 trials. The mean of total biomass amount (true value) and the mean of modelled biomass amounts in the pseudo region are 44.7 and 46.3 (t), respectively (Figure 2a). It shows an approximate 3.5% overestimate due to the scaling error. After compensating, around 94.8% of the scaling error is reduced on average. The average level illustrates that compensation rate tends to become stable with the increase of trials. When trials approach 4000 times, Figure 2b indicates a convergence. This implies that most of the scaling error should be removed using the compensator.
The parameters in Equation (7) can explain how the scaling error is removed for the VAR equation (y = rxk eux). Equation (7) contains two terms. The first is the same as the compensation rate for the typical allometric equation y = rxk, and the second is an addition for applying to the VAR equation. Both k and u affect the extent of error compensation. Because the VAR equation has a larger curvature than the typical allometric equation [35], the absolute value of η(μ, σ2) generally increases after adding the second term u(0.5u−1)σ2 to the first term in Equation (7). It means that the scaling error is larger using VAR than typical allometric equation. In extraordinary instances for some case, the volume-to-biomass relationship was considered as linear [33]. This makes u and k become zero and 1.0 respectively, and no scaling error will be caused in large-scale biomass estimation.

3.3. Comparison

The results of the case study exhibit the percentage of scaling error that may occur from a biomass estimation of provincial scale forest ecosystem. The volume-to-biomass equation is parameterized using field measured data (Figure 3) for this estimation. After compensating scaling error, mean biomass density declines by 3.6% and becomes 69.1 (Mg ha−1). This can result in a total forest biomass of 201.8 (106 Mg) for this species in the province (Table 2). The amount of reduced scaling error is 7.7 (106 Mg) using σ2 = 1721 m6 ha−2. The fine dotted line (Figure 3) denotes the curve from the power-law equation. It cannot properly fit the measurements especially when the volume exceeds 200 m3 ha−1.
The 3.6% overestimate is a remarkable percentage in this biomass estimation for the 2.92 (106 ha) forested area of Pinus yunnanensis. This error percentage is not necessarily for all species [22]. Pinus yunnanensis is a frigidly temperate coniferous forest species distributing from 23°1′ to 28°23′ N and 97°46′ to 105°54′ E in the province. It sustains drought and poor soils, widely adapts to extreme natural environmental conditions [41], and experiences fast growth on branch and leaves in middle-age [42,43]. This may result in its volume-to-biomass curve showing higher curvature than other species [23]. We found that the improved allometric VAR equation y = rxkeux is more feasible than the conventional power-law equation for Pinus yunnanensis. This finding implies that the stand condition is decisive for the curvature of the compensator in Equations (5) and (7), including parameters k and u. The values of parameters k (0.7623) and u (0.000536) are lower than the values for many species. This slowed growth under the poor site conditions is one of the possible reasons that the scaling error was reduced largely in the case study.

3.4. Uncertainty Analysis

As Equation 7 expressed, the estimate of variance σ2 is a key variable that affects the bias of the compensator. To evaluate the variance and the detailed situation is difficult for large-scale forests [38]. If the variance is unknown, it is necessary to assume a distribution of the stocking volume, for instance, uniform distribution. This would decrease the accuracy of compensation due to the bias of variance estimate. In the worst case, the variance bias could be Dmσ2, which can be derived as follows. We utilize σ2 = (xmax − xmin)2/12 for calculating the variance of uniform distribution, and rewrite Equation (8) as Dm = (μ − xmin)(xmaxμ). The relative error of variance becomes
(Dmσ2)/Dm = (2/3)(μ − xmin)/(xmaxμ),
As we have μ = (xmin + xmax)/2 for the uniform distribution, it leads to (Dmσ2)/Dm = 2/3. This implies that around 66% scaling error may be under-compensated by assuming a uniform distribution.
However, the worst condition is an extremity, in which the maximum variance Dm is only for the distribution having all stands located at either xmin or xmax. This stand structure does not exist in real forests. Practically, the distribution and forest structure are complex on tree height, stocking volume, and age group [39,44,45]. In our case study, we found that the distribution of stocking volume was close to an exponential curve on the area proportion (Table 2). Young stands have the largest area for Pinus yunnanensis [38]. Thus, the practical variance (σ2 = 1721) is less than the variance of uniform distribution (σ2 = 3158), which results in an estimate of 195.7 × 106 Mg. Comparing with the correct estimate (201.9 × 106 Mg) of total biomass on true distribution, an over-compensation could occur. This implies that the uniform distribution may not be a proper assumption for some regions in which large afforestation and reforestation have been promoted for long term.

4. Conclusions

A large-scale analysis of forest ecosystems is usually based on fine-scale allometric relationships or equations. Depending on the nonlinear strength of the equations, a scaling (or aggregation) error may result, where regional applications using the stand level volume-to-biomass relationship need to conduct error compensation to detect and remove the possible scaling error.
Our simulation indicates that, on average, approximately 94.8% of the scaling error can be reduced. We have shown that significant overestimates of biomass (~3.6%) can be removed from a provincial estimate for Pinus yunnanensis. However, various allometric equations will need different error compensators. We suggest assessing the effects of scaling errors on large-scale estimations of spatial and temporal distribution of forest biomass and carbon.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4907/10/7/602/s1, Supplementary Material S1: Supporting Information.

Author Contributions

C.Z. contributed to statistical analysis, simulation, and manuscript preparation. X.Z. conceived and designed experiment.

Funding

This work was supported by The National Key Research and Development Program of China (2017YFA0604401, 2016YFC0501101), Open Fund of State Key Laboratory of Remote Sensing Science (OFSLRSS201704) and Meteorology Scientific Research Fund in the Public Welfare of China (GYHY201506010).

Acknowledgments

The authors thank Joseph Carnegie and Jorge Potel for their helpful and valuable comments and suggestions. The authors thank Mingxia Yang for processing data and figures. The authors are also grateful to the three anonymous reviewers and the corresponding editors for their insightful comments and helpful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A1. Problem Background of Scaling Error

(a) Example

Supposition 1: A region includes two stands. The total area and volume are known for every stand (2 ha, 200 m3; 3 ha, 100 m3). If an allometric function is y = 2.7x0.7 e−0.0005x where x denotes stocking volume, we calculate the regional biomass as follows:
Y1 = 2[2.7(200/2)0.7e−0.0005(200/2)] + 3[2.7(100/3)0.7e−0.0005(100/3)] = 221.8 (Mg ha−1),
Supposition 2: Only total area and volume are known for the region (5 ha, 300 m3). The regional biomass may be calculated as
Y2 = 5[2.7(300/5)0.7e−0.0005(300/5)] = 230.1 (Mg ha−1),
This is an erroneous result, which shows an over-estimate by disregarding scaling error.

(b) Analysis

The improved allometric equation was presented as y = rxk eux [35] where r, k, and u are non-negative parameters. According to the actual situation, let x > 0 and k < 1.0. To apply the classical Jensen inequality for deriving scaling error, the function y must be a concave function. If letting g(x) = xkeux and y = rg(x), the second derivative of g(x) should be equal to or less than zero, i.e., u2x2 − 2kux – k + k2 ≤ 0. Then we can obtain the condition to apply the Jensen inequality: 0 ≤ x ≤ (k + k )/u.
The total biomass can be calculated through two different ways (Y1 and Y2):
Y1 = ∑aiyi = ∑ai(rxik euxi),
Y2 = Ak eux,
where ai represents area of ith stand, and xi denotes total stand stocking volume (both variables are unknown); A is region area including numerous stands, and μ represents mean total stocking volume of an entire region. A and μ are known variables that can be obtained easily from the forested area and total volume provided by the inventories.
For computing, the first way Y1 is precise and correct without scaling error since Y1 is added up from all stands. However, this expression is not applicable since ai and xi are unknown for each single stand in a large region. Therefore, the second way Y2 becomes the only choice to convert large-scale biomass directly from the mean stocking volume (m3 ha−1). We can rewrite Y1 and Y2 as follows:
Y 1 = a i y i = [ A ( a i A ) y i ] = [ Aw i   ( r x i k e u x i ) ] = A r [ w i ( x i k e u x i ) ] ,
Y2 = Ak eux = Ar(∑wixi)k euwixi,
where the random variable x is discrete with wi. According to convex inequalities, e.g., weighted Jensen inequality, we have ∑f(wixi) ≤ f(∑wixi) and Y1 ≤ Y2, i.e., an error (error = Y1 − Y2) exist. this error is caused by scaling forest biomass from stand-scale up to large-scale. To solve the error, let μg = ∑[wi(xikeuxi)], it leads to error = Ar (μgμk euμ). Hence the issue becomes to find the unknown μg using known μ and σ2. The step-by-step process is detailed in the results section.

Appendix A2. Derivation of Dm

Assuming n samples move from xmin to xmax one by one, the samples that are present at the point of xmin and xmax at any time are n − nz and nz respectively. Further, z is the percentage of the samples that have reached the point of xmax. Thus, we have
μ = [(nnz)xmin + nz xmax]/n,
μ(x2) = [(n − nz)(xmin)2 + nz(xmax)2]/n,
Simplifying them leads to:
μ = (1 − z)xmin + z xmax,
μ(x2) = (1 − z)(xmin)2 + z(xmax)2,
Then we convert Equations (A9)–(A11):
z = (μ − xmin)/(xmax − xmin),
Substitute Equations (A10) and (A11) into the variance formula σ2 = μ(x2)μ2, we have:
σ2 = −(μ − xmin)2 + (μ − xmin)(xmax − xmin),
Thus, the maximum variance is expressed as:
Dm = −(μ − xmin)2 + (μ − xmin)(xmax − xmin),
where Dm is the boundary of variance. represents the regional stocking volume. This is for any distribution where xmin and xmax are known.

References

  1. Odum, E.P.; Barrett, G.W. Fundamentals of Ecology, 5th ed.; Thomoson Brooks/Cole: Belmont, CA, USA, 2005. [Google Scholar]
  2. Huxley, J.S.; Teissier, G. Terminology of Relative Growth. Nature 1936, 137, 780–781. [Google Scholar] [CrossRef]
  3. Anderson, B.G.; Busch, H.L. Allometry in antennal. Biol. Bull. 1941, 81, 119–126. [Google Scholar] [CrossRef]
  4. Cowin, S.C. The Specific Growth Rates of Tissues: A Review and a Re-Evaluation. J. Biomech. Eng. 2011, 133, 41001. [Google Scholar] [CrossRef] [PubMed]
  5. Kittredge, J. Estimation of the amount of foliage of trees and stands. J. For. 1944, 42, 905–912. [Google Scholar]
  6. Niklas, K.J. Plant allometry: Is there a grand unifying theory? Biol. Rev. 2004, 79, 871–889. [Google Scholar] [CrossRef] [PubMed]
  7. Jorgensen, S.E.; Bendoricchio, G. Fundamentals of Ecological Modelling, 3rd ed.; Elservier Science: Amsterdam, The Netherlands, 2001. [Google Scholar]
  8. Searle, E.B.; Chen, H.Y.H. Persistent and pervasive compositional shifts of western boreal forest plots in Canada. Glob. Chang. Biol. 2017, 23, 857–866. [Google Scholar] [CrossRef] [PubMed]
  9. Reich, P.B.; Luo, Y.; Bradford, J.B.; Poorter, H.; Perry, C.H.; Oleksyn, J. Temperature drives global patterns in forest biomass distribution in leaves, stems, and roots. Proc. Natl. Acad. Sci. USA 2014, 111, 13721–13726. [Google Scholar] [CrossRef] [Green Version]
  10. Michaletz, S.T.; Cheng, D.; Kerkhoff, A.J.; Enquist, B.J. Convergence of terrestrial plant production across global climate gradients. Nature 2014, 512, 39–43. [Google Scholar] [CrossRef]
  11. Chu, C.; Bartlett, M.; Wang, Y.; He, F.; Weiner, J.; Chave, J.; Sack, L. Does climate directly influence NPP globally? Glob. Chang. Biol. 2016, 22, 12–24. [Google Scholar] [CrossRef]
  12. Bettencourt, L.M.A.; Lobo, J.; Helbing, D.; Kühnert, C.; West, G.B. Growth, innovation, scaling, and the pace of life in cities. Proc. Natl. Acad. Sci. USA 2007, 104, 7301–7306. [Google Scholar] [CrossRef] [Green Version]
  13. Li, X.; Wang, X.; Zhang, J.; Wu, L. Allometric scaling, size distribution and pattern formation of natural cities. Palgrave Commun. 2015, 1, 1–11. [Google Scholar] [CrossRef]
  14. 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]
  15. Komiyama, A.; Ong, J.E.; Poungparn, S. Allometry, biomass, and productivity of mangrove forests: A review. Aquat. Bot. 2008, 89, 128–137. [Google Scholar] [CrossRef]
  16. Cannell, M.G.R. World Forest Biomass and Primary Production Data; Academic Press: London, UK, 1982. [Google Scholar]
  17. Usoltsev, V.A. Forest Biomass and Primary Production Database for Eurasia, 2nd ed.; CD-ve.; Ural State Forest Engineering University: Yekaterinburg, Russia, 2013. [Google Scholar]
  18. Luo, Y.; Wang, X.; Zhang, X.; Lu, F. Biomass and Its Allocation of Forest Ecosystems in China; Chinese Forestry Publishing House Press: Beijing, China, 2013. [Google Scholar]
  19. Duursma, R.A.; Robinson, A.P. Bias in the mean tree model as a consequence of Jensen’s inequality. For. Ecol. Manag. 2003, 186, 373–380. [Google Scholar] [CrossRef]
  20. Gertner, G. Prediction bias and response surface curvature. For. Sci. 1991, 37, 755–765. [Google Scholar]
  21. Gardner, R.H.; Cale, W.G.; O’Neill, R.V. Robust Analysis of Aggregation Error. Ecology 1982, 63, 1771–1779. [Google Scholar] [CrossRef]
  22. Sambakhe, D.; Fortin, M.; Renaud, J.P.; Deleuze, C.; Dreyfus, P.; Picard, N. Prediction bias induced by plot size in forest growth models. For. Sci. 2014, 60, 1050–1059. [Google Scholar] [CrossRef]
  23. Rastetter, E.B.; King, A.W.; Cosby, B.J.; Hornberger, G.M.; O’Neill, R.V.; Hobbie, J.E. Aggregating fine-scale ecological knowledge to model coarser-resolution attributes of ecosystems. Ecol. Appl. 1992, 2, 55–77. [Google Scholar] [CrossRef]
  24. Kangas, A. On the bias and variance in tree volume predictions due to model and measurement errors. Scand. J. For. Res. 1996, 11, 281–290. [Google Scholar] [CrossRef]
  25. Zhou, X.; Lei, X.; Peng, C.; Wang, W.; Zhou, C.; Liu, C.; Liu, Z. Correcting the overestimate of forest biomass carbon on the national scale. Methods Ecol. Evol. 2016, 7, 447–455. [Google Scholar] [CrossRef]
  26. Woodall, C.W.; Heath, L.S.; Domke, G.M.; Nichols, M.C. Methods and Equations for Estimating Aboveground Volume, Biomass, and Carbon for Trees in the U.S. Forest Inventory, 2010; Tech. Rep. NRS-88.; USDA, Forest Service, Northern Research Station: Newtown Square, PA, USA, 2011; p. p. 30.
  27. Lei, X.; Tang, M.; Lu, Y.; Hong, L.; Tian, D. Forest inventory in China: status and challenges. Int. For. Rev. 2009, 11, 52–63. [Google Scholar] [CrossRef]
  28. Eggleston, H.S.; Buendia, L.; Miwa, K.; Ngara, T.; Tanabe, K. 2006 IPCC Guidelines for National Greenhouse Gas Inventories. In Forest land; IGES: Kanagawa, Japan, 2006. [Google Scholar]
  29. Brown, S.L.; Lugo, A. Biomass of Tropical Forests: A New Estimate Based on Forest Volumes. Science 1984, 223, 1290–1293. [Google Scholar] [CrossRef] [PubMed]
  30. Smith, J.E.; Heath, L.S.; Jenkins, J.C. Forest volume-to-biomass models and estimates of mass for live and standing dead trees of U.S. forests. USDA For. Serv. 2003, 1–62. [Google Scholar] [CrossRef]
  31. Guo, Z.; Fang, J.; Pan, Y.; Birdsey, R. Inventory-based estimates of forest biomass carbon stocks in China: A comparison of three methods. For. Ecol. Manag. 2010, 259, 1225–1231. [Google Scholar] [CrossRef]
  32. Henry, M.; Bombelli, A.; Trotta, C.; Alessandrini, A.; Birigazzi, L.; Sola, G.; Vieilledent, G.; Santenoise, P.; Longuetaud, F.; Valentini, R.; et al. GlobAllomeTree: International platform for tree allometric equations to support volume, biomass and carbon assessment. IForest 2013, 6, 1–5. [Google Scholar] [CrossRef]
  33. Fang, J.; Chen, A.; Peng, C.; Zhao, S.; Ci, L. Changes in forest biomass carbon storage in China between 1949 and 1998. Science 2001, 292, 2320–2322. [Google Scholar] [CrossRef]
  34. Boudewyn, P.; Song, X.; Magnussen, S.; Gillis, M.D. Model-Based, Volume-to-Biomass Conversion for Forested and Vegetated Land in Canada; Natural Resources Canada, Canadian Forest Service Pacific Forestry Centre Natural: Victoria, BC, Canada, 2007; p. 124. [Google Scholar]
  35. Ruark, G.; Bockheim, J.G.; Martin, G.L. Comparison of constant and variable allometric ratios for estimating Populus tremuloides biomass. For. Sci. 1987, 33, 294–300. [Google Scholar]
  36. Geron, C.D. Comparison of constant and variable allometric ratios for predicting foliar biomass of various tree genera. Can. J. For. Res 1988, 18, 1298–1340. [Google Scholar] [CrossRef]
  37. Berton, A.J.; Pregitzer, K.S.; Reed, D.D. Leaf area and foliar biomass relationships in northern hardwood forests located along an 800 km acid deposition gradient. For. Sci. 1991, 37, 1041–1059. [Google Scholar]
  38. Mather, A.S. Assessing the world’s forests. Glob. Environ. Chang. 2005, 15, 267–280. [Google Scholar] [CrossRef]
  39. Deng, X.; Huang, B.; Wen, Q.; Hua, C.; Tao, J.; Zheng, J. Dynamic of Pinus yunnanensis Forest Resources in Yunnan. J. Nat. Resour. 2014, 29, 1411–1419, (In Chinese with English abstract). [Google Scholar]
  40. Chinese Ministry of Forestry. Forest Resource Report of China for Periods 2003–2007; Department of Forest Resource and Management, Chinese Ministry of Forestry: Beijing, China, 2008. [Google Scholar]
  41. Lu, H. Preliminary estimation of forest carbon reserves of Pinus yunnanensis. For. Invent. Plan. 2010, 35, 91–93. [Google Scholar]
  42. Li, X. A review of researches on Pinus yunnanensis. J. Sichuan Agric. Univ. 1995, 13, 309–314. [Google Scholar]
  43. Wang, B.; Mao, J.F.; Zhao, W.; Wang, X.R. Impact of Geography and Climate on the Genetic Differentiation of the Subtropical Pine Pinus yunnanensis. PLoS ONE 2013, 8, e67345. [Google Scholar] [CrossRef] [PubMed]
  44. Dan, C.; Wu, Z. Studies on the biomass of Pinus yunnanensis forest. Acta Bot. Yunnanica 1991, 13, 59–64. [Google Scholar]
  45. Matasci, G.; Hermosilla, T.; Wulder, M.A.; White, J.C.; Coops, N.C.; Hobart, G.W.; Bolton, D.K.; Tompalski, P.; Bater, C.W. Three decades of forest structural dynamics over Canada’s forested ecosystems using Landsat time-series and lidar plots. Remote Sens. Environ. 2018, 216, 697–714. [Google Scholar] [CrossRef]
Figure 1. The percentage of scaling error compensated for Pinus yunnanensis Franch. η is the error compensation percentage given by Equation (7), in which this example uses parameters r = 2.7468, k = 0.7623 and, u = 0.000536. μ is the weighted mean of xi (stand stocking volume), and is also the mean regional volume. It is assumed to range from xmin (10 m3 ha1) to xmax (250 m3 ha1). σ2 is the variance of xi.
Figure 1. The percentage of scaling error compensated for Pinus yunnanensis Franch. η is the error compensation percentage given by Equation (7), in which this example uses parameters r = 2.7468, k = 0.7623 and, u = 0.000536. μ is the weighted mean of xi (stand stocking volume), and is also the mean regional volume. It is assumed to range from xmin (10 m3 ha1) to xmax (250 m3 ha1). σ2 is the variance of xi.
Forests 10 00602 g001
Figure 2. The efficiency of scaling error compensator based on statistic simulation. (a) Comparison between true value and modelled value of total amount of forest biomass in the pseudo region consisting of 100 stands (n = 100). (b) Percentage of compensated scaling error depending on trials. It becomes convergent after trials approach 4000 times. Note: the curve indicates the compensated scaling error in percentage, which is averaged from the trials.
Figure 2. The efficiency of scaling error compensator based on statistic simulation. (a) Comparison between true value and modelled value of total amount of forest biomass in the pseudo region consisting of 100 stands (n = 100). (b) Percentage of compensated scaling error depending on trials. It becomes convergent after trials approach 4000 times. Note: the curve indicates the compensated scaling error in percentage, which is averaged from the trials.
Forests 10 00602 g002
Figure 3. Stand stocking volume versus stand biomass density for parameterizing the volume-to-biomass equation y = rxk eux. The parameters estimated by nonlinear regression were r = 2.7468, k = 0.7623 and u = 0.000536 (m−3 ha) for Pinus yunnanensis. The fine dotted line is a conventional power-law curve (y = rxk, r = 2.6526, k = 0.7742, R2 = 0.7428). Under the condition of the Jensen inequality 0 ≤ x ≤ (k + k )/u, the range of x must be 0 ≤ x ≤ 3051 (m−3 ha). In practice, x is not greater than 500 (m−3 ha). The condition is met.
Figure 3. Stand stocking volume versus stand biomass density for parameterizing the volume-to-biomass equation y = rxk eux. The parameters estimated by nonlinear regression were r = 2.7468, k = 0.7623 and u = 0.000536 (m−3 ha) for Pinus yunnanensis. The fine dotted line is a conventional power-law curve (y = rxk, r = 2.6526, k = 0.7742, R2 = 0.7428). Under the condition of the Jensen inequality 0 ≤ x ≤ (k + k )/u, the range of x must be 0 ≤ x ≤ 3051 (m−3 ha). In practice, x is not greater than 500 (m−3 ha). The condition is met.
Forests 10 00602 g003
Table 1. Symbols and their descriptions.
Table 1. Symbols and their descriptions.
SymbolDescriptionUnit
ATotal forested area in a region, A = ∑ai.ha
aiThe area of ith stand.ha
C j s The combination. C j s = j!/[s!(js)!]. j represents the total number of elements, and s is the number of elements being chosen at a time.-
DmThe maximum variance, which is a function of expectation.m6 ha−2
EThe mathematical expectation operator.-
g(x)The function of random variable x, g(x) = xkeux.-
iStand number, i = 1, 2, …, n.-
kParameter, 0 < k ≤ 1.-
nThe number of stands.-
qDerivative order, q is set up to 4 in this study.-
rParameter, r > 0.-
uParameter, u > 0.m−3 ha
wiWeight or the probability of occurrence of ith stand, wi = ai/A.-
xStocking volume of a fine-scale area, e.g., a stand.m3 ha−1
xiStand stocking volume of ith stand.m3 ha−1
xmaxThe maximum possible value of x in the forest inventory.m3 ha−1
xminThe minimum possible value of x in the forest inventory.m3 ha−1
Y1Regional total biomass ideally supposed to accumulate from all stands.Mg
Y2Regional total biomass calculated from A and μ, Y2 = Akeux.Mg
yiStand biomass density of ith stand or sample plot, yi = r xik eux.Mg ha−1
zIntermediate parameter (0~1), the percentage of the samples at the point of xmax to total samples.-
Φ(μ, σ2)The compensator of scaling-up error.Mg ha−1
η(μ, σ2)The compensation rate; η(μ, σ2) = Φ(μ, σ2)/(Y2/A) = Φ(μ, σ2)/(keux).-
μ (xs)The expectation of xs, μ (xs) = E(xs), s = 1, 2.-
μThe expectation of x. It denotes regional stocking volume. μ = ∑wixi = E(x).m3 ha−1
μgThe expectation of g(x), μg = E[g(x)].-
νssth-order central moment, νs = E[(x − μ)s]= ∑ C j s μ(xj) (−μ)sj, sigma from j = 0 to s, s = 1, 2.-
σ2Variance of random variable x.m6 ha−2
Table 2. A case study of biomass estimation for Pinus yunnanensis at the provincial scale.
Table 2. A case study of biomass estimation for Pinus yunnanensis at the provincial scale.
Forest Age GroupsSum or Ave.
UnitYoungMiddleNear-MatureMatureOver-Mature
Total forested area A (data)106 ha1.0510.9550.4800.3410.0962.92
Total volume V (data)106 m341.1868.6245.6244.5622.45222.4
Area proportion (p)-0.3600.3270.1640.1170.0331.0
Mean stocking volume x and μm3 ha−139.271.995.1130.8233.876.1
p(x − 76.1)2 for calculating *m6 ha−2490.35.959.0348.7817.31721
Mean biomass Y Mg ha−1 71.7
Total biomass (71.7 × A)106 Mg 209.5
Corrected YMg ha−1 69.1
Removable scaling error % 3.6
Total biomass (69.1 × A)106 Mg 201.8
* Variance is calculated based on stocking volume of each age group (data of A and V refer to the reference [38]). Biomass density is calculated using y = rxkeux; parameters refer to the caption of Figure 3. Scaling error is calculated based on Equation (7).

Share and Cite

MDPI and ACS Style

Zhou, C.; Zhou, X. Removing the Scaling Error Caused by Allometric Modelling in Forest Biomass Estimation at Large Scales. Forests 2019, 10, 602. https://doi.org/10.3390/f10070602

AMA Style

Zhou C, Zhou X. Removing the Scaling Error Caused by Allometric Modelling in Forest Biomass Estimation at Large Scales. Forests. 2019; 10(7):602. https://doi.org/10.3390/f10070602

Chicago/Turabian Style

Zhou, Carl, and Xiaolu Zhou. 2019. "Removing the Scaling Error Caused by Allometric Modelling in Forest Biomass Estimation at Large Scales" Forests 10, no. 7: 602. https://doi.org/10.3390/f10070602

APA Style

Zhou, C., & Zhou, X. (2019). Removing the Scaling Error Caused by Allometric Modelling in Forest Biomass Estimation at Large Scales. Forests, 10(7), 602. https://doi.org/10.3390/f10070602

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