《Global Ecology and Conservation》：Distributional responses of benthic fish across contrasting dam diversion conditions and during a heatwave: associations with water temperature
编辑推荐：
专业术语翻译摘要部分：
引水坝分流（diversion dam）与热浪（heatwave）条件下底栖鱼类（benthic fish）的分布响应：与水温度的关联性研究
绕行段（bypass reach, BR）升温与冷水狭温性（cold-stenothermic
专业术语翻译摘要部分：
引水坝分流（diversion dam）与热浪（heatwave）条件下底栖鱼类（benthic fish）的分布响应：与水温度的关联性研究
绕行段（bypass reach, BR）升温与冷水狭温性（cold-stenothermic）鱼类分布相关，但与广温性（eurythermic）鱼类分布无关；研究人员发现，在引水停止两年后，野吕川杜父鱼（Cottus nozawae）的单位面积捕获量（catch per unit area, CPUA）相对于未受调控的参照河段出现显著的部分恢复；在热浪年份，尽管该物种与冷水栖息地相关联，其CPUA仍在广泛区域内上升；尾水渠（tailrace, TR）可能为C. nozawae提供热避难（thermal refuge）栖息地，值得进一步研究。
论文解读：
一、研究背景与科学问题
水温度是淡水生态系统结构与功能的核心驱动因素，影响初级生产力、溶解氧溶解度、食物网动态以及水生生物的分布、迁移和死亡率。淡水生态系统对水温变化尤为敏感，尤其是夏季升温，这源于水生变温动物（ectotherm）的生理特性、不成比例的高生物多样性以及污染、栖息地退化和破碎化等多重胁迫的叠加。近年来，河流水温受到人为因素的显著改变，主要包括大尺度气候变化和局部水坝建设。引水坝（diversion dam）通过将河流径流的大部分甚至全部引走用于水电、灌溉或城市供水，在下游形成径流大幅减少的绕行段（bypass reach, BR），同时引水经地下管道输送至下游非一体化的发电站后以尾水形式部分返回河流，形成尾水渠（tailrace, TR）。已有研究多聚焦于径流减少对鱼类丰度和群落结构的影响，而BR中热力条件的生态后果研究相对不足，尤其是BR与TR之间的热力异质性及其对鱼类分布的影响。
鱼类作为关键的水生变温脊椎动物，不同物种在热偏好和耐受限度上存在显著差异。冷水狭温性物种对升温尤为脆弱，而底栖鱼类通常具有活动性低、系统发育独特性高等特征，是检测热胁迫的理想指示物种。然而，以往研究在分类和功能上存在偏倚，主要集中于具有商业和休闲价值的移动性物种（如鲑科鱼类），对定居性底栖鱼类的响应研究严重不足。
2023年北半球经历了过去两千年中最温暖的副热带夏季，这一极端热浪事件与日本北海道丰平川（Toyohira River）流域的毛岩坝（Moiwa Dam）引水停止事件在时间上先后发生，为评估鱼类在复合水文与热力变化条件下的响应提供了难得的自然实验条件。该研究区域包含一个约12.85 km的BR，引水停止前BR水温显著升高，停止后恢复至接近自然状态。研究人员利用这一独特条件，针对两种热生态位截然不同的优势底栖鱼类——冷水狭温性的野吕川杜父鱼（Cottus nozawae）和相对广温性的石泥鳅（Barbatula oreas），提出两个核心问题：（1）物种的空间分布是否与引水及恢复水流相关的主河道热力梯度相关联？（2）这些分布在引水停止后两年内及热浪期间是否发生变化、如何变化？
二、主要技术方法
研究在2022年10月、2023年7月、2023年9月和2024年10月四个时点，采用三人小组背负式电捕鱼器（backpack electrofisher）在11个样点（8个主河道样点、3个支流样点）开展单次半定量电捕调查，每样点设3个重复（间距≥5 m），以单位面积捕获量（catch per unit area, CPUA）作为相对丰度指标。环境变量包括电导率（EC）、pH、溶解氧（DO）、水深、流速和底质特征等在重复内三点测量取均值。日均水温（daily mean water temperature, DMWT）通过一维稳态水流水力模型（HEC-RAS v5.0.7）模拟获取，支流样点水温基于站点特定的气温-水温关系估算。统计分析采用广义线性混合模型（GLMMs），包括时空分布分析（含双重差分（difference-in-differences, DiD）分析）、环境关联分析（含预测因子筛选、模型平均和交叉年模型）以及时间响应分析（基于基线环境-分布关系预测后续年份CPUA）。
三、研究结果
3.1 时空分布分析
对于C. nozawae，全局GLMM检测到显著的河段-时点交互作用，各时点特异性模型均确认显著的河段效应。2022年10月和2023年7月，D2（主河道-尾水交汇段）和TR的CPUA显著高于D1（下游末端）和BR；2023年9月和2024年10月，D2和TR仍保持最高CPUA，与D1的显著差异持续存在，而BR的CPUA逐步上升，至2024年10月BR显著高于D1。基于参照河段的DiD分析显示，BR在引水停止后两年内相对于未受调控的MK参照河段出现显著的CPUA增加（比例DiD：RRR = 7.97，校正p < 0.01；加性DiD：估计值 = 4.65 inds/100m2，校正p < 0.01）。对于B. oreas，河段效应不显著，未检测到河段-时点交互作用，表明其空间分布格局在时点间无明显变化。
3.2 环境关联分析
预测因子筛选为C. nozawae保留了DMWT、DO和Depth，为B. oreas保留了DMWT、DO和pH。对于C. nozawae，仲夏平均DMWT在三个秋季时点特异性模型中均与CPUA呈显著负相关；2023年7月时点，当年早夏平均DMWT不显著，而前一年（2022年）仲夏平均DMWT呈显著负相关且RMSE更低，表明其分布保留了与前一年仲夏热力格局的滞后关联。除DMWT外，仅DO在2023年9月显示显著正相关。秋季模型预测精度以2022年10月最高（RMSE最低），此后RMSE上升，表明环境-分布对应关系减弱。全局秋季GLMM检测到显著的时点-DMWT交互作用（LRT：χ2 = 6.27，df = 2，p = 0.044），DMWT斜率在2022年10月显著为负（β = -0.54，95% CI = -0.90至-0.17），2023年9月不显著（β = -0.17，95% CI = -1.05至0.72），2024年10月为负但窄幅不显著（β = -1.21，95% CI = -2.43至0.01）。对于B. oreas，所有时点均未检测到任何环境变量（包括DMWT）与CPUA的显著关联。
3.3 时间响应分析
基于2022年10月基线环境-分布关系的模型加权预测在2023年9月和2024年10月均系统性低估了主河道样点的观测CPUA，其中2023年9月DS-3样点的低估最为显著。这一模式与观测到的CPUA上升同时发生，尽管2023年9月所有样点及2024年10月非BR样点的仲夏平均DMWT均高于2022年，即出现了“温度升高但丰度增加”的时空格局，与环境中DMWT-CPUA的负向空间关联形成鲜明对比。体长频率分布显示，2023年9月的CPUA升高不仅限于最小个体，也涵盖较大体长等级。
四、讨论与结论
研究结果支持第一个假设：引水停止前，冷水狭温性C. nozawae的CPUA呈现显著的河段尺度空间变异，并与仲夏平均DMWT呈显著负相关，较高CPUA出现在较冷河段；而相对广温性的B. oreas未检测到显著的温度关联，空间结构弱。第二个假设也得到支持：引水停止和BR相对冷却后，C. nozawae在BR的CPUA相对于未受调控的MK参照河段显著增加，但恢复是部分的——BR与下游交汇段D2之间仍存在数值上的CPUA差异，且后续年份环境-分布对应关系减弱，观测CPUA在多个BR样点仍低于模型预测值。第三个假设未获支持：2023年热浪后C. nozawae的CPUA在大多数样点上升而非下降，超过基线环境-分布关系的预测，与其对DMWT的负向空间关联形成对比。
研究还发现，C. nozawae在超过其野外分布温度上限（约18°C）的BR中仍以低密度持续存在，且不限于特定体长等级，提示可能存在热避难所效应或资源缓冲等机制，但具体过程尚需遗传学、生理学、个体追踪和猎物资源评估等综合研究加以验证。DS-3样点预测偏差最大，可能与上游鱼道（约1 m高）将河道从约80 m缩窄至8 m、夏季平均流速增至约3 m/s的物理屏障效应有关，可能限制底栖鱼类的纵向连通性和向上游适宜热栖息地的移动。
总体而言，仲夏平均DMWT是C. nozawae分布的最一致环境关联因子，而水力和底质变量未显示显著关联。研究结果支持水温在观测的时空分布模式中发挥重要作用，同时表明引水停止后的种群恢复是多因素共同作用的渐进过程，需要更长时间尺度的监测来全面评估引水坝调控恢复措施对冷水狭温性底栖鱼类的长期生态效应。
