汾河流域植被覆蓋度時空特征與地形因子的關系

張會霞, 李永梅, 張 娜

(太原師范學院 地理科學學院, 山西 榆次 030619)

近幾年，自然環境問題在經濟發展中日益突出，對人民生活造成了很大的影響[1]。中國不斷加大生態文明建設，力求建設更加文明和諧更加美麗的中國。植被作為陸地生態系統最為重要的主體，不僅可以保留土壤中的水分，還可以防沙固土，減少水土流失，一定程度上反映出全球生態環境變化狀況[2]。植被覆蓋變化已是全球環境變化和陸表生態系統相關研究的熱門問題之一[3-4]，對研究植被變化情況及其程度則需引用植被覆蓋度(FVC)來進行反演，植被覆蓋度可將地表植被在遙感影像上量化，便于更好地進行研究分析[5]。

目前，植被覆蓋度提取方法主要以地表實測和遙感監測為主[5-6]，但多采用遙感監測。面對研究范圍較大的地區，在時間序列較長的情況下多數學者采用了高分辨率的MODIS數據產品[6]，以便得到更加精準的數據，如孫天瑤等[7]選取MODIS NDVI數據對塔里木河流域2000—2018年的植被覆蓋度時空格局進行研究，同時探討了氣候和人類活動對植被變化的影響；齊敬輝等[8]基于MODIS數據，采用像元二分模型、趨勢分析等方法對2000—2014年疏勒河流域植被覆蓋時空變化進行研究，得出疏勒河流域植被覆蓋度呈明顯改善的狀態；李曉蕾等[9]同樣基于MODIS遙感數據，利用趨勢分析、偏相關分析等方法對葉爾羌河流域植被覆蓋度時空變化進行研究，得出植被覆蓋度在16 a間呈現微弱的增長趨勢。

汾河作為黃河的主要支流，其生態環境保護狀況和沿線城市經濟構成黃河流域生態保護和高質量發展的重要組成部分和影響因素[10]。汾河流域地處黃土高原，土壤疏松、植被稀疏，在夏季多雨季節很容易造成水土流失，是嚴重缺水的區域。植被是地表生態系統中最為重要的組成部分[11]，是生態環境進行監測的基礎，因此對流域內的植被覆蓋度進行時空變化的研究具有一定的理論價值，結合遙感技術對植被進行反演，可為汾河流域的生態環境保護和經濟發展提供一定程度上的理論科學幫助。

本研究以汾河流域為研究區，選取2005—2020年中的5 a(2005，2009，2013，2017，2020年)MODIS NDVI數據，利用ArcGIS 10.2，ENVI 5.3等軟件對NDVI數據進行處理，對汾河流域的植被覆蓋度進行模型估算，從而研究其時空特征，同時探討不同地形因子與汾河流域植被覆蓋度關系。

1 研究區概況

汾河處于黃河的中游地段，其流域坐落于山西省內，地理位置35°13′04″—39°04′04″N，110°26′42″—113°26′56″E。汾河全長約為715 km，主要流經運城、臨汾、呂梁、晉中、太原、忻州等地級市中的43個縣市(區)，萬榮縣作為最終的流經地，在此處匯入黃河。汾河流域面積達到39 740 km2，超過山西全省總面積的四分之一。汾河，是山西最大的河流，汾河流域則是山西人口密集、工農業集中的重要區域，是山西省重點的發展保護地區。汾河流域地理位置見圖1。

圖1 汾河流域在山西省的位置

2 研究數據與方法

2.1 數據來源及處理

2.1.1 數據來源 本研究采用的MODIS NDVI遙感數據是從NASA官網(http:∥ecocast. arc. nasa. gov)下載的MOD13Q1數據產品，該產品的空間分辨率為250 m，時間分辨率為16 d，研究區域處在的行列號為h26v05和h27v05。其中選取了汾河流域2005—2020年中的5 a數據(2005，2009，2013，2017，2020年)來研究，以7—10月作為研究時段，此時段植物生長處于旺季，可更好的對研究區的植被覆蓋度進行時空特征的分析。DEM數據(數字高程模型)則是從地理空間數據云官網(http:∥www. gscloud. cn)下載所得，所用產品為90 m空間分辨率的SRTM DEM數據。

