劉 飛，梁智策，潘 帥，孫安俯，蔡杏偉，郭志強(qiáng)

南渡江下游海南長(zhǎng)臀鮠生物學(xué)參數(shù)估算

劉 飛1，梁智策1，潘 帥1，孫安俯1，蔡杏偉2，郭志強(qiáng)1

（1. 海南大學(xué)海洋學(xué)院//南海海洋資源利用國(guó)家重點(diǎn)實(shí)驗(yàn)室，?？?570228；2. 海南省海洋與漁業(yè)科學(xué)院，?？?571126）

【】開(kāi)展海南長(zhǎng)臀鮠（）的生物學(xué)及種群資源評(píng)估，為南渡江下游的海南長(zhǎng)臀鮠的種群保護(hù)和資源恢復(fù)提供基礎(chǔ)依據(jù)。本研究從2019年10月至2020年9月在南渡江下游進(jìn)行逐月采樣，共采集677尾海南長(zhǎng)臀鮠樣品。利用FiSAT Ⅱ軟件中的ELEFAN I 法對(duì)南渡江下游海南長(zhǎng)臀鮠的種群參數(shù)進(jìn)行估算，利用Beverton-Holt動(dòng)態(tài)模型對(duì)南渡江下游海南長(zhǎng)臀鮠的資源動(dòng)態(tài)進(jìn)行評(píng)估。南渡江下游海南長(zhǎng)臀鮠體長(zhǎng)，體質(zhì)量關(guān)系為：= 0.000 0073.16(2= 0.97；< 0.05)；海南長(zhǎng)臀鮠von Bertalanffy 生長(zhǎng)方程各參數(shù)為漸近體長(zhǎng)∞= 299.25 mm、生長(zhǎng)參數(shù)= 0.41，起點(diǎn)年齡0= -0.09 齡；根據(jù)生長(zhǎng)方程推算出生長(zhǎng)拐點(diǎn)年齡為3齡，拐點(diǎn)體長(zhǎng)約為226.38 mm，拐點(diǎn)體質(zhì)量為164.45 g。利用Pauly經(jīng)驗(yàn)公式計(jì)算出自然死亡系數(shù)（）為0.51，體長(zhǎng)變換漁獲曲線(xiàn)法分別計(jì)算出總死亡系數(shù)（）為1.23，捕撈死亡系數(shù)（）為0.72，捕撈開(kāi)發(fā)率（）為0.58；動(dòng)態(tài)綜合模型估算南渡江下游海南長(zhǎng)臀鮠的開(kāi)捕年齡為1.47 齡，對(duì)應(yīng)開(kāi)捕體長(zhǎng)為149.5 mm。海南長(zhǎng)臀鮠的漁業(yè)資源處于過(guò)度開(kāi)發(fā)狀態(tài)。

海南長(zhǎng)臀鮠；生物學(xué)參數(shù)；過(guò)度捕撈；南渡江

海南長(zhǎng)臀鮠（）隸屬于鲇形目（Siluriformes），長(zhǎng)臀鮠科（Cranoglanididae），長(zhǎng)臀鮠屬（），為海南島特有淡水魚(yú)類(lèi)，主要分布于海南島昌化江及南渡江水系[1]。近10 a以來(lái)，海南島淡水水域生態(tài)環(huán)境受到的影響不斷加劇，淡水水生生物的棲息環(huán)境受到破壞[2,3]。此外，尼羅羅非魚(yú)（）、云斑尖塘鱧（）等外來(lái)物種的入侵，也對(duì)海南長(zhǎng)臀鮠生存造成了巨大威脅[2]。目前，根據(jù)《中國(guó)脊椎動(dòng)物紅色名錄》，海南長(zhǎng)臀鮠的瀕危等級(jí)為瀕危（EN）[4]。有關(guān)海南長(zhǎng)臀鮠的研究大多集中在分類(lèi)及遺傳方面[5-7]，有關(guān)海南長(zhǎng)臀鮠生物學(xué)及資源動(dòng)態(tài)的研究未見(jiàn)報(bào)道。魚(yú)類(lèi)的生長(zhǎng)與死亡參數(shù)不僅是魚(yú)類(lèi)生物學(xué)研究的基礎(chǔ)，同時(shí)還是研究魚(yú)類(lèi)種群結(jié)構(gòu)變動(dòng)的重要依據(jù)[8]。魚(yú)類(lèi)的硬組織（鱗片、耳石等）被廣泛應(yīng)用于魚(yú)類(lèi)年齡與生長(zhǎng)的研究[9]，但熱帶地區(qū)的年輪存在不清晰的缺點(diǎn)，年齡鑒定較為困難[10]。而體長(zhǎng)頻率分析法在估算魚(yú)類(lèi)資源量方面，具有操作方便及無(wú)需年齡數(shù)據(jù)而被廣泛應(yīng)用于生長(zhǎng)與死亡參數(shù)的估算[11-13]。此外，魚(yú)類(lèi)準(zhǔn)確資源量評(píng)估是開(kāi)展?jié)O業(yè)資源保護(hù)的重要依據(jù)[14]。Beverton和Holt建立的單位補(bǔ)充量漁獲量的模型（B-H模型）被廣泛用來(lái)分析過(guò)度捕撈對(duì)漁業(yè)資源造成的許多影響[15,16]。

