海洋热浪压倒海洋保护区的缓冲效益：地中海柳珊瑚二十年的崩溃

《Global Change Biology》：Marine Heatwaves Overwhelm the Buffering Benefits of Marine Protected Areas: Two Decades of Collapse in a Mediterranean Gorgonian

　　

海洋热浪（MHW）正在全球范围内驱动沿海基础物种的大规模死亡，威胁其持久性及其所支持的生态系统。然而，这些影响在区域尺度上的长期种群证据仍然稀缺。因此，当前以海洋保护区（MPA）为主保护行动是否能够缓解MHW的影响仍是一个未解问题。研究人员通过分析西北地中海（38°–43°N，0°–8°E）49个浅水（<40 m深度）基础柳珊瑚Paramuricea clavata种群超过二十年的种群监测数据，并结合现场温度记录，填补了这一空白。利用层次贝叶斯模型，研究人员估算了区域尺度趋势，并评估了MHW暴露和保护状态的联合影响。研究人员发现，日益频繁的极端热事件已导致广泛的大规模死亡，过去二十年间区域尺度生物量中位数下降45%。在这一广泛衰退的背景下，保护通过维持更高的总体生物量暂时缓冲了MHW影响，尽管受保护种群经历了更大的绝对损失。综合来看，这些发现表明，仅靠局部保护无法抵消气候驱动的死亡率，并强调迫切需要将保护措施与全球气候行动相结合，以保护浅海生态系统。

**研究背景**
海洋热浪（Marine Heatwave, MHW）已成为全球沿海基础物种面临的最紧迫气候驱动威胁之一。珊瑚礁、海草床和海藻林等标志性海洋生态系统因物种热阈值的突破而急剧退化，基础物种的丧失会引发级联生态系统转变，甚至导致生态系统崩溃。然而，长期、大尺度的基础物种种群动态研究仍十分有限，阻碍了对物种持久性机制的理解。国际承诺到2030年保护30%海洋（CBD 2022）部分源于海洋保护区（Marine Protected Area, MPA）在气候适应和减缓中的作用，但MPA增强极端气候事件韧性的能力受到质疑。现有MPA有效性评估多集中于捕捞物种，而忽视生境形成物种。在气候韧性方面，基础物种的经验证据稀少，且通常涉及间接机制。相反，MHW后的大规模死亡事件（Mass Mortality Event, MME）在MPA内外均有报道，表明当前保护措施可能不足以缓解极端热异常的急性影响。这一问题在地中海尤为紧迫，该区域是全球变暖最快的区域之一，极端热事件已触发跨数千公里、深达45 m的盆地级MME，影响至少90种大型底栖生物。其中，红柳珊瑚（Paramuricea clavata）是受影响最严重的类群之一，该基础物种形成密集的森林状生境，其丧失对珊瑚藻群落（coralligenous assemblages）的生物多样性和生态系统功能产生级联影响。尽管热暴露与P. clavata死亡率之间的关联已明确，但长期、空间广泛的种群数据集仍然稀缺，尤其考虑到其缓慢的生活史（长寿、慢生长、晚熟）。随着MHW频率和复发增加，短期评估无法捕捉累积、延迟或非线性种群动态。因此，反复热暴露如何塑造长期种群轨迹，以及保护能否调节这些轨迹，仍缺乏经验证据。

**研究内容与结论**
为填补这些空白，研究人员利用涵盖20多年、49个浅水种群（<40 m深度）的独特数据集，分析了西北地中海P. clavata的长期种群趋势。通过五个种群描述符（受伤表面、高度、补充、密度、生物量），研究人员量化了盆地尺度趋势，并利用T-MEDNet网络的现场温度记录，聚焦生物量以关联重复MHW暴露。最后，基于MPA-guide框架（Grorud-Colvert et al. 2021）探索不同保护水平能否缓解P. clavata在观测MHW体制下的衰退。研究发现，日益频繁的极端热事件已导致广泛大规模死亡，过去二十年间区域尺度生物量中位数下降45%。在这一广泛衰退背景下，高保护水平暂时缓冲了MHW影响，维持了更高总体生物量，但受保护种群经历了更大的绝对损失。反事实预测表明，广泛的高保护可能使生物量中位数增加14%，部分抵消区域衰退。这些发现表明，仅靠局部保护无法抵消气候驱动的死亡率，迫切需要将保护措施与全球气候行动相结合。该研究发表在《Global Change Biology》。

