青藏高原多年冻土变化对水文过程的影响

赵 林 胡国杰 邹德富 吴晓东 马 露,3 孙 哲,3 原黎明,3 周华云,3 刘世博,3

1 南京信息工程大学 地理科学学院 南京 210044 2 中国科学院西北生态环境资源研究院 冰冻科学国家重点实验室 青藏高原冰冻圈观测研究站 兰州 730000 3 中国科学院大学 北京 100049

青藏高原及周边地区孕育着黄河、长江、恒河、澜沧江、印度河、萨尔温江和伊洛瓦底江等亚洲的重要河流，被称为“亚洲水塔”。其广泛分布的多年冻土通过独特的水分运移影响着区域水文和水循环过程。青藏高原现存多年冻土的总面积约为 106 万平方公里[1]。多年冻土是在青藏高原隆升过程中逐渐形成并扩张的，并经历了不同周期和尺度的气候冷暖波动。在多年冻土层形成的漫长过程中，反复的成冰作用将大量水分冻结并以固态形式长期储存于地下而形成了地下冰，现有地下冰储量约为 12.7 万亿立方米。

在气候变暖背景下，青藏高原多年冻土发生着广泛的退化[2-7]。最显著的特点是多年冻土温度升高和活动层增厚，浅表层多年冻土及其中的地下冰逐渐融化。活动层增厚引起的地表水分条件的改变是影响多年冻土区产汇流过程和生态过程的重要因素；而曾被固存并埋藏在多年冻土内部的地下冰的融化，则会导致更多的水分被释放并参与水分的年循环过程，同时还会引起地面沉降，对高原多年冻土区的水循环过程和气候产生影响[8-11]。多年冻土区地表的水热过程和变化在区域水文过程中也发挥着重要作用。与积雪、冰川和湖沼相比，冻融过程中的水分迁移和地下冰融化的产汇流过程对多年冻土区水循环过程的影响更为复杂。因此，明确青藏高原多年冻土变化及其影响对于认识多年冻土与区域水循环，以及生态和气候相互作用的关系具有重要意义。

多年冻土层的透水性能较差，作为一种大范围的区域性隔水层或弱透水层，在一定时空尺度上阻隔或显著减弱了大气降水、地表水同地下水之间的水力联系。因此，青藏高原大范围存在的多年冻土强烈影响着地表径流形成以及地下水的运移过程和分布格局。多年冻土的隔水作用可以提高流域融雪和降雨径流的产流量，而多年冻土退化会直接影响寒区地下水补给源和补给量、径流路径和排泄过程，以及地下水与地表水的交换等[14-19]。多年冻土的不均匀融化可导致多年冻土区冻土层的区域性稳定隔水作用不断减弱，冻结层上水水位随之下降，补给路径延长、加深，甚至可通过新形成局部融化的“天窗”（贯穿融区）直接补给冻结层下水或冻结层间水[20]，这对于局地到流域尺度的地下水循环，特别是地下水补、径、排过程都产生了深刻影响。多年冻土区的径流系数一般比非冻土区大，并随气候变暖、活动层增厚而明显减小。加拿大北极地区一个多年冻土小流域的观测结果显示，在某些气温较高、降水较多年份的实测径流量反而比正常年偏小很多[21]。此外，多年冻土退化及活动层增厚必将导致多年冻土中大量赋存的地下冰转化为液态水；被释放的液态水部分参与活动层的冻融过程，而其余部分将参与区域乃至全球的水循环过程[22]。因此，加深对多年冻土特征和变化的理解，可为进一步认识区域水循环过程提供理论基础。

很多报道表明，活动层增厚、地下冰融化已经引发了多年冻土区大量的热融现象，其中对局地水文循环过程影响极大的热融现象导致了热融湖塘的形成和扩展。在多年冻土退化初期，地下冰融化导致地表下陷形成积水洼地；积水洼地一旦出现，原有地表的水热平衡即被打破，多年冻土的融化也就不断加剧，热融洼地也随之不断扩张而形成热融湖塘[23]；当湖底的多年冻土全部融化之后，冻结层上水与冻结层下水贯通而形成贯通融区，湖水可能经由这个贯穿融区被迅速排干，之后在被疏干的湖盆底部开始形成新的多年冻土。可以看出，热融湖塘的形成—发展—消亡过程与地表水与地下水的转换直接相关，极大地影响局地的水文条件[24]。

基于前述认识，本文在长期观测数据的基础上，结合再分析数据产品，试图通过对青藏高原水文过程密切相关的多年冻土区气温、降水、地温和土壤含水量、地表融沉等特征基础资料的分析和综述，为“亚洲水塔”的变化及机理研究添砖垫瓦。

