北太平洋長(zhǎng)鰭金槍魚(yú)年齡與生長(zhǎng)初步研究

江建軍，許柳雄，2，3，朱國(guó)平，2，3，朱偉俊

（1.上海海洋大學(xué)海洋科學(xué)學(xué)院，上海 201306；2.大洋漁業(yè)資源可持續(xù)開(kāi)發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn)室，上海海洋大學(xué)，上海 201306；3.上海市高校大洋生物資源開(kāi)發(fā)和利用重點(diǎn)實(shí)驗(yàn)室，上海海洋大學(xué)海洋科學(xué)學(xué)院，上海 201306）

北太平洋長(zhǎng)鰭金槍魚(yú)年齡與生長(zhǎng)初步研究

江建軍1，許柳雄1，2，3，朱國(guó)平1，2，3，朱偉俊1

（1.上海海洋大學(xué)海洋科學(xué)學(xué)院，上海 201306；2.大洋漁業(yè)資源可持續(xù)開(kāi)發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn)室，上海海洋大學(xué)，上海 201306；3.上海市高校大洋生物資源開(kāi)發(fā)和利用重點(diǎn)實(shí)驗(yàn)室，上海海洋大學(xué)海洋科學(xué)學(xué)院，上海 201306）

北太平洋長(zhǎng)鰭金槍魚(yú)（Thunnus alalunga）是金槍魚(yú)漁業(yè)的主要捕撈對(duì)象之一，其資源的養(yǎng)護(hù)和管理日益受到各方關(guān)注，魚(yú)類年齡鑒定是漁業(yè)資源研究的重要內(nèi)容之一。根據(jù)2013～2014年北太平洋長(zhǎng)鰭金槍魚(yú)探捕所采集的227個(gè)脊椎骨樣本，用不同的染色方法分成4組進(jìn)行年齡的鑒定，選用線性模型、指數(shù)模型和冪指數(shù)模型對(duì)叉長(zhǎng)和脊椎骨半徑進(jìn)行擬合，最后選出結(jié)果最好的鑒定組和擬合模型，通過(guò)年齡鑒定估算長(zhǎng)鰭金槍魚(yú)的Von Bertalanffy生長(zhǎng)方程。研究結(jié)果表明，北太平洋長(zhǎng)鰭金槍魚(yú)的叉長(zhǎng)范圍為52～104 cm，優(yōu)勢(shì)叉長(zhǎng)組為81～90 cm，占總數(shù)的83.1%；年齡組成以5、6齡較多，占樣本總數(shù)的80.86%。最終得到北太平洋長(zhǎng)鰭金槍魚(yú)生長(zhǎng)方程為：Lt＝109.180 2［1－e－0.258783（t－1.99128）］，ttp＝2.25。

北太平洋；長(zhǎng)鰭金槍魚(yú)；脊椎骨；年齡鑒定；生長(zhǎng)

長(zhǎng)鰭金槍魚(yú)（Thunnus alalunga）具有高經(jīng)濟(jì)價(jià)值，其產(chǎn)量占全球金槍魚(yú)產(chǎn)量的7%左右，是海洋漁業(yè)資源的重點(diǎn)養(yǎng)護(hù)對(duì)象。長(zhǎng)鰭金槍魚(yú)主要分布于大西洋、印度洋和太平洋的熱帶和溫帶55°N～45°S之間的海域，但在10°N～10°S之間的表層水域很少出現(xiàn)。太平洋長(zhǎng)鰭金槍魚(yú)主要分為南部長(zhǎng)鰭（10°S以南）和北部長(zhǎng)鰭（10°N以北）兩個(gè)種群，而南北緯度10°之間基本上沒(méi)有出現(xiàn)，一般捕撈方式以延繩釣和流刺網(wǎng)等為主。太平洋長(zhǎng)鰭金槍魚(yú)是金槍魚(yú)漁業(yè)的主要捕撈對(duì)象之一，其資源的養(yǎng)護(hù)和管理日益受到各方關(guān)注。自1952年以來(lái)北太平洋長(zhǎng)鰭金槍魚(yú)的捕獲量大約為75 000 t/年，其中日本所占比例最大，約占72%左右［1］。

