劉寶柱, 方秀琴, 何祺勝, 榮祁遠(yuǎn)

(河海大學(xué)地球科學(xué)與工程學(xué)院，南京　210098)

?

基于MODIS數(shù)據(jù)和BFAST方法的植被變化監(jiān)測

劉寶柱, 方秀琴, 何祺勝, 榮祁遠(yuǎn)

(河海大學(xué)地球科學(xué)與工程學(xué)院，南京210098)

摘要：植被是聯(lián)結(jié)土壤、大氣和水分的自然“紐帶”，在全球氣候變化研究中具有“指示器”的作用。對歸一化植被指數(shù)(normalizeddifferencevegetationindex，NDVI)時間序列分析，可以為相關(guān)部門的工作和決策提供更好的支持。使用MODISNDVI數(shù)據(jù)結(jié)合BFAST(breaksforadditiveseasonalandtrend)方法實現(xiàn)對老哈河流域及周邊地區(qū)的植被變化監(jiān)測，并確定其NDVI時間序列出現(xiàn)突變點的時間節(jié)點。結(jié)合氣象數(shù)據(jù)以及數(shù)據(jù)本身的質(zhì)量作為影響因子，分析出現(xiàn)突變點的主要原因。研究結(jié)果表明，降水量、相對濕度、溫度、日照時數(shù)、流域蒸發(fā)量與NDVI變化趨勢呈正相關(guān)，風(fēng)速與NDVI變化趨勢相關(guān)性很小。降水量對NDVI變化的影響具有滯后性，滯后時間與降水量大小有關(guān)。關(guān)鍵詞：NDVI； 時間序列；BFAST； 變化監(jiān)測； 突變點

0引言

植被是地球系統(tǒng)生物圈的重要組成部分，包括森林、灌叢、草地、農(nóng)田和果園等，既是生態(tài)環(huán)境的重要組成部分，又是維持生態(tài)環(huán)境、發(fā)揮有效生態(tài)效能的功能體，是衡量自然生態(tài)環(huán)境狀況和性質(zhì)的主要指示物[1]。植被是聯(lián)結(jié)土壤、大氣和水分的自然“紐帶”，在全球氣候變化研究中具有“指示器”的作用[2]。植被覆蓋與氣候變化的關(guān)系是全球變化研究的重要內(nèi)容。研究地表植被覆蓋變化和氣候因子之間的關(guān)系，對區(qū)域乃至全球的能量、生物化學(xué)循環(huán)和水循環(huán)以及氣候變化響應(yīng)具有重要意義。

歸一化植被指數(shù)(normalizeddifferencevegetationindex，NDVI)能夠較好地反映植被的代謝強(qiáng)度及其季節(jié)性變化和年際變化特征，因而被廣泛應(yīng)用于植被的監(jiān)測與分類、物候分析、農(nóng)作物估產(chǎn)、植被-氣候關(guān)系分析等方面。利用NDVI時間序列數(shù)據(jù)研究不同類型植被對氣候變化的響應(yīng)，是當(dāng)前植被生態(tài)系統(tǒng)研究的熱點之一。最近幾年提出的BFAST(breaksforadditiveseasonalandtrend)方法能夠很好地監(jiān)測出NDVI時間序列中的季節(jié)和趨勢突變，特別是能夠指示突變發(fā)生的時空分布[3]。

老哈河流域是我國北方典型的半干旱區(qū)。本研究基于250m的MODISNDVI的時間序列數(shù)據(jù)，利用BFAST方法，監(jiān)測該區(qū)域的NDVI變化狀況，檢測突變信息，結(jié)合氣象數(shù)據(jù)分析突變的原因，以及產(chǎn)生的結(jié)果，為老哈河流域的生態(tài)安全評價、區(qū)域環(huán)境規(guī)劃以及生態(tài)建設(shè)提供參考依據(jù)。

1研究區(qū)概況及數(shù)據(jù)源

