馮君園，蔡強(qiáng)國，李朝霞，孫莉英

（1.中國科學(xué)院 地理科學(xué)與資源研究所 陸地水循環(huán)與地表過程重點(diǎn)實(shí)驗(yàn)室，北京100101；2.華中農(nóng)業(yè)大學(xué) 資源與環(huán)境學(xué)院，武漢430070）

我國凍土面積巨大，融水侵蝕是高海拔凍土區(qū)主要侵蝕類型。中國發(fā)生融水侵蝕的區(qū)域主要包括東北、內(nèi)蒙古東部和北部、新疆北部和西部、青藏高原等穩(wěn)定積雪區(qū)和冰川區(qū)，都屬于凍融作用強(qiáng)烈區(qū)域，據(jù)第2次全國土壤侵蝕遙感調(diào)查資料統(tǒng)計(jì)，全國凍土面積約215萬km2，占國土總面積的22.3%左右。景國臣［1］研究表明高海拔凍土區(qū)最主要的侵蝕類型就是融水徑流侵蝕。

融水徑流對(duì)下游河道流量的貢獻(xiàn)很大，融水侵蝕造成的水土流失對(duì)下游河道流量、輸沙有直接影響，危害嚴(yán)重。孫輝等［2］研究表明，青藏高原幾乎所有河流和湖泊在4—6月份解凍期間都比較渾濁。Swift等［3］分析了烏魯木齊河的英雄橋站水文資料，結(jié)果顯示：年徑流量的增加會(huì)使河流的平均含沙量、輸沙量明顯增加。河流泥沙直接影響河床變化，對(duì)河流水情及河流變遷有重大影響，河流輸沙增多引起水庫、渠道淤積，給防洪、灌溉、供水帶來困難。高寒區(qū)嚴(yán)重的融水侵蝕也給青藏鐵路、公路、輸油、輸氣管線的安全運(yùn)行帶來威脅［3－5］。

目前國內(nèi)外融水侵蝕研究較少，對(duì)融水侵蝕缺乏足夠的認(rèn)識(shí)，應(yīng)該重視融水侵蝕研究，包括侵蝕機(jī)理的研究和預(yù)報(bào)模型的開發(fā)。

1 融水侵蝕概念

融水侵蝕是季節(jié)性凍土在低溫環(huán)境反復(fù)的凍融作用下破碎、崩解，表層土壤隨著氣溫上升解凍后，在冰川和積雪融化形成的地表徑流沖刷作用下被剝離、搬運(yùn)和堆積的土壤侵蝕過程。

融水侵蝕是區(qū)別于融雪侵蝕和凍融侵蝕的典型侵蝕類型：（1）融雪侵蝕是融雪徑流作用下的侵蝕過程，主要是對(duì)土壤顆粒進(jìn)行沖刷和搬運(yùn)，同時(shí)凍融過程導(dǎo)致土壤性狀改變，進(jìn)而影響土壤抗蝕能力［6］。融雪侵蝕是融水侵蝕的一種形式，融雪侵蝕側(cè)重于融雪徑流，已有研究割裂了冰川、積雪的聯(lián)系，片面地針對(duì)高緯度地區(qū)融雪徑流作用而定義融雪侵蝕有失偏頗；融水侵蝕的定義包括冰川和積雪融水徑流，考慮我國高海拔寒區(qū)冰川積雪交錯(cuò)分布的情況，融水侵蝕的定義更為全面、合理。（2）目前，國內(nèi)外學(xué)者對(duì)于凍融侵蝕的定義與研究范疇的界定越來越清晰，大部分學(xué)者趨向于將由于溫度的頻繁變化造成的凍融交替所引起的土體或巖石機(jī)械破壞并在重力等作用下被搬運(yùn)、遷移、堆積的過程定義為凍融侵蝕［7］。凍融侵蝕與融水侵蝕關(guān)系緊密但區(qū)別明顯，二者都因前期凍融作用導(dǎo)致土壤破碎，為侵蝕提供豐富的物質(zhì)來源，但是融水侵蝕后期動(dòng)力源自冰川積雪融化后的坡面徑流沖刷，與凍融侵蝕的重力侵蝕力不同。融水侵蝕與凍融侵蝕時(shí)間上同時(shí)發(fā)生、空間上交錯(cuò)分布，兩種侵蝕類型都集中發(fā)生在溫度上升后的春季解凍期，隨著溫度的繼續(xù)升高，夏季仍是融水侵蝕的高發(fā)期；凍融侵蝕高發(fā)區(qū)也是融水侵蝕高發(fā)區(qū)，無法嚴(yán)格區(qū)分兩種侵蝕類型。

