海洋学报

影响因子	不分模态阶段	第1模态	第2模态	第3模态	第4模态
亲体量	99	34	86	100	98
SSH	76	41	26	23	32
SSS	20	46	27	24	22
SST	19	40	32	31	25
PDO	16	28	34	41	24
El Niño	14	24	49	23	33

影响因子	不分模态阶段	第1模态	第2模态	第3模态	第4模态
亲体量	99	34	86	100	98
SSH	76	41	26	23	32
SSS	20	46	27	24	22
SST	19	40	32	31	25
PDO	16	28	34	41	24
El Niño	14	24	49	23	33

影响因子	模型1	模型2	模型3	模型4	模型5
注：√表示该模型中包含的因子。
亲体量	√	√	√	√	√
SSS			√
PDO					√
SSH	√		√	√	√
El Niño
SST				√
PMP(MCMC)	36%	9%	9%	9%	7%
PMP(Exact)	36%	9%	9%	9%	7%

影响因子	模型1	模型2	模型3	模型4	模型5
注：√表示该模型中包含的因子。
亲体量	√	√	√	√	√
SSS			√
PDO					√
SSH	√		√	√	√
El Niño
SST				√
PMP(MCMC)	36%	9%	9%	9%	7%
PMP(Exact)	36%	9%	9%	9%	7%

模型	方案a	方案b	方案c
模型1	24	19	22
模型2	12	13	11
模型3	6	9	9
模型4	6	6	6
模型5	5	6	6
累计	52	51	53

模型	方案a	方案b	方案c
模型1	24	19	22
模型2	12	13	11
模型3	6	9	9
模型4	6	6	6
模型5	5	6	6
累计	52	51	53

不同模态下气候变化对智利竹筴鱼补充量的影响

PDF下载

张畅 ¹ , 李纲 ¹ , 陈新军 ¹^,²^,³^,⁴^,⁵^,^*

海洋学报 | 论文 2021,43(9): 48-58

收起

海洋学报 | 论文 2021, 43(9): 48-58

不同模态下气候变化对智利竹筴鱼补充量的影响

全屏

张畅¹, 李纲¹, 陈新军¹^,²^,³^,⁴^,⁵^,^*

作者信息

¹上海海洋大学海洋科学学院，上海 201306

²青岛海洋科学与技术试点国家实验室海洋渔业科学与食物产出过程功能实验室，山东青岛 266237

³上海海洋大学大洋渔业资源可持续开发教育部重点实验室，上海 201306

⁴上海海洋大学国家远洋渔业工程技术研究中心，上海 201306

⁵上海海洋大学农业农村部大洋渔业开发重点实验室，上海 201306

张畅（1994－），女，浙江省金华市人，博士研究生，主要研究方向为渔业资源。E-mail：itschang@foxmail.com

通讯作者:

陈新军（1967－），男，教授，浙江省义乌市人，从事渔业资源与渔场学等研究。E-mail：xjchen@shou.edu.cn

Impact of climate change on recruitment of Trachurus murphyi based on different regimes

Chang Zhang¹, Gang Li¹, Xinjun Chen¹^,²^,³^,⁴^,⁵^,^*

Affiliations

¹College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China

²Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China

³Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai Ocean University, Shanghai 201306, China

⁴National Engineering Research Center for Oceanic Fisheries, Shanghai Ocean University, Shanghai 201306, China

⁵Key Laboratory of Oceanic Fisheries Exploitation, Ministry of Agriculture and Rural Affairs, Shanghai 201306, China

出版时间: 2021-09-25 doi: 10.12284/hyxb2021130

文章导航

摘要

收起

智利竹筴鱼(Trachurus murphyi)是东南太平洋重要的经济鱼类之一，其资源量受补充量影响明显，了解补充量状况对智利竹筴鱼资源可持续利用和科学管理具有重要意义。本文基于模态分析将1971−2017年间智利竹筴鱼补充量划分不同模态，运用贝叶斯模型平均法，分析海表面温度、海表面盐度、海表面高度、厄尔尼诺和太平洋年代际振荡5个环境因子在不同模态中对补充量的解释能力，并探讨模态变动对补充量预测的影响。结果表明，第1模态(1971−1980年)更多的受捕捞因素的干扰；第2模态(1981−1990年)厄尔尼诺对补充量变动的解释概率最高；第3模态(1991−2001年)解释概率最高的环境因子为太平洋年代际振荡；第4模态(2002−2015年)厄尔尼诺为解释概率最高的环境因子。对比不考虑模态变动的分析结果，两者存在明显差异，基于不同模态的分析结果对智利竹筴鱼补充量变动的解释更为合理。研究认为，智利竹筴鱼补充量变动受到多个环境因子的影响，在不同模态时期起主导作用的环境因子也不同，推测年代际太平洋年代际振荡冷暖期交替与厄尔尼诺现象可能是诱发智利竹筴鱼补充量发生模态转变的重要因素。建议在未来智利竹筴鱼资源评估与管理中，应该考虑不同的模态变化及其影响因子。

关键词

智利竹筴鱼 / 贝叶斯模型平均法 / 模态变动 / 补充量 / 东南太平洋

Abstract

收起

