# Estimating the Aboveground Carbon Density of Coniferous Forests by Combining Airborne LiDAR and Allometry Models at Plot Level

^{1}College of Natural Resources and Environment, Northwest A&F University, Yangling, China^{2}College of Forestry, Northwest A&F University, Yangling, China^{3}College of Landscape Architecture and Arts, Northwest A&F University, Yangling, China

Forest carbon density is an important indicator for evaluating forest carbon sink capacities. Accurate carbon density estimation is the basis for studying the response mechanisms of forest ecosystems to global climate change. Airborne light detection and ranging (LiDAR) technology can acquire the vertical structure parameters of forests with a higher precision and penetration ability than traditional optical remote sensing. Combining top of canopy height model (TCH) and allometry models, this paper constructed two prediction models of aboveground carbon density (ACD) with 94 square plots in northwestern China: one model is plot-averaged height-based power model and the other is plot-averaged daisy-chain model. The correlation coefficients (*R*^{2}) were 0.6725 and 0.6761, which are significantly higher than the correlation coefficients of the traditional percentile model (*R*^{2} = 0.5910). In addition, the correlation between TCH and ACD was significantly better than that between plot-averaged height (AvgH) and ACD, and Lorey’s height (LorH) had no significant correlation with ACD. We also found that plot-level basal area (BA) was a dominant factor in ACD prediction, with a correlation coefficient reaching 0.9182, but this subject requires field investigation. The two models proposed in this study provide a simple and easy approach for estimating ACD in coniferous forests, which can replace the traditional LiDAR percentile method completely.

## Introduction

Forest carbon storage accounts for 82.5% of terrestrial vegetation carbon storage, which is the main component of the vegetation carbon sink (Cusack et al., 2014; Kauranne et al., 2017). Accurate calculations of forest carbon stocks are a hot topic in the field of forest carbon sink research. At present, large-scale estimations of forest carbon sinks are mainly realized by means of traditional optical remote sensing. Generally, the relationship between field survey data and remote sensing extraction parameters is established first and then extrapolated to the whole research scope; this technique is essentially remote sensing-assisted sampling surveying (Drake et al., 2003; Chi et al., 2017; Koju et al., 2018).

Traditional optical remote sensing can extract the spectral information and horizontal structure information of vegetation. However, with increasing biomass, saturation occurs easily, which affects the estimation accuracy of forest carbon storage (Zhao et al., 2016). Light detection and ranging (LiDAR) detects the distance between a sensor and target by emitting laser pulses and receiving reflections from the ground object. Thus, LiDAR can acquire high-precision three-dimensional information of the object. Furthermore, LiDAR has a certain penetrating ability and can obtain vertical structure information of forests, improving the estimation accuracy of forest height and structure and forest carbon storage (Dubayah and Drake, 2000; Naesset and Bjerknes, 2001; Hudak et al., 2002; Gwenzi and Lefsky, 2014).

With the big-data progress and increasing storage space in recent years, airborne LiDAR has become an important means of forest resource surveys and carbon storage research (Guo et al., 2017; Swetnam et al., 2017). Data processing methods are mainly divided into plot-based inversion and individual tree-based inversion. However, due to the large number of trees, complex spatial structure of forests and canopy shielding effect, single-tree segmentation algorithms are not yet mature. Therefore, developing plot-based inversion methods is indispensable (Ayrey et al., 2017; Dechesne et al., 2017).

There are two main approaches for the estimation of carbon density based on plots. One approach is the use of a variety of machine learning algorithms to establish the relationship between measured carbon density and LiDAR percentile metrics, which can make full use of the information contained in the point cloud to obtain increasing precision (Zhao et al., 2011; McRoberts et al., 2016). However, the modeling process is a black box operation, and the prediction results are difficult to explain. The other approach is the establishment of LiDAR inversion models directly based on allometry models (Mascaro et al., 2011; Asner and Mascaro, 2014). The premise of this method is that there is a similar allometric growth law for plot-level biomass and single-tree biomass. The key to this approach is finding the appropriate allometric growth model and the corresponding LiDAR extraction parameters.

A multiple linear regression model based on LiDAR percentiles is a popular method for estimating forest carbon density or biomass, which is widely used and has acceptable precision in different forest area (e.g., Boudreau et al., 2008; Zhao et al., 2009; Ferraz et al., 2016; Jimenez-Berni et al., 2018). Among them, Naesset and Gobakken (2008) explained 88% variation in aboveground biomass (AGB) and 85% belowground biomass using LiDAR derived variables in boreal forest (1395 circular sample plots with size 200–400 m^{2}); Levick et al. (2016) used similar methods to obtain the fitting accuracy of 92% in 1 ha plots and 68% in 0.05 ha plots in temperate forest; Cao et al. (2016) established two regression models for estimating the AGB using multi-temporal LiDAR data of subtropical forest with *R*^{2} of 0.74 and 0.79 in 0.09 ha plots, respectively; Dubayah et al. (2010) estimated AGB in tropical forest with an *R*^{2} of 0.90 in 0.5 ha plots. Previous studies have demonstrated that the accuracy and form of percentiles models are closely related to the LiDAR instruments (Naesset, 2009; Silva et al., 2017) and plot size (Maltamo et al., 2011; Mascaro et al., 2011) except for intrinsic characteristics of forest.

