黔中水利樞紐工程區(qū)植被覆蓋變化及其對地形與地貌的響應(yīng)

李 磊, 蘭安軍, 鐘九生

(貴州師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院, 貴陽 550025)

植被覆蓋度是指植被(包括枝、莖和葉)在地面的垂直投影面積在研究區(qū)總面積中的比例，是度量地表植被生長及分布狀況的一個重要指標(biāo)[1]。植被在物質(zhì)與能量循環(huán)、調(diào)節(jié)全球碳平衡和維護(hù)全球生態(tài)系統(tǒng)穩(wěn)定中具有無可取代的作用[2]。因此，在氣候變化加劇背景下，開展植被覆蓋度研究有重要現(xiàn)實意義[3]。

植被覆蓋度估算方法較多，傳統(tǒng)的方法是進(jìn)行地面實測，該方法雖精度高但范圍局限大且難以監(jiān)測長時間序列。遙感技術(shù)具有實時、連續(xù)和監(jiān)測范圍廣等優(yōu)勢[4]，與回歸模型[5]、像元二分模型等[6]方法的結(jié)合，可彌補地面測量的不足，降低人力物力和財力。其中植被指數(shù)和像元二分模型結(jié)合是測量植被覆蓋度常用的方法之一，成果豐碩，如氣候[7]、地形[8]、人類活動等[9]對植被覆蓋度變化響應(yīng)。地形因子(如高程、坡度及坡向等)對植被的分布起著基礎(chǔ)性作用，也是當(dāng)前的研究熱點。但是研究大多只結(jié)合地形因子的等級變化來分析植被覆蓋度的變化規(guī)律，對植被覆蓋變化類型在各地形因子上分布優(yōu)勢定量分析還較少。同時，植被覆蓋變化類型在不同地貌類型和區(qū)域有不同的分布特征和規(guī)律，尤其是在人類大型工程干預(yù)下，植被覆蓋度類型與地形地貌等級間呈現(xiàn)怎樣的空間格局，值得探究。

黔中水利樞紐工程于2009年動工，2016年第一期工程主體工程建設(shè)完成。在人類活動強(qiáng)烈干預(yù)下，原本工程性缺水嚴(yán)重，生態(tài)環(huán)境脆弱的區(qū)域植被覆蓋度變化趨勢是什么樣的，各個地形地貌等級間的植被覆蓋度分布特征如何？故選取此區(qū)域2001—2020年MODIS NDVI數(shù)據(jù)及DEM數(shù)據(jù)來開展研究，以期為區(qū)域水土保持和石漠化治理的科學(xué)實施提供參考，助力當(dāng)?shù)厣鷳B(tài)環(huán)境保護(hù)。

1 研究區(qū)概況

研究區(qū)位于貴州省中、西北部，經(jīng)緯度104°14′—107°10′E，25°28′—27°6′N，總面積16 630. 27 km2，地勢西高東低，海拔為372～2 889 m，大部分處于石漠化敏感區(qū)，生態(tài)系統(tǒng)較脆弱，可利用土地資源有限，開發(fā)強(qiáng)度不高[10]。黔中水利樞紐工程連接長江與珠江，是貴州省首個跨流域和地區(qū)的大型長距離水利調(diào)水工程，是黔中地區(qū)經(jīng)濟(jì)社會發(fā)展的生命線[11]。工程經(jīng)過貴州省貴陽市、安順市、六盤水市、畢節(jié)地區(qū)以及黔南州的部分地區(qū)[12]。地處亞熱帶季風(fēng)區(qū)，多年平均降水約為1 236.54 mm，屬于濕潤區(qū)，最熱月均溫22℃以上，最冷月均溫0～15℃，植被類型為亞熱帶常綠闊葉林，土壤主要為黃壤，兼有非地帶性石灰土[13]。研究區(qū)所在的黔中地區(qū)城市密集、交通通達(dá)度高，人口密度大、耕地資源集中、工業(yè)基礎(chǔ)好，是貴州省的中心地帶，區(qū)位優(yōu)勢十分明顯，發(fā)展?jié)摿薮骩13-15]。

2 材料與方法

2.1 數(shù)據(jù)來源與預(yù)處理

