PNAS: 用相变理论破解旱区生态突变

论文题目：Phase transition in a dryland shrub-grass ecosystem: Emergence through 

local trade-offs

论文链接：https://www.pnas.org/doi/epdf/10.1073/pnas.2618656123

发表时间：2026年11月11日

论文来源：PNAS

在复杂系统科学的语境中，一个永恒的问题是：宏观的有序与突变，是如何从微观个体的局部相互作用中自发涌现的？统计物理学中的相变理论（如伊辛模型）给出了经典范例：微观自旋的相互作用规则，在宏观尺度上演绎出磁化、临界涨落与系统相变。

传统上，生态学多采用案例描述的方法，或基于生物与环境之间的宏观反馈（如降水、温度等外部驱动）来研究系统动态。这种视角虽然能刻画系统的整体行为，但在解释宏观空间秩序究竟如何从个体行为中诞生时，往往缺乏自下而上的微观机制推演。这一困境实际上呼应了生态学史上著名的“机体论（Organicist）”与“个体论（Individualistic）”之争——前者将群落视为高度整合的超有机体，后者则强调其只是物种个体独立响应环境而形成的松散集合。而复杂系统的涌现理论，为化解这一对立提供了新的物理视角。

旱区灌-草生态系统，因其结构复杂度可控且存在明显的突变边界，成为了测试这一物理理论的理想试验场。近日，中国科学院地理科学与资源研究所余开亮课题组联合国内外学者在PNAS发表研究，将相变理论与基于交互核（Kernel）的元胞自动机模型结合，构建了一个自下而上的演化框架，并结合野外调查，清晰拆解了从微观植物策略到宏观系统相变的完整生成链条。

图1：旱区灌丛-草地-荒漠相变示意图。

微观权衡：

生态系统的“自旋相互作用”与底层搜索逻辑

在复杂的旱区环境中，植物无法移动，其生存本质上是一个在二维空间中进行资源搜索的最优空间拓扑问题。植物必须通过根系的物理延展，在不确定的环境下最大化资源捕获效率。这就是本研究的第一性起点——受能量约束的“广域-集约（Extensive-intensive）”资源获取权衡（EI trade-off），它本质上对应着两种经典的底层空间搜索算法：

• 灌木的“广度优先”逻辑（广域E型）：为了应对高度不确定的稀疏降水，灌木演化出庞大、稀疏且深扎的根系。这是一种全局覆盖策略，牺牲了单位面积的转化效率，换取了极高的环境容错率与空间接触面。

• 草本的“深度优先”逻辑（集约I型）：草本将根系高度集中于表层，进行高密度的局部抢占。这是一种极致的局部深挖策略，资源转化效率极高，但在长周期的极端干旱下缺乏缓冲。

图2：植物功能性状主成分分析。空间斑块特征投影于性状空间上，呈现出与性状特征的权衡关系。

这一基于空间拓扑的权衡，在野外观测中得到了植物功能性状与景观斑图特征的有力印证。 研究团队沿降水梯度测定了群落功能性状，并进行主成分分析（PCA）。在排除了降水梯度带来的协同变异后，性状空间清晰地分裂为“保守型”与“进取型”两个极端：灌木主要沿保守型轴分布，而草本则沿获取型轴分布。更具突破性的是，斑块大小和最近邻距离与这一“保守-获取”性状轴高度显著相关，其中大而稀疏的灌木斑块始终对应保守型性状，小而密集的草本斑块则对应获取型性状。这提示我们，传统生态学中备受关注的“保守-进取”经济谱，可能只是植物为了扩展空间觅食范围而不得不积累生物量的副产物；真正的底层驱动力，是以空间尺度为核心的EI权衡。

在模型中，这种空间拓扑策略的差异被严格形式化为物种特异的相互作用核（Interaction Kernel），并驱动了以下四个步骤的模拟规则：

1. 状态与扩散：二维晶格上的每个细胞被占据为灌木、草本或裸土。个体以一定概率产生后代并向全局随机扩散（模拟繁殖）。

