摘要
选择我国长江下游平原浅水湖泊升金湖86 cm沉积岩芯为研究对象,在可靠深度-年代模型基础上,通过沉积物质量累积速率(SAR)和总有机碳(TOC)含量计算升金湖过去千年有机碳埋藏速率变化;利用C/N与正构烷烃指标指示升金湖有机碳来源变化,并结合升金湖流域气候环境与人类活动历史记录探究该湖近千年有机碳埋藏的影响因素。研究结果表明:升金湖过去千年有机碳埋藏速率呈增加趋势,平均有机碳埋藏速率约为24.77 g/(m2·a),总有机碳储量约为3.33 Tg (1 Tg=1012 g)。根据升金湖有机碳埋藏速率变化特征,本研究将其划分为4个阶段:阶段Ⅰ:1000—1350 A.D.,有机碳埋藏速率较低但波动较大;阶段Ⅱ:1350—1600 A.D.,有机碳埋藏速率有所下降,为近千年来的最低水平;阶段Ⅲ:1600—1850 A.D.,升金湖有机碳埋藏速率迅速上升后,保持在较高水平;阶段Ⅳ:1850—2000 A.D.,升金湖有机碳埋藏速率继续上升,达到近千年来最高水平。进一步研究发现,升金湖有机碳埋藏速率主要受TOC含量的影响,与质量累积速率相关性较弱;阶段Ⅰ(1000—1350 A.D.)与阶段Ⅱ(1350—1600 A.D)沉积物有机质主要来源于内源水生植被,沉积环境以自然沉积为主;TOC含量增加推测与降水增多、温度下降有关,因此自然因素是此阶段湖泊有机碳埋藏速率变化的主控因素。阶段Ⅲ(1600—1850 A.D)与阶段Ⅳ(1850—2000 A.D.)沉积物有机质主要来源于外源陆生高等植被,沉积环境受人类活动影响更为强烈,如流域人口大量增长、耕地面积增加、化肥农药使用频繁等,给湖泊流域带来了大量的有机质输入,进而显著影响湖泊有机碳埋藏速率的变化。
Abstract
An 86-cm long sediment core, obtained from Lake Shengjin in the monsoon region of eastern China, served as the material for studying the mechanisms of organic carbon burial in the lake. Based on the reliable depth-age model, the sediment accumulation rate (SAR) and total organic carbon (TOC) were used to calculate the variation of organic carbon accumulation rate (OCAR) in the lake over the past millennium. The C/N ratio and n-alkanes proxies were used to indicate the sources of organic carbon. Then the influencing factors of the OCAR are discussed in terms of the organic matter source, historical climate change and human activities. Results reveal a general increasing trend of OCAR during the past millennium, with the average OCAR of approximately 24.77 g/(m2·a), and the total organic carbon storage of 3.33 Tg (1 Tg=1012 g). Specifically, the OCAR during the past millennium can be divided into four stages: Stage Ⅰ: 1000-1350 A.D., OCAR was generally low but fluctuated significantly; Stage Ⅱ: 1350-1600 A.D., OCAR decreased to the lowest level during the past millennium; Stage Ⅲ: 1600-1850 A.D., OCAR increased rapidly and maintained a relatively high level; Stage Ⅳ: 1850-2000 A.D., OCAR increased continuously and then reached the highest level. Our results further demonstrate that the sedimentary organic matter in the first two stages is mainly derived from endogenous aquatic plants. The increase in TOC content is speculated to be related to the increase in precipitation and the decrease in temperature, and the decrease in mass accumulation rate may be due to the stable hydrological conditions of the lake. Therefore, we suppose that the natural factors mainly controlled the variations of OCAR in the lake during Stages Ⅰ and Ⅱ. During Stage Ⅲ (1600-1850 A.D.) and Stage Ⅳ (1850-2000 A.D.), our result suggests that the sedimentary organic matter was mainly derived from exogenous terrestrial plants, and the sedimentary environment was more strongly influenced by human activities. The massive growth of population, the increase of cultivated land, as well as the frequent use of chemical fertilizers and pesticides, brought a large amount of organic matter into the lake basin. Therefore, we conclude that the variations of OCAR during these two stages were significantly affected by human activities, rather than natural factors.
Keywords
在全球变暖大背景下,碳排放日益增加,不同生态系统的碳循环过程受到众多学者和政府机构的关注[1-3]。IPCC最新发布的第六次评估(AR6)第一工作组报告指出,目前大气中CO2的浓度处于至少两百万年来的最高点,碳排放已显著加剧了全球气候变暖并且呈不可逆趋势[4]。由于湖泊面积较小,早期研究更关注海洋碳埋藏在全球碳循环中的作用,而忽视湖泊碳埋藏的潜力[5-7]。然而越来越多的研究表明,湖泊生态系统中也含有大量有机碳,在中小区域尺度碳循环中发挥着重要作用[8-16],根据Downing等[17]计算,全球湖泊总有机碳每年存储量约为0.05 Pg(1 Pg=1015 g)。而湖泊沉积物蓄积超过95%的湖泊碳,因此是湖泊碳的主要存储地,也是全球碳的重要源与汇[18]。而且,湖泊沉积物有着比其他碳库更高的固碳效率,Tranvik等估算全球湖泊碳埋藏速率约为0.6 Pg/a,是海洋碳埋藏速率的4倍[19]。由此可见,湖泊碳埋藏在全球及区域碳循环中的作用不容忽视。
过去千年是历史文献资料与地质沉积记录记载气候环境变化与人类活动信息相重叠的时期[20-22],包括中世纪暖期(MWP)、小冰期(LIA)和现代暖期(CWP)3个典型气候特征阶段。越来越多的研究表明,湖泊环境受人类活动影响愈发增强,因而其沉积物记录的碳埋藏过程也受到人类活动的显著影响[23-25]。因此,认识人类活动持续加强背景下湖泊碳埋藏变化的影响因素,了解较长时间尺度湖泊中碳的存储状态和变化规律[18],对深入探讨未来气候变化和人类活动对碳循环的影响具有重要意义。
过去全球的碳埋藏研究空间上主要聚焦于北美洲、欧洲和西伯利亚等高纬度地区;时间上主要集中于全新世和近百年两个时段[10-14]。我国湖泊碳埋藏研究起步相对较晚,但近十多年来学者们对我国湖泊碳埋藏开展了大量研究,使其迅速成为国内湖泊研究热点之一[26-31]。其中,长时间尺度主要集中于全新世以来湖泊碳埋藏速率与碳储量的估算,如Wang等通过综合我国82个湖泊碳累积记录及58个直接测定的湖泊有机碳埋藏数据,首次对我国湖泊的碳储量进行计算,并综述了我国自12000 cal a B.P.以来湖泊的碳埋藏速率及碳储量变化趋势[15],但时间精度相对较低;而短时间尺度主要集中于对人类活动影响下湖泊碳埋藏速率变化程度的研究[27-29]。长江下游平原地区气候环境优越,湖泊众多且湖泊生产力较高;同时该地区湖泊受人类活动影响显著,因此开展该地区历史时期湖泊碳埋藏变化及其与人类活动关系研究具有重要意义。本研究以长江下游平原地区典型浅水湖泊升金湖为研究对象,开展该湖过去千年高分辨率的碳埋藏重建工作,同时利用C/N、正构烷烃等指标重建湖泊有机碳来源的变化,并进一步结合周边气候变化及人类活动历史记录,分析影响升金湖有机碳埋藏的主要因素,探究人类活动影响前后有机碳埋藏速率的变化程度。本研究结果试图加深对我国长江下游平原地区浅水湖泊群碳埋藏的认识,并为深入认识湖泊有机碳碳埋藏变化对区域碳循环的影响提供理论参考。
1 湖区概况、样品采集与研究方法
1.1 湖区概况
升金湖(30°15′~30°28′N,116°58′~117°14′E)位于安徽省池州市东至县境内,濒临长江南岸,系长江下游河漫滩洼地经泥沙封淤积水而成,是保存较为完好的浅水河间洼地湖泊[32](图1)。湖泊东南岸为低山、丘陵,以灰岩、页岩为主;西北岸为平原,以亚黏土、砂砾为主;湖床自南向北倾斜[33]。升金湖集水面积可达1554 km2,平均水深为1.26 m,蓄水量为0.99×108 m3,自西向东分为上湖、中湖、下湖3部分(图1)。升金湖丰水期在每年的5—8月,湖泊面积超过140 km2,最高水位超过17 m;枯水期在每年的11月至次年4月,湖泊面积最小时仅为34 km2,最高水位不到10 m,有大量滩涂裸露[32,34-35]。湖泊流域属亚热带季风性气候,降水充足,地表径流丰富,四季分明,雨热同期;年平均日照时间约2000 h,年平均降水量和蒸发量分别为1600和757.5 mm,年平均气温为16.14℃,年平均无霜期约为240 d[33]。湖泊流域地带性土壤为黄红壤,沿湖四周河漫滩广泛分布潮土和水稻土,适宜农作物生长[36]。湖泊周边土地利用方式以农用地占比最高,其次为林地、草地;产业以农业为主,为当地农户提供日常所需,实现自给自足[37]。1960s升金湖下湖黄湓闸的修建,破坏了湖泊与长江的水文自然连通性,1986年始建升金湖自然保护区,1997年12月被国务院批准为国家级自然保护区[33]。

