留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

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

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

引用本文:
Citation:

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

    作者简介: 张 魁 (1987—),男,博士,助理研究员,从事渔业资源评估研究。E-mail: nedvedkui@163.com;
    通讯作者: 陈作志, zzchen2000@163.com
  • 中图分类号: S 937.3

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

    Corresponding author: Zuozhi CHEN, zzchen2000@163.com
  • CLC number: S 937.3

  • 摘要: 基于巴基斯坦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,处于过度捕捞状态。当前巴基斯坦海洋渔业资源已处于充分开发状态,多数经济类群已处于过度捕捞状态或者崩溃后的恢复状态,开发潜力较小,建议开发阿拉伯海鸢乌贼等大洋性渔业资源。
  • 图 1  1950—2016年巴基斯坦海洋渔业统计产量

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

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

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

    图 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  巴基斯坦24个海洋渔业类群2015年产量与最大可持续产量的比值

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

    表 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 results
    内禀增长率
    r
    最大可持续产量
    MSY/104 t
    可捕量/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
    最大可持续产量
    MSY/104 t
    可捕量/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. In: FAO Fisheries and Aquaculture Department[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]. FAO Species Identification Guide for Fishery Purposes. FAO, 2015.
    [4] FANNING L P, KHAN M W, KIDWAI S, et al. Surveys of the offshore fisheries resources of Pakistan-2010[R]. 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
  • [1] 刘莉莉 万荣 王熙杰 王四杰 王云中 . 基于系统动力学模型的渔业资源增殖放流效应分析. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2012.01.003
    [2] 李斌陈国宝1郭禹陈作志张俊1王东旭 . 南海中部海域渔业资源时空分布和资源量的水声学评估. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2016.04.004
    [3] 张俊张鹏陈作志陈国宝张魁许友伟孙铭帅 . 南海外海鲹科鱼类资源量及其分布. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2016.04.005
    [4] 晏磊谭永光杨炳忠张鹏李杰杨吝 . 基于张网渔业休渔前后的黄茅海河口渔业资源群落比较. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2016.06.001
    [5] 杨吝张旭丰谭永光张鹏 . 南海北部灯光罩网渔获组成及其对渔业资源的影响. 南方水产科学, doi: 10.3969/j.issn.1673-2227.2009.04.008
    [6] 袁梦汤勇徐姗楠陈作志杨玉滔江艳娥 . 珠江口南沙海域秋季渔业资源群落结构特征. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2017.02.003
    [7] 蔡研聪徐姗楠陈作志许友伟江艳娥杨长平 . 南海北部近海渔业资源群落结构及其多样性现状. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2018.02.002
    [8] 王跃中袁蔚文 . 南海北部底拖网渔业资源的数量变动. 南方水产科学,
    [9] 曾雷陈国宝于杰 . 南澳岛海域渔业资源声学评估与空间分布. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2018.02.004
    [10] 徐海龙谷德贤乔秀亭曹丹丹 . 渤海湾主要渔业资源长度与体质量关系分析. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2014.01.009
    [11] 龙鑫玲陈丕茂袁华荣冯雪余景舒黎明黎小国陈文静陈梓聪 . 深圳大鹏半岛沿岸海域渔业资源承载力评价. 南方水产科学, doi: 10.12131/20190024
    [12] 李健生程家骅 . 长江口水域主要渔业生物资源状况的分析. 南方水产科学,
    [13] 江艳娥林昭进黄梓荣 . 南海北部大陆架渔业生物多样性研究. 南方水产科学, doi: 10.3969/j.issn.1673-2227.2009.05.006
    [14] 蔡研聪陈作志徐姗楠张魁 . 北部湾二长棘犁齿鲷的时空分布特征. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2017.04.001
    [15] 杨吝张旭丰谭永光张鹏 . 南海北部灯光围网渔获组成分析. 南方水产科学, doi: 10.3969/j.issn.1673-2227.2009.06.012
    [16] 苏莹佳陈国宝周艳波马胜伟吴洽儿 . 2015—2017年南海海域伏季休渔制度实施效果评价. 南方水产科学, doi: 10.12131/20180149
    [17] 陈丕茂 . 渔业资源增殖放流效果评估方法的研究. 南方水产科学,
    [18] 陈丕茂袁华荣贾晓平秦传新蔡文贵余景舒黎明黎小国周艳波 . 大亚湾杨梅坑人工鱼礁区渔业资源变动初步研究. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2013.05.016
    [19] 李寇军邱永松王跃中 . 自然环境变动对北部湾渔业资源的影响. 南方水产科学,
    [20] 牛明香赵宪勇 . 卫星遥感和GIS技术在海洋渔业资源研究中的应用. 南方水产科学,
  • 加载中
图(4)表(3)
计量
  • 文章访问数:  556
  • HTML全文浏览量:  54
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-20
  • 录用日期:  2019-04-10
  • 网络出版日期:  2019-05-31

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

    作者简介:张 魁 (1987—),男,博士,助理研究员,从事渔业资源评估研究。E-mail: nedvedkui@163.com
    通讯作者: 陈作志, zzchen2000@163.com
  • 1. 中国水产科学研究院南海水产研究所,农业农村部外海渔业开发重点实验室,广东 广州 510300
  • 2. 上海海洋大学海洋科学学院,上海 201306
  • 3. 卡拉奇大学海洋生物卓越中心,巴基斯坦 卡拉奇 75270

摘要: 基于巴基斯坦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,处于过度捕捞状态。当前巴基斯坦海洋渔业资源已处于充分开发状态,多数经济类群已处于过度捕捞状态或者崩溃后的恢复状态,开发潜力较小,建议开发阿拉伯海鸢乌贼等大洋性渔业资源。

English Abstract

  • 巴基斯坦海域位于阿拉伯海北部,与伊朗和印度所属海域相邻,拥有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})$

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

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

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

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

      其中Θ为模型中的参数向量,[λ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中完成。

      类群
      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.

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

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

    • 当内禀增长率先验分布服从均匀分布[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] 情况下巴基斯坦海洋渔业资源最大可持续产量评估结果

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

      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,目前未处于过度捕捞状态。

      内禀增长率先验分布
      prior of intrinsic rate of increase
      评估结果 assessment results
      内禀增长率
      r
      最大可持续产量
      MSY/104 t
      可捕量/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

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

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

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

      类群
      fish group
      内禀增长率
      r
      最大可持续产量
      MSY/104 t
      可捕量/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.

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

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

      图  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  巴基斯坦24个海洋渔业类群2015年产量与最大可持续产量的比值

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

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

参考文献 (29)

目录

    /

    返回文章
    返回