劉志軍，曾業(yè)隆，嚴(yán) 蘭，譚 偉

(貴州大學(xué) 林學(xué)院，貴州 貴陽(yáng) 550025)

植被是覆蓋地表的植物群落的總稱，是全球生態(tài)系統(tǒng)的重要組成部分，其覆蓋狀態(tài)是生態(tài)系統(tǒng)良性循環(huán)的關(guān)鍵診斷參數(shù)，對(duì)全球物質(zhì)循環(huán)、能量平衡以及生態(tài)系統(tǒng)的水循環(huán)、碳循環(huán)、氣候變化、石漠化和水土流失方面起到了極其重要的作用[1-4]。植被指數(shù)是衡量地表植被狀況、指示生態(tài)環(huán)境變化的一個(gè)重要指標(biāo)[5]，在眾多植被指數(shù)中，歸一化植被指數(shù)(Normalized Difference Vegetation Index，簡(jiǎn)稱NDVI)是應(yīng)用最廣泛的植被指數(shù)[6]，利用NDVI不僅能較為真實(shí)地反映地表植被覆蓋狀況，也能大范圍長(zhǎng)時(shí)間序列地監(jiān)測(cè)地表植被覆蓋的變化[7-9]。因此，利用NDVI的時(shí)間序列數(shù)據(jù)研究植被的時(shí)空的變化對(duì)生態(tài)系統(tǒng)及其生態(tài)效應(yīng)具有重要意義。國(guó)內(nèi)外許多學(xué)者對(duì)不同區(qū)域的NDVI時(shí)空變化特征進(jìn)行了相關(guān)的研究。Park H S等[10]基于NDVI數(shù)據(jù)研究了東亞植被從1982年到2006年的時(shí)間變化特征；Piao S等[11]結(jié)合NDVI數(shù)據(jù)分析了1982～2006年歐亞大陸植被生長(zhǎng)變化的時(shí)空格局；劉靜等[12]利用NDVI數(shù)據(jù)采用轉(zhuǎn)移矩陣的方法，研究了毛烏素沙地1990～2007年間的植被覆蓋變化趨勢(shì)；戴聲佩等[13]利用GMMS/NDVI數(shù)據(jù)分析了我國(guó)西北地區(qū)1982～2006年植被的時(shí)空變化特征，研究結(jié)果認(rèn)為我國(guó)西北部地區(qū)近25年來(lái)NDVI呈增長(zhǎng)趨勢(shì)，且增速為0.5%/10a；進(jìn)一步地，謝曉華等[14]應(yīng)用像元二分原理，揭示了閩江流域植被覆蓋度的時(shí)空變化規(guī)律。在植被動(dòng)態(tài)變化趨勢(shì)方法上，劉亞龍等[15]利用非參數(shù)檢驗(yàn)方法，分析了膠東半島最近10年NDVI的變化趨勢(shì)，研究結(jié)果表明膠東半島的NDVI變化趨勢(shì)以衰減區(qū)域占主導(dǎo)地位；此外，王新明等[16]將R/S方法用于NDVI的時(shí)間序列分析，這種方法可以定量描述時(shí)間序列的持久性或反持久性，進(jìn)一步說(shuō)明了中國(guó)西部NDVI時(shí)間序列表現(xiàn)出更強(qiáng)的的持久性。王文輝等[17]采用逐像元趨勢(shì)分析，探究了福建長(zhǎng)汀植被覆蓋度的變化；袁麗華等[18]采用Theil-Sen Median趨勢(shì)分析和Mann-Kendall檢驗(yàn)，研究了黃河流域NDVI的變化趨勢(shì)特征。

