西藏东南部末次冰期早阶段冰川作用及其古气候意义

谢金明，周尚哲，许刘兵，杨太保

（1.兰州大学 资源环境学院，甘肃兰州 730000；2.华南师范大学 地理科学学院，广东广州 510631）

0 引言

西藏东南部念青唐古拉山东段，地貌特征主要表现为海拔6 000 m以上的夷平面以及深切的河谷。本区由于临近雅鲁藏布江大拐弯，南亚季风带来的暖湿气流可以伸入河谷，使之成为青藏高原上降水量最多的地区之一。丰沛的降水加之高大的地形，也使本区成为高原上现代海洋性冰川及第四纪冰川规模最大的地区之一。早在20世纪70年代，李吉均等［1］考察本区波堆藏布谷地古冰川遗迹，提出了“古乡冰期”和“白玉冰期”两次冰川作用。其中，白玉冰期（即末次冰期）分为早期和晚期两个阶段，相应的冰川沉积在波堆藏布主要支沟沟口（白玉沟、珠西沟等）保存最为完整，地貌上可划分为外冰碛垄和内冰碛垄，推测其形成年代分别对应于海洋氧同位素阶段（marine oxygen isotope stage，MIS）4和2。周尚哲等［2］应用宇宙成因核素（terrestrial in situ cosmogenic nuclide，TCN）10Be暴露测年手段，证实白玉冰期晚阶段冰碛垄的年代对应于MIS 2。目前，随着测年技术手段的发展及冰川数值年代的增加，波堆藏布流域的冰川演化序列日益完善，发生于MIS 12（甚至更早）［3］、末次冰消期［4-5］、新冰期以及小冰期［6］等时期的冰川作用也相继得到年代学的证实。然而，末次冰期早期的冰川作用年代还处于推测阶段，此次冰川作用究竟发生于MIS 4还是MIS 3，这一问题尚存在争议。因此，本文在野外地貌考察的基础上，运用光释光（optically stimulated luminescence，OSL）测年技术，对波堆藏布流域疑似形成于末次冰期早阶段的典型冰川/冰水沉积物进行定年，以期探讨末次冰期早阶段冰川作用时间及其古气候意义。

1 研究区概况

波堆藏布江发源于念青唐古拉山东段南坡，自西北向东南流100 km后汇入帕隆藏布江，流域面积约1.46×105km2（图1）。帕隆藏布江自东向西，至通麦汇合易贡藏布江，南流注入雅鲁藏布江。来自孟加拉湾的强大水汽，沿着雅鲁藏布江谷地深入帕隆藏布流域，带来丰沛的降水。根据波密气象站（29.87°N、95.77°E，海拔2 736 m）的观测数据，研究区年降水量约890 mm，其中6—9月的降水量占比超过了80%（1981—2011年，http：//data.cma.cn），而平衡线附近的降水量可达1 000～3 000 mm［7］。充沛降水加之高大地形的抬升作用为本区海洋性冰川发育提供了优越的条件。波堆藏布流域共有现代冰川482条，总面积974 km2［8］（图2）。流域内最高峰则普峰，海拔为6 364 m，其上发育了规模最大的现代冰川则普冰川，长达19.2 km，面积76.2 km2［9］。则普冰川同其他支沟冰川在第四纪期间多次进退于波堆藏布谷地，留下了丰富的古冰川作用遗迹。

图1 研究区地理位置及其DEMFig.1 Location of the study area and its DEM

