洪澤湖大銀魚和太湖新銀魚的生長、死亡參數(shù)及資源利用狀況

趙麗爽　程　飛　張　磊　王　欣　張勝宇　蔣功成　謝松光

(1. 中國科學(xué)院水生生物研究所中國科學(xué)院水生生物多樣性與保護(hù)重點(diǎn)實驗室, 武漢 430072; 2. 中國科學(xué)院大學(xué), 北京100049; 3. 中國科學(xué)院水生生物研究所淮安研究中心, 淮安 223002; 4. 江蘇省洪澤湖漁業(yè)管理委員會辦公室,淮安 223002; 5. 淮陰師范學(xué)院, 淮安 223300)

大銀魚(Protosalanx hyalocranius Abbott, 1901)和太湖新銀魚(Neosalanx taihuensis Chen, 1956)分別隸屬于鮭形目(Salmoniformes)銀魚科(Salangidae)的大銀魚屬和新銀魚屬, 主要分布在我國黃河、淮河和長江中下游及其附屬湖泊, 是我國重要的經(jīng)濟(jì)魚類[1]。20世紀(jì)70年代以來, 2種銀魚被引種移植到全國眾多湖泊和水庫[2]。大銀魚主要被移植到長江以北[3], 太湖新銀魚則多被移植到黃河以南[4]。2種銀魚的生命周期均為一年, 個體小、生長快、世代離散, 種群波動顯著[1,5]。

過度捕撈是影響魚類資源可持續(xù)利用的重要因素, 分為生長型過度捕撈和補(bǔ)充型過度捕撈[6]。由于銀魚世代離散, 補(bǔ)充型過度捕撈導(dǎo)致剩余產(chǎn)卵群體數(shù)量不足, 進(jìn)而使種群下一世代數(shù)量減少, 被認(rèn)為是其資源量下降的重要原因[7]。單位補(bǔ)充量產(chǎn)卵群體生物量(Spawning Biomass-recruitment,SBR)模型根據(jù)種群的生長、死亡參數(shù)和種群資源動態(tài), 估算捕撈后產(chǎn)卵群體的殘存量, 進(jìn)而評估種群的資源利用狀況[8,9]。當(dāng)SBR殘存量低于未開發(fā)狀態(tài)的20%時, 通常視為補(bǔ)充型過度捕撈, 可能導(dǎo)致種群崩潰[10,11]。平衡產(chǎn)量模型是漁業(yè)資源評估最常用的Beverton-Holt模型的簡化形式[6,12], 常用于分析開捕時間和捕撈強(qiáng)度對漁獲產(chǎn)量的影響, 并據(jù)此制定合理的漁業(yè)管理方案[13,14]。

洪澤湖(33°06—40′N, 118°10—52′E)是中國第四大淡水湖, 水域面積約2069 km2。大銀魚和太湖新銀魚是洪澤湖最主要的經(jīng)濟(jì)魚類之一, 2010年以來年漁產(chǎn)量最高曾達(dá)106kg[15], 但年際波動較大, 且近年來呈現(xiàn)下降趨勢(江蘇省洪澤湖漁業(yè)管理委員會辦公室, 未發(fā)表數(shù)據(jù))。洪澤湖銀魚捕撈采用銀魚拖網(wǎng)進(jìn)行, 每年1月1日至8月8日為禁捕期, 其中5月21日左右允許捕撈1周。據(jù)漁民反映, 在每年12月下旬捕撈期的最后階段2種銀魚均已很難捕撈、漁獲量極低。本研究對洪澤湖大銀魚和太湖新銀魚的周年生長和資源量動態(tài)進(jìn)行了分析, 評估種群捕撈利用狀況, 進(jìn)而優(yōu)化漁業(yè)管理方式, 以期為洪澤湖銀魚的資源可持續(xù)利用管理提供參考。

1　材料與方法

1.1　采樣水域與采樣方法

2015年7月至2016年6月每月中旬, 對洪澤湖大銀魚和太湖新銀魚進(jìn)行周年采樣（圖1）。采樣使用漁民用的銀魚拖網(wǎng)進(jìn)行, 拖網(wǎng)網(wǎng)口寬11.6 m, 高2.3 m; 網(wǎng)體長約25 m, 由4個網(wǎng)片組成, 側(cè)面網(wǎng)片網(wǎng)目為1.2 cm×1.2 cm, 上下網(wǎng)片網(wǎng)目為1.5 cm×1.5 cm; 末端網(wǎng)袋長1.3 m, 網(wǎng)目為0.3 cm×0.3 cm。采樣時在機(jī)動漁船左右兩側(cè)各設(shè)置一個拖網(wǎng), 拖網(wǎng)上沿恰好沒入水中, 漁船在固定起點(diǎn)(A, 33°17′19″N,118°47′17″E)和終點(diǎn)(B, 33°19′58″N, 118°46′50″E)間拖行5 km, 拖速為4.0 km/h。采樣在上午9：00—12：00間進(jìn)行。樣品采集后立即置于冰水保溫箱中, 帶回實驗室計數(shù), 進(jìn)行生物學(xué)測定。每次采樣2種銀魚生物學(xué)測定數(shù)量均在100尾以上, 當(dāng)樣品數(shù)量不足時, 增加拖行時間。

