高 惠, 滕麗微,2, 汪 洋, 王繼飛, 劉振生,2,*

1 東北林業(yè)大學(xué)野生動(dòng)物資源學(xué)院,哈爾濱 150040 2 國(guó)家林業(yè)局野生動(dòng)物保護(hù)學(xué)重點(diǎn)開(kāi)放實(shí)驗(yàn)室,哈爾濱 150040 3 寧夏賀蘭山國(guó)家級(jí)自然保護(hù)區(qū)管理局,銀川 750021

?

阿拉善馬鹿(Cervusalashanicus)生境適宜性評(píng)價(jià)

高 惠1, 滕麗微1,2, 汪 洋1, 王繼飛3, 劉振生1,2,*

1 東北林業(yè)大學(xué)野生動(dòng)物資源學(xué)院,哈爾濱 150040 2 國(guó)家林業(yè)局野生動(dòng)物保護(hù)學(xué)重點(diǎn)開(kāi)放實(shí)驗(yàn)室,哈爾濱 150040 3 寧夏賀蘭山國(guó)家級(jí)自然保護(hù)區(qū)管理局,銀川 750021

阿拉善馬鹿(Cervusalashanicus)目前僅分布于賀蘭山地區(qū),對(duì)該物種進(jìn)行生境適宜性的評(píng)價(jià)和分析是物種有效保護(hù)的前提和基礎(chǔ)。2013—2014年通過(guò)樣線調(diào)查及巡山資料查詢,確定阿拉善馬鹿出現(xiàn)位點(diǎn)86個(gè),結(jié)合13種環(huán)境變量數(shù)據(jù),利用最大熵(maximum entropy, MaxEnt)模型,并根據(jù)最大約登指數(shù)劃定適宜與不適宜生境區(qū),對(duì)賀蘭山地區(qū)阿拉善馬鹿的生境適宜性進(jìn)行評(píng)價(jià)。ROC曲線(receiver operating characteristic curve)檢測(cè)證明模型預(yù)測(cè)精度較高,研究結(jié)果表明：阿拉善馬鹿主要分布于賀蘭山東坡的中部和南部,以及西坡的中北部,適宜生境面積為667.87 km2,占研究區(qū)域面積的18.2%；礦區(qū)、坡度和海拔是影響阿拉善馬鹿分布的最主要環(huán)境變量,礦區(qū)對(duì)阿拉善馬鹿的影響最大,建議管理部門(mén)加大對(duì)此人為干擾的管控力度,控制和減少現(xiàn)有礦區(qū)的規(guī)模,以促進(jìn)該種群的發(fā)展。

阿拉善馬鹿；生境評(píng)價(jià)；最大熵模型；賀蘭山

生境適宜性評(píng)價(jià)是以動(dòng)物對(duì)生活環(huán)境的要求為基礎(chǔ),分析研究區(qū)域內(nèi)各生態(tài)因子的實(shí)際分布情況,從而明確研究范圍內(nèi)適合某動(dòng)物生存分布的區(qū)域[1]。對(duì)野生動(dòng)物進(jìn)行生境適宜性評(píng)價(jià)和分析對(duì)該物種的保護(hù)具有非常重要的作用。目前,利用3S(GIS,RS,GPS)技術(shù)結(jié)合生態(tài)位模型對(duì)生境適宜性進(jìn)行分析已越來(lái)越多的運(yùn)用到野生動(dòng)物生境評(píng)價(jià)之中[2- 5]。生態(tài)位模型的基本原理是根據(jù)每種生物特殊的生存環(huán)境,從目標(biāo)物種已知分布區(qū)出發(fā),利用數(shù)學(xué)模型歸納或模擬其生態(tài)位需求,將其投射到目標(biāo)地區(qū)即可得到目標(biāo)物種的適生區(qū)分布[6]。生態(tài)位模型中的最大熵(MaxEnt)模型是最具有代表性的物種分布模型之一[7- 9],近年來(lái)在國(guó)內(nèi)外物種地理分布的研究中得到廣泛應(yīng)用,且相較其他生態(tài)位模型表現(xiàn)出更好的預(yù)測(cè)能力[10- 12]。

