熊 姍，張海江，李 成 （1.中國(guó)科學(xué)院大學(xué)，北京 100049；2.中國(guó)科學(xué)院成都生物研究所，四川成都610041；.大自然保護(hù)協(xié)會(huì)，浙江 杭州 11118）

全球大約20%的陸棲脊椎動(dòng)物面臨著多種威脅并經(jīng)歷著種群快速下降和物種滅絕，包括25%的哺乳類、13%的鳥(niǎo)類、22%的爬行類以及41%的兩棲類[1]。對(duì)我國(guó)兩棲類的評(píng)估發(fā)現(xiàn)我國(guó)兩棲類受威脅物種達(dá)到43.1%[2]。顯然，兩棲動(dòng)物的生存境況極其嚴(yán)峻。監(jiān)測(cè)兩棲動(dòng)物種群動(dòng)態(tài)并確定其致危因素是應(yīng)對(duì)危機(jī)的關(guān)鍵步驟[3]。為此，我國(guó)分別于2011年由生態(tài)環(huán)境部南京環(huán)境科學(xué)研究所牽頭啟動(dòng)了“兩棲類示范觀測(cè)項(xiàng)目”[4]；在2014年由中國(guó)科學(xué)院成都生物研究所牽頭組織了“中國(guó)關(guān)鍵地區(qū)兩棲爬行動(dòng)物多樣性監(jiān)測(cè)與研究專項(xiàng)網(wǎng)”[5]。至此，我國(guó)的兩棲動(dòng)物監(jiān)測(cè)有了規(guī)?；捏w系，但目前在監(jiān)測(cè)中還存在如下幾個(gè)問(wèn)題沒(méi)有得到較好的解決：（1）兩棲動(dòng)物分布不集中，在制定監(jiān)測(cè)方案時(shí)需布設(shè)多條樣線或樣方[6]，人工成本較高。（2）物種的種群數(shù)量在季節(jié)以及年度間波動(dòng)較為顯著，氣候因素對(duì)觀測(cè)結(jié)果影響較大[7]，在制定監(jiān)測(cè)方案時(shí)需布置2～4個(gè)季節(jié)的調(diào)查[6]，時(shí)間成本高。（3）種群數(shù)量統(tǒng)計(jì)不準(zhǔn)確，兩棲動(dòng)物具有較強(qiáng)的季節(jié)遷徙習(xí)性，由于不能識(shí)別個(gè)體，只能采用平均個(gè)體發(fā)現(xiàn)率來(lái)代表種群數(shù)量。

被動(dòng)集成應(yīng)答器（passive integrated transponder，PIT）技術(shù)產(chǎn)生后，研究人員在1987年對(duì)于黑足鼬（Mustela nigripes）、蒙眼貂（（Mustela putorius furo）的研究中使用了PIT標(biāo)記，證明了PIT標(biāo)記對(duì)于實(shí)現(xiàn)個(gè)體識(shí)別的有效性[8]，隨后PIT標(biāo)記被廣泛運(yùn)用到野生動(dòng)物的標(biāo)記重捕研究中，包括哺乳類[9]、魚(yú)類[10]、鳥(niǎo)類[11]、爬行類[12]以及兩棲類[13]。FABER[14]首次采用PIT標(biāo)記技術(shù)研究了高山歐螈（Triturus alpestris），雌螈重捕率為74.6%，雄螈重捕率為82.2%，表明PIT標(biāo)記重捕法在兩棲動(dòng)物研究中的適用性，開(kāi)創(chuàng)了將PIT標(biāo)記技術(shù)運(yùn)用到兩棲動(dòng)物研究中的先河。此后PIT標(biāo)記重捕法廣泛運(yùn)用于兩棲動(dòng)物種群動(dòng)態(tài)、遷移模式等研究中。

該文介紹了一種利用繁殖點(diǎn)調(diào)查法與PIT標(biāo)記重捕法相結(jié)合，獲得短期兩棲動(dòng)物重捕數(shù)據(jù)后，基于開(kāi)放式種群標(biāo)記重捕模型的調(diào)查與分析方法。以四川省平武縣老河溝自然保護(hù)區(qū)內(nèi)的中華蟾蜍為代表物種，以5 d野外調(diào)查的標(biāo)記重捕數(shù)據(jù)，準(zhǔn)確快速地估算其種群數(shù)量。該方法提供了一種簡(jiǎn)潔有效的種群數(shù)量調(diào)查方法，對(duì)于珍稀物種的評(píng)估具有較高的價(jià)值。

