密度依赖性聚集驱动海岸沙丘空间自组织

《Proceedings of the National Academy of Sciences》：Density-dependent aggregation underlies spatial self-organization in coastal dunes

　　

摘要 在整个海岸生态系统中，栖息地形成生物与风和水的相互作用产生生物地貌反馈，从而形成空间格局化的景观。此类空间格局可通过不同的自组织机制涌现，这对生态系统韧性具有截然不同的影响。本研究检验了发育中的海岸沙丘的空间动态是否更符合尺度依赖反馈（scale-dependent feedback, SDF）还是密度依赖性聚集（density-dependent aggregation, DDA）——后者是一种较近期提出的空间自组织机制。研究人员利用高分辨率航空影像和耦合的高程时间序列，追踪了沙丘草斑块和沙粒堆积在早期沙丘发育过程中的空间动态。沙丘草斑块及其共同形成的沙丘体并未收敛于某一主导尺寸或间距，而是随时间粗化，其自相似尺寸分布向更大尺度偏移。沙粒堆积不成比例地集中在已经密集的植被区域，呈现出一种在临界阈值处改变符号的密度-生长关系。综合来看，这些发现与DDA而非SDF作为发育中海岸沙丘空间自组织的主导机制相一致。由此产生的密度-生长关系意味着对局部干扰的韧性具有阈值依赖性。一旦局部植被密度降至临界点以下，衰退将是自我强化而非自我抑制的，这对如何理解和管理沙丘干扰与恢复具有直接影响。

一、研究背景与科学问题

空间自组织是许多海岸生态系统（如海草床、贻贝床、红树林、盐沼、珊瑚礁和海岸沙丘）形成显著空间格局的核心过程。这些格局源于栖息地形成生物与环境动力（水流、风力）之间的双向反馈，即生物地貌反馈。长期以来，最被广泛研究的自组织机制是尺度依赖反馈（scale-dependent feedback, SDF），它由短距离正反馈（促进）和长距离负反馈（抑制）构成，能够使斑块趋向于一个主导尺寸和间距，形成规则的斑块、条纹或迷宫状格局，例如干旱区植被或贻贝床中的经典图案。

然而，近年来研究者提出了一种替代性机制——密度依赖性聚集（density-dependent aggregation, DDA）。DDA适用于“保守型”生态系统，其中一种移动资源（如沉积物、养分或消费者）在系统内部被重新分配，而净通量很小。在这种机制下，资源会向高生物活性斑块集中，导致大斑块通过吞并小斑块而生长合并，表现为持续粗化而非收敛于固定尺度，类似于材料科学中的相分离（phase-separation）动力学。这一区别至关重要，因为两种机制预测了截然不同的生态系统韧性：SDF驱动的系统在局部扰动后倾向于恢复原有格局，表现出较强的局部韧性；而DDA驱动的系统则表现出阈值依赖性韧性——一旦斑块密度低于临界值，衰退将自我强化，资源被其他斑块捕获，受损斑块难以恢复。

尽管已有建模和实证工作表明海岸沙丘是一个自组织系统，但其景观层面的主导自组织机制尚不明确。因此，本研究旨在通过分析发育早期海岸沙丘的空间格局和时间动态，检验DDA与SDF何者为更优的解释框架。

二、研究方法概述

研究人员选取荷兰沿岸三个发育中的海岸沙丘系统（阿默兰岛Ameland、特塞尔岛Texel和沙引擎Sand Engine）作为研究对象，利用国家航空调查提供的预处理正射镶嵌影像（分辨率0.1×0.1 m）提取沙丘草斑块。对于主研究区特塞尔岛，还使用了市政当局提供的立体摄影生成的年度数字高程模型（2014–2023年，分辨率0.05–0.07 m）。斑块检测采用基于地理空间深度学习平台Picterra的自定义目标检测与特征掩膜模型。数据分析包括：通过互补累积分布函数拟合斑块面积的最优理论分布模型，利用成对关联函数量化空间聚类特征，通过相对高程图（25米邻域半径）提取沙丘体并计算等效半径以追踪粗化动力学，以及将离散斑块转换为连续植被密度栅格（7米半径移动窗口均值）并与逐年高程变化进行耦合分析。

三、研究结果

沙丘草斑块的空间组织

跨站点比较显示，所有站点的斑块尺寸均最优拟合对数正态分布，呈现右偏长尾特征。成对关联函数揭示斑块在所有尺度上均呈聚集状态，未在任何尺度上出现分散信号，表明不存在主导斑块间距。在特塞尔主站点的时序分析中，斑块尺寸分布的互补累积分布函数随时间向右偏移，尾部增厚，表明大斑块的贡献持续增加；成对关联函数的短程聚集减弱且峰值向更大距离移动（从2.3米到4.8米），这与斑块生长和合并导致的粗化一致，而非收敛于固定波长。

沙丘体的粗化动力学

从相对高程图中提取的沙丘体等效半径同样呈现对数正态分布，互补累积分布函数随时间系统性右移，且最大斑块的变化最为显著。粗化速率随时间减缓，后续年份曲线间距缩小。当半径按每年均值归一化后，各年份互补累积分布函数塌缩为单一曲线，表明分布形状保持不变而仅特征尺度增长，即动态标度性。经总覆盖面积校正和不同分析尺度验证后，这些结果依然稳健。经验粗化指数接近相分离理论预期的1/3幂律增长。

海岸沙丘中的密度依赖性聚集

植被密度与高程变化的耦合分析表明，在最早的时间间隔内，高程变化随植被密度增加从净侵蚀转为净堆积，且在低密度区间存在临界转换点。交叉变异函数显示该关系不具有特定空间尺度。在多年累计效应中，最高密度类别的区域捕获了不成比例的沙量，而大量低密度区域始终低于系统平均增长率，空间对比持续增强。这一模式符合DDA的核心预测：移动资源（沙）优先被最密集的植被区域截留，以牺牲稀疏区域为代价。

四、讨论与结论

本研究表明，发育中海岸沙丘的空间自组织过程更符合DDA而非SDF。具体而言：沙丘草的斑块排列和尺寸分布呈现粗化而非波长选择；沙丘体的粗化显示出相分离的特征；沙积累不成比例地集中于高密度植被区，强化了高程对比。值得注意的是，虽然发育中的沙丘并非严格保守系统（沙不断从海滩输入），但其内部再分配速率远超过长期净堆积速率，形成了时间尺度分离，使得沙在最终固定前可被反复重新分配并向最具截留能力的区域集中，从而产生类似相分离的粗化。这暗示DDA的适用条件可扩展到“开放边界”系统。

关于生态韧性，经验密度-生长关系与DDA的理论预测高度吻合：生长率在中等密度阈值处改变符号，低于阈值的区域进入持续衰退。这意味着发育中沙丘的局部韧性是阈值依赖的，常见的人为扰动（如踩踏或小型风蚀坑）可能将植被密度推至临界水平以下，引发持续的衰退。然而，DDA也预示了尺度依赖的韧性：局部释放的资源可在系统别处被捕获，因此系统整体韧性可能仍然较高。

综上，研究呼吁建立针对开放边界生态系统的理论框架，其中生长限制性资源通过内部再分配被生物富集，即使系统总量随时间变化。这一框架不仅有助于解释空间格局的形成，更对指导沙丘修复、保护和韧性评估具有直接的实践意义。