Chilean jack mackerel (Trachurus murphyi) is one of the important economic fisheries in the Southeast Pacific Ocean. The stock abundance is inherently affected by recruitment. To find out the status of recruitment is of great significance to the sustainable yield and scientific management of resources. The recruitment of Chilean jack mackerel between 1971 and 2017 is divided into different regimes based on regime test. Then, the Bayesian model averaging is applied to the relationship between recruitment and five environmental factors (sea surface temperature, sea surface salinity, sea surface height, El Niño and Pacific decadal oscillation), and the impact of regime change in recruitment forecast is also discussed. The results show that the first regime (1971−1980) is more disturbed by fishing; in the second regime (1981−1990), El Niño has the highest probability of explaining for the change of recruitment; the environmental factor with the highest explanatory probability is Pacific Decadal Oscillation in the third regime (1991−2001); in the fourth regime (2002−2015), El Niño most likely explains the change in recruitment; comparing the analysis of regime shift and no regime shift models, there is a significant difference between the two. However, the results, with different regimes, are more reasonable to explain the change of recruitment. Recruitment abundance is affected by multiple environmental factors, but in different regimes, the dominant environmental factor changes. The regime shifts of recruitment are induced by the alternation of cold and warm periods of interdecadal Pacific decadal oscillation, and El Niño. It is suggested that different regime shift and their impact factors should be considered in the future assessment and management of Chilean jack mackerel resources.

Key words

Trachurus murphyi / Bayesian model averaging / regime shift / recruitment / Southeast Pacific Ocean

引用本文

张畅, 李纲, 陈新军. 不同模态下气候变化对智利竹筴鱼补充量的影响. 海洋学报, 2021 , 43 (9) : 48 -58 . DOI: 10.12284/hyxb2021130

Chang Zhang, Gang Li, Xinjun Chen. Impact of climate change on recruitment of Trachurus murphyi based on different regimes[J]. Haiyang Xuebao, 2021 , 43 (9) : 48 -58 . DOI: 10.12284/hyxb2021130

正文

收起

1 引言

收起

智利竹筴鱼（Trachurus murphyi）是东南太平洋重要的经济种类，于20世纪70年代开始大规模开发和利用，1995年达到历年最高产量，为496万t，其后多年渔业资源出现衰退，2013年产量仅35万t，为近30 a以来的最低产量，2019年捕捞产量又恢复到63.8万t^[1]。作为沿海国家和远洋渔业国家的重要捕捞对象之一，其资源评估与管理工作受到世界各国和南太平洋区域渔业管理组织的极大关注，补充量作为研究种群资源量变动的一个重要因素，也是资源评估管理中不可忽略的重要因素。

通常，渔业资源的补充量除了取决于资源本身的再生产机制，还受到生态环境的制约^[2]。在目前已有研究中，主要通过两个方面考虑环境对补充量的影响：一是使用一个或多个已知环境变量，构建短期的补充量预测模型^[3]；二是考虑长期环境气候变化对补充量的影响，即考虑模态变动^[4]。自然的模态变动可能是剧烈且突然的，与海洋气候联系密切，它有别与捕捞因素的影响，在渔业管理中有着十分重要的意义^[5]。近年来，基于贝叶斯模型平均法构建的亲体−补充量模型在渔业中得到较好应用^[6-7]，该方法不仅能分析各环境因子对资源补充量是否具有长期、稳健的解释能力，还能用于补充量预测。为此，本研究拟通过贝叶斯模型平均法，分析在不同模态中的环境因子对智利竹筴鱼的亲体补充量影响及其内在变化规律，为其资源评估和科学管理提供参考。

2 材料与方法

收起

2.1 渔业与环境数据

（1）渔业数据。1971−2017年智利竹筴鱼亲体量与补充量源自南太平洋区域渔业管理组织（South Pacific Reginal Fisheries Management Organization, SPRFMO）评估报告^[1,8]，以1龄鱼的资源量作为补充量，以参与繁殖的群体资源量作为亲体量。

（2）海洋环境数据。考虑智利竹筴鱼种群结构复杂，结合已有文献研究^[9-10]，选取目前已知的产卵期（11月至翌年2月）最大产卵场（35°～45°S，90°～75°W）（图1）的海表面温度(Sea Surface Temperature, SST)、海表面盐度（Sea Surface Salinity, SSS）、海表面高度（Sea Surface Height, SSH）的平均值作为影响补充量的潜在影响因子，SST数据来源于https://oceanwatch.pifsc.noaa.gov/erddap/griddap/index.html?page=1&itemsPerPage=1000，SSH与SSS数据来源于http://apdrc.soest.hawaii.edu/las_ofes/v6/dataset?catitem=71。海洋尼诺指数（Oceanic Niño Index, ONI）以Niño3.4区中每年11月至翌年2月的SSTA平均值来表示，时间范围为1971−2017年，数据来源于https://ggweather.com/ONI/oni.htm。太平洋年代际振荡(Pacific Decadal Oscillation, PDO）可直接造成太平洋及周边地区气候的年代际变化，对中上层鱼类资源造成影响，为此，本研究选取每年11月至翌年2月的平均太平洋年代际振荡，作为表征其资源量变化的环境因子，时间范围为1971−2017年，数据来源于http://research.jisao.washington.edu/pdo/PDO.latest。

2.2 模态检验

