石羊河流域上游植被時空變化及其對區域氣候的響應

唐 志 光,馬 金 輝,李 朝 奎,彭 煥 華,梁 繼

(1.湖南科技大學地理空間信息技術國家地方聯合工程實驗室，湖南 湘潭 411201；2.蘭州大學西部環境教育部重點實驗室，甘肅 蘭州 730000；3.蘭州大學遙感與地理信息系統研究所，甘肅 蘭州 730000)

石羊河流域上游植被時空變化及其對區域氣候的響應

唐 志 光1,馬 金 輝2,3,李 朝 奎1,彭 煥 華1,梁 繼1

(1.湖南科技大學地理空間信息技術國家地方聯合工程實驗室，湖南 湘潭 411201；2.蘭州大學西部環境教育部重點實驗室，甘肅 蘭州 730000；3.蘭州大學遙感與地理信息系統研究所，甘肅 蘭州 730000)

利用1999-2013年的SPOT-VEGETATION NDVI時間序列數據和降水、氣溫數據，對石羊河流域上游近15年的植被變化特征及其對降水和氣溫的響應關系進行了研究。結果表明：研究區NDVI的空間分布具有明顯的垂直地帶性；近15年中，81.3%的地區植被生長季NDVI呈增長趨勢；研究區平均NDVI在春、夏、秋季及植被生長季的增長幅度分別為11.75%、9.62%、5.98%和7.75%，植被生長季在延長。氣候變化對研究區植被的影響隨植被類型及其生長季節的不同而具有明顯差異，體現出較大的時空異質性；降水是影響研究區植被變化的主要因素。

石羊河流域；植被；NDVI；時空變化；氣候變化

0 引言

植被是連接土壤、大氣和水分的自然紐帶，在水土保持及區域氣候研究中扮演著重要角色[1,2]。植被與氣候間的相互作用主要表現在兩方面：1)氣候是決定地球上植被狀況的最主要因素，而植被是氣候最鮮明的反映和綜合標志；2)不同的植被類型通過影響植被與大氣間的物質(如CO2和水等)和能量(如太陽輻射、熱量等)交換影響氣候，改變后的氣候又通過植被與大氣間的這種物質和能量的交換對植被生長產生影響，最終影響全球的碳平衡格局。因此，監測植被的動態變化及其氣候響應關系，已成為全球變化研究的重要內容。

石羊河流域是我國西北干旱區頗受關注的水資源及生態環境脆弱地帶[3]，并且是甘肅省河西內陸河流域中人口最多、水資源開發利用程度最高、用水矛盾最突出、生態環境問題最嚴重、水資源對經濟社會發展制約性最強的地區[4]。石羊河流域上游祁連山區是其河流的發源地和水源形成區，更是水資源短缺的中下游地區的用水來源地。該流域的植被也大都集中在上游的祁連山區，作為石羊河流域陸地生態系統的代表，對石羊河流域的水源涵養起著重要的作用。因此，研究石羊河流域上游植被的時空變化及其氣候響應關系，對有效應對氣候變化，實現石羊河流域水資源的可持續發展及遏制荒漠化的發展態勢具有重要的現實意義，為研究和治理石羊河流域乃至干旱區內陸河流域的生態環境問題提供可靠的理論依據。

歸一化植被指數(NDVI)是檢測植被生長狀態、植被覆蓋度的重要遙感指標；NDVI遙感產品已被廣泛地應用于全球及區域尺度上的植被時空變化監測及其對氣候變化的響應關系研究。自1987年以來，我國西北地區正經歷著由暖干向暖濕轉變的氣候轉型，氣溫和降水均呈現出明顯的增長趨勢，而且這種氣候轉型的態勢將延續下去[5,6]。現有研究表明，在氣候變化的背景下，我國西北地區的植被也正發生顯著的變化[1,7-12]。雖然這些研究因研究范圍、選取數據以及時間跨度的不同，所得結果存在一定的差異，但多數研究[1,9-12]表明我國西北大多數地區的植被NDVI在過去近20年呈現增長的趨勢。在這樣的背景下，定量研究石羊河流域上游的植被變化趨勢及其氣候響應機制更具有必要性和迫切性。本研究利用長時間系列的SPOT-VEGETATION NDVI數據集、氣溫和降水數據，旨在分析近15年石羊河流域上游(由東到西包括大靖河、古浪河、黃羊河、雜木河、金塔河、西營河、東大河、西大河8個子流域，面積約10 567 km2)植被的時空變化特征及其對區域氣候變化的響應。

