富齿岩虫南移养殖与营养成分分析

李帅鹏, 徐斌, 王坤, 黄民巧, 杨大佐, 周一兵

李帅鹏, 徐斌, 王坤, 黄民巧, 杨大佐, 周一兵. 富齿岩虫南移养殖与营养成分分析[J]. 南方水产科学, 2021, 17(6): 93-100. DOI: 10.12131/20200222
引用本文: 李帅鹏, 徐斌, 王坤, 黄民巧, 杨大佐, 周一兵. 富齿岩虫南移养殖与营养成分分析[J]. 南方水产科学, 2021, 17(6): 93-100. DOI: 10.12131/20200222
LI Shuaipeng, XU Bin, WANG Kun, HUANG Minqiao, YANG Dazuo, ZHOU Yibing. Southwards cultivation and nutrition composition analysis of Marphysa maxidenticulata[J]. South China Fisheries Science, 2021, 17(6): 93-100. DOI: 10.12131/20200222
Citation: LI Shuaipeng, XU Bin, WANG Kun, HUANG Minqiao, YANG Dazuo, ZHOU Yibing. Southwards cultivation and nutrition composition analysis of Marphysa maxidenticulata[J]. South China Fisheries Science, 2021, 17(6): 93-100. DOI: 10.12131/20200222

富齿岩虫南移养殖与营养成分分析

基金项目: 国家重点研发计划项目 (2019YFC1407704);挪威研究理事会项目 (RCN: 280836);大连市杰出青年科学基金项目(2016RJ06);“兴辽计划”领军人才项目 (XLYC1805009)
详细信息
    作者简介:

    李帅鹏 (1992—),男,硕士研究生,研究方向为水产养殖。E-mail: 1345433665@qq.com

    通讯作者:

    杨大佐 (1979—),男,博士,研究员,从事多毛类动物养殖与资源恢复研究。E-mail: dzyang@dlou.edu.cn

  • 中图分类号: S 968.9

