山東省城鎮(zhèn)化與植被指數(shù)的時(shí)空關(guān)系研究

曹紀(jì)元 高劍飛 董超

摘? 要：為研究近20年城鎮(zhèn)化背景下，山東省的植被EVI時(shí)空變化受城鎮(zhèn)化進(jìn)程的影響，基于2001—2020年的MODIS-EVI數(shù)據(jù)檢驗(yàn)探討山東省植被的時(shí)空分布規(guī)律和植被變化趨勢(shì)，分析了植被變化與夜光數(shù)據(jù)的相關(guān)性。結(jié)果表明：山東省近20年整體植被呈增加趨勢(shì)，總增長(zhǎng)趨近年均值的11.80%。由城鎮(zhèn)擴(kuò)張引起的夜光指數(shù)變化對(duì)城市周邊植被覆蓋具有抑制作用，而遠(yuǎn)離城市則表現(xiàn)為一定的促進(jìn)作用。

關(guān)鍵詞：EVI;植被趨勢(shì);夜光數(shù)據(jù);城鎮(zhèn)化

中圖分類(lèi)號(hào)：Q948;X173? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：2096-4706（2020）24-0030-05

Spatiotemporal Relationship Study between Urbanization and Vegetation Index in Shandong Province

CAO Jiyuan，GAO Jianfei，DONG Chao

（College of Information Science and Engineering，Shandong Agricultural University，Taian? 271018，China）

Abstract：In order to study the situation that vegetation EVI spatiotemporal change is influenced by urbanization process in Shandong Province under the recent 20 years urbanization background，based on MODISEVI data validation from 2001 to 2020，the spatiotemporal distribution law and change trend of vegetation in Shandong Province are discussed，and the correlation between vegetation change and luminous data is analyzed. The results show that the overall vegetation in Shandong Province in recent 20 years present a rising trend，and the total growth gets close to 11.80%? of the annual average. The variation of luminous index caused by urban expansion will inhibit the vegetation coverage around the city，while far away from the city has a certain promotion to it.

Keywords：EVI;vegetation trend;luminous data;urbanization

0? 引? 言

植被作為自然界的重要組成部分，不僅是陸地生態(tài)系統(tǒng)中重要的一環(huán)，也是反映人與自然的一面鏡子[1]。近些年來(lái)，隨著我國(guó)城鎮(zhèn)化和工業(yè)化的迅速發(fā)展，人民對(duì)于擁有良好生態(tài)環(huán)境的美好生活也越來(lái)越向往[2]。而土地利用/覆被變化（LUCC）是引起全球環(huán)境變化的重要組成部分和主要原因之一，由于其對(duì)于引導(dǎo)我國(guó)土地資源合理利用具有非常重要的意義，我國(guó)學(xué)者開(kāi)展了一系列相關(guān)研究，以深化土地，環(huán)境，氣候，人類(lèi)活動(dòng)等因子對(duì)其產(chǎn)生的影響[3]。賈路等研究發(fā)現(xiàn)西安市降雨量與植被覆蓋度均值呈現(xiàn)出增加趨勢(shì)[4]。趙鴻雁等研究發(fā)現(xiàn)生長(zhǎng)季NDVI對(duì)不同氣候因子的響應(yīng)程度不同，氣溫最強(qiáng)，降水次之，太陽(yáng)輻射最弱[5]。劉啟興等通過(guò)分析氣候因素及人類(lèi)活動(dòng)對(duì)植被覆蓋的影響，發(fā)現(xiàn)氣候因素對(duì)黃河源區(qū)植被生長(zhǎng)起主要促進(jìn)作用，而人類(lèi)活動(dòng)則對(duì)黃河源區(qū)的生態(tài)環(huán)境產(chǎn)生一定的積極意義[6]。通過(guò)梳理相關(guān)文獻(xiàn)可知，局部地區(qū)的植被生長(zhǎng)對(duì)氣候因子的響應(yīng)程度存在顯著差異，人為因素對(duì)植被趨勢(shì)變化的影響也有利有弊。

