基于鱼类生物完整性的青龙湾水库生态系统健康评价

王宇星1,2,叶章义1,2,罗平1,2,李聪慧1,2,陈佳1,2,刘其根1,2,胡忠军1,2*

(1.上海海洋大学 农业农村部鱼类营养与环境生态研究中心,上海 201306;2.上海海洋大学 农业农村部淡水水产种质资源重点实验室,上海 201306)

摘要:为探究安徽省宁国市青龙湾水库生态系统的健康状况,于2023年6月—2024年1月,采用环境DNA宏条形码技术(environmental DNA metabarcoding,eDNA)按季节对青龙湾水库25个断面的鱼类组成开展调查。结果表明:调查共检测出了57种鱼类,隶属6目15科47属;综合营养状态显示,该水库方面处于中营养状态;调查选取29个候选指标进行分布范围、相关关系和判别能力分析,构建了以Shannon-Wiener多样性指数、鳅科鱼类百分比、中上层鱼类百分比、亲流型鱼类百分比、耐受性鱼类百分比、产沉性卵鱼类百分比等11个核心指标生物完整性指数;通过比较比值法和赋分法与核心指标的相关性分析结果,选取比值法构建了青龙湾水库生态系统健康评价标准;结果显示,青龙湾水库处于健康(F-IBI≥6.37)的断面占24%,亚健康(6.37>F-IBI≥4.78)的断面占56%,一般(4.78>F-IBI≥3.19)的断面占20%;相关性分析显示,F-IBI指数与水深(WD)、透明度(SD)、pH、水温(WT)呈显著负相关。研究表明,基于上述11个指标构建的青龙湾水库F-IBI体系可反映青龙湾水库的生态特征,较好地体现了青龙湾水库生态系统健康状况,可作为评价该水库生态系统健康的工具。

关键词鱼类生物完整性指数;鱼类群落结构;健康评价;青龙湾水库;eDNA

鱼类作为高级消费者,在维持水生态系统结构和功能方面具有重要作用[1]。生物完整性指生物群落维持自身平衡、保持结构稳定及适应环境变化的综合能力,其核心特征是与特定生境相适应的生物群落能够实现动态稳态[2],相关研究中,Karr[3]首先将鱼类作为指示生物构建鱼类生物完整性指数(the fish index of biological integrity,F-IBI),并利用它对美国中西部河流健康状况进行了评价。近年来,大量学者运用生物完整性指数评价了河流[4]、河口[5]和水库[6]等水体的健康状态,且指示生物类群也扩展到了底栖生物[7]和浮游生物[8]等其他生物类群中。

运用F-IBI开展水生态系统健康状况评价时,全面获取鱼类物种分布信息是重要基础。长久以来,通常采用传统网具如刺网和地笼等[9]开展鱼类调查,但某些鱼类具有数量少和擅隐蔽等特点,再加上网具的选择性原因导致一些物种难以捕获,并造成鱼类物种数量被低估和生态系统健康评价结果出现偏差等问题[10]。环境DNA宏条形码技术(environmental DNA metabarcoding,eDNA)的发展为生物监测提供了新途径[11]。eDNA方法通过提取环境样本中生物释放的遗传物质,经扩增、测序和生物信息学分析,从而检测量化各种生物出现与否及其相对丰度等信息[12]。与传统调查方法相比,eDNA技术在湖泊和河流等水域的鱼类多样性监测中展现出明显优势,在稀有物种检出率和群落物种组成解析方面尤为如此[13-14]。目前,采用eDNA方法所获数据构建鱼类生物完整性指数开展健康评价较少,且仅局限于湖泊[13]和河流[15]中,在水库中的应用尚不多见。青龙湾水库地处安徽省东南部宁国市,是该市主要饮用水水源地,水库内还设有青龙湖光倒刺鲃国家级水产种质资源保护区,主要保护物种为光倒刺鲃(Spinibarbus hollandi)和光唇鱼(Acrossocheilus fasciatus)等。随着生态旅游的持续开发,频繁的人类活动可能对该水库生态健康构成潜在威胁[16]。迄今为止,有关青龙湾水库的研究主要集中在浮游植物、水化学及鱼类增养殖技术等方面[17-18],鲜有关于其鱼类组成和生态系统健康评价等方面的报道。为此,2023年6月—2024年1月,本研究团队采用eDNA方法对青龙湾水库的鱼类物种组成开展了季节调查,利用物种组成数据构建了生物完整性指数评价体系,以期掌握该水库生态系统的健康状况,并为水库生态保护和水质管理提供科学依据。

1 材料与方法

1.1 采样点布置

根据青龙湾水库形态和水文等特点,在水库设置25个断面,在每个断面设置3个采样点,采样点覆盖了水库库区及其支流,分别于2023年6月、8月、11月和2024年1月开展鱼类eDNA样品采集(图1)。

图1 青龙湾水库采样点
Fig.1 Sampling sites in Qinglongwan Reservoir

1.2 方法

1.2.1 环境DNA的采集、提取、PCR及测序 环境DNA样品采集:在每个采样点,用10%次氯酸钠消毒过的采水器分别采集表层(0.5 m)、中层和底层(根据测量的实际水深分层)水样10 L,共采集30 L水样,充分混合后留取1 L装入无菌塑料瓶中。24 h内,使用消毒过的过滤装置和0.45 μm孔径的混合纤维素滤膜进行抽滤,每次抽滤时设置一个空白滤膜(过滤1 L纯净水)作为阴性对照。采用Fast DNA Spin Kit试剂盒(MPBio,USA)根据制造商提供的方案进行DNA提取。

