栖息地破碎化下鱼类种群密度响应的评估框架：五大湖支流低水头坝（low-head dams）上下游效应与 t50 滞后

《Oikos》：A framework for investigating fish population density responses to habitat fragmentation

　　

全球范围内，依赖河流的鱼类种群正受到河道内屏障（in-stream barriers）造成的破碎化，但负向种群响应对破碎化的普遍程度仍不清楚。研究人员开发了一个建模框架，用于估计鱼类在河道内屏障上游和下游的密度响应，并将其应用于了解低水头坝（low-head dams）如何影响 24 条 Great Lakes 支流的鱼类群落。在建模的 29 个物种中，研究人员发现强证据支持低水头坝上游栖息地中种群密度增加（n=6 种）和降低（n=6 种），而仅发现强证据支持坝下游种群增加（n=16 种）。坝诱导的鱼类种群密度变化完成 50% 所需时间，即 t₅₀，为数十年的量级；平均 t₅₀ 在低水头坝下游栖息地为 41 年，在上游栖息地为 21 年。t₅₀ 的估计提示负向效应比正向效应发生更快。然而，该比较的检验效能仍有限，因为用于推导 t₅₀ 的逻辑斯谛种群增长模型（logistic population growth model）仅对 29 种鱼中的 1 种具有强证据。研究人员提出不同机制解释观察到的响应：正向响应由群落结构（assemblage structure）变化驱动，负向响应由迁入（immigration）减少驱动。该研究展示了考虑集合种群动态（metapopulation dynamics）以及屏障建设与鱼类种群响应之间潜在时间滞后的框架，如何推进对河道内屏障如何构建河流中鱼类种群分布与多度的理解。

研究背景与问题：全球依赖河流的鱼类种群正受到河道内屏障（in-stream barriers）造成的栖息地破碎化（habitat fragmentation）影响，小型屏障包括低水头坝（low-head dams）、涵洞和闸门。鱼类种群密度对建坝的响应受生态背景和多种互作驱动因素影响，负向响应是否普遍仍不清楚。在五大湖流域，低水头坝与入侵海七鳃鳗控制密切相关，约 460 个最下游屏障用于限制成年海七鳃鳗进入上游产卵和育幼栖息地；但拆坝、建坝或鱼道改造政策需考虑对期望鱼类的附带影响。因此，研究人员需要构建基于集合种群动态（metapopulation dynamics）和灭绝债务（extinction debt）思想的采样与分析框架，评估低水头坝对五大湖支流鱼类种群密度的影响。

研究人员开展了什么研究及结论意义：研究人员开发了估计鱼类对栖息地破碎化密度响应的建模框架，并应用于 24 条 Great Lakes 支流低水头坝影响下的鱼类群落。目标是确定哪些鱼类的密度与破碎化强相关，量化屏障上下游的效应，并估计建坝到密度变化之间的滞后。结论表明，约一半受研究物种的密度与低水头坝强相关；上游同时出现正向和负向响应，下游强支持关联均为正向，松弛时间（relaxation time）对负向效应可能短于正向效应。该研究意义在于为拆坝和建坝权衡提供可纳入决策支持的工具，并表明河道内屏障会通过集合种群动态和时间滞后影响鱼类分布与多度。

关键技术方法概览：研究采用配对研究设计（paired study design），将 24 个含低水头坝溪流与未建坝或仅在高位建坝的参考溪流配对，共 47 条溪流、24 个配对；样本来自北美五大湖盆地（Great Lakes Basin），1996 年调查、1997 年部分重采样、2023 年按相同方法重采样，仅 4/274 站因高水失去样。采样使用背包电捕鱼，记录鱼种和总长；栖息地测宽、最大深度、Wentworth 底质、温度、电导率。以每单位努力捕获量（catch per effort, CPE）为密度指标，主成分分析（principal component analysis, PCA）提取 PC₁、PC₂；用负二项模型、贝叶斯 MCMC（Markov chain Monte Carlo, MCMC；Stan/RStan）、贝叶斯因子（Bayes factor, BF）和 bridge sampling 比较无效应 M₁、逻辑斯谛动态 M₂ 和阶跃 M₃。

