熊 飛，劉紅艷，段辛斌，劉紹平，陳大慶

1 江漢大學(xué)生命科學(xué)學(xué)院， 武漢 430056 2 中國水產(chǎn)科學(xué)研究院長江水產(chǎn)研究所， 武漢 430223

長江上游朱楊溪江段圓筒吻鮈種群參數(shù)和資源量

熊 飛1，劉紅艷1，段辛斌2, *，劉紹平2，陳大慶2

1 江漢大學(xué)生命科學(xué)學(xué)院， 武漢 430056 2 中國水產(chǎn)科學(xué)研究院長江水產(chǎn)研究所， 武漢 430223

圓筒吻鮈為長江上游特有種，三峽、向家壩、溪洛渡等梯級水利工程的建設(shè)，可能對其物種生存產(chǎn)生較大影響。為了解三峽工程蓄水后、金沙江一期工程蓄水前該物種的種群動態(tài)，2007—2009年對長江上游朱楊溪江段的圓筒吻鮈進(jìn)行了調(diào)查，利用體長頻率數(shù)據(jù)對其種群參數(shù)和資源量進(jìn)行了估算。長江上游朱楊溪江段圓筒吻鮈漁獲群體體長范圍為69.0—268.0 mm，體重范圍為3.9—230.4 g，平均體長為(181.4±26.9)mm，平均體重為(78.5±33.0)g；優(yōu)勢體長組為125.0—225.0 mm，約占總數(shù)的93.0%。圓筒吻鮈體長體重關(guān)系為：W=1.58×10-5L2.95(R2=0.94，P<0.01,n=401)，為勻速生長類型，可用Von Bertalanffy生長方程描述其生長規(guī)律。圓筒吻鮈生長參數(shù)為L∞=361 mm、k=0.21 a-1、t0=-0.68；死亡系數(shù)為Z=1.70、M=0.50、F=1.20。其開捕體長為163.4 mm，資源開發(fā)率為0.70，超過了其資源最適開發(fā)率(0.60)而接近最大開發(fā)率(0.75)。2007—2009年朱楊溪江段圓筒吻鮈年漁獲量分別為6716 尾(0.50 t)、22772 尾(1.87 t)和16139 尾(1.20 t)，平均值為15209 尾(1.19 t)。由體長股分析法估算出2007—2009年朱楊溪江段圓筒吻鮈年資源量分別為16361 尾/km(1.25 t/km)、13922 尾/km(1.74 t/km)和26836 尾/km(1.93 t/km)，均值為19040 尾/km(1.64 t/km)。目前圓筒吻鮈資源開發(fā)率偏高，建議提高開捕體長至193.5 mm，降低其開發(fā)率。建議逐步引導(dǎo)漁民轉(zhuǎn)產(chǎn)轉(zhuǎn)業(yè)，實施全年禁漁措施，加強(qiáng)筒吻鮈等特有魚類的保護(hù)。

種群參數(shù)；生長；死亡；資源評估；圓筒吻鮈

長江上游擁有豐富的珍稀、特有魚類資源，是我國生物多樣性保護(hù)的關(guān)鍵區(qū)域之一。三峽大壩等大型水利工程的建設(shè)，將阻礙魚類洄游通道、改變魚類棲息生境，對長江上游魚類群落產(chǎn)生顯著影響。三峽水庫蓄水后，使長江上游約600 km干流江段形成河谷型水庫，水流變緩，原來棲息于此的多種特有魚類將逐漸從庫區(qū)消失，其棲息生境大范圍萎縮[1]。溪洛渡、向家壩等更多梯級電站的建設(shè)，將加劇這種影響，長江上游特有魚類的種群生存令人擔(dān)憂，其種群動態(tài)頗受關(guān)注[2]。

