藏东南波堆藏布江流域古乡冰期冰川重建
1.
2.
Glacier reconstruction of the Guxiang Glaciation in the Bodui Zangbo Valley in southeastern Tibet
1.
2.
通讯作者:
编委: 武俊杰
收稿日期: 2021-06-01 修回日期: 2021-06-11 网络出版日期: 2021-07-29
基金资助: |
|
Received: 2021-06-01 Revised: 2021-06-11 Online: 2021-07-29
关键词:
Keywords:
本文引用格式
覃长雄, 许刘兵, 谢金明, 周尚哲.
QIN Changxiong, XU Liubing, XIE Jinming, ZHOU Shangzhe.
0 引言
冰川是寒冷气候的产物,对气温和/或降水变化十分敏感[1],古冰川相对现代冰川的物质平衡线高度(equilibrium-line altitude,ELA)、面积、长度、以及体积等特征参数的变化,对了解古气候和预测未来气候变化至关重要。其中,现代冰川的各项特征参数可根据实地观测、测量、遥感影像解译等方法直接估算得到,对古冰川而言,通常是根据现代冰川与当前气候环境要素(如气温、降水、地形等)的统计关系,再结合冰川地貌和年代学去反演过去不同时期冰川的各项关键参数。从目前的研究现状来看,对古冰川发育环境及特征参数的重建研究主要集中在末次冰期最盛期(Last Glacial Maximum,LGM)以来,主要原因是LGM以来的冰川地貌保存相对完好,地质年代也能准确的测定。
“古乡冰期”是20世纪80年代由李吉均等[2]提出来的冰期概念,同时提出的还有“白玉冰期”,这两个冰期名称在青藏高原第四纪冰川作用众多地方性冰期名称中非常具有代表性,数十年来一直被广泛沿用。其中,古乡冰期指代的是传统意义上的“倒数第二次冰期”,已有的宇宙成因核素10Be地表暴露测年结果[3]显示其发生于海洋氧同位素阶段(marine oxygen isotope stage,MIS)6。虽然古乡冰期概念已被众人接受,且发生年代也已界定,但古乡冰期命名地即藏东南波堆藏布江流域的冰川作用规模(如范围、冰厚及冰量等)和ELA等关键参数仍有待进一步考证。本文在多次详尽野外调查的基础上,结合高清遥感影像/卫星照片及定量测年结果,对波堆藏布江流域古乡冰期的冰川作用范围进行了重建,并进一步尝试运用冰面纵剖面(glacier surface profile,GSP)模型、积累区面积比率(accumulation area ratio,AAR)法和末端-源头高度比率(toe-to-headwall altitude ratio,THAR)法等方法对波堆藏布江流域“古乡冰期”的古冰川特征参数进行重建。
1 古乡冰期冰川遗迹及其年代
图1
古乡冰期指的是传统意义上的“倒数第二次冰期”,由李吉均等[2]于20世纪80年代提出,其命名地位于藏东南波堆藏布江与帕隆藏布江交汇处下游约17 km的古乡镇,即古乡冰期时,流域内支谷冰川汇入主谷(即波堆藏布江谷地),形成复合型山谷冰川,冰舌到达了古乡镇附近,遗存有较大规模的终碛-侧碛垄。
除古乡镇一带的冰碛垄外,波堆藏布江主谷两侧亦保存有古乡冰期的沉积物(图1)。保存在两侧谷坡、形态完好且规模宏大的侧碛垄从上游育仁一直延伸至白玉村附近,拔河高度从550 m降至260 m左右[图2(a)和(b)];自白玉向下游保存不甚理想,只见部分谷肩上残存有零星的冰碛物;但在波堆藏布江汇入帕隆藏布江处的卡达桥附近,尚完好地保留着拔河约200 m、长约3.2 km的高大侧碛垄,其上分布着大量的花岗岩漂砾[图2(c)]。周尚哲等[3]对卡达桥附近古乡冰期冰碛垄进行了10Be暴露测年,测年结果(约130 ka)显示古乡冰期发生于MIS 6。除波堆藏布江主谷外,古乡冰期冰碛垄在支谷亚龙藏布江两侧谷坡上也有保存,其中,亚龙藏布汇入波堆藏布处的侧碛垄保存得最为完整[图2(d)]。亚龙藏布江北岸的侧碛垄垄脊非常典型,从谷口往上游延伸约4.5 km,拔河高度由200 m升至320 m,而谷口南岸的侧碛垄因受到支沟的侵蚀,垄脊较为平坦。冰碛垄在上游也有断续保存,直到海拔3 700 m处,此处拔河约500 m。可见,古乡冰期时,波堆藏布江流域内支谷冰川都汇入了主谷,形成了主干长达100 km左右的复合型山谷冰川,冰川规模可能为现代的数倍[2]。
图2
图2
育仁(a)、林琼(b)、卡达桥(c)和亚龙藏布谷地(d)的古乡冰期冰碛垄照片
Fig.2
Photos of moraines deposited in the Guxiang Glaciation (Moraines in a, b, c, and d are located near Yuren, Linqiong, Kada Bridge, and Yalong Zangbo Valley, respectively)
2 研究方法
2.1 古乡冰期时冰川范围的划定
在多次详尽野外调查的基础上,结合已有冰川地貌年代学、高清遥感/卫星影像以及数字高程模型(digital elevation model,DEM),对古乡冰期冰川作用范围进行了重建(图3)。在划定古乡冰期冰川范围时,存在冰碛地貌的地方比较容易确定冰川作用范围,例如沿着波堆藏布江主谷从则普冰川至古乡镇这一段,以及亚龙藏布江与波堆藏布江交汇口附近,可以按照谷地两侧古乡冰期侧碛垄确定古冰川的作用范围(图1)。在古乡冰期冰碛地貌保留不理想甚至完全被侵蚀破坏的地方,如流域内的西南部、西北部以及中部,结合现代冰川的分布范围来大致划定古乡冰期冰川作用范围。而在流域内的东北部,既无明确的古乡冰期冰碛地貌,也没有现代冰川分布,主要依据形态清晰的冰斗和槽谷来划定古冰川作用范围。波堆藏布江流域内仍然存在还未定年或无法定年的冰碛垄,要准确地确定古乡冰期时的冰川作用范围,仍需进一步实地考察冰川沉积地貌以及进行相关的定年研究。因此,本文采取上述策略来界定古乡冰期冰川作用范围,可能与实际情况存在出入,但依据划定的古冰川范围估算出的古乡冰期ELA结果来看,与Zhou等[5]的研究结果基本一致(见3.2节)。
图3
图3
古乡冰期时波堆藏布江流域古冰川范围重建(彩色实线代表35条中流线分布,红色虚线代表流域边界)
Fig.3
Reconstruction of glacier outline in Guxiang Glaciation in the Bodui Zangbo River catchment (The solid colored lines represent glacier centerlines, and the dotted red line represents the boundary of the Bodui Zangbo River catchment)
2.2 古冰川模拟原理
将根据上文方法得到的古乡冰期冰川作用范围导入ArcGIS中,与分辨率为30 m的DEM数据(ALOS World 3D-30 m,
此模型基于的假设为冰川是理想的塑体,其形变由冰的重量及冰表面梯度引起的驱动应力(driving stress,
式中:ρ为冰川密度(900 kg∙m-3);g为重力加速度(9.81 m∙s-2);H为冰川厚度;h为冰表面高程;x为水平坐标(以冰川末端为原点,冰底和冰面高程为纵轴,距冰川末端的长度为横轴的坐标系)。有时用
但该式不适用于冰川坡度较大的地区。对于不规则的冰下地形,Benn等[7]指出必须沿着中流线按连续的点(编号1,2,…,i,i+1,i+2)计算,冰的表面坡度表示为
式中:Δx为每个点的间隔距离,本文取Δx=200 m。
因此,由
式中:H=h-B,B为冰川底部高程,此式可以解出冰面高程h。
为了解决冰舌处冰厚为0的问题,van der Veen[14]进一步推导得到
但该模型是以冰盖出发进行推导的,冰盖不受山谷两侧阻力的约束,其流动需要的驱动应力等于冰川底部剪切应力。而山谷冰川则需要考虑山谷两侧阻力的影响,因此引入形状因子f。
式中:
在冰川中流线处,
2.3 古冰川模拟过程
图4
图5
图5
根据冰川中流线E4模拟的冰川表面高程和其他参数
Fig.5
Characteristic parameters of the glacier reconstructed from the glacier centerline E4 (Figure a, b, c, and d represents basal shear stress, ice thickness, surface elevation, and shape factor, respectively)
使用上述模拟过程,对其他34条中流线进行模拟,将模拟结果分别用反距离加权(inverse distance weight,IDW)、泛克里金(universal Kriging,UK)和普通克里金(ordinary Kriging,OK)插值法进行插值,最终得到古冰面高程,而古冰川厚度可由模拟古冰面高程减去古冰川冰底高程,结果见图6。
图6
图6
古冰面高程(a,b,c)和古冰川厚度(d)模拟结果以及不同古冰面高程下的ELA重建结果(a,b,c中的黑线)(a,b,c分别为根据IDW、UK和OK插值法得到的古冰面高程)
Fig.6
The modeled glacier surface elevation (a,b,c) and glacier thickness (d) and reconstructed ELAs (black lines in a,b,c) based on different glacier surface elevation (Figures a,b,c show the glacier surface elevation using IDW, UK, and OK interpolation methods, respectively)
2.4 冰川ELA估算方法
由Meier[17]提出的AAR法被广泛用于计算古冰川和现代冰川的ELA,这一方法的关键是选择合适的AAR值。研究表明,分布于不同纬度的冰川,有不同的AAR值。中高纬度冰川的AAR值介于0.5~0.8[18]或0.55~0.65之间[19],而低纬地区则可能高达0.8[20]。AAR法假定冰川物质平衡随高度的变化趋近于一条直线,并且考虑了冰川的面积-高程信息。Kern等[21]认为冰川面积S与AAR值之间正相关,且呈对数变化关系:AAR=0.0648×lnS+0.483。对于一些数据缺乏的地区,当S为0.1~1 km²时,AAR=0.44±0.07;当S为1~4 km²时,AAR=0.54±0.07;当S大于4 km²时,AAR=0.64±0.04。Kaser等[22]认为热带地区无积雪覆盖的冰川最佳AAR范围在0.65~0.70之间,且这个范围对于热带有积雪覆盖的冰川同样适用,地貌证据也表明了这个取值范围同样适用于古冰川。
表1 第一次冰川编目[23]下波堆藏布江流域内5条现代冰川的AAR值及估算的第二次冰川编目[4]下各ELA值
Table 1
冰川名称 | 积累区面积/km² | 总面积/km² | 第一次冰川编目AAR | 第二次冰川编目ELA/m |
---|---|---|---|---|
平均值 | 0.73 | 4 451 | ||
朱西冰川 | 7.08 | 12.9 | 0.55 | 4 011 |
白玉冰川 | 16.75 | 24.4 | 0.69 | 4 048 |
则普冰川 | 59.44 | 72.4 | 0.80 | 4 390 |
日母冰川 | 12.52 | 15.2 | 0.80 | 4 849 |
关星冰川 | 37.62 | 47.6 | 0.79 | 4 956 |
THAR法也常用于估算古冰川和现代冰川的ELA,THAR值为ELA与冰川最低点高程之差和冰川最高点与最低点高程之差的比值。但该方法并未涉及到冰川的物质平衡或者冰川的面积-高程信息,THAR值的选取有一定的主观性[26-27]。Charlesworth[28]、Manley[29]曾使用冰川的中值高度来代替ELA。THAR值小于0.5时的ELA小于中值高度[30-33]。Meierding[32]在研究科罗拉多冰川时,认为在0.35~0.40之间的THAR值比较合适。使用该方法估算古冰川ELA时,古冰川的末端高程主要根据冰川终碛的高度计算,即采用沿终碛顶端的最低高度作为冰舌的最低高度[24];但最高高程不容易确定。为了选取合适的THAR值估算现代冰川和古冰川ELA,本文以第二次冰川编目[4]时观测的5条现代冰川(分别为朱西冰川、白玉冰川、则普冰川、日母冰川和关星冰川)为例(图1),将每条冰川上缘处的最高点作为最大高程,使用不同THAR值估算其ELA,估算结果见表2。结果显示,相同的THAR值下这5条冰川的ELA有着类似的变化趋势,即随着冰川海拔的增高其ELA增大,本文取这5条冰川ELA的平均值来代表整个流域的ELA。继而与第二次冰川编目下这5条冰川ELA(表1)的平均值4 451 m比较,发现THAR=0.35时二者比较接近(表2),因此,在重建古乡冰期ELA时,THAR值取0.35。
表2 波堆藏布江流域不同THAR值下的5条现代冰川ELA估算结果
Table 2
冰川名称 | ELA/m | ||||
---|---|---|---|---|---|
THAR=0.30 | THAR=0.34 | THAR=0.35 | THAR=0.36 | THAR=0.40 | |
平均值 | 4 345 | 4 435 | 4 458 | 4 480 | 4 570 |
朱西冰川 | 3 923 | 4 027 | 4 054 | 4 080 | 4 184 |
白玉冰川 | 3 996 | 4 087 | 4 110 | 4 133 | 4 225 |
则普冰川 | 4 246 | 4 362 | 4 392 | 4 421 | 4 538 |
日母冰川 | 4 744 | 4 814 | 4 831 | 4 848 | 4 917 |
关星冰川 | 4 817 | 4 884 | 4 901 | 4 918 | 4 985 |
3 结果与讨论
3.1 冰川规模变化
从中流线E4上的剪切力分布[图5(a)]发现,下游(距冰川末端的长度为0~45 km)剪切力变化较小,中游(距离冰川末端的长度为45~75 km)剪切力变化较大,上游(距离冰川末端的长度为75~120 km)剪切力较小,总体分布在7~120 kPa之间。经过计算,剪切力
3.2 古冰川ELA
冰面纵剖面模型只能模拟中流线上的冰面高程和冰川厚度,因此重建古冰川整个冰面高程时需要进行插值。本文使用了IDW、UK和OK等三种常用空间插值方法,得到了三个不同的古冰面高程分布(图6),结合AAR工具[12]和THAR法计算了相应的冰川ELA,结果显示:当采用AAR法时,古乡冰期时冰川的ELA分别为4 021 m、3 867 m、3 847 m。采用THAR法时,估算出古乡冰期时冰川ELA分别为3 720 m、3 649 m、3 649 m(表3)。从估算的结果来看,AAR法与THAR法得到的古冰川ELA相差较大,可能有以下几个方面的原因:①使用AAR法计算ELA时需要考虑冰川形状、类型和规模等特征[26],而THAR法没有考虑冰川形状或不同谷地的面积-高程变化[34];②采用AAR法时,可能会因为重建的古冰川表面面积不准确而产生误差;③使用THAR法时,若雪崩对冰川的积累影响较大,则会难以界定冰川的最高高程。因此,本文取两种估算方法得到的古冰川ELA平均值作为最终结果,这样能有效避免使用单一计算方法所带来的误差。
表3 波堆藏布江流域冰川ELA
Table 3
冰川作用时间 | ELA/m | 下降值/m | ||
---|---|---|---|---|
AAR=0.73 | THAR=0.35 | 平均值 | ||
现代 | 4 451* | 4 458 | 4 455 | |
古乡冰期(基于IDW) 古乡冰期(基于UK) | 4 021 3 867 | 3 720 3 649 | 3 871 3 758 | 584 697 |
古乡冰期(基于OK) | 3 847 | 3 649 | 3 748 | 707 |
不同的空间插值方法得到的古冰面高程对于重建古冰川冰量、ELA等重要参数影响较大。为选择合适空间插值方法,本文在中流线上的数千个插值点中随机抽取10个点,对比其不同插值结果与模拟值(表4中的实测高程)之间的误差,结果见表4。从误差指标可以看出,IDW插值的古冰面高程值的绝对误差小于5 m,与实测值之间的相对误差在0.01%~0.1%之间。相对地,UK和OK插值的绝对误差分别可达347 m和277 m,相对误差分别在0.002%~12%和0.002%~10%之间。由此可见,IDW插值法的古冰面高程插值精度明显好于UK和OK插值法。主要原因是IDW插值法对已知高程点设置权重值,并搜索插值点与已知点之间的最短距离来进行空间插值,而UK和OK插值法要求已知高程点的密度分布较大,对于密度较小的区域,空间插值精度较低。本文选取中流线上每个点之间的间隔为200 m,密度较小,已知点较多,因此UK和OK插值法的精度要低于IDW插值法。
表4 不同空间插值方法得到的古冰面高程精度对比
Table 4
编号 | 实测高程/m | 高程插值/m | 绝对误差/m | 相对误差/% | ||||||
---|---|---|---|---|---|---|---|---|---|---|
IDW | UK | OK | IDW | UK | OK | IDW | UK | OK | ||
1 | 2 678.742 | 2 678.072 | 3 025.929 | 2 956.566 | 0.67 | 347.19 | 277.82 | 0.025 | -12.961 | -10.371 |
2 | 2 705.931 | 2 706.710 | 2 923.167 | 2 882.967 | 0.78 | 217.24 | 177.04 | -0.029 | -8.028 | -6.543 |
3 | 2 718.123 | 2 719.080 | 2 834.148 | 2 811.306 | 0.96 | 116.03 | 93.18 | -0.035 | -4.269 | -3.428 |
4 | 2 734.359 | 2 733.184 | 2 911.702 | 2 881.433 | 1.18 | 177.34 | 147.07 | 0.043 | -6.486 | -5.379 |
5 | 2 753.043 | 2 753.531 | 2 753.087 | 2 753.087 | 0.49 | 0.04 | 0.05 | -0.018 | -0.002 | -0.002 |
6 | 2 776.828 | 2 776.500 | 2 980.947 | 2 946.629 | 0.33 | 204.12 | 169.80 | 0.012 | -7.351 | -6.115 |
7 | 2 783.399 | 2 783.711 | 2 971.017 | 2 940.639 | 0.31 | 187.62 | 157.24 | -0.011 | -6.741 | -5.649 |
8 | 2 791.636 | 2 790.393 | 2 871.972 | 2 859.999 | 1.24 | 80.34 | 68.36 | 0.044 | -2.878 | -2.449 |
9 | 2 813.275 | 2 813.493 | 2 920.051 | 2 909.309 | 0.22 | 106.78 | 96.03 | -0.007 | -3.796 | -3.413 |
10 | 3 889.426 | 3 885.117 | 3 991.450 | 3 989.136 | 4.31 | 102.02 | 99.71 | 0.111 | -2.623 | -2.564 |
4 结论
本文在冰川地貌外野实地考察和已有年代学为基础上,辅以高清遥感影像/卫星照片,对波堆藏布江流域古乡冰期时的冰川作用范围进行了重建;后进一步运用冰面纵剖面模型对波堆藏布江流域古乡冰期时古冰川进行了模拟;此外运用AAR和THAR方法对古乡冰期和现代冰川ELA进行了估算。古乡冰期时,波堆藏布流域冰川的平均厚度约为360 m,冰川总面积约为2 648 km²,体积约为953 km³,流域内的冰川覆盖率由古乡冰期的63%缩减至现在的22.4%;现代冰川ELA约4 455 m,古乡冰期时ELA约为3 871 m,较现代冰川下降了约584 m。
“古乡冰期”是李吉均先生主导提出的第四纪冰川重要冰期概念,谨以此文,纪念李吉均院士!
参考文献
Cosmogenic 10Be dating of Guxiang and Baiyu Glaciation
[J].
古乡冰期和白玉冰期的宇宙成因核素10Be定年
[J].
The second glacier inventory dataset of China (version
中国第二次冰川编目数据集(V
Glacial advances in southeastern Tibet during late Quaternary and their implications for climatic changes
[J].
Numerical reconstructions of valley glaciers and small ice caps
[M]// Denton G H, Hughes T J.
An ExcelTM spreadsheet program for reconstructing the surface profile of former mountain glaciers and ice caps
[J].
Reconstruction of Last Glacier Maximum paleo-glacier in the Qugaqie Valley in Nyainqentanglha Range
[D].
念青唐古拉山曲嘎切末次冰盛期古冰川重建
[D].
Cosmogenic 10Be surface exposure dating and glacier reconstruction for the Last Glacial Maximum in the Quemuqu Valley, western Nyainqentanglha Mountains, south Tibet
[J].
Kaluxung, south Tibet
[J].
藏南卡鲁雄峰枪勇冰川新冰期冰川发育探讨
[J].
Reconstruction past glacier since MIS
冈底斯中段冷布岗日嘎尔琼沟MIS3以来冰川重建及其平衡线高程估算
[J].
A GIS tool for automatic calculation of glacier equilibrium-line altitudes
[J].
The mechanics of glacier flow
[M].
Fundamentals of glacier dynamics
[M].
Automated modelling of spatially-distributed glacier ice thickness and volume
[J].
Proposed definitions for glacier mass budget terms
[J].
Recent variations in mass net budgets of glaciers in western North America
[C]//
Equilibrium line altitude of late Quaternary glaciers in the Southern Alps, New Zealand
[J].
Size specific steady-state accumulation-area ratio: an improvement for equilibrium-line estimation of small paleo-glaciers
[J].
Glacial and climate change in Boduizangbu Valley in Baiyu ice stage
[D].
波堆藏布谷地白玉冰期冰川与气候变化
[D].
Study on Quaternary glacial evolution in the southeastern Tibet and cosmogenic nuclides dating
[D].
西藏东南部第四纪冰川演化序列与宇宙成因核素测年研究
[D].
The methods for estimating the equilibrium line altitudes of a glacier
[J].
冰川物质平衡线的估算方法
[J].
A review on the reconstruction of paleoclimate by glacial landform in the Tibetan Plateau and adjacent mountains
[J].
青藏高原及毗邻山地利用冰川地貌重建古气候的研究综述
[J].
The late-glacial climate of north-west England
[J].
Modern and Wisconsinan snowlines in Alaska
[C]//
l
.: s.n.
Late Pleistocene glacial equilibrium-line in the Colorado Front Range: a comparison of methods
[J].
Himalayan glacial sedimentary environments: a framework for reconstructing and dating the former extent of glaciers in high mountains
[J].
Comparing reconstructed Pleistocene equilibrium-line altitudes in the tropical Andes of central Peru
[J].
/
〈 |
|
〉 |
