暖缘欧洲鳗玻璃鳗（Anguilla anguilla）入海口的体况与入湾动态：丰度与质量耦合的招募信号

《Estuarine, Coastal and Shelf Science》：Body condition and estuarine entry of glass eels ( Anguilla anguilla) at the European eel’s warm southern margin

　　

在物种温暖南缘的欧洲鳗（Anguilla anguilla）玻璃鳗入湾阶段，研究人员对个体体况与逐网次（haul-resolved）入湾密度进行了联合分析。研究于欧洲温暖边缘的一处河口采集了8,150尾玻璃鳗及251次网次样本。结果表明，体况随主要入湾脉冲达到峰值，并在季节两端回落。在冬季脉冲内部，体况自峰值密度到达期向脉冲尾部逐渐下降。涨潮、水温和盐度共同塑造了细尺度（fine-scale）的入湾密度格局。

论文解读

一、研究背景与科学问题

欧洲鳗（Anguilla anguilla）是一种兼性降海洄游、单次产卵的鱼类，其生命周期极为复杂：单一泛种群在萨加索海（Sargasso Sea）产卵，叶状幼体（leptocephalus）跨越大西洋后在大陆架变态为透明的玻璃鳗（glass eel），随后进入河口并上溯至淡水、半咸水和沿岸水域生长。河口过渡是招募（recruitment）的关键瓶颈，也是大陆种群招募指数估算与预测的核心环节。然而，该物种的招募已严重崩溃：以1960—1979年为参照期，2024年泛欧洲“其他欧洲地区”玻璃鳗指数仅为历史水平的约7%，“北海”指数约为1%。国际海洋考察理事会（ICES）将该种群评估为临界状态，并建议将所有人为死亡率降至零。欧洲鳗被列入IUCN红色名录“极危”等级及CITES附录II，并受欧盟专项恢复框架约束，但在物种南缘仍具有重要的文化与经济价值。

招募在原则上包含两个维度：到达的鳗鱼数量及其体况。然而，实际监测几乎完全依赖数量。玻璃鳗在变态期间及河口转运的大部分阶段不摄食，依赖叶状幼体阶段积累的有限能量储备。体况并非仅仅是海洋穿越的被动记录：体况较好的玻璃鳗更有可能维持上溯迁移，而能量耗竭则促使个体转向偏好咸水并提前定居。因此，入湾时的体况可能与丰度共变，同时携带超越数量的额外信息。这种耦合是否真实存在——即从同一次渔获中获得的体况描述符能否为已计数的数量增添信息——尚未在物种温暖南缘得到检验，这正是本研究的切入点。

此外，玻璃鳗入湾并非被动漂流，而是行为介导的迁移，选择性潮汐流输运（Selective Tidal Stream Transport, STST）是主导的近因机制：玻璃鳗在涨潮时上升至水柱、落潮时贴近河床，以低能量代价实现净上游位移。叠加于潮汐引擎之上，水温、盐度、淡水流量、降雨、浊度和光照共同调节玻璃鳗的活动性、垂直位置和可捕性，使入湾呈现强烈季节性并组织为离散脉冲。瓜达尔基维尔河口（Guadalquivir estuary，西班牙西南部）是欧洲鳗分布范围最南端的大型河口，是温暖边缘的前哨站点，但其已不再纳入协调监测网络——该站点向ICES/WGEEL招募指数提供科学玻璃鳗序列仅至2007年，目前物种大西洋南缘已无任何站点参与到达体长追踪。

二、研究设计与核心问题

研究人员提出了四项预测。H1：在单个招募季节内，平均体况呈非单调变化，在入湾峰值处达到最高，并向稀疏的秋季前兆和晚冬尾部两侧回落。H2：跨采样战役（campaign），平均体况与入湾密度正相关。H3：在脉冲内部，平均体况自起始向尾部下降，且该下降由给定体长下的质量承载而非由体长承载。H4：入湾密度在涨潮时高于落潮；在持续高浊度下昼夜主效应为零；潮汐相位与昼夜周期存在交互作用。H4在拟合前已明确设定，属确证性分析；H1—H3为探索性分析。研究人员将丰度与体况视为同一招募信号的两个侧面而非两个独立研究。

三、主要关键技术方法