采用鱼类通用引物MiFish-U/R(MiFish-U-F:5′GTCGGTAAAACTCGTGCCAGC 3′和MiFish-U-R:5′CATAGTGGGGTATCTAATCCCAGTTTG 3′)对样品12S高变区域进行扩增。PCR试验采用Pro Taq、反应体系20 μL:2×Pro Taq 10 μL、Forward Primer(5 μmol/L) 0.8 μL、Reverse Primer(5 μmol/L) 0.8 μL、Template DNA 10 ng/μL、补ddH2O至20 μL 。PCR采用ABI GeneAmp® 9700型。用ddH2O为模板进行PCR阴性对照,以此评判PCR扩增过程中是否出现污染。每份样品重复3次PCR扩增过程,并将同一样品的PCR产物混匀后,用20 g/L琼脂凝胶电泳检测,表明PCR产物为单条亮带且在过滤空白和阴性对照中均无扩增结果。扩增后经20 g/L琼脂糖凝胶电泳纯化,使用TruSeqTM DNA Sample Prep Kit构建Illumina文库,使用Illumina MiSeq PE300测序仪进行配对端测序。

1.2.2 生物信息学分析 采用fastp软件对测序得到的原始序列进行质控分析,采用FLASH软件进行拼接。区分样本后采用USEARCH v7.1软件按序列相似度≥97%聚类为可操作性分类单元(Operational Taxonomic Unit,OTU),从而获得代表性序列。为了解每个OTU对应的物种,采用 RDP classifier 贝叶斯算法对OTU 代表序列进行物种分类学注释,得到相应的OTU序列丰度表。利用每个物种的reads数除以所用物种的reads数和得到该物种的相对序列丰度。序列参考数据库为Nt数据(https://www.ncbi.nlm.nih.gov/)。

1.2.3 环境变量的测量 在eDNA取样的同时,进行每个断面的环境变量的测量,采用多参数水质分析仪(YSI ProDSS,美国)现场测量水温(WT,℃)、溶解氧(DO,mg/L)、总溶解固体(TDS,mg/L)、电导率(SPC,us/cm)、pH。使用赛氏盘和声纳深度计(美国SM-5A)测量透明度(SD,m)和水深(WD,m)。采集1 L水样,24 h内按照标准流程测定叶绿素a(Chl.a,μg/L)、高锰酸盐指数(CODMn,mg/L)、总磷(TP,mg/L)、溶解性总磷(DTP,mg/L)、可溶性磷酸盐总氮(TN,mg/L)、溶解性总氮(DTP,mg/L)、硝酸氮亚硝酸氮和氨态氮(NH3-N,mg/L)。

1.2.4 综合营养状态指数 综合营养状态指数(trophic level index,TLI)是用于评价水体营养状态的综合指标体系[19]。TLI的计算公式为

(1)

(2)

式中:TLI(Σ)为综合营养状态指数;TLI(j)为第j种参数的营养状态指数;Wj为第j种参数的营养状态指数的相关权重;rij为第j种参数与叶绿素a基准参数的相关系数;m为评价参数的个数。

根据TLI值将水库的营养状态分为贫营养(TLI(Σ)≤30)、中营养(30<TLI(Σ)≤50)、轻度富营养(50<TLI(Σ)≤60)、中度富营养(60<TLI(Σ)≤70)和重度富营养(TLI(Σ)>70)。

1.2.5 参照点的设立 根据以下原则筛选出参照点:1)Shannon-Wiener多样性指数H≥2;2)断面的所有水质指标均在Ⅱ类(GB 3838—2002《地表水 环境质量标准》)以上,其余断面为受损点[20]

1.2.6 鱼类生物完整性指标建立及筛选 根据青龙湾水库鱼类组成及生态学特征,挑选了种类组成和丰度、营养结构、耐受性、繁殖共位群和鱼类数量及健康情况5大类29个候选指标,用于构建青龙湾鱼类生物完整性体系[21-22](表1)。

表1 鱼类生物完整性评价候选指标
Tab.1 Candidate metrics for F-IBI assessment

候选指标分类 types of candidate metrics候选指标 candidate metrics对干扰响应 response to interference鱼类总物种数(M1)下降Shannon-Wiener多样性指数(M2)下降均匀度指数(M3)下降丰富度指数(M4)下降鲤科鱼类占总物种数的百分比(M5)下降鳅科鱼类占总物种数的百分比(M6)下降沙塘鳢科鱼类占总物种数的百分比(M7)上升虾虎鱼科鱼类占总物种数的百分比(M8)下降种类组成和丰度species composition and abundance鳢科鱼类占总物种数的百分比(M9)上升鲇科鱼类占总物种数的百分比(M10)上升鲿科鱼类占总物种数的百分比(M11)下降鮨科鱼类占总物种数的百分比(M12)上升上层鱼类物种数的百分比(M13)下降中上层鱼类物种数的百分比(M14)下降中下层鱼类物种数的百分比(M15)下降底层鱼类物种数的百分比(M16)下降亲流型鱼类物种数的百分比(M17)下降静水型鱼类物种数的百分比(M18)上升水流广适型鱼类物种数的百分比(M19)上升肉食性鱼类数量的百分比(M20)下降营养结构 trophic structure植食性鱼类数量的百分比(M21)下降杂食性鱼类数量的百分比(M22)上升耐受性 tolerance敏感性鱼类物种数的百分比(M23)下降耐受性鱼类物种数的百分比(M24)上升产沉性卵鱼类数量的百分比(M25)下降繁殖共位群 reproductive guild产黏性卵鱼类数量的百分比(M26)上升产漂浮型卵鱼类数量的百分比(M27)下降鱼类数量及健康情况fish number and healthy status外来鱼类物种个体数百分比(M28)上升广布种个体数量百分比(出现频率大于70%)(M29)上升