阿拉善馬鹿(Cervusalashanicus)原為馬鹿(Cervuselaphus)的一個(gè)亞種,為國(guó)家Ⅱ級(jí)重點(diǎn)保護(hù)野生動(dòng)物[13],《中國(guó)哺乳動(dòng)物多樣性及地理分布》(2015)將此亞種提升為種[14],依據(jù)此名錄目前中國(guó)境內(nèi)的馬鹿分為阿拉善馬鹿、四川馬鹿(C.macneilli)、西藏馬鹿(C.wallichii)、東北馬鹿(C.xanthopygus)和塔里木馬鹿(C.yarkandensis)。阿拉善馬鹿目前僅分布于寧夏和內(nèi)蒙古交界的賀蘭山中段,是我國(guó)唯一幸存的該亞種的有效種群[15- 17],劉振生等[18-19]和駱?lè)f等[20]已對(duì)賀蘭山地區(qū)阿拉善馬鹿的生境選擇進(jìn)行了較為詳細(xì)的研究,然而該地區(qū)阿拉善馬鹿的生境適宜性尚未見(jiàn)分析報(bào)道,并且阿拉善馬鹿在賀蘭山地區(qū)的種群數(shù)量一直處于較低水平,在此背景下,利用MaxEnt模型對(duì)賀蘭山地區(qū)的阿拉善馬鹿生境進(jìn)行評(píng)價(jià),明確阿拉善馬鹿在此區(qū)域的主要分布地,分析分布區(qū)與主要環(huán)境變量之間的關(guān)系,能夠?yàn)樵撐锓N的合理管護(hù)、促進(jìn)種群發(fā)展及進(jìn)一步探討同域分布草食性動(dòng)物的共存機(jī)制提供科學(xué)依據(jù)。

1 研究地區(qū)概況

賀蘭山位于寧夏回族自治區(qū)和內(nèi)蒙古自治區(qū)交界處(38°21′—39°22′N, 105°49′—106°42′E),由內(nèi)蒙古賀蘭山國(guó)家級(jí)自然保護(hù)區(qū)(賀蘭山西坡)和寧夏賀蘭山國(guó)家級(jí)自然保護(hù)區(qū)(賀蘭山東坡)組成。海拔一般為2000—3000 m,山勢(shì)陡峭、高差懸殊；為典型的大陸性氣候,全年干旱少雨,年均氣溫-0.9℃,年平均降水420 mm[15]；植被依據(jù)海拔的上升呈現(xiàn)明顯的垂直分異,自下而上依次為山地草原帶(1400—1600 m),山地疏林草原帶(1600—2000 m),山地針葉林帶(1900—3000 m),亞高山灌叢和草甸帶(3000—3556 m)[21]。

2 研究方法

2.1 阿拉善馬鹿出現(xiàn)點(diǎn)信息采集及處理

圖1 研究區(qū)域的阿拉善馬鹿出現(xiàn)點(diǎn)圖Fig.1 The occurrence points of Cervus alashanicus in study area

自2013年4月—2014年4月,在賀蘭山分4個(gè)季節(jié)對(duì)阿拉善馬鹿進(jìn)行樣線法調(diào)查,樣線根據(jù)賀蘭山的地形特征和溝系走向進(jìn)行布設(shè),共設(shè)定樣線25條,自溝口到山脊方向覆蓋該區(qū)域所有生境類型。利用GPS(集思寶 G330)對(duì)觀察到的阿拉善馬鹿實(shí)體、新鮮糞便、足跡等出現(xiàn)點(diǎn)信息進(jìn)行定位,每個(gè)季節(jié)調(diào)查1次。此外,查閱賀蘭山保護(hù)區(qū)2013—2014年的野外巡護(hù)記錄資料,提取阿拉善馬鹿出現(xiàn)位點(diǎn)信息。通過(guò)野外調(diào)查和資料查詢共獲得91個(gè)阿拉善馬鹿發(fā)現(xiàn)位點(diǎn),將其導(dǎo)入到ArcGIS中,去除出現(xiàn)在研究區(qū)域外的錯(cuò)誤坐標(biāo)點(diǎn),最后得到86個(gè)阿拉善馬鹿出現(xiàn)點(diǎn)信息(圖1)。將剩余位點(diǎn)轉(zhuǎn)換成MaxEnt軟件要求的.csv格式,字段設(shè)置為物種名稱、經(jīng)度和緯度。

2.2 環(huán)境變量信息獲取及處理