2 融水侵蝕過程機(jī)理

融水侵蝕主要受凍融作用及凍土層的影響。季節(jié)性凍土的凍融作用是高海拔寒區(qū)土壤侵蝕的重要?jiǎng)恿χ?，凍土層?duì)融水侵蝕有極大的促進(jìn)作用。

2.1 凍融作用的影響機(jī)理

凍融作用破壞土體結(jié)構(gòu)，影響土壤物理性質(zhì)，降低土壤抗蝕性。凍融作用對(duì)土壤物理性質(zhì)的影響，主要表現(xiàn)在水與冰轉(zhuǎn)化過程中水分體積脹縮對(duì)土壤團(tuán)聚結(jié)構(gòu)的破壞。當(dāng)氣溫下降到一定程度，土壤中的水分開始結(jié)冰。土壤孔隙中的冰晶增長(zhǎng)使土壤團(tuán)聚體受到擠壓或破裂。凍融交替也影響到土壤抗侵蝕性能，主要通過影響土壤結(jié)構(gòu)、土壤顆粒的粘聚力、水熱分布及其表面粗糙度等物理特性進(jìn)而影響土壤可蝕性和徑流侵蝕輸沙能力。凍融交替破壞后的土壤抗蝕性降低，極易發(fā)生融水侵蝕。Ferrick和Gatto［8］通過室內(nèi)試驗(yàn)驗(yàn)證凍融過程顯著增加坡地融凍土壤侵蝕的假設(shè)。結(jié)果表明，凍融循環(huán)顯著增加地表徑流，土壤侵蝕量是未凍土壤的2.4～5.0倍，特別是土壤含水量較高時(shí)更為嚴(yán)重。

凍融作用對(duì)土壤物理性質(zhì)的破壞程度受多因素影響。室內(nèi)試驗(yàn)表明凍融過程降低團(tuán)聚體穩(wěn)定性，破壞程度主要取決于土壤類型、含水量、團(tuán)聚體大小、凍結(jié)溫度及凍融循環(huán)次數(shù)等因素［9］。Formanek等［10］通過室內(nèi)測(cè)量顯示，粉砂壤土內(nèi)聚力強(qiáng)度經(jīng)歷第一次凍融循環(huán)后，比初始值的一半還要小，但是第二次和第三次凍融循環(huán)影響效果很小。Edwards等［11］的室內(nèi)模擬降雨試驗(yàn)結(jié)果顯示裸土融凍后土壤流失增加了90%，這些土壤流失主要是由于凍融土壤入滲降低引起的地表徑流量增加，融水徑流強(qiáng)烈沖刷所致。

2.2 凍土層作用機(jī)理

凍土層對(duì)融水侵蝕的影響是全方面的。土體在凍結(jié)過程中固結(jié)了大量水分，為融水侵蝕的發(fā)生準(zhǔn)備條件；土壤表層解凍后，下層的凍結(jié)層作為不透水層顯著促進(jìn)了融水侵蝕的發(fā)生；土壤解凍深度對(duì)融水侵蝕過程影響也很大。

凍土在形成過程中提高土壤含水量，為融水侵蝕做好前期準(zhǔn)備。凍土在凍結(jié)形成過程中貯存了大量固態(tài)水，土壤含水量顯著增大［12］。青藏高原多年凍土區(qū)10m深度以內(nèi)土層的平均重量含水量為18.1%。估計(jì)由于凍土變化平均每年從青藏高原多年凍土中由地下冰轉(zhuǎn)化成的液態(tài)水資源將達(dá)到50～110億 m3［13］，相當(dāng)于黃河蘭州站年徑流量的1／6～1／3，伴隨著融水徑流的侵蝕作用同樣不可忽視［14］。

