王志偉,王 茜,李世歌,王普昶,劉秀峰,謝彩云,史健宗,吳佳海,王小利,陸瑞霞,莫本田?剼?

(1.貴州省農(nóng)業(yè)科學(xué)院草業(yè)研究所，貴州 貴陽 550006；2.中國科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所 青藏高原冰凍圈觀測研究站，甘肅 蘭州 730020)

?

貴州喀斯特近30年植被生長特征分析

王志偉1,王 茜1,李世歌1,王普昶1,劉秀峰1,謝彩云1,史健宗2,吳佳海1,王小利1,陸瑞霞1,莫本田1?剼?

(1.貴州省農(nóng)業(yè)科學(xué)院草業(yè)研究所，貴州 貴陽 550006；2.中國科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所 青藏高原冰凍圈觀測研究站，甘肅 蘭州 730020)

貴州省地處南方喀斯特地區(qū)的中心，其生物學(xué)過程嚴(yán)重影響到環(huán)境變化。植物做為生物圈中的生產(chǎn)者，對陸地生態(tài)系統(tǒng)的貢獻(xiàn)同樣不容忽視。先進(jìn)型分辨率輻射儀的歸一化植被指數(shù)AVHRR GIMMS NDVI 3 g數(shù)據(jù)集作為一種有效的植被狀況監(jiān)測工具，可以提供全球范圍內(nèi)從20世紀(jì)80年代到21世紀(jì)初的資料。本研究使用該數(shù)據(jù)分析貴州省植被變化狀況，發(fā)現(xiàn)近30年來貴州省的植被整體上趨于生長增加的情形，其歸一化植被指數(shù)(NDVI)年增長率為0.002 1。不過，也存在兩個(gè)植被突變點(diǎn)(turing point，TP)，分別出現(xiàn)在1985年和1993年。在月NDVI數(shù)據(jù)中，各月整體變化率和TP與貴州省全省植被、栽培植物、草地和林地的NDVI變化趨勢則明顯不同。從空間分布角度分析，貴州省植被增長速度呈明顯的東南地區(qū)變化速率優(yōu)于西北地區(qū)，全球變暖是其中一方面的原因。同時(shí)，降水分布狀況也和植被變化狀況存在明顯的一致性。此外，人類活動對植被變化造成的影響也不容忽視。

植被；NDVI；喀斯特；石漠化；草地

植被作為陸地生態(tài)系統(tǒng)中最重要的組成部分，通過自身的光合作用和呼吸作用，以及對地表反照率和地表粗糙度的改變，在調(diào)節(jié)能量收支平衡和生物地球化學(xué)循環(huán)方面具有重要的意義[1]。其光合作用和呼吸作用的耦合，不僅可以影響地球氣候系統(tǒng)，同時(shí)還可以起到保持氣候穩(wěn)定的效果[1]。它不僅在全球物質(zhì)能量循環(huán)過程中起重要作用，而且對全球碳平衡的調(diào)節(jié)作用不可替代[2]。因此，分析植被特征在全球氣候變化狀況下的機(jī)理和規(guī)律，已經(jīng)成為眾多學(xué)者研究的熱點(diǎn)[3-7]。

有關(guān)植被狀況的監(jiān)測，可以借助衛(wèi)星遙感數(shù)據(jù)，該資料具有大范圍、重復(fù)觀測的特點(diǎn)，是長時(shí)間植被特征分析的有效工具[8-11]。最常用的方法是利用植被指數(shù)(vegetation indices，VI)數(shù)據(jù)，其中歸一化植被指數(shù)(normalized difference vegetation index，NDVI)的應(yīng)用最為廣泛[8,12]。NDVI是地物在電磁波譜范圍中紅波段與近紅外波段反射率的比值，與葉面積指數(shù)(leaf area index，LAI)、葉綠素濃度、光合有效輻射吸收比率(absorption of photosynthetically active radiation，fPAR)和總初級生產(chǎn)力(gross primary production quantity，GPP)的相關(guān)性很強(qiáng)[13]。目前，利用NDVI對植被機(jī)理的成熟研究很多[1,14-15]。在全球尺度上，利用1981-2006年的NDVI資料分析全球植被特征后表明：全球范圍內(nèi)遙感影像資料中，大部分的的像元值都在增長，特別是在北半球區(qū)域中植被的這種增長現(xiàn)象更加明顯[15]。在全國尺度上，對1982-1999年的NDVI資料計(jì)算后發(fā)現(xiàn)，18年來，中國植被四季平均NDVI均表現(xiàn)為上升趨勢，其中春季上升最為顯著，秋季上升趨勢最不顯著。

