20162018年长江下游安庆江段四大家鱼仔稚鱼资源调查分析*

2020-07-17丁隆强何晓辉李新丰方弟安徐东坡

湖泊科学 2020年4期

丁隆强，何晓辉，李新丰，方弟安，，徐东坡，

(1:上海海洋大学水产科学国家级实验教学示范中心，上海 201306)

(2:中国水产科学研究院淡水渔业研究中心，农业农村部长江下游渔业资源环境科学观测实验站，无锡 214081)

青鱼(Mylopharyngodonpiceus)、草鱼(Ctenopharynodonidellus)、鲢(Hypophthalmichthysmolitrix)和鳙(Aristichthysnobilis)统称为四大家鱼，是我国特有的经济鱼类，其产量曾占我国淡水鱼总产量的60%[1]. 长江干流是四大家鱼主要的繁殖场所，作为典型的产漂流性卵鱼类，其繁殖行为及早期生长与长江生态水文条件相适应. 水利工程的建设改变了长江原有的生态环境，使得现有的满足四大家鱼繁殖需求的江段减少，同时受过度捕捞等因素干扰，家鱼亲鱼群体破坏较为严重，致使长江流域家鱼数量逐渐减少[2]. 近年来长江中游监利江段和汉江丹江口水库上游早期资源调查结果均显示四大家鱼资源量显著下降，2007年监利江段四大家鱼卵苗径流量仅为1981年的1.34%，2014年丹江口水库上游家鱼卵苗径流量下降到1977年的0.25%[3-4]. 长江干流四大家鱼产卵场多集中在中上游江段，历史调查资料表明下游仅在湖口和彭泽江段发现四大家鱼产卵场[5]. 任鹏的研究表明长江下游由于江湖连通效应及支流的汇入，距离大坝越远的江段其水文形势越趋于自然状态，使得四大家鱼仔鱼资源量自中游向下游安庆江段有增多的趋势[6].

安庆江段沿岸有众多附属湖泊，分布有长江江豚自然生态保护区、长江刀鲚国家级水产种质资源保护区等重要生态功能区，是重要的渔业资源水域[7]. 位于上游的湖口江段与鄱阳湖形成特殊的江湖复合型生态系统，对四大家鱼的繁殖及育幼提供良好的生境，使得距离鄱阳湖较近的安庆江段家鱼仔鱼丰度高于南京和靖江江段[8]. 安庆市是作为长江三角洲重要的水运枢纽和化工城市，频繁的航运、化工污染和过度捕捞等因素导致安庆江段渔业资源日益衰退. 2015年春季安庆江段成鱼调查中四大家鱼仅有鲢和草鱼为常见种，且鲢个体呈现小型化现象，四大家鱼资源明显衰退[9]. 本研究以长江下游安庆段的家鱼仔稚鱼(初孵仔鱼至鳞片形成期[10])为研究对象，在繁殖高峰期对安庆江段进行持续监测，了解安庆江段四大家鱼仔稚鱼资源量现状，运用GAM模型解析该江段家鱼仔稚鱼丰度变化与水文环境因子之间的关系，为评估长江下游四大家鱼资源补充量、苗汛预测和保护提供基础数据.

1 材料与方法

1.1 采样时间和地点

通过查阅长江鱼类早期资源的历史文献[10]及现场勘察，2016-2018年4-8月(2016年：4月17日-8月15日；2017年：4-8月逐月采样一周；2018年：4月22日-8月10日，共计239天)于长江下游安庆段开展早期资源调查. 采样断面设置于长江干流与皖河口交汇处上游50 m处(图1)，断面设置北岸(30°29′5″N， 116°59′24″E)、南岸(30°28′48″N， 116°59′59″E)和江心(30°28′56″N， 116°59′40″E)3个采样点，南北两岸采样点距离岸边10 m左右.

图1 安庆江段鱼类早期资源采样点设置

1.2 调查方法和样品的保存鉴定