魚類的生長參數(shù)和死亡系數(shù)是漁業(yè)資源評估和管理的關(guān)鍵參數(shù)，是了解魚類種群數(shù)量變動的重要因素，亦是合理開發(fā)利用魚類資源及制定魚類資源保護(hù)策略的基礎(chǔ)。以往魚類種群參數(shù)和資源量估算主要依據(jù)魚類的年齡-體長數(shù)據(jù)，需要進(jìn)行年齡鑒定[3]。而利用體長頻率數(shù)據(jù)來進(jìn)行種群參數(shù)和資源量估算可以避免年齡鑒定，具有一定的優(yōu)勢[4]。隨著FiSATⅡ漁業(yè)評估軟件的不斷發(fā)展和完善，可以依據(jù)魚類的體長頻率數(shù)據(jù)估算魚類的生長參數(shù)、死亡系數(shù)和資源量等，從而提高了魚類資源評估的效率和精度[5]，該方法在魚類種群參數(shù)和資源量估算方面正在逐步得到推廣應(yīng)用[6-7]。我國這方面的研究和應(yīng)用主要集中在海洋魚類方面[8-9]，對河流魚類的研究相對較少[10]，遠(yuǎn)遠(yuǎn)不能滿足河流魚類資源管理和保護(hù)的需求。

圓筒吻鮈(Rhinogobiocylindricus)隸屬于鯉科(Cyprinidae)鮈亞科(Gobioninae)吻鮈屬(Rhinogobio)，主要分布于長江上游干、支流，是長江上游的特有種。目前已有一些關(guān)于圓筒吻鮈年齡和生長的報道[11-13]，但對其生長參數(shù)的估算主要依據(jù)年齡-體長數(shù)據(jù)。本研究根據(jù)2007—2009年長江上游江津朱楊溪江段的漁業(yè)調(diào)查資料，利用體長頻率數(shù)據(jù)對其種群參數(shù)和資源量進(jìn)行估算，以了解三峽工程蓄水后、金沙江一期工程蓄水前其種群動態(tài)和資源利用狀況，為其資源管理和物種保護(hù)提供依據(jù)，也可為大型水利工程的生態(tài)學(xué)效應(yīng)評估積累資料。

1 材料與方法

1.1 漁業(yè)調(diào)查

圖1 長江上游朱楊溪江段位置及樣點設(shè)置 Fig.1 Location of Zhuyangxi section in the upper Yangtze River and sampling sites

2007—2009年對長江上游江津朱楊溪江段的圓筒吻鮈進(jìn)行漁業(yè)資源調(diào)查與取樣，調(diào)查時間主要集中在5—12月。調(diào)查范圍為朱楊鎮(zhèn)上游至石門鎮(zhèn)約15 km的江段，主要調(diào)查站點包括朱楊、塘河口和羅灣壩(圖1)。調(diào)查漁具主要包括流刺網(wǎng)(網(wǎng)目3—6 cm)和百袋網(wǎng)(網(wǎng)目2 cm)。物種鑒定主要參考《四川魚類志》[14]。對圓筒吻鮈的體長和體重進(jìn)行測量，體長從吻端測量到尾部前端的最后一枚鱗片，精確到1 mm，體重測量精確到1 g。對調(diào)查江段的捕撈船只及單位捕撈努力量漁獲量進(jìn)行統(tǒng)計，根據(jù)作業(yè)江段捕撈的平均強(qiáng)度，流刺網(wǎng)和百袋網(wǎng)每天作業(yè)網(wǎng)次按10 網(wǎng)的標(biāo)準(zhǔn)計算單位捕撈努力量漁獲量。

1.2 種群參數(shù)估算

體長體重關(guān)系依據(jù)Keys公式W=aLb進(jìn)行擬合，用t檢驗法檢驗冪指數(shù)b值與3的偏離程度，判斷其生長是否為勻速生長類型[15-16]。利用Von Bertalanffy生長方程描述生長規(guī)律，漸近體長(L∞)、生長系數(shù)(k)和生長特征指數(shù)(φ)根據(jù)魚類體長頻率的時間系列數(shù)據(jù)借助FiSATⅡ漁業(yè)評估軟件中的體長頻率分析法(ELEFANⅠ)模塊進(jìn)行估算，總死亡系數(shù)(Z)、開捕體長(L50)由體長轉(zhuǎn)漁獲物曲線(Converted-length catch curve)模塊進(jìn)行估算[3]。理論生長起點年齡(t0)、自然死亡系數(shù)(M)根據(jù)Pauly經(jīng)驗公式獲得[17]：

ln(-t0)= -0.3922-0.2752ln(L∞)-1.038ln(k)

