贵州石笋记录的中晚全新世东亚夏季风变化

张振球 张伟宏 刘树双 翟秀敏 邵庆丰

关键词 石笋；氧同位素；全新世；4.2 kyr B.P.事件

0 引言

全新世发生了一系列百年尺度的气候突变事件[1?4]，表明全新世间冰期的气候并不稳定。其中，发生在4.2 kyr B.P.时期的气候突变过程被认为具有广泛的区域影响，被称为4.2 kyr B.P.事件[5]。近来，该事件也被作为中晚全新世的分界线[6]。因此，研究此事件的结构特征、区域气候响应不仅有助于理解季风突变事件的成因，而且对进一步完善气候地层学有着重要意义。

在东亚季风区，高分辨率的气候记录呈现了季风对百年尺度突变事件的响应[7?8]。Zhang et al.[9]通过综合分析我国季风区高精度定年和高分辨率记录，指出4.2 kyr B.P.事件具有时空差异性。在空间分布方面，我国季风区的北部和西南部区域在该时段整体表现为干旱状态，而在东南部区域则表现为湿润状态。在持续时间方面，北部和西南部记录的干旱4.2 kyr B.P.事件持续约500年（4.4～3.9 kyr B.P.），而东南部的湿润4.2 kyr B.P. 事件持续约300 年（4.2～3.9 kyr B.P.）[9]。近年来，位于黄土高原高分辨率的石笋记录表明，4.2 kyr B.P.事件的内部结构不同于其他发生在全新世的百年尺度干旱事件[10]。4.2 kyr B.P.事件期间，季风降水整体偏低，并呈现逐渐减少的趋势，内部发生有两次明显的强降水过程[10]。虽然重建记录表明我国的西南与北方地区在4.2 kyr B.P.事件都表现为干旱气候特征[9]，但是4.2 kyr B.P.事件的内部结构是否与北方一致，目前尚不明确。研究4.2 kyr B.P.事件内部结构的区域特征以及区域之间的异同，有助于明确4.2 kyr B.P.事件在不同区域是否存在持续时间方面的差异。此外，还有助于进一步理解4.2 kyr B.P.事件的发生原因。

贵州七星洞已经提供了晚更新世以来的气候环境变化历史记录[11?16]。这些成果展现了七星洞记录对不同时间尺度事件的敏感性，不仅标定了冰期Dansgaard-Oeschger和Heinrich事件的发生和持续时间[11?13]，还重现了全新世以来季风水循环的过程[14?16]。为了更详细地研究中晚全新世季风变化过程，本文提供了过去6 380～2 700 yr B.P.间高分辨率的七星洞石笋δ18O 记录，此记录与同区域董哥洞记录[17]在4.2 kyr B.P.事件上呈相同的振荡特征，并与北方石笋记录相似。此七星洞记录显示，在4.2 kyr B.P.事件内部，也有两次强降水过程，分别发生在4.15 kyr B.P.和3.93 kyr B.P.，但前者更为显著。4.2 kyr B.P.事件内部的这两次强季风过程成为对比不同区域4.2 kyr B.P.事件内部振荡的有效媒介。

1 材料与方法

编号为QX3的石笋（下文简称为QX3）发育于贵州省都匀市凯口镇凯酉村的七星洞（25°59' N，107°16' E，海拔1 020 m），距离董哥洞约100 km[17?20]。该洞长383 m，宽10～15 m，高约10～15 m。洞内年均温为16.1 ℃，相对湿度为95%。此洞穴发育于下石炭统摆佐组（C1b）灰色中厚层白云岩、白云质灰岩中，上覆基岩最大厚度约90 m。上覆土壤层相对较薄，普遍小于10 cm，局部基岩裸露。洞顶植被以灌木为主，草本植物茂盛，总体覆盖度较好。该地为亚热带季风气候，年均降水量约1 176 mm，77.5%的降水分布在4—9月[21]。现代洞穴监测结果表明，该洞洞穴滴水对外界降水响应敏感，洞穴滴水同位素具有明显的季节变化信号[22]。

