文章快速检索     高级检索
  生态与农村环境学报  2019, Vol. 35 Issue (9): 1172-1181   DOI: 10.19741/j.issn.1673-4831.2018.0784
0
骆马湖夏季浮游植物群落结构变化及其驱动因子分析
覃宝利 , 唐金玉 , 王宣朋 , 吴春 , 叶建勇 , 丁辰龙     
江苏省农业科学院宿迁农科所, 江苏 宿迁 223800
摘要:于2017年夏季对骆马湖浮游植物群落进行调查,探讨骆马湖浮游植物群落结构变化与环境因子的关系,以期为骆马湖生态保护提供科学依据。此次调查鉴定出浮游植物6门32属,其中,绿藻门属数最多,其后依次为蓝藻门和硅藻门,优势属主要为浮游蓝丝藻、微囊藻和小球藻。全湖浮游植物细胞丰度在2.63×105~2.85×107 cells·L-1之间,生物量在0.092~4.522 mg·L-1之间。全湖浮游植物Shannon-Wiener指数在0.60~2.60之间,平均值为1.75,且9月浮游植物Shannon-Wiener指数显著高于8月(P < 0.05),但不同点位之间多样性指数的差异未达显著水平(P>0.05)。与往期骆马湖调查结果相比,此次调查得到的全湖多样性指数差异不大,细胞丰度明显增加,优势属多集中在蓝藻门和绿藻门,尤其是浮游蓝丝藻明显增多。从近几年的优势属种类、藻类细胞密度变化来看,骆马湖富营养化程度依然在加剧。RDA分析结果表明,水温、溶解氧浓度和氮磷比的共同作用解释了骆马湖夏季浮游植物群落结构变化的28.16%。其中,骆马湖浮游植物群落结构的变化受采砂、围网养殖及夏季人类活动影响较大。因此,减少人为活动干扰对于保护骆马湖水质和南水北调工程水质势在必行。
关键词骆马湖    浮游植物    群落结构    驱动因子    
Analysis of Phytoplankton Community Structure Changes and the Driving Factors in Summer of Luoma Lake
QIN Bao-li , TANG Jin-yu , WANG Xuan-peng , WU Chun , YE Jian-yong , DING Chen-long     
Suqian Institute of Agricultural Sciences, Jiangsu Academy of Agricultural Sciences, Suqian 223800, China
Abstract: In order to provide a scientific basis for the ecological protection of Luoma Lake, northern Jiangsu Province, in the summer of 2017, the phytoplankton community in the lake was investigated to explore the relationship between phytoplankton community structure changes and environmental factors. A total of 6 phyla, including 32 genera of phytoplankton were identified. The number of Chlorophyta was the largest, followed by Cyanophyta and Diatoma, and the dominant genera were mainly Planktothrix, Microcystis and Chlorella. The phytoplankton abundance of the whole lake was ranged from 2.63×105 to 2.85×107 cells·L-1, and biomass was ranged from 0.092 to 4.522 mg·L-1. Shannon-Wiener index of phytoplankton in the whole lake ranged from 0.60 to 2.60, with an average of 1.75, moreover, the index in September was significantly higher than that in August (P < 0.05), but the difference in diversity index between different locations was no significant (P>0.05). Compared with the results of previous investigations in Luoma Lake, the diversity index of the whole lake did not have much difference, the cell abundance was significantly increased, and the dominant genera were concentrated in Cyanophyta and Chlorophyta, especially the number of Planktothrix increased distinctly. From the changes of dominant genera and algae cell density in recent years, the eutrophication degree of Luoma Lake has still been increasing. The results of RDA analysis show that the combination of water temperature, dissolved oxygen and nitrogen-phosphorus ratio explained 28.16% of the phytoplankton community structure change in the summer of Luoma Lake. Sand mining, purse seine culture and summer human activities had a comparatively large impact on this change. Therefore, it is imperative to reduce the disturbance of human activities to protect the water quality of Lake Luoma and the water quality of the South-to-North Water Transfer Project.
Key words: Luoma Lake    phytoplankton    community structure    driving factor    

骆马湖(34°00′~34°11′ N, 118°06′~118°18′ E)是淮河流域第3大湖泊, 江苏省第4大湖泊, 对南水北调东线工程具有重要的调蓄作用, 同时也是宿迁市的重要饮用水源地, 承担着宿迁市区供水功能[1]。骆马湖北面通过京杭大运河与山东南四湖相连, 南与洪泽湖相连, 继而与长江水系相通。入湖河流主要有沂河水系、南四湖水系和邳苍地区共40多条支流, 出流有3处, 分别经嶂山闸入新沂河, 经皂河闸入中运河, 经洋河滩闸入六塘河。据统计, 骆马湖年水位变幅为1.90~5.73 m, 年换水10次左右, 属典型过水性湖泊。近年来, 无序采砂、不合理养殖和入湖污染物增加等状况导致骆马湖湖泊生态系统遭到破坏, 水环境发生变化, 生物多样性降低[1], 使得骆马湖生态环境问题日益严重。

浮游植物作为湖泊生态系统中的初级生产者, 具有较短的生长周期和较强的扩散能力, 不同种类浮游植物对栖息地、营养物质和光照强度等的需求不同。种群密度对水质、水动力和气候条件等变化非常敏感, 故通常被用于指示水体环境变化, 特别是营养盐水平的变化[2]。欧盟水框架协议、海洋战略框架指令等国际政策法规[3-4]均提出将浮游植物群落作为水质评价的重要指标, 通过其生物量、群落组成和丰富度等指标来反映水体特征和水生态状况。水体的温度、光照、营养盐及水动力等环境因子又易影响浮游植物群落结构的变化[5]。因此, 进行浮游植物群落特征及其影响因子的研究, 对于了解整个湖泊生态系统的结构、功能及维护生态平衡具有重要意义。

