王彤新 李 莎 左 嬋 郭曉華 張庭康 張秀娟

1 長(zhǎng)江大學(xué)園藝園林學(xué)院 湖北荊州 434000

2 中南大學(xué)地球科學(xué)與信息物理學(xué)院 長(zhǎng)沙 410083

3 中國(guó)科學(xué)院地理科學(xué)與資源研究所生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測(cè)與模擬重點(diǎn)實(shí)驗(yàn)室 北京 100101

植被作為陸地生態(tài)系統(tǒng)中連接土壤、 大氣、水文等要素的樞紐, 在水土保持、 氣候調(diào)節(jié)、 碳氮循環(huán)等生態(tài)系統(tǒng)功能中發(fā)揮著不可或缺的作用[1－2]。 植被覆蓋動(dòng)態(tài)變化可用于反映植被生長(zhǎng)狀況和健康程度[3], 是評(píng)價(jià)土壤地表植被狀況的一項(xiàng)關(guān)鍵參數(shù)。 近年來(lái)衛(wèi)星遙感技術(shù)因不受自然和人為因素約束并且能夠在短期內(nèi)覆蓋大尺度植被監(jiān)測(cè)數(shù)據(jù)[4], 被廣泛應(yīng)用于監(jiān)測(cè)全球植被動(dòng)態(tài)變化。 歸一化植被指數(shù)(NDVI) 具有消除儀器定標(biāo)、 輻射誤差的特點(diǎn)以及時(shí)間序列長(zhǎng)、 運(yùn)算簡(jiǎn)便的特性, 是當(dāng)今植被指數(shù)評(píng)價(jià)中描述植被覆蓋和生物量情況的重要參數(shù), 被普遍應(yīng)用于植被動(dòng)態(tài)監(jiān)測(cè)和植被年際變化研究[5]。 目前已有研究大多是關(guān)于高原和流域地區(qū)植被NDVI 時(shí)空演變與氣候變化響應(yīng)等方面, 如黃土高原、 青藏高原、 喀斯特地區(qū)、 黃河流域、 長(zhǎng)江流域等植被覆蓋演變特征[6－13]。 面對(duì)城市植被覆蓋動(dòng)態(tài)變化與大氣污染間相關(guān)研究尚比較欠缺, 系統(tǒng)地探討植被覆蓋時(shí)空演變特征和空氣質(zhì)量的響應(yīng)對(duì)改善區(qū)域環(huán)境質(zhì)量、 科學(xué)管理植被配置和治理生態(tài)環(huán)境具有重要科學(xué)價(jià)值。

荊州市作為我國(guó)首批國(guó)家歷史文化名城之一,擁有內(nèi)陸最廣、 密度最高的水域和得天獨(dú)厚的生物資源優(yōu)勢(shì), 是國(guó)家森林城市和維護(hù)長(zhǎng)江中游生態(tài)系統(tǒng)的重要屏障。 在經(jīng)濟(jì)快速發(fā)展和城市化進(jìn)程中, 城市空氣污染已成為全球面臨的重大問(wèn)題。城市綠地作為城市的心肺, 在自然環(huán)境修復(fù)與城市人居環(huán)境改善中發(fā)揮重要作用, 自開(kāi)展長(zhǎng)江防護(hù)林建設(shè)、 山水林田湖草沙生態(tài)保護(hù)與修復(fù)、 荊江水生態(tài)大保護(hù)、 森林生態(tài)系統(tǒng)恢復(fù)等一系列生態(tài)修復(fù)工程以來(lái)[14－15], 植被覆蓋與區(qū)域環(huán)境發(fā)生巨大變化, 但是目前還未有系統(tǒng)的評(píng)估, 因此,本文以湖北省荊州市為研究對(duì)象, 利用MOD13A3 NDVI 數(shù)據(jù)源分析2000—2019 年荊州市植被NDVI動(dòng)態(tài)變化, 為推進(jìn)城市綠色生態(tài)發(fā)展和實(shí)現(xiàn)生態(tài)格局可持續(xù)發(fā)展提供參考。

