趙國慶 羅俊榮 唐峰華 樊 偉 宋學(xué)鋒 楊 超 張 衡

基于漁業(yè)數(shù)據(jù)的南極磷蝦48漁區(qū)漁場時空分布*

趙國慶1,2羅俊榮3唐峰華2樊 偉2宋學(xué)鋒2楊 超1,2張 衡2①

(1. 上海海洋大學(xué)海洋科學(xué)學(xué)院 上海 201306；2. 中國水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所 上海 200090；3. 福建正冠漁業(yè)有限公司 福建 福州 350500)

根據(jù)我國2010—2019年南極磷蝦()捕撈漁船的生產(chǎn)資料，分析了南極海域48漁區(qū)南極磷蝦漁場的分布特點，采用重心遷移軌跡模型和標(biāo)準(zhǔn)差橢圓(SDE)模型探討了南極磷蝦的漁場變動特征和規(guī)律。結(jié)果顯示，南極磷蝦捕撈量主要集中在48.1亞區(qū)，占比為70.30%，48.2亞區(qū)和48.3亞區(qū)的產(chǎn)量相差很小，占比分別為14.28%和15.42%；年間單位捕撈努力量漁獲量(CPUE)曲線上升，最小值為2012年，最大值為2019年；月間CPUE先增后降，最小值為1月，最大值為6月。48.1亞區(qū)的年間和月間漁場重心均往西南方向移動；48.2亞區(qū)年間的漁場重心東移，但移動范圍較小，月間規(guī)律不強；48.3亞區(qū)年間漁場重心南移，月間漁場重心向西北移動。經(jīng)SDE分析可知，48.1亞區(qū)漁場分布范圍最廣、離散程度最大，48.3亞區(qū)漁場方向性最強、向心力最明顯。48.1亞區(qū)漁場重心主要分布于布蘭斯菲爾德海峽，48.2亞區(qū)漁場重心分布于南奧克尼群島東側(cè)，48.3亞區(qū)漁場重心分布于南喬治亞群島東北側(cè)。聚類結(jié)果表明，48.1亞區(qū)年間漁場重心均較為集中，48.2和48.3亞區(qū)除2017年外，其他年間漁場重心較為集中。

南極磷蝦；漁場重心；標(biāo)準(zhǔn)差橢圓分析；聚類分析；時空分布；48漁區(qū)

南極磷蝦通常指南極大磷蝦(,以下稱磷蝦)，屬于甲殼類浮游動物，廣泛分布于南極大陸周圍，既是浮游植物的主要捕食者，又是魚類、頭足類、企鵝、海豹、鯨魚等的主要餌料，在食物網(wǎng)中起著承上啟下的作用，是整個南大洋生態(tài)系統(tǒng)物質(zhì)和能量流動的關(guān)鍵物種(Siegel, 2005; Atkinson, 2012; Stowasser, 2012)。在地球已知的多細(xì)胞動物中，磷蝦生物量最高，據(jù)估算，其資源量可達3.0×108～5.0×108t (Atkinson, 2009)。磷蝦脂類物質(zhì)含量很高，富含長鏈多不飽和脂肪酸(Ericson, 2018)，被稱為全球最大的動物蛋白質(zhì)庫，另外，還富含豐富的礦物元素和胡蘿卜素等(Wang, 2011)。因此，磷蝦具有十分巨大的開發(fā)價值。磷蝦捕撈業(yè)始于20世紀(jì)60年代早期，中國于1984年首次開展南極科學(xué)考察時即將磷蝦列為重點考察目標(biāo)，而正式商業(yè)性探捕開發(fā)始于2009—2010年度(國家海洋局極地專項辦公室, 2016)。目前，我國磷蝦捕撈技術(shù)和捕撈產(chǎn)量都取得了重大進展，由2009—2010年度4000 t上升到2019年50 000多t (CCAMLR, 2020)。

掌握磷蝦漁場時空分布對磷蝦商業(yè)性捕撈具有重大意義。目前，對磷蝦分布特征的研究，主要通過大型漁船的生產(chǎn)資料(Krafft, 2015)、科學(xué)調(diào)查資料(Marrari, 2008)、暫養(yǎng)實驗(Kawaguchi, 2010a)或遠(yuǎn)程監(jiān)測(Saenz, 2020)等進行，相關(guān)研究多結(jié)合環(huán)境數(shù)據(jù)(Atkinson, 2019)或借助聲學(xué)探測(Seok, 2018; 張吉昌等, 2012)來進行，側(cè)重于某一時空尺度下磷蝦的分布特征或已知因子如何影響磷蝦分布，并探討磷蝦分布是否存在統(tǒng)計意義上的聚集或離散性質(zhì)。Krafft等(2015)通過商業(yè)捕撈漁船對磷蝦豐度和垂直結(jié)構(gòu)等種群動態(tài)進行研究后認(rèn)為，使用商業(yè)捕撈漁船作為研究平臺，對了解漁業(yè)動態(tài)及量化生物資源量水平具有很大的優(yōu)勢。Patricio等(2020)通過2011—2016年磷蝦捕撈數(shù)據(jù)總結(jié)了智利磷蝦漁業(yè)現(xiàn)狀，展示了磷蝦捕撈深度、單位捕撈努力量漁獲量(CPUE)、網(wǎng)次地理分布等。張瑛瑛等(2020)利用2010—2015年1—6月的磷蝦捕撈數(shù)據(jù)研究了漁場重心總體上的軌跡變化，并探討了CPUE與海表溫度和葉綠素(Chl-)的關(guān)系。但是，目前基于長時間序列的磷蝦漁業(yè)數(shù)據(jù)對磷蝦漁場的變動時間、變動方向和變動格局的研究鮮有報道。因此，本研究依據(jù)我國2010—2019年磷蝦捕撈漁船的生產(chǎn)資料，探討了48區(qū)磷蝦漁場重心軌跡的遷移和分布格局情況，以期為磷蝦合理開發(fā)利用提供借鑒。

