洱海悬浮颗粒物和表层沉积物有机碳氮同位素来源特征及水质指示意义*

季宁宁，刘 永，王圣瑞

(1：北京大学环境科学与工程学院，国家环境保护河流全物质通量重点实验室，北京 100871)

(2：粤港水安全保障联合实验室，北京师范大学珠海校区水科学研究中心，珠海 519087)

(3：北京师范大学水科学研究院，北京 100875)

(4：云南省高原湖泊流域污染过程与管理重点实验室，昆明 650034)

近几十年来, 湖泊富营养化和生态系统退化现象频繁发生，氮、磷等营养盐的过量输入是主要因素之一[1]. 湖泊生态系统中，沉积物和悬浮颗粒物有机质是营养盐的重要蓄积库，是营养盐周转循环的关键环节[2]. 自然条件下，湖泊沉积物和悬浮颗粒物有机质主要来源于水体内源自生(如水生植物和浮游生物等)及外源输入(如陆生植物碎屑和土壤有机质等). 不同来源有机质记录有机质的生物活性及环境效应归宿[3]，影响着水生态系统中营养盐的生物利用度及其生物地球化学循环, 如细菌和微藻类来源的有机质更易降解，促进营养盐循环[4]; 降雨增加陆源有机质输入，从而改变湖泊有机质来源组分. 因此，鉴别沉积物和悬浮颗粒物有机质来源是理解与富营养化相关的营养盐内循环过程的基础，有利于从源头上控制营养盐负荷.

随着同位素示踪技术的发展，研究湖泊有机碳、氮稳定同位素已成为揭示水生态系统中碳、氮来源的强有力手段，它们能有效地识别湖泊生态环境的演变过程和人为影响因素[5]. 湖泊有机碳同位素能有效指示水生态系统初级生产力变化和有机碳来源及其变化规律[6]. 氮同位素不仅可以示踪有机质来源，而且能够反映氮来源及指示相关的生物地球化学循环过程[7]. 如倪兆奎等[8]运用稳定碳、氮同位素技术，研究了洱海不同河流水体悬浮颗粒物中有机碳、氮的来源，并探讨了其与流域环境和人类活动之间的关系. Xu等[9]研究了抚仙湖和星云湖颗粒态有机碳、氮同位素分布特征，揭示了贫、富营养湖泊碳、氮同位素的差异；王毛兰等[10]利用碳、氮同位素研究了鄱阳湖及其支流表层沉积物有机质来源. 此外，不同来源有机质中C/N比值差异明显，如多数细菌等微生物C/N比值分布在2.6～4.3之间，浮游植物C/N比值范围为7.7～10.1，而高等植物C/N比值高达50以上. 一般来讲，有机质中C/N比值较高表明陆源性来源，而比值较低表明内源藻类或者浮游生物有机残体来源[11]. 因此，湖泊悬浮颗粒物和沉积物有机碳、氮稳定同位素和C/N比值时空变化的综合信息对于有效地指示水生态系统有机碳与氮来源及营养盐循环具有重要意义. 然而，已有研究多为单独研究悬浮颗粒物或沉积物中有机碳与氮来源[12]，而对比分析悬浮颗粒物和沉积物有机碳与氮来源的时空差异，并揭示其与水质响应关系的研究鲜有报道. 特别是在藻华期，湖泊内源有机质贡献增大，加之湖泊水环境改变，悬浮颗粒物和沉积物有机碳与氮同位素受到一定的影响，其水质指示意义未见报道.

近20年来，随着经济快速发展，人类生产活动，特别是农业的快速发展引起流域面源污染负荷快速增加，导致洱海水质恶化[13]，水体透明度明显下降，水生植物严重退化，洱海处于由中营养向富营养转变的关键期. 由于洱海流域人类活动变化及其自身水生态环境演变，悬浮颗粒物和沉积物有机碳与氮来源结构也发生了一定的改变，这些变化可能对洱海水质及藻华暴发产生重要影响. 因此，对洱海悬浮颗粒物和沉积物有机碳与氮的溯源研究是进行营养盐负荷源头控制的关键步骤.