四大家鱼早期资源采集使用大型圆锥状浮游生物网进行主动采集. 网具为直径0.8 m、网深2.5 m、网目40目、孔径0.25 mm的圆锥浮游生物网，网后连接圆柱形鱼苗收集装置，网口固定流量计用于测量网口过水量，采样时使用机动船对江水上层进行拖拽. 调查期间，每日8：00-11：00之间对各采样点进行定量采集一次，根据采集到的鱼类早期资源总捕捞数确定每次采集时间(3～10 min)，采样同时记录每次拖网持续时间. 使用Dretec电子水温计测量水温，采用萨氏透明度盘测量透明度，使用德国HYDRO-BIOS公司数字网口流量计记录流量计转数来换算网口过水流量. 采样断面的水位和江水径流量数据来分别自安徽省水文信息网的安庆水文站和大通水文站(由于采样断面和铜陵市大通镇之间无大型支流汇入长江干流，故可以反映该断面的江水径流量情况).

主动网具采集的家鱼仔稚鱼样品由冰水保存，样品鉴定完成后使用75%中性酒精保存. 仔稚鱼种类鉴定按照曹文宣等[10]编写的《长江鱼类早期资源》，通过形态学方法在奥林巴斯解剖镜SZX16进行种类鉴定，并分类计数. 对于形态破损难以鉴定确认的物种使用分子方法鉴定，具体步骤是：提取样本的基因组DNA，基于COⅠ(Cytochorme C oxidase Subunit Ⅰ，线粒体细胞色素C氧化酶亚基Ⅰ)基因条形码，采用PCR扩增目的条带并测序，将所测结果与已经用成鱼建立的长江鱼类条形码数据库进行比对得到种类鉴定的结果[11].

1.3 数据处理与分析

1)仔稚鱼密度计算方法参考易伯鲁等[12]的计算方法，依据所采集四大家鱼仔稚鱼的数量、采集时间、流量计的始末差值和主动网口面积按下列公式处理：

qi=0.3Ci·a

(1)

Di=Ni/qit

(2)

式中，qi为第i次采集网具网口的过水流量(m3/s)，Ci为第i次采集流量计的流量差，a为主动网具网口面积(m2)，0.3为流量计的常数，t为每次采集仔稚鱼时间(s)，Di为第i次采集仔稚鱼的密度(ind./m3)，Ni为第i次采集仔稚鱼数量(ind.).

2)采样断面鱼苗平均密度与定点鱼苗密度相比系数：

(3)

3) 采集期间卵苗径流量：

Mi=Di·Qi·C

(4)

式中，Mi为第i次采集时段内通过该江段断面的鱼苗数(ind.)，Qi为第i次采集时的断面流量(m3/s).

非采集期间鱼苗径流量：

Mi,j+1=(Mi/ti+Mi+1/ti+1)ti,j+1/2

(5)

用相邻两次采集的径流量及其间隔时间进行插补计算. 式中：Mi，i+1为第i，i+1采集时间间隔内的鱼苗径流量(ind.)，ti，i+1为第i，i+1次采集时间间隔(s).

采集江段的鱼苗总径流量为：

M=∑Mi+∑Mi,j+1

(6)

依据刘建康对长江下游四大家鱼鱼苗垂直分布的研究结果，江水上层家鱼苗平均捕获率占总捕获率的43.5%[13]，估算本研究中繁殖期内通过安庆江段的四大家鱼仔稚鱼总径流量.

4)广义可加模型(generalized additive models，GAM)是一种非参数化或半参数化的回归分析方法，能很好地模拟应变量和一个或多个预测变量之间的关系[14]. 本研究运用GAM模型分析水文因子对家鱼仔稚鱼丰度变化的因子进行分析，其表达式为：

(7)

式中，g(ρ)为关联函数；ρ为仔稚鱼丰度(ind./m3)；α为适合函数中的截距；xj为解释变量，包括水温、透明度、水位上涨率和径流量上涨率；fi(xj)为解释变量关系的非参数函数；ε为误差项，与解释变量xj无关，E(ε)=0,ε=δ2模型采用样条平滑法，对数据图进行平滑处理.

数据处理使用Excel 2016软件，使用Origin 2018软件制作直方图等分析图. 数据分析使用SPSS 19.0软件，差异性分析使用单因素方差分析，相关性分析使用Pearson相关性分析. GAM模型的构建及其检验在R软件(R 3.5.3)的mgcv包内完成.

2 结果与分析

2.1 安庆江段四大家鱼仔稚鱼种类组成及径流量