1 研究區(qū)概況

荊州市 (111° 15′—114° 05′ E, 29° 26′—31°37′N(xiāo))地處湖北省中南部, 江漢平原腹地, 是長(zhǎng)江中游交通樞紐之一, 總面積14.1 萬(wàn)km2, 涵蓋2 個(gè)市轄區(qū)(荊州區(qū)、 沙市區(qū)), 4 個(gè)縣級(jí)市(石首市、 松滋市、 洪湖市、 監(jiān)利市), 2 個(gè)縣(公安縣、 江陵縣) (圖1)。 全市地勢(shì)略呈西高東低, 由低山丘陵向崗地、 平原逐漸過(guò)渡, 形成以平原崗地為主, 兼有少量丘陵、 低山的地形地貌。研究區(qū)屬亞熱帶季風(fēng)濕潤(rùn)氣候, 四季分明、 熱量豐富、 光照適宜、 雨水充沛、 雨熱同季、 無(wú)霜期長(zhǎng)。

圖1 荊州市行政區(qū)域

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

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

NDVI 數(shù)據(jù)來(lái)自美國(guó)國(guó)家航空航天局(National Aeronautics and Space Administration,NASA) 提供的2000—2019 年MOD13A3 數(shù)據(jù)產(chǎn)品(https:/ /neo.scigsfc. nasa. gov), 空間分辨率為1 km, 時(shí)間分辨率為16 天, 數(shù)據(jù)格式為HDF；DEM 數(shù)字高程數(shù)據(jù)來(lái)自地理空間數(shù)據(jù)云(https:/ /www.gscloud.cn/), 空間分辨率為90 m；荊州市行政邊界地圖來(lái)自國(guó)家基礎(chǔ)地理信息系統(tǒng)1 ∶400 萬(wàn)數(shù)據(jù)。 荊州市2015—2019 年空氣質(zhì)量數(shù)據(jù)來(lái)自空氣質(zhì)量在線(xiàn)監(jiān)測(cè)分析平臺(tái) (https:/ /www.aqistudy.cn/)。

首先使用NASA 提供的MODIS 重投影工具(MRT) 軟件對(duì)MOD13A3 數(shù)據(jù)進(jìn)行拼接和投影,將HDF 格式轉(zhuǎn)為T(mén)IF 格式, 并在ArcGIS 10.2 軟件中裁剪出研究區(qū)范圍, 采用TIMESAT3.2 軟件的S－G 濾波方法[16], 對(duì)數(shù)據(jù)進(jìn)行平滑和降噪處理； 然后用批量最大值合成法合成月尺度的NDVI數(shù)據(jù), 再對(duì)NDVI 數(shù)據(jù)進(jìn)行反演減少或消除云、大氣和太陽(yáng)高度角對(duì)遙感的影響, 處理后的240幅NDVI 分布影像采用均值法生成年NDVI 序列數(shù)據(jù)； 最后根據(jù)國(guó)家空氣質(zhì)量監(jiān)測(cè)平臺(tái)提供的2015—2019 年空氣質(zhì)量指數(shù)月統(tǒng)計(jì)歷史數(shù)據(jù)求均值獲得年均數(shù)據(jù)。

2.2 研究方法

2.2.1 最大值合成法

最大值合成法(MVC) 通常是指從幾個(gè)圖像中取最大像元值重新生成一副新的圖像, 可以有效減少或消除云層、 天氣等因素對(duì)結(jié)果的影響[17]。 本研究采用MVC 法得到2000—2019 年植被年尺度NDVI 年數(shù)據(jù), 代表荊州植被年內(nèi)最佳生長(zhǎng)狀態(tài)。

2.2.2 NDVI 計(jì)算

