冰川冻土, 2021, 43(3): 724-734 doi: 10.7522/j.issn.1000-0240.2021.0014

藏东南波堆藏布江流域古乡冰期冰川重建

覃长雄,1, 许刘兵,1, 谢金明2, 周尚哲1

1.华南师范大学 地理科学学院，广东 广州 510631

2.兰州大学 资源环境学院，甘肃 兰州 730000

Glacier reconstruction of the Guxiang Glaciation in the Bodui Zangbo Valley in southeastern Tibet

QIN Changxiong,1, XU Liubing,1, XIE Jinming2, ZHOU Shangzhe1

1.School of Geography，South China Normal University，Guangzhou 510631，China

2.College of Earth and Environmental Sciences，Lanzhou University，Lanzhou 730000，China

通讯作者: 许刘兵，教授，主要从事地貌与环境演变研究. E-mail： xuliubing234@163.com

编委: 武俊杰

收稿日期: 2021-06-01   修回日期: 2021-06-11   网络出版日期: 2021-07-29

基金资助: 国家自然科学基金项目.  41771065.  41271077

Received: 2021-06-01   Revised: 2021-06-11   Online: 2021-07-29

作者简介 About authors

覃长雄，硕士研究生，主要从事第四纪冰川研究.E-mail：1391789666@qq.com , E-mail：1391789666@qq.com

摘要

“古乡冰期”是青藏高原第四纪冰川作用中最具代表性的冰期之一，指代的是倒数第二次冰期，其命名依据来源于藏东南波堆藏布江谷地、保存于古乡一带的终碛-侧碛垄，已有的宇宙成因核素10Be暴露测年结果显示其发生于海洋氧同位素阶段（MIS）6。然而，古乡冰期时波堆藏布江流域冰川作用范围、冰量及平衡线高度（ELA）等关键信息，仍有待进一步研究。结合冰川地貌野外考察、高清遥感影像/卫星照片及绝对数值年龄，对波堆藏布江流域古乡冰期时的冰川作用范围进行了重建，在此基础上结合数字高程模型，运用冰面纵剖面模型对波堆藏布江流域古乡冰期时古冰川进行了数值模拟。结果显示：古乡冰期时流域内冰川的平均厚度约360 m，冰川总面积为2 648 km²，体积约为953 km³，流域内的冰川覆盖率由古乡冰期时的63%缩减至现在的22.4%。此外，运用积累区面积比率（AAR）法和末端-源头高度比率（THAR）法对波堆藏布江流域现代和古乡冰期时的冰川平衡线高度进行了估算。其中，现代冰川ELA为4 455 m，基于反距离加权法插值的古冰面高程计算出古乡冰期时ELA为3 871 m，比现代ELA低了584 m，古冰川ELA估算结果与前人结果基本一致。

关键词： 古乡冰期 ; 波堆藏布江 ; 冰川重建 ; 平衡线高度 ; 积累区面积比率 ; 末端-源头高度比率

Abstract

The Guxiang Glaciation is one of the most typical glaciations on the Tibetan Plateau during the Quaternary and has been widely used to represent the penultimate glaciation on the plateau. It was named from the frontal-lateral moraine at Guxiang， over 100 km down valley from the source area around the Guanxing Glacier， in the Bodui Zangbo River catchment in southeastern Tibet. Existing cosmogenic radionuclide 10Be surface exposure ages （~130 ka） show that the Guxiang Glaciation occurred during marine oxygen isotope stage 6， however， key information， such as ice volume， ice extent， and equilibrium line altitude （ELA） in the Bodui Zangbo River catchment during the Guxiang Glaciation， is an open question. In this study， we reconstruct the ice extent of the Guxiang Glaciation for the whole catchment on the basis of repeated field investigations， existing numerical dating results of moraines， and high-definition remote sensing images/satellite photos； estimate the ice volume and ice thickness using a glacier surface profile model run in Excel and ArcGIS software； and estimate the ELA of the catchment during the Guxiang Glaciation using the accumulation area ratio （AAR） and toe-to-headwall altitude ratio （THAR） methods. Our results show that during the Guxiang Glaciation， the total glacier area was 2 648 km²， the average glacier thickness was nearly 360 m， and the ice volume around 953 km³. The ice-cover rate has decreased from 63% in the Guxiang Glaciation to the present cover rate of 22.4%. Based on the paleo-glacier surface elevation interpolated with the inverse distance weighting method， the Guxiang Glaciation ELA is estimated to be 3 871 m， which is 584 m lower than the present ELA of 4 455 m. Our ELA estimation result is consistent with existing values.

