羅紫薇, 黎奕暉, 胡希軍,2, 杜心宇, 韋寶婧, 李芮芝, 曹詩怡,3, 胡 益

(1.中南林業(yè)科技大學(xué) 風(fēng)景園林學(xué)院, 長沙 410004; 2.中南林業(yè)科技大學(xué) 城鄉(xiāng)景觀生態(tài)研究所, 長沙 410004;3.湖南理工學(xué)院土木建筑工程學(xué)院, 湖南 岳陽 414006; 4.國家林業(yè)和草原局 中南調(diào)查規(guī)劃設(shè)計(jì)院, 長沙 410014)

加大自然保護(hù)區(qū)保護(hù)力度、改進(jìn)保護(hù)方式、維護(hù)森林生態(tài)系統(tǒng)穩(wěn)定性正是響應(yīng)生態(tài)文明建設(shè)目標(biāo)的重要舉措。近年來，湖南圍繞“健康森林、美麗濕地、綠色通道、秀美村莊”林業(yè)發(fā)展主題，強(qiáng)化森林資源保護(hù)管理，森林資源呈現(xiàn)持續(xù)穩(wěn)步增長態(tài)勢。金童山國家級自然保護(hù)區(qū)是2016年獲批湖南省第一個(gè)國家公園體制試點(diǎn)的重要片區(qū)之一，要回答如何實(shí)現(xiàn)人與自然的和諧共處這一根本問題，需要對保護(hù)區(qū)進(jìn)行全方位、多層次的數(shù)據(jù)收集和相關(guān)分析，為自然保護(hù)區(qū)管理體系的建設(shè)提供科學(xué)依據(jù)。

性質(zhì)、大小及形狀各異的森林景觀要素在空間上的分布形成了獨(dú)有的森林景觀格局[1]。景觀格局研究多從景觀結(jié)構(gòu)[2]、景觀類型轉(zhuǎn)移、景觀格局特征[3-5]、景觀格局模擬[6]和影響因子等[7]角度入手，從景觀格局現(xiàn)狀、內(nèi)部變化、未來預(yù)測和驅(qū)動力等方面進(jìn)行分析。針對景觀格局影響因子，現(xiàn)有研究中已有較多分析方法，例如：主成分分析[8]、多元線性回歸分析[9]、地理加權(quán)回歸分析[10]、除趨勢典范對應(yīng)分析等[11]對影響因子進(jìn)行定性或定量化研究。近年來，地理探測器被廣泛用于探測空間分異性、揭示分異背后驅(qū)動因素。地理探測器包括4個(gè)方面的內(nèi)容——分異及因子探測器、交互作用探測器、風(fēng)險(xiǎn)區(qū)探測器和生態(tài)探測器，其中，分異及因子探測器中的q值可代表各因子對整體空間分異的解釋程度，即探測全局驅(qū)動力和影響因子。地理探測器在分析區(qū)域不同類型空間分布影響因子上的功能強(qiáng)大，適用于進(jìn)行森林景觀格局影響因子方面的研究，在植被覆蓋及空間分布方面，王治國等[12]通過地理探測器計(jì)算q值分析自然因素、人文因素對植被覆蓋的影響，研究結(jié)果為當(dāng)?shù)馗纳浦脖桓采w環(huán)境提供理論依據(jù)；吳家榮等[13]運(yùn)用因子探測器識別影響丹霞梧桐分布的主要環(huán)境驅(qū)動因子；宋媛等[14]通過地理探測器從自然因素和人為因素兩個(gè)方面，識別甘肅省生態(tài)環(huán)境質(zhì)量影響因素強(qiáng)度。

本研究以金童山國家級自然保護(hù)區(qū)金童山片區(qū)為湖南森林景觀代表，通過對2005年、2009年、2015年三期遙感數(shù)據(jù)進(jìn)行景觀分類提取，通過景觀格局總體演變和景觀類型轉(zhuǎn)移矩陣分析2005—2015年金童山國家級自然保護(hù)區(qū)金童山片區(qū)的森林景觀格局演變狀況，并通過因子探測器定性分析影響森林景觀格局形成的內(nèi)外部影響因子、空間定量化給出指標(biāo)權(quán)重，并對影響因子權(quán)重與景觀類型轉(zhuǎn)移量之間的關(guān)系進(jìn)行解釋，可充分掌握保護(hù)區(qū)森林景觀的分布規(guī)律，為保護(hù)區(qū)生態(tài)調(diào)控、低干擾下的適度開發(fā)提供借鑒方法和科學(xué)素材。

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

