环境科学  2023, Vol. 44 Issue (11): 6116-6124   PDF    
基于宏条形码技术的白洋淀水华藻类识别及其驱动因子分析
陈婷1, 杜珣2,3, 陈义永2, 郭逍宇1, 熊薇2     
1. 首都师范大学资源环境与旅游学院, 北京 100048;
2. 中国科学院生态环境研究中心, 北京 100085;
3. 中国科学院大学, 北京 100049
摘要: 浮游藻类是引起水华暴发的主要原因.为筛选潜在水华藻类,评估白洋淀水华风险区域,于2020年8月对白洋淀373点位展开浮游藻类调查.利用宏条形码技术分析,解析水华藻类群落组成,同时采用显微镜计数法统计藻密度.根据总藻密度对白洋淀不同区域的水华程度进行评估,同时进一步针对水华藻类群落,耦合淀区水质条件,探究白洋淀不同区域水华藻类群落空间差异驱动因子,以甄别影响水华藻类群落结构关键环境因子.结果表明,95%以上采样区域无水华风险(藻类密度 < 2×106个·L-1),仅5个样点存在轻微水华风险.但水华藻类群落分析共检测到了90种水华藻类,其中优势水华藻种有20种,隶属于以绿藻门、蓝藻门和裸藻门为主.水华藻类群落结构在不同区域上具有显著空间异质性(P < 0.05).关键驱动因子解析结果表明,总磷(TP)、总氮(TN)和氨氮(NH4+-N)是造成水华藻类群落结构差异的关键因子.其中,门水平上,蓝藻门水华藻类与以上关键因子显著正相关;种水平上,硅藻门和绿藻门水华藻类与关键因子响应更显著.因此,水华藻类群落的研究结果表明,水华藻类的种类和丰度与水体营养盐浓度密切相关,建议控制白洋淀水体氮磷等营养盐的输入,防止水华暴发.
关键词: 白洋淀      宏条形码      水华藻类      群落结构      水污染      有害藻类     
Metabarcoding Profiling of Phytoplankton Communities Associated with Algal Blooms and Determining Related Drivers in Baiyangdian Lake
CHEN Ting1 , DU Xun2,3 , CHEN Yi-yong2 , GUO Xiao-yu1 , XIONG Wei2     
1. College of Resouurces Environment and Tourism, Capital Normal University, Beijing 100048, China;
2. Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China;
3. University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: Phytoplankton are the main cause of algal blooms. To identify bloom algae and assess the risks of the algal blooms in Baiyangdian Lake, a survey on 373 sites was conducted in August 2020. The phytoplankton were studied via both morphological-based density counting and metabarcoding profiling. Then, the bloom degree was classed according to algae density, and the relationship between the community of bloom algae and environmental variables were modeled to determine key factors constraining spatial variation in bloom algae communities. The results showed that more than 95% of the sampling sites were free from the risk of algal blooms(phytoplankton density < 2×106 cells·L-1), and only five sites had a slight risk of algal blooms. A total of 90 species with potential of algal blooming were detected, including 20 dominant species, which were mainly affiliated with Chlorophyta, Cyanophyta, and Euglenophyta. Communities of bloom algae significantly varied among different regions(P < 0.05). Total phosphorus(TP), total nitrogen(TN), and ammonia nitrogen(NH4+-N) were the key factors significantly affecting the spatial variation in algal bloom communities. At the phylum level, these key factors were significantly positively correlated with Chlorophyta, whereas at the species level, species in Bacillariophyta and Chlorophyta responded significantly to these key factors. Thus, our findings suggested that nutrient levels were significantly related to bloom algae communities, and we proposed that controlling the input of nutrients such as nitrogen and phosphorus and regulating the hydrological process of the lake would be effective management techniques to prevent algal blooms in Baiyangdian Lake.
Key words: Baiyangdian Lake      metabarcoding      algal blooms      community structure      water pollution      harmful algae     

水华是指由于水体富营养化等导致的藻类等浮游植物快速繁殖, 聚集在水体表面, 从而导致水体颜色呈绿色、蓝色或暗褐色的现象.水华是浅水湖泊重要的危害之一, 可造成水体质量急剧恶化, 严重危害淡水生态环境[1, 2].各国政府或组织已制定相关标准, 以加强对水华的监管和防控.例如:加拿大政府对娱乐用水限值为蓝藻密度不大于1×108个·L-1[3].世界卫生组织(World Health Organization, WHO)规定娱乐用水蓝藻密度<2×107个·L-1为低风险、2×107~10×107个·L-1为中等风险、10×107~100×107个·L-1为高风险水平[4].我国生态环境部2020年通过的《水华遥感与地面监测评价技术规范(试行)》(HJ 1098-2020)中根据总藻密度(个·L-1)将水华程度分为5类级别评价.水华暴发的实质是藻类等初级生产力的急剧增大的现象, 因此控制水体中初级生产力的营养物质是防控水华暴发的关键[5].另外, 大多数学者研究倾向于从众多复杂因子中寻找关键因子, 以更好地开展水华的预测和治理[6~8].

