多舛之海：哥伦比亚太平洋托图加斯湾仔稚鱼群落的 environmental drivers 研究

《Life》：Resilient Seas: Environmental Drivers of Larval Fish Assemblages in the Gulf of Tortugas, Colombian Pacific

　　

摘要
育幼栖息地通过支持鱼类早期生活阶段并缓冲环境变异性来维持鱼类种群，然而在热带东太平洋，这些栖息地在气候振荡下的功能仍知之甚少。研究人员分析了2010–2020年哥伦比亚托图加斯湾的浮游鱼卵和仔稚鱼群落，以评估其分类组成、时空变异性及环境驱动因子。利用丰富度估计量、多样性指数、卡斯韦尔中性模型（Caswell's neutral model）及多变量分析方法（PERMANOVA、PERMDISP、CAP、BIOENV），评价了ENSO循环和局地海洋学条件对育幼功能和恢复力的决定作用。群落以秘鲁鳀（Cetengraulis mysticetus）为主，石首鱼科（Sciaenidae）、仿石鲈科（Haemulidae）、银鲈科（Gerreidae）及鲹科（Carangidae）等具有小型渔业重要性的科亦有持续贡献。ENSO相位解释了年际变异性：厄尔尼诺年降低了多样性并削弱了滞留作用，而拉尼娜年则支持更为均衡的群落结构。空间异质性揭示了中心滞留区和周边多样性热点区。环境梯度，特别是叶绿素a、浮游动物生物量、透明度和pH，持续地塑造了仔稚鱼分布，反映了生产力、水质量和环流之间的耦合关系。这些发现表明，托图加斯湾发挥着双重育幼系统和生态韧性节点的功能，强调了在变化海洋中保护此类栖息地对基于生态系统的渔业管理和生物多样性保护的重要性。

**论文解读：哥伦比亚太平洋托图加斯湾仔稚鱼群落的育幼功能与生态韧性研究**

**1. 研究背景与科学问题**

海洋鱼类的早期生活史阶段——鱼卵与仔稚鱼（ichthyoplankton）——是海洋生物多样性的关键组成部分，其不仅贡献生物量、介导营养互作，还维系着生态系统的功能运转。仔稚鱼研究对于识别产卵场、估算成体种群规模及预测补充成功率具有重要的渔业科学价值。在全球变化背景下，理解生物多样性与生态系统功能如何响应气候变异和人为压力，使仔稚鱼研究日益受到关注。沿岸仔稚鱼群落与远洋群落存在显著差异，其结构受温度、盐度、海流、食物可得性及成体繁殖策略等复杂物理、化学和生物因素共同调控。在更宏观的尺度上，气候振荡，特别是厄尔尼诺与拉尼娜事件，通过改变生产力、环流及补充动态，对仔稚鱼群落产生深远影响。然而，关于全球过程与局地过程在塑造仔稚鱼群落中的相对重要性，学界仍存在争论。

育幼栖息地作为关键的生态节点，为海洋鱼类早期生活阶段提供庇护所、食物资源和适宜的水文条件，从而提升存活率和补充成功。在热带东太平洋（Eastern Tropical Pacific, ETP），沿岸仔稚鱼群落常由鳀科（Engraulidae）、石首鱼科、仿石鲈科及鲹科主导，反映了ENSO变率、上升流 regime 及河口-海洋相互作用的综合影响。尽管哥伦比亚的仔稚鱼研究可追溯至20世纪80年代，但长期连续性研究十分有限，主要原因在于时间跨度短、空间范围局限于单一河口或海湾、采样频率不规则及方法不一致。托图加斯湾（Gulf of Tortugas）地处哥伦比亚太平洋沿岸，受海洋 inflow、河流 discharge 及潮汐动力共同影响，其水文复杂性预示着潜在的育幼功能，但该区域缺乏长期的仔稚鱼群落评估。因此，理解局地生物物理驱动因子如何与大尺度气候振荡相互作用，对于评估气候变率下的育幼能力与韧性至关重要。本研究正是基于这一空白，通过十年监测数据，系统分析托图加斯湾仔稚鱼群落，检验三个核心假说：其一，海湾具有双重育幼功能，中心区域为滞留区（由少数物种主导），周边区域为多样性热点区；其二，局地生物物理因子调控仔稚鱼分布与群落结构；其三，ENSO变率通过海洋 Ni?o 指数（Oceanic Ni?o Index, ONI）相位驱动多样性、丰富度及分类组成的年际变化。

