麗江壩區(qū)植被覆蓋度時空變化驅(qū)動力研究

田瀟然，王 錦，余哲修

(1.西南林業(yè)大學(xué) 園林園藝學(xué)院，云南 昆明 650233；2.昆明理工大學(xué) 建筑與城市規(guī)劃學(xué)院，云南 昆明 650500；3.北京林業(yè)大學(xué) 林學(xué)院，北京 100083)

植被是地球生態(tài)系統(tǒng)的重要組成部分，為人類生存提供必需的食物資源和物質(zhì)環(huán)境，植被的生長狀態(tài)、分布特征與人類的生活生產(chǎn)、社會經(jīng)濟發(fā)展都密切相關(guān)[1]。植被覆蓋度(fractional vegetation cover，F(xiàn)VC)是植被生長狀況以及地表植被占比面積的重要量化指標，常用于生態(tài)評價和景觀生態(tài)分析的研究[2]。人為活動頻繁導(dǎo)致了植被的退化、景觀破碎化日益凸顯，人地間、人與生態(tài)環(huán)境間的矛盾加劇，協(xié)調(diào)城鎮(zhèn)快速擴張與植被[3]、生態(tài)環(huán)境之間的平衡成為當下的研究熱點。

近年來，國內(nèi)外學(xué)者針對植被覆蓋度分布變化特征與影響因素的研究有了較大進展，如楊旭超等[4]通過小波分析、趨勢分析法獲得昆明市呈貢區(qū)植被覆蓋度和土地的變化特征；賈路等[5]通過對西安市植被覆蓋度的變化特征、空間重心轉(zhuǎn)移進行分析，獲得城市植被覆蓋度的空間變化差異；達佤扎喜等[6]從城鎮(zhèn)化角度出發(fā)，對西寧市的景觀格局和植被覆蓋度的時空變化進行分析，探索城鎮(zhèn)發(fā)展對植被覆蓋度變化的影響；Y.H.Jinetal[7]構(gòu)建了半?yún)?shù)地理加權(quán)回歸模型，從縣域?qū)用娣治隽擞绊憦V東省植被覆蓋度變化的驅(qū)動因素。研究者采用不同方法探討了城市植被覆蓋度的變化特征和影響因素，然而針對由旅游業(yè)推動下的快速城鎮(zhèn)化與植被覆蓋度變化間的關(guān)聯(lián)性、驅(qū)動機制的研究則較少。

自1997年麗江古城列入世界文化遺產(chǎn)名錄以來，當?shù)芈糜螛I(yè)從探索、起步進入到快速發(fā)展階段[8]。旅游業(yè)拉動的經(jīng)濟與就業(yè)機會促進了當?shù)禺a(chǎn)業(yè)結(jié)構(gòu)的變更，同時也加速了城鎮(zhèn)化的進程。麗江地區(qū)于2003年撤地設(shè)市，于2009年麗江完成了城市建設(shè)的“三級跳”[9]，快速城鎮(zhèn)化使得當?shù)刂脖坏臅r空格局發(fā)生了明顯的變化。麗江作為旅游發(fā)展研究的典型案例，在城鎮(zhèn)化影響下植被變化趨勢、植被覆蓋度空間格局的探討則鮮有研究者關(guān)注。

本研究以麗江市的核心麗江壩區(qū)為對象，對其植被覆蓋度的計算結(jié)果進行同心圓空間梯度分析，結(jié)合回歸趨勢法，綜合探討植被覆蓋度時空變化特征，探析自然氣象變化、旅游經(jīng)濟發(fā)展、快速城鎮(zhèn)化、產(chǎn)業(yè)結(jié)構(gòu)調(diào)整等因素與植被生長狀態(tài)以及空間分布間的關(guān)聯(lián)性，為地區(qū)發(fā)展規(guī)劃、環(huán)境保護以及生態(tài)修復(fù)提供支撐和參考依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

