基于生態(tài)位模型的艾比湖國家級(jí)自然保護(hù)區(qū)馬鹿生境評(píng)價(jià)

艾科拜爾·木哈塔爾,熱木圖拉·阿卜杜克熱木,馬合木提·哈力克

新疆大學(xué)生命科學(xué)與技術(shù)學(xué)院,資源生物研究所,烏魯木齊 830046

?

基于生態(tài)位模型的艾比湖國家級(jí)自然保護(hù)區(qū)馬鹿生境評(píng)價(jià)

艾科拜爾·木哈塔爾,熱木圖拉·阿卜杜克熱木,馬合木提·哈力克*

新疆大學(xué)生命科學(xué)與技術(shù)學(xué)院,資源生物研究所,烏魯木齊 830046

生境評(píng)價(jià)和預(yù)測(cè)是對(duì)瀕危物種進(jìn)行有效保護(hù)的基礎(chǔ)。通過2013年9月和2014年10月對(duì)新疆艾比湖國家級(jí)自然保護(hù)區(qū)開展2次秋季野外調(diào)查共收集了92處馬鹿(Cervuselaphus)出現(xiàn)數(shù)據(jù),利用馬鹿出現(xiàn)數(shù)據(jù)作為分布點(diǎn)數(shù)據(jù),選取地形、植被類型和氣候因子3類23種因子作為生境變量,利用MAXENT生態(tài)位模型分析了新疆艾比湖國家級(jí)自然保護(hù)區(qū)馬鹿秋季生境適宜性分布特征和主要生境因子對(duì)馬鹿分布的影響。結(jié)果表明：模型預(yù)測(cè)結(jié)果較高,平均AUC (area under the curve,受試工作者曲線下面值)值為0.976；Jackknife檢驗(yàn)結(jié)果顯示：最熱月最高溫度對(duì)馬鹿生境分布的影響較大。植被類型和坡度對(duì)馬鹿生境分布的影響不大。海拔、年降雨量、氣溫日較差和最熱季平均溫度是影響馬鹿生境分布的主要生境因子。馬鹿秋季生境劃分為高適宜、次適宜、低適宜和不適宜4個(gè)等級(jí),馬鹿的高適宜生境區(qū)主要分布在研究區(qū)域的北部,次適宜及低適宜生境區(qū)則分布于高適宜生境區(qū)的邊緣,而不適宜生境區(qū)主要集中在西部和東部地區(qū)。研究不僅提供了馬鹿在艾比湖的實(shí)際分布狀況,也為馬鹿生境和生境因子的關(guān)系方面提供了一個(gè)重要的科學(xué)依據(jù)。

馬鹿；生態(tài)位模型；生境適宜性；Jackknife檢驗(yàn)

物種保護(hù)不但要維持種群數(shù)量,還要保護(hù)種群的生存環(huán)境,而了解物種種群的分布狀況是保護(hù)其生存環(huán)境的關(guān)鍵前提和基礎(chǔ)[1]。生境為野生動(dòng)物提供了生存、繁衍及種群發(fā)展所必須的資源[2],生境適宜度時(shí)是指生境能支持一個(gè)特定物種的潛在能力[3],是影響野生動(dòng)物生存和發(fā)展的重要因素,對(duì)物種進(jìn)行生境評(píng)價(jià)是對(duì)野生動(dòng)物進(jìn)行有效保護(hù)和管理的第一步[4]。

