王雅萌，陳新軍，2，3，4，5，陳 芃，陳子越

（1.上海海洋大學(xué)海洋科學(xué)學(xué)院，上海 201306；2.農(nóng)業(yè)部大洋漁業(yè)開發(fā)重點(diǎn)實(shí)驗(yàn)室，上海 201306；3.國(guó)家遠(yuǎn)洋漁業(yè)工程技術(shù)研究中心，上海 201306；4.大洋漁業(yè)資源可持續(xù)開發(fā)教育部重點(diǎn)實(shí)驗(yàn)室，上海 201306；5.農(nóng)業(yè)部大洋漁業(yè)資源環(huán)境科學(xué)觀測(cè)實(shí)驗(yàn)站，上海 201306）

黑潮是位于北太平洋副熱帶總環(huán)流系統(tǒng)中的西部邊界流，是世界海洋中第二大暖流，黑潮和親潮的廣泛交匯為海洋生物的生長(zhǎng)與發(fā)育帶來(lái)了豐富的餌料，使其附近相關(guān)區(qū)域成為世界海洋中漁業(yè)產(chǎn)量最高的海域之一［1］。據(jù)Seaaroundus網(wǎng)站（http：／／www.seaaroundus.org／）提供的數(shù)據(jù)，1950—2014年黑潮和親潮區(qū)域的漁獲物中，中上層魚種產(chǎn)量分別占到了總產(chǎn)量的67%～82%和22%～45%，主要經(jīng)濟(jì)魚種包括擬沙丁魚（Sardinops sagax）、太平洋褶柔魚（Todarodes pacificus）、白 海 鯡（Clupea pallasii pallasii）、日本竹筴魚（Trachurus japonicus）、日本鳀（Engraulis japonicus）、秋刀魚（Cololabis saira）、鰹（Katsuwonus pelamis）、日 本 鯖（Scomber japonicus）等。水溫變動(dòng)與親潮和黑潮海域的漁獲量變動(dòng)關(guān)系密切，例如沈惠明［2］認(rèn)為，1999年黑潮由于發(fā)生大彎曲，表層水溫偏高使得金槍魚種類出現(xiàn)10多年來(lái)未有的大豐產(chǎn)，而秋刀魚等則減產(chǎn)；ITOH等［3］通過(guò)回歸分析法發(fā)現(xiàn)，1978—2004年沙丁魚產(chǎn)卵數(shù)量與黑潮延伸區(qū)的SST具有線性負(fù)相關(guān)關(guān)系。TIAN 等［4］研究發(fā)現(xiàn)，1970—2000年大型秋刀魚CPUE與黑潮區(qū)域海表溫度呈正相關(guān)，與親潮區(qū)域海表溫度呈負(fù)相關(guān)。余為［5］研究發(fā)現(xiàn)，太平洋年代際震蕩（Pacific decadal oscillation，PDO）與柔魚（Ommastrephe sbartramii）的漁獲量變動(dòng)關(guān)系密切，PDO暖期會(huì)使柔魚漁獲量升高，冷期時(shí)使其下降。以往的研究都是對(duì)單一魚種漁獲量變動(dòng)及其影響因素進(jìn)行分析，尚未有對(duì)黑潮和親潮這2個(gè)大海洋生態(tài)系整體漁獲量變動(dòng)趨勢(shì)及與水溫變動(dòng)和相關(guān)氣候事件的關(guān)系進(jìn)行分析。因此，本文以歷年黑潮和親潮大海洋生態(tài)系中上層魚種漁獲量為研究對(duì)象，結(jié)合海表溫、Nino3.4指數(shù)、PDO等環(huán)境數(shù)據(jù)，利用廣義加性模型（generalized addictive model，GAM）分析漁獲量趨勢(shì)指標(biāo)與溫度變動(dòng)和氣候事件的關(guān)系，以期為未來(lái)黑潮和親潮大海洋生態(tài)系漁業(yè)資源合理開發(fā)和科學(xué)管理提供依據(jù)。

1 材料與方法

1.1 數(shù)據(jù)來(lái)源

