陽昌霞，阿的伍各，張 春

（1.成都理工大學(xué)地球科學(xué)學(xué)院，四川成都610059；2.鹽源縣林業(yè)和草原局，四川涼山615700）

地表的植被是遙感記錄和觀測的第一表層，遙感圖像中能夠反映出直觀的信息，是地球的重要組成之一，區(qū)域的植被特征與氣候、地貌及土壤條件密切相關(guān)[1]，能夠反映全球氣候變化，因此，利用有效的監(jiān)測方法正確監(jiān)測植被對人類活動具有幫助性意義.

遙感技術(shù)可從多時相、多波段提取地表覆蓋狀況[2]，具有大面積同步觀測，客觀性、時效性、綜合性和可比性，并且有經(jīng)濟性[3]，而傳統(tǒng)調(diào)查需要很多人力物力，花費很久的時間才能獲得地表大區(qū)域動態(tài)變化的數(shù)據(jù)，所以遙感技術(shù)不僅有時效且更經(jīng)濟.

國內(nèi)學(xué)者在植被覆蓋時空變化趨勢方面做了重點研究.曹旭娟等利用遙感技術(shù)的特性，利用空間數(shù)據(jù)對地表空間變化和定量化研究，定量化研究了區(qū)域生態(tài)系統(tǒng)的變化[4].李磊磊等分析了西藏自治區(qū)植被NDVI 與地表溫度上升密切相關(guān)，且西藏中部主要受降雨因子影響，西藏東部與西藏西部區(qū)域主要受到地表氣溫的影響，得出氣候變化是當(dāng)?shù)厣鷳B(tài)系統(tǒng)的響應(yīng)關(guān)系[5].劉振元等利用GIMMS/NDVI和MODIS/NDVI 遙感數(shù)據(jù)分析了青藏高原植被NDVI 變化規(guī)律及其影響因子[6].陸晴等利用GIMMS/NDVI 數(shù)據(jù)對青藏高原1982-2013 年高寒草地覆蓋時空變化及其對氣象因素的響應(yīng)進行研究，結(jié)果表明：青藏高原高寒草地生長季NDVI表現(xiàn)為從東南到西北逐漸減少的趨勢[7].楊志剛等采用Slope 趨勢分析方法和相關(guān)性分析方法，分析一江兩河地區(qū)的植被變化特征及其主要氣候驅(qū)動因素[8].樸世龍等對青藏高原草地植被覆蓋變化及其與氣候因子的關(guān)系進行研究[9].崔慶虎等進行了青藏高原草地退化原因述評[10].白淑英等進行了基于時間序列遙感數(shù)據(jù)的西藏山南地區(qū)植被覆蓋變化特征分析[11].高清竹等選擇草地植被覆蓋度作為評價草地退化的遙感監(jiān)測指標(biāo)，建立了藏北地區(qū)草地退化遙感監(jiān)測與評價的指標(biāo)體系[12].

國外學(xué)者Tucker 利用美國國家海洋和大氣管理局的氣象衛(wèi)星系列操作系列先進的極高分辨率輻射計傳感器的數(shù)據(jù)用于對非洲的土地覆蓋進行分類并監(jiān)測19 個月的非洲植被動態(tài)，研究綠葉植被密度和范圍的季節(jié)變化與熱帶輻合區(qū)運動相關(guān)的降雨模式之間存在對應(yīng)關(guān)系[13].Fensholt 等基于時間序列利用GIMMS與MODIS全球NDVI的變化監(jiān)測對全球植被變化進行了評估[14].

