彭晨阳，盛 煜，吴吉春，曹 伟，何彬彬

（1.中国科学院西北生态环境资源研究院冻土工程国家重点实验室，甘肃兰州730000；2.中国科学院大学，北京100049）

0 引言

多年冻土是指持续两年或两年以上温度处于0℃及以下岩土层［1］，高海拔多年冻土泛指高纬度多年冻土南界以南、一定海拔高度以上出现的多年冻土［2］。我国高海拔多年冻土面积居世界首位，总面积约173.20×104km²，其中约80%分布在青藏高原［3］。冻土在世界气象组织（WMO）的全球观测系统中被列为6个指示全球气候变化的冰冻圈因子之一［4］，其主要特征是低温且大部分地区的土层中含有丰富的地下冰。青藏高原被称为“亚洲水塔”［5］，其中分布的冰川和冻土是东亚地区重要的水资源来源。近五十年来青藏高原区域增温速率为全球平均值的两倍［6］，区域整体气候正逐步向着暖湿化方向发展。冻土对气候变化与人类活动反映极其敏感［7］，在气候变暖和人类活动双重影响下，冰川、冻土呈快速退化态势［8］。尽管多年冻土对气候变化响应相较于冰川表现不是很明显，但由于其存在范围广，多年冻土的变化对当地生态环境影响正在逐渐凸显，特别是对水资源的影响有着不可忽视的作用。

祁连山区位于青藏高原的东北缘，横跨甘肃青海两省，是我国甘青两省与河西地区主要内陆河流的发源地，是“亚洲水塔”的重要组成部分，是我国重要的生态屏障区之一［9］。祁连山区内多年冻土分布广泛，属于典型的高海拔多年冻土［10］。相关研究最早可以追溯到20世纪70年代，由中国科学院在祁连山木里煤矿地区建立了冻土长年观测站［1］。早年关于祁连山多年冻土空间分布的研究主要包含在青藏高原整体多年冻土分布研究的子集中［11-14］，或对祁连山区局部地区的多年冻土下界进行调查，如：王绍令等［15］对祁连山东西段局部地区多年冻土调查结果显示，西段喀克图地区多年冻土下界在3 950～4 000 m，东段达坂山垭口附近的多年冻土下界在3 780～3 830 m，年均地温偏高，多年冻土层较薄［16］。近年来，有许多学者对祁连山局部地区的多年冻土进行了研究，如吴吉春等［10，17］对祁连山东部冻土特征调查结果显示，该区冻土受垂直地带性影响显著，冻土分布呈现季节冻土-不连续多年冻土-大片连续冻土-连续冻土逐渐过渡分布的模式，且多年冻土正处于退化状态；王庆峰等［18］对黑河上游多年冻土考察结果显示，该区多年冻土下界在3 650～3 700 m；王生廷等［19］基于钻孔数据，对大通河源区多年冻土特征及变化进行了分析，结果显示源区冻土平均地温随海拔变化梯度约3.82℃·km-1，整体处于退化状态，年平均地温上升速率为0.0075℃·a-1。Wang等［20］基于模型响应的研究发现，受气候变影响，近几十年来祁连山区多年冻土面积减少了约2.63×104km2。

目前关于祁连山区多年冻土的研究，大多情况下都是伴随着区域经济开发活动展开的，对多年冻土研究的区域有比较明显的限制，多数研究都集中在东部地区，中西部地区几乎没有进行过系统调查。掌握祁连山区多年冻土的空间分布范围，是研究多年冻土退化对区内脆弱的高寒生态环境影响的前提。在第二次青藏高原综合科考中，对祁连山区不同区域的多年冻土下界进行了较为系统的调查，本文综合此次调查结果，结合部分钻孔资料以及前人的研究的成果，对祁连山区多年冻土空间分布规律进行了分析与总结。

1 研究区概况

本文研究区的地理范围介于93°30′～103°03′E、35°50′～39°50′N之间，基于美国太空总署（NASA）和国防部国家测绘局（NIMA）联合测量的SRTM（Shuttle Radar Topography Mission）数据，数据转化为WGS1984坐标系后，重采样成100 m的空间分辨率，采用以此数据海拔3 000 m作为界定祁连山区范围的基带，以高于此海拔的祁连山区作为本文研究对象（图1），面积约16.90×104km2。

图1 祁连山区地形及地理位置示意图Fig.1 Map showing the topography and geographical position of the Qilian Mountains