用于阿拉善馬鹿生境評(píng)價(jià)的環(huán)境變量分為地形因子、植被因子、水源及人為干擾因子。地形因子包括坡度、坡向和海拔,從國(guó)際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)網(wǎng)站(http:// datamirror.csdb.cn)下載30 m分辨率的數(shù)字高程模型(DEM)提取得到；植被因子劃分為高山草甸、針葉林、闊葉林、灌叢、草地和荒漠6種類型,從經(jīng)解譯的2007年賀蘭山TM遙感影像中提取獲得；礦區(qū)、道路和護(hù)林點(diǎn)3種人為干擾因子及水源的距離柵格圖層從賀蘭山1∶50000的矢量化地圖中提取,經(jīng)ArcGIS的空間分析工具進(jìn)行緩沖分析得到。所有環(huán)境變量在ArcGIS 10中利用空間分析工具進(jìn)行相關(guān)性檢測(cè),去除相關(guān)性過(guò)高(大于0.8)的因子,將剩余變量進(jìn)行邊界范圍和坐標(biāo)系的一致化,柵格大小為30 m×30 m,并輸出為MaxEnt軟件要求的格式。

2.3 MaxEnt模型預(yù)測(cè)

MaxEnt 3.3.3軟件從http://www.cs.princeton.edu/—schapire/maxent/網(wǎng)站下載得到。依次將阿拉善馬鹿出現(xiàn)點(diǎn)數(shù)據(jù)和環(huán)境變量數(shù)據(jù)導(dǎo)入后,將各環(huán)境因子設(shè)定為連續(xù)型(Continuous)變量或離散型(Categorical)變量,設(shè)定75%的數(shù)據(jù)用于模型的建立,剩余25%的數(shù)據(jù)用于模型的檢測(cè),用jackknife檢驗(yàn)各環(huán)境變量的重要性,其余設(shè)置成默認(rèn)值,結(jié)果以Logistic格式輸出并導(dǎo)入到ArcGIS 10中做進(jìn)一步處理。

2.4 模型驗(yàn)證

MaxEnt模型帶有自檢驗(yàn)功能,利用受試者工作特征曲線(Receiver operating characteristic curve, ROC曲線)對(duì)模型進(jìn)行檢驗(yàn),根據(jù)ROC曲線與橫坐標(biāo)圍成的面積即AUC(Area under ROC Curve)值對(duì)模型的準(zhǔn)確度進(jìn)行判定[22-24]。AUC值越大,表明環(huán)境變量與模型的相關(guān)性越大,預(yù)測(cè)的效果越好,評(píng)定標(biāo)準(zhǔn)設(shè)為：AUC值為0.5—0.6,不及格；0.6—0.7,較差；0.7—0.8,一般；0.8—0.9,良好；0.9—1.0,優(yōu)秀[25-26]。

2.5 阿拉善馬鹿適宜生境評(píng)定

將模型輸出的概率分布圖導(dǎo)入到ArcGIS 10中,依據(jù)最大約登指數(shù)(Youden′s index)選擇概率切斷點(diǎn)對(duì)MaxEnt模型預(yù)測(cè)的生境適宜圖進(jìn)行重分類,約登指數(shù)用來(lái)判斷ROC曲線的臨界值點(diǎn),指數(shù)越大表明預(yù)測(cè)效果越好[27-29]。將大于概率切斷點(diǎn)的區(qū)域作為阿拉善馬鹿的適宜生境,小于概率切斷點(diǎn)的區(qū)域作為阿拉善馬鹿不適宜生境,最終得到賀蘭山地區(qū)阿拉善馬鹿的生境適宜性分布等級(jí)圖。

3 結(jié)果

3.1 模型預(yù)測(cè)效果檢驗(yàn)

ROC曲線評(píng)價(jià)結(jié)果為：訓(xùn)練集(Training data)與驗(yàn)證集(Test data)的AUC值分別為0.965和0.898,表明MaxEnt模型的預(yù)測(cè)結(jié)果較好(圖2),所建模型可以用于分析賀蘭山地區(qū)馬鹿的分布情況。

3.2 環(huán)境因子對(duì)阿拉善馬鹿分布的影響

經(jīng)相關(guān)性分析表明,13個(gè)環(huán)境變量之間的空間相關(guān)性均小于0.8,都可作為阿拉善馬鹿生境適宜性評(píng)價(jià)因子,從Jackknife檢驗(yàn)結(jié)果看出(圖3),距礦區(qū)距離、坡度和海拔是對(duì)模型預(yù)測(cè)增益最大的3種因子,意味著阿拉善馬鹿的分布對(duì)這幾種因子的變化較敏感。

圖2 阿拉善馬鹿生境評(píng)價(jià)結(jié)果的ROC曲線驗(yàn)證Fig.2 ROC curve verification of Cervus alashanicus in the Helan Mountains, China

圖3 刀切法測(cè)定環(huán)境變量的重要性Fig.3 Jackknife test of environmental variable importance

各環(huán)境變量對(duì) MaxEnt模型的貢獻(xiàn)率表明：礦區(qū)(30.4%)、坡度(26%)、海拔(10.8%)和坡向(8.4%)這4種環(huán)境因子對(duì)模型的累積貢獻(xiàn)率達(dá)75.6%,對(duì)阿拉善馬鹿分布影響較大；其次為護(hù)林點(diǎn)(5.7%)、草地(4.2%)、道路(4.2%)、灌叢(4%)、高山草甸(3.8%),闊葉林(1.3%)和水源(0.9%)；而針葉林和荒漠這2個(gè)變量的貢獻(xiàn)率非常小。

各環(huán)境變量對(duì)阿拉善馬鹿分布預(yù)測(cè)的響應(yīng)曲線表明：阿拉善馬鹿的生境適宜性隨距礦區(qū)和距護(hù)林點(diǎn)距離的增加而明顯增加,在1000 m范圍內(nèi),隨距道路距離的增加適宜性也有增加趨勢(shì),隨著距水源距離的增加以及在超過(guò)30°的坡度范圍后適宜性都在逐漸降低,而隨著海拔的增加,適宜度出現(xiàn)兩個(gè)峰值,分別為1500—2000 m和2800—3000 m。

3.3 阿拉善馬鹿生境適宜性分布

圖4 阿拉善馬鹿生境適宜性分布圖 Fig.4 The distribution of habitat suitability of Cervus alashanicus in Helan Mountains, China

根據(jù)約登指數(shù)為“靈敏度+特異度-1”的求算方法,在MaxEnt模型運(yùn)算結(jié)果中,選取約登指數(shù)最大值時(shí)的概率切斷點(diǎn)0.31對(duì)阿拉善馬鹿分布預(yù)測(cè)圖進(jìn)行重分類,最終得到賀蘭山阿拉善馬鹿生境適宜性分布圖(圖4)。其中,適宜生境面積為667.87 km2,不適宜生境面積為3008.10 km2,分別占研究區(qū)域總面積的18.2%和81.8%。從分布圖中可以看出,阿拉善馬鹿的適宜生境主要分布于寧夏賀蘭山國(guó)家級(jí)保護(hù)區(qū)的中部和南部,內(nèi)蒙古賀蘭山國(guó)家級(jí)自然保護(hù)區(qū)中北部,賀蘭山北部的大部分區(qū)域都屬于阿拉善馬鹿的不適宜生境。

4 討論

阿拉善馬鹿的適宜生境主要分布于賀蘭山東坡的中部和南部以及西坡的中北部,這一結(jié)果可能是由于賀蘭山不同區(qū)段降水與植被覆蓋度不同造成的：保護(hù)區(qū)北段降水量約180 mm,區(qū)內(nèi)荒漠化嚴(yán)重,加上礦區(qū)較多對(duì)稀疏的植被產(chǎn)生較大的破壞；中部降水量約400 mm,1995年開(kāi)始的援德項(xiàng)目使中段的植被保存完好；南段降水約300 mm,由于1999年以前礦石開(kāi)采活動(dòng)嚴(yán)重,只在高山區(qū)域有稀疏森林與中段相連[21,30]。北段降水量少造成植被生長(zhǎng)較差,加上礦區(qū)對(duì)阿拉善馬鹿生存造成較強(qiáng)的人為干擾,導(dǎo)致賀蘭山北部大部分地區(qū)成為阿拉善馬鹿不適宜分布區(qū)。

研究結(jié)果表明海拔和坡度這兩個(gè)地形因子對(duì)阿拉善馬鹿具有較大影響,在超過(guò)30°的坡度后適宜度逐漸降低,在海拔1500—2000 m和2800—3000 m處適宜度出現(xiàn)兩個(gè)峰值,駱?lè)f[31]對(duì)阿拉善馬鹿的研究發(fā)現(xiàn)坡度是其生境選擇的一個(gè)重要限制因子,阿拉善馬鹿通常選擇坡度較緩的地方,同時(shí)發(fā)現(xiàn)在賀蘭山地區(qū)該物種存在明顯的季節(jié)性遷徙現(xiàn)象,冬春季多活動(dòng)于低海拔地區(qū)的山地草原帶,夏秋季偏愛(ài)選擇高海拔地區(qū)的山地針葉林地帶,與本次研究結(jié)果相對(duì)應(yīng)。這可能是由于冬春季氣溫較低,阿拉善馬鹿偏向于低海拔地區(qū)活動(dòng),而夏秋季天氣炎熱,阿拉善馬鹿喜歡選擇相對(duì)涼爽的高海拔地區(qū)生活。吳文等利用MaxEnt模型對(duì)小興安嶺南麓的馬鹿(Cervuselaphus)冬季適宜性進(jìn)行評(píng)價(jià),結(jié)果表明地形因子中的坡向貢獻(xiàn)率最高,海拔的貢獻(xiàn)率很低,與本研究結(jié)果不同,可能是由于兩個(gè)研究區(qū)域氣溫、風(fēng)向和植被分布存在差異,小興安嶺地區(qū)的馬鹿對(duì)西北坡表現(xiàn)出顯著的回避,偏愛(ài)在低海拔地區(qū)活動(dòng),相比較而言阿拉善馬鹿受海拔的影響要大于坡向；小興安嶺地區(qū)馬鹿的生境適宜度隨坡度(0—36°)的增大而不斷降低,與本研究結(jié)果相似,說(shuō)明兩個(gè)地區(qū)的馬鹿都喜歡平緩的生境[27,32]。

