Effects of morphological traits on body mass of Changfeng silver carp (Hypophthalmichthys molitrix) at different ages
-
摘要: 为研究长丰鲢 (Hypophthalmichthys molitrix) 形态特征对生长发育的影响，文章采用主成分分析和通径分析对其常规形态和框架结构进行了研究。结果表明，不同月龄长丰鲢各形态性状与体质量均呈极显著正相关 (P<0.01)；不同阶段体质量与形态性状的相关性有所不同。主成分分析结果显示6和12月龄表现为鱼体整体结构的协同生长发育，而36月龄主要是躯干部的增长。通径分析结果显示3个月龄中鳃盖背侧后缘至鳃盖腹侧后缘 (x8)、全长 (x11) 和体高 (x13) 与体质量的通径系数均达到显著水平 (P<0.05)，这3个性状在各阶段均对体质量有显著影响；同时建立了3个月龄组不同形态性状对体质量的最优回归方程。不同阶段各形态性状对体质量的影响有所不同，建议早期生长阶段将全长、体长和体高作为长丰鲢的主要选育目标，而成鱼阶段除了考虑全长，还需要重视鱼体宽度对体质量的影响。Abstract: To study the influence of morphological characteristics of Changfeng silver carp (Hypophthalmichthys molitrix) on its growth, we applied the principal component analysis and path analysis methods to analyze the conventional morphological measurements and truss network. The results show that the correlation coefficients between the morphological traits and body mass were extremely significant at different ages (P<0.01). Besides, the correlation coefficients between them were also different at different stages. Principal component analysis revealed that 6- and 12-month old H. molitrix showed a synergistic growth, as well as development of the overall structure, but it was mainly the growth of the trunk at 36-month age. Path analysis shows that the path coefficients of three morphological traits (x8, x11, x13) for the body mass were significant at three stages (P<0.05), and these traits had significant impacts on the body mass at all stages. The optimal regression equations of morphological traits to body mass were established for the three groups. The effects of morphological traits on the body mass were distinct at different stages. The results suggest that the total length, body length and body height should be selected as the breeding targets for H. molitrix at early growth stage. However, the target trait for selective breeding of adult fish should mainly be total length, and the synergistic breeding of fish body width should be emphasized.
-
据联合国粮农组织 (Food and Agriculture Organization of the United Nations, FAO) 报告，预计到2050年世界人口将达到90亿至100亿，全球粮食需求将增长70%~100% [1]，人们对富含蛋白质的水产品需求也随之增加。未来，海水养殖业将成为供给水产品的重要行业[2]，在不断变化的气候条件与环境因素下绘制物种水产养殖的潜在适宜区域，对海水养殖业的可持续发展至关重要。随着“向海经济”等战略思想[3]的提出，以及近岸养殖空间的饱和，海水养殖从近岸向离岸发展；同时，基于我国南海丰富的渔业资源，人们开始探索更多的适养品种尤其是金枪鱼等高端养殖品种。金枪鱼是全球公认的具有高营养价值和经济价值的高端海洋鱼类[4]，其消费市场广阔。世界金枪鱼养殖始于20世纪70年代，自1984年起FAO对全球金枪鱼养殖产量进行统计，当年北美洲首次养殖大西洋蓝鳍金枪鱼 (Thunnus thynnus)，其产量为50 t；2020年全球金枪鱼养殖产量约6.8×104 t，增长了1 368倍。金枪鱼适宜在水深50 m以内的海域养殖，可采用深远海网箱养殖方式[5]。因此，确定适宜的金枪鱼海水养殖区域已成为拓展金枪鱼养殖空间的关键问题。
海水养殖区域范围选址一般采用地理信息系统 (Geographic Information System, GIS)[6]和物种分布模型 (Species Distribution Model, SDM)[7]等方法，其影响因子涵盖物理、经济、社会因素等范畴。例如，基于GIS的多准则评价模型，将物理、环境和社会等多种因素纳入结构化模型，以此来确定日本扇贝水产养殖的潜力区域[8]；SDM基于环境生态理论[9]，利用算法来评估和预测物种在地理空间和时间上的分布，已广泛应用于预测物种海水养殖扩张的潜在区域[10]以及气候变化影响下的海水养殖生产潜力变化趋势[11]等。随着海水养殖面积的不断扩大和水产养殖设施的迭代发展，提前谋划布局中国南海金枪鱼类养殖潜在适宜区域具有一定的现实意义。
本研究选取南海金枪鱼类主要物种—黄鳍金枪鱼 (T. albacares) 和鲣 (Katsuwonus pelamis) 为研究对象，结合其出现的分布点和环境变量数据，运用SDM中使用广泛的最大熵 (Maximum Entropy, MaxEnt) 模型预测其潜在适生区的空间分布，并定量分析物种分布与主要环境因素之间的关系，以期为黄鳍金枪鱼和鲣资源的保护和利用提供参考。
1. 材料与方法
1.1 研究区域概况
中国南海海域是我国渔民历代作业的传统渔场，其南部海域鱼类资源丰富，浮游生物种类繁多[12]，也是我国重要的金枪鱼捕捞场所[13]。研究区域为我国十段线以内海域，105°E—124°E、2°N—25°N (图1)。
1.2 数据来源
1.2.1 物种出现数据
南海金枪鱼类部分出现数据从全球公开数据库Global Biodiversity Information Faculty (GBIF, www.gbif.org) 获取，其不仅包括黄鳍金枪鱼和鲣出现的记录次数，还有分布点位的经纬度。从1950年至今，南海共出现黄鳍金枪鱼30个位点、鲣34个位点。为降低采样偏差的影响并保留最大数量的有用信息，使用ArcGIS 10.8软件对数据集进行空间自相关处理，剔除重复或自相关数据，保留了25个黄鳍金枪鱼位点和29个鲣位点的发生数据和经纬度，并结合2010年3月—2013年2月南海金枪鱼探捕和生产监测到的87个黄鳍金枪鱼位点和39个鲣位点数据，最终共使用112个黄鳍金枪鱼位点和68个鲣位点，建立了相应的物种分布模型。
1.2.2 环境数据
影响鱼类自然种群分布模式的重要变量主要有水温、盐度、食物、溶解氧、叶绿素a和海水流速[14]。目前，国内外学者研究黄鳍金枪鱼和鲣渔场分布的主要环境因子包括表层水温 (Sea surface temperature, SST)、叶绿素a (Chl-a) 和表层盐度 (Sea surface salinity, SSS) 等[15-21]。因此，本研究以表层水温、底层水温 (Sea benthic temperature, SBT)、表层盐度、底层盐度 (Sea benthic salinity, SBS)、表层溶解氧 (Sea surface dissolved oxygen, SSDO)、底层溶解氧 (Sea benthic dissolved oxygen, SBDO) 以及叶绿素 a 7个环境因子为变量，建立南海金枪鱼类养殖潜在适生区物种分布模型。各环境因子数据主要下载于Bio-ORACLE v2.1 (https://www.bio-oracle.org/)，以5'×5' 表示经纬度网格：表层水温、表层溶解氧、表层盐度为目前平均值，底层水温、底层盐度、底层溶解氧为目前平均深度平均值；叶绿素 a浓度数据来自于哥白尼海洋环境监测服务中心 (Copernicus Marine Environment Monitoring Service, CMEMS, https://resources.marine. copernicus.eu)。
1.3 SDM
SDM是一种预测物种空间分布的模型，主要用于研究物种的地理分布与环境变量间的关系[22]，被广泛应用于区域生物多样性评估、种群变化、物种迁移、全球变化生物学、野生动物管理[23-24]等，以量化潜在海水养殖适宜区域[7]等方面的研究。SDM的计算方法包括回归分析、MaxEnt 模型、支持向量随机模型、生态位模型、决策树模型、随机森林模型等，各种方法在总结响应和预测变量之间的关系上存在差异；如在将物种分布转移到不同地理位置[25]或时间背景时，不同方法的预测结果可能存在显著性差异。而MaxEnt模型作为物种分布模型中的常用方法之一，其本质是一种机器学习技术，可用来预测和分析复杂的数据集；且基于信息论中的最大熵原理，可有效地提取数据中的模式，并可自动学习和更新模型，以适应新的数据。
1.4 模型构建
本研究运用了SDM中的MaxEnt模型，在0.5°×0.5° 网格上预测南海黄鳍金枪鱼和鲣的养殖潜在适宜性指数 (Suitability index, SI)。首先，使用MaxEnt 3.4.4软件分别将黄鳍金枪鱼和鲣的出现点样本随机分为2个数据集，将75%的记录设置为训练集，剩余25%的数据作为测试集。其次，采用Bootstrap抽样方法，重复运行10次，以消除随机性和重复数值。根据Logistic格式输出的结果，得到响应曲线，并运用环境贡献率与Jackknife检验结果分析主要环境因子对金枪鱼养殖适宜区分布的影响程度，促使选址更加合理科学可靠。受试者工作特征曲线下面积 (Area under the receiver operating characteristic curve, AUC) 检验值作为一种用于评估模型性能的指标，可反映模型的准确性和敏感性；因此，模型的准确性可以利用模型预测得出的AUC值来评估[19]。AUC值的取值范围为0~1，值越大则表明模型预测的准确度越高，0.5表示模型的预测能力为随机猜测水平，1表示具有极高的预测能力。最后，根据SI (0~1) 定量评估海水养殖适宜性，SI值越高，养殖地区的潜力就越大。研究将模型预测的概率结果以阈值为0.38、0.5、0.6进行划分，即SI小于0.38的区域被划分为非养殖区；SI介于0.38~0.5的区域被划分为低养殖区；SI介于0.5~0.6的区域被划分为中养殖区；SI大于0.6的区域被划分为高养殖区。
2. 结果
2.1 模型预测能力评价
最大熵模型检验结果如表1所示，两种金枪鱼的平均AUC检验值均大于0.848，表明黄鳍金枪鱼和鲣的物种分布模型具有较好的预测能力。模型预测的适生性模拟结果可用于进一步分析鱼类丰富度分布，如预测结果优异且性能良好，可用于研究南海金枪鱼类养殖适宜区的分布。
表 1 模型检验结果Table 1. Model test results种类
Species训练数据 (AUC 值)
Training data (AUC value)测试数据 (AUC 值)
Test data (AUC value)AUC 值标准差
Standard deviation of AUC value黄鳍金枪鱼 T. albacares 0.872 3 0.848 1 0.022 4 鲣 K. pelamis 0.880 3 0.872 1 0.028 8 2.2 环境变量重要性及响应曲线
2.2.1 环境变量重要性
Jackknife检验可以用来评估模型的准确性，以及探讨环境因子对金枪鱼栖息地分布的影响程度。根据Jackknife检验可知 (图2-a)，黄鳍金枪鱼在不包含某个环境因子得分中，不包含叶绿素 a浓度、底层溶解氧、表层水温的平均得分高于不包含底层盐度、表层溶解氧、表层盐度和底层水温这4项，这表明叶绿素 a浓度、底层溶解氧和表层水温是影响黄鳍金枪鱼栖息地分布的重要环境因子；而底层盐度、表层溶解氧、表层盐度和底层水温这4项对黄鳍金枪鱼栖息地分布的影响较小。鲣在不包含某个环境因子的得分中，不包含叶绿素 a浓度、底层溶解氧、表层水温的平均得分高于不包含表层盐度、底层盐度、底层水温和表层溶解氧，这表明叶绿素 a浓度、底层溶解氧、表层水温是影响鲣栖息地分布的重要环境因子，而表层盐度、底层盐度、底层水温以及表层溶解氧对鲣栖息地分布的影响相对较小。
![]()
图 2 Jackknife 检验不包含某个环境因子 (a) 或仅包含某个环境因子 (b) 的得分Figure 2. Score of Jackknife test without a certain environmental factor (a) or with only one environmental factor (b)在仅包含某个环境因子的得分中 (图2-b)，黄鳍金枪鱼仅包含表层溶解氧、底层水温和表层水温的平均得分高于仅包含叶绿素 a 浓度、表层盐度、底层溶解氧和底层盐度的得分，表明表层溶解氧、底层水温和表层水温是影响黄鳍金枪鱼栖息地分布的重要环境因子，而其余4种因素对黄鳍金枪鱼栖息地分布的影响较小。鲣中仅包含表层溶解氧、表层水温和表层盐度的平均得分高于仅包含叶绿素 a 浓度、底层水温、底层盐度和底层溶解氧的得分，表明表层溶解氧、表层水温和 表层盐度是影响鲣栖息地分布的重要环境因子，其余4种因素对鲣栖息地分布的影响较小。
结合Jackknife检验结果及7种环境因子的贡献率，分析得出表层溶解氧是影响黄鳍金枪鱼和鲣潜在适宜区分布最重要的环境变量，其环境因子贡献率分别为42.01%和29.87%。由图3可知，黄鳍金枪鱼中环境因子贡献率排名前3的为表层溶解氧、底层盐度、底层水温，累积贡献率占所有环境因子总贡献率的73.02%以上；置换重要值排前3名的分别是底层盐度、底层水温和表层溶解氧，累积贡献率为79.16%。对于鲣而言，环境因子贡献率排名前3的分别是表层溶解氧、底层盐度、表层水温，累积贡献率占所有环境因子总贡献率的69.41%以上；置换重要值排前3名的为底层盐度、表层溶解氧和表层盐度，累积贡献率为73.93%。可以看出，影响黄鳍金枪鱼潜在分布的重要环境变量主要为表层溶解氧、底层水温、底层盐度，而影响鲣潜在分布的重要环境变量主要为表层溶解氧、表层水温、底层盐度。
2.2.2 响应曲线
根据模型预测金枪鱼类潜在分布结果对环境变量的响应情况，黄鳍金枪鱼和鲣对表层溶解氧的响应曲线随算法的不同而变化。如图4-a所示，在表层溶解氧浓度为201.91~205.61 mol·m−3的海域，黄鳍金枪鱼和鲣的出现概率均较高 (大于0.5)。在表层溶解氧浓度大于201.52 mol·m−3时，黄鳍金枪鱼开始出现 (出现概率为0.38) 且随着浓度的增加其出现概率上升，当表层溶解氧浓度为202.30 mol·m−3时出现概率最大 (0.564)，之后随着浓度的增加出现概率呈下降趋势，在浓度为242.68 mol·m−3时达到临界值。因此，适宜黄鳍金枪鱼生长的表层溶解氧浓度介于201.52~242.68 mol·m−3。在表层溶解氧浓度大于200.83 mol·m−3时，鲣开始出现 (出现概率为0.38)且随着浓度的增加其出现概率上升，在浓度为203.07 mol·m−3时出现概率最大 (0.623)，之后随着浓度的增加出现概率呈下降趋势，在浓度为208.35 mol·m−3时达到临界值。由此可知，适宜鲣生长的表层溶解氧浓度介于200.83~208.35 mol·m−3。
![]()
图 4 重要环境因子对金枪鱼适生性指数的响应曲线Figure 4. Response curve of important environmental factors to suitable index of tuna适宜黄鳍金枪鱼生存的底层水温介于1.96~32.61 ℃，在底层水温介于29.80~32.61 ℃时，黄鳍金枪鱼的出现概率最大 (0.904，图4-b)；适宜鲣生存的表层水温介于19.71~28.96 ℃，在表层水温为25.16 ℃时，鲣的出现概率最大 (0.792，图4-c)。
图4-d显示，当底层盐度高于34.37‰ 时，黄鳍金枪鱼开始出现 (出现概率为0.38) ，其出现概率随底层盐度的增加而上升并在34.97‰~35.26‰时最大 (0.665)，因此底层盐度介于34.37‰~35.26‰时适宜黄鳍金枪鱼生长。鲣的底层盐度适宜范围介于34.30‰~35.26‰，在介于34.97‰~35.26‰ 时出现概率最大 (0.782)。
在以上重要环境因子中，黄鳍金枪鱼的最适宜范围分别为：表层溶解氧浓度介于201.52~242.68 mol·m−3，底层水温介于1.96~32.61 ℃，底层盐度介于34.37‰~35.26‰；鲣的最适宜范围为：表层溶解氧浓度介于200.83~208.35 mol·m−3，表层水温介于19.71~28.96 ℃，底层盐度介于34.30‰~35.26‰。
2.3 金枪鱼类养殖SI
在南海海域黄鳍金枪鱼和鲣的SI具有很强的区域变化。根据适宜SI对南海海域各网格单元的养殖概率进行定量评估，分布概率小于0.38为非适宜养殖区，分布概率大于0.6为潜在高适宜养殖区 (图5)。南海渔场与海域名参照《南海渔场作业图集》
1 。根据模型预测结果，黄鳍金枪鱼潜在适宜养殖区从南海海域中西部延伸至北部以及台湾岛东部海域，主要集中于台湾南部渔场的中西部、东沙渔场、海南岛东南部渔场的东南部、中沙东部渔场的中西部、西沙西部渔场的中东部、西中沙渔场、南沙西北部渔场的中东部、南沙中部渔场以及南沙中西部渔场的东北部，其潜在养殖高SI分布区域主要集中在东沙渔场、南沙中部渔场以及南沙中西部渔场的东北部；鲣潜在适宜养殖区从南沙群岛西北部海域延伸至台湾岛西南部，主要集中于台湾南部渔场的中西部、东沙渔场、海南岛东南部渔场的东南部、中沙东部渔场的中西部、西沙西部渔场的中东部、中沙渔场以及南沙西北部渔场的中东部，其潜在养殖高SI分布区域主要集中在台湾南部渔场的中西部以及东沙渔场。2.4 南海金枪鱼类潜在养殖丰富度分布
预测潜在金枪鱼类海水养殖面积因种类而异。在叠加金枪鱼类养殖潜在SI分布图的基础上，进一步绘制其潜在养殖丰富度分布图 (图6)。总体来说，黄鳍金枪鱼和鲣在我国南海海域养殖潜在区域分布较广，其中，台湾南部渔场的中西部、东沙渔场以及南沙中部渔场是其分布最为集中的区域。黄鳍金枪鱼和鲣在南海海域分布的高SI区域主要集中在台湾南部渔场的中西部海域。相比于黄鳍金枪鱼，鲣的潜在养殖区域面积较小，但处于潜在高适宜养殖区域的面积却更大 (图5)。
3. 讨论
本研究根据物种分布位点以及环境因子数据构建了SDM，对黄鳍金枪鱼和鲣的养殖潜在适宜区域进行了预测。结果表明，表层溶解氧、表层水温、底层水温、底层盐度是影响南海黄鳍金枪鱼和鲣养殖适宜区分布的重要环境变量。黄鳍金枪鱼潜在高适生区为南沙中西部渔场的东北部海域，鲣潜在高适生区为东沙渔场。自然环境因素作为养殖区域划分的基本要素，通过海洋中不同尺度的自然环境异质性来预测得到黄鳍金枪鱼和鲣的养殖潜在适宜区是可行的。
3.1 SDM在预测海水养殖潜力的适用性
当前和未来的海洋环境条件可能会在决定扩大海水养殖区域和物种方面发挥重要作用，温度、盐度、溶解氧作为影响金枪鱼类生长的要素，也是海水养殖能否成功的主要决定因素。研究物种适宜养殖区域是评估海水养殖可能性的第一步，SDM作为预测物种潜在海水养殖区域及分析环境适宜性结构不确定性的可靠工具[26]，已有研究表明，SDM适用于确定适宜养殖黄海鲑科的潜在海水养殖地[27-28]。本研究中，黄鳍金枪鱼和鲣的预测模型结果AUC值分别为0.848和0.872，均大于0.8，说明SDM对选择黄鳍金枪鱼和鲣的适宜养殖区域具有良好的预测能力。
在采用SDM选择金枪鱼类养殖区域时需要考虑以下因素：1) 离岸海水养殖物种通常被限制在一个有限区域内，即网箱养殖相比于近岸养殖而言，其技术与经济成本更高，考虑的因素更多。2) 我国南海不是金枪鱼的产卵海域，需要捕捞野生幼鱼和人工培育苗种进行养殖，其成活率有待考量；此外，养殖海域的环境因素随时间、空间动态变化 (如对养殖成功至关重要的温度、盐度等关键环境变量)，无法控制。3) 由于大多数养殖活动在相对较小的局部区域进行，设计海水养殖的SDM所需分辨率更高。
3.2 南海金枪鱼类潜在养殖区域的可行性
对黄鳍金枪鱼和鲣潜在养殖区域的分析表明，大部分海域在物理上适合其养殖、生长，并有可能扩大为新兴的离岸海水养殖业。本研究以潜在养殖SI为指标，定量评估了黄鳍金枪鱼和鲣的潜在养殖区，结果显示南海黄鳍金枪鱼和鲣的潜在高适生区较为集中，主要位于东沙渔场。目前一些水产养殖区域可能不是最优 (即小于高SI的区域)，但随着这些地区温度的上升，最适宜养殖区域面积将会扩大。而高SI区域面积因在不同海域中变化较大，所以需要密切跟踪其变化以确保准确性。此外，深远海养殖设施的发展、养殖技术的突破和国际市场的需求都是南海金枪鱼类养殖的驱动力。
3.3 限制因素
1) 气候变化。水产养殖对海洋环境和气候变化十分敏感，气候变化对全球渔业和水产养殖产生了直接或间接影响。面对极端气候变化，水产养殖遭遇了温度、pH、海平面上升、台风和海洋酸化的风险[29]。而要在南海东沙渔场海域实现黄鳍金枪鱼和鲣的养殖，台风是养殖过程中必须考虑的主要气候因素。利用可预测未来气候的区域情景方法，SDM可进一步评估未来气候因素的变化对海水养殖的潜在风险。
2) 模型预测能力。目前，生态位因子分析 (Ecological niche factor analysis, ENFA)、非参数概率生态位 (Non parametric probabilistic ecological niche, NPPEN)、随机森林 (Random forest, RF) 已被广泛用于估计水产养殖物种的潜在生产力。相比于采用SDM多集合方法来研究不同预测的变化差异，本文仅使用了SDMs集合中的最大熵模型，研究结果不能评估到所有潜在的适宜养殖区域。
3) 可持续承载能力。水产养殖区域的选取不仅要考虑环境因素的适宜性，还应在更广泛的养殖区域生态系统内考虑可持续承载能力。水产养殖可持续承载能力包括4个部分：物理承载能力，即物理上足以进行某种类型水产养殖的区域；生产承载能力，即从生物量或经济角度评估一个地区水产养殖生产能力的效用；生态承载能力，即在不产生难以接受的生态影响条件下，可能达到的最大生产水平；社会承载能力，即包括上述3种承载能力以及所有利益相关者之间的权衡，以满足人口(社会经济因素，如区域就业和传统渔业等) 和环境的需求[30]。本研究的重点是自然物理环境对物种水产养殖区域的影响，而在现实的水产养殖活动规划中往往需要运用多指标评价方法，充分考虑生产承载能力、生态承载能力和社会承载能力。
4) 养殖管理。水产养殖是一个复杂且长期的过程，养殖管理的便捷性与经济性成为降低生产成本的主要因素之一。随着大数据、物联网、人工智能等现代信息技术的广泛应用，水产养殖的精准化、精细化和智能化全过程生产管理更加普及，养殖管理水平和生产效率全面提升。此外，利用北斗导航系统可以预测养殖区域未来海况、潮流、风向等因子，特别是预防海域台风、赤潮等灾害造成的影响，以减少经济损失。但对生产者而言，前期的设施装备、智能管理系统以及技术研发等投入巨大，需要从技术、市场、资金等角度考虑金枪鱼养殖业发展的可行性。
5) 养殖设施。养殖装备设施作为保障深远海养殖可持续生产的基础要素，必须考虑金枪鱼类在养殖过程中的安全性。目前，根据养殖模式的划分，国内深远海渔业养殖设施装备主要包括自航养殖工船 (船舱养殖) 和网箱型养殖设施 (网箱养殖) 两种类型。自2017年以来，我国深远海养殖装备陆续投入使用，如山东青岛的“深蓝1号”、海南陵水的“5G+海洋牧场”及浙江舟山的“嵊海1号”等。但考虑到南海海域台风、风暴潮等影响和行业标准改进等因素，以及中国深远海养殖设施安全性与国际先进水平还存在一定差距，因此需要创新养殖装备核心技术，将海洋船舶工程装备与工业化养殖相融合，进一步加大对以大型养殖平台为核心的多元化养殖设备系统的研发力度，保障养殖全过程的安全性。
6) 养殖技术。养殖技术是物种水产养殖面积可持续扩大的基础条件。目前，我国已实现在外海养殖区开展黄鳍金枪鱼网箱养殖，建立了黄鳍金枪鱼深水网箱养殖技术体系，并开展了金枪鱼网箱养殖示范区建设[31]。但相较于其他国家，我国黄鳍金枪鱼深水网箱养殖技术仍处于起步阶段，扩大养殖区域需要考虑养殖技术的配套。而对于鲣来说，深水网箱养殖技术还需要更为深入的探索。
3.4 金枪鱼养殖建议
1) 采用SDM得出的金枪鱼潜在分布区域图，可为金枪鱼养殖区域选择提供指导。根据本研究结果，网箱可投放于离海南距离较近的东沙渔场，以减少运输成本。
2) 对养殖区域实施长期环境监测。环境变化对海水养殖业的影响巨大[32]，气候变化增加了水产养殖系统的复杂性和不确定性，如温度[33]、极端气候、溶解氧[34]及洋流[35]等，评估和预测环境因素对水产养殖的影响，并对在面临环境压力时，作出避免或减少损失的决策至关重要。
3) 评估和控制水产养殖对环境的影响。海水养殖对周围海域环境的影响复杂多样[36]。鱼类排出的粪便、饲料中含有的化学物质，甚至死鱼都可能污染环境[37]。例如，中国威海沿海地区长期的海水养殖活动导致pH呈现潜在下降趋势，整个研究区域的水质存在严重污染[38]。因此，研究应将物种养殖与周围环境的潜在影响结合起来。
4) 养殖区域要考虑海域使用权的立体多功能性及有限排他性。海域资源包括矿产资源、生物资源、海水资源等，同一海域空间的资源常以多种组合方式分布，使得海域具有立体开发的潜力。不同用海活动对海域自然属性的影响程度、影响方式和影响区域不同，不同用海活动间存在着兼容性和排他性，而同一海域空间的海域使用权是可以共存的。
5) 南海地理位置的特殊性对深远海养殖有诸多限制。南海作为重要的国际航道之一，拥有超过37条国际交通航线，也是我国极其重要的战略运输生命线，深远海养殖区域的选取要考虑航道通航因素。在国际主要航道分布图中，南海主要位于太平洋航线的东亚—马六甲国际主航道上，途经南沙中、西部渔场以及南沙西北部渔场的东部海域的船舶较多，而这些海域同时也是金枪鱼类的高适宜养殖区域，因此，在选择养殖区域时应充分考虑航区范围。
6) 海水养殖区域的选择，除考虑自然环境适宜度外，还需结合社会经济因素、科技水平及创新条件。金枪鱼养殖需要大规模的养殖工程系统，并依赖于高技术、高成本的陆基设施系统的开发、建设和维护。目前，我国金枪鱼养殖技术取得了新突破，如黄鳍金枪鱼深水网箱养殖示范区水体达6 600 m3，驯化成功率达到75.83%[31]；但同时也要加强养殖技术创新，日本作为全面掌握金枪鱼养殖技术的主要国家，其创新、技术、技能三位一体的基本经验[39]值得借鉴与学习，以促进我国金枪鱼养殖业的进一步发展。
4. 展望
根据南海黄鳍金枪鱼和鲣重要环境因子的最适宜范围，在今后研究中需要进一步分析其潜在高适生区重要环境因子的年际变动规律，从而更科学地分析高适生区的变化状况。SDM作为预测物种养殖适宜栖息地的有效工具，未来的水产养殖除了需要充分考虑中尺度环境异质性与气候变化之外，还需结合社会、经济、科技条件等因素，在理论生态适宜栖息地的基础上，得出南海金枪鱼潜在最佳适生区。在影响因素中不仅要全面考虑物理环境因素，还需要全面研究“生态适宜性”(即可持续承载力角度)，以精准预测适合黄鳍金枪鱼和鲣养殖的区域。此外，由于污染、海域内的人类活动和气候变化等因素，可能会导致目前适宜养殖的区域在未来变得不适宜。因此，未来的研究应考虑包括人类养殖活动对养殖区域环境的影响。
-
![]()
图 1 框架结构测量示意图
1. 头部腹侧末端；2. 吻端；3. 腹鳍起点；4. 头部背侧末端；5. 臀鳍起点；6. 背鳍起点
Figure 1. Schematic diagram of truss network measurements of H. molitrix
1. Ventral terminal of head; 2. Tip of snout; 3. Origin of pelvic fin; 4. Terminal of head back; 5. Origin of anal fin; 6. Origin of dorsal fin
表 1 表型数据测量表
Table 1 Phenotypic data measurements of H. molitrix
性状 Trait 性状 Trait x1 吻端至鳃盖腹侧后缘 (1—2) x9 头长 (HL) x2 鳃盖腹侧后缘至腹鳍起点 (1—3) x10 眼径 (ED) x3 吻端至鳃盖背侧后缘 (2—4) x11 全长 (TL) x4 鳃盖后上缘至腹鳍起点 (3—4) x12 体长 (BL) x5 腹鳍起点至臀鳍起点 (3—5) x13 体高 (BH) x6 鳃盖背侧后缘至背鳍起点 (4—6) x14 躯干长 (TL) x7 背鳍起点至腹鳍起点 (5—6) x15 尾柄长 (CPL) x8 鳃盖背侧后缘至鳃盖腹侧后缘 (1—4) x16 尾柄高 (CPH) 
下载: 导出CSV
表 2 不同月龄长丰鲢体质量与形态性状的描述性统计
Table 2 Descriptive statistics for body mass and morphometric traits of H. molitrix at different ages
性状
Trait6月龄 6 Months 12月龄 12 Months 36月龄 36 Months 均值
Mean标准差
SD变异系数
CV/%均值
Mean标准差
SD变异系数
CV/%均值
Mean标准差
SD变异系数
CV/%y/g 19.95 5.34 26.78 48.46 13.88 28.63 1 040.65 75.20 7.23 x1/cm 2.20 0.25 11.44 3.49 0.34 9.80 9.20 0.60 6.55 x2/cm 3.01 0.41 13.75 3.78 0.37 9.81 11.46 0.64 5.62 x3/cm 2.62 0.31 11.68 3.83 0.40 10.34 10.63 0.58 5.42 x4/cm 3.70 0.45 12.06 4.99 0.49 9.87 14.26 0.55 3.87 x5/cm 2.48 0.30 12.23 3.15 0.38 12.18 9.35 0.47 5.01 x6/cm 3.01 0.31 10.23 4.08 0.42 10.41 12.26 0.58 4.72 x7/cm 3.17 0.41 12.80 4.25 0.50 11.67 12.29 0.46 3.73 x8/cm 2.52 0.25 10.03 3.41 0.33 9.63 9.77 0.33 3.40 x9/cm 2.97 0.33 11.22 4.39 0.43 9.71 12.57 0.59 4.69 x10/cm 0.70 0.08 11.10 0.90 0.07 7.68 1.34 0.09 6.82 x11/cm 12.23 1.07 8.73 16.98 1.61 9.47 47.26 1.64 3.46 x12/cm 9.79 0.84 8.60 14.04 1.41 10.02 39.74 1.75 4.40 x13/cm 3.26 0.34 10.51 4.33 0.47 10.94 11.40 0.51 4.48 x14/cm 4.10 0.44 10.66 5.38 0.58 10.71 15.98 0.68 4.25 x15/cm 1.65 0.17 10.09 2.34 0.32 13.58 6.56 0.45 6.82 x16/cm 1.05 0.12 10.98 1.46 0.17 11.64 4.05 0.18 4.33
下载: 导出CSV
表 3 不同月龄长丰鲢形态性状的特征值和累计贡献率
Table 3 Eigenvalue and cumulative contribution rate of H. molitrix morphological traits at different ages
主成分
Principal component特征值 Eigenvalue 累计贡献率 Cumulative contribution rate/% 6月龄
6 Months12月龄
12 Months36月龄
36 Months6月龄
6 Months12月龄
12 Months36月龄
36 Months1 10.805 14.063 6.498 67.531 87.894 40.612 2 1.517 0.541 2.065 77.012 91.276 53.518 3 0.844 0.417 1.374 82.286 93.880 62.103 4 0.652 0.242 1.338 86.363 95.391 70.465
下载: 导出CSV
表 4 不同月龄入选的主成分特征向量
Table 4 Principal component eigenvector of H. molitrix at different ages
月龄 Month age 6 12 36 主成分 Principal component 1 2 3 4 1 2 3 4 1 2 3 4 x1 0.727 −0.290 −0.021 0.528 0.868 0.122 −0.431 −0.095 0.552 −0.486 0.209 0.420 x2 0.790 −0.247 0.188 −0.360 0.917 −0.117 0.236 −0.106 0.502 0.288 −0.160 −0.647 x3 0.831 −0.319 0.016 0.227 0.936 0.086 −0.184 −0.126 0.570 −0.377 0.354 −0.194 x4 0.903 −0.208 0.113 −0.177 0.972 −0.077 0.088 −0.047 0.860 −0.035 −0.036 −0.148 x5 0.828 −0.205 0.159 0.031 0.939 −0.154 −0.008 0.088 0.649 0.196 −0.328 0.090 x6 0.847 −0.243 0.269 0.087 0.944 −0.081 0.086 −0.007 0.790 −0.153 −0.200 0.049 x7 0.867 −0.154 0.122 −0.277 0.979 −0.080 0.080 −0.025 0.747 0.266 −0.089 −0.272 x8 0.910 0.010 −0.313 −0.048 0.970 0.029 −0.001 −0.013 0.808 −0.047 0.245 −0.086 x9 0.853 −0.374 0.074 0.024 0.956 0.075 −0.147 −0.156 0.626 −0.508 0.399 −0.137 x10 0.703 −0.098 −0.622 0.033 0.716 0.656 0.225 0.048 0.536 −0.057 0.186 0.398 x11 0.865 0.420 −0.002 0.020 0.983 −0.016 0.007 0.035 0.853 −0.044 −0.091 0.160 x12 0.865 0.431 0.027 0.052 0.994 −0.028 −0.009 −0.015 0.273 0.718 0.433 0.298 x13 0.839 0.360 −0.033 −0.004 0.975 −0.058 0.094 −0.046 0.137 0.758 0.492 0.091 x14 0.873 0.056 −0.319 −0.144 0.963 −0.130 0.064 0.062 0.726 0.333 −0.332 0.125 x15 0.643 0.517 0.292 0.208 0.886 0.014 −0.170 0.400 0.354 0.030 −0.468 0.496 x16 0.754 0.436 0.055 −0.072 0.965 −0.064 0.075 0.031 0.685 0.051 −0.024 −0.173
下载: 导出CSV
表 5 不同月龄长丰鲢性状对体质量的直接与间接通径系数
Table 5 Direct and indirect path coefficients of traits to body mass of H. molitrix at different ages
月龄
Month age性状
Trait相关系数
Correlation coefficient直接作用
Direct effect间接作用 Indirect effect ∑ x1 x3 x7 x8 x11 x12 x13 x15 x16 6 x8 0.812 0.105 0.708 − − − − 0.253 0.262 0.236 −0.043 − x11 0.950 0.330 0.621 − − − 0.081 − 0.333 0.269 −0.062 − x12 0.952 0.344 0.608 − − − 0.080 0.320 − 0.272 −0.064 − x13 0.920 0.312 0.609 − − − 0.079 0.284 0.301 − −0.055 − x15 0.668 −0.085 0.753 − − − 0.054 0.240 0.257 0.203 − − 12 x1 0.843 0.091 0.752 − − − 0.160 0.222 − 0.193 − 0.176 x8 0.966 0.191 0.774 0.076 − − − 0.257 − 0.230 − 0.211 x11 0.970 0.266 0.704 0.076 − − 0.184 − − 0.231 − 0.213 x13 0.964 0.241 0.722 0.073 − − 0.183 0.255 − − − 0.212 x16 0.962 0.222 0.739 0.072 − − 0.181 0.255 − 0.230 − − 36 x3 0.269 −0.146 0.416 − − 0.174 0.139 0.110 − −0.008 − − x7 0.708 0.462 0.246 − −0.055 − 0.133 0.124 − 0.044 − − x8 0.612 0.253 0.359 − −0.080 0.242 − 0.171 − 0.026 − − x11 0.592 0.252 0.340 − −0.064 0.227 0.172 − − 0.005 − − x13 0.342 0.174 0.168 − 0.007 0.117 0.038 0.007 − − − −
下载: 导出CSV
表 6 不同月龄长丰鲢形态性状对体质量的决定系数
Table 6 Determinant coefficient of morphometric attributes on body mass of H. molitrix at different ages
月龄
Month age性状
Traitx1 x3 x7 x8 x11 x12 x13 x15 x16 6 x8 − − − 0.011 − − − − − x11 − − − 0.053 0.109 − − − − x12 − − − 0.055 0.220 0.118 − − − x13 − − − 0.049 0.178 0.188 0.097 − − x15 − − − −0.009 −0.041 −0.044 −0.034 0.007 − 12 x1 0.008 − − − − − − − − x8 0.029 − − 0.036 − − − − − x11 0.040 − − 0.098 0.071 − − − − x13 0.035 − − 0.088 0.123 − 0.058 − − x16 0.032 − − 0.081 0.113 − 0.102 − 0.049 36 x3 − 0.021 − − − − − − − x7 − −0.051 0.213 − − − − − − x8 − −0.041 0.122 0.064 − − − − − x11 − −0.032 0.115 0.087 0.064 − − − − x13 − 0.002 0.041 0.013 0.002 − 0.030 − −
下载: 导出CSV
-
[1] 梁宏伟, 李忠, 罗相忠, 等. 长丰鲢与长江鲢形态差异与判别分析[J]. 水生生物学报, 2015, 39(5): 1059-1064. [2] 杨希, 王明涛, 白海锋, 等. 长丰鲢水花苗种在渭北地区培育效果对比分析[J]. 河北渔业, 2017(10): 32-33, 44. [3] YAKUBU A, MUHAMMED M M, ARI M M, et al. Correlation and path coefficient analysis of body weight and morphometric traits of two exotic genetic groups of ducks in Nigeria[J]. Bang J Anim Sci, 2015, 44(1): 1-9. doi: 10.3329/bjas.v44i1.23112
[4] 吴水清, 郑乐云, 罗辉玉, 等. 杂交石斑鱼 (斜带石斑鱼♀×赤点石斑鱼♂) 与其亲本形态性状比较研究[J]. 南方水产科学, 2017, 13(5): 47-54. [5] 吴水清, 罗辉玉, 张哲, 等. 不同月龄云龙石斑鱼表型性状的主成分与通径分析[J]. 大连海洋大学学报, 2019, 34(5): 680-687. [6] 李俊伟, 区又君, 温久福, 等. 室内循环水和池塘养殖四指马鲅的生长性能及形态性状与体质量的相关性研究[J]. 南方水产科学, 2020, 16(1): 27-35. [7] 黄建盛, 陈刚, 张健东, 等. 褐点石斑鱼不同月龄形态性状的主成分及通径分析[J]. 水产学报, 2017, 41(7): 1105-1115. [8] LAGO A A, REIS-NETO R V, REZENDE T T, et al. Quantitative analysis of black blotching in a crossbred red tilapia and its effects on performance traits via a path analysis methodology[J]. J Appl Genet, 2019, 60(3/4): 393-400.
[9] SHI D J. Correlation between main morphological traits and body weight of mandarin fish difference analysis of male and female morphological traits[J]. Fish Sci Aquacult, 2019, 1(1): 1-6.
[10] 方伟, 周胜杰, 赵旺, 等. 黄鳍金枪鱼5月龄幼鱼形态性状对体质量的相关性及通径分析[J]. 南方水产科学, 2021, 17(1): 52-58. [11] 黄建盛, 郭志雄, 陈刚, 等. 4月龄军曹鱼幼鱼形态性状与体质量的相关性及通径分析[J]. 海洋科学, 2019, 43(8): 72-79. [12] 张优平, 白海锋, 李引娣, 等. 微生物制剂对长丰鲢苗种培育效果的影响研究[J]. 水产养殖, 2020, 41(8): 39-40. [13] 杨希, 李涛, 白海锋, 等. 不同投喂模式培育长丰鲢的对比试验[J]. 水产养殖, 2017, 38(8): 13-15. [14] 韦玲静, 叶香尘, 梁克, 等. 广西本地鲢鱼和长丰鲢肌肉营养成分分析与品质评价[J]. 食品工业科技, 2020, 41(3): 285-290. [15] 全书月, 李忠, 梁宏伟, 等. 长丰鲢卵巢发育后期的肝脏、卵巢营养成分及脂肪酸变化[J]. 淡水渔业, 2018, 48(5): 79-85. [16] 郭华阳, 陈明强, 王雨, 等. 黄边糙鸟蛤野生群体主要经济性状间的相关性及通径分析[J]. 南方水产科学, 2013, 9(2): 1-8. [17] 叶乐, 王海山, 杨超杰, 等. 饥饿对眼斑双锯鱼 (Amphiprion ocellaris) 幼鱼形态性状与体质量相关性的影响[J]. 生态学报, 2021, 41(2): 1-8. [18] 苟盼盼, 王秀利, 窦冬雨, 等. 红鳍东方鲀不同家系群体的形态性状差异与相关性分析[J]. 大连海洋大学学报, 2019, 34(5): 674-679. [19] 王雪, 李莉, 刘元文, 等. 3月龄黑鲷形态性状对体质量的影响效果分析[J]. 大连海洋大学学报, 2019, 34(2): 233-238. [20] 李莉, 王雪, 菅玉霞, 等. 不同月龄大泷六线鱼形态性状与体质量的相关性及通径分析[J]. 上海海洋大学学报, 2019, 28(1): 58-66. [21] 窦亚琪, 梁旭方, 易提林, 等. 翘嘴鳜不同月龄性状的主成分与判别分析[J]. 中国水产科学, 2014, 21(6): 1116-1124. [22] 刘贤德, 蔡明夷, 王志勇, 等. 不同生长时期大黄鱼形态性状与体重的相关性分析[J]. 热带海洋学报, 2010, 29(5): 159-163. [23] 耿绪云, 马维林, 李相普, 等. 梭鱼 (Liza haematocheila) 外部形态性状对体重影响效果分析[J]. 海洋与湖沼, 2011, 42(4): 530-537. [24] 王新安, 马爱军, 许可, 等. 大菱鲆幼鱼表型形态性状与体重之间的关系[J]. 动物学报, 2008, 54(3): 540-545. [25] 张赛赛, 罗珺, 陈博锦, 等. 不同群体刀鲚形态差异研究[J]. 淡水渔业, 2020, 50(3): 26-33. [26] 何铜, 刘小林, 杨长明, 等. 凡纳滨对虾各月龄性状的主成分与判别分析[J]. 生态学报, 2009, 29(4): 2134-2142. [27] 李培伦, 刘伟, 王继隆, 等. 马苏大麻哈鱼不同月龄表型性状的主成分与判别分析[J]. 水产科学, 2017, 36(6): 707-713. [28] 韩慧宗, 姜海滨, 王斐, 等. 许氏平鲉不同月龄选育群体形态性状的主成分与通径分析[J]. 水产学报, 2016, 40(8): 1163-1172. [29] 汪月书, 李彩娟, 许郑超, 等. 梭鲈不同月龄性状的主成分分析与判别分析[J]. 水产养殖, 2016, 37(4): 16-22. [30] ALAKE C O. Traits associated with pod yield, seed Fe, and Zn contents in Okra landraces: a path analysis[J]. Int J Veg Sci, 2019, 26(1): 1-18.
[31] BOTELHO H A, LAGO A A, COSTA A L, et al. Application of morphometric measures in estimation of body weight and discrimination of Astyanax lacustris and Astyanax fasciatus[J]. Aquacult Res, 2019, 50(9): 2429-2436. doi: 10.1111/are.14196
[32] IKPEME E V, EKERETTE E E, UDENSI O U, et al. Assessment of morphological variation in wild and cultured populations of tilapia fish (Oreochromis niloticus)[J]. J Adv Biol Biotechnol, 2017, 13(2): 1-10.
[33] 魏海军, 邓正华, 陈明强, 等. 棕带仙女蛤数量性状的相关与通径分析[J]. 南方水产科学, 2019, 15(6): 34-40. [34] 区又君, 吉磊, 李加儿, 等. 卵形鲳鲹不同月龄选育群体主要形态性状与体质量的相关性分析[J]. 水产学报, 2013, 37(7): 961-969. [35] 张雷雷, 滕爽爽, 李腾腾, 等. 不同月龄青蛤形态性状对活体质量的影响分析[J]. 海洋科学, 2019, 43(12): 74-80. [36] 张进, 梁旭方, 郭小泽, 等. 翘嘴鳜形态性状对体质量的影响[J]. 华中农业大学学报, 2013, 32(5): 100-105. [37] 陈健, 郭丹, 翟子钦, 等. 西施舌形态性状与体质量性状的相关性分析[J]. 南方水产科学, 2021, 17(1): 45-51. [38] 耿平, 张魁, 陈作志, 等. 北部湾蓝圆鲹生物学特征及开发状态的年际变化[J]. 南方水产科学, 2018, 14(6): 1-9. [39] 冒树泉, 王秉利, 许鹏. 饥饿不同时间后投喂对许氏平鲥生长及体成分的影响[J]. 水产学杂志, 2017, 30(1): 26-31. -
期刊类型引用(10)
1. 肖辉裕，马爱军，孙志宾，王新安，孙伟恒，慈小静，俞兰良. 不同月龄膨腹海马表型形态性状与体重之间的关系. 水产学报. 2025(02): 139-153 . 
百度学术
2. 赵春暖，王晓梅，于涛，单秀娟，任利群，郑言鑫，李波，李致远，王喜文，陈洋. 膨腹海马不同生长阶段的形态特征研究. 南方农业学报. 2024(04): 1216-1224 . 百度学术
3. 段国庆，周华兴，汪焕，凌俊，胡玉婷，潘庭双，杨敏，邬凌云，江河. 60日龄瓦氏黄颡鱼选育群体形态性状对体重的影响分析. 安徽农业大学学报. 2023(01): 78-85 . 百度学术
4. 杨浩，蒲艳，高雷，段辛斌，刘绍平，陈大庆，李云. 基于几何形态测量学的磷酸三(1, 3-二氯-2-丙基)酯对鲢仔鱼形态影响研究. 水生生物学报. 2023(05): 819-826 . 百度学术
5. 张凯强，李刚，王庆龙，张龙，陈志方，温海深，李吉方，徐扬涛. 淡水养殖花鲈形态性状与体质量的相关性及通径分析. 中国海洋大学学报(自然科学版). 2023(09): 58-68 . 百度学术
6. 张涛，周剑光，张林，甘金华，何力. 大口黑鲈形态特征、染色体核型及眼晶状体乳酸脱氢酶的电泳分析. 西北农业学报. 2023(10): 1503-1511 . 百度学术
7. 孙以康，李美帅，陈四清，张岩，常青，李凤辉，葛建龙，边力. 不同日龄中华蛸形态性状对体质量的影响分析及曲线拟合. 海洋渔业. 2022(05): 568-576 . 百度学术
8. 李晓晖，邹桂伟，梁宏伟，何萱，沙航，冯翠，崔峰. 持续低氧胁迫对长丰鲢血液指标和组织抗氧化酶活性及相关基因表达的影响. 湖南农业大学学报(自然科学版). 2022(06): 723-729 . 百度学术
9. 夏威威，刘海洋，欧密，罗青，黄苏静，赵训金，熊炳源，吴国宏，陈昆慈，赵建. 斑鳢后备亲本体质量与形态性状间的关系分析. 南方农业学报. 2022(12): 3584-3592 . 百度学术
10. 童立豪，谭凡民，罗砚，倪孔平. 红树林人工湿地生态养殖黄鳍鲷形态性状对体质量的影响. 广西科学院学报. 2021(03): 248-255 . 百度学术
其他类型引用(1)
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 11
粤公网安备 44010502001741号
