楊成樂+張永福

摘要：利用2010～2014年Landsat 8數(shù)據(jù)，通過差值和相關(guān)性分析法，對(duì)呼圖壁縣植被指數(shù)（NDVI）的時(shí)空演變進(jìn)行分析，并對(duì)NDVI與氣候因子相應(yīng)的敏感性進(jìn)行分析。結(jié)果表明，呼圖壁縣植被覆蓋變化呈現(xiàn)明顯的空間差異，植被指數(shù)有明顯的地域特征，中部平原地區(qū)的NDVI最大，植被覆蓋度明顯較高，主要分布在林地、農(nóng)田；南部的高山、丘陵及北部荒漠地區(qū)的NDVI較低，植被覆蓋度也較低；呼圖壁縣具有明顯的季節(jié)性NDVI變化，NDVI增加最快在夏季和秋季，而冬季增長(zhǎng)緩慢，且NDVI對(duì)氣候因子（降水、氣溫）響應(yīng)存在時(shí)滯性；呼圖壁縣植被覆蓋NDVI與降水和氣溫呈極顯著相關(guān)關(guān)系，相關(guān)系數(shù)分別為0.843 1、0.870 2；復(fù)雜地形對(duì)植被指數(shù)有影響，農(nóng)田耕作區(qū)的NDVI變化更大程度上體現(xiàn)了人類活動(dòng)對(duì)植被指數(shù)的影響。隨著人口增加、城市化進(jìn)程加快、城鄉(xiāng)建設(shè)侵占大量土地，會(huì)減少植被覆蓋的面積，影響植被指數(shù)的分布。

關(guān)鍵詞：呼圖壁縣；植被指數(shù)；氣候因子；Landsat 8

中圖分類號(hào)：X144 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：0439-8114（2016）20-5229-04

DOI：10.14088/j.cnki.issn0439-8114.2016.20.016

Abstract： Based on Landsat 8 data during 2010～2014，temporal and spatial evolution of NDVI and its relationship with climate factors in Hutubi county were analyzed by using difference and correlation analysis method. The results showed that Hutubi county was covered by the change showed significant spatial difference vegetation index with obvious regional characteristics. Central Hutubi vegetation condition was better than the northern and southern desert region of mountains，hilly areas，the distribution of woodland and fields，vegetation index was significantly higher. The NDVI developed faster in summer and autumn than that in winter，but the responses of vegetation variation to precipitation and temperature had time lag. There was significant relation between monthly average NDVI and monthiy precipitation and temperature（P<0.01），with the corrlation coefficients 0.843 1 and 0.870 2. The impact of complex topography on vegetation indices，NDVI changes in agricultural farming area reflects the impact of human activities on vegetation indices greater extent. With population growth，urbanization，urban costruction occuied a lot of land， and it will reduce the area of the vegetation cover，vegetation index distribution influences.

Key words： Hutubi county； vegetation index； climate factor； Landsat 8

植被是地理環(huán)境的重要組成部分，植被的變化在全球變化研究中起著“指示器”的作用[1]。在植被覆蓋變化研究中，歸一化植被指數(shù)（NDVI，Normalized difference vegetation index） 能很好地反映植被覆蓋、生物量及生態(tài)系統(tǒng)參數(shù)的變化[2]。植被覆蓋變化是氣候因子、環(huán)境因子以及人類活動(dòng)共同作用的結(jié)果[3，4]。李曉兵等[5]通過研究得出從東南到西北，中國(guó)NDVI的變化與溫度、降水變化的相關(guān)系數(shù)逐漸增大；戴聲佩等[6]通過研究證實(shí)降水是中國(guó)西北干旱半干旱地區(qū)植被生長(zhǎng)的主要限制因子。國(guó)內(nèi)學(xué)者在分析全國(guó)范圍植被變化與氣候關(guān)系的研究中也發(fā)現(xiàn)植被指數(shù)對(duì)氣候的響應(yīng)存在著區(qū)域差異[7-9]。

1 研究區(qū)概況

呼圖壁縣地形復(fù)雜，地勢(shì)南高北低，交錯(cuò)分布著高山、丘陵、平原、沙漠和河流，屬于溫帶大陸性干旱半干旱氣候，自然環(huán)境復(fù)雜，氣候多樣，是北疆重要的農(nóng)業(yè)縣，以種植小麥、玉米、棉花為主。然而土地利用類型不盡合理。目前鮮有對(duì)該區(qū)域植被覆蓋變化與氣候因子的響應(yīng)關(guān)系研究。本研究采用Landsat 8數(shù)據(jù)研究呼圖壁縣植被覆蓋變化與氣候因子的響應(yīng)關(guān)系，旨在增進(jìn)對(duì)呼圖壁縣生態(tài)環(huán)境的認(rèn)識(shí)，從而為該縣農(nóng)林業(yè)生產(chǎn)和生態(tài)保護(hù)提供科學(xué)依據(jù)。

