1. 引言
鲟鱼是大约2亿年前出现的鲟形目鱼类,隶属于硬骨鱼纲辐鳍亚纲、软骨硬鳞总目鲟形目、鲟科、鲟属,是现存最古老的脊椎动物之一,有活化石之称 [1] [2]。在二十世纪,由于偷猎、环境污染、栖息地退化和建造水坝等人类活动干扰鲟鱼产卵迁徙,许多鲟鱼物种都处于极度濒危状态 [3]。野生鲟鱼物种越来越少,全世界都在努力恢复鲟鱼的数量。目前,鲟鱼的繁殖主要依赖于人工繁殖,以补充野生种群,自2000年以来,中国已成为世界上最大的鲟鱼养殖区 [4] [5] [6] [7]。本研究以北京市鲟龙种业养殖场为实验基地,该养殖场孵化全国80%的鲟鱼苗种,对水产养殖业的贡献无可比拟,为国家水产养殖做出巨大贡献,这与养殖场池塘区位、水源、地形、沟渠走向、养殖品种等有关。虽然目前已有较多对鲟鱼的野生环境和仿生态环境的研究,但尚未见对养殖鲟鱼生长环境中水生生物群落的报道。目前,基于水生生物群落特征与水环境关系来阐述水体生态状况,指导水体生态修复的研究主要集中于江、河、湖、水库等淡水生态系统中,对于集约化养殖水体的研究工作,仍有待深入。
本研究对鲟鱼集约化养殖水体浮游生物群落分布格局进行系统研究,基于多元统计分析,阐述鲟鱼集约化养殖水体浮游生物群落格局与水体环境因子的相互关系,为我国鲟鱼集约化养殖水生态系统水生生物群落提供基于浮游生物群落的基础生态学依据。
2. 材料与方法
2.1. 采样点布设
2021年9月15日,在北京市鲟龙种业养殖场共设置5个浮游生物采样点,分别为进水口、A池、B池、C池、D池;并记录不同池塘大小、养殖种类、养殖数目状况。浮游生物采样点分布如图1所示。
(A池:长22.5 m、宽8 m、水深1.4 m、池深1.8 m,达氏鳇4尾、俄罗斯鲟40尾、施氏鲟50尾、西伯利亚鲟350尾;B池:长17.5 m、宽8 m、水深1.5 m、池深1.8 m,达氏鳇50尾、西伯利亚鲟160尾;C池:长20 m、宽8.5 m、水深1.4 m、池深2 m,施氏鲟300尾;D池:长5.5 m、宽4.5 m、水深1.4 m、池深1.8 m,施氏鲟50尾。)
Figure 1. Location of sampling sites
图1. 采样点分布图
2.2. 样品采集与处理
2.2.1. 水环境因子测定
在各采样点使用多参数水质分析测定仪(YSI ProPlus)现场测定水体温度(WT)、溶解氧(DO)、pH、氨氮(NH3-N)、叶绿素a (Chl-a)。使用DR1900便携式分光光度计测定水体总氮(TN)、总磷(TP)、化学需氧量(COD)。
2.2.2. 浮游植物的采集、鉴定与计数
浮游植物样品的定性采集用25号浮游植物网采集,在选定的采样点于水面0.5 m深处以每秒20~30 cm的速度作∞形循回缓慢拖动,时长5~10分钟左右,打开活塞,收集浓缩样0~40 ml。带回实验室在10 × 40倍光学显微镜下观察分类;浮游植物样品的定量采集利用5000 mL采水器,在水深小于3 m的采样点,在水面下约0.5 m处,停留片刻,采集水样;在水深超过3 m的采样点,除了在水面下0.5 m处采集水样外,还需在离水底0.5 m处,再采集一次水样,取得混合水样1000 mL,然后,向采集到的水样中加入15 mL鲁哥氏液,进行固定,经过48 h静置沉淀后,用3~5 mm的橡皮管,虹吸抽掉上清液,余20~25 ml沉淀物转到30 ml的定量瓶,用上清液清洗容器壁,并放入定量瓶中,定容至30 ml。然后摇匀取出0.1 ml水样,利用浮游植物计数框在10 × 40倍光学显微镜下观察计数,计数时观察100个视野,每个样品观察2片,2片计数之差不大于15%,即视为有效结果,取其平均值,最后换算成每升水样藻类的细胞密度(个/L) [8]。藻类生物量采用细胞体积法推算,单细胞藻类的生物量主要根据藻类个体形状测量而得 [9]。浮游植物的分类鉴定均依照《中国淡水藻类—系统、分类及生态》 [10]。
2.2.3. 浮游动物的采集、鉴定与计数
浮游动物定性采用13 #浮游生物网(孔径为0.112 mm)在水面下0.5 m处划“∞”型大约5~10 min,将滤出液放入样品瓶中并加入4%的甲醛溶液进行现场固定;定量采集使用1 L有机玻璃采水器在不同水层采集10 L的水样,并用25#浮游生物网过滤、收集到样品瓶中,加入4%的甲醛溶液进行现场固定,带回实验室静置沉淀48 h后,浓缩至30 mL的样品瓶中,参照《原生生物学》、《淡水浮游生物图谱》 [11] [12] 进行样品分类鉴定及生物量换算。
2.3. 数据分析与评价标准
浮游生物群落结构参数计算运用Margalef物种丰富度指数(D)、Shannon物种多样性指数(H')、Pielou均匀度指数(J) (公式1~4)及浮游生物优势度等参数分析浮游生物群落结构特征 [13] [14] [15]。优势度值Y ≥ 0.02的浮游生物定为优势种。各参数计算公式如下:
(1)
(2)
(3)
(4)
式中:S为浮游生物总种数;N为浮游生物个体总数;Ni为第i种浮游生物个体数;fi为第i种浮游生物各采样点出现频率。
利用SPSS20.0软件进行相关性分析;利用Canoco5软件包对浮游生物和环境数据进行冗余分析(RDA)。浮游生物分级评价标准见表1。
![](Images/Table_Tmp.jpg)
Table 1. The evaluation standard of plankton biodiversity assessment
表1. 浮游生物分级评价标准
3. 结果与分析
3.1. 浮游植物群落特征
3.1.1. 物种组成及优势种
本次调查进水口共采集到浮游植物2门16种(表2),其中绿藻门种类最多,有10种(占62.50%);其次为硅藻门,有6种(占37.50%)。
养殖池共采集到浮游植物5门61种(表3),其中硅藻门种类最多,有36种(占59.02%);其次为绿藻门,有17种(占27.86%);蓝藻门种类4种(占6.56%);再其次为裸藻门,物种数为3种(占4.91%);甲藻门最少,为1种(占1.64%)。硅藻门、绿藻门合计占鉴定总数86.88%,是秋季鲟鱼养殖池浮游植物的主要组成类群。如图2所示,A号池浮游植物种类数最多,为34种,而B号池物种数较少,为24种。根据优势度的计算结果,鲟鱼养殖池浮游植物优势种主要有11种,分别为硅藻门的中型脆杆藻(Y = 0.107)、谷皮菱形藻(Y = 0.085)、短小舟形藻(Y = 0.227)、短肋羽纹藻(Y = 0.208)、尖针杆藻(Y = 0.229)、湖北小环藻(Y = 0.068)、近缘桥弯藻(Y = 0.208)、普通等片藻(Y = 0.127);绿藻门的短鼓藻(Y = 0.046)、小球藻(Y = 0.052);裸藻门的旋转囊裸藻(Y = 0.033)。
![](Images/Table_Tmp.jpg)
Table 2. The list of phytoplankton in intake
表2. 进水口浮游植物名录
![](Images/Table_Tmp.jpg)
Table 3. The list of phytoplankton in sturgeon ponds
表3. 鲟鱼养殖池浮游植物名录
![](//html.hanspub.org/file/6-2830186x13_hanspub.png?20220323082705087)
Figure 2. The phytoplankton species composition at different sampling sites
图2. 各采样点浮游植物种类组成
3.1.2. 丰度和生物量
本次调查进水口浮游植物丰度为7.8 × 105 cells/L,游植物总的生物量1.564 mg/L。鲟鱼养殖池各采样点浮游植物丰度为13.2 × 105~21.36 × 105 cells/L,浮游植物平均丰度为16.86 × 105 cells/L;浮游植物总的生物量为18.284 mg/L,平均生物量为4.571 mg/L。各采样点浮游植物丰度和生物量有显著差异(p < 0.01),如图3所示,浮游植物丰度最高点位于D池,生物量最高点位于A池。从浮游植物相对丰度和相对生物量来看,均是硅藻门占绝对优势,其次是绿藻门占优势。
![](//html.hanspub.org/file/6-2830186x14_hanspub.png?20220323082705087)
Figure 3. Abundance and biomass of phytoplankton in water intake and aquaculture pond
图3. 进水口、养殖池各采样点浮游植物丰度和生物量分布
3.1.3. 浮游植物群落多样性
本次调查进水口浮游植物Margalef物种丰富度指数(D)为3.59、Shannon物种多样性指数(H')为3.46、Pielou均匀度指数(J)为1.25。鲟鱼养殖池各采样点浮游植物Margalef物种丰富度指数(D)平均值约为5.81,介于4.89~6.48,最大值同样位于D池,最小值位于C池。Shannon-Wiener 指数平均值约为4.53,各点无明显差异;Pielou均匀度指数(J)平均值约为1.34,各点无明显差异(如图4)。
![](//html.hanspub.org/file/6-2830186x15_hanspub.png?20220323082705087)
Figure 4. Distribution of diversity index of phytoplankton in water intake and aquaculture pond
图4. 进水口、养殖池各采样点浮游植物多样性指数分布
3.2. 浮游动物群落特征
3.2.1. 物种组成及优势种
本次调查进水口鉴定出浮游动物3大类11种(表4),其中轮虫8种,占浮游动物总物种数的72.7%,原生动物2种,占18.2%,桡足类1种,占9.1%。养殖池鉴定出浮游动物2大类39种(表5),其中轮虫23种,占浮游动物总物种数的58.9%,原生动物16种,占41.1%。根据优势度计算得出,鲟鱼养殖池浮游动物优势种主要有5种,分别为螺形龟甲轮虫(Y = 0.08)、狭颈宿轮虫(Y = 0.06)和原生动物中的白钟虫(Y = 0.04)、尾草履虫(Y = 0.03)、郝氏皱甲轮虫(Y = 0.02)。
![](Images/Table_Tmp.jpg)
Table 4. The list of zooplankton in intake
表4. 进水口浮游动物名录
![](Images/Table_Tmp.jpg)
Table 5. The list of zooplankton in sturgeon ponds
表5. 鲟鱼养殖池浮游动物名录
3.2.2. 丰度和生物量
本次调查进水口浮游动物丰度为10.72 ind./L,生物量为0.0078 mg/L。养殖池各采样点浮游动物丰度为9.38~22.11 ind./L,浮游植物平均丰度为13.23 cells/L;浮游动物生物量范围为0.0133~0.0347 mg/L,平均生物量为0.0227 mg/L。如图5所示,浮游动物丰度的峰值出现在A池,其他各点无明显差异;其生物量峰值也出现在A池,最低值出现在进水口。各采样点丰度和生物量大致呈现出丰度越高,生物量也越高的趋势。从浮游动物相对丰度来看,进水口、B池、C池主要是轮虫占优势,而A池和D池中轮虫和浮游动物持平。从浮游动物相对生物量来看,进水口是以桡足类占主要优势,其余样点均为轮虫占主要优势。
![](//html.hanspub.org/file/6-2830186x16_hanspub.png?20220323082705087)
Figure 5. Abundance and biomass of zooplankton in water intake and aquaculture pond
图5. 进水口、养殖池各采样点浮游动物丰度和生物量分布
3.2.3. 浮游动物多样性
本次调查进水口浮游动物Margalef物种丰富度指数(D)为3.60、Shannon物种多样性指数(H')为3.63、Pielou均匀度指数(J)为1.51。鲟鱼养殖池各采样点浮游动物Margalef物种丰富度指数(D)平均值约为4.18,介于3.60~5.72,最大值位于A池,其它各点无明显差异;Shannon-Wiener指数平均值约为3.54,介于3.32~4.16,最大值位于A池,其它各点无明显差异;Pielou均匀度指数(J)平均值约为1.38,各点无明显差异(如图6)。
![](//html.hanspub.org/file/6-2830186x17_hanspub.png?20220323082705087)
Figure 6. Distribution of diversity index of zooplankton in water intake and aquaculture pond
图6. 进水口、养殖池各采样点浮游动物多样性指数分布
3.3. 环境因子
鲟鱼集约化养殖场调查过程中环境因子变化如表6所示。其中水温恒温,为16℃,该温度最适宜鲟鱼生长;各采样点整体偏碱性,养殖池pH值均大于进水口;叶绿素Chl-a进水口较低,四个养殖池平均值为1.975;进水口总氮值为1.2 mg/L,四个养殖池平均值为3.53 mg/L,均较高;进水口氨氮值为0.108 mg/L,四个养殖池平均值为0.357 mg/L;进水口总磷值为0.025 mg/L,四个养殖池平均值为0.197 mg/L;进水口溶解氧为8.01 mg/L,四个养殖池平均值为9.30 mg/L;化学需氧量C池最高,为7.4 mg/L,B池化学需氧量最低,为5.4 mg/L。PH、总氮、氨氮、总磷值、溶解氧均为养殖池大于进水口。
![](Images/Table_Tmp.jpg)
Table 6. The main environmental parameters in water intake and aquaculture pond
表6. 采进水口、养殖池各采样点的主要环境参数
3.4. 浮游生物与环境因子的RDA分析
选取鲟鱼集约化养殖场各样点16种浮游生物优势种(Y ≥ 0.02)与的8个环境因子首先进行DCA分析,结果显示最大梯度长度为1,小于3,因此对物种数据与环境因子数据作RDA分析,研究浮游生物丰度与环境因子间的关系。结果显示:浮游植物前两个轴的特征值分别为0. 485和0.266,共解释了群落变异程度的75.09%,物种和环境因子相关系数均达到1;浮游动物前两个轴的特征值分别为0.628和0.231,共解释了群落变异程度的85.86%,物种和环境因子相关系数均达到1。上述数据表明排序能够很好地反映浮游生物丰度与环境因子之间的关系。由图7 (左)可知,鲟鱼集约化养殖场浮游植物的大部分物种分布在第二排序轴的左侧,与TN、TP、NH3-N、PH、COD、叶绿素Chl-a呈正相关,与溶解氧呈负相关;谷皮菱形藻、短鼓藻、针尖杆藻、小球藻TN、TP、NH3-N、PH、叶绿素Chl-a呈负相关,与溶解氧呈正相关;普通等片藻与PH、TD呈正相关。浮游动物优势种中,白钟虫、狭颈宿轮虫、郝氏披甲轮虫与TN、TP、NH3-N、PH、叶绿素Chl-a呈正相关,与COD、DO具有负相关性;螺形龟甲轮虫与COD、DO具有正相关性,与TN、TP、NH3-N、PH、叶绿素Chl-a负相关;尾草履虫与DO呈正相关性。
4. 讨论
4.1. 鲟鱼集约化养殖水浮游生物群落分布特征
此次调查进水口共采集到浮游植物2门16种,其中绿藻门种类占多,其次为硅藻门;采集到浮游动物3大类11种,其中种类数轮虫 > 原生动物 > 桡足类。养殖池共采集到浮游植物5门61种,其中种类数硅藻门 > 绿藻门 > 蓝藻门 > 裸藻门 > 甲藻门;共采集到浮游动物2大类39种,其中轮虫种类占多,其次为原生动物。进水口与养殖池浮游生物种类组成差别较大,A、B、C、D四个鱼池的浮游生
(左图:sp1中型脆杆藻Fragilaria intermedia,sp2近缘桥弯藻Cymbella affinis,sp3短小舟形藻Navicula exigua,sp4旋转裸囊藻Trachelomonas volvocina,sp5短肋羽纹藻Pinnularia brevicostata,sp6湖北小环藻Cyclotella hubeiana,sp7普通等片藻Diatoma vulgare,sp8谷皮菱形藻Nitzschia palea,sp9短鼓藻Cosmarium abbreviatum,sp10针尖杆藻Synedra acus,sp11小球藻Chlorella vulgaris。右图:sp1白钟虫Vorticella alba,sp2狭颈宿轮虫Habrotrocha angusticollis,sp3螺形龟甲轮虫Keratella cochlearis,sp4尾草履虫Paramecium caudatum,sp5郝氏披甲轮虫Ploesoma hudsoni)
Figure 7. RDA sequence of phytoplankton (left), zooplankton (right) and environmental factor
图7. 浮游植物(左)、浮游动物(右)与环境因子的RDA排序图
物也具有差异,这说明养殖业对水质的影响较大,且养殖池水质与养殖物种、养殖数目、鱼池大小存在重要的关联。研究表明,中国河流水体浮游植物主要以硅藻、绿藻为主 [16] [17] [18] [19],本研究也得到同样的结果。进水口处于山涧溪流水生态系统,温度较低、透明度较高、水流速度较快,适合硅藻生存 [20]。养殖鱼池硅藻大多为尖杆藻、脆杆藻等有毒性藻类,鱼类滤食掉部分可食用藻类,造成硅藻数远大于绿藻的现象。养殖鱼池浮游动物群落以小型浮游动物轮虫和原生动物为主,大型浮游动物较少,表明下级营养级捕食压力较大,而轮虫和原生动物种外竞争压力较小 [21]。进水口存在少量的无节幼体,可能是水温较低造成的 [22],而养殖鱼池中不存在无节幼体,无节幼体主要以滤食细菌为主,低温不利于细菌的生长繁殖,且养殖鱼池定期进行消毒清杀,更无法为无节幼体提供丰富的食物来源。
4.2. 鲟鱼集约化养殖水浮游生物多样性与水质状况
生物多样性指数是通过浮游生物的种类、丰度、生物量等客观反映水体环境对生物产生的综合叠加效应 [22],一般来讲,生物多样性指数越高,代表群落结构越稳定、水质状况越好。进水口浮游生物Margalef物种丰富度指数评价显示水体呈现中度污染水平,浮游生物Pielou均匀度指数和Shannon-Wiener指数评价显示水体呈现清洁状态;而养殖鱼池Margalef物种丰富度指数评价显示水体整体呈现清洁–中度污染水平,浮游生物Pielou均匀度指数和Shannon-Wiener指数评价显示水体呈现清洁状态。养殖鱼池水虽为活水,但养殖生物高密度聚集,养殖生物可排泄大量的含氮化合物和含磷化合物,且饵料投喂过程中,只有20%~30%的氮、15%~30%的磷被鱼体吸收利用,其余的则直接排入水体 [23] [24] [25],而单因素理化指标显示养殖鱼池水质TN、TP均偏高,根据《地表水环境质量标准》 [26],养殖鱼池水质整体处于V类,从中足以看出不同的水质评价指标体系显示出不同的评价结果 [27],但反观其浮游植物物种组成,却发现大多为有害种,因此各个多样性指数较好只能说明其生态系统处于较稳定的状态,进行水质生态系统健康评价还需更系统、更全面的去考虑,需同时结合水体理化因子、浮游生物指示种、浮游生物生物量及丰度的评价结果 [28]。
4.3. 环境因子对浮游生物群落结构的影响
浮游生物群落结构受多种环境因子的综合影响,比如水温、pH、溶氧、透明度和是影响浮游动物群落的直接因子,N、P等营养盐类可通过影响浮游植物群落结构分布和种类组成发生改变而间接影响浮游动物,浮游动物对浮游植物的下行控制效应,对浮游生物的群落结构有重要的影响 [29] [30] [31] [32] [33]。本研究RDA分析表明硅藻种大多数种类与TN含量关系密切,硅藻种大多数种类为淡水水域生境广布种,耐受于中污性水体,常附着在水底岩石或有机碎屑表面,喜欢在营养盐相对较高的水体环境中生活,通常高浓度的氮含量有利于维持其的持续生长 [34] [35]。而谷皮菱形藻、针尖杆藻等与TN、TP呈负相关,相关研究表明,在一些淡水水体中,氮磷等营养盐在特定环境下对浮游动植物的生长也存在抑制作用 [36] [37]。从RDA排序图可知,大多数浮游动植物生物量与Chl-a关系密切,已有研究表明Chl-a的浓度变化可反映水体中浮游植物现存量的变化 [38]。TN、TP、NH3-N与白钟虫、狭颈宿轮虫、郝氏披甲轮虫呈正相关,螺形龟甲轮虫、尾草履虫呈负相关,说明营养盐类会直接影响浮游动物群落结构分布,且浮游动物滤食浮游植物,营养盐和Chl-a通过改变浮游植物群落结构,从而间接影响浮游动物的群落结构 [31] [39]。RDA排序图不仅可全面阐述物种和环境因子间的关系,还能反映各个物种在生活习性上的差异。短小舟形藻、旋转裸囊藻、短肋羽纹藻、湖北小环藻在排序图中位置十分接近,表明它们对生态环境的适应能力接近;而中型脆杆藻与谷皮菱形藻虽都为硅藻门,但是两者在排序图位置较远,即二者对生态环境的依赖差异较大 [40]。
5. 结论
通过对鲟鱼集约化养殖水浮游生物多样性调查与分析,得出以下结论:
1) 进水口共采集到浮游植物2门16种,养殖鱼池共采集到浮游植物5门61种,其中硅藻门和绿藻门在种类上占有优势,优势种种类丰富。
2) 进水口采集到浮游动物3大类11种,养殖鱼池共采集浮游动物2大类39种,其中轮虫种类占多,其次为原生动物。
3) 通过浮游动植物的指示作用以及利用生物多样性指数进行水质评价,得出水质处于清洁–中污状态,而单因素理化指标显示养殖鱼池水质整体处于V类,因此进行水质生态系统健康评价还需同时结合水体理化因子、浮游生物指示种、浮游生物生物量及丰度等进行系统的评价。
4) 影响鲟鱼集约化养殖水浮游动植物群落结构的环境因子包括TN、TP、NH3-N、PH、COD、叶绿素Chl-a、COD。综合考虑鲟鱼集约化养殖水浮游动植物与环境因子相互影响的关系,还需建立长期的监测系统,对维持养殖鱼池生态系统的平衡稳定、评定养殖物种的健康状况等具有重大意义。
基金项目
中央支持高水平人才项目《冷水鱼资源产业化可持续利用集成技术创新》,项目编号:GSP202014。
NOTES
*并列第一作者。
#通讯作者。