图1升金湖沉积物岩芯位置示意
Fig.1The location of sampling sediment core in Lake Shengjin
1.2 岩芯钻取与样品采集
课题组于2019年3月利用重力钻在升金湖浅水湿地区域采集一根长度为86 cm的岩芯,编号SJPC19-3,采样位置为30°20′N、117°0′E(图1)。岩芯采集后迅速运输至实验室并低温保存。由于岩芯9 cm以上部分为含水量很高的稀泥质沉积物,无法进行有效分样,因此本研究选用9~86 cm范围内的样品。岩芯9~30 cm范围内以0.5 cm为间隔进行分样,30~86 cm范围以1 cm间隔进行分样,共分得样品98个。将分得的样品用真空冷冻干燥机进行冻干处理后进行下一步分析。
1.3 研究方法
1.3.1 AMS14C年代测试
从SJPC19-3钻孔中分别挑选了深度为10、45、86 cm的3个样品,送往美国Beta Analytic公司进行AMS14C年代分析,测年材料包括植物残体和全岩有机碳,测试结果使用CALIB程序的INTCAL 13数据库进行年代校正。
1.3.2 总碳、总有机碳、总无机碳含量分析
对岩芯分样后得到的98个样品均开展总碳(TC)和总无机碳(TIC)含量的测试,其测试方法如下:利用Elementar Vario MAX cub型元素分析仪测试沉积物TC含量,依据高温催化氧化法将样品中所有C元素转换为CO2后进入非色散红外检测器(NDIR)测定,样品测试炉温为950℃,可保证样品所有C元素充分转换为CO2,提高数据准确率与测样效率[38-39];利用TOC-VCPH型总有机碳分析仪测试沉积物中TIC含量,利用加酸反应(H3PO4)将样品无机碳转化为CO2进入NDIR测定,TIC反应池温度为200℃,可保证样品中无机碳充分转化为CO2。总有机碳(TOC)含量为TC含量减去TIC含量获得。该方法经过标准样品验证,相对于提前酸化法直接测量样品中有机碳含量精度更高,测试结果也更准确可靠[38-39]。
1.3.3 正构烷烃处理与测试
将所采集的岩芯自顶部10 cm以下每隔1 cm取样进行正构烷烃指标的测试分析,共分析38个样品。分析方法如下:称取约10~15 g冷冻干燥研磨至80目(180 μm)以下的沉积物样品,利用Thermofisher ASE200型加速溶剂萃取仪提取总类脂物,萃取试剂为体积比为9∶1的二氯甲烷和甲醇混合试剂。将提取得到的总类脂物经过高纯氮气浓缩后,分别利用混合试剂(二氯甲烷∶异丙醇=2∶1,v/v)和混合试剂(乙醚∶乙酸=96∶4,v/v)作为洗脱剂过氨基硅胶柱(LC-NH2)分离中性组分和酸性组分。将分离后的中性组分分别利用混合试剂(正己烷∶二氯甲烷=9∶1,v/v)和混合试剂(二氯甲烷∶甲醇=1∶1,v/v)作为洗脱剂过硅胶柱分离正构烷烃组分和含有GDGTs的极性组分。将分离后的正构烷烃组分溶液定容至200 μL后进行GC/MS测定。GC/MS测试仪器为安捷伦HP6890N气相色谱仪,毛细柱型号为Agilent HP-1(30 m×0.32 mm×0.25 μm),毛细柱初始温度为50℃,随后以10℃/min的速率升至320℃并保留30 min,离子源温度250℃,氦气流速1.0 mL/min,采用特征碎片离子(M/Z=85)和相对保留时间定性,内标法定量。
升金湖有机质来源信息通过描述正构烷烃奇数碳链长与偶数碳链长相对丰度的碳优势指数CPI(carbon preference index)以及奇偶优势OEP(odd-even predominance)来进行解析,计算公式如下[40-42]:
(1)
(2)
1.3.4 升金湖碳埋藏速率计算
升金湖沉积物中有机碳埋藏速率(organic carbon accumulation rate,OCAR,g/(m2·a)),由升金湖每层的质量累积速率(sediment accumulation rate,SAR,g/(m2·a)),与每层的有机碳含量TOC(%)相乘得到,计算公式如下[43-45]:
(3)
式中,SR为线性累计速率(cm/a),由深度-年代模型计算得出;DBD(g/(m2·a))为沉积物干物质密度,计算方法为样品干重Md(g)除以干重体积V干(cm3)。
2 结果
2.1 年代-深度序列
表1升金湖SJPC19-3岩芯AMS14C测年结果
Tab.1 The AMS14C dating results of the core SJPC19-3 at Lake Shengjin