圖1　洪澤湖, 示大銀魚和太湖新銀魚采集時拖網(wǎng)拖行區(qū)間起點(diǎn)(A)和終點(diǎn)(B)Fig. 1　Hongze Lake, showing the trawling area (from A to B) for P. hyalocranius and N. taihuensis collection

1.2　數(shù)據(jù)獲取與分析

體長-體重關(guān)系和生長率在實驗室對大銀魚和太湖新銀魚的漁獲物進(jìn)行計數(shù)和稱重。測定個體體長(L, 0.1 mm), 稱體重(W, 0.01 g), 解剖鑒定性別, 繁殖季節(jié)肉眼確定雌性性腺發(fā)育分期[16,17]。體長-體重關(guān)系采用W=a×Lb進(jìn)行擬合。采用單因子方差分析(ANOVA)檢驗各月2種銀魚體長和體重差異顯著性, 并采用Bonferroni方法進(jìn)行多重比較。分別求出2種銀魚各月的平均體長和體重, 計算相鄰月份間體長和體重瞬時生長率(gl, gw)[18]：

式中, Ln和Wn分別為某月(n)大銀魚或太湖新銀魚的平均體長和體重; Ln+1和Wn+1則為下一月份的平均體長和體重。

最適生長方程洪澤湖大銀魚的繁殖季節(jié)為12至次年2月, 盛期為1月中旬; 太湖新銀魚春群的繁殖季節(jié)為3—4月, 盛期為3月(未發(fā)表數(shù)據(jù))。據(jù)此, 我們將1月25日和3月25日分別設(shè)為大銀魚和太湖新銀魚的出生時間, 即0齡, 并推算2種銀魚各月采樣時間對應(yīng)的年齡。采用Logistic生長方程(Logistic growth function, Logistic GF)、Gompertz生長方程(Gompertz growth function, Gompertz GF)和von Bertalanffy生長方程(Von Bertalanffy Growth Function, VBGF)擬合2種銀魚的年齡與體長生長關(guān)系, 使用最小二乘法估算各生長方程的參數(shù), 依據(jù)赤井信息準(zhǔn)則(Akaike’s information criterion,AIC)選擇最適生長方程[19]。3種生長方程分別為：

式中, t表示樣本年齡; t0為理論生長起點(diǎn), 即體長Lt=0時的起始年齡; L∞為漸近體長; K為生長系數(shù)。生長方程擬合采用R-3.3.1軟件進(jìn)行[20]。

種群死亡系數(shù)定義單位捕撈努力量漁獲量(Catch per Unit Effort, CPUE)為拖網(wǎng)拖行1 km2水域面積的漁獲物個體數(shù)量和生物量。使用CPUE逐月數(shù)量變動數(shù)據(jù), 結(jié)合洪澤湖當(dāng)前的銀魚捕撈管理方式(8月9日—12月31日為捕撈期, 同時5月21—28日可對銀魚資源進(jìn)行為期一周的捕撈), 求解以下方程組, 計算捕撈死亡系數(shù)F和自然死亡系數(shù)M[21]：

式中, Z為總死亡系數(shù); t2、t1分別為不同月份樣本的年齡差值, 單位為year;分別表示相應(yīng)年齡時對應(yīng)的漁獲物資源數(shù)量。在本研究中假定銀魚種群在一年內(nèi)的F和M恒定, 以年平均狀態(tài)表示,故用方程組求解的Z和M均值來表示種群的捕撈死亡系數(shù)和自然死亡系數(shù); F為Z與M的均值之差。

模型構(gòu)建使用平衡產(chǎn)量模型(The equilibrium output model)估算種群單位補(bǔ)充量漁產(chǎn)量隨開捕年齡和捕撈強(qiáng)度的變化[6,13]：

式中, Y為以重量表示的單位補(bǔ)充量的漁產(chǎn)量; N0為年齡t0時的假設(shè)個體數(shù), 假定N0=1000尾; r為初次開捕年齡與理論生長起點(diǎn)年齡之差, 即r=tc-t0; W∞為漸近體重

