王凱，孫永玲，孫世金，劉曉

（山東理工大學(xué)建筑工程學(xué)院，山東 淄博 255022）

0 引言

山地冰川是全球氣候變化的敏感指示器［1-2］，20世紀(jì)90年代以來(lái)，全球氣溫升高導(dǎo)致冰川加速消融，冰川表面流速是研究冰川動(dòng)力學(xué)的重要參數(shù)，因此監(jiān)測(cè)冰川流速變化變得至關(guān)重要［3-4］。但山地冰川地處偏遠(yuǎn)、自然條件惡劣、人力測(cè)量成本昂貴等使實(shí)地監(jiān)測(cè)變得困難，遙感技術(shù)的發(fā)展為冰川監(jiān)測(cè)提供了低成本、大范圍、高效率等有利監(jiān)測(cè)條件［5］，特別是光學(xué)遙感監(jiān)測(cè)具有視野廣闊、獲取信息量多、適應(yīng)性強(qiáng)、可用于動(dòng)態(tài)監(jiān)測(cè)等突出優(yōu)點(diǎn)，已廣泛應(yīng)用于全球范圍內(nèi)的冰川流速提?。?］。

冰川躍動(dòng)是一種特殊的冰川運(yùn)動(dòng)現(xiàn)象，是冰川不穩(wěn)定性的一種表現(xiàn)形式［7］。冰川在躍動(dòng)期一般表現(xiàn)出冰川表面流速比平靜期流速增大數(shù)十倍甚至數(shù)百倍［8］、冰川厚度在積累區(qū)減薄，消融區(qū)增厚和冰川末端前進(jìn)等特征［9］。冰川躍動(dòng)可能阻塞河流而引起水災(zāi)或冰川滑入阻塞湖引發(fā)冰湖潰決洪水，嚴(yán)重威脅下游區(qū)域人們的生命財(cái)產(chǎn)安全［10-11］。因此，研究躍動(dòng)冰川運(yùn)動(dòng)特征對(duì)預(yù)警冰川躍動(dòng)導(dǎo)致的冰川災(zāi)害具有重要的意義［12-13］。Copland等［14］確定音蘇蓋提冰川為躍動(dòng)型冰川，并發(fā)現(xiàn)該冰川含有大量折疊狀的冰磧物、冰川末端前進(jìn)5 km，該冰川第一次躍 動(dòng)發(fā) 生 在1973年 或 者更 早［15］。Sun等［16］根據(jù)Landsat ETM+數(shù)據(jù)獲取音蘇蓋提冰川1996—2006年的年平均冰川流速，發(fā)現(xiàn)音蘇蓋提冰川在1996—2003年流速保持在20～80 m·a-1，2004—2006年冰川流速整體變大，流速最大約為100 m·a-1。蔣宗立等［17］利用ALOS PALSAR數(shù)據(jù)估計(jì)2007—2009年音蘇蓋提冰川流速，發(fā)現(xiàn)音蘇蓋提冰川主冰川年平均流速約為100 m·a-1，而南支流比較活躍，冰川流速較大（流速最大約為350 m·a-1）。Rankl等［18］利用ALOS PALSAR和TerraSAR-X數(shù)據(jù)提取喀喇昆侖區(qū)域流速，發(fā)現(xiàn)音蘇蓋提冰川在2003—2009年年平均流速相對(duì)較低（流速最大約為109 m·a-1），而2009—2011年年平均流速開(kāi)始變大，該冰川在2011年發(fā)生躍動(dòng)，最大冰川流速約為547 m·a-1。

本文以中國(guó)境內(nèi)面積最大的音蘇蓋提冰川為研究區(qū)，選取了2018年1月1日—2021年11月16日20對(duì)Landsat-8影像，利用光學(xué)影像特征追蹤方法提取音蘇蓋提冰川表面流速，評(píng)估流速不確定性，并分析音蘇蓋提冰川流速的時(shí)空變化特征，此外，本文還利用ITS_LIVE數(shù)據(jù)分析了音蘇蓋提冰川北支流（西）、南支流、北支流（東）的流速時(shí)間變化特征及最大流速與其高程分布。

1 研究區(qū)域和數(shù)據(jù)

1.1 研究區(qū)域