馬鹿(Cervuselaphus)屬于哺乳綱、偶蹄目、反芻亞目、鹿科、鹿亞科、鹿屬[5],是國家Ⅱ級(jí)保護(hù)動(dòng)物。在中國主要分布于東北林區(qū)、寧夏賀蘭山、新疆北部、甘肅、青海、四川和西藏等地[6]。馬鹿在中國尚有一定數(shù)量,在黑龍江和吉林可能有近10萬只,但由于過量捕獵幼仔和生境的喪失,也逐漸產(chǎn)生危機(jī)。新疆是我國馬鹿的重要分布區(qū),有3個(gè)亞種：阿勒泰馬鹿(C.e.sibiricus)、天山馬鹿(C.e.songaricus)和塔里木馬鹿(C.e.yarkandensis)。塔里木的野生種群已經(jīng)由1.5萬只下降到4000—5000只；阿爾泰馬鹿由20世紀(jì)70年代的10萬只下降到4萬只左右；野生天山馬鹿則正以3000只/a左右的速度銳減[7- 10]。由于長期的亂捕濫獵和生境破壞導(dǎo)致馬鹿種群數(shù)量急劇下降,分布區(qū)也逐漸縮減。而對(duì)保護(hù)物種生境的研究,明確其主要影響因子和分布,是分析物種種群數(shù)量減少、瀕危原因的重要手段,同事也能為制定合理的保護(hù)對(duì)策提供依據(jù)[11]。

近年來MAXENT生態(tài)位模型在國內(nèi)外得到廣泛應(yīng)用,ISI Web of Science數(shù)據(jù)庫檢索結(jié)果表明：MAXENT生態(tài)位模型進(jìn)行預(yù)測(cè)研究的文獻(xiàn)數(shù)和每年的引文數(shù)逐年遞增[12]。有研究顯示,在預(yù)測(cè)物種潛在分布區(qū)時(shí),與同類的物種分布模型(species distribution models,SDMs)相比,MAXENT生態(tài)位模型表現(xiàn)出較好的預(yù)測(cè)能力[12]。在預(yù)測(cè)瀕危物種分布時(shí)也有明顯的優(yōu)勢(shì),如Pena等[13]利用3種生態(tài)位模型評(píng)估了巴西北部瀕危鳥類地理分布,結(jié)果表明MAXENT模型的預(yù)測(cè)精度最高；齊增湘等[14]利用MAXENT模型預(yù)測(cè)秦嶺山系黑熊(Ursusthibetanus)的生境適宜性,結(jié)果顯示模型精度達(dá)到0.991。徐衛(wèi)華等[15]采用MAXENT模型,對(duì)秦嶺地區(qū)的川金絲猴(Rhinopithecusroxllanae)的生境進(jìn)行評(píng)價(jià),結(jié)果顯示模型預(yù)測(cè)結(jié)果為0.979。本文應(yīng)用MAXENT模型對(duì)艾比湖馬鹿生境適宜性進(jìn)行評(píng)價(jià),分析影響馬鹿生存的主要生境因子并模擬其空間分布區(qū)域。

1 研究區(qū)域概況

艾比湖位于天山北側(cè)新疆西部,博爾塔拉蒙古族自治州境內(nèi),在行政區(qū)劃分上地跨博爾塔拉蒙古自治州精河縣、博樂市和阿拉山口口岸區(qū),總面積2670.85 km2。艾比湖為準(zhǔn)噶爾盆地最低洼地,也是新疆僅次于吐魯番艾丁湖的第二低點(diǎn),湖面海拔189 m。年平均氣溫6—8℃,年均降水量90.9 mm,年蒸發(fā)量3790 mm,年平均大風(fēng)天數(shù)達(dá)165 d。區(qū)內(nèi)典型的地帶性土壤為灰漠土、灰棕土和風(fēng)沙土。艾比湖湖底平坦,屬淺水湖；目前平均水深1.4 m,最深3 m。保護(hù)區(qū)內(nèi)目前的各種脊椎動(dòng)物有167種,其中兩棲動(dòng)物1目2科4種(蟾蜍科、蛙科),爬行動(dòng)物2目6科7屬10種,哺乳動(dòng)物6目12科32種,鳥類有16目28科111種,魚類2科10種。野生維管植物131科755屬3552種,種子植物53科191屬[16]。

2 研究方法

2.1 模型介紹

