可逆转的群落状态转变：气候驱动的大加勒比海委内瑞拉卡里亚科海盆(Cariaco Basin)浮游植物群落变迁

《Journal of Geophysical Research: Biogeosciences》：Reversible Regime Change: Climate-Driven Phytoplankton Community Shifts in the Cariaco Basin, Venezuela

　　

浮游植物群落是海洋生物地球化学循环的重要组成部分，也是对气候驱动的环境变化敏感的指示生物。长期时间序列能捕捉这种变异性，使研究人员得以解析环境变化对局地群落的影响。本研究利用1996—2016年月度水柱观测数据结合气候指数及气象资料，探讨委内瑞拉卡里亚科海盆浮游植物群落的变化。为考察大尺度气候变率的影响，研究人员考虑了与大加勒比海局地天气及上升流相关的北大西洋多年代际振荡(Atlantic Multidecadal Oscillation, AMO)指数及多元厄尔尼诺南方涛动指数(Multivariate El Ni?o Southern Oscillation index, MEI v.2)。聚类分析识别出2004年上升流由强转弱所分隔的两个明显群落状态；2014年条件转向更强上升流机制，群落组成开始向2004年前状态过渡。两个时期之间属(genus)丰富度或香农多样性(Shannon diversity)无显著差异，但低上升流期具显著更高的群落均匀度(Pielou's evenness)和更低的物种周转率(species turnover)。为量化气候及自下而上(bottom-up)驱动因子对浮游植物群落变化的影响，研究人员采用梯度森林分析(gradient forest analysis)，发现AMO指数是群落变化最强的预测变量，紧随其后的是与上升流相关的现场变量——海表温度(Sea Surface Temperature, SST)和硝酸盐浓度(nitrate concentration)，以及MEI v.2。本研究强调大尺度气候振荡（如AMO）对热带海岸生态系统浮游植物群落的影响，2014年起向2004年前状态的过渡暗示该生态系统在十年尺度气候强迫下具周期性行为。

《Journal of Geophysical Research: Biogeosciences》刊载论文解读——卡里亚科海盆浮游植物群落对气候驱动上升流变化的可逆性响应

一、研究背景与立项依据

卡里亚科海盆(Cariaco Basin)位于委内瑞拉北岸南加勒比海，是世界上最大的厌氧海盆之一，其初级生产主要由信风驱动的季节性上升流(upwelling)维持，是热带海域难得的拥有1995年起连续月度观测的CARIACO(CArbon Retention In A Colored Ocean)时间序列站点。已有前期研究发现2004年前后该海盆出现上升流减弱、海表温度(Sea Surface Temperature, SST)升高及叶绿素a(Chlorophyll a, Chl a)降低的现象，但此前分析截止于2010年且未涵盖大西洋多年代际振荡(Atlantic Multidecadal Oscillation, AMO)影响，也不清楚2004年后是否发生永久性状态转换(regime shift)，抑或具可逆性(reversible)。由于热带长期生态数据稀缺，且区分人为气候变化趋势与自然低频气候循环（如AMO）对浮游植物(phytoplankton)群落的影响尚存争议，研究人员利用完整1996—2016年CARIACO显微细胞计数及环境—气候数据，旨在：(a)明确浮游植物生物量深度分布与多样性在20年尺度的变化；(b)解析群落结构变迁及其与物理—生物地球化学变量的联系；(c)量化大尺度气候振荡〔AMO与多元厄尔尼诺南方涛动指数(Multivariate El Ni?o Southern Oscillation index, MEI v.2)〕对浮游植物群落变化的解释力，检验群落状态在十年尺度气候强迫下的可逆性假说。

二、主要技术方法与数据来源

研究人员基于CARIACO时间序列站(10.50°N, 64.67°W)1996—2016年月度Niskin瓶采水样与CTD剖面数据，包括上层100 m八标准深度硝酸盐(NO₃^-)、磷酸盐(PO₄^3-)、硅酸盐(SiO₄^4-)、盐度、SST及21℃等温线深度，浮游植物经中性硼砂缓冲甲醛固定后按Uterm?hl法镜检计数并归并至属(level of genus, 共223属、482种)，细胞丰度与Chl a深度积分至0—100 m。气候与气象数据取自ERA5再分析资料〔向东分量u10取负值为信风风速(wind speed)〕及NOAA MEI v.2、NSF NCAR提供的AMO指数(10年低通滤波)。主要统计方法含：基于二元Jaccard距离的完全联动层次聚类(hierarchical clustering)划分年际群落状态；非度量多维标度(Non-metric Multidimensional Scaling, NMDS)分析月分辨率群落矩阵(Bray-Curtis距离，立方根变换)；Wilcoxon秩和检验(Benjamini-Hochberg校正)比较各时段环境—生物多样性指标；置换多元方差分析(Permutational Multivariate Analysis of Variance, PERMANOVA)量化季节与年代际聚类对群落组成的解释率；梯度森林分析(gradient forest analysis，引入AMO 2年滞后、MEI v.2 4年滞后及原位变量时滞检验)识别驱动群落转折的自下而上(bottom-up)预测因子及阈值响应，仅纳入全时序中出现>10次的89个属作为响应变量。缺失月份(<6个月/年)依分析类型分别处理，2016年因样本不足不参与年均值计算但参与月尺度分析。

三、研究结果

3.1 Interannual Changes in the Environment and Ecosystem（环境与水生态系统的年际变化）

年标准化z值异常值显示：正信风异常对应SST与21℃等温线深度负异常（强上升流标志），AMO正异常期(2004—2013)伴随弱上升流（低风速、高SST、深等温线、低NO₃^-）及Chl a与功能类群细胞计数显著下降；AMO于2002年后转正而生物量响应滞后至2004年转折。Chl a约80%集中于0—50 m，强上升流季(12—5月，primary upwelling)表层生物量最高，弱上升流期叶绿素最大层(chlorophyll maximum)加深。属丰富度(genus richness)1996—1997年最高，1998—2009年下降，2011年后回升，与Chl a变化不完全同步；香农多样性(Shannon diversity, H′)与Pielou均匀度(Pielou's evenness, J′)呈类似但较平缓的年际波动，均匀度在低上升流期(Cluster 2)整体偏高。

3.2 Clustering Analysis Yields Contrasting Ecological Regimes（聚类分析揭示对比鲜明的生态状态）

基于年聚合二元Jaccard距离的层次聚类划分出两簇：Cluster 1为1996—2003(Early Cluster 1)与2014—2016(Late Cluster 1)，Cluster 2为2004—2013。年际属周转率(species turnover)2004年最高，Cluster 2期内较低且稳定，2014年再度升高指向向Cluster 1过渡。NMDS月分辨率排序沿NMDS1轴分离两簇，AMO与降水正异常关联Cluster 2，盐度与MEI v.2正异常关联Early Cluster 1，Late Cluster 1多数点落入Early Cluster 1空间但部分仍处Cluster 2区，表明处于过渡态。密度分布与配对Wilcoxon检验显示：Late Cluster 1与Early Cluster 1间SST、盐度、甲藻(dinoflagellates)及定鞭藻(haptophytes)无显著差异但均异于Cluster 2，支持部分恢复；硅藻(diatoms)与蓝细菌(cyanobacteria)在Late Cluster 1中仍与两时期均有差异，未完全恢复至1996—2003水平；属丰富度与Shannon指数三期间无显著差异，但Cluster 2的Pielou均匀度显著高于合并Cluster 1期(p_adj<0.05)，反映弱上升流下优势类群受抑、群落分配更均匀。

3.3 AMO and Upwelling-Related Variables Drive Community Changes（AMO及上升流相关变量驱动群落变化）

梯度森林分析含2年滞后AMO与4年滞后MEI v.2显示预测力排序：AMO > NO₃^-浓度 > SST > MEI v.2 > 21℃等温线深度 > 盐度 > 其余营养盐与气象因子。累积重要性曲线示AMO负异常(负相位，南移ITCZ、强化信风及上升流)引发硅藻、甲藻、定鞭藻强组成转折；MEI v.2正异常(El Ni?o)对定鞭藻与蓝细菌具阈值效应，负异常(La Ni?a)关联甲藻变动；NO₃^-在6—8 μmol·L^-1区间对应群落最强更替，浅21℃等温线(强上升流)主要影响硅藻响应。分季节探索性运行表明AMO全年均为首要或次位预测因子，主上升流季MEI v.2预测力升至首位。

四、讨论与结论总结

讨论指出2004年上升流转弱致硅藻与甲藻细胞计数降一个数量级、蓝细菌降近30倍，但Chl a降幅较小，结合HPLC色素数据推断群落向<20 μm显微检出限以小的细胞(小型/微微型浮游植物, nano- and picophytoplankton)偏移，与升温分层抑制大型细胞的理论预期一致。2014年起物理条件回归强上升流、NO₃^-与Chl a回升、分类组成(基于属水平有无)向1996—2003相似，虽硅藻等未完全复原但支持"可逆状态转变而非不可逆永久退化"的核心论点。局限含显微法漏检pico-/nanophytoplankton(如Synechococcus实际应旺发于分层期但镜检全为0)、2014—2016月覆盖不全及未纳入顶级下行(top-down)捕食压力(如摄食类浮游动物grazing)。研究者强调AMO是大加勒比海浮游植物群落变迁最强力的大尺度气候预测指标，突显长期连续观测对辨识非线性生态响应的价值。

结论(Conclusions)译文：

研究人员发现十年尺度气候变率——尤指AMO指数——强烈调控热带海岸卡里亚科海盆浮游植物群落组成。显微细胞计数揭示两分明群落状态：1996—2003年(Early Cluster 1)伴强上升流与AMO负异常；2004年后(Cluster 2)转为弱上升流与AMO正异常。两期总体多样性(genus richness与Shannon H′)无差异，但低上升流期具更高Pielou均匀度(J′)与更低物种周转率。此转变非永久性——2014年环境回归强上升流与AMO负异常，群落开始向2004年前状态过渡。多层次证据(NMDS空间位置、SST/盐度/定鞭藻与甲藻分布无差异于Early Cluster 1而异于Cluster 2)支持系统朝前期状态复原，部分类群(硅藻、蓝细菌)截至2016年尚未完全恢复，处过渡态。梯度森林分析确定AMO为群落变化最强预测变量，其次为NO₃^-浓度、SST及MEI v.2，证实大西洋与太平洋气候振荡经上升流动力学间接驱动该热带近海生态系统。研究提供证据显示十年气候强迫下浮游植物群落具向原初状态可逆的行为，强调长期生态时间序列及AMO作为区域生态状态指示因子的意义。