研究在瓜达尔基维尔河口下游主航道固定站点（距河口约8 km）开展，采用当地传统persiana渔具，以被动逆流方式过滤水柱。采样组织为每月新月（new moon）时段的战役，每战役含72小时连续采样（10—4月）或24小时（5—9月）；每个24小时周期产出四次网次——涨潮与落潮、昼与夜——各限定在平潮后最初两小时。入湾密度以每10⁵ m³过滤水体的个体数表示。生物测量方面，由单一训练有素的观察者对2013/14和2014/15两个季节的19个战役中阶段VB玻璃鳗进行体长（TL, mm）和湿重（精确至0.1 mg）测量，共8,150尾。体况以残差体况指数（Kr）量化，即对所有阶段VB个体合并拟合的log₁₀质量对log₁₀体长普通最小二乘回归的残差。入湾密度建模采用线性混合模型，随机截距为采样日嵌套于新月战役（1 | campaign/day），固定效应包括五个保留的z标准化连续预测变量和三个分类变量（潮汐相位、昼夜周期、风力等级）。体况—密度耦合在个体层面建模，脉冲内动态以事件（event）而非个体为复制单位进行回归分析。分析代码已完整存档于Zenodo仓库（https://doi.org/10.5281/zenodo.22766460）。

四、研究结果

4.1 体长、异速生长与体况指数

在9,201尾经色素分期的玻璃鳗中，8,150尾（88.6%）为阶段VB。阶段VB个体平均体长64.6 mm（SD 3.8），平均湿重0.248 g（SD 0.051）。质量—体长关系为log₁₀M = ?5.584 + 2.747·log₁₀L（R2 = 0.579），异速指数b = 2.747显著低于等速生长（t = ?9.76, P = 2 × 10^?22），即较长玻璃鳗按比例更轻。战役间质量—体长斜率存在异质性（2.20–3.52），但长度标准化指数Kr*与Kr在战役层面高度一致（Pearson r = 0.988），脉冲内分解在采用脉冲特异斜率时基本不变，表明该异质性不扭曲本研究推断。

4.2 战役间及季节内体况

体况在战役间变化强烈（ANOVA F_14,8103 = 276.8, P < 0.001, η2 = 0.32），两季均呈现晚冬最低、秋冬至初冬最高的重复模式。平均体长和平均质量在两季均呈平行季节过程：自秋季早期到达者下降至仲冬至晚冬的中期最低值（约62–64 mm和约0.20–0.22 g），随后向晚程恢复。平均Kr描述相同的季节内隆起，随主要入湾脉冲达到峰值并向稀疏秋季前兆和晚冬尾部回落，支持H1。将两季对齐至9月1日以来的天数后，2014/15队列在迁移程匹配点全程较小（log₁₀体长差?0.0122, P = 0.001；log₁₀湿重差?0.0309, P = 0.006），但体长差可完全解释质量差，Kr季节偏移为+0.0022（P = 0.83），无证据表明两队列体况存在差异。

4.3 体况与入湾密度

体况倾向于在入湾峰值处更高，支持H2的预测方向。个体层面分层模型中，体况随每log₁₀密度单位增加+0.038 Kr单位（95% CI 0.004–0.073; P = 0.031）。战役均值相关给出相同信号（Pearson r = 0.56, P = 0.030；Spearman ρ = 0.53, P = 0.047）。但该耦合被报告为提示性而非确立性：15个战役仅嵌套于两个招募季节内（滞后1残差自相关ρ? = 0.41，有效n ≈ 6.3），在有效样本上相关不再显著（P = 0.23）；置信区间0.07–0.83跨越从可忽略到强。高密度与低密度战役的个体Kr分布比较显示，第10、50、90百分位均在高密度战役中更高且幅度相近（+0.018、+0.012、+0.009），离散度不变（SD 0.060对0.061），Kolmogorov–Smirnov统计量0.107——这是位置平移而非截断的特征。在季节尺度上，平均Kr与到达顺序无单调关联（Spearman ρ = ?0.014, P = 0.96），体况追踪的是脉冲的幅度和相位而非日历。

4.4 脉冲内动态

在每个冬季脉冲内，平均体况自早期峰值密度战役向晚期尾部逐次下降。事件层面估计中，2013/14年下降显著（?0.132 Kr单位/100天，95% CI ?0.190至?0.073, P = 0.010），两冬季脉冲合并分析亦显著（?0.098/100天，95% CI ?0.146至?0.051, P = 0.003；8个事件），2014/15年呈一致趋势（?0.062/100天，P = 0.074）。逐事件剔除诊断显示斜率全程保持负值。精确分解ΔKr = Δlog₁₀M ? b·Δlog₁₀L表明，2013/14年变化几乎完全由给定体长下的质量承载（质量下降的115%），2014/15年约一半由体长缩短承载（54%对46%）；该归因在采用脉冲特异斜率时不变。体长在两脉冲内均无显著趋势，而质量下降（2014/15显著：?0.074 g/100天, P = 0.006；2013/14倾向：?0.061 g/100天, P = 0.12）。H3在预测方向和质量—体长归因上得到支持，但证据属探索性。