lnM=-0.0152-0.279ln(L∞)+0.6543lnk+0.463lnT

式中,T表示棲息江段的平均水溫，本研究中參考相關(guān)文獻(xiàn)取18.4 ℃[18]。魚類的資源開發(fā)率E=F/Z。采用Beverton-Holt動態(tài)綜合模型的相對單位補充量產(chǎn)量(Y′/R)和相對單位補充量生物量(B′/R)曲線評價圓筒吻鮈的資源利用情況，Emax表示能獲得最大產(chǎn)量時的開發(fā)率，E10表示資源群體邊際產(chǎn)量減少到理論原始資源邊際產(chǎn)量1/10時的開發(fā)率，被認(rèn)為是“最適”開發(fā)率，E50表示資源量下降到原始水平50%時的開發(fā)率[3]。

1.3 漁獲量計算

圓筒吻鮈年漁獲量(Y)由各種漁具的年漁獲量及其比例推算，計算公式為：

Y=∑(Yi×Pi)

式中,Pi為圓筒吻鮈在漁具i漁獲物中所占的數(shù)量百分?jǐn)?shù)，Yi為漁具i的年漁獲量，計算公式為：

Yi=xi×fi×ti

式中,xi為漁具i的單位捕撈努力量漁獲量(尾 船-1d-1)，fi為漁具i的作業(yè)船只數(shù)，ti為漁具i的年作業(yè)時間(d)，除去春季禁漁時間(2—4 月)和洪水季節(jié)(7—8 月)及節(jié)假日休息，本研究中各漁具的作業(yè)天數(shù)統(tǒng)一按200 d計算。

1.4 資源量估算

資源量估算借助FiSATⅡ中的體長股分析(Length-structured VPA)模塊進(jìn)行分析[5]。將圓筒吻鮈的年漁獲量按體長分組錄入模塊中，輸入L∞、k、M、a、b及Ft(最大體長組捕撈死亡系數(shù))等參數(shù)，運算后即可得各體長組的捕撈死亡系數(shù)和資源量。其中，F(xiàn)t采用迭代法確定，初始值取0.5[19]。各體長組資源量之和即為圓筒吻鮈的年資源量。

2 結(jié)果與分析

2.1 群體結(jié)構(gòu)

長江上游朱楊溪江段圓筒吻鮈漁獲群體體長范圍為69.0—268.0 mm，平均體長為(181.4±26.9)mm。優(yōu)勢體長組為125.0—225.0 mm，約占總數(shù)的93.0%，125.0 mm體長以下和225.0 mm體長以上的個體分別只占3.4%和3.7%(圖2)。圓筒吻鮈體重范圍為3.9—230.4 g，平均體重為(78.5±33.0)g。個體體重主要集中在130 g以下，約占個體總數(shù)的93.5%，150 g以上個體僅占3.7%。

圖2 長江上游朱楊溪江段圓筒吻鮈體長組成Fig.2 Length composition of Rhinogobio cylindricus in Zhuyangxi section of the upper Yangtze River

圖3 長江上游朱楊溪江段圓筒吻鮈長度和體重的關(guān)系 Fig.3 Relationship between length and weight of Rhinogobio cylindricus in Zhuyangxi section of the upper Yangtze River

2.2 生長參數(shù)

