本研究使用的HA购自天津光复试剂公司。 取1 g HA加入锥形瓶中, 加入100 mL去离子水配制样品。 一组在室温(25 ℃), 150 r· min^-1恒温振荡箱内, 根据恒温振荡箱的工作时间上限设置振荡时间为0~72 h, 间隔为3 h; 另一组在4 000 r· min^-1高速离心机内, 根据离心机最大工作时间设置的离心时间为0~60 min, 间隔为5 min。 静置4 h后, 将每个试样的上清液经0.45 μm水系微孔滤膜过滤得到DOM提取液, 稀释至UV₂₆₀< 0.02以去除内滤效应。

采用荧光分光光度计(F4600)测量样品的三维荧光光谱。 激发光源为氘灯。 激发波长范围为220~450 nm, 步长5 nm, 发射波长200~500 nm, 步长5 nm, 狭缝宽度5 nm, PMT电压400 V, 扫描速度1 200 nm· min^-1, 在1 cm石英荧光比色皿中测量。

利用MATLAB 2018对光谱数据(分析前以Milli-Q超纯水为空白扣除瑞利散射和拉曼散射的影响)进行FRI和PARAFAC解析。 通过PARAFAC识别荧光组分前, 将荧光强度归一化为激发波长350 nm处的拉曼信号强度(RU₃₅₀), 用最大荧光强度(F_max)值评估荧光峰强度。 通过FRI方法得出不同区域的积分标准体积代表各区域对应的DOM相对含量^[13], FRI所划分的五个区域见表1, 相关计算公式见式(1)— 式(3)。 利用Origin2018将处理好的光谱数据绘制成等高线谱图。

Фi, n=MFi∫∫exemI(λexλem)dλexdλem(1)

ФT, n=∑i=15Фi, n(2)

Pi, n=Фi, nФT, n×100%(3)

式中: Ф _{i, n}为荧光区域i的积分标准体积; Ф _i为荧光区域i的积分体积; λ _ex为激发波长; λ _em为发射波长; I(λ _exλ _em)为激发发射波长对应的荧光强度; Ф _{T, n}为总荧光区域积分标准体积; P_{i, n}为荧光区域i的积分标准体积占总积分标准体积比例; MF_i为倍增系数, 为荧光区域i的积分面积占总的荧光区域积分面积比例的倒数。

表1

Table 1

表1(Table 1)

表1 荧光积分区域划分范围 Table 1 Fluorescent integral area 区域 所代表有机物类型 激发波长/nm 发射波长/nm Ⅰ 类腐殖酸 250~450 380~550 Ⅱ 类富里酸 220~250 380~550 Ⅲ 类色氨酸 220~250 330~380 Ⅳ 类络氨酸 220~250 280~330 Ⅴ 微生物代谢产物 250~450 280~380

表1 荧光积分区域划分范围 Table 1 Fluorescent integral area

三维荧光光谱可以指示可溶性有机物的种类、 性质和相对含量等丰富信息, 从而反映提取样品中DOM的变化情况。 振荡和离心提取着HA水解的三维荧光光谱分别见图1、 图2。 结果表明, 随着提取时间的变化, DOM的光谱形状基本无明显变化, 但荧光峰强度呈现上升到下降的往复过程。

图1

Fig.1

	Figure Option ViewDownloadNew Window
	图1 离心提取下HA三维荧光光谱图Fig.1 3D fluorescence spectrogram of humic acid under centrifugal extraction

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 振荡提取下HA三维荧光光谱图Fig.2 3D fluorescence spectrogram of humic acid under Oscillating extraction

图3为离心和振荡提取下光谱区域体积积分变化情况。 荧光区域的体积积分(Ф _i)可间接表征不同区域特定结构荧光组分的相对含量。 结果表明, 与离心相比, 振荡各组分的占比变化幅度较小; 在离心和振荡处理下类腐殖质组分的积分占比(P1+P2)最大, 分别占比53.41%~70.67%, 74.55%~81.77%; 类蛋白质组分占比(P3+P4)分别在25.84%~42.65%, 14.85%~21.96%; 而微生物代谢产物的积分占比(P5)均在3%左右, 说明振荡处理提取效果较稳定。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 两种提取方式下荧光区域积分值变化情况 (a): 离心提取方式下荧光区域积分值变化情况; (b): 振荡提取方式下荧光区域积分值变化情况Fig.3 Variation of fluorescence region integration values under two extraction methods (a): Changes of regional integral values in centrifugal extraction; (b): Changes of regional integral values in oscillatory extraction