IPCC 第五次評估報告中指出，1880-2012 年期間全球地表平均氣溫增加了0.85 ℃，增溫速度超過了過去任何時期，并且全球地表平均溫度具有明顯的年代際和年際變化的特點[15]，其中根據(jù)遙感手段獲取的植被指數(shù)被廣泛應(yīng)用于植被監(jiān)測中，但是隨著遙感的技術(shù)在不斷地創(chuàng)新以及進步，傳統(tǒng)的獲取植被覆蓋度的方式已經(jīng)不能滿足如今實時性的需求，因此在做大尺度范圍的植被覆蓋動態(tài)變化分析時常采用低空間分辨率的NDVI 產(chǎn)品，其中MODIS是Terra 和Aqua 衛(wèi)星上搭載的主要傳感器之一，具有時間分辨率高等優(yōu)點，因此數(shù)據(jù)在大尺度植被分析上的應(yīng)用廣泛[16].

藏北地區(qū)植被變化對全球氣候變化的響應(yīng)具有顯著性，并且藏北地區(qū)的植被變化不但直接影響藏北牧區(qū)的社會經(jīng)濟發(fā)展，而且監(jiān)測預(yù)警草原沙化面積，水土流失，進一步為生態(tài)環(huán)境保護提供科學(xué)性的建議以及理論支撐，實現(xiàn)社會經(jīng)濟健康穩(wěn)定.本文以藏北地區(qū)植被變化對氣候變化的長期響應(yīng)為出發(fā)點，結(jié)合MODIS-NDVI 數(shù)據(jù)處理與分析，綜合分析了植被綠期與氣候因子的相互關(guān)系及其響應(yīng)特征以及對相關(guān)性特征的影響，旨在為研究區(qū)植被的科學(xué)管理及植被對氣候變化響應(yīng)的分析提供依據(jù).

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)藏北地區(qū)位于西藏自治區(qū)北部（圖1），主要包括那曲縣、安多縣、聶榮縣、比如縣，嘉黎縣、巴青縣、索縣、班戈縣、申扎縣、尼瑪縣共十個縣[16]，土地總面積約44.6萬km2（約66 901.2萬畝），是青藏高原的重要組成部分.

圖1 研究區(qū)藏北地區(qū)概況

研究區(qū)北面是昆侖山和唐古拉山，南面與南木林縣相鄰，東面與生達、八宿縣相鄰，西面與改則縣和措青縣相鄰.海拔高程大約在4 500 m 以上.藏北高原是我國長江、怒江、瀾滄江的發(fā)源地，對氣候變化以及人類的不合理開發(fā)和利用極為敏感，高寒草地分布廣泛，草原占自治區(qū)內(nèi)天然草地面積的59%.研究區(qū)內(nèi)冰川、雪山面積約0.91萬km2（約1 364.9萬畝），占全區(qū)土地總面積的2.04%.受高原地形的影響，氣候的突出特點是寒冷干燥，冬季漫長寒冷；受大氣環(huán)流和地形的影響，藏北地區(qū)的降水量從東向西、從南向北減少.藏北地區(qū)海拔高，光照充足，太陽能充足.

1.2 數(shù)據(jù)采集與預(yù)處理

（1）遙感數(shù)據(jù)：選取了2005 年至2015 年共11年的MODIS-NDVI 數(shù)據(jù).NDVI 數(shù)據(jù)分辨率為250 m，屬于MOD13Q1 數(shù)據(jù).MODIS-NDVI 數(shù)據(jù)來源于美國宇航局（National Aeronautics and Space Administration, NASA），經(jīng)過了水、云、重氣溶膠等處理.此數(shù)據(jù)按照統(tǒng)一算法開發(fā)的間隔16天的植被指數(shù)的最大值合成，能夠覆蓋整個研究區(qū)域.獲取數(shù)據(jù)后，將數(shù)據(jù)進行MRT 格式轉(zhuǎn)化，統(tǒng)一數(shù)據(jù)坐標(biāo)系統(tǒng)，轉(zhuǎn)換為WGS84 地理坐標(biāo)系，并將數(shù)據(jù)重采樣為( 1×1 )km分辨率.

