祁連山區多年凍土空間分布模擬

彭晨陽，盛 煜，吳吉春，曹 偉，何彬彬

（1.中國科學院西北生態環境資源研究院凍土工程國家重點實驗室，甘肅蘭州730000；2.中國科學院大學，北京100049）

0 引言

多年凍土是指持續兩年或兩年以上溫度處于0℃及以下巖土層［1］，高海拔多年凍土泛指高緯度多年凍土南界以南、一定海拔高度以上出現的多年凍土［2］。我國高海拔多年凍土面積居世界首位，總面積約173.20×104km2，其中約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×1010m3。近年來祁連山區中部降水明顯增加［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/)。