In this study, we attempt to find a simple plot-based LiDAR extraction parameter, establish allometry models of the aboveground carbon density (ACD) of the northern coniferous forest, and evaluate the accuracy of these models. The objectives of this study are (1) the selection of the best parameter for ACD prediction from the following three plot-based LiDAR extraction parameters: top of canopy height (TCH), AVG (plot-averaged height), and Lorey’s height (LorH); (2) the proposal of direct and indirect fitting models of TCH and ACD and comparison of their accuracy and (3) calculation of the ACD of the study area with the proposed models and comparisons of the results and spatial distribution characteristics.

## Materials and Methods

### Study Area

The study was conducted 50 km southwest of Zhangye City, Gansu Province, Northwest China (Figure 1). The study area is approximately 264 ha, and its centre is at 100°15′E, 38°32′N. The elevation ranges from 2700 to 3200 m, the annual rainfall is 200 to 500 mm, and the monthly average temperature is 5.4 to 19.6°C. The main vegetation is natural pure Qinghai spruce (*Picea crassifolia*) forest, which has both naturally renewed young trees and tall over-ripe old trees. The ecoregion classification is “cascade pure conifer forest.”

### LiDAR Data and Processing

The LiDAR data used in this study were acquired on June 2008 using a LiteMapper 5600 instrument that recorded up to five returns per pulse, along with their intensity. The average flight altitude was 3560 m, the relative height over ground was 760 m, and the flight speed was 227 km/h. The laser scanner adopted RIEGL LMS-Q560, and the wavelength was 1550 nm. The laser pulse width was 3.5 ns, and the laser pulse divergence angle was less than or equal to 0.5 mrad. The LiDAR point cloud used the WGS84 coordinate system and the UTM projection zone 47 in the northern hemisphere. To increase the point density, the flight was repeated seven times over the study area with a side overlap of approximately 90%. As a result, the average point cloud interval was decreased to 0.54 m, and the average point cloud density was 3.43/m^{2}.

Subsequently, a set of metrics (Table 1) was derived from the point cloud using the LAStools software package^{1}. The main processing steps were as follows: (1) the point cloud was filtered and classified to ground, trees and noise; (2) the normalized point cloud (NPC, also referred to as height above ground) was calculated; (3) height percentiles, density percentiles and canopy cover (CC) were derived from the NPC corresponding to each plot; (4) the digital surface model (DSM) and digital elevation model (DEM)were interpolated from the first echo and the last echo of the point cloud, respectively. The canopy height model (CHM) was the difference of the first two. (5) The TCH was extracted from the CHM based on each plot (the mean value of 400 pixels per plot).

In addition, in order to explore the effect of CHM pixel size on ACD prediction, we generated 10 CHMs from NPC, with pixel sizes from 1 to 10 m. When the pixel location of CHM corresponds to a laser point, the point’s height value is used as the pixel value. If the location corresponds to multiple laser points, the average value of height is used as the pixel value. For the pixels without corresponding laser point, inverse distance weighted (IDW) is adopted for interpolation, which can ensure smooth transition between the target pixel and surrounding pixels (Montealegre et al., 2015).

### Field Data and Processing

To calibrate and validate the models, the plot data were acquired simultaneously with the LiDAR data. A total of 94 square plots (20 m × 20 m), which included 5734 trees, were used. The four corners and the centre of each plot were measured using differential GPS (DGPS), and the error was less than 10 cm. For each tree with a diameter at breast height (DBH) greater than 5 cm, the tree type, diameter, height to crown base, crown width in cardinal directions, crown class, and crown transparency were measured. DBH was measured on all trees using a diameter tape, and the heights of all trees were measured using a laser ranging hypsometer with theoretical accuracy up to the decimeter level. Considering the canopy occlusion and human error, the average accuracy of the measured tree height was better than 0.5 m.

Using the species-specific allometry Eqs 1–4 in the study area (Wang et al., 1998; He et al., 2013), tree biomass components (stem, branch, foliage and fruit) were calculated from DBH and height. The AGB of each tree was equal to the sum of the AGB components and was then summed to obtain the AGB of each plot. These equations were constructed by destructive sampling and the fitting precision reached 0.9887, 0.9568, 0.8662, and 0.9340, respectively. Because the study area is a nature reserve and pure spruce forest with little human disturbance, it is believed that the field-estimated AGB of this study can also achieve such accuracy. Finally, the AGB was converted to the ACD using the conversion coefficient of 0.5034, which was attained using the potassium dichromate oxidation method on samples by Wang et al. (2000) in the same area. Additionally, Lorey’s height (LorH; Table 1) was also calculated based on each plot using Eq. 5, and the LorH values were compared with the TCH values extracted from LiDAR. For the same purpose, the average height (AvgH) of the plots was also calculated.

where *D* is DBH, *H* is tree height, BA* _{i}* and