1.1 研究區(qū)概況

湖南氣候、土壤適宜樹木生長，是我國南方重點(diǎn)林區(qū)省，森林資源豐富。根據(jù)第九次全國森林資源清查結(jié)果，湖南省森林面積為1 052.58萬hm2；森林蓄積量達(dá)40 715.73萬m3，森林覆蓋率達(dá)49.69%[15]，主要分布湘西、湘南、湘東的山地區(qū)域。湖南金童山國家級自然保護(hù)區(qū)位于湖南省西南邊陲的城步苗族自治縣境內(nèi)，地理坐標(biāo)為東經(jīng)110°07′48″—110°33′36″，北緯26°07′39″—26°20′06″，包括明竹老山片區(qū)和金童山片區(qū)，本研究僅對金童山片區(qū)進(jìn)行研究(圖1)。研究區(qū)域總面積為10 293 hm2，屬中山深谷地貌，氣候類型為中亞熱帶季風(fēng)濕潤氣候區(qū)山地氣候，年平均氣溫16.1℃，年降水量1 218.5 mm，全年日照時(shí)數(shù)在1 134.6～1 601.5 h，最高海拔1 784.1 m。保護(hù)區(qū)水平地帶性土壤為紅壤，母巖以花崗巖、砂巖等為主。保護(hù)區(qū)植物區(qū)系起源古老，屬華中、華南區(qū)系的過渡區(qū)域，具有典型的南嶺山地區(qū)域特征。

圖1 研究區(qū)地理區(qū)位

保護(hù)區(qū)建立歷程從20世紀(jì)八十年代開始，1981年，民族老山自然保護(hù)區(qū)成立；2009年，省級自然保護(hù)區(qū)建立；2013年，湖南金童山國家級自然保護(hù)區(qū)建立，并設(shè)立了湖南金童山國家級自然保護(hù)區(qū)管理處，對所轄的森林資源進(jìn)行保護(hù)和管理，可以說，近50 a來湖南金童山國家級自然保護(hù)區(qū)處于一定的保護(hù)狀態(tài)，人為干擾受到限制。

1.2 景觀分類體系構(gòu)建

根據(jù)研究目的，同時(shí)滿足遙感影像監(jiān)督分類條件，將森林景觀類型主要分為針葉林、常綠闊葉林、落葉闊葉林、常綠—落葉闊葉林、竹林、灌木林地、耕地、牧草地共8類?？紤]到研究區(qū)河流的形狀均為狹窄線狀，建設(shè)用地的類型均為房屋或道路，且被耕地所圍繞，因此，本研究遙感影像監(jiān)督分類時(shí)，綜合考慮影像分辨率、現(xiàn)狀特征及樹木的遮擋作用，將河流、道路、房屋三者信息納入耕地信息進(jìn)行監(jiān)督分類[16]，在完成遙感影像監(jiān)督分類后，結(jié)合相關(guān)數(shù)據(jù)，通過目視解譯將河流、道路、房屋單獨(dú)劃出，此3類數(shù)據(jù)在2005—2015年間均無明顯變化，因此不作為森林景觀類型參與相應(yīng)分析。

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

1.3.1 數(shù)據(jù)源 研究數(shù)據(jù)主要有研究區(qū)1∶10 000地形圖、森林資源2類調(diào)查數(shù)據(jù)(2015年)，2005年、2009年、2015年三期遙感數(shù)據(jù)。本研究選取了三期分辨率為30 m×30 m的Landsat遙感影像為數(shù)據(jù)源(地理空間數(shù)據(jù)云http:∥www.gscloud.cn/)。研究區(qū)內(nèi)用地多為林地，在選取成像時(shí)間時(shí)，多選取植被生長初期或末期質(zhì)量較好的月份。選取了數(shù)據(jù)類型級別為L1T、成像時(shí)間為2005年9月9日、2009年10月6日的TM影像以及2015年4月14日的OLI_TIRS影像，三期影像數(shù)據(jù)云覆蓋率分別為1%，0%，0.45%。其中，2005年9月9日數(shù)據(jù)存在有云區(qū)域，為使得3期數(shù)據(jù)研究區(qū)范圍內(nèi)云覆蓋率為0，用同年10月11日的數(shù)據(jù)進(jìn)行替換。

