《Ecology and Evolution》：Species Extinction Curves Under Emergent Spatiotemporal Land-Use Processes and the Role of Protected Areas
编辑推荐：
人为的土地利用变化是生物多样性丧失的主要驱动因素。预测由此产生的灭绝风险对于平衡发展与保护至关重要，但由于缺乏整合人类动态活动与物种空间分布的框架，这仍然充满挑战。在本研究中，研究人员开发了一个空间显式模型，该模型将土地利用进程与物种分布区大小分布相结合，以估
人为的土地利用变化是生物多样性丧失的主要驱动因素。预测由此产生的灭绝风险对于平衡发展与保护至关重要，但由于缺乏整合人类动态活动与物种空间分布的框架，这仍然充满挑战。在本研究中，研究人员开发了一个空间显式模型，该模型将土地利用进程与物种分布区大小分布相结合，以估算物种灭绝。该框架捕捉了随时间推移转化区域在随机涌现、空间扩张和合并的过程，同时考虑了物种的空间分布区结构。研究人员利用它比较了不同土地利用动态、地理梯度和保护区实施情景下的灭绝风险。通过使用特有区-面积关系（EARs）评估即时物种损失，研究发现新的、空间分散的土地利用斑块的频繁涌现，将物种损失转移到栖息地转化的后期阶段，产生了比从单一地点扩张的情景更陡峭的增长。然而，由于多点转化也加速了时间上的土地转化，这种转化模式可能对生物多样性构成更大的实时风险。研究人员还调查了保护区的作用，发现早期实施可以显著减少灭绝，尽管其有效性取决于土地利用模式。这些发现强调了将时空土地利用过程纳入灭绝模型的重要性，并强调了及时的保护策略对减轻生物多样性损失的紧迫性。
论文解读
一、 研究背景、现存问题与研究意义
人为的土地利用变化，如农业开垦和城市化，正以前所未有的速度改变地球生态系统，导致栖息地丧失和破碎化，成为全球生物多样性丧失的主要驱动力。准确预测由此引发的物种灭绝风险，是平衡经济发展与生物多样性保护、制定有效保护策略的关键前提。然而，现有的灭绝风险评估方法，尤其是基于物种-面积关系（Species-Area Relationship, SAR）和特有区-面积关系（Endemics-Area Relationship, EAR）的经典方法，通常采用简化的、静态的栖息地丧失假设。例如，许多研究假设栖息地从单一中心点向外单向扩张，或者维持固定的破碎化程度，未能充分考虑真实世界中土地利用过程的复杂性和动态性。
现实中的土地利用是一个复杂的时空过程：新的转化地块会在不同时间和地点随机涌现，已存在的转化地块会持续扩张甚至相互融合，导致栖息地的丧失和破碎化程度随时间动态变化。这种动态过程如何影响物种灭绝风险的累积曲线，以及保护区的建立能否有效减缓这一趋势，是当前生态学和保护生物学亟需解答的核心问题。本研究正是为了填补这一空白而开展的。
本研究的核心意义在于，开发了一个创新的空间显式模型框架，首次将基于Fisher-KPP反应-扩散方程的动态土地利用过程与基于宏观生态学理论的物种分布格局相耦合。该框架能够模拟土地利用斑块的随机涌现、空间扩张和融合，从而更真实地反映人为干扰下的景观演变。通过该模型，研究人员系统地比较了不同土地利用动态模式（单源扩张 vs. 多源并发扩张）、不同物种地理分布格局（随机分布 vs. 纬度梯度分布）以及不同保护区实施时机对物种灭绝曲线的影响。研究结果揭示了在多源并发土地利用情景下，物种灭绝风险会在后期急剧增加，凸显了传统简化模型可能低估灭绝风险的严重性，并为保护区尽早设立提供了强有力的理论依据，具有重要的科学价值和政策指导意义。
二、 关键技术方法概述
为开展此项研究，研究人员采用了以下主要技术方法：
- 1.
物种空间分布模型：基于Takashina等人提出的宏观生态学框架，通过空间点过程生成物种的分布区。物种分布区大小遵循对数尺度上左偏的经验分布（通过Beta分布生成），并考虑了两种空间配置：随机放置和沿y轴生成纬度梯度（Latitudinal Gradient, LG）。
- 2.
时空土地利用动态模型：采用Fisher-KPP型反应-扩散方程来描述土地利用强度u(x,y,t)的时空演变。该模型包含了土地利用强度的逻辑斯蒂增长、空间扩散以及新土地利用斑块的随机涌现。通过参数τ控制新土地利用事件的发生频率（τ=∞代表单源扩张，τ越小代表多源并发越频繁）。同时，设置了密度依赖系数K来控制转化斑块间的相互作用。
- 3.
物种灭绝评估方法：主要通过特有区-面积关系（EAR）来评估物种的即时灭绝，即当一个物种的整个地理分布区被完全转化为土地利用时，判定该物种灭绝。此外，也探讨了以丧失90%或95%分布区作为灭绝阈值的标准，并与物种-面积关系（SAR）进行了对比。
- 4.
保护区（PA）实施模拟：设定在特定区域被转化后（如20%或40%），在未转化区域中选取物种丰富度最高的网格作为保护区（占所有网格的5%）。保护区内不再允许新的土地利用活动，以此评估其缓冲作用。研究人员还使用成对相关函数g(r)量化了保护区的空间聚集程度。
三、 研究结果
3.1 独立土地利用扩张下的灭绝估计
- •
结果发现：在无密度依赖效应的独立扩张情景下，EAR曲线呈凸形。新土地利用事件发生频率越低（即τ越大，趋向于单源扩张），在同等转化面积下造成的物种即时损失更大。相反，高频率的多源并发扩张（τ较小）在转化初期物种损失较少，但在转化后期会出现更陡峭的上升，其EAR斜率在后期显著增大。这表明，多源并发扩张模式虽然看似延缓了基于面积的物种损失，但其带来的后期灭绝风险更为剧烈。
- •
时间维度分析：多源并发扩张（τ较小）会显著加快土地转化的实际速度，使得达到相同转化比例所需的时间更短。因此，基于面积的“延迟”并不意味着时间上的安全，反而可能导致物种损失在现实中更快发生。
- •
其他情景验证：上述定性趋势在非对称扩散、物种分布呈现纬度梯度等不同场景下均得到验证。当使用95%或90%分布区丧失作为灭绝标准时，也得到了相似的曲线形态，但在纬度梯度场景下，EAR斜率在高转化水平出现了峰值。SAR的结果则呈现出与EAR相反的镜像模式。
3.1.1 保护区在独立土地利用扩张下的影响
- •
结果发现：保护区的设立能有效“弯曲”灭绝曲线，降低其斜率，从而减少物种损失。保护区的引入时间越早，保护效果越好。在单源扩张（τ=∞）情景下，早期设立保护区的优势尤为明显。而在多源并发扩张（τ较小）情景下，由于物种损失本身进展较慢，早期和晚期设立保护区的效果差异相对减小。
- •
空间配置分析：通过成对相关函数g(r)分析发现，保护区的实际空间聚集程度受到土地利用模式的强烈影响。在单源扩张情景下，早期设立的保护区更倾向于聚集在未被侵占的高物种丰富度区域。
3.2 密度依赖土地利用扩张下的灭绝估计
- •
结果发现：当引入密度依赖效应（K=1）时，EAR曲线的总体定性趋势与独立扩张情景相似：多源并发扩张仍会导致后期更陡峭的物种损失。但是，不同τ值之间的差异变小，尤其是在τ=0.05和0.5之间。同时，由于转化面积可能出现非单调变化，导致EAR斜率的波动性增大。
3.2.1 保护区在密度依赖土地利用扩张下的影响
- •
结果发现：即使在密度依赖的扩张模式下，保护区的设立依然能有效降低物种灭绝风险，并且早期设立的效果优于晚期。保护区的空间配置同样受到土地利用模式的影响，但与独立扩张情景相比，其成对相关函数的交点向更短距离移动，表明密度依赖效应改变了保护区的最终空间布局。
四、 总结与讨论
讨论部分总结：本研究通过耦合物种地理分布与时空土地利用动态，证实了动态的土地利用过程会深刻改变物种灭绝曲线的形状。与传统的单源扩张假设相比，多源并发扩张模式会将物种损失推向转化后期，造成灭绝风险的急剧攀升。然而，由于多源扩张同时加快了土地转化的实际速度，这种基于面积的“延迟”并不代表时间上的安全，反而意味着生物多样性面临更快的现实威胁。研究强调了当前大量陆地景观已接近转化后期，物种灭绝风险可能正进入快速上升阶段。此外，保护区即使规模不大，也能有效“弯曲”灭绝曲线，挽救大量物种免于即时灭绝，其保护效益随着设立时间的提前而显著增强。但研究也指出，若考虑灭绝债务，则需要更大规模的保护区才能实现长期保护目标。论文呼吁，应尽快采取行动，正如《全球生物多样性框架》所强调的那样，及时实施基于区域的保护措施至关重要。
研究结论翻译：研究人员通过将动态土地利用过程纳入基于区域的灭绝框架，检验了物种灭绝曲线。在由EAR定义的即时灭绝情景下，观察到物种损失持续加速的趋势。当应用替代的灭绝定义时，出现了不同的模式，突显了灭绝预测对其潜在假设的敏感性。这些发现强调了与灭绝时间和债务相关的巨大不确定性，并指出了在栖息地破坏背景下，需要进一步研究以建立一致的、生态学上合理的灭绝标准。研究结果也强调了保护区在弯曲灭绝曲线方面可以发挥的重要作用。特别是，研究人员发现，保护区的早期实施在减轻即时物种损失方面具有显著益处。为了解决动态土地利用过程在基于区域的灭绝估计中整合有限的问题，研究人员开发了一个空间显式模型，该模型捕捉了景观转化的关键特征，包括转化区域的涌现、扩张和融合。虽然这个框架为分析土地利用变化下的灭绝提供了基础，但它仍然是一个简化的表述。未来的改进应包括纳入多样的土地利用类型和更清晰的灭绝阈值定义，这将增强研究人员理解和应对人为压力与生物多样性丧失之间复杂相互作用的能力。
