基于紅外相機技術和MaxEnt模型的黑麂(Muntiacus crinifrons)活動節律分析和潛在適生區預測

陳智強,趙增輝, 王遠飛, 韋 力,丁國驊, 林植華,*

1 麗水學院生態學院, 麗水 323000 2 浙江農林大學動物科技學院, 臨安 311300

全球氣候變化及人為干擾下的棲息地破壞,使生物多樣性面臨前所未有的威脅[1]。盡管有關氣候變化導致物種滅絕的證據相對有限,但在未來的幾十年中,氣候變化可能對全球范圍內的生物多樣性造成巨大威脅[2]。已被證明,其可導致物種的棲息地向高緯度和高海拔地區遷移以及生物季節活動時間的改變[3-4]。而隨著物種向高海拔地區遷移,適宜其棲息的土地面積的數量、質量下降以及當前和未來合適海拔帶之間的垂直間隙的擴大可能妨礙動物種群的擴散與交流,進而導致種群規模下降,最終導致物種滅絕[5-6]。山地物種的棲息地選擇具有高度的地方特異性[7-8]、狹窄溫度帶[6, 9-10]以及狹窄的海拔范圍[6, 11]。因此,了解山地動物的活動節律,預測其適生區分布,探究棲息地選擇的關鍵生態因子,有利于開展山地動物的保護與管理。

近些年,隨著紅外線觸發自動數碼相機陷阱技術(Infrared triggered cameras technology, ITCT, 以下簡稱紅外相機技術)的開發,越來越多的紅外相機技術運用于鳥類和獸類的調查中[12-13]。該技術與傳統的調查方法相比,具有可日夜連續監測、無損傷等優點,已成研究野生動物,尤其是珍稀瀕危野生動物的重要方法[14]。最大熵模型(MaxEnt)具有較好的物種分布預測能力[15-16],且被應用于珍稀瀕危野生動物的生境適宜性評價[17- 20],但對于種群稀少且行蹤隱蔽的動物難以有效評價。通過結合紅外相機技術所獲得的物種分布的準確經緯度信息,可提高生物適生區預測結果的準確性。

黑麂(Muntiacuscrinifrons)隸屬于偶蹄目鹿科麂屬,為典型的山地物種,是我國特有Ⅰ級重點保護動物,被《中國脊椎動物紅色名錄》列為瀕危[21]。黑麂棲息地的垂直分布較窄,主要棲息于1000m左右的山區常綠闊葉林以及常綠落葉混交林中的芒叢或灌木林中[22-23]。有關黑麂的研究主要涉及物種資源及資源量調查[22- 26]、生境選擇[27- 30]、食性[31-32]、繁殖[33- 36]、遺傳多樣性[37- 41]等內容。目前,采用紅外相機技術已在浙贛皖三省的9個自然保護區中陸續發現黑麂的活動[42- 51]。然而,有關黑麂適生區預測的研究多集中自然保護區內[30, 52-53],僅雷軍成等[54]在全國范圍內研究氣候變化對其分布區的潛在影響,側重于氣候變化對2050和2070年時黑麂適宜生境的影響,未對其現存潛在棲息地進行預測[54]。本研究以浙江遂昌牛頭山林場為研究區域,采用紅外相機技術對該區域內黑麂資源進行調查,分析黑麂的活動節律與環境溫度之間的關系。同時,采用MaxEnt模型對于浙閩贛皖四省的黑麂潛在棲息地進行預測,以期為現有黑麂自然種群資源的全面保護以及保護區的劃分和管理提供科學依據。

1 研究區域與方法

1.1 研究區域

