氣候變化下大興安嶺地區(qū)紫貂優(yōu)先保護(hù)區(qū)分析

蔣麗華，高俊琴①，萬(wàn)基中，2，劉吉平

(1.北京林業(yè)大學(xué)自然保護(hù)區(qū)學(xué)院，北京 100083；2.青海大學(xué)省部共建立三江源生態(tài)與高原農(nóng)牧業(yè)國(guó)家重點(diǎn)實(shí)驗(yàn)室，青海 西寧 810016；3.吉林師范大學(xué)旅游與地理科學(xué)學(xué)院，吉林 四平 136000)

近年來(lái)，氣候變化已經(jīng)引起物種空間分布格局發(fā)生巨大變化，物種分布范圍趨于向兩極或高海拔地區(qū)移動(dòng)[1]，生物多樣性也在下降[2]，給物種保護(hù)帶來(lái)更大挑戰(zhàn)。以往研究表明，氣候變化給某些物種潛在分布和適宜棲息地帶來(lái)很大影響[3-5]。ERASMUS等[3]認(rèn)為到2050年，南非179個(gè)物種中有78%的物種分布區(qū)將縮小，2%的物種將遭遇局地滅絕。在我國(guó)，東北地區(qū)白枕鶴(Grusvipio)當(dāng)前的潛在適宜生境將減少[5]，大熊貓(Ailuropodamelanoleuca)適宜生境到2070年代將消失52.9%～71.3%，且生境破碎化加劇[4]。

紫貂(Marteszibellina)屬于哺乳綱(Mammalia)食肉目(Carnivora)鼬科(Mustelidae)貂屬(Martes)，分布于俄羅斯、蒙古、中國(guó)、朝鮮和日本等地[6]，是我國(guó)一級(jí)保護(hù)動(dòng)物，被《中國(guó)生物多樣性紅色名錄——脊椎動(dòng)物卷》(2015)列為易危(VU)物種。大興安嶺地區(qū)是紫貂重要分布區(qū)，毗鄰俄羅斯，是野生動(dòng)物往返于西伯利亞和中國(guó)東北部地區(qū)的重要遷徙通道，對(duì)物種保護(hù)具有重要意義。大興安嶺地處寒溫帶氣候區(qū)，氣候變化敏感，其豐富的自然資源為紫貂提供了棲息環(huán)境和食物來(lái)源。但1950年代以來(lái)，我國(guó)東北地區(qū)紫貂的種群數(shù)量持續(xù)下降，其原因既有林木采伐和非法盜獵等人類干擾因素，也有氣候變化等自然因素[7]。近年來(lái)，我國(guó)大興安嶺地區(qū)大力開展天然林保護(hù)工程、林木禁伐限伐、生態(tài)移民工程等保護(hù)措施，使林區(qū)自然資源得到有效保護(hù)和恢復(fù)，人類干擾明顯減弱。在此基礎(chǔ)上，氣候變化所帶來(lái)的影響將成為紫貂保護(hù)工作的首要挑戰(zhàn)[7]。模擬分析不同氣候條件下紫貂的潛在分布，對(duì)其保護(hù)具有重要意義。最大熵(MaxEnt)模型被廣泛運(yùn)用于模擬物種潛在適宜生境分布格局[8]，在有限的痕跡點(diǎn)下其預(yù)測(cè)結(jié)果比其他生態(tài)位模型更為可靠[9]。同時(shí)，越來(lái)越多的研究者將物種分布模型與Zonation軟件相結(jié)合，開展物種保護(hù)規(guī)劃研究[10]。因此，以大興安嶺地區(qū)紫貂為研究對(duì)象，分析氣候變化對(duì)紫貂潛在適宜生境的影響，確定其優(yōu)先保護(hù)區(qū)域，并分析其現(xiàn)有保護(hù)空缺，可為大興安嶺地區(qū)紫貂的有效保護(hù)提供建議。

1 材料與方法

1.1 物種分布數(shù)據(jù)