（2）氣象數(shù)據(jù)：由于地理位置使得藏北地區(qū)的氣象站點較少（僅安多縣、班戈縣、申扎縣、那曲縣、索縣、比如縣6 個氣象站），因此，本文利用中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心獲取自1982 年以來逐年年降水量空間插值數(shù)據(jù)集以及逐年年均溫度空間插值數(shù)據(jù)集.將空間插值數(shù)據(jù)進行投影變化、重采樣、裁剪等處理，并通過ArcGIS 軟件下的柵格計算器計算11年間的平均溫度和平均降水量.

（3）DEM 數(shù)據(jù)：獲取了研究區(qū)內(nèi)90 m 分辨率的DEM 數(shù)據(jù)，統(tǒng)一數(shù)據(jù)坐標(biāo)系統(tǒng)，并進行裁剪獲取研究區(qū)的DEM數(shù)據(jù).

（4）除上述數(shù)據(jù)之外的其他主要數(shù)據(jù)：全國省級和縣級行政區(qū)劃.

1.3 研究方法

（1）最大值合成法：利用MVC（Maximum Value Composite）計算得到NDVI 的季節(jié)變化及NDVI 年份變化的最大值：

式中：Ii是第i 月的NDVI 值，Iij為第i 個月第j 旬的NDVI 值.最大值合成法（MVC）原理是對一旬的NDVI 值進行統(tǒng)計，選擇一旬內(nèi)天氣最晴朗、云量最少、對植被覆蓋影響最小的一天的最大值作為這一旬的NDVI 值，以降低衛(wèi)星與地表之間如云、大氣、太陽高度角等各種干擾程度.

（2）變異指數(shù)：利用變異指數(shù)CV（Coefficient of Variation）得到NDVI在時間序列上的穩(wěn)定性：

式中：xi為第i 年植被最大的NDVI 值，xˉ是這11 年間的NDVI 變化的平均值.變異系數(shù)CV 越大說明數(shù)據(jù)分布越離散，越不穩(wěn)定，NDVI 年季變化越大，反之則表明數(shù)據(jù)分布的較為緊湊，NDVI 時間序列數(shù)據(jù)較為穩(wěn)定[17].

（3）Slope 趨勢分析法：計算出每個象元間的變化情況，也就是可以反映出NDVI的變化趨勢：

式中：slope 是植被變化的斜率，n 為年份數(shù)，n=11，Ii是第i年的NDVI值.

（4）相關(guān)性分析法：通過要素對比，可得出對植被變化的主要影響因子，即固定其中一個變量來計算其余兩個變量之間的相關(guān)性大小，對比分析得主控因子.利用格網(wǎng)法計算出各個變量的元素屬性，同時在EXCEL中計算元素間的相關(guān)系數(shù)，并利用簡單相關(guān)系數(shù)值可看出相應(yīng)元素間的相關(guān)性；通過簡單相關(guān)系數(shù)進一步計算偏相關(guān)系數(shù)，由圖像可以直觀看出相關(guān)性分布，綜合分析降水，氣溫對植被變化的影響[17]，偏相關(guān)系數(shù)的計算公式如下：

其中：rab、rac、rbc分別表示變量a 與b、a 與c、b 與c 之間的簡單相關(guān)系數(shù)表示將變量c 固定后變量a 與b之間的偏相關(guān)系數(shù).則為正相關(guān)，反之則為負相關(guān)，偏相關(guān)系數(shù)的絕對值越大，則表示該象元處二者要素相關(guān)性越強.簡單的相關(guān)系數(shù)的計算方法如式（5）所示：

式中：Rxy為要素x 和y 之間的簡單相關(guān)系數(shù)，表示x與y 的之間的相關(guān)性程度有多高，其取值范圍介于[-1,1]之間；xi、yi分別表示為要素x和y在第i年該要素的值；xˉ和yˉ分別表示為兩個元素的年間平均值.