海南島地處熱帶，氣候高溫多雨，豐富水資源和復(fù)雜的生態(tài)環(huán)境使得海南島具有豐富內(nèi)陸?hù)~(yú)類(lèi)資源[1,17-18]。本研究以南渡江下游的海南長(zhǎng)臀鮠為研究對(duì)象，通過(guò)計(jì)算該物種體長(zhǎng)組成，體長(zhǎng)-體質(zhì)量關(guān)系、以及利用體長(zhǎng)頻率方法估算海南長(zhǎng)臀鮠的生長(zhǎng)及死亡參數(shù)，利用B-H模型對(duì)其種群資源進(jìn)行評(píng)估，為海南長(zhǎng)臀鮠的種群資源保護(hù)提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)域與樣品采集

本研究樣品采集于南渡江下游（圖 1），采樣時(shí)間為2019年10月至2020年9月，每月采樣時(shí)間定為10號(hào)至15號(hào)。每月利用3條多網(wǎng)目復(fù)合刺網(wǎng)（長(zhǎng)：30 m; 高：1.5 m；網(wǎng)目（2/cm）分別為1.0、6.0和12.5）和地籠（長(zhǎng)：4 m；網(wǎng)目：2= 0.4 cm）。將樣品收集起來(lái)，測(cè)量南渡江下游海南長(zhǎng)臀鮠的體長(zhǎng)（精確到1 mm）和體質(zhì)量（精確到0.01 g）等生物學(xué)參數(shù)。每月捕撈海南長(zhǎng)臀鮠數(shù)量見(jiàn)表1。

圖1 南渡江下游海南長(zhǎng)臀鮠采樣位點(diǎn)

1.2 研究方法

1.2.1 體長(zhǎng)與體質(zhì)量關(guān)系 利用=aL對(duì)樣品體長(zhǎng)與體質(zhì)量關(guān)系進(jìn)行擬合，關(guān)系式中為樣品體質(zhì)量（g），為肥滿(mǎn)度，為關(guān)系式系數(shù)，為魚(yú)體體長(zhǎng)（mm）。

1.2.2 生物學(xué)參數(shù)的計(jì)算 對(duì)樣品體長(zhǎng)以10 mm為組距進(jìn)行分組[19]，共50 ～ 59 mm，60 ～ 69 mm，270 ～ 279 mm，280 ～ 289 mm 等25組，并計(jì)算每組頻率。輸入海南長(zhǎng)臀鮠體長(zhǎng)頻率數(shù)據(jù)，利用FiSAT Ⅱ軟件中的 ELEFAN Ⅰ版塊估算海南長(zhǎng)臀鮠生長(zhǎng)系數(shù)（）和漸近體長(zhǎng)（∞）[20]。然后根據(jù)Pauly經(jīng)驗(yàn)公式，輸入和∞計(jì)算其理論生長(zhǎng)起點(diǎn)年齡0[21]：

ln (-0) = -0.3922 – 0.2752 ln∞- 1.038 ln，（1）

海南長(zhǎng)臀鮠總死亡系數(shù)（）利用體長(zhǎng)變換漁獲曲線(xiàn)法計(jì)算，自然死亡系數(shù)（）利用自然死亡系數(shù)估算模塊估算[19,22]。

表1 南渡江下游海南長(zhǎng)臀鮠樣本數(shù)量和規(guī)格

1.2.3 生長(zhǎng)特性指標(biāo) 基于von Bertalanffy 生長(zhǎng)方程[23]，輸入0（理論生長(zhǎng)起點(diǎn)年齡）、L（體長(zhǎng)）、m（體質(zhì)量）等參數(shù)，擬合海南長(zhǎng)臀鮠的體長(zhǎng)、體質(zhì)量生長(zhǎng)方程：

