麻秋云 薛 瑩 徐賓鐸 紀毓鵬 陳學剛韓東燕 孫遠遠
(1.中國海洋大學水產學院, 青島 266003; 2.青島市黃島區(qū)海洋與漁業(yè)局, 青島 266500)
耳石是硬骨魚類在生長過程中沉積在內耳中的結石, 由于其記錄的信息涵蓋了魚類生理生態(tài)學研究的諸多方面, 因此已逐漸成為國際上漁業(yè)科學和海洋環(huán)境科學研究的熱點[1]。由于耳石形態(tài)結構穩(wěn)定, 具有種的特征, 利用光學圖像處理軟件和傅立葉分析等統(tǒng)計學手段, 研究者將耳石形態(tài)學分析方法應用于魚類物種及群系識別領域, 取得了良好效果[2—4]。
在魚類胃含物分析的過程中, 許多餌料魚類在胃中多處于消化或半消化狀態(tài), 主要身體部位常常殘缺不全, 給種類鑒定和個體大小的測量帶來一定的難度[5]。由于耳石具有堅硬且不易被消化的特點,即使餌料魚類的消化程度很高, 其耳石仍然能夠較為完整的保留, 因此可以依據(jù)耳石的形態(tài)特征鑒別餌料魚類的種類, 這方面已經(jīng)在許多高營養(yǎng)級海洋動物(如海洋哺乳動物、海鳥、兇猛肉食性魚類等)的食性研究中得到了廣泛的應用[6—9]。除了種類鑒別, 通過建立耳石大小與魚體長度的關系式, 還可以根據(jù)胃含物中殘留耳石的個體大小來推算餌料魚類的體長等形態(tài)指標, 近年來國外學者在該領域開展許多相關的研究[10—14], 而我國在這方面的研究報道較少。
本文選擇膠州灣 12種常見的餌料魚類作為研究對象, 分析每種魚類的體長與耳石形態(tài)大小(包括耳石的長度、寬度、長半徑和短半徑)的關系, 并建立回歸方程, 旨在為深入研究高營養(yǎng)級魚類的攝食生態(tài)和食物網(wǎng)營養(yǎng)動力學提供基礎資料。
樣品采自2008年11月至2009年8月在膠州灣海域(36°06′—36°08′N, 120°12′—120°18′E)進行的逐月定點底拖網(wǎng)調查(圖1)。調查船為30 kW左右的單拖漁船, 平均拖速為 2.0 kn, 平均每站拖網(wǎng)時間為0.5h。調查網(wǎng)具網(wǎng)口高度1.6 m, 網(wǎng)口寬度12 m, 囊網(wǎng)網(wǎng)目 20 mm。調查共采集 12種餌料魚類的樣品1063尾, 各個魚種的樣品數(shù)目及體長范圍(表1)。

圖1 膠州灣調查站位Fig.1 Sampling stations in Jiaozhou Bay
魚類樣品帶回實驗室做進一步的分析, 鑒定種類, 并按照最新的《海洋調查規(guī)范》[15]進行生物學測定, 包括體長(mm)和體重(g)等生物學參數(shù)。將魚類的左、右矢耳石取出, 清洗干燥后, 編號記錄置于離心管中保存。采用尼康E800體視顯微鏡拍攝耳石的照片, 分類編號后應用Image-Pro Plus 6.0軟件測量耳石的長度、寬度、長半徑和短半徑。其中, 耳石長度是指平行于耳石主凹槽的耳石前后端最大長度, 耳石寬度是指垂直于主凹槽從耳石背側到腹側的最大寬度[12]; 耳石長半徑和短半徑分別為耳石中心核到耳石邊緣的最大距離和最小距離[16]。
根據(jù) Keys公式W=a*SLb計算各個魚種體長與體重的關系[17], 式中W為體重(g),SL為體長(cm),a為條件系數(shù),b為回歸系數(shù)。用t檢驗法檢驗回歸系數(shù)b值與3之間是否有顯著差異[18], 若差異不顯著則為等速生長, 反之則為異速生長[17]。運用最小二乘法對 12種魚類體長與耳石形態(tài)大小(包括耳石的長度、寬度、長半徑和短半徑)的關系進行回歸分析[10,12,14], 并應用方差分析(ANOVA)檢驗回歸方程的顯著性[18]。經(jīng)檢驗, 12種魚類左右矢耳石的形態(tài)參數(shù)均無顯著差異(P<0.05), 因此本文選取各個魚種的左耳石進行分析(圖2)。
本文所研究的12種餌料魚類分別隸屬于3目, 7科, 11屬。其中, 樣品數(shù)量最多的魚種是六絲鈍尾蝦虎魚Amblychaeturichthys hexanema(n=197), 樣品數(shù)量最少的魚種為Liza haematocheila(n=43)。結果表明, 12種魚類的體長與體重均呈顯著的冪函數(shù)關系(P<0.05), 回歸系數(shù)b值的范圍是 2.871—3.371,平均值為 3.159±0.110(95% CL), 所有魚種b值的95%的置信區(qū)間均在2.5—3.5的范圍內。b值最小的魚種是六絲鈍尾蝦虎魚, 最大的是細條天竺鯛Apogon lineatus。t檢驗表明,和箭魚銜Callionymus sagittal的b值與3無顯著差異(P<0.05), 屬等速生長; 李氏魚銜Callionymus richardsonii和六絲鈍尾蝦虎魚的b值顯著小于3(P<0.05), 屬負異速生長; 其余8種魚類的b值均顯著大于3(P<0.05), 為正異速生長。

