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 estimation, 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 be effectively achieved by existing globally covered, discrete footprint-based satellite platforms. To address this gap, this study was conducted using over 1117 km² of close-range Light Detection and Ranging (LiDAR) data, which enables the measurement of individual tree height in forest plots with high precision. Besides, 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 (h_a) and weighted mean height (h_w) at 30 m resolution across China. Due to limitations in obtaining basal area of individual tree within plots using UAV LiDAR data, this study calculated 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 large spatial scale, we also developed a machine learning-based mixed-effects model to map forest stand mean height across China. The results showed that the average h_a and h_w across China were 11.3 m and 13.3 m with standard deviations of 2.9 m and 3.3 m, respectively. The accuracy of mapped products was validated utilizing LiDAR and field measurement data. The correlation coefficient (𝑟) for h_a and h_w ranged from 0.603 to 0.906 and 0.634 to 0.889, while RMSE ranged from 2.6 to 4.1 m 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 h_a and h_w 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).