劉晶晶 馬春 蘇鵬程

摘要：冰湖潰決是冰川區(qū)突發(fā)的嚴(yán)重山地災(zāi)害之一。目前，全球氣候變暖，冰川普遍退縮，冰湖面積和潰決風(fēng)險(xiǎn)也隨之增大。因此，冰湖潰決預(yù)警方法的研究對(duì)冰川地區(qū)防災(zāi)減災(zāi)有著重要的實(shí)際意義?，F(xiàn)提供了一種綜合考慮氣溫與降雨條件影響的冰湖潰決預(yù)警方法，選擇預(yù)測日前期正積溫逐日增長速度值與預(yù)測日前期30 d累積降雨量作為氣溫和降雨的代表性指標(biāo)，通過對(duì)2010年前發(fā)生的21例潰決事件統(tǒng)計(jì)，結(jié)合鄰近21個(gè)氣象臺(tái)站的日均溫和日降雨量資料分析，建立了冰湖潰決的預(yù)警警戒線和冰湖潰決預(yù)警系統(tǒng)流程。其后，選擇2013年7月5日的然則日阿錯(cuò)冰湖潰決案例對(duì)建立的冰湖潰決預(yù)警系統(tǒng)進(jìn)行驗(yàn)證，發(fā)現(xiàn)預(yù)警效果良好。

關(guān)鍵詞：氣溫;降雨量;冰湖潰決;預(yù)警

中圖分類號(hào)：P694文獻(xiàn)標(biāo)志碼：A

開放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：[FK（W3。5*2][TP劉晶晶.TIF]

Early warning method of glacial lake outburst floods based on temperature and rainfall

LIU Jingjing1，2，3，MA Chun1，3，SU Pengcheng1，3

（1.Key Laboratory of Mountain Hazards and Surface Process & Institute of Mountain Hazards and Environment，Chinese Academy of Sciences & Ministry of Water Conservancy，Chengdu 610041，China;2.State Key Laboratory of Hydraulics and Mountain River Engineering，Sichuan University，Chengdu 610065，China;3.University of Chinese Academy of Sciences，Beijing 100049，China）

Abstract：

Glacial lake outburst floods （GLOFs） are serious disasters in glacial areas.At present，due to global warming，glaciers are generally diminishing;thus increasing the glacial lake area and the outburst risk.Therefore，the research on early warning methods for GLOFs is important to disaster prevention and relief in glacial areas.This paper proposes an early warning method using the temperature and rainfall as indices.We selected the daily growth rate of prior period positive cumulative temperature and the prior thirty-day cumulative rainfall as the representative indices of temperature and rainfall.Based on the statistics of the 21 outburst events before 2010 and the data from the 21 nearby meteorological stations，we built an outburst early warning threshold and a flow chart of the GLOF early warning system.The system was verified by the GLOF event in the Ranzeaco glacial lake on 2013-07-05.

Key words：

temperature;rainfall;glacial lake outburst floods;early warning

聯(lián)合國政府間氣候變化專門委員會(huì)（Intergovernmental Panel on Climate Change，IPCC）《第五次評(píng)估報(bào)告》[1]指出，大氣中溫室氣體尤其是二氧化碳的濃度持續(xù)增加導(dǎo)致1950年后全球溫度出現(xiàn)了明顯上升。在全球變暖的背景下，大部分高山區(qū)的冰川面積和體積有明顯減小，有些小規(guī)模冰川甚至消失，低緯度和中緯度的冰川縮減尤為明顯[2]。中國現(xiàn)代冰川面積59 406 km.2，冰儲(chǔ)量約為5 590 km.3[3]。自小冰期盛時(shí)以來，西部高山區(qū)平均升溫1.3 ℃，冰川的萎縮量相當(dāng)于現(xiàn)代冰川面積的20%， 預(yù)估2030年、2070年和2100年的升溫值分別為0.4～1.2 ℃、1.2～2.7 ℃和2.1～4.0 ℃，屆時(shí)冰川面積將分別減少12% 、28 % 和45 %[3]。

