察汗淖尔表层沉积物GDGTs分布特征及影响因素*

卢沛琪,薛 刚,郑艳红

(西北大学地质学系大陆动力学国家重点实验室,西安 710069)

重建大陆古气候与古环境一直是地质学研究的热点问题之一,这对更好地了解地球历史,预测未来气候变化具有重要意义。温度作为气候变化的重要因子,其变化显著地影响着气候带的分布、动植物的生存、人类健康及农业生产等。因此,厘清温度的定量变化至关重要。然而,当前依靠仪器记录温度变化在时间和空间上存在限制,并且近年来温度还广泛受到人类活动的影响,还不足以理解自然背景下温度变化的规律及其变化机制。

图1 GDGT化合物结构图Fig.1 Structure of GDGTs compounds

然而,湖泊沉积物brGDGTs来源复杂,既可以在原位水柱中产生[23-24],又可以接受土壤输入[25-28],将基于土壤的校正公式应用于以原位生产brGDGTs为主的湖泊沉积物时,重建的温度远低于实际值[29],因而在东非湖泊[30-31]、中国湖泊[5,25,28]乃至南北极湖泊[7,32]等都建立了适合湖泊原位生产的温度转换公式。因此,在进行湖泊brGDGTs古气候研究时,明确brGDGTs来源至关重要。

1 材料和方法

1.1 样品采集

图2 察汗淖尔采样点示意Fig.2 The sampling sites location of Lake Chahannaoer

图3 察汗淖尔月均温度(点线图)及降水量(柱状图)(来自1981-2010年中国气象局网站记录的数据集)Fig.3 Monthly mean temperature (Point-Line Chart) and precipitation (Histogram) of Chahannaoer from 1981 to 2010 (Dataset from the website of China Meteorological Administration)

1.2 pH值与含水量的测定

根据Weijers等的方法测量pH值,将样品与超纯水以1∶2.5比例充分混合,经离心机分离后提取上清液,再使用梅特勒托利多pH计进行pH值测定,共测定3次,取平均值,测量误差在±0.03之内。土壤含水量利用差量法进行测定,称量一定质量的样品(m1),将样品放置在冷冻干燥机干燥48 h,再次称量样品质量(m2),则样品含水量(SWC)=(m1-m2)/m2。

1.3 GDGTs抽提与测定

将冷冻干燥后的样品研磨成粉末(20 g左右),在室温下,使用二氯甲烷∶甲醇(DCM∶MeOH)混合溶剂(9∶1/V∶V)在快速溶剂萃取仪中抽取总脂质提取物,并通过旋转蒸发仪浓缩至1～2 mL。分别将正己烷(Hexane)和甲醇(MeOH)作为洗脱液,在硅胶柱内分离提取物的非极性组分和极性组分,收集后的组分在氮气流中进行干燥。使用正己烷溶解极性组分后通过0.45 μm聚四氟乙烯(PTFE)过滤器过滤,在氮吹仪下干燥过滤后的组分,最后定容至500或1000 μL,以待进行后续测试分析,前处理在西北大学地质学系完成。GDGTs测试在中国科学院地球环境研究所完成,使用的是岛津三重串联四极杆液相色谱质谱联用仪,单个样品总运行时间为2 h,进样量为50 μL,分析时加入10 μL的内标物(C46),在单离子扫描模式下对质子化离子进行扫描,通过测定各个化合物离子的峰面积对brGDGTs和isoGDGTs不同组分进行定量分析,根据内标物的丰度计算得到各组分的相对丰度。

1.4 计算公式

结合土壤类型和GDGTs丰度分布特征,本文选择了如下计算公式; 所有样品中IIIb和IIIc含量较少,因此使用Peterse等定义的MBT′和CBT指数计算GDGTs的甲基化和环化程度[9]:

MBT′=(Ia+Ib+Ic)/ (Ia+Ib+Ic+IIa+IIa′+IIb+IIb′+IIc+IIc′+IIIa+IIIa′)