魚(yú)類年齡鑒定是漁業(yè)資源研究的重要內(nèi)容之一，有關(guān)信息是研究魚(yú)類種群動(dòng)力學(xué)所需的基本生物學(xué)參數(shù)。利用金槍魚(yú)的脊椎骨對(duì)金槍魚(yú)進(jìn)行年齡鑒定在對(duì)魚(yú)類年齡鑒定中是一種比較常見(jiàn)的方法。國(guó)內(nèi)外許多研究大多采用金槍魚(yú)的第35節(jié)脊椎骨來(lái)作為鑒定其年齡的主要材料［1］。也有許多專家和學(xué)者采用耳石［1－2］和鰭條［2］對(duì)金槍魚(yú)進(jìn)行年齡鑒定。

迄今國(guó)內(nèi)學(xué)者關(guān)于北太平洋長(zhǎng)鰭金槍魚(yú)年齡與生長(zhǎng)的研究相對(duì)較少，而國(guó)外學(xué)者研究雖然相對(duì)較多，但發(fā)表時(shí)間都較久遠(yuǎn)。本文利用2013～2014年北太平洋長(zhǎng)鰭金槍魚(yú)資源漁場(chǎng)探捕調(diào)查數(shù)據(jù)，研究結(jié)果是北太平洋長(zhǎng)鰭金槍魚(yú)生物學(xué)信息的當(dāng)下體現(xiàn)，以期為該物種的資源評(píng)估和養(yǎng)護(hù)提供重要的基礎(chǔ)數(shù)據(jù)。

1 材料與方法

1.1 材料來(lái)源

調(diào)查漁船為大連遠(yuǎn)洋金槍魚(yú)釣有限公司“天祥16號(hào)”超低溫金槍魚(yú)延繩釣漁船。調(diào)查時(shí)間為2013年10月19日～2014年4月5日。調(diào)查海域?yàn)椋?0°～42°N、143°～164°W。調(diào)查海域?yàn)樯a(chǎn)海域，不影響船只的正常生產(chǎn)。觀察員在船上記錄釣獲的長(zhǎng)鰭金槍魚(yú)的體長(zhǎng)、體質(zhì)量和性別等數(shù)據(jù)，取下其第35節(jié)脊椎骨，并做冷凍處理后保存，隨船運(yùn)回上海海洋大學(xué)實(shí)驗(yàn)室。

1.2 樣本的染色

脊椎骨實(shí)驗(yàn)室預(yù)處理。將脊椎骨在水中煮5～6 min，當(dāng)附著的肌肉變軟時(shí)，用鋼絲球或者刷子將其擦除后進(jìn)行編號(hào)，然后用清水洗凈，室溫下風(fēng)干［2］。

按照文獻(xiàn)［2］對(duì)風(fēng)干的脊椎骨進(jìn)行染色，并對(duì)藥劑的選擇、染色步驟和染色時(shí)間進(jìn)行改進(jìn)，以更好地探索輪紋清晰呈現(xiàn)的方法。本文將染色方案分為4組，樣本數(shù)量分別為62個(gè)、54個(gè)、60個(gè)和51個(gè)，共計(jì)227個(gè)；其中雄性115個(gè)，雌性112個(gè)，具體染色步驟如表1所示。第一和第二試驗(yàn)組僅改變茜素紅的濃度，其余步驟完全相同；第三和第四試驗(yàn)組只是茜素紅的濃度做了改變，其余步驟完全相同。第一和第三試驗(yàn)組的茜素紅濃度相同，只是第一組的酒精濃度梯度做了一定的改變；第二和第四試驗(yàn)組的茜素紅濃度相同，只是第二組的酒精濃度梯度做了一定的改變，且和第一組相同。

圖1 樣本采集站位圖Fig.1 Sam pling stations

表1 4組試驗(yàn)的染色方法Tab.1 Dyeing method of four groups

