【摘要】【目的】探究生态系统服务权衡与协同耦合关系的驱动因素及其动态变化特征，对于明确生态恢复背景下社会-生态系统的互馈机制至关重要。【方法】本文以典型植被恢复区延安市为研究区，量化了2000—2020年碳固存、产水量、基流调节与土壤保持服务，采用生态系统服务综合指数（ESI）表征4项生态系统服务耦合关系；基于贝叶斯信念网络模型（BBNs）构建社会-生态驱动因素识别框架，并揭示了植被恢复以来研究区生态系统服务耦合关系的驱动因素及其变化特征。【结果】(1)2000—2020年延安市生态系统服务耦合关系显著改善，ESI高值区由中部向南北扩张，低值区集中在城镇周围以及北部边缘地区；但基流调节高值概率持续降低。(2)降水、植被覆盖度、实际蒸散发与气温对延安市生态系统服务耦合关系变化的综合贡献率达70%以上，且降水和植被覆盖度是主导驱动因素。(3)20年间，降水的贡献率波动减少，植被覆盖度的贡献率持续上升；概率推理表明，气候暖湿化将加剧基流调节下降的风险，而持续的植被恢复也将引起生态系统服务耦合关系优化的阈值效应。【结论】对植被恢复区开展长期监测、提高区域基流调节功能是维持植被恢复成效以应对气候变化和人类活动干扰的优先事项。