上海海洋大学学报  2020, Vol. 29 Issue (2): 268-279    PDF    
印度洋长鳍金枪鱼栖息地指数模型的构建与验证
张亚男1, 官文江1,2,3,4, 李阳东1,2,3,4     
1. 上海海洋大学 海洋科学学院, 上海 201306;
2. 农业农村部大洋渔业开发重点实验室, 上海 201306;
3. 农业农村部大洋渔业资源环境科学观测实验站, 上海 201306;
4. 上海市河口海洋测绘工程技术研究中心, 上海 201306
摘要:长鳍金枪鱼(Thunnus alalunga)因其经济价值、分布广泛,成为世界各国海洋渔业的主要捕捞对象之一。对长鳍金枪鱼进行渔场预报研究,可以提高捕捞效率和渔获产量,为渔业生产提供科学依据。选取了2006—2014年印度洋长鳍金枪鱼的生产数据,结合海洋表层的温度、盐度和叶绿素a浓度等3个环境因子,运用一元非线性指数模型按月份建立了各环境因子的印度洋长鳍金枪鱼单因子栖息地适应性指数(suitability index,SI)后,采用算术平均法获得综合栖息地适应性指数模型(habitat suitability index,HSI),并根据2016年印度洋渔业生产数据及相应海洋环境资料,基于地理信息系统(geographic information system,GIS)软件——ArcGIS对栖息地指数模型进行验证。研究表明该模型对渔场预报准确率约为90.56%,各HSI等级下平均预报准确率为87.46%,对于HSI等级4和5(IHSI>0.5)等较高HSI值渔区所代表的中心渔场,其平均准确率为71.82%,考虑到IHSI>0.5的产量平均比重达69.35%这一事实,说明所建立的HSI模型对印度洋长鳍金枪鱼具有较好的预报效果。
关键词栖息地指数模型    渔场预报    长鳍金枪鱼    印度洋    地理信息系统    

长鳍金枪鱼(Thunnus alalunga)是主要的经济性金枪鱼鱼种之一[1],在太平洋、大西洋和印度洋均有分布[2]。由于长鳍金枪鱼是大洋中上层高度洄游性鱼种,其空间分布与环境因子存在着密切联系,因此结合环境因子对该鱼种进行渔场预报是一种较为可靠的思路[3-5]。其中,栖息地指数(habitat suitability index,HSI)基于若干环境因子用来反映生物体对其栖息环境的适应程度,已经广泛地运用到渔场资源预报领域中[6-9]

目前,不少专家学者结合环境因子对太平洋等海域的长鳍金枪鱼渔场进行了研究,栖息地指数模型在东太平洋[10-12]和南太平洋[13-14]已经得到较好的研究和运用,但由于海洋环境的复杂性,栖息地指数模型在其他海域的研究还很匮乏。印度洋海域是世界上重要的金枪鱼生产海区,长鳍金枪鱼是该海域主要的金枪鱼渔获鱼种之一,主要的捕捞国家和地区包括日本、韩国以及我国台湾省,我国大陆在该海区的捕捞产量也在逐步上升[15-16]。关于印度洋长鳍金枪鱼渔场及资源量的研究,陈雪忠等[17]基于随机森林对长鳍金枪鱼渔场进行了预报,官文江等[2]对影响印度洋长鳍金枪鱼资源评估的因素进行了分析,然而相较于其他海域,印度洋长鳍金枪鱼资源量分布的相关研究还需要进一步加强。

为研究该海域长鳍金枪鱼渔场的分布情况及其与环境因子之间的关系,选取2006—2014年和2016年的印度洋金枪鱼渔业委员会(Indian Ocean Tuna Commission,IOTC)长鳍金枪鱼渔获数据,结合Argo海表温度(sea surface temperature,SST)、盐度(sea surface salinity,SSS)和卫星遥感叶绿素a浓度(chlorophyll a concentration,Chl.a)等因子数据,在地理信息系统(geographic information system,GIS)软件——ArcGIS支持下,建立并验证多环境因子的栖息地综合指数模型,探索渔场资源分布与环境因子之间的关系,为印度洋长鳍金枪鱼延绳钓渔业生产提供参考。

1 材料与方法 1.1 研究区域

研究区域为20°E~120°E,20°N~50°S的印度洋海域。图 1显示了根据IOTC延绳钓产量数据统计得到的印度洋海域2006—2014年和2016年共10年时间的5°×5°渔区累积产量分布。在这10年间,作业渔区遍布整个印度洋海域,高产量渔区主要分布在西南印度洋,最高的渔区累积产量约为2.7万t。