2.2 升金湖中长链正构烷烃分布特征及参数变化
升金湖SJPC19-3岩芯正构烷烃分布以碳数大于C26 的长链烷烃为主,中链正构烷烃相对较少,并具有明显的奇偶优势(图3)。长链正构烷烃n-C27~n-C35相对丰度变化于36.66%~90.51%之间,以奇数碳链n-C31相对丰度最高,平均为20.93%,其次为n-C29(平均为15.57%)和n-C33(平均为15.19%)。升金湖沉积物正构烷烃CPI值变化于2.11~7.14之间,平均值为4.19;OEP变化于1.95~6.03之间,平均值为3.27(附表Ⅰ)。
2.3 沉积岩芯TIC、TOC、SAR变化特征
过去千年来升金湖沉积物TIC含量变化于0.1%~1.0%之间,平均值为0.18%(图4)。1600 A.D.前大部分样品无法检测到TIC的存在,平均值仅为0.05%,有短暂时间节点突变。1600—1850 A.D.沉积物TIC含量剧烈波动上升,含量变化于0.10%~1.10%之间,平均值为0.49%,1850 A.D.后又重新下降,含量变化于0.1%~0.27%之间,平均值为0.09%。升金湖沉积岩芯TOC含量变化于0.65%~2.41%之间,平均值为1.35%,总体上呈上升趋势。TOC含量在1600 A.D.之前较低,但有微弱增加的趋势,变化范围为0.65%~1.46%,平均值为0.96%;1640—1670 A.D.期间TOC含量迅速上升,之后维持在高值范围,变化范围为0.93%~2.41%,平均值增加到1.83%。升金湖SAR变化于0.12~0.48 g/(m2·a)之间,平均值为0.17 g/(m2·a),其中在1000—1600 A.D.沉积岩芯SAR明显下降,变化于0.09~0.17 g/(m2·a),均值为0.26 g/(m2·a);从1600 A.D.之后沉积岩芯SAR逐渐趋于稳定,波动于0.16~0.42 g/(m2·a)之间,均值为0.15 g/(m2·a)。
2.4 沉积岩芯OCAR的变化特征
过去千年来升金湖OCAR变化于14.36~39.07 g/(m2·a)之间,平均值为24.77 g/(m2·a)(图4)。根据升金湖近千年OCAR的变化特征,可将其划分为以下4个阶段:阶段Ⅰ:1000—1350 A.D.,OCAR相对较低且波动较大,变化于17.47~33.62 g/(m2·a)之间,均值为24.02 g/(m2·a);阶段Ⅱ:1350—1600 A.D.,OCAR逐渐下降至最低水平,均值为20.41 g/(m2·a);阶段Ⅲ:1600—1850 A.D.,升金湖OCAR迅速升高并维持在较高水平,变化范围为15.38~38.63 g/(m2·a),均值达到27.11 g/(m2·a);阶段Ⅳ:1850—2000 A.D.,升金湖OCAR继续上升,达到了近千年来的最高水平,变化范围为24.27~39.07 g/(m2·a),均值达到29.53 g/(m2·a)。此外,将近千年来升金湖平均OCAR与升金湖面积及沉积年代相乘,得到近千年来升金湖有机碳储量约为3.33 Tg。早期岩芯OCAR变化与迅速下降的SAR变化趋势相一致(图4中黄色箭头);中后期,尤其在1600 A.D.前后岩芯OCAR变化与快速升高的TOC含量变化趋势一致(图4中红色箭头)。