通過(guò)MODIS 衛(wèi)星遙感獲取反映植被生長(zhǎng)狀態(tài)和空間分布密度的NDVI。 計(jì)算方法為近紅外波段反射值與紅光波段反射值之差與上兩者之和的比。公式如下[18]:

式(1) 中:Band2 為近紅外波段反射值；Band1 為紅光波段反射值；NDVI可以反映植被覆蓋度和植被生長(zhǎng)狀態(tài), 值域?yàn)椋?~1。NDVI<0 表示地面覆蓋為水系、 冰雪或云層等, 對(duì)可見(jiàn)光高反射；NDVI>0 表示地表有植被生長(zhǎng), 其值越高越趨于1, 植被覆蓋等級(jí)越好。

2.2.3 趨勢(shì)分析

采用一元線(xiàn)性回歸法分析變量間的線(xiàn)性關(guān)系及變化[19]。 本研究采用2000—2019 年荊州市NDVI 的年際變化表征荊州市植被在時(shí)間上的變化趨勢(shì), 借助F 檢驗(yàn)完成回歸方程的顯著性檢驗(yàn)。植被變化趨勢(shì)用Slope表示, 公式[20]為:

式(2) 中:n為年限(n＝20),i＝1, 2, 3,……,n；NDVIi代表第i年的數(shù)值；Slope代表植被NDVI 年際變化趨勢(shì)。 若Slope>0, NDVI 呈改善趨勢(shì)； 若Slope＝0, NDVI 在研究時(shí)段無(wú)增減；若Slope<0, NDVI 呈退化趨勢(shì)。

2.2.4 相關(guān)性分析

皮爾遜相關(guān)系數(shù)可以有效衡量?jī)蓚€(gè)變量間的相關(guān)程度。 為研究NDVI 與PM2.5、 NDVI 與PM10、NDVI 與CO、 NDVI 與SO2、 NDVI 與NO2、 NDVI與O3的相關(guān)性, 計(jì)算皮爾遜相關(guān)系數(shù)公式如下[21]:

式(3) 中:R為變量變化趨勢(shì)的相關(guān)性；xi代表第i年NDVI 值；yi代表對(duì)應(yīng)第i年P(guān)M2.5、PM10、 CO、 SO2、 NO2、 O3的年均值； 和分別表示要素樣本平均值；n為年限(n＝5)。R值域在－1~1, 若R>0 表示正相關(guān),R<0 表示負(fù)相關(guān),｜R｜值越大則相關(guān)性程度越高。

3 結(jié)果與分析

3.1 NDVI 時(shí)空分布特征

3.1.1 NDVI 時(shí)間變化特征

通過(guò)對(duì)荊州市植被NDVI 年際變化的分析可知(圖2), 2000—2019 年荊州市NDVI 年均變化呈現(xiàn)波動(dòng)上升, 從0.431 (2000 年) 增長(zhǎng)到0.474 (2019 年), 增加了0.043, 每年變化率為0.215%, 說(shuō)明研究時(shí)段荊州市植被NDVI 整體變化趨勢(shì)以改善為主。 2000—2019 年荊州市NDVI均值為0.454, 快速波動(dòng)變化集中在2000—2013年, 其中3 大上升階段的起點(diǎn)分別為2002、 2009和2012 年； 2013 年后NDVI 穩(wěn)定于0.470 且年際變化速率放緩。 各區(qū)域NDVI 均值(表1) 表明,與2000—2009 年相比, 2010—2019 年植被NDVI均增加, 增加范圍在4.1%~8.2%, 增量最高的區(qū)域?yàn)榻昕h。

表1 荊州市各區(qū)域NDVI 均值

圖2 2000—2019 年荊州市植被NDVI 變化趨勢(shì)

3.1.2 NDVI 空間變化特征

