八步沙林場防沙治沙區植被覆蓋度時空演變分析

胡曉娟，李 霞，王 訥

（1.甘肅省基礎地理信息中心，甘肅 蘭州 730000；2.八步沙林場，甘肅 古浪 733102）

植被覆蓋度是衡量地表植被狀況的重要指標，也是影響土壤侵蝕與水土流失的主要因子，對于區域環境變化和監測研究具有重要意義。隨著遙感技術的發展，植被覆蓋度的測量方法已由傳統的地面測量發展為遙感估算，這也為大面積植被覆蓋度的監測提供了可能。

位于騰格里沙漠南端的古浪縣是我國荒漠化重點監測縣，風沙線長達132 km，重點風口達20多個，沙害面積大，沙丘類型復雜，長期危害著周邊鄉鎮、耕地與道路。“八步沙六老漢”祖孫三代自1981年開始防沙治沙，使得區域植被覆蓋度得以明顯提升。為定量評估八步沙林場三代人治沙成效，本文以八步沙林場防沙治沙區域及其周邊鄉鎮為研究區，利用植被覆蓋度遙感估算方法探討了區域內近28 a來植被覆蓋變化的時空特征。

1 研究區概況

研究區位于甘肅省古浪縣北部，騰格里沙漠南緣，地理位置為103°28′E、37°36′N；海拔為1 620 m；年均氣溫為6.6℃，極端最高溫為37.2℃，極端最低溫為-26.4℃；年均降水量為207.7 mm，年均蒸發量為2 807.4 mm；最大風速為130 m/s，平均風速為3.8 m/s，年大風時間為182 d，年沙塵暴時間為30.5 d，多發生在冬、春季；無霜期為153 d；屬風力吹蝕堆積而成的風沙地貌，包括流動沙丘、固定半固定沙丘、平緩覆沙地等地貌類型。

2 數據來源與處理

2.1 數據來源

結合古浪縣八步沙林場防沙治沙的歷史階段與研究區植被生長旺盛期（6～9月），本文采用的數據為1991-06-25、1999-08-10、2005-07-25、2010-09-09、2015-08-14、2019-09-12等6景Landsat衛星遙感影像，軌道號為131/034，數據來源于地理空間數據云平臺（http://www.gscloud.cn）。

2.2 數據處理

1）影像預處理。以ENVI5.3為遙感圖像處理平臺，對下載的6景原始影像進行輻射定標、大氣校正等處理，并根據研究區邊界矢量文件對目標區域進行裁剪。由于下載的影像已進行初步的幾何糾正，具有較好的精度，因此本文無需對影像作幾何校正。

2）歸一化植被指數（NDVI）提取。NDVI是基于植物對紅光波段和近紅外波段的光譜通道來組合設計的，在提取植被信息方面具有較高的靈敏度，且能消除大部分的噪音、輻射、陰影等帶來的誤差。NDVI的計算公式為；

式中，Red為紅光波段反射值，對應Landsat4/5 TM影像第三通道，對應Landsat8 OLI影像第四通道；NIR為近紅外波段反射值，對應Landsat4/5 TM影像第四通道，對應Landsat8 OLI影像第五通道。

利用ENVI5.3中的NDVI工具計算各景影像的NDVI值，作為植被覆蓋度計算的基礎數據。

3）植被覆蓋度計算。植被覆蓋度是指植被（包括葉、莖、枝）在地面的垂直投影面積占統計區總面積的百分比。本文采用李苗苗[1]等在像元二分模型基礎上研究的模型，利用NDVI近似估算植被覆蓋度，即

式中，NDVIsoil為完全是裸土或無植被覆蓋區域的NDVI值；NDVIveg為完全被植被覆蓋區域的NDVI值，即純植被的NDVI值。

在沒有實測數據的情況下，在ENVI軟件中根據NDVI累計頻率表，提取累計頻率為2%的NDVI值作為NDVIsoil，累計頻率為99.5%的NDVI值作為NDVIveg。

在參考相關文獻的基礎上，結合研究區生態環境特點，本文將植被覆蓋度等級劃分為極高覆蓋度（VFC＞0.7）、高覆蓋度（0.5＜VFC≤0.7）、中覆蓋度（0.3＜VFC≤0.5）、低覆蓋度（0.1＜VFC≤0.3）和極低覆蓋度（VFC≤0.1）5個等級。