漁獲數(shù)據(jù)來(lái)自Seaaroundus網(wǎng)站（http：／／www.seaaroundus.org／），為黑潮和親潮大海洋生態(tài)系［6］（large marine ecosystem，LME，圖1）的所有魚種漁獲量數(shù)據(jù)，漁獲量單位為t，數(shù)據(jù)的時(shí)間范圍為1950—2014年，時(shí)間分辨率為年。

圖1 黑潮和親潮大海洋生態(tài)系分布圖Fig.1 Distribution map of Kuroshio current and Oyashio current

環(huán)境數(shù)據(jù)時(shí)間為1950—2014年，包括黑潮和親潮的表溫（SST）數(shù)據(jù)、PDO指數(shù)以及Nino3.4指數(shù)，分別來(lái)源自英國(guó)大氣數(shù)據(jù)中心（British Atmospheric Data Centre）網(wǎng) 站（https：／／www.metoffice.gov.uk／hadobs／hadisst／data／download.html）、華盛頓大學(xué)大氣與海洋聯(lián)合研究所網(wǎng)站（http：／／research.jisao.washington.edu／pdo／PDO.latest.txt）、NOAA國(guó)家天氣服務(wù)氣候預(yù)報(bào)中心（https：／／origin.cpc.ncep.noaa.gov／products／analysis_monitoring／ensostuff／ONI_v5.php），空間分辨率為0.5°×0.5°，時(shí)間分辨率均為月。

1.2 方法

1.2.1 數(shù)據(jù)的預(yù)處理

提取黑潮和親潮兩個(gè)生態(tài)系統(tǒng)的中上層魚種漁獲量數(shù)據(jù)，大多數(shù)魚類屬于狹溫性魚類［7］，結(jié)合fishbase網(wǎng) 站（https：／／www.fishbase.se／search.php）和fishdb網(wǎng)站（https：／／fishdb.sinica.edu.tw／mobi／home）上的魚種信息，將其劃分暖性、溫型、暖溫性和冷性等種類。分析發(fā)現(xiàn)，1950—2014年這兩個(gè)大海洋生態(tài)系的漁獲物均不包含冷性魚種，根據(jù)兩個(gè)大海洋生態(tài)系的產(chǎn)量大小，主要經(jīng)濟(jì)魚種及其種類屬性見表1。

表1 主要經(jīng)濟(jì)魚種及其種類屬性Tab.1 Main economic species and their attributes

將黑潮與親潮中上層魚種每年的漁獲量（單位：t）數(shù)據(jù)轉(zhuǎn)換為標(biāo)準(zhǔn)化常態(tài)分布值（Z值），利用Z值對(duì)漁獲量變動(dòng)趨勢(shì)進(jìn)行分析，其計(jì)算公式如下［8］：

式中：Ci為黑潮或親潮大海洋生態(tài)系的第i年漁獲量（單位：t），S為黑潮或親潮大海洋生態(tài)系1950—2014年漁獲量的標(biāo)準(zhǔn)差（單位：t）。

水溫的年間變動(dòng)采用兩大海洋生態(tài)系的SST年平均值表示，同時(shí)研究也考慮水溫變化趨勢(shì)，包含了SST時(shí)間上的變化程度和SST空間上的變化程度［9］。黑潮大海洋生態(tài)系SST時(shí)間上的變化程度（temporal variation of Kuroshio current’s SST，TVK）和親潮SST時(shí)間上的變化程度（temporal variation of Oyashio current’s SST，TVO），計(jì)算公式分別為：

式中：Ym為黑潮或親潮某一年第m月所有坐標(biāo)點(diǎn)的SST平均值，Xmi為某一年第m月空間第i個(gè)點(diǎn)的SST值，N為空間總點(diǎn)數(shù)，μT為黑潮或親潮某一年所有月份的SST平均值。

黑潮大海洋生態(tài)系SST空間上的變化程度（spatial variation of Kuroshio current SST，SVK）和親潮SST空間上的變化程度（spatial variation of Oyashio current SST，SVO）的計(jì)算公式分別為：

式中：Yi為黑潮或親潮某一年坐標(biāo)點(diǎn)的12個(gè)月平均SST值，Xmi為某一年第m月空間第i個(gè)點(diǎn)的SST值，μS為黑潮或親潮某一年所有月份的SST平均值，N為空間總點(diǎn)數(shù)。

