克里雅河流域植被覆蓋垂直地帶性分布與變化特征

阿不都拉·阿不力孜，依力亞斯江·努爾麥麥提，普拉提·莫合塔爾，白克拉木·孜克利亞

( 1．新疆大學旅游學院，烏魯木齊830049; 2．新疆大學綠洲生態教育部重點實驗室，烏魯木齊830046; 3.新疆大學資源與環境科學學院，烏魯木齊 830046；4.中央民族大學生命與環境科學學院，北京 100081)

全球氣候變化已在很大程度上影響著人類賴以生存的環境，加之人類日益頻繁的生產活動對自然界的干擾，加快了自然生態環境的演變速率，并導致了一系列環境惡化問題，如植被覆蓋的退化、土地沙漠化[1,2]。植被作為聯接土壤、大氣和水分的自然“紐帶”，受到水、土壤、熱、地形等環境要素的綜合影響，同時具有明顯的季節和年際變化及土地覆被變化的代表性，在氣候變化研究中起著“指示器”作用[3]。研究區域植被覆蓋與地形因子的關系對掌握該區域植被的垂直地帶性分布特征非常重要，給區域生態系統演變規律的探討以及氣候變化響應分析提供重要的科學依據。

植被覆蓋度作為主要的生態系統指標，區域植被蓋度變化體現了該區域自然和人類活動對生態環境的作用[4]。遙感具有大范圍連續觀測和數據獲取能力，是區域或全球植被覆蓋度估算的有效手段，其多時相、多樣化的空間和光譜分辨率特點，能提供不同尺度上的植被覆蓋信息[1]。隨著遙感技術的不斷發展，遙感已在諸多相關植被覆蓋研究中被廣泛應用，歸一化植被指數(NDVI)遙感數據集也越來越多，包括 AVHRR-NDVI、SPOT VGT-NDVI、MODIS-NDVI以及TM-NDVI等數據[5]。干旱區的植被覆蓋變化研究因該區域在全球變化研究中的重要性、區域水資源短缺性、年內降水分布不均勻性以及年際、季節溫差的高度變異性等特點成為研究的熱點區和重點區。目前，已有諸多學者利用多種遙感數據對植被覆蓋變化及其與氣溫和降水因子的關系進行了研究。Schmidt等[5]在以色列內蓋夫沙漠利用AVHRR-NDVI數據對其植被蓋度與降水量的關系進行了探討。信忠保等[6]利用GIMMS-NDVI和SPOT VGT-NDVI數據對黃土高原地區1981-2006年植被覆蓋的時空變化進行了研究，并從氣候變化和人類活動的角度分析了植被覆蓋變化的原因。馬明國等[7]利用1981-2001年的8 km AVHRR數據，采用均值法、差值法和一元線性回歸模擬法，分析了西北植被年內和年際變化特征，模擬了近21 a的 NDVI最大化值的變化趨勢線和NDVI增長幅度。吳云等[8]以MODIS-NDVI時間序列數據為基礎，利用像元二分模型對海河流域2000-2007年的植被覆蓋度進行了估算，分析了年最大植被覆蓋度的時空變化特征，并對植被覆蓋度與降雨量之間的響應關系進行了深入探討。袁麗華等[9]利用2000-2010 年的 250 m分辨率的MOD13Q1數據來研究黃河流域植被覆蓋區域的NDVI時空變化特征。徐興奎等[10]利用自然正交分解方法，基于NDVI數據，研究我國西北地區地表植被覆蓋的時空演變過程，并從氣候和人為兩大因素方面分析其變化原因。田源[11]基于Landsat TM/ETM+ 數據對新疆于田綠洲的植被覆蓋分布特征進行了研究，并分析了植被蓋度和溫度的關系。

流域是以水系為紐帶，將上、中、下游組成一個具有因果聯系的復合生態系統[12]。流域生態系統的植被覆蓋變化研究給氣候變化與生態系統演變相互作用關系研究提供基礎數據。克里雅河流域是典型的干旱區內陸河流域，對水資源短缺而開發強度極強的該流域來說，了解區域植被覆蓋度的分布特征以及演變規律，解析自然和人為因素對其影響是對區域生態環境的穩定性與農業可持續發展至關重要。本研究利用MODIS-NDVI時間序列數據和SRTM數據，基于像元二分模型對克里雅河流域2002-2013年間植被覆蓋度進行了估算，分析該流域植被覆蓋度時空變化特征及其與地形和氣象因子的關系。

1 材料與方法

1.1 研究區概況

