2000-2013年西安市植被覆蓋度時空演變

賈 路, 任宗萍, 李占斌,2, 徐國策, 時 鵬, 張譯心, 王 斌

(1.西安理工大學 省部共建西北旱區生態水利國家重點實驗室, 西安 710048;2.中國科學院 水利部 水土保持研究所 黃土高原土壤侵蝕與旱地農業國家重點實驗室, 陜西 楊凌 712100)

當前中國區域水土流失[1]與荒漠化、干旱等問題依然嚴重，急需進行生態環境治理。隨著工業化和城鎮化的迅速發展，廣大人民對我國生態環境美好的生活越來越向往[2]。植被在地表生態系統的能量交換、水循環、碳循環、生物地球化學循環和維持中發揮著重要作用，區域植被恢復情況的好壞，對區域生態環境變化具有重要意義。

在人類活動和氣候變化高強度共同作用下，區域水文過程發生了非常明顯的變化[3]。氣候變化、人類活動以及水文過程的變化直接或者間接對植被恢復等生態過程產生重要的影響[4]。標準化植被指數(NDVI)數據具有較高空間覆蓋和長時間序列等優點，已成為分析植被生長歷史、監測當前狀況和預測未來變化的最重要工具，通常被用于大尺度植被變化和空間密度分布的研究[5-7]。目前我國基于NDVI數據對很多地區的植被變化及其影響因素都進行分析。有研究表明，中國北方地區年NDVI呈輕微上升趨勢[8]。人類活動對植被恢復起著一定的積極的作用，在青海、西藏高原地區有近61.2%的草原地區實施了一定的工程措施進行植被恢復，其中氣候變化對56.7%的退化草原存在影響[9]。盡管以往的研究對不同地區植被覆蓋的變化和時空分布模式提供了十分有用的信息，但許多研究主要集中在植被覆蓋的時空分布特征的驅動因子和植被動態變化上，對于植被恢復的空間差異變化研究卻極為少見。因此，全面系統地對西安市植被覆蓋度時空變化及其差異進行分析研究是一個十分有意義的問題。

本文采用Mann-Kendall檢驗法對西安市2000—2013年植被覆蓋度進行趨勢分析，采用Hurst指數對西安市2000—2013年植被覆蓋度進行持續性分析，同時通過Pettitt檢驗法對西安市2000—2013年植被覆蓋度進行突變點分析，并結合相對發展率和重心轉移模型對2000—2013年西安市植被覆蓋度進行區域差異分析。研究成果有助于進一步深化對西安市植被覆蓋度變化及其影響因素的認識，有利于為西安市植被恢復等生態建設工程提供一定的科學依據。

1 研究區概況與方法

1.1 研究區概況

西安市地處陜西關中平原，位于107°—109°E，33°—34°N，共統轄9區4縣，總面積約9 983 km2，市區面積1 276 km2。地質狀況復雜多樣，地質、構造和巖相具有不同的發育特點，地形南高北低(圖1)，秦嶺位于西安市南部，渭河位于西安市北部。地貌類型主要為平原和基巖山地兩大類，平原所占比例為46%，基巖山地所占比例為54%。氣候類型為暖溫帶大陸性季風氣候，四季分明，雨熱同期，年內最高溫度在7月份，年內最低溫度在1月份，年平均降水為600 mm左右，由南向北逐漸遞減。秦嶺山地基本屬自然植被，渭河平原、驪東南丘陵與黃土臺塬屬栽培植被，兩大植被區域的分布與地貌區域范圍大體一致。境內的自然植被分布于秦嶺山區，隨海拔高度變化演替，依次出現落葉闊葉林、針闊葉混交林、針葉林、高山灌叢與草甸等植被類型。

圖1 西安市地理位置

1.2 數據來源

本文所用NDVI數據是一個長期的MODIS植被指數數據集。數據集由中國科學院計算機網絡信息中心國際科技數據鏡像網站提供(http:∥www.gscloud.cn)，數據時間為2000—2013年，空間分辨率為500 m，時間分辨率為旬，用最大值合成法(MVC)獲得NDVI年值，植被覆蓋率是根據以下公式計算的：

f=(NDVI-NDVImin)/(NDVImax-NDVImin)

(1)