在ArcGIS 軟件中采用自然間斷點(diǎn)分級(jí)法把荊州市2000—2019 年NDVI 從低值到高值依次聚類(lèi)為6 個(gè)等級(jí), 劃分標(biāo)準(zhǔn)為≤0.3、 0.3~0.4、0.4~0.5、 0.5~0.6、 0.6~0.7、 ≥0.7, 分別對(duì)應(yīng)裸地、 低植被覆蓋度、 較低植被覆蓋度、 中植被覆蓋度、 較高植被覆蓋度、 高植被覆蓋度。如圖3 所示, 荊州市整體植被覆蓋狀況因植被類(lèi)型和土地利用類(lèi)型不同在空間上呈現(xiàn)明顯地域差異, 大部分地區(qū)處于中、 較低植被覆蓋區(qū),東部為較低植被覆蓋區(qū), 西部和中部江陵縣為較高植被覆蓋區(qū)。 可見(jiàn), 荊州市NDVI 呈現(xiàn)“西高東低” 的空間分布特征, 較高植被覆蓋區(qū)主要分布在松滋市西部和江陵縣, 較低植被覆蓋區(qū)主要分布在洪湖市。 2010—2019 年相較于2000—2009 年, 荊州區(qū)、 松滋市、 石首市、 江陵縣和公安縣的較高植被覆蓋區(qū)增加, 監(jiān)利縣和洪湖市中植被覆蓋區(qū)增加, 其中江陵縣整體呈現(xiàn)較高植被覆蓋度。

圖3 2000—2019 年荊州市植被NDVI 空間分布

3.2 NDVI 動(dòng)態(tài)變化特征

2000—2019 年荊州市植被NDVI 呈現(xiàn)改善趨勢(shì)(圖4), 明顯改善的區(qū)域?yàn)樗勺淌形鞑康母吆０蔚貐^(qū)。 其中, 2000—2009 年荊州市植被NDVI呈現(xiàn)改善趨勢(shì), 大部分區(qū)域在0.5%~1%范圍變化, 主要位于荊州區(qū)、 公安縣和監(jiān)利市。 2010—2019 年荊州市植被NDVI 除在松滋市西部高海拔

圖4 2000—2019 年荊州市植被NDVI 變化趨勢(shì)及顯著性檢驗(yàn)

地區(qū)和其他零星區(qū)域呈現(xiàn)改善趨勢(shì)之外, 其余呈現(xiàn)退化趨勢(shì), 可能與全球變化下城市化進(jìn)程加快和人口密度增長(zhǎng)帶來(lái)的生態(tài)環(huán)境破環(huán)相關(guān)。 總體而言, 2000—2019 年荊州NDVI 變化趨勢(shì)呈現(xiàn)上升, 但增長(zhǎng)浮動(dòng)低于前10 年, 而后10 年可能受城鎮(zhèn)化的影響, 使荊州整體NDVI 出現(xiàn)大面積退化。

3.3 NDVI 對(duì)空氣質(zhì)量的響應(yīng)

3.3.1 空氣質(zhì)量年際變化特征

荊州市空氣污染物的年際變化濃度(圖5)表明, 2015—2019 年SO2和NO2濃度符合國(guó)家二級(jí)標(biāo)準(zhǔn)?環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)? (GB3095－2012),而可吸入顆粒物PM2.5和PM10未達(dá)到國(guó)家二級(jí)標(biāo)準(zhǔn), 說(shuō)明荊州市空氣質(zhì)量主要受到PM2.5和PM10的影響。

圖5 荊州市2015—2019 年P(guān)M2.5、 SO2、 PM10、 NO2年濃度年際變化

