巴基斯坦海洋渔业资源可捕量评估与开发现状

张魁, 耿平, 史登福, 许友伟, SHER KHANPanhwar, 陈作志

张魁, 耿平, 史登福, 许友伟, SHER KHANPanhwar, 陈作志. 巴基斯坦海洋渔业资源可捕量评估与开发现状[J]. 南方水产科学, 2019, 15(5): 1-8. DOI: 10.12131/20190065
引用本文: 张魁, 耿平, 史登福, 许友伟, SHER KHANPanhwar, 陈作志. 巴基斯坦海洋渔业资源可捕量评估与开发现状[J]. 南方水产科学, 2019, 15(5): 1-8. DOI: 10.12131/20190065
ZHANG Kui, GENG Ping, SHI Dengfu, XU Youwei, SHER KHAN Panhwar, CHEN Zuozhi. Assessment of allowable catches and exploitation state of marine fishery resources in Pakistani coastal waters[J]. South China Fisheries Science, 2019, 15(5): 1-8. DOI: 10.12131/20190065
Citation: ZHANG Kui, GENG Ping, SHI Dengfu, XU Youwei, SHER KHAN Panhwar, CHEN Zuozhi. Assessment of allowable catches and exploitation state of marine fishery resources in Pakistani coastal waters[J]. South China Fisheries Science, 2019, 15(5): 1-8. DOI: 10.12131/20190065

巴基斯坦海洋渔业资源可捕量评估与开发现状

基金项目: 国家自然科学基金项目 (31602157);中国水产科学研究院基本科研业务费专项资金 (2018GH03);2018年广东省促进经济发展专项资金 (海洋经济发展用途) 项目 (GDME-2018E004)
详细信息
    作者简介:

    张 魁 (1987—),男,博士,助理研究员,从事渔业资源评估研究。E-mail: nedvedkui@163.com

    通讯作者:

    陈作志 (1978—),男,博士,研究员,从事渔业资源和海洋生态研究。E-mail: zzchen2000@163.com

  • 中图分类号: S 937.3