不透水凍結(jié)層對(duì)融水侵蝕的影響在于限制水分入滲，從而使表層解凍土壤含水量增大，飽和條件下的土壤抗剪強(qiáng)度降低［15－16］。另外，凍結(jié)層也使融水更容易集流并具有沖刷能力。Sharratt等［17］試驗(yàn)表明，凍結(jié)層限制水的入滲，導(dǎo)致地表徑流和土壤含水量增加，從而增大了侵蝕量。Oygarden［18］通過研究挪威冬季的一次降雨，認(rèn)為未完全解凍層是導(dǎo)致大范圍片蝕和較小細(xì)溝侵蝕發(fā)生的主要原因。

土壤解凍深度對(duì)融水侵蝕過程有顯著影響。王榮軍等［19］分析了天山北坡軍墉湖流域春季融雪期季節(jié)性凍土濕度、雪深及流量數(shù)據(jù)，結(jié)果顯示表層10 cm以下季節(jié)性凍土的濕度劇烈變化造成融雪水下滲增加，削弱了洪峰峰值。范昊明等［20］有關(guān)草甸土融雪侵蝕模擬研究表明，解凍深度越小，前期融雪徑流侵蝕率越大，而后期侵蝕率緩慢增大，解凍深度越大，前期融雪徑流侵蝕率越小，后期細(xì)溝侵蝕出現(xiàn)后侵蝕率急劇增大。

綜上所述，筆者認(rèn)為凍融作用固結(jié)土壤水分，破壞土體，是融水侵蝕的準(zhǔn)備條件；隨著表層土壤解凍，下層不透水的凍結(jié)層限制水分入滲，增大表層土壤含水量，降低土壤抗蝕性，凍結(jié)層也使冰川積雪融水更容易集流并加強(qiáng)其沖刷能力，極大地增強(qiáng)了融水侵蝕強(qiáng)度；隨著解凍層深度增加，后期侵蝕率顯著增大。

3 融水侵蝕影響因素分析

考慮融水侵蝕發(fā)生過程，氣溫、輻射、降水、植被、風(fēng)和地形等是影響融水侵蝕的主要因素，其中氣溫、輻射和降水是影響凍融作用和融水徑流的直接因素，植被、風(fēng)和地形通過影響雪的分布及輻射而間接影響融水侵蝕過程。

融水徑流量是融水侵蝕的首要影響因素，冰川積雪的消融是一個(gè)復(fù)雜過程，從宏觀尺度看，季節(jié)變化與海拔高度是影響冰川積雪融化的主要因素。何元慶等［21］對(duì)玉龍雪山冰川進(jìn)行分析研究，認(rèn)為氣候變暖是玉龍雪山冰川退縮的主要原因。李巧媛［22］研究了升溫情景及降溫情景下青藏高原冰川的變化規(guī)律。研究發(fā)現(xiàn)，青藏高原冰川在氣候變暖情景下表現(xiàn)為不斷退縮并最終可能消亡的變化趨勢(shì)，其變化速度與氣溫變化速率直接相關(guān)。

在較小尺度上，輻射、地形、降雨、植被等是影響積雪融化的主要因素。李晶等［23］分析了2005年6—9月的野外觀測(cè)資料。結(jié)果表明，雪面最主要的能量來源是消融期凈輻射，感熱供熱也是積雪融化不可忽視的因素。焦劍等［24］對(duì)東北地區(qū)融雪期徑流及產(chǎn)沙特征進(jìn)行分析，結(jié)果表明：地貌對(duì)融雪期流域輸沙模數(shù)有顯著影響：丘陵漫崗區(qū)輸沙模數(shù)較大，平均是山區(qū)的29倍，融雪輸沙模數(shù)與流域面積有十分顯著的冪函數(shù)遞減關(guān)系。史彥江等［25］對(duì)伊犁河谷緩坡地融雪侵蝕特征進(jìn)行研究，結(jié)果表明：坡度與融雪強(qiáng)度呈正相關(guān)關(guān)系；融雪后，表層土壤含水量表現(xiàn)為：8°坡顯著低于5°，3°坡；日產(chǎn)沙量和泥沙侵蝕量與坡度呈正相關(guān)。Danny等［26］對(duì)美國俄勒岡州研究表明，在一次同時(shí)有降雨發(fā)生的融雪過程中，沒有林草覆蓋區(qū)域積雪融化雪水當(dāng)量是有較好林草覆蓋區(qū)域的5倍。林草覆蓋減緩了融雪速度，降雨對(duì)積雪融化具有重要的加速作用。Pomeroy等［27］對(duì)加拿大北部Yukon Territory地區(qū)灌木凍土帶研究發(fā)現(xiàn)，較高灌木區(qū)域冬季積雪數(shù)量雖然高出稀疏灌木區(qū)的147%，但融雪速率也快于稀疏灌木區(qū)47%，研究結(jié)果認(rèn)為是出露的灌木反射率低，灌木吸收太陽總能量高于積雪出露狀態(tài)，并把能量傳遞到地表積雪，從而加速融雪速率。灌木能很好實(shí)現(xiàn)這種能量傳遞，但高大喬木的能量傳遞較為有限。

