气候的变化、异常以及振荡驱动着海洋环境因子的变化，从而影响海洋生态系统的生命活动及过程^[1]。厄尔尼诺南方涛动 (El Niño Southern Oscillation, ENSO) (包括厄尔尼诺和拉尼娜事件) 起源于赤道太平洋，是太平洋海域气候变化的最强信号，具有年际变动周期^[2]，其基本特征是西太平洋暖池的地带性位移以及大气对流，导致了整个太平洋海盆物理性质和生态系统的改变；具体表现为对各种海洋环境参数的影响，包括海表面温度 (Sea surface temperature, SST)、叶绿素a浓度 (Chl-a) 等，以及对由海洋顶级捕食者 (如金枪鱼) 共同构成的相互关联的海洋生态系统的影响^[1,3-4]。

鲣 (Katsuwonus pelamis) 广泛分布于各大洋的热带、亚热带海域^[5]，中西太平洋热带海域是世界上最大的鲣渔场^[6-7]。因而，研究中西太平洋鲣资源的时空分布对掌握渔业生产规律和开发鲣资源尤为重要。鲣种群按集群特点进行划分，可分为自由鱼群 (Free swimming school, FSC) (亦称起水群)、流木鱼群 (Log school) 以及鲸豚附随群 (Marine mammal associated school) 3类^[8]；流木鱼群中亦包括人工集鱼装置 (Fish aggregation device, FAD) 鱼群。自20世纪80年代以来，FAD开始被广泛运用于鲣的围网作业中，是一种高效的捕捞模式^[9]，但捕捞副渔获物 [非目标种类的幼年大眼金枪鱼 (Thunnus obesus) 以及黄鳍金枪鱼 (T. albacares)等] 严重的问题随之而来^[10]。为此，中西太平洋海域的瑙鲁协定方 (Parties to nauru agreement, PNA) 制定了相关海洋管理政策以限制和减少FAD的使用^[11]。相较于前者，自由鱼群的个体较大，偏好聚集于浅水层^[12]，且围网渔获物中兼捕物种少，但由于其高速游动的特性，捕捞难度较大^[13]。因此，掌握围网鲣自由鱼群的时空分布特性，对于集中捕捉该种群以降低对漂浮物鱼群捕捉的依赖，从而并提高渔业生产效率具有积极作用，对保护金枪鱼的栖息环境与种群结构也具有重要意义^[14]。

海洋鱼类资源的时空分布极易受海洋环境因子影响，掌握其时空分布有利于鱼类资源的合理开发和利用。鲣的分布模式受海表面温度影响显著^[15]，已有研究普遍发现鲣资源主要分布在海表温度28~30 ℃的海域^[16-17]，因而该温度区间可作为研究鲣资源时空分布变动的一项参考指标^[18]。Lehodey等 ^[19]研究发现鲣作业渔场随暖池边缘(即与29 ℃等温线重合)发生了空间偏移。同时，气候变化在时空尺度上影响着海洋鱼类的分布模式，李政纬^[17]指出，29 ℃等温线东界会受厄尔尼诺与南方涛动影响，进而影响鲣围网渔场的经向分布。受ENSO影响，鲣渔场与暖池的时空分布变动具有相同趋势^[19-20]：即厄尔尼诺事件时，渔场重心随暖池东移；拉尼娜事件时，渔场重心随暖池西移^[21]。因此，将暖池变动作为研究围网鲣自由鱼群时空分布的参考指标，能够探究鲣渔场的时空分布特性，更好地进行集中、高效的渔业生产活动，为其资源开发提供科学依据。已有研究仅对历史渔场重心与极端气候指标(如南方涛动指数、海表温异常指数)的关系进行了探讨^[22-24]，而忽略了暖池相关物理海洋指标空间变动而导致的渔场重心变化，其中包括暖池重心处多出现高产量的鲣围网渔场^[25]，以及暖池右边缘经向扩展而具有的潜在东部栖息地^[26]。本研究以29 ℃等温线作为暖池边界，针对围网鲣自由鱼群的渔场重心与暖池相关指标进行相关性分析，并基于不同的气候模式下探究其时空分布变动之间的关系。

1.   材料与方法

1.1   数据来源

