甘肅省酒泉市土地沙漠化現狀及動態分析*

馬仲武，王新源， ※，王小軍，丁乾平，曲 浩

(1.甘肅省生態環境監測監督管理局，蘭州 730020； 2.中國科學院西北生態環境資源研究院，甘肅蘭州 730000)

0 引言

沙漠化是干旱、半干旱及部分半濕潤地區在多風干旱和疏松沙質地表條件下，由人類對土地過度開發利用和以風力吹蝕為主的自然營力等因素綜合作用下引起的土地退化過程，使原非沙漠地區呈現出以風沙活動、沙丘起伏為主要標志的類似沙漠景觀的地貌特征[1-2]。沙漠化由于其分布廣泛，危害嚴重，威脅人口眾多，已成為當今世界面臨的重大環境問題之一，引起國際社會的廣泛關注[3-4]。

甘肅省是我國受沙漠化脅迫最為嚴重的省份之一[5]，干旱的氣候、水資源的短缺以及不合理的人類活動引起以河西走廊為代表的部分區域土地沙漠化發展迅速，生態環境趨于脆弱和嚴酷[6]。近年來，針對甘肅省境內土地沙漠化動態變化與分析的研究已有不少，但主要集中于石羊河流域的民勤綠洲[7-9]、黑河中上游的張掖地區[10-11]以及黃河上游的甘南高原[12-13]，而關于酒泉市土地沙漠化現狀及動態分析的研究還較為欠缺。酒泉市地處甘肅省西北部，干旱少雨，水熱條件惡劣，生態環境脆弱[14]，是全省沙漠化最嚴重的區域之一。酒泉市同時也是河西走廊的重要節點城市和“一帶一路”經濟帶的重要黃金通道城市，具有較高的政治經濟地位。此外，酒泉市域范圍內保存有世界文化遺產莫高窟、全國重點文物保護單位玉門關在內的多處重要歷史文化遺跡。因此，掌握該地區沙漠化土地現狀及發展趨勢，明確驅動沙漠化正逆演替的關鍵因素，對區域內沙漠化防治與生態恢復、經濟社會有序持續發展乃至歷史文化遺跡的保存延續均具有重要意義。鑒于此，研究采用遙感(Remote Sensing)技術和人工現地調查相結合的方法，對酒泉市沙漠化現狀、動態變化規律及分布特征進行了監測和分析，旨在為該地區沙漠化防治和生態恢復與重建提供理論依據。

圖1 2009年酒泉市沙漠化土地類型分布

1 研究區概況

酒泉市位于甘肅省西北部，介于北緯38°19′～42°47′、東經92°45′～100°15′之間，面積19.2萬km2，占甘肅省總面積的42.0%，是甘肅省面積最大的市。酒泉市南接青海，西鄰新疆，北接內蒙古并與蒙古國接壤。地勢由東南向西北傾斜，地貌特征復雜多樣，北部為馬鬃山丘陵區，南部屬祁連山西麓，中部有河西走廊疏勒河谷地； 馬鬃山山地風高沙大，形成典型的戈壁景觀； 中部河西走廊疏勒河谷地地勢平坦、光熱充足； 南部祁連山山體巨大，海拔較高，部分山區終年積雪形成特有的冰川地貌。酒泉市大部屬典型的大陸性溫帶干旱區氣候，年日照時數3 033～3 317h，年平均氣溫3.9～9.3℃，多年平均降水量86mm，年蒸發量2 000～4 000mm，年日照時數3 000h以上，平均無霜期118～159d。常住人口110.83萬人， 2013年，全市生產總值達到642.7億元，財政收入85.44億元。據2009年酒泉市土地沙漠化監測結果，酒泉市沙漠化土地面積960.41萬hm2(圖1)。

2 數據影像資料及研究方法

2.1 數據資料

數據源采用甘肅省第五次荒漠化和沙化監測遙感影像數據，主要包括2013年SPOT5 2.5m高分影像、2013年資源一號與資源三號影像(成像季相選擇在植被生物量較大的6～9月份)。同時收集覆蓋全市1： 5萬最新地形圖、土地利用詳查圖、森林分布圖和土壤分布圖等數據資料用以信息提取和完善，對沙漠化動態變化采用文獻資料收集與氣象數據統計[15]相結合的方法進行分析。