1 材料與方法

1.1 數(shù)據(jù)來源及捕撈海域

漁業(yè)數(shù)據(jù)來源于我國磷蝦捕撈漁船(中國水產(chǎn)總公司、上海開創(chuàng)遠(yuǎn)洋漁業(yè)公司、遼寧遠(yuǎn)洋漁業(yè)公司等10艘漁船)的漁撈記錄資料(圖1)，包括作業(yè)日期(年、月、日)、作業(yè)位置(經(jīng)度和緯度)、漁獲產(chǎn)量、拖網(wǎng)次數(shù)、拖網(wǎng)時間、拖速、網(wǎng)具類型、拖網(wǎng)作業(yè)時的網(wǎng)口高度和網(wǎng)口水平擴張。時間跨度為2010—2019年，主要作業(yè)月份為每年的1—9月。作業(yè)區(qū)域主要分布在南極海域48.1區(qū)、48.2區(qū)和48.3區(qū)，空間分辨率為1°×1°(圖2)。

圖1 南極磷蝦拖網(wǎng)漁船數(shù)量和拖網(wǎng)次數(shù)統(tǒng)計

圖2 調(diào)查海域示意圖

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

1.2.1 漁獲數(shù)據(jù)處理 對漁獲數(shù)據(jù)進行標(biāo)準(zhǔn)化處理，采用標(biāo)準(zhǔn)化后的CPUE來表示磷蝦漁場資源豐度，單位為t/h，計算公式為：

1.2.3 標(biāo)準(zhǔn)差橢圓分析 標(biāo)準(zhǔn)差橢圓(standard deviational ellipse, SDE)是一種定量描述研究對象空間分布整體特征及時空演變過程的空間格局統(tǒng)計分析方法，能夠精確表達地理要素空間分布的整體特征(Lefever, 1926)，本研究采用產(chǎn)量作為權(quán)重字段，用其來揭示磷蝦漁場的變動方向和分布格局。SDE相關(guān)參數(shù)計算公式如下：

1.2.4 漁場重心聚類分析 利用歐式距離來計算年產(chǎn)量重心的距離，比較漁場重心的變化情況，公式如下(Carlo, 2021)：

式中，D為年和年產(chǎn)量重心之間的距離，、分別為第年的經(jīng)度和緯度，X、Y分別為第年的經(jīng)度和緯度。根據(jù)歐式距離，對2010—2019年各年的產(chǎn)量重心按照最小距離法進行聚類，分析比較年份和月份之間的差異性。

2 結(jié)果與分析

2.1 磷蝦產(chǎn)量和CPUE時空變化

將捕撈區(qū)域分為1°×1°若干單位小漁區(qū)，統(tǒng)計2010—2019年每單位漁區(qū)的磷蝦總產(chǎn)量，產(chǎn)量分布見圖3A，CPUE分布見圖3B。2010—2019年的產(chǎn)量集中于48.1亞區(qū)，48.2、48.3亞區(qū)的產(chǎn)量占比和有產(chǎn)量的單位小漁區(qū)分布數(shù)量相差較小，且產(chǎn)量多集中在各亞區(qū)的島嶼周圍。48.1、48.2和48.3亞區(qū)的產(chǎn)量分別占總產(chǎn)量的70.30%、14.28%和15.42%。各亞區(qū)CPUE的變化與產(chǎn)量的變化一致，各區(qū)域中平均CPUE值以48.1亞區(qū)最高，其次為48.2亞區(qū)，48.3亞區(qū)最低，分別為19.95、12.31和11.69 t/h。

圖3 2010—2019年磷蝦產(chǎn)量空間變化(A)和CPUE空間變化(B)

年間捕撈力量的投入存在較大差異，年間產(chǎn)量差距較大，因此，本研究僅對CPUE進行統(tǒng)計分析。從年間變化來看，CPUE呈上升趨勢，在2019年達到最大值，最小值發(fā)生在2012年。網(wǎng)次占比的變化情況基本與CPUE的變動趨勢相似，最小值和最大值分別發(fā)生在2010年和2016年(圖4A)。從月間變化來看，1—9月CPUE和網(wǎng)次占比均呈先增后降的趨勢，且均在9月達到最低值(圖4B)。

2.2 磷蝦產(chǎn)量重心的時空變化

從年間變化來看，48.1亞區(qū)的產(chǎn)量重心呈現(xiàn)西南–東南–西南的變動趨勢，產(chǎn)量重心總體上向西南方向移動(圖5A)；48.2亞區(qū)的產(chǎn)量重心在2010—2016年基本上在東西方向上來回移動，在2017年往南移動，2018年和2019年又基本移動到2013年重心的經(jīng)度上，但在緯度上東移，產(chǎn)量重心總體向東移動(圖5B)；48.3亞區(qū)的產(chǎn)量重心除2017年外總體上在南北方向上移動，且整體上往南移動(圖5C)。

圖4 2010—2019年磷蝦平均CPUE的年間(A)變化和月間(B)變化

圖5 48區(qū)磷蝦歷年產(chǎn)量重心年間(A、B和C)和月間(D、E和F)變化

A、B和C分別表示48.1、48.2和48.3亞區(qū)漁場重心的年間變動，D、E和F分別表示48.1、48.2和48.3亞區(qū)漁場重心的月間變動。

A, B and C represent the annual variation of fishing ground gravity center in sub region 48.1, 48.2 and 48.3 respectively; D, E and F represent the monthly variation of fishing ground gravity center in sub region 48.1, 48.2 and 48.3 respectively.

