趙 芳，林恭華，趙之重

(1.青海大學，青海 西寧 810016；2.中國科學院西北高原生物研究所，青海 西寧 810008)

三江源區(qū)位于青藏高原腹地，是青藏高原構成的主體[1]。三江源區(qū)是中國江河中下游地區(qū)和東南亞國家生態(tài)環(huán)境安全和區(qū)域可持續(xù)發(fā)展的生態(tài)屏障，同時也是中國生態(tài)系統(tǒng)最敏感和最脆弱的地區(qū)之一[2]。草地生態(tài)系統(tǒng)是三江源區(qū)的主體生態(tài)系統(tǒng)類型，對這一地區(qū)乃至全球氣候環(huán)境穩(wěn)定具有重要意義。三江源區(qū)以高寒植被為主體[3]，群落結(jié)構簡單，植被特征隨著時間和空間而變化。植被指數(shù)是對地表植被活動的簡單、有效的經(jīng)驗度量，分析三江源區(qū)植被指數(shù)時空變化對重要氣候因子的響應，有助于準確了解這一地區(qū)的植被現(xiàn)狀并為預測植被變化提供重要理論依據(jù)。

水熱條件(降水和氣溫)是表征氣候最重要的因子[4]，也是影響植被指數(shù)時空異質(zhì)性的重要因子[5]。當前，對于三江源區(qū)植被指數(shù)時空差異及其與水熱條件的關系研究，主要包括兩類方法。一類是以三江源區(qū)的整體或部分區(qū)域為空間單元，分析其植被指數(shù)在年際間的變化規(guī)律及其與水熱條件的關系[4,6-7]；另一類是以某一年或連續(xù)多年為時間單元，研究不同采樣點間植被指數(shù)的分布及其與水熱條件的關系[8-10]。然而，值得注意的是，三江源地區(qū)植被指數(shù)變化存在明顯的時空差異[11-12]，單獨對局部地區(qū)或三江源區(qū)整體進行時間動態(tài)分析可能會忽略空間水平上的特殊性，而單獨對個別年份進行空間差異分析則可能忽略時間水平上的異質(zhì)性細節(jié)。此外，值得一提的是，近10年來，在全球氣候變化的整體背景下，三江源區(qū)氣溫、降水特征有新的變化，對這一地區(qū)草地植被指數(shù)的時空變化可能存在新的影響。因此，利用最新資料對三江源區(qū)植被指數(shù)時空異質(zhì)性及其與水熱因素的關系進行系統(tǒng)分析是極其必要的。

大部分相關研究都是基于8 km精度的NOAA-AVHRR NDVI指數(shù)，2000年以來，精度更高(可達250 m)、對植被具有更高敏感度的MODIS Ⅵ產(chǎn)品(NDVI和EVI)逐漸成為草地監(jiān)測的主要基礎資料。然而，利用MODIS產(chǎn)品對三江源地區(qū)近10年來植被指數(shù)時空變化及其與水熱條件的關系研究，至今未見報道。本研究基于三江源區(qū)18個站點2000-2009共10年的氣溫、降水和MODIS(NDVI、EVI)數(shù)據(jù)，結(jié)合多種統(tǒng)計分析方法，重點探討三江源區(qū)植被指數(shù)的年際變化和空間變化對水熱條件變化的響應機制，為這一地區(qū)的草地形勢預測提供有價值的參考資料。

1 資料與分析方法

1.1衛(wèi)星遙感影像 采用NASA/MODIS產(chǎn)品MOD13Q1(MODIS/Terra Vegetation Indices 16-Day L3 Global 250m SIN Grid.)。因為三江源區(qū)牧草生長季為4月中旬-9月下旬，所以選用2000-2009年4月23日-9月30日合成的NDVI(Normalized Difference Vegetation Index)和EVI(Enhanced Vegetation Index)植被指數(shù)時空數(shù)據(jù)集。而且三江源區(qū)覆蓋兩景影像，所以每年共有20景影像。該數(shù)據(jù)可以從NASA’WIST數(shù)據(jù)網(wǎng)站下載，網(wǎng)址為https://wist.echo.nasa.gov/api/。

1.2氣象資料 收集了三江源區(qū)18個氣象臺站(圖1)2000-2009年年均溫和年均降水量觀測數(shù)據(jù)。

圖1 三江源區(qū)氣象臺站分布

1.3分析方法 首先，利用MRT(MODIS Reprojection Tool)對每年的原始影像進行同期鑲嵌處理，得到拼接后的10景NDVI和10景EVI影像；而后在ArcGIS 9.3中進行投影轉(zhuǎn)換，利用Spatial Analyst模塊中的Raster Calculator程序分別將10景NDVI和10景EVI影像合并為單景的NDVI和EVI均值影像(即對影像取平均值)。將18個臺站的經(jīng)緯度信息導入ArcGIS 9.3，使用Extract Values to Points工具，從每年NDVI和EVI均值影像中提取出每個臺站每年的NDVI像元值和EVI像元值。

