崔海濤， 鄭超剛， 張志剛，2，3， 張 宏， 郭 飛， 李奕曼， 王 建

（1. 南京師范大學 地理科學學院，江蘇 南京 210023； 2. 南京師范大學 虛擬地理環(huán)境教育部重點實驗室，江蘇 南京 210023；3. 江蘇省地理信息資源開發(fā)與利用協(xié)同創(chuàng)新中心，江蘇 南京 210023）

0 引言

青藏高原作為世界上海拔最高，我國面積最大的高原，是除南北極以外發(fā)育現(xiàn)代冰川最大的區(qū)域，也是我國保留第四紀冰川遺跡最豐富的區(qū)域［1］。由于所處的地理位置比較特殊，在亞洲甚至全球氣候變化中有著重要影響［2-4］。第四紀冰川遺跡詳細記錄了古冰川變化情況，是了解第四紀環(huán)境演化的重要載體［1］，對古冰川遺跡的空間分布及定量化研究，有助于認識冰期規(guī)模大小、恢復第四紀冰川演化史、了解古環(huán)境變化等［5-6］。

第四紀冰川地貌制圖作為記錄古冰川遺跡空間分布特征的重要載體［7］，其發(fā)展可分為定性研究和定量研究兩個階段。20 世紀50 年代到20 世紀末，是第四紀冰川地貌制圖定性研究階段，隨著前蘇聯(lián)學者馬爾科夫總結了地貌制圖的方法、形式，一些學者在青藏高原的不同區(qū)域制成了不同比例尺大小的冰川地貌圖［8-11］。21 世紀以來，隨著遙感監(jiān)測和無人機航拍在地貌學上的成熟運用，第四紀冰川地貌制圖從定性研究向定量研究轉(zhuǎn)變［12-16］。當前，已有學者對青藏高原地區(qū)的一些區(qū)域（青藏高原東北部巴顏喀拉山、青藏高原中部、青藏高原東南部沙魯里山、青藏高原東北部達里加山、念青唐古拉山中部）進行數(shù)字化冰川地貌制圖［17-22］，但是多數(shù)研究屬于半定量化，同時缺乏單個冰川地貌類型進行專題制圖。

位于青藏高原東南緣理塘縣和稻城縣之間的稻城古冰帽，第四紀時曾數(shù)次發(fā)育冰帽并保存了大量的古冰川遺跡，是開展綜合冰川地貌研究的理想?yún)^(qū)域。20世紀60年代起，一些學者陸續(xù)對該區(qū)的第四紀冰川遺跡分布進行研究［23］，姚檀棟等［24］根據(jù)巖盆、槽谷地貌分布特征，恢復了古冰帽范圍；李吉均等［25］簡要畫出了冰磧物、槽谷的地理位置以及三次冰期的作用范圍；鄭本興等［26］制作了冰川槽谷、冰蝕盆地分布圖，劃分了古冰帽侵蝕、堆積帶；Fu 等［19，27］在稻城古冰帽識別了冰磧壟、冰川谷等五種冰川地貌類型，描述了冰川地貌分布特征并建立了詳細的古冰川地貌圖。此外，部分學者基于冰川地貌年代學重建了該區(qū)第四紀冰川序列［28-36］。

目前，稻城古冰帽區(qū)第四紀冰川地貌數(shù)字化研究有待于進一步細化，冰川地貌類型專題圖有待補充完善。本文采用目視解譯和野外考察相結合的方法對該區(qū)冰川地貌進行識別和定量化研究，從而為該區(qū)古冰川遺跡在空間分布及數(shù)量獲得直觀上的認識，為第四紀冰期規(guī)模估算、古氣候重建提供基礎數(shù)據(jù)，也為當?shù)氐穆糜魏唾Y源調(diào)查提供一定的參考。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況