*H*

_{i}are the basal area and the height of the

*i*th tree, respectively, and

*a*and

*b*are regression fitting coefficients.

### Plot Allometry Models

The use of an allometry model is the main means of forest biomass calculation. This type of model is obtained by the regression of the sample forest harvesting and tree-measuring metrics and is a single-tree model for specific tree species in a specific region. The present work imitates the form of single-tree models at the plot level to find a suitable plot-level LiDAR metric to replace traditional tree-measuring metrics.

An idealized and simple tree allometry equation for special species is:

Since DBH is the most easily accessible and accurately measurable tree indicator, and there is an intrinsic relationship between the DBH and tree height, the model is widely used (Chave et al., 2005).

However, Eq. 6 cannot explain the variability of diameter and tree height growth caused by tree age, forest density, site conditions and management measures; the introduction of the tree height factor is necessary.

where *H* represents the tree height (m), and *a*, *b*1, and *b*2 are fitted coefficients.

The essence of LiDAR is ranging, which can directly estimate tree height. Therefore, this paper applies Eqs 6 and 7 to Eqs 8 and 9, which are plot-averaged height-based allometry models, and fits the equations as follows:

where ACD represents the aboveground carbon density (Mg C ha^{–1}), and *H* represents AvgH, LorH, or BA computed from field-measured data and TCH extracted from the CHM. BA’ is the fitted result in Eq. 10.

### Percentile Model

Light detection and ranging is highly sensitive to the three-dimensional structure of forest, because laser pulses can penetrate the canopy and then record all echo signals from the ground to the canopy surface. Therefore, a series of LiDAR metrics, such as height percentile, density percentile, variation coefficient, etc., have been successively extracted to capture key information of forest canopy (Nilsson, 1996; Lefsky et al., 2002; Naesset, 2002; Nelson et al., 2004).

In this study, the height and density percentiles extracted from the NPC were used to regression fit the ACD calculated from the field investigation data. The model is as follows, and the independent variables are described in Table 1. The prediction results of this model (Eq. 11) were compared with the prediction results of the allometry models (Eqs 8, 9) proposed in this paper.

### Model Fitting and Evaluation

All models in this study were fitted by the least squares (OLS) method (Meng et al., 2018). This method is simple and reliable, avoiding the algorithm differences of different fitting methods and making the fitting results more contrasting. All power models were first converted to a linear model by natural log transformation for regression fitting to correct both non-normality and heteroscedasticity. Then, the correlation of coefficients (*R*^{2}) and back-log root-mean-squared errors (RMSE) were employed to compare the performance of the models, and 10-fold cross-validation analysis was used to evaluate the stability of the models.

### Model Application

The study area was divided into a 20 m × 20 m grid using GIS software (Figure 2). The size and direction of the grid were the same as those of the field plots in order to reduce possible errors. Then, LiDAR parameters were extracted using each grid from the CHM and NPC and introduced into the allometry models (Eqs 8, 11) proposed in this paper for calculation. Thereby, ACD distribution maps of the study area were obtained, compared and evaluated. Since there was no point cloud at the boundary of the study area, grid incompleteness due to cropping did not affect the final prediction.

## Results

### Field Investigation Data Analysis

Using the comparison between the simple power-law model (Eq. 8) of ACD and the three plot-averaged metrics (AvgH, LorH, and BA)calculated from the field inventory, we found that the BA explains 91.8% of the variation in ACD, which is much higher than the 39.5% explained by AvgH and 10.1% explained by LorH (Table 2), and the convergence of BA is much better than that of AvgH and LorH (Figure 3), which indicates that BA is the optimal plot-averaged indicator for the inversion of ACD. It is not surprising that BA is a stronger predictor of AGB than height is, because BA can be measured with a lot better accuracy than height, and because DBH is weighted higher than tree height in Eq. 7. We also found that the exponent b is close to 1 (*b* = 1.0300), which indicates that the ACD is nearly linearly related to BA, and the cross-validated *R*^{2} value (0.9182) is reduced by 0.004 compared with the model-fitted *R*^{2} (0.9143), which indicates that the model using BA tends to be applicable and stable. This result means that, when we want to obtain the ACD of plots, we can discard the exhaustive field inventory data and only need to perform spatially explicit point-based measurements using the relascope or prism method.

**Figure 3.** The linear relationship between ln(AvgH) and ln(ACD) **(A)**, ln(LorH) and ln(ACD) **(B)**, ln(BA)and ln(ACD) **(C)**. The shadowed region shows the 95% confidence interval. **(D)** Boxplots of the field-surveyed BA and predicted BA using plot-averaged height (AvgH) and Lorey’s height (LorH).

Moreover, although AvgH and LorH are the most commonly used plot-averaged height indicators, when they were applied in Eq. 8 to predict ACD, the effect was poor, with *R*^{2} values of 0.3954 and 0.1014, respectively, and RMSE values of 13.5561 (Mg C ha^{–1}) and 15.9938 (Mg C ha^{–1}), respectively; the results with AvgH are slightly better than those with LorH (Table 2). This suggests that the plot-averaged height alone does not account for the variation in the ACD. Therefore, the plot-level ACD estimation (EACD) based on LiDAR should exclude the AvgH and LorH steps and directly fit ACD with LiDAR-extracted metrics (Figure 3D).