為了充分探討植被指數(shù)對水熱條件的響應規(guī)律，研究分別進行了分年(多臺站)和分臺站(多年)數(shù)據(jù)分析過程。在分析植被指數(shù)對水熱條件的分年動態(tài)響應時，將每年18個臺站的NDVI像元值和EVI像元值以及對應臺站的年均溫、年降水量在PASW18.0中進行相關分析、偏相關分析等。同樣，在分析植被指數(shù)對水熱條件的分臺站響應時，將各個臺站10年的NDVI像元值和EVI像元值與對應年的平均氣溫、降水量在PASW18.0中進行相關分析、偏相關分析等。

2 結(jié)果

2.1三江源區(qū)植被指數(shù)變化

2.1.1年際變化趨勢 從三江源區(qū)18個臺站在2000-2009年植被指數(shù)NDVI和EVI的變化趨勢圖(圖2)中可以發(fā)現(xiàn)，NDVI和EVI的走勢基本一致，不同的植被指數(shù)在反映三江源區(qū)植被變化時具有很高的相關性(R2=0.939)和同步性。從植被指數(shù)的線性趨勢線可以看出，三江源區(qū)近10年來植被指數(shù)總體上呈上升趨勢；值得注意的是，NDVI和EVI與年份的Spearman相關性較弱(RNDVI=0.430，REVI=0.321,P>0.200)，推測可能是三江源區(qū)植被指數(shù)的年際波動幅度較大所致。

2.1.2空間分布差異 三江源區(qū)植被指數(shù)存在很大的空間差異，值最大的清水河臺站和值最小的五道梁臺站的NDVI、EVI值分別相差4 380.30和3 232.60(圖3)。將18個臺站按地理位置分為東部、中部和西部地區(qū)臺站，其中東部、中部和西部地區(qū)NDVI均值分別為4 050.37、3 502.85和2 151.72，EVI均值分別為2 683.16、2 391.45和1 364.92，呈現(xiàn)出相同的規(guī)律，即由東向西逐漸降低。而在區(qū)域內(nèi)部，同樣存在明顯的空間差異，如中部地區(qū)的班瑪臺站其NDVI和EVI分別是興海臺站的2.48和2.22倍。

2.2三江源區(qū)植被指數(shù)變化對水熱條件的響應

2.2.1植被指數(shù)年際變化對水熱變化的響應 從2000-2009年18個臺站年均降水、氣溫與對應臺站NDVI、EVI值的直接相關和偏相關結(jié)果(表1)來看，NDVI、EVI與降水和氣溫在每一年的相關性顯著程度都不同，表現(xiàn)出較明顯的時間差異。NDVI、EVI在除了2002年的其他所有年份，與降水呈顯著的正相關關系，表明降水是影響近10年三江源區(qū)植被的主導因子。同時，NDVI、EVI在10年內(nèi)與氣溫的相關關系都不顯著，偏相關結(jié)果又顯示相關關系為負值，表明近10年來，氣溫對三江源區(qū)植被的影響較小，而且是負向影響。

圖2 三江源區(qū)NDVI、EVI在2000-2009年的變化趨勢

圖3 三江源區(qū)18個臺站NDVI、EVI的空間差異

表1 2000-2009年18個臺站年平均降水、氣溫和植被指數(shù)的相關性

2.2.2植被指數(shù)空間分布差異對水熱變化的響應 從18個臺站2000-2009年年均降水、氣溫與對應年的NDVI、EVI值的直接相關和偏相關結(jié)果(表2)來看，在三江源區(qū)這個大區(qū)域內(nèi)，各個臺站植被指數(shù)與水熱條件的相關顯著性各不相同，有的臺站相關性極顯著，有的則不顯著，其中分別有7個和5個臺站的NDVI或EVI與降水和氣溫呈顯著相關關系，表現(xiàn)出很大的空間差異。降水和氣溫對植被指數(shù)的影響仍然不同，多數(shù)臺站植被指數(shù)與降水的相關性好于氣溫(如瑪多)，也有臺站氣溫的相關性好于降水(如班瑪)，有的則是降水和氣溫的相關顯著性都較高(如曲麻萊)，總體而言，三江源區(qū)降水與NDVI、EVI的相關性好于氣溫。

表2 18個臺站2000-2009年年平均降水、氣溫和植被指數(shù)的相關性

3 討論與結(jié)論

三江源區(qū)植被指數(shù)在2000-2009年間呈明顯上升趨勢，降水和氣溫也是逐漸增加和升高，這表明，三江源區(qū)在2000-2009年間植被狀況逐漸恢復，植被指數(shù)總體上對水熱條件的響應很高。同時，三江源區(qū)植被指數(shù)存在明顯的空間差異。18個臺站的NDVI和EVI值各不相同，有的甚至相差5～6倍，區(qū)域之間也有明顯的差異，植被指數(shù)呈從東部往西部逐漸降低的趨勢，降水、氣溫均值也是逐漸降低，區(qū)域之間植被指數(shù)對水熱條件有很好的響應。因此認為，三江源區(qū)地處高寒地區(qū)，氣候條件惡劣，特別是進入21世紀以來，降水和氣溫的大幅度波動和不確定性對植被指數(shù)的時空差異影響很大。