（5）格網(wǎng)法：因為研究區(qū)面積約44.6 萬km2，以柵格作為評價單元進行統(tǒng)計分析的方法數(shù)據(jù)量過大，并且容易因數(shù)據(jù)的誤差導(dǎo)致分析結(jié)果不準(zhǔn)確，并且格網(wǎng)法可以減弱因研究區(qū)的氣象站點少而造成的研究誤差.因此建立格網(wǎng)，分析植被對氣候與地形因子的靜態(tài)響應(yīng)，以表格的形式記錄每一分析單元的平均值，用其中一個元素為橫坐標(biāo)，另一元素為縱坐標(biāo)，繪制出變量間的關(guān)系散點圖（質(zhì)心在網(wǎng)格內(nèi)的像元納入統(tǒng)計范圍）以及線性回歸線，并以此為依據(jù)進行分析討論.

2 結(jié)果與分析

2.1 植被覆蓋變化時空特征

2.1.1 月平均植被覆蓋變化特征

利用藏北地區(qū)2005-2015 年每年每個月的平均NDVI 值，得到這11 年間的每年的NDVI 變化的趨勢間的情況（如圖2）.從圖2 可知，藏北地區(qū)植被覆蓋在1-5月份基本不變，5月份之后NDVI呈現(xiàn)增加的趨勢至7 月份，且增加的幅度明顯，7 月份之后呈現(xiàn)不明顯下降至7月中旬上升至8月，在7-8月份之間NDVI 值達到最大，8 月份之后則快速下降至10月份，10月之后有回升狀態(tài)，但是上升幅度不明顯.

圖2 藏北地區(qū)2005-2015年月均NDVI值變化圖

年間的NDVI 呈現(xiàn)如此波動的原因是藏北地區(qū)的凍土從5 月初溫度升高時開始解凍，并且降雨量也增加，植被開始發(fā)芽生長，所以5 月開始NDVI 快速上升，植被快速增長到8 月份為植被長勢最好的階段，8 月份之后氣溫降低，降雨量也隨之減少，植被開始停止生長并向干枯，枯萎的趨勢發(fā)展，所以從8月下旬至9月初期階段NDVI值降低.根據(jù)藏北地區(qū)植被生長的實際情況，確定藏北地區(qū)植被生長季為5-8 月，5-10 月份的NDVI 對藏北地區(qū)的植被覆蓋具有代表性，變化趨勢為單峰型.

2.1.2 季平均植被覆蓋變化特征

利用藏北地區(qū)2005-2015年每月的NDVI，由最大值合成法（MVC）在ENVI 軟件中合成每年四季（理論意義上的四季，即春季為1-3 月，夏季為4-6月，秋季為7-9 月，冬季為10-12 月）的值，得到各個季節(jié)從2005-2015 年的植被變化趨勢，如圖3 所示.藏北地區(qū)2005-2015 年春季NDVI 值比較穩(wěn)定，基本為0.1，即植被覆蓋度低且隨著年份的變化趨勢不明顯.夏季變化明顯，NDVI 值在0.15～0.2 之間波動.秋季NDVI 波動明顯，可以看到2005-2007 三年間NDVI 值是呈下降的趨勢，到2008 年時NDVI 驟升.冬季相較秋季變化的相對穩(wěn)定，與春季相近，變化范圍在[0.08,0.13]之間.

圖3 藏北地區(qū)2005-2015季平均NDVI平均值變化圖

2.1.3 年平均植被覆蓋變化特征

為了分析藏北地區(qū)2005-2015 年間植被NDVI年際變化趨勢，對藏北地區(qū)這11 年月NDVI 進行最大值合成得到年NDVI 平均值，如圖4 所示，2005-2015 的11 年間藏北地區(qū)植被總體呈下降趨勢，除2009-2011年之外，波動比較小.

圖4 年均NDVI值變化趨勢

