形成与自组织13 分钟

相似降雨条件下,森林与稀树草原的边界由什么维持?

相似年降雨量不能独自决定森林与稀树草原的分界,季节性、土壤、火、种子来源和植被空间结构都参与其中。本文结合火处理研究、2023和2026年空间模型及遥感分析,用连通格网、边界传播和空间混合三个实际计算解释:覆盖比例相同不代表火的传播相同,局部双稳不保证已有边界静止,空间双峰也不直接证明同地点有两个稳态。反馈机制已有具体依据,但自然景观的恢复路径与阈值仍需地点和时间证据。

年雨量相似,树木经历的环境仍可能很不同

看一张年降雨地图,有些相近雨量的地区长着密闭森林,另一些却是草本与散生乔木共同构成的稀树草原。直觉上容易出现一个疑问:同样有这么多水,为什么不能都长成森林?

年总量只包含部分信息。同样的雨量可以分散在全年,也可以集中在短雨季;土层厚度、排水、地下水和坡位还会改变根系真正能利用的水。树木经历的干旱持续时间与草本燃料的干湿状态,都可能因此不同。1

与此同时,植被会反过来改变环境。树冠遮光,影响林下草本;草本形成可燃燃料;火改变树木的更新;树木能否长大,又改变下一轮火的传播条件。

所以,本题既要研究这些内部反馈,也要核对那些还没有被“相似降雨量”控制住的地点差异。这里的稀树草原是一个具有自身物种与功能的系统,不能仅凭树少,就预先解释为退化森林。

树能活下来,与森林能压住草火,是两个层次

在许多易燃草本丰富的景观中,火会反复烧掉幼树的地上部分。树未必彻底死亡,却可能不断从低处重新萌生,长期难以进入高大乔木阶段。

较长的无火时段,可以让部分树木积累树皮、增加高度或建立更强的恢复能力。但某棵树能够抵抗后来的火,和整个景观能够抑制草火,是两件不同的事。

后一个过程需要较大的树冠覆盖和相应林下环境。树冠形成后,遮光可能减少喜光草本;较湿的微环境与较低风速也会改变可燃性。燃料连续性下降,火更难进入或持续传播,树木又获得进一步生长的机会。24

另一条循环则可以维持开放景观:草本支持火,火限制部分树木长大,开放冠层继续支持草本。物种耐火性、资源和水分会改变每一环的强度。

这些循环能够提供历史依赖的机制。但“存在正反馈”还不等于已经证明两种长期稳态。反馈也可能只是让恢复变慢、让边界更清楚,或者放大一个持续存在的外部环境差异。

火的实验证据支持了什么

Gold等在2023年汇总了四个大陆30个火处理地点,包括20个稀树草原和10个森林地点。总体上,有火处理中的草本数量和草覆盖高于无火处理,而且响应程度随生物群落类型和气候而变。2

这不是单纯把卫星图上的“有火”和“有草”做相关。研究纳入了有燃烧与未燃烧对照的实验,最低持续三年;其中可以直接比较草覆盖的子集包括12个地点的239个样地。其他地点使用了生物量、覆盖率或密度等不同指标,再以响应比比较。

结果与火—草本反馈相容:抑制火之后,有些稀树草原的草本会减少;在原有森林中增加火,则未必产生同样大的草本响应。火的后果取决于原来的植被结构和气候。

但原研究也明确说明,它没有直接测试迟滞。要做到这一点,需要例如先让抑制火后的植被转变,再恢复原先的火条件,观察新状态是否仍能自行维持。持续改变火处理,本身就改变了外部条件,不能单凭最后不同的植被,证明同一组条件下存在两个稳态。

此外,研究排除了替代整个林分的树冠火,地点并非都位于同一气候与土壤环境。它没有覆盖所有可能的严重森林火灾,也没有把各种土壤、物种和微气候机制逐一分离。较可靠的结论是火能显著改变植被,且这种作用具有条件性;“整个景观必定有两个不可互相恢复的结局”仍需另外取证。

同样多的草,连成什么样也很重要

火需要接触到能继续燃烧的材料。只知道草占了多少面积,还不知道这些燃料是连成大片,还是被树木、湿地等分割开。

本篇用一个自选的60×60格网做了对照:两幅景观都含1620个草地格,占45%;一幅随机散布,另一幅连成一整片。暂时设定火只能沿上下左右相邻的草地传播,而且一旦接触草地就能继续传播。