中西太平洋围网鲣渔业生产数据来源于中西太平洋渔业委员会 (Western and Central Fisheries Commission,WCPFC)，数据包括作业日期、作业位置 (经纬度)、捕捞努力量、自由鱼群的渔获量等信息。其中时间跨度为1995—2019年，空间范围为120°E—150°W、20°S—20°N，空间分辨率为5°×5°，样本数共计7 721个。海洋环境数据为SST，来自哥伦比亚大学气候数据实验室网站 (http://iridl.ldeo.columbia.edu)，空间分辨率为0.5º×0.5°，时间分辨率为月。同时选取1995—2019 年美国国家海洋与大气管理局气候预测中心 (https://origin.cpc.ncep.noaa.gov) 的Nino 3.4区域 (170°W—120°W、5°N—5°S) 的平均海温距平作为海洋尼诺指数 (OceanicNiño Index, ONI)，其中ONI值介于 ±0.5为正常气候模式，大于0.5发生厄尔尼诺事件，小于 −0.5发生拉尼娜事件。

1.2   研究方法

1.2.1   重心计算

渔业研究中，资源丰度指数通常用于量化渔业中种群丰度的时空变化，如单位捕捞努力量渔获量 (Catch per unit effort, CPUE) 和渔获量^[27]。由于现代商业渔业的复杂性，常用的CPUE可能无法作为所有物种资源丰度指数的合理表征^[28]，因此将二者均纳入考虑来量化围网鲣自由鱼群的时空分布，进而与暖池相关指标进行比较分析，为围网鲣自由鱼群的渔业生产和资源开发提供借鉴。渔获量大小常用于直接表示渔业的资源丰度，而CPUE大小常被认为与渔业资源丰度成正比^[29]，因而可作为表示其局部资源丰度的指标^[30]。

本研究对1995—2019年间的渔场重心进行分析计算，作为研究其渔场时空分布变化的描述标准，通过使用渔场重心法来描述中西太平洋围网鲣自由鱼群 (以下简称自由鱼群) 的时空分布变化，环境因子SST的重心计算同理。由于鲣渔场与暖池变动间多为经向变化，因而仅对经向维度进行考虑，渔场重心计算范围与原始数据的空间范围 (120°E—150°W、20°S—20°N) 一致，同时为更准确地说明中西太平洋暖池的重心变化，在计算各项指标的重心时，本研究中的海域范围聚焦于110°E—80°W、45°S—45°N，公式如下^[31]：

式中：G_CPUE、G_Catch分别为自由鱼群渔场两种资源丰度指标的重心经度 (以下简称CPUE重心和Catch重心)；G_SST为 SST重心的经度 (以下简称暖池重心经度)；$ {{E}}_{i} $为单位捕捞努力量渔获量；$ {{C}}_{i} $为捕捞产量渔获重心；$ i $ 表示经度，$ {{L}}_{i} $为第 $ i $ 经度的重心经度；S_i为第i经度的温度重心。

暖池的右边缘通过自中西太平洋西部海盆起选取连续的29 ℃等温线右边缘的经度，记作R_SST。

1.2.2   统计方法

采用皮尔森相关性分析^[32]，对暖池指标 (暖池重心经度、暖池右边缘经度) 与渔场资源丰度重心经度以及气候环境指标ONI指数进行相关性分析。相关性在0.1~0.3为弱相关性，0.3~0.5间为中度相关性，大于0.5为强相关性；P<0.05表示相关，P<0.01为显著相关，P<0.001为极显著相关，P>0.05则表示相关性不显著。

1.2.3   不同气候模式下鲣资源空间分布

建立不同气候模式 (正常气候模式、厄尔尼诺年份、拉尼娜年份等) 下自由鱼群与暖池间的时空分布图，进一步说明不同气候模式下二者之间的空间相对位置及关系。

2.   结果

2.1   气候变化下暖池指标与鲣资源重心的关系

通过ONI指数对气候变化进行量化，并且采用皮尔森检验对暖池相关指标重心及右边缘、渔场资源丰度重心以及ONI指数进行相关性检验，发现其在统计学上均呈极显著相关性 (图1)。结果显示，渔场资源丰度重心指标，包括渔获量经向重心G_Catch与单位捕捞努力量渔获量经向重心G_CPUE，二者与暖池经向重心呈中度相关性，其中G_CPUE的相关性值略高于G_Catch。而在与暖池右边缘进行相关性检验时，G_Catch与G_CPUE同其呈弱相关性，G_CPUE表现仍稍高于G_Catch。在渔场资源丰度重心指标与海洋尼诺指数的相关性检验上，G_Catch、G_CPUE与之均呈中度相关性，此时G_Catch表现略高于G_CPUE。除主要研究渔场经向重心外，将暖池指标与海洋尼诺指数进行相关性检验，结果表明暖池经向重心G_SST和暖池右边缘R_SST与之呈强相关性，其中G_SST与之相关性值达0.79，R_SST为0.68。所有指标间的P均小于0.001，表现为极显著相关性。

图  1  皮尔森相关性检验

Figure  1.  Pearson correlation test

下载: 全尺寸图片 幻灯片

资源丰度的表征在渔业研究中尤为重要^[28]，经皮尔森相关性检验得出，CPUE与暖池指标间在相关性表现上相较Catch更优。此外，CPUE的计算避免了直接使用渔获量，而未纳入捕捞努力量的影响^[33]；因此，可选用CPUE作为围网鲣自由鱼群与暖池指标关系研究的资源丰度指标。

2.2   暖池指标与鲣资源重心的年际变化趋势

将各项指标的经向变化随年月变化进行对比，探究在不同气候条件下，渔场重心经向与暖池指标变动的时空分布变化 (图2)。结果表明，渔场重心指标间的经向年际变化基本一致，G_CPUE与暖池指标的相关性表现优于G_Catch，因而以下研究均以G_CPUE的经向变动进行探究。各项重心经向指标大多集中分布于140°E 以东、180°以西海域，其中G_CPUE与G_SST主要分布在160°E经线两端，而R_SST大多分布在180°以东 (图2-a)。G_CPUE、G_SST与R_SST随海洋尼诺指数的变动而发生变化 (图2-b)。多数情况下，当发生厄尔尼诺事件 (ONI>0.5)时，G_CPUE与G_SST东移越过160°E，同时R_SST向东移动，如1997年的6—10月等；发生拉尼娜事件 (ONI<−0.5)时，G_CPUE 与G_SST则分布在160°E以西海域，同时R_SST向西移动，如1995年8月—1996年3月等。暖池指标G_SST与R_SST间的经向变动关系随时间变化基本呈相同趋势。

图  2  各指标值随时间的变化

Figure  2.  Changes in values of each indicator along with time

下载: 全尺寸图片 幻灯片

G_CPUE与G_SST在经向上的变动趋势不完全相同。1995年1—4月、2003年1—4月、2007年1—5月和2016年1—5月等均为G_SST在G_CPUE以东，且经向变动趋势显著相反，而此时均为由发生厄尔尼诺事件向正常气候模式 (−0.5<ONI<0.5) 转变的月份。当G_SST在G_CPUE以西且经向变动趋势相反时，如1998年6—12月、1999年5—12月和a2010年6—12月等，此时正发生强烈的拉尼娜事件，受厄尔尼诺事件影响，G_SST多位于G_CPUE以东，但此时经向的变动趋势基本相同，如1997年6月—1998年4月、2015年6—12月等。G_CPUE与R_SST之间关系和G_SST一致。

2.3   不同气候模式下自由鱼群的分布

为探明上述不同气候模式下，自由鱼群与暖池间的时空分布关系，选取其中个别月建立时空分布 (图3)。拉尼娜事件选取2010年的11—12月，正常气候模式选取2013年11—12月，厄尔尼诺事件选取2015年11—12月，而厄尔尼诺事件转变为正常年份事件的特殊气候模式选取2016年4—5月。

图  3  不同气候模式下CPUE与暖池的时空分布

注：a、b为拉尼娜气候模式；c、d为正常气候模式；e、f为厄尔尼诺气候模式；g、h为厄尔尼诺气候模式向正常气候模式转变。

Figure  3.  Spatial-temporal distribution of CPUE and warm pool under different climate modes

Note: (a) and (b) are the La Niña climate modes; (c) and (d) are the normal climate modes; (e) and (f) are the El Niño climate modes; (g) and (h) are the climate modes shifting from the El Niño to the normal.

下载: 全尺寸图片 幻灯片

与ENSO相关的指标的位移发生在整个中西太平洋，发生拉尼娜事件时 (2010年11—12月)，暖池面积最小，右边缘东部抵达距离最短，接近于180°经线，此时CPUE的分布均位于180°经线以西，且集中于太平洋西部美拉尼西亚群岛海域。当正常年份时 (2013年11—12月)，暖池右边缘的位置相较拉尼娜事件时向东偏移，CPUE越过180°以东且少量分布。发生厄尔尼诺事件时 (2015年11—12月)，CPUE的分布随暖池右边缘位置的东移而向东扩展。此外，一个特殊的气候模式变化应当被纳入考虑，即厄尔尼诺事件转变为正常年份事件 (2016年4—5月)，此时暖池右边缘向西收缩，但CPUE的分布仍向东扩展，而相较于厄尔尼诺事件时在热带海域的均匀分布，此时仅在180°以东少量分布，且CPUE的大小和数量均下降。在所有的气候模式下，暖池右边缘均位于自由鱼群以东位置，自由鱼群基本均匀分布于暖池范围内，也同样说明了二者之间的紧密联系。

3.   讨论

不同类型的ENSO事件通过驱动暖池的空间位移从而改变区域内的海洋环境因子，进一步驱动海洋生态系统及其内部物种资源的响应^[34]。 位于中西太平洋的暖池是ENSO的基本构成要素^[15]，其表面积及位置随ENSO事件的变化而改变^[13,35]；即厄尔尼诺事件时，暖池的面积增大，东移到达太平洋中部海域；拉尼娜事件时，暖池的面积减小，西移到西太平洋海盆区域。从暖池指标(本研究中指暖池重心与暖池右边缘)的时空变化入手，可为研究围网鲣自由群对ENSO事件的响应提供依据。鲣的资源丰度以及空间分布受海洋环境因子影响^[36]，而海表温度对鲣资源状态的时空分布具有更显著的影响^[37]，在本研究中具体表现为资源丰度的大小及重心的空间分布随29 ℃海表温度场的影响而变化。由于ENSO驱动海表温度变化，受其影响，具有高度洄游的鲣自由鱼群会在沿海生态系统和公海海域之间移动^[38]。Williams和Ruaia^[39]指出，ENSO事件对鲣围网活动空间分布的影响表现为，在厄尔尼诺年份，捕捞活动通常进一步向东扩展，在拉尼娜期间，捕捞活动向西部地区收缩，与本研究的发现(图3)基本一致。捕捞作业活动的变动同暖池变动规律相同，而渔场重心可以表征捕捞作业活动空间的集中分布程度，但渔场重心与暖池间的协同关系也存在着特异性。受ENSO影响，暖池重心的空间变异程度 (变异系数Coefficient of variation，CV= 5.012) 要强于渔场重心 (CV = 2.523)的变化 (图2)，即表现为厄尔尼诺事件时暖池重心位于渔场重心东部，而当拉尼娜事件时暖池重心位于渔场重心西部。这可能是由于受捕捞作业方式影响，围网鲣自由鱼群的渔场重心相较于暖池受气候以及环境变化影响小^[40]。此外，自由鱼群多是体型较大的成年鲣，需寻求适宜的温度环境作为首选栖息地和合适的产卵地，厄尔尼诺事件会导致海表温度升高，温跃层变浅(温跃层是富含营养物质的混合层上下界之间交换的障碍)，营养物质更容易向混合层转移，从而促使深层营养盐在海表扩散，形成具有良好饵料的栖息环境^[41]。然而强烈厄尔尼诺事件时，温度过高导致暖池的长距离东扩，抑制了鲣在热带中西太平洋的产卵活动^[21]；但并未导致渔场范围相应东扩，反而此时在太平洋中部的岛国 (如基里巴斯群岛) 的围网渔获量更高，分布更集中。

而对于厄尔尼诺转变为正常年份这一特殊气候模式的变化，即暖池重心与渔场重心的变化趋势显著相反 (暖池重心向西移动而渔场重心仍向东)，可能由于暖池场的西移，导致了适宜的温度范围增大，以便于鲣适宜栖息地的形成。陈洋洋和陈新军 ^[23]研究得出，Nino3.4海区指数对CPUE的影响滞后0~2个月，这可能也是造成渔场重心和暖池重心之间变化不同步的原因。此外，暖池的东扩也与温跃层变浅，以及更强于平时的西太平洋信风有关，从而导致赤道西太平洋的初级生产力增加^[42]。因此，位于140°E—160°E的所罗门群岛和巴布亚新几内亚的资源丰度在厄尔尼诺现象结束后有所增加，以应对鲣捕捞量的增加以及栖息地收缩的现象，这也进一步解释了渔场重心与暖池重心之间变化的不同步。

本研究对围网鲣自由鱼群的渔场指标，以及暖池指标在不同气候模式下的时空分布变化进行探究发现，通过研究暖池重心的变化可以很好地探索及预测渔场重心变化的规律。而通过构建暖池场与自由鱼群资源丰度的时空分布关系发现，暖池右边缘的范围能够很好地与自由鱼群的空间分布联系起来，可为商业性捕捞围网鲣自由鱼群提供渔场边界的指示，也可为今后研究鲣渔场的分布范围提供一项有利参考。

本研究对不同资源丰度指标对围网鲣自由鱼群的表征的影响进行了研究，但选取的均为名义上的资源丰度指标 (名义的CPUE)，而未对CPUE进行标准化 (以消除人为捕捞等因素对资源丰度的影响)，可能会造成研究结果具有一定误差。且仅研究了自由鱼群对暖池指标的响应关系，未与由人工集鱼装置捕捞获得的流木鱼群进行对比。在未来研究中，将对不同集群种类的鲣对暖池响应的差异性进行更细致的研究，以期为中西太平洋的鲣围网作业提供更丰富详尽的科学依据。