2000-2015年青海省植被覆蓋的時空變化特征

代子俊，趙 霞，李冠穩(wěn)，王興春，龐龍輝

(青海師范大學(xué) 地理科學(xué)學(xué)院，中美合作(青海師范大學(xué)&普渡大學(xué))自然資源與環(huán)境建模實驗室，青海省自然地理與環(huán)境過程重點實驗室，青海 西寧 810008)

植被是連接各個生態(tài)系統(tǒng)的“紐帶”，對生態(tài)系統(tǒng)的穩(wěn)定起著重要作用[1]。歸一化植被指數(shù)(normalized difference vegetation index,NDVI)可以反映植被生長狀況及植被覆蓋變化特征，被用于研究小區(qū)域乃至全球的植被、生態(tài)環(huán)境變化，在植被指數(shù)中應(yīng)用最為廣泛[2-3]。

在氣候暖干化和人類活動綜合影響下，青海省植被覆蓋發(fā)生了重大變化。已有研究表明，青海省植被覆蓋整體呈增長趨勢[2]，局部區(qū)域有退化趨勢[4-5]，該區(qū)域植被覆蓋變化對氣候變化較為敏感[6-7]。劉櫟杉等[8]利用1 km MODIS/NDVI數(shù)據(jù)，采用線性趨勢分析、標(biāo)準(zhǔn)差等方法，研究了青海省2000-2009年“退耕還林還草”實施10年來植被覆蓋的時空變化特征，結(jié)果表明青海省植被覆蓋變化以輕微改善較多，明顯改善的地區(qū)主要集中在柴達木盆地東南邊緣、青海湖盆地、茶卡-共和盆地、河湟谷地及黃南低地；王莉雯等[9]基于1 km MODIS/NDVI數(shù)據(jù)，對青海省2001-2006年的地表植被覆蓋空間分布和時間序列變化進行了定量分析，發(fā)現(xiàn)5年中青海省植被覆蓋有所改善，植被覆蓋率增加了0.63%。然而，目前研究青海省植被NDVI變化所用影像的分辨率普遍較低[2,4-9]，采用的方法也多為線性趨勢統(tǒng)計方法，其規(guī)避誤差的能力較差,并且針對青海省植被NDVI未來變化趨勢的研究鮮見報道。

MODIS數(shù)據(jù)產(chǎn)品時間跨度長、時間分辨率高、覆蓋范圍廣且數(shù)據(jù)免費易獲取，因此被廣泛用于研究長時間跨度的植被覆蓋變化[10-11]；MODIS13Q1數(shù)據(jù)產(chǎn)品能夠詳盡地反映區(qū)域尺度上植被覆蓋的空間差異性，同時對植被覆蓋狀況敏感區(qū)域的生態(tài)環(huán)境監(jiān)測也具有明顯優(yōu)勢[12-14]。因此，本研究基于250 m分辨率的MODIS13Q1數(shù)據(jù)，通過MRT軟件拼接、投影轉(zhuǎn)換及最大值合成，得到2000-2015年青海省逐年植被生長季節(jié)NDVI值，采用規(guī)避誤差能力較強的Sen+Mann-Kendall趨勢變化分析方法，以像元為基本單元研究植被NDVI的空間變化特征，并通過Hurst指數(shù)分析，研究青海省未來植被NDVI變化趨勢；同時輔以青海省氣溫和降水?dāng)?shù)據(jù)，揭示植被NDVI變化對氣候變化的響應(yīng)規(guī)律，以期為青海省生態(tài)建設(shè)提供一定參考。

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

1.1 研究區(qū)概況

