卜亞勤，丁海勇

(1.南京信息工程大學(xué) 地理科學(xué)學(xué)院，南京 210044；2.南京信息工程大學(xué) 遙感與測(cè)繪工程學(xué)院，南京 210044)

0 引言

植被物候是指植被在生長(zhǎng)過(guò)程中受人類(lèi)活動(dòng)和環(huán)境因素影響而發(fā)生的發(fā)芽、展葉、開(kāi)花、結(jié)實(shí)、落葉等現(xiàn)象[1]，反映植被不斷適應(yīng)環(huán)境季節(jié)性變化的生長(zhǎng)發(fā)育規(guī)律。植被物候既是氣候變化的重要指示器[2]，也是外界環(huán)境變化的感應(yīng)器[3]。近年來(lái)，隨著遙感技術(shù)的進(jìn)步，基于遙感影像的植被物候監(jiān)測(cè)方法得到快速發(fā)展，相比于傳統(tǒng)的地面實(shí)測(cè)法，其不受時(shí)空限制，實(shí)現(xiàn)由小尺度向中大尺度的轉(zhuǎn)換。當(dāng)前，常用的遙感監(jiān)測(cè)數(shù)據(jù)主要有AVHRR、SPOT、MODIS和Landsat數(shù)據(jù)，其中，MODIS數(shù)據(jù)空間分辨率適宜，且具有穩(wěn)定的植被指數(shù)數(shù)據(jù)，被廣泛應(yīng)用于植被物候的研究中[4-7]。

城市是人類(lèi)活動(dòng)改變外界環(huán)境最明顯的場(chǎng)所，隨著城市化進(jìn)程的快速推進(jìn)，由城市熱島所引起的局地小氣候的改變，對(duì)生態(tài)系統(tǒng)的格局、過(guò)程和功能產(chǎn)生影響[8]。同時(shí)，城市景觀格局的改變，影響著景觀系統(tǒng)內(nèi)的物質(zhì)循環(huán)和能量流動(dòng)，這些都對(duì)城市植被的物候期產(chǎn)生影響。不少學(xué)者發(fā)現(xiàn)城市化導(dǎo)致城-鄉(xiāng)梯度上的物候差異，相對(duì)于農(nóng)村地區(qū)，城區(qū)植被生長(zhǎng)季開(kāi)始期更早、結(jié)束期更遲和長(zhǎng)度更長(zhǎng)[9-10]；Qiu等[11]采用不透水層百分比表示城市化水平，發(fā)現(xiàn)不透水層的增加導(dǎo)致植被生長(zhǎng)季開(kāi)始期提前、結(jié)束期推遲；Jeong等[12]發(fā)現(xiàn)韓國(guó)首爾地區(qū)人口密度的增加與春秋季物候期相關(guān)，導(dǎo)致該地區(qū)春季物候期提前、秋季物候期推遲。植被物候的改變對(duì)碳固定[13]、鳥(niǎo)類(lèi)的遷徙[14]、花粉過(guò)敏人群的健康[15]等產(chǎn)生影響，因此，研究城市化發(fā)展對(duì)植被物候的影響，對(duì)北京建設(shè)生態(tài)文明城市、保護(hù)生物多樣性等方面具有重要的現(xiàn)實(shí)意義。

本研究利用2001—2019年的MODIS EVI數(shù)據(jù)，提取北京的植被物候參數(shù)，即植被生長(zhǎng)季開(kāi)始期(start of season，SOS)、生長(zhǎng)季結(jié)束期(end of season，EOS)和生長(zhǎng)季長(zhǎng)度(length of season，LOS)，分析植被物候的時(shí)空變化特征及其對(duì)城市化的響應(yīng)。

1 研究區(qū)及數(shù)據(jù)

1.1 研究區(qū)概況