近年来, 众多学者对骆马湖水质及富营养化程度进行研究, 认为骆马湖处于中富营养化水平, 这些研究大多数集中于富营养化评价等方面[6-7], 有关骆马湖浮游植物群落结构变化及其驱动因子的报道不多。该文旨在通过对骆马湖进行采样调查, 分析浮游植物群落结构变化及其驱动因子的影响, 以期为骆马湖生态环境保护与管理提供数据支持。

1 研究方法 1.1 水样的采集和鉴定

于2017年8—9月年进行野外调查, 每月初采样1次, 共设置20个采样点(图 1)。采用GPS(集思宝G120BD)现场测定样点位置, 采用YSI 85溶氧仪测定水体温度(T)和溶解氧(DO)浓度, 采用YSI60 pH计测定pH, 采用塞氏盘测定透明度。用5 L采水器采集表层、中层和底层水样, 现场混匀后取5 L水样冷藏带回实验室分析总氮(TN)和总磷(TP)浓度等[8]水质指标。叶绿素a和藻蓝蛋白(phycocyanin, PC)浓度采用丙酮提取法测定[9]。每个采样点水样混合后另外保留1 000 mL, 立即加入鲁哥试剂固定用于分析浮游植物样品, 在稳定的实验台上固定24~36 h, 用细小虹吸管移出上清液, 最后定容到30 mL。取浓缩后的0.1 mL样品利用放大400倍的显微镜鉴定浮游植物种类并计数。计数方法为目镜视野法, 一般随机计数30~50个视野, 使得细胞数在300以上。由于浮游植物的相对体积质量接近1, 故可直接将浮游植物的体积换算成生物量(以湿重计), 即生物量为各自丰度乘以各自平均体积, 单位为mg·L-1, 单细胞生物量主要根据浮游植物个体形状测量得到。上述实验方法主要参见《淡水浮游生物研究方法》[10]

图 1 骆马湖采样点位分布 Fig. 1 Spatial distribution of sampling sites in Luoma Lake
1.2 数据分析

利用Excel 2003、Origin 8.5和SPSS 12软件制图并分析。利用Canoco 4.5软件的冗余分析(redundacy analysis, RDA)计算浮游植物组成与环境因子的相关性, RDA分析中环境因子和种群组成数据均经过lg (x+1)转换[11], 利用软件中的预筛选功能(forward selection)筛选对群落组成影响显著的环境因子, 采用499次非限制性蒙特卡罗排列(Monte Carlo permutations)对计算结果的显著性进行验证, 排序结果用表征种类、样品和环境因子关系的三序图表示。

浮游植物优势种用优势度(Y)来判定, 其计算公式为

$ Y = \left( {{n_i}/N} \right)\cdot{f_i}。$ (1)

式(1)中, fi为物种i在采样点的出现频率; ni为第i个物种的密度; N为所有物种总个数。物种优势度Y≥0.02即定为优势种[12]

浮游植物多样性采用SPADE(species prediction and diversity estimation)软件进行计算, 多样性指数[13]包括Shannon-Wiener指数(H)、Simpson指数(M)和Pielou指数(J), 其中, Shannon-Wiener指数受物种丰富度影响较大, Simpson指数较多地反映物种优势度, Pielou指数则是物种均一性的度量。HMJ计算公式为

$ H = - \sum\limits_{i = 1}^n {{P_i}{\rm{ln}}{P_i}} , $ (2)
$ M = \sum\limits_{i = 1}^n {{P_i}^2} , $ (3)
$ J = - \sum\limits_{i = 1}^n {{P_i}{\rm{ln}}{P_i}/{\rm{ln}}S}。$ (4)

式(2)~(4)中, Pi为第i个浮游植物种数量; S为每个采样点的浮游植物物种总数。

2 结果与分析 2.1 采样点位

根据环境保护部门在骆马湖设立的省控监测点位[14]和邹伟等[15]对骆马湖区生境特征的划分, 采样点划分情况如下:9~12#为戴场(围网养殖), 6~8、13#为三场(采砂区), 4~5#为嶂山闸, 14~18#为骆马湖乡(围网养殖), 1~3、19~20#为洋河滩闸(船只停靠、渔货销售、生活及游乐场所)。其中, 9~12#区域水深为3.4~3.9 m, pH为8.1~8.7, 透明度为60~105 cm, COD为10.9~14.1 mg·L-1; 6~8、13#区域水深为4.8~5.4 m, pH为8.0~8.5, 透明度为52~75 cm, COD为11.0~14.6 mg·L-1; 4~5#区域水深为3.6~4.3 m, pH为8.0~8.3, 透明度为65~80 cm, COD为9.0~13.0 mg·L-1; 14~18#区域水深为3.8~4.1 m, pH为8.5~8.8, 透明度为60~110 cm, COD为8.9~15.4 mg·L-1; 1~3、19~20#区域水深为2.8~4.3 m, pH为8.3~8.7, 透明度为70~100 cm, COD为9.4~14.3 mg·L-1。1~3、19~20#区域底质水草较多, 6~8、13#区域底质为砂质, 其他区域底质大部分为淤泥或黏土。

2.2 浮游植物群落组成及分布

此次调查鉴定到包括绿藻门、蓝藻门、硅藻门、裸藻门、甲藻门和隐藻门共计32属。其中, 绿藻门11属, 蓝藻门7属, 硅藻门6属, 裸藻门5属, 甲藻门2属, 隐藻门1属。此次调查鉴定到的绿藻门属数最多, 优势属(括号内数字为优势度)有盘星藻(0.04)、小球藻(0.03)、十字藻(0.02)和栅藻(0.01);硅藻门优势属有直链藻(0.01);蓝藻门优势属有浮游蓝丝藻(0.59)、微囊藻(0.03)、鱼腥藻(0.02)、颤藻(0.01)和螺旋藻(0.01);隐藻门优势属有隐藻(0.01)等, 优势属点位分布见表 1

