蔣宇威 李登秋

摘要 以浙江省杭州市臨安區(qū)森林植被作為研究對(duì)象，選擇火災(zāi)、泥石流和滑坡3種典型森林干擾類型，利用1984—2016年所有可用Landsat遙感影像，采用分段諧波模型擬合方法，對(duì)干擾發(fā)生前后植被的生長狀況進(jìn)行表征，對(duì)不同干擾類型導(dǎo)致的植被生長狀況及干擾后恢復(fù)過程進(jìn)行評(píng)估。結(jié)果表明：3種干擾事件均導(dǎo)致植被歸一化指數(shù)（NDVI）呈現(xiàn)顯著下降，火災(zāi)和泥石流導(dǎo)致植被的變化程度分別為22.7%和20.5%，而滑坡則導(dǎo)致植被的下降程度達(dá)42.7%。災(zāi)后植被均呈現(xiàn)快速恢復(fù)，但過程存在一定差異，火災(zāi)、泥石流和滑坡對(duì)植被所造成的破壞分別需要4、6和5年恢復(fù)到災(zāi)前植被覆蓋水平。時(shí)間序列Landsat數(shù)據(jù)能夠較好地表征干擾植被變化過程，可為監(jiān)測和評(píng)估不同干擾事件對(duì)植被生長的影響提供依據(jù)。

關(guān)鍵詞 地質(zhì)災(zāi)害;森林火災(zāi);植被干擾;NDVI;時(shí)間序列

中圖分類號(hào) S 771.8 ?文獻(xiàn)標(biāo)識(shí)碼 A ?文章編號(hào) 0517-6611（2021）14-0106-04

Abstract Using all available Landsat images between 1984 and 2016， this study selected three typical forest disturbance types （fire， mudslide and landslide） to assess their impacts on vegetation coverage in Linan District， Hangzhou City， Zhejiang Province. The piecewise harmonic model was used to fit the time series Landsat Normalized Difference Vegetation Index （NDVI）， which can characterize the vegetation coverage condition before and after the disturbance. The results showed that the NDVI decreased significantly after the three disturbance events occurred. The change degrees in the fire， mudflow and landslide areas were 22.7%， 20.5% and 42.7%， respectively. Vegetation recovered rapidly after the disasters， but there were some differences in the processes. The areas disturbed by fire， mudflow and landslide needed 4， 6 and 5 years to recover to the predisturbance coverage level， respectively. Dense time series Landsat data can represent the vegetation change process， and provide a basis for evaluating the impact of different disturbances on vegetation coverage.

Key words Geological disasters;Forest fire;Vegetation disturbance;NDVI;Time series

基金項(xiàng)目 國家自然科學(xué)基金（41701490）;國家重點(diǎn)研發(fā)計(jì)劃（2016YFC0503302）。

作者簡介 蔣宇威（1995—），男，浙江臺(tái)州人，從事災(zāi)害評(píng)估與時(shí)間序列遙感數(shù)據(jù)分析研究。*通信作者，副教授，博士，從事時(shí)間序列遙感森林變化監(jiān)測研究。

收稿日期 2020-11-06

干擾普遍存在于森林生態(tài)系統(tǒng)，是森林群落演替的驅(qū)動(dòng)力之一[1]。根據(jù)干擾起因，森林干擾可分為自然干擾和人為干擾2種類型[2]。地質(zhì)災(zāi)害與火災(zāi)是森林自然干擾的重要組成部分，地質(zhì)災(zāi)害是指由于自然或人為作用在地球表層比較強(qiáng)烈地破壞人類生命財(cái)產(chǎn)和生存環(huán)境的巖土體移動(dòng)事件，地質(zhì)災(zāi)害的發(fā)生往往會(huì)導(dǎo)致森林被掩埋毀壞。森林火災(zāi)是指失去人為控制，在林地內(nèi)自由蔓延和擴(kuò)展，直接減少森林面積，破壞森林結(jié)構(gòu)和森林環(huán)境，導(dǎo)致森林生態(tài)系統(tǒng)失衡、森林生物量下降、生產(chǎn)力減弱的自然災(zāi)害。準(zhǔn)確表征森林干擾的破壞程度以及植被恢復(fù)過程，對(duì)于保護(hù)和管理森林生態(tài)系統(tǒng)具有重要意義[3-4]。

