張鵬騫 ，胡理樂(lè)，白加德 *

1.北京麋鹿生態(tài)實(shí)驗(yàn)中心，北京 100076；2.北京生物多樣性保護(hù)研究中心，北京 100076；3.北京市園林綠化局，北京 100013；4.中國(guó)環(huán)境科學(xué)研究院，北京 100012

京津冀城市群是中國(guó)重要的政治和經(jīng)濟(jì)中心，推動(dòng)京津冀協(xié)同一體化發(fā)展是一個(gè)重大國(guó)家戰(zhàn)略。京津冀地區(qū)在快速城市化的過(guò)程中，實(shí)現(xiàn)了區(qū)域經(jīng)濟(jì)快速增長(zhǎng)，然而人類活動(dòng)強(qiáng)度的增加也給區(qū)域自然生態(tài)系統(tǒng)帶來(lái)了巨大的壓力。健康穩(wěn)定的自然生態(tài)系統(tǒng)和自然資源是支撐城市化發(fā)展的重要物質(zhì)基礎(chǔ)（閻世杰等，2019），京津冀地區(qū)內(nèi)有京津水源涵養(yǎng)生態(tài)功能保護(hù)區(qū)，在全國(guó)生態(tài)功能區(qū)劃中占有重要地位，因此加強(qiáng)京津冀地區(qū)的生態(tài)環(huán)境保護(hù)至關(guān)重要。近幾年，為響應(yīng)一體化協(xié)同可持續(xù)發(fā)展戰(zhàn)略，京津冀地區(qū)采取退耕還林、平原造林以及在城區(qū)建設(shè)綠地等一系列措施，生態(tài)環(huán)境已得到初步改善（孟丹等，2015），但生態(tài)安全形勢(shì)仍然嚴(yán)峻（劉軍會(huì)等，2018）。對(duì)京津冀地區(qū)生態(tài)環(huán)境狀況的評(píng)估以及如何有效地開(kāi)展生態(tài)保護(hù)和修復(fù)、如何進(jìn)行城市規(guī)劃，仍然是京津冀一體化協(xié)同發(fā)展的研究重點(diǎn)（鄧越等，2018；王喆等，2015）。

植被是生態(tài)環(huán)境中一個(gè)極其重要的組成部分，地表植被變化能直接反映區(qū)域的生態(tài)環(huán)境變化，尤其是對(duì)于京津冀地區(qū)這種因城市化而導(dǎo)致景觀破碎加劇的地區(qū)，地表植被變化可以很好地表征當(dāng)?shù)氐纳鷳B(tài)環(huán)境質(zhì)量（王靜等，2017）。歸一化植被指數(shù)（NDVI，normalization difference vegetation index）是監(jiān)測(cè)地區(qū)或全球植被變化的有效指標(biāo)之一，是植物生長(zhǎng)狀況及植被覆蓋度的最佳指示因子（孫智輝等，2010），并在全球范圍被廣泛應(yīng)用（Cai et al.，2014；Eastman et al.，2013；Tucker et al.，2005）。國(guó)內(nèi)學(xué)者也對(duì)遼西（王亭月，2020）、華北（孫艷玲等，2012）、陜北（程杰等，2020）、吉林（王景發(fā)，2020）、西南（曹云等，2020）等地區(qū)的NDVI時(shí)空變化特征以及NDVI對(duì)氣溫和降水的響應(yīng)（盧喬倩等，2020）進(jìn)行了研究，這些研究對(duì)評(píng)估當(dāng)?shù)厣鷳B(tài)環(huán)境狀況和制定可持續(xù)發(fā)展戰(zhàn)略提供了可靠支撐。

《京津冀協(xié)同規(guī)劃綱要》和《北京城市總體規(guī)劃（2016—2035年）》中將京津冀分為中部核心功能區(qū)、東部濱海發(fā)展區(qū)、南部功能拓展區(qū)和西北部生態(tài)涵養(yǎng)區(qū)4個(gè)功能區(qū)，每個(gè)功能區(qū)都有明確的空間范圍和發(fā)展重點(diǎn)。目前，已有學(xué)者基于遙感數(shù)據(jù)對(duì)京津冀地區(qū)的NDVI及影響因素進(jìn)行了研究（王靜等，2017；閻世杰等，2019），但在這些研究中，通常將京津冀地區(qū)作為一個(gè)整體分析，缺少各個(gè)功能區(qū)的對(duì)比分析。然而不同的功能區(qū)有不同的定位和發(fā)展模式，植被狀況也存在一定差異，分區(qū)分析京津冀地區(qū)的 NDVI，可反映在不同城市發(fā)展模式下植被狀況的變化，也可為各功能區(qū)制定具體發(fā)展措施提供科學(xué)參考。