3 研究結果分析

3.1 植被覆蓋度總體變化趨勢分析

本文通過計算得到1991年、1999年、2005年、2010年、2015年、2019年6期植被覆蓋度等級分布數據，如圖1所示，可以看出，1991-2019年區域內植被覆蓋度明顯提升，主要表現在兩個方面；①以裸露沙地為主的極低植被覆蓋度區域面積逐年減少，低植被覆蓋度及以上等級所占面積逐年增加，區域植被覆蓋由綠洲與沙漠接壤處向沙漠腹地延伸明顯；②植被覆蓋度等級提升明顯，主要集中在2015-2019年，中植被覆蓋度及以上等級所占面積顯著增加。

圖1 八步沙林場防沙治沙區域及其周邊鄉鎮植被覆蓋度分級圖

3.2 研究區植被覆蓋度時空演變分析

本文對相鄰年度植被覆蓋度數據進行矩陣運算，得到1991-1999年、1999-2005年、2005-2010年、2010-2015年、2015-2019年區域內各等級植被覆蓋度面積轉移矩陣，如表1～5所示。

表1 1991－1999年不同等級植被覆蓋度面積轉移矩陣/km2

由統計結果可知，1991-1999年區域內植被覆蓋度變化主要為低植被覆蓋度向極低植被覆蓋度轉換，轉換量為355.82 km2；極低植被覆蓋度向低植被覆蓋度的轉換量為162.93 km2；中植被覆蓋度向低植被覆蓋度的轉換量為143.82 km2；其他類型轉換面積相對較小，該時段內區域植被覆蓋度降低，主要原因為惡劣天氣影響，尤其是1993年的特大沙暴對區域內植被，尤其沙漠邊緣植被產生了毀滅性影響。1999-2005年區域內植被覆蓋度變化主要為極低植被覆蓋度向低植被覆蓋度的轉換，轉換量為309.02 km2，同時低植被覆蓋度向極低植被覆蓋度的轉換量為149.61 km2，極低植被覆蓋度的轉出量大于轉入量，因此極低植被覆蓋度區域面積呈減少趨勢；中植被覆蓋度總轉入量為117.62 km2，主要來源為低植被覆蓋度；其他類型之間轉換量較小。2005-2010年區域內植被覆蓋度變化主要為極低植被覆蓋度向低植被覆蓋度的轉換，轉換量為545.65 km2，同時極低植被覆蓋度的總轉入量為76.66 km2，遠小于轉出量，因此呈明顯減少趨勢；低植被覆蓋度向中植被覆蓋度的轉換量為69.07 km2，同時中植被覆蓋度向低植被覆蓋度的轉換量為55.79 km2，中植被覆蓋度轉入量大于轉出量；其他類型之間的轉換量較小。2010-2015年區域內植被覆蓋度變化仍以極低植被覆蓋度向低植被覆蓋度轉換為主，轉換量為883.84 km2，低植被覆蓋度區域面積呈顯著增加趨勢；低植被覆蓋度向中植被覆蓋度及以上等級的轉換量，以及中植被覆蓋度向高植被覆蓋度及以上等級的轉換量也開始明顯增加。2015-2019年區域內植被覆蓋度主要為低植被覆蓋度向中植被覆蓋度及以上等級的轉換，總轉換量為1 078.79 km2；其中中植被覆蓋度總轉入量最大，為791.98 km2，主要來源于低植被覆蓋度；高植被覆蓋度總轉入量為396.54 km2，主要來源于低植被覆蓋度和中植被覆蓋度；低植被覆蓋度總轉入量為298.36 km2，主要來源于極低植被覆蓋度；極高植被覆蓋度的轉入量也明顯增加，為123.05 km2，主要來源于中植被覆蓋度、低植被覆蓋度和高植被覆蓋度，該時段內植被覆蓋度轉換比其他時段劇烈，植被覆蓋度總體明顯提升。

表2 1999－2005年不同等級植被覆蓋度面積轉移矩陣/km2

表3 2005－2010年不同等級植被覆蓋度面積轉移矩陣/km2