使用SBR模型對洪澤湖當(dāng)前銀魚捕撈管理方式和依據(jù)平衡產(chǎn)量模型制定的新的管理方式進(jìn)行評估。假定2種方案下銀魚的捕撈死亡系數(shù)和自然死亡系數(shù)不變, 將模型中的時間步長設(shè)置為1周, 即0.0192 year[10], SBR模型公式為[8]：式中, tmax表示樣本的最大年齡; St和Gt分別表示刀刃式選擇假設(shè)下的網(wǎng)具和性成熟選擇性, 當(dāng)t＜50%性成熟體長時, Gt=0, 反之, Gt=1; 當(dāng)t＜初次開捕年齡時, St=0, 反之, St=1[22]; At對應(yīng)于t齡時是否處于漁業(yè)開捕狀態(tài), 若是, 則At=1, 否則At=0。

對于SBR模型中使用的50%性成熟體長, 針對2種銀魚1年生且1次性成熟的特點(diǎn), 以繁殖期性成熟的最小10尾雌性個體的體長均值Lm表示[23], 將其代入最適生長方程中, 計算對應(yīng)的性成熟年齡Tm。

2　結(jié)果

2.1　種群結(jié)構(gòu)與生長

共測定了洪澤湖1250尾大銀魚和1818尾太湖新銀魚的體長和體重。大銀魚體長和體重分別為29.0—182.6 mm和0.10—34.79 g; 太湖新銀魚體長和體重分別為20.4—82.7 mm和0.04—3.40 g。銀魚的體長-體重關(guān)系方程分別為：

太湖新銀魚：(n=1818;r2=0.951;P＜0.05)

大銀魚一個世代在漁獲物中出現(xiàn)的時間從4月到翌年2月, 9月以后大銀魚個體間體長差異逐漸增大、體長范圍逐漸增加, 12月后大個體在種群中所占比例下降(圖2a)。世代周期中存在4—6月和8—11月2個快速生長階段。這2個階段相鄰月份的體長和體重均差異顯著(P＜0.05)。體長和體重均在11月達(dá)到最大值。體長和體重瞬時生長率總體均呈下降趨勢, 12月后為負(fù)值(圖3、圖4)。

太湖新銀魚一個世代在漁獲物中出現(xiàn)的時間從5月到翌年4月, 為春季繁殖群體(春群)的后代;9月也采集到少量性成熟個體(秋群), 其后代于10月開始出現(xiàn)在漁獲物中, 但數(shù)量較少(圖2b)。本研究僅考慮春群。太湖新銀魚整個世代周期中持續(xù)生長, 無明顯的階段性。5—12月間相鄰月份間體長差異顯著(P＜0.05), 但體重差異不顯著(P＞0.05)。太湖新銀魚體長和體重瞬時生長率變化較大銀魚平穩(wěn), 總體呈下降趨勢, 12月后近乎于零(圖3、圖4)。

2.2　最適生長方程與性成熟大小

在3種生長方程中, VBGF對大銀魚生長擬合的AIC值最小, Logistic GF對太湖新銀魚生長擬合的AIC值最小(表1), 分別為2種魚類的最適生長方程(圖5)。大銀魚和太湖新銀魚卵巢發(fā)育至Ⅳ期的10尾最小個體的體長均值Lm分別為99.9和61.1 mm,代入對應(yīng)的生長方程, 得出大銀魚和太湖新銀魚的Tm分別為0.528 year和0.564 year。

2.3　CPUE逐月變動與種群死亡系數(shù)

單位水域面積(1 km2)中大銀魚和太湖新銀魚漁獲物數(shù)量均在5月最高, 分別為31809和78098 ind./km2。經(jīng)歷了5月中旬為期1周的捕撈后, 6月2種銀魚的資源數(shù)量均明顯下降; 6月至8月間下降相對平緩; 8月捕撈開始后下降明顯, 到12月時漁獲物數(shù)量均降至500 ind./km2以下。單位水域面積(1 km2)大銀魚和太湖新銀魚漁獲物生物量的最高值分別出現(xiàn)在6月和8月, 為73.07和49.33 kg/km2。大銀魚生物量在8—9月之間下降最快，太湖新銀魚生物量在9—10月間下降最快。

根據(jù)洪澤湖現(xiàn)行的捕撈管理方式, 5—6月間和8月以后, 2種銀魚的CPUE下降受捕撈和自然死亡的共同影響; 6—8月間, CPUE的下降只受自然死亡的影響。據(jù)此求出大銀魚F和M分別為8.583/year和3.292/year，太湖新銀魚F和M分別為7.006/year和1.146/year。

圖2　洪澤湖大銀魚(a)和太湖新銀魚(b)的體長頻率分布逐月變化Fig. 2　Monthly variations in body length frequency distribution for P. hyalocranius (a) and N. taihuensis (b) in Hongze Lake

圖3　洪澤湖大銀魚和太湖新銀魚(不包括秋群)的體長與其瞬時生長率的逐月變化(誤差線代表標(biāo)準(zhǔn)差; *表示兩個相鄰月份間體長差異顯著, P＜0.05)Fig. 3　Monthly fluctuations in body length and its instantaneous growth rate for P. hyalocranius and N. taihuensis in Hongze Lake(Error bars represent ± standard deviation; * indicates significant difference in body length between two adjacent months, P＜0.05)