北京位于我國(guó)華北平原西北邊緣，在39.5°N～41°N，115.4°E～117.5°E之間，背靠太行山余脈和燕山山脈，面向華北平原，東南距渤海150 km，土地面積為16 410.54 km2，2019年常住人口為2 153.6萬(wàn)人。北京以山地和平原為主，地勢(shì)總體上西北高、東南低，屬于暖溫帶半濕潤(rùn)大陸性季風(fēng)氣候，四季分明，夏季高溫多雨，冬季寒冷干燥，年平均溫度為10～12 ℃，年降水量為600 mm，夏季盛行東南風(fēng)，冬季盛行西北風(fēng)，其自然植被類(lèi)型包含落葉林、混交林、灌叢、草地等(圖1)。作為全國(guó)最發(fā)達(dá)的城市之一，北京城市熱島效應(yīng)明顯、景觀格局變化顯著，這些變化改變了植被的生長(zhǎng)環(huán)境，進(jìn)而影響著其植被物候。

圖1 研究區(qū)

1.2 數(shù)據(jù)源

本研究使用的遙感數(shù)據(jù)為MODIS EVI數(shù)據(jù)(MOD13Q1)和土地利用數(shù)據(jù)。MOD13Q1數(shù)據(jù)集的時(shí)間跨度為2001—2019年，從美國(guó)航空航天局(National Aeronautics and Space Administration，NASA)數(shù)據(jù)共享平臺(tái)(http：//reverb.echo.nasa.gov/)獲取，其時(shí)間分辨率為16 d，空間分辨率為250 m，包含NDVI和EVI數(shù)據(jù)集。與NDVI相比，EVI能夠減弱大氣和地面反射率的影響以及弱化植被冠層背景變化的影響[16-17]，被認(rèn)為更適用于監(jiān)測(cè)稀疏植被覆蓋的城市植被動(dòng)態(tài)[18-19]。北京建筑物密集，植被相對(duì)較為稀疏，因而本研究選擇EVI提取城市植被物候信息。土地利用數(shù)據(jù)為中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心(http：//www.resdc.cn)提供的2000、2005、2010、2015和2018年的土地利用/覆蓋數(shù)據(jù)，用來(lái)提取城市建成區(qū)，其空間分辨率為1 000 m。該數(shù)據(jù)集以Landsat TM/ETM遙感影像為主要數(shù)據(jù)源，經(jīng)人工目視解譯生成，是當(dāng)前我國(guó)精度最高的土地利用遙感監(jiān)測(cè)數(shù)據(jù)產(chǎn)品，在國(guó)家水文、生態(tài)、土地資源調(diào)查研究中發(fā)揮重要作用。統(tǒng)計(jì)數(shù)據(jù)來(lái)源于《北京市統(tǒng)計(jì)年鑒》，包含人口密度、人均GDP和建成區(qū)面積。

2 研究方法

2.1 時(shí)間序列重建

雖然MODIS EVI數(shù)據(jù)集經(jīng)過(guò)標(biāo)準(zhǔn)化大氣糾正處理，但仍不能有效去除氣溶膠、云與地物雙向性反射引起的噪聲[20]，因此需對(duì)時(shí)序數(shù)據(jù)進(jìn)行擬合重建，剔除云層、大氣的干擾。本研究選用Savizky-Glolay(S-G)濾波法對(duì)EVI時(shí)間序列數(shù)據(jù)進(jìn)行擬合，該方法有效地保留了原始植被的特征，且更加適用于人類(lèi)活動(dòng)干擾更強(qiáng)的區(qū)域[21]。

S-G濾波法通過(guò)局部多項(xiàng)式回歸模型來(lái)實(shí)現(xiàn)平滑時(shí)序數(shù)據(jù)[22]，其基本思想是：基于高階多項(xiàng)式，利用最小二乘法對(duì)滑動(dòng)窗口內(nèi)數(shù)據(jù)進(jìn)行最佳擬合。

2.2 植被物候參數(shù)提取

當(dāng)前常用的植被物候提取法有導(dǎo)數(shù)法、閾值法和滑動(dòng)平均法[23]。其中，閾值法操作簡(jiǎn)單，使用靈活，包含固定閾值法和動(dòng)態(tài)閾值法。相對(duì)于固定閾值法，動(dòng)態(tài)閾值法可消除因土壤背景值和植被類(lèi)型差異帶來(lái)的影響，故本研究使用該方法提取植被物候參數(shù)。動(dòng)態(tài)閾值法是根據(jù)植被的生理特征，將植被指數(shù)曲線(xiàn)在上升和下降階段達(dá)到曲線(xiàn)振幅一定比例的時(shí)間作為植被的生長(zhǎng)季開(kāi)始期和結(jié)束期。Zhou等[24]在研究中國(guó)32個(gè)城市時(shí)，選擇0.2作為閾值提取植被物候；胡召玲等[25]在研究中國(guó)東北城市的植被物候時(shí)，選擇0.25作為閾值。本研究經(jīng)過(guò)嘗試，參考相關(guān)文獻(xiàn)并結(jié)合實(shí)際情況，將提取生長(zhǎng)季開(kāi)始期的閾值設(shè)置為0.2，提取生長(zhǎng)季結(jié)束期的閾值設(shè)置為0.25(式(1)、式(2))。

λSOS=(EMAX-EMIN1)0.2

(1)

λEOS=(EMAX-EMIN2)0.25

(2)

式中：λSOS表示提取生長(zhǎng)季開(kāi)始期的EVI閾值；EMAX表示EVI最大值；EMIN1表示曲線(xiàn)上升階段EVI最小值；λEOS表示提取生長(zhǎng)季結(jié)束期的EVI閾值；EMIN2表示曲線(xiàn)下降階段EVI最小值。

2.3 建成區(qū)及緩沖區(qū)

根據(jù)中國(guó)縣級(jí)行政區(qū)劃圖，利用土地利用/覆蓋數(shù)據(jù)提取北京主城區(qū)建成區(qū)，其中，2000、2005、2010、2015和2018年分別代表2001—2002年、2003—2007年、2008—2012年、2013—2017年、2018—2019年[26]。文獻(xiàn)[27-28]在研究北美東部與中國(guó)城市化對(duì)植被物候的影響中，估算城市化對(duì)植被物候的影響足跡距城市周邊不到20 km，據(jù)此，本研究選擇城區(qū)外20 km范圍作為郊區(qū)，并利用ArcGIS沿城區(qū)邊界向外延伸作0～1 km、1～2 km、2～5 km、5～10 km、10～15 km和15～20 km的緩沖區(qū)，得到圖2，并計(jì)算城區(qū)和郊區(qū)以及各緩沖區(qū)內(nèi)植被物候的平均值。

圖2 城區(qū)及緩沖區(qū)

3 結(jié)果與分析

3.1 植被物候空間分布

基于連續(xù)的MODIS EVI產(chǎn)品數(shù)據(jù)，使用動(dòng)態(tài)閾值法提取北京的植被物候期(生長(zhǎng)季開(kāi)始期、生長(zhǎng)季結(jié)束期和生長(zhǎng)季長(zhǎng)度)，得到北京植被物候期的空間分布圖(圖3)。生長(zhǎng)季開(kāi)始期空間分布圖(圖3(a))中，植被生長(zhǎng)季開(kāi)始期多集中在90～105 d，整體從東南向西北推遲，與海拔有一定關(guān)系，其中，城區(qū)植被的生長(zhǎng)季開(kāi)始期最早，集中在75～90 d，城區(qū)周邊植被生長(zhǎng)季開(kāi)始期較晚，多在90～105 d之間。生長(zhǎng)季結(jié)束期空間分布圖(圖3(b))中，植被生長(zhǎng)季結(jié)束期多集中在290～300 d，東南部向西北部逐漸提前，其中，城區(qū)植被生長(zhǎng)季結(jié)束期最晚，集中在300～320 d，城區(qū)周邊植被生長(zhǎng)季結(jié)束期較早，多在290～300 d之間。生長(zhǎng)季長(zhǎng)度空間分布圖(圖3(c))中，植被生長(zhǎng)季長(zhǎng)度多集中在190～210 d，且分布規(guī)律與生長(zhǎng)季開(kāi)始期和結(jié)束期相同，東南向西北遞減，同時(shí)，城區(qū)植被生長(zhǎng)季長(zhǎng)度最長(zhǎng)，周邊稍短，分別集中在210～230 d和190～210 d。

圖3 北京植被物候期多年平均值

3.2 植被物候年際變化

以年份為自變量，以北京每年生長(zhǎng)季開(kāi)始期、結(jié)束期和長(zhǎng)度的平均值為因變量，引入一元線(xiàn)性回歸方程擬合北京2001—2019年植被物候參數(shù)的變化趨勢(shì)(圖4)。由圖4發(fā)現(xiàn)，北京出現(xiàn)植被生長(zhǎng)季開(kāi)始期提前、結(jié)束期推遲和長(zhǎng)度延長(zhǎng)的現(xiàn)象，平均植被生長(zhǎng)季開(kāi)始期、結(jié)束期和長(zhǎng)度分別為96.2、295.5和199 d。圖4(a)～圖4(c)中，生長(zhǎng)季開(kāi)始期平均每年提前0.69 d，波動(dòng)幅度較大，其中，2006年的生長(zhǎng)季開(kāi)始期最晚，超過(guò)110 d，2017年的生長(zhǎng)季開(kāi)始期最早，小于90 d；生長(zhǎng)季結(jié)束期平均每年延遲0.71 d，波動(dòng)幅度相對(duì)較為穩(wěn)定；生長(zhǎng)季長(zhǎng)度平均每年延長(zhǎng)1.4 d，其中，2006年生長(zhǎng)季長(zhǎng)度最短，2017年最長(zhǎng)。由生長(zhǎng)季開(kāi)始期和結(jié)束期的波動(dòng)幅度差異發(fā)現(xiàn)，生長(zhǎng)季開(kāi)始期對(duì)外界環(huán)境的變化更為敏感。

圖4 北京植被物候期年際變化趨勢(shì)

圖5為城區(qū)和郊區(qū)的植被物候期變化趨勢(shì)。相比于郊區(qū)的植被物候期，每年城區(qū)都表現(xiàn)出更早的生長(zhǎng)季開(kāi)始期、更晚的生長(zhǎng)季結(jié)束期和更長(zhǎng)的生長(zhǎng)季長(zhǎng)度。由圖5(a)可見(jiàn)，2001—2019年，城區(qū)和郊區(qū)的植被生長(zhǎng)季開(kāi)始期均呈現(xiàn)提前趨勢(shì)，但提前速度有所差異，城區(qū)植被的生長(zhǎng)季開(kāi)始期平均每年提前0.23 d，郊區(qū)提前速度較快，為0.68 d/a；圖5(b)中，城區(qū)與郊區(qū)均表現(xiàn)出植被生長(zhǎng)季結(jié)束期推遲的現(xiàn)象，城區(qū)的推遲速度較慢，平均每年推遲0.16 d，郊區(qū)較快，以平均每年0.72 d的速度推遲；圖5(c)中，城區(qū)和郊區(qū)的生長(zhǎng)季長(zhǎng)度都呈延長(zhǎng)趨勢(shì)，平均每年分別延長(zhǎng)0.42、1.42 d。整體來(lái)看，城區(qū)和郊區(qū)的植被物候期均呈現(xiàn)出相同的變化趨勢(shì)，但變化速度不同，郊區(qū)的植被物候期變化速度明顯快于城區(qū)。

圖5 北京城區(qū)和郊區(qū)植被物候期年際變化趨勢(shì)

3.3 城市化梯度上的植被物候變化

北京的植被物候期在城郊方向上表現(xiàn)出明顯的梯度現(xiàn)象(圖6)。由圖6可見(jiàn)，2001—2019年間，植被物候期與距城區(qū)的距離呈顯著的對(duì)數(shù)關(guān)系，城區(qū)的植被生長(zhǎng)季開(kāi)始期最早、結(jié)束期最晚和長(zhǎng)度最長(zhǎng)，分別為86.6 d、309.3 d和222 d。隨著距城區(qū)距離的增加，植被生長(zhǎng)季開(kāi)始期逐漸變晚、結(jié)束期逐漸提前、生長(zhǎng)季長(zhǎng)度逐漸縮短，同時(shí)，距離城區(qū)越遠(yuǎn)，植被物候期的變化幅度越小。總體來(lái)看，植被生長(zhǎng)季開(kāi)始期、結(jié)束期和長(zhǎng)度均在0～2 km范圍內(nèi)變化明顯，這歸因于城市近郊區(qū)的熱島強(qiáng)度和不透水層覆蓋率遠(yuǎn)大于遠(yuǎn)郊區(qū)，人類(lèi)活動(dòng)頻繁，城市化對(duì)植被物候期影響顯著。隨著距城區(qū)距離的增加，植被物候期變化的幅度越小，10 km以外的物候期變化趨勢(shì)趨于平緩，這表明該范圍外的植被物候期受城市化影響明顯減弱。

圖6 城郊梯度上的植被物候期變化趨勢(shì)

3.4 城市化特征分析

城市化也稱(chēng)城鎮(zhèn)化，是一個(gè)復(fù)雜的進(jìn)程，是社會(huì)由以農(nóng)業(yè)為主向以工業(yè)和服務(wù)業(yè)為主逐漸轉(zhuǎn)變的過(guò)程，包括人口城市化、經(jīng)濟(jì)城市化和土地城市化。本研究選擇人口密度表示人口城市化、人均GDP表示經(jīng)濟(jì)城市化、建成區(qū)面積占比表示土地城市化，以這三個(gè)指標(biāo)定義城市化進(jìn)程。2001—2019年北京的人口密度、人均GDP和建成區(qū)面積占比的年際變化趨勢(shì)如圖7所示。圖7(a)中，2001—2019年北京人口密度整體呈增長(zhǎng)趨勢(shì)，由2001年的824人/km2增長(zhǎng)至2019年的1 312人/km2，其年均增長(zhǎng)率為2.6%，其中，2001—2014年人口密度增速較快，年均增長(zhǎng)率為3.6%，2014—2017年趨于平穩(wěn)，因2017年后北京開(kāi)始整頓低端產(chǎn)業(yè)，人口稍有下降。2001—2019年人均GDP以10.3%的年均增長(zhǎng)率穩(wěn)步增長(zhǎng)，由2001年的28 097元增長(zhǎng)至2019年的165 220元，見(jiàn)圖7(b)，其中，2015年后增速最快。建成區(qū)面積占比呈波動(dòng)增長(zhǎng)趨勢(shì)，見(jiàn)圖7(c)，年均增長(zhǎng)率為3.5%，2009—2010年建成區(qū)面積占比有所下降，其他時(shí)間段均在增長(zhǎng)。整體而言，北京的城市化進(jìn)程明顯，其人口密度、人均GDP和建成區(qū)面積占比均呈現(xiàn)增長(zhǎng)趨勢(shì)。

圖7 城市化指標(biāo)年際變化趨勢(shì)

3.5 城市化對(duì)植被物候的影響

對(duì)北京城區(qū)和郊區(qū)的植被物候期與人口密度、人均GDP和建成區(qū)面積占比進(jìn)行皮爾遜(Pearson)相關(guān)性分析，得到表1。就城區(qū)來(lái)看，植被生長(zhǎng)季開(kāi)始期與城市化的三個(gè)指標(biāo)均無(wú)明顯相關(guān)關(guān)系(P>0.05)，植被生長(zhǎng)季結(jié)束期和長(zhǎng)度均與人均GDP存在顯著正相關(guān)關(guān)系(P<0.05)，人均GDP的增長(zhǎng)對(duì)植被生長(zhǎng)季結(jié)束期推遲和長(zhǎng)度延長(zhǎng)有重要作用。而郊區(qū)的植被物候期均與人口密度、人均GDP和建成區(qū)面積存在顯著關(guān)系，其中，生長(zhǎng)季開(kāi)始期與人口密度(P<0.05)、人均GDP(P<0.01)和建成區(qū)面積占比(P<0.05)均顯著負(fù)相關(guān)，生長(zhǎng)季結(jié)束期、生長(zhǎng)季長(zhǎng)度與人口密度、人均GDP和建成區(qū)面積占比均極顯著正相關(guān)(P<0.01)，城市化進(jìn)程對(duì)郊區(qū)的植被物候的影響顯著。由城區(qū)和郊區(qū)的植被物候期與各城市化因子的相關(guān)關(guān)系可以發(fā)現(xiàn)，郊區(qū)植被物候期受城市化影響更為劇烈，這可能因?yàn)槌鞘谢M(jìn)程中，北京城市擴(kuò)張明顯，郊區(qū)建筑物快速增加，且大型工廠多在郊區(qū)，這些都導(dǎo)致郊區(qū)水熱環(huán)境變化顯著，深刻影響到植被物候。

