利用耳石地球化学与产卵红巢计数估算濒危公牛鳟(Bull Trout)集合种群的产卵溪流贡献率

《Ecology of Freshwater Fish》：Estimates of Spawning Stream Contributions to an Imperilled Bull Trout Metapopulation Utilizing Otolith Geochemistry and Redd Counts

　　

摘要：迁徙鱼类物种在整个水域栖息地中经历了大幅衰退，但管理和保护迁徙物种的行动需要更好地理解局部种群对集合种群的相对贡献。这类种群统计信息对于优先在最高产的产卵支流进行保护和修复至关重要。本研究的目标是估算斯旺湖(Swan Lake)公牛鳟育肥群体中来自各产卵支流的贡献，并将该信息与历史红巢(redd)计数相关联，以探讨局部种群的生产力与种群统计动态。研究人员利用耳石地球化学与随机森林(Random Forest)分类法，对斯旺湖湖鳟(Lake Trout)抑制行动中兼捕的公牛鳟进行了产卵溪流起源鉴定。公牛鳟的溪流归属结果与数十年来调查获得的历史红巢计数数据进行了比较。结果证明，部分产卵支流对斯旺湖公牛鳟的贡献不成比例，且该流域内不同溪流具有差异化的相对生产力，而这一点仅通过检查红巢计数是无法显现的。了解斯旺河流域产卵支流间公牛鳟的种群统计动态，将有助于确定有效的缓解行动优先级，以辅助这一重要集合种群的恢复。

论文解读：利用耳石地球化学与红巢计数估算濒危公牛鳟集合种群的产卵溪流贡献

一、研究背景与意义

内陆淡水鱼类正日益受到气候变化、栖息地丧失及非本土物种相互作用的威胁。公牛鳟（Salvelinus confluentus）是北美太平洋西北部和加拿大西部的特有种，在美国受《濒危物种法》保护（列为“受威胁”），在加拿大被列为特别关注或受威胁。斯旺湖（Swan Lake）与斯旺河（Swan River）流域是公牛鳟的重要集合种群（metapopulation）栖息地。集合种群由潜在迁移连接的离散局部种群组成，其存续依赖于局部种群的动态平衡。

目前，资源管理部分常利用秋季产卵红巢（redd，雌鱼在砾石中挖掘的产卵坑）计数作为成体丰度的代理指标来监测公牛鳟趋势。然而，红巢计数存在观测误差，且雌鱼可能在一个季节内多次迁移并建造多个红巢，导致对局部种群生产力的估算存在偏差。与此同时，公牛鳟受非本土引入的湖鳟（Salvelinus namaycush）竞争与捕食压力巨大，斯旺湖的湖鳟抑制行动虽移除近6万尾，但公牛鳟红巢计数仍持续下降约90%。因此，理解不同产卵支流对湖泊育肥群体（adfluvial life history，江湖型生活史）的实际贡献，对精准保护至关重要。本研究发表于《Ecology of Freshwater Fish》，旨在结合耳石（otolith，鱼耳石）地球化学指纹与历史红巢数据，验证红巢计数是否能真实反映产卵溪流的贡献率。

二、主要关键技术方法

研究人员利用斯旺湖湖鳟抑制行动（2008–2016, 2022至今）的兼捕致死公牛鳟存档耳石（N=87）作为研究样本。为建立产卵溪流地球化学基线，研究人员采集了12条公牛鳟产卵溪流及斯旺河上游的水样，并以非洄游、家域小的黏滑杜父鱼（Slimy Sculpin, Cottus cognatus）作为替代种，采集其耳石（每溪约10尾）建立⁸⁷Sr/⁸⁶Sr、Sr/Ca及Ba/Ca比值的参考数据库。水样与耳石经激光剥蚀电感耦合等离子体质谱（LA-ICP-MS）及热电离质谱（TIMS）分析。研究人员从公牛鳟耳石核心外初次摄食检查点外延50μm区域提取幼体期地球化学信号，使用随机森林（Random Forest）分类模型（500棵树，70%训练/30%测试）将公牛鳟个体分配至其natal origin（natal spawning tributary，原生产卵溪流）。通过自举法（1000次）量化分配不确定性，并将分配的溪流贡献比例与2001–2019年对应的历史红巢计数比例进行皮尔逊相关性分析。

三、研究结果

3.1 水化学与杜父鱼耳石地球化学

斯旺河流域不同产卵溪流及干流位点的水体⁸⁷Sr/⁸⁶Sr同位素及元素/Ca比存在显著差异。Soup溪的⁸⁷Sr/⁸⁶Sr最低（0.72842±0.00105），Lion溪最高（0.73697±0.00086）。水体Sr/Ca在上斯旺河与南Woodward溪约为Elk溪的2倍以上，Ba/Ca在上斯旺河与Squeezer溪约为Woodward溪的7倍。