2.2 研究方法

2.2.1 遙感圖像處理與解譯標志

應用多波段遙感數據融合方法[16-17]，對各期影像按不同波段進行假色彩合成[3]，采用高斯一克呂格投影進行幾何精校正，從而獲得有利于沙漠化信息提取和識別的基礎影像。之后，依據不同程度沙漠化土地空間分布格局與光譜變化特征[12]，結合現地調查因子(植被、土壤、地形地貌和土地利用信息等)，建立沙漠化土地解譯標志。在此基礎上，根據現地調查獲得的沙漠化土地解譯標志與遙感影像，綜合2009年沙漠化監測結果等輔助信息，采用人—機交互判讀的方式(以SPOT5高分影像為主，資源一號、三號影像為補充參考)，通過ArcGIS9.3軟件對沙漠化信息進行判讀和提取，進而獲得可供分析的矢量化數據。最后，為了核定解譯工作的準確性，研究組于2014年9月5～29日對野外部分區域進行核查，隨機選取200個地塊與解譯結果進行對比，發現其解譯精度達96.5%，符合統計學要求，說明通過計算機軟件分析獲取的沙漠化土地面積、類型、程度和分布數據真實可靠。

2.2.2 沙漠化土地分類

以甘肅省2014年沙漠化監測技術細則(簡稱《細則》)[18]為依據，通過現地調查，綜合風蝕強度、植被蓋度、地貌特征及人類干預等多種因素，將沙漠化土地類型劃分為流動沙地(丘)、半固定沙地(丘)、固定沙地(丘)、露沙地、沙漠化耕地、非生物治沙工程地、風蝕殘丘、風蝕劣地和戈壁9種類型。之后將9種沙漠化土地類型與對應遙感影像的色調、紋理、形狀和分布等特征作對比分析，從而建立準確的解譯標志。

同時依照《細則》[18]將土地沙漠化程度劃分為輕度、中度、重度、極重度4個等級，并通過現地調查后建立對應的遙感解譯標志(表1)。

表1 酒泉市不同程度沙漠化土地景觀特征與解譯標志

2.2.3 年變化率

年變化率反映了研究期內沙漠化土地的變化情況。朱震達等[11]通過計算沙漠化的年增長率來度量沙漠化發展狀況，計算公式為：

(1)

式(1)中，R為沙漠化土地年變化率；n為測算相隔年限；D為測算終結年沙漠化土地面積；A為測算起始年沙漠化土地面積。

3 沙漠化土地分布現狀與變化特征

3.1 沙漠化土地類型分布現狀與變化特征

3.1.1 沙漠化土地類型分布現狀

根據2014年最新監測結果，酒泉市現有沙漠化土地面積956.34萬hm2。其中，流動沙地(丘)115.90萬hm2，占沙漠化土地面積的12.12%； 半固定沙地(丘)102.13萬hm2，占沙漠化土地面積的10.68%； 固定沙地(丘)95.09萬hm2，占沙漠化土地面積的9.94%； 風蝕劣地11.73萬hm2，占沙漠化土地面積的1.23%； 戈壁625.66萬hm2，占沙漠化土地面積的65.42%； 露沙地、沙漠化耕地、非生物治沙工程地和風蝕殘丘5.84萬hm2，占沙漠化土地面積的0.61%(表2)。

表2 酒泉市沙漠化土地類型現狀

沙漠化類型流動沙地半固定沙地固定沙地風蝕劣地戈壁其他類型面積(萬hm2)11590102139509117362566584占沙漠化土地面積比(%)121210689941236542061