圖4　洪澤湖大銀魚和太湖新銀魚(不包括秋群)的體重與其瞬時生長率的逐月變化(誤差線代表標(biāo)準(zhǔn)差; *表示兩個相鄰月份間體重差異顯著, P＜0.05)Fig. 4　Monthly fluctuations in body weight and its instantaneous growth rate for P. hyalocranius and N. taihuensis in Hongze Lake(Error bars represent ± standard deviation; * indicates significant difference in body weight between two adjacent months, P＜0.05)

表1　三種生長方程對大銀魚和太湖新銀魚生長擬合的方程參數(shù)Tab. 1　Estimated parameters of three different growth functions for P. hyalocranius and N. taihuensis in Hongze Lake

圖5　洪澤湖大銀魚(a)和太湖新銀魚(b)的最適生長方程曲線Fig. 5　The most fitted growth functions for P. hyalocranius (a)and N. taihuensis (b) in Hongze Lake

2.4　平衡產(chǎn)量模型分析

使用以上獲得的生長、繁殖和死亡等相關(guān)參數(shù)(表2), 分別構(gòu)建了大銀魚和太湖新銀魚單位補(bǔ)充量漁產(chǎn)量隨開捕年齡(圖6a)和捕撈死亡系數(shù)(圖6b)變化的曲線。結(jié)果顯示, 大銀魚和太湖新銀魚的單位補(bǔ)充量漁產(chǎn)量均隨開捕年齡的增大呈先升后降的倒U型曲線關(guān)系。當(dāng)大銀魚開捕年齡tc=0.525 year時, 能夠取得相應(yīng)最大的單位補(bǔ)充量漁產(chǎn)量Ymax=1.078 kg, 開捕年齡過早或過晚均不能獲得足夠的漁產(chǎn)量。就洪澤湖當(dāng)前漁業(yè)管理模式來說,大銀魚的開捕年齡約為0.538 year, 略大于tc, 此時的單位補(bǔ)充量漁產(chǎn)量Y=1.077 kg, 接近Ymax。當(dāng)太湖新銀魚開捕年齡tc=0.525 year時, 能夠取得相應(yīng)最大的單位補(bǔ)充量漁產(chǎn)量Ymax=0.684 kg。就當(dāng)前漁業(yè)管理模式來說, 太湖新銀魚的開捕年齡約為0.365 year, 此時的單位補(bǔ)充量漁產(chǎn)量Y=0.607 kg＜Ymax。將開捕時間向后推移20d (即8月29日), 即大銀魚和太湖新銀魚的開捕年齡分別為0.593 year和0.420 year時, 能夠獲得最大的總單位補(bǔ)充量漁產(chǎn)量1.708 kg。

表2　洪澤湖大銀魚和太湖新銀魚平衡產(chǎn)量模型和單位補(bǔ)充量產(chǎn)卵群體生物量模型參數(shù)Tab. 2　Growth and mortality parameters fitted into the equilibrium output model and SBR model for P. hyalocranius and N.taihuensis in Hongze Lake

大銀魚和太湖新銀魚的單位補(bǔ)充量漁產(chǎn)量均先隨捕撈死亡系數(shù)的增加而快速增高, 當(dāng)捕撈死亡系數(shù)增大到一定程度后, 漁產(chǎn)量上升速度減慢, 之后不再升高甚至出現(xiàn)下降趨勢。大銀魚捕撈死亡系數(shù)F=24.38/year時, 獲得最大單位補(bǔ)充量漁產(chǎn)量Ymax=1.114 kg; 當(dāng)前洪澤湖大銀魚捕撈死亡系數(shù)為8.583/year, 相應(yīng)漁產(chǎn)量Y=1.077 kg, 接近Ymax且尚未開發(fā)過度。太湖新銀魚捕撈死亡系數(shù)F=4.72/year時, 獲得最大的單位補(bǔ)充量漁產(chǎn)量Ymax=0.616 kg;而當(dāng)前洪澤湖太湖新銀魚捕撈死亡系數(shù)為7.006/year,已大于最適捕撈強(qiáng)度, 且相應(yīng)的漁產(chǎn)量Y=0.607 kg,小于Ymax。

圖6　洪澤湖大銀魚和太湖新銀魚(春群)單位水域面積(1 km2)漁獲物數(shù)量(a)和生物量(b)的逐月變動Fig. 6　Monthly variations in number of individuals (a) and biomass (b) for P. hyalocranius and N. taihuensis harvested across an area of 1 km2 by trawl in Hongze Lake

2.5　基于SBR模型分析的資源狀況評價

