Enhancing high-resolution forest stand mean height mapping in China through an individual tree-based approach with close-range lidar data

Abstract. Forest stand mean height is a critical indicator in forestry, playing a pivotal role in various aspects such as forest inventory, sustainable forest management practices, climate change mitigation strategies, monitoring of forest structure changes, and wildlife habitat assessment. However, there is currently a lack of large-scale, spatially continuous forest stand mean height maps. This is primarily due to the requirement of accurate measurement of individual tree height in each forest plot, a task that cannot effectively be achieved by existing globally covered, discrete footprint-based satellite platforms. To address this gap, this study was conducted using over 1117 km2 of close-range light detection and ranging (lidar) data, which enables the measurement of individual tree heights in forest plots with high precision. Apart from lidar data, this study incorporated spatially continuous climatic, edaphic, topographic, vegetative, and synthetic aperture radar data as explanatory variables to map the tree-based arithmetic mean height (ha) and weighted mean height (hw) at 30 m resolution across China. Due to limitations in obtaining the basal area of individual tree within plots using uncrewed aerial vehicle (UAV) lidar data, this study calculated the weighted mean height through weighting an individual tree height by the square of its height. In addition, to overcome the potential influence of different vegetation divisions at a large spatial scale, we also developed a machine-learning-based mixed-effects (MLME) model to map forest stand mean height across China. The results showed that the average ha and hw across China were 11.3 and 13.3 m with standard deviations of 2.9 and 3.3 m, respectively. The accuracy of mapped products was validated utilizing lidar and field measurement data. The correlation coefficient (r) for ha and hw ranged from 0.603 to 0.906 and 0.634 to 0.889, while the root mean square error (RMSE) ranged from 2.6 to 4.1 and 2.9 to 4.3 m, respectively. Comparing with existing forest canopy height maps derived using the area-based approach, it was found that our products of ha and hw performed better and aligned more closely with the natural definition of tree height. The methods and maps presented in this study provide a solid foundation for estimating carbon storage, monitoring changes in forest structure, managing forest inventory, and assessing wildlife habitat availability. The dataset constructed for this study is publicly available at https://doi.org/10.5281/zenodo.12697784 (Chen et al., 2024).