1 材料和方法

1.1 研究區(qū)域

老河溝縣級(jí)保護(hù)區(qū)地處四川省平武縣東部，東側(cè)毗鄰四川省青川縣唐家河國(guó)家級(jí)自然保護(hù)區(qū)，北側(cè)與甘肅省文縣白水江國(guó)家級(jí)自然保護(hù)區(qū)相鄰，南部為平武縣高村鄉(xiāng)，西邊與木皮鄉(xiāng)、木座鄉(xiāng)接壤，地理坐標(biāo)位為32°36'22″～32°25'52″N，104°32'42″～104°45'25″E。面積約 11 000 hm2。海拔 1 250～3 504 m，相對(duì)高差2 200 m左右，地形較復(fù)雜，平均坡度30°。該地屬于亞熱帶山地氣候，年均溫4～12℃，多年平均降水量950～1 300 mm。

1.2 研究對(duì)象

中華蟾蜍（Bufo gargarizans）隸屬于兩棲綱（Amphibia）、無(wú)尾目（Anura）、蟾蜍科（Bufonidae）、蟾蜍屬（Bufo），主要生活于海拔0～3 500 m的多種生態(tài)環(huán)境中。每年的3—5月，中華蟾蜍的卵產(chǎn)在臨時(shí)性、永久性靜水塘或流溪的洄水塘內(nèi)。卵數(shù)多，卵徑小，發(fā)育迅速，從卵孵化到蝌蚪變態(tài)上陸僅需2～3個(gè)月，季節(jié)性的來(lái)水即可成功繁殖。蝌蚪為靜水型蝌蚪：數(shù)量多，個(gè)體小，尾較短，尾肌弱，越冬的蝌蚪極為少見(jiàn)。比較適合水溫變異幅度大，容易干涸并具有豐富食物的靜水環(huán)境。成體為陸棲生活類型：體肥壯，后肢粗短，蹠突特別發(fā)達(dá)，有游離刃，皮膚厚，富于腺體或具瘰粒，生活在土洞內(nèi)或深埋在水底的泥沙內(nèi)。中華蟾蜍體型較大，雄蟾體長(zhǎng)63～90 mm，雌蟾體長(zhǎng)85～116 mm，是我國(guó)分布最廣的兩棲動(dòng)物物種，也是平原和低山丘陵的常見(jiàn)種[15]；在生活史的各個(gè)階段對(duì)環(huán)境的適應(yīng)能力極強(qiáng)[16]。種群數(shù)量在繁殖期與非繁殖期波動(dòng)較為顯著[17]；具有較強(qiáng)的季節(jié)性遷徙習(xí)性，無(wú)明顯的個(gè)體識(shí)別特征[18]。中華蟾蜍具有了兩棲動(dòng)物調(diào)查3個(gè)難點(diǎn)的代表性和典型性。

1.3 研究方法

1.3.1 繁殖點(diǎn)調(diào)查法

由于兩棲動(dòng)物是變溫動(dòng)物，生理活動(dòng)受到環(huán)境溫度的較大影響。多數(shù)的兩棲動(dòng)物具有明顯的繁殖期和非繁殖期，需要周期性地聚集到特定的繁殖場(chǎng)所進(jìn)行繁殖活動(dòng)。藉此特點(diǎn)，在兩棲動(dòng)物的繁殖季節(jié)，在特定的繁殖場(chǎng)所開(kāi)展集中連續(xù)調(diào)查，即為繁殖點(diǎn)調(diào)查法[19]。在溫帶地區(qū)，兩棲動(dòng)物的繁殖期主要集中在春天和初夏的幾周內(nèi)；在干旱地區(qū)，繁殖期集中在短促的雨季；而少數(shù)熱帶地區(qū)的兩棲動(dòng)物一年四季都有可能繁殖，繁殖期較不明顯。兩棲動(dòng)物在繁殖期的自然聚集，克服了其在非繁殖期分布不集中的難點(diǎn)，因此，繁殖點(diǎn)調(diào)查法的效率很高。該法可以用來(lái)監(jiān)測(cè)物種種群水平的變化、區(qū)系組成在不同時(shí)期的變遷，以及個(gè)體生態(tài)學(xué)研究。