由圖2分析可得，酒泉市沙漠化土地類型主要包括流動沙地(丘)、半固定沙地(丘)、固定沙地(丘)和戈壁4種，面積占比達總沙漠化土地面積的98.16%。其具體分布格局為，流動沙地(丘)主要分布于敦煌市西南部庫姆塔格沙漠、敦煌綠洲南部的鳴沙山地區、 金塔縣東南部巴丹吉林沙漠和環綠洲外圍、 阿克塞縣東南緣大哈勒騰河南岸、 肅北縣鹽池灣鎮東南部等。半固定沙地(丘)、固定沙地(丘)多集中于肅北縣和敦煌市境內、包括疏勒河兩岸的天然綠洲、人工綠洲(城鎮、農田等)與流動沙丘或戈壁的過渡區域。戈壁主要分布于南、北兩山至綠洲的緩沖區域，此外在祁連山北側沖積扇內，馬鬃山以北也有分布。風蝕劣地多分布在瓜州縣和阿克塞縣境內，占全市風蝕劣地總面積的99.62%。露沙地、沙漠化耕地、風蝕殘丘和非生物治沙工程地在酒泉市呈零星分布，僅占全市沙漠化土地面積的0.61%。

3.1.2 沙漠化土地類型變化特征

研究表明(表3)，相較2009年的監測數據，酒泉市土地沙漠化總面積減少4.07萬hm2，年變化率-0.08%。其中流動沙地(丘)減少10.48萬hm2，年變化率-1.72%； 半固定沙地(丘)增加6.27萬hm2，年變化率1.27%； 固定沙地(丘)增加2.82萬hm2，年變化率0.60%； 風蝕劣地減少2.22萬hm2，年變化率-3.41%； 戈壁減少3.89萬hm2，年變化率-0.12%； 露沙地增加1.07萬hm2，年變化率68.41%； 沙漠化耕地增加0.48萬hm2，年變化率35.34%； 風蝕殘丘增加1.94萬hm2，年變化率14.19%； 非生物治沙工程地減少0.05萬hm2，年變化率-8.74%。分析可知，酒泉市境內主要沙漠化土地類型中，流動沙地(丘)面積呈遞減趨勢，半固定沙地(丘)和固定沙地(丘)面積趨于增加，而戈壁面積變化輕微，基本處于動態平衡狀態，意味著酒泉市沙漠化土地生境狀況向利于物種生存和定居的方向發展，沙漠化土地面積處于緩慢下降過程。

表3 酒泉市沙漠化土地類型動態變化

沙漠化土地類型沙漠化總面積流動沙地半固定沙地固定沙地風蝕劣地戈壁其他類型面積變化(萬hm2)-407-1048627282-222-389343年變化率(%)-008-172127060-341-0121939

表4 酒泉市土地沙漠化程度現狀

沙漠化程度輕中重極重面積(萬hm2)3378122402049059527占沙漠化總面積比(%)353128021436224

3.2 沙漠化程度現狀與變化特征

3.2.1 沙漠化程度分布現狀

2014年的監測數據發現(表4)，輕度沙漠化土地33.78萬hm2，占沙漠化面積的3.53%； 中度沙漠化土地122.40萬hm2，占沙漠化面積的12.80%； 重度沙漠化土地204.90萬hm2，占沙漠化面積的21.43%； 極重度沙漠化土地595.27萬hm2，占沙漠化面積的62.24%。不同沙漠化發展階段的比例特征說明了酒泉市土地沙漠化程度以重度和極重度為主。

圖2 酒泉市沙漠化土地類型分布現狀 圖3 酒泉市土地沙漠化程度現狀

由圖3分析可得，輕度沙漠化土地主要分布于疏勒河兩岸及人口聚集的綠洲區，區內植被蓋度較高，土壤生態服務功能相對健全，是酒泉市土地生產力最高的區域； 中度沙漠化土地主要位于肅北縣和敦煌市境內，呈斑塊狀分布格局； 重度沙漠化土地主要分布于肅北北山(馬鬃山)以北區域、玉門市西北部和金塔縣西北部，主要的沙漠化類型為半固定沙丘和戈壁； 極重度沙漠化土地主要為流動沙丘(地)和戈壁，多分布于敦煌市西南部庫姆塔格沙漠，敦煌綠洲南部的鳴沙山地區以及南、北兩山至綠洲的緩沖區域內，具有面積大、延伸廣、集中連片的分布特征。

表5 酒泉市土地沙漠化程度動態變化

