崔树昆，蒋诗威，刘孝艳，涂路遥，周 鑫

(1：中国科学技术大学地球和空间科学学院，合肥 230026) (2：伯尔尼大学奥斯切尔中心气候变化研究所和地理学院，伯尔尼 3012)

小冰期(Little Ice Age，简称LIA)是全新世北大西洋冰筏事件最后一次发生的时期，也是距现代最近的全球性典型寒冷期，时间约为1400-1850 AD[1-3]. 小冰期期间北半球整体气温比器测记录的平均气温低约0.5℃，对人类社会的发展产生了重要的影响，是当前古气候研究的热点时段[4-5].

基于不同区域的地质记录，大量研究发现中国东部季风区小冰期内部存在明显的次一级干湿波动，但重建的降水记录和干湿变化有很大的区域差异. 在我国东北地区，呼伦湖孢粉数据[6]和镜泊湖沉积物记录[7]显示小冰期气候偏湿润；而小龙湾湖泊沉积物多指标分析[8]显示在1360-1450 AD、1590-1670 AD以及过去的150 a中存在显著的干旱期；内蒙古达里湖[9]在小冰期期间也出现低水位，说明季风降水较少，气候偏干旱；阿尔山天池[10]植被演化记录反映的近500 a气候变化显示1524-1680 AD、1743-1780 AD和1895-1945 AD 3个时期偏干旱，1680-1743 AD和1780-1895 AD则偏湿润. 在华北地区，北京石花洞石笋δ13C和δ18O记录[11]显示小冰期整体较中世纪暖期偏干旱，但其内部可分为两个时期，即1450-1750 AD气候相对干旱，1750-1900 AD气候偏湿润；白洋淀湖泊沉积物研究[12]显示1340-1580 AD气候偏干旱，1580-1800 AD气候干旱加剧，1800-1920 AD气候逐渐转向湿润. 在江淮地区，Cui等[13]基于长江中游地区孢粉记录提出小冰期可分为早期(气候轻度湿润)、中期(气候干旱)和晚期(气候湿润)3个阶段；江苏高邮湖沉积物记录[14]显示1600-1900 AD气候显著干旱；神农架大九湖泥炭记录[15]则显示整个小冰期气候相对湿润.

上述研究结果显示重建的中国东部季风区小冰期干湿状况区域差异较大，并且小冰期内部干湿波动的空间特征仍不清楚. 以往研究主要集中在中国北方，小冰期内部中国南方的地质记录较少，影响了对该区域小冰期期间季风降水演变的研究. 雁荡山地处中国南方的浙江省东南部，前人对该地区古气候的研究较少. 因此本文对雁荡山雁湖沉积物进行多指标分析，建立该区域水文变化历史，并与东亚季风区其他区域的指标序列进行对比，探讨小冰期期间季风降水的时空变化特征及其可能的驱动机制.

1 研究位置与方法

雁荡山位于浙江省乐清市东北境内，是全球典型性的中生代晚期复活破火山，处于环太平洋亚洲大陆边缘构造岩浆带中的中国东南沿海中生代火山岩带[16]，属于中低山丘陵区，地势西高东低，西部为低山丘陵，东部与乐清湾相接，为海积平原. 山脉多呈东北-西南向展布，平均海拔为500～600 m. 该地区属亚热带季风气候，雨量充沛，气候温暖湿润，水热资源丰富，年平均气温13.5℃，年平均降雨量为1935.6 mm，平均相对湿度77%. 植被类型主要为亚热带常绿阔叶林、常绿落叶阔叶混交林和针阔叶混交林[17]. 本研究具体采样地点雁湖是一处山顶湖盆区域，海拔为896 m，湖盆面积较小，长度约为70.7 m，宽度约为40 m. 研究组采样时期为枯水期，湖泊已干涸，据野外分析，当湖盆重新汇聚降水形成湖泊时，水位也较浅. 湖区流域面积较小，无河流汇入，其水体主要来源于径流作用和降水汇聚，周围植被较丰富，但湖区无水生植被生长.

YH沉积柱于2016年在浙江省雁荡山国家地质公园雁湖岗(28°21′25.93″N，121°0′58.58″E)钻取(图1)，长77 cm，在实验室对岩芯进行了对剖、照相与描述，按0.5 cm间距连续分割岩芯，共获得样品147个，样品经冷冻干燥处理后保存.