测年研究表明，流域内最早的冰川作用为尼通冰期，相当于MIS 12甚至更早，典型冰碛物主要遗存于波堆藏布主谷3处已相当残破的高谷肩上［3］。古乡冰期冰碛垄在主谷中表现得蔚为壮观，从上游育仁乡至白玉村沿着谷坡断续分布，另在波堆藏布汇流处卡达桥附近也保存了一处相对完好的高大侧碛垄（图2）。根据冰碛垄的保存位置，当时支谷冰川都汇入了主谷，形成长约100 km、规模为现代冰川6～8倍的复合型山谷冰川［1］。末次冰期最盛期时，波堆藏布主谷和亚龙藏布谷地均被古冰川所占据，冰川末端分别到达白玉村和达龙村附近，冰川规模为现代的4～7倍［1］。末次冰期时，白玉冰川和珠西冰川均前进至沟口，在沟口处留下了规模宏大的终碛垄（图2）。形成于末次冰期晚阶段的内终碛垄，可分为2～3道，冰碛垄形态相对完好，拔河高度介于200～350 m之间。周尚哲等［2］、Ou等［5］对白玉沟内冰碛垄进行了测年，其中漂砾的TCN10Be暴露年代介于（11.1±1.9）～（18.5±2.2）ka之间，冰碛物及冰水沉积物的石英单颗粒OSL年代则处于（8.0±2.5）～（14.4±4.2）ka之间，认为内冰碛垄为MIS 2期间冰川作用的产物。推测形成于末次冰期早阶段的外终碛垄，其形态大多残破，垄脊相对平缓，高度也大为降低，拔河约100～200 m（图2）。

除主要支沟外，疑似形成于末次冰期早阶段的冰碛垄在亚龙藏布谷地中也有保存（图2～3），李吉均等［1］推测冰期最盛时的古冰川长达75 km直到沟口附近。野外考察发现，亚龙藏布谷口南岸和北岸均保存有高低两套侧碛垄，其中，北岸高侧碛垄的规模最大。波堆藏布江于谷口处将冰碛垄切割成一个高达150 m的典型三角面，突起的垄脊从谷口往上游延伸约4.5 km，与松龙沟高大终碛垄相接，最后消失于谷坡上。相应地，冰碛垄的拔河高度则由200 m递增至320 m。南岸高侧碛垄因受冲沟的侵蚀作用，断续分布于谷坡上，垄脊的形态不明显。与高侧碛垄一样，低侧碛垄在北岸连续分布，而在南岸间断保存。北岸低侧碛垄垄脊形态明显，上游与松龙沟终碛垄相接，向下游延伸约4 km，末端高出地面不到2 m，最后尖灭于河流阶地沉积物中。保存于南岸的低侧碛垄，受河流改造作用明显，垄脊相对平坦，顶部最宽处可达70 m。野外调查发现，高低侧碛垄在上游也有相当程度的保存。高侧碛垄向上游一直断续分布，直到海拔约3 700 m处，此处拔河高度约500 m，而低侧碛垄仅在松龙沟上游附近有所保存（图2）。

图2 波堆藏布江流域现代冰川与古乡冰期和末次冰期冰碛垄分布Fig.2 Modern glaciers and moraines constructed during Guxiang Glaciation and Last Glaciation in the Bodui Zangbo River catchment

亚龙藏布谷口附近广泛发育了两级冰水阶地（图3）。这两级阶地的阶地面均相对平坦，T2阶地主要分布于波堆藏布东岸和亚龙藏布北岸，面积达4 km2，拔河高度20～50 m。T2阶地在亚龙藏布南岸的分布面积较小，为0.8 km2，拔河高度30～80 m。T1阶地在亚龙藏布南岸和北岸均有分布，两处的分布面积均约为1 km2，拔河高度3～10 m。从地貌上看，阶地发育时间要晚于冰碛垄，低侧碛垄末端部分被T2阶地沉积物所覆盖。

图3 亚龙藏布江与波堆藏布江汇流处附近的冰川地貌、河流阶地及采样点位置（a）和野外照片（b，c，d）Fig.3 Glacial landforms，river terraces，sampling sites（a）and field photos of moraines（b，c）and terrace T2（d）near the confluence of the Yalong Zangbo and Bodui Zangbo Rivers

2 研究方法

2.1 样品采集

