肖義發(fā),戴桂萍,王夢(mèng)君,黃杰,張小鵬,曾衛(wèi)
(1.國(guó)家林業(yè)和草原局西南調(diào)查規(guī)劃院,云南 昆明 650216;2.富寧馱娘江省級(jí)自然保護(hù)區(qū)管護(hù)局,云南 文山 663400)
景觀格局是指斑塊和其他組成單元的類(lèi)型、數(shù)目,以及空間分布與配置等,反映了各種生態(tài)過(guò)程在不同尺度上作用的結(jié)果[1-2]。景觀格局分析可以量化、分析景觀要素的結(jié)構(gòu)特征及其相互之間的空間分布關(guān)系,從而客觀地分析景觀要素類(lèi)型的分布規(guī)律、形成機(jī)制,為進(jìn)一步研究景觀功能和動(dòng)態(tài)提供基礎(chǔ)資料[3-5],對(duì)研究區(qū)域景觀結(jié)構(gòu)組成、功能及動(dòng)態(tài)等具有重要意義[6-8]。隨著國(guó)內(nèi)外對(duì)景觀格局研究成果的不斷豐富,目前景觀格局分析方法大致可以分為景觀格局指數(shù)法、空間統(tǒng)計(jì)學(xué)法、景觀格局模型與模擬法。從種群或森林斑塊組成結(jié)構(gòu)、異質(zhì)性分析、空間結(jié)構(gòu)等方面對(duì)景觀格局進(jìn)行研究[9-14],不僅可以根據(jù)幾期數(shù)據(jù)進(jìn)行景觀格局動(dòng)態(tài)分析,也可以根據(jù)參數(shù)建立模型,模擬未來(lái)景觀格局等[15-16]。森林是陸地生態(tài)系統(tǒng)的主體,森林植被景觀格局及其變化,如斑塊的大小和形狀會(huì)影響種群的生存能力和抗干擾能力等,不僅對(duì)森林生態(tài)系統(tǒng)內(nèi)的物質(zhì)循環(huán)和能量流動(dòng)產(chǎn)生影響,而且制約生態(tài)過(guò)程,乃至影響森林的演替[17]。
富寧馱娘江省級(jí)自然保護(hù)區(qū)(以下簡(jiǎn)稱(chēng)“保護(hù)區(qū)”)位于滇東高原的東南邊緣,是中國(guó)生物多樣性?xún)?yōu)先保護(hù)區(qū)域——桂西黔南石灰?guī)r生物多樣性保護(hù)優(yōu)先區(qū)和世界生物多樣性中心之一的北部灣中心,具有重要的保護(hù)與研究?jī)r(jià)值。目前,對(duì)保護(hù)區(qū)的研究主要集中在保護(hù)管理、動(dòng)物、植物、水文等方面[18-20],對(duì)景觀格局方面的研究尚處于空白階段。因此,開(kāi)展保護(hù)區(qū)景觀分布的研究,對(duì)認(rèn)識(shí)區(qū)域景觀類(lèi)型組成、分布特征及規(guī)律具有重要意義,可以為區(qū)域生物多樣性保護(hù)及自然資源管理提供參考。
保護(hù)區(qū)位于云南省東南部文山壯族苗族自治州富寧縣境內(nèi),地理坐標(biāo)為105°37′~106°07′E、23°35′~24°27′N(xiāo),包括剝隘片區(qū)、干南片區(qū)、鳥(niǎo)王山片區(qū)、谷拉河片區(qū),總面積19 711 hm2,詳見(jiàn)圖1。該區(qū)域?qū)倌蟻啛釒Ц咴撅L(fēng)氣候,年平均氣溫19.3 ℃,年均降水量1 194.5 mm,土壤類(lèi)型有磚紅壤、赤紅壤、紅壤、石灰土、紫色土等。
圖1 研究區(qū)地理位置
保護(hù)區(qū)主要保護(hù)對(duì)象有石灰山季雨林、河谷落葉季雨林、半常綠季雨林、季風(fēng)常綠闊葉林、落葉闊葉林等森林生態(tài)系統(tǒng),以及馱娘江濕地生態(tài)系統(tǒng)和珍稀、瀕危野生動(dòng)植物及其棲息地,主要功能是維護(hù)滇東生物多樣性穩(wěn)定,保護(hù)滇東獨(dú)特的自然資源和自然環(huán)境,使其免受人為干擾及破壞。區(qū)內(nèi)主要植被類(lèi)型包括石灰山季雨林、落葉季雨林、暖溫性針葉林、熱性稀樹(shù)灌木草叢、人工植被等;保護(hù)動(dòng)植物有德保蘇鐵(Cycas debaoensis)、叉孢蘇鐵(Cycas segmentifida)、 云豹(Neofelis nebulosa)、 金貓(Catopuma temmincki)、大靈貓(Viverra zibetha)等。
研究數(shù)據(jù)主要來(lái)源于馱娘江自然保護(hù)植被圖、1∶50 000 地形圖、2022 年Landsat8 30 m 分辨率的遙感影像和實(shí)地調(diào)查數(shù)據(jù)。對(duì)地形圖進(jìn)行數(shù)字化,運(yùn)用ENVI5.3 對(duì)Landsat 數(shù)據(jù)進(jìn)行校正、增強(qiáng)及裁剪處理。在遙感影像預(yù)處理的基礎(chǔ)上,依據(jù)遙感影像的光譜特征,結(jié)合地面實(shí)際調(diào)查數(shù)據(jù),參照《中國(guó)植被》《云南植被》分類(lèi)的原則,以植被類(lèi)型為基礎(chǔ),結(jié)合非植被因子確定研究區(qū)景觀類(lèi)型[21],將保護(hù)區(qū)景觀類(lèi)型劃分為22 個(gè)景觀類(lèi)型,包括飛機(jī)草-五節(jié)芒高草地灌草叢、楓香林、桄榔灌叢、桄榔-羊蹄甲灌叢、含楓香的中草草叢、含楓香-栓皮櫟中草草叢、含羽葉楸-千張紙的高草草叢、中國(guó)無(wú)憂(yōu)花-任木林、櫟類(lèi)灌叢、羅浮錐-小花紅花荷林、麻櫟-栓皮櫟林、芒種花-毛蕨菜灌叢、毛蕨菜草地、木棉-蒙自合歡林、木棉-楹樹(shù)林、人面子-野芒果林、栓皮櫟林、五節(jié)芒-蕨高草地18 個(gè)自然植被景觀和八角林、云南松林2個(gè)人工林植被景觀,以及農(nóng)地、水域2 類(lèi)非植被類(lèi)型景觀,詳見(jiàn)圖2。
圖2 景觀類(lèi)型分布
在ArcGIS10.8.2 平臺(tái)下,基于DEM 數(shù)據(jù)得到的海拔、坡度、坡向等地形因子數(shù)據(jù),并采用疊加分析法,將景觀類(lèi)型圖層與地形因子進(jìn)行疊加分析。本研究以150 m 為距,將海拔因子劃分為300 m 及以下、301~450 m、451~600 m、601~750 m、751~900 m、901~1 050 m、1 051~1 200 m、1 201~1 350 m、1 351 m 及以上9 個(gè)等級(jí);坡度因子設(shè)置為5°以下、6°~15°、16°~25°、26°~35°、36°~45°、46°及以上6 個(gè)等級(jí);坡向因子分為平坡、東、東北、東南、南、西南、西、西北、北9 個(gè)坡向。
參照郭晉平[22]研究,將景觀類(lèi)型斑塊規(guī)模等級(jí)劃分為小斑塊、中斑塊、中大斑塊、大斑塊、超大斑塊、巨斑塊6 類(lèi),并采用ArcGIS 10.8.2 分級(jí)統(tǒng)計(jì)分析。
從表2 可以看出,保護(hù)區(qū)景觀類(lèi)型組成以人面子-野芒果林、含楓香的中草草叢、楓香林、中國(guó)無(wú)憂(yōu)花-任木林4 種自然植被景觀為主,它們分別占景觀總面積的10.84%、12.84%、12.94%、13.11%;八角林、云南松林等人工林植被景觀分別占景觀總面積的2.02%、0.62%;農(nóng)地、水域等非植被景觀分別占景觀總面積的15.74%、3.97%。可見(jiàn),保護(hù)區(qū)是以熱帶、亞熱帶植被為基質(zhì)的森林植被景觀為主,但各景觀類(lèi)型分布不均衡,存在較大差異。
表1 景觀要素斑塊規(guī)模等級(jí)
表2 各景觀類(lèi)型特征指數(shù)
從圖3 可以看出,保護(hù)區(qū)景觀斑塊規(guī)模主要以小斑塊為主,其占景觀斑塊總數(shù)的85.78%,其余斑塊占比均在7%以下。在自然植被景觀方面,除櫟類(lèi)灌叢、麻櫟-栓皮櫟林、栓皮櫟林外,其他景觀類(lèi)型小斑塊占比均大于或等于50%;八角林、云南松林等人工林植被景觀也是以小斑塊為主,面積占比達(dá)92%以上;非植被類(lèi)型景觀則以中小斑塊為主。
圖3 各景觀類(lèi)型斑塊等級(jí)斑塊數(shù)分布
從圖4 可以看出,各景觀類(lèi)型斑塊大小等級(jí)的分布面積與斑塊分布數(shù)量存在較大的差異。從各景觀類(lèi)型斑塊大小等級(jí)的分布面積來(lái)看,保護(hù)區(qū)巨斑塊面積占比最大,達(dá)61.18%,其余斑塊大小等級(jí)分布的面積比例在6%~10%。在自然植被景觀方面,主要以超大斑塊和巨斑塊為主,說(shuō)明該區(qū)域內(nèi)自然植被景觀破碎化程度低;人工林植被景觀主要以中大斑塊、中斑塊和小斑塊為主,這類(lèi)景觀類(lèi)型主要是受人為干擾所影響;農(nóng)地、水域等非植被景觀主要以巨斑塊為主,該類(lèi)型景觀分布較為集中。
圖4 各景觀類(lèi)型斑塊大小等級(jí)面積分布
3.3.1 景觀類(lèi)型隨海拔變化分析
保護(hù)區(qū)地處低、中山區(qū)域,區(qū)內(nèi)海拔以700 m以下為主,占總面積的74.43%。其中300 m 以下占總面積的11.37%,301~450 m 占總面積的23.81%,451~600 m 占總面積的22.34%,601~750 m 占總面積的16.91%。
通過(guò)景觀類(lèi)型圖層與海拔圖層疊加分析,得到不同海拔段各景觀類(lèi)型分布,結(jié)果見(jiàn)圖5。從圖5 可以看出,自然植被景觀在各海拔段均有分布,其中飛機(jī)草-五節(jié)芒高草地灌草叢、栓皮櫟林主要分布在1 351 m 及以上的高海拔區(qū)域;楓香林主要分布在301~900 m 海拔段;桄榔-羊蹄甲灌叢主要分布在301~750 m 海拔段;含楓香的中草草叢主要分布在900 m 及以下海拔段;含楓香-栓皮櫟中草草叢主要分布在601~1 200 m 海拔段;含羽葉楸-千張紙的高草草叢主要分布在300 m 及以下海拔段;中國(guó)無(wú)憂(yōu)花-任木林在各海拔段均有分布,但主要分布在1 050 m 及以下海拔段;櫟類(lèi)灌叢主要分面于301~1 050 m 海拔段;羅浮錐-小花紅花荷林、麻櫟-栓皮櫟林、芒種花-毛蕨菜灌叢主要分布于1 050 m及以上海拔段;毛蕨菜草地主要分布在901 m 及以上海拔段;木棉-蒙自合歡林主要分布于301~1 050 m 海拔段;木棉-楹樹(shù)林主要分布于301~450 m 海拔段;人面子-野芒果林主要分布于600 m 及以下海拔段;五節(jié)芒-蕨高草地主要分布在750 m 及以下海拔段。人工林植被景觀方面,八角林主要分布在451~900 m 海拔段,其他海拔段僅有少量的分布;云南松林各海拔段均有分布。非植被景觀方面,農(nóng)地主要分布在900 m 及以下海拔段;水域主要分布在450 m 及以下海拔段的河谷地。
圖5 不同海拔景觀類(lèi)型分布
3.3.2 景觀類(lèi)型隨坡度變化分析
保護(hù)區(qū)平地較少,坡度5°以下的面積僅占2.71%;坡度6°以上的面積占比98.29%,其中以16°~25°坡度為主,面積占比達(dá)41.39%。
由圖6 可以得知,自然植被景觀從坡度分布面積上來(lái)看雖在各坡度均有分布,但主要分布在16°及以上的區(qū)域;受地形地貌及土地利用空間的限制,區(qū)域內(nèi)人工林植被景觀主要分布于16°~35°坡度;農(nóng)地、水域等非植被類(lèi)型景觀各坡度均有分布,但主要分布于坡度35°及以下區(qū)域,其中農(nóng)地主要分布坡度為6°~25°,以梯田為主;河流景觀主要為馱娘江,主要分布坡度為6°~15°。
圖6 不同坡度景觀類(lèi)型分布
3.3.3 景觀類(lèi)型隨坡向變化分析
由圖7 可知,保護(hù)區(qū)內(nèi)西南坡分布范圍較廣,占總面積的21.45%,主要分布有楓香林、含楓香的中草草叢、中國(guó)無(wú)憂(yōu)花-任木林、人面子-野芒果林等植被景觀;西北坡、平地較少,僅占保護(hù)區(qū)總面積的0.2%、0.26%;其余坡向面積分布較均勻,占比在9%~13%。自然植被景觀隨各坡向分布較為均勻,無(wú)明顯的分布規(guī)律;人工林植被則體現(xiàn)出一定的規(guī)律性,如八角林、云南松林等主要分布于陽(yáng)坡與半陽(yáng)坡上(東南、南、西、西南);農(nóng)地、河流等非植被類(lèi)型景觀無(wú)明顯的分布規(guī)律,各坡向均有分布,此類(lèi)景觀受海拔、坡度等地形地貌因子的影響較大。
圖7 不同坡向景觀類(lèi)型分布
(1)在景觀類(lèi)型組成方面,保護(hù)區(qū)以森林植被景觀類(lèi)型為主,其占保護(hù)區(qū)各景觀類(lèi)型總面積的80.29%,且天然植被占絕對(duì)優(yōu)勢(shì)。由此可見(jiàn),天然植被對(duì)維持保護(hù)區(qū)森林植被的穩(wěn)定性有著重要的控制作用,應(yīng)加強(qiáng)對(duì)天然植被保護(hù),減少人為干擾;同時(shí),應(yīng)加大對(duì)天然植被景觀破碎化區(qū)域的營(yíng)造林力度,促進(jìn)斑塊間的物種、物質(zhì)、能量和信息的交換,進(jìn)而保護(hù)物種的多樣性。
(2)保護(hù)區(qū)景觀斑塊規(guī)模主要以小斑塊為主,其占景觀斑塊總數(shù)的85.78%,其余斑塊數(shù)量占比均在3%以下;從斑塊面積來(lái)看,景觀以巨斑塊面積占比最大,達(dá)61.18%,其余斑塊大小等級(jí)分布的面積比例在6%~10%,各景觀類(lèi)型斑塊大小等級(jí)分布面積與斑塊數(shù)量分布存在較大的差異。
(3)保護(hù)區(qū)以低、中山為主,區(qū)內(nèi)700 m 及以下海拔占保護(hù)區(qū)總面積的74.43%,坡度16°~25°的面積占41.39%,除人工林植被景觀外,各景觀類(lèi)型在坡向方面無(wú)明顯分布差異。
(4)本研究?jī)H采用了一期實(shí)地調(diào)查數(shù)據(jù),結(jié)合Landsat8 影像對(duì)保護(hù)區(qū)景觀格局現(xiàn)狀及分布進(jìn)行了分析和評(píng)價(jià)。若要全面掌握保護(hù)區(qū)景觀格局的變化情況,還需對(duì)其開(kāi)展景觀格局動(dòng)態(tài)變化、景觀干擾因素,以及景觀預(yù)測(cè)模型等方面的研究,系統(tǒng)地掌握保護(hù)區(qū)內(nèi)景觀格局的變化和干擾因素,從而提出科學(xué)的干預(yù)手段,確保各景觀類(lèi)型健康、穩(wěn)定。