沙漠化程度輕度中度重度極重度面積變化(萬hm2)8841444-787-1947年變化率(%)626254-075-064

3.2.2 沙漠化程度變化特征

對比2014年和2009年兩期數據可得(表5)，輕度沙漠化土地面積增加8.84萬hm2，年變化率6.26%，中度沙漠化土地面積增加14.44 萬hm2，年變化率2.54%，重度沙漠化土地面積減少7.87萬hm2，年變化率-0.75%，極重度沙漠化土地面積減少19.47萬hm2，年變化率-0.64%。酒泉市沙漠化程度呈現出由極重度和重度向中度與輕度演變的動態規律，由于極重度和重度沙漠化土地的面積基數遠大于輕度和中度沙漠化土地，致使不同沙漠化程度土地間呈現出的轉化比率差異較大。

4 沙漠化土地變化原因分析

研究證實，酒泉市沙漠化土地面積在近5年內(2009～2014年)呈緩慢減少的變化特征，而王小軍等[19]通過影像判讀與現地調查，發現酒泉市沙漠化土地面積2004～2009年間整體遞減。綜合可見，自2004年起，在長達10年的時間梯度內，酒泉市沙漠化演替處于逆轉階段，沙漠化土地面積趨于減少，沙漠化程度逐步降低。

4.1 自然因素

水分是塑造和構建沙漠化土地植被群落結構與分布格局的主導要素[20]，而植被的生長分布狀況是表征沙漠化程度的重要因子[21]，因此，土壤水分含量在很大程度上影響著沙漠化的演替方向。在沙漠化地區，降水和地表徑流是土壤水分的主要補給源[22]，對此，分析近年來(2010～2014年)酒泉市降水量和主要河流徑流量的變化特征(表6)，有助于解讀該區域沙漠化土地面積總體減少，沙漠化程度趨于降低的客觀事實。

分析可得(表6)， 2010～2014年間酒泉市降水呈波動變化特征，整體表現出先增加后下降的規律， 2010年降水量最高，達148.5mm， 2014年最低，僅76.5mm，但5年間有3年降水量處于正距平狀態，且所有年份累計距平均為正值，表明該階段降水量相對較多，很大程度上提高了土壤水分有效性，有利于植被的定居繁衍。

表6 2010～2014年酒泉市降水量變化特征 mm

酒泉市主要的水系包括黑河水系和疏勒河水系，其中黑河的補給來源中，降水占70.2%，屬于典型的降水補給性河流[15]。水文數據顯示，黑河來水自2005年起多以豐水年為主，且李卓侖等[23]通過研究黑河徑流與氣象因子的響應關系，得出黑河徑流在未來短期內仍處于增加趨勢。疏勒河方面，類似于黑河水系， 2007年來徑流量總體處于偏豐狀態[24]，且藍永超等[25]對近50年疏勒河流域山區的氣候變化及其對出山徑流的影響研究顯示，流域近年來有連續降水量增加、徑流量偏豐的趨勢。

綜上發現，近年來(2009～2014年)，得益于酒泉市降水量和主要河流徑流量持續處于偏豐狀態，區域內生態系統普遍面臨的水分脅迫狀態得到緩解，改善了植被特別是對淺層土壤水分有較高依賴性的一年生草本植物的生存環境，一定程度上遏制了沙漠化的加劇與擴張。

4.2 人為因素

(1)生態工程

“十二五”期間，酒泉市陸續實施了重點公益林補償、“三北防護林”四期工程、退耕還林(草)等重點林業生態工程，截止2014年，酒泉市納入中央和省級財政補償的公益林面積達48.20萬hm2。2013年起，通過在敦煌市、金塔縣、玉門市實施沙漠化土地封禁保護補助試點工程，共完成固沙壓沙585.1hm2，減輕了人類活動對脆弱生態系統的不利影響，有效緩解了區域內土地沙漠化。同時在酒泉市境內的敦煌陽關、金塔攔河灣、玉門青山和阿克塞庫姆塔格設立4個國家級沙漠公園，這些公園在提高區域內生物多樣性水平、維持沙漠生態系統功能穩定性方面具有積極作用，同時有助于增強公眾的生態安全意識。此外，疏勒河流域生態移民、敦煌莫高窟崖頂風沙治理以及河西走廊農業綜合開發等工程的實施，亦對沙漠化的防治有所裨益。這些生態工程是人類對酒泉市境內嚴重退化生態系統(沙漠化)的有序合理干預，改善了下墊面的生境狀況，促進了沙生植被的生長、發育和繁衍，提高了生態系統的整體生產力，其直觀變化主要通過植被蓋度的顯著增加來反映，而植被蓋度是表征沙漠化程度的關鍵生態因子。說明生態工程以及政策的實施對沙漠化的遏制作用開始顯現。

