中西太平洋延绳钓黄鳍金枪鱼渔场时空分布与温跃层关系

杨胜龙，张忭忭，靳少非，樊伟*

(1. 中国农业部东海与远洋渔业资源开发利用重点实验室，上海 200090; 2. 中国水产科学研究院 渔业资源与遥感信息技术重点开放实验室，上海 200090; 3. 中国科学院大气物理研究所 东亚区域气候环境重点实验室，北京 100029)

中西太平洋延绳钓黄鳍金枪鱼渔场时空分布与温跃层关系

杨胜龙1，2，张忭忭1，靳少非3，樊伟1*

(1. 中国农业部东海与远洋渔业资源开发利用重点实验室，上海 200090; 2. 中国水产科学研究院 渔业资源与遥感信息技术重点开放实验室，上海 200090; 3. 中国科学院大气物理研究所 东亚区域气候环境重点实验室，北京 100029)

为了解热带中西太平洋延绳钓黄鳍金枪鱼(Thunnusalbacares)适宜的温跃层参数分布区间，采用Argo浮标温度信息和中西太平洋渔业委员会(The Western and Central Pacific Fisheries Commission，WCPFC)的黄鳍金枪鱼延绳钓渔获数据，绘制了热带中西太平洋月平均温跃层特征参数和月平均CPUE的空间叠加图，用于分析热带中西太平洋黄鳍金枪鱼中心渔场时空分布和温跃层特征参数间的关系。分析结果表明：热带中西太平洋温跃层上界深度、温度具有明显的季节性变化，而温跃层下界深度、温度季节性变化不明显，黄鳍金枪鱼中心渔场分布和温跃层季节性变化有关。全年中心渔场的位置分布在温跃层上界深度高值区域，随温跃层上界深度高值区域季节性南北移动。在新几内亚以东纬向区域(5°N～10°S，150°E～170°W)上界深度值全年都在70～100 m之间，全年都是延绳钓黄鳍金枪鱼中心渔场。中心渔场上界温度多在26℃以上，但是在上界温度超过30℃区域，CPUE值较小。中心渔场主要分布在温跃层下界深度两条高值带之间区域，在温跃层下界深度超过300 m和小于150 m区域，CPUE值均偏低。中心渔场主要分布在下界温度低于13℃区域，下界温度超过17℃难以形成中心渔场。频次分析和经验累积分布函数计算其适宜温跃层特征参数分布，得出中西太平洋黄鳍金枪鱼适宜的温跃层上界温度和深度分别是27～29.9℃和70～109 m；适宜的温跃层下界温度和深度分别是11～13.9℃和250～299 m 。文章初步得出中西太平洋黄鳍金枪鱼中心渔场温跃层各特征参数的适宜分布区间及季节变化特征，为我国金枪鱼实际生产作业提供技术支持。

黄鳍金枪鱼；中西太平洋；温跃层；Argo

1 引言

黄鳍金枪鱼(Thunnusalbacares)具有较高的商业价值，是我国远洋金枪鱼渔业的主要捕捞对象之一。标志放流和声学遥测等研究表明黄鳍金枪鱼具有高速游泳的能力，白天下潜到很深的水域觅食深水散射层生物(DSL)[1—7]，水温的垂直结构，尤其是温跃层分布特征在黄鳍金枪鱼渔场的形成中是极为重要的因素[8—10]。Zagaglia等[8]认为黄鳍金枪鱼这种高速移动，尤其是垂直方向的远涉会减少金枪鱼捕捞和海表温度关系。Lan等[9]研究表明相比海表温度，大西洋黄鳍金枪鱼延绳钓单位捕捞努力量渔获量(Catch per unit effort CPUE)和次表层的水温关系更密切，并推断较高的次表层水温会导致温跃层垂直分布更深，从而产生较高的黄鳍金枪鱼延绳钓CPUE。Maury等[10]研究认为，温跃层深度越深，日本黄鳍金枪鱼延绳钓CPUE值越大，而围网CPUE与温跃层深度分布成反比，与温跃层强度分布成正比。上述研究结果表明，温跃层分布对黄鳍金枪鱼的垂直分布有直接影响。中西太平洋在世界各渔区中金枪鱼产量最高，占全球产量近一半，也是我国金枪鱼主要作业渔区。因此分析中西太平洋温跃层时空分布，并结合延绳钓数据得出黄鳍金枪鱼适宜的垂直和水平空间分布非常重要。国内有关中西太平洋围网黄鳍金枪鱼渔场分布与海洋环境关系报道较多，而延绳钓黄鳍金枪鱼较少，而温跃层分布及其和黄鳍金枪鱼中心渔场分布关系国内外没有研究过。本文采用Argo浮标数据，绘制热带中西太平洋温跃层特征参数，分析热带中西太平洋温跃层时空分布和延绳钓黄鳍金枪鱼中心渔场时空分布关系，通过数值方法计算延绳钓黄鳍金枪鱼适宜的温跃层特征参数范围，为金枪鱼实际生产作业提供理论参考。