1.3.2 數(shù)據(jù)預(yù)處理 L1T數(shù)據(jù)產(chǎn)品已經(jīng)過輻射定標(biāo)，故數(shù)據(jù)預(yù)處理只需進(jìn)行大氣校正、正射校正、圖像融合和裁剪等步驟，空間投影坐標(biāo)系統(tǒng)使用WGS84UTM橫軸墨卡托投影。綜合實(shí)地調(diào)研采樣數(shù)據(jù)、Google Earth衛(wèi)星影像及森林資源二調(diào)數(shù)據(jù)，在ENVI中通過目視判讀建立訓(xùn)練樣本和驗(yàn)證樣本，并采用最大似然監(jiān)督分類法對三期影像數(shù)據(jù)進(jìn)行分類提取和結(jié)果驗(yàn)證，結(jié)果顯示Kappa系數(shù)分別為0.90，0.93，0.94，滿足分類提取要求。為提高分類準(zhǔn)確度，將分類結(jié)果進(jìn)行分類后處理與目視修改。將ENVI監(jiān)督分類數(shù)據(jù)導(dǎo)出為Shapefile格式，在Arc GIS軟件中與河流、道路、房屋Shapefile數(shù)據(jù)進(jìn)行聯(lián)合，得到研究區(qū)景觀類型分布圖，并將各景觀類型進(jìn)行編碼，轉(zhuǎn)為精度為10 m×10 m的柵格圖。

1.4 因子探測器

地理探測器(Geographical Detector)是探測空間分異性，以及揭示其背后驅(qū)動力的一組統(tǒng)計(jì)學(xué)方法，已在自然科學(xué)、社會科學(xué)、環(huán)境科學(xué)和人類健康等多領(lǐng)域應(yīng)用[17]。本研究采用地理探測器中的因子探測器，利用q值解釋各影響因子對景觀格局的影響程度。公式如下：

(1)

將所有自變量的q值歸一化處理后，可獲得影響景觀格局形成的各影響因子權(quán)重值。具體步驟為：(1)格點(diǎn)密度確定，空間離散點(diǎn)準(zhǔn)備。已知三期數(shù)據(jù)各景觀類型內(nèi)最小斑塊面積均大于10 000 m2，綜合考慮地理探測器計(jì)算精度和最大計(jì)算量，將研究范圍轉(zhuǎn)換為柵格大小100 m×100 m的柵格數(shù)據(jù)，即可保證在柵格轉(zhuǎn)點(diǎn)后，每個(gè)景觀斑塊都能在點(diǎn)數(shù)據(jù)集中表現(xiàn)。(2)景觀數(shù)據(jù)處理。將三期景觀類型分布數(shù)據(jù)重分類，用1～8代表8類景觀類型；提取各年各景觀類型點(diǎn)數(shù)據(jù)，在ArcGIS中進(jìn)行核密度分析，并采用自然間斷點(diǎn)法將進(jìn)行重分類分級。(3)影響指標(biāo)處理。一個(gè)區(qū)域內(nèi)景觀格局的形成往往受到自然地理因素和人為影響因素兩個(gè)方面的影響。本研究選取高程、坡度、坡向、水系、土壤以及交通因子(林區(qū)干道)、旅游因子(游步道)、村民生活干擾因子等8個(gè)影響指標(biāo)，水系、林區(qū)干道、游步道對周圍景觀格局的影響方式通過研究范圍中的點(diǎn)與其間的距離體現(xiàn)(表1)。在ArcGIS中采用以下劃分方法將上述影響指標(biāo)劃分為9個(gè)等級(僅土壤采用4分類)進(jìn)行離散化處理。

表1 影響指標(biāo)處理方法

