2. 农业部大洋渔业开发重点实验室, 上海 201306;
3. 国家远洋渔业工程技术研究中心, 上海 201306;
4. 农业部大洋渔业资源环境科学观测实验站, 上海 201306
鲣鱼(Katsuwonus pelamis)广泛分布于太平洋、大西洋、印度洋的热带、亚热带以及亚寒带海域, 在世界金枪鱼渔业中占有极其重要的地位[1]。世界主要金枪鱼产量的60%~70%来自太平洋海域[2],其中又有76%来自中西太平洋[3]和56%的产量来自围网捕捞[3-4],围网捕捞中的鲣鱼渔获比例约为70%~80%[5]。我国于2000年开始在该海域进行金枪鱼围网作业, 至2006年已发展到围网船只8艘、年产量约5万吨的规模[6],而近些年已达到围网船只20艘、年产量10万吨左右。国内外已有很多学者对中西太平洋鲣鱼资源渔场分布及其与海洋环境的关系做了大量研究[7-13]。研究表明,中西太平洋暖池与鲣鱼资源渔场分布关系密切。HAMPTON等[7-8]认为,鲣鱼分布会随厄尔尼诺—南方涛动现象(ENSO)产生相应的迁移。周甦芳[9]认为,厄尔尼诺发生时, 鲣鱼围网单位捕捞努力量渔获量经度重心较正常年份向东偏10~20个经度,拉尼娜年则向西偏10~20个经度。汪金涛等[10]认为,当Nino 3.4区海表温度异常值从低到高变化时,鲣鱼渔场重心也逐渐由西向东偏。以上研究均表明,渔场分布受气候环境变化的影响很大,与海表温度存在密切关系。在中西太平洋海域共有12个国家和地区,我国捕捞作业渔场基本上处在他国专属经济区管辖范围内[11-12]。现阶段我国金枪鱼围网渔业可持续发展面临着中长期金枪鱼资源与空间分布无法预测,且各岛国实行作业天数限制,导致入渔成本不断提高,因而影响金枪鱼围网的作业效率及产生的经济效益有所下降。中西太平洋是我国鲣鱼围网的重要作业海域,准确地预报该海域的中心渔场有利于提高捕捞效率。为此,本研究根据1995—2014年中西太平洋鲣鱼围网生产数据以及海洋环境数据,统计分析不同气候条件下渔业生产数据与气候因子之间的关系,以渔获量作为中心渔场的预报指标,将研究区域划分为22个5°×5°的海区,分别统计每个海区的渔获量与气候因子之间的关系并建立每个渔区的渔场预报模型,为准确快速寻找渔场提供科学依据。
1 材料与方法 1.1 材料来源中西太平洋鲣鱼围网渔获生产统计数据来源于南太平洋渔业委员会(http://www.spc.gov)。数据包括日本、韩国、中国、澳大利亚、美国、西班牙和南太平洋岛国等国家和地区,统计内容包括年、月、经度、纬度、投网次数以及渔获量。研究区域为5°S~5°N,125°E~180°E。时间为1995—2014年。空间分辨率为5°×5°,时间分辨率为月。
ENSO指标拟用Nino 3.4区海表温度距平值(SSTA)来表示,其数据来自美国NOAA气候预报中心(http://www.cpc.neep.noaa.gov/),时间单位为月。
1.2 研究方法采用美国NOAA气候预报中心的标准定义推断厄尔尼诺/拉尼娜事件[13-15]:Nino 3.4区的温度距平值连续3个月超过+0.5 ℃时, 则认为是发生了厄尔尼诺现象;Nino 3.4区的温度距平值连续3个月低于-0.5 ℃时, 则认为发生了拉尼娜现象;其他的情况则为正常情况。根据以上定义推断1995—2014年发生的异常环境事件(表 1),并根据表 1,将1995—2014年,以月为分辨率,从时间上将20年分为3个部分,分别为厄尔尼诺月份、拉尼娜月份、正常月份。
厄尔尼诺和拉尼娜现象会极大地影响鲣鱼的资源丰度[13],进而影响到渔场的分布及其渔获量。为此,统计不同区域渔获量的分布来研究鲣鱼的空间分布特征。
将1995—2014年渔业捕捞数据与环境数据(SSTA)相匹配,并根据以上对异常环境事件的定义以及不同月份的划分结果,将1995—2014年的统计数据分为3类,分别为厄尔尼诺月份数据、拉尼娜月份数据、正常月份数据。
将研究区域5°S~5°N,125°E~180°E以5°×5°为空间统计单位划分为22个海区进行统计。对每个海区分别统计不同SSTA范围(以0.5 ℃为间隔)与所对应的初值化渔获量的关系,并利用正态分布模型建立每个海区的渔场预报模型,所对应的y值为该海域占总海域产量的百分比。在3种气候条件下,分别选取80%的数据进行上述统计及建模,并利用剩余的20%数据进行验证,模型的优劣通过误差(变异系数,coefficient of variation,Cv)进行判断。变异系数Cv的计算方式如下:
式中:Cv为变异系数,Vstd为预测值标准差与实际值标准差之差,Vmean为预测值的平均值与实际值平均值之差。
2 结果 2.1 不同气候条件下产量情况统计不同气候条件下(厄尔尼诺、拉尼娜和正常月份)中西太平洋鲣鱼月总产量,结果发现:1995—2014年间厄尔尼诺月总产量为270.9万吨,拉尼娜月总产量为269.9万吨,正常月总产量为253.9万吨;从月平均产量来看,厄尔尼诺月均产量为2 997.4 t,拉尼娜月均产量为2 986.4 t,正常月均产量为2 809.4 t。总体来看,相比其他气候条件而言,正常年份产量要明显偏低(图 1)。
根据不同空间位置统计20年的鲣鱼产量,结果发现,不同空间以及不同气候条件下产量的有较大的差异。在130°E~145°E海域,北纬5°产量明显高于南纬5°产量,同时在这一海域,正常月份的产量均高于20万t,比其他极端气候的产量高;而在145°E~165°E以东海域,无论南北纬,产量均相对较高,且拉尼娜月份产量均高于20万t,比厄尔尼诺/正常月份的产量高;在165°E以东海域,产量相对较低,其中南纬5°海域中,拉尼娜月份产量明显比厄尔尼诺/正常月份的产量低(图 2)。
分析表明,在厄尔尼诺气候条件下,研究海区的SSTA和22个海区的渔获量之间均呈正态分布,相关系数均达到0.95以上(P < 0.01)(表 2)。所建立的模型均能很好地拟合SSTA与渔获量之间的关系。
分析表明,在拉尼娜气候条件下,研究海区的SSTA和22个海区的渔获量之间均呈正态分布,相关系数均在0.75以上,大多数均超过0.9(P < 0.01,表 3)。所建立的模型均能很好地拟合SSTA与渔获量之间的关系。
分析表明,在正常气候条件下,研究海区的SSTA和22个海区的渔获量之间均呈正态分布,相关系数均在0.95以上(P < 0.01,表 4)。所建立的模型均能很好地拟合SSTA与渔获量之间的关系。
根据80%的数据进行建模所得模型,利用剩余20%的数据进行模型验证,将不同海区的SSTA数据各自代入模型中,将获得的预报值与实际值进行比较。结果表明,3种不同气候条件下的预报结果与实际统计值之间均存在显著相关关系(P < 0.01,图 5)。同时,这3种不同气候条件下所建立的模型均有着较高的相关系数,超过0.6,拟合结果较好(图 5)。从误差(变异系数`Cv)来看,3种不同的气候条件下建立的模型Cv范围在4.9%~6.7%,均小于10%的误差区间。
陈新军等[6]认为,中西太平洋鲣鱼围网高产渔区空间位置主要集中在5°S~5°N、130°E~175°E海域,因此本研究中,以5°S~5°N,125°E~180°E为研究海域进行分析。厄尔尼诺等极端气候条件会极大地影响鲣鱼的渔获量[13],厄尔尼诺/拉尼娜月份的产量相对较高,这与中西太平洋暖池中心位置的变化有着密切的关系[16]。
在130°E~145°E海域,北纬产量明显高于南纬(图 2),这可能是由于北纬地区所占的海域面积较大,可捕捞的面积也较大,而南纬地区大部分为陆地,也直接导致北纬有较高的渔获量。145°E~165°E以东海域的产量较高,拉尼娜月份产量明显高于其他气候条件下的产量(图 2)。鲣鱼的主要栖息地会随着海表温的变化而变化[17]。根据渔获量重心计算可知,在拉尼娜条件下,暖池位置较正常年份向西偏移,鲣鱼栖息位置偏西北[10, 17]。赤道的上升流海域丰富的饵料生物会随着季风向西输送[18],因此在拉尼娜条件下比较适合鲣鱼的生长。郭爱等[18]的研究认为,气候正常年份CPUE值相对偏低,厄尔尼诺和拉尼娜年份CPUE值高于正常年份。在165°E以东海域,产量相对较低,且在拉尼娜年份下的渔获量也迅速下降,这可能是在拉尼娜条件下,海表温相对较低,更不利于鲣鱼的生长,由于不是鲣鱼理想的栖息场所,因此也直接导致了渔获量的下降。
在SSTA与渔获量的关系中可以发现,即使在不同的气候条件下,两者关系均符合正态分布(表 2~4)。相比较而言,拉尼娜年份的相关系数相对偏低,但是在其渔获量较高的145°E~160°E海域,对应的相关系数仍然较高(表 3)。这种变化与大尺度的环境变化有着密切的关系。沈建华等[17]和周甦芳等[9, 13]对于中西太平洋鲣鱼渔场及时空分布与ENSO之间关系进行研究发现,厄尔尼诺发生时,鲣鱼围网CPUE经度重心随着暖池的东扩而东移,拉尼娜发生时则随着暖池向西收缩而西移。鲣鱼的栖息地与中西太平洋暖池的范围有关:在厄尔尼诺发生时,暖池的东扩使得鲣鱼的适宜栖息范围变大,因此整个研究海域的温度变化可以较好地与鲣鱼渔获量建立关系;而拉尼娜发生时,暖池向西收缩,鲣鱼作为一种暖水性鱼类,对水温的变化较为敏感,这也会直接影响到其捕食活动,因此该海域的渔获量明显减少。
模型的验证结果表明,SSTA能较好地表现出渔获量的变化规律,预测值与实际值百分比有着较大的关联度(图 5)。SSTA反映出了海表温在一定时间内的变化,这可以很直观地分析出不同环境条件下鲣鱼资源量的变化。在上述分析中,暖池的移动直接导致了鲣鱼的渔获量在空间分布上发生了较大的变化,同时暖池的移动与相应的气候变化有着直接的关系,因此也可以认为鲣鱼渔获量在不同气候条件下发生了变化。因此,在对鲣鱼渔场预报的过程中,需要区分不同气候条件,以达到更为准确的预报结果。以往的相似研究主要对CPUE或捕捞努力量的变化进行分析[19],主要反映了相对资源量的丰度,在历年CPUE或捕捞努力量波动较大的情况下,更能反映渔场的真实情况。但从渔场预报的角度来说,直接获取的渔获量能够更加直观地表达该地区的资源量状况,尤其近年来在捕捞努力量已经趋于稳定的前提下,研究渔获量的变化对渔场预测更具有指导意义。
本文根据中西太平洋鲣鱼围网捕捞历史渔业生产数据,分析了渔获量在不同气候条件下的变化,研究了不同空间位置中3种气候条件鲣鱼渔获量波动规律,为预测渔场和相关预报工作提供了新的思路。之前的研究[20]主要针对不同气候条件下的鲣鱼资源丰度进行,为了更完整地体现整个变化规律,减少捕捞的盲目性从而提高捕捞效率,在后续研究中将进一步查明厄尔尼诺、拉尼娜对暖池变化的影响,引入捕捞努力量,并采用多种预报模型进行对比,包括BP神经网络模型[21-22]、栖息地指数模型[23-25]、灰色系统等[26]。进而更全面地对鲣鱼渔场进行预报,为提升中西太平洋鲣鱼的渔情预报精确度提供参考依据。
[1] |
COLLETTE B B, NAUEN C E. FAO species catalogue:vol. 2. Scombrids of the world. An annotated and illustrated catalogue of tunas, mackerels, bonitos and related species known to date[J]. FAO Fisheries Synopsis, 1983, 125(2): 83-86. |
[2] |
苗振清, 黄锡昌. 远洋金枪鱼渔业[M]. 上海: 上海科学技术文献出版社, 2003: 1-242. MIAO Z Q, HUANG X C. Distant water tuna fishery[M]. Shanghai: Shanghai Science and Technology Literature Press, 2003: 1-242. |
[3] |
谢营梁. 中西太平洋金枪鱼渔业的现状和展望(1)[J]. 现代渔业信息, 2002, 17(9): 18-20. XIE Y L. Status and perspectives of the western and central pacific ocean tuna fisheries (1)[J]. Modern Fisheries Information, 2002, 17(9): 18-20. DOI:10.3969/j.issn.1004-8340.2002.09.005 |
[4] |
齐建军, 王明彦. 世界金枪鱼围网渔业的现状与前景[J]. 海洋渔业, 2001, 23(1): 18-20. QI J J, WANG M Y. The present status and prospects of tuna purse seine fishery in the world[J]. Marine Fisheries, 2001, 23(1): 18-20. DOI:10.3969/j.issn.1004-2490.2001.01.005 |
[5] |
邹志来. 发展金枪鱼围网渔业的可行性研究报告[J]. 远洋渔业, 2000(6): 1-3. ZOU Z L. Feasibility study report of developing purse seine tuna fishery[J]. Distant Water Fishery, 2000(6): 1-3. |
[6] |
陈新军, 郑波. 中西太平洋金枪鱼围网渔业鲣鱼资源的时空分布[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 |
[7] |
HAMPTON J, LEWIS A, WILLIAMS P. The western and central Pacific tuna fishery: over view and status of stocks[R]. Oceanic Fisheries Programme SPC, 1999: 39.
|
[8] |
HAMPTON J. Estimation of tag reporting and tag-shedding rates in a large-scale tuna tagging experiment in the western tropical Pacific Ocean[J]. Fish Bulletin, 1997, 95: 68-79. |
[9] |
周甦芳. 厄尔尼诺-南方涛动现象对中西太平洋鲣鱼围网渔场的影响[J]. 中国水产科学, 2005, 12(6): 739-744. ZHOU S F. Impacts of the El Nino southern Oscillation on skipjack tuna purse-seine fishing grounds in the western and central Pacific Ocean[J]. Journal of Fishery Sciences of China, 2005, 12(6): 739-744. DOI:10.3321/j.issn:1005-8737.2005.06.012 |
[10] |
汪金涛, 陈新军. 中西太平洋鲣鱼渔场的重心变化及其预测模型建立[J]. 中国海洋大学学报, 2013, 43(8): 44-48. WANG J T, CHEN X J. Changes and prediction of the fishing ground gravity of skipjack (Katsuwonus pelamis) in western-central Pacific[J]. Periodical of Ocean University of China, 2013, 43(8): 44-48. |
[11] |
《世界大洋性渔业概况》编写组. 世界大洋性渔业概况[M]. 北京: 海洋出版社, 2011. World oceanic Fishery Status. World oceanic fishery status[M]. Beijing: China Ocean Press, 2011. |
[12] |
农业部渔业渔政管理局. 中国渔业年鉴2014[M]. 北京: 海洋出版社, 2014. The Ministry of Agriculture, Fisheries Fishery Administration. China fisheries yearbook[M]. Beijing: China Ocean Press, 2014. |
[13] |
周甦芳, 沈建华, 樊伟. ENSO现象对中西太平洋鲣鱼围网渔场的影响分析[J]. 海洋渔业, 2004, 26(3): 168-172. ZHOU S F, SHEN J H, FAN W. Impacts of the El Nino southern oscillation on skipjack tuna (Katsuwonus pelamis) purse-seine fishing grounds in the western and central Pacific Ocean[J]. Marine Fisheries, 2004, 26(3): 168-172. |
[14] |
张启龙, 翁学传, 侯一筠, 等. 西太平洋暖池表层暖水的纬向运移[J]. 海洋学报, 2004, 26(1): 33-39. ZHANG Q L, WENG X C, HOU Y J, et al. Zonal movement of surface warm water in the western pacific warm pool[J]. Acta Oceanologica Sinica, 2004, 26(1): 33-39. DOI:10.3321/j.issn:0253-4193.2004.01.004 |
[15] |
钟姗姗, 何金海, 刘宣飞, 等. ENSO十年际尺度变率的机制初探[J]. 热带海洋学报, 2004, 23(2): 28-36. ZHONG S S, HE J H, LIU X F, et al. Possible mechanisms of ENSO'S decadal variability[J]. Journal of Tropical Oceanography, 2004, 23(2): 28-36. DOI:10.3969/j.issn.1009-5470.2004.02.004 |
[16] |
LEHODEY P, BERTIGNAC M, HAMPTON J, et al. El Nino southern oscillation and tuna in the western Pacific[J]. Nature, 1997, 389(6652): 715-718. DOI:10.1038/39575 |
[17] |
沈建华, 陈雪冬, 崔雪森. 中西太平洋金枪鱼围网鲣鱼渔获量时空分布分析[J]. 海洋渔业, 2006, 28(1): 13-19. SHEN J H, CHEN X D, CUI X S. Analysis on spatial-temporal distribution of skipjack tuna catches by purse seine in the western and central Pacific Ocean[J]. Marine Fisheries, 2006, 28(1): 13-19. DOI:10.3969/j.issn.1004-2490.2006.01.003 |
[18] |
郭爱, 陈新军. ENSO与中西太平洋金枪鱼围网资源丰度及其渔场变动的关系[J]. 海洋渔业, 2005, 27(4): 338-342. GUO A, CHEN X J. The relationship between ENSO and tuna purse-seine resource abundance and fishing grounds distribution in the western and central Pacific Ocean[J]. Marine Fisheries, 2005, 27(4): 338-342. DOI:10.3969/j.issn.1004-2490.2005.04.015 |
[19] |
TIAN S Q, CHEN X J, CHEN Y, et al. Evaluating habitat suitability indices derived from CPUE and fishing effort data for Ommatrephes bratramii in the northwestern Pacific Ocean[J]. Fisheries Research, 2009, 95(2/3): 181-188. |
[20] |
陈洋洋, 陈新军. 厄尔尼诺/拉尼娜现象对中西太平洋鲣资源丰度的影响[J]. 上海海洋大学学报, 2017, 26(1): 113-120. CHEN Y Y, CHEN X J. Influence of El Nino/La Nina on the abundance index of skipjack in the western and central Pacific Ocean[J]. Journal of Shanghai Ocean University, 2017, 26(1): 113-120. |
[21] |
魏联, 陈新军, 雷林, 等. 西北太平洋柔鱼BP神经网络渔场预报模型比较研究[J]. 上海海洋大学学报, 2017, 26(3): 450-457. WEI L, CHEN X J, LEI L, et al. Comparative study on the forecasting models of squid fishing ground in the northwest Pacific Ocean based on BP artificial neural network[J]. Journal of Shanghai Ocean University, 2017, 26(3): 450-457. |
[22] |
苏新红, 金丰军, 杨奇志, 等. 基于BP神经网络模型的福建海域赤潮预报方法研究[J]. 水产学报, 2017, 41(11): 1744-1755. SU X H, JIN F J, YANG Q Z, et al. Red tide forecasting model based on BP neural network in Fujian sea area[J]. Journal of Fisheries of China, 2017, 41(11): 1744-1755. |
[23] |
蒋瑞, 陈新军, 雷林, 等. 秋冬季智利竹(竹夹)鱼栖息地指数模型比较[J]. 水产学报, 2017, 41(2): 240-249. JIANG R, CHEN X J, LEI L, et al. A comparative study on habitat suitability index of Trachurus murphyi during autumn and winter[J]. Journal of Fisheries of China, 2017, 41(2): 240-249. |
[24] |
王易帆, 陈新军. 中西太平洋鲣栖息地指数预报模型比较研究[J]. 上海海洋大学学报, 2017, 26(5): 743-750. WANG Y F, CHEN X J. Comparisons of habitat suitability index models of skipjack tuna in the Western and Central Pacific Ocean[J]. Journal of Shanghai Ocean University, 2017, 26(5): 743-750. |
[25] |
余为, 陈新军, 易倩. 不同气候模态下西北太平洋柔鱼渔场环境特征分析[J]. 水产学报, 2017, 41(4): 525-534. YU WEI, CHEN X J, YI Q. Analysis of variations in the environmental conditions on the fishing ground of neon flying squid (Ommastrephesbartramii) in the Northwestern Pacific Ocean under different climate modes[J]. Journal of Fisheries of China, 2017, 41(4): 525-534. |
[26] |
段丁毓, 陈芃, 陈新军, 等. 基于灰色系统的秘鲁鳀资源量预测模型的构建[J]. 上海海洋大学学报, 2018, 27(2): 284-290. DUAN D Y, CHEN P, CHEN X J, et al. The construction of biomass forecasting model for the anchoveta (Engraulisringens) by the grey system model[J]. Journal of Shanghai Ocean University, 2018, 27(2): 284-290. |
2. Key Laboratory of Oceanic Fisheries Exploration, Ministry of Agriculture, Shanghai 201306, China;
3. National Engineering Research Center for Oceanic Fisheries, Shanghai 201306, China;
4. Scientific Observing and Experimental Station of Oceanic Fishery Resources, Ministry of Agriculture, Shanghai 201306, China