本研究基于京津冀地區(qū)的功能分區(qū)，使用Google earth engine（GEE）遙感數(shù)據(jù)分析平臺(tái)和MOD13Q1-NDVI數(shù)據(jù)、Landsat Surface Reflectance（Landsat-SR）數(shù)據(jù)，對(duì)2000—2019年京津冀地區(qū)和北京市的NDVI時(shí)空變化特征及分布格局進(jìn)行了分析，以期反映京津冀各功能區(qū)的生態(tài)壞境狀況，為各功能區(qū)的生態(tài)保護(hù)修復(fù)和城市發(fā)展策略提供數(shù)據(jù)支撐與科學(xué)依據(jù)，對(duì)保護(hù)京津水源涵養(yǎng)生態(tài)功能保護(hù)區(qū)以及京津冀地區(qū)的綠色可持續(xù)發(fā)展有重要意義。

1 研究區(qū)與研究方法

1.1 研究區(qū)概況

京津冀地區(qū)地處華北平原北部（113°27′—119°50′E，36°03′—42°40′N），東臨渤海灣，西靠太行山，南為華北平原，北接壤燕山山脈。行政上包括北京和天津兩個(gè)直轄市以及河北省所轄的 13個(gè)地級(jí)市，區(qū)域土地總面積約2.18×105km2，常駐人口約1.1億，是北方經(jīng)濟(jì)規(guī)模最大的城市群。該地區(qū)東南部地勢(shì)平坦，西北部地勢(shì)較高，包含高原、山地、丘陵、平原、盆地等地貌類型，屬溫帶大陸性季風(fēng)氣候，春季干旱多風(fēng)，夏季高溫多雨，冬季寒冷干燥，年平均氣溫為 0—13 ℃，年降水量為300—800 mm，年無(wú)霜期為120—200 d。區(qū)域內(nèi)高原區(qū)以多年旱生草本植物為主，丘陵區(qū)以灌木為主，山間盆地區(qū)以喬木和灌草叢為主，平原區(qū)以耕地和人工植被為主。該地區(qū)作為中國(guó)的政治經(jīng)濟(jì)文化中心，具有重要的發(fā)展戰(zhàn)略地位，目前共建設(shè)有18個(gè)國(guó)家級(jí)自然保護(hù)區(qū)，區(qū)域北部從屬三北防護(hù)林的華北段，內(nèi)部還有京津水源涵養(yǎng)生態(tài)功能保護(hù)區(qū)，具有重大的生態(tài)安全意義（程琳琳等，2019）。根據(jù)《京津冀協(xié)同規(guī)劃綱要》和《北京城市總體規(guī)劃（2016—2035年）》，京津冀地區(qū)分為西北部生態(tài)涵養(yǎng)區(qū)、中部核心功能區(qū)、南部功能拓展區(qū)以及東部濱海發(fā)展區(qū)4個(gè)功能區(qū)（圖1）。

圖1 京津冀地區(qū)功能分區(qū)圖Fig.1 Functional zoning map of Beijing-Tianjin-Hebei region

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

1.2.1 數(shù)據(jù)來(lái)源

本研究分析了京津冀中部核心功能區(qū)、東部濱海發(fā)展區(qū)、南部功能拓展區(qū)和西北部生態(tài)涵養(yǎng)區(qū) 4個(gè)功能區(qū)和北京市平原區(qū)、淺山區(qū)、深山區(qū)3區(qū)近20年（2000—2019年）的NDVI變化特征，京津冀地區(qū)NDVI數(shù)據(jù)根據(jù)MOD13Q1影像計(jì)算得到，空間分辨率為 250 m；北京市 NDVI數(shù)據(jù)根據(jù)Landsat Surface Reflectance（Landsat-SR）影像計(jì)算得到，空間分辨率為30 m。影像的調(diào)用和處理均在Google Earth Engine（GEE）云計(jì)算地理信息處理平臺(tái)中實(shí)現(xiàn)。京津冀功能分區(qū)數(shù)據(jù)來(lái)自《北京城市總體規(guī)劃（2016—2035年）》。北京市平原區(qū)、淺山區(qū)和深山區(qū)三區(qū)數(shù)據(jù)來(lái)自《北京市淺山區(qū)保護(hù)規(guī)劃（2017—2035年）》，分別是市域內(nèi)海拔低于100 m、100—300 m和高于300 m的區(qū)域。京津冀市級(jí)行政區(qū)劃矢量數(shù)據(jù)來(lái)自中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http://www.resdc.cn）。