已有研究大多利用一種分析方法對(duì)某一區(qū)域進(jìn)行研究，相似研究很多，結(jié)果不盡相同，且研究區(qū)域局限，研究結(jié)果缺乏廣泛的適用性。融水侵蝕包括復(fù)雜的凍融過程和融水徑流作用，具體影響因素需要開展室內(nèi)模擬試驗(yàn)對(duì)比研究、分析各影響因素的作用機(jī)制。

4 融水侵蝕評(píng)價(jià)方法探討

融水侵蝕的評(píng)價(jià)方法比較少，參考已有研究成果，融水侵蝕評(píng)價(jià)可以借鑒GIS的定性評(píng)價(jià)方法和評(píng)價(jià)模型定量評(píng)估方法。

4.1 GIS定性評(píng)價(jià)

GIS評(píng)價(jià)方法簡(jiǎn)便易行，適用于多種侵蝕類型的宏觀評(píng)價(jià)。朱惇等［28］采用NEFCLASS模型通過地面監(jiān)測(cè)調(diào)查數(shù)據(jù)建立客觀、全面的評(píng)價(jià)標(biāo)準(zhǔn)，并結(jié)合GIS實(shí)現(xiàn)對(duì)研究區(qū)域土壤侵蝕的快速評(píng)價(jià)，經(jīng)驗(yàn)證結(jié)果具有較高可靠性。李成六［7］基于GIS采用因子比較法對(duì)三江源區(qū)的凍融侵蝕強(qiáng)度進(jìn)行了評(píng)價(jià)，結(jié)果顯示凍融侵蝕強(qiáng)烈的區(qū)域主要分布在地形起伏較大的山區(qū)。倪晉仁等［29］提出基于遙感監(jiān)測(cè)數(shù)據(jù)和最小圖斑，建立了土壤侵蝕快速評(píng)估模型。張建國［30］選取氣溫年較差等6個(gè)因子建立凍融侵蝕分級(jí)評(píng)價(jià)指標(biāo)體系，用加權(quán)加和的方法建立了適合西藏自治區(qū)的凍融侵蝕相對(duì)分級(jí)評(píng)價(jià)模型，并在地理信息系統(tǒng)軟件的支持下實(shí)現(xiàn)了西藏自治區(qū)凍融侵蝕相對(duì)分級(jí)。吳萬貞等［31］通過野外調(diào)查建立解譯標(biāo)志分析，結(jié)合凍融發(fā)生的條件對(duì)三江源地區(qū)凍融侵蝕動(dòng)力進(jìn)行分析并對(duì)其進(jìn)行了凍融侵蝕強(qiáng)度的分級(jí)，結(jié)果表明：三江源地區(qū)凍融侵蝕具有分布范圍廣，不同強(qiáng)度凍融侵蝕空間分布明顯，凍融侵蝕地區(qū)分布差異明顯等特點(diǎn)。張亦超［32］提出基于GIS技術(shù)構(gòu)建土壤風(fēng)蝕模型軟件的思路和方法，該軟件系統(tǒng)可以從時(shí)間和空間尺度上快速得到風(fēng)蝕模數(shù)及風(fēng)蝕強(qiáng)度分級(jí)結(jié)果。

綜上所訴，基于GIS的研究主要是采用因子比較法對(duì)侵蝕強(qiáng)度進(jìn)行分級(jí)評(píng)價(jià)，該方法選取因子主觀性太強(qiáng)，且各因子之間存在效果疊加作用，定性結(jié)果可信度不高。同時(shí)該方法需要實(shí)地調(diào)查資料作為結(jié)果驗(yàn)證的對(duì)照，需要大量的野外調(diào)查工作，耗費(fèi)大量的時(shí)間、人力和物力。

