基于ARIMA的海洋尼诺指数对中西太平洋黄鳍金枪鱼年际CPUE的影响

王啸, 刘文俊, 张健

王啸, 刘文俊, 张健. 基于ARIMA的海洋尼诺指数对中西太平洋黄鳍金枪鱼年际CPUE的影响[J]. 南方水产科学, 2023, 19(4): 10-20. DOI: 10.12131/20230007
引用本文: 王啸, 刘文俊, 张健. 基于ARIMA的海洋尼诺指数对中西太平洋黄鳍金枪鱼年际CPUE的影响[J]. 南方水产科学, 2023, 19(4): 10-20. DOI: 10.12131/20230007
WANG Xiao, LIU Wenjun, ZHANG Jian. Effect of Oceanic Niño index on interannual CPUE of yellowfin tuna (Thunnus albacares) in Western and Central Pacific Ocean based on ARIMA model[J]. South China Fisheries Science, 2023, 19(4): 10-20. DOI: 10.12131/20230007
Citation: WANG Xiao, LIU Wenjun, ZHANG Jian. Effect of Oceanic Niño index on interannual CPUE of yellowfin tuna (Thunnus albacares) in Western and Central Pacific Ocean based on ARIMA model[J]. South China Fisheries Science, 2023, 19(4): 10-20. DOI: 10.12131/20230007

基于ARIMA的海洋尼诺指数对中西太平洋黄鳍金枪鱼年际CPUE的影响

基金项目: 农业农村部全球重要鱼种资源动态监测评估项目;国家重点研发计划项目 (2020YFD0900803)
详细信息
    作者简介:

    王 啸 (1995—),男,硕士,研究方向为远洋渔业资源。E-mail: 457506879@qq.com

    通讯作者:

    张 健 (1979—),男,副教授,博士,研究方向为生态型渔具渔法。E-mail: j-zhang@shou.edu.cn

  • 中图分类号: S 931.3

