肖 戈，徐 博，張 衡，唐峰華，陳 峰，朱文斌

1. 中國水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所/農(nóng)業(yè)農(nóng)村部遠(yuǎn)洋與極地漁業(yè)創(chuàng)新重點(diǎn)實(shí)驗(yàn)室，上海 200090

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

3. 浙江省海洋水產(chǎn)研究所，浙江 舟山 316021

鳶烏賊 (Symplectoteuthis oualaniensis) 隸屬于頭足綱、槍形目、柔魚科，廣泛分布在印度洋和太平洋海域[1]，是這些海域重要的經(jīng)濟(jì)漁業(yè)資源之一，也是目前中國遠(yuǎn)洋漁船在西北印度洋公海即阿拉伯外海的主捕對象。21世紀(jì)初期中國開始對阿拉伯外海鳶烏賊開展試驗(yàn)性探捕和資源調(diào)查，發(fā)現(xiàn)阿拉伯海外海的北部鳶烏賊資源較豐富[2]，但隨后并未對其進(jìn)行商業(yè)型規(guī)模的開發(fā)和利用。直至2013年前后中國才有漁船在阿拉伯海外海進(jìn)行正式的商業(yè)性捕撈開發(fā)。目前中國在西北印度洋公海海域進(jìn)行商業(yè)捕撈的漁船近200艘，并且有逐年增加的趨勢。對阿拉伯海外海鳶烏賊資源進(jìn)行合理開發(fā)及可持續(xù)利用，有利于提升中國在印度洋海域漁業(yè)科學(xué)研究中的地位和影響力，同時(shí)對緩解近岸資源捕撈壓力具有積極意義。

近年來，印度洋鳶烏賊的研究主要集中在漁業(yè)生態(tài)學(xué)方面，如角質(zhì)顎和耳石等[3-6]，而環(huán)境因素對鳶烏賊的影響研究區(qū)域集中在中國南海[7-10]。有研究根據(jù)印度洋北部捕撈鳶烏賊3種方式 (燈光敷網(wǎng)、燈光罩網(wǎng)和魷釣) 的作業(yè)次數(shù)、漁獲量和單位捕撈努力量漁獲量 (Catch per unit effort, CPUE) 進(jìn)行統(tǒng)計(jì)，對漁場時(shí)空分布及其年間差異進(jìn)行分析[11]，但缺少對印度洋公海海域鳶烏賊漁場海洋環(huán)境變動的相關(guān)研究；且伴隨著《北太平洋公海漁業(yè)資源養(yǎng)護(hù)和管理公約》的實(shí)施，更多遠(yuǎn)洋漁船可能會從北太平洋選擇轉(zhuǎn)移至印度洋生產(chǎn)，因此從掌握該海域資源狀況及提升生產(chǎn)穩(wěn)定性等角度出發(fā)，對阿拉伯海外海鳶烏賊的時(shí)空分布和漁場環(huán)境信息的了解尤為必要。本研究利用2019年12月—2020年2月和2020年10—11月2個(gè)航次的調(diào)查資料，選取經(jīng)緯度、農(nóng)歷日、海表溫度 (Sea surface temperature,SST)、葉綠素a(Chla)濃度和海表溫度距平 (Sea surface temperature anomaly, SSTA)，開展上述要素在阿拉伯海外海鳶烏賊漁場的時(shí)空分布研究，旨為掌握印度洋鳶烏賊的漁汛規(guī)律、漁場的發(fā)生機(jī)理，為阿拉伯海外海漁場的精細(xì)化漁場預(yù)報(bào)提供科學(xué)參考和可靠依據(jù)。

1 材料與方法

1.1 調(diào)查要素

第1航次的調(diào)查漁船為燈光敷網(wǎng)漁船“欣海1202”，船總長50.95 m，船型寬8.2 m，型深4.2 m，總噸位 722 t，主機(jī)機(jī) 1 臺，功率 1 260 kW，副機(jī)600、400、250 kW各1臺；網(wǎng)具的網(wǎng)口周長約 420 m，網(wǎng)衣長度 250 m，最大網(wǎng)目尺寸 60 mm，最小網(wǎng)目尺寸30 mm；船上配有75盞集魚燈，功率為4 kW，另配有1盞2 kW的水下導(dǎo)魚燈。第2航次的調(diào)查漁船為燈光罩網(wǎng)漁船“興邦908”，船總長 61.64 m，船型寬 12.2 m，型深 5.2 m，總噸位 1 492 t，主機(jī) 1 臺，功率1 540 kW，副機(jī) 809、809、441 kW各1臺；網(wǎng)具的網(wǎng)口周長624 m，網(wǎng)衣長度約115 m，最大網(wǎng)目尺寸40 mm， 最小網(wǎng)目尺寸35 mm；船上配有300盞2 kW的水上集魚燈，8盞4 kW的水下集魚燈。