在亚龙藏布谷地低侧碛垄采集了4个（YL01-04）光释光测年样品（图3～4）。YL01的采样剖面，砾石大小混杂，无分选，但磨圆度较好，胶结程度高。样品采自一处长约2 m，宽20～30 cm的冰水透镜体中，冰水夹层为粉砂与细砾互层，以粉砂为主。YL02的采样冰碛层，砾石无分选，但有一定磨圆度，胶结程度高。样品采自砾石和粉砂互层的粉砂层中，透镜体周边出露有直径1 m的棱角状砾石。YL03的采样剖面以直径10 cm以内次棱角和次圆状的砾石及粉砂为主，胶结较好。新鲜剖面可见粉砂和细砾互层，采样物质为粉砂。YL04也采自附近粉砂与细砾互层处的粉砂层中。以上4个样品物质均为冰水夹层中的粉砂，来源于冰面湖或冰下河道沉积。

冰川作用时代可与相应的冰水阶地的发育年代相对应，因此在T2阶地剖面采了一个光释光年代样品（YLT18），用以验证冰碛垄的形成年代。YLT18的采样剖面出露高度约4 m，砾石直径大多介于5～30 cm之间，磨圆程度高，砾石之间充填砂及细砾，有分选。出露剖面的下部夹有一处长约2 m、宽约20 cm的粗砂层，样品采自这一粗砂层中。

采样时先挖开新鲜剖面，将一端塞了黑色塑料袋和棉花的不锈钢钢管（长20 cm，直径5 cm）砸入剖面中。取出钢管，用棉花和黑色塑料袋先后套住钢管的另一端，通过遮光的管道胶带将钢管的两端密封。

2.2 前处理

在红光实验室取出钢管两端少部分可能已经曝光的样品，用于含水量和年剂量的测定。

样品前处理之前先用300μm干筛将粗颗粒的砾石去除。将样品用10%稀盐酸和30%双氧水依次浸泡，去除样品中的碳酸盐和有机质。对于冰碛物样品，干筛提取38～63μm组分（中颗粒组分），用氟硅酸浸泡样品2周之后，通过10%稀盐酸去除样品与氢氟酸反应过程中产生的氟化物沉淀。对于冰水阶地样品，干筛提取90～125μm组分（粗颗粒组分），通过多钨酸钠重液分离出石英颗粒。用40%氢氟酸浸泡40 min，去除长石，随后用10%稀盐酸浸泡样品30 min，以去除氟化物。通过红外检测方法检测石英颗粒的纯度。若长石信号较高，则需再用氟硅酸或氢氟酸刻蚀样品，直到红外信号达到较低的水平（IRSL/OSL＜10%）为止。将纯石英样品用硅油均匀粘贴在直径0.97 cm的不锈钢测片中，样品的粘贴直径平均约为0.6 cm。

图4 采样点剖面图Fig.4 Photos showing the sampling sections

2.3 等效剂量与年剂量测试

等效剂量（equivalent dose，De）采用单片再生剂量法（single aliquot regeneration-dose，SAR）测试［10］。实验在华南师范大学地质年代学实验室完成，测试仪器为丹麦RisøTL/OSL DA-20全自动释光仪，内置90 Sr/90 Y人工β辐射源用于辐照。预热温度选定为260℃（10 s），试验剂量的预热温度为220℃（10 s）。测试石英释光信号的激发光源为波长（470±30）nm的蓝光二极管（强度90%），激发温度为120℃（40 s）。OSL信号通过7.5 mm的HoyaU-340滤光片进入9235QA光电倍增管内被接收并记录。

年剂量率计算中的U、Th和K的含量在北京第二核工业研究所测定。样品的实测含水量均在3%以内，考虑到样品运输及储存过程中的水分损失，结合当地降水量，石英中颗粒年代样品采用（10±5）%的估计值，粗颗粒采用（5±4）%的估计值［11-13］。宇宙射线的年剂量贡献根据样品的地理位置、海拔高度和采样深度计算［14］。剂量率根据Aitken的公式和参数进行计算［15］。

3 实验结果

3.1 释光特征

