观潮 · 模式预报与气象 AI 观察
六个大集合喂出的 CNN 判北极失冰风险:2007 年前有信号,2015 年起风险重新走高
> 快速失冰事件(RILE)指九月海冰在几年之内丢掉相当大一部分。气候模式里这类事件常见,观测记录太短、且 2005 年以来九月海冰没有统计显著的下降,所以真实北极到底有多脆弱很难判。作者用六个大集合(每个至少 20 个成员,历史与 SSP3-7.0 情景合并,1980—2090 年、只取九月)训练一个卷积神经网络,做「状态依赖的脆弱性」判别,再用留一法交叉验证到未见过的模式上。留一检验里实际 RILE 数随预测值单调上升,Spearman 秩相关 0.8 到 0.95;用到卫星时代观测上,2007 年前后那次下降之前有信号,2015 年前后至今风险走高。关键限定作者写在明处:CNN 能区分高风险与低风险状态,但给不出校准过的绝对概率;SHAP 可解释性分析假设特征可独立变化,而气候场空间自相关很强,这个假设并不成立。
缺的是哪一块
北极海冰预计本世纪会稳定进入季节性无冰状态,但下降速率的不确定性很大——内部变率强、模式物理不同、未来排放路径也不确定。
RILE 值得单列出来,是因为它的时间尺度是「几年」,比「本世纪内无冰」更直接影响生态、沿岸社区和航道。但它和「临界点」不是一回事:简单模式和复杂模式里,夏季海冰似乎都不适用不可逆的突变动力学;RILE 也不必然是失控反馈的结果,而是内部变率叠加强外强迫的产物。
麻烦在于两头都缺数据:不同模式里 RILE 的发生频率差别很大,而观测记录太短,2005 年以来九月海冰没有统计显著下降。模式分歧加观测稀少,直接评估真实系统的脆弱性非常困难。
数据与判据怎么定
项 · 设置
训练数据 · 六个气候模式大集合,每个在历史与 SSP3-7.0 情景下均至少 20 个成员
时段与月份 · 两情景合并,取 1980—2090 年,只用九月
输入变量 · 月平均位温 `thetao`、海冰密集度 `siconc`、海冰厚度 `sithick`
RILE 判据 · 九月海冰范围的五年滑动平均,其五年趋势连续至少四年低于阈值(沿用 Auclair and Tremblay 2018、Sticker et al. 2025 的定义)
平滑方式 · 改为**前向**滑动平均,避免训练时未来信息泄漏
验证方式 · 留一法交叉验证:留出整个模式不参与训练
观测侧数据 · PIOMAS 海冰厚度、EUMETSAT OSI SAF 海冰范围;次表层温度用 ASTE、EN4、ORAS5 三套互相独立的数据(多数取 1981—2025 年,ASTE 为 2002—2017 年)
判据本身还有一个副作用:用前向滑动平均,事件开始年会略早于最陡下降段,结束年也略早于海冰降到底。作者认为这不影响目标——他们关心的是「对 RILE 发生的脆弱性」,不是精确预测发生时刻。
六个模式在训练样本里的表现差别明显:CESM2、EC-Earth3、HadGEM3-GC3.1-MM 平均每个模拟大约出一次 RILE,而 MPI-ESM1-2-LR 大约每四个模拟才出一次。这些发生率低于 Sticker 等人 2025 年的结果,作者解释为该研究用的是更高排放的 SSP5-8.5 情景。
结果数字
检验项 · 结果
留一法区分能力 · 实际 RILE 数随 CNN 预测值单调上升,Spearman 秩相关系数 0.8—0.95
主要前兆因子 · 持续偏薄的海冰、剩余海冰范围、持续偏高的海洋上层热含量
SHAP 的时间含义 · 长期持续偏暖的边缘海对可预测的 RILE 贡献更大,而不是最近几个九月的异常增暖
观测期结果 · 2007 年前后那次海冰骤降之前出现 RILE 信号;2015 年前后至今风险上升
跨模型一致性 · 不同 SHAP 背景下的图像一致,说明网络学到的是 RILE 的通用前兆,而不只是拟合各个模式
训练样本构成会影响结果:把 EC-Earth3 和 HadGEM3-GC3.1-MM 从训练中剔除后,CNN 倾向于低估脆弱性;CanESM5 在留出时反而校准得更好。
原文图

图1 快速失冰事件的定义示例:HadGEM3 第 12 号成员在 2023 到 2028 年间 9 月海冰密集度断崖式下降(左、中两图),右图黑点标出事件时段,红线为前向滑动平均的趋势诊断。
作者承认的限制
限制 · 原文口径
概率未校准 · 总体看,CNN 能区分高风险与低风险状态,但**无法给出绝对意义上校准良好的 RILE 概率**。作者明确不要求留出模型在概率上校准,只要求能在时间上区分相对风险
SHAP 假设不成立 · SHAP 依赖「特征可独立变化」的假设,而海冰这类气候场空间自相关很强,该假设在气候系统中并不成立;作者说这只能给出网络工作方式的粗略图景
模式数量受限 · 研究必然受限于「成员数足够多的模式」,只能取到少数几个模式,走的是「用全部可用数据找通用规律」的路线,而不是按模拟观测趋势的能力筛模式
判据参数敏感 · 作者指出 RILE 定义对趋势阈值、滑动平均窗口、「连续四年」的要求都敏感,但为与已有研究保持一致,未做变参数试验
两个模型的一致性较差 · CanESM5 和 CESM2 在两种 SHAP 分析之间结果不一致,可能因为 CNN 对这两个模式的 RILE 预测本身不够清晰
一句话判断
CNN 在这里的用处不是预报「哪一年崩」,而是给当前北极状态打一个相对脆弱性分——它能可靠区分高低风险,但输出不能当概率用;2015 年以来风险曲线重新上扬这个信号,值得拿更长序列的观测继续盯。