候选指标的筛选通过以下3步进行:1)分布范围检测。对于候选指标中95%以上的值为 0的指标进行剔除。2)判别能力分析。绘制参照点和受损点每种指标的25%~75% 分位数箱型图,把箱体重叠情况记为IQ。当参照点与受损点箱体无重叠时,赋值IQ=3;当参照点与受损点箱体有重叠,但各自中位线均不在对方箱体内时,赋值IQ=2;当参照点与受损点箱体有重叠且仅有一方中位线在对方箱体内时,赋值IQ=1;当参照点与受损点箱体中位数值均在对方箱体内时,赋值IQ=0。选择 IQ≥2 的指标进入下一步冗余度分析。3)冗余度分析。为了检验候选指标所反映信息的独立性,计算上一步筛选出的指标之间的Spearman 相关性。当相关系数|r|>0.9 时,表示2个指标的独立性较差,保留其中包含信息量较大的指标[15]

1.2.7 赋值指标及评价标准 选用比值法和1、3、5赋分法对核心指标进行标准化。其中,比值法对受到干扰越强、分值越低的指标,按式(3)进行标准化赋分;对于受到干扰越强、分值越高的指标,按式(4)进行标准化赋分。赋分法将每个指标分为3等分,最差的等级记为1分、中间的等级记为3分,最好的等级记为5分;为了避免指标个数引起IBI总分值的不同,按式(5)换算IBI总分。

(3)

(4)

(5)

式中:为标准化后的参数值;Vi为标准化前的参数值;V5%为5%分位数值;V95%为95% 分位数值;Vmax为最大值;Pi为第i个指标的分数;n为指标个数 。

将标准化的分值加和即得到比值法的IBI 值。两种方法均以参照位点25% 分位数值为标准划分健康等级:健康、亚健康、一般、较差和极差。

1.3 数据处理

数据采用SPSS 26.0软件进行处理,以平均值±标准差(mean±S.D.)表示;显著性水平设为0.05;采用GraphPad Prism 6软件绘图。

2 结果与分析

2.1 青龙湾鱼类组成