如果等概率选一个草地格作为起点,在这一特定随机布局里,平均可达约22.7个草地格;在连成片的布局里,可达全部1620格。

自选格网的草地可达性。两种布局草地比例相同,但连通结构不同;这不是实地火灾概率,也未考虑风、湿度或跨树木传播。
自选格网的草地可达性。两种布局草地比例相同,但连通结构不同;这不是实地火灾概率,也未考虑风、湿度或跨树木传播。

这并不是说自然火灾会固定相差七十多倍。数值取决于格网、布局和传播规则。本例只把一个事实算清楚:面积比例不能替代连通性。

若继续允许火接触到的边缘树木有一定概率受损,森林损失还要看燃烧草地与树木接触了多少边界。一个草地团块大,不代表它的每一段周边都一样长;许多小团块也可能创造很长的树—草接触边界。

因此,单看总草量、总树量或总边界长度中的一个,都可能遗漏重要信息。

边缘既是扩张通道,也是受火接触的位置

Wuyts和Sieber在2023年把这一点落实到一个同时模拟森林、草地、燃烧和灰烬的空间模型中。3

模型里的树木可通过邻近扩散以及较弱的远距离补充建立新斑块;火则从草地开始,较容易经过草地,也能以较低概率进入森林边缘。快速的火与缓慢的植被变化相互塑造景观。

他们提出一个可解释的收支关系:森林向周边扩张的机会,随森林与草地接触的边缘增加;森林遭受火损失的机会,则还要按这些边缘连接到多大的草地团块加权。

这意味着,边缘长既可能帮助树木扩张,也可能增加火暴露。若森林把草地分成许多互不连通的团块,一次点火难以触及全部边缘;若草地连成一个大团块,单次火事件就可能接触更广的森林边界。

研究进一步展示,空间结构会改变宏观非线性关系。把树与草当作充分混合、仅用面积比例相乘的平均场近似,在一些接近均匀的状态下有用,但对聚集明显的景观和变化门槛可能产生较大偏差。

这种关系为实地研究提供了更具体的观测对象:除树冠面积之外,还应测草本连通、森林周长、边缘相连的燃料面积与种子来源。但模型提出的结构指标,仍需用真实火灾与恢复记录校准,不能仅凭一张地图就无条件读出恢复能力。

2026年的模型提醒:看几年,和看一场火,可能得到不同印象

Shen、Levin和Patterson在2026年把火、草本恢复与森林扩张放进具有空间传播核的随机跳跃模型,明确保留了燃烧与火后暂时缺乏燃料的状态。4

火可以在很短时间内传播并停止;草本在更长时间内重新形成可燃层;森林扩张和死亡则更慢。在同一过程中,短时间看见的可能是一场大火,按较长时间平均后却只剩下某种低树冠覆盖。

他们发现,空间模型可以长时间停留在与非空间平均模型明显不同的状态:草先恢复,燃料逐渐变得连续,之后又发生较大的火。把这些阶段平均掉,可能得到一个高灰烬覆盖的平衡,却难以描述实际模拟中的事件间隔与短期变化。

这个模型有几个需要保留的条件。它除了显式的局部火传播,还预设了随局部草量变化的S形可燃性函数;不能把其中所有非线性都称为从零生成。当前主要模拟也没有包含完整的降雨季节性、放牧、土壤差异和树木生活史。

模型之间的恢复边界也不同。2023年的模型包含外来种子或种子库带来的自发树木补充;2026年本文所用版本依赖已有树木向外扩展。如果这个模型内森林已经完全归零,又没有外来树木补充,就不会自行再生森林。这是转移规则带来的限制,不能自动解释为现实土地绝对不可恢复。

有限时长的随机模拟能显示长久维持的不同状态和转变路径,但不能仅凭两段轨迹就证明自然景观会永久停在两个状态中。

局部有两个稳态,边界就一定不动吗

不一定。两个相邻地区会交换种子、承受从另一侧来的火,也会改变彼此的微环境。孤立斑块的稳定性,与一条已经存在的边界能否停住,需要分别研究。

本篇用一个预设双稳结构的无量纲方程检查这个区别:

$$ \frac{\partial u}{\partial t} =\frac{\partial^2u}{\partial x^2}+u(1-u)(u-a), \qquad 0<a<1. $$

$u$表示示意树木状态,空间项概括邻近影响;它没有拟合某个真实森林。在没有空间作用时,$u=0$与$u=1$都稳定,中间的$u=a$不稳定。

