近15年黃土高原植被物候時空變化特征分析

李 強，張 翀，任志遠

?

近15年黃土高原植被物候時空變化特征分析

李 強1,2，張 翀3，任志遠1

（1陜西師范大學旅游與環境學院，西安 710119；2陜西學前師范學院環境與資源管理系，西安 710100；3寶雞文理學院陜西省災害監測與 機理模擬重點實驗室，陜西寶雞 721013）

【目的】黃土高原處于從濕潤向干旱過渡、從森林向草原過渡、從農業向牧業過渡的地區，是中國氣候變化與農業發展的敏感地帶，針對此區域的地表植被覆蓋物候特征研究，對于該地區農業生產、環境保護和生態建設具有重要的指示意義。通過分析不同時間序列、不同海拔高度與水熱條件下植被物候趨勢的差異，以期為黃土高原當前的農業生態環境建設與可持續發展提供有用的理論支撐與決策依據。【方法】基于1998—2012年SPOT VEGETATION旬值NDVI數據，并結合諧波分析法、線性趨勢等方法對黃土高原各年植被物候特征值進行了確定，并分析了物候的變化趨勢。【結果】（1）1998—2012年，生長季始期平均每年提前0.9 d。生長季末期平均每年推遲約0.8 d。在生長季始期提前和末期推遲的作用下，生長季長度平均每年延長1.7 d。（2）黃土高原水熱組合直接影響植被物候的空間差異性，植被生長的限制性氣溫為9℃，限制性降水分別為475 mm與540 mm，限制性海拔為1 750 m。（3）植被生長季長度趨勢與海拔和氣溫的空間偏相關系數分別為0.0591和0.0139，與降水的空間偏相關系數為-0.0174，三要素與生長季始期趨勢的相關程度較與生長季末期趨勢強。【結論】黃土高原植被物候特征趨勢明顯且穩定的區域主要分布在陜北高原與晉中北山地區。西北部干旱區荒漠草原區，物候變化主要受氣溫控制。半干旱地區農業與草原區，物候變化主要受到降水量控制。汾渭盆地農業區，物候變化受水熱共同作用。水熱的差異對秦巴山地闊葉林區植被物候的影響不明顯。海拔對黃土高原植被物候變化趨勢上的影響不明顯。生長季長度趨勢在生長季始期和生長季末期的共同作用下隨海拔和氣溫的增加的延長趨勢增強，隨降水增加的縮短趨勢增強。同種植被物候趨勢隨海拔、降水、氣溫的變化特征具有一致性，生長季始期的變化特征對于植被生長季長度變化影響較生長季末期強。

植被覆蓋；生長季；諧波分析；傅里葉插值；黃土高原

0 引言

【研究意義】黃土高原農業生態環境脆弱，水旱災害、氣象災害以及水土流失等農業災害比較頻繁和嚴重。植被作為陸地生態圈的重要組成，是氣候系統的重要元素[1]，對區域生態環境與自然災害有重要影響，對農業生產、環境保護和生態建設具有重要意義[2]。在這樣的背景下，掌握研究區地表植被覆蓋度年際間的變化規律與空間差異性特征，對評價區域農業生態系統的環境質量、調節生態過程、促進農業可持續發展具有重要的理論和實際意義[3]。【前人研究進展】歸一化植被指數（normalized difference vegetation index，NDVI）由紅波段與近紅外波段的反射率計算而來[4]，和植物的生產力密切相關[5]，并且NDVI趨勢可以用來衡量植被覆蓋的改善與退化[6]。NDVI趨勢被用在很多測算中，包括農業發展的生態響應[7]、物候變化[8]、作物狀況[9]、耕地利用[10]以及土壤肥力變化[11]。一系列研究表明在全球，特別是北半球，存在生長季的開始時刻提前以及生長季長度增大的趨勢[12]，生長季長度增大與氣溫的上升會加快地表水分蒸發，增加干旱脅迫與林火發生的可能[13]，并增加了固碳強度[14]，進而大范圍長時間地影響作物生產與農業發展。目前，許多學者對黃土高原植被物候變化特征進行了分析，謝寶妮等[15]分析認為黃土高原植被物候主要受氣溫影響，冬季和前年秋季氣溫上升是春季物候提前的主要驅動因子。韋振鋒等[16]認為陜甘寧多年物候同時與緯度和降水有關。張晗等[17]研究認為隨著緯度的升高陜西省植被返青期逐步推遲，枯黃期逐步提前。【本研究切入點】前人的研究均側重于單一要素對植被物候特征的影響，而植被生長并不是單個要素作用的結果，需要綜合分析多種因素的影響。本研究基于1998—2012年SPOT VEGETATION旬值NDVI數據，并結合GIS空間分析方法，將植被物候變化的生態過程與空間格局相耦合，對黃土高原近15年植被物候特征值進行了確定，并嘗試分析其時空差異性，進而深入探討海拔、氣溫和降水3種因素分別對植被物候時空變化的影響機制。【擬解決的關鍵問題】采用時間序列諧波分析法和線性趨勢等方法對黃土高原植被生長季時刻進行了確定，計算了黃土高原生長季長度及其趨勢，并進一步解釋不同海拔高度與水熱條件下植被物候趨勢的差異。