隨著遙感技術(shù)的發(fā)展，時(shí)間序列遙感已成為森林動(dòng)態(tài)監(jiān)測的重要手段[5-8]。目前，基于時(shí)間序列遙感數(shù)據(jù)對(duì)植被變化監(jiān)測已有大量研究[9]，如羅火錢等[10]采用多時(shí)相遙感影像評(píng)價(jià)了火燒跡地的森林恢復(fù)情況;梁超等[11]以災(zāi)害前后2期中分辨率成像光譜儀（moderateresolution imaging spectroradiometer，MODIS）歸一化植被指數(shù)（normalized difference vegetation index，NDVI）數(shù)據(jù)提取了災(zāi)害導(dǎo)致的植被破壞區(qū)，發(fā)現(xiàn)泥石流對(duì)植被的破壞程度最嚴(yán)重，崩塌次之，滑坡的破壞程度相對(duì)最低，且植被的破壞程度與海拔、坡度均有關(guān)系。侍昊等[12]基于2005—2008年MODIS NDVI構(gòu)建了森林災(zāi)害評(píng)估體系，以探測森林資源損失的空間分布范圍和災(zāi)害等級(jí)。不同災(zāi)害對(duì)植被干擾程度不同，植被覆蓋度對(duì)災(zāi)后損害響應(yīng)存在滯后性，受損植被恢復(fù)率與水系距離、海拔等相關(guān)，且存在明顯的空間差異[13-14]。長時(shí)間序列遙感能夠持續(xù)地監(jiān)測森林變化特征，為森林生產(chǎn)力與森林碳儲(chǔ)量評(píng)估提供有效的數(shù)據(jù)支撐[15]，但目前相關(guān)研究較多集中于單一的森林干擾，且多以中分辨率成像光譜儀（MODIS）等低分辨率或多期Landsat數(shù)據(jù)為主，對(duì)于我國亞熱帶森林地區(qū)小尺度、破碎化的干擾事件難以進(jìn)行有效監(jiān)測，且僅依靠多期遙感數(shù)據(jù)難以有效表征植被干擾和恢復(fù)過程。筆者基于歷史可用Landsat NDVI數(shù)據(jù)和分段諧波模型，分析滑坡、泥石流和火災(zāi)3種干擾類型對(duì)亞熱帶森林植被的干擾程度以及災(zāi)后植被恢復(fù)過程。

1 材料與方法

1.1 地質(zhì)災(zāi)害與火災(zāi)數(shù)據(jù)