本文選取老哈河流域作為研究區(qū)，其范圍如圖1所示。老哈河位于內(nèi)蒙古自治區(qū)東部，是西拉沐淪河的一級支流，流域面積約為1.9×104km2，位于N41°～43°，E117°～120°，流域內(nèi)地形較復(fù)雜，流域出口處高程約400m，源頭的七老圖山海拔2 000m左右。該流域?qū)僦袦貛О敫珊荡箨懶詺夂?，干燥少雨多風(fēng)沙，年降雨量和年平均氣溫分別從西北部山地的300mm和2 ℃變化到東南部低地的600mm和7 ℃。該地區(qū)的農(nóng)作物生長期為5—10月，農(nóng)耕期為4—11月[4]。老哈河流域下墊面多為黃土丘陵區(qū)，地表植被覆蓋率較低，水土流失較為嚴(yán)重。

圖1　研究區(qū)范圍

研究選用的遙感數(shù)據(jù)是2000年2月—2004年12月間Terra衛(wèi)星MOD13Q1的16d數(shù)據(jù)，空間分辨率為250m。這個數(shù)據(jù)可以頻繁地提供由人為因素造成的土地覆蓋變化信息[5]。MOD13Q1的16d數(shù)據(jù)是用受限視場角的最大NDVI值合成的[6]。影像在提取波段信息前，已經(jīng)經(jīng)過精幾何糾正和輻射校正。由于數(shù)據(jù)量較大，使用MRT(MODISreprojectiontool)軟件對MOD13Q1數(shù)據(jù)做拼接、裁剪、重投影以及NDVI波段提取的批處理。所用氣象數(shù)據(jù)來自區(qū)內(nèi)6個氣象站點(圖1)。選取的氣象數(shù)據(jù)，包括日照時數(shù)、相對濕度、平均溫度、降水量、風(fēng)速等主要用于與NDVI的突變結(jié)果做對比分析。其中NDVI值的提取是根據(jù)氣象站點的經(jīng)緯度信息，以3像元×3像元范圍作為興趣區(qū)，取其平均NDVI值作為該氣象站點的NDVI值，以消除個別像元可能存在的偶然誤差。

2研究方法

2.1突變點檢測

流程如圖2所示。根據(jù)BFAST方法對老哈河流域6個氣象站點處于2000—2004年年底期間NDVI突變點的檢測，并結(jié)合氣象數(shù)據(jù)分析其突變點產(chǎn)生的原因及最大影響因子。BFAST方法具體使用的是R語言開發(fā)的BFAST程序包對得到的NDVI數(shù)據(jù)進(jìn)行突變點檢測。BFAST程序包網(wǎng)址見http: //r-forge.r-project.org/R/?group_id=533。

圖2　突變點檢測流程圖

2.2分解模型

BFAST方法將集成的時間序列分解成趨勢、季節(jié)性和殘差3個部分[7]，根據(jù)一定的方法來檢測和表征趨勢和季節(jié)性部分中突然變化的部分。假設(shè)一個加法的分解模型可以迭代出與趨勢和季節(jié)性相匹配的分段線性模型。該模型的算法形式為[8]

Yt=Tt+St+et,t=1,…,n,

(1)

Tt=αj+βjt,

(2)

季節(jié)趨勢[9]可以表達(dá)為

(3)

式中： s為季節(jié)性的時長(如每年的觀察的次數(shù))； γi,j為第i個時段的影響因子。當(dāng)時間t為第i次觀測值時，dt,i=1，否則為0。因此，當(dāng)t為季節(jié)性第0個時段時，dt,i-dt,0=-1。對于t為其他季節(jié)性時段時，dt,i-dt,0=1。dt,i通常被認(rèn)為是季節(jié)性虛擬變量，它有2個可用的值為0和1來解釋回歸模型中的季節(jié)性時段。

2.3迭代算法

