观潮 · 模式预报与气象 AI 观察
CESM 40 个成员拆年际变率：中高纬下降主因是潜热输送随增暖变强，海冰退缩不是主导
> 增暖不只抬高平均温度，还会改变温度围绕平均值的年际摆幅：CMIP5/CMIP6 一致投影低纬变大、中高纬变小。这篇用一个加了天气噪声的随机湿能量平衡模型（SMEBM）去复现 CESM 大集合的纬向变率变化，并逐项拆开驱动。结论有两条反直觉的地方：一是中高纬变率下降的主因不是海冰退缩，而是时间平均增暖通过潜热输送的非线性把温度异常阻尼得更狠；二是低纬上升来自局地反馈参数变得不那么稳定，且把 ENSO 统计剔除后依然存在，所以不是 ENSO 主导。作者同时给出反向条件：如果时间平均增暖是空间均匀的，高纬变率反而会增大。
缺的是哪一块
年际温度变率怎么随温室强迫变化，投影结果本身分歧不大——CMIP5 和 CMIP6 在高排放情景（RCP8.5 或 SSP5-8.5）下都给出低纬增加、中高纬减小，其中北半球高纬的下降和陆地上低纬的上升跨模式最稳健。分歧在机制归属：
- 中高纬下降，多数研究归给海冰流失，理由是冰没了之后大气直接接触海洋的大热容（Stouffer 与 Wetherald，2007；Huntingford 等，2013；Brown 等，2017；Olonscheck 等，2021；Coquereau 等，2024）；另一条解释是极地放大增暖把经向温度梯度压平，通过热平流削弱变率（Screen，2014；Schneider 等，2015；Holmes 等，2016）。在年际尺度上，Coquereau 等（2024）认为梯度变平是次要的，Olonscheck 等（2021）和 Shi 等（2023）则认为它可能是重要贡献者。
- 低纬上升，此前被归因于陆面变干（Olonscheck 等，2021；Shi 等，2023）、平均云量带来的长波辐射增加（Shi 等，2023），或者 ENSO 本身的变化（Coquereau 等，2024）。
缺的是可开关的工具：已有机制工作几乎全在完备气候模式里做，过程互相纠缠，没法把某一项单独关掉看变率怎么变。Budyko-Sellers 型能量平衡模型虽然长于解释时间平均的强迫响应，但此前用于内部变率的是"干"版本（扩散近地面温度），湿静态能版本（MEBM）没有变率表示。
方法怎么搭的
做法是把外加天气噪声加进 MEBM，得到随机湿能量平衡模型（SMEBM，随机湿能量平衡模型），经向输送按湿静态能的 Fick 扩散参数化，同时显式区分感热与潜热输送；有效地表热容量和局地反馈参数不自由设定，而是从 CESM 拟合得到。
数据侧用 CESM 第一版大集合的 40 个成员：1920 到 2005 年历史强迫，之后转 RCP8.5 直到 2100 年。变率不用单个成员的时间序列算，而是取每一年 40 个成员之间的离散度作为该年的年际内部变率，依据是准各态历经假设（Olonscheck 与 Notz，2017）。早期时段取 2000—2025 年平均，末期时段取 2075—2100 年平均，两者相减即变率变化。
ENSO 的处理给了两套：一套是既有线性方法，一套是非线性、带滞后的机器学习方法，两套都用来把 ENSO 信号从地表温度里统计剔除，检验结论是否只是 ENSO 的影子。
隔离驱动的方式是参数替换实验：先用 21 世纪初的 CESM 参数跑出早期变率，再把参数逐项换成 21 世纪末的值重跑，比较差异。关键的一项是用乘积法则把输送辐合拆成三项——(i) 有效扩散系数作用于温度异常；(ii) 感热部分；(iii) 一个有效阻尼项，取决于本地湿静态能温度敏感度与邻近地区的相对高低。变异率对集合平均温度的依赖全部走潜热通道，感热通道不参与。
结果数字
环节 · 结果 · 出处位置
集合规模与情景 · CESM1 大集合 40 个成员；1920—2005 历史强迫，之后 RCP8.5 至 2100 · 方法节
变率口径 · 每年 40 个成员之间的标准差；早期 2000—2025，末期 2075—2100 · 方法节
早期变率型 · 极地最大，北半球大于南半球；赤道有次峰，滤波显示该峰主要由 ENSO 驱动 · 图 1a
末期变化型 · 副热带与热带低纬增加、中纬度减小；极地北半球减小、南半球增加，南大洋有局地减小 · 图 1b
低纬的相对量级 · 绝对变化很小，相对变化幅度与北半球高纬的变化同量级；最强的低纬增加落在陆地上 · 图 1c、附图 S1
ENSO 滤波 · 滤掉 ENSO 后变化型基本不变，只在热带与副热带有小影响 → 变化不是 ENSO 主导 · 结果节
SMEBM 复现度 · 大体捕捉到 CESM 的纬向变率变化型 · 摘要与结果节
中高纬驱动 · 海冰退缩抬高有效地表热容量 + 时间平均增暖调制的湿静态能辐合变化 · 摘要与结论节
海冰的地位 · 海冰流失**不是**变率变化的主导驱动，与 Coquereau 等（2024）相反 · 图 3b、讨论节
高纬阻尼机制 · 背景增暖抬高有效扩散系数，对温度异常阻尼更强；克劳修斯-克拉珀龙的非线性使该项在低纬增加更多 · 图 3c
极地通道 · 75°N 以北时间平均增暖远大于其他地区，使阻尼项由"放大"转为"削弱"，北半球高纬变率进一步下降 · 图 3c
反向条件 · 若时间平均增暖是空间均匀或仅弱极地放大，高纬变率会通过阻尼项**（iii）增加** · 附图 S4
低纬驱动 · 局地反馈参数变得不那么稳定，温度异常被阻尼得更弱 → 变率增加 · 结论节
早期极地变率峰值本身也是输送辐合造出来的：项（i）中，湿静态能对温度的敏感度越高，有效扩散系数越大，而克劳修斯-克拉珀龙的指数形式使热带该项远大于高纬，低纬异常因此被抹得更平；项（iii）在本地敏感度低于邻近时对异常起放大作用，极地敏感度低于邻近的中低纬，异常于是向极地汇集。作者提示热容和反馈参数的空间变化也参与形成这个峰值，但单靠输送机制已足以产生部分极地放大。
低纬那条结论有一个自洽性保证：反馈参数是在把 ENSO 从地表温度和大气顶通量两边都统计剔除之后估计的，因此其变化不太可能只是 ENSO 的间接影响。作者还指出，这个"低纬反馈变不稳定"与先前关于强迫响应中反馈演变的研究方向一致——那些研究把反馈演变归因于海表增暖型（或海洋吸热）改变低层稳定度、减少低云、使短波云反馈越来越正（Ceppi 与 Gregory，2017；Rose 与 Rayborn，2016）。本研究的含义是：这些此前只在时间平均强迫响应语境下讨论的反馈变化，可能同样作用于年际变率。
作者承认的限制
限制 · 原文口径
单模式 · 研究设计只解释 CESM 一个模式，识别出的过程未必推广到其他 CMIP 模式
理想化程度 · SMEBM 高度理想化，结果对不确定的拟合参数敏感，应作机制示意定性解读，不宜定量使用
反馈参数不确定 · 反馈参数是主要不确定源；回归只假设纯局地结构，相对同时考虑非局地反馈的做法会产生偏差
估计方法有偏风险 · 基于变率的反馈估计对滞后与时间尺度选择敏感（月平均对年平均）；所用零滞后回归在温度变率直接由随机大气顶通量驱动时会受回归稀释影响，因此估计可能有偏
部分缓解证据 · 21 世纪末低纬反馈变弱在强迫反馈参数中同样可见，后者不受同样的滞后与回归稀释问题影响
热平流机制未被排除 · 常数噪声强度假设按构造把热平流机制排除在 SMEBM 之外，因此本研究不否定该机制在年际尺度上的作用
结论依赖增暖型 · 结论依赖模拟出的时间平均增暖型；换成均匀增暖，高纬变率的符号会反过来
一句话判断
把年际变率的变化拆到"热容量 + 输送辐合 + 局地反馈"三项之后，中高纬的下降主要来自增暖本身让潜热输送的阻尼变强、而非海冰退缩——由于模型高度理想化、且结论依赖于极地放大增暖这一空间型，这套机制应当作定性解释来读，不宜当作定量结论。