Assessment of Pennahia macrocephalus stock in Beibu Gulf by length Bayesian biomass estimation method
-
摘要: 大头白姑鱼 (Pennahia macrocephalus) 是北部湾主要经济鱼类种群之一，研究其资源状况对北部湾渔业资源的可持续利用具有重要意义。该研究基于2006—2014年北部湾大头白姑鱼生物学数据，通过长度贝叶斯生物量估算法 (Length-based Bayesian biomass estimation method, LBB) 估算大头白姑鱼的渐进体长 (L∞)、相对自然死亡率 (M/k)、相对捕捞死亡率 (F/k)、相对总死亡率 (Z/k)、最适开捕体长 (Lc_opt)、开发率 (E) 和相对生物量开发指标 (B/B0)。结果表明，2006—2014年北部湾大头白姑鱼的平均L∞为265 mm，平均Lc_opt为157 mm，M/k、F/k和Z/k分别为1.39、3.58和4.96；Lc_opt呈上升趋势，而L∞呈下降趋势。目前大头白姑鱼E较高 (0.58)，B/B0 小于0.5，表明当前北部湾大头白姑鱼资源处于过度开发状态。Abstract: Pennahia macrocephalus, one of the major commercial fishes in the Beibu Gulf, is of great significance to study its stocks status. Based on the biological data of P. macrocephalus collected in the Beibu Gulf from 2006 to 2014, we estimated the asymptotic body length (L∞), relative natural mortality (M/k), relative fishing mortality (F/k), relative total mortality (Z/k), optimal body length-at-first-capture (Lc_opt), exploitation rate (E) and relative biomass index (B/B0) by the length-based Bayesian biomass estimation method (LBB). Results show that the average L∞ and Lc_opt were 265 mm and 157 mm, respectively, and the relative mortality rates of M/k, F/k and Z/k were 1.39, 3.58 and 4.96, respectively. The Lc_opt of P. macrocephalus showed an increasing trend while the L∞ showed a decreasing trend. At present, the E of P. macrocephalus is relatively high (0.58), indicating the resources have been in the state of overexploited. The B/B0 is less than 0.5, indicating the overexploitation of current stocks.
-
20世纪80年代以来，受海洋污染加重和过度捕捞的双重影响，近海渔业资源不断衰退，渔获物低龄化、小型化、低质化、营养级下降趋势明显，海洋生态系统逆向发展[1-2]。为保护海洋渔业资源，促进海洋渔业可持续发展，1995年我国正式实施海洋伏季休渔制度，实施区域为东、黄、渤海35°N以南至27°N以北海域[3]，1999年覆盖至南海海域[4]。近年来，伏季休渔制度也随着渔业资源情况的变化而进行调整。2017年初农业部颁布被称为“史上最严休渔”的新海洋伏季休渔制度，南海海域于5月1日至8月16日休渔，是制度实施以来调整最大的一次[5]。
目前关于海洋伏季休渔制度实施效果的定量研究多基于东海海域[6-9]，南海海域伏季休渔效果研究仅限部分海区与重点水域[6,9-12]，尚未见针对全海域伏季休渔效果评估的研究报道。本研究基于南海海洋捕捞网络数据和现场调研数据，综合不同作业类型渔船的日均产量产值、单位捕捞努力量渔获量(catch per unit effort，CPUE)、渔获率与渔获结构等指标，总体评价2015—2017年南海休渔制度实施效果，以期为进一步完善休渔制度、促进南海渔业资源可持续发展提供参考依据。
1. 材料与方法
1.1 数据来源
采用中国水产科学研究院南海水产研究所南海海洋捕捞生产结构调查资料，在南海三省区沿海各市主要渔港设立调查站点，信息员逐日调查登记站点各类渔船的渔获信息，包括船名、作业类型、渔获量、马力、航次产值、出海天数、作业渔场等。本文仅涉及2015—2017年休渔前后1个月南海主要作业类型(包括单拖、双拖、桁杆拖虾、围网、流刺网等)渔船的捕捞数据。调查站位见图1。
1.2 数据处理与统计分析
1.2.1 数据处理
日均产值产量计算公式为：
$${D_i} = \frac{{{Y_i}}}{{{T_i}}}$$ (1) $${C_i} = \frac{{{V_i}}}{{{T_i}}}$$ (2) 式中Di为某作业类型单艘渔船第i航次的日均产量(kg)，Ci为某作业类型单艘渔船i航次的日均产值(元)，Yi为该渔船i航次的总渔获量(kg)，Vi为该渔船i航次的总产值，Ti为该渔船的作业时间(d)。
CPUE表示特定时间内投入渔业捕捞生产工具的数量和强度[6]，计算公式为：
$$P = \overline {{P_i}} $$ (3) $${P_i} = \frac{{{Y_i}}}{{{E_i} \times {T_i}}}$$ (4) P为某作业类型渔船在一定时间段内的平均CPUE [kg·(kW·d)–1]，Pi为渔船第i航次的捕捞努力量，Ei为该船的主机功率(kW)。
按照统计习惯，选取7个拖网主要经济鱼种 [蓝子鱼(Siganus)、鲹类(Carangidae)[13]、二长棘鲷(Paerargyrops edita)、带鱼(Trichiurus haumela)、金线鱼(Nemipterus virgatus)、刺鲳(Psenopsis anomala)、鱿鱼(Loligo chinensis)] 作为研究对象。渔获率与渔获结构计算公式为：
$${P_a} = \frac{{\sum {{Y_a}} }}{{\sum {\left( {{E_i} \times {T_i}} \right)} }}$$ (5) $${R_a} = \frac{V}{Y}$$ (6) Pa表示a种鱼类的CPUE，Ya表示a种鱼类的渔获量 [kg·(kW·d)–1]，Ra为a种鱼类占渔获物的数量百分比[14]，V为某类型渔船一定时段该种鱼类的产量(kg)，Y为该类型渔船渔获总产量。
1.2.2 统计分析
捕捞生产结构数据。根据不同作业类型、时间分组，采用公式(1)、(2)和(4)计算各作业类型休渔前后1个月日均产量产值及CPUE。采用SPSS 20.0描述性统计计算不同作业类型单位渔船日均产值产量、各年休渔前后CPUE的“平均值±标准差(
$\overline X \pm {\rm SD}$ )”，结合渔业统计年鉴南海海域相应作业类型渔船数分析休渔前后的增长幅度。原始捕捞生产数据经SPSS K-S检验不符合正态分布，采用非参数检验-独立样本检验休渔前后短期与长期差异的显著性，显著性水平P值取0.05。综合增长率与显著性检验结果分析休渔前后日均产量产值及CPUE变化。渔获率与渔获结构分析。使用单拖渔船生产数据代入公式(5)和(6)获得单拖渔船的渔获率、渔获结构值，非参数检验-独立样本检验渔船休渔前后渔获数据差异的显著性，综合分析渔获变化情况。
2. 结果
2.1 渔船日均产量产值
自2007年后捕捞总产量变化趋于平缓(图2)，选取2015—2017年休渔前后1个月南海海域作业调查渔船的捕捞数据。虽存在捕捞站点差异，但捕捞渔场重合度较高，且调查渔船数量多，数据具有一定可比性。
![]()
图 2 南海海域20年间 (1997—2016年) 捕捞产量变化Figure 2. Fishing yield in South China Sea over past 20 years (1997−2016)2015年南海区双拖、单拖、桁杆拖虾3种作业日均产量休渔后比休渔前分别增加125.26%、35.86%和13.84%，仅双拖与单拖差异显著，产量变动幅度为双拖>单拖(P<0.05)；2016年桁杆拖虾、围网、单拖、双拖日均产量休渔后比休渔前分别增加125.98%、112.01%、80.64%、49.75%，仅单拖、双拖差异显著(P<0.05)，休渔前后产量变动幅度为单拖>双拖；2017年桁杆拖虾、围网、单拖日均产量休渔后比休渔前分别增加198.94%、118.47%、83.12%，且差异显著(P<0.05)，产量变动幅度为桁杆拖虾>围网>单拖(表1)。
表 1 2015—2017年休渔前后渔船日均产量统计结果Table 1. Statistical results of average daily yield before and after summer fishing moratorium during 2015−2017${\overline {\mathit{\boldsymbol{X}}}} \pm {\bf SD}$ 201505 201508 增长率/%
growth rate201605 201608 增长率/%
growth rate201704 201708 增长率/%
growth rate单拖 single trawler 626.74±674.67a 851.48±790.56b 35.86 634.10±726.81a 1 145.44±1 419.11b 80.64 655.14±845.97a 1 159.57±1 545.25b 77.00 双拖 bull trawler 2 773.50±1 833.75b 6 247.66±4 955.02a 125.26 4 512.46±4 257.73b 6 757.35±5 120.26a 49.75 3 211.07±2 709.62b 5 880.07±5 572.79b 83.12 桁杆拖虾 beam trawling 213.49±177.26b 243.03±218.68b 13.84 187.31±101.98b 423.28±501.54b 125.98 179.42±112.87b 536.36±429.28a 198.94 围网 purse seine 1 462.98±1 664.32ab 1 293.63±2 394.50bc –11.58 1 689.56±2 425.14abc 3 582.12±6 366.51a 112.01 1 933.52±4 208.51c 4 224.16±7 486.65a 118.47 流刺网 gill net 123.85±202.40a 101.45±173.27ab –18.09 115.05±216.32ab 78.75±201.03ab –31.55 45.97±55.52b 78.72±114.06a 71.24 钓具 fishing tackle 133.17±221.25ab 121.82± 190.21b – 8.52 137.02±196.59ab 169.47±251.34ab 23.68 129.08±182.61ab 185.04±236.86a 43.35 其他 other 102.92±123.89a 156.06±510.95a 51.63 79.35±91.59a 159.89±560.15a 101.50 121.08±145.45a 380.79±1 330.72a 214.49 注：不同上标字母表示某一作业类型的日均产量在不同时间差异显著(P<0.05)，下表同此 Note: Different lowercase letters indicate significant difference (P<0.05). The same case in the following tables. 2017年休渔制度将流刺网列入休渔范围，其日均产量休渔后比休渔前增长53.45%，差异显著(P<0.05)。钓具不休渔，因此其变化与休渔相关性较小。其他项包括定置张网、笼壶等作业渔船数量少，且各省休渔规定不同，但受休渔总体影响，日均产量休渔后比休渔前也有增长。3年数据总体对比，休渔对拖网影响较大，单拖、双拖、桁杆拖虾休渔前后差异显著(P<0.05)，增产幅度为桁杆拖虾>双拖>单拖(表1)。
综合所有作业类型分析，2015年、2016年、2017年休渔前后均差异显著且长期差异显著，渔船日均产量休渔后比休渔前分别增加19.21%、49.99%、114.64%，增产效果为2017年>2016年>2015年(表2)。
表 2 2015—2017年南海海域渔船日均产量、日均产值和CPUE均值Table 2. Average values of daily yield, daily production value and CPUE of fishing boats in South China Sea during 2015−2017201505 201508 增长率/%
growth rate201605 201608 增长率/%
growth rate201704 201708 增长率/%
growth rate日均产量/t daily yield 47 445.04c 56 561.28b 19.21 54 049.62c 81 070.71a 49.99 42 504.74d 91 233.06a 114.64 日均产值/百万元 daily production value/million Yuan 317.33c 318.60c 0.40 320.56c 372.03b 16.06 302.27c 450.80a 49.14 单位捕捞努力量渔获量/[kg·(kW·d)–1 CPUE 216 080.46c 231 468.31c 7.12 219 578.73c 380 524.81b 73.30 185 794.37d 391 835.35a 110.90 2015年南海区双拖、单拖、桁杆拖虾3种作业类型渔船休渔后日均产值比休渔前分别增加49.01%、34.41%、26.91%，仅单拖差异显著(P<0.05)；2016年单拖、双拖休渔后日均产值比休渔前分别增加34.11%、28.88%，且差异显著(P<0.05)，产值变动幅度为单拖>双拖。2017年桁杆拖虾、单拖、双拖日均产值休渔后比休渔前分别增加106.25%、53.67%、43.99%，差异显著(P<0.05)，产值增长幅度为桁杆拖虾>单拖>双拖。2017年休渔制度将流刺网列入休渔范围，其在休渔后产值增长，差异不明显。围网日均产值在休渔后减少，差异不显著。钓具与其他作业类型在休渔前后产值差异均不明显，受休渔制度影响较小。3年数据总体对比，单拖、双拖、桁杆拖虾差异较显著，休渔后产值增加，增产幅度为桁杆拖虾>双拖>单拖>围网。
综合所有作业类型分析，2017年、2016年休渔前后差异显著，2015年休渔前后差异不显著，3年变化不显著，2015年、2016年、2017年渔获物日均产值休渔后比休渔前分别增加0.40%、16.06%、49.14%，增值效果为2017年>2016年>2015年(表3)。
表 3 2015—2017年休渔前后渔船日均产值统计结果Table 3. Statistical results of average daily production value before and after summer fishing moratorium during 2015−2017${\overline {\mathit{\boldsymbol{X}}}} \pm {\bf SD}$ 作业类型
type of fishing vessels201505 201508 增长率/%
growth rate201605 201608 增长率/%
growth rate201704 201708 增长率/%
growth rate单拖 single trawler 5 763.49±2 860.79c 7 746.87±3 701.79b 34.41 6 230.31±2 470.67c 8 355.57±4 888.65b 34.11 6 286.32±3 257.39c 9 660.15±5 149.16a 53.67 双拖 bull trawler 14 055.21±10 375.48ab 20 943.16±15 242.59ab 49.01 15 082.40±9 090.37b 19 438.78±10 601.54a 28.88 16 375.87±8 845.72ab 23 580.17±16 977.61c 43.99 桁杆拖虾 beam trawling 3 392.19±1 973.87b 4 305.07±4 713.02b 26.91 3 092.06±1 777.26b 5 331.26±4 538.99b 72.42 3 037.33±1 701.58b 8 040.85±5 343.70a 164.73 围网 purse seine 8 348.29±8 310.18a 5 154.69±9 763.42b –38.25 9 428.75±14 645.69a 8 635.25±11 685.94a –8.42 11 934.76±20 451.73a 10 608.54±9 966.25a –11.11 流刺网 gill net 1 535.84±1 912.44ab 1 374.93±2 232.75a –10.48 1 543.25±1 915.26a 1 276.68±2 118.26a –17.27 1 248.45±1 740.08a 1 727.16±2 637.69a 38.34 钓具 fishing tackle 4 042.12±5 802.27ab 3 447.60±5 047.25a –14.71 4 105.74± 4 660.40a 4 732.23±6 306.55a 15.26 4 048.59±5 001.51a 5 757.97±6 766.10a 42.22 其他 other 1 163.73±2 418.94a 2 863.54±6 620.41a 146.07 2 800.96±3 783.53a 2 838.89±5 002.60a 1.35 3 381.70±6 550.63a 5 069.43±10 006.12a 49.91 2.2 捕捞努力量
2015年南海区双拖和单拖CPUE休渔后比休渔前分别增加81.34%和61.18%，且差异显著(P<0.05)，变动幅度为双拖>单拖。2016年围网、桁杆拖虾、单拖休渔后的CPUE分别增加194.24%、76.25%、62.28%，且差异显著(P<0.05)，变动幅度为围网>桁杆拖虾>单拖>双拖。2017年休渔制度调整后，围网、桁杆拖虾、单拖、流刺网和双拖的CPUE休渔后比休渔前分别增加184.08%、144.27%、72.30%、63.95%和72.30%，所有休渔类型渔船的CPUE均差异显著(P<0.05)，休渔前后变动幅度为围网>桁杆拖虾>单拖>流刺网>双拖。2017年流刺网被列入休渔范围，其CPUE在休渔后增长且差异显著(P<0.05)。就一年周期而言，休渔后多数作业渔船的CPUE增加，其中单拖、围网、桁杆拖虾、流刺网等差异显著。
综合所有作业类型分析，2015年休渔前后渔船CPUE差异不明显，2017年、2016年差异显著，2015年、2016年、2017年增长率分别为7.12%、73.30%、110.90%，CPUE增长幅度为2017年>2016年>2015年(表4)。
表 4 2015—2017年休渔前后各作业类型CPUE统计结果Table 4. Statistical results of CPUE of each type of fishing vessels before and after summer fishing moratorium during 2015−2017${\overline {\mathit{\boldsymbol{X}}}} \pm {\bf SD}$ 作业类型
type of fishing vessels201505 201508 增长率/%
growth rate201605 201608 增长率/%
growth rate201704 201708 增长率/%
growth rate单拖 single trawler 2.55±1.67b 4.11±2.67a 61.18 2.89±2.17b 4.69 ±4.08a 62.28 2.78±2.57b 4.79 ±4.64a 72.30 双拖 bull trawler 8.63±5.00c 15.65±8.46a 81.34 11.47±6.63bc 13.98±5.19ab –93.61 8.69± 6.51c 11.57± 3.98b 33.14 桁杆拖虾 beam trawling 1.89±1.51bc 2.01±1.32b 6.35 1.60±1.06bc 2.82±2.55a 76.25 1.31±0.97c 3.20±2.08a 144.27 围网 purse seine 8.23±8.85b 6.50±11.76c –99.99 9.37±11.00b 27.57±39.54a 194.24 10.49±21.72c 29.80±41.47a 184.08 流刺网 gill net 2.33±3.10b 2.27±2.62b –2.58 2.53±3.59b 2.08±2.01b –17.79 1.72±1.96b 2.82±3.47a 63.95 钓具 fishing tackle 1.93±1.91ab 2.03±1.55ab 5.18 1.92±1.38ab 2.28±1.82ab 18.75 1.63±1.23b 2.19±1.51a 23.00 2.3 单拖渔船的渔获结构和渔获率
单拖渔船为南海海域传统作业类型，且受休渔影响较为显著，以单拖渔获为例分析休渔变化情况。调查渔船基本信息见表5。
表 5 单拖渔船情况Table 5. Trawling survey boats information年份
year月份
month船次
Voy. No.马力
horsepower出海天数
days at sea2015 5 134 40 858 930 8 154 42 087 1 043 2016 5 126 35 537 686 8 140 41 247 906 2017 4 165 46 882 1 351 8 136 39 121 629 注：船次为所有调查渔船该月出海总次数，存在同艘渔船该月多次出海情况 Note: Voy. No. is the total number of fishing boats going out at sea that month, and there are some fishing boats going out several times within that month. 一年周期内，抽样调查船中主要经济鱼类休渔后渔获比例明显上升，其中以蓝子鱼最为明显，渔获量由2015年的1 265 kg增至84 385 kg，渔获比例由1.22%增至21.15%，且休渔前后差异显著(P<0.05)；鲹类、二长棘鲷、带鱼等重要经济鱼类渔获率增加，且在2016年、2017年休渔前后差异显著(P<0.05)。就长期效果而言，7种经济鱼类总渔获率由2015年的0.08 kg·(kW·d)–1增至2017年的1.16 kg·(kW·d)–1，增幅为1 335.62%，所占渔获比例从2015年的40.23%增至2017年的47.52%(表6)，种间结构出现恢复迹象。但3年来经济鱼类渔获率与渔获结构长期差异不显著(P>0.05)，休渔对渔获率与渔获结构影响仅局限于当年，长期效果难以巩固；且南海海域生产船次与总功率在休渔后迅速上升，甚至在2017年休渔后半个月便达到去年同期水平(表5)，失控的船数与功率形成了巨大的捕捞压力。
表 6 2015—2017年主要渔获种类渔获率与渔获结构Table 6. Fishing rate and fishing structure of main fish species during 2015−2017种类
species渔获率/kg·(kW·d)–1 catch rate 渔获结构/% fishing structure 2015年 2016年 2017年 2015年 2016年 2017年 5月 8月 5月 8月 4月 8月 5月 8月 5月 8月 4月 8月 蓝子鱼 Siganus 0.00 0.00 0.00 0.01 0.00 0.52 1.22 0.36 0.08 0.44 0.44 21.15 鲹类 Carangidae 0.00 0.16 0.01 0.92 0.02 0.14 1.45 11.67 0.97 53.17 2.79 5.67 二长棘鲷 Paerargyrops edita 0.01 0.06 0.00 0.08 0.00 0.07 2.69 4.72 0.47 4.58 0.38 3.02 带鱼 Trichiurus haumela 0.05 0.19 0.05 0.16 0.09 0.32 26.64 14.08 7.63 9.41 10.88 13.29 金线鱼 Nemipterus virgatus 0.00 0.13 0.01 0.03 0.00 0.02 2.03 9.37 1.49 1.70 0.22 0.78 刺鲳 Psenopsis anomala 0.01 0.03 0.00 0.04 0.00 0.02 4.65 2.15 0.76 2.04 0.10 0.69 鱿鱼 Loligo chinensis 0.00 0.05 0.01 0.03 0.01 0.07 1.54 3.50 1.65 1.55 0.66 2.92 3. 讨论
3.1 新调整伏季休渔制度的效果
2017年南海休渔时间延长为5月1日12时至8月16日12时，范围扩大至除钓具外的所有作业类型，捕捞辅助船首次同步休渔[5]。本研究发现，2017年伏季休渔效果优于2015年与2016年，南海海域休渔作业类型渔船日均产量、产值和CPUE增长率分别为114.64%、49.14%和110.90%，明显大于前2年，且休渔前后均差异显著。首年休渔的流刺网渔船日均产量、产值及CPUE增幅也效果明显。新休渔制度使渔民增产增收，且比2015年、2016年效果更好，这应该与休渔时间的延长、类型的扩展有关。至于新制度对资源种群结构的改善效果还需进一步研究。
3.2 经济效益
因各年休渔期差异，采用日均产量、日均产值的“平均值±标准差(
$\overline X \pm {\rm SD} $ )”和增长率来反映休渔影响幅度，通过非参数检验验证差异显著性 (P取0.05)，可为休渔效果验证提供理论依据。该方法因各省渔船功率及作业效率不同，同组数据离散，标准差较大，结果可能存在误差，但整体趋势不受影响。渔业资源极易受捕捞和环境变化影响，在经济效益上直观表现为产量呈较大波动。陈春亮等[15]和侯秀琼等[16]发现休渔后珠江口掺缯作业和桁杆拖虾的渔获量分别增加了10.8倍和5.6倍，程家骅[17]认为珠江口、北部湾的休渔带来了产量的暂时增长，这些研究均肯定了休渔具有增产的作用。本研究显示，2015年、2016年、2017年渔船日均产量休渔后比休渔前分别增加19.21%、49.99%、114.64%，增幅为2017年>2016年>2015年，休渔前后产量差异显著，且3年增产效果显著，但对产值影响有限，同年伏休后产值有所增长但长期效果不显著。张龙等[18]认为休渔后主要经济种类渔获个体质量大幅增加。Arendse等[19]认为对于易受捕捞干扰或洄游聚集繁殖而易被捕捞的鱼类，休渔可提高其增殖和增产效果。以上观点主要从抽样鱼类生长变化的生物学角度出发，而本文基于渔船的捕捞数据，从不同角度佐证了休渔是渔获数量与质量提高的主要原因。余景等[10]认为休渔期时大亚湾海域具有适宜的环境和相对充足的饵料，Davies等[20]发现禁渔区的幼鱼生存率有所提高，本文均进一步论证了本文上述观点。
针对产值变动，潘澎和李卫东[21]认为休渔缩短了作业时间，减少了燃油、人工、网具等的消耗；而徐莲莲和杨美丽[22]则认为休渔期渔船无法出海作业导致收益减少。以上观点从不同角度印证休渔确实能够影响产值，因国内相关定量研究极少，较难分析其中差异。本研究量化分析了休渔的增值效果，但受数据及方法限制无法进一步探究深层原因，且因渔船功率、船长及不同省区捕捞效率的差异，可能存在产值衡量标准不同。国内关于休渔制度下船队业绩、渔民收益方面的定量研究极少，而国外相关研究较多倾向于经济效益与资源保护的平衡[23]，值得进一步借鉴。
3.3 生态效益
本研究表明，一年周期内单拖、围网、桁杆拖虾、流刺网等渔船的CPUE增加且休渔前后差异显著。围网与刺网捕捞对象多为中上层鱼类，桁杆拖虾、单双拖捕捞对象多为底层和近底层鱼类，这些作业渔船的CPUE在休渔后增幅明显，也表明各水层的资源量有所增加[24]，资源结构改善。余景等[10]指出休渔后中上层及底层主要渔获CPUE增加明显；邹建伟等[11-12]研究发现2016年休渔后拖网、桁杆拖虾和灯光罩网CPUE比2015年分别增加了28.75%、65.14%和13.35%，拖网增幅最大，与本研究CPUE增幅2017年>2016年>2015年的结果相符。
休渔期是大部分鱼类的繁殖期与生长期[25]，有利于改善资源的种群结构。本研究发现蓝子鱼、带鱼、鲹类等拖网主要经济鱼类资源在休渔后有所恢复，部分重要经济鱼类的资源衰退情况也得到了一定的遏制。邹建伟等[26]也发现北部湾休渔后主要经济鱼类渔获率增加，本研究结论与之相似。晏磊等[27]指出从抽样鱼类生物学角度看，休渔后虽然渔获种类增加，但小型个体或大型鱼类幼体占多数；Watson等[28-29]发现禁渔区内的主要渔获种类的丰度、体长和生物量均明显占优。上述关于休渔后渔获种类增加的观点与本研究一致，但渔获体长的观点则与之相悖，究其原因可能是由各国休渔时长、环境的差异造成。本研究仅从渔获种类、渔获量角度分析，缺乏生物学数据，也无法探究休渔后鱼类群落结构特征。虽然资源密度与经济鱼类种类在休渔后增加，但体长、年龄结构、质量则会影响种群结构的稳定性与持久性，仍需进一步探究。
3.4 制度局限性
休渔制度很早便在国外实施用以保护渔业资源[30-31]，但各国在实践中均发现休渔难以遏制资源的衰退现状[32-34]。本研究表明，虽然休渔的短期增产效果较好，但对CPUE、渔获率与渔获结构的作用仅局限于当年，长期效果难以巩固，与邹建伟等[26]、蔡研聪等[35]关于资源量、资源密度的观点相符。王迎宾等[36]、刘凯等[37]采用鱼类群落的多样性指数与营养级提出休渔的保护成果到年末便不复存在；李忠炉等[38]、刘尊雷等[39]使用Ricker模型发现冬季鱼类资源密度跌至较低水平，个体体长比质量增长更快。以上学者多从渔获生物学特征变动或多样性角度考虑，而本研究基于3年捕捞渔船的CPUE、渔获率与渔获结构数据定量分析休渔的长期效果，能更全面地反映长期变化，但仅有休渔前后1个月的数据，难以探究休渔的持续性，需积累更多数据加以佐证。
中国的休渔制度实质上是一种“投入控制”制度，理论上只要保持捕捞量不超过渔业资源的再生量，渔业资源就不会被破坏。然而在实践中休渔制度很难、甚至不可能有效地控制渔业的捕捞力量[40]，无法解决渔船数量增加、马力增大和捕捞技术提高所带来的问题，且开捕后捕捞强度的增加也会抵消休渔效果[41]。本调查显示，休渔后南海海域渔船总船次与总功率呈迅猛增长态势，2017年8月开捕后仅半个月，渔船船次和总功率分别达2016年8月的97.14%和94.85%。在强大的捕捞压力下，几种经济鱼类资源量虽有所恢复，但渔获率与渔获结构的长期差异不显著，休渔难以改善资源衰退现状。
-
![]()
图 1 2006—2014年北部湾大头白姑鱼种群参数的变化趋势
Figure 1. Trends of population parameters of P. macrocephalus in Beibu Gulf from 2006 to 2014
![]()
图 2 大头白姑鱼捕捞强度 (F/M) 与资源利用状态
Figure 2. Fishing intensity (F/M) and stock utilization status of P. macrocephalus in Beibu Gulf from 2006 to 2014
表 1 采样航次、时间与样本量信息
Table 1 Information of sampling cruise, sampling month and sampling size
年份
Year月份
Month样本数
Sampling size采样航次
Sampling cruise2006 2 167 第1航次 8 695 第2航次 2007 2 937 第3航次 7 412 第4航次 2008 1 241 第5航次 7 280 第6航次 2009 2 328 第7航次 8 336 第8航次 2010 1 485 第9航次 8 335 第10航次 2011 1 441 第11航次 8 336 第12航次 2012 1 202 第13航次 8 298 第14航次 2013 1 383 第15航次 7 358 第16航次 2014 1 324 第17航次 7 34 第18航次 
下载: 导出CSV
表 2 2006—2014年大头白姑鱼群体结构
Table 2 Population structure of P. macrocephalus in summer and winter from 2006 to 2014
年份
Year冬季体长
Body length in winter/mm夏季体长
Body length in summer/mm范围
Range优势
Dominant平均
Mean范围
Range优势
Dominant平均
Mean2006 98~233 140~149 146 40~230 120~129 129 2007 50~278 150~159 140 39~215 100~109 110 2008 65~230 170~179 145 25~214 100~109 130 2009 68~245 140~149 132 69~212 110~119 121 2010 59~253 150~159 163 60~220 100~109 121 2011 45~290 100~109 114 55~220 90~99 116 2012 36~212 100~109 105 6.9~240 100~109 119 2013 47~232 130~139 129 28~240 110~119 130 2014 25~240 140~149 142 88~220 110~119 120 2006—2014 25~290 140~149 137 6.9~240 110~119 122
下载: 导出CSV
表 3 2006—2014年根据LBB估算的大头白姑鱼的种群参数
Table 3 Population parameters of P. macrocephalus estimated by LBB from 2006 to 2014
年份
Year月份
Month渐进体长
L∞/mm最适
开捕体长
Lc_opt/mm相对自然
死亡率
M/k最小值
Min最大值
Max相对捕捞
死亡率
F/k最小值
Min最大值
Max相对
总死亡率
Z/k最小值
Min最大值
Max开发率
E2006 2 275.00 142.00 1.31 1.07 1.46 0.37 0.21 0.64 1.67 1.57 1.77 0.22 8 279.00 169.00 1.36 1.13 1.57 2.34 2.13 2.72 3.7 3.53 3.92 0.63 2007 2 277.00 146.00 1.63 1.39 1.87 1.36 1.05 1.67 2.98 2.72 3.31 0.46 7 258.00 147.00 1.61 1.33 1.9 2.68 2.32 3.00 4.26 4.09 4.48 0.63 2008 1 252.00 136.00 1.50 1.21 1.75 1.13 0.81 1.39 2.59 2.35 2.75 0.44 7 225.00 122.00 0.91 0.82 1.02 0.12 0.06 0.22 1.04 0.95 1.12 0.12 2009 2 291.00 197.00 1.22 1.00 1.39 7.15 6.43 7.94 8.38 7.66 9.08 0.85 8 258.00 149.00 1.56 1.31 1.87 2.71 2.31 3.19 4.29 4.04 4.56 0.63 2010 1 303.00 172.00 1.57 1.27 1.85 2.30 1.97 2.81 3.88 3.61 4.16 0.59 8 269.00 160.00 1.55 1.31 1.88 3.55 3.14 3.84 5.11 4.85 5.32 0.69 2011 1 307.00 180.00 1.56 1.24 1.87 3.19 2.74 3.61 4.72 4.42 5.11 0.68 8 234.00 114.00 1.49 1.32 1.69 0.36 0.16 0.60 1.87 1.77 1.97 0.19 2012 1 213.00 145.00 1.22 0.80 1.49 7.93 7.01 9.12 9.09 8.27 10.30 0.87 8 285.00 156.00 1.50 1.18 1.71 1.34 1.06 1.66 2.83 2.70 2.95 0.47 2013 1 240.00 136.00 1.74 1.51 1.95 3.58 3.08 4.21 5.27 4.96 5.83 0.68 7 282.00 147.00 1.51 1.29 1.76 0.85 0.58 1.08 2.36 2.20 2.50 0.36 2014 1 242.00 226.00 0.15 −0.07 0.44 3.47 2.68 4.03 3.60 3.07 4.03 0.96 7 271.00 173.00 1.57 1.28 1.87 20.00 18.00 21.50 21.60 19.50 23.10 0.93 平均值 Mean 264.50 156.50 1.39 1.13 1.63 3.58 3.10 4.07 4.96 4.57 5.35 0.58 变异系数 CV 0.10 0.17 0.26 0.31 0.23 1.29 1.35 1.22 0.94 0.92 0.94 0.43
下载: 导出CSV
表 4 白姑鱼属种群参数的估算结果
Table 4 Estimation results of population parameters of Pennahia
种类
Species年份
Year渐进体长
L∞/mm最适开捕体长
Lc_opt/mm生长率
k相对自然死亡率
M/k开发率
E方法
Method文献
Reference白姑鱼 P. argentatus 1982 305.00 − 0.35 2.29 0.23 生长方程 [31] 1994 − − − − 0.50 白姑鱼 P. argentata 1983—1984 282.30 − 0.50 − − 生长方程 [32] 白姑鱼 P. argentata 1992—1993 382.00 211.00 0.42 2.21 0.74 ELEFAN I [19] 大头白姑鱼 P. macrocephalus 2006—2007 229.95 − 0.41 − − 最大似然法 [14] 斑鳍白姑鱼 P. pawak 2008—2009 220.32 146.53 0.58 − − 最大似然法 [20] 大头白姑鱼 P. macrocephalus 2006—2007 289.00 174.00 0.52 2.00 0.63 ELEFAN I [28] 大头白姑鱼 P. macrocephalus 2006—2014 264.50 156.50 − 1.49 0.58 LBB 本研究
下载: 导出CSV
-
[1] 史登福, 张魁, 陈作志. 基于生活史特征的数据有限条件下渔业资源评估方法比较[J]. 中国水产科学, 2020, 27(1): 12-24. [2] RANKIN P S, LEMOS R T. An alternative surplus production model[J]. Ecol Model, 2015, 313: 109-126. doi: 10.1016/j.ecolmodel.2015.06.024
[3] 冯波, 陈新军, 西田勤. 应用年龄结构产量模型评估印度洋黄鳍金枪鱼资源[J]. 生态学报, 2010, 30(13): 3375-3384. [4] 徐洁, 陈新军, 官文江. 适用于短生命周期种类资源评估模型的研究现状与展望[J]. 海洋湖沼通报, 2015(3): 113-124. [5] CHRYSAFI A, KUPARINEN A. Assessing abundance of populations with limited data: lessons learned from data-poor fisheries stock assessment[J]. Environ Rev, 2016, 24(1): 25-38. doi: 10.1139/er-2015-0044
[6] PRINCE J, VICTOR S, KLOULCHAD V, et al. Length based SPR assessment of eleven Indo-Pacific coral reef fish populations in Palau[J]. Fish Res, 2015, 171: 42-58. doi: 10.1016/j.fishres.2015.06.008
[7] BOWKER D W. The bias associated with the dispersion of length at age in the estimation of growth and mortality parameters from uni-modal length-frequency data[J]. J Fish Biol, 1996, 49(4): 743-747.
[8] 汪文婷, 田思泉, 戴小杰, 等. 中大西洋延绳钓渔业大眼金枪鱼体长频率时空分布[J]. 上海海洋大学学报, 2013, 22(5): 770-777. [9] FROESE R, WINKER H, CORO G, et al. A new approach for estimating stock status from length frequency data[J]. ICES J Mar Sci, 2018, 75(6): 2004-2015. doi: 10.1093/icesjms/fsy078
[10] BALDE B S, FALL M, KANTOUSSAN J, et al. Fish-length based indicators for improved management of the sardinella fisheries in Senegal[J]. Reg Stud Mar Sci, 2019, 31: 100801. doi: 10.1016/j.rsma.2019.100801
[11] 潘澎, 李纯厚, 罗家聪, 等. 中越北部湾捕捞力量变化问题研究[J]. 中国渔业经济, 2016, 34(5): 4-9. [12] 王雪辉, 邱永松, 杜飞雁, 等. 北部湾鱼类多样性及优势种的时空变化[J]. 中国水产科学, 2011, 18(2): 427-436. [13] 袁华荣, 陈丕茂, 贾晓平, 等. 北部湾东北部游泳生物资源现状[J]. 南方水产科学, 2011, 7(3): 31-38. [14] 颜云榕, 侯刚, 卢伙胜, 等. 北部湾大头白姑鱼生长特性及群体组成[J]. 中国海洋大学学报 (自然科学版), 2010, 40(6): 61-68. [15] 苏新红. 聚类分析在石首鱼类分类中的应用[J]. 生物数学学报, 2009, 24(4): 733-740. [16] YAMADA H, YAMADA U. Descriptive morphology of juvenile stages of two sciaenids, Miichthys miiuy and Pennahia macrocephalus, from the East China Sea[J]. Ichthyol Res, 1999, 46(1): 93-99. doi: 10.1007/BF02674952
[17] SASAKI K. Johnius aneus Bloch, a senior synonym of Pennahia macrophthalmus (Bleeker), with comments on the identity and status of an alleged lectotype of J. aneus (Sciaenidae: Perciformes)[J]. 魚類学雑誌, 1994, 40(4): 498-499. [18] 张其永, 林秋眠, 林尤通, 等. 闽南—台湾浅滩渔场经济鱼类食物网研究 (摘要)[J]. 福建水产科技, 1980(1): 17-18, 21. [19] 陈作志, 邱永松, 黄梓荣. 南海北部白姑鱼生长和死亡参数的估算[J]. 应用生态学报, 2005, 16(4): 712-716. [20] 颜云榕, 侯刚, 卢伙胜, 等. 北部湾斑鳍白姑鱼的年龄与生长[J]. 中国水产科学, 2011, 18(1): 145-155. [21] 何雄波, 陶雅晋, 侯刚, 等. 北部湾斑鳍白姑鱼群体结构特征与时空分布[J]. 广东海洋大学学报, 2015, 35(3): 35-42. [22] SIEGFRIED K I, SANSO B. Two Bayesian methods for estimating parameters of the von Bertalanffy growth equation[J]. Environ Biol Fish, 2006, 77(3): 301-308.
[23] FROESE R, WINKER H, GASCUEL D, et al. Minimizing the impact of fishing[J]. Fish Fish, 2016, 17(3): 785-802. doi: 10.1111/faf.12146
[24] PLUMMER M. JAGS: a program for analysis of Bayesian graphical models using Gibbs sampling[EB/OL]. [2020-10-12]. https://www.r-project.org/conferences/DSC-2003/Proceedings/Plummer.pdf.
[25] 联合国粮食与农业组织. 2018年世界渔业和水产养殖状况—实现可持续发展目标. ISBN 978-92-5-130691-8[R]. 罗马: 联合国粮食及农业组织, 2018. [26] 叶婷, 王迎宾, 周丛羽. 鱼类体长频率数据结构对生长参数估算的影响分析[J]. 水产科学, 2014, 33(5): 277-282. [27] WANG Y G, ELLIS N. Effect of individual variability on estimation of population parameters from length-frequency data[J]. Can J Fish Aquat Sci, 1998, 55: 2393-2401. doi: 10.1139/f98-134
[28] WANG X H, QIU Y S, DU F Y, et al. Population parameters and dynamic pool models of commercial fishes in the Beibu Gulf, northern South China Sea[J]. Chin J Oceanol Limnol, 2012, 30(1): 105-117. doi: 10.1007/s00343-012-1017-y
[29] RAMOS-MIRANDA J. Growth, mortality, maturity, and recruitment of the star drum (Stellifer lanceolatus) in the southern Gulf of Mexico[J]. Ceinc Mar, 2009, 35(3): 245-257.
[30] 李壮, 刘群. 应用渔业体长分析方法ELEFAN和SLCA估算鱼类生长参数的研究[J]. 海洋湖沼通报, 2007(3): 81-87. [31] 卢振彬, 戴泉水, 颜尤明. 福建近海主要底层经济鱼类的种群动态[J]. 台湾海峡, 1999, 18(1): 100-105. [32] 胡雅竹, 钱世勤. 白姑鱼年龄和生长的研究[J]. 海洋渔业, 1989(4): 158-162. [33] 陈丕茂. 南海北部主要捕捞种类最适开捕规格研究[J]. 水产学报, 2004, 28(4): 393-400. [34] 王雪辉, 邱永松, 杜飞雁, 等. 基于长度贝叶斯生物量法估算北部湾二长棘鲷种群参数[J]. 水产学报, 2020, 44(10): 1654-1662. [35] WANG Y B, WANG Y C, LIU S D, et al. Stock assessment using LBB method for eight fish species from the Bohai and Yellow Seas[J]. Front Mar Sci, 2020, 7: 164.
[36] GULLAND J A. Fish stock assessment: a manual of basic methods[M]. Chichester, UK: FAO/Wiley series on food and agriculture, 1983: 1-223.
[37] PONS M, COPE J M, KELL L T. Comparing performance of catch-based and length-based stock assessment methods in data-limited fisheries[J]. Can J Fish Aquat Sci, 2020, 77(6): 1026-1037. doi: 10.1139/cjfas-2019-0276
[38] WANG J Q, XU L L, LI B, et al. An evaluation of the effects of sample size on estimating length composition of catches from tuna longline fisheries using computer simulations[J]. Aqua Fish, 2020, 5(3): 122-130. doi: 10.1016/j.aaf.2019.09.001
[39] 李寇军, 邱永松, 王跃中. 自然环境变动对北部湾渔业资源的影响[J]. 南方水产, 2007, 3(1): 7-13. [40] 陈国宝, 李永振, 陈丕茂, 等. 鱼类最佳体长频率分析组距研究[J]. 中国水产科学, 2008, 15(4): 659-666. -
期刊类型引用(16)
1. 冯雪，戴小杰，袁华荣，成志健，陈丕茂. 外伶仃海洋牧场附近海域游泳生物群落结构及其与环境因子的关系. 上海海洋大学学报. 2024(01): 186-201 . 
百度学术
2. 田思泉，柳晓雪，花传祥，王寅，杜涣洋. 南海渔业资源状况及其管理挑战. 上海海洋大学学报. 2024(03): 786-798 . 百度学术
3. 樊菲，郑建丽，吴姗姗，张祝利. 我国伏季休渔制度效果分析与优化探讨. 海洋开发与管理. 2024(04): 152-160 . 百度学术
4. 许罕多，陈园园. 中国近海渔业资源衰退的动态评价研究. 海洋经济. 2024(05): 27-37 . 百度学术
5. 孙康，冯志鹏. 渔民对伏季休渔政策的满意度研究——基于黄、渤海渔村问卷调查的分析. 海洋开发与管理. 2023(01): 29-36 . 百度学术
6. 尹元一，章丽娟，廖东昇，李璐，谢文锋，冯波. 南海北部刺鲳最大可持续产量的评估. 渔业研究. 2023(04): 365-371 . 百度学术
7. 张磊，周艳波，马胜伟，黄应邦，高丽鹏，杜国昱，吴洽儿. 南海伏季休渔期秩序评估体系指标权重分析. 上海海洋大学学报. 2022(02): 491-501 . 百度学术
8. 蒋明峰，陈新军. 中东太平洋赤道公海海域鱿钓渔业自主休渔效果初步评价. 上海海洋大学学报. 2022(03): 670-676 . 百度学术
9. 郭迎香，杨李玲，许友伟，方艺菲，王萌，姜敬哲. 南海近岸硬骨鱼鳃组织的细菌群落及多样性分析. 南方水产科学. 2022(04): 44-53 . 本站查看
10. Qiang Wu，Qingpeng Han，Yue Jin，Tao Yang，Zhongyi Li，Xiujuan Shan，Xianshi Jin. Long-term and monthly changes in abundance, size composition and spatial distribution of the mantis shrimp Oratosquilla oratoria in the Bohai Sea. Acta Oceanologica Sinica. 2022(12): 80-86 . 必应学术
11. 郭琇，麦琳，张泽伟，林冲. 复合载体生物流化床在循环水养殖系统的调控效果. 广东化工. 2021(05): 111-114 . 百度学术
12. 张羽翔，刘维，赵海龙，吴程宏，陈敏. 海南岛近岸春秋季游泳动物多样性分析. 渔业研究. 2021(02): 133-142 . 百度学术
13. 吴程宏，张羽翔，刘维，赵海龙，陈敏. 基于上岸渔获调查的海南岛近海渔场伏季休渔效果评价. 渔业研究. 2021(02): 200-206 . 百度学术
14. 苏雪，陈新军. 渔业管理中投入和产出控制的比较研究. 海洋湖沼通报. 2021(03): 136-144 . 百度学术
15. 史登福，许友伟，孙铭帅，黄梓荣，陈作志，张魁. 广东海洋渔业资源可捕量评估. 海洋渔业. 2021(05): 521-531 . 百度学术
16. 廖维晓，李督. 南海渔业资源发展机制研究. 现代管理科学. 2019(12): 39-41 . 百度学术
其他类型引用(11)
计量
- 文章访问数: 817
- HTML全文浏览量: 266
- PDF下载量: 74
- 被引次数: 27
粤公网安备 44010502001741号
