基于GIS和RS的東北地區東北虎生境適宜性評價

李亞藏，馮仲科，黃季夏，楊　柳（.北京林業大學精準林業北京市重點實驗室，北京00083；2.安陽工學院，河南安陽455000）

?

基于GIS和RS的東北地區東北虎生境適宜性評價

李亞藏1,2，馮仲科1，黃季夏1，楊柳1
（1.北京林業大學精準林業北京市重點實驗室，北京100083；2.安陽工學院，河南安陽455000）

摘要：以中國東北地區（含東北三省和內蒙古東北部）為研究區域，依據穩定性、獨立性、主導性和綜合性原則，從植被、地理環境、人為干擾等3個層面選取植被類型、植被覆蓋度、海拔、坡度、坡向、水源、居民點、主要道路等為評價因子，對中國東北虎Panthera tigris altaica生境進行適宜性評價。在地理信息系統（GIS）和遙感（RS）技術的支持下，獲取各評價因子的定量信息。采用層次分析法（AHP）確定各評價因子的權重，并運用模糊數學模型對各單因子加權疊加，得到東北虎生境適宜性綜合評價分級圖并計算各分區面積。結果表明：①人為干擾和植被因素是影響東北虎分布的主要因素；②將研究區域劃分為東北虎最適宜分布區、次適宜分布區、一般適宜分布區和不適宜分布區，分別占研究區域總面積的21.97%，16.83%，22.02%，39.18%。最適宜區和次適宜區面積總和達到48.80%，說明東北地區適宜東北虎生存的潛在生境面積總量較為可觀。這為東北虎自然保護區的建立和人工飼養東北虎的放歸計劃提供決策依據。圖3表3參25

關鍵詞：動物學；地理信息系統（GIS）；遙感（RS）；東北虎；生境適宜性；數學模型

東北虎Panthera tigris altaica是世界上現存的5個虎Panthera tigris亞種中體形最大者，屬世界瀕危和區域旗艦物種［1］。由于捕獵和棲息地環境的破壞，中國境內的野生東北虎數量已不足20只［2-3］。雖然近些年中國對東北虎的保護做了大量工作，但仍然存在棲息地破碎化、保護區面積不足、人類生產經營活動逐漸向林區發展及生態廊道不暢等問題，嚴重限制了東北虎種群的恢復［4］。另外，隨著人工飼養東北虎數量的逐漸增加［5］，東北虎要放歸何處就成了一個值得思考的問題。目前，對東北虎的研究多集中在東北虎的生態行為學［6-9］、種群分布［2］、數量監測［2, 10-11］及其獵物種群［12-13］等方面，而對東北虎如此廣域生境篩選及其評價的研究較少。中國東北和俄羅斯遠東東北虎分布區的獵物種群密度相對較低，致使東北虎的家域面積較大，據研究，生活在該區1只成年雌虎需要的家域面積約為488 km2，雄虎的領域面積約為1 205 km2［14］；如果要保證種群的可持續繁衍大約需要20只左右成年雌虎，這就意味著至少需要8 000 km2的連續的生境面積，這也說明了生境研究的必要性［15］。本研究從東北虎最基本的生境條件入手，在掌握東北虎及其生境現狀的基礎上，以中國東北地區包括黑龍江、吉林、遼寧三省和內蒙古東北部為研究區域，分析影響東北虎分布的植被、地理環境和人為干擾三大因素，評價和預測東北虎的適宜生境，制定切實可行的規劃措施，以期為東北虎自然保護區的建立和放歸工作提供參考依據。

1　研究區概況與數據來源

1.1研究區概況