阿拉善馬鹿僅分布于賀蘭山地區(qū),從保護(hù)區(qū)建成以來(lái),加大了對(duì)野生動(dòng)物的保護(hù)力度,阿拉善馬鹿種群雖有所恢復(fù)但數(shù)量仍然較少。本次研究發(fā)現(xiàn),礦區(qū)這一人為干擾對(duì)該物種生境適宜度影響最大,20世紀(jì)80年代,賀蘭山地區(qū)有較多的工礦企業(yè),員工人數(shù)眾多,其生產(chǎn)和生活活動(dòng)對(duì)野生動(dòng)物造成了嚴(yán)重的干擾,如頻繁的運(yùn)輸車輛、機(jī)器的噪音以及多種因素導(dǎo)致的環(huán)境污染等[21],目前,保護(hù)區(qū)內(nèi)很多礦區(qū)已經(jīng)關(guān)閉,但仍存在的礦區(qū)干擾依舊較大,在調(diào)查期間發(fā)現(xiàn)離礦區(qū)較近的區(qū)域路面和植被都被采煤產(chǎn)生的粉塵染成黑色,人為活動(dòng)對(duì)當(dāng)?shù)刂脖坏钠茐募由宪囕v來(lái)往等產(chǎn)生的噪音嚴(yán)重降低了阿拉善馬鹿生境的質(zhì)量,礦區(qū)產(chǎn)生的干擾強(qiáng)度大這一結(jié)果在賀蘭山巖羊生境評(píng)價(jià)研究中也有體現(xiàn)[33]。研究表明,礦產(chǎn)開(kāi)發(fā)能夠?qū)е戮坝^格局的變化,使礦區(qū)周圍植被的種類和生物量降低,加劇野生動(dòng)物生境的破碎化[34]。張國(guó)鋼和張正旺[35]研究發(fā)現(xiàn)鐵礦開(kāi)采對(duì)褐馬雞(Crossoptilonmantchuricum)產(chǎn)生了屏障作用,降低其生境的連接度,使褐馬雞很少出現(xiàn)在礦區(qū)周圍?？锵容x[36]的研究也表明礦區(qū)等人為干擾因素對(duì)蜂桶寨大熊貓(Ailuropodamelanoleuca)的棲息地景觀格局產(chǎn)生明顯影響。本次研究結(jié)果可使相關(guān)部門(mén)有針對(duì)性的加大阿拉善馬鹿適宜生境區(qū)域的控制力度,特別是加強(qiáng)對(duì)礦區(qū)的管控,減少礦區(qū)的規(guī)模,重視礦區(qū)周圍植被和環(huán)境的恢復(fù),降低礦區(qū)等人為干擾對(duì)阿拉善馬鹿和其他野生動(dòng)物的影響,為阿拉善馬鹿種群的保護(hù)和發(fā)展提供有利條件。

