Ecological health assessment of Cyprinus carpio and Leiocassis longirostris national aquatic germplasm resources reserve in Huaihe River
-
摘要: 评价水域生态健康状况可判别水生态恢复成效,是水环境管理的有效手段。为探究淮河荆涂峡鲤 (Cyprinus carpio)、长吻鮠 (Leiocassis longirostris) 国家级水产种质资源保护区的水生态健康状况,根据2020年浮游植物群落结构的调查结果,结合形态功能群分类法,构建了基于浮游植物生物完整性指数的保护区生态健康评价体系。选择受干扰较小的采样点为参照点,对候选参数进行非参数检验和相关性分析,确定最终的评价参数,采用比值法确定分级评价标准,对保护区水体进行健康评价。结果显示,浮游植物优势功能群演变趋势为III+IV+V (春季)→III+IV+V+VI (夏季)→V+VI (秋季)→IV+V+VI (冬季),除核心区的优势功能群为III+IV+V+VI外,实验区、涡河和淮河干流优势功能群均为III+IV+V+VI。保护区68.06%的样点评价结果为“一般”及以上,总体水生态健康状况良好,具体表现为冬季优于秋季、夏季次之、春季最差的时间格局,以及淮河干流优于涡河、实验区优于核心区的空间格局。Abstract: Evaluation of water ecological health status helps to judge the effectiveness of water ecological restoration, which is an effective means of water environment management. In order to explore the water ecological health status of Cyprinus carpio and Leiocassis longirostris national aquatic germplasm resources reserve in Huaihe River, according to the survey results of phytoplankton community structure in 2020 and combined with the morphological functional group classification, we constructed the ecological health evaluation system of the reserve based on the biological integrity index of phytoplankton. Besides, we carried out the Mann-Whitney non-parametric test and correlation analysis of the candidate parameters to determine the final evaluation parameters, and used the ratio method to determine the grading evaluation standard for the evaluation of the health status of the reserve water. The results show that the evolution trend of dominant functional groups of phytoplankton was III+IV+V (Spring)→III+IV+V+VI (Summer)→V+VI (Autumn)→IV+V+VI (Winter). Except that the dominant functional groups in the core area were III+ IV+V+VI, the dominant functional groups of the experimental area, the Guo river and the Huaihe River were III+IV+V+VI. 68.06% of the sampling sites in the reserve were evaluated as "average" or above, indicating that the overall water ecological health status was good. Specifically, the temporal pattern of winter was better than autumn, followed by summer, and the worst was in spring; the spatial pattern of Huaihe River mainstream was better than Guohe River, and the experimental area was better than the core area.
-
硫酸盐还原菌(sulfate-reducing bacteria,SRB)在生态系统中是属于土著微生物类群,是兼性厌氧菌,广泛存在于各种环境中,是自然界硫循环中SO42-异化还原的专门承担者,将SO42-还原同时降解有机物获取能量,产物为对周围环境有害的有毒气体H2S。目前,国内外对海洋沉积物中SRB研究多集中在海洋工程、石油勘探、管道铺设等金属防腐方面[1-2]。近年来因病害使国内外海水养殖产量出现大幅度的滑坡,其主要原因是养殖环境恶化,尤其是底质的污染严重,大部分养殖海区投喂天然饵料,饵料系数较高,养殖年限较长,网箱底部堆积大量富含有机物的沉积物,为SRB大量繁殖提供了有利条件,SRB繁殖会产生大量有毒气体H2S,造成水体污染并毒害养殖生物,使养殖水体环境状况进入恶性循环状态。目前,国内外对网箱养殖环境中SRB的相关研究报道不多[3],因此,开展网箱养殖海区沉积物中SRB生态特性和相关防治是很有必要的。本文通过研究大鹏澳网箱养殖海域的SRB,了解此环境中SRB的生态特性、SRB数量变化,以及SRB与其他环境因子之间的关系,为以后修复养殖水体环境,提高养殖生物产量,实现海水网箱养殖业的可持续发展提供理论依据。
大鹏澳是大亚湾西南部的一个小内湾,水面面积约1 400 km2。20世纪90年代以来,海水网箱养殖得到了迅速发展。目前,网箱养殖区水面面积约有20 km2,养殖网箱约4 200箱,养殖区平均水深4.5 m,污染相对比较严重[4-5]。
1. 材料与方法
1.1 样品采集
2006年11月,在大亚湾大鹏澳网箱养殖海域利用GPS定位5个站点,其中1、2、3号站位于鱼类网箱养殖区,4号站位于浮筏贝类养殖区,5号站作为对照站,采样站点如图 1所示。在每个站用柱状采泥器采集柱状沉积物样品2管,将沉积物置于内径为5 cm的PVC管内并用橡胶塞密封好,其中1管用于测定硫化物、氧化还原电位(oxidation-reduction potential,ORP)和pH,另1管用于测定SRB含量,并清晰标记样品垂直方向,其中柱状沉积物SRB计数按照1 cm层、5 cm层、10 cm层深度进行取样。
1.2 SRB培养基配方和计数方法
培养基配方为K2HPO4 0.5 g · L-1,NH4Cl 1.0 g · L-1,MgSO4 2.0 g · L-1,Na2SO4 0.5 g · L-1,CaCI2 0.1 g · L-1,酵母膏1.0 g · L-1,乳酸钠4 mL ·L-1依次加入净化陈海水中溶解。调节pH值约为7.0~7.5[6]。密封后用蒸气压力灭菌器121~125℃灭菌15~20 min后冷却至室温。称取硫酸亚铁铵1.2 g,抗坏血酸0.4 g,在无菌箱(室)内均匀地摊在离紫外线灯30 cm处灭菌30 min。在无菌操作下,把硫酸亚铁铵和维生素C溶解于事先准备好的40 mL无菌水中。按每100 mL培养基各加入1.0 mL硫酸亚铁铵溶液和1.0 mL维生素C溶液。
硫酸盐还原菌的计数(本标准适用于工业循环冷却水中硫酸盐还原菌的测定,也适用于原水、生活用水及粘泥中硫酸盐还原菌的测定)采用我国国家标准《工业循环冷却水中硫酸盐还原菌的测定》(GB/T14643.5-1993)规定的MPN法[7]。
1.3 沉积物硫化物、ORP、pH的测定
1.3.1 硫化物的测定
沉积物中硫化物的测定采用《海洋监测规范》[8]中的硫离子选择电极法。取5.00 g混匀的湿样置于50 mL烧杯中,加入20 mL抗氧化络合剂使用液,充分搅拌5 min,静置,待沉积物完全沉淀后,用倾斜法将上清液倒入50 mL量瓶中;再重复2次用15 mL抗氧化络合剂使用液按上述步骤浸取残留的沉积物,合并3次上清液,再用抗氧化络合剂使用液定容至刻度。测定时将上清液移入50 mL烧杯中,插入参比电极和活化后的硫离子选择电极,一边搅拌至读数稳定后记录结果。每次测定后均需用蒸馏水洗涤电极并用滤纸吸干,以备下一样品测定。
1.3.2 ORP的测定
沉积物ORP的测定参照《海洋监测规范》[8]中的电位计法。
1.3.3 pH的测定
沉积物pH的测定按照《海洋监测规范》[8]规定采用复合电极方法。
上述分析仪器为上海雷磁PHSJ-4A实验室pH/ORP测定仪,配有硫离子选择电极、参比电极、氧化还原电极、感温电极和pH复合电极。
2. 数据处理
利用Microsoft Excel对SRB的检出率、含量、变化幅度等按照不同站位、不同层次进行分析,利用数理统计软件SPSS对SRB与硫化物、ORP、pH做相关分析。
3. 结果与讨论
3.1 沉积物中的SRB的数量
2006年11月大亚湾大鹏澳网箱养殖海域沉积物中SRB的数量如表 1所示,沉积物中SRB数量波动范围在900~110 000 ind · g-1之间。
表 1 大鹏澳网箱养殖海域沉积物中SRB的数量Table 1. Content of SRB in cage culture area sediment in Dapeng′ao Coveind · g-1 垂直深度/cm
vertical depth1号站
station 12号站
station 23号站
station 34号站
station 4对照站
control station1 45 000 45 000 110 000 25 000 30 000 5 30 000 4 500 45 000 15 000 15 000 10 11 500 2 500 4 500 900 4 500 图 2显示了各个站位表层沉积物SRB的分布情况,SRB平面分布特征是鱼类网箱区>对照区>浮筏贝类区,网箱区的数量明显高于对照区和浮筏贝类区,这可能是由于网箱区内的网箱设置比较密集,饵料主要是冰鲜小杂鱼,残饵和粪便形成的生物沉积,为SRB的生长提供了有利的环境;浮筏贝类区SRB数量最少,可能是由于此区域的养殖生物为太平洋牡蛎,无需投饵,故底部富含有机质的沉积物较少,SRB含量也较少。
各个站位沉积物中SRB的垂直分布如图 3所示,可以明显看出SRB数量的分布特征为1 cm层>5 cm层>10 cm层,呈现由上到下逐渐减少的趋势。验证了海洋沉积物中硫酸盐还原菌数量是由上向下逐渐降低的结论[9],这也与高爱国等[10]在北极楚科奇海调查得到的SRB分布规律相似。
3.2 沉积物中的硫化物
2006年11月大亚湾大鹏澳网箱养殖海域5个站位沉积物中硫化物的含量如表 2所示,其硫化物含量波动范围为372.48~703.98 μg · g-1。
表 2 表层沉积物中SRB与其他环境因子Table 2. The amount of SRB and other factors of surface sediment站位
stationSRB数量
content of SRB硫化物含量/μg·g-1
sulfide content氧化还原电位
ORPpH 1号站 station 1 45 000 529.44 -226.5 8.12 2号站 station 2 45 000 527.76 -270.3 8.18 3号站 station 3 110 000 703.98 -384.5 8.16 4号站 station 4 25 000 372.48 -171.5 8.21 对照站 control station 30 000 445.18 -199.5 8.17 由图 4可以明显看出,鱼类网箱养殖区(1、2、3号站)沉积物中硫化物含量要高于贝类区和对照区,其主要原因是鱼类网箱区需要大量投饵,底部富含有机质,这一环境利于SRB大量繁殖,异化还原硫酸盐形成硫化物;浮筏贝类养殖区的硫化物含量最低,这是因为贝类养殖过程中不需投饵,底部沉积物有机质含量较少;另外,该海域贝类区养殖年限较短,水体底部的溶解氧和pH值均较高,不利于SRB的生长,因此,硫化物含量较低。而对照区的5号站,其硫化物含量相对较高于贝类养殖区,这可能是因为5号站位于大鹏澳湾口,在海流作用下,受网箱区沉积物往湾口外迁移的影响要大于贝类区有关。
3.3 沉积物中的ORP
各采样站位沉积物的氧化还原电位如表 2所示,其波动范围在-171.5~-384.5 mV之间。图 5显示了2006年11月大亚湾大鹏澳网箱养殖海域5个站位沉积物的ORP平面分布特征。可以看出,鱼类网箱区的ORP要低于贝类区和对照区,其主要原因是由于网箱区残饵和排泄粪便长期积累,网箱区底部富含有机质,微生物分解大量有机质,产生还原性的物质(如:S2-、NH4+、Fe2+等),从而导致氧化还原电位较低。
3.4 沉积物中的pH
图 6显示了2006年11月大亚湾大鹏澳网箱养殖海域各个站位沉积物的pH值,其特征为网箱区低于对照区和贝类区,影响沉积物pH变化的主要因素是呼吸作用[11],因为网箱区沉积物的有机质含量高于对照区和贝类区,所以网箱区沉积物中的细菌呼吸作用和有机质在细菌参与下的分解作用相对活跃,导致pH下降。
3.5 SRB与其他因子的相关关系
2006年11月大亚湾大鹏澳网箱养殖海域表层沉积物中SRB含量、硫化物含量、ORP、pH值如表 2所示。利用SPSS对SRB与其他环境因子所做的相关分析结果显示,SRB与硫化物的含量为显著性正相关(R=0.96, P < 0.05,表 3),这与李培英等[12]研究的浙江至闽北陆架积物硫酸盐还原菌与硫化物的相关关系的结果相一致。由于SRB在沉积物环境中主要电子受体是硫酸盐,当SRB大量存在时就会有更多的还原态的硫化物生成。SRB含量与Eh的关系显示极显著性负相关关系(R=-0.97,P < 0.01,表 3),这与陈皓文等[13]对北部湾东侧沉积物SRB研究结果相一致。ORP的高低主要取决于沉积物中的DO含量,当沉积物处于严重缺氧时,有机物最终分解为还原态物质,Eh为负值,而厌氧的低氧化还原状态环境适合SRB的大量繁殖,所以这一分析结果与理论相一致。而SRB与pH的关系系数为R=-2.71(P>0.05,表 3),两者不具有显著性相关关系,由此可以看出,硫化物含量和ORP是反映SRB数量的2个非常重要的指标。
表 3 SRB与其他因子的相关系数Table 3. Correlation between SRB and other chemical factors对子 parallelism 相关系数R correlation coefficient n number of sample P significance factor SRB-硫化物 SRB-sulfide 0.96 5 0.011 SRB-氧化还原电位 SRB-ORP -0.97 5 0.006 SRB-pH -2.71 5 0.611 -
表 1 候选参数
Table 1 Candidate parameters
参数种类Parametertype 候选参数Candidate parameter 群落多样性Community diversity M1. 藻类物种数 M2. 蓝藻物种数 M3. 硅藻物种数 M4. 隐藻物种数 M5. 绿藻物种数 M6. 非硅藻物种数 M7. Shannon多样性指数 M8. Pielou均匀度指数 M9. Margalef丰富度指数 M10. Simpson多样性指数 群落丰富度Community richness M11. 藻类密度 M12. 蓝藻密度 M13. 硅藻密度 M14. 隐藻密度 M15. 绿藻密度 M16. 非硅藻密度 M17. 藻类生物量 M18. 蓝藻生物量 M19. 硅藻生物量 M20. 隐藻生物量 M21. 绿藻生物量 M22. 非硅藻生物量 群落均匀度Community evenness M23. 蓝藻密度/% M24. 硅藻密度/% M25. 隐藻密度/% M26. 绿藻密度/% M27. 非硅藻密度/% M28. 蓝藻生物量/% M29. 硅藻生物量/% M30. 隐藻生物量/% M31. 绿藻生物量/% M32. 非硅藻生物量/% 形态功能群 Morphology based functional groups M33. I类功能群密度 M34. II类功能群密度 M35. III类功能群密度 M36. IV类功能群密度 M37. V类功能群密度 M38. VI类功能群密度 M39. VII类功能群密度 M40. I类功能群生物量 M41. II类功能群生物量 M42. III类功能群生物量 M43. IV类功能群生物量 M44. V类功能群生物量 M45. VI类功能群生物量 M46. VII类功能群生物量 M47. I类功能群密度/% M48. II类功能群密度/% M49. III类功能群密度/% M50. IV类功能群密度/% M51. V类功能群密度/% M52. VI类功能群密度/% M53. VII类功能群密度/% M54. I类功能群生物量/% M55. II类功能群生物量/% M56. III类功能群生物量/% M57. IV类功能群生物量/% M58. V类功能群生物量/% M59. VI类功能群生物量/% M60. VII类功能群生物量/% 表 2 参数数据分析
Table 2 Parameter data analysis
季节 Season 参数 Parameter 最大值 Maximum value 最小值 Minimum value 5%分位数 5% Quantile 95%分位数 95% Quantile 参数赋分公式 Parameter assignment formula 对干扰响应 Response to interference 春季 Spring M7 2.254 0.979 1.026 2.246 M7/2.246 下降 M45 1.732 0.000 0.002 1.213 M45/1.213 下降 M58 0.822 0.029 0.086 0.717 0.822−M58/0.822−0.086 上升 夏季 Summer M28 0.677 0.001 0.027 0.668 M28/0.668 下降 M37 2 366 163.857 47 323.277 67 435.670 1 984 028.394 2 366 163.857−M37/2 366 163.857−67 435.670 上升 M50 0.365 0.018 0.019 0.349 0.365−M50/0.365−0.018 上升 秋季 Autumn M9 1.313 0.618 0.554 1.218 0.313−M9/1.313−0.505 上升 M11 4 571 428.571 1 220 940.550 890 150.843 4 321 561.668 4 571 428.571−M11/4 571 428.571−890 150.843 上升 M15 1 703 637.977 170 363.798 170 363.798 1 590 062.112 1 703 637.977−M15/1 703 637.977−170 363.798 上升 M44 13.294 0.003 0.056 6.643 13.294−M44/13.294−0.056 上升 冬季 Winter M8 0.918 0.294 0.314 0.910 M8/0.910 下降 M9 0.876 0.390 0.394 0.849 0.876−M9/0.876−0.394 上升 M22 0.559 0.006 0.019 0.476 0.559−M22/0.559−0.019 上升 M25 0.444 0.000 0.002 0.338 M25/0.338 下降 M26 0.667 0.082 0.099 0.575 M26/0.575 下降 M43 0.254 0.001 0.001 0.241 0.254−M43/0.254−0.001 上升 表 3 保护区P-IBI评价标准
Table 3 P-IBI evaluation criteria for reserve
季节 Season 健康 Health 亚健康 Sub-health 一般 General 较差 Worse 极差 Range 春季 Spring >2.64 1.98~2.64 1.32~1.98 0.66~1.32 <0.66 夏季 Summer >2.86 2.14~2.86 1.43~2.14 0.71~1.43 <0.71 秋季 Autumn >3.83 2.87~3.83 1.91~2.87 0.96~1.91 <0.96 冬季 Winter >5.95 4.46~5.95 2.97~4.46 1.49~2.97 <1.49 表 4 四季P-IBI与环境因子逐步回归
Table 4 Stepwise regression between P-IBI and environmental factors
季节Season 逐步线性回归方程 Stepwise linear regression equation R P F 春季 Spring yP-IBI=−1.464+0.951xTP−0.268xNH3-N 0.914 0.000 38.244 夏季 Summer yP-IBI=2.639−0.304xTP +0.775xNH3-N −1.417xNO2 −-N 0.958 0.000 52.548 秋季 Autumn yP-IBI=−3.044+0.697xTP 0.697 0.001 15.094 冬季 Winter yP-IBI=17.45−0.757xNH3-N 0.757 0.000 21.458 -
[1] 马金玉, 罗千里, 王文才, 等. 华阳河湖群表层沉积物磷形态及生物有效性[J]. 长江流域资源与环境, 2021, 30(12): 2962-2971. [2] 王晨, 陶孟, 李爱民, 等. 基于环境DNA宏条形码技术的秦淮河生物多样性研究[J]. 生态学报, 2022, 42(2): 611-624. [3] 马迎群, 曹伟, 赵艳民, 等. 典型平原河网区水体富营养化特征、成因分析及控制对策研究[J]. 环境科学学报, 2022, 42(2): 174-183. [4] 张宇航, 渠晓东, 王少明, 等. 浑河流域底栖动物生物完整性指数构建与健康评价[J]. 长江流域资源与环境, 2020, 29(6): 1374-1386. [5] 王旭, 王恒嘉, 王永刚, 等. 基于大型底栖动物完整性指数(B-IBI)的密云水库上游河流 (北京段) 水生态健康评价[J]. 生态与农村环境学报, 2022, 38(2): 157-167. [6] KARR R J. Assessment of biotic integrity using fish communities[J]. Fisheries, 2011, 6(6): 21-27.
[7] 田盼盼, 桑翀, 马徐发, 等. 基于周丛藻类群落结构的新疆额尔齐斯河生态健康评价[J]. 生态学报, 2022, 42(2): 778-790. [8] 张亚军, 秦秀东, 李恩军, 等. 基于浮游生物完整性指数的不同循环饥饿投饲模式黄颡鱼养殖塘水质评价[J]. 南方农业学报, 2020, 51(7): 1754-1763. [9] 刘园园, 阿依巧丽, 张森瑞, 等. 着生藻类和浮游藻类在三峡库区河流健康评价中的适宜性比较研究[J]. 生态学报, 2020, 40(11): 3833-3843. [10] 谢孟星, 钱新, 刘彤, 等. 基于微生物完整性指数的河流健康评价——以无锡市为例[J]. 环境科学学报, 2020, 40(3): 1112-1120. [11] 董贯仓, 冷春梅, 丛旭日, 等. 南水北调东线工程运行3年后东平湖浮游植物群落特征及环境驱动因子[J]. 湖泊科学, 2022, 34(1): 61-73. [12] BAEK S H, SON M, KIM D, et al. Assessing the ecosystem health status of Korea Gwangyang and Jinhae bays based on a planktonic index of biotic integrity (P-IBI)[J]. Ocean Sci J, 2014, 49(3): 291-311.
[13] WU N, SCHMALZ B, FOHRER N. Development and testing of a phytoplankton index of biotic integrity (P-IBI) for a German lowland river[J]. Ecol Indic, 2012, 13(1): 158-167.
[14] 顾毓蓉, 薛庆举, 万翔, 等. 基于P-IBI因子分析法评价生态修复后松雅湖水生态状况[J]. 应用与环境生物学报, 2020, 26(6): 1325-1334. [15] 黎明民, 骆鑫, 付家想, 等. 基于浮游植物生物完整性的北部湾生态健康评价[J]. 中国环境监测, 2018, 34(6): 113-121. [16] KRUK C, SEGURA A M. The habitat template of phytoplankton morphology-based functional groups[J]. Hydrobiologia, 2012, 698: 191-202.
[17] SALMASO N, NASELLI-FLORES L, PADISAK J. Functional classifications and their application in phytoplankton ecology[J]. Freshw Biol, 2015, 60(4): 603-619.
[18] 郭芳, 顾继光, 赵剑, 等. 类群划分方法对南亚热带水库夏季浮游植物群落与环境响应关系的影响[J]. 环境科学, 2020, 41(11): 5050-5059. [19] 尚聪聪, 肖慧, 朱琳, 等. 海州湾大竹蛏种质资源保护区生态环境质量分析与评价[J]. 渔业科学进展, 2021, 42(1): 11-17. [20] 邱阳凌, 林育青, 刘俊杰, 等. 淮河干流及主要支流夏季浮游植物群落生物多样性评价[J]. 环境科学学报, 2018, 38(4): 1665-1672. [21] 陈豪, 左其亭, 张永勇, 等. 淮河中上游浮游植物群落结构特征及营养状况评价[J]. 中国农村水利水电, 2019(2): 1-6, 10. [22] 赵秀侠, 侯冠军, 李静, 等. 2015年淮河干流安徽段浮游植物群落结构特征[J]. 湿地科学, 2017, 15(4): 497-504. [23] 周宇建, 张永勇, 花瑞祥, 等. 淮河中上游浮游植物时空分布特征及关键环境影响因子识别[J]. 地理研究, 2016, 35(9): 1626-1636. [24] 吕海深, 丁然, 朱永华. 淮河流域闸坝对径流量的影响分析[J]. 中国农村水利水电, 2020(4): 166-171,176. [25] 章宗涉, 黄祥飞. 淡水浮游生物研究方法[M]. 北京: 科学出版社, 1995: 18-356. [26] 周凤霞. 淡水微型生物与底栖动物图谱[M]. 北京: 化学工业出版社, 2011: 37-187. [27] 胡鸿钧, 魏印心. 中国淡水藻类——系统、分类及生态[M]. 北京: 科学出版社, 2006: 23-903. [28] 郭雯, 赵俞廷, 黎宇洋, 等. 红枫湖营养水平及其历史演化趋势[J]. 生态学杂志, 2020, 39(10): 3371-3378. [29] 汪琪, 黄蔚, 陈开宁, 等. 大溪水库浮游植物群落结构特征及营养状态评价[J]. 环境科学学报, 2020, 40(4): 1286-1297. [30] 王备新. 大型底栖无脊椎动物水质生物评价研究[D]. 南京: 南京农业大学, 2003: 86-89. [31] 蔡琨, 秦春燕, 李继影, 等. 基于浮游植物生物完整性指数的湖泊生态系统评价——以2012年冬季太湖为例[J]. 生态学报, 2016, 36(5): 1431-1441. [32] 肖善势, 郝雅宾, 刘金殿, 等. 应用生物完整性指数评价钱塘江流域—浙江段水生态系统健康[J]. 水产科学, 2021, 40(5): 740-749. [33] 张辉, 彭宇琼, 邹贤妮, 等. 新丰江水库浮游植物功能分组特征及其与环境因子的关系[J]. 中国环境科学, 2022, 42(1): 380-392. [34] 李璇, 王莲, 马卫星, 等. 不同河道型微污染水源预处理工程水质净化效能及水质稳定性评价[J]. 环境科学学报, 2020, 40(8): 2771-2784. [35] 胡韧, 蓝于倩, 肖利娟, 等. 淡水浮游植物功能群的概念、划分方法和应用[J]. 湖泊科学, 2015, 27(1): 11-23. [36] 杨毓. 西安地区湖库藻类演替特征及控藻与避藻技术的研究[D]. 西安: 西安建筑科技大学, 2017: 30-31. [37] 胡月敏. 浮游植物对高原水库水环境的指示及影响藻类生长的氮磷限制研究[D]. 贵州: 贵州师范大学, 2018: 32-33. [38] 刘凌, 朱燕, 李博韬, 等. 基于MBFG分类法的长江江苏段浮游植物生物完整性评价[J]. 水资源保护, 2020, 36(4): 13-20. [39] 郑晓宇, 顾詠洁, 金妍, 等. 不同氮、磷质量浓度下四尾栅藻的生长研究[J]. 生态环境学报, 2010, 19(11): 2663-2668. [40] 赵长森, 夏军, 王纲胜, 等. 淮河流域水生态环境现状评价与分析[J]. 环境工程学报, 2008(12): 1698-1704. [41] 何保, 燕亚东, 余苇, 等. 河流交汇口河床地貌演变研究进展[J]. 水资源保护, 2018, 34(6): 17-23. [42] 米潭, 姚建. Y型河流交汇区污染物扩散模拟研究[J]. 环境科学与技术, 2020, 43(11): 9-16. [43] 田鹏, 钱昶, 林佳宁, 等. 滦河流域大型底栖动物生物完整性指数健康评价[J]. 中国环境监测, 2019, 35(4): 50-58.