水华藻类是引起水华暴发的主要浮游藻类群落, 常见水华藻属有微囊藻属、鱼腥藻属、裸藻属和小环藻属等[9].对水华藻类的研究主要集中在蓝藻, 硅藻等具体优势类群上[10, 11], 针对全部水华藻类群落的研究鲜见.藻类鉴定是研究水华藻类的基础步骤, 传统方法通常借助显微镜对藻类进行鉴定和计数统计[12], 但是以多样性组成为目的的研究往往遇到形态难以区分, 费时费力的困境, 尤其针对大批量样品分析的情况.基于高通量测序的宏条形码技术(metabarcoding)通过多年发展, 具有高效、可靠和易于标准化等优点[13].目前, 基于宏条形码技术的浮游藻类的研究包括了提高浮游藻类鉴定的准确性和特异性的方法学研究[14], 以及浮游藻类群落多样性[15]、浮游藻类群落地理分布格局和影响因素[16]的理论性研究, 和评估淡水系统生态状况[17]应用性研究.

白洋淀位于河北省雄安新区, 属海河流域大清河南支水系湖泊, 是华北平原最大的淡水湖泊, 水域面积约366 km2.白洋淀是雄安新区的重要生态水体, 具有缓洪滞沥、调节局部气候、改善生态环境和保护生物多样性等多种生态功能[18].然而, 长期以来在人类活动和气候变化的综合作用下, 白洋淀湿地出现水源不足、水体污染、生物多样性减少和富营养化加剧等危机[19].据报道, 白洋淀近5年整体富营养化水平均为中度富营养, 未发生大面积大规模水华, 但小部分区域仍存在较高水华暴发风险[20].本研究利用传统显微镜计数法统计藻密度评估水体水华程度, 确定白洋淀潜在水华风险区域, 同时基于宏条形码技术筛选水华藻类, 结合淀区水质条件, 甄别影响水华藻类群落结构的关键环境驱动因子, 以期为白洋淀水华防控提供理论支撑.

1 材料与方法 1.1 样点布设

本研究于2020年夏季(8月10~18日)对白洋淀进行水质和浮游藻类样品采集.依据5大功能区特征共布设373个样点.具体如图 1所示, 包括北部入流区114个(N1~N114)、汇流混合区140个(H1~H140)、南部入流区69个(S1~S69)、西部入流区3个(W1~W3)和西南部入流区47个(WS1~WS47); 同时, 布设样点分布于3种不同的生境类型, 包括淀面(DM, 143个)、沟壕(GH, 64个)和水道(SD, 166个).

N.北部入流区, H.汇流混合区, S.南部入流区, W.西部入流区, WS.西南部入流区 图 1 全淀区浮游植物基线调查样点分布 Fig. 1 Sampling sites of phytoplankton in Baiyangdian Lake

1.2 样品采集与处理

浮游藻类:从底层到表层分别取3次水, 每次10 L, 混合后通过20 μm孔径浮游生物网过滤, 收集滤网内所有存留物, 加无水乙醇固定, 并定容到100 mL, 样品置于-20℃冰箱保存用于宏条形码实验.用于显微镜计数统计的定量样品则重复同样操作, 加入适量的鲁哥氏液固定, 并定容到100 mL.为避免样点之间交叉污染, 采样前用样点原位水彻底冲洗采水器及生物网.

水质样品:现场使用多参数水质检测仪(MYRON Company, USA)测定表层水(水下0.5 m处)的水温(T)、溶解氧(DO)、pH、总溶解性固体(TDS)和氧化还原电位(ORP)等.另外, 采集1 L表层水用于实验室测定总氮(TN)、总磷(TP)和氨氮(NH4+-N)等水质因子.实验室测定前, 所有样品均在4℃保存, 并于48 h内测定完成.水质测定方法参照《水和废水监测分析方法(第四版)(增补版)》.

1.2.1 藻密度分析方法和水华评价标准

针对鲁哥氏液保存样品, 在实验室内静置24~48 h, 用虹吸管缓慢吸走上清液, 保留50 mL浓缩样品作为镜检样品.藻密度参照《内陆水域浮游植物监测技术规程》(SL 733 2016)计算.藻密度结果参考《水华遥感与地面监测评价技术规范(试行)》(HJ 1098-2020)中水华评价分级标准(表 1)对白洋淀各样点水华风险进行评估.

表 1 水华评价分级标准 Table 1 Classification criteria for algal bloom evaluation

1.2.2 宏条形码分析实验