2001—2020年NDVI數(shù)據(jù)來源于美國國家航天局的MOD13Q1產(chǎn)品，空間和時間分辨率分別為250 m和16 d。首先利用MRT軟件將HDF格式的數(shù)據(jù)轉(zhuǎn)換為GeoTIFF，投影定義為UTM，其次采用最大值合成法(MVC，可更好地反映年植被覆蓋情況和消除數(shù)據(jù)中的誤差干擾[16])得到逐年的NDVI產(chǎn)品，最后以黔中水利樞紐工程區(qū)矢量邊界為掩膜進(jìn)行裁剪處理。

利用水文地質(zhì)和巖性圖，參考熊康寧等[17]研究成果，根據(jù)巖性類型[18]，將研究區(qū)劃分為純喀斯特區(qū)、亞喀斯特區(qū)與非喀斯特區(qū)3種地貌類型區(qū)域(圖1A)。

DEM(數(shù)字高程模型)數(shù)據(jù)來源于地理空間數(shù)據(jù)云(分辨率為30 m)，為降低數(shù)據(jù)誤差，保證像元尺度一致，將DEM數(shù)據(jù)的空間分辨率調(diào)整為250 m，再用研究區(qū)矢量邊界將其裁剪后獲得黔中水利樞紐工程區(qū)的DEM數(shù)據(jù)。高程、坡度及坡向數(shù)據(jù)通過DEM數(shù)據(jù)在ArcGIS軟件中提取得到。結(jié)合研究區(qū)實際將高程劃分為8級：≤500 m，500～800 m，800～1 100 m，1 100～1 400 m，1 400～1 700 m，1 700～2 000 m，2 000～2 300 m，>2 300 m(圖1B)；根據(jù)自然資源部頒布的《第三次全國土地調(diào)查技術(shù)規(guī)程》，將坡度劃分為5級：≤2°，2°～6°，6°～15°，15°～25°，>25°(圖1C)；坡向數(shù)據(jù)劃分為平地、東坡(45°～135°)、南坡(135°～225°)、西坡(225°～315°)和北坡(315°～45°)(圖1D)。

圖1 地貌分布和地形分級

2.2 研究方法

2.2.1 像元二分模型估算法 基于像元二分模型提取研究區(qū)植被覆蓋度(fc)，計算公式如下[19]：

fc=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)

(1)

式中:NDVIsoil為無植被覆蓋像元的NDVI值;NDVIveg為純植被覆蓋像元的NDVI值。參考相關(guān)研究并結(jié)合研究區(qū)實際情況[20],以累積頻率為5%的NDVI值為NDVIsoil,95%的NDVI值為NDVIveg。

2.2.2 植被覆蓋度變化類型提取模型 逐像元分析研究區(qū)植被覆蓋度變化，擬合其變化斜率，獲取植被覆蓋度的變化類型，計算公式如下[21]：

(2)

式中:n為監(jiān)測總年數(shù),在本研究中為20 a;fi為第i年的植被覆蓋度;slope為線性擬合斜率,可對黔中水利樞紐工程區(qū)20 a植被覆蓋度的變化趨勢進(jìn)行表征,其值為正,表示植被覆蓋度增加,反之則表示植被覆蓋度降低。

采用F檢驗對顯著性大小進(jìn)行檢驗，計算公式為：

(3)

2.2.3 分布指數(shù) 為明晰植被覆蓋度變化的空間分異特征及其對不同地形與地貌的響應(yīng)，引入分布指數(shù)來進(jìn)行計算，公式如下[22]：

P=(Sie/Si)/(Se/S)

(4)

式中:P為分布指數(shù),分布優(yōu)勢和P值兩者具有正相關(guān)性;Sie為在e區(qū)域內(nèi)i類型的面積;i為植被覆蓋變化的類型;Si為研究區(qū)內(nèi)i類型的總面積;Se表示e區(qū)域(本研究中指地形或地貌區(qū)域)的面積;S為研究區(qū)的總面積;當(dāng)P=1時,表明該類型在e區(qū)域上平穩(wěn)分布;P<1時,表明該類型在e區(qū)域為劣勢分布,反之為優(yōu)勢分布,且P值越大,分布優(yōu)勢越顯著。

3 結(jié)果與分析

3.1 年際變化分析