貴州省作為中國南方巖溶區(qū)的中心，是長江和珠江中上游最重要的生態(tài)屏障，具有獨(dú)特的喀斯特地貌環(huán)境單元[16]。其中，分布有大片連續(xù)的裸露碳酸鹽巖，這一地區(qū)生態(tài)環(huán)境脆弱，在資源開發(fā)和經(jīng)濟(jì)發(fā)展過程中存在嚴(yán)重的石漠化環(huán)境問題[17]。貴州巖溶區(qū)光照偏少、土層瘠薄、水土流失嚴(yán)重、雨水時(shí)空錯(cuò)位、旱澇災(zāi)害頻繁，影響籽實(shí)作物的品質(zhì)、產(chǎn)量及效益，但水熱條件好，尤其適合灌木、草本植物的生長，飼料資源豐富，適合草地畜牧業(yè)的發(fā)展[16]。發(fā)展草地畜牧業(yè)能保障巖溶山區(qū)農(nóng)民生計(jì)、促進(jìn)生態(tài)建設(shè)并維持其成果，調(diào)整農(nóng)村種植、產(chǎn)業(yè)及用地結(jié)構(gòu)，提高資源綜合利用效益，因此建立草地農(nóng)業(yè)生態(tài)系統(tǒng)是巖溶地區(qū)解決生態(tài)退化和石漠化的科學(xué)有效途徑[16-17]。

有關(guān)貴州省植被機(jī)理變化的研究主要在年尺度和季節(jié)尺度上展開[18-19]，時(shí)間范圍較短。而在30年時(shí)間尺度上，針對逐月數(shù)據(jù)分析巖溶中心貴州省NDVI變化狀況的研究和文獻(xiàn)，國內(nèi)外鮮見報(bào)道。為此擬采用1982-2011年的逐月NDVI數(shù)據(jù)，對貴州省的植被變化狀況進(jìn)行分析，研究貴州省的植被現(xiàn)狀，以期揭示30年來植被的動態(tài)變化規(guī)律，為喀斯特地區(qū)石漠化治理和生態(tài)環(huán)境保護(hù)提供數(shù)據(jù)支撐。

1 材料與方法

1.1 研究區(qū)概括

貴州喀斯特高原地處云貴高原東部[18]，位于我國西南部(圖1)，是青藏高原隆起時(shí)形成的一個(gè)海拔梯度大、地貌條件復(fù)雜、生態(tài)環(huán)境極為脆弱的獨(dú)特單元[16]，地勢起伏呈現(xiàn)明顯的東低西高趨勢。在碳酸鹽巖廣泛分布的地質(zhì)環(huán)境和溫暖濕潤季風(fēng)氣候的背景下，該區(qū)域海拔從90 m抬升至2 890 m，跨越24°-29° N，103°-109° E，總面積為1.76×105km2。其中，喀斯特地貌出露面積占73%，屬于我國乃至世界亞熱帶錐狀喀斯特分布區(qū)面積最大、發(fā)育最強(qiáng)烈的一個(gè)高原山區(qū)[19]。貴州省的植被類型圖(1∶1 000 000)來源于張新時(shí)[20]主編的《中華人民共和國植被圖》，電子地圖從中國科學(xué)院寒區(qū)旱區(qū)科學(xué)數(shù)據(jù)中心得到。通過該圖鑒計(jì)算得到貴州省植被類型主要包括栽培植被、灌叢、草叢、草甸、針葉林、闊葉林和針闊葉混交林，其占貴州省植被面積的比例分別為24.44%、35.24%、15.15%、2.70%、17.82%、4.64%和0.01%(圖1)。

圖1 研究區(qū)概況Fig.1 Geographical distribution of study area

1.2 數(shù)據(jù)源

本研究中使用的NDVI數(shù)據(jù)來源于先進(jìn)型分辨率輻射儀(advanced very high resolution radiometer，AVHRR)衛(wèi)星的全球庫存建模與制圖研究傳感器(global inventory modeling and mapping studies，GIMMS)。當(dāng)前，AVHRR NDVI產(chǎn)品已發(fā)展到第3代[21]，稱為AVHRR GIMMS NDVI 3 g產(chǎn)品，空間分辨率為(1/12)°。AVHRR GIMMS NDVI 3 g產(chǎn)品是唯一持續(xù)從1980s更新至2010s的NDVI數(shù)據(jù)集[14,22]。每半個(gè)月有一景NDVI產(chǎn)品，全年包括24景影像。去除1981年的半年數(shù)據(jù)，本研究選用的NDVI產(chǎn)品時(shí)間范圍為1982年1月到2011年12月，包含完整的30年數(shù)據(jù)，共720景影像。

為了去除云和雪，以及氣溶膠對NDVI數(shù)據(jù)的影響[21,23]，采用最大合成法，將每個(gè)月內(nèi)2景產(chǎn)品合成為1景月尺度數(shù)據(jù)，使得一年24景NDVI影像轉(zhuǎn)化成對應(yīng)于12個(gè)月的數(shù)據(jù)。同時(shí)，該產(chǎn)品已經(jīng)在部分地區(qū)經(jīng)過驗(yàn)證[5,12,24]，結(jié)果表明，GIMMS NDVI產(chǎn)品適宜于長時(shí)間序列的植被生長狀況分析[1]。

1.3 研究方法