在擬合分段線性模型和估算突變點之前，先測試是否有突變點出現(xiàn)在時間序列中。基于移動總和MOSUM(movingsum)測試的普通最小二乘OLS(ordinaryleastsquares)殘差被用來測試是否存在一個或多個突變點。如果測試顯著性變化(P<0.05),則使用Bai和Perron的方法來估算突變點，并根據(jù)貝葉斯信息準(zhǔn)則確定突變點的個數(shù)，以及突變點出現(xiàn)的時間和置信區(qū)間[10]。

2)趨勢系數(shù)αj,βj在區(qū)間j=1,…,m上，由基于M估算法的強(qiáng)健的回歸方程式(2)計算得到。趨勢估算為

(4)

4)季節(jié)性系數(shù)γi,j(i=1,…,s-1; j=1,…,m)可以從基于M估算方法的回歸方程中計算得出。季節(jié)性估算公式為

(5)

3結(jié)果與分析

3.1突變點的檢測結(jié)果

根據(jù)BFAST方法檢測突變點，6個氣象站點中發(fā)現(xiàn)赤峰氣象站和葉柏壽氣象站于2000—2004年年底期間出現(xiàn)季節(jié)性突變點，突變檢測結(jié)果分別如圖3所示。

(a) 赤峰氣象站 (b) 葉柏壽氣象站

圖3赤峰氣象站和葉柏壽氣象站NDVI時間序列分析

Fig.3ChifengmeteorologicalstationandYebaishoumeteorologicalstationNDVItimeseriesanalysis

輸出的突變點檢測結(jié)果見表1。表中括號內(nèi)數(shù)值為第i次觀測時間節(jié)點1a內(nèi)觀測的總次數(shù)。

表1　突變點檢測結(jié)果

從圖3(a)赤峰氣象站NDVI時間序列分析中可以看出，在季節(jié)性成分中出現(xiàn)了2個突變點，而在趨勢中沒有出現(xiàn)突變點，NDVI整體趨勢呈上升趨勢，但是上升趨勢并不明顯，其趨勢斜率為： 0.000 008 18。從整體數(shù)據(jù)可以看出，該點的植被生物量很低。從圖3(b)中可以看出，葉柏壽氣象站處的NDVI的季節(jié)性成分在2000年初至2004年底的時間序列中出現(xiàn)了2個突變點，趨勢成分中并沒有突變點的出現(xiàn)，且趨勢變化不大，其趨勢斜率為0.000 021 51，NDVI線性趨勢擬合始終維持在0.3左右。

3.2氣象數(shù)據(jù)與NDVI突變數(shù)據(jù)的對比分析

3.2.1赤峰氣象站數(shù)據(jù)對比分析

圖4為赤峰氣象站發(fā)生第一個季節(jié)性突變點處的NDVI趨勢圖和對應(yīng)時間點的日照時數(shù)、相對濕度、氣溫、降水量以及風(fēng)速圖的組合，便于進(jìn)行赤峰氣象站季節(jié)性成分突變點位于時間節(jié)點2002(21)數(shù)據(jù)的對比分析。

(a) 赤峰氣象站第一個突變點 (b) 日照時數(shù)

(c) 相對濕度(d) 平均氣溫

(e) 降水量(f) 風(fēng)速

圖4赤峰氣象站季節(jié)性成分突變點位于時間節(jié)點2002(21)的數(shù)據(jù)對比分析

Fig.4Chifengmeteorologicalstationseasonalcomponentbreakpointsdatacomparisonandanalysisat2002(21)

從圖4可以看出，風(fēng)速、降水量、日照時數(shù)、相對濕度與NDVI之間沒有明確的相關(guān)關(guān)系，僅氣溫與NDVI變化呈正相關(guān)關(guān)系，所以氣溫是導(dǎo)致NDVI時間序列季節(jié)性成分上產(chǎn)生突變信息主要因素。