2016-2018年共调查239天，采集家鱼仔稚鱼3178尾，其中以鲢数量最多，共2160尾，占家鱼仔稚鱼总捕捞数量的67.97%；其次是草鱼，共计732尾占23.03%；鳙279尾占8.78%；青鱼数量最少，仅采集到7尾，所占比例为0.22%. 2016-2018年四大家鱼仔稚鱼变化规律呈现出鲢所占比例逐年上升，草鱼比例不断下降，鳙所占比例较为稳定(表1). 根据采集到的四大家鱼仔稚鱼数量及断面流量估算，2016-2018年4-8月份通过安庆江段四大家鱼仔稚鱼径流量依次为21.70×108、14.62×108、12.05×108尾，占各年份仔稚鱼总径流量的2.43%、0.99%和0.65%.

表1 安庆江段四大家鱼仔稚鱼种类组成

2.2 四大家鱼仔稚鱼丰度的时空分布特征

2016年安庆江段家鱼仔稚鱼日平均丰度为8.57±9.03 ind./100 m3. 仔稚鱼自4月中旬呈波动性上升趋势至7月中旬达到峰值，集中出现在5月下旬至7月中旬，共出现3次苗汛，分别为5月21日、6月15日、7月4日，对应的仔稚鱼密度为26.50、37.83、43.23 ind./100 m3. 2017年家鱼日平均丰度为6.76±4.29 ind./100 m3，采集时间内最高丰度为13.77 ind./100 m3. 对应日期为7月8日. 2018年日平均丰度为7.49±12.52 ind./100 m3，仔稚鱼丰度随江段径流量变化呈现先升后降再上升的趋势，在5月27日和7月18日达到峰值，对应的密度分别为13.74和55.86 ind./100 m3(图2).

对断面北岸、南岸和江心3个采样点的家鱼密度数据进行显著性分析，验证家鱼密度随断面横向分布的差异性. Kruskal-Wallis检验表明，断面3个采样点的家鱼仔稚鱼漂流密度均存在显著性差异(2016年：P=0.026<0.05，n=57；2017年：P=0.041<0.05，n=5；2018年：P=0.046<0.05，n=49). 2016-2018年断面家鱼密度横向分布均呈现出北岸>南岸>江心，三年内北岸平均漂流密度依次为14.14、17.83、16.30 ind./100 m3；南岸平均漂流密度依次为7.63、4.79、3.66 ind./100 m3；江心平均漂流密度依次为3.97、3.81、2.53 ind./100 m3.

图2 安庆江段四大家鱼仔稚鱼丰度变化

2.3 四大家鱼仔稚鱼丰度与水文环境因子之间的关系

调查期间内，安庆江段水温变化范围为18.0～29.8℃，平均值为25.30±2.44℃；透明度为16～70 cm，平均值为46.78±15.64 cm；溶解氧浓度为6.17～8.07 mg/L，平均值为7.24±2.15 mg/L；pH 7.39～8.17，平均值为7.80±0.20；江水径流量20400～70200 m3/s，日上涨率在100～7000 m3/(s·d)；水位8.44～17.69 m，日上涨率在0.01～1.25 m/d之间. 通过对家鱼仔稚鱼丰度与水文环境因子进行的相关性分析发现，家鱼日丰度变化与水温呈显著正相关(P<0.05)，与透明度呈极显著负相关(P<0.01)，与江水径流量、水位均呈极显著正相关(P<0.01).

GAM模型构建过程中逐步加入各个水文因子，依据AIC准则来选择最优模型. 模型最终结果见表2，其总偏差解释率为69.8%，预测模型的校正决定系数R-sq为0.622>0.5，说明模型较高程度上揭示了水文因子对家鱼仔稚鱼丰度的影响.

应用正态Q-Q图和Shapiro-wilk检验验证模型残差的正态性，诊断GAM模型. 正态Q-Q图表明GAM模型基本符合正态分布(图3). Shapiro-wilk检验表明，模型残差符合正态分布(w=0.941，P=0.070>0.05).

表2 安庆江段四大家鱼仔稚鱼GAM模型方差分析

图3 四大家鱼仔稚鱼丰度的GAM模型Q-Q图

图4 四大家鱼仔稚鱼丰度与各因子之间关系的GAM分析