3.1.1 植被覆蓋度等級分布及數(shù)值年際變化 采用基于像元二分模型的植被覆蓋估算方法，對黔中水利樞紐工程區(qū)的植被覆蓋度進(jìn)行估算，利用研究時段內(nèi)逐年的最大植被覆蓋度值計算得出20 a的平均植被覆蓋度(圖2)。為便于進(jìn)一步分析研究，將所得到的20 a平均植被覆蓋度柵格數(shù)據(jù)在ArcGIS10.2軟件中進(jìn)行重分類，參考水利部頒布的《土壤侵蝕分類分級標(biāo)準(zhǔn)》(SL190—2007)，并根據(jù)研究區(qū)實際情況將植被覆蓋度劃分為低植被覆蓋區(qū)(fc<10%)、較低植被覆蓋區(qū)(10%≤fc<30%)、中植被覆蓋區(qū)(30%≤fc<50%)、較高植被覆蓋區(qū)(50%≤fc<70%)和高植被覆蓋區(qū)(fc≥70%)5個等級區(qū)間。由圖2可以看出，黔中水利樞紐工程區(qū)植被覆蓋分布具有顯著的差異，總體呈西高東低、南北高中部低的分布格局。黔中水利樞紐工程區(qū)以高和較高植被覆蓋度為主，其分別占研究區(qū)面積的70.28%，21.46%，中植被覆蓋區(qū)占5.73%，較低植被覆蓋區(qū)占1.64%，低植被覆蓋區(qū)面積占0.89%，植被覆蓋情況總體較好。

圖2 研究區(qū)2001-2020年平均植被覆蓋度空間分布

圖3為研究區(qū)植被覆蓋度均值的年際變化圖。20 a間研究區(qū)植被覆蓋度平均值為77.86%，年平均增長率為0.46%，總體趨勢表現(xiàn)為波動上升。年平均植被覆蓋度由2001年的75.05%上升至2020年的82.25%，兩年相差7.2%，有較大幅度的提升。其中植被覆蓋度最低的是2005年(73.08%)，低值可能與西南地區(qū)2005—2006年的干旱天氣有關(guān)[23]，2005—2012年植被覆蓋度顯著提升，與當(dāng)?shù)赝烁€林工程的實施有關(guān)，植被覆蓋度最高的是2020年(82.25%)。

黔中水利樞紐工程于2009年動工，2016年第一期工程主體工程建設(shè)完成，在人類大型工程的干預(yù)下，工程建設(shè)中期的2012—2013年植被覆蓋度有所下降，到2013年達(dá)到一個低值(75.88%)，水利樞紐工程建成以后，通過樞紐對研究區(qū)水資源的合理調(diào)配，研究區(qū)的植被覆蓋持續(xù)增加，說明水利工程的建設(shè)有利于研究區(qū)生態(tài)環(huán)境的良性發(fā)展。

圖3 2001-2020年研究區(qū)植被覆蓋度年際變化

3.1.2 研究區(qū)植被覆蓋度變化類型總體分布特征 為進(jìn)一步探究黔中水利樞紐工程區(qū)植被覆蓋度的變化趨勢，運用植被覆蓋變化的提取模型(公式2)提取2001—2020年的植被覆蓋度，得到研究區(qū)植被覆蓋度20 a變化斜率值(圖4A)；結(jié)合斜率值采用F檢驗法(公式3)，劃分研究區(qū)植被覆蓋度的變化類型(圖4B)。

F檢驗的結(jié)果表明，植被穩(wěn)定型占據(jù)研究區(qū)面積的47.82%，是研究區(qū)占比最大的植被覆蓋度變化類型，該類型分布較為均勻；其次為植被改善的區(qū)域，植被改善型占研究區(qū)面積的39.79%，主要分布在研究區(qū)西北部和南部，多為非喀斯特地區(qū)；植被退化型占研究區(qū)面積的12.39%，主要分布在研究區(qū)中部和東部，多為經(jīng)濟(jì)社會較為發(fā)達(dá)的地區(qū)，也是工程建設(shè)的主體區(qū)域，退化原因與水利工程建設(shè)等人類活動的干預(yù)有關(guān)，需得到足夠重視，應(yīng)當(dāng)大力實施封山育林，全面開展石漠化治理，使工程建設(shè)區(qū)周邊植被覆蓋率得以提升，有效遏制水土流失，確保研究區(qū)生態(tài)安全?？傮w來說，研究區(qū)植被覆蓋度20 a的變化情況較為穩(wěn)定，總體較好，且改善的比例大于退化的比例。

圖4 2001-2020年研究區(qū)植被覆蓋度年際變化斜率及變化類型

3.1.3 研究區(qū)植被變化類型區(qū)(縣)差異性分析 為更深入探索黔中水利樞紐區(qū)植被覆蓋變化類型的具體差異，統(tǒng)計28個區(qū)(縣)(僅統(tǒng)計研究區(qū)范圍內(nèi)涉及部分)不同植被覆蓋度變化類型在對應(yīng)區(qū)域的面積比(表1)。