2.1.2 數據處理 對于MODIS遙感影像數據，所獲取的數據格式為HDF，因汾河流域處在兩個行列號間(h26v05和h27v05)，則需利用MRT(MODIS Reprojection Tool)軟件進行拼接和格式轉化，地理坐標統一轉化為WGS-1984坐標系。在ENVI 5.3中利用波段運算去除NDVI中的無效值。在ArcGIS 10.2軟件中利用ArcMap根據汾河流域的矢量數據進行掩膜提取，從而得到汾河流域研究時段的所有影像數據。為了方便對每年的植被覆蓋度進行研究，本研究利用最大值合成法(MVC)將每年研究時段(7—10月)的MODIS NDVI數據合成最大NDVI數據集。

將DEM數據在ArcGIS 10.2軟件中進行地理坐標轉化為WGS-1984坐標系，再根據研究區進行掩膜提取，從而得到研究區的DEM數據，并在ArcGIS 10.2軟件中分別生成坡度圖、坡向圖、高程圖。

2.2 研究方法

2.2.1 像元二分模型 像元二分模型是一種利用遙感數據來進行估算的模型，它不僅簡單，而且很實用，它通常被用來估算某地區的植被覆蓋度[12]。植被覆蓋度(FVC)通常需要利用NDVI(歸一化植被指數)來輔助計算，利用植被在紅光波段具有高吸收，而在近紅外波段具有高反射的特點，采用改進的像元二分模型來進行估算植被覆蓋度，其公式為[6]：

式中：FVC表示植被覆蓋度； NDVIsoil表示裸土或沒有植被覆蓋地區植被指數； NDVIveg表示完全被植被覆蓋地區植被指數。本研究選取置信度為0.5%附近累計百分比作為NDVIsoil，同樣選取置信度在99.5%附近累計百分比作為NDVIveg，根據以上估算模型在ENVI 5.3中估算汾河流域每年植被覆蓋度。參考韓麗君等[6]對汾河流域研究將植被覆蓋度進行等間距分級,分級標準見表1。

表1 汾河流域植被覆蓋度劃分

2.2.2 線性趨勢分析 為了更好地對汾河流域2005—2020年的植被覆蓋度進行時空變化趨勢研究，植被覆蓋度隨時間變化研究采取了一元線性回歸分析法，基于最小二乘法來擬合植被覆蓋度的變化斜率[13]。具體公式為：

式中：θslope表示研究時段植被覆蓋度線性擬合變化趨勢斜率；n為總的研究年數(n=5)； FVCi表示第i年的植被覆蓋度(這里取植被覆蓋度均值)。其中θslope為正則表示植被覆蓋度改善,為負則代表植被覆蓋度減弱,為零則表示植被覆蓋度無變化。

3 結果與分析

3.1 汾河流域植被覆蓋度時空變化特征

3.1.1 植被覆蓋隨時間變化的特征 通過線性趨勢分析，再借助于ArcGIS 10.2軟件中的面積制表可得到汾河流域上、中、下游的植被覆蓋度變化程度。從整體上(表2)可以看出，2005—2020年汾河流域植被覆蓋度的變化程度主要為增加趨勢，其占據了研究區植被覆蓋總面積的70.19%；而植被覆蓋度不變的區域只有0.39%。其中，植被增加區域主要集中在汾河流域中下游地區，其在中游地區尤為顯著，面積比例為31.12%。同時，植被覆蓋度下降區域以中游地區為主，占據總面積的19.82%，植被覆蓋度不變的區域占比極少，但還是可以看出，汾河中游區域的植被變化程度不變的占比最多，但總體上來說，汾河流域的植被覆蓋度整體呈現改善的趨勢。