由圖4 可知，藏北地區(qū)年NDVI 范圍在[0.145,0.267]變動，多年平均值為0.243.2005-2008 年間，植被覆蓋自2005年有些許下降之后呈現(xiàn)為NDVI從0.25呈穩(wěn)定下降的狀態(tài)，2009年度有小幅度的上升，2009 年開始驟降至2010 年的0.22 之下，之后上升，而從2011 年到2015 年NDVI 值又一直呈現(xiàn)下降的趨勢，但植被覆蓋NDVI值仍≥0.25.可見藏北地區(qū)的植被覆蓋呈現(xiàn)下降的趨勢，草地大部分呈現(xiàn)退化的趨勢，這與高清竹對西藏草地呈嚴重退化趨勢的研究結(jié)果一致.

2.1.4 植被覆蓋空間分布特征

根據(jù)已有研究的植被覆蓋分級，并結(jié)合藏北地區(qū)的實際情況，將植被覆蓋分為5 個等級，分級標(biāo)準(zhǔn)如表1 所示.基于像元計算得到2005-2015 年藏北地區(qū)平均年度NDVI的空間分布（圖5）.

表1 植被覆蓋分級表

圖5 2005-2015研究區(qū)年平均NDVI植被覆蓋分布圖

由圖5 可知，藏北地區(qū)2005-2015 年間植被覆蓋的平均NDVI 為0.24，植被覆蓋呈現(xiàn)出從西向東增加的趨勢，NDVI 小于0.3 的低植被覆蓋區(qū)域分布于藏北地區(qū)的北部，其中還有些地方?jīng)]有植被的覆蓋，為無植被覆蓋區(qū).高植被覆蓋區(qū)主要分布于藏北東南部分.

圖5 表明，藏北地區(qū)主要以低植被覆蓋區(qū)與中植被覆蓋區(qū)為主，其中無植被覆蓋地區(qū)面積為8 191 km2，零星分布在藏北地區(qū)的中部和西北部地區(qū)；低植被覆蓋地區(qū)域面積約為30萬km2，主要分布于西北部、北部、申扎縣大部分地區(qū)；中植被覆蓋區(qū)域面積約為7.7 萬km2，主要分布于西南部、南部及東南部的大部分地區(qū)如安多縣、那曲縣、聶榮縣；高植被覆蓋區(qū)域面積為2.58 萬km2；極高植被覆蓋區(qū)域的面積最小，僅為385 km2，主要分布在氣候適宜，海拔較低的東南部地區(qū).

2005-2015 年藏北地區(qū)生長季中的植被呈東南向西北遞減的趨勢，由于藏北地區(qū)氣候梯度差異明顯，氣溫以及降水量由東南向西北隨海拔和緯度的升高而逐漸下降，導(dǎo)致藏北地區(qū)植被NDVI 表現(xiàn)出較強的空間抑制性.計算結(jié)果表明，藏北地區(qū)2005-2015 年11 年間的變異系數(shù)的平均值為0.05（圖6），變異系數(shù)較小，說明藏北研究區(qū)的NDVI 數(shù)據(jù)變異程度低，因此這11 年間的NDVI 值較為穩(wěn)定，數(shù)據(jù)可靠性高.

圖6 藏北研究區(qū)變異指數(shù)空間分布圖

2.1.5 植被覆蓋變化趨勢特征

利用Slope 趨勢分析方法對植被變化的趨勢進行分析計算，對藏北地區(qū)植被覆蓋退化趨勢進行分析.植被覆蓋退化分級的依據(jù)見表2.

表2 植被覆蓋退化分級表

對藏北地區(qū)退化趨勢進行分析得到如圖7 所示結(jié)果.可見藏北地區(qū)北部植被呈中度退化趨勢，面積約為11.3 萬km2；嚴重退化和輕微退化的地區(qū)的面積占了19 萬km2左右面積；基本不變的地區(qū)集中在西部和中部地區(qū)，面積約為7.1 萬km2；有改善的地區(qū)為南部和東南部地區(qū)，面積約為8.99萬km2.

圖7 藏北地區(qū)2005-2015年年均NDVI變化趨勢圖