**2. 研究材料与方法概述**

研究人员于2010–2020年在哥伦比亚太平洋托图加斯湾（3°20′–4°00′ N；77°00′–77°40′ W）开展了标准化仔稚鱼采样。采样网格在2010–2014年由12个等距站位组成，2015年扩展至15个站位，涵盖沿岸、河口及海洋过渡区。浮游动物采集采用双口浮游网斜拖，网口直径30 cm，网目300 μm与500 μm，两网目样品合并分析。样品经10%缓冲福尔马林固定，仔稚鱼依据哥伦比亚太平洋浮游鱼卵与仔稚鱼图谱及Moser图谱进行种类鉴定。现场使用多参数水质仪及CTD剖面仪测定温度、盐度、溶解氧和pH，透明度使用赛氏盘测定，叶绿素a浓度采用Lorenzen法分光光度测定，总中型浮游动物生物量作为食物可得性指标。表层环流（1 m）利用GPS漂流浮标评估，垂直稳定性以Brunt–V?is?l?频率（N²）量化，水平层化以势能异常（Φ）表征。数据分析涵盖非参数丰富度估计量（Chao1、Jackknife1、Bootstrap）、香农多样性（H′）、Pielou均匀度（J′）及Caswell中性模型统计量（V），多变量分析运用PERMANOVA、PERMDISP、CAP和BIOENV，ENSO影响通过ONI协变量纳入分析。

**3. 核心研究结果**

**3.1 分类组成与丰富度评估**
研究共采集7379尾仔稚鱼，鉴定出133个形态种，隶属于100属、34科、20目。秘鲁鳀占总数量的48%，石首鱼科（18%）、虾虎鱼科（13%）和仿石鲈科（3.3%）亦有显著贡献。丰富度估计显示观测丰富度接近预测值（Chao1 = 133种，Jackknife1 = 166种的80%，Bootstrap = 149种的89%）。98%的仔稚鱼处于前期弯曲阶段，表明对局地环流和滞留过程的依赖。Caswell中性模型偏离中性预测，与繁殖聚集和滞留趋势一致。

**3.2 表层环流格局**
GPS漂流浮标揭示了托图加斯湾表层环流的空间异质性：沿岸边缘为 inflow 区，中心区域为滞留 cell，海洋过渡区延伸出扩散 pathway。滞留区内浮标位移减小，趋向于仔稚鱼聚集；扩散 pathway 则连接生物量输出与多样性混合。这种环流模式的空间镶嵌表明海湾可能运行双重育幼系统。

**3.3 局地生物物理驱动因子与ENSO变率**
Spearman相关分析（n = 114）证实仔稚鱼密度与叶绿素a（r_s = 0.62, p < 0.01）和浮游动物生物量（r_s = 0.58, p < 0.01）显著正相关；多样性指数在透明度降低（r_s = ?0.47, p < 0.05）和pH下降（r_s = ?0.41, p < 0.05）时下降。BIOENV确定叶绿素a和透明度为群落组成最佳预测因子（r_s = 0.62），CAP排序强调生产力和水清澈度梯度为主要结构轴。垂直层化方面，El Ni?o年（2015–2016）上层20–40 m N²值升高，表明混合减弱、温跃层加深、营养通量减少；La Ni?a年（2010–2011）N²值较低，反映层化较弱、垂直混合增强。水平层化方面，El Ni?o年 inflow 减弱、滞留效率降低、水文梯度平缓；La Ni?a年 inflow 较强、梯度清晰，促进了滞留与扩散的共存。