2023—2024年4次采样分别获得了1 159 788、887 759、1 183 190、1 086 193条序列,经过注释后,共鉴定出了57种鱼类,隶属于6目15科47属(附表1,详见http://xuebao.dlou.edu.cn)。其中,鲤形目鱼类物种数最多,为40种,占全部物种数的70.2%;其次为鲈形目,鱼类物种数为8种,占全部物种数的14.0%。优势种方面:6月以翘嘴鲌(Culter alburnus,18.61%)、鲤(Cyprinus carpio,15.03%)、鲫(Carassius auratus,12.85%)的相对序列丰度最高;8月以鳙(Hypophthalmichthys nobilis,29.74%)、真吻虾虎鱼(Rhinogobius similis,20.52%)、间(Hemibarbus medius,10.34%)的相对序列丰度最高;11月以鳙(H.nobilis,40.21%)、鳜(Siniperca chuatsi,10.63%)、红鳍鲌(Chanodichthys erythropterus,7.05%)的相对序列丰度最高;1月以鳙(H.nobilis,53.11%)、鲫(C.auratus,12.56%)、(Hemiculter leucisculus,7.06%)的相对序列丰度最高。营养结构方面:肉食性鱼类17种、植食性鱼类2种、杂食性鱼类38种,分别占总鱼类物种数的29.8%、3.5%和66.7%。耐受性方面:敏感性鱼类14种、耐受性鱼类43种,分别占鱼类物种数的24.6%、75.4%。繁殖共位群方面:产沉性卵鱼类8种、产黏性卵鱼类39种、产漂浮型卵鱼类10种,分别占鱼类物种数的14.0%、68.4%和17.6%。外来入侵物种1种为蓝鳃太阳鱼(Lepomis macrochirus)。

2.2 营养状态评价

青龙湾水库各断面的综合营养状态指数(TLI)结果显示,除了S10处于轻度富营养状态外,其余断面均处于中营养状态(表5)。

2.3 鱼类生物完整性指标筛选

根据参照点的筛选标准,确定S2、S13和S21为参照点,其余为受损点。根据判别能力分析,M5、M7-M10、M12、M13、M18、M20-M22、M26-M28的IQ<2,因此从初选指标中删除。对于剩下的指标进行Spearman 相关性分析,结果表明,M2和M3的相关性为0.938,M29和M4的相关性为-0.946,M24和M23的相关性为-1.000,对于这些指标保留信息丰富的M2、M4和M24(表2)。最后得到青龙湾水库鱼类生物完整性评价的核心指标(图2)。

表2 青龙湾候选指标的Spearman相关性分析
Tab.2 Results of Spearman correlation analysis for candidate metrics in Qinglongwan Reservoir

指标metricsM2M3M4M6M11M14M15M16M17M19M23M24M25M29M2M30.938 M40.581 0.349 M60.528 0.469 0.631 M11-0.336 -0.302 -0.313 -0.334 M14-0.604 -0.442 -0.774 -0.684 0.269 M150.503 0.293 0.647 0.457 -0.364 -0.662 M16-0.081 -0.047 0.002 -0.019 0.506 -0.015 -0.336 M170.069 -0.028 0.396 0.405 -0.450 -0.438 0.313 -0.552 M19-0.373 -0.311 -0.572 -0.440 0.283 0.445 -0.254 0.183 -0.500 M230.029 0.109 -0.151 -0.059 -0.323 -0.122 0.056 -0.538 0.528 -0.413 M24-0.029 -0.109 0.151 0.059 0.323 0.122 -0.056 0.538 -0.528 0.413 -1.000 M250.438 0.456 0.221 0.239 -0.668 -0.269 0.083 -0.539 0.487 -0.397 0.510 -0.510 M29-0.563 -0.345 -0.964 -0.620 0.287 0.785 -0.652 -0.004 -0.324 0.473 0.195 -0.195 -0.230

(a)~(k)—核心指标;(l)~(m)为赋值方法;A—受损点;B—参照点。
(a)-(k)—core indicators;(l)-(m)—valuation methods;A—damage point;B——reference point.

图2 青龙湾鱼类生物完整性指数核心指标及赋值方法箱体图
Fig.2 Boxplot of core metrics and scoring methods for the fish-index of biological integrity in Qinglongwan Reservoir

2.4 青龙湾鱼类生物完整性评价

采用比值法和赋分法(表3),计算各采样点F-IBI值,得到青龙湾各断面健康状态。对两种赋值方法,按照参照点和受损点进行判别能力分析,发现两种方法均有良好的判别能力(IQ=3)。相关性分析结果显示,比值法和赋分法表现出良好的线性拟合(R2=0.78),表明两种方法具有方法学一致性。

表3 核心指标按比值和1、3、5赋分法计算标准
Tab.3 Calculation criteria for core indicators based on ratio and 1,3,5 scoring method

核心指标core metrics比值法统计分布ratio method for statistical distribution最大值maximum value5%分位数5% percentile95%分位数95% percentile比值法计算公式ratio method calculation formula赋分法统计分布scoring method for statistical distribution赋分标准scoring criteria25%分位数25% percentile75%分位数75% percentile135M22.630.542.57M2/2.571.351.96<1.351.35~1.96>1.96M432.0011.3031.70M4/31.7014.5022.00<14.5014.50~22.00>22.00M60.110.000.10M6/0.100.000.04<0.000.00~0.04>0.04M110.110.000.11M11/0.110.030.07<0.030.03~0.07>0.07M140.500.210.49M14/0.490.260.36<0.260.26~0.36>0.36M150.260.000.25M15/0.250.070.18<0.070.07~0.18>0.18M160.560.270.55M16/0.550.400.47<0.400.40~0.47>0.47M170.350.000.34M17/0.340.090.25<0.090.09~0.25>0.25M190.560.270.55(0.56-M19)/(0.56-0.27)0.350.45<0.350.35~0.45>0.45M240.850.680.84(0.85-M24)/(0.85-0.68)0.730.81<0.730.73~0.81>81.00M250.230.000.21M25/0.210.000.11<0.000.00~0.11>0.11

图3 比值法和赋分法评价结果的相关性分析
Fig.3 Correlation analysis of assessment outcomes by using the ratio method and the point-based scoring method

基于25%分数位划定比值法和赋分法的健康评价标准,分别计算两种方法对参照点和受损点评估的准确率(表4)。结果显示,比值法和赋分法对参照点评价准确率均为100%,但比值法对受损点评价的准确率为86.40%,高于赋分法的77.30%,基于以上结果选取比值法对青龙湾水库生态系统健康进行评价。

表4 青龙湾水库不同赋值方法的健康标准及准确率
Tab.4 Health standards and accuracy rates of different scoring methods for Qinglongwan Reservoir

分数位score position参照点 reference point比值法ratio method赋分法scoring method采样点sample准确率 accuracy rate/%比值法ratio method赋分法scoring method25%6.3748.125参照点reference point100100健康标准health criteria≥6.37≥48.125受损点damage point86.4077.30

按比值法计算结果,将参照点25%分数位划分为健康(F-IBI≥6.37),将25%分数位以下按照四分法划分为亚健康(6.37>F-IBI≥4.78)、一般(4.78>F-IBI≥3.19)、较差(3.19>F-IBI≥1.59)以及极差(F-IBI<1.59)。参照评价标准得到青龙湾水库各断面的健康状态(表5),其中处于亚健康的断面数最多(14个),占总断面数的56%,处在健康和一般的断面数分别有6个和5个,各占总断面数的24%和20%。青龙湾水库整体F-IBI指数为5.75,处于亚健康状态。

表5 青龙湾水库各采样点F-IBI得分、健康状况及营养状态
Tab.5 F-IBI scores,healthy status,and nutritional status of each sampling point in Qinglongwan Reservoir

断面sample比值法 ratio method赋分法 scoring methodF-IBI值 value of F-IBI健康状态 health statusF-IBI值F-IBI value 健康状态 health status营养状态nutritional statusS16.20亚健康48.13健康中营养S27.48健康59.13健康中营养S35.27亚健康37.13亚健康中营养S45.65亚健康42.63亚健康中营养S56.05亚健康48.13健康中营养S64.92亚健康39.88亚健康中营养S74.54一般34.38一般中营养S84.69一般37.13亚健康中营养S94.48一般37.13亚健康中营养S105.33亚健康45.38亚健康轻度富营养S118.07健康59.13健康中营养S127.09健康59.13健康中营养S138.38健康50.88健康中营养S144.86亚健康34.38一般中营养S156.06亚健康42.63亚健康中营养S164.23一般34.38一般中营养S174.66一般37.13亚健康中营养S185.83亚健康45.38亚健康中营养S195.85亚健康37.13亚健康中营养S205.24亚健康37.13亚健康中营养S216.37健康48.13健康中营养S225.60亚健康39.88亚健康中营养S234.95亚健康37.13亚健康中营养S246.34健康48.13健康中营养S255.57亚健康39.88亚健康中营养

2.5 鱼类生物完整性指数相关性分析

Spearman相关性分析结果显示,F-IBI指数与水深(WD)、透明度(SD)、pH、水温(WT)显著负相关,与综合营养状态指数(TLI)显著正相关(图4)。

图4 青龙湾水库F-IBI与环境因子Spearman相关性分析
Fig.4 Spearman correlation analysis between F-IBI and environmental factors in Qinglongwan Reservoir

3 讨论

3.1 青龙湾水库生态健康状况

作为水生态系统的关键生物类群,鱼类的物种多样性及其时空分布格局,可以有效表征水体健康状态和评估人为干扰强度[23]。生物完整性指数整合了多维生态指标,能系统反映水生态系统的结构和功能特征,因此可作为生态综合评价的有效工具。相较于传统调查方法,环境DNA技术展现出更高的物种检测灵敏度,如eDNA方法在英国三大湖泊中检测到的物种数显著高于刺网调查的历史记录[24],千岛湖的研究也证实其具有更优的物种多样性检测能力[25]。因此,本研究中采用eDNA技术全面解析了青龙湾水库的鱼类组成特征,为构建精准的生物完整性评价体系奠定了数据基础。

参照点、受损点的确定是构建F-IBI评价体系的关键。参照点的确立包括以下3种方法:1)采用水域的历史数据,选取人类干扰较小的断面作为参照点;2)采用定量标准化的方法;3)采用数据中最好的一项作为评价的期望值,通过对比选取参照点[26-28]。鉴于青龙湾水库历史资料匮乏,本研究中采用定量标准化方法确定参照基准。通过筛选得到了11种核心评价指标,其中Shannon-Wiener多样性指数和丰富度指数能表征青龙湾水库鱼类群落的稳定性;中上层、中下层及底层鱼类物种数占比达87.7%,可完整反映该水库鱼类的垂直分布情况;耐受性鱼类及鳅科(Cobitidae)和鲿科(Bagridae)等敏感类群的组成情况,可通过差异化环境响应机制指示生态系统受胁程度。以上多维度指标的综合应用,确保了青龙湾水库F-IBI评价体系的可靠性。

