基于MaxEnt模型模拟中西太平洋鲣自由鱼群栖息地的研究

汪伟松, 唐未, 龚一赫, 王学昉, 李玉伟

汪伟松, 唐未, 龚一赫, 王学昉, 李玉伟. 基于MaxEnt模型模拟中西太平洋鲣自由鱼群栖息地的研究[J]. 南方水产科学, 2023, 19(5): 11-21. DOI: 10.12131/20230011
引用本文: 汪伟松, 唐未, 龚一赫, 王学昉, 李玉伟. 基于MaxEnt模型模拟中西太平洋鲣自由鱼群栖息地的研究[J]. 南方水产科学, 2023, 19(5): 11-21. DOI: 10.12131/20230011
WANG Weisong, TANG Wei, GONG Yihe, WANG Xuefang, LI Yuwei. Modeling habitat of skipjack tuna of free swimming school in Western and Central Pacific Ocean based on MaxEnt model[J]. South China Fisheries Science, 2023, 19(5): 11-21. DOI: 10.12131/20230011
Citation: WANG Weisong, TANG Wei, GONG Yihe, WANG Xuefang, LI Yuwei. Modeling habitat of skipjack tuna of free swimming school in Western and Central Pacific Ocean based on MaxEnt model[J]. South China Fisheries Science, 2023, 19(5): 11-21. DOI: 10.12131/20230011

基于MaxEnt模型模拟中西太平洋鲣自由鱼群栖息地的研究

基金项目: 国家自然科学基金青年科学基金项目 (41506151);国家重点研发计划 “蓝色粮仓科技创新”专项 (2019YFD0901404);国家远洋渔业工程技术研究中心开放基金 (A1-2006-23-200204)
详细信息
    作者简介:

    汪伟松 (1996—),男,硕士研究生,研究方向为渔业海洋学。E-mail: weisong19961219@163.com

    通讯作者:

    李玉伟 (1984—),男,讲师,博士,研究方向为渔具力学。E-mail: ywli@shou.edu.cn

  • 中图分类号: S 932.4

