南水北调中线干渠浮游植物群落时空格局及其决定因子*

2021-05-10张春梅朱宇轩宋高飞米武娟毕永红王树磊梁建奎尚鸣宇

湖泊科学 2021年3期

张春梅，朱宇轩，宋高飞，米武娟，毕永红**,王树磊，梁建奎，尚鸣宇

(1:大连海洋大学，大连 116023) (2:中国科学院水生生物研究所，淡水生态与生物技术国家重点实验室，武汉 430072) (3:南水北调中线干线工程建设管理局，北京 100038)

浮游植物是水域生态系统中的初级生产者和基础环节，能对水体生态环境的变化做出快速直接的响应，生态环境变化也会直接或者间接影响浮游植物的群落结构[1]；一般认为，浮游植物种类数量的变化以及时空格局是营养盐、温度、光照、水动力条件等时空变化的结果[2-4]. 近年来，对于浮游植物群落结构的时空演变过程和关键影响因子的相关报道较多，丹江口水库不同区域的浮游植物优势群在同一季节具有显著差异，且随着季节变化其群落结构的演变趋势也有所不同[5]. 三峡水库的调查表明浮游植物具有明显的季节演变和空间差异[6]. 不同水体浮游植物群落结构特征不同，浮游植物因其对不同水环境的敏感响应和特异适应性而成为水环境监测的重要生物指示物种.

图1 南水北调中线干渠沿程采样点分布 (S1陶岔渠首、S2沙河渡槽进口、S3鲁山落地槽、 S4穿黄工程南岸、S5穿黄工程北岸、S6漳河北、 S7古运河暗渠、S8西黑山进口闸、S9天津外环河、 S10惠南庄北拒马河、S11团城湖)Fig.1 Distribution of sampling sites along the main channel of the middle route of South-to-North Water Diversion Project

南水北调中线工程解决了河南、河北、北京、天津4个省市的水资源短缺问题，具有重要的战略意义，发挥了巨大效益. 中线干渠为混凝土硬化渠道，与沿岸河流不交叉，水体特性不同于湖泊、河流和水库等自然水体，是一个全新封闭的人工系统[7]. 全面通水以来，在营养盐、水文水动力条件和生态因子的作用下，浮游植物群落结构处于动态演变过程中，特定时段还出现过某一特定类群生物量快速增殖的现象，成为影响输水水质的潜在生态问题. 2005年对中线工程水源地的研究[8]表明，陶岔渠首处无论在种类与数量上硅藻均占据绝对优势，春、秋季硅藻和绿藻较多，夏季蓝藻较为丰富. 2017年对中线工程水源地的研究[9]表明，群落演替由春季的硅藻、蓝藻向夏季的硅藻、蓝藻和隐藻以及秋、冬季的硅藻、隐藻和蓝藻转变. 尽管南水北调中线工程水源地浮游植物有一些相关研究报道，但迄今为止针对南水北调中线干渠水体浮游植物群落结构的研究资料较为匮乏.

鉴于中线干渠独特的渠道结构、较低的营养水平以及持续扰动掺混的流动状况，与丹江口水库来水中浮游植物的生境不同，因此本研究开展浮游植物群落结构与环境因子的同步调查，主要目的是验证以下假设：在跨流域长距离输水过程中浮游植物不仅仅是随水流的被动运输，可能存在沿程增殖生长行为；且浮游植物群落结构的时空异质性是由多因子共同决定的. 通过揭示干渠浮游植物群落结构特征及其演变的关键驱动因子，为深入理解南水北调中线干渠输水系统及其科学运行管理提供支撑.

1 材料与方法

1.1 研究区域与采样方法

南水北调中线干渠(32.68°～39.99°N, 111.72°～116.27°E)全长1432 km，多年平均调水量95亿m3，属于长距离跨气候带调水工程，南北纵跨8个纬度，四季气温变化分明. 本研究采样时间为2018年8月-2019年5月，按照夏季、秋季、冬季、春季每个季度月末采集一次，共计4次采样. 自河南南阳市淅川县陶岔渠首开始，到北京市团城湖共设置11个采样点(图1).

