环境科学  2019, Vol. 40 Issue (7): 3039-3048   PDF    
基于碳氮同位素的澜沧江水库TOC来源差异性分析
王从锋1,2, 胡子龙1, 杨正健1, 丹勇1, 赵萍1, 郝文超1, 聂小芬1, 聂睿1, 徐刘得1, 刘德富1,3     
1. 三峡大学水利与环境学院, 宜昌 443002;
2. 三峡库区生态环境教育部工程研究中心, 宜昌 443002;
3. 湖北工业大学土木建筑与环境学院, 河湖生态修复与藻类利用湖北省重点实验室, 武汉 430068
摘要: 水库建设所导致的生态问题正引发前所未有的深度思考,梯级水库在时间和空间尺度上的影响效应更加值得探究.为揭示新老水库有机碳分布特征及其沉积物TOC来源的差异,于2017年11月采集了苗尾、功果桥和大朝山这3个不同时期建设的水库的水样和沉积物样柱.测定了水体的温度(T)、溶解氧(DO)、总有机碳含量(TOC)、沉积物总有机碳(TOC)、氧化还原电位(ORP)、总氮(TN)和总磷(TP)等指标,并利用15 N和13C同位素,结合IsoSource软件,解析沉积物中TOC来源及其来源物质对相应水库沉积物中TOC的贡献量,从而探究其内在的碳循环机制和梯级水库演进模式.结果表明,苗尾、功果桥和大朝山水库水体有机碳质量浓度平均值分别为0.95、1.97和4.64 mg·L-1.对应水库沉积物中有机碳含量变化范围分别为4.41~81.63、18.30~28.42和9.16~14.46 g·kg-1.水库的梯级建设,使得新老水库的沉积物来源,周围补给等出现差异,使得新老水库TOC出现巨大差异.对于水体TOC,水体热力学状态和溶解氧的差异从而间接影响了水体中TOC的分布趋势.沉积物中主要考虑生源要素的影响,即沉积环境对有机质的保存能力是造成DCS、MV和GGQ沉积物垂向分布有差异的主要原因.而梯级水库的演进模式中,MV时空尺度上处于第一级段,是以累积上游来流TOC为主;GGQ处于第二级段,是以主要消耗分解上游来流TOC;DCS处于第三级段,是以主要积累水库周边TOC来源.
关键词: 沉积物      有机碳      碳氮同位素      差异分析      新老水库     
Analysis of Total Organic Carbon Source Differences Between New and Old Cascade Reservoirs using Carbon and Nitrogen Isotopes
WANG Cong-feng1,2 , HU Zi-long1 , YANG Zheng-jian1 , DAN Yong1 , ZHAO Ping1 , HAO Wen-chao1 , NIE Xiao-fen1 , NIE Rui1 , XU Liu-de1 , LIU De-fu1,3     
1. College of Hydraulic and Environmental Engineering, China Three Gorges University, Yichang 443002, China;
2. Engineering Research Center of Eco-environment in Three Gorges Reservoir Area Region, Yichang 443002, China;
3. Hubei Key Laboratory of Ecological Restoration of River-lakes and Algal Utilization, College of Civil Architecture and Environment, Hubei University of Technology, Wuhan 430068, China
Abstract: The ecological problems due to reservoir construction are causing unprecedented concern. To reveal the differences in organic carbon distribution characteristics and sediment sources of total organic carbon (TOC) between the old and new reservoirs, water samples, and sediment samples from reservoirs constructed in the three different periods of Miaowei, Gongguoqiao, and Dachaoshan were collected in November 2017. The temperature (T), dissolved oxygen (DO), TOC, redox potential (ORP), total nitrogen (TN), and total phosphorus (TP) of the water samples were measured. The isotopes 15N and 13C were used as indicators with IsoSource software to analyze the contributions of TOC sources and their source materials to the corresponding reservoir sediments, in order to explore the carbon cycle mechanism and evolution mode of reservoir. The results showed that the average concentrations of organic carbon in the waters of the Miaowei, Gongguoqiao, and Dachaoshan Reservoirs were 0.95 mg·L-1, 1.97 mg·L-1, and 4.64 mg·L-1, respectively. The range of organic carbon content in the corresponding sediments was 4.41-81.63 g·kg-1, 18.30-28.42 g·kg-1, and 9.16-14.46 g·kg-1, respectively. The cascade construction of the reservoirs resulted in a difference between the sediment sources of the new and old reservoirs and the surrounding recharge area, meaning that the TOC of the new and old reservoirs were significantly different. For the TOC of waterbodies, the difference between the thermodynamic state of water and dissolved oxygen indirectly affects the distribution trend of TOC. The sediments mainly reflect the influence of source elements, that is, the ability of the sedimentary environment to preserve organic matter was the main cause of the vertical distribution of DCS, MV, and GGQ sediments. In the evolution mode of cascade reservoir, the research shows that it can be preliminarily set as three stages. Firstly, due to the short age of MV, it is in the first stage and mainly accumulates the TOC from the upstream. GGQ is longer than the age of MV, and it is mainly used to decompose the upstream TOC, so it is defined as in the second stage. Finally, as an old reservoir, DCS mainly accumulates TOC sources around the reservoir, which can be regarded as the third stage.
Key words: sediment      organic carbon      carbon and nitrogen isotopes      difference analysis      new and old reservoir     