2019年2月份，美國(guó)航天局公布的衛(wèi)星數(shù)據(jù)表明，在2000年至2017年間，全球綠化面積增加了5.00%，其中約四分之一來(lái)自中國(guó)，貢獻(xiàn)比例居全球首位。分析顯示，中國(guó)的貢獻(xiàn)中42.00%來(lái)自植樹(shù)造林，32.00%來(lái)自集約農(nóng)業(yè)。而山東省作為我國(guó)的糧食大省，故對(duì)山東省城鎮(zhèn)化與植被指數(shù)的時(shí)空關(guān)系的研究具有重要意義。本文利用一元線(xiàn)性回歸、M-K趨勢(shì)檢驗(yàn)對(duì)研究區(qū)域進(jìn)行了EVI像元尺度的分析[7]，模擬山東省2001—2020年間的EVI變化趨勢(shì)，通過(guò)各個(gè)像元的變化趨勢(shì)反映整個(gè)山東省的變化規(guī)律，通過(guò)M-K趨勢(shì)檢驗(yàn)對(duì)2001年—2020年的MODIS數(shù)據(jù)進(jìn)行時(shí)序分析，得出山東省植被變化趨勢(shì)，以及利用NPP/VIRS夜光數(shù)據(jù)對(duì)城市區(qū)域進(jìn)行相關(guān)性分析，探討山東省夜光數(shù)據(jù)與EVI的關(guān)系。

1? 研究區(qū)域概況

山東省位于我國(guó)華東地區(qū)，介于北緯34°22.9′～38° 24.01′，東經(jīng)114°47.5′～122°42.3′之間，面積達(dá)15.79萬(wàn)平方千米，占全國(guó)總面積的1.64%。山東省作為我國(guó)的產(chǎn)糧大省，種植面積占到總面積的53.82%。山東中部山地突起，西南、西北低洼平坦，東部緩丘起伏，地形以丘陵為主，溫帶季風(fēng)氣候，降水集中，雨熱同季，春秋短暫，冬夏較長(zhǎng)，年均日照時(shí)數(shù)2 290～2 890小時(shí)，年平均降水量一般在550～950毫米。

2? 數(shù)據(jù)獲取與處理

2.1? 土地覆蓋數(shù)據(jù)及EVI數(shù)據(jù)源

研究采用的2001—2020年間的MOD13A2產(chǎn)品，所有數(shù)據(jù)在USGS網(wǎng)站（https：//ladsweb.nascom.nasa.gov/search）上獲取。MOD13A2以每像素為基礎(chǔ)，以1公里（km）的空間分辨率提供植被指數(shù)（VI）值。有兩個(gè)主要植被層。第一個(gè)是歸一化植被指數(shù)（NDVI），它被稱(chēng)為現(xiàn)有的美國(guó)國(guó)家海洋與大氣管理局-先進(jìn)超高分辨率輻射計(jì)（NOAA-AVHRR）得出的NDVI的連續(xù)性指數(shù)。第二植被層是增強(qiáng)植被指數(shù)（EVI），在高生物量地區(qū)具有更高的敏感性。該產(chǎn)品的算法會(huì)從16天的所有采集中選擇最佳的可用像素值。使用的標(biāo)準(zhǔn)是低云，低視角和最高NDVI/EVI值。由于NDVI指數(shù)對(duì)于高植被區(qū)域具有較低的靈敏度，所以采用MOD13A2增強(qiáng)植被指數(shù)EVI對(duì)山東省進(jìn)行植被趨勢(shì)分析，時(shí)間分辨率為16天，空間分辨率為1 km[8]。

原始數(shù)據(jù)通過(guò)MRT（MODIS Reprojection Tools）工具統(tǒng)一進(jìn)行投影拼接轉(zhuǎn)換，轉(zhuǎn)換為UTM投影，坐標(biāo)系為WGS84坐標(biāo)，并利用山東省矢量邊界對(duì)處理后的數(shù)據(jù)進(jìn)行裁剪[6]。

2.2? 夜光數(shù)據(jù)源

研究還使用了2014年—2020年的Soumi-NPP/VIIRS觀測(cè)夜間燈光數(shù)據(jù)表征城市面積變化。VIRS是高分辨率輻射儀AVHRR和地球觀測(cè)系列中分辨率成像光譜儀MODIS系列的拓展和改進(jìn)。VIIRS數(shù)據(jù)可用來(lái)測(cè)量云量和氣溶膠特性、海洋水色、海洋和陸地表面溫度、海冰運(yùn)動(dòng)和溫度、火災(zāi)、和地球反照率。通過(guò)對(duì)年際數(shù)據(jù)中值求解，表征夜光數(shù)據(jù)年際間變化。

3? 時(shí)空趨勢(shì)相關(guān)分析方法

本研究的數(shù)據(jù)分析方法主要采用了一元線(xiàn)性回歸和M-K趨勢(shì)檢驗(yàn)、相關(guān)性分析。