1.2 样品采集与处理

使用25#浮游生物网采集定性样品，用适量4%甲醛固定，定量样品用5 L玻璃采样器取表层水样置于瓶中，加入1%鲁哥氏液现场固定，样品静置48 h，沉降并浓缩至30 mL，充分摇匀后吸取0.1 mL至浮游生物计数框，在Olympus CX23下进行计数[11]. 藻类鉴定方法参照文献[12].

1.3 数据处理

优势度计算公式为：

(1)

式中，N表示同一样品中浮游藻类总个体数，ni表示第i种浮游藻类个体数，Y表示优势度，fi表示第i种浮游植物在样品中出现的频度，Y>0.02的物种定义为优势种[13].

图形绘制在Graphad prism 8软件中完成，采用ArcGIS 10.2软件作为空间分析工具，对沿线干渠各采样点浮游植物群落结构特征及细胞密度进行时空格局分析，应用IBM SPSS Statistics 25软件进行方差及相关性分析，采用Wilcoxon秩和检验方法检验不同季节穿黄工程前后细胞密度的差异性，使用R“vegan”包进行非度量多维尺度NMDS与ANOSIM差异性分析，使用R“UpSet”包中Upset函数进行样点间共有物种的维恩图分析.

采用Canoco 4.5软件[14]对物种数据与环境因子数据进行排序分析，通过对物种数据去趋势对应分析(DCA)确定排序模型，环境变量解释显著性由Monte Carlo测试检验，P值小于0.05则可用此方法进行检验排序，特征值越大，环境因子对藻类结构的影响越大，P值越小，相关性越显著，结果用物种-环境因子关系的双序图表示.

2 结果

2.1 浮游植物群落特征

2.1.1 物种组成与优势种调查期间共鉴定到浮游植物7门83属145种，以硅藻为主，共有68种，占总种类数的比例为46.89%；其次为绿藻，有47种，所占比例为32.41%；蓝藻16种，所占比例为11.03%；隐藻、裸藻、金藻、甲藻种类较少，所占比例分别为3.45%、0.27%、0.27%和0.68%.

从季节上看，浮游植物种类数排序为春季>秋季>冬季>夏季. 夏季绿藻种类最多，硅藻、蓝藻次之，其他季节种类数排序为硅藻>绿藻>蓝藻. 空间上陶岔渠首处(S1)物种数最低，从南到北呈递增趋势，沿程各样点均存在硅藻，且物种数呈现先增加后降低的波动性变化，以陶岔渠首处最低，穿黄工程北岸(S5)最高(图2).

图2 物种组成分布示意图：(a)空间分布；(b)各季节门水平变化；(c)各季节空间分布Fig.2 Distribution proportion of species: (a) spatial distribution; (b) species on phylum level varies from seasons; (c) spatial distribution varies from seasons

多集合维恩图(图3)显示，各样点浮游植物的共有种类有19种，主要为硅藻，包括桥弯藻(Cymbellasp.)、极小曲丝藻(Achnanthidumminutissimum)、内丝藻(Encyonopsissp.)、小环藻(Cyclotellasp.)、针杆藻(Synedrasp.)等. 蓝藻主要为束丝藻(Aphanizomenonsp.)、微囊藻(Microcystissp.)、假鱼腥藻(Pseudanabaenasp.)，隐藻主要为隐藻(Cryptomonassp.)、斜结隐藻(Plagioselmissp.). 陶岔处物种数目最少，穿黄北物种数目最多. 相比较于陶岔渠首处的物种种类，其他样点共有的新增加物种有4种，为空星藻(Coelastrumsphaericum)、维里纳优美藻(Delicataverena)、单角盘星藻具孔变种(Pediastrumsimplexvar.duodenarium)及瞳孔舟形藻(Naviculapupula). 各样点特有物种数分别为2、2、2、3、6、1、1、1、3、1、1种，特有物种较少，沿程物种变化较小.