QX3捡拾于洞穴中部，被发现时已碎裂为多段。可拼接起来的部分呈柱状，高780 mm，直径60～90 mm。QX3被切开抛光后，在剖面上未发现明显的沉积间断，但存在矿物相的差异（图1）。在225～255 mm和318～387 mm段，石笋沿生长轴方向呈疏松放射状，为白色不透明矿物，扫描电子显微镜下观察发现晶体呈纤维状（图1c），X射线衍射测试显示其为文石，具有高238U特征[23]。其余部分为白色透明致密矿物，对应的238U含量则相对较低且稳定，扫描电子显微镜下观察晶体均呈块状菱面体，未发现球粒状或纤维状结构，推测为纯方解石且未发生重结晶过程（图1b）。

为重建过去气候环境变化，对石笋进行稳定同位素和年龄测试。为获取稳定氧同位素数据，采用头部直径为0.3 mm 的钻头沿着生长轴钻取50～100 μg的粉末，取样间隔为1 mm，测试仪器为带有碳酸盐自动进样装置的MAT253同位素质谱仪，共测有779组氧同位素，分析误差优于0.05‰（±1σ），结果以Vienna Pee Dee Belemnite（VPDB）标准给出。对于年龄样品，采用头部直径为0.9 mm的钻头沿着石笋生长层钻取100～200 mg 的粉末，共测得14 件230Th 年龄，测试仪器为Neptune多接收电感耦合等离子质谱仪（MC-ICP-MS）。U和Th元素的分离纯化过程以及Neptune 测试分析方法分别见文献Shao et al.[24]和Cheng et al.[25]。以上相关测试工作均在南京师范大学同位素质谱实验室完成。

2 结果与讨论

2.1 年龄与氧同位素记录

表1呈现了QX3的14个230Th年龄结果，显示从石笋底部到顶部无倒序现象，符合地层层序规律。样品年代误差在10～75 yr之间，但在225～255 mm和318～387 mm 段的疏松放射状层段，误差为11～15yr，相对较小。在这两个层段的238U 含量较高，达7.7～11.9 μg/g，约是其余层段238U含量（0.3～0.4 μg/g）的20倍。在同一支石笋中238U含量突变，可能与组成石笋的矿物有关。一般来说，由于文石与方解石的结构不同，文石具有更强的离子位来捕获更大的阳离子，如易溶于水的铀酰根离子（UO22 +）[23，26]，因此U含量更高。在整个石笋沉积序列中，从底到顶230Th/238U的值逐渐增大，意味着铀元素衰变体系的封闭性。疏松的文石若发生U 元素的流失，230Th/238U 在此层段的值要较随后沉积的方解石值大。基于以上初步分析，QX3沉积时未发生U元素流失，U/Th年龄可靠。依据线性内插方法建立QX3氧同位素的时间标尺（图2a黑线），该时标序列得到MOD-AGE[27]模型模拟结果的证实（即内插时标序列均处于MOD-AGE 所模拟的置信区域内）。采用内插方法所建时标显示石笋QX3记录的持续时间为6 380～2 700 yr B.P.。

QX3 δ18O 记录的变化范围为-8.55‰～-5.48‰，平均值为-7.02‰，变化幅度为3.06‰，平均分辨率为4.7 yr。QX3 δ18O记录在6 380～2 700 yr B.P.时段呈逐步增加的趋势（图2b），并在长期趋势上叠加了一系列十年—百年尺度的次一级振荡，其中超出整个记录一倍标准偏差（1σ）的振荡主要发生在4 548～3 715 yr B.P.时段（图2b）。此时段对应于4.2 kyr B.P.事件。因此，QX3 δ18O记录对解析4.2 kyr B.P.事件的结构特征具有重要意义。

2.2 QX3石笋记录的我国西南气候变化

