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

4 期 王子琦，等：舟山渔场主要虾类资源状况及与环境因子的关系 50 卷

N
N
江苏省 W E
Jiangsu Province
S
31°30′
T1 T2 T18
! ! !
上海市
Shanghai City
T3 T4 T5 T6 T19 T20
! ! ! ! ! !
T7 T8 T9 T10 T21 T22
! ! ! ! ! !
舟山群岛
30°00′ Zhoushan Archipelago
T11 T12 T13 T14 T23 T24
! ! ! ! ! !
浙江省
Zhejiang Province T15 T16 T17
! ! !

! 站位 station
120°00′ 121°30′ 123°00′ 124°30′ 126°00′ E

图 1 舟山渔场桁杆拖虾作业渔船站位示意图

Fig. 1 Schematic diagram of beam trawl operation stations in the Zhoushan Fishing Ground
量渔获量 (t/h)。其次，通过 SE 反映均值估计的稳 2020 年 空 间 分 辨 率 采 用 0.125°×0.125°， 2021 —
定性与样本波动程度。 2024 年采用 0.083°×0.083°；SST 在 2020—2024 年
参照相关文献 [11-12] 的方法，通过 CPUE 数据， 空间分辨率采用 0.125°×0.125°。上述环境数据均
计算研究某时间段内虾拖网渔业各类渔获资源 按采样站位的经纬度位置进行空间匹配，并结合
的“重心”位置的月际、年际间时空变化规律。 月度及年度与渔业观测数据进行对应，用于后续
公式： 模型构建与环境响应分析。
k ) k )
∑ ( ∑ ( 为定量评价 GAM 的拟合表现，本研究以
CPUE ij ·x i CPUE ij ·y i
x = i ,y = i ln(CPUE+1) 作为响应变量进行建模。模型拟合优
∑ ∑
k k
2
CPUE ij CPUE ij 度与误差水平分别采用决定系数 R 和均方根误差
i i
2
式中，x 和 y 分别表示渔场重心的经度以及纬度； RMSE 进行衡量。R 反映模型对响应变量变异的
i 表示渔区编号，j 表示时间段，CPUE 则表示某 解释程度，公式：
j
i,
个渔区 i 在时间段 j 的单位捕捞努力量渔获量。将 2 SSE
R = 1−
每个渔区的 CPUE 值作为权重，计算所有渔区经 SST
纬度的加权平均值。CPUE 值越高，则该渔区对 式 中 ， SSE 为 残 差 平 方 和 ， SST 为 总 平 方 和 ；
重心位置影响越大。 RMSE 用于表征拟合误差的总体量级， y i 为观测
在软件 R 4.4.3 中利用 mgcv 数据包对相关数 值， ˆ y i为模型拟合值，公式：
据进行构建 GAM 模型，将主要利用虾类的资源 v
t
1 n ∑ 2
密度 (CPUE) 作为响应变量，将相关的环境以及时 RMSE = (y i − ˆy i )
n
空因子列为解释变量，分析它们之间的关系，并 i=1
探究各个因子同主要虾类 CPUE 的显著性关系 。 上述指标均在 ln(CPUE+1) 尺度下由观测值与
[10]
本研究选取了年份 (Year)、月份 (Month)、经度 拟合值计算获得。结果表明，两种虾类模型均获
(Lon)、纬度 (Lat) 为时空因子，DO、SST、SSS、 得较为一致的拟合效果。将所选因子放入模型中
Chl.a 为海洋环境因子进行分析。本研究所用环 进行拟合，通过赤池准则信息来对模型的拟合程
境因子数据的时空分辨率如下：Chl.a、DO 在 度进行检验，确定当 AIC 值最小，且偏差解释率
[13]
2020—2024 年空间分辨率为 0.25°×0.25°；SSS 在 为最大时，则为最优模型 。通过对相关变量进
中国水产学会主办 sponsored by China Society of Fisheries https://www.china-fishery.cn
3