圖5為赤峰氣象站發(fā)生第二個季節(jié)性突變點處的NDVI趨勢圖和對應(yīng)時間點的日照時數(shù)、相對濕度、平均氣溫、降水量以及風(fēng)速圖的組合，便于進(jìn)行赤峰氣象站季節(jié)性成分突變點位于時間節(jié)點2003(13)數(shù)據(jù)的對比分析。

(a) 赤峰氣象站第二個突變點 (b) 日照時數(shù)

(c) 相對濕度 (d) 平均氣溫

圖5-1赤峰氣象站季節(jié)性成分突變點位于時間節(jié)點2003(13)的數(shù)據(jù)對比分析

Fig.5-1Chifengmeteorologicalstationseasonalcomponentbreakpointsdatacomparisonandanalysisat2003(13)

(e) 降水量 (f) 風(fēng)速

圖5-2赤峰氣象站季節(jié)性成分突變點位于時間節(jié)點2003(13)的數(shù)據(jù)對比分析

Fig.5-2Chifengmeteorologicalstationseasonalcomponentbreakpointsdatacomparisonandanalysisat2003(13)

從圖5中可以看出，突變點處的NDVI值與下一個時間點的NDVI相差約0.2，說明是由突變點處的氣象變化而產(chǎn)生的。同樣，風(fēng)速與NDVI之間沒有明確的相關(guān)關(guān)系，從突變點的對應(yīng)時間點看，相對濕度和降水量發(fā)生了較大的變化。因此，該處的季節(jié)性成分突變點主要是由相對濕度和降水量的變化造成的。因為降水量、相對濕度對植物的光合作用速率、生長發(fā)育速度、物候期出現(xiàn)的早晚都有很大的影響。而植被的生長快慢和物候期的提前到來都會對季節(jié)性成分中的突變做出貢獻(xiàn)。

3.2.2葉柏壽氣象站數(shù)據(jù)對比分析

圖6為葉柏壽氣象站發(fā)生第一個季節(jié)性突變點處的NDVI趨勢圖和對應(yīng)時間點的日照時數(shù)、相對濕度、平均溫度、降水量以及風(fēng)速圖的組合，便于進(jìn)行葉柏壽氣象站季節(jié)性成分突變點位于時間節(jié)點2000(23)數(shù)據(jù)的對比分析。

(a) 葉柏壽氣象站第一個突變點 (b) 日照時數(shù)

(c) 相對濕度 (d) 平均氣溫

(e) 降水量 (f) 風(fēng)速

圖6葉柏壽氣象站季節(jié)性成分突變點位于時間節(jié)點2000(23)的數(shù)據(jù)對比分析

Fig.6Yebaishoumeteorologicalstationseasonalcomponentbreakpointsdatacomparisonandanalysisat2000(23)

從圖6中可以得出，平均溫度、降水量與NDVI變化呈正相關(guān)關(guān)系，而日照時數(shù)、相對濕度、風(fēng)速與NDVI相關(guān)性不明顯。從趨勢線中可以看出，平均溫度的趨勢線與NDVI趨勢線非常相似，則平均溫度為NDVI在20001001—20001230日期間變化的主要影響因子。植物蒸發(fā)蒸騰量通常主要包括植物蒸騰量和棵間土壤蒸發(fā)量。而平均溫度的降低導(dǎo)致植物的蒸騰量下降[11]。植物的蒸騰對植物的生長是有幫助的，蒸騰量越大，其生長速率越快。

圖7為葉柏壽氣象站發(fā)生第二個季節(jié)性突變點處的NDVI趨勢圖和對應(yīng)時間點的日照時數(shù)、相對濕度、平均溫度、降水量以及風(fēng)速圖的組合，便于進(jìn)行葉柏壽氣象站季節(jié)性成分突變點位于時間節(jié)點2001(15)數(shù)據(jù)的對比分析。

(a) 葉柏壽氣象站第二個突變點 (b) 日照時數(shù)

(c) 相對濕度 (d) 平均氣溫

(e) 降水量 (f) 風(fēng)速