1.2.2 主成分分析

主成分分析（PCA）的方法能夠?qū)Χ嗑S的變量進(jìn)行降維處理，提取能夠反映變量間相互關(guān)系的共同特征，在漁業(yè)生態(tài)領(lǐng)域有較廣泛的應(yīng)用［10］。分別提取黑潮、親潮中上層漁獲量主成分變化特征值（principal component score，PCs），進(jìn)行主成分年際變化特征分析。

1.2.3 GAM模型

廣義加性模型（generalized addictive model，GAM）是多元線性回歸的擴(kuò)展，其優(yōu)點(diǎn)是能直接處理響應(yīng)變量與多個(gè)解釋變量的非線性關(guān)系［11］。作為一種被廣泛使用的定量分析漁獲率與環(huán)境要素之間關(guān)系的方法，目前在漁業(yè)領(lǐng)域已有廣泛應(yīng)用。擬使用GAM模型分析黑潮和親潮大海洋生態(tài)系中上層魚種漁獲量標(biāo)準(zhǔn)化常態(tài)分布Z值與各變量因子之間的關(guān)系，其中Z值作為響應(yīng)變量，SST、PDO指數(shù)和Nino3.4指數(shù)作為解釋變量，GAM模型表達(dá)式為：

以赤池信息準(zhǔn)則（Akaike information criterion，AIC）來(lái)衡量GAM模型的擬合優(yōu)劣［11］，從而確定最優(yōu)模型。數(shù)據(jù)分析使用軟件為R3.4.2。

2 結(jié)果與分析

2.1 黑潮和親潮中上層魚種Z值變化情況

1）中上層魚種Z值變化?？傮w上，1950—1984年親潮大海洋生態(tài)系Z值大致高于黑潮大海洋生態(tài)系Z值，且親潮大海洋生態(tài)系Z值波動(dòng)更大。1958年親潮大海洋生態(tài)系Z值為峰值2.06，1971年為谷值-1.00，而黑潮大海洋生態(tài)系Z值在-1.50持續(xù)波動(dòng)。1984年以后，黑潮大海洋生態(tài)系Z值持續(xù)升高并在1987年以后一直高于親潮大海洋生態(tài)系Z值，1988年黑潮大海洋生態(tài)系Z值達(dá)到峰值1.80。1995年黑潮和親潮大海洋生態(tài)系Z值一起下降到谷值 -0.08、-2.20。2003和2005年黑潮和親潮大海洋生態(tài)系Z值相差較小，2010年親潮大海洋生態(tài)系Z值再次降低到谷值-1.90（圖2）。

圖2 1950—2014年黑潮和親潮中上層魚種Z值變化趨勢(shì)Fig.2 Variation of Z values of mid-pelagic fish species in Kuroshio current and Oyashio current from 1950 to 2014

2）暖性魚種Z值變化。1950—1972年黑潮和親潮大海洋生態(tài)系暖性魚種Z值變化均比較平緩（圖3），1975年前親潮Z值有較大增加，1973年達(dá)到峰值2.30。1975—1980年黑潮大海洋生態(tài)系暖性魚種Z值高于親潮大海洋生態(tài)系Z值。1984年親潮大海洋生態(tài)系暖性魚種Z值達(dá)到最大峰值2.90。1985—1999年親潮大海洋生態(tài)系暖性魚種Z值普遍低于黑潮。2001年親潮大海洋生態(tài)系暖性魚種Z值再次達(dá)到峰值2.60。2005年以后，黑潮和親潮大海洋生態(tài)系暖性魚種Z值變化均趨于平緩，都在-1.00（圖3）。

圖3 1950—2014年黑潮和親潮暖性魚種Z值變化趨勢(shì)Fig.3 Variation of Z values of warm fish species in Kuroshio current and Oyashio current from 1950 to 2014

