AN Zhiying, SUN Caizhi, HAO Shuai
Exploring the spatial heterogeneity of ecosystem services (ESs) and their driving factors under various interaction patterns is essential for informing sustainable development policies. Using Northeast China as a case study, this research investigates eight key ESs, including water yield (WY), carbon storage (CS), food provision (FP), habitat quality (HQ), soil conservation (SC), wind-break and sand-fixation (WS), water purification (WP) and aesthetic landscape (AL). The study examines the complexity of ESs from three dimensions: individual ES, ES pairs and ES bundles, and further evaluates their spatial heterogeneity and socio- ecological drivers. The results indicate that the spatial distribution of ESs remained relatively stable from 2000 to 2020. During this period, WY and FP increased significantly, CS and HQ remained relatively unchanged. SC, WS and AL followed an “increase-decrease-increase” trend, and WP exhibited a “decrease-increase” fluctuation. Overall, synergistic relationships among ES Pairs were more prevalent than trade-offs. Notably, CS showed trade-offs with over 70% of the other ESs, while HQ exhibited trade-offs with SC, WS, WP, and AL. The FP-HQ synergy bundle, primarily located in the Greater Hinggan Mountains and eastern regions, emerged as the dominant ES bundle. Ecological factors—such as solar radiation, temperature, slope, DEM, and NDVI—exerted a stronger influence on ES patterns than social factors like GDP and population density. Furthermore, these ecological drivers had a greater impact on individual ESs compared to ES pairs or ES bundles. These findings offer valuable insights for policymakers to understand the complex interrelationships among ESs and to design more effective and regionally tailored management strategies.