第1航次的調(diào)查時(shí)間為2019年12月—2020年2月，調(diào)查站位49個(gè)；第2航次的調(diào)查時(shí)間為2020年10—11月，調(diào)查站位54個(gè)。作業(yè)時(shí)每個(gè)站點(diǎn)燈誘3～5 h，記錄每站下網(wǎng)信息 (經(jīng)緯度、農(nóng)歷日) 和漁獲數(shù)據(jù)。調(diào)查海域位于阿拉伯海外海(60.18°E—65.0°E, 15.87°N—19.70°N) (圖1)。2 個(gè)航次的調(diào)查站點(diǎn)設(shè)計(jì)參考鳶烏賊的歷史中心漁場，并覆蓋了中心漁場的附近海域，以期開拓阿拉伯海外海鳶烏賊漁場范圍。

圖1 阿拉伯海外海鳶烏賊漁場調(diào)查站位分布Fig. 1 Distribution of survey stations in fishing grounds of S. oualaniensis in outer sea of Arabian Sea

1.2 環(huán)境要素與生產(chǎn)調(diào)查數(shù)據(jù)的采集

海洋環(huán)境要素的數(shù)據(jù)采集參照GB 17378.3—1998 《海洋監(jiān)測規(guī)范)》進(jìn)行，當(dāng)場采用CTD(Conductivity, Temperature, Depth system, Sea-Bird 911 plus CTD) 記錄 SST。生產(chǎn)調(diào)查數(shù)據(jù)包括每網(wǎng)次、每次下網(wǎng)的經(jīng)緯度、農(nóng)歷日、鳶烏賊漁獲量等。

1.3 遙感數(shù)據(jù)的獲取

