陳 慧

（1.國家測繪地理信息局第三航測遙感院，四川 成都 610100）

若爾蓋縣土地沙化遙感監(jiān)測

陳 慧1

（1.國家測繪地理信息局第三航測遙感院，四川 成都 610100）

針對沙化嚴重的若爾蓋縣地區(qū)，選取若爾蓋縣麥溪鄉(xiāng)為研究區(qū)，基于2009年和2014年兩期遙感影像數(shù)據(jù)，采用植被蓋度轉(zhuǎn)換法對研究區(qū)草地沙化信息進行提取，并根據(jù)國情普查數(shù)據(jù)檢驗2014期沙地提取精度。結(jié)果表明，麥溪鄉(xiāng)沙化類型以露沙地為主，兩期監(jiān)測中露沙地均超過地區(qū)沙化總面積的80%；2014年沙化面積為31 597．35 hm2，比2009年增加了1 266．91 hm2，流動沙地、固定沙地面積在減小，而半固定沙地和露沙地面積有所增加，沙地有此消彼長之態(tài)勢。

若爾蓋；土地沙化；遙感監(jiān)測；國情普查；植被蓋度轉(zhuǎn)換法

土地沙化是全球十大環(huán)境問題之首，若爾蓋縣作為長江和黃河上游重要的水源涵養(yǎng)地和重要的生態(tài)屏障，在長期自然和人為因素影響下，其沙化現(xiàn)象已日益嚴重，作為黃河重要水源，其水分涵養(yǎng)功能已被嚴重削弱[1-4]。近年來，遙感技術(shù)因其宏觀性、綜合性、周期性等特點，被廣泛應(yīng)用于生態(tài)環(huán)境動態(tài)監(jiān)測中[5]。針對土地沙化問題，國內(nèi)外相關(guān)研究者也提出了多種遙感監(jiān)測方法[6-8]。在自然諸要素中，植被是依賴性最大的一個要素，它對其他要素，如地形、地貌、土壤、水文條件、氣候等的變化最為敏感。因此，植被指標是區(qū)域尺度上沙化評價指標的首選。所以，利用遙感技術(shù)對沙化進行監(jiān)測也多通過植被指標反演[9-10]。

基于遙感影像和第一次地理國情普查數(shù)據(jù)，利用RS技術(shù)與GIS技術(shù)，監(jiān)測研究區(qū)土地沙化信息，并在此基礎(chǔ)上簡要探討該區(qū)域沙化形成的原因，并提出沙地治理措施，對區(qū)域內(nèi)沙化土地治理和維護區(qū)域生態(tài)安全均有重要現(xiàn)實意義。

1 研究區(qū)概況

四川省若爾蓋縣（33°05′N ～34°19′N，102°10′E～103°37′E）位于四川省西北部阿壩藏族羌族自治州（簡稱阿壩州）的最北部，幅員面積10 334 km2，是阿壩州畜牧業(yè)大縣，與甘肅省瑪曲縣、碌曲縣、卓尼縣、迭部縣和四川省阿壩縣、紅原縣、松潘縣、九寨溝縣接壤。自上個世紀60年代開始，若爾蓋地區(qū)開始出現(xiàn)沙化現(xiàn)象，到2000年的30多年間，其草地沙化面積增加307．7%，達到3 636 019 hm2，占區(qū)域總面積的7．25%[11]。原有沙地不斷擴張，一些非沙地向沙地發(fā)展，導(dǎo)致生態(tài)系統(tǒng)失去平衡，土地生產(chǎn)力不斷下降，生態(tài)環(huán)境日趨惡化。

