利用MaxEnt模型的舟山海域曼氏无针乌贼生境适宜性研究

张漫瑶, 石雯静, 桂峰, 曾宪敏, 徐开达, 赵晟

张漫瑶, 石雯静, 桂峰, 曾宪敏, 徐开达, 赵晟. 利用MaxEnt模型的舟山海域曼氏无针乌贼生境适宜性研究[J]. 南方水产科学, 2023, 19(5): 22-31. DOI: 10.12131/20230080
引用本文: 张漫瑶, 石雯静, 桂峰, 曾宪敏, 徐开达, 赵晟. 利用MaxEnt模型的舟山海域曼氏无针乌贼生境适宜性研究[J]. 南方水产科学, 2023, 19(5): 22-31. DOI: 10.12131/20230080
ZHANG Manyao, SHI Wenjing, GUI Feng, ZENG Xianmin, XU Kaida, ZHAO Sheng. Study of suitable habitats for Sepiella maindroni in Zhoushan sea areas based on MaxEnt model[J]. South China Fisheries Science, 2023, 19(5): 22-31. DOI: 10.12131/20230080
Citation: ZHANG Manyao, SHI Wenjing, GUI Feng, ZENG Xianmin, XU Kaida, ZHAO Sheng. Study of suitable habitats for Sepiella maindroni in Zhoushan sea areas based on MaxEnt model[J]. South China Fisheries Science, 2023, 19(5): 22-31. DOI: 10.12131/20230080

利用MaxEnt模型的舟山海域曼氏无针乌贼生境适宜性研究

基金项目: 国家重点研发计划“蓝色粮仓科技创新”专项 (2019YFD0901204-04);浙江省属高校基本科研业务费 (2021JD006);浙江省重点研发计划项目 (2019C02056)
详细信息
    作者简介:

    张漫瑶 (1997—),女,硕士研究生,研究方向为海域生境修复。E-mail: moyna1213@163.com

    通讯作者:

    赵 晟 (1973—),男,教授,博士,研究方向为海洋生态学。E-mail: zhaosh@zjou.edu.cn

  • 中图分类号: S 922.9+3