從月間變化來看，48.1亞區(qū)產(chǎn)量重心總體上逐漸向西南方向移動(圖5D)；48.2亞區(qū)的產(chǎn)量重心隨月份無明顯的變動規(guī)律，產(chǎn)量重心集中在46.20°W～46.75°W、60.90°S～61.55°S范圍內(nèi)，變動范圍較小(圖5E)；48.3亞區(qū)僅在6—9月有產(chǎn)量存在，產(chǎn)量重心總體上逐漸向西北移動(圖5F)。

2.3 磷蝦漁場變動方向和分布格局

為進一步揭示48區(qū)磷蝦漁場的變動方向和分布格局，采用SDE法對磷蝦的重心遷移軌跡進行了互補分析。如圖6A所示，對2010—2019年10年的數(shù)據(jù)進行匯總分析可得，從分布方向上看，磷蝦資源在48.1亞區(qū)基本上沿布蘭斯菲爾德海峽分布，在48.2亞區(qū)則分布在南奧克尼群島偏東方向，在48.3亞區(qū)分布在南喬治亞群島東北側(cè)。另外，通過橢圓的扁率、面積、長軸和短軸可知，磷蝦資源在48.1亞區(qū)的分布范圍最廣，離散程度最大；在48.3亞區(qū)的方向性最強，向心力最明顯。

圖6B、6C、6D和表1分別展示了48.1、48.2和48.3亞區(qū)2010—2019年每年磷蝦漁場的變動情況、產(chǎn)量分布SDE形狀的相關(guān)參數(shù)。48.1亞區(qū)磷蝦的漁場重心逐漸從西北往東南遷移到布蘭斯菲爾德海峽內(nèi)，年間擴張海域面積差距較大，在2017年達到最大值。除2011年外，漁場分布方位角均在54.05°～87.46°之間變動，表明磷蝦漁場分布呈東南–西北格局，2015、2017年的扁率較大，表明2015年和2017年的方向性較強。

圖6 2010—2019年磷蝦產(chǎn)量分布標(biāo)準(zhǔn)差橢圓

表1 2010—2019年磷蝦產(chǎn)量分布標(biāo)準(zhǔn)差橢圓形狀參數(shù)

Tab.1 Variation of parameters of standard deviational ellipse of production distribution of E. superba from 2010 to 2019

注：作圖時使用經(jīng)緯度作為度量單位，因此，短軸和長軸只用來提供距離比例和計算扁率。

Note: Latitude and longitude are used as the unit of measurement in drawing, so the minor axis and major axis only provide distance proportion and calculation of flatness.

48.2亞區(qū)除2011年和2018年外，漁場重心及擴張海域均在南奧克尼群島偏東方向，且結(jié)合圖6A可以得出，48.2亞區(qū)的磷蝦漁場擴張海域和方向基本受2011年和2017年支配，表明這2年的產(chǎn)量較其他年份占比大。漁場分布方位角的變動范圍為2.76°～ 103.35°，年間變動較大，不具有規(guī)律，2012、2015年的扁率較大，方向性較強。

48.3亞區(qū)除2017年外，其他年份漁場重心及擴張海域均在南喬治亞群島東北方向，且較為集中，在53°30?S～54°30?S、35°20?W～36°20?W之間(圖6E)，2017年的擴張海域最大，且漁場重心與其他年份均有較大差距，結(jié)合圖6A可知，漁場重心和擴張海域基本不受2017年影響，因此，2017年的產(chǎn)量占比較小。漁場分布方位角在90.27°～155.03°之間變動，表明漁場基本分布在東西方向，2017年的扁率最大，方向性最強。

2.4 產(chǎn)量重心聚類分析

2010—2019年48.1區(qū)不同年份產(chǎn)量重心分布差異較小，由歐式距離可得，年間產(chǎn)量重心的空間距離范圍在0.122 (2014年/2016年)～2.223 (2010年/2017年)之間，且歐氏距離超過2的僅有2010年/2017年和2010年/2019年，其他年份之間均小于2(表2)。經(jīng)聚類分析可得，若以空間距離2為閾值，則可將產(chǎn)量重心分為4類，即2010年為一組，2011年和2012年為一組，2013年、2014年、2015年、2016年和2018年為一組，2017年和2019年為一組(圖7A)。48.2亞區(qū)2017年產(chǎn)量重心分布與其他年份間差異較大，2017年與其他年份之間的歐式距離均>2，其他年份之間的歐式距離均<1，表明除了2017外，其他年份漁場重心較為集中(表3)。經(jīng)聚類分析可得，若以空間距離2為閾值，可以將年產(chǎn)量重心分為2類，即2017年為一組，其他年份為一組(圖7B)。48.3亞區(qū)2017年產(chǎn)量重心分布與其他年份差異較大，均>5，其他年份之間重心的歐式距離均<1，表明除2017年外，漁場重心較為集中(表4)。經(jīng)聚類分析可得，以空間距離2為閾值，可將年產(chǎn)量重心分為2類，即2017年為一組，其他年份為一組(圖7C)。

表2 48.1漁區(qū)磷蝦各年產(chǎn)量重心的空間距離

Tab.2 Distribution of annual catch gravity position of E. superba in the 48.1 fishing area

圖7 各漁區(qū)磷蝦年產(chǎn)量聚類結(jié)果

表3 48.2漁區(qū)磷蝦各年產(chǎn)量重心的空間距離

Tab.3 Distribution of annual catch gravity position of E. superba in the 48.2 fishing area

表4 48.3漁區(qū)磷蝦各年產(chǎn)量重心的空間距離

Tab.4 Distribution of annual catch gravity position of E. superba in the 48.3 fishing area

3 討論

3.1 磷蝦漁場時空變化特征