### TCH Models and the Comparisons

The plot-level LiDAR metric (TCH) was taken into the ideal simplest allometry (Eq. 8) and was subjected to log changes and linear fits. The result showed that TCH could explain 67.25% of the variation in ACD (Figure 4A), suggesting that TCH is a simple and effective predictor of ACD from LiDAR and that the classic allometry model (Eq. 6) can be extended from the tree-level to the plot-level scale.

**Figure 4.** **(A)** Linear relationship between ln(TCH) and ln(ACD). **(B)** Linear relationship between TCH and the basal area (BA). The shadowed region shows the 95% confidence interval. **(C)** ACD estimated using top of canopy height (TCH) from LiDAR compared to field-surveyed ACD. **(D)** ACD estimated using TCH and BA predicted from TCH compared to field-surveyed ACD. The black line is a 1:1 reference line.

Using regression by ordinary least squares, we modeled variation in BA to TCH for 94 plots, with resulting values of *R*^{2} = 0.6066 and RMSE = 5.1749 m^{2} ha^{–1} (Eq. 10, Figure 4B, and Table 3). By substituting this regression result into Eq. 9, the EACD could be generated without field inventory data. However, in a comparison of the fitting results, we found that the scatter plots were almost the same (Figures 4C,D); *R*^{2} only increased by 0.0036, and RMSE increased by 0.1163 Mg C ha^{–1} (Table 3). This strongly suggests that the daisy-chain method of TCH cannot achieve the same ACD prediction as the field-measured BA. Therefore, if we only use LiDAR-extracted TCH, the height-diameter model (Eq. 9) and height model (Eq. 8) have no essential difference in accuracy. Thus, this paper ultimately chose Eq. 8 as the final allometry model, and the parameter H used TCH.

### Percentile Model and a Comparison of Results

The following result (Eq. 12) was obtained by the multiple regression fitting of the surveyed ACD of 94 plots and the LiDAR percentile metrics listed in Table 1.

The ACD in the study area is closely related to h25 and d95. These two parameters can explain 59.1% of the ACD variation, with a RMSE of 11.6304 Mg C ha^{–1} (Figure 5A), and the prediction accuracy is lower than the 67.25% of Eq. 8 using TCH. It demonstrated that the height allometry model proposed in this paper can replace the traditional LiDAR percentile model with improved precision.

**Figure 5.** **(A)** Estimated ACD values of the percentile model versus the field investigation. **(B)** Boxplots of field-surveyed ACD and ACD estimated using the percentile model and TCH model.

Figure 5B illustrates the difference between the predicted values and the survey values of ACD. The median values of the MLR model and the TCH model are near 52 Mg C ha^{–1}, which is slightly lower than the measured ACD. Furthermore, the range of predicted values of the TCH model is slightly smaller than the surveyed value range, which is larger than the range of the MLR model. Therefore, compared with the MLR model, the TCH model has a wider prediction range and can represent larger and smaller values of ACD.

### Model Application

All grid values in the study area were calculated using our proposed TCH allometry model and percentile model, and then maps of ACD were produced. Figure 6 shows that the spatial distribution of the two maps is very similar. The high-density area of the map from the percentile model is slightly larger than that of the map from the TCH model (blue circle), and the low-density area demonstrates the opposite trend (blue rectangle). In addition, the density distribution percentages of the two maps are basically the same as those shown in the two pie charts. According to the grid statistics, the average ACD from the TCH model is 41.49 Mg C ha^{–1}, and the maximum value is 104.70 Mg C ha^{–1}, which is slightly larger than the values of 40.13 and 95.46 Mg C ha^{–1} from the percentile model. This resulted in an overall aboveground carbon reserve of the study area of 5535.54 Mg for the TCH model and 5433.06 Mg for the percentile model; the difference between the two models is only 1.89%. Although the accuracies of the TCH and the MLR models are not much different, the TCH model is much simpler and easier than the MLR model.

**Figure 6.** Study area carbon density map predicted using the TCH model **(A)** and percentile model **(B)**. Square is low density area, Round is high density area.

## Discussion

Our original purpose was to find a suitable plot-averaged LiDAR parameter and use existing allometry models to quickly and accurately predict the forest carbon density. The exponential model of TCH captures 67.25% of ACD changes (Table 3) and has a higher accuracy than the traditional percentile model of this study. We also realized that the accuracy of our prediction is relatively low. The possible reasons are (1) the plot size of 20 m × 20 m is relatively small, and the edge effect is obvious; (2) point cloud density is not enough, and the conical crown of spruce is not captured accurately and (3) the penetration rate of point cloud is insufficient, and the detection of lower wood is limited. However, for the research objective, we did effectively improve the accuracy of LiDAR’s prediction of ACD, simplify the prediction steps and solidify the form of the prediction model.