现在让左侧为高树木状态、右侧为低树木状态,并放入一条连接它们的边界。对这一方程,可以解析核对边界速度:

$$ c=\frac{1-2a}{\sqrt{2}}. $$

本篇实际数值求解了$a=0.3,0.5,0.7$三种条件,得到向右移动、静止和向左移动,与解析方向及速度一致。更细的网格和时间步长进一步减小了数值误差。

自选空间方程。左图显示相同类型的局部双稳结构可对应不同方向的边界运动;右图展示地点条件连续变化时,从两侧出发的边界趋近同一位置。
自选空间方程。左图显示相同类型的局部双稳结构可对应不同方向的边界运动;右图展示地点条件连续变化时,从两侧出发的边界趋近同一位置。

这个结果没有消除原方程中的两个均匀状态。它说明的是:一旦已经形成高低状态相邻的边界,就不能单凭“局部有两个吸引状态”判断这条边界一定静止。

如果让环境参数随位置变化,使左侧更利于高树木状态、右侧更利于低树木状态,边界可以被限制在环境梯度中的某个位置。本篇从两个不同初始边界出发,也实际算出了这种趋近。

真实森林的边界还可能受地形、河道、物种扩散距离和弯曲形状影响。上述方程预设了双稳反应项,只用来说明空间后果;其中的0.5和速度值不能换成真实世界的雨量或推进公里数。

为什么一张“双峰”分布还不能确定同地点双稳

一种常见证据是:在相似雨量下,低树冠覆盖和高树冠覆盖都很多,中间较少。它与双稳态相容,却不是双稳态的独有结果。

本篇另建了一个反例。每个地点都只有一个稳定覆盖值$m$,并按

$$ \frac{\mathrm du}{\mathrm dt}=m-u $$

向它恢复。一类地点的$m$接近0.25,另一类接近0.75。无论从高覆盖还是低覆盖开始,每个地点都回到它自己的唯一状态。

把这两类地点混在一起统计,覆盖分布就有两个峰。再对这个分布取负对数,也会出现类似“双谷”的图形。但模型中没有任何一个地点具有两个吸引状态。

自选空间混合反例。各地点只有一个恢复目标,合并统计仍能形成双峰及负对数密度的双谷;两类地点内部还包含连续的环境差异。
自选空间混合反例。各地点只有一个恢复目标,合并统计仍能形成双峰及负对数密度的双谷;两类地点内部还包含连续的环境差异。

这里的两类地点可以用来代表尚未控制的土壤、干季、坡位或人为影响差异。本例没有证明所有真实双峰都是这种混合,只是指出从双峰到同地点双稳,中间还需要证据。

Wuyts等在2017年对亚马逊地区进行空间分析,发现同时考虑季节性、土壤和距人为改造区的距离后,原本很宽的树冠双峰范围大幅缩小;靠近农业等区域的边缘效应可以解释相当一部分分布。5

但它也没有证明所有局部历史依赖都不存在。所谓较少受人类影响的区域,是按距农业或城市的空间标准划分的,不是被直接验证从未受过人类作用;“人类影响前森林”也来自重建产品。未纳入的环境差异、自然扰动与水文反馈,仍可能影响剩余重叠。

最新印度分析提供了什么,又不能直接提供什么

Dubey等发表于2025年底、列入2026卷期的研究,分析印度树冠分布与平均降雨、年际变动和季节性的关系,并使用气候情景估计未来变化。6

它很适合说明平均雨量不足以概括环境:同样雨量区间内,不同的年际变动和季节性,会对应不同的树冠分布。未来即使平均雨量增加,也不意味着森林的水分条件处处改善。

但该研究里的“稀树草原”和“森林”首先是树冠统计分类,分别采用10—40%和大于40%的覆盖阈值。这不等于已经鉴定了所有地点的物种组成与自然生境;从森林转为较稀疏的退化植被,也未必等于形成了有完整生态功能的天然稀树草原。

该研究将空间覆盖概率变成势景观,并据此讨论稳定性。原文明确依赖“用空间代替时间”的假设:把不同地点的分布,视作同一个地点长期可能经历的状态。这个转换只有在地点确实可比、过程条件与噪声等假设适当时,才能帮助推断动力学。上一节的反例展示了未控制异质性时的一个失效方式。

同样,气候模型预测某种树冠类别的概率较高,不自动等于它受扰动后恢复得更快、能承受更大的火灾,或已经测出了转变阈值。恢复能力需要连接到具体动力学与扰动。