1 多年冻土变化的基本特征

1.1 多年冻土区气温与降水变化

当前大部分研究基于青藏高原的各类数据资源揭示了青藏高原不同区域乃至整个青藏高原的气候变化背景[25-29]，但由于真正位于多年冻土区的气象台站极为稀少，导致现有结果对于高原海拔 4 000 m 以上地区气候变化背景的评估还存在极大不确定性。2004 年以来，中国科学院青藏高原冰冻圈观测研究站在青藏高原腹地多年冻土区沿青藏公路沿线的唐古拉、五道梁和西大滩 3 地建立了综合气象站，重点开展高原腹地多年冻土与气候相互关系监测研究。这 3 个气象站 2004—2016 年的年平均气温在 -5.8℃—-2.4℃ 间波动上升，气温升高趋势基本一致（图 1a），平均变化率约为0.05℃/a。其中，唐古拉气温增幅最大，约为 0.08℃/a；五道梁气温增幅最小，约为 0.02℃/a。2008 年的气温为观测期的最低值，平均值约为 -5℃；2010 年为观测期最高值，平均值约为 -3.3℃（图 1c）。

上述 3 个气象观测站点的年降水变化趋势基本一致，约在 210—580 mm 降水量区间呈明显的波动上升趋势（图 1b），平均变化速率为 7.49 mm/a。年平均降水量最小值和最大值出现在 2006 年和 2016 年，分别为 285 mm 和 436 mm（图 1c）。降水主要集中在每年5—9月，但年际降水集中期的分布有所差异。观测站点的降水量空间变化趋势与气温变化趋势相反，受纬度地带性规律和海拔的共同影响，降水量自北向南迅速增加。其中，唐古拉降水增幅最大，约为 10.68 mm/a；五道梁降水增幅最小，约为 4.81 mm/a。

图1 2004—2016 年唐古拉、五道梁和西大滩观测站点气象要素的年变化特征

可以看出，青藏公路沿线多年冻土区的气候在 10 多年来呈现变暖和变湿趋势。利用青藏高原气象站点和再分析资料研究结果表明，青藏高原自 1980 年以来气温变化速率在 (0.40℃—0.52℃)/10 a 之间，1998—2012 年的升温率为 0.20℃/10 a[29,30]。1979—2001年，青藏高原西北部降水增加速率为 3.99 mm/10 a，东南部增加速率约 16.84 mm/10 a[31,32]。总体而言，2004年以来，青藏高原腹地多年冻土区的气候变暖变湿趋势远较整个青藏高原的平均气候变化幅度要强烈。多年冻土区气候的这种变化无疑会导致多年冻土及多年冻土水资源特征发生显著变化。

1.2 活动层温度和含水量变化特征

1998 年以来，本项目组陆续在青藏公路沿线公路两侧的天然无破坏地表下建立了 9 个活动层观测站点（昆仑山垭口至两道河段），其观测内容是自地表到多年冻土上限之间厚度约为 1—4 m 不同深度处的土壤温度、未冻水含量、土壤热通量等。观测结果显示，这些观测站点的活动层厚度在 100—400 cm 之间，活动层厚度自 2004—2018 年呈现显著增加趋势（图 2a），平均变化率达到 21.7 cm/10 a（图 2b）。活动层底部温度变化范围在-3.2℃—0oC之间，从2004—2018 年呈现出明显的升温趋势，平均升温率为0.45oC/10 a（图 3）。

在施工期间应做好新旧路基的沉降过渡协调工作，保证新旧路基结合部分的协调控制，在提升沉降过渡控制效果的情况下，编制完善的计划方案，提升整体工作效果，达到预期的管理目的[7]。

大量研究表明，青藏高原多年冻土区活动层的融化过程开始于每年 5 月中下旬，结束于每年 9 月下旬—11月初[6,33-35]。受气温季节变化的影响，每年 9 月之后气候逐渐变冷，活动层自上而下的融化过程变得非常缓慢。因此，本文选取青藏公路沿线 5 个活动层观测站点 2004—2017 年不同深度土壤水分观测资料进行统计分析，探讨 2004 年以来活动层内部的水分变化。结果显示，活动层水分含量的年平均值和每年 9 月平均值均表现出自地表向下逐渐增大的趋势，活动层底部的土壤含水量较高。活动层冻结期土壤水分大多呈现冻结状态，土层中只含有较少量的未冻水，各站点 9 月平均土壤含水量介于 2.5%—39.7%，年平均值介于 2.1%—25.8%，年平均值明显低于 9 月平均值；5 个观测点活动层底部 9 月多年平均土壤含水量介于 21.7%—29.8%。通过趋势分析发现，2004—2018年，所有站点活动层底部的土壤含水量都呈现显著增加趋势（图 4），增加量为 11%—32%，平均每年增加 2% 以上；而活动层上部的土壤含水量却表现出降低趋势或者是没有显著变化，观测期内地表土壤含水量的降低幅度一般小于 2%。