式中：f為植被覆蓋度；NDVI為像素的植被指數；NDVImax，NDVImin為研究區NDVI的最大值和最小值。

降水數據來源于中國氣象數據服務中心中國地面氣候資料的年度數據集，包含西安市及其附近區域耀縣、漢中、武功、華山、商州、鎮安、佛坪7個氣象站點2000—2013年逐年降水量，通過算數平均法對7個雨量站年降雨量取平均值,將點雨量轉換為區域面降雨量。

1.3 研究方法

Mann-Kendall檢驗(M-K)是判斷數據序列變化趨勢的有效方法，被廣泛應用在氣象和水文序列的分析中[10]，本文采用Mann-Kendall檢驗進行年植被覆蓋度的變化趨勢分析。

Hurst指數(H)在定量表征時間序列的持續性或長期相關性具有廣泛的應用性，已被廣泛用于氣候和水文序列的分析中[11]。其中H代表Hurst指數。H的值在[0,1]的范圍內。當H=0.5時，表明植被覆蓋的時間序列是隨機序列，是不可持續的。當H>0.5時，植被覆蓋度的變化與目前的趨勢基本一致，表明植被的可持續性是正向的。H<0.5表示負可持續性，未來植被覆蓋變化將與當前趨勢相反。

Pettitt檢驗法[12]采用Mann-Whitney中Ut,N值檢驗同一總體中兩個樣本X1,…,Xt和Xt+1,…,XN，Pettitt檢驗的零假設為沒有變化點，當|Ut,N|取最大值時對應的Xt被認為是可能的突變點。當p≤0.05時認為數據中存在均值變異點。其顯著性水平可由下式計算：

(2)

相對發展率[13](NICH)可以用來分析研究區內，不同像元在研究時段初期與末期植被覆蓋度變化量與研究區整體研究時段初期與末期變化量的差異，計算公式如下：

(3)

式中：Y2i為2013年第i個柵格植被覆蓋度；Y1i為2000年第i個柵格植被覆蓋度；Y2為2013年西安市植被覆蓋度平均值；Y1為2000年西安市植被覆蓋度平均值。

空間重心轉移模型[14]對地理對象的空間分布進行描述的重要方法，可以用來研究植被覆蓋度空間重心隨時間的變化，本研究采用空間重心轉移模型來研究西安市2000—2013年植被覆蓋度的重心變化，計算公式如下：

(4)

(5)

式中：Xm和Ym為植被覆蓋度空間分布重心的經緯度；cmi為第i個柵格的植被覆蓋度值；xi為第i個柵格的經度；yi為第i個柵格的緯度。

本研究中所有數據計算通過ArcGIS 10.2軟件及R語言編程實現，使用Excel軟件進行基礎數據整理。

2 結果與分析

2.1 植被覆蓋度年際變化特征

2000—2013年西安市植被覆蓋度最小值空間變化范圍為0.33～0.97，最大值空間變化范圍為0.56～1.00，多年平均值空間變化范圍為0.44～0.98，標準差空間變化范圍為0.01～0.17，變異系數空間變化范圍為0.01～0.29。根據變異系數Cv的分類系統：弱變異Cv≤10%，中等變異10%

從圖2中可以看出，2000-2013年西安市植被覆蓋度變異程度中，中等程度變異面積占總面積10.02%，弱變異面積占總面積89.98%，總體來說，2000-2013年西安市植被覆蓋變異強度較弱，中等程度變異范圍較小，中等程度變異主要分布于西安市城區及臨潼區南部和藍田縣北部地區等，該區人類活動強烈，對植被變化影響較大。

圖2 2000-2013年西安市植被覆蓋度變異等級空間分布

2.2 植被覆蓋度的趨勢性與持續性

為了進一步分析2000—2013年西安市植被覆蓋度變化的趨勢特征和未來變化的持續性，結合Mann-Kendall檢驗和Hurst指數計算得到了2000—2013年西安市植被覆蓋度M-K趨勢檢驗圖和2000—2013年西安市植被覆蓋度Hurst指數分布圖。

如圖3所示，區域植被覆蓋度整體呈增加趨勢，其中呈增加趨勢的面積占區域總面積86.54%，顯著增加的面積占區域總面積的37.25%，顯著減少的面積占區域總面積的1.68%，顯著減少面積主要分布城區周邊，該區域人類活動強烈，對植被破壞較大。