將GIMMS NDVI 3 g原始數(shù)據(jù)旋轉(zhuǎn)180°，根據(jù)4個(gè)基準(zhǔn)點(diǎn)的地理定位坐標(biāo)分別為90.0-(1/24)°、-180.0+(1/24)°、-90.0+(1/24)°和180.0-(1/24)°，通過地理校準(zhǔn)方法(誤差控制在萬分之一以內(nèi))將NDVI數(shù)據(jù)轉(zhuǎn)化成投影類型為WGS-84的TIF文件，采用ArcGIS10.1處理。

NDVI數(shù)據(jù)中除NDVI值之外，還包括數(shù)據(jù)質(zhì)量控制f數(shù)值，計(jì)算方法如下：

f=NDVI-floor(NDVI/10)×10+1

(1)

NDVI=floor(NDVI/10)/1 000

(2)

式中：f表示質(zhì)量控制文件的數(shù)值，NDVI為數(shù)據(jù)值，floor()表示向下取整函數(shù)。f取值范圍從1到7，具體含義為，1-2、3-6和7分別代表該點(diǎn)象元值是優(yōu)良值、插值填補(bǔ)后數(shù)值(受云和傳感器影響)和缺失值。本研究中的NDVI數(shù)據(jù)全部選取f數(shù)值為1和2的象元值。

利用最小二乘線性回歸方法[1,5]分析貴州省1982年到2012年30年間的植被生長狀況，公式如下所示：

yt=a+bxt+ε

(3)

xt=1,2,3,…,30

式中：yt是NDVI序列，xt是第t年，a和b代表線性回歸方程的參數(shù)，a是截距，b是斜率，代表植被變化率，ε為方程誤差。

在長時(shí)間序列的NDVI數(shù)據(jù)分析中，通常會存在反映植被生長在時(shí)間尺度上產(chǎn)生突變的拐點(diǎn)(turing point，TP)，可以通過分段線性方程計(jì)算得到，方法如下：

(4)

式中：yt和xt的含義如式(3)中所示，j是GIMMS NDVI數(shù)據(jù)出現(xiàn)TP的年份，a0、b1和b2代表回歸系數(shù)，ε為誤差。在TP前的斜率是b1，在此之后斜率為b1+b2。在這里通過j∈ {4,6,…N-4} 來控制TP前后年份[1]。

2 結(jié)果與分析

根據(jù)貴州省已有的植被生長狀況研究結(jié)果可知[19,25]：無論在年尺度，還是在季節(jié)尺度，貴州省的植被生長狀況總體趨于更加旺盛的情形。不過，已有研究[19,25]的時(shí)間序列較短，對貴州省植被早期的分析[25]從1982年到1999年，近期的分析[19]從2000年到2006年。本研究將時(shí)間尺度擴(kuò)展為1982年到2011年，同時(shí)從月尺度上展開分析。

2.1 貴州省植被年變化的總趨勢

貴州省30年間的平均NDVI最大值為0.946 4，最小值為0.622 9，整體上而言，東部區(qū)域的植被生長狀況明顯好于西部區(qū)域(圖2)。其30年間的變化率顯示，大部分區(qū)域都呈現(xiàn)植被生長轉(zhuǎn)優(yōu)趨勢，僅畢節(jié)市和赫章縣及其周邊地區(qū)植被處于明顯轉(zhuǎn)劣的狀況。

貴州省的NDVI在30年間的變化中，存在兩個(gè)植被生長突變點(diǎn)(TP)，分別出現(xiàn)在1985年和1993年(圖3)，以上結(jié)果與在歐亞大陸[5]和亞洲大陸[1]得到的TP分別出現(xiàn)在1997年和2008年不同。不過，喜馬拉雅山脈克溪流域也發(fā)現(xiàn)兩個(gè)TP，分別出現(xiàn)在1994年和2000年[7]。

在1982年到2012年整個(gè)時(shí)間段內(nèi)，貴州省的植被NDVI數(shù)值，每年顯著增加0.002 1(P<0.001)。其中，1982年到1985年間，它每年的增長率為0.006 3(P<0.001)；1985年到1993年呈現(xiàn)相反的變化趨勢，每年存在0.004 9(P<0.001)的遞減率；1993年到2011年，再次出現(xiàn)每年0.003 9(P<0.001)的增長率。

不同植被類型具有不同的生物化學(xué)特征[26]，這些特性可能會導(dǎo)致不同的植被變化情況[2]。為了了解貴州省不同植被類型的變化規(guī)律和機(jī)理，根據(jù)我國1∶1 000 000植被類型圖[6,20]，貴州省的植被主要劃分為栽培植被、灌叢、草叢、草甸、針葉林、闊葉林和針闊葉混交林。其中，栽培植被在大類里被定義為栽培植被，其他6種定義為自然植被。

對栽培植被、草地和林地分析表明，其結(jié)果在變化趨勢上相似，在數(shù)值上存在差異(圖3)。在30年間整體上，栽培植被和林地的NDVI年增長速率一致，都為0.00 22(P<0.05)；草地略小于這兩類植被，NDVI增長率為0.002 1(P<0.05)。在1982年到1985年間，栽培植被和草地的NDVI增長率相似，分別為0.006 3(P<0.05)和0.006 6(P<0.05)；而林地的NDVI增長率較小，為0.005 4(P<0.05)。在1985年到1993年間，草地的NDVI減少率最大，為0.005 0(P<0.05)；之后依次為林地和栽培植被，分別為0.004 8(P<0.05)和0.004 6(P<0.05)。在1993年到2011年，栽培植被的NDVI增長率最高，為0.004 0(P<0.05)；草地和林地的NDVI增長率相同，為0.003 8(P<0.05)。