圓筒吻鮈體長(L)與體重(W)呈顯著冪函數(shù)關(guān)系(圖3)：W=1.58×10-5L2.95(R2=0.94，P<0.01，n=401)。體長與體重關(guān)系的冪指數(shù)與3均無顯著性差異(t=1.14

由ELEFANⅠ擬合圓筒吻鮈的生長參數(shù)，當(dāng)L∞=361 mm，k=0.21 a-1時擬合效果最佳，此時t0=-0.68，φ=4.44。拐點年齡為4.48 a，對應(yīng)的體長和體重為分別為238.8 mm和166.6 g。由ELEFANⅠ估計的圓筒吻鮈生長曲線見圖4，其體長和體重的Von Bertalanffy生長方程為：

Lt=361[1-e-0.21(t+0.68)]

Wt=564.8[1-e-0.21(t+0.68)]2.95

圖4 圓筒吻鮈的體長頻率序列數(shù)據(jù)及應(yīng)用ELEFANⅠ估計的生長曲線Fig.4 Length frequency of Rhinogobio cylindricus and growth curve estimated by ELEFANⅠ

2.3 死亡系數(shù)和開發(fā)率

圓筒吻鮈的長度轉(zhuǎn)漁獲物曲線見圖5，據(jù)此估算出其總死亡系數(shù)Z=1.70。其自然死亡系數(shù)M=0.50，捕撈死亡系數(shù)F=1.20，開發(fā)率E=0.70。不同年份的死亡系數(shù)和開發(fā)率見表1。由長度轉(zhuǎn)漁獲物曲線估算出圓筒吻鮈的開捕體長L50=163.4 mm。由Beverton-Holt動態(tài)綜合模型的相對補充量產(chǎn)量(生物量)曲線估算出，開捕體長為163.4 mm時，E10=0.60、E50=0.35、Emax=0.75(圖6)。圓筒吻鮈目前的資源開發(fā)率超過了其最適開發(fā)率(0.60)而接近其最大開發(fā)率(0.75)。

表1 長江上游朱楊溪江段圓筒吻鮈的死亡系數(shù)和開發(fā)率Table 1 Mortality rates and exploitation rates of Rhinogobio cylindricus in Zhuyangxi section of the upper Yangtze River

圖5 根據(jù)長度轉(zhuǎn)漁獲物曲線估算圓筒吻鮈的總死亡系數(shù)Fig.5 Length-converted catch curves of Rhinogobio cylindricus

圖6 圓筒吻鮈的相對單位補充量產(chǎn)量(Y′/R)和相對單位補充生物量(B′/R)曲線 Fig.6 Relatvie yield per recruit (Y′/R) and Relatvie biomass per recruit (B′/R) curve of Rhinogobio cylindricus

2.4 資源量

漁獲物分析表明，圓筒吻鮈在流刺網(wǎng)和百袋網(wǎng)漁獲物中的數(shù)量百分?jǐn)?shù)平均值分別為12.8%和2.6%。根據(jù)單位捕撈努力量漁獲量和圓筒吻鮈的比例估算出2007—2009年朱楊溪江段圓筒吻鮈年漁獲量分別為6716 尾(0.50 t)、22772 尾(1.87 t)和16139 尾(1.20 t)，年均值為15209 尾(1.19 t)。

利用Length-structured VPA模塊估算圓筒吻鮈的資源量，通過迭代法計算出2007—2009年其最大體長組的捕撈死亡系數(shù)Ft分別為0.1194、0.5030和0.1818。結(jié)果表明(圖7)，個體體長在160mm以下時，資源群體的損失以自然死亡為主，而在160mm以上時，資源群體的損失以捕撈死亡為主。2007—2009年朱楊溪江段圓筒吻鮈年資源量分別為245411、208836和402537 尾，年均值為285594 尾；以生物量計，年資源量分別為18.80、26.09和28.94t，年均值為24.61t。以單位長度河段來表示，2007—2009年圓筒吻鮈資源量分別為16361 尾/km(1.25t/km)、13922 尾/km(1.74t/km)和26836 尾/km(1.93t/km)，平均值為19040 尾/km(1.64t/km)。

圖7 2007—2009年圓筒吻鮈的體長結(jié)構(gòu)實際種群分析Fig.7 Length structured VPA for Rhinogobio cylindricus in 2007—2009

3 討論

3.1 種群參數(shù)和資源量估算

