基于不同方法天山西部山區(qū)雪線的確定及其變化規(guī)律分析

湯瑞　穆振俠

摘要：為了更好的指導(dǎo)天山西部山區(qū)冰雪水資源的開發(fā)利用，對(duì)其雪線變化規(guī)律的研究有至關(guān)重要的作用。借助2000-2016年MODIS數(shù)據(jù)、DEM數(shù)據(jù)以及1959-2015年氣象數(shù)據(jù)，采用遙感監(jiān)測(cè)雪線法、氣候雪線法、雪線場(chǎng)法，對(duì)喀什河流域雪線時(shí)空變化規(guī)律進(jìn)行分析研究，結(jié)果表明：（1）基于遙感監(jiān)測(cè)雪線法可以看出，研究區(qū)多年平均8日、年均雪線均呈上升趨勢(shì)；在季節(jié)尺度上，春季和冬季雪線呈上升趨勢(shì)，秋季和夏季雪線呈下降趨勢(shì)；結(jié)合氣候數(shù)據(jù)可知，年均雪線對(duì)降水有較好響應(yīng)，春季雪線對(duì)氣溫有較好響應(yīng)；（2）基于氣候雪線法可以看出，氣候雪線要高于遙感監(jiān)測(cè)雪線，僅從氣候角度分析，氣候變化導(dǎo)致雪線逐年呈上升趨勢(shì)；（3）采用雪線場(chǎng)分析可看出，雪線場(chǎng)的雪線分布主要受地形影響，以2 000 m高度為界限，其上下兩部分雪線場(chǎng)對(duì)降水有不同的響應(yīng)。

關(guān)鍵詞：MODIS；氣候雪線；雪線場(chǎng)；時(shí)空變化；氣候響應(yīng)

中圖分類號(hào)：P343.6文獻(xiàn)標(biāo)志碼：A文章編號(hào)：

16721683（2018）02010807

Abstract：

In order to better guide the development and utilization of ice and snow resources in western Tianshan mountains，it is of great importance to study the law of snowline change.With the help of the MODIS data from 2000 to 2016，the weather data from 1959 to 2015，and the DEM data，we studied the temporal and spatial variation of the snowline in Kashi river basin by using the remote sensing snowline method，the climate snowline method，and the fields of snowline altitude method.The results are as follows： （1） According to the remote sensing snowline method，the multiyear average （8 days） and the annual average snowline showed a rising tendency in the research area；on the seasonal scale，the snowline showed a rising tendency in spring and winter，and showed a declining tendency in autumn and summer.Based on the climatic data，the annual average snowline had a good response to precipitation，and the spring snowline had a good response to temperature.（2） According to the climate snowline method，the climate snowline was higher than the remote sensing snowline.From the climate perspective，the climate change caused the snowline to increase year by year.（3） It is shown from the analysis of the fields of snowline altitude that the distribution of the snowline is mainly influenced by topography.With the 2 000 m height as the dividing line，the two parts of the fields of snowline altitude have different responses to precipitation.

Key words：

MODIS；climate snowline；fields of snowline altitude；temporal and spatial variation；climate response

雪線作為冰凍圈內(nèi)對(duì)氣候最為敏感的要素之一[13]，其變化波動(dòng)直接反映著氣候的變遷[46]，長(zhǎng)時(shí)間序列的雪線也反映著冰川的消融和累積。雪線作為積雪覆蓋區(qū)和無雪區(qū)的交界線[6]，也同時(shí)反映著積雪面積的變化。而隨著全球氣候變暖，各高海拔山脈的積雪呈現(xiàn)出不同趨勢(shì)的減少，對(duì)區(qū)域水循環(huán)產(chǎn)生了極大的影響，所以動(dòng)態(tài)的雪線監(jiān)測(cè)對(duì)水資源利用保護(hù)和水庫的合理調(diào)洪調(diào)度可起到關(guān)鍵作用。特別在新疆這個(gè)氣候干燥，年內(nèi)溫差較大，降雨少、蒸發(fā)強(qiáng)烈的地區(qū)，雪線的研究對(duì)社會(huì)經(jīng)濟(jì)發(fā)展和自然環(huán)境保護(hù)有著更為重要的意義。近些年由于人類活動(dòng)影響，氣候變化劇烈，雪線這一極為敏感的氣候產(chǎn)物的研究得到許多學(xué)者們關(guān)注。目前雪線的研究包含了全球絕大數(shù)主要山脈及高原，研究主要發(fā)現(xiàn)，雪線不同的時(shí)空變化受山體基面[7]、經(jīng)度、緯度、氣溫、降水及下墊面等因素的影響。