评价结果显示,青龙湾水库整体处于亚健康状态,但该水库杂食性、耐受性及产黏性卵的鱼类占比较高,均在65%以上,表明其生态系统可能仍面临较大的外来干扰,这些干扰可能主要来自生态旅游业[29]带来的人类活动增加及种植业[30]排放的农业面源污染。

3.2 影响鱼类生物完整性的因素

S7~S9和S16~S17断面的健康状态处于一般水平、其他断面都处在亚健康及以上水平,这种健康程度的空间差异性可能与人类活动和生境异质性密切相关。有研究表明,人类活动会通过改变鱼类多样性维持机制影响生态系统健康[31]。S7~S9位于保护区(S1~S13)入口区域,频繁的巡逻活动可能会引发水域扰动,导致S7~S9三个断面健康程度较低。通常而言产沉性卵的鱼类需要稳定的底质和适度的水流等生境特征,以避免沉性卵被沉积物覆盖并保障其发育所需氧气[32]。S16和S17断面水体深度较大、水域面积宽广、水体流速缓慢,可能不利于亲流型鱼类(如光倒刺鲃)和产沉性卵鱼类的分布。这种生境适配性的缺失直接影响了评价体系中的相关核心指标,进而降低了这2个断面的健康等级。

本研究表明,F-IBI指数与pH、水深和透明度显著负相关,与TLI显著正相关。pH值会直接影响鱼类的分布,通常情况下,维持鱼类正常生活的pH值范围为6.5~8.5,过高的pH值下鱼类的存活率降低,影响鱼类的分布[33- 34]。青龙湾水库中健康状态处于一般水平断面的pH平均值为8.7,超出了鱼类正常生活的阈值,而其他处于亚健康及以上状态的断面的pH平均值(8.4)较低。有研究表明,在深水水库中鱼类群落存在着明显的分层现象,鱼类会根据生活习性选择适宜的水层,因此水深会直接影响鱼类的分布[35- 36]。如鳅科鱼类属于淡水底栖中小型鱼类,大多数属种生活在流水环境中。而青龙湾水库中F-IBI指数较低的断面,水深相对较深,影响了鳅科鱼类的分布,这可能导致了水深与F-IBI指数呈负相关。营养盐可以通过上行效应影响鱼类群落,TLI和透明度分布可以正向和反向表征初级生产力[37-38],生物多样性随着营养盐浓度的增加呈现单峰变化,营养盐达到某个阈值后多样性呈现下降趋势[39]。青龙湾水库处于中营养状态,表明水库的物种多样性及其他生物类群可能处在随营养盐增加而增加的阶段,因此该水库TLI越高的断面鱼类多样性可能越大,从而有助于提高F-IBI指数。另外,青龙湾水库中透明度低、TLI大的断面大多处在较浅的水域,这也会通过上述水深的影响提高F-IBI指数。

4 结论

1)本研究利用eDNA方法收集到的鱼类数据,构建了青龙湾水库鱼类生物完整性指数(F-IBI),结果表明,青龙湾水库整体处于亚健康状态,F-IBI指数与水深等负相关,与综合营养状态指数正相关。

2)本次调查构建的包含Shannon-Wiener多样性指数,以及鳅科鱼类、鲿科鱼类、中上层鱼类、中下层鱼类、亲流型鱼类、水流广适型鱼类、耐受性鱼类、产沉性卵鱼类及广布种个体数百分比等11个指标的生物完整性指数可以作为评价青龙湾水库生态系统健康的工具。