1 材料與方法

1.1 研究區概況

黃土高原西起祁連山支脈烏鞘嶺，東迄太行山，北到長城，南抵秦嶺。橫跨青、甘、寧、內蒙古、陜、晉、豫7省區，面積約63萬km2，約占全國陸地總面積的8%。是東南濕潤季風氣候向西北內陸干旱氣候過渡、暖溫帶落闊葉林向典型草原和荒漠草原過渡的過渡地帶。受溫帶大陸性氣候的控制，氣溫、降水量季相分明；從東南向西北，氣候依次為半濕潤氣候、半干旱氣候和干旱氣候；土壤依次為褐土、壚土、黃綿土和灰鈣土；植被依次出現森林草原、典型草原和荒漠草原。本區水土流失十分嚴重，是黃河泥沙的主要來源區。

1.2 數據來源

數據包括黃土高原1998—2012年SPOT VEGETATION旬值NDVI數據、131個臺站1998—2012年的年降水資料和年平均氣溫資料以及DEM數據。NDVI數據是由比利時佛萊芒技術研究所（Flemish Institute for Technological Research，Vito）VEGETATION影像處理中心（VEGETATION processing Centre，CTIV）負責預處理并提供免費下載，空間分辨率為1 000 m；降水和氣溫數據均來自于中國氣象數據網，采用經驗貝葉斯克里金法將其插值到空間上得到15年的年降水和年均溫柵格數據，插值柵格分辨率為1 000 m；數字高程模型為GDEMV2 30M分辨率數字高程模型數據（Digital Elevation Model，DEM），來源于地理空間云數據。

1.3 時間序列的諧波分析與物候特征提取

時間序列諧波分析法（harmonic analysis of time series，HANTS）是平滑和濾波兩種方法的綜合，它能夠充分利用遙感圖像存在時間性和空間性的特點，將其空間上的分布規律和時間上的變化規律聯系起來。首先，應用HANTS方法對原始的NDVI旬數據進行逐年平滑處理，然后利用傅里葉插值將平滑后的逐年NDVI旬數據插值成逐年NDVI天數據。

對平滑后的逐年NDVI天數據，計算植被生長季的逐年始期、末期與長度，以及多年平均始期、末期與長度。文中采用較為常用的最大比率法[18-19]，最大比率即為植被覆蓋最大變化速率。首先，計算逐年NDVI天數據相鄰時間點上的差分值（NDVI-NDVI，其中），差分可以反映出相鄰時間點的比率，利用差分值的最大值對應天數確定為生長季始期（start of growing season，SOG），最小值對應天數確定為生長季末期（end of growing season，EOG），EOG和SOG的之差為生長季長度（length of growing season，LOG）（圖1）。

圖1 NDVI 天數據最大相對變化值來確定物候特征

1.4 分析方法

文中采用Theil-Sen趨勢對植被物候時空特征進行趨勢分析，Theil-Sen趨勢既不受異常值影響，也不用服從一定的分布，結果更為科學、可信[20]，以剔除異常年份對物候特征造成的影響。計算公式如下：

其中，為趨勢，x，x為時間序列數據，t，t為時間序列對應的年份。

由于植被物候時空分布是多種因素作用的結果，文中分析了高程、降水和氣溫的影響，所以文中將3種要素進行等間距劃分，并分別統計出每種間隔內植被物候特征的均值，以反映不同高程帶、降水帶和氣溫帶上的時空變化差異。

2 結果

2.1 多年平均物候特征

2.1.1 多年平均物候空間分布 黃土高原植被物候多年均值空間分布見圖2。由黃土高原東南向西北方向，SOG逐漸推遲，EOG逐漸提前，LOG逐漸縮短，反映出氣候影響下的空間分異規律。黃土高原東南部，SOG主要集中在80—120 d之間，而在關中平原、汾河谷地與黃淮平原，SOG早于80 d，西北部則晚于150 d；黃土高原東南部，EOG主要晚于310 d，也表現出海拔上的差異，黃土高原西北部主要分布在280—300 d，后套平原、寧夏平原以及青海境內在260—270 d，晉中北、鄂爾多斯高原中部、隴中高原以北等山地區以及前套平原在270—280 d；黃土高原東南部，LOG主要大于210 d，而在西北部的平原區以及青海境內則小于120 d。

