河北壩上地區(qū)2000-2019年植被綠度動(dòng)態(tài)及其土地利用/覆被變化歸因分析

王彥芳, 裴宏偉, 趙 超

(1.河北地質(zhì)大學(xué) 河北省農(nóng)業(yè)干旱遙感監(jiān)測國際聯(lián)合研究中心/河北省高校生態(tài)環(huán)境地質(zhì)應(yīng)用技術(shù)研發(fā)中心, 河北 石家莊050031; 2.河北建筑工程學(xué)院 市政與環(huán)境工程系, 河北 張家口 075000; 3.河北省氣象科學(xué)研究所, 河北 石家莊 050021)

植被是陸地生態(tài)系統(tǒng)的重要組成部分，是連接大氣、水圈、生物圈和土壤圈的天然紐帶，通過光合作用、蒸散和地表反照率調(diào)節(jié)全球碳循環(huán)、水循環(huán)和能量交換等[1-4]。在全球范圍內(nèi)，植被活動(dòng)的變化主要受生物地球化學(xué)效應(yīng)和社會(huì)經(jīng)濟(jì)因素的影響，前者主要包括區(qū)域氣候變化、二氧化碳的“施肥效應(yīng)”和氮沉降等，人類活動(dòng)主要包括土地利用變化和社會(huì)經(jīng)濟(jì)因素[5]。研究區(qū)域植被的連續(xù)變化及探討其影響因素，能夠反映環(huán)境演變的趨勢和原因，對(duì)區(qū)域生態(tài)評(píng)價(jià)、生態(tài)安全和經(jīng)濟(jì)可持續(xù)發(fā)展具有重要意義[6]。

自然界的綠度與地表植被覆蓋率直接相關(guān)，在中分辨率的遙感真彩色合成影像上，植被覆蓋度越高，圖像顏色越綠[4]。隨著不同時(shí)間、空間及光譜分辨率的遙感衛(wèi)星資源的持續(xù)豐富，眾多研究者借助遙感手段來實(shí)現(xiàn)對(duì)地球表面的植物種類的識(shí)別以及地表植被的覆蓋度等的研究，如基于遙感反演能靈敏地反映植被生長特征(覆蓋度、葉面積指數(shù)、生物量等)[7-8]、生物物理生化特性(葉含水量、葉綠素含量等)[9-10]和生態(tài)系統(tǒng)參數(shù)(高度、長勢等)等[11]。其中，植被指數(shù)(vegetation index)選用對(duì)綠色植物強(qiáng)吸收的可見光紅光波段和對(duì)綠色植物高反射的近紅外波段構(gòu)建，成為衡量地表植被綠度的重要指標(biāo)[4]。在眾多植被指數(shù)中，歸一化植被指數(shù)(normalized difference vegetation index, NDVI)是植被生長狀態(tài)和植被綠度的最佳指示因子之一，在全球范圍內(nèi)廣泛應(yīng)用[5,12-13]。

已有的研究在區(qū)域及大陸尺度植被綠度變化規(guī)律、植被對(duì)氣候因子的響應(yīng)規(guī)律，以及氣候變化與人類活動(dòng)對(duì)植被變化的貢獻(xiàn)展開廣泛的研究[13-17]。研究表明，20世紀(jì)80年代以來，全球植被綠度不斷改善，尤其是北半球中高緯地區(qū)[13]，中國西北、華北、青藏高原和黃土高原等地區(qū)植被綠度也在不斷增長，不同區(qū)域植被變化受氣候變化和人類活動(dòng)影響程度不同[2,15-17]。土地是人類活動(dòng)的空間載體，改變土地利用方式是人類活動(dòng)影響植被綠度的重要手段，將土地利用數(shù)據(jù)與遙感NDVI相結(jié)合可以評(píng)估不同土地利用方式下區(qū)域植被生態(tài)狀況。例如，Chen等[14]研究發(fā)現(xiàn)中國和印度的陸地植被面積對(duì)全球綠化增量的貢獻(xiàn)約1/3，其中，中國占1/4，中國的綠化主要來自森林(42%)和農(nóng)田(32%)，而印度的綠化主要來自農(nóng)田(82%)，森林的貢獻(xiàn)較小(4.4%)。景振華等[18]分析了山西六大煤田NDVI對(duì)土地覆蓋類型的響應(yīng)特征；道日娜[19]研究認(rèn)為1995—2010年，三北地區(qū)NDVI減小的主要原因是大量的草地轉(zhuǎn)為耕地，NDVI增加與大量未利用地轉(zhuǎn)為耕地有關(guān)。但目前的研究，缺少結(jié)合土地利用方式變化的植被綠度變化的來源定量分析，明確區(qū)域植被綠度變化所發(fā)生的主要土地類型。因此，本研究以京津冀生態(tài)環(huán)境支撐區(qū)中重要的區(qū)域之一——河北壩上地區(qū)為例，采用2000—2019年的MODIS MOD13Q1 NDVI數(shù)據(jù)，結(jié)合2000，2010和2018年土地利用數(shù)據(jù)，分析壩上近20 a植被綠度變化及土地利用/覆被變化歸因，旨在為地區(qū)生態(tài)建設(shè)和京津冀生態(tài)環(huán)境支撐區(qū)的建設(shè)提供一定參考，為土地合理利用、生態(tài)環(huán)境保護(hù)和地區(qū)可持續(xù)發(fā)展提供科學(xué)依據(jù)和決策支持。