紫貂分布數(shù)據(jù)來(lái)源：(1)中國(guó)國(guó)家標(biāo)本資源平臺(tái)(http:∥www.nsii.org.cn/)；(2)GBIF數(shù)據(jù)庫(kù)(https:∥www.gbif.org/)；(3)保護(hù)地資源本底調(diào)查或者科學(xué)考察報(bào)告；(4)公開發(fā)表的文獻(xiàn)資料。為確保數(shù)據(jù)準(zhǔn)確性，紫貂分布信息按以下原則進(jìn)行篩選：(1)文獻(xiàn)優(yōu)先級(jí)由高到低依次為中國(guó)國(guó)家標(biāo)本資源平臺(tái)、保護(hù)地資源本底調(diào)查、GBIF數(shù)據(jù)庫(kù)和公開發(fā)表的文獻(xiàn)資料；(2)保留2000年以來(lái)縣級(jí)及縣級(jí)以下物種分布數(shù)據(jù)，去除重復(fù)分布點(diǎn)；(3)對(duì)來(lái)源于文獻(xiàn)資料的物種分布信息，按照所發(fā)表期刊影響因子大小和發(fā)表時(shí)間順序進(jìn)行篩選。最終篩選獲取紫貂分布信息及來(lái)源見表1。

表1 紫貂痕跡點(diǎn)信息Table 1 Information source of presences of Martes zibellina in Daxing′anling

基于紫貂分布數(shù)據(jù)和基礎(chǔ)空間數(shù)據(jù)建立柵格分辨率為10 km×10 km的物種分布數(shù)據(jù)庫(kù)。10 km×10 km的分辨率被廣泛運(yùn)用于區(qū)域尺度的物種分布建模[19]。由于獲取的紫貂分布均分布在縣級(jí)尺度以下，為防止高估紫貂分布區(qū)域，將紫貂分布與植被、高程、坡度和坡向疊加，保留包含紫貂棲息地類型的柵格為紫貂分布區(qū)，且每個(gè)柵格最多包含1個(gè)紫貂痕跡點(diǎn)。具體操作如下：(1)采用ArcGIS 10.0軟件建立大興安嶺地區(qū)10 km×10 km格網(wǎng)數(shù)據(jù)，并轉(zhuǎn)換為KML格式文件；(2)基于大興安嶺地區(qū)1∶100萬(wàn)植被分布圖提取適宜紫貂生存的植被類型，轉(zhuǎn)換為KML格式文件；(3)將1∶100萬(wàn)行政區(qū)劃圖(國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)共享服務(wù)平臺(tái)，http:∥www.geodata.cn/)轉(zhuǎn)換為KML格式文件；(4)將以上3個(gè)KML文件導(dǎo)入Google Earth軟件，根據(jù)物種分布區(qū)和生境信息確定有紫貂分布的柵格；(5)將Google Earth軟件建立的物種分布點(diǎn)導(dǎo)入ArcGIS 10.0軟件，獲得紫貂痕跡點(diǎn)(96個(gè))分布圖(圖1)。

1.2 環(huán)境變量選擇

1.2.1氣候變量

氣候變量選擇年均溫(Bio1)、平均氣溫日較差(Bio2)、最熱月最高溫(Bio5)和最冷月最低溫(Bio6)4個(gè)氣候因子，數(shù)據(jù)來(lái)源于世界氣候數(shù)據(jù)庫(kù)1.4版本(http:∥www.worldclim.org/)，分辨率為5.0弧分(約為10 km×10 km)。這4個(gè)因子在氣候變化條件下與紫貂分布顯著相關(guān)，相關(guān)系數(shù)∣r∣≥0.847[7]。RCP(representative concentration pathway)為未來(lái)人為溫室氣體排放濃度，RCP 4.5為中等溫室氣體排放情景，表示全球年度溫室氣體排放量在2040年左右達(dá)到排放峰值，之后下降；RCP 8.5為最高溫室氣體排放情景，表示排放量在整個(gè)21世紀(jì)持續(xù)上升。當(dāng)前氣候條件下的氣候變量采用1960—1990年全球多個(gè)站點(diǎn)觀測(cè)數(shù)據(jù)的平均值；未來(lái)氣候變量由BCC-CSM1-1、MIROC-ESM-CHEM和MIROC-ESM 3種常用全球氣候模型(global climate models，GCMs)對(duì)2041—2060和2061—2080年兩個(gè)階段分別在RCP 4.5和RCP 8.5兩種情景下預(yù)測(cè)得到。選取上述3種氣候模型進(jìn)行模擬預(yù)測(cè)，以減少單一模型對(duì)模擬結(jié)果的影響[20]。

