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

4 期 郭弘艺，等：长江口及邻近海域刀鲚的季节性时空分布特征 50 卷

N 0.05°×0.05°，按季节分别实施。
33°
1.4 空间分布特征分析
采样点 sampling sites
参照 Woillez 等 [24] 的方法，基于站位尺度资
江苏省
Jiangsu Province 源密度数据，计算 CG、惯性 (inertia)、惯性椭圆
32°
及各向同性指数 (isotropy)，并估算当年生群体与
非当年生群体间的全局共位指数 (global index of
崇明岛
长江口 collocation，GIC) 与空间分离指数 (separation)，以
Chongming Island
Yangtze Estuary
定量表征两类群体在不同季节的空间分布格局及
31° 上海市
Shanghai City 其相互关系。
东海
杭州湾 East China Sea 数据分组与筛选 空间指标按年份 (2024
Hangzhou Bay 年/2025 年) ×季节 (春季/夏季/秋季) ×年龄组 (当年
0 25 50 km 生群体/非当年生群体) 分组计算，以各站位经纬
30° `
121° 122° 123° 124° E 度及对应资源密度为基础数据。为避免零密度站
图 1 长江口及邻近海域采样站位分布图 位对位置和扩散估计造成偏差 ，仅纳入资源密
[24]
Fig. 1 Distribution of sampling sites in the Yangtze 度大于 0 的站位参与计算。若某组有效站位数不
River Estuary and adjacent coastal waters 足 3 个，则无法可靠构建惯性椭圆，该组指标记
为缺失值 (NA)，不纳入后续分析。
出现率 (frequency of occurrence, FO，%) 基于
几何分布指标计算 CG 计算公式：
站位尺度的有无数据计算：
∑ ∑
x i z i y i z i
N occ
FO = ×100%
N total CG x = ，CG y =
∑ ∑
式中， N occ 为刀鲚出现的站位数， N total 为调查总 z i z i
站位数。该指标用于表征刀鲚在不同季节的分布
式中， x i 和 分别为第 i 个站位的经度与纬度； z i
y i
普遍程度。
为该站位的资源密度。由于惯性相关指标以距离
对捕获个体进行生物学测定。采用精度为
为度量基础，在每个分组内以该组分布重心纬度
0.1 cm 的测量板测定体长 (SL)。完成外部形态测
CG 为局地近似基准，将经纬度偏差转换为公里
y
量后，解剖个体，通过目测性腺形态与色泽判定
尺度：
性别及性腺发育阶段，并依据 6 期成熟度分级系
统 (Ⅰ~Ⅵ期) 进行划分 [21-22] 。 ∆x km = (x i −CG x )×111×cosCG y
( )
年龄鉴定参照袁传宓等 [23] 的方法，以鳞片密 ∆y km = y i −CG y ×111
带数量为依据。鳞片密带形成于个体生长缓慢时 据此构建密度加权协方差矩阵 (惯性矩阵) S：
期，通常对应越冬阶段，年龄以冬龄表示，即个  ∑ 2 ∑ 
 z i ∆x z i ∆x km ∆y km  
  km 
体经历的冬季数。YOY 定义为 0 龄组，代表当年   ∑ ∑  
  z i  
 
  z i  
出生且尚未经历越冬的个体；非当年生个体 (non-     
S = ∑
 ∑ 
  z i ∆y 2  
YOY) 为≥1 龄组，即至少经历过 1 个冬季的个体。     z i ∆x km ∆y km km    
 ∑ ∑ 
 
 
z i z i
1.3 水文环境观测与空间插值
对矩阵 S 进行特征值分解，得最大与最小特
各采样站位在开展底拖网生物调查的同时，
征值 λ max 和 λ min ，据此计算惯性：
同步采用 CTD 剖面仪 (SBE25plus，Sea-Bird Elec-
tronics，美国) 现场测量底层水温 (T，℃) 和盐度 inertia = λ max +λ min
(S，PSU)。为表征研究区底层环境的空间分布特 惯性表征群体围绕分布重心的整体空间扩散
征，采用反距离加权法 (inverse distance weighting， 尺度。特征值同时用于定义惯性椭圆，其长、短
IDW) 对站位观测值进行空间插值，插值分辨率为 半轴分别为：
中国水产学会主办 sponsored by China Society of Fisheries https://www.china-fishery.cn
3