河流具备的陆-海二元纽带属性, 决定其本身也是陆地系统重要组成部分.传统的河流连续体[1]、洪水脉动[2]和河流生态系统四维结构[3]等生态理念, 如何与时俱进, 去解决社会、经济、人文和环境等多元矛盾的冲突.人类不得不面临这一选择困境.迫于对绿色能源的需求, 世界各国仍在大量修建大坝等水工建筑物[4].而水库“源-汇”效应在时-空尺度上都广泛影响世界河流的自然过程[5~8].目前“河流水库化”成为陆地水文系统的重要特点[9].碳在河流元素属性中不容忽视, 是底泥和水体重要组成部分.水库内部几乎所有生物、物理和化学变化都与其息息相关[10, 11].反映河流碳的生源要素指标中, 总有机碳(TOC)极其重要, 它能够反映水体中有机物污染程度[12].并且水体中有机碳是生物圈最大的活性有机碳库, 在生物、地质和化学过程发挥着主要作用[13].并且, 在检测方面, TOC测定方法比COD和BOD的测定更方便快捷和准确, 能全面反映水体含碳有机物污染状况, 因此对TOC的研究也越来越引起人们的重视[14, 15].

国内外对于水体中TOC的研究主要集中在海、河与浅水湖泊, 对于人造水库则鲜有报道[16].但人类活动修筑的大型水库或拦河坝, 已经影响到河流碳的生物地球化学过程[13, 17].其对生态环境的胁迫, 使得科研人员对于水库的TOC相关研究迫在眉睫.而目前针对与河流TOC相关研究表明其来源主要包括外源和内源[18].水库的建设, 使得河流这种自然模式受到影响.水库在时间层次上的更替需要很长的区间, 一般大型水库都是百年工程, 从而使得这种影响在时间尺度上极其深远[19].通过文献[20~26], 本文定义新水库库龄 < 10 a, 老水库库龄>10 a.正因为新老水库在时-空尺度的差异, 再加上水库本身不同的属性, 这些都使得梯级水库对河流的碳循环的影响程度和相关机制研究明显不足[20].

而本研究基于此出发点, 探究新老水库内在的相关碳循环机制.而在目前可利用的研究工具中, δ13C被广泛用作河流有机物质来源的示踪工具[21, 22], 并且常用C同位素的比值与标准同位素比值的千分偏差(δ13C)来表示稳定同位素13C, 如对密西西比河河口有机质的δ13C研究表明, 滨海沼泽和C3维管束植物是下游河水有机碳的主要来源[23].然而, 由于微生物降解对13C产生分馏作用, 影响有机质来源的指示作用[24].基于13C会产生分馏作用的不利影响, 因此本文综合运用δ15N、δ13C和TOC等, 去探究水库碳的相关机制.考虑到西南地区修建有大量的水库大坝.这些水库大部分属于特殊、狭长型河道梯级深水水库.水库的沉积物中碳滞留状况和来源属性, 以及对区域乃至全球碳循环的意义有多大等问题目前尚不清楚.由于有机碳与水体富营养化密切相关, 到目前为止, 苗尾(MV)、功果桥(GGQ)和大朝山(DCS)水库极少有这方面的报道, 且很少有学者和研究人员将梯级水库下的新老水库做综合研究, 因此本文综合前人的研究方法和理论[25, 26], 利用同位素技术及常规环境因子指标, 结合IsoSource软件, 分析其新老水库TOC的来源及分布, 以及相应的机制, 这不仅有助于研究澜沧江流域碳循环和有机质的贮藏量, 而且还有助于分析流域环境变化趋势, 以期为改善该流域生态环境提供参考依据.

1 材料与方法 1.1 研究区域与样点设置

澜沧江是湄公河上游在中国境内河段的名称, 是中国西南地区的大河之一.图 1所示, MV、GGQ和DCS水库均是澜沧江上的梯级水电站, MV水电站位于云南省大理白族自治州云龙县, 坝址距昆明544 km.地处横断山脉澜沧江纵谷地区, 是澜沧江干流水电基地中下游河段(古水至苗尾河段)一库七级开发方案中的最下游一个梯级电站, 该水电站于2016年左右蓄水, 并与该年建成投产, 距今只有短短2 a. GGQ位于云南省云龙县境内, 是澜沧江干流中下游“两库八级”梯级水电开发的最上游一级水电站, 该电站位于横断山纵谷区, 河谷深切, 于2011年左右蓄水, 距今7 a. DCS位于中国云南省云县与景东县交界的澜沧江中游, 那果河支流入口下游1.5 km处.在2001年左右蓄水, 于2003年建成投产, 距今17 a.通过本文定义, 可以初步认为DCS为老库, MV和GGQ为新库, 且MV库龄最小.

图 1 澜沧江采样点 Fig. 1 Sampling layout of the Lancang River

本研究于2017年11月, 分别对MV、GGQ和DCS进行采样.

1.2 样品采集与测试 1.2.1 样品采集与处理

本研究采用重力式柱状采泥器在样点中采集沉积物样品以及上覆水, 采泥器原理为依靠自身重力垂直采集沉积物样品, MV、GGQ与DCS沉积物按照2 cm一分层进行现场分样, 其他上游点位, 因沉积较少, 因此将所取沉积物混合后, 取表层, 并将沉积物用干净聚乙烯自封袋密封保存, 排尽袋中的空气, 上覆水立即装入300 mL聚乙烯采样瓶, 4密封保存后带回实验室[25].

根据澜沧江水文地质状况和流域分布情况, 选取了具有代表性的表层水采样地点, 使样品尽可能地覆盖整个流域.

1.2.2 同位素数据处理