祁连山区地处青藏高原东北前缘，因不断遭受构造隆升，逐渐成为高原的一部分［21］。祁连山区由一系列条块状、逆冲结构的造块体相互叠覆形成的NWW—SEE走向的7条大致平行的山脉和宽谷组成的，山脉东起于乌鞘岭，西止于当金山口与阿尔金山相接，南靠柴达木盆地，北临河西走廊［22］，东西长约1 000 km，南北宽约250 km，在一系列构造活动中谷地和盆地发育。

祁连山区水系丰富，西部海拔相对较高，大量冰碛堆积地貌，第四纪冰川发育；东部海拔相对较低，黄土堆积较厚［23］；河流以哈拉湖为中心，向四周辐射状外流。祁连山区东部水系主要是以黄河水系为主的大通河、湟水等；南部主要以流向柴达木盆地和青海湖的内陆河水系为主的巴音郭勒河、布哈河等；西部主要是党河、鱼卡河等；北部主要是疏勒河、北大河、黑河等。

祁连山区地处青藏高原、黄土高原和蒙新高原交汇地带，受大陆性气候和青藏高原气候综合影响，处于西风环流、青藏高原季风、东亚季风三个大气环流系统的交汇区［24］。张强等［25］对祁连山区空气水资源研究表明，祁连山区东西部降水差异较大，其中中部降水量最多，东部较少，东部高于西部，西部最少，平均每年水汽输送量为8.85×1010m³。近年来祁连山区中部降水明显增加［26］，贾文雄等［27］基于气象站的气温和降水资料的分析表明，祁连山区气候呈现暖湿化，与我国西北地区变化规律基本一致，对全球气候变化响应明显。曹广超等［28］基于祁连山区南坡气象站的气温数据表明，祁连山区多年平均气温呈上升趋势，增长率约0.35℃·（10a）-1。

2 数据与方法

目前我国对高海拔多年冻土的研究主要集中在青藏高原。研究方法主要是建立多年冻土空间分布与纬度、高程、平均地温、冻结指数等起决定性作用的宏观因素的统计关系。如吴青柏等［30］在对青藏公路沿线冻土分布研究中采用的年平均地温模型［29］、顶板温度模型；南卓铜等［31］在对青藏高原整体多年冻土分布研究中采用的地面冻结数模型。上述模型可以较好的在宏观上或局部地区模拟出多年冻土分布情况，但也存在一定的局限性，如年平均地温模型需要一定数量、较均匀分布的地温数据支持，顶板温度模型被证实在青藏高原不适用于-1.5℃以上的高温多年冻土分布模拟［30］。

祁连山位于多年冻土南界以南，其中的多年冻土属于高海拔多年冻土。本文尝试利用地带性规律对祁连山区多年冻土的空间分布规律进行探究。地带性规律是指程国栋提出的高海拔多年冻土空间分布遵从“三向地带性”规律，即由南北热量差异引起的纬度地带性［32］；由距海洋远近和大气环流造成的干燥度地带性；由水热条件随高度变化引起的垂直地带性。高海拔多年冻土与高纬度多年冻土分布相比，多年冻土分布虽然受纬度地带性、干燥度地带性影响，但主控因素是海拔，直接受垂直地带性控制［33］。吴吉春等［10］对祁连山区研究也表明，祁连山区多年冻土空间分布具有明显的垂直地带性。因此确定祁连山区地区的多年冻土分布面积，核心是查明祁连山区多年冻土分布下界。

Zhao等［34］提出多年冻土下界模型在局部大尺度的多年冻土分布研究中，多年冻土分布面积可能存在高估的现象，其中一部分原因可能与下界模型所采用的下界数据有关。传统的下界模型将发现多年冻土存在的最低海拔高度定义为该区域多年冻土分布下界，即将不连续多年冻土下界视为多年冻土下界［35］。而在多年冻土分布的边界区域的一定海拔范围内，真正存在多年冻土的区域十分有限。因此采用不连续多年冻土分布下界为依据的下界模型，在较大尺度的多年冻土分布研究中，一定程度上会高估多年冻土的分布面积，特别是在高原较为平坦区域，大量的非多年冻土被误判为多年冻土。而当达到大片连续多年冻土分布区的海拔高度时，多年冻土分布受局部因素影响减弱，只在少数区域或局部存在融区，其余区域一般均发育多年冻土，因此采用大片连续多年冻土分布区下界确定多年冻土分布范围显然更为合理。