自20世紀(jì)70年代磷蝦漁業(yè)商業(yè)性開發(fā)以來，發(fā)展迅速，磷蝦漁場的時空分布和資源量變動研究等越來越受重視。然而，有關(guān)磷蝦漁場的年間和季節(jié)性分布、資源量變化的研究結(jié)果多有分歧(Siegel, 2005)。磷蝦漁業(yè)活動一直處在變化中，漁場和漁季也隨之發(fā)生變化。20世紀(jì)70年代，磷蝦漁業(yè)遍布整個南極大陸周圍海域，80年代末，印度洋成為磷蝦捕撈的主要漁場，90年代開始逐漸向南大西洋轉(zhuǎn)移，自1996年以來，磷蝦漁業(yè)幾乎完全集中在南設(shè)得蘭群島、南喬治亞島和南奧克尼群島(圖8) (Krafft, 2015; Kawaguchi, 2007; CCAMLER, 2020)。自1973年以來，磷蝦漁場范圍逐漸縮小，捕撈量表現(xiàn)出較強的波動性。總體來看，20世紀(jì)90年代之前磷蝦的年均總產(chǎn)量要遠(yuǎn)大于90年代以后的年均產(chǎn)量，但90年代以后，其產(chǎn)量呈現(xiàn)逐年增加的趨勢(趙國慶, 2018)。磷蝦捕撈量的季節(jié)性變化也存在較大波動，20世紀(jì)80年代磷蝦捕撈產(chǎn)量主要集中在夏季，但自90年代以后，秋冬季磷蝦產(chǎn)量比例逐漸增大并占據(jù)優(yōu)勢(Kawaguchi, 2006)。趙國慶(2018)通過對以往相關(guān)文獻和資料總結(jié)分析后發(fā)現(xiàn)，磷蝦全年均可生產(chǎn)。目前，主要捕撈區(qū)域位于48區(qū)，其中，48.1區(qū)作業(yè)時間最長，1—6月產(chǎn)量較高；48.2區(qū)主要作業(yè)時間為3—8月；48.3區(qū)主要作業(yè)時間為5—10月，其中，5—8月產(chǎn)量較高。李顯森等(2015)對2012/2013漁季磷蝦的捕撈量進行統(tǒng)計分析后認(rèn)為，各漁區(qū)中CPUE均值以48.1區(qū)最高，48.3區(qū)最低；月均CPUE值在1—6月較為穩(wěn)定，7—9月逐漸下降。張瑛瑛等(2020)對中國2010—2015年磷蝦捕撈數(shù)據(jù)分析后發(fā)現(xiàn)，磷蝦漁場年際重心由南設(shè)德蘭群島周邊海域向東北方的南奧克尼群島、南喬治亞群島周邊海域推移，且大多數(shù)捕撈作業(yè)頻次分布于南設(shè)德蘭群島夏季時段(1—3月)。Francisco等(2018)通過對各個時空尺度上的捕撈漁船活動情況和磷蝦產(chǎn)量分析后發(fā)現(xiàn)，歷史上的漁場和漁獲量的空間分布已經(jīng)發(fā)生了很大變化，并確定了布蘭斯菲爾德海峽為主要磷蝦漁場。

據(jù)南極海洋生物資源養(yǎng)護委員會(Commission for the Conservation of Antarctic Marine Living Resources, CCAMLR)統(tǒng)計，自1996年以來，磷蝦的捕撈作業(yè)幾乎全部集中在48.1、48.2和48.3亞區(qū)(CCAMLR, 2020; CCAMLR Secretariat, 2020)，因此，本文對該海域的磷蝦進行時空分布的研究具有較大的現(xiàn)實意義。磷蝦廣泛分布于南極大陸周圍各個海域，棲息地范圍為19×106km2(Atkinson, 2009)，密集區(qū)基本上出現(xiàn)在大陸架邊緣、冰架邊緣、海峽以及島嶼周圍(Siegel, 2005; Atkinson, 2009)。Patricio等(2020)對智利2011—2016年磷蝦捕撈數(shù)據(jù)分析后發(fā)現(xiàn)，48.1區(qū)作業(yè)范圍逐年擴大；48.區(qū)的捕撈活動主要發(fā)生在南奧克尼群島西北部，捕撈頻次逐年減少；48.3區(qū)的捕撈活動主要集中在南喬治亞島東北部，捕撈范圍逐年縮小。這與本文所得出的結(jié)論相吻合，即磷蝦漁場重心在48.1亞區(qū)基本上沿布蘭斯菲爾德海峽分布，在48.2亞區(qū)則分布在南奧克尼群島東側(cè)，在48.3亞區(qū)分布在南喬治亞群島東北側(cè)，表明中國與智利的磷蝦捕撈區(qū)域有一定的重合。在此基礎(chǔ)上，本研究運用重心遷移軌跡模型和SDE模型，從更深層次挖掘磷蝦漁場重心在各個亞區(qū)的移動軌跡和分布格局。2010—2019年磷蝦產(chǎn)量和捕撈網(wǎng)次主要集中在48.1亞區(qū)，且斯菲爾德海峽的磷蝦產(chǎn)量大于48.2亞區(qū)和48.3亞區(qū)，基本與CCAMLR的統(tǒng)計結(jié)果相吻合(CCAMLR, 2020; CCAMLR Secretariat, 2020)。從本研究還可以看出，48.1亞區(qū)磷蝦漁場的分布范圍要遠(yuǎn)遠(yuǎn)大于48.2亞區(qū)和48.3亞區(qū)。

圖8 南極海洋生物資源養(yǎng)護委員會統(tǒng)計的磷蝦漁獲量空間分布(CCAMLR Secretariat, 2020)

3.2 磷蝦漁場變動的影響因素

磷蝦的商業(yè)性捕撈已近50年，磷蝦漁業(yè)多次被強調(diào)為世界上最后未被開發(fā)的漁業(yè)之一(Garcia, 2011)，具有很大的發(fā)展?jié)摿?，幾十年來，磷蝦漁場的主要海域、漁季都發(fā)生了巨大變化。