表1 植被物候期與人口密度、人均GDP和建成區(qū)面積占比的皮爾遜相關(guān)系數(shù)

4 討論

本研究使用人口密度、人均GDP和建成區(qū)面積占比表示城市化進(jìn)程，探討城市化進(jìn)程對(duì)植被物候的影響，發(fā)現(xiàn)植被物候期與城市化進(jìn)程相關(guān)性明顯，這與文獻(xiàn)[11]、文獻(xiàn)[12]的研究結(jié)果相一致。此外，植被物候期除受人類(lèi)活動(dòng)影響外，還受到氣候變化的影響，文獻(xiàn)[29]根據(jù)研究青藏高原東北部牧草物候發(fā)現(xiàn)，影響牧草物候期的主要因素是溫度，降水次之；文獻(xiàn)[11]對(duì)中緯度大城市的植被物候進(jìn)行研究，發(fā)現(xiàn)在大部分氣候區(qū)，不透水層的增加對(duì)植被生長(zhǎng)季開(kāi)始期提前、結(jié)束期推遲有重要作用，同時(shí)，氣候因素對(duì)植被物候變化起主導(dǎo)作用；Yuan等[30]發(fā)現(xiàn)黃河流域植被返青期受季前日最低溫度影響更大，且城市熱島造成城鄉(xiāng)植被物候的差異。