空间分辨率也需要看清:研究把数十公里尺度的气候网格赋给其中500米尺度的植被格。保留了细植被格,并没有凭空增加细尺度气候信息。未来部分使用的是特定排放情景下的条件推算,不是已经观察到的未来。

此外,旧气候关系对后来植被变化解释不了的部分,可以提示尚有别的驱动,却不能单独给人为干预、二氧化碳、自然恢复和分类误差逐项分账。它提供的是进一步归因的入口。

哪些观察最有助于判断一条真实边界

对真实景观,更有价值的是同时追踪边界两侧的条件与过程,而不只是在全国或全洲范围画一条树冠分布。

需要知道的包括:草本燃料能否连成片,火从哪里进入、在哪里停止,幼树是否反复失去地上部分,种子能否到达,以及边界在相近环境和扰动后向哪边移动。重复测量应足够长,才能把缓慢恢复与长期不同状态区分开。

如果声称存在历史依赖,还要检查目前的外部条件是否确实可比。长期维持不同燃烧制度造成的差异,与一次短暂扰动之后、恢复相同外部条件仍保持不同结果,提供的证据不同。

这些也有助于解释恢复为什么会停住。有时问题是持续的火或边缘扰动,有时是种子来源不足,有时是水土条件不支持此前设想的目标状态;对应的机制和验证方向并不相同。不能先把所有开放景观当成待植树的空位。

森林与稀树草原的边界,由植被生长、火的传播、空间连接和地点条件共同维持或移动。 现有研究已把若干反馈环节具体化,但由一张分布图到真实恢复路径,仍需要把空间和时间接起来。

可回查:草地团块与边界、两种布局计算汇总、边界运动轨迹、单稳地点轨迹、空间混合分布、原文草覆盖变换例子、计算结果与限制。