云貴高原是中國(guó)的第四大高原，中部是我國(guó)石漠化分布相對(duì)比較集中的地區(qū)，石漠化面積占石漠化總面積53.4%，且大部分區(qū)域?qū)儆诳λ固氐孛?，受特殊地質(zhì)地貌條件的影響，該區(qū)域生態(tài)環(huán)境十分脆弱，植被變化的穩(wěn)定性尚未明確。然而目前的研究成果主要集中在空間范圍較大的西南地區(qū)或縣級(jí)區(qū)域尺度，針對(duì)云貴高原整體區(qū)域NDVI時(shí)空變化的研究較少[19]；另外，很少有研究者對(duì)云貴高原NDVI的穩(wěn)定性進(jìn)行評(píng)價(jià)。鑒于上述情況此，本文利用2001～2014年的MODIS NDVI數(shù)據(jù)，結(jié)合趨勢(shì)分析和穩(wěn)定性分析方法，探討云貴高原近14年來(lái)植被NDVI時(shí)空變化特征及穩(wěn)定性，以期為該區(qū)域的生態(tài)環(huán)境保護(hù)和石漠化治理提供科學(xué)依據(jù)。

1 材料和方法

1.1 研究區(qū)概況

云貴高原是世界上喀斯特地貌發(fā)育最典型的地區(qū)之一，該區(qū)域水土流失強(qiáng)度大、石漠化嚴(yán)重、生態(tài)環(huán)境脆弱。云貴高原地區(qū)通常指云南省和貴州省，地處位于東經(jīng)100°～111°，北緯22°～30°之間，土地面積55.92萬(wàn)km2，占全國(guó)國(guó)土面積的5.8%(如圖1)。該區(qū)域地勢(shì)由北向南呈階梯式下降，西北部和東南部分別為最高和最低地帶，海拔介于400～3500 m，屬亞熱帶濕潤(rùn)區(qū)。此外，云貴高原氣候?yàn)閬啛釒Ъ撅L(fēng)氣候，氣候差別顯著，是中國(guó)森林植被類型最為豐富的區(qū)域，動(dòng)植物資源豐富。

1.2 數(shù)據(jù)來(lái)源與處理

NDVI數(shù)據(jù)為地理空間數(shù)據(jù)云得到的中國(guó)MODND1M月合成產(chǎn)品，來(lái)源于中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http：//www.gscloud.cn)；行政邊界數(shù)據(jù)為全國(guó)1∶400萬(wàn)基礎(chǔ)地理信息數(shù)據(jù)，來(lái)源于國(guó)家基礎(chǔ)地理信息中心。MODIS NDVI數(shù)據(jù)起止時(shí)間為2001年1月到2014年12月，空間分辨率為500 m，具有較高精度，時(shí)間分辨率為逐月，并用輻射定標(biāo)、大氣校正和去云等手段對(duì)原始數(shù)據(jù)進(jìn)行了預(yù)處理，數(shù)據(jù)使用前先對(duì)其進(jìn)行逐圖檢查，對(duì)異常值數(shù)據(jù)進(jìn)行替換處理。

1.3 研究方法

1.3.1均值法 為能有效反映出研究區(qū)年內(nèi)植被NDVI的總體情況，本研究通過(guò)計(jì)算年內(nèi)每月的NDVI平均值，代表當(dāng)年的NDVI值，公式如下：

圖1 云貴高原示意圖Fig.1 Map of the Yunnan-Guizhou plateau

(1)

此外，為能有效反映研究區(qū)14年植被NDVI的總體特征，本研究對(duì)2001-2014年植被的年累計(jì)平均值進(jìn)行計(jì)算。公式如下：

(2)

1.3.2一元線性回歸趨勢(shì)分析 在分析植被NDVI趨勢(shì)變化的方法中，一元線性回歸分析方法使用較為實(shí)用，該方法可以模擬每個(gè)柵格的變化趨勢(shì)，可以準(zhǔn)確真實(shí)地反映長(zhǎng)時(shí)間序列的植被NDVI變化特征[20]。張勃等[ 21]和Stow D等[22]利用此方法研究了植被覆蓋度的變化特征。本文利用一元線性回歸分析的方法，對(duì)研究區(qū)NDVI進(jìn)行逐像元線性回歸分析，得到逐個(gè)像元在時(shí)間內(nèi)的變化特征和趨勢(shì)。公式如下：

