《Biology》：Spatiotemporal Dynamics of Habitat Ecological Intactness and Responses of Wild Giant Panda Populations in Sichuan, China
编辑推荐：
中文摘要翻译
由人为干扰引起的栖息地退化是对全球生物多样性的主要威胁，然而大熊猫栖息地生态完整性的时空动态及相关种群响应仍缺乏量化研究。本研究聚焦于大熊猫国家公园四川片区及凉山山系，利用多源人为干扰数据及三次全国大熊猫调查数据，集成模糊叠加分析（fuzzy o
中文摘要翻译
由人为干扰引起的栖息地退化是对全球生物多样性的主要威胁，然而大熊猫栖息地生态完整性的时空动态及相关种群响应仍缺乏量化研究。本研究聚焦于大熊猫国家公园四川片区及凉山山系，利用多源人为干扰数据及三次全国大熊猫调查数据，集成模糊叠加分析（fuzzy overlay analysis）、Theil–Sen趋势分析（Theil–Sen trend analysis）、地理探测器（Geodetector）及Mann–Whitney U检验，分析了栖息地完整性的时空动态、干扰组分贡献及种群响应。2020年栖息地完整性均值为6.70；1980–2020年间完整性缓慢下降，73.7%的栖息地保持稳定。在完整性指数的干扰组分中，土地利用类型与道路建设对空间差异的贡献最大，且任意两个组分交互作用的联合贡献均超过单一组分。大熊猫出现位点的完整性高于系统背景点（六组调查—区域比较中有五组显著），且凉山山系的完整性高于国家公园。然而，完整性与种群密度无显著相关。结果表明，野生大熊猫集中于高完整性栖息地，尽管完整性单独并不能预测种群密度，该研究为以完整性为导向的保护管理提供了实证依据。
论文解读
一、研究背景与科学问题
人类活动如农业扩张、基础设施建设、资源开采和城市化已改变了全球超过70%的地表，造成栖息地丧失、景观破碎化和生物多样性下降。生态完整性指生态系统保留其自然基线且免受人为改造的程度，高完整性生态系统提供不可替代的生态系统服务，并作为濒危野生动物的重要庇护所。大熊猫作为中国西南山地特有的伞护种，经过数十年保护投入，野生种群已有所恢复，IUCN红色名录将其从濒危降为易危。然而，野生大熊猫目前分为33个局域种群，其中18个种群个体数不足10只，面临因栖息地退化和破碎化导致的较高灭绝风险。
已有生态完整性评估方法包括全球尺度的人类足迹指数、生物多样性完整性指数以及景观尺度生态完整性指数等。在熊猫分布区，已有研究评估了大熊猫国家公园及其邻近区域生态完整性的时空变化，但未涉及野生大熊猫种群对完整性梯度的响应，也未包含国家公园外山系的比较。凉山山系位于大熊猫分布区南缘，人地冲突突出，当地社区高度依赖森林资源，放牧压力大，人口增长迅速，家庭收入低，大熊猫种群小、栖息地破碎化严重、累积人为干扰强度高。该区域是否应纳入大熊猫国家公园体系，需要基于完整性的时空动态和种群响应提供科学依据。
本研究提出三个科学问题：1980–2020年间大熊猫国家公园四川片区和凉山山系栖息地生态完整性的时空格局如何？哪些人为干扰组分对完整性空间差异贡献最大及其交互作用如何？大熊猫出现位点和种群密度如何响应完整性空间差异？
二、研究设计与技术方法
研究区包括大熊猫国家公园四川片区（20177 km2）和凉山山系（15722 km2）。研究人员使用了三类数据集：第二次（1985–1988）、第三次（1999–2003）和第四次（2011–2014）全国大熊猫调查的分布数据，分别获得791、3098和1846个大熊猫出现位点；中国土地利用/覆被变化数据集（CNLUCC）30 m分辨率1980、1990、2000、2010和2020年五期数据；以及道路、水电水利工程和矿山等矢量人为干扰数据。
关键技术方法包括：模糊叠加分析整合多源干扰强度构建生态完整性指数；Theil–Sen中值斜率结合Mann–Kendall非参数检验逐像元分析1980–2020年完整性长期趋势；地理探测器（Geodetector）量化各干扰组分层对完整性空间异质性的相对贡献及交互作用；Mann–Whitney U检验比较出现位点与系统背景点的完整性差异；Pearson相关分析检验完整性均值与种群密度的关系。
三、研究结果
3.1 2020年栖息地生态完整性现状
2020年研究区平均栖息地完整性为6.70，72.0%区域完整性≥7.0，22.4%区域≤6.0，受到相对较强的人为改造。完整性随海拔升高而增加。国家公园四川片区平均完整性为7.04±1.61，高于凉山山系的6.26±2.34。3000 m以下国家公园完整性高于凉山山系，3000 m以上凉山山系略高。
3.2 1980–2020年栖息地完整性时空动态
1980–2020年间区域完整性总体缓慢下降，区域均值从1980年的6.73±2.06降至1990–2000年的最低值6.69±2.00，2010年小幅回升至6.70±2.00并保持至2020年，四十年间均值总变幅仅0.04。趋势分类显示73.7%栖息地保持稳定，0.4%呈下降趋势，0.5%呈上升趋势。国家公园内79.0%栖息地保持稳定，凉山山系为66.9%。五个时间点上国家公园完整性均高于凉山山系。
3.3 干扰组分对完整性空间差异的贡献
因子探测器显示各干扰组分对完整性空间变异的贡献排序为：土地利用类型（q=0.582，p<0.001）>道路建设（q=0.368，p<0.001）>水电建设（q=0.124，p<0.001）>水利建设（q=0.059，p<0.001）>矿山开采（q=0.018，p=0.359，不显著）。土地利用类型和道路建设是完整性指数的主导组分。交互探测器显示任意两个组分层的联合分层对完整性空间变异的解释力均超过单一组分，土地利用类型与道路建设的组合产生最高q值（0.768）。
3.4 大熊猫种群对栖息地完整性的响应
三次调查中，凉山山系大熊猫出现位点的完整性均显著高于国家公园内位点。国家公园内，第二次和第四次调查中熊猫出现位点完整性显著高于背景点，第三次差异不显著；凉山山系三次调查中均显著高于背景点。表明大熊猫出现与高完整性栖息地存在空间关联。但局域种群尺度（r=0.09，p=0.76，n=14）和保护区尺度（r=?0.26，p=0.15，n=33）上，完整性与种群密度均无显著相关。
四、讨论与结论
研究人员指出，高完整性栖息地并不必然满足大熊猫全部生存需求，种群密度还受食物资源、坡度、林分结构、保护区管理有效性、历史遗留干扰和种群恢复时滞效应等多因素共同调节。完整性可作为大熊猫栖息地评估的核心指标之一，但不应作为种群密度的唯一预测因子。
管理建议包括：大熊猫分布区基础设施规划应考虑干扰距离衰减关系，新建道路、水电和矿山项目应尽量避开核心栖息地；对关键分布区内低完整性斑块实施针对性生态修复并逐步拆除与保护目标不相容的设施；加强保护区一线执法以缓解放牧和过度采笋等弥散性干扰。若未来将凉山山系纳入国家公园体系，应优先配套社会经济政策并实施多干扰源协调管理。
研究结论：研究区平均栖息地完整性为6.70，72.0%区域≥7.0，总体水平较高但40年间呈缓慢下降趋势；土地利用类型和道路建设对完整性空间异质性贡献最大；大熊猫出现位点与高完整性栖息地相关联，该关联在干扰更重的凉山山系尤为突出；但完整性单独不能预测种群密度。
研究局限性包括：完整性指数依赖遥感可检测干扰，放牧、采笋、采药等弥散性活动未量化，可能导致凉山山系和高海拔区域完整性被高估；五期时间步长仅提供粗略的时间分辨率；三次调查空间精度不同且与最近年份完整性图层匹配，引入不确定性；地理探测器解释层为完整性指数组分本身，q值描述相对贡献而非独立驱动机制；背景点为系统网格中心而非可用性样本，熊猫—完整性关系应解释为空间关联而非已证实的栖息地选择。未来研究应将弥散性干扰纳入完整性评估体系，并将完整性指标纳入大熊猫国家公园及其周边山系的常规监测框架。