圖3 2000-2013年西安市植被覆蓋度M-K趨勢檢驗

由圖4可知，區域植被覆蓋度Hurst指數空間變化范圍為0.18～0.72，其中正向持續性的面積占區域總面積72.62%，負向持續性的面積占區域總面積27.38%。西安市35.30%面積植被恢復呈持續改善(H>0.5，且M-K統計量>0)。

圖4 2000-2013年西安市植被覆蓋度Hurst指數分布

2.3 植被覆蓋度突變點空間分布識別及突變年份分析

Pettit檢驗是一種常用地時間序列突變點的方法，可以簡單快速有效識別一個時間系列發生突變的時間，因此在突變點檢驗應用較為廣泛。為此，本研究對所研究時期內的植被覆蓋度進行Pettit檢驗，對西安市植被覆蓋度突變性進行分析，得到各西安市植被覆蓋度突變年份的時空分布情況。

如圖5所示，2000—2013年西安市植被覆蓋度突變年份均顯著發生在2004年、2005年、2006年、2007年、2008年，發生顯著突變年份的面積占總面積的17.58%，其中發生顯著突變年份在2004年、2005年、2006年、2007年、2008年的面積占總面積的比例依次是1.71%，4.85%，7.10%，2.86%，1.06%。本文將突變點前植被覆蓋度M-K檢驗統計量小于零且突變點后植被覆蓋度M-K檢驗統計量大于零的突變點定義為凹突變點，凹突變點表明外界環境對植被覆蓋度變化具有促進作用；將突變點前植被覆蓋度M-K檢驗統計量大于零且突變點后植被覆蓋度M-K檢驗統計量小于零的突變點定義為凸突變點，凸突變點表明外界環境對植被覆蓋度變化具有抑制作用。根據分析結果，2000—2013年西安市植被覆蓋度凹突變點占總突變點的20.67%，2000—2013年西安市植被覆蓋度凸突變點占總突變點的4.18%，因此，2000—2013年西安市植被恢復在2004年、2005年、2006年、2007年、2008年期間受到了外界環境的一定促進增長作用，發生突變的地點主要集中于灞橋區、藍田縣等人類活動較少的地區，該區域植被覆蓋度較高，植被覆蓋良好。

圖5 2000-2013年西安市植被覆蓋度突變點空間分布

2.4 植被覆蓋度變化的空間差異性及重心轉移

為了進一步分析在外界環境干擾下西安市植被覆蓋度變化的空間差異性，計算2000—2013年西安市植被覆蓋度相對發展率，得到2000—2013年西安市植被覆蓋度相對發展率空間分布圖。

對圖6進行分析得到，2000—2013年西安市植被覆蓋度相對發展率空間變化范圍為-9.07～7.49，相對發展率的負值占區域總面積的20.77%，相對發展率的負值主要出現在城區，說明城區植被覆蓋度變化與區域整體變化趨勢相反，人類活動的范圍和強度在不斷增大。相對發展率的最大值分布于灞橋區附近，說明該區域區植被覆蓋度增長速率優于區域整體植被恢復程度。

圖6 2000-2013年西安市植被覆蓋度相對發展率空間分布

由圖7可以看出，2000—2013年西安市植被覆蓋度空間重心呈現由西南逐漸向東北方向轉移的趨勢。2000—2001年植被覆蓋度重心從東北向西南轉移，2001—2002年植被覆蓋度重心從西南向東北轉移且轉移距離比較長，2002—2012年植被覆蓋度重心變化趨勢不明顯，2012—2013年植被覆蓋度重心從西南向東北轉移。2013年西安市各區縣人口數量從大到小排序依次為：長安區>雁塔區>臨潼區>碑林區>周至縣>藍田縣>蓮湖區>鄠邑區>未央區>灞橋區>新城區>高陵區>閻良區。整體來說，2000—2013年西安市植被覆蓋度增加趨勢逐漸從西南部的周至等人口密度較小的地區向臨潼區等人口密度較大的東北部區域轉移，說明隨著社會的進步，人們對植被等生態環境的保護更加注重。