致謝：寧夏賀蘭山國(guó)家級(jí)自然保護(hù)區(qū)管理局胡天華高級(jí)工程師和內(nèi)蒙古賀蘭山國(guó)家級(jí)自然保護(hù)區(qū)管理局王兆錠工程師幫助野外調(diào)查和數(shù)據(jù)收集，特此致謝。

[1] 蘭慧. 內(nèi)蒙古烏拉特荒漠鵝喉羚生境適宜性評(píng)價(jià)[D]. 北京: 中國(guó)林業(yè)科學(xué)研究院, 2015.

[2] Kanagaraj R, Wiegand T, Kramer-Schadt S, Anwar M, Goyal S P. Assessing habitat suitability for tiger in the fragmented Terai Arc Landscape of India and Nepal. Ecography, 2011, 34(6): 970- 981.

[3] 王蕓, 趙鵬祥. 青木川自然保護(hù)區(qū)大熊貓生境評(píng)價(jià). 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(1): 206- 212.

[4] Martin J, Revilla E, Quenette P Y, Naves J, Allainé D, Swenson J E. Brown bear habitat suitability in the Pyrenees: transferability across sites and linking scales to make the most of scarce data. Journal of Applied Ecology, 2012, 49(3): 621- 631.

[5] 李佳, 李言闊, 繆瀘君, 謝光勇, 袁芳凱. 江西桃紅嶺國(guó)家級(jí)自然保護(hù)區(qū)梅花鹿生境適宜性評(píng)價(jià). 生態(tài)學(xué)報(bào), 2014, 34(5): 1274- 1283.

