謝忠俍,肖聲華

(1. 陜西鐵路工程職業(yè)技術(shù)學(xué)院，陜西 渭南 714099; 2. 廣東省大寶山礦業(yè)有限公司，廣東 韶關(guān) 512128)

0 前 言

雪是地表最活躍的自然要素之一，雪的特征(如積雪面積、積雪分布、雪深等)是全球氣候、能量平衡、生態(tài)模型和水文的重要輸入?yún)?shù)。中國(guó)是一個(gè)積雪多并且雪災(zāi)頻繁的國(guó)家，穩(wěn)定的季節(jié)性積雪大概是420 km×104 km左右[1]，其主要分布在內(nèi)蒙古和青藏高原以及新疆北部和東北等地。復(fù)雜的氣候背景和地理環(huán)境及多樣的降雪天氣導(dǎo)致我國(guó)的積雪分布很不均勻，雪蓋特征多樣。積雪特征研究以往大多是基于氣象站和水文站直接觀測(cè)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，但是這種人工作業(yè)方式，容易受到歷史氣象資料缺失及站點(diǎn)數(shù)量稀少的影響，尤其是在西藏自然地理、氣象條件惡劣的情況下更使得地面觀測(cè)數(shù)據(jù)受到限制。近年來(lái)，隨著遙感技術(shù)的不斷發(fā)展，應(yīng)用遙感數(shù)據(jù)進(jìn)行積雪方面的研究取得了很大的進(jìn)步。美國(guó)國(guó)家海洋和大氣局(NOAA)早在1966年第一顆極軌衛(wèi)星發(fā)射以后就開(kāi)始制作北半球雪蓋產(chǎn)品，通過(guò)衛(wèi)星遙感監(jiān)測(cè)分析得到了全球積雪年際變化資料[2]。Frei等[3]通過(guò)遙感手段分析表明北美洲的積雪從20世紀(jì)30年代到80年代初在冬季一直呈現(xiàn)增加趨勢(shì)，隨后積雪呈減少趨勢(shì)。楊修群等[4]對(duì)北半球的雪蓋變化進(jìn)行了時(shí)空特征分析，數(shù)據(jù)采用NOAA提供的1988-1998年北半球逐周的雪蓋數(shù)據(jù)，研究發(fā)現(xiàn)1988-1998年這十年來(lái)，北半球的積雪變化最強(qiáng)烈的異常區(qū)域在我國(guó)的青藏高原。楊存建等[5]通過(guò)實(shí)驗(yàn)構(gòu)建新的積雪提取模型能有效提取積雪，并通過(guò)目視判讀檢驗(yàn),積雪提取效果理想。巴桑[6]等利用多源遙感數(shù)據(jù)分析西藏積雪變化趨勢(shì)，結(jié)果表明近年來(lái)西藏地區(qū)的積雪面積總體呈下降趨勢(shì)，同時(shí)存在明顯的季節(jié)性變化。楊倩[7]等利用MOD10A1通過(guò)積雪覆蓋指數(shù)研究了2000-2010年吉林省積雪時(shí)空變化特征及其與氣候的關(guān)系。楊志剛[8]等利用MODIS數(shù)據(jù)研究了近15年青藏高原積雪覆蓋時(shí)空變化。

西藏地區(qū)地處青藏高原的核心區(qū)域，是南亞和我國(guó)的許多大江大河的發(fā)源地。研究已表明該區(qū)域是我國(guó)積雪面積最大、時(shí)間最長(zhǎng)的地區(qū)之一，其對(duì)長(zhǎng)江流域夏季降水及東亞季風(fēng)及旱澇等災(zāi)害性天氣等會(huì)產(chǎn)生直接的影響[6,9]。文中利用MODIS積雪產(chǎn)品，采用月最大值合成的方式，獲取了西藏2002-2012年的積雪面積時(shí)空變化特征，并分析其積雪時(shí)空分布和變化特征，為西藏經(jīng)濟(jì)建設(shè)、環(huán)境保護(hù)及防災(zāi)減災(zāi)提提供科學(xué)依據(jù)。

2 研究區(qū)概況

西藏自治區(qū)土地面積120多萬(wàn)km2，平均海拔高于4 000 m，大小湖泊共計(jì)1 500多個(gè)。由于夏季西南季風(fēng)和冬季西風(fēng)的交替控制，西藏的氣候類(lèi)型復(fù)雜，干季和雨季的分別也非常明顯，干季一般在每年10月至第二年4月，雨季在每年的5月至9月。西藏的南部邊緣地區(qū)和藏北一線牧業(yè)生產(chǎn)受影響最嚴(yán)重的自然災(zāi)害就是雪災(zāi)。從20世紀(jì)80年代到現(xiàn)在，西藏曾發(fā)生多次重大雪災(zāi)，嚴(yán)重?fù)p害了牧民的生命財(cái)產(chǎn)安全。西藏的降雪南北多而中部少，南部邊緣降水量特點(diǎn)為雙峰型、那曲東北部和中西部則為單峰型[10]。同時(shí)由于高原地形復(fù)雜，基礎(chǔ)資料缺失及氣候獨(dú)特等原因，也給積雪研究帶來(lái)了許多困難。