来源与依据

  1. S1 研究综述与证据边界复评

    Reassessing the alternative ecosystem states proposition in the African savanna-forest domain

    Steven I. Higgins, Swarnendu Banerjee, Mara Baudena et al. · 2024

    New Phytologist243,1660–1669;2024 · 相关全文已核

    导论、Box1/2、地图/外部异质性与纵向/操控证据评价

    实际核理论定义与主要证据评价相关部分,不声称每个历史现场原件均已读。非洲AES理论合理、自然范围证据仍不充分,不升级为不存在任何火反馈。长期改变燃烧制度是外部条件改变,与同条件短扰动后保持不同状态区分;滞后和未控制协变量不能自动当迟滞。本文不继承Box1关于最大吸引盆等宽泛表述为严格数学定义;示例自行按局部稳定与前沿定义核对。

  2. S2 原始多地点火处理数据综合

    Herbaceous vegetation responses to experimental fire in savannas and forests depend on biome and climate

    Zachary J. Gold, Adam F. A. Pellegrini, Tyler K. Refsland et al. · 2023

    Ecology Letters26,1237–1246;2023-05-10在线 · 相关全文已核

    Methods纳入/分类/数据/统计;Results草本与草覆盖、Discussion限制

    71候选中30地点(20savanna/10forest),最低3年、排除stand-replacing crown fires及完全禁食草者;跨温带亚热带热带,不当全部热带林专有数据。239样地草覆盖为12地点子集(7savanna/5forest),其他指标以响应比整合。25地点三分类方式一致,5不一致按现场资料优先并作敏感性比较;非随机同环境大样本。模型含气候/火频率/群落及交互,另有site随机效应;不重复原模型统计。原文承认未逆转火处理检迟滞、未同地点相邻森林和savanna比较,不能直接证明同外部条件AES。打印平方根增量forest .974、savanna1.4实际反变换6例;如75起点用1.4得到101.21%,原文也承认模型越界,不用来建议真实火管理。原文一些图a/b标签与图注相反,本文不重复。数据Dryad10.5061/dryad.m37pvmd71元数据及5文件清单/哈希取得;API下载均401,本轮未取得原CSV,未声称重分析原样地。

  3. S3 原始空间随机模型与几何收支分析

    Emergent structure and dynamics of tropical forest-grassland landscapes

    Bert Wuyts, Jan Sieber · 2023

    PNAS120,e2211853120;2023-10-30在线、11-07卷期;PMC全文 · 相关全文已核

    FGBA规则/快慢分解、Eq4–9几何收支、平均场比较与Methods控制/近似条件

    格网4邻接、G/F/B/A连续时间过程,含外来种子/种子库beta、近邻alpha和自发死亡gamma;火草传播快、树慢。Eq9几何周长与草团块加权周长收支依赖pg接近1、pf小和快慢分离;森林边缘继续火传播忽略为小pf近似,未外推任意火灾。模型不预先设全景观火阈值,但微观规则和参数仍是假设。平均場准确性依赖聚集与状态,不否定所有平均場用途。unstable branch由人工反馈控制提取,非现场稳定状态观察;部分解析临界孔洞条件beta=gamma=0含两个吸收态,不能当所有森林恢复阈值。近似一维闭合和稳定性指标由模型支持,尚非任意地图直接量出自然韧性。主文相关Methods已读,SI未取得/重跑,本文不声称完整FGBA或bootstrap复现。与Wuyts2017作者重叠但问题和模型不同。

  4. S4 原始空间Markov跳跃及均场模型

    Spatial modeling of forest-savanna bistability: impacts of fire dynamics and timescale separation

    Kimberly Shen, Simon Levin, Denis Patterson · 2026

    Journal of Mathematical Biology92,44;2026-03-05;2025预印本2505.01149为同研究早期版 · 相关全文已核

    空间FGBA规则/传播核、均场近似与Section4/5暂态/边界

    实际核原始模型规则和主要空间/均场比较。火/草再生/森林尺度分离,2000随机位置常用有限仿真;森林从近邻种源长入,不含2023模型那种自发外来树补充,F=0没有内生回返路径是转移规则的后果。火模型既有显式传播核,也有局部草量S形cascade可燃项,不能说所有非线性完全从无阈值微观过程产生。无林GBA空间长暂态可偏离高灰烬均场,尚未完整分析暂态寿命;相同参数两种长时轨迹不等无限时稳态证明。当前未加入季节性、完整土壤、食草者和树种生活史,作者列为后续;phi为森林扩散率,不继承结尾把其写作fire spread的文字笔误。未重跑Gillespie或原完整分岔。

  5. S5 原始空间遥感/环境与反应扩散模型比较

    Amazonian forest-savanna bistability and human impact

    Bert Wuyts, Alan R. Champneys, Joanna I. House · 2017

    Nature Communications8,15519;2017-05-30;2018-02-21更正后HTML · 相关全文已核

    Results季节/土壤/人类边缘分层、模型与Discussion;相关SI图说明

    2018勘误10.1038/ncomms16179修参考文献23指向,不是新的独立数据,已读并保存。所谓自然区域按距离农业/城市>3km操作划分,不等真正无人影响。WCMC pre-human森林是重建产品;MAR与MSI单独均遗漏异质性,土壤及边缘作用可缩小广域双峰。局部剩余重叠、水文反馈与未测环境条件仍未决,不由此否定所有局部历史依赖。空间模型的移动/Maxwell边界适用其给定结构和初态,不误称局部双稳加任何空间作用都会成为唯一全局吸引子;我们的自选标量例子仅比较既有前沿。原遥感/空间模型未重跑;SI取得并核有关模型/图说明,不全复现。

  6. S6 原始树冠遥感统计与气候情景推算

    Forest-savanna stability in India under human interventions and changing climate

    Nivedita Dubey, Tejasvi A. Chauhan, Subimal Ghosh · 2026

    Communications Earth & Environment7,58;2025-12-10在线、2026卷期 · 相关全文已核

    树冠分类、GLM/未来情景、Methods67–85势密度/空间代时间、Discussion限制

    原主文相关结果/方法实际读。savanna10–40%、forest>40%为GMM树冠类,Methods一处forest<40%笔误不继承;不是完整物种/天然biome鉴定。2001/2020观察网格类型变更与森林身份/生态功能不同。气候0.5度复制给500m植被网格不新增细气候信息,k=5交叉验证没有报告空间分块,不视每像元独立外推验证。GLM状态概率不是直接扰动恢复率/阈值;空间PDF转势需所假设的一维SDE、稳定统计和可比地点/噪声,空间代时间与遍历性不能自动从横截面成立。负logdensity反例展示混合可伪造双谷,但不证明所有该类分析无效。旧模型与2020模型差异不唯一识别人为干预或CO2等,作者也留归因待做。2071–2100 SSP585为条件情景,不当已观察未来,未来稀树样退化状态不自动=天然savanna。SI已取得但未重拟合其表;Figshare30636614v1 API403,原图数据未得,不声称重算遥感势/预测。

最后更新 2026-09-10

来源