(2)節水型社會建設和關井壓田工程實施

解決酒泉市缺水問題，核心是提高用水效率，建設節水型社會[26]。2009年起，酒泉市通過積極推廣以滴灌、管灌、壟膜溝灌等為主的節水灌溉措施，使生產、生活乃至生態用水達到集約化利用的目的。截止2014年，已累計推廣高新節水工程面積1.60萬hm2，累計壓減傳統耕地0.34萬hm2，關閉機井318眼，治理重點風沙口0.23萬hm2，綜合治理沙漠化土地面積1.00萬hm2。節水型社會的建設和關井壓田在一定程度上緩解了沙漠化的進程。

(3)基礎設施建設

人類活動是推動沙漠化加劇的主要因素，但在特定條件下，人類活動能夠在一定程度上遏制沙漠化的加劇。隨著社會經濟的發展，酒泉市加大了基礎設施建設力度，各類公共交通設施，以及以1 000萬kW風電、100萬kW光電為主體的清潔能源基地的大規模建設[27]，使大量沙漠化土地轉變為非沙漠化土地。雖然這些工程建設面積有限，且未從根本上改變區域的氣候條件，但這種土地利用類型的改變，仍然是短期內減少沙漠化土地面積的重要途徑。

5 討論與結論

5.1 討論

研究表明， 2009年以來，酒泉市土地沙漠化呈“雙縮減”的動態特征，即土地沙漠化面積呈減少趨勢，沙漠化程度呈減輕態勢。但應當看到，至2014年，酒泉市土地沙漠化面積仍處于較高水平，達956.34萬hm2，年減少率僅為0.08%，土地沙漠化形勢依舊嚴峻，沙漠化防治仍是今后較長時間內酒泉市生態建設的核心工作。基于遙感技術的酒泉市沙漠化土地現狀與動態分析，在目視解譯過程中存在一定誤差，通過現地驗證與核實發現，這一誤差在統計學允許范圍內，因此，利用遙感技術來判定沙漠化土地現狀和演變趨勢是合理可靠的。該研究有助于酒泉市土地沙漠化時空動態數據庫的建立，為當地土地沙漠化的防治提供科學依據和決策支持，但由于數據的收集和提取過程較為粗糙，導致定性指標較多，定量化分析相對欠缺。對此，在今后研究中，在制定指標和提取數據時，應充分考慮其解釋力度和有效性，以更加精確地量化沙漠化土地的動態規律。

5.2 結論

(1)酒泉市現有沙漠化土地面積956.34萬hm2，沙漠化土地類型主要以流動沙地(丘)、半固定沙地(丘)、固定沙地(丘)和戈壁為主，其面積占到沙漠化土地總面積的98.16%。酒泉市土地沙漠化程度以重度、極重度為主，面積為800.14萬hm2，占沙漠化土地總面積的83.67%。

(2)與2009年監測數據對比，酒泉市沙漠化土地面積減少4.07萬hm2，年變化率-0.08%。沙漠化土地面積呈整體減少的動態變化。重度、極重度沙漠化土地面積年變化率分別為-0.75%和-0.64%，沙漠化程度呈減輕趨勢。

[1] 朱震達.中國北方沙漠化現狀及發展趨勢.中國沙漠， 1985， 5(3)： 3～11

[2] 王濤， 朱震達.我國沙漠化研究的若干問題—Ⅰ.沙漠化的概念及其內涵.中國沙漠， 2003, 23(3): 209～214