脊椎骨輪紋讀取由2人完成。1號(hào)為主要讀者，一般具有讀取輪紋的經(jīng)驗(yàn)；2號(hào)為輔助讀者［2］。1號(hào)讀取兩次，2號(hào)讀取1次。每次讀取輪紋時(shí)，不向讀者提供長(zhǎng)鰭金槍魚(yú)的叉長(zhǎng)（fork length，F(xiàn)L）和其它相關(guān)信息。1號(hào)讀者兩次讀取輪紋的時(shí)間間隔至少為2周，避免前一次對(duì)輪紋的記憶而產(chǎn)生影響。如果某個(gè)樣本3次讀取的年齡都不一樣，則舍棄；如果兩次讀數(shù)相同，另一次讀數(shù)的誤差在1齡以內(nèi)，則以兩次相同的讀數(shù)為準(zhǔn)，反之則舍棄［3］。

1.3 數(shù)據(jù)的分析與處理

1.3.1 叉長(zhǎng)頻次分析

長(zhǎng)鰭金槍魚(yú)叉長(zhǎng)以5 cm組間距表示，采用t檢驗(yàn)分析其差異性；應(yīng)用協(xié)方差分析（analysis of covariance，ANCOVA）檢驗(yàn)雌、雄性長(zhǎng)鰭金槍魚(yú)的叉長(zhǎng)與脊椎骨半徑（圖2）（vertebral radius，VR）之間是否呈現(xiàn)顯著性差異。

圖2 脊椎骨半徑和輪紋示意圖Fig.2 Radius of vertebrae and wheel lines diagram

1.3.2 各組平均百分比誤差指數(shù)的差異與對(duì)比

平均百分比誤差指數(shù)［3］是指不同鑒定者之間年齡鑒定的重復(fù)性指標(biāo)，當(dāng)平均百分比誤差指數(shù)小于10%時(shí)，表明年輪鑒定精確度較高，其結(jié)果是可信的，否則需重新鑒定或者作廢。計(jì)算公式如下：

式（1）中：N為第j樣本重復(fù)計(jì)算的次數(shù)；xij為第j樣本第i次計(jì)數(shù)的輪紋數(shù)目；xj為第j樣本重復(fù)計(jì)數(shù)的平均輪紋數(shù)。

由于本文將染色分為4組進(jìn)行，所以各組的染色效果有可能是明顯不同的。因此，對(duì)比其平均百分比誤差指數(shù)來(lái)判定最終的染色理想程度是一個(gè)重要的指標(biāo)。

1.3.3 叉長(zhǎng)和脊椎骨半徑關(guān)系

分別通過(guò)冪函數(shù)、線性方程及指數(shù)方程擬合叉長(zhǎng)和脊椎骨半徑之間的關(guān)系［4］。

式（2）～（4）中：FL為長(zhǎng)鰭金槍魚(yú)的叉長(zhǎng)（cm），R為脊椎骨半徑（mm），a1、a2、a3、b1、b3均為回歸系數(shù)，b2為截距。

選用赤遲信息量標(biāo)準(zhǔn)（AIC）［5］對(duì)3種模型結(jié)果進(jìn)行檢驗(yàn)，以尋求最合適的模型。AIC值最適合相對(duì)較小的模型，可簡(jiǎn)化表達(dá)式為：

式（5）中：n為樣本數(shù)目，RSS為殘差平方和［6］，即規(guī)劃求解后的最小值，k為模型參數(shù)個(gè)數(shù)。

1.3.4 叉長(zhǎng)的逆算

選取AIC值最小的模型進(jìn)行長(zhǎng)鰭金槍魚(yú)逆算叉長(zhǎng)計(jì)算。其中，如果指數(shù)模型最適合，可以直接通過(guò)指數(shù)回歸關(guān)系計(jì)算，而線性模型和冪函數(shù)模型則需要運(yùn)用Fraser-Lee’s法和Monastyrsky法計(jì)算［7］，其公式分別為：

Fraser-Lee’s法：