2.2 貴州省植被的月尺度變化規(guī)律

在1982年到2011年間，貴州省植被年內(nèi)NDVI月變化規(guī)律如圖4所示。1、2月份植被生長位于最低點(diǎn)，之后開始顯著增長，在7、8月份達(dá)到峰值，從9月份開始顯著下降，一直持續(xù)到12月份。同時(shí)，12月份的植被生長狀況要優(yōu)于1、2月份。

圖2 貴州省NDVI在1982-2011年間的分布狀況Fig.2 Geographical distribution of NDVI in Guizhou province from 1982 to 2011

圖3 NDVI在1982-2011年的年變化率曲線Fig.3 Interannual variations of annual mean NDVI from 1982 to 2011 for the entire Guizhou Province

表1展示了1982年到2011年貴州省年度內(nèi)各月的NDVI整體變化律，其各月30年間整體變化率與圖2中貴州省全省、栽培植物、草地和林地的NDVI變化狀況存在明顯差異。1月份每年有0.002 5的顯著降低趨勢(P<0.05)，4月和5月每年分別有0.005 0和0.004 4的顯著增長，同時(shí)11月和12月每年分別有0.002 4和0.001 5的增長(P>0.05)，其它月份在30年間的變化率并不明顯。

植被生長突變點(diǎn)(TP)在各月數(shù)據(jù)中出現(xiàn)的年份也明顯不同。存在兩個(gè)TP點(diǎn)數(shù)據(jù)的月份包括6個(gè)月，分別為1月(1994年和1999年)、2月(1987年和2000年)、4月(1990年和1994年)、6月(1989年和2008年)、9月(1986年和1993年)和10月(1989年和2003年)。存在一個(gè)TP的月份包括4個(gè)月份，分別為3月(1998年)、5月(2005年)、8月(1993年)和12月(1990年)。7月和11月不存在TP。

圖4 貴州省1982-2011年NDVI月變化規(guī)律Fig.4 NDVI distribution in different months of one year in Guizhou Province

注：誤差線大小為SE，樣本數(shù)為30，通過0.05顯著性檢驗(yàn)。

Note：Values are means±SE (n=30) and bars indicate SE. Different lower case letters indicate significant difference at 0.05 level. (Duncan’s test).

表1 1982-2011年貴州省年均各月NDVI整體變化率Table 1 The 30-year series of monthly mean NDVI date from 1982 to 2011 over Guizhou Province

注：以“*”表示在0.05水平上顯著相關(guān)；“**”表示在0.01水平上顯著相關(guān)。TP1為1985年的NDVI突變點(diǎn)，TP2為1993年的NDVI突變點(diǎn)。

Note：“*” and “**” indicate significant correlation at 0.05 and 0.01 level, respectively；TP1and TP2are the turning point of NDVI in 1985 and 1993, respectively.

3 討論

3.1 數(shù)據(jù)中的不穩(wěn)定因素

已有研究[27]表明，長時(shí)間序列的NDVI中誤差來源主要包括以下幾個(gè)方面：1)傳感器的不一致；2)天氣變化造成的干擾(如云、灰塵、氣溶膠和積雪覆蓋);3)地表變化引發(fā)的干擾(如火災(zāi)和植物病、蟲害)。

本研究中所使用的AVHRR GIMMS NDVI 3 g產(chǎn)品本身的時(shí)間覆蓋范圍十分廣泛，因此完全可以克服傳感器不一致的問題[1]。同時(shí)，很多學(xué)者[12,22,24,28-30]也指出NDVI最大合成方法能夠有效克服天氣原因造成的NDVI誤差。有關(guān)地表變化造成的NDVI變化，在本研究中不做具體分析。因?yàn)?，一方面相比無植被分布區(qū)域其變化并非十分劇烈[1]；另一方面GIMMS NDVI的像元分辨率較低，一個(gè)像元點(diǎn)覆蓋的面積約為64 km2，在其影像中會對一些局部區(qū)域做出“平滑”處理。

3.2 氣溫對植被變化的影響