1 數據與方法

1.1 數據及預處理

1.1.1 NDVI數據 本文所利用的植被NDVI數據集為1999年1月—2013年12月的SPOT-VEGETATION 逐旬NDVI產品，來源于http://free.vgt.vito.be/home.php網站。該數據由瑞典的Kiruna 地面站負責接收，并經過了嚴格的質量控制和預處理，最終生成逐日1 km 的全球數據[13]。其預處理包括輻射定標、幾何校正、大氣校正，以及利用清除云和大氣干擾的最大化合成法(MVC)[14]生成10 d 最大化合成的NDVI數據。該數據集已被廣泛應用于植被動態變化研究中[8,15-18]。

1.1.2 氣象數據與處理 本研究中的氣象數據主要包括1998-2013年研究區域內及周邊的39個氣象臺站的降水和氣溫資料。為了方便逐像元分析研究區NDVI變化對降水、氣溫的響應關系，利用與回歸模型相結合的克里格插值方法對氣象臺站的氣溫和降水數據均進行了空間插值處理，插值后的空間分辨率為1 km(與NDVI數據一致)。通過建立降水/氣溫的回歸模型初步反演研究區降水/氣溫的空間分布,并利用克里格法對降水/氣溫回歸模型的殘差值進行修正。降水/氣溫回歸模型的基礎數據為研究區DEM(90 m)數據提取的地形因子,結合氣象站點的降水/氣溫分別建立其多元線性回歸方程。將各氣象站點的降水/氣溫觀測值減去多元線性回歸方法的預測數據得到殘差值,并對殘差進行普通克里格空間插值;最后用殘差空間插值結果修正回歸模型估算的降水/氣溫數據，得到修正后的研究區降水/氣溫空間分布數據。通過交叉驗證的方法，得出研究區降水和氣溫空間插值結果具有較高的精度，其平均誤差分別為9.25%和7.73%。

在研究區的高海拔地區，由于氣象臺站數據稀少，降水/氣溫的插值結果可能存在較大的偏差；然而，這并不會影響后文中的分析結果，因為在后文中這些高海拔的非植被覆蓋區未參與相關計算。