GAM模型结果显示，根据偏差解释率，各因子相对重要性依次为水位日上涨率>透明度>水温>径流量日上涨率. 水文因子与家鱼仔稚鱼丰度之间关系的分析结果图表明(图4)，当水温在19.0～27.3℃时，家鱼仔稚鱼大量出现，其丰度随水温升高整体呈现上升趋势；在超过27.3℃后，其丰度逐渐下降. 透明度范围主要集中在24～57 cm之间，家鱼仔稚鱼密度与透明度呈现显著负相关性，即江水透明度下降时家鱼仔稚鱼密度上升. 水位和径流量的日上涨率与家鱼仔稚鱼丰度均呈现极显著正相关性，说明家鱼丰度随水位和径流量的日上涨率增大有所上升. 其中家鱼高峰期对应的水温、水位日上涨率和径流量日上涨率分别为27.3℃、0.23 m/d和900 m3/(s·d).

3 讨论

3.1 安庆江段四大家鱼仔稚鱼资源现状

三峡大坝运行以来，长江中上游四大家鱼卵苗种类组成出现以滤食浮游生物的鲢鱼苗数量迅速增加，草鱼所占比例不断下降，鳙和青鱼处于较低水平[2，15]，本研究结果与中上游相类似. 长江下游安庆江段2016-2018年四大家鱼组成中鲢以67.97%的比例占绝对优势，且所占比例有一定上升趋势；其次是草鱼，呈现出逐年下降的趋势；鳙和青鱼所占比例较低. 与2013年[6]结果相比青鱼资源量占比大幅度下降，可能原因是由于长江下游悬移质粒径较小，泥沙和有机质等营养物质在水流搬运过程中不易沉积，不利于大型底栖动物的生存和建群[16]. 作为底栖生物食性的青鱼，饵料的匮乏直接影响其生长及育幼过程中正常的摄食活动，使得青鱼早期资源量得不到有效补充. 长江干流中青鱼资源量的减少应引起重视.

由于沿江支流的汇入及水流作用，距离大坝较远的江段其水文情势越接近自然状态，大坝对鱼类资源的不利影响就会越小[17]. 2016年安庆江段四大家鱼仔稚鱼日平均密度为8.57 ind./100 m3，显著高于同期中游黄石江段的卵苗日均密度[18]，表明四大家鱼资源量自长江中游部分江段向下游安庆江段有增多的趋势. 长江下游距离三峡大坝较远，水文特征受影响较小，同时鄱阳湖与长江的连通关系不仅有利于家鱼亲鱼出湖繁殖，也有利于仔鱼入湖育幼，鄱阳湖和长江形成的江湖复合生态系统可能是安庆江段家鱼资源量较大的原因.

近年来对长江流域早期资源的调查结果均呈现出四大家鱼资源量处于不断衰退的趋势. 2014-2015年刘明典等对长江中游早期资源调查结果显示大坝下泄的不饱和水流改变了原有的水文环境，使得宜昌江段四大家鱼产卵场规模呈减少趋势[19]. 受大量涉水工程和围湖造田的影响，安庆江段和沿江湖泊的连通关系均发生变化，导致江湖洄游性鱼类自然繁殖受到影响. 本研究中估算2016-2018年通过安庆江段的四大家鱼仔稚鱼径流量分别为21.70×108尾、14.62×108尾、12.05×108尾，呈现出逐年下降的趋势. 受江湖阻隔的影响，安庆江段鱼类群落结构表现出洄游性鱼类减少的趋势[20]，长江干流四大家鱼繁殖亲本数量的减少，使得仔稚鱼径流量逐年下降. 此外，作为典型的产漂流性卵鱼类，四大家鱼产卵规模同江水涨水过程密切相关[21]. 2018年江段径流量较往年相比上涨持续时间较短，上涨幅度较低. 涨水过程的缩减一定程度上影响四大家鱼的正常繁殖活动.

3.2 影响四大家鱼仔稚鱼时空分布差异的因素

鱼类早期资源的时间动态能直接反映鱼类繁殖时间的变化特征[22]. 安庆江段四大家鱼仔稚鱼集中出现在5-7月，高峰期多发生在6月下旬和7月，与长江中游四大家鱼繁殖时间大致相同. 2016-2018年安庆江段四大家鱼苗高峰期丰度有所上升，但其主汛期存在推迟且持续时间有缩短的迹象. 2016和2017年主汛期均处于6月下旬至7月上旬，2018年则推迟至7月中下旬. 可能是由于2018年5月份江水涨幅较低，不满足家鱼繁殖所需的水文条件，使得四大家鱼集中在7月中下旬江水上涨过程时进行繁殖活动.