2. 卷积构建相互作用场：所有植物的分布与各自的“核kernel”进行空间卷积叠加，生成一个连续的相互作用场。灌木的核表现为大半径、低强度的长程负反馈（广域资源获取）和一定的局部正反馈（正相互作用）；草本的核则表现为小半径、高强度的短程负反馈（集约资源获取）。

3. 存活阈值判定：在相互作用场中，每个元胞根据所在位置的效应强度决定存活。当净效应低于物种特异的生存阈值时，面临生存压力（以一定概率死亡）；草本的存活阈值低于灌木，反映了其高资源利用效率带来的竞争优势。

4. 空间更新与涌现：死亡的元胞释放空间，存活者继续迭代。正是由于元胞繁殖过程（正反馈）与长程的资源抑制（长程负反馈）在数学上构成了经典的尺度依赖反馈机制，系统无需任何外部模板，通过局部规则的迭代便自发破缺对称性，最终涌现出图灵斑图（Turing patterns）的宏观相结构。

图3：灌木-草本相互作用核对比与元胞自动机模拟流程图。

相变的两种面貌：外场驱动与无序驱动

当外部环境（如降水梯度/干扰梯度）缓慢改变并逼近临界点时，系统展现了两种截然不同的物理相变形态。

断崖式跳跃：外场驱动的一级不连续相变

在年均降水 180-250mm 的极窄窗口内，旱区灌草系统发生骤然切换。从野外观测来看，总植被盖度虽随降水呈线性上升，但灌木与草本的组成比例却发生了剧烈的非线性骤变，斑块大小和最近邻距离也出现不连续重组；同时，斑块空间分布的偏度（Skewness）在过渡带内急剧上升，这是系统逼近临界状态时发出的典型空间预警信号。

其物理实质在于强烈的局部互斥：灌木通过大半径、低强度的长程根系网络，在斑块间消耗水分，排除了草本的建立；而草本则依靠小半径、高强度的短程密集抢占，抑制了灌木幼苗的更新。当降水（类似于外部磁场H）越过特定阈值时，这种微观层面的拮抗作用累积至临界点，导致系统发生系统级的分岔。这种突变伴随着强烈的双稳态与迟滞效应（Hysteresis）。由于植物斑块一旦形成，其内部的局域促进作用会“锁定”当前状态，系统必须跨越一个极高的能量壁垒才能重新洗牌。这意味着生态修复的严重路径依赖：一旦灌丛锁定景观，即便降水恢复，也很难自发退回草地。

图4:灌丛-草地相变模拟结果与野外观测对比，均表现出相似的灌-草盖度突变与相格局的变化。

温水煮青蛙：无序驱动的二级连续相变

在沙漠边缘，风沙扰动（类似于 Ising 模型中的温度 T 带来的热涨落）正不断侵蚀植被的空间秩序。与降水驱动的灌草切换不同，这种退化并非瞬间的断崖式跳跃，而是一个连续渐变的过程。其微观物理实质在于正反馈修复能力与植被盖度的深度绑定。植被内部存在着由局域促进作用（如遮阴、固沙、水分汇聚）构成的短程正反馈，这种机制原本负责“修复”被破坏的空间秩序。然而，随着沙埋扰动强度不断加剧，植物被掩埋致死，空间秩序被直接破坏；更致命的是，修复能力本身依赖于植被的总覆盖度，盖度越低正反馈修复力越弱，系统抗扰动能力越差，从而陷入恶性循环。当扰动强度逼近临界阈值D_c时，植被盖度呈现出典型的幂律衰减 C～(D_c-D)^β。这标志着系统发生了无序驱动的二级连续相变（类似磁铁在居里点逐渐失去磁性），表征着空间秩序在物理规律支配下连续、渐进的崩塌。