遙感數(shù)據(jù)包括 Chla濃度和 SSTA。Chla濃度數(shù)據(jù)來自美國航空航天局 (National Aeronautics and Space Administration, NASA) 網(wǎng)站 (https://www.nasa.gov/) 提供的中分辨率成像光譜儀 (Moderateresolution Imaging Spectroradiometer, MODIS) 衛(wèi)星傳感器；SSTA來自美國國家海洋和大氣局 (National Oceanic and Atmospheric Administration,NOAA) 的 Ocean Watch 網(wǎng)站 (https://oceanwatch.pifsc.noaa.gov/)。時(shí)間分別率均 1 d，其中 Chla濃度的空間分辨率為 1/24°×1/24° (約 4 km)，SSTA的空間分別率為5 km。將遙感數(shù)據(jù)進(jìn)行0.25°×0.25°的處理后與生產(chǎn)調(diào)查數(shù)據(jù)進(jìn)行對照匹配。

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

鳶烏賊的CPUE為每個(gè)站點(diǎn)的各網(wǎng)次的漁獲量(t·網(wǎng)?1)；根據(jù)農(nóng)歷日和實(shí)際生產(chǎn)設(shè)定：月光夜為農(nóng)歷日的初十至廿二，月黑夜為農(nóng)歷日的廿三至三十、初一至初九；設(shè)定以0.5 ℃為1個(gè)SST區(qū)間，以 0.1 mg·m?3為 1 個(gè) Chla區(qū)間，以 0.2 ℃ 為1個(gè)SSTA區(qū)間。

式中：ycatch為每日站位的總產(chǎn)量 (t)；nnet為每日站位的放網(wǎng)次數(shù)。

廣義可加模型 (Generalized Additive Model,GAM) 是 1986 年 Hastie 和 Tibshirani[12]提出的一種廣義線性模型的非參數(shù)化擴(kuò)展，能直觀地對多個(gè)解釋變量和相應(yīng)變量之間進(jìn)行非線性關(guān)系的處理。本研究將可能影響響應(yīng)變量 (阿拉伯海外海鳶烏賊CPUE) 的解釋變量分為時(shí)空因素和海洋環(huán)境因素，其中時(shí)空因素包括經(jīng)緯度和農(nóng)歷日；海洋環(huán)境要素包括SST、Chla和SSTA。綜上，構(gòu)建的GAM模型如下：

式中：為防止零值出現(xiàn)，采用CPUE加上1后進(jìn)行對數(shù)化處理的方法；s()為樣條平滑 (Smoothing spline)，可用于單變量平滑及多變量的各項(xiàng)同性平滑；lon為經(jīng)度；lat為緯度；lc為農(nóng)歷日；SST (℃)；Chla (mg·m?3)；SSTA (℃)。

利用赤池信息準(zhǔn)則 (Akaike's Information Criterion, AIC) 檢驗(yàn)加入解釋變量后模型的擬合度，選擇最優(yōu)的解釋變量，AIC的值越小，表示模型的擬合效果越好[13]。利用F檢驗(yàn)評估各因子的顯著性。構(gòu)建GAM模型和計(jì)算采用R 4.0.5軟件包的mgcv 函數(shù)庫。繪圖采用 ArcGIS 10.2、Excel 2019和R 4.0.5軟件進(jìn)行。

1.5 最佳模型的篩選

模型中逐步加入各解釋變量后的詳細(xì)變化情況見表1。方差解釋率與決定系數(shù)不斷增大，且最終各因素均對模型存在顯著性影響。lg(CPUE+1)～s(lon, lat)+s(lc)+s(SST)+s(Chla)+s(SSTA) 為最佳模型，解釋率為83.7 %。殘差分布圖 (圖2) 表明最終選擇的最優(yōu)模型能較好地解釋阿拉伯海外海鳶烏賊CPUE的變化。

表1 最優(yōu)廣義可加模型分析結(jié)果Table 1 Analysis results of optimal Generalized Additive Model

圖2 鳶烏賊最佳廣義可加模型殘差圖Fig. 2 Residual plot of optimal Generalized Additive Model

2 結(jié)果

2.1 調(diào)查概況

2個(gè)航次的阿拉伯海外海鳶烏賊資源調(diào)查情況見表2。調(diào)查期間，共捕獲鳶烏賊277.44 t，平均CPUE 為 2.69 t·網(wǎng)?1；第 2 航次的總產(chǎn)量和 CPUE平均值均大于第1航次；2個(gè)航次的SST介于24.56～28.60 ℃，Chla質(zhì)量濃度介于 0.16～1.49 mg·m?3，SSTA 介于?0.45～1.16 ℃。

表2 2個(gè)航次的阿拉伯海外海鳶烏賊資源調(diào)查概況Table 2 Survey of S. oualaniensis resources by two voyages in outer sea of Arabian Sea

2.2 CPUE 分布情況

圖3為2個(gè)航次鳶烏賊CPUE的分布情況。第 1 航次的調(diào)查范圍為 60°E—63.50°E，15.50°N—18.50°N，總漁獲量 113.20 t，平均 CPUE 為 2.31 t·網(wǎng)?1，高 CPUE 區(qū)域主要集中在 60°E—61°E，16°N 和 62°E，16°N—18°N 附近 (圖3-a)；第 1 航次的調(diào)查范圍為 61.50°E—65°E，16°N—20°N，總漁獲量 164.24 t，平均 CPUE 為 3.04 t·網(wǎng)?1，高 CPUE 區(qū)域主要集中在 61.50°E—63°E，17°N—17.75°N 附近 (圖3-b)。

圖3 阿拉伯海外海鳶烏賊單位捕撈努力量漁獲量分布注：a. 第一航次；b. 第二航次；圖6、圖8同此。Fig. 3 CPUE distribution of S. oualaniensis in outer sea of Arabian SeaNote:a. 1st voyage; b. 2nd voyage. The same case in Fig.6 and Fig. 8.

2.3 主要海洋環(huán)境因子與 CPUE 的疊加分布

漁場的SST為24.5～25.5 ℃時(shí)，鳶烏賊的總漁獲量為60.01 t，占總漁獲量的53.0%，CPUE為3.14 t·網(wǎng)?1，是平均 CPUE 的 1.16 倍 (圖4-a)；在26.5～27.5 ℃組內(nèi)，鳶烏賊的漁獲量占總漁獲量的85.6%，共140.7 t，CPUE呈現(xiàn)下降趨勢，最高值出現(xiàn)在 26.5～27.0 ℃ 組內(nèi)，為 8.31 t·網(wǎng)?1(圖4-b)。

圖4 不同海表溫度組的鳶烏賊漁獲量與單位捕撈努力量漁獲量分布Fig. 4 Catch and CPUE of S. oualaniensis in different SST groups

Chla質(zhì)量濃度為 0.3～0.4 mg·m?3時(shí)，漁獲量為85.73 t，占總漁獲量的75.73 %，其對應(yīng)的CPUE為 2.38 t·網(wǎng)?1，CPUE 最高值是 3.73 t·網(wǎng)?1，對應(yīng)的Chla為 0.2～0.3 mg·m?3(圖5-a)，CPUE 與 Chla的疊加圖 (圖6-a) 顯示，CPUE高值區(qū)顏色較深，說明其對應(yīng) Chla范圍小于 0.5 mg·m?3；漁獲量和CPUE最高值均落在 Chla質(zhì)量濃度0.4～0.5 mg·m?3范圍內(nèi)，分別為 140.7 t和 4.69 t·網(wǎng)?1(圖5-b)。因?yàn)槲丛?Chla質(zhì)量濃度為 0.5～0.7 mg·m?3的區(qū)域進(jìn)行調(diào)查，因此該范圍內(nèi)的CPUE與漁獲量均為0。CPUE 高值區(qū)域集中在 0.4～0.8 mg·m?3范圍內(nèi)(圖6-b)。

圖5 不同葉綠素a組的鳶烏賊漁獲量與單位捕撈努力量漁獲量分布Fig. 5 Catch and CPUE of S. oualaniensis in different Chl a groups

圖6 阿拉伯海外海鳶烏賊單位捕撈努力量漁獲量分布與葉綠素 a疊加圖Fig. 6 CPUE distribution and stack distribution of Chl a of S. oualaniensis in Arabian Sea

SSTA介于0.2～0.6 ℃時(shí)，兩組的漁獲量約占總漁獲量一半，CPUE分別為 3.82和 4.22 t·網(wǎng)?1(圖7-a)；漁獲量和CPUE的最高值分別為84.2 t和 9.78 t·網(wǎng)?1，是平均 CPUE 的 3.21 倍 (圖7-b)。2個(gè)航次的高CPUE主要集中在圖中黃色區(qū)域附近，即SSTA為0.2～1.0 ℃的區(qū)域 (圖8)。

圖7 不同海表溫度距平組的鳶烏賊漁獲量與單位捕撈努力量漁獲量Fig. 7 Catch and CPUE of S. oualaniensis in different SSTA groups

圖8 阿拉伯海外海鳶烏賊單位捕撈努力量漁獲量分布與海表溫度距平疊加圖Fig. 8 CPUE distribution and stack distribution of SSTA of S. oualaniensis in Arabian Sea

2.4 最佳模型顯示要素對 CPUE 的影響作用

通過選擇最佳GAM模型分析得到經(jīng)緯度在所有變量中的方差解釋率最高 (34.5%)。圖9-a中黑色實(shí)線的范圍表明阿拉伯海外海鳶烏賊CPUE較高的區(qū)域有 3 處，分別是 60°E—61.50°E，15.90°N—16.50°N；60.80°E—64°E，17°N—17.80°N；62°E—64.8°E，18.8°N—19.10°N；圖9-b顯示，在經(jīng)緯度較小范圍的區(qū)域，隨著經(jīng)緯度變大，對鳶烏賊CPUE的影響逐漸下降，并且在 62°E—63°E，17°E—18°N附近有一穹頂狀突起，對應(yīng)了圖3-a中的1個(gè)CPUE高值密集區(qū) (62°E、16°N—18°N) 和圖 3-b 的CPUE高值密集區(qū)。

圖9 GAM模型中各因子對鳶烏賊單位捕撈努力量漁獲量的影響作用Fig. 9 GAM-derived effects of each factors on CPUE of S. oualaniensis

農(nóng)歷日的方差解釋率為20.0%。圖9-c顯示農(nóng)歷日對鳶烏賊CPUE的影響呈現(xiàn)先升后小幅下降的趨勢，并在農(nóng)歷日為15—23日呈上升趨勢，之后下降。GAM模型最適農(nóng)歷日介于17—25日。2個(gè)航次中，月黑夜的漁獲量、作業(yè)次數(shù)和CPUE均大于月亮夜，第1和第2航次月黑夜的產(chǎn)量分別為月亮夜的3.01和2.00倍，但作業(yè)次數(shù)差別不大(表3)。

表3 農(nóng)歷日鳶烏賊漁獲量和單位捕撈努力量漁獲量情況Table 3 Catch and CPUE of S. oualaniensis on lunar calendar day

SST的方差解釋率為20.0%。影響曲線呈先降低、回升后再降低的趨勢，波峰位于27.0～28.0 ℃ (圖9-d)，置信區(qū)間表明最適SST介于26.5～27.5 ℃。Chla的方差貢獻(xiàn)率最低 (2.9%)，影響曲線呈倒“V”形。當(dāng) Chla為 0.30～0.50 mg·m?3時(shí)，置信度較高，CPUE與Chla呈現(xiàn)正相關(guān)關(guān)系(圖9-e)；當(dāng) Chla大于 0.50 mg·m?3時(shí)，則為負(fù)相關(guān)關(guān)系。SSTA的方差解釋率為15.9%，系環(huán)境因子中解釋率最高的因子。當(dāng)SSTA為0～0.4 ℃時(shí)，影響曲線呈遞增趨勢；當(dāng)SSTA＞0.4 ℃時(shí)，則CPUE與SSTA呈負(fù)相關(guān) (圖9-f)。置信區(qū)間表明最適SSTA范圍在0～0.43 ℃。

3 討論

3.1 阿拉伯海外海鳶烏賊 CPUE 時(shí)空變化

第2航次的漁獲總量和CPUE平均值高于第1航次 (表2)，相較于第1航次的燈光覆網(wǎng)漁船，燈光罩網(wǎng)漁船在作業(yè)水深與作業(yè)方式上有所不同，導(dǎo)致了2次探捕的結(jié)果不同，燈光罩網(wǎng)的作業(yè)水深較燈光覆網(wǎng)的深，因此燈光罩網(wǎng)漁獲物的種類、年齡結(jié)構(gòu)等較燈光覆網(wǎng)更豐富[14]，且季節(jié)性漁獲物也存在差異。

鳶烏賊CPUE高值的集中區(qū)域主要有3個(gè)，分別是 60°E—61°E，16°N；62°E，16°N—18°N；61.50°E—63°E，17°N—17.75°N (圖2)。陳新軍和葉旭昌[15]通過對 2003 年印度洋西北海域 (58°E—65°E，2°N—18°N) 探捕生產(chǎn)情況研究發(fā)現(xiàn)，10 月下旬鳶烏賊中心漁場分布在61°E，15°N—16°N附近海域，這與本研究的結(jié)果相似。錢衛(wèi)國等[16]發(fā)現(xiàn)在2004年冬季調(diào)查期間60°E，14°30'N附近海域鳶烏賊的CPUE最高，該研究采用的數(shù)據(jù)為魷釣漁船數(shù)據(jù)，因此造成與本研究結(jié)果有差異，本研究采用的是燈光圍網(wǎng)漁船數(shù)據(jù)，在燈光集魚效果、作業(yè)深度、漁獲選擇性等方面與魷釣作業(yè)方式的調(diào)查有較大區(qū)別，并且該研究的調(diào)查時(shí)間在2003年左右，與本研究相距時(shí)間較長，而印度洋西北海域的季風(fēng)氣候的情況復(fù)雜，在秋冬季節(jié)為東北季風(fēng)而夏季為西南季風(fēng)[17]，導(dǎo)致每年鳶烏賊和澳洲鮐 (Scomber australasicus)[18]等中上層魚類的漁場分布和產(chǎn)量會產(chǎn)生年際間和季節(jié)性波動。本研究調(diào)查持續(xù)時(shí)間較短且不連續(xù)，存在樣本量不足的缺點(diǎn)，故暫不適宜進(jìn)行季節(jié)性的具體研究，有待后續(xù)補(bǔ)充調(diào)查。

根據(jù)GAM模型結(jié)果，農(nóng)歷日 (月亮夜與月光夜) 對阿拉伯海外海鳶烏賊的CPUE具有極其顯著的影響 (P≤0.001)。影響曲線在15 d后呈現(xiàn)升高趨勢 (圖9-c)，說明在滿月附近的時(shí)段，漁船為了減少作業(yè)成本會減少每天的作業(yè)次數(shù)，而在滿月之后，月光照度開始下降，漁船會進(jìn)入正常生產(chǎn)狀況；月黑夜總作業(yè)次數(shù)是月亮夜的1.38倍，月黑夜的總漁獲量較月亮夜增加了108.21% (表3)，基本與晏磊等[19]、招春旭等[20]、李杰等[21]研究結(jié)果一致，進(jìn)一步表明農(nóng)歷日對鳶烏賊分布的影響顯著。另外，有研究表明，除了月光會導(dǎo)致水上集魚燈的效果降低，還會對水溫、海流等外部環(huán)境條件造成影響[22]；同時(shí)，不同物種的趨光特性也會有所不同，錢衛(wèi)國等[23]研究證明印度洋鳶烏賊適宜光照強(qiáng)度低于智利外海莖柔魚 (Dosidicus gigas)。因此依據(jù)主捕漁獲物的生物學(xué)習(xí)性，選擇不同類型不同功率的集魚燈，適當(dāng)配置水上集魚燈和水下集魚燈的數(shù)量，有利于漁船降低生產(chǎn)成本，提高經(jīng)濟(jì)效益，減少漁獲物中的兼捕物種，保護(hù)該海域的生物資源。

3.2 主要環(huán)境要素對鳶烏賊 CPUE 的影響

上升流是一種海水垂直向上運(yùn)動的現(xiàn)象，上升的水體通常來自數(shù)百米深處，雖然上升流海域的面積僅占海洋總面積的1‰，但漁獲量卻占世界海洋總漁獲量的50%[7,24-25]。上升流是印度洋鳶烏賊漁場形成與變動的內(nèi)在驅(qū)動[26]，表層海水與低溫的深層海水交匯，SST發(fā)生變化，且因?yàn)樯仙鲙由詈５讓拥臓I養(yǎng)鹽來到表層，表層藻類生長旺盛，導(dǎo)致Chla濃度上升，由此形成了阿拉伯海外海漁場。SSTA能反映漁場的水溫變化狀況，與上升流的強(qiáng)弱具有一定的相關(guān)性[27]。因此，本文利用SST、Chla和SSTA作為海洋環(huán)境因子，研究其與鳶烏賊CPUE的關(guān)系是可行的。本研究中SST、Chla濃度、SSTA對阿拉伯海外海鳶烏賊CPUE均存在影響 (P≤0.05)。