图3 样点间浮游植物种类多集合维恩图(黑色点表示该位置有数据，点线图表示样点间存在交集，柱状图上的数字表示样点间的共有物种数，左边柱形图表示各样点的物种数)Fig.3 Multi-set venn diagram of phytoplankton species between sampling sites(Black dots indicated the existence of data at the sites, the dot plots indicated the intersection of the different sites, the numbers on the bar chart indicate the number of species in common between different sites, the bar chart on the left showed the number of species at each site)

浮游植物优势种共19种(表1)，不同季节优势类群存在差异：夏季以硅藻、绿藻、蓝藻为优势类群，空星藻优势度最大. 秋、冬、春季以硅藻、绿藻、蓝藻、隐藻为优势类群，分别以针杆藻、尖针杆藻(Synedraacus)、桥弯藻、斜结隐藻占据较大优势. 桥弯藻和极小曲丝藻为全年优势种.

2.1.2 群落的时空格局不同季节浮游植物群落组成存在一定差异(图4)，细胞密度主要是由硅藻、绿藻、蓝藻贡献，硅藻在全年占据优势. 春季以硅藻为主，其相对丰度沿程递增，隐藻、金藻主要出现在渠首(S1). 夏季绿藻相对丰度较高，蓝藻、硅藻次之，隐藻主要出现在渠首(S1)，穿黄工程前(S1～S3)与穿黄工程后(S5～S11)均以硅藻、绿藻为主，但所占比例存在差异. 秋、冬季为硅藻>绿藻>蓝藻，与夏季相比，蓝藻比例下降，隐藻的相对丰度在穿黄工程前各样点增加，冬季硅藻与绿藻相对丰度沿程大体呈相反的变化. 陶岔渠首处以硅藻、绿藻、隐藻及蓝藻为主要优势类群，冬季蓝藻所占比例较小.

表1 南水北调中线干渠浮游植物优势种及优势度

2.1.3 细胞密度与叶绿素a浓度浮游植物细胞密度的季节变化范围为2.16×106～1.92×107cells/L，平均密度为(6.79±7.15)×106cells/L，其中夏季>春季>秋季>冬季，单因素方差分析显示夏季显著高于其他季节(P<0.05)，其他各季节间无显著差异(图5a). 从空间尺度上分析，沿程各样点浮游植物细胞密度呈现先增加后降低、再增加再降低的“M”型变化趋势，穿黄工程南北(S4～S5)与古运河暗渠(S7)处均有较高值，春、夏季穿黄工程后显著高于穿黄工程前，夏、秋季古运河暗渠处均存在峰值(图5b).

Chl.a浓度变化范围为1.09～8.85 μg/L，平均值为(4.09±2.89)μg/L，其中夏季>秋季>春季>冬季，单因素方差分析显示夏季显著高于其他季节，与细胞密度变化趋势一致(图5a).

2.2 环境因子

2.3 群落结构与环境因子的关系

环境因子与样点的RDA分析结果显示(图6)，前两个排序轴累计解释了54.60%的物种变异和74%累计解释方差. 所有样本按照季节排序，分为春季、夏季、秋季、冬季. 确认水温和营养盐是影响浮游植物的主要因子，其次是流量和浊度.

对主要优势类群蓝藻、绿藻与硅藻进行多元线性分析(表3)可知，硅藻相对丰度与水温呈显著负相关，绿藻的相对丰度与水温呈显著正相关，蓝藻细胞密度与流速呈显著负相关，绿藻与硅藻的细胞密度与 NH3-N 浓度、流量呈显著负相关.

表2 南水北调中线干渠全年水环境因子平均值

图5 浮游植物细胞密度与叶绿素a浓度的时空格局Fig.5 Spatial and temporal changes of phytoplankton cell density and chlorophyll-a concentration