4.2 評(píng)價(jià)模型定量評(píng)估

融水侵蝕評(píng)價(jià)模型研究較少，沒有針對(duì)的預(yù)報(bào)模型，可以對(duì)其他模型進(jìn)行修訂、校正實(shí)現(xiàn)融水侵蝕定量評(píng)估，常用的修訂模型有 USLE（Universal Soil Loss Equation）和 WEPP（Water Erosion Prediction Project）模型。

USLE模型研究較多，高峰等［33］基于GIS和USLE模型對(duì)欽江流域土壤侵蝕進(jìn)行定量評(píng)估，并得出欽江流域年均土壤侵蝕模數(shù)為2 608.87t／（km2·a）。國內(nèi)外有針對(duì)USLE模型適用于冬季過程的修訂研究以實(shí)現(xiàn)融雪侵蝕的評(píng)估［32］，但是由于經(jīng)驗(yàn)?zāi)Ｐ突静豢紤]融雪侵蝕的物理過程，經(jīng)過反復(fù)修訂后可能在本地應(yīng)用較好，但若進(jìn)行其他區(qū)域應(yīng)用，則需進(jìn)行大量修訂、校正，沒有廣泛適用性［1］。

WEPP模型起步較晚，但應(yīng)用廣泛。代龍華等［34］基于WEPP模型研究紫色土坡面水蝕，結(jié)果表明該模型不僅預(yù)測(cè)較為準(zhǔn)確，還可估算土壤侵蝕的時(shí)空分布及全坡面或坡面任意一點(diǎn)的凈土壤流失量及隨時(shí)間的變化。WEPP模型對(duì)融雪侵蝕的預(yù)報(bào)則通過空氣溫度、太陽輻射、蒸發(fā)傳熱、降水過程幾個(gè)方面的能量平衡計(jì)算，實(shí)現(xiàn)對(duì)雪的累積、雪密度、積雪是否融化、融雪徑流發(fā)生、土壤凍結(jié)與融化等方面的模擬，進(jìn)而對(duì)融雪侵蝕進(jìn)行以小時(shí)為單位的模擬預(yù)報(bào)［35］。雖然WEPP模型對(duì)融雪侵蝕物理過程的模擬考慮相對(duì)完備，但仍有不足之處，如模型中并未考慮凍融作用對(duì)土壤水分、熱量遷移而造成的土壤可蝕性變化等問題，而這種變化將直接影響到侵蝕量預(yù)測(cè)的精度［1］。

綜上所述，筆者認(rèn)為融水侵蝕評(píng)價(jià)要參考GIS的定性評(píng)價(jià)結(jié)合我國凍融侵蝕的普查標(biāo)準(zhǔn)，進(jìn)行流域的坡面和溝道劃分，建立產(chǎn)流匯水模型；利用已建立的坡面和溝道侵蝕產(chǎn)沙模型實(shí)現(xiàn)流域侵蝕產(chǎn)沙模擬；實(shí)現(xiàn)從坡面到溝道到流域的侵蝕產(chǎn)沙預(yù)報(bào)。

5 研究展望

高海拔寒區(qū)冰川融水、融雪與土壤融凍交互作用，實(shí)際上是一個(gè)水熱耦合、多因素綜合作用的復(fù)雜侵蝕過程，融水侵蝕機(jī)理等各方面研究都需要進(jìn)一步加強(qiáng)和補(bǔ)充。

（1）融水侵蝕的機(jī)理研究較少，也缺乏融水侵蝕模型的開發(fā)研究，參考已有模型的研究，借鑒 WEPP等土壤侵蝕模型方法，以融水坡面侵蝕過程模型為基礎(chǔ)，結(jié)合已建立的氣象模型和土壤融凍模型，采取室內(nèi)模擬試驗(yàn)的侵蝕參數(shù)，利用有限元計(jì)算方法進(jìn)行數(shù)值計(jì)算，進(jìn)行坡面和溝道侵蝕產(chǎn)沙過程的模擬。以GIS為基礎(chǔ)，進(jìn)行流域的坡面和溝道劃分，建立產(chǎn)流匯水模型；利用已建立的坡面和溝道侵蝕產(chǎn)沙模型實(shí)現(xiàn)流域侵蝕產(chǎn)沙模擬；實(shí)現(xiàn)從坡面到溝道到流域的侵蝕產(chǎn)沙預(yù)報(bào)。