在这一连续相变的边缘，系统演化出了一种特殊的抗逆策略——边界效应（Boundary effect）。在植被与荒漠的交界处，边缘植物通过增强局域促进（如克隆繁殖与地下根系资源共享），同时也会扩大资源获取范围，从而在强烈扰动下维持最后的空间斑块。从相变视角看，这是系统在无序边缘通过重构局部规则来延缓崩塌。这种非线性特征为荒漠化预警提供了明确的物理判据：当植被盖度的下降在逼近临界点时呈现出幂律加速，正是系统即将彻底滑向无结构荒漠相的明确信号，这为在“温水煮青蛙”式的退化中抓住干预的最后窗口期提供了理论依据。

图5：遥感与无人机展示边界效应，表现为面对流沙干扰时植被斑块的扩大与植被盖度的非线性下降。元胞自动机模型通过在kernel 中引入增强的正相互作用与扩大的负相互作用模拟出相同的现象——斑块扩大与盖度的幂律衰减。

相变边缘的脆弱平衡：钉扎效应与精灵圈的出现

在一级相变的临界阈值附近，系统往往会展现出“亚稳态（Metastability）”，如同物理学中零度以下却依然保持液态的“过冷水”。在旱区灌草过渡带观察到的“精灵圈（Fairy circles）”，正是这种亚稳态的生态具象化。而解释其长期存在的物理机制，正是“钉扎效应（Pinning effect）”。在凝聚态物理中，当相变发生（如磁畴翻转）时，材料内部的缺陷或杂质会阻碍相界的移动，将相界“钉”在某处，使系统在去除外力后仍能局部维持旧相。在旱区生态系统中，这种“缺陷”就是灌木斑块中心的极端竞争。

模型推演揭示了精灵圈的动力学机理：当降水增加，草本表层水分利用加剧，灌木斑块的资源获取半径被迫收缩。斑块中心因竞争最激烈而“空心化”，但这恰恰形成了一个空间“缺陷”——空心环汇集了稀缺的径流水，刚好能维持外围灌木环的存活。这个空心结构就像一个物理上的“钉子”，将灌木相界死死钉扎在草本基质中，阻止了系统完全滑向草地相。只有当外部环境提供足够大的驱动力，越过钉扎势垒时，这个脆弱的亚稳态才会最终崩溃，系统才彻底完成相变。

图6：精灵圈模拟与亚稳态示意图。中国北方旱区首次记录到精灵圈格局。

范式跨越：走向严密的交叉科学

这项研究的核心并非仅仅借用了“相变”的词汇，而是在数学结构与物理机制上，将统计物理的经典 伊辛Ising 模型与旱区生态系统完成了一次严密的跨尺度映射。这种映射跨越了物理实体与生物个体的边界，确立了微观生态规则与宏观相变动力学之间的同构关系：

 图7：伊辛模型与灌草相变模型同构性示意图。

• 微观单元（Microscopic unit）：Ising 模型中的微观“自旋（Spin，向上或向下）”，对应着二维晶格上占据一个位置的植物状态（灌木、草本或裸地）。

• 微观相互作用（Microscopic interaction）：物理上的自旋对齐（Spin alignment），即相邻自旋倾向于同向排列，在生态系统中对应着局部竞争排斥导致的优势固化。具体而言，当灌木形成局部优势时，其广域根系的水分消耗会强烈抑制斑块内的草本生长；当草本形成优势时，其高密度的空间抢占则会直接抑制灌木的更新。正是这种微观层面的互斥，驱动系统涌现出了类似铁磁畴（ferromagnetic domains）的宏观空间斑图（Turing patterns）。

• 宏观序参量（Order parameter）：伊辛模型中的“磁化强度M”（平均自旋朝向），在生态模型中对应着灌木的优势度/盖度。它度量了系统宏观相的有序程度。

• 控制参量（Control parameters）：Ising 模型中的外部磁场 H 对应着降水，它作为一种外部偏好力，驱动系统在灌丛与草地之间切换；伊辛模型中的温度 T 则对应着外部的风沙扰动 D，它作为热涨落，不断破坏空间结构的秩序。