各區(qū)(縣)的面積比例變化差異較大，其中25個區(qū)(縣)的植被改善面積比例大于植被退化面積比例，開陽、清鎮(zhèn)、關(guān)嶺、紫云、赫章、納雍、威寧、織金、惠水、貞豐、長順、六枝、水城、鐘山的改善與退化面積比例差較為顯著，說明這些區(qū)(縣)20 a來的生態(tài)建設(shè)有了顯著成效；而白云區(qū)、觀山湖區(qū)、修文縣的植被退化面積比例超過了植被改善的面積比例，說明這3個區(qū)(縣)(研究區(qū)范圍內(nèi)涉及部分)在經(jīng)濟(jì)發(fā)展的同時，存在植被退化風(fēng)險，需要加強(qiáng)對生態(tài)環(huán)境的治理保護(hù)。

3.2 研究區(qū)植被覆蓋度變化類型的地形分布特征

將研究區(qū)植被變化類型分別與各地形因子(高程、坡度、坡向)進(jìn)行空間分析，運用分布指數(shù)計算不同植被覆蓋度變化類型在各地形因子上的分布情況，探究和分析其分布特征。

3.2.1 植被覆蓋度變化類型隨高程分布特征 隨著高程的變化，不同高度上接受的輻射量、土壤性質(zhì)與氣流狀況等自然要素也會發(fā)生變化[24]。從圖5A可以看出，隨高程上升，植被覆蓋度變化類型也有不同程度的變化，其總體變化趨勢為：植被改善型先增后減，之后再增再減；植被穩(wěn)定型先增再減后又回升；植被退化型的P值則是先減后增，后再次減少后又增加。植被穩(wěn)定型先增再減后又回升。植被改善型分別在700～1 000 m和大于1 700 m的區(qū)域呈優(yōu)勢分布，其余高程區(qū)域呈劣勢分布，其中在高程為2 300 m時植被改善最為顯著，大于2 300 m時呈下降趨勢，但未降至1以下。植被退化型在500 m以下時，分布優(yōu)勢明顯，500～800 m區(qū)間P值迅速降低至1以下，當(dāng)高程處于800～1 400 m時又有所回升，大于1 400 m再次下降，其中在1 300～1 600 m區(qū)域內(nèi)植被退化優(yōu)勢分布，原因在于農(nóng)村居民點在此區(qū)域分布最集中，受生活和工程建設(shè)等人類活動的影響較大[25]，導(dǎo)致此區(qū)域的植被退化情況較為顯著。植被穩(wěn)定型在高程小于800 m時，P值有所上升，趨近平穩(wěn)分布，在800～1 700 m的高程區(qū)域上為優(yōu)勢分布，1 700～2 300 m區(qū)域內(nèi)為劣勢分布，當(dāng)高程超過2 300 m后，P值又逐漸回升。

表1 研究區(qū)范圍各區(qū)(縣)植被覆蓋度變化情況所占該區(qū)域比例 %

圖5 植被覆蓋度變化類型P值隨高程的變化趨勢和面積比

從圖5B可以看出，研究區(qū)各高程區(qū)域中，植被改善型和植被穩(wěn)定型均占有較大比例，其合計占比達(dá)77%以上，在這其中又以植被穩(wěn)定型占比更多，占比均值為46.6%。植被退化型的占比在各高程區(qū)域上都不高，以高程≤500 m的22.99%為最高，2 300 m的1.93%為最小。由上可知，在≤500 m和800～1 600 m的高程上，植被退化明顯，這與該區(qū)域內(nèi)人類活動如人居生活、生產(chǎn)建設(shè)等的干擾有關(guān)，在其他高程上，植被穩(wěn)定和植被改善明顯，原因是該區(qū)域地勢陡峭，人類活動的干擾較少，同時，研究區(qū)生態(tài)保護(hù)及生態(tài)修復(fù)工程使得植被得以保留。

3.2.2 植被覆蓋度變化類型隨坡度分布特征 坡度大小對地表能量轉(zhuǎn)換和物質(zhì)流動有直接影響[26]。由圖6A可看出，隨著坡度的增加，植被改善型先增加，后增速逐漸減緩；植被退化型持續(xù)降低，坡度大于15°后降速逐漸減緩；植被穩(wěn)定型則是先增后減，但幅度不大。其中，植被改善型在15°～25°時，P值逐漸趨近并達(dá)到優(yōu)勢分布，但超過25°時增速逐漸減緩；植被退化型表現(xiàn)為兩極分化，坡度在15°以下時為優(yōu)勢分布，坡度超過15°時為劣勢分布；植被穩(wěn)定型在各坡度區(qū)域上均平穩(wěn)分布(P值圍繞1上下波動，幅度不大)。