3 研究數(shù)據(jù)與方法

3.1 研究數(shù)據(jù)

目前，NASA面向全球用戶共提供了陸地、大氣以及海洋三類(lèi)共計(jì)44種的MODIS標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品。在MODIS的各種標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品里，MOD10(Terra提供)和MYD10(Aqua提供)代表的是陸地2級(jí)或者3級(jí)的標(biāo)準(zhǔn)雪蓋產(chǎn)品，它們是由美國(guó)國(guó)家雪冰數(shù)據(jù)中心(NSIDC)的分布式主動(dòng)存檔中心(DAAC)面向全球用戶進(jìn)行免費(fèi)發(fā)布。本文收集的覆蓋西藏自治區(qū)，時(shí)間段為2002-2012年的空間分辨率為500 m的雪蓋產(chǎn)品MOD10A2影像數(shù)據(jù)，其影像生成頻率為每8天一次，它是利用每8景 MOD10A1 影像合成1期得到的產(chǎn)品，這樣合成能較大程度的降低云覆蓋對(duì)積雪識(shí)別的影響。

3.2 研究方法

1)數(shù)據(jù)預(yù)處理

利用 Arcgis10.1、Envi5.4等專(zhuān)業(yè)數(shù)據(jù)處理軟件對(duì) MOD10A2積雪產(chǎn)品進(jìn)行鑲嵌和坐標(biāo)變換、投影轉(zhuǎn)換，選用最鄰近方法進(jìn)行重采樣,再利用西藏矢量邊界對(duì)其進(jìn)行裁剪，最后的數(shù)據(jù)結(jié)果以GEOTIFF格式輸出。

2)影像合成、積雪數(shù)據(jù)提取

選擇將研究區(qū)進(jìn)行每月4景的MOD10A2圖像融合，為了去除云污染對(duì) MOD10A2雪蓋產(chǎn)品分類(lèi)所造成的影響,根據(jù)MOD10A2產(chǎn)品編碼及意義，通過(guò)設(shè)定地表類(lèi)型的優(yōu)先級(jí)(積雪>湖冰>水體>陸地云層無(wú)意義)的方法來(lái)對(duì)積雪面積進(jìn)行合成處理，其合成規(guī)則[11-12]見(jiàn)表1。對(duì)像元數(shù)據(jù)重分類(lèi)并賦值(積雪為1，其他為0)，統(tǒng)計(jì)不同時(shí)期積雪像元數(shù)，按月積雪面積最大合成規(guī)則獲得研究區(qū)積雪面積[13]。

表1 MOD10A2產(chǎn)品合成規(guī)則

注：*0、1、4、11、254和255分別代表傳感器數(shù)據(jù)丟失、不確定、數(shù)據(jù)錯(cuò)誤、夜晚或傳感器停止工作或在極地區(qū)域、傳感器數(shù)據(jù)飽和及填充數(shù)據(jù)。

3)積雪覆蓋時(shí)間指數(shù)計(jì)算

根據(jù)積雪覆蓋信息提取的結(jié)果統(tǒng)計(jì)得到西藏地區(qū)10年內(nèi)每年的月最大積雪面積、月最小積雪面積，據(jù)此得到西藏積雪的時(shí)間變化特征，同時(shí)對(duì)時(shí)間變化特征進(jìn)行年內(nèi)變化和年際變化分析，為分析西藏積雪覆蓋的空間分布特征。根據(jù)文獻(xiàn)[14-15]提出的積雪覆蓋時(shí)間指數(shù) 對(duì)西藏近10年的各像元平均積雪指數(shù)進(jìn)行提取、計(jì)算，具體見(jiàn)公式1。

(1)

4 結(jié)果分析

4.1 積雪的年內(nèi)變化時(shí)間變化特征

在積雪覆蓋的年內(nèi)變化上，西藏積雪面積隨著季節(jié)的更替有明顯的變化，存在有積雪和融雪兩個(gè)時(shí)期，年內(nèi)積雪變化有明顯的兩個(gè)波峰和一個(gè)波谷(見(jiàn)圖1)。第1個(gè)波峰出現(xiàn)于2月，其平均積雪面積達(dá)到58.6萬(wàn)km2，最大積雪面積也達(dá)到81.1萬(wàn)km2。2月以后，隨著春季的來(lái)臨，積雪開(kāi)始慢慢消融，雖然4月最大積雪面積有所上升，但平均積雪面積及最小積雪面積都在2月后逐漸減少。至5月份，隨著夏季的來(lái)臨，積雪面積在快速的減少，直到7月達(dá)到最少，月平均積雪面積僅27.1萬(wàn)km2。8月份后，隨著氣溫的下降，積雪面積又開(kāi)始快速上升，到10月份達(dá)到第2個(gè)波峰，其平均面積達(dá)到47.3萬(wàn)km2，最大積雪面積達(dá)到65.8萬(wàn)km2。10月份以后，積雪面積開(kāi)始小幅下滑。這是由于10月以后西藏進(jìn)入冬季，氣溫極低，降水不充分，在一定程度上限制了積雪的堆積。