磷蝦漁業(yè)往往跟隨磷蝦群進行捕撈作業(yè)，因此磷蝦漁場的選擇與其集群行為具有很大關(guān)系。集群是磷蝦的主要生態(tài)行為，在它們的整個生命史中都發(fā)揮著重要作用，它們通過集群在某些區(qū)域高度集中，但是在某些區(qū)域幾乎完全沒有分布(Murphy, 1988)，有時磷蝦集合成群時面積可達100萬km，內(nèi)含200萬t磷蝦(Siegel, 2005)。磷蝦的集群和分布是多種因素互相作用的結(jié)果，很少有研究能夠證明單一因素諸如溫度、鹽度、溶解氧濃度、溶解有機質(zhì)、海水理化性質(zhì)、洋流、鋒面系統(tǒng)和經(jīng)緯度等與磷蝦密度或分布具有明顯的可預(yù)測關(guān)系(Siegel, 2016)。磷蝦往往因為洋流和水深的共同影響集中分布于陸架邊緣和峽谷，可能會密集聚集形成熱點區(qū)域(Bernard, 2017)。本文與相關(guān)研究(Patricio, 2020)發(fā)現(xiàn)，近年來，磷蝦捕撈多集中于島嶼周圍(南奧克尼群島和南舍得蘭群島)和海峽(布蘭斯菲爾德海峽)內(nèi)。據(jù)研究，磷蝦生物量的聚集也依賴于潮汐(Bernard, 2017)、風(fēng)向(Krafft, 2015)的情況，在晝夜潮汐期間，安弗斯島附近水域的生物量顯著更高，當(dāng)西風(fēng)為主和局地混合潮處于日變化狀態(tài)時，大密度聚集更頻繁發(fā)生。但不得不提的是，磷蝦的集群分布難以預(yù)測，到目前為止，仍無法準(zhǔn)確掌握磷蝦的集群機制(Krafft, 2015)。

磷蝦漁業(yè)是一種利益驅(qū)動的經(jīng)濟性活動，漁船通常以利潤最大化為目標(biāo)(Krafft, 2015)，磷蝦漁場的變遷、漁獲量的增減除了受各種物理環(huán)境因素的影響外，還與磷蝦的商業(yè)用途、市場的需求、捕撈漁船的增減和磷蝦管理政策等有很大關(guān)系(Nicol, 2012)。影響漁船收入的因素有很多，但磷蝦集群的密度具有預(yù)測性。Krafft等(2015)指出，聲學(xué)得到磷蝦集群密度最高的深度與實際捕撈深度一致，并表示漁船更傾向于在垂直方向上瞄準(zhǔn)高密度磷蝦群。但有時磷蝦漁場的變動可能并不完全由磷蝦集群的密度決定。Hewitt等(2004)分析了南設(shè)德蘭群島周圍海域磷蝦捕撈隨時間的變化模式后指出，該海域磷蝦密度較高的區(qū)域主要有象島東端、象島和喬治國王島嶼之間和希里夫角附近，但前二者均不是磷蝦的主要漁場。磷蝦捕撈量的時空波動與船隊組成也有一定關(guān)聯(lián)。前蘇聯(lián)自1972年開始在南極水域進行商業(yè)捕撈以來一直占據(jù)主導(dǎo)地位，年捕撈量超過20萬t，主要的捕撈地為南極洲東部海域和南大西洋(Nicol, 2011)。20世紀(jì)90年代初期前蘇聯(lián)解體后，磷蝦產(chǎn)量大幅下降，主要來自于如本和韓國，后來隨著挪威、中國、韓國等國家相繼加入磷蝦的捕撈，其產(chǎn)量有較大提升，挪威逐漸占據(jù)主導(dǎo)地位，磷蝦漁場逐漸向南大西洋遷移(陳雪忠等, 2009; 趙憲勇等, 2016; CCAMLR Secretariat, 2020)。

因為磷蝦漁船的捕撈活動只發(fā)生在無冰區(qū)，海冰是影響捕撈漁船作業(yè)的主要因素，捕撈熱點區(qū)域位置與海冰覆蓋度關(guān)系密切，這也導(dǎo)致磷蝦漁場季節(jié)性明顯(Kawaguchi, 2010b)。另外，海冰為磷蝦提供了優(yōu)良的棲息環(huán)境，為越冬的磷蝦群提供了優(yōu)良的索餌環(huán)境，也減少了被捕撈和獵食的風(fēng)險(Brierley, 2002)。據(jù)研究，48漁區(qū)磷蝦夏季的磷蝦豐度與上一季度冬季海冰的面積呈正比(Hewitt, 2004)，而磷蝦漁業(yè)中，CPUE與海冰總面積年間變化呈顯著負(fù)相關(guān)關(guān)系(戴立峰等, 2012)。近年來，全球氣候變暖導(dǎo)致海冰形成延遲、無冰區(qū)域擴大，進而使磷蝦捕撈的時間跨度變得越來越長，曾只在夏季的捕撈活動現(xiàn)已持續(xù)到秋季(Kawaguchi, 2006)。而冰川融化對磷蝦的生存具有不容忽視的威脅，據(jù)研究，氣候引起的冰川融水導(dǎo)致海洋中懸浮顆粒增加，大量的巖性顆粒影響了磷蝦的攝食、吸收能力和交配行為，這可能是南極洲波特灣磷蝦大量死亡的原因(Fuentes, 2016)。隨著全球氣候變暖，在過去的90年里，西南大西洋磷蝦的分布逐漸向南收縮，在分布的北部極限區(qū)域密度急劇下降，密度重心越來越靠近南極冰架(Atkinson, 2019)，這可能也是導(dǎo)致48.1亞區(qū)磷蝦漁場重心逐漸向西南方向移動的主要原因。據(jù)研究，在未來幾十年里南大洋將會持續(xù)變暖(Michael, 2021)，而這可能進一步對磷蝦的棲息地、分布等產(chǎn)生重要影響。Veytia等(2020)模擬了磷蝦棲息地質(zhì)量的季節(jié)性變化，指出未來80年里海洋溫度和初級生產(chǎn)力的改變可能會導(dǎo)致磷蝦棲息地萎縮，漁場不斷向高緯度地區(qū)轉(zhuǎn)移。