1.1.3 土地利用、DEM數據 研究中的植被類型劃分是以甘肅省2000年1∶10萬土地利用數據集(來源于寒區旱區科學數據中心，http://westdc.westgis.ac.cn)為基礎，選用2008-2011年的Landsat5-TM遙感數據并結合高分辨率的Google Earth圖像，采用目視解譯的方法，最終按照研究需要(探討不同植被類型NDVI對氣候變化響應的差異性)將研究區土地利用類型劃分為林地、草地、耕地和非植被類型，其中非植被類型為高海拔的冰雪、裸巖等無植被覆蓋區,不參與植被變化的分析。還采用90 m分辨率的SRTM DEM數據，主要用于提取地形因子輔助氣溫、降雨量數據的插值處理等。

1.2 主要分析方法

一元線性回歸可模擬每個像元的變化趨勢，對區域內的逐像元模擬則可綜合反映區域的時空格局演變。該方法也常應用于植被變化的定量分析中[8,19,20]。本研究采用該方法定量計算研究區NDVI的變化趨勢，公式如下：

(1)

式中：ΘSlope表示線性擬合回歸方程的趨勢線斜率，n為參與計算的總年份數，i為年序號，NDVIi表示第i年的春季平均NDVI。在本研究中通過式(1)計算得到的變化趨勢圖反映的是1999-2013年研究區的NDVI的總體變化趨勢。當ΘSlope>0 時，則表示NDVI在15年間呈增加的趨勢。

計算了近15年研究區NDVI的百分比變化幅度(η)，以反映近15年的NDVI變化量，公式如下：

(2)

對于植被NDVI與氣溫、降水之間相互關系的研究主要采用計算其皮爾遜相關系數的方法。

2 結果與討論

2.1 研究區NDVI空間分布特征

根據研究區1999年1月-2013年12月的逐旬SPOT-VEGETATIONNDVI數據集，計算得到石羊河流域上游近15年的年均NDVI、年最大值NDVI分布圖(圖1,見封3)。為了進一步分析研究區植被覆蓋與海拔高程間關系，將各柵格年最大值NDVI與其對應的海拔高程作散點圖(圖2)。由圖1和圖2可以看出，植被NDVI的分布呈現明顯的垂直地帶性。從東北向西南隨著海拔高度的增加NDVI值先增大后減少，植被呈現東南-西北方向的帶狀分布，年均NDVI小于0.1的地區主要分布在高海拔的冰雪、裸巖等無植被覆蓋區。在3 200m以下的地區，植被NDVI隨海拔的升高而增大，而在海拔3 200m以上，植被NDVI隨海拔的升高而減小。植被覆蓋較高的地區(年最大值NDVI>0.6)主要在海拔2 800～3 500m的地區(圖3)。植被覆蓋呈現這種垂直地帶性規律的原因可能有：在海拔2 800m以下的地區降雨量較少，氣溫相對較高，水分蒸發嚴重，因而植被覆蓋較差，而且隨著海拔升高降水量增加，植被覆蓋也增強；在海拔2 800～3 500m的地區降雨量較豐富，氣溫也適合植被生長，因而植被較茂盛；而在海拔3 500m以上，隨著海拔的升高氣溫逐漸降低，低氣溫成為制約植被生長的主要因素，因而植被覆蓋逐漸減少；到了海拔4 000m以上植被覆蓋稀疏或基本無植被，主要是冰雪、裸巖及其他未利用地。

圖2 年最大值NDVI與海拔高程關系Fig.2 The relationship of annual maximum NDVI and altitude

圖3 研究區年均NDVI、年降水量及年均氣溫的年際波動Fig.3 Interannual variations of annual mean NDVI,annual precipitation and annual mean temperature in the study area

2.2 研究區氣候及NDVI年際波動特征

從1999-2013年研究區年均NDVI、年降水量及年均氣溫的變化(圖3)可知近15年研究區NDVI呈增加趨勢，總體年均增長速度為0.002(R2=0.41),該植被NDVI增加的趨勢與我國西北地區植被活動的相關研究結論[9,12]一致。研究區年降水量和年均氣溫也呈增加趨勢，增長速度分別為3.138 mm(R2=0.19)和0.0194℃(R2=0.23)。其中，年均NDVI與年降水量的年際波動特征大致接近，體現出明顯的相關性。年均NDVI與年降水的相關系數(R)為0.74，而與年均氣溫的相關系數(R)僅為0.15。

2.3 研究區NDVI的時空變化特征

為了反映研究區NDVI變化的時空格局特征，分別對春、夏和秋季以及整個生長季的NDVI進行逐像元的趨勢分析與百分比變化幅度計算。根據研究區NDVI的年內波動特征，本研究中的植被生長季定義為3-11月，即包含春、夏、秋季。圖4顯示了1999-2013年研究區平均生長季NDVI的變化趨勢及其變化幅度。按照不同植被類型分區，統計了近15年各植被類型在不同季節的NDVI變化趨勢及變化幅度(表1)。

從圖4(見封3)可以看出，大部分地區的平均生長季NDVI體現出增加的趨勢，說明在近15年中大部分地區的植被覆蓋增強，植被狀況得到良好的改善。植被增強的區域面積比例占81.3%，其中研究區北部的大部分區域(包括東大河、金塔河以及西營河流域的下游等)植被改善較明顯，NDVI的變化趨勢斜率大于0.004(即年均NDVI的增長值大于0.004)，近15年變化幅度大于10%，并且 NDVI增長幅度大于20%的地區占研究區總面積的25%。植被退化的區域面積比例為19.7%，主要分布在植被覆蓋度低的高海拔地區。但是，從整個研究區看，近15年石羊河流域上游植被NDVI的總體趨勢是在增強，植被狀況得到改善。

表1 15年不同植被類型的NDVI變化趨勢(Θslope)及百分比變化幅度(η)Table 1 Annual increment of NDVI and increase rate of NDVI for different vegetation types during 1999-2013

從表1可以看出，近15年研究區平均生長季NDVI的增長趨勢及變化幅度分別為0.0025 a-1、7.75%。從不同植被類型看，林地、草地和耕地的平均生長季NDVI均呈增加的趨勢，其中，耕地增長最快，其近15年的平均生長季NDVI的增長趨勢為0.003 a-1，變化幅度為10.64%。從不同季節看，研究區NDVI增長最快的為夏季，達0.0048 a-1；春季增幅最大，近15年增長了11.75%；秋季增長最慢，增長趨勢和變化幅度分別為0.0019 a-1和5.98%。

不同植被類型的NDVI在各生長季均體現出增加趨勢，但其增長趨勢及幅度差異較大。其中，夏季林地的NDVI增長最快(0.0052 a-1)，春季耕地NDVI增幅最大(19.99%)。研究區春、秋季的NDVI增長趨勢表明，植被生長周期的開始日期在提前、結束日期在推后，植被生長季在延長。這與“在全球變暖背景下植被生長季在延長”的相關結論一致[21,22]。

2.4 NDVI對氣溫和降水變化的響應

為探討研究區NDVI變化對氣候變化的響應關系，逐像元計算了NDVI與氣溫、降水之間的相關系數，以構建研究區氣候變化與植被覆蓋之間的關系(表2、表3)。

表2 不同類型植被NDVI與降水的相關系數Table 2 Correlation coefficients between seasonal mean NDVI and precipitation for different vegetation types

注：*表示通過0.05的置信度檢驗。

表3 不同類型植被NDVI與氣溫的相關系數Table 3 Correlation coefficients between seasonal mean NDVI and temperature for different vegetation types

注：**和*分別表示通過0.01和0.05的置信度檢驗。

從整個研究區生長季的NDVI看，其與降水呈明顯正相關(R=0.472,P=0.112)，而與氣溫的相關性較弱(R=0.054,P=0.857)；不同植被類型均體現出相似的特征。這說明本研究區，降水對植被生長季NDVI的影響要大于氣溫，是影響平均生長季NDVI的主要因素(表2、表3)。

無論是從生長季還是不同季節看，三種植被類型的NDVI與降水的相關系數基本為耕地>草地>林地，說明降水對林地、草地和耕地的植被生長的促進作用依次增大。這可能是由于耕地集中在海拔相對低的地區，降水量較少，氣溫較高，蒸發量較大，因而降水是影響其植被狀況的主要因素，致使NDVI與降水的相關性大。而林地主要分布在海拔相對高的地區，降水量豐富，氣溫較低，蒸發量較小，降水量對植被狀況的影響相對較小，因而林地的NDVI與降水量的相關性小。各植被類型NDVI與降水的相關性在夏季較大(R均大于0.45)，夏季降水是促進植被生長的重要因素；而春、秋季NDVI與降水呈弱相關性(表2)。

雖然生長季NDVI與氣溫之間相關性較弱，但是研究區春季NDVI與氣溫呈現強正相關性，草地、耕地和林地與氣溫的相關系數分別為0.532、0.575(P<0.01)和0.596(P<0.01)，說明春季氣溫對植被的生長具有明顯的促進作用。春季氣溫的上升可能是導致研究區春季平均NDVI上升的主要原因，是植被對全球變暖響應的重要信號。夏季氣溫與三種植被的夏季NDVI均呈負相關，體現了夏季氣溫對三種植被夏季生長的抑制作用，這種抑制作用由大至小依次為耕地、林地和草地。夏季氣溫與植被NDVI的負相關性，主要是由于研究區夏季氣溫的升高加速地表水分的蒸發，從而抑制植被的生長(表3)。

3 結論

本文的主要研究結論如下：1)研究區NDVI的空間分布呈現出明顯的垂直地帶性。從東北向西南隨著海拔高度的增加NDVI的值先增大后減小，植被呈東南-西北方向帶狀分布。2)近15年，研究區大部分地區植被活動在增強，總體年均NDVI與生長季NDVI均呈增加趨勢，年均增長量分別為0.002和0.0025。生長季NDVI增強的區域占研究區面積的81.3%，主要分布在北部區域。不同類型植被的生長季NDVI以及春、夏、秋季的NDVI都在上升，植被生長季在延長。3)研究區年均NDVI、生長季NDVI變化都與降水量有較大的相關性，降水量是影響植被狀況的主要因素；氣溫對植被的影響主要體現在春季的促進作用和夏季的抑制作用，春季氣溫的上升可能是導致研究區春季平均NDVI上升的主要原因，是植被對全球變暖響應的重要信號；不同類型植被體現出對降水和氣溫不同程度的敏感性，降水對三種植被類型的促進作用為耕地>草地>林地，春季氣溫的促進作用為林地>耕地>草地，夏季氣溫的抑制作用為耕地>林地>草地。