Study of suitable habitats for Sepiella maindroni in Zhoushan sea areas based on MaxEnt model

  • 摘要:

    曼氏无针乌贼 (Sepiella maindroni) 具有较高的经济和营养价值,近年因过度捕捞、生境破坏等导致其资源衰退。为了定量评估多元生态因子对曼氏无针乌贼分布的影响,以曼氏无针乌贼的重要产卵地和栖息地—舟山海域为研究区域,基于2017—2021年浙江近海的底拖网调查数据、2021—2022年自主航次的调查数据和全球生物多样性信息机构 (Global Biodiversity Information Facility, GBIF) 数据,利用最大熵 (Maximum Entropy, MaxEnt) 模型分析曼氏无针乌贼的生境适宜性,探索其修复潜力和适生区的空间分布格局。结果表明,海表温度最大值、海表盐度最小值和pH是影响曼氏无针乌贼分布的主要生态因子。曼氏无针乌贼偏向于在海洋表面温度较高、盐度较大的海域活动,且能接受弱碱性环境。当舟山海域内海表温度最大值介于25.5~28.5 ℃、海表盐度最小值介于20‰~30‰、pH介于7.8~8.3时,曼氏无针乌贼存在的概率较大。曼氏无针乌贼在舟山海域的适生区面积分布较广 (30 213.40 km2),占研究区域总面积的75.11%。其中高适生区总体分布在中街山列岛、嵊泗列岛和马鞍列岛,面积为4 600.78 km2。次高适生区、中适生区和低适生区的面积分别为8 205.52、9 723.31和7 683.78 km2。从修复潜力上看,曼氏无针乌贼在舟山海域有较大的恢复潜力。

    Abstract:

    Cuttlefish (Sepiella maindroni) has high economic and nutritional value, but its resources have declined in recent years due to overfishing and habitat destruction. To quantitatively assess the impact of multiple ecological factors on the distribution of S. maindroni, we selected the Zhoushan sea area which is an important spawning ground and habitat for S. maindroni as the research area, and used the MaxEnt (Maximum Entropy) model to analyze the habitat suitability of S. maindroni as well as to explore its restoration potential and spatial distribution of suitable habitats, based on the bottom trawl survey data from the Zhejiang coastal area from 2017−2021, the survey data from autonomous cruises in 2021−2022, and data from the Global Biodiversity Information Facility (GBIF). The research results show that the maximum sea surface temperature (SST), minimum sea surface salinity (SSS) and pH were the main ecological factors affecting the S. maindroni distribution. S. maindroni tended to be active in areas with higher sea surface temperatures and higher salinity, and could tolerate a weakly alkaline environment. In the Zhoushan sea area, S. maindroni was more likely to exist with maximum SST of 25.5−28.5 ℃, minimum SSS of 20‰−30‰, and pH value of 7.8−8.3. The suitable habitat area in the Zhoushan sea area was quite extensive, covering 30 213.40 km2, which accounted for 75.11% of the total area of the research area. Among these, the highly suitable habitats were mainly distributed in the sea areas of Zhongjieshan Islands, Shengsi Islands and Ma'an Islands, with an area of 4 600.78 km2. The areas of the sub-highly suitable, moderately suitable and lowly suitable habitats were 8 205.52, 9 723.31 and 7 683.78 km2, respectively. In terms of restoration potential, S. maindroni has significant recovery potential in the Zhoushan sea area.

  • 增殖放流是目前有效恢复渔业资源的方法之一[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模型对曼氏无针乌贼进行生境适宜性评价,以期为曼氏无针乌贼的增殖放流选址提供参考,并为其资源养护提供科学支撑。

    构建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  曼氏无针乌贼舟山海域的分布点
    Figure  1.  Scatter plot of S. maindroni in Zhoushan sea area

    海洋环境数据集来源于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
    下载: 导出CSV 
    | 显示表格

    导入数据至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格式的栅格数据输出。

    参考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)。

    本文使用受试者工作特征曲线 (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  ROC 曲线检验模型结果
    Figure  2.  Model results tested by ROC curve

    表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
    下载: 导出CSV 
    | 显示表格

    生态因子重要性结果显示 (图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时,曼氏无针乌贼存在的概率较大。

    图  4  生态因子响应曲线
    Figure  4.  Response curves for ecological factors

    在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升高呈下降趋势。

    利用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
    下载: 导出CSV 
    | 显示表格

    为检验模型结果,采取访问调查的方式走访舟山市各区县典型渔村,开展曼氏无针乌贼捕捞出现海域调研。范围包括定海区 (金塘镇、长白乡)、普陀区 (沈家门、朱家尖、六横岛、桃花镇、东极岛)、岱山县 (衢山镇、长涂镇) 和嵊泗县 (嵊山镇、枸杞乡、洋山镇)。据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),该结果与模型预测结果总体相符。

    图  7  渔民捕捞曼氏无针乌贼海域的调研结果
    Figure  7.  Results of a survey on fishing area of S. maindroni

    海表温度和海表盐度与头足类分布联系紧密[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的变化敏感。因此后续对曼氏无针乌贼潜在适生区进行保护时需注意海洋水质及其酸碱度。

    MaxEnt模型预测的物种适生区分布是模型在一定的研究区域和时间尺度下根据已有的约束条件 (目标物种的存在数据及其所处环境信息) 而运算出的物种时空分布拟合[49],具有科学性和直观性,有利于为物种保护区规划研究和物种增殖放流选址提供参考。根据曼氏无针乌贼的适生区分布图及面积大小,可获得其在研究区域的分布特征。首先,曼氏无针乌贼在研究区域内的高适生区总体分布以嵊泗列岛、马鞍列岛和中街山列岛为核心。该特征与曼氏无针乌贼的繁殖特点有关,浙江北部的曼氏无针乌贼种群在需要进行交配和产卵时会游向浅海区和岛屿周围。舟山海域内岛屿众多,诸如嵊泗列岛、马鞍列岛和中街山列岛等岛屿为曼氏无针乌贼繁殖创造了优良条件。其次,曼氏无针乌贼适生区分布面积广,占研究区总面积的75.11%,大部分高适生区集中在舟山本岛的东北方海域,从马鞍列岛到舟山本岛东侧海域呈南北方向延伸,具有一定的连通性。且高适生区与次高适生区包含该海域的3个海洋特别保护区,其中在高适生区内包含2个国家级海洋特别保护区,分别为嵊泗马鞍列岛海洋保护区和普陀中街山列岛海洋生态特别保护区。次高适生区则包含了舟山市东部省级海洋保护区。利用该特征,加强对舟山海域内海洋特别保护区的保护,维持其生态系统平衡,有利于保护曼氏无针乌贼潜在适生区,从而间接促进其资源恢复。除自然条件外,资金和政策的支持也同样给予了曼氏无针乌贼资源恢复优势。2017年舟山市为海洋保护区提供法律支撑,出台了第一部国家级海洋特别保护区地方性法规《舟山国家级海洋特别保护区管理条例》。2018—2022年,为生态修复累计投入资金超过1亿元,陆续投放了近50万立方米的人工鱼礁[50]。因此,鉴于曼氏无针乌贼在舟山海域内适生区分布面积广、适生区保护程度高、资金投入充足及相关保护政策完善,可得出结论:舟山海域曼氏无针乌贼资源恢复的潜力较大。

    根据适生区分布情况可采取相应策略进行曼氏无针乌贼增殖放流,从而促进其资源恢复。首先,可在曼氏无针乌贼产卵期选择高适生区投放曼氏无针乌贼成体,该措施有利于较大程度提高曼氏无针乌贼产量。原因在于,与其他海洋生物相比,头足类动物的生命周期相对较短,物种的数量变化受产卵数量和孵化数量的影响较大。并且曼氏无针乌贼产卵期较短,一般为4月中旬至6月下旬,亲体繁殖后不久便死亡[51]。其次,次高适生区、中适生区和低适生区3个适生等级的面积大小相近。从生境保护的角度出发,提高曼氏无针乌贼资源有两种策略,一是扩大曼氏无针乌贼高适生区面积,二是将次高适生区、中适生区和低适生区向更高一级的适生区转化。在面对次高适生区、中适生区和低适生区面积大小相近的情况时,选择次高适生区进行生境修复获得的效益较大。

    MaxEnt模型在头足类中的应用逐渐增加,本文以头足类曼氏无针乌贼为研究对象展现了MaxEnt模型在多元生态因子条件下模拟预测物种适生区分布及物种生境需求分析上的适用性。同时,与ArcGIS软件相结合获得了较为直观的结果,可为曼氏无针乌贼增殖放流选址提供参考。但该模型仍具有一定的局限性,在适生区等级划分上仍然具有一定的主观性。此外,头足类生物活动范围广,其分布不具有明确的分界线,而MaxEnt模型结果的准确性受样本量和研究范围的影响[52],因此需考虑研究区的边界范围。综上,应用MaxEnt模型分析头足类时,应注重模型参数优化和模型预测结果的精度检验。大部分头足类具有洄游分布的特点,根据其洄游规律划分时间尺度,再利用MaxEnt模型进行物种生境适宜性分析可获得更为精准的物种生境适宜区域分布情况。由于调查过程中受疫情影响,部分月份收集统计到的曼氏无针乌贼分布点较少。为了保证模型的预测精度,本文未将时间尺度划分至月份后再进行曼氏无针乌贼的生境适宜性分析。在今后的研究中,可进一步补充相关数据,探索不同时空尺度下曼氏无针乌贼的生境适宜性,或选择不同物种分布模型与MaxEnt模型的预测结果进行对比,多方位、多角度地为曼氏无针乌贼的资源养护提供科学参考和理论支持。

  • 图  1   曼氏无针乌贼舟山海域的分布点

    Figure  1.   Scatter plot of S. maindroni in Zhoushan sea area

    图  2   ROC 曲线检验模型结果

    Figure  2.   Model results tested by ROC curve

    图  3   刀切法检验生态因子对曼氏无针乌贼分布影响的重要程度

    Figure  3.   Importance of ecological factors on distribution of S. maindroni by Jackknife test

    图  4   生态因子响应曲线

    Figure  4.   Response curves for ecological factors

    图  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

    图  7   渔民捕捞曼氏无针乌贼海域的调研结果

    Figure  7.   Results of a survey on fishing area of S. maindroni

    表  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
    下载: 导出CSV

    表  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
    下载: 导出CSV

    表  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
    下载: 导出CSV
  • [1]

    RICHARDSON L E, LENFANT P, CLARKE L J, et al. Examining current best-practices for the use of wild post-larvae capture, culture, and release for fisheries enhancement[J]. Front Mar Sci, 2023, 9: 1058497. doi: 10.3389/fmars.2022.1058497

    [2] 钟明, 侍昊, 安树青, 等. 中国野生动物生境适宜性评价和生境破碎化研究[J]. 生态科学, 2016, 35(4): 205-209. doi: 10.14108/j.cnki.1008-8873.2016.04.028
    [3] 吕国敏, 吴进锋, 陈利雄. 我国头足类增养殖研究现状及开发前景[J]. 南方水产, 2007, 3(3): 61-66.
    [4]

    HEIKKINEN R K, LUOTO M, ARAÚJO M B, et al. Methods and uncertainties in bioclimatic envelope modelling under climate change[J]. Prog Phys Geogr, 2006, 30(6): 751-777. doi: 10.1177/0309133306071957

    [5] 徐晓萱, 谢玉, 刘姝含, 等. 基于GAM模型的浙江近海曼氏无针乌贼时空分布研究[J]. 浙江海洋大学学报 (自然科学版), 2022, 41(5): 400-407.
    [6] 魏裙倚, 崔国辰, 玄文丹, 等. 海表面温度及叶绿素a浓度对西北印度洋鸢乌贼时空分布的影响[J]. 中国水产科学, 2022, 29(3): 388-397.
    [7] 李杰, 张鹏, 晏磊, 等. 南海中南部海域鸢乌贼CPUE影响因素的GAM分析[J]. 中国水产科学, 2020, 27(8): 906-915.
    [8] 杨洋. 基于四种模型的肉苁蓉潜在地理分布预测及空间格局变化分析[D]. 西安: 陕西师范大学, 2017: 25-28.
    [9] 图雅. 基于物种分布模型的中国针茅属植物潜在分布区及其影响因子分析[D]. 北京: 北京林业大学, 2020: 72-73.
    [10]

    MEROW C, SILANDER J A. A comparison of Maxlike and Maxent for modelling species distributions[J]. Methods Ecol Evol, 2014, 5(3): 215-225. doi: 10.1111/2041-210X.12152

    [11] 陈新美, 雷渊才, 张雄清, 等. 样本量对MaxEnt模型预测物种分布精度和稳定性的影响[J]. 林业科学, 2012, 48(1): 53-59. doi: 10.11707/j.1001-7488.20120110
    [12]

    TENG S, SU N, LEE M, et al. Modeling the habitat distribution of Acanthopagrus schlegelii in the coastal waters of the Eastern Taiwan Strait using MaxEnt with fishery and remote sensing data[J]. J Mar Sci Eng, 2021, 9(12): 1442. doi: 10.3390/jmse9121442

    [13]

    HARTE J, NEWMAN E A. Maximum information entropy: a foundation for ecological theory[J]. Trends Ecol Evol, 2014, 29(7): 384-389. doi: 10.1016/j.tree.2014.04.009

    [14] 杨继超, 董民星, 种衍飞, 等. 最大熵模型在海洋生物适生区预测中的应用[J/OL]. 应用海洋学学报. https://kns.cnki.net/kcms/detail/35.1319.P.20230302.1347.002.html.
    [15] 周海涛, 那晓东, 臧淑英, 等. 最大熵 (MaxEnt) 模型在物种栖息地研究中的应用[J]. 环境科学与管理, 2016, 41(3): 149-151. doi: 10.3969/j.issn.1673-1212.2016.03.035
    [16] 张嘉容, 杨晓明, 田思泉. 基于最大熵模型的南太平洋长鳍金枪鱼栖息地预测[J]. 中国水产科学, 2020, 27(10): 1222-1233.
    [17] 龚彩霞, 陈新军, 高峰, 等. 水温变暖对西北太平洋柔鱼潜在栖息地分布的影响[J]. 海洋学报 (中文版), 2022, 44(7): 95-102.
    [18]

    PIRTLE J L, SHOTWELL S K, ZIMMERMANN M, et al. Habitat suitability models for groundfish in the Gulf of Alaska[J]. Deep-Sea Res II, 2017, 165: 303-321.

    [19] 曹子豪, 迟长凤, 刘慧慧, 等. 不同地理群体曼氏无针乌贼肌肉营养成分分析比较与评价[J]. 食品科学, 2015, 36(4): 101-105. doi: 10.7506/spkx1002-6630-201504019
    [20] 倪正雅, 徐汉祥. 浙江近海乌贼资源评估及乌贼渔业管理[J]. 海洋渔业, 1986(2): 51-54.
    [21] 李继姬, 郭宝英, 吴常文. 浙江海域曼氏无针乌贼资源演变及修复路径探讨[J]. 浙江海洋学院学报 (自然科学版), 2011, 30(5): 381-385, 396.
    [22] 覃涛, 俞存根, 陈全震, 等. 舟山渔场及邻近海域头足类 (Cephalopod) 种类组成和数量分布[J]. 海洋与湖沼, 2011, 42(1): 124-130. doi: 10.11693/hyhz201101019019
    [23] 陈伟峰, 叶深, 余玥, 等. 浙南近海头足类种类组成及生态位分析[J]. 水生生物学报, 2021, 45(2): 428-435. doi: 10.7541/2021.2018.262
    [24]

    PANG Y M, TIAN Y J, FU C H, et al. Variability of coastal cephalopods in overexploited China seas under climate change with implications on fisheries management[J]. Fish Res, 2018, 208: 22-33. doi: 10.1016/j.fishres.2018.07.004

    [25] 张义浩, 王志铮, 吴常文, 等. 舟山群岛定生海藻种类组成、生态分布及区系特征研究[J]. 浙江海洋学院学报 (自然科学版), 2002, 21(2): 98-105.
    [26] 刘连为, 隋宥珍, 徐开达, 等. 曼氏无针乌贼两个地理群体繁殖生物学比较研究[J]. 海洋湖沼通报, 2022, 44(4): 57-62. doi: 10.13984/j.cnki.cn37-1141.2022.04.008
    [27]

    WARREN D L, MATZKE N J, CARDILLO M, et al. ENMTools 1.0: an R package for comparative ecological biogeography[J]. Ecography, 2021, 44(4): 504-511. doi: 10.1111/ecog.05485

    [28] 尹飞, 王春琳, 宋微微. 曼氏无针乌贼幼体生态因子耐受性的研究[J]. 湛江海洋大学学报, 2005, 25(4): 39-43.
    [29] 黄伟卿, 陈宇光, 张艺, 等. 盐度胁迫对曼氏无针乌贼胚胎发育和生长性能的影响[J]. 渔业现代化, 2022, 49(6): 100-107. doi: 10.3969/j.issn.1007-9580.2022.06.013
    [30] 郑美丽, 肖金华, 郑微云, 等. 曼氏无针乌贼的趋光特性[J]. 厦门大学学报 (自然科学版), 1980(3): 91-99.
    [31]

    GUISAN A, ZIMMERMANN N E. Predictive habitat distribution models in ecology[J]. Ecol Model, 2000, 135(2): 147-186.

    [32] 谢平, 赵羽西, 桑燕芳, 等. 基于相关系数的水文周期变异分级方法及验证[J]. 水力发电学报, 2018, 37(12): 33-43. doi: 10.11660/slfdxb.20181204
    [33]

    COBOS M E, PETERSON A T, BARVE N, et al. kuenm: an R package for detailed development of ecological niche models using MaxEnt[J]. PeerJ, 2019, 7: e6281. doi: 10.7717/peerj.6281

    [34]

    SHCHEGLOVITOVA M, ANDERSON R P. Estimating optimal complexity for ecological niche models: a jackknife approach for species with small sample sizes[J]. Ecol Model, 2013, 269: 9-17. doi: 10.1016/j.ecolmodel.2013.08.011

    [35] 阮欧, 刘绥华, 陈芳, 等. 基于多源遥感的贵州草海国家级自然保护区黑颈鹤生境适宜性评价[J]. 生态学报, 2022, 42(5): 1947-1957.
    [36]

    MASTRANDREA M D, FIELD C B, STOCKER T F, et al. Guidance Note for Lead Authors of the IPCC Fifth Assessment Report on Consistent Treatment of Uncertainties[C]. The IPCC Panel, Intergovernmental Panel on Climate Change (IPCC), CA, USA, IPCC Cross-Working Group Meeting on Consistent Treatment of Uncertainties Jasper Ridge, 2010: 2-3.

    [37]

    PETERSON A T, PAPE M, SOBERÓN J. Rethinking receiver operating characteristic analysis applications in ecological niche modeling[J]. Ecol Model, 2017, 213(1): 63-72.

    [38] 王运生, 谢丙炎, 万方浩, 等. ROC曲线分析在评价入侵物种分布模型中的应用[J]. 生物多样性, 2007, 15(4): 365-372. doi: 10.3321/j.issn:1005-0094.2007.04.005
    [39] 曹杰, 陈新军, 刘必林, 等. 鱿鱼类资源量变化与海洋环境关系的研究进展[J]. 上海海洋大学学报, 2010, 19(2): 232-239.
    [40] 肖戈, 徐博, 张衡, 等. 阿拉伯海外海鸢乌贼渔场时空分布与海洋环境要素的研究[J]. 南方水产科学, 2022, 18(4): 10-19. doi: 10.12131/20210217
    [41] 陈峰, 李楠, 方舟, 等. 浙江近岸海域春夏季剑尖枪乌贼栖息地分布变化规律[J]. 上海海洋大学学报, 2021, 30(5): 847-855. doi: 10.12024/jsou.20201203245
    [42] 吴天, 梁君, 周永东, 等. 中街山列岛海域曼氏无针乌贼栖息地偏向性研究[J]. 浙江海洋大学学报 (自然科学版), 2022, 41(5): 408-417.
    [43] 刘姝含, 刘连为, 徐开达, 等. 浙江中北部近海曼氏无针乌贼资源现状[J]. 浙江海洋大学学报 (自然科学版), 2022, 41(4): 286-293.
    [44] 袁健美, 张虎, 贲成恺, 等. 江苏近海曼氏无针乌贼增殖放流效果评估[J]. 水产养殖, 2023, 44(2): 12-16. doi: 10.3969/j.issn.1004-2091.2023.02.003
    [45] 张建设. 曼氏无针乌贼养殖生物学特性和血细胞免疫功能研究[D]. 厦门: 厦门大学, 2007: 155.
    [46]

    YIN F, SUN P, PENG S M, et al. The respiration, excretion and biochemical response of the juvenile common Chinese cuttlefish, Sepiella maindroni at different temperatures[J]. Aquaculture, 2013, 402/403: 127-132. doi: 10.1016/j.aquaculture.2013.03.018

    [47]

    PENG S M, YIN F, SHI Z H, et al. Optimum water temperature for the growth of juvenile common Chinese cuttlefish, Sepiella maindroni (de Rochebrune 1884)[J]. J Shellfish Res, 2011, 30(2): 205-209. doi: 10.2983/035.030.0202

    [48] 雷舒涵, 吴常文, 高天翔, 等. 金乌贼和曼氏无针乌贼胚胎发育及其盐度耐受能力的比较研究[J]. 中国水产科学, 2011, 18(2): 350-359.
    [49]

    HALVORSEN R, MAZZONI S, BRYN A, et al. Opportunities for improved distribution modelling practice via a strict maximum likelihood interpretation of MaxEnt[J]. Ecography, 2015, 38(2): 172-183. doi: 10.1111/ecog.00565

    [50] 朱海洋. 舟山渔场: 从“赶海”到“养海”[N]. 农民日报, 2023-02-16(6).
    [51] 宋旻鹏, 汪金海, 郑小东. 中国经济头足类增养殖现状及展望[J]. 海洋科学, 2018, 42(3): 149-156. doi: 10.11759/hykx20180128001
    [52] 季乾昭, 王荣兴, 黄志旁, 等. 样本量与研究范围变化对MaxEnt模型准确度的影响: 以黑白仰鼻猴为例[J]. 兽类学报, 2019, 39(2): 126-133.
  • 期刊类型引用(3)

    1. 吴春梅,杨柳明,王维奇,罗欣婕. 闽江河口湿地围垦后土壤微生物残体碳特征. 环境科学学报. 2024(12): 360-368 . 百度学术
    2. 洪妍,杨平,仝川,赵光辉,李玲,唐晨,张彦,谭莹莹,阮蔓菁. 亚热带河口区对虾养殖池塘浮游植物初级生产力变化. 湖泊科学. 2022(03): 881-893 . 百度学术
    3. 林于蓝,陈钰,尹晓雷,曾庆淞,王维奇. 围垦养殖与退塘还湿对闽江河口湿地土壤铁碳结合特征的影响. 环境科学学报. 2022(07): 466-477 . 百度学术

    其他类型引用(4)

图(7)  /  表(3)
计量
  • 文章访问数:  396
  • HTML全文浏览量:  136
  • PDF下载量:  78
  • 被引次数: 7
出版历程
  • 收稿日期:  2023-04-16
  • 修回日期:  2023-05-17
  • 录用日期:  2023-06-05
  • 网络出版日期:  2023-06-09
  • 刊出日期:  2023-10-04

目录

/

返回文章
返回