西北太平洋柔魚(yú)漁場(chǎng)分布與渦動(dòng)能變化的相關(guān)關(guān)系

劉瑜，鄭全安，李曉峰*

( 1. 上海海洋大學(xué) 海洋科學(xué)學(xué)院，上海 201306；2. 馬里蘭大學(xué) 大氣海洋科學(xué)系，馬里蘭州 科利吉帕科 20742)

1 引言

柔魚(yú)(Ommastrephes bartramii)為我國(guó)西北太平洋公海重要經(jīng)濟(jì)魚(yú)類(lèi)之一，通常隨著黑潮進(jìn)行南北索餌、繁殖洄游，并在黑潮與親潮交匯區(qū)形成漁場(chǎng)[1]。在交匯區(qū)冷暖水匯流，海水垂直運(yùn)動(dòng)加強(qiáng)，營(yíng)養(yǎng)鹽富集并伴有浮游植物大量繁殖，從而成為高生產(chǎn)力區(qū)，吸引浮游動(dòng)物和魚(yú)類(lèi)來(lái)此索餌與繁衍[2]。黑潮延伸體是北太平洋中尺度渦動(dòng)能（Eddy Kinetic Energy，EKE）最大的海域。研究表明黑潮延伸體大彎曲以及高生產(chǎn)力區(qū)域均與渦旋活動(dòng)有相互聯(lián)系[3–4]。渦旋運(yùn)動(dòng)過(guò)程會(huì)影響環(huán)流、海溫及葉綠素的垂直和水平分布，從而影響資源豐度和漁場(chǎng)分布。以往研究多為黑潮形態(tài)或者黑潮流量對(duì)柔魚(yú)漁場(chǎng)分布的影響[5–6]，或者將渦動(dòng)能作為海洋環(huán)境要素之一，與其他環(huán)境因子如海表溫度、葉綠素、海面高度等建立模型預(yù)測(cè)漁業(yè)資源豐度與漁場(chǎng)分布[7–8]，而黑潮延伸體渦動(dòng)能對(duì)于柔魚(yú)漁場(chǎng)分布的影響研究闕如。因此，本文根據(jù)中國(guó)魷釣船在西北太平洋海域的柔魚(yú)生產(chǎn)數(shù)據(jù)和衛(wèi)星高度計(jì)遙感數(shù)據(jù)，分析柔魚(yú)漁場(chǎng)的渦動(dòng)能與柔魚(yú)的單位捕撈努力量漁獲量（Catch Per Unit Effort，CPUE）的相關(guān)關(guān)系，同時(shí)計(jì)算黑潮延伸體渦動(dòng)能的時(shí)空變化及其對(duì)柔魚(yú)漁場(chǎng)分布的影響，以期尋找柔魚(yú)漁場(chǎng)適宜的渦動(dòng)能范圍，探索西北太平洋柔魚(yú)漁場(chǎng)分布對(duì)黑潮延伸體渦動(dòng)能變化的響應(yīng)，為該漁業(yè)資源的評(píng)估與科學(xué)開(kāi)發(fā)管理提供依據(jù)。

2 材料與方法

2.1 數(shù)據(jù)來(lái)源和預(yù)處理

柔魚(yú)生產(chǎn)數(shù)據(jù)取自中國(guó)遠(yuǎn)洋漁業(yè)協(xié)會(huì)上海海洋大學(xué)魷釣技術(shù)組，數(shù)據(jù)內(nèi)容包括作業(yè)日期、經(jīng)度、緯度、日產(chǎn)量（t）和作業(yè)船數(shù)等，數(shù)據(jù)時(shí)間為2010–2016年的7–11 月，漁船作業(yè)范圍主要分布在38°～46°N，147°～170°E。CPUE 定義為每艘船每天的捕撈產(chǎn)量，單位為t/d，作為漁業(yè)資源豐度的指標(biāo)。采用漁場(chǎng)分布重心作為資源量的空間分布和變化指標(biāo)[9]，計(jì)算公式為

