微生物降解及自团聚作用对鄱阳湖流域瑶湖水中不同粒径胶体荧光物质变化的影响*

黄 娴,聂明华,丁明军,晏彩霞,刘慧慧,吕泽兰

(江西师范大学地理与环境学院, 鄱阳湖湿地与流域研究教育部重点实验室, 南昌 330022)

天然胶体(natural colloids, NCs)通常是指一维粒径介于1 nm～1 μm之间的物质[1]. NCs广泛存在于水环境中，由于其颗粒小、比表面积大、活性基团丰富，因此NCs能够在环境中长距离迁移，并深刻影响着环境中重金属、药物等污染物的赋存分配、迁移转化等环境行为[2-4]. 湖泊水体中的NCs主要为外源或内源输入，其主导因素分别为人类活动以及湖泊水体中微生物活动、水动力条件等[5]. 相较于天然湖泊，城市内湖一般水深较浅，因此NCs受水体扰动影响更大，湖泊底泥再悬浮作用、湖水自身高生产力等使得城市内湖水体中NCs丰度更高[6-7]. 加之城市内湖水体通常污染较重，受人为影响较大[8-9]. 因此，研究城市内湖水中NCs对了解湖水环境中污染物的迁移转化意义重大.

受气候、水动力条件、外源输入等的影响，不同湖泊中的NCs化学组成差异较大[10-11]. 并且，NCs物质组成的差异会进一步影响其在湖泊水环境中的持久性. 如有研究发现，随着水流滞留时间的增加，外源性NCs降解较快，但内源NCs却能够长期存在[10-11]. 但亦有研究报道，在湖泊中，由于内源NCs比外源NCs的C/N比高、碳水化合物更多，从而更容易被降解并参与到微食物网中[12-13]. 湖泊水体中的NCs主要通过微生物降解、光化学降解以及絮凝沉淀作用发生转化[14-16]. 由于光化学降解主要发生在水质较好的水体或典型有色溶解性有机质中[17-18]，因此在水质相对较差的城市内湖中微生物降解则是NCs在湖水中转化的重要途径. 目前国内相关研究主要关注太湖等富营养湖泊中微生物对NCs降解的影响[13,19]. 这类湖泊中藻类、水生植物生长旺盛，水中NCs的类蛋白比例较高，其NCs降解过程必然不同于类腐殖质比例较高的城市内湖水中NCs. 然而，目前有关城市内湖水体中微生物对NCs物质组成降解特性影响的研究还较少，尤其是有关不同粒径NCs荧光特性的研究更待开展. 除微生物降解过程，自团聚过程导致的絮凝沉淀是湖泊NCs转化的另一重要途径之一，该过程对碳循环、养分动态变化等起到重要作用，近年来引起广泛关注[20-21]. 但已有研究较少关注该过程对NCs理化特性(如粒径、组成及不同粒级分配关系等)的影响，并且在研究过程中常常与微生物降解过程分开考虑，这显然与事实不符[22]. 因此，综合考虑微生物降解(生物降解)和自团聚过程(非生物降解)对湖泊水体中NCs理化性质(如不同粒径NCs荧光特性)的影响，对于深入了解水环境中NCs迁移转化等环境行为具有重要意义.

瑶湖是南昌市区内最大的封闭型浅水湖，受水产/畜禽养殖、城镇化过程等人为活动影响，水体污染严重[23-24]. NCs是环境中污染物的重要载体，对瑶湖水中NCs环境行为的研究有助于深入了解瑶湖水污染问题. 然而，目前有关瑶湖水质的研究主要集中于氮、磷等常规污染物的检测[25]，而针对水中NCs的迁移转化及其影响下的污染物环境行为等相关研究严重不足. 基于此，本文选取瑶湖水体及注入瑶湖的一条黑臭河流水体为研究对象，利用切向超滤技术及三维荧光光谱法(3D-EEM)，并结合平行因子分析法，对比不同粒径NCs荧光物质在灭菌和有微生物存在下各自的黑暗降解情况，以探究微生物降解和自团聚过程对湖泊水体中不同粒径NCs荧光物质降解的影响，从而为湖泊水环境中NCs迁移转化及其对污染物环境行为的影响研究提供数据支持.

1 材料与方法

1.1 水样的采集与处理