Hjd為痕跡點(diǎn)，Bhd為保護(hù)地，Dxal為大興安嶺地區(qū)，Xj為縣界，Sj為省區(qū)界，Gj為國(guó)界。審圖號(hào)：GS(2019)3621號(hào)圖1 大興安嶺地區(qū)保護(hù)地及紫貂痕跡點(diǎn)分布Fig.1 Distribution of protected areas and presences of Martes zibellina in Daxing′anling

1.2.2地形變量

地形變量選擇海拔、坡度和坡向3個(gè)地形因子，數(shù)據(jù)來(lái)源于美國(guó)地質(zhì)調(diào)查局(http:∥reverb.echo.nasa.gov/reverb/redirect/wist)數(shù)字高程模型，分辨率為30 m×30 m，采用ArcGIS 10.0軟件得到坡度和坡向數(shù)據(jù)。

1.2.3植被分布

1∶100萬(wàn)植被分布數(shù)據(jù)來(lái)源于國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)共享服務(wù)平臺(tái)(http:∥www.geodata.cn/)。大興安嶺植被類型主要為落葉松林、樟子松林、白樺林和山楊林等，植被數(shù)據(jù)可用于物種痕跡點(diǎn)的確定和物種潛在分布的模擬。

1.3 保護(hù)地空間數(shù)據(jù)

通過(guò)生態(tài)環(huán)境部網(wǎng)站(http:∥www.mee.gov.cn/)得到大興安嶺地區(qū)63個(gè)自然保護(hù)區(qū)名錄。共獲得55個(gè)大興安嶺地區(qū)保護(hù)地功能區(qū)劃(分區(qū))或邊界圖，其中14個(gè)自然保護(hù)區(qū)功能區(qū)劃圖來(lái)源于中國(guó)林業(yè)網(wǎng)(http:∥www.forestry.gov.cn/)，23個(gè)自然保護(hù)區(qū)邊界圖從各保護(hù)地管理局獲取，18個(gè)大興安嶺地區(qū)國(guó)家濕地公園分區(qū)圖從濕地保護(hù)管理中心獲取。另有3個(gè)省級(jí)和5個(gè)縣級(jí)及縣級(jí)以下自然保護(hù)區(qū)未能獲取邊界圖，根據(jù)地理位置標(biāo)出(圖1)。

1.4 模型模擬

1.4.1模型運(yùn)行

采用MaxEnt 3.4.1軟件(http:∥biodiversityinformatics.amnh.org/open_source/maxent/)模擬紫貂潛在適宜生境。采用ArcGIS 10.0和DIVA-GIS 7.5.0(http:∥www.diva-gis.org/)軟件，分別建立4個(gè)氣候變量、3個(gè)地形變量和1個(gè)植被分布變量柵格文件，統(tǒng)一所有圖層坐標(biāo)系、邊界和分辨率(5.0弧分)，將紫貂痕跡點(diǎn)和環(huán)境變量數(shù)據(jù)輸入MaxEnt模型。繪制環(huán)境變量響應(yīng)曲線，選擇輸出數(shù)值為邏輯斯蒂值(Logistic)。隨機(jī)選擇75%的物種分布點(diǎn)作為訓(xùn)練組，用以獲取建立最大熵模型所需參數(shù)；剩余25%的物種分布點(diǎn)作為測(cè)試組，用以驗(yàn)證模型準(zhǔn)確性[8]。重復(fù)次數(shù)為4，最大迭代次數(shù)為500，其他參數(shù)為默認(rèn)值。