以浙江省杭州市臨安區(qū)為調(diào)查區(qū)域，臨安區(qū)森林覆蓋率達(dá)78.2%，森林資源豐富，且氣候季節(jié)性明顯，是森林火災(zāi)易發(fā)區(qū)域。臨安區(qū)以中低山丘陵為主要地貌，雨量充沛季節(jié)，易發(fā)生滑坡、泥石流等地質(zhì)災(zāi)害，尤以山核桃種植區(qū)清涼峰鎮(zhèn)與昌化鎮(zhèn)為主要易發(fā)區(qū)域?；谝巴庹{(diào)查及歷史數(shù)據(jù)確定3個(gè)典型干擾區(qū)域，包括1個(gè)火災(zāi)區(qū)域、1個(gè)滑坡區(qū)域及1個(gè)泥石流區(qū)域（表1）。結(jié)合野外調(diào)查數(shù)據(jù)、Google Earth以及Landsat歷史影像，明確每個(gè)干擾事件的影響范圍（圖1），用于計(jì)算每期Landsat影像干擾區(qū)域NDVI平均值。

1.2 時(shí)間序列Landsat NDVI數(shù)據(jù)

Landsat衛(wèi)星是由美國航空航天局（NASA）發(fā)射的用于探測地球資源與環(huán)境的系列地球觀測衛(wèi)星系統(tǒng)，自1972年7月23日以來，已發(fā)射8顆衛(wèi)星。該研究通過USGS下載了覆蓋研究區(qū)的1984—2016年云量小于80%的所有Landsat TM/ETM+/OLI數(shù)據(jù)（圖2），共收集1 045幅Landsat L1T影像數(shù)據(jù)（其中，軌道P119/R039 483幅，軌道P118/R041 563幅）。Landsat L1T產(chǎn)品已經(jīng)經(jīng)過大氣校正、幾何校正、云和陰影的檢測，采用式（1）對(duì)每景數(shù)據(jù)計(jì)算NDVI，來表征植被生長狀況。NDVI能夠綜合反映植物生長、覆蓋、生物量和植被種類等情況，通過植被NDVI值的變化，可以分析區(qū)域內(nèi)植被的生長及其變化特征，并能監(jiān)測干擾后的植被恢復(fù)情況。

NDVI=（NIR- R ）/（NIR+ R ）（1）

式中，NIR為近紅外波段反射率， R 為紅波段反射率。

1.3 研究方法

根據(jù)每個(gè)干擾事件的發(fā)生范圍和1984—2016年每期Landsat NDVI數(shù)據(jù)，利用R語言提取出每景影像中干擾區(qū)域的NDVI平均值，采用時(shí)間序列諧波模型對(duì)每個(gè)干擾事件發(fā)生前后時(shí)間序列NDVI進(jìn)行擬合[式（2）]，模型系數(shù)由普通最小二乘法估計(jì)，采用公式（3）～（5）評(píng)估干擾對(duì)森林植被的影響。在時(shí)間序列模型中， βi是表征植被生長的總體特征，ai和bi 代表了物候?qū)W和日照角差異引起的NDVI的年內(nèi)變化， ci 捕獲NDVI的年際變化。

其中， βi和ci 代表截距和年際變化（斜率）系數(shù)， t start， i 和 t end， i 代表干擾前（后）時(shí)間序列NDVI起始點(diǎn)（儒略日）和終止點(diǎn)（儒略日）。NDVIstart， i 和NDVIend， i 代表干擾前（后）時(shí)間序列起始點(diǎn)和終止點(diǎn)剔除季節(jié)影響后的NDVI。 A 表示干擾導(dǎo)致的NDVI凈變化， Ar 表示干擾導(dǎo)致的NDVI相對(duì)變化程度。

2 結(jié)果與分析

2.1 火災(zāi)對(duì)植被的干擾過程

時(shí)間序列Landsat NDVI和諧波模型擬合結(jié)果表征了火災(zāi)對(duì)植被的影響過程（圖3），火災(zāi)發(fā)生前，該區(qū)域以杉木林為主，NDVI處于較穩(wěn)定的季節(jié)波動(dòng)狀態(tài)，擬合曲線的斜率接近于0，表明在未遭受明顯干擾的情況下，植被變化主要表現(xiàn)為季節(jié)波動(dòng)。火災(zāi)后發(fā)生后（2008年），時(shí)間序列NDVI表現(xiàn)出明顯的轉(zhuǎn)折點(diǎn)，NDVI在火災(zāi)發(fā)生后顯著降低，由災(zāi)害前的0.61降低為0.47（表2），植被在受到干擾后恢復(fù)過程明顯，NDVI呈現(xiàn)顯著上升趨勢，且2008—2016年一直處于快速恢復(fù)狀態(tài)，在2012年基本恢復(fù)到干擾前水平，但干擾后植被NDVI的波動(dòng)幅度強(qiáng)于干擾發(fā)生前，主要原因在于火災(zāi)發(fā)生后，部分森林轉(zhuǎn)變?yōu)檗r(nóng)田、苗木，使得植被季節(jié)變化增強(qiáng)。

2.2 地質(zhì)災(zāi)害對(duì)植被的干擾過程