對(duì)海南長(zhǎng)臀鮠的體長(zhǎng)，體質(zhì)量生長(zhǎng)方程的進(jìn)行一階和二階求導(dǎo)，計(jì)算其體長(zhǎng)、體質(zhì)量的生長(zhǎng)速度和加速度曲線(xiàn)方程。

1.2.4 魚(yú)類(lèi)單位補(bǔ)充量漁獲量

1.2.4.1 估算平均選擇體長(zhǎng) 利用體長(zhǎng)變換曲線(xiàn)對(duì)總死亡系數(shù)（）進(jìn)行計(jì)算：ln (/△) =+t，其中為第個(gè)體長(zhǎng)組的數(shù)量（尾），△為其他體長(zhǎng)組與第個(gè)體長(zhǎng)組的年齡間隔，為常數(shù)，t為第個(gè)體長(zhǎng)組的相對(duì)年齡，-即為總死亡系數(shù)（）[24]。然后計(jì)算各點(diǎn)觀(guān)察與期望的比值及其累計(jì)率，當(dāng)累計(jì)率達(dá)50.0%時(shí)所對(duì)應(yīng)的體長(zhǎng)即為海南長(zhǎng)臀鮠的開(kāi)捕體長(zhǎng)（）[25]。

1.2.4.2 相對(duì)單位補(bǔ)充量漁獲量分析 Beverton和Holt提出動(dòng)態(tài)綜合模型（dynamic pooled model）以計(jì)算魚(yú)類(lèi)種群參數(shù)變化和漁獲量之間的關(guān)系[15]。Pauly對(duì)模型進(jìn)行了改進(jìn)以分析相對(duì)/（單位補(bǔ)充量漁獲量）和相對(duì)單位補(bǔ)充量生物量/（相對(duì)單位補(bǔ)充量生物量）[26]。在FiSAT Ⅱ軟件中，開(kāi)發(fā)率一般有三種：1）max，是最大漁獲量時(shí)的開(kāi)發(fā)率[27]；2）0.5，為資源量下降到原始水平50.0%時(shí)的開(kāi)發(fā)率[27]；3）0.1，為/邊際增長(zhǎng)減少10.0%時(shí)的開(kāi)發(fā)率[27]。

2 結(jié)果

2.1 體長(zhǎng)分布

本研究在南渡江下游共采集海南長(zhǎng)臀鮠677 尾，體長(zhǎng)分布為56 ～ 285 mm，平均體長(zhǎng)為 194.73 mm，南渡江下游海南長(zhǎng)臀鮠的優(yōu)勢(shì)體長(zhǎng)為200 ～ 290 mm，占全部體長(zhǎng)組的84.47%（圖 2）。

圖2 海南長(zhǎng)臀鮠體長(zhǎng)分布

2.2 體長(zhǎng)與體質(zhì)量關(guān)系

通過(guò)對(duì)南渡江下游海南長(zhǎng)臀鮠的體長(zhǎng)和體質(zhì)量進(jìn)行擬合，海南長(zhǎng)臀鮠體長(zhǎng)和體質(zhì)量擬合的關(guān)系式為：

= 0.000 0073.16(2= 0.97；< 0.05;= 677)，（4）

其體長(zhǎng)體質(zhì)量關(guān)系的冪指數(shù)為3.16，與3具有顯著性差異（< 0.05），因此海南長(zhǎng)臀鮠為正異速生長(zhǎng)（圖 3）。

圖3 海南長(zhǎng)臀鮠體長(zhǎng)與體質(zhì)量關(guān)系式

2.3 生長(zhǎng)方程及生長(zhǎng)參數(shù)

利用FISAT Ⅱ軟件中的ELEFAN Ⅰ方法中掃描功能計(jì)算海南長(zhǎng)臀鮠的漸近體長(zhǎng)（∞）為299.25 mm，生長(zhǎng)參數(shù)為0.41。然后根據(jù)Pauly提出的用于計(jì)算理論生長(zhǎng)起點(diǎn)年齡的經(jīng)驗(yàn)公式計(jì)算出0為-0.09，因此海南長(zhǎng)臀鮠的齡時(shí)的體長(zhǎng)與體質(zhì)量理論生長(zhǎng)方程分別為：

L= 299.25 [1 - e-0.41 (t + 0.09)]， （5）

m= 467.05 [1 - e-0.41 (t + 0.09)]3.16， （6）