本文按照不连续多年冻土分布区和大片连续多年冻土分布区分类建立多年冻土下界统计模型，将大片连续多年冻土分布区视为祁连山区多年冻土分布区、不连续多年冻土区视为多年冻土与季节冻土之间存在的过渡区，对祁连山区多年冻土的空间分布规律进行分析。本文所采用的多年冻土下界数据主要来源于三个方面：①青藏高原第二次综合科学考察实地踏勘获取的数据。②早年科学考察中基于道路勘察、冻土监测的祁连山区多年冻土区数据。③部分前人文献中所提及的祁连山区多年冻土下界。

2.1 多年冻土下界数据

在第二次青藏高原综合科考实地考察中，在祁连山区不同流域的溯源踏勘中依据冻土环境，将多年冻土下界按照不连续多年冻土分布区和大片连续多年冻土分布区分类统计（图2）。

由于存在某些对多年冻土的发育影响较大的局部因素，在一定的海拔高度范围内，多年冻土存在于某些最利于多年冻土发育的区域，如阴坡山沟出口的缓坡等地带，但由于不同区域的地形地貌等条件差异较大，多年冻土呈零星或岛状分布，出现这类多年冻土的区域最低海拔值，即可代表不连续多年冻土区下界；随着海拔升高，多年冻土在某些不利于多年冻土发育的区域也开始发育，如平坦的谷地等，局部因素对多年冻土发育影响减弱，多年冻土呈大片连续乃至连续分布，出现此类现象的最低海拔即可代表大片连续多年冻土区分布下界。由于不连续多年冻土分布区，存在大量的非多年冻土，多年冻土表现为零星或岛状分布，因此可以将该区域视为大片连续冻土分布区和季节冻土分布区之间的过渡区。

另一部分数据来源于多年冻土分布边界区域的钻孔点，主要是在疏勒河、德令哈、大通河等地或道路两侧人工钻取的钻孔点。钻孔点数据包括钻孔点的经纬度、冻土类型、上限深度、体积含冰量等信息，通过对设立在多年冻土分布边缘地带的钻孔点及周边环境信息分析，也可以获得相应地区不连续分布及大片连续分布的多年冻土区下界值。

图2 祁连山区科考踏勘路线Fig.2 Map showing the scientific investigation route in the Qilian Mountains

此外，祁连山区的多年冻土一直受到相关研究者的关注，不同的研究者开展了有关多年冻土不同方面的研究和调查工作，其中有些工作涉及到了祁连山区多年冻土下界。程国栋等［33］、吴吉春等［10］、王绍令等［15］、王庆锋等［36］在文献中提到多年冻土下界，这些工作大部分完成时间较早，虽然多年冻土退化的速度在加快［37］，活动层厚度以0.15 m·（10a）-1的速度增加［38］，但这对近几十年的尺度上的下界变化影响不是特别显著，所以这些早期调查的下界对目前多年冻土分布仍然具有指导意义。通过结合实地考察结果，经过甄别和具体分析后，也可以获得相应地区不连续分布及大片连续分布的多年冻土下界值。

将上述来源所获得的下界资料，分为大片连续分布与不连续多年冻土分布区下界数据进行整理，由于不连续多年冻土分布区多年冻土分布受局地地形影响较大，在不连续多年冻土分布区下界数据添加坡向分类，坡向基于数据点的经纬度坐标，在已有的高分辨率数字高程模型（DEM）中通过ESRI公司出版的ArcMap10.3软件进行提取，结果如下所示（表1和表2）。

表1 大片连续多年冻土分布区下界数据表Table 1 Longitude，latitude and the lower bound of large tracts of continuous permafrost in the Qilian Mountains

表2 不连续多年冻土分布区下界数据表Table 1 Longitude，latitude and the lower bound of discontinuous permafrost in the Qilian Mountains

2.2 多年冻土下界统计模型

本文利用线性回归方法，建立多年冻土下界统计模型，基于经纬度分析祁连山不同区域的多年冻土下界，探究祁连山区多年冻土的空间分布规律。由于大片连续多年冻土分布区海拔较高，多年冻土呈大片乃至连续分布，局部不利因素对多年冻土发育影响较弱。因此在建立回归方程时，只将受局部不利因素影响较大的不连续多年冻土分布区，按照阳坡阴坡分类建立回归方程，对于大片连续多年冻土分布区，只考虑经纬度变化对其下界变化的影响。通过对经纬度与下界的回归分析，得出：

大片连续多年冻土分布区下界统计模型：

不连续多年冻土分布区下界统计模型（阴坡）：

不连续多年冻土分布区下界统计模型（阳坡）：

式中：Long代表经度；Lat代表纬度。

