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

4 期 水 产 学 报 50 卷

将鱼卵及仔稚鱼丰度进行数据标准化，以其 环境因子与全部调查月份的鱼卵或仔稚鱼优势种
丰度 (Abundance) 即单位滤水量采集到的鱼卵、仔 和重要种生态密度数据进行环境变量矩阵与群落
稚鱼数 [3] 构建 CLUSTER 及 NMDS 中 Bray-Curtis 距离矩阵的相关性分析时，设定显著性水平为
相似度矩阵的密度数据 (尾/m )，公式： P<0.05。相关系数 r>0 表示二者呈正相关，r<0 表
3
M 示呈负相关。
Abu = (6)
R×q×π×r 2
2 结果
式中，M 为鱼卵、仔稚鱼数量 (尾)；R 为网口流
量计转数差值 (r)；q 为网口流量计距离转数比
2.1 水体基础环境因子的时空动态
(m/r)；r 为大型浮游生物水平拖网的网口半径 (m)。
对鱼卵与仔稚鱼资源密度数据矩阵，分别采 2024 年 3 月—2025 年 2 月，四十里湾西部海
用四次方根及 lg(X+1) 转换进行转化处理，以降低 域 SST 呈显著季节变化 (图 2)，于 8 月达到年内
优势种和稀有种对分析结果的影响。在此基础上， 峰值，随后逐步下降，至次年 2 月降至全年最低。
以月份为尺度单位进行聚类分析 (cluster analysis) SBT 变化趋势与 SST 一致但变幅较小，空间上总
和非度量多维标度分析 (NMDS)，以进一步揭示 体表现为近岸高、外海低；全年各月 SST 与 SBT
鱼卵与仔稚鱼群落在时间尺度上的结构变化规律， 差异不显著 (U=18 712.5, P=0.055, Mann-Whitney U-
NMDS 结果的可靠性通过胁强系数 (stress) 进行检 test)(图 3-a)。
验，当 stress 值小于 0.05 时，表明分析结果具有 3—5 月 SSS 空间差异不明显 (图 4)，6 月呈
较 好 的 代 表 性 [18-19] 。 同 时 ， 采 用 相 似 性 分 析 东低西高趋势，7 月表现为近岸低、外海高，9 月
(ANOSIM)，以 Global R 值和 P 值作为评估分组结 及次年 2 月受鲁北低盐沿岸水入侵影响，盐度梯
构显著性的指标。进一步通过 SIMPER 方法对聚 度显著增大并呈现自西北湾口向东南湾底部递增
类分析所得组别进行解析，以探明不同群落间的 格局。SBS 与 SSS 变化基本一致，仅 8 月表现为
差异，并依据平均相异性贡献率 (累积至 90%) 自外海向近岸递减；全年各月 SSS 与 SBS 无显著
与平均相似性贡献率 (累积至 90%) 确定组内主要 差异 (U=19 015.0, P=0.132>0, Mann-Whitney U-test)。
特征物种及组间的主要分歧物种。 四十里湾西部海域主要渔业生物产卵期为
[3]
1.4 鱼类早期资源优势种和重要种类时空分布 5—7 月 ，表层五项溶解无机营养盐浓度呈明显
月际差异，整体表现为 6 月浓度普遍偏高。空间
与环境因子关系分析
上 ， 高 值 区 主 要 集 中 于 近 岸 及 逛 荡 河 入 海 口
鉴于鱼卵和仔稚鱼主要分布于海水表层，表 +
(S23 站) 等陆源输入明显的区域。其中，NH 在
4
层环境因子对其影响更显著，故最终选取 SST、
6 月整体浓度最高，且呈现由近岸向离岸递减的
SSS、Depth、DO、pH 和 WW 作为解释变量用于 3− − − 2−
梯度分布特征；PO 、NO 、NO 及 SiO 的浓
3
4
3
2
后续 Mantel 检验和 CCA 分析。
度峰值亦出现在 S23 站附近 (图 5)。
将筛选后环境因子与各月鱼卵或仔稚鱼优
2.2 群落结构与物种多样性
势种及重要种生态密度数据进行去趋势对应分
析 (DCA) 以揭示环境变量与物种分布的关系。 2024 年 3 月—2025 年 2 月，四十里湾西部海
根据 DCA 结果中物种变量矩阵第一排序轴梯度 域共采集鱼卵 346 608 粒，仔稚鱼 1 071 尾。经鉴
长度确定适用的排序分析方法：若梯度长度大 定共记录鱼卵、仔稚鱼 45 种，隶属于 11 目 27 科。
于 4.0，采用典范对应分析 (CCA)；若小于 3.0， 其中鱼卵 21 种 (7 目 15 科)，以鲈形目为主；仔稚
则采用冗余分析 (RDA)；若介于 3.0 至 4.0，则两 鱼 33 种 (8 目 18 科)，同样以鲈形目占优势。IRI
[20]
种方法均适用 。本研究鱼卵和仔稚鱼第一排序 结果显示，鱼卵群落 3 月以石鲽 (Kareius bicol-
轴的梯度分别为 5.43 和 7.29，均选用 CCA 开展 oratus) 为 重 要 种 ； 4 月 牙 鲆 (Paralichthys oli-
后续分析。 vaceus) 为优势种，鳀和弯棘䲗(Callionymus curvic-
Mantel 检验适用于检验两个矩阵之间的相关 ornis) 为重要种；5—6 月鳀为优势种，分别伴随
[21]
性，该非参数统计方法对数据分布无特定要求 ， 黑鲷 (Acanthopagrus schlegelii)、弯棘䲗及少鳞鱚
已广泛应用于群落结构研究 [22-24] 。在对筛选后的 (Sillago japonica) 和叫姑鱼 (Johnius grypotus) 等重

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