先说清楚"一样大"指的是哪一个量
2026年8月发表的一组CESM气候模式实验给出一个反直觉的对照。三次模拟都从同一个预工业参照态出发，都把二氧化碳浓度线性升到预工业值的4.07倍，唯一的差别是用多长时间升上去：+0.5 ppm/yr用了1750个模式年，+2.5 ppm/yr用了350年，+5.0 ppm/yr用了175年。慢升那次，全球平均近地表气温升到+5.5 °C，大西洋经向翻转环流（AMOC）始终维持；两次快升的，AMOC分别在+2.02 °C和+1.99 °C开始崩溃。1
先把比较口径钉住：这里"一样大"的是二氧化碳终值，不是升温幅度。慢升那次最终暖得更多，却没有崩溃。作为参照，2015—2025年Mauna Loa观测到的二氧化碳年增速平均为+2.6 ppm/yr，落在两档快升之间。1
如果结局只由终点决定，这三次实验应该收敛到同一类状态。它们没有。所以问题不是"阈值是几度"，而是"什么条件下终点决定不了结局"。
"越过临界"其实是三件不同的事
同一次突然翻转，至少可以由三种机制产生。2012年的一篇工作把它们分开命名：控制参数越过某个值、原有吸引态失去稳定性或不再存在，叫分岔诱导（B-tipping）；固定参数下靠涨落跳出吸引域，叫噪声诱导（N-tipping）；参数变化太快、系统"未能追随连续变化的准静态吸引子"，叫速率诱导（R-tipping）。3后续工作把R-tipping的判定形式化：它由非自治的阈值与边缘态刻画，一般不能由冻结参数后的分岔来解释。4
在本篇这种"单调升温到某个终值后长期保持"的设定下，三者给出的可检验预期不同：B-tipping只认终点，R-tipping在终点固定时仍随速率改变结局，N-tipping给的是概率而不是一刀切的门槛。本篇在一个可完全掌控的模型里把三者分开，再回头看那三条气候模式曲线。
本篇的最小对象：一个五箱环流模型
CESM那类模式跑一次要几百万核时，不能用来扫参数。入口论文因此另配了一个五箱大西洋环流模型（南、热带、热带南缘、北、深层五个箱，动态求解各箱盐度、温度与温跃层深度），并把代码和数据一并发表。12本库按其发表脚本中的方程、参数与初始化文件另写了一份实现（不调用原代码），本篇全部计算在这份实现上完成。
模型里的"AMOC强度"是一个诊断量：
北方箱比热带南缘箱重，环流就强；$D$是温跃层深度。实际进入输运方程的是$H(q_N)\,q_N$（$H$为阶跃函数）：$q_N$作为密度诊断量取负时，只表示北向下沉分支关闭，并不表示出现一支反向环流。北大西洋的额外淡水输入$E_A$使北方箱变淡、$\rho_{\mathrm n}$下降，环流减弱；环流减弱又减少向北的盐输运，使北方箱更淡——这就是盐输运反馈，也是这类模型出现双稳态的原因。5
先验证实现，再用它下结论
把发表稳态文件中$E_A=0.40$ Sv那一组数值当初值，本库实现算出的$q_N$为13.877746191869 Sv，与发表值13.877746191866 Sv一致到12位有效数字；把该状态代回本库方程，按年归一化的残差为$7\times10^{-15}$量级。这说明重写的是同一个模型，不是另一个模型。
再做延拓：固定副极地大气增暖$\Delta T_{\mathrm n}^{a}$，以$E_A$为参数追踪"开态"分支，并用扩展牛顿系统直接求鞍结点。本库求得的鞍结点与发表延拓结果在四个可对照的增暖水平上最大偏差为$1.1\times10^{-8}$ Sv。
增暖使门槛左移：$\Delta T_{\mathrm n}^{a}=0$时鞍结点在$E_A=0.4861$ Sv，$+5$ °C时降到$0.3420$ Sv，近似线性，斜率$k=0.0288$ Sv/°C。所以在固定淡水强迫0.40 Sv下持续增暖，门槛迟早降到强迫值以下——终点本身就是越界的。
顺带查出来的一件事：开态在鞍结点之前已经失稳
延拓时逐点算雅可比谱，发现一件入口论文全文未提及的事（该文把双稳区描述为由鞍结点界定）：一对复共轭特征值在鞍结点之前就越过零轴。$\Delta T_{\mathrm n}^{a}=0$时Hopf点在$E_A=0.47984$ Sv，比鞍结点0.48612 Sv低0.0063 Sv（1.3%），振荡周期约847年；各增暖水平上这个间隔都在0.006 Sv上下。这不是差分步长的假象：$10^{-6}$、$10^{-7}$、$10^{-8}$三种步长在$E_A=0.483$（两者之间）都给出$+0.0027$/yr的最大实部。从该处的精确不动点出发不加扰动可以一直停住，加$10^{-6}$的相对扰动后轨迹缓慢外扩，第2824年崩溃；在$E_A=0.470$（Hopf之下）同样的扰动4000年内始终跟随原态，最大偏离0.0009 Sv。
这不推翻入口论文的结论——间隔只有1.3%，双稳区宽度几乎不变——但它对机制讨论有用："平衡点仍在、而且局部稳定"与"任何扰动都能恢复"不是一回事。
这里要把话停在算过的地方。本库只确认了一对复特征值过零，即Hopf型线性失稳；没有算第一Lyapunov系数，没有延拓周期轨道，也没有测吸引域的几何。另一个AMOC五箱模型（Wood等，用FAMOUS标定；其箱为北大西洋、热带大西洋、印度—太平洋、南大洋与底层，与本篇的南/热带/热带南缘/北/深层不是同一套划分）里，2019年的一项工作确认了开态经亚临界Hopf失稳、吸引域在分岔前已经有界，于是"上升得足够快"本身就能把状态甩出吸引域。5那是一条与本篇主线不同的"追不上"路径，也是本篇观察到的Hopf窗口的一个候选解读——但它在那个模型里被证实，在本篇模型里只是候选。后面用对照实验看这条路径的行为特征是否出现。
让门槛移动的，是一个需要时间长出来的反馈
只有门槛移动还不够。若强迫固定、只有门槛下移，结局仍完全由终点决定。入口论文在CESM里识别出的关键是另一件事：变暖本身会减少北大西洋的淡水输入。慢升过程中蒸发增强、海冰后退（既减少绝热盖效应，也减少融冰淡水），表面浮力通量向稳定方向变化，北大西洋盐化，深层北大西洋水把盐分输出盆地，盐输运反馈被削弱。这套稳定化按百年尺度调整（受深层盐度异常输出的时间控制），而变暖造成的减弱按年到十年尺度发生——两个时间尺度不同，正是速率能决定结局的原因。1
箱模型里，这一层写成一个有效淡水强迫：
$\overline{E_A}$是外加的固定强迫，$\Delta T_{\mathrm n}$是副极地海洋箱相对初始态的升温，$\gamma$是温度—淡水耦合强度。作者由CESM慢升实验估得参照值约0.069 Sv/°C（34°S—65°N），示例中取$\gamma=0.0325$ Sv/°C。1
关键在于$\Delta T_{\mathrm n}$不会瞬间跟上大气。本库对副极地箱做阶跃响应实验（$\gamma=0$，$\overline{E_A}=0.40$ Sv，阶跃+0.25/+0.5/+1.0 °C），测得稳态增益$c=0.511$（海洋箱升温约为大气异常的一半），达到最终响应$1-1/e$的时间为60年。但这条曲线不是单指数：它先在5年内冲到最终值的60%（局地海气恢复很快），再被随之减弱的环流带回到53%附近（20年）、70%（100年），最后由盆地尺度的调整在约400年补齐。60年这个数字掩盖了这三段结构。后面会看到，快段与慢段分别管住两件事——快段决定极快强迫下还能不能自救，慢段决定慢强迫下还剩多少余量。
三条路径：同样升到 +5 °C，结局分两类
固定$\overline{E_A}=0.40$ Sv，把副极地大气异常线性升到$+5$ °C后保持，跑3000年：
- $\gamma=0$、慢升0.01 °C/yr：有效强迫一直是0.40 Sv，门槛在$\Delta T_{\mathrm n}^{a}\approx3.0$ °C（第300年）降到0.40以下，不久崩溃。这是标准的分岔诱导：终点处已经没有开态，任何速率都到不了别的结局。
- $\gamma=0.0325$、慢升0.01 °C/yr：海洋箱有时间跟上，最终$\Delta T_{\mathrm n}=3.04$ °C，有效强迫降到0.3012 Sv，低于$+5$ °C时的门槛0.3420 Sv，环流一路追住，3000年后$q_N=12.26$ Sv。原文对这一条给出的数字是$\Delta T_{\mathrm n}=3.04$ °C与0.30 Sv。1
- $\gamma=0.0325$、快升0.05 °C/yr：外加背景强迫（0.40 Sv）、耦合规则与最终大气异常（$+5$ °C）都与上一条相同，只是100年就升到位。海洋箱还没暖起来，有效强迫仍在0.39 Sv附近，而门槛已降到0.342，轨迹越过移动中的鞍结点，环流崩溃。
注意最后一条与上一条的有效强迫并不相同（0.3898对0.3012 Sv）：$E_A$依赖海洋状态，一条追住、一条崩溃，状态不同，有效强迫自然不同。相同的是外部设定与反馈规则，不是系统最终所处的位置。这也是为什么冻结$E_A$画出的阈值图只能当作参照系，不能当作完整耦合系统的稳定性判据。
真正隔离速度作用的是后两条：外部设定与反馈规则完全相同，只有升温速率不同，结局分成两类。第一条不是速度对照，而是拆掉反馈后的对照，用来说明没有这个反馈时结局由终点决定。后两条的差别就是"未能追随一个仍然存在的稳定态"。
为什么"外部强迫变快"本身不够——以及什么时候够
一个容易滑过去的问题：是不是只要外部参数变得够快，系统就会追不上？答案取决于起点离边界有多远。
关掉耦合（$\gamma=0$），把固定强迫放在门槛以内，速率从0.0003扫到5 °C/yr：
- $\overline{E_A}=0.30$ Sv（$+5$ °C时鞍结点0.3420 Sv、Hopf点0.3360 Sv，余量0.036 Sv）：最慢的0.0003 °C/yr与最快的5 °C/yr都没有崩溃；主脚本另扫的24个速率（0.002–5 °C/yr）同样0个崩溃。速率跨了四个数量级，结局完全一样。
- $\overline{E_A}=0.33$ Sv（离$+5$ °C的Hopf点只有0.006 Sv）：最慢的0.0003 °C/yr存活，5 °C/yr崩溃，二分得到临界速率0.00132 °C/yr——存在，但比耦合那条路径的临界速率小约21倍。主脚本以0.002 °C/yr起步的24个速率因此24个全崩，并不是“与速率无关”，而是全部快于这个很小的临界速率。
两个终点在冻结系统里都是稳定平衡（主特征值分别为$-0.00216$与$-0.00200$/yr），整条升温路径也都没有越过移动中的Hopf点。所以0.33那一档的崩溃既不是分岔诱导（终点仍是局部稳定的平衡点），也不是"终点守不住"（冻结主特征值为负），而是在到达终点的路上被甩了出去。最自然的读法是吸引域在Hopf点附近已经变小、升温带来的滞后足以把状态推出去；这与另一个AMOC箱模型中已经证实的亚临界Hopf图像一致5，但本篇没有直接测量吸引域，只能说它是与观察相容的解读，不是本篇证明的结论。
于是"快"要变成危险，需要两件事之一：平衡态朝吸引域边界移动（本篇$\gamma$耦合那条），或者吸引域本身已经很小（贴近Hopf那条）。离边界足够远时（0.30那一档），单调的外部变化再快也追得上——所以"变化越快越危险"不是一条无条件的规则，它要问离边界多远、边界在往哪走。
临界速率能不能算出来
R-tipping的经典判据很朴素：平衡态在移动，实际状态落后一段距离，落后量超过到边界的余量就脱离。线性近似下瞬时滞后等于调整时间尺度乘以平衡态的漂移速度，$L=M^{-1}r$。3
写到本问题上。令门槛$E_c(\Delta T_{\mathrm n}^{a})\approx E_{c0}-k\,\Delta T_{\mathrm n}^{a}$，海洋箱以时间尺度$\tau$、增益$c$滞后地跟随大气，则线性升温$\Delta T_{\mathrm n}^{a}=rt$下的余量为
这里把海洋响应理想化成了单指数，而实测曲线并不是；下一节会把这一步换成实测响应。第二项是"门槛下降"与"反馈补偿"的净较量，第三项是滞后代价，正比于$r\tau$。当$\gamma c<k$时余量在升温段单调下降，最危险的时刻是升温刚结束那一刻；令$m=0$并略去指数项，得到只用四个可测量量写成的临界速率：
代入本库测得的$m_0=0.0861$ Sv、$k=0.0288$ Sv/°C、$c=0.511$、$\tau=60$年，在$\gamma=0.0325$处给出$r_{\mathrm c}\approx0.0250$ °C/yr；箱模型实测0.0282 °C/yr，发表数组0.0282 °C/yr，同一量级。
但要看清这条式子的结构。把它展开：
本篇$m_0-k\Delta T_{\max}=-0.0580$ Sv为负，所以$r_{\mathrm c}$随$\gamma$上升（不是反比），并向上限$\Delta T_{\max}/\tau\approx0.083$ °C/yr靠拢。两个后果：一是"反馈越强越难被甩掉"这个方向被判据正确给出；二是它结构上永远给不出"找不到临界速率"，因为它有一个有限上限。
误差方向也不固定：$\gamma=0.0325$处偏低11%，而$\gamma=0.03$处给0.0202、箱模型只有0.0085，高估2.4倍——在那里它会高估可承受的速率。原因之一是它把鞍结点余量当成全部余量，而本篇已经查出完整系统先在Hopf点失稳。把它当作标定点附近的量级启发即可，不要当作安全线。
把两处理想化换掉：真实的滞后形状，以及逃逸也要时间
上面那条一行判据有两处理想化：把海洋响应当成单指数，把"越界"直接当成结局。实际上稳定态消失后状态是以有限速度离开的，门槛若很快回到它上方，系统还能被救回来。两处一起换掉：取鞍结点标准形描述逃逸，滞后部分直接用实测的阶跃响应做线性卷积，
其中$S$是实测的单位阶跃响应形状，$t_{\mathrm e}$是升温段时长。除逃逸时间常数$\tau_x$外全部参数都已测得；$\tau_x$只用$\gamma=0.0325$一个点标定，其余各点是预测。
标定后的对照（单位 °C/yr）：γ=0.0300：箱模型0.0085／最小模型0.0211／发表0.0085；γ=0.0325：箱模型0.0282／最小模型0.0285／发表0.0282；γ=0.0350：箱模型0.0461／最小模型0.0380／发表0.0463；γ=0.0400：箱模型0.0966／最小模型0.0770／发表0.0973；γ=0.0450：箱模型0.2120／最小模型在测试范围内未找到／发表0.2155。在两者都给出有限值的4个点上，最小模型与箱模型相差0.80–2.47倍；最小模型在$10^{-3}$—$10^{3}$ °C/yr的测试范围内，从$\gamma=0.0450$起找不到临界速率；箱模型（最高测到3 °C/yr）的同一转变发生在0.045与0.05之间；发表数组在$\gamma\approx0.0475$以上取其扫描上限值100 °C/yr，同样表示范围内未找到。三者都只说明"测试范围内未见翻转"，不是数学意义上的对所有速率成立——但这个转变本身被复现了，位置早了一档。标定只用$\gamma=0.0325$一个点定出$\tau_x=0.47$，其余各点是预测。
顺带回答强耦合下测试范围内为什么连越界都不发生：把$\gamma$提到0.06、升温速率提到1 °C/yr（5年升完+5 °C），有效强迫与移动门槛之间的最小间距仍有0.025 Sv，从未越界，末态$q_N=14.6$ Sv。原因就在上面那条阶跃曲线的快段：5年内海洋箱已完成60%的响应，0.06 Sv/°C乘上这一段就足以把有效强迫压到门槛之下。
有一个下限：耦合太弱时，终点本身就守不住
比较发表的临界速率数组时，本库先得到一处不一致：发表数组在$\gamma\le0.0275$处为缺测（按其脚本，缺测表示该耦合下所有测试速率都崩溃），而本库最初的扫描在$\gamma=0.0275$、升温速率0.001 °C/yr下"没有崩溃"。回查发现是本库自己的问题：升温段结束后只跑了2000年。把同一条件延长，$q_N$在斜坡后第2000年仍是11.9 Sv，到第8000年已崩溃。
顺着这条线做了一个本来没计划的计算：把耦合写进方程后直接求闭环终态（$\Delta T^{a}=+5$ °C、$\overline{E_A}=0.40$ Sv、$E_A=\overline{E_A}-\gamma\Delta T_{\mathrm n}$）并算其特征值。结果是$\gamma\le0.0275$时闭环终态虽然作为定态存在，却是振荡不稳定的（$\gamma=0.0275$时主特征值实部$+0.00031$/yr，周期约462年）；$\gamma\ge0.03$时主特征值转为负实数（$-0.0021$/yr），终态可以长期保持。
所以"临界速率"只在终态守得住时才有意义：耦合更弱时不是"升得太快"，而是"这个终点根本守不住"。所测各点中$\gamma\le0.0275$不稳定、$\gamma\ge0.028$稳定（$\gamma=0.028$的主特征值实部为$-0.00060$/yr），分界落在这两点之间，本篇没有细定位。还要说明：线性不稳定只表明这个定态不能被稳健追随，并不排除存在别的非定常开态吸引子；本篇是结合长时间轨迹判定它确实走向崩溃的，没有做全吸引子搜索。发表的图6d把这一片标成灰区（所有测试速率都失败），本库的闭环计算给出了灰区的一个具体原因。这也说明为什么一次数值实验的"没崩溃"必须附上运行时长——本篇后续的临界速率判定统一在升温段后再跑4000年，并做了尾段长度敏感性检查（尾段由2000年延长到12000年后，四个测试点中只有「γ=0.0275、速率0.001 °C/yr」的判定改变，由存活变为崩溃，与该耦合下终点本身不稳定一致；γ=0.0325的两个点判定不变）。
三条机制留下的不同痕迹
终点 vs 速率。 离边界足够远时（$\gamma=0$、$\overline{E_A}=0.30$ Sv）结局只随终点变，四个数量级的速率同命；$\gamma=0.0325$时终点固定而结局随速率变。这是B与R的分界。但要注意关掉耦合本身并不保证与速率无关——上一节紧贴Hopf的0.33那一档就有临界速率，只是小得多。
一刀切 vs 概率。 给有效强迫加上标准差0.01 Sv的白噪声：无噪声时崩溃比例是一条竖线（0.006:0%、0.010:0%、0.016:0%、0.022:0%、0.030:100%、0.042:100%、0.060:100%、0.090:100%），加噪声后变成0.006:0%、0.010:5%、0.016:2%、0.022:5%、0.030:85%、0.042:100%、0.060:100%、0.090:100%（每档40个成员）。噪声没有把方向倒过来，它把确定性的门槛抹成一段概率过渡——这正是N-tipping与R-tipping叠加时应有的样子，也是为什么单看一次崩溃很难反推机制。
升得快危险 vs 撤得慢危险。 把强迫做成一次超调脉冲：20年内从0.40升到0.50 Sv（越过0.4861的鞍结点），再以不同速度回到0.40。峰值固定、上升段固定，只改回撤时间：回撤40年以内全部存活，60年以上全部崩溃。峰值只比鞍结点高0.0139 Sv，越过门槛后的停留时间才是关键。所以"速率"在两处的作用方向相反：上升段快是危险的（追不上），回撤段快是安全的（停留短）。只说"变化越快越危险"会说反其中一半。
比较时刻。 一项独立的模式研究提醒了第四种混淆。在另一个模式中比较不同CO₂上升速率时，更快的速率造成更严重的暂态减弱，且该效应在0.5%/yr以上饱和；但更早的研究曾报告"慢升在CO₂倍增那一刻反而更弱"——同一批实验，量在不同时刻读，速率依赖的符号可以翻转。6所以"速率决定结局"这句话必须绑定在什么时刻、用什么变量判定。
条件改变时，这套解释还剩多少
模型层面。 本篇的全部计算在一个五箱概念模型里完成。它能重现入口论文的稳态与分岔位置，但只模仿了CESM中被识别出的一条速率相关机制（表面淡水输入随变暖减少），不含海冰绝热效应，也不含深对流位置变化——CESM慢升实验里第800模式年前后拉布拉多海开启深对流、AMOC增强约3 Sv那种事，箱模型里没有对应物。1另一个箱划分不同的AMOC五箱模型（Wood等，FAMOUS标定）里，失稳则由亚临界Hopf主导。5箱模型能证明"这类机制足以产生速率依赖"，不能证明真实海洋按哪一条走。
气候模式层面。 那三条CESM曲线来自单一模式、单一背景hosing（0.45 Sv）。该模式的AMOC崩溃需要约0.5 Sv的额外淡水输入，约为当前格陵兰融水速率的65倍；作者明确把背景hosing当作对已知模式偏差的订正而非物理输入，参照态的AMOC约13 Sv，也弱于观测的约17 Sv。1作者对辐射强迫方向刻意避免"临界点"一词，因为是否存在辐射强迫下的鞍结点仍未确定，他们只称"崩溃"。
观测层面。 26°N的RAPID直接观测自2004年才开始，序列长度不足以判定分岔距离；基于再分析的淡水输运指标$F_{\mathrm{ovS}}$趋势为负、指向不稳定方向，但同一批作者写明现有再分析不足以充分估计到倾覆的距离，而且该指标只在接近平衡时可解释——快速强迫下正好不适用。7本篇不给出任何崩溃年份，也不认为箱模型的临界速率可以换算成真实世界的安全升温速率。
这套想法有多老。 1997年就有简化耦合模式报告过：同样升到750 ppmv，100年内升上去导致环流永久关闭，升得更慢则只是减速；当时给出的原因也是表层浮力向深层的传递需要时间。8（本库未取得该文全文，只按摘要引用。）二十九年后的新处，是用一个完整耦合模式做到同样的对照，并把稳定化机制拆到收支层面。
还没有解决的
- 真实系统的$\gamma$与$\tau$是多少。 入口论文由CESM慢升实验估出约0.069 Sv/°C的参照值，但这是一个模式、一种情景下的诊断。临界速率对这两个量敏感：按上面的展开式，$m_0-k\Delta T_{\max}$为负时$r_{\mathrm c}$随$\gamma c$上升并趋于$\Delta T_{\max}/\tau$，$\tau$越长两项都越小。而$\gamma$与$\tau$都还没有独立的观测约束。
- 辐射强迫下到底有没有鞍结点。 淡水强迫方向的双稳态证据较强，辐射强迫方向仍未确定；如果没有，"越过移动门槛"这个几何叙述就要换一种说法。
- 两条"追不上"路径在真实系统里谁更相关。 稳定化反馈来不及长出来（本篇主线）与吸引域在Hopf点附近变小（本篇只观察到线性失稳，未测吸引域）会给出不同的早警信号：前者应表现为浮力/淡水收支项的符号变化，后者应表现为百年尺度振荡的阻尼下降。这两类信号能否在同一批观测里分开，本篇没有回答。
- 速率阈值能不能像温度阈值那样被报告。 入口论文建议把变暖阈值改写成变暖速率阈值，同时承认临界速率依赖背景气候态与稳定化机制强度。一个依赖状态的阈值还能不能当作政策口径的"安全线"，是这条路线要面对的问题。