本文试图分析洱海悬浮颗粒物与表层沉积物有机碳与氮稳定同位素的时空分布特征，建立端元混合模型并结合C/N比值定性和定量分析悬浮颗粒物与表层沉积物有机碳与氮来源及贡献率差异，探究其水质指示作用，并结合环境因子等信息，利用随机森林回归分析悬浮颗粒物与表层沉积物有机碳与氮来源特征对水质的影响，以期为有效地控制洱海污染负荷提供理论指导.

1 材料与方法

1.1 研究区域概况与样品采集

洱海 (25°36′～25°58′N，100°06′～100°18′E)是中国西南地区第二大淡水湖，是云南省大理市重要的饮用水源[14]. 其流域面积2565 km2，水面面积249 km2(图1)，平均水深为10.5 m，最大水深为20.7 m. 该流域属亚热带西南高原季风气候带，年均气温15.1℃，年均降水量1000 mm，95%降雨集中于5-10月，称为雨季，其余月份为旱季. 洱海流域入湖河溪共计117条，由弥苴河、罗时江、永安江构成的“北三江”水系，南部波罗江水系和西部“苍山十八溪”是洱海主要入湖通道，洱海唯一出湖河流为西洱河[15]，水流总体由北向南流动. 由于河流流入量高于流出量，且平均停留时间为2.75年，这意味着外源负荷和内部循环对水质变化至关重要[16]. 北部和中西部的农业及南部区域的人类活动是洱海污染的主要来源. 近10年来，洱海水质总体为《地表水环境质量标准》(GB 3838-2002)Ⅱ～Ⅲ类，但每年9-10月水华高风险期，洱海流域污染物随降雨集中入湖，在适宜的温度和光照等条件下，蓝藻水华时有发生[17]，特别在2013年暴发大规模水华，全湖平均叶绿素a(Chl.a)浓度高达39.6 μg/L，洱海处于富营养化初期阶段和治理的关键期.

图1 洱海流域及采样点

根据湖底地形可将洱海分为北部 (康廊-海舌以北)、中部 (挖色湖心为中心，龙龛-海舌)和南部 (龙龛以南) 3个湖区[18]. 按照代表不同时空特征的原则[19]，分别于2013年4月、2013年7月、2013年10月和2014年1月，利用GPS定位，同时采集北部、中部和南部的表层水体样品(水面下0.5 m)及表层沉积物样品(0～10 cm)(图1). 1月和4月为旱季，7月和10月为雨季. 采集洱海水样储存于聚乙烯瓶中低温避光保存，并于48 h内运回实验室分析. 所有水样通过预先500℃下灼烧 4 h 左右的玻璃纤维滤膜 (Whatman GF/F, 0.7 μm, 47 mm)过滤收集悬浮颗粒物，冷冻干燥，在4℃黑暗中保存直至分析. 利用抓斗式底泥采样器现场采集表层沉积物(0～10 cm)，放自封袋中标号低温保存，带回实验室后去除植物碎屑、碎石和贝壳等，放冷冻干燥机中真空冷干，然后取适量样品研磨过筛(孔径为0.075 mm)后置于自封袋中备用.

1.2 样品的处理与分析

1.2.2 碳氮元素分析及碳氮同位素测定 研磨过筛后的沉积物样品和冷干称重后悬浮颗粒物滤膜滴加1 mol/L盐酸完全浸湿，并置于含浓盐酸的干燥器中熏蒸48 h，去除样品中的无机碳，再用超纯水反复洗涤，直至pH值为7左右，再将其冷冻干燥[20]. 经上述处理后的样品在中国农业科学院农业环境与可持续发展研究所用元素分析仪(Flash EA 2000HT)和同位素质谱仪(Thermo Fisher Scientific Inc., USA)联用测定总有机碳、总氮含量及稳定同位素组成[21]，计算公式分别为:

δ13C (‰)=[(Rsample-Rstandard)/Rstandard]×1000,R=13C/12C

(1)

δ15N (‰)=[(Rsample-Rstandard)/Rstandard]×1000,R=15N/14N

(2)

式中，Rsample为待测样品同位素比值；Rstandard为标准参比样品同位素比值；13C/12C和15N/14N分别对应国际标准Vienna PDB与大气中的氮标准. 样品C/N比值可由测定的TOC与TN含量计算得到. δ13C和δ15N的平均分析误差约为0.13‰和0.1‰.