音蘇蓋提冰川（35°55′～36°15′N，75°55′～76°21′E），位于中國(guó)新疆喀喇昆侖山脈喬戈里峰的北坡（圖1），總長(zhǎng)約42 km，冰舌長(zhǎng)約4 200 m，面積達(dá)390 km2，最高海拔約7 050 m，冰川末端海拔約4 000 m，是中國(guó)面積、儲(chǔ)冰量最大的山谷冰川［19-20］。地形上呈NWW-SEE走向，該冰川是由4條主要支流和10余條小支流匯合而成的樹(shù)枝狀山谷冰川，該區(qū)域?qū)儆诘湫痛箨懶蜌夂颍瑲鉁剌^低，降水稀少［21］。

圖1 音蘇蓋提冰川位置（黑色實(shí)線為音蘇蓋提冰川邊界，白色實(shí)線為冰川中心線，白色虛線是冰川橫剖面線，底圖是2018年4月7日獲取的Landsat-8影像）Fig.1 Location of Skamri Glacier（The black solid line is the border of Skamri Glacier.The white solid line is the central line of Skamri Glacier.The white dotted line is the glacier cross section.The background is Landsat-8 image from April 7，2018）

1.2 數(shù)據(jù)

本文選取L1T（Level 1T）完全覆蓋音蘇蓋提冰川的Landsat-8 OLI全色波段影像數(shù)據(jù)，數(shù)據(jù)全部來(lái)源于美國(guó)地質(zhì)調(diào)查局（USGS）（https：//earthexplorer.usgs.gov/）。OLI傳感器的全色波段數(shù)據(jù)空間分辨率為15 m，衛(wèi)星每16天可以實(shí)現(xiàn)一次全球覆蓋［22］。本文選取云覆蓋率小10%的影像，以消除云覆蓋對(duì)冰川流速提取的影響，共計(jì)有20對(duì)Landsat-8影像用于提取2018年1月1日到2021年11月16日音蘇蓋提冰川表面流速（表1）。

表1 覆蓋音蘇蓋提冰川的Landsat-8數(shù)據(jù)Table 1 Landsat-8 images covering Skamri Glacier

SRTM（Shuttle Radar Topography Mission）是美國(guó)奮進(jìn)號(hào)航天飛機(jī)利用成像雷達(dá)于2000年2月11日至22日成功獲取60°N至56°S間80%的陸地面積，對(duì)原數(shù)據(jù)處理后，獲取平面精度±20 m，高程精度±16 m的全球數(shù)字高程模型（DEM）［23］。公開(kāi)發(fā)布SRTM數(shù)字高程產(chǎn)品由美國(guó)國(guó)家航空航天局（NASA）（https：//earthexplorer.usgs.gov/）提供。由SRTM獲取的DEM數(shù)據(jù)有30 m和90 m兩種分辨率的產(chǎn)品［24］。本文選取的是SRTM 30 m DEM數(shù)據(jù)，用以獲取音蘇蓋提冰川的表面高程。

2 方法

2.1 冰川邊界提取

冰川邊界是描述冰川位置和規(guī)模的重要特征，也是監(jiān)測(cè)冰川變化信息的重要參數(shù)。目前冰川邊界提取有目視解譯和計(jì)算機(jī)解譯兩種方法。本文結(jié)合了全球冰雪中心（http：//www.glims.org/）發(fā)布的RGI 6.0（Randolph Glacier Inventory）［25］冰川邊界數(shù)據(jù)為基礎(chǔ)，結(jié)合音蘇蓋提冰川區(qū)域2018年1月1日—2021年11月16日的Landsat-8影像數(shù)據(jù)，采用目視解譯方法對(duì)冰川邊界進(jìn)行調(diào)整，確定音蘇蓋提冰川邊界。

2.2 冰川流速提取