表 1 骆马湖夏季浮游藻类优势属空间分布 Table 1 Distribution of dominant phytoplankton species in summer in Luoma Lake

8—9月骆马湖全湖浮游植物丰度在2.63×105~2.85×107 cells·L-1之间, 平均为1.50×108 cells·L-1, 其中, 蓝藻门丰度占总丰度的82.2%, 为绝对优势门, 其后依次为绿藻门和硅藻门, 分别占15.3%和1.2%。9月浮游植物总丰度为8月的87.4倍, 蓝藻门总丰度是8月的9.5倍。从时间分布来看, 9月浮游植物总生物量显著高于8月(P < 0.05)(图 2), 8月绿藻门最多(11.803 mg·L-1), 9月蓝藻门最多(19.134 mg·L-1)(图 2)。

图 2 不同月份骆马湖浮游植物总生物量、主要门类生物量和优势物种生物量分布 Fig. 2 Phytoplankton distribution of total biomass, main species biomass and dominant species biomass in different months in Luoma Lake

8月小球藻生物量(0.613 mg·L-1)>浮游蓝丝藻(0.473 mg·L-1)>微囊藻(0.314 mg·L-1), 优势属间生物量相差不大, 但9月优势属生物量相差较大, 浮游蓝丝藻生物量最高(14.91 mg·L-1), 小球藻次之(2.882 mg·L-1), 微囊藻最少(0.074 mg·L-1)(图 2)。从空间分布来看, 8月全湖浮游植物分布差异较大, 其中, 3~15#点位绿藻门、硅藻门生物量较高, 2、16#点位微囊藻属生物量明显高于其他点位。9月绿藻门、硅藻门生物量分布差异不大, 蓝藻门生物量大体上为1~3、7~8、14~20#点位较高, 4~6、9~13#点位较低, 浮游蓝丝藻属生物量分布与之类似。

2.3 浮游植物多样性指数及NMDS分析

骆马湖8—9月浮游藻类Shannon-Wiener指数在0.60~2.60之间, 平均值为1.75, Simpson指数在0.26~0.78之间, 平均值为0.74, Pielou指数在0.37~0.96之间, 平均值为0.79(图 3)。

图 3 骆马湖夏季藻类多样性时空分布 Fig. 3 Temporal and spatial distribution of algal diversity in summer in Luoma Lake

9月Shannon-Wiener指数显著高于8月(P < 0.05), Simpson指数个别点位显著高于8月(P < 0.05), Pielou指数与8月相比无显著差异(P>0.05)。全湖各采样点位之间Shannon-Wiener指数、Simpson指数和Pielou指数均无显著差异(P>0.05), 但3者之间均呈显著相关(P < 0.05)。

图 4显示, 8、9月骆马湖各点位浮游藻类组成差异明显, 9月各点位藻类组成相似性较大, 分布较集中, 8月各点位藻类组成相似性更小, 分布较分散。

圆圈表示8月点位,三角形表示9月点位。 图 4 骆马湖夏季藻类NMDS样本得分排序图 Fig. 4 NMDS sample score sorting map in summer in Luoma Lake
2.4 浮游植物与环境因子的关系 2.4.1 骆马湖夏季主要环境因子差异分析

通过方差分析和多重比较发现, 骆马湖8月(28.9 ℃)、9月(25.4 ℃)平均水温有显著差异(P < 0.05), 但空间上无显著差异。ρ(DO)在时间和空间上均无显著差异。9月透明度(75.25±12.92)显著低于8月(88.45±22.96)(P < 0.05)。ρ(TN)平均值在8、9月之间无显著差异, 在空间上1~6#点位显著低于7~20#点位(P < 0.05)(图 5)。9月ρ(TP)平均值〔(0.38±0.12) mg·L-1〕显著高于8月〔(0.23±0.12) mg·L-1〕(P < 0.05), 在空间上无显著差异(图 5)。

图 5 夏季骆马湖不同点位总氮(TN)和总磷(TP)浓度 Fig. 5 Total nitrogen (TN) and total phosphorus (TP) contents at different points in summer in Luoma Lake
2.4.2 浮游植物藻蓝蛋白(PC)和TP浓度及总氮总磷比〔ρ(TN):ρ(TP)〕的变化关系

浮游植物PC与氮磷比和TP相关性见图 6

图 6 夏季骆马湖藻蓝蛋白与ρ(TN):ρ(TP)比值和TP浓度的相关性 Fig. 6 Correlation between the changes of phycocyanin and ρ(TN):ρ(TP) ratio, and ρ(TP) in summer in Luoma Lake

图 6所示, 8月全湖不同点位PC浓度分别与ρ(TN):ρ(TP)比值和TP浓度呈显著正相关(P < 0.05)和显著负相关(P < 0.01), 9月不同点位PC浓度分别与ρ(TN):ρ(TP)比值和TP浓度呈显著负相关(P < 0.01)和显著正相关(P < 0.05)。

2.5 浮游植物群落与环境因子的RDA分析

通过RDA分析(图 7)发现, 在选取的主要环境因子中, 经蒙特卡罗检验, 有3个因子与浮游植物分布显著相关(P < 0.05)。3个显著的环境因子(箭头表示)共同解释了浮游植物群落变化的28.16%。RDA1和RDA2分别解释了浮游植物群落变化的18.9%和5.1%。

图 7 骆马湖夏季环境因子与浮游植物组成的RDA排序图 Fig. 7 RDA sorting map of environmental factors and phytoplankton in summer in Luoma Lake
3 讨论 3.1 骆马湖浮游植物群落组成与生物量变化

