中国西南喀斯特地区自然与人为因素对生产-生活-生态功能相互作用的制约

《Agriculture》：Natural and Anthropogenic Factors Condition Production–Living–Ecology Function Interactions in the Karst Region of Southwest China

　　

协调生产-生活-生态（Production–Living–Ecology, PLE）土地利用需要准确了解PLE功能之间如何相互作用，然而大多数研究将这些相互作用视为静态属性，并孤立地报告其阈值。研究人员转而将它们视为由其驱动因素所制约的动态状态。针对中国西南喀斯特地区，研究人员利用地理探测器（Geodetector）识别了PLE功能的主要自然因素和主要人为因素，通过配对级规则为每个功能对固定一个自然和一个认为制约轴，并沿这两个梯度重建了2010年和2019年每个功能对的帕累托前沿（Pareto frontier）；曲线转折点区分了反转型阈值和缓冲型阈值，调节分区由此直接得出。生产-生态和生活-生态在拟合范围内是权衡关系，而生产-生活在大部分生产范围内边界上升。等权重参考最优值下的总效益在所有三个功能对中均上升（0.552至0.636、0.701至0.726以及0.823至0.872），在1000次分块自助法（block bootstrap）重采样中无一下降。边界在约束函数上升和下降的区域之间也存在差异，但这两组在海拔和人口密度方面不平衡，因此这一对比仅作为描述性比较报告。沿海拔梯度，生活-生态相互作用每年显示一个反转型阈值，分别在595米和756米处，其置信区间定义了一个498至939米的预防性管理包络线；沿人口密度梯度，在子流域尺度上没有阈值幸存。该框架为脆弱区域的分异调控提供了可操作的基础。

论文解读：自然与人为因素对中国西南喀斯特地区生产-生活-生态功能相互作用的制约

**研究背景与问题**
土地资源的多功能性和多宜性使生产、生活、生态（Production–Living–Ecology, PLE）功能共存于同一地块。快速工业化与城镇化导致三类用地相互转换频繁，供需平衡被打破，功能间权衡加剧，整体效益下降，生态安全风险升高。已有研究多从驱动因素或相互作用两个侧面展开：地理探测器（Geodetector）用于量化自然与人为因子对功能的解释力，而统计关联、空间相关与聚类分析用于识别权衡或协同。然而，这些方法只能提供空间平均的线性关联，无法刻画关系随功能水平升降的非线性转变，也难以识别关系反转的临界点。同时，多数研究将相互作用视为功能对的内在静态属性，孤立报告阈值，导致诊断与空间调控之间脱节。中国西南贵州—广西喀斯特地区是中国重要生态屏障和石漠化严重区，地形破碎、土壤浅薄、生产生活用地极度稀缺，人类活动与自然背景均呈现强梯度，系统韧性低，一旦相互作用被破坏难以自我恢复。因此，本研究将PLE功能相互作用视为由驱动因素制约的动态状态，构建连接诊断与调控的分析框架，为脆弱区域的分异调控提供依据。论文发表在《Agriculture》。

**研究方法概述**
研究区涵盖贵州9个地级市和广西中西部7个地级市，总面积约2.02×10⁵ km²，基于250 m分析网格共获得3,235,552个有效分析单元，并划分为1017个子流域。首先空间化11个PLE子功能，采用最小-最大归一化统一量纲；利用地理探测器识别每个功能的优势自然因子和优势人为因子，并通过配对级规则确定每个功能对的分析轴；随后从函数值视角提取帕累托前沿（Pareto frontier）曲线，采用分块自助法（block bootstrap）估计不确定性；再从驱动因子视角，在子流域尺度上计算功能对相关系数，沿海拔和人口密度梯度拟合前沿曲线，通过曲线转折点识别反转型与缓冲型阈值；最后进行分段线性回归和多种敏感性检验。

**研究结果**
**PLE功能的驱动因素**
通过地理探测器因子检测，海拔（Ele）在两年中均为生产和生态功能的优势自然因子，也是2010年生活功能的优势自然因子；2019年增强植被指数（EVI）取代海拔成为生活功能的自然轴。人口密度（Pop）在两年中均为三个功能的优势人为因子。所有因子对交互均为增强型，其中海拔与人口密度的交互解释力最强，表明自然轴主要与人为干扰共同发挥作用。

**PLE功能相互作用的效益分析**
从函数值视角，生产-生态和生活-生态在拟合范围内呈权衡关系，而生产-生活的边界在大部分生产范围内上升，不构成权衡。等权重参考最优值下，三个功能对的总效益在2010—2019年间均上升：生产-生活从0.552升至0.636，生产-生态从0.701升至0.726，生活-生态从0.823升至0.872；配对检验显示该增加在1000次重采样中无一下降，且方向对估计设置和权重不敏感。其中生活-生态在2019年出现唯一的内部最优点。

**方向对比分析**
以第三功能为控制，分别对约束函数上升区和下降区拟合可达边界。结果显示，三个功能对的两条边界均发生交叉：上升组在约束函数低端高于下降组，在高端低于下降组。生产-生活的分离幅度最大（0.242），生产-生态和生活-生态分别仅为0.040和0.032。但两组在海拔和人口密度上不匹配，因此该对比仅作为描述性报告，不解释为方向变化的因果效应。

**驱动因子制约下的相互作用分析**
沿海拔梯度，生产-生活始终为协同，但可达协同水平随海拔升高而下降；生产-生态变化幅度极小（0.070—0.075）；生活-生态变化最大，低海拔为明显权衡，随海拔上升转为弱协同，约1600 m以上又回到弱权衡。生活-生态在2010年和2019年各有一个反转型阈值，分别为595 m和756 m，其95%置信区间合并为498—939 m的预防性管理包络线。生产-生活在1037—1547 m出现结构转变，但类型在两年间不稳定，未赋予管理类型。人口密度梯度上，所有10个候选阈值均未通过保留规则，因为曲线变化处位于子流域平均密度范围之上，无法划分两侧子流域。

**讨论与结论**
讨论部分指出，将相互作用重构为受驱动因素制约的动态状态，而非静态标签，是连接诊断与调控的关键。自然因子（海拔）与人为因子（人口密度）以相反方向制约相互作用：随海拔上升，人类活动减弱，生产、生活与生态的冲突缓解；而人口密度上升则加剧生活-生态权衡。阈值类型学为调控提供了直接操作依据：反转型阈值界定关键调整区，缓冲型阈值界定关键保护区。研究还明确了几项不确定性，包括仅有两个时间点、石漠化数据单一时期、道路和POI数据替代，以及人口密度驱动的阈值在子流域尺度无法估计等。结论部分总结了三条主要发现：第一，生产-生态和生活-生态为权衡，而生产-生活边界上升，总效益在等权重参考最优下三年间均上升；第二，海拔和人口密度是主要的自然与人为制约轴，相互作用的方向和强度随驱动因子值系统性变化；第三，海拔驱动的生活-生态反转阈值（595 m和756 m）合并为498—939 m的管理包络线，可作为关键调整区；人口密度阈值未保留。该框架为喀斯特及其他生态脆弱区协调PLE用地和分异国土空间调控提供了可操作工具。