繁殖點(diǎn)調(diào)查法可以調(diào)查成體，也可以調(diào)查幼體。然而，幼體的密度受氣候和捕食者影響很大，在短期內(nèi)也會(huì)出現(xiàn)很大的變化[20]；READING[21]標(biāo)記了7 259只剛變態(tài)的幼體，僅有48只成年后回到出生的池塘繁殖；因此，幼體的密度并不能預(yù)測(cè)成體種群的大小。

在老河溝保護(hù)區(qū)，中華蟾蜍一般在3月開(kāi)始進(jìn)入水塘或河流緩流處繁殖，山區(qū)水塘集雨而成，干涸較快，卵和蝌蚪死亡很多。2013年在保護(hù)區(qū)內(nèi)修建一處人工池塘（32.483 95°N，104.71 967°E，1 373 m），池塘呈梯形，面積為101.25 m2。池塘內(nèi)引入山溪水，水位常年保持在20～30 cm。池內(nèi)雜草叢生，雌雄抱對(duì)后共同移向合適的產(chǎn)卵地點(diǎn)，卵包裹在管狀的膠質(zhì)卵帶內(nèi)，纏繞在水池的水草上。蝌蚪孵化后主要攝食水池中的水草及落葉。由于面積大、常年有水、幼體食物豐富，吸引了大量的中華蟾蜍和中國(guó)林蛙（Rana chensinensis）等兩棲動(dòng)物進(jìn)入繁殖，因此該池塘成為保護(hù)區(qū)內(nèi)重要的兩棲動(dòng)物繁殖點(diǎn)之一。2014年3月12—16日，對(duì)該繁殖點(diǎn)的中華蟾蜍成體種群數(shù)量開(kāi)展了短期調(diào)查。

圖1 中華蟾蜍的繁殖池塘Fig.1 Breeding ponds of Bufo gargarizans

1.3.2 PIT標(biāo)記重捕法

兩棲動(dòng)物的多數(shù)成體在繁殖期間只有較短時(shí)間會(huì)在繁殖場(chǎng)出現(xiàn)[22-23]。要準(zhǔn)確地評(píng)估物種的種群數(shù)量，則需要采用標(biāo)記重捕法進(jìn)行調(diào)查。傳統(tǒng)的兩棲動(dòng)物標(biāo)記多為損傷性的個(gè)體識(shí)別技術(shù)，如剪指（趾）、烙印、植入性熒光標(biāo)記、趾環(huán)、彩漆、色斑等[24-27]。BROWN[28]探究了 PIT 標(biāo)記對(duì)兩棲動(dòng)物的影響，結(jié)果表明，PIT標(biāo)記對(duì)于兩棲動(dòng)物個(gè)體的存活率及身體狀況的負(fù)面影響可忽略不計(jì)。此后，該方法被廣泛應(yīng)用于兩棲動(dòng)物監(jiān)測(cè)和研究，是實(shí)現(xiàn)兩棲動(dòng)物標(biāo)記和個(gè)體識(shí)別的最佳方法[29]。

PIT標(biāo)簽是一種體積小的玻璃封裝的微型芯片，尺寸僅有1.4 mm×8.4 mm，質(zhì)量為0.06 g。每個(gè)標(biāo)簽有唯一的編號(hào)，通過(guò)將標(biāo)簽注射進(jìn)動(dòng)物體內(nèi)，實(shí)現(xiàn)對(duì)動(dòng)物的個(gè)體識(shí)別，通過(guò)多次重捕，快速準(zhǔn)確地評(píng)估野生兩棲動(dòng)物的種群數(shù)量。PIT標(biāo)記重捕法和繁殖點(diǎn)調(diào)查法相結(jié)合[30]，可以大幅提高傳統(tǒng)監(jiān)測(cè)數(shù)據(jù)的準(zhǔn)確率。