自20世纪以来, 骆马湖富营养化程度正在加剧[6], 浮游植物丰度在逐年增加(表 2)。浮游植物细胞密度越高, 表明湖泊营养化程度越高。优势类群占比可以用于指示湖泊的污染状态[16], 如中等营养水平湖泊中蓝藻、绿藻、硅藻和隐藻都占有一定比例, 富营养型湖泊则常以绿藻、蓝藻类占优势[13]

表 2 近年来骆马湖浮游植物调查结果比较 Table 2 Comparison of phytoplankton survey of Luoma Lake in recent years

此次调查得到的浮游植物生物密度均高于以往调查结果(表 2[6-7,17-19]), 藻类优势属种类数也与以往调查结果不同, 如2014年骆马湖夏、秋两季优势属为裸藻、伪鱼腥藻、纤维藻和隐藻[6], 此次调查优势属为浮游蓝丝藻、微囊藻、小球藻和十字藻(Y≥0.02), 且丝状蓝藻明显增多。因此, 与以往调查结果中优势属和浮游植物细胞密度做对比, 结果表明骆马湖富营养化程度仍然在加剧。

骆马湖8月总生物量、主要门类生物量及优势属生物量分布没有明显规律。9月总生物量呈西南区域明显高于东北区域的趋势, 与杨士建[14]的调查结果相似。9月浮游植物总丰度为8月的87.4倍, 这主要是由于9月浮游蓝丝藻大量增殖所致, 且此次调查表明丝状蓝藻占比较以往调查结果明显增加。浮游蓝丝藻是水华蓝藻的重要常见种属。浮游蓝丝藻属(又称浮丝藻属)是1988年ANAGNOSTIDIS等从颤藻(Oscillatoria)中分离出来的, 以阿氏浮丝藻(Planktothrix agardhii)〔阿氏颤藻(Oscillatoria agardhii)〕作为模式种[20]。而在中国一直未见有关浮游蓝丝藻属的报道, 直到2003年陈宇炜等[21]介绍了包括Planktothrix的3个浮游蓝藻属的国际最新分类系统, 并将Planktothrix称为浮游蓝丝藻属。在欧洲的很多水体, 浮游蓝丝藻常年处于优势地位, 且产生很高毒素, 带来众多环境问题, 是目前欧洲学者蓝藻研究中的重点之一。关于浮游蓝丝藻的报道主要集中在欧洲, 且研究重点在其毒素方面[20]。近年来, 在我国各地也有发现, 很多水体中浮游蓝丝藻占据优势地位。目前, 这些水体的毒素监测报道很少。由于骆马湖是宿迁市的重要饮用水源地, 承担着宿迁市区供水功能。此次调查到的骆马湖中浮游蓝丝藻是否产毒, 还需后续研究。

3.2 骆马湖浮游植物群落多样性

通过方差分析和多重比较, 发现骆马湖夏季生物多样性指数不存在显著的空间差异(P>0.05)。9月骆马湖Shannon-Wiener指数显著高于8月(P < 0.05), 说明在夏季骆马湖生物多样性增加, 种类明显增多, 丰富度明显增加。此与NMDS分析结果一致。Simpson指数和Pielou指数均变化不大, 且两者相关系数为0.817, 说明骆马湖夏季全湖浮游植物种类分布比较均匀。Shannon-Wiener指数与Simpson指数和Pielou指数的相关系数分别为0.945和0.630, 由此可以推测, 骆马湖夏季物种多样性指数的变动, 可能更多地是由于种类数的变动而非均匀度的变化。

典型浮游植物的Shannon-Wiener指数基本在2.4~2.6之间[22]。此次调查中骆马湖Shannon-Wiener指数处于0.60~2.60之间, 说明有些点位浮游植物Shannon-Wiener指数明显偏低。该调查结果与彭凯等[6]对骆马湖的调查结果相似, 与田家怡等[23]于2005年的调查结果相比偏高, 与张小明等[7]于2015年对骆马湖北部的调查结果相比偏低, 可能与调查点位不同有关。与其他湖泊相比, 与邻近的南四湖[24]相似, 与云贵高原13个湖泊的Shannon-Wiener指数[13]相差较远, 比洪泽湖浮游植物的Shannon-Wiener指数[6]偏低。

3.3 环境因子对骆马湖浮游植物组成的影响

影响浮游藻类种群的因素有很多, 如水流速度、水温、营养盐、生物因素等。在河流、湖泊和水库中, 水力滞留时间的长短决定了浮游藻类的种类组成和生物量[18]。此次调查结果表明8月骆马湖湖区北部绿藻门、硅藻门占比较高, 其中, 绿藻门盘星藻、空球藻和实球藻较多, 硅藻门直链藻属较多, 而相对体积质量较小的蓝藻门占比较低, 彭凯等[8]也发现类似现象。骆马湖来水主要为上游的沂河水系, 皂河闸由于常年处于关闭状态, 出水以新沂河-嶂山闸为主, 南部地区水流相对较为平缓[8]。由于硅藻个体较大, 适合生长在水力滞留时间短、水流流速快的水体中, 而蓝藻适合在水力滞留时间长的水体中生长[18]。因此, 8月骆马湖藻类群落构成主要与水体水力滞留时间短和流速较大有关。9月主要是蓝藻门的分布差异较大, 其中, 1~3#(洋河闸)、7~8#(采砂区)、14~20#(围网养殖)点位蓝藻门生物量较高, 说明这些区域水体营养物质含量较高, 水质较差。采砂导致沉积物中蓄积的营养物质释放至水体中, 养殖投饵也会增加水体的营养负荷, 而洋河闸水质较差可能受捕鱼船只停靠、渔货销售生活区污染的影响; 4~6#(嶂山闸)、9~13#(围网养殖、采砂区)点位蓝藻门生物量较低。13#点位位于主采砂区, 其蓝藻生物量较低, 可能是采砂导致水体过于浑浊而抑制了包括蓝藻在内的所有藻类的生长, 其他点位较低可能是调查期间恰好处于开闸放水期, 带走了大量营养物质。骆马湖为典型的过水性湖泊, 其水质变化是对区域内人类活动强烈程度的反馈[25]。因此, 骆马湖现阶段水质状况与夏季人类活动加剧有关。