式（6）～（7）中：FL為長(zhǎng)鰭金槍魚(yú)被釣獲時(shí)測(cè)定的叉長(zhǎng)（mm），Ln為第n輪年輪形成時(shí)的逆算叉長(zhǎng)（mm），R為脊柱骨半徑，Rn為各年齡所對(duì)應(yīng)的半徑（mm），b1為脊柱骨半徑-叉長(zhǎng)線性回歸關(guān)系式中的截距，b2為脊椎骨半徑-叉長(zhǎng)的冪函數(shù)回歸系數(shù)。

1.3.5 長(zhǎng)鰭金槍魚(yú)生長(zhǎng)方程

將所得最佳模型的逆算叉長(zhǎng)值代入Von Bertalanffy生長(zhǎng)方程，利用Excel2013中的規(guī)劃求解功能估算其生長(zhǎng)參數(shù)的最佳值。Von Bertalanffy生長(zhǎng)方程是描述金槍魚(yú)類生長(zhǎng)最重要的模型，很多學(xué)者都用此方法來(lái)描述金槍魚(yú)的生長(zhǎng)［7］。Von Bertalanffy生長(zhǎng)方程為：

式（8）中：Lt為t齡時(shí)的叉長(zhǎng)；L∞為極限體長(zhǎng)；K為瞬時(shí)相對(duì)生長(zhǎng)速率；t為年齡；t0為叉長(zhǎng)為0時(shí)對(duì)應(yīng)的理論年齡。

長(zhǎng)鰭金槍魚(yú)的Von Bertalanffy生長(zhǎng)曲線的拐點(diǎn)年齡為：

式（9）中：ttp為生長(zhǎng)拐點(diǎn)；lnb為參數(shù)；k為瞬時(shí)相對(duì)生長(zhǎng)速率；t0為叉長(zhǎng)為0時(shí)對(duì)應(yīng)的理論年齡。

2 結(jié)果與分析

2.1 叉長(zhǎng)頻次

實(shí)驗(yàn)采集的有效樣本209個(gè)，叉長(zhǎng)范圍為52～104 cm，平均叉長(zhǎng)為（84.79±4.951）cm，優(yōu)勢(shì)叉長(zhǎng)為81～90 cm，占總數(shù)的83.1%，其中雌性叉長(zhǎng)范圍為63～104 cm，平均叉長(zhǎng)為（84.43± 5.274）cm，優(yōu)勢(shì)叉長(zhǎng)為81～90 cm，占雌性總數(shù)的81.82%；雄性叉長(zhǎng)范圍為52～96 cm，平均叉長(zhǎng)為（85.23±4.642）cm，優(yōu)勢(shì)叉長(zhǎng)為81～90 cm，占雄性總數(shù)的84.47%。

圖3 北太平洋長(zhǎng)鰭金槍魚(yú)叉長(zhǎng)組頻次分布Fig.3 Fork length frequency distribution of the north Pacific albacore tuna

2.2 年齡組成和平均百分比誤差指數(shù)分析

對(duì)染色后的227個(gè)脊椎骨樣本進(jìn)行了鑒定，其中達(dá)到鑒定要求的209個(gè)。在沒(méi)有達(dá)到要求的18個(gè)樣品中，第一組5個(gè)，第二組2個(gè)，第三組6個(gè)，第四組5個(gè)。各組平均百分比誤差指數(shù)如表2所示，APE2＜APE1＜APE4＜APE3，表明第一試驗(yàn)組和第二試驗(yàn)組的染色效果明顯好于和第四試驗(yàn)組。

表2 各試驗(yàn)組平均百分比誤差指數(shù)Tab.2 Average percentage error index

在達(dá)到鑒定要求的209個(gè)樣本中，雄性樣本106個(gè)，雌性樣本103個(gè)；樣本叉長(zhǎng)范圍52～104 cm，年齡組成在1～9a，其中5～6a共169個(gè)，占總樣本數(shù)的80.86%（圖4，圖5，表3）。協(xié)方差檢驗(yàn)結(jié)果表明，2013年和2014年雌、雄長(zhǎng)鰭金槍魚(yú)叉長(zhǎng)與脊椎骨半徑的關(guān)系均無(wú)顯著性差異（P＞0.05），故分析時(shí)將上述兩年的雌、雄長(zhǎng)鰭金槍魚(yú)樣本數(shù)據(jù)合并。