针对每个乙醇保存样品, 首先, 采用试剂盒DNeasy Blood and Tissue Kit(Qiagen, 69506)提取浮游藻类总DNA.然后, 采用藻类通用18S rDNA引物1389F:TTGTACACACCGCCC和1510R: CCTTCY GCAGGTTCACCTAC以及16S rDNA引物515F: GTGCCAGCMGCCGCGGTAA和806R: GGACTAC HVGGGTWTCTAAT作为分子标记分别对真核和原核藻类进行PCR扩增, 将每个样品的扩增引物加上特异性8碱基标签进行混合建库.扩增体系和程序设置参考相似研究[21].每个样品进行3次重复PCR扩增以避免扩增偏差.利用2%的琼脂糖凝胶电泳检测扩增结果, 并采用SanPrep Column PCR产物纯化试剂盒对PCR产物进行纯化.最后, 利用紫外分光光度计NanoDrop 2000进行PCR产物纯化质量和浓度测定, 并按照等摩尔量原则将PCR产物混合建库, 采用Illumina NovaSeq平台PE150和PE250双端测序法对18S和16S文库进行测序.

1.3 生物信息学分析

生物信息学分析主要通过本地服务器平台Galaxy(http://mem.rcees.ac.cn:8080/)完成.首先, 根据特异性的标签将获得的原始序列分配到每个样品, 然后进行序列拼接.以每5个碱基为一个操作窗口(windows), 删除:①序列标签(tags)或引物有错配的序列, ②平均Q值(Phred scores)<20的序列, ③包含任何不确定碱基的序列, 18S和16S保留序列长度区间分别为115~125 bp和248~258 bp.采用UNOISE算法对滤错后的数据进行聚类, 得到ZOTUs和每个ZOTU对应的序列信息, 即ZOTU表格.最后利用Seed软件, 真核藻类用PR2, 原核用SILVA数据库为参照, 对得到每个OTU的物种信息进行分类注释.其中, 比对结果筛选标准为:E值<10-80, 最小覆盖度(coverage)>85%, 序列相似度(identity)>90%[22, 23].结合水华藻类历史资料筛选出水华藻类物种, 从而揭示各研究点位的水华藻类群落.

1.4 生态学与统计学分析 1.4.1 环境差异性分析

为了分析不同区域及生境类型采样点的环境异质性, 首先采用主成分分析(principal component analysis, PCA)对不同样点之间水质环境数据进行聚类.然后, 采用相似性分析(analysis of similarity, ANOSIM)对不同区域类型及生境类型间的差异程度进行显著性检验.在统计分析之前, 除pH值外, 其余环境变量采用log10(x+1)转换.所有分析均基于所有环境变量的欧式距离计算.PCA分析在R语言环境由rda函数实现[24], ANOSIM分析在Primer 5.0软件中实现[25].

1.4.2 浮游植物群落多样性分析

采用Primer 5.0软件计算水华藻类群落Shannon-Wiener多样性指数和Pielou均匀度指数[26, 27], 并用Origin 2021绘制α多样性指数箱线图.利用Excel计算各水华藻类优势度以确定淀区内优势水华藻种.

1.4.3 环境与生物群落相互关系分析

采用冗余分析(redundancy analysis, RDA), Mantel相关分析和Pearson相关性分析多角度探究水华藻类群落与环境因子之间的关系.在RDA分析过程中, 采用方差膨胀因子(variance inflation factor, VIF)检验去除共线性强的环境因子, 前向选择(forward selection)和蒙特卡罗检验(1 000次置换检验)解析影响水华藻类群落的关键因子.前向选择、VIF检验和蒙特卡罗检验和RDA分析均在软件Canoco 5.0实现.Mantel相关性分析和Pearson相关性分析在R语言环境由ggcor包实现.

2 结果与分析 2.1 白洋淀水体环境的空间异质性

结合环境因子PCA排序分析和ANOSIM相似性分析结果表明, 不同区域之间环境条件显著不同, 且不同生境类型之间环境条件显著不同.PCA分析中, 主成分PC1方差贡献率为28%, 其中TP、TN和NH4+-N因子载荷比重较大, 且与PC1正相关; 同理主成分PC2方差贡献率为21.1%, 其中DO和pH因子载荷比重较大, 且与PC2正相关.ANOSIM相似性分析结果显示(表 2), 不同区域之间的水体环境具有显著差异, 两两区域间差异程度不同, 其中H与N区水体环境差异最大; 同理, 不同生境之间的环境差异显著, DM与GH和SD水体环境均差异极显著.

表 2 白洋淀环境因子相似性的显著性检验结果1) Table 2 Significance test results for ANOSIM in Baiyangdian Lake

2.2 水华风险区域评估

依据藻密度的水华风险区域评估结果表明, 95%以上样点无水华风险(藻密度<2×106个·L-1), 仅在淀面(点位编号:DM135, 位于N区; DM166, 位于H区)与水道(点位编号:SD224, 位于N区; SD17和SD245, 位于WS区)共5个点位检测到藻类密度为2×106~1×107个·L-1, 达到轻微水华程度.