根據(jù)體長頻率數(shù)據(jù)進(jìn)行魚類種群參數(shù)和資源量估算，可以避免年齡鑒定的主觀誤差而影響參數(shù)評估的準(zhǔn)確性，特別適合年齡鑒定比較困難的評估對象，在漁業(yè)資源評價中取得了較好的應(yīng)用效果[20- 22]。ELEFAN技術(shù)估算魚類生長參數(shù)主要根據(jù)魚類體長頻率的時間系列數(shù)據(jù)進(jìn)行體長頻率樣品重構(gòu)，將體長頻率形成的波峰和波谷放大，以擬合最佳生長曲線，從而估算出合適的生長參數(shù)，該方法對采樣頻率和樣品容量要求較高[23]。本研究由于春季禁漁和魚類種群數(shù)量偏少等原因，采樣頻率和樣品容量都受到了一定限制，但圓筒吻鮈體長數(shù)據(jù)(69.0—268.0mm)基本包括了從小到大的各個體長段的個體，具有一定的代表性。以往對圓筒吻鮈生長參數(shù)的估算主要依據(jù)體長-年齡數(shù)據(jù)，估算結(jié)果表明，長江上游圓筒吻鮈的L∞為336—389mm，k為0.18—0.21a-1，φ為4.37—4.44。本研究利用FiSATⅡ軟件中的ELEFANⅠ方法估算出L∞、k、φ分別為361mm、0.21a-1和0.44，L∞值位于以往研究結(jié)果范圍之內(nèi)，k和φ值與以往研究結(jié)果接近，兩種方法估算的結(jié)果較為一致(表2)。自然死亡系數(shù)M值估算結(jié)果的準(zhǔn)確性受評估模型的影響較大，考慮到本研究中已由ELEFAN技術(shù)估算出了L∞和k值，直接采用了比較通用的Pauly經(jīng)驗公式進(jìn)行估算[17]。一般認(rèn)為M/k值應(yīng)在1.5—2.5范圍內(nèi)[3]，本研究中M/k值為2.4，在合理范圍之內(nèi)。

表2 圓筒吻鮈種群參數(shù)估算結(jié)果比較Table 2 Comparison of estimated population parameter of Rhinogobio cylindricus

魚類資源量評估有多種方法，如標(biāo)志放流法、初級生產(chǎn)力估算法、魚卵仔魚調(diào)查法、水聲學(xué)探測法、年齡結(jié)構(gòu)股分析法和體長結(jié)構(gòu)股分析法等[3]，其中，體長股分析法不需要鑒定研究對象的年齡，具有采樣工作量小和數(shù)據(jù)結(jié)構(gòu)簡單等優(yōu)點而受到青睞[24-25]。傳統(tǒng)的體長股分析中, 各體長組的Ft值是根據(jù)最大體長組的開發(fā)率估算的，不存在迭代過程，一般將最大體長組開發(fā)率定為0. 5，但實際上最大體長最的開發(fā)率可能與0. 5差異較大，從而影響估算結(jié)果的準(zhǔn)確性；而通過FiSATⅡ中的體長股分析模塊可以輕松使用迭代法確定最大體長組的捕撈死亡系數(shù)Ft，從而提高了魚類資源量估算的效率和精度[10,19]。本研究通過該方法，在種群參數(shù)估算的基礎(chǔ)上，對長江上游圓筒吻鮈的資源量進(jìn)行了評估，結(jié)果表明，2007—2009年圓筒吻鮈資源量分別為16361 尾/km(1.25 t/km)、13922 尾/km(1.74 t/km)和26836 尾/km(1.93 t/km)，該結(jié)果可為今后評估圓筒吻鮈的種群數(shù)量長期變化規(guī)律提供基礎(chǔ)。

3.2 資源利用與保護(hù)

開捕體長和捕撈死亡系數(shù)是影響漁業(yè)資源量和漁獲量的兩個重要因素，也是制定漁業(yè)管理措施的兩個主要依據(jù)[1]。本研究表明，圓筒吻鮈生長拐點對應(yīng)的體長為238.8 mm，即生長速度在該體長以前逐漸增加，之后呈下降趨勢，而圓筒吻鮈的開捕體長為163.4 mm，表明其還未充分生長就被捕撈。Froese和Binohlan提出能獲得最大相對漁獲量的最適捕撈體長(Lopt)可由其最小性成熟體長(Lm)估算獲得：logLopt=1.053Lm-0.0565[26]，圓筒吻鮈最小性成熟年體長約為168.0 mm[11]，由此估算出最適開捕體長為193.5 mm。因此，將圓筒吻鮈的開捕體長由目前的163.4 mm提高到193.5 mm，可讓圓筒吻鮈充分生長，獲得較高的相對漁獲量。