青海省位于中國西北部內(nèi)陸腹地，地理位置為89°35′～103°04′E和31°39′～39°19′N，屬于高原大陸性氣候，氣溫普遍較低，區(qū)域間氣溫差異顯著且垂直變化明顯。根據(jù)青海省地形特征、氣候差異和人類活動等因素，將研究區(qū)劃分為柴達木盆地、環(huán)湖和祁連山區(qū)、東部農(nóng)業(yè)區(qū)及青南牧區(qū)4個子區(qū)。東部農(nóng)業(yè)區(qū)年平均氣溫2～9 ℃，年降雨量250～550 mm，降雨較多且主要集中于7-9月；柴達木盆地年平均氣溫2～5 ℃，年降雨量近200 mm，氣候干燥，年日照時長達3 000 h以上；東北部的祁連山區(qū)及青南高原年均溫度較低，年降雨量一般為100～500 mm。受地形和氣候條件限制，青海省植被類型多以高寒灌叢、高寒草甸及高寒草原為主，荒漠和山地草原次之，森林植被則比較少。青海省北部降水量由東向西逐漸減少，植被類型由東北部溫性草原植被向柴達木盆地東部草原荒漠、中部荒漠、西北部風(fēng)蝕殘丘、流沙過渡；南部的青南高原處在西北-東南近似東西走向的東昆侖山和唐古拉山之間，植被類型由東南部的寒溫性針葉林、高寒灌叢向中部瑪沁-玉樹高寒灌叢、高寒草甸及西北黃河-長江上游高寒草甸過渡[15]。在1∶400萬中國植被圖基礎(chǔ)上，通過青海省行政區(qū)邊界裁剪，得到青海省植被類型圖。青海省子區(qū)劃分及植被類型分布見圖1。

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

遙感數(shù)據(jù)為2000-2015年的MODIS13Q1數(shù)據(jù)集，來源于美國航天航空(NASA)數(shù)據(jù)中心，網(wǎng)址為https://wist.echo.nasa.gov/api。該數(shù)據(jù)集為16 d合成的L3級數(shù)據(jù)產(chǎn)品，空間分辨率為250 m×250 m。青海省覆蓋了h25v5、h26v5兩幅MOD13Q1數(shù)據(jù)。借鑒趙健赟等[2]對青海省年內(nèi)植被生長變化的研究結(jié)果，青海省植被生長季節(jié)一般在5-9月，故本研究選取5-9月的NDVI數(shù)據(jù)，每月2景，生長期內(nèi)共160景。首先利用MRT(MODIS ReProjection Tools)軟件進行影像拼接和投影轉(zhuǎn)換，然后采用最大值合成法(maximum value composite,MVC)獲得逐年植被生長季節(jié)5-9月最大合成NDVI數(shù)據(jù)。氣象數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn)。

1.3 研究方法

1.3.1 NDVI年際變化 利用MVC法[16-17]得到青海省2000-2015年逐年每月的NDVI數(shù)據(jù)，公式如下:

(1)

研究區(qū)所有像元2000-2015年逐年生長季節(jié)平均NDVI值計算公式為：

(2)

1.3.2 Sen+Mann-Kendall趨勢變化分析 將能反映趨勢變化的Theil-Sen Median方法和檢驗趨勢變化顯著性的Mann-Kendall方法結(jié)合[18-20]，與常見的植被年際變化一元線性回歸方法相比，該方法對數(shù)據(jù)誤差的規(guī)避能力較強，可以減少異常值的影響，使分析檢驗結(jié)果更準(zhǔn)確可靠[21]。本研究應(yīng)用MATLAB 2016生成Theil-Sen Median和Mann-Kendall函數(shù)。Theil-Sen Median趨勢分析計算公式如下:

2000≤i

(3)

Mann-Kendall檢驗公式如下：

(4)

其中,

方差V(S)=n(n-1)(2n+5)/18。統(tǒng)計量Z通過下式計算：

1.3.3 Hurst指數(shù)分析 R/S分析是一種基于長程相關(guān)思想的時間序列分析方法[22]，由H.E.Hurst于1965年最先提出，近期被用于研究植被NDVI變化的可持續(xù)性判斷[23-25]。為使結(jié)果準(zhǔn)確可靠，應(yīng)用IDL8.5生成Hurst函數(shù)，其基本步驟如下：

給定時間序列{NDVI(t)}，t=1,2,3, …，n。對于任意正整數(shù)u，定義該時間均值序列：

(5)

1≤t≤u。

(6)

(7)

(8)

若存在R/S∝uH，則表示時間序列{NDVI(t)}存在Hurst現(xiàn)象，Hurst指數(shù)用H值表示。根據(jù)H值可判斷植被NDVI變化是完全隨機還是存在持續(xù)性。當(dāng)H=0.5時，表示NDVI時間序列自相關(guān)性為0，即時間序列前后變化是完全隨機的；當(dāng)0.5

Ⅰ.柴達木盆地；Ⅱ.環(huán)湖和祁連山區(qū)；Ⅲ.東部農(nóng)業(yè)區(qū)；Ⅳ.青南牧區(qū)Ⅰ.Qaidam basin;Ⅱ.Central Qinghai lake and Qilian mountains;Ⅲ.Eastern agricultural zone;Ⅳ.Qingnan pastoral area圖1 青海省子區(qū)劃分及植被類型分布Fig.1 Subdivision of Qinghai province and the distribution of vegetation types

