陳曉寧，蔣好忱

(國家測繪地理信息局第一地形測量隊(duì)，陜西 西安 710054)

植被在全球變化研究中起著“指示器”的作用，對植被的動(dòng)態(tài)監(jiān)測可以在一定程度上反映氣候的變化趨勢，是生態(tài)環(huán)境研究的熱點(diǎn)之一[1-4]。歸一化植被指數(shù)(normalized difference vegetation index,NDVI)是反映植被所吸收的有效輻射比例的重要指標(biāo)，對植被的長勢和生長量非常敏感，在一定程度上代表地表植被覆蓋變化[5]。因此，NDVI常用來估測土地覆蓋面積大小、植被光合作用能力、植被生物量、區(qū)域蒸散量、土壤水含量及旱情災(zāi)害等[6-13]。

許多學(xué)者研究了不同地區(qū)NDVI與主要?dú)夂蛞蜃拥年P(guān)系。文獻(xiàn)[14]利用SPOT VGT NDVI數(shù)據(jù)和氣象資料，分析了1998—2007年我國東部植被在不同季節(jié)對氣溫、降水變化的時(shí)空響應(yīng)特征；文獻(xiàn)[15]利用GIMMS/NDVI數(shù)據(jù)分析了我國西北地區(qū)1982—2006年植被NDVI時(shí)空變化特征及其影響因子；文獻(xiàn)[16]利用變化斜率法、主成分分析法和相關(guān)分析法得出流域NDVI與降雨的良好相關(guān)性；文獻(xiàn)[17]使用統(tǒng)計(jì)學(xué)方法，探討了NDVI時(shí)空變化及與氣候因子的關(guān)系。本文對秦嶺山區(qū)的多年NDVI變化情況進(jìn)行分析，同時(shí)使用相關(guān)分析法和趨勢分析法對不同植被類型下影響NDVI的氣候因子進(jìn)行研究，為氣候變化背景下制定環(huán)境資源保護(hù)和利用策略提供參考。

1 數(shù)據(jù)與方法

1.1 數(shù)據(jù)基礎(chǔ)

以陜西秦嶺地區(qū)為研究樣區(qū)，植被數(shù)據(jù)資料采用AVHRR-NDVI資料，空間分辨率為8 km，將數(shù)據(jù)進(jìn)行預(yù)處理拉伸并標(biāo)準(zhǔn)化至-1～1之間。選取安康、佛坪、商州、華山4個(gè)氣象站點(diǎn)(見表1)，獲取1982—2006年的氣溫和降水量數(shù)據(jù)，資料來自于中國氣象局氣象信息中心氣候資料室。研究范圍及氣象站點(diǎn)的空間位置分布如圖1所示，秦嶺地區(qū)植被覆蓋類型如圖2所示。

表1 氣象站點(diǎn)

1.2 研究方法

采用線性傾向估算法和相關(guān)系數(shù)法對秦嶺地區(qū)植被指數(shù)的時(shí)空變化及其與氣候因子的關(guān)系進(jìn)行分析。構(gòu)建樣本與時(shí)間之間的線性回歸方程，回歸系數(shù)與回歸常數(shù)的最小二乘估計(jì)見式(1)。利用回歸系數(shù)b與相關(guān)系數(shù)r之間的關(guān)系，求出時(shí)間ti與變量xi之間的相關(guān)系數(shù)，見式(2)。

(1)

(2)

NDVI與氣候因子的相關(guān)系數(shù)可通過式(3)得到

(3)

式中，n為樣本數(shù)量。

2 結(jié)果與分析

2.1 NDVI變化分析

對安康、佛坪、商州、華山4個(gè)氣象站的氣溫和降水?dāng)?shù)據(jù)進(jìn)行計(jì)算，求出1982—2006年的年平均氣溫、年總降水量和年平均植被指數(shù)。分別對1982—2006年的年平均氣溫、年總降水量及植被指數(shù)進(jìn)行擬合，結(jié)果如圖3—圖5所示。

1982—2006年間，秦嶺地區(qū)年均溫最低為10.7℃，最高為12.8℃，整體呈不斷波動(dòng)上升趨勢，25年間上升幅度達(dá)到了2.1℃。

由圖4可知，雖然各年的總降水量值大小不一，但就變化趨勢來看變化趨于平穩(wěn)，平均在790 mm左右。年總降水量的變化趨勢不明顯，約為-7.172 6 mm/a。

1982—2006年，NDVI的多年變化趨勢有所上升，反映了研究區(qū)域的植被有所增加，植被的保護(hù)相對較好，每年增幅約為0.000 7。