GAM模型表明，阿拉伯海外海鳶烏賊的最適SST介于26.5～27.5 ℃。林東明和陳新軍[28]根據(jù)在印度洋西北部海域2004年9月—2005年4月的生產(chǎn)調(diào)查資料，發(fā)現(xiàn)印度洋鳶烏賊的漁場主要分布在 60°E—61.50°E，15°N—16.50°N，作業(yè)漁場的適宜SST為25～29 ℃，該結(jié)果與本研究范圍基本一致；另外周艷波等[29]認(rèn)為SST、Chla濃度和海面高度 (Sea surface height, SSH) 與南海外海鳶烏賊分布有關(guān)，其中最適宜SST為25～31 ℃；徐紅云[30]研究發(fā)現(xiàn)1—12月南海外海鳶烏賊漁場的SST介于25.6～29.3 ℃并呈現(xiàn)由南到北遞減的趨勢；李杰等[21]研究表明29.00～29.49 ℃適合南海中南部海域鳶烏賊的捕獲，以上相關(guān)南海海域的研究結(jié)果有所差異，是因?yàn)檎{(diào)查海域的具體位置及調(diào)查時(shí)間不同，并且研究阿拉伯海外海鳶烏賊漁汛期為10月—翌年4月，調(diào)查期間SST介于27～29 ℃，而南海外海鳶烏賊全年均有作業(yè)，漁汛期間的SST介于25～31 ℃[29]，這導(dǎo)致南海外海鳶烏賊的SST范圍較阿拉伯海外海的范圍寬。

Chla濃度決定了浮游動物如橈足類、糠蝦類等的分布，阿拉伯海外海海域豐富的浮游動物資源是鳶烏賊的重要餌料[16]，因此Chla濃度能在一定程度上影響鳶烏賊的分布。本研究中，阿拉伯海外海鳶烏賊的適宜 Chla質(zhì)量濃度介于 0.30～0.50 mg·m?3。這與楊曉明等[2]以及余為和陳新軍[31]的結(jié)論有一定差異，筆者推測是上述研究中所用漁具漁法為釣鉤，因此漁獲組成和個(gè)體差異等與燈光圍網(wǎng)的漁獲會存在一定差異。此外，由于阿拉伯海外海的海洋環(huán)境復(fù)雜，同時(shí)受到索馬里海流和赤道海流的影響[32]，Chla濃度范圍變化較大，2個(gè)航次的調(diào)查研究時(shí)間較短，而沒有覆蓋完整的鳶烏賊漁汛期。同時(shí)，通過對比南海外海鳶烏賊的相關(guān)研究，發(fā)現(xiàn)阿拉伯海外海鳶烏賊的最適Chla濃度范圍高于南海外海[7,10,29-30]，也表明不同海域的相同物種會因?yàn)樯巢煌a(chǎn)生一定差異。