参考文献:

[1] HE J Y,WU Z Q,HUANG L L,et al.Diversity,distribution,and biogeography of freshwater fishes in Guangxi,China[J].Animals,2022,12(13):1626.

[2] 金洪宇,李雷,张泽鹏,等.基于鱼类完整性指数评价倭肯河流域健康状况[J].水产学杂志,2023,36(6):109-117,126.JIN H Y,LI L,ZHANG Z P,et al.Assessment of ecosystem health of the Woken River on fish diversity and fish-index of biological integrity[J].Chinese Journal of Fisheries,2023,36(6):109-117,126.(in Chinese)

[3] KARR J R.Biological integrity:a long-neglected aspect of water resource management[J].Ecological Applications,1991,1(1):66-84.

[4] 李文,滕航,田辉伍,等.应用鱼类生物完整性指数评价黑水河河流健康状况[J].水生态学杂志,2024,45(2):73-83.LI W,TENG H,TIAN H W,et al.River health assessment of Heishui River using the fish-based index of biotic integrity[J].Journal of Hydroecology,2024,45(2):73-83.(in Chinese)

[5] 谢志伟,杨刚,庄平,等.基于鱼类完整性指数(F-IBI)的长江口生态健康状况评价[J].海洋渔业,2024,46(2):140-50.XIE Z W,YANG G,ZHUANG P,et al.Ecological health assessment of the Yangtze Estuary based on fish index of biological integrity(FIBI) [J].Marine Fisheries,2024,46(2):140-50.(in Chinese)

[6] 惠筠,方冬冬,贺海战,等.基于鱼类生物完整性指数的丹江口水库水生态健康评价[J].河南水产,2024(4):32-35,44.HUI J,FANG D D,HE H Z,et al.Water ecological health evaluation of Danjiangkou Reservoir based on fish biological integrity index[J].Henan Fisheries,2024(4):32-35,44.(in Chinese)

[7] 李文忠,赵子群,支翔宇,等.基于大型底栖动物生物完整性的北京山区小流域水生态健康评价[J].干旱区资源与环境,2024,38(12):118-128.LI W Z,ZHAO Z Q,ZHI X Y,et al.Assessment of water ecological health in small watersheds in the mountainous areas of Beijing based on the biological integrity of macrozoobenthos[J].Journal of Arid Land Resources and Environment,2024,38(12):118-128.(in Chinese)

[8] 郜星晨,张琪,苏巍,等.基于浮游植物生物完整性指数的金沙江下游河流生态健康评价[J].水生态学杂志,2025,46(1):11-19.GAO X C,ZHANG Q,SU W,et al.Ecosystem health assessment of lower Jinsha River based on the phytoplanktonic index of biotic integrity (P-IBI)[J].Journal of Hydroecology,2025,46(1):11-19.(in Chinese)

[9] DUKAN N,CORNELIS I,MAES S,et al.Vertical and horizontal environmental DNA (eDNA) patterns of fish in a shallow and well-mixed North Sea area[J].Scientific Reports,2024,14(1):16748.

[10] WANG S P,YAN Z G,HNFLING B,et al.Methodology of fish eDNA and its applications in ecology and environment[J].Science of the Total Environment,2021,755:142622.

[11] STOECKLE M Y,ADOLF J,CHARLOP-POWERS Z,et al.Trawl and eDNA assessment of marine fish diversity,seasonality,and relative abundance in coastal New Jersey,USA[J].ICES Journal of Marine Science,2021,78(1):293-304.

[12] ZHU M Y,KUROKI M,KOBAYASHI T,et al.Comparison of fish fauna evaluated using aqueous eDNA,sedimentary eDNA,and catch surveys in Tokyo Bay,Central Japan[J].Journal of Marine Systems,2023,240:103886.

[13] 郭昌陇,张鹏,郑伟,等.基于eDNA鱼类群落监测的武汉市及周边城市湖泊水生态健康评价[J/OL].湖泊科学,1-16[2025-05-26].http://kns.cnki.net/kcms/detail/32.1331.P.20250121.1749.006.html.GUO C L,ZHANG P,ZHENG W,et al.Aquatic ecosystem health assessment of Urban Lakes in Wuhan and surrounding cities based on eDNA fish community monitoring[J/OL].Journal of Lake Sciences,1-16[2025-05-26].http://kns.cnki.net/kcms/detail/32.1331.P.20250121.1749.006.html.(in Chinese)

[14] 言柯程,李建超,田永军,等.基于环境DNA metabarcoding和底拖网调查的南黄海西部鱼类多样性比较[J].中国海洋大学学报(自然科学版),2023,53(5):71-81.YAN K C,LI J C,TIAN Y J,et al.Comparison of fish biodiversity in western south Yellow Sea based on environmental DNA metabarcoding and trawl survey[J].Periodical of Ocean University of China(Natural Science Edition),2023,53(5):71-81.(in Chinese)

[15] 郭宁宁.基于eDNA技术研究赤水河流域鱼类分布及影响因素[D].北京:中国环境科学研究院,2023.GUO N N.Study on fish distribution and influencing factors in Chishui River basin based on eDNA technology[D].Beijing:Chinese Academy of Environmental Sciences,2023.(in Chinese)

[16] 李丽玲.可持续旅游发展与生态文明规划理论与实践[J].黑龙江环境通报,2025,38(1):38-40.LI L L.Theory and practice of sustainable tourism development and ecological civilization planning[J].Heilongjiang Environmental Journal,2025,38(1):38-40.(in Chinese)