雪線的監(jiān)測(cè)主要有兩種途徑：一是由雪線所在位置的氣象要素建立經(jīng)驗(yàn)公式推導(dǎo)雪線高度[8]，或是利用野外數(shù)據(jù)采集提取雪線高度[10]。二是利用遙感影像結(jié)合地理信息系統(tǒng)提取雪線高度[1112]，隨著遙感技術(shù)的發(fā)展，在地面觀測(cè)站稀少的地區(qū)應(yīng)用氣象衛(wèi)星進(jìn)行積雪研究得到廣泛關(guān)注[12]，國(guó)內(nèi)外學(xué)者基于MODIS衛(wèi)星、HJ1/CCD衛(wèi)星和Landsat衛(wèi)星等光學(xué)遙感對(duì)粒雪線高度進(jìn)行提取[1319]，也有基于SAR等微波遙感對(duì)雪線高度進(jìn)行提取[12，20]。在我國(guó)西北地區(qū)地廣人稀，地面監(jiān)測(cè)站網(wǎng)密度極為稀疏，并且國(guó)內(nèi)利用遙感衛(wèi)星的雪線研究大部分集中在青藏高原地區(qū)[21]，新疆地區(qū)內(nèi)的研究還較為稀少，其中研究大多使用單一的方法對(duì)雪線高度進(jìn)行確定，沒有體現(xiàn)出不同方法提取出的雪線高度之間的差異。本文分別借助遙感和經(jīng)驗(yàn)公式的方法從不同角度對(duì)天山西部山區(qū)喀什河流域的雪線進(jìn)行提取，并分析其時(shí)空變化規(guī)律。

1研究區(qū)概況

喀什河流域位于新疆中段天山山脈西部，東經(jīng)81°40′-85°10′、北緯43°25′-44°15′處，集水總面積達(dá)9 590 km2。地勢(shì)由東向西逐漸降低，山峰平均海拔高于2 000 m，最高海拔至4 615 m。流域?yàn)闇貛Т箨懶詺夂?，水汽主要來源于大西洋暖濕空氣的輸送，由于其地處高寒干旱山區(qū)，多數(shù)山峰均分布著季節(jié)性積雪和冰川，導(dǎo)致流域內(nèi)汛期徑流主要源于融雪徑流補(bǔ)給。流域控制斷面為烏拉斯臺(tái)水文站，站點(diǎn)位于東經(jīng)83°7′、北緯43°48′處，海拔為1 440 m，站點(diǎn)處多年平均降雨量為3823 mm，多年平均氣溫為64 ℃，干流主體為SENW走向，下游向西至與鞏乃斯河匯合。根據(jù)天山自然地理分帶，將喀什河劃分為：丘陵帶（海拔<1 600 m）、山地森林帶（海拔1 600～2 000 m）、高山草甸帶和亞高山草甸帶（2 200～3 300 m）以及冰雪帶和終年積雪帶（>3 300 m）[20]。

2數(shù)據(jù)來源及研究方法

為更加合理的分析研究區(qū)雪線變化規(guī)律，指導(dǎo)雪線相關(guān)方面的研究，采用遙感監(jiān)測(cè)雪線法、氣候雪線法以及雪線場(chǎng)法[22]對(duì)流域雪線時(shí)空分布規(guī)律進(jìn)行研究。

2.1數(shù)據(jù)來源

（1）MODIS數(shù)據(jù)。數(shù)據(jù)來源于美國(guó)航天局NASA（https：//reverb.echo.nasa.gov/reverb/）的八日合成MOD10A2積雪產(chǎn)品數(shù)據(jù)，分辨率為500 m×500 m。

（2）DEM數(shù)據(jù)。數(shù)據(jù)來源于美國(guó)航天局NASA（http：//srtm.csi.cgiar.org/SELECTION/inputCoord.asp）的DEM數(shù)據(jù)，分辨率為90 m×90 m。

（3）實(shí)測(cè)氣象數(shù)據(jù)。數(shù)據(jù)來源于新疆水文年鑒中尼勒克氣象站的實(shí)測(cè)數(shù)據(jù)，選取時(shí)間序列為1959-2015年的月均降水和氣溫?cái)?shù)據(jù)。

2.2研究方法

（1）遙感監(jiān)測(cè)雪線法。

基于所選用的分辨率為500 m×500 m的8日合成積雪數(shù)據(jù)MOD10A2以及90 m×90 m DEM數(shù)據(jù)，經(jīng)過投影轉(zhuǎn)換和去除邊界處理以及疊加DEM數(shù)據(jù)，得到研究區(qū)不同時(shí)期8日雪線，即積雪像元、陸地像元的邊界線及其高程數(shù)據(jù)。流域平均雪線的采用如下計(jì)算為：

h=[SX（]∑[DD（]n[]i=1[DD）]hi[]n[SX）][JY] （1）

式中：h為流域平均雪線高度（m）；hi為單個(gè)雪線像元的高度（m）；n為雪線像元個(gè)數(shù)。

（2）氣候雪線法。