使用预先处理过的WhatmanGF/FGF-F玻璃纤维滤膜(0.47 μm, 预先在马弗炉500℃下燃烧4 h以除去无机物, 直至滤膜恒质量)收集水体表层悬浮颗粒物, 具体参考相关人员方法[26].

沉积物则将已经分层的样品, 首先称取土样0.5~1.5 g装入离心管中, 用0.1 mol·L-1稀HCl反复洗涤去掉土样中的无机碳, 再用去离子水反复洗涤, 直到离心管中的溶液呈中性, 置于50℃下烘干的土样放入玛瑙研钵碾磨, 过筛200目, 利用精度为百万分之一电子天平称取少量碾磨过的土样, 连滤膜一起用锡箔纸包好并放入冰箱于-20℃下冷冻避光保存.包于锡杯中, 最后将包好的样品逐个放入仪器里进行测试.

1.2.3 环境因子数据处理

水温(T)和溶解氧(DO)这两个指标均用HydrolabDS5多参仪(美国)现场测定, 沉积物氧化还原电位(ORP)由中国科学院南京土壤研究所研发QX-6530便捷式氧化还原电位仪测定. TP用钼锑抗比色法, TN用凯氏定氮法测定.

1.2.4 沉积物与水体TOC分析

水体中TOC采用型号为德国耶拿MutiN/C3100的TOC分析仪, 实验时取10 mL样品, 加入0.5 mL 2 mol·L-1盐酸后上机检测, 主要测定原理是非分散性红外线技术(NDIR).底泥中的TOC采用型号为意大利EA300元素分析仪, 其利用的专利为flash燃烧技术, 实验时取10 g经真空冷冻干燥过的底泥样品上机检测.

1.3 数据分析

基于IsoSource软件的物源解析方法, 其基本原理是根据稳定同位素质量守恒, 以质量守恒混合模型为基础, 开发用于计算源头值>(n+1)(n为多种同位素值)潜在贡献比例.该模型在一定增量范围内, 通过测试的同位素值信息, 使用标准线性混合模型来模拟每一种可能污染比例, 其总值为100%, 从而确定混合物各部分的比例[26].

用该软件依据同位素各端元值.算出贡献率结果与沉积物含量依次数乘后, 即得沉积物中TOC的具体含量值.

(1)

式中, Q为该层次的TOC含量(g·kg-1), P为来源百分比(%), F为该层TOC含量(g·kg-1).

2 结果与分析 2.1 澜沧江新老水库水体-沉积物TOC分布规律

为证明新老水库TOC存在差异性, 阐明其具体表征.通过研究发现, 图 2可看出新老水库水体TOC质量浓度在垂向上呈现出先减小后增大的趋势(DCS除外); DCS和MV水体中TOC质量浓度最大值均出现在表层水体中, 而GGQ水体中TOC质量浓度最大值出现在上覆水体中.相对于MV和GGQ水库, DCS水体表层TOC质量浓度最大, 其值为4.94 mg·L-1.上覆水中, MV水库TOC含量最小, 其值为1.08 mg·L-1.总体上, 各层水体TOC质量浓度表现为:DCS水库水体TOC质量浓度>GGQ水库水体TOC质量浓度>MV水库水体TOC质量浓度.从箱线图可以看出, MV与DCS的TOC无明显摆动幅度, 而GGQ变化明显, 分布区间广.沉积物中TOC变化趋势如图 2所示, 在0~4 cm处MV水库TOC含量≫GGQ水库TOC含量>DCS水库TOC含量, 其最大值81.63g·kg-1, 在4~6 cm处, MV水库沉积物TOC含量急剧递减, 而GGQ和DCS沉积物TOC含量增加, 在6~10 cm处, GGQ水库沉积物TOC含量>DCS水库TOC含量>MV水库TOC含量, 且三水库TOC含量均呈现递减的趋势.在0~6 cm处, MV和DCS水库TOC含量呈现出先减后增, 而GGQ水库TOC含量呈现递增趋势; 在6~16 cm处, MV和GGQ水库TOC含量变化趋势相同, 变化趋势为先减后增, 而DCS在6~8 cm处TOC含量表现为递增, 8~12 cm处含量表现为递减.其中MV和GGQ水库沉积物TOC含量在4~6 cm处均达到最大值, 其含量分别为16.55g·kg-1与28.43g·kg-1.整体上, GGQ水库沉积物TOC含量较MV和DCS都大.且从箱线图可看出, MV沉积物变化剧烈, 而DCS, GGQ相对无明显变化.

图 2 TOC分布及其箱线图 Fig. 2 Total organic carbon distribution and the associated box plot

2.2 澜沧江流域有机质同位素溯源特性

为说明澜沧江新老水库TOC量的属性与其来源有关, 本研究试图通过溯源分析去解决这个问题.因此通过表 1分析, 可看出MV、GGQ和DCS沉积物同位素的分布趋势, MV、GGQ和DCS的N同位素分别处于5.28‰~6.24‰、4.61‰~7.60‰和2.70‰~4.52‰.而其C同位素MV处于-16.34‰~-24.60‰区间; GGQ处于-16.96‰~-26.48‰区间; DCS处于-23.65‰~-27.94‰区间.

表 1 沉积物同位素值1) Table 1 Isotopic values of sediments