最大熵原理(MAXENT)起源于信息科學(xué),1957年最早由Jaynes提出,在許多學(xué)科中有廣泛應(yīng)用,近年來引起了生態(tài)學(xué)家們的關(guān)注[12]。MAXENT是一種基于生態(tài)位原理的模型,根據(jù)物種“出現(xiàn)點(diǎn)”(出現(xiàn)馬鹿糞便、食跡、足跡、臥跡等位置進(jìn)行定點(diǎn)記錄的數(shù)據(jù))的生境變量特征得出約束條件,探尋此約束條件下最大熵的可能分布,以此來預(yù)測(cè)目標(biāo)物種在研究地區(qū)的生境分布[17]。2004年,Phillips等人開發(fā)了MAXENT軟件,近來廣泛用于物種生境適宜區(qū)的預(yù)測(cè)和評(píng)價(jià),表現(xiàn)出了良好的預(yù)測(cè)能力[18]。該模型提供了自檢驗(yàn)功能,可以自動(dòng)生成ROC曲線進(jìn)行模型的模擬預(yù)測(cè)自檢,且在對(duì)動(dòng)物生境進(jìn)行評(píng)價(jià)與預(yù)測(cè)時(shí),只需動(dòng)物“出現(xiàn)點(diǎn)”的數(shù)據(jù),并且具有較高的精度[19- 21]。

2.2 數(shù)據(jù)來源與處理

2.2.1 馬鹿分布點(diǎn)數(shù)據(jù)

圖1 艾比湖馬鹿出現(xiàn)點(diǎn)分布圖Fig.1 Species occurrence date

馬鹿分布點(diǎn)數(shù)據(jù)來源于2013年9月和2014年10月通過走訪當(dāng)?shù)鼐用窈拖驅(qū)?duì)新疆艾比湖國家級(jí)自然保護(hù)區(qū)開展2次秋季野外巡邏和檢測(cè)以及艾比湖國家級(jí)自然保護(hù)區(qū)管理站的野外檢測(cè),用GPS對(duì)所觀察到的馬鹿實(shí)體或馬鹿出現(xiàn)點(diǎn)位置進(jìn)行定點(diǎn)記錄,在艾比湖與馬鹿同域分布的有蹄類還有鵝喉羚(Gazellasubgutturosa),可根據(jù)體型大小和形態(tài)特征對(duì)它們進(jìn)行區(qū)分,第一次和第二次野外檢測(cè)中分別得到了55和37處馬鹿出現(xiàn)點(diǎn),共得到了92處包含經(jīng)緯度的馬鹿出現(xiàn)點(diǎn)(包括糞便、足跡、食跡、臥跡等)(圖1)。

2.2.2 生境變量數(shù)據(jù)

(1)氣候因子。

從世界氣象數(shù)據(jù)庫(http://www.worldclim.org/)中下載全球19個(gè)生物氣候因子的1 km2分辨率的柵格文件1950年至2000年的平均值,然后用ArcGIS10.1軟件來提取艾比湖區(qū)域的氣候變量作為模型預(yù)測(cè)的生境變量(表1)。

表1 19個(gè)生物氣候變量

(2)地形因子。

地形選取3個(gè)因子,包括海拔(altitude)、坡度(slope)和坡向(aspect)。采用來源于中國科學(xué)院科學(xué)數(shù)據(jù)庫分辨率30m×30m的數(shù)字高程地圖(DEM)數(shù)據(jù),在ArcGIS10.1中進(jìn)行空間分析獲得。

(3)植被數(shù)據(jù)。

先從谷歌剪切器下載艾比湖地區(qū)的地形圖,然后通過目視解析法做出植被分布圖,把艾比湖地區(qū)的植被類型分為泡泡刺(NitrariasphaerocarpaMaxim),草地(Meadow),檉柳(Tamarixspp.),胡楊林(Populuseuphratica),沼澤(Swamp),湖泊(Lake),宜林地(Suitable land for forest),枸杞(Lyciumchinense),蘆葦(Phragmitesaustralis),柳樹(salixbabylonica),梭梭(Haloxylonammodendron),沙漠(Desert),等12類(圖2)。

圖2 研究區(qū)域的植被分布圖Fig.2 Vegetation distribution map of the study area

2.3 MAXENT模型的分析方法

2.3.1 MAXENT模型的數(shù)據(jù)

物種分布點(diǎn)數(shù)據(jù)：將物種分布點(diǎn)的數(shù)據(jù)在Excel中保存為MAXENT軟件可以識(shí)別的以逗號(hào)分隔的CSV數(shù)據(jù)格式,字段包括物種名稱、經(jīng)度和緯度。

生境變量數(shù)據(jù)：以ArcGIS10.1為平臺(tái),分別建立各種生境變量的柵格文件,將所有生境變量數(shù)據(jù),包括氣候因子(BIO1-BIO19)、地形因子(DEM)、植被類型在ArcGIS10.1軟件下轉(zhuǎn)換為ASCII格式,并且按照MAXENT軟件要求將所有圖層統(tǒng)一邊界和坐標(biāo)系統(tǒng)。

2.3.2 影響馬鹿?jié)撛诜植缄P(guān)鍵生境變量的篩選