Effect of Oceanic Niño index on interannual CPUE of yellowfin tuna (Thunnus albacares) in Western and Central Pacific Ocean based on ARIMA model

  • 摘要: 黄鳍金枪鱼 (Thunnus albacares) 为高度洄游的大洋性鱼类,有较高的生态和经济价值,中西太平洋 (Western and Central Pacific Ocean, WCPO) 是全球金枪鱼捕捞产量最高的海区。为了解和预测中西太平洋黄鳍金枪鱼不同渔业对气候变化的反应,根据1990—2020年世界各国在中西太平洋的围网和延绳钓作业以及海洋尼诺指数 (Oceanic Niño index, ONI) 数据,分析了常规自回归积分滑动平均模型 (Autoregressive Integrated Moving Average Model, ARIMA) 和加入ONI标准差为协变量的动态ARIMA模型在渔业资源量研究中的适用性,以及ONI对中西太平洋黄鳍金枪鱼年际单位捕捞努力量渔获量 (Catch per unit effort, CPUE) 的影响。结果表明:1) 常规ARIMA模型能够充分考虑中西太平洋黄鳍金枪鱼年CPUE的变化特征,可用于黄鳍金枪鱼年CPUE的长期拟合;2) 相比常规ARIMA模型,动态ARIMA模型的拟合度更好,拟合值和真实值的相关性更高,同时平均绝对误差、均方根误差更小;3) ONI对中西太平洋赤道南北海域黄鳍金枪鱼的年CPUE影响不同,相对而言,在赤道以北,ONI的影响因素更关键,模型的拟合度更高;4) ONI对中西太平洋不同渔业的黄鳍金枪鱼的年CPUE影响有差别,对中西太平洋黄鳍金枪鱼延绳钓渔业存在滞后1~2年的影响,而在强厄尔尼诺和强拉尼娜现象时,对围网渔业的影响速度较快,不存在滞后。
    Abstract: As a highly migratory pelagic fish, yellowfin tuna (Thunnus albacares) has high ecological and economic value. The Western and Central Pacific Ocean (WCPO) is the sea area with the highest tuna production of all oceans. In order to understand and predict the response of yellowfin tuna to climate change at different life stages in WCPO, we used the catch data of yellowfin tuna in purse seining and pelagic longlining and Oceanic Niño index (ONI) data from 1990 to 2020 in the WCPO to validate the applicability of general ARIMA (Autoregressive integrated moving average) model and dynamic ARIMA model, so as to explore the influence of the ONI on the interannual CPUE (Catch per unit effort) of yellowfin tuna. The results show that: 1) General ARIMA models could be used for long-term fitting of annual CPUE of yellowfin tuna in the WCPO, taking full account of the variability characteristics of annual CPUE of yellowfin tuna. 2) Compared with the general ARIMA model, the dynamic ARIMA model provided a better fit and a higher correlation between the fitted and true values, as well as smaller mean absolute and root mean square errors. 3) The influence of the ONI on the annual CPUE of yellowfin tuna differed between the northern and southern equatorial regions of the WCPO, with the ONI being a more critical factor and a better model fit relatively north of the equator. 4) The ONI had different impacts on the annual CPUE of yellowfin tuna in different fisheries in the WCPO, with a 1–2 years' lag in the ONI for the yellowfin tuna longline fishery in the WCPO, and a faster impact on the purse seine fishery during strong El Niño and strong La Niña events, without a lag.
  • 钝缀锦蛤 (Tapes dorsatus) 隶属软体动物门、瓣鳃纲、帘蛤目、帘蛤科、缀锦蛤属,俗称沙包螺,主要分布于西南太平洋热带、亚热带海区,常见于我国东海、南海和台湾海峡[1]。其主要栖息于海区潮间带中、低潮区至潮下带及浅海3~15 cm的泥沙底质中,营埋栖生活,在潮间带下部珊瑚块和碎石泥沙中亦有栖息[2-3]。钝缀锦蛤为中大型经济贝类,成体可达100 g以上,具有个体大、生长快、适应性强、品质优等特点,且味道鲜美、营养丰富。目前钝缀锦蛤养殖业主要兴起于北部湾浅海区,以底播养殖为主[4]

    国内外有关钝缀锦蛤的研究报道主要集中在生态习性[2-6]、分类学[7-8]、繁殖生物学[9-12]、育苗[13-14]等方面。对钝缀锦蛤的养殖实验仅见于不同养殖方式和底质对澳大利亚钝缀锦蛤生长与存活的影响[15]以及钝缀锦蛤在防城港海区中培期筐式养殖防风浪模式实验[16]。随着钝缀锦蛤养殖规模的逐渐扩大,传统底播养殖方式出现养殖成活率低、采捕时贝壳破碎率高、劳动量大、回捕率低等问题。探索高效且便于管理的养殖模式是解决这些问题的有效途径之一。已有研究表明,菲律宾蛤仔 (Ruditapes philippinarum)[17]、扁平钳蛤 (Isognomon ephippium)[18]等筏式吊养模式比底播养殖模式具有成活率高、生长速度快等优点。目前尚未见有关钝缀锦蛤在广西北部湾海域不同养殖密度及养殖模式的比较研究。本实验采用筐式养殖开展钝缀锦蛤中培期和养成期不同养殖密度的生长比较研究,并在防城港白龙珍珠湾首次进行网笼吊养,探索钝缀锦蛤在广西北部湾3个海区更为高效且方便管理的养殖模式和各海区的适宜养殖密度,为其在广西北部湾海域的规模化及高效养殖提供参考。

    实验所用的钝缀锦蛤稚贝于2021年3月1日购自广东省雷州市覃斗苗场,平均壳长 (2.69±0.51) mm、平均壳高 (1.10±0.18) mm、平均壳宽 (2.21±0.45) mm、平均体质量 (2.17±0.17) mg,2021年3月1日—5月8日在防城港市广西富群海水种苗繁殖有限公司室内标粗68 d后,取大小均一、壳形完整、体质健康无损伤、活力好的钝缀锦蛤贝苗进行海区筐式中培实验,中培实验结束后继续进行海区养成实验。

    海区中培和养成所用的实验筐均为聚乙烯塑料筐 (长45 cm×宽35 cm×高15 cm),并搭配孔径1 cm的塑料硬盖 (长45 cm×宽35 cm) (图1-a1-b)。塑料筐四周的上部分有3排透水口,筐下部分及底部密闭。实验时为方便筐内外进行水流、物质交换,在筐底部钻9~12个孔径为0.50 cm的圆形透水口。

    图  1  中培、养成阶段筐具和吊养网笼
    Fig. 1  Plastic baskets and hanging net cage used for nursery and adult culture stage

    网笼吊养实验使用聚乙烯胶丝编织而成的单层网笼 (底径45 cm、高40 cm、网孔3 cm×3 cm),如图1-c所示。网笼内置养殖网袋 (底径35 cm、高30 cm、网孔1 cm×1 cm),网袋内装入1 kg浮石颗粒基质 (火山石比重为0.45 g·cm−3、孔隙率75%~80%、基质粒径1.50~2.50 cm)。

    钝缀锦蛤海区中培实验于2021年5月15日开始,至2021年7月14日结束,实验周期为60 d。海区中培实验前随机取室内中培稚贝30粒测量其初始规格:平均壳长 (6.01±0.56) mm、平均壳高 (2.22±0.18) mm、平均壳宽 (4.61±0.37) mm、平均体质量 (0.04±0.01) g。中培实验分别在广西钦州市犀牛脚镇渔港海区 (108.73°E, 21.63°N)、北海市营盘镇沙虫坪海区 (109.28°E, 21.26°N)、防城港市白龙珍珠湾海区 (108.22°E, 21.52°N) 潮间带低潮区近潮下带滩涂进行,大潮潮高约4~5 m。海区中培实验设计5个密度组:d1 (1 000 粒·筐−1)、d2 (1 500 粒·筐−1)、d3 (2 000 粒·筐−1)、d4 (2 500 粒·筐−1)、d5 (3 000 粒·筐−1),每个密度设置12个平行组。养殖实验场地见图2。实验过程中将塑料筐半埋入滩涂,筐上部的四周孔洞外露,并加盖孔径1 cm的塑料硬盖。

    图  2  3 个海区钝缀锦蛤养殖实验场地
    Fig. 2  Experimental mariculture site of T. dorsatus in three sea areas

    海区中培实验期间,每15 d退潮时对每个密度的3个平行组进行不重复取样测量,每个平行组随机取30粒稚贝测量壳长、壳宽、壳高以及体质量 (精确到0.01 g),并计算成活率。测量当日实验海区的水温、滩温、盐度、pH、透明度和溶解氧 (DO),分别用水银温度计、盐度计、CT-6023高精度pH计、塞氏盘、美国YSI ProPlus便携式水质分析仪测量。每次采样测量提前进入养殖实验场地,当海水还未完全退去时在实验场地四周及中心5个点使用1 L的玻璃采水器采集滩涂底层水样,混合均匀后装于1.5 L聚乙烯塑料瓶中,藻样采集后用3%~5% (体积分数) 甲醛固定,经藻类沉降器沉降,尼康倒置显微镜 (TS100) 观察藻属种类及数量。待滩涂露出后将温度计插入养殖实验场地滩涂表层下5~10 cm处测量滩温。

    此外,对实验海区的滩涂底质粒径组成进行调查,分别取3个实验海区滩涂表层下约10 cm处的底质200 g,烘干机烘干,以标准检验筛 (GB/T 6003.1—2012) 筛分,参考赵东波等[19]和刘志刚等[20]的海砂粒级分类标准,按优势粒组命名法分类:底质粒径2.00~0.50 mm为粗砂,0.50~0.25 mm为中砂,0.25~0.063 mm为细砂,粒径<0.063 mm为粉砂质黏土。

    钝缀锦蛤海区养成实验于2021年7月14日开始,至2022年7月19日结束,实验周期为370 d。海区养成前随机取30粒中培幼贝测量初始规格:平均壳长 (17.42±1.63) mm、平均壳宽 (12.01±1.16) mm、平均壳高 (6.95±0.73) mm、平均体质量 (0.87±0.24) g。海区养成与中培实验在相同地点,以筐式进行养殖。实验过程中塑料筐半埋入滩涂,筐上部的四周孔洞外露,并加盖孔径1 cm的塑料硬盖。海区养成实验设置4个密度组:D1 (20 粒·筐−1)、D2 (30 粒·筐−1)、D3 (40 粒·筐−1)、D4 (50 粒·筐−1),每个密度设25个平行组。实验期间,每2个月退大潮时对每个密度的3个平行组进行不重复取样测量,测量时将该平行组所有样本全部筛出,随机取30粒样本 (样本数量不足30粒时取该组的全部贝苗),测量壳长、壳宽、壳高和体质量,并计算成活率。每次取样测量时对当日该实验海区的环境因子进行调查,测量项目同海区中培期。

    钝缀锦蛤海区网笼吊养实验于2022年7月26日开始,至2023年1月20日结束,实验周期为178 d。海区网笼吊养前随机取30粒钝缀锦蛤1.5龄贝测量初始规格:平均壳长 (44.65±2.83) mm、平均壳宽 (18.27±1.25) mm、平均壳高 (30.06±1.82) mm、平均体质量 (14.48±2.84) g。海区网笼吊养实验地点在防城港市白龙珍珠湾海区,海区网笼吊养实验设置3个密度组: C1 (30 粒·笼−1)、C2 (40 粒·笼−1)、C3 (50 粒·笼−1),按密度梯度放入网笼内网袋基质中,每个密度设置12个平行组。网笼通过吊绳垂挂在珍珠贝养殖浮排上进行垂下式养殖,保持养殖网笼吊养深度为距水面约1.5 m (水深1.5 m)、笼间距20 cm。实验期间,每2个月退大潮时对每个密度的3个平行组进行不重复取样测量,测量方法及项目同海区养成。每次取样测量时对当日该实验海区的环境因子进行调查,项目同海区中培期。

    实验数据采用Excel 2019软件整理,结果以“平均值±标准差 ($ \overline { x}\pm s$)”表示,使用SPSS 26.0软件对数据进行单因素方差分析 (One-way ANOVA),通过Duncan's多重比较法检验不同组间的差异显著性 ,P<0.05表示差异显著。

    钝缀锦蛤海区中培实验结束时在3个不同海区的生长指标和成活率见图3。北海海区d1组除壳宽外的生长指标和成活率均为最高,其中壳长为 (16.51±1.07) mm、壳高为 (11.64±1.11) mm、体质量为 (0.72±0.18) g、成活率为 (78.23±0.92)%。北海海区各密度组的成活率显著高于钦州和防城港海区 (P<0.05),除北海海区d5组的生长指标 (壳长、壳宽及体质量) 小于防城港和钦州海区外,其他各密度组生长指标均大于防城港和钦州海区。北海海区d1—d4密度组的壳长、壳宽、壳高和体质量均高于钦州海区相应组,北海海区d3组各生长指标显著高于钦州和防城港海区 (P<0.05)。

    图  3  钝缀锦蛤在广西北部湾3个海区中培期的生长和成活率
    注:同一密度、不同海区大写字母不同表示组间差异显著(P<0.05),相同表示组间差异不显著 (P>0.05);同一海区、不同密度,小写字母不同表示组间差异显著(P<0.05),相同表示组间差异不显著 (P>0.05)。图4同此。
    Fig. 3  Growth and survival rate of nursery culture T. dorsatus in three sea areas of Beibu Gulf, Guangxi
    Note: Different capital letters for the same density and different sea areas represent significant differences between groups (P<0.05), while the same capital letters represents no significant differences between groups (P>0.05). Different lowercase letters for the same sea area, different densities represent significant differences between groups (P<0.05), while the same lowercase letters represents no significant differences between groups (P>0.05). The same case in Fig. 4.

    在海区水温27.50~30.30 ℃、滩温28.00~31.00 ℃、盐度28‰~29‰ 的条件下,中培实验结束时,北海海区d3组的成活率为 (76.30±1.01)%,与d1、d2组无显著性差异 (P>0.05),显著大于d4和d5组 (P<0.05)。d3组的壳长、体质量显著大于d4和d5组 (P<0.05) (图3)。结果表明,钝缀锦蛤在北海海区的适宜中培密度为2 000 粒·筐−1

    在海区水温27.50~31.70 ℃、滩温28.50~32.00 ℃、盐度27‰~28‰ 的条件下,防城港海区钝缀锦蛤各生长指标和成活率随着中培密度的增加呈下降趋势,中培实验结束时,d2组 (1500 粒·筐−1) 成活率为 (69.31±2.63)%,与d1组差异不显著 (P>0.05),但显著高于其余各密度组 (P<0.05)。d2组壳长、体质量显著高于d4和d5组 (P<0.05) (图3)。结果表明,钝缀锦蛤在防城港海区的适宜中培密度为1 500 粒·筐−1

    在海区水温29.50~30.00 ℃、滩温30.00~30.50 ℃、盐度23‰~27‰ 的条件下,中培实验结束时,钦州海区d2组成活率为66.31%,与d1组无显著性差异 (P>0.05),显著大于d3、d4和d5组 (P<0.05),且d2组壳长、体质量显著大于d4和d5组 (P<0.05) (图3)。结果表明,钝缀锦蛤在钦州海区的适宜中培密度为1 500 粒·筐−1

    3个海区结束养成时钝缀锦蛤各生长指标和成活率见图4,养成期钝缀锦蛤的体质量生长曲线见图5

    图  4  北部湾海区钝缀锦蛤养成期的生长和成活率
    Fig. 4  Growth and survival rate in adult culture of T. dorsatus in Beibu Gulf sea area
    图  5  海区养成期钝缀锦蛤各密度组体质量生长曲线
    Fig. 5  Body mass growth curves of each density group of adult culture of T. dorsatus in sea area

    北海D1组钝缀锦蛤的各生长指标和成活率均最高,其中壳长为 (55.65±3.98) mm、壳宽为 (22.09±1.55) mm、壳高为 (36.68±2.19) mm、体质量为 (27.98±4.98) g,成活率为 (95.00±5.00)%。Duncan's法多重比较分析表明,北海和防城港海区D1组的成活率无显著性差异 (P>0.05),但均显著大于钦州海区 (P<0.05),其余3个密度组的成活率在3个海区间差异显著 (P<0.05)。钝缀锦蛤在3个海区不同养殖密度下的壳长、壳宽、壳高、体质量差异均达到显著水平 (P<0.05),不同养殖密度下钝缀锦蛤各生长指标和成活率为北海>防城港>钦州。

    在海区水温19.50~27.90 ℃、滩温20.80~29.80 ℃、盐度27‰~31‰ 的条件下,北海海区养殖密度为D3时,产量最大 (平均861.03 g·筐−1)。D3组的壳长、体质量、成活率显著高于D4组 (P<0.05)。北海海区各密度组的体质量增长速度变化趋势相同,第6—第10个月水温较低 (19.50~25.30 ℃),体质量增长率从0.09 g·d−1降至0.04 g·d−1,第10—第12个月水温回升 (25.30~27.20 ℃),体质量增长速度加快,由0.04 g·d−1上升至0.08 g·d−1。结果表明,北海海区适宜的养成密度为40 粒·筐−1

    在海区水温17.80~29.30 ℃、滩温19.00~30.00 ℃、盐度26‰~27‰ 的条件下,防城港海区养殖密度为D3 (40 粒·筐−1) 时,平均产量最大 (390.15 g·筐−1),壳长为 (44.36±2.15) mm、体质量为 (14.45±2.55) g、成活率为 (67.50±6.61)%。Duncan's法多重比较分析表明,D3组的成活率显著大于D4组 (P<0.05)。防城港海区第8个月水温最低 (17.80 ℃),第8—第10个月各密度组的体质量增长率降低,由0.05 g·d−1降为0.01 g·d−1,第10—第12个月水温升高,体质量增长速度加快,由0.01 g·d−1上升至0.04 g·d−1。结果表明,防城港海区适宜的养成密度为40 粒·筐−1

    在海区水温16.70~29.80 ℃、滩温17.30~30.40 ℃、盐度23‰~26‰ 的条件下,钦州海区D1组的产量最高,平均122.65 g·筐−1,且成活率显著大于D2、D3和D4组 (P<0.05)。钦州海区第6—第10个月水温降低,4个密度组的体质量增长速度由0.03 g·d−1下降至0.01 g·d−1,第10—第12个月水温回升,体质量增长率回升,从0.01 g·d−1上升至0.03 g·d−1。结果表明,钦州海区钝缀锦蛤适宜的养成密度为20 粒·筐−1

    在水温18.50~28.90 ℃、盐度26‰~27‰条件下,钝缀锦蛤在防城港白龙珍珠湾海区网笼吊养模式的各生长指标和成活率见表1。C1组各生长指标均最大,C3组成活率最高、产量最大;各个密度组的壳长、壳宽、壳高、体质量及成活率无显著性差异 (P>0.05)。结果表明,网笼吊养钝缀锦蛤成活率均在90%以上,且产量随着密度的增加而升高。

    表  1  防城港海区钝缀锦蛤网笼吊养模式的生长和成活率
    Table  1  Growth and survival rate of net cage culture of T. dorsatus in Fangchenggang sea area
    密度 Density/
    (grains·cage−1)
    壳长
    Final SL/mm
    壳宽
    Final SW/mm
    壳高
    Final SH/mm
    体质量
    Final BG/g
    成活率
    Survival rate/%
    产量
    Yield/g
    C1 (30) 52.00±2.72a 20.67±1.23a 34.34±1.78a 22.13±3.73a 95.56±3.85a 634.42
    C2 (40) 51.97±3.95a 20.39±1.64a 33.69±2.61a 21.77±4.96a 94.17±1.44a 901.20
    C3 (50) 51.57±2.39a 20.67±1.29a 33.66±1.94a 21.98±3.46a 96.00±2.00a 1020.96
    注:同列上标小写字母不同表示组间差异显著(P<0.05),相同表示组间差异不显著 (P>0.05)。 Note: Different and the same lowercase letters within the same column represent significant differences (P<0.05) and insignificant differences (P>0.05) between groups, respectively.
    下载: 导出CSV 
    | 显示表格

    为进一步分析钝缀锦蛤在广西北部湾3个养殖海区中培及养成期的生长存活差异,对3个养殖海区的水温、滩温、盐度、透明度、pH、溶解氧 (DO)、滩涂底质以及浮游植物藻属丰度进行周年取样调查。3个海区理化因子调查结果见图6图7图8。实验期间北海海区水温19.50~30.30 ℃ (平均26.10 ℃)、滩温20.80~31.00 ℃ (平均27.40 ℃)、盐度27‰~31‰ (平均28.9‰)、pH 7.88~8.41 (平均8.17)、透明度1.50~1.70 m (平均1.61 m)、DO 6.84~8.33 mg·L−1 (平均7.55 mg·L−1);防城港海区水温17.80~31.70 ℃ (平均25.50 ℃)、滩温19.00~32.00 ℃ (平均26.30 ℃)、盐度26‰~28‰ (平均26.7‰)、pH 7.81~8.34 (平均8.07)、透明度1.50~1.65 m (平均1.55 m)、DO 6.73~8.13 mg·L−1 (平均7.33 mg·L−1);钦州海区水温16.70~30.00 ℃ (平均25.10 ℃)、滩温17.30~30.50 ℃ (平均26.10 ℃)、盐度23‰~27‰ (平均25‰)、pH 7.65~8.32 (平均7.97)、透明度1.35~1.50 m (平均 1.41 m)、DO 5.28~7.64 mg·L−1 (平均7.08 mg·L−1)。北海海区的年平均水温、滩温、盐度、透明度、pH、DO略高于防城港和钦州海区。北海海区滩温、水温、pH的变化最小。结果表明,北海海区水质环境变化较稳定,防城港次之,钦州较不稳定。

    图  6  3 个海区水温和滩温的周年变化
    Fig. 6  Annual changes of water temperature and beach temperature in three sea areas
    图  7  3 个海区盐度和溶解氧的周年变化
    Fig. 7  Annual changes of salinity and dissolved oxygen in three sea areas
    图  8  3个海区pH和透明度的周年变化
    Fig. 8  Annual changes of pH and transparency in three sea areas

    3个海区的滩涂底质调查结果见表2。钦州海区滩涂底质粒径以细砂为优势粒径组 (77.33%);北海海区滩涂底质粒径以粗砂 (49.00%) 和中砂 (29.05%) 为优势粒径组;防城港海区滩涂底质粒径以中砂 (59.28%) 和细砂 (36.17%) 为优势粒径组。

    表  2  不同海区滩涂底质粒径组成比例
    Table  2  Composition of sand beach sediment particles size in different sea areas %
    海区
    Sea area

    Gravel
    粗砂
    Coarse sand
    中砂
    Medium sand
    细砂
    Fine sand
    粉砂质黏土
    Silty clay
    > 2.00 mm2.00~0.50 mm0.50~0.25 mm0.25~0.063 mm< 0.063 mm
    钦州 Qinzhou 3.70 10.15 6.35 77.33 2.45
    北海 Beihai 0.75 49.00 29.05 21.00 0.20
    防城港 Fangchenggang 0.32 3.51 59.28 36.17 0.72
    下载: 导出CSV 
    | 显示表格

    3个海区浮游植物藻属种类和生物量见图9。在北海海区的藻属种类共26种,其中24种属于硅藻门、2种属于甲藻门;防城港海区藻属种类23种,其中20种属于硅藻门、1种属于甲藻门、2种属于蓝藻门;钦州海区的藻属种类22种,其中20种属于硅藻门、1种属于甲藻门、1种属于蓝藻门。北海、防城港、钦州海区的浮游藻类生物量分别为 (13 775±613.90)、 (6 700±1 644.69)和 (4 825±601.56) 个·L–1。结果表明,海区浮游植物藻属的种类和数量均为北海>防城港>钦州。

    图  9  3 个不同养殖海区浮游植物藻属的种类和数量
    Fig. 9  Species of algal genus and quantity of phytoplankton in three different mariculture areas

    已有研究表明,养殖密度是影响贝类生长和存活的重要因素。养殖密度越大,贝类排泄物越多,更容易造成养殖环境污染,影响贝类生存空间。高密度也会加大贝类对空间和饵料的竞争、降低成活率;此外高密度养殖还会降低贝类抵御环境变化及病原侵袭的能力,但养殖密度过小则会降低空间利用率,减少经济效益。刘志刚等[20]研究表明皱肋文蛤 (Meretrix lyrate) 中培密度对种苗生长速度和成活率存在显著性影响 (P<0.05),即密度越低,生长速度越快,成活率越高。李永仁等[21]研究表明天津沿海菲律宾蛤仔 (Ruditapes philippinarum) 幼贝壳长大于22 mm时,高密度组生长缓慢,差异显著 (P<0.05),成贝养成期当壳长达到30 mm以上时,4 000 粒·m−2 密度组生长速度显著低于其他低密度组,死亡率高于50%,其他低密度组死亡率介于15%~22%。潘渊博等[22]对福建定海湾菲律宾蛤仔的研究表明,养殖密度影响菲律宾蛤仔的产量、形态构建和生长规律。于佐安等[23] 研究表明虾夷扇贝 (Pationopecten yessoensis) 低密度养殖组 (10和15枚·层−1) 的平均壳高大于高密度养殖组 (20、25和30枚·层−1),而累积死亡率则低于高密度养殖组。刘德经等[24] 研究表明西施舌 (Coelomactra antiquata) 稚贝的生长与放养密度关系密切,最终平均体长、平均体质量均随放养密度的增加而降低。何毛贤等[25]对马氏珠母贝 (Pinctada martensii) 的研究表明,在一定实验范围内养殖密度对体质量的增长有显著影响 (P<0.05),高密度组 (200只·笼−1) 的生长较低密度组慢。

    本研究表明,北海、防城港和钦州海区养成密度为40、40、20 粒·筐−1时产量达到峰值,北海和防城港海区产量随着筐式养成密度的增加呈先增后降的峰值变化趋势。栗志民等[26]研究养成期墨西哥湾扇贝 (Argopecten irradians concentricus) 养殖密度对产量的影响,结果表明大板笼养殖密度为60个·层−1、小板笼养殖密度为50个·层−1时两种笼具单层产量最大,且随着密度的升高产量呈先增后降的峰值变化趋势。

    本研究表明,当养殖密度超过筐式养殖最适密度时,钝缀锦蛤的生长和成活率均显著降低 (P<0.05),且在适宜筐式养殖密度范围内,养殖环境对钝缀锦蛤的生长和存活均产生影响,当养殖密度大于适宜养殖密度时,养殖密度成为制约钝缀锦蛤生长和存活的主要因素。廖锐等[27]研究了养殖密度对鱼类福利的影响,表明降低成活率是为了避免未来更为剧烈的种间竞争而采取的自我保护性调节措施。

    本实验中培期北海海区适宜的中培密度为2000 粒·筐−1,钦州和防城港海区的为1500 粒·筐−1;养成期北海和防城港海区适宜的养成密度为40 粒·筐−1,钦州海区为20 粒·筐−1。结果表明,养殖密度对钝缀锦蛤中培和养成期的生长和存活产生显著影响 (P<0.05),即随着密度的增加,其生长和成活率均呈下降趋势。3个实验海区的适宜中培和养成密度有所不同,可能与其生长的环境因子不同有关。结合本实验的海区环境调查结果表明,北海海区水温、滩温和pH波动较小,水质环境更稳定,浮游植物藻属种类及数量占比相对钦州、防城港海区更多,底质粒径较大,更适宜钝缀锦蛤的生长存活。

    目前广西北部湾海域钝缀锦蛤多以开放式底播养殖为主,养殖过程存在贝苗逃逸、敌害袭扰等问题,且采捕方式多为人工耙取,劳动强度大,采捕过程容易出现贝壳破碎等现象,采捕效率低、损耗大。为寻找便于管理且高效的养殖模式,本研究以筐式养殖模式开展中培、养成实验,并首次采用网笼进行钝缀锦蛤吊养实验,探讨了2种养殖模式对其生长和存活的影响。

    本实验表明,筐式养殖钝缀锦蛤在北海海区中培期成活率可达78.23%,海区养成成活率可达95.00%,并且筐式养殖采收方便。杨凌等[16]采用筐式中培钝缀锦蛤+吊养马氏珠母贝防风浪养殖模式在150 m2的滩涂上养殖60 d获经济效益达3 260.93元。筐式养殖模式在其他底播贝类中已有研究报道,王超奇等[28]使用塑料筐 (长510 mm×宽350 mm×高170 mm) 和塑料箱 (长500 mm×宽350 mm×高230 mm) 进行施氏獭蛤 (Lutraria sieboldii) 海区中培实验,使用塑料筐 (长580 mm×宽420 mm×高290 mm ) 进行海区养成实验,结果表明,中培期北海和防城港海区箱式养殖表现更佳,钦州海区筐式养殖模式表现更佳。养成期适宜在北海进行,适宜养成密度为40 粒·筐−1。蔡德建等[29]使用长方形塑料筐 (长50 cm×宽35 cm×高15 cm)进行大獭蛤 (Lutmria maxima) 人工浅海沉筐式养殖,表明浅海沉筐中间培育生长速度快,成活率高。因此,筐式养殖模式可作为北部湾海域钝缀锦蛤的高效海区中培和养成模式。

    滩涂调查表明,北海海区大部分滩涂平缓,退大潮几乎全部露空,北海网箱养殖一般都设在浅海区,不在滩涂低潮区;钦州竹筏式吊养牡蛎的海区为半咸淡水,盐度偏低,不适合钝缀锦蛤生长;防城港白龙珍珠湾滩涂低潮区常年浮排网笼吊养马氏珠母贝,吊养水深合适且盐度高,适宜钝缀锦蛤生长。本研究表明,网笼吊养可降低风浪以及退潮露空时间对钝缀锦蛤摄食及生长的影响,与筐式养殖相比,网笼吊养使贝类摄食、呼吸等不受退潮露空时间影响,可以缩短养殖时间,且其排泄物更容易扩散,避免了排泄物沉积对其生长环境的污染;同时采用轻便的浮石作为基质替代海沙,钝缀锦蛤加工时不需要进行吐沙处理,提高了钝缀锦蛤的品质。本研究结果表明,网笼吊养钝缀锦蛤成活率均在90%之上,最高可达96%。田园等[17]研究表明筏式养殖菲律宾蛤仔的成活率在90%以上,秋、冬季的成活率显著高于底播养殖 (P<0.05),这与本研究结果相似。陈瑞芳等[30]使用定制PVC养殖筐 (长80 cm×宽60 cm×高40 cm) 进行织锦巴非蛤 (Paphia textile) 海区悬挂养殖,结果表明海区悬挂模式可提高织锦巴非蛤的成活率,由30%左右提升至50%以上。

    海水的理化性质对水生生物的生长分布起着至关重要的作用[31]。本实验在广西北部湾海域钦州犀牛脚、北海沙虫坪、防城港白龙珍珠湾3个海区同时开展钝缀锦蛤的生长对比实验,结果表明,海区中培和养成期相同密度条件下钝缀锦蛤的生长和成活率均为北海>防城港>钦州。

    1) 水质因子。实验期间钦州海区周年水温16.70~30.00 ℃、滩温17.30~30.50 ℃、盐度23‰~27‰;北海海区周年水温19.50~30.30 ℃、滩温20.80~31.00 ℃;盐度27‰~31‰;防城港海区周年水温17.80~31.70 ℃、滩温19.00~32.00 ℃、盐度26‰~28‰。已有研究表明钝缀锦蛤生长的半致死低温为 (9.20±0.15) ℃、高温为 (34.25±0.09) ℃,适宜盐度为23.58‰~45.15‰,最适盐度为33.90‰[5]。本实验3个养殖海区的水温、滩温和盐度均在其适宜的存活范围。本研究表明,北海海区的水质环境相对较钦州和防城港海区更稳定,钝缀锦蛤在水质较稳定的环境中生长更快,成活率更高。

    2) 底质粒径。底质粒径会影响贝类的生长、存活及分布[20,32]。本实验滩涂底质分析结果表明,北海滩涂底质以粗砂为主 (49.00%),防城港以中砂为主 (59.28%),钦州海区以细砂为主 (77.33%),中培期在北海和防城港实验海区效果较好,养成期在北海海区效果较好。此外,海区网笼吊养使用比重为0.45 g·cm−3、孔隙率为75%~80%、颗粒粒径为1.50~2.50 cm的浮石作为基质,钝缀锦蛤的成活率始终保持在90%以上。本实验表明钝缀锦蛤适宜在底质粒径为0.20 mm 以上的滩涂进行海区中培,在底质粒径为0.50 mm以上的滩涂进行海区养成。刘亦涵等[33]研究表明缀锦蛤 (Tapes literatus) 具有明显的底质选择性,大粒径 (15~30 目) 的基底对缀锦蛤的生长有促进作用,在致密的细砂粒径 (80~100目) 基质几乎没有生长。该研究结果与本实验结果相似。

    3) 藻类的种类和数量。本研究表明,海区浮游植物的数量对钝缀锦蛤的生长具有显著性影响 (P<0.05)。3个养殖海区的周年浮游植物藻属种类和数量分布表明,北海养殖海区藻属种类最为丰富,钦州与防城港养殖海区相近。在藻属数量上北海海区占优势,其次为防城港海区,钦州海区最少。这与黄亚楠等[34]研究墨西哥湾扇贝选育系F7 在广西北部湾海域的生长比较,对广西北部湾海域的浮游植物数量调查结果相同。

    4) 干露胁迫。邱清波等[35]研究了不同干露胁迫时间对缀锦蛤耗氧率的影响,结果表明缀锦蛤在干露72 h后耗氧率低于正常水平。本实验筐式养殖模式各实验地点位于潮间带低潮区,由于地理位置不同,退大潮后露空时间不同,北海沙虫坪海区退潮后露空时间约6 h,防城港白龙珍珠湾海区约7 h,钦州犀牛脚渔港海区约8 h。采用网笼吊养模式在浮排上进行垂下养殖,吊笼距离水面1.50 m,不存在露空时间,但浮排吊养模式抗风浪能力相对较差。滩涂筐式养殖模式的抗风浪能力相对较强,但其成活率相对较低,可能与退大潮时间在白天时有关,太阳直晒导致水温、滩温较高,露空胁迫影响了其摄食、生长和呼吸。黄洋等[2]研究湛江雷州半岛西部乌石及遂溪草潭海区的野生钝缀锦蛤生态习性,表明耐干露能力与气温成反比,气温越高,耐干露能力越弱,这与本实验结果相似。

    本研究表明,钝缀锦蛤在广西北部湾3个海域低潮区适宜采用筐式进行海区中培和养成,其中北海海区适宜的中培密度为2 000 粒·筐−1、养成密度为40 粒·筐−1;防城港海区适宜的中培密度为1 500 粒·筐−1、养成密度为40 粒·筐−1;钦州海区适宜的中培密度为1 500 粒·筐−1、养成密度为20 粒·筐−1。北海中培和养成效果优于防城港和钦州,此外,在防城港白龙珍珠湾海区开展大规格钝缀锦蛤网笼吊养是可行的,各个密度组生长指标及成活率无显著性差异 (P>0.05) 且均在90%以上。本研究比较研究了筐式中培、养成和网笼吊养2种模式及其养殖效果,可为钝缀锦蛤在广西北部湾海区大规模养殖提供参考。

  • 图  1   一阶差分后的自相关函数 (ACF) 和偏相关函数 (PACF) 图

    注:a—b. 延绳钓赤道以北;c—d. 延绳钓赤道以南;e—f. 围网赤道以北;g—h. 围网赤道以南;图2—图4同此。

    Figure  1.   ACF and PACF diagrams after first difference

    Note: a–b. North of the equator for longline; c–d. South of the equator for longline; e–f. North of the equator for purse seine; g–h. South of the equator for purse seine. The same case in Fig. 2–Fig. 4.

    图  2   常规ARIMA模型残差序列自相关函数 (ACF) 和偏相关函数 (PACF) 图

    Figure  2.   ACF and PACF diagrams of general ARIMA model residual sequence

    图  3   动态 ARIMA模型拟合残差检验 (ACF 和 Q-Q 图)

    Figure  3.   Residual test of dynamic ARIMA model (ACF and Q-Q diagrams)

    图  4   常规 ARIMA 模型和动态 ARIMA 模型的拟合优度比较

    Figure  4.   Comparison of goodness of fit between general ARIMA model and dynamic ARIMA model

    图  5   1990—2020 年中西太平洋不同渔业黄鳍金枪鱼年 CPUE 模型拟合值与实际值关系

    a. 延绳钓赤道以北;b. 延绳钓赤道以南;c. 围网赤道以北;d. 围网赤道以南。

    Figure  5.   Relationship between fitted and actual values of annual CPUE model for yellowfin tuna from different fisheries in Western and Central Pacific during 1990−2020

    a. North of the equator for longline; b. South of the equator for longline; c. North of the equator for purse seine; d. South of the equator for purse seine.

    表  1   不同区域渔业的水平时间序列和一阶拆分序列的单位根检验表

    Table  1   Unit root test of horizontal sequence and difference sequence for various regional fisheries

    时间序列    
    Time series    
    数据处理
    Data processing
    单位根检验
    Unit root test
    显著性水平 Significance level
    1%5%10%
    延绳钓赤道以北
    North of equator for longline
    水平序列 −1.296 −2.62 −1.95 −1.61
    一阶差分序列 4.922 −4.15 −3.50 3.18
    延绳钓赤道以南
    South of equator for longline
    水平序列 −0.586 −2.62 −1.95 −1.61
    一阶差分序列 5.571 −4.15 −3.50 3.18
    围网赤道以北
    North of equator for purse seine
    水平序列 −0.746 −2.62 −1.95 −1.61
    一阶差分序列 4.818 −4.15 −3.50 3.18
    围网赤道以南
    South of equator for purse seine
    水平序列 −0.808 −2.62 −1.95 −1.61
    一阶差分序列 4.818 −4.15 −3.50 3.18
    注:粗体字表示该序列为平稳时间序列。 Note: The boldface numbers indicate that the series is a stationary time series.
    下载: 导出CSV

    表  2   ARIMA 时间序列全部可能模型以及参数评价

    Table  2   All possible models and evaluation of parameters of ARIMA time series

    时间序列
    Time series
    一般 ARIMA 模型
    General ARIMA model
    模型赤池信息量准则
    AIC
    小样本校正的 AIC
    AICc
    贝叶斯信息准则
    BIC
    延绳钓赤道以北
    North of equator for longline
    (0,1,0) 81.60 81.74 83.00
    (0,1,1) 83.06 83.51 85.87
    (1,1,0) 83.29 83.73 86.09
    (1,1,1) 83.69 84.62 87.9
    延绳钓赤道以南
    South of equator for longline
    (0,1,1) 98.58 99.03 101.39
    (0,1,4) 97.14 99.64 104.15
    (1,1,1) 100.57 101.49 104.77
    (1,1,4) 98.51 102.16 106.92
    围网赤道以北
    North of equator for purse seine
    (0,1,0) 116.51 116.65 117.91
    (0,1,1) 110.06 110.51 112.86
    (1,1,0) 109.95 110.39 112.75
    (1,1,1) 111.62 112.54 115.82
    围网赤道以南
    South of equator for purse seine
    (2,1,0) 103.42 104.35 107.63
    (0,1,2) 102.8 103.72 107.00
    (2,1,2) 103.93 106.43 110.94
    (0,1,0) 110.24 110.39 111.64
    注:粗体字表示该模型为最优模型。 Note: The boldface numbers indicate the optimal model.
    下载: 导出CSV
  • [1]

    WU Y L, LAN K W, TIAN Y. Determining the effect of multiscale climate indices on the global yellowfin tuna (Thunnus albacares) population using a time series analysis[J]. Deep Sea Res II, 2020, 175: 104808-104819. doi: 10.1016/j.dsr2.2020.104808

    [2]

    MCCLATCHIE S. Sardine biomass is poorly correlated with the Pacific decadal oscillation off California[J]. Geophysical Res Lett, 2012, 39(13): 1-6.

    [3]

    BRUCE F P, MONICA P R. Climate change impacts on fisheries and aquaculture: a global analysis[M]. New York: John Wiley & Sons, 2017: 523-568.

    [4] 王雅萌, 陈芃, 陈新军. 气候变化下西北太平洋大海洋生态系海表面温度特征分析[J]. 上海海洋大学学报, 2021, 30(5): 875-883. doi: 10.12024/jsou.20200603074
    [5] 郭爱, 陈新军, 范江涛. 中西太平洋鲣鱼时空分布及其与ENSO关系探讨[J]. 水产科学, 2010, 29(10): 591-595. doi: 10.3969/j.issn.1003-1111.2010.10.006
    [6] 周为峰, 陈亮亮, 崔雪森, 等. 异常气候下温跃层及时空因子中西太平洋黄鳍金枪鱼渔场分布的影响[J]. 中国农业科技导报, 2021, 23(10): 192-201.
    [7]

    SYAMSUDDIN M L, SAITOH SI, HIRAWAKE T, et al. Effects of El Niño-Southern Oscillation events on catches of bigeye tuna (Thunnus obesus) in the Eastern Indian Ocean off Java[J]. Fish Bull, 2013, 111(2): 175-188. doi: 10.7755/FB.111.2.5

    [8]

    ERAUSKIN-EXTRAMIANA M, ARRIZABALAGA H, HOBDAY A J, et al. Large-scale distribution of tuna species in a warming Ocean[J]. Global Change Biol, 2019, 25(6): 2043-2060. doi: 10.1111/gcb.14630

    [9]

    Food and Agriculture Organization of the United Nations. FAO STAT statistics database[R]. Roma: FAO, 2018: 284-290.

    [10]

    Western and Central Pacific Fisheries Commission. Tuna fishery yearbook[R]. The Federated States of Micronesia: WCPFC, 2019: 139-144.

    [11]

    OLSON R J, YOUNG J W, MENARD F, et al. Bioenergetics, trophicecology, and niche separation of tunas[J]. Adv Mar Biol, 2016, 74: 199-344.

    [12] 杨胜龙, 靳少非, 吴祖立, 等. 太平洋金枪鱼渔场关键次表层环境变量的季节变化[J]. 海洋科学, 2015, 39(5): 36-46.
    [13] 朱国平, 许柳雄. 东太平洋金枪鱼延绳钓大眼金枪鱼渔场与表层温度之间的关系[J]. 海洋环境科学, 2007(4): 333-336. doi: 10.3969/j.issn.1007-6336.2007.04.009
    [14] 曹晓怡, 周为峰, 樊伟, 等. 大眼金枪鱼渔场与环境关系的研究进展[J]. 海洋渔业, 2008, 30(2): 176-182. doi: 10.3969/j.issn.1004-2490.2008.02.012
    [15] 杨胜龙, 张胜茂, 蒋兴伟, 等. 热带大西洋大眼金枪鱼和黄鳍金枪鱼渔场温跃层的时空变化特征[J]. 应用海洋学学报, 2013, 32(3): 349-357. doi: 10.3969/J.ISSN.2095-4972.2013.03.007
    [16] 李鹏, 许柳雄, 周成, 等. 中西太平洋金枪鱼围网自由鱼群渔场重心变动及其与南方涛动指数的关系[J]. 南方水产科学, 2020, 16(2): 20-27.
    [17] 宋利明, 李玉伟, 高攀峰. 帕劳群岛附近海域延绳钓渔场大眼金枪鱼 (Thunnus obesus)的环境偏好[J]. 海洋与湖沼, 2009, 40(6): 768-776. doi: 10.3321/j.issn:0029-814X.2009.06.015
    [18]

    SOYKAN C U, EGUCHI T, KOHIN S, et al. Prediction of fishing effort distributions using boosted regression trees[J]. Ecol Appl, 2014, 24(1): 71-83. doi: 10.1890/12-0826.1

    [19] 周永东, 徐汉祥, 潘国良, 等. 东海区鲳鱼、小黄鱼资源量及其持续渔获量的估算[J]. 浙江海洋学院学报(自然科学版), 2013, 32(1): 1-5.
    [20] 宋大德, 汪金涛, 陈新军, 等. 时间序列分析模型在黄海南部小黄鱼资源量预测中的应用[J]. 海洋学报, 2020, 42(12): 26-33.
    [21]

    MONLLOR-HURTADO A, PENNINO M G, SANCHEZ-LIZASO J L. Shift in tuna catches due to ocean warming[J]. PLoS One, 2017, 12(6): 178-196.

    [22] ROBERT H S, DAVID S S. 时间序列分析及其应用[M]. 李红成, 译. 北京: 机械工业出版社, 2020: 202-230.
    [23] 李颖若, 韩婷婷, 汪君霞, 等. ARIMA 时间序列分析模型在臭氧浓度中长期预报中的应用[J]. 环境科学, 2021, 42(7): 3119-3126.
    [24] 乔宇, 闫振飞, 冯承莲, 等. 几种典型模型在物种敏感度分布中的应用和差异分析[J]. 环境工程, 2021, 39(10): 85-95. doi: 10.13205/j.hjgc.202110012
    [25]

    FOSTER R G, ROENNEBERG T. Human responses to the geophysical daily, annual and lunar cycles[J]. Curr Biol, 2008, 18(17): R784-R794. doi: 10.1016/j.cub.2008.07.003

    [26] 袁红春, 高子玥, 张天蛟. 基于改进的XGboost模型预测南太平洋长鳍金枪鱼资源丰度[J]. 海洋湖沼通报, 2022, 44(2): 112-120.
    [27] 杨旭, 黄进, 秦泽宇, 等. 基于多尺度特征融合的人群计数算法[J]. 计算机系统应用, 2022, 33(1): 226-235.
    [28] 章贤成, 汪金涛, 陈新军. 基于BP神经网络的阿根廷滑柔鱼资源CPUE标准化研究[J]. 渔业科学进展, 2022, 43(2): 11-20.
    [29] 杨彩莉, 杨晓明, 朱江峰. 不同类型厄尔尼诺事件中环境因子对中西太平洋金枪鱼围网鲣分布响应[J]. 南方水产科学, 2021, 17(3): 8-18. doi: 10.12131/20210014
    [30] 王少琴, 许柳雄, 朱国平, 等. 中西太平洋金枪鱼围网的黄鳍金枪鱼CPUE时空分布及其与环境因子的关系[J]. 大连海洋大学学报, 2014, 29(3): 15-18.
    [31]

    WHITNEY J E, AL-CHOKHACHY R, BUNNELL D B, et al. Physiological basis of climate change impacts on North American inland fishes[J]. Fisheries, 2016, 41(7): 332-345. doi: 10.1080/03632415.2016.1186656

    [32]

    LEHODEY P, ALHEIT J, BARANGE M, et al. Climate variability, fish, and fisheries[J]. J Climate, 2006, 19(20): 5009-5030. doi: 10.1175/JCLI3898.1

    [33]

    BRANDER K M. Global fish production and climate change[J]. Proc Nat Acad Sci USA, 2007, 104(50): 19709-19714. doi: 10.1073/pnas.0702059104

    [34]

    LEHODEY P, SENINA I, CALMETTES B, et al. Modelling the impact of climate change on Pacific skipjack tuna population and fisheries[J]. Climatic Change, 2012, 119(1): 95-109.

    [35] 石肖飞, 王啸, 王佚兮, 等. 热带中西太平洋海域黄鳍金枪鱼的摄食生物学特性[J]. 南方水产科学, 2022, 18(1): 43-51. doi: 10.12131/20210140
    [36]

    TREMBLAY-BOYER L, MCKECHNIE S, PILLING G, et al. Stock assessment of yellowfin tuna in the Western and Central Pacific Ocean, WCPFC-SC13-2017/SA-WP-06[R]. Rarotonga, Cook Islands: WCPFC, 2017: 20-32.

    [37] 代丹娜, 刘洪生, 戴小杰, 等. ENSO现象与东太平洋黄鳍金枪鱼围网CPUE时空分布的关系[J]. 上海海洋大学学报, 2011, 20(4): 571-578.
    [38]

    SUMAILA U R, CHEUNG W W, LAM V W, et al. Climate change impacts on the biophysics and economics of world fisheries[J]. Nat Climate Change, 2011, 1(9): 449-456. doi: 10.1038/nclimate1301

    [39]

    SUN C, WANG W R, YEH S. Reproductive ecology of yellowfin tuna in the Central and Western Pacific Ocean, WCPFC-SC1 BI WP-1[R]. Noumea New Caledonia: WCPFC, 2005: 15-26.

    [40]

    MUHLING B A, LAMKIN J T, ALEMANY F, et al. Reproduction and larval biology in tunas, and the importance of restricted area spawning grounds[J]. Fish Pacific, 2017, 27(4): 697-732.

    [41] 蔡格菁, 傅海彬, 蒋仁斌, 等. 基于ARIMA模型的渔业经济预测及其优化[J]. 计算机与现代化, 2019(4): 88-91.
    [42] 李勇, 冯家成, 李娜, 等. 太湖水体 Chl-a 预测模型 ARIMA 在引排水方案优化中的应用[J]. 环境工程, 2022, 40(10): 72-79.
    [43]

    FREE C M, THORSON J T, PINSKY M L, et al. Impacts of histori-cal warming on marine fisheries production[J]. Science, 2019, 363(6430): 979-983. doi: 10.1126/science.aau1758

    [44]

    PLAGANYI E. Climate change impacts on fisheries[J]. Science, 2019, 363(6430): 930-931. doi: 10.1126/science.aaw5824

    [45]

    HAZEN E L, JORGENSEN S, RYKACZEWSKI R R, et al. Predicted habitat shifts of Pacific top predators in a changing climate[J]. Nat Climate Change, 2013, 3(3): 234-238. doi: 10.1038/nclimate1686

  • 期刊类型引用(6)

    1. 陈颖,赵菊花,杨羽,王洁,杜秀君,刘海霞. 茶多酚调节AMPK/SIRT1/PGC-1信号通路对光老化小鼠模型氧化损伤的影响. 西部医学. 2025(02): 205-211 . 百度学术
    2. 宋玉琼,胡晓,李来好,陈胜军,相欢,杨贤庆,吴燕燕. 食物蛋白源酪氨酸酶抑制肽的研究进展. 陕西师范大学学报(自然科学版). 2023(03): 16-28 . 百度学术
    3. 姜红伶. 强脉冲光与CO_2点阵激光治疗面部皮肤光老化疗效对比分析. 中国医疗美容. 2022(02): 25-29 . 百度学术
    4. 许锋,陈得科,陈晓刚,聂金梅,马志浩,孙恢礼,陈忻. 坛紫菜防晒剂MAAs提取、鉴定及在化妆品体系的稳定性研究. 南方水产科学. 2022(04): 156-162 . 本站查看
    5. 罗艳,黄权新,蔡捷. 牡蛎酶解产物的种类、生物活性及应用研究进展. 中国食物与营养. 2022(11): 49-53 . 百度学术
    6. 张雪妍,秦小明,林海生,曹文红,郑惠娜,高加龙,章超桦. 牡蛎酶解超滤组分对TM4小鼠睾丸支持细胞的氧化损伤保护作用. 南方水产科学. 2021(05): 118-125 . 本站查看

    其他类型引用(10)

