环境科学  2023, Vol. 44 Issue (4): 2052-2061   PDF    
典型湖泊有机聚集体时空特征及驱动因子
谢贵娟1,2,3, 龚伊2, 朱富成1, 刘昌利1, 卢宝伟1, 邓辉1, 汤祥明2,3     
1. 皖西学院生物与制药工程学院, 六安 237012;
2. 中国科学院南京地理与湖泊研究所湖泊与环境国家重点实验室, 南京 210008;
3. 中国科学院大学资源与环境学院, 北京 100049
摘要: 有机聚集体(OA)是水体中物质循环和能量流通的重要枢纽,但在不同营养水平湖泊中的比较研究仍相当匮乏.于2019~2021年间,采集了不同季节贫营养抚仙湖、中营养天目湖、中度富营养太湖和重度富营养星云湖表层水样,利用扫描电镜、表面荧光显微镜和流式细胞仪等研究了4个典型湖泊中OA及其附着细菌的时空变化特征.结果表明:①抚仙湖、天目湖、太湖和星云湖中,OA丰度的年均值分别为1.4×104、7.0×104、27.7×104和16.0×104ind·mL-1;OA附着细菌丰度的年均值分别为0.3×106、1.9×106、4.9×106和6.2×106 cells·mL-1;附着细菌占总细菌的比值分别为30%、31%、50%和38%.②夏季OA丰度显著高于秋冬季;但夏季附着细菌占总细菌的比值为26%,显著低于其它3个季节.③湖泊营养状态是影响OA及其附着细菌丰度的最重要环境因子,可单独解释OA及其附着细菌丰度时空变化的50%和68%.④OA上富含营养盐和有机质,尤其是星云湖,其OA中颗粒态磷和氮占比高达69%和59%,有机质占比达79%.未来气候变化和湖泊藻类水华扩张背景下,藻源性OA在有机质矿化降解及营养盐循环中的作用可能会进一步加大.
关键词: 有机聚集体(OA)      附着细菌      湖泊营养状态      藻类水华      沉积物再悬浮      营养盐循环      驱动因子     
Spatio-temporal Characteristics of Organic Aggregates and the Driving Factors in Typical Lakes
XIE Gui-juan1,2,3 , GONG Yi2 , ZHU Fu-cheng1 , LIU Chang-li1 , LU Bao-wei1 , DENG Hui1 , TANG Xiang-ming2,3     
1. College of Biology and Pharmaceutical Engineering, West Anhui University, Lu'an 237012, China;
2. State Key Laboratory of Lake Science and Environment, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing 210008, China;
3. College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: Organic aggregates (OA) are the important circulation hub of matter and energy in aquatic ecosystems. However, the comparison studies on OA in lakes with different nutrient levels are limited. In this study, spatio-temporal abundances of OA and OA-attached bacteria (OAB) in oligotrophic Lake Fuxian, mesotrophic Lake Tianmu, middle-eutrophic Lake Taihu, and hyper-eutrophic Lake Xingyun were investigated in different seasons during 2019-2021 using a scanning electron microscope, epi-fluorescence microscope, and flow cytometry. The results showed that: ① the annual average abundances of OA in Lake Fuxian, Lake Tianmu, Lake Taihu, and Lake Xingyun were 1.4×104, 7.0×104, 27.7×104, and 16.0×104 ind·mL-1, whereas the annual average abundances of OAB in the four lakes were 0.3×106, 1.9×106, 4.9×106, and 6.2×106 cells·mL-1. The ratios of OAB: total bacteria (TB) in the four lakes were 30%, 31%, 50%, and 38%, respectively. ② OA abundance in summer was significantly higher than that in autumn and winter; however, the ratio of OAB: TB in summer was approximately 26%, which was significantly lower than that in the other three seasons. ③ Lake nutrient status was the most important environmental factor that affected the abundance variations of OA and OAB, accounting for 50% and 68% of the spatio-temporal variations in OA and OAB abundances. ④ Nutrient and organic matters were enriched in OA, especially in Lake Xingyun; the proportions of particle phosphorous, particle nitrogen, and organic matters in this lake were as high as 69%, 59%, and 79%, respectively. Under the circumstance of future climate change and the expansion of lake algal blooms, the effects of algal-originated OA in the degradation of organic matters and nutrient recycling would be increased.
Key words: organic aggregates(OA)      particle-attached bacteria      lake nutrient status      algal blooms      sediment resuspension      nutrient recycling      driving factor     

有机聚集体(organic aggregates, OA)是水生态系统中的重要组分之一[1, 2].它主要由浮游生物残体及其代谢产物、各种有机碎屑、无机颗粒、自养和异养微生物等, 通过物理、化学或生物过程而粘附在一起的不规则多孔颗粒物[3].OA及其周际会形成一个营养丰富的微环境, 吸引大量细菌定殖其中.水体中细菌可分为浮游细菌和颗粒附着细菌, 浮游细菌营自由生活, 而附着细菌一般在颗粒物表面营附着生活[4].有研究表明, OA附着细菌的数量比同体积水体中浮游细菌的数量高2~3个数量级[2, 5], 其代谢活性和生长速率也比浮游细菌高数倍[6~8].OA形成以后, 仍和周围水体的微生物、溶解性物质和其他颗粒物进行相互作用[9, 10].在大量高活性附着微生物的作用下, OA上的颗粒有机物(POM)能够被快速分解并释放溶解性有机物(DOM)至水体中, 供应浮游细菌生长[11, 12].OA也会被部分小型食碎屑动物或原生动物牧食, 从而将有机质向更高一级的食物链传递.由于其在有机质降解和食物网中的重要作用, OA被视为水体中营养盐矿化、能量流通、元素生物地球化学循环的热点[13~18], 和微生物相互作用的枢纽[19, 20].