（2）融水侵蝕過程是一個(gè)多因素共同影響的過程，了解各因素對(duì)融水侵蝕的影響效應(yīng)和影響機(jī)理是對(duì)融水侵蝕進(jìn)行預(yù)報(bào)和防治的基礎(chǔ)。目前融水侵蝕影響因素缺乏系統(tǒng)的研究，應(yīng)該加強(qiáng)野外監(jiān)測(cè)、室內(nèi)試驗(yàn)等，定量地給出各因素對(duì)融水侵蝕影響的效應(yīng)。

（3）在深入分析融水侵蝕過程機(jī)理和影響因素的基礎(chǔ)上，合理制定融水侵蝕評(píng)價(jià)標(biāo)準(zhǔn)，對(duì)我國融水侵蝕區(qū)侵蝕強(qiáng)度合理分級(jí)，并采取相應(yīng)的防治措施減輕融水侵蝕的危害。

［1］ 景國臣.凍融侵蝕及其形式探討［J］.黑龍江水利科技，2003，31（4）：111－112.

［2］ 孫輝，秦紀(jì)洪，吳楊.土壤凍融侵蝕交替生態(tài)效應(yīng)研究進(jìn)展［J］.土壤，2008，40（4）：505－509.

［3］ Swift D A，Nienow P W，Hoey T B.Basal sediment evacuation by subglacial meltwater：Suspended sediment transport from Haut Glacierd＇Arolla，Switzerland［J］.Earth Surface Processes and Landforms，2005，30（7）：867－883.

［4］ Haritashya U K，Singh P，Kumar N.Suspended sediment from the Gangotri Glacier：Quantification，variability and associations with discharge and air temperature［J］.Journal of Hydrology 2006，321（1）：116－130.

［5］ Singh P，Ramasatri K S，Kumar N.Suspended sediment transport from the dokriani glacier in the Garhwal Himalayas［J］.Nordic Hydrology，2003，34（3）：221－244.

［6］ 范昊明，武敏，周麗麗，等.融雪侵蝕研究進(jìn)展［J］.水科學(xué)進(jìn)展，2013，24（1）：146－152.

［7］ 李成六，馬金輝，唐志光，等.基于GIS的三江源區(qū)凍融侵蝕強(qiáng)度評(píng)價(jià)［J］.中國水土保持，2011（4）：41－43.

［8］ Ferrick M G，Gatto L W.Quantifying the effect of a freeze－thaw cycle on soil erosion：Laboratory experiments［J］.Earth Surface Processes and Landforms，2005，30（10）：1305－1326.

［9］ Six J，Bossuyt H，Degryze S，et al.A history of research on the link between（micro）aggregates soil biota，and soil organic matter dynamics［J］.Soil ＆ Tillage Research，2004，79（1）：7－31.

［10］ Formanek G E，McCool D K，Papendick R I.Freezethaw and consolidation effects on strength of a wet silt loam［J］.Transactions of the ASAE，1984，27（6）：1749－1752.

［11］ Edwards L M，Burney J R.Soil erosion losses under freeze／thaw and winter ground cover using a laboratory rainfall simulator［J］.Canadian Agricultural Engineering，1987，29（2）：109－115.

［12］ 常學(xué)向，王金葉，金博文，等.祁連山林區(qū)季節(jié)性凍土凍融規(guī)律及其水文功能研究［J］.西北林學(xué)院學(xué)報(bào)，2001，16（S1）：26－29.

［13］ 丁永建，潘家華.氣候與環(huán)境變化對(duì)生態(tài)和社會(huì)經(jīng)濟(jì)影響的利弊分析［M］∥秦大河，陳宜瑜，李學(xué)勇.中國氣候與環(huán)境演變（下卷）.北京：科學(xué)出版社，2005.