蓝藻中微囊藻明显占优势的点位为8月的2#和16#, 其他大部分区域均为浮游蓝丝藻占优势, 尤其是9月占绝对优势。8月TP浓度与PC浓度呈显著负相关, 2#和16#点位TP浓度在全湖中较低。而微囊藻具有较高的磷吸收最大摄取速率, 比其他藻类有更强的储存磷的能力[26]。因此, 2#和16#点位微囊藻占优势的结果可能是原有TP浓度较高所致。另外, 浮游动物对浮游藻类的捕食作用具有选择性, 不同类型藻类受捕食作用的影响是不同的。研究表明, 捕食性微小动物僧帽触口虫对孟氏浮游蓝丝藻有明显的控制作用[27]。由于僧帽触口虫对浮游蓝丝藻的捕食, 反而使微囊藻生长更旺盛。因此, 浮游动物牧食其他藻类而非微囊藻也可能是2#和16#点位微囊藻占优势的原因之一。微囊藻最佳生长温度高于其他藻类, 如太湖微囊藻最适生长温度为30~35 ℃[26]。金相灿等[28]研究表明浮游蓝丝藻比微囊藻更适合低温生长。浮游蓝丝藻属是从颤藻属中分离出来的新属, 很多特性与颤藻相似, 在水温约为25 ℃时最易暴发[29]。而骆马湖9月平均水温为25.4 ℃, 更适合浮游蓝丝藻的生长, 这可能是9月浮游蓝丝藻占绝对优势的原因之一。

磷是浮游植物生长繁殖所必需的营养元素, 是影响浮游植物群落组成的重要环境因子之一。湖泊的富营养化, 尤其是水体中磷浓度的增加, 通常会导致水体中浮游植物种群组成朝着形成水华的蓝藻演替。如太湖梅梁湾6月浮游蓝丝藻属占据主导地位时, TP浓度高于其他月份[30]。当磷浓度很高时, 比较容易出现浮游动物对绿藻的捕食, 遏制绿藻的生长, 促使蓝藻快速生长形成优势种[31]。当氮、磷营养盐浓度严重超标时, 较小的氮磷比〔ρ(TN):ρ(TP)≈10〕有利于形成拟浮丝藻水华[29]。拟浮丝藻属和浮游蓝丝藻属为相近属[20]。骆马湖9月ρ(TN):ρ(TP)约为12.5:1。因此, 9月浮游蓝丝藻生物量的增加可能与ρ(TN):ρ(TP)比值较低及TP浓度增加有关。

另外, 薛迪[32]调查结果显示南四湖夏、秋季水体中蓝藻门的伪鱼腥藻占据优势地位, 此次调查中的浮游蓝丝藻与其都是丝状蓝藻。在水体中丝状蓝藻藻丝游离漂浮, 利用自身的悬浮机制, 能够获得更多的营养和光照[32]。流速较高的水体不利于相对体积质量较小的细胞停留, 一般在流水系统中蓝藻门与硅藻门较多[33], 但微囊藻的悬浮机制与水体的稳定性密切相关[34], 在水力滞留时间短的水体浮游植物组成中会失去竞争优势[35], 而丝状体在流水和浑浊水体更占优势[36]。骆马湖是过水性湖泊, 径流量较大, 流速较快, 而且常年采砂作业导致骆马湖每遇大风之后水质浑浊, 且笔者研究中9月骆马湖水体透明度与PC浓度呈显著负相关(P < 0.05, 数据未列出)。因此, 水体流速较快和透明度较低也可能是浮游蓝丝藻占优势的原因之一。综上所述, 9月骆马湖浮游蓝丝藻占优势, 可能是由其本身的生理特点以及温度、营养盐、水体流速、透明度及其他生物等诸多环境因素所引发。

RDA分析主要针对环境因子与浮游植物群落结构变化的关系, 在自然条件下, 多个环境因子之间往往是相互作用, 共同对浮游植物群落结构起作用。蒋伊能[37]通过RDA分析发现营养盐、水温和离子类环境数据共同作用时解释了抚仙湖浮游植物群落结构变化的25%, 未解释的为75%。邻近的南四湖夏季主要受水温和总磷浓度的影响, 秋季主要受总氮浓度、总磷浓度和pH的影响[32]。笔者研究中RDA分析结果表明, 氮磷比、水温和溶解氧浓度的共同作用解释了骆马湖夏季浮游植物群落结构变化的28.16%, 未解释的为71.84%。由于骆马湖生态环境的复杂性, 其浮游植物群落结构的变化可能受到此次调查中未检测的其他水质因子、气象条件、水动力以及生物因素等的驱动。