2 材料与方法

2.1 研究区域

黄鳍金枪鱼是暖水性鱼类，延绳钓作业主要在热带区域，选择(25°S～25°N，130°E～130°W)作为研究区域(见图1)。

2.2 研究数据

2.2.1 Argo浮标数据

Argo是“全球海洋观测网”计划，该计划设想在全球大洋中每隔3个经纬度布放一个卫星跟踪浮标，组成一个由3 000个浮标构成的庞大的Argo全球海洋观测网，至2007年正式完成[11]。因此本文采用2007—2011年Argo浮标剖面水温和深度数据进行温跃层特征分析，数据来自于中国Argo实时数据中心(http://www.argo.org.cn/)。

2.2.2 渔获量数据

采用中西太平洋渔业委员会(the Western and Central Pacific Fisheries Commission,WCPFC)2007—2011年按年、月份统计的黄鳍金枪鱼渔获量和捕捞努力量数据，空间分辨率采用国际金枪鱼官方统计精度5°×5°。渔业数据有投放钩数、渔获产量、渔获尾数和作业日期、地点(经度、纬度)等参数。按5°×5°统计方格内的渔获率(CPUE，单位：尾/千钩)计算公式为：

(1)

式中，CPUE(i,j)，Nfish(i,j)，Nhook(i,j)分别是第i个经度、第j个纬度处方格的月平均CPUE，月总渔获尾数和月总投钩数。公式(1)可以消除投影后低纬度和高纬度网格大小不同带来的影响[12]。

2.3 研究方法

2.3.1 温跃层等值线数据计算

采用Zhou等[13]的温跃层判别方法，取大洋温跃层强度最低标准值为0.05℃/m，对温度剖面逐层判断。把连续满足跃层标准的作为一个跃层段；对不连续者，如果跃层段之间的间隔小于10 m(当上界深度小于50 m)或小于30 m(当上界深度大于50 m)，则将两段合并进行跃层标准值判定。合并后，如果温度梯度大于或等于标准值，则合并为一个温跃层段；否则以上界深度50 m为界，分别在50 m以浅、以深，选取跃层强度强者，如强度相等，则选跃层厚度厚者为温跃层段。要求合并后的跃层厚度不小于10 m(当上界小于50 m)或不小于20 m(上界大于50 m)。具体的网格化计算参考文献[14]。按照Levitus[15]的季节划分，将北半球的季节划分如下：1—3月为冬季，4—6月为春季，7—9月为夏季，10—12月为秋季。

2.3.2 黄鳍金枪鱼中心渔场定义

按公式(1)计算2007—2011年各月网格内CPUE月平均值，共1 839个。计算1 839个CPUE的平均值、均方差和四分位数(Q1～Q4)。大于Q3的CPUE称为高值CPUE，即认为CPUE较高，而其所属渔区定义为黄鳍金枪鱼中心渔场。

2.3.3 适宜温跃层参数空间分析

把CUPE数据按月分别和温跃层上界深度、温度和温跃层下界深度、温度进行匹配，在空间上进行数据叠加，绘制CPUE和温跃层特征参数空间叠加后的月分布图，并分析CPUE、温跃层特征参数时空分布特征。最后定量分析黄鳍金枪鱼渔场和温跃层特征参数关系，找出黄鳍金枪鱼中心渔场温跃层参数变化范围。

2.3.4 适宜温跃层参数数值计算

黄鳍金枪鱼最适温跃层特征参数分别通过频次分析和经验累积分布函数(ECDF)得到[16]。计算与中心渔场CPUE对应温跃层特征参数的平均值和均方差，及适宜温跃层特征参数区间；采用非数统计K-S(Kolmogorov-Smirnov)检验方法进行显著性检验，同时采用K-S检验，检验高值CPUE对应的4个温跃层变量是否为正态分布。

2.3.5 软件工具

空间分析、数值计算和图片绘制均采用Matlab2010(a)软件及相关工具包。

3 结果与分析

3.1 CPUE时空分布

2007—2011年中西太平洋黄鳍金枪鱼延绳钓月平均CPUE均值是2.34尾/千钩(SD=2.86，n=1 839)，Q3是3.66尾/千钩。从空间上看，中西太平洋延绳钓黄鳍金枪鱼中心渔场呈现集中到离散，再定集中的空间变化趋势。4-8月份中心渔场分布比较集中，分布在澳大利亚和新几内亚周围外海，这段时期正是太平洋黄鳍金枪鱼捕捞盛渔期。从9月份开始，渔场开始向四周扩散，中心渔场分布零散，到第二年3月份再集中。全年在新几内亚外海CPUE都很高。

3.2 温跃层上界深度

中西太平洋温跃层上界深度月平均空间分布图(图1)表明，温跃层月平均上界深度大致呈纬向带状分布，且在南北半球存在明显不同；而且其季节变化特征也十分明显，呈冬深、夏浅的分布趋势。1—4月份，大致在10°S以北纬向大部分海域，温跃层上界深度超过90 m；而在10°S以南海域，其上界深度却均小于60 m，CPUE值几乎都小于Q3。7—10月份恰好相反，5°N以南海域的上界深度大于90 m，而5°N以北海域的上界深度小于60 m。中间5、6月和11、12月是中西太平洋海域的季节转换月份，5、6月温跃层月平均上界深度在90 m，11、12月温跃层月平均上界深度在100 m。在新几内亚以东纬向区域(5°N～10°S，150°E～170°W)上界深度值全年都在70 m以上，1—3月份最深达110 m。上述区域CPUE值全年都较高，全年都是延绳钓黄鳍金枪鱼中心渔场。在澳大利亚外海(10°S以南)在3—7月份有季节性中心渔场出现，4—6月份最多，中心渔场对应的温跃层上界深度值变化大，从3月份的50 m深至7月份100 m。从7月份开始，赤道以南区域温跃层上界深度值变大以后，上述区域中西渔场逐渐消失，但仍有零星渔场出现。在新几内亚外海(5°N～10°S，130°E～150°W)全年也有中心渔场出现，在6—9月份CPUE值达到全年最高，其他月份CPUE值少有回落。上述区域温跃层上界深度值全年变化不大，在70～90 m之间。在15°N以北三角区域，温跃层上界深度从3月份开始变小，4月份上述区域大部分地方上界深度值低于40 m，到9月份深度值逐渐上升。上述区域全年CPUE值都小，尤其在深度值低的月份，CPUE更小。

图1 温跃层上界深度和CPUE空间叠加图Fig.1 The overlay map of CPUE and the upper boundary depth of thermocline

3.3 温跃层上界温度

中西太平洋温跃层上界温度月平均空间分布图(图2)同样表现出季节性变化特征。1—4月份在15°N以北区域，温跃层上界温度出现低温区域(小于26℃)，5月份开始，温跃层上界温度高温区域向北移动，15°S以南区域，温跃层上界温度低于26℃。全年黄鳍金枪鱼中心渔场高值区域，温跃层上界温度多在26℃以上，低于25℃，CPUE普遍较小。但是在上界温度超过30℃区域，CPUE多小于Q3。 7月份到12月份，15°～25°S澳大利亚外海，温跃层上界温度下降到25℃以下，中心渔场渐渐消失，CUPE值变小。从3月份开始，上述区域上界温度上升到26℃以上，重新开始形成中心渔场。在新几内亚外海及以东纬向区域，CPUE值全年都很高，上界温度在29℃左右。

图2 温跃层上界温度和CPUE空间叠加图Fig.2 The overlay map of CPUE and the upper boundary temperature of thermocline

图3 温跃层下界深度和CPUE空间叠加图Fig.3 The overlay map of CPUE and the lower boundary depth of thermocline

图4 温跃层下界温度和CPUE空间叠加图Fig.4 The overlay map of CPUE and the lower boundary temperature of thermocline

3.4 温跃层下界深度

中西太平洋温跃层下界深度月平均空间分布图(见图3)表明，在10°～20°N，0°～20°S纬度带是跃层最深的区域，且从西到东两条高值带逐渐向高纬移动，跃层最深可达350 m 以深。在0°～10°N ，跃层深度相对较浅，从西太平洋向东逐渐向高纬拓展，其深度为200 m 左右。在热带跃层深度高值以外至20°为跃层深度较浅区域。全年在赤道北部的高值区域，下界深度浅于250 m，虽有部分区域有中心渔场形成，但CPUE值普遍偏低。在中心渔场集中的新几内亚外海及以东纬向区域，温跃层下界深度常年在250～300 m。在研究区域东部，跃层下界深度较浅，该区域没有渔获数据。

3.5 温跃层下界温度

图4呈现了研究海域温跃层月平均下界温度的分布。可以看到，温跃层下界温度分布似乎没有明显的季节性差异，温度分布在10～20℃之间。主要分布特征为，15°N以北和15°S以南海域，温跃层下界温度值相对较大；期间(即在15°N和15°S之间的纬向带内)呈现一片低值区，且从西到东两条高值带逐渐向高纬移动。受纬向区域影响，中西渔场对应的温跃层下界温度跨度较大。在澳大利亚外海(10°S以南)3—7月份季节性中心渔场区域，温跃层下界温度较高，在16～19℃。在新几内亚外海及以东纬向区域，温跃层下界温度低于13℃。

3.6 中心渔场适宜温跃层参数范围

2007—2011年黄鳍金枪鱼延绳钓中心渔场所在区域，上界温度分布在22～29.9℃之间(见图5a)，83.88%的中心渔场CPUE分布在27～29.9℃之间，中心渔场CPUE趋向于集中在28～28.9℃。中心渔场的下界温度分布在11～19.9℃之间(见图6b)，67.9%的中心渔场分布在11～13.9℃之间，中心渔场CPUE趋向于集中在11℃(见图5b)。中心渔场的上界深度分布在30～139 m之间，73%的中心渔场分布在60～109 m之间，中心渔场趋向于集中在90 m(见图5c)。中心渔场的下界深度在150～299 m之间(图5d)，70.5的中心渔场分布在250～299 m之间。K-S单样本检验表明，高值CPUE对应的温跃层上界深度频数服从正态分布(α=0.05)。

ECDF分析结果显示检验结果是显著的(p<0.01)，表明中心渔场CPUE和温跃层上界温度、深度，以及温跃层下界深度、温度有密切关系。结果表明CPUE和4个变量关系密切的区域分别是27.75～30.25℃(29±1.25，见图7a)、10.38～15.62℃(13±2.62，见图7b)、68.11～111.89 m(90±21.89，见图7c)和237.45～302.55 m(270±32.55，见图7d)。

综合频次分析和ECDF分析结果，本文取中西太平洋黄鳍金枪鱼适宜的温跃层上界温度和深度分别是27～29.9℃和70～109 m；适宜的温跃层下界温度和深度分别是11～13.9℃和250～299 m。

4 讨论

4.1 中西太平洋延绳钓黄鳍金枪鱼空间分布与温跃层关系

标志放流和声学调查等研究表明黄鳍金枪鱼垂直分布呈现明显的昼夜不同，晚上都在混合层或温跃层以上水域，白天大部分时间在温跃层以上，但会游到温跃层一下水域觅食。如东太平洋黄鳍金枪鱼在夜晚超过93%的时间分布在温跃层(20℃等温线)以上水域，多浅于50 m[2—5]。太平洋琉球群岛附近黄鳍金枪鱼一天的大部分时间在温跃层以上区域，56%的时间在低于表层水温1℃浅水区域[1]。在中西太平洋，黄鳍金枪鱼幼鱼夜间分布都浅于100 m，比东太平洋黄鳍金枪鱼垂直分布深，全天多数时间都在温跃层之上水域，少数时间游到温跃层以下[17]。黄鳍金枪鱼是暖水性鱼类，栖息和产卵需要在一定的水温之上。如东北太平洋黄鳍金枪鱼纬向空间分布受18℃等温线影响[2-5]，中西太平洋金枪鱼围网黄鳍金枪鱼在29℃左右产量最高[18]。本文研究表明中心渔场多分布在27～29.9℃，集中在28～28.9℃，在温跃层上界温度高于30℃和低于25℃难以形成中心渔场；同时中心渔场趋向于集中在上界深度90 m区域。图1和图2中显示，上界深度低和温度值过低地方不容易形成中心渔场。图3和图4显示在温跃层下界温度高，同时下界深度值大的地方，虽有渔获，但难以形成好的中心渔场，如在15°N以北三角区域，这个和大西洋黄鳍金枪鱼空间分布类似[18]。在上述区域温跃层下界温度高，同时下界深度值大，黄鳍金枪鱼追捕的DSL也会下潜的更深的冷水区域寻求庇护使得黄鳍金枪鱼的觅食变得困难。同时，金枪鱼类是视觉和机会捕食者，深度大的地方光线很弱，这降低了黄鳍金枪鱼的捕食能力，相比新几内亚以东海域，黄鳍金枪鱼在该区域生成环境相比恶劣。此外黄鳍金枪鱼群体觅食时分布水层较深，而且垂直分布不集中，而延绳钓捕捞主要针对该群体，因而渔获率较低。

4.2 中西太平洋延绳钓黄鳍金枪鱼垂直分布和温跃层关系

黄鳍金枪鱼游动速度快，全天虽然多数分都在温跃层之上水域，但会经常下游到温跃层以下水域捕食DSL，黄鳍金枪鱼在觅食DSL时有能力突破温跃层到深层冷水区域索饵。标志放流研究表明黄鳍最深可以快速下潜到水下1 000多米[1—5]。不同的水温垂直结构，或者可以捕食的生物组成都会影响金枪鱼的水层分布模式[17]，因此不同种类、年龄和大小的金枪鱼，在不同的个体垂直水层分布、停留时间会不同，从而影响渔船捕捞效率。在中西太平洋，温跃层要垂直分布要比东太平洋分布明显更深。标志放流表明，在珊瑚海和中西太平洋，大眼金枪鱼游动的最大深度在400～450 m[20]，明显比热带东太平洋225～270 m深。东太平洋成年黄鳍金枪鱼频繁进入到150～250 m水层觅食DSL，周边水温约12℃[2—5]。在新几内亚附件幼体黄鳍金枪鱼索饵的水层深度在200～350 m，平均在200～250 m[17]。与大眼金枪鱼分布类似，黄鳍金枪鱼在太平洋西部比东太平洋成年黄鳍金枪鱼垂直分布深，这种分布趋势和温跃层深度分布相似。

印度洋黄鳍金枪鱼适宜的温跃层下界深度、温度区间是140～200 m和13～16℃[21]；大西洋黄鳍金枪鱼适宜的温跃层下界深度、温度区间是150～249 m和12～15℃[19]。本文得出的中西太平洋黄鳍金枪鱼适宜的温跃层下界温度和深度分别是11～13.9℃和250～299 m，比印度洋和大西洋要深。延绳钓调查表明，印度洋黄鳍金枪鱼高渔获率水温为16℃[22]，大西洋在13℃[23]。在太平洋，标志放流研究指出东北太平洋成年黄鳍金枪鱼频繁进入12水域℃[2-5]，中西太平洋黄鳍金枪鱼进入12℃水域，甚至可能更冷的水层[17]。杨胜龙等[17]通过比较印度洋和大西洋黄鳍金枪鱼的高渔获率水层和温跃层关系认为，影响印度洋和大西洋黄鳍金枪鱼垂直分布的环境因子是温跃层下界变量。认为因为热带大西洋和印度洋存在不同的温跃层下界深度和温度[24]，才导致了两大洋延绳钓黄鳍金枪鱼高渔获率的水层和水温各不相同。可以推测在中西太平洋黄鳍金枪鱼垂直分布水层也受温跃层下界深度和温度影响，因为中西太平洋温跃层比热带大西洋和印度洋分布深，使得该海域黄鳍金枪鱼觅食时索饵分布更深。

Schaefer等[2—5]指出，黄鳍金枪鱼有能力突破温跃层进入深水层，甚至超过1 000 m，黄鳍金枪鱼垂直分布不受温跃层影响。但同时指出，黄鳍金枪鱼频繁进入深水层是为了觅食DSL，相关的调查研究也证实这一结论。而太平洋DSL调查研究表明，DSL的垂直分布和温跃层有关，DSL白天分布在温跃层以下区域，晚上在温跃层以上。当温跃层变浅时，DSL垂直分布也变浅[25]。Marchal等[26]调查指出，大西洋声音散射层(Sound Scattering Layer SSL，和前面DSL都表示海洋中层生物)的密度和金枪鱼捕捞分布一致，夜间SSL生物分布与温跃层深度分布一致，垂直分布在温度梯度最大处。大眼金枪鱼觅食DSL下层生物，相比大眼金枪鱼，黄鳍金枪鱼因垂直游动水层浅，冷水持续游动时间段，所以黄鳍金枪鱼觅食的是DSL上层生物[27]。因此我们推断，温跃层通过影响DSL昼夜垂直分布以及水层分布，进而影响黄鳍金枪鱼昼夜垂直分布和水层分布。

4.3 环境变量的选择

本文首次采用Argo数据重构了中西太平洋温跃层特征参数，分析得出中西太平洋黄鳍金枪鱼适宜的温跃层分布区间。以往研究太平洋黄鳍金枪鱼与海洋环境因子关系，多采用海表温度(SST)，表层盐度(SSS)、叶绿素浓度a(Chla)和海面高度(SSH)[18]。中西太平洋黄鳍金枪鱼大部分时间在表层100 m以浅，但很少活动在10 m以浅[17]，所以采用温跃层温度比SST跟合理。黄鳍金枪鱼处于食物链顶部，不是直接食用Chla，而SSS被证实是对延绳钓黄鳍金枪鱼影响最弱的环境因子[10]，Arrizabalaga等采用海平面高度变异替代温跃层[28]。研究证实，相比海表表层环境变量，温跃层对延绳钓黄鳍金枪鱼CPUE影响更大，因为延绳钓捕捞的黄鳍金枪鱼分布水层深，而且这种垂直分布会减弱与海洋表层变量的关系[8—10]。因此本文采用温跃层分析黄鳍金枪鱼生态环境栖息属性比较合适，适宜的温跃层上界温度和深度分别是27～29.9℃和70～109 m；适宜的温跃层下界温度和深度分别是11～13.9℃和250～299 m。可以参考本文绘制的月平均温跃层参数空间分布图寻找中心渔场，并选择投钩的深度，以提高黄鳍金枪鱼的捕捞效率。溶解氧是黄鳍金枪鱼分布一个重要环境因子。在太平洋，溶解氧在中西部分布深，东部分布浅。溶解氧2 mL/L等深线在太平洋西部低于300 m，太平洋中部为200～300 m，太平洋东部在100 m。而在澳大利亚东海岸的珊瑚海，1.5 mL/L溶解氧深度分布在500～600 m[20]。和大西洋一样，中西太平洋溶解氧分布深，不影响黄鳍金枪鱼垂直分布[10]，本文没有研究溶解氧对黄鳍金枪鱼分布。

[1] Matsumoto T，Kitagawa T，Kimura S. Vertical behavior of juvenile yellowfin tunaThunnusalbacaresin the southwestern part of Japan based on archival tagging[J]. Fisheries Science，2013，79(3): 417-424.

[2] Schaefer K M，Fuller D W，Block B A. Movements，behavior，and habitat utilization of yellowfin tuna (Thunnusalbacares) in the northeastern Pacific Ocean，ascertained through archival tag data[J]. Marine Biology，2007，152(3): 503-525.

[3] Schaefer K M，Fuller D W，Block B A. Vertical movements and habitat utilization of skipjack (Katsuwonuspelamis)，yellowfin (Thunnusalbacares)，and bigeye (Thunnusobesus) tunas in the equatorial eastern Pacific Ocean，ascertained through archival tag data[C]//Nielsen J L，Arrizabalaga H，Fragoso N，et al. Reviews: Methods and Technologies in Fish Biology and Fisheries，Tagging and Tracking of Marine Animals with Electronic Devices. Berlin: Springer，2009，9: 121-144.

[4] Schaefer K M，Fuller D W，Block B A. Movements，behavior，and habitat utilization of yellowfin tuna (Thunnusalbacares) in the Pacific Ocean off Baja California，Mexico，determined from archival tag data analyses[J]. Fisheries Research，2011，112(1/2): 22-37.

[5] Schaefer K M，Fuller D W，Aldana G. Movements，behavior，and habitat utilization of yellowfin tuna (Thunnusalbacares) in waters surrounding the Revillagigedo Islands Archipelago Biosphere Reserve，Mexico[J]. Fisheries Oceanography，2014，23(1): 65-82.

[6] Yokawa K，Saito H，Kanaiwa M，et al. Vertical distribution pattern of CPUE of atlantic billfishes and associated species estimated using longline research data[J]. Bulletin of Marine Science，2006，79(3): 623-634.

[7] Weng K C，Stokesbury M J W，Boustany A M，et al. Habitat and behaviour of yellowfin tunaThunnusalbacaresin the Gulf of Mexico determined using pop-up satellite archival tags[J]. Journal of Fish Biology，2009，74(7): 1434-1449.

[8] Zagaglia C R，Lorenzzetti J A，Stech J L. Remote sensing data and longline catches of yellowfin tuna (Thunnusalbacares) in the equatorial Atlantic[J]. Remote Sensing of Environment，2004，93(1/2): 267-281.

[9] Lan Kuowei，Lee M A，Lu H J，et al. Ocean variations associated with fishing conditions for yellowfin tuna (Thunnusalbacares) in the equatorial Atlantic Ocean[J]. ICES Journal of Marine Science，2011，68(6): 1063-1071.

[10] Maury O，Gascuel D，Marsac F，et al. Hierarchical interpretation of nonlinear relationships linking yellowfin tuna (Thunnuslbacares) distribution to the environment in the Atlantic Ocean[J]. Canadian Journal of Fisheries and Aquatic Sciences，2001，58(3): 458-469.

[11] China Argo Real-time Data Center. Argo newsletter[R]. Hangzhou: The second Institute of Oceanography，State Oceanic Administration. 2012，28(4): 7-9.

[12] Lee P F，Chen I C，Tzeng W N. Spatial and temporal distribution patterns of bigeye tuna (Thunnusobesus) in the Indian Ocean[J]. Zoological Studies，2005，44(2): 260-270.

[13] Zhou Yanxia，Li Binglan，Zhang Yijun，et al. World Oceanic Thermocline Characteristics in Winter and Summer[J]. Marine Science Bulletin，2002，21(1): 16-22.

[14] Yang Shenglong，Ma Junjie，Wu Yumei，et al. Seasonal variability of temperature profile inThunnusobesusandThunnusalbacaresfishing ground in the Indian Ocean[J]. Marine Sciences，2012，36(7): 97-103.

[15] Levitus S. Climatological Atlas of the world ocean[J]. Eos，Transactions American Geophysical Union，1983，64(49): 962-963.

[16] Zainuddin M，Saiton K，Saiton S I. Albacore (Thunnusalalunga) fishing ground in relation to oceanographic conditions in the western North Pacific Ocean using remotely sensed satellite data[J]. Fish Oceanogr，2008，17(2): 61-73.

[17] Leroy B，Itano D G，Nicol S. Preliminary analysis and observations on the vertical behaviour of WCPO Skipjack，Yellowfin and Bigeye Tuna in association with anchored FADs，as indicated by acoustic and archival tagging data[C]. Western & Central Pacific Fisheries Commission，2007.

[18] 王少琴，许柳雄，朱国平，等. 中西太平洋金枪鱼围网的黄鳍金枪鱼CPUE时空分布及其与环境因子的关系[J]. 大连海洋大学学报，2014，29(3): 303-308.

Wang Shaoqin，Xu Liuxiong，Zhu Guoping，et al. Spatial-temporal profiles of CPUE and relations to environmental factors for yellowfin tunaThunnusalbacoresfrom purse-seine fishery in Western and Central Pacific Ocean[J]. Journal of Dalian Fisheries University，2014，29(3): 303-308.

[19] 杨胜龙，马军杰，张禹，等. 大西洋中部延绳钓黄鳍金枪鱼渔场时空分布与温跃层的关系[J]. 生态学报，2013，33(19): 6345-6353.

Yang Shenglong，Ma Junjie，Zhang Yu，et al. Relationship between the temporal-spatial distribution of longline fishing grounds of yellowfin tuna (Thunnusalbacares) and the thermocline characteristics in the Central Atlantic Ocean[J]. Acta Ecologica Sinica，2013，33(19): 6345-6353.

[20] Evans K，Langley A，Clear N P，et al. Behaviour and habitat preferences of bigeye tuna (Thunnusobesus) and their influence on longline fishery catches in the western Coral Sea[J]. Canadian Journal of Fisheries and Aquatic Sciences，2008，65(11): 2427-2443.

[21] 杨胜龙，张禹，张衡，等. 热带印度洋黄鳍金枪鱼渔场时空分布与温跃层的关系[J]. 生态学报，2012，32(3): 671-679.

Yang Shenglong，Zhang Yu，Zhang Heng，et al. The relationship between the temporal-spatial distribution of fishing ground of yellowfin tuna (Thunnusalbacares) and themocline characteristics in the tropic Indian Ocean[J]. Acta Ecologica Sinica，2012，32(3): 671-679.

[22] Mohri M，Nishida T. Consideration on distribution of adult yellowfin tuna (Thunnusalbacares) in the Indian Ocean based on Japanese tuna longline fisheries and survey information[J]. IOTC Proc，2000(3): 276-282.

[23] 宋利明，陈新军，许柳熊. 大西洋中部黄鳍金枪鱼(Thunnusalbacares)的垂直分布与有关环境因子的关系[J]. 海洋与湖沼，2004，35(1): 64-68.

Song Liming，Chen Xinjun，Xu Liuxiong. Relationship between vertical distribution of yellowfin tunas’ (Thunnusalbacares) and the concerned environmental factors in the Central Atlantic Ocean[J]. Oceanologia et Limnologia Sinica，2004，35(1): 64-68.

[24] 杨胜龙，化成君，靳少非，等. 热带印度洋黄鳍金枪鱼水平-垂直分布空间分析[J]. 海洋与湖沼，2014，45(2): 281-290．

Yang Shenglong，Hua Chengjun，Jing Shaofei，et al. Spatial analysis the horizontal and vertical distribution of yellowfin tunaThunnusalbacaresin the tropical Indian Ocean[J]. Oceanologia et Limnologia Sinica，2014，45(2): 281-290.

[25] Hazen E L，Johnston D W. Meridional patterns in the deep scattering layers and top predator distribution in the central equatorial Pacific[J]. Fisheries Oceanography，2010，19(6): 427-433.

[26] Marchal E，Gerlotto F，Stéquert B. On the relationship between scattering layer，thermal structure and tuna abundance in the Eastern Atlantic equatorial current system[J]. Oceanologica Acta，1993，16(3): 261-272.

[27] Brill R W，Bigelow K A，Musyl M K，et al. Bigeye tuna (Thunnusobesus) behavior and physiology and their relevance to stock assessments and fishery biology[J]. Collective Volume of Scientific Papers ICCAT，2005，57(2): 142-161.

[28] Arrizabalaga H，Pereira J G，Royer F，et al. Bigeye tuna (Thunnusobesus) vertical movements in the Azores Islands determined with pop-up satellite archival tags[J]. Fisheries Oceanography，2008，17(2): 74-83.

Relationship between the temporal-spatial distribution of longline fishing grounds of yellowfin tuna (Thunnusalbacares) and the thermocline characteristics in the Western and Central Pacific Ocean

Yang Shenglong1，2，Zhang Bianbian1，Jin Shaofei3，Fan Wei1

(1.KeyLaboratoryofEastChinaSea&OceanicFisheryResourcesExploitationandUtilization，MinistryofAgriculture，Shanghai200090，China;2.KeyandOpenLaboratoryofRemoteSensingInformationTechnologyinFishingResource，EastChinaSeaFisheriesResearchInstitute，Shanghai200090,China;3.KeyLaboratoryofRegionalClimate-EnvironmentforTemperateEastAsia，InstituteofAtmosphericPhysics，ChineseAcademyofSciences，Beijing100029，China)

We evaluated the isoline distribution of thermocline characteristics (the upper and lower boundary temperatures and depths) in the yellowfin tuna (Thunnusalbacares) fishing grounds in Western and Central Pacific Ocean. We plotted the thermocline characteristics contour on a spatial overlay map using data collected on a monthly basis from Argo buoys and monthly CPUE (catch per unit effort) from yellowfin tuna longlines from the Western and Central Pacific Fisheries Commission (WCPFC). In addition，frequency analysis and the empirical cumulative distribution function (ECDF) were used to calculate the optimum ranges for the thermocline characteristics of the central fishing grounds. Our analysis suggested that there were significant seasonal variations in the upper boundary temperature and depth of the thermocline in the central fishing grounds，which significantly influenced the temporal and spatial distribution of the yellowfin tuna population. However，the lower boundary temperature and depth of thermocline has little seasonal variation. The overlay maps suggest that the central fishing grounds were observed in areas where the upper boundary depth of thermocline was deeper，and south-north moved follow by thermocline. The values between 70 and 100 m in the east zonal areas of New Guinea，and the central fishing grounds were found all year. The fishing grounds distributed where the upper boundary temperature of the thermocline was higher than 26℃，but CPUE was lower than Q3 while temperature higher than 30℃. The fishing grounds located between the two high value shape of the lower boundary depth of thermocline，if the depth was more than 300 m or less than 150 m，the CPUE tended to be low. The lower boundary temperature of the thermocline in the fishing grounds was lower than 13℃ all year in the equatorial zone. Conversely，if the temperature was higher than 17℃，the hooking rates are very low. Frequency analysis and the empirical cumulative distribution function (ECDF) were used to calculate the optimum range of thermocline characteristics. The optimum depth range of the upper boundary temperature and depth were 27-29.9℃ and 70-109 m .The optimum depth range of the lower boundary and the lower boundary temperature were 250-299 m and 11-13.9℃，respectively. The results were confirmed by using the Kolmogorov-Smirnov test. A preliminarily indication of the optimum distribution interval and seasonal change characteristics of each thermocline in the longline central fishing ground of yellowfin tuna was obtained in this paper. The conclusion could provide a reference for improving the efficiency of yellowfin tuna longline fishing and aid tuna resource management in Western and Central Pacific Ocean．

Thunnusalbacares; Western and Central Pacific Ocean; thermocline; Argo

10.3969/j.issn.0253-4193.2015.06.008

2014-09-22；

2014-12-15。

上海市自然科学基金(14ZR1449900)；科技支撑计划项目(2013BAD13B01)。

杨胜龙(1982—)，男，江西省九江市人，主要从事渔场次表层环境和金枪鱼渔场变动研究。E-mail:ysl6782195@126.com

*通信作者：樊伟，研究员。E-mail:fanwee@126.com

S931.3

A

0253-4193(2015)06-0078-10

杨胜龙，张忭忭，靳少非，等. 中西太平洋延绳钓黄鳍金枪鱼渔场时空分布与温跃层关系[J]. 海洋学报，2015，37(6)：78—87，

Yang Shenglong，Zhang Bianbian，Jin Shaofei，et al. Relationship between the temporal-spatial distribution of longline fishing grounds of yellowfin tuna (Thunnusalbacares) and the thermocline characteristics in the Western and Central Pacific Ocean[J]. Haiyang Xuebao，2015，37(6)：78—87，doi:10.3969/j.issn.0253-4193.2015.06.008