式中，X、Y 分別為漁場(chǎng)作業(yè)重心的經(jīng)度和緯度；Ci為作業(yè)點(diǎn)i 的CPUE；Xi和Yi分別為作業(yè)點(diǎn)i 的緯度和經(jīng)度；n 為作業(yè)總次數(shù)。

環(huán)境數(shù)據(jù)為衛(wèi)星高度計(jì)測(cè)得的海表面高度異常(Maps Sea Level Anomaly，MSLA)數(shù)據(jù)，取自法國(guó)空間局AVISO 網(wǎng)站（http://www.aviso.oceanobs.com），其中包括海面高度數(shù)據(jù)和地轉(zhuǎn)流數(shù)據(jù)。該數(shù)據(jù)采用墨卡托投影方式插值，空間分辨率（1/4）° × （1/4）°，時(shí)間分辨率為天。渦動(dòng)能的計(jì)算公式為[10]：

式中，EKE 為渦動(dòng)能，單位為cm2/s2；=u?=v?和分別是緯向速度u 和經(jīng)向速度 v 的年平均值。

2.2 研究方法

如圖1 所示，西北太平洋柔魚(yú)漁場(chǎng)主要分布在黑潮延伸體以北位置。本文以A 區(qū)（32°～38°N，138°～168°E）代表黑潮延伸體區(qū)域，B 區(qū)（31°～34°N，138°～142°E）代表黑潮大彎曲發(fā)生的主要區(qū)域，再按經(jīng)度范圍（142°～166°E）等分為C 區(qū)（142°～150°E）、D 區(qū)（150°～158°E）和E 區(qū)（158°～166°E）3 個(gè)子區(qū)域。參考馮劭華等[11]提出的面積指數(shù)，本文定義黑潮延伸體各區(qū)域內(nèi)EKE 大于3 000 cm2/s2所有點(diǎn)的加權(quán)和作為EKE 的強(qiáng)度（仍記作EKE）；然后對(duì)不同子區(qū)域EKE大于3 000 cm2/s2的數(shù)據(jù)點(diǎn)進(jìn)行線性回歸擬合，求出一條直線代表該區(qū)域的位置，并計(jì)算作業(yè)點(diǎn)到該直線的垂直距離，記作DIST。

廣義可加模型（Generalized Additive Models,GAM）能夠直接處理響應(yīng)變量與多個(gè)解釋變量之間的非線性關(guān)系，也可以獨(dú)立得出解釋變量對(duì)響應(yīng)變量的影響效應(yīng)[12]。本文計(jì)算每個(gè)漁撈作業(yè)點(diǎn)的EKE，按月平均匹配組成樣本集，利用頻率分析法和GAM 模型，估算適宜的EKE 范圍；然后將CPUE 和黑潮延伸體各區(qū)域的EKE 強(qiáng)度、漁場(chǎng)重心到EKE 位置的距離，按天和月平均數(shù)據(jù)組成樣本集，采用回歸分析法和相關(guān)性分析法，分析黑潮延伸體EKE 的強(qiáng)度變化及其對(duì)漁場(chǎng)空間分布和資源豐度的影響。

3 結(jié)果

3.1 柔魚(yú)CPUE 與作業(yè)漁場(chǎng)EKE 的關(guān)系分析

圖2 示出2010?2016 年7?11 月的西北太平洋柔魚(yú)的月均CPUE 和EKE 的變化，得出7 年月平均CPUE 和EKE 呈負(fù)相關(guān)，相關(guān)系數(shù)為0.78（P ＜ 0.01）。每年的月均CPUE 均呈先增加后遞減趨勢(shì)，7 月為1.16 t/d，而2012 年7 月最低為0.68 t/d；8 月和9 月增高，分別為2.40 t/d 和2.63 t/d，2016 年9 月最高為4.48 t/d；10 月和11 月降低，分別為2 t/d 和1.82 t/d。月均EKE 與CPUE 的季節(jié)變化呈相反趨勢(shì)，7?11 月的平均EKE 分別為264 cm2/s2、206 cm2/s2、100 cm2/s2、145 cm2/s2和173 cm2/s2，其中除了2015 年7 月、8 月和2016 年7 月之外，月平均EKE 低于300 cm2/s2，最低在2011 年9 月為30 cm2/s2，而該月平均CPUE 為3.09 t/d，為月平均CPUE 第二高月份。

