西北干旱区典型水库浮游植物群落结构特征及驱动因子*

杨宋琪，高兴亮，王丽娟，祖廷勋，王丹霞，罗光宏**

(1：河西学院，甘肃省微藻技术创新中心，甘肃省河西走廊特色资源利用重点实验室，张掖 734000)(2：西南大学，三峡库区生态环境教育部重点实验室，重庆 400715)

随着社会和经济的快速发展及人口的增加，河流、湖泊受到严重污染甚至退化，人类饮用水、农业灌溉和水产养殖等越来越依赖水库供水[1-3]，尤其在自然湖泊较少的干旱地区，水库发挥着极其重要的功能. 在水库人为调度过程中，蓄水排水使其成为高度动态而又复杂的生态系统[4]，使得水生高等植物等初级生产者难以得到稳定的水环境进行生长繁殖，致使浮游植物演替成为水库优势类群[1].

浮游植物是水生态系统初级生产力的主要贡献者，其组成和多样性变化直接影响水生态系统的结构与功能，同时，因浮游植物对水环境因子变化非常敏感而被广泛应用于湖泊、水库的日常监测中[5]. Han等研究发现，水库因水位波动使浮游植物群落的组成和动态存在较大个体差异[6]，而氮、磷、溶解氧、温度等亦是水库浮游植物物种分布及群落演替的主要驱动因子[7-9]. 近年来，面源污染及生活污水的排放使得水库水体持续富营养化[10-11]，水华频发[12-13]，尤其水华蓝藻大量增殖后释放的藻毒素及异味物质会对农作物、鱼虾贝类及人体健康造成不可逆的损害[14-15]，致使水库水生态安全面临严重威胁，生态服务价值丧失[16].

我国西北干旱地区降雨稀少，天然湖泊稀缺，河流径流量小，因此，水库水资源在饮用水供应及农业灌溉方面起着不可替代的作用. 那么面临极端的气候条件和水文过程的影响，干旱区水库浮游植物组成和丰度会呈现怎样的特征呢？以及主导因素是什么？目前针对西北干旱区水库的类似研究并不多. 为此，笔者以位于干旱气候区河西走廊中段的张掖为例，对其境内分布的主要水库水体环境因子和浮游植物群落结构进行调查分析，综合评估水库水质状况，并进一步明确驱动水库浮游植物物种分布及群落结构变化的关键环境因子，从而为水库水生态系统精准的监测和评估及管理提供科学依据.

1 研究区域与方法

1.1 张掖自然概况及采样点设置

张掖地处甘肃省河西走廊中段，属温带大陆性干旱-半干旱气候区，年均气温6～8℃，年均降雨量104～328 mm，且降雨多集中于夏季. 2017年，对张掖范围内不同分布区域的8座典型水库进行水样的采集(表1，图1)，其中高海拔高度大于2000 m的水库3座，分别为李桥水库、海潮坝水库和鹦鸽嘴水库；低海拔区域水库5座，分别为小海子水库、马郡滩水库、刘家深湖水库、马尾湖水库和天城水库，其中小海子水库为主要的渔业养殖场，刘家深湖水库为螃蟹养殖区.

图1 张掖采样水库位置分布

表1 被调查水库的基本信息

1.2 样品的采集与分析

依据水库面积和周围地理特征在每个水库选择具有代表性的3个位点，分别于2017年3月、6月、9月中旬和11月底采样，浮游植物定性分析样品用25#浮游生物网采集表层水样并立即用Lugol溶液固定保存. 使用5 L有机玻璃采水器于水体表层0.5 m深处分采集水样装入1 L塑料瓶中，立即加入15 mL Lugol溶液进行固定并摇匀保存，样品带回实验室静置48 h后，用虹吸管移去上清液浓缩至30 mL进行计数. 计数时，将浓缩液充分混匀后立即取0.1 mL于浮游植物计数框计数，藻类的鉴定参考《中国淡水藻类》、《中国西藏硅藻》等文献[17-20].

1.3 水体理化指标的测定

1.4 多样性指数及优势度的计算

1.4.1 多样性指数 Shannon-Wiener多样性指数(H′)：

(1)

Simpson指数(D)：

(2)

Pielou均匀度指数(E)：

(3)

Margalef指数(d)：

d=(S-1)/lnN

(4)

式中，N为总个体数，pi为第i个种类密度与总密度的比值，S为藻类的物种数.

1.4.2 优势度 Manaughton优势度指数(Y)：

(5)

式中，Ni为样品中某种藻的数量，N为采集样品中浮游植物总数量，fi为样品中某浮游植物出现的频率. 当物种优势度大于0.02时，该种即为优势种群[22].

1.5 数据处理