由魚類的死亡系數(shù)可以計算出魚類的開發(fā)率E，該參數(shù)反映了種群資源的利用程度。本研究表明，圓筒吻鮈現(xiàn)有開發(fā)率為0.70，與熊星等的研究結(jié)果一致[13]。一般認(rèn)為魚類的最適開發(fā)率為0.5[27]，若按此標(biāo)準(zhǔn)，圓筒吻鮈資源處于過度利用狀態(tài)。在目前的開捕體長(L50=163.4 mm)下，圓筒吻鮈的開發(fā)率已接近其最大開發(fā)率(Emax=0.75)，理論上在此開發(fā)率下可獲得較大的產(chǎn)量，但從生物學(xué)和經(jīng)濟(jì)學(xué)等多方面考慮，應(yīng)該以最適產(chǎn)量和最適開發(fā)率為管理目標(biāo)，可考慮降低其現(xiàn)有開發(fā)率至最適開發(fā)率(E10=0.60)以下。

隨著三峽、向家壩和溪洛渡等大型梯級水利工程的建設(shè)，長江上游魚類棲息生境已支離破碎，適應(yīng)流水生境的圓筒吻鮈的棲息生境萎縮，其壩上和壩下種群也得不到有效交流。圓筒吻鮈為產(chǎn)漂流性卵魚類[28]，其產(chǎn)出的卵隨水漂流發(fā)育，必須要保證足夠的漂流距離才能發(fā)育成具有主動游泳能力的魚苗，這些梯級電站的建設(shè)，將影響其產(chǎn)卵場的分布和卵、苗的發(fā)育，從而影響其種群生存。因此，應(yīng)采取補救措施，加強(qiáng)對圓筒吻鮈的保護(hù)力度。目前，江津至宜賓干流江段仍保持著適宜圓筒吻鮈等特有魚類生存的流水生境，在長江上游特有魚類保護(hù)中具有重要意義。從珍稀、特有物種資源保護(hù)的角度考慮，應(yīng)逐步引導(dǎo)漁民轉(zhuǎn)產(chǎn)轉(zhuǎn)業(yè)，實施全年禁漁措施，加強(qiáng)長江上游珍稀特有魚類國家級自然保護(hù)區(qū)的建設(shè)和管理，促進(jìn)對圓筒吻鮈等特有物種的有效保護(hù)。

[1] Park Y S, Chang J B, Lek S, Cao W X, Brosse S. Conservation strategies for endemic fish species threatened by the Three Gorges Dam. Conservation Biology, 2003, 17(6): 1748- 1758.

[2] 曹文宣. 長江上游特有魚類自然保護(hù)區(qū)的建設(shè)及相關(guān)問題的思考. 長江流域資源與環(huán)境, 2000, 9(2): 131- 132.

[3] 詹秉義. 漁業(yè)資源評估. 北京: 中國農(nóng)業(yè)出版社, 1995.

[4] Wang Y G, Ellis N. Maximum likelihood estimation of mortality and growth with individual variability from multiple length-frequency data. Fishery Bulletin, 2005, 103: 380- 391.

[5] Gayanilo F C Jr, Sparre P, Pauly D. FAO-ICLARM Stock Assessment Tools II (User′s Guide). Worldfish Center, FAO, Rome, 2005.

[6] Ghosh S, Sivadas M, Abdussamad E M, Rohit P, Koya K P S, Joshi K K, Chellappan A, Rathinam M M, Prakasan D, Sebastine M. Fishery, population dynamics and stock structure of frigate tunaAuxisthazard(Lacepede, 1800) exploited from Indian waters. Indian Journal of Fisheries, 2012, 59(2): 95- 100.

[7] Johnson M G, Tamatamah A R. Length frequency distribution, mortality rate and reproduction biology of Kawakama (Euthynnusaffinis-Cantor, 1849) in the coastal waters of Tanzania. Pakistan Journal of Biological Sciences, 2013, 16(21): 1270- 1278.

[8] 陳作志, 邱永松, 黃梓榮. 南海北部白姑魚生長和死亡參數(shù)的估算. 應(yīng)用生態(tài)學(xué)報, 2005, 16(4): 712- 716.

[9] 劉元文, 薛瑩, 魏邦福, 任一平. 海州灣5種主要魚類種群參數(shù)估算. 中國水產(chǎn)科學(xué), 2014, 21(1): 125- 133.