1 研究區(qū)概況

壩上地區(qū)指京津冀地區(qū)北部，系內(nèi)蒙古高原的南緣，其中本研究涉及的河北壩上地區(qū)包括張家口市的4縣兩區(qū)(張北縣、康?？h、尚義縣、沽源縣、察北管理區(qū)、塞北管理區(qū))和承德市的豐寧滿族自治縣和圍場滿族蒙古族自治縣一部分。該區(qū)域總面積近1.90×104km2，占河北省總面積約10%，海拔多在1 200～1 800 m，地形為丘陵和平原，呈東南高、西北低格局(圖1)；河北壩上人口約112萬人，占河北省總?cè)丝诘?.7%[20]。氣候條件方面，壩上系大陸性季風(fēng)氣候，自然條件惡劣，表現(xiàn)在寒冷、多風(fēng)、干旱等，壩上年均氣溫1～2 ℃，無霜期90～120 d，年降水量400 mm左右[20]。壩上以閃電河為界，東部屬于壩上外流區(qū)，降雨較多，地表覆蓋以森林草原和草甸草原為主；西部屬于封閉的內(nèi)流區(qū)，氣候干旱，以干草原為主，歷史上曾是一個(gè)原生態(tài)的自然生態(tài)系統(tǒng)，曾經(jīng)河網(wǎng)密布，水淖豐富[21]。由于受氣候變化和人類開灘種地和過度放牧等活動(dòng)的影響，草場逐步退化，湖淖不斷萎縮，地下水位下降明顯。

圖1 研究區(qū)范圍及地形

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù)來源

本研究所用的主要數(shù)據(jù)包括遙感NDVI數(shù)據(jù)和土地利用數(shù)據(jù)。NDVI數(shù)據(jù)為MODIS (moderate-resolution Imaging Spectroradiometer， MODIS)MOD13Q1數(shù)據(jù)，時(shí)間序列為2000—2019年，空間分辨率為250 m。為了提高數(shù)據(jù)精度，該數(shù)據(jù)采用最大值合成方法(MVC)減少了云和氣溶膠等的影響，合成后數(shù)據(jù)的時(shí)間分辨率為16 d。土地利用數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心(http:∥www.resdc.cn/)，包括2000，2010和2018年的3期壩上地區(qū)土地利用數(shù)據(jù)，該數(shù)據(jù)來源于國家尺度1∶10萬比例尺多時(shí)期土地利用/土地覆蓋遙感監(jiān)測數(shù)據(jù)庫(CNLUCC)，以美國陸地衛(wèi)星Landsat遙感影像數(shù)據(jù)作為主信息源，通過人工目視解譯獲得，數(shù)據(jù)空間分辨率為30 m[22]。土地利用類型根據(jù)土地資源及利用屬性包括耕地、林地、草地、水域、居民地和未利用土地6個(gè)1級(jí)類型，分類結(jié)果的綜合評(píng)價(jià)精度在94.3%以上[22]。

2.2 分析方法

2.2.1 植被綠度 NDVI可以很好地反映區(qū)域植被綠度。根據(jù)植被綠度季節(jié)性變化差異，本研究同時(shí)采生長季平均NDVI和年最大NDVI分別評(píng)價(jià)區(qū)域年平均和最大綠度。由于受氣候條件的影響，該區(qū)域的植被生長主要集中在4—10月。本研究以該時(shí)期內(nèi)平均NDVI表示年生長季平均植被綠度。首先將MOD13Q1 NDVI每年16 d間隔的23個(gè)波段合成一個(gè)文件，生長季平均NDVI即波段范圍為7～20的平均值。2000—2019年逐像元生長季平均NDVI和年最大NDVI計(jì)算公式為：