圖7 2000-2013年西安市植被覆蓋度重心轉移

3 討 論

3.1 氣候因素及對植被恢復的影響

降雨作為影響區域植被變化的主要氣候因素[15]，降水變化對區域徑流變化和土壤水的變化具有直接的影響，從而影響到植被對水分的吸收，進而影響到植被的生長變化。近些年來全球氣候發生顯著變化，氣溫升高，降水增多，對區域下墊面變化具有巨大影響[16]。

從圖8可以看出，2000—2013年西安市面降雨量最小值為570.18 mm，最大值為1 002.74 mm，多年平均值為734.00 mm，變異系數為0.17，變異強度為中等程度。2000—2013年西安市植被覆蓋度均值變化的最小值為0.83，最大值為0.90，多年平均值為0.87，變異系數為0.02，變異強度為弱變異。2000—2013年西安市面降雨量和區域植被覆蓋度均值均呈現增加趨勢，根據M-K趨勢檢驗結果，降雨量變化呈不顯著增加趨勢，植被覆蓋度均值呈顯著增加趨勢。降雨量的Hurst指數為0.51，植被覆蓋度均值的Hurst指數為0.66，降雨量與植被覆蓋度均值未來呈現增加趨勢，降雨量與植被覆蓋度均值皮爾遜相關性系數為0.47(p<0.09)，呈不顯著正相關。根據轉移重心模型，2000年西安市降雨量空間分布重心坐標為108.72°E，33.83°N，2013年西安市降雨量空間分布重心坐標為108.54°E，33.79°N，2000—2013年西安市降雨量空間分布重心呈現從東北向西南轉移的趨勢，與植被變化呈現相反的趨勢。

3.2 不同土地利用變化對植被覆蓋的影響

由圖9可知，2000—2013年西安市各類土地利用植被覆蓋度均有增加趨勢，2000年西安市耕地、林地、草地、建設用地、未利用地的植被覆蓋度均值依次為0.80，0.91，0.86，0.76，0.83，其中林地植被覆蓋度最大，建設用地植被覆蓋度最小。2013年西安市耕地、林地、草地、建設用地、未利用地的植被覆蓋度均值依次為0.87，0.94，0.93，0.76，0.87，其中林地植被覆蓋度最大，與2000年相比草地與耕地的植被覆蓋度增加最多，建設用地植被覆蓋度最小，與2000年相比基本未發生變化。

圖8 2000-2013年西安市植被覆蓋度與降雨時間序列

圖9 2000-2013年西安市各類型土地利用植被覆蓋度變化

以上研究說明，2000—2013年西安市降雨對植被變化存在一定程度促進作用，但主要是土地利用等人類活動引起的植被覆蓋度變化。

4 結 論

(1) 2000—2013年西安市植被覆蓋度變異程度整體較弱，中等程度變異面積占總面積10.02%，弱變異面積占總面積89.98%，主要分布于西安市城區及臨潼區南部和藍田縣北部地區等。

(2) 2000—2013年西安市植被覆蓋度整體呈增加趨勢，其中呈增加趨勢的面積占區域總面積86.54%，正向持續性的面積占區域總面積72.62%，負向持續性的面積占區域總面積27.38%，35.30%面積植被恢復呈持續改善。

(3) 2000—2013年西安市植被覆蓋度突變年份均顯著發生在2004年、2005年、2006年、2007年、2008年，發生顯著突變年份的面積占總面積的17.58%。

(4) 2000—2013年西安市植被覆蓋度相對發展率空間變化范圍為-9.07～7.49，相對發展率的負值占區域總面積的20.77%，城區植被覆蓋度變化與區域整體變化趨勢相反，人類活動的范圍和強度在不斷增大，西安市植被覆蓋度空間重心呈現由西南逐漸向東北方向轉移的趨勢。

(5) 2000—2013年西安市面降雨量和區域植被覆蓋度均值均呈現增加趨勢，降雨量與植被覆蓋度均值皮爾遜相關性系數為0.47(p<0.09)，呈不顯著正相關。降雨量空間分布重心呈現從東北向西南轉移的趨勢，與植被變化呈現相反的趨勢，西安市降雨對植被變化存在一定程度促進作用，但主要起作用的是非降雨因素。