(3)

式中：θslope表示變化趨勢(shì)的斜率；n表示總的年數(shù)；NDVIi表示第i年的年均NDVI值；θslope>0表示植被覆蓋呈增加趨勢(shì)，θslope<0表示植被覆蓋呈減少趨勢(shì)，θslope=0表示植被覆蓋無(wú)變化趨勢(shì)。本文采用F檢驗(yàn)對(duì)趨勢(shì)進(jìn)行顯著性檢驗(yàn)[23]，根據(jù)檢驗(yàn)結(jié)果將變化趨勢(shì)結(jié)果劃分為五個(gè)等級(jí)：極顯著減少(θslope<0，p<0.01)；顯著減少(θslope<0，0.010.05)；顯著增加(θslope>0，0.010，p<0.01)。

1.3.3變異系數(shù)法 變異系數(shù)法可以分析植被NDVI的波動(dòng)情況，劉洋等[24]和Milich L等[25]利用此方法分析了植被覆蓋度的波動(dòng)規(guī)律，李卓等[26]利用此方法分析了京津冀地區(qū)植被覆蓋度的穩(wěn)定性。本文通過(guò)逐像元計(jì)算云貴高原2001～2014年NDVI年均值的變異系數(shù)，可以得到NDVI在時(shí)間序列上的穩(wěn)定性。公式如下：

(4)

2 結(jié)果與分析

2.1 云貴高原NDVI時(shí)空分布格局

為了研究云貴高原植被NDVI隨時(shí)間變化的特點(diǎn)，采用均值法對(duì)2001～2014年14a年間各年、各月求平均值，能有效反映出研究區(qū)年和月的變化。圖2是用均值法求出2001～2014年每年平均NDVI值的年際變化情況。圖2表明，在2001～2014年間，云貴高原植被年平均NDVI值整體呈現(xiàn)上升趨勢(shì)，年均值在0.63～0.73之間波動(dòng)，數(shù)值變化范圍較小，增速為0.00395/a。2002～2005年NDVI值整體呈現(xiàn)下降的趨勢(shì)，在0.672～0.631之間波動(dòng)，2005年達(dá)到14a間的最低值；2005年之后，除2010年、2012、2014年有較大的下降外，NDVI值整體呈現(xiàn)上升趨勢(shì)，2013年達(dá)到最大值。查閱相關(guān)文獻(xiàn)[27]發(fā)現(xiàn)，云貴高原在2005、2010、2012、2014年發(fā)生了不同程度的干旱，由此推測(cè)干旱災(zāi)害對(duì)植被NDVI變化造成影響。因此，從年際變化來(lái)看，14a間云貴高原的植被覆蓋總體向好的方向發(fā)展，植被覆蓋得到很大的提高，說(shuō)明云貴高原的植被逐漸變好。

圖3是用均值法求出2001～2014年每個(gè)月平均NDVI值的變化情況。圖3表明，云貴高原植被NDVI在年內(nèi)分布整體呈單峰曲線，分布不均勻，月均值在0.599～0.777之間波動(dòng)。1～3月曲線呈緩慢下降的趨勢(shì)，NDVI值較均勻，出現(xiàn)一年中的最低值0.599；從3月份開(kāi)始，隨溫度升高，降雨量增加，植被進(jìn)入發(fā)芽期，NDVI值逐漸增大，9月達(dá)到一年中的最大值，接近0.777；9月之后，氣溫開(kāi)始下降，降雨量降少，植被進(jìn)入休眠期，NDVI值開(kāi)始逐漸下降，減小的幅度較大。因此從14年總體來(lái)看，3～9月植被覆蓋度逐漸增加，9月份達(dá)到最大；8～11月植被覆蓋度沒(méi)有較大的變化，維持在較高值；11月到第二年的3月植被覆蓋度下降幅度最大，降到一年中最低值。

