Heterogeneity of mantle length-body mass relationship in different Ommastrephes bartramii populations based on linear mixed model
-
摘要:
为研究不同群体柔鱼 (Ommastrephes bartramii) 生长的异质性,文章根据2016—2018年在北太平洋采集的柔鱼样本,对该物种胴长和体质量的性别和月间差异进行了分析。构建了广义线性模型 (Generalized linear model, GLM) 和9个线性混合模型 (Linear mixed model, LMM),用于研究不同群体柔鱼的胴长和体质量关系 (m=aLb) 在性别及时间上的差异。结果显示,东、西部群体不同性别柔鱼胴长、体质量均存在极显著差异 (P<0.01);各月份间也均存在极显著差异 (P<0.01)。东部群体异速生长参数 (b) 的固定值 (2.763 4) 小于3,表示其为负异速增长,肥满度随胴长增加而减少,胴体趋于细长;西部群体b的固定值 (3.117 19) 大于3,表示其为正异速增长,肥满度随胴长增加而增大,胴体趋于宽厚。研究表明,性别和月份对不同群体柔鱼胴长和体质量关系均有显著影响,线性混合模型能将性别和月份的异质性通过随机效应在单个模型中更准确、方便地体现出来,进一步证实了此模型在个体生长异质性研究中的优势。
Abstract:We had collected neon flying squid (Ommastrephes bartramii) samples from 2016 to 2018 in North Pacific, and studied their gender and monthly differences of mantle length-body mass relationship, so as to understand the heterogeneity in the growth of different O. bartramii populations. One generalized linear (GLM) model and nine linear mixed models (LMM) were applied to study the difference in gender and time for mantle length-body mass relationship (m=aLb). The results show that there was a significant difference in the mantle length and body mass between female and male O. bartramii of eastern and western populations (P<0.01). Moreover, there was a significant difference between different months (P<0.01). The fixed value of b (2.763 4) was less than 3 in the eastern population, which indicates that O. bartramii was growing at negative allometry, and the degree of fullness was negatively related to mantle length, and the body tended to be slender. The fixed value of b (3.117 19) was more than 3 in the western population, which indicates that O. bartramii was growing at positive allometry, and the degree of fullness was positively related to the mantle length, and the body tended to be thick. The study indicates that month and gender have a significant effect on the mantle length-body mass relationship for different O. bartramii populations. The LMM reflects the heterogeneity of gender and month more accurately and conveniently in a single model through random effects, thereby further confirming the advantage of this model in the study of heterogeneity in individual growth.
-
河口盐沼湿地初级生产力高、生境异质性强的生态环境为众多水生生物提供了丰富的饵料资源和栖息生境[1],因此被认为是鱼类重要的育幼场所[2-4]。长江口湿地属于河口型冲击沙洲,在陆海物质交汇、咸淡水混合、径流和潮汐的相互作用下,形成了特殊而又多样的生态环境,主要有潮沟、盐沼植被和光滩3种亚生境类型[5],它们为水生生物提供了产卵场、育幼场、索饵场和洄游通道[6],且对维持流域和近海鱼类生物多样性、保障渔业资源有效补充方面发挥着重要作用。近年来长江口盐沼受到了众多人为因素的干扰并已造成严重影响,如水域污染引起鱼类死亡、渔获量下降以及饵料生物减少;滩涂围垦等各种水利工程建设使栖息生境和育幼场所破碎化;过度捕捞和毁灭性渔业生产方式对亲鱼和幼鱼资源造成了严重破坏[7]。中国花鲈 (Lateolabrax maculatus) 也是深受其害的鱼类之一。
中国花鲈广泛分布于我国沿海及河口,是长江口的重要经济鱼类[7],对河口进行选择性利用。夏、秋季活跃于河口近岸水域,11月至翌年1月到外海产卵[8]。国内外学者已对中国花鲈的分类学、人工繁育、养殖、营养学、食性等进行了大量相关研究[9-14],但有关其如何选择、利用育幼生境的研究尚不多见。
生态位宽度是指一个种群在一个群落中所利用的各种资源的总合[15]。时间生态位宽度反映物种在时间尺度上对资源利用的情况;空间生态位宽度反映物种在空间维度上的分布范围,其大小反映物种对不同空间资源利用的能力大小[16]。本研究根据在长江口九段沙湿地3种亚生境的调查资料,分析中国花鲈的生态位、栖息分布以及种群特征的变化情况,旨为丰富中国花鲈对长江口盐沼湿地的利用提供基础资料,也为中国花鲈资源保护和可持续利用提供科学依据,并为长江口盐沼湿地生态修复提供新视角。
1. 材料与方法
1.1 样品采集与测量
研究地点位于上海九段沙湿地国家级自然保护区 (121°50'E—122°05'E、31°03'N—31°38'N,图1),受亚热带季风气候的影响,四季分明,冬季寒冷干燥,夏季暖热多雨,年均降水量为1 145 mm,年平均气温为15.7 ℃,1月最低 (平均4.2 ℃),7月最高 (平均27.3 ℃)。九段沙附近潮汐属于非正规半日浅海潮,平均潮差为2.67 m,最大潮差为4.62 m。冬季为枯水季节 (11月—翌年1月),被高盐度咸水控制,夏季为丰水季节 (7—9月),被低盐度江水控制[17]。
本研究于2019年5月 (春季)、8月 (夏季)、11月 (秋季) 和2020年1月 (冬季) 利用插网在光滩、海三棱藨草 (Scirpus mariqueter) 植被区和潮沟中各设一站点进行调查。插网由网目为15 mm的网翼和网囊组成,其中网翼长5 m、高1 m,网口1 m×1 m,网囊长6 m。每次调查插网连续收集3 d渔获物,当场冰上冷藏并带回实验室,分类鉴定到种。对渔获物进行生物学测量,包括全长、体长 (精确至1 mm)、体质量 (精确至0.01 g) 等。根据孙帼英等[8]研究并结合解剖结果辨别中国花鲈的生长发育情况。
1.2 数据分析
采用Pinkas相对重要性指数 (Index of relative importance, IRI) 来分析优势种[18],公式为:
$$ {\rm{IRI}}=(N{\text{%}}+W{\text{%}})\times F{\text{%}} $$ (1) 式中
$ N{\text{%}} $ 为某种类的尾数在渔获物总尾数所占百分比;$ W{\text{%}} $ 为某种类的质量在渔获总质量所占百分比;$ F{\text{%}}$ 为某种类出现的站位数占调查站位总数百分比,将IRI>500定义为优势种[19]。应用Shannon-Winner公式计算生态位宽度[20]:
$$\begin{split} &\qquad \qquad \qquad \qquad \qquad \;\\ &{B}_{{i}}=-{\sum }_{j=1}^{r}\left({P}_{i_{\dot{J}}}\mathrm{ln}{P}_{i{j}}\right) \end{split}$$ (2) 式中
${P}_{i_{\dot{J}}}$ =$ {n}_{ij}/{N}_{ij} $ 和$ {P}_{kj}={n}_{kj}/{N}_{kj} $ 分别代表种$ i $ 和$ k $ 在第$ j $ 个资源状态下的个体数占该种所有个体数的比例,$ j $ 表示某个调查航次。生态位宽度反映物种对资源的利用能力,空间生态位宽度可反映某物种分布区域广泛与否,时间生态位宽度可反映物种时间上出现的长短[21]。采用摄食等级来描述摄食强度[22-23]。2. 结果
2.1 研究区域的环境特征
调查期间,九段沙3种亚生境的水温随季节发生明显变化,水温夏季最高、冬季最低,3种亚生境间水温无显著性差异 (P>0.05);盐度则相反,夏季最低,冬季最高,3种亚生境间盐度无显著性差异 (P>0.05);3种亚生境最大淹没水深在夏季均最高,其他月份相差不大,潮沟比海三棱藨草、光滩的更高且差异显著 (P<0.05,表1)。春季至夏季海三棱藨草生长旺盛,秋季至冬季海三棱藨草逐渐枯萎、消亡。
表 1 九段沙湿地不同亚生境水温、盐度和水深的季节变化Table 1. Seasonal variation of water temperature, salinity and water depth in different subhabitats of Jiuduansha Wetland季节
Season水温
Water temperature/℃盐度
Salinity最大水深
Maximum water depth/m潮沟
Tidal creek海三棱藨草
S. mariqueter光滩
Mudflat潮沟
Tidal creek海三棱藨草
S. mariqueter光滩
Mudflat潮沟
Tidal creek海三棱藨草
S. mariqueter光滩
Mudflat春季 Spring 13.5 16.1 13.6 1.83 1.46 1.49 3.5 1.6 1.7 夏季 Summer 28.2 28.3 28.5 1.04 0.43 0.81 4.3 2.1 2.2 秋季 Autumn 14.7 14.8 14.5 7.24 7.12 7.77 3.4 1.5 1.6 冬季 Winter 7.6 8.1 7.9 14.88 12.34 14.51 3.1 1.4 1.6 2.2 长江口九段沙盐沼湿地主要种及其生态位
调查期间共采集鱼类样本8目12科,29种,1 346尾。3种亚生境数量百分比N>1%的鱼类有5种:中国花鲈、斑尾刺虾虎鱼 (Acanthogobius ommaturus)、鮻 (Liza haematocheila)、小黄鱼 (Larimichthys polyactis)、黄姑鱼 (Nibea albiflora)。不同物种在3种亚生境的优势度有差异,中国花鲈、斑尾刺虾虎鱼在3种亚生境中均为优势种 (IRI>500),鮻、小黄鱼、黄姑鱼在部分亚生境中为优势种 (表2)。
表 2 九段沙3种亚生境主要种的相对重要性指数Table 2. Relative importance index of main species in three subhabitats of Jiuduansha Wetland种类
Species潮沟
Tidal creek海三棱藨草
S. mariqueter光滩
MudflatN/% W/% F/% IRI N/% W/% F/% IRI N/% W/% F/% IRI 中国花鲈 L. maculatus 26.5 24.8 100 5 130 35.9 40.5 100 7 645 17.7 43.5 100 6 119 斑尾刺虾虎鱼 A. ommaturus 23.8 63.2 100 8 698 8.3 7.9 75 1 216 2.8 9 100 1 180 鰉 L. haematocheila 1.5 4.9 75 480 20.5 40.7 75 6 117 6.1 28.3 100 3 444 小黄鱼 L. polyactis 40 1.6 25 1 040 12.8 0.2 25 326.6 17.1 0.8 50 897 黄姑鱼 N. jordan 1.6 0.1 25 43 1.9 0.02 25 48.6 22.3 5.7 25 700 中国花鲈的空间生态位宽度和时间生态位宽度变化较大,空间生态位宽度介于0.16~0.88,在春、夏、秋季依次递减,且春季和夏季时空间生态位宽度比其他优势种的大,春、夏季中国花鲈在3种亚生境的分布较秋、冬季更均匀;中国花鲈的时间生态位宽度介于0.39~0.78,在时间尺度上中国花鲈对潮沟的利用程度大于海三棱藨草植被区和光滩亚生境 (表3)。
表 3 九段沙3种亚生境主要鱼种生态位宽度Table 3. Niche breadth of main fish species in three subhabitats of Jiuduansha Wetland种类
Species空间生态位宽度
Spatial niche breadth时间生态位宽度
Temporal niche breadth春季
Spring夏季
Summer秋季
Autumn冬季
Winter潮沟
Tidal creek海三棱藨草
S. mariqueter光滩
Mudflat中国花鲈 L. maculatus 0.88 0.52 0.16 0.26 0.78 0.66 0.39 斑尾刺虾虎鱼 A. ommaturus 0.17 0.41 0.45 0.51 0.97 0.22 0.18 鰉 L. haematocheila 0.57 0.21 0.45 0.89 0.27 0.90 0.72 小黄鱼 L. polyactis 0.57 0 0.01 0 0.16 0.14 0.28 黄姑鱼 N. albiflora 0 0.57 0 0 0.29 0.11 0.17 2.3 长江口中国花鲈分布特征
本研究共采集到中国花鲈449尾,3种亚生境数量差异明显,其中潮沟266尾,海三棱藨草植被区121尾,光滩62尾。各亚生境中国花鲈的数量在不同季节间变化显著,春、夏、秋和冬季各亚生境中国花鲈的数量先逐渐减少、后增加,春季最多,秋季最少 (图2)。
中国花鲈在3种亚生境的数量与总数量的比例随季节发生明显变化,潮沟中的数量占比在春、夏、秋、冬季逐渐减小;海三棱藨草和光滩的数量占比在春、夏、秋、冬季逐渐增加 (图3)。
2.4 中国花鲈的体长特征
中国花鲈性成熟时体长约350 mm,本研究捕获的中国花鲈体长最大为339 mm,最小为20 mm,经对比和解剖发现所捕获的中国花鲈均为幼鱼。体长小于70 mm的个体最多,占比94.24%,且该体长规格的个体数在不同亚生境中的比例存在差异,占比大小依次为潮沟>海三棱藨草>光滩。
中国花鲈的平均体长在3种亚生境均表现出春、夏、秋季逐渐增大的趋势,海三棱藨草亚生境中的平均体长在春季与夏季之间的变化幅度明显大于光滩和潮沟 (图4)。
各季节各亚生境中国花鲈平均体长的标准偏差较大,个体间的体长差异大。为进一步分析各亚生境间的体长差异,本研究将中国花鲈的体长数据进行合理分组。结果显示,春季中国花鲈体长主要集中在10~50 mm体长组,对该体长区间的3种亚生境体长数据进行单因子方差分析发现,潮沟、海三棱藨草和光滩亚生境中国花鲈的平均体长分别为34.13、29.27和32.15 mm,海三棱藨草亚生境中国花鲈的体长与其他亚生境有显著性差异 (P<0.05),潮沟亚生境与光滩亚生境的中国花鲈体长无显著性差异 (P>0.05);夏季各亚生境体长为91~130 mm的个体占比均最大,此体长组潮沟、海三棱藨草和光滩亚生境的平均体长分别为101.85、113.17和112.4 mm,潮沟亚生境中的中国花鲈体长与海三棱藨草亚生境和光滩亚生境差异显著 (P<0.05);秋季各亚生境体长分布较为分散,未发现体长小于90 mm的个体,各亚生境中的中国花鲈体长无显著性差异 (P>0.05);冬季中国花鲈体长集中在131~210 mm,未发现体长大于251 mm的个体 (图5)。
春、夏、秋季除了占比较大的小规格个体外,还存在大于251 mm的大规格个体,且此规格的个体在海三棱藨草和光滩亚生境出现的比例更高,但冬季3种亚生境中均无此规格的个体出现。
2.5 各亚生境中国花鲈平均摄食强度变化
中国花鲈的平均摄食等级在潮沟亚生境呈现从春季到冬季逐渐增大的趋势;在光滩亚生境呈现从春季到秋季逐渐增大而冬季骤减的趋势;在海三棱藨草亚生境呈现夏季最大、秋冬季逐渐减小的趋势。
中国花鲈平均摄食等级在春、夏、秋季表现为海三棱藨草>潮沟>光滩;冬季较其他季节特殊,平均摄食等级表现为潮沟>海三棱藨草>光滩 (图6)。
3. 讨论
3.1 优势种与生态位
本研究发现长江口盐沼湿地优势种为中国花鲈、斑尾刺虾虎鱼、鮻和小黄鱼,且中国花鲈优势度比其他优势种大,说明中国花鲈在九段沙湿地的群落中起主导和控制作用。本研究与杨刚等[24]在长江口的调查结果——优势种为中国花鲈、鮻和窄体舌鳎 (Cynoglossus gracilis)及庄平等[25]的调查结果——优势种为长蛇 (Saurogobio dumerili)、中国花鲈、窄体舌鳎、刀鲚 (Coilia nasus)、鮻和鲻 (Mugil cephalus) 略有差异。导致这种差异主要有2个原因:1) 时间。本研究进行4次采样,对应春、夏、秋、冬4个季节,较全面地反映了中国花鲈在九段沙湿地的动态变化,杨刚等[24]采样时间是2010年5月、8月、11月,庄平等[25]采样时间是2004—2006年每年的5—8月,与本研究采样时间略有差异。2) 生境。本研究采样地点包括3种不同的亚生境,能更客观地反映九段沙盐沼湿地鱼类群落特征,而杨刚等[24]、庄平等[25]对采样地点未作详细划分,可能因某些物种的选择偏好性而低估了盐沼湿地鱼类群落的物种丰度[26]。
本研究中中国花鲈的空间生态位宽度在春、夏、秋季逐渐减小,秋季的较小是因为光滩亚生境中国花鲈的数量偏少,为潮沟亚生境和海三棱藨草亚生境中数量的1/2,说明中国花鲈随季节变化在3种亚生境中的分布差异增大,这可能是因为季节变化引起水温、水位、盐度等环境因素的变化,进而影响各亚生境中的环境和饵料生物。春、夏季中国花鲈的空间生态位较秋、冬季大,说明春、夏季中国花鲈在3种亚生境的分布较秋、冬季更均匀,此时中国花鲈对生境的选择利用差异较小,侧面反映了春、夏季更适合中国花鲈在河口盐沼湿地育幼和生长。
时间生态位宽度反映物种在时间尺度上对资源利用的情况[16,27],中国花鲈时间生态位宽度在潮沟、海三棱藨草、光滩亚生境间逐渐减小。光滩亚生境中国花鲈的时间生态位最小,说明光滩亚生境中的中国花鲈在时间上较其他亚生境的分布差异大。潮沟亚生境中国花鲈的时间生态位宽度较其他亚生境大,说明潮沟亚生境中的中国花鲈在时间尺度上较其他亚生境的分布差异小,这可能是因为潮沟亚生境环境复杂,相较其他亚生境,季节变化对潮沟的干扰较小。
3.2 九段沙湿地中国花鲈的时空分布差异
河口洄游性鱼类一般受到自身生活史及生物学特性的影响,其分布具有明显的时空交替变化特征。本研究结果表明中国花鲈作为近海咸淡水鱼类,在不同亚生境和不同季节均存在数量差异 (潮沟的数量>海三棱藨草>光滩),春季数量最多。中国花鲈为肉食性鱼类,具有主动捕食的特征,其季节分布差异受自身生物学特性和饵料生物的影响[28-29]。作为洄游性鱼类,中国花鲈在秋季至冬季性成熟,洄游到水位深、盐度高的水域产卵,随后仔、稚鱼在早春游到近岸低盐度水域索饵[7],因此,本研究调查到的数量表现出春季最多、随后逐渐减小、冬季有所回升的趋势。春季,大量中国花鲈幼体集中游向河口盐沼湿地索饵、躲避敌害,此区域的低盐度和丰富的适口饵料生物适合其在此栖息[30-31];夏季,中国花鲈逐渐长大,游泳能力增强,食性也发生转变,为满足生长需求会游向其他水域以获取更多的食物[14,32-33];秋季,成体和亚成体需游向深水位、高盐度水域为繁育做准备,因此秋季在盐沼湿地捕到的中国花鲈个体很少;冬季,产过卵的中国花鲈会游到河口盐沼湿地索饵,这可能是冬季捕获个体数量多于秋季的原因。有研究表明长江口潮间带是中国花鲈的育肥场所[8],中国花鲈在长江口3种亚生境的空间分布差异,可能与不同亚生境的生境特质有关。潮沟和海三棱藨草植被区比光滩亚生境有更丰富的底栖生物,3种亚生境浮游生物优势种和丰度也存在差异[34],且本研究发现3种亚生境中中国花鲈的摄食等级有明显差异,因此,中国花鲈数量在长江口盐沼湿地的时空变化可能是基于其自身生物学特性以及生境间的差异、季节性利用河口育幼功能的结果。
3.3 中国花鲈的体长特征
国内外的一些研究表明,盐沼湿地是一些近岸渔业物种的重要栖息生境[4,35]。张衡和朱国平[36]通过在长江口的调查推测长江口盐沼湿地是重要的育幼场,金斌松[37]在崇明东滩和九段沙湿地调查发现潮沟中存在大量的中国花鲈幼鱼,平均体长为24 mm。本研究春季中国花鲈的平均体长为43.6 mm,经解剖并与孙帼英等[8]研究对比发现,本研究所捕中国花鲈均为幼鱼;3种亚生境中国花鲈体长小于50 mm数量占比均达86%以上,幼体比例较高,说明长江口盐沼湿地是中国花鲈的育幼场。中国花鲈在3种亚生境的体长存在差异,尤其在春季和夏季。相较于潮沟,大规格个体在海三棱藨草、光滩出现的比例更高,这可能因为海三棱藨草、光滩亚生境较为开阔,更利于中国花鲈追捕游泳动物。
-
图 4 柔鱼不同月份胴长体质量分布图
a. 东部群体胴长分布图;b. 东部群体体质量分布图;c. 西部群体胴长分布图;d. 西部群体体质量分布图
Figure 4. Distribution of mantle length and body mass of O. bartramii in different months
a. Distribution of mantle length for eastern population; b. Distribution of body mass for eastern population; c. Distribution of mantle length for western population t; d. Distribution of body mass for western population
图 5 不同群体柔鱼体质量-胴长关系
a. 东部群体雌性体质量-胴长关系;b. 东部群体雄性体质量-胴长关系;c. 西部群体雌性体质量-胴长关系;d. 西部群体雄性性体质量-胴长关系
Figure 5. Mantle length-body mass relationship for different O. bartramii populations
a. Relationship between BM and ML of female for eastern population; b. Relationship between BM and ML of male for eastern population; c. Relationship between BM and ML of female for western population; d. Relationship between BM and ML of male for western population
表 1 不同月份柔鱼样本基本信息
Table 1 Basic information of O. bartramii samples in different months
群体
Population年份
Year月份
Month样本数量
Number of sample胴长范围
Range of mantle length/mm东部群体 Eastern population 2016—2017 5 394 163~500 6 456 212~516 西部群体 Western population 2016—2018 7 124 187~290 8 283 195~357 9 331 199~386 10 165 190~377 11 105 215~421 表 2 不同群体柔鱼胴长-体质量关系模型的拟合效果指标
Table 2 Fitting effect index (AIC) of mantle length-body mass relationship models for different O. bartramii populations
种群
Population模型缩写
Model abbreviation模型
Model对数转化
Log-transformed赤池信息量准则
AIC均方根误差
RMSE东部群体
Eastern populationLMR m = a×Lb ln(m)=ln(a)+b×ln(L) −466.402 0.183 265 RSS m=a×L(b+RSS) ln(m)=ln(a)+(b+RSS)×ln(L) −472.338 0.181 617 RSI m=a×exp(RSI)×Lb ln(m)=[ln(a)+RSI]+b×ln(L) −472.27 0.181 641 RSI&S m=a×exp(RSI)×L(b+RSS) ln(m)=[ln(a)+RSI]+(b+RSS)×ln(L) −469.057 0.181 428 RMS m=a×L(b+RMS) ln(m)=ln(a)+(b+RMS)×ln(L) −489.107 0.179 784 RMI m=a×exp(RMI)×Lb ln(m)=[ln(a)+RMI]+b×ln(L) −492.772 0.179367 RMI&S m=a×exp(RMI)×L(b+RSI) ln(m)=[ln(a)+RMI]+(b+RMS)×ln(L) −595.137 0.166 468 RS&MS m=a×L(b+RSS+RMS) ln(m)=ln(a)+(b+RSS + RMS)×ln(L) −494.452 0.178 558 RS&MI m=a×exp(RSI)×exp(RMI)×Lb ln(m)=[ln(a)+RSI + RMI]+b×ln(L) −497.694 0.178 211 RS&MI&S m=a×exp(RSI)×exp(RMS)×L(b+RSI+RSI) ln(m)=[ln(a) + RSI + RMS]+(b+RSI + RMI)×ln(L) −596.138 0.165 394 西部群体
Western populationLWR m = a×Lb ln(m)=ln(a)+b×ln(L) −1 795.63 0.099 003 RSS m=a×L(b+RSS) ln(m)=ln(a)+(b+RSS)×ln(L) −1 797.15 0.098 614 RSI m=a×exp(RSI)×Lb ln(m)=[ln(a)+RSI]+b×ln(L) −1 797.31 0.098 603 RSI&S m=a×exp(RSI)×L(b+RSS) ln(m)=[ln(a)+RSI]+(b+RSS)×ln(L) −1 796.72 0.098 512 RMS m=a×L(b+RMS) ln(m)=ln(a)+(b+RMS)×ln(L) −1 807.04 0.097 626 RMI m=a×exp(RMI)×Lb ln(m)=[ln(a)+RMI]+b×ln(L) −1 807.41 0.097 615 RMI&S m=a×exp(RMI)×L(b+RSI) ln(m)=[ln(a)+RMI]+(b+RMS)×ln(L) −1 807.59 0.097 215 RS&MS m=a×L(b+RSS+RMS) ln(m)=ln(a)+(b+RSS + RMS)×ln(L) −1 808.14 0.097 347 RS&MI m=a×exp(RSI)×exp(RMI)×Lb ln(m)=[ln(a)+RSI + RMI]+b×ln(L) −1 808.56 0.097 333 RS&MI&S m=a×exp(RSI)×exp(RMS)×L(b+RSI+RSI) ln(m)=[ln(a) + RSI + RMS]+(b+RSI + RMI)×ln(L) −1 808.09 0.096 976 注:第一列为柔鱼的不同群体,第二列为第三、四列模型的缩写 RSS、RSI、RSI&S、RMS、RMI、RMI&S、RS&MS、RS&MI 和 RS&MI&S 分别为斜率 b 的性别随机效应、截距 ln(a) 的性别随机效应、截距 ln(a) 与斜率 b的性别随机效应、斜率 b 的月份随机效应、截距 ln(a) 的月份随机效应、截距 ln(a) 与斜率 b 的月份随机效应、斜率 b的性别与月份随机效应、截距 ln(a) 的性别与月份随机效应,以及截距 ln(a) 与斜率 b 的性别与月份随机效应Note: The first column shows different O. bartramii populations. The second column shows the abbreviations of models detailed in the third and fourth columns. RSS, RSI, RSI&S, RMS, RMI, RMI&S, RS&MS, RS&MI和RS&MI&S indicate random effects on slope (b) in sexual difference, random effects on intercept [ln (a)] in sexual difference, random effects on intercept [ln (a)] and slope (b) in sexual difference, random effects on slope (b) in months, random effects on intercept [ln (a)] in months, random effects on intercept [ln (a)] and slope (b) in months, random effects on slope (b) in sexual difference and months, random effects on intercept [ln (a)] fin sexual difference and months, and random effects on intercept [ln (a)] and slope (b) in sexual difference and months, respectively. 表 3 不同月份东西部群体柔鱼胴长与体质量关系
Table 3 Mantle length-body mass relationship for different O. bartramii populations in different months
群体
Population月份
Month胴长与体质量关系
Mantle length-mass relationship样本数量
Number of samples判定系数R2
Coefficient of determinationP 东部群体
Eastern population5 m =1.9×10−3L2.295 7 394 0.670 6 0.001 6 m =1×10−5L3.160 1 456 0.973 8 0.001 西部群体
Western population7 m = 9×10−6L3.205 3 124 0.889 5 0.001 8 m = 2×10−5L3.090 9 283 0.934 5 0.001 9 m = 2×10−5L3.059 7 331 0.952 0 0.001 10 m = 1×10−5L3.178 3 165 0.972 8 0.001 11 m = 5×10−5L2.902 8 105 0.947 7 0.001 表 4 不同性别东西部群体柔鱼胴长与体质量关系
Table 4 Mantle length-body mass relationship of O. bartramii of different sex for different populations
群体
Population性别
Sex胴长与体质量关系
Mantle length-mass relationship样本数量
Number of samples判定系数R2
Coefficient of determinationP 东部群体
Eastern population雌性 m=3×10−5L2.981 8 744 0.892 9 0.001 雄性 m=2×10−4L2.621 0 106 0.794 3 0.001 西部群体
Western population雌性 m=1×10−5L3.159 9 512 0.968 1 0.001 雄性 m =2×10−5L3.103 5 496 0.941 5 0.001 -
[1] 王尧耕, 陈新军. 世界大洋性经济柔鱼类资源及其渔业[M]. 北京: 海洋出版社, 2005: 124-137. [2] MURATA M. Oceanic resources of squids[J]. Mar Freshw Behav Phy, 1990, 18(1): 19-71. doi: 10.1080/10236249009378779
[3] BOWER J R, ICHII T. The red flying squid (Ommastrephes bartramii): a review of recent research and the fishery in Japan[J]. Fish Res, 2005, 76(1): 39-55. doi: 10.1016/j.fishres.2005.05.009
[4] 陈新军. 渔业资源与渔场学[M]. 北京: 海洋出版社, 2004: 29-212. [5] YATSU A, TANAKA H, MORI J. Population structure of the neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean[J]. Jpn Mar Fish Resour Res Cent, 1998: 31-48.
[6] FANG Z, LIU B, LI J, et al. Stock identification of neon flying squid (Ommastrephes bartramii) in the North Pacific Ocean on the basis of beak and statolith morphology[J]. Sci Mar, 2014, 78(2): 239-248. doi: 10.3989/scimar.03991.06A
[7] 方舟, 陈新军, 瞿俊跃, 等. 北太平洋柔鱼角质颚形态及生长年间差异[J]. 上海海洋大学学报, 2020, 29(1): 109-120. doi: 10.12024/jsou.20181202485 [8] WATANABE H, KUBODERA T, ICHII T, et al. Feeding habits of neon flying squid Ommastrephes bartramii in the transitional region of the central North Pacific[J]. Mar Ecol Prog Ser, 2004, 266(1): 173-184.
[9] CHEN C S, CHIU T S. Variations of life history parameters in two geographical groups of the neon flying squid, Ommastrephes bartramii, from the North Pacific[J]. Fish Res, 2003, 63(3): 349-366. doi: 10.1016/S0165-7836(03)00101-2
[10] YATSU A. Distribution of flying squid, Ommastrephes bartramii, in the north Pacific based on a jigging survey, 1976−1983[J]. Bull Natl Res Inst Far Seas Fish, 1992(29): 13-37.
[11] YATSU A, MORI J. Early growth of the autumn cohort of neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean[J]. Fish Res, 2000, 45(2): 189-194. doi: 10.1016/S0165-7836(99)00112-5
[12] SUN M C, ZHANG S H, QIAN W G. Approach on biological characteristics of Ommastrephes bartramii on the central and eastern squid fishing grounds in the North Pacific[J]. Mar Tieher, 2001, 23(1): 21-24.
[13] YATSU A, MIDORIKAWA S, SHIMADA T, et al. Age and growth of the neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean[J]. Fish Res, 1997, 29(3): 257-270. doi: 10.1016/S0165-7836(96)00541-3
[14] MCINTYRE A D. Cephalopod life cycles. Vol. 1. Species accounts[J]. J Anim Ecol, 1986, 55(1): 395-396.
[15] GRANADOS-AMORES J, SALINAS-ZAVALA C A, FLORES-ORTEGA J R, et al. Length-weight relationship and condition factor for 7 loliginid squid species in Mexican waters[J]. Cienc Mar, 2019, 45(4): 175-180. doi: 10.7773/cm.v45i4.3015
[16] VASCONCELOS P, PEREIRA F, CARVALHO A N, et al. Weight-length relationships and relative growth of the cuttlefish (Sepia officinalis): causes and effects of hypoallometry[J]. Thalassas: Int J Mar Sci, 2018, 32(4): 323-331.
[17] 衷思剑, 麻秋云, 刘淑德, 等. 基于线性混合效应模型的黄鮟体长体重关系的时空差异[J]. 中国水产科学, 2018, 25(6): 152-160. [18] XU H, SUN Y J, WANG X J, et al. Linear mixed-effects models to describe individual tree crown width for China-fir in Fujian province, southeast China[J]. PLoS One, 2015, 10(4): e0122257. doi: 10.1371/journal.pone.0122257
[19] MA Q Y, JIAO Y, REN Y P. Linear mixed-effects models to describe length-weight relationships for yellow croaker (Larimichthy spolyactis) along the north coast of China[J]. PLoS One, 2017, 12(2): e0171811. doi: 10.1371/journal.pone.0171811
[20] 瞿俊跃, 方舟, 陈新军. 基于线性混合模型的蓝点马鲛叉长体重生长关系的月间及性别差异[J]. 中国水产科学, 2020, 8: 953-960. [21] 廖为耕. 浅谈北太平洋之海洋环境及赤鱿之洄游[J]. 中国水产(台刊), 1980, 457: 39-47. [22] LIPINSKI M R, UNDERHILL L G. Sexual maturation in squid: quantum or continuum[J]. S Afr J Mar Sci, 1995, 15: 207-223. doi: 10.2989/02577619509504844
[23] 李春喜, 邵云, 姜丽娜. 生物统计学[M]. 北京: 科学出版社, 2008: 85-110. [24] GUEORGUIEVA R, KRYSTAL J H. Move over anova: progress in analyzing repeated-measures data and its reflection in papers published in the archives of general psychiatry[J]. Arch Gen Psychiat, 2004, 61(3): 310-317. doi: 10.1001/archpsyc.61.3.310
[25] 王斌会. 多元统计分析及R语言建模[M]. 广州: 暨南大学出版社, 2010: 157-158. [26] KEYS A B. The weight-length relationship in fishes[J]. P Natl A Sci USA, 1928, 14(12): 922-925. doi: 10.1073/pnas.14.12.922
[27] BURNHAM K P, ANDERSON D R. Model selection and multi-model inference: a practical information theoretic approach[M]. 2nd ed. New York: Springer Press, 2002: 484-488.
[28] BERNAL-RUSIEL J L, GREVE D N, REUTER M, et al. Statistical analysis of longitudinal neuroimage data with Linear Mixed Effects models[J]. NeuroImage, 2013, 66(2): 249-260.
[29] 黄洪亮, 郑元甲, 程家骅. 北太平洋海区柔鱼生物学特征研究[J]. 海洋渔业, 2003, 25(3): 21-24. [30] ICHII T, MAHAPATRA K, OKAMURA H, et al. Stock assessment of the autumn cohort of neon flying squid (Ommastrephes bartramii) in the North Pacific based on past large-scale high seas driftnet fishery data[J]. Fish Res, 2006, 78(2): 286-297.
[31] KATUGIN O N. Patterns of genetic variability and population structure in the North Pacific squids Ommastrephes bartramii, Todarodes pacificus, and Berryteuthis magister[J]. Bull Mar Sci, 2002, 71(1): 383-420.
[32] HILE R. Age and growth of the cisco, Leucichthys artedi (le Sueur), in the lakes of the Northeastern highlands, Wisconsin[J]. Bull Bur Fish, 1936, 48(11): 211-317.
[33] KEITH A, BIGELOW, KEVIN C, et al. Hatch dates and growth of Ommastrephes bartramii paralarvae from Hawaiian waters as determined from statolith analysis[J]. Rec Adv Fish Biol, 1993: 15-24.
[34] 唐峰华, 伍玉梅, 樊伟. 2010年北太平洋海区巴特柔鱼(Ommastrephe bartrami)的生物学特征[J]. 海洋科学进展, 2012, 30(2): 257-265. doi: 10.3969/j.issn.1671-6647.2012.02.013 [35] FORSYTHE J W. Accounting for the effect of temperature on squid growth in nature: from hypothesis to practice[J]. Mar Freshw Res, 2004, 55(4): 331-339. doi: 10.1071/MF03146
[36] ICHII T, MAHAPATRA K, SAKAI M, et al. Life history of the neon flying squid: effect of the oceanographic regime in the North Pacific Ocean[J]. Mar Ecol Prog Ser, 2009, 378: 1-11. doi: 10.3354/meps07873
[37] 章寒, 郑基, 虞聪达, 等. 东太平洋赤道海域茎柔鱼主要生物学特性比较研究[J]. 中国水产科学, 2019, 26(4): 745-755. [38] DOMINGUES P M, SYKES A, ANDRADE J P. The effects of temperature in the life cycle of two consecutive generations of the cuttlefish Sepia officinalis (Linnaeus, 1758), cultured in the Algarve (South Portugal)[J]. Aquacult Int, 2002, 10(3): 207-220. doi: 10.1023/A:1022148802078
[39] DOMINGUES P M, POIRIER R, DICKEL L, et al. Effects of culture density and live prey on growth and survival of juvenile cuttlefish, Sepia officinalis[J]. Aquacult Int, 2003, 11(3): 225-242. doi: 10.1023/A:1024803802486
[40] 陈锦辉, 王学昉, 田思泉, 等. 基于混合效应模型的印度洋大眼金枪鱼生长特征异质性分析[J]. 中国水产科学, 2020, 27(5): 570-578.