圖4 北太平洋長(zhǎng)鰭金槍魚(yú)叉長(zhǎng)與年齡的關(guān)系Fig.4 Relationship between fork length and age of north Pacific albacore tuna

圖5 3～9a各年齡段的染色樣本Fig.5 Dyeing sam p les of all ages

表3 各年齡階段的叉長(zhǎng)范圍和樣本個(gè)數(shù)Tab.3 Fork length com position by age of the north Pacific albacore tuna

2.3 叉長(zhǎng)和脊椎骨半徑的關(guān)系

對(duì)于達(dá)到鑒定要求的209個(gè)樣本，利用圓規(guī)和直尺逐個(gè)測(cè)量其脊柱骨半徑和各年輪對(duì)應(yīng)的半徑，然后分別通過(guò)冪函數(shù)、線性方程和指數(shù)方程擬合長(zhǎng)鰭金槍魚(yú)叉長(zhǎng)和脊椎骨半徑的關(guān)系，選擇相關(guān)系數(shù)R2值最大和AIC值最小的方程作為最適方程（表4）。

2.4 叉長(zhǎng)的逆算

4組實(shí)驗(yàn)結(jié)果發(fā)現(xiàn)指數(shù)模型的AIC值最小，表明選取指數(shù)模型擬合長(zhǎng)鰭金槍魚(yú)叉長(zhǎng)和脊椎骨半徑之間的關(guān)系最合適，根據(jù)脊椎骨半徑可直接利用指數(shù)模型對(duì)長(zhǎng)鰭金槍魚(yú)叉長(zhǎng)直接進(jìn)行逆運(yùn)算，見(jiàn)圖6。

表4 叉長(zhǎng)和脊椎骨半徑的關(guān)系Tab.4 Relationship between fork length and radius of vertebrae of north Pacific albacore tuna

圖6 各試驗(yàn)組叉長(zhǎng)和脊椎骨半徑的實(shí)測(cè)值和擬合值Fig.6 M easured values and fitted values of fork length and radius of backbones of north Pacific albacore tuna for each experimental group

表5 第二試驗(yàn)組各年齡長(zhǎng)鰭金槍魚(yú)的平均逆算叉長(zhǎng)Tab.5 Average inverse calculated fork length of north Pacific albacore by age based on the second experimental group

單因子方差分析（ANOVA）表明，4個(gè)試驗(yàn)組在鑒定長(zhǎng)鰭金槍魚(yú)年齡時(shí)不存在顯著差異（P＞0.05）。然而，第二試驗(yàn)組的指數(shù)模型AIC（24.736 91）明顯小于其它3個(gè)試驗(yàn)組（表4），表明該試驗(yàn)組比其它3個(gè)試驗(yàn)組能更好地?cái)M合長(zhǎng)鰭金槍魚(yú)年齡與叉長(zhǎng)關(guān)系?；诘诙囼?yàn)組各年齡的長(zhǎng)鰭金槍魚(yú)平均逆算叉長(zhǎng)如表5所示。

2.5 生長(zhǎng)方程

根據(jù)第二試驗(yàn)組所得的1～8齡長(zhǎng)鰭金槍魚(yú)逆算叉長(zhǎng)數(shù)據(jù)，應(yīng)用Von Bertalanffy生長(zhǎng)模型，通過(guò)規(guī)劃求解擬合了北太平洋長(zhǎng)鰭金槍魚(yú)的年齡與生長(zhǎng)，確定最終的生長(zhǎng)方程如下：

拐點(diǎn)年齡ttp為2.25。

3 討論

3.1 脊椎骨的染色