2 結(jié)果與分析

2.1 青海省平均NDVI的空間分布特征

利用2000-2015年的逐年生長季節(jié)NDVI數(shù)據(jù)，取平均值后得到青海省柴達木盆地、環(huán)湖和祁連山區(qū)、東部農(nóng)業(yè)區(qū)及青南牧區(qū)4個子區(qū)16年間平均NDVI空間分布圖(圖2)。

Ⅰ.柴達木盆地；Ⅱ.環(huán)湖和祁連山區(qū)；Ⅲ.東部農(nóng)業(yè)區(qū)；Ⅳ.青南牧區(qū)Ⅰ.Qaidam basin; Ⅱ.Central Qinghai lake and Qilian mountains; Ⅲ.Eastern agricultural zone; Ⅳ.Qingnan pastoral area圖2 2000-2015年青海省平均NDVI的空間分布Fig.2 Spatial distribution of mean NDVI in Qinghai from 2000 to 2015

由圖2可以看出，青海省植被NDVI空間分布特征整體呈現(xiàn)從西北向東南逐漸增大的趨勢，植被NDVI值較好(0.3以上)的區(qū)域占青海省面積的53.7%。青海省北部降水從東向西逐漸減少，干旱現(xiàn)象逐漸加重，植被類型由東北部溫性草原植被向柴達木盆地東部草原荒漠、柴達木盆地中部荒漠、柴達木盆地西北部風(fēng)蝕殘丘、流沙過渡，因此西部柴達木盆地具有較低的NDVI值；而東部農(nóng)業(yè)區(qū)海拔相對較低，雨雪豐富，植被生長旺盛，同時也是青海省主要農(nóng)作物種植區(qū)，NDVI值也較高；青海省南部的青南牧區(qū)，自東南向西北，地形由山地河谷、峽谷向高原過渡，植被類型由東南部的寒溫性針葉林、高寒灌叢向中部瑪沁-玉樹高寒灌叢、高寒草甸及西北黃河-長江上游高寒草甸過渡[15]，因此NDVI也由東南向西北逐漸降低。

2.2 青海省NDVI的年際變化

2.2.1 青海全省及各子區(qū)的變化 根據(jù)2000-2015年逐年生長季節(jié)NDVI值，基于整幅影像取平均值(式(2))，制作年際NDVI變化圖(圖3-a)。由圖3-a可以看出，16年來青海省NDVI整體呈增長趨勢，增速為每10年1.5%，并且通過了逐年P(guān)=0.05水平的顯著性檢驗(表1)；NDVI值在0.37～0.42波動，其中2003-2005年、2008-2010年NDVI快速上升，2008年NDVI值最低，處于波谷。

為了更準(zhǔn)確地分析青海省16年來植被NDVI的變化規(guī)律，通過分區(qū)方法，統(tǒng)計各子區(qū)NDVI的變化差異。從圖3-b-e可以看出，東部農(nóng)業(yè)區(qū)、青南牧區(qū)、環(huán)湖和祁連山區(qū)以及柴達木盆地NDVI均呈增加趨勢，各子區(qū)NDVI變化規(guī)律與青海省整體NDVI變化規(guī)律一致，不同之處在于各子區(qū)NDVI變化速率及相對變化幅度稍有差異。其中東部農(nóng)業(yè)區(qū)增長最快，增速為每10年2.3%(表1)，且年際差異顯著，該區(qū)2000-2002年持續(xù)增長，但2008-2010年相對變化幅度低于青海省整體；環(huán)湖和祁連山區(qū)及青南牧區(qū)增速分別為每10年2.2%和1.3%，但均未通過顯著性檢驗；柴達木盆地增速為每10年1.4%，并且通過了顯著性檢驗(P<0.01)。柴達木盆地NDVI變化明顯，主要是由于近年來積極實施植樹造林，改變了原有的荒漠生態(tài)系統(tǒng)格局，植被覆蓋明顯改善；青南牧區(qū)NDVI整體呈增加趨勢且相對變化幅度較大，這與國家實施的三江源區(qū)生態(tài)保護與建設(shè)工程密切相關(guān)，同時也與局部地區(qū)過度放牧、亂挖蟲草及氣候變化有關(guān)。

