周 勇,苗潤(rùn)澤,任 毅,楊秀冬,黎 勇,劉 庚
(1.吉林省林業(yè)科學(xué)研究院,吉林 長(zhǎng)春 130033;2.野生生物保護(hù)學(xué)會(huì)(美國(guó))北京代表處,北京 100012;3.吉林龍灣國(guó)家級(jí)自然保護(hù)區(qū)管理局,吉林 通化 135102;4.柳河縣八里哨林場(chǎng),吉林 通化 135321)
生境是野生動(dòng)物進(jìn)行繁衍和生命活動(dòng)的場(chǎng)所,是物種生存和發(fā)展所需生物和非生物因子的集合[1].隨著環(huán)境變化和人類活動(dòng)干擾的加強(qiáng),日益嚴(yán)重的物種生境破碎化問題[2-3]導(dǎo)致物種棲息地適宜生境面積減少、質(zhì)量下降,影響物種的擴(kuò)散和遷移,干擾物種種群間的基因交流,減少物種多樣性[4-5].東北虎(Pantheratigrisaltaica)和東北豹(Pantheraparduaorientalis)是全球生物多樣性保護(hù)的旗艦物種,在維持健康生態(tài)系統(tǒng)功能中占據(jù)不可替代的重要地位.近年來,人類活動(dòng)對(duì)東北虎豹天然生態(tài)廊道產(chǎn)生嚴(yán)重干擾,天然生態(tài)廊道消失殆盡,保護(hù)東北虎豹首先就要加強(qiáng)棲息地生境保護(hù)與生境連通[6-7],連接島狀分布的適宜、潛在生境,識(shí)別東北虎豹潛在的棲息、繁衍、遷徙和擴(kuò)散的重要生境廊道[8],這對(duì)保障種群向內(nèi)陸擴(kuò)散、增加基因交流、擴(kuò)大種群數(shù)量、維持和提高物種多樣性具有重要意義[9-10].
目前,有關(guān)生境廊道的研究、定義有很多,于成龍等、戴璐等[11-12]基于最小累積阻力模型,將提取的最小成本路徑視為生態(tài)廊道;李瑾等[13]利用水文分析原理,將累積耗費(fèi)距離最短的山谷線稱為生態(tài)廊道;Wang Chenxu等[14]運(yùn)用ArcGIS中的Linkage Mapper插件提取生態(tài)廊道,并給出了廊道范圍內(nèi)起關(guān)鍵作用的生態(tài)節(jié)點(diǎn);張遠(yuǎn)景等[15]通過“斑塊-廊道-基質(zhì)”理論,分析了哈爾濱市生態(tài)網(wǎng)絡(luò)連接度強(qiáng)弱的空間分布.
在吉林省東北虎豹主要分布區(qū)琿春市,國(guó)道(G331)穿境而過,是東北虎豹棲息地阻隔最重要的干擾因子,也對(duì)棲息地、居民點(diǎn)和有蹄類分布產(chǎn)生持續(xù)影響.但是,有關(guān)這種干擾因子是如何影響東北虎豹生境整體連通性的研究尚未見報(bào)道,因此,本文以琿春市東北虎豹分布區(qū)域?yàn)閷?shí)例,運(yùn)用廊道設(shè)計(jì)模型(Linkage Mapper)提取生態(tài)廊道和生態(tài)夾點(diǎn),分析識(shí)別了琿春市國(guó)道周邊東北虎豹的廊道區(qū)域,以為東北虎豹棲息地生境恢復(fù)和種群西遷提供科學(xué)依據(jù).
吉林省琿春市(130°03′21″~130°18′33″E,42°25′20″~43°30′18″N)位于吉林省東南部圖們江下游地區(qū)的延邊朝鮮族自治州境內(nèi)(見圖1),地處中、朝、俄三國(guó)交界地帶,轄區(qū)面積5 145 km2,主要包括中山、低山和丘陵河谷地貌3種地貌類型.林地面積為4 370 km2,森林覆蓋率達(dá)85%,以針葉林、闊葉林和闊葉混交林為主;氣候溫和濕潤(rùn),屬于溫帶海洋性氣候,冬暖夏涼,年均氣溫5.65 ℃,平均降水量617.9 mm,無霜期140~160 d,8月份平均氣溫21.2 ℃.境內(nèi)的東北虎豹國(guó)家公園以琿春林業(yè)局青龍臺(tái)林場(chǎng)東部為起點(diǎn),向吉林省汪清縣和黑龍江東寧市輻射,是東北虎豹種群向內(nèi)陸擴(kuò)散的主要通道;轄區(qū)內(nèi)國(guó)道(G331)穿境而過,全長(zhǎng)188 km,北通黑龍江東寧市、南抵中朝俄邊界的防川村,沿途分布市區(qū)、鄉(xiāng)鎮(zhèn)和農(nóng)村等人類活動(dòng)區(qū)..
土地利用數(shù)據(jù)來源于國(guó)家基礎(chǔ)地理信息中心(http:∥www.ngcc.cn/ngcc/)2020年全球30 m地表覆蓋數(shù)據(jù),根據(jù)琿春市行政邊界矢量裁剪獲取的森林、水體、耕地、建設(shè)用地、草本沼澤等6類土地利用數(shù)據(jù);DEM數(shù)據(jù)來源于地理空間數(shù)據(jù)云(http:∥www.gscloud.cn/)ASTERGDEM 30 m分辨率數(shù)字高程數(shù)據(jù);河流、道路、居民點(diǎn)、琿春市行政邊界矢量數(shù)據(jù)來自國(guó)家地理信息中心的1∶250 000基礎(chǔ)地理信息數(shù)據(jù);郁閉度數(shù)據(jù)來自琿春市森林資源檔案;有蹄類數(shù)據(jù)為2021年有蹄類野生動(dòng)物野外調(diào)查發(fā)現(xiàn)的點(diǎn)位數(shù)據(jù),調(diào)查對(duì)象為狍(Capreoluspygargus)和野豬(Susscrofa).
研究方法包括4部分:生態(tài)源地篩選、阻力面構(gòu)建、生態(tài)廊道綜合識(shí)別和夾點(diǎn)、障礙點(diǎn)識(shí)別.
1.3.1 生態(tài)源地的篩選
根據(jù)琿春市森林資源生境類型和面積將生態(tài)源地斑塊進(jìn)行處理整合.一是斑塊生境類型的選擇:根據(jù)東北虎豹生存活動(dòng)適宜生境選擇針葉林、闊葉林和針闊混交林3種生境類型;二是國(guó)道單側(cè)生態(tài)源地選取:根據(jù)研究目的,并考慮居民點(diǎn)、道路和河流對(duì)源地的影響進(jìn)行選取;三是細(xì)碎斑塊和獨(dú)立斑塊選取:刪除鄰近市區(qū)且相隔距離較遠(yuǎn)的獨(dú)立斑塊,選取相隔距離近且適宜野生動(dòng)物遷徙的細(xì)碎斑塊.
1.3.2 阻力面構(gòu)建
(1) 阻力因子選取.根據(jù)琿春市自然資源狀況、人為主導(dǎo)的景觀分布情況,通過考察東北虎豹適宜生境,與相關(guān)研究、保護(hù)人員開展訪談、調(diào)查,選定土地利用類型、森林郁閉度、坡度、海拔、道路、河流、居民點(diǎn)和有蹄類距離等8類阻力因子,并構(gòu)建各因子生態(tài)阻力基面.
(2) 阻力值確定及綜合阻力面構(gòu)建.利用ArcGis軟件對(duì)阻力因子進(jìn)行分級(jí)和相關(guān)性分析,結(jié)合研究區(qū)域?qū)嶋H,運(yùn)用“鄰域分析”工具對(duì)道路、河流和有蹄類距離3類阻力因子做緩沖區(qū)確定分析,最后設(shè)定道路中林間土路、一般硬化路面、主要公路和國(guó)道(G331)緩沖區(qū)距離分別為10,15,25和25 m;河流中山澗溪流、村邊河流、支流和主要河流緩沖區(qū)距離分別為3,10,30和50 m;有蹄類距離中近、較近和較遠(yuǎn)緩沖區(qū)距離分別為500,1 000和2 000 m.
參考已有的研究[16],確定8類阻力因子的分級(jí)標(biāo)準(zhǔn)及各級(jí)別阻力值.為更好滿足Linkage Mapper插件數(shù)據(jù)要求,將棲息地因子中各阻力值限定在1~100區(qū)間,其中1代表最適宜、100代表不適宜;將食源因子、有蹄類緩沖區(qū)距離設(shè)定為-100~-1區(qū)間,其中-100代表最適宜、-1代表不適宜.阻力因子與阻力值設(shè)定結(jié)果見表1.為避免各阻力因子之間高度的自相關(guān)對(duì)模型分析產(chǎn)生干擾,運(yùn)用ArcGis軟件中的“波段集統(tǒng)計(jì)”分析工具,對(duì)8類阻力因子進(jìn)行相關(guān)性分析.
表1 阻力因子與阻力值
依據(jù)專家打分法,通過相關(guān)研究專家和保護(hù)管理人員的評(píng)分,確定各阻力因子權(quán)重,結(jié)果見表2.利用“地圖代數(shù)”分析工具對(duì)8類阻力因子阻力面加權(quán)疊加構(gòu)建綜合阻力面.各阻力基面柵格數(shù)據(jù)輸出過程中,各柵格數(shù)據(jù)的波段數(shù)位為1、像素深度為8位;除道路中林間土路、一般硬化路面和主要公路柵格數(shù)據(jù)像元大小為10 m×10 m,其余阻力因子?xùn)鸥駭?shù)據(jù)像元大小均為30 m×30 m.
表2 阻力因子權(quán)重
1.3.3 生態(tài)廊道綜合識(shí)別
利用Linkage Pathways Tool工具,通過導(dǎo)入的東北虎豹生態(tài)源地和綜合阻力面數(shù)據(jù),計(jì)算東北虎豹在源地間遷徙經(jīng)過各阻力基層面所克服的累積阻力值,遷移擴(kuò)散的最小成本路徑(Least-cost path method,LCP)即為生態(tài)廊道路徑.本研究中截?cái)嗌鷳B(tài)走廊的成本加權(quán)距離閾值設(shè)置為200 000[17].
1.3.4 夾點(diǎn)、障礙點(diǎn)識(shí)別
Linkage Mapper插件通過電路理論和障礙分析來檢測(cè)最小成本路徑中的夾點(diǎn)和障礙點(diǎn),以確定重要的核心區(qū)域和走廊;Circuitscape程序是基于電路理論對(duì)不同景觀進(jìn)行連接建模[18],從而將電路理論與運(yùn)動(dòng)生態(tài)學(xué)相聯(lián)系.本文運(yùn)用Pinchpoint Mapper工具識(shí)別生態(tài)夾點(diǎn),在綜合阻力基柵格數(shù)據(jù)輸出過程中,柵格數(shù)據(jù)的波段數(shù)位為1、像素深度為32位、像元大小為10 m×10 m;選擇“all-to-one”迭代運(yùn)算模式,調(diào)用Circuitscape程序?qū)?3個(gè)東北虎豹生態(tài)源地斑塊之間的電流進(jìn)行組合,識(shí)別廊道內(nèi)生態(tài)夾點(diǎn)位置.鑒于廊道寬度不會(huì)影響夾點(diǎn)位置與區(qū)域的連通性,本文設(shè)置廊道寬度加權(quán)成本距離為300,以電流值評(píng)估夾點(diǎn)對(duì)于維持整個(gè)景觀連接性的重要性[19],提取的高電流強(qiáng)度區(qū)域即為夾點(diǎn)區(qū)域.
本文運(yùn)用Barrier Mapper工具識(shí)別生態(tài)障礙點(diǎn).在綜合阻力基柵格數(shù)據(jù)輸出時(shí),柵格數(shù)據(jù)的波段數(shù)位為1、像素深度為32位、像元大小為10 m×10 m,此時(shí)最小檢測(cè)半徑為20 m、最大檢測(cè)半徑為80 m、半徑步長(zhǎng)值為20 m.不勾選“改進(jìn)得分相對(duì)于最小成本路徑百分比選項(xiàng)”,得到最大改進(jìn)分?jǐn)?shù),以改進(jìn)分?jǐn)?shù)評(píng)估修復(fù)障礙點(diǎn)對(duì)整個(gè)景觀連通性的重要性,提取的高改進(jìn)分?jǐn)?shù)區(qū)域即為障礙點(diǎn)區(qū)域.
1.3.5 廊道寬度識(shí)別
國(guó)道(G331)上的廊道寬度通常定義為有一定障礙但具有較好景觀連通性的道道的寬度[20].運(yùn)用Barrier Mapper工具,在檢測(cè)半徑20,40,60和80 m條件下,計(jì)算耗費(fèi)累計(jì)阻力值最小的廊道寬度,利用SPSS 19.0軟件對(duì)4類檢測(cè)半徑下的12個(gè)廊道寬度進(jìn)行成對(duì)樣本T檢驗(yàn).
依據(jù)生境類型及斑塊分布,將研究區(qū)域內(nèi)的針葉林、闊葉林和針闊混交林劃分為23個(gè)生態(tài)源地,結(jié)果見圖2,生態(tài)源地總面積4 462.32 km2.
圖2 琿春市主要生態(tài)源地分布
在選定的土地利用類型、森林郁閉度、坡度類型、海拔類型、道路類型、河流類型、居民點(diǎn)類型和有蹄類距離8類阻力因子中,棲息地因子為東北虎豹遷移擴(kuò)散的重要阻礙.土地利用類型中針闊混交林、闊葉林相對(duì)于針葉林更適宜東北虎豹的生存,大水域阻隔了有蹄類動(dòng)物的遷徙,耕地、建設(shè)用地對(duì)東北虎豹遷徙造成強(qiáng)烈的人為干擾,草本沼澤對(duì)東北虎豹的遷徙耗能、食源豐富度均有影響.森林郁閉度影響東北虎豹自身的隱藏、對(duì)食物的追蹤.適宜的坡度可減少東北虎豹能量的消耗.東北虎豹在中山區(qū)域的生存阻礙小于低山和丘陵區(qū)域.硬化路面對(duì)生境斑塊的阻隔作用要強(qiáng)于林間土路,屬極大的人為干擾.村邊河流阻礙遷徙主要因?yàn)檩^大的人為干擾,支流、山澗溪流則是遷徙過程中的重要水源.居民點(diǎn)是居民的生產(chǎn)生活區(qū),存在較大的人為干擾,遠(yuǎn)離居民點(diǎn)則更有利于東北虎豹的擴(kuò)散[21].食源因子可促進(jìn)遷徙擴(kuò)散,與有蹄類野生動(dòng)物位置越近則更利于東北虎豹的棲息、遷徙.
根據(jù)阻力因子的阻力值及權(quán)重,對(duì)全部因子進(jìn)行相關(guān)性分析,結(jié)果見表3.由表3可見,各阻力因子相關(guān)系數(shù)均小于0.553,所以各因子之間不存在相關(guān)性[22].
表3 阻力因子相關(guān)性分析
綜合阻力面計(jì)算結(jié)果中,柵格顏色越深表示阻力越大,即東北虎豹種群擴(kuò)散消耗的能量越多,越不利于種群擴(kuò)散.由研究區(qū)綜合阻力面計(jì)算結(jié)果(見圖3)可見,東北虎種群擴(kuò)散阻力主要分布在國(guó)道及兩側(cè)的連接部分;此外城區(qū)及周邊區(qū)域顏色最深,表示該區(qū)域?qū)|北虎豹擴(kuò)散影響最大.
圖3 綜合阻力面分布
在國(guó)道(G331)上,東北虎豹的潛在通道主要由廊道寬度和最短路徑長(zhǎng)度兩部分標(biāo)識(shí).本文通過Linkage Pathways Tool工具識(shí)別廊道12條(見表4);基于Pinchpoint Mapper和Barrier Mapper工具識(shí)別生態(tài)節(jié)點(diǎn)13個(gè),包括6個(gè)夾點(diǎn)區(qū)域,此區(qū)域?qū)Ψ乐箺⒌赝嘶嘘P(guān)鍵作用,應(yīng)重點(diǎn)保護(hù)或優(yōu)先修復(fù);識(shí)別7個(gè)障礙點(diǎn)區(qū)域,結(jié)合土地利用類型應(yīng)對(duì)此區(qū)域進(jìn)行一定的修復(fù)以增強(qiáng)區(qū)域內(nèi)景觀的連通性.對(duì)4類檢測(cè)半徑下12個(gè)廊道寬度進(jìn)行的成對(duì)樣本T檢驗(yàn)結(jié)果表明,60和80 m步長(zhǎng)情況下長(zhǎng)度差異顯著(t=-2.448,P=0.032);排除80 m步長(zhǎng)情況,分別比較20,40 和60 m步長(zhǎng)類別下,距離與平均距離(16 317 m)間的差異為1 698,153 和1 545 m.以基于40 m步長(zhǎng)的距離為廊道寬度,識(shí)別廊道寬度12個(gè).生境廊道、夾點(diǎn)和障礙點(diǎn)示意圖見圖4—6.
表4 廊道基本信息
圖4 生境廊道影像示意圖
圖5 夾點(diǎn)示意圖
圖6 障礙點(diǎn)示意圖
(1) 本研究主要是基于貫穿琿春市域的國(guó)道(G331)進(jìn)行生態(tài)廊道的提取,在生態(tài)源地的識(shí)別中,考慮到研究區(qū)域的高森林覆蓋率及森林資源的連續(xù)性,以國(guó)道為界將東北虎豹棲息地分為東西兩部分,分別合并生境類型相同或相似的斑塊,從而在南北方向上將本來連通的區(qū)域劃分為多處源地,這在一定程度上對(duì)生態(tài)廊道的提取存在影響,但這種影響也僅是表現(xiàn)在提取廊道的重要性方面.
(2) 本研究在考慮東北虎豹的食源時(shí),一是以狍和野豬代表東北虎豹的有蹄類食物,對(duì)于主要獵物的馬鹿(Cervuscanadensis)和梅花鹿(Cervusnippon)沒有進(jìn)行全面調(diào)查[23];二是通過距有蹄類距離替代不同區(qū)域內(nèi)有蹄類野生動(dòng)物的種群密度以表征食物對(duì)東北虎豹的吸引,這兩點(diǎn)有待加強(qiáng).
(3) 本研究基于國(guó)道(G331)的生態(tài)廊道、夾點(diǎn)、障礙點(diǎn)的研究思路和方法,可為人造工程與野生動(dòng)物遷徙保護(hù)之間的矛盾找出科學(xué)的解決辦法;同時(shí),也可為其他物種遷徙及保護(hù)提供科學(xué)參考.
(4) Hodgson等[24]的研究表明,很少有哺乳動(dòng)物在穿越棲息地斑塊中使用預(yù)測(cè)的廊道,因此,為使預(yù)測(cè)的廊道結(jié)果更為準(zhǔn)確,可通過在預(yù)測(cè)廊道區(qū)域架設(shè)紅外相機(jī)進(jìn)行驗(yàn)證.
本文在結(jié)合了琿春市自然及社會(huì)資源情況下,基于國(guó)道(G331)開展了東北虎豹生境廊道的研究,結(jié)果識(shí)別東北虎豹潛在生境廊道12條、夾點(diǎn)6個(gè)和障礙點(diǎn)7個(gè),這對(duì)促進(jìn)東北虎豹在國(guó)道兩側(cè)的遷徙連通,并進(jìn)一步優(yōu)化提升東北虎豹的生存空間具有一定的現(xiàn)實(shí)意義.在主要阻力因子中,土地利用類型、森林郁閉度、有蹄類距離和坡度對(duì)綜合阻力面有較大影響;對(duì)于提取的生態(tài)節(jié)點(diǎn)、夾點(diǎn)區(qū)域?qū)Ψ乐箺⒌赝嘶嘘P(guān)鍵作用,應(yīng)重點(diǎn)保護(hù)或優(yōu)先修復(fù),對(duì)障礙點(diǎn)區(qū)域要結(jié)合土地利用類型進(jìn)行修復(fù)以增強(qiáng)區(qū)域內(nèi)景觀的連通性.本文選定40 m步長(zhǎng)距離為廊道寬度,只是在最低程度上保障東北虎豹的遷徙,隨著選定步長(zhǎng)的變化,夾點(diǎn)區(qū)域要保護(hù)或修復(fù)的面積可能會(huì)隨之改變.