MAXENT模型預(yù)測(cè)中必須用到的數(shù)據(jù)有兩個(gè),分別是物種空間分布點(diǎn)的數(shù)據(jù)和生境變量數(shù)據(jù),其中物種分布點(diǎn)的數(shù)據(jù)是必不可少的,而生境變量數(shù)據(jù)因?yàn)橛卸鄠€(gè)變量,并不是每個(gè)變量對(duì)物種潛在分布預(yù)測(cè)都是必需的,因此需要剔除掉一些對(duì)MAXENT模型預(yù)測(cè)結(jié)果貢獻(xiàn)較小的生境變量,只留下對(duì)模型預(yù)測(cè)結(jié)果貢獻(xiàn)較大的關(guān)鍵限制因子,以提高模型的預(yù)測(cè)精度。在MAXENT模型中軟件本身提供了實(shí)現(xiàn)這一步驟的模塊,即刀切法(Jackknife檢驗(yàn))測(cè)定變量的重要性,通過比較每個(gè)生境變量對(duì)模型結(jié)果的影響程度來來剔除影響較小的變量。

2.3.3 馬鹿生境模型構(gòu)建

將物種的分布點(diǎn)的數(shù)據(jù)和生境變量數(shù)據(jù)導(dǎo)入MAXEN模型中,隨機(jī)選取75%的分布點(diǎn)的數(shù)據(jù)用于建立模型,剩下25%的分布點(diǎn)的數(shù)據(jù)用于模型檢驗(yàn),其余參數(shù)設(shè)為模型默認(rèn)值,構(gòu)建馬鹿分布模型進(jìn)行模擬,取10次模擬結(jié)果的平均值作為最終模擬結(jié)果[22]。

2.3.4 MAXENT模型的檢驗(yàn)方法

圖3 MAXENT模型的ROC曲線及AUC面積Fig.3 ROC Curve and AUC values of the MAXENT model

模型的檢驗(yàn)采用受試者工作特征曲線,即ROC曲線(receiver operating characteristic curve),該方法最早用于雷達(dá)信號(hào)接受能力的評(píng)估[23],后被廣泛用于評(píng)價(jià)醫(yī)學(xué)診斷試驗(yàn)試驗(yàn)性能,ROC曲線以假陰性率為橫坐標(biāo),真陽性率為縱坐標(biāo),繪制而成,ROC曲線下的面積即為AUC值(area under curve)[24]。AUC值越大,表示生境變量與預(yù)測(cè)物種地理分布模型之間相關(guān)性越大,越能將該物種有分布和無分布判別開,預(yù)測(cè)效果也就越好[25- 26]。評(píng)價(jià)標(biāo)準(zhǔn)為：AUC值為0.50—0.60,失敗(fail)；0.60—0.70,較差(poor)；0.70—0.80,一般(fair)；0.80—0.90,好(good)；0.90—1.0,非常好(excellent)[27- 28](圖3)。

3 結(jié)果與分析

3.1 MAXENT預(yù)測(cè)結(jié)果檢測(cè)