于2018年3月采集江西省南昌市最大的淡水湖泊——瑶湖水体(L)及注入瑶湖的一条黑臭河流水体(R)，采样点分布见图1. 水样采自表层30 cm内水体，采样的同时测定了水体中温度、溶解氧等基本理化参数(表1). 水样运回实验室后立即用经过400℃煅烧后的1 μm和0.45 μm玻璃纤维滤膜(PAUL, 美国)先后过滤分别得到含1 μm和0.45 μm粒径NCs的水样. 接着将含0.45 μm粒径NCs的水样利用切向超滤技术得到真溶解态水样(<1 kDa) (内置1 kDa超滤膜，聚醚砜滤膜，5×10-3m2，Vivaflow 50，赛多利斯，德国). 具体超滤方法参照作者前期研究[3]，超滤倍数(cf)为10. 因此，通过上述过滤、超滤处理，得到L和R含4种粒径NCs的水样，分别记为W原水、W1 μm、W0.45 μm和W1 kDa.

图1 瑶湖水体样品采集位置分布Fig.1 The distribution of sampling sites of waters from Lake Yaohu

表1 水样理化参数

1.2 培育实验

所有样品的培育过程均于暗室中进行以去除光降解对实验的影响. 首先，将L和R水样含4种粒径NCs的水体(即W原水、W1 μm、W0.45 μm和W1 kDa，各100 mL)分为2组，Ⅰ组为加入0.2 g/L NaN3的灭菌组，Ⅱ组为未做灭菌处理的对照组，实验设置2组平行样. 灭菌组反映的是自团聚过程对体系中荧光特性的影响，而对照组反映的是微生物降解及自团聚过程共同作用对体系中荧光特性的影响. 待获得含所有粒径NCs的水样后，立即测定其3D-EEMs(初始值). 接着加入NaN3，并开始计时，分别于第3、8、15、30天取样测定3D-EEMs. 除了测定3D-EEMs图谱外，同时还测定每一时间间隔对照组中W1 μm及W1 kDa水样的细菌丰度及溶解性有机碳(DOC)浓度[3,26].

1.3 三维荧光光谱测定

所有水样的三维荧光光谱(3D-EEMs)数据使用FLS 980荧光光谱检测器测定，仪器配备了450 W氙灯(爱丁堡, 英国). 激发(Ex)和发射(Em)波长范围分别为200～450、250～600 nm，增量均为5 nm. 狭缝宽度为4 nm，滞留时间为0.1 s. 考虑到内滤波效应，利用Milli-Q超纯水对R水样及部分L水样进行了不同浓度稀释，稀释倍数依据水样在254 nm处的吸光度是否小于0.3而定[27].

1.4 数据处理

