冰岛—苏格兰溢流水（ISOW）在过去80万年间持续双稳态及其对AMOC稳定性的长期约束

《SCIENCE ADVANCES》：Persistent bistability of Nordic Seas overflow waters over 800,000 years

　　

大西洋经向翻转环流（Atlantic Meridional Overturning Circulation, AMOC）的稳定性是未来气候预估中的主要不确定性。然而，轨道时间尺度上AMOC临界行为（tipping behavior）的直接证据仍然稀少。在此，研究人员将来自北大西洋沉积岩芯的12个流速记录与数据驱动的混合模型相结合，重建了过去800,000年间冰岛—苏格兰溢流水（Iceland-Scotland Overflow Water, ISOW）的强度；ISOW是深部大西洋环流的关键下支。重建结果揭示弱（冰期）与强（间冰期）环流模式持续交替、两种状态内普遍存在千年尺度变率，以及弱—强恢复与强—弱崩塌之间的系统性不对称。阈值和相空间分析表明，状态转换发生在一个狭窄的大气CO₂和海平面窗口内，并常遵循类滞后（hysteresis-like）回路。尽管双稳态在过去的800,000年间持续存在，但其敏感性和时序随时间演化。这些发现为北欧海溢流水的双稳态及其在变暖世界中对AMOC恢复力（resilience）的意义提供了长期、高分辨率、数据驱动的约束。

论文解读：该研究关注大西洋经向翻转环流（Atlantic Meridional Overturning Circulation, AMOC）在轨道时间尺度上的稳定性。AMOC稳定性是未来气候预估的重要不确定性，但轨道尺度上AMOC临界行为的直接观测证据不足；双稳态、滞后和阈值行为主要来自概念箱式模型与耦合气候模型，且不同模型给出的阈值窗口差异较大，缺乏连续、高分辨率的深部AMOC强度重建。ISOW是北大西洋深层水（North Atlantic Deep Water, NADW）的关键下支，但过去冰期—间冰期贡献约束不足；基于²³¹Pa/²³⁰Th的记录时限较短且解释存在争议，已有ISOW重建多集中于近期。为此，研究人员从Reykjanes Ridge沿ISOW路径的两个沉积岩芯提取ISOW强度指数，构建过去380 ka的ISOWdata，并用混合模型延伸至800 ka，识别阈值、滞后和敏感性演化。论文发表在《SCIENCE ADVANCES》。其重要意义在于为北欧海溢流水双稳态提供长期、高分辨率、数据驱动的约束，并为评估气候模式中的稳定性边界提供基准。

关键方法概况（不超过250字）：研究使用Reykjanes Ridge沿ISOW路径的两个北大西洋沉积岩芯MD03-2685-84Cq（约2000 m）和MD03-2673-74Cq（约2800 m），保留AICC2012年代，选取12条流速代理记录并每100年插值，采用蒙特卡洛主成分分析（principal component analysis, PCA；10,000次）构建ISOWdata。随后以线性回归（ln(IRD)与NPS）加XGBoost机器学习（machine learning, ML）的混合模型，结合CO₂、底栖有孔虫δ¹⁸O、冰筏碎屑（Ice-Rafted Debris, IRD）、NPS及轨道参数，经100次蒙特卡洛和Optuna调参外推至800 ka，并用SHAP归因；阈值检测采用7-ka平滑、高斯混合分离器、持续交叉与CO₂-全球平均海平面（global mean sea level, GMSL）核密度；谱分析使用Pyleoclim小波和小波相干。

研究结果：ISOWdata stack over the last 380,000 years：通过两岩芯12条记录合成，过去四个冰期—间冰期旋回中ISOW在间冰期强、冰期弱，转换突然且无持久中间平台；去除7-ka滑动平均后，强弱状态内均有幅度和时序可比的千年尺度波动，说明千年尺度变率不局限于某一状态。ISOWhybrid stack over the last 800,000 years：混合模型将记录延伸至80万年，交叉验证和敏感性试验表明重建稳健；SHAP显示IRD和副极地北大西洋温度主导，其次冰体积和CO₂，轨道参数影响较弱但不可忽略。80万年间ISOW持续呈强、弱两模式，转换尖锐，弱—强恢复约2000年，比强—弱崩塌约5000年快约2.5倍；两模式均含亚轨道变率，说明千年尺度变率是系统内在特征。Forcing thresholds of ISOW regime shifts：数据驱动阈值检测将每次转换映射到CO₂-GMSL空间，发现弱—强转换约在CO₂ 242–265 ppm、GMSL ?55至?39 m，强—弱转换集中在CO₂ 201–208 ppm、GMSL ?57至?50 m；阈值呈窗口而非单值，反映多驱动交互。Asymmetric loops on glacial-interglacial timescales：80万年识别23次状态转换、11个冰期—间冰期旋回；相空间显示重复滞后环，同一CO₂和GMSL下ISOW可处两态，体现历史记忆；部分旋回为完整环，部分仅在CO₂或GMSL轴表现滞后，差异与背景气候和强迫幅度有关。

讨论部分总结：ISOW strength as a constraint on the deep North Atlantic circulation方面，三条独立证据支持ISOWdata反映盆尺度深部北大西洋环流而非纯局地：两岩芯不同水深相位振幅高度一致；沉积物源显示南北输送持续至全冰期；与多站点NADW磁学stack千年尺度一致，时序偏移≤800年。与Bermuda Rise ²³¹Pa/²³⁰Th记录互补；ISOW不是总AMOC直接量，而是北欧海溢流和深对流的稳健约束。Persistent ISOW bistability over the last 800,000 years方面，数据与混合stack显示80万年间两种模式持续交替，转换集中于受限CO₂-GMSL窗口并具滞后环，支持深大西洋环流多平衡态和路径依赖；千年变率叠加于双稳态，挑战深部不稳定仅限冰期或中间态的观点，稳定模式并非静止平衡。Mechanisms of ISOW regime shifts方面，低CO₂和中等冰量下海冰扩张、淡水增加和层化增强削弱深对流，盐平流反馈锁定弱态；较高CO₂和中等冰量下，退冰期淡水层延迟恢复，海冰退缩后盐平流反馈快速锁定强态，解释恢复快于崩塌；阈值集中于约?60至?40 m GMSL，为气候模式提供基准。Evolving sensitivities of Nordic Seas overflow waters方面，Mid-Brunhes Event（MBE）前后ISOW强度无系统变化，但100 ka频带相位关系改变：MBE前几乎同步跟踪边界强迫，MBE后滞后所有四类驱动5000–10000年；说明双稳态机制仍在，但有效敏感性和响应时间随背景气候重组。

研究结论翻译：总之，这些结果表明ISOW双稳态在中—晚更新世持续存在，但并非北欧海溢流系统的固定属性；其有效敏感性和响应时间取决于背景气候状态。展望未来，研究结果将关注点从北欧海溢流是否能发生突然、滞后转变转向其稳定性边界如何响应人为强迫。随着大气CO₂水平目前升高至更新世范围之外，一个开放问题是深对流将保持在双稳态的强模式内，还是其稳定性景观将被根本重塑，从而可能开启气候模式实验所提示的崩塌路径。