近100年來，全球氣候逐漸變暖，特別是在20世紀(jì)末的20 多年中，增溫異常迅速，甚至在20世紀(jì)90年代時(shí)全球平均氣溫出現(xiàn)有觀測資料以來的最高值[31]。受全球變暖的影響，貴州省不同區(qū)域的植被也表現(xiàn)出不同的增長趨勢。同樣，已有研究[22]指出，在北半球緯度35°到75°的地域內(nèi)植被呈現(xiàn)顯著增長的狀況。P在歐亞大陸，1982年到2006年NDVI的年增長率為0.000 5[5]。同時(shí)，植被生長主要分為兩個(gè)階段。在1982年到1997年之間NDVI表現(xiàn)為增長，年增長率為0.001 8；在1997年到2006年之間NDVI表現(xiàn)為減少，年減少率為0.001 3。在全國范圍，指出在1982年到1999年這18年間，NDVI在春季上升最快、最顯著，平均每年以1.3%的速度增長[2]。青藏高原地區(qū)NDVI在1982-2003年的年增長率為0.000 5[32]。在更小的局部區(qū)域，如克溪地區(qū)的NDVI在1998-2011年的年增長率為0.000 8[7]。在1994年前NDVI的增長率為0.001 9，1994年到2000年NDVI減少了0.005 8，2000年之后NDVI的增長率為0.003 4。在甘南州夏河縣?？撇菰希瑢ι锪康难芯恐赋?，該區(qū)域的生物量干質(zhì)量在2000年到2013年的年增加趨勢為50.57 kg·hm-2(R2=0.590 7，P<0.001)[33]。同時(shí)，其它研究[34]也指出植物群落豐富度增加同氣候變暖的相關(guān)性很高。

3.3 其它因素對植被變化的影響

降雨與NDVI存在良好和顯著的線性關(guān)系，通常是影響植被生產(chǎn)力的最關(guān)鍵因素之一[35]。研究表明，貴州省的植被具有明顯的變好趨勢，同時(shí)降水也存在不同程度的升高，從區(qū)域上看，降水量呈現(xiàn)明顯的從東南部向西北部遞減的趨勢[36]，與本研究中NDVI變化率的分布結(jié)果一致。

地表植被變化是人類活動和氣候因素共同作用的結(jié)果[37]。貴州省許多地區(qū)存在人口密度高、人均受教育年限低、經(jīng)濟(jì)社會發(fā)展落后的特點(diǎn)，而且居民耕種方式落后、環(huán)保意識薄弱，進(jìn)而促使當(dāng)?shù)馗貑萎a(chǎn)低下[38]。傳統(tǒng)觀念中的小農(nóng)意識和“多子多?！庇^點(diǎn)直接導(dǎo)致當(dāng)?shù)厝丝谠鲩L速度明顯高于全國平均人口增速，當(dāng)現(xiàn)有人口增長超過土地資源本身的承載力時(shí)，在落后的生產(chǎn)方式條件下，會造成“越窮越生，越生越窮”的惡性循環(huán)。貧困，也不可避免的使得人們不斷的毀林開荒，來滿足對耕地和糧食的需求，進(jìn)而導(dǎo)致植被生態(tài)環(huán)境趨于惡化。

4 結(jié)論

基于AVHRR GIMMS NDVI 3 g產(chǎn)品，利用最小二乘線性回歸方法和分段線性方程可以獲取貴州省長時(shí)間序列的植被特征狀況。相比其它地區(qū)，貴州省在1982年到2011年間的植被生長具有更加明顯的增長趨勢，其NDVI年增長率為每年0.002 1。植被可以增加土壤入滲率，減小地表的水土流失[39]。因此，構(gòu)建良好的植被生長環(huán)境是治理土地石漠化的有效手段。特別是在喀斯特地區(qū)，地勢多變，山體眾多，更適合以小尺度種草養(yǎng)畜的方式來改善植被生長環(huán)境。同時(shí), 開發(fā)畜產(chǎn)品支撐我國畜產(chǎn)品有效供給，可以通過對草地畜牧業(yè)的整體規(guī)劃，來有效保障巖溶山區(qū)農(nóng)民就業(yè)、促進(jìn)生態(tài)環(huán)境建設(shè)、調(diào)整農(nóng)村種植產(chǎn)業(yè)及用地結(jié)構(gòu)，以此來更加合理、有效、科學(xué)的改良石漠化問題。

References:

[1] Chen B Z,Xu G,Coops N C,Ciais P,Innes J L，Wang G，Myneni R B，Wang T，Krzyzanowski J，Li Q，Cao L，Liu Y.Changes in vegetation photosynthetic activity trends across the Asia-Pacific region over the last three decades.Remote Sensing of Environment，2014，144:28-41.

[2] 樸世龍,方精云.1982-1999年我國陸地植被活動對氣候變化響應(yīng)的季節(jié)差異.地理學(xué)報(bào),2003,58(1):119-125.

Piao S L,Fang J Y.Seasonal changes in vegetation activity in response to climate changes in China between 1982 and 1999.Acta Geographica Sinica,2003,58(1):119-125.(in Chinese)

[3] Zhang G L,Dong J W,Xiao X M,Hu Z M,Sheldon S.Effectiveness of ecological restoration projects in Horqin Sandy Land,China based on SPOT-VGT NDVI data.Ecological Engineering,2012,38(1):20-29.

[4] Delbart N,Le Toan T,Kergoat L,Fedotova V.Remote sensing of spring phenology in boreal regions:A free of snow-effect method using NOAA-AVHRR and SPOT-VGT data (1982-2004).Remote Sensing of Environment,2006,101(1):52-62.

[5] Piao S L,Wang X H,Ciais P,Zhu B,Wang T,Liu J.Changes in satellite-derived vegetation growth trend in temperate and boreal Eurasia from 1982 to 2006.Global Change Biology,2011,17(10):3228-3239.