图2SJPC19-3岩芯年代-深度序列
Fig.2Age-depth model of the core SJPC19-3

图3升金湖SJPC19-3岩芯中长链正构烷烃平均相对丰度
Fig.3Averaged relative abundance of sedimentary n-alkanes in core SJPC19-3 from Lake Shengjin
3 讨论
3.1 升金湖OCAR与其他区域对比
相比于长时间尺度的湖泊碳埋藏研究,升金湖近千年平均OCAR(约24.77 g/(m2·a)),高于全新世以来全球大于5000 km2大型湖泊OCAR(约5 g/(m2·a));低于全球小型湖泊(<100 km2)的OCAR(约72 g/(m2·a))[46],与我国大暖期湖泊OCAR(约22.3 g/(m2·a))较为接近[9];相较于短时间尺度的湖泊碳埋藏研究,高于长江下游平原的大型湖泊太湖(约14.4 g/(m2·a))[47]、巢湖(约7.8 g/(m2·a))[48]等近百年的碳埋藏速率,但远低于一些富营养化的小型湖泊碳埋藏速率(148~17000 g/(m2·a))[49-50]。整理并对比近千年相关湖泊的平均碳埋藏速率表明,升金湖OCAR高于同位于长江下游平原的天目湖(14.10 g/(m2·a))[51]、华中地区的东湖(17.48 g/(m2·a))[52],而与淮河下游浅水湖泊南四湖(21.32 g/(m2·a))接近[53],低于西南地区星云湖(47.14 g/(m2·a))[54]。综上,升金湖OCAR介于长江下游大型浅水湖泊与小型湖泊之间,显著低于同样受人类活动影响较为剧烈的西南地区湖泊。