1.4.1 平行因子分析(PARAFAC) 在进行PARAFAC处理前，对所有3D-EEMs数据利用Matlab 2012a进行瑞利和拉曼散射校正[28]. 同时利用Milli-Q超纯水在Ex=350 nm下的拉曼峰进行水拉曼单位校正(r.u.)[29]. PARAFAC模型通过Matlab 2012a软件及DOMFluor工具箱对选取的80个(2种水样×4种粒径×5个时间段×2个平行样)3D-EEMs数据进行分析[30]. 分析的可靠性通过对半分析及残差、载荷分析得以保证，模型逐步从3个组分数递增到6个组分进行检验，结果表明4组分数能较好地通过对半检验、残差及载荷分析. 模型输出的最大荧光强度(Fmax)值即为各样品对应的各荧光组分强度[31]. 对输出的PARAFAC组分通过在线OpenFluor数据库平台进行物质判定(http://www.openfluor.org)[32].

通过上一分级粒径水样的荧光强度扣除下一级粒径水样的荧光强度得到相应粒径NCs中荧光组分的荧光强度，以0.45～1 μm 粒径NCs为例(公式(1))：

Intensity(0.45～1 μm)=Intensity (<1 μm)-Intensity (<0.45 μm)

(1)

1.4.2 荧光指数 以Fn(280)及Fn(355)作为水样中类蛋白及类腐殖质的浓度[33]. 具体计算方法如下：Fn(280)为Ex=280 nm、Em=340～360 nm时的最大荧光强度值；Fn(355)为Ex=355 nm、Em=440～470 nm时的最大荧光强度值. 以上所有荧光强度均为散射校正、超纯水空白校正后的结果.

腐殖化指数(humification index,HIX)最早用于评价土壤中溶解性有机质的成熟度时引入[34].HIX定义为Em=254 nm、Ex为300～345 nm和435～480 nm波长下峰面积的比值，可表示有机质组分的腐殖化程度随降解过程的变化情况.

1.4.3 降解动力学 水样中荧光物质的降解率(R, %)的计算方法为：

R=C/C0×100%

(2)

式中,C和C0分别为荧光组分的荧光强度在某一降解时刻值及初始值.

含不同粒径NCs水样的PARAFAC组分及荧光指数随时间的降解可采用准一级反应动力学进行模拟，公式为：

ln (C/C0)=-k·t

(3)

式中,k为降解速率(d-1),t为降解时间(d).

2 结果与讨论

2.1 PARAFAC组分特征分析

通过PARAFAC模型，共解析出4个荧光组分. PARAFAC组分的3D-EEM图谱及最大Ex、Em波长、物质组成及信息描述等见图2及表2. C1(335/420 nm)和C4(290/340 nm)在OpenFluor数据库中的匹配度较高，说明C1和C4为环境中常见荧光物质. 其中C1为可见富里酸[35]，与表层沉积物及悬浮颗粒物中水溶性有机质的人为源类腐殖物质相似[36]；C4组分与河流沉积物来源的类蛋白物质相似[37]，为典型海洋中微生物源类蛋白物质[38-39]. 由检验模型匹配度的收敛系数(TCC, Tucker Convergence Coefficient)可知，数据库中可与C2(370/455 nm)和C3(395/515 nm)匹配的模型数较少(TCC值较小)，其中C2为土壤富里酸[40-41]，C3为类腐殖质[42].

第二，旅游危机事件本身是网络舆情主体作用的客体。旅游危机事件是旅游危机事件网络舆情的核心构成部分。旅游危机事件网络舆情的演进首先取决于旅游危机事件本身的发展和处置进展，因此旅游危机事件是决定线上网络舆情传播的重要因素。

图2 PARAFAC组分3D-EEM图及其最大Ex/Em位置Fig.2 3D-EEM and the maximum excitation/emission wavelength positions of PARAFAC components

表2 OpenFluor数据库中PARAFAC组分的最大TCC值*

*时间：2019年4月14日(UTC).

瑶湖水体(L)及其支流水体(R)初始时刻含不同粒径NCs水样的荧光强度分布见图3. 两种水样的荧光强度均随粒径减小而减小，这说明所有粒径分级水样中均含有荧光物质. 所有粒径分级水样均以类腐殖物质为主(C1～C3)，类蛋白物质较少(C4)，尤其是R水样尤为明显. 对于L水样，含不同粒径NCs水体的C1～C4组分荧光强度相差不大，但随粒径减小各组分荧光强度减小的情况不同. W原水→W1 μm、W1 μm→W0.45 μm各组分荧光强度变幅相差不大(均值0.03 r.u.)，但W0.45 μm→W1 kDa各组分荧光强度的变幅较大(均值0.045 r.u.). R水样中，总荧光强度明显高于L水样，说明支流水体中荧光物质较多，可能是瑶湖水体荧光物质的来源之一. 但所有粒径分级水样C1～C4组分荧光强度的变化规律与L水样一致，即W原水→W1 μm、W1 μm→W0.45 μm各组分荧光强度变幅较小(0.01～0.02 r.u.)，W0.45 μm→W1 kDa水样则高达0.08 r.u.. 以上结果说明大部分的荧光物质主要分布在1 kDa～0.45 μm粒径范围内，而这一粒径范围正是典型胶体物质部分，这与前人研究结果一致[43].

图3 初始时刻瑶湖水体(L)和支流水体(R)含不同粒径NCs水样荧光组分的分布Fig.3 Distribution of fluorescence intensity in water samples with different sizes of NCs from Lake Yaohu(L) and its tributaries(R) samples at the initial moment

2.2 荧光组分的变化

2.2.1 不同粒径NCs中荧光组分荧光强度的变化 L及R水样不同粒径NCs中C1～C4组分随时间的变化如图4所示. 对于L水样，无论有无微生物存在，在培育的第3天各荧光组分总荧光强度均出现小幅上升，随后开始明显下降(第3～15天)，但培育15天后荧光强度下降速率减慢. 对于R水样，除未出现在培育第3天荧光强度小幅上升的情况，其他变化规律与L水样一致. 总体来说，R水样对照组(有微生物)中荧光组分的减少比灭菌组多.

微生物降解对不同粒径NCs中荧光物质变化的影响在L水样不明显，但在R水样影响较大(除C1组分外). C1～C4组分在不同粒径NCs中的减少以小粒径NCs(<1 kDa和1 kDa～0.45 μm)较多、大粒径NCs较少(0.45～1 μm和>1 μm)，在R水样中甚至出现中、小粒径(1 kDa～0.45 μm或0.45～1 μm) NCs荧光强度增加的情况(R水样对照组C2、C3、C4组分). 该现象在胶体尺度的相关研究中亦有发现，即颗粒态(>0.45 μm、10 kDa～0.45 μm、1～10 kDa)荧光物质的荧光强度随时间不仅没有减少，反而呈增加趋势，而真溶解态(<1 kDa)中荧光强度明显减少，该研究认为这可能是因为微生物影响下真溶解态水样中的荧光物质转移至胶体甚至是微米级颗粒中所致[44]. 该现象同时说明支流水体中微生物对不同粒径NCs荧光物质的微生物降解作用比湖泊水体显著，这可能与支流水体较高的细胞丰度有关(图5)，同时也说明了支流水体微生物可能是湖泊水体微生物的重要来源，但相关结果有待进一步验证.

图4 含不同粒径NCs瑶湖水体(L)和支流水体(R)在微生物存在及灭菌条件下C1～C4组分的降解情况 (图中红色虚线下方表示相应组分荧光强度的增加)Fig.4 Degradation of fluorescence component in water samples with different size of NCs from Lake Yaohu(L) and its tributaries(R) samples with and without microorganisms

2.2.2 不同粒径NCs荧光组分的降解动力学 L和R水样灭菌组和对照组在不同培育时间下的C1～C4组分荧光强度C与初始值C0的比值(降解率)随时间变化趋势如图6所示. 总体来说，L水样的荧光物质在实验过程中共减少了50%左右(灭菌组与对照组差别不大)，而R水样的荧光物质在灭菌组与对照组中分别减少了50%和70%. 这说明L水样中自团聚作用是导致其荧光物质减少的主要原因，而R水样中荧光物质的减少则是自团聚作用及微生物降解共同作用的结果(自团聚作用占50%，微生物降解占20%). 此外，类腐殖质的降解率(即C1、C2和C3组分降解率的平均值)与总体降解率(C1～C4所有组分降解率的平均值)相近并略低于总体降解率，这说明L和R水样中荧光物质的减少主要表现为类腐殖物质的减少. 而类蛋白质的降解率在灭菌组中与类腐殖质的降解率相差不大，但在对照组中则明显高于类腐殖质的降解率，且随粒径减小二者的降解率差值逐渐减小，即W原水为20%、W1 μm为15%～18%、W0.45 μm为10%、W1 kDa则小于5%. 这说明类蛋白物质主要分布在大粒径水体中，且主要受微生物降解的影响.

L水样中所有粒径荧光组分在培育第3天时均出现了一个明显的小幅上升(该结果与前文分析一致，图4). 这可能与降解初期大粒径(W1 μm)水体及小粒径(W1 kDa)水体中的细胞数均显著增加有关，约为初始值的1.5倍(图5折线图)，这与文献报道结果一致[44]. 含不同粒径NCs水样中各组分在灭菌组中的降解差别不大，说明自团聚过程在所有粒径水样中均存在，且差别不大. 但对照组中各荧光组分的变化较明显，其中以W1 kDa降解较缓慢、W0.45 μm次之，这与对照组中W1 kDa的DOC浓度降解较慢的趋势一致(图5柱状图). 结合灭菌组结果可知，大粒径水体中除了有自团聚过程导致的荧光组分减少外，还有微生物降解的作用，而小粒径水体中有机碳及其荧光物质的减少则主要为自团聚作用占主导.

图5 瑶湖水体(L)和支流水体(R)对照组中W1 μm及W1 kDa细胞丰度及DOC浓度的变化过程Fig.5 Changes of the bacterial abundance and dissolved organic carbon concentration in <1 μm and <1 kDa water sample from Lake Yaohu(L) and its tributaries(R) samples with microorganisms

图6 瑶湖水体(L)和支流水体(R)含不同粒径NCs水样荧光组分 在微生物存在及灭菌条件下的降解率Fig.6 Degradation rates of fluorescence components in water samples with different size of NCs from Lake Yaohu(L) and its tributaries(R) samples with and without microorganisms

进一步分析发现，大部分荧光物质的降解符合准一级动力学模型(除R水样W1 kDa的C2、C3、C4组分外)(表3)，且拟合R2值有随水样粒径减小而减小的趋势. 总体来说，含相同粒径NCs的水样、同一荧光物质的降解以对照组中的速率更快，但W1 kDa差异相对较小. 由于细菌、病毒等微生物的粒径基本上大于1 kDa[1]，因此W1 kDa可认为是不含微生物的“灭菌”水体. 由此可知，大粒径水体中荧光物质的减少受微生物降解的影响比小粒径大，而小粒径水体中荧光物质的降解主要是自团聚作用占主导，这与前文分析结果一致. 从细节来看，L水样灭菌组随着粒径减小，其荧光组分的降解速率(k)呈微小波动变化；对照组中则是随粒径减小k值随之减小(表3)，这是因为微生物降解作用随粒径减小而减小. 并且，灭菌组和对照组中均以C4组分的k值最高，尤其是对照组更显著，这说明C4组分代表的类蛋白组分降解更快. 对于R水样由于W1 kDa中除C1组分外，均不符合准一级动力学模型(R2<0.5)，因此未给出相应k值. 但根据已有数据可知其变化规律与L水样一致，即随水体粒径减小荧光物质的降解速率逐渐降低. 与L水样类似，R水样在对照组中C4组分的降解速率较高，说明微生物对类蛋白物质的降解较快. 这与上文C4组分荧光强度的减小比例分析结果一致，与胶体尺度范围内类蛋白质、碳水化合物等大分子物质主要是微生物降解占主导、类腐殖质等小分子物质主要是自团聚作用占主导的结论亦一致，研究表明这一现象与其分子量及碳组成有关[44].

表3 瑶湖水体(L)和支流水体(R)含不同粒径NCs水样在微生物存在及 灭菌条件下各参数降解速率k值及R2值

“-”表示R2<0.5.

2.3 类蛋白及类腐殖物质浓度的变化

图7 瑶湖水体(L)和支流水体(R)含不同粒径NCs水样中Fn(280)和 Fn (355)在微生物存在及灭菌条件下的降解率Fig.7 Degradation rates of Fn (280) and Fn (355) in water samples with different sizes of NCs from Lake Yaohu(L) and its tributaries(R) samples with and without microorganisms

2.4 腐殖化指数HIX值的变化

随培育过程的进行，L和R水样含不同粒径NCs水体的腐殖化程度总体来说呈下降趋势(图8). 并且相较于灭菌组，对照组的腐殖化程度变化较大，即中、大粒径(W原水、W1 μm、W0.45 μm)水体的腐殖化程度随时间呈先增加后减小的趋势，而小粒径(W1 kDa)水体的腐殖化程度随时间始终呈减小趋势. 这可能是因为在降解初期由于微生物量的增加，增加了将脂肪碳转化为芳香结构的能力，进而使得腐殖化程度增加[45-46]. 两种水样含不同粒径NCs水体HIX值随时间的变化过程均符合准一级动力学模型(除R水样对照组W原水外)(表3). 对于L水样，HIX值的降解速率k值在含各粒径NCs水体中的差异较小. 但R水样中明显的存在随水体粒径减小，HIX值的k值增加的趋势，表明小粒径水体的腐殖化程度减小较快、大粒径水体较慢，这可能与大粒径水体在降解初期腐殖化程度略有增加有关.

图8 瑶湖水体(L)和支流水体(R)含不同粒径NCs水样中HIX在微生物存在及灭菌条件下的降解率Fig.8 Degradation rates of HIX in water samples with different sizes of NCs from Lake Yaohu(L) and its tributaries(R) sample with and without microorganisms

3 结论

1)通过PARAFAC模型，瑶湖(L)及其支流(R)水样含不同粒径NCs水体荧光物质随时间变化共解析出4个荧光组分，即C1～C3为类腐殖物质、C4为类蛋白物质.

2)在初始时刻，两种水样均以类腐殖物质较多、类蛋白物质较少，且均以小粒径NCs(1 kDa～0.45 μm)中的荧光强度更高，但总体来说以R水样荧光组分的荧光强度较高.

3)两种水样含各粒径NCs水体荧光物质随时间的变化过程表明：①天然水体中荧光物质随时间的减少是微生物降解和自团聚作用共同作用的结果；②由于真溶解态中荧光物质的降解主要受自团聚作用影响，因此在灭菌和对照组中的差别较小，但在中、大粒径水体中荧光物质变化是微生物降解和自团聚作用共同作用结果，因此总体降解速率更快，尤其是类蛋白组分.

4)荧光参数分析结果表明：①类蛋白物质的降解主要受微生物作用，尤其是大粒径水体中的降解速率较快，并出现类蛋白物质由大粒径向小粒径NCs转移的现象；②类腐殖物质的降解主要是自团聚作用占主导，并可能出现由小粒径向大粒径NCs的转移；③腐殖化指数随降解过程总体呈减小趋势，微生物作用下大粒径水体腐殖化程度呈先增后减趋势，小粒径水体则以较快速率始终呈减小趋势.