Moreover, since the TCH is derived from the mean of the CHM based on the plots, the TCH is also subject to the pixel size. We extracted 10 CHMs from the LiDAR point cloud, with pixel sizes from 1 to 10 m, and then extracted the corresponding TCHs to fit the ACD. As the pixel size increased, *R*^{2} continually decreased, and the RMSE continually increased (Figure 7). This result indicated that the smaller the CHM pixel is, the better the fitting effect of TCH will be. This study was limited to a point cloud density of 3.43/m^{2}, so the minimum pixel size was 1 m. In addition, we found that when the pixel sizes were 5 and 7 m, the fitting effect fluctuated slightly, but this fluctuation did not affect the overall law. The reason for this finding requires further study. Similarly, fitting accuracy is also limited by the size of the plot and the number of samples. A larger plot size means a smaller boundary effect, and a larger plot number means a smaller outlier influence (Ni et al., 2014; Gwenzi and Lefsky, 2017). However, larger plots are more expensive and time-consuming than smaller plots, so finding an optimal plot size and number in coniferous forest will be an important task for future study.

**Figure 7.** Fitting trends of TCH and ACD under different CHM precisions. **(A)** Declining *R*^{2}; **(B)** Increasing RMSE.

The three-dimensional visualization of the point cloud in plot no. 1 (Figure 8A) suggests that the forest point cloud includes the crown, some of the lower layer, some of the trunks and the woodland gap. Therefore, the TCH data derived from the point cloud also contain the above information. However, the average height of the plot (AvgH) ignores the forest gap and is therefore slightly higher than the TCH (Figure 8B). Although LorH is widely used for the estimation of forest biomass (Mitchard et al., 2012; Gwenzi and Lefsky, 2016), LorH is mainly used to evaluate site quality and mostly reflects the largest trees in the forest; therefore, its value is larger than that of TCH and AvgH (Figure 8B), and LorH is not applicable for fitting the ACD of irregular and mixed forests. This explains why TCH is the optimal ACD predictor.

**Figure 8.** **(A)** 3D scatter plot of the normalized point cloud (NPC) in plot 1. **(B)** Boxplots of TCH, AvgH and LorH. White bullet is mild outliers of boxplot.

We also recognized a flaw in the ACD prediction at the plot scale. Whether in the field measurement phase, the plot-based TCH extraction phase, or the final ACD prediction phase, our resolution is fixed at 20 m × 20 m. This inevitably leads to the conversion of the continuous ACD distribution in nature into a discontinuous distribution, which may cause a large jump phenomenon at the boundary. Therefore, selecting the appropriate interpolation algorithm to restore the continuity of the ACD will help improve the prediction accuracy of our proposed models (Loquin and Strauss, 2010; DeWitt et al., 2017). In addition, we only adopted linear regression fitting based on the least squares method in this paper, and although this method is simple and practical, it is not necessarily the best method. With the rise of machine learning in LiDAR research (Zhou et al., 2017; Jin et al., 2018; Lin et al., 2018), it will be necessary to compare various machine learning algorithms in future research to find the best way to fit the allometry models.

Finally, we must emphasize that although our proposed TCH-based allometric approach is an efficient LiDAR-assisted ACD prediction method, the allometry model used for plot calculation is generally targeted to a specific region and species (Picard et al., 2015; Duque et al., 2017), so it is necessary to re-select an appropriate allometry model for other tree species and ecological regions when our method is used. Moreover, developing a general ACD prediction model based on LiDAR for forests across ecological regions and species will be the focus of our future research.

## Conclusion

Using the traditional allometry growth model theory, this paper proposed two models based on TCH extracted from LiDAR data. The first model was a simple power model (only using TCH) based on the diameter allometry, and the second model was a daisy-chain model (TCH → BA′ → ACD) based on diameter-height allometry. A comparison of the results suggested there was little difference in the fitting accuracy and error distribution between models. In addition, this paper compared the traditional LiDAR percentile method with the proposed method and found that the latter method had a higher precision, fewer parameters, more concise steps and more stable forms than the former method. Furthermore, the implicit hypothesis in our study, the traditional allometry model of individual trees can be extrapolated to the plot scale, was confirmed. The LiDAR-assisted ACD estimation method proposed in this study will accelerate the application of airborne LiDAR technology in forest carbon density measurements and provide an accurate data basis for forest ecosystem research.

## Data Availability

Publicly available datasets were analyzed in this study. This data can be found here: http://www.heihedata.org/data/.

## Author Contributions

HH and WL: conceptualization, methodology, and writing-original draft preparation. HH: software, resources, and data curation. XZ and PZ: validation. HH and XZ: formal analysis. HH, WL, and XZ: investigation. QC: writing, review, and editing. HH and PZ: visualization. QC and PZ: supervision, project administration, and funding acquisition.

## Funding

This research was funded by the National Key Research and Development Program of China (Grant No. 2016YFD0600203), the Introducing the International Advanced Agricultural Science and Technology Program (948 Program) of China (Grant No. 2014-4-27), the National Natural Science Foundation of China (Grant No. 41701398), and Chinese Universities Scientific Fund (Grant No. 2452017108).

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