(1)

NDVImax,k=max(NDVIk,i)

(2)

2.2.2 線性傾向估計(jì) 本研究利用線性傾向定量的估算各柵格植被綠度變化趨勢。對(duì)于樣本量為n的變量xi來說，線性傾向估計(jì)通過建立變量序列xi和時(shí)間ti之間的一元線性回歸關(guān)系[4]，即

(3)

式中：回歸系數(shù)b表示變量x的趨勢傾向，其絕對(duì)值大小上升或者下降的速率，符號(hào)表示變化方向，b可以通過最小二乘法擬合得到：

(4)

(5)

(6)

因此，每種土地利用類型對(duì)于區(qū)域植被綠度變化的貢獻(xiàn)量為Rl

(7)

年最大綠度變化及其定量分析與生長季平均值相同。

3 結(jié)果和分析

3.1 河北壩上地區(qū)植被綠度時(shí)空變化特征

3.1.1 壩上地區(qū)植被綠度的空間分布特征 不同土地利用狀況下，地區(qū)的植被綠度存在明顯差異。2000—2019年壩上生長季多年平均植被綠度和年最大NDVI的多年平均態(tài)如圖2a,2b所示。由圖2可知，壩上生長季平均NDVI為0.39，年最大值均值為0.63，空間上呈現(xiàn)壩東高于壩西。植被最差的區(qū)域主要分布在干涸的湖淖區(qū)，其次是康?？h北部，尚義縣南部造林初期區(qū)域植被。NDVI變異系數(shù)可以體現(xiàn)植被穩(wěn)定性，2000—2019年NDVI變異系數(shù)結(jié)果顯示，壩東植被較壩西穩(wěn)定，壩西地區(qū)變異系數(shù)浮動(dòng)在15%之內(nèi)，植被綠度變異最大為壩西湖淖區(qū)，其次是尚義南部和北部地區(qū)(圖2c)。結(jié)合土地利用狀況，除水域區(qū)，植被變異較大的區(qū)域主要是人類活動(dòng)影響較大的區(qū)域，如耕地植被和建設(shè)用地的植被，其次是草地和未利用地，林地最為穩(wěn)定。統(tǒng)計(jì)2000，2010和2018年不同土地類型的平均NDVI植被生長曲線(圖2d)，結(jié)果表明，林地植被綠度最高，其次草地、耕地、未利用地，除水域區(qū)域外，建設(shè)用地區(qū)域植被綠度最低。從5月中旬起，人工植被即耕地和建設(shè)用地NDVI增長速率明顯高于林地、草地等自然植被。

圖2 壩上地區(qū)植被綠度空間分布特征

3.1.2 近20 a壩上地區(qū)植被綠度變化趨勢 2000—2019年間，河北壩上地區(qū)植被綠度(生長季均值或年最大值)呈現(xiàn)整體增加的趨勢，其中，年生長季均值增速為0.044/10 a，年最大值平均增速更大，達(dá)到0.063/10 a，植被綠度增加的區(qū)域占總面積的90%以上(圖3)。增幅較大的區(qū)域主要分布在尚義縣中南部、康保縣和張北縣西部、沽源縣大部地區(qū)以及豐寧和圍場的中北部地區(qū)，這些區(qū)域增加趨勢的顯著性水平也較高?？当？h北部零星區(qū)域?yàn)橄陆第厔?，但并不顯著(圖4)。

圖3 河北壩上地區(qū)2000-2019年生長季平均NDVI和年最大NDVI的變化趨勢

圖4 河北壩上地區(qū)2000-2019年生長季平均NDVI和年最大NDVI變化趨勢的顯著性分析

生長季平均綠度呈顯著(p<0.05)增加趨勢的比例為83%，不顯著上升趨勢的為15%，不顯著下降的比例為1.6%，顯著下降的區(qū)域面積不足0.5%，主要為建設(shè)用地?cái)U(kuò)張。年最大植被綠度變化趨勢和生長季平均狀況大致相同，但程度不同，其中，綠度顯著上升的面積占比為60%，不顯著上升趨勢的面積比例次之為34.8%，不顯著下降區(qū)域的占比為4.8%，顯著下降趨勢不足0.4%。因此，無論是生長季平均綠度還是最大綠度，近20 a均呈現(xiàn)上升趨勢，顯著上升區(qū)占比在60%～83%。