圖3 示出2010?2016 年7?11 月CPUE 在各EKE 區(qū)間內(nèi)的累積頻率分布，得出CPUE 主要分布在EKE 范圍為0～100 cm2/s2，該范圍的累積頻率分別為32%、58%、71%、62%和46%，而在0～200 cm2/s2范圍內(nèi)累積頻率分別為53%、82%、89%、78% 和69%。作業(yè)位置主要集中在EKE 較低的區(qū)域。結(jié)合EKE 與CPUE 關(guān)系的GAM 分析圖（圖4），在95%置信度下，CPUE 受到EKE 的影響顯著（P ＜ 0.01），EKE 對(duì)CPUE 的有效影響范圍為0～1 500 cm2/s2。在0～500 cm2/s2范圍內(nèi)，CPUE 隨著EKE 增加而減少，尤其在0～200 cm2/s2呈顯著線性負(fù)相關(guān)趨勢(shì)。在500～1 200 cm2/s2范圍內(nèi)，隨EKE 升高CPUE 增大，呈非線性相關(guān)。GAM 分析結(jié)果顯示，柔魚(yú)漁場(chǎng)的最佳EKE 范圍為25～150 cm2/s2。

圖 1 研究區(qū)域和漁場(chǎng)重心位置（黃色圓圈）Fig. 1 Study area and the locations of fishing ground centroids(yellow circles)

3.2 柔魚(yú)CPUE 與黑潮延伸體EKE 的關(guān)系分析

圖5 示出2010?2016 年7?11 月的黑潮延伸體EKE_A 與CPUE 的相關(guān)分析，可見(jiàn)年平均EKE_A 與CPUE 呈負(fù)相關(guān)，二者的相關(guān)系數(shù)為0.81（P ＜ 0.05），這意味著黑潮延伸體中尺度EKE 越活躍的年份，柔魚(yú)的CPUE 越低。7 年整體月平均EKE_A 與CPUE呈負(fù)相關(guān)，相關(guān)系數(shù)為0.35（P ＜ 0.05）。圖5 還顯示，不同年份的月平均變化趨勢(shì)有所不同。經(jīng)計(jì)算，2010–2011 年，月平均EKE_A 與CPUE 為正相關(guān)，相關(guān)系數(shù)為0.80（P ＜ 0.05）；2012?2016 年，二者為負(fù)相關(guān)，相關(guān)系數(shù)為0.51（P ＜ 0.05）；2015–2016 年負(fù)相關(guān)系數(shù)達(dá)0.78（P ＜ 0.05），其中2016 年9 月的EKE_A 最低為108×104cm2/s2，而該月CPUE 高達(dá)4.47 t/d，為研究期間最高值。對(duì)于2010 年和2011 年，盡管二者月平均為正相關(guān)，但年平均EKE_A 均不高，并且月平均EKE_A的變化幅度較低。這可能與大尺度氣候年際變化有關(guān)，同時(shí)也受到漁場(chǎng)的作業(yè)位置的影響。

圖 2 2010?2016 年7?11 月柔魚(yú)的月平均CPUE 和EKE 變化Fig. 2 Variations of the monthly average CPUE and EKE from July to November during 2010?2016

圖 3 2010?2016 年7?11 月CPUE 在各EKE 區(qū)間內(nèi)的累積頻率分布Fig. 3 Distribution of CPUE frequencies in relation to EKE from July to November during 2010?2016