ROC曲線評(píng)價(jià)結(jié)果顯示：AUC值為0.976,表明MAXENT模型預(yù)測(cè)結(jié)果非常好。所建模型可用于艾比湖地區(qū)馬鹿生境評(píng)價(jià)。

3.2 生境變量篩選結(jié)果

由圖4刀切法測(cè)定生境變量的重要性中可以看出,植被類型、坡度對(duì)其分布影響最小,因此,在選擇最優(yōu)秀MAXENT預(yù)測(cè)模型時(shí)剔除這兩個(gè)影響最小的因子來構(gòu)建最優(yōu)秀模型進(jìn)行預(yù)測(cè)。

從圖中可以看出最熱月最高溫度(BIO5)對(duì)艾比湖馬鹿分布影響最大,其次是海拔、年降雨量(BIO12)、氣溫日較差(BIO2)、最熱季平均溫度(BIO1O)、最冷月最低溫度(BIO6)和最干季降水量(BIO17)。

圖4 MAXENT 模型刀切法測(cè)定的各種生境變量的重要性 Fig.4 MAXENT Jackknife measuring of habitat variable importance

3.3 各種生境變量對(duì)MAXENT模型預(yù)測(cè)的影響

圖5列出了全部生境變量的反饋曲線,顯示了每個(gè)生境變量對(duì)模型預(yù)測(cè)的影響。從圖中可以看到,植被類型和坡度的變化對(duì)MAXENT模型預(yù)測(cè)的結(jié)果不明顯,海拔、最熱月最高溫度(BIO5)、年降雨量(BIO12)、溫日較差(BIO2)的變化對(duì)MAXENT模型預(yù)測(cè)的結(jié)果很明顯,而其他變量的變化對(duì)模型預(yù)測(cè)的結(jié)果沒有較大的影響。

圖5 MAXENT 模型中研究物種對(duì)各種生境變量的響應(yīng)曲線反饋曲線Fig.5 Response curves of habitat variables in MAXENT models

3.4 艾比湖馬鹿生境適宜性分布

利用MAXENT模型對(duì)艾比湖馬鹿的生境適宜性進(jìn)行評(píng)價(jià),將模型輸出的ASCII文件導(dǎo)入到ArcGIS10.1中轉(zhuǎn)換成浮點(diǎn)型柵格數(shù)據(jù),按照專家經(jīng)驗(yàn)法將生境適宜圖重新分類4個(gè)適宜等級(jí):0—0.05為不適宜區(qū),0.05—0.2為低適宜區(qū),0.2—0.5為次適宜區(qū),0.5—1為高適宜區(qū),最終得到艾比湖地區(qū)馬鹿的生境適宜性分布圖(圖6)

圖6 艾比湖馬鹿生境分布Fig.6 Habitats distribution of Cervuselaphus in the Ebinur Reserve

從圖中可以看出,艾比湖馬鹿的高適宜生境區(qū)主要分布在研究區(qū)域的北部,次適宜及低適宜生境區(qū)則分布于高適宜生境區(qū)的邊緣,而不適宜生境區(qū)主要集中在西部和東部地區(qū),其中適宜生境區(qū)和不適宜生境區(qū)占總面積的9.4%和90.6%。

4 討論