洞穴石笋同位素分馏平衡是其进行古气候环境重建的基本前提。除了Hendy检验[28]，同一或不同洞穴记录的重复性检验也是常用的验证方法。本文QX3 δ18O记录与此洞另一支已发表QX1石笋[29]在重叠时段有着良好的对应关系（图3）。尽管两记录在局部有微小差异，但在百年—数十年尺度上呈现相似的振荡特征。两记录良好的重现性一方面说明QX3沉积时未发生显著的动力分馏，处于同位素平衡分馏状态；另一方面也暗示QX3在沉积过程中受岩溶渗流系统、蒸发作用等因素影响较小，QX3 δ18O记录可以用作过去气候环境变化的代用指标。中国石笋δ18O指标在大气水循环和全球对比领域有突出优势，但有些研究认为其指示的气候意义仍不明确，存在较大争议[30?34]。Cheng et al.[35]和Zhang et al.[36]通过综合分析我国石笋记录，认为石笋δ18O记录在轨道和千年尺度指示了季风强度的平均变化状态。最近，石笋δ18O记录与其他反映降水量指标记录（如黄土、湖泊等）之间的“矛盾”逐渐得到澄清，即它们可能分别响应于东亚季风系统的不同方面，石笋δ18O主要响应于大尺度季风环流和季风雨带的移动，可指示东亚夏季风总体强度变化[37?38]。但百年尺度各地的δ18O记录则存在差异[15]，甚至呈反向的变化[39]。有些洞穴附近的器测降水与降水δ18O记录有很好的对应关系，明确了石笋δ18O 记录可指示季风降水变化[39?40]。在我国西南区域，石笋δ18O记录在百年—数十年尺度上被认为主要反映了季风降水的变化，δ18O越小，指示降水量较大，反之亦然[41]。已有的研究表明贵州七星洞洞穴滴水对大气降水响应迅速[42]，δ18O的振荡主要反映了与降雨量变化相联系的不同季风降水云团对该地区降雨贡献的剧烈变化[14]。为进一步验证这一解释，我们对比了QX3 234U/238U 初始值（234U/238U）0 与δ18O记录的异同（图2）。对比结果显示，（234U/238U）0 与δ18O记录在数百年尺度上的变化基本一致。他们的对应关系可解释为：当季风整体加强时（δ18O值减小），增加的降水导致土壤中234U优先淋滤作用减弱，土壤水中234U减少，从而导致石笋234U/238U值减小[43]。七星洞石笋δ18O与234U/238U初始值有着很好的对应关系，说明当地降水量的重要贡献。但是上游降水对石笋δ18O的影响程度，还需更深入的研究。

QX3 δ18O记录显示亚洲季风平均强度在6 380～2 700 yr B.P.逐渐减弱。该减弱过程与北半球夏季太阳辐射逐渐降低相对应，支持了太阳辐射变化在轨道尺度对季风变化的控制作用。然而，有研究发现季风在中晚全新世减弱的开始时间（即适宜期结束时间）随着纬度的增加而延迟[44]。附近的董哥洞记录的全新世季风开始减弱的时间为7.0 kyr B.P.[17]，QX3在生长期内δ18O持续偏正的变化支持了董哥洞记录的季风演化过程。太阳辐射变化调控着亚洲夏季风的整体变化，这些亚轨道及千年尺度的振荡表现为其他气候系统如北大西洋[45]、太平洋[44]扰动了太阳辐射对季风环流的控制。

叠加在季风逐渐减弱长期趋势上最显著的百年尺度振荡主要发生在4 548～3 715 yr B.P.时段。此外，较为显著的季风振荡还发生在3 485～3 230 yr B.P.和5 800～5 560 yr B.P.时段，持续时间约250 yr的季风减弱事件结构简单，均表现为季风先逐渐减弱后又增强的特征（图4 灰色阴影部分①④）。其中，5 800～5 560 yr B.P.时段的弱季风事件内部存在一次数十年尺度的次级季风增强过程，4 548～3 715 yr B.P.时段的季风振荡对应于具有广泛区域影响的4.2 kyr B.P.事件。QX3 δ18O记录显示，在4.2 kyr B.P.干旱事件内部约4.15 kyr B.P.时段有一次显著的季风强降雨事件，奠定了该弱季风事件两（干）谷夹一（湿）峰的整体结构。同时，在3.93 kyr B.P.也出现一次降水增强的过程，但其幅度和持续时间较之于4.15 kyr B.P.强降水事件则相对较弱。