通过拆半分析、 残差分析、 平方误差和、 解释方差和核一致性分析, 确定两种提取方式DOM平行因子分析模型的最佳因子数为3, 其荧光光谱图如图5。 表2为所有组分的激发和发射波长最大值的位置, 在Open Fluor数据库中查询到三个组分的多个强匹配项(tucker congruence coefficient, TCC> 0.95), 提供了荧光成分的近似描述。 C3组分的荧光峰最大值较C1组分发生了红移, 表现出更复杂的性质。 同C1相比, C3生物生产、 降解和吸附的强度更低^[14]。

表2

Table 2

表2(Table 2)

表2 HA的DOM荧光组分特征 Table 2 Characteristics of fluorescent component in humic acid DOM 组分 Ex/Em(nm) 物质及描述 文献中数值/nm C1 270(340)/440 生物可降解类腐殖质[14, 15] 265(360)/450 C2 245/410 微生物类富里酸类物质[16, 17] 235/408 C3 290 (405)/495 未定义类腐殖酸类物质[18, 19] (285/395)/497

表2 HA的DOM荧光组分特征 Table 2 Characteristics of fluorescent component in humic acid DOM

两种提取方式下F_max变化情况如图4所示, 振荡提取样品的F_max为C1> C2> C3组分, 离心提取样品的F_max为C2> C1> C3组分。 C1对整个模拟荧光得分的平均贡献率为44.03%, C2平均贡献率44.28%之间, C3平均贡献率为11.68%。

图4

Fig.4

	Figure Option ViewDownloadNew Window
	图4 两种提取方式下Fmax值变化情况 (a): 离心提取方式下Fmax值变化情况; (b): 振荡提取方式下Fmax值变化情况Fig.4 Variation of Fmax values under two extraction methods (a): Changes of regional integral values in centrifugal extraction; (b): Changes of regional integral values in oscillatory extraction

图5

Fig.5

	Figure Option ViewDownloadNew Window
	图5 基于PARAFAC解析出两种提取方式下HA的DOM荧光组分Fig.5 Fluorescence components of humic acid under the two extraction methods by PARAFAC model

提取过程中振荡处理的区域体积积分值和F_max值均大于离心, 且其占比的变化幅度小于离心处理。 说明振荡处理提取量充分且提取的组分相对较稳定。 因为高速离心机工作时间最大为60 min, 恒温振荡箱工作时间上限为4 500 min, 可以通过较长时间的振荡提取分子结构比较复杂的组分。 离心是利用转子高速旋转产生的强大离心力, 在较短时间内迫使这些微粒克服扩散产生沉降运动^[20], 多用来分离不同密度的混合物。

两种提取方法的区域积分值和F_max变化趋势基本一致, 呈现先上升后下降又恢复上升趋势, 直至上升幅度趋于平缓。 图6(a)显示, 类络氨酸与类腐殖质、 微生物代谢产物组分呈负相关, 与类蛋白质组分呈正相关。 可能是由于各组分之间存在复杂的相互作用, 导致类络氨酸变化规律与其他组分存在差异。 在提取前期DOM被不断释放含量上升; 出现下降时可能是出现絮凝现象, 当释放量大于絮凝量会出现上升, 反之则会下降。 同时, 大分子物质的絮凝能力大于小分子物质^[20], 这一点在振荡提取后期(72 h)中荧光区域Ⅲ 体积积分(Ф ₃)、 区域Ⅳ 体积积分(Ф ₄)较初期上升, 而区域Ⅰ 体积积分(Ф ₁)、 区域Ⅱ 体积积分(Ф ₂)则出现下降能够得到证实。

图6

Fig.6

	Figure Option ViewDownloadNew Window
	图6 HA DOM组分的相关性分析 注: * 表示存在显著性差异(p< 0.05), * * 表示存在极显著性差异(p< 0.01) (a): FRI方法下区域积分值相关性分析; (b): PARAFAC方法Fmax值相关性分析Fig.6 Correlation analysis of DOM components in humic acid (a): Correlation analysis of regional integral values by FRI method; (b): Correlation analysis of Fmax values by PARAFAC method