基于不同年份的相關分析結(jié)果顯示，三江源區(qū)植被指數(shù)與水熱條件的相關性存在時間上的差異。降水在某些年份與植被指數(shù)呈極顯著相關，而在有的年份卻不顯著，因此認為，有些研究用某一年或一個時間段的數(shù)據(jù)分析水熱條件對該區(qū)植被的影響，可能會影響結(jié)論的準確性。與此同時，將每年18個臺站的氣溫和降水取均值后與對應年的年均NDVI、EVI值進行直接相關和偏相關分析，結(jié)果表明相關性都不顯著；可見，植被指數(shù)對水熱條件響應規(guī)律的時間異質(zhì)性也隨空間地域的不同而變化，如若將三江源區(qū)大面積范圍內(nèi)的植被指數(shù)看作一個整體進行分析，也可能發(fā)生誤判的情況。分臺站的相關分析結(jié)果表明，植被指數(shù)與水熱條件的相關性在三江源區(qū)18個臺站間各不相同。如果片面地選取整個三江源區(qū)或某個局部地區(qū)進行分析，很容易導致錯誤的結(jié)論。將每個臺站10年的降水、氣溫取均值后與NDVI、EVI進行相關性分析，相關性不顯著，可見，空間異質(zhì)性也有時間段的差異。綜上所述，鑒于三江源區(qū)植被指數(shù)及其對水熱條件響應規(guī)律具有嚴重的時空差異，在這一地區(qū)進行相關分析時應充分考慮研究對象的時空特性。

唐紅玉等[4]和陳瓊等[7]研究顯示，氣溫是影響三江源區(qū)植被生長的主導因子。而本研究表明，降水對三江源區(qū)植被時空變化的影響大于溫度，即降水才是影響三江源區(qū)植被變化的主導因子。由于之前的相關研究所用的植被指數(shù)為8 km的NDVI數(shù)據(jù)，而本研究結(jié)合250 m的NDVI和EVI數(shù)據(jù)同時進行分析，精度和數(shù)據(jù)量都高于前者，結(jié)果應更能反映真實規(guī)律。值得一提的是，不論是分年份，還是分臺站，都發(fā)現(xiàn)降水總體上對植被指數(shù)呈正向影響，而氣溫總體上呈負向影響，具體原因有待進一步的分析。

[1]陳桂琛.三江源自然保護區(qū)生態(tài)保護與建設[M].西寧:青海人民出版社,2007:1-42.

[2]樊江文,邵全琴,劉紀遠,等.1988-2005年三江源草地產(chǎn)草量變化動態(tài)分析[J].草地學報,2010,18(1):5-10.

[3]王堃,洪紱曾,宗錦耀,等.“三江源”地區(qū)草地資源現(xiàn)狀及持續(xù)利用途徑[J].草地學報,2005,13:28-47.

[4]唐紅玉,肖風勁,張強,等.三江源區(qū)植被變化及其對氣候變化的響應[J].氣候變化研究進展,2006,2(4):177-180.

[5]李曉東,傅華,李鳳霞,等.氣候變化對西北地區(qū)生態(tài)環(huán)境影響的若干進展[J].草業(yè)科學,2011,28(2):286-295.

[6]徐維新,古松,趙新全.氣候持續(xù)變暖引起三江源地區(qū)植被出現(xiàn)階段性變化新特點[C].北京:第26屆中國氣象學會年會氣候變化分會,2009:943-955.

[7]陳瓊,周強,張海峰,等.三江源地區(qū)基于植被生長季的NDVI對氣候因子響應的差異性研究[J].生態(tài)環(huán)境學報,2010,19(6):1284-1289.

[8]于龍.基于“3S”技術的三江源區(qū)果洛州草地動態(tài)空間分析[D].西寧:中國科學院西北高原生物研究所,2009.

[9]馬軒龍.基于3S技術對青海省草地資源生產(chǎn)力的監(jiān)測[D].蘭州:蘭州大學,2008.

[10]米兆榮,張耀生,趙新全,等.NDVI和EVI在高寒草地牧草鮮質(zhì)量估算和植被動態(tài)監(jiān)測中的比較[J].草業(yè)科學,2010,27(6):13-19.

[11]楊建平,丁永建,陳仁升,等.長江黃河源區(qū)高寒植被變化的NDVI記錄[J].地理學報,2005,60(3):467-478.

[12]張鐿鋰,劉林山,擺萬奇,等.黃河源地區(qū)草地退化空間特征[J].地理學報,2007,61(1):3-14.