冰湖是冰川運(yùn)動(dòng)的產(chǎn)物，通常由冰川挖蝕成的洼坑和冰磧物堵塞冰川槽谷積水而形成，以冰川作用或以冰川融水為主要補(bǔ)給源[4-5]。冰湖受氣候影響顯著，隨著氣候變暖，冰川大規(guī)模退縮，冰湖水位上升，補(bǔ)給充足的大型冰湖庫容量增大，而相應(yīng)的小型冰湖在數(shù)量上減少[6]。隨著冰湖面積的擴(kuò)張，冰湖發(fā)生潰決的頻率和規(guī)模也隨之增大[7]。冰湖潰決是冰川區(qū)突發(fā)的嚴(yán)重山地災(zāi)害，其導(dǎo)致的洪水和泥石流突發(fā)性強(qiáng)、流量大、破壞性強(qiáng)、波及范圍大，在沿程又會(huì)激發(fā)其他次生災(zāi)害，從而造成人員傷亡與經(jīng)濟(jì)損失，其災(zāi)害損失比一般的暴雨激發(fā)的洪水和泥石流要嚴(yán)重得多[8-9]。根據(jù)姚曉軍等的研究[8]和筆者后期調(diào)查，發(fā)現(xiàn)1935年以來，西藏地區(qū)發(fā)生的且已有記錄的冰湖潰決事件共30起，給冰湖下游的城市鄉(xiāng)村、居民點(diǎn)、公路設(shè)施、水利電力等帶來了巨大的影響，甚至波及周邊國家[10-11]。

鑒于冰湖潰決后帶來的系列風(fēng)險(xiǎn)，冰湖潰決方面的研究受到越來越多研究者的關(guān)注[12]。冰湖潰決的誘發(fā)因素包括冰崩、雪崩、冰滑坡、冰雪強(qiáng)烈消融、高強(qiáng)度持續(xù)降水等外部誘因[13]，以及冰磧壩內(nèi)核冰融化、冰磧壩管涌等內(nèi)部誘因[14]。目前關(guān)于冰湖潰決方面的研究，多是利用遙感手段[15]，提取下墊面特征（包括冰湖特征、其后冰川特征，終磧堤特征等），進(jìn)行冰湖潰決的風(fēng)險(xiǎn)評(píng)估[16-17]。盡管研究者們很早就認(rèn)識(shí)到氣候變化的水熱組合對(duì)于冰湖潰決和泥石流活動(dòng)起著控制作用[18]，例如研究者發(fā)現(xiàn)1988年光謝錯(cuò)冰湖潰決就是因?yàn)榍捌诘某掷m(xù)高溫和冰川融水的潛蝕作用[19]。并且前期研究中，也發(fā)現(xiàn)西藏冰湖潰決事件都發(fā)生于日均溫>0 ℃的日子，且海拔越高的冰湖發(fā)生潰決的時(shí)間越晚[20-21]。但由于缺乏長期連續(xù)的氣候資料，使得目前的研究只能單純的討論某次事件下的氣候背景，而沒法弄清氣候變化對(duì)冰湖潰決的作用機(jī)制，也沒法依據(jù)氣候條件對(duì)冰湖潰決事件進(jìn)行前期預(yù)警。

冰湖潰決的影響因素包括冰川積累面積和厚度、冰湖大小及儲(chǔ)水量、湖堤結(jié)構(gòu)與穩(wěn)定性等，這些影響因素都與氣溫、降雨有著密切關(guān)系，所以本文認(rèn)為氣溫和降雨是觸發(fā)冰湖潰決的影響因素。與氣候波動(dòng)事件相比，冰湖潰決的發(fā)生雖然是一種低概率事件，但基于冰湖潰決發(fā)生的自然特征規(guī)律可以確定，對(duì)于既定的某一危險(xiǎn)冰湖，雖然不一定存在所謂某一確定的能夠觸發(fā)冰湖潰決的氣溫指標(biāo)值與降雨指標(biāo)值，但一定存在一個(gè)與冰湖潰決發(fā)生高度相關(guān)的以某一或某些氣溫與降雨指標(biāo)為變量的冰湖潰決危險(xiǎn)臨界線。超過該臨界線，冰湖潰決發(fā)生的概率高，低于該臨界線，潰決發(fā)生的概率低，由此可以解決冰湖潰決預(yù)測預(yù)警的技術(shù)問題?；诖耍疚膶?萬多組氣溫與降雨長期觀測數(shù)據(jù)及歷次冰湖潰決事件進(jìn)行關(guān)聯(lián)性分析，試圖從眾多涉及氣溫與降雨的氣象指標(biāo)中篩選出具有代表性的指標(biāo)，確定其與冰湖潰決事件之間存在的符合規(guī)律的函數(shù)關(guān)系，為冰湖潰決的預(yù)警工作提供基礎(chǔ)。