為了探索使輪紋得以清晰呈現(xiàn)的更好方法，本文根據(jù)藥劑選擇、染色步驟、染色時(shí)間等，分四組對(duì)脊椎骨進(jìn)行染色。鑒于染色淺了不利觀察，染色深了則會(huì)將一些年輪覆蓋掉，因此染色過(guò)程中尤其是在對(duì)酒精濃度和茜素紅濃度的掌握上有很高的要求。為了更好地對(duì)比實(shí)驗(yàn)結(jié)果，實(shí)驗(yàn)采取了改變單一變量的方法。文章采用了冪函數(shù)、線性和指數(shù)函數(shù)三種方法來(lái)擬合長(zhǎng)鰭金槍魚(yú)叉長(zhǎng)和脊椎骨的關(guān)系，通過(guò)對(duì)比AIC的值最終選擇了指數(shù)模型來(lái)作為逆算叉長(zhǎng)公式，因?yàn)橹笖?shù)函數(shù)的AIC要明顯小于其它兩種模型，且指數(shù)模型的相關(guān)系數(shù)最大。另外，不論在平均百分比誤差指數(shù)還是在AIC上，第二組試驗(yàn)的結(jié)果都要明顯優(yōu)于其它三組。所以實(shí)驗(yàn)二的染色方法相對(duì)其它三組更加適合對(duì)長(zhǎng)鰭金槍魚(yú)進(jìn)行染色，達(dá)到的染色效果也更加明顯。

脊椎骨經(jīng)染色的最大優(yōu)點(diǎn)就是，經(jīng)染色以后的脊椎骨輪紋清晰明顯，輪紋完整，十分適合進(jìn)行年齡的鑒定。在染色過(guò)程中也遇到了許多問(wèn)題，染色的前期處理比較麻煩，再就是在對(duì)編號(hào)的處理上，由于整個(gè)染色過(guò)程比較長(zhǎng)，大約會(huì)持續(xù)一周的時(shí)間［2］，因此，在這個(gè)過(guò)程中標(biāo)記的繩索容易脫落，這樣會(huì)導(dǎo)致在后面鑒定年齡時(shí)找不到對(duì)應(yīng)的編號(hào)，造成樣品的浪費(fèi)，以后可以在該方面進(jìn)行完善。

3.2 叉長(zhǎng)和脊椎骨半徑關(guān)系

在對(duì)叉長(zhǎng)和脊椎骨半徑關(guān)系進(jìn)行分析的時(shí)候，文章采用了線性模型、指數(shù)模型以及冪函數(shù)模型來(lái)擬合叉長(zhǎng)和脊椎骨半徑的關(guān)系，最終得出指數(shù)模型為最合適的函數(shù)。由于不同種類、不同區(qū)域以及不同樣本的魚(yú)類都會(huì)對(duì)模型的結(jié)果產(chǎn)生直接性的影響。所以，為了找到最合適的模型來(lái)模擬長(zhǎng)鰭金槍魚(yú)的生長(zhǎng)，應(yīng)該盡可能多地應(yīng)用更多的模型來(lái)分析。文章在最初對(duì)叉長(zhǎng)和脊椎骨關(guān)系進(jìn)行分析時(shí)，也曾經(jīng)采用過(guò)對(duì)數(shù)模型和多項(xiàng)式模型，并且發(fā)現(xiàn)多項(xiàng)式模型得到的結(jié)果是最優(yōu)的。但在最后進(jìn)行叉長(zhǎng)逆運(yùn)算時(shí)，經(jīng)多方面查找也沒(méi)有找到合適的逆運(yùn)算方程，國(guó)內(nèi)外相關(guān)領(lǐng)域的專家也沒(méi)有使用過(guò)這兩種模型的，所以最終放棄了這兩種模型。

3.3 生長(zhǎng)方程

國(guó)內(nèi)外許多專家學(xué)者對(duì)長(zhǎng)鰭金槍魚(yú)的生長(zhǎng)進(jìn)行過(guò)研究（表6），研究結(jié)果表明長(zhǎng)鰭金槍魚(yú)極限叉長(zhǎng)最小值為102.9 cm，最大的為146.5 cm，本研究的極限叉長(zhǎng)為109.1 cm。生長(zhǎng)曲線的平均曲率最大為0.4，最小為0.134，本研究為0.258。t0最大為1.999，最小為－2.273，本研究為－1.991，本研究的所有參數(shù)都在范圍之內(nèi)。NOSE等［8］、OTSU等［9］對(duì)金槍魚(yú)的年齡鑒定過(guò)大，導(dǎo)致最后研究結(jié)果t0的值大于零。CLEMENS等［10］實(shí)驗(yàn)樣本太少，且叉長(zhǎng)范圍太小，導(dǎo)致研究結(jié)果的極限叉長(zhǎng)過(guò)大，方程到高齡部分不能進(jìn)行有效的收斂。由于樣本原因，可能造成不同的結(jié)果。本研究結(jié)果與BELL［11］在1962年通過(guò)耳石對(duì)北太平洋長(zhǎng)鰭金槍魚(yú)的研究結(jié)果比較相近，在一定程度上說(shuō)明本研究結(jié)果是有理論依據(jù)的，研究結(jié)果是可信的。