氣候雪線是反映大氣固態(tài)降水零平衡的雪線，也稱之為理論雪線[2324]。其在雪線確定過程中只考慮降水和氣溫這兩種氣候因素，不受地形、坡向等下墊面因素影響，故采用水文氣象數(shù)據(jù)進(jìn)行研究區(qū)氣候雪線的確定。據(jù)吳錫浩[23]等學(xué)者研究表明，我國(guó)西部地區(qū)的雪線年均溫（Ta）與年降水量（Pa）、年均溫（Ta）與年固態(tài)降水率（Psa）以及年均溫（Ta）與年降雪量（Sa）之間的相關(guān)關(guān)系如下[2324]：

Ta=-9.22+0.00407Pa[JY] （2）

Ta=-9.22+0.00407Psa[JY] （3）

Ta=-9.67+0.00604Sa[JY] （4）

Sa=Pa×Psa[JY] （5）

據(jù)蔣復(fù)初[24]研究認(rèn)為，當(dāng)山峰高于氣候雪線出現(xiàn)多年積雪時(shí)，積雪的消融和累積主要決取決夏季氣溫，且我國(guó)西北地區(qū)夏季7月均溫最高，約74 ℃，故認(rèn)為達(dá)到氣候雪線區(qū)域的條件為Ta≤0 ℃，Pa≥50 mm，Psa≥60%，T7<7 ℃。

（3）基于空間尺度進(jìn)行雪線變化規(guī)律分析。

為合理的確定雪線在空間的變化規(guī)律，以提取得到的雪線像元為依據(jù)，選用目前應(yīng)用較多、效果較好的Kriging空間插值法進(jìn)行雪線的插值處理，以得到研究區(qū)消除地形影響的空間雪線場(chǎng)。

3結(jié)果與分析

31遙感監(jiān)測(cè)雪線分析

考慮到雪線的高度隨著時(shí)間推移在年內(nèi)和年際間發(fā)生不同變化，[JP+1]基于喀什河流域2000-2016年八日合成MOD10A2積雪產(chǎn)品數(shù)據(jù)，采用遙感監(jiān)測(cè)雪線法得到研究區(qū)遙感監(jiān)測(cè)的雪線數(shù)據(jù)，并進(jìn)行處理分析，得到多年平均8日的最高、平均及最低雪線高度圖（圖2）和融雪期和累積期的期初期末雪線高度圖（圖3）以及年均雪線與降水、氣溫關(guān)系圖（圖4）。

從圖2可以看出，最高雪線、最低雪線、平均雪線都呈先升高后降低的趨勢(shì)。每年1月到3月為積雪穩(wěn)定期，在此時(shí)間段中雪線高度較為穩(wěn)定，平均雪線維持在1 7544 m左右，但由于逆溫層的存在，導(dǎo)致最高雪線和最低雪線在2月底有突變現(xiàn)象，使雪線高度有明顯的抬升；融雪期從每年3月中旬開始，此時(shí)平均雪線為1 6501 m，之后氣溫上升，積雪消融，雪線高度持續(xù)抬升，直到8月底積雪停止消融達(dá)到最高平均雪線，雪線高度為3 4356 m，在此高度以上為流域的永久積雪覆蓋區(qū)；累積期從每年9月初開始，此時(shí)平均雪線為3 2486 m，之后氣溫下降，高海拔地區(qū)積雪積累，雪線高度持續(xù)下降，直到11月中旬積雪累積達(dá)到最大，此時(shí)平均雪線高度為1 5871 m；11月中旬之后，又為積雪穩(wěn)定期，平均雪線穩(wěn)定的維持在1 7528 m。分析最高雪線和最低雪線的差值，得出夏季的雪線高度平均變化范圍為6386 m，而冬季的雪線高度最大變化范圍為1 7754 m，表明冬季雪線高度的變化波動(dòng)要大于夏季雪線高度的變化波動(dòng)。

由圖3可以看出，融雪期初雪線在1 200～2 000 m范圍波動(dòng)，雪線呈上升趨勢(shì)，表明春季融雪開始時(shí)間提前，最高雪線點(diǎn)出現(xiàn)在2008年，雪線高度為1 912 m，反映了2008年融雪開始時(shí)間較其他年份更早；融雪期末雪線高度在3 000～4 000 m范圍波動(dòng)，雪線呈微弱上升趨勢(shì)，表明流域的最高雪線稍有抬升；累積期末雪線高度在700～2 000 m范圍波動(dòng)，雪線呈下降趨勢(shì)，表明積累期結(jié)束時(shí)間提前，累積期末最低雪線點(diǎn)出現(xiàn)在2016年，雪線高度為693 m，反映了2016年累積期較其他年份更早的結(jié)束。累積期末和融雪期初的雪線高度大致相同，在這一階段雪線高度保持穩(wěn)定。