图1 采样点位置(左图黑色圆点为本研究对比位置点，分别为： 1湖光岩玛珥湖；2仙山；3天目山；4华北平原；5公海)Fig.1 Site location (The black dots in the left map are sites used for comparison in this study: 1 Huguangyan Maar Lake;2 Xianshan Mountain ; 3 Tianmu Mountain; 4 North China Plain; 5 Lake Gonghai)

选取6个样品送往美国Beta实验室进行AMS14C测年，测年材料包括全岩有机碳和炭屑. 测年结果使用Calib7.0程序的IntCal13数据库[18]进行校正，并使用2σ校正年龄(表1). 总有机碳(TOC)测试在中国科学技术大学理化科学实验中心进行，具体步骤为：(1) 使用球形磨样机对样品进行研磨；(2) 取1.5 g左右的样品，加入稀盐酸多次搅拌，去除样品中的无机碳酸盐；(3) 用中性去离子水水洗并且多次离心至中性(pH = 7)，置于恒温箱烘干后保存；(4) 称取一定量的被测样品，用锡舟紧密包裹后送入氧化炉中，由自动进样器送入Elementary vario EL cube型元素分析仪，在960℃的氧化条件下测定TOC的百分含量. 色度测试在中国科学技术大学进行，将样品装入石英杯使用CR-700R测色仪测试. 粒度分析在泰山学院土壤与环境变化实验室进行，具体方法如下：取少量样品加10～20 mL 30%的H2O2，摇匀后在150℃条件下加热以除去有机质，之后加入10 mL 10%的盐酸除去碳酸盐，静置48 h后加入10 mL 10%的分散剂(六偏磷酸钠)，超声10 min后，用Mastersizer 2000型激光粒度仪进行粒度测试. 仪器测量范围为0.02～2000 μm，粒级分辨率为0.01Φ，重复测试的相对误差<2%.

2 结果

2.1 年代

YH钻孔6个样品的AMS14C测年数据如表1所示. 由于测年误差，个别数据存在倒转现象，但整体上符合沉积规律. 为了更好地处理年代-深度关系，我们基于贝叶斯公式的“WinBacon2.2”年代模型[19]，使用R软件[20]为不同的年代控制点选择合适的函数，从而建立了YH钻孔的年代框架(图2).

年代-深度模型(图2)显示，YH钻孔在41 cm以下沉积速率较低，41～57 cm之间沉积速率为0.011 cm/a，57～77 cm为0.012 cm/a. 41 cm以上沉积速率较高，0～13.0 cm之间沉积速率为0.058 cm/a，13.0～23.5 cm为0.058 cm/a，23.5～35.0 cm为0.15 cm/a，35.0～37.5 cm为0.093 cm/a，37.5～41.0 cm为0.026 cm/a.

表1 YH钻孔AMS14C测年结果

图2 YH钻孔年代-深度模型Fig.2 The age-depth model of the YH Core

2.2 岩性和代用指标

由YH钻孔沉积剖面(图3a)可以看出，整个沉积剖面有明显的分层现象，上部(41 cm以上)颜色较深，有机质含量高；底部(41 cm以下)为浅黄色，颗粒偏粗，有机质含量低. YH钻孔的各项环境指标数据分布如图3b所示，L*为亮度，在明亮(此时L*=100)与黑暗(此时L*=0)之间变化，a*值表示颜色从绿色(-a*)到红色(+a*)之间变化，a*/L*为通过亮度校正后的红度值. 钻孔粒度数据按照黏土(<4 μm)、粉砂(4～64 μm)和砂(>64 μm)分类. 结合YH钻孔岩性和各项环境指标数据，整个沉积序列可分为3段:

(1) 0～23 cm岩芯为浅褐色，TOC含量平均值为1.83%，有机质含量较高，粒度显示该段岩芯黏土含量较低，粉砂含量较高，中值粒径在8.79～12.3 μm之间，平均值为9.47 μm，沉积物颗粒中等偏细；(2) 23～41 cm岩芯为深褐色，TOC含量平均值为2.55%，有机质含量为最高，该段沉积黏土含量较高，粉砂含量较低，中值粒径在7.02～8.93 μm之间，平均值为7.72 μm，说明沉积物颗粒偏细；(3) 41～77 cm岩芯为黄褐色，TOC含量平均值为0.67%，有机质含量为最低，黏土和粉砂含量均较低，而砂含量较高，中值粒径在8.16～19.3 μm之间，平均值为10.22 μm，沉积物颗粒偏粗.