該研究區域為中國東北地區，包括黑龍江、吉林、遼寧三省及內蒙古東部地區。地理坐標范圍：38°43′29″～53°41′52″N，111°0′23″～135°53′50″E。全區面積為157.67萬hm2。東北地區自南向北跨中溫帶與寒溫帶，屬溫帶季風氣候，四季分明，夏季溫熱多雨，冬季寒冷干燥。自東南向西北，年降水量自1 000 mm降至300 mm以下，從濕潤區、半濕潤區過渡到半干旱區。該區東、南和北面被水域環繞，西面為陸界，內側是大興安嶺和長白山系的高、中、低山及丘陵，平均海拔為1 000 m；中部為低海拔的東北平原，平均海拔約200 m。東北地區森林覆蓋率高，森林植被以紅松Pinus koraiensis，落葉松Larix olgensis，云杉Picea asperata，山楊Populus davidiana，白樺Betula platyphylla，黃蘗Phellodendron amurense，水曲柳Fraxinus mandschurica，紫椴Tilia amurensis等樹種為主。棲居的獸類有東北虎，豹Panthera pardus，貉Nyctereutes procyonoides，黑熊Ursus thibetanus，紫貂Martes zibellina，水獺Lutra lutra，黃鼬Muastela sibirca等［2］。

1.2數據來源

1.2.1遙感影像遙感影像采用空間分辨率為250 m的MODIS中國區域歸一化植被指數（NDVI）產品和空間分辨率為30 m×30 m的美國陸地衛星Landsat-5遙感影像。NDVI植被指數產品預處理包括：反演、拼接、切割、投影轉換、單位換算等步驟，遙感影像處理在ENVI軟件支持下實現，過程包括輻射定標和大氣校正、幾何精校正及配準、圖像鑲嵌和裁剪、組合波段的選擇與確定、圖像增強處理等步驟。

1.2.2植被圖收集《1∶1 000 000中國植被圖集》［16］，對其進行掃描，以ArcGIS 10.2為平臺，進行矢量化數據和屬性數據錄入，然后利用已有數據進行校正和配準。形成包括針葉林、闊葉林、灌叢、栽培植被、草原、草甸等類型的植被矢量圖。

1.2.3數字高程模型（DEM）DEM數據下載于中國科學院計算機網絡信息中心地理空間數據云，選用GDEM 30 m的產品。

1.2.4東北各省矢量圖省市界、區縣界、道路、河流等矢量圖均來自國家1∶100 000矢量數據集（下載于中國國家基礎地理信息中心）。

2　研究方法

基于遙感（RS）和地理信息系統（GIS）技術，在野外生態調查和收集圖件資料以及相關文本資料［17］的基礎上，依據穩定性、獨立性、主導性和綜合性原則［18］從植被、地理環境、人為干擾3個層面，選取植被類型、植被覆蓋度、海拔、坡度、坡向、水源、居民點、道路8個因子作為生境適宜性評價因子。植被類型結合矢量化的1∶1 000 000中國植被圖，對Landsat-5遙感影像分類得到；植被覆蓋度采用歸一化植被指數NDVI來表現。海拔、坡度和坡向來自30 m空間分辨率的DEM影像數據；水源主要來自河流、湖泊矢量數據，河流級別包括5級及以上；居民點包括地級市和縣級市；道路最小至省道。上述數據在地理信息系統軟件ArcGIS 10.2的支持下，進行數據編輯，并將其轉化為和國家1∶100 000基礎地形數據庫一致的投影坐標系統得到評價所需的各評價因子信息。最后將不同的空間數據離散為90 m×90 m的柵格。再運用適宜性評價模型和GIS的空間分析功能分別計算出每個柵格的生境適宜性綜合評價值，然后制定分級標準，最后得到東北地區東北虎生境適宜性分級圖并計算各級面積。東北地區東北虎生境適宜性評價技術路線如圖1所示。

圖1　東北地區東北虎生境適宜性評價技術路線Figure 1 Technical route of Amur tigers habitat suitability evaluation

3　適宜性評價

3.1單因子生境適宜性評價

依據東北虎生活對棲息地的要求，以及人類活動與自然環境的關系，并綜合考慮前人的研究成果［3,17,19-21］，分析確定東北虎單因子生境適宜度的分級標準（表1），并得到各單因子適宜性評價分級圖。

表1　東北虎單因子生境適宜性分級標準Table 1 Classification standard of Amur tigers habitat suitability based on single factor

3.2評價因子權重的確定

采用層次分析法確定權重［22-23］。首先將指標分為目標層、準則層和方案層3個層次，建立層次結構（圖2），通過專家對判斷矩陣表的填寫構建判斷矩陣，求解矩陣，得到特征向量和最大特征根及最大特征向量，最后對指標進行一致性檢驗，結果為CR=0.015 68＜0.1，具有滿意的一致性。再參考文獻［24］的方法計算出8個指標的組合權重值（表2）。