图 3为澜沧江流域C和N同位素分布值, 其端元简化后, 将其划分为这5个端元, 可知周边土壤C的同位素区间在-21‰~-27‰, N的同位素在0‰~8‰; 污水C的同位素区间在-21‰~-52‰, N的同位素在1‰~5‰; 悬浮有机质C的同位素区间在-21‰~-32‰, N的同位素在3‰~19‰; 枯枝落叶C的同位素区间在-19‰~-37‰, N的同位素在0‰~7‰; 上游来流C的同位素区间在-3‰~-27‰, N的同位素在2‰~9‰.表 1可知MV沉积物C和N同位素范围分别在-16.35‰~-24.60‰和5.28‰~6.24‰; GGQ沉积物C, N同位素范围分别在-16.96‰~-26.48‰和4.61‰~7.60‰; DCS沉积物C和N同位素范围分别在-23.65‰~-27.94‰和2.70‰~4.52‰. GGQ、DCS和MV的沉积物都落在多端元交汇处, 且DCS的悬浮有机质也落入多端元交汇处, MV和GGQ的悬浮有机质只落入其本体区间内.

图 3 同位素端元 Fig. 3 Isotope end element diagram

图 4为各来源对新老水库相应TOC贡献量, MV水库各层沉积物TOC的相关贡献量, 数据证明在各有机质来源中, 沉积物受上游来流的影响较大, 特别是0~2 cm和2~4 cm处沉积物TOC含量分别达39.59g·kg-1和26.26g·kg-1, 岸边泥土、耕地有机质和悬浮有机质对沉积物TOC贡献量占总体的50%;排污污水对TOC贡献量较少, 各层均小于8%的贡献量.不考虑小范围波动, 总体垂向上各来源随沉积物深度的加深, 上游来流对TOC的贡献量逐渐减小.对GGQ水库各层沉积物TOC的贡献量, 从表层至沉积物16 cm处, 上游来流对TOC的贡献量随沉积物深度的增加, 从4.62g·kg-1增加到11.11 g·kg-1, 且均为总体所占TOC的20%.从表层至沉积物12 cm处, 上游来流和排污污水对DCS沉积物各层TOC的贡献量总体上要小于其上游的MV和GGQ水库, 但上游来流对TOC的贡献量随沉积物深度的增加而增加, 而排污污水对TOC的贡献量随深度增加呈现出先减小后增加的趋势.

图 4 TOC来源示意 Fig. 4 Schematic diagram of total organic carbon sources

2.3 澜沧江新老库水体-沉积物环境因子特征

环境因子是水库内在属性的外在表现, 通过环境因子分析, 可以找出新老水库的同-异性.因此, 通过图 5可知, 水体温度, MV<GGQ<DCS, 在该时段内所测水体温度, DCS无水温分层, 而GGQ和MV存在弱水温分层, GGQ水体温度在9.42~10.61℃区间上变化, MV水体温度在8.39~9.80℃区间上变化, DCS稳定在16.74℃.其水体溶解氧含量可知MV>GGQ>DCS, 其刚好和水体温度相反.且MV是先增加后降低, 其在9.70~10.07 mg·L-1的区间上变化, 而GGQ于0~15 m处逐渐增加到9.63 mg·L-1, 后在15~35 m处于稳定, 且基本保持在9.65 mg·L-1; DCS只逐渐从8.09减小至7.48 mg·L-1.关于ORP:MV在37~-53 mV的范围内变化, GGQ在-17~-176 mV的范围内变化, 且在5 cm处有小范围波动DCS在-44.3~-257.3 mV的范围内变化, 其变化速率MV>GGQ>DCS.沉积物间隙水中TN,DCS较MV与GGQ波动幅度明显,DCS区间在4.04 ~9.43 mg·L-1, 而MV和DCS区间分别在1.77~2.75 mg·L-1和2.55 ~5.21 mg·L-1.沉积物间隙水中TP也与TN属性同步,DCS也较MV与GGQ波动幅度明显,DCS区间在0.03 ~0.12 mg·L-1.对于浊度来说, MV的浊度随着深度的增加, 不断增加, 从2.09NTU到8.8NTU; GGQ的浊度随着深度的增加, 不断减少, 然后在35 m出现增大趋势, 其变化区间在11~5.91NTU; DCS的浊度无明显变化, 在54 m左右有小范围波动, 其在1~2NTU内变化.

图 5 环境因子指标 Fig. 5 Environmental factor indicators

3 讨论

关于水库碳循环如图 6所示, 受上游来流输入、大气碳库以及水库周围碳源影响.而水库内部有机碳和无机碳系统的相互作用, 在不同生物作用、化学作用与物理作用下, 使得梯级水库TOC来源在时空尺度出现巨大差异性.碳来源的多途径性包括: ①在水库蓄水后淹没的原区域的碳[27~30]; ②从上游输入到水库中的有机物[31]; ③水库系统中生长和死亡各种生物; ④消落带的植被[32]; ⑤流域内受侵蚀土壤所带来的碳[33].

图 6 水库碳循环示意 Fig. 6 Schematic of the reservoir carbon cycle

3.1 水库建设对水体-沉积物TOC垂向分布的影响机制