本文利用光學(xué)影像特征追蹤方法提取音蘇蓋提冰川表面流速，其方法的基本原理是對(duì)于不同時(shí)期冰川表面的遙感影像進(jìn)行相關(guān)性分析。由于冰川表面特征點(diǎn)是隨冰川一起運(yùn)動(dòng)的，計(jì)算兩景影像特征點(diǎn)移動(dòng)的位移，即可獲取冰川表面流速［26-27］。圖2為提取音蘇蓋提冰川年平均流速的流程圖。由于本文選取的L1T數(shù)據(jù)已經(jīng)進(jìn)行過(guò)輻射校正和地面控制點(diǎn)幾何校正，所以本文中的Landsat-8影像數(shù)據(jù)均不需要做正射校正，只需做影像配準(zhǔn)和互相關(guān)計(jì)算。本文利用ENVI軟件對(duì)獲取時(shí)間較早的Landsat-8數(shù)據(jù)為參考影像，并以獲取時(shí)間較晚的影像為匹配影像進(jìn)行配準(zhǔn)，并利用COSI-Corr（Co-registration of Optically Sensed Images and Correlation）軟件采用頻率域互相關(guān)算法對(duì)配準(zhǔn)后的影像對(duì)進(jìn)行互相關(guān)計(jì)算，獲取音蘇蓋提冰川東西向和南北向偏移量和信噪比［28-29］。在互相關(guān)計(jì)算中，首先以初始窗口進(jìn)行互相關(guān)計(jì)算，得到像素級(jí)的冰川偏移量，再以最終窗口進(jìn)行互相關(guān)計(jì)算得到冰川亞像素級(jí)偏移量。本文初始搜索窗口設(shè)為128×128（像元），最終搜索窗口為64×64（像元），步長(zhǎng)設(shè)為2（像元），最終獲得30 m分辨率的東西向和南北向偏移量。根據(jù)光學(xué)影像特征追蹤方法計(jì)算出偏移量的結(jié)果存在粗大誤差，首先將結(jié)果中大于3倍標(biāo)準(zhǔn)差的異常值剔除，此外，利用COSI-Corr軟件中的非局部均值濾波對(duì)冰川流速結(jié)果進(jìn)行濾波處理，非局部均值濾波對(duì)圖像中具有相同性質(zhì)的區(qū)域進(jìn)行分類并加權(quán)平均達(dá)到降噪效果，本文將搜索區(qū)域的窗口大小設(shè)為21×21像元，將鄰域塊大小設(shè)為7×7像元，噪聲參數(shù)H設(shè)為1.6倍的標(biāo)準(zhǔn)差。最后根據(jù)濾波后的東西向和南北向偏移量以及影像對(duì)間隔時(shí)間獲取音蘇蓋提冰川年平均流速。

圖2 光學(xué)影像特征方法提取冰川表面流速流程圖Fig.2 Flow chart showing how to derive the glacier surface flow velocity by using optical image feature tracking

由于音蘇蓋提冰川位于地形復(fù)雜的高海拔偏遠(yuǎn)山區(qū)，受到惡劣的自然條件限制，很難獲得實(shí)測(cè)數(shù)據(jù)來(lái)對(duì)遙感影像所獲取的表面流速進(jìn)行驗(yàn)證。理論上，進(jìn)行互相關(guān)計(jì)算時(shí)，冰川穩(wěn)定區(qū)域（冰川內(nèi)和周邊基巖露頭）不會(huì)產(chǎn)生偏移量，因此，可利用冰川穩(wěn)定區(qū)域的偏移量來(lái)評(píng)估冰流結(jié)果的不確定性。其公式如下：

式中：eoff為冰川穩(wěn)定區(qū)偏移量的誤差；MED為冰川穩(wěn)定區(qū)偏移量的平均值；SE為冰川穩(wěn)定區(qū)偏移量的標(biāo)準(zhǔn)誤差，其具體計(jì)算

式中：STDV為冰川穩(wěn)定區(qū)偏移量平均值的標(biāo)準(zhǔn)差；Neff為空間非相關(guān)像元個(gè)數(shù)，其計(jì)算公式如下：

式中：Ntotal為冰川穩(wěn)定區(qū)像元的總數(shù)；PS為像元分辨率；D為空間相關(guān)的最大距離，空間相關(guān)的最大距離一般是像元分辨率的20倍［30］。通過(guò)上述方法，本文提取的音蘇蓋提冰川表面流速不確定性約為12 m·a-1。

3 結(jié)果分析

