金岳,方舟,葉旭昌、2、3,李鳳瑩
(1.上海海洋大學(xué) 海洋科學(xué)學(xué)院,上海201306;2.國家遠洋漁業(yè)工程技術(shù)研究中心,上海201306;3.上海海洋大學(xué) 大洋漁業(yè)資源可持續(xù)開發(fā)省部共建教育部重點實驗室,上海201306)
沙丁魚Sardina pilchardus 為鯡形目的典型物種,是重要的洄游性經(jīng)濟魚類,廣泛分布于西大西洋各海域及地中海[1-2]。作為近海暖水性魚類,沙丁魚一般不出現(xiàn)在外海和大洋。它們游速快,棲息于近海中上層,但秋冬季由于水溫較低會移動到較深海域[3]。根據(jù)沙丁魚在夜間易被亮光吸引的特性,捕撈者一般采用燈光誘捕法進行捕撈。摩洛哥地處非洲西北部,西臨大西洋,海岸線長達2446 km。由于該海域受加那利寒流和涌升的影響,營養(yǎng)物質(zhì)豐富,餌料充足,因此魚類種類繁多,其中以中上層的沙丁魚聞名于世[4]。據(jù)統(tǒng)計,2007年世界沙丁魚漁獲量為99 萬t,而摩洛哥沿岸的漁獲量為57 萬t,占全世界沙丁魚總漁獲量的57%,近幾年基本穩(wěn)定在60 萬t 左右,沙丁魚產(chǎn)量位居世界第一[5]。本研究中,依據(jù)2011年6—8月中國資源探捕漁船在摩洛哥專屬經(jīng)濟區(qū)海域采集的沙丁魚樣本,對其基礎(chǔ)生物學(xué)進行初步分析,旨在為今后更好地利用沙丁魚資源奠定基礎(chǔ)。
調(diào)查船為上海蒂爾遠洋漁業(yè)有限公司的“HUYU855 號”漁船,該船屬冷海水拖網(wǎng)漁船,總長為39.20 m,型寬為7.38 m,設(shè)計吃水為3.5 m,總噸位為413 t,主機功率為1491 kW。所使用的網(wǎng)具為四片式中層拖網(wǎng),為荷蘭MARITIEM 公司設(shè)計的Pelagic trawl 650 MTR 型漁具,網(wǎng)板采用6 m2立式V 型網(wǎng)板,總重約1300 kg,采用雙手綱連接。囊網(wǎng)最后部分采用PES 聚乙烯網(wǎng)線,其他部分采用不同粗度的尼龍編繩或編線,囊網(wǎng)由單層網(wǎng)衣組成,前部網(wǎng)衣網(wǎng)目內(nèi)徑為44 mm,末端網(wǎng)目尺寸為55 mm。
沙丁魚樣本的采集海域為摩洛哥專屬經(jīng)濟區(qū)海域(21°00' ~26°00'N、15°30' ~20°00'W),調(diào)查站位的設(shè)置如圖1所示。
圖1 摩洛哥專屬經(jīng)濟區(qū)沙丁魚采樣站點示意圖Fig.1 Sampling stations of sardine Sardina pilchardus at EEZ waters of Morocco
作業(yè)漁法為單拖作業(yè),拖曳時間根據(jù)漁場魚群和漁獲物數(shù)量而定,時間不定,可從幾十分鐘至6 ~7 h 不等,拖速基本維持在4 kn 以上。探捕時間為2011年6—8月,隨機取樣后共獲得2036尾樣本,冷凍后帶回實驗室進行測量與分析。
對沙丁魚生物學(xué)的測定內(nèi)容包括叉長、體質(zhì)量、性別、性腺成熟度等。叉長用量魚板測量,精確至1 mm。性別的確定主要基于對性腺的形狀、外觀和結(jié)構(gòu)的分析。根據(jù)黃海水產(chǎn)研究所改編的《海洋水產(chǎn)資源調(diào)資手冊》(1981年),將性腺成熟度劃分為Ⅰ~Ⅵ期,并分為性未成熟(Ⅰ~Ⅲ期)、性成熟(Ⅳ~Ⅴ期)和繁殖后(雄性為排精后,雌性為產(chǎn)卵后)(Ⅵ期)3 個等級。根據(jù)陳新軍[6]對魚類食性研究中所用的直接觀察法進行描述,將沙丁魚的胃級分為0 ~4 五個等級,其中0級為空胃,1級為胃含物不足1/2,2級為胃含物占胃的1/2,3級為胃含物超過胃的1/2,4級為滿胃。
(1)利用頻度分析法分析沙丁魚樣品的叉長和體質(zhì)量組成,組間距分別為20 mm和30 g。
(2)采用線性回歸法,求得雌雄叉長與體質(zhì)量之間的關(guān)系:
其中:W 為體質(zhì)量(g);L 為叉長(mm);a、b為估算參數(shù)。
(3)利用方差檢驗(ANOVA),對不同性腺成熟度個體的叉長和體質(zhì)量進行比較,差異特別顯著的,用最小顯著性差異法(LSD)進行組間多重比較,顯著性水平設(shè)為0.05。
(4)采用Logistic 曲線,推算不同性別的沙丁魚個體達初次性成熟時的叉長方程為
其中:Pi為成熟個體占組內(nèi)樣本的比例;Li為各組叉長(mm);c、d 為估算參數(shù)。
沙丁魚個體初次性成熟時的叉長L50%=-c/d。
對叉長分析的結(jié)果表明:摩洛哥海域沙丁魚樣本中,雌性個體的叉長范圍為150 ~270 mm,優(yōu)勢叉長為190 ~230 mm,占雌性總數(shù)的81.51%;雄性個體的叉長范圍為150 ~250 mm,優(yōu)勢叉長為190 ~230 mm,占雄性總數(shù)的84.47%(圖2)。雌性個體叉長在較小個體(<190 mm)和較大個體(>230 mm)范圍均有一定量的分布,而雄性個體叉長較為集中(190 ~230 mm)。雌性的叉長較雄性大,不同性別間存在顯著性差異(t=4.856,P=0.000<0.01)。
對體質(zhì)量分析的結(jié)果表明:沙丁魚雌性個體的體質(zhì)量范圍為60 ~210 g,優(yōu)勢體質(zhì)量為90 ~150 g,占雌性總數(shù)的85.46%;雄性個體的體質(zhì)量范圍為60 ~210 g,優(yōu)勢體質(zhì)量為90 ~150 g,占雄性總數(shù)的93.44%(圖2)。雌、雄個體的優(yōu)勢體質(zhì)量相差不大,但雄性優(yōu)勢體質(zhì)量的比例相對較高。
圖2 沙丁魚叉長和體質(zhì)量的頻率分布Fig.2 Distribution histogram of fork length and body weight in sardine Sardina pilchardus
通過對體質(zhì)量與叉長的對應(yīng)值進行分析,得到二者之間的相關(guān)性。由于雌、雄個體的叉長和體質(zhì)量均存在顯著性差異,為此,本研究中將雌、雄分開進行討論,結(jié)果如圖3所示。雌、雄沙丁魚的體質(zhì)量(W)與叉長(L)的關(guān)系方程分別為
W雌性=7×10-6L3.1076,R2=0.8483,P<0.001,
W雄性=9×10-6L3.0653,R2=0.8140,P<0.001。
圖3 沙丁魚體質(zhì)量與叉長的關(guān)系Fig.3 Relationship between fork length and body weight in sardine Sardina pilchardus
經(jīng)過對2036 尾沙丁魚樣本進行鑒別,樣本中雌性958 尾(占47.05%),雄性1026 尾(占50.39%),性別未定52 尾(占2.56%),雌雄性比為0.93∶1。從圖4可見:叉長≤170 mm和叉長>230 mm 的樣本中,雌雄性比均較大,分別為2.67∶1和3.3∶1,叉長為171 ~230 mm 的樣本中雌雄性比為(0.66 ~1.12)∶1。
統(tǒng)計分析發(fā)現(xiàn):雌性沙丁魚個體中未成熟個體、成熟個體、繁殖后個體分別占76.47%、19.64%、3.89%,其中,Ⅰ~Ⅵ期個體分別占2.21%、38.13%、36.13%、8.40%、11.24%和3.89%,大部分樣本都處在Ⅱ期和Ⅲ期,未成熟個體占樣本的大多數(shù);對于雄性個體,未成熟個體、成熟個體、繁殖后個體分別占69.04%、30.08%、0.88%,其中,Ⅱ~Ⅵ期個體分別占47.46%、21.58%、18.26%、11.82%和0.88%,繁殖后的個體非常少,未成熟個體所占比例比雌性大。無論雌性還是雄性,沙丁魚的叉長和體質(zhì)量均隨性腺成熟度的增大而增大(表1)。由于雌性Ⅰ期樣本較少,而雄性Ⅰ期沒有樣本以及Ⅵ期個體的樣本較少,因此,本研究中未對這些發(fā)育期的數(shù)據(jù)進行分析討論。
圖4 不同叉長范圍沙丁魚的雌雄性比組成變化Fig.4 Variation in sexual ratio in sardine Sardina pilchardus in different fork length groups
經(jīng)ANOVA 分析發(fā)現(xiàn),雌性沙丁魚不同性腺成熟度期的叉長(F=70.39,P=0.000<0.01)和體質(zhì)量(F=33.782,P=0.000<0.01)均存在極顯著性差異。經(jīng)LSD 分析表明,雌性Ⅱ期和Ⅲ期之間以及與其他等級的叉長和體質(zhì)量均存在極顯著性差異(P=0.000<0.01),Ⅳ期、Ⅴ期和Ⅵ期之間的叉長和體質(zhì)量均無顯著性差異(P>0.05)(表1)。
雄性沙丁魚不同性腺成熟度期的叉長(F=28.428,P=0.000<0.01)和體質(zhì)量(F=12.084,P=0.000<0.01)均存在極顯著性差異。經(jīng)LSD 分析表明,雄性Ⅱ期的叉長與其他等級均存在極顯著性差異(P=0.000<0.01),Ⅲ期與Ⅴ期的叉長存在極顯著性差異(P=0.000<0.01),而Ⅲ期與Ⅳ期以及Ⅳ期與Ⅴ期的叉長均無顯著性差異(P >0.05);雄性Ⅱ期的體質(zhì)量與Ⅳ期和Ⅴ期以及Ⅲ期與Ⅴ期的體質(zhì)量均存在極顯著性差異(P=0.001<0.01),Ⅱ期與Ⅲ期的體質(zhì)量存在顯著性差異(P=0.049<0.05),而Ⅲ期與Ⅳ期以及Ⅳ期與Ⅴ期的體質(zhì)量均無顯著性差異(P>0.05)(表1)。
對同一叉長組的雌性個體比較發(fā)現(xiàn)(圖5),在叉長≤170 mm 的樣本中,未出現(xiàn)Ⅳ期和Ⅵ期性腺,出現(xiàn)最多的是Ⅱ期性腺,Ⅰ期性腺次之;在叉長范圍為171 ~190 mm 的樣本中,各期性腺均有出現(xiàn),最多的是Ⅱ期性腺,Ⅲ期性腺次之;在叉長范圍為191 ~210 mm 的樣本中,Ⅰ期性腺僅有1尾,出現(xiàn)最多的是Ⅱ期性腺,Ⅲ期性腺次之;在叉長范圍為211 ~230 mm 的樣本中,未出現(xiàn)Ⅰ期性腺,出現(xiàn)最多的是Ⅲ期性腺,Ⅱ期性腺次之;在叉長范圍為231 ~250 mm 的樣本中,未出現(xiàn)Ⅰ期性腺,出現(xiàn)最多的是Ⅴ期性腺,Ⅳ期性腺次之。對同一叉長組的雄性個體比較發(fā)現(xiàn)(圖5),在叉長≤170 mm 的樣本中只出現(xiàn)了Ⅱ期性腺;在叉長范圍為171 ~190 mm 的樣本中,也未出現(xiàn)Ⅵ期性腺,出現(xiàn)最多的是Ⅱ期性腺,其他各期分布較均勻;在叉長范圍為191 ~210 mm 的樣本中,各期性腺都有出現(xiàn),最多的是Ⅱ期性腺,Ⅲ期性腺次之;在叉長范圍為211 ~230 mm 的樣本中,Ⅵ期性腺個體僅為6 尾,其他各期性腺數(shù)量基本一致;在叉長范圍為231 ~250 mm 的樣本中,Ⅱ、Ⅲ、Ⅳ、Ⅴ期性腺均有出現(xiàn),出現(xiàn)最多的是Ⅴ期性腺。
表1 不同性腺成熟度期叉長和體質(zhì)量的分布Tab.1 Distribution of fork length and body weight in sardine Sardina pilchardus with different maturity
圖5 不同叉長組沙丁魚性成熟度的組成Fig.5 The profile of maturity stage of sardine Sardina pilchardus in different fork length groups
利用Mann-Whitney U 檢驗發(fā)現(xiàn),雌雄之間性腺成熟度的分布存在顯著性差異(Z=2.594,P=0.000<0.05)。用列聯(lián)表法(contingency tables)對不同叉長組的雌雄個體性成熟度進行分析,結(jié)果表明:叉長≤170 mm(χ2=11.691,df=3,P=0.009<0.05)、叉長為171 ~190 mm(χ2=31.275,df=4,P=0.000<0.05)和叉長為211 ~230 mm(χ2=25.521,df=5,P=0.000<0.05)3 個叉長組雌雄個體性腺成熟度的分布均存在顯著性差異,而叉長為191 ~210 mm(χ2=7.77,df=5,P=0.169 >0.05)和231 ~250 mm(χ2=5.863,df=4,P=0.210>0.05)兩個組則不存在顯著性差異。
經(jīng)過Logistic 曲線的擬合,雌性和雄性的半數(shù)個體達初次性成熟時的叉長曲線如圖6所示,方程分別為
半數(shù)個體達初次性成熟時的叉長為L50%:雌性個體L50%=228.14 mm,雄性個體L50%=219.78 mm。
圖6 不同性別沙丁魚初次性成熟時的叉長(L50%)Fig.6 The fork length(L50%)at first maturity in females and males of sardine Sardina pilchardus
對所取樣本進行攝食等級分析,發(fā)現(xiàn)胃級為0級和1級的樣本最多,且所占比例相近(不分雌雄),而并未出現(xiàn)胃級為4級的樣本。對雌雄個體分別進行分析(圖7),發(fā)現(xiàn)雌性個體中,胃級為0級、1級、2級、3級的樣本分別占44.50%、42.03%、11.53%、1.94%,未出現(xiàn)胃級為4級的樣本;雄性個體中,胃級為0級、1級、2級、3級的樣本分別占46.76%、43.02%、9.21%、1.01%,也未出現(xiàn)胃級為4級的樣本。
圖7 沙丁魚攝食等級分析Fig.7 Analysis of food item proportion in sardine Sardina pilchardus
沙丁魚屬于較小個體的上層魚類,生長迅速,年齡最高可達到9 齡[7-8]。本研究中發(fā)現(xiàn),雌雄個體的優(yōu)勢叉長均在190 ~230 mm 范圍內(nèi),雌性個體的叉長大于雄性個體,不同性別之間的叉長分布和體質(zhì)量分布均存在顯著性差異。摩洛哥海域沙丁魚的生長在不同性別中存在著差異[9],這在對其他臨近海域不同群體沙丁魚的研究中也有所報道[10-12]。從叉長與體質(zhì)量的關(guān)系中也可以發(fā)現(xiàn),雌性個體的生長較快,這就更有利于沙丁魚的繁殖,對物種本身也是一種保護策略[5]。
摩洛哥海域沙丁魚的性別比例在不同的叉長范圍內(nèi)有著不同的變化。在較小的叉長范圍(<170 mm)和較大的叉長范圍(>231 mm)內(nèi),雌性比例明顯高于雄性,在其他叉長范圍內(nèi)雄性略多于雌性。Amenzoui等[5]對摩洛哥海域沙丁魚的繁殖生物學(xué)進行研究時也發(fā)現(xiàn),在較小的叉長范圍內(nèi),雌性個體較多,這可能是因為該地區(qū)為沙丁魚的繁殖場,該地聚集了大量性成熟的雌性沙丁魚所致。這種現(xiàn)象在摩洛哥北岸的大西洋海域[10]和阿爾沃蘭海海域[13]也有發(fā)現(xiàn)。Parrish等[14]根據(jù)地理分布和生物學(xué)特征,將北大西洋的沙丁魚分為4 個種群:分布于北海(57°N)至西班牙北岸(43°N)的北歐大西洋種群(Septentrional Atlantic stock),分布于葡萄牙沿岸至直布羅陀海峽(43°N ~36°N)的南歐大西洋種群(Meridional Atlantic group),分布于摩洛哥沿岸(36°N ~28°N)的摩洛哥種群(Moroccan stock),以及分布于毛里塔尼亞沿岸(西撒哈拉地區(qū)沿岸)(28°N ~21°N)的撒哈拉種群(Saharian stock)。從地理分布上看,本研究中的樣本屬于撒哈拉種群,屬于大個體。不同種群在不同時期的性別比例上也有一定的差異。
本研究中發(fā)現(xiàn),未成熟的雌雄個體均超過樣本的65%以上。沙丁魚全年都產(chǎn)卵,但是產(chǎn)卵的高峰期主要在每年的10月至翌年2月[5,15],本研究中的采樣時間并不在產(chǎn)卵高峰期內(nèi),因此也是樣本中未成熟個體較多的主要原因。叉長和體質(zhì)量都隨著性腺成熟度的上升而增加,從叉長和體質(zhì)量的增加趨勢來看,都可以發(fā)現(xiàn)未成熟個體(Ⅰ期、Ⅱ期和Ⅲ期)與成熟個體(Ⅳ期、Ⅴ期和Ⅵ期)之間存在顯著性差異,同時Ⅰ期、Ⅱ期和Ⅲ期之間也存在顯著性差異。這說明沙丁魚在未成熟時,生長迅速,各個時期均有一定的差異;而到性成熟階段,生長速度減緩,各個時期的差異不大。魚類在生長過程中,在初期生長很快,大量的能量僅供給魚體的生長;而在后期,由于大量的能量還要供給性腺的發(fā)育,因此個體生長的速度就相對變慢[16]。在對不同叉長的個體進行分析時發(fā)現(xiàn),雌雄之間的性成熟度分布存在顯著性差異,這與不同性別沙丁魚的性腺成熟時間不同有關(guān)。而在較小個體時(<190 mm),雌性需要迅速生長以便早于雄性,在叉長為211 ~230 mm 階段,雌性需要大量的能力供給卵巢生長,所以這兩個階段雌性的性腺成熟度分布與雄性有很大的不同。叉長為191 ~210 mm 時沙丁魚處在生長的中間期,叉長為231~250 mm 時大部分沙丁魚已經(jīng)達到性成熟期,因此,這兩個階段的雌雄差異并不大。另外,一些生態(tài)原因,如可利用的食物和溫度等都會對沙丁魚的性成熟產(chǎn)生影響[17-18]。
本研究中,對沙丁魚達到初次性成熟時的叉長進行分析,結(jié)果發(fā)現(xiàn),雌性和雄性的初次性成熟時的叉長比較接近,但是與其他學(xué)者的研究結(jié)果有所差異[5,11,19-20]。撒哈拉種群屬于較大群體,因此,初次性成熟時的叉長會較其他地區(qū)的種群大。而樣本中沒有發(fā)現(xiàn)滿胃個體,基本上都處于0 ~1級,這可能與未成熟個體較多和捕撈的時間有關(guān)。
沙丁魚主要是歐洲沿岸國家和北非國家的捕撈對象,中國首次派船隊對該海域進行探捕。本研究基于2011年6—8月中國探捕船采集的樣本,對沙丁魚的生物學(xué)特征進行了初步分析。中國學(xué)者之前也對中國近岸的沙丁魚類進行過研究[21-23],但目前尚未對其他海域不同種類的沙丁魚進行過調(diào)查研究,因此,在今后的研究中應(yīng)該更科學(xué)地確定樣本采集時間,并結(jié)合相關(guān)的海洋環(huán)境數(shù)據(jù),如溫度[24]、葉綠素a[25]等,盡量認定沙丁魚的年齡[26],以便對沙丁魚的年齡生長和生活史有更深入地了解。
[1]Furnestin J,F(xiàn)urnestin M L.La reproduction de la sardine et de l'anchois des c?tes atlantiques du Maroc(saisons et aires de ponte)[J].Rev Trav Inst Pêches Marit,1959,23(1):79-104.
[2]Ettahiri O.étude de la phase planctonique de la sardine,Sardina pilchardus(Walb.),et de l'anchois,Engraulis encrasicolus(L.)des c?tes atlantiques marocaines,thése[D].Brest:Universitéde Bretagne Occidentale,1996:262.
[3]Muzinic R.Migrations of adult sardines in the central Adriatic[J].Netherland Journal of Sea Research,1973,7:19-30.
[4]Sinovcic G.Stock size assessment of sardine,Sardina pilchardus(Walb.)population from the cenatral eastern Adriatic on the basis of VPA method[J].Acta Adriatica,1991,32(2):869-884.
[5]Amenzoui K,F(xiàn)erhan-Tachinante F,Yahyaoui A,et al.Analysis of the cycle of reproduction of Sardina pilchardus(Walbaum,1792)off the Moroccan Atlantic coast[J].C R Biologies,2006,329(10):892-901.
[6]陳新軍.漁業(yè)資源與漁場學(xué)[M].北京:海洋出版社,2004:27-32.
[7]Sinovcic G.Estimation of growth,mortality,production and stock size of sardine,Sardina pilchardus(Walb.),from the middle Adriatic[J].Acta Adriatica,1986,27(1/2):67-74.
[8]Sinovcic G.Responsible exploitation of the sardine Sardina pilchardus(Walb.),population in the coastal region of the Adriatic Sea[J].Periodicum Biologorum,2000,102(Suppl.1):47-54.
[9]Turner S C,Grimes C B,Able K W.Growth,mortality,and age/size structure of the fisheries for the tilefish,Lopholatilus chaemaelonticeps,in the middle Atlantic-Southern New England region[J].Fisheries Bulletin,1983,81:751-763.
[10]Belvèze H.Biologie et dynamique des populations de sardine(Sardina pilchardus)peuplant les c?tes atlantiques et proposition pour un aménagement des pêcheries,thèse d'état[D].Bretagne:Universitéde Brest Occidentale,1984:531.
[11]Abad R,Giraldez A.Reproduccion,factor de condicion y talla de primer madurez de la sardina,Sardina pilchardus(Walb.),del litoral de Malaga,mar de Alboran(1989-1992)[J].Bol Inst Esp Oceanogr,1993,9(1):145-155.
[12]Boely T.étude du cycle sexuel de la sardinelle ronde(Sardinella aurita Val.1847)au Sénégal,Cah.ORSTOM[J].Ser Oceanogr,1982,17(1):3-13.
[13]Garcia A,F(xiàn)ranco C,Sold A,et al.Sardina(Sardina pilchardus,Walb.)daily egg production of the Galicia,Cantabrian and Bay of Biscay Waters[J].ICES C M,1991,H(37):15.
[14]Parrish R H,Serra R,Grant W S.The monotypic sardines,Sardina and Sardinops:their taxonomy,distribution,stock structure,and zoogeography[J].Can J Fish Aquat Sci,1989,46(11):2019-2036.
[15]Muzinic R.Contributonà l' étude de l'oecologie de la sardina(Sardina pilchardus Walb.)dans l'Adriatique orientale[J].Acta Adriatica,1954,5(10):1-219.
[16]Sinovcic G,Kec V C,Zorica B.Population structure,size at maturity and condition of sardine,Sardina pilchardus(Walb.,1792),in the nursery ground of the eastern Adriatic Sea(Krka River Estuary,Croatia)[J].Estuarine,Coastal and Shelf Science,2008,76(4):739-744.
[17]Nikolsky G V.The Ecology of Fishes[M].London and New York:Academic Press,1963:352.
[18]Hempel G.On the importance of larval survival for the population dynamics of marine food fish[J].California Cooperative Oceanic Fisheries Investigations,1965,10:13-23.
[19]Villamil M M,Nespereira J M L,Pajuelo J M G,et al.Periodo reproductory madurez sexual de la sardina Sardina pilchardus(Walbaum,1792)en aguas de Gran Canaria(islas Canarias)[J].Bol Inst Esp Oceanogr,1997,13(1/2):47-55.
[20]Perez N,Porteiro C,Alvarez F.Contribucion al conocimiento de la biologia de la sardina de Galicia[J].Bol Inst Esp Oceanogr,1985,2(3):27-37.
[21]王啟華,姚善成,邱兆星.黃海區(qū)遠東擬沙丁魚的年齡與生長[J].海洋通報,1988,7(1):56-62.
[22]楊圣云,丘書院.臺灣海峽南部金色小沙丁魚的年齡和生長[J].臺灣海峽,1984,3(1):100-111.
[23]高天翔,楊天燕.中、日青鱗小沙丁魚形態(tài)學(xué)比較分析[J].中國海洋大學(xué)學(xué)報:自然科學(xué)版,2008,38(2):201-206.
[24]Ramzi A,Hbid M L,Ettahiri O.Larval dynamics and recruitment modeling of the Moroccan Atlantic coast sardine(Sardina pilchardus)[J].Ecological Modelling,2006,197(3/4):296-302.
[25]Ganias K.Linking sardine spawning dynamics to environmental variability[J].Estuarine,Coastal and Shelf Science,2009,84:402-408.
[26]Dulcic J.Estimation of age and growth of sardine,Sardina pilchardus(Walbaum,1792),larvae by reading daily otolith increments[J].Fisheries Research,1995,22(3):265-277.