Keywords： Guxiang Glaciation ; Bodui Zangbo River ; glacier reconstruction ; equilibrium-line altitude （ELA） ; accumulation area ratio （AAR） ; toe-to-headwall altitude ratio （THAR）

PDF (8899KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

覃长雄, 许刘兵, 谢金明, 周尚哲. 藏东南波堆藏布江流域古乡冰期冰川重建[J]. 冰川冻土, 2021, 43(3): 724-734 doi:10.7522/j.issn.1000-0240.2021.0014

QIN Changxiong, XU Liubing, XIE Jinming, ZHOU Shangzhe. Glacier reconstruction of the Guxiang Glaciation in the Bodui Zangbo Valley in southeastern Tibet[J]. Journal of Glaciology and Geocryology, 2021, 43(3): 724-734 doi:10.7522/j.issn.1000-0240.2021.0014

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 古乡冰期冰川遗迹及其年代

波堆藏布江发源于念青唐古拉山东段南坡（29°89′~30°67′ N、95°08′~95°96′ E，图1），流域面积约1.46×104 km²，邻近雅鲁藏布江大拐弯，处于南亚夏季风进入高原的水汽通道上。流域内的高山（最高峰则普海拔6 364 m）迎风面降水丰富，使得该区域发育了规模庞大的第四纪和现代冰川，冰川ELA和冰川末端高度基本处于同纬度的最低值［2］。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   波堆藏布江流域位置及流域内现代冰川［4］和古乡冰期冰碛垄分布图

Fig.1   Map showing location of the Bodui Zangbo River catchment， modern glaciers［4］， and moraines of the Guxiang Glaciation within the catchment


古乡冰期指的是传统意义上的“倒数第二次冰期”，由李吉均等［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

新窗口打开| 下载原图ZIP| 生成PPT

图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

新窗口打开| 下载原图ZIP| 生成PPT

图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，https：//www.eorc.jaxa.jp/ALOS/en/aw3d30/index.htm）叠加，可得到古冰川的面积。

后运用冰面纵剖面模型模拟古乡冰期冰川的冰面形态和冰厚等特征参数。冰面纵剖面模型由Schilling等［6］提出，后经Benn等［7］改进，现已被广泛运用于各地（含青藏高原及周边山地）古冰川形态的重建［8-12］。

此模型基于的假设为冰川是理想的塑体，其形变由冰的重量及冰表面梯度引起的驱动应力（driving stress，τd）和冰川底部的屈服应力（yield stress，τy）的变化引起。当驱动应力大于屈服应力时，冰川发生运动；当驱动应力小于屈服应力时，就不会发生冰川的运动，但冰川会变厚从而增加驱动应力。为了保持冰川的稳定状态，即τy=τd，冰川的厚度会不断进行调整，从而引起驱动应力的变化。基于此得到了下列表达式。

τy=τd=ρgHdhdx

式中：ρ为冰川密度（900 kg∙m-3）；g为重力加速度（9.81 m∙s-2）；H为冰川厚度；h为冰表面高程；x为水平坐标（以冰川末端为原点，冰底和冰面高程为纵轴，距冰川末端的长度为横轴的坐标系）。有时用sinα代替dhdx，α为冰川表面坡度［13］。根据以上等式可求出沿着中流线的冰川厚度H，结果为抛物线。

Hx2=2τyρgx

但该式不适用于冰川坡度较大的地区。对于不规则的冰下地形，Benn等［7］指出必须沿着中流线按连续的点（编号1，2，…，i，i+1，i+2）计算，冰的表面坡度表示为

dhdx=(hi+1-hi)/Δx

式中：Δx为每个点的间隔距离，本文取Δx=200 m。

因此，由式（1）和式（3）得

hi+1=hi+(τyH)iΔxρg

式中：H=h-B，B为冰川底部高程，此式可以解出冰面高程h。

为了解决冰舌处冰厚为0的问题，van der Veen［14］进一步推导得到

hi+12-hi+1(Bi+Bi+1)+hi(Bi+1-Hi)-2Δxτ¯y/ρg=0

式中：τ¯y为相邻两点之间的平均屈服应力。该式是一个标准的二次方程，其标准模式为ay2+by+c=0，其中a=1，b=-（Bi+Bi+1），c=hi（Bi+1-Hi）-2Δxτ¯y/ρg。Benn等［7］将式（5）做成了一个基于Excel的计算冰川表面轮廓的ExcelTM表格程序。此程序只需要输入Δx、ρ、g、τy以及B就可以得到式（5）的解（h<B的解没有实际物理意义），即中流线上每个点的冰面高程以及冰厚数据。τy取值范围一般为50~150 kPa［13］，后续可用来检验该模型的可用性。

但该模型是以冰盖出发进行推导的，冰盖不受山谷两侧阻力的约束，其流动需要的驱动应力等于冰川底部剪切应力。而山谷冰川则需要考虑山谷两侧阻力的影响，因此引入形状因子f。

τB=τy=fτd

式中：τB为底部剪切应力。

在冰川中流线处，τB=τd，f为

f=AHp

式中：A为冰川横截面面积；H为冰厚度；p为横截面的“冰川周长”［15］；f的取值范围一般为0.4~1［8］。此时的驱动应力大于同等条件下冰盖流动所需要的驱动应力，在计算时，需要把式（4）中的τy用τy/f替换。

2.3　古冰川模拟过程

使用ArcGIS软件并结合DEM数据以及研究区的谷歌地球（Google Earth）影像图，在古冰川范围上提取出中流线，为了使中流线上插值点分布范围更全面，减小插值误差，人工选取了23条分布在面积较大的现代冰川上的中流线和12条分布在主谷位置的中流线（合计35条冰川中流线，图3）。在ArcGIS中，按等距离（Δx=200 m）对中流线进行分割，把分割得到的线段转化为点，最后提取每个点的高程数据，古冰川模拟流程如图4。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   古冰川模拟流程图

Fig.4   Flow chart of paleo-glacier simulation


在此ExcelTM表格程序中进行古冰川表面高程模拟时，步长（Δx=200 m）、冰的密度（ρ=900 kg∙m-3）、重力加速度（g=9.81 m∙s-2）是已知的，而B为冰底地形高程，存在现代冰川的地方需要减掉其厚度。输入的剪切应力τy一般取值100 kPa，根据谷歌地球影像图，可求得终碛-侧碛垄的高程数据作为“目标高程”，将其与中流线上相应点的高程对应起来。下面将以中流线E4（图3，古乡至关星冰川北部顶端）为例子，模拟结果见图5。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图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）


模拟过程中，从下游（古乡）至上游（关星）沿着中流线作出几道横剖面，计算每一段的形状因子［图5（d）］并插值，之后通过不断调整中流线上各点的剪切力，使冰面高程与“目标高程”（终碛-侧碛垄高度）接近一致［图5（c）］。

使用上述模拟过程，对其他34条中流线进行模拟，将模拟结果分别用反距离加权（inverse distance weight，IDW）、泛克里金（universal Kriging，UK）和普通克里金（ordinary Kriging，OK）插值法进行插值，最终得到古冰面高程，而古冰川厚度可由模拟古冰面高程减去古冰川冰底高程，结果见图6。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图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之间，且这个范围对于热带有积雪覆盖的冰川同样适用，地貌证据也表明了这个取值范围同样适用于古冰川。

本文以朱西冰川、白玉冰川、则普冰川、日母冰川和关星冰川为例，在已有冰川编目数据的基础上，计算这5条冰川AAR的平均值，并将得到的AAR平均值用来估算古乡冰期时冰川的ELA。具体操作如下：①根据第一次冰川编目［23］这5条冰川的积累区面积和总面积计算AAR值，取此AAR的平均值来估算古乡冰期时的冰川ELA；②根据得到的AAR值和第二次编目的冰川范围来计算各冰川的ELA（表1），其平均值用于下文中选取THAR值。考虑到这些冰川上部宽、下部窄，粒雪盆（积累区）与冰舌（消融区）都发育正常［24-25］，因此用这5条冰川AAR的平均值（0.73）来恢复古乡冰期时的冰川ELA。

表1   第一次冰川编目［23］下波堆藏布江流域内5条现代冰川的AAR值及估算的第二次冰川编目［4］下各ELA值

Table 1  The five modern glaciers AAR values in the Bodui Zangbo River catchment in glacier inventory of China［23］ and their estimated ELA values in second glacier inventory of China［4］

冰川名称积累区面积/km²总面积/km²第一次冰川编目AAR第二次冰川编目ELA/m
平均值0.734 451
朱西冰川7.0812.90.554 011
白玉冰川16.7524.40.694 048
则普冰川59.4472.40.804 390
日母冰川12.5215.20.804 849
关星冰川37.6247.60.794 956

新窗口打开| 下载CSV


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 estimates of five modern glaciers under different THAR values in the Bodui Zangbo River catchment

冰川名称ELA/m
THAR=0.30THAR=0.34THAR=0.35THAR=0.36THAR=0.40
平均值4 3454 4354 4584 4804 570
朱西冰川3 9234 0274 0544 0804 184
白玉冰川3 9964 0874 1104 1334 225
则普冰川4 2464 3624 3924 4214 538
日母冰川4 7444 8144 8314 8484 917
关星冰川4 8174 8844 9014 9184 985

新窗口打开| 下载CSV


3 结果与讨论

3.1　冰川规模变化

从中流线E4上的剪切力分布［图5（a）］发现，下游（距冰川末端的长度为0~45 km）剪切力变化较小，中游（距离冰川末端的长度为45~75 km）剪切力变化较大，上游（距离冰川末端的长度为75~120 km）剪切力较小，总体分布在7~120 kPa之间。经过计算，剪切力τy平均值为61.3 kPa，平均形状因子为0.75，等效平均剪切力τy/f为81.7 kPa，在剪切力合理取值范围内（50~150 kPa）［13］。图5（b）是冰面高程与“目标高程”相吻合时的冰川厚度分布，由图5（c）中冰面高程与冰底地形高程相减所得。可以从图5（b）和5（c）看出，从下游到上游，冰川厚度总体上呈增大的趋势，其最大值分布在上游，但临近冰川陡峭的后壁时厚度迅速减小。另外，下游地形较为平缓而中上游地形起伏较大。一般来说，地形平缓的下游更容易积累冰川，但在下游主谷内冰川厚度较小，从白玉村往上游［图5（c），距离末端50 km以上］冰川厚度逐渐增大，则普村之后冰川厚度基本维持着很大的值，这种趋势在图6（d）中也很好地体现了出来。从图6（a）~（c）中可以得知当时的ELA大致分布在则普村（图1则普冰川末端）附近，因此冰厚的这种分布是合理的。

通过对其余34条中流线的模拟，可以重建古乡冰期时古冰川的表面高程和冰川厚度（图6），大致估算出古乡冰期时波堆藏布江谷流域内冰川的形态参数，即流域内平均冰川厚度约360 m，最大厚度约893 m，面积约2 648 km²，体积约953 km³。古乡冰期时，流域内约63%的区域被冰川所覆盖，而根据第二次冰川编目统计［4］的数据，现代冰川的覆盖率约为22.4%。

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  ELAs in the Bodui Zangbo River catchment

冰川作用时间ELA/m下降值/m
AAR=0.73THAR=0.35平均值
现代4 451*4 4584 455

古乡冰期（基于IDW）

古乡冰期（基于UK）

4 021

3 867

3 720

3 649

3 871

3 758

584

697

古乡冰期（基于OK）3 8473 6493 748707

注：*为第二次冰川编目ELA平均值。

新窗口打开| 下载CSV


不同的空间插值方法得到的古冰面高程对于重建古冰川冰量、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  Precision comparison of paleo-glacier surface elevation estimated from different spatial interpolation methods

编号实测高程/m高程插值/m绝对误差/m相对误差/%
IDWUKOKIDWUKOKIDWUKOK
12 678.7422 678.0723 025.9292 956.5660.67347.19277.820.025-12.961-10.371
22 705.9312 706.7102 923.1672 882.9670.78217.24177.04-0.029-8.028-6.543
32 718.1232 719.0802 834.1482 811.3060.96116.0393.18-0.035-4.269-3.428
42 734.3592 733.1842 911.7022 881.4331.18177.34147.070.043-6.486-5.379
52 753.0432 753.5312 753.0872 753.0870.490.040.05-0.018-0.002-0.002
62 776.8282 776.5002 980.9472 946.6290.33204.12169.800.012-7.351-6.115
72 783.3992 783.7112 971.0172 940.6390.31187.62157.24-0.011-6.741-5.649
82 791.6362 790.3932 871.9722 859.9991.2480.3468.360.044-2.878-2.449
92 813.2752 813.4932 920.0512 909.3090.22106.7896.03-0.007-3.796-3.413
103 889.4263 885.1173 991.4503 989.1364.31102.0299.710.111-2.623-2.564

新窗口打开| 下载CSV


最终选择IDW插值得到的古冰面高程，结合AAR法和THAR法计算古乡冰期时波堆藏布江流域内冰川ELA，结果取平均值3 871 m（图6和表3），相较于现代冰川ELA（4 455 m）下降了约584 m（表3）。Zhou等［5］曾对古乡冰期时波堆藏布江流域内冰川ELA进行估算，结果显示，古乡冰期时冰川ELA约为3 800 m，相较于现代冰川ELA（4 600 m）下降约800 m，本文古乡冰期时冰川ELA的平均值与此前研究之间仅相差约71 m，下降值相差约216 m。

4 结论

本文在冰川地貌外野实地考察和已有年代学为基础上，辅以高清遥感影像/卫星照片，对波堆藏布江流域古乡冰期时的冰川作用范围进行了重建；后进一步运用冰面纵剖面模型对波堆藏布江流域古乡冰期时古冰川进行了模拟；此外运用AAR和THAR方法对古乡冰期和现代冰川ELA进行了估算。古乡冰期时，波堆藏布流域冰川的平均厚度约为360 m，冰川总面积约为2 648 km²，体积约为953 km³，流域内的冰川覆盖率由古乡冰期的63%缩减至现在的22.4%；现代冰川ELA约4 455 m，古乡冰期时ELA约为3 871 m，较现代冰川下降了约584 m。

“古乡冰期”是李吉均先生主导提出的第四纪冰川重要冰期概念，谨以此文，纪念李吉均院士！

参考文献

Shi Yafeng， Cui Zhijiu， Su Zhen. The Quaternary glaciations and environmental variations in China［M］. Shijiazhuang： Hebei Science and Technology Publishing House， 2006.

[本文引用: 1]

施雅风， 崔之久， 苏珍. 中国第四纪冰川与环境变化［M］. 石家庄： 河北科学技术出版社， 2006.

[本文引用: 1]

Li Jijun， Zheng Benxing， Yang Xijin， et al. Glaciers of Xizang （Tibet）［M］. Beijing： Science Press， 1986： 130-140.

[本文引用: 4]

李吉均， 郑本兴， 杨锡金， 等. 西藏冰川［M］. 北京： 科学出版社， 1986： 130-140.

[本文引用: 4]

Zhou Shangzhe， Xu Liubing， Colgan P M， et al.

Cosmogenic 10Be dating of Guxiang and Baiyu Glaciation

［J］. Chinese Science Bulletin， 2007， 52（10）： 1387-1393.

[本文引用: 2]

周尚哲， 许刘兵， Colgan P M， 等.

古乡冰期和白玉冰期的宇宙成因核素10Be定年

［J］. 科学通报， 2007， 52（8）： 945-950.

[本文引用: 2]

Guo Wanqin， Xu Junli， Liu Shiyin， et al.

The second glacier inventory dataset of China （version

1.0）［DS］. Lanzhou： Cold and Arid Regions Science Data Center， 2014.

[本文引用: 6]

郭万钦， 许君利， 刘时银， 等.

中国第二次冰川编目数据集（V

1.0）［DS］. 兰州： 寒区旱区科学数据中心， 2014.

[本文引用: 6]

Zhou Shangzhe， Wang Jie， Xu Liubing， et al.

Glacial advances in southeastern Tibet during late Quaternary and their implications for climatic changes

［J］. Quaternary International， 2010， 218（1）： 58-66.

[本文引用: 2]

Schilling D H， Hollin J T.

Numerical reconstructions of valley glaciers and small ice caps

［M］// Denton G H， Hughes T J. The last great ice sheets. New York： Wiley， 1981： 207-220.

[本文引用: 1]

Benn D I， Hulton N R J.

An ExcelTM spreadsheet program for reconstructing the surface profile of former mountain glaciers and ice caps

［J］. Computers and Geosciences， 2010， 36（5）： 605-610.

[本文引用: 3]

Jin Yao.

Reconstruction of Last Glacier Maximum paleo-glacier in the Qugaqie Valley in Nyainqentanglha Range

［D］. Beijing： University of Chinese Academy of Sciences， 2013.

[本文引用: 2]

金耀.

念青唐古拉山曲嘎切末次冰盛期古冰川重建

［D］. 北京： 中国科学院大学， 2013.

[本文引用: 2]

Dong Guocheng， Xu Xiangke， Zhou Weijian， et al.

Cosmogenic 10Be surface exposure dating and glacier reconstruction for the Last Glacial Maximum in the Quemuqu Valley， western Nyainqentanglha Mountains， south Tibet

［J］. Journal of Quaternary Science， 2017， 32（5）： 639-652.

Liu Jinhua， Yi Chaolu， Li Yingkui. Reconstruction of the Neoglacial glacier in the Qiangyong Valley， Mt.

Kaluxung， south Tibet

［J］. Quaternary Sciences， 2018， 32（2）： 348-354.

刘金花， 易朝路， 李英奎.

藏南卡鲁雄峰枪勇冰川新冰期冰川发育探讨

［J］. 第四纪研究， 2018， 32（2）： 348-354.

Wang Youqi， Zhou Weijian， Wu Yubin.

Reconstruction past glacier since MIS

3 and estimate its area， volume and equilibrium-line altitudes in Gaerqiong Valley， Lunpo Gangri central Gangdese［J］. Journal of Earth Environment， 2020， 11（2）： 169-180.

王友琪， 周卫健， 毋宇斌.

冈底斯中段冷布岗日嘎尔琼沟MIS3以来冰川重建及其平衡线高程估算

［J］. 地球环境学报， 2020， 11（2）： 169-180.

Pellitero R， Rea B R， Spagnolo M， et al.

A GIS tool for automatic calculation of glacier equilibrium-line altitudes

［J］. Computers and Geosciences， 2015， 82： 55-62.

[本文引用: 2]

Nye J F.

The mechanics of glacier flow

［M］. 2nd ed. Journal of Glaciology， 1952， 2（12）： 82-93.

[本文引用: 3]

van der Veen C J.

Fundamentals of glacier dynamics

［M］. Rotterdam， the Netherlands： Balkema， 1999： 462.

[本文引用: 1]

Paterson W S B. The physics of glaciers［M］. 3rd ed. London： Butterworth-Heinemann， 1994： 480.

[本文引用: 1]

James W H M， Carrivick J L.

Automated modelling of spatially-distributed glacier ice thickness and volume

［J］. Computers and Geosciences， 2016， 92： 90-103.

Meier M F.

Proposed definitions for glacier mass budget terms

［J］. Journal of Glaciology， 1962， 4（33）： 252-263.

[本文引用: 1]

Meier M F， Post A S.

Recent variations in mass net budgets of glaciers in western North America

［C］// Variations of the regime of existing glaciers： Symposium of Obergurgl 1962. Wallingford， UK： International Association of Hydrological Sciences， 1962： 63-77.

[本文引用: 1]

Porter S T.

Equilibrium line altitude of late Quaternary glaciers in the Southern Alps， New Zealand

［J］. Quaternary Research， 1975， 5： 27-47.

[本文引用: 1]

Kaser G， Osmaston H. Tropical glaciers［M］. Cambridge， UK： Cambridge University Press， 2002： 149-192.

[本文引用: 1]

Kern Z， Laszlo P.

Size specific steady-state accumulation-area ratio： an improvement for equilibrium-line estimation of small paleo-glaciers

［J］. Quaternary Science Reviews， 2010， 29： 2781-2787.

[本文引用: 1]

Kaser G， Osmaston H. Tropical glaciers［M］. Cambridge， UK： Cambridge University Press， 2002： 207.

[本文引用: 1]

Mi Desheng， Xie Zichu， Luo Xiangrui. Glacier inventory of China XI： the Ganga drainage basin and Indus River drainage basin［M］. Xi’an： Xi’an Cartographic Publishing House， 2002.

[本文引用: 3]

米德生， 谢自楚， 罗祥瑞. 中国冰川编目XI： 恒河水系和印度河水系［M］. 西安： 西安地图出版社， 2002.

[本文引用: 3]

Wang Xiaoli.

Glacial and climate change in Boduizangbu Valley in Baiyu ice stage

［D］. Lanzhou： Lanzhou University， 2006.

[本文引用: 2]

王孝理.

波堆藏布谷地白玉冰期冰川与气候变化

［D］. 兰州： 兰州大学， 2006.

[本文引用: 2]

Wang Jie.

Study on Quaternary glacial evolution in the southeastern Tibet and cosmogenic nuclides dating

［D］. Lanzhou： Lanzhou University， 2007.

[本文引用: 1]

王杰.

西藏东南部第四纪冰川演化序列与宇宙成因核素测年研究

［D］. 兰州： 兰州大学， 2007.

[本文引用: 1]

Cui Hang， Wang Jie.

The methods for estimating the equilibrium line altitudes of a glacier

［J］. Journal of Glaciology and Geocryology， 2013， 35（2）： 345-354.

[本文引用: 2]

崔航， 王杰.

冰川物质平衡线的估算方法

［J］. 冰川冻土， 2013， 35（2）： 345-354.

[本文引用: 2]

Cui Hang， Cao Guangchao， Chen Kelong， et al.

A review on the reconstruction of paleoclimate by glacial landform in the Tibetan Plateau and adjacent mountains

［J］. Journal of Glaciology and Geocryology， 2021， 43（1）： 254-262.

[本文引用: 1]

崔航， 曹广超， 陈克龙， 等.

青藏高原及毗邻山地利用冰川地貌重建古气候的研究综述

［J］. 冰川冻土， 2021， 43（1）： 254-262.

[本文引用: 1]

Charlesworth J K. The Quaternary era［M］. London： Arnold， 1957.

[本文引用: 1]

Manley G.

The late-glacial climate of north-west England

［J］. Liverpool and Manchester Geological Journal， 1959， 2： 188-215.

[本文引用: 1]

Pewe T L， Reger R D.

Modern and Wisconsinan snowlines in Alaska

［C］// Proceedings of the 24th International Geological Congress.

[本文引用: 1]

［S.

l

.： s.n.］， 1972： 187-197.

[本文引用: 1]

Andrews J T. Glacial systems： an approach to glaciers and their environments［M］. North Scituate， MA， USA： Duxbury Press， 1975.

Meierding T C.

Late Pleistocene glacial equilibrium-line in the Colorado Front Range： a comparison of methods

［J］. Quaternary Research， 1982， 18： 289-310.

[本文引用: 1]

Benn D I， Owen L A.

Himalayan glacial sedimentary environments： a framework for reconstructing and dating the former extent of glaciers in high mountains

［J］. Quaternary International， 2002， 97： 3-25.

[本文引用: 1]

Ramge J M， Smith J A， Rodbell D T.

Comparing reconstructed Pleistocene equilibrium-line altitudes in the tropical Andes of central Peru

［J］. Journal of Quaternary Science， 2005， 20： 777-788.

[本文引用: 1]

/

〈 〉

陇ICP备2021001824号-10
　 甘公网安备 62010202000676号
版权所有 © 《冰川冻土》编辑部
地址：兰州市天水中路8号　
邮编：730000　电话：0931-8260767　
E-mail：edjgg@lzb.ac.cn
总访问量 　 今日访问 　 在线人数
技术支持：北京玛格泰克科技发展有限公司