摘要
为兼顾调查成本和调查效果并实现多目标远洋渔业资源调查,根据2021—2023年西北太平洋渔业资源调查数据,以澳洲鲐、日本鲭和远东拟沙丁鱼的相对资源量以及总相对资源量为调查目标,使用Kriging插值法模拟调查目标的空间分布,利用计算机模拟方法比较了简单随机抽样(Design 1)、系统抽样(Design 2)、基于层面积分配方法的分层随机抽样(Design 3)和基于层资源量分配方法的分层随机抽样(Design 4)等4种采样设计对调查目标的评估效果,以相对误差(Relative estimation error, REE)和相对偏差(Relative bias, RB)衡量评估效果的精度,采用准确度变化率(Accuracy change rate, ACR)对调查站点数进行选择。结果显示,4种采样设计对调查目标的评估精度不同,表现为Design1<Design2<Design3<Design4。Design1、Design3和Design4的REE值随站点数量增加而降低,Design2的REE值随站点数量增加波动较大。Design4对日本鲭、远东拟沙丁鱼和总相对资源量的评估精度更高,但高估了日本鲭的相对资源量。当采用Design4进行模拟采样时,综合考虑所有调查目标,28个站位为最适调查站位数量。结果表明,分层随机抽样(Design4)可优先作为西北太平洋渔业资源调查采样设计方案。
关键词
渔业资源调查是研究鱼类生长特
西北太平洋区域是世界粮农组织(Food and Agriculture Organization of the United Nations,FAO)统计的所有区域中产量最高的区
数据来源于2021—2023年“淞航”号调查船在西北太平洋开展的渔业资源综合调查,调查时间为2021—2023年6—8月。调查海域为148°E~164°E、35°N~45°N,共设70个综合调查站点,调查站点栅格经纬度大小为2°×1°,其中2021年渔业资源拖网站点42个,2022年36个,2023年39个(

图1 2021—2023年西北太平洋调查站点空间分布
Fig.1 Sampling stations in Northwest Pacific Ocean from 2021 to 2023
本研究通过Kriging插值法获得调查目标的资源分布情况作为模拟抽样的总体,以70个综合调查站点为总体潜在调查站点(
由于在2021年和2022年的调查中,未将获取的鲐进行日本鲭和澳洲鲐的区分,同时2021—2023年的总相对资源量和远东拟沙丁鱼的相对资源量空间分布特征较为相似(

图2 2021—2023年总相对资源量和远东拟沙丁鱼相对资源量空间分布
Fig.2 Spatial distribution of total relative abundance and the relative abundance of Sardinops malanostictus in 2021-2023
本研究比较了SRS、SYS和StS在估计各目标鱼种相对资源量及总相对资源量的表现。
SRS:每次抽样从70个潜在调查站点中随机抽取n个站点,进行不放回抽样,相应采样设计为Design 1。
SYS:按调查船航行路线对潜在调查站点进行编号,随机抽取站点m作为起点,再按照固定间隔k=70/n从总体中抽取n个站点(若k不为整数,则对k向下取整),进行环形系统抽样,相应采样设计为Design 2。
StS:根据各目标鱼种相对资源量及总相对资源量的空间分布特征,按照层间异质性大、层内异质性小的分层原则将调查水域划分为3层(相对资源量由高到低分别为A、B、C,

图3 目标鱼种相对资源量及总相对资源量预测分布及分层示意图
Fig.3 Stratified schematic diagram illustrating the predicted distribution of the total relative abundance and the relative abundance of target species
采样站点按照每层包含的潜在站点数量比例进行分配,即层内的调查站点数与层面积呈正比,相应采样设计为Design 3。
采样站点按照每层的相对资源量比例进行分配,即层内的调查站点数与层相对资源量呈正比,若分配的站点数超过该层总站点数量,则将超出部分的站点数按照相对资源量比例分配给其他层,相应采样调查设计为Design 4。
采用相对误差RE
REE与RB的计算公式:
REE=100% | (1) |
RB=100% | (2) |
式中:V为采样设计第i次模拟抽样调查的估计值均值;V为目标鱼种相对资源量和总相对资源量“真实”值的均值;N为模拟次数,本研究中N=1 000。
基于1 000次重抽样的REE均值来计算准确度变化率(Accuracy change rate,ACR),以此作为衡量REE变化量与站点数量变化的比
ACR= | (3) |
式中:Sn、Smin为不同站点数,Sn大于Smin,且Smin为站点数在[3,n-1]区间内REE最小值对应的站点数;REEn、REEmin为不同站点数对应的1 000次REE均值,选取ACR值为0.2作为确定REE显著降低的标准线。
总体来看,简单随机抽样和分层随机抽样对4种调查目标总体估计的REE变化趋势一致,均随站点数量的增加而降低,且2种分层随机抽样的REE小于简单随机抽样(

图4 不同采样设计对各调查目标抽样相对误差(REE)变化趋势
Fig.4 REE values of the different sampling designs based on different objectives
对于不同调查目标,4种采样设计的表现不同。当以澳洲鲐相对资源量为调查目标时,分层随机抽样的REE值较小,当抽样站点数<35时,Design 3的REE值最小,当抽样站点数>35时,Design 3和Design 4的REE值趋于一致(
当以澳洲鲐、远东拟沙丁鱼的相对资源量和总相对资源量为调查目标时,4种采样设计的RB值在-5%~5%的数值范围内围绕0值上下波动,均为无偏估计(

图5 不同采样设计对各调查目标抽样相对偏差(RB)变化趋势
Fig.5 RB values of the different sampling designs based on different objectives
Design 1的ACR均大于0,并随抽样站点数量的增加而降低(图

图6 不同采样设计的ACR 值随样本量的变化
Fig.6 Changes in accuracy change rate (ACR) of different sampling designs with different sampling stations
当对不同调查目标模拟抽样时,Design 4的REE值较小且最为稳定,最优采样设计可选择为Design 4。当Design 4对不同的调查目标进行模拟抽样时,最适调查站点数量存在不同(
调查目标 Survey objectives | REE范围 Range of REE | 最适调查站点数 Optimal sampling stations | 最适调查站点数时的REE REE of optimal sampling stations |
---|---|---|---|
澳洲鲐 Scomber australasicus | 1.20~6.51 | 28 | 3.01 |
日本鲭 Scomber japonicus | 12.27~31.97 | 26 | 13.33 |
远东拟沙丁鱼Sardinops melanostictus | 0.06~13.16 | 18 | 7.85 |
总相对资源量 Total relative abundance | 0.71~10.44 | 20 | 5.68 |
本研究以多个鱼种相对资源量以及总相对资源量为调查目标,比较了4种采样设计在不同调查目标下的抽样表现。研究发现,相比简单随机抽样和系统抽样,分层随机抽样的模拟抽样表现较好。简单随机抽样是最简单的抽样方法,其他抽样方法在某种程度上包含简单随机抽样。本研究中简单随机抽样的REE、RB随抽样站点数的增加而降低,但抽样精度相对不高,这与相关研
系统抽样可以覆盖调查区域范围,并能兼顾随机性,常被作为初始调查方案,但在某些调查努力量范围内,其调查精度存在波
分层随机抽样的估计效果比简单随机抽样和系统抽样精度高,这与相关研
调查目标 Survey objectives | A层 Stratum-A | B层 Stratum-B | C层 Stratum-C | 总体 Total |
---|---|---|---|---|
澳洲鲐 Scomber australasicus | 0.047 0 | 0.074 3 | 0.012 8 | 0.558 7 |
日本鲭 Scomber japonicus | 897.091 2 | 22.064 4 | 0.391 0 | 813.812 4 |
远东拟沙丁鱼 Sardinops melanostictus | 1 830.057 0 | 993.998 2 | 0.112 4 | 6 941.720 5 |
总相对资源量 Total relative abundance | 2 248.528 0 | 1 083.006 0 | 27.739 1 | 8 572.700 8 |
以日本鲭相对资源量为调查目标时,Design4对日本鲭存在高估现象,这可能与日本鲭的空间分布较为集中有
在各目标鱼种及总相对资源量真值获取上,本研究使用了Kriging插值法,并未考虑环境及其他因素对物种分布的影响,存在一定的局限性。联合物种分布模型和广义加性模型可以分析不同物种之间的相互作用及对其分布产生的影响,并提供更准确的预测结
参考文献
柳昊, 张文超, 李文甲, 等. 黄海中部夏季鳀鱼卵、仔稚鱼分布与早期生长特征[J]. 上海海洋大学学报, 2024, 33(1): 172-185. [百度学术]
LIU H, ZHANG W C, LI W J, et al. Distribution and early growth characteristics of anchovy eggs, larvae and juveniles in the central Yellow Sea in summer[J]. Journal of Shanghai Ocean University, 2024, 33(1): 172-185. [百度学术]
王淼娣, 王雪辉, 杜飞雁, 等. 基于LBB方法估算北部湾竹䇲鱼种群参数[J]. 上海海洋大学学报, 2022, 31(1): 212-222. [百度学术]
WANG M D, WANG X H, DU F Y, et al. Estimation of the population parameters of Trachurus japonicus in the Beibu Gulf based on the LBB method[J]. Journal of Shanghai Ocean University, 2022, 31(1): 212-222. [百度学术]
鲁红月, 汪金涛, 张畅. 基于产卵场海表面温度的太平洋褶柔鱼冬生群资源评估[J]. 上海海洋大学学报, 2022, 31(3): 650-657. [百度学术]
LU H Y, WANG J T, ZHANG C. Stock assessment of winter-spawning cohort for Todarodes pacificus based on the sea surface temperature in the spawning ground[J]. Journal of Shanghai Ocean University, 2022, 31(3): 650-657. [百度学术]
孙菁煜, 戴小杰, 朱江峰, 等. 淀山湖鱼类多样性分析[J]. 上海水产大学学报, 2007, 16(5): 454-459. [百度学术]
SUN J Y, DAI X J, ZHU J F, et al. Analysis of the fish species diversity in Dianshan Lake[J]. Journal of Shanghai Fisheries University, 2007, 16(5): 454-459. [百度学术]
ZHAO J, YANG K E, MA J. Optimization of sampling effort for different fishery groups in the Yangtze river estuary, China[J]. Marine and Coastal Fisheries, 2022, 14(4): e10214. [百度学术]
STARR R M, CARR M, MALONE D, et al. Complementary sampling methods to inform ecosystem-based management of nearshore fisheries[J]. Marine and Coastal Fisheries, 2010, 2(1): 159-179. [百度学术]
HOWARD R A, CIANNELLI L, WAKEFIELD W W, et al. Comparing fishery-independent and fishery-dependent data for analysis of the distributions of Oregon shelf groundfishes[J]. Fisheries Research, 2023, 258: 106553. [百度学术]
SCHEIRER K, CHEN Y, WILSON C. A comparative study of American lobster fishery sea and port sampling programs in Maine: 1998-2000[J]. Fisheries Research, 2004, 68(1/3): 343-350. [百度学术]
ZHAO J, CAO J, TIAN S Q, et al. Evaluating sampling designs for demersal fish communities[J]. Sustainability, 2018, 10(8): 2585. [百度学术]
赵静, 章守宇, 林军, 等. 不同采样设计评估鱼类群落效果比较[J]. 应用生态学报, 2014, 25(4): 1181-1187. [百度学术]
ZHAO J, ZHANG S Y, LIN J, et al. A comparative study of different sampling designs in fish community estimation[J]. Chinese Journal of Applied Ecology, 2014, 25(4): 1181-1187. [百度学术]
SHI Y C, ZHANG X M, HE Y R, et al. Stock assessment using length-based Bayesian evaluation method for three small pelagic species in the northwest Pacific ocean[J]. Frontiers in Marine Science, 2022, 9: 775180. [百度学术]
刘思源, 张衡, 杨超, 等. ENSO事件下西北太平洋远东拟沙丁鱼和日本鲭栖息地协同变化特征[J]. 海洋学报, 2024, 46(1): 39-52. [百度学术]
LIU S Y, ZHANG H, YANG C, et al. Exhibit covariation characteristics in the habitat changes of Sardinops melanostictus and Scomber japonicus in the northwestern Pacific Ocean under ENSO event[J]. Haiyang Xuebao, 2024, 46(1): 39-52. [百度学术]
蔡锴, 周雨霏, 麻秋云, 等. 西北太平洋鲐的资源养护和管理研究[J]. 中国水产科学, 2024, 31(4): 439-453. [百度学术]
CAI K, ZHOU Y F, MA Q Y, et al. Research on the conservation and management for chub mackerel (Scomber japonicus) in the Northwest Pacific[J]. Journal of Fishery Sciences of China, 2024, 31(4): 439-453. [百度学术]
马金. 基于计算机模拟的河口渔业资源抽样调查设计比较与优化——以长江口为例[D]. 上海: 上海海洋大学, 2020. [百度学术]
MA J. Simulation-based evaluation of sampling design for different fishery groups in estuary: a case study in the Yangtze estuary[D]. Shanghai: Shanghai Ocean University, 2020. [百度学术]
LIU Y, CHEN Y, CHENG J H. A comparative study of optimization methods and conventional methods for sampling design in fishery-independent surveys[J]. ICES Journal of Marine Science, 2009, 66(9): 1873-1882. [百度学术]
CHEN Y. A Monte Carlo study on impacts of the size of subsample catch on estimation of fish stock parameters[J]. Fisheries Research, 1996, 26(3/4): 207-223. [百度学术]
CAO J, CHEN Y, CHANG J H, et al. An evaluation of an inshore bottom trawl survey design for American lobster (Homarus americanus) using computer simulations[J]. Journal of Northwest Atlantic Fishery Science, 2014, 46: 27-39. [百度学术]
MA J, TIAN S Q, GAO C X, et al. Evaluation of sampling designs for different fishery groups in the Yangtze River estuary, China[J]. Regional Studies in Marine Science, 2020, 38: 101373. [百度学术]
MA J, LI B, ZHAO J, et al. Environmental influences on the spatio-temporal distribution of Coilia nasus in the Yangtze River estuary[J]. Journal of Applied Ichthyology, 2020, 36(3): 315-325. [百度学术]
张国晟, 王晶, 张崇良, 等. 基于多目标鱼种资源量指数估计的不同采样设计比较[J]. 水产学报, 2021, 45(5): 700-715. [百度学术]
ZHANG G S, WANG J, ZHANG C L, et al. Comparison of sampling designs of fishery-independent survey in estimating abundance indices of multiple target species[J]. Journal of Fisheries of China, 2021, 45(5): 700-715. [百度学术]
金勇进, 杜子芳, 蒋妍. 抽样技术[M]. 2版. 北京: 中国人民大学出版社, 2008. [百度学术]
JIN Y J, DU Z F, JIANG Y. Sampling technique[M]. 2nd ed. Beijing: China Renmin University Press, 2008. [百度学术]
刘勇. 渔业资源评估抽样调查方法的理论探讨与研究[D]. 上海: 华东师范大学, 2012. [百度学术]
LIU Y. Theoretical study on the sampling methods of survey for fishery stock estimation[D]. Shanghai: East China Normal University, 2012. [百度学术]
刘思源, 张衡, 杨超, 等. 基于最大熵模型的西北太平洋远东拟沙丁鱼和日本鲭栖息地差异[J]. 上海海洋大学学报, 2023, 32(4): 806-817. [百度学术]
LIU S Y, ZHANG H, YANG C, et al. Differences in habitat distribution of Sardinops melanostictus and Scomber japonicus in the northwest Pacific based on a maximum entropy model[J]. Journal of Shanghai Ocean University, 2023, 32(4): 806-817. [百度学术]
周茉, 方星楠, 余为, 等. 厄尔尼诺和拉尼娜事件下西北太平洋柔鱼栖息地时空分布差异[J]. 上海海洋大学学报, 2022, 31(4): 984-993. [百度学术]
ZHOU M, FANG X N, YU W, et al. Difference of spatio-temporal distribution of neon flying squid Ommastrephes bartramiii in the northwest Pacific Ocean under the El Niño and La Niña events[J]. Journal of Shanghai Ocean University, 2022, 31(4): 984-993. [百度学术]
潘邵媛, 王学昉, 田思泉, 等. 不同采样设计对长江口海洋环境监测效果影响的评价[J]. 海洋环境科学, 2021, 40(3): 435-441. [百度学术]
PAN S Y, WANG X F, TIAN S Q, et al. Evaluation of the effect of different sampling designs on marine environmental monitoring in the Yangtze River Estuary[J]. Marine Environmental Science, 2021, 40(3): 435-441. [百度学术]
潘邵媛, 王学昉, 田思泉, 等. 长江口中华鲟保护区海洋环境监测浮标站点的优化设计[J]. 海洋学报, 2021, 43(4): 55-64. [百度学术]
PAN S Y, WANG X F, TIAN S Q, et al. The design of the stations of marine environmental monitoring buoys in the Chinese sturgeon nature reserve in the Changjiang River Estuary[J]. Haiyang Xuebao, 2021, 43(4): 55-64. [百度学术]
刘思源, 张衡, 杨超, 等. 西北太平洋远东拟沙丁鱼与日本鲭种群动态特征及其与环境因子关系研究进展[J]. 大连海洋大学学报, 2023, 38(2): 357-368. [百度学术]
LIU S Y, ZHANG H, YANG C, et al. Relationship between stock dynamics and environmental variability for Japanese sardine (Sardinops sagax) and chub mackerel (Scomber japonicus) in the Northwest Pacific Ocean: a review[J]. Journal of Dalian Ocean University, 2023, 38(2): 357-368. [百度学术]
刘彦琳. 基于GAM和MaxEnt的西北太平洋公海日本鲭渔场形成机制及潜在栖息地分布研究[D]. 上海: 上海海洋大学, 2023. [百度学术]
LIU Y L. Based on GAM and MaxEnt the formation mechanism and potential habitat distribution of Scomber japonicus fisheries in the high seas of the Northwest Pacific Ocean[D]. Shanghai: Shanghai Ocean University, 2023. [百度学术]
WANG J, XU B D, ZHANG C L, et al. Evaluation of alternative stratifications for a stratified random fishery-independent survey[J]. Fisheries Research, 2018, 207: 150-159. [百度学术]
ZHANG C L, XU B D, XUE Y, et al. Evaluating multispecies survey designs using a joint species distribution model[J]. Aquaculture and Fisheries, 2020, 5(3): 156-162. [百度学术]
张涛, 赵天亚, 栾静, 等. 基于物种分布模型分析环境因子对海州湾偶见种资源分布的影响[J]. 海洋学报, 2023, 45(7): 69-78. [百度学术]
ZHANG T, ZHAO T Y, LUAN J, et al. Analysis of the influence of environmental factors on the distribution of occasional species in the Haizhou Bay based on species distribution model[J]. Haiyang Xuebao, 2023, 45(7): 69-78. [百度学术]
XU B D, REN Y P, CHEN Y, et al. Optimization of stratification scheme for a fishery-independent survey with multiple objectives[J]. Acta Oceanologica Sinica, 2015, 34(12): 154-169. [百度学术]
王雅萌, 陈芃, 陈新军. 气候变化下西北太平洋大海洋生态系海表面温度特征分析[J]. 上海海洋大学学报, 2021, 30(5): 874-883. [百度学术]
WANG Y M, CHEN P, CHEN X J. Characteristics of sea surface temperature for large marine major ecosystems in the northwest Pacific under climate change[J]. Journal of Shanghai Ocean University, 2021, 30(5): 874-883. [百度学术]