图2 青藏公路沿线各观测点活动层厚度2004—2018 年变化趋势

图3 青藏公路沿线2004—2018 年各观测点活动层底部温度及变化趋势

图4 青藏高原不同深度土壤水分在2004—2017 年间变化特征（每年9 月平均土壤含水量）

研究表明，青藏高原多年冻土退化、活动层增厚会导致表层土壤含水量减少[36]。上述结果显示，2004年以来，青藏公路沿线活动层内部的水分状况整体呈现下部增大，表层基本不变或略有降低的趋势，这与同期降水量呈显著增加趋势相吻合。换言之，降水量的增加一定程度缓减了地表土壤含水量的减少趋势。

2 多年冻土区降水和土壤含水量时空变化趋势

为更好地探求青藏高原多年冻土变化的水资源效应，亟待查明多年冻土区空间尺度的降水量和土壤储水量之间的协同变化特征。本文选取中国气象信息中心中国地面降水 0.5°×0.5° 格点数据集和第五代全球大气再分析资料（ERA5）土壤水分再分析数据产品，以探求青藏高原多年冻土区的降水和土壤含水量空间变化格局。

2.1 多年冻土区降水时空变化趋势

中国气象信息中心中国地面降水 0.5°×0.5° 格点数据集是相对较为可信的一套数据，对青藏高原降水量空间插值，尤其是对降水量变化的描述较好。对该数据的分析表明，青藏高原多年冻土区的年平均降水量存在明显的区域差异，自西北向东南从几十至上千毫米递增；其中西北部降水量偏低，东南部降水量较大（图 5a）。从年际变化来看，1980—2018 年，除青藏高原东南部区域的年降水量呈明显减少趋势外，青藏高原大部分区域的降水量呈明显增加趋势，中西部尤为明显，增加速率范围在 2.0—4.4 mm/a，降水变化速率最大值约出现在唐古拉山南坡附近。其中，比较干旱的西北部多年冻土区年降水量在 99—400 mm 之间，但该区域的年降水变化速率最大（1.5—3.3 mm/a）；而在东南部降水量较大的部分区域呈现微弱的减小趋势（图 5b）。总体上，青藏高原多年冻土降水增加率为29.5 mm/10 a。

2.2 多年冻土区土壤含水量时空变化趋势

ERA5 土壤水分再分析数据产品在青藏高原具有较好的表现[37]。本文利用该数据分析了青藏高原多年冻土区 1980—2018 年 6— 8月0—100 cm 土壤含水量均值变化特征。从图 6a 中可以看出，多年冻土区0—100 cm 土壤含水量主要分布在 0.1—0.4 m3/m3，其中青藏高原东南部土壤含水量较高，最大值为 0.69 m3/m3；而在青藏高原西南部和腹地土壤含水量较低，最小值为 0.06 m3/m3，多年冻土区土壤含水量平均值为 0.30 m3/m3。从变化趋势上来看，如图 6b 所示，多年冻土区土壤含水量整体呈现增加趋势，每 10 年增加率为 0.005 m3/m3；从区域来看，在青藏高原东北部土壤含水量较高的区域有略微减小的趋势，每 10 年变化率为 -0.01—0 m3/m3；而在青藏高原西南部和腹地大部分区域呈现增加趋势，每 10 年变化率为 0.01—0.03 m3/m3。可以看出，土壤含水量较低的区域增加趋势明显，而在土壤含水量较高的区域呈现减少的趋势，这与前述站点监测结果非常一致。

图5 1980—2018 年青藏高原多年冻土区多年降水量变化的空间分布

图6 青藏高原多年冻土区1980—2018 年土壤含水量变化（6—8 月）

3 多年冻土变化的影响

3.1 多年冻土变化对区域水文过程的影响

从前述结果可以看出，2004 年以来，青藏高原多年冻土区的气温在升高，降水量在增加，这种变化幅度远大于青藏高原气候变化的平均幅度，也大于 1980 年以来青藏高原气候整体的变化幅度。在这样的气候变化背景下，多年冻土区活动层厚度在增加，温度在升高，而土壤含水量整体在增加[38-41]（图 4）。