[1] PIAO S,MOHAMMAT A,FANG J,et al.NDVI-based increase in growth of temperate grasslands and its responses to climate changes in China[J].Global Environmental Change,2006,16(4):340-348.

[2] SUN H,WANG C,NIU Z.Analysis of the vegetation cover change and the relationship between NDVI and environmental factors by using NOAA time series data[J].Journal of Remote Sensing,1998,2:210-216.

[3] 馬金珠,王剛,李丁,等.民勤盆地沙漠化,貧困化與水資源開發利用[J].自然資源學報,2006,21(4):551-558.

[4] 任建民,忤彥卿,貢力.人類活動對內陸河石羊河流域水資源轉化的影響[J].干旱區資源與環境,2007,21(8):7-11.

[5] 施雅風,沈永平,胡汝驥.西北氣候由暖干向暖濕轉型的信號、影響和前景初步探討[J].冰川凍土,2002,24(3):219-226.

[6] SHI Y,SHEN Y,KANG E,et al.Recent and future climate change in northwest China[J].Climatic Change,2007,80(3-4):379-393.

[7] PIAO S,FANG J,ZHOU L,et al.Interannual variations of monthly and seasonal normalized difference vegetation index(NDVI)in China from 1982 to 1999[J].Journal of Geophysical Research:Atmospheres(1984-2012),2003,108(D14):4401-4413.