麗江壩區(qū)位于麗江市中部(100°10′0″-100°20′0″E，26°45′0″-27°05′0″N)，為南北長、東西較短的山間盆地(圖1)。壩區(qū)內(nèi)地勢平坦，四周山體環(huán)繞，土地面積194.39 km2。壩區(qū)內(nèi)下轄古城區(qū)的大研街道、西安街道、祥和街道、束河街道、金山白族鄉(xiāng)，以及玉龍納西自治縣黃山鎮(zhèn)與白沙鎮(zhèn)的部分地區(qū)。屬于暖溫帶高原山地季風(fēng)氣候，年平均降水量1 108.8 mm，全年降水不均，集中于夏季，年平均氣溫13.46 ℃。麗江壩區(qū)是麗江市政治、經(jīng)濟、文化的核心區(qū)，早期壩區(qū)內(nèi)以農(nóng)業(yè)經(jīng)濟為主，在大力開發(fā)旅游業(yè)后，大研古鎮(zhèn)、束河古鎮(zhèn)、白沙古鎮(zhèn)等成為全市旅游熱點區(qū)域[10]，旅游業(yè)逐漸替代農(nóng)業(yè)成為當?shù)氐闹еa(chǎn)業(yè)。

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

研究數(shù)據(jù)為1987-2018年連續(xù)32 a的Landsat影像，來源于美國地質(zhì)勘探局(United States geological survey，USGS)。植被生長旺季(6-9月)為麗江地區(qū)的雨季，影像云量較大，因此選用非雨季1-5月的影像作為研究數(shù)據(jù)。影像在大氣校正、裁剪等預(yù)處理后，計算每期影像的NDVI值，并形成NDVI時間序列。由于NDVI時間序列數(shù)據(jù)存在噪聲，因此運用局部加權(quán)回歸法(locally weighted regression，LOWESS)進行平滑處理去除噪聲干擾，重新擬合時間曲線。平滑處理后的NDVI時序曲線中出現(xiàn)了1991、1994、2010年3個主要突變點，結(jié)合研究的起止年份，確定1987、1991、1994、2010、2018年為5個研究時間節(jié)點(圖2、表1)。

表1 研究區(qū)遙感影像信息

1.3 研究方法

1.3.1 植被覆蓋度估測 運用像元二分模型計算植被覆蓋度。該方法假設(shè)1個像元的地表有植被覆蓋和無植被覆蓋兩部分所組成，植被覆蓋度(FVC)為植被面積在像元中所占的比例[11]。其計算公式為

(1)

式中：NDVI為植被指數(shù)(NDVI)的值，NDVIveg為被綠色植被覆蓋區(qū)域的NDVI值，NDVIsoil為沒有被植被覆蓋地區(qū)的NDVI值。參考相關(guān)研究成果[12-13]，提取像元NDVI值累積分布處于95%分位處的像元NDVI均值為NDVIveg，提取累積分布處于5%分位處的像元NDVI均值為NDVIsoil，估算麗江壩區(qū)內(nèi)的FVC。

1.3.2 趨勢分析 采用一元線性回歸分析法獲得麗江壩區(qū)連續(xù)32 a植被覆蓋度的變化趨勢，探討其變化特征。計算公式為

(2)

1.3.3 空間梯度分析 采用同心圓空間梯度模型進行植被覆蓋度的梯度分析，根據(jù)麗江壩區(qū)的整體形態(tài)確定其幾何中心，以幾何中心為圓心、1 km為間距，將壩區(qū)劃分為18個同心圓采樣環(huán)帶(圖3)，對每個采樣環(huán)帶進行景觀格局指數(shù)計算以及梯度分析[14-15]。

1.3.4 景觀格局指數(shù)選擇和計算 植被在快速城鎮(zhèn)化的影響下逐漸呈現(xiàn)出不規(guī)則、破碎化的分布狀態(tài)，本研究在參考前人成果的基礎(chǔ)上結(jié)合麗江壩區(qū)的發(fā)展具體情況，選擇類型水平指數(shù)中的NP(斑塊數(shù)量)、PD(斑塊密度)、AMFRAC(面積加權(quán)平均分維數(shù))、COHESION(連接度指數(shù))、AI(聚集度指數(shù))5個指數(shù)量化植被覆蓋度的時空分布特征以及破碎化程度。選用的景觀指數(shù)定義及計算公式見文獻[16-18]。