1.2.2 歸一化植被指數(shù)（NDVI）數(shù)據(jù)提取

NDVI是植物生長(zhǎng)狀態(tài)和植被空間分布的指示因子，與地表植被的覆蓋率呈正比關(guān)系，對(duì)于同一種植被，NDVI值越大，表明地表植被覆蓋率越高，植被長(zhǎng)勢(shì)越好。本研究使用Landsat ETM+數(shù)據(jù)，提取北京市2000—2019年NDVI數(shù)據(jù)；使用Modis植被指數(shù)產(chǎn)品提取京津冀地區(qū)2000—2019年NDVI數(shù)據(jù)。

式中，ρNIR為近紅外波段反射率；ρRED為紅光波段反射率。

使用研究區(qū)域夏季（6—8月）NDVI均值來(lái)表征該地區(qū)當(dāng)年的植被覆蓋狀況。公式如下：

1.2.3 NDVI變化趨勢(shì)分析

為去除NDVI數(shù)據(jù)中因云和大氣影響而產(chǎn)生的誤差和噪聲，使用S-G濾波算法通過(guò)在滑動(dòng)窗口內(nèi)對(duì)給定高階多項(xiàng)式進(jìn)行最小二乘擬合，對(duì)時(shí)間序列數(shù)據(jù)進(jìn)行重構(gòu)，以提高數(shù)據(jù)質(zhì)量（黃耀歡等，2009）。采用最小二乘法線性擬合模擬每個(gè)像元年NDVI值的變化趨勢(shì)，分析植被在研究時(shí)段的時(shí)空變化規(guī)律。計(jì)算歸一化植被指數(shù)序列與時(shí)間序列的像元回歸斜率Slope，若Slope為正則表示NDVI隨時(shí)間變化呈上升趨勢(shì)，為負(fù)則表示NDVI隨時(shí)間呈下降趨勢(shì)，Slope的絕對(duì)值越大，植被覆蓋狀況變化越明顯。其計(jì)算公式為：

式中，Slope為像元NDVI變化趨勢(shì)斜率；j為監(jiān)測(cè)時(shí)間段內(nèi)的年序號(hào)；NDVIj為第j年的夏季NDVI均值；n為監(jiān)測(cè)時(shí)間段的年數(shù)。本研究將近20年分為2000—2009年和2010—2019年兩個(gè)時(shí)間段，分別計(jì)算近20年NDVI總的變化趨勢(shì)以及分段變化趨勢(shì)。

2 結(jié)果與分析

2.1 京津冀地區(qū)總體NDVI時(shí)空變化特征

京津冀地區(qū)的NDVI在前4年（2000—2004年）增長(zhǎng)明顯，隨后基本保持不變，略有波動(dòng)（圖2），其中在2004、2008、2013、2016、2018年為峰值，2006、2009、2014、2017年為谷值。20年間（2000—2019年）NDVI的變化區(qū)間為0.529—0.632，最小值出現(xiàn)在2000年，最大值出現(xiàn)在2008年。整體來(lái)看，NDVI波動(dòng)增長(zhǎng)，增速為0.025/(10 a)，植被狀況逐年變好。

圖2 京津冀地區(qū)總體NDVI年變化Fig.2 Annual change of overall NDVI in Beijing-Tianjin-Hebei region

京津冀地區(qū)NDVI空間分布差異明顯（圖3），分布特征為西北低東南高，和京津冀地區(qū)的海拔分布恰好相反，海拔高的地區(qū)，NDVI普遍較低。除此之外，北京、天津等城市的中心城區(qū)以及張家口市NDVI指數(shù)明顯低于其他地區(qū)，植被狀況較差，但也在逐年改善。

圖3 京津冀地區(qū)2000、2009、2019年NDVI空間分布Fig.3 Spatial distribution of NDVI in Beijing-Tianjin-Hebei region in 2000, 2009 and 2019

2.2 京津冀4個(gè)功能區(qū)NDVI時(shí)空變化特征