表2 2005-2020年汾河流域上、中、下游趨勢分析結果

利用像元二分模型，對汾河流域植被覆蓋度2005—2020年整合的5期遙感數據進行估算統計分析，以此得到每年不同植被覆蓋度的分布類型的面積和比例情況(圖2)。從圖2中可以看出，從2005—2020年整體上看，研究區的植被覆蓋呈現改善狀態，植被覆蓋度主要由低、中低、中向中高、高變化，其中：中、中高植被覆蓋度變化最為顯著，中覆蓋度區域面積減少了10.93%，中高覆蓋度區域面積增加了15.56%。

圖2 2005-2020年汾河流域不同植被覆蓋度區域面積比例

據圖2可知，2005—2009年汾河流域的植被覆蓋度變化主要體現在中高覆蓋度上，即從2005年的32.31%增加到2009年的38.18%，增加了5.87%，低、中低植被覆蓋度均下降。

2013年的植被覆蓋度顯著增高，從中可以看出中高、高植被覆蓋度占據汾河流域總植被覆蓋區的79.35%，其增加的主要原因是于2009年開始，為加大水資源管理力度并改善汾河流域的水污染狀況，山西省政府開始實行綠化政策工作，將流域的綠化與植樹造林相結合起來，組織各級政府進行人工植樹造林并加大其力度。

從2013—2017年高植被覆蓋區明顯降低，由2013年的17 836.13 km2減少到2017年的8 656.69 km2，但中、中高植被覆蓋區顯著增加，由2013年的18 668.75 km2增加到2017年的26 544.56 km2，在此期間由于汾河流域遭遇了旱情，造成了植被的退化，且在2017年山西省天氣冷暖異常，7月高溫天氣歷史最多，最高氣溫為41.9 ℃。

2017—2020年低、中低、中植被覆蓋區顯著減少，中高、高植被覆蓋度顯著上升，高覆蓋度達到了27.44%。2017年6月，山西省全面開啟了汾河流域生態保護綠色發展修復工作，先后出臺了各種修復方案和條例，截止目前，已初步實現了“水量豐富起來、水質好起來”。2020年汾河流域的植被覆蓋度出現上升的趨勢。

3.1.2 植被覆蓋空間分布及變化特征 在ArcGIS 10.2中將5期汾河流域植被覆蓋度圖進行分級，結果如圖3。由圖3可知，高植被覆蓋度主要分布在汾河流域邊界的縣市，沿著呂梁山脈和太行山脈分布，而處于兩條山脈之間的太原盆地植被覆蓋度受人文因素的影響則主要以中高植被覆蓋度為主，其周圍區域是低覆蓋度，但太原市主城區的植被覆蓋度明顯一直處于低覆蓋度的情況，忻定盆地的植被覆蓋度則一直處于中度覆蓋以下。

2005年，處于汾河下游地區的植被覆蓋度主要的植被覆蓋度類型很明顯為低植被覆蓋度，植被覆蓋率極低，是因為人們長期對水資源的過渡開采利用與污染。2005—2009年，植被覆蓋明顯改善的區域為汾河下游以南的區域地段，主要為從低植被覆蓋度轉變為中植被覆蓋度，對于忻州盆地以北部分區域，植被覆蓋度由低覆蓋轉變為中度覆蓋。2013年，植被覆蓋度顯著大幅改善，低、中低植被覆蓋度幾乎都轉變為中高、高植被覆蓋度，尤其表現在臨汾盆地和運城盆地，使汾河流域的高植被覆蓋度達到最大化，可以看出省政府對汾河流域的重視。2013—2017年，沿太行山脈西北區域的縣市、汾河流域下游以南地區，植被覆蓋度由原來的高植被覆蓋度變為中及中高植被覆蓋度，汾河流域一半以上的地區植被覆蓋度達到60%。2017—2020年，汾河流域管轄區主要對低植被覆蓋度的地區(植被覆蓋度在60%以下)進行大力整改，沿太行山脈的區域植被再次得到改善。