土壤含水量的增大与多年冻土的隔水作用密切相关。一方面，在冻结过程中未冻水的迁移导致更多的未冻水向多年冻土上限迁移。而活动层在融化过程中，融雪、降水和土壤水分在重力、温度和渗透梯度共同作用下迁移到多年冻土上限之上，导致土壤水分聚集在活动层底部[42]。另一方面，因为活动层厚度增加，导致底部厚层地下冰融化释放出大量冻结水。以可可西里观测站点（QT01）为例，在观测期内（2004—2017 年），活动层厚度从 160 cm 增加到 176 cm。例如，在 2004—2006 年观测深度 165 cm处尚处于冻结状态（图 2），随着多年冻土融化，地下冰释放大量液态水，导致活动层底部含水量增加。研究表明，活动层加深、融化期延长增加的多年冻土层上水的冬季补给也不可小觑；活动层增厚也可以通过增加储水空间而减少夏季径流[43]。多年冻土退化会导致蒸发、径流、水位等水循环要素变化[44-46]；并对径流、水位和水力联系有直接影响[47,48]。由于多年冻土的弱透水性，活动层底部增加的这部分冻结层上水有可能通过侧向流在低洼处形成径流，进而增加青藏高原多年冻土区江河湖泊的径流量。

图7 青藏高原活动层储水量与降水变化趋势

此外，冻土退化会导致流域内更多的地表水入渗变成地下水，造成流域地下水储水量加大[49]，导致冬季径流增加[50,51]。但是通过分析唐古拉观测点（QT04）和西大滩观测站点（QT09）站点活动层储水量和降水量的变化趋势发现，降水量显著增加，而活动层储水量却没有明显变化（图 7）。这除了与表层强烈的蒸发和侧向流有关，其原因是否与青藏高原多年冻土区径流在增加有关还需要进一步研究。高纬度地区的研究表明，地下冰融化、活动层加深和融化期延长可能会导致冬季基流和暖季径流显著增加[52,53]；在不连续冻土区，地下冰融化显著增加了冬季的河川基流，冻土退化使得径流的季节分布更加平缓[54-56]；多年冻土退化对径流分配的影响与多年冻土覆盖率密切相关，只有在多年冻土覆盖率高的流域，多年冻土退化才会引起产流过程的较大变化，而对于多年冻土低覆盖率流域，多年冻土退化的影响则较小[57]。然而，多年冻土对气候变暖的响应是一个缓慢的过程，因此冻土退化对水文过程的影响也是一个渐变过程。随着青藏高原多年冻土区的变暖变湿（图 5 和 6），多年冻土的退化对水文过程影响的程度究竟有多大，具体的影响过程如何还需要进一步探究。

3.2 多年冻土变化对湖泊面积变化可能的影响

1990—2010 年青藏高原湖泊面积增加了26%[58,59]。湖泊水位升高和水量增大现象归因于区域降水量增加[60]。黄河源区湖泊扩张、数量增多的主要原因正是同期降水量增大、蒸发量减少和冻土退化加大了地下冰融化水补给量的直接反应[61]；班戈错湖泊水位持续上涨，湖面扩大的原因之一是由于该流域上游在海拔约 4 800 m 部位存在一片面积较大的多年冻土湿地，而该多年冻土区的边缘有热融湖塘发育。近 20 年来该地区的多年冻土的退化及其中部分地下冰的融化可能是导致班戈错湖泊水位上升的主要原因之一[62]。

相关研究表明，“亚洲水塔”主要河川径流和湖泊短期内可能增加[63]，并且这一趋势可能将在未来（21 世纪中期）得以持续[64]，其中印度河、恒河、雅鲁藏布江、澜沧江上游等径流增加的幅度相对较大[65]。但随着冰川规模不断萎缩，冰雪融水的贡献率将会不断下降；而多年冻土区活动层增厚，土壤蓄水容量将增大，从而导致降雨下渗量增加和基流量增大[66,67]；多年冻土退化时地下冰的融化对湖泊有一定的补给作用，但是对湖泊的补给作用的大小取决于很多因素，如气候变化的方式、幅度、冻土融化速率、深度，以及局地岩性、排水和水文地质条件等[68]。除了多年冻土区湖泊上升之外，也有研究认为区域地下水位下降的主因是冻土退化，其导致冻结层上水水位持续下降甚至消失，从而可能导致区域产流减少[69-72]。因此，多年冻土退化和水循环的相互作用需要进一步研究。

3.3 多年冻土变化对地下冰释放量及地表形变的影响