SSTA，又稱為SST異常，常用于捕捉SST在的異常變化特征，相關(guān)研究發(fā)現(xiàn)，當(dāng)阿拉伯海外海漁場附近有上升流發(fā)生，會出現(xiàn)較大的SST梯度和Chla梯度變化，導(dǎo)致魚群匯集從而形成漁場，而當(dāng)上升流消失時(shí)，漁場消失[2]。環(huán)境因子中，SSTA對鳶烏賊CPUE的影響極其顯著 (P≤0.001)，且關(guān)聯(lián)度大于SST和 Chla。當(dāng)SSTA小于0 ℃或大于0.4 ℃時(shí)，整體趨勢與鳶烏賊CPUE呈現(xiàn)負(fù)相關(guān)關(guān)系，這與鳶烏賊作為頭足類動物對SST的異常變動較為敏感這一生物學(xué)特性相一致，異常的溫度變動不僅會影響頭足類自身的生長[33]，還可能對餌料生物的豐度產(chǎn)生影響[16]。后期研究計(jì)劃利用SSH分布和海洋環(huán)流等物理要素推測漁場的上升流情況，結(jié)合漁業(yè)生產(chǎn)調(diào)查資料以探究阿拉伯海外海鳶烏賊漁場的內(nèi)在驅(qū)動因子，有利于今后合理地開發(fā)和利用阿拉伯海外海漁業(yè)資源。