圖 4 EKE 與CPUE 關(guān)系 的GAM 分析圖Fig. 4 GAM analysis of the relationship between EKE and CPUE

本研究進(jìn)一步計(jì)算黑潮延伸體各子區(qū)域EKE 與CPUE 的相關(guān)關(guān)系，為發(fā)現(xiàn)不同區(qū)域的EKE 變化情況，本研究計(jì)算了黑潮延伸體B、C、D 和E 子區(qū)域2010?2016 年7?11 月逐天EKE 變化，圖6 結(jié)果顯示，黑潮延伸體35 個(gè)月EKE 總和，E 區(qū)為244×106cm2/s2，遠(yuǎn)低于C 區(qū)的652×106cm2/s2和D 區(qū)的648×106cm2/s2，可見(jiàn)黑潮延伸體EKE 強(qiáng)度的年際變化為由西向東逐漸遞減。其中，C 區(qū)和D 區(qū)不同年份的變化規(guī)律略有不同，2013 年、2014 年、2016 年為C 區(qū)EKE 高于D區(qū)，其他年份相反。對(duì)于不同區(qū)域，EKE_B以斜率為0.003 2（R2=0.01）遞增，在2010 年、2013 年和2015 年9 月出現(xiàn)3 個(gè)峰值，分別為24×104cm2/s2、35×104cm2/s2和40×104cm2/s2；EKE_C、EKE_D、EKE_E 呈遞減趨勢(shì)，斜率分別為?0.005 8（R2=0.01）、?0.042 5（R2=0.24）和?0.016 9（R2=0.18）。EKE_C 峰值在2012 年、2013 年7 月和2014 年8 月，分別為125×104cm2/s2、148×104cm2/s2和113×104cm2/s2，EKE_D 在2012 年7 月峰值為150×104cm2/s2；EKE_E 峰值在2012 年、2014 年10 月和2013 年9 月，分別為53×104cm2/s2、62×104cm2/s2和55×104cm2/s2。研究發(fā)現(xiàn)，C、D、E 3 個(gè)區(qū)域出現(xiàn)峰值的月份對(duì)應(yīng)的CPUE 則較低，7?10 月CPUE 最低的年份依次為2012 年、2014 年、2013 年和2012 年，而當(dāng)某區(qū)域出現(xiàn)峰值的年份，該區(qū)域的EKE 總強(qiáng)度也要高于其他年份。以上分析顯示，柔魚(yú)的月平均CPUE變化響應(yīng)黑潮延伸體的不同子區(qū)域的EKE 強(qiáng)度變化，這與該月漁場(chǎng)位置距離動(dòng)力活躍區(qū)域的遠(yuǎn)近有關(guān)，因此有必要進(jìn)一步分析不同區(qū)域EKE 強(qiáng)度與漁場(chǎng)位置的相關(guān)性。

