Study of impacts of environmental factors on stock and recruitment relationship of the anchovy stock in the Yellow Sea
-
摘要:
鱼类年际资源量的波动可以归因于年间环境条件的变化和该种鱼类亲体数量的变化。文章根据1990~2001年间黄海中南部鳀鱼声学调查评估结果，以及黄海千里岩海区在此年间的表层水温和营养盐统计数据，以Ricker模型(R=αSe-βS)为基础对黄海鳀鱼(Engraulis japonicas)亲体-补充量关系进行了初步研究。结果表明，黄海千里岩水域表层水温、磷酸盐浓度等环境条件因素对补充量有重要影响。
Abstract:Changes in fish year-class strength have been attributed to year-to-year variability in environmental conditions, such as sea surface temperature (SST) and phosphate concentration, and spawning stock biomass(SSB).In this study, we examined the relationship between environmental factors, SSB and recruitment for the anchovy (Engraulis japonicas) stock in the Yellow Sea using a traditional stock-recruitment Ricker model. The recruitment fluctuation of the anchovy stock in the Yellow Sea was successfully reproduced by the Ricker stock-recruitment model with environmental factors. The findings suggest that the impacts of SST and phosphate concentration observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea on the anchovy stock in the Yellow Sea are very important.
-
Keywords:
- stock-recruitment model /
- environmental condition index /
- AIC /
- BIC
-
鳀鱼(Engraulis japonicas)属鳀科(Engranlidae)、鳀属(Engraulis)，是一种小型中上层游泳鱼类[1]。广泛分布于北起库页岛，南至台湾海峡的太平洋西北部水域。我国主要分布在渤海、黄海和东海海域[2]。作为一种小型中上层鱼类，鳀鱼的数量变动容易受到各种环境因素的影响。
亲体-补充量关系(stock and recruitment relationship，SRR)是渔业资源评估和管理的基础[3-4]之一。依据SRR模型也可以计算和确定许多重要的渔业管理目标值，如最大持续产量(maximal sustainable yield，MSY)，和维持MSY所需的亲体量等。然而，由于补充量的剧烈波动、亲体量的测量误差以及时间序列的相关性和空间分布上的变动性等常常使参数的估计及模型的选择存在很大的不确定性[4]。SRR的估计是渔业资源评估中最困难的任务之一。
渔业生物的亲体数量对种群的补充量大小具有重要的影响，BEVERTON-HOLT[5]和RICKER[6]先后确立了SRR，建立了相应的鱼类繁殖模型，这些模型长期为渔业生物学家在不同类型的渔业种群中应用[7]。此文选用对环境因子响应比较显著的Ricker模型对资源数据进行拟合[8]。
在以前对黄海鳀鱼的亲体-补充量(stock and recruitment，SR)研究当中，研究者都将各种环境因素综合到一个误差项中进行研究。此文尝试将各个环境因子独立出来，再通过多元回归方法进行分析，探讨其对于SRR的影响程度。
1. 材料与方法
1.1 数据来源
1.1.1 黄海鳀鱼SR数据
以ZHAO等[8]研究中1987~2002年的渔业声学调查数据和渔获量统计数据及所估算的自然死亡系数和成时间序列的SR数据作为此文的数据来源(图 1)。选取1990~2001年间的12组SR数据作为观测数据。
![]()
图 1 1987~2001年黄海鳀鱼亲体-补充量数据散点图图中连线连接相邻年份的补充数据Figure 1. Scatter plot of the stock and recruitment data for the anchovy (E.japonicas) stock in the Yellow SeaThe line joins recruitment values with the year indicated.以1龄鱼的个体数量作为补充量的量度指标。用前一年度冬季鳀鱼声学调查测得的越冬群体生物量与前一年度鳀鱼1/2的年渔获量之差作为亲体量的度量指标。其中前一年度1月份鳀鱼越冬群体的生物量则是在相应生物学评估结果基础上，经过以VPA返算的补充量对其1龄鳀鱼的生物量进行调整而得[8]。假设雌雄个体比例为1 : 1，每个成熟个体的生殖力与体重呈简单的正比函数关系。
1.1.2 水温、营养盐数据
黄海水文状况以国家海洋局黄海千里岩海洋环境监测站的观测数据为指标。采用1990~2001年间黄海千里岩海区的表层水温(sea surface temperature，SST)(图 2)和磷酸盐浓度(phosphate concentration，PC)数据(图 3)作为此研究的原始数据。
![]()
图 2 黄海千里岩1990~2001年4~7月份平均日表层水温数据(折线图)与黄海鳀鱼补充量数据(柱状图)变化趋势图Figure 2. Variations of mean SST from April to July (line) and recruitment (bar) observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea in the period of 1990~2001![]()
图 3 黄海千里岩1990~2001年平均磷酸盐浓度数据(折线图)与黄海鳀鱼补充量数据(柱状图)变化趋势图Figure 3. Variations of mean phosphates concentration (line) and recruitment (bar) observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea in the period of 1990~20011.2 方法
1.2.1 SRR模型
选用通常使用的研究亲体数量(spawning stock biomass，SSB)和补充量(recruitment，R)之间关系的Ricker模型[8]，并使用多元回归的分析方法筛选了黄海千里岩海域逐月的SST、年平均SST、年平均PC等环境数据。
Ricker的SRR模型表达式为：R=αSe-βS
其中α、β是模型中的待定参数。
用环境条件指数$\alpha_i=\alpha_0+\sum\limits_{i=1}^n \alpha_i X_i$取代模型中的参数α，Ricker模型的表达形式变为：
$$ R=\left(\alpha_0+\sum\limits_{i=1}^n \alpha_i X_i\right) S_t \exp \left(-\beta S_t\right) $$ α0，α1，……，αi为环境条件指数回归方程中的参数，X1，……，Xi的含义如表 1。
表 1 黄海千里岩水域环境因子与黄海鳀鱼补充量数据的相关系数Table 1. The correlation of environmental factors and recruitment data observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea环境因子
environmental factorsXi 相关系数
correlation coefficient年平均磷酸盐浓度
mean phosphates concentrationX1 0.4091 平均表层水温
average sea surface temperature1月 X2 0.2656 2月 X3 0.3707 3月 X4 -0.0761 4月 X5 -0.6213 5月 X6 -0.6614 6月 X7 -0.4987 7月 X8 -0.5290 8月 X9 -0.0451 9月 X10 -0.4223 10月 X11 -0.3174 11月 X12 -0.1593 12月 X13 -0.1838 全年 X14 -0.4679 4~7月 X15 -0.8096 1.2.2 AIC(akaike information criterion)与BIC(bayesian information criterion)
AIC与BIC在统计学上，作为模型选择的标准已经被广泛应用到很多领域[9]。此文采用AIC和BIC标准对SRR模型进行比较。
AIC=-2ln (L)+2m
BIC=-2ln(L)+mln(n)
其中m是模型中参数的个数，n是观测数据的个数，L表示最大似然值。取得最小AIC或BIC的模型被选为最适模型。
2. 结果
根据图 1所示的1990~2001年黄海鳀鱼亲体数量(×106 t)与补充量(×1011尾)资料，用Ricker模型进行拟合，求得各项参数值为α=1.5931，β=0.3058。AIC与BIC值分别为37.2062，38.176。
表 1显示了黄海鳀鱼补充量与黄海千里岩水域逐月平均、全年平均、4~7月平均SST及PC数据的相关系数。
这里，设$\alpha_t=\frac{R_t}{S_t} \exp \left(\beta S_t\right)$是一个对黄海鳀鱼补充量有重要综合影响的可变环境因子[10-11]，其与补充量的关系如图 4所示。为了对各项环境数据进行统计计算，设定Xi代表各项环境因子，其对应关系如表 1所示。
![]()
图 4 黄海千里岩环境条件指数αt和黄海鳀鱼补充量指数R变化趋势图Figure 4. Variations of environmental conditions index and the recruitment observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea有$\alpha_i=\alpha_0+\sum\limits_{i=1}^n \alpha_i X_i$[10-11]，根据表 1中显示的相关系数，采用多元线性回归方法拟合对αt有重要影响的环境因子，年平均PC、9月平均和4~7月平均SST，即有αt=α0+α1X1+α10X10+α15X15。求得各项参数值α0=7.974；α1=-0.0144；α10=0.216；α15=-0.713。
Rt=(7.974-0.0144Xt+0.216X10-0.713X15)Stexp (-0.3058St)，AIC与BIC的值分别为9.3159和11.7404。将由此求得的Rt的计算值与补充量R的观测值进行作图对比(图 5)，较之单纯使用Ricker模型，加入环境指数的模型拟合的效果更好。
![]()
图 5 黄海鳀鱼补充量观测值和计算值比较图Figure 5. Comparison of observed and calculated recruitment of anchovy (E.japonicas) stock in the Yellow Sea3. 讨论
3.1 环境因子与补充量的相关性
逐月分析黄海千里岩海域平均月SST与补充量的关系，可以看出，4~7月份的SST对黄海鳀鱼的补充量有着显著的影响，因为黄海鳀鱼性腺发育、生殖和仔稚鱼育肥的主要时期是每年的4~10月份[8]，由表 1也可以看出4~7月的SST与黄海鳀鱼补充量的相关系数最大。
3.2 环境条件指数
从图 4中可以看到，1996年(包括1996年)之前的环境条件指数与补充量的变化趋势以及波动程度都极其相似，1996年之后的环境因子与补充量的变化趋势虽然比较相似，但是变化程度有了明显的不同。由ZHAO等[8]研究的年捕捞量数据可知，1997~2000年间的黄海鳀鱼产量逐年明显高于其它年，捕捞力量的加剧可能导致了这一现象的出现。
黄海鳀鱼SRR的主要影响因素在1997年之后发生了变化，1990~1997年间，黄海鳀鱼补充量在波动中保持相对的稳定，1999年之后黄海鳀鱼资源大幅衰退，补充量随之显著降低。由此可见，只有在亲体量变化水平相当的情况下，环境因子的影响才相对显著。
3.3 模型结果比较分析
表 2是黄海鳀鱼SRR模型参数的估计值，以及AIC、BIC的计算值。从中可以看出，加入环境条件因子的Ricker模型对补充量的估计效果更好。其AIC与BIC的值都远远小于单纯用Ricker模型进行拟合的结果。说明环境条件对黄海鳀鱼SRR的影响是比较明显的，对环境条件因子的研究还是有一定意义的。
表 2 黄海鳀鱼SRR模型参数的估计值，AIC、BIC计算值Table 2. Estimated parameters and AIC, BIC in stock recruitment models (SRR)模型
model参数 estimated parameters AIC BIC α β α0 α1 α10 α15 Ricker模型
Ricker model1.5931 0.3058 - - - - 37.2062 38.176 加入αt的Ricker模型
Ricker model with αt- 0.3058 7.974 -0.0144 0.216 -0.713 9.3159 11.7404 鳀鱼资源的衰退，将对鳀鱼渔业产生巨大的负面影响。从图 1可以看出，黄海鳀鱼种群正在迅速向原点位置衰退，为了保证黄海鳀鱼渔业的可持续发展和利用，预防性管理措施的实施是极为必要的。
-
![]()
图 1 1987~2001年黄海鳀鱼亲体-补充量数据散点图
图中连线连接相邻年份的补充数据
Figure 1. Scatter plot of the stock and recruitment data for the anchovy (E.japonicas) stock in the Yellow Sea
The line joins recruitment values with the year indicated.
![]()
图 2 黄海千里岩1990~2001年4~7月份平均日表层水温数据(折线图)与黄海鳀鱼补充量数据(柱状图)变化趋势图
Figure 2. Variations of mean SST from April to July (line) and recruitment (bar) observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea in the period of 1990~2001
![]()
图 3 黄海千里岩1990~2001年平均磷酸盐浓度数据(折线图)与黄海鳀鱼补充量数据(柱状图)变化趋势图
Figure 3. Variations of mean phosphates concentration (line) and recruitment (bar) observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea in the period of 1990~2001
![]()
图 4 黄海千里岩环境条件指数αt和黄海鳀鱼补充量指数R变化趋势图
Figure 4. Variations of environmental conditions index and the recruitment observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea
![]()
图 5 黄海鳀鱼补充量观测值和计算值比较图
Figure 5. Comparison of observed and calculated recruitment of anchovy (E.japonicas) stock in the Yellow Sea
表 1 黄海千里岩水域环境因子与黄海鳀鱼补充量数据的相关系数
Table 1 The correlation of environmental factors and recruitment data observed at the Qianliyan Marine Environment Monitoring Station in the Yellow Sea
环境因子
environmental factorsXi 相关系数
correlation coefficient年平均磷酸盐浓度
mean phosphates concentrationX1 0.4091 平均表层水温
average sea surface temperature1月 X2 0.2656 2月 X3 0.3707 3月 X4 -0.0761 4月 X5 -0.6213 5月 X6 -0.6614 6月 X7 -0.4987 7月 X8 -0.5290 8月 X9 -0.0451 9月 X10 -0.4223 10月 X11 -0.3174 11月 X12 -0.1593 12月 X13 -0.1838 全年 X14 -0.4679 4~7月 X15 -0.8096 
下载: 导出CSV
表 2 黄海鳀鱼SRR模型参数的估计值，AIC、BIC计算值
Table 2 Estimated parameters and AIC, BIC in stock recruitment models (SRR)
模型
model参数 estimated parameters AIC BIC α β α0 α1 α10 α15 Ricker模型
Ricker model1.5931 0.3058 - - - - 37.2062 38.176 加入αt的Ricker模型
Ricker model with αt- 0.3058 7.974 -0.0144 0.216 -0.713 9.3159 11.7404
下载: 导出CSV
-
[1] 张春霖, 成庆泰, 郑傈珊, 等. 黄渤海鱼类调查报告[M]. 北京: 科学出版社, 1955: 51-52. [2] 朱德山, IVERSEN S A. 黄、东海鳀鱼及其他经济鱼类资源声学评估的调查研究[J]. 海洋水产研究, 1990, 11: 18-31. [3] RICKER W E. Computation and interpretation of biological statistics of fish populations[J]. Bull Fish Res Board Can, 1975, 191(18): 1-382. https://www.semanticscholar.org/paper/Computation-and-interpretation-of-biological-of-Ricker/5a8d0d7094c356e3b851fd66bd929ed0e56aabfd
[4] HILBORN R, WALTERS C J. Quantitative fisheries stock assessment: Choice, dynamics and uncertainty[M]. New York: Chapman & Hall, 1992. doi: 10.1007/978-1-4615-3598-0
[5] BEVERTON R J H, HOLT S J. On the dynamics of exploited fish populations[J]. Fish Inverst, 1957, 19(2): 1-533. doi: 10.1007/bf00044132
[6] RICKER W E. Stock and recruitment[J]. J Fish Res Bd Can, 1954, 108(11): 599-623. doi: 10.1139/F54-039
[7] QUINN Ⅱ T J, DERISO R B. Quantitiative fish dynamics[M]. New York: Oxford University Press, 1999. doi: 10.1093/oso/9780195076318.001.0001
[8] ZHAO X, HAMRE J, LI F, et al. Recruitment, sustainable yield and possible ecological consquences of the sharp decline of the anchovy (Engraulis japonicus) stock in the Yellow Sea in the 1990s[J]. Fish Oceanogr, 2003, 12(4): 495-501. doi: 10.1046/j.1365-2419.2003.00262.x
[9] WANG Yanjun, LIU Qun, YE Zhenjiang. A bayesian analysis on the anchovy stock (Engraulis japonicus) in the Yellow Sea[J]. Fish Res, 2006, 82(1/3): 87-94. doi: 10.1016/j.fishres.2006.08.007
[10] 詹秉义. 渔业资源评估[M]. 北京: 中国农业出版社, 1995: 208-215. https://book.douban.com/subject/2048040/ [11] 邓景耀, 叶昌臣. 渔业资源学[M]. 重庆: 重庆出版社, 2000: 124-131. https://www.dushu.com/book/10041543/
计量
- 文章访问数: 5383
- HTML全文浏览量: 157
- PDF下载量: 3466
粤公网安备 44010502001741号