前文所述的发生在225～255 mm和318～387 mm层段的文石矿物分别对应于3 856～3 731 yr B.P.和4 548～4 250 yr B.P.时段，分别位于4.2 kyr B.P.事件的双谷结构中（图4灰色阴影部分②③）。此双谷所夹的位于4.15 kyr B.P.的强降水信号（δ18O负偏）为方解石记录。此方解石为原生方解石，并非由文石转变而成的转化方解石，真实地记录了当时的季风变化信号。证据如下：1）4.2 kyr B.P.事件内部的方解石与其余层段的方解石铀含量相当，说明石笋中的文石和方解石都在石笋沉积时形成，不存在文石向方解石转化的现象。2）4.2 kyr B.P.事件内部方解石与其他部分（包括文石）的230Th年龄线性相关，未发生年龄倒转现象，符合地层层序规律。在整个石笋样品中，从底到顶，230Th/238U值逐渐增加，无乱序现象。文石一般发育在高温或者干旱的环境中[47?48]，相对来说洞内温度较为稳定，文石的发育可能与大气降水减少有关。以方解石为主体的QX3样品，在生长过程中多处于滴水长期淋滤状态，即降水较多引起持续滴水的湿润环境。而4.2 kyr B.P.事件发生时，可能因降水较少造成持续性滴水减少，滴水速率变慢，渗流水在上覆岩层中滞留时间变长，造成七星洞上覆围岩（如白云岩、白云质灰岩）中较多Mg元素被淋滤，滴水中Mg/Ca比值增加。当滴水中Mg/Ca 浓度增加到一定值时，高浓度的Mg2+可能抑制方解石的形成，促成文石产生[26，47?50]。因而，4.2 kyr B.P.事件内部的两个文石层发育，表明此时期气候较为干旱。相反，4.2 kyr B.P.事件内部方解石层的发育，则主要由于降水增多造成持续滴水的湿润环境，再次形成方解石矿物[47]。此方解石层的发育，印证了4.2 kyr B.P.事件内部存在季风增强的过程。

QX3记录的三次百年尺度季风干旱事件均得到附近董哥洞和莲花洞记录的支持[17，46]。有区别的是湖南莲花洞石笋在约3.32 kyr B.P.附近显示季风在约70 yr时间内快速减弱，并在弱季风状态下持续了约300 yr，并不同于董哥洞与七星洞记录，此差异可能源于局地气候的分异。但三者在4.2 kyr B.P.事件内部结构的一致性，再次验证了QX3揭示的4.2 kyr B.P.弱季风事件结构的可靠性。

已有研究显示，我国西南季风区与北方石笋都记录了4.2 kyr B.P. 事件发生时的干旱过程[10，17，46]。QX3 δ18O记录的4.2 kyr B.P.干旱事件内部两次季风降水增强过程，得到了同区域董哥洞记录的支持[17]。这两次季风增强事件在甘肃乌鸦洞、陕西九仙洞和祥龙洞记录中均有表现（图5）[10，44，51]，说明南北方在4.2 kyr B.P.干旱事件的背景下，都发生有两次季风降水增强波动。除QX3记录外，其他记录均显示，发生在4.15 kyr B.P.的强季风过程比3.93 kyr B.P.事件持续时间长且强度大。QX3 δ18O 记录的发生在4.15 kyr B.P.的强降水事件在整个记录中也较为突出（图5）。季风雨带的南北移动和停留时间的不同导致了南北方地区在4.2 kyr B.P.事件内的降水差异[52]。我国北方和西南地区记录的4.2 kyr B.P.事件内部结构一致的两次强降水过程说明两地在该事件发生时季风雨带的西侧可能停留于七星洞之南。或者4.2 kyr B.P.事件发生时，受北高纬反气旋影响，季风雨带整体向东退缩。

2.3 4.2 kyr B.P. 事件时亚洲季风与热带和南半球气候的联系