图4升金湖沉积岩芯碳埋藏指标随时间变化:(a)TOC,(b)OCAR,(c)SAR,(d)TIC
Fig.4Time series of carbon accumulation proxies in sediment of Lake Shengjin: (a) TOC, (b) OCAR, (c) SAR, (d) TIC
3.2 升金湖OCAR与SAR、TOC含量的关系
湖泊的OCAR因TOC含量与SAR变化而改变,因此对各阶段升金湖岩芯OCAR、SAR和TOC含量进行相关性分析。结果表明(图5),升金湖岩芯OCAR与TOC含量呈现显著正相关性,而岩芯SAR与OCAR没有明显的相关性,说明沉积物TOC含量变化是升金湖OCAR的主要决定因素,因此分析湖区近千年来TOC含量的变化原因,是探究升金湖OCAR的关键所在。
3.3 升金湖有机碳来源分析
湖泊沉积物中的有机质来源主要分为内源有机质和外源有机质,其中内源主要包括水生细菌、藻类和大型水生植物等,而外源主要为湖泊流域内输入的陆生高等植物和土壤有机质等。本研究为定量研究升金湖内源和外源有机质对湖泊的相对贡献量,根据课题组前期已发表C/N数据(摩尔比,根据质量比乘以1.167计算得到),采用二元模式对有机质内外源进行估算,计算公式如下[8,12,55]:

图5升金湖OCAR与其他指标的相关性热图
Fig.5Correlation heatmap between OCAR and other proxies from Lake Shengjin
(4)
式中,C(i)和N(i)分别代表了第i层沉积物中的TOC和TN含量,Cl(i)和Nl(i)分别代表了第i层沉积物中来自外源的TOC和TN含量,Cu(i)和Nu(i)分别为第i层沉积物样品中来自内源的TOC和TN含量,Rl(i)和Ru(i)分别代表了第i层沉积物中来自陆源和内源的碳氮比。由上述公式可得:
(5)
式中,Rl(i)和Ru(i)根据零级近似可分别赋予Rl=Rl(i),Ru=Ru(i),通过统计得到升金湖的C/N变化范围为7.67~14.80,因此,结合前人研究经验,本文分别赋值升金湖的Rl和Ru为7和15[8],由此得出升金湖沉积物有机质中内外源的TOC含量(图6d)。结果表明,1000—1350 A.D.,升金湖沉积物有机质外源显著低于内源输入量。在1350—1600 A.D.期间升金湖沉积物有机质外源输入略有增加,达到与内源相当水平,但输入量仍较低。1600 A.D.之后,升金湖沉积物有机质外源输入量迅速增加,推测流域在1600 A.D.后有机质以外源输入为主(图6d)。