[10] 吳金明, 婁必云, 趙海濤, 黎良, 陳永祥, 王劍偉. 赤水河魚類資源量的初步估算. 水生態(tài)學(xué)雜志, 32(3): 99- 103.

[11] 馬惠欽. 長江干流圓筒吻鮈生物學(xué)的初步研究 [D]. 重慶: 西南大學(xué), 2001.

[12] 王美榮, 楊少榮, 劉飛, 黎明政, 但勝國, 劉煥章. 長江上游圓筒吻鮈年齡與生長的研究. 水生生物學(xué)報, 2012, 36(2): 262- 269.

[13] 熊星. 圓筒吻鮈(RhinogobiocylindricusGünther)基礎(chǔ)生物學(xué)及資源量研究 [D]. 重慶: 重慶師范大學(xué), 2013.

[14] 丁瑞華. 四川魚類志. 成都: 四川科學(xué)技術(shù)出版社, 1994.

[15] Richer W E. Computation and Interpretation of Biological Statistics of Fish Populations (Bulletin 191 of the Fisheries Research Board of Canada). Caldwell, New Jersey: Blackburn Press, 1975.

[16] Pauly D. Fish Population Dynamics in Tropical Waters: A Manual for Use with Programmable Calculators (ICLARM Studies and Reviews 8). Manila, Philippines: International Center for Living Aquatic Resources Management, 1984.

[17] Pauly D. On the interrelationships between natural mortality, growth parameters, and mean environmental temperature in 175 fish stocks. ICES Journal of Marine Science, 1980, 39(2): 175- 192.

[18] 田輝伍, 段辛斌, 熊星, 羅宏偉, 劉紹平, 陳大慶. 長江上游長薄鰍生長和種群參數(shù)的估算. 長江流域資源與環(huán)境, 2013, 22(10): 1305- 1312.

[19] 吳斌, 方春林, 賀剛, 傅培峰. FiSAT II軟件支持下的體長股分析法探討. 南方水產(chǎn)科學(xué), 2013, 9(4): 94- 98.

[20] 陳國寶, 李永振, 陳丕茂, 舒黎明. 魚類最佳體長頻率分析組距研究. 中國水產(chǎn)科學(xué), 2008, 15(4): 659- 666.

[21] Wang X H, Qiu Y S, Du F Y, Lin Z J, Sun D R, Huang S L. Population parameters and dynamic pool models of commercial fishes in the Beibu Gulf, northern South China Sea. Chinese Journal of Oceanology and Limnology, 2012, 30(1): 105- 117.

[22] Liu Q, Xu B D, Ye Z J, Ren Y P. Growth and mortality of small yellow croaker (Larimichthyspolyactis) inhabiting Haizhou Bay of China. Journal of Ocean University of China, 2012, 11(4): 1672- 5182.

[23] 李壯, 劉群. 應(yīng)用漁業(yè)體長分析方法ELEFAN和SLCA估算魚類生長參數(shù)的研究. 海洋湖沼通報, 2007, (3): 81- 87.

[24] 劉軍, 王劍偉, 苗志國, 曹文宣. 長江上游宜賓江段長鰭吻鮈種群資源量的估算. 長江流域資源與環(huán)境, 2010, 19(3): 276- 280.

[25] 周永東, 徐漢祥, 潘國良, 賀舟挺. 東海區(qū)鯧魚、小黃魚資源量及其持續(xù)漁獲量的估算. 浙江海洋學(xué)院學(xué)報: 自然科學(xué)版, 2013, 32(1): 1- 5.

[26] Froese R, Binohlan C. Empirical relationships to estimate asymptotic length, length at first maturity and length at maximum yield per recruit in fishes, with a simple method to evaluate length frequency data. Journal of Fish Biology, 2000, 56(4): 758- 773.

[27] Gulland J A. The Fish Resources of the Ocean. London: Fishing News (Books), Ltd., 1971.

[28] 李修峰, 黃道明, 謝文星, 謝山, 常秀嶺, 楊漢運, 張友謙, 何家慶. 漢江中游產(chǎn)漂流性卵魚類產(chǎn)卵場的現(xiàn)狀. 大連水產(chǎn)學(xué)院學(xué)報, 2006, 21(2): 105- 111.