[6] Zhang G L,Zhang Y J,Dong J W,Xiao X.Green-up dates in the Tibetan Plateau have continuously advanced from 1982 to 2011.Proceedings of the National Academy of Sciences of the United States of America,2013,110(11):4309-4314.

[7] Zhang Y,Gao J,Liu L,Wang Z,Ding M,Yang X.NDVI-based vegetation changes and their responses to climate change from 1982 to 2011:A case study in the Koshi River Basin in the middle Himalayas.Global and Planetary Change,2013,108:139-148.

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

[9] 王勇,于信芳,莊大芳.近30年北京市高爾夫球場時(shí)空變化特征.草業(yè)學(xué)報(bào),2015,24(8):188-198.

Wang Y,Yu X F,Zhuang D F.Temporal and spatial characteristics of Beijing Golf Courses over the past 30 years.Acta Prataculturae Sinica,2015,24(8):188-198.(in Chinese)

[10] 趙鳳杰,王正浩,王慧萍,吳慧慧,劉航偉,王廣君,張澤華.不同密度短星翅蝗危害后羊草的高光譜變化及對產(chǎn)草量的影響.草業(yè)學(xué)報(bào),2015,24(3):195-203.

Zhao F J,Wang Z H,Wang H P,Wu H H,Liu H W,Wang G J,Zhang Z H.The effects of hyper spectal change on grassland biomass after damage byCalliptamusabbreviatespopulations of different densities.Acta Prataculturae Sinica,2015,24(3):195-203.(in Chinese)

[11] Che M,Chen B,Innes J L,Wang G,Dou X,Zhou T,Zhang H,Yan J,Xu G,Zhao H.Spatial and temporal variations in the end date of the vegetation growing season throughout the Qinghai-Tibetan Plateau from 1982 to 2011.Agricultural and Forest Meteorology,2014,189-190:81-90.

[12] Zeng F W,Collatz G,Pinzon J,Ivanoff A.Evaluating and quantifying the climate-driven interannual variability in global inventory modeling and mapping studies (GIMMS) normalized difference vegetation index (NDVI3g) at global scales.Remote Sensing,2013,5(8):3918-3950.

[13] Goetz S J,Prince S D.Modelling terrestrial carbon exchange and storage:Evidence and implications of functional convergence in light-use efficiency.Advances in Ecological Research,1999，28：57-92.

[14] Chen B,Zhang X,Tao J,Wu J,Wang J,Shi P,Zhang Y,Yu C.The impact of climate change and anthropogenic activities on alpine grassland over the Qinghai-Tibet Plateau.Agricultural and Forest Meteorology,2014,189-190:11-18.

[15] de Jong R,de Bruin S,de Wit A,Schaepman M E,Dent D L.Analysis of monotonic greening and browning trends from global NDVI time-series.Remote Sensing of Environment,2011,115(2):692-702.

[16] 高漸飛,蘇孝良,熊康寧,周瑋.貴州巖溶地區(qū)的草地生態(tài)環(huán)境與草地畜牧業(yè)發(fā)展.草業(yè)學(xué)報(bào),2011,20(4):279-286.

Gao J F,Su X L,Xiong K N,Zhou W.Grasslands ecoenvironment and stockbreeding development in the Karast areas of Guizhou Province.Acta Prataculturae Sinica,2011,20(4):279-286.(in Chinese)

[17] 熊康寧,李晉,李明忠.典型喀斯特石漠化治理區(qū)水土流失特征與關(guān)鍵問題.地理學(xué)報(bào),2012,67(7):878-888.

Xiong K N,Li J,Li M Z.Features of soil and water loss and key issues in demonstration areas for combating karst rocky desertification.Acta Geographica Sinica,2012,67(7):878-888.(in Chinese)

[18] 高江波,吳紹洪,蔡運(yùn)龍.區(qū)域植被覆蓋的多尺度空間變異性——以貴州喀斯特高原為例.地理研究,2013,32(12):2179-2188.

Gao J B,Wu S H,Cai Y L.Investigating the spatial heterogeneity of vegetation cover at multi scales:A case study in Karst Guizhou Plateau of China.Acta Geographica sinica,2013,32(12):2179-2188.(in Chinese)

[19] 王冰,楊勝天.基于NOAA/AVHRR的貴州喀斯特地區(qū)植被覆蓋變化研究.中國巖溶,2006,25(2):157-162.

Wang B,Yang S T.Change in vegetation cover in Guizhou karst regions based on NOAA/AVHRR.Carsologica sinica,2006,25(2):157-162.(in Chinese)

[20] 張新時(shí).中國植被類型圖.北京:地質(zhì)出版社,2007.

Zhang X S.Chinese Vegetation Type Map.Beijing:Geological Publishing House,2007.(in Chinese)

[21] Pinzon J,Tucker C.A Non-Stationary 1981-2012 AVHRR NDVI3g Time Series.Remote Sensing,2014,6(8):6929-6960.

[22] Beck H E,McVicar T R,van Dijk A I J M,Schellekens J,de Jeu R A M,Bruijnzeel L A.Global evaluation of four AVHRR NDVI data sets:Intercomparison and assessment against Landsat imagery.Remote Sensing of Environment,2011,115(10):2547-2563.