圖2 2001～2014年云貴高原NDVI年際變化Fig.2 Inter-annual variation of the NDVI in Yunnan-Guizhou plateau from 2001 to 2014

圖3 2001～2014年云貴高原NDVI月變化Fig.3 Monthly variation of the NDVI in Yunnan-Guizhou plateau from 2001 to 2014

云貴高原NDVI季節(jié)空間變化具有明顯的變化規(guī)律。統(tǒng)計(jì)得出云貴高原NDVI季節(jié)空間變化具有明顯的變化特點(diǎn)，從圖4可以看出，變化主要集中在中部和東部地區(qū)，植被NDVI隨季節(jié)的變化呈現(xiàn)由東向西先增加后減少的變化規(guī)律，西南地區(qū)植被NDVI隨季節(jié)變化不明顯。夏季植被NDVI值最高，其次是秋季，冬季的植被NDVI值最低。西南地區(qū)主要以常綠植被為主，對(duì)季節(jié)的響應(yīng)不明顯，東南部植被覆蓋較高，有部分常綠植被，受季節(jié)影響較小，中部喀斯特地區(qū)，植被稀疏，受季節(jié)變化影響最大。

圖4 2001～2014年云貴高原NDVI季節(jié)變化Fig.4 Seasonal variation of the NDVI in Yunnan-Guizhou plateau from 2001 to 2014

根據(jù)2001～2014年的年均NDVI數(shù)據(jù)，運(yùn)用ArcGIS軟件逐像元計(jì)算14a平均值得到云貴高原平均NDVI空間分布圖(圖5)。由于受海拔和特殊地貌的影響，云貴高原的NDVI空間分布具有顯著性的差異性，如圖5云貴高原的NDVI在空間整體分布呈西南、東南高，西北、中部低的空間格局，均值為0.68。高值區(qū)主要分布在西南、東南地區(qū)，其最高值在南部的西雙版納地區(qū)，均值為0.80，主要原因是該地區(qū)氣候?yàn)闊釒Ъ撅L(fēng)氣候，氣候溫暖濕潤(rùn)，有大量的熱帶雨林，其次是東南部的黔東南地區(qū)和東北部的赤水地區(qū)，主要原因是黔東南地區(qū)原始生態(tài)保存完好，境內(nèi)有雷公山、云臺(tái)山、佛頂山等原始森林，赤水境內(nèi)有亞熱帶常綠闊葉林原生植被帶。低值區(qū)主要分布在中部和西北部，中部的畢節(jié)地區(qū)、曲靖市、六盤水市是云貴高原的中心，是中國(guó)石漠化分布相對(duì)較集中的地區(qū)，嚴(yán)重的石漠化導(dǎo)致植被稀疏，西北部的迪慶藏族自治州是云貴高原向青藏高原的過(guò)渡帶，這里海拔較高是北半球緯度最低的雪山群，高寒植被區(qū)。圖中紅色部分是NDVI值小于0.1的無(wú)植被區(qū)，主要是城市區(qū)、滇池、洱海、撫仙湖、程海、梅里雪山、玉龍雪山等。

圖5 2001～2014年云貴高原年均NDVI空間分布Fig.5 Spatial distribution of average annual DNVI in Yunnan-Guizhou plateau from 2001 to 2014

2.2 云貴高原NDVI趨勢(shì)分析