作为自然水体中最具生物活性的微型“岛屿”[20], OA形状类似不规则多孔性的雪花团, 有“海洋雪花”[15, 21]、“湖泊雪花”[2]和“河流雪花”[22]之称.目前, 有关OA的研究主要集中在海洋和深水湖泊, 因为在这些环境中, 由于沉降距离长, 容易形成大型OA(粒径>500 μm), 便于收集观测.而水动力扰动强烈的浅水湖泊中主要以小型OA(粒径范围为5~500 μm)为主, 相关研究主要集中在太湖[23, 24].由于浮游藻类是OA的主要来源之一[25], 影响藻类丰度的因子, 如水体营养盐水平和季节变化等, 均直接或间接影响水体中OA的来源、组成、理化特征及其附着微生物的群落演替[26, 27].当前在气候变化背景下, 水体富营养化对水生态系统的影响加剧, 全球尺度上湖泊藻类水华呈扩张趋势[28, 29], 深刻影响着湖泊生态系统的结构和功能.随着营养水平的增高, 这种源自藻类的有机聚集体“微岛”的丰度, 预期也将增高, 其生态学作用和环境效应不容忽视.然而, 在湖泊生态系统中, 目前仍缺乏不同营养水平和不同水深湖泊中OA及其附着微生物时空差异的对比研究.鉴于此, 本文选择4个典型湖泊为研究对象, 即贫营养深水湖(云南抚仙湖)、中营养人工水库型湖泊(江苏天目湖)、富营养大型浅水湖泊(江苏太湖)和重度富营养中等水深湖泊(云南星云湖), 通过对典型湖泊中OA及其附着细菌4个季节的调查分析, 探讨营养盐水平和季节变化对OA及其附着细菌时空特征的影响, 旨在为未来气候变化背景下, 以有机聚集体为核心载体的湖泊营养盐循环及其生态环境效应研究提供基础数据和新的视角.

1 材料与方法 1.1 研究区概况

太湖是我国第三大淡水湖泊, 水面面积约2 338 km2, 蓄水量约44.3亿m3, 平均水深1.9 m, 最大水深3 m.太湖水动力扰动强烈, 入湖河流多, 其中位于湖西的大浦河为入湖水量最大的河流[图 1(a)].由于地处经济发达的江苏平原河网地区, 外源污染负荷长期通过西北部的河流输入, 导致太湖富营养化和蓝藻水华问题十分突出[30~32].

图 1 4个不同营养水平湖泊的采样位点示意 Fig. 1 Sampling sites of the four typical lakes with different trophic levels

天目湖(沙河水库)水面面积12 km2[图 1(b)], 蓄水量约1.1亿m3, 平均水深7 m, 最大水深14 m[33].天目湖自1961年建成蓄水以来, 一直是周边溧阳市的饮用水水源地, 水体处于中营养水平[34, 35].

抚仙湖是我国最大的高原深水湖, 水面面积约216.6 km2[图 1(c)], 在水位为1 721 m时, 最大水深155 m, 平均水深89.6 m[36].它是我国蓄水量最大的淡水湖, 常年蓄水量为189亿m3.2011~2020年间, 抚仙湖水质整体为Ⅰ类, 处于贫营养状态[37].

星云湖是云南省九大高原湖泊之一, 和抚仙湖同属珠江流域南盘江水系, 水面面积约34.7 km2[图 1(c)], 蓄水量约1.84亿m3, 平均水深5.3 m, 最大水深11 m[36].有研究显示星云湖目前水污染严重, 水质为劣Ⅴ类, 是重度富营养湖泊, 常年有蓝藻水华暴发[38, 39].为减小对抚仙湖水质的影响, 当地政府已于2002年在连接星云湖和抚仙湖的隔河建立闸控, 因此目前这两个湖已不连通.

1.2 样点布设和采样

在上述4个湖泊中各设置一个断面, 进行样点布设.在太湖, 从最大入湖河流大浦河的入口处, 至湖东部漫山岛一线设置了12个采样点, 在其它3个湖泊各设置10个采样点(图 1).分别于2019年11月(秋季)、2020年4月(春季)、2020年7月(夏季)和2021年1月(冬季)进行采样.在每个采样点表层约50 cm处, 用采水器采集5~20 L水样, 置预先清洗过的塑料桶中, 运至实验室, 进行水体理化参数的分析.现场取47.5 mL原水, 置于灭过菌的50 mL离心管中(管中事先加入40%的无颗粒甲醛2.5 mL, 使甲醛终浓度为2%), 放入含有冰袋的保温箱中, 带回实验室4℃保存, 用于OA和细菌的检测; 藻类计数的原水(40 mL)带回实验室后保存在-20℃.后续分析在一周内完成.

1.3 环境参数分析

在每个取样点水下约50 cm处, 用多参数水质分析仪(YSI EXO2, Yellow Springs Instruments Inc., USA)测量水温(WT, ℃); 用Secchi盘测量水体透明度(SD).总氮(TN)、总溶解氮(TDN)、总磷(TP)、总溶解磷(TDP)、叶绿素a(Chla)、化学需氧量(COD)、溶解性有机碳(DOC)、总悬浮固体(SS)和烧失量(LOI)的测量, 参照标准方法进行[40].用TN和TDN的差值计算颗粒态氮(PN); 用TP和TDP的差值计算颗粒态磷(PP).用以下公式计算SS中有机物(OM)的质量分数(COM):

根据文献[41], 采用SD、TN、TP、Chla和COD这5个参数计算综合营养状态指数(TLI), 以表征各湖泊营养状态的时空变化.TLI≤30, 为贫营养; 30<TLI≤ 50, 为中营养; TLI>50, 为富营养.其中, 50<TLI≤ 60, 为轻度富营养; 60<TLI≤ 70, 为中度富营养; TLI>70, 为重度富营养.