3）溫性魚種Z值變化。1950—1975年黑潮和親潮大海洋生態(tài)系溫性魚種Z值變化趨勢(shì)基本相同（圖4）。1976—1980年親潮大海洋生態(tài)系Z值高于黑潮大海洋生態(tài)系且相差較大。1981—1996年黑潮和親潮大海洋生態(tài)系Z值變化趨勢(shì)也大致相同，但1996年黑潮大海洋生態(tài)系Z值達(dá)到最大峰值3.50，而親潮大海洋生態(tài)系Z值只有1.60。1997—2002年黑潮大海洋生態(tài)系Z值高于親潮，且1998年為親潮大海洋生態(tài)系Z值達(dá)到谷值-0.8。2004—2009年親潮大海洋生態(tài)系Z值高于黑潮，且2007年為親潮大海洋生態(tài)系Z值達(dá)到最大峰值2.40（圖4）。

圖4 1950—2014年黑潮和親潮溫性魚種Z值變化趨勢(shì)Fig.4 Variation of Z values of temperate fish species in Kuroshio current and Oyashio current from 1950 to 2014

4）暖溫性魚種Z值變化。1950—1989年黑潮和親潮大海洋生態(tài)系暖溫性魚種Z值變化趨勢(shì)均是增加，1984年親潮大海洋生態(tài)系Z值達(dá)到峰值2.47，黑潮Z值在1989年達(dá)到峰值2.05。1989年以后，黑潮和親潮大海洋生態(tài)系暖溫性魚種Z值均呈下降趨勢(shì)，且1998—2014年親潮Z值一直高于黑潮Z值（圖5）。

圖5 1950—2014年黑潮和親潮暖溫性魚種Z值變化趨勢(shì)Fig.5 Variation of Z values of warm-temperate fish species in Kuroshio current and Oyashio current from 1950 to 2014

2.2 黑潮和親潮大海洋生態(tài)系SST變化

1950—2014年親潮大海洋生態(tài)系SST在5.9～7.3℃之間波動(dòng)，1985—1990年SST呈持續(xù)的增長(zhǎng)趨勢(shì)并于1990年達(dá)到峰值，為7.2℃；2004—2014年親潮大海洋生態(tài)系SST在較高值7℃波動(dòng)（圖6）。1950—1986年黑潮大海洋生態(tài)系SST在21.5～22.5℃波動(dòng)且趨勢(shì)較穩(wěn)定，1986年下降至谷值21.7℃，1998年達(dá)到最大峰值23.1℃，之后不斷波動(dòng)下降（圖6）。SST空間上的變動(dòng)，總體上親潮大海洋生態(tài)系小于黑潮大海洋生態(tài)系，波動(dòng)范圍分別為1.47～2.34℃和3.5～4.4℃（圖7）。SST時(shí)間上的變動(dòng)，總體上親潮大海洋生態(tài)系大于黑潮大海洋生態(tài)系，波動(dòng)范圍分別為3.9～4.9℃和3.2～4.0℃（圖8）。

圖6 1950—2014年黑潮和親潮SST變化Fig.6 Variation of SST in Kuroshio current and Oyashio current from 1950 to 2014

圖7 1950—2014年黑潮和親潮SST空間變化趨勢(shì)Fig.7 Spatial variation of SST in Kuroshio current and Oyashio current from 1950 to 2014

圖8 1950—2014年黑潮和親潮SST時(shí)間變化趨勢(shì)Fig.8 Temporal variation of SST in Kuroshio current and Oyashio current from 1950 to 2014

2.3 主成分年際變化特征

主成分分析結(jié)果可以反映漁獲物組成變化的主要特征，對(duì)黑潮的漁獲物產(chǎn)量進(jìn)行主成分分析后得出，前兩位的主成分方差累計(jì)解釋率達(dá)到60.0%，黑潮的第一、第二主成分的解釋率分別為29.2%、26.8%，并且在1987年左右均出現(xiàn)得分轉(zhuǎn)換。對(duì)親潮的漁獲物產(chǎn)量進(jìn)行主成分分析后得出，前兩位的主成分方差累計(jì)解釋率達(dá)到62.3%，親潮的第一、第二主成分的解釋率分別為35.5%、26.8%，并且第一主成分在1987年出現(xiàn)峰值，第二主成分在1987年出現(xiàn)得分轉(zhuǎn)換。