1數(shù)據(jù)來源與遴選方法

1.1研究區(qū)與過去冰湖潰決事件

通過對(duì)文獻(xiàn)資料的整理和野外考察搜集，這里統(tǒng)計(jì)了1930后發(fā)生于我國西藏地區(qū)的冰湖潰決事件共30起，詳見表1。

1.2氣溫和降雨量數(shù)據(jù)來源

本文中涉及到的氣溫和降雨數(shù)據(jù)來自于中國氣象數(shù)據(jù)網(wǎng)（http：//data.cma.cn），包括中國地面氣候資料年值數(shù)據(jù)集、中國地面氣候資料日值數(shù)據(jù)集和中國地面國際交換站氣候資料日值數(shù)據(jù)集（V3.0）（更新時(shí)間截止到2015年5月）。

西藏的氣象臺(tái)站密度較內(nèi)陸低，全區(qū)共有71個(gè)縣，但目前只設(shè)立了39個(gè)氣象觀測站點(diǎn)，且氣象臺(tái)站建立時(shí)間晚，缺測數(shù)據(jù)序列長。本方法所需要日均溫與日降雨量觀測數(shù)據(jù)，本應(yīng)當(dāng)選擇冰湖所在區(qū)域的氣象臺(tái)站提供的數(shù)據(jù)。但在實(shí)際工作中，冰湖位于高寒高海拔區(qū)域，這些區(qū)域由于硬件條件限制，均無配套的自身區(qū)域氣象臺(tái)站，無法提供冰湖所在地的長期日均溫觀測數(shù)據(jù)與日降雨量數(shù)據(jù)。

本次分析中氣象數(shù)據(jù)網(wǎng)數(shù)據(jù)提供至2015年5月，所以無名湖的冰湖潰決事件未納入此次分析。并且由于局部氣象臺(tái)站數(shù)據(jù)序列不全和冰湖潰決發(fā)生時(shí)間記錄不清的原因，塔阿錯(cuò)、窮比嚇瑪錯(cuò)、魯惹錯(cuò)、桑旺錯(cuò)、章藏布、次仁瑪錯(cuò)、[JP2]給曲冰湖的潰決事件也未納入本次分析。如下，僅對(duì)其余22例冰湖潰決事件進(jìn)行分析。經(jīng)過39個(gè)站點(diǎn)的對(duì)比和分析，選擇鄰近臺(tái)站（d<200 km）作為冰湖潰決氣溫和降雨數(shù)據(jù)來源，具體冰湖選擇的氣象站點(diǎn)見表2。

1.3氣溫和降雨量數(shù)據(jù)修正

由于冰湖與氣象臺(tái)站之間存在距離和海拔的差異，所以來自于氣象臺(tái)站的數(shù)據(jù)不能直接體現(xiàn)冰湖當(dāng)?shù)氐臍夂驐l件。在進(jìn)行分析之前，需對(duì)各氣象臺(tái)站的數(shù)據(jù)進(jìn)行修正。

1.3.1個(gè)別缺測數(shù)據(jù)處理

個(gè)別測站會(huì)出現(xiàn)個(gè)別日數(shù)據(jù)缺測的情況，這里采用趨勢預(yù)測法對(duì)其進(jìn)行補(bǔ)差。若個(gè)別測站長時(shí)間缺測，則選用其他測站數(shù)據(jù)修正替代或者直接棄用該站點(diǎn)數(shù)據(jù)。例如阿亞錯(cuò)冰湖事件的分析中，發(fā)現(xiàn)1970年定日氣象站數(shù)據(jù)整年缺測，則只能棄用定日測站數(shù)據(jù)，僅采用聶拉木和日喀則測站數(shù)據(jù)。

1.3.2數(shù)據(jù)距離修正

氣象臺(tái)站通常都距冰湖一定距離，所以來自于氣象臺(tái)站的氣溫和降雨量數(shù)據(jù)都需進(jìn)行距離修正。這里主要選擇距離反比權(quán)重法（IDW）對(duì)數(shù)據(jù)進(jìn)行距離修正。距離反比權(quán)重法（IDW）是一種中間插值方法，利用插值點(diǎn)和樣點(diǎn)間的距離為權(quán)重，進(jìn)行加權(quán)平均[22]，公式如下：

1.3.3數(shù)據(jù)海拔修正