F检验（F-test）通常是用来分析用了超过一个参数的统计模型，以判断该模型中的全部或一部分参数是否适合用来估计母体，可以较好的体现检验方程整体的显著性。对式（1）、式（2）与式（3）做回归方程显著性检验，由结果（表3）可知P＞F概率远小于0.05，结果表现显著，回归方程显著性较好，统计模型具有一定的可信度。

表3 回归方程显著性检验Table 3 Results of the significance test of regression equation

2.3 多年冻土分布模拟制图

在ESRI公司的ArcMap10.3版本中，将重采样的DEM数据转换成带经纬度及海拔数据（H）的点，并赋以基于DEM数据提取的坡向数据。依据每个点所在的经纬度坐标，通过式（1）计算出大片连续多年冻土分布区下界值H1，再依据坡向数据，分阴坡阳坡通过式（2）和式（3）计算出阴坡不连续多年冻土分布区下界值H2与阳坡不连续多年冻土分布区下界值H3。再将其海拔（H）与H1、H2、H3值对比，当若H＞H1，则认为其所在海拔高于大片连续多年冻土分布区下界值，即存在大片连续分布多年冻土，将其视为多年冻土区；若H≤H2或H≤H3，则认为不存在多年冻土，将其视为季节冻土区；若H≤H1，但H＞H2或H＞H3，即认为其所在海拔低于大片连续多年冻土分布区下界值，但高于不连续多年冻土分布区下界值，即存在不连续多年冻土，将其视为多年冻土区与季节冻土区之间的过渡区，由此获得祁连山区多年冻土空间分布模拟（图3）。此外，由于过渡区存在着大量的非多年冻土区，比较合理的估算祁连山区多年冻土分布面积的方法是只考虑大片连续多年冻土分布区所代表的多年冻土区的面积，将过渡区中存在的部分多年冻土视作大片连续多年冻土区中非多年冻土区部分的补偿，最终得出祁连山区多年冻土区面积为8.03×104km2，从多年冻土区向季节冻土区过渡的区域面积约1.43×104km2。

图3 祁连山区多年冻土空间分布模拟Fig.3 Spatial distribution of the simulated permafrost in the Qilian Mountains

3 验证与讨论

3.1 与实地钻孔点的交叉验证

本文所采用的钻孔点主要完成时间在2004—2011年，分布于道路两侧或部分河流流域。每个钻孔数据都包括钻孔点的经纬度、冻土类型、上限深度、体积含冰量等信息。钻孔数据共有548个，其中季节冻土的钻孔数为117个，多年冻土的钻孔数为431个，具体分布情况如下所示（表4、图4）。

钻孔点数据能很好的体现钻孔点附近多年冻土分布的真实情况，将下界统计模型计算得到的模拟结果与已有的钻孔数据做交叉验证（表5），结果表明：共有374个钻孔在模拟结果为多年冻土区内，其中有337个钻孔为多年冻土钻孔，多年冻土钻孔占比90.1%；共有169个钻孔点位于模拟结果为过渡区中，其中有92个钻孔为多年冻土钻孔，多年冻土钻孔占比54.4%。钻孔在钻取时存在一定的人为的干扰因素，这种干扰因素表现为在钻取钻孔时人为的倾向选择多年冻土存在概率较高的区域。由于多年冻土分布边界区域存在着大量的非多年冻土，这种人为干扰因素造成的影响，在过渡区对验证结果的影响尤为显著。因此可以推断，基于钻孔数据的验证的结果表明：模拟结果为多年冻土的区域，实际多年冻土面积占比应略低于90.1%，验证结果在一定程度上证明了所得到的多年冻土分布区具有较好的准确性；而在模拟结果为过渡区的区域，实际多年冻土面积占比应远低于54.4%，验证结果也表明，该区域存在着大量的季节多年冻土区。另一方面，由于处在较低海拔的季节冻土区钻孔数量较少，基于实地钻孔数据交叉验证的结果对于季节冻土区的验证不具有代表性。

表4 祁连山区已有钻孔点信息表Table 4 Information of the drilling points in the Qilian Mountains

图4 祁连山区钻孔点位置示意图Fig.4 Map showing the distribution of drill points in the Qilian Mountains

表5 多年冻土模拟与钻孔点交叉验证结果Table 5 Results of cross-validation between permafrost simulation result and borehole observation

3.2 与黑河上游冻土图的交叉验证