圖6 植被覆蓋度變化類型P值隨坡度的變化趨勢和面積比

從圖6B可以看出，研究區(qū)各坡度范圍內(nèi)，植被改善型和植被穩(wěn)定型占比較大，其合計占比均在61.48%以上，其中又以植被穩(wěn)定型占比較大，占比均值為46.12%。隨著坡度的升高，植被退化型占比持續(xù)下降，其在坡度≤2°的區(qū)域內(nèi)比例最大，為38.52%，坡度大于15°后降至1%以下，在坡度大于25°的區(qū)域占比最小，僅0.33%。植被覆蓋度變化受坡度的影響較大[27]，坡度小于6°時，地勢起伏小，人類活動干擾較其他區(qū)域大。坡度在6°～15°時，植被改善和植被穩(wěn)定均呈現(xiàn)上升趨勢，但未達(dá)到優(yōu)勢分布，植被改善并不顯著，說明此坡度區(qū)域內(nèi)植被處于不穩(wěn)定狀態(tài)，需加強(qiáng)生態(tài)保護(hù)，避免植被發(fā)生退化。坡度大于15°時地勢起伏大，人類活動影響小，植被改善類型在此區(qū)域內(nèi)呈優(yōu)勢分布。

3.2.3 植被覆蓋度變化類型隨坡向分布特征 不同坡向的土壤含水量和太陽輻射量存在較大的差異，對植被的分布特征也有較大影響[8]。由圖7A可看出，當(dāng)北坡轉(zhuǎn)換到東坡時，植被改善型的P值上升，植被退化和植被穩(wěn)定型的P值均呈下降趨勢；當(dāng)西坡轉(zhuǎn)換到南坡時，植被改善型和植被穩(wěn)定型的變化趨勢均保持平穩(wěn)，植被退化型稍有上升，但變化不大。其中，植被改善型呈優(yōu)勢分布的坡向為南坡和西坡；植被退化型呈優(yōu)勢分布的坡向為平地、北坡和東坡，在其中又以平地優(yōu)勢最大，北坡和東坡優(yōu)勢略差；植被穩(wěn)定型在平地劣勢顯著，除平地外分布平穩(wěn)。

圖7 植被覆蓋度變化類型P值隨坡向的變化趨勢和面積比

從圖7B可以看出，各坡向上，植被改善型和植被穩(wěn)定型的占比均較大，合計占比達(dá)75.56%以上。隨著坡向的轉(zhuǎn)變，植被覆蓋變化類型也有不同程度的變化，北坡轉(zhuǎn)向東坡時，植被退化和植被穩(wěn)定型的面積占比略有下降，植被改善型占比稍有上升；西坡向南坡轉(zhuǎn)變時，植被退化輕微上升，變化不大。植被改善型占比最高的坡向為南坡，占比為41.27%，植被退化型占比最高的坡向是平地，為24.44%。

由上可知，各植被變化類型在平地上具有顯著分異特征，其余坡向則差異微弱，在陰坡(北坡和西坡)上，植被退化類型北優(yōu)西劣，植被改善類型北劣西優(yōu)；在陽坡(東坡和南坡)上，植被改善型在南坡略有優(yōu)勢，在東坡趨于穩(wěn)定，植被退化型在東坡優(yōu)勢稍顯，在南坡趨于穩(wěn)定。當(dāng)坡向由陰坡轉(zhuǎn)為陽坡時，植被穩(wěn)定和退化類型的分布指數(shù)及面積比略有下降，植被改善型則略有上升，但變化幅度都不大。原因是與陰坡相比，陽坡的水分更合適、太陽輻射更充足，更適宜植被的生長發(fā)育。

3.3 植被覆蓋度變化類型的地貌分布特征