[23] Xu L,Myneni R B,Chapin Iii F S,Callaghan T V,Pinzon J E,Tucker C J,Zhu Z,Bi J,Ciais P,Tommervik H,Euskirchen E S,Forbes B C,Piao S L,Anderson B T,Ganguly S,Nemani R R,Goetz S J,Beck P S A,Bunn A G,Cao C,Stroeve J C.Temperature and vegetation seasonality diminishment over northern lands.Nature Climate Change,2013,3(6):581-586.

[24] Tucker C,Pinzon J,Brown M,Slayback D,Pak E,Mahoney R,Vermote E,El Saleous N.An extended AVHRR 8-km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data.International Journal of Remote Sensing,2005,26(20):4485-4498.

[25] 丹利,謝明.基于MODIS資料的貴州植被葉面積指數(shù)的時(shí)空變化及其對氣候的響應(yīng).氣候與環(huán)境研究,2009,14(5):455-464.

Dan L,Xie M.The spatiotemporal variation of leaf area index in Guizhou and its response to climate based on MODIS data.Climatic and Environmental Research,2009,14(5):455-464.(in Chinese)

[26] 梁天剛,馮琦勝,黃曉東,任繼周.草原綜合順序分類系統(tǒng)研究進(jìn)展.草業(yè)學(xué)報(bào),2011,20(5):252-258.

Lang T G,Feng Q S,Huang X D,Ren J Z.Review in the study of comprehensive sequential classification system of grassland.Acta Prataculturae Sinica,2011,20(5):252-258.(in Chinese)

[27] Forkel M,Carvalhais N,Verbesselt J,Mahecha M,Neigh C,Reichstein M.Trend change detection in NDVI time series:Effects of inter-annual variability and methodology.Remote Sensing,2013,5(5):2113-2144.

[28] Sobrino J A,Julien Y.Global trends in NDVI-derived parameters obtained from GIMMS data.International Journal of Remote Sensing,2011,32(15):4267-4279.

[29] Tian H,Cao C,Chen W,Bao S,Yang B,Myneni R B.Response of vegetation activity dynamic to climatic change and ecological restoration programs in Inner Mongolia from 2000 to 2012.Ecological Engineering,2015,82:276-289.

[30] Eckert S,Hüsler F,Liniger H,Hodel E.Trend analysis of MODIS NDVI time series for detecting land degradation and regeneration in Mongolia.Journal of Arid Environments,2015,113:16-28.

[31] 郭鈮,朱燕君,王介民,鄧朝平.近22年來西北不同類型植被NDVI變化與氣候因子的關(guān)系.植物生態(tài)學(xué)報(bào),2008,32(2):319-327.

Guo N,Zhu Y J,Wang J M,Deng C P.The relationship between NDVI and climate elements for 22 years in different vegetation areas of Northwest China.Journal of Plant Ecology,2008,32(2):319-327.(in Chinese)

[32] Peng J,Liu Z H,Liu Y H,Wu J,Han Y.Trend analysis of vegetation dynamics in Qinghai-Tibet Plateau using Hurst Exponent.Ecological Indicators,2012,14(1):28-39.

[33] 孟寶平,陳思宇,崔霞,馮琦勝,梁天剛.基于多源遙感數(shù)據(jù)的高寒草地生物量反演模型精度——以夏河縣桑科草原試驗(yàn)區(qū)為例.草業(yè)科學(xué),2015,32(11):1730-1739.

Meng B P,Chen S Y,Cui X,Feng Q S,Liang T G.The accuracy of grassland vegetation biomass estimated model based on multi-source remote sensing data——A case of experimental area in Sangke grassland in Xiahe County.Pratacultural Science,2015,32(11):1730-1739.(in Chinese)

[34] 任國華,鄧斌,后源.黃河源區(qū)沼澤濕地退化過程中植物群落特征的變化.草業(yè)科學(xué),2015,32(8):1222-1229.

Ren G H,Deng B,Hou Y.Changes of community characteristics in the degradation process of the alpine swamp wetland in the Yellow River Source area.Pratacultural Science,2015,32(8):1222-1229.(in Chinese)

[35] 陳操操,謝高地,甄霖,耿艷輝,冷允法.涇河流域植被覆蓋動態(tài)變化特征及其與降雨的關(guān)系.生態(tài)學(xué)報(bào),2008,28(3):925-938.

Chen C C,Xie G D,Zhen L,Geng Y H,Leng Y F.Analysis of Jinghe watershed vegetation dynamics evaluation of its relation to precipitation.Acta Ecologica Sinica,2008,28(3):925-938.(in Chinese)

[36] 鄭有飛,牛魯燕,吳榮軍,吳戰(zhàn)平,駱興江,蔡子穎.1982―2003年貴州省植被覆蓋變化及其對氣候變化的響應(yīng).生態(tài)學(xué)雜志,2009,28(9):1773-1778.