1.2.3 有机碳与氮来源端元贡献率估算 由于湖泊有机碳、氮同位素来源的复杂性，本研究利用端元平衡混合模型(Iso Source)，计算每种来源端元的贡献率(%). 端元混合模型是以质量守恒混合模型为基础，开发用于计算源头值超过n+1(n为多种同位素值)的潜在贡献比值[22]. 该模型在一定增量范围内，通过实测同位素信息，使用标准线性混合模型模拟不同端元物质所占比例[23]. 该模型的质量守恒等式为：

δm=fa·δa+fb·δb+fc·δc

(3)

1=fa+fb+fc

(4)

式中，δm代表样品中的同位素值；δa、δb、δc分别代表a、b、c 3种来源端元同位素分布范围的中间值作为端元值；fa、fb、fc分别代表a、b、c 3种端元物质相应的贡献百分比.

1.3 数据处理与分析

1.3.2 随机森林 随机森林是以决策树为基础，对数据进行分类和回归的新工具[26]. 随机森林的分析方法是通过对数据的每一个解释变量分别随机置换以得出每个解释变量对被解释变量的贡献程度，以均方误差(MSE)表示. 本研究利用随机森林回归模型(random forest analysis) 探究悬浮颗粒物和表层沉积物有机碳、氮同位素以及水环境因子对多营养盐循环指标和Chl.a浓度影响的贡献程度，MSE值越高表示该变量越重要. 利用R软件中“random forest” 数据包进行随机森林模型预测; 利用“rf Utilities”数据包对模型整体检验; 利用“rf Permute”数据包对随机森林中每个变量对模型贡献程度的显著性检验. 最后，使用Sigmaplot 14.0 软件进行科学绘图.

1.3.3 其他分析 利用Spearman相关系数评估有机碳、氮稳定同位素和水体营养盐指标之间关系. 使用非参数Mann-Whitney U检验比较营养盐及有机碳、氮同位素的时空差异性. 采用SPSS 22.0进行数据描述性统计分析. 采用OriginPro 2017等相关软件进行制图.

2 结果与讨论

2.1 洱海水质时空变化特征分析

研究区水体基本理化性质见表1，旱季ρ(DO)(均值为7.42 mg/L±0.08 mg/L)显著高于雨季(6.23 mg/L±0.29 mg/L)(P<0.05)，其在北、中、南湖区差异性不显著(P>0.05). ORP值在旱季为99.9～138.8 mV(均值为122.5 mV±14.2 mV)，氧化性较强；在雨季为-70.7～14.7 mV(均值为-4.3 mV±27.1 mV)，还原性较强. ORP影响水体营养盐的转化和有效性，较高的ORP有利于好氧微生物的生长，而较低的ORP则利于厌氧菌的繁殖，进而使水质变坏[27]. 洱海水体pH在旱、雨季差异不显著(P>0.05)，其变化范围为8.53～8.76，呈弱碱性. WT在雨季(21.9℃±2.2℃)显著高于旱季(12.4℃±1.2℃)(P<0.05)，其区域性差异不显著(P>0.05). 旱季陆源输入负荷较低[16], WT较低，浮游植物生物量较低，水生植物占优势，其通过光合作用消耗了CO2，促使水体pH和DO浓度升高，为氧化环境；而雨季陆源输入负荷增加，加之WT升高，藻类快速繁殖，藻类等浮游植物占优势，透明度下降[28]，水生植物退化，水生和浮游植物残体及陆源有机质分解，消耗DO，释放CO2，促使水体DO浓度和pH值下降，为还原环境.

2.2 洱海悬浮颗粒物有机碳与氮同位素特征及来源