研究结果：模型支持方面，29 个物种中，无屏障效应模型 M₁ 仅对 8 个物种有正支持；15 个物种仅含屏障效应模型获支持，其中 M₃ 对 12 个物种获支持，M₂ 对 1 个物种 white sucker 获支持，6 个物种证据不确定。上下游效应方面，下游效应量 α>1 见于 97% 的鱼（28/29），其中 16/29 的 90% 可信区间（credible intervals, CI）不含零；最大下游密度增加为 white sucker Catostomus commersonii（均值 40.8）、blacknose shiner Notropis heterolepis（24.3）和 American brook lamprey Lethenteron appendix（21.5）。上游效应变异更大：α>1 见于 59%（17/29），90% CI 不含零者 6/29；α<1 见于 41%（12/29），90% CI 不含零者 6/29。上游最大增加为 American brook lamprey（22.0）、blacknose shiner（14.0）、central mudminnow Umbra limi（9.9）；最大下降为 round goby Neogobius melanostomus（0.03）、yellow perch Perca flavescens（0.04）、logperch Percina caprodes（0.16）。t₅₀ 方面，对 19 个最可信含屏障效应模型的物种进行探索性估计：下游平均 t₅₀ 为 41 年，上游平均为 22 年；下游最低 25 年，上游有 41% 物种（7 种）t₅₀ 小于 10 年，且上游快速响应者均为负关联（α<1）。细尺度栖息地关联方面，72%（21/29） modeled 物种的密度与细尺度栖息地测量强关联；38%（11/29）与 PC₁ 强关联，48%（14/29）与 PC₂ 强关联，14%（4/29）与二者均强关联；PC₁ 强关联者多呈正关联（73%），PC₂ 强关联者强烈偏向负关联（93%）。

讨论部分总结：研究人员认为，约一半受研究鱼类的密度与低水头坝强关联；上游既有正向也有负向响应，下游强支持关联全为正向，松弛时间对负向效应短于正向效应，但该比较为探索性且逻辑斯谛模型对多数物种支持较弱。上游密度下降可能与鱼类在上下游间迁移受阻有关，如 round goby 存在季节性迁移，低水头坝阻断产卵成鱼上溯并改变幼鱼补充分布。上游密度增加机制未被研究明确识别；原文提出捕食释放和种间竞争释放可能是原因，American brook lamprey 的上游增加则可能与针对入侵海七鳃鳗的化学农药非目标死亡降低有关。下游未支持关键栖息地丧失导致下降的假设；可能原因是所研究低水头坝并非完全迁移屏障，关键栖息地分布于坝上下游，或源-汇动态（source–sink dynamics）掩盖了补充失败。研究人员假设参考栖息地鱼群处于平衡密度，并假设由未建坝到建坝平衡的过渡遵循逻辑斯谛种群增长；若种群本身增长或下降，破碎化效应可能被掩盖。未来可放宽过渡动态函数形式，采用广义逻辑斯谛增长模型，发展 n-state 集合种群模型和发生率函数模型（incidence function model, IFM），并进行重复纵向采样。研究还存在选择偏差，仅分析常见种和可涉水栖息地，未代表全部五大湖河流鱼类群落。

研究结论翻译：研究人员开发并应用了一个估计五大湖支流低水头坝对鱼类种群密度响应的框架，为权衡拆坝或建坝提供潜在工具。例如，这些估计可纳入优化区域屏障移除策略的决策支持工具，或用于特定溪流和决策背景的详细定量模型。研究还表明，为限制入侵海七鳃鳗上游迁移而设计建造的屏障限制了 round goby 的分布，而 round goby 是威胁五大湖支流本地鱼类群落的另一入侵鱼。研究人员强调，分析未代表全部本地和归化五大湖鱼类群落；对稀有或受威胁物种、春秋迁移期以及可涉水和不可涉水栖息地的针对性采样，可建立对小尺度河道内屏障影响五大湖鱼类的全面理解。框架应用表明，在没有建坝前后鱼类群落响应数据的应用中，该框架可有效解决鱼类种群对河道内屏障响应的不确定性。基于模拟的检验和在其他系统中的应用，可明确框架有用决策背景范围，以及稳健估计河道内屏障效应和松弛时间所需采样要求。支持生态理论定量检验的方法学进展，对理解鱼类种群密度如何响应改变支流鱼类移动的保护行动之变异性和一致性至关重要。