聚类分析(Cluster)和多维尺度转换排序(non-metric multidimensional scaling，nMDS)分析使用Primer 5软件. 采用R软件对浮游植物密度及环境因子进行Spearman相关性分析和物种与环境的CCA分析. 所有物种在满足出现的频度>12.5%且物种在至少一个样点的相对密度≥1%的基础上[23]，再选择优势种进行排序，物种矩阵及环境数据(除pH值外)均经过lg(x+1)转换[24]. 对物种信息进行去趋势对应分析(DCA)，排序轴梯度长度(LGA)均大于3，因此选择基于单峰模型的典范对应分析(CCA).

2 结果分析

2.1 水体理化因子变化

表2 张掖8座典型水库环境因子年均值及变化范围*

2.2 浮游植物群落结构特征

张掖8座典型水库共计检出浮游植物8门106属294种(表3). 其中硅藻门36属137种，占总种类数的46.59%，绿藻门38属81种，占总种类数的27.55%，蓝藻门17属44种，占总种类数的14.97%，甲藻门、隐藻门、金藻门和黄藻门种类占比共计为10.89%. 各水库检出物种差异较大，其中K2和K3种类数最少，且硅藻门种类占比均达到65%以上；但K4、K6、K7和K8检出物种数量显著增加，均达到140种以上，同时蓝藻门和绿藻门物种数量比例也相应上升.

2.3 浮游植物密度大小及优势种

调查期间张掖8座典型水库浮游植物密度范围在0.3×104～4.38×107cells/L之间(图2). 前3个季度K6、K7和K8密度均显著高于其他水库(P<0.05)；而在冬季，K2和K4水库藻密度最高. 浮游植物的演替方面，春季除K1水库外，各水库优势种以硅藻门种类为主；夏、秋季，K2和K3水库为硅藻门种类占优，但K4、K6、K7和K8水库蓝藻门种类藻密度急剧增加；冬季，随着温度的降低，各水库硅藻门种类密度均达到最大.

图2 不同水库各季节浮游植物密度

优势种分析可知(表4)，4个季节中K1水库优势种演替明显，而K2全年则主要以冷水型硅藻门种类微小曲壳藻(Achnanthesminutissima)、延长等片藻细弱变种(Daitomaelongatumvar.tenuis)、弯曲真卵形藻(Eucocconeisflexella)和扇形藻(Meridiumsp.)等占据优势；K3水库中硅藻门尖针杆藻(Synedraacus)在4个季度均占据明显优势；K4水库在前三个季度优势种主要集中于蓝藻门、绿藻门，仅在冬季硅藻门的尖针杆藻、尖端菱形藻(Nitzschiaacula)占优；K5水库优势种主要为硅藻门的尖针杆藻、谷皮菱形藻(Nitzschiapalea)和金藻门的分歧锥囊藻(Dinobryondivergens)等；K6、K7和K8水库则呈现典型的春、冬季硅藻门占优势、夏、秋季蓝藻占优势的特点，包括史密斯微囊藻(Microcystissmithii)、惠氏微囊藻(Microcystiswesenbergii)、细小平裂藻(Merismopediaminima)、假鱼腥藻(Pseudanabaenasp.)等.

表4 张掖8座典型水库不同季节浮游植物优势种

2.4 浮游植物多样性指数

2017年度，张掖8座水库浮游植物多样性指数H′、D、E和d年均值分别为2.26、0.77、1.20和0.86(表5). 其中，K7和K8水库浮游植物的H′、D和d均显著高于其他6座水库(P<0.05)，E最高值则出现在K3水库.

表5 张掖8座典型水库浮游植物多样性指数

2.5 8座水库间浮游植物群落相似性分析

以各位点的浮游植物密度组成原始数据矩阵，运用Primer 5进行聚类分析(Cluster)和多维尺度转换排序(non-metric Multidimensional Scaling,nMDS)分析，使得浮游植物群落结构相似性较高的聚成一组. 为更直观地观察不同水库间浮游植物组成的相似关系，在nMDS序图上勾画了聚类组. 结果表明，8个典型水库可以分为3组，其中K2独立一组，K1和K5为一组，其余位点为一组(图3). 由nMDS排序图也可看出，多维标度分析结果与聚类分析结果一致.

图3 8座典型水库聚类和多维尺度转换排序图

图4 浮游植物总密度与环境因子的Spearman相关性分析

2.6 浮游植物与环境因子之间的相关性