圖2　東北地區東北虎生境適宜性綜合評價指標體系Figure 2 Evaluation index system of Amur tigers habitat suitability

表2　各評價因子權重值Table 2 Weight of each evaluation factor

3.3生境適宜性綜合評價

3.3.1評價模型的建立應用適宜性評價的模糊數學模型進行綜合評價值的計算。計算公式如下：

式（1）中：i為柵格編號；k為評價因子編號；n為評價因子總數；si為第i個柵格的綜合評價值；wk為第k個因素的權重；ci（k）為第i個柵格的第k個評價因子適宜度評價。

3.3.2單因子加權疊加獲得綜合適宜度通過上面分析已經對各單因子適宜性進行了評價，并確定了各影響因子的權重?，F根據適宜性為各評價因子賦值，以最適宜為標準賦值為1，次適宜、一般適宜和不適宜依次賦值為2，3和4。然后，根據層次分析法求得的權重按照式（1）進行加權求和得到各個柵格的綜合評價值。經過加權疊加得到每個柵格的綜合評分值范圍為1.000 00～3.890 79，若該柵格的所有因子均在最適宜分布區，疊加后的值為最小值1.000 00，說明此處的環境最適宜東北虎生存；當各因子均在不適宜分布區，疊加后為最大值3.890 79；同理，各因子均為次適宜分布區，疊加后的總值是2；均在一般適宜分布區，疊加后的總值為2.968 01。但因為各因子的權重相差懸殊，居民點、道路、植被蓋度和植被類型4個因子的權重總和為0.829 30。因此，為了更為準確地預測東北虎的最適宜分布區，本研究將分區的臨界值做了適當調整。將3個權重最大的因子（居民點、道路和植被蓋度）的最適宜分布區賦值與植被類型的次適宜分布區和其他4個因子的一般適宜分布區賦值的總和作為劃分最適宜區和次適宜區的分界點（1.442 65）；將權重最大的5個因子次適宜區賦值總和作為次適宜區與一般適宜區分界點（1.781 58）；將權重最大的3個因子一般適宜區賦值的總和作為劃分一般適宜區和不適宜區的分界點（2.184 15），重新分類，將研究區分類為最適宜分布區（1～1.442 65），次適宜分布區（1.442 65～1.781 58），一般適宜分布區（1.781 58～2.184 15）和不適宜分布區（2.184 15～3.890 79）。評價結果如圖3所示。利用ArcGIS對4個適宜性等級的面積進行統計，結果如表3所示。

圖3　東北地區東北虎生境適宜性綜合評價分區Figure 3　Comprehensive evaluation distribution of the Amur tigers habitat suitability

表3　東北地區東北虎生境適宜性綜合評價分區面積統計表Table 3　The partition area of comprehensive evaluation of Amur tigers habitat suitability

4　結果與討論

本研究采用層次分析法確定了各評價因子的權重，認為居民點、道路和植被蓋度所占權重最大；其次為植被類型；海拔、坡度、坡向和水源對東北虎的分布影響不大。由此可見，人為干擾和植被因素是影響東北虎分布的主要因素，這與前人的研究結果相一致［1］。以后在進行有關東北虎研究和保護過程中要加強該方面的管理。本研究發現，東北地區東北虎最適宜分布區主要分布于長白山脈、張廣才嶺、小興安嶺、千山山脈和大興安嶺的部分地區，其面積為346 388.4 hm2，占研究區總面積的21.97%。該區植被類型豐富，林分郁閉度高，林分類型多以針、闊混交林或純林為主體，林內有豐富的食物分布，另外該地區城鎮密度小，且靠近邊境地區，人為干擾小。次適宜區面積為265 464.8 hm2，占研究區面積的16.83%，所占比例最小，主要為最適宜區的一些緩沖區域。一般適宜區主要分布在研究區的西南部，該區域居民點稀少，影響東北虎分布的主要因素為植被因素，其植被類型以灌草為主。不適宜區主要集中在中部和南部的東北平原，該區域主要為城鎮分布區，人口密集，植被覆蓋度低，植被類型主要以農田為主，不適宜東北虎生存。張明海等［25］研究認為，已知的野生東北虎出現點主要分布在黑龍江東部及東南部，吉林東部及東北部，均在本研究評價結果范圍內，驗證了評價結果的可靠性。