京津冀4個(gè)功能區(qū)的NDVI均值年變化有一定差異（圖4），西北部生態(tài)涵養(yǎng)區(qū)NDVI均值最高，其次是南部功能拓展區(qū)和中部核心功能區(qū)，東部濱海發(fā)展區(qū)NDVI均值最低。整體來(lái)看，前4年（2000—2004年），4個(gè)功能區(qū)的NDVI均明顯增長(zhǎng)。隨后，NDVI在西北部生態(tài)涵養(yǎng)區(qū)波動(dòng)增長(zhǎng)，于2007、2009、2014、2017年出現(xiàn)谷值，于 2019年達(dá)到最大值；在其他3個(gè)功能區(qū)波動(dòng)變化，沒(méi)有增長(zhǎng)，在南部功能拓展區(qū)于2006、2009、2014年出現(xiàn)谷值，在中部核心功能區(qū)于2006、2010、2014年出現(xiàn)谷值，在東部濱海發(fā)展區(qū)于 2001、2014、2017年出現(xiàn)谷值。2008年，NDVI在西北部生態(tài)涵養(yǎng)區(qū)、南部功能拓展區(qū)和中部核心功能區(qū)均出現(xiàn)明顯的峰值，在東部濱海發(fā)展區(qū)出現(xiàn)較小的峰值。

圖4 京津冀不同區(qū)域NDVI年變化Fig.4 Annual changes of NDVI in different regions of Beijing-Tianjin-Hebei

NDVI在4個(gè)功能區(qū)的空間分布有明顯差異（圖3）。2000年，各區(qū)NDVI差異較大，但4個(gè)功能區(qū)中均有NDVI的高值和低值集中的地區(qū)。NDVI較高的區(qū)域分布在西北部生態(tài)涵養(yǎng)區(qū)中的承德市、北京市北部和西部、保定市西部和秦皇島市北部區(qū)域，其次是南部功能拓展區(qū)中的衡水市、邢臺(tái)市和邯鄲市大部分地區(qū)。NDVI較低的區(qū)域分布在西北部生態(tài)涵養(yǎng)區(qū)中的張家口市、承德市西北部，東部濱海發(fā)展區(qū)中天津市、滄州市的濱海地區(qū)，以及中部核心功能區(qū)的大部分地區(qū)。2009年，NDVI低值地區(qū)整體減少。西北部生態(tài)涵養(yǎng)區(qū)中除張家口市NDVI較低外，其他區(qū)域植被覆蓋情況良好。中部核心功能區(qū)的NDVI較2000年相比整體有所增加，其中北京市和天津市城市中心的NDVI明顯低于其他區(qū)域。東部濱海發(fā)展區(qū)中滄州市的NDVI有了明顯增長(zhǎng)，但濱海區(qū)的NDVI指數(shù)仍明顯低于其他地區(qū)。南部功能拓展區(qū)植被狀況整體較好，NDVI低值地區(qū)最少。2019年，西北部生態(tài)涵養(yǎng)區(qū)和其他城市中心的NDVI明顯提高，植被狀況明顯改善，其他地區(qū)NDVI的空間分布較之前沒(méi)有明顯變化。

近20年（2000—2019年），NDVI的變化趨勢(shì)在4個(gè)功能區(qū)表現(xiàn)為整體改善，局部退化，改善的區(qū)域明顯多于退化的區(qū)域（圖5，表1）。其中南部功能拓展區(qū)、中部核心功能區(qū)和東部濱海發(fā)展區(qū)的改善區(qū)域占比都在90%以上，改善趨勢(shì)比西北部生態(tài)涵養(yǎng)區(qū)明顯。植被退化的區(qū)域主要分布在西北部生態(tài)涵養(yǎng)區(qū)的承德市西北部、張家口市東北部以及北京市中心地區(qū)，植被改善明顯的區(qū)域主要分布在南部功能拓展區(qū)和中部核心功能區(qū)的南部。將2000—2019年分為兩個(gè)時(shí)間段，前10年（2000—2009年）各功能區(qū)植被改善區(qū)域的占比都比后 10年（2010—2019年）多，但差異較小，除西北部生態(tài)涵養(yǎng)區(qū)外，差異都在1%以內(nèi)。后10年，張家口市的植被改善區(qū)域明顯增加，但承德市北部增加了一些植被退化的區(qū)域，其他3個(gè)功能區(qū)的NDVI變化趨勢(shì)分布較上一個(gè)10年沒(méi)有明顯變化。

圖5 京津冀地區(qū)NDVI變化趨勢(shì)空間分布Fig.5 Spatial distribution of NDVI change trend in Beijing-Tianjin-Hebei region