圖7 顯示2010?2016 年7?11 月柔魚(yú)CPUE 漁場(chǎng)分布重心主要范圍為39°～44°N，151°～160°E。經(jīng)計(jì)算，月平均CPUE 與緯度重心呈正相關(guān)，相關(guān)系數(shù)為0.65（P＜0.01）。而月平均CPUE 與經(jīng)度重心無(wú)顯著相關(guān)，故采用CPUE 緯度重心來(lái)反映漁場(chǎng)位置。 由于柔魚(yú)漁場(chǎng)的變動(dòng)，緯度重心分別受不同子區(qū)域的影響，表1 為2010?2016 年柔魚(yú)CPUE 緯度重心與黑潮延伸體EKE 的相關(guān)性分析，相關(guān)系數(shù)取該月相關(guān)最高的區(qū)域。結(jié)果顯示，7 月CPUE 緯度重心對(duì)EKE_D 有明顯響應(yīng)，而2011 年、2014 年和2015 年7 月經(jīng)度分布重心變化幅度大且較長(zhǎng)時(shí)間位于高經(jīng)度區(qū)域，因此，對(duì)EKE_E 響應(yīng)明顯；8 月，除2014 年響應(yīng)EKE_C，其余年份均響應(yīng)EKE_B，且相關(guān)系數(shù)均超過(guò)0.70（P ＜0.01）；9 月，2014?2016 年響應(yīng)EKE_C，2012 年響應(yīng)EKE_E，其他年份與EKE 相關(guān)不顯著，其中，2012 年9 月的經(jīng)度重心為7 年里最高月份；10 月，2010 年響應(yīng)EKE_D，2011?2014 年響應(yīng)EKE_C，而2015 年和2016 年無(wú)顯著相關(guān)；11 月，為漁船返程期，均與EKE無(wú)顯著相關(guān)。由此可見(jiàn)，黑潮延伸體EKE 越活躍，柔魚(yú)的緯度分布重心越偏北。研究發(fā)現(xiàn)，除了與對(duì)應(yīng)的區(qū)域EKE 活躍相關(guān)外，也與漁場(chǎng)作業(yè)位置與活躍區(qū)的距離有關(guān)。

圖 5 2010?2016 年7?11 月CPUE 與黑潮延伸體EKE_A 相關(guān)變化圖Fig. 5 Variations of the monthly average CPUE and EKE_A of the Kuroshio Extension from July to November during 2010?2016

圖 6 2010?2016 年7?11 月黑潮延伸體各子區(qū)域的EKE 變化Fig. 6 Variations of the EKE in sub-areas of the Kuroshio Extension from July to November during 2010?2016

由于7 月和11 月漁場(chǎng)分布經(jīng)度變化幅度大，故取盛漁期8?10 月CPUE 與各區(qū)DIST 進(jìn)行GAM 分析，結(jié)果如圖8 所示。在95%置信度下，GAM 模型顯示A、B、D、E 區(qū)域的DIST 對(duì)CPUE 影響極為顯著（P ＜ 0.01），C 區(qū)影響不顯著（P＞0.05），D 區(qū)變化與A 區(qū)的變化趨勢(shì)最接近。A 區(qū)和D 區(qū)，在DIST 為800～1 000 km 范圍內(nèi)，CPUE 與DIST 呈正相關(guān)，小于800 km和大于1 000 km 置信度較低，二者沒(méi)有顯著相關(guān)性；B 區(qū)由于距離作業(yè)漁區(qū)較遠(yuǎn)，CPUE 與DIST 在600～800 km 和1 300～1 600 km 范圍內(nèi)呈負(fù)相關(guān)，在800～1 000 km 范圍內(nèi)呈正相關(guān)。C 區(qū)DIST 主要分布在700 km 左右范圍內(nèi)，與CPUE 沒(méi)有顯著相關(guān)性。總的來(lái)說(shuō)，漁場(chǎng)與黑潮延伸體的距離過(guò)近（＜ 600 km）或者過(guò)遠(yuǎn)（＞ 1 000 km），則DIST 對(duì)CPUE 沒(méi)有影響，而DIST 在800～1 000 km 范圍內(nèi)，CPUE 隨DIST 增加而增大，其中最適宜的DIST 范圍為850～950 km。

4 討論與分析