退化地區(qū)主要是尼瑪縣的北部地區(qū)、班戈的北部地區(qū)；基本不變的地區(qū)主要是申扎縣大部分地區(qū)和尼瑪縣南部地區(qū)以及班戈縣南部地區(qū)，同時還有安多縣的西北部地區(qū)；植被覆蓋改善的有安多縣東部地區(qū)、聶榮縣、那曲縣、嘉黎縣、巴青縣、比如縣以及索縣.

2.2 植被NDVI對氣候變化的響應(yīng)

2.2.1 氣溫變化特征

藏北地區(qū)因其特殊地理位置，年平均溫度介于-2.8 ℃～1.6 ℃，冬季寒冷且時間長；2005-2015年這11 年間的溫度介于-11.24 ℃～8.65 ℃之間，平均氣溫為-1.85 ℃（圖8）.藏北地區(qū)高海拔地區(qū)尼瑪縣和班戈縣的北部多年平均溫度低，東部低海拔的索縣最高，呈現(xiàn)出北部高海拔地區(qū)溫度低，南部低海拔地區(qū)溫度高的趨勢，藏北地區(qū)的溫度變化大，溫度總體呈現(xiàn)上升的變化（圖9）.

藏北地區(qū)在11 年間年均氣溫出現(xiàn)過兩次明顯的波動，分別為2008 年和2011 年，氣溫急劇降低使得藏北地區(qū)凍土面積增大而不利于植被的生長，2010 年溫度的明顯升高會加速冰川和永久性積雪的融化，增加當(dāng)?shù)氐目衫盟Y源，對農(nóng)牧發(fā)展產(chǎn)生有利的影響，但是在帶來有利影響的同時對會改變藏北地區(qū)的水域環(huán)境，對生態(tài)脆弱的地區(qū)造成一定的影響，在降水量不變的時候可能會對植被的生長產(chǎn)生不好的影響，即抑制植被的生長.

圖8 藏北地區(qū)年平均溫度分布圖

圖9 藏北地區(qū)2005-2015年溫度變化折線圖

2.2.2 降水量變化特征

藏北地區(qū)由于其特殊的地理位置，平均年降雨量在247.3 mm～513.6 mm 之間，年蒸發(fā)總量大于自然降水量[14]，2005-2015 年的年平均降水量為425 mm.降水量的增加，使得藏北地區(qū)的濕潤指數(shù)增加，也可使表土層水分增加，對藏北地區(qū)的草地退化和地表沙化有一定的抑制作用，可以減緩草地的退化速度，尤其是干旱半干旱的西部地區(qū)，降水的增加更能減緩草地退化的狀況[18].圖10 的年均降水量分布顯示，降水量由東到西依次遞減，班戈縣北部高海拔地區(qū)、尼瑪縣西部與申扎縣等高海拔草原無人區(qū)年均降雨量最少，可以明顯看到藏北地區(qū)東南部聶榮縣、那曲縣、比如縣、嘉黎縣、巴青縣以及索縣等低海拔地區(qū)的年降水量較大.

將2005-2015 年間藏北地區(qū)的降水量作圖分析得到圖11，可見藏北地區(qū)的降水量處于波動下降的趨勢，由散點圖和擬合線可知，擬合程度高.其中這11 年來藏北地區(qū)降水量有明顯的3 次降水量最大，分別為2005 年、2008 年、2011 年，其值為500 mm、497 mm、503 mm，最小出現(xiàn)在2015 年（266.25 mm）.氣候暖濕向變化不僅對該地區(qū)的動植物生長提供了有利的條件，也能夠使草地的退化速度減慢維持生態(tài)系統(tǒng)的穩(wěn)定.由于研究區(qū)冰川多，常年凍土的面積大，溫度升高會改變研究區(qū)的水系狀況，在帶來有利的生態(tài)環(huán)境同時，長期的溫度升高也會加劇藏北地區(qū)的生態(tài)系統(tǒng)脆弱性，帶來不利的影響.

圖10 藏北地區(qū)年均降水量分布圖

圖11 藏北地區(qū)2005-2015降水量變化折線圖

2.2.3 植被NDVI與氣候因子的相關(guān)性分析