[17] 贾俊威.港口湾水库斑点叉尾鮰网箱养殖技术[J].现代农业科技,2012(4):324-325.JIA J W.Cage culture technology of channel catfish in Gangkouwan reservoir[J].Modern Agricultural Science and Technology,2012(4):324-325.(in Chinese)

[18] 赵秀侠,卢文轩,李静,等.青龙湖冬季浮游植物群落结构及其与环境因子的关系[J].水生态学杂志,2020,41(2):83-90.ZHAO X X,LU W X,LI J,et al.Winter phytoplankton community structure and its relationships with environmental factors in Qinglong Lake,Anhui Province[J].Journal of Hydroecology,2020,41(2):83-90.(in Chinese)

[19] 齐荣斌.河北衡水湖营养状态变化特征与预测[J].水利技术监督,2025,33(2):166-168.QI R B.Variation characteristics and prediction of nutritional status in Hengshui Lake,Hebei Province[J].Technical Supervision in Water Resources,2025,33(2):166-168.(in Chinese)

[20] 潘丽波,黄雪妍,刘晶晶,等.基于鱼类生物完整性的赤水河健康评估研究[J].环境科学研究,2023,36(8):1532-1542.PAN L B,HUANG X Y,LIU J J,et al.Chishui river health assessment based on fish biotic integrity index[J].Research of Environmental Sciences,2023,36(8):1532-1542.(in Chinese)

[21] 李慧清,罗遵兰,孙光,等.应用鱼类生物完整性指数评价北运河河流健康状况[J].水生态学杂志,2023,44(1):34-40.LI H Q,LUO Z L,SUN G,et al.Health assessment of north canal river based on the fish-index of biotic integrity[J].Journal of Hydroecology,2023,44(1):34-40.(in Chinese)

[22] 郜星晨,张琪,苏巍,等.基于鱼类生物完整性的金沙江黑水河健康评估及修复模式探讨[J].长江流域资源与环境,2022,31(1):104-112.GAO X C,ZHANG Q,SU W,et al.Environmental quality assessment based on fish biotic integrity index and eco-restoration mode discussion for the Heishui River,Jinsha River[J].Resources and Environment in the Yangtze Basin,2022,31(1):104-112.(in Chinese)

[23] CHEN D Q,XIONG F,WANG K,et al.Status of research on Yangtze fish biology and fisheries[J].Environmental Biology of Fishes,2009,85(4):337-357.

[24] HNFLING B,HANDLEY L L,READ D S,et al.Environmental DNA metabarcoding of lake fish communities reflects long-term data from established survey methods[J].Molecular Ecology,2016,25(13):3101-3119.

[25] HAO L,GU K D,ZHOU Y,et al.Comparing diversity and structure of freshwater fish assemblages using environmental DNA and gillnetting methods:a case study of a large deep reservoir in East China[J].Ecological Indicators,2024,166:112538.

[26] 凌文翠,武亚林,方瑶瑶,等.基于鱼类生物完整性的北京市山区河流生态健康评价[J].生态毒理学报,2024,19(4):131-142.LING W C,WU Y L,FANG Y Y,et al.Ecosystem health evaluation based on fish-index of ecological integrity in Beijing mountainous area rivers[J].Asian Journal of Ecotoxicology,2024,19(4):131-142.(in Chinese)

[27] 牛明香,王俊,左涛,等.基于鱼类生物完整性指数的黄河口近岸渔业水域健康评价[J].水生态学杂志,2023,44(6):45-52.NIU M X,WANG J,ZUO T,et al.Health assessment of coastal fishery waters in the Yellow River Estuary based on the fish index of biotic integrity[J].Journal of Hydroecology,2023,44(6):45-52.(in Chinese)

[28] 胡兴坤,阙延福,李伟涛,等.基于种类数基准值优化设定的汉江中下游水生生物完整性指数评价[J/OL].水生态学杂志,1-14[2025-05-26].https://doi.org/10.15928/j.1674-3075.202407050258.HU X K,QUE Y F,LI W T,et al.Evaluation of index of aquatic biological integrity in the middle and lower reaches of Hanjiang River based on optimal setting of species number reference value[J/OL].Journal of Hydroecology,1-14[2025-05-26].https://doi.org/10.15928/j.1674-3075.202407050258.(in Chinese)

[29] 朱金宜.安徽宁国:“落羽杉红了,乡村游火了” [N].人民日报(海外版),2023-10-26(11).ZHU J Y.Ningguo,Anhui:“The feather fir is red,and the countryside is on fire”[N].People’s Daily (Overseas Edition),2023-10-26(11).(in Chinese)

[30] 邓佳良,黄黄,王耀,等.安徽省宁国市山核桃种植区典型滑坡机制与稳定性分析:以石岭村滑坡为例[J].中国地质调查,2024,11(1):93-104.DENG J L,HUANG H,WANG Y,et al.Mechanism and stability analysis of typical landslide in pecan planting area of Ningguo City in Anhui Province:a case study of Shilingcun landslide[J].Geological Survey of China,2024,11(1):93-104.(in Chinese)

[31] BURGAD A A,ADAMS G L,ADAMS R.Temporal and spatial dynamics of fish community structure during watershed alteration in two Ouachita River systems[J].Ecology of Freshwater Fish,2019,28(3):459-472.

[32] 朱其广,张琪,杨志,等.三峡库区支流磨刀溪产粘沉性卵鱼类早期资源现状[J].水生态学杂志,2023,44(1):101-107.ZHU Q G,ZHANG Q,YANG Z,et al.Status of early stage resources of fish with sticky and demersal eggs in modaoxi river,a tributary of the Three Gorges Reservoir[J].Journal of Hydroecology,2023,44(1):101-107.(in Chinese)