冰湖海拔高度大多都高于所在區(qū)域氣象站的海拔，在以氣溫作為指標(biāo)時(shí)，必須考慮海拔高度的差異對(duì)氣溫的影響[23]。由于目前降雨量的海拔修正，無明確的認(rèn)可方法，故僅對(duì)氣溫?cái)?shù)據(jù)進(jìn)行海拔修正，這里引入青藏高原地區(qū)的氣溫垂直遞減率，對(duì)數(shù)據(jù)進(jìn)行修正。依據(jù)青藏高原及其周邊地區(qū)的氣溫垂直遞減率等值線分布圖[24]，對(duì)中國地面氣候資料的溫度進(jìn)行修正，溫度修正公式如下：

式中：TH為修正后冰湖所在海拔的溫度值;T0表示氣象站實(shí)際測得溫度值（℃）;

1.3.4五日滑動(dòng)平均

對(duì)氣溫?cái)?shù)據(jù)進(jìn)行了前期處理后，為了消除不穩(wěn)定的波動(dòng)，顯示出溫度變化的平穩(wěn)性，充分利用熱量資源，減小氣溫突變誤差，進(jìn)一步采用五日滑動(dòng)平均，遴選出第1日均溫≥0 ℃的日期。

T..*1=（T-4+T-3+…+T1）/5，…，T..*i=（Ti-4+Ti-3+…+Ti）/5[JY]（3）

式中：T1是通過五日滑動(dòng)平均法求算出的穩(wěn)定通過界限溫度0 ℃的起始日;i是起始日到潰決日的天數(shù)。

1.4氣溫和降雨量代表性指標(biāo)

1.4.1氣溫的代表性指標(biāo)—正積溫逐日增長速度值（TV）

在前期研究中，發(fā)現(xiàn)冰湖發(fā)生潰決的當(dāng)日，其日均溫都≥0 ℃，即冰湖潰決都發(fā)生在消融日，這與冰湖潰決的觸發(fā)機(jī)制有關(guān)[20]。所以，認(rèn)為日均溫≥0 ℃的時(shí)期，對(duì)于冰湖潰決的發(fā)生是有效時(shí)期?？紤]到積溫的總量和升溫的過程，對(duì)于冰湖潰決觸發(fā)都有積極的作用。因此提出了采用正積溫逐日增長速度值TV作為氣溫核心預(yù)警指標(biāo)，其計(jì)算方法如下。

（1）獲取日均溫?cái)?shù)據(jù)序列{Tn}，以{Tn}為基礎(chǔ)，篩選確定第1個(gè)日均溫≥0 ℃的日期D1，統(tǒng)計(jì)序列{Tn}中自日期D1至預(yù)測日的天數(shù)i。

（2）依式TDC=∑[DD（]n[]i=1[DD）]Ti計(jì)算預(yù)測日前期正積溫累計(jì)值TDC，式中，Ti為日均溫值（℃）;i為D1至預(yù)測日的天數(shù)（d）。

（3）采用冪函數(shù)曲線TDC=Ai.β擬合，得到正積溫累積過程曲線方程（即正積溫逐日增長的過程線方程），并求得過程曲線方程導(dǎo)數(shù)T′DC=（Ai.β）′=Aβiβ-1，其中常數(shù)部分Aβ即為預(yù)測日前期正積溫逐日增長速度值TV，即Aβ=TV。

1.4.2降雨量的代表性指標(biāo)—前期30 d累積降雨量（RDC）

對(duì)于前期降雨量指標(biāo)的選取，分別對(duì)時(shí)間窗口為10 d、20 d、30 d、40 d、50 d、60 d、70 d、80 d與90 d的前期降雨量指標(biāo)進(jìn)行關(guān)聯(lián)性分析，確定前期30 d累積降雨量與冰湖潰決事件間具有更良好的影響關(guān)系，由此最終選取其作為降雨核心預(yù)警指標(biāo)，其計(jì)算方法如下。

（1）獲取日降雨數(shù)據(jù)序列{Rn}、日均溫?cái)?shù)據(jù)序列{Tn};以{Tn}為基礎(chǔ)，篩選確定第1個(gè)日均溫≥0 ℃的日期d，日期d是前期30 d累積降雨量計(jì)算起始日，則以日期d前第30 d為前期30 d降雨量計(jì)算起始日d-30，統(tǒng)計(jì)自d-30至預(yù)測日的前期30日降雨量天數(shù)n。