海南長(zhǎng)臀鮠年齡與體長(zhǎng)、體質(zhì)量擬合曲線(xiàn)見(jiàn)圖4。

圖4 海南長(zhǎng)臀鮠von Bertalantly生長(zhǎng)曲線(xiàn)

利用體長(zhǎng)變換漁獲物曲線(xiàn)法計(jì)算總死亡系數(shù)（），并對(duì)其進(jìn)行擬合（圖5）：

圖5 體長(zhǎng)變換漁獲曲線(xiàn)估算南渡江下游海南長(zhǎng)臀鮠總死亡系數(shù)（Z）

ln() = -1.23+ 9.482 (2= 0.96)， （7）

式中，為相對(duì)年齡，為數(shù)量。

擬合曲線(xiàn)的系數(shù)為-1.23。因此，海南長(zhǎng)臀鮠的總死亡系數(shù)（）為1.23，總死亡系數(shù)（）的95%置信區(qū)間為0.99 ～ 1.47。

將南渡江下游海南長(zhǎng)臀鮠的漸近體長(zhǎng)∞（299.25 mm）和生長(zhǎng)參數(shù)（0.41）及采樣點(diǎn)年平均水溫 26.8 ℃代入Pauly[22]經(jīng)驗(yàn)公式計(jì)算出0為-0.09，。利用體長(zhǎng)變換曲線(xiàn)計(jì)算自然死亡系數(shù)為0.51、總死亡系數(shù)為1.23、捕撈死亡系數(shù)得為0.72，開(kāi)發(fā)率為為0.58。

2.4 生長(zhǎng)特性與拐點(diǎn)年齡

通過(guò)對(duì)南渡江下游海南長(zhǎng)臀鮠的生長(zhǎng)方程進(jìn)行求導(dǎo)得海南長(zhǎng)臀鮠的生長(zhǎng)速度及生長(zhǎng)加速度，由圖 6、7可見(jiàn)海南長(zhǎng)臀鮠的生長(zhǎng)變化過(guò)程。

體長(zhǎng)生長(zhǎng)速度dL/d= 122.69 e-0.41 (t +0.09)。

體長(zhǎng)生長(zhǎng)加速度d2L/d2= -50.30 e-0.41 (t +0.09)。

圖6 海南長(zhǎng)臀鮠體長(zhǎng)生長(zhǎng)速度與加速度變化曲線(xiàn)

體質(zhì)量生長(zhǎng)速度dm/d= 574.47 e-0.41 (t + 0.09)×[1- e-0.41(t + 0.09)]2。

體質(zhì)量生長(zhǎng)加速度d2m/d2= 235.53 e-0.41 (t +0.09)×[1- e-0.41 (t + 0.09)] × [3e-0.41 (t + 0.09)- 1]。

圖7 海南長(zhǎng)臀鮠體質(zhì)量生長(zhǎng)速度與加速度變化曲線(xiàn)

南渡江下游海南長(zhǎng)臀鮠體長(zhǎng)的增長(zhǎng)隨年齡的增長(zhǎng)而減小，在高齡時(shí)期逐漸平穩(wěn)趨于0（圖 6）。海南長(zhǎng)臀鮠的體質(zhì)量的生長(zhǎng)速度呈現(xiàn)出先上升后下降的趨勢(shì)。其中在3齡之前，海南長(zhǎng)臀鮠的體質(zhì)量的生長(zhǎng)速度隨年齡的增長(zhǎng)而增長(zhǎng)，3齡后隨年齡增長(zhǎng)而下降，因此推測(cè)南渡江下游的海南長(zhǎng)臀鮠的體質(zhì)量的最大生長(zhǎng)速度出現(xiàn)在2 ～ 3齡之間。通過(guò)計(jì)算得海南長(zhǎng)臀鮠的生長(zhǎng)拐點(diǎn)年齡為3 齡（體長(zhǎng)226.38 mm，體質(zhì)量164.45 g）。

2.5 單位補(bǔ)充量漁獲量

2.5.1 估算平均選擇體長(zhǎng)利用FISAT Ⅱ軟件對(duì)開(kāi)捕體長(zhǎng)進(jìn)行估算，當(dāng)累計(jì)率達(dá)到50%時(shí)對(duì)應(yīng)的體長(zhǎng)為（c）為149.5 mm，開(kāi)捕體長(zhǎng)（c）對(duì)應(yīng)的年齡為1.47。

