周先峰
汶川極重災區(qū)影響植被分布因子是多元的。選取2002-2020 年MYD13Q1-NDVI 數(shù)據(jù),利用地理加權回歸模型分析氣候和其他四個因子的植被變化驅(qū)動力。研究結(jié)果表明:GWR 模型認為年均氣溫對NDVI 的影響最大,其次是年降水、道路歐氏距離、水系歐氏距離。
全球環(huán)境日益漸變并逐漸威脅到人類和社會的發(fā)展,植被作為大氣環(huán)境、土壤環(huán)境和水分循環(huán)的關鍵連接點,植被是自然變化的“顯示器”,研究植被變化趨勢,并探討植被變化的驅(qū)動因子。研究區(qū)地理位置特殊,自然災害頻發(fā),本文利用2002-2020 年MODIS-NDVI數(shù)據(jù)分析極重災區(qū)植被覆蓋變化,從氣候和其他四個因子分析NDVI 與各個因子的耦合關系,分析植被變化的驅(qū)動力。
研究區(qū)位于屬于青藏高原與四川盆地交界處,本文以汶川極重災區(qū)青川縣、平武縣以及汶川縣等共十個縣市為研究區(qū)(圖1),東西橫跨約207 km,南北長約260 km,所轄國土面積約26090km2。

圖1 研究區(qū)位置示意圖
本文采用中分辨率成像光譜儀MYD13Q1-NDVI 數(shù)據(jù),收集2002-2020 年16d 合成產(chǎn)品,利用MRT 軟件進行預處理。從中國氣象網(wǎng)下載2000-2019 年的氣溫、降水數(shù)據(jù)。從全國地理信息資源目錄服務系統(tǒng)獲取2017年水系和道路數(shù)據(jù)生成250m 空間分辨率柵格數(shù)據(jù)。
地理加權回歸(GWR)它是將空間地理位置作為參數(shù)的一種普通線性回歸,為消除數(shù)據(jù)噪音影響,采用最大值合成法合成數(shù)據(jù)。本文應用GWR4 軟件運行GWR。Brunsdon 等對GWR 的原理和算法作了較為詳細的介紹,GWR 模型如公式1 所示:

氣候因子具有長期性和累積性效應,年均氣溫53.61%的區(qū)域呈正相關;年降水22.26%的區(qū)域呈正相關(圖2)。降水較高區(qū)域?qū)DVI 變化影響較高,較低區(qū)域?qū)DVI 變化影響較低。

圖2 氣候因子地理加權回歸估計系數(shù)


圖3 其他因子地理加權回歸估計系數(shù)
道路歐氏距離63.79%的區(qū)域呈負相關,道路主要對植被改善總體產(chǎn)生抑制作用;水系84.85%的區(qū)域呈負相關(圖3),研究區(qū)的水系和水資源對人類活動產(chǎn)生重要的影響。
本文基于MODIS-NDVI 數(shù)據(jù)四個因子,分析2002-2020 年極重災區(qū)十個縣、市NDVI 與影響因子之間的關系,得到以下結(jié)論:GWR 模型認為年均氣溫對NDVI 的影響最大,其次是降水、道路、水系。氣候因子中,對NDVI的影響隨著降水量的增加而增大,介于降水量的最值之間氣溫對NDVI 的影響最大;其他因子中,道路和水系主要對植被改善總體上是抑制作用,人類活動對NDVI 的影響不可忽視。