生境評(píng)價(jià)模型目前可以分為3種類型：機(jī)理模型、回歸模型和生態(tài)位模型。機(jī)理模型并不需要物種分布點(diǎn)數(shù)據(jù),而是根據(jù)生境因子對(duì)物種分布的影響建立相應(yīng)的評(píng)價(jià)準(zhǔn)則,進(jìn)而模擬得到物種的適宜性生境,但因其沒有考慮生境的可達(dá)性,且在對(duì)因子劃分等級(jí)及確定權(quán)重上存在主觀性,因此具有一定局限性[4]。相比回歸模型,生態(tài)位模型只需動(dòng)物“出現(xiàn)點(diǎn)”的數(shù)據(jù)而不需要“未出現(xiàn)點(diǎn)”的數(shù)據(jù),根據(jù)在艾比湖野外考察的實(shí)際經(jīng)驗(yàn),發(fā)現(xiàn)在某個(gè)地點(diǎn)未發(fā)現(xiàn)馬鹿的蹤跡并不代表馬鹿未曾在此地出現(xiàn)過,因此記錄的“出現(xiàn)點(diǎn)”的數(shù)據(jù)并非十分準(zhǔn)確,此外,物種“出現(xiàn)點(diǎn)”數(shù)據(jù)的多少對(duì)模型預(yù)測(cè)結(jié)果也有很大影響,通常物種分布點(diǎn)的數(shù)據(jù)越多,模型預(yù)測(cè)結(jié)果越可靠,但該模型在“出現(xiàn)點(diǎn)”數(shù)據(jù)較少的情況下,也具有較高的預(yù)測(cè)精度,如Silva等[29]在研究蜥蜴的1個(gè)稀有種的潛在分布區(qū)時(shí),僅僅利用17個(gè)分布點(diǎn)的信息構(gòu)建模型,得到的模型精度達(dá)到0.99,在獲得少量分布點(diǎn)的信息情況下模型精度仍達(dá)到優(yōu)秀水平。該方法避免了機(jī)理模型的主觀性和回歸模型的輸入數(shù)據(jù)不易獲取等缺點(diǎn),更適用于分布點(diǎn)的信息匱乏的野生動(dòng)物生境分布研究[30]。

本研究基于MAXENT生態(tài)為模型,結(jié)合ArcGIS10.1,以地形因子、植被類型和氣候因子作為生境變量因子,與分布點(diǎn)的數(shù)據(jù)結(jié)合對(duì)新疆艾比湖國家級(jí)自然保護(hù)區(qū)里面的馬鹿的生境進(jìn)行評(píng)價(jià),模型分析表明：最熱月最高溫度(BIO5)、年降雨量(BIO12)、氣溫日較差(BIO2)、最熱季平均溫度(BIO1O)、最冷月最低溫度(BIO6)和最干季降水量(BIO17)等氣候因子對(duì)艾比湖馬鹿的分布影響起重要的作用,在地形因子方面,海拔最為重要,其次是坡向。李言闊等[31]的研究結(jié)果也指出海拔和坡向是黑龍江省完達(dá)山地區(qū)馬鹿生境分布的主要影響因子。本研究中坡度因子對(duì)馬鹿生境分布的影響很小,主要是因?yàn)轳R鹿行動(dòng)敏捷、活動(dòng)范圍較大,對(duì)坡度要求較低,同時(shí),本研究區(qū)地勢(shì)比較平緩,幾乎沒有坡度很大的地方。植被類型對(duì)艾比湖馬鹿的生境分布沒有太大的影響,主要原因是雖然把植被分為不同的幾個(gè)類型,但是艾比湖國家級(jí)自然保護(hù)區(qū)里面大部分地區(qū)的植被類型是混交林,同一區(qū)域范圍會(huì)有各種植被共存,這是導(dǎo)致植被類型對(duì)艾比湖馬鹿生境分布沒有明顯的影響。

研究過程中沒有加上水源因素,是因?yàn)榘群壹?jí)自然保護(hù)區(qū)大部分區(qū)域總體上類似,以前的水源都不存在了,即艾比湖里面不同的區(qū)域離水源距離沒有明顯的差異,而且吳文等[32]利用MAXENT生態(tài)為模型對(duì)小興安嶺南麓馬鹿冬季適宜生境進(jìn)行研究,結(jié)果發(fā)現(xiàn),水源對(duì)馬鹿生境分布的影響很低,所以以艾比湖的實(shí)際情況出發(fā)本研究沒有考慮了水源因素。

致謝：新疆博州艾比湖國家級(jí)自然保護(hù)區(qū)員工幫助采樣,新疆大學(xué)資環(huán)與環(huán)境科學(xué)學(xué)院許仲林老師幫助軟件使用、數(shù)據(jù)分析及討論,新疆大學(xué)資環(huán)與環(huán)境科學(xué)學(xué)院阿不都拉·阿不力孜博士幫助寫作,特此感謝。