2.4 環(huán)境因子和氣候因子相關(guān)性分析

分析發(fā)現(xiàn)（表2），黑潮海域Nino3.4與TVK、PDO在統(tǒng)計(jì)上分別具有負(fù)相關(guān)和正相關(guān)（P＜0.05），PDO與SST空間變化程度在統(tǒng)計(jì)上具有正相關(guān)（P＜0.05）。分析也發(fā)現(xiàn)（表3），親潮海域Nino3.4與TVO、PDO在統(tǒng)計(jì)上分別具有負(fù)相關(guān)和正相關(guān)（P＜0.05），PDO與SST、SST與TVO在統(tǒng)計(jì)上均具有正相關(guān)（P＜0.05）。

表2 黑潮海域環(huán)境因子和氣候因子相關(guān)性系數(shù)Tab.2 Correlation coefficient between Kuroshio environmental factors and climate factors

表3 親潮海域環(huán)境因子和氣候因子相關(guān)性系數(shù)Tab.3 Correlation coefficient between Oyashio environmental factors and climate factors

2.5 GAM 模型結(jié)果

依據(jù)AIC最小值確定了以下最優(yōu)擬合模型（表4）：黑潮中上層魚種Z值的最優(yōu)擬合變量組合為SST＋PDO ＋SVK，模型偏差解釋率為53.0%，其中，SVK貢獻(xiàn)率最大，為28.0%，其次為SST和PDO，分別為14.9%和10.4%。黑潮中上層暖性魚種Z值的最優(yōu)擬合變量組合為PDO＋Nino3.4，模型偏差解釋率為31.0%，其中，PDO為17.6%，多于Nino3.4指數(shù)的13.7%。黑潮中上層溫性魚種Z值的最優(yōu)擬合變量為SST，其模型偏差解釋率為25.0%。黑潮中上層暖溫性魚種Z值的最優(yōu)擬合變量組合為PDO＋SST，模型偏差解釋率為40.0%，其中，PDO為24.5%，多于SST的15.4%。親潮中上層魚種Z值的最優(yōu)擬合變量組合為SST＋SVO，其中，SST為44.5%，多于SVO的16.9%。親潮中上層暖性魚種Z值的最優(yōu)擬合變量組合為PDO＋Nino3.4＋SVO，其中，PDO貢獻(xiàn)率最大，為19.59%，其次為SVO和Nino3.4，貢獻(xiàn)率分別為15.7%和11.01%。親潮中上層溫性魚種Z值的最優(yōu)擬合變量為SST，其模型偏差解釋率為27.1%。親潮中上層暖溫性魚種Z值的最優(yōu)擬合變量組合為PDO＋SST＋Nino3.4＋SVO，模型偏差解釋率為60.8%，其中，PDO貢獻(xiàn)率最大，為24.3%，其次為SST和Nino3.4，分別為23.4%和9.1%。

表4 各Z值最優(yōu)擬合模型參數(shù)Tab.4 The optimal model parameters of each Z value

2.5.1 黑潮和親潮大海洋生態(tài)系中上層漁獲量Z值

黑潮大海洋生態(tài)系SST和PDO指數(shù)對(duì)其中上層魚種漁獲量Z值的影響極顯著，隨SST、SVK和PDO指數(shù)的增加，其中上層魚種漁獲量Z值也顯著增加（圖9）。親潮大海洋生態(tài)系SST對(duì)其Z值的影響極其顯著，SST與親潮大海洋生態(tài)系中上層魚種漁獲量Z值之間有正相關(guān)關(guān)系，當(dāng)SST為6.1℃和6.8℃時(shí)，Z值為峰值，當(dāng)SST為6.6℃和7.2℃時(shí)，Z值為谷值（圖10）。親潮大海洋生態(tài)系SST在空間上的變化程度對(duì)其Z值的影響也顯著，兩者間具有顯著的正相關(guān)關(guān)系（圖10）。

圖9 黑潮海域中上層魚種Z值的GAM 模型Fig.9 GAM model of Z value of mid-pelagic fish species in Kuroshio current

圖10 親潮海域中上層魚種Z值的GAM 模型Fig.10 GAM model of Z value of mid-pelagic fish species in Oyashio current