本文將一元線性回歸趨勢(shì)分析與F檢驗(yàn)結(jié)合起來(lái)，用來(lái)反映云貴高原2001～2014年每個(gè)柵格的植被NDVI變化趨勢(shì)。圖6是把一元線性回歸趨勢(shì)分析的結(jié)果和F檢驗(yàn)的結(jié)果相結(jié)合，得到逐像元的2001～2014年云貴高原植被NDVI的變化趨勢(shì)的空間分布特征。從植被NDVI變化趨勢(shì)的空間分布可以看出：2001～2014年云貴高原植被NDVI值主要以增加為主，且增加的區(qū)域遠(yuǎn)遠(yuǎn)大于減少的區(qū)域。植被NDVI極顯著增加的區(qū)域主要分布在云貴高原中部的畢節(jié)地區(qū)、昭通地區(qū)、曲靖市、文山地區(qū)，以及云貴高原東北部的遵義地區(qū)的東北部和西北部的臨滄市，主要原因是這些區(qū)域植被面積增加并大部分得到改善，以喀斯特地區(qū)石漠化分布集中的畢節(jié)、六盤水、黔西南為例，通過(guò)退耕還林、封山育林以及水土保持工程對(duì)石漠化區(qū)域進(jìn)行治理，使得植被面積明顯增加。植被NDVI顯著性增加的區(qū)域主要分布在云貴高原的東部和北部的大部分地區(qū)，中部主要呈現(xiàn)零散分布，主要原因是，人類活動(dòng)近14a來(lái)未對(duì)植被進(jìn)行破壞，保護(hù)環(huán)境意識(shí)加強(qiáng)。變化不明顯的區(qū)域主要分布在云貴高原的南部西雙版地區(qū)、東部的黔東南地區(qū)和遵義地區(qū)的中部，以及西部的大理市、怒江、迪慶、德宏地區(qū)，以上這些區(qū)域，植被覆蓋度較好，并沒(méi)有遭到嚴(yán)重破壞，如西雙版納地區(qū)的熱帶雨林，黔東南的原始森林，赤水地區(qū)的原生植被帶。植被NDVI顯著減少的區(qū)域主要分布在昆明市、貴陽(yáng)市、曲靖市、大理市、保山市、銅仁市的城市邊緣，紅河地區(qū)的東南部也有減少區(qū)域，同時(shí)在其它地州市也有零星的減少區(qū)域。因?yàn)檫@些區(qū)域城市化擴(kuò)展導(dǎo)致植被面積減少，特別地，在喀斯特地區(qū)，生態(tài)脆弱，受人為干擾強(qiáng)度大，導(dǎo)致植被覆蓋度減少程度明顯。

圖6 2001～2014年云貴高原年均NDVI變化趨勢(shì) Fig.6 Spatial distribution of the coefficient of variation of the inter-annual NDVI and precitation in Yunnan-Guizhou plateau from 2001 to 2014

表1是云貴高原NDVI的變化趨勢(shì)統(tǒng)計(jì)，可以看出14a間全區(qū)域的NDVI值變化趨勢(shì)以極顯著增加為主，極顯著和顯著減少的比例較少，各類變化趨勢(shì)的占總面積的比例：極顯著增加>變化不顯著>顯著增加>極顯著減少>顯著減少。統(tǒng)計(jì)表明，有占總面積59.9%的區(qū)域植被覆蓋呈現(xiàn)顯著性增長(zhǎng)的趨勢(shì)，有占總面積0.7%的區(qū)域植被覆蓋度呈現(xiàn)顯著性減少的趨勢(shì)，其余占總面積39.4%的區(qū)域植被變化不顯著。

表1 NDVI變化趨勢(shì)顯著統(tǒng)計(jì)Tab.1 Significance of the trend of NDVI

以上結(jié)果說(shuō)明，云貴高原植被NDVI整體為增加的趨勢(shì)，增加的區(qū)域主要是石漠化較集中的畢節(jié)地區(qū)、六盤水、文山地區(qū)，顯著減少的趨勢(shì)較少，主要分布在人口相對(duì)頻繁的城市周邊和喀斯特生態(tài)脆弱的地區(qū)。

表2 NDVI變異系數(shù)統(tǒng)計(jì)Tab.2 Variation coefficient of NDVI

圖7 2001～2014年云貴高原年均NDVI變異程度 Fig.7 Spatial distribution of coefficient of variation of the inter-annual NDVI in Yunnan-Guizhou plateau from 2001 to 2014

2.3 云貴高原植被NDVI空間穩(wěn)定性評(píng)價(jià)