為了研究植被覆蓋對氣候變化的動態(tài)響應(yīng)，本文計算了每個像元的NDVI與年均氣溫與年均降水量間的簡單相關(guān)系數(shù)以及偏相關(guān).基于拓展時間序列進行逐像元相關(guān)分析可知，通過格網(wǎng)法得到藏北地區(qū)2005-2015 年的11 年平均NDVI和年平均溫度以及平均降水量的散點圖.

對植被NDVI與溫度的偏相關(guān)性分析（圖12）表明，植被NDVI與溫度間的系數(shù)為-0.27～0.32，在散點圖上可以看出藏北地區(qū)的植被與溫度呈連續(xù)的線性相關(guān)，即溫度升高植被NDVI也增加.

圖12 藏北地區(qū)年平均NDVI與年均氣溫間相關(guān)性散點圖

植被NDVI 與溫度整體表現(xiàn)出一次函數(shù)的關(guān)系，其斜率k=0.0423，相關(guān)系數(shù)R2=0.3005，且根據(jù)ArcGIS計算得出年均NDVI與年均溫度間的偏相關(guān)系數(shù)的平均值為-0.11，可見植被NDVI 與溫度的相關(guān)性不是很高，原因是整個藏北地區(qū)的溫度相對都較低，2005-2015 年平均溫度僅為-1.85 ℃，常年低溫，植被的變化與氣溫的變化相關(guān)性不是很明顯.其中偏相關(guān)性分布與分區(qū)分別如圖13a、13b 所示，藏北地區(qū)尼瑪和班戈縣北部因為海拔較高，氣溫較低、寒冷等特點以及山岳冰川和廣泛發(fā)育的多年凍土使得北部地區(qū)的生態(tài)系統(tǒng)脆弱，造成草地退化，因此北部地區(qū)植被NDVI 與氣溫的變化的相關(guān)性較強；東南部索縣、比如縣、嘉黎縣海拔較低，氣溫較高，植被長勢較好，植被的生長依賴氣溫的穩(wěn)定，所以藏北地區(qū)東南部地區(qū)植被NDVI與氣溫的變化正相關(guān)性較強.整個藏北地區(qū)的植被NDVI 與年均氣溫間的相關(guān)性分布中，呈負相關(guān)的地區(qū)約占71.3%，呈正相關(guān)的地區(qū)占28.7%，其中主要呈負相關(guān)的地區(qū)是位于南部地區(qū)的申扎縣和東南部地區(qū)的那曲縣.

NDVI 植被變化與降水量的偏相關(guān)系數(shù)為-0.35～0.46（圖14a），由散點圖15 可知，NDVI 值隨著降水量的升高呈連續(xù)的線性變化，整體表現(xiàn)出一次函數(shù)y=0.0009x-0.1521 的關(guān)系，其中斜率k=0.0009，相關(guān)系數(shù)為R2=0.6292，散點與線性的擬合度較強，基于軟件計算偏相關(guān)性可知平均偏相關(guān)系數(shù)為0.0025，由此可知植被覆蓋變化與降水量間呈正相關(guān)性，且相關(guān)性強，即隨著降水量的增加，植被長勢越好，植被覆蓋度越高.其中藏北地區(qū)的西北部地區(qū)屬于不敏感和輕度敏感的土壤侵蝕低敏感地區(qū)，降水少，起伏度小，即坡度較緩，海拔高，因此降水量的變化多少對植被覆蓋變化的影響以及相關(guān)性不明顯；東南部申扎、聶榮、那曲、安多、比如等縣降水量高，地勢起伏度大，屬于土壤侵蝕敏感區(qū)域，植被覆蓋與該地區(qū)的降水量的變化之間的相關(guān)性明顯且呈正相關(guān)，植被覆蓋及長勢隨著降水量的增加而增加；中部地區(qū)尼瑪和班戈縣的中部為東西部過渡區(qū).如圖14b 所示，其中整個藏北地區(qū)的植被NDVI與年均氣溫間的相關(guān)性分布，呈負相關(guān)的地區(qū)約占38.8%，呈正相關(guān)的地區(qū)占61.2%，因此整個藏北地區(qū)的植被NDVI與降水量總體呈正相關(guān).