[8] 宋怡,馬明國.基于 SPOT VEGETATION 數據的中國西北植被覆蓋變化分析[J].中國沙漠,2007,27(1):89-93.

[9] DAI S,ZHANG B,WANG H,et al.Vegetation cover change and the driving factors over northwest China[J].Journal of Arid Land,2011,3(1):25-33.

[10] 李秀花,師慶東,常順利,等.中國西北干旱區 1981-2001 年 NDVI 對氣候變化的響應分析[J].中國沙漠,2009,29(5):989-994.

[11] PIAO S,FANG J,JI W,et al.Variation in a satellite-based vegetation index in relation to climate in China[J].Journal of Vegetation Science,2004,15(2):219-226.

[12] ZHAO X,TAN K,ZHAO S,et al.Changing climate affects vegetation growth in the arid region of the northwestern China[J].Journal of Arid Environments,2011,75(10):946-952.

[13] MAISONGRANDE P,DUCHEMIN B,DEDIEU G.VEGETATION/SPOT:An operational mission for the Earth monitoring;presentation of new standard products[J].International Journal of Remote Sensing,2004,25(1):9-14.

[14] HOLBEN B N.Characteristics of maximum-value composite images from temporal AVHRR data[J].International Journal of Remote Sensing,1986,7(11):1417-1434.

[15] VERBESSELT J,SOMERS B,VAN AARDT J,et al.Monitoring herbaceous biomass and water content with SPOT VEGETATION time-series to improve fire risk assessment in savanna ecosystems[J].Remote Sensing of Environment,2006,101(3):399-414.