1.4 有机聚集体、细菌和藻类的检测

有机聚集体丰度的检测和扫描电镜形态学观测参照文献[23, 42].除抚仙湖水样用原水外, 其他湖水样进行适当稀释后, 取出5 mL加100 μL的DAPI(4′, 6-二脒基-2-苯基吲哚, 100 μg·mL-1), 染色10 min; 用手持真空泵在较低压力下抽滤至0.2 μm黑色滤膜(PoreticsTM)上, 加无荧光浸镜油(Immersol® 518F, ZEISS)制片后, 用表面荧光显微镜(Axio Imager.A2, Zeiss)观察拍照, 每个样品至少随机拍摄20张照片, 包含至少100个OA.用Image J 1.52软件计算各个OA的表面积, 本研究中将表面积大于5 μm2, 且有明显细菌附着的颗粒归为OA.OA丰度的计算公式如下:

式中, AOA为有机聚集体丰度, ind·mL-1; WA为滤膜的有效过滤面积, μm2; FA为在所用的放大倍数下照片的实际面积, μm2; N为在所用的放大倍数下照片中的有机聚集体个数; F为水样的稀释系数; V为所过滤的水样的体积, mL.

用流式细胞仪(Becton Dickinson)进行藻类丰度(Aalgae)和细菌丰度的计数, 具体步骤见文献[43].总细菌(total bacteria, TB)丰度(ATB)和浮游细菌(free-living bacteria, FLB)丰度(AFLB)分别用原水和5 μm滤膜过滤后的滤液进行检测[44].OA附着细菌(OA-attached bacteria, OAB)丰度:(AOAB)=ATB-AFLB.

1.5 统计分析和绘图

