《Research in Cold and Arid Regions》：Spatiotemporal Evolution and Driving Mechanisms of Ecosystem Service “Supply-Demand-Matching” in the Gansu-Qinghai Contiguous Region of the Upper Yellow River
编辑推荐：
位于黄河上游的甘青毗连区作为国家级战略水源涵养与生态安全屏障，其生态系统极为脆弱，并面临气候变暖和人类活动加剧的双重压力。研究生态系统服务（Ecosystem Services, ES）的供需动态对于循证生态管理和可持续资源治理至关重要。本研究整合了InVES
位于黄河上游的甘青毗连区作为国家级战略水源涵养与生态安全屏障，其生态系统极为脆弱，并面临气候变暖和人类活动加剧的双重压力。研究生态系统服务（Ecosystem Services, ES）的供需动态对于循证生态管理和可持续资源治理至关重要。本研究整合了InVEST模型和RUSLE方程，量化了2000至2023年间五种ES——产水（Water Yield, WY）、食物生产（Food Production, FP）、碳储存（Carbon Storage, CS）、土壤保持（Soil Retention, SR）和生境质量（Habitat Quality, HQ）——的年际供给、需求及供需比（Ecosystem Service Supply-Demand Ratio, ESDR）。研究人员应用包括ArcGIS、GeoDa和Canoco在内的空间分析工具，考察了其时序轨迹、空间匹配模式及驱动机制。主要发现如下：WY供给表现出显著的年际波动，整体呈微弱上升趋势，但在2022–2023年急剧下降；FP供给在研究期内增长超过一倍；CS、SR和HQ供给总体保持稳定。WY、CS和HQ的需求呈上升趋势，集中在东部河谷廊道和主要城市走廊，而FP需求持续处于低位。不同服务的ESDR轨迹各异：WY的ESDR波动剧烈；FP的ESDR在2001年骤降，随后约从2004年开始恢复，但到2023年仍略低于2000年水平；CS的ESDR逐渐下降；SR的ESDR总体保持稳定；HQ的ESDR维持高覆盖率，但在城市群区域存在持续赤字。双变量局部空间自相关（Bivariate LISA）分析揭示了结构一致的西-东空间对比：生态完好的高地地区呈现高-低盈余区，而东部河谷廊道则表现为低-高错配区。冗余分析（Redundancy Analysis, RDA）证实了一个双重驱动框架，其中水文气候和地形因子主导供给侧，而人口密度和国内生产总值（Gross Domestic Product, GDP）则支配需求侧，且在整个研究期间，驱动结构从气候-自然主导型演变为自然-社会二元配置型。这些发现为黄河流域差异化的生态修复分区、水资源治理及空间定向生态补偿政策提供了科学依据。
研究背景与现存问题
黄河上游甘青毗连区地处青藏高原与黄土高原交汇地带，是国家重要的水源涵养区和生态安全屏障。然而，该区域生态系统结构简单、抗干扰能力弱，在气候变暖和人类活动加剧的双重压力下，局部出现草地退化、荒漠化和水土流失等生态问题。生态系统服务（Ecosystem Services, ES）研究虽已从概念探讨转向供需评估与错位诊断，但仍存在四方面不足：第一，多数研究分别识别影响供给和需求的驱动因素，却鲜有研究厘清这些因素如何共同塑造供需平衡维度及其贡献随时间空间的演化；第二，许多研究基于横截面统计相关性推断驱动因素，难以捕捉非线性效应与时变交互作用；第三，干旱区生态系统长期、高分辨率时空数据集的匮乏制约了对耦合社会-生态驱动力演化轨迹的追踪；第四，现有研究在供给-需求关系的量化方法、匹配评估指标（如指数法与空间叠置法或聚类方案间的差异）、以及山地流域过渡条件下服务跨空间单元流动性的考量等方面仍存在分歧。为此，研究人员开展了此项研究。
研究内容与结论
研究人员采用InVEST模型和RUSLE方程，结合ArcGIS、GeoDa和Canoco等空间分析与统计工具，系统量化了2000至2023年间五种关键ES——产水（WY）、食物生产（FP）、碳储存（CS）、土壤保持（SR）和生境质量（HQ）——的供给、需求及供需比（ESDR）。研究旨在揭示以下核心问题：各ES的供需关系在时间上如何演变？供给与需求的空间关系如何表现，可识别出何种独特的空间匹配模式？哪些关键自然和社会经济因素解释了这些供需关系的内在时空动态？
研究发现：WY供给呈现强年际波动，2022–2023年急剧下降；FP供给翻倍以上增长；CS、SR和HQ供给基本稳定。WY、CS和HQ需求呈上升趋势且集中于东部河谷与城市走廊。WY的ESDR波动剧烈，FP的ESDR在2001年骤降后逐步恢复但仍未达2000年水平，CS的ESDR逐渐下降，SR的ESDR总体稳定，HQ的ESDR覆盖率高但在城市群区域持续赤字。双变量LISA分析显示稳定的西-东空间分异：生态完好的西部高地多为高-低盈余区，而东部河谷廊道则为低-高错配区。RDA证实了双重驱动框架：水文气候与地形因子主导供给侧，人口密度与GDP支配需求侧，且驱动结构从气候-自然主导向自然-社会二元配置演变。这些结论为差异化生态修复分区、水资源治理及空间定向生态补偿政策提供了科学基础。
关键技术方法
- 1.
生态系统服务量化模型：采用InVEST模型的产水模块和碳储存模块分别估算WY和CS的供给；利用修正通用土壤流失方程（RUSLE）计算SR供给；基于归一化植被指数（NDVI）与县级作物产量统计数据，按像元NDVI占比进行空间降尺度分配FP供给；运用InVEST的生境质量模块评估HQ供给。
- 2.
供需匹配度评估：通过计算ESDR定量表征供给与需求的平衡状态，其中ESDR > 0表示供给盈余，ESDR < 0表示供给赤字。
- 3.
空间自相关分析：运用双变量局部空间自相关（Bivariate LISA）识别四种空间匹配类型：高高集聚（HH）、低低集聚（LL）、低高错配（LH）和高低盈余（HL）。
- 4.
冗余分析（RDA） ：以五种ES的供给、需求和ESDR为响应变量，选取年均降水量（PRE）、年潜在蒸散量（PET）、数字高程模型（DEM）、坡度（SLOPE）、人口密度（POP）、GDP、建设用地比例（CON）、耕地比例（CUL）和NDVI共九个因子为解释变量，通过CANOCO 5.0软件揭示驱动机制。
研究结果
3.1 生态系统服务供需及供需比关系的时空特征
3.1.1 产水服务（WY）：WY供给呈现强年际波动，整体略有上升但末期明显下降，空间上呈南高、东西低的格局。WY需求经历三阶段变化，2010年左右达峰后下降并趋于稳定，东部城市化地区需求较高。ESDR呈“高-低-高-低”阶段性波动，2023年为研究期最低点，东西部边缘地区多呈供给赤字。
3.1.2 食物生产服务（FP）：FP供给稳步上升，总量翻倍以上，高供给区沿农业主带和河谷低地分布。FP需求持续低迷且空间变化小。ESDR呈波动上升趋势，主要由供给侧发展驱动，空间格局与供给一致，南部高海拔高原持续低值。
3.1.3 碳储存服务（CS）：CS供给极其稳定，高供给集中在东南森林带，低值区分布于谷底和黄土丘陵沟壑区。CS需求总体偏低但高需区逐渐扩大，2020年后集中于中部北部城市走廊。ESDR紧密跟随供给空间格局但呈减弱趋势，原因是需求增长压缩了碳汇盈余空间。
3.1.4 土壤保持服务（SR）：SR供给保持稳定，高值聚集于东南山区和上游林带，河岸廊道和黄土丘陵区较弱。SR需求空间异质性更大，高需区集中于北部山麓、主河道和坡耕地。ESDR总体有利，但坡耕地区和下游河谷景观存在持续性低值错配。
3.1.5 生境质量服务（HQ）：HQ供给极其稳定，高值集中于东南山地森林，低值沿河谷和城市密集区分布。HQ需求持续偏低，高需斑块沿城市走廊稳定存在。ESDR大面积高值，但城市群区域存在持续性赤字。
3.2 生态系统服务供需空间匹配特征：双变量LISA分析揭示了五种ES在2000–2023年间的四种匹配类型。总体上，WY受自然地理控制，西北干旱区以LL为主，东南山地多为HL盈余，最尖锐的LH紧张区沿主河道带状分布。CS总体有利，HL在东南林草区占优且逐步扩展。SR和HQ呈现持久的南北对立，北部黄土坡耕地区LH稳定并可扩张，东南山林区HL持久。FP由南部高原和山地景观中的LL主导，年际动态最强。
3.3 生态系统服务供给、需求及供需比的影响因素分析：RDA结果显示，六个代表年份（2000、2005、2010、2015、2020、2023）中轴1和轴2的累计解释力稳定在45%至52%之间。供给侧，PRE、PET和NDVI始终与WY、HQ、CS、SR的供给及其ESDR正相关，构成基础性决定因素；NDVI的贡献从2015年起逐步增强。DEM与SR和CS供给负相关。FP供需向量独立于其他ES。需求侧，POP、GDP、CUL和CON与各需求指标保持正向驱动关系。社会经济与自然驱动向量的结构对立揭示了供需失衡的核心机制。时间上，驱动框架呈现出从气候-自然主导向自然-社会二元配置的渐进轨迹。
讨论与结论
讨论：研究区域的ES供给-需求关系在时空上呈现分异，产生了持续的供需失衡热点。CS和HQ供给的稳定以及SR供给的相对平稳符合缓慢的生物物理动态和植被介导缓冲，而WY对年际水文气候变化的响应更敏感，FP则对长期农业投资和灌溉扩张做出反应。需求增长，尤其是东部和中北部城市及河谷走廊的结构性社会经济集中，增加了生态威胁和消费压力。因此，尽管大面积区域的ES状况改善或稳定，但东部河谷和城市中心的低ESDR区域仍需持续关注。与湿润流域相比，研究区的WY供给具有更强的年际波动性；与黄土高原相比，土壤保持供需动态存在结构相似性但FP供给对灌溉基础设施依赖度更高；与干旱内陆盆地相比，研究区因受季风影响的东南边缘而拥有更多样化的ES供给组合。本研究的空间错配模式与多个流域尺度的发现存在共性，但驱动机制因旅游依赖性流域的季节性脉动与本区域长期结构性因素的差异而不同。基于比率的ESDR指标隐含的服务自由流动假设在地形障碍显著的山区流域存在局限，而基于双变量LISA的空间自相关方法部分弥补了这一不足。
结论：从管理角度看，西部生态完好的高供给区应加强保护并建立生态补偿机制；东部高需求区需严格控制城市扩张并优化产业能源强度；对于持续供需错配的区域，跨区域生态产品交易平台（如水权、碳排放交易）是缓解结构性生态压力的可行机制。