The authors would like to thank Northwest Institute of Eco-Environment and Resources (CAS) for support in LiDAR data acquisition and also thank Gansu Province Qilian Water Resource Conservation Forest Research Institute for support in field investigation.

## Footnotes

## References

Asner, G. P., and Mascaro, J. (2014). Mapping tropical forest carbon: calibrating plot estimates to a simple LiDAR metric. *Remote Sens. Environ.* 140, 614–624. doi: 10.1016/j.rse.2013.09.023

Ayrey, E., Fraver, S., Kershaw, J. A., Kenefic, L. S., Hayes, D., Weiskittel, A. R., et al. (2017). Layer stacking: a novel algorithm for individual forest tree segmentation from LiDAR point clouds. *Can. J. Remote Sens.* 43, 16–27. doi: 10.1080/07038992.2017.1252907

Boudreau, J., Nelson, R. F., Margolis, H. A., Beaudoin, A., Guindon, L., and Kimes, D. S. (2008). Regional aboveground forest biomass using airborne and spaceborne LiDAR in Quebec. *Remote Sens. Environ.* 112, 3876–3890. doi: 10.1016/j.rse.2008.06.003

Cao, L., Coops, N. C., Innes, J. L., Sheppard, S. R. J., Fu, L. Y., Ruan, H. H., et al. (2016). Estimation of forest biomass dynamics in subtropical forests using multi-temporal airborne LiDAR data. *Remote Sens. Environ.* 178, 158–171. doi: 10.1016/j.rse.2016.03.012

Chave, J., Andalo, C., Brown, S., Cairns, M. A., Chambers, J. Q., Eamus, D., et al. (2005). Tree allometry and improved estimation of carbon stocks and balance in tropical forests. *Oecologia* 145, 87–99. doi: 10.1007/s00442-005-0100-x

Chi, H., Sun, G., Huang, J., Li, R., Ren, X., Ni, W., et al. (2017). Estimation of forest aboveground biomass in changbai mountain region using ICESat/GLAS and Landsat/TM Data. *Remote Sens.* 9:707. doi: 10.3390/rs9070707

Cusack, D. F., Axsen, J., Shwom, R., Hartzell-Nichols, L., White, S., and Mackey, K. R. (2014). An interdisciplinary assessment of climate engineering strategies. *Front. Ecol. Environ.* 12, 280–287. doi: 10.1890/130030

Dechesne, C., Mallet, C., Le Bris, A., and Gouet-Brunet, V. (2017). Semantic segmentation of forest stands of pure species combining airborne lidar data and very high resolution multispectral imagery. *ISPRS J. Photogramm. Remote Sens.* 126, 129–145. doi: 10.1016/j.isprsjprs.2017.02.011

DeWitt, J. D., Warner, T. A., Chirico, P. G., and Bergstresser, S. E. (2017). Creating high-resolution bare-earth digital elevation models (DEMs) from stereo imagery in an area of densely vegetated deciduous forest using combinations of procedures designed for lidar point cloud filtering. *Gisci. Remote Sens.* 54, 552–572. doi: 10.1080/15481603.2017.1295514

Drake, J. B., Knox, R. G., Dubayah, R. O., Clark, D. B., Condit, R., Blair, J. B., et al. (2003). Above-ground biomass estimation in closed canopy Neotropical forests using lidar remote sensing: factors affecting the generality of relationships. *Glob. Ecol. Biogeogr.* 12, 147–159. doi: 10.1046/j.1466-822X.2003.00010.x

Dubayah, R. O., Sheldon, S. L., Clark, D. B., Hofton, M. A., Blair, J. B., Hurtt, G. C., et al. (2010). Estimation of tropical forest height and biomass dynamics using lidar remote sensing at La Selva, Costa Rica. *J. Geophys. Res.* 115:G00E09. doi: 10.1029/2009jg000933

Duque, A., Saldarriaga, J., Meyer, V., and Saatchi, S. (2017). Structure and allometry in tropical forests of Choco. Colombia. *For. Ecol. Manage.* 405, 309–318. doi: 10.1016/j.foreco.2017.09.048

Ferraz, A., Saatchi, S., Mallet, C., Jacquemoud, S., Goncalves, G., Silva, C. A., et al. (2016). Airborne lidar estimation of aboveground forest biomass in the absence of field inventory. *Remote Sens.* 8:653. doi: 10.3390/rs8080653

Guo, X., Coops, N. C., Tompalski, P., Nielsen, S. E., Bater, C. W., and Stadt, J. J. (2017). Regional mapping of vegetation structure for biodiversity monitoring using airborne lidar data. *Ecol. Inform.* 38, 50–61. doi: 10.1016/j.ecoinf.2017.01.0051574-9541

Gwenzi, D., and Lefsky, M. A. (2014). Modeling canopy height in a savanna ecosystem using spacebome lidar waveforms. *Remote Sens. Environ.* 154, 338–344. doi: 10.1016/j.rse.2013.11.024