[1] Wilson J W, Sexton J O, Jobe R T, Haddad N M. The relative contribution of terrain, land cover, and vegetation structure indices to species distribution models. Biological Conservation, 2013, 164: 170- 176.

[2] 王金亮, 陳姚. 3S技術(shù)在野生動(dòng)物生境研究中的應(yīng)用. 地理與地理信息科學(xué), 2004, 20(6): 44- 47.

[3] 王秀磊. 普氏原羚生境的景觀動(dòng)態(tài)與適宜性評(píng)價(jià)研究[D]. 北京: 中國林業(yè)科學(xué)研究院, 2004.

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

[5] 吳建平. 哺乳動(dòng)物分類講義. 哈爾濱: 東北林業(yè)大學(xué), 2005:120- 133.

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

[7] 許慶翔, 張明海, 路秉信. 黑龍江省野生馬鹿種群資源調(diào)查研究. 經(jīng)濟(jì)動(dòng)物學(xué)報(bào), 2004, 4(1): 57- 62.

[8] 高行宜,姚軍. 北天山喀拉烏成山的馬鹿資源調(diào)查. 地方病通報(bào), 2007, 22(2): 6- 8.

[9] 艾尼瓦爾·吐米爾, 馬合木提·哈力克. 塔里木河中游塔里木馬鹿種群數(shù)量及其棲息地現(xiàn)狀的初步分析. 新疆農(nóng)業(yè)科學(xué), 2008, 45(4): 743- 749.

[10] 張顯理, 李志剛, 李正, 馬勇璽, 張鐵師, 翟浩. 寧夏賀蘭山馬鹿春季種群數(shù)量與種群動(dòng)態(tài)研究. 寧夏大學(xué)學(xué)報(bào): 自然科學(xué)版, 2006, 27(3): 263- 265.

[11] 歐陽志云, 劉建國, 肖寒, 譚迎春, 張和民.臥龍自然保護(hù)區(qū)大熊貓生境評(píng)價(jià). 生態(tài)學(xué)報(bào), 2001, 21(11): 1869- 1874.

[12] 邢丁亮, 郝占慶. 最大熵原理及其在生態(tài)學(xué)研究中的應(yīng)用. 生物多樣性, 2011, 19(3): 295- 302.

[13] 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.

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

[15] 徐衛(wèi)華, 羅翀. MAXENT模型在秦嶺川金絲猴生境評(píng)價(jià)中的應(yīng)用. 森林工程, 2010, 26(2): 1- 3, 26- 26.

[16] 陳蜀江, 侯平, 李文華, 李虎. 新疆艾比湖濕地自然保護(hù)區(qū)綜合科學(xué)考察. 烏魯木齊: 新疆科學(xué)技術(shù)出版社, 2006.

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

[18] 羅翀, 徐衛(wèi)華, 周志翔, 歐陽志云, 張路. 基于生態(tài)位模型的秦嶺山系林麝生境預(yù)測(cè). 生態(tài)學(xué)報(bào), 2011, 31(5): 1221- 1229.

[19] Stachura-SkierczyńskaK, Tumiel T, SkierczyńskiM. Habitat prediction model for three-toed woodpecker and its implications for the conservation of biologically valuable forests. Forest Ecology and Management, 2009, 258(5): 697- 703.

[20] Haegeman B, Etienne R S. Entropy maximization and the spatial distribution of species. The American Naturalist, 2010, 175(4): E74-E90.

[21] 李明陽, 巨云為, Kumar S, Stohlgren T J. 美國大陸外來入侵物種斑馬紋貽貝(Dreissenapolymorpha)潛在生境預(yù)測(cè)模型. 生態(tài)學(xué)報(bào), 2008, 28(9): 4253- 4258.

[22] Behdarvand N, Kaboli M, Ahmadi M, Nourani E, Mahini A S, Aghbolaghi M A. Spatial risk model and mitigation implications for wolf-human conflict in a highly modified agroecosystem in western Iran. Biological Conservation, 2014, 177: 156- 164.

