溫利紅,張 衡,方 舟,3,4,5,6*,陳新軍,3,4,5,6

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

引 言

鳶烏賊(Sthenoteuthisoualaniensis)，隸屬頭足綱(Cephalopoda)，柔魚(yú)科(Ommastrephidae)，鳶烏賊屬(Sthenoteuthis)。該種類(lèi)生命周期短，繁殖能力強(qiáng)，資源量豐富，是世界重要的經(jīng)濟(jì)性頭足類(lèi)，廣泛分布于印度洋、太平洋的赤道和亞熱帶海域[1-2]。目前印度洋西北部海域和南海是我國(guó)捕撈鳶烏賊的主要漁場(chǎng)[3-4]。印度洋沿岸國(guó)漁業(yè)資源豐富，但除少數(shù)發(fā)達(dá)國(guó)家外，大多數(shù)國(guó)家的經(jīng)濟(jì)發(fā)展水平不高，漁業(yè)資源也尚未得到開(kāi)發(fā)，相對(duì)于其他公海鳶烏賊漁場(chǎng)，開(kāi)發(fā)潛力較大，日益受到人們重視[4-5]。由于鳶烏賊具有晝夜垂直遷移習(xí)性和趨光性，目前我國(guó)針對(duì)該物種的主要捕撈方式為燈光敷網(wǎng)、燈光罩網(wǎng)和魷釣[6-8]。

合理分析環(huán)境因素與漁獲量的關(guān)系，是了解和掌握鳶烏賊資源量及分布規(guī)律，合理開(kāi)發(fā)和利用鳶烏賊漁業(yè)資源的重要前提之一，但兩者之間的關(guān)系往往不是簡(jiǎn)單的線性關(guān)系。單位捕撈努力量漁獲量(catch per unit effort,CPUE)通常被假定與漁業(yè)資源量之間成正比關(guān)系，是衡量漁業(yè)資源密度的一個(gè)重要指標(biāo)，被廣泛應(yīng)用于漁業(yè)資源評(píng)估與管理中[9-10]。由于在實(shí)際漁業(yè)生產(chǎn)過(guò)程中，受漁場(chǎng)時(shí)空分布、海洋環(huán)境、氣候變化及作業(yè)漁船的捕撈方式等多種因素的影響，CPUE與資源量之間存在較為復(fù)雜的關(guān)系，不能準(zhǔn)確表達(dá)出不同作業(yè)方式之間的捕撈努力量差異，無(wú)法全面、真實(shí)地反映漁業(yè)資源量及分布情況[11]。在此情況下，需要對(duì)CPUE進(jìn)行標(biāo)準(zhǔn)化處理，去除影響因素，減小誤差，從而能夠科學(xué)的衡量漁業(yè)資源量[10]。目前，研究漁業(yè)資源與環(huán)境之間的關(guān)系模型方法有：廣義線性模型(generalized linear model,GLM)和廣義可加性模型(generalized additive model,GAM)，這兩種模型方法已被國(guó)內(nèi)外廣泛到漁業(yè)中[12]。其中GLM模型主要解決響應(yīng)變量和解釋變量之間的線性問(wèn)題，由于CPUE受多個(gè)變量的影響，且與各影響變量之間的關(guān)系較為復(fù)雜，可能存在非線性關(guān)系，而GAM模型是GLM模型的非參數(shù)化拓展，是一種非線性關(guān)系模型，能夠較好地處理非線性問(wèn)題，更加有效地探討漁業(yè)資源與影響變量之間的關(guān)系，提高研究的準(zhǔn)確性[13-14]。為此，本文根據(jù)2017—2019年中國(guó)遠(yuǎn)洋漁業(yè)協(xié)會(huì)魷釣技術(shù)組和公海圍拖網(wǎng)技術(shù)組提供的生產(chǎn)統(tǒng)計(jì)數(shù)據(jù)，利用GLM模型篩選關(guān)鍵時(shí)空環(huán)境因子，以GAM模型分析印度洋北部鳶烏賊漁場(chǎng)與時(shí)空環(huán)境因子間的關(guān)系，并以此對(duì)進(jìn)行CPUE標(biāo)準(zhǔn)化，為我國(guó)對(duì)印度洋鳶烏賊漁業(yè)資源的合理利用和漁情預(yù)報(bào)提供相關(guān)的參考依據(jù)。