柔魚(yú)為短生命周期魚(yú)類(lèi)，對(duì)海洋環(huán)境變化較為敏感，并且具有晝夜垂直洄游的習(xí)性[7]。EKE 活躍的區(qū)域，海面高度變化大，海水常發(fā)生強(qiáng)烈的垂直運(yùn)動(dòng)，并引起溫度梯度急劇變化，超過(guò)柔魚(yú)適應(yīng)范圍，而溫躍層成了魚(yú)類(lèi)垂直運(yùn)動(dòng)時(shí)難以逾越的障礙，并且海流流速大時(shí)不利于魚(yú)群集聚。研究發(fā)現(xiàn)，柔魚(yú)漁場(chǎng)的EKE與CPUE 季節(jié)變化呈顯著負(fù)相關(guān)，相關(guān)系數(shù)為0.78（P ＜0.01），其中9 月柔魚(yú)CPUE 最高，而對(duì)應(yīng)的最低EKE。根據(jù)CPUE 的累積頻率分布，漁場(chǎng)EKE 范圍主要分布在0～200 cm2/s2內(nèi)，其中9 月70%以上EKE 范圍為0～100 cm2/s2，CPUE 在各EKE 區(qū)間范圍內(nèi)累積頻率為逐漸遞減趨勢(shì)。GAM 分析顯示CPUE 受到EKE的影響顯著，有效影響范圍為0～1 500 cm2/s2，最適宜EKE 范圍為25～150 cm2/s2，本研究得出的以上適宜范圍與以往研究結(jié)果基本一致[7]。盡管柔魚(yú)喜好EKE較低的海域，但柔魚(yú)并不集聚在EKE 小于20 cm2/s2的水域，這可能是因?yàn)楹Ａ骱蜏囟忍荻茸兓粼隰~(yú)類(lèi)的適應(yīng)范圍內(nèi)，魚(yú)群也需要由渦旋引起的海水垂直運(yùn)動(dòng)帶來(lái)的餌料。

北太平洋柔魚(yú)通常進(jìn)行季節(jié)性南北洄游，在長(zhǎng)距離的洄游過(guò)程中，黑潮和親潮勢(shì)力的此消彼長(zhǎng)以及黑潮延伸體的環(huán)境條件對(duì)柔魚(yú)的豐度與分布、死亡率以及補(bǔ)充量等產(chǎn)生很大影響[13–14]?？偨Y(jié)前人關(guān)于黑潮延伸體形態(tài)對(duì)柔魚(yú)的影響主要為黑潮延伸體發(fā)生大彎曲，呈大的蛇形擺動(dòng)時(shí)，柔魚(yú)漁場(chǎng)的水溫會(huì)降低，資源豐度降低，漁場(chǎng)分布重心向北偏移[4–6,12]。通常在黑潮大彎曲發(fā)生的年份，黑潮延伸體路徑復(fù)雜多變不穩(wěn)定，常伴有渦旋產(chǎn)生并從黑潮延伸體主軸脫落[3]，從而導(dǎo)致EKE 增加。本研究發(fā)現(xiàn)，黑潮延伸體EKE_A 與CPUE 年平均呈顯著負(fù)相關(guān)關(guān)系，EKE 越活躍的年份，柔魚(yú)的CPUE 越低，與黑潮延伸體大彎曲產(chǎn)生的結(jié)果相一致。但在年際月均變化中，2010 年和2011 年兩年為正相關(guān)，其余年份為負(fù)相關(guān)，其中2015 年和2016 年呈顯著負(fù)相關(guān)，這可能是因?yàn)楹诔毖由祗w大彎曲主要位于B 區(qū)，而由于渦旋和周?chē)Ｋ嗷プ饔卯a(chǎn)生的湍流耗散，黑潮延伸體EKE 表現(xiàn)為由西向東遞減的變化趨勢(shì)，每年遞減的時(shí)間有所區(qū)別，因此，月平均CPUE變化響應(yīng)黑潮延伸體的不同區(qū)域EKE 強(qiáng)度的變化。黑潮延伸體EKE 強(qiáng)度在春、夏季較大，這種季節(jié)變化機(jī)制有多種可能的解釋，如風(fēng)場(chǎng)的季節(jié)性變化引起海洋正壓響應(yīng)等[15]。

表 1 CPUE 緯度重心與黑潮延伸體EKE 的相關(guān)關(guān)系Table 1 Correlations of the latitudinal centroids of CPUE to EKE of the Kuroshio Extension

圖 7 2010?2016 年7?11 月CPUE 漁場(chǎng)分布重心變化Fig. 7 Fishing ground centroids from July to November, 2010?2016