呼圖壁縣位于新疆維吾爾自治區(qū)，轄6個(gè)鎮(zhèn)、 1個(gè)鄉(xiāng)，總面積9 721.6 km2。地處東經(jīng)86°5′～87°8′、北緯43°7′～45°20′。呼圖壁縣地勢(shì)南高北低，自東南向西北傾斜，縣內(nèi)地貌類型復(fù)雜多樣，南部為高山、丘陵，平均海拔2 400 m，占總面積的31.6%；中部為沖積平原，海拔460～700 m，是農(nóng)作物種植區(qū)，占總面積的43.2%；北部為沙漠地，海拔360～460 m，占總面積的25.2%。縣內(nèi)有兩大水系分別是呼圖壁河、雀爾溝河。

呼圖壁縣在中緯度西風(fēng)帶控制下，位于歐亞大陸中心帶，屬中溫帶大陸性氣候，冬長(zhǎng)夏短，冬寒夏熱，晝夜溫差大。境內(nèi)海拔高度相差較大，全縣南北地區(qū)的氣候也有著明顯的差異，低山、平原和沙漠地區(qū)屬中溫帶，南部中山和高山地區(qū)屬寒溫帶。溫度由北向南逐漸降低。

2 數(shù)據(jù)來(lái)源及預(yù)處理

Landsat 8數(shù)據(jù)是由中國(guó)科學(xué)院地理空間數(shù)據(jù)云平臺(tái)下載得到的，影像分辨率是1 km，時(shí)間分辨率是10 d。在ENVI5.0和ArcGIS10.0軟件下經(jīng)過大氣校正、輻射校正等圖像預(yù)處理，得到研究區(qū)的NDVI數(shù)據(jù)。

氣象數(shù)據(jù)來(lái)自于中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)提供的中國(guó)農(nóng)業(yè)基本氣象資料月值數(shù)據(jù)集，該數(shù)據(jù)集包含了研究區(qū)范圍內(nèi)及周邊的21個(gè)標(biāo)準(zhǔn)氣象臺(tái)站的月均值數(shù)據(jù)，對(duì)這21個(gè)臺(tái)站月均值數(shù)據(jù)通過克里金插值，對(duì)插值得到的溫度降水?dāng)?shù)據(jù)進(jìn)行投影轉(zhuǎn)換、圖像裁剪等處理，得到研究區(qū)所需氣象數(shù)據(jù)。

3 研究方法

本研究中通過一元線性回歸方法研究各變量的變化趨勢(shì)，皮爾森相關(guān)系數(shù)研究植被NDVI和溫度降水的相關(guān)關(guān)系，并且運(yùn)用顯著性檢驗(yàn)方式來(lái)驗(yàn)證研究結(jié)果。

1）NDVI均值法。利用NDVI柵格數(shù)據(jù)和經(jīng)過插值裁剪得到的氣象數(shù)據(jù)對(duì)研究區(qū)的年均溫、年均降水和年均NDVI進(jìn)行統(tǒng)計(jì)，以氣象部門季節(jié)劃分標(biāo)準(zhǔn)為依據(jù)。其中，春季為3～5月，夏季為6～8月，秋季為9～11月，冬季為12月至次年2月。

2）差值分析法。采用差值分析法即兩個(gè)時(shí)相NDVI相對(duì)差值衡量NDVI的空間變化大小[10]，計(jì)算每個(gè)柵格上的相對(duì)變化差，即

其中，X為相對(duì)差，j為前期，i為后期。

3）相關(guān)分析法。計(jì)算各月、季節(jié)、年份的平均NDVI，將其與相應(yīng)時(shí)段內(nèi)的降水量及氣溫因素進(jìn)行相關(guān)性分析，通過皮爾森公式分別計(jì)算出NDVI-溫度與NDVI-降水之間的相關(guān)系數(shù)，計(jì)算公式如下：

其中，rxy表示x、y之間的相關(guān)系數(shù)，xij表示第i年第j月的NDVI，x表示NDVI的月平均值、季度平均值或年平均值，yij表示第i年第j月的降水或溫度值，y表示溫度或降水的月平均值、季度平均值或年平均值。r介于-1～1之間，當(dāng)0

如果r通過了顯著性檢驗(yàn)P值（P<0.05），則NDVI與溫度以及降水的相關(guān)性都具有很好的統(tǒng)計(jì)學(xué)意義。

4 結(jié)果與分析

4.1 植被指數(shù)時(shí)空變化特征