(4)數(shù)據(jù)集準(zhǔn)備。在Arc GIS中將上述各年的景觀類型分布數(shù)據(jù)、各景觀類型密度分級數(shù)據(jù)和影響指標(biāo)數(shù)據(jù)采用多值提取指點(diǎn)的方式提取至整體網(wǎng)格點(diǎn)；將各屬性表分別轉(zhuǎn)出至Excel參與地理探測器運(yùn)算。將因子探測器結(jié)果歸一化處理。

2 結(jié)果與分析

2.1 景觀格局演變分析

本研究從森林景觀格局的外部總體特征、內(nèi)部景觀類型轉(zhuǎn)移兩個(gè)方面，利用統(tǒng)計(jì)分析、轉(zhuǎn)移矩陣分析法對景觀格局演變進(jìn)行多角度分析。

2.1.1 景觀格局演變總體分析 研究區(qū)總面積為10 288.94 hm2，由2005年、2009年、2015年研究區(qū)森林景觀類型分布情況及各景觀類型面積統(tǒng)計(jì)表可以得知(圖2，表2)：(1)從景觀類型的平面分布來看，常綠闊葉林主要分布于研究區(qū)北部；落葉闊葉林主要分布于西部和南部；耕地主要分布于道路兩旁；針葉林主要為人工種植的杉木林，故而分布比較集中，主要分布于保護(hù)區(qū)中部高程較低的地帶；牧草地主要分布于高海拔區(qū)域，且覆蓋范圍較小；灌木林與竹林的分布較為零散；(2)從景觀類型的規(guī)模來看，研究區(qū)90%以上區(qū)域均為森林景觀，以常綠闊葉林、落葉闊葉林為其景觀基質(zhì)，占總面積超60%，面積占比較大的景觀類型往往分布較為集中；其次為針葉林，牧草地占比最??；(3)從各景觀類型面積變化情況來看，常綠闊葉林面積整體呈現(xiàn)減少狀態(tài)，落葉闊葉林取代常綠闊葉林成為研究區(qū)最主要的景觀類型；灌木林地和竹林呈現(xiàn)增長趨勢，其中灌木林地增長幅度較大；針葉林、常綠落葉闊葉林、牧草地占比變化幅度較小。

表2 研究區(qū)各年景觀類型面積

圖2 各年景觀類型分布

2.1.2 景觀類型轉(zhuǎn)移 在Arc GIS中利用面積制表工具，對2005年、2009年、2015年三期數(shù)據(jù)的景觀類型面積轉(zhuǎn)化情況進(jìn)行統(tǒng)計(jì)，分別得到2005—2009年、2009—2015年研究區(qū)域景觀類型轉(zhuǎn)移矩陣、并對轉(zhuǎn)入、轉(zhuǎn)出量進(jìn)行計(jì)算(表3—4)。

由表3—4可知、2005—2009年景觀轉(zhuǎn)移面積總量為4 116.95 hm2，占保護(hù)區(qū)總面積的40.01%，保護(hù)區(qū)景觀類型動態(tài)變化較為頻繁；2009—2015年景觀轉(zhuǎn)移面積總量為3 700.58 hm2，占比35.97%，動態(tài)變化情況減緩。落葉闊葉林在2009年、2015年有80.04%，70%的部分未發(fā)生景觀類型轉(zhuǎn)變，最為穩(wěn)定，而常綠闊葉林在2005—2009時(shí)期段的未轉(zhuǎn)變量為48.01%，主要向落葉闊葉林、針葉林產(chǎn)生大量轉(zhuǎn)移，2009—2015年為70.13%，其變化逐漸趨于穩(wěn)定；耕地、針葉林的穩(wěn)定性也較為良好。竹林、灌木林地、牧草地等景觀類型由于斑塊面積較小，發(fā)生轉(zhuǎn)移情況普遍較為頻繁。耕地和牧草地主要與灌木林地之間發(fā)生相互轉(zhuǎn)換。

表3 2005-2009年研究區(qū)域景觀類型轉(zhuǎn)移矩陣 hm2

表4 2009-2015年研究區(qū)域景觀類型轉(zhuǎn)移矩陣 hm2