推荐阅读
芮氏刻肋海胆胚胎及幼体发育过程
余炜欢 et al., 南方水产科学, 2025
基于基尔霍夫近似模型的南海大黄鱼声学目标强度研究
王文卓 et al., 南方水产科学, 2024
基于轴向特征校准和时间段网络的鱼群摄食强度分类模型研究
徐波 et al., 南方水产科学, 2024
海洋牧场建设效益评价研究进展与展望
袁华荣 et al., 南方水产科学, 2024
海洋中尺度涡对南太平洋长鳍金枪鱼延绳钓渔获量的影响
ZHENG Chunwen et al., JOURNAL OF SHANGHAI OCEAN UNIVERSITY, 2025
大西洋黄鳍金枪鱼资源丰度不同预报模型的比较分析
丁鹏 et al., 大连海洋大学学报, 2024
Body roundness index improves the predictive value of cardiovascular disease risk in hypertensive patients with obstructive sleep apnea: a cohort study
Cai, Xintian et al., CLINICAL AND EXPERIMENTAL HYPERTENSION, 2023
Enhanced impact of the aleutian low on increasing the central pacific enso in recent decades
Chen, Shangfeng et al., NPJ CLIMATE AND ATMOSPHERIC SCIENCE, 2023
Intra-annual flow distribution of the rivers in the yesil river basin
WATER CONSERVATION AND MANAGEMENT, 2024
Use of ephemeroptera as a bioindicator for assessing surface water quality in indonesia
WATER CONSERVATION AND MANAGEMENT, 2024
Powered by
图(5)  /  表(2)
计量
  • 文章访问数:  498
  • HTML全文浏览量:  139
  • PDF下载量:  59
  • 被引次数: 16
出版历程
  • 收稿日期:  2023-01-14
  • 修回日期:  2023-02-18
  • 录用日期:  2023-03-01
  • 网络出版日期:  2023-03-07
  • 刊出日期:  2023-08-04

目录

/

返回文章
返回