稻城古冰帽位于青藏高原東南緣（29°02′～30°08′ N，99°48′～100°50′ E）（圖1），坐落于橫斷山脈中沙魯里山中部高原區(qū)。介于德格—鄉(xiāng)城斷裂帶和甘孜—理塘斷裂帶之間，出露巖性以印支期花崗巖為主。研究區(qū)北部與毛埡壩盆地和理塘盆地相鄰，東部與無量河斷裂帶和甲洼—康嘎盆地相接，南部以稻城河斷裂帶為界，西部與稻城河谷和希曲毗鄰［26，29-30，37］。稻城古冰帽主體部分是由甘孜州南部稻城縣和理塘縣之間的海子山夷平面及其殘丘組成［30］，長約135 km，呈西北—東南走向（NWSE），面積達3 600 km2，平均海拔在4 500 m左右。

圖1 稻城古冰帽位置示意圖Fig. 1 The location of paleo-Daocheng ice cap

目前整個稻城古冰帽區(qū)已無現(xiàn)代冰川，第四紀時該區(qū)曾多次發(fā)育冰帽并保留了包括冰川槽谷、冰蝕盆地、古冰斗、冰磧壟、羊背巖等在內(nèi)的各種冰川侵蝕堆積地貌［24-25］。根據(jù)冰川地貌空間分布特征，稻城古冰帽被劃分為冰帽中心微弱侵蝕和堆積帶、中等侵蝕和堆積帶、強烈侵蝕和冰蝕槽谷帶以及古冰帽邊緣強烈堆積帶［26］（圖2），記錄了冰川作用的整個過程。

圖2 稻城古冰帽冰川地貌特征圖（據(jù)鄭本興等［26］改繪）Fig. 2 The map of glacial geomorphic zones in the palco-Daocheng ice cap region （Modified after Zheng Benxing et al.［26］）

1.2 研究方法

本文遙感影像數(shù)據(jù)來源于Google Earth Pro 在線地圖，數(shù)字高程模型（DEM）數(shù)據(jù)來源于SRTMDEM 12.5 m，下載網(wǎng)址：ftp：//ftp.glcf.umiacs.umd.edu/glcf/SRTM/Degree Tiles/。

采用目視解譯和野外考察相結合的方法識別稻城古冰帽第四紀冰川地貌類型，根據(jù)解譯人的專業(yè)知識和經(jīng)驗，對遙感影像進行判讀與分析，通過影像上物體的形狀、位置、色調(diào)等特征，識別影像上目標物體及獲取分布特征信息并對部分冰川地貌進行現(xiàn)場驗證。首先，在稻城古冰帽區(qū)多次野外考察的基礎上，并結合多位學者在此研究古冰川地貌的成果，初步認識該研究區(qū)第四紀冰川侵蝕堆積地貌類型，深刻了解該區(qū)冰川湖、冰川谷、羊背巖、冰磧壟四種冰川地貌成因、形態(tài)特征、不同地貌之間的聯(lián)系及可能形成的區(qū)域位置（表1），這為接下來目視解譯奠定了基礎。其次，根據(jù)該區(qū)冰川地貌特征（圖3）并參考其他學者識別冰川地貌的方法及原理［17-19］，對Google Earth Pro 在線影像圖進行仔細判讀，獲取該區(qū)第四紀冰川地貌類型及空間分布。再次，通過野外考察對目視解譯結果進一步驗證。最后，對識別的冰川地貌進行分類，導出KML 格式文件（. kml）并依次加載到ArcGIS 10.5 中。另外，依據(jù)冰川地貌類型所屬規(guī)范色彩進行上色處理，區(qū)分不同冰川地貌類型，并選取WGS 1984 UTM Zone 47N 坐標進行投影處理，在屬性表中計算不同冰川地貌的面積、數(shù)量等指標數(shù)據(jù)，基于該區(qū)數(shù)字高程模型制成對應的冰川地貌專題圖。

表1 冰川地貌特征Table 1 Glacial landform features

圖3 稻城古冰帽區(qū)冰川地貌影像Fig. 3 The map of glacial geomorphic imagery in the palco-Daocheng ice cap area： glacial lake （a）， glacial valley （b），roche moutonnee （c）， glacial moraine （d）

2 結果與分析

2.1 冰川湖