由图2可知，洱海悬浮颗粒物POM浓度在旱、雨季差异显著(P<0.05)，雨季(8.6 mg/L±1.5 mg/L)是旱季(2.5 mg/L±0.6 mg/L)的3.4倍；空间分布特征为南部(9.9 mg/L)>中部(8.4 mg/L)>北部(7.6 mg/L). POM浓度的时空特征与Chl.a浓度一致，说明其与Chl.a具有同源性. 悬浮颗粒物的δ13C在旱、雨季差异显著(P<0.05)，其变化范围分别为-31.75‰～-18.21‰(均值为-25.34‰±4.14‰)和-23.8‰～-14.7‰(均值为-20.2‰±3.3‰). 雨季δ13C空间异质性显著(P<0.05)，偏负程度为北部(-22.8‰±3.14‰)>中部(-20.9‰±2.14‰)>南部(-16.9‰±4.14‰). 如图3所示，根据悬浮颗粒物δ13C值变化范围初步判定有机碳来源于陆源植物(如C3、C4植物)、土壤有机质、水生植物和浮游植物等[8,29]. 旱季以陆源C3植物、土壤有机质和水生植物混合为主，雨季以水生植物、浮游植物和陆源C4植物混合为主.

由图2可知，δ15N在旱、雨季差异显著(P<0.05)，变化范围分别为4.9‰～7.4‰(均值为6.4‰±1.3‰)、7.4‰～10.8‰(均值为9.3‰±1.8‰). δ15N结果显示水生植物(5.53%±5.12‰)和浮游植物(7.36%±3.89%)并不存在显著差异(P>0.05)，因此，水生植物和浮游植物合并为一个氮来源端元. 根据δ15N值初步判断其来源主要为土壤流失氮、大气干湿沉降、工业及生活污水、化肥、水生植物和浮游植物等混合输入[30]. 有机质C/N比值较高代表陆源植物来源，一般大于15[31]；而比值较低代表内源菌藻类或者浮游生物有机残体，介于3～10之间[11]；且随C/N比值降低分解程度增加[32]. 由图3所示，旱季悬浮颗粒物有机质C/N比值(9.1～16.9(均值为13.3±2.7))显著高于雨季(4.6～8.9(均值为7.1±1.6))，表明悬浮颗粒物有机质来源在雨季以内源菌藻降解产物为主，易被降解；在旱季陆源贡献相对增加，不易降解.

图2 洱海悬浮颗粒物δ13C、δ15N、C/N和POM特征值的时空分布特征

图3 洱海悬浮颗粒物δ13C(a)、δ15N(b)和C/N(c)的来源分布特征[12,33]

2.3 洱海表层沉积物有机碳与氮同位素特征及来源

如图4、5所示，表层沉积物有机质含量w(TOM)在雨季(81.5 g/kg±26.2 g/kg)高于旱季(70.6 g/kg±19.4 g/kg). 全年w(TOM)的空间异质性表现为：北部(85.7 g/kg±31.2 g/kg)>南部(77.1 g/kg±14.9 g/kg)>中部(53.3 g/kg±12.7 g/kg). 可能由于北三江为洱海的主要水源，携带大量污染物在北部沉积，且北部湖湾水草茂盛，近年来水生植物退化腐烂导致大量内源有机碎屑沉积，使洱海北部w(TOM)较高. 中部区域水深较深，且由于环流作用，有机质以小颗粒沉积为主，且随着深度增加，DO浓度较低，有机质在厌氧微生物作用下不断矿化分解，从而使w(TOM)在中部呈下降趋势. 洱海表层沉积物δ13C在旱、雨季差异不显著(P>0.05)，变化范围分别为-25.1‰～-8.6‰(均值为-18.5‰±5.5‰)和-24.3‰～-13.0‰(均值为-18.9‰±4.2‰)；δ13C空间异质性显著(P<0.05)，δ13C偏负程度为中部(-22.9‰±1.5‰)>南部(-17.9‰±1.3‰)>北部(-14.8‰±0.3‰). 因此，初步判断表层沉积物有机碳来源主要受陆源C4、C3植物、土壤有机质和水生植物影响. 北部以陆源C4植物和水生植物混合为主，中部以陆源C3植物、土壤有机质和水生植物混合为主，南部以陆源C4植物为主. 洱海表层沉积物δ15N在旱、雨季差异显著(P<0.05)，变化范围分别为1.9‰～4.9‰(均值为3.6‰±1.5‰)和0.7‰～7.8‰(均值为4.2‰±1.8‰)；北、中、南部分别为3.8‰±0.1‰、4.5‰±0.5‰、4.2‰±0.2‰. 相比3大富营养化湖泊太湖(4.5‰～15.2‰)、巢湖(3.3‰～10.4‰)和滇池(5.2‰～11.0‰)的沉积物δ15N范围[30,34-35]，洱海表层沉积物δ15N变化范围较小，说明洱海沉积物受外源氮污染负荷影响相对较小. 根据δ15N值范围初步判断洱海表层沉积物氮受土壤流失、化肥、陆源植物、工业及生活污水、水生植物和浮游植物等来源共同影响. 洱海表层沉积物有机质C/N在旱、雨季无显著差异(P>0.05)，变化范围分别为8.5～25.6(均值为12.8±5.3)和8.5～17.1(均值为11.4±3.2)；区域性差异显著(P<0.05)，表现为北部(15.0±0.8)>南部(12.1±0.4)>中部(9.6±0.7). 先前研究认为沉积物有机质C/N>8，即被认为是受2种物源的影响，而陆源有机质所占的比例越高，C/N值越大，说明洱海表层沉积物有机质来源有一定的复杂性[20].