通過(guò)一元線(xiàn)性回歸獲取EVI指數(shù)的變化趨勢(shì)，對(duì)EVI指數(shù)進(jìn)行M-K趨勢(shì)檢驗(yàn)，從統(tǒng)計(jì)學(xué)的角度揭曉植被EVI指數(shù)隨時(shí)間的變化趨勢(shì)，利用相關(guān)性分析方法對(duì)夜光指數(shù)年均值數(shù)據(jù)和植被EVI進(jìn)行相關(guān)性定量解算。

3.1? 線(xiàn)性回歸和M-K趨勢(shì)分析

本文采用線(xiàn)性回歸分析法定量研究EVI的變化趨勢(shì)，該方法可以模擬每個(gè)柵格的變化趨勢(shì)，以單個(gè)像元時(shí)間變化特征反映整個(gè)空間變化規(guī)律，綜合反映時(shí)空格局演變[9]。該方法通過(guò)對(duì)預(yù)處理后的MOD13A2 EVI數(shù)據(jù)集進(jìn)行逐像元應(yīng)用，獲取EVI數(shù)據(jù)的時(shí)序變化特征，計(jì)算公式如式（1）所示：

（1）

Slope為2001—2020年間時(shí)間分辨率為16 dpi的EVI數(shù)據(jù)的變化趨勢(shì);n為年份總數(shù)，i為年序數(shù)，EVIi為第i年的EVI值。當(dāng)Slope為正時(shí)，說(shuō)明植被狀況趨向正向趨勢(shì)發(fā)展;反之，當(dāng)Slope為負(fù)時(shí)，則說(shuō)明植被狀況趨于退化[10]。

M-K趨勢(shì)檢驗(yàn)是一種常用于綠色植被的非參數(shù)檢驗(yàn)方法。具有不受少量異常值干擾，定量化程度高，計(jì)算方便等優(yōu)點(diǎn)。通過(guò)對(duì)2001年—2020年的時(shí)間序列進(jìn)行M-K趨勢(shì)檢

驗(yàn)，判別EVI數(shù)據(jù)隨時(shí)間的變化趨勢(shì)。計(jì)算公式如式（2）～式（5）所示：

（2）

（3）

（4）

（5）

Zmk為單個(gè)像元的M-K標(biāo)準(zhǔn)化檢測(cè)統(tǒng)計(jì)量。S統(tǒng)計(jì)量是通過(guò)前后影像的EVI數(shù)據(jù)計(jì)算得到的，Var（S）為S統(tǒng)計(jì)量的方差，求解公式可簡(jiǎn)化為式（5）。本文將其置于0.05置信水平上進(jìn)行分析，即Zmk>1.96。

3.2? 相關(guān)性分析

通過(guò)對(duì)夜光數(shù)據(jù)進(jìn)行相關(guān)分析，從人文活動(dòng)方面進(jìn)行驅(qū)動(dòng)因子分析。采用皮爾斯相關(guān)系數(shù)方法，計(jì)算如式（6）所示：

（6）

ra為夜光數(shù)據(jù)的相關(guān)系數(shù)，ai，bi分別為第i年的夜光數(shù)值和EVI數(shù)據(jù)，，為n年間的夜光數(shù)值和EVI均值。

4? 結(jié)果與分析

4.1? 山東省植被覆蓋度分析

對(duì)2001年—2020年間的山東省EVI數(shù)據(jù)進(jìn)行均值操作，得到山東省的植被覆蓋情況，如圖1所示，可以發(fā)現(xiàn)，山東省EVI多年均值的高值主要分布于魯西北地區(qū)（EVI>0.3），而低值則主要分布在沿海及城市周邊區(qū)域（EVI<0.1）。圖2所示為山東省近20年間的EVI均值年際變化情況。結(jié)果顯示，山東省20年間植被覆蓋呈現(xiàn)出增長(zhǎng)趨勢(shì)——2001—2020的EVI總增長(zhǎng)趨近年均值的11.80%，20年間全省EVI增長(zhǎng)率達(dá)到23.30%。

通過(guò)對(duì)20年間山東省EVI影像進(jìn)行一元線(xiàn)性回歸分析，得出每個(gè)柵格的變化趨勢(shì)，如圖3所示。從空間分布上可以發(fā)現(xiàn)山東省城市區(qū)域及沿海區(qū)域的EVI呈現(xiàn)出明顯的下降趨勢(shì)，而EVI存在明顯增長(zhǎng)區(qū)域主要分布在魯西地區(qū)。總體結(jié)果表明，山東省近20年間植被指數(shù)顯著呈增強(qiáng)態(tài)勢(shì)，而圖中藍(lán)色部分明顯集中于城鎮(zhèn)周邊，表明城市周邊區(qū)域的植被覆蓋出現(xiàn)了退化現(xiàn)象。