圖7葉柏壽氣象站季節(jié)性成分突變點位于時間節(jié)點2001(15)的數(shù)據(jù)對比分析

Fig.7Yebaishoumeteorologicalstationseasonalcomponentbreakpointsdatacomparisonandanalysisat2001(15)

從圖7中可以看出，在2001年8月13日之前降水頻繁，降水量最高達(dá)到45mm，距突變時間點最近的一次降水發(fā)生2001年8月7日。降水量對NDVI的影響滯后了至少6d，降水量過大，導(dǎo)致物候期的提前到來，從而產(chǎn)生了季節(jié)性突變點。

3.3源數(shù)據(jù)精度分析

季節(jié)性突變點的產(chǎn)生很有可能是因為所采用的MODIS數(shù)據(jù)的本身精度就不夠高，或者源數(shù)據(jù)中有云的存在等等，導(dǎo)致突變點的產(chǎn)生不完全是因為氣候變化導(dǎo)致的。在ENVI中提取赤峰氣象站和葉柏壽氣象站處的數(shù)據(jù)精度信息，得到赤峰氣象站時間節(jié)點為2002(21)突變點的VI質(zhì)量如表2所示。

表2　2002(21)突變點的VI質(zhì)量

VI(2 189)： 很有可能有云像素，VI有用性為最低質(zhì)量，氣溶膠含量低，不是云邊，沒有進(jìn)行大氣二向性校正，不是混合云，檢測區(qū)域為陸地，無冰和雪覆蓋，無陰影。

VI(2 266)： 生產(chǎn)VI數(shù)據(jù)，但是由其他數(shù)據(jù)對其檢驗的，VI有用性不是很高，氣溶膠含量高，不是云邊，未進(jìn)行大氣二向性校正，無混合云，檢測區(qū)域為陸地，無冰和雪覆蓋，無云陰影。

VI(18 508)： 產(chǎn)生VI，且質(zhì)量好，VI有用性為最低質(zhì)量，氣溶膠含量為中，不是云邊，未進(jìn)行大氣二向性校正，不是混合云，檢測區(qū)域為暫時性有水區(qū)域，有冰和雪覆蓋，無云陰影。

從上面的數(shù)據(jù)質(zhì)量結(jié)果可以看出，VI有用性普遍都很低，且不是云邊，未進(jìn)行大氣二向性校正，不是混合云，在時間節(jié)點為2002(19)和2002(21)的檢測區(qū)域為陸地，而在2002(22)時間節(jié)點的檢測區(qū)域為暫時性水域，說明該處發(fā)生的突變情況主要由降雨造成； 此外，3個時間節(jié)點均未被冰或雪覆蓋。結(jié)合之前圖表對比分析，可以得出結(jié)論： 該季節(jié)性突變點出現(xiàn)的原因主要是由降雨造成的。

赤峰氣象站時間節(jié)點為2003(13)突變點的VI質(zhì)量如表3所示。

表3　2003(13)突變點的VI質(zhì)量

VI(2 257)： 產(chǎn)生的數(shù)據(jù)很有可能有云像素，VI有用性為質(zhì)量較低，氣溶膠含量較低，不是云邊，未進(jìn)行大氣二向性校正，無混合云，檢測區(qū)域為陸地，無冰和雪覆蓋，無陰影。

VI(34 842)： 產(chǎn)生VI，但是需要檢驗其他QA，VI有用性為質(zhì)量較低，氣溶膠含量具有氣候性，不是云邊，未進(jìn)行大氣二向性校正，無混合云，檢測區(qū)域為暫時性水域，無冰和雪覆蓋，有云陰影。

VI(2 189)： 很有可能有云像素，VI有用性為最低質(zhì)量，氣溶膠含量低，不是云邊，沒有進(jìn)行大氣二向性校正，不是混合云，檢測區(qū)域為陸地，無冰和雪覆蓋，無陰影。