克里雅河流域位居塔克拉瑪干沙漠南緣，昆侖山中段北麓，受大陸性干旱氣候和山盆相間的地貌格局影響，流域中部發育了典型的綠洲-荒漠生態系統[13]。流域內海拔最高點6 962 m，海拔高程5 000 m以上的山峰終年冰雪覆蓋，北部沙漠海拔高程在1 000 m以上。克里雅河流域自南向北可分為5個自然景觀帶，即山帶、低山丘陵帶、山前戈壁平原帶、洪積扇平原帶和沙漠帶[14](見圖1)。

圖1 克里雅河流域自然景觀帶分布Fig.1 Landscape zonal distribution in Keriya River Basin

四季分明、晝夜溫差大、降水稀少、蒸發量大、春夏多風沙和浮塵天氣是該流域的氣候特點，屬典型的極端干旱區[15]。平原綠洲區年均降水僅有14 mm左右，蒸發量則高達2 500 mm左右，綠洲主要依靠山區冰雪融水和部分地下水灌溉[11]。流域范圍內土壤鹽漬化和沙漠化現象共存，嚴重制約著植被生長分布和綠洲農業的發展，生態環境十分脆弱。流域氣候干旱、土壤貧瘠、多風沙的特點制約著植被的生長發育和分布，南部山區由于地貌和氣溫的影響，植被分布面積較小，區系較單一；而流域中部的洪積扇平原綠洲區植被受少降水、多風沙、浮塵為特點的極端干旱氣候的沖擊，除了人工培育的多種糧食、經濟作物和防護林外，還有蘆葦和檉柳等自然植被分布；北部沙漠區的主要植被種類有胡楊、檉柳、蘆葦、甘草、駱駝刺和羅布麻等[14]。

1.2 數據來源

本文所用數據包括遙感影像數據、氣象數據和縣志統計數據，均為在官方網站上正式公布的可靠數據，數據詳細清單如下：

(1)所用遙感影像數據為2002-2013年的MODIS-NDVI時間序列數據，空間分辨率為500 m × 500 m和30 m × 30 m的SRTM數據。其中，MODIS-NDVI數據來自于美國航空航天局(NASA)提供的MOD13A1數據集(http:∥glovis.usgs.gov)。SRTM高程數據從中國科學院計算機網絡信息中心(CNIC)地理空間數據云平臺下載(www.gscloud.cn)。為進行定量信息提取和分析，在ENVI和ArcGIS軟件環境中對遙感數據進行裁剪和投影轉換等預處理工作。

(2)研究區的氣象數據為1956-2012年的該區域降水、氣溫數據，從中國氣象局(CMA)國家氣象信息中心的中國氣象科學數據共享服務網絡平臺獲取(http:∥cdc.cma.gov.cn)。

(3)研究區相關的自然、地形地貌和地表覆被數據來源于于田縣志和中國知網(www.cnki.net)。

1.3 植被覆蓋度遙感估算

本文利用像元二分法模型來估算克里雅河流域的植被覆蓋度。首先，對MODIS-NDVI時間序列數據進行投影轉換，再利用于田縣行政邊界數據對其進行裁剪，最后用公式(1)對研究區植被覆蓋度進行估算并進行相關分析。

像元二分法模型[16]是線性混合像元分解模型中最簡單的模型，它假設像元值由植被與非植被覆蓋地表2部分構成[1]。模型原理是遙感傳感器觀測到的光譜信息S由綠色植被的光譜信息Sveg與裸土的光譜信息Ssoil2個組分因子線性加權合成。歸一化植被指數(NDVI)是植被生長狀況的最佳指示因子，并與植被覆蓋度具有良好的相關性。根據像元二分法，一個像元的NDVI值可表示為有植被覆蓋區域的NDVIveg和裸土部分NDVIsoil的線性加權。植被覆蓋度(fc)的計算公式(1)如下：

fc=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)

(1)

式中：NDVIveg為純植被覆蓋像元的NDVI值；NDVIsoil為裸土或無植被像元的NDVI值。在沒有實測數據的情況下，NDVIveg可以取遙感影像中NDVI的最大值；NDVIsoil可以取影像中NDVI的最小值。

2 結果與分析

在克里雅河流域2002-2013年的植被覆蓋度數據基礎上，統計分析該流域的平均植被覆蓋度的年際變化趨勢，并結合海拔高度差異進一步探討了不同區域的植被蓋度變化。

2.1 克里雅河流域植被覆蓋隨時間的變化特征

2002-2013年克里雅河流域的年平均植被覆蓋度隨時間的變化可以反映該流域的植被覆蓋年際變化(見圖2)。近12 a來，年平均植被覆蓋度值為25.67%～33.22%，其中，2002-2005年，植被覆蓋度呈現微小的減少趨勢，而在2005-2008年，植被覆蓋度變化幅度比較大，呈現一個周期的正弦式變化特征，2009-2012年，植被覆蓋度逐年上升，從2009年的28.80%提高到2012年的31.01%，然而在2013年，植被覆蓋度出現了猛減現象，下降到25.67%。除了2006年、2008年和2013年外，其余年份植被覆蓋度值年際變化為28.67%～31.01%，其幅度較小。

