自然·通讯：提前一年预警冰川跃动与雨林生态临界点

导语

当我们需要预测地球系统是否即将发生灾难性突变时，一个基本前提是：我们所依赖的数据必须是“干净”的，没有趋势，没有季节性。然而，现实中的冰川流速、植被指数无不带有强烈的周期性波动，传统的预处理方法稍有不慎，就可能将真正的危机信号掩埋于人为的偏差之中。有没有一种方法，可以在不“清理”数据的前提下，直接评估周期性系统的稳定性？2026年2月6日发表于Nature Communications的一项研究给出了肯定的答案。这项研究基于Floquet乘数与动态模态分解，提出了一套无需数据预处理的稳定性评估框架，并在合成模型、冰川跃动与亚马逊雨林等多个系统中验证了其有效性。这一方法不仅对数据噪声具有高度鲁棒性，还能从时空网格中提取出系统失稳的空间模式，为预测复杂系统的临界转变提供了全新的分析路径。

关键词：临界慢化（Critical Slowing Down，CSD），Floquet乘数，系统稳定性，冰川跃动，亚马逊雨林，动态模态分解（Dynamic Mode Decomposition (DMD)）

王璇丨作者

赵思怡丨审校

论文题目：Predicting instabilities in transient landforms and interconnected ecosystems 论文链接：https://www.nature.com/articles/s41467-026-68944-w 发表时间：2026年2月6日 论文来源：Nature Communications

传统方法的困境：当“噪音”成为障碍

自然系统（如亚马逊雨林、山地冰川）可能发生难以预测的状态突变，对生态与气候影响深远。科学家常借助临界慢化（CSD）理论捕捉突变预警信号，其核心是通过统计指标异常判断系统接近临界点，但该方法需依赖平稳数据，必须先去除趋势和季节性波动。

这种预处理充满不确定性，尤其在冰川、雨林等季节性波动显著且不规则的系统中，易引入虚假信号或掩盖真实危机。例如冰川速度受季节影响剧烈，传统方法需去除“平均季节变化”，但冰川季节规律本身不断变化，残留偏差会干扰判断。雨林研究中，对各区域单独去季节化、去趋势的处理，不仅计算量大，还可能割裂系统整体性。

面对这一困境，2026年2月Nature Communications的这项研究提供了新解决方案，巧妙地将Floquet乘数与动态模态分解结合，无需预处理即可处理周期性数据，还能扩展至时空网格，为系统稳定性评估提供了新工具。

新方法的核心逻辑：绕开“噪音”，直击“核心”

既然季节性本身就是系统的一个重要组成部分，为何要费力地将其剔除？能不能将季节性视为系统的一种“周期性轨道”，然后去分析系统围绕这条轨道的稳定性呢？这正是Floquet乘数的用武之地。

Floquet乘数的概念源于对周期系数线性微分方程的研究。对于一个具有已知周期（比如一年）的系统，我们可以构建一个“单值矩阵”，这个矩阵描述了系统状态在经过一个完整周期后的演化。矩阵的特征值，就是Floquet乘数。它们的模长直接告诉我们系统的稳定性：如果所有乘数的模都小于1，那么系统就是稳定的，周期轨道会像一块磁铁一样，将扰动拉回原位，一旦有一个乘数的模接近并超过1，就意味着系统开始“漂移”，即将离开当前的周期性轨道，奔向另一个新状态。

图1. Floquet乘数

关键在于，Floquet乘数的计算无需去除季节性。研究者们利用动态模态分解算法，可以从时序数据中直接估计出这个单值矩阵。动态模态分解就像一位擅长提取主要特征的专家，它能够从看似混乱的、包含噪声的数据中，捕捉到最关键的、相干性最强的动态模式。它将系统的状态向量（例如，一系列时间点的冰川速度）通过一个线性算子联系起来，而这个算子的特征值就是我们需要的Floquet乘数。

在实际操作中，研究者会面对一组特征值。其中总有一个或多个模长接近1的特征值，它们代表了系统稳定的季节性周期。而我们需要关注的，是那些模长在逐渐增大、并向1逼近的“非季节性”特征值。正是这些“不安分”的特征值，预示着系统内在恢复力的丧失和临界点的迫近。通过这种方式，方法成功地将周期信号和临界信号分离开来，避免了传统方法中因去除季节性而引入的偏差。

冰川预警：从“点”到“面”的飞跃

为了验证方法的有效性，研究者选取了两个典型冰川：一个位于美国阿拉斯加，另一个位于中亚的喀喇昆仑山脉。这两个冰川的数据展示了其表面流速在2014至2025年间的剧烈变化。在传统的CSD框架下，这些数据因其巨大的季节性振幅而难以处理。然而，通过计算Floquet乘数，研究者成功地从中提取出了清晰的预警信号。