2.2.2 不同植被類型的變化 青海省植被類型以草原與稀樹灌木草原、草甸與草本沼澤、灌叢與萌生矮林為主，其次是荒漠、旱地、針葉林。不同植被類型的NDVI變化規(guī)律見圖4。

圖4 2000-2015年青海省各植被類型NDVI的年際變化Fig.4 Inter-annual variability of NDVI in different vegetation types in Qinghai from 2000 to 2015

由圖4可知，所有植被類型NDVI均呈增加趨勢，但增加幅度略微不同。由表2可以看出，旱地、草原與稀樹灌木草原的NDVI增速最快，分別為每10年2.7%和2.3%，其中草原與稀樹灌水草原通過了顯著性檢驗(P<0.05)；針葉林、荒漠的NDVI增速分別為每10年1.5%和 1.1%，但荒漠的年際變化幅度較小，這是由于荒漠地區(qū)人類活動干擾較少，因此氣候變化是該區(qū)植被改善的主要原因；草甸與草本沼澤、灌叢與萌生矮林的NDVI增速均為每10年1.0%，呈改善趨勢，且兩者的相對變化幅度較大，這與青海省近幾年實施退耕還林還草政策有關(guān)。

表2 2000-2015年青海省各植被類型NDVI的變化率Table 2 Changes of NDVI in different vegetation types in Qinghai from 2000 to 2015

2.3 青海省NDVI的變化趨勢

對16年來青海省植被生長季節(jié)NDVI進行Theil-Sen Median趨勢變化分析，采用Mann-Kendall方法檢驗趨勢變化的顯著性，得到基于像元尺度的2000-2015年青海省NDVI變化趨勢空間分布特征。再根據(jù)前人研究[25]，將Theil-Sen Median植被NDVI趨勢分析的分級結(jié)果與Mann-Kendall趨勢顯著性檢驗的分級結(jié)果在Arcgis10.3中疊加，可將NDVI變化結(jié)果劃分為5種類型(表3)。由表3可以看出，16年間青海省植被覆蓋改善區(qū)域(60.51%)明顯大于退化區(qū)域(17.87%)，其中植被覆蓋明顯改善區(qū)域占青海省總面積的21.26%，輕微改善區(qū)域比重最大，占39.25%，穩(wěn)定不變的區(qū)域占21.62%，輕微退化區(qū)域占15.75%，嚴重退化區(qū)域占2.12%。

表3 2000-2015年青海省植被NDVI變化趨勢統(tǒng)計分析Table 3 Statistical analysis of vegetation NDVI variation in Qinghai from 2000 to 2015

將植被變化趨勢和顯著性檢驗疊加，結(jié)果(圖5)顯示，16年間青海省植被改善的區(qū)域分布較集中，而植被輕度退化區(qū)域分布較零散，且植被改善區(qū)域的面積顯著大于退化區(qū)域。祁連山西部及柴達木盆地周邊地區(qū)由于近年來積極實施植樹造林，改變了原有的荒漠生態(tài)系統(tǒng)格局，植被覆蓋明顯改善；玉樹州中西部由于采取圍欄封育、造林種草、生態(tài)移民等措施，植被覆蓋也得到明顯改善；祁連山東部及三江源區(qū)南部，則由于過度放牧、亂挖亂采，導(dǎo)致植被覆蓋輕微退化；柴達木盆地中部鹽湖和沖積扇的人工植被區(qū)以及湟水谷地，因受農(nóng)田開墾、旅游、道路修建等人類活動影響較大，導(dǎo)致植被退化嚴重；柴達木盆地中部半荒漠帶和柴達木盆地西北部荒漠帶，因氣候干燥、人類活動較少，所以植被覆蓋基本保持不變。

2.4 青海省NDVI變化的可持續(xù)性分析

Hurst指數(shù)可以反映序列自相關(guān)性及序列發(fā)展的相關(guān)強度[25]。將植被覆蓋變化結(jié)果與Hurst指數(shù)疊加，結(jié)果(圖6)顯示，青海省未來植被持續(xù)改善區(qū)域主要分布在柴達木盆地周邊、東昆侖山和共和盆地,持續(xù)性穩(wěn)定不變區(qū)域主要分布在柴達木盆地西北部荒漠帶,持續(xù)性退化區(qū)域零星分布在湟水谷地和柴達木盆地湖泊周邊等人類活動頻繁的地方,未來植被由退化轉(zhuǎn)為改善趨勢的區(qū)域零星分布在青南牧區(qū)，未來變化趨勢不確定區(qū)域主要分布在柴達木盆地中部半荒漠帶，未來植被由改善轉(zhuǎn)為退化趨勢的區(qū)域零星分布在青海湖盆地西部以及柴達木盆地周邊。