图3 YH钻孔岩性图(a)和指标序列(b)(▲为AMS14C测年点，红色虚线为岩性分界)Fig.3 The lithology (a) and proxies (b) of YH Core (Triangles show AMS14C ages and the red dotted lines show distinctive lithological changes)

由于：1) 本研究所采取的6个AMS14C测年样品多数集中在顶部至41.0 cm处，底部测年的数据较少；2) 底部沉积物TOC含量较低，可能导致全岩有机碳测试年龄存在误差，另外据钻孔剖面岩性初步判断(1)、(2)段岩芯为湖相沉积，而(3)段岩芯可能为早期岩石的风化产物而非湖泊沉积相，因此本文将着重探讨41.0 cm(大约1400 AD)以上沉积序列所反映的气候变化信息.

3 讨论

3.1 1400 AD以来研究区环境指标的变化特征与湖区古气候重建

湖泊沉积物的TOC含量常被用于反映沉积时期生物量的多少以及氧化还原环境条件[21]. 通常在温暖湿润气候时期，湖泊内水生生物及藻类发育，且流域内植被茂盛，沉积物中TOC含量较高[22-23]. 研究区雁湖为山顶面积较小、水深较浅的湖泊，同时雁湖属于亚热带季风区且地处较高海拔位置，相比于温度的变化，该地区的降水量在季节-年代际范围都具有较大波动，影响湖区及周围流域的生物量，因此认为降水量可能是控制雁湖沉积物中TOC变化的主要因素[24]. 当然，对于浅湖来说，其他因素也会影响沉积物中TOC含量，因此本文结合了TOC、粒度和色度进行多指标对比分析.

湖泊沉积物颗粒的粗细变化能够反映水体的扩张和收缩，从而指示沉积环境的干湿变化[25-26]. 雁湖属于封闭湖泊且水位波动较大，其沉积物来源单一，沉积物颗粒大小与水位高低密切相关. 在气候干旱期，湖水退却，采样点距离湖岸较近，粗颗粒物质易于到达；反之，在气候湿润期，湖泊扩张，采样点距离湖岸较远，粗颗粒物质难以到达，该位置的沉积颗粒较细[27-28]. 因此，雁湖沉积物颗粒的粗细变化能够反映湖区的干湿变迁，即：细颗粒沉积物标志着湿润气候，粗颗粒沉积物反映干旱气候.

湖泊沉积物的色度也经常用于气候变化研究，其常用指标包括亮度L*和红度a*. 亮度L*主要受控于沉积物中碳酸盐和有机质含量变化，碳酸盐含量越高，则沉积物的亮度值就越大，而有机质含量的增高则使沉积物变暗[29-30]. 红度a*与沉积物中氧化铁含量有关[31-32]，沉积物有机质含量和湖泊水位高低会影响氧化铁等致色矿物的形成，即高水位、有机质含量丰富的偏还原环境抑制了致色矿物形成，而低水位、有机质含量低的偏氧化环境促进了矿物形成.

图4 1400 AD以来雁湖沉积记录的气候 变化序列(虚线为文献记录时间)Fig.4 Time-series of proxies in sediments of Lake Yanhu since 1400 AD (The dotted lines correspond to the ages in historical documents)

雁湖地区1400 AD以来气候变化可分为4个阶段(图4). 阶段Ⅰ: 1400-1600 AD期间，沉积物TOC含量逐渐增加，说明湖区及周围流域生物量较高，进而反映了较多的季风降水. 这段时期L*逐渐降低反映了有机质含量高使得沉积物偏暗，a*/L*整体偏低表明沉积环境偏还原条件，抑制了赤铁矿等致色矿物的形成. 另外中值粒径在该时期也逐渐减小，沉积物颗粒偏细. 各项指标均说明，1400-1600 AD期间，研究区降水量较丰富，气候偏湿润，湖区集水发育，水位较高；阶段Ⅱ: 1600-1650 AD期间，TOC含量大幅降低，L*和a*/L*相应升高，中值粒径显示沉积物颗粒变粗，说明这段时期研究区降水量减少，气候转向干旱，湖区水位下降；阶段Ⅲ: 1650-1750 AD期间，TOC含量出现回升，L*和a*/L*降低，但该时期中值粒径并无显著变化，仍显示沉积物颗粒偏粗，说明该时期研究区降水有小幅度增加，但湖区水位仍然较低；阶段Ⅳ: 1750-2000 AD期间，TOC含量逐渐降低且之后一直较低，L*和a*/L*升高且保持在较高水平，中值粒径偏大，各项指标均说明研究区气候一直处于干旱状态，降水偏少.