2.5.2 計(jì)算相對(duì)單位補(bǔ)充量漁獲量 通過(guò)FISAT Ⅱ軟件中的B-H模型計(jì)算，南渡江下游海南長(zhǎng)臀鮠單位補(bǔ)充量漁獲量（/）與開(kāi)發(fā)率（）的關(guān)系如圖所示（圖8）。其中，= 1.24，c/L= 0.50，估算出max= 0.689，0.1= 0.604，0.5= 0.363。

圖8 海南長(zhǎng)臀鮠相對(duì)單位補(bǔ)充量漁獲量、生物量與開(kāi)發(fā)率的關(guān)系

通過(guò)分析南渡江下游海南長(zhǎng)臀鮠的相對(duì)單位補(bǔ)充量等值曲線(xiàn)，可得/與開(kāi)捕體長(zhǎng)c和開(kāi)發(fā)率之間的動(dòng)態(tài)變化趨勢(shì)（圖 9）。

圖9 海南長(zhǎng)臀鮠相對(duì)單位補(bǔ)充量與開(kāi)發(fā)率和開(kāi)捕體長(zhǎng)的關(guān)系

如圖9所示，當(dāng)前南渡江下游海南長(zhǎng)臀鮠漁業(yè)點(diǎn)為A點(diǎn)，且開(kāi)發(fā)率（）為 0.58，c/L為 0.50[即開(kāi)捕年齡（c）為 1.47，對(duì)應(yīng)的開(kāi)捕體長(zhǎng)（c）為 149.5 mm]。南渡江下游海南長(zhǎng)臀鮠理論上的最佳開(kāi)發(fā)狀態(tài)為B點(diǎn)，在/為 1.24 不變的情況下，單位補(bǔ)充量漁獲量/達(dá)到最大值 0.077，其對(duì)應(yīng)的開(kāi)發(fā)率約為1.00，c/L約為 0.80。若從A點(diǎn)調(diào)到B 點(diǎn)，單位補(bǔ)充量漁獲量/將達(dá)到最大值，其對(duì)應(yīng)的開(kāi)發(fā)率增大72%及c/L增大60%。目前海南長(zhǎng)臀鮠開(kāi)捕體長(zhǎng)為149.5 mm，而生長(zhǎng)拐點(diǎn)對(duì)應(yīng)的體長(zhǎng)為226.38 mm，若增加開(kāi)捕體長(zhǎng)，/呈現(xiàn)先增大后減小的趨勢(shì)。當(dāng)c/L為0.7，即開(kāi)捕體長(zhǎng)為209.45 mm時(shí)，/達(dá)到最大值0.06。

3 討論

3.1 參數(shù)估算與資源量評(píng)估

有關(guān)魚(yú)類(lèi)年齡的鑒定的常用材料包括鱗片、耳石等硬組織材料[28]。海南長(zhǎng)臀鮠體表無(wú)鱗片，因此無(wú)法利用鱗片進(jìn)行年齡的鑒定。而利用耳石進(jìn)行年齡的鑒定存在耗時(shí)長(zhǎng)的問(wèn)題，且在進(jìn)行年齡的鑒定過(guò)程中需要較高技能，不同的鑒定者對(duì)于同一個(gè)耳石材料得出的年齡的鑒定結(jié)果一致性的概率小于90%[29]。而長(zhǎng)度頻率法是利用魚(yú)類(lèi)體長(zhǎng)頻率對(duì)生長(zhǎng)參數(shù)進(jìn)行估算，避免利用硬組織鑒定年齡過(guò)程中的主觀(guān)性誤差，因此在漁業(yè)資源的評(píng)估中得到了廣泛應(yīng)用[29-30]。但是由于對(duì)于不同研究對(duì)象和采樣數(shù)據(jù)誤差會(huì)影響模型本身的適用性。因此在利用體長(zhǎng)頻率法（LCA）進(jìn)行評(píng)估時(shí)需要大量的樣品以及足夠的月份，且盡可能獲得大個(gè)體樣本[31]。而本研究利用地籠以及多網(wǎng)目復(fù)合刺網(wǎng)進(jìn)行逐月采樣，所采集樣本體長(zhǎng)頻率分布合理，且大個(gè)體樣本較多（圖2），保證了體長(zhǎng)頻率法估算海南長(zhǎng)臀鮠的準(zhǔn)確性。

