Page 17 - 《水产学报》2026年第04期

4 期 韩青鹏，等：黄渤海小黄鱼种群动态及其驱动因素 50 卷

征，其路径系数仅为 –0.03，表明在当前模型架构
600
下，海水污染对该种群动态呈微弱负影响。综上，
580 驱动因子的相对重要性依次为捕捞压力、管理措
东向/km eastward 560 施、气候变化和水质环境 (图 7)。
540
520 种群动态潜变量与产卵潜力比 (载荷= 0.94)
和相对产卵群体生物量 (载荷=0.77) 均高度相关，
500
说明这两个指标能很好地反映种群的健康状态。
2000 2005 2010 2015 2020 2025
相比补充量的载荷较低 (0.25)，表明补充量的随机
年份/年
year 波动性较大，不能单独作为衡量种群状态的稳定
(a) 性指标。其他潜变量与观测变量之间具有较高的
载荷，例如，管理措施潜变量 (载荷>0.9) 和气候
变化潜变量 (载荷>0.8) 对应的观测指标具有极高
3 880
的解释力，证明了模型构建的合理性。
3 860
北向/km northward 3 840 3 讨论

3 820
3 800 3.1 小黄鱼种群的时空动态与空间管理

2000 2005 2010 2015 2020 2025 在已有研究 [3,15-16] 基础上，分析了黄渤海小黄
年份/年 鱼冬季时空分布和资源丰度的长期变化，综合分
year
析了密度热点、分布重心、空间各向异性、时间
(b)
变异性 (线性趋势) 及聚类结果，揭示了小黄鱼越
冬群体清晰且一致的空间分布特征；进一步明确
有效分布面积/ln(km²) effective area occupied 11.0 的海域亦是其重要越冬场。热点分布分析发现，
11.5
了越冬场的位置及其变化，黄海北部水深大于
50 m
其主要在水深
以深的海域越冬，少量分布
50 m
海域，极少分布在
以浅海域。空
在
30 m
30~50 m
10.5
间各向异性结果进一步表明，种群分布并非各向
同性，而是沿特定方向 (西北-东南) 受环境梯度约
2000 2005 2010 2015 2020 2025 束；这种方向性结构为密度热点的拉伸形态及其
年份/年 与主要水文特征 (水深和黄海暖流) 的空间一致性
year
提供了机制层面的解释，即小黄鱼在越冬场的分
(c)
布与物理海洋过程强耦合，受控于以等深线走向、
图 3 2001—2025 年越冬场小黄鱼的东向和北向重心
[37]
以及有效分布面积 潮流及水团锋面等为主导空间定向过程 。韩青
(a) 东向重心，(b) 北向重心，(c) 有效分布面积。 鹏等 [15] 研究同样发现，黄渤海小黄鱼越冬群体核
心密度区年间迁移和黄海暖流与大陆沿岸冷水系
Fig. 3 Eastward and northward center of gravity and
effective area occupied of the L. polyactis in the overwin- 强弱变动有关。
tering ground in each year of the period 2001-2025 单秀娟等 [17] 发现，小黄鱼分布重心受气候
(a) eastward COG, (b) northward COG, (c) effective area occupied. (PDO 的延迟效应)、捕捞和种群数量变动的影响，
持续向中国近岸迁移。捕捞压力和气候变化是决
健康起到了积极作用。
定黄渤海鱼类分布和种群数量变动的重要影响因
气候变化同样表现为正面影响，路径系数为 子 [10,38-39] ，这一结论在春、夏季渤海小黄鱼时空动
0.27。该潜变量主要由 PDO 及其滞后项 PDO_lag1 态研究中同样得到验证 。PDO 与东亚大气环流
[40]
定义，说明当前气候条件在一定程度上有利于该 及中国气候年代际变化关系密切。冬季 PDO 暖位
种群的生存或补充。水质环境以海水污染变量表 相期，阿留申低压与蒙古高压增强，长江流域至

中国水产学会主办 sponsored by China Society of Fisheries https://www.china-fishery.cn
7