3.2 与历史文献记录的对比

文献资料中记录了大量的古气候信息，雁湖地区的水文状况在多篇历史文献中也有反映. 1590s，明代陈仁锡在其所著《潜确类书》中记载到：“雁荡山在东清县山顶有一湖，方可十里，水常不涸，春雁归时都宿此，故名”，说明此时期雁湖湖泊面积较大，水位较高，对应于本研究得出的1550-1600 AD研究区降水量较多，气候偏湿润的阶段. 明代地理学家徐霞客分别在1613 AD和1632 AD(图4中黑色虚线标注)两次游历雁荡山并且留下两篇游记，第1篇(1613 AD)记载到：“湖中草满，已成芜田”，第2篇(1632 AD)记载到：“洼中积水成芜草之地，青青弥望满眼，视野所及处，所称雁湖也”，说明该时期雁湖水位较低，在年代误差范围内，对应于本研究中1600-1650 AD湖区降水量减少，气候偏干旱的阶段. 另外，岩性显示钻孔约23 cm(约1613 AD)以上沉积物颜色逐渐变浅，对照YH钻孔TOC含量降低和粒度变粗说明了湖泊沉积物中有机质含量减少，气候较为干旱的状态. 综上，本研究各项环境指标分析结果与历史文献记录的气候变化特征较吻合，说明雁湖沉积物小冰期以来的沉积序列年代框架和环境指标所反映的气候变化信息较准确.

图5 南北地区干湿记录以及与太平洋年代 际振荡的对比(橙色阴影代表小冰期前期，蓝色 阴影代表小冰期后期)：(a) YH钻孔TOC含量； (b) 湖光岩玛珥湖TOC含量；(c) 福建仙山泥炭 α-纤维素δ13C；(d) 浙江天目山泥炭Rb/Sr； (e) 公海磁化参数S-300mT；(f) 华北平原干湿指数； (g) 孢粉记录重建的北方年平均降水量； (h) 太平洋年代际振荡指数(PDO指数)Fig.5 Comparison of records in eastern China with the Pacific Decadal Oscillation (orange shadows represent the early LIA, blue shadows represent the late LIA): (a) TOC content from YH Core; (b) TOC content from Huguangyan Maar Lake; (c) the percentage of cellulose δ13C from the Xianshan Peat in Fujian Province; (d) the Rb/Sr ratio of the Tianmu Mountain in Zhejiang Province; (e) the variations of S-300mT in Lake Gonghai; (f) dry-wet index in North China Plain; (g) the pollen-based PANN anomaly in northern China; (h) Pacific Decadal Oscillation index (PDO index)

3.3 季风降水的区域对比

为了进一步确定小冰期期间东亚季风区南部的降水变化特征，本文将雁湖沉积物的指标与被用于反映降水变化的湖光岩玛珥湖沉积物TOC含量[33]、福建仙山泥炭α-纤维素δ13C记录[34]和浙江天目山千亩田泥炭Rb/Sr值[35]进行对比. 分析发现(图5a～d)：在年代误差范围内，上述地区的降水变化趋势较为一致, 总体上在小冰期前期(1400-1600 AD)气候偏湿润，而小冰期后期(1600-1850 AD)气候偏干旱. 天目山、仙山、湖光岩玛珥湖和雁荡山同处于中国东南季风区，因此4个地区的气候变化序列可能反映了中国东南地区的干湿变化趋势.

为深入研究小冰期内部中国东部季风区南北降水时空变化特征，我们将上述南方古气候记录与北方地区(图5e～g)降水重建记录(包括公海沉积物磁性特征记录[36]、基于历史文献重建的华北干湿指数[37]和孢粉记录重建的北方年平均降水量[38])进行对比. 结果显示：与南方地区记录反映的干湿状态所不同，北方古气候记录显示在小冰期前期气候偏干旱，而在小冰期后期气候偏湿润. 上述对比说明在百年时间尺度上，南北方地区降水气候变化趋势相反，小冰期期间中国东部季风区降水存在“南涝北旱”和“北涝南旱”的空间变化特征.