圖4中年均雪線總體趨勢(shì)較為穩(wěn)定，多年平均雪線高度為2 4861 m。逐年雪線高度表現(xiàn)出較為明顯的波動(dòng)特征，波動(dòng)變化可分為五個(gè)階段，前五年為穩(wěn)定階段，雪線維持在2 330～2 530 m之間，之后每三年為一個(gè)階段表現(xiàn)出波動(dòng)上升或下降趨勢(shì)，波動(dòng)幅度隨時(shí)間推移越來越劇烈。結(jié)合年降水資料可以看出，降水量和雪線高度呈明顯反相關(guān)關(guān)系，降水量增多，雪線下降，反之降水減少，雪線抬升。但2001、2012年雪線和降水表現(xiàn)出良好的同步性，結(jié)合年均氣溫發(fā)現(xiàn)，2001年均氣溫（74 ℃）較2007年（67 ℃）高07 ℃，2012年均氣溫（69 ℃）較2011年（68 ℃）高01 ℃，年均氣溫的上升加劇了全年積雪的消融。通過分析得出，年均雪線的變化主要受降水制約，年氣溫也對(duì)年均雪線有影響作用。年均雪線高度最高點(diǎn)出現(xiàn)在2008年，高度為2 6521 m，結(jié)合2008年的年降水總量2855 mm得出，2008年為2000-2016年之間的特枯水年，較少的降水可能直接阻礙到了積雪的累積；最低點(diǎn)出現(xiàn)在2016年，高度為2 1245 m，2016年為特豐水年，僅到6月底降水量就達(dá)3059 mm、7月28日到8月2日降水量達(dá)766 mm，較多的降水推進(jìn)積雪的累積，并結(jié)合MODIS數(shù)據(jù)可以看出全流域積雪覆蓋時(shí)段較長(zhǎng)，導(dǎo)致2016年雪線高度最低。

為了分析不同季節(jié)的雪線高度變化規(guī)律，再依據(jù)季節(jié)尺度進(jìn)行劃分，得到各季節(jié)雪線的高度圖（圖5）。由圖5可以看出，夏季雪線呈微弱下降趨勢(shì)，秋季雪線基本呈穩(wěn)定趨勢(shì)，夏季平均雪線高度為3 3162 m，秋季為2 4543 m，春季雪線和冬季雪線呈明顯上升趨勢(shì)，春季平均雪線高度為2 3439 m，冬季為1 7606 m，冬季雪線上升幅度最大，表明春季積雪消融量逐年增大，冬季低海拔地區(qū)積雪的累積有明顯減少趨勢(shì)?？紤]到春季為高寒山區(qū)融雪的主要時(shí)期，融雪期在分析河流補(bǔ)排關(guān)系及洪水預(yù)報(bào)中有重要意義，故對(duì)春季雪線和氣候關(guān)系進(jìn)行分析研究。從圖6可以看出，春季雪線和春季均溫表現(xiàn)為良好的正相關(guān)關(guān)系，即隨著氣溫上升或下降，雪線呈現(xiàn)抬升或下降的趨勢(shì)，但2008年均氣溫和雪線呈相反性，結(jié)合2007年和2008年的降水分析，2008年降水量（951 mm）較2007年降水量（186 mm）低了909 mm，春季降水量的減少影響到了積雪的消融，使得氣溫和雪線表現(xiàn)為反相關(guān)。降水量和雪線部分年份呈正相關(guān)，而部分年份雪線高度和降水量表現(xiàn)為反相關(guān)。通過綜合分析得出，在研究區(qū)春季雪線的高低主要受氣溫制約，同時(shí)降水也對(duì)其有影響作用。

3.2氣候雪線分析

為分析氣候雪線的具體高度，以100 m高度為

間隔，借鑒已有資料參考研究區(qū)的降水和氣溫在垂直方向上的變化規(guī)律插值得到不同高程帶的年均溫、年降水、年固態(tài)降水率及7月均溫以尋求達(dá)到條件的氣候雪線高度[24]。具體以氣溫遞減率以065 ℃/（100m）為標(biāo)準(zhǔn)[25]，降水計(jì)算公式為：

Pi+1=Pi+Hi×Q[JY]（6）

式中：Pi+1為下100 m區(qū)域高度的降水量（mm）；Pi為此100 m區(qū)域高度的降水量（mm）；Q為降水遞減率（mm/（100m））；Hi為降水遞減率修正系數(shù)。