Modeling habitat of skipjack tuna of free swimming school in Western and Central Pacific Ocean based on MaxEnt model

  • 摘要:

    由于漂流人工集鱼装置 (Fish aggregating device, FAD) 的大量使用对金枪鱼种群带来的负面效应,金枪鱼围网渔业转向捕捞自由群成为发展趋势,因此开展鲣 (Katsuwonus pelamis) 自由群的栖息地利用研究非常必要。使用2016—2020年中西太平洋渔业委员会 (Western and Central Pacific Fisheries Commission, WCPFC) 统计的月度金枪鱼渔业数据和不同层的水温 (SST、Temp200)、海表盐度 (SSS)、溶解氧浓度 (DO0、DO50、DO200)、东西向海流速度 (EV)、南北向海流速度 (NV)、混合层深度 (MLD)、叶绿素 a 浓度 (CHL0、CHL50、CHL100、CHL200) 共13个环境变量,通过最大熵模型 (Maximum Entropy Model, MaxEnt) 模拟鲣自由群的栖息地分布及其月变化规律。结果表明:模型测试和训练集的AUC值及灵敏度值均大于0.90,真实技巧统计值大于0.80,模型具有很强的预测能力,可用于鲣的栖息地适宜性模拟;SST和DO200是影响鲣自由群栖息地偏好的关键因子,最适范围分别为30~31  ℃、114~153 mmol·m−3。研究期内,鲣自由群高适宜栖息地主要靠近巴布亚新几内亚和所罗门群岛海域,不同时期向东延伸的范围有较大变化,经度差达到6°。研究结果可为中国金枪鱼围网船队进行鲣自由群中心渔场的预报提供参考。

    Abstract:

    Due to the negative effects of extensive use of drifting artificial fish aggregating devices (FADs) on tuna stocks, tuna purse seine fishing has become a development trend towards catching free swimming school, so it is necessary to specifically study the habitat use of free swimming school of skipjack tuna (Katsuwonus pelamis). In this study, we used monthly tuna fishery data from the Western and Central Pacific Fisheries Commission (WCPFC) from 2016 to 2020, different layers of water temperature (SST, Temp200), sea surface salinity (SSS), dissolved oxygen concentration (DO0, DO50, DO200), east-west current velocity (EV), north-south current velocity (NV), mixed layer depth (MLD), chlorophyll a concentration (CHL0, CHL50, CHL100, CHL200), and a total of 13 environmental variables by Maximum Entropy (MaxEnt) model to simulate the habitat distribution of the free swimming school and their monthly variation patterns. The results show that the AUC and sensitivity values of both the test and training set of the model were greater than 0.90, and the true skill statistics values were greater than 0.80, indicating that the model has strong predictive ability and can be used for the habitat suitability modeling of skipjack tuna. SST and DO200 were the key factors affecting the habitat preference of free swimming school, with the optimal ranges of 30−31  ℃ and 114−153 mmol·m−3, respectively. During the survey period, the highly suitable habitat for free swimming school was mainly near the waters of Papua New Guinea and Solomon Islands, with a large variation in the range extending eastward in different periods, and the difference in longitude reached 6°. The results provide references for the prediction of the central fishing ground of free swimming school of skipjack tuna by Chinese tuna purse seine fleet.

  • 根据中西太平洋渔业委员会 (Western and Central Pacific Fisheries Commission, WCPFC) 2021年金枪鱼渔业年鉴统计的主要金枪鱼物种 [鲣 (Katsuwonus pelamis)、黄鳍金枪鱼 (Thunnus albacares)、大眼金枪鱼 (T. obesus)和长鳍金枪鱼 (T. alalunga)] 的全球渔获量数据表明,鲣在全球总产量中约占56%,是金枪鱼围网渔业最为重要的捕捞对象[1]。金枪鱼类所捕食的饵料生物分散,为寻找位置动态变化的饵料需要快速迁徙游动或聚集成群,根据集群特点通常将金枪鱼群分为随附群 (Associated school),如流木群 (Log school) 及人工集鱼装置 (Fish aggregating device, FAD) 随附群,和非随附群如自由群 (Free swimming school)[2]。其中FAD随附群较为稳定,投网成功率较高,被围网渔业广泛使用,但FAD的大量使用,导致海表的物理环境被人为改变,兼捕过多的金枪鱼类幼鱼,会对目标鱼种的种群健康和海洋中上层生态系统造成严重的负面影响[3]。为此WCPFC于2008年开始实施禁渔期措施并要求购买作业天数,以期引导围网渔业转向自由群,金枪鱼围网渔业转向金枪鱼自由鱼群捕捞成为发展趋势[4]

    目前,已有不同模型应用于鲣栖息地环境偏好及栖息地的选择研究,如Mugo等[5]使用广义加性模型 (Generalized Additive Model) 预测了北太平洋西部鲣的栖息地分布状况;Dueri等[6]使用APECOSM-E数值模型 (Apex Predator Ecosystem Model-estimation) 评估了三大洋气候变化对鲣生理、空间分布和丰度的影响;Adam和Sibert[7]使用平流扩散反应模型 (Advection Diffusion Reaction Model) 分析了马尔代夫海域鲣栖息地的时空变化;Zhou等[8]开发了综合栖息地指数评估了中西太平洋鲣的栖息地质量和气候的关系;Mugo和Saitoh[9]使用支持向量机、随机森林、最大熵模型 (Maximum Entropy Model, MaxEnt) 等多种机器学习算法对比分析了西北太平洋鲣的栖息地分布。其中,MaxEnt是一种基于仅出现数据预测未知分布的机器学习方法,不需要以丰度数据作为物种信息,已广泛应用于多种海洋鱼类、鲸豚类、底栖生物等栖息地的分布模拟中[10-11]

    目前,对于鲣栖息地利用的研究大多局限于海表环境因子,如海表温度、叶绿素a浓度、海表盐度,很少将垂直环境因子纳入栖息地建模考虑,这可能会影响模型的预测性能[12]。另一方面,金枪鱼围网对自由鱼群的瞄准捕捞具有很高的空网率,大量围捕失败产生的零值导致以单网次渔获量表征的鱼群丰度被大幅低估[13],所以以鱼群出现作为自由群的物种发生信息,可能更适用于金枪鱼围网渔业数据的栖息地建模。因此,本研究以中西太平洋海域鲣的自由鱼群为研究对象,尝试利用MaxEnt确定影响其栖息地适宜性的环境因子和分布特征的月变化规律,旨在为中国金枪鱼围网船队准确地预测该海域鲣自由群的中心渔场提供参考,并提高捕捞效率,降低生产成本。

    鲣自由鱼群的渔业统计数据来自WCPFC,研究范围为120°E—140°W、35°N—35°S的中西太平洋海域,围网作业区域覆盖除4个“口袋型”公海禁渔区以外的所有海域[14] (图1),研究时段为2016—2020年。渔业生产数据包括以月为时间分辨率、1°×1°为空间分辨率的作业时间、作业位置、分鱼种和鱼群类型的产量统计等信息。本研究共提取鲣自由群的空间出现信息8 645条,并按月输入MaxEnt中建模分析。

    图  1  2016—2020 年中西太平洋金枪鱼围网捕捞自由鱼群的作业位置分布
    Figure  1.  Distribution of operating positions for tuna purse seine fishing of free fish schools in Western and Central Pacific Ocean from 2016 to 2020

    鲣的主要栖息深度为海表到200 m水层[15],其活动受各水层温度、盐度、溶解氧、叶绿素a浓度、海流速度、混合层深度等多种海洋环境因子的影响[16-18]。因此,本研究选取0、50、100、150、200 m水层的温度、溶解氧、叶绿素a浓度、盐度、东西及南北向的海流速度、混合层深度等环境因子作为建模的环境变量[19-22]。所有环境数据均下载于哥白尼海洋环境监测服务中心 (https://resources.marine.copernicus.eu),时间分辨率为月,空间分辨率均为0.25°×0.25°。

    使用MAXENT 3.4.4 (http://biodiversityinform-atics.amnh.org/open source/maxent/) 软件建模。输入层包括各月鲣自由鱼群的物种分布数据,利用ArcGIS 10.2软件将各月环境数据转换成ASC II栅格数据后输入环境图层。运行前随机选择70%的分布数据作为构建模型的训练数据,其余30%作为验证模型精度的测试数据,循环该过程20次,消除可能存在的随机性,其余参数均使用默认设置。

    由于环境因子之间存在自相关和多重共线性问题[23-24],采用方差膨胀因子 (Variance inflation factor, VIF) 对0、50、100、150和200 m水层的海水温度、溶解氧、叶绿素a浓度、海水盐度以及东西及南北向海流速度、混合层深度等环境因子进行筛选。VIF<3时的环境变量作为建模变量[25-27],其中海表温度 (SST)、200 m水温 (Temp200)、海表盐度 (SSS)、混合层深度 (MLD)、东西向的海流速度 (EV)、南北向的海流速度 (NV)、表层溶解氧浓度 (DO0)、50 m溶解氧浓度 (DO50)、200 m溶解氧浓度 (DO200)、表层叶绿素a浓度 (CHL0)、50 m叶绿素a浓度 (CHL50)、100 m叶绿素a浓度 (CHL100)、200 m叶绿素a浓度 (CHL200) 等13个环境变量被保留建模。

    分别以测试集AUC、训练集AUC、灵敏度、真实技巧统计 (True skill statistics, TSS) 4个指标共同评估模型的性能。AUC指受试者工作特征曲线 (Receiver operating characteristic curve, ROC) 下的面积,取值范围介于0~1[28];一般认为AUC值小于0.7时模型不具备良好的预测能力,0.7~0.9表示预测效果良好,大于0.9表示模型性能表现优异,预测存在与否时具有很高的准确性[29]。灵敏度是模型正确预测出现点的概率,无固定判定标准,其值越接近于1表明模型预测能力越佳[30]。TSS是评估物种存在−不存在时的衡量指标,范围介于−1~1,其中1表示与预测完全一致,0或更小值表示模型性能不优于随机模型[31]

    利用MaxEnt计算得出的环境因子年平均贡献率并用刀切法 (Jackknife) 评估影响中西太平洋鲣自由群潜在分布的环境变量的重要性。贡献率表示在所有环境因子的共同作用下,单一因子对于整个模型的贡献[32]。刀切法反映“仅有此变量”和“没有此变量”时模型的得分。“仅有此变量”时得分越高,表明该变量有较高预测能力,对物种潜在分布贡献越大。“没有此变量”则表示除此变量外其他所有变量贡献的得分,越低表明除去此变量对物种分布的影响越大[33]

    根据2016—2020年测试集AUC显示,所有模型的AUC均值都大于0.97,最高值和最低值分别出现在2019年 (0.98) 和2017年 (0.96) (图2-a)。模型训练集AUC均值与测试集AUC均值相同,均大于0.97,最高值和最低值分别出现在2019年 (0.98) 及2018年 (0.97),2018年模型AUC值变化明显高于其他年份 (图2-b)。模型的灵敏度均值均大于0.92,最高值和最低值分别出现在2019年 (0.96) 和2016年 (0.90),2016年的变化范围显著高于其他年份 (图2-c)。模型的TSS均值均大于0.85,最高值和最低值分别出现在2019年 (0.92) 和2016年 (0.81) (图2-d)。评估模型性能的4个指标均表明该模型具有很强的预测能力,可用于模拟中西太平洋鲣自由群栖息地分布。

    图  2  2016—2020 年模型性能的综合评价
    Figure  2.  Comprehensive evaluation of model performance from 2016 to 2020

    从各环境因子年平均贡献率可知 (表1),SST年平均贡献率最高,均超过55%,CHL50次之,均在10%以上,DO200年平均贡献率均超过5%,其他环境因子如Temp200、SSS的年平均贡献率很少超过5%。SST、CHL50、DO200 3个因子对模型的累计贡献率超过75%。根据Jackknife检验不同环境变量的正则化测试增益也支持这一结果 (图3)。“仅有此变量”的得分前两名分别为SST (2.13) 和DO200 (1.50),明显高于其他变量;随后较为重要的变量是浅水层的DO和CHL。综合环境因子年平均贡献率和Jackknife检验结果,得出影响中西太平洋鲣自由群栖息地偏好最关键的环境因子为SST和DO200

    表  1  2016—2020 年环境因子年平均贡献率
    Table  1.  Average annual contribution rate of environmental factors from 2016 to 2020
    环境因子
    Environmental factor
    年平均贡献率 Average annual contribution rate/%
    20162017201820192020均值 Mean
    海表温度 SST 59.35 61.11 68.61 63.15 55.47 61.54
    50 m 叶绿素 a 浓度 CHL50 11.41 10.59 10.94 12.07 13.86 11.77
    200 m 溶解氧浓度 DO200 6.14 9.90 7.12 5.50 5.38 6.81
    200 m 水温 Temp200 4.73 3.79 3.80 5.45 7.44 5.04
    表层叶绿素 a 浓度 CHL0 2.90 3.81 1.18 4.25 6.48 3.72
    海表盐度 SSS 4.77 3.75 3.06 2.22 3.61 3.48
    海流速度 EV 2.68 1.13 1.01 1.03 1.69 1.51
    表层溶解氧浓度 DO0 1.77 1.84 1.13 1.09 1.63 1.49
    混合层深度 MLD 1.77 1.18 1.06 2.35 0.87 1.45
    200 m 叶绿素 a 浓度 CHL200 1.24 1.79 0.92 1.41 1.69 1.41
    100 m 叶绿素 a 浓度 CHL100 2.09 0.39 0.62 0.72 0.60 0.88
    50 m 溶解氧浓度 DO50 0.84 0.39 0.43 0.46 0.96 0.62
    南北向的海流速度 NV 0.32 0.33 0.11 0.32 0.34 0.28
    下载: 导出CSV 
    | 显示表格
    图  3  2016—2020 年使用 Jackknife 检验得到的测试增益
    Figure  3.  Test gain obtained using Jackknife test from 2016 to 2020

    图4图5分别展示了研究期内中西太平洋海域鲣各月栖息地适宜性指数 (Habitat Suitability Index, HSI) 随SST和DO200的变化趋势。在SST低于28  ℃的海域,鲣自由群出现的概率很低,而后随着水温的升高HSI会出现快速上升,一般在30~31  ℃达到顶峰,然后再次下降,但不同月的情况具有较大差异 (图4)。在DO200高于50 mmol·m−3的海域,自由群才有可能出现,其分布概率均小于0.2,此时的环境不适宜自由群活动,随着DO浓度增加,HSI上升明显,DO200介于114~153 mmol·m−3是鲣最适宜的栖息环境,随后HSI快速下降,一般超过200 mmol·m−3时就不适合鲣生存,不同月份间的阈值差异不显著 (图5)。

    图  4  海表温度对中西太平洋鲣栖息地适宜性的响应曲线
    Figure  4.  Response curves of SST to habitat suitability of skipjack tuna in Western and Central Pacific Ocean
    图  5  200 m 水层溶解氧浓度对中西太平洋鲣栖息地适宜性的响应曲线
    Figure  5.  Response curves of DO200 to habitat suitability of skipjack tuna in Western and Central Pacific Ocean

    图6展示了2016—2020年鲣自由群栖息地适宜性分布的变动情况。在研究期内,自由群的高适宜栖息地均位于赤道海域附近 (8°N—10°S),核心区域从巴布亚新几内亚和所罗门群岛向东延伸,但是不同月份间延伸的位置有较大变化,从而导致高适宜性区域面积发生消长。高适宜区域向东延伸最远的月份为9月,位于171°W,最小月份为6月,位于177°W,两者经度差可达到6°。各月高适宜栖息地面积会发生一定的改变,向东高适宜性区域面积远小于核心区域,最大在9月,最小在12月,各月栖息地范围变动总体呈现先向北迁移而后向东再向北迁移的趋势。

    图  6  2016—2020 年中西太平洋鲣自由鱼群潜在栖息地适宜性指数分布
    注:图中白色区域表示无数据。
    Figure  6.  Distribution of potential habitat suitability index for free swimming school of skipjack tuna in Western and Central Pacific Ocean from 2016 to 2020
    Note: The white area in the figure indicates no data.

    模型准确性评估一般采用阈值无关指标或阈值相关指标,AUC是一种常见的阈值无关指标,用于评估真实存在与不存在数据的匹配程度,其中测试集与训练集的AUC值在模型预测能力评估中可作为良好的指标[34]。灵敏度及TSS则作为一种阈值相关指标,用于模型预测物种出现点的概率,其值越接近于1说明模型预测能力越佳[35]。通常,会选择结合多种指标评估栖息地模型的综合性能,譬如West等[36]在探究入侵物种栖息地适宜性时,使用测试集AUC和灵敏度评估了MaxEnt模型的性能,认为AUC大于0.9时模型具有高预测性能,灵敏度越大表示模型预测结果越接近于真实情况。本研究使用MaxEnt计算得出的AUC、灵敏度值均大于0.9,TSS值均大于0.8,表明模型预测结果有很高的可信度,具有较强的预测能力,可用于中西太平洋鲣自由群栖息地适宜性评价的研究。

    事实上,MaxEnt自2004年提出之后迅速在国内外得到广泛应用,在生态学领域发展迅速[37-38]。MaxEnt根据物种分布与环境因子的不完全信息,探寻熵最大 (即随机变量不确定性的度量) 的随机分配,预测并估算该物种的潜在栖息地分布,即使用不完整的分布点数据也可实现准确预测,是一种基于“当前存在”预测“未知分布”的机器学习方法[39]。目前,MaxEnt是模拟海洋物种栖息地研究中使用频率最高的方法,在鸟类、鱼类、植物等物种栖息地研究中应用广泛[38, 40-43]。在金枪鱼栖息地的研究中,Mugo和Saitoh[9]使用MaxEnt对西北太平洋鲣的适宜栖息地分布进行模拟,Nursan等[41]应用MaxEnt对印度尼西亚海域鲣潜在栖息地进行建模分析;结果均表明MaxEnt在探究热带金枪鱼栖息地分布与环境因子关系的研究中表现优异。

    鲣自由群的渔场分布特征会随月份产生变化,已有相关研究对此做了阐述。陈世泳[43]使用1994—2004年台湾省金枪鱼围网渔业统计资料研究中西太平洋鲣不同群性 (自由群、天然流木群、人工集鱼装置群) 渔场重心的月际变动,并结合南方涛动指数 (Southern Oscillation Index, SOI) 及SST进行分析;发现自由群在厄尔尼诺事件 (El Niño) 发生时,表征中西太平洋暖池东边界的29  ℃等温线向东太平洋推移,促使鲣向东迁移,在拉尼娜 (La Niña) 时期则向西移动,El Niño及La Niña交替出现,造成海水温度异常的周期循环而使鲣向不同方向迁徙。在时间上,金枪鱼围网渔场重心会沿赤道作季节性移动,6月以前集中于140°E—165°E,之后向东移动至165°E—170°E,9—12月逐渐向西返回,在纬度上渔场重心集中于5°N—5°S,纬度重心在各月间并无明显差异。Zhou等[44]使用2013—2017年的中国金枪鱼围网船队统计数据并结合SST和DO等海洋环境因子研究了中西太平洋鲣自由群的栖息地偏好,在其研究期内也发生了El Niño及La Niña现象,发现栖息地的适宜性分布会沿着暖池位移方向做季节性移动,进而带动鲣自由群进行迁徙,其中La Niña事件对鲣丰度及栖息地的影响更为显著。

    本研究的区域范围相比上述研究更加广泛,研究期内2016、2018和2019年发生的El Niño事件及2017和2020年发生的La Niña事件均能促使暖池位置移动,进而导致自由群的高适宜栖息地在东西方向上发生位移。另一方面,暖池东边界的位置不仅在东西方向变化,也会随月份自南向北发生转变[45]。由图6可知,1—4月鲣自由群栖息地主要位于6°N—11°S,自5月栖息地开始自南向北移动,至9月在南北方向上栖息地范围变化不大,随后栖息地向南移动,于12月到达最南端的11°S附近海域。综上所述,可以获知中西太平洋海域ENSO循环周期会驱动鲣自由群的高适宜性栖息地范围发生相应变动,最终导致高度洄游鱼类的种群重心不断改变。另外需要注意的是,由于4个“口袋型”公海内禁止金枪鱼围网作业 (图1),所以无法对鲣自由群在这些区域的真实分布进行采样,目前整个研究区域的潜在适宜区估计来自于采样区域调查结果的外推预测 (图6),其可靠性需要今后通过使用独立于渔业的数据进行进一步验证。

    鲣的自由群和漂流物随附群具有完全不同的集群机制,对于栖息地的选择和环境偏好也可能具有很大差异,本文专门针对自由群进行研究。本研究发现,SST、DO200是决定中西太平洋鲣自由群栖息地适宜性的关键环境因子;其中SST是决定鲣分布最重要的因素,可作为寻找鲣渔场的水温指标[46]。已有大量研究证实鲣的分布会受SST的强烈影响,SST被认为是解释温暖水域鲣出现及丰度变化的主因[47]。在适宜范围方面,已往研究认为中西太平洋金枪鱼围网捕获的金枪鱼大多分布于水温28  ℃以上水域,其中鲣的主要栖息水温为28~31  ℃[48]。在本研究中,当SST低于28  ℃时,2016—2020年各月鲣的出现概率普遍低于0.2,表明鲣在SST 28  ℃以下的海域难以栖息。高适宜性栖息地的SST范围为30~31  ℃,但很多月份在SST高于32  ℃的情况下,鲣的分布概率仍然很高 (图4),这说明在气候变化背景下不同时间情景的适宜性偏好或环境耐受力会有所区别,这种现象在不区分月际模型时难以探测。

    相对于SST,目前关于DO浓度对自由群栖息地选择影响的研究仍较少,特别是考虑次表层DO的研究。鲣作为一种高度洄游性鱼类,在游泳时肌肉中的蛋白质转换需要消耗大量氧气 (O2)[49]。Bushnell和Brill[50]认为DO浓度会制约鲣的游泳速度,从而影响其分布。鲣一天中有60%~80%的时间在温跃层 (200 m) 以上活动[51],生理和觅食需求是改变鲣空间分布的关键因素,并有助于形成鲣的高适宜栖息地[52]。标记放流实验则表明鲣虽然主要在温跃层活动,但由于身体的热惯性增强,偶尔会短暂地潜入温跃层以下,并对一定深度的DO具有高耐受性[53]。因此,在海洋的中上层生境中,DO浓度在不同时空尺度上制约着鲣在垂直方向上的分布[54]。本研究发现200 m水层的DO浓度对于中西太平洋鲣自由群的分布相当重要,该指标过低 (小于50 mmol·m−3) 或过高 (大于200 mmol·m−3) 均会制约鲣的分布,其最适范围介于114~153 mmol·m−3。但需要指出的是,区别于SST适宜性的叠加曲线,不同月份间DO200曲线的变化趋势十分接近,这说明次表层DO对中西太平洋鲣的分布具有稳定的约束作用 (图5)。

    本研究利用MaxEnt探究了中西太平洋鲣自由群栖息地适宜性的时空分布规律,和对其分布具有重要作用的海洋环境因子。目前,本研究仅根据相关研究选取了常见的环境因子进行建模分析,建议在未来纳入更多具有潜在影响的环境变量进行建模,如温跃层深度[55]、涡动能[56]等。另外,本研究仅揭示了鲣在自由栖息状态下的栖息地环境偏好和适宜性格局,但由于中西太平洋海域FAD的大量使用,人为地改变了大洋表层的栖息环境,会对鲣的栖息地选择造成显著干扰,这需要在今后的研究中进行进一步分析。

  • 图  1   2016—2020 年中西太平洋金枪鱼围网捕捞自由鱼群的作业位置分布

    Figure  1.   Distribution of operating positions for tuna purse seine fishing of free fish schools in Western and Central Pacific Ocean from 2016 to 2020

    图  2   2016—2020 年模型性能的综合评价

    Figure  2.   Comprehensive evaluation of model performance from 2016 to 2020

    图  3   2016—2020 年使用 Jackknife 检验得到的测试增益

    Figure  3.   Test gain obtained using Jackknife test from 2016 to 2020

    图  4   海表温度对中西太平洋鲣栖息地适宜性的响应曲线

    Figure  4.   Response curves of SST to habitat suitability of skipjack tuna in Western and Central Pacific Ocean

    图  5   200 m 水层溶解氧浓度对中西太平洋鲣栖息地适宜性的响应曲线

    Figure  5.   Response curves of DO200 to habitat suitability of skipjack tuna in Western and Central Pacific Ocean

    图  6   2016—2020 年中西太平洋鲣自由鱼群潜在栖息地适宜性指数分布

    注:图中白色区域表示无数据。

    Figure  6.   Distribution of potential habitat suitability index for free swimming school of skipjack tuna in Western and Central Pacific Ocean from 2016 to 2020

    Note: The white area in the figure indicates no data.

    表  1   2016—2020 年环境因子年平均贡献率

    Table  1   Average annual contribution rate of environmental factors from 2016 to 2020

    环境因子
    Environmental factor
    年平均贡献率 Average annual contribution rate/%
    20162017201820192020均值 Mean
    海表温度 SST 59.35 61.11 68.61 63.15 55.47 61.54
    50 m 叶绿素 a 浓度 CHL50 11.41 10.59 10.94 12.07 13.86 11.77
    200 m 溶解氧浓度 DO200 6.14 9.90 7.12 5.50 5.38 6.81
    200 m 水温 Temp200 4.73 3.79 3.80 5.45 7.44 5.04
    表层叶绿素 a 浓度 CHL0 2.90 3.81 1.18 4.25 6.48 3.72
    海表盐度 SSS 4.77 3.75 3.06 2.22 3.61 3.48
    海流速度 EV 2.68 1.13 1.01 1.03 1.69 1.51
    表层溶解氧浓度 DO0 1.77 1.84 1.13 1.09 1.63 1.49
    混合层深度 MLD 1.77 1.18 1.06 2.35 0.87 1.45
    200 m 叶绿素 a 浓度 CHL200 1.24 1.79 0.92 1.41 1.69 1.41
    100 m 叶绿素 a 浓度 CHL100 2.09 0.39 0.62 0.72 0.60 0.88
    50 m 溶解氧浓度 DO50 0.84 0.39 0.43 0.46 0.96 0.62
    南北向的海流速度 NV 0.32 0.33 0.11 0.32 0.34 0.28
    下载: 导出CSV
  • [1]

    Food and Agriculture Organization of the United Nations. The State of World Fisheries and Aquaculture, 2021[R]. Rome: FAO, 2021: 127-148.

    [2]

    TRYGONIS V, GEORGAKARAKOS S, DAGORN L, et al. Spatiotemporal distribution of fish schools around drifting fish aggregating devices[J]. Fish Res, 2016, 177: 39-49. doi: 10.1016/j.fishres.2016.01.013

    [3]

    HALLIER J P, GAERTNER D. Drifting fish aggregation devices could act as an ecological trap for tropical tuna species[J]. Mar Ecol Prog Ser, 2008, 353: 255-264. doi: 10.3354/meps07180

    [4]

    TSENG C T, SUN C L, YEH S Z, et al. Spatio-temporal distributions of tuna species and potential habitats in the Western and Central Pacific Ocean derived from multi-satellite data[J]. Int J Remote Sens, 2010, 31(17/18): 4543-4558.

    [5]

    MUGO R, SAITOH S E I I, NIHIRA A, et al. Habitat characteristics of skipjack tuna (Katsuwonus pelamis) in the western North Pacific: a remote sensing perspective[J]. Fish Oceanogr, 2010, 19(5): 382-396. doi: 10.1111/j.1365-2419.2010.00552.x

    [6]

    DUERI S, BOPP L, MAURY O. Projecting the impacts of climate change on skipjack tuna abundance and spatial distribution[J]. Glob Chang Biol, 2014, 20(3): 742-753. doi: 10.1111/gcb.12460

    [7]

    ADAM M S, SIBERT J R. Population dynamics and movements of skipjack tuna (Katsuwonus pelamis) in the Maldivian fishery: analysis of tagging data from an advection-diffusion-reaction model[J]. Aquat Living Resour, 2002, 15(1): 13-23. doi: 10.1016/S0990-7440(02)01155-5

    [8]

    SALAZAR J E, BENAVIDES I F, CABRERA C V P, et al. Generalized additive models with delayed effects and spatial autocorrelation patterns to improve the spatiotemporal prediction of the skipjack (Katsuwonus pelamis) distribution in the Colombian Pacific Ocean[J]. Reg Stud Mar Sci, 2021, 45: 101-115.

    [9]

    MUGO R, SAITOH S I. Ensemble modelling of skipjack tuna (Katsuwonus pelamis) habitats in the western north Pacific using satellite remotely sensed data: a comparative analysis using machine-learning models[J]. Remote Sens, 2020, 12(16): 2591-2605. doi: 10.3390/rs12162591

    [10]

    ALABIA I D, SAITOH S I, MUGO R, et al. Seasonal potential fishing ground prediction of neon flying squid (Ommastrephes pastrami) in the western and central North Pacific[J]. Fish Oceanogr, 2015, 24(2): 190-203. doi: 10.1111/fog.12102

    [11]

    CHEN B Y, HONG Z, HAO X Q, et al. Environmental models for predicting habitat of the Indo-Pacific humpback dolphins in Fujian, China[J]. Aquat Conserv, 2020, 30(4): 787-793. doi: 10.1002/aqc.3279

    [12]

    YOSHINAGA D H, FRANK H A. Histamine-producing bacteria in decomposing skipjack tuna (Katsuwonus pelamis)[J]. Appl Environ Microb, 1982, 44(2): 447-452. doi: 10.1128/aem.44.2.447-452.1982

    [13]

    DAI L, WANG X, STAPLES K W, et al. Factors influencing successful fishing of tuna free-swimming schools in the equatorial western Pacific Ocean[J]. Turk J Fish Aquat Sc, 2019, 20(5): 341-350.

    [14]

    SIBERT J, SENINA I, LEHODEY P, et al. Shifting from marine reserves to maritime zoning for conservation of Pacific bigeye tuna (Thunnus obesus)[J]. P Natl Acad Sci USA, 2012, 109(44): 18221-18225. doi: 10.1073/pnas.1209468109

    [15]

    PRINCE E D, GOODYEAR C P. Hypoxia-based habitat compression of tropical pelagic fishes[J]. Fish Oceanogr, 2006, 15(6): 451-464. doi: 10.1111/j.1365-2419.2005.00393.x

    [16]

    LEHODEY P, BERTIGNAC M, HAMPTON J, et al. El Niño Southern Oscillation and tuna in the western Pacific[J]. Nature, 1997, 389(6652): 715-718. doi: 10.1038/39575

    [17]

    DI Y U, CHANG F C, BIN W, et al. Characterization of acid-and pepsin-soluble collagens from spines and skulls of skipjack tuna (Katsuwonus pelamis)[J]. Chin J Nat Medicines, 2014, 12(9): 712-720. doi: 10.1016/S1875-5364(14)60110-2

    [18]

    ROGER C. The plankton of the tropical western Indian Ocean as a biomass indirectly supporting surface tunas (yellowfin, Thunnus albacares and skipjack, Katsuwonus pelamis)[J]. Environ Biol Fish, 1994, 39: 161-172. doi: 10.1007/BF00004934

    [19]

    LEROY B, ITANO D G, USU T, et al. Vertical behavior and the observation of FAD effects on tropical tuna in the warm-pool of the western Pacific Ocean[M]//NIELSEN J L, ARRIZABALAGA H, FRAGOSO N, et al. Tagging and tracking of marine animals with electronic devices. Reviews: methods and technologies in fish biology and fisheries, vol 9. Dordrecht: Springer, 2009: 161-179.

    [20]

    KIYOFUJI H, AOKI Y, KINOSHITA J, et al. Northward migration dynamics of skipjack tuna (Katsuwonus pelamis) associated with the lower thermal limit in the western Pacific Ocean[J]. Prog Oceanogr, 2019, 175: 55-67. doi: 10.1016/j.pocean.2019.03.006

    [21]

    SCHAEFER K M. Vertical movement pattern of skipjack tuna (Katsuwonus pelamis) in the eastern equatorial Pacific Ocean, as revealed with archival tags[J]. Fish Bull, 2007, 105: 379-389.

    [22] 杨胜龙, 周甦芳, 周为峰, 等. 基于Argo数据的中西太平洋鲣渔获量与水温、表层盐度关系的初步研究[J]. 大连水产学院学报, 2010, 25(1): 34-40.
    [23]

    GARCIA C B, GARCIA J, LOPEZ MARTIN M M, et al. Collinearity: revisiting the variance inflation factor in ridge regression[J]. J Appl Stat, 2015, 42(3): 648-661. doi: 10.1080/02664763.2014.980789

    [24]

    IGARASHI H, SAITOH S I, ISHIKAWA Y, et al. Identifying potential habitat distribution of the neon flying squid (Ommastrephes bartramii) off the eastern coast of Japan in winter[J]. Fish Oceanogr, 2018, 27(1): 16-27. doi: 10.1111/fog.12230

    [25] 赵静, 柳晓雪, 吴建辉, 等. 零膨胀模型在珍稀鱼类资源时空分布预测中的应用: 以长江口刀鲚为例[J]. 生态学杂志, 2020, 39(9): 3155-3163. doi: 10.13292/j.1000-4890.202009.028
    [26]

    Sagarese S R, Frisk M G, Cerrato R M, et al. Application of generalized additive models to examine ontogenetic and seasonal distributions of spiny dogfish (Squalus acanthias) in the Northeast (US) shelf large marine ecosystem[J]. Can J Fish Aquat Sci, 2014, 71(6): 847-877. doi: 10.1139/cjfas-2013-0342

    [27] 马金, 黄金玲, 陈锦辉, 等. 基于GAM的长江口鱼类资源时空分布及影响因素[J]. 水产学报, 2020, 44(6): 936-946.
    [28]

    PHILLIPS S J, ANDERSON R P, SCHAPIRE R E. Maximum entropy modeling of species geographic distributions[J]. Ecol Model, 2006, 190(3/4): 231-259.

    [29]

    MANEL S, WILLIAMS H C, ORMEROD S J. Evaluating presence absence models in ecology: the need to account for prevalence[J]. J Appl Ecol, 2001, 38(5): 921-931. doi: 10.1046/j.1365-2664.2001.00647.x

    [30]

    ZU T P, KANG R, WEN M L, et al. Belief reliability distribution based on maximum entropy principle[J]. IEEE Access, 2017, 6: 1577-1582.

    [31]

    ALLOUCHE O, TSOAR A, KADMON R. Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic (TSS)[J]. J Appl Ecol, 2006, 43(6): 1223-1232. doi: 10.1111/j.1365-2664.2006.01214.x

    [32]

    MEYSMAN F J R, BRUERS S. Ecosystem functioning and maximum entropy production: a quantitative test of hypotheses[J]. Philosophical Transactions of the Royal Society B:Biol Sci, 2010, 365(1545): 1405-1416. doi: 10.1098/rstb.2009.0300

    [33]

    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

    [34]

    ELITH J H, GRAHAM C P, ANDERSON R, et al. Novel methods improve prediction of species distributions from occurrence data[J]. Ecography, 2006, 29(2): 129-151. doi: 10.1111/j.2006.0906-7590.04596.x

    [35]

    ZU T, KANG R, WEN M, et al. Belief reliability distribution based on maximum entropy principle[J]. IEEE Access, 2017, 6: 1577-1582.

    [36]

    WEST A M, KUMAR S, BROWN C S, et al. Field validation of an invasive species MaxEnt model[J]. Ecol Inform, 2016, 36: 126-134. doi: 10.1016/j.ecoinf.2016.11.001

    [37]

    KUMAR P. Assessment of impact of climate change on Rhododendrons in Sikkim Himalayas using Maxent modelling: limitations and challenges[J]. Biodivers Conserv, 2012, 21(5): 1251-1266. doi: 10.1007/s10531-012-0279-1

    [38]

    MORENO R, ZAMORA R, MOLINA J R, et al. Predictive modeling of microhabitats for endemic birds in South Chilean temperate forests using Maximum Entropy (MaxEnt)[J]. Ecol Inform, 2011, 6(6): 364-370. doi: 10.1016/j.ecoinf.2011.07.003

    [39]

    LIU S J, YANG J. Modeling spatial patterns of forest fire in Heilongjiang Province using Generalized Linear Model and Maximum Entropy Model[J]. Chin J Ecol, 2013, 32(6): 1620-1628.

    [40] 梁阳阳, 陈康, 崔凯, 等. 气候变化情景下须鳗鰕虎鱼在中国的潜在地理分布[J]. 大连海洋大学学报, 2022, 37(5): 739-746.
    [41]

    NURSAN M, YONVITNER Y, AGUS S B. Distribution of skipjack (Katsuwonus pelamis) fishing areas using purse seine fishing equipment in WPP 573[J]. J Trop Fish Manag, 2022, 6(1): 11-20.

    [42]

    ZHANG K L, YAO L J, MENG J S, et al. Maxent modeling for predicting the potential geographical distribution of two peony species under climate change[J]. Sci Total Environ, 2018, 634: 1326-1334. doi: 10.1016/j.scitotenv.2018.04.112

    [43] 陈世泳. 中西太平洋正鲣渔场推移与表明水温变化之关系[D]. 基隆: 台湾海洋大学, 2006: 24-26.
    [44]

    ZHOU C, HU Y Y, CAO J, et al. Comparison of nominal and standardized catch per unit effort data in quantifying habitat suitability of skipjack tuna in the equatorial Pacific Ocean[J]. Acta Oceanol Sin, 2022, 41(3): 1-10. doi: 10.1007/s13131-021-1922-z

    [45]

    PICAUT J, IOUALALEN M, MENKÈS C, et al. Mechanism of the zonal displacements of the Pacific warm pool: implications for ENSO[J]. Science, 1996, 274(5292): 1486-1489. doi: 10.1126/science.274.5292.1486

    [46]

    ELY B, VIÑAS J, ALVARADO BREMER J R, et al. Consequences of the historical demography on the global population structure of two highly migratory cosmopolitan marine fishes: the yellowfin tuna (Thunnus albacares) and the skipjack tuna (Katsuwonus pelamis)[J]. BMC Evol Biol, 2005, 5(1): 1-9. doi: 10.1186/1471-2148-5-1

    [47]

    LEHODEY P. The pelagic ecosystem of the tropical Pacific Ocean: dynamic spatial modelling and biological consequences of ENSO[J]. Prog Oceanogr, 2001, 49(1/2/3/4): 439-468.

    [48]

    BOYCE D G, TITTENSOR D P, WORM B. Effects of temperature on global patterns of tuna and billfish richness[J]. Mar Ecol Prog Ser, 2008, 355: 267-276. doi: 10.3354/meps07237

    [49]

    MATHIEU-COSTELLO O, BRILL R W, HOCHACHKA P W. Structural basis for oxygen delivery: muscle capillaries and manifolds in tuna red muscle[J]. Comp Biochem Physiol A, 1996, 113(1): 25-31. doi: 10.1016/0300-9629(95)02059-4

    [50]

    BUSHNELL P G, BRILL R W. Responses of swimming skipjack (Katsuwonus pelamis) and yellowfin (Thunnus albacares) tunas to acute hypoxia, and a model of their cardiorespiratory function[J]. Physiol Zool, 1991, 64(3): 787-811. doi: 10.1086/physzool.64.3.30158207

    [51]

    SCHAEFER K M, FULLER D W. Vertical movement pattern of skipjack tuna (Katsuwonus pelamis) in the eastern equatorial Pacific Ocean, as revealed with archival tags[J]. Fish Bull, 2007, 105: 379-989.

    [52]

    ZAINUDDIN M, KIYOFUJI H, SAITOH K, et al. Using multi-sensor satellite remote sensing and catch data to detect ocean hot spots for albacore (Thunnus alalunga) in the northwestern North Pacific[J]. Deep-Sea Res II, 2006, 53(3/4): 419-431. doi: 10.1016/j.dsr2.2006.01.007

    [53]

    BERNAL D, BRILL R W, DICKSON K A, et al. Sharing the water column: physiological mechanisms underlying species-specific habitat use in tunas[J]. Rev Fish Biol Fisher, 2017, 27(4): 843-880. doi: 10.1007/s11160-017-9497-7

    [54]

    DEARY A L, MORET-FERGUSON S, ENGELS M, et al. Influence of central Pacific oceanographic conditions on the potential vertical habitat of four tropical tuna species1[J]. Pac Sci, 2015, 69(4): 461-475. doi: 10.2984/69.4.3

    [55]

    ANDRADE H A. The relationship between the skipjack tuna (Katsuwonus pelamis) fishery and seasonal temperature variability in the south western Atlantic[J]. Fish Oceanogr, 2003, 12(1): 10-18. doi: 10.1046/j.1365-2419.2003.00220.x

    [56]

    LLOPIZ J K, HOBDAY A J. A global comparative analysis of the feeding dynamics and environmental conditions of larval tunas, mackerels, and billfishes[J]. Deep-Sea Res II, 2015, 113: 113-124. doi: 10.1016/j.dsr2.2014.05.014

图(6)  /  表(1)
计量
  • 文章访问数:  371
  • HTML全文浏览量:  176
  • PDF下载量:  89
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-02-04
  • 修回日期:  2023-05-08
  • 录用日期:  2023-05-24
  • 网络出版日期:  2023-06-04
  • 刊出日期:  2023-10-04

目录

/

返回文章
返回