2. 大洋渔业资源可持续开发教育部重点实验室, 上海 201306;
3. 国家远洋渔业工程技术研究中心, 上海 201306;
4. 农业农村部大洋渔业开发重点实验室, 上海 201306
阿根廷滑柔鱼(Illex angentinus)是大洋性浅海种[1-2], 主要分布在西南大西洋的35°S~52°S海域, 是目前世界头足类中最为重要的种类之一[3]。马尔维纳斯群岛(Malvinas Islands)海域是阿根廷滑柔鱼的重要作业渔场之一, 该海域阿根廷滑柔鱼平均年渔获量为11.9万t, 在高产年份, 马尔维纳斯群岛为世界提供了鱿鱼总量的10%[4]。阿根廷滑柔鱼资源极易受到海洋环境和气候变化的影响, 其资源丰度会出现剧烈的年间变化[5]。已有研究表明, 海表面温度(sea surface temperature, SST)[6]、海平面高度(sea surface height, SSH)[7]、叶绿素a质量浓度(chlorophy a, Chl.a)[8]等环境因子均会影响阿根廷滑柔鱼资源量。此外, 南极涛动(Antarctic oscillation, AAO)[9]等气候因子也影响着阿根廷滑柔鱼资源。国内外学者对阿根廷滑柔鱼渔场时空分布[10]、中心渔场和资源丰度预报模型[11-12]等都做了初步研究, 但研究区域多集中在阿根廷专属经济区和公海区域, 王越等[13]研究了马尔维纳斯群岛海域阿根廷滑柔鱼资源丰度与海表面温度之间的关系, 但目前未见有该海域阿根廷滑柔鱼资源丰度年度灾变预测的研究。
灰色系统理论由我国著名学者邓聚龙教授于1982年创立[14], 包括了灰色预测、灰色关联等研究方法, 已在渔业科学领域进行了很好的应用[15], 如资源丰度预测[16]、渔业灾变年份预测[17]、渔场空间分布[18]等。也有学者[19]应用灰色系统理论对太平洋褶柔鱼(Todarodes pacificus)秋生群体丰歉年进行灾变预测, 取得了较好的预测效果。因此, 本文采用灰色系统理论中的灰色GM(1, 1)模型, 对马尔维纳斯群岛海域阿根廷滑柔鱼丰歉年进行科学预测, 为该海域阿根廷滑柔鱼渔业资源管理与可持续开发提供参考。
1 材料与方法 1.1 研究区域和数据来源阿根廷滑柔鱼的渔业生产数据来自于渔业统计年报[20-23], 统计时间为1995—2019年, 统计信息为每年的渔获量和作业船数, 以单船年捕捞产量(CPUE, t/船)表征阿根廷滑柔鱼资源丰度。其作业渔场为47°S~57°S, 52°W~64°W海域(图 1)[20]。
环境数据为南极海冰覆盖范围, 来自美国国家冰雪数据中心(https://nsidc.org/data/seaice_index/), 时间尺度为天, 单位为百万平方千米。
1.2 分析方法利用GM(1, 1)模型对马尔维纳斯群岛海域阿根廷滑柔鱼资源丰度进行灾变预测[24], 取马尔维纳斯群岛海域阿根廷滑柔鱼CPUE作为原始数列, 记为
取序列的第3分位数和中位数分别为灾变阈值的上限ξ1和下限ξ2, 将原数列X中大于ξ1或小于ξ2的数认为是异常值, 将其对应的时刻组成新的灾变日期序列, 即
再将Q(0)日期序列进行累加, 建立GM(1, 1)模型进行灾变预测。灰色预测模型GM(1, 1)的计算方法参考文献[15]。其级比λ(0)(k)为
其级比偏差为
模型的验证采用残差检验、后验差检验和级比偏差值检验[25], 残差越小, 模型的精度越高, 方差比C越小越好,小误差概率P越大越好。模型精度检验等级见表 1[14]。其中当发展系数|a|≤0.3时, 模型可用于中长期预测; 当0.3<|a|≤0.5时可用于短期预测。若所有的ρ(k)均小于0.2, 则认为模型可行, 其值越小越好; 若所有的ρ(k)均小于0.1, 则认为模型达到较高要求[26]。
根据统计分析显示(图 2), 1995—2019年25年间, 马尔维纳斯群岛海域阿根廷滑柔鱼产量波动较大且有着明显的年间变化, 其年平均产量为99404.28t。2009年阿根廷滑柔鱼产量最低, 仅为3t。此后几年产量开始上升, 到2015年达到最高, 为332863t, 2016年产量急剧下降, 只有2 297t, 之后几年的产量都处于低水平。
1995至2019年期间, 阿根廷滑柔鱼CPUE年平均为1118.807t/船, CPUE变化趋势与产量变化趋势基本一致(图 2), 2009年的CPUE值最小, 仅为0.143t/船。此后几年CPUE开始上升, 2015年CPUE达到最大值, 为3140.217t/船, 2016年CPUE急剧下降, 只有22.086t/船, 2016年以后其CPUE都处于较低水平。由图 2可知, 2007、2008年的CPUE出现异常, 虽然年产量较低, CPUE却处于较高水平, 原因可能与当年作业渔船的减少有关。
2.2 模型预测构建根据灾变值上下限的计算结果(上限灾变值为1615t/船, 下限灾变值为784t/船), 确定丰年年份为1999、2000、2007、2008、2014和2015年, 其对应的时间序列为5、6、13、14、20和21;歉年年份为1995、1996、2002、2004、2005、2009、2010、2016、2017、2018和2019年, 对应的时间序列为1、2、8、10、11、15、16、22、23、24和25。为此, 分别构建丰年和歉年的预测模型:
丰年预测的GM(1, 1)模型为X(t)=33.689exp(0.234t)-28.689
歉年预测的GM(1, 1)模型为X(t)=57.527exp(0.149t)-56.527
灾变模型预测值与真实值之间的相对误差见表 2, 其中丰年灾变模型的相对误差为1.200%~47.924%, 平均相对误差为16.212%;歉年灾变模型的相对误差为0.731%~282.915%, 平均相对误差为35.656%(表 2)。
模型的级比偏差计算值如表 3所示, 丰年灾变模型平均级比偏差为0.193, 歉年灾变模型平均级比偏差为0.199, 均小于0.2, 因此上述构建的灰色预测达到统计要求。其灰色预测模型的有关参数如表 4所示, 丰年灾变模型a=-0.234, 歉年灾变模型a=-0.149, 因此均可用于中长期预测。同时, 其方差比均小于0.35, 小概率误差P均为1.000, 均大于0.95(表 4), 构建模型的精度达到Ⅰ级。
根据建立的灾变预测模型, 对其丰年和歉年进行预测。最后一个丰年为2015年, 下3次出现超过灾变点阈值对应的序号分别为28.57、36.10和45.61(表 4), 因此未来其资源丰度出现丰年的年份为2021、2028和2038年。
同样, 最后一个歉年为2019年, 下3次出现超过灾变点阈值对应的序号分别为33.92、39.36和45.68(表 4), 可得未来其资源丰度出现歉年的年份为2024、2029和2036年。
3 讨论与分析 3.1 阿根廷滑柔鱼资源丰度与气候环境变化的关系阿根廷滑柔鱼具有生命周期短、生长快、终生一次繁殖等特点, 其每一代的资源量取决于上一代亲体产生的补充量[27]。同时, 其资源丰度变化受到多个方面因素的影响, 陈新军等[28]根据2002年1—6月西南大西洋阿根廷滑柔鱼生产数据, 对阿根廷滑柔鱼渔场、CPUE分布及其与表温的关系进行了分析, 认为2004年低产的重要原因是因为作业海域海表温异常偏高, 比往年高出3~4℃。NIGMATULLIN等[29]认为2003—2004年阿根廷滑柔鱼低产的原因可能是2003年10月至12月, 马尔维纳斯海流的加剧导致巴西海流向东偏移, 大部分鱿鱼幼体向东移动, 幼体死亡率大大提升, 同时研究还认为前几年高强度捕捞也是导致其低产的原因。叶绿素a质量浓度一般用于估算海洋生产力, 可用于表征浮游植物量的多少, 通过郑丽丽等[30]的研究发现, 2000年与2007年的叶绿素a质量浓度较高, 有利于形成较好的阿根廷滑柔鱼渔场。
南极海冰覆盖范围指南极地区海冰浓度超过15%的区域面积, 可用于作为指示南极海冰变化的一种大尺度气候变化指标, 它也会对海洋环境和海洋生物产生直接或间接影响。刘赫威等[31]通过研究西南大西洋阿根廷滑柔鱼CPUE与海冰覆盖范围之间的关系, 认为南极海冰的覆盖范围会影响阿根廷滑柔鱼渔场内不同水层的水温, 进而影响其资源丰度和渔场的分布。根据美国国家冰雪数据中心(https://nsidc.org/data/seaice_index/)提供的南极海冰覆盖范围数据, 并结合阿根廷滑柔鱼CPUE丰年和歉年的年份分析(图 3), 发现海冰覆盖范围下降时, 容易出现低产年份, 而海冰覆盖范围上升时, 易出现高产年份。同时在海冰覆盖范围明显波动的年份(2002—2003年, 2011—2012年), 其灾变年份类型仍会维持一定时间, 这说明海冰覆盖范围对阿根廷滑柔鱼资源的影响有着一定的滞后性。另外, 阿根廷滑柔鱼资源丰度与大尺度的大气变化存在一定关系, CHANG等[9]发现南极涛动(Atarctic oscillation, AAO)通过影响其他生物或非生物因素, 在影响阿根廷滑柔鱼资源丰度上具有滞后性。张忠等[32]研究发现, 南极海冰的变化会影响阿根廷滑柔鱼栖息地垂直水温的变化, 导致其栖息地向阿根廷近海移动, 马尔维纳斯群岛海域内阿根廷滑柔鱼资源量下降。此外, 南极海冰的变化会引起海水含盐量变化, 间接影响阿根廷滑柔鱼的资源分布。
通过对马尔维纳斯群岛海域1995—2019年阿根廷滑柔鱼CPUE原始数据的分析, 并运用灰色灾变模型较好地预测其丰歉年, 建立的丰年和歉年模型的平均相对误差分别为16.212%和35.656%, 级比值检验与后验差比检验结果均显示模型精度良好, 可用于中长期预测。这一结果可为渔业生产提供科学指导。据有关资料报道, 2021年马尔维纳斯群岛海域阿根廷滑柔鱼的捕捞产量达到17.2万t, 平均CPUE超过上限灾变值为1615t/船, 是一个阿根廷滑柔鱼资源丰年, 与本模型的预测结果一致, 因此本研究划定的灾变上下限值是可信的。
影响灾变预测模型构建的因素较多, 比如灾变上下限值的划定, 以及预测时间的长短和环境因素的影响等。在进行年份灾变预测时, 数据序列的长短也会影响模型的精度, GM(1, 1)模型一般适用于样本数量少(比如20个以内)时使用, 样本过多难以达到所要求的序列满足的光滑度, 样本过少会导致结果失真, 都会导致预测失败。灾变上下限值的划定, 以及预测时间的长短的选择与优化, 可以通过建立不同灾变预测模型来比较和确定, 从而解决这一难题。同时, 在后续的灾变预报模型建立中, 可考虑影响阿根廷滑柔鱼资源丰度的气候和环境因子, 建立多个GM(1, h)模型, 并与没有考虑影响因子的GM(1, 1)模型进行比较分析。如王言丰等[12]选取灰色关联度较高海区的SST作为预报因子, 构建了GM(0, N)模型较好地拟合了西南大西洋阿根廷滑柔鱼资源丰度的状况。目前大范围的海洋环境数据都来自于卫星遥感, 往往不能满足渔情预报中对数据的及时性和完整性的要求, 会在一定程度上影响到模型的构建[19]。
阿根廷滑柔鱼资源丰度会受到生物与非生物环境因子的影响, 本研究仅对其资源丰度进行了灾变年份的预测, 若要更加准确地预测其资源丰度变化, 在后续的灾变预报模型建立中, 可考虑影响阿根廷滑柔鱼资源丰度的气候和环境因子, 建立多个GM(1, h)模型并进行优化比较, 也可以尝试人工智能等方法, 结合更多的因子来进行预测, 筛选出最优模型对阿根廷滑柔鱼资源丰度进行长期预报。本研究中仅对马尔维纳斯群岛海域阿根廷滑柔鱼资源丰度进行了丰歉年的灾变预测, 结果显示灰色灾变预测在样本量少、信息不完全的情况下, 可以达到较好的结果, 适用于阿根廷滑柔鱼的短期预测, 预测结果具有说服力, 可为马尔维纳斯群岛海域阿根廷柔滑鱼的开发和利用提供相关参考。
[1] |
BRUNETTI N E, ELENA B, ROSSI G R, et al. Summer distribution, abundance and population structure of Illex argentinus on the Argentine shelf in relation to environmental features[J]. South African Journal of Marine Science, 1998, 20(1): 175-186. DOI:10.2989/025776198784126386 |
[2] |
IVANOVIC M L, BRUNETTI N E. Food and feeding of Illex argentinus[J]. Antarctic Science, 1994, 6(2): 185-193. DOI:10.1017/S0954102094000295 |
[3] |
陈新军. 世界头足类资源开发现状及我国远洋鱿钓渔业发展对策[J]. 上海海洋大学学报, 2019, 28(3): 321-330. CHEN X J. Development status of world cephalopod fisheries and suggestions for squid jigging fishery in China[J]. Journal of Shanghai Ocean University, 2019, 28(3): 321-330. |
[4] |
王晓晴. 福克兰(马尔维纳斯)群岛海域鱿钓渔业以及主要许可鱿钓船队发展分析[J]. 渔业信息与战略, 2018, 33(3): 206-214. WANG X Q. Analysis on the development of squid jigging fishery and main licensed jigging fleets in Falkland (Malvinas) Islands waters[J]. Fishery Information & Strategy, 2018, 33(3): 206-214. |
[5] |
刘赫威, 余为, 陈新军. 西南大西洋阿根廷滑柔鱼资源及其对环境响应的研究进展[J]. 中国水产科学, 2020, 27(10): 1254-1265. LIU H W, YU W, CHEN X J. A review of Illex argentinus resources and the responses to environmental variability in the southwest Atlantic Ocean[J]. Journal of Fishery Sciences of China, 2020, 27(10): 1254-1265. |
[6] |
陈新军, 赵小虎. 西南大西洋阿根廷滑柔鱼产量分布与表温关系的初步研究[J]. 大连水产学院学报, 2005, 20(3): 222-228. CHEN X J, ZHAO X H. The relationship between the distribution of production of squid Illex argentinus and sea surface temperature in the Southwest Atlantic Ocean[J]. Journal of Dalian Fisheries University, 2005, 20(3): 222-228. DOI:10.3969/j.issn.1000-9957.2005.03.011 |
[7] |
陈芃, 陈新军. 基于最大熵模型分析西南大西洋阿根廷滑柔鱼栖息地分布[J]. 水产学报, 2016, 40(6): 893-902. CHEN P, CHEN X J. Analysis of habitat distribution of Argentine shortfin squid (Illex argentinus) in the southwest Atlantic Ocean using maximum entropy model[J]. Journal of Fisheries of China, 2016, 40(6): 893-902. |
[8] |
CRESPI-ABRIL A C, BARÓN P J. Revision of the population structuring of Illex argentinus (Castellanos, 1960) and a new interpretation based on modelling the spatio-temporal environmental suitability for spawning and nursery[J]. Fisheries Oceanography, 2012, 21(2/3): 199-214. |
[9] |
CHANG K Y, CHEN C S, WANG H Y, et al. The Antarctic Oscillation index as an environmental parameter for predicting catches of the Argentine shortfin squid (Illex argentinus)(Cephalopoda: Ommastrephidae) in southwest Atlantic waters[J]. Fishery Bulletin, 2015, 113(2): 202-212. DOI:10.7755/FB.113.2.8 |
[10] |
陈新军, 陆化杰, 刘必林, 等. 利用栖息地指数预测西南大西洋阿根廷滑柔鱼渔场[J]. 上海海洋大学学报, 2012, 21(3): 431-438. CHEN X J, LU H J, LIU B L, et al. Forecasting fishing ground of Illex argentinus by using habitat suitability model in the southwest Atlantic[J]. Journal of Shanghai Ocean University, 2012, 21(3): 431-438. |
[11] |
胡贯宇, 陈新军, 汪金涛. 基于不同权重的栖息地指数模型预报阿根廷滑柔鱼中心渔场[J]. 海洋学报, 2015, 37(8): 88-95. HU G Y, CHEN X J, WANG J T. Forecasting fishing ground of Illex argentinus based on different weight habitat suitability index in the southwestern Atlantic[J]. Acta Oceanologica Sinica, 2015, 37(8): 88-95. DOI:10.3969/j.issn.0253-4193.2015.08.008 |
[12] |
王言丰, 陈新军, 陈芃, 等. 基于灰色系统西南大西洋阿根廷滑柔鱼资源丰度预测模型的构建[J]. 海洋学报, 2019, 41(4): 64-73. WANG Y F, CHEN X J, CHEN P, et al. Prediction of abundance index of Argentine shortfin squid in the Southwest Atlantic Ocean based on grey system model[J]. Acta Oceanologica Sinica, 2019, 41(4): 64-73. |
[13] |
王越, 汪金涛, 陈新军, 等. 福克兰海域阿根廷滑柔鱼资源丰度年间变化与其栖息地海表温度的关系[J]. 海洋湖沼通报, 2021, 43(5): 109-114. WANG Y, WANG J T, CHEN X J, et al. Interannual change of Illex argentinus abundance and its relationship with sea surface temperature in Falkland Waters[J]. Transactions of Oceanology and Limnology, 2021, 43(5): 109-114. |
[14] |
解明阳, 陈新军. 基于文献计量学的灰色系统理论在渔业科学中的应用研究进展[J]. 海洋湖沼通报, 2019(5): 117-126. XIE M Y, CHEN X J. Advances in the application of bibliometrics-based grey system theory in fisheries science[J]. Transactions of Oceanology and Limnology, 2019(5): 117-126. |
[15] |
陈新军. 灰色系统理论在渔业科学中的应用[M]. 北京: 中国农业出版社, 2003: 1-32. CHEN X J. Grey system theory in fisheries science[M]. Beijing: China Agriculture Press, 2003: 1-32. |
[16] |
张硕, 李莉, 陈新军. 基于海表温因子的太平洋褶柔鱼冬生群资源丰度预测模型比较[J]. 水产学报, 2018, 42(5): 704-710. ZHANG S, LI L, CHEN X J. A comparative study on forecasting model of the stock abundance index for the winter-spawning cohort of Todarodes pacificus in the Pacific Ocean based on the factor of SST[J]. Journal of Fisheries of China, 2018, 42(5): 704-710. |
[17] |
谢骏, 黄樟翰, 肖学铮, 等. 台湾鳗鲡苗种丰欠的灰色年灾变预测[J]. 水利渔业, 1998(1): 21-22. XIE J, HUANG Z H, XIAO X Z, et al. Grey catastrophe year prediction for the abundance of Taiwan eel seed[J]. Journal of Hydroecology, 1998(1): 21-22. |
[18] |
陈新军, 郑波. 中西太平洋金枪鱼围网渔业鲣鱼资源的时空分布[J]. 海洋学研究, 2007, 25(2): 13-22. CHEN X J, ZHENG B. Study on the temp-spatial distribution of skipjack resources for tuna purse seine in the west-central Pacific Ocean[J]. Journal of Marine Sciences, 2007, 25(2): 13-22. DOI:10.3969/j.issn.1001-909X.2007.02.002 |
[19] |
张磊, 陈新军, 汪金涛, 等. 基于灰色系统的太平洋褶柔鱼秋生群资源丰度灾变预测[J]. 大连海洋大学学报, 2020, 35(4): 607-611. ZHANG L, CHEN X J, WANG J T, et al. Catastrophe prediction of abundance of autumn cohort of pleated squid Todarodes pacificus in Japan Sea and East China Sea based on gray system theory[J]. Journal of Dalian Ocean University, 2020, 35(4): 607-611. |
[20] |
Falkland Islands Government Fisheries Department. Fisheries statistics, volume 21, 2016[EB/OL]. [2021-11-13]. http://www.fis.com/falklandfish/html/publications.html.
|
[21] |
Falkland Islands Government Fisheries Department. Fisheries statistics, volume 22, 2017[EB/OL]. [2021-11-13]. http://www.fis.com/falklandfish/html/publications.html.
|
[22] |
Falkland Islands Government Fisheries Department. Fisheries statistics, volume 23, 2018[EB/OL]. [2021-11-13]. http://www.fis.com/falklandfish/html/publications.html.
|
[23] |
Falkland Islands Government Fisheries Department. Fisheries statistics, volume 24, 2019[EB/OL]. [2021-11-13].
|
[24] |
陈业华, 邱菀华. 灰色灾变预测模型及其应用[J]. 北京航空航天大学学报, 1998, 24(1): 4. CHEN Y H, QIU W H. Grey model of calamity forecast and its application to the calamity forecast[J]. Journal of Beijing University of Aeronautics and Astronautics, 1998, 24(1): 4. |
[25] |
李燕. 灰色预测模型的研究及其应用. 杭州: 浙江理工大学, 2012. LI Y. Research and application of the grey forecast model. Hangzhou: Zhejiang Sci-Tech University, 2012. |
[26] |
刘思峰, 邓聚龙. GM(1, 1)模型的适用范围[J]. 系统工程理论与实践, 2000, 20(5): 121-124. LIU S F, DENG J L. The range suitable for GM(1, 1)[J]. Systems Engineering-Theory & Practice, 2000, 20(5): 121-124. DOI:10.3321/j.issn:1000-6788.2000.05.023 |
[27] |
林东明. 阿根廷滑柔鱼的产卵策略研究. 上海: 上海海洋大学, 2015. LIN D M. Spawning strategy of Argentine shortfin squid, Illex argentinus(Cephalopoda: Ommastrephidae) in the southwest Atlantic. Shanghai: Shanghai Ocean University, 2015. |
[28] |
陈新军, 刘金立. 巴塔哥尼亚大陆架海域阿根廷滑柔鱼渔场分布及与表温的关系分析[J]. 海洋水产研究, 2004, 25(6): 19-24. CHEN X J, LIU J L. Preliminary analysis on the relationship between the distribution of fishing ground of Illex argentius and SST in the Patagonian shelf[J]. Marine Fisheries Research, 2004, 25(6): 19-24. |
[29] |
NIGMATULLIN C M, ZIMIN A V, SUNDAKOV A Z. The stock and fishery variability of the Argentine squid Illex argentinus in 1982-2004 related to environmental conditions. ICES CM 2004/CC: 10, 2004: 1-21.
|
[30] |
郑丽丽, 伍玉梅, 樊伟, 等. 西南大西洋阿根廷滑柔鱼渔场叶绿素a分布及其与渔场的关系[J]. 海洋湖沼通报, 2011(1): 63-70. ZHENG L L, WU Y M, FAN W, et al. The distribution of chlorophyll-a and its relationship with the Illex argentinus fishing ground of Southwest Atlantic Ocean[J]. Transactions of Oceanology and Limnology, 2011(1): 63-70. DOI:10.3969/j.issn.1003-6482.2011.01.009 |
[31] |
刘赫威, 余为, 陈新军, 等. 南极海冰变动对西南大西洋阿根廷滑柔鱼资源丰度及空间分布的影响[J]. 水产学报, 2021, 45(2): 187-199. LIU H W, YU W, CHEN X J, et al. Influence of Antarctic sea ice variation on abundance and spatial distribution of Argentine shortfin squid Illex argentinus in the Southwest Atlantic Ocean[J]. Journal of Fisheries of China, 2021, 45(2): 187-199. |
[32] |
张忠, 刘赫威, 余为, 等. 南极海冰变化驱动的海水降温效应对西南大西洋阿根廷滑柔鱼栖息地的影响[J]. 上海海洋大学学报, 2021, 30(6): 1101-1112. ZHANG Z, LIU H W, YU W, et al. Effects of Antarctic sea ice variations-driven seawater cooling on habitat pattern of Argentine shortfin squid Illex argentinus in the southwest Atlantic Ocean[J]. Journal of Shanghai Ocean University, 2021, 30(6): 1101-1112. |
2. Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai 201306, China;
3. National Engineering Research Center for Oceanic Fisheries, Shanghai 201306, China;
4. Key Laboratory of Oceanic Fisheries Exploration, Ministry of Agriculture and Rural Affairs, Shanghai 201306, China