依據(jù)尼勒克站氣象數(shù)據(jù)，借助式（2）至式（6），依照上述四個(gè)氣候雪線取值參考條件，得出喀什河流域氣候雪線高度的年際變化趨勢(shì)（圖7），并對(duì)比同時(shí)段的遙感雪線，分析兩種雪線高度的變化規(guī)律。由圖7可以看出，氣候雪線總體呈上升趨勢(shì)，多年平均雪線高度為3 0621 m。低值點(diǎn)出現(xiàn)在1972年，雪線高度為2 706 m，高值點(diǎn)出現(xiàn)在2015年，高度為3 506 m，兩者之間差值達(dá)1 400 m，結(jié)合氣溫?cái)?shù)據(jù)得出，不同年份的相同高度處氣溫相差1 ℃，氣候雪線要相差200～300 m。在1972-1973年間及2003-2004年間有較大的雪線高度差，差值均為400 m，1972年氣候雪線高度2 706 m處7月均溫為64 ℃，在1973相同高度2 706 m處7月均溫為91 ℃，相差27 ℃；2003年氣候雪線高度2 806 m處7月均溫為69 ℃，在2004年相同高度2 806 m處7月均溫為92 ℃，相差23 ℃；表明在1972-1973年和2003-2004年這兩組相鄰年份間氣溫變化較其他年份劇烈。對(duì)比2000-2015年的氣候雪線和遙感年均雪線發(fā)現(xiàn)，2000-2015年間氣候雪線和遙感雪線兩者趨勢(shì)相同，均呈上升趨勢(shì)，表明研究區(qū)的氣候變化導(dǎo)致雪線高度上升，積雪消融量增加。對(duì)比兩種雪線的高度，氣候雪線高度要明顯高于遙感雪線，兩者最大相差1 00575 m，結(jié)合遙感雪線高度為實(shí)際觀測(cè)影像數(shù)據(jù)處理得出這一原因，表明在流域中植被、山體坡度及地形等下墊面因素阻礙了氣候變化導(dǎo)致的融雪過程，使得受綜合因素影響的雪線高度要低于只受氣候因素影響的雪線，且兩種雪線各年份差值均不相同也說明下墊面等因素對(duì)雪線的影響隨時(shí)間發(fā)生著動(dòng)態(tài)變化。

3.3雪線場(chǎng)分析

為分析雪線在空間尺度上變化規(guī)律，基于Kriging插值法[26]生成喀什河流域雪線高度場(chǎng)（圖8）。研究區(qū)總體雪線高度趨勢(shì)由西向東逐漸升高，由南向北逐漸升高，沒有表現(xiàn)出緯度經(jīng)度的地帶性，而是表現(xiàn)出明顯山體地形高度分布對(duì)雪線高度抬升的影響，具體地區(qū)的雪線受山體的坡度、坡向等下墊面因素及水汽的運(yùn)輸?shù)挠绊懕憩F(xiàn)出較為復(fù)雜的分布。多年平均流域雪線高度分布在800～4 200 m之間，最低值出現(xiàn)在流域西南部分，最高值出現(xiàn)在流域東部地區(qū)，東西跨度上雪線高度差3 600 m左右；南北跨度上在流域西部差值在1 000～1 600 m之間，雪線高度由南向北抬升。在流域東部差值基本小于1 000 m，雪線高度表現(xiàn)出流域上下邊界處高、中部低的趨勢(shì)。

依天山自然地理分帶對(duì)區(qū)域雪線分析得出：丘陵帶主要分布在流域西南部分，占流域西部大部分面積，面積達(dá)2 51649 km2，此帶地勢(shì)平緩，雪線高度分布在800～1 600 m，雪線等值線總體較稀疏，表明分帶內(nèi)雪線總體上升趨勢(shì)平緩，但分帶南邊界處有一密集等值線區(qū)，雪線高度達(dá)1 800 m。山地森林帶主要分布在流域中南部和西北部分，跨度面積較丘陵帶小，面積達(dá)2 23663 km2，此帶雪線高度分布在1 600～2 200 m，總體雪線變化梯度較丘陵帶大，分帶南部雪線高度達(dá)2 400 m。高山草甸和亞高山草甸帶主要分布在流域北側(cè)和南側(cè)居中部分，所占流域面積最大，面積達(dá)3 29798 km2，在此分帶內(nèi)雪線等值線也是全流域最為密集的，雪線分布在2 200～3 200 m；流域南邊界處附近，為海拔較低的山峰，雪線高度等值線出現(xiàn)較多同心圓狀，此處雪線高度為分帶內(nèi)最高；靠近流域東部，雪線高度變化梯度較大，呈現(xiàn)南北高中部低的趨勢(shì)，變化幅度在2 000～2 200 m之間。冰雪帶和終年積雪帶主要分布在流域最東處和東北邊界處，此分帶所占流域面積最小，僅為1 5759 km2，此帶雪線高度最高，在3 300 m以上，雪線高度等值線密度較高山草甸帶稀疏，分帶東部多為流域高海拔山峰，雪線等值線呈現(xiàn)同心圓狀，雪線高度達(dá)全流域最高，在3 400 m以上。分帶中部?jī)?nèi)有一低值區(qū)，雪線高度在 3 000～3 200 m之間，表明此處為較低高程的河谷地區(qū)。

由圖9及圖10對(duì)比圖8可以看出多年平均雪線場(chǎng)在流域北部較年雪線場(chǎng)更為平滑，進(jìn)一步的消除了地形對(duì)雪線高度場(chǎng)的影響。2003和2013兩年的雪線場(chǎng)總體變化不大，均與地形高度有著相同的變化趨勢(shì)，具體部分表現(xiàn)出在2 000 m海拔以下豐水年雪線高度較高于枯水年的高度，豐水年平均雪線高度為1 61609 m，枯水年為1 48317 m，兩個(gè)年份的雪線高度差為13292 m；在2 000 m以上豐水年雪線高度較低于枯水年的，豐水年平均雪線高度為2 97174 m，枯水年為3 0527 m，高度差為8096 m。結(jié)合氣溫?cái)?shù)據(jù)得出，2 000 m處年均溫在0 ℃左右，為積雪或融雪的臨界溫度。表明在2 000 m以下降水對(duì)積雪起到消融作用，2 000 m以上降水對(duì)積雪的累積起到推進(jìn)作用。