表6 國(guó)內(nèi)外學(xué)者關(guān)于長(zhǎng)鰭金槍魚(yú)生長(zhǎng)方程的研究結(jié)果Tab.6 Research results about albacore tuna grow th equation by fisheries researchers hom e and abroad

3.4 不足與展望

本研究采用的總樣本數(shù)為227個(gè)，由于分成四組進(jìn)行研究，每組的樣本數(shù)量就相對(duì)較少，在對(duì)每一組進(jìn)行數(shù)據(jù)分析時(shí)，顯得略微不足。另外，高齡魚(yú)與低齡魚(yú)的數(shù)量相對(duì)較少，這在進(jìn)行逆算叉長(zhǎng)和生長(zhǎng)方程的分析時(shí)會(huì)有一定的局限性，不能從統(tǒng)計(jì)學(xué)上完全地反映其最適結(jié)果，這也是國(guó)內(nèi)外在年齡與生長(zhǎng)的研究方面共同存在的缺陷，望以后在做漁業(yè)調(diào)查時(shí)可以彌補(bǔ)這方面的缺陷。

［1］ WELLSR JD，KOHIN S，TEO S L H，et al.Age and growth of North Pacific albacore（Thunnusalalunga）：Implications for stock assessment［J］.Fisheries Research 2013：55－62.

［2］ 李鳳瑩，印度洋中南部大眼金槍魚(yú)年齡與生長(zhǎng)研究［D］.上海海洋大學(xué)，2014.

LI F Y.The age and growth of bigeye tuna（Thunnusobesus）in the South-central Indian Ocean［D］.ShanghaiOcean University，2014.

［3］ BEAMISH R J，F(xiàn)OURNIER D A.A method for comparing the precision of a set of age determinations［J］.Canadian Journal of Fisheries and Aquatic Sciences，1981，38（8）：982－983.

［4］ 周 成，朱國(guó)平，陳錦淘，等.印度洋南部大眼金槍魚(yú)年齡鑒定及其與生長(zhǎng)的關(guān)系［J］.中國(guó)水產(chǎn)科學(xué)，2012，19（3）：536－544.

ZHOU C，ZHU G P，CHEN J T，et al.Relationship between age and growth of bigeye tuna（Thunnusobesus）in the southern Indian Ocean［J］.Journal of Fishery Sciences of China，2012，19（3）：536－544.

［5］ AKAIKE H.Information measures and model selection［J］.The International Statistical Institute，1983，50（1）：277－291.

［6］ CHENY，JACKSONDA，HARVEYH H.A comparison of Von Bertalanffy and polynomial functions in modelling fish growth data［J］.Canadian Journal of Fisheries and Aquatic Sciences，1992，49（6）：1128－1235.

［7］ SUN C，WANG S，YEH S.Age and growth of the swordfish（XiphiasgladiusL.）in the waters around Taiwan determined from anal-fin rays［R］.Fishery bulletin，2002，100（4）：822－835.

［8］ NOSEY，KAWATSUH，HIYAMA Y.Age and growth of Pacific tunas by scale reading［J］.Tokyo University Press，1967：701－716.

［9］ OTSU T.Albacoremigration and growth in the North Pacific Ocean as estimated from tag recoveries［J］.Pacific Science，1960：14（3）：257－266.

［10］ CLEMENS H B.The migration age and growth of Pacific albacore（Thunnus germo）［J］.State of California Department of Fish and Game Fish Bulletin，1961：115－128.