2015—2019 年荊州市空氣污染變化趨勢(shì)(表2) 表明, PM2.5變化率為－6.325 μg?m－3?a－1(R2＝0.916), PM10變化率為－7.350 μg?m－3?a－1(R2＝0.833), SO2變化率為－4.308 μg?m－3?a－1(R2＝0.997), CO 變化率為－0.059 μg?m－3?a－1(R2＝0.703), 說(shuō)明在湖北省大氣污染防治攻堅(jiān)幫扶和長(zhǎng)江經(jīng)濟(jì)帶生態(tài)環(huán)境保護(hù)行動(dòng)下, 荊州市空氣質(zhì)量正在逐步好轉(zhuǎn)。

表2 2015—2019 年空氣污染物變化趨勢(shì)

3.3.2 NDVI 與空氣質(zhì)量相關(guān)分析

采用2015—2019 年NDVI 年際變化與同時(shí)段PM2.5、 PM10、 CO、 NO2、 SO2、 O3的年度變化趨勢(shì)開(kāi)展相關(guān)分析, 結(jié)果(表3) 表明: NDVI 與PM2.5、 CO 呈顯著負(fù)相關(guān) (P<0.05), 5 年內(nèi)PM2.5濃度降低23.92 μg?m－3, CO 濃度降低0.28μg?m－3； NDVI 與PM10、 NO2呈極顯著負(fù)相關(guān)(P<0.01), 5 年內(nèi)PM10濃度降低26.17 μg?m－3,NO2濃度降低4.42 μg?m－3； 而NDVI 與SO2、 O3相關(guān)性不顯著, 5 年內(nèi) SO2濃度降低17.2 μg?m－3, O3濃度降低4.42 μg?m－3。 PM10與NDVI 負(fù)相關(guān)顯著性最高,R為－0.992, 其次分別是 NO2( － 0.963)、 PM2.5( － 0.953)、CO (－0.880)。 由此可見(jiàn), 植被NDVI 的增加可以有效地改善空氣質(zhì)量, 降低空氣中PM2.5、 PM10、CO 和NO2的含量, 特別是降低PM10的含量。

表3 2015—2019 年NDVI 與空氣污染物變化的相關(guān)性

4 討論

4.1 近20 年NDVI 變化原因

本研究結(jié)果顯示, 2000—2019 年荊州市植被NDVI 以每年變化0.215%的速率增加, 2000—2009 年的NDVI 年際增長(zhǎng)變化趨勢(shì)高于2010—2019 年。 空間范圍內(nèi)荊州市大多為中植被覆蓋區(qū)域, 較高植被覆蓋區(qū)分布在松滋市、 江陵縣, 較低植被覆蓋區(qū)分布在洪湖市。 近20 年荊州市植被的改善得益于地方部門(mén)開(kāi)展的一系列生態(tài)保護(hù)措施, 包括長(zhǎng)江防護(hù)林建設(shè)、 碼頭治理、 退垸還湖、山水林田湖草沙生態(tài)修復(fù)和長(zhǎng)江大保護(hù)工程等[22－24], 這些工程的實(shí)施促進(jìn)荊州植被NDVI 的增加。 湖北省生態(tài)文明建設(shè)水平持續(xù)改善, 環(huán)境—經(jīng)濟(jì)—社會(huì)系統(tǒng)的協(xié)調(diào)發(fā)展水平日益提高, 荊州市在2009 年開(kāi)展國(guó)家園林城市活動(dòng)后, 2012年授予國(guó)家園林城市稱(chēng)號(hào), 這可以解釋植被NDVI在2009 年后呈上升趨勢(shì), 并且在2012 后植被變化趨勢(shì)較為穩(wěn)定。 荊州市作為長(zhǎng)江大保護(hù)中游的生態(tài)屏障, 以創(chuàng)建生態(tài)文明先行示范區(qū)為目標(biāo),在監(jiān)利市開(kāi)展長(zhǎng)江防護(hù)林建設(shè)、 江陵縣開(kāi)展退田還湖工程、 公安縣打造植樹(shù)10.67 萬(wàn)hm2(16 萬(wàn)畝), 現(xiàn)如今荊州入圍湖北長(zhǎng)江荊江段及洪湖流域山水林田湖草沙一體化保護(hù)和修復(fù)工程項(xiàng)目,綠色發(fā)展的步伐從未停止。