圖7是由ArcGIS軟件求出2001-2014年14a間云貴高原逐個(gè)像元的c.vNDVI值，根據(jù)以往研究結(jié)果，本文參照以往劉洋等[24]和李卓等[26]的研究結(jié)果將云貴高原c.vNDVI變異值分為三個(gè)等級(jí)：非常穩(wěn)定(c.vNDVI≤0.05)、穩(wěn)定(0.050.1)。從圖7和表2可以看出，2001～2014年云貴高原植被NDVI值總體處于中等波動(dòng)變化，高波動(dòng)變化所占的比例較少。如圖低波動(dòng)區(qū)域占總面積的37.3%，主要分布在北部地區(qū)，以西雙版納、大理、德宏分布較集中，其次是東南部的黔東南、銅仁中部地區(qū)，以上這些地區(qū)水熱條件和原生植被較好，故植被波動(dòng)變化較小；中等波動(dòng)變化區(qū)域占總面積的59.0%，主要分布在畢節(jié)中部以東地區(qū)、六盤水、曲靖北部地區(qū)，由于該地區(qū)植被主要以灌木為主，植被稀疏。高波動(dòng)區(qū)域占總面積的3.7%，主要分布石漠化嚴(yán)重的畢節(jié)西部地區(qū)、昭通西部和高海拔地區(qū)的迪慶和怒江，該地區(qū)植被受海拔和降水的影響較大，故波動(dòng)性較高。

3 結(jié)論與討論

本研究分析了云貴高原近14年來(lái)植被NDVI時(shí)空分布格局、變化趨勢(shì)和穩(wěn)定性等特征，對(duì)NDVI空間分布研究，得到空間分布上具有明顯的差異性，東南部與西南部明顯高于中部，西北部地區(qū)較低；對(duì)NDVI趨勢(shì)分析，得到云貴高原植被整體逐漸變好，特別是石漠化嚴(yán)重地區(qū)的植被呈增加趨勢(shì)，說(shuō)明喀斯特地區(qū)石漠化治理工程和生態(tài)恢復(fù)工程對(duì)植被的快速恢復(fù)起到重要作用；對(duì)近14年NDVI年均值分析，得到年均值在0.63～0.73之間波動(dòng)，整體呈現(xiàn)上升的趨勢(shì)，但在2005年達(dá)到14年間的最低值，分析發(fā)現(xiàn)2005年主要是云南遭受近50年來(lái)最大干旱和貴州發(fā)生嚴(yán)重夏旱[28]，說(shuō)明極端干旱氣候事件對(duì)植被NDVI的變化會(huì)產(chǎn)生一定影響。此外對(duì)空間穩(wěn)定性分析，得到處于低波動(dòng)變化占比為37.3%，主要分布在自然保護(hù)區(qū)和原始生態(tài)保持較好地區(qū)；中等波動(dòng)變化占比為59.0%，主要分布在以灌木為主的植被稀疏地區(qū)；高波動(dòng)變化占比為3.7%，主要分布在石漠化嚴(yán)重地區(qū)和高海拔地區(qū)，說(shuō)明NDVI穩(wěn)定性受植被類型和人類活動(dòng)及氣候影響較大。