2014年3月12日20點(diǎn)開(kāi)始，工作人員捕捉進(jìn)入繁殖場(chǎng)的中華蟾蜍；記錄中華蟾蜍的性別、抱對(duì)情況、測(cè)量頭體長(zhǎng)、體重、將PIT標(biāo)簽用注射器注射到胯部皮膚與肌肉之間，用醫(yī)用膠水封閉注射孔，防止電子標(biāo)簽掉落；同時(shí)，由于PIT標(biāo)簽有丟失的情況[31-32]，為了正確估算標(biāo)簽丟失率，研究中需要采用2種獨(dú)立的標(biāo)記方法進(jìn)行標(biāo)記[33]，所以該研究采用PIT為主要標(biāo)記和剪趾法輔助的綜合標(biāo)記方法，用醫(yī)用酒精消毒后剪取每只成體最長(zhǎng)的第4趾的一節(jié)趾骨作為輔助標(biāo)記，剪趾后使用液體創(chuàng)可貼；全部標(biāo)記后統(tǒng)一釋放回繁殖場(chǎng)。其后的4晚重復(fù)相同的工作；運(yùn)用掃描儀掃描所有個(gè)體的PIT編號(hào)，記錄已標(biāo)記個(gè)體的PIT編號(hào)，對(duì)未標(biāo)記的個(gè)體注射PIT標(biāo)簽進(jìn)行標(biāo)記。

標(biāo)記重捕法的準(zhǔn)確性主要取決于樣本量，當(dāng)重捕幾率小時(shí)樣本量要大于200[34]；該次調(diào)查共用PIT標(biāo)記了212只中華蟾蜍。在5 d的調(diào)查期間，僅有1枚標(biāo)簽脫落（表1）。

表1 中華蟾蜍的標(biāo)記重捕資料Table1 Capture-Mark-RecapturedataofBufogargarizans

1.3.3 分析方法

等捕性檢驗(yàn)：個(gè)體間具等捕性是標(biāo)記重捕模型估算種群參數(shù)的前提條件。所有的標(biāo)記重捕模型都是基于以下假設(shè)：個(gè)體抽樣是完全隨機(jī)的，標(biāo)記個(gè)體與未標(biāo)記個(gè)體有著相同的捕獲概率。如果種群個(gè)體間不具等捕性，則應(yīng)用模型估算的種群參數(shù)均不可靠[35-36]。因此，應(yīng)用標(biāo)記重捕模型分析前首先采用Leslie法對(duì)中華蟾蜍個(gè)體等捕性進(jìn)行檢驗(yàn)[35]。如再捕獲數(shù)與再捕獲期望值差異顯著即表明個(gè)體間不具等捕性。

模型選擇：經(jīng)檢驗(yàn)，中華蟾蜍個(gè)體間滿足等捕性條件后，采用R軟件（R Core Team 2017）中的RMark包[37]，對(duì)中華蟾蜍標(biāo)記重捕資料進(jìn)行分析。RMark包提供了100余個(gè)利用標(biāo)記重捕資料估算種群參數(shù)的模型[38]。

中華蟾蜍在繁殖期不斷地有新的個(gè)體進(jìn)入繁殖場(chǎng)，同時(shí)也有個(gè)體離開(kāi)繁殖場(chǎng)，屬開(kāi)放式種群。因此，選擇適用于開(kāi)放式種群的Jolly-Seber模型[39-40]。Jolly-Seber模型有很多適用于短期或長(zhǎng)期標(biāo)記重捕資料的子模型，POPAN模型是Jolly-Seber模型中的一個(gè)子模型[41-42]，適用于基于連續(xù)的標(biāo)記重捕資料，估算個(gè)體進(jìn)入繁殖場(chǎng)的遷入率（t）、遷出率（1-φ）、存活率（φ）、重捕率（p）、每晚種群數(shù)量（N）、種群總數(shù)（NGross）。MADSEN 等[42]在 1987—1991年采用繁殖點(diǎn)調(diào)查法對(duì)大蟾蜍（Bufo bufo）進(jìn)行了剪趾標(biāo)記重捕，運(yùn)用POPAN模型分別估計(jì)了每年的種群數(shù)量。該文擬采用該模型適于短期標(biāo)記重捕資料的計(jì)算優(yōu)勢(shì)，估算中華蟾蜍繁殖種群數(shù)量。