图5展示了样品YL02的光释光衰减曲线和生长曲线。从衰退曲线可以看出，OSL信号在蓝光激发下的最初几秒衰减很快，表明OSL信号以快速组分为主。该样品的循环比在0.95～1.05之间，说明试验剂量对感量变化的校正较为理想。样品的热转移值介于0.15%～0.76%之间，表明测试过程中热转移效应对等效剂量的贡献较低，对年代结果产生的影响很小。图6为所有样品的De散点图。

图5 样品YL02（虚线）和标样（实线）天然剂量的衰减曲线（小图为样品YL02的生长曲线）Fig.5 OSL decay curves of natural dose of Sample YL02（dashed line）and standard sample（solid line）（Theinserted figure shows the growth curve of Sample YL02）

图6 样品D e散点图Fig.6 D e scatter plots for the samples

3.2 预热坪与剂量恢复实验

为选择合适的SAR程序的预热温度，本文选择样品YL02进行预热坪和剂量恢复实验。对样品制备20个样品测片用于预热坪实验，测片分为5组，分别在220℃、240℃、260℃、280℃、300℃的条件下预热10 s，试验剂量的预热温度则为220℃，预热时间为10 s。每组预热温度的De平均值如图7（a）所示，结果相差不大，且各个温度条件下的循环比和热转移值都在合理范围内［图7（c）］。参考前人相关研究［16-17］，本文选择260℃作为预热温度，实验剂量的预热温度为220℃，预热时间均为10 s。

对于剂量恢复实验，先准备10个样品测片，将其释光信号清除后，对每个测片给定200 Gy的辐照剂量。将测片分为5组，每组分别在220℃、240℃、260℃、280℃、300℃的条件下预热，测试程序如上。每组测片的De平均值介于198.78～203.59 Gy之间，与给定剂量的比值均在0.9～1.1之间［图7（b）］，且各个温度条件下的样品循环比和热转移值均介于合理阈值范围内［图7（d）］，说明以上测试程序可以恢复已知的剂量。

图7 样品YL02预热坪和剂量恢复实验结果及其循环比和热转移Fig.7 Results of preheat plateau test（a），recycling ratio，and recuperation（c）of Sample YL02，and results of dose recovery test（b），recycling ratio，and recuperation（d）of Sample YL02

4 年代结果

OSL年代结果见表1。De值采用中值年代模型（central age model，CAM）计算得出。

5 讨论

5.1 测年结果分析

如表1所示，冰碛物YL01-04这4个样品的CAM年代分别为（56.4±4.2）ka、（60.3±3.5）ka、（65.9±3.9）ka和（31.2±2.1）ka，河流阶地沉积物样品YLT18的CAM年代为（65.2±6.0）ka。根据野外调查及地貌地层关系，T2阶地的形成时代应该接近或稍晚于冰碛垄的建造时代，基于此，YL04的年代应当剔除。YL04年代结果偏年轻可能是由于后期曝光所致。其余3个冰碛物样品（YL01、YL02和YL03）的年代介于（56.4±4.2）～（65.9±3.9）ka之间，这一年代结果可与保存于白玉村南侧谷坡附近的侧碛垄上的8个漂砾TCN10Be暴露年代结果相对应［（57.7±5.7）～（72.0±7.4）ka（未发表）］。这两处冰川沉积物的年代均对应于MIS 4，表明冰碛垄为MIS 4冰期冰川作用的产物。松龙沟高大冰碛垄上的漂砾TCN10Be暴露年代位于（14.8±0.9）～（19.5±1.2）ka之间（未发表），表明其形成于MIS 2，根据地层地貌接触关系，末次冰期早阶段（MIS 4）的冰川规模要大于晚阶段（MIS 2），与波堆藏布主谷的情况一致。此外，基于以上测年结果及地层地貌关系，推测高侧碛垄可能形成于古乡冰期。

表1 光释光测年结果Table 1 OSL dating results

5.2 青藏高原及周边山地MIS 4冰川遗迹