圖2 2002-2013年克里雅河流域平均植被覆蓋度年際變化特征曲線Fig.2 Annual mean vegetation coverage curve in Keriya River Basin from 2002 to 2013

變異系數[17]是總體樣本標準差與總體均值之比，即Cv=σ/μ，是衡量結構的風險程度或穩定程度的可靠性指標，是一個無量量綱。為進一步分析克里雅河流域的年平均植被覆蓋度隨時間變化的穩定性，對其進行變異系數的計算，其結果見表1。

克里雅河流域的年均植被覆蓋度隨時間的變異系數在一定程度上可以說明該流域植被群落與植被覆蓋變化的穩定性(見圖3)。在2002-2006年，植被覆蓋度變異系數由0.29減少為0.24，植被群落和覆蓋度變化趨于穩定。而從2006-2013年，植被覆蓋度變異系數由2006年0.24不斷增大為2013年0.37，植被群落與植被覆蓋的穩定性下降。

表1 2002-2013年克里雅河流域植被覆蓋度均值、標準差及變異系數Tab.1 Annual mean value, standard deviation and variation coefficient of vegetation coverage in Keriya River Basin from 2002 to 2013

圖3 2002-2013年克里雅河流域平均植被覆蓋度變異系數變化特征曲線Fig.3 Coefficient of vegetation coverage variation curve in Keriya River Basin from 2002 to 2013

2.2 克里雅河流域植被覆蓋度垂直地帶性分布特征與變化趨勢

為進一步分析克里雅河流域植被隨海拔高度的覆蓋分布特征，在不同海拔高度上對其2002-2013年的植被覆蓋度值進行均值處理，其結果見圖4。海拔高度1 000～1 500 m范圍內的平均植被覆蓋度為39.07%，提升到1 500～2 000 m時，其平均值下降到31.08%；2 000～3 500 m區域內，其平均植被覆蓋度依次增加，在3 000～3 500 m高度范圍內均值最高，其值為54.32%；在3 500～7 000 m區域內，區域平均植被覆蓋度呈現隨高度增加而減小的趨勢，其中6 500～7 000 m高度區域的平均植被覆蓋度值為19.39%，為最低植被覆蓋高度區(見圖4)。

圖4 克里雅河流域植被覆蓋度垂直地帶性分布特征Fig.4 Vertical zonal distribution characteristics of vegetation coverage in Keriya River Basin

不同海拔高度的植被覆蓋變化幅度會有差異，因此，對研究區不同海拔高度的平均植被覆蓋度在2002-2013年的最大值、最小值和變化幅度進行了統計分析(見表2)。

表2 2002-2013年克里雅河流域不同海拔高度植被平均、最大、最小覆蓋度及變化幅度Tab.2 Annual mean, maximum, minimum values of vegetation coverage and variations in Keriya River Basin from 2002 to 2013

圖5中可以看出，在過去12 a內，3 000～3 500 m海拔高度范圍的平均植被覆蓋度的變化差異最大，其變化幅度值為21.32%；其次，在1 000～1 500和2 500～3 000 m高度區的平均植被覆蓋度變化較大，其變化幅度值分別為12.01%和12.73%；其他海拔高度區域的植被雖存在平均覆蓋度的差異，但其年際變化幅度相當，其值為6.60%～8.30%(見表2)。

圖5 克里雅河流域不同海拔高度植被覆蓋度年際變化Fig.5 Annual vegetation coverage change at different altitude in Keriya River Basin

選取克里雅河流域上游高山帶、低山丘陵帶、中游沖積扇平原帶和下游沙漠帶為4個不同海拔高度區域植被覆蓋變化的典型代表，對其植被覆蓋度變化趨勢進行了分析(見圖6)。上游高山帶、中游沖積扇平原帶和下游沙漠帶的植被覆蓋變化趨勢線斜率為負值，植被覆蓋度值近12 a來呈下降趨勢，植被覆蓋逐漸退化。低山丘陵帶植被覆蓋變化趨勢線斜率為正值，植被覆蓋度值近12 a里有上升趨勢，植被覆蓋有所改善。其中，植被覆蓋度下降的3個區域的趨勢線斜率絕對值均高于低山丘陵帶區域趨勢線斜率的絕對值，說明流域內的植被覆蓋退化程度高于植被改善程度。

圖6 克里雅河流域上游高山帶、低山丘陵帶、中游沖積扇平原帶和下游沙漠帶的植被覆蓋度年際變化Fig.6 Annual change of vegetation coverage at upstream high mountain area, low hilly area, midstream alluvium and diluvium plain area, and downstream desert area of Keriya River Basin

3 植被覆蓋變化的影響因子分析

