賴(lài)云森 束志文 牛 地 周玲美

安徽省地質(zhì)礦產(chǎn)勘查局332地質(zhì)隊(duì)，安徽 黃山 245000

0 引言

植被通常是指地球表面某一區(qū)域所覆蓋的植物群落，在調(diào)節(jié)氣候和維持生態(tài)平衡方面發(fā)揮著關(guān)鍵作用，是陸地生態(tài)系統(tǒng)的重要組成部分［1-2］。植被作為陸地生態(tài)系統(tǒng)對(duì)全球氣候變化響應(yīng)的重要“指示器”［3-4］，探究植被覆蓋度的時(shí)空變化對(duì)全球氣候變化的響應(yīng)，可為應(yīng)對(duì)全球氣候變化提供理論依據(jù)［5］，在全球氣候變化研究中具有重要意義［6］。歸一化植被指數(shù)（Normal Difference Vegetation Index，NDVI）作為表征植被生長(zhǎng)狀態(tài)和植被覆蓋度的最佳指標(biāo)［7］，與植被分布密度呈線性相關(guān)［8-9］。近年來(lái)，越來(lái)越多的學(xué)者基于NDVI 數(shù)據(jù)探究不同區(qū)域及尺度的植被覆蓋變化［2］。如張玉婷等［10］基于MODIS-NDVI 數(shù)據(jù)，利用最大值合成法、均值法和趨勢(shì)分析法對(duì)2011—2020 年甘肅省武威市植被覆蓋度時(shí)空變化進(jìn)行分析，發(fā)現(xiàn)2011—2020 年武威市植被NDVI整體呈增加趨勢(shì)；涂又等［11］基于AVHRR NDVI3g 遙感數(shù)據(jù)探究了1982—2015 年中國(guó)植被NDVI 時(shí)空變化及驅(qū)動(dòng)因子，結(jié)果表明大多數(shù)植被呈持續(xù)穩(wěn)定變化狀態(tài)，且中國(guó)植被變化主要受到氣候變化和人類(lèi)活動(dòng)的共同驅(qū)動(dòng)；劉珊珊等［12］采用2001—2015 年MODIS-NDVI 數(shù)據(jù)分析了云南省植被覆蓋變化，結(jié)果顯示云南省植被覆蓋度呈現(xiàn)基本穩(wěn)定或增加的趨勢(shì)。安徽省位于我國(guó)華東地區(qū)，地形地貌復(fù)雜，隨著經(jīng)濟(jì)的不斷發(fā)展及城鎮(zhèn)化的快速發(fā)展，植被覆蓋也發(fā)生著巨大的變化，對(duì)人類(lèi)社會(huì)和生態(tài)環(huán)境的影響愈發(fā)明顯［1］。因此，筆者基于MODIS-NDVI 數(shù)據(jù)分析2000—2021 年安徽省NDVI 的時(shí)空變化特征，以期為安徽省的生態(tài)環(huán)境保護(hù)和社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展提供理論依據(jù)。

1 研究區(qū)概況

安徽省位于我國(guó)華東地區(qū)，地理坐標(biāo)為北緯29°41′～34°38′、東經(jīng)114°54′～119°37′，總面積為14.01 萬(wàn)km2［13］。安徽省地勢(shì)西高東低，南高北低，地貌以山地、丘陵和平原為主，其中山地約占總面積的29.52%，丘陵占29.01%，平原占24.82%。安徽省屬暖溫帶與亞熱帶過(guò)渡區(qū)，受季風(fēng)影響較為顯著，四季變化明顯，為植被的生長(zhǎng)提供了良好的氣候條件。

2 數(shù)據(jù)來(lái)源與研究方法

2.1 數(shù)據(jù)來(lái)源與預(yù)處理