4結(jié)語

為了更好的指導(dǎo)高寒山區(qū)冰雪水資源的開發(fā)，雪線的研究起到至關(guān)重要的作用。因此，借助研究區(qū)遙感數(shù)據(jù)及水文氣象數(shù)據(jù)，采用多種方法對(duì)喀什河流域雪線的時(shí)空變化規(guī)律進(jìn)行分析，得到以下結(jié)論。

（1）借助遙感監(jiān)測(cè)雪線法可以看出，每年1月到3月為積雪穩(wěn)定期；融雪期從每年3月中旬開始，至8月底結(jié)束，達(dá)到流域最高雪線；累積期末結(jié)束在11月中旬。在季節(jié)尺度上，夏季雪線呈微弱下降趨勢(shì)，秋季雪線基本呈穩(wěn)定趨勢(shì)，春季雪線和冬季雪線呈明顯上升趨勢(shì)；在年際尺度上，逐年雪線高度表現(xiàn)出較為明顯的波動(dòng)特征，前五年為穩(wěn)定階段，之后每三年為一個(gè)階段表現(xiàn)出波動(dòng)上升或下降趨勢(shì)。年均雪線對(duì)降水有很好的響應(yīng)，表現(xiàn)為反相關(guān)關(guān)系；春季雪線對(duì)氣溫有很好的響應(yīng)，表現(xiàn)為正相關(guān)關(guān)系。

（2）借助氣候雪線法可以看出，2000-2015年間氣候雪線和遙感雪線均呈上升趨勢(shì)。氣候雪線高度要明顯高于遙感雪線，兩者最大相差1 00575 m，說明在流域中植被、山體坡度及地形等下墊面因素阻礙了氣候變化導(dǎo)致的融雪過程，且兩種雪線各年份差值均不相同也說明下墊面等因素對(duì)雪線的影響隨時(shí)間發(fā)生著動(dòng)態(tài)變化。

（3）通過空間雪線場(chǎng)可以看出，雪線場(chǎng)的總體高度趨勢(shì)由西向東逐漸升高，由南向北逐漸升高，最低值出現(xiàn)在流域西南部分，最高值出現(xiàn)在流域東部地區(qū)，且由自然地理分帶分析得出，在研究區(qū)雪線高度分布主要受地形和下墊面影響。對(duì)比豐水年和枯水年的空間雪線場(chǎng)，以2 000 m為界限，界限以下降水對(duì)積雪起到?jīng)_刷消融作用，界限以上降水加大了積雪的積累。

盡管采用了多種數(shù)據(jù)，借助了多種方法分析了雪線的時(shí)空分布規(guī)律，但由于遙感數(shù)據(jù)分辨率較低，山區(qū)水文氣象站點(diǎn)一般布設(shè)在出山口，沒有足夠的代表性，從而使得對(duì)雪線變化規(guī)律的認(rèn)識(shí)存在不足。隨著高空水文氣象數(shù)據(jù)和更高精度的遙感數(shù)據(jù)等數(shù)據(jù)源的出現(xiàn)，能為以后更準(zhǔn)確的研究雪線分布規(guī)律的起到了重要作用。[HJ1.8mm]

參考文獻(xiàn)（References）：

[1]沈永平，蘇宏超，王國(guó)亞，等.新疆冰川、積雪對(duì)氣候變化的響應(yīng)（Ⅰ）：水文效應(yīng)[J]，冰川凍土，2013，35（3）：531527.（SHENG Y P，SU H C，WANG G Y，et al.The responses of glaciers and snow cover to climate change in Xinjiang（Ⅰ）：Hydrological Effect[J].Journal of Glaciology and Geocryology，2013，35（3）：531527.（in Chinese）） DOI：10.7522/j.issn.10000240.2013.0061.

[2]張國(guó)飛，李忠勤，王文彬，等.天山烏魯木齊河源1號(hào)冰川19592009年物質(zhì)平衡變化過程及特征研究[J].冰川凍土，2012（6）：13011309.（ZHANG G F，LI Z Q，WANG W B，et al.Change processes and characteristics of mass balance of the Urumqi glacier no.1 at the headwaters of the Urumqi River，Tianshan Mountains，During 19592009[J].Journal of Glaciology and Geocryology，2012（6）：13011309.（in Chinese））

[3]張連成，胡列群，李帥，等.基于GIS的天山南北坡雪線分布特征及其影響因素分析[J].水土保持研究，2015，22（3）：259263.（ZHANG L C，HU L Q，LI S，et al.GISBased analysis on the snowline distributional characteristics of south and north slopes of Tianshan mountain and their influencing factors[J].Research of Soil and Water Conservation，2015，22（3）：259263.（in Chinese）） DOI：10.13869/j.cnki.rswc.2015.03.046.