[3] 段翰晨， 王濤，薛嫻，等.基于RS與GIS的科爾沁沙地沙漠化時空演變.中國沙漠， 2013, 33(2): 470～477

[4] Millennium Ecosystem Assessment Council.Ecosystems and Human Well-being：Desertification Synthesis.Washington，D C：World Resources Instiute， 2005： 1～22

[5] 王濤. 干旱區綠洲化、荒漠化研究的進展與趨勢.中國沙漠， 2009, 29(1): 1～9

[6] 王根緒， 程國棟，沈永平.近50年來河西走廊區域生態環境變化特征與綜合防治對策.自然資源學報， 2002, 17(1): 78～86

[7] 陳翔舜， 高斌斌，王小軍，等.甘肅省民勤縣土地荒漠化現狀及動態.中國沙漠， 2014, 34(4): 970～974

[8] 魏懷東， 周蘭萍，徐先英，等.2003～2008年甘肅民勤綠洲土地荒漠化動態監測.干旱區研究， 2011, 28(4): 572～579

[9] 何俊林， 吉亞輝，柳建平.甘肅民勤沙漠化防治創新機制研究.開發研究， 2008， (2)： 39～42

[10]齊善忠， 王濤，羅芳，等.黑河中游張掖地區沙漠化土地動態變化.山地學報， 2005, 23(2): 153～157

[11]劉蔚， 王濤，鄭航，等.黑河流域不同類型土地沙漠化驅動力分析.中國沙漠， 2008, 28(4): 634～641

[12]逯軍峰， 董治寶，胡光印，等.甘肅省瑪曲縣土地沙漠化發展及其成因分析.中國沙漠， 2012, 32(3): 604～609

[13]胡夢珺， 潘寧惠，左海玲，等.基于RS和GIS的瑪曲高原土地沙漠化時空演變研究.生態學報， 2017, 37(3): 922～931

[14]張龍生， 李萍，張建旗.甘肅省生態環境脆弱性及其主要影響因素分析.中國農業資源與區劃， 2013, 34(3): 55～59

[15]甘肅省地礦局第三水文地質工程地質隊. 甘肅省北大河流域地下水動態監測報告.2012

[16]胡光印， 董治寶，逯軍峰，等.黃河源區1975～2005年沙漠化時空演變及其成因分析.中國沙漠， 2011, 31(5): 1079～1086

[17]吳薇. 土地沙漠化監測中TM影像的利用.遙感技術與應用， 2001, 16(2): 86～90

[18]甘肅省林業廳. 甘肅省荒漠化和沙漠化監測技術細則.2014

[19]王小軍， 陳翔舜，魏金平，等.甘肅省2004～2009年土地沙漠化時空變化分析.中國沙漠， 2013, 33(1): 33～37

[20]李新榮， 張志山，黃磊，等.我國沙區人工植被系統生態—水文過程和互饋機理研究評述.科學通報， 2013，58(Z1)： 397～410

[21]趙哈林， 周瑞蓮，趙學勇，等.科爾沁沙地沙漠化正、逆過程的地面判別方法.中國沙漠， 2008, 12(1): 8～15

[22]Dunkerley D L.Intra-event intermittency of rainfall：An analysis of the metrics of rain and no-rain periods.Hydrological Processes， 2015， (29)： 3294～3305

[23]李卓侖， 王乃昂，李育，等.近50年來黑河出山徑流對氣候變化的響應.水土保持通報， 2012, 32(2): 7～11

[24]李計生， 胡興林，黃維東，等.河西走廊疏勒河流域出山徑流變化規律及趨勢預測.冰川凍土， 2015, 37(3): 803～810

[25]藍永超， 仵彥卿，康爾泗，等.祁連山北麓出山徑流對氣候變化的響應.蘭州大學學報(自然科學版)， 2001, 37(4): 125～132

[26]羅志楨. 酒泉市水資源現狀及面臨的問題與對策.中國農業資源與區劃， 2007, 28(7): 13～16

[27]李玉忠， 胡秉安.酒泉清潔能源產業園區建設研究.中國農業資源與區劃， 2011, 32(5): 16～21