2.2 影響因子分析

本研究主要從自然地理因素和人為影響因素兩個(gè)方面進(jìn)行定量分析，從社會及政策角度進(jìn)行原因分析及補(bǔ)充。遵循區(qū)域特殊性和數(shù)據(jù)可獲取性等原則，為反映10 a間影響森林景觀格局的因子權(quán)重變化，選擇高程、坡度、坡向、水系、土壤等5類自然地理因素因子；鑒于研究區(qū)的用地性質(zhì)為國家級自然保護(hù)區(qū)，選用交通因子(林區(qū)干道)、旅游因子(游步道)、村民生活干擾因子等3類人為影響因素。

影響因子分析從景觀級別和景觀類型級別兩個(gè)方面分別進(jìn)行描述。

2.2.1 景觀級別影響因子分析 為探究各影響因子對整體景觀格局影響的權(quán)重大小，在地理探測器中，分別將三期整體景觀類型分布作為因變量，當(dāng)年8類影響因子作為自變量進(jìn)行運(yùn)算，q值歸一化后的結(jié)果(圖3)顯示，各指標(biāo)因素的影響力排序?yàn)?005年：村民生活干擾>高程>坡向>土壤>林區(qū)干道>坡度>水系>游步道，2009年：高程>坡向>土壤>坡度>林區(qū)干道>村民生活干擾>游步道>水系，2015年：高程>坡度>村民生活干擾>土壤>坡向>林區(qū)干道>水系>游步道。

圖3 2005-2015年整體景觀格局影響因子權(quán)重分布

(1)因子屬性上，自然地理因素與人為影響因素均對區(qū)域景觀格局形成產(chǎn)生重要影響。自然地理因素占優(yōu)勢，三年占比均超65%，這是由于研究區(qū)以山地為主，內(nèi)部環(huán)境復(fù)雜，小范圍的立體小氣候特點(diǎn)，結(jié)合金童山較為清晰的土壤垂直分布規(guī)律，影響著植被分布，進(jìn)而形成多變的景觀格局。人為影響因素方面，由于村民生活干擾因子在2005年影響強(qiáng)度最大，2009年強(qiáng)度最小，其他因子在10 a間變化程度不明顯，總體人為影響強(qiáng)度變化呈相同規(guī)律；(2)總體權(quán)重占比上，2005—2015年，自然地理因素中高程影響各年整體景觀格局的力度始終最大，水系影響最小，3類人為因素中村民生活干擾的影響力度最大；(3)權(quán)重時(shí)間變化上，高程、坡向、土壤因素對景觀格局的影響權(quán)重呈波動趨勢，坡度的影響權(quán)重呈上升趨勢；由于2009年湖南城步金童山省級自然保護(hù)區(qū)成立，隨后《關(guān)于設(shè)置城步苗族自治縣金童山自然保護(hù)區(qū)管理處的批復(fù)》、《湖南城步金童山省級自然保護(hù)區(qū)管理辦法》等政策的發(fā)布，保護(hù)區(qū)林區(qū)干道、村民生活干擾影響受到有效控制，影響權(quán)重減小。

2.2.2 景觀類型級別影響因子分析 將各影響因子作為自變量，分別將三年各景觀類型分布作為因變量，載入地理探測器運(yùn)算，獲得影響各景觀類型分布的影響因子權(quán)重大小及三期數(shù)據(jù)的變化情況(圖4—5)。

圖4 各景觀類型兩類影響因素權(quán)重占比

