刘 熠,任 鹏,杨习文,李新丰,徐东坡**,方弟安**
(1:中国水产科学研究院淡水渔业研究中心, 农业部长江下游渔业资源环境科学观测实验站, 无锡 214081)(2:上海海洋大学水产科学国家级实验教学示范中心, 上海 201306)
刀鲚(Coilianasus)属鲱形目(Clupeiformes)、鳀科(Engraulidae)、鲚属(Coilia),俗称刀鱼,分布在中国、日本、朝鲜沿海和与之相通的河流、附属湖泊,其成鱼脂肪丰满,肉质鲜美,是长江中下游重要的经济鱼类之一[1-3]. 近年来,随着长江下游城市经济的迅速发展,长江下游干流江段沿岸密集分布着许多工厂建设,引起水质污染加剧,渔业方面的捕捞力度也随需求量增加而加大,种种原因导致水域生境受到的压力过大,渔业资源明显下降[4-6]. 受上述因素影响,刀鲚的产卵场和育幼场不断遭到破坏,洄游距离也逐渐缩短,1970s鄱阳湖分布着的大量洄游性刀鲚产卵场已经极少[7-8];长江口的刀鲚年产量历史最高3945 t,到2009年已经下滑至年均产量仅为86.2 t;长江定居型刀鲚的资源也不容乐观,种群的小型化现象较为严重[7,9-10]. 产卵场和育幼场的破坏是导致刀鲚补充群体减少的重要原因,鱼类栖息地是鱼类产卵和育幼的主要场所,为鱼类种群提供了必要生存基础[11]. 因此,对刀鲚栖息地的保护是展开刀鲚保护工作的首要方向,针对刀鲚栖息地的保护需对其早期资源的进行全面研究.
历代学者对刀鲚早期资源有较系统的探究,如在其资源动态监测、食性、刀鲚人工繁殖、年龄与生长等方面[12-18],但对长江下游各江段刀鲚仔稚鱼的研究较少. 长江下游江段是指湖口至东海入海口的一段长约940 km的干流江段,途径江西、安徽、江苏和上海等多个省市[19]. 其中湖口江段连通鄱阳湖,是鄱阳湖与长江下游进行生态联系的首要江段, 其水域生境复杂, 为多种洄游型鱼类提供了关键洄游通道[20];安庆和靖江江段均是已知刀鲚的重要繁殖场,靖江江段刀鲚产量非常高,其刀鲚年际产量能占长江下游总产量的28%[21]. 可见,系统地对长江下游刀鲚早期资源的时空格局进行探究,有利于对长江下游刀鲚栖息地进行多方面的保护. 为此,本研究利用鱼类早期资源调查方法,探究长江下游湖口、安庆和靖江江段刀鲚仔稚鱼的时空格局,对比不同江段刀鲚仔稚鱼的丰度情况,分析刀鲚仔稚鱼丰度的变化规律及其与主要环境因子的关系,以期为长江下游刀鲚资源的保护提供科学依据.
鱼类早期资源调查方法依据曹文宣所著《长江鱼类早期资源》[22]. 长江刀鲚的繁殖期从4月下旬开始,到11月份才结束,据历史调查结果,江西和湖南江段繁殖期在5-7月,安徽江段在4-8月,江苏江段在4-10月,崇明江段在4-11月,总的繁殖高峰期集中在5-8月[7,14]. 据此,本研究于2016年5月18日至8月15日(刀鲚繁殖高峰期内)在长江下游湖口、安庆、靖江3个江段同时开展刀鲚早期资源调查.
采样点设置如图(图1)所示,湖口、安庆、靖江江段各取一个断面,在断面南北两岸和江段中心各取一个位点,按顺序分别对其编号:湖口南岸(A1:29°47′12″N,116°16′30″E)、中心(A2:29°47′29″N,116°16′20″E)、北岸(A3:29°47′41″N,116°16′7″E);安庆南岸(B1:30°28′48″ N,116°59′59″E)、中心(B2:30°28′56″N,116°59′40″E)、北岸(B3:30°29′5″N,116°59′24″E);靖江南岸(C1:31°56′19″N,120°6′9″E)、中心(C2:31°56′52″N,120°6′30″E)、北岸(C3:31°57′26″N,120°6′40″E).
鱼类早期资源的采样使用主动网具,网具规格为0.5 mm网目的圆锥状浮游生物网,网长2 m, 网口面积为0.50 m2,锥形网尾部连有一个圆柱型的鱼苗收集装置,网口固定一个螺旋桨式的流量计用于测量过水流量,采样时使用机动船进行拖拽. 在调查期间, 每日8:00至12:00之间在每个断面设置的三个采样点进行定量采集,根据采集到总的鱼类早期资源数量(300~500尾)确定每次采集时间,采样同时记录每次拖网持续的时间.
图1 长江下游刀鲚早期资源采样点设置Fig.1 Sampling distributions on Coilia nasus resources of early life history stages in the lower reaches of the Yangtze River
环境因子的采集使用Dretec电子水温计测量水温,采用萨氏透明度盘测量透明度,使用德国HYDRO-BIOS公司数字网口流量计记录流量计转数来换算过水流量,在各个采样点的每次采样均进行环境因子数据的统计. 浮游生物网中所采集到的仔稚鱼样品使用5%的福尔马林溶液进行保存带回实验室在体式解剖镜(Olympus SZX 16)下依据曹文宣[22]的方法进行仔稚鱼种类的鉴定与统计计数,其中刀鲚卵黄囊期至脊索弯曲期的仔鱼鉴定主要基于刀鲚仔鱼的肌节数目,结合其细长的体态、较长的消化道来进行准确鉴定,发育至脊索弯曲后期后的刀鲚体态逐渐接近于成鱼,可基于其体态特征与鳍条位置进行鉴定. 所有鱼类早期资源样品鉴定完成后使用75%的中性酒精保存,将刀鲚仔稚鱼样品单独分开存留.
在一个江段的调查中,将所有点位均进行了鱼类早期资源采样的日期作为一次有效采样天数. 仔稚鱼密度计算方法参考曹文宣[22]和易伯鲁等[23]的计算方法,依据所采集仔稚鱼的数量、采集时间、流量计流量差、圆锥状浮游生物网网口面积按下列公式处理:
Qi=(0.3Ci·a)/t
(1)
Di=Ni/(Qi·t)
(2)
式中,Qi为第i次采集圆锥状浮游生物网网口的过水流量(m3/s),Ci为第i次采集流量计的流量差,a为主动网具网口面积(m2),0.3为流量计的常数,t为每次采集仔稚鱼的时(s),Di为第i次采集仔稚鱼的密度(ind./m3),Ni为第i次采集仔稚鱼数量(ind.).
使用PRIMER 5.0软件分析刀鲚仔稚鱼的时间动态. 将调查期间各江段不同月份作为样本, 刀鲚仔稚鱼密度为变量, 建立1个仔稚鱼密度的原始矩阵, 再根据Bray-Curtis相似性系数, 构建相似性三角矩阵. 根据不同月份的相似性指数, 使用组平均连接法进行等级聚类分析(Cluster analysis)和非参数多维尺度分析(non-Metric Multidimensional scaling, MDS), 并通过相似性分析(ANOSIM)分析对不同类群间差异显著性进行检验[24-27]. 检验 MDS分析结果的优劣用胁强系数(stress)来衡量,通常认为当 stress<0.2时,其图形有一定的解释意义;当stress<0.1时,可以认为是一个好的排序;当stress<0.05时,说明具有很好的代表性[28]. 数据处理使用Excel 2016软件,使用Origin 2018软件制作直方图、折线图、断层图等各种分析图. 数据分析使用SPSS 19.0软件处理, 差异性分析使用单因素方差分析,相关性分析采Pearson相关性检验.
由于鱼类早期资源调查是在野外环境下展开的生态学实验,容易受恶劣天气等环境影响导致不能保证每天有效采样. 在2016年5月18日至8月15日设定调查时间内,湖口江段受环境因素影响较小,有效采样时间为87天,共采集刀鲚仔稚鱼837尾;安庆江段少数时间受大风影响,有效采样时间为78天,共采集刀鲚仔稚鱼1239尾;靖江江段靠近长江口,江面较宽且易受潮汐影响,导致风浪较大不易过江,有效采样时间为40天,共采集刀鲚仔稚鱼1281尾.
调查发现,在采样期间湖口江段平均密度为1.56 ind./100 m3,安庆江段平均密度为10.46 ind./100 m3,靖江江段平均密度为29.94 ind./100 m3. 3个江段中刀鲚仔稚鱼密度在7月初之前的整体趋势均为在波动中逐渐上升,湖口江段在6月20日达到密度最高峰(50.37 ind./100 m3),随后安庆江段在7月1日达到密度最高峰(96.05 ind./100 m3),靖江江段在6月26日达到密度最高峰(320.02 ind./100 m3). 湖口江段的刀鲚仔稚鱼密度达到峰值后恢复至之前较低水平. 安庆江段和靖江江段达到峰值后则在波动中不断下降,其中安庆江段在7月12-20日一直保持较高水平,并在7月18日出现小高峰56.92(ind./100 m3);靖江江段在6月29日和7月26日出现了2次小高峰,峰值分别为172.01和39.98 ind./100 m3(图2a1、b1、c1).
刀鲚仔稚鱼密度变化在空间上的对比结果如图2的a2、b2、c2所示. 湖口江段密度在中心点达到的峰值时最大,为130.15 ind./100 m3;南岸点其次,为19.23 ind./100 m3;北岸点最小,为6.77 ind./100 m3. 对湖口江段3个采样点进行差异性分析结果表明,南北两岸刀鲚仔稚鱼密度存在极显著性差异(n=88,P<0.01). 安庆江段密度在南岸点达到的峰值最大,为173.96 ind./100 m3;北岸点其次,为113.01 ind./100 m3;中心点最小,为42.26 ind./100 m3. 对安庆江段3个采样点进行差异性分析结果表明,南北两岸刀鲚仔稚鱼密度差异不显著性,南北两岸和中心点均存在显著性差异(n=78,P<0.05). 靖江江段密度在南岸点达到的峰值最大,为930.24 ind./100 m3;北岸点其次,为318.03 ind./100 m3;中心点最小,为101.54 ind./100 m3. 对靖江江段3个采样点进行差异性分析结果表明,南岸和中心点存在显著性差异(n=40,P<0.05).
调查期间,湖口江段平均水温为25.76℃,水温变化范围为21.47~30.20℃,变化幅度达到8.73℃;安庆江段平均水温为25.78℃,水温变化范围为21.07~30.00℃,变化幅度达到8.93℃;靖江江段平均水温为26.02℃,水温变化范围为20.87~30.03℃,变化幅度达到9.16℃. 3个江段在平均水温和温差变化等方面相差不大,且均在6月中旬至7月中旬出现小高峰后回调了一段时间,各江段刀鲚仔稚鱼的密度最高峰均出现在这个时期(图3a1、b1、c1). 对各点水温和刀鲚仔稚鱼密度的相关性分析均显示无显著性相关(P>0.05).
通过对各江段透明度的统计,湖口江段平均透明度为42.03 cm,最高峰值为55.67 cm,7月2日前透明度一直在34.7~55.7 cm之间波动,7月4日透明度下降到整个调查期间最低值(22.6 cm),然后缓慢回升. 安庆和靖江江段也有一个透明度骤然下降的过程,安庆平均透明度为29.15 cm,最高峰值为41.33 cm,在7月9日降至最低(16 cm);靖江平均透明度为21.48 cm,最高峰值为30.33 cm,在7月8日降至最低10 cm后再波动回升(图3a2、b2、c2). 对各点透明度和刀鲚仔稚鱼密度的相关性分析均显示无显著性相关(P>0.05).
通过对各江段过水流量的统计,湖口江段呈波动性上升状态,刀鲚平均密度最高时过水流量为0.52 m3/s,调查期间平均过水流量为0.44 m3/s;安庆江段呈波动性下降趋势,刀鲚平均密度最高时过水流量为0.64 m3/s,调查期间平均过水流量为0.65 m3/s;靖江江段过水流量整体呈波动性下降趋势,平均密度最高时过水流量为0.40 m3/s,调查期间平均过水流量为0.35 m3/s(图3a3、b3、c3). 对安庆江段过水流量和刀鲚仔稚鱼密度的相关性分析显示两者显著相关(P<0.05),其它2个江段无无显著性相关(P>0.05).
对5-8月的刀鲚仔稚鱼进行月份尺度上的空间分布分析. 结果表明(图4),5月份安庆江段密度分布较为平均,其中北岸点密度较高,为0.40 ind./100 m3;湖口与靖江江段未采集到刀鲚仔稚鱼. 6月份湖口江段密度主要集中在中心点,为6.26 ind./100 m3;安庆江段南北两岸密度相差不大,分别为18.58和19.21 ind./100 m3;靖江江段刀鲚仔稚鱼主要集中分布在南岸点,为119.53 ind./100 m3,且各采样点的密度均高于其它江段. 7月份湖口江段刀鲚仔稚鱼主要分布在南岸点,为2.47 ind./100 m3;安庆江段主要分布在南北两岸,分别为21.14和17.79 ind./100 m3,中心点仔稚鱼密度最低;靖江江段仔稚鱼主要集中分布在南岸点,密度为59.17 ind./100 m3. 8月份湖口江段刀鲚仔稚鱼集中在中心点,为1.10 ind./100 m3;安庆江段主要分布在南岸点,为0.36 ind./100 m3;靖江江段主要集中分布在南北两岸,分别为6.03和5.13 ind./100 m3,且各采样点密度均高于其它江段.
图2 湖口、安庆和靖江江段刀鲚仔稚鱼密度的逐日变化(a1、b1和c1分别为湖口、安庆、靖江江段刀鲚仔稚鱼平均密度变化;a2、b2和c2分别为湖口、安庆、靖江江段刀鲚仔稚鱼在不同采样点的密度变化)Fig.2 The daily density change of larvae and juvenile C. nasus in Hukou, Anqing and Jingjiang sections (a1, b1 and c1 showed the average density of larvae and juvenile C. nasus in Hukou, Anqing and Jingjiang sections, respectively; a2, b2 and c2 showed the density of larvae and juvenile C. nasus in different sampling sites of Hukou, Anqing and Jingjiang sections, respectively)
图3 湖口、安庆和靖江江段刀鲚仔稚鱼密度与环境因子的变化关系(a1、b1和c1分别为湖口、安庆、靖江江段刀鲚仔稚鱼密度与水温的关系;a2、b2和c2分别为湖口、安庆、靖江江段刀鲚仔稚鱼密度与透明度的关系;a3、b3和c3分别为湖口、安庆、靖江江段刀鲚仔稚鱼密度与过水流量的关系)Fig.3 Changes between environmental factors with larvae and juvenile C. nasus density in Hukou, Anqing and Jingjiang sections(a1, b1 and c1 showed the relationships between C. nasus larvae and juvenile density and water temperature in the Hukou, Anqing and Jingjiang sections, respectively; a2, b2 and c2 showed the relationships between C. nasus larvae and juvenile density and transparency in the Hukou, Anqing and Jingjiang sections, respectively; a3, b3 and c3 showed the relationships between C. nasus larvae and juvenile density and water flow in the Hukou, Anqing and Jingjiang sections, respectively)
图4 5-8月长江下游各江段刀鲚仔稚鱼密度的空间分布Fig.4 The spatial distribution of density of larvae and juvenile Coilia nasus in each section of the lower reaches of Yangtze River from May to August
将6-8月各江段刀鲚仔稚鱼密度进行聚类分析,由各月份聚类分析树状图可以看出,以相似度22.01为界限可将9个月份分为两组(图5):安庆8月单独聚为一组,其余月份聚为一组. 其中,湖口6-8月、安庆6-7月、靖江6-8月以相似度44.47为界限可分为两个小组:湖口6月、7月和8月相似度较高,聚为一组;靖江6月、7月、8月和安庆6月、7月相似度较高,聚为一组.
基于各时期间Bray-Curtis相似性指数的9个月份刀鲚仔稚鱼密度MDS平面图(图6)的胁强系数为0.01 (Stress<0.05),能很好地解释各月份间的聚类结果, MDS图具有很好的代表性. 相似性分析(Analysis of similarity, ANOSIM)对组间差异显著性进行检验,结果显示R=0.85,P(level%) =11.1%, 说明两组间差异性不显著.
图5 6-8月各江段刀鲚仔稚鱼密度聚类分析Fig.5 The density cluster analysis of larvae and juvenile Coilia nasus in each section of the Yangtze River from June to August
图6 6-8月各江段刀鲚仔稚鱼密度非参数多维尺度分析Fig.6 The density in nonparametric multidimensional scaling analysis of larvae and juvenile Coilia nasus in each section of the Yangtze River from June to August
黎雨轩[14]在长江洄游性刀鲚繁殖生物学的研究中,对安庆、靖江和崇明3个江段刀鲚成熟个体的性体指数(GSI)和性腺发育分期进行了时空变化上的分析,发现刀鲚在4月份不进行繁殖,从长江下游往上游,5、6月份在河口附近进行繁殖,6、7月份在靖江繁殖,7、8月则在安庆繁殖. 王生等[29]在鄱阳湖湖口对刀鲚幼鱼出现时间的研究,推测出湖口江段洄游刀鲚成熟个体产卵时间大概在5-6月左右. 本研究发现,刀鲚仔稚鱼在湖口、安庆和靖江江段在6月下旬至7月上旬之间均出现一段高峰期,刀鲚密度的高峰期出现时间较短,这可能和其产卵类型有关,徐钢春等对刀鲚卵巢的研究发现刀鲚属于和鳗鲡等鱼类相似的同步型产卵类型,一年产一次卵,但一生中可多次繁殖活动[17, 30]. 本研究中刀鲚仔稚鱼大量出现在湖口和靖江江段主要发生在6月下旬,安庆江段发生在7月上旬,这和上述学者研究刀鲚成熟个体在这几个江段产卵的时期大致相同.
现在普遍观点认为长江中刀鲚具有两种生态类型:淡水定居型与溯河洄游型. 原来直接将短颌鲚划分为淡水定居型,长颌鲚划分为溯河洄游型,但近年徐钢春等[31]通过耳石微化学证明了短颌鲚也存在溯河洄游型个体,姜涛等[3]也通过对短颌鲚和长颌鲚矢耳石形态学的对比鉴定发现了长颌鲚中存在淡水定居型个体,传统的通过上颌骨长短来区分刀鲚生态类型的方法已经有一定的不确定性. 湖口江段靠近鄱阳湖,其鱼类组成中定居型刀鲚较多[20]. 从刀鲚洄游的角度分析各江段暴发期出现的时间,刀鲚自下而上洄游,游至湖口产卵的时间应在安庆之后,但本研究发现湖口刀鲚仔稚鱼密度高峰期出现在安庆江段之前,基于刀鲚洄游产卵的时期先后和湖口江段的鱼类组成结构,推断这种现象是由于湖口江段定居型刀鲚产卵造成的. 安庆江段衔接湖口江段,也是长江下游鱼类重要的洄游通道之一,鱼类仔稚鱼时期活动能力较弱,处于被动漂流阶段,一般顺水漂流而下,7月份安庆江段刀鲚仔稚鱼的大量出现可能有部分来自湖口江段[32].
各月份刀鲚仔稚鱼密度的空间分布和聚类分析结果显示,安庆和靖江江段具有一定的相似性,湖口江段则被分隔开来,刀鲚仔稚鱼密度总体在湖口江段较少,但8月份湖口密度高于安庆,表明8月份湖口水域产卵群体的数量较安庆多. 究其原因,这可能与湖口江段地理位置的特殊性对刀鲚发育的影响有关. 基于刘熠等[33]调查研究发现,2016年7月中下旬湖口江段发生了江水倒灌现象. 湖口江段连通长江和鄱阳湖,作为一个必要的环境单元结构参与形成江湖复合生态系统,而江水倒灌则是这个复杂生态系统的表现机制之一,它对许多鱼类的繁殖和补充起着重要的作用[34-35]. 李云良等[36]对长江倒灌鄱阳湖的水动力研究表明,贯穿整个湖区的主河道的流速是受倒灌影响最显著的区域. 徐钢春等[37]研究发现,刀鲚性腺的完全发育成熟需要一定程度的江水流速刺激,长江倒灌正好为鄱阳湖中的刀鲚性腺发育提供了相应条件,使得湖口水域刀鲚产卵群体增多.
刀鲚自然条件下产卵较喜好平缓的水流,少敌害而浮游动物丰富的水域,如长江沿岸水流较缓的通江支流和湖泊等[38]. 本研究中南岸、北岸和中心3个采样点刀鲚密度之间差异在各江段的表现不尽相同,湖口江段南北两岸差异性显著,安庆江段中心点与两岸差异性显著,靖江江段南岸与中心点差异性显著,这可能是由于沿岸地形和水体动力学差异所产生的现象. 湖口采样点处位于鄱阳湖与长江交汇的不远处,由于江湖水的汇合产生,中心地带的流速加快,造成南北两岸分布上的差异[39]. 有研究表明,洄游性刀鲚产卵群体通常在靠近江心的水域快速上溯寻找产卵场,而仔稚鱼和幼鱼则选择栖息在沿岸一带育肥和移动[40]. 安庆江段江面相对较窄,水域环境较为平和,因此,在环境较良好的情况下,刀鲚仔稚鱼更多地聚集在两岸一带. 靖江江段靠近长江口,常有沙洲将江水分流,加上潮汐的影响也易造成空间上的差异[21]. 地理环境的影响在对各月份刀鲚密度的空间分布结果中也有体现,湖口江段在大部分月份都集中在中心点,而安庆南北两岸在各月分布都较为均匀,靖江江段南岸在各月均保持优势.
刀鲚在自然条件下产卵,除了自身内部生殖腺需逐渐步入成熟阶段,体内囤积大量的脂肪为洄游和性腺发育提供大量能量外,达到繁殖条件还需要外界环境的刺激,如水温、盐度和水流量等因素[7]. 历史对刀鲚产卵场基本需求的研究得出,刀鲚喜好选择产卵水深约0.5~3.0 m,水流需要平缓,流速约0.057~0.075 m/s,水温在15.0~27.5℃,水质良好无污染且富有浮游动物的水域[38]. 本研究在湖口、安庆和靖江3个江段刀鲚仔稚鱼高峰期时的水流量均处于较低水平,温度分别为25.87、25.27和21.67℃,均达到了刀鲚产卵所需条件. 何为等[41]在南通和安庆江段对刀鲚性腺发育规律进行研究,发现性腺的成熟和水温及洄游时间密切相关;徐钢春等[42]通过培育的刀鲚性腺的抽样检查对自然状态下产卵的规律分析得到刀鲚产卵的最低和最高水温分别为20.0和30.5℃,本研究的3个江段温度变化范围为20.87~30.03℃,符合性腺在自然条件下在发育成熟的温度条件.
本研究在刀鲚繁殖期内对长江下游湖口、安庆和靖江江段的刀鲚仔稚鱼时空格局进行了系统的研究,结果发现,在调查期间,靖江江段刀鲚仔稚鱼平均密度最高;由于刀鲚生活史的同步型产卵,刀鲚仔稚鱼在3个江段均出现一段时间较短的暴发期,而湖口刀鲚仔稚鱼的暴发时间早于安庆和靖江江段的现象可能是由于湖口江段定居型刀鲚产卵造成;湖口江段8月产卵群体较多,推测可能是由于长江江水倒灌鄱阳湖的特殊环境现象刺激了刀鲚性腺的发育所致;各江段断面空间上的刀鲚仔稚鱼密度差异可能是沿岸地形和水体动力学差异所造成的. 本研究为探索资源变动规律与资源保护策略的关系提供了科学依据,但后期仍需进行许多更为细致的工作,如增加调查江段并对刀鲚早期不同发育阶段组成的时空变化进行研究等,今后将持续系统调查刀鲚早期资源的分布,尤其是产卵场的确定和特征分析,以期能够更加清晰地认识长江下游刀鲚仔稚鱼的时空变化特征.
致谢:感谢上海海洋大学的叶昆同学和刘思磊同学在野外鱼类早期资源采样中的帮助.