1.4.2模型準(zhǔn)確性評(píng)估

當(dāng)物種痕跡點(diǎn)有限時(shí)，MaxEnt模型對(duì)一定區(qū)域內(nèi)特定物種分布區(qū)域模擬結(jié)果較優(yōu)[9]。受試者工作特征曲線(receiver operator characteristic curves，ROC)下的面積，即AUC值，被廣泛用于物種分布模型(species distribution model，SDM)模擬能力的評(píng)估。AUC值越高，模型模擬結(jié)果則越準(zhǔn)確。按AUC值可將模型評(píng)估結(jié)果分為5個(gè)等級(jí)：預(yù)測(cè)較差(0.5～0.6)、預(yù)測(cè)一般(>0.6～0.7)、預(yù)測(cè)較準(zhǔn)確(>0.7～0.8)、預(yù)測(cè)很準(zhǔn)確(>0.8～0.9)和預(yù)測(cè)極準(zhǔn)確(>0.9～1.0)。當(dāng)以25%的物種分布點(diǎn)作為測(cè)試組時(shí)，要求AUC值>0.7，模型模擬的物種潛在適宜生境才能用于進(jìn)一步研究[21]。

1.4.3潛在適宜生境

MaxEnt模型輸出結(jié)果為連續(xù)物種分布適宜性指數(shù)(habitat suitability index，HSI)，范圍為[0，1]。采用ArcGIS 10.0軟件分別計(jì)算未來(lái)氣候情景下3種氣候模型模擬結(jié)果平均值，作為紫貂在該氣候情景下相應(yīng)時(shí)段的潛在分布情況。HSI值越大，生境適應(yīng)性就越好。HSI值≥0.5的區(qū)域?yàn)樽硝踹m宜生境，HSI值<0.5的區(qū)域?yàn)椴贿m宜生境。

1.5 保護(hù)規(guī)劃