牛頭山屬武夷山系仙霞嶺山脈分支,地處浙江省麗水市遂昌縣東北部和金華市武義縣西南部。遂昌牛頭山林場擁有鄭崗嶺、大西坑、荷花地、天堂等7個林區,面積約為7000 hm2,絕大多數山地為中山區,少數山地為低山丘陵,相對高差1300 m,主峰海拔1560 m。該區域雨量充沛,氣候溫暖濕潤,屬于中亞熱帶濕潤季風區,一年中氣候有明顯的季節變化。年平均氣溫為10—16 ℃,年總降水量為1800—2400 mm[55]。林場全境植物群落多樣,區系組成復雜,包括常綠闊葉林、落葉—常綠闊葉林混交林、落葉闊葉林、竹林等多種植被類型。分布于該區域的瀕危珍稀動物主要有黃腹角雉(Tragopancaboti)、白頸長尾雉(Syrmaticusellioti)、白鷴(Lophuranycthemera)、黑麂、鬣羚(Capricornismilneedwardsii)和獼猴(Macacamulatta)等(數據未發表)。

1.2 紅外相機數據獲取

根據遂昌牛頭山林場的地形、植被類型、獸類活動路徑等多種因素,采用樣線法進行紅外相機布設,兩相機間的布設距離不小于200 m。2018年1月在主峰下和桃花嶺樣區各選擇12個位點,西肩嶺和大西坑樣區各9個位點,高山濕地樣區9個位點,九盤里樣區6個位點,共計57個位點(圖1),監測時長為1年。監測位點涵蓋竹林、針葉林、常綠闊葉林、灌木林和針闊混交林5種植被類型,海拔寬度為648—1413 m。

圖1 已報道黑麂的自然保護區位點及遂昌牛頭山紅外相機位點Fig.1 Map of the nature reserve site of M. crinifrons have been reported and infrared cameras installation sites for studying M. crinifrons in Suichang Niutou mountain, Zhejiang ProvinceGNJ: 牯牛降國家級自然保護區; JLF: 九龍峰省級自然保護區; AQLF: 安徽清涼峰國家級自然保護區; ZQLF: 浙江清涼峰國家級自然保護區; GTS: 古田山國家級自然保護區; JLS: 九龍山國家級自然保護區; FYS: 鳳陽山—百山祖國家級自然保護區; WYL:烏巖嶺國家級自然保護區; JXWY: 江西武夷山國家級自然保護區; NTS: 遂昌牛頭山。A: 西肩嶺; B: 大西坑; C: 高山濕地; D: 九盤里; E: 主峰下; F: 桃花嶺

每個位點布設1臺 Ltlacorn 感應紅外相機。相機多固定于具有明顯動物活動痕跡,如糞便、獸徑旁的大型喬灌木上,距地面約0.5 m。紅外相機統一設置,包括拍攝模式(照片+視頻)、連拍(3張)、時間間隔(1 s)和靈敏度(中)等；詳細記錄相機布設地點的地理坐標、海拔、植被類型等。監測期間,及時檢查與更新32G SD卡和電池。

1.3 紅外相機數據讀取與分析

紅外相機數據包括照片和視頻。通過照片分析及視頻輔助對所拍攝到的物種進行鑒定。根據《浙江動物志: 獸類》[56]的描述進行黑麂照片判定,提取獨立有效照片,以降低同一物種的照片相關性。獨立有效照片的判定標準:(1)相同或不同物種的不同個體連續照片和視頻；(2)相同物種的連續照片和視頻之間時間間隔大于 30 min；(3)相同物種的不連續照片或者視頻。符合以上任意一條標準,視為獨立有效照片[57-58]。每張獨立有效照片視為一次有效探測,并對拍攝時間、環境溫度和相機位點等信息進行匯總。

1.4 文中所涉及的術語定義

i:記錄時間點數目,時間間隔為1 h,全天24 h拍攝黑麂的日活動模式；ti:第i個數據點的時間值,如6:00—7:00時間值為6,11:00—12:00時間值為11；N:觀察期間相機工作日天數。1臺相機連續正常工作24 h記為1個相機工作日；A:觀察期間樣地內黑麂的總有效探測次數；pi:在ti時刻拍攝到的黑麂有效探測次數；Ti:在i℃溫度拍攝到的黑麂有效探測次數；Si:在某個季節下的黑麂有效探測次數(3—5月定義為春季；6—8月定義為夏季；9—11月定義為秋季；12—翌年2月定義為冬季)；Ei:在某個海拔區段內的黑麂有效探測次數(將相機海拔區段分為600—800 m、800—1000 m、1000—1200 m和>1200 m)；Vi:在某個植被類型下的黑麂有效探測次數；ri:ti時刻的相對活動強度,定義為ti時刻內拍攝到的黑麂有效檢測次數占樣地內A的比例即ri=pi/A。