魚(yú)類(lèi)生長(zhǎng)參數(shù)和自然死亡系數(shù)是漁業(yè)資源評(píng)估中的兩個(gè)重要參數(shù)[32]。自然死亡系數(shù)受到環(huán)境水溫、魚(yú)類(lèi)生長(zhǎng)狀況等諸多因素的影響[33]。魚(yú)類(lèi)生長(zhǎng)參數(shù)與魚(yú)類(lèi)的漸近體長(zhǎng)L具有負(fù)相關(guān)關(guān)系，與魚(yú)類(lèi)自然死亡系數(shù)具有正相關(guān)關(guān)系[26]。當(dāng)e-< 1時(shí)，表明von Bertalanffy生長(zhǎng)方程能夠反映該魚(yú)類(lèi)的生長(zhǎng)特性[33-34]。而本研究中= 0.41，e-= 0.66 < 1，表明von Bertalanffy可以用來(lái)擬合海南長(zhǎng)臀鮠的生長(zhǎng)。根據(jù)Gulland[35]提出資源最佳開(kāi)發(fā)率為0.5，本研究中南渡江下游海南長(zhǎng)臀鮠的資源開(kāi)發(fā)率為0.58，表明海南長(zhǎng)臀鮠處于過(guò)度捕撈狀態(tài)。根據(jù)Pauly[21]等提出最佳開(kāi)發(fā)率水平（= 0.4），即= 0.16，同樣表明南渡江下游海南長(zhǎng)臀鮠處于過(guò)度開(kāi)發(fā)狀態(tài)。

3.2 南渡江下游海南長(zhǎng)臀鮠資源利用與保護(hù)

魚(yú)類(lèi)種群大小，開(kāi)捕年齡及捕撈強(qiáng)度等因素會(huì)影響漁業(yè)資源的狀態(tài)，進(jìn)而影響漁業(yè)資源的可利用性[36]。其中，捕撈死亡系數(shù)與初次開(kāi)捕體長(zhǎng)是漁業(yè)管理中兩個(gè)重要的人為控制方式，通過(guò)對(duì)這兩個(gè)因素的控制可以實(shí)現(xiàn)最大可持續(xù)產(chǎn)量[19]。本研究通過(guò)B-H動(dòng)態(tài)綜合模型對(duì)南渡江下游海南長(zhǎng)臀鮠漁業(yè)狀態(tài)進(jìn)行分析，結(jié)果表明南渡江下游的海南長(zhǎng)臀鮠的野生資源利用具有不合理性（圖 9），主要為南渡江下游海南長(zhǎng)臀鮠的捕撈強(qiáng)度偏大，而開(kāi)捕年齡則較小。由圖8 可看出保持當(dāng)前開(kāi)發(fā)率= 0.58不變，增大開(kāi)捕體長(zhǎng)至209.45 mm，相對(duì)單位補(bǔ)充量漁獲量可達(dá)到0.06，/增長(zhǎng)1%，與最佳狀態(tài)下的0.077基本持平；而降低南渡江下游的海南長(zhǎng)臀鮠捕撈強(qiáng)度，則會(huì)導(dǎo)致/降低。通過(guò)對(duì)增大開(kāi)捕體長(zhǎng)與增大捕撈強(qiáng)度兩種方式比較，可以得出增大開(kāi)捕體長(zhǎng)比減小捕撈強(qiáng)度更能達(dá)到保護(hù)資源的同時(shí)又合理利用資源的目的。

[1] 中國(guó)水產(chǎn)科學(xué)研究院珠江水產(chǎn)研究所. 海南島淡水及河口魚(yú)類(lèi)志[M].廣州：廣東科技出版社, 1986: 1-406.

[2] 申志新, 李高俊, 蔡杏偉, 等. 海南省淡水野生魚(yú)類(lèi)多樣性演變及保護(hù)建議[J]. 中國(guó)水產(chǎn), 2018(11): 56-60.

[3] 蔡杏偉, 申志新. 海南長(zhǎng)臀鮠[J]. 生物資源, 2020, 42(3): 366.

[4] 蔣志剛, 江建平, 王躍招, 等. 中國(guó)脊椎動(dòng)物紅色名錄[J]. 生物多樣性, 2016, 24(5): 501-551.

[5] 楊家堅(jiān). 海南長(zhǎng)臀鮠人工繁殖與胚胎發(fā)育的研究[J]. 上海水產(chǎn)大學(xué)學(xué)報(bào), 2008, 17(6): 684-688.