[6] 蔡靜蕓, 張明明, 粟海軍, 張海波. 生態(tài)位模型在物種生境選擇中的應(yīng)用研究. 經(jīng)濟(jì)動(dòng)物學(xué)報(bào), 2014, 18(1): 47- 52, 58.

[7] Phillips S J, Anderson R P, Schapire R E. Maximum entropy modeling of species geographic distributions. Ecological Modelling, 2006, 190(3/4): 231- 259.

[8] Phillips S J, Dudík M. Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography, 2008, 31(2): 161- 175.

[9] De Castro Pena J C, Kamino L H Y, Rodrigues M, Mariano-Neto E, De Siqueira M F. Assessing the conservation status of species with limited available data and disjunct distribution. Biological Conservation, 2014, 170: 130- 136.

[10] 吳文浩, 李明陽(yáng). 基于生態(tài)位模型的松材線蟲(chóng)潛在生境預(yù)測(cè)方法研究. 林業(yè)調(diào)查規(guī)劃, 2009, 34(5): 33- 38.

[11] 曹向鋒, 錢(qián)國(guó)良, 胡白石, 劉鳳權(quán). 采用生態(tài)位模型預(yù)測(cè)黃頂菊在中國(guó)的潛在適生區(qū). 應(yīng)用生態(tài)學(xué)報(bào), 2010, 21(12): 3063- 3069.

[12] Elith J, Leathwick J R. Species distribution models: ecological explanation and prediction across space and time. Annual Review of Ecology, Evolution, and Systematics, 2009, 40(1): 677- 697.

[13] 昶野, 張明明, 劉振生, 胡天華, 李志剛. 賀蘭山同域分布巖羊和馬鹿的夏季食性. 生態(tài)學(xué)報(bào), 2010, 30(6): 1486- 1493.

[14] 蔣志剛, 馬勇, 吳毅, 王應(yīng)祥, 周開(kāi)亞, 馮祚建, 劉少英, 羅振華, 李春旺. 中國(guó)哺乳動(dòng)物多樣性及地理分布. 北京: 科學(xué)出版社, 2015.

[15] 王小明, 李明, 唐紹祥, 劉志霄, 李元廣, 盛和林. 賀蘭山偶蹄類動(dòng)物資源及保護(hù)現(xiàn)狀研究. 動(dòng)物學(xué)雜志, 1999, 34(5): 26- 29.

[16] 盛和林. 中國(guó)鹿類動(dòng)物. 上海: 華東師范大學(xué)出版社, 1992: 1- 215.

[17] 李明, 王小明, 盛和林, 玉手英利, 增田隆一, 永田純子, 大泰司紀(jì)之. 馬鹿四個(gè)亞種的起源和遺傳分化研究. 動(dòng)物學(xué)研究, 1998, 19(3): 177- 183.