本研究中除采样点位图用ArcMap 10.5绘制外, 其他统计分析和可视化在R 3.5.3(https://www.r-project.org)平台和RStudio 1.1.463交互界面完成.由于所测得的数据部分不符合正态分布, 因此, 利用“boxplerk.R”脚本[45], 通过非参数Kruskal-Wallis检验来比较不同样品间差异的显著性; 参照文献[46, 47], 利用mgcv 1.8-26包[48], 通过广义可加模型(GAM)分析影响AOAAOAB的最优环境因子组合, 利用hier.part 1.0-6包采用层次分割理论来分解每个环境因子的相对贡献率[49].

2 结果与分析 2.1 不同湖泊主要环境参数比较

抚仙湖、天目湖、太湖和星云湖TLI的均值分别为24、46、62和73, 分别符合贫营养、中营养、中度富营养和重度富营养的标准, 非参数统计检验表明各个湖泊间TLI差异显著(图 2).除WT外, 所测得其它参数在4个湖泊间差异均显著(P<0.001).星云湖pH最高, 均值达9.0, 显著高于抚仙湖(8.8)、太湖(8.3)和天目湖(8.2).水体中TN、TP、TDP、COD、DOC和Chla在4个湖泊间变化趋势相似, 从贫营养的抚仙湖至重度富营养的星云湖, 其均值依次增高.TDN、Aalgae和SS从抚仙湖到太湖均值依次增高, 但星云湖低于太湖且高于其它两个湖泊, 随着营养程度的提高总体呈左偏单峰曲线.

箱线图内标识Kruskal-Wallis秩和检验总体结果:***表示P<0.001, ns表示差异不显著, P>0.05; 小写字母表示组间差异的显著性, 相同字母表示差异不显著(P>0.05), 不同字母表示差异显著(P<0.05); 箱体内圆表示平均值, 横线表示中位数, 箱体外圆表示异常值 图 2 4个湖泊主要环境参数比较 Fig. 2 Comparisons of main environmental parameters among the four lakes

颗粒物上, 贫营养抚仙湖的COM最高, 均值达86%; 太湖最低, 均值为23%; 星云湖COM也较高, 均值为79%(图 2).PN占比在不同湖泊间差异显著.星云湖最高, 均值为59%; 其次是太湖, 均值为30%; 抚仙湖和天目湖的均值较低, 分别为22%和20%.星云湖、太湖和天目湖中PP占比依次为69%、64%和64%, 三者差异不显著, 但均显著高于抚仙湖的40%.

2.2 典型有机聚集体形态特征

表面荧光显微镜和扫描电镜观测结果显示, 有机聚集体形态组成主要分为两大类, 一类是絮状的源自有机碎屑和透明胞外物质聚合而成的聚集体; 另一类主要由衰亡藻类及其表面附着细菌构成的藻源性有机聚集体(图 3).需要说明的是, 处于旺盛生长阶段的健康藻类颗粒, 其表面很少能观测到细菌的附着.在本研究中, 这种表面无细菌附着的藻类颗粒不被计数为OA.

(a)和(b)为絮状碎屑有机聚集体的表面荧光显微镜和扫描电镜图, (c)和(d)为藻源性有机聚集体的表面荧光显微镜和扫描电镜图, (a)和(c)中白色杆状或球状亮点均为附着在有机聚集体上的大量细菌 图 3 典型絮状碎屑有机聚集体和藻源性有机聚集体表面荧光显微镜和扫描电镜图 Fig. 3 Epifluorescence and scanning electron micrographs showing the typical flocculent debris and algal-originated organic aggregates

2.3 不同湖泊有机聚集体基本特征的差异

不同营养水平的湖泊之间, AOA差异极显著[图 4(a1)].抚仙湖中AOA最低, 均值为1.4×104ind·mL-1; 天目湖AOA的均值为7.0×104ind·mL-1; 太湖AOA最高, 均值为27.7×104ind·mL-1, 其变化范围也最大, 为(4.5~73.8)×104ind·mL-1; 星云湖AOA的均值为16.0×104ind·mL-1, 略低于太湖.星云湖OA表面积均值为178 μm2, 显著高于其它3个湖泊[图 4(a2)].抚仙湖、天目湖和太湖OA的平均表面积分别为74、67和90 μm2, 三者差异不显著(P>0.05).

箱线图内标识Kruskal-Wallis秩和检验总体结果:***表示P<0.001, **表示P<0.01; 小写字母表示组间差异的显著性, 相同字母表示差异不显著(P>0.05), 不同字母表示差异显著(P<0.05); 箱体内圆表示平均值, 横线表示中位数, 箱体外圆表示异常值; (a1)~(a4)为不同湖泊间差异, (b1)~(b4)不同季节间差异 图 4 有机聚集体丰度、表面积、附着细菌丰度和占比在不同湖泊和不同季节间差异的箱线图 Fig. 4 Boxplots showing the differences in AOA, surface area, AOAB, and the ratio of AOAB: ATBamong the four lakes and four seasons

星云湖AOAB最高, 均值为6.2×106 cells·mL-1; 太湖次之, 均值为4.9×106 cells·mL-1[图 4(a3)].两个湖泊之间附着细菌丰度差异不显著, 但均显著高于天目湖和抚仙湖.附着细菌占总细菌的比例以太湖最高, 均值为50%, 显著高于其它3个湖泊[图 4(a4)].星云湖、天目湖和抚仙湖附着细菌占比分别为38%、31%和30%.

2.4 不同季节有机聚集体基本特征的差异

将4个湖泊数据汇总分析, OA丰度、表面积、OA附着细菌丰度和占比在不同季节间差异显著.夏季和春季AOA最高, 均值分别为17.8×104和16.7×104ind·mL-1, 两者差异不显著; 夏季显著高于秋季(11.0×104ind·mL-1)和冬季(9.4×104ind·mL-1)[图 4(b1)].秋季OA平均表面积最大, 为138 μm2; 冬季次之, 均值为115μm2; 夏季均值为92 μm2; 3个季节间差异不显著; 秋冬季显著高于春季的均值61 μm2[图 4(b2)].

AOAB秋季均值为4.4×106cells·mL-1, 显著高于夏季(2.3×106cells·mL-1)和冬季(1.7×106cells·mL-1), 但与春季(5.3×106cells·mL-1)差异不显著[图 4(b3)].附着细菌占总细菌的比例夏季最低, 均值为26%, 显著低于冬季(44%)、秋季(43%)和春季(38%)[图 4(b4)].

2.5 有机聚集体及其附着细菌丰度与环境因子的关系

利用GAM模型分析了影响OA和OAB丰度的最优环境因子组合, 以及各环境因子对OA和OAB丰度变化的相对解释率.通过剔除共线性参数和不显著的参数, OA丰度最优环境因子GAM模型为:

附着细菌丰度最优环境因子GAM模型为:

最优模型中各个变量的平滑曲线, 显示环境因子与OA和OAB丰度的关系基本上都是非线性的[图 5(a)图 6(a)].GAM结果显示, TLI、SS、DOC和WT这4个因子合计解释了lgAOA变化的94.8%, 其中前3个变量的相对贡献分别为53%、28%和16%[图 5(b)]; TLI、DOC和WT这3个因子合计解释了lgAOAB变化的92.9%, 其中前两个因子的相对贡献分别为73%和23%[图 6(b)].4个湖TLI、SS、DOC和WT这4个因子的季节变化显著, 特别是WT, 季节差异大(图 7).

(a)AOA与环境因子的GAM结果, (b)各因子对模型的相对贡献; 实线表示响应变量方程的拟合曲线, 虚线表示95%的置信区间; 横坐标上的短黑线为边际地毯, 表示实测样本值的分布情况, 下同 图 5 有机聚集体丰度和环境因子的GAM平滑函数拟合结果 Fig. 5 Results of GAM analysis between AOA and the environmental parameters

图 6 有机聚集体附着细菌丰度和环境因子的GAM平滑函数拟合结果 Fig. 6 Results of GAM analysis between AOAB and the environmental parameters

箱线图内标识Kruskal-Wallis秩和检验结果:***表示P<0.001, **表示P<0.01, *表示P<0.05; 小写字母表示组间差异的显著性, 相同字母表示差异不显著(P>0.05), 不同字母表示差异显著(P<0.05); 箱体内圆表示平均值, 横线表示中位数, 箱体外圆表示异常值 图 7 4个湖泊TLI、SS、DOC和WT不同季节间差异的箱线图 Fig. 7 Boxplots showing the differences in TLI, SS, DOC, and WT among the four seasons in the studied lakes

3 讨论 3.1 影响有机聚集体丰度的主要环境因子

本研究发现:影响OA丰度变化的主要环境因子之一是湖泊的营养水平, 可单独解释AOA时空变化的50%[94.8%×53%; 图 5(b)].从贫营养抚仙湖至富营养太湖, 随着营养水平的升高, OA丰度有明显增高趋势.这和晁建颖等[24]在太湖不同营养水平湖区中OA丰度的研究结果类似.随着营养水平的增加, OA表面积也有增长的趋势, 且重富营养星云湖OA平均表面积显著高于其它3个湖[图 4(a2)].造成星云湖的OA平均面积较大的原因, 可能与其湖水中存在更多的藻源性OA有关.这和星云湖水体中Chla含量显著高于其它湖泊是一致的(图 2).有模拟实验显示, 营养盐显著促进了OA数量和面积的增长[26].

虽然重富营养星云湖OA平均表面积最高, 但其OA丰度却低于中富营养太湖, 说明除营养盐外, 还有其它环境因子影响OA丰度的变化.GAM结果显示SS也是影响OA丰度变化的主要因子之一, 对OA丰度变化的相对贡献为28%[图 5(b)].SS主要来源于水体中的藻类和表层沉积物再悬浮颗粒.太湖OA丰度显著高于星云湖, 可能主要和表层沉积物再悬浮有关[42, 50].太湖是大型浅水湖泊, 平均水深只有1.9 m, 水体扰动频繁, 沉积物发生再悬浮的概率和频度都要远大于平均水深为5.3 m的星云湖, 致使大量小型OA悬浮在水体中.这也符合本研究的观测结果, 即太湖SS浓度显著高于其它3个湖泊(图 2).

DOC和WT也显著影响OA丰度变化.OA上的颗粒有机物在微生物降解过程中向水体释放大量DOC, 吸引有趋化能力的细菌利用DOC, 进行碳代谢和能量流通[51].DOC对OA丰度变化的影响可能是OA上颗粒有机碳(POC)代谢的结果, 而非原因[52].在藻源性OA丰富的星云湖中, 其水体中DOC的含量也显著高于其它3个湖泊(图 2).水温的作用可能更多地影响OA丰度的季节变化.夏季水温高, 藻类生长旺盛, 藻源性OA的丰度高[53].这可能是夏季OA丰度显著高于秋冬两季的原因之一[图 4(b1)].

3.2 影响有机聚集体附着细菌丰度的主要环境因子

本研究中, 湖泊营养水平也是影响OA附着细菌丰度的主导因子, 可单独解释AOAB时空变化的68% [92.9%×73%; 图 6(b)].从贫营养抚仙湖至重富营养星云湖, 附着细菌平均丰度从0.3×106 cells·mL-1增加到6.2×106 cells·mL-1, 增加了约19倍[图 4(a3)].在营养盐缺乏的海洋环境中, OA附着细菌的丰度在(0.1~2.9)×106 cells·mL-1之间[54, 55]; 贫营养深水湖中, 其丰度范围在(0.1~0.5)×106 cells·mL-1之间[56]; 而在生产力较高营养丰富的河口生态系统中, 其丰度可高达250×106 cells·mL-1[22].

DOC和WT也是影响AOAB时空变化的主要环境因子.细菌的生长代谢除了要利用氮磷营养外, 还需要利用碳源作为合成菌体所必需的原料, 同时也作为细菌代谢的主要能量来源.DOC主要来自POM的降解释放, 可以间接反映颗粒物上碳源的多少.因此, 水体中DOC和氮磷含量高时, 总细菌和OAB丰度也相应增加.水温可能在OAB丰度的季节变化中起作用.冬季水温低时, OAB的平均丰度也最低[图 4(b3)].需要指出的是, 当夏季水温最高时, OAB的丰度及其占比均较低, 这可能和湖泊中藻类的生理特性有关.夏季藻类生长旺盛, 活性高, 竞争排斥作用导致健康的藻类颗粒上附着细菌数量少[图 3(d)].而秋季, 随着水温的降低, 衰亡藻类的相对比例升高, 附着细菌可能受到化学信号诱导[57], 从而吸附至藻类表面并大量增殖, 以分解有机质[18].因此, 秋季OAB的丰度显著高于冬季和夏季[图 4(b3)].

3.3 营养水平对有机聚集体潜在生态环境效应的影响

本研究表明, 颗粒态营养盐和有机质是湖泊中重要组分(图 2).PP∶ TP随着湖泊营养水平的提高呈上升趋势, 在重富营养星云湖中比值高达69%. PN∶ TN从中营养天目湖至重富营养星云湖也显著升高, 在星云湖中比值高达59%. COM的值在深水贫营养抚仙湖中最高, 但在重富营养星云湖中也高达79%.OA是水体中颗粒态营养盐和有机质的主要贮存库.在Maggiore湖中, 作为小型OA的透明胞外聚合颗粒物, 其碳的相对含量可占总有机碳的54%[58].水体中营养水平的提高可导致藻类的大量生长, 严重时可形成大规模的藻类水华.而内源性藻类残体作为水体中OA的重要组成部分, 因含有大量易被微生物利用的有机质[59], 可吸附大量附着细菌并在短时间内大量增殖, 使营养盐得以快速矿化[24, 60], 以正反馈形式维持藻类水华的再次暴发, 从而形成营养盐在“藻源性OA-矿化释放-藻类再快速利用-形成水华-衰亡形成藻源性OA”之间的快速流通通路.附着细菌是OA上最重要的活性组分之一, 本研究观测到随着营养水平的提高, OAB的丰度也呈显著上升的趋势, 这与上述营养盐在OA和藻类水华间快速流通的假设是可以相互印证的.需要说明的是, 湖泊营养水平的提高促进OA上有机质矿化的同时, 也会释放大量DOC和氮磷至水体中, 支撑水体中浮游细菌的增殖, 从而提高湖泊整体营养盐内源释放的速率, 使得富营养化湖泊的治理更加艰难.此外, 有研究指出, 相比水体中浮游态, 有机颗粒物上栖息的病原菌种类更丰富, 耐药基因更多, 对人类健康的负面影响也更大[20, 61].在全球气候变化导致湖泊藻类水华扩张的背景下[28, 29], 富营养化湖泊中藻源性OA在有机质矿化中的核心作用及其生态环境效应需要进一步研究.

4 结论

(1) 湖泊OA主要来源于藻类残体和沉积物再悬浮.OA上富含营养盐和有机质, 其中重富营养星云湖中颗粒态磷和氮占比高达69%和59%, 有机质占比达79%.

(2) 贫营养抚仙湖、中营养天目湖、中度富营养太湖和重度富营养星云湖中OA丰度的年均值分别为1.4×104、7.0×104、27.7×104和16.0×104ind·mL-1.星云湖OA表面积均值为178 μm2, 显著高于其它3个湖泊.季节上, 夏季OA丰度显著高于秋冬季; 但秋季OA的平均表面积最大.湖泊营养状态、悬浮物浓度、DOC浓度和水温是影响OA丰度时空变化的主要环境因子.

(3) 有机聚集体上附着细菌的丰度随着湖泊营养水平的升高呈显著上升趋势.从贫营养至重度富营养, 四个湖泊中附着细菌丰度的年均值分别为0.3×106、1.9×106、4.9×106和6.2×106 cells·mL-1.大型浅水富营养化太湖中, 附着细菌占总细菌的比值最高, 均值达50%.季节上, 春季和秋季OA上附着细菌丰度显著大于冬季; 夏季附着细菌占总细菌的比值为26%, 显著低于其它3个季节.湖泊营养状态、DOC浓度和水温是影响附着细菌丰度时空变化的主导因子.

致谢: 感谢田伟、瞿猛、李冬、薛静琛和张雨晴等在采样和理化数据测试分析中的帮助.

参考文献
[1] Alldredge A L, Cole J J, Caron D A. Production of heterotrophic bacteria inhabiting macroscopic organic aggregates (marine snow) from surface waters[J]. Limnology and Oceanography, 1986, 31(1): 68-78. DOI:10.4319/lo.1986.31.1.0068
[2] Grossart H P, Simon M. Limnetic macroscopic organic aggregates (lake snow): occurrence, characteristics, and microbial dynamics in Lake Constance[J]. Limnology and Oceanography, 1993, 38(3): 532-546. DOI:10.4319/lo.1993.38.3.0532
[3] Simon M, Grossart HP, Schweitzer B, et al. Microbial ecology of organic aggregates in aquatic ecosystems[J]. Aquatic Microbial Ecology, 2002, 28(2): 175-211.
[4] Grossart H P. Ecological consequences of bacterioplankton lifestyles: changes in concepts are needed[J]. Environmental Microbiology Reports, 2010, 2(6): 706-714. DOI:10.1111/j.1758-2229.2010.00179.x
[5] Grossart H P, Hietanen S, Ploug H. Microbial dynamics on diatom aggregates in Øresund, Denmark[J]. Marine Ecology Progress Series, 2003, 249: 69-78. DOI:10.3354/meps249069
[6] Grossart HP, Ploug H. Bacterial production and growth efficiencies: direct measurements on riverine aggregates[J]. Limnology and Oceanography, 2000, 45(2): 436-445. DOI:10.4319/lo.2000.45.2.0436
[7] Smith D C, Simon M, Alldredge A L, et al. Intense hydrolytic enzyme-activity on marine aggregates and implications for rapid particle dissolution[J]. Nature, 1992, 359(6391): 139-142. DOI:10.1038/359139a0
[8] Worm J, Søndergaard M. Dynamics of heterotrophic bacteria attached to Microcystis spp. (Cyanobacteria)[J]. Aquatic Microbial Ecology, 1998, 14(1): 19-28.
[9] Baumas C M J, Le Moigne F A C, Garel M, et al. Mesopelagic microbial carbon production correlates with diversity across different marine particle fractions[J]. The ISME Journal, 2021, 15(6): 1695-1708. DOI:10.1038/s41396-020-00880-z
[10] Enke T N, Leventhal G E, Metzger M, et al. Microscale ecology regulates particulate organic matter turnover in model marine microbial communities[J]. Nature Communications, 2018, 9(1). DOI:10.1038/s41467-018-05159-8
[11] Cho B C, Azam F. Major role of bacteria in biogeochemical fluxes in the ocean's interior[J]. Nature, 1988, 332(6163): 441-443. DOI:10.1038/332441a0
[12] Zhang Y, Xiao W, Jiao N Z. Linking biochemical properties of particles to particle-attached and free-living bacterial community structure along the particle density gradient from freshwater to open ocean[J]. Journal of Geophysical Research: Biogeosciences, 2016, 121(8): 2261-2274. DOI:10.1002/2016JG003390
[13] Azam F, Long R A. Sea snow microcosms[J]. Nature, 2001, 414(6863): 495-498. DOI:10.1038/35107174
[14] Azam F, Smith D C, Steward G F, et al. Bacteria-organic matter coupling and its significance for oceanic carbon cycling[J]. Microbial Ecology, 1994, 28(2): 167-179. DOI:10.1007/BF00166806
[15] Alldredge A L, Silver M W. Characteristics, dynamics and significance of marine snow[J]. Progress in Oceanography, 1988, 20(1): 41-82. DOI:10.1016/0079-6611(88)90053-5
[16] Grossart H P, Simon M. Bacterial colonization and microbial decomposition of limnetic organic aggregates (lake snow)[J]. Aquatic Microbial Ecology, 1998, 15(2): 127-140.
[17] Stocker R, Seymour J R, Samadani A, et al. Rapid chemotactic response enables marine bacteria to exploit ephemeral microscale nutrient patches[J]. Proceedings of the National Academy of Sciences of the United States of America, 2008, 105(11): 4209-4214. DOI:10.1073/pnas.0709765105
[18] Mori J F, Ueberschaar N, Lu S P, et al. Sticking together: inter-species aggregation of bacteria isolated from iron snow is controlled by chemical signaling[J]. The ISME Journal, 2017, 11(5): 1075-1086. DOI:10.1038/ismej.2016.186
[19] Bižić-Ionescu M, Ionescu D, Grossart HP. Organic particles: heterogeneous hubs for microbial interactions in aquatic ecosystems[J]. Frontiers in Microbiology, 2018, 9. DOI:10.3389/fmicb.2018.02569
[20] Lyons M M, Ward J E, Gaff H, et al. Theory of island biogeography on a microscopic scale: organic aggregates as islands for aquatic pathogens[J]. Aquatic Microbial Ecology, 2010, 60(1): 1-13.
[21] Thiele S, Fuchs B M, Amann R, et al. Colonization in the photic zone and subsequent changes during sinking determine bacterial community composition in marine snow[J]. Applied and Environmental Microbiology, 2015, 81(4): 1463-1471. DOI:10.1128/AEM.02570-14
[22] Böckelmann U, Manz W, Neu T R, et al. Characterization of the microbial community of lotic organic aggregates ('river snow') in the Elbe River of Germany by cultivation and molecular methods[J]. FEMS Microbiology Ecology, 2000, 33(2): 157-170.
[23] Tang X M, Gao G, Qin B Q, et al. Characterization of bacterial communities associated with organic aggregates in a large, shallow, eutrophic freshwater lake (Lake Taihu, China)[J]. Microbial Ecology, 2009, 58(2): 307-322. DOI:10.1007/s00248-008-9482-8
[24] 晁建颖, 高光, 汤祥明, 等. 生态系统结构对太湖有机聚集体理化性质的影响[J]. 环境科学, 2009, 30(11): 3196-3202.
Chao J Y, Gao G, Tang X M, et al. Effects of ecosystem structure on the organic aggregates physical and chemical parameters of Lake Taihu[J]. Environmental Science, 2009, 30(11): 3196-3202. DOI:10.13227/j.hjkx.2009.11.021
[25] Lundgreen R B C, Jaspers C, Traving S J, et al. Eukaryotic and cyanobacterial communities associated with marine snow particles in the oligotrophic Sargasso Sea[J]. Scientific Reports, 2019, 9. DOI:10.1038/s41598-019-45146-7
[26] Grossart H P, Berman T, Simon M, et al. Occurrence and microbial dynamics of macroscopic organic aggregates (lake snow) in Lake Kinneret, Israel, in fall[J]. Aquatic Microbial Ecology, 1998, 14(1): 59-67.
[27] Allgaier M, Grossart HP. Seasonal dynamics and phylogenetic diversity of free-living and particle-associated bacterial communities in four lakes in northeastern Germany[J]. Aquatic Microbial Ecology, 2006, 45(2): 115-128.
[28] Ho J C, Michalak A M, Pahlevan N. Widespread global increase in intense lake phytoplankton blooms since the 1980s[J]. Nature, 2019, 574(7780): 667-670. DOI:10.1038/s41586-019-1648-7
[29] Hou X J, Feng L, Dai Y H, et al. Global mapping reveals increase in lacustrine algal blooms over the past decade[J]. Nature Geoscience, 2022, 15(2): 130-134. DOI:10.1038/s41561-021-00887-x
[30] 朱广伟, 秦伯强, 张运林, 等. 2005-2017年北部太湖水体叶绿素a和营养盐变化及影响因素[J]. 湖泊科学, 2018, 30(2): 279-295.
Zhu G W, Qin B Q, Zhang Y L, et al. Variation and driving factors of nutrients and chlorophyll-a concentrations in northern region of Lake Taihu, China, 2005-2017[J]. Journal of Lake Sciences, 2018, 30(2): 279-295.
[31] 秦伯强, 王小冬, 汤祥明, 等. 太湖富营养化与蓝藻水华引起的饮用水危机——原因与对策[J]. 地球科学进展, 2007, 22(9): 896-906.
Qin B Q, Wang X D, Tang X M, et al. Drinking water crisis caused by eutrophication and cyanobacterial bloom in Lake Taihu: cause and measurement[J]. Advances in Earth Science, 2007, 22(9): 896-906.
[32] Xu H, Mccarthy M J, Paerl H W, et al. Contributions of external nutrient loading and internal cycling to cyanobacterial bloom dynamics in Lake Taihu, China: implications for nutrient management[J]. Limnology and Oceanography, 2021, 66(4): 1492-1509. DOI:10.1002/lno.11700
[33] Xie G J, Tang X M, Shao K Q, et al. Bacterial diversity, community composition and metabolic function in Lake Tianmuhu and its dammed river: effects of domestic wastewater and damming[J]. Ecotoxicology and Environmental Safety, 2021, 213. DOI:10.1016/j.ecoenv.2021.112069
[34] 祁红娟, 龙麟, 周俊. 天目湖水生态环境变化趋势及影响因素研究[J]. 环境保护与循环经济, 2017, 37(11): 44-47.
[35] 杨文斌, 段文秀, 崔扬, 等. 天目湖沙河水库水生态安全状况长期变化及影响因素[J]. 环境科学, 2021, 42(10): 4739-4752.
Yang W B, Duan W X, Cui Y, et al. Long-term changes and drivers of ecological security in Shahe Reservoir, China[J]. Environmental Science, 2021, 42(10): 4739-4752.
[36] 王苏民, 窦鸿身. 中国湖泊志[M]. 北京: 科学出版社, 1998.
[37] 杨春艳, 李天翠, 施艳峰, 等. 抚仙湖近十年水质变化特征及保护措施初探[J]. 人民长江, 2021, 52(11): 47-55.
Yang C Y, Li T C, Shi Y F, et al. Study on variation characteristics of water quality and protection measures for Fuxian Lake in past decade[J]. Yangtze River, 2021, 52(11): 47-55.
[38] 郑田甜, 赵祖军, 赵筱青, 等. 云南星云湖水质变化及其人文因素驱动力分析[J]. 湖泊科学, 2018, 30(1): 79-90.
Zheng T T, Zhao Z J, Zhao X Q, et al. Water quality change and humanities driving force in Lake Xingyun, Yunnan Province[J]. Journal of Lake Sciences, 2018, 30(1): 79-90.
[39] 刘绍俊, 施艳峰, 翟竟余, 等. 星云湖微囊藻密度周年变化及其与水质指标的关系[J]. 环境化学, 2021, 40(7): 2064-2072.
Liu S J, Shi Y F, Zhai J Y, et al. Annual variations of Microcystis density and their relationships with water quality indices in Xingyun Lake[J]. Environmental Chemistry, 2021, 40(7): 2064-2072.
[40] 金相灿, 屠清瑛. 湖泊富营养化调查规范[M]. (第二版). 北京: 中国环境科学出版社, 1990.
[41] 王明翠, 刘雪芹, 张建辉. 湖泊富营养化评价方法及分级标准[J]. 中国环境监测, 2002, 18(5): 47-49.
Wang M C, Liu X Q, Zhang J H. Evaluate method and classification standard on lake eutrophication[J]. Environmental Monitoring in China, 2002, 18(5): 47-49.
[42] Tang X M, Gao G, Chao J Y, et al. Dynamics of organic-aggregate-associated bacterial communities and related environmental factors in Lake Taihu, a large Eutrophic shallow lake in China[J]. Limnology and Oceanography, 2010, 55(2): 469-480.
[43] Gong Y, Tang X M, Shao K Q, et al. Dynamics of bacterial abundance and the related environmental factors in large shallow Eutrophic Lake Taihu[J]. Journal of Freshwater Ecology, 2017, 32(1): 133-145.
[44] Xie G J, Tang X M, Gong Y, et al. How do planktonic particle collection methods affect bacterial diversity estimates and community composition in Oligo-, Meso- and Eutrophic lakes?[J]. Frontiers in Microbiology, 2020, 11. DOI:10.3389/fmicb.2020.593589
[45] Borcard D, Gillet F, Legendre P. Numerical ecology with R (2nd ed.)[M]. Cham: Springer, 2018.
[46] 邓建明, 秦伯强, 王博雯. 广义可加模型在R中的快捷实现及蓝藻水华预测分析中的应用[J]. 生态学杂志, 2015, 34(3): 835-842.
Deng J M, Qin B Q, Wang B W. Quick implementing of generalized additive models using R and its application in blue-green algal bloom forecasting[J]. Chinese Journal of Ecology, 2015, 34(3): 835-842.
[47] Xie G J, Tang X M, Shao K Q, et al. Spatiotemporal patterns and environmental drivers of total and active bacterial abundances in Lake Taihu, China[J]. Ecological Indicators, 2020, 114. DOI:10.1016/j.ecolind.2020.106335
[48] Wood S N. Generalized additive models: an introduction with R (2nd ed.)[M]. Boca Raton: CRC Press, 2017.
[49] Walsh C, Nally R M. hier. part: hierarchical partitioning. version 1.0-6[EB/OL]. https://CRAN.R-project.org/package=hier.part, 2020-03-03.
[50] Kiørboe T, Jackson G A. Marine snow, organic solute plumes, and optimal chemosensory behavior of bacteria[J]. Limnology and Oceanography, 2001, 46(6): 1309-1318.
[51] Zhou L, Zhou Y Q, Tang X M, et al. Resource aromaticity affects bacterial community successions in response to different sources of dissolved organic matter[J]. Water Research, 2021, 190. DOI:10.1016/j.watres.2020.116776
[52] Cai J, Bai C R, Tang X M, et al. Spatiotemporal variations in seston C∶N∶P stoichiometry in a large eutrophic floodplain lake (Lake Taihu): do the sources of seston explain stoichiometric flexibility?[J]. Water, 2020, 12(1). DOI:10.3390/w12010036
[53] Grossart H P, Simon M, Logan B E. Formation of macroscopic organic aggregates (lake snow) in a large lake: the significance of transparent exopolymer particles, plankton, and zooplankton[J]. Limnology and Oceanography, 1997, 42(8): 1651-1659.
[54] Ghiglione J F, Conan P, Pujo-Pay M. Diversity of total and active free-living vs. particle-attached bacteria in the euphotic zone of the NW Mediterranean Sea[J]. FEMS Microbiology Letters, 2009, 299(1): 9-21.
[55] Rink B, Martens T, Fischer D, et al. Short-term dynamics of bacterial communities in a tidally affected coastal ecosystem[J]. FEMS Microbiology Ecology, 2008, 66(2): 306-319.
[56] Lemarchand C, Jardillier L, Carrias J F, et al. Community composition and activity of prokaryotes associated to detrital particles in two contrasting lake ecosystems[J]. FEMS Microbiology Ecology, 2006, 57(3): 442-451.
[57] Urvoy M, Labry C, L'Helguen S, et al. Quorum sensing regulates bacterial processes that play a major role in marine biogeochemical cycles[J]. Frontiers in Marine Science, 2022, 9. DOI:10.3389/fmars.2022.834337
[58] Callieri C, Hernández-Avilés J S, Eckert E M, et al. Transparent exopolymer particles (TEP), phytoplankton and picocyanobacteria along a littoral-to-pelagic depth-gradient in a large subalpine lake[J]. Journal of Limnology, 2021, 80(3). DOI:10.4081/jlimnol.2021.2021
[59] Keskitalo K H, Bröder L, Jong D, et al. Seasonal variability in particulate organic carbon degradation in the Kolyma River, Siberia[J]. Environmental Research Letters, 2022, 17(3). DOI:10.1088/1748-9326/ac4f8d
[60] Tang X M, Chao J Y, Chen D, et al. Organic-aggregate-attached bacteria in aquatic ecosystems: abundance, diversity, community dynamics and function[A]. In: Marcelli M(Ed.). Oceanography[M]. Rijeka: InTech-Open Access Company, 2012.
[61] He Y K, Bai M H, He Y D, et al. Suspended particles are hotspots for pathogen-related bacteria and ARGs in coastal beach waters of northern China[J]. Science of the Total Environment, 2022, 817. DOI:10.1016/j.scitotenv.2022.153004
典型湖泊有机聚集体时空特征及驱动因子
谢贵娟, 龚伊, 朱富成, 刘昌利, 卢宝伟, 邓辉, 汤祥明