冰川湖常見有冰蝕湖、冰磧湖，因稻城古冰帽范圍較大，且歷經(jīng)多次冰川作用，湖泊成因難以嚴格區(qū)分，故本文統(tǒng)稱冰川湖。在稻城古冰帽區(qū)，共識別1 096 個冰川湖，面積約60 km2，占古冰帽1.7%，平均密度達0.3個.km-2。主要分布在海子山夷平面上（圖4）。冰川湖按面積大小分為6個等級，其中面積小于0.01 km2的冰川湖占冰川湖總量的46.08%，面積在0.02 km2以上的冰川湖占冰川湖總量的34.21%，0.01～0.02 km2之間的冰川湖占總量的19.71%（表2）。

表2 冰川湖面積、數(shù)量統(tǒng)計Table 2 Statistics on the area and quantity of glacial lakes

圖4 稻城古冰帽區(qū)冰川湖泊分布圖與實景圖Fig. 4 Distribution map and scene map of the glacial lakes of the paleo-Daocheng ice cap area： distribution map of the glacial lakes of the paleo-Daocheng ice cap area （a），photo of Xingyicuo lake （b）

其中，古冰帽區(qū)最大的冰川湖——興伊措，由夷平面上的兩列冰磧壟堵住出口而形成，面積達6.1 km2，湖深3～4 m，最深處可達10 m，是稻城河的主要源頭，分布在海子山中部，受氣候變暖影響，冰川湖面積縮?。▓D4）。

2.2 冰川谷

稻城古冰帽區(qū)的古冰斗因后期長時間受到外動力作用影響，而且該區(qū)曾多次冰帽發(fā)育，致使形成的冰斗形態(tài)特征不明顯，因此，我們把冰川槽谷（“U”形谷）和冰斗統(tǒng)稱為冰川谷進行地貌識別。

關于研究區(qū)的冰川谷地貌，共識別約370條，約1 166 km2，占古冰帽區(qū)32.4%。分布上，主要存在于古冰帽區(qū)邊緣地帶，其中，南部和北部的山地冰川谷發(fā)育顯著，而中部地區(qū)發(fā)育較少（圖5）。

圖5 稻城古冰帽區(qū)冰川谷地分布圖、實景圖及冰川谷剖面圖Fig. 5 Distribution map， scene map and glacial valley profile of the glacial valley of the paleo-Daocheng ice cap area：distribution map of the glacial valley of the paleo-Daocheng ice cap area （a）， glacial valley of Kuzhaori （b）， glacial valley of Xionggu （c）， cross section profiles of glacial valleys （x，y） and fluvial valley （z） （d）

2.3 羊背巖

羊背巖，是經(jīng)歷冰川侵蝕作用后形成于地面上的石質(zhì)小丘。形態(tài)特征表現(xiàn)為：頂部渾圓，形似羊背，具有卵形的基部，羊背巖具有迎冰坡緩、背冰坡陡的特點，這是由于迎冰坡以磨蝕作用為主，導致坡面比較平緩光滑；背冰坡以挖蝕作用為主，導致坡面比較陡峻粗糙（破碎）。對海子山夷平面考察時，發(fā)現(xiàn)數(shù)個羊背巖，它們占地面積較大，海拔較高。而且羊背巖長軸延伸方向（從坡緩處向坡陡處）與冰川運動方向一致，故羊背巖可作為指示冰川運動方向的標志物（圖6）。

圖6 稻城古冰帽區(qū)羊背巖分布圖與實景圖Fig. 6 Distribution map and scene map of the roche moutonnees of the paleo-Daocheng ice cap area： distribution map of the roche moutonnees of the paleo-Daocheng ice cap area （a），scene map of the roche moutonnee 1 （b）， scene map of the roche moutonnee 2 （c）

依據(jù)野外考察經(jīng)驗和羊背巖形態(tài)特征，在研究區(qū)共識別41 個羊背巖，面積約14 km2，占古冰帽面積0.4%。其中，最大羊背巖面積達到1.07 km2，但面積在0.2～0.4 km2的數(shù)量最多，占比高達51.22%，面積達0.6 km2以上的數(shù)量較少（表3）。從分布上看，羊背巖主要形成于古冰帽中心地區(qū)，主體部分集中在海子山夷平面向南部山地的過渡區(qū)域，多數(shù)分布在S217公路北側。