图5 主要优势种与环境因子之间的典范对应分析(C1席藻(Phormidium sp.)，C2细小平裂藻(Merismopedia minima)，C3束丝藻(Aphanizomenon sp.)，C4史密斯微囊藻(Microcystis smithii)，C5棒胶藻(Rhabdogloea sp.)，C6假鱼腥藻(Pseudanabaena sp.)，C7惠氏微囊藻(Microcystis wesenbergii)，G1四尾栅藻(Scenedesmus quadricauda)，G2球衣藻(Chlamydomonas globosa)，D1舟形藻(Navicula sp.)，D2微小曲壳藻(A. minutissima)，D3延长等片藻细弱变种(Daitoma elongatum var. tenuis)，D4等片藻(Diatoma sp.)，D5扇形藻(Meridium sp.)，D6尖针杆藻(Synedra acus)，D7尖端菱形藻(Nitzschia acula)，D8谷皮菱形藻(Nitzschia palea)，P1飞燕角甲藻(Ceratium hirundinella)，Chr1分歧锥囊藻(Dinobryon divergens))

3 讨论

水体中浮游植物群落结构和动态是多个环境因子在时空序列上综合作用的结果[25]，一般而言，浮游植物随季节变化呈现一定的演替规律，即春、冬季硅藻和金藻为主，夏、秋季蓝藻和绿藻占优，而冬季浮游植物种类数和密度均较低[17]，K2水库全年硅藻占优除外，其余7座水库演替规律一致. Stomp等对美国540个湖泊调查发现，地理格局及海拔对浮游植物多样性及丰度产生显著影响[26]. 然而，水库受人为水位调度和污染物排放等影响，浮游植物群落的演替亦会发生显著改变[27-28]. 2017年度，对西北干旱气候区8座典型水库进行调查，共计检出浮游植物8门106属294种，以硅藻门、绿藻门和蓝藻门种类为主(占比分别为48.35%、26.64%和14.47%)，浮游植物种类及密度随季节和空间地理分布呈现显著差异(P<0.05). 聚类和nMDS分析结果表明，8座水库主要分为3组，其中处于祁连山区高海拔的海潮坝水库(K2)独立一组，水库水源主要来自高山冰雪融水，浮游植物主要以喜低温和贫营养环境的硅藻如微小曲壳藻、延长等片藻细弱变种和扇形藻等为主要优势种；其次，李桥水库(K1)和马郡滩水库(K5)主要以分歧锥囊藻、飞燕角甲藻、尖针杆藻和谷皮菱形藻等物种为主要优势种而聚为一组；其余5座水库聚为一组，其中刘家深湖水库(K6)、马尾湖水库(K7)和天城水库(K8)水库位于黑河张掖段下游，承载了黑河及周围农田面源污染汇入的氮、磷等营养物质，尤其在夏、秋季高水温环境中蓝藻门种类快速繁殖并占据明显优势，而K4受渔业养殖的影响，水体富营养化程度高.

浮游植物多样性指数的变化能够直观地评判水体营养状况及其对浮游植物群落结构所产生的影响[29]，水体受到污染越严重时，浮游植物多样性指数值越小[30]. 2017年张掖8座水库Shannon-Wiener指数、Simpson指数、Pielou均匀度指数和Margalef指数年均值分别为2.26、0.77、1.20和0.86，可知水库属β中污型. 然而，浮游植物多样性与水质关系较复杂，易受水体类型、计数方法和鉴定种类的影响[31]，因此，再结合况琪军等[29]利用浮游植物细胞密度和物种组成对水质的评价标准，由于李桥水库(K1)、海潮坝水库(K2)、鹦鸽嘴水库(K3)和马郡滩水库(K5)浮游植物密度年均值小于106cells/L，且群落结构主要以硅藻、甲藻和金藻门种类为主，因此可判定上述水体为贫营养；而其余4座水库藻密度年均值都大于106cells/L(刘家深湖水库和马尾湖水库超过107cells/L)，同时浮游植物以蓝藻、绿藻占优，因而判定水体为中营养-中富营养.

4 结论

1)2017年张掖8座典型水库共计检出浮游植物8门106属294种，以硅藻门、绿藻门和蓝藻门种类为主，占比分别为48.35%、26.64%和14.47%，浮游植物种类及密度随季节和空间地理分布呈现显著差异(P<0.05).

2)聚类和nMDS分析结果表明，8座典型水库主要分为3组，其中高海拔区域的海潮坝水库(K2)独立一组，以喜低温和贫营养环境的硅藻为主要优势种；第二组为较高海拔区域李桥水库(K1)和鹦鸽嘴水库(K5)，以分歧锥囊藻和谷皮菱形藻、尖针杆藻等为主要优势种；其余5座富营养化程度相对较高的水库聚为一组，尤其在夏、秋季，蓝藻门种类占据优势明显.

3)Spearman相关性分析和CCA分析发现，TN、Spc、T和海拔是水库浮游植物密度大小及优势种分布的主要影响因子，低海拔区域水库中细小平裂藻、史密斯微囊藻和惠氏微囊藻等水华蓝藻随着水温升高而占据绝对优势，尤其史密斯微囊藻密度在夏季超过水华阈值(107cells/L)，值得引起注意.

致谢：感谢中国科学院海洋研究所徐沙博士对本文写作的帮助.