以往研究发现, 水体中TOC在藻类水华暴发期以及水华过后容易积累, 浮游植物是产生TOC的主要来源[32].浮游植物在生长期间会排泄出有机物质, 据研究, 这部分可占到光合作用产物的30%~40%, 而且光合作用产物中的很大一部分转化为TOC释放到水体中的过程是非常快的, 几乎与光合作用同时进行[34, 35].这意味着浮游植物初级生产力高, 产生的TOC也高. TOC含量在GGQ、MV和DCS表、中和底层水体差异明显, 其可能原因很多. DCS与其他2座水库的主要差异与浮游植物有关.而不同的水体温度和溶解氧等, 将影响水库浮游植物的垂直分布[36, 37].适当的水体温度能促进浮游植物的新陈代谢和生长, 温度对水体的影响, 表现在间接地影响表层水体浮游植物的悬浮, 以及浮游植物垂直分布和演替[38].水体温度, MV<GGQ<DCS; 而水体溶解氧MV>GGQ>DCS, 水体热力学状态的差异和溶解氧的差异从而间接影响了水体中TOC的分布趋势.相对于MV和GGQ, 不考虑沉积物的影响, 影响DCS水体TOC可能还是主要由于周围枯枝落叶等有机质来源丰富, 随雨水或者地表径流流入河中.被水中浮游植物和微生物利用分解, 从而影响了水体中的TOC分布. MV的水体理化因子与功果桥的无明显差异, 其主要影响因素可能是沉积物的保存状态差异, 和外部环境的影响.沉积物经过分解矿化等过程后, GGQ和MV底部都呈现出增大的趋势, 而外部环境有机质输入影响使得其在表层和中层, 有着截然不同的浊度分布. MV底层浊度较高, 说明沉积物再悬浮或矿化等分解作用的趋势强, 而受周围环境的影响小, 但受上游来流影响大, 水体的TOC来源基本是分解水库内部TOC和水库截留的上游来流有机质.而GGQ的浊度表层高, 说明水库周围的有机质输入, 其逐渐沉降, 从而使得水体TOC呈现逐渐减少的趋势, 最终使得水体中GGQ的TOC大于MV, 变化趋势也有所差别.水体浊度是一种光学效应, 它反映了光线在透过水层时受到的阻碍.从浊度的分布来看, DCS的浊度偏小, 且无明显变化, 说明在水体中悬浮有机质较少的条件下, 光照充足, 进而水体中的悬浮有机质大多被分解成可溶态的碳, 所以尽管水体的浊度低, 但是TOC仍然很高.这些都说明MV、GGQ和DCS在水体TOC分布上存在本质区别.

沉积物中TOC含量取决于沉积物中有机质的输入量和沉积环境对有机质的保存能力[39], 从TOC等指标的相关性分析中, 看出GGQ、DCS和MV有着明显的区别.有机质的输入量来源与水生植物、流域侵蚀带来的陆源植物碎屑和人类污染排放物[40], 而MV相对于DCS、GGQ其表层有很高的TOC, 是由于水库蓄水后淹没大量储存在植被和土壤中的碳及外源泥沙有机质尚未矿化造成的.而造成GGQ、DCS和MV在沉积物TOC存在明显差异与其两者不同的理化环境和成库状态有关, 致使其沉积物处于不同的有机物赋存状态, 水库淹没陆地生态系统后, 改变被淹没土壤的化学性质, 导致其与之前形成截然不同的TOC分布和来源, 即沉积环境对有机质的保存能力不同.一方面, 可能是在不同层次的氧化还原环境氛围中, 使得TOC出现差异; 另一方面, 也可能是生源要素的含量差异, 使得TOC的来源与分布有着显著的差异.在不同层次的氧化还原环境氛围中, 其TOC的来源与分布应该都有显著的差异.从Eh中可以看出MV表层处于弱还原, 其余处于还原状态环境, 而GGQ和DCS整体处于还原环境.可以发现其三者基本都处于还原环境, 因此主要考虑生源要素的影响.而TN和TP变化剧烈, 且含量上的差异, 整体上TN趋势上是MV<GGQ<DCS, 而TP是GGQ<MV<DCS.这种差异影响了相关底栖动物和微生物等的生物活性, 导致不同程度的有机质滞留与分解, 从而使这些物质沉积在沉积物中量的差异性.沉积环境对有机质的保存能力, 即生源要素量的差异, 造成GGQ、DCS和MV沉积物垂向分布有如此差异的主要原因.

3.2 澜沧江不同水库沉积物来源差异分析

图 7所示, 水库的梯级建设, 使得新老水库存在一定的拦截模式. DCS、GGQ和MV分别距今17、7和2 a. MV建库时间短, 但由于梯级水库拦截作用, 使得GGQ和DCS上游来流补充的TOC减少. MV库龄短, 因此只有表层0~4 cm TOC沉积的含量高, 在此区段内MV水库TOC含量GGQ水库TOC含量>DCS水库TOC含量, 其最大值为81.63 g·kg-1, 底部在7 cm之后沉积的TOC含量就相对于GGQ和DCS来说较小, 且各种来源都很小.这符合相关人员研究的成果[20], 其研究表明水库存在拦截效应.而GGQ上游来流在总量上占比小, 且随深度增加逐渐升高.可以推测MV建成后, 其拦截作用, 使得表层补给速率小于分解速率, 从而表层TOC通过矿化等作用被分解, 释放到水体当中, 而中下层滞留于水库底部. GGQ建库时间相对比DCS短, 在上游来流的TOC分解还不完全的前提下, 周围植被土壤, 岸边泥, 枯枝落叶等物质的综合影响, 使得GGQ的TOC来源不同于MV.结合水库对河流碳循环的有机碳输入过程影响的相关研究, 对Verde的河流-水库系统研究表明, 水库拦截了72%的颗粒有机碳, 水库出口处TOC含量明显高于入口处, 而水库输出的TOC中, 约1/3来自于水库[41].并且在中国西洪家渡水库以及猫跳河流域水库研究也表明, 河流经水库作用后水体无机碳体系也发生了显著变化[42].这些学者和研究人员的研究均表明了水库通过影响河流的物理因子、化学反应、生物过程、从而造成了河流碳循环的异常.

图 7 梯级水库拦截模式 Fig. 7 Cascade reservoir interception mode