Ⅰ.柴達木盆地；Ⅱ.環(huán)湖和祁連山區(qū)；Ⅲ.東部農(nóng)業(yè)區(qū)；Ⅳ.青南牧區(qū)Ⅰ.Qaidam basin；Ⅱ.Central Qinghai lake and Qilian mountains；Ⅲ.Eastern agricultural zone；Ⅳ.Qingnan pastoral area圖5 2000-2015年青海省NDVI的變化趨勢Fig.5 Changing trend of vegetation NDVI in Qinghai from 2000 to 2015

Ⅰ.西寧市；Ⅱ.海西州；Ⅲ.海北州；Ⅳ.海東市；Ⅴ.海南州；Ⅵ.黃南州；Ⅶ.果洛州；Ⅷ.玉樹州Ⅰ.Xining city;Ⅱ.Haixi prefecture;Ⅲ.Haibei prefecture; Ⅳ.Haidong city;Ⅴ.Hainan prefecture;Ⅵ.Huangnan prefecture;Ⅶ.Guoluo prefecture;Ⅷ.Yushu prefecture圖6 青海省不同地區(qū)植被覆蓋未來變化趨勢Fig.6 Future trend of vegetation cover in different regions of Qinghai

將表3植被NDVI變化分析結(jié)果與Hurst指數(shù)結(jié)果疊加，得到青海省未來植被NDVI變化持續(xù)性的分級結(jié)果，可劃分為6種情況，見表4。由表4可知，青海省植被覆蓋未來改善的區(qū)域占整個區(qū)域面積的62.23%，其中持續(xù)性改善區(qū)域占38.01%，由退化轉(zhuǎn)為改善的區(qū)域占24.22%；持續(xù)性穩(wěn)定不變區(qū)域占整個區(qū)域面積的8.43%；植被覆蓋未來有退化趨勢的區(qū)域占整個區(qū)域面積的18.3%，其中持續(xù)性退化區(qū)域占7.74%，由改善變?yōu)橥嘶瘏^(qū)域占10.56%；植被覆蓋未來變化不確定區(qū)域占整個區(qū)域面積的11.04%?？傮w而言，青海省大部分植被未來呈改善趨勢。

表4 青海省植被NDVI變化與可持續(xù)性疊加分析Table 4 Analysis on change and sustainability of vegetation NDVI in Qinghai

2.5 青海省植被覆蓋變化對氣象因子及人類活動的響應(yīng)

由2000-2015年青海省主要地區(qū)年均降雨量及年均氣溫數(shù)據(jù)可知，降雨量較大的區(qū)域主要集中在青海省東南部，東北部的祁連山地區(qū)次之，可可西里及西北的柴達木盆地最少，總趨勢是由東南向西北遞減，與植被覆蓋分布情況較一致；祁連山中西部及青南高原年均溫較低，而中部的湟水谷地及柴達木盆地年均溫較高，總趨勢是南北低中間高。年際變化上，對逐年最大合成NDVI數(shù)據(jù)、青海省內(nèi)35個氣象站點逐年平均氣溫和降水?dāng)?shù)據(jù)(5-9月)進行相關(guān)分析(圖7)發(fā)現(xiàn)，植被覆蓋與降水量、氣溫的相關(guān)系數(shù)分別為0.61和0.62，均通過了顯著性檢驗(P<0.05)。

圖7 2000-2015年青海省NDVI與氣溫、降水及人類活動的關(guān)系Fig.7 Relationship between vegetation NDVI and temperature,precipitation and human activities inQinghai from 2000 to 2015

由圖7-a，b可以看出，2000-2015年氣溫和降水整體呈增加趨勢，2010年以后降水量波動與NDVI波動一致，2008年NDVI最低，這與全省在2008年前后氣溫較低同步。

植被覆蓋變化與人類活動息息相關(guān)，尤其與國家退耕還林還草等大型生態(tài)恢復(fù)政策相關(guān)。由圖7-c可以看出，自1999年國家實行退耕還林還草生態(tài)恢復(fù)政策以來，青海省人工造林面積不斷攀升，同時NDVI整體也呈增加趨勢，說明近年來的生態(tài)恢復(fù)政策對植被覆蓋變化起到了積極作用。但同時青海省的生態(tài)建設(shè)也面臨著嚴重的威脅，首先是草原面積不斷退化，載畜量逐漸下降[7]；其次土壤侵蝕程度日益擴大，水土流失嚴重[26]。青海省肉用畜出欄數(shù)逐年增加，天然草場過度放牧，優(yōu)良牧草因被反復(fù)踐踏、食用，得不到繁衍生息的機會，造成草原不同程度的退化(圖7-d)。