4 結(jié)論

本文采用GAM模型分析了時(shí)空及海洋環(huán)境因子對阿拉伯海外海鳶烏賊CPUE的影響，初步探究了鳶烏賊CPUE變化與時(shí)空因子和SST等海洋環(huán)境因子的關(guān)系，并得出了相關(guān)結(jié)論，但仍有所不足：1) 本研究中雖然SST和Chla對CPUE有顯著影響，但模型解釋率并不高，這可能與經(jīng)緯度在模型中偏差解釋率較大有關(guān)[27]；2) 影響阿拉伯海外海鳶烏賊CPUE的環(huán)境因子可能不僅有SST、Chla和SSTA，也有其他表征海流或物理海洋過程的因子未被考慮；3) 本研究采用的是2個(gè)不同航次的探捕數(shù)據(jù)，雖然可以保證數(shù)據(jù)的真實(shí)性和準(zhǔn)確性，但調(diào)查的經(jīng)緯度范圍未能覆蓋整個(gè)西北印度洋漁場區(qū)域，調(diào)查時(shí)間也未能包含整個(gè)漁汛期，存在數(shù)據(jù)覆蓋面窄的缺點(diǎn)。綜上，在今后的研究中，建議將多年份及大面積的調(diào)查數(shù)據(jù)進(jìn)行整合，包含整個(gè)西北印度洋的漁汛期和漁場，結(jié)合海洋環(huán)境數(shù)據(jù)開展深入的調(diào)查與重點(diǎn)分析，為今后該區(qū)域漁場系統(tǒng)的分析提供科學(xué)依據(jù)。