《Journal of Theoretical Biology》：Vegetation pattern formation induced by local growth outpacing susceptibility to non-local competition
编辑推荐：
干旱和半干旱生态系统中规则的植被斑图是一种广泛存在的现象，通常被认为是植物自组织的结果。传统模型多关注于非局部竞争引发的排斥区机制。本文作者提出了一个包含局部生长、非局部竞争和生物量扩散的通用框架，并进行了系统性分析。他们识别出两种可引发图灵不稳定性的机制：文献中已记载的、由植被斑块间累积竞争压力导致的排斥区机制，以及一种新颖的、在均匀平衡点附近局部增长超越竞争敏感性的机制。通过线性稳定性分析和两个基准模型的数值模拟验证，该研究不仅建立了不稳定性发生的条件，还提供了清晰的生态学解释。这项研究拓展了我们对纯竞争驱动下植被空间自组织模式形成机制的理解，为更准确地预测生态系统动态提供了理论工具。
在广袤的干旱与半干旱地区，大自然常以令人惊叹的规律性“作画”——疏密相间的植被条带、间隙均匀的植被斑块，如同大地的指纹。这些被称为“植被斑图”的现象，不仅是生态学中一个迷人的谜题，更是生态系统功能与稳定性的关键指示器。数十年来，科学家们一直试图解开这些图案背后的形成之谜。主流理论认为，这是植物之间“尺度依赖反馈”的结果，即近距离的促进（如遮荫保水）与远距离的竞争（如根系争夺水分）共同作用，驱动了空间自组织。基于此，反应-扩散模型等数学工具被成功用于再现这些斑图。
然而，一个核心问题依然悬而未决：在纯粹的竞争性相互作用下，是否还存在传统“排斥区”机制以外的、能够驱动植被形成规则斑图的其他生态“引擎”？传统上，许多模型依赖于非局部竞争核函数的傅里叶变换出现负值（这通常对应着植被斑块周围明确的“排斥区”）来产生图灵不稳定性。但这是否是唯一途径？为了系统性地探索竞争驱动下斑图形成的全部可能性，并明确其生态学含义，来自荷兰莱顿大学数学研究所的Jelle Van Der Voort、Ricardo Martinez-Garcia和Arjen Doelman开展了一项研究。他们的成果以“Vegetation pattern formation induced by local growth outpacing susceptibility to non-local competition”为题，发表在《Journal of Theoretical Biology》上，为我们理解植被空间自组织提供了新的视角。
为了回答上述问题，研究人员主要运用了理论建模与数学分析相结合的方法。首先，他们构建了一个通用的单组分建模框架，该框架整合了局部增长函数g(u)、非局部竞争敏感函数s(u)、竞争压力函数c(u)、表征竞争空间范围的内核函数φ(x)以及描述生物量扩散的扩散项。通过对该模型进行线性稳定性分析，推导出均匀稳态解发生图灵不稳定性的充要条件。研究特别关注了具有紧支集（如方顶、抛物线、余弦、三角核）的非负竞争内核。此外，他们还通过数值模拟，对两个具体的基准模型（非局部Fisher-KPP模型和一个具有非线性竞争函数的模型）进行了计算，以验证理论分析结果并展示超越线性近似范围的斑图动力学。
2. 一个通用的单组分模型
本研究提出了一个通用框架来描述植被的空间动力学，其核心方程包含了局部增长、非局部自限竞争和生物量扩散。模型中的函数g(u), s(u), c(u)被赋予了明确的生态学意义，分别对应局部生长过程、对竞争压力的敏感度以及施加的竞争压力。竞争的空间结构由非负的核函数φ(x)描述，并假设其具有紧支集和对称性。生物量扩散（如种子传播）则用扩散项简化表示。该框架概括了文献中许多非局部模型。
3. 图灵斑图的萌生
3.1. 线性稳定性分析
研究首先分析了模型均匀稳态解的存在性与稳定性。通过傅里叶模态展开和线性化，得到了扰动增长率ω(k)的色散关系。图灵不稳定性要求均匀态在时空均匀扰动下稳定（ω(0)<0），但在特定空间波长扰动下失稳（存在k>0使ω(k)>0）。分析导出了失稳的临界条件，该条件联系了局部增长与竞争敏感性的相对变化率[d/du (g/s)]、竞争压力的变化率c‘以及核函数傅里叶变换φ?(k)的最小值。当扩散系数D小于一个由模型参数决定的最大值Dmax时，空间不稳定性就会发生。
3.2. 图灵不稳定性的机制
线性稳定性分析的关键不等式揭示，图灵不稳定性可以通过两种不同的机制实现：
机制I（排斥区机制）：当竞争核的傅里叶变换φ?(k)在某些波数下为负值时，即使d/du (g/s)为正值（即局部增长相对于敏感性占优），只要满足特定不等式，不稳定性仍可发生。这对应于生态学上熟知的“排斥区”，即植被斑块周围因强竞争形成的裸地区域。
机制II（局部增长超越敏感性机制）：这是一种新发现的机制。当核函数的傅里叶变换φ?(k)始终非负时（例如三角核），图灵不稳定性仍可能发生，但此时要求d/du (g/s)必须为负值。这意味着在均匀平衡态附近，局部增长速率g(u)相对于竞争敏感性s(u)增长得更快。从生态学角度解读，即植被的局部增长能力“超越”了其对外部竞争压力的敏感度，从而在空间上获得了优势。
4. 基准模型的数值模拟
为验证理论并展示机制，研究者对两个模型进行了数值模拟。
4.1. 非局部 Fisher-KPP 模型
该模型是通用框架的一个特例，其中g(u)=u(1-u), s(u)=u, c(u)=u。使用方顶核进行模拟，结果表明系统经历了一个超临界的图灵分岔。当扩散系数D降低穿过临界值Dc时，均匀态失稳，形成稳定、鲁棒的空间周期斑图，其波长与线性理论预测的临界波数kc一致。这证实了机制I在传统模型中的作用。
4.2. 具有非线性竞争函数的模型
此模型设定g(u)=u, s(u)=u/(1+u), c(u)=u，其竞争压力函数c(u)为线性，但敏感函数s(u)饱和。研究发现，即使采用傅里叶变换恒为非负的三角核，系统也能产生图灵斑图。线性稳定性分析证实，此例中正是机制II在起作用，因为d/du (g/s)在平衡点处为负值。数值模拟再次展示了超临界分岔和稳定斑图的形成，有力支持了新机制的存在性。
本研究通过构建并分析一个高度概括的植被动力学模型，系统性地探寻了纯非局部竞争驱动下图灵斑图形成的潜在机制。研究不仅重新确认了经典的、与排斥区相关的机制（机制I），更重要的是，首次揭示并阐释了一种全新的机制（机制II）：即当局部增长在均匀态附近超越竞争敏感性时，同样可以引发空间不稳定性并导致规则植被斑图的形成。这一发现突破了“核函数傅里叶变换必须变号”的传统认知，拓宽了竞争驱动自组织的理论边界。
该研究的核心意义在于其普适性和解释力。通过使用未指定具体形式的函数，研究者建立了一个可用于检验各类生态假设的框架。明确的生态学解读（如将s(u)和c(u)与植物的效应和响应性状关联）使得理论条件能够转化为可检验的生态学预测。这有助于弥合数学模型与实证生态学之间的鸿沟。例如，新机制提示我们，在某些生态系统（如竞争敏感度随生物量增加而饱和的植被）中，即使竞争影响在空间上平滑衰减（无明确排斥区），也可能观察到自组织斑图。这鼓励生态学家在野外调查和实验设计中，更细致地量化植物的生长、竞争敏感度和竞争压力性状，以甄别不同斑图形成机制的主导作用。
总之，这项研究深化了我们对干旱区生态系统空间自组织基础原理的理解，为预测气候变化或人类干扰下植被斑图的动态响应提供了更坚实的理论基础，并在理论生态学与数学生物学领域贡献了一个关于模式形成机理的深刻案例。