冰川沉积物定量测年研究表明，MIS 4冰期冰川作用在青藏高原及周边山地广泛存在（图8）。在青藏高原东北部，东昆仑东吉康纳湖地区［18］和巴颜喀拉山［19］等地的冰碛物10Be暴露年代结果显示，MIS 4期间两地的古冰川曾大范围扩张。MIS 4冰期冰川作用遗迹在青藏高原中部的唐古拉山［20-21］以及南部的念青唐古拉山西段［21］和申扎杰岗日［22］等地区也有发现，这些地区冰碛物的10Be暴露年代结果揭示当时的古冰川规模要大于后续冰期冰川作用规模。在喜马拉雅山地区，如中段南坡Gorkha地区［23］，西段北坡纳木那尼峰地区［24］，南坡加瓦尔西北部［25］，以及西北段印度西北部地区［26-29］，这些地区的漂砾10Be暴露年代或冰碛物OSL年代皆证实了MIS 4冰期冰川作用的存在。形成于MIS 4期间的冰川遗迹在青藏高原西部，如阿里地区查西康米仁农坝河谷［30］，西昆仑［31］等地区也有保存。相较于以上地区，主要受西风影响的青藏高原西北部及天山等地区MIS 4冰期冰川遗迹数量较多。冰川沉积物10Be暴露年代结果表明，喀喇昆仑山南部熊彩岗日地区［32］，东南段拉达克山脉［33］、洪扎谷地［34］、塔什库尔干谷地［35］均保存有此次冰川作用的遗迹。在帕米尔高原地区，如东帕米尔瓦恰谷地［36］，帕米尔（塔吉克斯坦）雅什库尔湖［37-39］等地区也发现有MIS 4冰期冰川作用遗迹。天山地区，如东天山的博格达峰［40］、天格尔峰［41］和喀尔里克山［42］等地区，中天山的木扎尔特河谷［43］，以及西天山的阿莱山脉科克苏谷地［44］、特斯基阿莱太山脉［45-46］、阿拉巴什盆地［47-48］、吉萨尔山脉（塔吉克斯坦）［49］、伊尼尔切克谷地［50］等地区，这些地区冰碛物的OSL年代、ESR年代或10Be暴露年代结果也证实古冰川于MIS 4期间曾大规模扩张。此外，横断山脉东侧中段螺髻山［51］及横断山脉腹地白马雪山［52］等地区也保存有4阶段冰川作用遗迹。

图8 青藏高原及其周边山地有测年数据的MIS 4冰川遗迹分布Fig.8 Dated glacial landforms formed during MIS 4 on the Tibetan Plateau and itssurrounding mountains：1.Luoji Mountain［51］，2.Baimaxue Shan［52］，3.eastern Nyainqêntanglha Range（thisstudy），4.Macha Khola Valley，Gorkha Himal［23］，5.Lake Donggi Cona，eastern Kunlun Shan［18］，6.Bayan Har Shan［19］，7.Tanggula Shan［20-21］，8.western Nyainqêntanglha Range［21］，9.Mount of Jaggang，the Xainza Range［22］，10.Gurla Mandhata（Naimon’anyi）［24］，11.NW Garhwal，Central Himalayas［25］，12.NW Himalaya of north-west India［26-27，29，41］，13.Miren Nongba Valley，Chaxikang［30］，14.Mawang Gangri Range［32］，15.Ladakh Range［33］，16.Hunza Valley，Karakoram Mountains［34］，17.Tashkurgan Valley，southeast Pamir［35］，18.West Kunlun，Karakax［31］，19.Waqia Valley，Chinese Pamir［36］，20.Lake Yashilkul，Pamir，Tajikistan［37-39］，21.Karlik Range，easternmost Tian Shan［42］，22.Bogeda Peak area，Tianshan Range［40］，23.Ala Valley of the Tianger Range，eastern Tian Shan［41］，24.Muzart River Valley，central Tianshan Range，China［43］，25.Inylchek Valley，eastern Kyrgyz Tian Shan［50］，26.Teskey Ala-Too Range，Kyrgyz［45-46］，27.Ala Bash Basin［47］，28.Koksu Valley of the Alay Range［47］，29.Gissar Range，Tajikistan［49］