[33] 肖莹洁,张洁,徐鹏,等.赣江中游鱼类群落结构及多样性[J].江西科学,2024,42(3):490-497.XIAO Y J,ZHANG J,XU P,et al.Community structure and diversity of fish in the middle reaches of Ganjiang River[J].Jiangxi Science,2024,42(3):490-497.(in Chinese)

[34] 刘成凤,严云志,韩胜盼,等.基于eDNA宏条形码技术的巢湖夏季鱼类群落结构及多样性[J/OL].水生态学杂志,1-11[2025-05-26].https://doi.org/10.15928/j.1674-3075.202305040116.LIU C F,YAN Y Z,HAN S P,et al.Fish community structure and diversity in summer of Chaohu Lake based on eDNA metabacroding technology[J/OL].Journal of Hydroecology,1-11[2025-05-26].https://doi.org/10.15928/j.1674-3075.202305040116.(in Chinese)

[35] MARTIN B T,CZESNY S J,WAHL D H.Vertical distribution of larval fish in pelagic waters of southwest Lake Michigan:implications for growth,survival,and dispersal[J].Journal of Great Lakes Research,2011,37(2):279-288.

[36] XIE H Y,CAI Y C,LI C H,et al.Unveiling vertical migration patterns of the mesopelagic fish Ceratoscopelus warmingii in cold seep ecosystems:a trophic niche perspective[J].Ecological Indicators,2024,166:112533.

[37] 代培,阎明军,周游,等.太湖五里湖沿岸带浮游植物群落生态特征(2014—2015年)[J].长江流域资源与环境,2018,27(10):2348-2357.DAI P,YAN M J,ZHOU Y,et al.Ecological characteristics of phytoplankton community structure in the littoral zone of Lake Wuli,lake Taihu in 2014-2015[J].Resources and Environment in the Yangtze Basin,2018,27(10):2348-2357.(in Chinese)

[38] 袁晨,裴瑞华,郑梓瑶,等.广东东莞松木山水库鱼类群落结构及其与环境因子的关系[J].上海海洋大学学报,2022,31(4):940-949.YUAN C,PEI R H,ZHENG Z Y,et al.Structural characteristics of fish assemblage and its relationship with environmental factors in Songmushan Reservoir,Dongguan City of Guangdong Province[J].Journal of Shanghai Ocean University,2022,31(4):940-949.(in Chinese)

[39] CHEN L J,LIU Q,PENG Z R,et al.Rotifer community structure and assessment of water quality in Yangcheng Lake[J].Chinese Journal of Oceanology and Limnology,2012,30(1):47-58.

Assessment of ecosystem health in Qinglongwan Reservoir based on the Fish Index of Biological Integrity

WANG Yuxing 1,2,YE Zhangyi 1,2,LUO Ping 1,2,LI Conghui 1,2,CHEN Jia 1,2,LIU Qigen 1,2,HU Zhongjun 1,2*

(1. Ministry of Agriculture and Rural Affairs Fish Nutrition and Environmental Ecology Research Center,Shanghai Ocean University,Shanghai 201306,China;2. Ministry of Agriculture and Rural Affairs Freshwater Aquatic Germplasm Resources Key Laboratory,Shanghai Ocean University,Shanghai 201306,China)

AbstractTo study the health status of the ecosystem in Qinglongwan Reservoir, Ningguo Cith,Anhui Province, we investigated the fish composition in 25 sites in the Qinglongwan Reservoir seasonally from June 2023 to January 2024. Using eDNA methods, we identified 57 fish species,which were classified into six orders, 15 families, and 47 genera. The overall nutritional status indicated that the Qinglongwan reservoir is in a mesotrophic state. We selected 29 candidate indicators for analysis of distribution range, correlation, and discriminative ability and ultimately constructed a fish index of biological integrity (F-IBI) based on 11 core indicators, including the Shannon-Wiener diversity index, percentage Cobitidae fish species, percentage of mid-upper layer fish, percentage of rheophilic fish, percentage of tolerant fish, and percentage of fish that produce demersal eggs. Considering comparative analysis of correlation between the ratio method,scoring method and the core indicators, we chose to use the ratio method to establish the ecological health assessment standard for the Qinglongwan Reservoir. The results showed that 24% of the sites in the Qinglongwan Reservoir were healthy (F-IBI≥6.37), 56% were sub-healthy (6.37>F-IBI≥4.78), and 20% were general (4.78>F-IBI≥3.19). Correlation analysis indicated that the F-IBI index was significantly negatively correlated with water depth, transparency, pH, and water temperature. The results suggest that the F-IBI system constructed based on the 11 indicators reflects the ecological characteristics of the Qinglongwan Reservoir and adequately represents the health status of the reservoir’s ecosystem. Therefore, it can serve as a tool to evaluate the healthy status of the reservoir’s ecosystem.

Key wordsfish index of biological integrity; fish community structure; health assessment; Qinglongwan Reservoir; environmental DNA

中图分类号X 826

文献标志码:A

收稿日期2025-02-19

基金项目国家重点研发计划(2023YFD2400901)

作者简介王宇星 (1998—),男,硕士研究生。E-mail:WangYX981223@163.com

通信作者胡忠军 (1975—),男,博士,副教授。E-mail:zjhu@shou.edu.cn

DOI:10.16535/j.cnki.dlhyxb.2025-031

文章编号:2095-1388(2025)05-0843-11