隨著磷蝦漁業(yè)的發(fā)展，磷蝦漁場的掌控能力和捕撈技術(shù)都有了長足的發(fā)展。自2007年以來，挪威將連續(xù)泵吸技術(shù)應(yīng)用到磷蝦捕撈，磷蝦漁獲產(chǎn)量持續(xù)增加，已達到季節(jié)性配額，導(dǎo)致漁場提前關(guān)閉(Francisco, 2018)。目前，磷蝦漁業(yè)多集中在西南大西洋(48區(qū))，CCAMLR對48.1、48.2和48.3亞區(qū)的捕撈出發(fā)限額分別為15.5萬t、27.9萬t和27.9萬t (Nicol, 2012)。48.1亞區(qū)作為捕撈熱點區(qū)域，已多次達到觸發(fā)水平，對磷蝦捕撈時長產(chǎn)生了一定影響(CCAMLR Secretariat, 2020)。目前，CCAMLR以亞區(qū)為單位進行捕撈限額的制定，而并未考慮磷蝦的聚集特征。南極半島周邊水域，特別是布蘭斯菲爾德海峽，作為磷蝦重要的產(chǎn)卵場、索餌場和越冬場，大量磷蝦聚集在此，捕撈區(qū)域和捕撈季節(jié)與企鵝和鯨魚的捕食活動有較大的重合(Weinstein, 2017)。Francisco等(2018)指出，應(yīng)在較小的空間單元中重新分配預(yù)防性捕撈限制，以盡量減少對依賴磷蝦的捕食者的不良影響。大約20年前，CCAMLR探索過這種依據(jù)小單元進行管理的方法，目前仍有待實施(Hewitt, 2004)。在局部漁獲產(chǎn)量不斷提升的前提下，這種管理方法很可能會得到采用，而這將進一步對磷蝦漁場的時空分布產(chǎn)生影響。

磷蝦漁場的選擇需要多方面考慮，在進行安全捕撈的前提下，既要嚴(yán)格遵守CCAMLR的管理措施，又需要實現(xiàn)利益最大化，特別是目前氣候變化可能已經(jīng)對南極產(chǎn)生了不可逆的影響，這就要求我們要實時掌握磷蝦資源和其生存環(huán)境的相關(guān)信息，并做出迅捷反應(yīng)，掌握更先進的泵吸捕撈技術(shù)，提升漁場的掌控能力，在世界磷蝦漁業(yè)中占據(jù)有利地位。

ATKINSON A, SIEGEL V, PAKHOMOV E A,. A re-appraisal of the total biomass and annual production of Antarctic krill. Deep-Sea Research Part Ⅰ, 2009, 56(5): 727–740

ATKINSON A, NICOL S, KAWAGUCHI S,. Fittinginto southern ocean foodweb models: A review of data sources and their limitations. CCAMLR Science, 2012, 19: 219–245

ATKINSON A, HILL S L, PAKHOMOV E A,. Krill () distribution contracts southward during rapid regional warming. Nature Climate Change, 2019, 142–149

BERNARD K S, CIMINO M, FRASER W,. Factors that affect the nearshore aggregations of Antarctic krill in a biological hotspot. Deep Sea Research Part Ⅰ: Oceanographic Research Papers, 2017, 126: 139–147

BRIERLEY A S, FERNANDES P G, BRANDON M A,. Antarctic krill under sea ice: Elevated abundance in a narrow band just south of ice edge. Science, 2002, 295(5561): 1890–1890

CARLO R. From the Euclidean distance to compositional dissimilarity: What is gained and what is lost. Acta Oecologica, 2021, 111: 1–3

CCAMLR. CCAMLR statistical bulletin volume 32. Tasmania, Australia: Commission for the conservation of Antarctic marine living resources, 2020

CCAMLR Secretariat. Fishery report 2020:in Area 48. Tasmania, Australia: Commission for the conservation of Antarctic marine living resources, 2020

CHEN M F, LIU F, ZHAO Y L,. The expansion and gradient evolution of impervious surface within the Guangdong-Hong Kong-Macao Greater Bay Area from 1987 to 2017. Journal of Geo-Information Science, 2021, 23(3): 443–455 [陳明發(fā), 劉帆, 趙耀龍, 等. 1987—2017年粵港澳大灣區(qū)城市群的不透水地表擴張與梯度演化. 地球信息科學(xué)學(xué)報, 2021, 23(3): 443–455]

CHEN X Z, XU Z L, HUANG H L,. Development strategy on Antartic krill resource utilization in China. Journal of Fishery Sciences of China, 2009, 16(3): 451–458 [陳雪忠, 徐兆禮, 黃洪亮, 等. 南極磷蝦資源利用現(xiàn)狀與中國的開發(fā)策略分析. 中國水產(chǎn)科學(xué), 2009, 16(3): 451–458]

DAI L F, ZHANG S M, FAN W. The abundance of Antarctic krill and its relationship to sea ice and sea surface temperature. Chinese Journal of Polar Research, 2012, 24(4): 352–360 [戴立峰, 張勝茂, 樊偉. 南極磷蝦資源豐度變化與海冰和表溫的關(guān)系. 極地研究, 2012, 24(4): 352–360]

ERICSON J A, HELLESSEY N, NICHOLS P D,. Seasonal and interannual variations in the fatty acid composition of adultDana, 1850 () samples derived from the Scotia Sea krill fishery. Journal of Crustacean Biology, 2018, 38(6): 1–11