1.3.5 驅(qū)動力分析

1.3.5.1 驅(qū)動因子選取 根據(jù)麗江壩區(qū)植被覆蓋度變化特征及城鎮(zhèn)化發(fā)展進程，從歷年《云南統(tǒng)計年鑒》《云南年鑒》和《麗江年鑒》中篩選出能夠代表自然生態(tài)、人口、社會經(jīng)濟變化的3類因素。由于影響植被覆蓋度時空變化的因子眾多[19]，其中自然與人為因素的影響程度較為突出。自然因素中的降水、日照、氣候變化等直接關(guān)系到植被的生長狀態(tài)[20-22]，同時城市的擴張、人為活動都能在短期內(nèi)改變植被的分布狀態(tài)以及空間格局[23-25]，而城市擴張與人口數(shù)量、國民生產(chǎn)總值、財政收入等社會經(jīng)濟、城鎮(zhèn)化水平因素密切相關(guān)。因此，選取19個因子進行驅(qū)動力分析，并對其進行分類(表2)。

表2 麗江壩區(qū)植被覆蓋度變化驅(qū)動因子

1.3.5.2 偏最小二乘回歸 偏最小二乘回歸(partial least squares regression，PLSR)集合了多元線性回歸分析、主成分分析和典型相關(guān)分析的基本特征和優(yōu)點，在樣本數(shù)量較小時能夠有效反映出多個自變量或因變量之間的相關(guān)性，同時能排除系統(tǒng)內(nèi)的干擾，多次成分提取不剔除變量[26-28]，保證了信息的完整性。偏最小二乘回歸分析在土地利用、景觀格局、綠色空間變化等領(lǐng)域的研究中都得到了應(yīng)用[29]。

將自然、人口、社會經(jīng)濟三方面選取出來的因子作為自變量X，植被覆蓋度作為因變量Y，用偏最小二乘回歸定量分析自然、人口、社會經(jīng)濟因素與植被覆蓋度變化之間的關(guān)系，同時從定性角度解析32 a間城鎮(zhèn)化建設(shè)、旅游開發(fā)、政府決策等對壩區(qū)植被覆蓋度變化的影響。

2 結(jié)果與分析

2.1 植被覆蓋度時空變化特征

對1987-2018年的植被覆蓋度逐年分析，獲得32 a連續(xù)變化趨勢計算結(jié)果，并通過四分位法分為4個等級(表3)，繪制變化趨勢分布圖(圖4)。

表3 植被覆蓋度變化趨勢分級閾值標準

從結(jié)果可以看出，4個等級所占的面積都較為均衡，其中明顯改善的植被面積略高。壩區(qū)中主城區(qū)與玉龍新城內(nèi)，以及村落周邊區(qū)域的植被發(fā)生了明顯的退化現(xiàn)象，輕微退化與輕微改善的區(qū)域主要交錯分散于主城區(qū)、玉龍新城外圍以及壩區(qū)南北兩端的耕地中，而遠離城市的壩區(qū)邊緣的植被則得到明顯改善。植被覆蓋度整體變化特征呈現(xiàn)出：人為活動干擾越大退化程度越強，遠離城市的半自然環(huán)境中的植被得到恢復(fù)。

2.2 同心圓空間梯度變化分析

結(jié)合壩區(qū)植被生長狀況和遙感計算結(jié)果將植被覆蓋度劃分為5個等級，用以對植被覆蓋度時空變化的分布特征進行定量分析(表4)。

表4 植被覆蓋度分級閾值標準

由圖5可知，5個等級的斑塊數(shù)量整體變化趨勢都較為相似，從1～6環(huán)帶為快速增加，6環(huán)帶斑塊數(shù)量達到最高，6～18環(huán)帶都呈現(xiàn)下降特征，說明主城區(qū)與玉龍新城之間的耕地內(nèi)植被斑塊數(shù)量較多且分布零散，并未連成片。