基于这种严密的同构映射，系统在临界点附近展示了两种截然不同的物理相变动力学：

• 外场驱动的一阶相变：当外部降水越过临界阈值时，系统发生不连续的“磁化翻转”。模型从两种不同初始状态出发，成功复现了真实的迟滞现象（Hysteresis）——系统被锁定在当前状态，必须跨越极大的势垒才会发生状态切换。

• 无序驱动的二阶相变：当风沙扰动加剧时，类似于温度升高破坏自旋排列，植被盖度在逼近临界点时呈现出典型的幂律衰减  ，对应着向无结构荒漠相的连续崩塌。

面对生态学难以进行全局控制实验的固有局限，这项研究展示了复杂系统范式的一种严密路径：通过确立受约束的微观规则，借助形式化模型进行演绎，并与野外多尺度观测相互印证，同样能对宏观生态突变给出确切的机制解释。更重要的是，这一自下而上的框架为生态学中经典的“机体论与个体论”之争提供了一个基于物理机制的出口。宏观的“机体式”整体秩序（如规则斑图、系统双稳态与突变阈值），并非源于某种神秘的超有机体力量，而是完全由“个体式”的微观植物策略在空间相互作用中，于临界点附近被放大而涌现出的集体行为。

从生命底层的搜索策略，到植物斑块的空间分布，再到宏观的系统相变，相变理论在此成为连通微观与宏观的物理语言。这不仅为旱区生态系统的临界转变与早期预警提供了可推演、可检验的分析框架，揭示了精灵圈与边界效应的动力学本质，也为复杂系统科学跨越实体边界、向其他生物学聚合系统乃至社会生态系统延伸提供了极具参考价值的范式。它表明，在能量约束与空间相互作用的底层逻辑上，生物与非生物系统共享着同样深邃的数学与物理法则。

作者简介：

孙靖尧（论文一作），中国科学院地理科学与资源研究所特别科研助理，研究方向：旱区生态系统动力学、斑图动力学、生态相变理论。学术理念：致力于将复杂系统科学的“自下而上”建构思想引入生态学，探究微观尺度的局部相互作用如何驱动宏观生态格局的涌现与临界突变。代表性成果发表于 PNAS、Soil Biology and Biochemistry、Ecological Modelling 等期刊。

余开亮（通讯作者），中国科学院地理科学与资源研究所研究员，国家引才计划青年项目入选者，研究方向：全球变化生态学、旱区生态系统稳定性、生态系统碳循环。学术简介：长期从事生态系统生态学研究，关注气候变化与人类活动对生态系统结构与功能的影响。相关研究成果以第一或通讯作者身份发表于PNAS、Nature water、Nature Communications、Science Advances 等国际顶级期刊。

生态学本身就是一种充满隐喻的、整体化的自然科学方法论，对人类社会的组织、社群、各类系统都提供了观摩自然界的“模型”和“理论”照鉴。理论生态学是理论中的理论、方法中的方法，因此对这一领域的熟悉和梳理，对更系统地走进生态学的思想宝库、熟练分析和研究工具，具有基础而重要的意义。分享理论生态学框架与前沿、对理论生态学进行祛魅和进一步发展其与相关领域的连接，是我们发起举办这次读书会的初衷。“生态学欢迎所有人”——人生而来源于大自然，生活在各种的生态系统之中，终将归宿于大自然，取材自然、问道天地，遵循和应用自然的秩序和规律。相信理论生态学的读书会会给所有参与者带来一次很有启示意义的学习和共研旅程，从方法论的角度深度理解生态系统更本质的规律和特征。

集智俱乐部联合北京林业大学大学副教授李周园，普利高津奖章得主、Towson大学Brian D. Fath教授以及北京大学理论生态学课题组博士研究生于越共同发起「理论生态学读书会」，旨在深入探讨理论生态学的基础思想与前沿进展，通过分享经典文献与最新研究，促进对生态学复杂性、共存机制及生态系统动态的理解，推动理论生态学与实际生态问题的连接与创新。