本研究的結(jié)果與華維等[29]和張勃等[20]研究的西南地區(qū)植被覆蓋變化特征結(jié)果相似，同時(shí)與劉世梁等[30]研究云南省植被NDVI時(shí)空變化特征和李仕容等[31]研究貴州省植被覆蓋度時(shí)空變化特征與趨勢(shì)的結(jié)果相似，說(shuō)明本研究結(jié)果具有較大利用價(jià)值。此外我們分析發(fā)現(xiàn)不同數(shù)據(jù)源監(jiān)測(cè)得到的結(jié)果會(huì)存在差異，本文研究使用的MODIS NDVI數(shù)據(jù)與前人研究使用的SPOT VGT數(shù)據(jù)不同，并且在空間分辨率和時(shí)間上也不同，因此不同數(shù)據(jù)源對(duì)NDVI變化趨勢(shì)和穩(wěn)定性的影響有待于進(jìn)一步探討。同時(shí)地形地貌對(duì)植被NDVI時(shí)空分布具有一定影響，下一步研究可以考慮喀斯特和石漠化對(duì)NDVI的變化和穩(wěn)定性有何影響。此外，人類活動(dòng)和極端氣候?qū)χ脖籒DVI的變化起到重要作用，如何科學(xué)評(píng)估生態(tài)工程和干旱氣候的對(duì)植被NDVI的影響，是今后云貴高原生態(tài)恢復(fù)和生態(tài)保護(hù)的重要參考，也是未來(lái)研究的重點(diǎn)。

參考文獻(xiàn)：

[1] 張 萍，曾信波.植被蓄水保土功能研究[J]. 山地農(nóng)業(yè)生物學(xué)報(bào)，1999(5)：300-304.

[2] 李玉山.黃土高原森林植被對(duì)陸地水循環(huán)影響的研究[J]. 自然資源報(bào)，2001(05)：427-432.

[3] 王德?tīng)t，朱守謙，黃寶龍.石漠化過(guò)程中土壤理化性質(zhì)變化的初步研究[J]. 山地農(nóng)業(yè)生物學(xué)報(bào)，2003，22(3)：204-207.

[4] 呂 蒙，鐘悅之.NDVI在生態(tài)環(huán)境研究中的應(yīng)用進(jìn)展[J]. 農(nóng)村經(jīng)濟(jì)與科技，2011，22(06)：11-14.

[5] 張喜旺，吳炳方.基于中高分辨率遙感的植被覆蓋度時(shí)相變換方法[J]. 生態(tài)學(xué)報(bào)，2015，35(04)：1155-1164.

[6] Gutman G，Ignatov A.The derivation of the green vegetation fraction from NOAA/AVHRR data for use in numerical weather prediction models[J].InternationalJournalofRemoteSensing，1998，19(8)：1533-1543.

[7] 田慶久，閔祥軍.植被指數(shù)研究進(jìn)展[J]. 地球科學(xué)進(jìn)展，1998(04)：10-16.

[8] 顧祝軍，曾志遠(yuǎn).遙感植被蓋度研究[J]. 水土保持研究，2005(02)：18-21.

[9] 孫紅雨，王長(zhǎng)耀，牛 錚，等.中國(guó)地表植被覆蓋變化及其與氣候因子關(guān)系──基于NOAA時(shí)間序列數(shù)據(jù)分析[J]. 遙感學(xué)報(bào)，1998(03)：204-210.

[10] Park H S，Sohn B J.Recent trends in changes of vegetation over East Asia coupled with temperature and rainfall variations[J].JournalofGeophysicalResearchAtmospheres，2010，115(D14)：1307-1314.

[11] Piao S，Wang X，Ciais P，etal. Changes in satellite-derived vegetation growth trend in temperate and boreal Eurasia from 1982 to 2006[J].GlobalChangeBiology，2011，17(10)：3228-3239.

[12] 劉 靜，銀 山，張國(guó)盛，等.毛烏素沙地17年間植被覆蓋度變化的遙感監(jiān)測(cè)[J]. 干旱區(qū)資源與環(huán)境，2009，23(07)：162-167.

[13] 戴聲佩，張 勃，王海軍.中國(guó)西北地區(qū)植被NDVI的時(shí)空變化及其影響因子分析[J]. 地球信息科學(xué)學(xué)報(bào)，2010，12(03)：315-321.

[14] 謝曉華，賴日文，李永實(shí).基于RS技術(shù)的閩江流域植被覆蓋度時(shí)空變化分析[J]. 貴州大學(xué)學(xué)報(bào)：自然科學(xué)版，2008，25(5)：536-539.