Zheng Y F,Niu L Y,Wu R J,Wu Z P,Luo X J,Cai Z Y.Vegetation cover change in Guizhou of Southwest China in 1982―2003 in response to climate change.Chinese of Ecology,2009,28(9):1773-1778.(in Chinese)

[37] 周洪建,王靜愛,岳耀杰,李睿.人類活動對植被退化/恢復(fù)影響的空間格局——以陜西省為例.生態(tài)學(xué)報(bào),2009,29(9):4847-4856.

Zhou H J,Wang J A,Yue Y J,Li R.Research on spatial pattern of human induced vegetation degradation and restoration:A case study of Shaanxi Province.Acta Ecologica Sinica,2009,29(9):4847-4856.(in Chinese)

[38] 龔勤林,陳說.西南喀斯特地區(qū)資源環(huán)境保護(hù)與生態(tài)文明建設(shè)研究.西部發(fā)展評論,2014(7):44-53.

Gong Q L,Chen S.The study of resources and environment protection and ecological civilization construction in the southwest karst area.The Western Region Development Reviews,2014(7):44-53.(in Chinese)

[39] 高洋,王根緒,高永恒.長江源區(qū)高寒草地土壤有機(jī)質(zhì)和氮磷含量的分布特征.草業(yè)科學(xué),2015,32(10):1548-1554.

Gao Y,Wang G X,Gao Y H.Distribution characteristics of soil organic matter and nitrogen,phosphor content in alpine grassland ecosystem in upper Yangtze River.Pratacultural Science,2015,32(10):1548-1554.(in Chinese)

(責(zé)任編輯 王芳)

The growth characteristics of vegetation in the Karst regions of Guizhou from 1981 to 2012

Wang Zhi-wei1, Wang Qian1, Li Shi-ge1, Wang Pu-chang1, Liu Xiu-feng1, Xie Cai-yun1,Shi Jian-zong2, Wu Jia-hai1, Wang Xiao-li1, Lu Rui-xia1, Mo Ben-tian1

(1.Guizhou Institute of Prataculture, Guizhou Academy of Agricultural Sciences, Guiyang 550006, China;2.Cryosphere Research Station on the Qinghai-Tibet Plateau, Cold and Arid Regions Environmental and Engineer Research Institute, Chinese Academy of Sciences, Lanzhou 730020, China)

Guizhou province is located in the center area of Karst region, its ecological environment has been changed by biological process significantly. It is noted that the plants play an important role in the terrestrial ecosystem as the producer of the biosphere. Therefore, the updated Global Inventory Modeling and Mapping Studies (GIMMS) third generation Normalized Difference Vegetation Index (NDVI) dataset by using global satellite Advanced Very High Resolution Radiometer (AVHRR) was used to investigate the changes of vegetation in the Guizhou Province from 1982 to 2011. The results demonstrated that the overall NDVI increased from 1982 to 2011 with an average rate of 0.002 1 per year (P<0.001). However, there still were two turning points (TP) in 1985 and 1993, respectively. There are large variations between the NDVI trend and TPs in the whole area and in different vegetation types. From the spatial distribution of vegetation growth in Guizhou, the rate of vegetation increased in southeastern was superior to northwestern obviously. One reason was global warming, and meanwhile the vegetation change consisitent with the spatial distribution of precipitation. In addition, the impact of human activities on the changes of vegetation should be paid more attentions.

vegetation; NDVI; Karst region; rocky desertification; grassland

Mo Ben-tian E-mail:gzcymbt163@sohu.com

10.11829/j.issn.1001-0629.2015-0662

2015-11-30 接受日期：2016-04-05

貴州省農(nóng)業(yè)科學(xué)院科研項(xiàng)目([2016]貴州省近三十年來植被生長特征分析)；貴州省農(nóng)業(yè)科學(xué)院專項(xiàng)資金([2016]032號)；貴州省科技創(chuàng)新人才團(tuán)隊(duì)項(xiàng)目(黔科合平臺人才[2016]5617號)

王志偉(1983-)，男，陜西府谷人，助理研究員，博士，研究方向?yàn)椴莸剡b感。E-mail:wzw1206@163.com共同

第一作者：王茜(1984-)，女，甘肅蘭州人，助理研究員，博士，研究方向?yàn)椴輼I(yè)科學(xué)。E-mail：snoopy0729@163.com

莫本田(1967-)，男(布依族)，貴州獨(dú)山人，研究員，學(xué)士，研究方向?yàn)椴輼I(yè)科學(xué)。E-mail：gzcymbt163@sohu.com

S812;Q945.3

A

1001-0629(2016)11-2188-09*

王志偉,王茜,李世歌,王普昶,劉秀峰,謝彩云,史健宗,吳佳海,王小利,陸瑞霞,莫本田.貴州喀斯特近30年植被生長特征分析.草業(yè)科學(xué),2016,33(11)：2180-2188.

Wang Z W,Wang Q,Li S G,Wang P C,Liu X F,Xie C Y,Shi J Z,Wu J H,Wang X L,Lu R X,Mo B T.The growth characteristics of vegetation in the Karst regions of Guizhou from 1981 to 2012.Pratacultural Science，2016,33(11)：2180-2188.