**3.4 双重育幼功能与整合模式**
PERMDISP检测到不同航次间群落离散度显著差异（F = 4.78, p(perm) = 0.001），PERMANOVA显示航次、站位和区域层面均存在显著年际变化（p < 0.01）。中心滞留区以秘鲁鳀的繁殖聚集为特征，丰富度与均匀度较低；周边扩散区支持包括鲹科、虾虎鱼科、石首鱼科和鳀科在内的更广泛科类，作为多样性热点区。H′和Hill数（N1, N2）在周边站位较高，而Margalef丰富度（d）和Pielou均匀度（J′）在中心区域较低。

**3.5 整合结果总结**
十年监测显示托图加斯湾仔稚鱼群落在分类组成上以秘鲁鳀为绝对优势，同时石首鱼科、仿石鲈科、银鲈科和鲹科的持续性贡献指示了河口-陆架连通性。拉尼娜后年份往往呈现较高丰富度和较均衡群落，厄尔尼诺阶段则对应多样性降低和强优势度。局地因子（透明度、叶绿素a、浮游动物生物量和pH）与生产力、水质量和环流的耦合效应共同塑造了仔稚鱼分布。生态学上，中心滞留区维持鳀鱼生物量，周边扩散区充当多样性庇护所，二者共同构成生物量 reservoir 与生物多样性 refugia 的双重功能，使海湾成为哥伦比亚热带东太平洋中的韧性热点。

**4. 讨论与结论**

讨论部分指出，托图加斯湾的育幼功能并非由单一因素控制，而是生物群落、物理环流和气候振荡相互作用的结果。中心滞留区通过 estuarine–shelf 连通性集中关键渔业种生物量，周边扩散区则通过维持分类广度增强系统韧性。El Ni?o期间秘鲁鳀占比超过60%，与 inflow 减弱、滞留效率降低和生产力下降相伴随，表明气候振荡可通过将生物量集中于少数优势种而破坏育幼功能稳定性；La Ni?a期间丰富度较El Ni?o年增加25–30%，H′几乎翻倍，说明周边扩散区发挥了缓冲作用，维持了生态异质性。Caswell中性模型统计量在La Ni?a期间趋近中性，暗示更均衡的补充过程。垂直和水平层化的交互作用在不同气候阶段产生了不同的机制：El Ni?o期间高N²值与深温跃层、低营养通量相关，而平缓的Φ梯度限制了扩散；La Ni?a期间低N²值和清晰的水平梯度促进了生产力和多样性的提升。在管理层面，中心滞留区作为生物量 reservoir，其易受El Ni?o抑制的特性要求采取气候响应式管理策略；周边扩散区作为生物多样性 refugia，其缓冲作用表明应保护栖息地 mosaic 而非仅关注生物量丰富区域。将环流动态、局地生物物理驱动因子和ENSO变率纳入基于生态系统的渔业管理，能够制定兼顾生物量与生物多样性保护的空间和时间显式策略。

研究结论明确：托图加斯湾的育幼功能由分类组成、环流动态、局地生物物理驱动因子和ENSO变率共同塑造。中心滞留区作为少数物种主导的生物量 reservoir，支撑补充潜力；周边扩散区作为高分类丰富度的生物多样性 refugia，促进多样性输出和生态连通性。二者的互补作用平衡了聚集与扩散过程，通过缓冲环境变率增强了系统韧性。El Ni?o通常关联多样性降低，La Ni?a则带来部分恢复，体现了系统对气候振荡的脆弱性与适应能力。尽管采样月份和站位数量存在差异，通过按局地汇总和标准化指标最小化了这些影响。保护这两种栖息地类型对于维持渔业生产力和生物多样性至关重要，将局地生物物理调节与全球气候影响纳入生态系统管理框架，能够制定预期变率并强化韧性的策略。在“迈向韧性海洋”的更广阔背景下，托图加斯湾作为保护育幼栖息地以增强海洋生物多样性和生态功能的范例，为变化海洋中的沿岸生态系统提供了韧性发展路径。