由圖4可知，自然、人文兩類因素對不同森林景觀類型分布的影響力度存在差異。在8類景觀類型中，針葉林、常綠闊葉林、落葉闊葉林和常綠、落葉闊葉林受自然因素影響較大，各年權(quán)重均在0.6以上，由景觀分布圖可知，此類景觀類型斑塊面積較大，意味著同一類型景觀覆蓋面廣、分布集中，人為影響程度有限，因此，受到地形、地貌及小氣候等自然地理因素的影響較大；人為因素對耕地分布的影響明顯，這是由于研究區(qū)內(nèi)道路和村莊等人為影響因子與耕地均主要分布在中、低海拔區(qū)域，游客、農(nóng)民活動較多，同時(shí)，金童山片區(qū)的試驗(yàn)區(qū)毗鄰南山牧場，政府鼓勵(lì)周邊村民利用山林、草場資源積極發(fā)展林下種植和養(yǎng)殖業(yè)，轉(zhuǎn)換經(jīng)濟(jì)發(fā)展模式，減少對森林資源的依賴，因此在10 a間耕地面積發(fā)生增長，故耕地對兩類影響因素的響應(yīng)程度相當(dāng)，甚至在2009年人為影響因素的影響強(qiáng)度超過了50%；同樣由于一系列政策、措施的實(shí)行，竹林、灌木林地的兩類影響因素權(quán)重波動較大。

由圖5可知，針葉林、常綠闊葉林、落葉闊葉林分布受高程影響最大，統(tǒng)計(jì)各景觀點(diǎn)數(shù)據(jù)中高程分布，可見此3類景觀分別集中分布在海拔1 000～1 400 m，800～1 200 m，1 000～1 600 m，其他高程地帶少有分布，常綠闊葉林、落葉闊葉林同時(shí)受土壤分布的影響；常綠、落葉闊葉林主要受到高程、土壤和游步道的影響，竹林對高程、游步道的影響反映明顯，牧草地主要受到高程、土壤的較大影響；而對于中低海拔且斑塊較破碎的灌木林地和耕地而言，更多的受到高程、游步道、村民生活干擾的影響，與權(quán)重結(jié)果相對應(yīng)；坡度、坡向、水系對各景觀類型分布影響較小。

圖5 2005-2015年各景觀類型影響因子權(quán)重分布

2.2.3 影響因子權(quán)重與各類型轉(zhuǎn)移量關(guān)系 由上述影響因子分析可以對景觀格局演變分析中的景觀類型轉(zhuǎn)移進(jìn)行解釋：(1)針葉林、常綠闊葉林、落葉闊葉林是研究區(qū)內(nèi)的主要景觀類型，面積較大，地理探測器結(jié)果顯示其自然地理因素占比較大，面積轉(zhuǎn)移以自然演替為主；(2)常綠、落葉闊葉林因人為影響而發(fā)展受限，景觀斑塊大小和位置發(fā)生變動，被以落葉闊葉林為主的其他景觀類型替代，同時(shí)其他類型也經(jīng)過自然演替形成了常綠落葉闊葉林，因此其總面積呈增長趨勢；(3)竹林因其無性繁殖能力及蔓延生長特點(diǎn)，對研究區(qū)典型森林景觀進(jìn)行侵入，占據(jù)大量其他景觀類型面積，故面積轉(zhuǎn)移矩陣中出現(xiàn)針葉林、常綠闊葉林、落葉闊葉林向竹林轉(zhuǎn)移的情況；(4)一方面牧草地經(jīng)自然演替轉(zhuǎn)變?yōu)楣嗄玖值兀硪环矫婀嗄玖值厥茌^大的人為干擾影響而發(fā)生逆向演替形成了牧草地，2005—2009年由于人為干擾影響發(fā)生的速率大于自然演替，故在此間牧草地面積增加，隨著保護(hù)力度的加大，人為干擾減少，牧草地逐漸向更高層景觀類型轉(zhuǎn)變；(5)耕地景觀的分布主要受人為影響因素的干擾，耕地的增加主要是由于人為的開墾，而棄耕又會導(dǎo)致耕地向灌木林演替，兩個(gè)原因的作用范圍存在差異，導(dǎo)致了耕地面積的增加及其景觀格局的形成。

2.3 討 論