3.2 河北壩上地區(qū)土地利用/覆被變化特征

2000—2018年的河北壩上地區(qū)土地利用分布和轉(zhuǎn)移情況如圖5和表1所示，自“游牧區(qū)”轉(zhuǎn)變?yōu)椤鞍朕r(nóng)半牧區(qū)”以來，壩上地區(qū)耕地面積不斷增加，目前耕地面積已達(dá)到區(qū)域土地面積的1/2以上，成為河北壩上地區(qū)主要的土地利用類型，其次是草地和林地，3種類型占到區(qū)域總面積的91%。但耕地主要分布在中西部地區(qū)，林草地主要分布在壩東地區(qū)。2000—2018年壩上耕地變幅不大(圖5)，耕地面積轉(zhuǎn)移主要表現(xiàn)在轉(zhuǎn)為林地、草地和建設(shè)用地，變化的原因主要是由于農(nóng)田開墾、退耕還林退耕還草、耕地撂荒和城鎮(zhèn)擴(kuò)張(表1)等。與草地轉(zhuǎn)換密切的土地類型主要有耕地、林地和未利用地，但總體來說，草地面積近20 a減少了10%。同樣減少的還有水域面積，但壩上河流和湖泊多屬于季節(jié)性水域，由于數(shù)據(jù)來源的季節(jié)性差異導(dǎo)致土地利用數(shù)據(jù)的水域面積提取具有較大不確定性，但從實(shí)地調(diào)查和多源遙感數(shù)據(jù)看，近20 a河北壩上水域面積大幅萎縮。壩上林地以增加為主，主要來源于耕地和草地的退耕還林、植樹造林活動(dòng)。

圖5 河北壩上地區(qū)2000，2010和2018年土地利用狀況

表1 河北壩上地區(qū)2000-2018年土地利用轉(zhuǎn)移矩陣 hm2

建設(shè)用地增幅最大，耕地和草地為主要的轉(zhuǎn)入類型，轉(zhuǎn)出類型主要為耕地，主要是對(duì)壩上部分空心化嚴(yán)重村莊，實(shí)施整體拆遷和易地扶貧搬遷，對(duì)“空心村”拆遷后的土地進(jìn)行復(fù)墾，從轉(zhuǎn)出類型看大部分轉(zhuǎn)移為耕地，少部分轉(zhuǎn)為草地。與未利用地轉(zhuǎn)移較多的土地類型為耕地和草地，主要表現(xiàn)為農(nóng)業(yè)開墾和草地沙化，另外部分湖淖干涸形成鹽堿地?？偟膩碚f，2000年以來河北壩上地區(qū)土地利用變化以政策性引導(dǎo)為驅(qū)動(dòng)，隨著國家“三北”防護(hù)林建設(shè)工程、退耕還林退耕還草等生態(tài)項(xiàng)目的落地實(shí)施和不斷推進(jìn)，包括草地、林地等生態(tài)用地小幅增加，生態(tài)環(huán)境有所改善。例如尚義縣南部，自2002年實(shí)施退耕還林以來，甲石河鄉(xiāng)和下馬圈鄉(xiāng)共退耕地3 000 hm2以上[23]，2017年林草覆蓋率達(dá)70%以上。

但同時(shí)，近十幾年蔬菜的規(guī)?；⒓s化種植加大對(duì)區(qū)域水資源的消耗，耕地面積雖然變化不大，但是耕作強(qiáng)度增加，從雨養(yǎng)農(nóng)業(yè)轉(zhuǎn)為灌溉農(nóng)業(yè)，加大區(qū)域生態(tài)安全威脅。

3.3 區(qū)域植被綠度動(dòng)態(tài)的土地利用/覆被變化歸因