亚洲—澳大利亚季风系统是全球季风的组成部分。冰消期及近1 000年的时段内，千年—百年尺度的气候事件（如Heinrich 1， Bolling-Allerod， YoungerDrays， Little Ice Age等）在亚洲和澳—印季风中呈反相关系[53?55]。对比两区域石笋记录发现，这种反相关系可能也存在于4.2 kyr B.P.事件及其内部更短的时间尺度上。重现性良好的印度尼西亚石笋LR06-B1[56]和LR06-B3[53]δ18O记录显示4.2 kyr B.P.事件发生时，澳—印季风整体增强（图6a）。七星洞和董哥洞石笋记录的4.2 kyr B.P.事件内部的两次强降雨过程在澳—印季风记录中也有表现，但为季风降水减弱事件。甚至，在事件幅度上也与亚洲季风系统记录相对应并呈反相位关系。其中发生在早期的干旱事件（图6中标号为1）的幅度与亚洲季风记录一致，都强于晚期季风波动（图6中标号为2）。这种对应关系说明亚洲与澳—印季风在4.2 kyr B.P.事件时的紧密耦合关系，证实了这种反相位的耦合关系可能不受气候边界条件影响。也进一步证明亚洲与澳—印季风变化的动力联系[57]。4.2 kyr B.P.事件与其他全新世如8.2 kyr B.P.事件类似，一般都认为是发生于北大西洋地区或者由北大西洋温盐环流将此事件的信号传向全球[2，9，17，58]。近期，南北两极冰芯研究发现北大西洋温盐环流将北高纬信号传向南极将需约200 年[59]。4.2 kyr B.P.事件时，亚洲和澳—印季风在百年尺度的良好对应关系说明大气环流在传输南北半球信号中的重要作用。Denton et al.[60]指出冰消期发生时大气系统将北高纬冰消信号传向南半球，促使环南极流增强，导致冰期储存于南大洋的CO2释放到大气中，最终完成全球冰消过程。按照此理论，4.2 kyr B.P.事件发生时，北半球西风带南移，亚洲季风减弱，热带辐合带南移，此时南半球的澳—印季风增强[53]。通常，热带辐合带的南北移动可解释南北半球季风系统的反相位关系[61?62]，但其变化的根本原因是否与北大西洋温盐环流直接相关，还需要进一步研究[63]。

亚洲季风水汽主要来源于热带海洋，明确亚洲季风在4.2 kyr B.P.时的变化是否受ENSO影响，对完整理解4.2 kyr B.P.事件的发生原因具有重要意义。厄瓜多尔湖泊沉积物记录[64]显示在4.2 kyr B.P.事件时ENSO与QX3 δ18O记录在百年尺度上有很好的对应关系（图7）。当亚洲季风水循环增强时，La Nina事件的发生频率增强，反之亦然。如在4.15 kyr B.P.处的季风增强事件对应于La Nina的多发期（图7阴影部分），此时，沃克环流上升支的气流活动增强，输送到东亚的水汽也增多。Tan[65]指出ENSO可通过西太平洋副热带高压来影响季风降水在百年尺度的变化。如在La Nina态时，西太平洋副热带高压减弱东缩。这样的变化有利于热带海洋水汽向大陆输送，导致季风降水增多。但是，石笋δ18O在4.2 kyr B.P.事件的内部振荡是否源于ENSO控制下的印度洋与太平洋不同水汽比例的变化，目前还不清楚。若我国东部地区在4.2 kyr B.P.事件内部也明确存在两次百年尺度季风增强事件，则不能排除ENSO对大尺度季风环流的影响。ENSO 变化与QX3 δ18O 记录在4.2 kyr B.P.事件内部百年尺度的对应性，说明热带太平洋对东亚季风4.2 kyr B.P.事件的发生具有重要调控作用。需要指出的是，QX3 δ18O记录与ENSO变化记录仅在4.2 kyr B.P.事件时有明显的对应关系。而其他时段两者关系并不明显，这可能与不同气候系统的振荡频率及驱动因素有关。

3 结论

基于高精度230Th定年和高分辨率的贵州七星洞石笋氧同位素数据，重建了6 380～2 700 yr B.P.时段的东亚夏季风演化序列。与我国其他洞穴石笋记录一致，该石笋δ18Ｏ记录呈现逐渐偏正的长期变化趋势，表明中全新世以来东亚夏季风持续减弱，这可能响应于北半球夏季太阳辐射逐渐降低的变化。QX3δ18Ｏ记录的4 548～3 715 yr B.P.时段季风快速突变过程，对应于4.2 kyr B.P.事件。高分辨率石笋δ18Ｏ记录捕捉到4.2 kyr B.P.时期季风整体干旱的背景下发生两次显著的季风强降雨期。此两次强降雨期得到我国北方记录的进一步支持。此外，还发现亚洲季风区4.2 kyr B.P.事件的结构特征与澳—印季风区的石笋记录呈反相位的耦合关系。在东亚和澳—印季风系统中，4.2 kyr B.P事件内部结构在持续时间和振荡幅度方面的一致性，证实了在数十年—百年的更短尺度上季风系统间的成因联系。同时，4.2 kyr B.P.事件的结构特征与ENSO变化相对应，在La Nina事件发生频率升高期，东亚夏季风增强，暗示了低纬热带太平洋可能在一定程度上对东亚夏季风4.2 kyr B.P.事件的发生具有重要控制作用。