A Review of Sampling and Modeling Techniques for Forest Biomass Inventory

Heng Wu 1,2信封 纯色填充 and Hui Xu 1,*信封 纯色填充

1    College of Forestry, Southwest Forestry University, Kunming 650224, China

2    Southwest Survey and Planning Institute, National Forestry and Grassland Administration, Kunming 650031, China

*Author to whom correspondence should be addressed.

A&R 2023, Vol. 1, No. 1, 0002; https://doi.org/10.59978/ar01010002

Received: 29 January 2023; Revised: 22 March 2023; Accepted: 17 May 2023; Published: 24 May 2023

Copyright © 2023
This is an Open Access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC BY 4.0) (
https://creativecommons.org/licenses/by/4.0/)

Abstract: Forest biomass is the energy base and material source of forest ecosystem cycle, which is expressed by the dry matter weight or energy accumulated per unit area and time. It is also an important index to study the structure and function of forest ecosystem, and is the premise and basis of scientific management of forest ecosystem. In this paper, the concept, development history, and research status of forest biomass were reviewed. The sampling methods, model construction methods of forest biomass survey were analyzed. Finally, the research prospects and summaries of key technologies of forest biomass inventory and monitoring were put forward.

Keywords: biomass sampling; estimation model; remote sensing; biomass inventory

1. History and Methodologies

1.1. History

According to the time line clustering view of forest biomass studies, the first biomass studies were conducted by Ebermayer (1876) in Germany to measure the amount of branches and leaves and the weight of wood in several forests, and the results have been cited by geochemists in calculating the chemical elements in the biosphere for more than 50 years. Before the 1950s, the study of forest biomass received little attention (Boysen, 1910; Kitterge, 1944). It was not until the 1950s that researchers from various countries around the world began to pay attention to the study of forest biomass. Researchers in Japan and the United States successively carried out the study of forest productivity, and began to conduct actual investigation and data collection on the biomass and productivity of major forest ecosystems. In the early 1970s, with the implementation of the International Geosphere and Biosphere Program (IGBP), the study of vegetation biomass and productivity introduced the perspective of ecosystem, grasped the process of ecosystem material production from an overall level and combined with environmental factors, and greatly developed the study of forest biomass and productivity.

After the 1980s, some scholars established regression equations of biomass using easy measurement factors of stand to study the biomass in different areas (Brandeis et al., 2006; Montes et al., 2000; Nascimento & Laurance, 2002). Others have studied the distribution patterns of forest structure and biomass in Northern Europe and America, and the dynamics of forest biomass changes after disturbance (de Wit et al., 2006; Giese et al., 2003;  Kauffman et al., 2003; Lehtonen et al., 2004). With the mature application of satellite remote sensing technology in geography, some scholars have used the TM, ETM+ remote sensing data (Dong et al., 2003; Labrecque et al., 2006; Suganuma et al., 2006) and satellite Radar data (Austin et al., 2003; Hyde et al., 2007; Lucas et al., 2006) for studying forest biomass in a different area. With the accumulation of time, the remote sensing data volume is increasing, which brings new opportunities and challenges to the study of biomass estimation by remote sensing. the initial research on forest biomass started from the forest ecosystem and forest succession, and a large number of studies on forest biomass and carbon density were carried out from the 1990s to 2000. In the first decade of the new century, a large number of scholars carried out research on remote sensing and carbon sequestration, which are closely related to continuous inventory and management of forest resources, as shown in Figure 1.

Figure 1. Time line clustering view of forest biomass studies.

1.2. Methodologies

The traditional methods of biomass inventory and estimation are still dominant. These traditional methods could be divided into carbon dioxide balance method, micrometeorological field method, direct harvest method and expansion factor method (Wu et al., 2023). The direct harvest method is the most accurate survey method and the most practical method for terrestrial forests. Direct harvest method can be divided into average wood method, clear cutting method, and relative growth method. The traditional methods of biomass investigation and estimation are heavy in workload, complicated in process, poor in representativeness, and have not formed a system of measurement techniques. Therefore, the traditional methods cannot timely reflect the quick changes of macroscopical ecosystems and ecological environment conditions, which could not meet the practical requirements. With the development of “3S” (GIS, RS, and GPS) technology, studies on vegetation productivity and biomass based on remote sensing technology have developed from the traditional ground measurement on a small scale and two-dimensional scale to the estimation of a large scale and multidimensional space-time, so that forest biomass at different spatial scales from stand to region can be estimated quickly (Hill et al., 2018), accurately.

2. Sampling Techniques

Sampling is a method based on probability theory. The random selection of samples can ensure the representativeness of samples, avoid human interference and deviation, and estimate the sampling error (Wu et al., 2023). Different sampling methods should be used for different purposes (Hou et al., 2021). In practical problems, a specific sampling scheme is mostly composed of a variety of basic sampling methods (Bagaram & Tóth, 2020). The methods of systematic sampling, stratified sampling and random sampling are used in biomass inventory at various levels. As a classical sampling method, equal probability sampling has been widely used (Hawbaker et al., 2009), and methods such as unequal probability sampling, remote sensing based sampling and sparse population distribution sampling (Lei & Tang, 2007; Sterba, 2009) are more targeted in biomass inventory. The study of forest biomass sampling mainly focuses on two keywords, measurement and estimation (Perez-Cruzado et al., 2021), which are closely related to forest carbon storage. The study of forest biomass sampling is gradually extended to survey design and estimation model, as shown in Figure 2.

Figure 2. Cluster view of forest biomass sampling studies.

2.1. Equal Probability Sampling

2.1.1. Random Sampling

The sample units obtained by random sampling are scattered, which is not conducive to the development of actual forest resources investigation (Yu, 1974), and the accuracy depends on the number of sample units (Jin & Zhao, 2001; Meng et al., 1995). When the standard deviation of biomass between the sample units of systematic sampling is large, even if the samples are organized according to the method of systematic sampling, the estimated result is closer to the random sampling. To ensure the accuracy of sampling inventory results, it is necessary to increase the sample units. At present, the most common approach to estimating provincial forest biomass is by using data from continuous forest resource inventory plots data (Qin et al., 2017). However, the biomass estimation of this method is summarized from the level of individual trees to the level of plot, and then estimated the total biomass of the region. In the process, there are a large number of errors from different sources (McRoberts et al., 2013) and uncertainties, such as measurement errors, model uncertainties, and young trees that do not measure, which could lead to underestimation of forest biomass data (Poudel et al., 2015). The error sources of forest biomass estimation are diverse, and there are interaction and error conduction effects among them, so scientific sampling design is particularly important (Montesano et al., 2015).

2.1.2. Stratified Sampling

Compared with random sampling, stratified sampling often has significant efficiency for estimation. The effect of stratified sampling depends on the accuracy of prior information. The key point is that according to the classification of attribute characteristics, the variance between sample units in the same layer should be as small as possible, while the variance between sample units in different layers should be as large as possible. When the standard deviation of total biomass is large and the standard deviation between subpopulations is small, the estimation accuracy of using stratified sampling is greater than that of the random sampling (Zeng et al., 1995). Stratified sampling is often used to investigate forest volume in forest resource planning and design investigation (Yang, 1993). Cluster sampling divides the survey population into disjoint groups and then surveys the whole. By increasing the number of investigation units of adjacent samples to improve the accuracy and sampling efficiency and has been widely used in some forest resource surveys distributed in clusters (Li et al., 2019; Shi, 2012).

2.1.3. Cluster Sampling

Since adaptive sampling design was proposed in the 1990s (Thompson, 1991), sample organization method, estimation method and practical application have all been improved and perfected. Since traditional sampling methods do not consider the inconsistent contributions of different clusters, traditional methods such as systematic sampling may lead to the reduction of sampling efficiency and estimation accuracy (Gao & Gao, 2018; Hua et al., 2014; Huang., 2018; Zeng et al., 2018). Adaptive sampling design for sparse distribution population relying on the correlation between sample units can effectively reduce the number of sampling units (Hero et al., 2013; Lin et al., 2009); Traditional cluster sampling for aggregate distribution may lead to estimation bias (Xiao, 2004). The combination method of applying the results of biomass spatial distribution pattern analysis to adaptive sampling is unequal probability sampling (Holmberg & Lundevaller, 2015). Different distribution patterns should have different sampling probabilities. Adaptive sampling combined with unequal probability sampling can fully consider the spatial distribution difference of biomass, the status difference of different clusters, to carry out adaptive unequal probability sampling.

2.2. Unequal Probability Sampling

The advantage of unequal probability sampling is to improve the estimation accuracy and reduce the sampling error. The process of preparing the sampling frame is more complicated. PPES (sampling with probability proportional to an estimate of size) and PPP (probability proportional prediction) in forest resource inventory is the practical application of unequal probability sampling theory. In the 1970s, the forestry inventory began to use an unequal probability sampling design (Shi et al., 2009). Angular gauge for tree measurement is a typical application of unequal probability sampling in forest resource investigation. Different sample organization forms and unequal probability sampling combine to form different sampling design schemes (Good et al., 2001; Zhou & Sun, 2004). The sample sizes in cluster sampling are often different, the sampling probabilities of groups of different sizes are specified, and sampling with unequal probabilities can ensure the accuracy of sampling estimation under the requirement of fewer sampling units (Peng, 1998). Unequal probability sampling is more complicated than equal probability sampling in sampling frame design and sampling probability determination, but it can effectively improve the sampling efficiency (Li, 2000).

2.3. Other Sampling Methods

Unequal probability sampling is applicable to situations where the status of sampling units in the population is inconsistent or the units surveyed are inconsistent with those of the sampled population. The prerequisite of unequal probability sampling is that the sampling probability of sample units can be determined by the auxiliary information of each unit. Based on the prior information of the distribution pattern of forest biomass, the original sampling proportion was set, and the sampling adaptation was conducted according to the sampling probability of neighborhood sample units, so as to improve sampling efficiency. The Ripley’s K(d) analysis, aggregation index analysis, nearest neighbor analysis and spatial autocorrelation analysis are the main spatial pattern analysis methods (Jiang et al., 2009). The nearest neighbor analysis selects and measures the distance between each base unit and its nearest base unit in turn, then calculates the mean value of the nearest neighbor distance of all base units in the region, and compares it with the expected mean value of random distribution. Aggregation index is used to calculate and analyze the distance between adjacent sample units to describe the spatial distribution and aggregation state of biomass. The spatial autocorrelation analysis among sample units can effectively reduce sample redundancy and reduce investigation cost (Zhao et al., 2022).

3. Modeling Techniques

Stand growth model refers to a mathematical function or a group of mathematical functions describing the relationship between stand growth and site conditions, which is used to estimate the development process of a stand under certain conditions (Gao et al., 2017). The main uses of model include updating forest biomass data, evaluating the benefits of different forest management measures, evaluating the impact of disturbance activities on forest ecosystem, predicting the yield of forest sustainable management. Models can be classified into various types according to their purpose of use, model structure, and object of reflection. The correlation between forest carbon storage and aboveground biomass is still the main research of forest biomass model. Random forest and support vector machines are widely used in forest biomass model research, and forest measurement is combined, as shown in Figure 3.

Figure 3. Cluster view of forest biomass modeling studies.

3.1. Empirical Model

3.1.1. Tree Level

The biomass model of individual trees is a kind of model based on simulating the dry matter weight of each component (trunk, branch, leaf, bark, root, etc.) of each tree in the stand. More than 3000 biomass models have been established globally, involving more than 100 tree species (Basuki et al., 2009; Case & Hall, 2008; Chojnacky, 2002; Jenkins et al., 2003; Muukkonen, 2006). In the modeling process, it is necessary to measure the biomass of a certain number of sample trees as the basic data. Once the model is established, the continuous forest inventory data can be used to estimate the biomass of the whole stand in the same kind of stands with a precision. Therefore, the research of biomass model of individual trees has always been a hot spot (Kleinn et al., 2020).

According to the number of independent variables, tree level biomass models can be divided into unitary, binary and multivariate biomass models (Zeng & Tang, 2010). Models can be divided into linear, nonlinear and polynomial models, according to the form of tree level biomass models. Parameter estimation methods for biomass model construction of tree level biomass model mainly include traditional regression method (Basuki et al., 2009; Vallet et al., 2006), nonlinear likelihood uncorrelated regression method (Bi et al., 2004), linear or nonlinear joint estimation method (Tang et al., 2000), dummy variable method (Fu et al., 2012), mixed effects model method (Fehrmann et al., 2008), measurement error model method (Zeng & Tang, 2012), spatial regression method (Ou et al., 2014), and so forth.

In the research and application of biomass models, the component model is incompatible with the aggregate model (Parresol, 1999; Xu & Zhang, 2002). Tang et al. (2000) proposed a compatible tree level biomass model and estimation method combined with the current forest volume resources inventory method, namely the nonlinear joint estimation method. Compatible biomass models mainly include two types, namely the biomass model compatible with volume, and the biomass model compatible with total and component.

3.1.2. Stand Level

With the proposal of the IGBP, the results of previous ecological research systems and biomass data at stand level have been extended to landscape, regional and even global scales in the study of global climate change (Xu & Zhang, 2002). According to the independent variables, stand level biomass model can be divided into the model based on stand factor and the model based on volume. The biomass model based on volume can be subdivided into biomass expansion factor (BEF) and continuous biomass conversion functions (CBCF).

The dependent variable of the stand biomass model is usually the biomass per unit area of each organ of the stand, while the independent variable is the total basal area, dominant height, and average height of the stand to construct the stand biomass model of each dimension (Luo et al., 2009). Because the measurement of biomass at stand level rarely uses clear-cutting to obtain measured data, mostly obtains the biomass data based on the calculation of biomass model of individual trees or the measurement of standard trees (Dong & Li, 2016). Uncertainty analysis in the process of conversion from individual tree to stand level is also an important aspect of stand biomass research (Qin et al., 2017).

Scholars began to study to improve the accuracy of forest biomass estimation, and thus proposed a series of research methods (Fang, Chen, et al., 2001). Among them, BEF is a method to obtain the total biomass of stand by multiplying the average biomass of stand by the total volume of the forest types (Kauppi et al., 1992). The shortcomings of the BEF method are mainly reflected in the conversion factors, such as wood density and conversion ratio of total biomass and aboveground biomass, which are taken as constants. Fang, Wang et al. (1998) pointed out that stand biomass and volume were related to forest type, age, site conditions, stand density and other factors, and that the use of constant biomass conversion factor could not accurately estimate forest biomass. The continuous function method of conversion factors changes the constant average conversion factor into an age-graded conversion factor to achieve a more accurate estimation of forest biomass at the national or regional scale (Edgar et al., 2019; Kauppi et al., 1992; Turner et al., 1995). The conversion model between biomass and volume is a hot topic in recent years, which needs to be further verified from the region and tree species, and the model relationship should be established comprehensively and systematically.

3.2. Nonparametric Model

Traditional statistical regression methods could not effectively describe the complex nonlinear relationship between forest biomass and measured data in certain situation, as well as practical problems such as high dimension, and the derived relationship is usually only applicable to this region. Although the accuracy can be improved by the learning method, the “black box” operation only shows their complex action process through the simulation of some training data sets, which is difficult to reflect the mechanism between biomass and remote sensing parameters. In order to improve the nonlinear prediction ability of biomass model, data mining and machine learning methods were applied to forest biomass estimation, including decision tree, K-NN method, support vector machine and artificial neural network.

3.2.1. Artificial Neural Network

Artificial neural network (ANN) is based on spectral information, vegetation index and texture characteristics as input variables of neural network, forest aboveground biomass of sample plot survey as output variables, select part of the sample data input neural network system for training to obtain a model algorithm, and then estimate the forest aboveground biomass according to the model algorithm (Foody et al., 2003; Xu et al., 2011).

3.2.2. Decision Tree

Decision tree (DT) is a method to approximate discrete value function, which can be regarded as a tree prediction model. The basic algorithms include random forest and gradient lifting decision tree. Decision tree integration methods can remove noise interference well, training complexity is low, prediction is accurate, and the model is easy to display, but there may be problems with over-fitting of training data (Carreiras, et al., 2012).

3.2.3. K-NN Method

K-nearest neighbor (K-NN) classification algorithm is also known as the reference plot method. The forest aboveground biomass value of a certain pixel of remote sensing image is obtained by weighted value of K measured sample points closest to the pixel in feature space (Chirici et al., 2008; Tuominen, et al., 2010), and then monitored forest aboveground biomass according to the sample plot data. The K-NN method can estimate the forest biomass and maintain the heterogeneity and similarity characteristics of the spatial distribution of carbon density, but its estimation results are often higher than those estimated by using the sample plot data (Labrecque et al., 2006).

3.2.4. Support Vector Machine

The principle of support vector machine (SVM) can be summarized, the nonlinear transformation defined by the inner product function is used to transform the input space into a high-dimensional space, and then the optimal classification surface is obtained in this space. Each intermediate node corresponds to a support vector, and the output is a combination of nodes (Zhang, 2000). Support vector regression machine is a special form of SVM and a kernel theory of regression analysis and equation approximation (Englhart et al., 2011). It overcomes the defects of insufficient data and over-learning of traditional prediction methods and has unique advantages in solving small sample and high-dimensional problems. However, improper selection of kernel function would cause errors in estimation results.

3.3. Remote Sensing Inversion Model

3.3.1. Based on Tree Structure Parameters

As for the structural parameters of trees, it is difficult to obtain the DBH by satellite remote sensing because the trunk is heavily shielded by the crown. As the vertical structure parameter of trees, height can be accurately measured by active lidar remote sensing. The height can be calculated based on the time interval between the echo signal received by the lidar and the echo signal from the ground. Then, the biomasses were calculated by the tree height or stand average height. Many space-borne full-waveform tree height inversion models have been developed (Lefsky et al., 2005).

3.3.2. Based on Vegetation Feature Index

Spectral information analysis of optical remote sensing images of forests can reflect the biophysical characteristics of forests. A variety of vegetation indices can be applied for modeling based on forest optical remote sensing data. Commonly used vegetation index is usually a linear or nonlinear combination of spectral reflectance ratio of two or more bands (Zeng et al., 2022), which can compress multi-dimensional spectral information into one index channel. By statistical regression of various vegetation indices, structural parameters such as canopy density and leaf area index can be further estimated. Based on the statistical regression method, the relationship between the spectral reflectance provided by optical remote sensing images and the vegetation biomass can be established, and then the regional vegetation biomass can be estimated.

3.3.3. Based on Physical Mechanism

The mechanism model (or process model) is used to describe the vegetation growth process at different spatial-temporal scales, such as photosynthesis, respiration, decomposition and oxygen cycle, etc. It simulates the process of solar energy conversion into chemical energy, and the process of plant body and soil water loss accompanied by canopy evapotranspiration and photosynthesis according to the principles of plant physiology, to achieve the estimation of forest productivity. The mechanism model is incorporated into models of global change and nutrient cycling, with biomass being only one of the model’s output variables, taking CENTURY, CARAIB, TEM and so forth as Examples. The disadvantage of these models are often too complicated and need more input variables, the application of the model often depends on the quality of the data. The mechanism model emphasizes more on the description of various action processes within the ecosystem, and the estimation results are generally more reliable (Xu & Cao, 2006).

4. Research Prospects

The biomass inventory sampling should make full use of the existing spatial data for distribution pattern analysis, and carry out sampling design based on the results of spatial distribution pattern analysis (Wu & Xu, 2021). Spatial sampling technique can carry out with unequal probability sampling based on the spatial distribution pattern of forest biomass, which can effectively improve the investigation efficiency (Li et al., 2009; Wu et al., 2004). At the same time, in order to monitor large scale macroscopic changes and development trends of forest biomass, remote sensing data were used as sampling basic data and auxiliary survey data, and the method of aerial sampling was used to arrange sample plots on remote sensing maps to estimate forest biomass (Liu, 2001; Liu, 2016) to meet the practical requirements of biomass estimation at different regional scales and spatial distribution characteristics (Hetzer et al., 2020; Zhu et al., 2020).

In recent years, scholars have presented the biomass models of various tree species, and carried out the study of forest biomass at multiple scales, such as individual, population, community, ecosystem, region and biosphere (McRoberts, 2001). It is necessary to study the biomass of individual tree deeply, including different geographical provenances, different development stages, and different natural zones, so as to establish a weight index model of biomass to achieve a more accurate estimation of biomass in different stand types. At present, there is still a lack of research on multi-level stand models. Only the biomass of standing trees with a certain diameter in forest was estimated, while the biomass of trees with a smaller diameter, understory shrubs and herbaceous plants was ignored. Therefore, the relationship between the total stand biomass and the biomass of living trees was clarified, and the biomass model of understory vegetation was established to solve the problem that the existing model neglected the biomass of smaller trees, understory shrubs, and herbs. In the future, model development, the effects of forest biological factors and non-biological factors on forest biomass should be considered comprehensively, especially the effects of stand volume, age, and climate factors on forest biomass estimation.

Remote sensing has the characteristics of macroscopical, comprehensive, dynamic, rapid and repeatable, and its band information has a certain correlation with forest biomass, so it has become the main method for estimating regional forest biomass (Wirasatriya et al., 2022). Each remote sensing data has certain limitations in spatial, spectral and radiative resolution, which affects the ability of remote sensing technology to estimate forest biomass. These factors result in the instability of the accuracy of estimating ground forest biomass with different remote sensing data (Yu et al., 2022). Combining ecological factors, topographic factors, environmental factors, and remote sensing data to build a forest biomass estimation model with multi-source data can inhibit the influence of these factors. When the optical remote sensing data are used to estimate the areas with high forest biomass, the problem of remote sensing information saturation comes out. As a result, the changes of biomass cannot be accurately reflected, which becomes the bottleneck of forest biomass estimation by remote sensing method. Steininger (2000) found that LandSatTM images had a saturation problem when estimating biomass. Its saturation threshold was 15kg/m2. Lu (2005) also found this problem when estimating the Amazon basin in Brazil. Similar problems exist in estimating forest biomass with radar data (Wang et al., 2006). Multi-sensor data integration and estimation methods can solve this problem to some extent. Multi-source remote sensing image by multi-sensor to estimate forest biomass has become a developing trend to address the saturation threshold issue (Zhu et al., 2020). Combining remote sensing data of different sensors, different time, different spatial resolution different spectral resolution, and selecting the optimal information to estimate forest biomass is a problem that needs to be studied in data assimilation, as well as a problem faced by contemporary remote sensing development (Zeng et al., 2022). There is a certain gap between field sampling data and remote sensing image data, which is also a problem faced by remote sensing data assimilation.

5. Summary

The change of ecological environment caused by climate warming is becoming more and more obvious. It is affecting the pattern of ecological system and the sustainable development of human society. It has become a global environmental problem recognized by the international community. Forest is one of the main terrestrial ecosystems, which has an important carbon sink function. Forest carbon reserves account for about 80% of the above-ground part of the land and 40% of the underground part of the land, and the annual carbon fixed amount accounts for more than 2/3 of the whole terrestrial ecosystem. Forest in China will become a large carbon sink, which will play a positive role in mitigating the rise of global atmospheric carbon dioxide concentration.

It is a research hotspot and technical difficulty in the field of natural resources investigation and monitoring that efficiently and accurately monitors the annual dynamic change of forest biomass. National forest and grassland ecological comprehensive monitoring in 2021 and 2022 has optimized survey organization. The number of sample units to be investigated by equiprobability sampling is still quite large, under the specified accuracy and reliability. The estimation of forest biomass below the provincial level takes the regional total carbon storage as the control number, and divides them into small populations according to the principle of hierarchical control to produce the biomass data. Due to the lack of precision and the inability to sensitively reflect the dynamic changes of forest biomass at different scales, the application of the monitoring results in forest carbon sink accounting and other work is limited.

Based on the multi-stage sampling framework, researchers adopted the three-phase sampling method, combined with the characteristics of the regional distribution pattern of forest biomass, and carried out multi-scale unequal probability spatial stratified complex sampling design. The model sampling inference assisted by randomization inference was used to solve the problems of the limited sample size and small area estimation, data missing and measurement error under the specified precision, and the complicated sampling design uncertainty analysis and reliability evaluation were carried out. It is expected to form a set of complex sampling design and data inference technology for annual monitoring of forest biomass under hierarchical control, so as to meet the low-cost, rapid and accurate annual counting demand of forest biomass to be applied in practical work.

CRediT Author Statement: Heng Wu: Conceptualization, Methodology, Data curation, Writing – original draft, Software, Validation and Writing – review & editing; Hui Xu: Conceptualization, Methodology, Visualization, Investigation, Supervision and Writing – review & editing.

Data Availability Statement: Not applicable.

Funding: This research was funded by the National Natural Science Foundation of China (31560209 and 31760206) and the Academician Workstation of Yunnan Province of China (2018IC066).

Conflicts of Interest: The authors declare no conflict of interest.

Acknowledgments: We would like to acknowledge all the people who have contributed to this paper.

References

Austin, J. M., Mackey, B. G., & Van Niel, K. P. (2003). Estimating forest biomass using satellite radar: an exploratory study in a temperate Australian eucalyptus forest. Forest Ecology and Management, 176(1–3), 575–583. https://doi.org/10.1016/s0378-1127(02)00314-6 

Bagaram, M. B., & Tóth, S. F. (2020). Multistage sample average approximation for harvest scheduling under climate uncertainty. Forests, 11(11), 1230. https://doi.org/10.3390/f11111230

Basuki, T. M., van Laake, P. E., Skidmore, A. K., & Hussin, Y. A. (2009). Allometric equations for estimating the above-ground biomass in tropical lowland dipterocarp forests. Forest Ecology and Management, 257(8), 1684–1694. https://doi.org/10.1016/j.foreco.2009.01.027

Bi, H., Turner, J., & Lambert, M. J. (2004). Additive biomass equations for native eucalypt forest trees of temperate Australia. Trees, 18(4), 467–479. https://doi.org/10.1007/s00468-004-0333-z

Boysen, J. (1910). Studier over skovtraernes forhold til lyset Tidsskr. F Skorvaessen, 22, 1116.

Brandeis, T. J., Delaney, M., Parresol, B. R., & Royer, L. (2006). Development of equations for predicting Puerto Rican subtropical dry forest biomass and volume. Forest Ecology and Management, 233(1), 133–142. https://doi.org/10.1016/j.foreco.2006.06.012

Carreiras, J. M. B., Vasconcelos, M. J., & Lucas, R. M. (2012). Understanding the relationship between aboveground biomass and ALOS PALSAR data in the forests of Guinea-Bissau (West Africa). Remote Sensing of Environment, 121, 426–442. https://doi.org/10.1016/j.rse.2012.02.012

Case, B. S., & Hall, R. J. (2008). Assessing prediction errors of generalized tree biomass and volume equations for the boreal forest region of west-central Canada. Canadian Journal of Forest Research, 38(4), 878–889. https://doi.org/10.1139/x07-212

Chirici, G., Barbati, A., Corona, P., Marchetti, M., Travaglini, D., Maselli, F., & Bertini, R. (2008). Non-parametric and parametric methods using satellite images for estimating growing stock volume in alpine and Mediterranean forest ecosystems. Remote Sensing of Environment, 112(5), 2686–2700. https://doi.org/10.1016/j.rse.2008.01.002

Chojnacky, D. C. (2002). Allometric scaling theory applied to FIA biomass estimation. U.S. Department of Agriculture, Forest Service, North Central Research Station. https://www.nrs.fs.usda.gov/pubs/gtr/gtr_nc230/gtr_nc230_096.pdf

de Wit, H. A., Palosuo, T., Hylen, G., & Liski, J. (2006). A carbon budget of forest biomass and soils in southeast Norway calculated using a widely applicable method. Forest Ecology and Management, 225(1–3), 15–26. https://doi.org/10.1016/j.foreco.2005.12.023

Dong, J., Kaufmann, R. K., Myneni, R. B., Tucker, C. J., Kauppi, P. E., Liski, J., Buermann, W., Alexeyev, V., & Hughes, M. K. (2003). Remote sensing estimates of boreal and temperate forest woody biomass: Carbon pools, sources, and sinks. Remote Sensing of Environment, 84(3), 393–410. https://doi.org/10.1016/s0034-4257(02)00130-x

Dong, L., & Li, F. (2016). Comparison of three stand-level biomass estimation methods. Chinese Journal of Applied Ecology, 27(12), 3862–3870. https://doi.org/10.13287/j.1001-9332.201612.030

Ebermayer, E. (1876). Die gesammte lehre der waldstreu mit rücksicht auf die chemische statik des waldbaues. Unter zugrundlegung der in den Königl. Staatsforsten bayerns angestellten untersuchungen. Springer. https://doi.org/10.1007/978-3-642-91491-1

Edgar, C. B., Westfall, J. A., Klockow, P. A., Vogel, J. G., & Moore, G. W. (2019). Interpreting effects of multiple, large-scale disturbances using national forest inventory data: A case study of standing dead trees in east Texas, USA. Forest Ecology and Management, 437, 27–40. https://doi.org/10.1016/j.foreco.2019.01.027

Englhart, S., Keuck, V., & Siegert, F. (2011). Aboveground biomass retrieval in tropical forests—The potential of combined X- and L-band SAR data use. Remote Sensing of Environment, 115(5), 1260–1271. https://doi.org/10.1016/j.rse.2011.01.008

Fang, J., Chen, A., Peng, C., Zhao, S., & Ci, L. (2001). Changes in forest biomass carbon storage in China between 1949 and 1998. Science, 292(5525), 2320–2322. https://doi.org/10.1126/science.1058629

Fang, J., Wang, G., Liu, G., & Xu, S. (1998). Forest biomass of China: An estimate based on the biomass-volume relationship. Ecological Applications, 8(4), 1084–1091. https://doi.org/10.1890/1051-0761(1998)008[1084:Fbocae]2.0.Co;2

Fehrmann, L., Lehtonen, A., Kleinn, C., & Tomppo, E. (2008). Comparison of linear and mixed-effect regression models and a knearest neighbour approach for estimation of single-tree biomass. Canadian Journal of Forest Research, 38(1), 1–9. https://doi.org/10.1139/x07-119

Foody, G. M., Boyd, D. S., & Cutler, M. E. J. (2003). Predictive relations of tropical forest biomass from Landsat TM data and their transferability between regions. Remote Sensing of Environment, 85(4), 463–474. https://doi.org/10.1016/s0034-4257(03)00039-7

Fu, L., Zeng, W., Tang, S., Guo, Q., Li, Z., & Cheng, X. (2012). Using linear mixed model and dummy variable model approaches to construct compatible single-tree biomass equations at different scales—A case study for Masson pine in Southern China. Journal of Forest Science, 58(3), 101–115. https://doi.org/10.13275/j.cnki.lykxyj.2011.03.011

Gao, H., Dong, L., & Li, F. (2017). Modeling variation in crown profile with tree status and cardinal directions for planted larix olgensis henry trees in northeast. Forests, 8(5), 139. https://doi.org/10.3390/f8050139

Gao, Y., & Gao, C. (2018). Forest biomass extraction based on ZY-3 remote sensing image. Beijing Surveying and Mapping, 32(10), 1186–1191. https://doi.org/10.19580/j.cnki.1007-3000.2018.10.017

Giese, L. A. B., Aust, W. M., Kolka, R. K., & Trettin, C. C. (2003). Biomass and carbon pools of disturbed riparian forests. Forest Ecology and Management, 180(1–3), 493–508. https://doi.org/10.1016/s0378-1127(02)00644-8

Good, N. M., Paterson, M., Brack, C., & Mengersen, K. (2001). Estimating tree component biomass using variable probability sampling methods. Journal of Agricultural Biological and Environmental Statistics, 6(2), 258–267. https://doi.org/10.1198/108571101750524599

Hawbaker, T. J., Keuler, N. S., Lesak, A. A., Gobakken, T., Contrucci, K., & Radeloff, V. C. (2009). Improved estimates of forest vegetation structure and biomass with a LiDAR-optimized sampling design. Journal of Geophysical Research-Biogeosciences, 114. https://doi.org/10.1029/2008jg000870

Hero, J. M., Castley, J. G., Butler, S. A., & Lollback, G. W. (2013). Biomass estimation within an Australian eucalypt forest: Mesoscale spatial arrangement and the influence of sampling intensity. Forest Ecology and Management, 310, 547–554. https://doi.org/10.1016/j.foreco.2013.08.062

Hetzer, J., Huth, A., Wiegand, T., Dobner, H. J., & Fischer, R. (2020). An analysis of forest biomass sampling strategies across scales. Biogeosciences, 17(6), 1673–1683. https://doi.org/10.5194/bg-17-1673-2020

Hill, A., Mandallaz, D., & Langshausen, J. (2018). A double-sampling extension of the German national forest inventory for design—based small area estimation on forest district levels. Remote Sensing, 10(7), 1052. https://doi.org/10.3390/rs10071052

Holmberg, H., & Lundevaller, E. H. (2015). A test for robust detection of residual spatial autocorrelation with application to mortality rates in Sweden. Spatial Statistics, 14, 365–381. https://doi.org/10.1016/j.spasta.2015.07.001

Hou, Z., Domke, G. M., Russell, M. B., Coulston, J. W., Nelson, M. D., Xu, Q., & McRoberts, R. E. (2021). Updating annual state- and county-level forest inventory estimates with data assimilation and FIA data. Forest Ecology and Management, 483, 118777. https://doi.org/10.1016/j.foreco.2020.118777

Hua, W., Qiu, Y., Xu, B., Jiang, X., & Qiu, T. (2014). Estimated forest carbon sequestration of the Fujian province based on growth and biomass model. Journal of Southwest Forestry University, 34(6), 35–43.

Huang, Y. (2018). Estimating forest carbon storage of talin forest farm with forest inventory data, northeast China. Journal of North-East Forestry University, 46(5), 12–16,37. https://doi.org/10.13759/j.cnki.dlxb.2018.05.003

Hyde, P., Nelson, R., Kimes, D., & Levine, E. (2007). Exploring LIDAR-RaDAR synergy - predicting aboveground biomass in a southwestern ponderosa pine forest using LiDAR, SAR and InSAR. Remote Sensing of Environment, 106(1), 28–38. https://doi.org/10.1016/j.rse.2006.07.017

Jenkins, J., Chojnacky, D. C., Heath, L. S., & Birdsey, R. A. (2003). National scale biomass estimators for United States tree species. Forest Science, 49(1), 12–35.

Jiang, C., Wang, J., & Cao, Z. (2009). A review of geo-spatial sampling theory. Acta Geographica Sinica, 64(3), 368–380.

Jin, L., & Zhao, R. (2001). Determination method of optimum quadrat size for simple random sampling of pupa population of Pinus tabulaeformis. Liaoning Forestry Science and Technology, (2), 42–44.

Kauffman, J. B., Steele, M. D., Cummings, D. L., & Jaramillo, V. J. (2003). Biomass dynamics associated with deforestation, fire, and, conversion to cattle pasture in a Mexican tropical dry forest. Forest Ecology and Management, 176(1–3), 1–12. https://doi.org/10.1016/s0378-1127(02)00227-x

Kauppi, P. E., Mielikäinen, K., & Kuusela, K. (1992). Biomass and carbon budget of European forests, 1971 to 1990. Science, 256(5053), 70–74. https://doi.org/10.1126/science.256.5053.70

Kitterge, J. (1944). Estimation of amount of foliage of trees and shrubs. J. Forest, 42, 905–912.

Kleinn, C., Magnussen, S., Nolke, N., Magdon, P., Alvarez-Gonzalez, J. G., Fehrmann, L., & Perez-Cruzado, C. (2020). A comparison of four methods to map biomass from Landsat-TM and inventory data in western Newfoundland of field inventory estimation of aboveground biomass through an alternative view on plot biomass. Forest Ecosystems, 7(1), 57. https://doi.org/10.1186/s40663-020-00268-7

Labrecque, S., Fournier, R. A., Luther, J. E., & Piercey, D. (2006). A comparison of four methods to map biomass from Landsat-TM and inventory data in western Newfoundland. Forest Ecology and Management, 226(1–3), 129–144. https://doi.org/10.1016/j.foreco.2006.01.030

Lefsky, M. A., Harding, D. J., Keller, M., Cohen, W. B., Carabajal, C. C., Del Bom Espirito-Santo, F., Hunter, M. O., & de Oliveira, R. C., Jr. (2005). Estimates of forest canopy height and aboveground biomass using ICESat. Geophysical Research Letters, 32(22). https://doi.org/10.1029/2005gl023971

Lehtonen, A., Mäkipää, R., Heikkinen, J., Sievänen, R., & Liski, J. (2004). Application of factors (BEFs) for Scots pine, Norway spruce and birch according to stand age for boreal forests. Forest Ecology and Management, 188(1–3), 211–224.

Lei, Y., & Tang, S. (2007). Application of adaptive cluster sampling in multi-resource inventory. Scientia Silvae Sinicae, 43(11), 132–137.

Li, C., Chen, Q., & Tang, B. (2009). Large-scale forest resources monitoring by means of spatial sampling of hi-resolution remote sensing images. Forest Resources Management, 2, 106–110. https://doi.org/ 10.13466/j.cnki.lyzygl.2009.02.013

Li, J. (2000). Research on equal generalization of unequal probability problem in sampling. Statistical Research, 17(4), 54–57. https://doi.org/10.19343/j.cnki.11-1302/c.2000.04.011

Li, Y., Chen, Z., Lei, J., Chen, X., Yang, Q., & Wu, T. (2019). Study on spatial balance sampling of forest resources survey in Haikou. Forest Resources Management, (2), 47. https://doi.org/10.13466/j.cnki.lyzygl.2019.02.007

Lin, G., Wen, X., & She, G. (2009). Selection and analysis of auxiliary factors using sampling method of sub-compartment. Journal of Nanjing Forestry University (Natural Sciences Edition), 33(1), 121–123.

Liu, H. (2001). Sampling method with remote sensing for monitoring of cultivated land changes on large scale. Nongye Gongcheng Xuebao (Transactions of the Chinese Society of Agricultural Engineering), 17(2), 168–171.

Liu, L. (2016). Inversion of forest carbon density in Changning country based on GF-1 remote sensing image [Master’s thesis, Central South University of Forestry and Technology].

Lu, D. (2005). Aboveground biomass estimation using Landsat TM data in the Brazilian Amazon. International Journal of Remote Sensing, 26(12), 2509–2525. https://doi.org/10.1080/01431160500142145

Lucas, R. M., Cronin, N., Lee, A., Moghaddam, M., Witte, C., & Tickle, P. (2006). Empirical relationships between AIRSAR backscatter and LiDAR-derived Forest biomass, Queensland, Australia. Remote Sensing of Environment, 100(3), 407–425. https://doi.org/10.1016/j.rse.2005.10.019

Luo, Y., Zhang, X., Wang, X., Zhu, J., Hou, Z., & Zhang, Z. (2009). Forest biomass estimation methods and their prospects. Scientia Silvae Sinicae, 45(8), 129–134. https://doi.org/10.3321/j.issn:1001-7488.2009.08.023

McRoberts, R. E. (2001). Imputation and model-based updating techniques for annual forest inventories. Forest Science, 47(3), 322–330.

McRoberts, R. E., Naesset, E., & Gobakken, T. (2013). Inference for lidar-assisted estimation of forest growing stock volume. Remote Sensing of Environment, 128, 268–275. https://doi.org/10.1016/j.rse.2012.10.007

Meng, F., Xu, S., & Yang, X. (1995). Application of random sampling survey in non-productive consumption survey of forest resources. Forest Resources Management, (4), 35–37.

Montes, N., Gauquelin, T., Badri, W., Bertaudiere, V., & Zaoui, E. H. (2000). A non-destructive method for estimating above-ground forest biomass in threatened woodlands. Forest Ecology and Management, 130(1–3), 37–46. https://doi.org/10.1016/s0378-1127(99)00188-7

Montesano, P. M., Rosette, J., Sun, G., North, P., Nelson, R. F., Dubayah, R. O., Ranson, K. J., & Kharuk, V. (2015). The uncertainty of biomass estimates from modeled ICESat-2 returns across a boreal forest gradient. Remote Sensing of Environment, 158, 95–109. https://doi.org/10.1016/j.rse.2014.10.029

Muukkonen, P. (2006). Forest inventory-based large-scale forest biomass and carbon budget assessment: new enhanced methods and use of remote sensing for verification. Dissertationes Forestales, 2006(30). https://doi.org/10.14214/df.30

Nascimento, H. E. M., & Laurance, W. F. (2002). Total aboveground biomass in central Amazonian rainforests: a landscape-scale study. Forest Ecology and Management, 168(1–3), 311–321. https://doi.org/10.1016/s0378-1127(01)00749-6

Ou, G., Wang, J., Xiao, Y., & Xu, H. (2014). Modeling individual biomass of Pinus Kesaya var. langbianensis natural forests by geographically weighted regression. Forest Research, 27(2), 213–218. https://doi.org/10.13275/j.cnki.lykxyj.2014.02.012

Parresol, B. R. (1999). Assessing tree and stand biomass: A review with examples and critical comparisons. Forest Science, 45(4), 573–593. 

Peng, N. (1998). Discussion on sampling technology with unequal probability. Statistical Research, 15(3), 3–5.

Perez-Cruzado, C., Kleinn, C., Magdon, P., Alvarez-Gonzalez, J. G., Magnussen, S., Fehrmann, L., & Noelke, N. (2021). The horizontal distribution of branch biomass in European beech: A model based on measurements and tls based proxies. Remote Sensing, 13(5), 1041. https://doi.org/10.3390/rs13051041

Poudel, K. P., Temesgen, H., & Gray, A. N. (2015). Evaluation of sampling strategies to estimate crown biomass. Forest Ecosystems, 2. https://doi.org/10.1186/s40663-014-0025-0

Qin, L., Zhang, M., & Zhong, S. (2017). Model uncertainty in forest biomass estimation. Acta Ecologica Sinica, 37(23), 7912–7919. https://doi.org/10.5846/stxb201609281973

Shi, J. (2012). The study on spatial distributions and adaptive cluster sampling for Dongzhaigang mangrove in Hainan [Doctoral dissertation, Beijing Forestry University].

Shi, J., Lei, Y., & Zhao, T. (2009). Progress in sampling technology and methodology in forest inventory. Forest Research, 22(1), 101–108.

Steininger, M. K. (2000). Satellite estimation of tropical secondary forest above-ground biomass: data from Brazil and Bolivia. International Journal of Remote Sensing, 21(6–7), 1139–1157. https://doi.org/10.1080/014311600210119

Sterba, S. K. (2009). Alternative model-based and design-based frameworks for inference from samples to populations: From polarization to integration. Multivariate Behavioral Research, 44(6), 711–740. https://doi.org/10.1080/00273170903333574

Suganuma, H., Abe, Y., Taniguchi, M., Tanouchi, H., Utsugi, H., Kojima, T., & Yamada, K. (2006). Stand biomass estimation method by canopy coverage for application to remote sensing in an and area of Western Australia. Forest Ecology and Management, 222(1–3), 75–87. https://doi.org/10.1016/j.foreco.2005.10.014

Tang, S., Zhang, H., & Xu, H. (2000). Study on establish and estimate method of compatible biomass model. Scientia Silvae Sinicae, 36(S1), 19–27. 

Thompson, S. K. (1991). Adaptive cluster sampling: Designs with primary and secondary units. Biometrics, 47(3), 1103. https://doi.org/10.2307/2532662

Tuominen, S., Eerikainen, K., Schibalski, A., Haakana, M., & Lehtonen, A. (2010). Mapping biomass variables with a multi-source forest inventory technique. Silva Fennica, 44(1), 109–119, 458. https://doi.org/10.14214/sf.458

Turner, D. A., Koerper, G. J., Harmon, M. E., & Lee, J. E. (1995). A carbon budget for forests of the conterminous United States. Ecological Applications, 5(2), 421–436. https://doi.org/10.2307/1942033

Vallet, P., Dhote, J. F., Le Moguedec, G., Ravart, M., & Pignard, G. (2006). Development of total aboveground volume equations for seven important forest tree species in France. Forest Ecology and Management, 229(1–3), 98–110. https://doi.org/10.1016/j.foreco.2006.03.013

Wang, H., Ouchi, K., Watanabe, M., Shimada, M., Tadono, T., Rosenqvist, A., Romshoo, S. A., Matsuoka, M., Moriyama, T., & Uratsuka, S. (2006). In search of the statistical properties of high-resolution polarimetric sar data for the measurements of forest biomass beyond the rcs saturation limits. IEEE Geoscience and Remote Sensing Letters, 3(4), 495–499. https://doi.org/10.1109/lgrs.2006.878299  

Wirasatriya, A., Pribadi, R., Iryanthony, S. B., Maslukah, L., Sugianto, D. N., Helmi, M., Ananta, R. R., Adi, N. S., Kepel, T. L., Ati, R. N. A., Kusumaningtyas, M. A., Suwa, R., Ray, R., Nakamura, T., & Nadaoka, K. (2022). Mangrove above-ground biomass and carbon stock in the karimunjawa-kemujan islands estimated from unmanned aerial vehicle-imagery. Sustainability, 14(2), 706. https://doi.org/10.3390/su14020706

Wu, H., & Xu, H. (2021). A review of the application of sampling techniques in forest biomass inventory. Journal of Southwest Forestry University (Natural Science), 41(3), 183–188. https://doi.org/10.11929/j.swfu.202007059

Wu, H., Xu, H., Tian, X., Zhang, W., & Lu, C. (2023). Multistage sampling and optimization for forest volume inventory based on spatial autocorrelation analysis. Forests, 14(2), 250. https://doi.org/10.3390/f14020250

Wu, Q., Yang, B., Pei, Z., & Wang, F. (2004). Influence of small features on crop area estimation at a using remote sensing and a doubler sampling method. Transactions of the Chinese Society of Agricultural Engineering, 20(3), 130–133.

Xiao, Y. (2004). Application of “3S” technology and sampling techniques to dynamic monitoring of forest resources. Journal of Southwest Forestry University, 24(2), 60–64.   

Xu, H., & Zhang, H. (2002). Study on forest biomass model. Yunnan Science and Technology Press.

Xu, X., & Cao, M. (2006). An analysis of the applications of remote sensing method to the forest biomass estimation. Geo-Information Science, 8(4), 122–128.

Xu, X., Du, H., Zhou, G., Ge, H., Shi, Y., Zhou, Y., Fan, W., & Fan, W. (2011). Estimation of aboveground carbon stock of Moso bamboo (Phyllostachys heterocycla var. pubescens) forest with a Landsat Thematic Mapper image. International Journal of Remote Sensing, 32(5), 1431–1448. https://doi.org/10.1080/01431160903551389

Yang, D. (1993). Selection of the best sampling control method for the total stock of forest resource secondary survey. Anhui Forestry Science and Technology, (3), 34–36.

Yu, P. (1974). Simple sampling survey method of firewood forest volume. Forest Resources Management, (3), 31–36.

Yu, Y., Pan, Y., Yang, X., & Fan, W. (2022). Spatial scale effect and correction of forest aboveground biomass estimation using remote sensing. Remote Sensing, 14(12), 2828. https://doi.org/10.3390/rs14122828

Zeng, P., Zhang, W., Li, Y., Shi, J., & Wang, Z. (2022). Forest total and component Above-Ground Biomass (AGB) estimation through C- and L-band polarimetric SAR Data. Forests, 13(3), 442. https://doi.org/10.3390/f13030442

Zeng, W., Chen, X., Pu, Y., Li, M., & Tang, S. (2018). Comparison of different methods for estimating forest biomass and carbon storage based on national forest inventory data. Forest Research, 31(1), 66–71. https://doi.org/10.13275/j.cnki.lykxyj.2018.01.008

Zeng, W., Luo, Q., & Peng, C. (1995). Study on efficiency of two-stage cluster sampling in forest inventory. Forest Research, (5), 483–488.

Zeng, W., & Tang, S. (2010). Using measurement error modeling method to establish compatible single-tree biomass equations system. Forest Research, 23(6), 797–803. https://doi.org/10.13275/j.cnki.lykxyj.2010.06.004

Zeng, W., & Tang, S. (2012). Modeling compatible single-tree aboveground biomass equations for masson pine (Pinus massoniana) in southern China. Journal of Forestry Research, 23(4), 593–598. https://doi.org/10.1007/s11676-012-0299-4

Zhang, X. (2000). Introduction to statistical learning theory and support vector machines. Acta Automatica Sinica, 26(1), 32–42. https://doi.org/10.16383/j.aas.2000.01.005

Zhao, J., Zhao, L., Chen, E., Li, Z., Xu, K., & Ding, X. (2022). An improved generalized hierarchical estimation framework with geostatistics for mapping forest parameters and its uncertainty: a case study of forest canopy height. Remote Sensing, 14(3), 568. https://doi.org/10.3390/rs14030568

Zhou, C., & Sun, Q. (2004). Causes and correction methods of sample structural deviation in hierarchical multi-stage unequal probability sampling. Journal of Statistics and Information, 19(6), 16–18. https://doi.org/10.3969/j.issn.1007-3116.2004.06.004

Zhu, Y., Feng, Z., Lu, J., & Liu, J. (2020). Estimation of forest biomass in Beijing (China) using multisource remote sensing and forest inventory data. Forests, 11(2), 163. https://doi.org/10.3390/f11020163