本文来自微信公众号:中国工程院院刊(ID:CAE-Engineering),本文选自中国工程院院刊《Engineering》2021年第2期,作者:陈求稳、张建云、陈宇琛、莫康乐、王骏、唐磊、林育青、陈磊、高勇、姜伟、张与馨,原文标题:《建坝河流鱼类繁殖的流速调控丨Engineering》,题图来自:视觉中国


保护濒危和重要经济鱼类以及控制入侵鱼类是全世界水利工程的巨大挑战。当前,已有共识是流速会影响鱼类在河流中的产卵,但仍缺乏可靠的科学支持和相关机制,以综合考虑产卵、受精、孵化、存活等因素,建立流速与鱼类繁殖之间的定量关系。


中国工程院张建云院士科研团队在中国工程院院刊《Engineering》2021年第2期刊发《建坝河流鱼类繁殖的流速调控》,通过实验室控制实验和野外现场实验,量化了流速与家鱼繁殖率之间的关系。


文章研究结果表明,促发家鱼产卵需要一个最小的速度,并且需要一定的速度范围以维持其产卵过程。但是,实验发现,随着流速的增加,家鱼的胚胎孵化和仔鱼发育会受到抑制。


综合考虑产卵和孵化以及仔鱼存活需求,文章确定了河流中家鱼繁殖的优化流速过程。该发现对改进水库运行以构建合理、精确的生态流量来调控鱼类繁殖具有重要意义,并在三峡工程应用中取得了有前景的效果。


一、引言


地球上大多数大河都被密集地筑坝,以进行水力发电、洪水管理、供水和航行,并且大坝的数量仍在快速地增长。筑坝不仅改变了自然水文状况、泥沙通量和水质,而且还对鱼类群落产生了不良的生态影响。大坝的运行会影响本地鱼类的产卵活动,从而影响其数量和种群,这是世界各地的河流都会有的现象,如马德拉河、科罗拉多河和澜沧江-湄公河。


鱼的产卵活动很可能是由多种环境因素触发的,包括流速、水温、溶解氧、浊度和初级产物。众所周知,流速对于产漂流卵的鱼类的产卵特别重要。


中国四大家鱼包括鲢(Hypophthalmichthys molitrix)、草鱼(Ctenopharyngodon idellus)、青鱼(Mylopharyngodon piceus)和鳙(Aristichthys nobilis),是中国渔业养殖和捕捞的重要鱼类。在20世纪,长江全部淡水鱼的捕捞量中有60%是家鱼。


但是,在过去的几十年中,由于三峡大坝的建设引起了水文条件的显著变化,鱼的产卵活动受到了严重影响。2003年,三峡水库(TGR)开始运行时,鱼卵和仔鱼的数量仅为2002年的10%。