至今有諸多學者對全球或區域植被覆蓋變化及其影響因素進行了研究[3,4,9,18-20]。植被生長分布及其覆蓋率通常受到土壤類型、土壤水分、氣溫和降水條件以及高程等多種自然環境因子的制約，其中氣溫和降水是影響植被覆蓋的顯著性因子[3]。

3.1 克里雅河流域年均氣溫與年降水量變化

根據于田縣志記載[14]，克里雅河流域年平均氣溫東西方向變化小，西部略高于東部；南北之間的溫差比較大，北部地勢低，氣溫較高，由北向南隨著地勢的升高，氣溫逐漸降低。由圖7可知，在1956-1997年，克里雅河流域的年均氣溫微弱降低，而從1998年以后快速上升。在過去的57 a間，克里雅河流域的年均氣溫整體上呈現上升趨勢[0.13 ℃/(10 a)，R為0.137 6]。克里雅河流域的年降水量由北向南遞增。北部沙漠帶極其干燥，年降水量僅14 mm左右；中部的洪積扇平原農區為45.8 mm左右，南部低山丘陵帶達113 mm，海拔高度3 300～4 000 m地區的降水量可達150 mm。該流域的年降水量表現出增加的趨勢，然而年間的變化幅度較大。由圖8中的年降水量和降水距平變化曲線可知，在過去的57 a內，年降水量的總體增加率為3.30 mm/(10 a)，其相關系數R為0.027 9。

圖7 克里雅河流域年均氣溫及其距平變化Fig.7 Annual change of mean temperature and its anomaly in Keriya River Basin

圖8 克里雅河流域年降水量及其距平變化Fig.8 Annual change of rainfall and its anomaly in Keriya River Basin

3.2 植被覆蓋與氣象因子變化的響應關系

從克里雅河流域平均氣溫和年降水量變化分析結果可知，年均氣溫和年降水量均有微小的增加趨勢。在克里雅河流域4個不同景觀帶的植被覆蓋度變化數據的基礎上，其對氣溫和降水變化的響應進行分析和探討。

整體上克里雅河流域植被覆蓋與氣溫呈負相關，而與降水呈正相關。區域植被覆蓋變化以退化為趨勢的上游高山帶、中游沖積扇平原帶和下游沙漠帶的植被覆蓋度變化趨勢與降水的變化趨勢相對一致，與氣溫的變化呈現負相關關系。低山丘陵帶的植被覆蓋變化幅度較大，與降水的變化趨勢存在很強的正相關性，其趨勢相當一致，而對氣溫呈負相關關系(見圖9)。其中，中游沖積扇平原帶和下游沙漠帶的植被覆蓋度變化趨勢與氣溫和降水的相關度較差，可能是因為人類生產活動干擾引起的。

圖9 克里雅河流域植被覆蓋、氣溫和降水量的距平變化Fig.9 Anomaly change of vegetation coverage, temperature and rainfall in Keriya River Basin

4 結 語

本文以MODIS-NDVI時間序列數據和SRTM高程數據為基礎，基于像元二分模型對克里雅河流域2002-2013年的植被覆蓋度進行估算，并分析了該流域植被覆蓋垂直地帶性分布特征，探討了植被覆蓋度時空變化趨勢及其與氣象因子的相互關系。本文得出以下幾點結論。

(1)2002-2013年克里雅河流域的年平均植被覆蓋度值為25.67%～33.22%，呈減少趨勢，植被群落與植被覆蓋的穩定性下降。

(2)流域植被覆蓋存在明顯的垂直地帶性分布特征，其海拔為3 000～3 500 m范圍的植被覆蓋度均值最高，其值為54.32%，而6 500～7 000 m的平均植被覆蓋度值為19.39%，為最低植被覆蓋高度區。在3 500～7 000 m區域內，區域平均植被蓋度呈現隨高度的增加而減小的趨勢。在過去12 a內3 000～3 500 m海拔高度的平均植被覆蓋度的變化差異最大，其變化幅度值為21.32%。

(3)上游高山帶、中游沖積扇平原帶和下游沙漠帶的趨勢線斜率為負值，植被覆蓋度值近12 a來呈下降趨勢，植被覆蓋逐漸退化。低山丘陵帶趨勢線斜率為正值，有上升趨勢，植被覆蓋逐漸改善。

(4)總體上，流域植被覆蓋與氣溫變化呈負相關關系，而與降水變化表現出正相關。然而，上游高山帶、中游沖積扇平原帶和下游沙漠帶的植被覆蓋度與降水的變化趨勢相對一致，與氣溫的變化呈現負相關關系。低山丘陵帶的植被覆蓋變化幅度較大，與降水的變化趨勢存在很強的正相關性，其趨勢相當一致，而對氣溫呈負相關關系。