2.1.2 多年平均物候的氣候分異特征 黃土高原植被物候與降水量關系十分密切（圖3-(a-c)）。降水量每增加100 mm，SOG提前13 d，EOG推遲3 d，LOG增大16 d。降水量在200—475 mm，物候隨降水量變化不明顯；降水量高于475 mm以上區域，物候隨降水量變化的波動性增強，尤其是EOG；降水量在475—670 mm，隨降水量增大，SOG提前最為明顯，而EOG波動性推遲，LOG波動性增大；降水量在670 mm以上的區域，SOG開始推遲，EOG先推遲再提前，LOG減小。降水量在670 mm以上的區域主要分布在黃土高原東南部，植被覆蓋較高，且人口密集，主要分布為種植作物，植被物候與降水量的關系不明顯可能受人類活動影響導致，比如一年兩熟作物的種植導致物候特征提取精度降低。

圖3-(d-f)反映出，黃土高原氣溫與降水量相比，隨著氣溫的升高，EOG推遲以及LOG增大強度均高于隨降水量變化強度，而SOG相對較小。黃土高原氣溫每升高2℃，SOG提前12 d，EOG推遲9 d，LOG增大21 d。可以看出在氣溫低于6℃的地區，物候隨氣溫變化不明顯。隨著氣溫升高，EOG推遲明顯，其次是LOG的增大，而氣溫在6—9℃之間SOG甚至表現為推遲，高于9℃，SOG明顯提前，LOG增大明顯。

圖3 黃土高原植被物候多年均值與水熱之間的關系

圖4顯示的是每1 mm降水量帶，每0.1℃氣溫帶上物候多年均值變化。可以看出降水量在475 mm與540 mm處，以及氣溫在9℃處，SOG、EOG和LOG存在明顯變化。在降水量小于475 mm的區域，隨著降水量的增加，物候變化不明顯；隨著氣溫的升高，SOG逐漸推遲，EOG逐漸推遲，LOG變化不明顯。主要是因為該區域屬于干旱半干旱地區，降水量在小于475 mm的范圍內變化，植被生長不明顯，而氣溫的上升會加快地表水分蒸發，從而抑制生長初期植被生長并延緩了生長末期時間。降水量在475—540 mm之間的區域，當氣溫低于9℃，物候變化不明顯；當氣溫高于9℃，SOG逐漸提前，EOG逐漸推遲，LOG逐漸增大。該區域降水量大于475 mm，滿足植被生長需求，并且氣溫高于9℃時，便會促進植被生長。在降水量大于540 mm的區域，隨著降水量的增加，除降水量在640—690 mm之間的區域SOG有所提前外，SOG變化不明顯，EOG與LOG波動性較強，規律性較差；隨著氣溫的變化，SOG變化不明顯，EOG與LOG也呈較強波動性，主要是由于該地區農作物生長季末期受人類活動（收割與耕種等）影響造成的。但是在降水量大于700 mm的地區，EOG與LOG隨著氣溫的變化不明顯，隨降水量變化波動性加強，該區域主要集中在秦嶺山地區，植被覆蓋高，土壤含水量豐富，降水量和氣溫的變化對植被生長作用不明顯。

2.1.3 多年平均物候的海拔分異特征 圖5中海拔每升高1 000 m，SOG推遲13 d，EOG提前9 d，LOG縮短23 d。在海拔1 750 m以下地區，物候隨海拔變化波動比較強，隨著海拔升高，SOG主要變現為推遲，EOG變化不明顯，LOG主要表現為減小；在海拔1 750 m以上地區，隨海拔升高，SOG不斷推遲，EOG在1 750—2 350 m區域明顯提前，大于2 350 m則變化不明顯，LOG不斷減小。海拔大于1 750 m的地區主要分布在隴中高原與青海境內，另外晉北山地區也有分布，這些地區海拔高，氣溫低，并且降水量相對較少，主要以自然植被為主；而海拔小于1 750 m的地區地形復雜，降水與氣溫空間差異較大，且容易受到人類活動影響。1 750 m為黃土高原植被生長隨海拔明顯分異的界限，高于1 750m，植被生長相對緩慢，低于1 750 m則植被生長長度延長。

2.2 物候年際變化特征