FRANCISCO S C, BILLY E, ARATA J A,. Spatial and temporal dynamics of the Antarctic krill fishery in fishing hotspots in the Bransfield Strait and South Shetland Islands. Fisheries Research, 2018, 208: 157–166

FUENTES V, ALURRALDE G, MEYER B,. Glacial melting: An overlooked threat to Antarctic krill. Scientific Reports, 2016, 1–6

GARCIA S M, ROSENBERG A A. Food security and marine capture fisheries: Characteristics, trends, drivers and future perspectives. Philosophical Transactions of the Royal Society B, 2010, 365(1554): 2869–2880

HEWITT R P, WATKINS J, NAGANOBU M,. Biomass of Antarctic krill in the Scotia Sea in January/February 2000 and its use in revising an estimate of precautionary yield. Deep-Sea ResearchⅡ, 2004, 51: 1215–1236

KAWAGUCHI S, NICOL S, TAKI K,. Fishing ground selection in the Antarctic krill fishery: Trends in patterns across years, seasons and nations. CCAMLR Science, 2006, 13: 117–141

KAWAGUCHI S, NICOL S. Learning about Antartic krill from the fishery. Antarctic Science, 2007, 19(2): 219–230

KAWAGUCHI S, CANDY S G. Quantifying movement behaviour of vessels in the Antarctic krill fishery. CCAMLR Science, 2009, 16: 131–148

KAWAGUCHI S, KING R, MEIJERS R,. An experimental aquarium for observing the schooling behaviour of Antarctic krill (). Deep-Sea Research PartⅡ, 2010a, 57(7–8): 683–692

KAWAGUCHI S, NICOL S, PRESS A J. Direct effects of climate change on the Antarctic krill fishery. Fisheries Management and Ecology, 2010b, 16: 424–427

KRAFFT B A, SKARET G, KNUTSEN T. An Antarctic krill () hotspot: Population characteristics, abundance and vertical structure explored from a krill fishing vessel. Polar Biology, 2015, 38(10): 1687–1700

LEFEVER D W. Measuring geographic concentration by means of the standard deviational ellipse. American Journal of Sociology, 1926, 32(1): 88–94

LEHODEY P, BERTIGNAC, HAMPTON J,. El Nino Southern Oscillation and tuna in the western Pacific. Nature, 1997, 389(6652): 715–718

LI X S, ZHANG J C, ZHAO X Y,. Spatiotemporal distribution ofin CCAMLR Area 48 during 2012/2013 fishing season. Progress in Fishery Sciences, 2015, 36(4): 1–9 [李顯森, 張吉昌, 趙憲勇, 等. 2012/2013漁季CCAMLR 48區(qū)南極磷蝦()資源時空分布. 漁業(yè)科學(xué)進展, 2015, 36(4): 1–9]

MARRARI M, DALY K L, HU C. Spatial and temporal variability of SeaWiFS chlorophylldistributions west of the Antarctic Peninsula: Implications for krill production. Deep-Sea Research PartⅡ, 2008(55): 377–392

MICHAEL K, SUBERG LA, WESSELS W,. Facing Southern Ocean warming: Temperature effects on whole animal performance of Antarctic krill (). Zoology, 2021, 146: 1–31

MURPHY E J, MORRIS D J, WATKINS J L,. Scales of interaction between Antarctic krill and the environment. Springer Berlin Heidelberg, 1988

NICOL S, FOSTER J. Recent trends in the fishery for Antarctic krill. Aquatic Living Resources, 2003, 16(1): 42–45

NICOL S, FOSTER J, KAWAGUCHI S. The fishery for Antarctic krill: Recent developments. Fish and Fisheries, 2011, 13(1): 30–40

PATRICIO M A, RENZO R, áLVARO D C. Chilean Antarctic krill fishery (2011–2016). Latin American Journal of Aquatic Research, 2020, 48(2): 179–196

Polar Special Office of State Oceanic Administration. Investigation and assessment of krill and other biological resources in the waters around Antarctica. Beijing: Ocean Press, 2016, 3 [國家海洋局極地專項辦公室. 南極周邊海域磷蝦等生物資源考察與評估. 北京: 海洋出報社, 2016, 3]

SAENZ B T, AINLEY D G, DALY K L,. Drivers of concentrated predation in an Antarctic marginal-ice-zone food web. Scientific Reports, 2020, 10(1): 1–17

SEOK G C, EUN A Y, DOO H A,. Characterization of frequency and aggregation of the Antarctic krill () using acoustics. Ocean Science Journal, 2018, 53(4): 1–11

SIEGEL V. Distribution and population dynamics of: Summary of recent findings. Polar Biology, 2005, 29(1): 1–22

SIEGEL V. Biology and ecology of Antarctic krill. Advances in Polar Ecology. Springer International Publishing, 2016, 22– 28

STOWASSER G, ATKINSON A, MCGILL R,. Food web dynamics in the Scotia Sea in summer: A stable isotope study. Deep Sea Research PartⅡ, 2012, 59(1): 208–221

VEYTIA D, CORNEY S, MEINERS K M,. Circumpolar projections of Antarctic krill growth potential. Nature Climate Change, 2020, 10(6): 1–8

WANG L Z, XUE C H, WANG Y M,. Extraction of proteins with low fluoride level from Antarctic krill () and their composition analysis. Journal of Agricultural and Food Chemistry, 2011, 59(11): 6108–6112

WEINSTEIN B G, DOUBLE M, GALES N,. Identifying overlap between humpback whale foraging grounds and the Antarctic krill fishery. Biological Conservation, 2017, 210: 184–191