图4 洱海表层沉积物δ13C、δ15N、C/N及TOM的时空分布特征

图5 洱海表层沉积物δ13C(a)、δ15N(b)和C/N(c)的来源分布特征

2.4 洱海悬浮颗粒物和表层沉积物有机碳与氮来源特征差异

利用δ13C、δ15N和C/N 3种不同组合综合分析悬浮颗粒物和表层沉积物有机质来源(图6)，不断缩小特征范围，降低同位素重合的部分，从而精确悬浮颗粒物和表层沉积物有机质来源特征[12,35]. 结果表明，洱海悬浮颗粒物有机质旱季以水生植物和土壤有机质混合来源为主；雨季以浮游植物来源为主. 表层沉积物来源复杂，在旱、雨季均以水生植物、土壤有机质和陆源植物混合来源为主. 因此，悬浮颗粒物与表层沉积物之间联系并不紧密, 雨季藻类快速繁殖对悬浮颗粒物有机质贡献较大，其沉降过程易降解，导致表层沉积物有机质浮游植物贡献较低，以物理侵蚀产物或长期沉积的悬浮物为主.

图6 洱海旱季(a、b、c)和雨季(d、e、f)悬浮颗粒物和表层沉积物δ13C、δ15N和C/N之间的关系(A：浮游植物；B：水生植物 ；C：土壤有机质 ；D：陆源植物)

利用Iso Source模型分别计算悬浮颗粒物和表层沉积物有机碳、氮来源端元贡献率(图7). 结果表明：悬浮颗粒物有机碳来源贡献在旱季为陆源C3植物(38.2%～51.2%)>水生植物(13.0%～32.9%)>陆源C4植物(2.5%～20.4%)>土壤有机质(8.8%～13.4%)>浮游植物(4.2%～12.8%)；在雨季为浮游植物(38%～50%)>土壤有机质(24%～27%)>陆源C3(4%～19%)>C4植物(7%～17%)>水生植物(12%～14%). 悬浮颗粒物氮来源贡献在旱季以陆源植物(34%～47%)和工业及生活污水(18%～22%)为主，在雨季以水生和浮游植物(32%～45%)以及化肥(18%～26%)为主. 说明从旱季到雨季，藻类等浮游植物来源的有机碳和氮贡献增大；悬浮颗粒物氮来源由外源向内源、由点源向面源转变. 表层沉积物有机碳在旱、雨季来源均以陆源C4植物为主(48.2%±19.1%)，氮来源旱季以陆源植物为主(44.3%±10.1%)转变为雨季以化肥为主(30.3%±6.8%). 先前研究表明[36]，2013-2014年洱海流域种植作物主要为C3植物，占总种植面积的50%以上，包括稻谷、大麦、蚕豆、蔬菜等. C4植物主要为玉米，种植距离洱海水体较近，约占总种植面积的20%，集中于洱海流域北部(上关镇、挖色镇、双廊镇、右所镇)和南部(凤仪镇)(图1)，因此，洱海北部和南部湖区表层沉积物C4植物来源的有机碳较高. 北部水生植物生物量和覆盖度较大，导致表层沉积物北部水生植物来源有机碳(23%～32%)高于中部和南部(4%～19%)；南部藻类生物量(26.9 μg/L±2.7 μg/L)高于北部(19.4 μg/L±3.3 μg/L)和中部(22.6 μg/L±4.7 μg/L)(表1)，导致南部悬浮颗粒物和表层沉积物有机碳浮游植物来源贡献(14%～50%)高于北部(1%～41%)和中部(5%～38%). 表层沉积物氮来源在旱季以陆源植物为主(33%～52%)，在雨季以化肥为主(25%～38%)，特别是在北部(38%)高于中部和南部(25%～28%)，证实雨季洱海北部受流域农业面源污染影响增大，北部流域农田密布，长期不合理施用化肥，导致大量农田氮素流失[37].