(1)

CBT=-lg[(Ib+IIb+IIb′)/(Ia+IIa+IIa′)]

(2)

MAAT=0.81-5.67CBT+31.0MBT′

(3)

pH=7.90-1.97CBT

(4)

(5)

CBT′=lg[(Ic+IIa′+IIb′+IIc′+IIIa′+IIIb′+IIIc′)/(Ia+IIa+IIIa)]

(6)

(7)

pH=7.15+1.59CBT′

(8)

此外,我们还应用了Yang等针对中国干旱地区土壤提出的温度校正公式[41],以及Wang等建立的适合北方土壤的区域校正[42]:

MAAT(Yang)=20.9-13.4 f(IIa+IIa′)-17.2 f(IIIa+IIIa′)-17.5 f(IIb)+11.2 f(Ib)

(9)

(10)

pH(Wang)=1.65CBT′+6.93

(11)

根据Sinninghe Damsté等提出的公式计算四甲基化、五甲基化和六甲基化brGDGTs丰度[43],用于分析brGDGTs来源:

Tetramethylated brGDGTs=Ia+Ib+Ic

(12)

Pentamethylated brGDGTs=IIa+IIb+IIc+IIa′+IIb′+IIc′

(13)

Hexamethylated brGDGTs=IIIa+IIIb+IIIc+IIIa′+IIIb′+IIIc′

(14)

使用De Jonge等定义的CI指数以评估生物群落对brGDGTs分布的影响[44]:

CI=Ia/(Ia+IIa+IIIa)

(15)

2 结果

2.1 IsoGDGTs组成与丰度

所有样品均检测到丰富的isoGDGTs和brGDGTs,总含量在1.50～29.22 ng/g之间。如图4所示,湖泊表层沉积物与土壤呈现相似的分布特征,isoGDGTs含量分别在0.14～2.45、0.07～11.14 ng/g之间,主要以crenarchaeol为主(平均占比为48%),其次是GDGT-0(29%)、GDGT-2(7%)、GDGT-1(6%)、crenarchaeol′(5%)和GDGT-3(4%)。

2017年扬中市依托中兴物联网项目，深化智慧停车设备系统的开发管理，进一步推进了智慧停车工程建设发展。

图4 察汗淖尔湖泊及其周围表层土壤GDGTs分布(误差棒为标准偏差)Fig.4 Distribution of GDGTs in surface sediments and surrounding soils of Lake Chahannaoer (Error bars indicate the standard deviation)

2.2 BrGDGTs组成与丰度

大部分样品brGDGTs丰度较高,平均占GDGTs总量的80%以上。湖泊沉积物中brGDGTs含量在1.60～13.27 ng/g 之间,周围土壤中总含量是0.24～18.08 ng/g。brGDGTs中五甲基化brGDGTs含量最多(图4),其次是六甲基化和四甲基化brGDGTs。IIa及IIa′占湖泊沉积物brGDGTs的36%～41%,但在表层土壤中占32%～51%。

在所有样品中,IIIa′丰度普遍高于其异构体IIIa。而对于五甲基化brGDGTs,在湖泊表层沉积物中,IIa′含量(分别为34%)高于IIa(5%),但土壤中IIa′ (20%)低于IIa (23%)。所有样品中其他化合物相对丰度较低,多数样品无法检测到IIIc、IIIc′。综上所述,表层沉积物和土壤brGDGTs分布基本相似,主要以Ⅱ型为主,然而,IIa和IIa′在沉积物和土壤中的优势地位不同。

2.3 BrGDGTs相关指标

图5 湖泊表层沉积物和土壤中和CBT′、MBT′ 和CBT指标交叉图Fig.5 Cross plots of and CBT′, MBT′ and CBT for surface sediments and soils

3 讨论

3.1 察汗淖尔湖泊沉积物及其周围土壤GDGTs的来源