城區(qū)與郊區(qū)的植被物候變化程度有所不同，郊區(qū)的植被物候變化更為劇烈，這可能與北京城區(qū)人口流出、郊區(qū)人口暴增、郊區(qū)發(fā)展迅速有關(guān)。2010—2019年間，北京主城區(qū)(東城區(qū)、西城區(qū)、朝陽(yáng)區(qū)、豐臺(tái)區(qū)、石景山區(qū)、海淀區(qū))人口增加445萬(wàn)人，而近郊區(qū)(房山區(qū)、通州區(qū)、順義區(qū)、昌平區(qū)、大興區(qū)、門(mén)頭溝區(qū))人口增加546萬(wàn)人。其中，2000—2019年，主城區(qū)流出48萬(wàn)人，近郊區(qū)暴增約223萬(wàn)人，郊區(qū)人口增長(zhǎng)更加明顯，同時(shí)，北京城市化發(fā)展不斷向郊區(qū)擴(kuò)張，這些導(dǎo)致郊區(qū)植被生長(zhǎng)的外部環(huán)境變化顯著，植被物候變化劇烈。

5 結(jié)束語(yǔ)

本研究基于2001—2019年MOD13Q1植被指數(shù)產(chǎn)品，利用S-G濾波算法和動(dòng)態(tài)閾值法提取了北京的植被生長(zhǎng)季開(kāi)始期、結(jié)束期和長(zhǎng)度，研究其植被物候期的時(shí)間變化趨勢(shì)和空間分布特征，并探究城市化對(duì)城區(qū)和郊區(qū)植被物候的影響，研究結(jié)論如下。

1)北京出現(xiàn)植被生長(zhǎng)季開(kāi)始期提前、結(jié)束期推遲和長(zhǎng)度延長(zhǎng)的現(xiàn)象。分區(qū)域看，城區(qū)和郊區(qū)也均表現(xiàn)相同的現(xiàn)象，但城區(qū)和郊區(qū)的植被物候期變化速度有所差異，郊區(qū)植被物候變化更為明顯。

2)植被物候期在城郊梯度上表現(xiàn)出明顯的空間異質(zhì)性，且植被物候期與距城區(qū)的距離呈顯著的對(duì)數(shù)關(guān)系，城區(qū)的生長(zhǎng)季開(kāi)始期最早、結(jié)束期最遲、長(zhǎng)度最長(zhǎng)。

3)植被物候期與城市化進(jìn)程相關(guān)性明顯。分區(qū)域看，郊區(qū)的植被物候期與城市化因子(人口密度、人均GDP和建成區(qū)面積占比)相關(guān)性更強(qiáng)，城區(qū)稍弱。