研究采用的遙感數(shù)據(jù)主要為MODIS-NDVI數(shù)據(jù)，來(lái)源于美國(guó)航空航天局提供的MOD13Q1 產(chǎn)品。該數(shù)據(jù)空間分辨率為250 m，時(shí)間分辨率為16 d，時(shí)間范圍為2000—2021 年。筆者利用MRT 工具將MODISNDVI 數(shù)據(jù)進(jìn)行拼接及空間投影轉(zhuǎn)換，投影方式為通用橫軸墨卡托投影（Universal Transverse Mercator Projection，UTM），坐標(biāo)系為GCS-WGS-84 坐標(biāo)系，并利用安徽省矢量文件對(duì)其進(jìn)行裁剪，得到最終的研究區(qū)域NDVI 數(shù)據(jù)。采用最大值合成法（Maximum Value Composition，MVC）將各旬?dāng)?shù)據(jù)合成逐年的月數(shù)據(jù)，然后利用均值法計(jì)算得到逐年的月平均數(shù)據(jù)和年平均數(shù)據(jù)。

研究采用的數(shù)字高程模型（Digital elevation model，DEM）數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心（https：//www.resdc.cn/），其空間分辨率為90 m，利用Arcgis軟件對(duì)其重采樣分辨率至250 m。

2.2 研究方法

2.2.1 Theil-Sen趨勢(shì)分析

采用Theil-Sen 趨勢(shì)分析對(duì)2000—2021 年安徽省NDVI的時(shí)空變化趨勢(shì)進(jìn)行分析。Theil-Sen趨勢(shì)分析是一種非參數(shù)統(tǒng)計(jì)的趨勢(shì)分析方法，其優(yōu)點(diǎn)在于樣本數(shù)據(jù)不受少數(shù)極端值的影響，Theil-Sen 趨勢(shì)分析計(jì)算公式為

式（1）中：Mi和Mj分別表示第i年和j年的NDVI值。如果β＞0，表示NDVI 呈增加趨勢(shì)，反之則表示呈下降趨勢(shì)。

2.2.2 Mann-Kendall統(tǒng)計(jì)檢驗(yàn)

采用Mann-Kendall 統(tǒng)計(jì)檢驗(yàn)方法對(duì)NDVI 變化趨勢(shì)顯著性進(jìn)行判斷。該方法作為一種非參數(shù)檢驗(yàn)方法，不要求所使用的樣本數(shù)據(jù)服從特定的分布［14］，而且不易受到異常值干擾，計(jì)算也較為簡(jiǎn)便，計(jì)算公式為

式（2）至（5）中：S表示統(tǒng)計(jì)量，V（S）表示S的方差，Z表示標(biāo)準(zhǔn)化后的統(tǒng)計(jì)量，n表示時(shí)間序列的長(zhǎng)度，sgn表示符號(hào)函數(shù)。在給定α=0.05置信水平判定NDVI變化趨勢(shì)的顯著性，|Z|≤Z1-α/2，說(shuō)明數(shù)據(jù)變化趨勢(shì)不顯著；|Z|＞Z1-α/2，說(shuō)明數(shù)據(jù)變化趨勢(shì)顯著［15］。

3 結(jié)果與分析

3.1 2000—2021年安徽省NDVI時(shí)間變化特征

3.1.1 月均NDVI變化特征

由圖1 可知，2000—2021 年安徽省月均NDVI 整體呈現(xiàn)先上升后下降的趨勢(shì)。具體來(lái)說(shuō)，1—5 月安徽省NDVI 呈上升趨勢(shì)，6 月安徽省NDVI 略微下降，7—8 月安徽省NDVI 再次呈上升趨勢(shì)，在8 月達(dá)最高值（0.77），9—12 月安徽省NDVI 呈下降趨勢(shì)。安徽省作為我國(guó)主要的農(nóng)業(yè)生產(chǎn)區(qū)，省內(nèi)大部分地區(qū)種植的糧食產(chǎn)物可以達(dá)到一年兩熟，且交替種植。在6 月安徽省NDVI出現(xiàn)輕微的下降，是由于6月安徽省各地種植的玉米、大豆、水稻等農(nóng)作物被收割，導(dǎo)致NDVI 有所下降。7—8 月，安徽省氣候條件較好，高溫多雨，有利于植被的生長(zhǎng)，所以該段時(shí)期植被覆蓋度較高。9—10月為農(nóng)作物的收獲期，安徽省NDVI 急劇下降，11—12月植被逐漸枯萎，安徽省NDVI 呈逐漸下降趨勢(shì)。整體而言，2000—2021 年安徽省月均NDVI 變化趨勢(shì)表現(xiàn)為雙峰特征，這與鄧凱等［1］的研究結(jié)果一致。