從上面的數(shù)據(jù)質(zhì)量結(jié)果可以看出，VI有用性普遍都很低，且不是云邊，未進(jìn)行大氣二向性校正，不是混合云，在時間節(jié)點為2003(12)和2003(15)的檢測區(qū)域為陸地，而在2003(13)時間節(jié)點的檢測區(qū)域為暫時性水域，說明該處發(fā)生的突變情況主要是降雨造成，此外，3個時間節(jié)點均未被冰或雪覆蓋，且2003(15)時間節(jié)點處有云陰影，使得數(shù)據(jù)質(zhì)量較低，但是突變點為2003(13)，影響不大。結(jié)合之前圖表對比分析，可以得出結(jié)論： 該季節(jié)性突變點出現(xiàn)的原因主要是由降雨造成的。

葉柏壽氣象站時間節(jié)點為2000(23)突變點的VI質(zhì)量見表4。

表4　2000(23)突變點的VI質(zhì)量

VI(2 120)： 產(chǎn)生VI，質(zhì)量好，VI有用性為質(zhì)量較低，氣溶膠含量為中，不是云邊，未進(jìn)行大氣二向性校正，無混合云，檢測區(qū)域為陸地，無冰/雪覆蓋，無云陰影。

從上面質(zhì)量檢查結(jié)果可以看出，該時間節(jié)點的數(shù)據(jù)質(zhì)量較好，產(chǎn)生突變點的情況與數(shù)據(jù)本身無關(guān)，而從之前的圖表對比分析中得出，產(chǎn)生該突變的原因主要是與流域蒸發(fā)量有關(guān)?？梢缘贸鼋Y(jié)論為： 2000(23)時間節(jié)點處的季節(jié)性突變主要是因為流域蒸發(fā)量的下降。

葉柏壽氣象站時間節(jié)點為2001(15)突變點的VI質(zhì)量見表5。

表5　2001(15)突變點的VI質(zhì)量

VI(2 116)： 產(chǎn)生VI，質(zhì)量好，VI有用性為質(zhì)量較低，氣溶膠含量為中，不是云邊，未進(jìn)行大氣二向性校正，無混合云，檢測區(qū)域為陸地，無冰/雪覆蓋，無云陰影。

從上面質(zhì)量檢查結(jié)果可以看出，該時間節(jié)點的數(shù)據(jù)質(zhì)量較好，產(chǎn)生突變點的情況與數(shù)據(jù)本身無關(guān)，而從之前的圖表對比分析中得出，產(chǎn)生該突變的原因主要是與降水量有關(guān)。降水量最大值為45mm左右，且降水量對NDVI的影響具有滯后性，滯后時間為6d。所以可以得出結(jié)論為： 2001(15)時間節(jié)點處的季節(jié)性突變主要是因為降水量突然增多原因造成的。

4結(jié)論與展望

本文使用了BFAST方法將NDVI時間序列數(shù)據(jù)分成3個部分，分別為季節(jié)性成分、趨勢成分和殘差，并對季節(jié)性成分和趨勢成分做回歸分析，檢測并記錄突變點的信息； 結(jié)合氣象數(shù)據(jù)分析產(chǎn)生突變點的原因并得出結(jié)論。目前國內(nèi)較少見做該方面研究的報道。本文結(jié)合了氣象數(shù)據(jù)對突變點產(chǎn)生的原因進(jìn)行分析，更加詳細(xì)地解釋了產(chǎn)生突變的驅(qū)動因素和變化機(jī)制。研究得出以下結(jié)論：

1)在2000年初至2004年底，對6個氣象站點處的NDVI時間序列突變進(jìn)行了監(jiān)測，分別在赤峰氣象站監(jiān)測到季節(jié)性突變點，時間節(jié)點為2002(21)和2003(13)以及葉柏壽氣象站監(jiān)測季節(jié)性突變點，時間節(jié)點為2000(23)和2001(15)。

