Optical remote sensing data have been widely used for estimating forest aboveground biomass (AGB). However, the use of optical images is often restricted by the saturation of spectral reflectance for forests that have multilayered and complex canopy structures and high AGB values and by the effect of spectral reflectance from underlayer shrub, grass, and bare soil for young stands. This usually leads to overestimations and underestimations for smaller and larger values, respectively, and makes it very challenging to improve the estimation accuracy of forest AGB. In this study, a novel methodology was proposed by incorporating stand age as a dummy variable into four models to improve the estimation accuracy of the Pinus densata forest AGB in Yunnan of Southwestern China. A total of eight models, including two parametric models (LM: linear regression model and LMC: LM with combined variables), two nonparametric models (RF: random forest and ANN: artificial neural network) without the age dummy variable, and four corresponding models with the age dummy variable (DLM, DLMC, DRF, and DANN), were compared to estimate AGB. Landsat 8 Operational Land Imager (OLI) images and 147 sample plots were acquired and utilized. The results showed that (1) compared with the two parametric models, the two nonparametric algorithms resulted in significantly greater estimation accuracies of Pinus densata forest AGB, and the increases of accuracy varied from 8% to 32% for 100 modeling plots and from 12% to 35% for 47 test plots based on root mean square error (RMSE); (2) compared with the models without the age dummy variable, the models with the age dummy variable greatly reduced the overestimations for the plots with AGB values smaller than 70 Mg/ha and the underestimations for the plots with AGB values larger than 180 Mg/ha and, thus, significantly improved the overall estimation accuracy by 14% to 42% for the modeling plots and by 32% to 44% for the test plots based on RMSE; and (3) the texture measures derived from the Landsat 8 OLI images contributed more to improving the estimation accuracy than the original spectral bands and other transformations. This implied that two nonparametric models, coupled with the use of the age dummy variable and texture measures, offered a great potential for improving the estimation accuracy of Pinus densata forest AGB.

Uncertainties in forest aboveground biomass (AGB) estimates resulting from over- and underestimations using remote sensing data have been widely studied. The uncertainties may occur due to the spatial effects of the plot data. In this study, we collected AGB data from a total of 147 Pinus densata forest sample plots in Yunnan of southwestern China and analyzed the spatial effects on the estimation of AGB. An ordinary least squares (OLS) and four spatial regression methods were compared for the estimation using Landsat 8-OLI images. Through the spatial analysis of AGB and residuals of model predictions, it was found that the spatial autocorrelation and heterogeneity of the plot data could not be ignored. Compared with the OLS, the impact of the spatial effects on AGB estimation could be reduced slightly by the spatial lag model (SLM) and the spatial error model (SEM) and greatly reduced by the linear mixed effects model (LMM) and geographically weighted regression (GWR) based on the distributions of prediction residuals, global Moran’s I, and Z score. The spatial regression models had better performance for model fitting and prediction because of the reduction in overestimations and underestimations for the forests with small and large AGB values, respectively. However, the reductions in the overestimations and underestimations varied depending on the spatial regression models. The GWR provided the most accurate predictions with the largest R² (0.665), the smallest root mean square error (34.507), and mean relative error (−9.070%) by greatly reducing the AGB interval for overestimations occurring and significantly increasing the threshold of AGB from 150 Mg/ha to 200 Mg/ha for underestimations. Thus, GWR offered the greatest potential of improving the estimation of Pinus densata forest AGB in Yunnan of southwestern China.

It is important to achieve estimates at the minimum cost, with no greater uncertainty than that which is appropriate for the objectives of the inventory. The aim of this study was to estimate the forest volume efficiently and accurately by sampling and analyzing the existing forest survey data, which is also a technical challenge. In this work, we used the spatial statistics tools in the ArcGIS software to analyze spatial autocorrelations with the data from the sixth to ninth continuous forest inventories (CFI) of Sichuan Province from 2002, 2007, 2012, and 2017. Based on the sampling framework of the CFI, we divided the sampling units into five groups using different methods to create the second-stage samples. Combined with the spatial autocorrelation analysis results, we selected certain samples from the collection of second-stage samples through stratified sampling to form the third-stage sampling units. We applied the sampling ratio, sampling accuracy, workload, and costs as the evaluation indexes for the sampling efficiency analysis. The main results are as follows: Before conversion, the forest volume density had a positively skewed distribution. There was substantial positive spatial autocorrelation, and its intensity was affected by the distance scale, especially at 187.3 km, where the spatial processes of clustering were most pronounced. At the significance level of α = 0.01, the high-volume stands were mainly concentrated in the Aba Prefecture, Garze Prefecture, and Liangshan Prefecture, while the low-volume stands were mainly concentrated in the Sichuan Basin region. The heterogeneous gatherings were staggered between the high-volume areas and low-volume areas, while the transition zone between the three prefecture regions and basin region was randomly distributed. With 95% reliability, the average estimation accuracy of the systematic sampling, random sampling, and cluster sampling in the second stage was 94.09%, which is less accurate than the CFI estimation accuracy. The mean correlation coefficients (R) between the estimated value of the forest volume and the observations of the systematic sampling, random sampling, and cluster sampling in the second stage were 0.95, 0.98, and 0.96, respectively. The relative differences (RD%) were −0.52, −0.39, and −0.36, respectively. The spatial stratified sampling in the third stage, which is based on spatial distribution pattern information, significantly reduced the sampling ratio to 1.68 per 10,000, compared with the average ratios of the CFI sampling and second-stage sampling, which were 13.73 per 10,000 and 2.75 per 10,000, respectively. With 95% reliability, the mean accuracy of the spatial stratified sampling in the third stage was 93.05%, the R was 0.94, and the RD% was −0.09. Spatial stratified sampling is more in line with the actual work conducted in annual surveys because it effectively reduces the sample size using prior spatial information, which can better meet the requirements of the annual output.