對(duì)山東省2001年—2020年EVI指數(shù)做M-K趨勢(shì)分析，如圖4所示。山東省2001—2020年EVI趨勢(shì)分析表明：UF統(tǒng)計(jì)量在2001年—2003年間小于0，表明山東EVI指數(shù)整體呈減少趨勢(shì);在2004年—2020年間，UF統(tǒng)計(jì)量一直大于0，且于2011年UF與UB相交于一點(diǎn)，且相交點(diǎn)超過(guò)顯著性α小于0.05的臨界值1.96，可以由此判斷出自2007年后增長(zhǎng)趨勢(shì)明顯。

4.2? 夜光數(shù)據(jù)與植被指數(shù)分析

對(duì)2014年—2020年夜光數(shù)據(jù)進(jìn)行均值處理，將夜光數(shù)據(jù)與EVI進(jìn)行皮爾斯相關(guān)系數(shù)處理，分別如圖5，圖6所示。

統(tǒng)計(jì)結(jié)果如表1所示，表明EVI與夜光指數(shù)年均值呈負(fù)相關(guān)的區(qū)域面積占總面積的23.33%，而EVI與夜光指數(shù)年均值呈正相關(guān)的區(qū)域面積占總面積的76.67%，表明山東省植被生長(zhǎng)與夜光指數(shù)年均值關(guān)系較大。結(jié)合圖6，我們可以發(fā)現(xiàn)，研究區(qū)EVI與夜光指數(shù)年均值呈正相關(guān)的地區(qū)主要分布在魯中山地丘陵地區(qū)（R>0.4）;而對(duì)于植被的抑制作用則主要集中于城市周邊地區(qū)（R<-0.4）。因此，城鎮(zhèn)化對(duì)周邊植被覆蓋有一定的抑制作用。主要表現(xiàn)在城市擴(kuò)張過(guò)程中所引起的土地利用/覆蓋變化等。而在遠(yuǎn)離城市地區(qū)，植被EVI指數(shù)呈現(xiàn)出正增長(zhǎng)的態(tài)勢(shì)，其中原因仍需要進(jìn)一步的探討。故，城鎮(zhèn)化進(jìn)程對(duì)于植被覆蓋具有一定的促進(jìn)作用。

為進(jìn)一步探討人類(lèi)活動(dòng)對(duì)EVI的影響，在2014—2020年的夜光指數(shù)年均值的基礎(chǔ)上，對(duì)青島市進(jìn)行剖面分析，如圖7所示。分別對(duì)夜光指數(shù)與植被指數(shù)相關(guān)系數(shù)、青島市土地覆蓋類(lèi)型、夜光指數(shù)剖面圖進(jìn)行剖面分析，如圖8，9，10所示，對(duì)根據(jù)圖10可以發(fā)現(xiàn)，青島市城市內(nèi)夜光指數(shù)顯著高于城市外側(cè)的農(nóng)業(yè)用地，夜光指數(shù)在城市內(nèi)部達(dá)到峰值，遠(yuǎn)離城市則表現(xiàn)為低值。結(jié)合圖8圖9，可以發(fā)現(xiàn)在青島市土地覆蓋類(lèi)型發(fā)生改變時(shí)，往往也是夜光指數(shù)與植被指數(shù)相關(guān)系數(shù)發(fā)生較大波動(dòng)的時(shí)候。特別是在青島市城鄉(xiāng)交界處，夜光指數(shù)與植被指數(shù)呈負(fù)相關(guān)系數(shù)關(guān)系，而在城市內(nèi)部相關(guān)系數(shù)以正相關(guān)系數(shù)為主，這可能與青島市近些年來(lái)注重城市綠化有關(guān)，如2011年青島市政府通過(guò)《青島市城市綠化條例》，使城市植被覆蓋度有了較大的提升。因此，城鎮(zhèn)化對(duì)植被覆蓋具有抑制作用，主要集中于城鄉(xiāng)交接的土地利用/覆蓋變化區(qū)域;而對(duì)于遠(yuǎn)離城市區(qū)域則以促進(jìn)作用為主，至于原因仍需要進(jìn)行相關(guān)研究。

5? 結(jié)? 論