河流鱼类早期阶段受非生物和生物因素影响会存在一定的漂流规律. 非生物因子包括洪水的水位涨幅变化、河床形态、光照强度、水温水文条件等，生物因子有鱼类的发育阶段、栖息地偏好及其种间关系等[23]. 徐薇等研究发现长江鱼类在早期阶段存在水平分布上的显著性差异，断面空间卵苗密度与流速存在负相关性关系[24]. 本文在分析采样断面四大家鱼仔稚鱼丰度空间分布时，发现江心丰度显著低于南北两岸. 江心水流量大，流速较两岸快，不利于仔稚鱼的栖息和生存，使得具有一定游泳和主动选择能力的仔鱼聚集在沿岸缓流回水区. 此外其分布还受生物等其他因素的影响，采样断面3个点位的年际变化呈现南岸丰度均值逐年接近江心，北岸家鱼仔稚鱼丰度有一定的增长趋势. 采样断面属于顺直型江段，两岸江水流速无较大差异，其南北两岸家鱼仔稚鱼丰度可能受生境差异等因素影响. 北岸拥有较高的植被覆盖率，相比南岸航运及农耕等人类活动所造成的污染可以提供更为丰富的生物饵料，有利于仔稚鱼早期的摄食活动.

3.3 水文环境因子对四大家鱼仔稚鱼丰度的影响

四大家鱼繁殖特性相似，其产卵所需的水文条件也相似，不同种类在相同产卵场进行繁殖活动，常被视作一个类群进行研究. 作为典型的产漂流性卵鱼类，亲鱼需要在特定的水温条件下，同时有适宜水位涨幅才能促使排卵[25]. 水温被认为是影响鱼类产卵的重要因素，适宜四大家鱼产卵的水温一般在18～24℃. 长江中游监利断面家鱼卵苗高峰期在6月下旬至7月初，水温范围在20～24℃[26]. 本研究中安庆江段同期水温在22.4～25.2℃，略高于监利江段. 推测水温的差异是导致家鱼资源量长江中游向下游有增多趋势的原因. GAM模型结果中安庆江段水温范围在19～27.3℃之间，随着水温的上升，家鱼仔稚鱼丰度整体呈上升趋势.

透明度随江水径流量的变化被认为是鱼类产卵起始的重要因素. 张晓敏等研究2004-2007年汉江中下游洪峰过程中发现，当水温和涨水条件均满足家鱼繁殖需求的情况下，仅在透明度低于14.5 cm时监测到家鱼卵[27]. 因此洪水过程中透明度的高低和家鱼自然繁殖存在一定的相关性. 本文结果表明家鱼仔稚鱼密度与透明度呈显著负相关关系，随着透明度下降，家鱼密度升高. 即在繁殖季节，透明度的降低提高了家鱼产卵的可能性，同时有利于仔稚鱼躲避敌害的能力.

此外，水流刺激也是四大家鱼产卵的必备条件. 易雨君等探究长江中游家鱼栖息地适宜度曲线的结果表明水位的变幅是家鱼繁殖活动的必备条件，鱼苗径流量与洪峰呈正相关关系，长江中游家鱼产卵的理想水位日上涨率为0.3 m/d[28]. 本研究发现安庆江段四大家鱼仔稚鱼丰度伴随着江水水位和径流量上涨有显著的上升趋势，证明水位和径流量的上涨是刺激四大家鱼繁殖的重要因子. GAM模型结果显示安庆江段水位日上涨率达到0.23 m/d或径流量日上涨率达到900 m3/(s·d)时，仔稚鱼丰度开始显著上升.

3.4 对长江下游四大家鱼资源保护的建议

长江下游城市经济发展迅速，人类活动频繁，为保护和恢复四大家鱼资源，基于以往和本文研究结果建议：(1)坚持贯彻落实长江流域禁渔制度，保证家鱼亲本的正常繁殖. (2)动态监测鱼类资源，以便制定和调整增殖保护对策. 针对下游鳙和青鱼数量稀少的情况，应通过开展科学增殖放流，增加繁殖群体数量，促进渔业资源恢复. (3)应注意保护洲滩及岸线的自然属性，尤其是实施航道整治和护岸等工程时. (4)严控各类污染物的排放，保证鱼类栖息地环境质量. (5)确保通江湖泊与长江干流的连通性，并逐步恢复其他沿江湖泊及支流与长江干流的连通性，充分发挥江湖复合生态系统的作用，以此来保护和恢复长江下游家鱼资源.