图6 样方与环境因子的RDA二维排序图Fig.6 Biplot diagram for RDA of the relationship between environment variable and sample sites

表3 多元线性回归分析

2.4 时空异质性

NMDS分析显示所有样本按照季节分为4组(图7a)，有明显的季节变化. ANOSIM结果显示(图7b)组间差异显著大于组内(P<0.01，R=0.767)，即浮游植物群落结构组成在时间上的差异显著大于空间上的差异，且不同季节间两两比较均存在显著差异(P<0.01).

基于UPGMA算法构建聚类树状结构，结果如图8所示：距离相近的样点群落组成相似性较高，渠首处(S1)与相邻样点间群落相似性较低. 穿黄工程前(S1～S3)与穿黄工程后(S5～S11)可大致分为两组，Wilcoxon秩和检验显示：春季穿黄工程前与穿黄工程后的绿藻、硅藻的细胞密度存在显著差异(P<0.05)，穿黄工程后显著高于穿黄工程前，夏季穿黄工程后的蓝藻、绿藻、硅藻细胞密度均显著高于穿黄前. 秋季穿黄后的硅藻细胞密度显著高于穿黄前. 冬季绿藻、金藻的细胞浓度存在显著差异.

图7 浮游植物群落组成的时间变异分析：(a)NMDS ；(b)ANOSIMFig.7 Temporal variability analysis of phytoplankton community composition: (a) NMDS; (b) ANOSIM

图8 基于Bray-curtis距离矩阵的聚类样点图Fig.8 Bray-curtis dissimilarity based on distance matrix showing the clustering of sampling sites

3 讨论

3.1 浮游植物群落结构及其演变

浮游植物季节演变的一般规律为：春、秋季以喜低温的硅藻、金藻为主，夏季以喜高温的蓝藻[15]和绿藻为主，冬季藻类种类与数量都相对较少[12]. 千岛湖[16]冬、春季浮游植物主要类群为硅藻、隐藻，夏季为绿藻、硅藻，秋季为蓝藻、绿藻和硅藻. 富春江水库[17]春季绿藻、隐藻、硅藻占优，夏季蓝藻及硅藻占优，秋季硅藻和隐藻占据较大优势. 中线干渠浮游植物群落季节变化与千岛湖及富春江类似，全年以硅藻为主，春季为硅藻型，夏季为绿藻-蓝藻-硅藻型，秋、冬季为硅藻-绿藻-蓝藻型. 这些水体均属于清洁水体，群落特征与水体类型具有明显的对应关系. 另外，不同水体的浮游植物群落因环境过滤作用而显示明显的差异，黑河张掖段的上游营养水平较低，主要优势类群为普生性的硅藻清洁指示种，下游营养水平较高，主要优势类群为蓝、绿藻中的中-富营养型指示种[18]. 广州地表水浮游植物的演变显示，休耕期间营养水平较低，硅藻为主要优势类群，播种期间营养水平高，蓝、绿藻取代硅藻成为主要优势类群[19]. 中线干渠全年水质较好，所有断面均可达到Ⅱ类及以上标准，营养水平为贫-中营养型，硅藻种类最多，为68种，春、秋、冬季均以硅藻种类为主要类群，这种以硅藻为主的浮游植物群落是水体营养水平较低的特征之一.

从NMDS与ANOSIM分析中(图7)可以看出，浮游植物群落结构变化的时间效应比空间效应更加明显，这可能是渠道内水体连通性以及快速流动导致的有效掺混削弱了空间距离对浮游植物群落组成的影响，而南北两地的显著季节差异和温度特征对浮游植物群落的影响更明显.

3.2 环境因子对浮游植物群落的影响