圖13 藏北地區(qū)2005-2015年NDVI值與年均溫度的相關(guān)性分布與分區(qū)圖

圖14 藏北地區(qū)2005-2015年NDVI值與年均降水量的偏相關(guān)性分布與分區(qū)圖

圖15 藏北地區(qū)年平均NDVI與年均降水量間的相關(guān)性散點圖

雖然植被覆蓋變化與溫度和降水量都變現(xiàn)出正相關(guān)關(guān)系，但由R2可知，植被覆蓋的變化對于降水量的相關(guān)性明顯大于年均氣溫的相關(guān)性，由此可見藏北地區(qū)的植被覆蓋變化是溫度以及降水量共同作用影響的結(jié)果.

3 討論與結(jié)論

本文根據(jù)MODIS-NDVI 數(shù)據(jù)、氣象數(shù)據(jù)、DEM 數(shù)據(jù)，采用MVC 方法，通過變異系數(shù)的計算、Slope 趨勢分析、簡單相關(guān)系數(shù)的計算以及偏相關(guān)系數(shù)的計算方法，研究了藏北地區(qū)2005-2015 年間的植被覆蓋月變化、季節(jié)間的變化以及年份間的變化，進行了植被變化特征以及變化趨勢的分析，研究植被變化與氣候之間的響應(yīng)關(guān)系，得到如下結(jié)論：

（1）根據(jù)MODIS-NDVI 數(shù)據(jù)分析得到藏北地區(qū)植被覆蓋呈降低的趨勢，根據(jù)月NDVI 確定藏北地區(qū)植被生長季為5-8月，變化趨勢為單峰型，根據(jù)季NDVI 變化可知植被對季節(jié)的變化敏感，有特定的生長季，夏秋兩季的植被覆蓋明顯有起伏.年平均NDVI 為-0.23～0.90，平均值為0.24.植被覆蓋變化從研究區(qū)的西北部向研究區(qū)的東南部呈現(xiàn)增加的趨勢，總體植被覆蓋朝著良好的改善方向走，改善的區(qū)域為東南部海拔較低的溫暖濕潤地區(qū)，中部地區(qū)的植被沒有明顯的改善情況，西北部地區(qū)植被處于退化趨勢.

（2）藏北地區(qū)2005-2015 年的年平均氣溫和年降水量東南部地區(qū)明顯高于西北部地區(qū)，氣溫上升趨勢，降水量呈小幅度下降趨勢，總體氣候呈暖濕方向發(fā)展.同時藏北地區(qū)近11 年來溫度具有上升趨勢，溫度波動明顯，降水量小幅下降.

（3）植被覆蓋變化與氣候因子（溫度和降水量）呈正相關(guān)性.可以看出藏北地區(qū)東南部植被變化成覆蓋較好的情況，西北部邊緣、北部邊緣、中部地區(qū)植被覆蓋低.藏北地區(qū)的植被覆蓋變化不是受一個單一的因子的影響，而是由多個因子共同作用，但是在總體驅(qū)動中，氣候因子響應(yīng)中降水量的影響大于氣溫的影響.

在選擇時間序列上選取了11 年的數(shù)據(jù)，研究時間尺度較小，在對氣候變化的研究中，氣候因子只考慮了年均氣溫及降水量，沒有對其他因子如濕度、蒸發(fā)量、太陽的日照情況等進行研究討論.研究植被覆蓋變化對氣候因子的響應(yīng)，利用的是相關(guān)性和偏相關(guān)性方法進行分析，而氣候因子等各個因子間的相關(guān)性并不局限于簡單的線性相關(guān)，今后的研究過程中可考慮更多的方法對植被與氣候因子間的響應(yīng)研究.