图6升金湖沉积岩芯OCAR与有机质来源变化指标对比:(a)TOC;(b)OCAR;(c)C/N;(d)内、外源TOC含量;(e)正构烷烃CPI;(f)正构烷烃OEP
Fig.6Comparison between organic carbon accumulation and source proxies: (a) TOC; (b) OCAR; (c) C/N; (d) TOC content of endogenous and exogenous sources; (e) n-alkane CPI; (f) n-alkane OEP
正构烷烃广泛存在于藻类、水生植物和高等植物等各种生物体中,具有结构稳定、不易降解的特点,描述正构烷烃奇数碳链与偶数碳链相对比例的碳优势指数CPI和奇偶优势OEP能够记录有机质的来源信息[40]。由于浮游藻类与水生植物奇偶优势不明显,而陆生高等植物因表皮叶蜡含有大量奇偶优势明显的n-C27~n-C33的高碳数正构烷烃,其CPI值一般大于5,因此CPI与OEP值越大,往往表明陆生高等植物对湖泊沉积有机质贡献量越多[41-42]。升金湖岩芯正构烷烃CPI值变化趋势结果表明,1600 A.D.前岩芯正构烷烃CPI和OEP处于低值阶段,1600 A.D.后岩芯正构烷烃CPI和OEP值迅速升高,CPI基本维持在5以上(图6e、f),这表明1600 A.D.前湖区以水生植物为主,1600 A.D.后湖区陆生高等植物比例迅速增加。因此,通过上述指标变化,推测升金湖流域1600 A.D.前沉积物有机质主要来源于内源水生植物,1600 A.D.后沉积物有机质主要来源于外源陆生高等植被。同时将升金湖TOC含量与外源进行相关性分析表明,二者呈现高度正相关性(图5),即过去千年来升金湖TOC含量随着外源有机质输入的增加而呈上升趋势(图6中箭头)。
3.4 升金湖OCAR影响因素分析
上述研究表明湖泊OCAR受TOC含量、SAR变化的影响,而二者又因气候(温度、降水等)、人类活动等因素影响而发生改变[8-16],因此本研究结合周边气候变化及历史人类活动记录,分析过去千年影响升金湖沉积物OCAR变化的主要因素。
3.4.1 气候因素
有机物经过复杂的风化、搬运沉积等地质过程,在湖泊中持续累积并形成沉积有机碳,因此OCAR与气候变化密切相关。前人研究表明,气候对OCAR的影响主要表现在温度和降水两个方面,温度和降水可以直接影响湖泊水体及沉积物的水热条件与微生物活性,从而影响有机碳的矿化分解速率,进而影响到沉积物有机碳埋藏[27];还能够通过影响湖泊流域的植被类型、生物量、湖泊特征(水位、水动力、面积等)进而影响湖泊OCAR[56-57]。
(1)降水
诸多研究表明,当流域降水增加时,植被生长繁茂使得流域土壤碳储量增加,同时流域地表径流输入增强,河流中的溶解有机碳搬运入湖增多,进而造成湖泊OCAR的增加[58]。但也有研究认为,湿润环境下植被覆盖率提高、水动力减弱、湖泊水位上升、面积扩张,也可能造成流域内的侵蚀作用减弱以及单位面积湖泊水体TOC含量减少,进而造成OCAR下降[59]。因此降水对湖泊OCAR的影响仍存在较大争议。
本研究通过沉积物化学蚀变指数(chemical index of alteration,CIA)以及课题组前期研究成果中Ti元素、平均粒径等指标来揭示升金湖流域水文信息[33,60],进而探究湖区水文环境与OCAR间的关系(图7)。前人研究表明,湿润环境下,Ti元素容易富集而含量相对升高,而CIA则通过元素活动组分与惰性组分的关系来指示源区风化程度,CIA值越高表明源区化学风化越强、水热条件优越;反之,CIA值越低反映化学风化作用较弱、水热条件较差[61-62]。而平均粒径则可指示采样点在湖泊中的位置,进而反映湖泊水位与水动力情况,平均粒径较小时表明采样点距湖岸较远、湖泊水位较高、水动力较弱;反之表明采样点距湖岸较近、湖泊水位较低、水动力较强[63]。
结果表明,在1600 A.D.前,升金湖Ti、CIA、平均粒径等指标与TOC、TIC含量的相关性较好,即当湖区水热条件较好、水位较高时,湖区TOC含量往往会增加,而水热条件较差、水位下降则会导致TOC含量减少、TIC含量增加,说明此时段升金湖沉积物中TOC、TIC含量与流域水文环境关系较为密切。从阶段Ⅰ(1000—1350 A.D.)到阶段Ⅱ(1350—1600 A.D.)沉积物Ti、CIA值明显升高,平均粒径与TIC含量则逐渐降低,反映湖区向湿润化转变、湖泊水位增高、采样点附近水动力减小的水文环境,该环境有利于水生植物的繁衍以及径流输入有机质的增加,进而造成湖泊水体中有机质的潜在来源增多,同时水位升高也减少了采样点位置在氧气中的暴露时间,从而减缓了有机质的分解与氧化。这与湖区TOC含量在这一时期明显增加的趋势相一致。值得注意的是,阶段Ⅰ前期沉积物SAR显著高于其他时段,较高的SAR与较低的TOC含量,导致SAR、OCAR出现协同变化趋势,可能的原因是湖泊流域水热环境波动剧烈,且表现为高温干旱趋势,导致采样点附近侵蚀作用较强、沉积物干密度较高,有利于SAR的提升;而之后流域的水热环境趋于稳定。
阶段Ⅲ(1600—1850 A.D.),岩芯Ti元素、CIA等指标显著降低,平均粒径与TIC含量剧烈波动上升,指示流域水文条件较为干旱;阶段Ⅳ(1850—2000 A.D.),各指标指示流域水文环境开始转湿,与周边湖泊降水重建结果较为一致[64-65](图7g、h),推测湖泊水体中有机质潜在来源在这两个阶段先减少后增加。但本研究结果显示岩芯TOC含量表现为持续快速升高趋势,可能说明水文环境并非是该阶段湖泊TOC含量迅速增加的主导因素。值得注意的是,本研究对象升金湖属于浅水湖泊且岩芯钻孔地点位于地形平坦的滨岸湿地,阶段Ⅲ气候相对干燥、水位下降,可能导致采样点位置露出水面,陆生或挺水植被大量生长,从而造成TOC含量增加。
(2)温度
研究发现,温度升高,有机碳的矿化速率加剧,导致有机碳埋藏效率下降;温度降低,有机碳矿化分解速率变慢,导致TOC含量的增加[27-28]。升金湖所处华东地区的温度距平[66]重建结果显示(图7i),该地区温度在约1000—1850 A.D.期间呈总体降温趋势。已有研究也表明,我国各区域的温度变化在近千年均表现出长期的下降趋势,这反映中世纪暖期与现代暖期两个温暖时期中间存在一段持续时间较长、相对较冷的小冰期[67-68]。因此除了上述流域降水量因素外,1000—1850 A.D.升金湖沉积物中TOC含量的增加还可能受到了温度下降导致的矿化度减弱的影响。

