Page 15 - 《水产学报》2026年第04期
P. 15

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

              选取多个指标评估模型与数据的适配性，判定标                            心区域多呈现负载荷，暗示其对整体变化的响应
              准：χ 检验     P > 0.05、均方根误差近似值< 0.08、              较弱或存在非同步性。2001—2025            年小黄鱼密度
                    2
              比较拟合指数 > 0.9、标准化残差均方根< 0.08；                     在较多区域 (尤其是密度热点核心区) 的                PC1 载荷
              计算克朗巴赫系数 (Cronbach’s α)，α > 0.7       表示观        为正，表明越冬场核心区域的密度变化趋势具有
              测变量对潜变量的代表性良好；P < 0.05。结果显                       一致性。线性趋势分析显示密度年际变化存在空
              示，本研究构建的         SEM  模型通过了上述检验。                 间异质性，研究区域西部主要呈正斜率，表明近
                   上述模型运行与结果绘制均使用              R  语言实现。        年来该海域小黄鱼资源丰度增加；西南部呈较大

                                                               负斜率 (< –0.2) 表明，该海域小黄鱼种群恢复态
               2    结果
                                                               势不佳或越冬场的迁移。这表明小黄鱼生物量密
                                                               度长期变化在越冬场具有空间异质性。
               2.1    小黄鱼在越冬场的时空分布模式
                                                                    空间聚类分析  通过 K–means 与层次聚类
                    分布模式与指标变化  本研究构建的时空                        确定的小黄鱼越冬场的空间分区。两种聚类结果
              模型满足收敛标准。研究期间首尾和中间年份的                            较为一致，小黄鱼越冬群体具有清晰的空间结构；
              小黄鱼密度和热点分布图显示，2001                年，黄海中          最优聚类数为       2，即越冬场可划分为           2  块区域：
              部、东南部存在较为广泛的小黄鱼高密度分布区，                           ①北黄海、中黄海西部和南黄海大部分海域；②
              黄海北部也是小黄鱼重要的栖息地。小黄鱼种群                            黄海东部和黄海东南部海域，后者也是高生物量
              密度自    2001  年后总体呈下降趋势，密度热点数量                    密度的主要分布区 (图         5)。
              减少、收缩至黄海中部，近年来开始增加，2025
                                                                2.2    小黄鱼资源状况的长期演变
              年黄海南部密度增加程度高于黄海北部海域 (图                     2-
              a, b)。空间各向异性分析表明，小黄鱼相遇概率                             运行结果显示，LIME          模型满足收敛标准。
              的空间相关尺度显著大于正渔获率，且两者均表                            小黄鱼捕捞死亡系数         F  在  2001—2005  年处于较低
              现出明显的方向性特征，即             2001—2025  年小黄鱼         水平，随后呈波动上升趋势，在               2015—2021  年严
              在越冬场沿西北-东南方向分布 (图              2-c)。             重超过目标参考点         F targe t  并超过限制参考点    F limit
                   小黄鱼种群分布东向重心 (图             3-a) 在初期位        处于过度捕捞水平 (图         6)。近年来，F     呈现急剧下
              于研究区域偏东位置，之后总体呈现显著西移 (P                          降，在    2023—2025  年低于目标参考点         F target ，处
              < 0.05)，至  2025  年主要分布在研究区域偏西位置，                 于可持续捕捞压力水平。SPR             与  F  呈现明显相反
              表明整体分布范围在经向上发生了显著的向西移                            趋势，在    2015—2021  年处于不可持续的资源衰退
              动。相比之下，北向重心在             2001—2010  年往北部         状 态 ， 之 后 在   2023  年 开 始 恢 复 到 目 标 参 考 点
                                                                   0
              迁移，但在整个研究期间无纬度方向的显著性                             SPR 以上。
                                                                  4
              移 动   (P  > 0.05)。 小 黄 鱼 种 群 有 效 分 布 面 积 在           小黄鱼补充量表现出高度的年际波动性，在
              2001—2025  年呈现先上升后下降的趋势，与生物                      2001、2007—2008   和  2024  年都观察到明显的“脉
              量密度热点趋势相反。                                       冲式”补充高峰。近年来，小黄鱼补充量总体维持
                    密度空间分布和种群动态时序差异  2001—                     在中低水平，少数强补充年份与较长时间的低补
              2025  年，黄海中部和东南部越冬场小黄鱼的中位                        充期交替出现，说明小黄鱼补充量具有较强的不
              数密度分布较高 (图        4-a)。中部大部分区域显示出                稳定性，可能受环境和亲体生物量的共同影响。
              较低的    CV (图  4-b)，表明这些区域个体数量持续                  相对产卵群体生物量指数轨迹显示了典型的资源
              较高，在时间上较为稳定。中位数密度低的区域                            衰退过程。该指数从           2001  年的高值迅速下降，
              总体也呈现相对低的          CV。                           2010  年后长期维持在低水平，近年来呈现有所恢
                   时间序列的第一主成分 (PC1) 代表小黄鱼密                     复的迹象。模型预测的平均体长 (绿线) 与观测数据
              度变化的主导时间模式，解释了密度总变异的                             (黑点) 拟合良好。平均体长在初期较高，随后随
              82.8%。PC1 的正载荷主要与高密度核心区重叠                        着捕捞强度的增加而下降，反映了种群中大个体
              (图  4-c)，说明这些区域具有高度同步性的年际变                       被移除的现象 (截断效应)。渔具选择性显示，当
              化，并在时间尺度变异中起主导作用；分布非核                            前  S 和 0  S 分别为  9.79  和  12.96 cm。
                                                                        5
                                                                        9
                                                                  5
              中国水产学会主办  sponsored by China Society of Fisheries                          https://www.china-fishery.cn
                                                            5
   10   11   12   13   14   15   16   17   18   19   20