4.5 入湾的环境过滤

在单次潮汐尺度上（n = 246网次，92天，29战役），两个预测变量在所有规格中均稳健。涨潮为最强效应（β = +0.474, 95% CI 0.312–0.636, P < 0.001），涨潮网次密度显著高于落潮——选择性潮汐流输运的特征，确证H4第一部分。热—季节复合效应强烈为负（β = ?0.439, P < 0.001），因水温和日历季节高度共曲（最坏情况共曲度≈0.93），该系数不能解读为独立热效应。盐度为负（β = ?0.192, P < 0.001），指示集中于淡水影响区。15天降雨为正但未显著（β = +0.159, P = 0.079），且规格依赖，被报告为提示性相关。流速、7天风速和昼夜主效应均不显著；强风为临界（β = +0.217, P = 0.066）。模型解释R2m = 0.625、R2c = 0.797。方差主要集中于战役尺度（σ2campaign = 0.131对σ2day = 0.016）。

4.6 潮汐×昼夜结构

尽管昼夜主效应为零，预设的潮汐×昼夜交互作用显著（LRT χ2? = 8.71, P = 0.003）：涨潮增强在白天更大，落潮渔获在夜间相对较高（单元均值：涨潮/昼1.74 > 涨潮/夜1.60 > 落潮/夜1.10 > 落潮/昼0.77）。这确证H4的其余两部分。浊度全程较高（中位数126.5 NTU），但既非密度的显著预测变量，也非昼夜效应的显著调节变量，因此浊度被处理为合理但未经证实的解释。

五、讨论与结论

本研究的核心组织性结果是：丰度与体况耦合但非冗余——大部分个体在脉冲峰值以优于平均的体况进入，但每个脉冲内部到达鳗鱼的给定体长质量下降。季节内体况隆起调和了表面矛盾的观测：丰度—体况正耦合、高密度时整个体况分布上移、脉冲内下降、以及季节尺度到达顺序效应缺失——体况追踪脉冲的幅度和相位而非日历。丰度—体况耦合被解读为对共同强迫的共享响应而非因果联系，可能机制为高密度和高体况共同标记近期输送事件，新到达且海洋储备完整的鳗鱼以数量优势抵达采样站；脉冲间期则由储备已耗竭的残留河口种群主导。脉冲内给定体长质量下降是最具信息量的结果，但样本为连续横断面而非同一鳗鱼的重复测量，无法单独区分河口内耗竭与内在更轻、更晚出生的海洋队列到达，也不能区分选择性定居、差异停留时间等组成变化；但该分析全程限定于阶段VB，排除了色素分期推进的混淆。质量—体长下降与已发表的未摄食阶段VB玻璃鳗能量消耗速率一致（18°C下40天损失约39%体质量），但Kr为形态测量指标，本研究未测量能量含量，所描述的是与耗竭一致的关联而非已证实的河口内过程。

环境模型提供了解读体况的局部背景。涨潮主导是STST的特征，负热—季节系数在温暖边缘窗口的下降支上得到解释，负盐度效应与淡水携带嗅觉线索的上游引导一致。潮汐×昼夜交互作用是唯一真正新颖的环境发现：昼夜主效应为零而涨潮增强在白天更大、落潮渔获在夜间较高，与持续高浊度下夜间光照线索缺失的假设一致，但同样可解读为新月采样下潮汐相位主导、昼夜效应较弱。

管理相关信息是：丰度和体况作为招募的描述符耦合但非冗余，仅数量监测遗漏了质量维度——到达个体的形态体况，以及推断的后续迁移潜力——而这些信息同一渔获中已经包含。这是潜力的展示而非验证的指标：将Kr确立为监测描述符需要更多季节、跨河口复制，并理想地以直接能量测量校准。本研究证明该信号存在、内部一致，且无需额外采样即可从单观察者、单阶段残差体况设计中获得。此外，这些数据涉及监测空白：瓜达尔基维尔向ICES/WGEEL指数提供科学序列仅至2007年，目前大西洋南缘已无站点纳入协调招募集，本研究为恢复该站点的监测代表性提供了方法学基础。