ZHANG J C, ZHAO X Y, WANG X L,. Numericalization of pictorial image of Antarctic krill echogram taken from a commercial echo-sounder. Progress in Fishery Sciences, 2012, 33(4): 64–71 [張吉昌, 趙憲勇, 王新良, 等. 商用探魚儀南極磷蝦聲學(xué)圖像的數(shù)值化處理. 漁業(yè)科學(xué)進展, 2012, 33(4): 64–71]

ZHANG H, ZHANG S M, CUI X S,. Spatiotemporal dynamics in the location of the fishing grounds and catch per unit effort (CPUE) for Chilean jack mackerel (Nichols, 1920) from Chinese trawl fleets on the high seas of the Southeast Pacific Ocean, 2001–2010. Journal of Applied Ichthyology, 2015, 31(4): 646–656

ZHANG Y Y, XU B, ZHANG H,. Interannual and monthly variations of catch per unit effort and the relation with sea surface temperature and chlorophyll concentration in fishing grounds (Fishing Area 48) of Antarctic krill. Chinese Journal of Ecology, 2020, 39(5): 1685–1694 [張瑛瑛, 徐博, 張衡, 等. 南極磷蝦漁場(48漁區(qū))CPUE的年、月變化及其與海表溫度、葉綠素濃度的關(guān)系. 生態(tài)學(xué)雜志, 2020, 39(5): 1685–1694]

ZHAO G Q. Study on the spatial-temporal distribution of Antarctic krill () fishing ground and its relationship with environmental factors. Master′s Thesis of Yantai University, 2018, 6 [趙國慶. 南極磷蝦在漁場中的分布及與環(huán)境因子之間關(guān)系研究. 煙臺大學(xué)碩士研究生學(xué)位論文, 2018, 6]

ZHAO X Y, ZUO T, LENG K L,. Engineering science and technology challenges in the Antarctic krill fishery. Engineering Sciences, 2016, 18(2): 85–90 [趙憲勇, 左濤, 冷凱良, 等. 南極磷蝦漁業(yè)發(fā)展的工程科技需求. 中國工程科學(xué), 2016, 18(2): 85–90]

Temporal and Spatial Distribution of Antarctic Krill in 48 Fishing Areas Based on Fishery Data

ZHAO Guoqing1,2, LUO Junrong3, TANG Fenghua2, FAN Wei2, SONG Xuefeng2, YANG Chao1,2, ZHANG Heng2①

(1. College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China; 2. East China Sea Fisheries Research Institute, Chinese Academy of Fishery Science, Shanghai 200090, China;3. Fujian Zhengguan Fishery Development Coperation, Fuzhou, Fujian 350500, China)

Antarctic krill () is a vital part of the food web, supporting substantial fishery in the Southern Ocean. Krill represent the world′s largest animal protein pool, rich in lipids, nutritional long-chain polyunsaturated fatty acids, mineral elements, and carotene, with great potential for development. It is of great significance to understand the changes in fishing grounds, as they profoundly affect the krill population. The distribution characteristics of Antarctic krill fishing grounds in 48 areas were analyzed based on the production statistics of Antarctic krill fishing vessels in China from 2010 to 2019. The gravity center migration trajectory model and standard deviational ellipse (SDE) model were used to analyze the trends and characteristics of Antarctic krill fishing grounds. The results showed that the production of Antarctic krill was mainly concentrated in sub-region 48.1, accounting for 70.30% of total production. The yields of the 48.2 and 48.3 sub-regions showed little difference, accounting for 14.28% and 15.42%, respectively. The annual catch per unit effort (CPUE) curve increased over time, with a minimum value in 2012 and a maximum value in 2019. The monthly CPUE first increased and then decreased, with the minimum value in January and the maximum value in June. The annual and monthly gravity center of Antarctic krill in sub-region 48.1 moved southwest. The annual gravity center of Antarctic krill in sub-region 48.2 moved to the east, but the range of movement was small, and the monthly trends were weak. The annual gravity center of Antarctic krill in sub-region 48.3 moved to the south, and the monthly gravity center moved to the northwest. According to SDE analysis, krill fishing grounds in sub-region 48.1 had the widest distribution and the greatest dispersion, with the strongest directivity and the most obvious centripetal force in sub-region 48.3. The center of gravity of the fishing ground was distributed along the Bransfield Strait in sub-region 48.1, east of the South Orkney Islands in sub-region 48.2, and northeast of the South Georgia Islands in sub-region 48.3. Clustering results showed that the annual fishing ground gravity centers were concentrated in sub-region 48.1. The annual fishing ground gravity centers of sub-regions 48.2 and 48.3 were relatively concentrated in all years, except in 2017.

Antarctic krill; Fishing ground gravity center; Standard deviation ellipse analysis; Cluster analysis; Spatial and temporal distribution; 48 fishing areas

ZHANG Heng, E-mail: zhangziqian0601@163.com

10.19663/j.issn2095-9869.20210407004

S931

A

2095-9869(2022)04-0081-12

*國家重點研發(fā)計劃(2018YFC1406802)資助 [This work was supported by National Key Research and Development (2018YFC1406802)].趙國慶，E-mail: zgq617717@163.com

張 衡，副研究員，E-mail: zhangziqian0601@163.com

2021-04-07,

2021-06-09

http://www.yykxjz.cn/

趙國慶, 羅俊榮, 唐峰華, 樊偉, 宋學(xué)鋒, 楊超, 張衡. 基于漁業(yè)數(shù)據(jù)的南極磷蝦48漁區(qū)漁場時空分布. 漁業(yè)科學(xué)進展, 2022, 43(4): 81–92

ZHAO G Q, LUO J R, TANG F H, FAN W, SONG X F, YANG C, ZHANG H. Temporal and spatial distribution of Antarctic krill in 48 fishing areas based on fishery data. Progress in Fishery Sciences, 2022, 43(4): 81–92

(編輯 馮小花)