Gwenzi, D., and Lefsky, M. A. (2016). Plot-level aboveground woody biomass modeling using canopy height and auxiliary remote sensing data in a heterogeneous savanna. *J. Appl. Remote Sens.* 10:016001. doi: 10.1117/1.jrs.10.016001

Gwenzi, D., and Lefsky, M. A. (2017). Spatial modeling of lidar-derived woody biomass estimates collected along transects in a heterogeneous savanna landscape. *IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens.* 10, 372–384. doi: 10.1109/Jstars.2016.2582148

He, Q. S., Chen, E. X., An, R., and Li, Y. (2013). Above-ground biomass and biomass components estimation using LiDAR data in a coniferous forest. *Forests* 4, 984–1002. doi: 10.3390/f4040984

Hudak, A. T., Lefsky, M. A., Cohen, W. B., and Berterretche, M. (2002). Integration of lidar and Landsat ETM plus data for estimating and mapping forest canopy height. *Remote Sens. Environ.* 82, 397–416. doi: 10.1016/S0034-4257(02)00056-1

Jimenez-Berni, J. A., Deery, D. M., Rozas-Larraondo, P., Condon, A. T. G., Rebetzke, G. J., James, R. A., et al. (2018). High throughput determination of plant height, ground cover, and above-ground biomass in wheat with LiDAR. *Front. Plant Sci.* 9:237. doi: 10.3389/fpls.2018.00237

Jin, S., Su, Y., Gao, S., Wu, F., Hu, T., Liu, J., et al. (2018). Deep learning: individual maize segmentation from terrestrial lidar data using faster R-CNN and regional growth algorithms. *Front. Plant Sci.* 9:866. doi: 10.3389/fpls.2018.00866

Kauranne, T., Joshi, A., Gautam, B., Manandhar, U., Nepal, S., Peuhkurinen, J., et al. (2017). LiDAR-assisted multi-source program (LAMP) for measuring above ground biomass and forest carbon. *Remote Sens.* 9:154. doi: 10.3390/Rs9020154

Koju, U. A., Zhang, J., Maharjan, S., Zhang, S., Bai, Y., Vijayakumar, D. B. I. P., et al. (2018). A two-scale approach for estimating forest aboveground biomass with optical remote sensing images in a subtropical forest of Nepal. *J. For. Res.* 31, 1–18.

Lefsky, M. A., Cohen, W. B., Parker, G. G., and Harding, D. J. (2002). Lidar remote sensing for ecosystem studies. *Bioscience* 52, 19–30.

Levick, S. R., Hessenmoller, D., and Schulze, E. D. (2016). Scaling wood volume estimates from inventory plots to landscapes with airborne LiDAR in temperate deciduous forest. *Carbon Balance Manag.* 11:7. doi: 10.1186/s13021-016-0048-7

Lin, Y., Jiang, M., Pellikka, P., and Heiskanen, J. (2018). Recruiting conventional tree architecture models into state-of-the-art LiDAR mapping for investigating tree growth habits in structure. *Front. Plant Sci.* 9:220. doi: 10.3389/fpls.2018.00220

Loquin, K., and Strauss, O. (2010). “Convolution filtering and mathematicalmorphology on an image: a unified view,” in *Linear Filtering and Mathematical Morphology on an Image: A Bridge. ICIP: International Conference on Image Processing*, Cairo, 3965–3968.

Maltamo, M., Bollandsas, O. M., Naesset, E., Gobakken, T., and Packalen, P. (2011). Different plot selection strategies for field training data in ALS-assisted forest inventory. *Forestry* 84, 23–31. doi: 10.1093/forestry/cpq039

Mascaro, J., Detto, M., Asner, G. P., and Muller-Landau, H. C. (2011). Evaluating uncertainty in mapping forest carbon with airborne LiDAR. *Remote Sens. Environ.* 115, 3770–3774. doi: 10.1016/j.rse.2011.07.019

McRoberts, R. E., Domke, G. M., Chen, Q., Naesset, E., and Gobakken, T. (2016). Using genetic algorithms to optimize k-nearest neighbors configurations for use with airborne laser scanning data. *Remote Sens. Environ.* 184, 387–395. doi: 10.1016/j.rse.2016.07.007

Meng, R., Dennison, P. E., Zhao, F., Shendryk, I., Rickert, A., Hanavan, R. P., et al. (2018). Mapping canopy defoliation by herbivorous insects at the individual tree level using bi-temporal airborne imaging spectroscopy and LiDAR measurements. *Remote Sens. Environ.* 215, 170–183. doi: 10.1016/j.rse.2018.06.008

Mitchard, E. T. A., Saatchi, S. S., White, L. J. T., Abernethy, K. A., Jeffery, K. J., Lewis, S. L., et al. (2012). Mapping tropical forest biomass with radar and spaceborne LiDAR in Lope National Park, Gabon: overcoming problems of high biomass and persistent cloud. *Biogeosciences* 9, 179–191. doi: 10.5194/bg-9-179-2012

Montealegre, A., Lamelas, M., and Riva, J. (2015). Interpolation routines assessment in ALS-derived digital elevation models for forestry applications. *Remote Sens.* 7, 8631–8654. doi: 10.3390/rs70708631