使用以上獲得的生長、繁殖和死亡等相關(guān)參數(shù)(表2), 分別構(gòu)建了2種不同漁業(yè)管理方案下大銀魚和太湖新銀魚種群的SBR隨捕撈死亡系數(shù)變動的曲線：

方案1： 洪澤湖現(xiàn)行的漁業(yè)管理方式。8月9日—12月31日為捕撈期, 同時5月可對銀魚資源進(jìn)行為期1周的捕撈。

方案2： 假定對兩種銀魚的捕撈強(qiáng)度不發(fā)生變化(即F恒定), 將開捕時間向后推遲20d, 8月29日—12月31日為捕撈期, 同時取消5月為期1周的捕撈。

大銀魚在當(dāng)前漁業(yè)管理方案1下, 相對于未開發(fā)狀態(tài)(100%), SBR殘存量為20.234%, 即3.326 g,接近于補(bǔ)充型過度捕撈評價閾值20%, 種群有可能因產(chǎn)卵親體不足而導(dǎo)致資源衰退; 在方案2下, SBR殘存量相當(dāng)于未開發(fā)狀態(tài)的36.721%, 即6.036 g, 遠(yuǎn)大于補(bǔ)充型過度捕撈評價閾值20%, 大銀魚漁業(yè)資源和產(chǎn)卵群體得到了較好的保護(hù)(圖7)。在2種不同管理方案下, 太湖新銀魚SBR殘存量分別相當(dāng)于未開發(fā)狀態(tài)的7.503%和12.862%, 即0.972 g和1.666 g(圖8)。在優(yōu)化的管理方案2下, 太湖新銀魚的產(chǎn)卵親體資源雖然得到了一定程度的保護(hù), 但其漁業(yè)資源仍處于補(bǔ)充型過度捕撈水平。

3　討論

3.1　大銀魚和太湖新銀魚的階段性生長差異

圖7　平衡產(chǎn)量模型估算的洪澤湖大銀魚和太湖新銀魚單位補(bǔ)充量漁產(chǎn)量隨開捕年齡(a)和捕撈死亡系數(shù)(b)變化的曲線Fig. 7　Relationships of yield per-recruitment to the age at first capture (a) and fishing mortality coefficient (b) estimated from the equilibrium output mode for P. hyalocranius and N. taihuensis in Hongze Lake

圖8　兩種不同捕撈管理方案下洪澤湖大銀魚(a)和太湖新銀魚(b)單位補(bǔ)充量產(chǎn)卵群體生物量隨捕撈死亡系數(shù)F變化的曲線Fig. 8　Relationships of spawner biomass per-recruitment to fishing mortality coefficients for P. hyalocranius (a) and N.taihuensis (b) under two different fishery management schemes in Hongze Lake

魚類的生長具有階段性, 生活史的不同階段,個體的生理特征、食性以及環(huán)境因素等差異, 會引起生長速度的改變[24,25]。本研究發(fā)現(xiàn), 洪澤湖大銀魚存在2個快速生長階段, 分別在4—6月和8—11月間。魚類一般在早期生活史階段和性成熟前體長生長較快[26], 這一階段的快速生長可能是魚類維持種群數(shù)量的一種生態(tài)適應(yīng)方式, 有利于降低其被捕食的壓力、快速達(dá)到性成熟體長[27,28]。洪澤湖大銀魚繁殖期在12—2月(作者未發(fā)表數(shù)據(jù)), 4—6月間的第一個快速生長階段, 在個體性成熟之前; 同時,該時期水體溫度逐步上升, 初級生產(chǎn)力升高, 水體中浮游動物等大銀魚的主要餌料生物豐富, 是這一階段大銀魚快速生長的營養(yǎng)保障[29]。8—11月間第二階段快速生長可能主要源于大銀魚的食性轉(zhuǎn)化。唐富江發(fā)現(xiàn)興凱湖大銀魚的體長大于60 mm以后, 逐漸由攝食浮游動物轉(zhuǎn)為攝食魚蝦類[30]; 本文作者也發(fā)現(xiàn), 自8月起蝦類和魚類逐漸取代枝角類和橈足類成為洪澤湖大銀魚的主要食物(作者未發(fā)表數(shù)據(jù))。魚蝦類相對于浮游動物而言具有體積大、能量高、動物蛋白豐富等特點(diǎn), 更易于被魚類吸收和利用, 作為生長和繁殖的能量供應(yīng)[31]。因此,大銀魚在8—11月的快速生長很可能是因為有更多的能量來源投入到生長中。12月后, 大銀魚種群體長和體重均表現(xiàn)為負(fù)增長, 這可能與其進(jìn)入繁殖期后大個體較早的繁殖和死亡有關(guān)[32], 需進(jìn)一步研究確證。