Southwards cultivation and nutrition composition analysis of Marphysa maxidenticulata

  • 摘要: 该研究通过南移养殖的方法,测定了富齿岩虫 (Marphysa maxidenticulata) 在广西北海地区的生长状况,并比较了南北方虫体营养成分。结果显示,经240 d的养殖,虫体平均体长达 (12.66±1.52) cm,平均体质量达 (2.30±0.86) g。试验前期,海泥底质组虫体的特定生长率显著高于粗砂底质组;而粗砂底质组的虫体在养殖后期表现出较高的特定生长率。营养组成分析结果显示,广西北海养殖和山东乳山野生虫体均含有18种氨基酸,但其基本营养组成存在显著差异 (P<0.05)。北海养殖虫体氨基酸总量 (TAA)、呈味氨基酸 (DAA)、非必需氨基酸 (NEAA)、半必需氨基酸 (SEAA)、必需氨基酸 (EAA) 含量均高于乳山野生虫体。在脂肪酸组成方面,2个地区的富齿岩虫均含有27种脂肪酸,且以棕榈酸和油酸为主要脂肪酸;乳山野生虫体饱和脂肪酸含量和多不饱和脂肪酸含量与广西北海养殖虫体无显著差异;北海养殖虫体单不饱和脂肪酸总含量显著高于乳山野生虫体 (P<0.05)。
    Abstract: In this paper, we investigated the growth condition and nutritional composition of Marphysa maxidenticulata cultured in Beihai by using the "northern worm and southern cultivation" method. The results show that the average body length of M. maxidenticulata grew up to (12.66±1.52) cm and the average body mass was (2.30±0.86) g within 240 d. The specific growth rate of M. maxidenticulata in mud sediments was significantly higher than that of the coarse sands at the early phase of culture, but the late phase of culture was just opposite. The results of nutrition composition analysis show that the nutritional composition contained 18 species amino acids in both cultured M. maxidenticulata in Beihai and wild M. maxidenticulata in Rushan, but their contents were significantly different (P<0.05). The contents of total amino acid, delicious amino acid, non-essential amino acid, semi-essential amino acid and essential amino-acid of cultured M. maxidenticulata in Beihai were higher than those of wild M. maxidenticulata in Rushan. In terms of fatty acid composition, 27 fatty acids had been determined in both cultured and wild M. maxidenticulata, and palmitic acid and oleic acid were the main fatty acids. There was no significant difference in the contents of saturated fatty acids and polyunsaturated fatty acids between cultured M. maxidenticulata in Beihai and wild M. maxidenticulata in Rushan. But the content of total monounsaturated fatty acids of cultured M. maxidenticulata in Beihai was significantly higher than that of wild M. maxidenticulata in Rushan (P<0.05).
  • 深水金线鱼 (Nemipterus bathybius) 属鲈形目、金线鱼科、金线鱼属,为暖水性底层经济鱼类,在我国主要分布于南海北部和台湾周边海域,是底拖网渔业的主要渔获物之一[1-2]。体型较大的深水金线鱼个体通常以甲壳类动物、小型鱼类和头足类为食,较小的个体以桡足类、介形虫和片脚类为食。深水金线鱼生长速度快且性成熟较早,具有分批产卵、生命周期短等生物学特征,可以承受一定的捕捞压力[3]。随着南海北部近海渔业资源捕捞压力的持续上升,深水金线鱼资源状况也开始呈现过度开发现象,近年来其渔获物多为幼鱼,种群结构以低龄化个体为主[4-6]。南海幅员辽阔,但我国在南海大部分的捕捞活动和产量都集中在北部陆架区。北部陆架区是重要的经济鱼类产卵场,也是重要的渔场[7]。改革开放以来,随着沿海地区经济的快速发展,人民群众对海产品的需求量持续上升,海洋渔船数量和捕捞努力量不断增加[8],导致了南海北部近海渔业资源的衰退,特别是典型的半封闭海湾[9-10]。Zhang等[11]对北部湾9种主要经济鱼类进行资源评估,结果表明9种鱼类在不同时期均遭受了过度捕捞。《中国渔业统计年鉴》的全国海洋捕捞鱼类产量数据显示,2018—2022年南海金线鱼类产量平均为31.4万t,约占南海海洋鱼类总捕捞产量的16.2%,该类群经济价值高,其中深水金线鱼是主要的优势种类之一。以往对南海北部深水金线鱼种群的资源评估多针对生长、死亡和性成熟等生物学特征,且评估数据多为2010年以前的。对北部湾深水金线鱼种群参数的评估研究[4,11]表明,其生物学参数年际变化明显,种群个体虽生长速度有所增加,但小型化现象突出,捕捞个体多为幼鱼。

    渔业资源评估是制定合理管理策略的基础,科学有效的评估结果有助于实现渔业资源的可持续开发。模型是进行渔业资源评估的重要工具,数据缺乏条件下的资源评估模型主要分为单鱼种评估方法、多鱼种评估方法和生态系统模型,目前国际上常用的是单鱼种的数据缺乏方法[12]。相对于其他渔业资源数据,长度数据的收集操作简便、结果读取直接、采集成本低,因此基于长度数据的渔业资源评估方法近些年发展较为迅速。本研究根据2014—2019年南海北部底拖网调查的渔业生物学资料,使用基于长度的贝叶斯生物量评估 (Length-based Bayesian biomass estimation method, LBB) 和基于长度的繁殖潜力比 (Length-based spawning potential ratio, LBSPR) 2种模型,对南海北部深水金线鱼种群资源现状进行评估,以期为该种群的科学管理和可持续利用提供技术支撑。

    本研究所用样品数据均来自2014—2019年中国水产科学研究院南海水产研究所在南海北部海域底拖网调查资料,调查站点根据渔业资源按水深分布[13]的特点设计(图1)。所有航次的调查采样与生物学测定均按照GB/T 12763.3—2007《海洋调查规范》进行。采集的深水金线鱼样本冷冻保存,带回实验室进行生物学测量和鉴定,包括体长、体质量、性成熟度、胃饱满度等。每个站位深水金线鱼样本多于50尾时随机选取50尾进行测定,样本不足50尾时则全部测定,体长和体质量的测定分别精确到1 mm和1 g。12个航次共测定深水金线鱼3 059尾,样品信息见表1

    图  1  南海北部近海研究区域及站点分布
    Figure  1.  Study area and station distribution in offshore northern South China Sea
    表  1  南海北部深水金线鱼的样品信息
    Table  1.  Sample information of N. bathybius in northern South China Sea
    年份
    Year
    月份
    Month
    样本数
    Sampling
    number
    体长 Body length/mm
    范围
    Range
    均值±标准差
    Mean±Standard
    deviation
    2014 7 293 93~193 121.2±16.3
    10 225 77~184 119.9±21.6
    2015 1 120 73~172 121.8±18.7
    4 330 75~189 111.4±19.2
    2016 7 304 84~200 121.8±19.4
    10 78 74~178 122.1±23.8
    2017 1 71 72~176 109.3±18.7
    4 197 72~158 115.0±15.3
    2018 3 220 75~176 113.8±19.0
    8 87 78~182 126.0±24.2
    2019 3 524 73~176 106.6±16.8
    9 610 65~192 117.3±19.0
    下载: 导出CSV 
    | 显示表格

    采用幂指数拟合深水金线鱼体长与体质量关系,计算公式为:

    $$ W {\mathrm{=}}a{L}^{b} $$ (1)

    式中:W为体质量 (g);L为体长 (mm);a为生长的肥满度指数;b为幂指数。

    深水金线鱼的生长关系可通过von Bertalanffy生长方程[4]进行拟合:

    $$L_{{t}}{\rm{=}}L_{\rm {inf }}\left\{1{\text{−}}{\mathrm{exp}} \left[{\text{−}}K\left(t{\text{−}}t_0\right)\right]\right\} $$ (2)

    式中:Ltt龄的体长;Linf为渐近体长;K为生长系数;t0为体长为0时的理论年龄。本研究使用R软件包TropFishR[14]中的ELEFAN (Electronic length frequency analysis) 技术估算深水金线鱼的渐近体长 (Linf) 和生长系数 (K)。参数t0采用以下Pauly经验方程[15]进行计算:

    $$ \lg\left({\text{−}}t_0\right){\mathrm{=}}{\text{−}}0.392\,2{\text{−}}0.275 \lg L_{\rm {inf }}{\text{−}}1.038 \lg K $$ (3)

    根据以下经验公式[16]计算生长性能指数 (Φ):

    $$ \Phi{\mathrm{=}}\lg K{\mathrm{+}}2 \lg L_{\rm {inf }} $$ (4)

    自然死亡率 (M) 用以下2种经验公式[17-18]估算:

    $$ M_1{\mathrm{=}}\mathrm{e}^{{\text{−}}0.015\,2{\mathrm{+}}0.654\,3 \times \ln ({K}){\text{−}}0.279 \times \ln \left(L_{\text {inf }}\right){\mathrm{+}}0.463\,4 \times \ln (T)}$$ (5)
    $$M_2{\mathrm{=}}4.118 \times K^{0.73} \times L_{\text {{\mathrm{inf}} }}^{{\text{−}}0.33} $$ (6)

    式中:T为环境平均温度,本研究采用27 ℃[19];取M1M2的平均值,以获得南海北部深水金线鱼自然死亡率(M)相对准确的值。

    统计深水金线鱼繁殖期雌性个体的性腺发育状况,按照I—VI期标准进行划分,性腺发育达到IV期或以上,则认为该个体为性成熟个体[4]。50%性成熟体长(L50)采用最小二乘法 (Least square method, LSM)拟合以下逻辑斯蒂模型得到[20]

    $$ P_i{\mathrm{=}}\frac{1}{1{\mathrm{+}}{\mathrm{exp}} \left[{\text{−}}r\left(L_i{\text{−}}L_{50}\right)\right]} $$ (7)

    式中:Pi为体长组Li对应的成熟个体比例;r为模型参数。95%性成熟体长(L95)由1.1×L50估算[21]

    LBB是近年来新开发的一种根据长度-频率数据来评估鱼类种群开发状况的方法,该方法使用贝叶斯蒙特卡罗马尔科夫链 (MCMC)来计算死亡率参数和相对种群数量,适用于可以终生持续生长的物种[22],该方法可以估算目标物种的渐近体长 (Linf)、50%个体可被渔具捕获的体长值 (Lc)、最适开捕体长 (Lc_opt)、相对自然死亡率 (M/K)、捕捞死亡率相对于自然死亡率 (F/M)、已经开发的生物量相对于未开发的生物量 (B/B0) 和相对种群大小 (B/BMSY,其中MSY为最大可持续产量)。LBB方法只需要长度-频率数据便可进行计算,但也可以通过设置M/KLinfL50等先验参数,以获得更准确的结果[23] (表2)。

    表  2  LBB和LBSPR方法所需数据目录
    Table  2.  Data for LBB and LBSPR methods
    参数
    Parameter
    估算值
    Estimated value
    基于长度的繁殖潜力比模型
    Length-based spawning potential ratio
    (LBSPR) assessment model
    基于长度的贝叶斯估计模型
    Length-based Bayesian biomass
    (LBB) assessment model
    长度-频率数据
    Length–frequency data
    /
    渐近体长 (Linf)
    Asymptotic length/cm
    23.7 cm
    生长系数 (K)
    Growth coefficient
    0.326
    自然死亡率 (M)
    Natural mortality
    0.762
    相对自然死亡率M/K 2.33
    50%性成熟体长 (L50)
    Body length at 50% sexual maturity/cm
    11.76
    95%性成熟体长 (L95)
    Body length at 95% sexual maturity/cm
    12.94
    渐近体长的变异系数
    Coefficient of variation of Linf
    0.1
    长度-体质量关系参数 (a)
    Length-body mass relationship parameter
    6.854 0×10−5
    长度-体质量关系参数 (b)
    Length-body mass relationship parameter
    2.797 4
    注:“√” 代表必需数据;“〇”代表可选数据。 Note: "√" represents required data; "〇" represents optional data.
    下载: 导出CSV 
    | 显示表格

    根据B/BMSYLc/Lc_opt的评估结果,确定深水金线鱼的种群开发现状。如果一个种群的B/BMSY 1,可认为该种群处于健康状态;当0.8 ≤ B/BMSY < 1时,认为该种群处于轻度过度开发状态;当0.5 B/BMSY < 0.8时,认为该种群正处于过度捕捞状态;当0.2≤B/BMSY < 0.5时,认为该种群处于严重过度捕捞状态;当B/BMSY < 0.2时,认为该种群资源已经崩溃[24]。如果估计的Lc/Lc_opt< 1,认为种群遭受了“生长过度捕捞”(当鱼的平均捕捞尺寸小于每次捕捞可产生最大产量的尺寸时),其中Lc_opt为最适开捕长度,表示可保持种群处于健康状态且捕捞量接近MSY的体长值。LBB方法在2018年由Forese等[25]提出,本研究使用的模型R语言代码(LBB-33a)来自基尔亥姆霍兹海洋研究中心 (http://oceanrep.geomar.de/44832/)。

    在基于长度-频率的数据缺乏分析法中,繁殖潜力比 (Spawning potential ratio, SPR) 是最常用的参考点之一,可以为数据缺乏的渔业管理提供具有决策性的建议[26],SPR定义为在一定捕捞压力下(F>0),当前种群与未被捕捞种群的繁殖潜力比值[27]。繁殖潜力指在不考虑密度制约的情况下,补充量群体的潜在产卵量。LBSPR是基于SPR这一生物学参考点建立的,广泛应用于数据缺乏种群的资源评估[21,28-29]。LBSPR假设被捕捞鱼类种群的SPR是相对捕捞死亡率 (F/M) 和2个生活史比率 (M/KL50/Linf) 的函数,与LBB不同的是,LBSPR的运行除长度-频率外还需要其他参数 (即M/KLinfL50L95,如表2所示),这些参数的估计可通过R软件包TropFishR、M的经验公式和逻辑斯蒂模型计算获得。

    在LBSPR模型的运行过程中,使用最大似然估计法来估算50%和95%选择性下的平均体长($S_{L_{50}} $和$S_{L_{95}} $)和相对捕捞死亡率(F/M),若SPR≤0.2,说明该种群繁殖潜力比过低,种群资源濒临崩溃;当0.2<SPR≤0.35时,说明种群处于充分开发状态;当0.35<SPR<0.4时,则说明种群处于MSY水平[21,26,29]。LBSPR方法在2015年由Hordyk等[30]提出,本研究使用R语言LBSPR软件包完成(https://mirrors.pku.edu.cn/CRAN/web/packages/LBSPR/index.html)。

    本文对2个评估模型的LinfM/KL50 3种参数进行不确定性分析,以判断不同先验参数在不同区间的敏感性。LinfM/KL50根据LBB对3种参数所估算的95%置信区间的数值进行上下浮动,来保证先验信息相对准确,其中Linf分为6组 (数值如表3中第2—第7组所示),体长数值间隔为2 mm;L50分为6组 (数值如表3中第8—第13组所示),体长数值间隔为2 mm;M/K分为6组 (数值如表3中第14—第19组所示),数值间隔为0.02。

    表  3  LBB和LBSPR在不同先验参数下的评估结果
    Table  3.  Estimated results of LBB and LBSPR parameters with different priors
    组号
    Group No.
    先验参数
    Prior parameter
    LBSPR结果
    Results of LBSPR
    LBB结果
    Results of LBB
    Linf/mm L50/mm M/K F/M SPR B/BMSY F/M Lc/Lc_opt
    1 237 117.6 2.33 1.71 0.19 0.89 1.00 0.85
    2 232 1.73 0.20 0.99 0.90 0.87
    3 234 1.76 0.20 0.94 0.97 0.86
    4 236 1.71 0.19 0.90 1.00 0.85
    5 238 1.84 0.18 0.87 1.10 0.84
    6 240 1.73 0.19 0.84 1.10 0.83
    7 242 1.94 0.17 0.81 1.10 0.82
    8 237 112.6 2.33 1.71 0.21 0.89 1.00 0.85
    9 114.6 1.71 0.21 0.89 1.00 0.85
    10 116.6 1.71 0.20 0.89 1.00 0.85
    11 118.6 1.71 0.19 0.89 1.00 0.85
    12 120.6 1.71 0.18 0.89 1.00 0.85
    13 122.6 1.71 0.17 0.89 1.00 0.85
    14 237 117.6 2.28 1.77 0.19 0.84 1.10 0.83
    15 2.30 1.75 0.19 0.87 1.00 0.84
    16 2.32 1.73 0.19 0.87 1.00 0.84
    17 2.34 1.70 0.20 0.89 1.00 0.84
    18 2.36 1.68 0.20 0.92 1.00 0.86
    19 2.38 1.66 0.20 0.92 0.99 0.86
    下载: 导出CSV 
    | 显示表格

    2014—2019年采集的深水金线鱼样品体长范围为6.5~20.0 cm,体质量范围为9.3~202.0 g。雄性、雌性和整体的体长-体质量关系 (图2) 分别为:

    图  2  南海北部雄性、雌性、整体深水金线鱼体长-体质量关系
    Figure  2.  Length-body mass relationship among male, female and all N. bathybius in northern South China Sea
    $$ \text { 雄性}: W{\rm{=}}6.911\,8 \times 10^{{\text{−}}5} L^{2.794\,8}\left(r^2{\rm{=}}0.920\,1\right) $$ (8)
    $$ \text { 雌性 } :W{\mathrm{=}}8.211\,0 \times 10^{{\text{−}}5} L^{2.761\,0}\left(r^2{\mathrm{=}}0.911\,5\right) $$ (9)
    $$ \text { 整体}: W{\mathrm{=}}6.854\,0 \times 10^{{\text{−}}5} L^{2.797\,4}\left(r^2{\mathrm{=}}0.920\,3\right) $$ (10)

    在R软件包TropFishR中使用ELEFAN方法提供的von Bertalanffy生长参数估计,评估的Linf为23.7 cm,K为0.326,t0为 −0.54,估算的生长性能指数(Φ)为2.26。von Bertalanffy生长函数为:

    $$ L_t{\mathrm{=}}23.7 \times\{{\text{−}}{\mathrm{exp}} [{\text{−}}0.33 \times(t{\mathrm{+}}0.54)]\} $$ (11)

    使用2种经验公式估算的深水金线鱼自然死亡率M分别为0.908和0.616,本研究使用2个结果的平均值0.762。因此,LBB和LBSPR的先验信息参数M/K为2.33 (表2)。

    根据非线性回归拟合的LSM 逻辑斯蒂曲线 (图3),估算的深水金线鱼50%性成熟体长(L50)为11.76 cm,95%性成熟体长(L95)为1.1×L50=12.94 cm。

    图  3  南海北部深水金线鱼性成熟比例与体长拟合的逻辑斯蒂曲线
    Figure  3.  Logistic model fitted for relationship between body length and percentage of mature N. bathybius in northern South China Sea

    根据LBB评估结果(图4)显示,南海北部深水金线鱼的渐近体长Linf估算值为23.8 cm (95%置信区间为23.4~24.2 cm),相对自然死亡率M/K为2.35 (95%置信区间为2.22~2.49),相对生物量B/B0为0.31 (95%置信区间为0.23~0.38),种群资源量水平B/BMSY为0.89 (95%置信区间为0.66~1.10),Lc/Lc_opt为0.85,表明深水金线鱼处于轻度过度开发状态和生长性过度捕捞状态,估算的相对捕捞死亡率F/M为1.00 (95%置信区间为0.81~1.20),表明捕捞死亡率和自然死亡率基本一致,开发率 (E) 可由F/(F+M)计算,结果为0.5,Lc_opt为11.1 cm,捕获的约60%的深水金线鱼个体超过了这一长度。

    图  4  南海北部深水金线鱼基于LBB方法的分析结果
    注:左图为模型对体长数据的拟合曲线,右图为LBB方法的预测;Lc是50%个体被渔具捕获的体长;Lopt为未开发种群获得最大生物量的捕捞体长。
    Figure  4.  Estimated values for N. bathybius in northern South China Sea by length-based Bayesian biomass assessment model
    Note: The left figure shows the fits of the model to length data, and the right figure shows the predictions of the LBB analysis; Lc. Length of 50% of the individuals captured by the gear; Lopt. Fishing body length for obtaining maximum biomass in undeveloped populations.

    根据LBSPR的输出结果 (图5),南海北部深水金线鱼的相对死亡率F/M为1.71,开发率E为0.63,高于LBB计算出的开发率,SPR值为0.19,表明深水金线鱼处于过度捕捞状态。选择性和成熟度曲线显示性成熟时的长度低于首次捕获时的长度 (图5-b)。在50%和95%选择性 ($S_{L_{50}} $和$S_{L_{95}} $)下,预估的体长分别为9.9 和11.9 cm (图5-c)。假设本次的长度-频率数据处于稳定状态[26],目前南海深水金线鱼的SPR值为0.19 (图5-c),与其阈值水平相比 (SPR=0.4),该种群目前所捕获样品比预期的个体小,尤其是大体长个体数量明显低于预期 (图5-d)。

    图  5  南海北部深水金线鱼基于LBSPR方法的分析结果
    注:a. 长度-频率分布(柱状图),预测的种群大小组成(折线图);b. LBSPR模型拟合的性成熟和选择性曲线,其中L50=11.76 cm和L95=12.94 cm时;c. 选择性参数($S_{L_{50}}$和 $S_{L_{95}} $)、捕捞死亡率与自然死亡率的比值(F/M)和繁殖潜力比的估计值;d. 在目标SPR=0.4时,观测到的长度-频率数据与预期的大小组成。
    Figure  5.  Length-based spawning potential ratio for N. bathybius in northern South China Sea
    Note: a. Length-frequency distribution (Pillars), and predicted fished size composition (Line chart); b. Maturity and selectivity curves from the fitted LBSPR model when L50 was 11.76 cm and L95 was 12.94 cm; c. Distribution of mean selectivity parameters ($S_{L_{50}} $ and $S_{L_{95}} $), fishing mortality to natural mortality (F/M), and spawning potential ratio; d. Observed length-frequency data against an expected size composition at a target SPR of 0.4.

    LBSPR和LBB 模型的估算结果对参数Linf的设置极为敏感,对参数M/K的设置较为敏感,对参数L50的设置不敏感 (表3)。从第2—第6组Linf的数值变化2 mm,即可引起2种方法估算所得的各个结果的明显变化,第1、第4、第5组参数Linf的数值相差1 mm,各个估算结果同样变化明显,可判断LBSPR和LBB模型的估算结果对参数Linf的设置极为敏感。根据第8—第13组计算结果、第8组和第13组参数L50差值为10 mm,LBB的估算结果并无明显变化,表明LBB对参数L50的设置不敏感;LBSPR的结果则出现了2种情况,根据第8—第13组计算结果,F/M的估算结果未随参数L50的变化发生明显改变,SPR的值在参数L50变化为2 mm条件下变化较小,根据第1、第9、第12组参数L50变化为3 mm条件下的结果对照来看,SPR计算结果变化较为明显,SPR对参数L50设置较为敏感,2种模型的另外4个估算结果无明显变化,LBSPR和LBB 模型的估算结果对参数L50的设置不敏感。根据第14—第19组参数M/K的数值相差0.02的设定,结果表明,2种模型的估算结果变化较为明显,根据第1、第16、第17组,参数M/K变化为0.01的情况下,LBSPR模型估算的F/M结果和LBB模型估算的B/BMSY结果变化明显,另外3个值的估算结果变化小;根据第1、第14、第19组,参数M/K的变化为0.05的情况下,LBSPR和LBB模型的估算结果产生了明显变化,这2种模型的估算结果对参数M/K设置的变化比较敏感。通过对先验参数的不确定性分析,发现LBB和LBSPR模型对参数Linf设置的变化极为敏感,对参数M/K设置的变化比较敏感,对参数L50设置的值不敏感,故在使用这2种基于长度的数据缺乏模型时对参数LinfM/K应谨慎设置。

    本研究使用LBB和LBSPR 2种基于体长的数据缺乏评估方法,对2014—2019年南海北部深水金线鱼进行评估,LBB模型评估结果显示种群资源量水平(B/BMSY)、Lc/Lc_optE分别为0.89、0.85和0.5,LBSPR模型评估结果显示繁殖潜力比 (SPR)为0.19,E为0.63。一般认为SPR=0.4是比较合理的目标参考点,种群在此水平下可长期维持可持续发展,而SPR=0.2常被设定为最低限制参考点,低于此值则表明种群遭受过度捕捞,种群补充量不足,难以维持种群稳定,必须采取严格的管理措施来恢复种群资源量[26,31-32],LBB模型估算结果表明深水金线鱼正处于过度捕捞状态。Gulland[33]提出鱼类的最适开发率E为0.5,大于0.5为过度开发,2种方法的评估结果略显差异。同一种群,使用不同的数据缺乏评估方法可能产生不同的结果。例如,Xu等[34]使用LBB和LBSPR 2种基于体长的数据缺乏评估方法,对2012—2014年、2019—2021年南海细鳞圆鲹 (Decapterus macarellus) 2个种群进行资源评估,2种方法的结果均表明2019—2021年该种群已经被过度捕捞,而这2种方法在对2012—2014年南海细鳞圆鲹种群状况评估中出现了分歧:LBB模型评估结果显示为健康状态,而LBSPR模型评估结果显示为过度捕捞状态。Liao等[35]结合基于产量 (Monte Carlo Catch Maximum Sustainable Yield Model, CMSY)、长度 (LBB)和丰度指数 (Abundance-based method, AMSY) 3种数据缺乏评估方法对2005年南大西洋长鳍金枪鱼(Thunnus alalunga)进行资源评估,基于产量(CMSY)和基于长度(LBB) 的模型估算结果表明在2005年南大西洋长鳍金枪鱼种群处于健康(从过度捕捞状态中恢复健康)状态,而基于丰度指数 (AMSY) 的模型估算结果表明该种群处于过度捕捞状态。本研究中,LBB模型计算结果表明南海北部深水金线鱼处于轻度过度开发状态和生长性过度捕捞状态,而LBSPR模型计算结果表明其正处于过度捕捞状态。通过结合不同模型,能更完整地展现出种群的实际情况,为渔业管理提供更准确的科学依据。

    大多数鱼类的表型特征在经过较强的捕捞压力胁迫后会产生适应性响应,其中个体规格的变化较为明显[36]。根据陈作志等[4]对北部湾深水金线鱼种群参数的研究 (表4,本研究的种群参数在表4中记为2017年),在持续多年的捕捞胁迫下,深水金线鱼的群体结构发生了较为明显的变化,渐近体长(Linf)和50%性成熟体长(L50)在减小后均有所回升,拐点分别为1997年和2009年,Linf从1992年242 mm减小至1997年220 mm后,1997—2009年增长到235 mm,2009—2017年略有增长,为237 mm;L50在1992—2009年呈减小趋势,从128.6 mm减小至110.5 mm,2009—2017年增长到117.6 mm。生长系数(K)在1992—2009年呈增大趋势,2009—2017年K减小至0.326,生长加快的现象有所缓解,这一变化趋势与北部湾蓝圆鲹 (D. maruadsi)[37]和南海北部带鱼(Trichiurus japonicus)[38]类似,这也与陈作志等[4]提到的环境胁迫理论相同,即过大的捕捞强度使得资源群体的密度稀疏,剩余群体的饵料丰度和易得性大幅增加,从而改善其生长环境,使得生长加快,这种变化在停止捕捞或捕捞强度下降时会出现逆转。

    表  4  不同时期深水金线鱼种群参数
    Table  4.  Population parameters of N. bathybius in different periods
    参数
    Parameter
    年份 Year
    1992 1997 2009 2017
    渐近体长 Asymptotic length (Linf)/mm 242 220 235 237
    生长系数Growth coefficient (K) 0.39 0.45 0.48 0.326
    50%性成熟体长Body length at 50% sexual maturity (L50)/mm 128.6 116.5 110.5 117.6
    下载: 导出CSV 
    | 显示表格

    Liao等[35]研究表明,在使用数据缺乏的方法时,应注意模型参数先验的设置和不同类型数据中隐含的不同趋势。对于LBSPR模型,Linf估计的不确定性反映在估计的F/M和SPR上,其中Linf输入越高,F/M越高,SPR越低,反之亦然(表3)[39]

    LBSPR模型结果的准确性依赖于以下主要假设能否满足[26]:1) 渔具选择性曲线符合逻辑斯蒂型选择性曲线;2) 鱼类在整个生命周期中体型持续增加;3) 体长频率数据具有良好的代表性,能代表所在海区内种群的体长结构特征。

    LBB 是一种简单快速的估算鱼类相对种群参数的方法,与其他独立的渔业资源评估方法相比,它不需要年龄、性成熟度、补充量、生长和死亡参数、捕捞努力量等信息,只需要能代表评估对象整个世代的长度频率数据。长度频率数据是渔业生物学测定中最易获取和测量误差最小的数据,不受野外船只晃动和内脏有无的影响[40]。在长度-频率能代表资源各开发阶段的长度组成时,LBB 估算的结果能较好地反映其真实情况[41]。Zhang等[11,42]研究表明LBB模型估计参数(LinfLc/Lc_optZ/KB/B0B/BMSY)对样本量不敏感,电子体长频率分析 (ELEFAN) 方法的先验信息对LBB法有效。与ELEFAN相比,LBB可以估计更多的相关结果 (如B/BMSYLc/Lc_opt)。因此,建议使用ELEFAN和LBB方法来拟合数据贫乏的鱼类种群的长频数据,因为它们在估计管理参考点方面是互补的[11]。根据表3,通过对先验参数的不确定性分析,发现LBB和LBSPR模型对参数Linf设置的变化极为敏感,对参数M/K设置的变化较为敏感,对参数L50设置的值不敏感,故在使用这2种基于长度的数据缺乏模型时对LinfM/K 2种参数应谨慎设置。本研究只针对M/KL50Linf 3个参数对LBB和LBSPR的结果影响进行分析,而渔业数据的变化多种多样,在数据采集时,并非只有一种先验参数会发生变化,在以后的研究中将进一步探讨不同误差的情况,为数据有限条件下的渔业资源评估方法优化提供依据。

  • 图  1   试验期内海水盐度与温度变化

    Figure  1.   Change in sea water salinity and temperature during experiment

    图  2   不同底质组富齿岩虫体长、体宽、体质量和特定生长率随时间变化的情况

    Figure  2.   Body length, body width, body mass and special growth rate of M. maxidenticulata with different substrates along with time

    表  1   北海养殖与乳山野生的富齿岩虫基本营养组成 (干质量, N =30)

    Table  1   Basic nutritional composition of M. maxidenticulata cultured in Beihai and wild M. maxidenticulata in Rushan (dry mass) %

    组别
    Group
    水分
    Moisture
    灰分
    Ash
    粗蛋白质
    Crude protein
    粗脂肪
    Crude fat
    北海养殖 Cultured in Beihai 76.09±0.01a 7.31±0.05a 60.96±0.01a 24.13±0.01a
    乳山野生 Wild in Rushan 74.28±0.01b 9.75±0.50b 64.24±0.01b 19.71±0.01b
    注:同列中标有不同小写字母者表示两组间有显著性差异 (P<0.05)。 Note: The values with different letters within the same column are very significantly different (P<0.05).
    下载: 导出CSV

    表  2   北海养殖与乳山野生的富齿岩虫氨基酸组成

    Table  2   Amino acid compositionof M. maxidenticulata cultured in Beihai and wild M. maxidenticulata in Rushan %

    氨基酸 
    Amino acids 
    北海养殖
    Cultured in Beihai
    乳山野生
    Wild in Rushan
    天冬氨酸# Asp 5.39±0.23a 5.14±0.01a
    谷氨酸# Glu 7.12±0.21a 7.09±0.00a
    甘氨酸# Gly 4.82±0.06a 4.87±0.02a
    丙氨酸# Ala 3.19±0.10a 2.95±0.02a
    丝氨酸 Ser 2.03±0.03a 2.01±0.02a
    酪氨酸 Tyr 1.84±0.06a 1.79±0.00a
    脯氨酸 Pro 1.72±0.03a 1.78±0.06a
    半胱氨酸 Cys 1.06±0.01a 1.02±0.04a
    组氨酸** His 1.48±0.33a 1.08±0.02a
    精氨酸** Arg 3.82±0.10a 3.63±0.02a
    苏氨酸* Thr 2.15±0.05a 2.13±0.01a
    缬氨酸* Val 2.43±0.09a 2.38±0.02a
    蛋氨酸* Met 1.00±0.00a 0.75±0.01b
    异亮氨酸* Ile 2.20±0.09a 2.18±0.01a
    亮氨酸* Leu 3.81±0.14a 3.84±0.00a
    赖氨酸* Lys 4.24±0.22a 3.90±0.04a
    苯丙氨酸* Phe 2.04±0.04a 2.06±0.01a
    色氨酸* Trp 0.50±0.02a 0.50±0.05a
    氨基酸总量 TAA 49.76±1.75a 48.09±0.04a
    必需氨基酸 EAA 18.37±0.67a 17.75±0.01a
    半必需氨基酸 SEAA 5.30±0.42a 4.72±0.00a
    非必需氨基酸 NEAA 26.09±0.66a 25.63±0.03a
    呈味氨基酸 DAA 20.51±0.60a 20.05±0.05a
    EAA/TAA 36.91±0.47a 36.91±0.02a
    SEAA/TAA 10.64±0.54a 9.80±0.02a
    NEAA/TAA 52.46±0.24a 53.29±0.14a
    DAA/TAA 41.23±0.24a 41.69±0.14a
    注:#. 呈味氨基酸;*. 必需氨基酸;**. 半必需氨基酸;同行中标有不同小写字母者表示组间有显著性差异 (P<0.05);表4同此。 Note: #. Delicious amino acid (DAA); *. Essential amino acid (EAA); **. Semi-essential amino acid (SEAA); the values with different letters within the same row are very significantly different from each other (P<0.05). The same case in Table 4.
    下载: 导出CSV

    表  3   北海养殖与乳山野生的富齿岩虫氨基酸营养评价

    Table  3   Amino acids nutritional evaluation of M. maxidenticulata cultured in Beihai and wild M. maxidenticulata in Rushan

    必需氨基酸
    Essential amino acid
    FAO/WHO 评分模式
    FAO score mode/(mg·g−1)
    鸡蛋蛋白
    Egg protein/(mg·g−1)
    北海养殖
    Cultured in Beihai
    乳山野生
    Wild in Rushan
    氨基酸评分
    AAS
    化学评分
    CS
    氨基酸评分
    AAS
    化学评分
    CS
    赖氨酸 Lys 340 441 1.28 0.99 1.12 0.86
    缬氨酸 Val 310 411 0.80 0.61 0.75 0.56
    亮氨酸 Leu 440 534 0.89 0.73 0.85 0.70
    异亮氨酸 Ile 250 331 0.90 0.68 0.85 0.64
    苏氨酸 Thr 250 292 0.88 0.75 0.83 0.71
    苯丙氨酸+酪氨酸 Phe+Tyr 380 565 1.05 0.70 0.98 0.66
    色氨酸 Trp 60 99 0.85 0.52 0.82 0.50
    必需氨基酸指数 EAAI 69.91 65.28
    下载: 导出CSV

    表  4   北海和乳山野生的富齿岩虫脂肪酸含量

    Table  4   Fatty acid compositionof M. maxidenticulata cultured in Beihai and wild M. maxidenticulata in Rushan %

    脂肪酸
    Fatty acids
    北海养殖
    Cultured in Beihai 
    乳山野生
    Wild in Rushan
    豆蔻酸 C14:0 1.15±0.02a 0.97±0.03a
    十五烷酸 C15:0 0.44±0.12a 0.40±0.06a
    棕榈酸 C16:0 15.96±0.16a 14.73±0.58a
    十七碳酸 C17:0 0.95±0.01a 1.22±0.04a
    硬脂酸 C18:0 5.48±0.05a 4.93±0.02a
    花生酸 C20:0 0.27±0.00a 0.15±0.01a
    二十一烷酸 C21:0 1.81±0.02a 1.78±1.62a
    山嵛酸 C22:0 0.44±0.01a 0.60±0.01a
    二十三烷酸 C23:0 3.83±0.01a 6.74±0.16b
    木蜡酸 C24:0 0.14±0.00a 0.16±0.01a
    饱和脂肪酸总量 ∑SFA 30.47±0.05a 31.70±2.53a
    棕榈一烯酸 C16:1 2.88±0.01a 2.42±0.12a
    十七碳一烯酸 C17:1 ND 0.52±0.23
    油酸 (反) C18:1 n-9t 1.00±0.04a 0.22±0.06b
    油酸 (顺) C18:1 n-9c 20.42±0.07a 13.84±0.17b
    花生一烯酸 C20:1 n-9 4.50±0.02a 4.35±0.16a
    芥子酸 C22:1 n-9 0.43±0.00a 1.99±0.06b
    神经酸 C24:1 n-9 4.19±0.02a 7.20±0.16b
    单不饱和脂肪酸总量 ∑SUFA 33.41±0.06a 30.54±0.85b
    亚油酸 (反) * C18:2 n-6t 0.10±0.04a 0.19±0.04a
    亚油酸 (顺) * C18:2 n-6c 15.59±0.03a 6.53±0.43b
    α-亚麻酸* C18:3 n-6 0.19±0.00a 4.80±4.64b
    γ-亚麻酸# C18:3 n-3 0.34±0.00a 0.85±0.03a
    花生二烯酸* C20:2 9.53±0.05a 7.76±0.28b
    花生三烯酸* C20:3 n-6 1.14±0.00 ND
    花生四烯酸* C20:4 n-6 3.29±0.02a 3.93±0.10a
    二十二碳二烯酸* C22:2 0.13±0.00a 2.07±0.13b
    花生五烯酸# EPA C20:5 n-3 1.71±0.00a 2.31±0.05a
    二十二碳五烯酸# DPA 1.22±0.00a 4.40±0.08b
    二十二碳六烯酸# DHA C22:6 n-3 2.88±0.00a 4.92±0.12b
    多不饱和脂肪酸总量 ∑PUFA 36.12±0.00a 37.76±3.38a
    ∑n-3 PUFA 6.14±0.01a 12.48±0.27b
    ∑n-6 PUFA 29.98±0.01a 25.28±3.65b
    EPA+DHA 4.58±0.00a 7.23±0.17b
    ∑n-3 PUFA/∑n-6 PUFA 0.20±0.00a 0.51±0.08b
    (EPA+DHA)/∑n-3 PUFA 0.75±0.00a 0.58±0.00b
    (EPA+DHA)/∑n-6 PUFA 0.15±0.00a 0.29±0.05b
    注:*. n-6系列多不饱和脂肪酸;#. n-3系列多不饱和脂肪酸;ND. 未测出。 Note: *. n-6 poly unsaturated fatty acids (PUFA); #. n-3 poly unsaturated fatty acids (PUFA); ND. Not detected.
    下载: 导出CSV
  • [1]

    WANG Z, ZHANG Y, QIU J W. A new species in the Marphysa sanguinea complex (Annelida, Eunicidae) from Hong Kong[J]. Zool Stud, 2018, 57(48): 1-13.

    [2]

    MOLINA-ACEVEDO I C, CARRERA-PARRA L F. Reinstatement of three species of the Marphysa sanguinea complex (Polychaeta: Eunicidae) from the Grand Caribbean Region[J]. Zootaxa, 2015, 3925(1): 37-55. doi: 10.11646/zootaxa.3925.1.3

    [3]

    LAVESQUE N, DAFFE G, BONIFÁCIO P, et al. A new species of the Marphysa sanguinea complex from French waters (Bay of Biscay, NE Atlantic) (Annelida, Eunicidae)[J]. ZooKeys, 2017, 716(10): 1-17.

    [4]

    ZANOL J, SILVA T S C, HUTCHINGS P. Marphysa (Eunicidae, Polychaete, Annelida) species of the Sanguinea group from Australia, with comments on pseudo-cryptic species[J]. Invertebr Biol, 2016, 135(4): 328-344. doi: 10.1111/ivb.12146

    [5]

    LAVESQUE N, DAFFE G, GRALL J, et al. Guess who? On the importance of using appropriate name: case study of Marphysa sanguinea (Montagu, 1813)[J]. ZooKeys, 2019, 2(859): 1-15.

    [6]

    PREVEDELLI D, N'SIALA G M, SIMONINI A R. Life cycle of Marphysa sanguinea (Polychaeta: Eunicidae) in the Venice Lagoon (Italy)[J]. Mar Ecol, 2010, 28(3): 384-393.

    [7] 毛昕. 岩虫 (Marphysa sanguinea) 消化道形态学、组织学和组织化学的研究[D]. 青岛: 中国海洋大学, 2004: 3-5.
    [8] 杨大佐, 周一兵, 陈爱华, 等. 岩虫室内人工繁育的初步研究[J]. 水产科学, 2011, 30(9): 572-574. doi: 10.3969/j.issn.1003-1111.2011.09.012
    [9]

    THI THU E V, RAHMAN M, PHOO W W, et al. Salinity effects on growth and survival of the polychaete rockworm Marphysa sanguinea (Montagu, 1813) juveniles and adults[J]. J Aquacult Res Dev, 2019, 10(2): 1-7.

    [10]

    PREVEDELLI D, SIMONINI R. Life cycles in brackish habitats: adaptive strategies of some polychaetes from the Venice lagoon[J]. Oceanologica Acta, 2003, 26(1): 77-84. doi: 10.1016/S0399-1784(02)01232-X

    [11]

    LI S L, CHEN Y, ZHANG M, et al. Complete mitochondrial genome of the marine polychaete, Marphysa sanguinea (Polychaeta, Eunicida)[J]. Mitochondr DNA, 2014, 27(1): 42-43.

    [12]

    ZANOL J, HALANYCH K M, STRUCK T H, et al. Phylogeny of the bristle worm family Eunicidae (Eunicida, Annelida) and the phylogenetic utility of noncongruent 16S, COI and 18S in combined analyses[J]. Mol Phylogenet Evol, 2010, 55(2): 660-676. doi: 10.1016/j.ympev.2009.12.024

    [13]

    ZHAO H, WANG Y X, YANG D Z, et al. An analysis of genetic diversity in Marphysa sanguinea from different geographic populations using ISSR polymorphisms[J]. Biochem Sys Ecol, 2016(64): 65-69.

    [14] 张文龙, 刘声宇, 刘晓威, 等. 温度和体质量对岩虫碳收支的影响[J]. 大连海洋大学学报, 2018, 33(2): 95-99.
    [15]

    YANG D, CHEN F, ZHOU Y, et al. Diel variation in metabolism and ammonia excretion of Marphysa sanguinea (Polychaeta: Eunicidae)[J]. Chin J Oceanol Limnol, 2016, 34(6): 1209-1217. doi: 10.1007/s00343-016-4340-x

    [16]

    JANG M, SHIM W J, HAN G M, et al. Formation of microplastics by polychaetes (Marphysa sanguinea) inhabiting expanded polystyrene marine debris[J]. Mar Poll Bull, 2018, 131: 365-369.

    [17]

    MDAINI Z, CAFSI M E, TREMBLAY R, et al. Spatio-temporal variability of biomarker responses and lipid composition of Marphysa sanguinea, Montagu (1813) in the anthropic impacted lagoon of Tunis[J]. Mar Poll Bull, 2019, 144: 275-286. doi: 10.1016/j.marpolbul.2019.04.065

    [18]

    LIU Y B, HUTCHINGS P, KUPRIYANOVA E. Two new species of Marphysa quatrefages, 1865 (Polychaeta: Eunicida: Eunicidae) from northern coast of China and redescription for Marphysa orientalis Treadwell, 1936[J]. Zootaxa, 2018, 4377(2): 191-215. doi: 10.11646/zootaxa.4377.2.3

    [19]

    LEELATANAWIT R, UAWISETWATHANA U, KHUDET J, et al. Effects of polychaetes (Perinereis nuntia) on sperm performance of the domesticated black tiger shrimp (Penaeus monodon)[J]. Aquaculture, 2014, 433(1): 266-275.

    [20]

    LEELATANAWIT R, UAWISETWATHANA U, KLANCHUI A, et al. Transcriptomic analysis of male black tiger shrimp (Penaeus monodon) after polychaete feeding to enhance testicular maturation[J]. Mar Biotechnol, 2017, 19(2): 125-135. doi: 10.1007/s10126-017-9738-8

    [21]

    LYTLE J S, LYTLE T F, OGLE J T. Polyunsaturated fatty acid profiles as a comparative tool in assessing maturation diets of Penaeus vannamei[J]. Aquaculture, 1990, 89(3): 287-299.

    [22]

    FAO/WHO. Energy and protein requirements[R]. Rome: FAO, 1973: 52.

    [23] 中国预防医学科学院营养与食品卫生研究所. 食物成分表(全国代表值)[M]. 北京: 人民卫生出版社, 1991: 30-82.
    [24]

    PELLETT P L, YONG V R. Nutritional evaluation of protein foods[M]. Tokyo: United National University Press, 1980: 26-29.

    [25] 毛东东, 张凯, 欧红霞, 等. 投喂皇竹草和配合饲料对草鱼生长及肌肉营养成分的影响[J]. 大连海洋大学学报, 2018, 33(1): 7-13.
    [26] 马双, 郝淑贤, 李来好, 等. 几种鱼卵营养成分对比分析[J]. 南方水产科学, 2019, 15(4): 113-121. doi: 10.12131/20180266
    [27] 李钧, 林祥志, 谢双如, 等. 贝类养殖技术之一皱纹盘鲍南北接力养殖技术研究[J]. 中国水产, 2007(11): 46-47. doi: 10.3969/j.issn.1002-6681.2007.11.027
    [28] 蔡含筠, 马明辉. 岩虫的个体发育[J]. 辽宁师范大学学报(自然科学版), 1994, 17(4): 328-335.
    [29] 于海志. 岩虫 (Marphysa Sanguinea) 的性腺发育和生殖周期[D]. 青岛: 中国海洋大学, 2004: 5.
    [30] 吕峰, 巩杰, 倪建忠, 等. 不同饵料投喂对南美白对虾亲虾性腺成熟及E75基因表达的影响[J]. 黑龙江畜牧兽医, 2018(8): 198-202.
    [31] 姚卫军, 黄翔鹄, 李活. 不同天然饵料对凡纳滨对虾亲虾性腺发育的影响[J]. 广东海洋大学学报, 2009, 29(4): 84-88. doi: 10.3969/j.issn.1673-9159.2009.04.019
    [32] 刘天红, 于道德, 李红艳, 等. 东营养殖双齿围沙蚕营养成分分析及膳食营养评价[J]. 水产科学, 2017, 36(2): 160-166.
    [33] 曹启猛, 刘清兵, 于跃芹, 等. 山东沿海养殖和野生双齿围沙蚕营养成分比较及饵料安全性评价[J]. 中国水产科学, 2016, 23(5): 1164-1172.
    [34] 汪春玲. 不同投喂频率和饲料中DHA/EPA比例对凡纳滨对虾亲虾繁殖性能及后代质量的影响[D]. 上海: 上海海洋大学, 2018: 6-10.
    [35] 田相利. 变温对中国对虾(Fenneropenaeus chinensis OSBECK)生长的影响及其生物能量学机制[D]. 青岛: 中国海洋大学, 2004: 2-6.
    [36] 段志刚, 吴金英, 李文笙, 等. 低温对罗非鱼类影响的相关研究进展[J]. 南方水产科学, 2011, 7(6): 79-84.
    [37] 孔祥会, 王桂忠, 李少菁. 低温驯化锯缘青蟹鳃抗氧化防护、ATPase及膜脂肪酸组成变化[J]. 水生生物学报, 2007, 31(1): 61-68.
图(2)  /  表(4)
计量
  • 文章访问数:  575
  • HTML全文浏览量:  164
  • PDF下载量:  37
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-10-21
  • 修回日期:  2021-04-19
  • 录用日期:  2021-05-29
  • 网络出版日期:  2021-09-09
  • 刊出日期:  2021-12-05

目录

/

返回文章
返回