本研究發(fā)現(xiàn)這種變化更多影響著漁場(chǎng)的分布，由于渦旋在運(yùn)動(dòng)過(guò)程中能量不斷耗散[16]，因此其影響具有一定的范圍，這導(dǎo)致響應(yīng)區(qū)域與漁場(chǎng)的距離遠(yuǎn)近相關(guān)，研究結(jié)果顯示，8?10 月盛漁期DIST 在800～1 000 km范圍內(nèi)，CPUE 隨DIST 增加而增大，其中最適宜的DIST 范圍為850～950 km。前文計(jì)算給出的CPUE 2013 年、2014 年和2015 年連續(xù)3 年對(duì)黑潮延伸體的響應(yīng)特征卻各不相同，本文將這3 年盛漁期月份的漁場(chǎng)重心與月均EKE 分布圖疊加，以更好的說(shuō)明響應(yīng)區(qū)域的特征。從圖9 中可以看出，2013 年黑潮延伸體路徑變化復(fù)雜，黑潮延伸體主軸南北側(cè)渦旋較多，這意味著8 月B 區(qū)EKE 增強(qiáng)，9 月B 區(qū)和C 區(qū)EKE 均較強(qiáng)，因此，2013 年8 月和9 月漁場(chǎng)位置分別響應(yīng)B 區(qū)和C 區(qū)的EKE 強(qiáng)度。在2014 年，黑潮延伸體路徑總體平順，渦旋較少，僅C 區(qū)在黑潮延伸體南側(cè)出現(xiàn)一個(gè)明顯的渦旋，該渦旋持續(xù)發(fā)展到10 月，因此，2014 年漁場(chǎng)位置盛漁期期間一直響應(yīng)C 區(qū)的EKE 強(qiáng)度變化。而2015 年黑潮路徑雖然總體彎曲較多，但渦旋數(shù)量不多，在B 區(qū)的一個(gè)環(huán)狀的渦旋持續(xù)到10 月稍減弱，因此2015 年漁場(chǎng)位置盛漁期在8 月份響應(yīng)B 區(qū)的EKE 強(qiáng)度變化。而漁場(chǎng)與黑潮延伸體的距離的響應(yīng)表現(xiàn)為，2014 年的經(jīng)緯度均較高，2015 年的經(jīng)緯度均為7 年中最低的年份，而CPUE 最高月份為10 月，與其他年份的盛漁期時(shí)間相反。結(jié)合表1，研究結(jié)果顯示，CPUE 的緯度重心響應(yīng)該月EKE 強(qiáng)度最高的子區(qū)域，并與其響應(yīng)子區(qū)域的EKE 呈正相關(guān)，因此當(dāng)EKE 強(qiáng)度增大，漁場(chǎng)位置越偏北。

圖 8 不同區(qū)域DIST 與CPUE 關(guān)系的GAM 分析Fig. 8 GAM analysis of the relationship between DIST of sub-areas and CPUE

由此可見(jiàn)，黑潮及黑潮延伸體的EKE 變化不僅影響柔魚(yú)漁場(chǎng)分布，而且影響著漁場(chǎng)與黑潮延伸體的距離，從而影響柔魚(yú)的洄游路線和漁場(chǎng)變動(dòng)，而且大尺度氣候變化也影響著漁場(chǎng)的變動(dòng)[17]，今后需要進(jìn)一步開(kāi)展黑潮及黑潮延伸體與其他環(huán)境因子如溫度、海表面高度和葉綠素等對(duì)柔魚(yú)資源豐度和漁場(chǎng)分布變動(dòng)影響的綜合分析研究。

圖 9 2013?2015 年8?10 月漁場(chǎng)重心與月平均EKE 分布圖（黃色圓圈表示漁場(chǎng)重心）Fig. 9 Maps of fishing ground centroids (yellow circles) with backgrounds of monthly average EKE from August to October, 2013?2015