本文通過(guò)對(duì)山東省近20年來(lái)的增強(qiáng)型植被指數(shù)進(jìn)行一元線(xiàn)性回歸與M-K趨勢(shì)檢驗(yàn)，結(jié)合土地覆蓋圖，通過(guò)分析夜光遙感影像，對(duì)近20年來(lái)的植被情況進(jìn)行了初步的探討，并分析了人類(lèi)活動(dòng)對(duì)其所造成的影響。

山東省EVI指數(shù)呈現(xiàn)出西北高東部低的空間分布。近20年間山東省植被覆蓋度呈現(xiàn)出明顯增長(zhǎng)態(tài)勢(shì)，且從2007年起增長(zhǎng)趨勢(shì)顯著，2020年同2001年相比，EVI增長(zhǎng)了23.30%，近20年的EVI總增長(zhǎng)也趨近于年均值的11.80%。

城鎮(zhèn)化對(duì)植被覆蓋度起到了一定程度的促進(jìn)作用。促進(jìn)作用主要分布于遠(yuǎn)離城鎮(zhèn)區(qū)域，其中以魯中區(qū)域貢獻(xiàn)較大（R>0.4）;魯中區(qū)域由于多為山地丘陵地貌，受人類(lèi)活動(dòng)負(fù)面影響較弱;而大部分城市的邊緣區(qū)域植被呈現(xiàn)顯著的減少趨勢(shì)，且與夜光指數(shù)具有顯著負(fù)相關(guān)關(guān)系，表明城市的發(fā)展和擴(kuò)張對(duì)于周邊植被會(huì)產(chǎn)生不利影響，但本文只考慮了人類(lèi)活動(dòng)中的夜光指數(shù)，對(duì)人類(lèi)活動(dòng)中的開(kāi)墾、污染等人類(lèi)活動(dòng)因子未加考慮，對(duì)深刻揭示城鎮(zhèn)化與植被覆蓋的關(guān)系的研究還有待繼續(xù)。

參考文獻(xiàn)：

[1] 毛德華，王宗明，羅玲，等.基于MODIS和AVHRR數(shù)據(jù)源的東北地區(qū)植被NDVI變化及其與氣溫和降水間的相關(guān)分析 [J].遙感技術(shù)與應(yīng)用，2012，27（1）：77-85.

[2] 翟紹果，諶基東.共建美好生活的時(shí)代蘊(yùn)意、內(nèi)涵特質(zhì)與實(shí)現(xiàn)路徑 [J].西北大學(xué)學(xué)報(bào)（哲學(xué)社會(huì)科學(xué)版），2017，47（6）：20-26.

[3] 劉新衛(wèi)，陳百明，史學(xué)正.國(guó)內(nèi)LUCC研究進(jìn)展綜述 [J].土壤，2004（2）：132-135.

[4] 賈路，任宗萍，李占斌，等.2000-2013年西安市植被覆蓋度時(shí)空演變 [J].水土保持研究，2019，26（6）：274-279.

[5] 趙鴻雁，陳英，周翼，等，甘肅中東部植被生長(zhǎng)季NDVI時(shí)空變化及其對(duì)氣候因子的響應(yīng) [J].干旱區(qū)地理，2019，42（6）：1427-1435.

[6] 劉啟興，董國(guó)濤，景海濤，等，2000-2016年黃河源區(qū)植被NDVI變化趨勢(shì)及影響因素 [J].水土保持研究，2019，26（3）：86-92.

[7] 徐建華.現(xiàn)代地理學(xué)中的數(shù)學(xué)方法 [M].北京：高等教育出版社，2002.

[8] 賀忠華，張育慧，何月，等，浙江省近20年植被變化趨勢(shì)及驅(qū)動(dòng)因子分析 [J]，生態(tài)環(huán)境學(xué)報(bào)，2020，29（8），1530-1539.

[9] 趙維清，李經(jīng)緯，褚琳，等.近10年湖北省植被指數(shù)時(shí)空變化特征及其驅(qū)動(dòng)力 [J].生態(tài)學(xué)報(bào)，2019，39（20）：7722-7736.

[10] 何彬方，馮妍，吳文玉，等.安徽省近十年植被指數(shù)時(shí)空變化特征 [J].生態(tài)學(xué)雜志，2010，29（10）：1912-1918.

作者簡(jiǎn)介：曹紀(jì)元（2000—），男，漢族，浙江金華人，本科在讀，研究方向：遙感科學(xué)與技術(shù);高劍飛（2000—），男，漢族，山東臨沂人，本科在讀，研究方向：測(cè)繪程序與地理信息系統(tǒng);董超（1984—），男，漢，山東昌邑人，講師，博士，研究方向：土地資源利用。