［11］ BELL R R.Age determination of the Pacific albacore of the California coast［J］.California Fish and Game，1962（48）：39－48.

［12］ YABUTA Y，YUKINAWA M.Growth and age of albacore［J］.Report of Nankai Reginoal Fisheries Research Laboratory，1963（17）：111－120.

［13］ LABELLE，HAMPTON J，BAILEY K，et al.Determination of age and growth of South Pacific albacore（Thunnus alalunga）using threemethodologies［J］.Fishery Bulletin，1993（91）：649－663.

［14］ SANTIAGO J，ARRIZIBALAGA H.An integrated growth study for North Atlantic albacore（Thunnus alalungaBonn.1788）［J］.ICESJournalofMarine Science，2005（62）：740－749.

［15］ FARLEY J，CLEAR N.Preliminary study of age growth and spawning activity of albacore in Australia’s Eastern Tuna and Billfish Fishery［C］.Western and Central Pacific Fisheries Commission， 2008：36.

［16］ WILLIAMS A J，F(xiàn)ARLEY JH，HOYLE S D，et al.Spatial and sex-specific variation in growth of albacore tuna across the South Pacific Ocean.Plos One，2012，7（6）：9063－9070.

［17］ CHEN K S，SHIMOSE T，TANABE T，etal.Age and growth of albacore（Thunnusalalunga）in the North Pacific Ocean［J］.Journal of Fish Biology，2012（80）：2328－2344.

Prelim inary study on age and grow th of north Pacific albacore Thunnus alalunga

JIANG Jian-jun1，XU Liu-xiong1，2，3，ZHU Guo-ping1，2，3，ZHUWei-jun1
（1.College of Marine Sciences，Shanghai Ocean University，Shanghai201306，China；2.Key Laboratory of Oceanic Fisheries Resources Exploitation of Shanghai Education Commission，Shanghai Ocean University，Shanghai201306，China；3.Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources，Ministry of Education，Shanghai Ocean University，Shanghai201306，China）

The north Pacific albacore Thunnusalalunga is one of themain target species of tuna fisheries，the conservation and management of its resources has caught increasing attention by regional tuna management organizations.The fish age identification is an important contentof fishery resources research.There is a need for a simple practical and accurate agingmethod to document the population structure，age and growth of the north Pacific albacore.We collected 227 vertebrae samples of albacore in the North Pacific between 2013 and 2014，divided them into four groups for age identification in the experiment by different dyeingmethods and then used linear model，exponentialmodel and power index model to fit the relationship between the fork length and radius of vertebrae，the best historical and fittingmodelwas chosen to calculate the Von Bertalanffy growth equation for the North Pacific albacore.Results indicated that the fork length of albacore ranged from 52 cm to 104 cm，with the average fork length of84.79±4.951 cm，the dominant fork length（83.1%）was between 81－90 mm.The fork length of female albacore ranged from 63 cm to 104 cm，with the average fork length 84.43±5.274 cm and the dominant length（81.82%）was between 81－90 mm.The fork length of male albacore ranged from 52 cm to 96 cm，with the average fork length 85.23±4.642 cm and the dominant fork length（80.86%）was between 81－90 mm.The 5 and 6 age class accounted for 80.86%of the total number sampled.The single factor analysis of variance showes that there is no significant difference among four experiments in the identification of albacore age.However，the AIC indexmodel of experimental group is lower than other three test groups，and it shows that the experimental group can better fit age-fork length relationship for albacore tuna.The Von Bertalanffy growth function is，Lt＝109.180 2［1－e－0.258783（t－1.99128）］，ttp＝2.25.

north Pacific；albacore tuna；vertebrae；age identification；growth

S 932.4

A

1004－2490（2016）01－0001－08

2015－07－14

農(nóng)業(yè)部財(cái)政部北太平洋長(zhǎng)鰭金槍魚(yú)探捕項(xiàng)目

江建軍，男，碩士研究生，主要從事遠(yuǎn)洋漁業(yè)資源研究。E-mail：827324606＠qq.com

許柳雄，教授。E-mail：lxxu＠shou.edu.cn