[6] 謝少林, 呂子君, 周愛(ài)國(guó), 等. 基于線(xiàn)粒體Cytb和COⅠ基因?qū)χ榻?、海南和紅河長(zhǎng)臀鮠群體遺傳結(jié)構(gòu)分析[J]. 生態(tài)科學(xué), 2016, 35(3): 65-72.

[7] 樂(lè)小亮, 趙爽, 劉海林, 等. 基于線(xiàn)粒體細(xì)胞色素b基因全序列的海南長(zhǎng)臀鮠南渡江種群遺傳變異分析[J]. 生態(tài)科學(xué), 2010, 29(3): 247-250.

[8] GILPIN J A, CHADWICK N E. Life-history traits and population structure of pederson cleaner shrimps[J]. The Biological Bulletin, 2017, 233(3): 190-205.

[9] ASSEFA W W, WONDIE A, ENYEW B G. Population dynamics and exploitation patterns ofin Lake Tana, northwest Ethiopia[J]. Lakes & Reservoirs: Research & Management, 2019, 24(4): 344-353.

[10]LOWE-MCCONNELL R H. Ecological studies in tropical fish com-munities[M]. London: Cambridge University Press, 1987.

[11] 高春霞, 麻秋云, 田思泉, 等. 浙江南部近海小黃魚(yú)生長(zhǎng)、死亡和單位補(bǔ)充量漁獲量[J]. 中國(guó)水產(chǎn)科學(xué), 2019, 26(5): 925-937.

[12] 郭峻宏, 李軍, 沈忱, 等. 閩江口龍頭魚(yú)()的生物學(xué)參數(shù)及資源評(píng)估[J]. 廣東海洋大學(xué)學(xué)報(bào), 2019, 39(5): 56-64.

[13] 葉金清, 柳圭澤, 袁秀堂, 等. 基于生物學(xué)特征的遼河口鮻資源開(kāi)發(fā)與保護(hù)[J]. 大連海洋大學(xué)學(xué)報(bào), 2020, 35(2): 302-307.

[14] 吳金明, 婁必云, 趙海濤, 等. 赤水河魚(yú)類(lèi)資源量的初步估算[J]. 水生態(tài)學(xué)雜志, 2011, 32(3): 99-103.

[15] BEVERTON R J H, HOLT S J. On the dynamics of exploited fish populations[M]. Dordrecht: Springer Netherlands, 1993.

[16] RICKER W E. Handbook of computations for biological statistics of fish populations [J]. Can Fish Res Board Bull, 1958: 119-300.

[17] 中國(guó)水產(chǎn)科學(xué)研究院珠江水產(chǎn)研究所, 華南師范大學(xué), 暨南大學(xué), 上海水產(chǎn)大學(xué), 湛江水產(chǎn)學(xué)院. 廣東淡水魚(yú)類(lèi)志[M]. 廣東: 廣東科技出版社, 1991: 1-622.

[18] 鄭慈英. 珠江魚(yú)類(lèi)志[M]. 北京: 科學(xué)技術(shù)出版社, 1989: 1-452.

[19] 詹秉義. 漁業(yè)資源評(píng)估[M]. 北京: 中國(guó)農(nóng)業(yè)出版社, 1995: 1-353.

[20] GAYANILO F C, SPARRE P, PAULY D. FAO-ICLARM Stock Assessment Tools Ⅱ (Fisat Ⅱ): User’s Guide (Revised Version) [M]. Rome: FAO Computerized Information Series (Fisheries), 2005: 1-152.

[21] PAULY D. Length-converted catch curves and the seasonal growth of fishes[J]. Fishbyte, 1990, 8(3): 24-29.

[22] PAULY D, Morgan G. Length-based Methods in Fisheries Research [M]. World Fish1987: 5-89.

[23] VON BERTALANFFY L. A quantitative theory of organic growth (Inquires on growth laws.Ⅱ) [J]. Human Biology, 1938, 10: 181-213.

[24] SPARRE P, WILLMANN R. BEAM 4, a bio-economic multi-species, multi-fleet, multi-plant, multi-area extension of the traditional forecast model[M]. Copenhagen: ICES, Statistical Committee. 1992.

[25] 何寶全, 李輝權(quán). 珠江河口棘頭梅童魚(yú)的資源評(píng)估[J]. 水產(chǎn)學(xué)報(bào), 1988, 12(2): 125-134.

[26] PAULY D, SORIANO M L. Some practical extensions to Beverton and Holt’s relative yield-per-recruit model[J]. First Asian Fisheries Forum, 1986(5): 26-31.