Naesset, E. (2002). Predicting forest stand characteristics with airborne scanning laser using a practical two-stage procedure and field data. *Remote Sens. Environ.* 80, 88–99. doi: 10.1016/s0034-4257(01)00290-5

Naesset, E. (2009). Effects of different sensors, flying altitudes, and pulse repetition frequencies on forest canopy metrics and biophysical stand properties derived from small-footprint airborne laser data. *Remote Sens. Environ.* 113, 148–159. doi: 10.1016/j.rse.2008.09.001

Naesset, E., and Bjerknes, K. O. (2001). Estimating tree heights and number of stems in young forest stands using airborne laser scanner data. *Remote Sens. Environ.* 78, 328–340. doi: 10.1016/S0034-4257(01)00228-0

Naesset, E., and Gobakken, T. (2008). Estimation of above- and below-ground biomass across regions of the boreal forest zone using airborne laser. *Remote Sens. Environ.* 112, 3079–3090. doi: 10.1016/j.rse.2008.03.004

Nelson, R., Short, A., and Valenti, M. (2004). Measuring biomass and carbon in Delaware using an airborne profiling LIDAR. *Scand. J. For. Res.* 19, 500–511. doi: 10.1080/02827580410019508

Ni, X. L., Park, T., Choi, S. H., Shi, Y. L., Cao, C. X., Wang, X. J., et al. (2014). Allometric scaling and resource limitations model of tree heights: Part 3. Model optimization and testing over continental china. *Remote Sens.* 6, 3533–3553. doi: 10.3390/rs6053533

Nilsson, M. (1996). Estimation of tree heights and stand volume using an airborne lidar system. *Remote Sens. Environ.* 56, 1–7. doi: 10.1016/0034-4257(95)00224-3

Picard, N., Rutishauser, E., Ploton, P., Ngomanda, A., and Henry, M. (2015). Should tree biomass allometry be restricted to power models? *For. Ecol. Manage.* 353, 156–163. doi: 10.1016/j.foreco.2015.05.035

Silva, C. A., Hudak, A. T., Vierling, L. A., Klauberg, C., Garcia, M., Ferraz, A., et al. (2017). Impacts of airborne lidar pulse density on estimating biomass stocks and changes in a selectively logged tropical forest. *Remote Sens.* 9:1068. doi: 10.3390/Rs9101068

Swetnam, T. L., Gillan, J. K., Sankey, T. T., Mcclaran, M. P., Nichols, M. H., Heilman, P., et al. (2017). Considerations for achieving cross-platform point cloud data fusion across different dryland ecosystem structural states. *Front. Plant Sci.* 8:2144. doi: 10.3389/fpls.2017.02144

Wang, J., Che, K., Fu, H., Chang, X., Song, C., and Ke, H. (1998). Study on biomass of water conservation forest in qilian mountains. *J. Fujian Coll. For.* 18:5. doi: 10.13324/j.cnki.jfcf.1998.04.009

Wang, J., Che, K., and Jang, Z. (2000). Study on carbon balance of qinghai spruce forest in qilian mountains. *J. Northwest For. Univ.* 15, 9–14. doi: 10.3969/j.issn.1001-7461.2000.01.002

Zhao, K. G., Popescu, S., Meng, X. L., Pang, Y., and Agca, M. (2011). Characterizing forest canopy structure with lidar composite metrics and machine learning. *Remote Sens. Environ.* 115, 1978–1996. doi: 10.1016/j.rse.2011.04.001

Zhao, K. G., Popescu, S., and Nelson, R. (2009). Lidar remote sensing of forest biomass: a scale-invariant estimation approach using airborne lasers. *Remote Sens. Environ.* 113, 182–196. doi: 10.1016/j.rse.2008.09.009

Zhao, P., Lu, D., Wang, G., Wu, C., Huang, Y., and Yu, S. (2016). Examining spectral reflectance saturation in landsat imagery and corresponding solutions to improve forest aboveground biomass estimation. *Remote Sens.* 8:469. doi: 10.3390/rs8060469

Keywords: LiDAR, AGB, ACD, allometry model, NPC, CHM, TCH

Citation: Hao H, Li W, Zhao X, Chang Q and Zhao P (2019) Estimating the Aboveground Carbon Density of Coniferous Forests by Combining Airborne LiDAR and Allometry Models at Plot Level. *Front. Plant Sci.* 10:917. doi: 10.3389/fpls.2019.00917

Received: 07 March 2019; Accepted: 28 June 2019;

Published: 10 July 2019.

Edited by:

Cyril Pommier, INRA Centre Versailles-Grignon, FranceReviewed by:

András Zlinszky, Institute of Ecology Research Center (MTA), HungaryEric R. Casella, Forestry Commission England, United Kingdom

Copyright © 2019 Hao, Li, Zhao, Chang and Zhao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Qingrui Chang, changqr@nwsuaf.edu.cn; Pengxiang Zhao, zhaopengxiang@nwsuaf.edu.cn