北京植被物候時空變化及其對城市化的響應

卜亞勤，丁海勇

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

0 引言

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

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

本研究利用2001—2019年的MODIS EVI數據，提取北京的植被物候參數，即植被生長季開始期(start of season，SOS)、生長季結束期(end of season，EOS)和生長季長度(length of season，LOS)，分析植被物候的時空變化特征及其對城市化的響應。

1 研究區及數據

1.1 研究區概況

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

圖1 研究區

1.2 數據源

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

2 研究方法

2.1 時間序列重建

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

S-G濾波法通過局部多項式回歸模型來實現平滑時序數據[22]，其基本思想是：基于高階多項式，利用最小二乘法對滑動窗口內數據進行最佳擬合。

2.2 植被物候參數提取

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

λSOS=(EMAX-EMIN1)0.2

(1)

λEOS=(EMAX-EMIN2)0.25

(2)

式中：λSOS表示提取生長季開始期的EVI閾值；EMAX表示EVI最大值；EMIN1表示曲線上升階段EVI最小值；λEOS表示提取生長季結束期的EVI閾值；EMIN2表示曲線下降階段EVI最小值。

2.3 建成區及緩沖區

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

圖2 城區及緩沖區

3 結果與分析

3.1 植被物候空間分布

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

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

3.2 植被物候年際變化

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

圖4 北京植被物候期年際變化趨勢

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

圖5 北京城區和郊區植被物候期年際變化趨勢

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

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

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

3.4 城市化特征分析

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

圖7 城市化指標年際變化趨勢

3.5 城市化對植被物候的影響

對北京城區和郊區的植被物候期與人口密度、人均GDP和建成區面積占比進行皮爾遜(Pearson)相關性分析，得到表1。就城區來看，植被生長季開始期與城市化的三個指標均無明顯相關關系(P>0.05)，植被生長季結束期和長度均與人均GDP存在顯著正相關關系(P<0.05)，人均GDP的增長對植被生長季結束期推遲和長度延長有重要作用。而郊區的植被物候期均與人口密度、人均GDP和建成區面積存在顯著關系，其中，生長季開始期與人口密度(P<0.05)、人均GDP(P<0.01)和建成區面積占比(P<0.05)均顯著負相關，生長季結束期、生長季長度與人口密度、人均GDP和建成區面積占比均極顯著正相關(P<0.01)，城市化進程對郊區的植被物候的影響顯著。由城區和郊區的植被物候期與各城市化因子的相關關系可以發現，郊區植被物候期受城市化影響更為劇烈，這可能因為城市化進程中，北京城市擴張明顯，郊區建筑物快速增加，且大型工廠多在郊區，這些都導致郊區水熱環境變化顯著，深刻影響到植被物候。

表1 植被物候期與人口密度、人均GDP和建成區面積占比的皮爾遜相關系數

4 討論

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

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

5 結束語

本研究基于2001—2019年MOD13Q1植被指數產品，利用S-G濾波算法和動態閾值法提取了北京的植被生長季開始期、結束期和長度，研究其植被物候期的時間變化趨勢和空間分布特征，并探究城市化對城區和郊區植被物候的影響，研究結論如下。

1)北京出現植被生長季開始期提前、結束期推遲和長度延長的現象。分區域看，城區和郊區也均表現相同的現象，但城區和郊區的植被物候期變化速度有所差異，郊區植被物候變化更為明顯。

2)植被物候期在城郊梯度上表現出明顯的空間異質性，且植被物候期與距城區的距離呈顯著的對數關系，城區的生長季開始期最早、結束期最遲、長度最長。

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