Assessment of allowable catches and exploitation state of marine fishery resources in Pakistani coastal waters

  • 摘要:

    基于巴基斯坦1950—2015年海洋渔业统计产量数据,利用Catch-MSY模型对巴基斯坦海洋渔业总可捕量以及24个重要经济类群的最大可持续产量 (maximum sustainable yield, MSY) 和可捕量进行了评估。结果显示,内禀增长率(r)先验分布对MSY的评估结果影响不大,巴基斯坦海洋渔业MSY为40.53×104 t,总可捕量为36.47×104 t。2015年海洋渔业产量为36.10×104 t,目前未处于过度捕捞状态。24个重要经济类群的评估结果显示,有8个类群 (军曹鱼、魣类、石斑鱼类、鲳类、鲯鳅、宝刀鱼、马鲹和白带鱼) 2015年产量超过MSY,处于过度捕捞状态。当前巴基斯坦海洋渔业资源已处于充分开发状态,多数经济类群已处于过度捕捞状态或者崩溃后的恢复状态,开发潜力较小,建议开发阿拉伯海鸢乌贼(Sthenoteuthis oualaniensis)等大洋性渔业资源。

    Abstract:

    According to the catch statistics data during 1950−2015, we applied a Catch-MSY model to estimate the allowable catches and maximum sustainable yields (MSYs) for 24 commercial fish groups and total allowable catch (TAC) in Pakistani coastal waters. Results show that MSY and TAC in Pakistani coastal waters were 40.53×104 t and 36.47×104 t, respectively. The total catches in 2015 were less than MSY, suggesting that the fisheries are not overfished. The MSY results show that eight fish groups (cobia, Barracudas, Groupers, Butterfishes, common dolphinfish, dorab wolf-herring, crevalle jack and largehead hairtail) have been overfished at present since the catches in 2015 exceeded MSY. At present, the marine fishery resources in Pakistani coastal waters are fully exploited. Most commercial fish groups are overfished or recovering after collapse and their development potential is relatively small. It is recommended to exploit oceanic fishery resources such as Sthenoteuthis oualaniensis in Arabian Sea in future fisheries development.

  • 巴基斯坦海域位于阿拉伯海北部,与伊朗和印度所属海域相邻,拥有1 100 km的海岸线以及24×104 km2的专属经济区[1]。海洋渔业在国民经济中起着重要作用,体现在提供动物蛋白、解决就业和出口创汇等方面。海洋渔业产量自20世纪60年代到20世纪末为高速增长阶段,产量从1961年的5.38×104 t增加到1999年的47.47×104 t,1993年产量最高,为49.91×104 t,近10余年产量维持在34.80×104 t左右 (图1)。2016年,巴基斯坦海洋渔业产量为37.63×104 t,占捕捞水产品的73.3%和全部水产品的56.2%,渔业产品出口额为3.52亿美元,占商品出口贸易总额的1.7%[2]。巴基斯坦的海洋渔业资源捕捞技术比较落后,其海洋渔业主管部门积极推动了多种渔业发展模式,包括开放专属经济区与他国渔船合作开发海洋渔业资源[3-4],为我国远洋渔业发展和“一带一路”倡议的实施带来机遇[5]

    图  1  1950—2016年巴基斯坦海洋渔业统计产量
    Figure  1.  Statistic data of marine fisheries catches in Pakistani waters from 1950 to 2016

    目前,巴基斯坦海洋渔业资源中仅有少数种类存在资源评估[1-2, 6-8],大部分经济种类缺乏最大可持续产量 (maximum sustainable yield, MSY) 和可捕量等的评估,难以进行有效管理。另外,巴基斯坦缺乏渔业调查数据,已有MSY评估研究多基于剩余产量模型[6, 9-10],单位捕捞努力量渔获量 (catch per unit effort, CPUE) 数据未经过标准化处理,评估结果存在一定的不确定性[11]。对于渔业数据缺乏,国际上通常采用基于产量统计数据的评估模型,如depletion-corrected average catch (DCAC)[12]、depletion-based stock reduction analysis (DB-SRA)[13]等,但是以上方法多针对长寿命鱼类[14],对于巴基斯坦近海多数经济鱼类不适用。本研究采用一种基于统计产量数据和种群参数先验信息的Catch-MSY模型[15]对巴基斯坦海洋渔业总可捕量及24个重要经济类群的可捕量和MSY进行评估,以期为巴基斯坦海洋渔业管理和我国在阿拉伯海的远洋渔业提供科学依据。

    用于模型分析的1950—2015年巴基斯坦海洋渔业产量数据由卡拉奇大学 (University of Karachi) 巴方专家Sher Khan Panhwar等提供,经过联合国粮农组织(FAO)统计数据库(http://www.fao.org/statistics/databases/en/) 校准。评估的24个重要经济类群分别为鳀类(Anchovies,英文名称均按照FAO数据库中类群名称,下同)、魣类 (Barracudas)、乌鲳 (Formio niger)、鲳类 (Butterfishes)、鲹类 (Carangids)、军曹鱼 (Rachycentron canadum)、鲯鳅 (Coryphaena hippurus)、石首鱼类 (Croakers)、宝刀鱼 (Chirocentrus dorab)、石斑鱼类 (Groupers)、石鲈类 (Grunts)、长头小沙丁鱼 (Sardinella longiceps)、马鲹 (Caranx hippos)、白带鱼 (Trichiurus lepturus)、银纹笛鲷 (Lutjanus argentimaculatus)、鲻类 (Mullets)、康氏马鲛 (Scomberomorus commerson)、海鳗 (Muraenesox cinereus)、鲷类 (Porgies)、鳐类 (Rays)、海鲇 (Arius thalassinus)、金线鱼类 (Threadfin breams)、舌鳎类 (Tonguefishes)和大甲鲹 (Megalaspis cordyla)。

    评估模型采用一种简化的产量模型Catch-MSY模型,评估过程不需要CPUE数据,而以资源量水平和内禀增长率 (r) 的先验分布代替[15],产量数据和模型参数均可按照类群得到,符合巴基斯坦渔业资源特征以及数据现状。Catch-MSY模型为:

    $${B_t} = {\lambda _0}k\exp ({\nu _t})$$ (1)
    $${B_{t + 1}} = [{B_t} + r{B_t}(1 - {B_t}/k) - {C_t}]\exp ({\nu _t})$$ (2)

    其中Btt年的资源量,k为环境容量,Ctt年的渔获量;假定过程误差符合对数正态分布,因此νt为均值为0,方差为σ2的标准正态分布;λ0为起始资源量水平B1/k

    采用如下伯努利分布作为似然函数:

    $$\begin{gathered} L(\Theta |{C_t}) = 1{\rm{ }}\qquad{\lambda _{01}} \leqslant {B_{n + 1}}/k \leqslant {\lambda _{02}} \\ \qquad\quad \; {\rm{ = 0 }}\qquad{\lambda _{01}} > {B_{n + 1}}/k > {\lambda _{02}} \\ \end{gathered} $$ (3)

    其中Θ为模型中的参数向量,[λ01, λ02] 为最终年份资源量水平的先验分布区间。这样的似然函数可以保证r-k参数组合可以得到种群状态的有效解[5]

    研究采用重要性重抽样 (SIR)[16]方法计算参数的后验分布,每次计算的迭代次数50 000次。利用得到的r-k联合后验分布计算MSY,MSY=0.25rk,采用MSY的90%作为可捕量的设置标准[15]

    模型中参数内禀增长率和资源量水平的先验分布均采用均匀分布形式[14-15]。其中内禀增长率的先验分布根据Fishbase数据库的鱼类恢复力分级法[14,17]确定。资源量水平的先验分布则根据评估对象的开发状态以及产量与数据中最大产量的比值确定[14-15]。如1954年巴基斯坦渔业开发皆为小型木船,石首鱼类捕捞努力量和产量都很低,因此初始年份1954年资源量水平设置为0.5~0.9,石首鱼类的产量在2002年达到顶峰后有所下滑,最终年份2015年资源量水平设置为0.3~0.7。

    由于巴基斯坦渔业种类繁多,难以确定其海域渔业资源综合内禀增长率的先验分布,本研究参考中国南海区综合种群内禀增长率评估结果[14],设置了4种不同的先验分布(0.6~1.5;0.4~1.5;0.6~1.7;0.4~1.7),以评估不同内禀增长率先验分布区间对MSY评估结果的影响。研究中评估类群的产量数据序列以及参数先验分布设置见表1。建模和数据分析都在R语言3.3.1中完成。

    表  1  评估对象产量数据序列及参数先验分布设置
    Table  1.  Catch data series and prior distribution for fish groups in stock assessment
    类群
    fish group
    数据序列
    data series
    先验分布 prior distribution
    内禀增长率
    r
    起始年份资源量水平
    first year B1/k
    最终年份资源量水平
    final year Bn/k
    鳀类 Anchovies 1988—2015 [0.6, 1.5] [0.5, 0.9] [0.01, 0.4]
    魣类 Barracudas 1978—2015 [0.05, 0.5] [0.5, 0.9] [0.3, 0.7]
    乌鲳 Formio niger 1962—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    鲳类 Butterfishes 1985—2015 [0.2, 1] [0.5, 0.9] [0.3, 0.7]
    鲹类 Carangids 1962—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    军曹鱼 Rachycentron canadum 1962—2015 [0.05, 0.5] [0.5, 0.9] [0.3, 0.7]
    鲯鳅 Coryphaena hippurus 1983—2015 [0.6, 1.5] [0.5, 0.9] [0.3, 0.7]
    石首鱼类 Croakers 1954—2015 [0.2, 1] [0.5, 0.9] [0.3, 0.7]
    宝刀鱼 Chirocentrus dorab 1971—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    石斑鱼类 Groupers 1962—2015 [0.05, 0.5] [0.5, 0.9] [0.3, 0.7]
    石鲈类 Grunts 1962—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    长头小沙丁鱼 Sardinella longiceps 1950—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    马鲹 Caranx hippos 1981—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    白带鱼 Trichiurus lepturus 1980—2015 [0.05, 0.5] [0.5, 0.9] [0.01, 0.4]
    银纹笛鲷 Lutjanus argentimaculatus 1962—2015 [0.05, 0.5] [0.5, 0.9] [0.01, 0.4]
    鲻类 Mullets 1970—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    康氏马鲛 Scomberomorus commerson 1950—2015 [0.2, 1] [0.5, 0.9] [0.3, 0.7]
    海鳗 Muraenesox cinereus 1978—2015 [0.05, 0.5] [0.5, 0.9] [0.01, 0.4]
    鲷类 Porgies 1970—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    鳐类 Rays 1962—2015 [0.05, 0.5] [0.5, 0.9] [0.01, 0.4]
    海鲇 Arius thalassinus 1950—2015 [0.05, 0.5] [0.5, 0.9] [0.01, 0.4]
    金线鱼类 Threadfin breams 1985—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    舌鳎类 Tonguefishes 1950—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    大甲鲹 Megalaspis cordyla 1987—2015 [0.2, 1] [0.5, 0.9] [0.3, 0.7]
    注:[ ]. 均匀分布区间 Note: Square brackets represent uniform distributions.
    下载: 导出CSV 
    | 显示表格

    当内禀增长率先验分布服从均匀分布[0.6, 1.5]情况下的模型评估结果见图2。结果显示,巴基斯坦海洋渔业在1992—2002年间产量超过MSY,处于过度捕捞状态,近年来产量有所下降,但是仍在MSY附近 (图2-a)。内禀增长率 (r) 环境容量 (k) 具有明显的负相关关系 (图2-b),ln(r) 和ln(k) 则具有明显的线性关系 (图2-c)。MSY的后验概率密度分布要比参数rk的相对范围窄,即变异系数CV要低 ( 图2-d~f)。

    图  2  内禀增长率参数先验分布为 [0.6, 1.5] 情况下巴基斯坦海洋渔业资源最大可持续产量评估结果
    a. 1950—2015年统计产量与MSY评估结果(实线为几何平均数,虚线为正负2倍标准差,下同;b. 参数rk的先验均匀分布,黑色点为后验组合;c. 为ln(r) 和 ln(k)的线性关系以及MSY的几何平均数 (实线);d~f. 为参数rk以及MSY的后验概率密度
    Figure  2.  Model outputs for fisheries in Pakistani waters with r of [0.6, 1.5]
    a. catch history from 1950 to 2015 with MSY estimation (solid line) ± 2SD (dash line); b. prior uniform distribution of rk, and the black dots are posterior combinations; c. the relationship between ln(r) and ln(k) with geometric mean MSY (solid line) ± 2SD (dash line); d~f. posterior densities of r, k and MSY

    4种不同内禀增长率先验分布评估的巴基斯坦海洋渔业MSY分别为40.88×104 t、40.30×104 t、40.13×104 t和40.80×104 t (表2),相差不大,取平均值40.53×104 t;海洋渔业总可捕量为36.47×104 t。2015年巴基斯坦海洋渔业产量为36.10×104 t,目前未处于过度捕捞状态。

    表  2  4种不同内禀增长率先验分布评估的巴基斯坦海洋渔业最大可持续产量和可捕量评估结果
    Table  2.  Estimated MSY and allowable catch under four prior levels of intrinsic rate of increase for marine fisheries in Pakistani waters
    内禀增长率先验分布
    prior of intrinsic rate of increase
    评估结果 assessment result
    内禀增长率
    r
    最大可持续产量/104 t
    MSY
    可捕量/104 t
    allowable catch
    [0.6, 1.5] 0.976 (CV=0.31) 40.88 (CV=0.023) 36.79
    [0.4, 1.5] 0.637 (CV=0.37) 40.30 (CV=0.027) 36.27
    [0.4, 1.7] 0.654 (CV=0.39) 40.13 (CV=0.026) 36.11
    [0.6, 1.7] 0.905 (CV=0.28) 40.80 (CV=0.022) 36.72
    注:CV. 变异系数 Note: CV. coefficient of variation
    下载: 导出CSV 
    | 显示表格

    巴基斯坦24个重要经济类群的评估结果显示,所有类群的内禀增长率介于0.11~0.85,营养级较高的鱼类如鳐类、银纹笛鲷、海鳗和石斑鱼类内禀增长率较低,在0.2以内 (表3)。有8个类群 (军曹鱼、魣类、石斑鱼类、鲳类、鲯鳅、宝刀鱼、马鲹和白带鱼) 2015年产量超过MSY,处于过度捕捞状态,其中军曹鱼、魣类、石斑鱼类和鲯鳅产量超过MSY 20%以上,过度捕捞状态严重 (图3)。其他16个类群2015年产量小于MSY,未处于过度捕捞状态,其中金线鱼类、鳀类、鳐类和鲹类2015年产量低于MSY的40%以上。24个类群中可捕量超过104 t的类群有鳀类、石首鱼类、石斑鱼类、长头小沙丁鱼、白带鱼、鳐类和海鲇 (表3)。

    表  3  巴基斯坦24个重要经济类群评估结果
    Table  3.  Assessment results of 24 important commercial fish groups in Pakistani waters
    类群
    fish group
    内禀增长率
    r
    最大可持续产量/104 t
    MSY
    可捕量/104t
    allowable catch
    2015年产量/104t
    catch in 2015
    鳀类 Anchovies 0.73 14.37 (0.09) 12.93 6.45
    魣类 Barracudas 0.26 4.46 (0.12) 4.01 6.69
    乌鲳 Formio niger 0.26 2.91 (0.13) 2.62 2.29
    鲳类 Butterfishes 0.54 4.37 (0.08) 3.93 4.93
    鲹类 Carangids 0.28 10.09 (0.06) 9.08 2.61
    军曹鱼 Rachycentron canadum 0.22 2.31 (0.03) 2.08 3.89
    鲯鳅 Coryphaena hippurus 0.85 3.58 (0.07) 3.22 4.61
    石首鱼类 Croakers 0.41 18.77 (0.11) 16.89 16.16
    宝刀鱼 Chirocentrus dorab 0.26 2.68 (0.13) 2.41 2.78
    石斑鱼类 Groupers 0.20 11.46 (0.17) 10.31 15.98
    石鲈类 Grunts 0.30 4.84 (0.09) 4.36 4.01
    长头小沙丁鱼 Sardinella longiceps 0.38 35.98 (0.13) 32.38 27.78
    马鲹 Caranx hippos 0.32 6.70 (0.05) 6.03 7.32
    白带鱼 Trichiurus lepturus 0.26 12.46 (0.06) 11.21 12.66
    银纹笛鲷 Lutjanus argentimaculatus 0.13 1.56 (0.03) 1.40 1.26
    鲻类 Mullets 0.29 9.82 (0.09) 8.84 9.53
    康氏马鲛 Scomberomorus commerson 0.52 9.89 (0.18) 8.90 9.79
    海鳗 Muraenesox cinereus 0.14 3.16 (0.21) 2.84 2.76
    鲷类 Porgies 0.29 2.90 (0.14) 2.61 1.96
    鳐类 Rays 0.11 13.92 (0.03) 12.53 4.72
    海鲇 Arius thalassinus 0.33 21.57 (0.22) 19.41 19.29
    金线鱼类 Threadfin breams 0.35 5.06 (0.08) 4.55 2.89
    舌鳎类 Tonguefishes 0.30 1.42 (0.09) 1.28 1.12
    大甲鲹 Megalaspis cordyla 0.47 3.75 (0.10) 3.38 3.63
    注:括号内为变异系数CV值 Note: Value in brackets are coefficient of variations.
    下载: 导出CSV 
    | 显示表格
    图  3  巴基斯坦24个海洋渔业类群2015年产量与最大可持续产量的比值
    虚线代表比值为1
    Figure  3.  Ratios of catch in 2015 to MSY for 24 important commercial fish groups in Pakistani waters Dash line represent that the ratio is equal to 1.

    渔业的量化管理是保证渔业可持续发展的重要举措[18-20],然而全球仅有不足1%的渔业种类进行过资源评估[21],大多数渔业资源由于缺乏调查数据,难以用传统的方法评估MSY和可捕量[22]。基于产量的评估模型可以利用统计产量数据和渔业种群的生活史特征参数等对MSY进行评估,不再需要调查数据,其中比较常用的几种如DCAC、DBSRA和only reliable catch stocks (ORCS) 等[23],但是这些方法对评估对象的种群生活史参数等有所限制[14]。巴基斯坦产量多以类群进行统计,并且许多种类的自然死亡系数大于0.2[1,6,9],无法利用以上方法进行评估。本研究采用的Catch-MSY模型对渔业种群生物学特征没有选择性,模型所需的统计产量数据允许按类群给出,更适合用于当前巴基斯坦海域的资源评估。

    内禀增长率先验分布的设置是模型评估过程中不确定性的重要来源,本研究采用Fishbase数据库对每个类群的内禀增长率进行分级,主要依靠von Bertallanffy生长方程[24]中的生长参数、性成熟年龄、最大年龄以及繁殖力等参数确定[17]。但是巴基斯坦海洋渔业资源综合种群内禀增长率缺乏相关研究,因此本研究以与巴基斯坦海域相近纬度的南海近海评估结果为依据[14],设置了4种不同的内禀增长率先验分布评估其不确定性。结果显示,4种不同的内禀增长率先验分布评估的MSY分别为40.88×104 t、40.30×104 t、40.13×104 t和40.80×104 t,变动幅度在2%以内,说明内禀增长率的先验分布对MSY的评估结果影响不大,这一结果与南海的评估结果一致[14]。因此,利用Catch-MSY模型进行评估时,评估结果更受制于统计产量数据的精确度,本研究采用的统计产量数据来自巴基斯坦官方统计并在FAO数据库中逐条验证。

    模型评估的巴基斯坦海洋渔业MSY为40.53×104 t,总可捕量为36.47×104 t,而2015年巴基斯坦海洋渔业产量为36.10×104 t,虽然当前的捕捞产量未超过MSY,但是已经达到总可捕量,需要进行适当的管理,控制捕捞产量。24个重要经济类群的评估结果显示,8个类群处于过度捕捞状态,但是16个未处于过度捕捞状态的类群中,康氏马鲛、鲻类和大甲鲹2015年的产量超过了可捕量,需要进行管理控制产量。另外,金线鱼类、鳀类、鳐类和鲹类等4个类群虽然2015年产量低于MSY的40%以上,但是从产量时间序列来看(图4),这4个类群捕捞产量均在过去一段时期内超过MSY,过度捕捞后资源崩溃导致目前的产量处于低位。因此,这4个类群虽然产量低于MSY,但是资源仍未恢复,需要实行禁捕等措施,待资源恢复。该模型对与巴基斯坦海域相近纬度的南海近海海域渔业资源进行了评估,得到了总可捕量和11个重要类群可捕量的评估结果[14]。结果显示,南海近海石斑鱼类等低内禀增长率、高营养级鱼类的过度捕捞比较严重,而巴基斯坦中较低内禀增长率的石斑鱼类、军曹鱼和魣类等也有类似结果。低内禀增长率的渔业种群对应着其低繁殖力和低恢复力,一旦种群崩溃便很难恢复[25]

    图  4  巴基斯坦金线鱼类、鳐类、鳀类和鲹类等4个海洋渔业类群的统计产量与最大可持续产量评估值
    实线为平均值,虚线为95%置信区间
    Figure  4.  Statistical catches and estimated MSYs of threadfin breams, rays, anchovies and carangids fisheries in Pakistani waters
    Solid lines are average values and dash lines are 95% confidence intervals.

    评估结果显示,当前巴基斯坦海洋渔业资源已处于充分开发状态,多数经济类群已处于过度捕捞状态或者崩溃后的恢复状态,开发潜力较小。与我国南海深海区类似,阿拉伯海现有丰富的渔业资源,如鸢乌贼、金枪鱼和鲣等[26-27],尤其是鸢乌贼最具开发潜力[28]。我国多次对阿拉伯海鸢乌贼资源进行了探捕[28],南海区也已形成了以灯光罩网渔船为主的捕捞作业船队[29];因此,以阿拉伯海的鸢乌贼渔业作为切入点,充分利用中巴全天候战略合作伙伴关系加强合作,与巴基斯坦共同开发印度洋公海渔业资源,对“一带一路”倡议的推动,中国和巴基斯坦渔业的可持续发展,以及中巴双方渔业领域科学家的交流与合作均有重要意义。

  • 图  1   1950—2016年巴基斯坦海洋渔业统计产量

    Figure  1.   Statistic data of marine fisheries catches in Pakistani waters from 1950 to 2016

    图  2   内禀增长率参数先验分布为 [0.6, 1.5] 情况下巴基斯坦海洋渔业资源最大可持续产量评估结果

    a. 1950—2015年统计产量与MSY评估结果(实线为几何平均数,虚线为正负2倍标准差,下同;b. 参数rk的先验均匀分布,黑色点为后验组合;c. 为ln(r) 和 ln(k)的线性关系以及MSY的几何平均数 (实线);d~f. 为参数rk以及MSY的后验概率密度

    Figure  2.   Model outputs for fisheries in Pakistani waters with r of [0.6, 1.5]

    a. catch history from 1950 to 2015 with MSY estimation (solid line) ± 2SD (dash line); b. prior uniform distribution of rk, and the black dots are posterior combinations; c. the relationship between ln(r) and ln(k) with geometric mean MSY (solid line) ± 2SD (dash line); d~f. posterior densities of r, k and MSY

    图  3   巴基斯坦24个海洋渔业类群2015年产量与最大可持续产量的比值

    虚线代表比值为1

    Figure  3.   Ratios of catch in 2015 to MSY for 24 important commercial fish groups in Pakistani waters Dash line represent that the ratio is equal to 1.

    图  4   巴基斯坦金线鱼类、鳐类、鳀类和鲹类等4个海洋渔业类群的统计产量与最大可持续产量评估值

    实线为平均值,虚线为95%置信区间

    Figure  4.   Statistical catches and estimated MSYs of threadfin breams, rays, anchovies and carangids fisheries in Pakistani waters

    Solid lines are average values and dash lines are 95% confidence intervals.

    表  1   评估对象产量数据序列及参数先验分布设置

    Table  1   Catch data series and prior distribution for fish groups in stock assessment

    类群
    fish group
    数据序列
    data series
    先验分布 prior distribution
    内禀增长率
    r
    起始年份资源量水平
    first year B1/k
    最终年份资源量水平
    final year Bn/k
    鳀类 Anchovies 1988—2015 [0.6, 1.5] [0.5, 0.9] [0.01, 0.4]
    魣类 Barracudas 1978—2015 [0.05, 0.5] [0.5, 0.9] [0.3, 0.7]
    乌鲳 Formio niger 1962—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    鲳类 Butterfishes 1985—2015 [0.2, 1] [0.5, 0.9] [0.3, 0.7]
    鲹类 Carangids 1962—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    军曹鱼 Rachycentron canadum 1962—2015 [0.05, 0.5] [0.5, 0.9] [0.3, 0.7]
    鲯鳅 Coryphaena hippurus 1983—2015 [0.6, 1.5] [0.5, 0.9] [0.3, 0.7]
    石首鱼类 Croakers 1954—2015 [0.2, 1] [0.5, 0.9] [0.3, 0.7]
    宝刀鱼 Chirocentrus dorab 1971—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    石斑鱼类 Groupers 1962—2015 [0.05, 0.5] [0.5, 0.9] [0.3, 0.7]
    石鲈类 Grunts 1962—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    长头小沙丁鱼 Sardinella longiceps 1950—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    马鲹 Caranx hippos 1981—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    白带鱼 Trichiurus lepturus 1980—2015 [0.05, 0.5] [0.5, 0.9] [0.01, 0.4]
    银纹笛鲷 Lutjanus argentimaculatus 1962—2015 [0.05, 0.5] [0.5, 0.9] [0.01, 0.4]
    鲻类 Mullets 1970—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    康氏马鲛 Scomberomorus commerson 1950—2015 [0.2, 1] [0.5, 0.9] [0.3, 0.7]
    海鳗 Muraenesox cinereus 1978—2015 [0.05, 0.5] [0.5, 0.9] [0.01, 0.4]
    鲷类 Porgies 1970—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    鳐类 Rays 1962—2015 [0.05, 0.5] [0.5, 0.9] [0.01, 0.4]
    海鲇 Arius thalassinus 1950—2015 [0.05, 0.5] [0.5, 0.9] [0.01, 0.4]
    金线鱼类 Threadfin breams 1985—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    舌鳎类 Tonguefishes 1950—2015 [0.2, 1] [0.5, 0.9] [0.01, 0.4]
    大甲鲹 Megalaspis cordyla 1987—2015 [0.2, 1] [0.5, 0.9] [0.3, 0.7]
    注:[ ]. 均匀分布区间 Note: Square brackets represent uniform distributions.
    下载: 导出CSV

    表  2   4种不同内禀增长率先验分布评估的巴基斯坦海洋渔业最大可持续产量和可捕量评估结果

    Table  2   Estimated MSY and allowable catch under four prior levels of intrinsic rate of increase for marine fisheries in Pakistani waters

    内禀增长率先验分布
    prior of intrinsic rate of increase
    评估结果 assessment result
    内禀增长率
    r
    最大可持续产量/104 t
    MSY
    可捕量/104 t
    allowable catch
    [0.6, 1.5] 0.976 (CV=0.31) 40.88 (CV=0.023) 36.79
    [0.4, 1.5] 0.637 (CV=0.37) 40.30 (CV=0.027) 36.27
    [0.4, 1.7] 0.654 (CV=0.39) 40.13 (CV=0.026) 36.11
    [0.6, 1.7] 0.905 (CV=0.28) 40.80 (CV=0.022) 36.72
    注:CV. 变异系数 Note: CV. coefficient of variation
    下载: 导出CSV

    表  3   巴基斯坦24个重要经济类群评估结果

    Table  3   Assessment results of 24 important commercial fish groups in Pakistani waters

    类群
    fish group
    内禀增长率
    r
    最大可持续产量/104 t
    MSY
    可捕量/104t
    allowable catch
    2015年产量/104t
    catch in 2015
    鳀类 Anchovies 0.73 14.37 (0.09) 12.93 6.45
    魣类 Barracudas 0.26 4.46 (0.12) 4.01 6.69
    乌鲳 Formio niger 0.26 2.91 (0.13) 2.62 2.29
    鲳类 Butterfishes 0.54 4.37 (0.08) 3.93 4.93
    鲹类 Carangids 0.28 10.09 (0.06) 9.08 2.61
    军曹鱼 Rachycentron canadum 0.22 2.31 (0.03) 2.08 3.89
    鲯鳅 Coryphaena hippurus 0.85 3.58 (0.07) 3.22 4.61
    石首鱼类 Croakers 0.41 18.77 (0.11) 16.89 16.16
    宝刀鱼 Chirocentrus dorab 0.26 2.68 (0.13) 2.41 2.78
    石斑鱼类 Groupers 0.20 11.46 (0.17) 10.31 15.98
    石鲈类 Grunts 0.30 4.84 (0.09) 4.36 4.01
    长头小沙丁鱼 Sardinella longiceps 0.38 35.98 (0.13) 32.38 27.78
    马鲹 Caranx hippos 0.32 6.70 (0.05) 6.03 7.32
    白带鱼 Trichiurus lepturus 0.26 12.46 (0.06) 11.21 12.66
    银纹笛鲷 Lutjanus argentimaculatus 0.13 1.56 (0.03) 1.40 1.26
    鲻类 Mullets 0.29 9.82 (0.09) 8.84 9.53
    康氏马鲛 Scomberomorus commerson 0.52 9.89 (0.18) 8.90 9.79
    海鳗 Muraenesox cinereus 0.14 3.16 (0.21) 2.84 2.76
    鲷类 Porgies 0.29 2.90 (0.14) 2.61 1.96
    鳐类 Rays 0.11 13.92 (0.03) 12.53 4.72
    海鲇 Arius thalassinus 0.33 21.57 (0.22) 19.41 19.29
    金线鱼类 Threadfin breams 0.35 5.06 (0.08) 4.55 2.89
    舌鳎类 Tonguefishes 0.30 1.42 (0.09) 1.28 1.12
    大甲鲹 Megalaspis cordyla 0.47 3.75 (0.10) 3.38 3.63
    注:括号内为变异系数CV值 Note: Value in brackets are coefficient of variations.
    下载: 导出CSV
  • [1]

    KALHORO M A, LIU Q, VALINASSAB T, et al. Population dynamics of greater lizardfish, Saurida tumbil from Pakistani waters[J]. Pak J Zool, 2015, 47(4): 921-931.

    [2]

    FAO. Fishery and Aquaculture Country Profiles. Pakistan (2009). Country Profile Fact Sheets.[EB/OL]. (2017-12-01).http://www.fao.org/fishery/facp/PAK/en.

    [3]

    PSOMADAKIS P N, OSMANY H B, MOAZZAM M. Field identification guide to the living marine resources of Pakistan[R]. Rome: FAO, 2015.

    [4]

    FANNING L P, KHAN M W, KIDWAI S, et al. Surveys of the offshore fisheries resources of Pakistan-2010[R]. Rome: FAO Fisheries and Aquaculture Circular, 2011: 1065.

    [5] 陈峰, 蒋日进, 朱文斌, 等. 巴基斯坦海洋渔业现状与合作开发对策分析[J]. 海洋开发与管理, 2016, 33(12): 13-18. doi: 10.3969/j.issn.1005-9857.2016.12.003
    [6] KALHORO M A. Population dynamics of nemipteridae and synodontidae fish families from Pakistani waters[D]. 青岛: 中国海洋大学, 2014: 17-43.
    [7]

    KALHORO M A, LIU Q, WARYANI B, et al. Growth and mortality of brushtooth lizardfish, Saurida undosquamis, from Pakistani waters[J]. Pak J Zool, 2014, 46(1): 139-151.

    [8]

    MEMON K H, LIU Q, KALHORO M A, et al. Maximum sustainable yield estimates of barramundi Lates calcarifer fishery from Pakistani waters[J]. Indian J Geo-Mar Sci, 2015, 44(6): 825-832.

    [9]

    KALHORO M A, LIU Q, MEMON K H, et al. Maximum sustainable yield of greater lizardfish Saurida tumbil fishery in Pakistan using the CEDA and ASPIC packages[J]. Acta Oceanol Sin, 2015, 34(2): 68-73. doi: 10.1007/s13131-014-0463-0

    [10]

    PANHWAR S K, LIU Q, AMIR S A, et al. Performance comparison between logistic and generalized surplus-production models applied to the Sillago sihama fishery in Pakistan[J]. J Ocean U China, 2012, 11(3): 401-407. doi: 10.1007/s11802-012-1930-x

    [11] 张魁, 刘群, 廖宝超, 等. 渔业数据失真对两种非平衡剩余产量模型评估结果的影响比较[J]. 水产学报, 2018, 42(9): 1378-1389.
    [12]

    MACCALL A D. Depletion-corrected average catch: a simple formula for estimating sustainable yields in data-poor situations[J]. ICES J Mar Sci, 2009, 66(9): 2267-2271.

    [13]

    DICK E J, MACCALL A D. Depletion-based stock reduction analysis: a catch-based method for determining sustainable yields for data-poor fish stocks[J]. Fish Res, 2011, 110(2): 331-341. doi: 10.1016/j.fishres.2011.05.007

    [14] 张魁, 廖宝超, 许友伟, 等. 基于渔业统计数据的南海区渔业资源可捕量评估[J]. 海洋学报, 2017, 39(8): 25-33. doi: 10.3969/j.issn.0253-4193.2017.08.003
    [15]

    MARTELL S, FROESE R. A simple method for estimating msy from catch and resilience[J]. Fish Fish, 2013, 14(4): 504-514. doi: 10.1111/j.1467-2979.2012.00485.x

    [16]

    HADDON M. Modelling and quantitative methods in fisheries (Second Edition) [M]. London: Chapman and Hall , 2001: 76-98.

    [17]

    FROESE R, PAULY D. Estimation of life history key facts. In: FishBase 2000: concepts, design and data sources. (eds Froese R and Pauly D) [R]. Philippines: ICLARM, 2000.

    [18] 慕永通. 渔业管理: 以基于权利的管理为中心[M]. 青岛: 中国海洋大学出版社, 2006: 1-25.
    [19] 易传剑, 周梅芳, 俞存根. 政府规制理论在我国近海渔业管理中应用的探讨[J]. 水产学报, 2012, 36(5): 787-793.
    [20] 郭文路, 黄硕琳. 总可捕量制度不同实施方式的比较分析[J]. 海洋湖沼通报, 2001(4): 61-69. doi: 10.3969/j.issn.1003-6482.2001.04.012
    [21]

    RICARD D, MINTO C, JENSEN O P, et al. Examining the knowledge base and status of commercially exploited marine species with the RAM Legacy Stock Assessment Database[J]. Fish Fish, 2012, 13(4): 380-398. doi: 10.1111/faf.2012.13.issue-4

    [22]

    CARRUTHERS T R, PUNT A E, WALTERS C J, et al. Evaluating methods for setting catch limits in data-limited fisheries[J]. Fish Res, 2014, 153(5): 48-68.

    [23]

    BERKSON J, THORSON J T. The determination of data-poor catch limits in the United States: is there a better way?[J]. ICES J Mar Sci, 2015, 72(1): 237-242. doi: 10.1093/icesjms/fsu085

    [24] 詹秉义. 渔业资源评估[M]. 北京: 中国农业出版社, 1995: 25-31.
    [25]

    ALLISON G. The influence of species diversity and stress intensity on community resistance and resilience[J]. Ecol Monogr, 2004, 74(1): 117-134. doi: 10.1890/02-0681

    [26] 韩保平. 巴基斯坦海洋渔业概况[J]. 海洋渔业, 1999(2): 90-91.
    [27] 李善勋. 巴基斯坦的海洋渔业[J]. 中国水产, 1983(9): 28-29.
    [28] 陈新军, 钱卫国, 田思泉. 阿拉伯海北部公海海域鸢乌贼资源密度及其分布[J]. 海洋科学进展, 2006, 24(3): 360-364. doi: 10.3969/j.issn.1671-6647.2006.03.014
    [29] 张鹏, 张俊, 李渊, 等. 秋季南海中南部海域的一次灯光罩网探捕调查[J]. 南方水产科学, 2016, 12(2): 67-74. doi: 10.3969/j.issn.2095-0780.2016.02.010
  • 期刊类型引用(6)

    1. 张静,戴佳玥,来新昊,刘旭祥,张浩,王学锋,汤保贵. 卵形鲳鲹应对流速胁迫的代谢组学分析. 海洋学报. 2023(05): 53-63 . 百度学术
    2. 段鹏飞,田永胜,李振通,李子奇,陈帅,黎琳琳,王心怡,王林娜,刘阳,李文升,王晓梅,李波. 棕点石斑鱼(♀)×蓝身大斑石斑鱼(♂)杂交后代与棕点石斑鱼低氧耐受能力初步研究. 中国水产科学. 2022(02): 220-233 . 百度学术
    3. 逯云召,于燕光,薄其康,马超,宓慧菁,孙晓旺. 大泷六线鱼幼鱼的摄食节律研究. 渔业现代化. 2021(02): 35-39 . 百度学术
    4. 李志成,江飚,钟志鸿,李诗钰,何润真,唐嘉嘉,李安兴. 硫酸铜治疗卵形鲳鲹淀粉卵涡鞭虫病的研究. 南方水产科学. 2021(03): 108-114 . 本站查看
    5. 韩明洋,周胜杰,杨蕊,胡静,马振华. 温度胁迫下卵形鲳鲹仔鱼骨骼组织病理及分子表征. 南方农业学报. 2021(11): 3147-3156 . 百度学术
    6. 黄小林,戴超,虞为,杨洁,杨育凯,李涛,林黑着,黄忠,孙莘溢,舒琥. 丁香酚对卵形鲳鲹幼鱼的麻醉效果. 广东海洋大学学报. 2020(04): 124-131 . 百度学术

    其他类型引用(3)

图(4)  /  表(3)
计量
  • 文章访问数:  5345
  • HTML全文浏览量:  2174
  • PDF下载量:  99
  • 被引次数: 9
出版历程
  • 收稿日期:  2019-03-19
  • 修回日期:  2019-04-09
  • 录用日期:  2019-05-24
  • 网络出版日期:  2019-05-30
  • 刊出日期:  2019-10-04

目录

/

返回文章
返回