Analysis of Influencing Factors of Terrestrial Carbon Sinks in China Based on LightGBM Model and Bayesian Optimization Algorithm

With accelerating climate change and urbanization, regional carbon balance faces increasing uncertainty. Terrestrial carbon sinks play a crucial role in advancing China’s sustainable development under the dual-carbon strategy. This study quantitatively modeled China’s terrestrial carbon sink capacity and analyzed the multidimensional relationships between impact factors and carbon sinks. After preprocessing multi-source raster data, we introduced kernel normalized the difference vegetation index (kNDVI) to the Carnegie–Ames–Stanford approach (CASA) model, together with a heterotrophic respiration (Rh) empirical equation, to simulate pixel-level net ecosystem productivity (NEP) across China. A light gradient-boosting machine (LightGBM) model, optimized via Bayesian algorithms, was trained to regress NEP drivers, categorized into atmospheric components (O3, NO2, and SO2) and subsurface properties (a digital elevation model (DEM), enhanced vegetation index (EVI), soil moisture (SM)), and human activities (land use/cover change (LUCC), POP, gross domestic product (GDP)). Shapley Additive Explanation (SHAP) values were used for model interpretation. The results reveal significant spatial heterogeneity in NEP across geographic and climatic contexts. The pixel-level mean and total NEP in China were 268.588 gC/m2/yr and 2.541 PgC/yr, respectively. The north tropical zone (NRZ) exhibited the highest average NEP (828.631 gC/m2/yr), while the middle subtropical zone (MSZ) and south subtropical zone (SSZ) demonstrated the most stable NEP distributions. LightGBM achieved high simulation accuracy, further enhanced by Bayesian optimization. SHAP analysis identified EVI as the most influential factor, followed by SM, NO2, DEM, and POP. Additionally, LightGBM effectively captured nonlinear relationships and variable interactions.