4.1.1 空間植被指數(shù) 從宏觀上看，呼圖壁縣植被指數(shù)有明顯的地域特征，中部平原地區(qū)的NDVI最大，植被覆蓋度明顯較高，主要分布在林地、農(nóng)田。南部的高山、丘陵及北部荒漠地區(qū)的NDVI較低，植被覆蓋度也較低，如圖1所示。

4.1.2 季節(jié)植被指數(shù) 選取2013年3月至2014年4月各季節(jié)的NDVI平均值，從季節(jié)上分析NDVI的變化。由圖2可知，呼圖壁縣植被覆蓋隨季節(jié)呈現(xiàn)推移變化趨勢(shì)，季節(jié)變化明顯。

夏季的植被覆蓋情況最好，特別是呼圖壁縣中部平原地區(qū)，0.8

北部荒漠地區(qū)一年四季的NDVI基本沒有太大的變化，0.1

中部為沖積平原，種植農(nóng)作物，NDVI隨季節(jié)變化明顯?？梢钥闯觯S著月份逐漸增加，NDVI從0.1增加到0.8，再下降到0.1～0.2。

NDVI在季節(jié)變化中，春到夏季呈上升趨勢(shì)，表明這兩個(gè)季節(jié)促進(jìn)了植被活動(dòng)。春季增加較為明顯，而夏季增加相對(duì)微弱。夏季的植被覆蓋度高。秋季農(nóng)作物成熟收獲，收獲后植被覆蓋有所降低。

4.1.3 月植被指數(shù) 由呼圖壁縣各月NDVI的變化曲線（圖3）得知，在NDVI數(shù)月變化中，7、8、9月值較高。其中，8月NDVI達(dá)到最大，為0.8，由于夏季雨量充足，溫度適宜，植被生長(zhǎng)達(dá)到一年中的峰值，植被覆蓋達(dá)到一年中的最大范圍。1～3月由于冬季溫度過低，抑制植被生長(zhǎng)，NDVI較低，進(jìn)入4月后溫度逐漸上升，植被加快生長(zhǎng)，NDVI也逐漸增大，而進(jìn)入9～12月，溫度又降低，植被生長(zhǎng)速度變慢，NDVI逐漸減小。

4.2 植被指數(shù)變化的氣候原因分析

4.2.1 NDVI與月平均溫度的相關(guān)性分析 統(tǒng)計(jì)2010～2014年研究區(qū)內(nèi)各氣象站點(diǎn)相應(yīng)年份的月平均溫度，將該月平均溫度與相應(yīng)月份的NDVI進(jìn)行相關(guān)性分析。通過皮爾森公式計(jì)算得到NDVI與溫度的相關(guān)關(guān)系，從月平均溫度與NDVI關(guān)系曲線圖4來(lái)看，NDVI與溫度有高度的相關(guān)關(guān)系，相關(guān)系數(shù)為0.870 2。植被NDVI隨著溫度的增加而增加，這主要是因?yàn)橹脖簧L(zhǎng)需要一定的熱量條件。

根據(jù)2013～2014年各月平均溫度，統(tǒng)計(jì)得出一年內(nèi)月平均溫度，將其與各月平均NDVI進(jìn)行相關(guān)性分析，得到年內(nèi)NDVI與溫度的關(guān)系。年內(nèi)NDVI與溫度基本都呈正態(tài)分布，且兩者的相關(guān)性較好，NDVI隨著溫度的增加而增加，隨著溫度的降低而降低。7月溫度達(dá)到最高，而NDVI則在8月達(dá)到最高。出現(xiàn)了成長(zhǎng)期NDVI隨溫度變化存在滯后現(xiàn)象，有可能是NDVI受前期溫度的影響，或者是夏季溫度過高蒸發(fā)量變大。植物可能由于受到水分的限制導(dǎo)致生長(zhǎng)受到抑制。

但從總體上看，溫度與NDVI相關(guān)性較高，波動(dòng)較一致。在同一溫度范圍內(nèi)，植被類型是影響二者相關(guān)性的主要因子。

4.2.2 NDVI與月平均降水的相關(guān)性分析 統(tǒng)計(jì)2010～2014年研究區(qū)內(nèi)各氣象站點(diǎn)相應(yīng)年份的月平均降水，將該月平均降水與相應(yīng)月份的NDVI進(jìn)行相關(guān)性分析。通過皮爾森公式計(jì)算得到 NDVI與降水的相關(guān)關(guān)系。由月平均降水量與NDVI關(guān)系曲線（圖5）可以看出，NDVI與降水有高度的相關(guān)關(guān)系，相關(guān)系數(shù)為0.843 1，相對(duì)于NDVI與月平均溫度而言，NDVI與平均降水的相關(guān)性要大一些。降水量在0～150 mm時(shí)，NDVI隨著降水量的增加而增大，降水量超過150 mm時(shí)，NDVI隨著降水量的增加有下降趨勢(shì)，或因降水量增加但是溫度不適宜時(shí)，NDVI增加不明顯，冬季植被覆蓋度低，NDVI也低。