结果显示，在冰川进入高速跃动期前至少一年，代表非周期性不稳定的特征值就开始显著上升，并越过了代表失稳的阈值“1”。这证明了该方法能够提前预报冰川跃动的发生，这对于预测冰川对海平面上升的贡献以及冰川湖溃决洪水风险具有重要的实际意义。

图2. 全冰川稳定性

论文进一步指出，仅分析冰川上单个点的数据，其预警信号的强弱和时机很大程度上依赖于观测点的选择，这无疑引入了主观性。为此，研究者将分析对象从单个点扩展到整个冰川的面。他们将冰川表面数百个点的速度场数据，构建成一个时空数据立方体，并运用相同的动态模态分解框架进行分析。这时，整个冰川作为一个相干系统，其整体的稳定性变化被揭示出来。结果显示，对全冰川的分析能提供更稳定、更可靠的预警信号。更重要的是，通过分析动态模态分解输出的空间模态，研究者能够可视化出是冰川的哪些区域正在“蠢蠢欲动”，驱动着整个系统走向不稳定。这些区域往往是局部激活的，而非整个冰川均匀加速，这为理解冰川跃动的物理机制提供了宝贵的空间信息。

亚马逊雨林：空间上的不稳定模式

长期以来，科学家们担忧，由于森林砍伐和气候变化，亚马逊雨林可能正在接近一个临界点，从茂密的雨林转变为稀树草原。之前的研究基于传统的CSD指标（如时间自相关），已经暗示雨林的恢复力在过去几十年中有所下降。但这项新研究提供了一个更宏大、更具空间洞察力的视角。

图3. 亚马逊植被稳定性

通过分析该特征值对应的空间模态，研究者能够清晰地指出，整个亚马逊地区的不稳定性主要集中在南部地区。这片区域正是过去几十年里森林砍伐最为严重的地区，人类活动的影响也最为突出。研究进一步发现，这一不稳定性模式与森林覆盖变化、人类足迹指数、火灾频率、降水以及干旱指数等多种驱动因素都存在复杂的空间关联，这表明雨林的不稳定性是多种因素在不同时空尺度上共同作用的结果。这种方法不再依赖于对成千上万个独立像素点的时间序列进行“点对点”的预警信号分析，而是直接从系统的时空演替中提取出最具决定性的不稳定性模式。

优势与局限并存：未来的应用与挑战

这项研究成功地将基于特征值的稳定性分析方法拓展到了两类最具挑战性的数据领域：周期性数据和高维时空数据。通过巧妙地引入Floquet乘数，它绕过了令人头疼的数据预处理难题，使得对冰川、冻土、季节性河流等系统的CSD研究成为可能。而通过应用动态模态分解，它将分析维度从“点”提升到了“面”，使我们能够捕捉到系统失稳的空间模式，这无疑会加深我们对许多复杂系统突变机制的理解。为我们提供了一套更强大、更灵活的分析框架。从监测冰川跃动、预警河流决口，到追踪海洋环流的变化、揭示生态系统空间自组织背后的不稳定机制，该方法有望在未来的地球系统科学研究中扮演关键角色。它向我们展示了，当我们不再费力地“清理”自然信号中的“噪音”，而是尝试理解并分析这种周期性的动态本身时，或许更能看清，这些看似平静的系统深处，正在悄然发生着怎样的变化。

地球系统科学读书会

世界气象组织《2023年全球气候状况》报告确认2023年是有观测记录以来最暖的一年。气候变化正在以高温、干旱、洪水、野火和沙尘暴等极端天气的形式吸引人们的广泛关注。世界经济论坛《2024全球风险报告》将气候变化作为首要值得关注的风险。地球作为一个多要素、非线性的开放复杂系统，要素间相互作用关系复杂，往往牵一发而动全身。在人类活动深刻影响下，我们该如何理解并有效应对正在面临的气候变化以及其带来的社会经济等一系列议题，实现人类与地球的可持续发展？为了能够深入理解人类世背景下地球系统各要素之间复杂的相互作用与演化机制，并为人类应对未来的地球系统科学重大挑战提供一套科学的认知框架，集智俱乐部联合清华大学讲席教授陈德亮、北京师范大学教授樊京芳、东莞理工学院特聘副研究员陈爱芳、南开大学副教授戴启立老师和爱尔兰都柏林大学博士生班崭共同发起，将组织大家从新的研究范式出发梳理相关文献，并深入研读其中涉及的理论与模型。

详情请见：

1.

2.

3.

4.

5.

6.

7.

8.

9.