圖1 西藏積雪面積特征值年內(nèi)變化

4.2 積雪的年際變化時(shí)間變化特征

在積雪覆蓋的年際變化上，西藏積雪面積在總體上呈現(xiàn)雙峰波動(dòng)的特征，但是波動(dòng)幅度不大(見(jiàn)圖2)，10年內(nèi)西藏的平均積雪面積在39.9萬(wàn)km2左右，約占研究區(qū)總面積的32.2%。年平均積雪面積的兩個(gè)峰值分別出現(xiàn)在2005年和2008年。在2008年，年平均積雪面積達(dá)到42.1萬(wàn)km2，這與2008年冬天我國(guó)出現(xiàn)大范圍冰雪災(zāi)害天氣這個(gè)事實(shí)相符。年平均積雪面積的兩個(gè)波谷出現(xiàn)在2007年和2010年，這兩年的積雪面積僅為38.2萬(wàn)km2和35.5萬(wàn)km2。近10年來(lái)，積雪面積的年最大值出現(xiàn)的月份大多在1月或者2月，而積雪面積年最小值只出現(xiàn)在7月和8月，其中有8年都出現(xiàn)在7月(見(jiàn)表2)。

圖2 西藏積雪面積特征值年際變化表2 西藏積雪面積特征值的年際變化

從研究區(qū)10年的積雪覆蓋率的時(shí)間序列的變化來(lái)看(見(jiàn)圖3)，西藏的積雪覆蓋率周期性變化顯著，每一年的積雪覆蓋率都呈現(xiàn)出年初和年末高，中間時(shí)段低。在10年里，西藏地區(qū)每月的平均積雪覆蓋率基本上都在10%～60%之間，每年的月平均最小覆蓋率基本上都在7月份，約10%左右。積雪覆蓋率的年內(nèi)變化波動(dòng)范圍較大，而年際變化波動(dòng)范圍較小。

圖3 西藏積雪覆蓋率的時(shí)間序列變化

4.3 積雪的空間分布特征

西藏地區(qū)積雪在空間上的分布受地域的影響差異明顯(圖4)。高積雪區(qū)(積雪時(shí)間指數(shù)大于0.6)主要分布在西藏的林芝的大部分地區(qū)、那曲縣東部、昌都西部以及噶爾縣北部、日喀則南部邊緣等地。這樣的分布很大程度上是因?yàn)槲鞑乇辈刻幱谖黠L(fēng)帶的上升區(qū)、南部處于孟加拉灣與印度的洋暖濕氣流的作用區(qū)，而這些地區(qū)的降水很豐富，這就為積雪的發(fā)展提供了有利的條件。而西藏的中部由于受到高山地形阻隔等的影響，水汽輸送較少，且該地區(qū)的海拔較低，溫度較高，不利于積雪的發(fā)展與堆積，所以中部的積雪覆蓋時(shí)間指數(shù)低。

圖4 西藏積雪10年平均積雪時(shí)間指數(shù)圖

5 結(jié) 論

本文主要通過(guò)收集2002-2012年共10年的西藏MOD10A2雪蓋產(chǎn)品，將數(shù)據(jù)進(jìn)行4期合成一景，得到每個(gè)月的積雪覆蓋影像。分析出西藏10年來(lái)積雪面積的年內(nèi)、年際及空間分布特征。其結(jié)論如下：

1)在積雪覆蓋的年內(nèi)變化上，西藏積雪面積隨著季節(jié)的更替有明顯的變化，存在有積雪和融雪兩個(gè)時(shí)期，年內(nèi)積雪變化有明顯的兩個(gè)波峰和一個(gè)波谷，波峰出現(xiàn)于2月和7月，波谷出現(xiàn)于7月。

2)在積雪覆蓋的年際變化上，西藏積雪面積在總體上呈現(xiàn)雙峰波動(dòng)的特征，但是波動(dòng)幅度不大。10年來(lái)，積雪面積的年最大值出現(xiàn)的月份大多在1月或者2月，而積雪面積的年最小值只出現(xiàn)在7月和8月。

3)空間上，西藏地區(qū)積雪受地域影響差異明顯。高積雪區(qū)主要分布在西藏的林芝的大部分地區(qū)、那曲縣東部、昌都西部以及噶爾縣北部、日喀則南部邊緣等地。