图7升金湖OCAR与自然影响因素分析:(a)TOC;(b)OCAR;(c)TIC;(d)岩芯平均粒径[60];(e)化学蚀变指数[60];(f)Ti[60];(g)基于泥炭重建的大九湖降水量变化[64];(h)基于色度A*/L*重建的女山湖降水量变化[65];(i)华东温度距平[66]
Fig.7Analysis of natural influencing factors of organic carbon accumulation rate of Lake Shengjin: (a) TOC; (b) OCAR; (c) TIC; (d) average grain size; (e) CIA; (f) Ti; (g) precipitation change of Dajiuhu Peat based on pollen analysis; (h) Lake Nvshan precipitation change based on A*/L* reconstruction; (i) temperature anomaly of eastern China
3.4.2 人类活动
人类活动主要通过改变土地利用方式影响湖泊碳循环,例如农业活动会导致流域内化肥施用量增加、围垦与耕地面积扩张,不仅能够为湖泊输入更多的营养盐,造成入湖氮磷等元素的增加,还可能破坏地表植被,造成入湖泥沙量增多,进一步导致土壤碳向湖泊输入增多、SAR提高,最终造成OCAR增加[23-25,69]。
本研究将升金湖OCAR变化指标与流域内历史人类活动指标进行对比(图8)。对比结果发现,阶段Ⅰ(1000—1350 A.D.)与阶段Ⅱ(1350—1600 A.D.)安徽省人口数量一直维持在低值状态[59],说明当时人类活动相对较弱,因此对区域环境的影响相对较小(图8g)。而阶段Ⅲ(1600—1850 A.D.)特别是明末清初时段,安徽省人口数快速增长[70],这与升金湖OCAR、TOC含量在此时段的快速上升呈现出高度相关性。升金湖流域以农为本,明清时期人口的快速增长给土地带来巨大的压力,人均耕地面积逐年减小,人口所带来巨大的生存压力只能通过扩大耕地面积来解决,进而导致流域内土地利用方式的变化[71]。大量研究也证明,明清时期皖江流域圩田兴起并发展进入全盛时期,圩田地域分布广泛,导致水土流失加剧[72-73]。为进一步理解人类农业活动对升金湖流域OCAR的影响,本研究通过历史文献资料粗略整理出中国耕地面积与池州府围垦耕地面积的相关数据[74-77](图8h,i),以及利用沉积岩芯重金属砷(As)的含量来反映流域内农业生产和化肥农药使用情况[78](图8f)。指标对比结果表明:阶段Ⅲ(1600—1850 A.D.),中国耕地面积与池州府耕地面积快速增加,岩芯As含量显著上升,并保持在较高水平,这可能反映出农业生产活动在这一时间段内迅速加强,土地利用方式发生改变以及生产过程中化肥使用开始增加。同时,对比湖泊氮磷硫[33]等元素变化特征也发现,此时段流域内较高的耕地面积可能导致氮磷硫等营养元素快速增加,为湖泊输入了大量的营养盐。以上指标变化可说明,自阶段Ⅲ开始人类农业生产已显著影响到升金湖的TOC沉积过程。阶段Ⅳ(1850—2000 A.D.),升金湖OCAR、TOC含量均继续升高,多指标也指示该阶段人类活动进一步加剧,可能反映出人类活动对湖区碳埋藏的影响进一步扩大。值得注意的是,1600 A.D.受人类活动显著影响后,升金湖平均OCAR约为(28.06 g/(m2·a)),较1600 A.D.之前(平均OCAR约22.62 g/(m2·a))增长约21%。相较于其他湖区受人类活动显著影响的湖泊,增长幅度较慢,这也与前人研究结果相一致[16,44],主要原因是长江下游浅水湖泊群由于气候、浅水等性质而发生强烈的矿化和分解作用,因此导致该湖区在人类活动强烈、有机碳输入较高的背景下累积量相对较少,但不可否认的是该湖区碳埋藏依然是一个巨大的区域性碳汇。