[18] 劉振生, 曹麗榮, 翟昊, 胡天華, 王小明. 賀蘭山區(qū)馬鹿對(duì)冬季生境的選擇性. 動(dòng)物學(xué)研究, 2004, 25(5): 403- 409.

[19] 劉振生, 張明明, 李志剛, 胡天華, 翟昊. 賀蘭山馬鹿冬季取食和臥息生境選擇. 獸類學(xué)報(bào), 2009, 29(2): 133- 141.

[20] 駱?lè)f, 張明明, 劉振生, 李志剛, 胡天華, 翟昊. 賀蘭山馬鹿冬春季生境的選擇. 生態(tài)學(xué)報(bào), 2009, 29(5): 2757- 2763.

[21] 狄維忠. 賀蘭山維管植物. 西安: 西北大學(xué)出版社, 1987.

[22] Vanagas G. Receiver operating characteristic curves and comparison of cardiac surgery risk stratification systems. Interactive Cardiovascular and Thoracic Surgery, 2004, 3(2): 319- 322.

[23] 王運(yùn)生, 謝丙炎, 萬(wàn)方浩, 肖啟明, 戴良英. ROC曲線分析在評(píng)價(jià)入侵物種分布模型中的應(yīng)用. 生物多樣性, 2007, 15(4): 365- 372.

[24] 李明陽(yáng), 席慶, 徐海根. 基于WEB數(shù)據(jù)庫(kù)和Maxent模型的白頭葉猴潛在生境評(píng)價(jià)研究. 林業(yè)資源管理, 2010, (4): 31- 36, 57.

[25] 齊增湘, 徐衛(wèi)華, 熊興耀, 歐陽(yáng)志云, 鄭華, 甘德欣. 基于MAXENT模型的秦嶺山系黑熊潛在生境評(píng)價(jià). 生物多樣性, 2011, 19(3): 343- 352.

[26] Araújo M B, Pearson R G, Thuiller W, Erhard M. Validation of species-climate impact models under climate change. Global Change Biology, 2005, 11(9): 1504- 1513.

[27] 吳文, 李月輝, 胡遠(yuǎn)滿, 陳龍, 李悅, 李澤鳴, 聶志文, 陳探. 小興安嶺南麓馬鹿冬季適宜生境評(píng)價(jià). 生物多樣性, 2016, 24(1): 20- 29.

[28] Jiménez-Valverde A, Lobo J M. Threshold criteria for conversion of probability of species presence to either-or presence-absence. Acta Oecologica, 2007, 31(3): 361- 369.

[29] 侯寧, 戴強(qiáng), 冉江洪, 焦迎迎, 程勇, 趙成. 大相嶺山系泥巴山大熊貓生境廊道設(shè)計(jì). 應(yīng)用與環(huán)境生物學(xué)報(bào), 2014, 20(6): 1039- 1045.

[30] 張明春. 賀蘭山巖羊的種群動(dòng)態(tài)及警戒行為研究[D]. 上海: 華東師范大學(xué), 2012.

[31] 駱?lè)f. 賀蘭山巖羊(Pseudoisnayaur)和馬鹿(Cervuselaphusalxaicus)的食性及生境選擇比較研究[D]. 哈爾濱: 東北林業(yè)大學(xué), 2011.

[32] 李言闊, 張明海, 蔣志剛. 基于生境可獲得性的完達(dá)山地區(qū)馬鹿(Cervuselaphusxanthopygus)冬季生境選擇. 生態(tài)學(xué)報(bào), 2008, 28(10): 4619- 4628.

[33] 劉振生, 高惠, 滕麗微, 蘇云, 王曉勤, 孔芳毅. 基于MAXENT模型的賀蘭山巖羊生境適宜性評(píng)價(jià). 生態(tài)學(xué)報(bào), 2013, 33(22): 7243- 7249.

[34] 康薩如拉, 牛建明, 張慶, 韓硯君, 董建軍, 張靖. 草原區(qū)礦產(chǎn)開(kāi)發(fā)對(duì)景觀格局和初級(jí)生產(chǎn)力的影響—以黑岱溝露天煤礦為例. 生態(tài)學(xué)報(bào), 2014, 34(11): 2855- 2867.

[35] 張國(guó)鋼, 張正旺. 山西五鹿山地區(qū)鐵礦開(kāi)采對(duì)褐馬雞活動(dòng)的影響. 北京師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 1999, 35(2): 266- 269.