考慮土地利用變化，采用公式(5)—(7)分析不同時(shí)期各土地利用/覆被的面積變化和強(qiáng)度變化對(duì)區(qū)域植被綠度的貢獻(xiàn)(表2)。從表2可以看出，除了2000—2010年未利用地的貢獻(xiàn)為負(fù)之外，其他時(shí)段和類型的貢獻(xiàn)均為正值。2000—2018年，壩上地區(qū)生長季平均NDVI的變化量中，1/2以上發(fā)生在耕地，包括耕地開墾、撂荒、退耕還林還草等措施導(dǎo)致的耕地面積的變化，以及種植結(jié)構(gòu)變化、旱地水澆地變化和氣候變化引起的耕地植被綠度的變化；其次，林地和草地的貢獻(xiàn)率分別為28.62%和15.07%，水域、建設(shè)用地和未利用地的總貢獻(xiàn)率在7%左右。在不同發(fā)展時(shí)期，各土地利用方式的貢獻(xiàn)有一定差異，2000—2010年耕地的影響明顯大于其他土地類型，耕地總面積和耕地綠度都有明顯增加，對(duì)區(qū)域生長季平均NDVI的貢獻(xiàn)量為0.034，貢獻(xiàn)率達(dá)到65.74%，耕地和林地面積和綠度變化的影響占區(qū)域植被變化的90%以上。2010年后，雖然耕地植被綠度仍呈增加趨勢，但由于面積出現(xiàn)下降，因此，耕地的影響逐漸下降。相反，在退耕還林退耕還草的政策下，林草面積不斷上升，耕地、林地和草地影響趨于均等(表2)。

表2 河北壩上地區(qū)不同時(shí)期各土地利用/覆被變化對(duì)區(qū)域植被綠度變化的貢獻(xiàn)率 %

耕地對(duì)年最大NDVI的影響較生長季平均NDVI更大，結(jié)果表明，2000—2018年的年最大NDVI的變化中57.22%來自于耕地，高于生長季平均值的50.52%，其次，仍然是林地和草地，貢獻(xiàn)率分別為21.73%和14.41%，較對(duì)生長季平均NDVI的影響略小。不同階段，耕地的貢獻(xiàn)率差異明顯，2000—2010年耕地的綠度貢獻(xiàn)率達(dá)67.29%，2010—2018年降至43.74%，林地的影響較為穩(wěn)定，草地的貢獻(xiàn)率上升明顯。水域、建設(shè)用地和未利用地由于面積占比較少，對(duì)區(qū)域植被綠度總趨勢影響不大，但呈現(xiàn)明顯的增加趨勢，從前10 a的2.18%～3.18%增加到2010—2018年的9.62%～10.44%。

4 討 論

作為首都和京津冀重要生態(tài)環(huán)境支撐區(qū)和水源涵養(yǎng)區(qū)，河北壩上地區(qū)植被狀況對(duì)于其生態(tài)功能的發(fā)揮起著重要作用，因此對(duì)壩上地區(qū)的植被綠度的分析具有重要意義。本研究以MODIS MOD13Q1 NDVI遙感數(shù)據(jù)為反映植被綠度的數(shù)據(jù)源，結(jié)合Landsat土地利用數(shù)據(jù)，分析2000—2019年壩上植被綠度的年際變化趨勢，并定量分析其變化的土地利用/覆被變化歸因。植被綠度空間特征與時(shí)間序列變化方面，本研究的結(jié)果與前人研究結(jié)果類似[24-28]，壩上植被綠度(覆蓋)和穩(wěn)定性均表現(xiàn)為從壩西到壩東依次遞增，2000—2019年壩上地區(qū)的植被綠度呈現(xiàn)明顯的增加趨勢，60%以上的區(qū)域呈顯著增加趨勢[27]，土地沙化程度減輕，生態(tài)狀況好轉(zhuǎn)，相應(yīng)地，區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值呈增加態(tài)勢，2000—2015年生態(tài)系統(tǒng)服務(wù)總價(jià)值共增加334.83億元[24]。

影響因素方面，通常研究主要分析氣候變化和人類活動(dòng)的影響，尤其是氣候條件的波動(dòng)。例如陳輝等[28]研究表明1987—2000年冀北地區(qū)土地利用類型變化、錯(cuò)季蔬菜生產(chǎn)等人為因素對(duì)NDVI的變化作用不明顯，降水、特別是生長季降水分布特征是影響NDVI變化的重要因素；類似地，2000—2013年壩上地區(qū)年均NDVI也與年降水的相關(guān)性較大[29]。壩上地區(qū)作為典型的水分限制區(qū)，植被對(duì)降水的變化較為敏感。