2.2.1 泥石流對(duì)植被的干擾過程。

圖4表征了泥石流對(duì)植被的影響及其恢復(fù)過程，災(zāi)害發(fā)生前，植被覆蓋度相對(duì)較好，NDVI圍繞0.7呈現(xiàn)季節(jié)性波動(dòng)，該地點(diǎn)由于突降特大暴雨造成山體塌方，由此引發(fā)泥石流，泥石流導(dǎo)致NDVI顯著下降，由災(zāi)害前的0.70降低為0.56（表2），與火災(zāi)干擾相類似，泥石流發(fā)生后植被呈現(xiàn)快速恢復(fù)，于2011年左右恢復(fù)至災(zāi)前水平，與災(zāi)害發(fā)生前相比NDVI季節(jié)性波動(dòng)無明顯差異。

2.2.2 山體滑坡對(duì)植被的干擾過程。

圖5為滑坡對(duì)植被的影響過程，該地區(qū)地勢較陡（坡度達(dá)到30～50°），土層厚度達(dá)3～8 m，地質(zhì)環(huán)境條件脆弱;存在人類工程活動(dòng)影響，切坡修路改變了坡面的原始狀況，降低了山體的穩(wěn)定性，加上受“莫拉克”臺(tái)風(fēng)影響導(dǎo)致的持續(xù)降雨，使土體水分處于過飽和狀態(tài)，從而引發(fā)滑坡。相比于火災(zāi)發(fā)生范圍，滑坡發(fā)生區(qū)域面積較小，在Landsat上僅表現(xiàn)為2個(gè)像元，滑坡發(fā)生前NDVI均值處于0.5左右，NDVI呈現(xiàn)一定的增加趨勢，滑坡導(dǎo)致NDVI下降較為顯著，由災(zāi)害前的0.59降低為0.33，變化率達(dá)42.7%（表2），滑坡發(fā)生后植被恢復(fù)速度顯著，災(zāi)害發(fā)生后NDVI于2014年恢復(fù)到災(zāi)前水平。

3 結(jié)論與討論

該研究通過長時(shí)間序列密集遙感數(shù)據(jù)表征植被長期變化及干擾事件對(duì)植被生長過程的影響，基于諧波模型擬合能夠較好地消除時(shí)間序列異常值和季節(jié)波動(dòng)對(duì)植被變化監(jiān)測的不利影響，通過對(duì)火災(zāi)、泥石流和滑坡3種典型干擾事件影響下的時(shí)間序列NDVI分析，發(fā)現(xiàn)3種干擾事件均導(dǎo)致NDVI顯著降低，其中火災(zāi)和泥石流導(dǎo)致植被的變化程度分別為22.7%和20.5%，滑坡則導(dǎo)致植被的下降程度達(dá)42.7%。災(zāi)后植被均呈現(xiàn)快速的恢復(fù)，但恢復(fù)過程存在一定差異，火災(zāi)、泥石流和滑坡對(duì)植被所造成的破壞分別需要4、6和5年恢復(fù)到災(zāi)前植被覆蓋水平。

該研究中火災(zāi)所導(dǎo)致的干擾面積較大，而泥石流和滑坡則僅發(fā)生在有限的幾個(gè)像元內(nèi)，且存在一定的混合像元，可能會(huì)低估地質(zhì)災(zāi)害對(duì)植被的破壞能力。泥石流發(fā)生的區(qū)域，雖然NDVI表明植被覆蓋良好，但由于處于山核桃種植區(qū)，林下植被多被清理，水土流失嚴(yán)重，是發(fā)生泥石流的重要原因。在火災(zāi)發(fā)生區(qū)域，部分像元災(zāi)后土地利用類型發(fā)生了變化，由森林轉(zhuǎn)為農(nóng)田、苗木和建設(shè)用地，使得整個(gè)火災(zāi)區(qū)域的植被恢復(fù)更為復(fù)雜，更為詳細(xì)地評(píng)估火災(zāi)對(duì)植被的影響需要在像元尺度開展。此外，部分像元干擾發(fā)生后，NDVI呈現(xiàn)持續(xù)下降的趨勢，表明干擾發(fā)生后植被呈現(xiàn)持續(xù)惡化，這些特征需在具體的災(zāi)害評(píng)估中進(jìn)一步明確其影響因素。