模型結(jié)果評(píng)估：目前最優(yōu)的方法是通過(guò)比較赤池信息準(zhǔn)則（akaike information criterion，AICc）[43]和權(quán)重（AICc weight）選擇對(duì)數(shù)據(jù)支持度最高的模型；其次，△AICc為該模型與所有候選模型中最小的AICc值之間的差值，△AICc越小，模型越好。當(dāng)其他模型△AICc＜2時(shí)需要對(duì)所有△AICc＜2的模型計(jì)算結(jié)果進(jìn)行模型平均[44]。

2 研究結(jié)果

2.1 捕獲量和標(biāo)記量

2014年3月12—16日，每晚捕獲量為40～64只，共標(biāo)記了212只中華蟾蜍個(gè)體，重捕45次，新標(biāo)記個(gè)體數(shù)隨著總標(biāo)記個(gè)體數(shù)的增加逐漸減少（表1）。平均重捕率22.5%，重捕率最低為12.0%，最高為31.7%（表2）。

表2 各取樣期的重捕率Table 2 Recaptures of the toads PIT tagged in the date

2.2 中華蟾蜍等捕性

利用Leslie法對(duì)中華蟾蜍個(gè)體等捕性進(jìn)行檢驗(yàn)，即以取樣時(shí)點(diǎn)i前個(gè)體的被捕獲次數(shù)（j=0，1，2，…），把在該時(shí)點(diǎn)捕獲釋放的Si只個(gè)體分組，然后利用卡方檢驗(yàn)（χ2-test）比較i+1時(shí)點(diǎn)上述各組的再捕獲數(shù)（b）與其期望值（E）的差異，差異顯著即表明調(diào)查種群個(gè)體不具等捕性，差異不顯著則表明調(diào)查種群個(gè)體具等捕性。Leslie法對(duì)中華蟾蜍個(gè)體等捕性檢驗(yàn)結(jié)果表明，在各取樣期不同捕獲次數(shù)的個(gè)體再捕率均無(wú)顯著差異（χ2=0.79，χ2＜χ0.052=14.07，自由度d?=7，P＞0.05），表明中華蟾蜍個(gè)體間具等捕性，可滿足標(biāo)記重捕模型估算種群參數(shù)的前提條件。

2.3 中華蟾蜍種群數(shù)量

遷入率（t）、存活率（φ）、重捕率（p）是利用POPAN模型計(jì)算種群數(shù)量的3個(gè)重要參數(shù)，將3個(gè)參數(shù)分別定義為常數(shù)值（～1）或隨時(shí)間變動(dòng)的值（～time），共產(chǎn)生8個(gè)隨機(jī)模型，然后用標(biāo)記重捕數(shù)據(jù)分別進(jìn)行擬合。擬合結(jié)果顯示，在8個(gè)隨機(jī)模型中，φ和p為常數(shù)、遷入率隨時(shí)間變化的模型2的AICc值最低，權(quán)重值為0.52，為最優(yōu)模型；其次是模型1，其△AICc值為1.87（表3）。由于2個(gè)模型的△AICc＜2，需對(duì)2個(gè)的模型計(jì)算結(jié)果進(jìn)行平均，則中華蟾蜍每晚在繁殖場(chǎng)的種群數(shù)量為252～377只，繁殖種群總數(shù)為（641±86）只（表4）。

3 討論

3.1 蟾蜍的適應(yīng)能力和瀕危現(xiàn)狀

全球兩棲動(dòng)物共有7 968種，其中蟾蜍科有617種（www.amphibaweb.org），是全世界分布最廣的科級(jí)兩棲動(dòng)物，占據(jù)著多種不同的棲息地。VAN BOCXLAER等[45]對(duì)蟾蜍科的進(jìn)化史進(jìn)行了重構(gòu)，發(fā)現(xiàn)7個(gè)性狀推動(dòng)了蟾蜍的全球擴(kuò)張。7個(gè)性狀包括：（1）適應(yīng)缺水的環(huán)境；（2）有毒腺，防止被捕食；（3）體內(nèi)存儲(chǔ)了大量脂肪，能量足；（4）個(gè)體相對(duì)較大，體表面積比值小，失水少；（5）繁殖場(chǎng)所靈活多樣；（6）卵的數(shù)量大；（7）蝌蚪屬外部營(yíng)養(yǎng)型，自主覓食，不需要母親提供能量。大多數(shù)原始蟾蜍都有這些特性，從而使該群體擁有了擴(kuò)大分布區(qū)的能力。