1987-2018年，5個等級的斑塊密度總體出現(xiàn)了明顯的波動，1～2環(huán)帶快速下滑，2～16環(huán)帶變化緩慢，17～18環(huán)帶變化各不相同(圖6)。L-FVC中，2018年1～9環(huán)帶斑塊密度最高。LM-FVC、M-FVC、HM-FVC和H-FVC變化中，2010年2～6環(huán)帶的斑塊密度總體最低，由于1994-2010年主城區(qū)面積快速擴張，其周邊的耕地轉(zhuǎn)為建設(shè)用地，導(dǎo)致植被大面積減少。HM-FVC、H-FVC中，2018年9～18環(huán)帶斑塊密度較高。麗江市從2006年啟動創(chuàng)建國家級園林城市工作開始至2012年“創(chuàng)園”成功，綠化建設(shè)和生態(tài)環(huán)境受到重視，植被得以逐步恢復(fù)，因此5個等級中2018年的斑塊密度總體高于2010年，但植被斑塊密度的增加也表明了破碎度呈現(xiàn)加劇趨勢。

面積加權(quán)平均斑塊分維數(shù)的計算中(圖7)，5個等級都呈現(xiàn)出先快速上升，后平緩變化，再快速下降的特征。L-FVC中，2018年2～18環(huán)帶中整體下降趨勢最顯著；LM-FVC中，2～8環(huán)帶內(nèi)2018年的斑塊形狀較為復(fù)雜；M-FVC中，2～8環(huán)帶內(nèi)1991年與1994年的斑塊形狀較規(guī)則，2018年最為復(fù)雜。HM-FVC的變化中，1987年1～7環(huán)帶的斑塊形狀更復(fù)雜；而H-FVC中2018年1～7環(huán)帶的斑塊形狀更規(guī)則。說明1987-2018年主城區(qū)內(nèi)及周邊區(qū)域的LM-FVC、M-FVC的植被斑塊受到建設(shè)活動干擾形狀逐漸趨于復(fù)雜，隨著環(huán)境改造、城市綠化的加強，MH-FVC、H-FVC的斑塊形狀逐漸規(guī)則，人工景觀特征更為顯著。

由圖8可知，5個等級的景觀連接度變化特征各不相同，L-FVC中，2018年斑塊的連接度整體低于其他4個時間點，8～18環(huán)帶中尤為突出。LM-FVC、M-FVC變化中，1987-2018年3～16環(huán)帶整體呈現(xiàn)緩慢升高的趨勢，說明主城區(qū)內(nèi)的植被斑塊連接度低于郊區(qū)，同時2018年的連接度整體較高。主城區(qū)內(nèi)的綠化種植以點狀、線狀方式分散布局，主城區(qū)周邊的耕地在城市擴張、道路新建的過程中被割裂，植被連接度降低，HM-FVC、H-FVC中，2018年2～8環(huán)帶的斑塊連接度總體較低。

由圖9可知，5個等級植被覆蓋度的聚集度變化各不相同。L-FVC中，2018年的斑塊聚集度整體低于其他4個時間節(jié)點。LM-FVC、M-FVC中，1987-2018年1～18環(huán)帶的聚集度總體呈現(xiàn)逐漸增加的趨勢，其中2010、2018年1～9環(huán)帶的聚集度高于1987、1991、1994年，說明2010、2018年主城區(qū)及周邊區(qū)域LM-FVC、M-FVC的斑塊分布較為密集。HM-FVC、H-FVC中，1987、1991、1994、2010年2～8環(huán)帶的斑塊聚集度整體高于2018年，說明隨著城市建筑、道路、人工設(shè)施不斷增加，植被覆蓋度較高的斑塊分布變得更加松散。

2.3 壩區(qū)植被覆蓋度時空變化驅(qū)動力分析

