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

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

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

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

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

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

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

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

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

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

1 材料與方法

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

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

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

1.2 種群參數(shù)估算

體長(zhǎng)體重關(guān)系依據(jù)Keys公式W=aLb進(jìn)行擬合，用t檢驗(yàn)法檢驗(yàn)冪指數(shù)b值與3的偏離程度，判斷其生長(zhǎng)是否為勻速生長(zhǎng)類型[15-16]。利用Von Bertalanffy生長(zhǎng)方程描述生長(zhǎng)規(guī)律，漸近體長(zhǎng)(L∞)、生長(zhǎng)系數(shù)(k)和生長(zhǎng)特征指數(shù)(φ)根據(jù)魚類體長(zhǎng)頻率的時(shí)間系列數(shù)據(jù)借助FiSATⅡ漁業(yè)評(píng)估軟件中的體長(zhǎng)頻率分析法(ELEFANⅠ)模塊進(jìn)行估算，總死亡系數(shù)(Z)、開捕體長(zhǎng)(L50)由體長(zhǎng)轉(zhuǎn)漁獲物曲線(Converted-length catch curve)模塊進(jìn)行估算[3]。理論生長(zhǎng)起點(diǎn)年齡(t0)、自然死亡系數(shù)(M)根據(jù)Pauly經(jīng)驗(yàn)公式獲得[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動(dòng)態(tài)綜合模型的相對(duì)單位補(bǔ)充量產(chǎn)量(Y′/R)和相對(duì)單位補(bǔ)充量生物量(B′/R)曲線評(píng)價(jià)圓筒吻鮈的資源利用情況，Emax表示能獲得最大產(chǎn)量時(shí)的開發(fā)率，E10表示資源群體邊際產(chǎn)量減少到理論原始資源邊際產(chǎn)量1/10時(shí)的開發(fā)率，被認(rèn)為是“最適”開發(fā)率，E50表示資源量下降到原始水平50%時(shí)的開發(fā)率[3]。

1.3 漁獲量計(jì)算

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

Y=∑(Yi×Pi)

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

Yi=xi×fi×ti

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

1.4 資源量估算

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

2 結(jié)果與分析

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

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

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

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

2.2 生長(zhǎng)參數(shù)

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

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

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

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

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

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

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

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

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

圖6 圓筒吻鮈的相對(duì)單位補(bǔ)充量產(chǎn)量(Y′/R)和相對(duì)單位補(bǔ)充生物量(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模塊估算圓筒吻鮈的資源量，通過迭代法計(jì)算出2007—2009年其最大體長(zhǎng)組的捕撈死亡系數(shù)Ft分別為0.1194、0.5030和0.1818。結(jié)果表明(圖7)，個(gè)體體長(zhǎng)在160mm以下時(shí)，資源群體的損失以自然死亡為主，而在160mm以上時(shí)，資源群體的損失以捕撈死亡為主。2007—2009年朱楊溪江段圓筒吻鮈年資源量分別為245411、208836和402537 尾，年均值為285594 尾；以生物量計(jì)，年資源量分別為18.80、26.09和28.94t，年均值為24.61t。以單位長(zhǎng)度河段來表示，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年圓筒吻鮈的體長(zhǎng)結(jié)構(gòu)實(shí)際種群分析Fig.7 Length structured VPA for Rhinogobio cylindricus in 2007—2009

3 討論

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

根據(jù)體長(zhǎng)頻率數(shù)據(jù)進(jìn)行魚類種群參數(shù)和資源量估算，可以避免年齡鑒定的主觀誤差而影響參數(shù)評(píng)估的準(zhǔn)確性，特別適合年齡鑒定比較困難的評(píng)估對(duì)象，在漁業(yè)資源評(píng)價(jià)中取得了較好的應(yīng)用效果[20- 22]。ELEFAN技術(shù)估算魚類生長(zhǎng)參數(shù)主要根據(jù)魚類體長(zhǎng)頻率的時(shí)間系列數(shù)據(jù)進(jìn)行體長(zhǎng)頻率樣品重構(gòu)，將體長(zhǎng)頻率形成的波峰和波谷放大，以擬合最佳生長(zhǎng)曲線，從而估算出合適的生長(zhǎng)參數(shù)，該方法對(duì)采樣頻率和樣品容量要求較高[23]。本研究由于春季禁漁和魚類種群數(shù)量偏少等原因，采樣頻率和樣品容量都受到了一定限制，但圓筒吻鮈體長(zhǎng)數(shù)據(jù)(69.0—268.0mm)基本包括了從小到大的各個(gè)體長(zhǎng)段的個(gè)體，具有一定的代表性。以往對(duì)圓筒吻鮈生長(zhǎng)參數(shù)的估算主要依據(jù)體長(zhǎng)-年齡數(shù)據(jù)，估算結(jié)果表明，長(zhǎng)江上游圓筒吻鮈的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)確性受評(píng)估模型的影響較大，考慮到本研究中已由ELEFAN技術(shù)估算出了L∞和k值，直接采用了比較通用的Pauly經(jīng)驗(yàn)公式進(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

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

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

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

由魚類的死亡系數(shù)可以計(jì)算出魚類的開發(fā)率E，該參數(shù)反映了種群資源的利用程度。本研究表明，圓筒吻鮈現(xiàn)有開發(fā)率為0.70，與熊星等的研究結(jié)果一致[13]。一般認(rèn)為魚類的最適開發(fā)率為0.5[27]，若按此標(biāo)準(zhǔn)，圓筒吻鮈資源處于過度利用狀態(tài)。在目前的開捕體長(zhǎng)(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)以下。

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

[3] 詹秉義. 漁業(yè)資源評(píng)估. 北京: 中國(guó)農(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] 陳作志, 邱永松, 黃梓榮. 南海北部白姑魚生長(zhǎng)和死亡參數(shù)的估算. 應(yīng)用生態(tài)學(xué)報(bào), 2005, 16(4): 712- 716.

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

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

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

[12] 王美榮, 楊少榮, 劉飛, 黎明政, 但勝國(guó), 劉煥章. 長(zhǎng)江上游圓筒吻鮈年齡與生長(zhǎng)的研究. 水生生物學(xué)報(bào), 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] 田輝伍, 段辛斌, 熊星, 羅宏偉, 劉紹平, 陳大慶. 長(zhǎng)江上游長(zhǎng)薄鰍生長(zhǎng)和種群參數(shù)的估算. 長(zhǎng)江流域資源與環(huán)境, 2013, 22(10): 1305- 1312.

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

[20] 陳國(guó)寶, 李永振, 陳丕茂, 舒黎明. 魚類最佳體長(zhǎng)頻率分析組距研究. 中國(guó)水產(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è)體長(zhǎng)分析方法ELEFAN和SLCA估算魚類生長(zhǎng)參數(shù)的研究. 海洋湖沼通報(bào), 2007, (3): 81- 87.

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

[25] 周永東, 徐漢祥, 潘國(guó)良, 賀舟挺. 東海區(qū)鯧魚、小黃魚資源量及其持續(xù)漁獲量的估算. 浙江海洋學(xué)院學(xué)報(bào): 自然科學(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] 李修峰, 黃道明, 謝文星, 謝山, 常秀嶺, 楊漢運(yùn), 張友謙, 何家慶. 漢江中游產(chǎn)漂流性卵魚類產(chǎn)卵場(chǎng)的現(xiàn)狀. 大連水產(chǎn)學(xué)院學(xué)報(bào), 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

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

2014- 03- 10;

日期:2015- 04- 20

10.5846/stxb201403100406

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

熊飛，劉紅艷，段辛斌，劉紹平，陳大慶.長(zhǎng)江上游朱楊溪江段圓筒吻鮈種群參數(shù)和資源量.生態(tài)學(xué)報(bào),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.