對植被在1982—2006年的變化作趨勢分析，結(jié)果如圖6所示。植被指數(shù)在研究區(qū)域的中西部和東北部都有明顯的增加(黑色和灰色所示)。在南部和西北部的部分地區(qū)植被指數(shù)值一度處于低值，甚至出現(xiàn)負(fù)值，最低值為-0.002 3。針對植被指數(shù)的變化情況，對1982—2006年植被指數(shù)變化的顯著性作了檢驗(yàn)，結(jié)果如圖7所示，植被指數(shù)的變化在中西部和東北部比較顯著，而在南部和西北大面積區(qū)域不顯著。圖7與圖6所示情況基本一致，說明在中西部和東北部植被指數(shù)的變化較為明顯。

將研究區(qū)植被類型數(shù)據(jù)與NDVI變化數(shù)據(jù)進(jìn)行疊加分析，進(jìn)一步研究25年來不同植被類型的變化情況(見表2)。從表2可以看到，增加最為明顯的植被為溫帶、亞熱帶、熱帶常綠灌叢矮林，平均值為0.001 2，而針葉林的變化最小。

表2 不同植被類型NDVI變化情況

2.2 不同區(qū)域NDVI與氣候因子相關(guān)性

將1982—2006年逐年NDVI值與同期的逐年年均溫進(jìn)行相關(guān)分析，相關(guān)系數(shù)如圖8所示。研究區(qū)域中西部和東北的大部分地區(qū)植被指數(shù)變化與氣溫之間呈明顯的正相關(guān)關(guān)系，相關(guān)系數(shù)最大值達(dá)到0.72，表明中西部和東北大部分地區(qū)的植被指數(shù)變化受氣溫的影響較大；南部和西北地區(qū)呈明顯的負(fù)相關(guān)關(guān)系(白色區(qū)域)，最大負(fù)相關(guān)系數(shù)達(dá)到-0.41，由此可見該區(qū)域植被指數(shù)變化受氣溫的負(fù)面影響較大；此外，局部小范圍地區(qū)(灰色所示)的植被幾乎不受氣溫的影響。對植被指數(shù)與氣溫之間關(guān)系的顯著性作分析，結(jié)果如圖9所示，氣溫對植被指數(shù)影響明顯的區(qū)域也主要分布于中西部和東北區(qū)域。

分析氣溫對不同植被覆蓋類型的影響程度(見表3)，相互關(guān)系最為顯著的植被類型為果樹園、經(jīng)濟(jì)林，平均值達(dá)到0.52，說明果樹園、經(jīng)濟(jì)林受氣溫影響大，針葉林則受氣溫的影響最小。

表3 不同植被類型下氣溫與植被指數(shù)的相關(guān)關(guān)系

采用與上述相同的方法，分析NDVI與年總降水量的相關(guān)關(guān)系(如圖10、圖11所示)。從圖中可以看出僅有南部與北部的小部分區(qū)域的植被指數(shù)與降水呈正相關(guān)(黑色所示)，并且最大正相關(guān)系數(shù)只有0.473；而其他大部分地區(qū)植被指數(shù)變化與降水呈明顯的負(fù)相關(guān)(白色所示)，最大負(fù)相關(guān)系數(shù)達(dá)到-0.516，因此大部分地區(qū)的植被指數(shù)受降水的負(fù)面影響較大。除個(gè)別地區(qū)外，絕大部分地區(qū)的植被指數(shù)與降水之間的關(guān)系不顯著。

由表4可見，降水與植被的相互關(guān)系最顯著的植被類型為闊葉林，平均值為-0.13，說明闊葉林受降水影響較大；相互關(guān)系最不顯著的植被類型為果樹園、經(jīng)濟(jì)林，平均值僅為-0.07，受降水的影響較小。

表4 不同植被類型下降水與植被指數(shù)的相關(guān)關(guān)系

3 結(jié)論與展望

NDVI是反映植被覆蓋的一個(gè)重要指標(biāo)，其變化是區(qū)域植被覆蓋狀況的體現(xiàn)，分析其與氣候因子的相互關(guān)系，具有重要的理論和實(shí)際意義。本文所得結(jié)論如下：

(1) 陜西秦嶺地區(qū)年均氣溫在1982—2006年間呈不穩(wěn)定性上升，上升幅度為2.1℃，這主要是由于近幾十年來，隨著溫室氣體的不斷排放，全球氣溫不斷上升；年總降水量的變化趨勢趨于平穩(wěn)，平均值在756 mm左右；該區(qū)域的植被指數(shù)有略微的上升。

(2) 研究發(fā)現(xiàn)，1982—2006年秦嶺地區(qū)的NDVI在中西部和東北部與氣溫呈正相關(guān)，且比較顯著，在南部和西北部呈明顯的負(fù)相關(guān)，植被指數(shù)受氣溫的影響比較大。北部和南部的部分地區(qū)的植被指數(shù)與降水呈正相關(guān)，其他區(qū)域與降水呈明顯的負(fù)相關(guān)，植被指數(shù)與降水的關(guān)系不顯著，受降水的影響較小。總的來說，氣溫是影響該地區(qū)植被指數(shù)變化的主要因素。

由于受到數(shù)據(jù)空間分辨率的限制，精度更高的分析結(jié)果有待進(jìn)一步的完善。