相對多度指數(relative qbundance index, RAI),用以代表動物種群的相對數量, 指數大小與動物種群相對數量呈正相關[59],公式如下:RAI=(A/N)×1000

晝行性指數(diurnal-nocturnal index, DNI)用于判別黑麂的習性屬于晝行性還是夜行性[60]。若DNI數值接近13/24(約為0.54),表示黑麂活動沒有明顯的晝夜差別；DNI>13/24,表明黑麂活動以晝行性為主；DNI<13/24,表明黑麂活動夜行性為主。公式如下:

時間段活動強度指數(time-period activity intensity index, TAII):衡量黑麂日活動節律差異的指標[61]。

TAII=(pi/A)×100

季節活動強度指數(seasonal activity intensity index, SAII):衡量黑麂季節活動節律差異的指標。

SAII=(Si/A)×100

溫度活動強度指數(temperature activity intensity index, TEAII):分析黑麂對環境溫度的選擇[61]。

TEAII=(Ti/A)×100

海拔活動強度指數(elevation activity intensity index, EAII):分析黑麂對海拔的選擇。

EAII=(Ei/A)×100

植被類型活動強度指數(vegetational activity intensity index, VAII):分析黑麂對植被類型的選擇。

VAII=(Vi/A)×100

1.5 環境溫度的讀取

環境溫度受海拔影響且在山地小氣候影響下,每臺相機布設點的氣溫在同一時刻下也可能存在差異。Welbourne等[62]研究顯示,紅外相機的溫度傳感裝置較為敏感,其工作原理是根據被拍攝對象的熱輻射與周圍環境的熱量差異而觸發相機拍攝。因此,紅外相機所記錄的環境溫度具有一定的參考意義。

1.6 黑麂分布區預測

1.6.1黑麂分布點信息

通過遂昌牛頭山紅外相機監測、文獻檢索搜集到黑麂現代分布記錄28條。參照王茹琳等[63]的方法確定標本地理坐標并對全部分布記錄進行篩選,最終篩選獲得16個分布點。

1.6.2環境變量數據的獲取與選擇