太湖新銀魚個體相對較小, 其體長在5—10月呈快速持續(xù)生長狀態(tài), 無明顯階段性。太湖新銀魚終生攝食浮游動物, 與大銀魚不同, 沒有明顯的食性轉(zhuǎn)化[33]。太湖新銀魚這一食性特點(diǎn)可能是其生長過程沒有明顯階段性的主要原因。11月以后, 太湖新銀魚個體體長無明顯的增長趨勢。洪澤湖太湖新銀魚繁殖季節(jié)主要在3—4月, 11月之后, 繁殖季節(jié)前, 是性腺發(fā)育能量積累的主要階段[34]。因此,能量投入到性腺發(fā)育中, 可能是太湖新銀魚11月后體長停止增長的重要原因。

3.2　洪澤湖銀魚資源利用評估與管理

魚類種群可提供漁獲量的多少取決于種群數(shù)量、生長、存活、開捕年齡和捕撈強(qiáng)度等, 其中僅有捕撈強(qiáng)度(捕撈死亡系數(shù))和初次捕撈年齡(初次開捕體長)能夠被人為控制, 因此漁業(yè)管理中常通過調(diào)整這2個因素以獲取最大可持續(xù)漁產(chǎn)量[6]。依據(jù)平衡產(chǎn)量模型分析, 洪澤湖大銀魚和太湖新銀魚在當(dāng)前捕撈強(qiáng)度下的漁產(chǎn)量均已接近其資源群體所能提供的最大漁產(chǎn)量, 調(diào)整捕撈強(qiáng)度對進(jìn)一步提高漁產(chǎn)量的影響不大; 而調(diào)整2種銀魚的初次開捕年齡, 更有利于獲得較高的總漁產(chǎn)量。

在洪澤湖當(dāng)前的捕撈管理方式下, 大銀魚和太湖新銀魚的SBR殘存量分別為未開發(fā)狀態(tài)的20.23%和7.50%, 表明對2種銀魚的捕撈均處于補(bǔ)充型過度捕撈水平, 且太湖新銀魚捕撈后的SBR殘余量已遠(yuǎn)低于種群崩潰的警戒線。由此看來, 較大的捕撈強(qiáng)度和補(bǔ)充型過度捕撈, 可能是導(dǎo)致近年來洪澤湖大銀魚和太湖新銀魚資源波動明顯的重要原因。因此, 很有必要對洪澤湖當(dāng)前銀魚資源捕撈管理方式進(jìn)行優(yōu)化, 避免過度捕撈的發(fā)生。

本研究構(gòu)建平衡產(chǎn)量模型時發(fā)現(xiàn), 將洪澤湖銀魚開捕時間推遲20d可獲得大銀魚和太湖新銀魚的最大總漁產(chǎn)量。在漁業(yè)管理中, 這種開捕時間的推遲也將影響最終捕撈后的產(chǎn)卵親體生物量[8,9]。對洪澤湖的銀魚捕撈管理方式進(jìn)行優(yōu)化后, 大銀魚的SBR殘存量提高至未開發(fā)狀態(tài)的36.72%, 既滿足了獲得最高漁產(chǎn)量的需求, 又不會帶來太大的漁業(yè)風(fēng)險; 優(yōu)化管理后太湖新銀魚的SBR殘存量亦提高至未開發(fā)狀態(tài)的12.86%。可見, 優(yōu)化的漁業(yè)管理方式實現(xiàn)了獲取最高漁產(chǎn)量和保護(hù)產(chǎn)卵親體資源上的“雙贏”。同時也應(yīng)注意到, 優(yōu)化的漁業(yè)管理方案對大銀魚產(chǎn)卵群體的保護(hù)作用更為顯著; 在未改變當(dāng)前捕撈強(qiáng)度的情況下, 盡管太湖新銀魚在捕撈后的產(chǎn)卵親體資源量也得到了一定程度的提升, 但依然低于補(bǔ)充型過度捕撈評價閾值20%, 故還需減小對太湖新銀魚的捕撈強(qiáng)度。因此, 針對洪澤湖太湖新銀魚的漁業(yè)資源管理, 在延遲開捕時間的基礎(chǔ)上,建議還可以通過限制漁獲物的起捕體長、增大捕撈網(wǎng)具的網(wǎng)目規(guī)格等方式控制捕撈強(qiáng)度, 同時可采取人工放流等措施提高其下一代的資源補(bǔ)充能力,以避免種群資源的進(jìn)一步衰退。

參考文獻(xiàn)：

[1]Xie Y H, Xie H. Classification, distribution and population ecology of Salangidae fishes [J]. Chinese Journal of Fisheries, 1997, 10(2)： 11—19 [解玉浩, 解涵. 銀魚科魚類的分類分布和種群生態(tài). 水產(chǎn)學(xué)雜志, 1997, 10(2)：11—19]