[27] 田輝伍, 岳興建, 陳大慶, 等. 怒江東方墨頭魚(yú)資源量和死亡參數(shù)的估算[J]. 生態(tài)學(xué)雜志, 2012, 31(1): 235-240.

[28] 張學(xué)健, 程家驊. 魚(yú)類(lèi)年齡鑒定研究概況[J]. 海洋漁業(yè), 2009, 31(1): 92-99.

[29] 高春霞, 田思泉, 戴小杰. 淀山湖刀鱭的生物學(xué)參數(shù)估算及其相對(duì)單位補(bǔ)充量漁獲量[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2014, 25(5): 1506-1512.

[30] CUNHA N L, CATELLA A C, KINAS M A. Growth parameters estimates for a small fish of the Pantanal, Brazil:(Characiformes; Characidae)[J]. Brazilian Journal of Biology, 2007, 67(2): 293-297.

[31] HOENING, J. M., CSIRKE, J. SAN DE RS, M. J, et al. Data acquisition for length-based stock assessment: report of Working Group I [C]. Italy, International Centre for Living Aquatic Resources Management, 1987:343-352.

[32] 朱國(guó)平, 宋旗. 運(yùn)用脂褐素鑒定甲殼類(lèi)年齡的研究進(jìn)展[J]. 生態(tài)學(xué)雜志, 2016, 35(8): 2225-2233.

[33] STERGIOU K I, VOULGARIDOU P. Trends in various biological parameters of the European sardine,(Walbaum, 1792), in the Eastern Mediterranean Sea[J]. Scientia Marina, 2003, 67(S1): 269-280.

[34] 費(fèi)鴻年, 張?jiān)娙? 水產(chǎn)資源學(xué)[M]. 北京: 中國(guó)科學(xué)技術(shù)出版社, 1990: 1-787.

[35] GULLAND J A. fish stock assessment: a manual of basic methods [M]. Chichester: Wiley Interscience, 1983.56-76.

[36] 王繼隆, 劉偉, 魯萬(wàn)橋, 等. 黑龍江中游烏蘇里白鮭資源現(xiàn)狀評(píng)估[J]. 生態(tài)學(xué)雜志, 2019, 38(6): 1824-1829.

Estimation of Biological Parameters ofin the Downstream of Nandu River

LIU fei1, LIANG Zhi-ce1, PAN Shuai1, SUN An-fu1, CAI Xing-wei2, GUO Zhi-qiang1

（1.//,570228,; 2.,571126,）

【】To perform the biology and population resource assessment ofand to provided basic information for the population protection and resource recovery ofin the downstream of Nandu River. 【】A total of 677 fish samples were collected from the downstream of Nandu River from October 2019 to September 2020 monthly. The population parameters ofwere calculated by the ELEFAN I part of FiSAT Ⅱ software, and the population dynamic ofwas estimated by Beverton-Holt dynamic model.【】The length()-mass() relationship ofin the downstream of Nandu River was=0.000 0073.16(2= 0.97;< 0.05); Von Bertalanffy growth equation parameters ofwere asymptotic length∞= 299.25 mm, curvature parameter= 0.41,0(theoretical age when the theoretical length of fish was 0) =-0.09years. Calculated by growth equation, the inflection age was 3 years, inflection body length was226.38 mm, inflection mass was164.45 g for.The natural mortality calculated by Pauly’s empirical function was 0.51, the total mortality() coefficient, fishing mortality () and exploitation rate() estimated by length-converted catch curve method were 1.23，0.72, and 0.58,respectively. The capture age and body length ofcalculated by dynamic pool model were 1.47 years and 149.5 mm, respectively.【】in the downstream of Nandu River is overexploited.

; biological parameters; overfishing; Nandu River

S932.4

A

1673-9159（2021）06-0131-07

10.3969/j.issn.1673-9159.2021.06.016

劉飛，梁智策，潘帥，等. 南渡江下游海南長(zhǎng)臀鮠生物學(xué)參數(shù)估算[J]. 廣東海洋大學(xué)學(xué)報(bào)，2021，41（6）：131-137.

2021-06-30

2020年海南省重點(diǎn)研發(fā)計(jì)劃（ZDYF2020193）；國(guó)家自然基金（31872574, 41866003, 31960731）

劉飛（1999－），男，碩士研究生，主要研究方向?yàn)闈O業(yè)資源保護(hù)。E-mail: liufei@hainanu.edu.cn

郭志強(qiáng)，E-mail: guozq@hainanu.edu.cn

（責(zé)任編輯：劉嶺）