［14］ 秦大河，丁永建.冰凍圈變化及其影響研究：現(xiàn)狀、趨勢(shì)及關(guān)鍵問題［J］.氣候變化研究進(jìn)展，2009，5（4）：187－195.

［15］ 劉海昆，黃樹祥，王慧.論凍土對(duì)土壤水分動(dòng)態(tài)的影響［J］.黑龍江水利科技，2002，30（3）：87.

［16］ 王金葉，康爾泗，金博文.黑河上游林區(qū)凍土的水文功能［J］.西北林學(xué)院學(xué)報(bào)，2001，16（S1）：30－34.

［17］ Sharratt B S，Lindstrom M J，Benoit G R，et al.Runoff and soil erosion during spring thaw in the northern US corn belt［J］.Journal of Soil and Water Conservation，2000，55（4）：487－494.

［18］ Oygarden L.Rill and gully development during an extreme inter runoff event in Norway［J］.Catena，2003，50（2）：217－242.

［19］ 王榮軍，劉志輝，李誠志，等.融雪期季節(jié)性凍土濕度變化對(duì)融雪洪水的影響［J］.水土保持研究，2013，20（4）：277－281.

［20］ 范昊明，武敏，周麗麗，等.草甸土近地表解凍深度對(duì)融雪侵蝕影響模擬研究［J］.水土保持學(xué)報(bào)，2010，24（6）：28－31.

［21］ 何元慶，章典.氣候變化是玉龍雪山冰川退縮的主要原因［J］.冰川凍土，2004，26（2）：230－231.

［22］ 李巧媛.不同氣候變化情景下青藏高原冰川的變化［D］.長(zhǎng)沙：湖南師范大學(xué)，2011.

［23］ 李晶，劉時(shí)銀，張勇.天山難破科契卡爾巴西冰川消融期雪面能量平衡研究［J］.冰川凍土，2007，29（3）：366－373.

［24］ 焦劍，謝云，林燕，等.東北地區(qū)融雪期徑流及產(chǎn)沙特征分析［J］.地理研究，2009，28（2）：333－344.

［25］ 史彥江，宋鋒惠，羅青紅，等.伊犁河谷緩坡地融雪侵蝕特研究［J］.新疆農(nóng)業(yè)科學(xué)，2009，46（5）：1111－1116.

［26］ Danny M，John K，Dave T，et al.The sensitivity of snowmelt process to climate conditions and forest cover during rain－on－snow：A case study of the1996Pacific Northwest flood［J］.Hydrol Process，1998，12：1569－1587.

［27］ Pomeroy J W，Bewley D S，Essery R L H，et al.Shrub tundra snowmelt［J］.Hydrol Process，2006，20（4）：923－941.

［28］ 朱惇，王天巍，蔡崇法，等.基于模糊神經(jīng)系統(tǒng)與GIS的區(qū)域土壤侵蝕快速評(píng)價(jià)［J］.農(nóng)業(yè)工程學(xué)報(bào)，2010，26（1）：103－107.

［29］ 倪晉仁，李秀霞.基于最小圖斑的土壤侵蝕快速評(píng)估方法［J］.應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2007，15（4）：425－434.

［30］ 張建國，劉淑珍，楊思全，等.西藏凍融侵蝕分級(jí)評(píng)價(jià)［J］.地理學(xué)報(bào)，2006，61（9）：911－918.

［31］ 吳萬貞，劉峰貴.三江源地區(qū)凍融侵蝕動(dòng)力分析及其分布特點(diǎn)［J］.青海師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2010，26（1）：57－61.

［32］ 張亦超，史明昌，岳德鵬，等.基于GIS的土壤風(fēng)蝕模型軟件構(gòu)建［J］.中國水土保持科學(xué)，2013，11（1）：69－74.

［33］ 高峰，華璀，盧遠(yuǎn)，等.基于GIS和USLE的欽江流域土壤侵蝕評(píng)估［J］.水土保持研究，2014，21（1）：18－22.

［34］ 代龍華，曹叔尤，劉興年，等.基于 WEPP模型的紫色土坡面水蝕預(yù)報(bào)［J］.中國水土保持科學(xué)，2008，6（2）：60－65.

［35］ 劉寶元，謝云，張科利，等.土壤侵蝕預(yù)報(bào)模型［M］.北京：中國科學(xué)技術(shù)出版社，2001.