1 材料與方法

1.1 材料來(lái)源

印度洋北部鳶烏賊生產(chǎn)統(tǒng)計(jì)資料來(lái)自中國(guó)遠(yuǎn)洋漁業(yè)協(xié)會(huì)魷釣技術(shù)組和公海圍拖網(wǎng)技術(shù)組，該統(tǒng)計(jì)資料包含燈光敷網(wǎng)、燈光罩網(wǎng)和魷釣3種作業(yè)方式。統(tǒng)計(jì)內(nèi)容包括作業(yè)日期、作業(yè)次數(shù)、作業(yè)經(jīng)度、作業(yè)緯度和漁獲量?？臻g分辨率為1°×1°，本文研究的海域范圍為10°～25°N，50°～75°E(圖1)。時(shí)間為2017—2019年。

圖1 印度洋北部海域范圍(注：黑色線框內(nèi)為研究區(qū)域)

海表面溫度(Sea Surface Temperature,SST)、海表面高度(Sea Surface Height,SSH)、海表鹽度(Sea Surface Salinity,SSS)、光合有效輻射(Photosynthetically Active Radiation,PAR)、風(fēng)速(Wind Speed,WS)以及流速(Current Speed,U)被認(rèn)為是影響大洋性魷魚(yú)漁場(chǎng)分布的重要環(huán)境因子[15-16]，因此在本研究中納入CPUE標(biāo)準(zhǔn)化的影響因子。海表溫度和光合有效輻射數(shù)據(jù)來(lái)源于美國(guó)大氣與海洋局(National Ocean and Atmosphere Administration,NOAA)中太平洋觀測(cè)網(wǎng)點(diǎn)(https://oceanwatch.pifsc.noaa.gov/erddap/index.html)；海表鹽度、海表面高度、風(fēng)速和流速數(shù)據(jù)來(lái)源于夏威夷大學(xué)網(wǎng)站(http://apdrc.soest.hawaili.edu/data/data.php)。時(shí)間分辨率為月，空間分辨率為0.5°×0.5°。

1.2 數(shù)據(jù)預(yù)處理

1.2.1 計(jì)算

單位捕撈努力量漁獲量(catch per unit effort, CPUE)可以作為表征鳶烏賊資源密度的指標(biāo)，計(jì)算公式如下[17]：

(1)

式中，CPUE單位為t/次；C表示一艘漁船一天的產(chǎn)量；E表示其對(duì)應(yīng)的作業(yè)次數(shù)，燈光敷網(wǎng)和燈光罩網(wǎng)以網(wǎng)次計(jì)算、魷釣是按照每天作業(yè)位置的變化次數(shù)來(lái)計(jì)算。

1.2.2 數(shù)據(jù)匹配

通過(guò)Excel和Matlab軟件的相關(guān)程序包，將下載的海表溫度(SST)、海表鹽度(SSS)、海表面高度(SSH)、光合有效輻射(PAR)、風(fēng)速(WS)和流速(U)等環(huán)境因子，利用克里金插值法與漁業(yè)數(shù)據(jù)(經(jīng)緯度、產(chǎn)量、標(biāo)準(zhǔn)化CPUE以及捕撈方式)進(jìn)行匹配[6-7]，使得環(huán)境數(shù)據(jù)與漁業(yè)數(shù)據(jù)一一對(duì)應(yīng)，均調(diào)整為0.5°×0.5°的分辨率。

1.3 變量的相互獨(dú)立性檢驗(yàn)

本研究根據(jù)鳶烏賊漁業(yè)特點(diǎn)選擇解釋變量為年、月、經(jīng)度、緯度、海表面溫度、海表鹽度、海表面高度、光合有效輻射、風(fēng)速、流速和作業(yè)方式(fishing types, FT)。利用方差膨脹系數(shù)(VIF)和Spearman相關(guān)系數(shù)對(duì)解釋變量進(jìn)行相互獨(dú)立性檢驗(yàn)。當(dāng)VIF<10時(shí)，表明變量之間不存在多重共線性，VIF>10時(shí)，表明變量之間存在嚴(yán)重的多重共線性[18]。根據(jù)獨(dú)立性檢驗(yàn)，結(jié)果顯示，各個(gè)解釋變量的VIF均小于10，表明解釋變量之間不存在嚴(yán)重多重共線性問(wèn)題(表1)。

表1 解釋變量間方差膨脹因子

1.4 研究方法

目前CPUE標(biāo)準(zhǔn)化方法主要采用廣義線性模型(GLM)和廣義加性模型(GAM)。國(guó)內(nèi)外較多研究表明，利用GLM模型可獲得影響CPUE的主要因子及其貢獻(xiàn)度,并可考慮多因子的交互效應(yīng)；而GAM模型在CPUE標(biāo)準(zhǔn)化中的效果優(yōu)于GLM模型[12,19-20]。因此本研究首先利用廣義線性模型(GLM)對(duì)影響鳶烏賊資源密度的因子進(jìn)行顯著性檢驗(yàn)，篩選出具有顯著性影響的因子，隨后利用廣義加性模型(GAM)對(duì)顯著性因子和鳶烏賊資源密度關(guān)系進(jìn)行研究，并利用該模型嘗試對(duì)不同時(shí)間尺度(年和月)鳶烏賊CPUE進(jìn)行標(biāo)準(zhǔn)化。

GLM模型假設(shè)響應(yīng)變量的期望值與解釋變量呈線性關(guān)系[21]：

(2)

式中，g為鏈接函數(shù);μi=E(Yi),Yi為第i個(gè)響應(yīng)變量;Xi為第i個(gè)響應(yīng)變量的解釋變量;β為模型估計(jì)參數(shù)。本研究假設(shè)CPUE服從對(duì)數(shù)正態(tài)分布，則GLM模型表達(dá)式為：

ln(CPUE+1)=Year+Month+Latitude+Longitude+SST+SSS+PAR+WS+U+FT+ε

(3)

式中，CPUE為每艘船每次作業(yè)的捕撈產(chǎn)量，對(duì)CPUE加上常數(shù)1是為解決CPUE出現(xiàn)零值的現(xiàn)象；ε為u誤差項(xiàng)，假設(shè)其服從正態(tài)分布。GLM模型中，將時(shí)間(年、月)、空間(經(jīng)度、緯度)、環(huán)境(SST、SSS、PAR、WS、SSH、U)和作業(yè)方式因子作為解釋變量，其中變量年、月、經(jīng)度、緯度和作業(yè)方式作為離散變量，其他變量(SST、SSS、PAR、WS、SSH、U)作為連續(xù)變量。

GAM是GLM模型的非線性拓展，通過(guò)將函數(shù)與相應(yīng)變量進(jìn)行移動(dòng)的變化，將基于指數(shù)分布的回歸與一般線性回歸進(jìn)行整合[22]，即：

g(μi)=α+∑i=1fi(xi)+ε

(4)

式中，g為鏈接函數(shù);μi=E(Yi)，Yi為第i個(gè)響應(yīng)變量;Xi為第i個(gè)響應(yīng)變量的解釋變量;ε為模型估計(jì)參數(shù);fi為平滑函數(shù)。根據(jù)GLM模型篩選出的顯著性解釋變量依次加入GAM模型，則GAM模型的表達(dá)式為：

ln(CPUE+1)～s(Year)+s(Month)+s(Latitude)+s(Longitude)+
s(SST)+s(SSS)+s(PAR)+s(WS)+s(U)+factor(FT)+ε

(5)

式中，對(duì)CPUE+1是為防止響應(yīng)變量出現(xiàn)零值，再進(jìn)行了對(duì)數(shù)化處理；s為自然立方樣條平滑(natural cube spline smoother)；s(year)為年效應(yīng)；s(month)為月效應(yīng)；s(latitude)為緯度效應(yīng)；s(longitude)為經(jīng)度效應(yīng)；s(SST)為海表溫度效應(yīng)；s(SSS)為海表鹽度效應(yīng)；s(PAR)為光合有效輻射效應(yīng)；s(WS)為風(fēng)速效應(yīng)；s(U)為流速效應(yīng)；由于作業(yè)方式的數(shù)值組成太少，因此將作業(yè)方式以因子factor(FT)形式處理。

根據(jù)赤池信息量準(zhǔn)則(AIC)值，選取最佳模型[23]。AIC值計(jì)算如下：

AIC=-2lnl(p1,p2,…,pm,σ2)+2m

(6)

式中，m為模型中參數(shù)的個(gè)數(shù)。

本研究使用R(V3.2.2)進(jìn)行統(tǒng)計(jì)分析與處理。

2 結(jié)果與分析

2.1 解釋變量ln(CPUE+1)的統(tǒng)計(jì)分布檢驗(yàn)

首先驗(yàn)證ln(CPUE+1)是否服從正態(tài)分布，經(jīng)K-S檢驗(yàn)，ln(CPUE+1)的數(shù)據(jù)點(diǎn)在正態(tài)Q-Q圖中基本形成一條直線(圖2-B)，ln(CPUE+1)服從正態(tài)分布(圖2-A)，這說(shuō)明本文研究中關(guān)于ln(CPUE+1)服從正態(tài)分布的假設(shè)是合理的，可以運(yùn)用GLM模型和GAM模型進(jìn)行分析。

圖2 2017—2019年印度洋北部鳶烏賊ln(CPUE+1)的頻次分布及其檢驗(yàn)

2.2 模型分析

2.2.1 GLM模型分析

通過(guò)GLM模型對(duì)各因子進(jìn)行顯著性檢驗(yàn)(見(jiàn)表2)。本文以P<0.05來(lái)確認(rèn)變量是否具有顯著性，結(jié)果表明，年、月、緯度、SST、WS、U以及作業(yè)方式均為顯著性變量，且除因子SST對(duì)CPUE的影響為顯著性外，其他因子都為極顯著性(P<0.01)；經(jīng)度、SSS、PAR及SSH為不顯著性變量(P>0.05)，對(duì)CPUE的影響不明顯。因此，選擇7個(gè)顯著性解釋變量(年、月、緯度、SST、WS、U以及作業(yè)方式)納入GAM模型對(duì)CPUE進(jìn)行標(biāo)準(zhǔn)化。

表2 GLM模型自變量顯著性檢驗(yàn)

2.2.2 GAM模型分析

隨后將上述具有顯著性的解釋變量年、月、緯度、SST、WS、U以及作業(yè)方式逐一加入GAM模型中，進(jìn)行運(yùn)算分析，發(fā)現(xiàn)所選取的因子建立的GAM模型，AIC值最小，其擬合效果最好(見(jiàn)表3)，最終得到最佳的GAM模型為：

ln(CPUE+1)～s(Year)+s(Month)+s(Latitude)+s(SST)+s(WS)+s(U)+factor(FT)

(7)

并且根據(jù)P值，得出年、月、緯度、SST、WS、U及作業(yè)方式均為顯著性變量，對(duì)CPUE的影響都為極顯著性(P<0.01)。模型對(duì)CPUE的總偏差解釋為21.6%，其中作業(yè)方式變量對(duì)CPUE的影響最大，解釋了5.6%的總偏差，說(shuō)明作業(yè)方式對(duì)CPUE的影響最大；隨后影響由大到小依次是年(5.22%)、緯度(4.94%)、月(3.94%)、SST(1.2%)、U(0.6%)、WS(0.1%)。

表3 GAM模型統(tǒng)計(jì)結(jié)果

2.3 各因子對(duì)CPUE的影響

從時(shí)間因子來(lái)看，年為極顯著性變量，對(duì)CPUE具有較大影響；CPUE總體在2017—2018年呈增長(zhǎng)趨勢(shì)，至2018年達(dá)到峰值，隨后呈小幅度下降趨勢(shì)，2019年較為集中。月變化對(duì)CPUE的影響也極顯著，總體變化幅度較大，具有明顯的季節(jié)變化，其中1—7月呈下降趨勢(shì)，并在7月達(dá)到最低值，此后隨月份逐漸上升，并在12月份達(dá)到最大值。從空間因子來(lái)看，緯度為極顯著性變量，變化趨勢(shì)總體呈波浪式起伏，且變化幅度較大，首先在10°N～12°N處呈下降趨勢(shì)，并在12°N附近達(dá)最小值，后隨著緯度的增加呈上升趨勢(shì)，于18°N附近達(dá)到最大值，隨后又呈下降趨勢(shì)，其中CPUE在14°N～17°N范圍分布較為集中(圖3)。

圖3 基于GAM模型的時(shí)空與環(huán)境因子效應(yīng)對(duì)印度洋北部鳶烏賊CPUE的影響

從環(huán)境因子來(lái)看，SST為極顯著性變量，由圖2可知，印度洋北部鳶烏賊作業(yè)區(qū)域的SST分布范圍為23～31℃，CPUE主要集中分布在26～28.5 ℃之間，CPUE在23～28 ℃范圍內(nèi)呈小幅度下降趨勢(shì)，至28 ℃達(dá)到最小，隨后逐漸上升，于30.5 ℃達(dá)到最大值。CPUE隨WS增大而緩慢下降趨勢(shì)，主要集中分布在5～7 m/s之間。CPUE隨U的增大呈波浪式下降趨勢(shì)，在0.05 m/s附近達(dá)到最大值，主要集中分布在0.01～0.05 m/s之間(圖3)。

2.4 名義CPUE和標(biāo)準(zhǔn)化CPUE比較

從圖4可知，名義CPUE和基于GAM模型標(biāo)準(zhǔn)化處理后的CPUE變化趨勢(shì)隨著年份的增加呈上升趨勢(shì)，于2019年達(dá)到最高值；整體上，標(biāo)準(zhǔn)化CPUE均略高于名義CPUE，且標(biāo)準(zhǔn)化后的CPUE變化趨勢(shì)與名義CPUE變化趨勢(shì)一致(圖4a)。

2017年9—12月，2018年4—5月、9月、11—12月，2019年1—4月名義CPUE高于標(biāo)準(zhǔn)化CPUE；其余時(shí)間段內(nèi)標(biāo)準(zhǔn)化CPUE均高于或接近對(duì)應(yīng)的名義CPUE。其中2018年11月，名義CPUE和標(biāo)準(zhǔn)化CPUE都達(dá)到最大值；2017年2月，名義CPUE達(dá)到最小值，2018年7月，標(biāo)準(zhǔn)化CPUE達(dá)最小值。總體上名義CPUE的變化波動(dòng)幅度較大，而標(biāo)準(zhǔn)化CPUE變化幅度較小，且名義CPUE與標(biāo)準(zhǔn)化后的CPUE變化趨勢(shì)基本一致(圖4b)。

圖4 基于GAM模型的印度洋北部鳶烏賊年(a)和月(b)平均標(biāo)準(zhǔn)化CPUE

3 討論與分析

3.1 時(shí)空因子對(duì)CPUE的影響

從時(shí)間因子來(lái)看，GAM模型結(jié)果表明，年對(duì)CPUE的貢獻(xiàn)率為5.22%。其中，2018年CPUE的影響效應(yīng)比其余年份高，主要是因?yàn)?018年的生產(chǎn)總量和頻次數(shù)據(jù)相對(duì)較多，因此所統(tǒng)計(jì)的總產(chǎn)量較其他年份多，對(duì)資源密影響相對(duì)較大。同時(shí)月間波動(dòng)性較大，10月到次年1月對(duì)CPUE的影響較大。鳶烏賊為一年生群體，全年可捕撈作業(yè)，一方面，夏季受西南季風(fēng)的高強(qiáng)度且持久性的影響，索馬里沿岸上升流不斷將海底豐富的營(yíng)養(yǎng)鹽帶到海表層，使得浮游動(dòng)植物大量繁殖，當(dāng)上升流到達(dá)索馬里半島外的大渦旋處，在季風(fēng)和大渦旋的共同作用下產(chǎn)生強(qiáng)大的離岸流，將近海岸豐富的營(yíng)養(yǎng)鹽和葉綠素以及沿海獨(dú)特的浮游植物運(yùn)輸?shù)桨⒗Ｖ胁亢Ｓ?，此過(guò)程存在一定的滯后性，所以秋冬季季節(jié)生產(chǎn)力也較高[24-25]；另一方面，9月以后，隨著夏季風(fēng)的消退、相關(guān)的海洋和氣候環(huán)境的變化，包括云量減少、風(fēng)速減弱、降雨減少等因素的綜合影響下，秋冬季的初級(jí)生產(chǎn)力和次級(jí)生產(chǎn)力顯著增加，生物量也顯著升高[26]。此外，夏季受西南季風(fēng)的影響，印度洋北部海域風(fēng)浪較大，不利于捕撈作業(yè)，同時(shí)由于2017年6—7月漁場(chǎng)轉(zhuǎn)移其其他海域，導(dǎo)致印度洋北部海域無(wú)生產(chǎn)數(shù)據(jù)，造成該兩月對(duì)CPUE的影響程度較低。

從空間因子來(lái)看，GAM模型結(jié)果表明，緯度對(duì)CPUE的貢獻(xiàn)率為4.94%，空間范圍在14°N～18°N、59°E～63°E尤其對(duì)CPUE影響較大。該海域受索馬里海流和季風(fēng)海流的影響，營(yíng)養(yǎng)鹽豐富，表層藻類(lèi)繁盛，葉綠素濃度增高，提供了豐富的餌料，是鳶烏賊的一個(gè)重要漁場(chǎng)所在地[27-28]。根據(jù)本文統(tǒng)計(jì)數(shù)據(jù)可知，漁船作業(yè)位置主要集中在16°N～18°N，61.5°E～63°E海域，該海域是鳶烏賊的中心漁場(chǎng)，資源量豐富[29]。陳新軍等[30]研究也認(rèn)為印度洋鳶烏賊漁場(chǎng)在12°N～18°N、58°E～61°E具有較高的CPUE，本文研究結(jié)果與其結(jié)論較一致。

3.2 環(huán)境因子對(duì)CPUE的影響

在所研究的環(huán)境因子中，SST對(duì)CPUE的影響最大。鳶烏賊作為暖水性的大洋性頭足類(lèi)，其漁場(chǎng)分布與SST有著高度的關(guān)聯(lián)性[31]。陳新軍等[3]研究認(rèn)為，印度洋中心漁場(chǎng)的水溫為26.4～29.0 ℃；余為等[32]在印度洋西北海域鳶烏賊的棲息地研究中，得出鳶烏賊的最適SST范圍為27～29 ℃。本文研究中，作業(yè)溫度集中在26～28.5 ℃，與上述研究結(jié)果較一致。與此同時(shí)，較大尺度氣候因子也影響著SST的變化[33]。其中印度洋偶極子(Indian Ocean Dipole,IOD)是影響印度洋氣候環(huán)境變化的重要特殊氣候現(xiàn)象：當(dāng)IOD發(fā)生時(shí)，印度洋處于海溫異常狀態(tài)，主要通過(guò)?！?dú)庀嗷プ饔糜绊懹《妊蠛Ｅ韬Ｑ蟓h(huán)境，IOD處于正相位時(shí)期，西印度洋海溫高于東印度洋，且印度洋季風(fēng)環(huán)流減弱；相反，處于負(fù)相位時(shí)期，西印度洋海溫低于東印度洋，季風(fēng)環(huán)流增強(qiáng)[34]。Wei等[35]研究表明，赤道印度洋海域的印度洋偶極子指數(shù)(DMI)從2017—2019年經(jīng)歷了正負(fù)波動(dòng)到最大正值的過(guò)程，IOD從負(fù)相位轉(zhuǎn)變?yōu)檎辔?，西印度洋海表面溫度逐漸升高。同時(shí)Anderson等[36]的研究結(jié)果也表明這種偶極子分布形態(tài)存在于海表溫度的變化中。

另外海流和風(fēng)速對(duì)印度洋鳶烏賊CPUE也有著較大影響，這在本研究的GAM模型結(jié)果里也有所體現(xiàn)。印度洋北部被非洲大陸、亞洲大陸三面環(huán)抱，使得海洋環(huán)境和氣候變化深受大陸的影響，成為世界典型的季風(fēng)海域，季風(fēng)是驅(qū)動(dòng)印度洋環(huán)境變化和海洋生態(tài)動(dòng)力過(guò)程的最主要因素之一[37]。印度洋北部受氣壓帶和風(fēng)帶位置的季節(jié)性移動(dòng)，以及海陸熱力性質(zhì)差異的影響，海域多為風(fēng)成海流和上升流，進(jìn)而海域表層環(huán)流促使深層和表層的海水進(jìn)行交換，將海水深層豐富的營(yíng)養(yǎng)鹽帶到海表層，帶來(lái)較高的初級(jí)生產(chǎn)力，使得浮游動(dòng)植物生物量大量富集，并呈現(xiàn)出明顯的區(qū)域差異和季節(jié)變化[38]，這些因素都會(huì)對(duì)鳶烏賊資源量和空間分布產(chǎn)生潛在影響。此外，由于作業(yè)方式的差異，風(fēng)速的大小對(duì)魷釣作業(yè)方式影響很大，風(fēng)速大會(huì)使釣具上下浮動(dòng)頻率變大，影響捕撈效率[39]。但對(duì)燈光敷網(wǎng)和燈光罩網(wǎng)捕撈作業(yè)影響較小。全球氣候變化對(duì)印度洋海洋環(huán)境變化產(chǎn)生重要影響，影響的因素也是多方面的。因此，可以考慮將U和WS作為突出的潛在因子應(yīng)用于后續(xù)的漁情預(yù)報(bào)中。

3.3 不同作業(yè)方式對(duì)CPUE的影響

以往的研究中，由于作業(yè)方式較為單一，往往會(huì)忽略其他因素對(duì)CPUE所產(chǎn)生的影響，而在印度洋鳶烏賊漁業(yè)中，存在三種作業(yè)方式，不同作業(yè)方式的CPUE與單位漁船數(shù)量、單位漁船主機(jī)功率、漁具作業(yè)性能、捕撈作業(yè)技術(shù)水平和作業(yè)時(shí)間等因素有關(guān)[39-40]。本研究中燈光敷網(wǎng)和燈光罩網(wǎng)作業(yè)方式較相似，因此兩者之間差異較小。然而，魷釣作業(yè)對(duì)CPUE的影響與前兩者相比有著較大的差異，GAM模型中的結(jié)果也說(shuō)明了這一問(wèn)題。燈光敷網(wǎng)和燈光罩網(wǎng)為主動(dòng)捕撈，主要的捕撈對(duì)象是趨光性的魷魚(yú)和中上層魚(yú)類(lèi)，漁獲選擇性較低，作業(yè)位置轉(zhuǎn)移較為靈活，捕撈努力量較高，且捕撈方式受到海洋環(huán)境變化的影響較小，因此捕撈效率的變化較小，CPUE起伏波動(dòng)較?。欢溼灥淖鳂I(yè)方式相對(duì)被動(dòng)，且捕撈位置相對(duì)固定，漁獲選擇性高，易受捕撈對(duì)象資源量變動(dòng)的制約，且釣具易受海流、風(fēng)浪的影響，因此導(dǎo)致捕撈效率變化較大，CPUE起伏波動(dòng)較大。因此，后續(xù)的相關(guān)研究中，不僅需要考慮時(shí)空環(huán)境因子對(duì)印度洋鳶烏賊CPUE的影響，也不應(yīng)忽視捕撈方式所帶來(lái)的影響。

4 結(jié)論

本文根據(jù)2017—2019年印度洋北部鳶烏賊漁業(yè)生產(chǎn)數(shù)據(jù)，通過(guò)使用GLM和GAM模型，分析時(shí)空因子和重要的環(huán)境因子對(duì)印度洋北部鳶烏賊CPUE變化的影響，以期為科學(xué)管理和利用鳶烏賊漁業(yè)資源提供理論依據(jù)。在后續(xù)的研究中，應(yīng)考慮因子間的綜合效應(yīng)以及加強(qiáng)時(shí)間序列的樣本采集工作；同時(shí)在比較不同作業(yè)方式的差異時(shí)，還應(yīng)該考慮到作業(yè)方式對(duì)資源量反映的準(zhǔn)確性，結(jié)合多種方法綜合考慮資源量和環(huán)境因子間的關(guān)系。