本文利用獲取的2018—2021年20對(duì)完全覆蓋音蘇蓋提冰川遙感影像數(shù)據(jù)，采用頻率域互相關(guān)算法，得到了年平均流速圖（圖3～5）。從圖中看出，冰川流速?gòu)谋ㄉ嫌蔚奖ㄏ掠沃饾u減小，冰川中心流速比冰川邊緣流速大，符合冰川運(yùn)動(dòng)規(guī)律。冰川末端由于氣溫高、海拔低和坡度小的原因?qū)е履┒吮趶?qiáng)烈，幾乎監(jiān)測(cè)不到末端流速，所以冰川末端流速接近零［21］。如圖所示，2018年1月到2019年6月，音蘇蓋提冰川南支流（圖1的B-B'）年平均流速（最大達(dá)到441 m·a-1）大于北支流（西）（圖1的AA'）年平均流速（最大達(dá)到138 m·a-1），而在2019年6月到2021年11月，北支流（西）年平均流速（最大達(dá)到1 008 m·a-1）大于南支流年平均流速（最大達(dá)到413 m·a-1）。音蘇蓋提冰川北支流（東）在2021年7月流速明顯增大（流速最大達(dá)到280 m·a-1），可能發(fā)生躍動(dòng)。此外，通過(guò)流速圖發(fā)現(xiàn)，主冰川流速?zèng)]有明顯變化。

圖3 2018年1月—2019年6月音蘇蓋提冰川年平均流速（m·a?1），底圖為SRTM DEM山影效果圖Fig.3 Annual mean velocity（m·a?1）of Skamri Glacier from January 2018 to June 2019 overlaid on the hillshade image of SRTM DEM

為了進(jìn)一步分析冰川流速的時(shí)間變化特征，本文提取了音蘇蓋提冰川北支流（西）中心線流速隨時(shí)間的變化（圖6）。音蘇蓋提冰川北支流（西）流速在2018年1月—2019年6月相對(duì)比較穩(wěn)定，流速最大約為138 m·a-1，而從2019年6月，北支流（西）冰川流速突然增大，流速最大達(dá)到333 m·a-1，流速持續(xù)增大到2020年7月，流速最大達(dá)到1 008 m·a-1［距北支流（西）冰川末端（圖1的A'）4 500 m處］，是2018年1月份流速的7倍。隨后流速開(kāi)始減小，到2021年11月流速減小到300 m·a-1，但此時(shí)流速仍是2018年1月份流速的2倍。根據(jù)北支流（西）中心線流速變化，可推斷音蘇蓋提冰川北支流（西）在2019年6月末開(kāi)始發(fā)生躍動(dòng)，至2021年11月仍處于躍動(dòng)期。

圖6 2018—2021年音蘇蓋提冰川北支流（西）中心線流速圖（中心線位置見(jiàn)圖1的A-A'剖面線）Fig.6 Annual velocity profiles along the centerline from terminus to up-glacier of the north（west）Skamri Glacier from 2018 to 2021（The location of the central line is marked in Fig.1 A-A'）

為了分析研究時(shí)間段內(nèi)冰川末端變化，本文利用Landsat-8影像展示了音蘇蓋提冰川北支流（西）冰川末端變化（圖7）。在2018年8月—2020年8月期間，北支流（西）的冰川末端沒(méi)有明顯變化，在此之后，該支冰川末端緩慢前進(jìn)，并出現(xiàn)了支流末端變寬的現(xiàn)象，至2021年9月，北支流（西）冰川末端較2018年8月前進(jìn)了約320 m。

圖7 2018—2021年音蘇蓋提冰川北支流（西）冰川末端變化Fig.7 The changes of terminus in the north（west）of Skamri Glacier from 2018 to 2021

圖8為音蘇蓋提冰川南支流的中心線流速隨時(shí)間變化圖。如圖所示，2018年1—5月，南支流冰川流速明顯增大，最大流速?gòu)?42 m·a-1增大到441 m·a-1，隨后流速開(kāi)始減小，至2020年1月，流速減小到340 m·a-1，而從2020年1月—2021年2月，冰川流速持續(xù)增大，流速最大約為413 m·a-1，之后流速又開(kāi)始減小，到2021年11月，流速約為299 m·a-1，音蘇蓋提冰川南支流在研究時(shí)段內(nèi)流速一直比較活躍。

圖8 2018—2021年音蘇蓋提冰川南支流中心線流速圖（中心線位置見(jiàn)圖1的B-B'剖面線）Fig.8 Annual velocity profiles along the centerline from terminus to up-glacier of the south Skamri Glacier from 2018 to 2021（The location of the central line is marked in Fig.1 B-B'）