（2）以{Rn}為基礎(chǔ)，依式RDC=∑[DD（]30[]n=1[DD）]Rn計(jì)算自日期d-30起至預(yù)測日的前期30 d累積降雨量RDC。式中，Rn為日降雨量值（mm）。

2以氣溫和降雨量為指標(biāo)的冰湖潰決預(yù)警方法

2.1[JP2]基于氣溫和降雨量的預(yù)警警戒線的建

通過對(duì)2010年前發(fā)生的21例潰決事件（然則日阿錯(cuò)冰湖潰決事件作為驗(yàn)證事件），涉及17個(gè)冰湖約5萬組日均溫?cái)?shù)據(jù)和降水量數(shù)據(jù)進(jìn)行分析，統(tǒng)計(jì)21例潰決事件的前期正積溫逐日增長速度值（TV）和前期30 d累積降雨量（RDC），見表4。

基于上述的統(tǒng)計(jì)結(jié)果，對(duì)比潰決事件和無潰決事件的前期正積溫逐日增長速度值（TV）和前期30 d累積降雨量（RDC），描點(diǎn)繪制圖1，設(shè)定預(yù)警警戒線。圖1表明我國境內(nèi)記錄的17個(gè)冰湖的21例潰決事件中，有16例潰決事件位于警戒線TV=-0.0193RDC+3.0018之上。這說明使用本發(fā)明提供的預(yù)警方法可以實(shí)現(xiàn)對(duì)危險(xiǎn)冰磧湖發(fā)生的76.2%的潰決事件預(yù)警，其預(yù)警效率很高。但由于冰磧湖是復(fù)雜多因素耦合的現(xiàn)象，所以在氣候波動(dòng)劇烈的情況下，本方法仍存在約20%左右的漏報(bào)率。

2.2冰湖潰決預(yù)警系統(tǒng)流程

依據(jù)前述分析的過程和建立的冰湖潰決預(yù)警線，建立了冰湖潰決預(yù)警系統(tǒng)的流程圖（圖2）。本方法所需要的日均溫與日降雨量觀測數(shù)據(jù)，應(yīng)當(dāng)選擇冰湖所在區(qū)域的氣象臺(tái)站提供的長期數(shù)據(jù)。但在實(shí)際工作中，冰湖位于高寒高海拔區(qū)域（例如我國青藏高原的冰湖海拔均超過3 800 m）這些區(qū)域由于硬件條件限制，幾乎均無配套的自身區(qū)域氣象臺(tái)站，無法提供冰湖所在地點(diǎn)的長期日均溫觀測數(shù)據(jù)與日降雨量數(shù)據(jù)。因此，依照氣象學(xué)常規(guī)方法，本方法中所需要日均溫與日降雨量觀測數(shù)據(jù)通常從鄰近氣象臺(tái)站獲取。鄰近氣象臺(tái)站的選取原則與方法適用相關(guān)的氣象學(xué)原則與常規(guī)方法，獲得的氣象數(shù)據(jù)經(jīng)常規(guī)氣象學(xué)數(shù)據(jù)修正方法修正后（通常是海拔高差修正、水平距離修正等）適用于本發(fā)明方法。[JP2]通常情況下，選擇兩個(gè)及以上鄰近氣象臺(tái)站作為數(shù)據(jù)來源。為保證預(yù)警結(jié)果的有效性和后期的數(shù)據(jù)篩選的可能性，控制中心調(diào)用的氣象臺(tái)站長期溫度與降雨觀測數(shù)據(jù)應(yīng)當(dāng)具備足夠的規(guī)模，即需要獲取足夠時(shí)間長度的日均溫?cái)?shù)據(jù)和日降雨量數(shù)據(jù)形成日均溫?cái)?shù)據(jù)序列和日降雨量數(shù)據(jù)序列。需至少包括120個(gè)觀測日的日均溫?cái)?shù)據(jù)與120個(gè)觀測日的日降雨量數(shù)據(jù)，才能保證后期分析處理的可靠性，且數(shù)據(jù)序列越長，其可靠度越高。

依據(jù)圖2所示流程，可以對(duì)任一冰湖的潰決可能進(jìn)行預(yù)警。如下，選擇了2013年發(fā)生的然則日阿錯(cuò)冰湖潰決事件，參照冰湖潰決預(yù)警系統(tǒng)流程，進(jìn)行事件驗(yàn)證分析。

