Spatiotemporal Dynamics and Future Climate Change Response of Forest Carbon Sinks in an Ecologically Oriented County

Research on forest carbon sinks is crucial for mitigating global climate change and achieving carbon peaking and neutrality. However, studies at the county level remain relatively limited. This study utilized multi-source remote sensing data and the Carnegie–Ames-Stanford Approach (CASA) and soil respiration models to estimate the forest net ecosystem productivity (NEP) in Taoyuan County from 2000 to 2023. The spatiotemporal differentiation was analyzed using seasonal Mann–Kendall tests, Theil–Sen slope estimation, and standard deviation ellipses. The forest NEP for 2035 was predicted under multiple climate scenarios (SSP1-2.6, SSP2-4.5, SSP5-8.5) by applying a discrete coupling of the Patch-generating Land Use Simulation (PLUS) model, incorporating territorial spatial planning policy, and using the CASA model. The results indicated that the Taoyuan County forest NEP exhibited a fluctuating upward trend from 2000 to 2023, with higher (lower) values in the west/south (east/north). Under future warming and humidification, the overall forest NEP in Taoyuan County was projected to decrease by 2035, with predicted NEP values across scenarios ranking as SSP5-8.5 > SSP1-2.6 > SSP2-4.5. The findings offer practical insights for improving local forest management, optimizing forest configuration, and guiding county-level “dual-carbon” policies under future climate and land use change, thereby contributing to ecological sustainability.