[36] 匡先輝. 基于RS的蜂桶寨大熊貓棲息地景觀格局變化研究[D]. 雅安: 四川農(nóng)業(yè)大學(xué), 2012.

Habitat assessment of red deer (Cervusalashanicus) in the Helan Mountains, China

GAO Hui1, TENG Liwei1,2, WANG Yang1, WANG Jifei3, LIU Zhensheng1,2,*

1CollegeofWildlifeResources,NortheastForestryUniversity,Harbin150040,China2KeyLaboratoryofConservationBiology,StateForestryAdministration,Harbin150040,China3HelanMountainsNationalNatureReserveofNingxia,Yinchuan750021,China

Cervusalashanicusonly occurs in the Helan Mountains. Evaluating the habitat suitability ofCervusalashanicusis the premise and foundation of effective protection for this population. We used a Geographic Information Systems (GIS) and the Maximum Entropy (MaxEnt) model to predict habitat suitability forCervusalashanicusin the Helan Mountains. The data needed in this research, included occurrence records ofCervusalashanicusand environmental factors. We collected 86 GPS coordinates forCervusalashanicusoccurrance obtained by line transect surveys and patrol data from the Helan Mountain National Nature Reserve from 2013 to 2014. Environmental factors were grouped into three categories, including topographic data (elevation, slope, and aspect) extracted from the Digital Elevation Model (DEM), six vegetation factors extracted from the Thematic Mapper (TM) image of the Helan Mountains (2007), and water resources and human interference factors, including roads, mining, and forest protection points extracted from the 1∶50000 vector maps of the Helan Mountains. GIS was used to produce the data needed in the model. The distribution ofCervusalashanicuswas predicted using the model of MaxEnt. The area under the curve (AUC) values of a Receiver Operating Characteristic (ROC) curve was used to assess the accuracy of the MaxEnt model, and the Youden Index was applied to determine the threshold value for habitat classification in the MaxEnt model. The simulated habitat was divided into two categories, including suitable and unsuitable habitat using a threshold value of 0.31. Finally, the distribution map of habitat suitability ofCervusalashanicusin the Helan Mountains was generated from the data and software, and the AUC value for training data and test data were 0.956 and 0.898 respectively, which indicated that the result of the assessment was excellent. The results showed that suitable habitat forCervusalashanicuswas primarily distributed in the central and southern regions of the Ningxia Helan Mountain National Nature Reserve (eastern part) and to a lesser extent in the northern region of the Inner-Mongolia Helan Mountain National Nature Reserve (western part), with a total area of 667.87 km2, accounting for 18.2% of the total assessment area. Most northern areas of Helan Mountains were unsuitable habitat forCervusalashanicusbecause of lower precipitation and more mining. Jackknife tests and the contribution of environmental variables indicated that mining (30.4%), slope (26%), altitude (10.8%), and aspect (8.4%) were the main factors affecting habitat selection ofCervusalashanicus, followed by the forest protection points (5.7%), grasslands (4.2%), roads (4.2%), shrubs (4%), alpine meadows (3.8%), broad-leaved forests (1.3%), and water resources (0.9%). Coniferous forest and desert had little effect on the habitat selection ofCervusalashanicus. The response curves showed that habitat suitability increased gradually with the increase of distance from mine and road, and decreased when the slop was more than 30°. Results from our study suggested that mining was the most serious human interference for theCervusalashanicus. Managers need to take more regulatory actions to protect the suitable habitat areas ofCervusalashanicusin order to promote the conservation of this population.

Cervusalashanicus; habitat suitability assessment; maximum entropy model; Helan Mountains

中央高?；究蒲袠I(yè)務(wù)費(fèi)資助項(xiàng)目(2572014CA03,DL13EA01)；國(guó)家自然科學(xué)基金資助項(xiàng)目(31372221)

2016- 03- 25; 網(wǎng)絡(luò)出版日期:2017- 02- 22

10.5846/stxb201603250534

*通訊作者Corresponding author.E-mail: zhenshengliu@163.net

高惠, 滕麗微, 汪洋, 王繼飛, 劉振生.阿拉善馬鹿(Cervusalashanicus)生境適宜性評(píng)價(jià).生態(tài)學(xué)報(bào),2017,37(11):3926- 3931.

Gao H, Teng L W, Wang Y, Wang J F, Liu Z S.Habitat assessment of red deer (Cervusalashanicus) in the Helan Mountains, China.Acta Ecologica Sinica,2017,37(11):3926- 3931.