2.3 白洋淀水华藻类组成及分布特征

根据宏条形码分析结果, 结合相关文献及历史资料, 在白洋淀整个研究区域共检出90种水华藻类, 隶属于7个门, 包括金藻门、隐藻门、裸藻门、甲藻门、硅藻门、绿藻门和蓝藻门.其中, 从藻类丰度看, 绿藻门水华藻类相对丰度最高(38%), 裸藻门次之(27%), 隐藻门最少(0.04%).从种类上看, 水华藻群落以绿藻门和蓝藻门水华种类为主, 其中绿藻门水华藻类30种, 占33.33%, 包括1纲、1目、4科、4属和30种; 蓝藻门24种, 占26.67%, 包括2纲、2目、3科、9属和24种.因此, 无论从丰度还是物种数来看, 白洋淀最主要的水华藻类为绿藻门.同时, 裸藻门和蓝藻门也是白洋淀重要的水华藻类.

对群落α多样性的分析发现不同区域水华藻类α多样性存在差异(图 2), 基于不同区域, N区水华藻类群落Shannon-Wiener指数显著低于H和S区(P<0.01), S区域Evenness指数显著高于H区(P<0.001).基于不同生境, Shannon-Wiener和Evenness指数均值显示, GH生境水华藻类的多样性略高于DM和SD.但是, 不同生境类型水华藻类α多样性差异不显著(P>0.05).进一步ANOSIM相似性分析表明不同区域之间的水华藻类群落结构具有显著差异(P=0.001, R=0.185).但是不同生境类型之间的水华藻类群落结构差异性不显著(P=0.43, R=0.002).

Kruskal-Wallis F检验, *表示P<0.05, **表示P<0.01, ***表示P<0.001 图 2 不同区域和不同生境α多样性指数的箱线图 Fig. 2 Boxplots of α diversity index in different regions and habitats

优势度计算结果表明90种水华藻类中共有20种优势水华藻种, 其中包括绿藻门与裸藻门藻类各6种, 蓝藻门藻类5种, 硅藻门藻类2种, 甲藻门1种(表 3).H区优势水华藻种最多, 有14种, 占所有优势种70%.此外, 5种水华藻种在所有分区内均为优势物种, 分别是活泼裸藻(Euglena agilis)、Cyclotella scaldensis、冠盘藻(Stephanodiscus hantzschii)、Chlorella sp.和Chlamydomonas sp.

表 3 白洋淀不同区域水华浮游藻类优势种1) Table 3 Dominant species in different regions of Baiyangdian Lake

2.4 水华藻类群落结构与环境因子的关系

水华藻类群落与白洋淀环境因子间进行RDA分析.先通过方差膨胀因子(VIF)分析, 剔除共线性显著的TDS和COND(VIF>10).通过前向选择和

蒙特卡罗检验, 筛选出TP(R2=0.141, P=0.001)、TN(R2=0.118, P=0.001)和NH4+-N(R2=0.163, P=0.001)这3个环境因子是造成藻华藻类群落结构差异的主导因子.根据筛选出的关键环境因子构建RDA模型(图 3), 其中所示的环境因子共解释了水华藻类群落结构差异的42.3%, RDA排序图前两轴累计解释率达到98.26%, RDA1为93.41%, RDA2为4.85%.图 3显示了丰度前10的水华藻种与水体关键环境因子的响应规律, 可以发现不同的水华藻类与3个关键因子间响应规律呈现出差异, 以NH4+-N为例, 常温衣藻(Chlamydomonas monadina)和冠盘藻(Stephanodiscus hantzschii)等与NH4+-N呈正相关, 而Gymnodinium sp.和Euglena sp.则呈负相关(图 3).

图 3 白洋淀不同区域水华藻类群落与水体环境因子的RDA分析 Fig. 3 RDA biplot of algal bloom communities and water environmental factors in different areas and habitats of Baiyangdian Lake

为了进一步研究水华藻类与环境因子间的响应的关系, 在门水平上, 选取丰度, 种类数和优势度排名前3的蓝藻门, 裸藻门和绿藻门与环境因子进行Mantel相关性分析[图 4(a)], 结果发现蓝藻门水华藻类与TP、TN和NH4+-N显著正相关.在种水平上, 将全部20种优势水华藻类分别对关键环境因子做Pearson相关性分析[图 4(b)], 结果发现硅藻门和绿藻门所有优势种均与TP、TN和NH4+-N呈现不同程度的显著关系.以上结果表明, 氮磷等营养元素是影响水华藻类群落差异的关键因子, 且对不同水华藻类的影响程度呈现出差异.

1.TN, 2.TP, 3.NH4+-N, 4.T, 5.DO, 6.pH, 7.ORP, 8.COND, 9.NTU; *表示P<0.05, **表示P<0.01, ***表示P<0.001 图 4 水华藻类与环境因子相关性分析 Fig. 4 Correlation analysis between algal blooms and environmental factors