蟾蜍雖然具有強(qiáng)大的擴(kuò)散和適應(yīng)能力，然而其生存狀況卻不容樂(lè)觀。對(duì)全球5 743種兩棲動(dòng)物的評(píng)估結(jié)果顯示[46]，32.5%（1 856種）的兩棲動(dòng)物受到滅絕威脅；而在評(píng)估的461種蟾蜍科物種中，45.6%的蟾蜍受到滅絕威脅。中國(guó)的蟾蜍科共20種，7種（35.0%）受到滅絕威脅[2]。

中華蟾蜍是我國(guó)分布最廣的兩棲動(dòng)物，數(shù)量大，食性廣，是農(nóng)田和森林里最重要的生態(tài)益蟲(chóng)之一[47]。由于其分布面積大、數(shù)量多、遷移能力強(qiáng)，只能通過(guò)繁殖期的集中調(diào)查評(píng)估其種群數(shù)量；而在非繁殖期，由于中華蟾蜍長(zhǎng)距離的遷徙、分散覓食，調(diào)查難度非常大。通過(guò)對(duì)中華蟾蜍繁殖期的種群數(shù) 量評(píng)估對(duì)于其他兩棲動(dòng)物調(diào)查也有啟示意義。

表3POPAN模型篩選結(jié)果Table 3 Model selection results from the POPAN model

表4 基于POPAN模型的中華蟾蜍種群數(shù)量Table 4 Population size estimates from the POPAN model

3.2 復(fù)合型調(diào)查方法的優(yōu)勢(shì)和不足

監(jiān)測(cè)方法對(duì)監(jiān)測(cè)指標(biāo)和數(shù)據(jù)質(zhì)量具有重要影響。通常兩棲動(dòng)物監(jiān)測(cè)中使用的監(jiān)測(cè)方法有樣線法、樣方法、圍欄陷阱法等[5-6]。以物種多樣性監(jiān)測(cè)為例，對(duì)佛羅里達(dá)州的22種兩棲類進(jìn)行了調(diào)查方法比較[48]，結(jié)果發(fā)現(xiàn)：遇見(jiàn)法發(fā)現(xiàn)12種，樣方法發(fā)現(xiàn)6種，圍欄陷阱法發(fā)現(xiàn)7種。

以種群密度和動(dòng)態(tài)監(jiān)測(cè)為例，在研究卵齒蟾屬（Eleutherodactylus）種群密度和種群動(dòng)態(tài)時(shí)發(fā)現(xiàn)[34]：遇見(jiàn)法與樣線法得到的種群密度較低，而且對(duì)種群動(dòng)態(tài)的預(yù)測(cè)結(jié)果偏差大；標(biāo)記重捕法得到的種群密度也較低，但對(duì)種群動(dòng)態(tài)的預(yù)測(cè)結(jié)果較準(zhǔn)確。遇見(jiàn)法的預(yù)測(cè)準(zhǔn)確率為0.10～0.52，樣線法是0.27，標(biāo)記重捕法的預(yù)測(cè)準(zhǔn)確率為0.35～0.61。研究表明，標(biāo)記重捕法的評(píng)估結(jié)果更精確，偏差更小，而且能夠及時(shí)有效地監(jiān)測(cè)種群數(shù)量和動(dòng)態(tài)[49]。