2.5.2 黑潮和親潮大海洋生態(tài)系暖性魚種

PDO指數(shù)對(duì)黑潮大海洋生態(tài)系中上層暖性魚種Z值的影響極其顯著，兩者具有正相關(guān)關(guān)系，Nino 3.4指數(shù)與黑潮大海洋生態(tài)系中上層暖性魚種Z值間則具有顯著的負(fù)相關(guān)關(guān)系（圖11）。

圖11 黑潮海域中上層暖性魚種Z值的GAM 模型Fig.11 GAM model of Z value of mid-pelagic warm fish species in Kuroshio current

PDO和Nnio3.4指數(shù)對(duì)親潮大海洋生態(tài)系中上層暖性魚種Z值的影響極其顯著，SST與親潮大海洋生態(tài)系Z值間有負(fù)相關(guān)關(guān)系，PDO與親潮大海洋生態(tài)系Z值間有正相關(guān)關(guān)系。親潮大海洋生態(tài)系SST在空間上的變化程度也會(huì)顯著影響其Z值，SST在空間上的變化程度與親潮大海洋生態(tài)系中上層暖性魚種Z值間具有負(fù)相關(guān)關(guān)系（圖12）。

圖12 親潮海域中上層暖性魚種Z值的GAM 模型Fig.12 GAM model of Z value of mid-pelagic warm fish species in Oyashio current

2.5.3 黑潮和親潮大海洋生態(tài)系溫性魚種

黑潮海域SST對(duì)其中上層溫性魚種Z值的影響顯著，兩者具有正相關(guān)關(guān)系（圖13）。親潮大海洋生態(tài)系SST和PDO指數(shù)與其溫性魚種Z值間具有顯著的正相關(guān)關(guān)系（圖14）。

圖13 黑潮海域中上層溫性魚種Z值的GAM 模型Fig.13 GAM model of Z value of mid-pelagic temperate fish species in Kuroshio current

圖14 親潮海域中上層溫性魚種Z值的GAM 模型Fig.14 GAM model of Z value of mid-pelagic temperate fish species in Oyashio current

2.5.4 黑潮和親潮大海洋生態(tài)系暖溫性魚種

PDO指數(shù)對(duì)黑潮海域中上層暖溫性魚種Z值的影響極其顯著，兩者間具有正相關(guān)關(guān)系；黑潮大海洋生態(tài)系SST與其Z值間則具有顯著的負(fù)相關(guān)關(guān)系（圖15）。

圖15 黑潮海域中上層暖溫性魚種Z值的GAM 模型Fig.15 GAM model of Z value of mid-pelagic warm-temperature fish species in Kuroshio current

親潮SST、PDO指數(shù)、Nino3.4指數(shù)以及SST在空間上的變化均會(huì)顯著影響親潮中上層暖溫性魚種Z值，SST的增加對(duì)親潮大海洋生態(tài)系Z值的影響波動(dòng)較大，SST為6.8℃時(shí)，Z值為峰值，在SST為6.55℃和7.05℃時(shí)，Z值為谷值；PDO指數(shù)與親潮大海洋生態(tài)系暖溫性魚種Z值間存在正相關(guān)關(guān)系（圖16）。

圖16 親潮海域中上層暖溫性魚種Z值的GAM 模型Fig.16 GAM model of Z value of pelagic warm-temperature fish species in Oyashio current

3 討論

本研究通過(guò)對(duì)黑潮和親潮大海洋生態(tài)系中上層暖性、溫性以及暖溫性魚種標(biāo)準(zhǔn)常態(tài)分布值Z的分析，發(fā)現(xiàn)兩大海洋生態(tài)系的中上層魚種總漁獲量Z值變動(dòng)在時(shí)間上的變化有一定差異（圖2），但是各種類型的漁獲量Z值變化有著較為一致的趨勢(shì)，造成這種差異的原因可能主要在于黑潮和親潮大海洋生態(tài)系獲得的各種類型漁獲物的比例存在不同。根據(jù)歷年漁獲量數(shù)據(jù)的分析，黑潮的主要漁獲物以暖溫種類為主，親潮則以溫性種類為主。