[2]Qin W. Discussion on the present situation and development for icefishes transplanting to lakes and reservoirs in China [J]. Fisheries Science and Technology Information,1994, 21(4)： 156—158 [秦偉. 我國湖泊水庫銀魚移植現(xiàn)狀及漁業(yè)發(fā)展的探討. 水產(chǎn)科技情報, 1994, 21(4)：156—158]

[3]Ge Z Q, Zhu J, Shen Q Z, et al. A study on status and prospective of Protosalanx hyalocranius (Abbott) transplanting technique [J]. Modern Fisheries Information,2001, 16(7)： 13—16 [戈志強(qiáng), 朱江, 沈其璋, 等. 關(guān)于大銀魚 Protosalanx hyalocranius (Abbott) 移植現(xiàn)狀及其發(fā)展前景的探討. 現(xiàn)代漁業(yè)信息, 2001, 16(7)： 13—16]

[4]Gu L W. The status of Taihu Lake icefish transplant In China [J]. Modern Fisheries Information, 1998, 13(9)：23—25 [顧良偉. 我國太湖銀魚移植現(xiàn)狀. 現(xiàn)代漁業(yè)信息, 1998, 13(9)： 23—25]

[5]Qin J H, Xie P, Xu J. The growth features of Neosalanx taihuensis in Lake Fuxian and its fisheries utilization [J].Acta Hydrobiologica Sinica, 2009, 33(2)： 352—354 [覃劍暉, 謝平, 徐軍. 撫仙湖太湖新銀魚的生長特性及其漁業(yè)利用. 水生生物學(xué)報, 2009, 33(2)： 352—354]

[6]Zhan B Y. Fish Stock Assessment [M]. China Agriculture Press. 1995, 141—323 [詹秉義. 漁業(yè)資源評估. 中國農(nóng)業(yè)出版社. 1995, 141—323]

[7]Huang Z L, Chang J B. Can icefish (Salangidae) production be predicted [J]? Acta Ecologica Sinica, 2001, 21(1)：86—93 [黃真理, 常劍波. 銀魚的產(chǎn)量能預(yù)報嗎?. 生態(tài)學(xué)報, 2001, 21(1)： 86—93]

[8]Booth A J, Buxton C D. Management of the panga Pterogymnus laniarius (Pisces： Sparidae) on the Agulhas Bank,South Africa using per-recruit models [J]. Fisheries Research, 1997, 32(1)： 1—11

[9]Govender A, Al-Oufi H, McIlwain J L, et al. A per-recruit assessment of the kingfish (Scomberomorus commerson) resource of Oman with an evaluation of the effectiveness of some management regulations [J]. Fisheries Research, 2006, 77(2)： 239—247

[10]Caddy J. A Short Review of Precautionary Reference Points and Some Proposals for Their Use in Data-poor Situations [M]. Food and Agriculture Origination. 1998,8—22

[11]Feng B, Hou G, Wang X F. Research on length-based yield per recruit and spawning-biomass per recruit of Trichiurus lepturus in Beibu Gulf [J]. Journal of Tropical Organisms, 2011, 1(4)： 362—366 [馮波, 侯剛, 王學(xué)鋒. 基于體長的北部灣帶魚單位補(bǔ)充量親體量研究. 熱帶生物學(xué)報, 2011, 1(4)： 362—366]

[12]Ricker W E. Computation and Interpretation of Biological Statistics of Fish Populations [M]. Dept. of Fisheries and Oceans, Canada. 1987, 175—178

[13]Xie C X, Wang M X, Jin H. Study on the growth characteristics of Culter erythopterus in Baoan Lake [J]. Resources Environment and Fishery Ecological Management of Macrophytic Lakes, 1995, 281： 266—271 [謝從新, 王明學(xué), 金暉. 保安湖紅鰭原鲌生長特性的研究. 草型湖泊資源, 環(huán)境與漁業(yè)生態(tài)學(xué)管理 (一), 1995, 281：266—271]

[14]Lin L S, Cheng J H, Ling J Z, et al. First capture sizes of major commercial fishes in the East China Sea Region[J]. Journal of Fishery Sciences of China, 2006, 13(2)：250—256 [林龍山, 程家驊, 凌建忠, 等. 東海區(qū)主要經(jīng)濟(jì)魚類開捕規(guī)格的初步研究. 中國水產(chǎn)科學(xué), 2006,13(2)： 250—256]

[15]Lin M L, Zhang T L, Ye S W, et al. Status of fish resources, historical variation and fisheries management strategies in Hongze Lake [J]. Acta Hydrobiologica Sinica, 2013, 37(6)： 1118—1127 [林明利, 張?zhí)昧? 葉少文,等. 洪澤湖魚類資源現(xiàn)狀、歷史變動和漁業(yè)管理策略.水生生物學(xué)報, 2013, 37(6)： 1118—1127]

