1. Introduction
With the intensification of global warming and the greenhouse effect, the carbon cycle has become a hot spot in global climate change research [
1,
2] and forests, as the largest carbon reservoir in terrestrial ecosystems [
3] store more than 80% of carbon [
4]. Since forest aboveground biomass (AGB) is a key biophysical parameter for measuring carbon and is generally used to quantify the contribution of forests to the global carbon cycle and the forest aboveground carbon stock (AGC) is an important parameter for assessing carbon sequestration capacity and carbon balance above the forest soil layer, it is crucial to understand the role of forests in carbon cycling and climate change.
Forest AGCs have traditionally been measured by field measurements, and although the traditional field measurement method is highly accurate, it is destructive to forests, laborious, and excessively expensive, mainly for small samples [
5], and is not suitable for spatially continuous AGC estimation [
6]. How to improve the accuracy of spatially continuous forest AGC estimation is still an active research area [
7].
Remote sensing can provide accurate and rapid information on vegetation cover over large areas and is now widely used to estimate forest AGB and AGCs [
8]. The vegetation indices generated from multispectral images can be used to estimate AGB and AGCs [
9]. However, the inability of optical images to penetrate the canopy and provide information on the vertical structure of the forest may lead to uncertain estimates in areas with dense canopies [
10,
11].
Light detection and ranging (LiDAR) can quickly and accurately acquire 3D information on vegetation [
12] and horizontal and vertical structural information on the forest canopy surface and can overcome the optical remote sensing saturation problem, and its emitted laser beam can penetrate the forest canopy. The acquired point cloud can be used to accurately estimate tree height, diameter at breast height, canopy size, and other structural attributes [
13,
14] and is therefore widely used to estimate AGCs. These LiDAR-derived structural attributes can be used with field measurements to estimate forest AGB through different models. Lu et al. [
15] used LiDAR data to extract individual tree structural parameters and input them into an anisotropic growth model to obtain sample-scale AGB. Chen et al. [
16] used a multiple regression model to estimate forest AGB and found that the accuracy of the nonlinear model was better than that of the linear model. Luo S et al. [
17,
18] analyzed the data fusion method of LiDAR and optical images from the perspective of the three-dimensional structure of forest stands, but it is only applicable to data estimation in small areas because of the small coverage area, which limits the utilization of remote sensing data. However, the variables used vary depending on the type of tree species in the study area; determining the best variables remains challenging, and due to the lack of forest canopy spectral information, the accuracy of tree species classification from LiDAR data is limited in complex vegetation conditions [
19].
To further improve the accuracy of forest AGB and AGC estimation, a data fusion method of optical images and LiDAR data was proposed [
7], which not only provides spectral information but also forest structure information and therefore can improve the accuracy of forest AGC estimation [
20,
21]. Popescu, S. et al. [
22,
23] showed that adding a large number of field measurements for modeling based on data fusion can also improve the estimation of forest AGB and AGC. Kim et al. [
24] indicated that combining spectral information with attributes derived from LiDAR data is more suitable for assessing AGB and AGC estimates than using optical images or LiDAR data alone. The above study was based on a large amount of field measurement data, and the variables were entered into the regression model to estimate the accuracy. However, there are challenges in finding a method to build a regression model to accurately estimate forest AGCs with relatively few field measurements. We hypothesized that forest-related attributes (tree height, diameter at breast height) combined with remote sensing spectral indices through anisotropic relationships could be used to accurately estimate forest AGCs at the sample plot scale with fewer field measurements.
In this study, we aim to develop a new method for forest AGC estimation by combining forest structural properties and spectral information using the anisotropic relationship assumed above. To improve the estimation of forest AGCs, in this study, we added horizontal structure variables from multispectral images to the measured data and vertical structure variables from LiDAR data to establish a multiple regression model and to explore the effect of combining multisource remote sensing variables on model accuracy. The accuracy was evaluated by comparing and analyzing the multivariate linear model and multivariate power model to estimate the forest AGCs at the sample plot scale. The spatial distribution of AGC density was mapped. This approach enables a more accurate estimation of forest AGCs and provides data support for the carbon cycle and sustainable forest management.
4. Discussion
4.1. Potential of LiDAR and Multispectral Image Synergy for Forest AGC Estimation
Compared to manual surveying, the backpack LiDAR allows for accurate scanning and real-time data integration while on the move, providing a more flexible and efficient way to collect data for forest inventory [
36]. The backpack LiDAR requires only one surveyor to carry the equipment across the survey site during data acquisition, significantly reducing time and costs and increasing efficiency [
37]. As shown in
Table 8, when collecting point cloud data of 10 m × 40 m samples, the traditional measurement method requires 3–4 people to collect the data at the same time, while the backpack LiDAR only requires 1 person to complete the collection; the traditional manual measurement takes approximately 36 min to finish measuring a sample area, while the backpack LiDAR only takes approximately 5 min; the traditional measurement method requires preprocessing of the collected data (inputting the data on the record sheet into Excel), which takes approximately 14 min, while the backpack LiDAR takes approximately 10 min to preprocess the data. The time required for preprocessing the collected data (inputting the data on the record sheet into Excel) is approximately 14 min, while the preprocessing time for the backpack LiDAR point cloud data is approximately 10 min, and the internal data processing time depends on the size of the data set and the computer configuration; overall, the time spent by the backpack LiDAR after collecting a sample plot is approximately 30 min faster than the traditional measurement method. The above description illustrates the time efficiency of backpack LiDAR. For optical data, acquiring airborne multispectral images in good weather conditions improves efficiency and reduces costs to some extent. Thus, the use of combined optical imagery and LiDAR further reduces the cost of assessing forest abatement. It makes it possible to map near real-time carbon stocks over large areas [
38].
Optical images have been applied in earlier studies to estimate forest AGB and AGC, but the results show that the penetration of optical signals is weak. Spectral images mainly record the horizontal structure of the forest and cannot record the vertical structure information of the woods. However, LiDAR can penetrate the forest canopy and record the vertical structure information. It is good to make up for the deficiency of optical images. In this study, there are two main reasons why the improvement with the addition of multispectral information is slight. The first reason could be that when visible light from multispectral data is saturated in dense forest areas [
39], the accuracy is lower in complex forest structures, resulting in deviations between the estimated and measured AGC of NDVI. Another reason could be that LiDAR forest structure properties are strongly correlated with AGC, and adding multispectral information does not improve much. Although these improvements are not significant, the novel multisensor earth observation approach combining satellite-based LiDAR data using machine learning techniques for optical data has enabled accurate measurements of carbon stocks and provided adequate data support for forest mitigation [
40,
41,
42,
43]. For example, Jiao et al. [
38] proposed a practical framework for assessing forest abatement using the fusion of optical satellite images and spaceborne LiDAR data. Shen et al. combined Landsat TM/ETM + and ALOS l-band SAR images of Guangdong Province to map AGB data of subtropical forests. The results showed a good correlation of AGB based on multisensor photos [
44]. Our results further suggest that combining LiDAR data and multispectral data is essential to improve the accuracy of AGB and AGC estimation.
4.2. Analysis of the Major Challenges and Uncertainties in Estimating Forest AGC
To address the challenges in vegetation biomass and carbon stock estimation (i.e., whether the simultaneous acquisition of large-scale forest structural and spectral information can improve the analysis of biomass and carbon stocks [
14,
45]), in this study, we combined forest structural attributes and spectral data to estimate forest AGC at the sample plot scale. (1) Although it is difficult to accurately capture changes in forest AGCs using only structural and spectral information for forest stands with complex structures, our proposed combined modeling approach with multisource data improved the accuracy of AGC estimation from 90.29% to 90.6% because the information on canopy spectral heterogeneity was provided by multispectral images. (2) In addition, we selected the best regression model to fit the AGC from multiple linear regression models and multiple power regression models. Multiple power regression models had higher AGC estimation accuracy than multiple linear regression models (
Table 5 and
Table 7). This suggests that the
L. gmelinii in our study area is consistent with a power anisotropy relationship. The use of the power anisotropic relationship can improve the accuracy of forest AGC estimation, and this relationship based on forest structural properties and spectral information is a new approach to improving forest AGC estimation. Based on forest structural attributes and spectral information, this method can be used to explore the relationship between tree metabolism and biomass [
46], and this relationship may be more stable in similar landscapes [
7].
There are still some uncertainties in this study. First, the backpack LiDAR laser beam cannot penetrate the lower canopy in dense forest structures [
47]. Second, the backpack LiDAR has challenges capturing the tops of trees in the upper canopy due to the shading caused by trees in the lower canopy, thus leading to significant differences between the height estimated by LiDAR and the measured height in this study. All the experimental results and conclusions of this study are currently valid only for coniferous forests with relatively simple stand structures, and their validity in broadleaf forests, mixed forests, or other forest types with more complex stand structures needs to be verified based on more forest sample plots and remote sensing data.
4.3. Estimating the Late Stage of Forest AGC Research and Outlook
This study focuses on the theme of combining multispectral images and LiDAR data for estimating regional-scale forest AGC, from field sample measurements to preprocessing such as atmospheric correction, radiometric correction, and geometric correction of multispectral images and LiDAR data cropping, resampling, denoising, filtering, near-ground classification, and normalized point clouds, to constructing a forest AGC estimation model for complex terrain conditions and then performing spatial extension of forest AGC at the regional scale. However, due to the shortage of measured data and the complexity of mountainous terrain conditions, the accuracy of regional forest AGC estimation by combining multisource remote sensing data is not currently accurate. A series of studies need to be continued.
At the current stage, calibration and validation still require high-quality field real-world data. Due to the complexity of mountainous terrain conditions, more accessible locations were selected to collect field inventories, leading to spatial discontinuity and discrete problems in producing regional forest AGC density spatial distribution maps. In future studies, we will try to select spatially continuous sample sites to collect data. Due to time constraints, limited samples were collected, and more minor sample data were only applicable to single wood or regional forest AGC estimation. It could not represent the forest AGC stock in the whole Dural Forest.
Poor GPS signals during data acquisition by backpack LiDAR can directly affect the quality of the track files, leading to point cloud solution failure or point cloud data solution failure with significant absolute coordinate errors. In addition, acquiring point cloud data with high-precision absolute coordinates is crucial for localizing individual trees in the sample area. Therefore, how to efficiently and accurately acquire absolute georeferenced point cloud data in a dense forest without GPS signals by backpack LiDAR is the focus of future research. In addition, the L. gmelinii-like ground has dense branches. To avoid scratching the backpack LiDAR instrument, the branches must be cut off in advance along the design route to ensure the safe operation of the backpack LiDAR. Therefore, the timing and quality of backpack LiDAR data acquisition due to forest stands and the accuracy of image data acquisition need to be further verified in more operating environments. Due to the small amount of spectral band data acquired from the multispectral data used in this study, the calculated vegetation spectral indices correlate less with the forest AGC. In future studies, the ability of UAVs with hyperspectral imagers at different flight altitudes to acquire regional forest vegetation spectra and combine them with LiDAR data for regional forest carbon stock inversion needs to be explored.
In general, the combination of LiDAR data and traditional remote sensing data can better complement each other’s data sources, which will help the acquisition and classification of feature information, improve the accuracy of estimating various parameters of ecosystems, and enhance the overall function of ecological monitoring and simulation. How to effectively combine multiple remote sensing data sources for environmental research is a hot issue.
5. Conclusions
In this study, based on the measured data, LiDAR vertical structure variables, and the addition of multispectral image horizontal structure variables to establish a multiple regression model, the following conclusions were obtained by comparing and analyzing the multivariate linear model and the multivariate power model to estimate the forest AGC at the sample scale.
A multivariate model was developed to predict the AGC and was tested for accuracy using LiDAR-estimated DBH and tree height as independent variables and the measured AGC as the dependent variable. The highest accuracy of the estimated AGC was found for the multiplicative power regression model based on LiDAR-estimated DBH (R2 = 0.903, RMSE = 10.91 Pg). The AGC values predicted by LiDAR ranged from 4.1–279.12 kg C. The accuracy of the LiDAR-estimated diameter at breast diameter was much higher than that of the tree height. This result may be due to the tall vegetation cover in the study area and the narrow beam of ground-based LiDAR, which makes it very difficult to search for targets in space due to the influence of occlusions and directly affects the interception probability and detection efficiency of the targets.
LiDAR data combined with the multispectral estimation of the AGC and determination of the accuracy showed that the multiplicative power regression model with the highest accuracy included the DBH-predicted AGC estimated by NDVI combined with LiDAR data (R2 = 0.906, RMSE = 10.87 Pg); the predicted AGC values were in the range of 3.93–449.07 kg C. The accuracy of AGC inversion using NDVI alone was extremely low (R2 = 0.056, RMSE = 36.25). This is because multispectral optical image data cannot accurately reflect vegetation height information, and spatial data effects are lacking.
Multiple regression analysis modeling demonstrated the potential of estimating the AGC from multisource remote sensing data. The model’s prediction accuracy was high (R
2 = 0.87–0.90) compared to the prediction accuracy in other studies [
48,
49]. The results showed that the addition of multispectral image variables to the predictive model for LiDAR estimation explained the variation in AGC estimation improvement. When LiDAR data and multispectral data are combined to estimate the AGC, LiDAR data are both accurate and include the spectral characteristics of multispectral optical images.
In general, the AGC is related not only to the structural features of trees that can be extracted from LiDAR data but also to the carbon coefficients that can be reflected in the multispectral information. Therefore, if both LiDAR and multispectral data are available, the fusion of LiDAR with multispectral data is the best method to accurately estimate forest AGC. This study could provide a valuable resource for researchers and forest managers to obtain more accurate AGC values.