本研究聚焦區(qū)域植被綠度發(fā)生的空間，即通過不同土地利用類型的綠度和面積變化的分解，定量得到每一種土地利用/覆被類型的變化(植被綠度和面積)對(duì)于整個(gè)區(qū)域植被綠度的影響。其中，某種土地利用方式的植被綠度變化，既有氣候因素的影響也有人類活動(dòng)的影響。研究發(fā)現(xiàn)，2000—2019年壩上區(qū)域植被綠度變化的1/2以上的貢獻(xiàn)發(fā)生在耕地，其次是林地和草地，貢獻(xiàn)率分別為21.73%～28.62%，14.41%～15.07%。結(jié)合土地利用強(qiáng)度與土地利用類型轉(zhuǎn)移特征，可以得出耕地的貢獻(xiàn)量中，除了氣候原因之外人類活動(dòng)影響最大的是農(nóng)業(yè)種植結(jié)構(gòu)變化，如傳統(tǒng)雨養(yǎng)耕地到水澆地的轉(zhuǎn)變導(dǎo)致利用強(qiáng)度的轉(zhuǎn)變，另一方面在于耕地開墾、退耕還林還草工程、易地搬遷舊宅基地騰退復(fù)墾項(xiàng)目等引起的耕地面積和使用強(qiáng)度的變化。林地和草地植被綠度的變化類似，除了受降雨、氣溫等自然要素因素的影響之外，植樹造林工程、退耕還林還草、人工林撫育等項(xiàng)目都是影響林地和草地面積和綠度變化的關(guān)鍵因素。2010年后，耕地的影響逐漸下降，耕地、林地和草地影響趨于均等。為建設(shè)好首都水源涵養(yǎng)功能區(qū)和生態(tài)環(huán)境支撐區(qū)，壩上地區(qū)仍需加快生態(tài)環(huán)境恢復(fù)，加大恢復(fù)力度，進(jìn)一步增加林草生態(tài)用地面積，對(duì)“空心村”治理后的原宅基地及耕地復(fù)墾后實(shí)施還林還草，避免流轉(zhuǎn)作為水澆地。

本研究結(jié)果定量的指出區(qū)域植被綠度發(fā)生的主要土地利用/覆被類型，掌握了區(qū)域綠度變化來源，為地區(qū)生態(tài)建設(shè)和京津冀生態(tài)環(huán)境支撐區(qū)的建設(shè)提供一定參考，為京津冀地區(qū)可持續(xù)發(fā)展、土地合理利用及生態(tài)環(huán)境保護(hù)提供科學(xué)依據(jù)和決策支持。但該研究仍具有一定的局限性和不確定性。一方面未能對(duì)歸因進(jìn)行更細(xì)致的分析，如2000—2018年耕地引起的區(qū)域植被綠度50.51%～57.22%的變化中，來自于耕地面積的變動(dòng)和耕地綠度變化所占比例各為多少，耕地與其他土地利用類型之間的轉(zhuǎn)移各對(duì)區(qū)域植被綠度產(chǎn)生多大的影響，還缺少深層次的分析。以后的研究可以結(jié)合具體的研究內(nèi)容和應(yīng)用領(lǐng)域，可以針對(duì)某地類，展開更加深入的分析。

5 結(jié) 論

(1) 2000—2019年河北壩上地區(qū)植被生長季平均NDVI為0.39，年最大值均值為0.63，空間上呈現(xiàn)壩東高于壩西，壩東植被穩(wěn)定性較壩西高，林地植被綠度最高，其次草地、耕地、未利用地，除水域區(qū)域外，建設(shè)用地區(qū)域植被綠度最低。

(2) 2000—2019年河北壩上地區(qū)植被綠度呈現(xiàn)增加趨勢，其中，生長季平均NDVI增速為0.044/10 a，年最大值增速更大，達(dá)到0.063/10 a，植被綠度增加的區(qū)域占總面積的95%以上，其中顯著增加趨勢占60%～83%。

(3) 土地利用/覆被變化歸因分析表明，2000—2018年，壩上地區(qū)植被綠度變化量中，耕地貢獻(xiàn)了50.51%～57.22%，林地和草地分別貢獻(xiàn)21.73%～28.62%和14.41%～15.07%，水域、建設(shè)用地和未利用地的總貢獻(xiàn)率在6%左右。其中，2000—2010年耕地植被變化對(duì)區(qū)域植被綠度變化的貢獻(xiàn)率在65%以上，林地和草地植被貢獻(xiàn)不足1/3；2010—2018年，耕地的貢獻(xiàn)在降低，其他土地利用類型的貢獻(xiàn)呈現(xiàn)逐步上升趨勢，耕地、林地和草地的貢獻(xiàn)趨于均等。