研究利用GAM模型對(duì)影響Z值變動(dòng)的環(huán)境因素進(jìn)行了分析。研究發(fā)現(xiàn)，總體上親潮大海洋生態(tài)系的GAM模型解釋率要高于黑潮大海洋生態(tài)系的，影響其漁獲量Z值變化的顯著指標(biāo)也多于黑潮，原因可能為親潮區(qū)域流經(jīng)冷水性海流，流速較慢，海洋環(huán)境特征穩(wěn)定，環(huán)境特征的較小變化及氣候變化就會(huì)更為直接地對(duì)大海洋生態(tài)系中的種類分布與資源狀況造成影響；而黑潮區(qū)域海流為暖水性的，流速較快，且包含了較多的與黑潮交匯的區(qū)域，產(chǎn)生渦流、鋒面等不同的海洋動(dòng)力學(xué)要素，是多種魚類的索餌場(chǎng)所在。漁場(chǎng)的環(huán)境指標(biāo)不僅僅是溫度，區(qū)域內(nèi)部的海面高度、渦動(dòng)能（渦流），流速（黑潮），水色如葉綠素濃度（代表了交匯區(qū)索餌場(chǎng)的鋒面及餌料情況）等環(huán)境指標(biāo)可能更為重要［13-15］。此外，本研究采用應(yīng)變量式漁獲量數(shù)據(jù)，捕撈努力量在這兩個(gè)區(qū)域的年間差異可能會(huì)對(duì)模型的結(jié)果造成影響，因此后續(xù)的分析需要收集相關(guān)捕撈數(shù)據(jù)對(duì)模型進(jìn)行補(bǔ)充。

研究認(rèn)為，影響黑潮大海洋生態(tài)系總體漁獲量Z值的顯著性指標(biāo)是PDO，親潮大海洋生態(tài)系是SST和SVO，這可能與它們的漁獲物組成有關(guān)，中上層魚類種群不同的棲息地分布、產(chǎn)卵繁殖以及營(yíng)養(yǎng)級(jí)特點(diǎn)決定了對(duì)環(huán)境變化的響應(yīng)存在差異［10］。黑潮大海洋生態(tài)系的漁獲物主要是暖溫性魚種，親潮大海洋生態(tài)系主要是溫性魚種。研究表明，PDO能夠顯著影響到黑潮區(qū)域的水文特征，如溫度鋒、渦動(dòng)能、海表面高度等，由此帶來(lái)了漁場(chǎng)環(huán)境的改變［16-18］，例如，PDO負(fù)位相時(shí)期的西北太平洋SST降低，冬季日本鳀漁場(chǎng)重心隨等溫線向南移動(dòng)［19］，北太平洋大海洋生態(tài)系沙丁魚的資源量受到SST的影響，其CPUE與PDOI呈正相關(guān)［20］。相反，親潮來(lái)源于白令海峽和鄂霍茨克海，為從北向南的冷水性海流，包含了SST和SVO的影響，這可以用前述親潮水體海洋環(huán)境相對(duì)穩(wěn)定來(lái)解釋，水溫的變暖能夠較好的反映到漁獲物的變化上。當(dāng)氣候發(fā)生變化時(shí)，由于海洋自身的熱慣性，特定區(qū)域的海洋生態(tài)系統(tǒng)響應(yīng)規(guī)模較氣候系統(tǒng)小且表現(xiàn)出時(shí)空非同步性［21］，魚類資源豐度可能存在滯后的響應(yīng)，此外，李婷等［22］發(fā)現(xiàn)，凈初級(jí)生產(chǎn)力變化對(duì)漁場(chǎng)資源豐度具有顯著影響，其中相關(guān)機(jī)制有待后續(xù)進(jìn)一步研究。

本文初步分析了SST和PDO、Nino3.4指數(shù)對(duì)黑潮和親潮大海洋生態(tài)系中上層暖性、溫性以及暖溫性魚種的影響，可為今后環(huán)境和氣候事件對(duì)黑潮和親潮大海洋生態(tài)系經(jīng)濟(jì)魚種資源量及其影響機(jī)理的探究提供一定依據(jù)。