5 小結(jié)與討論

1）在空間變化尺度上，呼圖壁縣植被覆蓋變化呈現(xiàn)明顯的空間差異，植被指數(shù)有明顯的地域特征，中部平原地區(qū)的NDVI最大，植被覆蓋度明顯較高，主要分布在林地、農(nóng)田。南部的高山、丘陵及北部荒漠地區(qū)的NDVI較低，植被覆蓋度也較低。

2）在植被季節(jié)變化方面，春到夏季呈上升趨勢(shì)，表明這兩個(gè)季節(jié)植被活動(dòng)得到促進(jìn)。春季增加較為明顯，而夏季增加相對(duì)微弱。夏季的植被覆蓋度高。秋季農(nóng)作物成熟收獲，收獲后植被覆蓋有所降低。

3）NDVI與平均溫度和平均降水量有較高的相關(guān)性，NDVI與月平均溫度的相關(guān)系數(shù)為0.870 2，NDVI與月平均降水量的相關(guān)系數(shù)為0.843 1?？梢姕囟群徒邓渴怯绊懼脖桓采w的重要因素，但是NDVI與降水和溫度存在一定時(shí)滯性，夏季溫度過高蒸發(fā)量變大。植物可能由于受到水分的限制導(dǎo)致生長(zhǎng)受到抑制，降水量超過150 mm時(shí)，NDVI隨著降水量的增加有下降趨勢(shì)，或因降水量增加但是溫度不適宜時(shí)，NDVI增加不明顯。即NDVI是由溫度和降水共同作用的結(jié)果。

4）呼圖壁縣大體上處于地形較為平坦的地區(qū)，但是南部地區(qū)有高山和丘陵，所以也不能忽略復(fù)雜地形對(duì)植被指數(shù)的影響，農(nóng)田耕作區(qū)的NDVI變化更大程度上體現(xiàn)了人類活動(dòng)對(duì)植被指數(shù)的影響。隨著人口增加，城市化進(jìn)程加快，城鄉(xiāng)建設(shè)侵占大量土地，會(huì)減少植被覆蓋的面積，影響著植被指數(shù)的分布。

參考文獻(xiàn)：

[1] 孫紅雨，王長(zhǎng)耀，牛 錚，等.中國(guó)地表植被覆蓋變化及其與氣候因子的關(guān)系：基于NOAA時(shí)間序列數(shù)據(jù)分析[J].遙感學(xué)報(bào)，1998，2（3）：204-210.

[2] 戴聲佩，張 勃，王海軍.中國(guó)西北地區(qū)植被NDVI的時(shí)空變化及其影響因子分析[J].地球信息科學(xué)學(xué)報(bào)，2010，12（3）：315-321.

[3] 趙玉萍，張憲洲，王景升，等.1982年至2003年藏北高原草地生態(tài)系統(tǒng)NDVI與氣候因子的相關(guān)分析[J].資源科學(xué)，2009，31（11）：1988-1998.

[4] 信忠保，許炯心，鄭 偉.氣候變化和人類活動(dòng)對(duì)黃土高原植被覆蓋變化的影響[J].中國(guó)科學(xué)D輯：地球科學(xué)，2007，37（11）：1504-1514.

[5] 李曉兵，史培軍.基于NOAA/AVHRR數(shù)據(jù)的中國(guó)主要植被類型NDVI變化規(guī)律研究[J].植物學(xué)報(bào)，1999，41（3）：314-324.

[6] 戴聲佩，張 勃，王海軍，等.中國(guó)西北地區(qū)植被覆蓋變化驅(qū)動(dòng)因子分析[J].干旱區(qū)地理，2010，33（4）：636-643.

[7] 何 月，樊高峰，張小偉，等.浙江省植被NDVI動(dòng)態(tài)及其對(duì)氣候的響應(yīng)[J].生態(tài)學(xué)報(bào)，2012，32（14）：4352-4362.

[8] 高志強(qiáng)，劉紀(jì)遠(yuǎn).基于遙感和GIS的中國(guó)植被指數(shù)變化的驅(qū)動(dòng)因子分析及模型研究[J].氣候與環(huán)境研究，2000，5（2）：155-164.

[9] 劉 燦，高陽(yáng)華，李月臣，等.基于NDVI的重慶市植被覆蓋變化及其對(duì)氣候因子的響應(yīng)[J].長(zhǎng)江流域資源與環(huán)境，2013，22（11）：1514-1520.

[10] 劉憲鋒，任志遠(yuǎn).西北地區(qū)植被覆蓋變化及其與氣候因子的關(guān)系[J].中國(guó)農(nóng)業(yè)科學(xué)，2012，45（10）：1954-1963.