表3 羊背巖面積、數(shù)量統(tǒng)計Table 3 Statistics on the area and quantity of roche moutonnees

2.4 冰磧壟

冰磧壟是在冰川消融過程中，由冰川搬運的物質(zhì)堆積而成的壟狀地形，是典型的第四紀冰川遺跡。本研究區(qū)保留了大量冰磧壟地貌，以終磧壟和側磧壟為主，終磧壟是指在冰川末端由冰磧堆積形成的弧形壟狀地形，側磧壟是指在冰川兩側邊緣形成的壟狀地形。

曾有多位學者在稻城古冰帽不同區(qū)域的冰磧壟上采集漂礫、冰磧碎屑等暴露測年樣品（圖7），不僅為冰磧壟的識別提供了地理位置信息，而且進一步驗證了冰磧壟識別的精確性。此外，依據(jù)在Google Earth Pro 數(shù)據(jù)中預設采樣點位置與實際野外采集暴露測年樣品時用GPS 定點位置信息相吻合，表明制圖精度較高。在該區(qū)，共識別約1 268 列冰磧壟，但受后期冰川、流水、風化等侵蝕作用影響，一些冰磧壟遭受破壞，形態(tài)特征不再顯著甚至消亡。地形上，冰磧壟主要保留在海子山夷平面、山麓地帶和槽谷末端；空間上，呈現(xiàn)南多北少，其中夷平面上有3 組清晰可見的冰磧壟發(fā)育，表面分布大量排列有序的冰川漂礫。其余則分布在南北兩端的冰川谷邊緣地帶，以庫照日冰磧壟序列最多，形態(tài)特征最為典型，是重建青藏高原第四紀冰期序列的重要區(qū)域（圖8）。

圖7 稻城古冰帽區(qū)冰磧壟上采樣點及對應冰期分布圖Fig. 7 Distribution map of sampling sites and corresponding glacial stage on moraine of the paleo-Daocheng ice cap area

圖8 稻城古冰帽區(qū)冰磧壟分布圖與實景圖Fig. 8 Distribution map and scene map of the glacial moraine of the paleo-Daocheng ice cap area： distribution map of the glacial moraine of the paleo-Daocheng ice cap area （a），sampling site 1 on moraine in July 2021 （b）， sampling site 2 on moraine in July 2021 （c）

2.5 稻城古冰帽冰川地貌圖

在稻城古冰帽區(qū)已識別冰川湖、冰川谷、羊背巖、冰磧壟四種第四紀冰川地貌類型，冰川湖主要分布于冰川古道內(nèi)和夷平面上；羊背巖形成于冰川侵蝕區(qū)的海子山夷平面中心區(qū)域；冰磧壟環(huán)繞冰川谷尾端發(fā)育。根據(jù)第四紀冰川地貌圖可知，不同的冰川地貌類型，其形成的地域有所差異；在位置關系上，彼此之間存在一定的聯(lián)系（圖9）。

圖9 稻城古冰帽區(qū)冰川地貌分布圖Fig. 9 Distribution map of the glacial landform of the paleo-Daocheng ice cap area

3 討論

從野外考察到影像判讀，共識別1 096 個冰川湖、370 條冰川谷、41 個羊背巖和1 268 列冰磧壟。稻城古冰帽區(qū)作為研究青藏高原第四紀冰川地貌最早區(qū)域之一，姚檀棟等［24］統(tǒng)計稻城古冰帽區(qū)域巖盆1 145 個，平均密度0.3 個·km-2?；诒ê梢虻囊巴饪疾煅芯堪l(fā)現(xiàn)，海子山夷平面上最大冰川湖——興伊措，主要是冰川前進過程中侵蝕形成的巖盆，后期冰川后退過程中形成的冰磧壟或冰磧物堵住了巖盆出口，積水形成的冰川湖泊［26］。分布上契合了古冰帽冰川湖以海子山夷平面和冰川谷為主特點，成因上有待進一步研究，需加強對野外深入考察。

