Biological study of reproductive characteristics of Acanthopagrus schlegeli in Daya Bay
-
摘要:
为开展大亚湾黑鲷(Acanthopagrus schlegeli)的增殖和保护，于2018年1月、3月、4月和8—12月共采集了721尾黑鲷样本，研究了大亚湾黑鲷繁殖生物学的周年变化特征。结果表明，1)繁殖期的雌性全长为175~487 mm，均值为 (297.92±56.61) mm，体质量为90.96~1 907.12 g，均值为 (511.36±307.53) g，雄性全长为110~440 mm，均值为 (237.61±61.05) mm，体质量为21.92~1 370.00 g，均值为 (277.76±210.33) g，雌性个体全长显著大于雄性个体全长(P<0.001)；2)性腺成熟系数和性腺发育表明大亚湾黑鲷繁殖期为11月至翌年3月；3)具有雌雄同体个体存在，表明该物种会发生性逆转，且表现为雄性个体先成熟再向雌性转变；4) 50%性成熟雌性全长达到248.52 mm，雄性为148.03 mm；5)卵径频率分布表明黑鲷为分批产卵型鱼类，绝对繁殖力为84 100~2 655 168粒，均值为 (728 732±499 547) 粒；全长相对繁殖力为380.54~8 271.55粒·mm−1，均值为 (2 731.44±1 630.93) 粒·mm−1，体质量相对繁殖力为277.58~3 592.83粒·g−1，均值为 (1 443.87±766.82) 粒·g−1，绝对繁殖力与全长、体质量及性腺质量均呈显著线性相关。
Abstract:We collected 721 individuals of Acanthopagrus schlegeli in Daya Bay in January, March, April, August−December, 2018, and then investigated the annual variation of reproductive characteristics of A. schlegeli in Daya Bay so as to provide a theoretical basis for its artificial breeding and conservation. The results show that: 1) the spawning period falls between November and March. The total length ranged 175−487 mm for females [averge (297.92±56.61) mm]. The body mass ranged 90.96−1 907.12 g for females [average (511.36±307.53) g]. The total length ranged 110−440 mm for males [average (237.61±61.05) mm]. The body mass ranged 21.92−1 370.00 g for males [average (277.76±210.33) g]. The full lengths of females were significantly longer than those of males (P<0.001). 2) The gonad maturity coefficient and gonad development indicate that the reproductive period of A. schlegeli in Daya Bay was from November to March. 3) The existence of androgynous individuals indicates that the species will undergo sexual reversal, and the males will mature first and then change to females. 4) The total lengths (TL) at 50% sex maturity for females and males were 248.52 mm and 148.03 mm, respectively. 5) A. schlegeli is a multiple spawner with indeterminate fecundity by oocyte size-frequency distributions. The number of the individual fertility ranged from 84 100 to 2 655 168 granules (average 728 732±499 547). The relative fertility of total length was 380.54−8 271.55 eggs·mm−1 [(2 731.44±1 630.93) egg·mm−1]. The relative fertility of body mass was 277.58−3 592.83 eggs·g−1 [average (1 443.87±766.82) eggs·g−1]. There was a linear correlation among absolute fecundity and total length, body mass and gonadal mass.
-
Keywords:
- Acanthopagrus schlegeli /
- sex reversal /
- gonadosomatic index /
- spawns in batches /
- fertility
-
增殖放流是目前有效恢复渔业资源的方法之一[1]，掌握放流种类的生态学特征是实施有效增殖放流的基本要求，选择适宜放流海区是确保放流物种成活率的关键。其中，物种的生境适宜性分析是保障物种增殖放流获得最大效益的重要环节。该过程需分析目标物种与所处环境之间的关系，分析目标物种生境中关键因子的需求，构建目标物种栖息地选择模型，最后将其潜在适生地理分布区直观地表达在数字化地形图或遥感图上[2]。头足类具有市场价值高、生长周期短和生长速度快等优点[3]，对头足类增殖放流时进行物种的生境适宜性分析有利于进一步实现渔业资源可持续发展和蓝色经济增长的目标。
最大熵 (Maximum Entropy, MaxEnt) 模型是基于物种已知的分布数据和相关环境因子对物种分布及其生态需求进行模拟，从而分析得出物种的潜在适生区的物种分布模型 (Species Distribution Models, SDMs)[4]。在头足类研究领域，目前物种分布模型的应用以广义加性模型(Generalized Additive Model, GAM) 为主[5-7]。但与GAM模型相比，MaxEnt模型的预测结果较为准确[8-9]。MaxEnt模型的预测精度在一定程度上受所选生态因子的影响，但通过筛选生态因子和参数优化后，模型可在研究样本量较小的情况下获得较高的预测精度[10-11]。该模型不仅能有效处理环境因子间的相互作用，还可定量揭示环境因子对物种分布的影响[12]。适用于不同的生境类型、空间尺度和分类类群[13]，在海洋生物保护方面极具发展前景[14]，常被用于分析物种栖息地需求、气候变化对物种分布的影响、物种入侵监测和自然保护区的选择和保护[15]。例如，MaxEnt模型被应用于预测南太平洋长鳍金枪鱼 (Thunnus alalunga) 的栖息地[16]，模拟水温变化对太平洋柔鱼 (Ommastrephes bartramii) 潜在栖息地的影响[17]，分析底层鱼类在阿拉斯加湾的生境适宜性等[18]。
曼氏无针乌贼 (Sepiella maindroni) 是我国重要的头足类生物之一，隶属于软体动物门、头足纲、二鳃亚纲、十腕目、乌贼科、无针乌贼属。其可食部分比重大 (约占总体质量的92%)、蛋白质含量高且肉质富含钙 (Ca)、磷 (P)、铁 (Fe) 等元素[19]，具有较高的营养和经济价值。曼氏无针乌贼曾是东海的优势物种，在浙江省的年产量曾高达6万吨[20]，但受过度捕捞、亲体过度开发以及产卵场破坏等影响，其产量自1981年后持续下降[21]，直至20世纪90年代初基本绝迹[22]。开展大规模增殖放流后，其产量在浙江海域呈上升趋势，至2008年增殖效果初显[21,23]，但目前其产量仍不能满足人们的需求[24]。舟山海域岛屿众多，适宜曼氏无针乌贼卵附着的藻类丰富[25]，是其重要的栖息地和产卵地。该海域内曼氏无针乌贼种群胴长范围为29~141 mm，体质量范围为5.1~341.8 g[26]，目前对该海域内曼氏无针乌贼生境适宜性的研究较少。本文以舟山海域曼氏无针乌贼为研究对象，收集其分布数据和舟山海域的海洋环境数据，结合ArcGIS 10.8软件和MaxEnt模型对曼氏无针乌贼进行生境适宜性评价，以期为曼氏无针乌贼的增殖放流选址提供参考，并为其资源养护提供科学支撑。
1. 材料与方法
1.1 物种分布数据
构建MaxEnt模型需要物种分布点的经纬度数据。本研究的物种分布数据来源于2017—2021年浙江近海的底拖网调查数据、2021—2022年自主航次的底拖网调查数据和全球生物多样性信息机构GBIF (Global Biodiversity Information Facility, https://www.gbif.org/)。2017—2021年调查范围为120°E—131°E、23°N—33°N，站位采用均匀网格设置，每隔经纬度15'设置1个采样站点。2021—2022年自主航次调查范围为122°30'E—123°20'E、29°36'N—30°30'N。调查船为“浙普渔68712”单船底拖网渔船，使用单拖网，网目尺寸为80 mm，网囊最小网目尺寸为20 mm，每网拖曳1 h，平均拖速3 kn。调查时记录捕获曼氏无针乌贼时的经纬度信息。GBIF是一个由世界各国政府资助的国际开放基础数据库，数据来源包括物种名录、书籍、相关论文、动植物考察和观察资料以及数字化的在线数据集，目前已有5 000多万条分布点的数据记录。为扩大样本量，利用该数据库及相关文献资料，在底拖网调查数据的基础上补充曼氏无针乌贼的分布点数据。为避免因群集效应造成的取样偏差，利用ENMTools (http://purl.oclc.org/enmtools) 剔除冗余数据[27]。最终共得到曼氏无针乌贼有效分布样点167个 (图1)，通过Excel 2020软件转换为csv格式。采用ArcGIS 10.8 软件绘制地图，数据来源于全国地理信息资源目录服务系统 [https://www.webmap.cn/，1∶100万公众版基础地理信息数据 (2021)]。
1.2 海洋环境数据
海洋环境数据集来源于Bio-ORACLE (https://bio-oracle.org) 和全球海洋环境数据库 (Global Marine Environment Datasets, GMED)。研究表明，海表面温度 (Sea surface temperature, SST)、海表面盐度 (Sea surface salinity, SSS)、pH、光照均可影响曼氏无针乌贼的生长和活动[28-30]。在此基础上补充选取影响海洋生物分布的相关生态因子 (海洋流速、初级生产力、硝酸盐含量等) 共同构成模型预实验的生态因子集合。
相关性过高的生态因子会影响模型的运行结果[31]，为提高预测结果的准确性，对预实验中的生态因子进行筛选。利用SPSS 20.0软件中的相关性分析法和MaxEnt模型中的刀切法 (Jackknife)，以贡献率小且相关性绝对值 |r|>0.8为标准，剔除易导致模型过拟合的生态因子[32]，最终获得与曼氏无针乌贼相关的关键生态因子 (表1)。
表 1 筛选出的关键生态因子Table 1. Selected key ecological factors生态因子
Ecological factor单位
Unit海流流速最大值 Maximum currents velocity m·s−1 海流流速最小值 Minimum currents velocity m·s−1 溶解氧质量浓度最小值 Minimum dissolved oxygen mg·L−1 底层光照强度最大值 Maximum light at bottom μmol·(m2·s)−1 硝酸盐质量浓度最大值 Maximum nitrate mg·L−1 磷酸盐质量浓度最大值 Maximum phytoplankton mg·L−1 初级生产力最小值 Minimum primary productivity g·(m3·day)−1 海表盐度最小值 Minimum sea surface salinity ‰ 海表温度最大值 Maximum sea surface temperature ℃ 海表温度最小值 Minimum sea surface temperature ℃ 酸碱度 pH 1.3 模型参数优化及运行
导入数据至MaxEnt模型，运用由Cobos等[33]发明的kuenm软件包进行MaxEnt模型的参数优化。其利用了R语言和MaxEnt模型的功能，帮助实现校准、创建和评估最终的模型。为最大化获得参数优化效果，设置调控倍频 (Regularization multiplier, RM) 区间为0.5~4.0，间距为0.5 (0.5、1.0、1.5、2.0、2.5、3.0、3.5、4.0)，其他设置为系统默认值。优先选择AIC (Akaike information criterion) 值最小和遗漏率最低的调控系数。AIC值是衡量统计模型拟合优良性的一种标准，能权衡所估计模型拟合数据的优良性，遗漏率是用于所用的数据在特定阈值上的评估，可分析建模时的数据是否存在空间自相关性[34]。根据运行结果，设置调控系数为3.0。
导入处理后的数据至MaxEnt 3.4.4 (https://biodiversityinformatics.amnh.org/open_source/maxent/)运行，勾选MaxEnt模型中Jackknife工具，获得主要生态因子对曼氏无针乌贼在舟山海域分布影响的权重值。随机选取75%的分布点用于建立模型，其余25%的分布点用于模型验证，且为了保证模型稳定性，采用10次重复交叉验证[35]。模型其他参数设置为默认值，结果以Logistic格式的栅格数据输出。
1.4 划分物种适生区等级
参考IPCC报告有关评估可能性的划分方法[36]，根据曼氏无针乌贼分布的概率值，将其适宜区分为5个等级，分别为P<0.2、0.2≤P<0.4、0.4≤P<0.6、0.6≤P<0.8和0.8≤P<1，依次为非适生区 (Non-adaptive areas)、低适生区 (Low-adaptive areas)、中适生区 (Middle-adaptive areas)、次高适生区 (Sub-superior suitable area) 和高适生区 (High-adaptive areas)。
2. 结果
2.1 模型结果精度及检验
本文使用受试者工作特征曲线 (Receiver operating characteristic curve, ROC) 验证最大熵模型的预测精度。ROC曲线下面积 (Area under curve, AUC) 是评估预测物种分布的模型准确性的主要工具[37]。AUC值越高、越接近1则表示模型的性能越好，当AUC值大于0.9时，表示结果准确度较高[38]。应用表1中的生态因子进行建模，AUC值为0.995 (图2)。证明选择表1中的生态因子建立MaxEnt模型进行预测的结果准确度较高，可用于曼氏无针乌贼在舟山海域的生境适宜性分析。
2.2 物种对生态因子的响应
表2反映了生态因子对MaxEnt模型预测的相对贡献率，其中SST最大值、SSS最小值和pH是影响曼氏无针乌贼在舟山海域分布的主要生态因子。三者对模型的相对贡献率分别为27.8%、25.5%和20.3%，均超过20%，并且其相对贡献率之和占所有生态因子相对贡献率之和的一半以上。初级生产力最小值、磷酸盐质量浓度最大值、SST最小值和硝酸盐质量浓度最大值对模型的相对贡献率介于1%~15%，对曼氏无针乌贼分布的影响小于SST最大值、SSS最小值和pH。溶解氧质量浓度最小值、底层光照强度最大值、最大海流流速和最小海流流速的相对贡献率较小，均小于1%。
表 2 生态因子对模型的相对贡献率Table 2. Percent contributions of ecological factors to MaxEnt model生态因子
Ecological factor相对贡献率
Relative
contribution
rate/%海表温度最大值 Maximum sea surface temperature 27.8 海表盐度最小值 Minimum sea surface salinity 25.5 酸碱度 pH 20.3 初级生产力最小值 Minimum primary productivity 12.2 磷酸盐质量浓度最大值 Maximum phytoplankton 8.7 海表温度最小值 Minimum sea surface temperature 2.5 硝酸盐质量浓度最大值 Maximum nitrate 2.2 溶解氧质量浓度最小值 Minimum dissolved oxygen 0.4 底层光照强度最大值 Maximum light at bottom 0.2 海流流速最大值 Maximum currents velocity 0.2 海流流速最小值 Minimum currents velocity 0.1 生态因子重要性结果显示 (图3)，仅包含SST最大值、SSS最小值和pH等生态因子建模时，模型获得的增益值较大。排除相对贡献率较大的生态因子 (如：SST最大值、SSS最小值和pH) 后，分别利用除这类生态因子外的其他生态因子建模，此时模型所获得的增益值较小。表明模型预测过程中SST最大值、SSS最小值和pH包含有价值的信息，无法被其他生态因子替代。
![]()
图 3 刀切法检验生态因子对曼氏无针乌贼分布影响的重要程度Figure 3. Importance of ecological factors on distribution of S. maindroni by Jackknife test生态因子响应曲线可反映生态因子与曼氏无针乌贼存在概率之间的关系 (图4)。因相对贡献率较高的生态因子对曼氏无针乌贼分布的影响较大，故重点分析曼氏无针乌贼对SST最大值、SSS最小值、pH等生态因子的响应。以0.5概率值作为分界点，结果显示，舟山海域内SST最大值介于25.5~28.5 ℃，SSS最小值介于20‰~30‰，pH介于7.8~8.3时，曼氏无针乌贼存在的概率较大。
在SST最大值介于25.5~26.6 ℃时，曼氏无针乌贼的存在概率随温度升高而增加，温度达到26.6 ℃后其存在概率随温度升高而减少。SSS最小值响应曲线呈先增加后减少的变化，SSS最小值达25‰时曼氏无针乌贼物种存在的概率最大。在pH的响应曲线中，pH为8.2时曼氏无针乌贼的存在概率达到最大值。pH介于7.1~8.2时，其存在概率随pH升高呈上升趋势；pH介于8.2~8.4时，存在概率随pH升高呈下降趋势。
2.3 物种的适生区分布及检验
利用ArcGIS 10.8软件对模型结果进行重分类，获得曼氏无针乌贼在舟山海域的适生区分布图 (图5)。经计算，曼氏无针乌贼在舟山海域的适生区面积为30 213.40 km²，占研究区域总面积的75.11%。适生区分布范围较广，呈南北方向分布，北至马鞍列岛海域，南至舟山本岛东侧海域。高适生区面积为4 600.78 km²，主要分布在嵊泗列岛、马鞍列岛和中街山列岛海域，少量分散在马鞍列岛东侧。该海域内曼氏无针乌贼次高适生区面积为8 205.52 km²，中适生区面积为9 723.31 km²，低适生区面积为7 683.78 km²，三者面积相近 (表3)。高适生区、次高适生区分布范围内包含了舟山市内的3个海洋特别保护区 (图6)。
![]()
图 5 曼氏无针乌贼在舟山海域的适生区分布Figure 5. Distribution of suitable habitat of S. maindroni in Zhoushan Islands![]()
图 6 舟山海域内海洋特别保护区地理位置图Figure 6. Geographical location map of marine special reserves in Zhoushan sea area表 3 适生区面积及其占海域面积百分比Table 3. Area of suitable habitat and its percentage of sea area区域
Area概率
Probability面积
Square/km2占比
Percentage/%高适生区 High-adaptive area 0.8≤P<1 4 600.78 11.44 次高适生区 Sub-superior suitable area 0.6≤P<0.8 8 205.52 20.40 中适生区 Middle-adaptive area 0.4≤P<0.6 9 723.31 24.17 低适生区 Low-adaptive area 0.2≤P<0.4 7 683.78 19.10 非适生区 Non-adaptive area P<0.2 10 007.89 24.88 研究区域 Research area 40 221.29 100.00 为检验模型结果，采取访问调查的方式走访舟山市各区县典型渔村，开展曼氏无针乌贼捕捞出现海域调研。范围包括定海区 (金塘镇、长白乡)、普陀区 (沈家门、朱家尖、六横岛、桃花镇、东极岛)、岱山县 (衢山镇、长涂镇) 和嵊泗县 (嵊山镇、枸杞乡、洋山镇)。据2021年舟山市统计年鉴，2020年定海区、普陀区、岱山县和嵊泗县渔业人口分别为12 575、73 596、56 961和44 007人。综合考虑人力物力等客观因素及调研结果准确度，2022年共访问调查了1 285人，其中定海区87人、普陀区505人、岱山县391人和嵊泗县302人。所调查渔民中70.3%为专业从业人员，19.2%为兼业从业人员，10.5%为临时从业人员。调查统计结果表明，较多渔民在中街山列岛海域、嵊泗列岛海域和浪岗山列岛海域捕获曼氏无针乌贼，渔民数量分别占调查总体的32.84%、24.51%和15.02% (图7)，该结果与模型预测结果总体相符。
3. 讨论
3.1 影响曼氏无针乌贼分布的主要生态因子
海表温度和海表盐度与头足类分布联系紧密[39]，因此在头足类分布与生态因子的关系研究中常被作为重要生态因子[40-41]。同时，溶解氧、叶绿素、深度、悬浮物和混合层深度也常被作为分析头足类分布的重要环境变量。目前针对曼氏无针乌贼分布与环境变量的关系，主要探讨了水温、盐度、叶绿素、水深和溶解氧与其栖息地选择之间的偏好关系。pH、初级生产力、营养盐含量、海流流速以及光照强度也与海洋生物分布关系密切。但针对这部分生态因子对曼氏无针乌贼分布影响的研究较少，因此本文在其他学者研究[42-44]的基础上进行补充研究。该过程涉及的生态因子较多，为避免主观选取生态因子导致模型预测结果产生误差，文章采取相关性分析法和刀切法筛选生态因子。模型运行结果表明，SST最大值、SSS最小值和pH对模型的贡献率较大，是影响曼氏无针乌贼的主要生态因子，该结果与目前已有研究的结果总体一致[5,42]。初级生产力、营养盐含量对曼氏无针乌贼分布具有一定影响，而溶解氧浓度、底层光照强度和海流流速对其分布的影响较小。
SST最大值、SSS最小值和pH生态因子的响应曲线表明，曼氏无针乌贼偏向在海表温度较高、盐度较大的环境活动，且能接受弱碱性的环境，舟山海域内SST最大值介于25.5~28.5 ℃、SSS最小值介于20‰~30‰、pH介于7.8~8.3时，曼氏无针乌贼存在的概率较大。该结论与实验室中观察不同温度、盐度、pH下曼氏无针乌贼活力的实验结果相符[45]，同时进一步验证了MaxEnt模型结果的可靠性。在舟山海域中，适宜曼氏无针乌贼生长的SST最大值介于25.5~28.5 ℃，且温度达到 (26±0.5) ℃时其存在概率可达0.75以上。符合在最佳温度范围内曼氏无针乌贼具有高生长速度和高成活率的特点[46]。Peng等[47]发现水温在25~27 ℃时，曼氏无针乌贼的瞬时生长速率、摄食速率、食物转化率和体长增加长度最佳。排除Peng等[47]在实验过程中所选种群 (浙南—闽东种群) 与舟山海域内曼氏无针乌贼 (浙北种群) 之间的差异，模型对曼氏无针乌贼在海表温度阈值的分析较为准确。SST最小值对模型的贡献率小于SST最大值 (表2)，表明曼氏无针乌贼分布受高温变化的影响较大，且在低温条件下曼氏无针乌贼的活动和摄食频率降低，即曼氏无针乌贼更适宜在温度较高的环境下生存。盐度有利于海洋生物维持渗透压，对曼氏无针乌贼幼体而言，稳定的盐度环境有利于其发育和成熟[48]。该特性解释了SSS最小值对曼氏无针乌贼分布的影响较大，且在海水盐度较低的条件下曼氏无针乌贼存在概率较小的原因。此外，SSS最小值达到25‰ 时，曼氏无针乌贼存在概率最大，表明曼氏无针乌贼对盐度的耐受性较强。曼氏无针乌贼存在概率与pH变化关系表明，曼氏无针乌贼对pH的耐受范围较窄，仅在pH为7.8~8.3时的存在概率较大。面对敌害或外界环境条件突变时，曼氏无针乌贼会出现喷墨现象，该现象在pH过高时也会发生[28]。结合该现象和物种在pH达8.2后其存在概率降幅较大的现象可知，曼氏无针乌贼对pH的变化敏感。因此后续对曼氏无针乌贼潜在适生区进行保护时需注意海洋水质及其酸碱度。
3.2 曼氏无针乌贼适生区分布及其修复策略
MaxEnt模型预测的物种适生区分布是模型在一定的研究区域和时间尺度下根据已有的约束条件 (目标物种的存在数据及其所处环境信息) 而运算出的物种时空分布拟合[49]，具有科学性和直观性，有利于为物种保护区规划研究和物种增殖放流选址提供参考。根据曼氏无针乌贼的适生区分布图及面积大小，可获得其在研究区域的分布特征。首先，曼氏无针乌贼在研究区域内的高适生区总体分布以嵊泗列岛、马鞍列岛和中街山列岛为核心。该特征与曼氏无针乌贼的繁殖特点有关，浙江北部的曼氏无针乌贼种群在需要进行交配和产卵时会游向浅海区和岛屿周围。舟山海域内岛屿众多，诸如嵊泗列岛、马鞍列岛和中街山列岛等岛屿为曼氏无针乌贼繁殖创造了优良条件。其次，曼氏无针乌贼适生区分布面积广，占研究区总面积的75.11%，大部分高适生区集中在舟山本岛的东北方海域，从马鞍列岛到舟山本岛东侧海域呈南北方向延伸，具有一定的连通性。且高适生区与次高适生区包含该海域的3个海洋特别保护区，其中在高适生区内包含2个国家级海洋特别保护区，分别为嵊泗马鞍列岛海洋保护区和普陀中街山列岛海洋生态特别保护区。次高适生区则包含了舟山市东部省级海洋保护区。利用该特征，加强对舟山海域内海洋特别保护区的保护，维持其生态系统平衡，有利于保护曼氏无针乌贼潜在适生区，从而间接促进其资源恢复。除自然条件外，资金和政策的支持也同样给予了曼氏无针乌贼资源恢复优势。2017年舟山市为海洋保护区提供法律支撑，出台了第一部国家级海洋特别保护区地方性法规《舟山国家级海洋特别保护区管理条例》。2018—2022年，为生态修复累计投入资金超过1亿元，陆续投放了近50万立方米的人工鱼礁[50]。因此，鉴于曼氏无针乌贼在舟山海域内适生区分布面积广、适生区保护程度高、资金投入充足及相关保护政策完善，可得出结论：舟山海域曼氏无针乌贼资源恢复的潜力较大。
根据适生区分布情况可采取相应策略进行曼氏无针乌贼增殖放流，从而促进其资源恢复。首先，可在曼氏无针乌贼产卵期选择高适生区投放曼氏无针乌贼成体，该措施有利于较大程度提高曼氏无针乌贼产量。原因在于，与其他海洋生物相比，头足类动物的生命周期相对较短，物种的数量变化受产卵数量和孵化数量的影响较大。并且曼氏无针乌贼产卵期较短，一般为4月中旬至6月下旬，亲体繁殖后不久便死亡[51]。其次，次高适生区、中适生区和低适生区3个适生等级的面积大小相近。从生境保护的角度出发，提高曼氏无针乌贼资源有两种策略，一是扩大曼氏无针乌贼高适生区面积，二是将次高适生区、中适生区和低适生区向更高一级的适生区转化。在面对次高适生区、中适生区和低适生区面积大小相近的情况时，选择次高适生区进行生境修复获得的效益较大。
3.3 MaxEnt模型在头足类中的应用
MaxEnt模型在头足类中的应用逐渐增加，本文以头足类曼氏无针乌贼为研究对象展现了MaxEnt模型在多元生态因子条件下模拟预测物种适生区分布及物种生境需求分析上的适用性。同时，与ArcGIS软件相结合获得了较为直观的结果，可为曼氏无针乌贼增殖放流选址提供参考。但该模型仍具有一定的局限性，在适生区等级划分上仍然具有一定的主观性。此外，头足类生物活动范围广，其分布不具有明确的分界线，而MaxEnt模型结果的准确性受样本量和研究范围的影响[52]，因此需考虑研究区的边界范围。综上，应用MaxEnt模型分析头足类时，应注重模型参数优化和模型预测结果的精度检验。大部分头足类具有洄游分布的特点，根据其洄游规律划分时间尺度，再利用MaxEnt模型进行物种生境适宜性分析可获得更为精准的物种生境适宜区域分布情况。由于调查过程中受疫情影响，部分月份收集统计到的曼氏无针乌贼分布点较少。为了保证模型的预测精度，本文未将时间尺度划分至月份后再进行曼氏无针乌贼的生境适宜性分析。在今后的研究中，可进一步补充相关数据，探索不同时空尺度下曼氏无针乌贼的生境适宜性，或选择不同物种分布模型与MaxEnt模型的预测结果进行对比，多方位、多角度地为曼氏无针乌贼的资源养护提供科学参考和理论支持。
-
![]()
图 1 Logistic方程拟合雌、雄全长区间性成熟百分比求得的初次性成熟全长 (L50)
Figure 1. L50 for female and male A. schlegelii by ratio of interval maturity of full length fitted by Logistic equation
![]()
图 2 黑鲷各性别鱼在繁殖期内的体质量、全长分布频率
Figure 2. Distribution frequency of body mass and total length of female and male A. schlegelii during reproductive period
![]()
图 4 大亚湾黑鲷雌、雄鱼性成熟系数 (平均值±标准差) 的周年变化
数字代表样本个数，上面一列为雌性，下面一列为雄性
Figure 4. Annual variation of gonadosomatic indices (
${\overline X}\pm{\rm{SD}}$ ) for female and male A. schlegeli in Daya BayThe number indicates the quantity of samples, female above and male below.
![]()
图 5 大亚湾黑鲷雌、雄鱼每月性腺阶段的变化
数字代表样本个数
Figure 5. Monthly variation in gonad stages for female and male A. schlegeli in Daya Bay
The number indicates the quantity of samples.
-
[1] 吴仁协, 刘静, 樊冀蓉, 等. 黑棘鲷的命名和分类地位探究[J]. 海洋科学, 2011, 35(5): 117-119. [2] 林金, 陈涛, 陈琳, 等. 大亚湾黑鲷标志放流技术[J]. 水产学报, 2001, 25(1): 79-83. [3] 江世贵, 苏天凤, 夏军红, 等. 中国近海鲷科鱼类种质资源及其利用[M]. 北京: 海洋出版社, 2012: 17-18. [4] 国家海洋局第三研究所. 大亚湾海洋生态文集(Ⅰ)[C]. 北京: 海洋出版社, 1989: 1-5. [5] 国家海洋局第三研究所. 大亚湾海洋生态文集(Ⅱ)[C]. 北京: 海洋出版社, 1990: 1-42. [6] 邱永松, 曾小光, 陈涛, 等. 南海渔业资源与渔业管理[M]. 北京: 海洋出版社, 2008: 115-116. [7] 广东省海岛资源综合调查大队, 广东省海岸带和海涂资源综合调查领导小组办公室. 广东省海岛资源综合调查报告[R]. 广州: 广东科技出版社, 1995: 121-122. [8] 柯东胜, 彭晓娟, 吴玲铃, 等. 大亚湾典型生态系统状况调查与分析[J]. 海洋环境科学, 2009, 28(4): 421-425. doi: 10.3969/j.issn.1007-6336.2009.04.016 [9] 王雪辉, 杜飞雁, 邱永松, 等. 1980-2007年大亚湾鱼类物种多样性、区系特征和数量变化[J]. 应用生态学报, 2010, 21(9): 2403-2410. [10] 王家玥, 林清, 周睿东, 等. 广东省环境水体中的氚活度浓度水平[J]. 辐射防护, 2017, 37(4): 265-273. [11] 黎广媚. 疏浚泥海洋倾倒对大亚湾口海域环境影响的研究[D].青岛: 中国海洋大学, 2004: 1-62. [12] 孔祥昆. 由大亚湾核电站外海面油污事故引发的思考[J]. 船海工程, 2007, 36(5): 136-138. doi: 10.3963/j.issn.1671-7953.2007.05.042 [13] 杨玉峰, 梁浩亮, 范开文, 等. 惠州海域沉积物重金属污染与潜在生态风险[J]. 海岸工程, 2018, 37(1): 80-88. doi: 10.3969/j.issn.1002-3682.2018.01.009 [14] 梁庆洋, 齐占会, 巩秀玉, 等. 大亚湾鱼类深水网箱养殖对环境的影响[J]. 南方水产科学, 2017, 23(5): 25-32. doi: 10.3969/j.issn.2095-0780.2017.05.004 [15] 黄洪辉, 林钦, 王文质, 等. 大鹏澳海水鱼类网箱养殖对水环境的影响[J]. 南方水产, 2005, 1(3): 9-17. doi: 10.3969/j.issn.2095-0780.2005.03.002 [16] 李磊, 蒋玫, 沈新强, 等. 0#柴油和原油水溶性成分对黑鲷肝脏7-乙氧基异吩噁唑酮-脱乙基酶活性及其mRNA表达的影响[J]. 应用与环境生物学报, 2014, 20(6): 1027-1032. [17] 吕敢堂, 王志铮, 邵国洱, 等. 4种重金属离子对黑鲷幼鱼的急性毒性研究[J]. 浙江海洋学院学报(自然科学版), 2010, 29(3): 206-210. doi: 10.3969/j.issn.1008-830X.2010.03.002 [18] 蒋玫, 林钦, 沈新强. 航道疏浚泥悬浮液对黑鲷胚胎的发育及仔鱼摄食生长的影响[J]. 生态环境学报, 2009, 18(5): 1674-1678. doi: 10.3969/j.issn.1674-5906.2009.05.013 [19] 殷名称. 鱼类生态学[M]. 基隆: 水产出版社, 1998: 105-128. [20] LAW C, MITCHESON Y. Reproductive biology of black seabream Acanthopagrus schlegelii, threadfin porgy Evynnis cardinalis and red pargo Pagrus major in the northern South China Sea with consideration of fishery status and management needs[J]. J Fish Biol, 2017, 91(1): 101-125. doi: 10.1111/jfb.2017.91.issue-1
[21] 铃木克美, 熊国强. 硬骨鱼类的雌雄同体现象[J]. 水产科技情报, 1990(5): 151-154. [22] 孙远东, 谭立军, 唐新科, 等. 鱼类性别决定及其基因的研究进展[J]. 安徽农业科学, 2008, 36(18): 7691-7692, 7695. doi: 10.3969/j.issn.0517-6611.2008.18.083 [23] WU G C, CHANG C F. Oocytes survive in the testis by altering the soma fate from male to female in the protandrous black porgy, Acanthopagrus schlegeli[J]. Biol Reprod, 2013, 88(1): 1-9.
[24] 陈淑吟, 张志勇, 刘海林, 等. 池养黑鲷性逆转组织学观察[J]. 海洋科学, 2017, 41(7): 51-58. [25] WOONTTON R J. Ecology of teleost fishes[M]. London: Chapman an Hall, 1991: 159-195.
[26] 林昭进, 王雪辉, 江艳娥. 大亚湾鱼卵数量分布及种类组成特征[J]. 中国水产科学, 2010, 17(3): 543-550. [27] 邓利, 张为民, 林浩然, 等. 黑鲷生长激素及其受体的季节变化[J]. 水产学报, 2001, 25(3): 203-208. [28] 胡兴华. 鲷鱼人工繁殖研究-嘉腊Chrysophrys major及黑鲷Acanthopagrus schlegeli之探讨[C].台湾省水产研究所澎湖分所实验报告汇集, 1983(3): 1-48. [29] 郑镇安. 黑鲷人工繁殖及育苗的研究[J]. 渔业研究, 1983(3): 3-16. [30] 蔡兴邦. 黑鲷自然产卵及苗种培育的初步观察[J]. 江苏海洋水产, 1983(5): 63-71. [31] 张孝威. 黑鲷卵子及仔, 稚, 幼鱼的形态观察[J]. 动物学报, 1980, 26(4): 331-336. [32] 李加儿, 区又君, 许波涛, 等. 黑鲷生物学及其种苗生产技术研究[J]. 南海水产研究, 1990(2): 36-44. [33] 岳敏娟, 尤永隆, 林丹军. 温度对鲫鱼性腺分化的影响[J]. 动物学杂志, 2009, 44(1): 9-16. doi: 10.3969/j.issn.0250-3263.2009.01.002 [34] 施兆鸿. 盐度对黑鲷卵巢发育的影响[J]. 水产学报, 1996, 20(4): 357-360. [35] CAMBRAY J A, BRUTON M N. The reproductive strategy of a barb, Barbus anoplus (Pisces: Cyprinidae), colonizing a man-made lake in South Africa[J]. J Zool, 2010, 204(2): 143-168.
[36] 严云志. 抚仙湖外来鱼类生活史对策的适应性进化研究[D]. 武汉: 中国科学院研究生院(水生生物研究所), 2005: 21-105. [37] 中国农业百科全书总编辑委员会水产业卷编辑委员会. 中国农业百科全书:水产业卷上[M]. 北京: 中国农业出版社, 1994: 216-217. [38] 郭禹. 花尾胡椒鲷与黑鲷的超声波标志跟踪研究[D]. 大连: 大连海洋大学, 2016: 18-31. [39] 韩兴勇. 渔业资源保护和渔民的意识关系[J]. 中国渔业经济, 2000, 18(6): 29-30. doi: 10.3969/j.issn.1009-590X.2000.06.014 [40] 唐议, 李富荣, 庄会柏, 等. 我国渔民弱势群体问题与对策分析[J]. 中国渔业经济, 2006, 24(4): 3-7. doi: 10.3969/j.issn.1009-590X.2006.04.001 [41] 邹玉芹, 张东芝, 王为璋. 黑鲷网箱养殖技术[J]. 齐鲁渔业, 2000, 17(3): 21-22. [42] 朱德芬. 黑鲷人工养殖技术讲座: 第一讲黑鲷生物学特性及增养殖概况[J]. 水产养殖, 1996(1): 30-31. [43] 钱国英. 鱼类的性别及人工控制[J]. 动物学杂志, 2000, 35(1): 47-52. doi: 10.3969/j.issn.0250-3263.2000.01.019 [44] 邹记兴, 陶友宝, 向文洲, 等. 人工诱导点带石斑鱼性逆转的组织学证据及其机制探讨[J]. 高技术通讯, 2003, 13(6): 81-86. doi: 10.3321/j.issn:1002-0470.2003.06.019 [45] 储张杰. 黄鳝性逆转研究进展[J]. 水利渔业, 2006, 26(6): 19-23. doi: 10.3969/j.issn.1003-1278.2006.06.007 -
期刊类型引用(3)
1. 吴春梅，杨柳明，王维奇，罗欣婕. 闽江河口湿地围垦后土壤微生物残体碳特征. 环境科学学报. 2024(12): 360-368 . 
百度学术
2. 洪妍，杨平，仝川，赵光辉，李玲，唐晨，张彦，谭莹莹，阮蔓菁. 亚热带河口区对虾养殖池塘浮游植物初级生产力变化. 湖泊科学. 2022(03): 881-893 . 百度学术
3. 林于蓝，陈钰，尹晓雷，曾庆淞，王维奇. 围垦养殖与退塘还湿对闽江河口湿地土壤铁碳结合特征的影响. 环境科学学报. 2022(07): 466-477 . 百度学术
其他类型引用(4)
下载:
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 7
粤公网安备 44010502001741号