PC是蓝藻的特征色素蛋白之一, 通常被用来表征蓝藻生物量。大量研究表明, 高磷负荷会导致较高的浮游植物生物量。但在笔者研究中, 8月蓝藻PC浓度和TP浓度变化未显示出预期的正相关关系, 可能是物理因素而不仅仅是磷决定了浮游植物密度和生物量。这种解释与ZHANG等[38]的研究结论相似, 其研究表明抚仙湖浮游植物密度的主要驱动因子是物理变量, 特别是在水下光照条件下, 尽管总磷浓度(1.89~21.8 μg·L-1)很低。XIE等[39]在武汉东湖的围隔实验证明氮磷比无论小于还是大于29都会发生蓝藻(微囊藻)水华。笔者研究中8月ρ(TN)和ρ(TP)平均值分别为3.02和0.23 mg·L-1, 已经充分满足了蓝藻生长的需要[40], 说明ρ(TN):ρ(TP)比值已经不是影响8月蓝藻生长的主要因子。研究表明, 较低的ρ(TN):ρ(TP)比值是蓝藻水华产生的结果[39], 但笔者研究中8月蓝藻生物量较低, 可能没有足够的生物量来引发ρ(TN):ρ(TP)比值降低。当氮、磷浓度严重超标时, 拟浮丝藻大规模的爆发还会使ρ(TN):ρ(TP)比值下降[29], 因而9月PC浓度与ρ(TN):ρ(TP)比值、TP浓度的关系可能受浮游蓝丝藻大量增殖的影响较大。因此, 蓝藻PC与TP浓度及ρ(TN):ρ(TP)比值的关系在8月和9月相反可能主要是由于8月蓝藻不占优势及其他物理因素所致。

4 结论

通过对骆马湖夏季浮游植物的调查, 分析骆马湖水质现状及水体营养状态, 探讨骆马湖浮游植物群落结构变化与环境因子的关系, 得出以下结论:

(1) 与往期骆马湖调查结果相比, 骆马湖夏季多样性指数差异不大, 细胞丰度明显增加, 优势种较集中于蓝藻门和绿藻门, 尤其是浮游蓝丝藻明显增多。从近几年优势属种类、藻类细胞密度变化来看, 骆马湖富营养化程度依然在加剧。

(2) RDA分析结果表明, 水温、溶解氧浓度和氮磷比的共同作用解释了骆马湖夏季浮游植物群落结构变化的28.16%。其中, 骆马湖浮游植物群落结构的变化受采砂、围网养殖及夏季人类活动影响较大。骆马湖夏季总体为浮游蓝丝藻占优势, 可能是由浮游蓝丝藻本身的生理特点以及温度、营养盐、水体流速、透明度及其他生物等诸多环境因素所引发。而比浮游蓝丝藻生长速率高的微囊藻几乎不占优势, 说明骆马湖在调查期间可能属于微囊藻生长受抑水体, 而骆马湖水体浮游蓝丝藻水华暴发几率会增大。