[23] Leshowitz B. Comparison of ROCCurvesfrom one-and two-interval rating-scale procedures. The Journal of the Acoustical Society of America, 1969, 46(2B): 399- 402.

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

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

[26] Catry F X, Rego F C, Ba??oF L, Moreira F. Modeling and mapping wildfire ignition risk in Portugal. International Journal of Wildland Fire, 2009, 18(8): 921- 931.

[27] Swets J A. Measuring the accuracy of diagnostic systems. Science, 1998, 240(4857): 1285- 1293.

[28] 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.

[29] De Novaes e Silva V, Pressey R L, Machado R B, VanDerWal J, Wiederhecker H C,Werneck F P, Colli G R. Formulating conservation targets for a gap analysis of endemic lizards in a biodiversity hotspot. Biological Conservation, 2014, 180: 1- 10.

[30] 陳新美, 雷淵才, 張雄清, 賈宏炎. 樣本量對(duì)MaxEnt模型預(yù)測(cè)物種分布精度和穩(wěn)定性的影響. 林業(yè)科學(xué), 2012, 48(1): 53- 59.

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

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

Assessing habitat suitability forCervuselaphusin the Ebinur Lake National Nature Reserve

Akbar Muhtar,Rahmutulla Abdukerim,Mahmut Halik*

CollegeofLifeScienceandTechnology,XinjiangUniversity,Urumqi830046,China

Habitat assessing and prediction is important when attempting to protect endangered species. In the Ebinur Lake National Nature Reserve, Xinjiang, A total of 92 recorded GPS coordinates showed thatCervuselaphuswas present in the reserve in September 2013 and October 2014. We grouped 23 habitat predictor variables into three classes—terrain, vegetation structure, and climatic factors. We used the MAXENT model to predict the potentially suitable autumn habitat distribution forC.elaphusand determined the contribution each habitat factor made to the distribution characteristics. The accuracy of our prediction models was accessed by the area under the curve (AUC) values for a receiver operating characteristic (ROC) curve, and the validation showed that the results had high average AUC of 0.976. The results of a Jackknife test indicated that the maximum temperature in the warmest month had the strongest influence on autumnC.elaphushabitat suitability, followed by altitude, annual precipitation, mean diurnal range, and mean temperature during the warmest quarter. Vegetation type and slope had little effect on habitat distribution. The simulated habitat was divided into four classes—the most suitable habitat, moderate suitability habitat, low suitability habitat, and unsuitable habitat. The most suitable habitats forC.elaphuswere mainly distributed in the northern parts of the study area, and habitats with moderate and low suitability were in the marginal areas of the most suitable habitats. In contrast, the western and eastern parts of the study area were classified as unsuitable habitats forC.elaphus. This study provides information on the actual distribution ofC.elaphusin the Ebinur Lake National Nature Reserve and identifies important indicators that could be investigated by future studieson this species and its habitat.

Cervuselaphus; niche model;habitat suitability; Jackknife test

國家自然科學(xué)基金項(xiàng)目(31360266,31560600)；新疆大學(xué)自治區(qū)重點(diǎn)學(xué)科資助項(xiàng)目(動(dòng)物學(xué)科)；中華人民共和國全球環(huán)境基金艾比湖流域可持續(xù)管理與生物多樣性保護(hù)項(xiàng)目(zx- 3- 1)

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

10.5846/stxb201603180488

*通訊作者Corresponding author.E-mail: mahmuthalik@xju.edu.cn

艾科拜爾·木哈塔爾,熱木圖拉·阿卜杜克熱木,馬合木提·哈力克.基于生態(tài)位模型的艾比湖國家級(jí)自然保護(hù)區(qū)馬鹿生境評(píng)價(jià).生態(tài)學(xué)報(bào),2017,37(11):3919- 3925.

Akbar Muhtar,Rahmutulla Abdukerim,Mahmut Halik.Assessing habitat suitability forCervuselaphusin the Ebinur Lake National Nature Reserve.Acta Ecologica Sinica,2017,37(11):3919- 3925.