GDGT-0是isoGDGTs中来源最广泛的化合物,它可由奇古菌(Thaumarchaeota)、泉古菌(Crenarchaeota)和广古菌(Euryarchaeota)等古菌群落合成[45],其中产甲烷菌是GDGT-0的主要贡献源[46],而crenarchaeol及其异构体是氨氧化奇古菌产生的特征标志物[47],因此GDGT-0/crenarchaeol比值可判断isoGDGTs的主要来源[46]。结果显示,所有样品GDGT-0/crenarchaeol比值在0.2～2之间,说明该地区奇古菌占主要地位。目前已培育出两组奇古菌,一种是在海洋水柱中占主导的GroupⅠ.1a[48-50];另一种是在土壤中可产生更多crenarchaeol′的GroupⅠ.1b[51-52],故可利用crenarchaeol/crenarchaeol′比值区分环境中奇古菌类型。当crenarchaeol/crenarchaeol′比值>25时,奇古菌种类主要是GroupⅠ.1a型,如果比值<25,则为GroupⅠ.1b型奇古菌[6]。察汗淖尔所有样品crenarchaeol/crenarchaeol′比值<25,表明研究区内古菌主要为GroupⅠ.1b型奇古菌。

图6 察汗淖尔样品四甲基化、五甲基化和六甲基化brGDGTs丰度三元图,并与中国北方土壤和青海湖等表层沉积物数据比较[42,55]Fig.6 Ternary diagram of the fractional abundances of tetramethylated, pentamethylated, and hexamethylated brGDGTs from the different sample sets[42,55]

3.2 细菌支链GDGTs分布的影响因素

多项研究表明pH值、含水量是影响湖泊沉积物brGDGTs分布的重要因素[8,34,36-37,47,56-60],可利用RDA分析检验这些环境参数对察汗淖尔brGDGTs分布的贡献。结果如图7所示,前两个RDA轴记录了brGDGTs分布受土壤含水量影响较小(P=0.06),主要由土壤pH值控制(P=0.02)。土壤pH值与IIIa′、IIa、Ia丰度呈正相关,与IIIa、IIa′呈负相关,造成了CBT与土壤pH值之间的相关性较弱(R=0.1,图7)。这一现象也与其他研究相符,如青藏高原的表层土壤碱性越强,pH值与CBT的相关关系越趋于平缓[61];对于中国北方的碱性土壤(pH>8),CBT与pH值的正相关性不显著[42]。本文所有样品pH在7.9～9.2之间变化,碱性较强,这可能是导致pH值与CBT不具有相关性的原因之一。

图7 BrGDGTs冗余分析(RDA)图,显示了现有环境因子(红色箭头)与brGDGTs化合物及相关指数(蓝色箭头)之间的关系Fig.7 Redundancy analysis (RDA) showing relationships between existing environment variables (red arrows) and compounds and indexes of brGDGTs (blue arrows)

3.3 基于土壤MBT′/CBT指数重建温度和土壤pH值

察汗淖尔湖泊表层沉积物及其周围土壤brGDGTs与北方表层土壤分布较为一致,因此在重建温度和pH值时使用基于土壤的校准公式更合适 (图8)。考虑到IIa和IIa′在湖泊沉积物和土壤中优势地位不同,需将两者的温度和pH值分别进行计算。在湖泊沉积物中,当使用原始全球土壤校准公式 (公式(3))时,重建的大气温度范围在1.0～4.9℃,平均为2.7℃,接近年均温度4.0℃;重建的土壤pH值(公式(4))为5.8～6.8,平均是6.3,低于平均测量pH值(8.9)。使用修订后的全球土壤校准公式(公式(7))的计算结果介于9.8～17.5℃,平均温度为13.4℃,同样高于预期值;计算的pH值(公式(8))范围在7.7～8.2之间,平均为7.9,与实际测量值相比偏低。当应用干旱地区土壤校准时(公式(9)),重建温度为7.5～11.9℃,平均是9.8℃,显然高于观测温度。运用北方区域温度校准公式(公式(10))得到的结果为 4.1～12.1℃,平均温度是7.8℃,略低于器测温度(4℃);重建的土壤pH值(公式(11))在7.5～8.0之间,均低于测量值。