由于城鎮過度擴張，人工建設大量開發，導致野生東北虎的棲息地遭到嚴重破壞，這使得適宜東北虎生存的潛在生境（面積僅占21.97%）的破碎度較高，即東北虎分布區島嶼化相當嚴重，可能造成了東北虎近年來有向俄羅斯遠東地區遷徙的趨勢。因此，對現有保護區的最適宜分布區建立生態廊道，連接破碎化區域，是當前保護區管理部分急需解決的首要問題。鑒于俄羅斯遠東地區有較多東北虎出沒，建議將中國東北虎保護區與俄羅斯遠東地區的保護區聯合，形成跨國界的保護區網絡。同時，中國東北地區要改變林業經營模式，建立長期的生態環境監測體系。

由于受數據來源的限制，本研究未能將影響東北虎分布的獵物因子納入評價因子當中，可能會對預測結果造成一定程度的偏差，這將是我們今后的工作重點。另外，在利用層次分析法和模糊評價模型進行評價過程中，帶有較強的主觀性，今后考慮能否結合實際情況進行其他評價方法的選擇或模型的改進，為野生東北虎的生存和繁衍以及人工飼養東北虎的放歸計劃提供更為科學的參考。

5　參考文獻

［1］李鐘汶，鄔建國，寇曉軍，等.東北虎分布區土地利用格局與動態［J］.應用生態學報, 200，20（3）：713 - 724.LI Zhongwen, WU Jianguo, KOU Xiaojun, et al.Land use pattern and its dynamic changes in Amur tiger distribution region［J］.Chin J Appl Ecol, 2009, 20（3）: 713 - 724.

［2］周紹春，孫海義，張明海，等.黑龍江省東北虎分布區域及其數量動態［J］.獸類學報, 2008, 28（2）: 165 - 173.ZHOU Shaochun, SUN Haiyi, ZHANG Minghai, et al.Regional distribution and population size fluctuation of wild A-mur tiger（Panthera tigris altaica）in Heilongjiang Province［J］.Acta Theriol Sin, 2008，28（2）: 165 - 173.

［3］曲藝，欒曉峰.基于最小費用距離模型的東北虎核心棲息地確定與空缺分析［J］.生態學雜志，2010，29（9）：1866 - 1874.QU Yi, LUAN Xiaofeng.Core habitat identification and gap analysis for Amur tiger（Panthera tigris altaica）based on least-costdistance model［J］.Chin J Ecol，2010, 29（9）: 1866 - 1874.

［4］李冰，張恩迪，張振華，等.吉林琿春自然保護區東北虎種群的初步監測［J］.獸類學報，2008，28（4）: 333 -341.LI Bing, ZHANG Endi, ZHANG Zhenhua, et al.Prelmiinary monitoring of Amur tiger population in Jilin Hunchun National Nature Reserve［J］.Acta Theriol Sin, 2008，28（4）: 333 - 341.

［5］郭玉榮，劉穎，許魯，等.東北虎林園虎糞的資源化利用研究［J］.野生動物學報，2014，35（1）：15 - 18.GUO Yurong, LIU Ying, XU Lu, et al.Resource value and utilization of tiger feces in Siberian Tiger Park［J］.Chin J Wildlife, 2014, 35（1）: 15 - 18.

［6］劉振生，滕麗微，李楓，等.籠養和半散放東北虎行為時間分配及活動規律差異的比較［J］.生態學雜志，2004，23（1）：20 - 23.LIU Zhensheng, TENG Liwei, LI Feng, et al.Differences of time budget and activity rhythm of captive and semi-free Amur tiger（Panthera tigris altaica）［J］.Chin J Ecol, 2004, 23（1）: 20 - 23.

［7］滕麗微.籠養條件下東北虎行為模式的研究［D］.哈爾濱：東北林業大學，2001.TENG Liwei.A Study on the Behavioral Models of Amur Tiger in Captivity［D］.Harbin: Northeast Forestry University，2001.