3然則日阿錯(cuò)案例驗(yàn)證

冰湖然則日阿錯(cuò)（30°28′07″N， 93°31′55″E）位于西藏自治區(qū)嘉黎縣忠玉鄉(xiāng)，于2013年7月5日發(fā)生潰決，此次潰決造成忠玉鄉(xiāng)238戶1 160人不同程度的受災(zāi)，其中49戶房屋被徹底沖毀，部分人員和牲畜失蹤，農(nóng)田遭到破壞，通往嘉黎縣唯一的1座水泥橋被沖毀，損失估計(jì)達(dá)2.7億元[25]。

根據(jù)然則日阿錯(cuò)冰湖的地理位置，選擇的數(shù)據(jù)主要來自于嘉黎、林芝、比如氣象站，參照冰湖

潰決預(yù)警系統(tǒng)流程對(duì)數(shù)據(jù)進(jìn)行修正后納入分析，其中涉及到自2013年1月1日至7月5日的186個(gè)氣象觀測日，經(jīng)正積溫分析首先確定折麥錯(cuò)冰湖自5月16日起存在潰決風(fēng)險(xiǎn)，將風(fēng)險(xiǎn)范圍縮小至51 d，其后計(jì)算求得2013年5月16日至2013年7月5日期間的正積溫逐日增長速度值（TV）和前期30 d累積降雨量（RDC），見表5。

將表5中所列的正積溫逐日增長速度值（TV）和前期30 d累積降雨量（RDC），與前期設(shè)定的潰決警戒線TV=-0.0193RDC+3.0018對(duì)比繪制圖3。從圖3中，可以看到在初篩后納入分析的51 d中，通過預(yù)警線的設(shè)置，將風(fēng)險(xiǎn)日集中到了11個(gè)報(bào)警日。最終結(jié)果顯示冰湖在11個(gè)報(bào)警日之一（2013年7月5日）發(fā)生了潰決。

4結(jié)論

本文提供了一種綜合考慮氣溫與降雨條件影響的冰湖潰決預(yù)警方法，選擇預(yù)測日前期正積溫逐日增長速度值與預(yù)測日前期30 d累積降雨量作為氣溫和降雨的代表性指標(biāo)，通過對(duì)2010年前發(fā)生的21例潰決事件統(tǒng)計(jì)，結(jié)合鄰近21個(gè)氣象臺(tái)站的日均溫和日降雨量資料分析，建立冰湖潰決的預(yù)警警戒線TV=-0.0193 RDC+3.0018和冰湖潰決預(yù)警系統(tǒng)流程。其后，選擇2013年7月5日的然則日阿錯(cuò)冰湖潰決案例對(duì)建立的冰湖潰決預(yù)警系統(tǒng)進(jìn)行驗(yàn)證，在186個(gè)預(yù)測日中劃定出11個(gè)風(fēng)險(xiǎn)日，其中1日發(fā)生潰決，其預(yù)警效果良好。

[11]呂儒仁.西藏泥石流與環(huán)境[M].成都：成都科技大學(xué)出版社，1999：106-136.（LYU R R.Debris flow and environment in Tibet[M].Chengdu：Chengdu Science and Technology University Press，1999：106-136.（in Chinese））

[12]姚治君，段瑞，董曉輝，等.青藏高原冰湖研究進(jìn)展及趨勢[J].地理科學(xué)進(jìn)展，2010，29（1）：10-14.（YAO Z J，DUAN R，DONG X H，et al.The progress and trends of Glacial Lakes research on Qinghai-Tibet Plateau[J].Progress in Geography，2010，29（1）：10-14.（in Chinese）） DOI：10.11820/dlkxjz.2010.01.002.

[13]CLAGUE J J，EVANS S G.A review of catastrophic drainage of moraine-dammed lakes in British Columbia[J].Quaternary Science Reviews，2000，19（17-18）：1763-1783.DOI：10.1016/S0277-3791（00）00090-1.

[14]劉晶晶，程尊蘭，李泳，等.西藏終磧湖潰決形式研究[J].地學(xué)前緣，2009，16（4）：372-380.（LIU J J，CHENG Z L，LI Y，et al.A study of the outburst form of the End-Moraine Lake in Tibet[J].Geoscience Frontiers，2009，16（4）：372-380.（in Chinese）） DOI：10.3321/j.issn：1005-2321.2009.04.037.