冰川谷主要分布于稻城古冰帽邊緣地帶，可能與海子山夷平面地勢起伏小有關，第四紀冰川發(fā)育以冰帽形式覆蓋整個夷平面，致使難以形成冰川谷。古冰帽邊緣區(qū)，因構造斷裂和流水侵蝕作用而形成地勢起伏較大的地形，冰川重力作用加上地形陡峭，古冰川運動速度較快，侵蝕原有的河谷河床而形成冰川谷，其中以南北兩端的山谷分布最為顯著。這與其他學者描述冰流槽谷分布位置一致［24-26］。凸顯了稻城古冰帽是以海子山夷平面為中心的冰帽區(qū)和以邊緣區(qū)為主的山谷冰川構成。

關于羊背巖，最早由李吉均等［25］在橫斷山二次考察時發(fā)現(xiàn)，羊背巖主要形成于稻城古冰帽中心區(qū)域，其特征是迎冰面比較平緩、光滑，背冰面陡峭而破碎，部分羊背巖表面有冰磧堆積物。鄭本興等［26］同樣在稻城古冰帽中心區(qū)域發(fā)現(xiàn)了羊背巖存在。在對稻城古冰帽區(qū)野外考察發(fā)現(xiàn)，羊背巖發(fā)育于海子山夷平面中心區(qū)，該區(qū)域基巖裸露而且冰川以侵蝕作用為主［26］，是塑造羊背巖的理想?yún)^(qū)。其長軸指向指示了古冰川的前進方向與實際中古冰川前進方向相同（NE-SW），間接表明了該工作的可靠性。

稻城古冰帽冰磧物主要以終磧壟和側磧壟形式保存，是研究第四紀冰期的重要標志。李吉均等［25］在海子山夷平面上發(fā)現(xiàn)了多列冰磧壟，并認為它們是同一冰期作用產(chǎn)物，此外在庫照日地區(qū)發(fā)現(xiàn)了具有歷史意義的紅色風化殼覆蓋的終磧壟。鄭本興等［26］研究表明，冰磧壟主要分布于古冰帽邊緣槽谷地帶，特別是稻城河北側的山麓地帶冰磧壟發(fā)育顯著。21 世紀以來，關于冰磧壟的研究主要集中于冰川年代學，分別在古冰帽西緣冰磧壟、海子山夷平面冰磧壟、海子山北端冰磧壟以及古冰帽東緣冰磧壟都有詳細描述與記載［29-36］。Fu 等［33］研究發(fā)現(xiàn)，海子山夷平面上清晰分布著3道主線冰磧壟，形成于末次冰消期。其他冰磧壟主要分布在山谷山麓地帶，主要是冰進作用而成，表明在第四紀冰期，海子山夷平面上發(fā)育形成冰帽。

綜上所述，通過與其他學者研究結果相比，驗證了本文冰川地貌數(shù)據(jù)化研究的可靠性、真實性和精確性。此外加強野外考察與影像判讀之間的聯(lián)系，分析不同冰川地貌之間存在的位置關系，從而提高冰川地貌數(shù)字化研究的準確性。

4 結論

以第四紀冰川遺跡豐富的稻城古冰帽為研究區(qū)，利用目視解譯和現(xiàn)場考察相結合的方法，在Google Earth Pro 遙感影像上識別了冰川湖、冰川谷、羊背巖和冰磧壟四種冰川地貌類型并制成對應的稻城古冰帽區(qū)冰川地貌分布圖。

初步統(tǒng)計，冰川湖約1 096 個，主要分布在海子山夷平面上；冰川谷約370條，以南北兩端的山谷為主；羊背巖約41個，發(fā)育于稻城古冰帽中心區(qū)域；冰磧壟約1 268 列，主要分布在東西兩側的冰川谷尾端、山麓地帶以及海子山夷平面上。

從地貌類型和分布上表明稻城古冰帽由以南部山地為主形成的山谷冰川和以海子山夷平面為主形成的冰帽構成，而且很好地指示了古冰川的前進方向和冰川性質(zhì)變化，為后面的古冰川地貌特征研究提供參考。