圖1 2000—2021年安徽省月均NDVI變化趨勢(shì)

3.1.2 年均NDVI變化特征

由圖2 可知，2000—2021 年安徽省年均NDVI 整體呈顯著增加趨勢(shì)（P＜0.05），增率為0.003 4/a，在2001年安徽省NDVI最?。?.52），在2021年安徽省NDVI達(dá)最大值（0.61），這表明2000—2021 年安徽省植被狀況整體呈持續(xù)改善的趨勢(shì)。

圖2 2000—2021年安徽省年均NDVI變化趨勢(shì)

3.2 2000—2021年安徽省NDVI空間變化特征

3.2.1 2000—2021年安徽省NDVI空間分布特征

2000—2021年安徽省NDVI空間分布存在明顯的空間異質(zhì)性，總體上呈現(xiàn)南高北低、西高東低的趨勢(shì)。安徽省NDVI高值區(qū)域主要分布在皖南山區(qū)和皖西大別山地區(qū)，是因?yàn)檫@些地區(qū)地勢(shì)較高，人類(lèi)活動(dòng)較少，且氣候條件較好，較為適宜植被的生長(zhǎng)，植被覆蓋率較高；NDVI 低值區(qū)域集中分布在淮河流域、長(zhǎng)江沿線和巢湖等流域；NDVI 中值區(qū)域集中分布在淮北平原、江淮丘陵地區(qū)，是因?yàn)檫@些地區(qū)地勢(shì)較為平坦，人類(lèi)活動(dòng)比較頻繁，大部分地區(qū)為農(nóng)作物區(qū)，植被覆蓋率較低。

3.2.2 2000—2021年安徽省NDVI空間變化趨勢(shì)

利用Theil-Sen 趨勢(shì)分析方法，分析了2000—2021 年安徽省NDVI 空間變化趨勢(shì)。2000—2021 年，安徽省大部分地區(qū)NDVI呈增加趨勢(shì)，其中NDVI呈增加趨勢(shì)的區(qū)域約占研究區(qū)總面積的89.33%，而NDVI呈降低趨勢(shì)的區(qū)域僅占研究區(qū)總面積的10.67%。Mann-Kendall 顯著性檢驗(yàn)分析結(jié)果表明：安徽省3.99%的區(qū)域NDVI呈顯著降低趨勢(shì)，主要集中在沿江平原區(qū)域長(zhǎng)江兩側(cè)地帶、皖中巢湖地區(qū)及淮河以北部分地區(qū)；6.68%的區(qū)域NDVI呈降低趨勢(shì)，但不顯著，主要集中在沿江平原和江淮丘陵地區(qū)；16.13%的區(qū)域NDVI 呈增加趨勢(shì)，但不顯著，主要集中在淮北平原和淮河以南地區(qū)；NDVI 呈顯著增加趨勢(shì)的區(qū)域占研究區(qū)總面積的73.20%，主要集中在淮河平原北部和東部、皖西大別山地區(qū)及皖南山區(qū)。

4 結(jié)論

筆者基于MODIS-NDVI 數(shù)據(jù)分析了2000—2021年安徽省NDVI時(shí)空變化特征，得出如下結(jié)論。

第一，在時(shí)間上，2000—2021 年安徽省月均NDVI總體呈先上升后下降的趨勢(shì)，在8 月NDVI 達(dá)最高值。2000—2021年安徽省年均NDVI整體呈顯著增加趨勢(shì)，增率為0.003 4/a，2021 年NDVI 達(dá)最高值（0.61），表明2000—2021年安徽省的植被生長(zhǎng)狀況得到了極大的改善。

第二，在空間上，2000—2021年安徽省NDVI高值區(qū)域主要分布在皖南山區(qū)和皖西大別山地區(qū)，NDVI低值區(qū)域主要分布在淮北平原和江淮丘陵地區(qū)。2000—2021 年安徽省大部分區(qū)域NDVI 呈增加趨勢(shì)，其中NDVI 顯著增加區(qū)域占比高達(dá)73.20%，顯著降低區(qū)域僅占3.99%。