2.5 耳石數據及年齡成長分析
2.5.1 日輪判讀
早期生活史日齡的沉積,在過去許多鰻形目鰻鱺科的魚類中已被證實,無論 是野生溫帶種類 (Tsukamoto, 1989; Martin, 1995; Cieri and McCleave, 2001)、熱帶種 類 (Arai et al., 2000a; Sugeha et al., 2001),或是人工飼養實驗 (Umezawa et al., 1989;
Shinoda et al., 2004) 皆有報告,鰻魚初期生活史的日齡計算也已經普遍應用在不同 的物種上,因此本研究沿用此結果,假設小鰭鰭尾鯙早期生活史耳石之沉積亦是 一天形成一輪。
耳石日週輪的形成過程中,是由白晝時成長快速的連續帶與夜晚時成長較慢 的不連續帶所組成,在穿透光下,連續帶明亮,不連續帶則呈現深色,而一圈日 週輪的定義是包含一輪明亮的連續帶與一圈深色的不連續帶。計算從 FFC (first feeding check) 到不同記號位置的日輪數目 (圖九),TMC、TGC及TMG (dsys)。圖中,
MC (metamorphosis check) 代表輪寬從最窄處開始變寬前的最後一輪;GC (growth check) 代表 MC 後之一深刻記號,記號外之輪紋結構忽然以倍數急遽增加,甚至 模糊不清;而MG 則代表 MC 與 GC 之間輪寬微增加之區域。耳石日輪數目以 SAS 套裝軟體 (learning edition 4.1) 進行 one – way ANOVA 檢定是否有顯著差異,並以 Duncan's multiple range test 進行事後檢定。
2.5.2 日輪輪寬及平均耳石成長率計算
以影像分析軟體 (ImageJ®),統一沿最長軸 (anterior) 計算從 FFC 到 GC 間每 輪的輪寬 (increment width),同一標本以每三輪做一平均輪寬,並繪製成圖表以觀 察輪寬的時序變化。同時測量從FFC 到 MC、GC,以及 MG 的距離 (圖九),RMC、 RGC及 RMG (μm)。再以此距離除以相對應之日輪天數,便可得到平均耳石成長率 (average otolith growth rate),GMC、GGC及GMG,其單位為μm / day ( e.g., RGC / TGC
26
= GGC)。平均耳石成長率之統計檢定方法與日齡相同。
2.5.3 年輪判讀
雖然過去沒有關於小鰭鰭尾鯙年輪形成有效性驗證的文章,但張 (2009) 對臺 灣地區同為鯙科魚類的雲紋裸胸鯙 (Gymnothorax chilospilus) 進行耳石邊緣成長 率 (marginal increment ratio, MIR) 分析,發現年輪為一年形成一輪,且形成的高峰 位於10 – 1 月間; Matić-Skoko et al. (2011) 對地中海的澤生鯙 (Muraena helena) 進行MIR 分析,也得到了相似的結果。因此,本研究假設小鰭鰭尾鯙之年輪亦是 一年形成一輪,在顯微鏡下觀察5 % EDTA 腐蝕後的耳石,可見明顯之明帶與暗帶,
以暗帶之計算作為年齡的標準,並參考張 (2009)、ICES 2009 年對歐洲鰻鱺及美洲 鰻鱺年齡判讀之Workshop 及 Matić-Skoko et al. (2011) 之讀法,判讀小鰭鰭尾鯙之 年齡。在年齡的判讀上,以實際觀察到的耳石情況,配合捕捉月份、假設之出生 日與暗帶形成時間來估算年齡至月的精確度,詳細讀輪方法如 (圖二十) 所示。
2.5.4 體長體重關係式
將臺灣地區樣點小鰭鰭尾鯙測量所得之體長與體重資料,繪製成曲線圖並用 公式來描述之,如下:
BW = a TL
式中,a 與 b 是係數;BW 代表體重 (g);TL 代表全長 (cm)。
2.5.5 年齡長度關係
根據臺灣地區耳石年輪判讀的結果對應全長資料,以每 2 公分為一間距,建 構小鰭鰭尾鯙之年齡長度關係,用以了解臺灣地區小鰭鰭尾鯙的年齡組成結構。
27
2.5.6 成長方程式
將臺灣地區小鰭鰭尾鯙的年齡及全長帶入 von Bertalanffy growth equation (VBGE) 中,並以 R 進行非線性回歸得到另外三未知參數 k、t 0及L∞。VBGE 的 公式如下:
L = L 1 − e
式中,Lt為t 歲時的體長;L∞為理論之極限體長;k 為成長係數;t 0為體長為零時 之理論年齡。
2.5.7 讀輪精準度分析
所有耳石樣本在讀取前皆會以亂數重新排序,並給予新編號,年輪及日輪皆 判讀兩次以上,每次判讀至少間距兩個禮拜。若兩次判讀結果相符則採用結果;
若不相同,則按照前兩次的判讀結果,進行第三次判讀,並採用第三次的判讀結 果進行後續分析。
為了討論耳石讀輪的精準度,另計算平均百分誤差 (average percent error, APE;
Beamish and Fournier, 1981) 與變異係數 (coefficient of variation, CV; Campana, 2001) 加以評估,並繪製年齡別誤差圖 (age-bias plot; Campana, 2001)。APE 及 CV 的公式如下:
APE = 1 −
× 100%
CV =
∑ −
− 1 × 100%
式中,R 代表重複讀輪的次數;Xij代表在第 i 次判讀時,第 j 個樣本的輪數;Xj
代表j 個樣本在重複計算後之平均輪數。
28