而在其他区域,家鱼已成为具有威胁性的入侵物种之一。由于它们较高的繁殖和运动能力,其种群迅速扩大,已成功入侵了许多河流和湖泊。家鱼在美国的入侵和种群激增导致本地鱼类种群锐减。在澳大利亚许多河流中,家鱼维持着高生存量(高达 3144 kg·hm–2


在中国,鳙和鲢的引入对高原湖泊中的许多本地或特有物种造成了破坏性影响,导致了不可逆转的生态变化。尽管已经采取了多种措施来控制家鱼,澳大利亚甚至提出了一种高风险的想法,通过释放病毒以防止家鱼种群扩张,但收效甚微。


调控繁殖过程,包括产卵、受精、孵化和胚胎发育,可能是维持种群或控制家鱼入侵的有效措施。但是,到目前为止,研究仅集中在产卵上,这些研究表明家鱼产卵发生在浑浊的河流中,最低水温要求为18 ℃,还有特定的流速要求。综合考虑产卵、受精、孵化、仔鱼存活等繁殖过程,建立流速与家鱼繁殖之间的定量关系,目前还缺乏坚实的科学依据。


本研究以家鱼为目标物种,通过物理模型验证两个假设:①流速与产卵行为以及家鱼胚胎发育之间存在定量关系;②存在一个最佳的流速来平衡产卵行为和胚胎发育的要求。同时,开展了多年的家鱼产卵野外试验,与物理模型实验结果相互印证。本研究通过对河流水利工程适应性管理,为提升或控制鱼类繁殖提供坚实的基础。


二、材料与方法


所有实验均在南京水利科学研究院当涂实验中心进行,并获得研究机构动物管理与使用委员会(IACUC)的许可。


1. 野外试验


在2012—2016年的野外试验性生态调度过程中,我们在长江干流宜都市附近的固定断面处采集了漂浮性鱼卵样品(图1)。在繁殖季节,经常会在该断面监测到家鱼的仔鱼和流卵。


据报道,从宜昌到杨家嘴的断面上游约10 km是产卵场。该河段宽近1 km,既没有易蚀的浅滩,也没有狭窄的地段,这些特殊的形态特征创造了有利于家鱼产卵的水流条件。


图1 长江中游鱼卵取样的研究河段河床标高和横截面(CS1~CS4)
图1 长江中游鱼卵取样的研究河段河床标高和横截面(CS1~CS4)


在为期五年的试验性生态调度期间,采样于每年的五月中旬开始,当时的水温高于18 ℃的阈值。连续两周采样,直到没有发现鱼卵为止。在采样期间,每天的8:00—9:00和14:00—15:00收集鱼卵。用网(直径为 0.5 m,网眼为500 μm)收集漂流的卵。网被悬挂在水面上,并逆流向上游定向。同时,采用旋转式流速仪用于测量流速。研究河段上游边界处的流量(m3 ·s–1和下游边界处的水位(m)从长江水利委员会获得。


收集到的鱼卵会被立即检查。直径为4~5 mm的卵被认为是潜在的家鱼卵。收集时都确定了鱼卵的发育阶段和受精后的时间。将鉴定出的鱼卵在实验室中进行培养,直到确认它们的种类为止,从而最终估算出每个采样的鱼卵密度(每1000 m3 鱼卵数量)


为获得整个河段的流场,建立了二维非恒定流水动力模型。通过获取的流量和水位数据以及在CS1~CS4 (图1)四个断面处测得的速度对模型进行率定和验证,附录A第1节中介绍了模型基本信息、率定和验证等信息。然后使用校准后的模型来计算整个河段的流场。


根据流速和鱼卵胚胎发育时间(漂移时间),计算至鱼卵收集段面的漂移距离,确定上游产卵场。然后,从流场模拟结果中提取确定的产卵场的流速,将其用于统计分析。


2. 实验室控制实验


实验室控制实验分别于2017年和2018年5月15日至6月20日进行。实验装置是一个大型的环形水槽[图 2(a)]。安装了两台水下推流器和四台水下水泵提升流速。水槽顺直部分被加宽加深,以方便安置水下推流器的叶轮,并且在两侧都安装了渔网,以防止鲢进入叶轮区域。


在水槽周围布置了22个测量点,并通过旋转式流速仪(南京水利水文自动化研究所,LS20B)测量了0.5 m深度处的流速。水槽周围有6个水下摄像头,以观察家鱼的产卵行为。实验过程中,水深设定为0.9 m,水质与驯化池的水质一致。使用多功能水质分析仪(美国YSI Inc.,YSI 650MDS)测量水质参数。实验设计了6种具有不同理论最大流速的方案(见附录A中的表 S1),每个方案每年重复两次。


性成熟的雄性[长度:(37±1.3) cm;体重:(3.6± 0.4) kg;数量:72尾]和雌性[长度:(43±2.1) cm;体重: (4.5±0.6) kg;数量:72尾]家鱼来源于湖北省黄冈市的家鱼孵化中心。实验之前,将这些家鱼在驯养池中驯养,每天两次提供商业饲料以满足其饮食的需求。7天后随机选择6尾雄性和6尾雌性家鱼作为实验组。


在移入环形水槽之前,向选定的家鱼注射一次促黄体生成激素释放激素拮抗剂(LHRHA2)。雌性家鱼注射剂量为2 μg·kg–1,雄性家鱼注射剂量减少一半,注射位置在胸鳍的底部。这些家鱼要经过三天的测试,或者直到产卵为止。


每次测试后,使用亚致死剂量的NaOH中和的 3-氨基苯甲酸乙酯甲磺酸盐(MS-222)麻醉雌性家鱼,进行解剖,根据卵巢饱和度确认发生过产卵行为的雌性的数量。


为了研究流速对胚胎发育的影响,建立了一个小的环形槽[图2(b)],它使用小型水泵来控制流速。于 2018年6月测试了三种方案,最大流速分别为0.4 m·s–1、 0.8 m·s–1和1.2 m·s–1。每个方案处理重复三次。鱼卵孵化后24 h新孵化的仔鱼的总孵化率和存活率用于表征胚胎发育。


以与上述相同的方式监测水质参数。在小的环形水槽中,水位保持恒定,每天更换一半的水。如果产卵受阻,则从大水槽中选出200颗鱼卵,以评估每次测试中胚胎的发育。


为了获得整个环形槽的流速,在研究中使用了计算流体动力学(CFD)方法,该方法在CFX 11.0(美国 ANSYS公司)中完成。通过从22个测量点观察到的流速数据对模型进行校准和验证。附录A中的第3节介绍了模型原理、率定和验证的详细信息。


图2 产卵试验深度为1.8 m的大环形水槽(a),以及胚胎实验深度为0.2 m的小环形水槽(b)<br label=图片备注 class=text-img-note>
图2 产卵试验深度为1.8 m的大环形水槽(a),以及胚胎实验深度为0.2 m的小环形水槽(b)


3. 实验数据的统计分析


通过单方差分析(ANOVA)研究了来自实验室和野外试验的鱼卵数据以及与之对应的流速数据,以确定流速对家鱼产卵和胚胎发育的影响。使用Duncan多差距测试和Tukey-t 检验对显著性进行了分析,显著性值(P)设为0.05。数据表示为平均值±标准误差(SE)。所有分析均使用SigmaPlot 11.0(Systat Software Inc.,美国)进行。


三、结果


1. 长江中游家鱼繁殖野外试验


通过五年的野外试验,监测到20次产卵事件(虚线框),其中发现了11个峰值产卵事件(较高的虚线框)(图3)。对每个产卵事件,都计算了持续时间(图3中每个粉红色块的宽度)。然后获得平均持续时间,大约为4.5 d。


图3 长江中游2012–2016年的产卵事件及相应的流速和流速上升率。黑色曲线表示产卵点的流速。点是观察到的鱼卵密度。红色虚线表示产卵事件。蓝点表示在触发水流速度下的鱼卵密度,该速度定义为发生在产卵事件中的第一个产卵动作时的流速。红点是在合适的速度下发生显著产卵的密度。蓝色箭头指示产卵事件期间的流速上升率(a)<br label=图片备注 class=text-img-note>
图3 长江中游2012–2016年的产卵事件及相应的流速和流速上升率。黑色曲线表示产卵点的流速。点是观察到的鱼卵密度。红色虚线表示产卵事件。蓝点表示在触发水流速度下的鱼卵密度,该速度定义为发生在产卵事件中的第一个产卵动作时的流速。红点是在合适的速度下发生显著产卵的密度。蓝色箭头指示产卵事件期间的流速上升率(a)


在采样期间,水温处于合适的水平,范围从18.5 ℃到24.0 ℃。鱼卵的鉴定表明,收集的卵中约55%属于鲢,38%属于草鱼,5%属于黑鱼,只有1%属于鳙鱼。结果还表明,不同家鱼的鱼卵比例每年差别不大,这表明五年之间比例没有显著变化。


明确了鱼胚的发育阶段,以确定鱼卵的漂移时间。被检验的鱼卵主要在囊胚发育的最后一个阶段,占总数的 72.42%,而在其他发育阶段的鱼卵占27.58%。


流体动力学模型的率定和验证结果如图S3所示。可以看出,验证过程中流速和水位的相对均方根误差分别在0.7%~4.1%和0.15%~0.23%范围内,这表明该河流段的模拟流场是可靠的。测算的产卵地点表明,产卵事件的地点在胭脂坝和红花套之间[图4(a)],与以前的研究一致。


从水动力模型结果中提取出每个定位地点产卵事件的相应流速,并与图3中的卵密度数据结合分析。统计结果(图4)显示,在所调查的河段中,家鱼产卵的触发流速为(1.30±0.20)m·s–1 [图4(b)],而家鱼产卵的适宜流速为(1.40±0.23)m·s–1 [图4(c)]


图4(a)研究河段的流速和产卵事件的产卵地点;在田间试验中,在长江中游的触发速度(n = 20)(b)和优先流速(n = 11)(c)。SD:标准偏差;n:产卵事件的数量<br label=图片备注 class=text-img-note>
图4(a)研究河段的流速和产卵事件的产卵地点;在田间试验中,在长江中游的触发速度(n = 20)(b)和优先流速(n = 11)(c)。SD:标准偏差;n:产卵事件的数量


2. 水槽模拟实验中家鱼的繁殖


图5给出了不同实验情景下大型环形水槽的流场。由于在最大速度分别为0.8 m·s–1和1.8 m·s–1的实验场景中未观察到产卵行为,因此未进一步分析其实验结果。 6台水下摄像机(图2)的视频显示,所有产卵活动(请参见补充视频“繁殖鱼和产卵”)都发生在弯道部分。


因此,从CFD模拟流场中提取了这些弯道处的流速,以进行进一步分析。实验期间水槽中的水质条件示于表S4中。


图5 不同实验情况下大环形水槽的模拟流场<br label=图片备注 class=text-img-note>
图5 不同实验情况下大环形水槽的模拟流场


对提取的各实验场景流速和产卵雌性比例进行统计分析[图6(a)],发现该比例随着流速的增加先显著增加,直到1.4 m·s–1,然后在1.6 m·s–1时急剧下降。在1.4 m·s–1流速下观察到产卵雌性的比例最高,达到 67% [图6(a)]


但是,鱼卵的最大受精率出现在流速为1.2 m·s–1的情况下[图6(b)],并且显著高于流速为 1.0 m·s–1和1.6 m·s–1情况。受精率在1.2~1.4 m·s–1之间时数据没有显著差异。


随着流速的增加,家鱼胚胎的孵化率(孵化的卵/ 实验卵总数)下降[图6(c)]。在1.2 m·s–1的速度下,孵化率低至15%。幼鱼在24 h内的存活率(存活的幼鱼/总实验卵)显示出与孵化率相似的趋势。


在静态水中的存活率约为12%,在1.2 m·s–1时,存活率降至4%。孵化率和存活率低主要是由于从大水槽收集的实验卵的低受精率,因为产卵实验期间大水槽中的高流速导致了低受精率。


图6 在实验室实验中,不同流速下家鱼的繁殖。(a)产卵雌鱼与流速的比率(n = 4);(b)受精率与流速的关系(n= 4);(c)受精卵的孵化率和幼鱼的存活率(n = 3)。使用一维高斯函数拟合曲线<br label=图片备注 class=text-img-note>
图6 在实验室实验中,不同流速下家鱼的繁殖。(a)产卵雌鱼与流速的比率(n = 4);(b)受精率与流速的关系(n= 4);(c)受精卵的孵化率和幼鱼的存活率(n = 3)。使用一维高斯函数拟合曲线


3. 室内实验和野外试验一致性分析


在水槽实验中,最高鱼卵密度出现在1.4 m·s–1的流速下,最低鱼卵密度出现在1.6 m·s–1的流速下[图7 (a)]野外试验产卵事件数据的分析表明,速度范围为 1.05~1.60 m·s–1适合家鱼在天然水体中产卵[图7(b)]


采用Vismara等的方法对室内和野外试验数据进行归一化处理,使用高斯回归拟合归一化的产卵数据(S)和流速(v),发现对于水槽实验,S = exp [–(v – 1.3062)2 / 0.0707](R2 = 0.9289,P < 0.0001),野外试验S = exp [–(v – 1.3143)2 / 0.0968](R2 = 0.9828,P < 0.0001)(图7),实验室和野外试验的结果高度一致,可以归一化写为S = exp [–(v – 1.31)2 / 0.097]。因此,确定在河流中家鱼产卵的最佳流速为1.31 m·s–1


图7 鱼类繁殖与流速的关系。在水槽实验中不同速度水平下的归一化鱼卵密度及其高斯拟合曲线(a);天然产卵场产卵活动的频率直方图和高斯拟合曲线(b)<br label=图片备注 class=text-img-note>
图7 鱼类繁殖与流速的关系。在水槽实验中不同速度水平下的归一化鱼卵密度及其高斯拟合曲线(a);天然产卵场产卵活动的频率直方图和高斯拟合曲线(b)


4. 研究结果的工程试验性应用


研究结果在指导三峡水库适应性运行方面产生了显著的效果。根据最低速度限值和逐日上升速度、最大速度限值以及宜都四大家鱼产卵场的断面几何形状,在产卵事件开始的第一天,建议流量从12 500 m3 ·s–1(流速为1.05 m·s–1开始,在产卵达到高峰时的第4天逐渐增加到18 600 m3 ·s–1(流速为1.40 m·s–1,然后迅速下降以支持孵化和仔鱼生存[图8(a)和附录A中的第4节]


在家鱼产卵季节(三峡水库下游的5月中旬至6月中旬),重复上述流量过程。在2017—2019年的测试期间,按照与野外试验中采集鱼卵相同的程序监测,以评估效率。结果表明,建议的三峡水库生态调度在增加家鱼产卵方面取得了持续的显著效果[图8(b)]


图8 在三峡大坝下游的作业中应用流速工程的建议流程,以改善家鱼的繁殖(a);自流程实施以来,三峡大坝下游家鱼产卵的变化(b)。数据来自中国三峡公司和中国环境监测中心<br>
图8 在三峡大坝下游的作业中应用流速工程的建议流程,以改善家鱼的繁殖(a);自流程实施以来,三峡大坝下游家鱼产卵的变化(b)。数据来自中国三峡公司和中国环境监测中心


四、讨论


通过野外试验和实验室控制实验,量化了流速对家鱼(一种产漂流性鱼卵的物种)产卵和胚胎发育的影响。研究结果直接表明一定的流速是鲢产卵的先决条件。在水槽实验中,成功产卵亲鱼的百分比随着流速的增加而增加,并在流速为1.4 m·s–1时达到最大值(图3)


但是,当流速高达1.6 m·s–1时,产卵行为将受到严重抑制,这是由于家鱼一直处于游泳状态、体力消耗过多所致。对于长江自然产卵场,在1.05~1.60 m·s–1的较宽速度范围内均观察到了鲢鱼卵。天然水体中适宜产卵的流速范围较宽,归因于复杂的水文条件和生境特征。除了家鱼产卵的适宜流速外,触发流速是另一个重要指标。


在 0.8 m·s–1的速度下,水槽中未发现产卵行为,确定的触发速度约为1.0 m·s–1(图4)。在天然河段,产卵速率为 0.90~1.90 m·s–1(图4)。不同鱼类产卵所需的流速存在巨大差异,雌性大西洋鲑鱼在平均流速为0.53 m·s–1的区域产卵,而在流水条件下,青鳉鱼产卵会减少,当水流静止时又恢复产卵。本研究表明,尽管鲢多在相对静态的水域中生活和生长,但其繁殖仍需要一定的水流速度。


流速不仅触发了家鱼的产卵行为,而且还影响了鱼卵的受精率。在1.2 m·s–1观察到最高的受精率,但在 1.2~1.4 m·s–1之间受精率有显著差异(图3)。由于鲢是体外受精,精卵融合通过水流混合在水中发生。


由于精子的游动能力在时间和空间上都受到限制,因此高于1.6 m·s–1的流速会导致精子和卵子相遇失败,从而使受精率迅速下降。由于鲢鱼卵是漂浮的,并且在生命的早期阶段必须保持悬浮状态,因此在其鱼卵受精后需要在一定的流速下孵化。


Garcia等研究了鲢鱼卵在不同流速下的沉降,发现0.2 m·s–1可以支持其鱼卵的悬浮。然而,在本研究中,随着流速的增加,鲢鱼卵的孵化率和存活率均显著下降。结果表明,静态水域更有利于鱼卵早期的生长发育。


根据Murphy和Jackson的研究结果,低至0.15~0.25 m·s–1的流速足以让鲢鱼卵悬浮在大湖区的四个支流中。鱼类生命的早期阶段,特别是在快速成长的胚胎和幼鱼阶段,最适合感应环境压力,因为胚胎和幼鱼比成年鱼更容易受到环境压力的影响。


环境因素的细微变化可能会导致表型发生重大变化,从而对后期个体发育和生存产生重大影响。过高的流速会导致鱼卵的机械损伤,这可能是鱼卵孵化和仔鱼存活同时受到抑制的主要原因。受精后,鱼卵向下游漂流,胚胎发育在相对缓慢的水流中进行。


虽然较高流速可以刺激鱼类产卵,但由于在快速水流下幼鱼的存活率较低,因此在较高流速下鱼类的数量并不会增加。水流除了维持鱼卵悬浮的主要功能外,在较高的流速下水交换和溶解氧的含量也会增加,这有利于鱼卵孵化。在我们的水槽实验中,水深较浅,鱼卵的密度也较低,因此氧气条件不是限制其早期发育的因素。


众所周知,硬骨鱼的繁殖性能受垂体激素,尤其是促性腺激素(GtH)的调节,而垂体激素受促性腺激素释放(GnRH)的控制。人们认为环境条件,如增加流速、长光照周期和暖和的温度是控制硬骨鱼繁殖的重要因素。


通过刺激性腺活性,观察到促性腺激素释放激素的活性上调。据报道,流速上升是许多硬骨鱼(包括鲢)产卵的关键。与普通家鱼(如锦鲤)不同,由于成熟的雌性家鱼卵巢发育停滞,鲢无法在静态水体中自然产卵,人工注射促黄体生成激素释放激素II 已被证明可有效诱导鲢产卵。


本实验中促黄体生成激素的应用证明了产卵过程中流速的重要性。流速为 0.8 m·s–1时,即使注入促黄体生成激素也不会引起鲢产卵活动,只有当流速达到1.0 m·s–1时才发生产卵。


在中国,长江是鲢的主要栖息地。据报道,长江上游水坝给家鱼繁殖造成了较大影响。新产卵场的形成将极大影响种群维持。本研究提供了家鱼产卵所需的明确流速范围,为大坝的合理运行以及家鱼保护提供了重要依据。


与长江家鱼种群规模减少相反,家鱼作为入侵物种正在许多国家引起严重的生态灾难。鲢已被引入至少88个国家或地区进行水产养殖、控制藻华或作为入侵物种,其中有超过三分之一已建立了自我维持的种群。鲢相对于本地物种的竞争优势以及对浮游生物群落的深远影响已造成严重后果。


据报道,2008年伊利诺伊河拉格兰奇河段鲢的生物量占鱼类总捕捞量的 51%。繁殖性状的高可塑性是成功建立新种群的原因。同时,一系列物理和非物理措施,包括水坝隔离、电流屏障、声屏障或生物学方法,使用锦鲤疱疹病毒、绝育或致死雌性个体,都已被用于抑制鲢大量繁殖。然而这些措施未能成功解决问题,甚至更糟的是,一些措施引发了意想不到的生态风险。


由于产卵和生命早期发育在入侵物种种群建立中具有重要作用,因此控制产卵过程对控制入侵至关重要。鲢产卵对流速的要求很严格,但在生命史其余阶段则要灵活得多。因此,通过水坝控制流速可能是抑制家鱼入侵的一种经济有效的方式。


考虑到世界上多数大河都建有水坝,这些发现为繁殖季节调整水库运行以操控流速来提升(1.2 <v <1.4 m·s–1或抑制(v <0.8 m·s–1或v >1.8 m·s–1家鱼种群提供了坚实的基础。


五、结论


近期研究过度强调了水坝建设的负面影响,特别是对鱼类的影响。家鱼在中国粮食安全和蛋白质获取中发挥着重要作用;然而,由于大坝造成水文情势发生显著变化,因此家鱼在长江的自然种群急剧缩小。流速是鱼类早期发育的关键因素,对鱼类种群起决定性作用。


本研究首次通过实验室和野外试验,同时考虑了产卵和胚胎发育需求,定量确定了家鱼繁殖的适宜流速范围。该发现提供了家鱼自然种群管控的非工程方法,如通过水坝调节河流流量。未来,研究将集中于考虑水文不确定性下梯级水库生态多目标优化调度方面。


改编原文:

Qiuwen Chen, Jianyun Zhang, Yuchen Chen, Kangle Mo, Jun Wang, Lei Tang, Yuqing Lin, Lei Chen, Yong Gao, Wei Jiang, Yuxin Zhang.Inducing Flow Velocities to Manage Fish Reproduction in Regulated Rivers[J].Engineering,2021,7(2):178-186.


本文来自微信公众号:中国工程院院刊(ID:CAE-Engineering),作者:陈求稳, 张建云, 陈宇琛, 莫康乐, 王骏, 唐磊, 林育青, 陈磊, 高勇, 姜伟, 张与馨