表1 京津冀地區(qū)NDVI變化趨勢(shì)面積比例統(tǒng)計(jì)Table 1 Statistics on the proportion of NDVI change area in Beijing-Tianjin-Hebei region

2.3 北京NDVI時(shí)空變化特征

《北京市淺山區(qū)保護(hù)規(guī)劃（2017—2035年）》中將北京市分為平原區(qū)、淺山區(qū)、深山區(qū)3區(qū)，本研究使用 30 m的 Landsat Surface Reflectance（Landsat-SR）影像集對(duì)北京市3區(qū)的NDVI數(shù)據(jù)進(jìn)行了統(tǒng)計(jì)分析（圖6）。北京市深山區(qū)的NDVI均值最高，其次是淺山區(qū)，平原區(qū)最低，這3個(gè)區(qū)域的NDVI時(shí)間變化趨勢(shì)略有不同，整體呈波動(dòng)上升趨勢(shì)。其中，淺山區(qū)NDVI增長(zhǎng)最明顯（k=0.120/(10 a)），20年平均增長(zhǎng)了 51.95%，深山區(qū)和平原區(qū)的增速分別為0.106/(10 a)和0.073/(10 a)。NDVI的波動(dòng)趨勢(shì)在3個(gè)區(qū)域基本一致，于2003年和2012年均出現(xiàn)明顯的谷值。2007年，平原區(qū)和淺山區(qū)的NDVI出現(xiàn)了谷值，而深山區(qū)的NDVI在這一年略有增加；2015年，深山區(qū)的NDVI出現(xiàn)明顯的谷值，平原區(qū)和淺山區(qū)的NDVI基本保持不變。

圖6 北京不同區(qū)域NDVI年變化Fig.6 Annual changes of NDVI in different regions of Beijing

3 討論

京津冀地區(qū)近20年來(lái)NDVI增加是因?yàn)?000年以來(lái)采取的林草地植被保護(hù)、小流域生態(tài)治理和高原地區(qū)水土保持等措施取得了成效（吳建生等，2015）。尤其是2008年為籌備奧運(yùn)會(huì)，京津冀地區(qū)圍繞北京市開(kāi)展綠化美化等一系列生態(tài)工程來(lái)改善生態(tài)環(huán)境質(zhì)量，這使得2008年NDVI明顯增加，植被狀況明顯改善（孫曉鵬等，2012）。2009年NDVI出現(xiàn)明顯的下降是因?yàn)樵摰貐^(qū)遭受了接近歷史極限的干旱災(zāi)害，對(duì)植被造成了嚴(yán)重的影響（陳權(quán)亮等，2010；李卓等，2017）。京津冀地區(qū) NDVI較高的區(qū)域主要分布西北部生態(tài)涵養(yǎng)區(qū)，該區(qū)從屬三北防護(hù)林的華北段，作為阻止風(fēng)沙入侵京津兩中心城市的重要屏障，該區(qū)域采取了建造木草結(jié)合的防護(hù)網(wǎng)等一系列防護(hù)措施，故近 20年該區(qū)域大部分地區(qū)的 NDVI穩(wěn)定在較高值（孟丹等，2015）。但該區(qū)域內(nèi)的高原地區(qū)水熱條件較差，仍存在土地荒漠化和水土流失等問(wèn)題（李卓等，2017），這可能是張家口市和承德市西北部植被退化的原因。中部核心功能區(qū)、南部功能拓展區(qū)和東部濱海發(fā)展區(qū)位于京津冀平原地區(qū)，海拔較低，氣候溫暖潮濕，良好的氣候因素以及平原造林和沿海防護(hù)林工程使這3個(gè)區(qū)域大部分地區(qū)的植被得到改善，生態(tài)環(huán)境狀況明顯好轉(zhuǎn)（潘梅等，2020）。除張家口市外，NDVI較低的區(qū)域主要集中在北京市、天津市、石家莊市、衡水市和邯鄲市等城市中心區(qū)域，此類區(qū)域城市化進(jìn)程較快，大量林地和農(nóng)用耕地轉(zhuǎn)化為城市用地，導(dǎo)致城市及周邊地區(qū)植被覆蓋減少，低NDVI區(qū)域的面積擴(kuò)大，影響到生態(tài)環(huán)境質(zhì)量（王靜等，2014；Xin et al.，2008）。近20年，政府重視生態(tài)文明建設(shè)，在這些城市中心建設(shè)了大量綠地來(lái)吸收大氣污染和降低熱島效應(yīng)，同時(shí)植被狀況得到改善（Wu，2013；Zhou et al.，2014），這也表明政府對(duì)生態(tài)與城市化和諧發(fā)展的重視對(duì)生態(tài)環(huán)境改善有正面影響并且成效明顯。