地表温度对冻土的空间分布有着重要的影响。葛社民［39］基于黑河上游源区内外九个站点0 cm处的地表温度，通过空间插值，结合冻土模拟方法，获得了黑河上游山区冻土模拟分布图（下称：黑河上游冻土图），制图成果将冻土分为多年冻土与季节冻土两类，其成果范围大致在本文所模拟的范围内（图5）。将其在本文模拟范围内的部分截取，用于与本文模拟的结果做叠加分析，以进一步验证模拟结果的准确性。

图5 黑河上游冻土图位置示意图Fig.5 Map showing the permafrost in the upper reaches of the Heihe River

将截取所获得的部分与本文模拟结果做叠加分析（图6），按照本文所模拟的结果，将叠加分析的结果划分为多年冻土区、过渡区和季节冻土区三个范围。

图6 与黑河上游冻土图的叠加分析结果Fig.6 Superposition analysis results with the permafrost map in the upper reaches of the Heihe River

叠加分析（图6）结果表明：模拟结果截取的部分与黑河上游冻土图有较好的一致性，差异主要集中在多年冻土发育的边界区域。在本文模拟结果为多年冻土区的范围内，黑河上游冻土图中整体上为多年冻土，约占该范围面积的96.46%；在本文模拟结果为季节冻土区的范围内，黑河上游冻土图中未出现多年冻土。由于所截取的部分黑河上游冻土图的制图使用的是九个站点的实测地表温度数据，具有一定的准确性，且制图方法与本文不同，因此尽管该冻土图只占本文所模拟区域的一小部分，但检验结果仍旧具有一定的代表性。叠加分析结果在一定程度上验证了本文模拟结果中多年冻土区和季节冻土区的准确性。过渡区中存在着大量的季节冻土，叠加分析结果也有较好的体现。在过渡区范围内，黑河上游冻土图中多年冻土主要集中在本文模拟结果为多年冻土区的边界区域，多年冻土面积约占过渡区总面积的24.65%，验证结果对本文模拟结果有较好的支撑。

3.3 与其他制图成果的比较

本文模拟结果得出祁连山区多年冻土区面积为8.03×104km2，占研究区总面积的47.5%，从多年冻土区向季节冻土区过渡的区域面积约1.43×104km2，占研究区总面积的8.5%。研究区内多年冻土在空间分布上表现为以哈拉海为中心向四周扩散分布的格局，哈拉湖至祁连山区外围山区多年冻土呈大片连续分布，外围山区多年冻土受山体、河谷的地形控制逐渐呈不连续分布，直至过渡到季节冻土区。将模拟结果与《青藏高原冻土图》［13］、《中国冰川冻土沙漠图》［12］以及青藏高原新绘制冻土分布图（2017）［40］祁连山区部分比较，结果如图7所示。

受海拔因素影响，祁连山区多年冻土发育直接受垂直地带性规律控制。由上图可以直观得出，本文基于祁连山区实地踏勘数据，利用冻土下界由高程之间的统计关系，建立的空间模型所模拟的结果整体上与《青藏高原冻土图》《中国冰川冻土沙漠图》以及青藏高原新绘制冻土分布图（2017）中祁连山部分相似，差异主要集中在多年冻土边界区域。

《青藏高原冻土图》同样是基于冻土下界由高程的统计关系所模拟的结果。研究区在《青藏高原冻土图》中多年冻土面积为9.02×104km2，该值接近与本文模拟结果中多年冻土区面积与过渡区面积之和9.46×104km2。造成与本模拟结果边界区域差异的主要原因可能与前者所依据的下界数据较低有关。即《青藏高原冻土图》在多年冻土分布边界区域，将大量季节冻土判断为多年冻土，在一定程度上，高估了祁连山区多年冻土的分布面积。

《中国冰川冻土沙漠图》主要基于青藏铁路沿线钻孔获取的平均地温数据，建立地温经验公式，从而对青藏高原整体多年冻土分布进行模拟。研究区在《中国冰川冻土沙漠图》中多年冻土分布整体上与本文模拟结果差异较小，差异较大的区域集中在祁连山区东部多年冻土分布边界区域，特别是大通河以及布哈河流域。在流域区域，地温受地表植被、土壤含水率影响较大，因此较以地温经验公式为基础的《中国冰川冻土沙漠图》，本文的模拟结果具有更高的可靠性。

图7 《青藏高原冻土图》（a）、《中国冰川冻土沙漠图》（b）、青藏高原新绘制冻土分布图（2017）（c）中的祁连山区多年冻土分布及本文模拟结果（d）Fig.7 Permafrost distribution over the Qilian Mountains，quoted from：the map of permafrost distribution over the Tibetan Plateau（1996）（a），the map of glacier，permafrost and deserts in China（2005）（b），a new map of permafrost distribution on the Tibetan Plateau（2017）（c）and the result of this paper（d）