[16]Sun G Y. On the maturity of the female gonads and the type of spawning of the icefish [J]. Journal of Fisheries of China, 1985, 9(4)： 363—368 [孫幗英. 大銀魚卵巢的成熟期和產(chǎn)卵類型. 水產(chǎn)學(xué)報, 1985, 9(4)： 363—368]

[17]Wu L. Comparative study on life history strategy of introduced Neosalanx taihuensis at different latitudes [D].Thesis for Doctor of Science, Institute of Hydrobiology,Chinese Academy of Sciences, Wuhan. 2012 [吳朗. 太湖新銀魚不同緯度移殖種群生活史對策比較研究. 博士學(xué)位論文, 中國科學(xué)院水生生物研究所, 武漢. 2012]

[18]Yin M C. Ecology of Fishes [M]. China Agriculture Press. 1995, 53 [殷名稱. 魚類生態(tài)學(xué). 中國農(nóng)業(yè)出版社.1995, 53]

[19]Katsanevakis S. Modelling fish growth： model selection,multi-model inference and model selection uncertainty[J]. Fisheries Research, 2006, 81(2)： 229—235

[20]Guan W J. The Application of R in Marine Fishery [M].China Ocean Press. 2015, 164—173 [官文江. R語言在海洋漁業(yè)中的應(yīng)用. 海洋出版社. 2015, 164—173]

[21]Jensen A L. Comparison of catch-curve methods for estimation of mortality [J]. Transactions of the American Fisheries Society, 1985, 114(5)： 743—747

[22]Taylor B M, McIlwain J L. Beyond abundance and biomass： effects of marine protected areas on the demography of a highly exploited reef fish [J]. Marine Ecology Progress Series, 2010, 411(12)： 243—258

[23]Sparre P, Venema S C. Introduction to Tropical Fish Stock Assessment-Part 1： Manual [M]. Food and Agriculture Origination. 1998, 1—283

[24]Marchant T A, Peter R E. Seasonal variations in body growth rates and circulating levels of growth hormone in the goldfish, Carassius auratus [J]. Journal of Experimental Zoology Part A: Ecological Genetics and Physiology, 1986, 237(2)： 231—239

[25]Bacon P J, Gurney W S C, Jones W, et al. Seasonal growth patterns of wild juvenile fish： partitioning variation among explanatory variables based on individual growth trajectories of Atlantic salmon (Salmo salar) parr[J]. Journal of Animal Ecology, 2005, 74(1)： 1—11

[26]Fuiman L A, Werner R G. Fishery Science： the Unique Contributions of Early Life Stages [M]. Blackwell Science. 2002, 33—63

[27]Dannewit J, Petersson E. Association between growth,body condition and anti-predator behavior in maturing and immature brown trout parr [J]. Journal of Fish Biology, 2001, 59(4)： 1081—1091

[28]Stallings C D, Coleman F C, Koenig C C, et al. Energy allocation in juveniles of a warm-temperate reef fish [J].Environmental Biology of Fishes, 2010, 88(4)： 389—398

[29]Meekan M G, Carleton J H, McKinnon A D, et al. What determines the growth of tropical reef fish larvae in the plankton： food or temperature [J]? Marine Ecology Progress Series, 2003, 256(1)： 193—204

[30]Tang F J, Liu W, Wang J L, et al. Diet composition and transition of clearhead icefish (Protosalanx hyalocranius)in Lake Xingkai [J]. Zoological Research, 2013, 34(5)：493—498 [唐富江, 劉偉, 王繼隆, 等. 興凱湖大銀魚食物組成與食性轉(zhuǎn)化. 動物學(xué)研究, 2013, 34(5)：493—498]

[31]Jauncey K. The effects of varying dietary protein level on the growth, food conversion, protein utilization and body composition of juvenile tilapias (Sarotherodon mossambicus) [J]. Aquaculture, 1982, 27(1)： 43—54

[32]Zhang K X, Zhuang D D, Zhang L, et al. Biology study on the Protosalanx hyalocranius and its propagation in Hongze Lake [J]. Journal of Fisheries of China, 1981,5(1)： 29—39 [張開翔, 莊大棟, 張立, 等. 洪澤湖所產(chǎn)大銀魚生物學(xué)及其增值的研究. 水產(chǎn)學(xué)報, 1981, 5(1)：29—39]

[33]Liu Z W, Zhu S Q. Food and feeding behavior of ice-fish(Neosalanx taihuensis Chen) in Dianchi Lake [J]. Acta Zoologica Sinica, 1994, 40(3)： 253—261 [劉正文, 朱松泉. 滇池產(chǎn)太湖新銀魚食性與攝食行為的初步研究. 動物學(xué)報, 1994, 40(3)： 253—261]

[34]McBride R S, Somarakis S, Fitzhugh G R, et al. Energy acquisition and allocation to egg production in relation to fish reproductive strategies [J]. Fish and Fisheries, 2015,16(1)： 23—57