图7 不同端元物质对样品有机碳与氮的贡献率

悬浮颗粒物和表层沉积物有机碳、氮来源特征差异性(图8)表现为：悬浮颗粒物有机碳旱季以陆源C3植物(46%)和雨季以浮游植物(44%)为主转变为表层沉积物有机碳在旱、雨季均以陆源C4植物(48%～58%)和土壤有机质(17%～26%)为主，表明悬浮颗粒物在沉降过程中浮游植物和陆源C3植物来源有机碳易降解，参与水体营养盐循环，而陆源C4植物和土壤有机质来源有机碳易沉积，对水质影响相对较小. 雨季悬浮颗粒物氮以水生和浮游植物来源(40%)为主转变为表层沉积物氮以化肥(30%)、土壤流失(22%)和陆源植物(21%)来源为主，推断悬浮颗粒物在沉降过程中水生和浮游植物来源氮易降解，而陆源植物、土壤流失和化肥氮易沉降累积. 此外，C/N是能量-物质传输效率的一个关键指标[38].悬浮颗粒物沉降过程中有机氮优先矿化降解导致表层沉积物C/N(11.4±3.2)高于悬浮颗粒物C/N(7.1±1.6). 表层沉积物C/N升高可能降低跨营养水平能量传递效率[39]，在一定程度上控制沉积物营养盐释放，对水质影响较小. 表层沉积物δ15N值明显低于悬浮颗粒物也表明15N消耗是悬浮颗粒物沉降过程中N循环过程的重要组成部分[40-41]. 先前研究也证实洱海表层沉积物有机质矿化程度低，以结构复杂的腐殖质为主[42].

悬浮颗粒物与表层沉积物δ13C、δ15N、C/N没有显著相关性(表2)，进一步证实表层沉积物来源于降解转化后长期沉积的悬浮物或者物理侵蚀的产物[20]，局部环境的改造作用也使悬浮颗粒物和表层沉积物有机质δ13C、δ15N、C/N产生差异[39]. 因此，悬浮颗粒物和表层沉积物有机碳和氮来源特征差异性揭示不同来源端元有机碳和氮迁移的转化和降解特征(图8)，记录和反映其生物地球化学循环.

图8 悬浮颗粒物和表层沉积物有机碳与氮来源贡献差异分析

2.5 洱海悬浮颗粒物和表层沉积物有机碳与氮来源特征对水质的指示意义

表2 洱海悬浮颗粒物和表层沉积物有机碳、氮来源特征与水质的相关关系