為了分析各支流對(duì)主冰川流速的影響，本文提取了音蘇蓋提冰川主冰川橫剖面線流速（圖9）。通過(guò)T1、T2橫剖面線流速變化發(fā)現(xiàn)，北支流（西）發(fā)生躍動(dòng)之前，主冰川在北支流（西）匯入之前流速?zèng)]有增大反而減小，最大流速?gòu)?5 m·a-1減小到14 m·a-1［圖9（a）］，而當(dāng)北支流（西）發(fā)生躍動(dòng)后，主冰川流速整體變大，并且橫剖面T2流速（最大流速為100 m·a-1）大于T1的流速（最大流速為76 m·a-1）［圖9（c）］，顯然，北支流（西）的躍動(dòng)影響了主冰川的流速；從圖9（a）和9（b）中的T3、T4橫剖面線流速變化發(fā)現(xiàn)，北支流（東）匯入前后對(duì)主冰川流速?zèng)]有明顯變化，但2021年10月，北支流（東）匯入主冰川之后，主冰川最大流速?gòu)?3 m·a-1增大到73 m·a-1，可能跟北支流（東）發(fā)生躍動(dòng)有關(guān)；通過(guò)T5、T6橫剖面線流速變化發(fā)現(xiàn)，南支流匯入后主冰川流速增大，由此可推斷，音蘇蓋提冰川主冰川流速受南支流影響，此外，北支流躍動(dòng)會(huì)導(dǎo)致主冰川流速變大。

圖9 音蘇蓋提冰川主冰川橫剖面線圖，流速為2018年11月—2019年3月（a）、2020年11月—2021年11月（b）、2021年10月—2021年11月（c），T1、T2、T3、T4、T5和T6橫剖面線位置見(jiàn)圖1Fig.9 The cross-sectional diagram of Skamri main glacier shows the annual average velocity from November 2018 to March 2019（a），from November 2020 to November 2021（b），and from October 2021 to November 2021（c）.The positions of T1，T2，T3，T4，T5 and T6 cross-sectional lines are shown in Fig.1

表2為音蘇蓋提冰川主冰川和各支冰川在研究時(shí)間段內(nèi)的流速最大值及其高程分布。主冰川在2018年1—3月最大流速高程分布在4 980 m處，2018年3月—2020年1月最大流速高程上升到5 410 m處，隨后最大流速高程又降低到4 980 m處，在2021年11月，最大流速高程分布約在5 410 m處；南支流在2018年1月—2021年11月最大流速的高程分布相對(duì)比較穩(wěn)定，高程分布在4 860～4 960 m之間；北支流（西）在2018年1月—2020年7月最大流速的高程從5 058 m逐漸降低到5 001 m，隨后又開(kāi)始上升，在2021年11月，最大流速的高程約為5 106 m；北支流（東）在2018年1月—2021年11月，最大流速的高程分布在5 150～5 300 m之間。

表2 音蘇蓋提冰川主冰川和各支冰川流速最大值及其高程分布Table 2 Maximum velocity and elevation distribution of Skamri main glacier and its branches

圖4 2019年6月—2021年2月音蘇蓋提冰川年平均流速（m·a?1），底圖為SRTM DEM山影效果圖Fig.4 Annual mean velocity（m·a?1）of Skamri Glacier from June 2019 to February 2021 overlaid on the hillshade image of SRTM DEM

圖5 2021年2月—2021年11月年音蘇蓋提冰川年平均流速（m·a?1），底圖為SRTM DEM山影效果圖Fig.5 Annual mean velocity（m·a?1）of Skamri Glacier from February 2021 to November 2021 overlaid on the hillshade image of SRTM DEM

在空間上，主冰川和各支流冰川最大流速高程分布也存在差異。在2018年1—3月和2019年11月—2021年10月期間，北支流（東）最大流速的高程高于主冰川，南支流以及北支流（東）；而在2018年3月—2019年11月和2021年10—11月期間，主冰川最大流速的高程高于其他三條支流冰川。對(duì)于該區(qū)域流速最大值高程分布的時(shí)空差異的原因需要將來(lái)更多的研究去探討。

4 討論