物種的種群數(shù)量和動(dòng)態(tài)是確定瀕危等級(jí)，制定保護(hù)物種名錄的重要依據(jù)[50]。謝鋒[51]采用繁殖點(diǎn)調(diào)查法和標(biāo)記重捕法調(diào)查了國(guó)家二級(jí)保護(hù)動(dòng)物鎮(zhèn)海棘螈（Echinotriton chinhaiensis）的種群數(shù)量。繁殖點(diǎn)調(diào)查法發(fā)現(xiàn)1997、1998和1999年雌性個(gè)體總數(shù)分別是50、88、89只；對(duì)發(fā)現(xiàn)的部分個(gè)體剪指（趾）標(biāo)記，1997、1998和1999年分別剪指（趾）標(biāo)記了31、53、54只，重捕后經(jīng)Petersen估計(jì)模式推算，1998年的雌性種群數(shù)量為（156±30）只，1999年的雌性種群數(shù)量為（181±11）只。正是基于對(duì)鎮(zhèn)海棘螈種群數(shù)量的準(zhǔn)確評(píng)估結(jié)果，2000年國(guó)家林業(yè)局啟動(dòng)了鎮(zhèn)海棘螈保護(hù)小區(qū)的規(guī)劃和建設(shè)，為我國(guó)的瀕危物種保護(hù)建立了科學(xué)示范。

通過(guò)標(biāo)記來(lái)實(shí)現(xiàn)個(gè)體識(shí)別，通過(guò)多次重捕以提高預(yù)測(cè)準(zhǔn)確率，標(biāo)記重捕法是一種較好的物種多樣性和種群動(dòng)態(tài)監(jiān)測(cè)方法[52]。謝鋒等[53]特別指出，每種監(jiān)測(cè)方法的適用類群不同，監(jiān)測(cè)指標(biāo)各異，在野外工作中同時(shí)采用2～3種方法可以更好地監(jiān)測(cè)物種多樣性和種群的變化。在標(biāo)記重捕研究中，應(yīng)用PIT標(biāo)記有以下優(yōu)點(diǎn)：（1）保留率高：監(jiān)測(cè)期間，僅有一個(gè)個(gè)體標(biāo)簽脫落；（2）識(shí)別精確：標(biāo)簽編號(hào)具有唯一性，不會(huì)出現(xiàn)誤讀的情況，而剪指（趾）法、烙印等標(biāo)記隨著標(biāo)記時(shí)間的增加而錯(cuò)誤識(shí)別的幾率也會(huì)增加；（3）負(fù)擔(dān)?。弘娮訕?biāo)簽重量?jī)H0.06 g；（4）使用壽命長(zhǎng)，適用于中長(zhǎng)期標(biāo)記重捕研究。但PIT標(biāo)記在對(duì)兩棲動(dòng)物標(biāo)記時(shí)有著以下局限：（1）標(biāo)簽會(huì)在動(dòng)物體內(nèi)移動(dòng)，可能會(huì)導(dǎo)致誤認(rèn)為標(biāo)簽丟失；（2）不適用于體型過(guò)小的兩棲動(dòng)物（體重小于2 g）[54]，尤其是兩棲動(dòng)物幼體[29]。

繁殖點(diǎn)調(diào)查法與標(biāo)記重捕法相結(jié)合的復(fù)合型調(diào)查方法有如下優(yōu)點(diǎn)：繁殖點(diǎn)調(diào)查法利用了兩棲動(dòng)物集群繁殖的生物學(xué)習(xí)性，克服了個(gè)體分布不集中的難點(diǎn)，捕捉難度小，人工成本低；在較短的、較固定的繁殖期進(jìn)行調(diào)查，以繁殖種群為調(diào)查對(duì)象，減小了不同季節(jié)的種群數(shù)量的波動(dòng)，調(diào)查時(shí)間也大大縮短；PIT標(biāo)記實(shí)現(xiàn)了個(gè)體識(shí)別，對(duì)動(dòng)物個(gè)體的傷害較小，可以得到較長(zhǎng)期和較準(zhǔn)確的種群動(dòng)態(tài)數(shù)據(jù)，以及生長(zhǎng)、運(yùn)動(dòng)遷移、健康狀況等信息。該方法適用于在臨時(shí)性的池塘、或者永久性的池塘、或溪流中的兩棲動(dòng)物的繁殖種群調(diào)查；尤其對(duì)于繁殖期明顯，短期內(nèi)聚集，爆發(fā)式繁殖的廣布種非常有效；對(duì)于狹域分布的、數(shù)量較少的珍稀物種的快速評(píng)估也具有很高的價(jià)值。