表1 膠州灣12種餌料魚類的樣品數(shù)量和體長范圍Tab.1 Sample size and size range of twelve prey fish species from Jiaozhou Bay

圖2 膠州灣12種餌料魚類矢耳石的圖像(左耳石)Fig.2 Images of sagittal otolith for 12 prey fish species in Jiaozhou Bay (left otolith)
方差分析表明, 膠州灣12種餌料魚類體長與耳石長度和寬度之間都存在極顯著的線性回歸關系(P<0.01)。在體長與耳石長度的回歸方程中, 有9種魚類的R2值大于0.8, 回歸方程的擬合效果較好, 其余3種魚類(鳀Engraulis japonicus、李氏魚銜和箭魚銜)的R2值相對較小(0.7
在12種餌料魚類中, 有 10種魚類的體長與耳石長度和耳石寬度回歸方程之間的R2值比較接近;而其余 2種魚類(鳀和箭魚銜)2個回歸方程的R2值相差較大, 擬合效果存在明顯差異, 鳀其中體長與耳石長度(R2=0.778)的擬合效果要優(yōu)于體長與耳石寬度(R2=0.643), 而箭魚銜體長與耳石寬度(R2=0.923)的擬合效果要明顯優(yōu)于體長與耳石長度(R2=0.742)(表 3)。

表2 膠州灣12種餌料魚類體長與體重回歸方程的參數(shù)及其決定系數(shù)(R2)Tab.2 Length-weight relationship parameters and determination coefficient for twelve prey fish species from Jiaozhou Bay
膠州灣 12種餌料魚類體長與耳石長半徑和短半徑之間都存在極顯著的線性回歸關系(P<0.01)。在體長與耳石長半徑的回歸方程中, 有9種魚類的R2值大于0.8, 回歸方程的擬合效果較好, 其余3種魚類(鳀、李氏魚銜和箭魚銜)的R2值相對較小(0.7
在12種餌料魚類中, 有8種魚類的體長與耳石長半徑和耳石短半徑回歸方程之間的R2值比較接近;其余4種魚類(鳀、、細條天竺鯛和箭魚銜)2個回歸方程的R2值相差較大, 擬合效果存在較大的差異,其中箭魚銜的體長與耳石短半徑(R2=0.930)的擬合效果要優(yōu)于體長與耳石長半徑(R2=0.750), 而其余 3種魚類的體長與耳石長半徑的擬合效果均優(yōu)于體長與耳石短半徑(表3)。
膠州灣位于黃海之濱, 山東半島的南岸, 灣內平均水深7 m, 通過一個狹窄的灣口(約3.1 km)與黃海相通, 是一個典型的半封閉型淺海灣[19]。由于膠州灣沿岸有大量的淡水流入, 帶來豐富的營養(yǎng)鹽, 使得灣內水域生產力高, 餌料生物豐富, 是多種經(jīng)濟魚、蝦、蟹類繁殖、育幼和索餌的場所[19]。本文所研究的 12種魚類均是目前膠州灣海域數(shù)量較為豐富的小型餌料魚類[20,21], 它們既是小型低營養(yǎng)級餌料生物的捕食者, 又被許多高營養(yǎng)級肉食性魚類所捕食,在海洋生物群落和食物網(wǎng)中占有重要的地位[22—24]。

表3 膠州灣12種餌料魚類體長與耳石長度、寬度、長半徑和短半徑的回歸方程及其決定系數(shù)(R2)Tab.3 Regression equations and determination coefficient for fish standard length (SL) with otolith length (OL), width (OW), long radius (LR) and short radius (SR) for twelve prey fish species from Jiaozhou Bay
研究發(fā)現(xiàn), 膠州灣12種餌料魚類的體長與耳石長度、寬度、長半徑和短半徑之間均呈顯著的線性相關, 這與國內外很多學者的研究結果相似[10,14,16,25]。在所有 48個線性回歸方程中, 有 37個回歸方程的R2值均大于0.8, 擬合效果較好; 有9個回歸方程的R2值在0.7—0.8之間; 其余2個回歸方程的R2值小于0.7, 擬合效果略差。其中, 各方程擬合效果都最好的是六絲鈍尾蝦虎魚, 其R2值均接近于 1; 各方程擬合效果最差的是 鳀 , 其R2值均小于0.8。各個魚種的擬合效果存在差異, 可能與不同魚類及其耳石的生長速率和增長模式不同有關[16,24,26]。
通過比較同種魚類各個回歸方程的決定系數(shù)R2值可知, 在這12種餌料魚類中, 多數(shù)魚種體長與耳石長度和寬度的擬合效果要略好于體長與耳石半徑的擬合效果, 而且與耳石半徑相比, 耳石的長度和寬度更容易測量, 誤差也較小, 因此在實際應用中,可以通過測量耳石的長度和寬度來推算餌料魚類的體長等形態(tài)參數(shù)。另一方面, 由于多數(shù)魚類的耳石比較脆, 尤其是耳石的前端更容易破損, 難以準確測量耳石的長度, 因此, 與耳石長度相比, 體長與耳石寬度的回歸方程更加適合于實際應用, Pietro,et al.[14]也得出過相同的結論。
需要說明的一點是, 由于耳石在捕食者的胃中會經(jīng)歷一定程度的化學腐蝕或機械磨損[11,27], 因此通過測量胃含物中耳石的大小來推算魚體長度時可能會存在一定程度的偏差, 因此, 需要進行相應的校對以減小誤差。另外, 由于不同海域或不同種群的魚類, 其體長和耳石的生長模式會存在一定的差異[28], 因此在將本文的研究結論應用于其他海域時,需要考慮到上述因素的影響。
[1] Dou S Z.An introduction to fish otolith research: techniques and applications [J].Studia Marina Sinica, 2007, 48:93—113 [竇碩增.魚類的耳石信息分析及生活史重建—理論、方法與應用.海洋科學集刊, 2007, 48: 93-113]
[2] Petursdottir G, Begg G A, Marteinsdottir G.Discrimination between Icelandic cod (Gadus morhuaL.) populations from adjacent spawning areas based on otolith growth and shape[J].Fisheries Research, 2006, 80(2-3): 182—189
[3] Xiang J G, Xiang J, Wang X L,et al.Somatic and lapillus otolith ontogenetic growth and development in larval and juvenileElopichthys bambusarichardson in the Xiang River[J].Acta Hydrobiologica Sinica, 2011, 35(5): 817—822 [向建國, 向勁, 王星璐, 等.鳡湘江仔稚魚個體和耳石生長發(fā)育研究.水生生物學報, 2011, 35(5): 817—822]
[4] Zhang G H, Dan S G, Miao Z G,et al.Otolith morphology of six cyprinid species with the use in species and stock discriminations [J].Acta Hydrobiologica Sinica, 1999, 23(6):683—688 [張國華, 但勝國, 苗志國, 等.六種鯉科魚類耳石形態(tài)以及在種類和群體識別中的應用.水生生物學報,1999, 23(6): 683—688]
[5] Potier M, Ménard F, Benivary H D,et al.Length and weight estimates from diagnostic hard part structures of fish, crustacea and cephalopods forage species in the western Indian Ocean [J].Environmental Biology of Fishes, 2011, 92(3):413—423
[6] Fitch J E, Brownell R L.Fish otoliths in cetacean stomach and their importance in interpreting feeding habits [J].Journal of the Fisheries Research Board of Canada, 1968, 25(12):2561—2574
[7] Pierce G J, Boyle P R.A review of methods for diet analysis in piscivorous marine mammals [J].Oceanography and Marine Biology Annual Review, 1991, 29: 409—486
[8] Tollit D J, Steward M J, Thompson P M,et al.Species and size difference in the digestion of otoliths and beaks: placations for estimates of pinniped diet composition [J].Canadian Journal of Fisheries and Aquatic Sciences, 1997, 54(1):105—119
[9] Watanabe H, Kubodera T, Ichii T,et al.Feeding habits of neon flying squidOmmastraphes bartramiiin the transitional region of the central North Pacific [J].Marine Ecology Progress Series, 2004, 266: 173—184
[10] Gamboa D A.Otolith size versus weight and body-length relationships for eleven fish species of Baja California,Mexico [J].Fishery Bulletin, 1991, 89(4): 701—706
[11] Granadeiro J P, Silva M A.The use of otoliths and vertebrae in the identification and size-estimation of fish in predator-prey studies [J].Cybium, 2000, 24(4): 383—393
[12] Harvey J T, Loughlin T R, Perez M A,et al.Relationship Between Fish Size and Otolith Length for 63 Species of Fishes from the Eastern North Pacific Ocean [M].Seattle:NOAA/ National Marine Fisheries Service.2000, 150
[13] Waessle J A, Lasta C A, Favero M.Otolith morphology and body size relationships for juvenile Sciaenidae in the Río de la Plata estuary (35-36°S) [J].Scientia Marina, 2003, 67(2):233—240
[14] Battaglia P, Malara D, Romeo T,et al.Relationships between otolith size and fish size in some mesopelagic and bathypelagic species from the Mediterranean Sea (Strait of Messina, Italy) [J].Scientia Marina, 2010, 74(3): 605—612
[15] China State Bureau of Technical Supervision.Specifications for oceanographic survey—Part 6: Marine biological survey GB/T 12763.6-2007 [S].Beijing: China Standard Press.2007[中國國家技術監(jiān)督局.海洋調查規(guī)范第6部分: 海洋生物調查GB/T 12763.6-2007.北京: 中國標準出版社.2007]
[16] Ou Y J, Liao R, Li J E,et al.Otolith growth characteristics ofCollichthys lucidusfrom the Pearl River Estuary [J].Chinese Journal of Zoology, 2012, 47(1): 88—95 [區(qū)又君, 廖銳, 李加兒, 等.珠江口魚棘頭梅童耳石的生長特性.動物學雜志, 2012, 47(1): 88—95]
[17] Anderson R O, Neumann R M.Length, weight and associated structural indices [A].In: Murphy B R, Wills D (Eds.),Fisheries Techniques [C].Bethesda: American Fisheries Society.1996, 447—481
[18] Du R Q.Biostatistics [M].Beijing: China Higher Education Press.1997, 69—211 [杜榮騫.生物統(tǒng)計學.北京: 高等教育出版社.1997, 69—211]
[19] Liu R Y.Ecology and Living Resources of Jiaozhou Bay[M].Beijing: Science Press.1992, 2—3 [劉瑞玉.膠州灣生態(tài)學和生物資源.北京: 科學出版社.1992, 2—3]
[20] Xu B D, Zhang F, Mei C,et al.Characteristics of fish community structure in the central Jiaozhou Bay in spring and summer [J].Chinese Journal of Applied Ecology, 2010,21(6): 1558—1564 [徐賓鐸, 張帆, 梅春, 等.膠州灣中部海域春、夏季魚類群落結構特征.應用生態(tài)學報, 2010,21(6): 1558—1564]
[21] Mei C, Xu B D, Xue Y,et al.Fish community structure and species diversity during autumn and winter in the central waters of Jiaozhou Bay [J].Journal of Fishery Sciences of China, 2010, 17(1): 110—118 [梅春, 徐賓鐸, 薛瑩, 等.膠州灣中部海域秋、冬季魚類群落結構及其多樣性研究.中國水產科學, 2010, 17(1): 110—118]
[22] Wei S, Jiang W M.Study on food web of fishes in the Yellow Sea [J].Oceanologia et Limnologia Sinica, 1992, 23(2):182—192 [韋晟, 姜衛(wèi)民.黃海魚類食物網(wǎng)的研究.海洋與湖沼, 1992, 23(2): 182—192]
[23] Tang Q S.Strategies of research on marine food web and trophodynamics between high trophic levels [J].Marine Fisheries Research, 1999, 20(2): 1—11 [唐啟升.海洋食物網(wǎng)與高營養(yǎng)層次營養(yǎng)動力學研究策略.海洋水產研究,1999, 20(2): 1—11]
[24] Zhang B.Preliminary studies on marine food web and trophodynamics in China coastal seas [D].Thesis for Doctor of Science.Ocean University of China, Qingdao.2005 [張波.中國近海食物網(wǎng)及魚類營養(yǎng)動力學關鍵過程的初步研究.博士學位論文, 中國海洋大學, 青島.2005]
[25] Liu W, Zhan P R, Zhao C G,et al.Study on morphological development of fall chum salmon (Oncorhynchus keta) otoliths in Heilongjiang River [J].Acta Hydrobiologica Sinica,2010, 34(6): 1069—1076 [劉偉, 戰(zhàn)培榮, 趙春剛, 等.黑龍江秋大麻哈魚耳石形態(tài)發(fā)育研究.水生生物學報, 2010,34(6): 1069—1076]
[26] Xie Y H, Li B, Fu L J,et al.Daily growth increment in otolith and growth for fry—young fish of bighead (Aristichthys nobilis) [J].Journal of Fishery Sciences of China, 1995, 2(2):34—42 [解玉浩, 李勃, 富麗靜, 等.鳙仔-幼魚耳石日輪與生長的研究.中國水產科學, 1995, 2(2): 34—42]
[27] Jobling M, Breiby A.The use and abuse of fish otoliths in studies of feeding habits of marine piscivores [J].Sarsia,1986, 71(3—4): 265—274
[28] Campana S E, Casselman J M.Stock discrimination using otolith shape analysis [J].Canadian Journal of Fisheries and Aquatic Sciences, 1993, 50(5): 1062—1083