随机森林预测变量重要性估计表明(图9)，影响水体多营养循环指标的主导因子是悬浮颗粒物δ13C和C/N，重要度分别为6.9%和5.7%，其次是水环境因子中ORP和WT，重要度分别为4.3%和4.0%；影响Chl.a浓度最重要的指标是WT，重要度为10.7%，其次是POM和C/N，重要度分别为9.5%和9.3%；表层沉积物有机质及其碳、氮同位素对水体多营养循环指标和Chl.a影响程度最小，说明POM来源特征结合环境因子显著影响水质，而表层沉积物有机质来源对水质无显著影响. 洱海POM来源变化特征与水质波动具有良好的响应关系，其水质的指示意义为：从旱季到雨季，洱海POM浓度增加，且有机碳与氮来源由以陆源植物为主转变为以浮游植物为主，δ13C、δ15N均呈现增大趋势，这是因为藻类大量繁殖，当水体中溶解性CO2和硝酸盐无法满足藻类需求时，浮游植物开始富集13C、15N，随着水体中碳酸盐和硝酸盐的减少，有机质δ13C、δ15N逐渐偏正[47]. 悬浮颗粒物与表层沉积物有机碳与氮来源特征差异表明浮游植物来源的POM在沉降过程易降解释放溶解态营养盐，是重要的营养盐源，导致水体ρ(DON)、ρ(TOC)、ρ(TN)和ρ(TP) 增大(表1)，ρ(SRP)和ρ(NH3-N)降低是由于SRP和NH3-N被藻类大量吸收利用，从而影响水体营养盐循环；由于POM矿化作用消耗DO，ORP值降低，处于还原性环境(表1)，利于厌氧菌的繁殖[27]，导致水体恶性循环. 先前研究证实POM中有机氮通过矿化作用释放氨氮，同时O2消耗刺激碱性磷酸酶活性促进磷释放，从而驱动藻华暴发[48]. 有机氮优先矿化导致C/N比值增加，所以POM矿化分解影响Chl.a浓度[4]. 由此可见，湖泊浮游植物来源POM增加驱动了水体营养盐循环，造成严重的水质问题，对藻华暴发具有重要指示意义，需引起关注. 洱海表层沉积物主要来源于长期沉积降解转化后的悬浮物、难降解的陆源植物或者土壤侵蚀物，有机质结构较稳定，对水质影响相对较小. 先前研究也证实洱海沉积物有机质比长江中下游湖泊有机质含量高，但其结构稳定，不易降解，生物有效性较低[49-50].

图9 基于随机森林模型的多营养盐循环指标(a)和叶绿素a(b)的主要预测因素(*表示0.05水平下变量对模型的贡献显著. 悬表示悬浮颗粒物，沉表示表层沉积物)

根据洱海有机碳和氮来源端元贡献的时空差异，对洱海保护提出针对性管理措施. 在加强控制外源负荷的同时(如提高洱海流域植被覆盖率，减少土壤流失；重点削减北部农业面源及周边村落生活污染源等)，雨季还需重点采取除藻措施，不仅有效地控制藻生物量，而且对阻断藻-水质内源恶性循环，防止“二次水华”暴发也具有积极作用.

3 结论

洱海悬浮颗粒物和表层沉积物有机碳与氮同位素时空变化特征提供了其来源信息，研究结果为洱海水华期水质变化的解析提供了新的视角，并为提出针对性管理措施提供了理论依据. 主要研究结论如下：

1)悬浮颗粒物和表层沉积物有机碳与氮来源差异明显. 悬浮颗粒物有机碳来源旱季以陆源C3植物为主(46.0%±6.9%)转变为雨季以浮游植物为主(43.3%±6.1%)，氮来源旱季以陆源植物为主(40.7%±6.5%)转变为雨季以水生和浮游植物为主(39.9%±6.6%)；表层沉积物有机碳在旱、雨季来源均以陆源C4植物为主(48.2%±19.1%)，氮来源旱季以陆源植物为主(44.3%±10.1%)转变为雨季以化肥为主(30.3%±6.8%). 两者有机碳与氮来源差异揭示水生和浮游植物来源的有机碳、氮易降解，对水质影响较大，而陆源C4植物和土壤来源的有机碳和化肥氮易沉积，对水质影响相对较小.

2)POM来源差异(δ13C、C/N)、ORP和WT是影响水体多营养循环指标的重要因子，POM、C/N和WT是影响Chl.a浓度的重要因子，说明POM来源特征结合水环境因子显著影响水质，对水华暴发具有重要指示意义，需引起关注. 而表层沉积物有机碳与氮来源特征对水质影响不显著. 从旱季到雨季，WT升高，浮游植物快速生长，浮游植物来源POM增加，δ13C偏负程度降低，C/N值下降，其易降解，成为重要的营养盐源，进而影响水体营养盐循环，促进藻类“二次”暴发；由于POM矿化降解耗氧，改变水体ORP，造成水质恶性循环.

3)对于洱海保护，根据洱海有机碳与氮来源贡献的时空差异，提出针对性的管理措施. 除加强外源负荷控制，雨季还应重点控制浮游植物生物量，避免“二次”水华发生.