对PARAFAC和FRI分析得到的不同提取时间下的区域积分值和F_max进行相关性分析和主成分分析, 结果如图6、 图7所示。 DOM组分间的相关关系表明它们来源或其结构具有某种联系。 相关性结果显示, 离心处理下荧光区域积分Ф ₁, Ф ₂和Ф ₅均与总量呈现显著正相关(p< 0.05), Ф ₃, Ф ₄与总量的关系较弱; C1, C2和C3均与总量呈现显著正相关(p< 0.05)。 类腐殖质组分与总量的提取情况相似, 类蛋白质与总量的提取情况存在差异, 说明各组分之间的相互作用不同, 总量的提取情况无法概括组分的情况, 需要对组分和总量分别设置提取时间。 PARAFAC和FRI解析DOM的主成分1对总方差的解释分别达到67.9%和66.0%, 第2主成分对总方差的解释分别达到32.0%和25.7%, 累积分别达到99.9%和91.7%, 说明主成分分析法对原数据的损失不大。 两种解析方法下主成分1的主要影响因子均是类腐殖质组分(Ф ₂, C1), 表明类腐殖质组分与总量的提取情况相似。

图7

Fig.7

	Figure Option ViewDownloadNew Window
	图7 HA DOM 组分的主成分分析 (a): PARAFAC方法得到的Fmax值主成分分析; (b): FRI方法下区域积分值主成分分析Fig.7 Principle component analysis of DOM components in humic acid (a): Principle component analysis of Fmax values by PARAFAC method; (b): Principle component analysis of regional integral values under FRI method

分别根据区域积分值、 F_max以及经济角度考虑, 筛选出的最佳提取时间结果见表3。 以往研究中^[21], DOC的提取过程符合一级动力学模型, 与本研究中得到结果存在差异, 可能是本研究中提取时间较长, 存在分段符合动力学模型的情况。 由于本研究中每次组分提取量趋于稳定时的峰值相差较小, 从经济角度, 选择区域积分值和F_max初期出现峰值的时间。 此外, 本研究针对总量提取时间的选取, 为体现组分的提取效果, 以各组分占比变化趋于稳定为前提, 筛选区域积分值和F_max初期出现峰值的时间。 离心处理的组分最佳提取时间同理可得。 同时, 根据相关性分析和主成分分析发现, 类富里酸组分(区域Ⅱ 、 C2)与总量的相关性较高且是主成分1的主要影响因子, 二者具有相同的提取时间。 说明总量可与影响因子较大且相关性较强的组分选择相同的提取时间。 综上, 两种提取方法各有优点, 振荡提取的分子量大于离心处理, 提取组分比例较为稳定。 离心提取下总量与组分最佳提取时间可选为一致。

表3

Table 3

表3(Table 3)

表3 两种提取方法下最佳提取时间结果 Table 3 Best extraction time results under two extraction methods FRI PARAFAC 组分名称 离心/min 振荡/h 组分名称 离心/min 振荡/h P1 45 21 C1 45 21 P2 45 12 C2 45 12 P3 45 12 C3 45 39 P4 45 12 总量 45 12 P5 45 12 总量 45 12

表3 两种提取方法下最佳提取时间结果 Table 3 Best extraction time results under two extraction methods

PARAFAC分析和FRI方法对DOM组分的解析情况存在差异^[22], 本研究根据以上两种方法识别出的组分类别存在差异, 但是通过分析得到的最佳时间呈现较高的一致性。 虽然从表3来看, 振荡处理中PARAFAC分析识别出的C3为类腐殖酸组分, 与FRI方法中的区域Ⅰ 相对应, 但是二者的最佳提取时间不同。 这是由于FRI方法把一部分类腐殖酸类物质划分入了类富里酸物质^[23], 同时C3的发射波长也较同为类腐殖酸组分的C1发生红移, 结构更为复杂, 需要延长提取时间, 这一点与DOC在提取过程中得出的结论一致^[21]。

同时还必须强调本研究的局限性, 因为本研究只针对这一种外源有机物的最佳提取时间, 不同类别的外源有机物组分类别及分布存在差异, 组分与总量的关系是否存在差异, 需要后续进行比较研究。 此外, 在提取方案中对组分的最佳提取剂种类、 样品浓度等可参考本文开展进一步研究。