2.2.1 物候區域尺度上的年際變化 黃土高原植被LOG變化趨勢最為明顯，其次為SOG，而EOG趨勢最弱。SOG呈提前趨勢，平均每年平均提前0.9 d。EOG呈推遲趨勢，但趨勢相對較弱，平均每年推遲約0.8 d。LOG變化趨勢較強，呈上升趨勢，平均每年的變化幅度為1.7 d。SOG趨勢高于全球的0.3 d·a-1、歐亞大陸的0.33 d·a-1、北半球的0.31 d·a-1、中國的0.79 d·a-1，以及北半球的0.02 d·a-1。EOG的推遲趨勢均高于全球的0.05 d·a-1、歐亞大陸的0.61 d·a-1、北半球的0.23 d·a-1和中國的0.37 d·a-1，以及北半球的 0.25 d·a-1。LOG的增加趨勢遠遠高于全球的0.38 d·a-1，高于北半球的0.56 d·a-1、歐亞大陸的1.33 d·a-1和中國的1.16 d·a-1，以及北半球的0.28 d·a-1[21-25]。

2.2.2 物候年際變化的空間分布特征 采用Sen趨勢計算了黃土高原植被物候年際趨勢空間差異，如圖6所示。73.71% 的區域SOG呈提前趨勢，26.29% 的區域有推遲趨勢，其中提前的區域占24.27%，主要分布在陜北高原與晉中北山地區，其次隴中高原以西也有分布。

38.73%的區域EOG具有提前趨勢，61.27%的區域具有推遲趨勢，其中提前的區域約為4.34%，推遲的區域約為16.27%。從圖6-b可以看出，大部分地區EOG呈推遲趨勢，推遲趨勢集中在0—2.0 d·a-1，推遲的區域主要分布于高原西北部邊緣和南部。

從LOG變化來看，80.77%的區域具有增長趨勢，19.23%的區域具有縮短趨勢，其中增長的區域約為30.06%，縮短的區域約為4.38%，增長與縮短的空間分布與SOG提前與推遲的分布極為相似，說明黃土高原生長季長度主要由SOG的變化引起。

圖5 黃土高原植被物候多年均值與海拔之間的關系

圖6 黃土高原植被物候變化趨勢的空間分布

2.2.3 生長季物候年際變化的海拔分異特征 從SOG 來看（圖7），850—4 000 m的地區呈提前趨勢；小于850 m的區域呈推遲趨勢；大于4 000 m的區域，表現為推遲趨勢，但顯著性隨海拔升高逐漸減小。從EOG來看，低于2 500 m的地區，呈推遲趨勢，顯著性隨海拔減小，但在1 350—2 500 m顯著性逐漸回升；高于2 500 m的地區，呈提前趨勢，但是顯著性隨海拔升高在逐漸減小。LOG主要呈增長趨勢，顯著性隨海拔升高整體上是減小的，在小于850 m或大于4 000 m的地區表現為縮短趨勢。

黃土高原海拔小于850 m的地區為汾渭盆地農業區，占總面積的11.41%；海拔在850—2 500 m之間的區域，植被物候變化趨勢相對最為明顯，占總面積的82.93%；海拔在2 500—4 000 m之間包括祁連山森林與高寒草原區以及江河源區-甘南高寒草甸草原區，面積占5.34%；大于4 000 m的區域占0.33%。所以海拔對黃土高原植被物候變化趨勢上的影響相對于水熱條件不明顯。

2.2.4 生長季物候年際趨勢的水熱條件的分異特征 從SOG與降水量的關系來看（圖8），隨降水量變化SOG基本上呈提前趨勢，在260—560 mm的區域，隨降水量的增大，提前幅度有所減小，由1.3 d·a-1減小到0.4 d·a-1，但是顯著性在升高；560—720 mm的地區，推遲趨勢明顯；大于720 mm的地區，提前幅度增強，由0增加到1.0 d·a-1，而顯著水平先減后增。從EOG 來看，絕大部分研究區呈推遲趨勢，僅在630 m左右呈提前趨勢；推遲幅度隨降水量增大在緩慢升高，顯著水平先增后減。LOG主要呈增長趨勢，在350—560 mm地區上增長，在600—660 mm的地區表現為減小趨勢。

從SOG與氣溫的關系來看（圖9），氣溫在小于8℃地區上提前幅度緩慢增大，8—12℃區域上呈提前趨勢。從EOG來看，氣溫在大于5℃地區上，隨氣溫升高表現為推遲趨勢，趨勢幅度與顯著性均在波動性增加。LOG主要呈增長趨勢，在小于10℃地區上，增長幅度與顯著性隨氣溫升高波動性增大，而在大于10℃地區上，幅度與顯著性在不斷減小。

圖7 黃土高原植被物候年際變化特征在海拔上的分異特征（a：SOG；b：EOG；c：LOG）

圖8 黃土高原植被物候年際變化特征在降水上的分異特征（a：SOG；b：EOG；c：LOG）

圖9 黃土高原植被物候年際變化特征在氣溫上的分異特征（a：SOG；b：EOG；c：LOG）

2.3 植被物候年際趨勢的海拔、水熱空間分異特征

