Climate and land use changes impact the trajectories of ecosystem service bundles in an urban agglomeration: Intricate interaction trends and driver identification under SSP-RCP scenarios.

The delivery of ecosystem services (ESs), particularly in urban agglomerations, faces substantial threats from impending future climate change and human activity. Assessing ES bundles (ESBs) is critical to understanding the spatial allocation and interactions between multiple ESs. However, dynamic projections of ESBs under various future scenarios are still lacking, and their underlying driving mechanisms have received insufficient attention. This study examined the Beijing-Tianjin-Hebei urban agglomeration and proposed a framework that integrates patch-generating land use simulation into three shared socioeconomic pathway (SSP) scenarios and clustering analysis to assess spatiotemporal variations in seven ESs and ESBs from 1990 to 2050. The spatial trajectories of ESBs were analyzed to identify fluctuating regions susceptible to SSP scenarios. The results indicated that (1) different scenarios exhibited different loss rates of regulating and supporting services, where the mitigation of degradation was most significant under SSP126. The comprehensive ES value was highest under SSP245. (2) Bundles 1 and 2 (dominated by regulating and supporting services) had the largest total proportion under SSP126 (51.92 %). The largest total proportion of Bundles 4 and 5 occurred under SSP585 (48.96 %), with the highest provisioning services. The SSP126 scenario was projected to have the least ESB fluctuation at the grid scale, while the most occurred under SSP585. (3) Notably, synergies between regulating/supporting services were weaker under SSP126 than under either SSP245 or SSP585, while trade-offs between water yield and non-provisioning services were strongest. (4) Forestland and grassland proportions significantly affected carbon sequestration and habitat quality. Climatic factors (precipitation and temperature) acted as the dominant drivers of provisioning services, particularly water yield. Our findings advocate spatial strategies for future regional ES management to address upcoming risks.