根據(jù)潛在適宜生境數(shù)據(jù)采用Zonation 4.0軟件(https:∥www.helsinki.fi/en/r￣e￣s￣e￣a￣r￣c￣h￣g￣r￣o￣u￣p￣s/m￣e￣t￣a￣p￣o￣p￣u￣l￣a￣t￣ion-research-centre/software#section-14300/)計(jì)算各區(qū)域優(yōu)先等級(jí)，以確定在當(dāng)前和未來(lái)氣候條件下紫貂優(yōu)先保護(hù)區(qū)域。計(jì)算時(shí)采用核心區(qū)移除規(guī)則(core-area cell removal rule)，翹曲因子(warp factor)設(shè)為1，其他參數(shù)為默認(rèn)值。

根據(jù)“愛知目標(biāo)”中提出的到2020年實(shí)現(xiàn)陸地和內(nèi)陸水域保護(hù)地覆蓋率達(dá)17%的目標(biāo)[22]，將優(yōu)先等級(jí)≥0.83的優(yōu)先保護(hù)區(qū)域作為紫貂目標(biāo)優(yōu)先保護(hù)區(qū)域，將其與保護(hù)地進(jìn)行疊加分析，分析氣候變化下紫貂優(yōu)先保護(hù)區(qū)域受保護(hù)情況。

2 結(jié)果與分析

2.1 模型預(yù)測(cè)結(jié)果評(píng)估

通過(guò)查閱資料，共獲取紫貂痕跡點(diǎn)96個(gè)，MaxEnt模型模擬結(jié)果的AUC值為0.852，表明模擬結(jié)果準(zhǔn)確。由表2可知，對(duì)模型貢獻(xiàn)率較大的變量為最冷月最低溫(42.5%)、植被類型(31.5%)、年均溫(9.8%)和平均氣溫日較差(9.1%)，其他指標(biāo)貢獻(xiàn)率較小。分析貢獻(xiàn)率較大的最冷月最低溫和植被類型與紫貂適宜性指數(shù)的反應(yīng)曲線發(fā)現(xiàn)，紫貂適宜生境(HIS值≥0.5)中，最冷月最低溫應(yīng)低于-37 ℃，植被優(yōu)勢(shì)種為興安落葉松。

表2 環(huán)境變量及其貢獻(xiàn)率Table 2 Description and percent contribution of 8 environmental variables

2.2 紫貂潛在適宜生境

由圖2可知，在當(dāng)前氣候條件下，紫貂潛在適宜生境面積為56 043.44 km2，占研究區(qū)總面積的13.74%，集中分布于大興安嶺中部及西北部地區(qū)。由圖3～4可知，至2041—2060和2061—2080年，RCP 4.5氣候情景下紫貂潛在適宜生境面積占研究區(qū)總面積比例由13.74%分別降至6.31%和6.52%，RCP 8.5氣候情景下則分別降至為4.05%和0.02%。未來(lái)氣候變化對(duì)紫貂潛在適宜生境影響極大，紫貂潛在適宜生境總體呈縮減趨勢(shì)，RCP 8.5氣候情景下的縮減程度大于RCP 4.5氣候情景。潛在適宜生境減少區(qū)域主要分布于大興安嶺北部，即漠河縣和額爾古納市北部區(qū)域。位于根河市與呼瑪縣交界地帶和額爾古納市北部的潛在適宜生境在未來(lái)氣候情景下均保持穩(wěn)定，而位于根河市中部地區(qū)的潛在適宜生境呈現(xiàn)擴(kuò)張趨勢(shì)。

審圖號(hào)：GS(2019)3621號(hào)圖2 當(dāng)前氣候條件下大興安嶺地區(qū)紫貂潛在適宜生境分布Fig.2 Distribution of potential suitable habitats of Martes zibellina in Daxing′anling relative to current climate scenario

2.3 紫貂優(yōu)先保護(hù)區(qū)

大興安嶺地區(qū)紫貂潛在適宜生境優(yōu)先保護(hù)等級(jí)見圖5。由于RCP 8.5氣候情景下2061—2080年紫貂潛在適宜生境面積僅占研究區(qū)總面積的0.02%，因此未計(jì)算該條件下優(yōu)先保護(hù)區(qū)域。在當(dāng)前氣候條件下，大興安嶺地區(qū)紫貂優(yōu)先保護(hù)區(qū)域主要分布于額爾古納市東北部、漠河縣西南部、呼瑪縣西部和根河市東部，其中漠河縣西南部、呼瑪縣西部和根河市東部的交界地帶是我國(guó)東北地區(qū)最需要優(yōu)先保護(hù)的區(qū)域之一[23]。RCP 4.5氣候情景下，紫貂優(yōu)先保護(hù)區(qū)域主要分布于根河市中東部及北部。RCP 8.5氣候情景下，紫貂優(yōu)先保護(hù)區(qū)域主要分布于根河市中部。

將當(dāng)前氣候條件下紫貂目標(biāo)優(yōu)先保護(hù)區(qū)域與保護(hù)地疊加，得到當(dāng)前紫貂保護(hù)空缺(圖6)。在當(dāng)前氣候條件下紫貂目標(biāo)優(yōu)先保護(hù)區(qū)域保護(hù)率最大，達(dá)46.19%，紫貂保護(hù)空缺主要分布于漠河縣西南部與額爾古納市交界地帶(區(qū)域a)、根河市東北部與呼瑪縣交界地帶(區(qū)域b)。保護(hù)空缺區(qū)域主要包含落葉針葉林、落葉闊葉林和草本濕地3種土地覆被類型，沒(méi)有居住地和水田(圖6)。區(qū)域a位于內(nèi)蒙古烏瑪與阿魯自然保護(hù)區(qū)之間，且與兩個(gè)保護(hù)區(qū)相連接，但有交通用地穿插其中；區(qū)域b位于內(nèi)蒙古阿魯和汗馬以及黑龍江嶺峰和呼中4個(gè)自然保護(hù)區(qū)之間，但未與阿魯自然保護(hù)區(qū)相連接。

審圖號(hào)：GS(2019)3621號(hào)圖3 RCP 4.5氣候情景下大興安嶺地區(qū)紫貂潛在適宜生境變化Fig.3 Variation of potential suitable habitats of Martes zibellina in Daxing′anling relative to RCP 4.5 climate scenario

審圖號(hào)：GS(2019)3621號(hào)圖4 RCP 8.5氣候情景下大興安嶺地區(qū)紫貂潛在適宜生境變化Fig.4 Variation of potential suitable habitats of Martes zibellina in Daxing′anling relative to RCP 8.5 climate scenario

3 討論

3.1 氣候變化背景下紫貂的潛在分布區(qū)變化

筆者研究發(fā)現(xiàn)在RCP 4.5和RCP 8.5氣候情景下紫貂潛在適宜生境面積占大興安嶺地區(qū)總面積比例將減少7.21～13.72百分點(diǎn)，縮減區(qū)域集中分布于大興安嶺地區(qū)北部。

大興安嶺地區(qū)北部為興安落葉松集中分布區(qū)，興安落葉松分布區(qū)具有溫度低、氣溫年較差大、太陽(yáng)輻射和年積溫較小等特點(diǎn)[24]，RCP 4.5和RCP 8.5氣候情景下大興安嶺北部地區(qū)增溫最明顯[25]，這將導(dǎo)致該區(qū)域興安落葉松分布面積大幅減少。以往研究表明紫貂偏好棲息于以興安落葉松為優(yōu)勢(shì)種的成熟林和老林地[16]，而大興安嶺地區(qū)北部興安落葉松分布區(qū)域的縮減可能導(dǎo)致紫貂潛在適宜生境減少。

此外，有研究表明，2010—2016年紫貂分布范圍較1950年代縮減超過(guò)50%，這表明隨氣候變化，紫貂在大興安嶺的適宜生境呈一直減少的趨勢(shì)[7]。不同氣候變化情景對(duì)紫貂潛在適宜生境影響的程度也不同。RCP 8.5氣候情景下，紫貂潛在適宜生境面積減少程度大于RCP 4.5氣候情景。

審圖號(hào)：GS(2019)3621號(hào)A為當(dāng)前氣候條件；B和C分別為2041—2060和2061—2080年RCP 4.5氣候情景下；D為2061—2080年RCP 8.5氣候情景下。圖5 大興安嶺地區(qū)紫貂潛在適宜生境優(yōu)先保護(hù)等級(jí)Fig.5 Distribution of priority protection grade of potential suitable habitats of Martes zibellina in Daxing′anling

a為漠河縣西南部與額爾古納市交界地帶，b為根河市東北部與呼瑪縣交界地帶。審圖號(hào)：GS(2019)3621號(hào)圖6 當(dāng)前氣候條件下大興安嶺地區(qū)紫貂優(yōu)先保護(hù)區(qū)域和保護(hù)空缺Fig.6 Distribution of priority protection area and protection gaps of Martes zibellina in Daxing′anling relative to current climate scenario

一方面，這可能與不同氣候情景下最冷月最低溫不同有關(guān)，MaxEnt模型模擬結(jié)果顯示，紫貂潛在適宜生境受最冷月最低溫的影響極大，且其生境適宜性指數(shù)隨著該指標(biāo)的升高而降低。另一方面，氣候變暖程度的增加將加深物種的受威脅程度[26]，氣溫上升幅度的增加可能加大對(duì)紫貂分布的負(fù)面影響，RCP 8.5氣候情景下2081—2100年升溫幅度大于RCP 4.5氣候情景，紫貂可能難以應(yīng)對(duì)較大的變暖幅度，適宜生境縮減程度也將加劇。這與大興安嶺如駝鹿和黑嘴松雞等其他物種潛在適宜生境變化模擬結(jié)果較為一致[27-28]。此外，將增溫控制在1.5 ℃以下，動(dòng)植物的氣候避難所相比增溫2 ℃條件下將增加5.5%～14%[26]，可減輕變暖對(duì)物種遷移和擴(kuò)散的壓力。因此，物種分布模擬研究應(yīng)綜合考慮多種氣候變化情景，為生物多樣性保護(hù)工作提供更全面的依據(jù)。

3.2 紫貂的保護(hù)空缺

在當(dāng)前氣候條件下，紫貂目標(biāo)優(yōu)先保護(hù)區(qū)域僅有總面積的46.19%受到保護(hù)地保護(hù)。截至2018年，大興安嶺地區(qū)已經(jīng)建立面積為45 704.10 km2的保護(hù)地，占大興安嶺地區(qū)總面積的24.08%，但保護(hù)空缺仍然存在。紫貂保護(hù)空缺區(qū)域主要為根河市中東部、額爾古納市與漠河縣交界地帶，保護(hù)空缺區(qū)域所含斑塊與內(nèi)蒙古烏瑪和阿魯自然保護(hù)區(qū)均有連接，但部分斑塊之間有道路分布，這可能成為紫貂遷移和活動(dòng)的障礙，建議加強(qiáng)對(duì)紫貂運(yùn)動(dòng)路徑的監(jiān)測(cè)和跟蹤，在阻隔紫貂運(yùn)動(dòng)的道路區(qū)建設(shè)生態(tài)廊道；建議適當(dāng)增加保護(hù)地范圍，特別是大興安嶺北部烏瑪國(guó)家級(jí)自然保護(hù)區(qū)與根河市中東部間的區(qū)域；加強(qiáng)現(xiàn)有保護(hù)地保護(hù)能力建設(shè)，提高對(duì)保護(hù)地外圍區(qū)域的關(guān)注。對(duì)于保護(hù)空缺所在地——額爾古納市、根河市、漠河縣和呼瑪縣，各縣(市)相關(guān)部門應(yīng)當(dāng)加強(qiáng)交流與合作，促進(jìn)紫貂綜合保護(hù)。未來(lái)在開展生物多樣性保護(hù)工作中，應(yīng)綜合考慮大興安嶺地區(qū)優(yōu)先保護(hù)區(qū)域變化趨勢(shì)及當(dāng)前保護(hù)地建設(shè)狀況，適當(dāng)調(diào)整現(xiàn)有保護(hù)地體系，為瀕危物種的有效保護(hù)提供保障。

3.3 模型評(píng)估

筆者研究所采用的MaxEnt模型簡(jiǎn)單易行，依靠物種痕跡點(diǎn)數(shù)據(jù)及環(huán)境變量就可開展模擬分析，通過(guò)痕跡點(diǎn)對(duì)應(yīng)的環(huán)境變量預(yù)測(cè)物種的潛在分布；在眾多物種分布模型中，MaxEnt模型預(yù)測(cè)準(zhǔn)確性較高[8]，被廣泛應(yīng)用于物種潛在分布模擬與評(píng)估研究。ZHANG等[7]模擬分析了1950—2016年紫貂在我國(guó)東北地區(qū)分布變化，發(fā)現(xiàn)2000—2016年紫貂在大興安嶺地區(qū)潛在分布區(qū)域位于大興安嶺中部及北部地區(qū)，筆者研究中紫貂當(dāng)前潛在適宜生境分布與之一致。然而，全球氣候變化不僅體現(xiàn)在氣溫和降水條件的改變，還表現(xiàn)為干旱、極端降水和熱浪等極端氣候事件增多，極端氣候事件是物種調(diào)節(jié)種群動(dòng)態(tài)和分布的關(guān)鍵驅(qū)動(dòng)因子[29]，這使得氣候變暖對(duì)物種的影響可能比模型預(yù)測(cè)結(jié)果更加強(qiáng)烈。同時(shí)，紫貂分布還會(huì)受到人類活動(dòng)影響[13]。未來(lái)研究需要進(jìn)一步考慮人類活動(dòng)和極端氣候事件對(duì)紫貂分布的影響。