3.4 驱动机制分析

太平洋年代际振荡(Pacific Decadal Oscillation，PDO)是北太平洋海温多年代际变化的主要模式[39-40]. 在PDO暖相位时期，北太平洋中西部海温偏冷，而东部海温则高于往年平均温度，PDO冷相位时期情况则相反[41-43]. 现代气象观测数据显示我国东部夏季降水的空间格局和时间变异性与太平洋海温异常具有较好的相关性[44-45]. 例如，一些研究认为PDO相位的移动与华北降水的年代际变化密切相关，在PDO暖期，降水低于正常水平，华北地区干旱加剧，反之则情况相反[46-48]. Zhu等发现：与1979-1999年相比，2000-2008年黄淮流域夏季降水增加，长江流域降水减少是由1999年前后PDO由暖向冷的相变引起的[49]. Yang等利用1951-2015年160个台站月降水数据研究发现：在太平洋海温模式暖期(冷期)，中国东部地区表现为“南涝北旱”(“北涝南旱”)[50].

本文将中国东部南北地区重建的1400 AD以来的干湿记录与PDO指数[51](图5h)进行对比，分析发现：小冰期前期(1400-1600 AD)，PDO处于暖相位，中国南方地区降水量相对较高，气候湿润，北方地区降水量则相对较低，气候干旱；小冰期后期(1600-1850 AD)，PDO处于冷相位，中国南方地区降水量减少，气候偏干旱，而北方地区则呈现相反趋势，其降水量增加，气候偏湿润.

关于中国东部旱涝分布变化与PDO的内在联系，有研究[52-53]认为太平洋海温异常可以通过影响西太平洋副热带高压(西太副高)的强度和位置进而改变中国东部降水空间格局. 一般而言，西太副高会在每年的6月中旬和7月中下旬出现向北迁移的过程，伴随着西太副高的迁移，中国东部地区雨带地理位置也会发生相应变化[54,55]. 在PDO暖相位时期，西太副高增强并向西南扩展，使得雨带在中国南方地区停留较长时间并且导致更多降水，而北方地区降水则减少；相反在PDO冷相位时期，西太副高减弱且向东北退缩，导致中国南方地区提前结束了梅雨期，而北方地区则迎来更多降水[56]. 本研究中地质记录对比也显示小冰期前期受PDO暖事件影响，中国东部地区降水表现为“南涝北旱”；小冰期后期受PDO冷事件影响，东部地区降水表现为“北涝南旱”. 综上分析，小冰期以来，PDO对中国东部季风区南-北“偶极型”空间降水特征具有重要影响作用.

4 结论

本研究对YH钻孔进行了多指标对比分析，并与中国东部南北地区季风降水指标序列以及PDO指数进行对比分析，初步得到以下结论：

1)1400 AD以来研究区的气候变化大致可分为4个阶段：a)1400-1600 AD，气候偏湿润，湖区集水发育；b)1600-1650 AD，降水量减少，湖区水位下降；c)1650-1750 AD，研究区气候出现回转，降水量有小幅增加，但湖区水位仍较低；d)1750-2000 AD，研究区气候一直处于偏干旱状态. 其中1500-1650 AD期间研究区的干湿变化与历史文献记录能够很好地对应.

2)在百年时间尺度上，雁湖沉积记录与天目山泥炭、福建仙山泥炭以及湖光岩玛珥湖沉积记录的干湿变化同步，指示了中国东南地区气候波动的一致性.

3)小冰期以来中国东部季风区南北降水记录反映的气候变化状况大体上呈反相位关系，说明在百年时间尺度上，东部地区降水存在着“南涝北旱”和“北涝南旱”的“偶极型”特征. 这种空间降水特征与PDO具有密切相关性：小冰期前期(1400-1600 AD)受PDO暖相位事件影响，中国东部地区表现为“南涝北旱”，而小冰期后期(1600-1850 AD)受PDO冷相位事件影响，表现为“北涝南旱”的特征.

致谢：感谢黄开建在采样过程中的帮助.