由于植被物候年際趨勢受到降水、氣溫、海拔的綜合影響，上述分異特征中無法剔除其他兩種要素的影響，所以對像元某一鄰域內的所有柵格進行偏相關分析，其中鄰域內所有柵格的物候特征年際趨勢作為因變量，領域內三要素作為自變量，剔除其他兩種要素得到物候趨勢與某一要素的偏相關系數，以反映植被物候在3種影響因素上的空間分異特征（圖10—12）。

圖10 黃土高原植被SOG與海拔、降水、氣溫的空間偏相關系數

圖11 黃土高原植被EOG與海拔、降水、氣溫的空間偏相關系數

圖12 黃土高原植被LOG與海拔、降水、氣溫的空間偏相關系數

整體上，黃土高原植被物候年際趨勢隨DEM、降水和氣溫表現出明顯的空間分異特征。其中SOG隨著海拔增加呈提前趨勢，隨降水量增加呈推遲趨勢，隨氣溫增加呈提前趨勢；EOG趨勢空間分異特征正好與SOG相反；LOG在SOG和EOG的共同作用下隨DEM和氣溫的增加呈延長趨勢，隨降水增加呈縮短趨勢（表1）。

表1 黃土高原植被物候趨勢與影響因素的空間偏相關系數

空間上，植被物候趨勢隨DEM的空間分異特征明顯，隨降水和氣溫的變化主要體現在較小區域范圍，大范圍上差異性不明顯（圖10—12），植被類型的空間差異性產生的影響不容忽視。因此，利用植被類型對空間偏相關系數進行統計，并基于空間偏相關系數利用K-means法對各個植被類型進行聚類，如表2所示。

從聚類結果可以看出，鹽地落葉灌叢、荒漠化草原、強旱生小灌木-半灌木-灌木與沙生植被歸屬為第1類，該類植被SOG隨DEM增加呈提前趨勢，EOG隨DEM增加呈推遲趨勢，LOG隨DEM增加呈延長趨勢，溫室氣體的增加以及全球變暖會致使各地區植被物候始期提前和末期推遲，從而導致生長期延長，隨著氣溫的升高，高海拔植被感溫效果較為明顯，因此各種植被物候趨勢隨海拔變化具有較強的一致性；其中，鹽地落葉灌叢、荒漠化草原與沙生植被SOG隨降水增加呈提前趨勢，強旱生灌木SOG卻隨降水增加呈推遲趨勢；各個植被的EOG趨勢與降水量的空間相關性較SOG趨勢與降水量的相關性強烈，僅沙生植被EOG隨降水量增加呈推遲趨勢，因此其LOG隨降水增加呈延長趨勢，其他3種植被EOG均在降水量增加下不斷提前，從而導致其LOG趨勢隨降水增加而縮短；該類植被LOG隨氣溫上升呈延長趨勢，其中荒漠化草原、強旱生灌木與沙生植被主要是SOG隨氣溫上升的提前趨勢引起的，而鹽地落葉灌叢是由EOG隨氣溫上升的推遲趨勢引起的。

竹林、強旱生雜類草與沼澤植被各自成一類。竹林對水熱條件要求均較高，其SOG隨DEM與降水的增加均呈提前趨勢，EOG均呈推遲趨勢，從而竹林LOG隨DEM與降水的增加呈延長趨勢；氣溫的上升雖然使EOG推遲趨勢增強，但由于SOG的推遲更甚，所以竹林LOG隨氣溫上升呈縮短趨勢；強旱生雜類草LOG隨DEM增加呈延長趨勢，由SOG的提前趨勢引起，LOG隨降水和氣溫的增加呈縮短趨勢，均由EOG的提前趨勢引起；沼澤植被LOG隨降水增加呈縮短趨勢，主要由EOG提前趨勢引起，LOG均隨DEM和氣溫增加呈延長趨勢，由SOG提前與EOG推遲共同決定。

寒溫性針葉林、高寒灌叢、高寒草原與（亞）高山草甸物候變化特征類似。除過（亞）高山草甸LOG隨氣溫升高呈延長趨勢外，其他植被LOG均隨DEM、降水與氣溫增加呈縮短趨勢，主要是由于SOG的推遲趨勢和EOG的提前趨勢共同作用的結果。