2)季節(jié)性突變點2002(21)產(chǎn)生的原因主要是由降雨造成的； 季節(jié)性突變點2003(13)產(chǎn)生的原因主要是由降雨造成的； 季節(jié)性突變點2000(23)產(chǎn)生的原因主要是由流域蒸發(fā)量下降造成； 季節(jié)性突變點2001(15)產(chǎn)生的原因主要是由降水量突然增多造成。

3)在2000年初至2004年底出現(xiàn)突變點的原因與數(shù)據(jù)本身質(zhì)量無關(guān)，數(shù)據(jù)本身質(zhì)量較好。

4)NDVI變化趨勢與降水量、流域蒸發(fā)量、氣溫、相對濕度、日照時數(shù)呈正相關(guān)趨勢，與風(fēng)速相關(guān)性不大。降水量對NDVI變化的影響具有滯后性，滯后時間與降水量大小有關(guān)。

5)不同時期NDVI變化與降水量、流域蒸發(fā)量、氣溫、相對濕度、日照時數(shù)及風(fēng)速等相關(guān)性都不一樣，影響程度也不同。

本文研究可以為自然和農(nóng)業(yè)土地生產(chǎn)評估、農(nóng)作物狀況以及物候變化研究提供參考依據(jù)。本研究還存在一些不足之處，只對氣象站點處的突變信息進(jìn)行提取和分析，尚未實現(xiàn)整個面上的突變信息的提取和分析，在后續(xù)研究中，會結(jié)合全球同化數(shù)據(jù)對面上的突變點進(jìn)行分析。

參考文獻(xiàn)(References)：

[1]張飛,塔西甫拉提·特依拜,丁建麗,等.新疆典型鹽漬區(qū)植被覆蓋度遙感動態(tài)監(jiān)測——以渭干河-庫車河三角洲綠洲為例[J].林業(yè)科學(xué),2011,47(7):27-35.

ZhangF,TashpolatT,DingJL,etal.DynamicallymonitoringvegetationcoverbyremotesensinginthetypicalsalinizationregionofXinjiang:AcasestudyindeltaoasisofWeiganandKuqaRivers[J].ScientiaSilvaeSinicae,2011,47(7):27-35.

[2]劉蕾.2001—2005年新疆植被覆蓋動態(tài)變化及原因分析[J].干旱環(huán)境監(jiān)測,2007,21(3):146-148.

LiuL.AnalysisonXinjiangvegetationcoveragedynamicchangesanditscausesfrom2001to2005[J].AridEnvironmentalMonitoring,2007,21(3):146-148.

[3]VerbesseltJ,HyndmanR,NewnhamG,etal.Detectingtrendandseasonalchangesinsatelliteimagetimeseries[J].RemoteSensingofEnvironment,2010,114(1):106-115.

[4]方秀琴,任立良.西遼河的老哈河流域土地利用遙感動態(tài)監(jiān)測[J].地球信息科學(xué)學(xué)報,2009,11(1):125-131.

FangXQ,RenLL.DetectionoflandusechangeinLaohaheRiverBasinofWestLiaoheRiverbasedonremote-sensingimages[J].JournalofGeo-InformationScience,2009,11(1):125-131.

[5]TownshendJRG,JusticeCO.Selectingthespatialresolutionofsatellitesensorsrequiredforglobalmonitoringoflandtransformations[J].InternationalJournalofRemoteSensing,1988,9(2):187-236.

[6]HueteA,DidanK,MiuraT,etal.OverviewoftheradiometricandbiophysicalperformanceoftheMODISvegetationindices[J].RemoteSensingofEnvironment,2002,83(1/2):195-213.

[7]HueteA,JusticeC,vanLeeuwenW.ModisVegetationIndex(MOD13):AlgorithmTheoreticalBasisDocument[EB/OL].(1999-04-30).http://eospso.gsfc.nasa.gov/ftp_ATBD/REVIEW/MODIS/ATBD-MOD-13/atbd-mod-13.pdf.