本研究表明, 不同库龄水库TOC来源差异化现象明显.建库时间最长的DCS, 一方面上游来流对其有机碳的补充量不到10%, 占总量少.另一方面周边补充量大于80%, 说明DCS已经使上游来流TOC基本分解, 而其累积TOC是后期补充的周围碳. DCS的TOC含量很高的现象, 可以证明老水库在上游来流很小的前提下, 除水库初期沉积的TOC和浮游生物产生的有机物外, 这些周边补充的有机质分解(比如污水, 植被腐蚀质以及冲蚀岸边的土壤等)也会成为TOC含量升高的来源.研究表明, 在库龄60 a的水库中仍然可以找到完好无损的树枝、树干、树根和下层土壤中的腐殖质[43].同样, Lokka水库蓄水30 a后, 库底仍然进行着强烈的厌氧分解, 而浅水区对土壤的侵蚀也会为水库提供有机碳[44~46].这也就解释了为何DCS比GGQ和MV建库时间都长, 而沉积物中仍然有如此高含量的TOC[47].一方面, 空间尺度上, 随着梯级水库的不断建设, 使得上游来流的碳逐渐被拦截; 另一方面, 时间尺度上, 随着时间的推移, 其水库早期蓄积的有机质被分解消耗, 而后期补充的是周围的碳源, 且梯级水库建设对于DCS、GGQ和MV的影响程度上DCS>GGQ>MV.因此如何定义临界值, 在发挥河流的自然调节能力的同时, 又使其避免生态失衡, 这有待进一步研究.

4 结论

(1) 影响澜沧江梯级水库新老水库水体TOC的主要原因是由于水体热力学状态和溶解氧的差异.

(2) 沉积物中主要考虑生源要素的影响, 即沉积环境对有机质的保存能力是造成DCS、MV和GGQ沉积物垂向分布有差异的主要原因.

(3) 在新库阶段, 澜沧江水库TOC主要与上游来流补充有关; 在老库阶段, 其TOC主要与周围补充有关.

(4) MV在时空尺度上处于第一级段, 以累积上游来流TOC为主; GGQ处于第二级段, 主要消耗分解上游来流TOC; DCS处于第三级段, 主要积累水库周边TOC来源.

