Improving Pinus densata Carbon Stock Estimations through Remote Sensing in Shangri-La: A Nonlinear Mixed-Effects Model Integrating Soil Thickness and Topographic Variables

Abstract

Forest carbon sinks are vital in mitigating climate change, making it crucial to have highly accurate estimates of forest carbon stocks. A method that accounts for the spatial characteristics of inventory samples is necessary for the long-term estimation of above-ground forest carbon stocks due to the spatial heterogeneity of bottom-up methods. In this study, we developed a method for analyzing space-sensing data that estimates and predicts long time series of forest carbon stock changes in an alpine region by considering the sample’s spatial characteristics. We employed a nonlinear mixed-effects model and improved the model’s accuracy by considering both static and dynamic aspects. We utilized ground sample point data from the National Forest Inventory (NFI) taken every five years, including tree and soil information. Additionally, we extracted spectral and texture information from Landsat and combined it with DEM data to obtain topographic information for the sample plots. Using static data and change data at various annual intervals, we built estimation models. We tested three non-parametric models (Random Forest, Gradient-Boosted Regression Tree, and K-Nearest Neighbor) and two parametric models (linear mixed-effects and non-linear mixed-effects) and selected the most accurate model to estimate Pinus densata’s above-ground carbon stock. The results showed the following: (1) The texture information had a significant correlation with static and dynamic above-ground carbon stock changes. The highest correlation was for large-window mean, entropy, and variance. (2) The dynamic above-ground carbon stock model outperformed the static model. Additionally, the dynamic non-parametric models and parametric models experienced improvements in prediction accuracy. (3) In the multilevel nonlinear mixed-effects models, the highest accuracy was achieved with fixed effects for aspect and two-level nested random effects for the soil and elevation categories. (4) This study found that Pinus densata’s above-ground carbon stock in Shangri-La followed a decreasing, and then, increasing trend from 1987 to 2017. The mean carbon density increased overall, from 19.575 t·hm−2 to 25.313 t·hm−2. We concluded that a dynamic model based on variability accurately reflects Pinus densata’s above-ground carbon stock changes over time. Our approach can enhance time-series estimates of above-ground carbon stocks, particularly in complex topographies, by incorporating topographic factors and soil thickness into mixed-effects models.