3 讨论 3.1 白洋淀水体环境与水华风险区域评估

白洋淀作为雄安新区生态建设的重要依托, 由于人类活动产生的污染物长期积累导致其治理措施效果不明显, 多处点位水质仍为Ⅳ~劣Ⅴ类水平, 不能达到水功能区划的要求[28~30].水质主要受降雨和水文循环等自然因素及生产生活和人工调水等人为因素的影响[31].夏季强降雨和高温导致此时淀区水质最差[32], 淀区环境因子PCA排序分析显示, 营养盐浓度对淀区内水质差异贡献最大, 且淀区水体环境存在空间异质性.以生境类型来看, DM的ρ(TN)和SD的ρ(TP)都处于较高水平, 平均值分别为1.37 mg·L-1和0.07 mg·L-1.以区域类型来看, N、H和WS区是营养盐平均浓度排名前3的区域.

本研究选取2020年夏季利用藻密度调查评估淀区水华风险区域, 结果更具有代表性.结果显示2020年夏季白洋淀95%以上地区无水华风险(藻类密度<2×106个·L-1), 仅在DM(点位编号:DM135, 位于N区; DM166, 位于H区)与SD(点位编号:SD224, 位于N区; SD17和SD245, 位于WS区)共5个点位检测到藻类密度为2×106~1×107个·L-1, 达到轻微水华程度.DM风险点位具体位于南刘庄(DM135)和光淀张庄(DM166), SD风险点具体位于小白洋淀(SD224和SD245)和烧车淀(SD17), 相关研究表明以上片区氮磷污染严重, 其中南刘庄水质最差, 属于重度富营养状态[33~35].

3.2 白洋淀水华藻类群落结构差异及其驱动因子

基于宏条形码分析, 共计得到90种白洋淀潜在水华藻类, 隶属于7个门, 以绿藻门和蓝藻门种类为主, 占比达60%.蓝藻和绿藻在温暖环境中更易繁殖, 同时绿藻门在富营养化环境中也更具优势[36, 37], 因此绿藻门物种丰度和数量在夏季白洋淀水华藻类群落中占据优势.通过优势度确定了淀区不同分区内20种水华优势藻种, 以绿藻门和裸藻门为主.这与针对已经暴发水华的淡水系统藻类研究有区别, Xin等[38]发现汉江中下游水华优势种为小环藻属(Cyclotella sp.), 还有一些相关研究指出蓝藻是水华暴发优势群落[39, 40], 针对白洋淀区域浮游植物优势种的调查研究也与本研究结果呈现出差异[41].造成上述差异除了与地理气候等自然因素有关, 也可能是研究对象(水华藻类群落与浮游植物群落)的不同.

结合群落α多样性和ANOSIM相似性分析发现, 与不同生境类型相比, 水华藻类群落差异在不同区域上更显著(图 2), 同时, 白洋淀水环境因子在不同区域也具有显著的空间异质性(表 2).淀区的水华藻类群落和水体环境都展现出较强的空间异质性, 说明淀区环境复杂, 且不同区域群落结构不同, 需分别关注和分析.RDA分析探究白洋淀水华藻类群落和水环境因子之间的关系, 筛选出TP、TN和NH4+-N这3个环境因子是造成水华藻类群落结构差异的关键因子(P<0.05), 且3个环境因子之间呈极显著正相关(图 3).国内学者研究白洋淀浮游植物和环境因子关系, 发现水温(T)、pH和溶解氧(DO)也会不同程度地影响浮游植物群落[42, 43].而本研究中氮磷等营养元素是影响水华藻类群落差异的关键因子, 说明特定浮游植物群落与环境响应关系也更为集中和具体.其他湖泊也有类似结果, Wang等[44]在研究滇池北部赤潮时, 发现TN和TP是关键因素.Tong等[45]发现湖泊营养盐浓度年内剧烈变化, 为蓝藻在一年内持续更长时间的水华创造了有利条件.基于水华藻类群落门水平和优势种水平分别与环境因子做相关性分析发现(图 4), 不同分类水平对关键环境因子的响应规律不同, 门水平上, 蓝藻与营养盐浓度响显著正相关; 种水平上, 则呈现出硅藻门和绿藻门优势种与营养盐浓度响应更显著, 可能是与高级分类单元相比, 低级分类单元情况更为复杂, 涉及到个体物种间的差异和相互作用等.

基于以上对白洋淀水华藻类群落的调查分析和水华暴发风险区域评估, 结合白洋淀实际情况, 提出以下水华防控意见:建议重点控制淀面和水道等宽阔水域的营养元素的输入, 降低水体富营养化程度; 同时保证白洋淀流域水华风险区域的基础水位, 有助于维持白洋淀高生物多样性, 降低水华暴发风险.

4 结论