The aboveground biomass (AGB) of a forest is an important indicator of the forest’s terrestrial carbon storage and its relation to climate change. Due to the advantage of extensive spatial coverage and low cost, coarse-resolution remote sensing data is the main data source for wall-to-wall mapping of forest AGB at the regional scale. Despite this, improving the accuracy and efficiency of forest AGB estimation is a major challenge. In this study, two optical imageries, Moderate Resolution Imaging Spectroradiometer (MODIS) 500 m imagery and Fengyun-3C Visible and Infrared Radiometer (FY-3C VIRR) 1000 m imagery, were used and compared for forest AGB estimation in Yunnan Province, southwest China. One parametric approach, multiple linear regression (MLR), and two nonparametric approaches, k-nearest neighbor (KNN) and random forest (RF), were applied for the two imagery datasets, respectively. We evaluated the performance of the combination of remote sensing data and modeling approaches by comparing the accuracies and also explored the potential of FY-3C imagery data in forest AGB estimation at the regional scale as it was used for this purpose for the first time. We found that the machine learning models KNN and RF provided better results than MLR. From the three approaches for both MODIS and FY-3C imagery, RF performed best with R² values of 0.84 and 0.81 and RMSE of 23.18 and 23.43, respectively. Estimation of forest AGB based on MODIS was marginally better than the estimation based on FY-3C. FY-3C imagery could therefore be an additional optical remote sensing data source of coarse spatial resolution, comparable to MODIS data which has been widely used for regional forest AGB estimation. Indices related to forest canopy moisture levels from both types of imagery were sensitive to forest AGB. The RF model and MODIS imagery were then applied to map the spatial variation of forest AGB of Yunnan Province. As a result of our study, we determined that Yunnan Province has a total forest AGB of 2123.22 Mt, with a mean value of 58.05 t/ha for forestland in 2016.

The estimation of forest aboveground biomass (AGB) using Landsat 8 operational land imagery (OLI) images has been extensively studied, but forest aboveground biomass (AGB) is often difficult to estimate accurately, in part due to the multi-level structure of forests, the heterogeneity of stands, and the diversity of tree species. In this study, a habitat dataset describing the distribution environment of forests, Landsat 8 OLI image data of spectral reflectance information, as well as a combination of the two datasets were employed to estimate the AGB of the three common pine forests (Pinus yunnanensis forests, Pinus densata forests, and Pinus kesiya forests) in Yunnan Province using a parametric model, stepwise linear regression model (SLR), and a non-parametric model, such as random forest (RF) and support vector machine (SVM). Based on the results, the following conclusions can be drawn. (1) As compared with the parametric model (SLR), the non-parametric models (RF and SVM) have a better fitting performance for estimating the AGB of the three pine forests, especially in the AGB segment of 40 to 200 Mg/ha. The non-parametric model is more sensitive to the number of data samples. In the case of the Pinus densata forest with a sample size greater than 100, RF fitting provides better fitting performance than SVM fitting, and the SVM fitting model is better suited to the AGB estimation of the Pinus yunnanensis forest with a sample size of less than 100. (2) Landsat 8 OLI images exhibit superior accuracy in estimating the AGB of the three pine forests using a single dataset. Variables, such as texture and vegetation index variables, which can reflect the comprehensive reflection information of ground objects, play a significant role in estimating AGBs, especially the texture variables. (3) By incorporating the combined dataset with characteristics of tree species distribution and ground object reflectance spectrum, the accuracy and stability of AGB estimation of the three pine forests can be improved. Moreover, the employment of a combined dataset is also effective in reducing the number of estimation errors in cases with AGB less than 100 Mg/ha or exceeding 150 Mg/ha.

It is essential to analyze the spatial autocorrelation and heterogeneity of aboveground biomass (AGB). But it is difficult to accurately describe due to the lack of data in clear-cutting plots. Thus, measuring the AGB directly in a clear-cutting plot can provide a reference for accurately describing the spatial variation. Therefore, a 0.3-hectare clear-cutting sample plot of Pinus kesiya var. langbianensis natural forest was selected, and the AGB was calculated by each component. The intra-group variance was quantitatively described in terms of spatial heterogeneity, and the spatial autocorrelation was explored by global and local Moran’s I. The results indicated that (1) there was different spatial heterogeneity for the different trees and organs. The intra-group variance tended to be stable after 20 m for P. kesiya var. langbianensis (PK) and other upper trees (UPs) and after 10 m for the other lower trees (LTs). (2) The spatial autocorrelation of AGB and wood biomass was similar, while the bark biomass and foliage biomass were consistent. PK and other UPs also exhibited strong spatial autocorrelation, with maximum Moran’s I values of 0.1537 and 0.1644, respectively. (3) There was spatial heterogeneity in the different components except for the bark of PK. The lowest spatial heterogeneity was found for LT.