其他植被類型雖被聚為一類，但是其物候隨三因素的空間變化特征差異性較大，經過對其進行重新聚類，聚為7類，說明這些植被類型物候變化趨勢受三因素影響的空間差異性較大。其中，落葉闊葉林與溫性落葉灌叢物候趨勢特征相似，其LOG均隨DEM增加呈縮短趨勢，主要由于SOG推遲趨勢導致；LOG均隨氣溫上升呈延長趨勢，由SOG的提前引起；LOG受降水影響特征不一致。雜木林、草叢與草本栽培植被物候隨三要素變化特征高度一致，其LOG隨DEM與氣溫增加均呈延長趨勢，是SOG提前趨勢與EOG推遲趨勢共同作用的結果；LOG隨降水增加呈縮短趨勢，主要為SOG推遲趨勢造成的，另外還有EOG提前趨勢的作用。典型草原與河-湖灘地草甸LOG均隨三要素的增加呈延長趨勢，且SOG與EOG均隨三要素呈提前趨勢，而SOG提前趨勢更強。綜上，植被物候趨勢隨DEM、降水、氣溫的變化表現出明顯的自然分帶特征，SOG的變化特征對于植被LOG變化影響較EOG強烈。

表2 黃土高原各個植被類型物候趨勢與影響因素的空間偏相關系數及其聚類結果

SOG-DEM為不同植被類型SOG與DEM空間偏相關系數的均值，其他以此類推；聚類結果中1、2等表示大類，31、32等表示第3類重聚類的結果，個位數為重聚類類別

SOG-DEM is the spatial mean coefficient of partial correlation of SOG & DEM in different vegetation types, and so on. 1, 2, etc are big class. 31, 32, etc are re-cluster results, where single digits are re-cluster classes

3 討論

黃土高原植被物候特征變化穩定且趨勢明顯的區域主要分布在陜北高原與晉中北山地區，表現為生長季長度增長。該地區屬半干旱地區，降水量相對較少，加之全球處于升溫階段，生長季長度的增加必然會加劇土壤水分的減小，近15年來，該區域生長季長度的增長主要是由人類活動的影響造成的，三北防護林建設以及退耕還林還草政策的實施，人工林草的種植在很大程度上改變了地表覆蓋類型，并增加了該地區植被覆蓋，從而改變了黃土高原植被的物候特征。

目前，植被物候對氣候、地形等要素的響應受到了國內外學者的廣泛關注與研究[15, 26-27]。然而在許多影響因素中存在著信息的冗余，影響因素之間相互影響、相互依存、互為因果，且植被物候是多種因素綜合作用的結果，所以如何剔除其他影響因素的影響，分析植被物候與某因素之間的關系至關重要。本文選擇鄰域計算方法，在鄰域內利用偏相關系數表達植被物候趨勢與某因素的相關程度，既剔除了其他因素影響，也保證了結果是基于像元計算的，并統計出不同植被類型的偏相關系數，基于空間偏相關系數利用K-means法對各個植被類型進行聚類，發現植被物候趨勢隨DEM、降水、氣溫的變化表現出明顯的自然分帶特征，而同一地域不同植被物候變化的驅動機制差異性較大，因此不同類型植被物候變化的驅動力分析值得關注。同時還需關注年內各因子對植被物候的影響。

另外，本文僅分析了黃土高原植被物候趨勢在水熱條件與海拔高度上的分異特征，并沒有分析這種趨勢的產生機制，而不同地域物候的變化應該主要受當地氣候要素的變化影響，比如，降水量滿足植被生長的區域，如果氣溫處于增長趨勢，則可能導致生長季始期的提前與末期的推遲，從而導致生長季長度的增長。由于篇幅限制，本文未做更進一步探討，筆者在以后的研究中將對這方面進行更深層次的分析。

4 結論

黃土高原植被物候多年均值由東南向西北，生長季始期（SOG）逐漸推遲，生長季末期（EOG）逐漸提前，生長季長度（LOG）逐漸縮短，反映出氣候的空間分異規律。1998—2012年間SOG呈提前趨勢，平均每年提前0.9d，EOG平均每年推遲約0.8d，LOG在SOG和EOG共同作用下變化率達1.7 d·a-1。黃土高原生長季長度變化主要由SOG的變化引起，物候特征變化趨勢明顯且穩定的區域主要分布在陜北高原與晉中北山地區，其次隴中高原以西也有分布。

黃土高原植被物候年際趨勢隨DEM、降水和氣溫表現出明顯的空間分異特征。其中SOG趨勢隨海拔增加的提前趨勢增強，隨降水量增加的推遲趨勢增強，隨氣溫增加的提前趨勢增強；EOG趨勢空間分異特征正好與SOG相反；LOG在SOG和EOG的共同作用下隨DEM和氣溫的增加的延長趨勢增強，隨降水增加的縮短趨勢增強。同種植被物候趨勢隨DEM、降水、氣溫的變化特征具有一致性，SOG的變化特征對于植被LOG變化影響較EOG強烈。

References

[1] Foley J A, Levis S, Costa M H, Cramer W, Pollard D. Incorporating dynamic vegetation cover within global climate models., 2000, 10: 1620-1632.