2.3.1 驅(qū)動力主成分分析 經(jīng)偏最小二乘回歸分析，根據(jù)投影重要性(variable important projection，VIP)大小，同時結(jié)合精度分析，1～4個主成分的貢獻率分別為：71.70%、7.50%、4.30%、7.40%。截止到前4個主成分，其累計貢獻率達到90.90%(表5)，已經(jīng)涵蓋了19個原始因子較為充足的信息。因此選取前4個主成分進行偏最小二乘回歸分析。

表5 主成分與研究項精度分析

提取的1～4個主成分分別用字母A、B、C、D表示，各主成分與原始驅(qū)動因子的表達式為

A=-0.182X1+0.052X2-0.112X3-0.282X4+0.294X5+0.313X6+0.054X7+0.366X8+0.240X9+0.211X10+0.240X11+0.157X12+0.240X13+0.209X14+0.290X15+0.171X16+0.233X17+0.240X18+0.202X19

(3)

B=0.147X1+0.303X2-0.128X3+0.865X4-0.180X5-0.225X6+0.022X7-0.382X8+0.012X9+0.063X10+0.013X11+0.183X12+0.008X13-0.005X14-0.148X15+0.171X16+0.025X17+0.013X18+0.052X19

(4)

C=0.172X1-0.021X2+0.494X3+0.142X4-0.089X5-0.240X6+0.263X7-0.475X8+0.139X9+0.143X10+0.139X11+0.309X12+0.119X13-0.103X14-0.462X15+0.301X16+0.114X17+0.132X18+0.283X19

(5)

D=-0.223X1+0.428X2+0.247X3+0.061X4+0.203X5-0.160X6+0.725X7-0.354X8+0.077X9+0.032X10+0.139X11-0.004X12-0.012X13-0.089X14-0.828X15+0.255X16-0.011X17+0.052X18+0.215X19

(6)

由式(3)～式(6)可知，主成分A主要受到人口密度、非農(nóng)人口、總?cè)丝凇⒐房瓦\量的正向影響，受到日照時數(shù)和降水量的負向影響；主成分B主要受到日照時數(shù)的正向影響，受到人口密度較大負向影響；主成分C主要受到相對濕度、全社會固定資產(chǎn)投資總額、財政收入的正向影響，受到人口密度、公路客運量的負向影響；主成分D主要受到農(nóng)業(yè)人口、平均氣溫的正向影響，而受到公路客運量的負向影響顯著。

主成分A是一個由城鎮(zhèn)化發(fā)展程度來衡量的指標，其中人口密度、非農(nóng)人口有著顯著的正向推動作用，人口向城市聚集，促進了公路客運量的上升以及非農(nóng)人口的增加，因此，在主成分A的驅(qū)動下麗江壩區(qū)的植被經(jīng)歷了由農(nóng)業(yè)植被不斷向城市人工綠化發(fā)展的變化過程[30]。主成分B受到平均氣溫、日照時數(shù)的正向影響十分突出，由此看出主成分B是由自然因素主導(dǎo)的驅(qū)動力。主成分C中相對濕度、全社會資產(chǎn)投資總額的正向影響明顯，人口密度和公路客運量負向影響突出，是人口、社會經(jīng)濟驅(qū)動的綜合指標[31-32]。在C的驅(qū)動下，壩區(qū)植被從受到自然因素影響逐步成為受社會經(jīng)濟干擾的人工景觀。平均氣溫、農(nóng)業(yè)人口對主成分D有著較強的正向推動作用，由此可見在農(nóng)業(yè)的驅(qū)動下植被成為一個受到人工勞作干擾的半自然景觀。

2.3.2 植被覆蓋度變化的驅(qū)動分析 經(jīng)過偏最小二乘回歸計算，得到分析結(jié)果表達式(7)

Y=0.097X1+0.223X2+0.248X3+0.463X4-0.089X5-0.273X6+0.315X7-0.502X8+0.066X9+0.079X10+0.082X11+0.200X12+0.031X13-0.093X14-0.518X15+0.258X16+0.037X17+0.056X18+0.188X19

(7)

由式(7)可知，人口、社會經(jīng)濟因素對麗江壩區(qū)植被覆蓋度變化的影響程度要高于自然因素，其中農(nóng)業(yè)人口、全社會固定資產(chǎn)投資總額、財政收入對植被覆蓋度變化起到較為顯著的正向影響，人口密度、公路客運量、非農(nóng)人口對變化的負向影響突出。