［8］KERLEY L L, GOODRICH J M, MIQUELLE D G, et al.Reproductive parameter of wild female Amur（Siberian）tigers（Pantheratighs altaica）［J］.J Mammal, 2003, 84（1）: 288 - 298.

［9］GOODRICH J M.MIQUELLE D G.Transloeation of problem Amur tigers Panthera tigris altaica to alleviate tiger-humanconflicts［J］.Oryx, 2005, 39（4）: 454 - 457.

［10］張常智，張明海，姜廣順.中國野生東北虎數量監測方法有效性評估［J］.生態學報，2012，32（19）：5943 -5952.ZHANG Changzhi, ZHANG Minghai, JIANG Guangshun.Assessment of monitoring method for population abundance of Amur tiger in Northeast China［J］.Acta Ecol Sin, 2012, 32（19）: 5943 - 5952.

［11］樸正吉，睢亞臣，崔志剛，等.長白山自然保護區貓科動物種群數量變化及現狀［J］.動物學雜志，2011，46 （3）：78 - 84.PIAO Zhengji, SUI Yachen, CUI Zhigang, et al.The history and current status of felid population in Changbai Mountain Nature Reserve［J］.Chin J Zool, 2011, 46（3）: 78 - 84.

［12］張常智，張明海.黑龍江省東完達山地區東北虎獵物種群現狀及動態趨勢［J］.生態學報，2011，31（21）：6481 - 6487.ZHANG Changzhi, ZHANG Minghai.Population status and dynamic trends of Amur tiger’s prey in Easter Wandashan Mountain, Heilongjiang Province［J］.Acta Ecol Sin, 2011, 31（21）: 6481 - 6487.

［13］周紹春，張明海，孫海義，等.黑龍江省完達山東部林區東北虎獵物生物量［J］.生態學報, 2011, 31（1）：145 - 153.ZHOU Shaochun, ZHANG Minghai, SUN Haiyi, et al.Prey biomass of the Amur tiger（Panthera tigris altaica）in the eastern Wanda Mountains of Heilongjiang Province［J］.Acta Ecol Sin, 2011, 31（1）: 145 - 153.

［14］李穎.揮春地區基于東北虎生境選擇的生態廊道研究［D］.吉林：延邊大學，2011.LI Ying.The Research on Ecology Corridor of Amur Tiger Based on Habitat Choice Research in Hunchun Area［D］.Jilin: Yanbian University, 2011.

［15］羅述金.中國虎的概況［J］.生物學通報，2010，45（1）：1 - 5.LUO Shujin.Tigers in China［J］.Bull Biol, 2010, 45（1）：1 - 5.

［16］中國科學院中國植被圖編輯委員會.1∶1 000 000中國植被圖集［M］.北京：科學出版社，2001.

［17］周紹春.東北虎及其獵物的種群大小、生境選擇與評價研究［D］.哈爾濱：東北林業大學，2011.ZHOU Shaochun.Relationship Between Amur Tiger and Prey: Population Size, Habitat Seleetion and Evaluation in Eastem Wandashan Mountains of Heilongjiang Province［D］.Harbin: Northeast Forestry University, 2011.

［18］周建飛，曾光明，黃國和，等.基于不確定性的城市擴展用地生態適宜性評價［J］.生態學報，2007，27（2）：774 - 783.ZHOU Jianfei, ZENG Guangming, HUANG Guohe, et al.The ecological suitability evaluation on urban exoansion land based on uncertainties［J］.Acta Ecol Sin, 2007, 27（2）: 774 - 783.

［19］李振新.中國長白山地區東北虎潛在棲息地研究［M］.北京：中國林業出版社，2010.

［20］馬天.黑龍江省野生東北虎調查管理信息系統的研建［D］.北京：北京林業大學，2013.MA Tian.Research and Construction of Investigation Information Management System of Panthera tigris altaica in Heilongjiang Province［D］.Beijing: Beijing Forestry University, 2013.

［21］杜婧.黑龍江省野生東北虎生境適宜度評價［D］.北京：北京林業大學，2014.DU Jing.The Suitability Evaluation of the Wild Amur Tiger Habitat in Heilongjiang Province［D］.Beijing: Beijing Forestry University, 2014.