图 1 印度洋长鳍金枪鱼渔场10年产量累计分布图 Fig. 1 Cumulative distribution map of Indian Ocean albacore tuna fishing grounds for ten years
1.2 数据来源及处理

主要利用了渔业生产数据、SST、SSS和Chl.a等数据。其中,选取2006—2014年的渔业数据和相关环境因子数据建立HSI模型,而2016年数据则用于对该模型进行精度验证。渔业生产数据来源于IOTC(http://www.iotc.org/)的长鳍金枪鱼生产统计数据,这些数据包括经纬度、作业日期、放钩数、渔获尾数与渔获产量等要素,数据的时间尺度为月,水平空间尺度为1°×1°或5°×5°。其中,研究所需数据是按年分月统计的产量数据,将该数据集两种空间尺度数据进行质量控制(剔除0和1低质量代码的数据)之后,按月份统一网格化成5°×5°渔区产量数据。

SST和SSS数据均来源于中国Argo实时资料中心(http://www.argo.org.cn/)的逐年逐月全球海洋Argo网格数据集产品(BOA-Argo),其水平空间分辨率是1°×1°,垂向共58层。利用ArcGIS 10.3软件的多维工具分别提取温度与盐度数据,并将2006—2014年和2016年数据按月份通过平均值法重采样到与渔区产量数据一致的5°×5°空间分辨率。

Chl.a数据来源于美国国家海洋和大气管理局(National Oceanic and Atmospheric Administration,NOAA)OceanWatch网站(http://oceanwatch.pifsc.noaa.gov/index.html)的Aqua MODIS卫星的月均产品,其空间分辨率为0.05°×0.05°。利用ArcGIS 10.3软件将2006—2014年和2016年数据按月份通过平均值法重采样到与渔区产量数据一致的5°×5°空间分辨率。

1.3 栖息地指数模型的构建

利用2006—2014年的各月渔获产量数据分别与相应时期的SST、SSS和Chl.a数据构建单一因子的适应性指数(Suitability Index,SI)模型。各月单因子SI构建方法如下:

将SST按0.5 ℃、SSS按0.1、Chl.a按0.02 mg/m3组距进行分类,每个组距间隔统计其对应的长鳍金枪鱼总产量。然后按照公式(1)计算各月单因子栖息地指数ISI[18]

    (1)

式中:ISI为各环境因子的适应性指数;Cii分类间隔的捕捞产量;Cmax为各分类间隔中的最大捕捞产量。这里假定分类间隔中最高渔获产量Cmax为长鳍金枪鱼资源最适宜的因子区间,认定其ISI为1;当分类间隔渔获量为0时,则认定长鳍金枪鱼资源最不适宜的因子区间,其ISI为0[19]

使用各月SST、SSS、Chl.a因子各分类间隔的中间值与各间隔相应ISI之间的关系构建一元非线性模型,运用DPS 7.5软件求解一元非线性回归方程参数。一元非线性模型公式为:

    (2)

式中:ISI为各环境因子的适应性指数;X为各环境因子分类间隔所对应中间值;ab为估算参数。

当针对SST、SSS和Chl.a 3种环境因子分别构建单一因子适应性指数模型后,采用算术平均法构建HSI模型[20]IHSI在0~1即不适宜到最适宜之间变化。计算公式为:

    (3)

式中:IHSI为综合栖息地适应性指数,ISI-SSTISI-SSSISI-Chl.a分别代表SST、SSS和Chl.a各环境因子的单因子适应性指数。

1.4 栖息地指数模型的验证

为检验上述HSI模型的预报结果空间分布及精度情况,选用2016年海洋环境数据对长鳍金枪鱼HSI分布进行预报,并利用IOTC实际产量数据对预报结果进行验证。首先,利用2016年5°×5°的月平均SST、SSS和Chl.a 3种环境因子,基于ArcGIS 10.3软件栅格计算功能,使用公式(2)和公式(3)计算覆盖印度洋海域各渔区的IHSI值,即得到2016年各月的HSI预报结果。然后,将2016年各月的HSI预报结果和实际产量数据按自然断裂法(Natural Break-Jenks)[21]从低到高自然划分为5个等级。其中,IHSI值5个等级分别为:等级1(0≤IHSI < 0.1)、等级2(0.1≤IHSI < 0.3)、等级3(0.3≤IHSI < 0.5)、等级4(0.5≤IHSI < 0.7)、等级5(0.7≤IHSI≤1)。最后,将计算及分类所得的2016年各月HSI及其实际产量等级叠加可视化,以分析和对比他们的吻合程度。

为进一步定量和定性分析与检验渔场预报结果,根据上述印度洋长鳍金枪鱼预报HSI分布及其实际产量分级数据,利用ArcGIS 10.3软件提取各渔场分类后的HSI等级,并与实际产量等级做差。对于同一个5°×5°作业渔区,如果预报HSI等级与实际产量等级差小于或等于2,则认为该模型计算的HSI能够准确预测渔场的适宜度; 如果等级差大于2,则认为模型不能正确预测该渔场的适宜度[22]

2 结果 2.1 SI模型的建立

图 2表示了1—12月SST、SSS和Chl.a 3种环境因子的适应性指数散点图及其基于一元非线性指数模型(公式2)建立的适应性指数SI拟合曲线。由图 2中各环境因子散点图和对应的SI拟合曲线可知,印度洋长鳍金枪鱼3种环境因子的适应性均呈季节性变化。

图 2 1—12月SST、SSS和Chl.a因子印度洋长鳍金枪鱼适应性曲线 Fig. 2 Monthly suitability index fitting curves of Indian Ocean albacore tuna to SST, SSS and Chl.a

春季(3—5月)和夏季(6—8月)印度洋长鳍金枪鱼最适SST范围均为17.5~19 ℃。秋季(9—11月)最适范围变化较大:9月为17.5~18.5℃,10月为23.5~25℃,11月为25~27℃。冬季(12—2月)最适范围在25~28.5℃。春夏季(3—8月)的最适SSS范围为35.3~35.7。秋季的9月最适盐度范围为35.5~35.7,10—11月以及冬季的12—2月的最适范围降低为34.8~35.3。

春季最适Chl.a范围在0.06~0.22 mg/m3,夏季为0.12~0.16 mg/m3。秋季最适范围总体在0.04~0.2 mg/m3:9月为0.08~0.2 mg/m3,10月为0.05~0.12 mg/m3,11月为0.04~0.06 mg/m3。冬季最适范围为0.02~0.07 mg/m3

表 1-3为1—12月SST、SSS和Chl.a 3种环境因子对应的SI模型的估算参数、相关系数R2、及各参数相应的P值、95%置信区间。由表 1~表 3可见,除3月和9月外,其余月份SST对应的SI模型R2均大于0.65,拟合系数均通过显著性检验(P < 0.001);除1—3月外,SSS对应的SI模型R2均大于0.70,拟合系数均通过显著性检验(P < 0.001);除8月外,Chl.a对应的SI模型R2均大于0.72,拟合系数均通过显著性检验(P<0.001)。

表 1 SST适应性指数模型统计量 Tab.1 Statistics of SST suitability index model
表 2 SSS适应性指数模型统计量 Tab.2 Statistics of SSS suitability index model
表 3 Chl.a适应性指数模型统计量 Tab.3 Statistics of Chl.a suitability index model
2.2 HSI模型渔场预报及其验证

将计算及分类所得的2016年各月HSI及其实际产量等级叠加可视化,以分析和对比预报结果和实际产量分布的吻合程度,叠加结果如图 3所示。由图 3可知,印度洋海域长鳍金枪鱼实际产量等级>3的作业渔区主要分布在预报HSI等级4和5(IHSI>0.5)较高HSI值渔区,较高HSI分布范围略大于实际高产量分布范围。春季(4—5月)、夏季(6—8月)实际高产量等级与预报的高HSI等级均往印度洋中南部偏移,主要呈现带状分布在50°E~95°E、30°S~40°S附近海域的渔区;秋季(9—11月)渔场开始转移,9月的实际高产量等级与预报高HSI等级仍然主要分布在45°E~70°E、10°S~25°S附近海域的渔区,10—11月的高产量等级与高HSI等级已转移至45°E~65°E、15°S~25°S附近海域的渔区;冬季(12—2月)的实际高产量等级与预报高HSI等级分布在45°E~65°E、15°S~25°S附近海域的渔区。图 3中,仅3月份作为冬、春季过渡月,两种等级分布差异较大:预报较高HSI等级分布在印度洋南部50°E~95°E、35°S~40°S海区,而实际作业高产量等级在45°E~65°E、15°S~20°S与50°E~95°E、35°S~40°S带状海域附近均有分布。

图 3 2016年各月印度洋长鳍金枪鱼栖息地指数HSI及其实际产量等级分布图 Fig. 3 Distribution of monthly HSI and its actual production grades for Indian Ocean albacore tuna in 2016

表 5表 6分别表示渔场总体预报结果和各HSI等级下渔场预报准确率。由表 5可知,2016年1—12月印度洋长鳍金枪鱼渔场平均准确率达90.56%。由表 6可知,低HSI等级(等级1~3)准确率普遍较高;高HSI等级(等级4~5)通常为中心渔场,其准确率明显低于低HSI等级的准确率。表 7统计了2016年各月不同HSI等级下渔获产量占该月总产量的比重。由表 7可知,低HSI等级(等级1~3)下的渔获产量仅占全年总渔获产量的30.65%;高HSI等级(等级4、5)下的渔获产量占全年总渔获产量的69.35%。综合表 5~表 7可知,HSI等级4和5下的平均预报准确率为71.82%,考虑到预报HSI等级4和5的产量占总产量之比为69.35%(29.65%+39.70%),说明HSI模型能够较好地预报印度洋长鳍金枪鱼中心渔场的分布。

表 5 2016年各月渔场预报结果统计 Tab.5 Statistics of monthly fishing ground forecast results in 2016
表 6 2016年各月不同HSI等级下渔场预报准确率统计 Tab.6 Statistics of accuracy for monthly fishing ground forecast according to each HSI grade in 2016
表 7 2016年各月不同HSI等级下渔获产量比重 Tab.7 The proportion of monthly catch production according to each HSI grade in 2016
3 分析与讨论 3.1 各环境因子适应性指数分析

图 2的模型曲线图和表 1~3模型拟合精度可知,印度洋长鳍金枪鱼与SST、SSS和Chl.a 3种环境因子存在显著的指数关系(P均小于0.001),表明基于一元非线性指数模型能够较好地反映印度洋长鳍金枪鱼在各环境因子下的适应性变化趋势。

另外,从图 2的散点图及其拟合的模型曲线可知,印度洋长鳍金枪鱼最适宜温度存在明显的季节性变化,冬季(12—2月)最适温度分布在27.5~28.5 ℃;3月为冬春过渡月最适温度分布在18.5~19 ℃;春、夏季(4—8月)最适温度分布在17.5~18.5 ℃;秋季温度适应性变化较大,9月的最适宜温度接近夏季,10月最适温度分布在23.5~25℃,11月最适温度分布在25~27 ℃,与冬季(1—2月)温度范围趋近。在已有的一些针对南太平洋和东太平洋等其他海域SST与渔场关系的研究中,长鳍金枪鱼最适温度随季节变化不明显。如范永超等和闫敏等在南太平洋长鳍金枪鱼最适SST研究中,范围均在27~29 ℃[23-24]。王震等在东太平洋长鳍金枪鱼研究中作业渔场多分布在SST为24~29 ℃的海域[25]。这些研究中其研究区域一般较小且偏向热带海域,SST季节性变化相对较小,因而研究得出的最适温度随季节变化不明显。而本研究选取的印度洋海域研究区域较大,纬度范围为20°N~50°S,横跨温带和热带海洋,温度随季节性变化明显,从而导致印度洋长鳍金枪鱼在各季节存在较为明显的最适宜温度差异。

盐度对大部分的鱼类分布没有直接影响,主要通过海流和水团等间接影响长鳍金枪鱼的分布[8]。本研究中印度洋长鳍金枪鱼总体适宜盐度大约为34.8~35.7,但也呈现微小的季节性变化,其中秋冬季前后(10—12月和1—2月)偏低一些,范围在34.8~35.3;而春夏季前后(3—9月)偏高一些,范围在35.3~35.7,这与其他海域长鳍金枪鱼盐度影响程度类似[14]。范永超等在南太平洋长鳍金枪鱼HSI模型研究中最适宜盐度为35.2~35.4[23],林显鹏等在所罗门群岛海域长鳍金枪鱼的盐度与环境因子的关系研究中,最适盐度范围在35.6~35.9[26]

浮游植物作为海洋食物链的初级生产者,为整个生态系统提供物质和能量来源,是维持海洋生物生存和发展的基础。其中,Chl.a能够反映浮游植物生物量,基于食物链的原理使以其为食物的浮游动物的生物量发生变动,进而使以浮游动物为食的海洋性鱼种数量和分布发生改变[27]。本研究中,长鳍金枪鱼Chl.a最适宜范围同样受季节的影响,总体范围分布在0.03~0.22 mg/m3,冬季最适宜浓度范围0.06~0.16 mg/m3;4~5月最适宜浓度在0.13~0.22 mg/m3;11-12月最适范围在0.04~0.06 mg/m3。相较于其他海域的长鳍金枪鱼最适宜Chl.a范围大体保持了一致。闫敏等在南太平洋长鳍金枪鱼渔场分布与环境因子的研究中最适Chl.a浓度为0.02~0.08 mg/m3,与本研究相比总体略低[24]

3.2 栖息地指数HSI与渔场分布的关系

使用2016年SST、SSS和Chl.a环境因子数据,基于ArcGIS平台对2016渔场的HSI分布进行了可视化预报及验证,从其他海域长鳍金枪鱼HSI模型仅从统计角度进行验证其精度来看,本文工作更具实际参考意义。印度洋长鳍金枪鱼渔场分布呈季节性变化,模型拟合出的HSI可视化预报结果与实际渔场分布基本一致。1—2月与10—12月主要分布在45°E~70°E、10°S~25°S海域的热带渔区;4—9月主要呈现带状分布在30°S~40°S附近海域的温带渔区。其中,春夏过渡的3月份渔场在温带和热带海域均有分布,结合图 2表 1~3可知,SST、SSS和Chl.a三因子SI模型精度以及模型曲线、拟合散点图均未出现双峰即双适宜度范围的情况,但散点图和其适应性曲线表明:3月份最适宜温度与冬季(1—2月)相差较大,偏向春夏季最适宜温度;最适宜盐度和叶绿素断层虽总体不明显,但同样偏向春夏季最适宜范围。长鳍金枪鱼为高度洄游性大洋鱼种,特别是季节过渡月,一旦因海流、水团、温跃层等其他未知因子的影响,极易造成中心渔场预报的偏移或不准确[14]

从模型统计精度来看,栖息地模型在印度洋区域各HSI等级下月平均准确率为87.46%,结合空间可视化预报结果可知其整体精度略优于其他海域长鳍金枪鱼整体预报率,能够反映该海域整体渔业资源空间分布;然而,高HSI预测准确率上小于低HSI预报准确率,且高HSI等级下渔获量占全年比重约为69.35%,这与其他海域长鳍金枪鱼预报准确率基本类似[8, 13-14]

4 结论

选取2006—2014年印度洋长鳍金枪鱼渔获数据,结合Argo数据集中的SST、SSS和卫星遥感Chl.a等环境因子数据,在基于一元指数模型构建单因素适应性模型的基础上建立了多环境因子的栖息地综合指数模型,并在GIS软件支持下使用2016年实际渔捞产量数据对其进行了验证。研究结果表明:3种环境因子拟合出的单因素适应性模型均通过显著性检验,能够较好反映印度洋长鳍金枪鱼各月的适应性情况,各环境因子的适应性均呈季节性变化;该栖息地指数模型对渔场预报准确率约为90.56%,各HSI等级下平均预报准确率为87.46%,IHSI>0.5所代表的中心渔场平均准确率为71.82%,同时考虑到IHSI>0.5的产量平均比重达69.35%,说明基于SST、SSS和Chl.a 3种环境因子的综合栖息地指数模型具有较好的预报效果。

然而,选取IOTC的5°×5°空间分辨率渔业数据属于较大尺度范围的建模,因此,在印度洋局部小区域范围的准确预报存在局限性,还需更加精细的生产数据进行验证。另外,长鳍金枪鱼渔场分布还与其他海洋环境因子关系密切,例如海流、水团、温跃层、海面高度等[28],都是需要探索的方向。金枪鱼为大洋性中上层高度洄游性鱼种,具有垂直分布的现象[29],其他海域已有结合Argo剖面数据探索金枪鱼垂向分布情况[30-32],但印度洋海域的长鳍金枪鱼垂直分布情况仍需进一步地探索和研究。因此,在以后的研究中需综合考虑其他环境因子、垂向相关因素以及结合其他的函数模型,以完善长鳍金枪鱼栖息地模型。

参考文献
[1]
CHEN I C, LEE P F, TZENG W N. Distribution of albacore (Thunnus alalunga) in the Indian Ocean and its relation to environmental factors[J]. Fisheries Oceanography, 2010, 14(1): 71-80.
[2]
官文江, 朱江峰, 高峰. 印度洋长鳍金枪鱼资源评估的影响因素分析[J]. 中国水产科学, 2018, 25(5): 1102-1114.
GUAN W J, ZHU J F, GAO F. Analysis of influencing factors on stock assessment of the Indian Ocean albacore tuna (Thunnus alalunga)[J]. Journal of Fishery Sciences of China, 2018, 25(5): 1102-1114.
[3]
宋利明, 谢凯, 赵海龙, 等. 库克群岛海域海洋环境因子对长鳍金枪鱼渔获率的影响[J]. 海洋通报, 2017, 36(1): 96-106.
SONG L M, XIE K, ZHAO H L, et al. Effects of environmental variables on catch rates of Thunnus alalunga in waters near Cook Islands[J]. Marine Science Bulletin, 2017, 36(1): 96-106.
[4]
YEN K W, LU H J, CHANG Y, et al. Using remote-sensing data to detect habitat suitability for yellowfin tuna in the Western and Central Pacific Ocean[J]. International Journal of Remote Sensing, 2012, 33(23): 7507-7522.
[5]
毛江美, 陈新军, 余景. 基于神经网络的南太平洋长鳍金枪鱼渔场预报[J]. 海洋学报, 2016, 38(10): 34-43.
MAO J M, CHEN X J, YU J. Forecasting fishing ground of Thunnus alalunga based on BP neural network in the South Pacific Ocean[J]. Haiyang Xuebao, 2016, 38(10): 34-43.
[6]
陈程, 陈新军, 雷林, 等. 基于栖息地指数模型的毛里塔尼亚头足类底拖网渔场研究[J]. 上海海洋大学学报, 2017, 26(6): 939-945.
CHEN C, CHEN X J, LEI L, et al. Study on fishing ground of bottom trawl based on the habitat suitability model in Mauritania[J]. Journal of Shanghai Ocean University, 2017, 26(6): 939-945.
[7]
王易帆, 陈新军. 中西太平洋鲣栖息地指数预报模型比较研究[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.
[8]
杨嘉樑, 黄洪亮, 宋利明, 等. 基于分位数回归的库克群岛海域长鳍金枪鱼栖息环境综合指数[J]. 中国水产科学, 2014, 21(4): 832-851.
YANG J L, HUANG H L, SONG L M, et al. An integrated habitat index for albacore tuna (Thunnus alalunga) inthe waters near the Cook Islands based on the quantile regression method[J]. Journal of Fishery Sciences of China, 2014, 21(4): 832-851.
[9]
佟佳琦, 陈锦辉, 高春霞, 等. 基于栖息地适应性指数的长江口刀鲚时空分布特征[J]. 上海海洋大学学报, 2018, 27(4): 584-593.
TONG J Q, CHEN J H, GAO C X, et al. Temporal-spatial distribution of Coilia nasus in the Yangtze River Estuary based on habitat suitability index[J]. Journal of Shanghai Ocean University, 2018, 27(4): 584-593.
[10]
任中华, 陈新军, 方学燕. 基于栖息地指数的东太平洋长鳍金枪鱼渔场分析[J]. 海洋渔业, 2014, 36(5): 385-395.
REN Z H, CHEN X J, FANG X Y. Forecasting fishing grounds of Thunnus alalunga in the eastern Pacific based on habitat suitability index[J]. Marine Fisheries, 2014, 36(5): 385-395.
[11]
沈智宾, 陈新军, 汪金涛. 基于海表温度和海面高度的东太平洋大眼金枪鱼渔场预测[J]. 海洋科学, 2015, 39(10): 45-51.
SHEN Z B, CHEN X J, WANG J T. Forecasting of Bigeye tuna fishing ground in the Eastern Pacific Ocean based on sea surface temperature and sea surface height[J]. Marine Sciences, 2015, 39(10): 45-51.
[12]
张艳波.东南太平洋长鳍金枪鱼垂直分布与栖息地指数研究[D].上海: 上海海洋大学, 2015.
ZHANG Y B. Vertical distribution and integrated habitat index for albacore (Thunnus alalunga) in the southeast Pacific Ocean[D]. Shanghai: Shanghai Ocean University, 2015.
[13]
范江涛, 陈新军, 钱卫国, 等. 南太平洋长鳍金枪鱼渔场预报模型研究[J]. 广东海洋大学学报, 2011, 31(6): 61-67.
FAN J T, CHEN X J, QIAN W G, et al. Models of forecasting fishing ground of Thunnus alalunga in the South Pacific Ocean[J]. Journal of Guangdong Ocean University, 2011, 31(6): 61-67.
[14]
郭刚刚, 张胜茂, 樊伟, 等. 基于表层及温跃层环境变量的南太平洋长鳍金枪鱼栖息地适应性指数模型比较[J]. 海洋学报, 2016, 38(10): 44-51.
GUO G G, ZHANG S M, FAN W, et al. Comparative research on habitat suitability index models of albacore tuna (Thunnus alalunga) based on surface and thermocline environmental variables in the South Pacific Ocean[J]. Acta Oceanologica Sinica, 2016, 38(10): 44-51.
[15]
朱江峰, 戴小杰, 官文江. 印度洋长鳍金枪鱼资源评估[J]. 渔业科学进展, 2014, 35(1): 1-8.
ZHU J F, DAI X J, GUAN W J. Stock assessment of albacore Thunnus alalunga in the Indian Ocean[J]. Progress in Fishery Sciences, 2014, 35(1): 1-8.
[16]
马璐璐, 朱江峰, 耿喆, 等. 运用生物量动态模型评估印度洋长鳍金枪鱼资源[J]. 上海海洋大学学报, 2018, 27(2): 259-264.
MA L L, ZHU J F, GENG Z, et al. Stock assessment of albacore (Thunnus alalunga) in the Indian Ocean usingbiomass dynamics model[J]. Journal of Shanghai Ocean University, 2018, 27(2): 259-264.
[17]
陈雪忠, 樊伟, 崔雪森, 等. 基于随机森林的印度洋长鳍金枪鱼渔场预报[J]. 海洋学报, 2013, 35(1): 158-164.
CHEN X Z, FAN W, CUI X S, et al. Fishing ground forecasting of Thunnus alalunga in Indian Ocean based on random forest[J]. Acta Oceanologica Sinica, 2013, 35(1): 158-164.
[18]
官文江, 高峰, 雷林, 等. 多种数据源下栖息地模型及预测结果的比较[J]. 中国水产科学, 2015, 22(1): 149-157.
GUAN W J, GAO F, LEI L, et al. Comparisons of the habitat suitability index models developed by multi-source data and forecasting[J]. Journal of Fishery Sciences of China, 2015, 22(1): 149-157.
[19]
SONG L M, ZHOU Y Q. Developing an integrated habitat index for bigeye tuna (Thunnu sobesus) in the Indian Ocean based on longline fisheries data[J]. Fisheries Research, 2010, 105(2): 63-74.
[20]
CHEN X J, LI G, FENG B, et al. Habitat suitability index of Chub mackerel (Scomber japonicus) from July to September in the East China Sea[J]. Journal of Oceanography, 2009, 65(1): 93-102.
[21]
JENKS G F. Optimal data classification for choroplethmaps[M]. Kansas: University of Kansas, Occasional Paper 2, 1977.
[22]
高峰, 陈新军, 范江涛, 等. 西南大西洋阿根廷滑柔鱼中心渔场预报的实现及验证[J]. 上海海洋大学学报, 2011, 20(5): 754-758.
GAO F, CHEN X J, FAN J T, et al. Implementation and verification of intelligent fishing ground forecasting of Illex argentinus in the Southwest Atlantic[J]. Journal of Shanghai Ocean University, 2011, 20(5): 754-758.
[23]
范永超, 陈新军, 汪金涛. 基于多因子栖息地指数模型的南太平洋长鳍金枪鱼渔场预报[J]. 海洋湖沼通报, 2015(2): 36-44.
FAN Y C, CHEN X J, WANG J T. Forecasting central fishing ground of Thunnus alalunga based on multi-factors habitat suitability index in the South Pacific Ocean[J]. Transactions of Oceanology and Limnology, 2015(2): 36-44.
[24]
闫敏, 张衡, 樊伟, 等. 南太平洋长鳍金枪鱼渔场CPUE时空分布及其与关键海洋环境因子的关系[J]. 生态学杂志, 2015, 34(11): 3191-3197.
YAN M, ZHANG H, FAN W, et al. Spatial-temporal CPUE profiles of the albacore tuna (Thunnus alalunga) and their relations to marine environmental factors in the South Pacific Ocean[J]. Chinese Journal of Ecology, 2015, 34(11): 3191-3197.
[25]
王震, 陈新军, 雷林. 东太平洋长鳍金枪鱼栖息地指数模型的比较[J]. 广东海洋大学学报, 2017, 37(1): 58-64.
WANG Z, CHEN X J, LEI L. Comparison of Thunnus alalunga in the Eastern Pacific based on habit atsuitability index model[J]. Journal of Guangdong Ocean University, 2017, 37(1): 58-64.
[26]
林显鹏, 郭爱, 张洪亮, 等. 所罗门群岛海域长鳍金枪鱼的垂直分布与环境因子的关系[J]. 浙江海洋学院学报(自然科学版), 2011, 30(4): 303-306.
LIN X P, GUO A, ZHANG H L, et al. Relationship between vertical distribution of albacore tuna Thunnus alalunga and the concerned environmental factors in the Solomon Sea[J]. Journal of Zhejiang Ocean University (Natural Science), 2011, 30(4): 303-306.
[27]
于杰, 陈国宝, 张魁, 等. 南海中部海域夏季叶绿素a浓度垂向分布特征[J]. 南方水产科学, 2016, 12(4): 1-8.
YU J, CHEN G B, ZHANG K, et al. Vertical distribution of summer chlorophyll a concentrationin the middle South China Sea[J]. South China Fisheries Science, 2016, 12(4): 1-8.
[28]
宋婷婷, 樊伟, 伍玉梅. 卫星遥感海面高度数据在渔场分析中的应用综述[J]. 海洋通报, 2013, 32(4): 474-480.
SONG T T, FAN W, WU Y M. Review on the fishery analysis with satellite remote sensing sea surface height data[J]. Marine Science Bulletin, 2013, 32(4): 474-480.
[29]
HOLLAND K N. Horizontal and vertical movements of Yellowfin and Bigeye tuna associated with fish aggregating devices[J]. Journal of Cell Biology, 1990, 148(5): 957-970.
[30]
杨胜龙, 张禹, 张衡, 等. 热带印度洋黄鳍金枪鱼渔场时空分布与温跃层的关系[J]. 生态学报, 2012, 32(3): 671-679.
YANG S L, ZHANG Y, ZHANG H, et al. The relationship between the temporal-spatial distribution of fishing ground of yellowfin tuna (Thunnus albacares) and themocline characteristics in the tropic Indian Ocean[J]. Acta Ecologica Sinica, 2012, 32(3): 671-679.
[31]
ASHLEY J W, ALLAIN V, NICOL S J, et al. Vertical behavior and diet of albacore tuna (Thunnu salalunga) vary with latitude in the South Pacific Ocean[J]. Deep Sea Research Part Ⅱ:Topical Studies in Oceanography, 2015, 113: 154-169.
[32]
郭刚刚, 张胜茂, 樊伟, 等. 南太平洋长鳍金枪鱼垂直活动水层空间分析[J]. 南方水产科学, 2016, 12(5): 123-130.
GUO G G, ZHANG S M, FAN W, et al. Spatial analysis of vertical active layer of albacore tuna (Thunnus alalunga) in the South Pacific Ocean[J]. South China Fisheries Science, 2016, 12(5): 123-130.
Construction and verification of a habitat suitability index model for the Indian Ocean albacore tuna
ZHANG Yanan1, GUAN Wenjiang1,2,3,4, LI Yangdong1,2,3,4     
1. College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China;
2. Key Laboratory of Oceanic Fisheries Exploration, Ministry of Agriculture and Rural Affairs, Shanghai 201306, China;
3. Scientific Observing and Experimental Station of Oceanic Fishery Resources, Ministry of Agriculture and Rural Affairs, Shanghai 201306, China;
4. Shanghai Engineering Research Center of Estuarine and Oceanographic Mapping, Shanghai 201306, China
Abstract: Because of its economic value and wide distribution, albacore tuna (Thunnus alalonga) has become one of the main fishing targets in the world's marine fisheries. Prediction of albacore tuna can improve fishing efficiency and yield, and provide scientific basis for fishery production.The production data of albacore tuna in the Indian Ocean from 2006 to 2014 and three environmental factors of temperature, salinity and chlorophyll a concentration in the surface layer of the ocean were used in this study. The single-factor Suitability Index (SI) of Indian Ocean albacore tuna with various environmental factors was established monthly by using the single-variable non-linear index model. Then the arithmetic average method was used to obtain the comprehensive Habitat Suitability Index (HSI) model. Using the 2016 Indian ocean albacore tuna production data and the corresponding marine environment data, furthermore, the HSI model was verified based on the ArcGIS platform. The results show that the accuracy of the monthly fishing ground forecast is about 90.56% and the overall forecast accuracy for each HSI grade is 87.46%. Moreover, for the central fishing ground with IHSI>0.5, the average accuracy rate is 71.82%. Considering that the average yield ratio of IHSI>0.5 is 69.35%, it can be concluded that the established HSI model had a promising forecast effect for the Indian Ocean albacore tuna.
Key words: habitat suitability index model     fishing ground forecast     albacore tuna(Thunnus alalunga)     Indian Ocean     geographic information system