[15]王欣，劉時(shí)銀，姚曉軍，等.我國喜馬拉雅山區(qū)冰湖遙感調(diào)查與編目[J].地理學(xué)報(bào)，2010，65（1）：29-36.（WANG X，LIU S Y，YAO Z J，et al.Glacier lake investigation and inventory in the Chinese Himalayas Based on the remote sensing data[J].Journal of Geographical Sciences，2010，65（1）：29-36.（in Chinese）） DOI：10.11821/xb201001004.

[16]王世金，秦大河，任賈文.冰湖潰決災(zāi)害風(fēng)險(xiǎn)研究進(jìn)展及其展望[J].水科學(xué)進(jìn)展，2012，23（5）：735-742.（WANG S J，QIN D H，REN J W.Progress and prospect in risk assessment of hazards from glacier lake outbursts[J].Advances in Water Science，2012，23（5）：735-742.（in Chinese）） DOI：10.14042/j.cnki.32.1309.2012.05.024.

[17]汪宙峰，張廷山，王成武.西藏喜馬拉雅山地區(qū)冰湖潰決的預(yù)測模型及其應(yīng)用研究[J].冰川凍土，2016，38（2）：388-394.（WANG Z F，ZHANG T S，WANG C W.Prediction model and its application for glacial lake outburst in the Himalayas area，Tibet[J].Journal of Glaciology and Geocryology，2016，38（2）：388-394.（in Chinese）） DOI：10.7522/j.issn.1000-0240.2016.0042.

[18]呂儒仁，李德基.西藏波密冬茹弄巴的冰雪融水泥石流[J].冰川凍土，1989，11（2）：148-160.（LV R R，LI D J.Ice-snow-melt water debris flows in the Dongru Longba（Gully）Bomi County，Xizang（Tibet）[J].Journal of Glaciology and Geocryology，1989，11（2）：148-160.（in Chinese））

[19]李德基，游勇.西藏波密米堆冰湖潰決淺議[J].山地研究，1992，10（4）：219-224.（LI D J，YOU Y.Bursting of the Midui Moraine Lake in Bomi，Xizang[J].Mountain Research，1992，10（4）：219-224.（in Chinese））

[20]LIU J J，CHENG Z L，LI Y.The 1988 glacial lake outburst flood in Guangxieco Lake，Tibet，China[J].Natural Hazards and Earth System Sciences，2014，14（11）：3065-3075.DOI：10.5194/nhess-14-3065-2014.

[21]LIU J J，CHENG Z L，SU P C.The relationship between air temperature fluctuation and Glacial Lake Outburst Floods in Tibet，China[J].Quaternary International，2014，321（2）：78-87.DOI：10.1016/j.quaint.2013.11.023.

[22]PRICE D T，MCKENNEY D W，NALDER I A，et al.A comparison of two statistical methods for spatial interpolation of Canadian monthly mean climate data[J].Agricultural and Forest meteorology，2000，101（2）：81-94.DOI：10.1016/s0168-1923（99）00169-0.

[23]ROBESON S M，JANIS M J.Comparison of temporal and unresolved spatial variability in multiyear time-averages of air temperature[J].Climate Research，1998，10（1）：15-26.DOI：10.3354/cr010015.

[24]李巧媛，謝自楚.高原區(qū)氣溫垂直遞減率的分布及其特點(diǎn)分析——以青藏高原及其周邊地區(qū)為例[J].石河子大學(xué)學(xué)報(bào)（自然科學(xué)版），2006，24（6）：719-723.（LI Q Y，XIE Z C.Analyses on the characteristics of the vertical lapse rates of temperature-Taking Tibetan Plateau and its adjacent area as an example[J].Journal of Shihezi University（Natural Science），2006，24（6）：719-723.（in Chinese）） DOI：10.3969/j.issn.1007-7383.2006.06.016.

[25]孫美平，劉時(shí)銀，姚曉軍，等.2013 年西藏嘉黎縣 “7.5” 冰湖潰決洪水成因及潛在危害[J].冰川凍土，2014，36（1）：158-165.（SUN M P，LIU S Y，YAO X J，et al.The cause and potential hazard of glacial lake outburst flood occurred on July 5，2013 in Jiali County，Tibet[J].Journal of Glaciology and Geocryology，2014，36（1）：158-165.（in Chinese）） DOI：10.7522/j.issn.1000-0240.2014.0020.