[8]VerbesseltJ,HyndmanR,ZeileisA,etal.Phenologicalchangedetectionwhileaccountingforabruptandgradualtrendsinsatelliteimagetimeseries[J].RemoteSensingofEnvironment,2010,114(12):2970-2980.

[9]LinderholmHW.Growingseasonchangesinthelastcentury[J].AgriculturalandForestMeteorology,2006,137(1/2):1-14.

[10]BaiJS,PerronP.Computationandanalysisofmultiplestructuralchangemodels[J].JournalofAppliedEconometrics,2003,18(1):1-22.

[11]屈艷萍,康紹忠,張曉濤,等.植物蒸發(fā)蒸騰量測定方法述評[J].水利水電科技進(jìn)展,2006,26(3):72-77.

QuYP,KangSZ,ZhangXT,etal.Areviewofmethodsformeasurementofevaportranspirationfromplants[J].AdvancesinScienceandTechnologyofWaterResources,2006,26(3):72-77.

(責(zé)任編輯： 李瑜)

Monitoring the changes of vegetation based on MODIS data and BFAST methods

LIU Baozhu, FANG Xiuqin, HE Qisheng, RONG Qiyuan

( School of Earth Science and Engineering, Hohai University, Nanjing 210098, China)

Abstract:Vegetationisanatural“l(fā)ink”whichlinkssoil,airandwaterandan“indicator”inglobalclimatechangeresearch.Usingnormalizeddifferencevegetationindex(NDVI)time-seriesanalyses,wecanprovidebettersupportfortherelevantresearchesanddecision-making.UsingMODISNDVIdatabindingwithBFAST(breaksforadditiveseasonalandtrend)method,theauthorsimplementedmonitoringvegetationdynamicsintheLaohaheRiverBasinandthesurroundingareas,andidentifieditsNDVItime-seriesabruptchangepointsoccurringintime.Themeteorologicaldataandthequalityofthedataitselfwerealsousedasaninfluencefactoranalysisofthemainreasonforthebreakpoints.Itisfoundthatprecipitation,relativehumidity,temperature,sunshineandwaterevaporationarepositivelycorrelatedwithNDVItrends,whilewindspeedislesscorrelatedwithNDVItrends.What’smore,theprecipitationandsunshinehourimpactonNDVIchangehasacertainlag.

Keywords:NDVI;time-series;BFAST;changemonitoring;breakpoints

doi:10.6046/gtzyyg.2016.03.23

收稿日期：2015-04-14；

修訂日期：2015-08-28

基金項目：水利部公益性行業(yè)科研專項經(jīng)費項目“東北灌區(qū)節(jié)水灌溉生態(tài)與增產(chǎn)效應(yīng)評估研究”(編號： 201401001)、國家自然科學(xué)基金項目“復(fù)雜地形條件下多源遙感數(shù)據(jù)森林生物量協(xié)同反演研究”(編號： 41101308)及“北方半干旱區(qū)典型土地利用變化的水文效應(yīng)”(編號： 41201027)共同資助。

中圖法分類號：TP79

文獻(xiàn)標(biāo)志碼：A

文章編號：1001-070X(2016)03-0146-08

第一作者簡介：劉寶柱(1993-)，男，碩士研究生，主要從事定量遙感研究。Email:bzliu@hhu.edu.cn。

通信作者：何祺勝(1981-)，男，博士，副教授，主要從事定量遙感、GIS建模與開發(fā)等方面的研究。Email:heqis@hhu.edu.cn。

引用格式： 劉寶柱,方秀琴,何祺勝,等.基于MODIS數(shù)據(jù)和BFAST方法的植被變化監(jiān)測[J].國土資源遙感,2016,28(3):146-153.(LiuBZ,FangXQ,HeQS,etal.MonitoringthechangesofvegetationbasedonMODISdataandBFASTmethods[J].RemoteSensingforLandandResources,2016,28(3):146-153.)