4.2 NDVI 變化與空氣質(zhì)量的關(guān)系

荊州市NDVI 的增加對(duì)改善空氣質(zhì)量具有重要意義。 荊州市2015—2019 年大氣污染物均呈現(xiàn)下降趨勢(shì), 但PM2.5和PM10濃度年均值未達(dá)國(guó)家二級(jí)標(biāo)準(zhǔn), 是荊州市主要的大氣污染源, 同時(shí)也是湖北省主要大氣污染物[25]。 相關(guān)性分析結(jié)果表明NDVI 與大氣污染物PM2.5、 PM10、 CO 和NO2均呈現(xiàn)顯著負(fù)相關(guān), 負(fù)相關(guān)程度最大的因素為PM10, 其次分別為NO2、 PM2.5、 CO。 本研究指出降低PM2.5和PM10是改善荊州空氣質(zhì)量的關(guān)鍵,綠地能有效吸附空氣污染物, 植被NDVI 提升能有效改善空氣質(zhì)量, 與已有研究結(jié)果一致[26－27]。實(shí)施相關(guān)生態(tài)修復(fù)的長(zhǎng)江防護(hù)林建設(shè)、 森林生態(tài)系統(tǒng)恢復(fù)等工程, 一方面可以提升荊州市植被覆蓋和空氣環(huán)境質(zhì)量； 另一方面能牢筑長(zhǎng)江中上游生態(tài)安全屏障, 確保一江清水永續(xù)東流、 構(gòu)建和諧共生生態(tài)格局[28]。

4.3 研究局限

本文研究荊州市植被NDVI 時(shí)空演變特征與其對(duì)空氣質(zhì)量的影響, 對(duì)優(yōu)化城市綠地系統(tǒng)規(guī)劃、城市綠地建設(shè)和實(shí)現(xiàn)區(qū)域長(zhǎng)江岸線(xiàn)生態(tài)環(huán)境保護(hù)具有重要意義。 研究中未考慮到季相或者極端天氣等對(duì)植被NDVI 變動(dòng)的影響, 未來(lái)應(yīng)充分探索植被覆蓋動(dòng)態(tài)變化與多環(huán)境要素的關(guān)系； 空氣質(zhì)量除了受植被NDVI 時(shí)空演變特征影響外, 還受到工業(yè)、 交通、 人為活動(dòng)和氣象因素影響, 因此本研究可能會(huì)夸大植被提升改善城市空氣質(zhì)量的效果； 本文重點(diǎn)研究了荊州市域和空氣質(zhì)量的關(guān)系, 暫未考慮各縣市綠地空間分布對(duì)空氣質(zhì)量的影響, 今后對(duì)此應(yīng)進(jìn)一步深入研究。

5 結(jié)論

1) 2000—2019 年荊州市植被NDVI 呈增加趨勢(shì), 增速每年為0.215%, 植被覆蓋程度較好。 在空間分布上呈現(xiàn)“西高東低” 的分異特征, 江陵縣植被覆蓋情況最好, 松滋市次之, 荊州區(qū)植被增長(zhǎng)比率最高。

2) 荊州市2000—2019 年NDVI 整體動(dòng)態(tài)變化呈現(xiàn)改善趨勢(shì), 植被覆蓋增加區(qū)域大于減少區(qū)域,但增長(zhǎng)幅度相較于2000—2009 年植被覆蓋增長(zhǎng)趨勢(shì)放緩。

3) 荊州市主要空氣污染物為PM2.5和PM10。NDVI 與大氣污染物PM2.5、 PM10、 CO、 NO2呈顯著負(fù)相關(guān), NDVI 對(duì)PM10的影響最為顯著, 表明植被覆蓋度增加有助于降低可吸入顆粒物濃度。