参考文献
[1] Vannote R L, Minshall G W, Cummins K W, et al. The river continuum concept[J]. Canadian Journal of Fisheries and Aquatic Sciences, 1980, 37(1): 130-137. DOI:10.1139/f80-017
[2] Junk W J, Bayley P B, Sparks R E. The flood pulse concept in river-floodplain systems[A]. In: Dodge D P (Ed.). Proceedings of the International Large River Symposium[C]. Ontario: Canadian Special Publications Fisheries and Aquatic Sciences, 1989.
[3] Ward J V. The four-dimensional nature of lotic ecosystems[J]. Journal of the North American Benthological Society, 1989, 8(1): 2-8. DOI:10.2307/1467397
[4] Hooper R P, Aulenbach B T, Kelly V J. The national stream quality accounting network:a flux-based approach to monitoring the water quality of large rivers[J]. Hydrological Processes, 2001, 15(7): 1089-1106. DOI:10.1002/(ISSN)1099-1085
[5] Meybeck M, Helmer R. The quality of rivers:from pristine stage to global pollution[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 1989, 75(4): 283-309. DOI:10.1016/0031-0182(89)90191-0
[6] Vörösmarty C J, Sharma K P, Fekete B M, et al. The storage and aging of continental runoff in large reservoir systems of the world[J]. Ambio, 1997, 26(4): 210-219.
[7] Kelly V J. Influence of reservoirs on solute transport:a regional-scale approach[J]. Hydrological Processes, 2001, 15(7): 1227-1249. DOI:10.1002/(ISSN)1099-1085
[8] Humborg C, Danielsson Å, Sjöberg B, et al. Nutrient land-sea fluxes in oligothrophic and pristine estuaries of the Gulf of Bothnia, Baltic Sea[J]. Estuarine, Coastal and Shelf Science, 2003, 56(3-4): 781-793. DOI:10.1016/S0272-7714(02)00290-1
[9] 王雨春, 朱俊, 马梅, 等. 西南峡谷型水库的季节性分层与水质的突发性恶化[J]. 湖泊科学, 2005, 17(1): 54-60.
Wang Y C, Zhu J, Ma M, et al. Thermal stratification and paroxysmal deterioration of water quality in a canyon-reservoir, Southwestern China[J]. Journal of Lake Sciences, 2005, 17(1): 54-60. DOI:10.3321/j.issn:1003-5427.2005.01.009
[10] Schlesinger W H, Melack J M. Transport of organic carbon in the world's rivers[J]. Tellus, 1981, 33(2): 172-187. DOI:10.3402/tellusa.v33i2.10706
[11] Ittekkot V. Global trends in the nature of organic matter in river suspensions[J]. Nature, 1988, 332(6163): 436-438. DOI:10.1038/332436a0
[12] Ni H G, Lu F H, Luo X L, et al. Riverine inputs of total organic carbon and suspended particulate matter from the Pearl river delta to the coastal ocean off south China[J]. Marine Pollution Bulletin, 2008, 56(6): 1150-1157. DOI:10.1016/j.marpolbul.2008.02.030
[13] 韩志伟, 刘丛强, 吴攀, 等. 大坝拦截对河流水溶解组分化学组成的影响分析——以夏季乌江渡水库为例[J]. 长江流域资源与环境, 2009, 18(4): 361-367.
Han Z W, Liu C Q, Wu P, et al. Damming effect on water dissolved chemistry components-a case study of Wujiangdu reservoir water in summer, southwest China[J]. Resources and Environment in the Yangtze Basin, 2009, 18(4): 361-367. DOI:10.3969/j.issn.1004-8227.2009.04.011
[14] Yu H Y, Rezaee R, Wang Z L, et al. A new method for TOC estimation in tight shale gas reservoirs[J]. International Journal of Coal Geology, 2017, 179: 269-277. DOI:10.1016/j.coal.2017.06.011
[15] Liu E F, Shen J, Hang E L, et al. A geochemical record of recent anthropogenic nutrient loading and enhanced productivity in lake Nansihu, China[J]. Journal of Paleolimnology, 2010, 44(1): 15-24. DOI:10.1007/s10933-009-9382-z
[16] 滕明德, 迟峰, 高庚申. 万峰湖总有机碳空间分布特征[J]. 中国环境监测, 2014, 30(2): 35-39.
Teng M D, Chi F, Gao G S. Spatial characteristics of total organic carbon in Wanfeng Reservoir[J]. Environmental Monitoring in China, 2014, 30(2): 35-39. DOI:10.3969/j.issn.1002-6002.2014.02.007
[17] Kim B, Choi K, Kim C, et al. Effects of the summer monsoon on the distribution and loading of organic carbon in a deep reservoir, lake Soyang, Korea[J]. Water Research, 2000, 34(14): 3495-3504. DOI:10.1016/S0043-1354(00)00104-4
[18] 张永领. 河流有机碳循环研究综述[J]. 河南理工大学学报(自然科学版), 2012, 31(3): 344-351.
Zhang Y L. The review of the research of the riverine organic carbon cycle[J]. Journal of Henan Polytechnic University (Natural Science), 2012, 31(3): 344-351. DOI:10.3969/j.issn.1673-9787.2012.03.020
[19] dos Santos M A, Matvienko B, Rosa L P, et al. Gross greenhouse gas emissions from Brazilian Hydro reservoirs[A]. In: Tremblay A, Varfalvy L, Roehm C, et al (Eds.). Greenhouse Gas Emissions-Fluxes and Processes: Hydroelectric Reservoirs and Natural Environments[C]. Berlin, Heidelberg: Springer, 2005.267-291.
[20] Shi W Q, Chen Q W, Yi Q T, et al. Carbon emission from cascade reservoirs:spatial heterogeneity and mechanisms[J]. Environmental Science & Technology, 2017, 51(21): 12175-12181.
[21] 赵海娟, 肖琼, 吴夏, 等. 漓江地表水体有机碳来源[J]. 环境科学, 2017, 38(8): 3200-3208.
Zhao H J, Xiao Q, Wu X, et al. Sources of organic carbon in the surface water of Lijiang river[J]. Environmental Science, 2017, 38(8): 3200-3208.
[22] Stepanauskas R, Moran M A, Bergamaschi B A, et al. Sources, bioavailability, and photoreactivity of dissolved organic carbon in the Sacramento-San Joaquin river delta[J]. Biogeochemistry, 2005, 74(2): 131-149. DOI:10.1007/s10533-004-3361-2
[23] Wang X C, Chen R F, Gardner G B. Sources and transport of dissolved and particulate organic carbon in the Mississippi river estuary and adjacent coastal waters of the northern Gulf of Mexico[J]. Marine Chemistry, 2004, 89(1-4): 241-256. DOI:10.1016/j.marchem.2004.02.014
[24] Waterson E J, Canuel E A. Sources of sedimentary organic matter in the Mississippi river and adjacent Gulf of Mexico as revealed by lipid biomarker and δ13CTOC analyses[J]. Organic Geochemistry, 2008, 39(4): 422-439. DOI:10.1016/j.orggeochem.2008.01.011
[25] 望雪, 程豹, 杨正健, 等. 澜沧江流域沉积物间隙水-上覆水营养盐特征与交换通量分析[J]. 环境科学, 2018, 39(5): 2126-2134.
Wang X, Cheng B, Yang Z J, et al. Differences in diffusive fluxes of nutrients from sediment between the natural river areas and reservoirs in the Lancang river basin[J]. Environmental Science, 2018, 39(5): 2126-2134.
[26] 唐咏春, 徐飘, 杨正健, 等. 澜沧江流域水体悬浮颗粒物δ15 N空间差异及成因分析[J]. 环境科学, 2018, 39(11): 4964-4971.
Tang Y C, Xu P, Yang Z J, et al. Spatial difference and causes analysis of the δ15 N of suspended particulate matter in the Lancang river basin[J]. Environmental Science, 2018, 39(11): 4964-4971.
[27] Delmas R, Galy-Lacaux C, Richard S. Emissions of greenhouse gases from the tropical hydroelectric reservoir of Petit Saut (French Guiana) compared with emissions from thermal alternatives[J]. Global Biogeochemical Cycles, 2001, 15(4): 993-1003. DOI:10.1029/2000GB001330
[28] Rosa L P, Aurelio Dos Santos M, Matvienko B, et al. Biogenic gas production from major Amazon reservoirs, Brazil[J]. Hydrological Processes, 2003, 17(7): 1443-1450. DOI:10.1002/hyp.v17:7
[29] Kelly C A, Rudd J W M, Bodaly R A, et al. Increases in fluxes of greenhouse gases and methyl mercury following flooding of an experimental reservoir[J]. Environmental Science & Technology, 1997, 31(5): 1334-1344.
[30] Galy-Lacaux C, Delmas R, Kouadio G, et al. Long-term greenhouse gas emissions from hydroelectric reservoirs in tropical forest regions[J]. Global Biogeochemical Cycles, 1999, 13(2): 503-517. DOI:10.1029/1998GB900015
[31] Abril G, Guérin F, Richard S, et al. Carbon dioxide and methane emissions and the carbon budget of a 10-year old tropical reservoir (Petit Saut, French Guiana)[J]. Global Biogeochemical Cycles, 2005, 19(4): GB4007.
[32] Minor E C, Simjouw J P, Mulholland M R. Seasonal variations in dissolved organic carbon concentrations and characteristics in a shallow coastal bay[J]. Marine Chemistry, 2006, 101(3-4): 166-179. DOI:10.1016/j.marchem.2006.02.005
[33] Fearnside P M. Do hydroelectric dams mitigate global warming? The case of Brazil's Curuá-una dam[J]. Mitigation and Adaptation Strategies for Global Change, 2005, 10(4): 675-691. DOI:10.1007/s11027-005-7303-7
[34] 彭兴跃, 洪华生, 商少凌. 台湾海峡真光层有机碳动力学研究——DOC大幅度的日变化[J]. 海洋学报, 1997, 19(3): 57-65. DOI:10.3321/j.issn:0253-4193.1997.03.001
[35] Wood A M, Rai H, Garrnier J, et al. Practical approaches to algal excretion[J]. Marine Microbial Food Webs, 1992, 61: 21-38.
[36] Ahn C Y, Chung A S, Oh H M. Rainfall, phycocyanin, and N:P ratios related to cyanobacterial blooms in a Korean large reservoir[J]. Hydrobiologia, 2002, 474(1-3): 117-124.
[37] Han B P, Armengol J, Garcia J C, et al. The thermal structure of Sau reservoir (NE:Spain):a simulation approach[J]. Ecological Modelling, 2000, 125(2-3): 109-122. DOI:10.1016/S0304-3800(99)00176-3
[38] 李秋华, 商立海, 李广辉, 等. 万峰湖浮游植物群落的时空分布[J]. 生态学杂志, 2011, 30(5): 1031-1038.
Li Q H, Shang L H, Li G H, et al. Temporal and spatial characteristics of phytoplankton community in Wanfeng Reservoir[J]. Journal of Ecology, 2011, 30(5): 1031-1038.
[39] 马龙, 吴敬禄. 安固里淖湖积物中总有机碳含量及其碳同位素的环境意义[J]. 自然资源学报, 2009, 24(6): 1099-1104.
Ma L, Wu J L. Environmental significance from organic carbon and its isotope of Angulinao lake sediment[J]. Journal of Natural Resources, 2009, 24(6): 1099-1104. DOI:10.3321/j.issn:1000-3037.2009.06.017
[40] Kendall C, Silva S R, Kelly V J. Carbon and nitrogen isotopic compositions of particulate organic matter in four large river systems across the United States[J]. Hydrological Processes, 2001, 15(7): 1301-1346. DOI:10.1002/(ISSN)1099-1085
[41] Parks S J, Baker L A. Sources and transport of organic carbon in an Arizona river-reservoir system[J]. Water Research, 1997, 31(7): 1751-1759. DOI:10.1016/S0043-1354(96)00404-6
[42] 李干蓉, 刘丛强, 陈椽, 等. 猫跳河流域梯级水库夏-秋季节溶解无机碳(DIC)含量及其同位素组成的分布特征[J]. 环境科学, 2009, 30(10): 2891-2897.
Li G R, Liu C Q, Chen C, et al. Dissolve inorganic carbon and its carbon isotope composition in cascade reservoir of the Maotiao river during summer and autumn[J]. Environmental Science, 2009, 30(10): 2891-2897. DOI:10.3321/j.issn:0250-3301.2009.10.013
[43] Van Coillie R, Visser S A, Campbell P G C, et al. Evaluation de la degradation du bois de conifères immergés durant plus d'un demi-siècle dans un reservoir[J]. Annales de Limnologie, 1983, 19(2): 129-134. DOI:10.1051/limn/1983014
[44] Huttunen J T, Väisänen Hellsten S K, et al. Fluxes of CH4, CO2, and N2 O in hydroelectric reservoirs Lokka and Porttipahta in the northern boreal zone in Finland[J]. Global Biogeochemical Cycles, 2002, 16(1): 1003.
[45] 李思亮, 刘丛强, 丁虎, 等. δ13C-DIC在河流风化和碳生物地球化学过程中的应用[J]. 地球环境学报, 2012, 3(4): 929-935.
Li S L, Liu C Q, Ding H, et al. The application of δ13C-DIC on weathering and biogeochemical processes of carbon in rivers[J]. Journal of Earth Environment, 2012, 3(4): 929-935.
[46] 杨玉雪, 向鹏, 卢玮琦, 等. 贵州乌江渡水库沉积速率及碳氮埋藏通量估算[J]. 地球与环境, 2017, 45(1): 66-73.
Yang Y X, Xiang P, Lu W Q, et al. The sedimentation rate and burial fluxes of carbon and nitrogen in Wujiangdu reservoir, Guizhou, China[J]. Earth and Environment, 2017, 45(1): 66-73.
[47] 李昶.不同水库淤积形态对总有机碳、总氮埋藏通量的影响[D].上海: 华东师范大学, 2018.
Li C. The impacts of different types of reservoir siltation on the burial fluxes of total organic carbon and total nitrogen[D]. Shanghai: East China Normal University, 2018.