[2] IPCC. 4th Assessment report of the intergovernmental panel on climate change., 2007: 52.

[3] 穆少杰, 李建龍, 陳奕兆. 2001-2010年內蒙古植被覆蓋度時空變化特征. 地理學報, 2012, 67(9): 1255-1268.

MU S J, LI J L, CHEN Y Z. Spatial differences of variations of vegetation coverage in Inner Mongolia during 2001-2010., 2012, 67(9): 1255-1268. (in Chinese)

[4] Tucker C J. Red and photographic infrared linear combinations for monitoring vegetation., 1979, 8(2): 127-150.

[5] Prince S D, Tucker C J. Satellite remote sensing of rangelands in Botswana II. NOAA AVHRR and herbaceous vegetation., 1986, 7(11): 1555-1570.

[6] Alcaraz-Segura D, Chuvieco E, Epstein H E, Kasischke E S, Trishchenko A. Debating the greening vs. browning of the North American boreal forest: Differences between satellite datasets., 2009, 16(2): 760-770.

[7] Pettorelli N, Vik J O, Mysterud A, Gaillard J M, Tucker C J, Stenseth N C. Using the satellite-derived NDVI to assess ecological responses to environmental change., 2005, 20(9): 503-510.

[8] White M A, de Beurs K M, Didan K, Inouye D W, Richardson A D, Jensen O P, O’KEEFE J, ZHANG G. Intercomparison, interpretation, and assessment of spring phenology in North America estimated from remote sensing for 1982-2006., 2009, 15: 2335-2359.

[9] Tottrup C, Rasmussen M S. Mapping long-term changes in savannah crop productivity in Senegal through trend analysis of time series of remote sensing data., 2004, 103(3): 545-560.

[10] Hüttich C, Herold M, Schmullius C, Egorov V, Bartalev S A. Indicators of Northern Eurasia's land-cover change trends from SPOT-VEG ETATION time-series analysis 1998-2005., 2007, 28: 4199-4206.

[11] Symeonakis E, Drake N. Monitoring desertification and land degradation over sub-Saharan Africa., 2004, 25(3): 573-592.

[12] Sparks T H, Aasa A, Huber K, Wadsworth R. Changes and patterns in biologically relevant temperatures in Europe 1941-2000., 2009, 39: 191-207.

[13] Zhang K, Kimball J S, Mu Q, Jones L A, Goetz S J, Running S W. Satellite based analysis of northern ET trends and associated changes in the regional water balance from 1983 to 2005., 2009, 379(1): 92-110.

[14] White M A, Running S W, Thornton P E. The impact of growing-season length variability on carbon assimilation and evapotranspiration over 88 years in the eastern US deciduous forest., 1999, 42(3): 139-145.

[15] 謝寶妮, 秦占飛, 王洋, 常慶瑞. 基于遙感的黃土高原植被物候監測及其對氣候變化的響應. 農業工程學報, 2015, 31(15): 153-160.

Xie B N, Qin Z F, Wang Y, Chang Q R. Monitoring vegetation phenology and their response to climate change on Chinese Loess Plateau based on remote sensing., 2015, 31(15): 153-160. (in Chinese)

[16] 韋振鋒, 任志遠, 張翀. 1999-2010年陜甘寧黃土高原區氣候對植被物候的影響. 水土保持通報, 2014, 34(5): 232-236.

Wei Z F, Ren Z Y, Zhang C. Impact of climate on vegetation phenology on Loess Plateau in Shaanxi-Gansu-Ningxia Region during 1999-2010., 2014, 34(5): 232-236. (in Chinese)

[17] 張晗, 任志遠. 基于Whittaker濾波的陜西省植被物候特征. 中國沙漠, 2015, 35(4): 901-906.

Zhang H, Ren Z Y. Remote sensing analysis of vegetation phenology characteristics in Shaanxi province based on Whittaker Smooth method., 2015, 35(4): 901-906. (in Chinese)

[18] Yu H Y, Luedeling E, Xu J C. Winter and spring warming result in delayed spring phenology on the Tibetan Plateau., 2010, 107(51): 22151-22156.

[19] Zhang X Y, Friedl M A, Schaaf C B, Strahler A H, Hodges J C F, Gao F, Reed B C，Huete A. Monitoring vegetation phenology using MODIS., 2003, 84: 471-475.

[20] Kendall M G.. New York: Hafner Publishing Company, 1962.

[21] Julien Y, Sobrino J A. Global land surface phenology trends from GIMMS database., 2009, 30: 3495-3513.

[22] Zhou L M, Tucker C J, Kaufmann R K, MYNENI R B. Variation in northern vegetation activity inferred from satellite data of vegetation index during 1981 to 1999., 2001, 106: 20069-20083.