青藏高原新绘制冻土分布图（2017）是利用中分辨率成像光谱仪（MODIS）地表温度（LSTs）的冻融指数及多年冻土顶部温度（TTOP）模型模拟生成的。由于成图时间较晚，分辨率较高，且综合了多种遥感数据，因此较《青藏高原冻土图》与《中国冰川冻土沙漠图》具有更高的可靠性。除去资料缺失区域，本文模拟结果在青藏高原新绘制冻土分布图（2017）中差异较小，多年冻土分布结果基本一致。但在大通河流域、塔塔棱河流域以及党河流域等区域，对于多年冻土分布边界区域模拟有一定差异，这可能与该区土壤含水率较高，植被发育较好，从而对地表温度遥感解译产生了影响有关。此外，由于本文模拟结果基于地带性规律，对于祁连山区多年冻土模拟没有出现资料缺失区域，因此本文模拟结果可以更好的反映祁连山区多年冻土分布的情况。

综上，本文模拟结果与早年出版的青藏高原多年冻土分布图中祁连山部分，具有相似的空间分布特征，与近年新绘制的青藏高原多年冻土分布图中祁连山部分具有较好的一致性。由于本文依据的多年冻土下界点是依据实地勘探所获取的，在祁连山区域综合了更多的实地资料，具有更高的可靠性。此外，较已出版的三幅图件，还给出了位于多年冻土与季节冻土之间的过渡区的范围，因此在祁连山地区本文模拟结果具有更好的适用性。

4 结论与展望

多年冻土分布面积一般来自模型估算，而估算的准确性取决于模型的合理性和关键数据量的支持，祁连山区多年冻土分布主要控制因素是海拔，多年冻土下界是确定其多年冻土分布的关键参数。本文基于青藏高原第二次综合科考，结合冻土现象分析和核实冻土资料缺失区域多年冻土下界，获取了较为详实的关键数据量，并建立了多年冻土分布的统计模型，绘制了祁连山区多年冻土空间分布模拟图，并通过已有的钻孔点资料和前人所做的多年冻土图进行了验证与比较。根据下界统计模型及制作完成的祁连山区多年冻土空间分布图，对祁连山区多年冻土空间分布特征的基本认识归纳如下：

（1）祁连山区多年冻土在空间上呈现出以哈拉湖为中心向四周扩散的分布格局。哈拉湖至祁连山区外围山区多年冻土大片连续分布，外围山区多年冻土受山体、河谷的地形控制逐渐呈不连续分布，直至过渡到季节冻土区。

（2）在多年冻土区与季节冻土区之间存在着过渡区，由于过渡区存在着大量的非多年冻土区，比较合理的估算祁连山区多年冻土分布面积的方法是只考虑大片连续多年冻土分布区所代表的多年冻土区的面积，将不连续多年冻土中存在的部分多年冻土视作大片连续多年冻土区中非多年冻土区部分的补偿，根据统计模型分析，最终得出祁连山区多年冻土区面积为8.03×104km2，从多年冻土区向季节冻土区过渡的区域面积约1.43×104km2。

（3）祁连山区多年冻土分布具有良好的地带性规律，多年冻土下界随经纬度增加而降低。在祁连山区东北部多年冻土下界最低，如冷龙岭一带大片连续分布多年冻土下界约为3 540 m，不连续多年冻土分布的下界约3 410 m；在西南部多年冻土下界最高，如大柴旦附近，大片连续多年冻土分布下界达到了4 050 m，不连续多年冻土分布下界也达到了3 900 m。整体上下界值呈西南＞西北＞东南＞东北的规律。

本文以实地考察资料为基础，估算了祁连山区多年冻土空间分布，所获得的下界统计模型模拟结果较好，为祁连山区多年冻土空间分布提供了一定的参考。另一方面，受条件限制，祁连山区西部以及季节冻土区的钻孔点较少，尚需进一步的工作以对目前已取得的模拟结果进行进一步验证和改良。

致谢：《黑河上游山区冻土模拟分布图》数据来源于“国家青藏高原科学数据中心”(http：//data.tpdc.ac.cn)，《青藏高原冻土图》与《中国冰川冻土沙漠图》数据来源于“中国西部环境与生态科学数据中心”(http：//westdc.westgis.ac.cn/)。