Spatial-Temporal Variation and Driving Forces of Carbon Storage at the County Scale in China Based on a Gray Multi-Objective Optimization-Patch-Level Land Use Simulation-Integrated Valuation of Ecosystem Services and Tradeoffs-Optimal Parameter-Based Geographical Detector Model: Taking the Daiyun Mountain’s Rim as an Example

Exploring and predicting the spatiotemporal evolution characteristics and driving forces of carbon storage in typical mountain forest ecosystems under land-use changes is crucial for curbing the effects of climate change and fostering sustainable, eco-friendly growth. The existing literature provides important references for our related studies but further expansion and improvements are needed in some aspects. This study first proposed an integrated framework comprising gray multi-objective optimization (GMOP), Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST), the Patch-level Land Use Simulation Model (PLUS), and optimal parameter-based geographical detector (OPGD) models to further expand and improve on existing research. Then, the integrated model was used to analyze the spatial–temporal variation in land-use pattern and carbon storage at the county scale in China’s Daiyun Mountain’s Rim under four scenarios in 2032, and analyze the driving force of spatial differentiation of carbon storage. The results indicated that (1) land-use change primarily involves the mutual transfer among forest, cultivated, and construction land, with approximately 7.2% of the land-use type area undergoing a transition; (2) in 2032, the natural development scenario projects a significant reduction in forest land and an expansion of cultivated, shrub, and construction lands. Conversely, the economic priority, ecological priority, and economic–ecological coordinated scenarios all anticipate a decline in cultivated land area; (3) in 2032, the natural development scenario will see a 2.8 Tg drop in carbon stock compared to 2022. In contrast, the economic priority, ecological priority, and economic–ecological coordinated scenarios are expected to increase carbon storage by 0.29 Tg, 2.62 Tg, and 1.65 Tg, respectively; (4) the spatial differentiation of carbon storage is jointly influenced by various factors, with the annual mean temperature, night light index, elevation, slope, and population density being the key influencing factors. In addition, the influence of natural factors on carbon storage is diminishing, whereas the impact of socioeconomic factors is on the rise. This study deepened, to a certain extent, the research on spatiotemporal dynamics simulation of carbon storage and its driving mechanisms under land-use changes in mountainous forest ecosystems. The results can serve to provide scientific support for carbon balance management and climate adaptation strategies at the county scale while also offering case studies that can inform similar regions around the world. However, several limitations remain, as follows: the singularity of carbon density data, and the research scope being confined to small-scale mountainous forest ecosystems. Future studies could consider collecting continuous annual soil carbon density data and employing land-use simulation models (such as PLUS or CLUMondo) appropriate to the study area’s dimensions.