青藏高原多年冻土年均地下冰融化量总计可达 12.7×103km3水量[73]，活动层加深 25 cm 引起的地下冰融水量相当于欧亚大陆北部主要河川径流的增加量。青藏公路沿线多年冻土上限附近地下冰主要由其上部的活动层水进行补给（59%—87%），其次是下部的多年冻土水（13%—41%）[74]；另外，下垫面的差异会显著影响活动层水对上限附近地下冰的补给比例。比如，高寒草甸区域的活动层水对其下部上限附近地下冰的补给比例（59%—69%）小于高寒草原所在区域的活动层水补给比例（70%—87%）。

降水作为活动层水的主要补给源，对地下冰的发育和形成也有着间接的贡献[75]，地表径流对浅层地下冰的形成也有显著的补给作用[76]。而针对多年冻土区地下冰的释放量到底去向何处，目前研究还较少。在北麓河热融湖塘地区，浅层地下冰融水对当地热融湖塘的补给比例最大能达到60%左右[77]；在昆仑山垭口地区，地下冰融水对地面径流也有着一定的贡献，浅层地下冰的融水对当地地表径流的补给比例达到 37.4%，仅次于冰川融水 56.7% 的补给比例[78]；在黄河源区，地下冰融水对当地地表径流的贡献比例为13.2%—16.7%[79]。本文的研究结果反映了活动层厚度增加导致的地下冰融化对活动层底部土壤含水量有显著的影响；然而，针对多年冻土退化对地下冰的释放量有多大影响，地下冰释放量到底对区域水循环的影响程度有多大，到目前为止还没有被量化。

气候变暖背景下，多年冻土退化还会造成多年冻土地面发生长期的沉降形变，这主要是由于进入多年冻土内部的能量过剩时多年冻土上限处地下冰融化所导致。研究表明，年间地表形变的高值区主要分布在地下冰含量较高的地区。例如，青藏高原的五道梁、北麓河地区和俄罗斯的勒拿河三角洲地区，多年冻土上限附近地下冰的融化导致了较高的年间沉降量。在五道梁地区年间沉降量可达 10.28 mm[80]。多年冻土的年间形变也有明显的空间差异性——北麓河高山地区多年冻土存在较为稳定，年平均形变量为 -1.78 mm[81]；而稀疏植被区域多年冻土地表形变有着明显的不同，其年形变量在 -16—0 mm 之间，退化更为明显[82]。

多年冻土地下冰含量较少的区域年间形变量较小，主要在青藏高原的西北部和西藏当雄县附近[83,84]。随着气候变暖的加剧，地下冰融化速率加快，多年冻土年间地表下沉的速率也随之增加。有研究表明，基于 ERS1/2 及 Envisat 数据获得的 1995—1 9 9 9 年格陵兰岛东北部地区的年平均沉降速率为 -0.3—-2.4 mm/a，而 2006—2009 年地表沉降速率增加到 -0.8—-2.7 mm/a，沉降速率可以很好地反映气候在变暖、多年冻土在退化[85]。因此，多年冻土退化导致的地面形变均与多年冻土内部地下冰直接相关，所以加强对地表形变的监测可为多年冻土区地下冰储量的变化提供可靠的依据。

4 结论

多年冻土退化给青藏高原及周边地区的水循环过程和水资源时空分布带来显著影响。具体而言，多年冻土层消融可释放水分补充地下水，从而影响区域水文过程。但这部分水分到底对区域水资源调节作用有多大尚待进一步深入研究。本文以观测事实为基础，结合再分析数据分析了青藏高原多年冻土变化及可能的影响，主要得到的 3 条结论。

（1）近 10 年来，青藏高原多年冻土区呈现变暖变湿的特征，多年冻土呈现显著的退化趋势。其地温明显升高、活动层增厚和活动层底部含水量增加，活动层底部温度变化率平均为 0.45oC/10 a，活动层厚度变化率达到了 21.7 cm/10 a。

（2）1980—2017 年青藏高原多年冻土区降水呈微弱的增加趋势，增加率为 29.5 mm/10 a。多年冻土区在1980—2018 年 6—8 月的 0—100 cm 土壤含水量平均值为 0.30 m3/m3；并且多年冻土区降水整体呈现增加趋势，每 10 年的增加率为 0.005 m3/m3。

（3）多年冻土退化将对青藏高原水文过程、湖泊面积变化等会有明显的影响，尤其是多年冻土退化导致的地下冰融化对于水循环会产生显著影响，并可能引起多年冻土地面发生长期的沉降形变；但多年冻土退化对水循环究竟影响的程度有多大需要进一步深入探索。