參考文獻(xiàn)

[1] 于立忠，朱教君，張艷紅，等.森林干擾度評(píng)價(jià)[J].生態(tài)學(xué)雜志，2009，28（5）：976-982.

[2] 高偉，韓孟孟，辛秀，等.干擾理論及其對(duì)森林資源的影響[J].中國城市林業(yè)，2009，7（5）：16-18.

[3] 劉孝富，王文杰，李京，等.災(zāi)后生態(tài)恢復(fù)評(píng)價(jià)研究進(jìn)展[J].生態(tài)學(xué)報(bào)，2014，34（3）：527-536.

[4] MENG Y Y，LIU X N，WU L，et al.Spatiotemporal variation indicators for landscape structure dynamics monitoring using dense normalized difference vegetation index time series[J/OL].Ecological indicators，2019，107[2020-04-28].https：//doi.org/10.1016/j.ecolind.2019.105607.

[5] 李洛晞，沈潤平，李鑫慧，等.基于MODIS時(shí)間序列森林?jǐn)_動(dòng)監(jiān)測指數(shù)比較研究[J].遙感技術(shù)與應(yīng)用，2016，31（6）：1083-1090.

[6] MUNROE D K，SOUTHWORTH J，TUCKER C M.Modeling spatially and temporally complex landcover change：The case of western Honduras[J].Professional geographer，2004，56（4）：544-559.

[7] KENNEDY R E，COHEN W B，SCHROEDER T A.Trajectorybased change detection for automated characterization of forest disturbance dynamics[J].Remote sensing of environment，2007，110（3）：370-386.

[8] POTAPOV P V，TURUBANOVA S A，HANSEN M C，et al.Quantifying forest cover loss in Democratic Republic of the Congo，2000-2010，with Landsat ETM + data[J].Remote sensing of environment，2012，122：106-116.

[9] ZHU Z，ZHANG J X，YANG Z Q，et al.Continuous monitoring of land disturbance based on Landsat time series[J/OL].Remote sensing of environment，2020，238[2020-04-28].https：//doi.org/10.1016/j.rse.2019.03.009.

[10] 羅火錢，陳文惠.火燒跡地森林恢復(fù)的植被指數(shù)時(shí)間序列特征譜研究[J].四川林勘設(shè)計(jì)，2012（1）：23-30.

[11] 梁超，趙廷寧，史常青，等.基于NDVI的汶川大地震前后北川縣次生地質(zhì)災(zāi)害區(qū)植被破壞評(píng)估[J].中國水土保持科學(xué)，2013，11（4）：86-92.

[12] 侍昊，王笑，薛建輝，等.基于MODIS/NDVI時(shí)間序列的森林災(zāi)害快速評(píng)估方法：以貴州省為例[J].生態(tài)學(xué)報(bào)，2012，32（11）：3359-3367.

[13] 李京忠，曹明明，邱海軍，等.汶川地震區(qū)災(zāi)后植被恢復(fù)時(shí)空過程及特征：以都江堰龍溪河流域?yàn)槔齕J].應(yīng)用生態(tài)學(xué)報(bào)，2016，27（11）：3479-3486.

[14] 王愛愛，臧淑英，王翠珍，等.重建NDVI時(shí)間序列及火后森林恢復(fù)時(shí)空動(dòng)態(tài)分析[J].哈爾濱師范大學(xué)自然科學(xué)學(xué)報(bào)，2017，33（4）：54-61.

[15] 沈文娟，李明詩.基于長時(shí)間序列Landsat影像的南方人工林干擾與恢復(fù)制圖分析[J].生態(tài)學(xué)報(bào)，2017，37（5）：1438-1449.