黏滑杜父鱼耳石的⁸⁷Sr/⁸⁶Sr、Sr/Ca及Ba/Ca在不同溪流间均存在极显著异质性（Kruskal-Wallis, p≤0.001），且耳石值与对应溪流的水化学值高度线性相关（⁸⁷Sr/⁸⁶Sr, R2=0.75; Sr/Ca, R2=0.56; Ba/Ca, R2=0.74），证明杜父鱼耳石可有效记录当地水文地球化学信号，适合作为公牛鳟幼体期生境的替代基线。

3.2 公牛鳟年龄结构

斯旺湖采集的公牛鳟年龄跨度为2–7龄，以3、4、5龄为主（约占87%），其中5龄鱼最频（32.6%）。年龄与体长呈显著正相关（R2=0.62, p<0.001）。

3.3 集合种群统计与溪流贡献

随机森林模型整合⁸⁷Sr/⁸⁶Sr、Sr/Ca、Ba/Ca时预测准确率最高（整体测试集准确率88%，袋外误差OOB=17.98%，Kappa=0.87）。部分溪流（Goat, Upper Swan River, Lion）训练集分类误差为0，而Squeezer与Woodward溪误差较高。

基于模型对存档公牛鳟耳石的分配，近58%的个体原生源于3条溪流：Elk溪（中位数贡献0.22, 95% CI: 0.12–0.28）、Jim溪（0.18, 0.07–0.32）与Lion溪。Upper Swan River无分配记录。Goat溪分配比例最低（0.02, 0.0–0.06）。

将分配比例与同期（2001–2019）红巢计数比例对比，二者呈正相关（Pearson r=0.60, p=0.086）。Lost溪与Elk溪的红巢比例与耳石分配比例基本一致。Jim溪红巢占比约9%，却贡献了21%的耳石分配个体；反之，Goat溪红巢常年占约10%，耳石分配仅2%。Piper、Lost、Soup等溪红巢计数低且分配比例也低，显示其小种群维持了与产出相称的向集合种群的输入。

四、讨论总结与结论翻译

讨论部分总结

研究人员发现，大部分溪流（如Lost、Elk）的红巢计数与耳石推算的后期生活史贡献相符，表明红巢可大致反映这些局部种群的存续与输入。Soup、Lost、Piper等小种群虽红巢少，但维持了成比例贡献，在优势种群（Elk、Lion）衰退时可能成为集合种群动态的重要缓冲，值得保护。

中型种群（Woodward、Squeezer、Lion溪）红巢计数陡降且在耳石分配中代表不足，可能正在失去未来繁殖群，需加强监测。

Jim溪与Goat溪是主要例外：Jim溪红巢占比低但耳石分配高（模型对该溪误分最多，置信区间宽），Goat溪红巢稳定但耳石分配极低（模型分类准且误分率低）。这可能与不同生活史表达（如幼鱼离溪后转为河流型fluvial）、迁移距离、密度依赖生长，或非本土Brook Trout（Salvelinus fontinalis）的竞争捕食导致幼体存活率差异有关。若Goat溪红巢多但存活到湖泊的成鱼少，说明限制发生在孵化后到入湖前的阶段，保护重点应是提高幼鱼存活而非单纯改善产卵栖息地；反之，在Soup、Lost、Piper溪，产卵栖息地的提升可能直接增加混合群体。

与既往微卫星遗传研究（DeHaan et al., 2011）对比，耳石分配显示最远的Elk溪贡献最大，而非近湖溪流，部分差异可能源于遗传方法依赖生殖隔离假设及样本构成不同。研究局限包括使用替代种耳石基线、幼鱼可能在其他支流重养导致信号偏移、以及部分溪流（Woodward, Squeezer）分类误差较高。未来可结合遗传与地球化学，或试用Brook Trout作替代种提升置信度。管理上，应在Elk等高产源种群保护生境，在Goat等产卵多但存活低的溪流针对幼鱼瓶颈（如控制非本土鳟竞争）行动；同时继续优化斯旺湖湖鳟抑制策略（如针对未转食鱼类的阶段、提升刺网效率、减少公牛鳟兼捕），以协同恢复集合种群生产力。

研究结论部分翻译

耳石地球化学与历史红巢计数数据的结合，为深入理解斯旺河流域公牛鳟集合种群中各产卵溪流的相对贡献提供了有效途径。相关已有研究同样利用耳石中记录的物理与环境信息来描述种群动态，而本研究的附加红巢数据使研究人员能够将局部种群的相对生产力与斯旺河流域的公牛鳟现状进行对比。该方法增强了研究人员在原生产卵与育肥溪流中确定保护与修复优先序的能力。