3 結(jié) 論

1)青海省植被覆蓋空間分布整體從西北向東南逐漸增加，NDVI值較低地區(qū)主要集中在柴達木盆地，較高區(qū)域集中在環(huán)湖和祁連山區(qū)、東部農(nóng)業(yè)區(qū)以及青南牧區(qū)東南部，這種格局與青海省氣候特征及地理格局相吻合，NDVI值達0.3以上的區(qū)域占整個區(qū)域的53.7%。

2)16年來青海省植被整體呈增長趨勢，增長率為每10年1.5%；年均NDVI值在0.37～0.42波動，其中2007-2008年NDVI急劇下降，2003-2005年、2008-2010年NDVI快速上升；各分區(qū)植被變化規(guī)律與整體變化規(guī)律基本一致，東部農(nóng)業(yè)區(qū)NDVI增速最快，環(huán)湖和祁連山地區(qū)、柴達木盆地次之，青南牧區(qū)增速最慢；就植被類型而言，旱地、草原與稀樹灌木草原NDVI增速最快，針葉林、荒漠次之，草甸與草本沼澤、灌叢與萌生矮林增速較慢。

3)16年間青海省植被覆蓋改善區(qū)域(60.51%)明顯大于退化區(qū)域(17.87%)，其中植被覆蓋明顯改善的區(qū)域占總面積的21.26%，主要集中在祁連山西部和柴達木盆地周邊；植被覆蓋輕微改善區(qū)域面積最大，占39.25%，除柴達木盆地中部和西北部荒漠區(qū)外均有分布；植被覆蓋狀況基本不變的區(qū)域占21.62%，主要集中在柴達木盆地中部和西北部荒漠區(qū)；植被輕微退化區(qū)域占15.75%，零星分布在青南牧區(qū)東南部、環(huán)湖和祁連山區(qū)東南部及東部農(nóng)業(yè)區(qū)東北部；植被嚴重退化區(qū)域占2.12%，主要分布在柴達木盆地中部鹽湖及沖積扇的人工植被區(qū)、河湟谷地。

4)植被覆蓋變化與Hurst指數(shù)疊加結(jié)果表明，青海省植被覆蓋未來改善的區(qū)域占整個區(qū)域面積的62.23%，其中持續(xù)性改善占38.01%，主要分布在柴達木盆地周邊、東昆侖山、共和盆地，由退化轉(zhuǎn)為改善的區(qū)域占24.22%，零星分布在青南牧區(qū)；穩(wěn)定不變區(qū)域為8.43%；未來變化不確定區(qū)域占整個區(qū)域面積的11.04%，主要分布在柴達木盆地中部半荒漠及西北部荒漠帶；未來有退化趨勢的區(qū)域占整個區(qū)域面積的18.3%，其中持續(xù)性退化的占7.74%，主要分布在柴達木盆地中部的鹽湖及沖積扇的人工植被區(qū)、河湟谷地，由改善變?yōu)橥嘶厔莸膮^(qū)域占10.56%，主要分布在青海湖盆地及柴達木盆地周邊?？傮w而言，青海省大部分地區(qū)植被覆蓋未來呈改善趨勢。

青海省是青藏高原的重要組成部分，近年來，由于人類活動和氣候變化等因素的直接或間接影響，使得青海省各地區(qū)植被覆蓋發(fā)生了不同程度的變化。本研究所用方法較常規(guī)回歸斜率方法能更準(zhǔn)確地反映青海省植被覆蓋變化的趨勢及變化程度，且基于像元尺度描繪了青海省未來植被變化趨勢的可持續(xù)性，對青海省“退耕還林還草”政策的有效實施具有一定的指導(dǎo)意義。但同時研究仍存在很多局限性，首先，植被覆蓋變化具有一定的時效性且氣候變化具有一定階段性，本研究時間跨度較短，對研究結(jié)果有一定的影響；其次，植被覆蓋受氣溫、降水、地形、海拔、土地利用方式及人類活動等多種因素影響，后續(xù)研究需要結(jié)合引起NDVI變化的眾多因子進行綜合分析。