本文共選取環境因子20個,包括氣候因子19個(BIO1—BIO19)、地形因子1個(Altitude)。氣候與地形數據均來自于WorldClim環境數據庫(http://www.Worldclim.org/),以當代時間段(1970—2000年)的生物環境因子均值為結果,空間分辨率為30″(約等于1 km2)。用DIVA—GIS軟件提取16個分布點的19個當代氣候變量數值與海拔數值。為降低因子間多重共線性對于模型的干擾,本文參考Yang等[64]提出的方法,運用多重共線性分析來檢驗環境因子間的相關性,對于兩因子之間相關性絕對值|r|>0.85,只選用一個因子進入模型。根據因子相關性分析結果,最終篩選6個環境變量,包括ALT(海拔)、BIO3(等溫性)、BIO6(最冷月份最高溫)、BIO9(最干季度均溫)、BIO12(年均降雨量)、BIO14(最干月降雨量)。

1.6.3MaxEnt模型預測與精度驗證

采用MaxEnt 3. 3. 3 k軟件(htp://www.cs.princeton.edu/-schapire /MaxEnt /)對黑麂的潛在分布區進行預測。MaxEnt模型的運行屬性設置為:模型重復運行5次；每一次運行都使用分布樣點中75%的隨機樣本為訓練集,剩余25%的樣本為測試集；使用“Do jackknife to measure variable importance”功能對各氣候變量進行Jackknife檢驗以評價其重要性；使用“Create response curves”功能創建各氣候變量的響應曲線以評估其對模型預測結果的影響；采用“cross validate”的重復運行類型；以Logistic格式輸出結果。構建黑麂分布模型進行模擬,以5次模擬結果的平均值作為最終結果[65]。

AUC因其不易受到物種分布率與閾值的影響,評價精度穩定,常被用于物種分布模型的精度驗證[66]。本文MaxEnt 模型采用受試者工作特征曲線(receiver operating characteristic curve, ROC曲線)分析法驗證模型精度。ROC曲線下的面積為AUC值,AUC 值越大,表示環境變量與預測物種地理分布模型間的相關性越大,越容易分辨該物種有無分布,預測效果也就越好。評價標準為:AUC值為0.5—0.6,失敗；0.6—0.7,較差；0.7—0. 8,一般；0.8—0.9,好；0.9—1.0,非常好[67]。

1.6.4適生區劃分

利用ArcGIS 10.0的轉換工具將5次模型平均預測結果的ASCII編碼文件結果轉化為柵格格式按照人工分級法進行重分類,根據黑麂的生境適宜性分成4個等級,即非適宜區(<0.3)、低適宜區(0.3—0.5)、中適宜區(0.5—0.7)、高適宜區(>0.7),最終得到黑麂在浙閩贛皖四省的潛在適生區空間分布圖。同時,利用ArcGIS 10.0軟件的空間統計功能,計算不同潛在適生區面積及比例。

1.7 數據分析

全部數據用SPSS 18.0進行分析。用G檢驗和線性相關分析對相應數據進行統計分析。所有統計值用平均值±標準誤表示,顯著性水平設置為0.05。

2 結果

2.1 黑麂分布區及遂昌牛頭山紅外相機拍攝概況

結合已報道的黑麂野外種群相關資料,其主要棲息地位于北緯27.5°—30.2°,東經117.4°—119.7°范圍內(圖1)。本研究中,遂昌牛頭山的紅外相機有效工作日為13655個；在月份上,5月(1410)最多,12月(900)最少；在季節上,春季(4069個)最多,冬季(2971)最少。紅外相機視頻與照片文件約 89176份,從中獲取到4個樣區的黑麂活動監測數據,總有效探測次數為26次。4個樣區中監測到黑麂的相機布設平均海拔范圍為668—1356 m,植被類型包括針葉林、常綠闊葉林、灌木林和竹林,有效監測平均次數范圍為1—4次；各樣區RAI指數在1.27—3.85范圍內(表1)。另外,2018年6月與7月在桃花嶺樣區3個不同位點均記錄到成對黑麂活動。圖2為本研究中紅外相機拍攝到4張具有代表性的黑麂照片。

表1 遂昌牛頭山4個監測到黑麂的樣區情況

2.2 日活動節律

本研究中,黑麂的DNI指數為0.96。紅外相機在4個樣區監測到黑麂的日活動發生在4:00—20:00期間,且活動頻次在不同時間段存在顯著差異(G=26.45,df=16,P<0.05)。TAII指數在8：00—10:00和16:00—18:00時最高,占總活動頻次的55.62%(圖3)。

不同季節的TAII指數在同一時間段內并不一致。春季TAII指數在8:00—10:00和14:00—16:00時出現峰值(圖3)；夏季在4:00—6:00、8:00—10:00和16:00—18:00時出現峰值(圖3)；秋季則在6:00—10:00、12:00—14:00和16:00—18:00時出現峰值(圖3)。春夏季的TAII指數均在最后時間段內最大,而秋季的TAII指數在三個時間段內差異不大(圖3)。

圖3 遂昌牛頭山黑麂全年、春季、夏季和秋季的日活動節律Fig.3 Daily activity rhythm of M. crinifrons during one year, spring, summer and fall in Suichang Niutou Mountain, Zhejiang Province

2.3 季節活動節律

黑麂的季節活動頻次在不同季節間存在顯著差異(G=20.9,df=3,P<0.001),SAII指數在夏季最高,春季次之,秋季最低,冬季未監測到黑麂的活動(圖4)。環境溫度在不同季節間存在顯著差異(F3, 27=42.8,P<0.001),進一步用Tukey多重比較,兩兩季節間均存在顯著差異(allP<0.03),環境溫度從高到低依次為夏季、春季、秋季和冬季。SAII指數與季節平均溫度呈顯著正相關(r=0.962,F1, 2=24.6,P<0.05)(圖4)。

圖4 遂昌牛頭山黑麂季節活動強度指數與季節溫度情況. 圖中給出兩者關系圖和回歸方程、回歸系數Fig.4 Seasonal activity intensity index (SAII) of M. crinifrons and seasonal temperature in Suichang Niutou Mountain, Zhejiang Province. The relationship of SAII and seasonal temperature and regression equation and regression coefficient are given in the Figure

2.4 環境溫度、海拔和植被類型的選擇

本研究中,監測到黑麂活動的平均環境溫度為(21±1.22) ℃(范圍為(6.5—32.8) ℃)。以每5 ℃為間隔,黑麂的活動頻次在不同環境溫度區間存在顯著差異(G=16.36,df=5,P<0.01),TEAII指數呈偏峰分布,在21—25 ℃范圍內TEAII指數高達38.46%(圖5)；在不同海拔區段間活動頻次存在顯著差異(G=16.5,df=3,P<0.001),EAII指數在1000—1200 m海拔最高(圖5)；在不同植被類型下活動頻次存在顯著差異(G=16.8,df=3,P<0.001),VAII指數在針葉林時最高(圖5)。

圖5 遂昌牛頭山黑麂的溫度、海拔和植被類型活動強度指數情況Fig.5 Temperature, elevational and vegetational activity intensity index of M. crinifrons in Suichang Niutou Mountain, Zhejiang Province

2.5 潛在適生區預測

圖6為浙閩贛皖四省黑麂的棲息地適應性分布圖,其潛在適生區主要有浙贛皖潛在分布區(黃山-天目山-白際山-懷玉山-千里崗)、浙贛潛在分布區(武夷山北段-仙霞嶺)、洞宮山潛在分布區和括蒼山潛在分布區等4個區域。ROC曲線評價結果顯示MaxEnt模型的預測結果達到優秀水平(AUC平均值±標準差=0.976±0.015; 圖6),說明黑麂潛在適生區預測結果可靠且精度較高。

黑麂適生區面積統計表明,黑麂在浙閩贛皖四省的潛在適生區面積約為25980. 62 km2,約占四省總面積的4.87%。其中,高適宜區面積約為3700.73 km2,所占四省總面積的比重約為0.69%；中適宜區面積約為7552.51 km2,所占四省總面積的比重約為1.41%；低適宜區面積約為14727.39 km2,所占四省總面積的比重最大,約為2.76%。

從MaxEnt模型刀切法分析結果可以看出(圖6):使用單環境因子建模時,BIO12、ALT、BIO14、BIO6對模型增益效果最為明顯,表明這4個環境因子的貢獻較大；而BIO3、BIO9對模型增益效果不明顯,表明這2個環境因子的貢獻較小。

圖6 黑麂在浙閩贛皖的潛在適生區分布情況、MaxEnt模型對黑麂預測結果的ROC曲線驗證及基于刀切法的環境變量重要性分析Fig.6 Potential suitable habitats of M. crinifrons in four provinces, ROC curve verification of prediction process for M. crinifrons by Maxent model and importance analysis of environmental factors by Jackknife methodALT:海拔Altitude;BJ03、BJ06、BJ09、BJ12、BJ14分別代表等溫性、最冷月份最高溫、最干季度均溫、年均降雨量、最干月降雨量

MaxEnt模型運行結果中包含生態因子響應曲線,即存在概率與生態因子閾值范圍的關系(圖7)。根據黑麂的響應曲線,劃分6種環境因子的適宜值范圍(存在概率大于0.5的范圍)。BIO12響應曲線顯示,年均降雨量在大于1750 mm,為黑麂分布的適宜降雨量范圍,且降雨量在2800—3200 mm時黑麂分布概率最大。ALT響應曲線顯示,海拔小于2000 m時,分布概率隨著海拔的增加而增加,且在大于750 m時,其存在概率大于0.5。BIO14響應曲線顯示,黑麂的適宜范圍在40—140 mm。BIO6響應曲線顯示,黑麂的適宜范圍在-3—1 ℃。BIO3響應曲線顯示,黑麂的適宜范圍在19.77—26。BIO9響應曲線顯示,黑麂的適宜范圍在-3—7 ℃。

圖7 海拔、BIO3、BIO6、BIO9、BIO12和BIO14等6個生態因子響應曲線Fig.7 Responses curve of six ecological factors including altitude, BIO3, BIO6, BIO9, BIO12 and BIO14

3 討論

3.1 遂昌牛頭山林場黑麂的種群現狀

根據RAI指數大小比較可知,在監測到黑麂的4個樣區中,以桃花嶺樣區的種群數量最多,接近安徽清涼峰國家級自然保護區的種群數量[49],甚至超過浙江清涼峰國家級自然保護區的種群數量[43]。2018年6月與7月,布設在桃花嶺樣區的3個不同紅外相機均記錄到兩只黑麂成對活動。這表明,在遂昌牛頭山林場存在一定數量的黑麂種群。但遂昌牛頭山林場屬于國有林場,與其他已報道黑麂活動的9個自然保護區相比,保護力度較小、人為干擾大,黑麂種群受威脅程度大。

3.2 遂昌牛頭山林場黑麂的活動溫度選擇

環境溫度是影響動物活動節律的關鍵氣候因子[68]。哺乳動物可通過自身體內氧化代謝產熱在較大的環境溫度范圍內維持高而相對恒定的體溫[69]。而這種恒定體溫的維持是以高能量消耗為代價的[70]。當環境溫度適宜時,動物體內溫度與體外溫度的差距較少,體溫調節所需能耗大幅降低[71]。研究結果顯示,在一定環境溫度范圍內(5—30 ℃),黑麂SAII指數與環境溫度呈正相關關系,且21—25 ℃的溫度范圍為黑麂的活動高峰。因此,我們推測21—25 ℃是黑麂活動的適宜環境溫度范圍,在該適溫區內黑麂的體溫調節能耗較低。

3.3 遂昌牛頭山林場黑麂的活動節律

本研究通過計算黑麂的DNI指數與TAII指數,推斷黑麂為以晨昏型為主的晝行性動物,這與章書聲等[72]的結果一致。一天之內,黑麂一般有2或3個活動高峰期,夏秋兩季為3個,春季為2個。在古田山自然保護區中,黑麂的日活動高峰期為2個[72],這可能與兩地的氣候、海拔、主要植被類型等環境因子的差異有關。夏季由于日出時間較早,其最早的活動高峰在4:00—8:00內出現；同時,夏秋兩季的晚活動高峰較春季有時間滯后。這種活動高峰的季節變化可能為吻合日出日落的季節變化。黑麂的晨昏型活動節律且活動高峰時間隨季節變化而變化的習性,可能是為了通過進食嫩葉上的露水來提供水分,進而降低搜尋水源時被捕食的風險[73]。

當外界環境條件(如溫度,光照等)變化時,動物會通過調整自身行為以形成特定條件下的活動節律模式[74]。本研究表明,溫度是影響遂昌牛頭山黑麂季節活動節律的重要因子。監測期間,研究樣區的冬季平均溫度為6.78 ℃,處于黑麂的非適宜活動溫度內,這也是本研究未監測到黑麂的重要原因之一。在夏季,黑麂的活動頻次占全年總活動頻次的53.85%,且除晨昏時間段外,其日活動節律也在8:00—10:00內出現高峰。黑麂喜好在夏季活動可能原因有(1)監測期間,研究樣區的夏季平均溫度為23.18 ℃,環境溫度適宜黑麂活動；(2)黑麂全年以灌木為主要食物,食性更接近嫩食者[32],夏季為多數灌木的生長旺盛期,嫩葉較多,食物資源豐富；(3)確保產仔期間及新生后代的安全是動物在長期進化上所面對的問題,黑麂雖為全年繁殖[33],但在具有更好的隱蔽條件的夏季進行繁殖更為安全[75-76]。

3.4 遂昌牛頭山林場黑麂的海拔與植被類型的選擇

黑麂偏好棲息在海拔為1000—1200 m之間的針葉林中,其次為常綠闊葉林。三尖杉(Cephalotaxusfartunei)、光葉菝葜(Smilaxglabra)、矩形鼠刺(Iteachinensis)、南五味子(Kadsuralongipedunculata)、馬銀花(Rhododendronovatum)5種黑麂主要食物資源[22, 32]主要分布于針葉林、常綠闊葉林與針闊混交林。同時,Althoff等[77]認為植被郁閉度越低,林下小生境溫度就越高。針葉林的郁閉度較低,利于棲息地環境積溫[27]。棲息在以中山區和種植經濟樹種占絕大多數的牛頭山林場中,黑麂的棲息地選擇偏好可以獲得一定的熱量、食物資源以及隱蔽條件,這與黑麂在九龍山國家級自然保護區內的棲息地選擇結果類似[27, 29]。

3.5 黑麂的分布區預測

棲息地破碎化和喪失是導致生物多樣性減少和物種滅絕的重要原因[78]。而有效開展保護工作的前提是對于受脅物種的潛在分布區進行預測,充分了解其現有生境情況[79]。本研究預測結果顯示,黑麂現有適生區十分狹小,僅占四省總面積的4.87%,且其棲息地呈島嶼狀分布于四省的部分山地丘陵。這與盛和林研究結果大體相符[22],但存在一定的差異,即在洞宮山為黑麂分布的第四個區域。

黑麂是中國特有種,目前預測僅存在浙閩贛皖四省。由于棲息地的破碎化和喪失、人為干擾以及氣候變化,黑麂種群面臨著滅絕的危險。因此,基于本研究的模型預測結果提出以下建議:(1)重點識別和建立各適宜生境之間的遷徙廊道。黑麂現有4個地理分布區域呈分隔狀態,分布區內部也存在一定的水系或地理阻隔。如在洞宮山潛在分布區和括蒼山潛在分布區,因水系阻隔,其為孤立的潛在分布區,且潛在分布區內部適生區呈破碎的斑塊狀,棲息于此的黑麂滅絕風險大。因此,要保護每個分布區域的高適生區,重點識別和建立各適宜生境之間的遷徙廊道,以加強不同分布區間的不同黑麂種群的基因交流。(2)加強對黑麂潛在分布區內的種群資源進行調查。黑麂分布區狹窄、適宜棲息地破碎化嚴重,例如千里崗地區為黑麂潛在高適生區,但未有黑麂種群資源相關的報道。應加強對黑麂潛在分布區內的種群資源進行調查,以充分了解其種群現狀,為在浙閩贛皖四省進一步開展全面保護措施提供科學依據。(3)建立黑麂保護網絡。黑麂高適生區內有存在諸多國家級自然保護區。馬頭山國家級自然保護區-江西武夷山國家級自然保護區-九龍山國家級自然保護區貫穿武夷山北段-仙霞嶺形成一個野生動物生境廊道,有利于黑麂的浙贛潛在分布區種群的擴散和基因交流。但適生區不只存在于現有的自然保護區,在保護區周邊、諸多林場也分布有黑麂的潛在適生區。因此,在黑麂適生區內的林業局、保護區管理局以及村委會之間應加強宣傳,建立起黑麂保護網絡。