[16] CAO X,CHEN X,BAO A,et al.Response of vegetation to temperature and precipitation in Xinjiang during the period of 1998-2009[J].Journal of Arid Land,2011,3(2):94-103.

[17] FAN J,HUANG J,ZHANG M.Retrieval of cropping index in China using time series of SPOT vegetation NDVI[J].Sensor Letters,2013,11(6-7):1134-1140.

[18] HOU X,GAO S,NIU Z,et al.Extracting grassland vegetation phenology in North China based on cumulative SPOT-VEGETATION NDVI data[J].International Journal of Remote Sensing,2014,35(9):3316-3330.

[19] STOW D,DAESCHNER S,HOPE A,et al.Variability of the seasonally integrated normalized difference vegetation index across the north slope of Alaska in the 1990s[J].International Journal of Remote Sensing,2003,24(5):1111-1117.

[20] 田義超,梁銘忠.基于 SPOT-VEGETATION 的北部灣經濟區植被變化特征[J].農業現代化研究,2014(4):465-471.

[21] PIAO S,NAN H,HUNTINGFORD C,et al.Evidence for a weakening relationship between interannual temperature variability and northern vegetation activity[J].Nat Commun,2014,5.doi:10.1038/ncomms6018.

[22] 王宏,李曉兵,李霞,等.基于 NOAA NDVI 和 MSAVI 研究中國北方植被生長季變化[J].生態學報,2007,27(2):504-515.

Spatiotemporal Changes of Vegetation and the Responses to Temperature and Precipitation in Upper Shiyang River Basin

TANG Zhi-guang1,MA Jin-hui2,3,LI Chao-kui1,PENG Huan-hua1,LIANG Ji1

(1.National-LocalJointEngineeringLaboratoryofGeo-spatialInformationTechnology,HunanUniversityofScienceandTechnology,Xiangtan411201;2.KeyLaboratoryofWestChina′sEnvironmentalSystem,MinistryofEducation,LanzhouUniversity,Lanzhou730000;3.InstituteofRSandGeographicalInformationSystem,LanzhouUniversity,Lanzhou730000,China)

The Shiyang River Basin is a typical arid inland river basin in northwestern China.It is also a well-known fragile area of water resource and ecological environment,and a sensitive area of climate change.The vegetation in Shiyang River Basin is mostly concentrated in its upstream,and plays an important role in ecological environment of this watershed.Therefore,studies on the spatiotemporal variation of vegetation in the upper Shiyang River Basin and the impacts of changing temperature and precipitation on the vegetation growth are of significance to dealing with climate change,realizing sustainable development of ecological environment in Shiyang River Basin.In this study,the spatiotemporal changes of vegetation growth in upper Shiyang River Basin together with their responses to climate changes were investigated using SPOT-VEGETATION Normalized Difference Vegetation Index(NDVI)and climate datasets from 1999 to 2013.Results reveal that there is an obvious vertical zonal characteristic of the spatial distribution of NDVI.About 81.3% of the study area shows an increasing trend in growing season NDVI.The average NDVI of the study area increases at rates of 7.75% for the growing season(March-November),11.75% for spring(March-May),9.62% for summer(June-August),and 5.98% for autumn(September-November)over the study period.The increase of NDVI in spring and autumn suggests the growing season of the vegetation in this study area has been prolonged.The effects of climate changes on vegetation growth vary with the types of vegetation and seasons,which shows a large spatial and temporal heterogeneity.As compared with temperature,precipitation is the dominant climatic factor affecting the interannual variations of vegetation.

Shiyang River Basin;vegetation;NDVI;spatiotemporal changes;climate change

2015-11-13；

2016-02-26

國家自然科學基金項目(41501070、31400409、41101319);湖南科技大學校級科研項目(E51520)

唐志光(1985-)，男，博士，講師，主要從事環境遙感與GIS應用方面的研究。E-mail：tangzhg11@lzb.ac.cn

10.3969/j.issn.1672-0504.2016.03.021

P935.1;P46

A

1672-0504(2016)03-0116-05