景觀格局是受到多種因素相互作用下的復(fù)合結(jié)果。本研究以湖南金童山國家級自然保護(hù)區(qū)金童山片區(qū)為湖南森林景觀代表，從自然地理、人為兩個(gè)角度提取了高程、坡度、坡向、水系、土壤、交通(林區(qū)干道)、旅游(游步道)、村民生活干擾因子等共計(jì)8個(gè)因子進(jìn)行其森林景觀格局形成的影響因素研究，其中水系、主路和次路的影響通過距離遠(yuǎn)近表示，3個(gè)距離數(shù)據(jù)是在TerrSet軟件下，運(yùn)用Distance工具獲得的平面距離分布圖，與山地環(huán)境中的實(shí)際距離存在差異，可能在一定程度上影響權(quán)重的結(jié)果。同時(shí)，本研究區(qū)復(fù)雜的小氣候效應(yīng)更加劇了全面揭示影響因子的難度，為尋求更為細(xì)致和本土化的景觀格局影響因子，未來可提高研究精度，考慮生物因素，對局部環(huán)境指標(biāo)進(jìn)行監(jiān)測。

本研究創(chuàng)新性地分析了影響因子權(quán)重與各類型轉(zhuǎn)移量之間的關(guān)系，但仍處于一個(gè)定性的認(rèn)識階段，今后可嘗試通過定量的方法將影響因子與轉(zhuǎn)移量之間的內(nèi)在關(guān)系進(jìn)行量化，剖析景觀類型發(fā)生轉(zhuǎn)移以及轉(zhuǎn)移量大小的深層原因。

針對景觀格局模擬，可利用現(xiàn)有數(shù)據(jù)，選擇模型對景觀格局進(jìn)行模擬預(yù)測，并考慮運(yùn)用地理探測器因子探測的方法對適宜性圖集中所用的影響因子圖集進(jìn)行權(quán)重修正，避免人為主觀判定，提高適宜性圖集的客觀性、科學(xué)性和本土適用性。同時(shí)，針對森林景觀格局的特殊性，可嘗試通過植物遙感反演結(jié)合地面實(shí)測的方法，運(yùn)用編程手段反演、計(jì)算植被的各項(xiàng)指標(biāo)，并對未來植被生長情況進(jìn)行預(yù)測，從而達(dá)到對未來整個(gè)景觀格局模擬預(yù)測、提出科學(xué)干預(yù)手段、保證森林景觀格局穩(wěn)定且健康的最終目的。

3 結(jié) 論

(1)湖南金童山自然保護(hù)區(qū)金童山片區(qū)90%以上區(qū)域均為森林景觀，以常綠闊葉林、落葉闊葉林為其景觀基質(zhì)，牧草地面積占比最?。怀＞G闊葉林面積減少量較大，2009年被落葉闊葉林超越成為研究區(qū)最主要的景觀類型。2005—2009年景觀類型動態(tài)變化較為頻繁，2009—2015年情況減緩。常綠闊葉林向落葉闊葉林、針葉林產(chǎn)生大量轉(zhuǎn)移，落葉闊葉林轉(zhuǎn)移率小，最為穩(wěn)定；竹林、灌木林地、牧草地等景觀類型轉(zhuǎn)移情況較為頻繁。

(2)整體上，自然地理因素是主要影響區(qū)域整體景觀格局形成的因素，除2005年村民生活干擾對森林景觀的影響力度最大外，2009年、2015年高程的影響力度始終最大，游步道及水系的影響最?。桓骶坝^類型上，人為因素對耕地的分布影響明顯，常綠闊葉林和落葉闊葉林內(nèi)部自然條件情況復(fù)雜，受自然地理因素影響較大，其中，受高程影響最大；常綠、落葉闊葉林主要受到高程、土壤和游步道的影響，竹林對高程、游步道的影響反映明顯，牧草地主要受到高程、土壤的較大影響；坡度、坡向、水系對各景觀類型分布影響較小。

(3)影響因子分析可對景觀類型轉(zhuǎn)移進(jìn)行解釋：常綠闊葉林、落葉闊葉林面積轉(zhuǎn)移以自然演替為主；常綠落葉闊葉林因人為影響導(dǎo)致景觀斑塊大小和位置發(fā)生變動，與其他景觀類型發(fā)生相互轉(zhuǎn)換；竹林對典型森林景觀進(jìn)行侵入，占據(jù)大量其他景觀類型面積；自然演替與人為影響的雙重作用導(dǎo)致牧草地與灌木林地、耕地與灌木林地間發(fā)生景觀類型轉(zhuǎn)變?？梢?，研究區(qū)森林景觀格局的形成是各影響因子對各景觀類型綜合作用的結(jié)果。