图8 基于不同土壤校准公式计算的大气温度(a)和土壤pH误差值(b)(ΔpH是CBT-pH计算值与测量值的差值)Fig.8 Comparison of reconstructed temperature based on different soil calibrations (a) and the differences(ΔpH) between CBT-inferred pH and the measured pH (b)

而对于土壤,使用原始MBT′/CBT校准(公式(3))的计算结果为-2.1～9.1℃,平均为3.4℃,略低于器测记录温度(4℃);重建的土壤pH值(公式(4))为5.7～7.4,均低于平均测量值(pH=8.5)。应用修订后的土壤校准公式(公式(7))重建的温度为-3.8～19.2℃,平均是5.8℃;计算的土壤pH值(公式(8))范围在6.2～8.2之间,普遍低于实际测量值。利用Yang等建立的干旱地区校准公式(公式(9))得到的结果是7.5～12.0℃,平均温度为9.8℃,远高于年均温度4℃[41]。北方区域校准公式(公式(10))重建的温度范围较大,介于-10.0～14.0℃之间,平均为-0.1℃;重建的土壤pH值(公式(11))为6.0～8.1,与测量值相比偏低。

综上所述,使用全球MBT′/CBT土壤校准[9]得到的结果较为接近记录的年均温度,但仍存在一定偏差,区域校准结果并不比全球土壤校准更精确。值得注意的是,使用针对干旱环境或消除了5-Me和6-Me异构体影响的最新校准获得的温度结果偏高,这与其他干旱、寒冷地区情况一致[58,62],可能是由于在这种特殊的寒冷地区,brGDGTs生产细菌倾向于在夏季活跃[28,63]。

为了进一步探究重建温度产生偏移的原因,本研究利用CI指数评估细菌群落对brGDGTs分布的影响[44]。De Jonge等认为土壤brGDGTs分布随微生物群落变化而变化,并依据冷簇和暖簇细菌群落脂质将全球土壤数据划分为暖集群(CI>0.64)和冷集群(CI<0.64)[44]。经计算,我们发现察汗淖尔湖泊表层沉积物和土壤细菌群落基本与全球干旱土壤一致,均属于冷簇 (图9)。以往研究表明,多数冷簇土壤细菌brGDGTs分布与温度相关性较弱或不显著[44],而主要受土壤pH值控制[42,58],这与察汗淖尔结果一致。显然,细菌群落组成和土壤pH值可能引起上述重建结果出现偏差。

图9 察汗淖尔和全球干旱地区土壤的群落指数(CI)箱线图Fig.9 Boxplot of the community index (CI) values for the sediments and soils from Lake Chahannaoer and global arid regions

4 结论

本研究对察汗淖尔现代表层样品进行了GDGTs分析,结果显示湖泊沉积物及其周围土壤与北方表层土壤GDGTs分布相似,isoGDGTs以crenarchaeol 和GDGT-0为主,brGDGTs中五甲基化brGDGTs(IIa、IIa′)丰度最高,但IIa和IIa′在沉积物和土壤中优势地位不同。GDGT-0/crenarchaeol比值介于0～2之间,crenarchaeol与其异构体比值小于25,说明isoGDGTs主要由GroupⅠ.1b型奇古菌产生。土壤pH值是影响brGDGTs分布的主要环境因素,CI指数显示研究区产brGDGTs细菌属于冷簇,与温度相关性较弱或不显著,这些造成了应用多种温度校准公式的计算结果与器测记录不一致,重建的土壤pH值偏低。因此,在这种特殊的寒冷干旱环境中,GDGTs分布较为复杂,需谨慎选择校正公式以获得准确的气候信息。

致谢:感谢中国科学院地球环境研究所对实验工作提供的支持。