参考文献
[1]
申霞, 洪大林, 谈永锋, 等. 骆马湖生态环境现状及其保护措施[J]. 水资源保护, 2013, 29(3): 39-43, 50. [ SHEN Xia, HONG Da-lin, TAN Yong-feng, et al. Ecological Environment of Luoma Lake and Protection Measures[J]. Water Resources Protection, 2013, 29(3): 39-43, 50. DOI:10.3969/j.issn.1004-6933.2013.03.009] (0)
[2]
BADRUZZAMAN M, PINZON J, OPPENHEIMER J, et al. Sources of Nutrients Impacting Surface Waters in Florida:A Review[J]. Journal of Environmental Management, 2012, 109: 80-92. DOI:10.1016/j.jenvman.2012.04.040 (0)
[3]
DOMINGUES R B, BARBOSA A, GALVO H. Constraints on the Use of Phytoplankton as a Biological Quality Element Within the Water Framework Directive in Portuguese Waters[J]. Marine Pollution Bulletin, 2008, 56(8): 1389-1395. DOI:10.1016/j.marpolbul.2008.05.006 (0)
[4]
PADISÁK J, HAJNAL É, NASELLI-FLORES L, et al. Convergence and Divergence in Organization of Phytoplankton Communities Under Various Regimes of Physical and Biological Control[J]. Hydrobiologia, 2010, 639(1): 205-220. DOI:10.1007/s10750-009-0021-5 (0)
[5]
蒋伊能, 张涛, 陈丽, 等. 抚仙湖冬-春季浮游植物群落结构变化及其影响因子识别[J]. 海洋湖沼通报, 2016(6): 87-98. [ JIANG Yi-neng, ZHANG Tao, CHEN Li, et al. Ldentification of Phytoplankton Community Changes and Environmental Factors During the Winter-Spring Transition in Fuxian Lake[J]. Transactions of Oceanology and Limnology, 2016(6): 87-98.] (0)
[6]
彭凯, 李太民, 刘利, 等. 苏北骆马湖浮游植物群落结构及其水质生物评价[J]. 湖泊科学, 2018, 30(1): 183-191. [ PENG Kai, LI Tai-min, LIU Li, et al. Community Structure of Phytoplankton and Bio-Assessment of Water Quality in Lake Luoma, Northern Jiangsu Province[J]. Journal of Lake Sciences, 2018, 30(1): 183-191.] (0)
[7]
张小明, 李超, 李倩. 骆马湖浮游藻类调查与水质评价[J]. 治淮, 2016(1): 29-30. DOI:10.3969/j.issn.1001-9243.2016.01.017 (0)
[8]
EBINA J. Simultaneous Determination of Total Nitrogen and Total Phosphorus in Water Using Peroxodisulfate Oxidation[J]. Water Research, 1983, 17(12): 1721-1726. DOI:10.1016/0043-1354(83)90192-6 (0)
[9]
阎荣, 孔繁翔, 韩小波. 太湖底泥表层越冬藻类群落动态的荧光分析法初步研究[J]. 湖泊科学, 2004, 16(2): 163-168. [ YAN Rong, KONG Fan-xiang, HAN Xiao-bo. Analysis of the Recruitment of the Winter Survival Algae on the Sediments of Lake Taihu by Fluorometry[J]. Journal of Lake Science, 2004, 16(2): 163-168. DOI:10.3321/j.issn:1003-5427.2004.02.010] (0)
[10]
章宗涉, 黄祥飞. 淡水浮游生物研究方法[M]. 北京: 科学出版社, 1991: 232-246. (0)
[11]
MUYLAERT K, VAN DER GUCHT K, VLOEMANS N, et al. Relationship Between Bacterial Community Composition and Bottom-Up Versus Top-Down Variables in Four Eutrophic Shallow Lakes[J]. Applied and Environmental Microbiology, 2002, 68(10): 4740-4750. DOI:10.1128/AEM.68.10.4740-4750.2002 (0)
[12]
JIANG Y J, HE W, LIU W X, et al. The Seasonal and Spatial Variations of Phytoplankton Community and Their Correlation With Environmental Factors in a Large Eutrophic Chinese Lake (Lake Chaohu)[J]. Ecological Indicators, 2014, 40: 58-67. DOI:10.1016/j.ecolind.2014.01.006 (0)
[13]
张民, 于洋, 钱善勤, 等. 云贵高原湖泊夏季浮游植物组成及多样性[J]. 湖泊科学, 2010, 22(6): 829-836. [ ZHANG Min, YU Yang, QIAN Shan-qin, et al. Phytoplankton Community Structure and Biodiversity in Summer Yunnan-Guizhou Plateau Lakes[J]. Journal of Lake Sciences, 2010, 22(6): 829-836.] (0)
[14]
杨士建. 骆马湖富营养化发生机制与防治途径初探[J]. 环境导报, 2003(9): 11-12. (0)
[15]
邹伟, 李太民, 刘利, 等. 苏北骆马湖大型底栖动物群落结构及水质评价[J]. 湖泊科学, 2017, 29(5): 1177-1187. [ ZOU Wei, LI Tai-min, LIU Li, et al. Macrozoobenthic Community Structure and Water Quality Assessment of Lake Luoma, Jiangsu Province, China[J]. Journal of Lake Sciences, 2017, 29(5): 1177-1187.] (0)
[16]
庞科, 姚锦仙, 王昊, 等. 额尔古纳河流域秋季浮游植物群落结构特征[J]. 生态学报, 2011, 31(12): 3391-3398. [ PANG Ke, YAO Jin-xian, WANG Hao, et al. Community Structure Characteristics of Phytoplankton in Argun River Drainage Area in Autumn[J]. Acta Ecologica Sinica, 2011, 31(12): 3391-3398.] (0)
[17]
桂佳, 胡韧, 辛艳萍, 等. 一株来源于水库的浮游蓝丝藻SK3株的生长及产毒特性[J]. 湖泊科学, 2010, 22(3): 400-404. [ GUI Jia, HU Ren, XIN Yan-ping, et al. Growth and Toxin-Producing Characteristics of Planktothrix agardhii SK3 From a Reservoir[J]. Journal of Lake Sciences, 2010, 22(3): 400-404.] (0)
[18]
代玲玲.三峡水库富营养化藻类特征及环境因素影响研究[D].重庆: 重庆大学, 2007. [DAI Ling-ling.Study on Euthophication Algal Characters and Environmental Factors Effect of Three Gorge Reservoir[D].Chongqing: Chongqing University, 2007.] (0)
[19]
杨刚, 胡瑾.洪泽湖、骆马湖藻类监测及富营养化评价[R].淮河流域综合治理与开发科技论坛.:, 2010. [YANG Gang, HU Jin.Algae Monitoring and Eutrophication Assessment in Hongze Lake and Luoma Lake[R].Science and Technology Forum on Comprehensive Management and Development of Huaihe River Basin.:, 2010.] (0)
[20]
吴忠兴, 余博识, 彭欣, 等. 中国水华蓝藻的新记录属:拟浮丝藻属(Planktothricoides)[J]. 武汉植物学研究, 2008, 26(5): 461-465. [ WU Zhong-xing, YU Bo-shi, PENG Xin, et al. Planktothricoides, a Newly Recorded Genus of Water Bloom Forming Cyanophyta in China[J]. Journal of Wuhan Botanical Research, 2008, 26(5): 461-465. DOI:10.3969/j.issn.2095-0837.2008.05.007] (0)
[21]
陈宇炜, 李朋富, Martin DOKULIL. 浮游藻类三个常见属(颤藻属、直链硅藻属和针杆藻属)学名变更的解释[J]. 湖泊科学, 2003, 15(1): 85-94. [ CHEN Yu-wei, LI Peng-fu, DOKULIL M. Clarification for Modifying 3 Common Phytoplankton Genera Names[J]. Journal of Lake Science, 2003, 15(1): 85-94. DOI:10.3321/j.issn:1003-5427.2003.01.013] (0)
[22]
KILHAM S S.Phytoplankton Ecology: Structure, Function and Fluctuation[M]//HARRIS G P.Trends in Ecology & Evolution.London: Chapman & Hall, 1987: 380.DOI: 10.1016/0169-5347(87)90143-1. (0)
[23]
田家怡, 田静, 申保忠, 等. 南水北调东线工程输水沿线浮游植物分布与多样性[J]. 山东科学, 2008, 21(1): 28-34. [ TIAN Jia-yi, TIAN Jing, SHEN Bao-zhong, et al. Phytoplankton Distribution and Diversity Along the East Route of South-to-North Water Transfer Project[J]. Shandong Science, 2008, 21(1): 28-34.] (0)
[24]
冷春梅, 王亚楠, 董贯仓, 等. 南四湖浮游植物群落结构特征及季相变化[J]. 山东农业大学学报(自然科学版), 2014, 45(2): 193-197. [ LENG Chun-mei, WANG Ya-nan, DONG Guan-cang, et al. Community Structure Characteristics and the Seasonal Variation of Phytoplankton in Lake Nansi[J]. Journal of Shandong Agricultural University(Natural Science Edition), 2014, 45(2): 193-197. DOI:10.3969/j.issn.1000-2324.2014.02.07] (0)
[25]
DOU M, ZHANG Y, LI G Q. Temporal and Spatial Characteristics of the Water Pollutant Concentration in Huaihe River Basin From 2003 to 2012, China[J]. Environmental Monitoring and Assessment, 2016, 188(9): 522. DOI:10.1007/s10661-016-5503-0 (0)
[26]
孔繁翔, 高光. 大型浅水富营养化湖泊中蓝藻水华形成机理的思考[J]. 生态学报, 2005, 25(3): 589-595. [ KONG Fan-xiang, GAO Guang. Hypothesis on Cyanobacteria Bloom-Forming Mechanism in Large Shallow Eutrophic Lakes[J]. Acta Ecologica Sinica, 2005, 25(3): 589-595. DOI:10.3321/j.issn:1000-0933.2005.03.028] (0)
[27]
庞燕, 金相灿, 储昭升, 等. 几种捕食性微小动物对铜绿微囊藻及孟氏浮游蓝丝藻生长的抑制作用[J]. 环境科学学报, 2008, 28(1): 56-61. [ PANG Yan, JIN Xiang-can, CHU Zhao-sheng, et al. Degradation of Growing Microcystis aeruginosa and Planktothrix mougeotii by Several Predacious Microanimals[J]. Acta Scientiae Circumstantiae, 2008, 28(1): 56-61. DOI:10.3321/j.issn:0253-2468.2008.01.009] (0)
[28]
金相灿, 储昭升, 杨波, 等. 温度对水华微囊藻及孟氏浮游蓝丝藻生长、光合作用及浮力变化的影响[J]. 环境科学学报, 2008, 28(1): 50-55. [ JIN Xiang-can, CHU Zhao-sheng, YANG Bo, et al. Effects of Temperature on Growth, Photosynthesis and Buoyancy Regulation of the Cyanobacteria Microcystis flos-aquae and Planktothrix mougeotii[J]. Acta Scientiae Circumstantiae, 2008, 28(1): 50-55.] (0)
[29]
税永红, 唐亚, 李巧巧, 等. 成都犀池拟浮丝藻水华成因探析[J]. 四川环境, 2010, 29(1): 17-23. [ SHUI Yong-hong, TANG Ya, LI Qiao-qiao, et al. Formation Causes of Water Bloom of Planktothricoides in the Xichi Lake of Chengdu[J]. Sichuan Environment, 2010, 29(1): 17-23. DOI:10.3969/j.issn.1001-3644.2010.01.004] (0)
[30]
宋晓兰, 刘正文, 潘宏凯, 等. 太湖梅梁湾与五里湖浮游植物群落的比较[J]. 湖泊科学, 2007, 19(6): 643-651. [ SONG Xiao-lan, LIU Zheng-wen, PAN Hong-kai, et al. Phytoplankton Community Structure in Meiliang Bay and Lake Wuli of Lake Taihu[J]. Journal of Lake Sciences, 2007, 19(6): 643-651. DOI:10.3321/j.issn:1003-5427.2007.06.004] (0)
[31]
王小冬.蓝藻水华形成的实验模拟[D].南京: 中国科学院南京地理与湖泊研究所, 2010. [WANG Xiao-dong.Experimental Simulation of Cyanobacterial Bloom Formation[D].Nanjing: Nanjing Institute of Geography and Lakes, Chinese Academy of Sciences, 2010.] (0)
[32]
薛迪.南四湖浮游植物群落结构时空变化特征[D].济南: 山东大学, 2016. [XUE Di.Spatial and Temporal Variations of Phytoplankton Community Structure in Nansi Lake[D].Jinan: Shandong University, 2016.] (0)
[33]
WEHR J D, DESCY J P. Use of Phytoplankton in Large River Management[J]. Journal of Phycology, 1998, 34(5): 741-749. DOI:10.1046/j.1529-8817.1998.340741.x (0)
[34]
OLIVER R L. Floating and Sinking in Gas-Vacuolate Cyanobacteria1[J]. Journal of Phycology, 1994, 30(2): 161-173. DOI:10.1111/j.0022-3646.1994.00161.x (0)
[35]
REYNOLDS C S. The Ecology of Phytoplankton[M]. Cambridge: Cambridge University Press, 2006: 1-435. (0)
[36]
REYNOLDS C S. The Long, the Short and the Stalled:On the Attributes of Phytoplankton Selected by Physical Mixing in Lakes and Rivers[J]. Hydrobiologia, 1994, 289(1/2/3): 9-21. (0)
[37]
蒋伊能.抚仙湖浮游植物群落与多样性的时空变化特征及其影响因子研究[D].昆明: 云南师范大学, 2017. [JIANG Yi-neng.Spatio-Temporal Distribution and Environmental Drivers of Phytoplankton Community and Diversity in Fuxian Lake[D].Kunming: Yunnan Normal University, 2017.] (0)
[38]
ZHANG X, XIE P, CHEN F Z, et al. Driving Forces Shaping Phytoplankton Assemblages in Two Subtropical Plateau Lakes With Contrasting Trophic Status[J]. Freshwater Biology, 2007, 52(8): 1463-1475. DOI:10.1111/j.1365-2427.2007.01776.x (0)
[39]
XIE L Q, XIE P, LI S X, et al. The Low TN:TP Ratio, a Cause or a Result of Microcystis Blooms?[J]. Water Research, 2003, 37(9): 2073-2080. DOI:10.1016/S0043-1354(02)00532-8 (0)
[40]
许海, 朱广伟, 秦伯强, 等. 氮磷比对水华蓝藻优势形成的影响[J]. 中国环境科学, 2011, 31(10): 1676-1683. [ XU Hai, ZHU Guang-wei, QIN Bo-qiang, et al. Influence of Nitrogen-Phosphorus Ratio on Dominance of Bloom-Forming Cyanobacteria (Microcystis aeruginosa)[J]. China Environmental Science, 2011, 31(10): 1676-1683.] (0)