[15] 劉亞龍，王 慶，畢景芝，等.基于Mann-Kendall方法的膠東半島海岸帶歸一化植被指數(shù)趨勢(shì)分析[J]. 海洋學(xué)報(bào)(中文版)，2010，32(03)：79-87.

[16] 王新明，王長(zhǎng)耀，牛 錚.應(yīng)用R/S方法分析NDVI時(shí)間序列[J]. 地理與地理信息科學(xué)，2005(05)：20-23，48.

[17] 王文輝，馬祥慶，田 超，等.福建長(zhǎng)汀植被覆蓋度變化的主要驅(qū)動(dòng)影響因子及影響力分析[J]. 福建農(nóng)林大學(xué)學(xué)報(bào)(自然科學(xué)版)，2017，46(03)：277-283.

[18] 袁麗華，蔣衛(wèi)國(guó)，申文明，等.2000～2010年黃河流域植被覆蓋的時(shí)空變化[J]. 生態(tài)學(xué)報(bào)，2013，33(24)：7798-7806.

[19] 靖娟利，王永鋒.西南巖溶區(qū)NDVI時(shí)空變化及其與氣候因子的關(guān)系[J]. 水土保持研究，2016，23(5)：169-174.

[20] 張 君，延軍平.1982-2013年陜西不同植被類型NDVI變化特征分析[J]. 干旱區(qū)資源與環(huán)境，2017，31(04)：86-92.

[21] 張 勃，王 東，王桂鋼，等.西南地區(qū)近14a植被覆蓋變化及其與氣候因子的關(guān)系[J]. 長(zhǎng)江流域資源與環(huán)境，2015，24(6)：956-964.

[22] Stow D，Daeschner S，Hope A，etal. Variability of the Seasonally Integrated Normalized Difference Vegetation Index Across the North Slope of Alaska in the 1990s[J].InternationalJournalofRemoteSensing，2003，24(5)：1111-1117.

[23] 穆少杰，李建龍，陳奕兆，等.2001～2010年內(nèi)蒙古植被覆蓋度時(shí)空變化特征[J]. 地理學(xué)報(bào)，2012，67(09)：1255-1268.

[24] 劉 洋，李誠(chéng)志，劉志輝，等.1982～2013年基于GIMMS-NDVI的新疆植被覆蓋時(shí)空變化[J]. 生態(tài)學(xué)報(bào)，2016，36(19)：6198-6208.

[25] Milich L，Weiss E.GAC NDVI interannual coefficient of variation (CoV) images：ground truth sampling of the Sahel along north-south transects[J].InternationalJournalofRemoteSensing，2000，21(2)：235-260.

[26] 李 卓，孫然好，張繼超，等.京津冀城市群地區(qū)植被覆蓋動(dòng)態(tài)變化時(shí)空分析[J]. 生態(tài)學(xué)報(bào)，2017(22)：1-9.

[27] 韻 婕，任福民，李憶平，等.1960～2010年中國(guó)西南地區(qū)區(qū)域性氣象干旱事件的特征分析[J]. 氣象學(xué)報(bào)，2014，72(2)：266-276.

[28] 韓蘭英，姚玉璧，李憶平，等.近60年中國(guó)西南地區(qū)干旱災(zāi)害規(guī)律與成因[J]. 地理學(xué)報(bào)，2014，69(5)：632-639.

[29] 華 維，范廣洲，李洪權(quán)，等.西南地區(qū)近21年來(lái)NDVI變化特征分析[J]. 成都信息工程學(xué)院學(xué)報(bào)，2008，23(1)：91-97.

[30] 劉世梁，田韞鈺，尹藝潔，等.云南省植被NDVI時(shí)間變化特征及其對(duì)干旱的響應(yīng)[J]. 生態(tài)學(xué)報(bào)，2016，36(15)：4699-4707.

[31] 李仕蓉.2000～2013年貴州省植被覆蓋時(shí)空變化特征與趨勢(shì)[J]. 貴州農(nóng)業(yè)科學(xué)，2014，42(07)：202-205.