北京市近20年的NDVI波動(dòng)上升，空間上表現(xiàn)為深山區(qū)＞淺山區(qū)＞平原區(qū)。深山區(qū)主要分布在北京市的西北部，屬溫帶半濕潤(rùn)季風(fēng)氣候，海拔較高，受人類活動(dòng)影響相對(duì)較小，且 2000年以來(lái)采取風(fēng)沙源治理等多種措施（石莎等，2009），良好的自然條件和人工的修復(fù)使NDVI明顯增加。淺山區(qū)是深山區(qū)和平原區(qū)的過(guò)渡地帶，也是北京中心城市建設(shè)的主要生態(tài)屏障（胡艷霞等，2020），區(qū)域內(nèi)分布著潮白河水系、懷柔和密云水庫(kù)等水源地，有較強(qiáng)的水源涵養(yǎng)功能（楊文杰等，2020）。但淺山區(qū)也是北京市自然生態(tài)環(huán)境與人類活動(dòng)相互影響最強(qiáng)烈的區(qū)域，建議在此對(duì)大面積灌木和林地區(qū)域進(jìn)行治理，優(yōu)化林分組成，提高森林生態(tài)系統(tǒng)穩(wěn)定性（劉陽(yáng)等，2020）。平原區(qū)的NDVI最低，該地區(qū)為北京城市中心，人口密度高，植被主要以園林綠地和四旁樹(shù)為主。近 20年來(lái)，得益于平原造林、天然林保護(hù)、退耕還林等綠化工程以及城市小氣候的改善（賈寶全，2013；謝勝金等，2020），平原區(qū)NDVI不斷增長(zhǎng)，由此可見(jiàn)，協(xié)調(diào)好城市建設(shè)布局與生態(tài)綠化的關(guān)系是城市快速發(fā)展和生態(tài)文明建設(shè)的良好途徑。

4 結(jié)論

京津冀地區(qū)近20年的NDVI波動(dòng)變化，略有增長(zhǎng)，增速為0.025/(10 a)。其中前4年（2000—2004年）為明顯增長(zhǎng)階段，2004—2019年為波動(dòng)變化階段，NDVI相對(duì)較高的年份為2008年，明顯下降的年份為 2009年和 2014年。京津冀不同功能區(qū)的NDVI差異明顯，表現(xiàn)為西北部生態(tài)涵養(yǎng)區(qū)＞南部功能拓展區(qū)＞中部核心功能區(qū)＞東部濱海發(fā)展區(qū)。近20年，植被改善的區(qū)域占京津冀地區(qū)總面積的88.09%，退化的區(qū)域僅占11.91%；前10年（2000—2009年）植被改善的區(qū)域占京津冀地區(qū)總面積的95.13%，退化的區(qū)域面積僅占4.87%。后10年（2010—2019年）植被改善的區(qū)域略有減少，占總面積的92.88%，退化的區(qū)域增加，占7.12%。植被改善的區(qū)域主要分布在南部功能拓展區(qū)和中部核心功能區(qū)，退化的區(qū)域主要分布在西北部生態(tài)涵養(yǎng)區(qū)的張家口市和承德市，張家口市的 NDVI明顯低于其他地區(qū)，但在逐年改善。北京市的 NDVI整體呈波動(dòng)增長(zhǎng)趨勢(shì)，其中深山區(qū)NDVI最高，其次是淺山區(qū)和平原區(qū)。

整體來(lái)看，得益于平原造林、天然林保護(hù)、退耕還林等綠化工程，近 20年京津冀大部分地區(qū)的植被狀況得到了改善，但張家口市的NDVI仍然明顯低于其他地區(qū)，承德市北部植被覆蓋高但出現(xiàn)了植被退化現(xiàn)象，張家口市和承德市都位于西北部生態(tài)涵養(yǎng)區(qū)內(nèi)，對(duì)京津冀地區(qū)綠色可持續(xù)發(fā)展有重要意義，應(yīng)繼續(xù)加強(qiáng)張家口市的綠化工程，增加植被覆蓋，優(yōu)化承德北部的植被組成，提高植被穩(wěn)定性。