综上，无论是西风区还是季风区，MIS 4期间古冰川均曾大规模扩张，表明4阶段冰川作用在青藏高原及其周边山地的出现具有普遍性。前人研究认为，青藏高原及其边缘山地MIS 4期间的冰川活动主要受控于中纬度西风［33，42，53-54］，冰期时西风带来的低温和降水增加了冰川物质积累量，进而驱动冰川前进。特别是在中纬度西风为主要水汽来源的地区，如天山［41，47，49］、帕米尔地区［35，44］，MIS 4冰期为末次冰期最盛期，冰川作用范围为末次冰期以来最大。

5.3 藏东南MIS 4冰期冰川作用发生的气候条件

从MIS 5a开始，北半球低纬度太阳辐射持续减弱，至71 ka左右达到最低值［55］［图9（b）］。青藏高原东部若尔盖盆地RM孔140 ka以来的湖泊碳、氧同位素记录显示，MIS 4期间气候较冷［56］，表明随着太阳辐射减弱，青藏高原季风区的气温也逐渐降至低值。根据云南腾冲青海湖泊沉积岩中所揭示的低纬陆地高分辨率年均温记录，MIS 4期间气温要比MIS 5a低2～3℃，而这样的温度变幅基本达到了季风区冰期-间冰期的变幅［57］［图9（e）］。因此，藏东南地区山地冰川于MIS 4期间大范围扩张正是对冰期气温下降的响应，与青藏高原及其边缘山地的部分研究结果相符［42，49，53-54］。深海钻孔中的δ18O记录也表明，MIS 4时期是大陆冰量最大时期之一［58］［图9（c）］。

这一时期的石笋δ18O记录指示了偏弱的季风强度［59］［图9（d）］，表明藏东南地区在该时期由西南夏季风带来的降水减少，此阶段增强的中纬度西风可能成其主要水汽来源［33，53-54］。到了MIS 3时期，太阳辐射量的增加及季风的增强导致冰川开始后退［33，53］。随着冰后期降水量及冰川融水量的增多，河流流量增加及搬运能力增强，大量的冰水沉积物首先堆积在河谷中，后期河流下切形成T2阶地。

图9 西藏东南部MIS 4冰期冰川作用与亚洲季风区的古气候记录对比Fig.9 Climatic archives relative to the glacial activities during MIS 4 across the monsoonal Asia：probability density estimate of the OSL ages（Black dots show age of moraine，and the blue one indicates that of river terrace）（a），summer（June-July-August）mean isolation at 30°N［55］（b），benthic δ18O records［58］（c），speleothemδ18O data from Dongge，Hulu，and Sanbao Cavesin China［59］（d），reconstructed annual air temperature at Lake Tengchongqinghai［57］（e）

6 结论

（1）基于前人研究和野外地貌调查，运用光释光测年技术，对西藏东南部波堆藏布流域疑似形成于末次冰期早阶段的冰碛垄进行测年，年代介于（56.4±4.2）～（65.9±3.9）ka之间，对应于MIS 4冰期。这一年代结果与青藏高原其他地区及边缘山地同时期的冰川测年结果具有可比性，表明MIS 4冰期冰川作用发生的普遍性。

（2）通过对比北半球低纬度地区夏季太阳辐射和亚洲季风区古气温和古降水指标记录，藏东南地区MIS 4期间的冰川扩张可能是对北半球低纬度地区太阳辐射减弱及气温下降的响应，与季风降水无关。MIS 4之后，印度夏季风增强，气温回升及降水增加导致冰川后退。

“白玉冰期”是李吉均先生于二十世纪七八十年代主导提出的第四纪冰川重要冰期概念，谨以此文，纪念李吉均院士！