Relationship between center of gravity and environmental factors of main catches of purse seine fisheries in North Pacific high seas based on GAM
-
摘要:
北太平洋公海灯光围网渔业是近年来我国新兴的远洋渔业项目，其中日本鲭 (Scomber japonicus)、远东拟沙丁鱼 (Sardinops sagax)、巴特柔鱼 (Ommastrephes bartrami) 和秋刀鱼 (Cololabis saira) 是主要渔获组成，厘清关键鱼种的渔场重心与环境因子的关系对研究渔场的形成机制及指导渔业生产具有重要意义。该研究利用2016—2017年北太平洋灯光围网渔业渔获数据及海表温 (Sea surface temperature, SST)、叶绿素 (Chl-a) 等环境因子，通过广义可加模型等分析渔获量的季节性变化及其与环境的关系。结果表明，北太平洋渔场重心整体上在渔汛初期向东北方向偏移，8—9月后期往西南折返。研究发现日本鲭渔场集中在148°E—154°E、40°N—42°N，远东拟沙丁鱼渔场集中在149°E—153°E、40°N—42°N，巴特柔鱼渔场集中在150°E—154°E、40°N—42°N，秋刀鱼渔场集中在150°E—153°E、41°N—42°N；2016—2017年北太平洋公海日本鲭的单位捕捞努力量渔获量 (Catch per unit effort, CPUE) 最适SST介于13~20.18 ℃，最适Chl-a介于0.1~0.8 mg·m−3；远东拟沙丁鱼CPUE最适SST介于12~19.1 ℃，最适Chl-a介于0.3~0.88 mg·m−3；巴特柔鱼 CPUE 最适 SST 介于 14~22 ℃，最适Chl-a介于0.2~1.1 mg·m−3；秋刀鱼CPUE 最适 SST 介于 14~21 ℃，最适Chl-a介于0.24~0.98 mg·m−3。
Abstract:Light seine fishery in North Pacific high seas is a new offshore fishery project in China in recent years. Scomber japonicus, Sardinops sagax, Ommastrephes bartramii and Cololabis saira are the main economic species in that area. It is important to clarify the relationship between the center of gravity and environmental factors for the study on the formation mechanism of fishery and the guidance of fishery production. Based on the 2016−2017 North Pacific Light Seine Fishery catch data and environmental factors such as sea surface temperature (SST) and chlorophyll a (Chl-a), we analyzed the seasonal variation of catch and its relationship with environment were analyzed by generalized additive model (GAM). On the whole, the center of gravity of North Pacific fisheries shifted northeast at the beginning of the fishing flood, and turned southwest during August to late September. The results show that S. japonicas fishing grounds were mainly at 148°E−154°E, 40°N−42°N; S. melanostictus fishing grounds were mainly at 149°E−153°E, 40°N−42°N; O. bartramii fishing grounds were mainly at 150°E−154°E, 40°N−42°N; C. saira fishing grounds were mainly at 150°E−153 °E, 41°N−42°N. Moreover, the optimal SST of S. japonicus CPUE was 13−20.18 ℃, and the optimal Chl-a was 0.1−0.8 mg·m−3. The optimal SST of S. sagax CPUE was 12−19.1 ℃, and the optimal Chl-a was 0.3−0.88 mg·m−3. The optimal SST of O. bartramii CPUE was 14−22 ℃, and the optimal Chl-a was 0.2−1.1 mg·m−3. The optimal SST of C. saira CPUE was 14−21 ℃, and the optimal Chl-a was 0.24−0.98 mg·m−3.
-
Keywords:
- Light seine /
- Commercial species /
- Environmental factors /
- GAM /
- North Pacific high seas
-
三疣梭子蟹 (Portunus trituberculatus) 隶属甲壳纲、十足目、梭子蟹科、梭子蟹属，广泛分布于西北太平洋及东南亚沿海，处于食物链较底层，生长速度快、适应能力强、分布范围广，是我国近海海洋生态系统中的优势物种，是海洋生态系统中物质传递和能量流动的重要载体，决定着海洋生态系统的许多重要生态过程。近年来，有关三疣梭子蟹的研究大多集中在生物学特性与养殖[1-5]、基因表达与免疫[6-9]等方面，对其年龄鉴定的相关研究甚少。但年龄鉴定一直是甲壳类生物研究的重要内容，有利于了解及掌握其种群结构、开展资源评估，为制订渔业管理措施提供依据，因此提供更加精确的直接年龄鉴定方法十分必要。
甲壳类生物因其特殊的生长方式 (蜕壳)，其年龄鉴定尤其困难。长期以来大多数学者采用标志重捕法[10]和体长频率法[11]等间接方法鉴定其年龄，但因方法的局限性，鉴龄的精确度不高；放射性同位素分析法[12]虽可较为精准鉴龄，但成本较高。最新研究表明甲壳类硬组织结构 (眼柄和胃磨) 能够作为直接鉴定其年龄的方法[13]。Leland等[14-15]利用甲壳动物的胃磨研磨截面微结构进行年龄鉴定，为甲壳动物年龄鉴定研究提供了新方向。Kilada等[13,16]先后发现美国螯龙虾 (Homarus americanus) 与斑纹黄道蟹 (Cancer irroratus) 两者的胃磨生长纹与年龄有着一定关系。现有研究中对三疣梭子蟹进行年龄鉴定的较少，蒋瑞等[17]发现三疣梭子蟹眼柄生长纹与其生长有关，但未肯定宽带个数代表其年龄。本研究通过对三疣梭子蟹胃磨微结构进行观察，探讨利用胃磨微结构对三疣梭子蟹进行年龄鉴定的可行性，分析胃磨微结构与基础生物学数据及年龄间的关系，以期为今后深入开展三疣梭子蟹年龄鉴定和其资源评估提供参考依据。
1. 材料与方法
1.1 实验样本的采集
三疣梭子蟹样本共185只，分别采集于三疣梭子蟹南通养殖基地 (100只)、江苏南通近海 (45只) 和上海芦潮港附近 (40只)。
1.2 胃磨的制片
测量取样个体生物学数据 [甲长 (Carapace length，CL)、甲宽 (Carapace width，CW)、体质量 (Body mass，m) ] 后取出胃磨，放置于70%的乙醇溶液中运送回实验室，在实验室中将胃磨转移至乙醇甘油混合溶液中 [V (乙醇)∶V (甘油)∶V (水)=6∶3∶1]；实验时使用解剖针剔除胃磨中附着的有机物，并使用蒸馏水冲洗干净。采用国内外胃磨结构划分的常规方法[18]，可将胃磨结构分为5个部分，即前幽门骨、外幽门骨、尾贲门骨、轭贲门骨和翼贲门骨 (图1)。以往的甲壳类鉴龄研究表明，外幽门骨和尾贲门骨经过预处理并切片后能在显微镜下清晰地观察到类似鱼类耳石或甲壳类眼柄中生长纹结构[15]，因此本研究采用这两种结构进行切片观察。
![]()
图 1 三疣梭子蟹胃磨结构展开示意图Figure 1. Schematic diagram of gastric mill structure of P. trituberculatus采用蒋瑞等[19]眼柄的制片方法，对于胃磨中外幽门骨和尾贲门骨在进行包埋前先使用无水乙醇脱水，再放入包埋盒中使用调配好的冷埋树脂溶液固定 (固化剂与亚克力粉按照体积比1.5∶1的比例调配，固化剂黏度为60 mPas)。包埋好的树脂样本放在阴凉通风处24 h待其硬化，后将包有胃磨的模具依次用120、240、600、1 200和2 500目的水磨砂纸打磨，在打磨过程中不断在显微镜下观察，直至能清晰观察到胃磨结构的平面，之后对此树脂薄片抛光，再清除其表面杂质，用吸水纸擦干表面水分。然后将抛光过的截面朝下，用热熔胶粘在载玻片上，使用同样的方法对另一面进行打磨，直至在显微镜下模具截面的生长纹清晰可见，再抛光得到切片。最后将做好的切片用密封袋密封保存，以备观察微结构时使用。
1.3 生长纹读取方法
生长纹的计数以胃磨中“暗纹”为计数标准。使用目视方法确定角质层边界，并且以该边界为计数起点。将胃磨切片照片依据照片的可读性分为3个等级。1) 清晰：沿着单个平面读数，可以连续且快速地读取所有生长纹的切片；2) 基本清晰：无法在单个平面内连续读取所有生长纹，但是通过平面的平移可以将生长纹全部读取的切片；3) 无法读取：无法辨别生长纹分界、明显脱钙、生长纹不完全及有明显缺陷的切片，这切片的读数不计入后续的数据分析。
每个样本的胃磨切片照片分别由3名不同实验室人员在相互独立的环境下判读其轮纹数量。采用变异系数 (CV) 以校对每个样本读数的准确性，公式为:
$$ {\rm{CV}} = \frac{{\sqrt {\frac{{{{\left( {{R_1} - R} \right)}^2} + {{\left( {{R_2} - R} \right)}^2} + {{\left( {{R_3} - R} \right)}^2}}}{2}} }}{R} \times 100{\text{%}} $$ (1) 其中R1、R2、R3分别表示3名不同人员对同一样本的读数，R表示3位人员读数的平均值，根据其他研究的经验[20-21]，当每个样本读数的变异系数值不高于10%时认为该次读数可信。
1.4 蜕壳次数的估算
通过周期性的蜕壳实现生长是甲壳动物最明显的生长特点。掌握甲壳类动物的蜕壳规律是研究其年龄与生长等信息的重要依据。高天龙等[22]认为三疣梭子蟹蜕壳周期时间随着期数的增加而增大，但其个体生命史中只存在有限数量的蜕壳次数，同时引入了单位阶梯函数来表示蜕壳次数f(m) 与体质量m的关系，具体函数为：
当m <0，u(m) =0；当m >0，u(m) =1；
$$ \begin{array}{c} f\left( m \right) = u\left( m \right) + u\left( {m - 0.014} \right) + u\left( {m - 0.033} \right) + u( {m }-\\ 0.104 ) + u\left( {m - 0.411} \right) + u\left( {m - 1.4} \right) + u\left( {m - 3.89} \right) +\\ u\left( {m - 11.56} \right) +u\left( {m - 32.77} \right) + u\left( {m - 66.28} \right) + \\ u\left( {m - 114.02} \right) + u\left( {m - 177.84} \right) + u\left( {m - 256.61} \right) \end{array} $$ (2) 该阶梯函数对于蜕壳阶段在5—11期的三疣梭子蟹判对率高达91.38%，因此本文采用该阶梯函数来推算取样三疣梭子蟹的蜕壳次数，并以此来探索蜕壳次数与胃磨生长纹数量之间的关系。
1.5 生长方程模型
由于缺乏对其的精准鉴龄，有关其生长方程的研究也较少。故本研究采用目前推算水生生物常用的7种生长方程模型：线性模型、幂函数模型、指数模型、对数模型、Logistic模型、von Bertalanffy模型 (v-B模型) 和Gompertz模型来拟合三疣梭子蟹的生长方程，并开展最适模型的选择，各模型公式为：
线性方程，
$$ f\left( t \right) = {\rm{a}} + {\rm{b}}t $$ (3) 幂函数方程，
$$ f\left( t \right) = {\rm{a}}{t^{\rm{b}}} $$ (4) 指数方程，
$$ f\left( t \right) = {\rm{a}}{{\rm{e}}^{{\rm{b}}t}} $$ (5) 对数函数方程，
$$ f\left( t \right) = {\rm{a}}\ln t + {\rm{b}} $$ (6) Logistic生长方程，
$$ f\left( t \right) = \frac{{{{{L}}_\infty }}}{{1 + {{\rm{e}}^{ - {\rm{K}}\left( {t - {{\rm{t}}_0}} \right)}}}} $$ (7) von Bertalanffy生长方程，
$$ f\left( t \right) = {L_\infty } \times \left[ {1 - {{\rm{e}}^{ - {\rm{K}}\left( {t - {t_0}} \right)}}} \right] $$ (8) Gompertz生长方程，
$$ f\left( t \right) = {L_\infty } \times {{\rm{e}}^{\left[ {1 - {{\rm{e}}^{ - {\rm{K}}\left( {t - {t_0}} \right)}}} \right]}} $$ (9) 其中f (t) 表示三疣梭子蟹的甲长 (CL)、甲宽 (CW) 或体质量 (m)，t表示胃磨生长纹个数，a、b、K、L∞和t0表示对应公式中常数项。对于Logistic方程、v-B方程和Gompertz方程，L∞表示无限年龄的三疣梭子蟹渐进的甲长、甲宽或体质量，t0表示理论上当甲长、甲宽或体质量为0时三疣梭子蟹的生长纹数量，在一般情况下，该值为负值。K值则表示常数。为了更加快速地拟合各公式，本研究中在迭代过程中起始值默认为所有样本中甲长、甲宽或体质量的最大值，而为了避免出现理论上的错误，t0和K的起始值则为0，并约束−t0和K值为正值。
在选取最适函数时，使用赤池信息量准则(Akaike information criterion, AIC)，AIC为最小时，则认为该模型在常用7个生长模型当中为最适模型，选定最适函数后使用回归系数 (R2) 作为判断数据对该模型的拟合程度。
AIC计算方法为：
$$ {\rm{AIC}} = 2k + n{\rm{ln}}\left( {\frac{{{\rm{RSS}}}}{n}} \right) $$ (10) 式中k为模型的参数数量，n为样本数量，RSS为所有样本量的残差平方和。
2. 结果
2.1 胃磨的微结构
胃磨位于虾蟹类的消化系统中，大部分虾蟹类食道通入一个能够收缩的三角形囊状胃，胃内部角质层发达，与胃壁共同形成的嵴构成胃磨。由于胃磨形状的不规则性及胃磨由多块相对独立的部分组合而成，因此将胃磨在显微镜下所观察的微结构分为两部分：1) 在肉眼下能观察到类似与动物牙齿的齿板部分；2) 不包含齿板的剩余白色部分 (柄部及翼部)。生长纹主要分布于胃磨各部分的柄部和翼端，本研究中胃磨微结构特指能够观察到生长纹的柄部切面微结构。
通过显微镜下观察三疣梭子蟹胃磨微结构后发现 (图2)，胃磨微结构由4个部分组成：蜡质层、外角质层、钙化层和膜层。蜡质层紧邻外角质层，外角质层位于胃磨靠向外一侧，边缘规则且光滑，但因内部含有少量色素透光度低，但随着切片厚度减小，透光度有所提升，主要由几丁质构成，在显微镜下观察不到任何的周期性生长纹结构。膜层在胃磨内侧与肌肉或其他组织接触部分，在显微镜下可以观察到黑色部分，透光度较高，钙化层是胃磨微结构中占比最大的一部分，钙化程度较高，在显微镜下能观察明暗相间的生长纹结构，与膜层为平移关系，同时钙化层中还存在垂直生长纹。
![]()
图 2 外幽门骨 (a)、尾贲门骨 (b) 和眼柄 (c) 微结构及生长纹Figure 2. Microstructure and growth bands of explyloric ossicle (a), urocardiac ossicl (b) and eyestalk (c)与同样作为甲壳类直接年龄鉴定材料的眼柄微结构 (图2-c) 相比，胃磨微结构与其略有不同。眼柄切片中色素层充满了红色色素，而胃磨中表层较薄，色素层呈现黑色，且色素层占比比眼柄切片小。胃磨钙化层所占比例要大于眼柄中钙化层。在胃磨和眼柄切片中均可观察到密集且间距均匀的细纹，还可以观察到颜色较深、排列稀疏的宽带，一个宽带内可以看到包含数条或十数条细纹。
2.2 生长纹读数的精确性
生长纹读数的准确度会直接影响到最终的实验结果，因此在利用生长纹个数拟合生长方程之前，要判断对样本生长纹读数的精确度。Kilada等[23]利用眼柄和胃磨直接鉴定南极磷虾 (Euphausis superba) 年龄发现，当变异系数CV小于10%，该组数据可被使用。在本研究中，由3名实验人员单独对同一样本的胃磨生长纹读数，经计算得出平均变异系数为8.6% (n=185)，因此该生长纹数据准确可靠，可用于拟合有关生长方程的研究。
2.3 生长模型建立
2.3.1 甲长、甲宽与体质量生长模型
通过基本生物学测量分析后，185只样本中包含雌蟹120只，雄蟹65只。总体样本中甲长为20~99 mm [平均值 (57.27±1.30) mm]，甲宽介于43~190 mm [平均值 (116.25±2.22) mm]，体质量介于3.70~365.89 g [平均值 (105.86±5.92) g]。雌雄三疣梭子蟹基础生物学数据见表1，在平均值上雌蟹数据要略高于雄蟹。分别对所有的三疣梭子蟹甲长、甲宽和体质量数据进行分析，得出以下三疣梭子蟹甲长CL与甲宽CW的关系 (图3、图4、图5)。
表 1 三疣梭子蟹样本个体基础生物学数据Table 1. Basic biological data of gastric abrasion growth samples of P. trituberculatus项目
Item雄蟹 Male 雌蟹 Female 甲长 CL/mm 甲宽 CW/mm 体质量 m/g 甲长 CL/mm 甲宽 CW/mm 体质量 m/g 样本量 Sample number 65 65 65 120 120 120 最大值 Maximum 99 190 365.89 97 189 363.17 最小值 Minimum 20 43 3.7 24 60 7.29 平均值 Mean 50.05 105.47 85.33 61.18 122.08 166.97 标准误 Standard error 2.27 4.17 10.16 1.47 2.67 7.10 ![]()
图 3 雄性 (a) 与雌性 (b) 三疣梭子蟹甲长与甲宽关系图Figure 3. Relationship between carapace length and width of male (a) and female (b) P. trituberculatus![]()
图 4 雄性 (a) 与雌性 (b) 三疣梭子蟹甲宽与体质量关系图Figure 4. Relationship between carapace width and body mass of male (a) and female (b) P. trituberculatus![]()
图 5 185只三疣梭子蟹甲长与甲宽、甲宽与体质量关系图Figure 5. Relationship between carapace length and carapace width, carapace width and body mass of 185 individuals of P. trituberculatus雄性三疣梭子蟹：
$$ {\rm{CW}} = 1.77{\rm{CL}} + 16.94\;(n = 65,{R^2} = 0.93,P < 0.01) $$ (11) 雌性三疣梭子蟹：
$$ {\rm{CW}} = 1.72{\rm{CL}} + 16.93\;(n = 120,{R^2} = 0.89,P < 0.01) $$ (12) 185只三疣梭子蟹：
$$ {\rm{CW}} = 1.72{\rm{CL}} + 17.93\;(n = 185,{R^2} = 0.89,P < 0.01) $$ (13) 三疣梭子蟹甲宽CW与体质量m的关系：
雄性三疣梭子蟹：
$$ m = 8.28 \times {10^{ - 5}}{\rm{C}}{{\rm{W}}^{2.91}}(n = 65,{R^2} = 0.87,P < 0.01) $$ (14) 雌性三疣梭子蟹：
$$ m = 4.92 \times {10^{ - 5}}{\rm{C}}{{\rm{W}}^{3.02}}(n = 120,{R^2} = 0.92,P < 0.01) $$ (15) 185只三疣梭子蟹：
$$ m = 6.31 \times {10^{ - 5}}{\rm{C}}{{\rm{W}}^{2.97}}(n = 185,{R^2} = 0.92,P < 0.01) $$ (16) 2.3.2 甲长、甲宽、体质量与生长纹的关系
根据最小AIC法原则 (表2)，185个样本中，甲长、甲宽与胃磨生长纹个数的关系均适合v-B生长模型来表示 (图6)，关系式分别为：
表 2 7个生长模型函数拟合结果Table 2. Seven growth model function fitting results模型
Model甲长-生长纹个数
CL-t甲宽-生长纹个数
CW-t体质量-生长纹个数
m-tR2 AIC R2 AIC R2 AIC 线性 Linear 0.884 473.77 0.922 539.78 0.868 710.89 幂函数 Power 0.834 504.01 0.886 571.06 0.910 678.42 指数 Exponent 0.865 486.36 0.906 555.32 0.913 675.90 对数 Logarithm 0.841 500.07 0.897 562.54 0.698 780.75 逻辑斯蒂 Logistic 0.897 465.69 0.953 499.09 0.907 683.98 v-B 0.910 449.19 0.995 4 498.26 0.914 677.15 冈珀茨 Gompertz 0.870 484.63 0.949 504.56 0.922 668.93 注：加粗表示拟合结果最好 Note: Bold indicates the best fit result. ![]()
图 6 三疣梭子蟹甲长、甲宽、体质量与生长纹数目的关系Figure 6. Relationship between carapace length, carapace width, body mass and growth increments of P. trituberculatus$$ {\rm{CL}}=103.42 \times ( {1-{{\rm{e}}^{ - 0.027\;7t-1.149}}} )( {n=185,{R^2}=0.91,P < 0.01} ) $$ (17) $$ {\rm{CW}} = 202.96 \times ( {1 - {{\rm{e}}^{ - 0.032\;2t}}} )(n = 185,{R^2} = 0.90,P < 0.01) $$ (18) 体质量与胃磨生长纹个数的关系适合使用Gompertz模型来表示 (图6)，且关系式为：
$$ m = 356.27 \times {{\rm{e}}^{\left[ {1 - {{\rm{e}}^{ - 0.027\;1\left( {t - 59.6} \right)}}} \right]}}\left( {n = 185,{R^2} = 0.92,P < 0.01} \right) $$ (19) 同样，雄性三疣梭子蟹甲长、甲宽、体质量与其胃磨生长纹数目的关系均适于用Gompertz模型描述，表达式分别为：
$$ {\rm{CL}} = 51.296 \times {{\rm{e}}^{\left[ {1 - {{\rm{e}}^{0.028\left( {t - 26.014} \right)}}} \right]}}( {n = 65,{R^2} = 0.98,{{P}} < 0.01} ) $$ (20) $$ {\rm{CW}} = 87.853 \times {{\rm{e}}^{\left[ {1 - {{\rm{e}}^{ - 0.033\left( {t - 18.471} \right)}}} \right]}}( {n = 65,{R^2} = 0.99,{{P }}< 0.01} ) $$ (21) $$ m = 2{\rm{}}512.410 \times {{\rm{e}}^{\left[ {1 - {{\rm{e}}^{0.015\left( {t - 129.841} \right)}}} \right]}}( {n = 65,{R^2} = 0.96,{{P}}< 0.01} ) $$ (22) 对于雌性三疣梭子蟹，其甲宽以及体质量与胃磨生长纹数目的关系均适于用幂函数生长模型来描述，表达式分别为：
$$ {\rm{CW}} = 17.704{t^{0.573}}(n = 120,{R^2} = 0.89,P < 0.01) $$ (23) $$ m = 0.357{t^{1.684}}(n = 120,{R^2} = 0.81,P < 0.01) $$ (24) 雌性甲长与胃磨生长纹数目的关系则适于用v-B生长模型描述，表达式为：
$$ {\rm{CL}}=98.007\times\left[ {1-{{\rm{e}}^{ - 0.037{\rm{}}\left( {t - 1.201} \right){\rm{}}}}} \right](n=120,{R^2}=0.84,P < 0.01) $$ (25) 2.4 生长纹与预估蜕壳次数关系
在本研究中选取70个胃磨切片清晰的样本，利用单位阶梯函数的判定，得到其样本的预估蜕壳次数。通过对胃磨生长纹和预估蜕壳次数的拟合，两者存在幂函数关系 (图7)，关系式为u=0.207 1t2.105 2(n=70，R2=0.77)。
![]()
图 7 三疣梭子蟹预估蜕壳次数与胃磨生长纹关系拟合图Figure 7. Relationship between estimate molting preriod and growth increments of P. trituberculatus3. 讨论
Leland等[14]在研究红鳌鳌虾 (Cherax quadricarinatus) 胃磨时发现胃磨研磨截面展示出与鱼类耳石相近的同心生长纹，这开启了虾蟹类年龄鉴定研究的新方向，而后又使用钙黄绿素染色追踪红鳌鳌虾胃磨生长后，在胃磨切片中能明显看到生长痕迹，这也是目前证明胃磨能够作为甲壳类直接年龄鉴定最有力的证据[15]。Sheridan等[24]和Kilada等[13]在实验中采用外幽门骨或尾贲门骨作为甲壳类年龄鉴定，结果显示两者组织切片在显微镜下能够观察到清晰的生长纹，并且其生长纹与年龄信息显著相关。Sheridan等[24]利用硬组织生长纹计数法对普通黄道蟹 (Cancer pagurus) 鉴龄时，其轭贲门骨中生长纹清晰可见，但未在眼柄中观察到清晰的生长纹，这不排除是样本制备所导致。本研究结果显示，胃磨微结构中存在与眼柄相似的生长纹，且胃磨微结构中钙化层比例较大，生长纹清晰可见。对于甲壳类动物来说，目前为止没有哪一种方法可以独立对其进行精确鉴龄，因此可以同时对眼柄及胃磨生长纹分析，两者互为辅助，并与多种鉴龄方法相结合，交叉验证，以提高鉴龄的精确性。对甲壳类年龄鉴定最大的问题就是所使用的材料是否会在甲壳类独特的生长过程 (周期性蜕壳) 中被蜕去。如果在蜕壳过程中所使用的年龄鉴定材料随着老化的外壳丢失，那么该鉴龄材料则只能记录两次蜕壳间的生长信息。但是，Kilada等[13]在对雪蟹 (Chionoecetes opilio) 蜕壳的观察中发现，其眼柄不会随着蜕壳过程丢失，这使得眼柄可以成为甲壳类鉴龄材料。同样地，Leland等[15]也证明了在水晶蟹 (Chaceon albus) 蜕壳过程中，其胃磨也不会随着蜕壳过程而丢失。Krafft等[25]研究发现，胃磨当中的生长纹与已知能够作为甲壳类年龄鉴定材料的眼柄生长纹数量呈函数关系。因此胃磨中生长纹可以作为年龄鉴定材料。三疣梭子蟹通过蜕壳实现生长，在整个生长周期内，蜕壳次数有限。王小刚[26]认为三疣梭子蟹的龄期 (每蜕壳1次，称作1个龄期) 与其蜕壳次数呈正相关关系。故本研究引入高天龙等[22]建立的通过体质量判断期龄的判定函数，以此来探索蜕壳次数与胃磨生长纹数量之间的关系。结果发现，蜕壳次数与生长纹数目之间呈幂函数关系。即在幼蟹时期，三疣梭子蟹每次蜕壳间胃磨所累积的生长纹数目要远少于成年后三疣梭子蟹每次蜕壳间胃磨所累积的生长纹数目。
国内外渔业生物学研究者在对甲壳类生物研究中，各甲壳类甲宽 (头胸甲长) 与体质量关系的研究均采用幂函数，其中幂函数系数b均介于2.5~3.5。通过对所采集的样本分析，本研究中总样本的幂指数系数为3.05，雄蟹为2.91，雌蟹为3.02。即可认为本研究所使用样本符合大部分三疣梭子蟹的生长规律，即等速生长，说明三疣梭子蟹个体在整体生活史中体型随着质量增加等比增大。同时本研究中雌雄三疣梭子蟹幂指数系数相差不大，说明三疣梭子蟹甲宽与体质量关系函数在性别上没有显著差异。与王小刚等[27]使用v-B生长模型所计算得出的渐近甲宽 (L∞=210 mm) 相比，本研究中通过胃磨生长纹代替实际年龄所拟合出的函数其渐近甲宽 (L∞=202.96 mm) 要小。同时从以往的研究来看，在20世纪后期林景祺[28]在中国海洋渔业资源中指出中国近海三疣梭子蟹甲宽范围在30~250 mm，优势甲宽为120~190 mm；张健等[29]研究表明，在中国东海近海的三疣梭子蟹甲宽范围在63.3~201.5 mm，优势甲宽仅为110~160 mm。说明中国近海部分三疣梭子蟹偏向于小型化，且相对于早期个体有明显差距，这也是本研究中渐近甲宽、甲长和体质量要略小于早几年研究的原因。三疣梭子蟹甲宽与甲长之间为线性关系，这也和大多数学者研究结果一致。此外，逻辑斯蒂生长方程、v-B生长模型和Gompertz方程均可用于描述水生生物生长，而在本研究中对于三疣梭子蟹胃磨生长纹个数与其基础生物学数据的关系中，甲长与甲宽和胃磨生长纹数目关系适于用v-B生长模型描述，而体质量和胃磨生长纹数目则更适于用Gompertz模型进行描述。Kilada等[30]在雌性欧洲龙虾 (Homarus gammarus) 和雌性斑纹黄道蟹 (Cancer irroratus) 的轭贲门骨切片中分别观察到了4个和9个宽生长纹，且证实宽带个数与两者实际年龄相近。Kilada和Acuna[16]也在黄蹲龙虾 (Curvimunida johni) 的中位骨切片中观察到了宽带，且认为此宽带是以年为单位形成的。Ismail等[31]研究得出榄绿青蟹 (Scylla oliviacea) 轭贲门骨切片中的宽带数目可以作为估计其年龄的依据。本研究在三疣梭子蟹外贲门骨切片中也观察到了宽带，依据前人研究结果，推测可以用宽带数目估计其年龄。
本研究表明，采用胃磨对三疣梭子蟹鉴龄是可行的，但本实验采用的均为未知年龄的三疣梭子蟹样本，因此只能依据前人研究成果推测胃磨切片中宽带个数与实际年龄相近，后续应使用已知年龄的样本对该推测进行验证。在探讨三疣梭子蟹胃磨生长纹与其年龄的相关性时，笔者尚无法对胃磨生长纹的形成机制作出合理解释。未来的研究重点应落在以下两个方向：1) 延长实验样本 (已知年龄) 采集周期并缩短采集间隔，以验证更加精确的周期性实验，提高三疣梭子蟹胃磨鉴龄的精确性；2) 采用材料学方法 (电子显微镜、荧光染色和同位素分析等) 探究三疣梭子蟹胃磨的形成机制。
-
![]()
图 1 北太平洋公海灯光围网作业区域
Figure 1. Light purse seine fishing area in North Pacific high seas
![]()
图 2 2016—2017年北太平洋公海主要渔获物月平均产量和捕捞努力量渔获量
Figure 2. Average monthly catches and CPUE of main species in North Pacific high seas during 2016−2017
![]()
图 3 2016—2017年北太平洋公海主要渔获品种的渔场重心季节性变化
左上方灰线所围区域为专属经济区
Figure 3. Seasonal variation of fishing ground gravity center of main species catches in North Pacific during 2016−2017
The area encircled by gray line on the upper left is Exclusive Economic Zone.
![]()
图 4 2016—2017年各月份海表温度与渔获量分布叠
Figure 4. Stack distribution of sea surface temperature and CPUE during 2016−2017
![]()
图 5 2016—2017年各月份叶绿素a浓度与渔获量分布叠加
Figure 5. Stack distribution of Chl-a concentration and CPUE during 2016−2017
![]()
图 6 2016—2017年北太平洋公海灯光围网渔场主要品种渔获量与环境因子的GAM模拟结果
Figure 6. GAM simulation results of main species cathes and environmental factors in fishing grounds of North Pacific high seas during 2016−2017
表 1 2016—2017年北太平洋渔场环境因子与渔获量GAM模拟的检验系数值和显著性
Table 1 Inspection coefficient value and significance for GAM simulation of environmental factors and catches in North Pacific fishing grounds during 2016−2017
鱼种
Fish species环境因子
Environment factor有效自由度
Effective degrees of freedomF P 日本鲭 S. japonicas 经度 Lon 7.12 8.76 1.48×10−11 *** 纬度 Lat 5.47 8.16 5.54×10−9 *** 海表温 SST 3.43 7.23 5.20×10−6 *** 叶绿素a Chl-a 1.00 11.65 6.79×10−4 *** 沙丁鱼 S. sagax 经度 Lon 1.00 0.99 0.32 纬度 Lat 1.00 8.34 4.11×10−3 ** 海表温 SST 3.65 2.53 3.46×10−2 * 叶绿素a Chl-a 1.00 1.22 0.27 柔鱼 O. bartrami 经度 Lon 4.59 3.88 1.52×10−3 ** 纬度 Lat 1.00 12.80 4.38×10−4 *** 海表温 SST 2.74 4.42 2.93×10−3 ** 叶绿素a Chl-a 3.13 5.18 5.88×10−4 *** 秋刀鱼 C. saira 经度 Lon 1.00 11.45 9.54×10−4 *** 纬度 Lat 1.00 16.73 7.60×10−5 *** 海表温 SST 1.00 9.93 2.12×10−3 ** 叶绿素a Chl-a 1.00 1.24 0.27 注：*. P<0.05; **. P<0.01; ***. P<0.001 
下载: 导出CSV
-
[1] FAO. State of world fisheries and aquaculture (SOFIA) 2018. [EB/OL]. (2018-06-01) [2019-10-1]. http://www.fao.org/fisheries/en/
[2] YU W, GUO A, ZHANG Y, et al. Climate-induced habitat suitability variations of chub mackerel Scomber japonicus in the East China Sea[J]. Fish Res, 2018, 207: 63-73. doi: 10.1016/j.fishres.2018.06.007
[3] 李纲, 陈新军. 东海鲭鱼资源和渔场时空分布特征的研究[J]. 中国海洋大学学报 (自然科学版), 2007, 37(6): 921-926. [4] 程家骅, 林龙山. 东海区鲭鱼生物学特征及其渔业现状的分析研究[J]. 海洋渔业, 2004, 26(2): 73-78. doi: 10.3969/j.issn.1004-2490.2004.02.001 [5] YATSU A, WATANABE T, ISHIDA M, et al. Environmental effects on recruitment and productivity of Japanese sardine Sardinops melanostictus and chub mackerel Scomber japonicus with recommendations for management[J]. Fish Oceanogr, 2005, 14(4): 263-278. doi: 10.1111/j.1365-2419.2005.00335.x
[6] KAMIMURA Y, TAKAHASHI M, YAMASHITA N, et al. Larval and juvenile growth of chub mackerel Scomber japonicus in relation to recruitment in the western North Pacific[J]. Fish Sci, 2015, 81(3): 505-513. doi: 10.1007/s12562-015-0869-4
[7] 由上龍嗣, 西嶋翔太, 井須小羊子ら. 平成30 (2018) 年度マサバ太平洋系群の資源評価[R]. 中央水産研究所, 2018: 4-8. [8] 渡邊千夏子. マサバ太平洋系群の繁殖特性の変化とその個体群動態への影響[J]. 水産海洋研究, 2010, 74: 46-50. [9] 西田宏, 川端淳, 目黒清美, ら. マサバとゴマサバの分布と回遊1−幼魚[J]. 水産海洋研究, 2001, 65(4): 201. [10] 目黒清美, 梨田一也, 三谷卓美. マサバとゴマサバの分布と回遊 (成魚)(総特集 マサバとゴマサバ太平洋系群−−漁業·資源·管理の現状と将来展望)[J]. 海洋, 2002, 34(4): 256-260. [11] 川崎健. マサバ太平洋系群未成魚の生態について[J]. 東海区水産研究所研究報告, 1968(55): 59-114. [12] 飯塚景記. 東北海区におけるマサバ未成魚の生態−3−八戸沖で越冬すると推定される魚群について[J]. 東北区水産研究所研究報告, 1974(33): 37-49. [13] WATANABE C, YATSU A. Long-term changes in maturity at age of chub mackerel (Scomber japonicus) in relation to population declines in the waters off northeastern Japan[J]. Fish Res, 2006, 78(2-3): 323-332. doi: 10.1016/j.fishres.2006.01.001
[14] 黒田一紀. 日本の太平洋沿岸域におけるさば属魚類の産卵期, 産卵場及び産卵量水準の動向[J]. 水産海洋研究, 1992, 56: 65-72. [15] 加藤充宏, 渡邊千夏子. マサバとゴマサバの成熟·産卵および食性 (総特集 マサバとゴマサバ太平洋系群——漁業·資源·管理の現状と将来展望)[J]. 海洋, 2002, 34(4): 266-272. [16] 庄之栋, 马超, 刘勇, 等. 北太平洋公海秋季鲐鱼生物学特性的初步研究[J]. 海洋湖沼通报, 2018(06): 125-131. [17] 薛嘉伦, 樊伟, 唐峰华, 等. 基于最大熵模型预测西北太平洋平洋公海鲐潜在栖息地分布[J]. 南方水产科学, 2018, 14(1): 92-98. doi: 10.3969/j.issn.20950780.2018.01.012 [18] 田思泉, 陈新军. 不同名义CPUE计算法对 CPUE 标准化的影响[J]. 上海海洋大学学报, 2010, 19(2): 240-245. [19] 唐峰华, 崔雪森, 杨胜龙, 等. 海洋环境对中西太平洋金枪鱼围网渔场影响的 GIS 时空分析[J]. 南方水产科学, 2014, 10(2): 18-26. doi: 10.3969/j.issn.2095-0780.2014.02.003 [20] 化成君, 张衡, 伍玉梅, 等. 中东太平洋金枪鱼延绳钓中心渔场的时空变化[J]. 生态学杂志, 2014, 33(5): 1243-1247. [21] 郑波, 陈新军, 李纲. GLM 和 GAM 模型研究东黄海鲐资源渔场与环境因子的关系[J]. 水产学报, 2008, 32(3): 379-386. [22] AUSTIN M P, MEYERS J A. Current approaches to modelling the environmental niche of eucalypts: implication for management of forest biodiversity[J]. Forest Ecol Manag, 1996, 85(1/2/3): 95-106. doi: 10.1016/S0378-1127(96)03753-X
[23] 周甦芳, 樊伟, 崔雪森, 等. 环境因子对东海区帆式张网主要渔获物渔获量影响[J]. 应用生态学报, 2004(9): 1637-1640. doi: 10.3321/j.issn:1001-9332.2004.09.030 [24] BELLIDO J M, PIERCE G J, WANG J. Modelling intra-annual variation in abundance of squid Loligo forbesi in Scottish waters using generalised additive models[J]. Fish Res, 2001, 52(1/2): 23-39. doi: 10.1016/S0165-7836(01)00228-4
[25] SWARTZMAN G, HUANG C, KALUZNY S. Spatial analysis of Bering Sea groundfish survey data using generalized additive models[J]. Can J Fish Aquat Sci, 1992, 49(7): 1366-1378. doi: 10.1139/f92-152
[26] Burnham K,Anderson D. Model Selection and Multimodel Inference: a practical information theoretic approach [M]. 2nd ed. New York: Springer, 2002: 267-347.
[27] 戴澍蔚, 唐峰华, 樊伟. 北太平洋平洋公海日本鲭资源分布及其渔场环境特征[J]. 海洋渔业, 2017, 39(4): 372-382. doi: 10.3969/j.issn.1004-2490.2017.04.002 [28] 朱国平, 李纲, 郑晓琼, 等. 东海鲐鱼资源时空分布特征[J]. 生态科学, 2011, 30(l): 1-7. [29] 唐峰华, 樊伟, 伍玉梅, 等. 北太平洋柔鱼渔场资源与海洋环境关系的季节性变化[J]. 农业资源与环境学报, 2015(3): 242-249. [30] 苏杭, 陈新军, 汪金涛. 海表水温变动对东、黄海鲐鱼栖息地分布的影响[J]. 海洋学报, 2015, 37(6): 88-96. [31] 袁小楠, 陈新军, 李纲. 黄海日本鲐灯光围网渔场重心年际变化及其与环境因子关系[J]. 广东海洋大学学报, 2016, 36(3): 110-114. doi: 10.3969/j.issn.1673-9159.2016.03.018 [32] 王良明, 李渊, 张然, 等. 西北太平洋日本鲭资源丰度分布与表温和水温垂直结构的关系[J]. 中国海洋大学学报, 2019, 49(11): 29-38. -
期刊类型引用(5)
1. 郭美莎，李钰静，梁延鹏，马雲，王海琴，车佳祥，宋晓红. 对羟基苯甲酸丙酯(PrP)对雌性食蚊鱼的组织损伤及相关基因表达的影响. 生态毒理学报. 2024(03): 373-386 . 
百度学术
2. 蒋永生，周珊珊，周永东，徐开达，镇小蔓，焦礼诗，张秋宏，渠遥. 运输水温对曼氏无针乌贼成体死亡率及生理指标的影响. 南方水产科学. 2024(04): 107-115 . 本站查看
3. 何静怡，郑伟，黄卉，岑剑伟，赵永强，王田，魏涯，郝淑贤，杨少玲，陈琛. 不同温度、盐度条件对草鱼暂养及应激保活的影响. 大连海洋大学学报. 2024(04): 597-605 . 百度学术
4. 廖嘉仪，熊梓彤，李志力，刘其根，孙佳敏. 急性低温和复温对青田田鱼鳃组织代谢的影响. 南方农业学报. 2024(09): 2813-2823 . 百度学术
5. 李欣阳，艾力扎提·哈力夫，高居杰，姜志强. 急性低温胁迫对许氏平鲉鳃组织抗氧化能力和组织结构的影响. 河北渔业. 2024(11): 16-20 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 2329
- HTML全文浏览量: 961
- PDF下载量: 57
- 被引次数: 7
粤公网安备 44010502001741号