［22］黃輝玲，吳次芳，張守忠.黑龍江科技學院建筑黑龍江省土地整治規劃效益分析與評價［J］.農業工程學報，2012，28（6）：240 - 246.HUANG Huiling, WU Cifang, ZHANG Shouzhong.Benefits analysis and evaluation on land consolidation lanning in Heilongjiang province［J］.Trans Chin Soc Agric Eng, 2012, 28（6）: 240 - 246.

［23］向欣，羅煜，程紅勝，等.基于層次分析法和模糊綜合評價的沼氣工程技術篩選［J］.農業工程學報，2014，30（18）：205 - 212.XIANG Xin, LUO Yu, CHENG Hongsheng, et al.Biogas engineering technology screening based on analytic hierarchy process and fuzzy comprehensive evaluation［J］.Trans Chin Soc Agric Eng, 2014, 30（18）: 205 - 212.

［24］鮑艷，胡振琪，王建峰，等.層次分析法在土地開發中的適宜性評價［J］.西安科技大學學報，2005，25（2）：179 - 183.BAO Yan, HU Zhenqi, WANG Jianfeng, et al.Analytical hierarchy process based on GIS assessment of suitability in land exploitation and consolidation［J］.J Xi’an Univ Sci Technol, 2005, 25（2）: 179 - 183.

［25］張明海，馬建章.中國野生東北虎現狀及其保護愿景展望［J］.動物學雜志，2010，45（1）：165 - 168.ZHANG Minghai, MA Jianzhang.Current status and protection vision of wild northeast tiger in China［J］.Chin J Zool, 2010, 45（1）: 165 - 168.

Habitat suitability evaluation for Amur tigers in northeast China based on GIS and RS

LI Yacang1,2, FENG Zhongke1, HUANG Jixia1, YANG Liu1
（1.Institute of GIS, RS, and GPS, Beijing Forestry University, Beijing 100083, China; 2.Anyang Institute of Technology, Anyang 455000, Henan, China）

Abstract:This article researched the habitat suitability evaluation for Amur tigers in northeast China, which found the potential distribute area and provided decision-making for the establishment of Amur tiger nature reserve and the plan of getting back to nature of artificial breeding Amur tigers.The study area in northeast China included the three provinces of northeast China and the northeast region of Inner Mongolia.According to the principles of stability, independence, dominance, and comprehensiveness, the habitat suitability of Amur tigers in northeast China was assessed using the three aspects of vegetation, geographical environment, and human disturbance with vegetation type, vegetation coverage, elevation, gradient, slope direction, water source, residential areas, and main roads selected as evaluation variables.Data for each variable was collected quantitatively using geographic information systems（GIS）and remote sensing（RS）.Then the Amur tiger habitat suitability evaluation classification figure was obtained using the Analytic Hierarchy Process（AHP）method to determine the weight of each evaluation factor, and using the fuzzy mathematical model for the single factor weighted superposition.Afterward, the partitioning of areas was conducted.Results showed that human disturbance and vegetation were the main factors affecting the Amur tiger distribution.The study area was classified as most appropriate（21.97%）, second suitable（16.83%）, generally appropriate（22.02%）, and inappropriatebook=266,ebook=87distribution areas（39.18%）.Thus, the living space that could be appropriated for Amur tigers compared to the total potential habitat area in the Northeast China was considerable approximately 48.8%.［Ch, 3 fig.3 tab.25 ref.］

Key words:zoology; geographic information system（GIS）; remote sensing（RS）; Amur tiger; habitat suitability; mathematical model

作者簡介：李亞藏，副教授，從事地理信息系統在資源環境上的應用研究。E-mail：liyacang@sohu.com。通信作者:馮仲科，教授，博士生導師，從事精準林業與3S技術研究。E-mail：fengzhongke@126.com

基金項目：“十二五”國家科技支撐計劃項目（2012BAH34B01）

收稿日期：2015-03-19；修回日期：2015-09-01

doi:10.11833/j.issn.2095-0756.2016.02.011

中圖分類號：S759.9；Q959.9

文獻標志碼：A

文章編號：2095-0756（2016）02-0265-07