[23] Jeong S J, Ho C H, Gim H J, BROWN M E. Phenology shifts at start vs. end of growing season in temperate vegetation over the Northern Hemi-sphere for the period 1982-2008., 2001, 17: 2385-2399.

[24] Piao S L, Fang J Y, Zhou L M, Ciais P, Zhu B. Variations in satellite-derived phenology in China’s temperate vegetation., 2006, 12: 672-685.

[25] Zhang X Y, Friedl M A, Schaaf C B, Strahler A H, Hodges J C F, Gao F, Reed B C, Huete A. Monitoring vegetation phenology using MODIS., 2003, 84: 471-475.

[26] 俎佳星, 楊健. 東北地區植被物候時序變化. 生態學報, 2016, 36(7): 2015-2023.

ZU J X, YANG J. Temporal variation of vegetation phenology in northeastern China., 2016, 36(7): 2015-2023. (in Chinese)

[27] Wang C, Cao R Y, Chen J, Rao Y，Tang Y. Temperature sensitivity of spring vegetation phenology correlates to within-spring warming speed over the Northern Hemisphere., 2015, 50(3): 62-68.

（責任編輯 楊鑫浩）

Analysis of Temporal and Spatial Variation of Vegetation Phenology in the Loess Plateau

LI Qiang1,2, ZHANG Chong3, REN Zhi-yuan1

(1Tourism and Environment College of Shaanxi Normal University, Xi’an 710119;2Department of Environment and Resource Management, Shaanxi Xueqian Normal University, Xi’an 710100;3Shaanxi Key Laboratory of Disaster Monitoring and Mechanism Modeling, Baoji University of Arts and Sciences, Baoji 721013, Shaanxi)

【Objective】The Loess plateau is in the transitional region from wetness to dryness, from forest to grassland, from farming to animal husbandry, being the sensitive zone in climate change and agricultural development in China, the research on phenological feature of surface vegetation coverage in the region is of directive significance to agricultural production, environmental protection and ecological construction. Analysis of difference in phenological trends of vegetation in different time series and altitudes and hydrothermal conditions shall offer theoretical support and decision basis for current agricultural ecological environment improvement and sustainable development on loess plateau.【Method】Phenological feature values of vegetation on Loess Plateau every year were determined and phenological change trend was analyzed based on NDVI of ten-day values of SPOT VEGETATION from 1998 to 2012 and combined with harmonic analysis method and linear trend method.【Result】(1) From 1998 to 2012, the start of growing season advanced by 0.9 d on average every year and the end of growing season delayed by about 0.8 d on average every year, the length of growing season every year extended by 1.7 d on average under the joint action of advance at the start of growing season and delay at the end of growing season. (2) Hydrothermal condition on Loess Plateau has an immediate impact on phenological difference, the restrictive temperature for vegetation growth is 9℃，with the restrictive precipitation of 475 mm and 540 mm, respectively and restrictive altitude of 1 750 m. (3) The spatial partial correlation coefficients between the length trend of growing season of vegetation and altitude and air temperature are 0.0591and 0.0139 respectively, the spatial partial correlation coefficients between the length trend of growing season of vegetation and precipitation is -0.0174, therefore, the degree of correlation between three factors and the trend at the start of growing season is higher than that at the end of growing season.【Conclusion】The zones showing a significant and stable trend of phenological feature of vegetation on loess plateau are primarily distributed on plateau in the north of Shaanxi and Beishan in the middle of Shanxi. The phenological change in arid area and desert and grassland area in northwest is mainly subject to control by air temperature. The phenological change in semi-arid area and farming and grassland areas is mainly subject to control by precipitation. The phenological change in agricultural area in Fenwei Basin is subject to hydrothermal condition. Difference in hydrothermal condition has an insignificant impact on vegetation phenology in broad leaved forest zone. Altitude has an insignificant impact on change trend of vegetation phenology on loess plateau. The extension trend of growing season is on the increase with the increase of altitude and air temperature under the joint action of the start and end of growing reason, the shortening trend of growing season is on the increase with the increase of precipitation, the change characteristics of the same type of vegetations in terms of phenological trend based on the change of altitude, precipitation and air temperature are consistent, the change characteristics at the start of growing season have a greater influence on length change of growing season compared with that at the end of growing season.

vegetation coverage; growth season; harmonic analysis; Fourier interpolation; loess plateau

2015-11-20；接受日期：2016-09-30

國家自然基金青年科學基金（41301618）、陜西省2015年度科學技術研究發展計劃項目（陜西省青年科技新星項目2015KJXX-45）、2015年度陜西省博士后科研項目二等資助

李強，E-mail：liqiangis@163.com。通信作者任志遠，E-mail：renzhy@snnu.edu.cn