ITS_LIVE數(shù) 據(jù)［31］（The inter-mission time series of land ice velocity and elevation）來(lái)源于NASA“研究環(huán)境制作地球系統(tǒng)數(shù)據(jù)記錄2017”項(xiàng)目（https：//its-live.jpl.nasa.gov/），該數(shù)據(jù)包含高亞洲地區(qū)1985—2018年逐年冰川流速數(shù)據(jù)，本文選取2013—2018年流速數(shù)據(jù)來(lái)研究音蘇蓋提冰川表面運(yùn)動(dòng)。圖10（a）～（f）是基于ITS_LIVE數(shù)據(jù)獲取的2013—2018年的年均流速圖，圖10（g）～（i）是基于Landsat-8數(shù)據(jù)提取的2019—2021年的年平均流速圖。根據(jù)ITS_LIVE流速結(jié)果和年平均流速結(jié)果發(fā)現(xiàn)2013—2016年北支流（西）流速相對(duì)穩(wěn)定，最大流速介于55～74 m·a-1之間，從2017年流速緩慢增大，到2020年流速達(dá)到最大（流速約為811 m·a-1）。通過(guò)年平均流速變化可發(fā)現(xiàn)，北支流（西）在2019年3月—2021年11月期間發(fā)生躍動(dòng)。劉星月等［32］利用2015—2016年Sentinel-1數(shù)據(jù)提取音蘇蓋提冰川表面流速，發(fā)現(xiàn)北支流（西）上部積累區(qū)在2016年5月至6月期間速度突增，最高流速達(dá)534 m·a-1，而本文發(fā)現(xiàn)音蘇蓋提冰川北支流（西）在2016年整體年平均流速都較?。蹐D10（d）］，因此推斷北支流（西）在2016年5月至6月可能是由于雪崩導(dǎo)致的流速增大。

圖10 音蘇蓋提冰川2013—2018年ITS_LIVE流速（a）～（f），2019年3月—2021年11月年平均流速（g）～（i）Fig.10 ITS_LIVE Velocity（a）～（f）from 2013 to 2018，and mean annual velocity（g）～（i）from March 2019 to November 2021 of Skamri Glacier

此外，音蘇蓋提冰川主冰川流速在2013年較大，最大流速達(dá)到172 m·a-1，隨后流速開(kāi)始減小，一直減小到2015年（最大流速約為133 m·a-1），2015年之后流速?zèng)]有明顯變化，流速介于108～133 m·a-1之間。音蘇蓋提冰川北支流（東）在2013—2019年3月流速較穩(wěn)定，最大流速介于33～63 m·a-1之間，從2019年3月流速開(kāi)始增大，2021年11月流速達(dá)到之前的兩倍（最大流速為121 m·a-1），因此可判斷北支流（東）發(fā)生躍動(dòng)。音蘇蓋提冰川南支流在2013—2021年期間流速比較活躍，但沒(méi)有明顯變化，最大流速介于349～358 m·a-1之間。

續(xù)表2

5 結(jié)論

本文基于頻率域互相關(guān)算法提取了音蘇蓋提冰川2018年1月1日—2021年11月16日20期冰川表面流速，分析了音蘇蓋提冰川流速的時(shí)空變化特征，根據(jù)研究結(jié)果得出以下結(jié)論：

（1）音蘇蓋提冰川流速存在空間差異，在2018年1月—2019年6月期間，音蘇蓋提冰川南支流流速遠(yuǎn)大于其北支流（西）流速，而2019年6月北支流（西）發(fā)生躍動(dòng)后，改變了空間流速分布，呈現(xiàn)北支流（西）流速大于南支流流速。

（2）根據(jù)冰川流速隨時(shí)間的變化，本文首次通過(guò)冰川流速發(fā)現(xiàn)北支流（西）和北支流（東）發(fā)生躍動(dòng)，北支流（西）冰川在2019年6月—2021年11月處于躍動(dòng)期，且該支流冰川末端在2020年8月開(kāi)始緩慢向主冰川推進(jìn)，至2021年9月，約向主冰川推進(jìn)320 m；而北支流（東）在2021年7—11月可能發(fā)生躍動(dòng)。

（3）通過(guò)研究時(shí)間段內(nèi)冰川橫剖面線流速變化發(fā)現(xiàn)，音蘇蓋提冰川主冰川流速在南支流匯入后明顯增大，此外，兩條北支流的躍動(dòng)也會(huì)導(dǎo)致主冰川流速增大。

（4）音蘇蓋提冰川主冰川和各支流冰川最大流速的高程分布在時(shí)間上和空間上均存在差異，將來(lái)需要更多的研究去探討時(shí)空差異的原因。