Population parameters and population abundance ofRhinogobiocylindricusin Zhuyangxi section of the upper Yangtze River

XIONG Fei1, LIU Hongyan1, DUAN Xinbin2, *, LIU Shaoping2, CHEN Daqing2

1SchoolofLifeSciences,JianghanUniversity,Wuhan430056,China2YangtzeRiverFisheriesResearchInstitute,ChineseAcademyofFisheryScience,Wuhan430223,China

The upper Yangtze River supports many rare and endemic fish species and plays a critical role in biodiversity conservation in China. However, the fish fauna is under serious threat from large cascaded hydroelectric projects, including damming in Three Gorges, Xiangjiaba, and Xiluodu. The state of the endemic species is a major ecological concern. Estimates of population abundance and demography are needed for assessing fish population dynamics and effective fisheries management. A traditional method of estimating population parameters is based on the age-length relationship, which requires detailed age determination. A simpler alternative method is based on length-frequency data. We employed the second method and used the software FiSATⅡ (FAO-ICLARM Stock Assessment Tools) to estimate the growth rate and mortality ofRhinogobiocylindricus, an endemic species of the upper Yangtze River, based on a survey in the Zhuyangxi section conducted in 2007—2009, a period after the Three Gorges Reservoir was impounded, but before the Xiangjiaba and Xiluodu reservoirs were. Growth and mortality parameters of the fish were estimated based on length-frequency data.R.cylindricusranged from 69.0 to 268.0 mm in length and 3.9 to 230.4 g in weight, with an average length of (181.4±26.9) mm and an average weight of (78.5±33.0) g. The length group of 125.0—225.0 mm dominated the catches (93.0% of the total number). The length-weight relationship was well fitted with a power function,W= 1.58 × 10-5L2.95(R2= 0.94,P< 0.01,n= 401), a result that means that the growth of the fish is isometric and adequately described by the Von Bertalanffy growth equation. Three growth-related population parameters,L∞,k, andt0, were estimated as 361 mm, 0.21 a-1, and -0.68, respectively. The three mortality-related parameters,Z,M, andF, were estimated as 1.70, 0.50, and 1.20, respectively. The minimum catchable size was 163.4 mm, and the exploitation rate observed was 0.70, higher than the estimated optimum exploitation rate (0.60), but lower than the maximum exploitation rate (0.75). The annual yield ofR.cylindricusin Zhuyangxi section was 6716 ind. (0.50 t) in 2007, 22772 ind. (1.87 t) in 2008, and 16139 ind. (1.20 t) in 2009, with an average of 15209 ind. (1.19 t). Population abundance was estimated by applying Length-structured Virtual Population Analysis to length-frequency and annual yield data. The estimated annual-average abundance ofR.cylindricusin Zhuyangxi was 16361 ind./km (1.25 t/km) in 2007, 13922 ind./km (1.74 t/km) in 2008, and 26836 ind./km (1.93 t/km) in 2009, respectively, with an average of 19040 ind./km (1.64 t/km). The estimated growth-related population parameters ofR.cylindricusby length-frequency data in this study are consistent with previous studies based on age- length data, which indicates that the result is reliable. Considering the current high exploitation rate forR.cylindricus, we recommend a minimum catchable length of 193.5 mm to decrease its exploitation rate. We further suggest that management authorities establish a closed period for fishing and assist commercial fishermen with alternative career options to reduce harvesting pressure on the endemic species.

population parameters; growth; mortality; stock assessment;Rhinogobiocylindricus

國家自然科學(xué)基金(51109091，51310105036)

2014- 03- 10;

日期:2015- 04- 20

10.5846/stxb201403100406

*通訊作者Corresponding author.E-mail: duan@yfi.ac.cn

熊飛，劉紅艷，段辛斌，劉紹平，陳大慶.長江上游朱楊溪江段圓筒吻鮈種群參數(shù)和資源量.生態(tài)學(xué)報,2015,35(22):7320- 7327.

Xiong F, Liu H Y, Duan X B, Liu S P, Chen D Q.Population parameters and population abundance ofRhinogobiocylindricusin Zhuyangxi section of the upper Yangtze River.Acta Ecologica Sinica,2015,35(22):7320- 7327.