基于连续t检验的模态变动分析(Sequential T-test Analysis of Regime Shifts, STARS），进行补充量模态分析^[11]。当模态改变时，可通过模态跃变指数（Regime Shift Index, RSI）反映，其计算公式为

(1)$dif f = t\sqrt {\sigma _L^2/L} ,$

(2)${{\rm{RSI}}_{i,j}} = \sum\limits_{i = 1}^{j + m} {\frac{{x_i^*}}{{L{\sigma _L}}}} ,\quad m = {\rm{0}},{\rm{1}}, \cdots ,L - 1,$

式中，diff为两个连续且存在显著差异的格局模态平均值的差值；j为新模态的起始年份；L为所选取的截断长度；t为置信水平p下，2L−2自由度的t分布值；$ {\sigma }_{L}^{2} $为方差；RSI为模态跃变指数；当模态向上转变，则$ {x}_{i}^{\mathrm{*}}={{x}}_{{i}}{-}{{\overline{x}}_{{R}2}'} $，当模态向下转变，则$ {x}_{i}^{\mathrm{*}}={{\overline{x}}_{{R}2}'}{-{x}}_{{i}} $，则$ {{\overline{x}}_{{R}2}'}= $ $ {{\overline{x}}_{{R}1}}\pm dif f $；x_i第i个变量；$ {{\overline{x}}_{{R}1}} $为L年内变量x的平均值；${{\overline{x}}_{{R}2}'} $为新模态中变量x的平均值。根据相关文献和参数计算设定标准^[11-12]，设置显著性检验置信水平p等于0.1，截断长度L等于20，huber权重参数H等于2。

2.3 补充量模型构建

选用Ricker模型构建基于环境因子的亲体−补充量模型^[13]。其公式为

(3)${{R}} = {{\alpha S}}{{\rm{e}}^{{{\rm{\beta }}_0}{{s}} + \sum\limits_{i = 1}^n {{{\rm{\beta }}_i}} {X_i}}},$

模型两边取对数后可转化为线性形式：

(4)$\ln (R/S) = \ln {\rm{\alpha }} + {{\rm{\beta }}_0}{{S}} + \sum\limits_{i = 1}^n {{{\rm{\beta }}_i}} {X_i},$

式中，${R}/{S} $为繁殖率，又称为再生长率；α、β、β₀为模型中的待定参数；R为补充量；S为亲体量；X表示影响因子；n为影响因子的数量。本次研究中变量因子包括亲体量、SST、SSS、SSH、ONI和PDO，共6个变量，可构建k=2⁶个模型，所有模型构成模型集合M={M₁，M₂，···，M_k}。

2.4 贝叶斯模型平均法

2.4.1 模型参数估计

贝叶斯模型平均法（Bayesian Model Averaging，BMA）的基本原理为，根据已有的主观信息指定先验分布，解释变量的后验包含概率（Posterior Inclusion Probabilities，PIP）大小可作为选择解释变量的客观标准，以模型后验概率（Posterior Model Probabilities，PMP）为权重对2.3节中所有单项模型进行加权平均，形成包含所有信息的集合M，用于被解释变量的预测^[14]。

在BMA分析框架下，参数$ \beta $不再是固定常数，其与模型M一样服从某一特定的统计分布，所以得到的参数$ \beta $的后验均值也是加权平均的形式。统计分析使用基于R中的BMS包执行^[15-16]。

为了不对模型的优劣加入主观色彩，模型先验设为均匀先验分布（即相同的模型先验概率），参数先验选择单位信息先验。利用马尔科夫链蒙特卡洛（Markov Chain Monte Carb，MCMC）方法在模型空间中抽样，计算迭代10 000次，抛弃前2 000次的计算结果。

模型后验概率PMP表达式为

(5)$\begin{split} {{p}}\left( {{M_\gamma }\big|y,X} \right)& = \frac{{{{p}}\left( {y\big|{M_\gamma },X} \right)p\left( {{M_\gamma }} \right)}}{{{{p}}\left( {y\big|X} \right)}} \\& = \frac{{{{p}}\left( {y\big|{M_\gamma },X} \right)p\left( {{M_\gamma }} \right)}}{{\displaystyle\sum\limits_{i = 1}^k {{p}} \left( {y\big|{M_i},X} \right)p\left( {{M_i}} \right)}},\end{split}$

各因子的后验包含概率PIP的表达式为

(6)${{p}}\left( {\beta \big|y,X} \right) = \sum\limits_{i = 1}^k {p\left( {{M_i}\big|y,X} \right)} ,$

参数$\, \beta $的后验分布为

(7)${{p}}\left( {\beta \big|y,X} \right) = \sum\limits_{i = 1}^k {p\left( {\beta \big|{M_i},y,X} \right)p\left( {{M_i}\big|y,X} \right)} ,$

式中，y为ln(R/S)；X表示影响因子；$ {M}_{\gamma }\in M $；$ p\left({M}_{\gamma }\right) $为模型$ {M}_{\gamma } $的先验分布概率；$\displaystyle\sum\limits_{i = 1}^k {{p}} \left( {y\big|{M_i},X} \right)p\left( {{M_i}} \right)$为观测数据的概率总和。

2.4.2 补充量预测

贝叶斯模型平均法不仅可用于解释变量的选择，还可用于预测被解释变量。不同模态中，不同环境因子以不同的组合形式拟合多个亲体补充量模型，以PMP为权重对所有可能的模型进行加权平均，通过BMA综合多个模型进行集合预报，求得最后两年的预报结果^[17]。为便于考虑模态变动对亲体−补充量的影响，将根据模态分析结果，设计多个基于BMA的亲体−补充量预测模型用于对比分析。

(8)${y_t} = \sum\limits_{i = 1}^k {{y_{t,i}}} {{p}}\left( {{M_i}\big|y,X} \right),$

式中，$ {y}_{t} $为t年ln(R/S)预测值；$ {y}_{t,i} $为t年第i个模型的ln(R/S)预测值；$ {p}\left({M}_{i}\big|y, X\right) $为第i个模型的后验概率。

3 结果

收起

3.1 模态检验

1971−2017年，智利竹筴鱼补充量发生了显著的年际变化，大体上呈现先上升再降低，再上升降低的连续波动过程，具体表现为补充量在1971−1985年处于波动上升阶段，且在1985年达到最高值，在随后的几年里其补充量持续下降，1995年后再次上升，于2000年达到第2个高峰值，随后补充量下降，在2002−2015年维持在较低的范围内（图2）。STARS分析结果表明，1971−2017年间智利竹筴鱼资源补充量可分为5个模态：第1模态为1971−1980年，该模态中补充量处于上升阶段，补充量平均水平相对较低，多个国家和地区对该渔业进行了开发和利用；第2模态为1981−1990年，该模态中的补充量平均水平最高；第3模态为1991−2001年，该阶段补充量明显的起伏变化，但其平均值仍保持在较高的水平；第4模态为2002−2015年，该模态中的补充量平均水平最低；第5模态为2016−2017年，为新的模态阶段，与前一个模态相比，该阶段补充量水平有明显的上升，但该模态时间序列短，故在后续分析中不予讨论（图2，图3）。

3.2 气候变化特征

由图4可知，在第1模态中，厄尔尼诺累计值持续下降。在第2模态中，正负序列（不同气候事件）交替的频次高，变动趋势不稳定。在第3模态中，连续4 a为正数序列，厄尔尼诺累计值持续上升，经历连续两年减少，1997年迎来强厄尔尼诺，厄尔尼诺指数高达1.38，其后4 a为负数序列，累计值持续下降。在该模态中，厄尔尼诺指数累计值有明显的上升下降过程，变动趋势明显。在第4模态中，2002−2006年间整体趋势上升，2007−2014年整体趋势下降，但连续增长或下降的过程没有第3模态明显。分析发现，发生模态变化的年份都为厄尔尼诺年（图4）。

由图5可知，在第1模态中，PDO指数累计值持续下降。在第2模态中，1981−1988年PDO指数累计值持续上升，其后两年持续下降。在第3模态中，1991年PDO指数骤降至−1.78，整体趋势变动不明显，但增减变化的频率更高。在第4模态中，2003年PDO指数上升至1.86，其后小幅增长，2008−2014年整体趋势下降，2015年又上升至2.25。PDO指数的累计值在3个模态中都经历了增长再下降的过程，但在第2模态、第4模态中该过程更明显、更规律，在第3模态中增减变化的频率更高。分析发现，发生模态变动的年份前后，PDO指数都有明显的变化（图5）。

3.3 不同模态中影响因子的筛选分析

在不考虑模态变动的亲体−补充量关系中，亲体量的PIP为99%，在环境因子中，SSH的PIP最高，为76%。第1模态中，亲体量的PIP占34%，PIP最高的环境因子为SSS，其PIP值为46%；第2模态中，亲体量的PIP为86%，El Niño的PIP值为49%，为该模态中PIP值最高的环境因子；第3模态中，亲体量的PIP为100%，PIP最高的环境因子为PDO，其值为41%；第4模态中，亲体量的PIP为98%，PIP最高的环境因子为El Niño，其值为33%（表1）。

3.4 预测模型

在不考虑模态变动的亲体−补充模型中，6个因子各自组合形成多个模型，表2列举PMP最高的前5个模型。模型1为基于SSH的亲体−补充量模型，后验概率最高，为36%。5个模型集合预报的累计后验概率达到了70%。

考虑模态变动对亲体−补充量的影响，依次去除第1模态、第2模态、第3模态，设计3个基于BMA的亲体补充量模型(表3)：（1）方案a，由第2模态、第3模态、第4模态数据（1981−2015年）构建亲体补充量模型；（2）方案b，基于第3模态、第4模态（1991−2015年）构建亲体补充量模型；（3）方案c，由第4模态（2002−2015年）构建亲体补充量模型。3个方案中，方案a中最佳模型的后验概率最高为24%，各个方案的累计后验概率基本相同(表3)。

由图6可知，从预测结果看，BMA可通过概率分布定量评价预报的不确定性（图6）。不区分模态的预测结果最接近真实值，优于方案a、b、c，方案c最为保守，预测值明显小于其他方案。但其他方案的预测概率约为60%～80%，方案c的预测概率大于80%，优于其他方案。

4 分析与讨论

收起

4.1 气候变化与补充量模态变动之间的联系

南太平洋智利竹筴鱼分布极其广泛，东起智利、秘鲁和厄瓜多尔沿岸，西至新西兰、澳大利亚沿海，在30°～60°S、78°～177°W广阔海域均有密集的智利竹筴鱼群体，该区域也称为“竹筴鱼带”^[18]。该海域大尺度环流、中尺度海流以及各种海洋环境要素的变化活跃，智利竹筴鱼的补充量受到环境要素、气候事件的影响明显^[19-20]。变动的生态环境也可能会引起竹筴鱼补充量的模态变动，从而给资源评估和管理带来不可忽略的影响^[21]。

本研究基于STARS对智利竹筴鱼补充量进行模态分析，结果表明其存在多个模态，且在不同模态中，占主导地位的环境气候因子并不相同。在第1模态中，亲体量对补充量变动的解释概率明显偏低，各影响因子的解释概率也明显存在异常，可能的原因是该阶段渔业资源变化更多的受到捕捞因素的影响^[22]，并不适用于分析补充量与环境的关系。其余模态与传统环境因子分析（不区分模态变动）相比，在环境因子的解释概率排序上有明显区别。在传统环境因子分析中，SSH占有极高的解释概率，而考虑模态变动的分析结果认为，各环境因子在不同的时期存在不同的排序。从统计结果分析，SSH的PIP明显偏大且优于其他环境因子，极少有环境因子的解释概率超过60%^[23]，从目前已有的智利竹筴鱼与环境关系的研究结果分析，也更偏向于区分模态的分析结果^[24-25]。传统环境因子分析往往将长时间序列的资源变动归因于一个关键因子上，而忽略了资源变动受到多个环境因子的共同影响，在不同时期占主导地位的环境因子也会发生改变。要深入了解渔业资源对环境气候的响应以及鱼类动态变动及生态系统调节机制等内容，模态分析是不可忽略的。

在不同环境时期，智利竹筴鱼的补充量水平会发生改变，基于这些不同时期的平均补充量水平来拟合亲体−补充量关系和制定管理方案，可能不利于渔业的可持续发展^[26]。在智利竹筴鱼的评估报告^[1]中，选取两个基于不同时间段亲体−补充量关系（1970−2015年、2000−2015年），并分别放入渔业资源评估模型，其结果存在较大的区别。研究认为，基于2000−2015年数据的分析结果更为保守，可用于最终的资源评估与管理。此外，SPRFMO评估报告^[1]中选取的时间段（2000−2015年）与本研究中的第4模态相近，基于第4模态的预测结果也明显比其余方案更保守（图6）。但从预测结果看，4种预测方案中，不考虑模态变动的方案最接近真实值，基于第4模态的结果偏差最大，其原因可能是用于预测的2016年和2017年补充量水平高于第4模态补充量水平，与2015年以前有较大的差异（图2），因此只基于第4模态（低水平）进行预测其结果自然会偏小。目前，在SPRFMO智利竹筴鱼评估报告中，数据段的选取更多的依赖于人为经验，但本文基于模态检验的科学方法将智利竹筴鱼划分为多个模态进行分析，能为数据段的选取及资源评估提供了科学依据。

4.2 厄尔尼诺和PDO事件对智利竹筴鱼补充量的影响

在第2模态、第4模态中，厄尔尼诺皆为解释概率最高的环境因子，其指数的正负序列（不同气候事件）交替频次高，持续处于某一气候事件的时间短，与第3模态中的表现有明显的区别。这可能是由于不同气候事件高频的交替发生，这种不稳定的气候环境对智利竹筴鱼资源补充量产生了更大的影响。研究认为，智利竹筴鱼在面对不同气候事件时，存在不同的应对机制。以厄尔尼诺和拉尼娜为例，在厄尔尼诺期间，海表面温度、海表面高度、风速、盐度等众多环境要素都会发生改变，智利竹筴鱼仔稚鱼分布、洄游路径等多方面产生变化^[27]。同时，由于捕食者减少、不同性质海水扩展等原因，智利竹筴鱼栖息范围会有所扩展^[9]。Arcos等^[28]研究认为，厄尔尼诺对智利竹筴鱼仔稚鱼肥育场的影响往往会持续3～4 a的时间。在拉尼娜期间，气候条件与厄尔尼诺存在显著不同，海水的降温可能导致智利附近海域的幼鱼进入厄瓜多尔和智利的专属经济区^[29]。

PDO是一种反映从年际到年代际的气候变率强信号，是太平洋重要的气候系统之一，对太平洋乃至全球的大气环流有着重要影响，也是气候预测的关键因子^[30-31]。除了对各地区气候的影响，PDO对海洋鱼类种群的兴衰也起着调节作用^[32]。PDO与许多中上层鱼类的补充量有着直接的联系，PDO正负相位时期、冷暖时期对应的补充量都有明显的变化^[33]。Dioses等^[34]研究认为，年际间PDO变化与智利竹筴鱼的纬度隔离有一定关系。分析认为，1977−1999年为PDO暖期，1999−2015年为PDO冷期^[35]。冷暖期交替发生在本研究中的第3模态时期，这也解释了该模态中PDO年际间变动更为不稳定，而不稳定的年际间PDO变化对智利竹筴鱼补充量产生了更大的影响。

本研究认为，在发生模态变动年份前后，均发生了厄尔尼诺事件且PDO指数有明显的跃变，这表明智利竹筴鱼补充量与厄尔尼诺和PDO气候事件存在紧密联系，PDO年代际冷暖期交替与厄尔尼诺现象等可能是诱发智利竹筴鱼补充量发生模态转变的重要因素。Espino^[36]研究也表明，PDO、厄尔尼诺、智利竹筴鱼三者间的变动存在紧密联系。一方面，不同尺度间的气候变化存在交互作用，年代际PDO变化会影响厄尔尼诺现象，而物理环境中这些低频的变化可能会导致海洋生态系统的改变^[37-38]，从而导致智利竹筴鱼补充量的模态改变。另一方面，PDO和厄尔尼诺会以不同的方式和不同的量级对渔业资源造成影响^[14]。在今后的研究中，需要结合物理海洋、气候学等多学科，加强了解智利竹筴鱼补充量对气候环境变化的响应机制和机理。

4.3 其他环境因素对智利竹筴鱼补充量的影响

海表面高度是用于检测和表征中尺度涡流的重要物理量^[39-40]。涡流系统对智利海域叶绿素的聚集起着重要作用，冬季高浓度叶绿素海域多数是因此形成^[41]。同时，涡流和上升流之间的耦合使得高初级生产力区域扩张，持续时间延长，浮游生物随之聚集，为智利竹筴鱼仔鱼提供连续的食物来源，有利于其生长^[42-43]。此外，涡流的垂直影响深度完全覆盖仔鱼生活水层，部分仔鱼被夹带入涡流内会被运到外海，在适当条件下可能会促进产生新的仔、稚鱼栖息地^[44-45]。在早期生活史阶段，涡流等中尺度海流变化过程对智利竹筴鱼空间分布、生长以及资源补充量等有着不可忽略的影响。

在本研究中，温盐的PIP基本维持在20%～30%，这表明两者对智利竹筴鱼补充量也具有一定的影响。东南太平洋海域拥有众多不同温盐特质的水团，占主导地位的水团改变及强弱改变都会导致竹筴鱼分布情况的改变^[46]。当高温（温度大于18.5℃）、高盐（盐度大于34.9）亚热带表层水团南移并覆盖近岸海域时，鱼群分布更为集中，渔获量增加；当低温（约12℃）、低盐（约34.25）的亚南极中层水团强于赤道表层水时，智利竹筴鱼分布在更深的水层中，渔获量明显减少。秘鲁海域的智利竹筴鱼常聚集在表温为11～28℃的水域，智利和公海海域智利竹筴鱼聚集在表温为9～20℃的水域, 当海表温上升时，智利竹筴鱼适宜栖息地有明显的南移趋势。与其他鱼类不同，盐度对智利竹筴鱼的影响更为直接^[47]，张敏等^[29]研究认为，较稳定的盐度将有利于智利竹筴鱼集群的稳定。

模态检验方法的参数设置不仅需要遵循模型使用的前提条件，同时还需考虑智利竹筴鱼渔业资源本身变动规律及其特性，任意降低截断时间长度可能会出现更多的模态^[48]，时间序列过短的模态不利于渔业资源评估和管理，选择一个恰当合理的时间序列长度在不同模态的渔业资源评估是极为重要的。气候环境变化对智利竹筴鱼早期生活史阶段的影响是不可忽略的，合理的划分模态，有助于探究不同阶段的环境与补充量的影响模式。其次，利用模态变动的方法将长期气候变化对补充量的影响纳入考虑，有助于提高对鱼类动态及其生态系统调节机制的理解，制定更合理的管理策略^[49]。同时，贝叶斯模型平均法将补充量影响机制归因于多个环境因子共同作用的结果，分析每个因子不同的解释能力，还能以概率分布的形式描述补充量预报的准确性。模态分析与贝叶斯模型平均法的结合适用于多信息集合，可在今后的研究中对更多的环境气候因素进行分析，促进物理海洋、气候、生物等多学科交叉，建立起科学的智利竹筴鱼补充量预报模型，为后续的渔业资源评估和管理提供基础。

基金

收起

国家重点研发计划（2019YFD0901404）；上海市科技创新行动计划（19DZ1207500）

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

South Pacific Regional Fisheries Management Organisation. 7th Scientific Committee Meeting Report[R]. Cuba: SPRFMO, 2019.

Kuparinen A, Keith D M, Hutchings J A. Increased environmentally driven recruitment variability decreases resilience to fishing and increases uncertainty of recovery[J]. ICES Journal of Marine Science, 2014, 71(6): 1507−1514.

徐冰, 陈新军, 陆化杰, 等. 秘鲁外海茎柔鱼资源丰度和补充量与海表温度的相关关系[J]. 海洋渔业, 2013, 35(3): 296−302.

Xu Bing, Chen Xinjun, Lu Huajie, et al. Relationship between sea surface temperature and abundance index, recruitment of Dosidicus gigas in the waters off Peru[J]. Marine Fisheries, 2013, 35(3): 296−302.

Steele J H. Regime shifts in fisheries management[J]. Fisheries Research, 1996, 25(1): 19−23.

Beamish R J, Noakes D J, McFarlane G A, et al. The regime concept and natural trends in the production of Pacific salmon[J]. Canadian Journal of Fisheries and Aquatic Sciences, 1999, 56(3): 516−526.

Brodziak J, Legault C M. Model averaging to estimate rebuilding targets for overfished stocks[J]. Canadian Journal of Fisheries and Aquatic Sciences, 2005, 62(3): 544−562.

王亮. 贝叶斯模型平均方法研究综述与展望[J]. 技术经济与管理研究, 2016(3): 19−23.

Wang Liang. Overview and prospect of Bayesian model averaging[J]. Technoeconomics & Management Research, 2016(3): 19−23.

South Pacific Regional Fisheries Management Organisation. 6th Scientific Committee Meeting Report[R]. Cuba: SPRFMO, 2018.

Hintzen N T, Corten A, Gerlotto F, et al. Hydrography and Jack Mackerel stock in the South Pacific[R]. Netherlands: IMARES, 2013.

Muck P, Castillo O S, Carrasco S. Abundance of sardine mackerel and horse mackerel eggs and larvae and their relationship to temperature, turbulence and anchoveta biomass off Peru. In the Peruvian anchoveta and its upwelling ecosystem: three decades of change[M]. Peru: ICLARM, 1990: 268−275.

Rodionov S, Overland J E. Application of a sequential regime shift detection method to the Bering Sea ecosystem[J]. ICES Journal of Marine Science, 2005, 62(3): 328−332.

袁兴伟, 刘尊雷, 程家骅, 等. 气候变化对冬季东海外海中下层游泳动物群落结构及重要经济种类的影响[J]. 生态学报, 2017, 37(8): 2796−2808.

Yuan Xingwei, Liu Zunlei, Chen Jiahua, et al. Impact of climate change on nekton community structure and some commercial species in the offshore area of the northern East China Sea in winter[J]. Acta Ecologica Sinica, 2017, 37(8): 2796−2808.

詹秉义. 渔业资源评估[M]. 北京: 中国农业出版社, 1995, 291−298.

Zhan Bingyi. Fish Stock Assessment[M]. Beijing: China Agriculture Press, 1995: 291−298.

Amini S, Parmeter C F. Bayesian model averaging in R[J]. Journal of Economic and Social Measurement, 2011, 36(4): 253−287.

Draper D. Assessment and propagation of model uncertainty[J]. Journal of the Royal Statistical Society: Series B (Methodological), 1995, 57(1): 45−70.

Stefan Z, Martin F. Bayesian model averaging employing fixed and flexible priors: The BMS package for R[J]. Journal of Statistical Software, 2015, 68(4): 1−37.

刘尊雷, 袁兴伟, 杨林林, 等. 亲体量和环境对东海小黄鱼补充成功率的影响[J]. 应用生态学报, 2015, 26(2): 588−600.

Liu Zunlei, Yuan Xingwei, Yang Linlin, et al. Effect of stock abundance and environmental factors on the recruitment success of small yellow croaker in the East China Sea[J]. Chinese Journal of Applied Ecology, 2015, 26(2): 588−600.

张敏, 邹晓荣. 大洋性竹筴鱼渔业[M]. 北京: 中国农业出版社, 2011.

Zhang Min, Zou Xiaorong. Oceanic Jack Mackerel Trachurus murphyi Fishery[M]. Beijing: China Agriculture Press, 2011.

Pineda J, Hare J A, Sponaugle S. Larval transport and dispersal in the coastal ocean and consequences for population connectivity[J]. Oceanography, 2007, 20(3): 22−39.

Kasai A, Komatsu K, Sassa C, et al. Transport and survival processes of eggs and larvae of jack mackerel Trachurus japonicus in the East China Sea[J]. Fisheries Science, 2008, 74(1): 8−18.

Ianelli J N, Hollowed A B, Haynie A C, et al. Evaluating management strategies for eastern Bering Sea walleye pollock (Theragra chalcogramma) in a changing environment[J]. ICES Journal of Marine Science, 2011, 68(6): 1297−1304.

South Pacific Regional Fisheries Management Organization. Report of the science working group[R]. Cuba: SPRFMO, 2011.

张畅, 陈新军. 海洋环境因子对澳洲鲐亲体补充量关系的影响—基于贝叶斯模型平均法的研究[J]. 海洋学报, 2019, 41(2): 99−106.

Zhang Chang, Chen Xinjun. The impact of environmental factors on stock-recruitment relationship of spotted mackerel—Based on Bayesian model averaging method[J]. Haiyang Xuebao, 2019, 41(2): 99−106.

Li Gang, Zou Xiaorong, Chen Xinjun, et al. Standardization of CPUE for Chilean jack mackerel (Trachurus murphyi) from Chinese trawl fleets in the high seas of the Southeast Pacific Ocean[J]. Journal of Ocean University of China, 2013, 12(3): 441−451.

常永波. 几种人工神经网络模型在智利竹筴鱼渔场渔情预报中的比较研究[D]. 上海: 上海海洋大学, 2016.

Chang Yongbo. Fishery forecasting comparative study of Trachurus murphyi in the southeast pacific fishing grounds based on artificial neural network model[D]. Shanghai: Shanghai Ocean University, 2016.

Haltuch M A, Punt A E, Dorn M W. Evaluating the estimation of fishery management reference points in a variable environment[J]. Fisheries Research, 2009, 100(1): 42−56.

Blanco J L, Thomas A C, Carr M E, et al. Seasonal climatology of hydrographic conditions in the upwelling region off northern Chile[J]. Journal of Geophysical Research: Oceans, 2001, 106(C6): 11451−11467.

Arcos D F, Cubillos L A, Núñez S P. The jack mackerel fishery and El Niño 1997–98 effects off Chile[J]. Progress in Oceanography, 2001, 49(1/4): 597−617.

张敏, 邹晓荣, 季星辉, 等. 东南太平洋公海水域智利竹筴鱼探捕及其商业开发前景探讨[J]. 水产学报, 2005, 29(3): 386−391.

Zhang Min, Zou Xiaorong, Ji Xinghui, et al. Discussion on exploratory fishing of Chilean jack mackerel (Trachurus murphyi) in the open sea of the southeast Pacific and prospect of its commercial exploitation[J]. Journal of Fisheries of China, 2005, 29(3): 386−391.

吕俊梅, 琚建华, 张庆云, 等. 太平洋年代际振荡冷、暖背景下ENSO循环的特征[J]. 气候与环境研究, 2005, 10(2): 238−249.

Lü Junmei, Ju Jianhua, Zhang Qingyun, et al. The characteristics of ENSO cycle in different phases of Pacific decadal oscillation[J]. Climatic and Environmental Research, 2005, 10(2): 238−249.

杨修群, 朱益民, 谢倩, 等. 太平洋年代际振荡的研究进展[J]. 大气科学, 2004, 28(6): 979−992.

Yang Xiuqun, Zhu Yimin, Xie Qian, et al. Advances in studies of Pacific decadal oscillation[J]. Chinese Journal of Atmospheric Sciences, 2004, 28(6): 979−992.

Mantua N J, Hare S R. The Pacific decadal oscillation[J]. Journal of Oceanography, 2002, 58(1): 35−44.

Zwolinski J P, Demer D A. Environmental and parental control of Pacific sardine (Sardinops sagax) recruitment[J]. Ices Journal of Marine Science, 2014, 71(8): 2198−2207.

Dioses T. Abundance and distribution patterns of Jack mackerel Trachurus murphyi in Peru[J]. Revista Peruana de Biologia, 2013, 20(1): 67−74.

Yu Wei, Wen Jian, Chen Xinjun, et al. Trans-Pacific multidecadal changes of habitat patterns of two squid species[J]. Fisheries Research, 2021, 233: 105762.

Espino M. The Jack mackerel Trachurus murphyi and the environmental macro-scale variables[J]. Revista Peruana De Biología, 2013, 20(1): 9−20.

Wang Chunzai, Fiedler P C. ENSO variability and the eastern tropical Pacific: A review[J]. Progress in Oceanography, 2006, 69(2/4): 239−266.

Lorenzo E D, Ohman M D. A double-integration hypothesis to explain ocean ecosystem response to climate forcing[J]. Proceedings of the National Academy of Sciences of the United States of America, 2013, 110(7): 2496−2499.

Chaigneau A, Gizolme A, Grados C. Mesoscale eddies off Peru in altimeter records: Identification algorithms and eddy spatio-temporal patterns[J]. Progress in Oceanography, 2008, 79(2/4): 106−119.

Wang Xin, Du Yunyan, Zhou Chenghu, et al. An improved, SSH-based method to automatically identify mesoscale eddies in the ocean[J]. Journal of Tropical Oceanography, 2013, 32(2): 15−23.

Correa-Ramirez M A, Hormazábal S, Yuras G. Mesoscale eddies and high chlorophyll concentrations off central Chile (29°–39°S)[J]. Geophysical Research Letters, 2007, 34(12): L12604.

Govoni J J, Hare J A, Davenport E D, et al. Mesoscale, cyclonic eddies as larval fish habitat along the southeast United States shelf: a Lagrangian description of the zooplankton community[J]. ICES Journal of Marine Science, 2010, 67(3): 403−411.

Morales C E, Torreblanca M L, Hormazabal S, et al. Mesoscale structure of copepod assemblages in the coastal transition zone and oceanic waters off central-southern Chile[J]. Progress in Oceanography, 2010, 84(3/4): 158−173.

Parada C, Núñez S, Correa M, et al. Advances in biophysical modeling of Chilean jack mackerel in the South Pacific[J]. ICES CM, 2010, 20: 32.

Parada C, Mullon C, Roy C, et al. Does vertical migratory behaviour retain fish larvae onshore in upwelling ecosystems? A modelling study of anchovy in the southern Benguela[J]. African Journal of Marine Science, 2008, 30(3): 437−452.

Serra R. Important life history aspects of the Chilean jack mackerel, Trachurus symmetricus murphyi[J]. Investigacion Pesquera (Chile), 1991, 36: 67−83.

陈春光. 东南太平洋智利竹筴鱼渔场空间分布及其栖息地适宜性指数的研究[D]. 上海: 上海海洋大学, 2015.

Chen Chunguang. Study on fishing ground spatial distribution and habitat suitability index (HSI) model of Chilean Jack Mackerel (Trachurus murphyi) in the Southeast Pacific Ocean[D]. Shanghai: Shanghai Ocean University, 2015.

Rodionov S N. A sequential algorithm for testing climate regime shifts[J]. Geophysical Research Letters, 2004, 31(9): L09204.

Wayte S E. Management implications of including a climate-induced recruitment shift in the stock assessment for jackass morwong (Nemadactylus macropterus) in south-eastern Australia[J]. Fisheries Research, 2013, 142: 47−55.

2021年第43卷第9期

PDF下载

128

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2021130

接收时间：2021-01-28
首发时间：2026-02-26
出版时间：2021-09-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2021-01-28
修回日期：2021-04-29

基金

国家重点研发计划（2019YFD0901404）；上海市科技创新行动计划（19DZ1207500）

作者信息

¹上海海洋大学海洋科学学院，上海 201306

²青岛海洋科学与技术试点国家实验室海洋渔业科学与食物产出过程功能实验室，山东青岛 266237

³上海海洋大学大洋渔业资源可持续开发教育部重点实验室，上海 201306

⁴上海海洋大学国家远洋渔业工程技术研究中心，上海 201306

⁵上海海洋大学农业农村部大洋渔业开发重点实验室，上海 201306

通讯作者:

陈新军（1967－），男，教授，浙江省义乌市人，从事渔业资源与渔场学等研究。E-mail：xjchen@shou.edu.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.12284/hyxb2021130

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT