馬有成，朱國(guó)平，張 健,3，王 嘯，張鴻霖，石建高

1.上海海洋大學(xué) 海洋科學(xué)學(xué)院，上海 201306

2.宏?yáng)|漁業(yè)股份有限公司，福建 福州 350015

3.國(guó)家遠(yuǎn)洋漁業(yè)工程技術(shù)研究中心/大洋漁業(yè)資源可持續(xù)開(kāi)發(fā)教育部重點(diǎn)實(shí)驗(yàn)室/農(nóng)業(yè)農(nóng)村部大洋漁業(yè)開(kāi)發(fā)重點(diǎn)實(shí)驗(yàn)室，上海 201306

4.上海開(kāi)創(chuàng)遠(yuǎn)洋漁業(yè)有限公司，上海 200082

5.中國(guó)水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所，上海 200090

鰹 (Katsuwonuspelamis) 是廣泛分布于世界各大洋熱帶和亞熱帶海域的中上層魚(yú)類，具有高度洄游性、漁業(yè)作業(yè)類型多、管理區(qū)域復(fù)雜等特點(diǎn)，是金槍魚(yú)漁業(yè)的主要目標(biāo)種類之一[1]。鰹種群分布模式不是隨機(jī)變化的，而是受到海洋環(huán)境因子和氣候因素變化的影響。海洋環(huán)境因子，如海表面溫度(Sea surface temperature,SST)[2-3]、海表面鹽度 (Sea surface salinity,SSS)[4]和葉綠素a濃度 (Chlorophyllaconcentration,Chl-a)[5]不僅是鰹種群分布的限制性因素，甚至影響其生存。而氣候現(xiàn)象，如厄爾尼諾和拉尼娜現(xiàn)象，也會(huì)對(duì)其產(chǎn)生消極或積極的影響[6-8]。因此，掌握鰹資源豐度的變化規(guī)律，對(duì)準(zhǔn)確預(yù)測(cè)其種群分布、評(píng)估資源量、合理開(kāi)發(fā)和持續(xù)利用資源具有重大意義[9]。

由于氣候和海洋環(huán)境變化，Putri 和Zainuddin[10]預(yù)計(jì)鰹的豐度將出現(xiàn)小幅度下降，因此，全球?qū)灊O業(yè)資源的關(guān)注進(jìn)一步提升。在太平洋[11-12]、印度洋[13-14]和大西洋[15]各海域，圍繞多種作業(yè)方式，各國(guó)學(xué)者均開(kāi)展了氣候變化和海洋環(huán)境變動(dòng)對(duì)鰹資源豐度的影響相關(guān)研究，然而鮮有文獻(xiàn)提及或關(guān)注到降水等因素對(duì)于鰹等大洋性高度洄游物種種群分布的影響。大洋中降水量的變化是對(duì)厄爾尼諾和拉尼娜現(xiàn)象影響的直觀表現(xiàn)，異常的溫度變化會(huì)改變極端降水事件的發(fā)生頻率、強(qiáng)度和持續(xù)時(shí)間[16-18]，而降水對(duì)不同海洋生物的繁殖[19-20]、生活史變化[21]、遷徙行為[22-23]等均具有不同程度的影響，因此忽視降水對(duì)此類生物資源豐度的影響，可能無(wú)法全面、真實(shí)地反映漁業(yè)資源量及其與多種海洋環(huán)境因子的關(guān)系。

為了解和掌握降水對(duì)不同海域鰹群體分布的影響，本文使用毛里塔尼亞海域雙拖漁業(yè)數(shù)據(jù)、中西太平洋圍網(wǎng)和延繩釣漁業(yè)數(shù)據(jù)，結(jié)合降水、SSS、SST 等主要環(huán)境因子，分析不同作業(yè)方式下鰹單位捕撈努力量漁獲量 (Catch per unit effort,CPUE) 的時(shí)空分布特征，以及降水等海洋環(huán)境因子對(duì)不同作業(yè)方式下鰹CPUE 的影響，以期初步探明降水對(duì)不同鰹群體資源分布及變動(dòng)的影響機(jī)制，為今后進(jìn)一步開(kāi)展鰹漁業(yè)資源評(píng)估和養(yǎng)護(hù)，實(shí)現(xiàn)漁業(yè)可持續(xù)發(fā)展提供科學(xué)依據(jù)。

1 材料與方法

1.1 漁業(yè)數(shù)據(jù)

雙拖鰹數(shù)據(jù)來(lái)源于宏?yáng)|國(guó)際 (毛塔) 漁業(yè)發(fā)展有限公司雙拖網(wǎng)中上層作業(yè)船隊(duì)統(tǒng)計(jì)的2017—2019 年的生產(chǎn)漁撈日志，涉及生產(chǎn)漁船共11 艘，漁船作業(yè)范圍為20°W—10°W，16°N—21°N，漁撈日志數(shù)據(jù)包含漁船作業(yè)時(shí)間、作業(yè)位置、放網(wǎng)數(shù)量、鰹漁獲量等。圍網(wǎng)鰹數(shù)據(jù)來(lái)源于上海開(kāi)創(chuàng)遠(yuǎn)洋漁業(yè)有限公司圍網(wǎng)漁船統(tǒng)計(jì)的2018—2020 年的生產(chǎn)漁撈日志，涉及生產(chǎn)漁船共7 艘，漁船作業(yè)范圍為155°E—180°E、180°W—165°W，10°S—6°N，漁撈日志數(shù)據(jù)包含漁船放網(wǎng)時(shí)間、收網(wǎng)時(shí)間、作業(yè)位置、鰹漁獲量等。延繩釣鰹數(shù)據(jù)來(lái)源于美洲熱帶金槍魚(yú)委員會(huì)網(wǎng)站http://www.iattc.org，時(shí)間為2010—2020 年，生產(chǎn)漁船為日本延繩釣漁船，作業(yè)范圍為77°W—148°W，37.5°S—32.5°N；空間分辨率為5°×5°，時(shí)間分辨率為月；數(shù)據(jù)內(nèi)容包括作業(yè)時(shí)間、作業(yè)位置、鉤數(shù)、鰹漁獲量。3 種鰹捕撈方式作業(yè)范圍如圖1 所示。

圖1 鰹漁業(yè)作業(yè)范圍Fig.1 Scope of operation of skipjack tuna fishery

1.2 環(huán)境數(shù)據(jù)

海洋環(huán)境因子中海表面溫度、海表面鹽度數(shù)據(jù)來(lái)自哥白尼數(shù)據(jù)中心網(wǎng)站https://data.marine.copernicus.eu/，空間范圍為180°W—180°E，33°S—33°N，空間分辨率為0.25°×0.25°，時(shí)間分辨率為月；降水?dāng)?shù)據(jù)來(lái)自NCEP/NCAR (美國(guó)國(guó)家環(huán)境預(yù)測(cè)中心/國(guó)家大氣研究中心) 的再分析數(shù)據(jù)集網(wǎng)站https://psl.noaa.gov/，空間范圍為180°W—180°E，90°S—90°N，空間分辨率為2.5°×2.5°，時(shí)間分辨率為月；利用克里金空間插值法[24]對(duì)降水?dāng)?shù)據(jù)進(jìn)行了不同空間分辨率重采樣，根據(jù)已知的大尺度數(shù)據(jù)點(diǎn)賦權(quán)重預(yù)測(cè)小尺度未知點(diǎn)的值。

1.3 數(shù)據(jù)處理

1.3.1 CPUE 標(biāo)準(zhǔn)化

使用CPUE 作為指示鰹豐度指標(biāo)[11]。對(duì)于3 種捕撈方式鰹數(shù)據(jù)計(jì)算月平均CPUE 以指示資源相對(duì)豐度，同時(shí)為了匹配環(huán)境數(shù)據(jù)，使用一個(gè)漁區(qū)單元(0.25°×0.25°、0.25°×0.25°、5°×5°) 內(nèi)月度總漁獲量除以捕撈努力量來(lái)分別計(jì)算雙拖、圍網(wǎng)和延繩釣鰹CPUE，即：

式中：Yy,m,i表示單元漁區(qū)內(nèi)的月度鰹CPUE；∑Cy,m,i表示單元漁區(qū)內(nèi)的月度總漁獲量；∑Ey,m,i表示單元漁區(qū)內(nèi)月度總捕撈努力量；y為年份；m為月份；i為漁區(qū)單元。

1.3.2 統(tǒng)計(jì)方法

廣義加性模型 (Generalized Additive Model,GAM) 是一種非參數(shù)化的廣義多元線性回歸方法，可以直接處理響應(yīng)變量與多個(gè)解釋變量之間的非線性關(guān)系，有效地提升研究精度[25]。利用廣義加性模型分析3 種捕撈作業(yè)方式鰹CPUE 與降水等環(huán)境因子的關(guān)系，即以鰹CPUE 作為響應(yīng)變量，以年度(y) 和月份 (m) 作為時(shí)間因子解釋變量，以漁區(qū)經(jīng)度 (Llon) 和緯度 (Llat) 作為空間因子解釋變量，降水量 (Pprec,mm?d-1)、海表面溫度 (TSS,℃)、海表面鹽度 (SSS,‰) 作為環(huán)境因子解釋變量，建立模型如下：

式中：為防止零值的情況，采用YCPUE加上1 再進(jìn)行對(duì)數(shù)化處理，S表示自然立方樣條平滑；ε表示模型誤差。

利用方差膨脹因子 (Variance inflation factor,VIF) 對(duì)時(shí)空和環(huán)境因子進(jìn)行多重共線性檢驗(yàn)，且VIF 閾值小于10 時(shí)認(rèn)為影響因子間不存在多重共線[26]。采用逐步加入解釋變量的方法，結(jié)合解釋變量擬合的顯著性水平和赤池準(zhǔn)則信息 (Akaike information criterion,AIC) 對(duì) GAM 完整模型進(jìn)行解釋變量篩選；AIC 值越小，模型擬合效果越好，以此獲得最優(yōu)模型[27]。

以上模型的構(gòu)建在R 4.1.3 軟件中實(shí)現(xiàn)，通過(guò)R 語(yǔ)言mgcv 包中的gam 函數(shù)實(shí)現(xiàn)模型擬合和car 軟件包的VIF 函數(shù)判斷因子的多重共線性。

2 結(jié)果

2.1 時(shí)空和環(huán)境因子的共線性

時(shí)空和環(huán)境因子的共線性分析結(jié)果如表1 所示，所選取的時(shí)空和環(huán)境因子 VIF 值均小于10，表明各因子之間不存在共線性問(wèn)題，均可作為解釋變量用于GAM 模型分析。

表1 影響因子之間的共線性分析Table 1 Collinearity analysis among influencing factors

2.2 最優(yōu)GAM 模型

將解釋變量逐步加入模型，分別構(gòu)建得到關(guān)于3 種作業(yè)方式鰹CPUE 與各解釋變量的GAM 模型，獲得GAM 模型統(tǒng)計(jì)參數(shù)如表2 所示。雙拖作業(yè)方式最優(yōu)模型保留所有變量；圍網(wǎng)作業(yè)方式最優(yōu)模型首先剔除不顯著解釋因子月份，再進(jìn)行模型構(gòu)建；延繩釣作業(yè)方式首先剔除不顯著解釋因子經(jīng)度和降水，再進(jìn)行模型構(gòu)建。

表2 GAM 模型擬合及最優(yōu)模型篩選Table 2 GAM model fitting and optimal model screening

2.3 雙拖捕撈作業(yè)鰹CPUE 與各因子關(guān)系

2.3.1 雙拖鰹CPUE 分布時(shí)空特征

雙拖捕撈作業(yè)方式捕撈到鰹的GAM 模型結(jié)果如表3 和圖2 所示。由表3 可見(jiàn)，時(shí)空因子 (月份、經(jīng)度和緯度) 均與鰹CPUE 呈極顯著關(guān)系 (P<0.001)。由圖2 可見(jiàn)，鰹CPUE 在10 月到翌年3 月增長(zhǎng)速度較快，在3 月達(dá)到最高值，之后開(kāi)始下降至5 月，而后5—10 月出現(xiàn)先上升再下降的趨勢(shì)；在經(jīng)度上，鰹CPUE 在17.5°W—17.0°W 上，自西向東呈遞減的趨勢(shì)；而后出現(xiàn)上升趨勢(shì)，在16.5°W時(shí)達(dá)到最大，隨后出現(xiàn)自西向東遞減趨勢(shì)；緯度與鰹CPUE 存在顯著正相關(guān)關(guān)系，即其CPUE 自南向北逐漸上升。

表3 各影響因子對(duì)3 種作業(yè)方式鰹CPUE 的影響Table 3 Effects of environmental factors on CPUE of three fishing methods for skipjack

圖2 雙拖鰹CPUE 與各因子關(guān)系Fig.2 Relationship between CPUE and various factors for skipjack of pair trawl

2.3.2 雙拖鰹CPUE 與環(huán)境因子關(guān)系

環(huán)境因子 (降水、SSS 和SST) 均與鰹CPUE 呈極顯著性相關(guān) (P<0.001)。鰹漁場(chǎng)降水范圍主要介于0～10 mm?d-1，降水為0～5 mm?d-1時(shí)，鰹CPUE與其呈正相關(guān)；介于5～7 mm?d-1時(shí)，降水與CPUE呈負(fù)相關(guān)；而后隨著降水的增加鰹CPUE 升高(圖2-d)。鰹漁場(chǎng)SSS 范圍主要介于35‰～37‰，鰹CPUE 隨著SSS 的增加而上升，SSS 為36‰ 時(shí)CPUE 最高，其后開(kāi)始下降 (圖2-e)；SST 在17～18 ℃范圍內(nèi)，鰹CPUE 隨著溫度的升高而降低，在18～20 ℃范圍內(nèi)，鰹CPUE 與SST 呈正相關(guān)關(guān)系，而超出此范圍時(shí)，鰹CPUE 隨著SST 的增加而稍微下降。

2.4 圍網(wǎng)捕撈作業(yè)鰹CPUE 與各因子關(guān)系

2.4.1 圍網(wǎng)鰹CPUE 分布時(shí)空特征

用圍網(wǎng)捕撈的鰹的GAM 模型結(jié)果如表3 和圖3 所示。時(shí)間因子 (月份) 對(duì)鰹CPUE 無(wú)影響，空間因子 (經(jīng)緯度) 對(duì)鰹CPUE 有顯著性影響 (P<0.05)。在經(jīng)度梯度上，鰹CPUE 呈現(xiàn)出自西向東逐漸降低的趨勢(shì) (圖3-a)；在緯度梯度上，鰹CPUE自南向北先升高再小幅度下降并趨于平緩，再小幅度增加，最后小幅度降低 (圖3-b)。

圖3 圍網(wǎng)鰹CPUE 與各因子的關(guān)系Fig.3 Relationship between CPUE and various factors for skipjack of seine

2.4.2 圍網(wǎng)鰹CPUE 與環(huán)境因子關(guān)系

降水和SST 對(duì)鰹CPUE 的影響顯著 (P<0.05)，SSS 對(duì)鰹CPUE 的影響極其顯著 (P<0.001，表3)。漁場(chǎng)主要降水范圍為0～15 mm?d-1，介于0～3 mm?d-1時(shí)，鰹CPUE 與降水呈正相關(guān)；介于3～6 mm?d-1時(shí)，鰹CPUE 與降水呈負(fù)相關(guān)；超出此范圍，則鰹CPUE 與降水呈正相關(guān) (圖3-c)。SSS 介于34‰～34.5‰ 時(shí)鰹CPUE 較為穩(wěn)定，超過(guò)此范圍后，鰹CPUE 隨著SSS 的增加而升高，而后在34.5‰～35.3‰ 時(shí)達(dá)到最高值，之后逐漸降低 (圖3-d)。SST 與鰹CPUE 呈負(fù)相關(guān)關(guān)系，鰹CPUE 較高的SST 范圍介于23.5～24.5 ℃ (圖3-e)。

2.5 延繩釣捕撈作業(yè)鰹CPUE 與各因子關(guān)系

2.5.1 延繩釣鰹CPUE 分布時(shí)空特征

用延繩釣捕撈的鰹的GAM 模型結(jié)果如表3 和圖4 所示。月份與鰹CPUE 呈極顯著相關(guān)關(guān)系 (表3)，在1 月鰹CPUE 最高，隨之開(kāi)始下降，5—7 月稍有回升而后又繼續(xù)下降，10—12 月有上升的趨勢(shì)(圖4-a)；空間因子緯度對(duì)鰹CPUE 有顯著性影響(表3)，整體上由低緯至高緯，鰹CPUE 逐漸升高(圖4-b)。

圖4 延繩釣鰹CPUE 與各因子關(guān)系Fig.4 Relationship between CPUE and various factors for skipjack of longline

2.5.2 延繩釣鰹CPUE 與環(huán)境因子關(guān)系

降水對(duì)鰹CPUE 無(wú)顯著性影響 (P>0.05)，SSS 和SST 對(duì)鰹CPUE 有顯著性影響 (表3)。鰹漁場(chǎng)SSS 范圍主要為34‰～36‰，鰹CPUE 與SSS 呈正相關(guān)關(guān)系 (圖4-e)；鰹漁場(chǎng)SST 范圍主要為15～29 ℃，鰹CPUE 先隨著溫度升高而增加，當(dāng)超過(guò)22 ℃時(shí)鰹CPUE 有下降趨勢(shì) (圖4-d)。

3 討論

3.1 SSS 與SST 對(duì)3 種作業(yè)方式鰹CPUE 的影響

魚(yú)類通常對(duì)較小范圍的鹽度變化反應(yīng)劇烈，鹽度通過(guò)影響魚(yú)體的滲透壓調(diào)節(jié)來(lái)改變其分布[28]。作為一種暖水性魚(yú)類，金槍魚(yú)的分布和豐度表現(xiàn)出與鹽度變化顯著的相關(guān)關(guān)系[29]，鰹屬于狹鹽性魚(yú)類，喜生活于高鹽海域，其種群分布與SSS 之間存在很強(qiáng)的相關(guān)性[30]。楊勝龍等[31]研究指出鰹生活的鹽度范圍為34.0‰～35.5‰，本研究結(jié)果與其相似，從GAM 模型圖可以看出，鰹喜好的鹽度較高，漁場(chǎng)分布集中在SSS 為34‰～36‰ 范圍內(nèi)，雙拖作業(yè)鰹CPUE 較高時(shí)的SSS 約為36‰，這可能與其生存區(qū)域的海洋環(huán)境有關(guān)，毛里塔尼亞沿岸有加那利寒流自北向南海岸走向流動(dòng)，帶來(lái)北部的冷水，與由南向北流動(dòng)的幾內(nèi)亞暖水在布蘭科海峽區(qū)域混合，形成顯著的涌升流[32]，故此區(qū)域內(nèi)的SSS 高于其他海域。

魚(yú)類的生理活動(dòng)有其最佳的SST 范圍，SST 是影響魚(yú)類繁殖、生長(zhǎng)及游泳能力等的一個(gè)重要因素[29]。已有研究證明，金槍魚(yú)物種的垂直和水平分布由SST 變化所引起[33]，且因生理差異，金槍魚(yú)在不同生命階段對(duì)SST 的敏感程度也不同，幼魚(yú)隨著年齡的增長(zhǎng)對(duì)SST 會(huì)更加敏感。本研究發(fā)現(xiàn)，SST 對(duì)3 種作業(yè)方式鰹的影響極為顯著(P<0.001)，且鰹CPUE 較高時(shí)的SST 不同，雙拖作業(yè)鰹CPUE 在SST 為20 ℃時(shí)最高，圍網(wǎng)作業(yè)鰹CPUE 在SST 為24 ℃時(shí)最高，延繩釣作業(yè)鰹CPUE在SST 為22 ℃時(shí)最高；3 種作業(yè)方式鰹CPUE 有著不同的SST 適宜范圍，但當(dāng)SST 超過(guò)最適區(qū)間后其CPUE 均隨著SST 的升高而降低。說(shuō)明SST 的變化對(duì)不同海域的鰹資源有著正面和負(fù)面的影響，當(dāng)赤道區(qū)溫度升至33～34 ℃時(shí)，將超過(guò)鰹的耐受范圍 (其耐受范圍介于20～32 ℃)[14]，此時(shí)SST 對(duì)鰹的生存和分布產(chǎn)生了負(fù)面影響。

3.2 降水對(duì)3 種作業(yè)方式鰹CPUE 的影響

鰹全年可以進(jìn)行產(chǎn)卵，繁殖潛力大，每批產(chǎn)卵約100 萬(wàn)枚[34]，其繁殖所需的能量主要來(lái)自食物攝入[35]，因此鰹的特點(diǎn)是生長(zhǎng)快速和與環(huán)境密切相關(guān)的產(chǎn)卵行為。鰹遷移距離大，水平遷移距離平均超過(guò)1 000 n mile[36]，與其他金槍魚(yú)物種相比，鰹的垂直運(yùn)動(dòng)有限，僅限于中上層水層，因?yàn)槠鋵?duì)低溶解氧和低溫耐受性有限[37]；其體溫調(diào)節(jié)能力也較有限，僅適合棲息于相對(duì)溫暖、含氧量高的水層，因此溫躍層深度變化會(huì)對(duì)鰹的捕撈產(chǎn)生影響；海洋環(huán)流的變化可能會(huì)影響海洋中餌料，從而影響鰹的小規(guī)模分布[38]。鰹對(duì)海洋環(huán)境條件變化十分敏感[6]，種群分布極易受其影響，降水可以通過(guò)影響海洋初級(jí)生產(chǎn)中的有機(jī)物，導(dǎo)致養(yǎng)分輸送和初級(jí)生產(chǎn)者物種組成產(chǎn)生變化[39-40]，通過(guò)改變淡水輸送的幅度和時(shí)間，影響有機(jī)物的輸送和海域處理有機(jī)物的停留時(shí)間和環(huán)境條件 (光照、分層、微生物群落等)[41-42]，還可能通過(guò)改變沉積物、營(yíng)養(yǎng)物質(zhì)和有機(jī)物的輸送以及改變分層和溶解氧濃度等影響海洋的物理和生物地球化學(xué)環(huán)境[43]，從而影響鰹的繁殖、遷徙和攝食來(lái)影響其資源的分布。

本研究發(fā)現(xiàn)降水對(duì)3 種作業(yè)方式鰹CPUE 的影響不同。降水對(duì)雙拖和圍網(wǎng)作業(yè)鰹CPUE 的影響顯著，且對(duì)鰹CPUE 的影響較為一致，這2 種作業(yè)方式的鰹均存在適宜生存的降水范圍；而降水對(duì)延繩釣作業(yè)鰹CPUE 的影響不顯著 (表3)，這可能是因?yàn)榻邓畬?duì)不同作業(yè)方式區(qū)域內(nèi)鰹的棲息地、行為和海洋環(huán)境等影響不同所致。本研究中，雙拖作業(yè)方式捕撈的鰹在毛里塔尼亞沿岸海域，存在著幾內(nèi)亞暖流和加那利寒流相交匯的情況，通常寒暖流的交匯會(huì)使得該區(qū)域各種漁業(yè)資源較為豐富[44]，這些環(huán)境條件有利于幼魚(yú)存活。降雨的變化會(huì)影響沿海地區(qū)的鹽度、水渾濁度和殺蟲(chóng)劑等污染物的濃度，進(jìn)而影響生物的繁殖行為，其對(duì)高繁殖力和高生長(zhǎng)速度的鰹有顯著性影響[45]。本研究結(jié)果顯示，中西太平洋圍網(wǎng)鰹CPUE 與降水存在顯著相關(guān)關(guān)系，而東太平洋延繩釣鰹CPUE 與降水無(wú)顯著相關(guān)關(guān)系，這可能由鰹的生物學(xué)特性所致，因鰹在洄游過(guò)程中會(huì)受到海洋壞境變化的性影響[7]，對(duì)于中西太平洋熱帶海域，因SST 常年在28 ℃以上，形成了西太平洋暖池，其提供了太平洋最大比例的金槍魚(yú)捕撈量，是鰹理想的索餌場(chǎng)[46]，且在生物生產(chǎn)力方面顯示出顯著的動(dòng)態(tài)變化；而對(duì)于東太平洋海域，除了最強(qiáng)厄爾尼諾事件期間外，初級(jí)生產(chǎn)力和“下游”金槍魚(yú)餌料豐度在空間和時(shí)間上都很高且穩(wěn)定。與東太平洋相比，中西太平洋暖池的次級(jí)生產(chǎn)力較低但充足，可補(bǔ)償溫度、氧氣和水透明度等環(huán)境條件的不足[47]，造成降水對(duì)兩個(gè)海域鰹CPUE的影響有差異。

降水對(duì)3 種作業(yè)方式鰹CPUE 影響的差異，還有可能受作業(yè)方式的影響，如捕撈效率差異、捕撈鰹生長(zhǎng)發(fā)育階段差異、漁業(yè)空間分辨率差異以及捕撈海域餌料差異等。雙拖和圍網(wǎng)漁具捕撈鰹效率高，以鰹作為主捕目標(biāo)，在空間和時(shí)間分布上較為有限[48]；延繩釣漁具捕撈鰹效率較低，鰹作為兼捕物種，其CPUE 值較低，因此數(shù)據(jù)的穩(wěn)定性不高，可能對(duì)本研究中降水與鰹CPUE 的相關(guān)性產(chǎn)生了影響。本文采用延繩釣漁業(yè)數(shù)據(jù)，由于不同作業(yè)方式的模型處理和CPUE 計(jì)算方法一致，且延繩釣漁業(yè)可提供大空間尺度和時(shí)間尺度相對(duì)一致的數(shù)據(jù)[49]，因此仍是優(yōu)良的數(shù)據(jù)來(lái)源；日本延繩釣船隊(duì)主要捕撈對(duì)象為各個(gè)魚(yú)種中體型較大的個(gè)體，因此本文使用的是鰹生長(zhǎng)發(fā)育成熟階段的數(shù)據(jù)，預(yù)測(cè)的資源分布也應(yīng)被視為成熟鰹的分布變化，此階段鰹的游泳和覓食能力較強(qiáng)。有研究發(fā)現(xiàn)叉長(zhǎng)大于45 cm 的鰹才可以在黑潮充分覓食后，遷移到黑潮最北端[50]，圍網(wǎng)和雙拖作業(yè)方式捕撈的鰹集群性較高，且捕撈的幼魚(yú)對(duì)降水影響的響應(yīng)較高，這可能導(dǎo)致降水對(duì)延繩釣作業(yè)方式鰹CPUE 的影響不顯著；本文使用的雙拖和圍網(wǎng)作業(yè)鰹數(shù)據(jù)均為海上實(shí)測(cè)漁業(yè)數(shù)據(jù)，空間分辨率較小 (0.25°)，而延繩釣作業(yè)鰹數(shù)據(jù)為美洲熱帶金槍魚(yú)委員會(huì)網(wǎng)站記錄的組織數(shù)據(jù)，空間分辨率較大 (5°)，隨著空間尺度的增大，各個(gè)小漁區(qū)所包括的海域面積增加，漁區(qū)內(nèi)所包含的環(huán)境信息也增多[51]，可能會(huì)影響降水與鰹CPUE 的相關(guān)性；在延繩釣漁業(yè)作業(yè)過(guò)程中，用釣線結(jié)縛裝餌料的釣鉤來(lái)引誘金槍魚(yú)吞食[52]，會(huì)造成捕撈海域的鰹餌料較多，從而導(dǎo)致鰹CPUE 受降水的影響不顯著。

4 小結(jié)

本研究使用GAM 模型對(duì)降水、SSS、SST 海洋環(huán)境因子與鰹CPUE 關(guān)系進(jìn)行分析，結(jié)果表明降水與常規(guī)海洋環(huán)境因子一樣會(huì)對(duì)鰹資源變動(dòng)產(chǎn)生影響。建議在今后分析環(huán)境因子對(duì)漁業(yè)CPUE 的影響效果時(shí)，將降水納入常規(guī)海洋環(huán)境因子中。后續(xù)研究需要探討更多的環(huán)境因子，如葉綠素a濃度、海表面高度等，以及大尺度氣候變化對(duì)漁業(yè)資源的影響，從而更全面地解釋物種資源豐度變動(dòng)與環(huán)境因子間的關(guān)系。