研究期內(nèi)麗江產(chǎn)業(yè)結(jié)構(gòu)發(fā)生了深刻變化，三大產(chǎn)業(yè)由1990年的“一三二”的結(jié)構(gòu)轉(zhuǎn)變?yōu)?018年的“三二一”[33]，2018年全市經(jīng)濟增長主要由第二、三產(chǎn)業(yè)支撐，第三產(chǎn)業(yè)的貢獻率達到32.17%，其中旅游業(yè)和服務(wù)業(yè)都呈現(xiàn)逐年增長的趨勢。同時人口密度與公路客運量的增長加快了城鎮(zhèn)化的進程，因此快速城鎮(zhèn)化、社會經(jīng)濟發(fā)展成為引發(fā)植被覆蓋度變化的重要因素。

2.3.1.1 自然因子對植被覆蓋度的影響 降水量、平均氣溫、相對濕度、日照時數(shù)都是植被生長的重要決定因素[34]。由表6可知，自然因子對植被覆蓋度的變化都有正向的影響，其中日照時數(shù)的正向作用與VIP值均最大，降水量的影響遠小于其他3個因子。麗江市從2009年末開始經(jīng)歷了3年的連續(xù)干旱，較少的降雨量導(dǎo)致壩區(qū)內(nèi)黑龍?zhí)度和耆珨嗔?，各景區(qū)、公園用水緊張，農(nóng)田灌溉困難，因此降水量在植被生長、植被覆蓋度的增加中起到正向影響作用有限。

表6 自然因子與植被覆蓋度的偏最小二乘回歸分析結(jié)果

2.3.1.2 人口因子對植被覆蓋度的影響 4個因子中只有農(nóng)業(yè)人口對植被覆蓋度起到了正向影響，總?cè)丝?、非農(nóng)人口、人口密度都是起到負向的影響作用，其中人口密度負向作用最大，VIP值最高，說明人口密度對植被覆蓋度的解釋力度最大(表7)。人口的增長加速了城鎮(zhèn)化進程，2007、2013、2018年麗江城鎮(zhèn)化率分別為25%、32.63%、40.44%。人口密度的上升需要通過增加建筑密度、擴大城市面積等來滿足居民不斷提升的居住、生活、娛樂的需求，植被的生長空間便受到擠壓，因此人口聚集、城鎮(zhèn)化加快與植被覆蓋的減少有著較高的相關(guān)性。

表7 人口因子與植被覆蓋度的偏最小二乘回歸分析結(jié)果

2.3.1.3 經(jīng)濟因子對植被覆蓋度的影響 由表8中可知，經(jīng)濟因子對植被覆蓋度的變化都起到一定作用(VIP>0.5)，財政收入、全社會固定資產(chǎn)投資總額、旅游收入的正向影響較大，公路客運量對植被覆蓋度的負向影響最為突出。

表8 經(jīng)濟因子與植被覆蓋度的偏最小二乘回歸分析結(jié)果

麗江蓬勃發(fā)展的旅游業(yè)，以及地方財政收入與全民固定資產(chǎn)投資額的上漲，都讓城市綠化與生態(tài)環(huán)境建設(shè)得到了更多資助，城市綠化環(huán)境不斷得到改善。市區(qū)內(nèi)的綠化覆蓋率從2006年的27%提升到2008年的32%，人均公共綠地面積由10 m2增至15 m2。隨著綠化建設(shè)的不斷加強，2018年市區(qū)綠地面積達到559.34萬m2，人均公園綠地面積17.5 m2，綠化覆蓋率達到34.57%[35]。