**主要技术方法**
研究人员利用层次贝叶斯模型（Hierarchical Bayesian Models）框架，包括层次广义加性模型（Hierarchical Generalized Additive Model, HGAM）估计区域共有趋势，以及层次线性模型（Hierarchical Linear Model, HLM）评估MHW暴露和保护对生物量的影响。MHW检测基于Hobday et al. (2016)定义，使用“heatwaveR”包，并采用固定基线（记录前十年）及阈值23°C（基于热耐受实验和现场观测）。现场温度数据来自T-MEDNet网络（Bensoussan et al. 2019），覆盖43个种群（至少10年数据）。保护状态按MPA-guide框架分为“高保护”（完全保护+高度保护）和“轻保护”（轻度保护）。样本来源：西北地中海9个研究区域，49个种群，监测时间跨度1999–2023年，共282个采样时间点，记录51,175个菌落。

**研究结果**
**3.1 红柳珊瑚在西北地中海的共有种群趋势**
通过层次贝叶斯模型分析，受伤表面、密度和生物量呈现明显长期方向性趋势。受伤表面从1999年的5.8%（95% CrI: 2.9%–10.7%）增至2023年的27.5%（22.7%–32.9%），呈现五倍增加。密度从1999年的41.4 colonies/m²（34.3–49.9）持续下降至2019年的18.6（16.8–20.5），随后略有回升至23.4（20.9–26.3）。生物量从1999年的356.3 g/m²（280.5–473.9）降至2023年的196.5 g/m²（165.2–234.3），中位数减少45%（66% CrI: 34%–54%, 95% CrI: 20%–63%）。高度和补充呈温和变化，但高度在受伤高峰期间下降，补充在死亡高峰后增加。

**3.2 生物量趋势、MHW暴露和保护效应**
基于层次线性模型，43个种群中38个种群下降概率>0.5，28个>0.75，18个>0.90。高保护种群下降概率平均0.85，轻保护种群0.64。MHW暴露（累积天数，0–164天）与生物量负相关：暴露越高，下降速率越快，中等暴露（下四分位数）即导致下降（?0.17, 95% CrI: ?0.37至0.05）。时间×暴露交互作用为负（?0.33, ?0.66至?0.01），表明MHW暴露加速生物量下降。保护状态调节该效应：高保护下，低暴露与高暴露的对比更强，但参数不确定性较大（95% CrI跨零）。高保护主效应为正（0.34, ?0.36至1.03），83%后验分布>0；时间×暴露×保护三向交互为负（?0.27, ?0.77至0.23），86%后验分布<0。反事实分析显示，若所有种群均为高保护，生物量中位数增加14%（66% CrI: ?0.8%至31%），82%后验分布为正。

**讨论与结论**
讨论部分总结了三个关键点。**4.1 重复大规模死亡事件、气候避难所的侵蚀与恢复限制**：两个受伤高峰（1999/2003年和2016年后）对应极端MHW事件，早期未受影响区域后来也出现严重冲击，表明浅水气候避难所正被侵蚀。死亡高峰之间受伤水平恢复（通过死枝断裂或菌落脱落），但生物量持续下降，因为小菌落对总生物量贡献小。重复MHW导致大面积老菌落丧失，破坏了缓慢生活史策略，恢复潜力极低。**4.2 保护对海洋热浪的缓冲有限**：MHW暴露是长期生物量轨迹的稳健预测因子，高暴露导致更陡下降。高保护种群在无MHW时维持更高生物量，但暴露时损失更大，形成“保护悖论”（protection paradox）。反事实分析表明高保护可部分抵消损失，但无法增强恢复，因为MHW复发频率远超种群恢复速率。**4.3 对地中海温带礁保护的意义**：缺乏未保护站点可能掩盖保护益处；与热带珊瑚礁不同，地中海八放珊瑚恢复有限，缺乏功能冗余和响应多样性，温带礁对气候变化脆弱。红柳珊瑚森林消失预计导致生物多样性下降、生态系统功能受损及生态系统服务减少。保护深水种群（>40 m）可能维持生态角色，但其分布和局部压力暴露尚不明确。

**研究结论**：本研究提供了迄今最全面的P. clavata衰退估计，揭示过去二十年间该基础物种生物量在西北地中海浅水区域下降34%–54%。在MHW日益频繁和强烈的背景下，所有监测种群均出现温度驱动的死亡率，多数为中度至重度影响，且恢复证据有限。在这一广泛衰退背景下，高保护种群维持了更高生物量，尽管暴露于MHW时遭受更大绝对损失。反事实预测表明，广泛的高保护可能使生物量中位数增加14%，部分抵消区域衰退。这些发现提供了经验证据：虽然高保护带来重要益处，但仅靠保护无法阻止种群崩溃，当气候驱动的死亡率超过恢复时间时，这种脆弱性在生长缓慢、长寿物种中被放大。因此，减缓未来气候轨迹仍是最紧迫的保护优先事项，同时需将保护措施与全球气候行动相结合，以保护浅海生态系统。