[4]唐志光，王建，梁繼，等.基于MODIS的青藏高原雪線高度遙感監(jiān)測(cè)[J].遙感技術(shù)與應(yīng)用，2015，30（4）：767774.（TANG Z G，WANG J，LIANG J，et al.Monitoring of snowline altitude over the Tibetan Plateau based on MODIS Data[J].Remote Sensing Technology And Application，2015，30（4）：767774.（in Chinese）） DOI：10.11873/j.issn.10040323.2015.4.0767.

[5]楊存建，趙梓健，倪靜，等.基于MODIS數(shù)據(jù)的川西積雪時(shí)空變化分析[J].中國(guó)科學(xué)：地球科學(xué)，2011，41（12）：17431750.（YANG C J，ZHAO Z J，NI J，et al.Temporal and spatial analysis of changes in snow cover in Western Sichuan Based on MOIDS images[J].Science China Earth Sciences，2011，41（12）：17431750.（in Chinese）） DOI：10.1007/s1143001143365.

[6]趙軍，黃永生，師銀芳，等.2000—2012年祁連山中段雪線與氣候變化關(guān)系[J].山地學(xué)報(bào)，2015，33（6）：683689.（ZHAO J，HUANG Y S，SHI Y F，et al.Relationship between snow line change and climate change in the Middle of Qilian Mountains during 20002012[J].Mountain Research，2015，33（6）：683689.（in Chinese）） DOI：10.16089/j.cnki.10082786.000082.

[7]韓芳，張百平，譚靖，等.山體基面高度對(duì)青藏高原及其周邊地區(qū)雪線空間分布的影響[J].地理研究，2014，33（1）：2330.（HAN F，ZHANG B P，TAN J，et al.The Effect of mountain basal elevation on the distribution of snowline with different mountain basal elevations in Tibetan Plateau and its Surrounding Areas[J].Geographical Research，2014，33（1）：2330.（in Chinese）） DOI：10.11821/dlyj201401003.

[8]VERDHEN A，CHAHAR B R，GANJU A，et al.Modeling Snow Line Altitudes in the Himalayan Watershed[J].Journal of Hydrologic Engineering，2016，21（1）：04015056.DOI：10.1061/（ASCE）HE.19435584.0001255.

[9]張連成，胡列群，李帥，等.基于遙感影像的雪線高程提取方法對(duì)比分析——以天山北坡為例[J].干旱區(qū)地理，2015，38（04）：788796.（ZHANG L C，HU L Q，LI S，et al.Comparison of snowline elevation extraction method based on remote sensing image：A case of north slope of Tianshan Mountains[J].Arid Land Geography，2015，38（04）：788796.（in Chinese）） DOI：10.13826/j.cnki.cn651103/x.2015.04.017.

[10]劉耕年，陳藝鑫，張梅，等.喜馬拉雅山佩枯崗日冰川地貌的年代學(xué)、平衡線高度和氣候研究[J].冰川凍土.2011，33（05）：959970.（LIU G N，CHEN Y X，ZHANG M，et al.Glacial landform chronology and environment reconstruction of Peiku Gangri，Himalayas[J].Journal of Glaciology and Geocryology，2011，33（05）：959970.（in Chinese））

[11]史建橋，白淑英，高吉喜，等.近10a西藏高原雪線時(shí)空變化及其與氣象因素關(guān)系分析[J].沙漠與綠洲氣象.2015，9（03）：713.（SHI J Q，BAI S Y，GAO J X，et al.Spatial and temporal variation of snow line and its relationships with meteorological factors in Tibet Plateau in recent 10 years[J].Desert and Oasis Meteorology，2015，9（03）：713.（in Chinese）） DOI：10.3969/j.issn.10020799.2015.03.002.

[12]HUANG L，LI Z，TIAN B，et al.Classification and snow line detection for glacial areas using the polarimetric SAR image[J].Remote Sensing of Environment，2011，115（7）：17211732.doi：10.1016/j.rse.2011.03.004.

[13]KRAJCI P，HOLKO L，PERDIGAO R A P，et al.Estimation of Regional Snowline Elevation （RSLE） from MODIS Images for Seasonally Snow Covered Mountain Basins[J].Journal of Hydrology，2014，519（B）：17691778.DOI：10.1016/j.jhydrol.2014.08.064.

[14]KRAJ P，HOLKO L，PARAJKA J.Variability of snow line elevation，snow cover area and depletion in the main slovak basins in winters 20012014[J].Journal of Hydrology and Hydromechanics，2016，64（1）：1222.DOI：10.1515/johh20160011.

[15]楊俊華，秦翔，吳錦奎，等.祁連山老虎溝流域春季積雪屬性的分布及變化特征[J].冰川凍土，2012，34（05）：10911098.（YANG J H，QING X，WU J K，et al.Distribution and Variation of Spring Snow Cover in Laohugou Watershed of the Qilian Mountains[J].Journal of Glaciology and Geocryology，2012，34（05）：10911098.（in Chinese））