由圖8A可以看出，植被改善型在非喀區(qū)為優(yōu)勢分布，在亞喀區(qū)域和純喀區(qū)為劣勢分布；植被退化型在非喀表現(xiàn)為劣勢分布，亞喀和純喀區(qū)表現(xiàn)為優(yōu)勢分布；植被穩(wěn)定型在各地貌區(qū)平穩(wěn)分布，波動不大。非喀區(qū)優(yōu)勢分布的是植被改善型，劣勢分布的是植被退化型，說明非喀地貌區(qū)植被改善較為顯著；純喀區(qū)優(yōu)勢分布的是植被退化型，劣勢分布的是植被改善型，且又以亞喀區(qū)較為顯著，說明純喀和亞喀地貌區(qū)的植被退化比非喀區(qū)更為顯著。

圖8 植被覆蓋度變化類型P值在不同地貌下的分布指數(shù)變化趨勢和面積比

由圖8B可知，各地貌區(qū)，植被改善型和植被穩(wěn)定型的占比均較大，其合計占比均在85%以上，又以植被穩(wěn)定型占比較大，均在47%以上。植被退化型在各地貌類型中面積占比均較低，亞喀斯特地貌的14.14%最高，非喀區(qū)的8.23%最小。在純喀區(qū)和亞喀區(qū)植被退化較為明顯，原因是喀斯特區(qū)域生態(tài)系統(tǒng)的脆弱性，導(dǎo)致喀斯特地貌區(qū)更容易遭到破壞，且破壞后恢復(fù)困難。非喀區(qū)內(nèi)植被改善和植被穩(wěn)定顯著。

4 討論與結(jié)論

4.1 討 論

研究中采用的像元二分模型，計算較為簡單且不受地域限制，所提取的植被覆蓋度能夠較為全面地反映研究區(qū)植被覆蓋的生長及分布情況，但該方法僅以置信度區(qū)間來截取NDVI上下閾值，具有一定局限性。此外，本研究中所使用MODIS數(shù)據(jù)MOD13Q1級NDVI產(chǎn)品已排除云和大氣的干擾，且與其他各項數(shù)據(jù)保證了像元尺度的一致性，能夠較為準(zhǔn)確地反映研究區(qū)植被覆蓋度的時空分異特征，但眾多的中間過程難免導(dǎo)致植被覆蓋變化的模擬結(jié)果有一定誤差，盡管如此，本文力求使用科學(xué)合理的數(shù)據(jù)與模型方法，來分析黔中水利樞紐工程區(qū)的植被覆蓋變化分異特征，不足之處仍有待進(jìn)一步探索和研究。同時，本文只分析了研究區(qū)不同地形、地貌的植被覆蓋變化分異特征，而沒有對其相互之間的關(guān)系進(jìn)行更深入的分析。研究區(qū)地形、地貌條件復(fù)雜，未來需要進(jìn)一步深入研究不同地形與地貌的植被覆蓋變化的具體差異原因，以便更為精確、更加全面地得到黔中水利樞紐工程區(qū)植被生長的分異特征和影響因子，為研究區(qū)水土保持和石漠化治理及大型調(diào)水工程的科學(xué)實施提供更加完備而科學(xué)的數(shù)據(jù)支撐。

4.2 結(jié) 論

黔中水利樞紐工程區(qū)2001—2020年，年平均植被覆蓋狀況較好，不同地形(高程、坡度、坡向)和地貌(純喀斯特、亞喀斯特、非喀斯特)條件下，植被覆蓋空間分布及變化存在顯著的分異規(guī)律，主要有以下特征：

(1) 20年間，研究區(qū)植被覆蓋情況總體較好。年平均增長率為0.46%，呈波動上升趨勢，尤其在工程干預(yù)之后，原本工程性缺水嚴(yán)重，生態(tài)環(huán)境脆弱的區(qū)域，生態(tài)環(huán)境得到不斷改善。在20 a平均植被覆蓋度分級中，低、較低、中、較高和高植被覆蓋區(qū)分別占工程區(qū)總面積的0.89%，1.64%，5.73%，21.46%，70.28%。

(2) 研究區(qū)植被穩(wěn)定型占比最高，達(dá)47.82%；植被改善型次之，為39.79%；植被退化型最低，僅為12.39%。除白云區(qū)、觀山湖區(qū)、修文縣3個區(qū)(縣)外，其余區(qū)(縣)的植被改善面積均超過植被退化面積。

(3) 隨高程和坡度上升，植被退化型的分布優(yōu)勢降低。隨坡向變化，平地區(qū)植被退化顯著，其余坡向差異不明顯；從地貌上看，非喀區(qū)植被改善比純喀和亞喀區(qū)顯著，純喀和亞喀地貌區(qū)較非喀區(qū)植被退化突出。