城市主干道香格里拉大道及其延伸段分別在1997、1999年建成通車，2013年大理麗江高速路通車，暢通的交通帶來了更多的游客和務(wù)工人員，麗江游客量由1987年5 000人次升至2018年的4 643.3萬人次。逐年上漲的游客量激發(fā)了運輸業(yè)、餐飲業(yè)、酒店業(yè)、建筑業(yè)以及地產(chǎn)業(yè)的繁榮[36-38]，主城區(qū)面積快速增加，其周邊大量耕地被征用。因此，在上升的客運量、快速城鎮(zhèn)化的負向影響下植被快速消減，同時建設(shè)活動頻繁、車流量增加等因素加劇了環(huán)境的污染，植被生長也受到影響。

3 結(jié)論與討論

3.1 結(jié)論

1987-2018年，麗江壩區(qū)中有占總面積24.46%地區(qū)的植被覆蓋度出現(xiàn)了明顯退化現(xiàn)象，集中于主城區(qū)與玉龍新城內(nèi)以及村落周邊。人為干擾程度的強弱對植被覆蓋度的變化有著直接的影響，人口流動以及社會經(jīng)濟活動越頻繁的地區(qū)植被退化程度就越顯著，反之遠離城市的壩區(qū)邊緣的植被得到逐步恢復(fù)。

在研究時段內(nèi)，主城區(qū)內(nèi)植被覆蓋地區(qū)的破碎化程度高于郊區(qū)。1987-2018年，主城區(qū)內(nèi)中低、中等植被覆蓋度區(qū)域的斑塊形狀逐漸復(fù)雜，中高、高覆蓋地區(qū)的斑塊形狀越發(fā)規(guī)則；主城區(qū)內(nèi)中低、中等植被覆蓋度地區(qū)的景觀連接度增加，而中高、高覆蓋度地區(qū)降低；主城區(qū)及其周邊地區(qū)的中低、中等植被覆蓋度的斑塊聚集度高，中高、高覆蓋地區(qū)的聚集度逐漸降低。日益顯著的破碎化、斑塊形狀的變化以及植被的連接度的降低都與城鎮(zhèn)化建設(shè)活動有較強的關(guān)聯(lián)性。

麗江產(chǎn)業(yè)結(jié)構(gòu)的變更成為壩區(qū)植被覆蓋度變化的誘因，人口、社會經(jīng)濟是導(dǎo)致變化的主要驅(qū)動力，其中人口密度、公路客運量起到突出的負向作用，投影重要性分別為1.172、1.198。自然因子中日照時數(shù)對植被覆蓋度的變化有著顯著的正向作用。

3.2 討論

已有研究表明[39-41]，由于在人口、經(jīng)濟的驅(qū)動作用下城鎮(zhèn)化進程加速，會直接引起植被的空間分布格局的變化。本研究運用同心圓空間梯度模型對快速城鎮(zhèn)化的旅游地區(qū)植被覆蓋度的空間變化特征進行了詳細描述。然而，麗江壩區(qū)整體為南北長、東西短的不規(guī)則帶狀形態(tài)，在采用同心圓空間梯度分析中存在著每個環(huán)帶面積不均、分析結(jié)果不能明確指出植被覆蓋度變化的具體方位等問題，之后的研究將結(jié)合多方向采樣模式進行綜合分析[42]，為植被空間分布特征描述提供更為準確的量化依據(jù)。植被變化是一個較為復(fù)雜的過程，本研究中僅選擇了類型水平中的景觀指數(shù)分析了植被的分布特征，今后的研究可以從景觀水平上選擇多樣性指數(shù)進一步探討人為活動干擾下植被的多樣性變化特征[43]。本研究構(gòu)建了麗江壩區(qū)植被變化的驅(qū)動力模型，發(fā)現(xiàn)旅游經(jīng)濟影響下的城鎮(zhèn)化驅(qū)動作用最為突出，但影響植被變化的驅(qū)動因子眾多[44-45]，且因子間的關(guān)聯(lián)程度也較為復(fù)雜，在驅(qū)動力因子的選擇中不可避免地會受到主觀經(jīng)驗的干擾，從而產(chǎn)生信息遺漏或冗余。今后的研究可以采用機器學(xué)習(xí)方法篩選驅(qū)動因子和回歸模型[46]，為研究分析提供客觀依據(jù)。