該復(fù)合型調(diào)查方法不能準(zhǔn)確監(jiān)測(cè)的情況如下：（1）繁殖期較不明顯的物種，如熱帶地區(qū)一年內(nèi)多次繁殖的物種；（2）繁殖點(diǎn)較分散的物種，如水樹(shù)蛙屬（Aquixalus）的物種；（3）個(gè)體太小不能攜帶PIT標(biāo)簽的物種，如水蛙屬（Hylarana）、浮蛙屬（Occidozy?ga）、姬蛙屬（Microhyla）、水樹(shù)蛙屬（Aquixalus）、纖樹(shù)蛙屬（Gracixalus）的物種。

3.3 POPAN模型用于種群數(shù)量分析的優(yōu)越性

目前，常用的標(biāo)記重捕模型有Lincoln-Index模型、Jolly-Seber模型、POPAN模型、Robust Design模型等。Lincoln-Index模型是最簡(jiǎn)單的標(biāo)記重捕模型，只需要一次標(biāo)記、一次重捕所獲得的數(shù)據(jù)就可以估計(jì)種群數(shù)量，但只適用于無(wú)出生、死亡、遷入、遷出個(gè)體的封閉式種群，并且所獲得的結(jié)果誤差較大[55]。而野外種群大多是開(kāi)放式種群，Jolly-Seber法被認(rèn)為是研究開(kāi)放式種群最好的方法之一[56]。POPAN模型就是屬于Jolly-Seber法的一個(gè)模型，能夠基于≥4 d的標(biāo)記重捕數(shù)據(jù)估算種群數(shù)量等參數(shù)；而Jolly-Seber模型的其余子模型僅能估算存活率、遷入率等參數(shù)，無(wú)法估算種群大小。Robust Design模型則結(jié)合了開(kāi)放式種群和封閉式種群的特點(diǎn)，但在估算種群數(shù)量等參數(shù)時(shí)需要重捕周期發(fā)生在不同時(shí)間尺度（年、月、日）上，每一次重捕周期內(nèi)有≥2次重捕[57-59]。因此，在重捕周期僅有一個(gè)時(shí)間尺度上的情況下，要獲得種群數(shù)量等參數(shù)時(shí)，POPAN模型是最適用的標(biāo)記重捕模型。

CITTA 等[40]在研究白鯨（Delphinapterus leucas）的種群動(dòng)態(tài)時(shí)，模擬了不同遷移率水平下的標(biāo)記重捕數(shù)據(jù)，然后運(yùn)用這些數(shù)據(jù)跟POPAN模型進(jìn)行了擬合。結(jié)果表明，POPAN模型所估算的種群數(shù)量結(jié)果95%的置信區(qū)間都包含了模擬的種群數(shù)量，表明POPAN模型適用于標(biāo)記重捕數(shù)據(jù)，擬合的結(jié)果較為準(zhǔn)確。但在該研究中，采用的是遺傳標(biāo)記，重捕時(shí)出現(xiàn)了識(shí)別模糊的情況，對(duì)種群數(shù)量估算有一定的影響。而這次對(duì)中華蟾蜍采用了PIT和剪趾的復(fù)合標(biāo)記方法，標(biāo)記的212只個(gè)體在5 d的重捕周期中僅有一個(gè)個(gè)體出現(xiàn)標(biāo)簽丟失的情況；并且每個(gè)PIT標(biāo)簽都有唯一的電子編號(hào)，沒(méi)有識(shí)別錯(cuò)誤；因此，運(yùn)用POPAN模型對(duì)中華蟾蜍種群數(shù)量的估算結(jié)果因PIT標(biāo)記實(shí)現(xiàn)了對(duì)個(gè)體的準(zhǔn)確識(shí)別，對(duì)中華蟾蜍種群數(shù)量的估算更準(zhǔn)確。

致謝：四川省平武縣老河溝保護(hù)區(qū)的陳祥輝、蔡興貴、舒成、張海江、強(qiáng)國(guó)權(quán)、胡小梅、向遵忠、杜惠、吳紅等協(xié)助調(diào)查工作，中國(guó)科學(xué)院成都生物研究所陳有華研究員指導(dǎo)模型選擇和分析，Janak Raj Khatiwada博士校訂英文摘要，特此致謝。