2005—2020年，汾河流域植被覆蓋類型分布如圖3所示，太原市主城區的植被覆蓋度雖在不同程度上得到改善，但幅度較小，因此是需要重點關注。

圖3 2005-2020年汾河流域各期植被覆蓋度分布

3.2 地形因子與植被覆蓋度的關系分析

地形因子的不同會造成水熱狀況的不同，從而對植被覆蓋度的分布格局具有不同的影響程度[14-15]。本研究對汾河流域的坡度、坡向、高程的劃分參考白建偉[5]對汾河流域的分級(表3)。2005—2020年汾河流域高程分布特征如圖4所示。

表3 2005-2020年汾河流域地形因子劃分指標

圖4 汾河流域高程分布特征

3.2.1 高程與汾河流域植被覆蓋度的關系 通過對汾河流域的DEM數據進行分級劃分，從而得到汾河流域的高程(圖4)。植被覆蓋度在不同的高程(海拔)有著不同的分布。

在ArcGIS 10.2中將轉換為矢量數據的高程圖與各時期的植被覆蓋度進行分區統計，計算得到各期不同高程下的植被覆蓋度均值(圖5)。由圖5可見，隨著高程的不斷增高，汾河流域的平均植被覆蓋度在不斷增加，且在大于2 000 m的高程地區汾河流域的平均植被覆蓋度達到最大，超過了85%，同時在1 500～2 000 m的地區平均植被覆蓋度達到了80%左右。得益于在高海拔地區人類無法干預植被的生長。在1 000～1 500 m的高程地區，汾河流域的植被覆蓋度超過60%但不足80%。而在小于1 000 m的地區，由于人類活動頻繁，對于植被覆蓋具有一定的破壞程度，其植被覆蓋度不超過60%，覆蓋類型主要為中低植被覆蓋度。

圖5 2005-2020年汾河流域各時期

2005—2020年汾河流域的平均植被覆蓋度在逐年增高，自山西省政府出臺《山西省汾河流域水資源管理和水環境保護條例》以來，汾河流域實施綠化工作，2013年汾河流域實施大量的人工造林和退耕還林政策，此年不同高程的平均植被覆蓋度略高于其余研究年份。

3.2.2 坡度與汾河流域植被覆蓋度的關系 利用DEM數據，依據表3的分級指標在ArcGIS 10.2下生成坡度圖(圖6)。坡度代表了地表單元的陡緩程度。在很大程度上，坡度會大幅度影響水分再分配的過程，在再分配過程中對土壤中的有機物和水分進行控制，從而影響植被的生長與空間分布[16]。

圖6 汾河流域坡度分布特征

在ArcGIS 10.2中將轉換為矢量數據的坡度圖與各時期的植被覆蓋度進行分區統計，計算得到各期不同坡度下的植被覆蓋度均值(圖7)。從圖7中可以看出，隨著坡度的不斷增大，汾河流域的平均植被覆蓋度呈現微弱的上升趨勢，但上升的程度不大。在坡度大于12°的區域，汾河流域的平均植被覆蓋度達到了最大，超過了70%，但不足80%，并沒有達到高植被覆蓋度。在9°～12°的坡度范圍內，其平均植被覆蓋度超過了60%。綜合來看，在坡度超過9°的區域，汾河流域的植被覆蓋度處于中高覆蓋度，而且在2013，2020年不同坡度的平均植被覆蓋度均超過了60%，2005—2020年汾河流域的植被呈現緩慢的改善趨勢。

圖7 2005-2020年汾河流域各時期不同坡度下的平均植被覆蓋度

3.2.3 坡向與汾河流域植被覆蓋度的關系 通過DEM數據，依據表3的分級指標在ArcGIS 10.2下生成坡向圖(圖8)。坡向就是坡面法線在水平面上的投影的方向 (也可以通俗理解為由高及低的方向)，不同的坡面朝向會使太陽輻射量及土壤水分的多少存在著明顯的差異，進而造成植被分布特征的不同[17]。