(1) 基于藻密度评估白洋淀水华风险区域显示, 2020年夏季白洋淀95%以上地区无水华风险(藻类密度<2×106个·L-1), 仅5个样点存在轻微藻华风险:DM135, 位于N区; DM166, 位于H区; SD224, 位于N; SD17和SD245, 位于WS区.上述区域水质差, 氮磷污染严重, 建议控制白洋淀淀面和水道等宽阔水域氮磷等营养盐的输入, 监管维持水华风险区域的基础水位, 抑制水华藻类生长, 防止水华暴发.

(2) 通过宏条形码技术白洋淀淀区共检测到90种水华藻类, 隶属于7个门, 以绿藻门和蓝藻门为主, 其中优势水华藻类20种.淀区的水华藻类群落和水体环境均在不同区域上具有更强的空间异质性.TP、TN和NH4+-N这3个环境因子是造成水华藻类群落结构差异的关键因子, 水华藻类群落对关键环境因子的响应程度和规律在门分类和种分类水平上不同.

参考文献
[1] 王菁晗, 何吕奇姝, 杨成, 等. 太湖、巢湖、滇池水华与相关气象、水质因子及其响应的比较(1981—2015年)[J]. 湖泊科学, 2018, 30(4): 897-906.
Wang J H, He L Q S, Yang C, et al. Comparison of algal bloom related meteorological and water quality factors and algal bloom conditions among Lakes Taihu, Chaohu, and Dianchi(1981-2015)[J]. Journal of Lake Sciences, 2018, 30(4): 897-906.
[2] 张虎军, 宋挺, 朱冰川, 等. 太湖蓝藻水华暴发程度年度预测[J]. 中国环境监测, 2022, 38(1): 157-164.
Zhang H J, Song T, Zhu B C, et al. Annual forecast of the extent of cyanobacteria bloom in Taihu lake[J]. Environmental Monitoring in China, 2022, 38(1): 157-164. DOI:10.19316/j.issn.1002-6002.2022.01.15
[3] Health Canada. Guidelines for Canadian recreational water quality(3rd ed.)[S]. Ottawa, Canada: Health Canada, 2012.
[4] World Health Organization. Guidelines for safe recreational water environments. Volume 1, coastal and fresh waters[M]. France: World Health Organization, 2003.
[5] Yu P X, Gao R, Zhang D Z, et al. Predicting coastal algal blooms with environmental factors by machine learning methods[J]. Ecological Indicators, 2021, 123. DOI:10.1016/J.ECOLIND.2020.107334
[6] Wang C, Wang Z Y, Wang P F, et al. Multiple effects of environmental factors on algal growth and nutrient thresholds for harmful algal blooms: application of response surface methodology[J]. Environmental Modeling & Assessment, 2016, 21(2): 247-259.
[7] Xia R, Zhang Y A, Wang G S, et al. Multi-factor identification and modelling analyses for managing large river algal blooms[J]. Environmental Pollution, 2019, 254. DOI:10.1016/j.envpol.2019.113056
[8] Zhou Z X, Yu R C, Zhou M J. Resolving the complex relationship between harmful algal blooms and environmental factors in the coastal waters adjacent to the Changjiang River estuary[J]. Harmful Algae, 2017, 62: 60-72. DOI:10.1016/j.hal.2016.12.006
[9] 陈修康. 水环境治理中藻华的危害及其防控技术[J]. 绿色环保建材, 2021(10): 27-28. DOI:10.16767/j.cnki.10-1213/tu.2021.10.014
[10] 张民, 阳振, 史小丽. 太湖蓝藻水华的扩张与驱动因素[J]. 湖泊科学, 2019, 31(2): 336-344.
Zhang M, Yang Z, Shi X L. Expansion and drivers of cyanobacterial blooms in Lake Taihu[J]. Journal of Lake Sciences, 2019, 31(2): 336-344.
[11] 李建, 尹炜, 贾海燕, 等. 汉江中下游硅藻水华研究进展与展望[J]. 水生态学杂志, 2020, 41(5): 136-144.
Li J, Yin W, Jia H Y, et al. Research progress on diatom blooms in the middle and lower Hanjiang River: review and advances[J]. Journal of Hydroecology, 2020, 41(5): 136-144.
[12] 孟溪. 藻类鉴定及分类方法研究进展[J]. 生物化工, 2019, 5(2): 102-104.
Meng X. The progress of the methods to identify and classify phytoplankton[J]. Biological Chemical Engineering, 2019, 5(2): 102-104.
[13] 张军毅, 孙蓓丽, 朱冰川, 等. 基于分子标记的藻类鉴定研究进展[J]. 湖泊科学, 2021, 33(6): 1607-1625.
Zhang J Y, Sun B L, Zhu B C, et al. Research progress in the taxonomic identification of algae on the basis of molecular markers[J]. Journal of Lake Sciences, 2021, 33(6): 1607-1625.
[14] Gao Y C, Fang H D, Dong Y H, et al. An improved method for the molecular identification of single dinoflagellate cysts[J]. PeerJ, 2017, 5. DOI:10.7717/peerj.3224
[15] Dzhembekova Ni, Moncheva S, Ivanova P, et al. Biodiversity of phytoplankton cyst assemblages in surface sediments of the Black Sea based on metabarcoding[J]. Biotechnology & Biotechnological Equipment, 2018, 32(6): 1507-1513.
[16] Gao Y C, Chen Y Y, Xiong W, et al. Distribution patterns of dinoflagellate communities along the Songhua River[J]. PeerJ, 2019, 7. DOI:10.7717/peerj.6733
[17] Hanžek N, Udovič M G, Kajan K, et al. Assessing ecological status in karstic lakes through the integration of phytoplankton functional groups, morphological approach and environmental DNA metabarcoding[J]. Ecological Indicators, 2021, 131. DOI:10.1016/j.ecolind.2021.108166
[18] 杨苗, 龚家国, 赵勇, 等. 白洋淀区域景观格局动态变化及趋势分析[J]. 生态学报, 2020, 40(20): 7165-7174.
Yang M, Gong J G, Zhao Y, et al. Analysis of dynamic changes and trends in the landscape pattern of the Baiyangdian Region[J]. Acta Ecologica Sinica, 2020, 40(20): 7165-7174.
[19] 刘世存, 王欢欢, 田凯, 等. 白洋淀生态环境变化及影响因素分析[J]. 农业环境科学学报, 2020, 39(5): 1060-1069.
Liu S C, Wang H H, Tian K, et al. Ecological environmental changes and influencing analysis in Baiyangdian Lake[J]. Journal of Agro-Environment Science, 2020, 39(5): 1060-1069.
[20] 段茂庆, 许维, 苑飞燕, 等. 白洋淀近5年水生态环境质量变化趋势与营养状态分析[J]. 环境化学, 2022, 41(6): 1988-2000.
Duan M Q, Xu W, Yuan F Y, et al. Analysis on the trend and nutritional status of water environment quality in Baiyangdian Lake in recent 5 years[J]. Environmental Chemistry, 2022, 41(6): 1988-2000.
[21] Esenkulova S, Sutherland B J G, Tabata A, et al. Comparing metabarcoding and morphological approaches to identify phytoplankton taxa associated with harmful algal blooms[J]. Facets, 2020, 5(1): 784-811.
[22] Yang Y Z, Gao Y C, Chen Y Y, et al. Interactome-based abiotic and biotic impacts on biodiversity of plankton communities in disturbed wetlands[J]. Diversity and Distributions, 2019, 25(9): 1416-1428.
[23] Machado K B, Antunes A M, Targueta C P, et al. DNA metabarcoding reveals the responses of prokaryotes and eukaryotes microbiota to warming: Are the patterns similar between taxonomic and trophic groups?[J]. Ecological Indicators, 2020, 115. DOI:10.1016/j.ecolind.2020.106452
[24] Oksanen J, Kindt R, Legendre P, et al. The vegan package, community ecology package[M]. 2010.
[25] Clarke K R, Gorley R N. Primer v5: user manual/tutor3[M]. Plymouth: Primer-E Ltd, 2001.
[26] Ludwig J A, Reynolds J F. Statistical ecology: a primer on methods & computing[M]. New York: John Wiley & Sons, Inc., 1988.
[27] Pielou E C. Ecological diversity[M]. New York: John Wiley & Sons, 1975.
[28] 任晓强, 管孝艳, 陶园, 等. 白洋淀流域水环境风险评估综述[J]. 中国农村水利水电, 2021(1): 22-27.
Ren X Q, Guan X Y, Tao Y, et al. Review of water environment risk assessment in Baiyangdian Basin[J]. China Rural Water and Hydropower, 2021(1): 22-27.
[29] 王淑娟, 龚艳伟, 刘新军, 等. 白洋淀流域水体新污染物研究现状综述[J]. 环境与发展, 2022, 34(8): 102-108.
Wang S J, Gong Y W, Liu X J, et al. Review on the research status of emerging contaminants in Baiyangdian Lake[J]. Environment and Development, 2022, 34(8): 102-108.
[30] 高秋生, 焦立新, 杨柳, 等. 白洋淀典型持久性有机污染物污染特征与风险评估[J]. 环境科学, 2018, 39(4): 1616-1627.
Gao Q S, Jiao L X, Yang L, et al. Occurrence and ecological risk assessment of typical persistent organic pollutants in Baiyangdian Lake[J]. Environmental Science, 2018, 39(4): 1616-1627.
[31] Galbraith L M, Burns C W. Linking land-use, water body type and water quality in southern New Zealand[J]. Landscape Ecology, 2007, 22(2): 231-241.
[32] 王欢欢, 白洁, 刘世存, 等. 白洋淀近30年水质时空变化特征[J]. 农业环境科学学报, 2020, 39(5): 1051-1059.
Wang H H, Bai J, Liu S C, et al. Spatial and temporal variations in the water quality of Baiyangdian Lake in the recent 30 years[J]. Journal of Agro-Environment Science, 2020, 39(5): 1051-1059.
[33] 王子铭, 杨丽虎, 宋献方. 雄安建设初期白洋淀水质时空差异及其对降雨和土地利用的响应[J]. 环境科学, 2023, 44(7): 3820-3834.
Wang Z M, Yang L H, Song X F. Spatial-temporal variation of water quality and its response to precipitation and land use in Baiyangdian Lake in the early stage of the construction of Xiong'an new area[J]. Environmental Science, 2023, 44(7): 3820-3834.
[34] 赵黔伟, 汪敬忠, 魏浩, 等. 白洋淀水体氮、磷和重金属的变化特征及影响因素[J]. 水资源与水工程学报, 2020, 31(6): 103-108.
Zhao Q W, Wang J Z, Wei H, et al. Variation characteristics of nitrogen, phosphorus and heavy metals in waters from Baiyangdian Lake and the influencing factors[J]. Journal of Water Resources & Water Engineering, 2020, 31(6): 103-108.
[35] 朱金峰, 周艺, 王世新, 等. 白洋淀湿地生态功能评价及分区[J]. 生态学报, 2020, 40(2): 459-472.
Zhu J F, Zhou Y, Wang S X, et al. Ecological function evaluation and regionalization in Baiyangdian Wetland[J]. Acta Ecologica Sinica, 2020, 40(2): 459-472.
[36] 侯秀丽, 苑春刚, 李学平, 等. 滇池氮磷浓度变化对蓝、绿、硅藻年际变化的影响[J]. 水生态学杂志, 2018, 39(1): 16-22.
Hou X L, Yuan C G, Li X P, et al. Effect of nitrogen and phosphorus concentrations on the planktonic algae dynamics in Dianchi Lake[J]. Journal of Hydroecology, 2018, 39(1): 16-22.
[37] Jensen J P, Jeppesen E, Olrik K, et al. Impact of nutrients and physical factors on the shift from cyanobacterial to chlorophyte dominance in shallow Danish lakes[J]. Canadian Journal of Fisheries and Aquatic Sciences, 1994, 51(8): 1692-1699.
[38] Xin X K, Zhang H, Lei P, et al. Algal blooms in the middle and lower Han River: characteristics, early warning and prevention[J]. Science of the Total Environment, 2020, 706. DOI:10.1016/j.scitotenv.2019.135293
[39] Kim J, Jones J R, Seo D. Factors affecting harmful algal bloom occurrence in a river with regulated hydrology[J]. Journal of Hydrology: Regional Studies, 2021, 33. DOI:10.1016/j.ejrh.2020.100769
[40] Zhu C M, Zhang J Y, Wang X, et al. Responses of cyanobacterial aggregate microbial communities to algal blooms[J]. Water Research, 2021, 196. DOI:10.1016/j.watres.2021.117014
[41] 陈家琪, 董丽, 麻晓梅, 等. 基于eDNA技术的白洋淀微型生物群落监测[J]. 农业环境科学学报, 2021, 40(8): 1773-1786.
Chen J Q, Dong L, Ma X M, et al. Microbial community monitoring in Baiyangdian Lake based on eDNA technology[J]. Journal of Agro-Environment Science, 2021, 40(8): 1773-1786.
[42] 李娜, 周绪申, 孙博闻, 等. 白洋淀浮游植物群落的时空变化及其与环境因子的关系[J]. 湖泊科学, 2020, 32(3): 772-783.
Li N, Zhou X S, Sun B W, et al. Spatiotemporal variation of phytoplankton community and its relationship with environmental factors in the Lake Baiyangdian[J]. Journal of Lake Sciences, 2020, 32(3): 772-783.
[43] 金磊, 李林钰, 周杨, 等. 白洋淀三大典型水域浮游植物群落及水质评价[J]. 河北大学学报(自然科学版), 2017, 37(3): 329-336.
Jin L, Li L Y, Zhou Y, et al. Phytoplankton community and water quality analysis of three typical waters in Baiyangdian Lake[J]. Journal of Hebei University(Natural Science Edition), 2017, 37(3): 329-336.
[44] Wang J H, Li C, Xu Y P, et al. Identifying major contributors to algal blooms in Lake Dianchi by analyzing river-lake water quality correlations in the watershed[J]. Journal of Cleaner Production, 2021, 315. DOI:10.1016/j.jclepro.2021.128144
[45] Tong Y D, Xu X W, Qin M, et al. Lake warming intensifies the seasonal pattern of internal nutrient cycling in the eutrophic lake and potential impacts on algal blooms[J]. Water Research, 2021, 188. DOI:10.1016/j.watres.2020.116570