大量研究表明，水温、营养物质、水动力条件等多种因素可共同影响浮游植物群落的结构和功能[22]；南水北调中线干渠作为输水渠道，并不是天然河流，调水过程每年按照相关计划进行，不同时段的调水规程存在差异，本研究的全年监测数据经RDA与多元线性回归分析均确认水温、TP浓度、NH3-N浓度、流速和流量是影响中线干渠浮游植物群落结构的重要因子，该研究结果与以往其他水体的研究类似，说明中线干渠具有某些自然水体的属性.

水温对浮游植物种类和季节变化的重要影响[23]主要是由于浮游植物细胞的诸多代谢过程为酶促反应，这些生化反应过程受酶活性的影响，而温度对酶活性有决定作用. 水温升高加速藻类种群的生长增殖，且不同藻类在水温适应性上存在物种差异，这导致优势种随着水温变化而变化，水温通过调控细胞活性近而成为浮游植物群落结构变化的主要控制因子[24]. 研究表明蓝、绿藻的最适温度要高于硅、隐藻[23]. 本研究中干渠水温为5.10～34.30℃，四季变化明显且存在显著差异，RDA结果显示水温是影响群落结构的关键因子，多元线性回归分析显示水温与硅藻的相对丰度呈显著负相关，与蓝藻、绿藻的细胞密度呈显著正相关，低温条件下，蓝绿藻相对丰度下降，硅藻相对丰度占据优势. 如空星藻的最适生长温度为30～35℃[25]，干渠夏季水温较高，适合空星藻生长，因此它是夏季的优势种之一；隐藻常被认为在水温较低时种群密度较高[26]，干渠冬、春季水温较低，这使得隐藻成为冬、春季的优势种. 因此水温影响和决定着干渠浮游植物细胞密度以及优势类群的演变.

营养物质是浮游植物生长繁殖必需的，也是浮游植物季节变化的重要环境影响因素[27]. TN、TP是影响浮游植物细胞密度和群落结构的重要因子，有研究显示早春时节藻类大量增殖的原因主要是对氮的快速吸收利用[22]；本研究也显示NH3-N浓度与细胞密度呈显著负相关，表明浮游植物生长会消耗大量NH3-N. 同时，由于硅藻生长的最适磷浓度较低(0.002～0.01 mg/L)，随着水体中磷浓度的降低，蓝、绿藻相对丰度降低，硅藻逐渐占据优势[28]. 干渠TN浓度为0.39～2.13 mg/L，TP浓度为0～0.01 mg/L，氮浓度高，磷浓度较低，N/P比值远大于16，表现为磷限制，这种营养条件更适合硅藻的生长.

流量、流速可显著影响藻类群落组成以及优势类群的演变[29]. 一般来说，蓝藻和绿藻较适合在静止或弱紊动水体中生长，硅藻在流动或紊动环境下可获得较大的增殖优势[30]. 中线干渠水体持续流动高度掺混，流速范围为0.1～2.4 m/s，平均流速为0.95 m/s，相比较于其他水体流速较大[31]. 本研究中，蓝藻、绿藻的细胞密度与流速呈显著负相关，硅藻的相对丰度与流速呈正相关，在中线干渠这种动态混合良好的生境中，高流速使得硅藻的优势地位增强，环境过滤与适者生存的选择作用更利于硅藻的增殖，特别是那些比表面积大且对低光有适应性的硅藻，如桥弯藻、舟形藻、曲丝藻等[32]. 极小曲丝藻具有良好的环境适应性，易在高流速水体中成为优势类群[33]，因此，它成为中线干渠的全年优势种. 本研究结果显示，流量与硅藻、绿藻的细胞密度呈显著负相关，流量加大可能会对水体中藻类密度具有一定的稀释作用. 对西江下游浮游植物的研究也表明，流量加大有利于种类丰富度的增加，但不利于细胞密度的增加[34]. 可见，流量加大在一定程度上可以改变藻类细胞密度，杜绝大量藻类短时间内快速生长增殖或衰老死亡对水体营养盐的急剧影响并可能导致水体异味物质含量升高等问题，可保障水质的稳定.