王君兰,邓 伟,张平宇,张家武
(1:中国科学院兰州文献情报中心,兰州730000)
(2:甘肃省地质调查院,兰州730000)
(3:兰州大学西部环境教育部重点实验室,西部环境与气候变化研究院,兰州730000)
湖泊沉积物中的有机质含量是反映湖泊或流域生产力和湖泊沉积环境的常规指标,其碳同位素(δ13Corg)可以反映湖泊或流域内植被类型的变化,从而反映气候环境的变化.因此有机质含量及其δ13Corg常被用来恢复区域的古环境、古植被变化情况[1-7].常见的有机质含量的测定方法有烧失法(LOI)、滴定法、元素分析仪测定法(TOC)等,前人的研究结果表明,尽管不同方法测定的值存在差异,但变化趋势相同[8-13],均能代表沉积物中有机质含量的变化情况.碳同位素δ13Corg的测定一般选择全岩样品,并需要去除样品中较大的植物残体及无机碳酸盐,避免由此带来样品的不均一性,同时较好地保留有机质组分,才能使实验结果较准确地反映古气候环境的变化[14-20].但对于有机质含量变化较大的湖泊沉积物,如出现较多植物残体时,目视挑选去除植物残体(全样不过筛,研磨足够细且均匀即可)和筛除掉植物残体的细颗粒组分(如过120目筛和360目筛等)对样品均一性的影响、其同位素组成的差别以及对于古环境解释的影响尚无详细研究和评估.本文选取内蒙古额济纳旗黑河下游嘎顺诺尔GXN剖面的样品,通过选用不同大小颗粒组分的样品进行有机质含量及有机碳同位素的对比实验,讨论不同组分有机质碳同位素结果的差异,为湖泊沉积物有机质碳同位素分析提供优选方案.
本文采用的沉积物样品于2009年11月在嘎顺诺尔东部(42°27'03″N,100°54'05″E)采集,海拔 890 m.四周是干涸的湖底,在湖底上挖掘探坑剖面(GXN)至地下水出露为止,探坑深2 m(图1).在剖面上以1 cm间隔采样,共采集样品200个,样品装入自封袋密封后运回室内,经真空冷冻干燥仪干燥后待分析.
图1 嘎顺诺尔GXN剖面图及其岩性Fig.1 Photo of the GXN section and lithological profile
1.2.1 烧失量法 将坩埚用稀盐酸清洗后放在预先烧至550℃的马弗炉中灼烧4~6 h,恒重后称重,称约1 g干样(精确到0.0001 g)放入坩埚中,将坩埚放入105℃烘箱干燥12 h,称重后在马弗炉中550℃灼烧2~3 h,关闭电源,打开炉门,待部分热量散失后,用长钳取出坩埚置于干燥器中冷却.待温度降至室温时在分析天平上称重.烧失部分的重量占样品干重的百分比即为烧失量(LOI)[21].样品的有机质含量计算公式为:LOI(%)=(550℃时烧失量/干重)×100%.
1.2.2 有机质含量及有机碳同位素的测定 1)将经冷冻干燥过的样品分为3份,一份挑取掉对土壤有机质有影响的大的植物残体,未经过筛直接将全样用玛瑙研钵研磨均匀;另一份样品挑取掉植物残体后用蒸馏水分层湿筛,上层筛孔孔径为120目(125 μm),下层孔径为360目(45 μm),将通过下层45 μm孔径的细粒部分烘干后研磨均匀;另外挑选含植物残体较多的样品(91~115 cm)过120目筛,烘干后用玛瑙研钵研磨均匀;
2)透水坩埚置于马弗炉中550℃灼烧4 h左右,取出置于干燥器中冷却;
3)取上述适量样品(参照烧失量结果视有机质的多少而定)放在透水坩埚中,将透水坩埚置于培养皿中,用10%的盐酸滴入透水坩埚中,滴加至反应缓和时,加入过量的盐酸充分反应24 h,目的是去除样品中的碳酸盐;
4)将透水坩埚放在能将其支起并透水的架子上,在透水坩埚中加入蒸馏水洗至中性(用pH试纸测试),将透水坩埚放在烘箱中,75℃烘24 h,将充分干燥后的样品混合均匀,用玛瑙研钵研磨至120目左右,装入小试管中编号保存备用;
5)用万分之一天平称取样品,并以干净的特制小锡舟紧密包好(包样过程尽量将锡舟中空气排挤干净,否则空气中的CO2会影响到样品测定值的准确性),然后放入旋转式自动载样盘中进行上机(元素分析仪-同位素比例质谱仪-EA-IRMS)测试,误差<0.2‰.对于每个样品,必须予以准确称量,燃烧释放的CO2在通过色谱柱后与标准样品的CO2气体的峰积分面积相比较,从而得到待测样品的碳含量(以元素C的质量百分数形式给出),即总有机碳(TOC)[22].
以上实验均在兰州大学西部环境教育部重点实验室完成.
准确测定沉积物中总有机碳含量,需要有效地将有机碳和无机碳区分开来.烧失法采用特定温度下(550℃)的烧失量来表示有机质含量;研究表明,高镁碳酸盐在400℃以下会发生分解,而一些难降解的有机物要在1050℃以上才能被氧化,所以该方法在测量样品时产生误差较大,且不易控制.而元素分析仪法采用无氧化性的盐酸溶液溶解掉沉积物中的无机碳酸盐,然后上机测试有机碳含量.实验过程中也会造成偏差,比如去除碳酸盐过程可能将可溶性有机物溶解,但相对而言更易控制[13].
从图2可以看出,代表有机质含量的3条曲线(a、b和c)形态相似,说明不同组分不同分析方法测得的有机质含量变化趋势一致.整个剖面有机碳含量较低,仅在85~115 cm段有机质含量较高,与野外相应位置观察到该段有大量芦苇的植物残体相一致.对太湖现代沉积物的研究[23]表明,在生有芦苇等水生植物的水域有机质含量高,烧失量也较高;在敞水区或湖心区,有机质含量低,烧失量亦低,烧失量与有机质含量之间存在密切的关系.剖面中植物残体含量高段烧失量大致为全样有机碳含量的2倍,这与Dean等[9]的研究结果一致.
全样样品在85~115 cm(植物残体含量高)段的TOC(图2b)较细颗粒组分的TOC(图2a)含量高,其余层位全样比细颗粒组分有机质含量稍低.这表明在有机质含量较低的湖泊沉积物中,不溶解的有机质颗粒主要为较细的颗粒,过筛去除了比重较大的沉积物颗粒后,该部分有机质所占的比例相应增加.而含植物残体多的湖泊沉积物中,过筛过程中多数粗颗粒有机质被筛除,过筛后的细颗粒组分TOC含量较全样低.尽管实际研究中烧失量因测定起来快速简便、经济的特点而经常用来代表样品有机质含量,但因前处理不同(未去除碳酸盐),在有机质含量较低的沉积物中,烧失量波动较大(图2c),可能与不同样品碳酸盐含量、可溶有机组分不同以及有机质烧失过程中部分碳酸盐分解有关.因此用烧失量来粗略代表有机质含量时,LOI曲线与其他指标对比时,其细节变化不具有对比意义,这在研究中应引起注意.
烧失法和元素分析仪法测得的TOC含量的相关性分析(图3)表明,同一样品(如全样,图3a)TOC(元素分析)和LOI尽管统计结果相关(n=200),但部分数据点仍较为离散.过筛后的TOC与全样的LOI(图3b)相关性最低,这与上述有无前处理过程有关;元素分析仪测得的全样和过筛的细颗粒组分的TOC相关性最好(图3c),在有机质含量低时,二者较为接近,有机质含量高时,二者数据也较为离散,说明取不同组分进行元素分析测量时,样品的均一性对测量结果影响较大.过筛的细颗粒组分在有机质含量低时段出现一定的波动(图2a),而在有机质含量高的时段波动反而较小.这反映了过筛过程对有机质的分选作用:细颗粒组分相对富集,粗颗粒的有机质(如残体)得到有效去除.
嘎顺诺尔GXN剖面沉积物全样和细颗粒组分的δ13Corg值曲线形态相似(图2d、e),细颗粒组分有机碳同位素变化范围为 -25.51‰ ~ -19.23‰,平均值为 -23.03‰;全样有机碳同位素值在 -24.88‰ ~-18.36‰之间,平均值为 -22.63‰.全样比细颗粒组分的有机碳同位素值偏正,两者的差值(D-value,图2f)变化范围为-1.51~3.13,平均值为0.39.两者之间的差值也可能与有机质组分赋存的颗粒范围有关.为了进一步分析不同粒径的样品测得的碳同位素值结果,选取植物残体含量较多的91~115 cm段样品,增加过120目筛组分的样品进行测试,并与全样和过360目筛的细颗粒组分进行对比,结果如图4.
图2 嘎顺诺尔GXN剖面全样和细颗粒组分的有机质含量及δ13Corg值曲线Fig.2 TOC content and δ13Corgof the bulk sample and the fine-grained components of the GXN section
图3 烧失法和元素分析仪法测得有机质含量结果的相关关系Fig.3 Relationship of organic contents by LOI and elemental analysis methods
同一样品不同粒径组分的有机碳δ13Corg值趋势基本一致(图4),但全样的δ13Corg值(图4a)波动最大,过120目筛组分的δ13Corg值(图4c)次之,过360目筛组分的δ13Corg值(图4b)波动最小.对同一样品不同组分的δ13Corg值进行相关分析,结果如图5.全样的δ13Corg值与过不同孔径筛子组分的δ13Corg值相关程度均不高(图5a、b),但过120目筛和360目筛后的沉积物δ13Corg值相关性十分显著(R2=0.94)(图5c),证明过筛的过程对有机质进行了分选,使样品更加均匀,同时分选出的有机质颗粒可能更集中于某一相同来源的有机质(如陆源的C3植物),使过360目筛的最细颗粒的δ13Corg值偏负一些.颗粒粒径对固体样品的均匀混合也会产生一定影响.一般而言,颗粒粒径越小,样品越容易混合均匀,相应的分析结果离散度越小[24].如图5显示全样比细颗粒组分的离散度要大.因此本文的实验结果表明,随着湖泊沉积物全样被不同孔径筛分后,有机质含量及其δ13Corg值均发生变化,随着过筛的孔径变细,有机质中的细颗粒碎屑相对富集,残体影响减小,同位素结果波动减小.因此使用湖泊沉积物有机质 δ13Corg值进行区域对比时,应查明所使用样品的处理方式,相同处理方式下的δ13Corg结果更具有可比性.波动较大的δ13Corg结果可能是样品不均匀所致.对比结果表明,过360目筛后的样品,细颗粒有机质相对富集,样品混合更均匀,有利于燃烧更为充分,反应更加完全,更有利于分析区域有机质δ13Corg的稳定变化趋势,从而排除剧烈波动给指标对比和区域对比带来的干扰.此外,嘎顺诺尔GXN剖面沉积物全样和细颗粒组分(过360目筛)的TOC与 δ13Corg值的相关性(图6)均很低,说明有机质碳同位素值与有机质含量关系不大.
图4 嘎顺诺尔GXN剖面沉积物全样和细颗粒组分的δ13Corg值对比曲线Fig.4 Comparison of the δ13Corgof the bulk sample and the fine-grained components of the GXN section
图5 嘎顺诺尔GXN剖面沉积物全样和细颗粒组分的δ13Corg值的相关关系Fig.5 Relationship of δ13Corgvalues between bulk sample and fine-grained components of the GXN section
通过分析嘎顺诺尔湖泊沉积物不同组分有机质含量及其有机碳同位素值,得到以下结论:
1)湖泊沉积物有机质含量较高时,烧失量LOI大致相当于元素分析仪法得到的全样TOC的2倍;有机质含量较低时,LOI与TOC差别较大,因为TOC分析前处理时去除了可溶的有机质组分;全样的TOC较细粒TOC略低,显示出过筛过程中对有机质颗粒的分选作用.两种方法存在很好的相关性,比较而言,烧失量测定方法仅能粗略地用来代表有机质含量,与其他指标对比时细节变化可能不代表环境意义;用元素分析仪测定得到的TOC结果更准确.
2)嘎顺诺尔GXN剖面沉积物全样和细颗粒组分的δ13Corg值对比表明,全样的有机碳同位素值存在较大波动,随着过筛孔径变小,δ13Corg值趋于稳定,表明δ13Corg值与样品均匀程度有关,且过360目筛的细颗粒组分的碳同位素值较全样或过120目筛后组分的δ13Corg偏负.这一差异与有机质组分赋存的颗粒范围有关,筛后分选(富集)了低有机质含量段的细颗粒有机质,也剔除了大颗粒残体的影响,因此我们认为过360目筛的细颗粒组分更有利于充分反应,且可获得更有利于区域对比的同位素值.不同湖泊沉积物δ13Corg值对比时,应注意研究使用样品的前处理方式是否相同,相同处理方式下的δ13Corg值更具有对比意义.
图6 嘎顺诺尔GXN剖面沉积物全样/细颗粒组分的TOC与δ13Corg值的相关关系Fig.6 Relationship of TOC and δ13Corgvalue of bulk sample or fine-grained components of the GXN section
[1]吴敬禄,王苏民.湖泊沉积物中有机质碳同位素特征及其古气候.海洋地质与第四纪地质,1996,16(2):103-109.
[2]沈 吉,王苏民,羊向东.湖泊沉积物中有机碳稳定同位素测定及其古气候环境意义.海洋与湖沼,1996,27(4):400-403.
[3]沈 吉,吴瑞金,安芷生.大布苏湖沉积剖面有机碳同位素特征与古环境.湖泊科学,1998,10(3):8-12.
[4]余俊清,王小燕,李 军等.湖泊沉积有机碳同位素与环境变化的研究进展.湖泊科学,2001,13(1):72-78.
[5]张恩楼,沈 吉,夏威岚等.青海湖沉积物有机碳及其同位素的气候环境信息.海洋地质与第四纪地质,2002,22(2):105-108.
[6]马 龙,吴敬禄.安固里淖湖积物中总有机碳含量及其碳同位素的环境意义.自然资源学报,2009,24(6):1099-1104.
[7]吴敬禄,王苏民.湖泊沉积物有机质所揭示的环境气候信息.湖泊科学,1996,8(2):113-118.
[8]王浩然,陈业材.湖泊沉积物有机碳分析的最优化方法.贵州环保科技,1996,2(4):21-24.
[9]Jr Dean WE.Determination of carbonate and organic matter in calcareous sediments and sedimentary rocks by loss on ignition:Comparison with other methods.Journal of Sedimentary Research,1974,44:242-248.
[10]Heiri O,Lotter AF,Lemcke G.Loss on ignition as a method for estimating organic and carbonate content in sediments:reproducibility and comparability of results.Journal of Paleolimnology,2001,25:101-110.
[11]Santisteban JI,Mediavilla R,Lopez-Pamo E et al.Loss on ignition:a qualitative or quantitative method for organic matter and carbonate mineral content in sediments?Journal of Paleolimnology,2004,32:287-299.
[12]Beaudoin A.A comparison of two methods for estimating the organic matter content of sediments.Journal of Paleolimnology,2003,29:387-390.
[13]卢凤艳,安芷生.鹤庆钻孔沉积物总有机碳、氮含量测定的前处理方法及其环境意义.地质力学学报,2010,16(4):393-401.
[14]曹蕴宁,刘卫国,宁有丰等.土壤有机碳同位素样品制备过程的影响因素讨论.地球化学,2005,34(4):395-404.
[15]王金权.样品预处理对有机碳同位素分析结果的影响.古生物学报,2005,44(3):472-477.
[16]吴 夏,黄俊华,白 晓等.沉积岩总有机质碳同位素分析的前处理影响.地球学报,2008,29(6):677-683.
[17]刘 强,刘嘉麒,陈晓雨等.18.5kaB.P.以来东北四海龙湾玛珥湖全岩有机碳同位素记录及其古气候环境意义.第四纪研究,2005,25(6):711-721.
[18]Midwood AJ,Boutton TW.Soil carbonate decomposition by acid has little effect on δ13C of organic matter.Soil Biology and Biochemist,1998,10(11):1301-1307.
[19]Schubert CJ,Nielsen B.Effects of decarbonation treatments on δ13C values in marine sediments.Marine Chemistry,2000,72(1):55-59.
[20]Lorrain A,Savoye N,Chauvaud L et al.Decarbonation and preservation method for the analysis of organic C and N contents and stable isotope ratios of low-carbonated suspended particulate material.Analytica Chimica Acta,2003,491(2):125-133.
[21]张家武,金 明,陈发虎等.青海湖沉积岩芯记录的青藏高原东北部过去800年以来的降水变化.科学通报,2004,49(1):10-14.
[22]饶志国,陈发虎,曹 洁等.黄土高原西部地区末次冰期和全新世有机碳同位素变化与C3/C4植被类型转换研究.第四纪研究,2005,25(1):107-113.
[23]袁旭音,陈 俊,季俊峰等.太湖现代沉积物的物质组成和形成条件分析.南京大学学报:自然科学版,2002,38(6):756-765.
[24]喻 涛,李春园.盐酸温度时间及粒径对海洋沉积物碳酸盐去除的影响.热带海洋学报,2006,25(6):33-38.