若爾蓋縣氣候?qū)俑咴瓉喓畮О霛駶櫞箨懶约撅L(fēng)氣候，氣候寒冷，常年無夏，無絕對無霜期；1月平均氣溫為-10．6 ℃，7月平均氣溫為10．8 ℃，年平均氣溫為0．6～1．1 ℃；年降水量為650～750 mm，年內(nèi)降水量分配不均，有70%的降水集中在5～7月[1]。區(qū)域冷濕的氣候有利于沼澤形成和發(fā)育，是我國沼澤率最高的地方。區(qū)內(nèi)主要河流有黃河、黑河、熱曲和白龍江。四川省第一次地理國情普查成果顯示，若爾蓋縣地表覆蓋狀況是以草地為主，占69．1%，其次是林地占27．2%，未利用地（荒漠與裸露地表）占1．2%，見圖1。

圖1 2015年若爾蓋縣地表覆蓋分布圖

麥溪鄉(xiāng)位于若爾蓋縣以西110 km，面積約804 km2。是著名的若爾蓋濕地核心區(qū)，同時也是若爾蓋縣沙化最嚴重的鄉(xiāng)鎮(zhèn)，從圖1也可以看出，麥溪鄉(xiāng)境內(nèi)有大面積的未利用土地。因此本文選取麥溪鄉(xiāng)為研究區(qū)。

2 研究方法

2.1 數(shù)據(jù)源及預(yù)處理

本文數(shù)據(jù)源主要是2009年TM（thematic mapper）影像和2014年OLI（opera tional land imager）影像，時相均為7月，軌道號為131/36和131/37，空間分辨率為30 m，影像質(zhì)量清晰；若爾蓋縣第一次全國地理國情普查數(shù)據(jù)（簡稱國情普查數(shù)據(jù)）現(xiàn)勢性為2014年。

對收集到的遙感影像進行波段合成、幾何精糾正和圖像鑲嵌/裁剪等處理，得到若爾蓋縣兩期遙感鑲嵌裁剪整體影像圖。

2.2 沙地提取

沙化作為土地退化的一種極端形式，其外在表現(xiàn)為植被退化和土壤退化兩種類型。所以，利用遙感技術(shù)對沙化土地進行監(jiān)測也多通過植被或土壤的指標。本文基于兩期遙感影像數(shù)據(jù)，采用蓋度轉(zhuǎn)換法（又稱植被指數(shù)法）估算植被蓋度，再根據(jù)植被蓋度作為沙化分類指標提取相應(yīng)的沙化土地信息。蓋度轉(zhuǎn)換法，即選取與植被覆蓋度有良好相關(guān)性的植被指數(shù)－NDVI（歸一化植被指數(shù)），并通過NDVI與植被覆蓋度的關(guān)系估算植被覆蓋度信息。其中，蓋度轉(zhuǎn)換法公式如下：

式中，F(xiàn)表示植被蓋度；NDVIi表示被求像元的NDVI值；NDVIν和NDVIs分別表示純植被像元和無植被像元的NDVI值[5]。

根據(jù)國情普查數(shù)據(jù)中的解譯樣本信息選取純植被和裸沙進行定位，提取對應(yīng)像元的植被指數(shù)，得到純植被的植被指數(shù)值和裸沙的植被指數(shù)值，然后通過ArcGIS軟件求取研究區(qū)內(nèi)各點的植被指數(shù)，最后代入公式（1）進行計算得到各點植被蓋度值，根據(jù)蓋度分級得到沙地分類信息。

3 沙地分析

沙化土地是在各種氣候條件下，由于各種因素形成的地表呈現(xiàn)以沙（礫）物質(zhì)為主要標志的退化土地。按照§2．2給出的植被蓋度計算方法，得到區(qū)域內(nèi)植被蓋度信息。按照國家林業(yè)局統(tǒng)一的土地沙化分類系統(tǒng)，將兩期沙化土地根據(jù)植被蓋度進行分類，得到流動沙地、半固定沙地、固定沙地、露沙地和沙化耕地，面積統(tǒng)計結(jié)果見表1，空間分布信息見圖2。

3.1 沙地面積變化分析

統(tǒng)計得出區(qū)域內(nèi)2個時期沙地的面積，并計算出它們的增減變化（表1）。

表1 若爾蓋縣麥溪鄉(xiāng)2009～2014年沙地變化/hm2

從表1可以看出，麥溪鄉(xiāng)沙地類型包括流動沙地、半固定沙地、固定沙地和露沙地，沒有沙化耕地。沙地類型以露沙地為主，占全部沙化土地的80%以上。其余沙化類型面積相對較少。露沙地屬于沙化的初級階段，這類沙地植被蓋度較低，在有風(fēng)等外力影響下，地表沙化土就會隨風(fēng)飄散，處于沙漠化的前期階段。

2009年研究區(qū)沙地面積為30 330．45 hm2，至2014年增加到31 597．35 hm2，比2009年增加了4．1%，5 a內(nèi)凈增加面積為1 266．91 hm2。沙地面積從2009年占研究區(qū)總面積的37．7%增加到2015年的39．3%，沙地面積呈現(xiàn)出增加的趨勢。

從不同類型沙地上看，流動沙地在監(jiān)測時段內(nèi)減少601．83 hm2，比2009年減少了28．1%。流動沙地土地為完全沙質(zhì)特征，不具備植物生長的基本土壤條件，是沙化治理中最難以治理的類型，在監(jiān)測時段內(nèi)流動沙地大面積減少表明區(qū)域內(nèi)流動沙地治理初顯成效。半固定沙地增加高達1 212．66 hm2，比2009年增加了將近1倍，遠大于流動沙地減少量，反映出研究區(qū)內(nèi)相當部分沙化土地呈現(xiàn)惡化趨勢，有轉(zhuǎn)為流動沙地的可能性。

3.2 沙地空間格局分析

麥溪鄉(xiāng)沙化土地主要分布在黃河和黑河兩岸河谷地帶、興錯湖和幕措干湖區(qū)域以及沼澤邊緣退化地帶見圖2。

從圖2可以看出，麥溪鄉(xiāng)沙化土地呈現(xiàn)出此消彼長的態(tài)勢。例如，位于麥溪鄉(xiāng)東北的幕措干湖區(qū)域，2009年該湖泊已經(jīng)全部干涸，湖床裸露成為新的沙源地，到2014年幕措干湖區(qū)已實現(xiàn)部分區(qū)域蓄水，沙地面積明顯減少；而位于麥溪鄉(xiāng)西南地區(qū)的興措湖地區(qū)由2009年的固定沙地演變成2014年的半固定和流動沙地，危及周邊草場安全。

圖2 若爾蓋縣麥溪鄉(xiāng)沙化土地空間分布

3.3 精度檢驗

為了驗證植被蓋度法提取麥溪鄉(xiāng)沙化土地信息的準確性，現(xiàn)利用國情普查數(shù)據(jù)對沙地提取結(jié)果進行精度檢驗。計算機自動選取地理國情普查數(shù)據(jù)中沙地點43個，對2014期提取出的沙化土地進行精度驗證，整體平均誤差為15．6%，即分類精度為84．4%，表明該方法在該研究區(qū)提取沙地信息可行。

4 結(jié) 語

本文選取若爾蓋縣沙化土地最為嚴重的麥溪鄉(xiāng)為研究區(qū)，對該區(qū)域5 a來沙化土地的時空變化進行了分析。主要結(jié)論如下：

1）麥溪鄉(xiāng)土地沙化現(xiàn)象嚴重，其沙化面積超過區(qū)域總面積的37%，明顯高于若爾蓋縣平均水平，且土地沙化面積有明顯增加的趨勢；

2）麥溪鄉(xiāng)沙化土地以露沙地為主，占所有沙化土地的80%以上，表明該區(qū)域土地沙化處于初級階段，絕大部分沙地能夠通過科學(xué)治理得以恢復(fù)；

3）流動沙地在監(jiān)測區(qū)內(nèi)大面積減少，共601．83 hm2，比2009年減少近1/3，說明該區(qū)域在治理流動沙地方面獲得了一定的成效；

4）麥溪鄉(xiāng)半固定沙地在監(jiān)測時段內(nèi)顯著增加，表明區(qū)域內(nèi)土地有惡化趨勢，需引起重視；

5）麥溪鄉(xiāng)沙化土地主要分布在河谷地區(qū)、干涸湖區(qū)以及沼澤邊緣退化地帶，且在監(jiān)測時段內(nèi)呈現(xiàn)出此消彼長的態(tài)勢。

針對麥溪鄉(xiāng)嚴重的土地沙化及沙化擴張問題，建議從以下幾方面加以遏制和治理：

1）加強相關(guān)現(xiàn)代高科技知識的學(xué)習(xí)，加大對草地植被修復(fù)技術(shù)的研究和草原蟲鼠害科學(xué)治理等方面的投入，以便為該地區(qū)生態(tài)環(huán)境恢復(fù)提供技術(shù)支持；

2）減少畜牧業(yè)占比，逐步向生態(tài)農(nóng)業(yè)和生態(tài)旅游業(yè)轉(zhuǎn)變，同時啟動生態(tài)補償機制；

3）加大立法保護，做到有法必依、違法必究，以達到減少人為破壞的目的。

[1] 蔣錦剛,李愛農(nóng),邊金虎,等．1974～2007年若爾蓋縣濕地變化研究[J]．濕地科學(xué)，2012(10)：318-326

[2] 吳玉．基于RS與GIS的四川省若爾蓋縣生態(tài)環(huán)境狀況評價[D]．成都：成都理工大學(xué),2010

[3] 王文麗,董治寶,胡光印,等．若爾蓋高原近30年沙地年變化趨勢分析[J]．中國沙漠,2008,28(4)：617-621

[4] 蒙嘉文,左林,蔡應(yīng)君,等．若爾蓋縣土地沙化現(xiàn)狀分析及治理對策研究[J]．四川林業(yè)科技,2013,34(4)：42-46

[5] 陳慧,杜耘,肖飛,等．漢江中游河谷平原植被指數(shù)時空變化及其與沙化土地動態(tài)的關(guān)聯(lián)關(guān)系[J]．長江流域資源與環(huán)境,2013,22(9)：1 221-1 226

[6] 王曉慧．沙化土地遙感監(jiān)測機理和方法研究[D]．北京：中國林業(yè)科學(xué)研究院,2007

[7] 李金亞．科爾沁沙地草原沙化時空變化特征遙感監(jiān)測及驅(qū)動力分析[D]．北京：中國農(nóng)業(yè)科學(xué)院,2014

[8] 徐瑤．藏北草地退化遙感監(jiān)測與生態(tài)安全評價[D]．成都：成都理工大學(xué),2014

[9] 王正興,劉闖．植被指數(shù)研究進展：從AVHRR-NDVI到MODIS-EVI[J]．生態(tài)學(xué)報,2003,23(5)：979-987

[10] 杜子濤,占玉林,王長耀．基于NDVI序列影像的植被覆蓋變化研究[J]．遙感技術(shù)與應(yīng)用,2008,23(1)：47-50

[11] 雍國偉,石承蒼,邱鵬飛．川西北高原若爾蓋草地沙化及濕地萎縮動態(tài)遙感監(jiān)測[J]．山地學(xué)報,2003,21(6)：758-762

P237

B文章編號：1672-4623（2017）06-0031-03

10.3969/j.issn.1672-4623.2017.06.009

2016-03-04。

項目來源：數(shù)字制圖與國土信息應(yīng)用工程國家測繪地理信息局重點實驗室開放研究基金資助項目（DM2014SC06）。

陳慧，碩士研究生，主要從事GIS、RS應(yīng)用及災(zāi)害監(jiān)測方面研究。