Fig.8Relationship between organic carbon accumulation rate changes in Lake Shengjin and historical human activities: (a) TOC; (b) TN; (c) TP; (d) S; (e) OCAR; (f) As; (g) The historical population of Anhui Province[70]; (h) The estimated historical cultivated land areas in China; (i) The estimated historical cultivated land areas of Chizhou Prefecture
本研究探讨了流域内气温降水以及人类活动对升金湖OCAR变化的可能影响,依据上述探讨分析内容,总结升金湖过去千年来碳埋藏变化的4个阶段大致可分为两类模式(图9):
阶段Ⅰ(1000—1350 A.D.)与阶段Ⅱ(1350—1600 A.D.),升金湖流域的OCAR整体处于低水平,有机质含量较低。此时段人类活动强度较低,OCAR受气候条件的影响更为强烈。温度下降、降水增多可引起有机碳矿化分解速率变慢以及有机质潜在来源增加;气候干旱则会引起TIC含量上升,TOC含量下降;而水热环境波动剧烈则可能引起SAR的增加。综上,可认为此时段升金湖OCAR以自然沉积为主(图9)。阶段Ⅲ(1600—1850 A.D.)与阶段Ⅳ(1850—2000 A.D.),升金湖流域的OCAR显著提高,TOC含量迅速升高,使得升金湖OCAR整体达到高水平。此时段流域的水文条件、气温等对升金湖OCAR的响应相对较弱,而人口快速增加、耕地面积扩大、农业施肥增多等可能是沉积物TOC含量与SAR快速上升的主要原因。综上,可认为此时段升金湖OCAR受到人类活动的影响更为强烈(图9)。

图9升金湖OCAR变化机制示意(箭头粗细代表了影响程度,(+)代表积极效应,(-)代表消极影响)
Fig.9The sketch of organic carbon accumulation rate mechanisms of Lake Shengjin (arrow thickness represents impact degree, (+) represents positive effect, (-) represents negative effect)
4 结论
1)升金湖过去千年OCAR约为24.77 g/(m2·a),总有机碳储量约为3.33 Tg。其中阶段Ⅰ(1000—1350 A.D.)与阶段Ⅱ(1350—1600 A.D.),升金湖OCAR处于低水平;从阶段Ⅲ(1600—1850 A.D.)开始迅速上升并在阶段Ⅳ(1850—2000 A.D.)达到峰值。此外升金湖OCAR的变化主要是由岩芯TOC含量的变化引起,而受质量累积速率的影响较小。
2)升金湖流域1600 A.D.前沉积物有机质主要来源于内源水生植物,1600 A.D.后沉积物有机质主要来源于外源陆生高等植被,因此1600 A.D.后升金湖TOC含量的迅速增加可能与湖区陆生高等植物大量输入有关。
3)升金湖流域有机碳埋藏沉积环境存在两种模式,其中阶段Ⅰ与阶段Ⅱ以自然沉积为主,湖区OCAR变化推测与降水增多、温度下降有关。阶段Ⅲ与阶段Ⅳ则受人类活动影响更为强烈,湖区人口大量增长、耕地面积增加、施肥农药频繁等给湖泊流域带来大量的有机质,进而影响到湖泊OCAR变化。
5 附录
附表Ⅰ见电子版(DOI: 10.18307/2025.0251)。