[16]陳安安，陳偉，吳紅波，等.20002013年木孜塔格冰鱗川冰川粒雪線高度變化研究[J].冰川凍土，2014，36（5）：10691078.（CHENG A A，CHENG W，WU H B，et al.The variations of firn line altitude on the Binglinchuan Glacier，Ulugh Muztagh during 20002013[J].Journal of Glaciology and Geocryology，2014，36（5）：10691078.（in Chinese）） DOI：10.7522/j.issn.10000240.2014.0129.[HJ1.7mm]

[17]李帥，侯小剛，鄭照軍，等.基于2001—2015年遙感數(shù)據(jù)的天山山區(qū)雪線監(jiān)測(cè)及分析[J].水科學(xué)進(jìn)展，2017，28（3）：364372.（LI S，HOU X G，ZHEN Z J，et al.Monitoring and analysis of snowline on Tianshan Mountains based on remotely sensed data from 2001 to 2015[J].Advances in Water Science，2017，28（03）：364372.（in Chinese）） DOI：10.14042/j.cnki.32.1309.2017.03.006.

[18]張其兵，康世昌，張國(guó)帥.念青唐古拉山脈西段雪線高度變化遙感觀測(cè)[J].地理科學(xué)，2016，36（12）：19371944.（ZHANG Q B，KANG S C，ZHANG G S，et al.Changes of snow line altitude for glaciers on Westernnyainqentanglha range observed by remote sensing[J].Scientia Geographica Sinica，2016，36（12）：19371944.（in Chinese）） DOI：10.13249/j.cnki.sgs.2016.12.020.

[19]鮑偉佳，劉時(shí)銀，吳坤鵬，等.一種基于MODIS積雪產(chǎn)品的雪線高度提取方法[J].冰川凍土，2017，39（2）：259272.（BAO W J，LIU S Y，WU K P，et al.A method for extracting snow line altitude based on MODIS snow product[J].Journal of Glaciology and Geocryology，2017，39（2）：259272.（in Chinese）） DOI：10.7522/j.issn.10000240.2017.0030.

[20]ARIGONYNETO J，SAURER H，SIM？ES J C，et al.Spatial and temporal changes in drysnow line altitude on the Antarctic Peninsula[J].Climatic Change，2009，94（12）：1933.DOI：10.1007/s1058400995501.

[21]鄧育武.西藏南部冰川系統(tǒng)雪線場(chǎng)的建立及其空間分布特征[D].長(zhǎng)沙：湖南師范大學(xué)，2005.（DENG Y W.Creating fields of equilibrium line altitude and analyzing their space distribution characteristics in the Southern Tibet Glacier System[D].Changsha：Hunan Normal University，2005.（in Chinese））

[22]郭忠明，顧祝軍，吳紅波，等.冰川雪線高度研究進(jìn)展[J].遙感技術(shù)與應(yīng)用，2016，31（4）：645652.（GUO Z M，GU Z J，WU H B，et al.Research progress of glacier snowline altitude[J].Remote Sensing Technology and Application，2016，31（4）：645652.（in Chinese）） DOI：10.11873/j.issn.10040323.2016.4.0645.

[23]吳錫浩，祝一志.中國(guó)氣候雪線的高程計(jì)算及變化規(guī)律[J].科學(xué)通報(bào)，1990，35（6）：451453.（WU X H，ZHU Y Z.Elevation calculation and variation law of climate snowline in China[J].Chinese Science Bulletin，1990，35（6）：451453.（in Chinese））

[24]蔣復(fù)初，吳錫浩，王書兵，等.中國(guó)氣候雪線空間分布特征[J].地質(zhì)力學(xué)學(xué)報(bào)，2002，8（4）：289296.（JIANG F C，WU X H，WANG S B，et al.Characteristics of space distribution of the climatic snowline in China[J].Journal of Geomechanics，2002，8（4）：289296.（in Chinese））

[25]穆振俠.高寒山區(qū)降水垂直分布規(guī)律及融雪徑流模擬研究[D].烏魯木齊：新疆農(nóng)業(yè)大學(xué)，2010.（MU Z X.Research on vertical distribution law of precipitation and snowmelt runoff simulation in high cold alpine areas[D].Urumqi：Xinjiang Agricultural University，2010.（in Chinese））

[26]王利平，謝自楚，丁良福，等.基于GIS的羌塘高原冰川系統(tǒng)雪線場(chǎng)的建立及其空間分布特征[J].干旱區(qū)地理，2010，33（5）：692702.（WANG L P，XIE Z C，DIN L F，et al.Establishing the field of equilibrium line altitude of the glacier system in the Chiangtang Plateau based on geographic information system and its distribution characteristics[J].Arid Land Geography，2010，33（5）：692702.（in Chinese） DOI：10.13826/j.cnki.cn651103/x.2010.05.005.