圖8 汾河流域坡向分布特征

在ArcGIS 10.2中將轉換為矢量數據的坡向圖與各時期的植被覆蓋度進行分區統計，計算得到各期不同坡向下的植被覆蓋度均值(圖9)。由圖9可知，坡向對于汾河流域植被覆蓋的分布沒有特別明顯的影響，平均植被覆蓋度基本都超過了60%。但在無坡向地區，汾河流域植被覆蓋度明顯低于其他坡向，陰坡地區的平均植被覆蓋度明顯高于其他坡向，相反在半陰坡、陽坡、半陽坡平均植被覆蓋度呈現稍下降趨勢。是因為坡向對植被覆蓋的主要影響體現在溫度和降水上[18]，陰坡的太陽輻射小，土壤水分大幅度保留，造就了植被的相對茂盛。

圖9 2005-2020年汾河流域各時期不同坡向下的平均植被覆蓋度

在各坡向的地區內，變化起伏不明顯，但總體上呈現小幅度上升的趨勢。2013年，各坡向的植被覆蓋度達到最大值，除無坡向外，其他坡向平均植被覆蓋度均超過了70%，且最大值為73.19%。綜上所述，坡向對汾河流域的植被覆蓋度的影響程度不大。

4 討論與結論

4.1 結 論

2005—2020年，汾河流域植被覆蓋度總體呈現增加趨勢，增幅相對比較明顯，2013年汾河流域的平均植被覆蓋度明顯高于其他年份，較為突出。2005—2020年，汾河流域植被覆蓋度主要以中高覆蓋度為主，其面積比例均超過了30%，其次為中高、高覆蓋度，低覆蓋度所占比例極低，表明汾河流域的植被覆蓋度狀況整體呈現良好的狀態。

(1) 從空間分布上來看，高植被覆蓋度主要分布在汾河流域周邊的縣市，沿著呂梁山脈和太行山脈分布，中高植被覆蓋度分布較為均勻，低、中低植被覆蓋度主要分布在臨汾、運城、太原市的主城區。從總體特征上來看，越靠近山區的區域植被覆蓋度越高，而處于盆地的區域植被覆蓋度不高，可能受人文因素的影響較為大。

(2) 隨著高程的不斷增高，汾河流域植被覆蓋度呈現明顯的上升趨勢，在大于2 000 m的區域，植被覆蓋度達到最高。

(3) 隨著坡度增大，汾河流域平均植被覆蓋度不斷上升，但上升的程度不大，在坡度大于12°時，植被覆蓋度最高。

(4) 坡向對汾河流域植被覆蓋度影響不明顯，但在陰坡處，植被覆蓋度較高。

4.2 討 論

汾河流域的植被覆蓋度在時間跨度上總體呈現逐步增加的趨勢，這與白建偉[5]、韓麗君等[6]對汾河流域植被覆蓋度的研究相符合。其次，汾河流域植被覆蓋度與高程和坡度存在著明顯的相關性，但相對于坡向而言，對植被造成影響的主要是太陽的輻射及土壤中水分，而汾河流域的坡向相對均勻一些，所以坡向對汾河流域的植被覆蓋度影響并不明顯。本文對汾河流域植被覆蓋度的研究存在著一定的局限性和不足，如對汾河流域植被覆蓋度的分布變化并未逐年進行研究分析，人類活動及氣候變化對植被覆蓋度的具體影響劃分研究等等。本研究對汾河流域植被覆蓋度的研究內容相對較少，其土地利用類型也可能為汾河流域植被覆蓋度的重要影響因素，但受數據等各方面的限制，影響因子只是結合地形因子進行探討，且因汾河流域相對較大，對研究數據的獲取和處理存在一定的限制。應在往后的研究中加入更多的影響因素分析，采用更便捷的數據處理方法，進而實現汾河流域最大程度的綠色發展。