表4 2010－2015年不同等級植被覆蓋度面積轉移矩陣/km2

表5 2015－2019年不同等級植被覆蓋度面積轉移矩陣/km2

由此可見，區域內植被覆蓋度時空變化包括2005年和2015年兩個轉折點，2005年之前，區域內植被覆蓋度轉換相對緩慢，表現為極低植被覆蓋度與低植被覆蓋度的相互轉換、沙漠與低植被覆蓋度植被的相互轉變，區域內植被覆蓋度規模小且不穩定，易受沙漠侵害；2005年以后，極低植被覆蓋度區域開始明顯減少，表現為植被面積快速增加、沙漠面積急速減少，區域內植被規模逐漸穩定且快速增加；2015之前，區域內植被覆蓋度主要表現為極低植被覆蓋度向低植被覆蓋度的轉換，即沙漠向植被轉換；2015年以后，區域內植被覆蓋度變化主要表現為植被覆蓋度等級間的梯度遞增。

3.3 研究區植被覆蓋度變化影響因素分析

1）自然因素。氣溫與降水是影響區域植被覆蓋度變化的主要自然因素，尤其是降水量對沙漠區域天然植被的生長具有決定性影響。由收集的古浪縣1991-2019年年平均氣溫（圖2）和年累積降水量資料（圖3）可知，受全球氣候變暖影響，監測區氣溫和降水量均呈波動上升趨勢，2010年以后，降水量波動趨勢減緩，整體處于研究時段的高水平，并在2019年達到最高值，這與研究區植被覆蓋度時空演變的趨勢一致，因此氣溫與降水的增加是植被覆蓋度等級提升的重要影響因素。

圖2 研究區年均氣溫趨勢圖

圖3 研究區降水量趨勢圖

2）人為因素。獨特的地理位置和地貌特征決定了古浪縣的生態屏障地位。多年來，古浪縣始終堅持把防沙治沙作為一項長期性、全局性的戰略任務來抓，努力在風沙線上構建百里生態屏障；將八步沙林場探索形成的因地制宜、科學防沙治沙經驗在全縣范圍內進行推廣實施，并實施了國家重點生態功能區轉移支付、沙化土地封禁保護區建設、省級防沙治沙、三北防護林等多項國家重點生態建設工程；將防沙治沙范圍從城區機關擴大到全縣各鄉鎮和社會團體，在沙區大規模開展義務壓沙。為推動和鼓勵全民參與治沙，將境內需要治理的宜林沙荒地以及治理后形成的灌木林地、未成林造林地實行承包經營。由植被覆蓋度變化的區域特征可知，在重點防沙治沙、封育管護的區域植被覆蓋度變化相對較顯著。由此可見，大規模實施的防沙治沙和生態環境保護措施以及社會公眾的積極響應與參與，是研究區植被覆蓋面積持續擴展、生態環境改善的重要因素。

4 結 語

1）1991-2019 年研究區植被覆蓋面積顯著提升。極低植被覆蓋度區域面積呈減少趨勢，低植被覆蓋度及以上等級區域面積呈增加趨勢。

2）從植被覆蓋度的時空演變特征來看，1991-2005年區域內植被覆蓋度變化主要為極低植被覆蓋度與低植被覆蓋度之間的相互轉換；2005-2015年區域內植被覆蓋度變化主要為極低植被覆蓋度向低植被覆蓋度的轉換；2015-2019年區域內植被覆蓋度變化主要為低植被覆蓋度向中植被覆蓋度及以上等級的轉換，植被覆蓋度的轉換量較其他時段更劇烈。

3）從植被覆蓋度變化的區域來看，2015年之前，研究區內植被覆蓋度變化以面積擴展為主，由綠洲與沙漠接壤處向騰格里沙漠腹地持續延伸；2015年以后植被覆蓋度以等級提升為主，在重點防沙治沙和封育管護區域，植被覆蓋度等級提升較顯著。

4）研究區植被覆蓋度的改善受自然和人為雙重因素影響。氣溫、降水量的增加以及大規模防沙治沙與生態環境保護政策的實施、社會公眾的積極參與，是研究區植被覆蓋增加、生態環境改善的重要因素。