网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

西北太平洋渔业资源调查采样设计比较  PDF

  • 张烨 1
  • 赵静 1
  • 叶旭昌 1
  • 刘志良 2
  • 马金 1,3,4
1. 上海海洋大学 海洋生物资源与管理学院,上海 201306; 2. 上海海洋大学 船舶管理中心,上海 201306; 3. 国家远洋渔业工程技术研究中心,上海 201306; 4. 大洋渔业资源可持续开发教育部重点实验室,上海 201306

中图分类号: S 932.4

最近更新:2024-07-12

DOI: 10.12024/jsou.20240404526

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

为兼顾调查成本和调查效果并实现多目标远洋渔业资源调查,根据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)可优先作为西北太平洋渔业资源调查采样设计方案。

渔业资源调查是研究鱼类生长特

1、种群分2、生物多样性以及渔业资源变3-4,掌握渔业资源现状的重要途5,对渔业资源管理措施或政策的制定与调整有重要意6。渔业资源调查可分为渔业资源科学调查和渔业生产调7,其中,渔业资源科学调查易受调查成本、调查方法、调查范围、调查周期等因素影8-9,远洋渔业资源科学调查尤为突出。远洋渔业资源调查通常距离远、站点跨度大、费用高,优化调查设计方案,提升调查的精确度和准确度具有重要实践意10

西北太平洋区域是世界粮农组织(Food and Agriculture Organization of the United Nations,FAO)统计的所有区域中产量最高的区

11。澳洲鲐(Scomber australasicus)、日本鲭(Scomber japonicus)和远东拟沙丁鱼(Sardinops melanostictus)是西北太平洋重要经济鱼种,且日本鲭和远东拟沙丁鱼也是北太平洋渔业委员会(The North Pacific Fisheries Commission, NPFC)的优先管理物11。近年来,我国在西北太平洋海域开展了渔业资源科学调查,并分析了主要经济鱼种资源变动12-13,但对该海域科学调查站位设计方案尚未开展研究。本研究以西北太平洋澳洲鲐、日本鲭和远东拟沙丁鱼等的相对资源量以及总体相对资源量为调查目标,运用计算机模拟方法进行不同采样设计的抽样调查精度比较,以期为优化现行西北太平洋渔业资源调查方案和远洋渔业资源调查设计提供研究参考。

1 材料与方法

1.1 数据来源

数据来源于2021—2023年“淞航”号调查船在西北太平洋开展的渔业资源综合调查,调查时间为2021—2023年6—8月。调查海域为148°E~164°E、35°N~45°N,共设70个综合调查站点,调查站点栅格经纬度大小为2°×1°,其中2021年渔业资源拖网站点42个,2022年36个,2023年39个(图1)。“淞航”号船长85.00 m,宽14.96 m,平均拖速4~6 kn,每个站点拖网1~2 h。调查网具为四片式中层单囊拖网,网口周长434 m,网具总长度97.1 m。渔获物进行种类鉴别后冷冻保存,带回实验室进行生物学分析,按照拖网时间1 h和拖速5 kn进行数据标准化处理,得到各调查站点目标鱼种的相对资源量及总相对资源量(kg/h),以目标鱼种相对资源量和总相对资源量为调查目标进行模拟研究。

图1  2021—2023年西北太平洋调查站点空间分布

Fig.1  Sampling stations in Northwest Pacific Ocean from 2021 to 2023

1.2 模拟过程

本研究通过Kriging插值法获得调查目标的资源分布情况作为模拟抽样的总体,以70个综合调查站点为总体潜在调查站点(图1),采用简单随机抽样(Simple random sampling, SRS)、系统抽样(Systematic sampling, SYS)和分层随机抽样(Stratified sampling, StS)对“总体”进行重抽样(1 000次)。选择相对误差(Relative estimation error, REE)和相对偏差(Relative bias, RB)评价各抽样方法对调查目标的抽样估计准确度和偏离情况,使用REE的均值计算准确度变化率(Accuracy change rate, ACR),通过比较评价指标的数值区间分布来评价各采样设计的抽样表

14

1.2.1 抽样总体的获取

由于在2021年和2022年的调查中,未将获取的鲐进行日本鲭和澳洲鲐的区分,同时2021—2023年的总相对资源量和远东拟沙丁鱼的相对资源量空间分布特征较为相似(图2),因此将2023年4种调查目标的空间分布作为抽样总体。采用Kriging插值法获得澳洲鲐、日本鲭、远东拟沙丁鱼的相对资源量及总相对资源量在调查海域的“真实”分布情况,并假设该总体能够代表调查时该海域调查目标的“真实”分布,以此作为“真值”进行模拟抽样,比较不同采样设计的抽样精确度与准确度。

图2  2021—2023年总相对资源量和远东拟沙丁鱼相对资源量空间分布

Fig.2  Spatial distribution of total relative abundance and the relative abundance of Sardinops malanostictus in 2021-2023

1.2.2 采样设计

本研究比较了SRS、SYS和StS在估计各目标鱼种相对资源量及总相对资源量的表现。

SRS:每次抽样从70个潜在调查站点中随机抽取n个站点,进行不放回抽样,相应采样设计为Design 1。

SYS:按调查船航行路线对潜在调查站点进行编号,随机抽取站点m作为起点,再按照固定间隔k=70/n从总体中抽取n个站点(若k不为整数,则对k向下取整),进行环形系统抽样,相应采样设计为Design 2。

StS:根据各目标鱼种相对资源量及总相对资源量的空间分布特征,按照层间异质性大、层内异质性小的分层原则将调查水域划分为3层(相对资源量由高到低分别为A、B、C,图3),为了保证具有统计意义,每层最少分配3个采样站点。

图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。

1.2.3 评价指标

采用相对误差REE

15-16评价采样估计值的准确度和精密度,REE可以反映相对资源量估计值和“真实值”之间的差异情况。采用相对偏差RB17评价采样估计值的偏差,RB可以反映估计值与“真实值”的偏离情况,可以反映采样估计值对总体均值的高估和低估情况。REE值越小,RB值越趋近于0,采样设计效果越16

REE与RB的计算公式:

REE=i=1N(Viestimated-Vtrue)2NVtrue×100% (1)
RB=i=1NViestimatedN-VtrueVtrue×100% (2)

式中:Viestimated为采样设计第i次模拟抽样调查的估计值均值;Vtrue为目标鱼种相对资源量和总相对资源量“真实”值的均值;N为模拟次数,本研究中N=1 000。

基于1 000次重抽样的REE均值来计算准确度变化率(Accuracy change rate,ACR),以此作为衡量REE变化量与站点数量变化的比

15。由于传统ACR计算方法为使用某一站点数对应的准确度统计量值与上一站点数对应的准确度统计量值进行计算以确定准确度的变化,考虑到系统抽样的REE存在随站点数量增加而呈现波动的情18-19,导致在某些站点数量虽然ACR的值较大,但实际其准确度可能并未显著增加。因此,本研究ACR的计算方式在参考相关研究的基础上进行了优1520

ACR=REEn-REEminSmin-Sn (3)

式中:SnSmin为不同站点数,Sn大于Smin,且Smin为站点数在[3,n-1]区间内REE最小值对应的站点数;REEn、REEmin为不同站点数对应的1 000次REE均值,选取ACR值为0.2作为确定REE显著降低的标准线。

1.2.4 抽样过程

(1)对抽样总体分别进行不同采样设计的重抽样,站点数量为3~50(分层随机抽样为9~50),进行连续抽样。根据模拟数据计算不同采样设计下的各评价指标数值。(2)重复上述过程1 000次,计算不同采样设计的REE和RB,以此比较不同采样设计的表现。(3)在1 000次循环抽样结果基础上,使用1 000次REE的均值计算ACR,以确定各采样设计最优站点数。

2 结果

2.1 不同采样设计的相对误差

总体来看,简单随机抽样和分层随机抽样对4种调查目标总体估计的REE变化趋势一致,均随站点数量的增加而降低,且2种分层随机抽样的REE小于简单随机抽样(图4)。总体上,系统抽样对4种调查目标总体估计的REE呈降低趋势,但在个别抽样站点数量(如23、35)存在陡增现象,当抽样间隔k相同时,站点数量越多,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值趋于一致(图4a)。当以日本鲭相对资源量为调查目标时,Design 4的REE值最小,当抽样站点数达到20后,REE值变化幅度降低;Design 3的REE值呈阶段式下降趋势,当抽样站点数达到30后,REE值略高于Design 1(图4b)。当以远东拟沙丁鱼相对资源量和总相对资源量为调查目标时,分层随机抽样的REE值均较小,且Design 4的REE值小于Design 3(图4c、4d)。

2.2 不同采样设计的相对偏差

当以澳洲鲐、远东拟沙丁鱼的相对资源量和总相对资源量为调查目标时,4种采样设计的RB值在-5%~5%的数值范围内围绕0值上下波动,均为无偏估计(图5)。除Design 2外,其余3种采样设计的RB值分布区间随站点数的增加逐渐减小,Design2的RB值则在-5%~5%的数值范围内波动变化(图5a、5c、5d)。当以日本鲭的相对资源量为调查目标时,Design 1和Design 2的RB值在-10%~10%的数值范围内围绕0值上下波动,为无偏估计,分层随机抽样的RB值均在8%~16%的数值范围内围绕12%上下波动,表明分层随机抽样高估了日本鲭的相对资源量,且随站点数量的增加RB值分布区间逐渐减小,Design 4的RB值总体小于Design 3(图5b)。

图5  不同采样设计对各调查目标抽样相对偏差(RB)变化趋势

Fig.5  RB values of the different sampling designs based on different objectives

2.3 准确度变化率

Design 1的ACR均大于0,并随抽样站点数量的增加而降低(图6a、6b、6c、6d),表明Design 1的抽样精度随抽样站点数量的增加而提升;Design 1在估计澳洲鲐相对资源量且抽样站点数量为30时,ACR达到标准值0.2,在估计其他调查目标时,ACR均未达到0.2。Design 2的ACR在抽样站点数量达到35时,ACR达到0.2,并且在抽样站点大于35后ACR均为负值(图6e、6f、6g、6h),表明在站点数达到35后,继续增加站点数量并不能达到更好的抽样效果。Design 3在估计澳洲鲐相对资源量和总相对资源量时,在站点数为20和37时,ACR达到0.2,表明在站点数量达到20和37时,REE下降趋势减缓(图6i、6l);在估计日本鲭和远东拟沙丁鱼相对资源量时,ACR值趋近于0值(图6j、6k),表明随抽样站点数量增加Design 3的抽样精度呈现平滑上升趋势。除澳洲鲐外,Design 4在估计日本鲭、远东拟沙丁鱼相对资源量和总相对资源量时,ACR随抽样站点数增加呈下降趋势(图6m、6n、6o、6p),在抽样站点数量分别为26、18和20时ACR首次达到0.2,表明在抽样站点数量分别达到26、18和20时Design 4的抽样精度上升趋势减缓。

图6  不同采样设计的ACR 值随样本量的变化

Fig.6  Changes in accuracy change rate (ACR) of different sampling designs with different sampling stations

2.4 最优采样设计

当对不同调查目标模拟抽样时,Design 4的REE值较小且最为稳定,最优采样设计可选择为Design 4。当Design 4对不同的调查目标进行模拟抽样时,最适调查站点数量存在不同(表1),当以ACR=0.2为评价标准并综合考虑REE变化规律时,在保证一定采样精确度和准确度的前提下,最适调查站点数量可选择为28。

表1  Design 4的REE范围及最优站点数
Tab.1  REE range of Design 4 and the optimal sampling stations for each survey objective

调查目标

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

3 讨论

本研究以多个鱼种相对资源量以及总相对资源量为调查目标,比较了4种采样设计在不同调查目标下的抽样表现。研究发现,相比简单随机抽样和系统抽样,分层随机抽样的模拟抽样表现较好。简单随机抽样是最简单的抽样方法,其他抽样方法在某种程度上包含简单随机抽样。本研究中简单随机抽样的REE、RB随抽样站点数的增加而降低,但抽样精度相对不高,这与相关研

21结果一致。由于大洋环境变化剧烈,鱼类分布受到温度、海流等环境因素的影响而呈现集群分19,简单随机抽样的完全随机性会导致估计误差较大,因此在多数情况下并不适22

系统抽样可以覆盖调查区域范围,并能兼顾随机性,常被作为初始调查方案,但在某些调查努力量范围内,其调查精度存在波

151720,本研究系统抽样也发现相同结果。系统抽样调查精度波动与抽样间隔存在关系,系统抽样采用连续抽样,抽样间隔k采取向下取整,当k减小时,REE出现急剧增加的现象;当k不变时REE随站点数量的增加呈现下降趋势。例如当抽样站点数达到23和35后,系统抽样REE出现陡增的波动,抽样站点数23和35则分别对应k为3和2时最大站点数。鱼类会根据环境变化选择合适的栖息地导致其时空分布的不均23-24,当k减小,系统抽样的空间覆盖率下降,抽样站点分布与鱼类分布存在空间错位,导致其抽样精度下降;当k相同时,随站点数量的增加,站点空间覆盖率增加,其抽样精度升高。

分层随机抽样的估计效果比简单随机抽样和系统抽样精度高,这与相关研

25-26结果一致。两种分层随机抽样在不同调查目标下的估计效果并不相同,根据一般抽样理论,基于层面积大小分配方式的采样设计往往表现更好。在本研究中,日本鲭、远东拟沙丁鱼相对资源量和总相对资源量各站点差异较大,总体方差均较大。以日本鲭、远东拟沙丁鱼相对资源量和总相对资源量为调查目标时,基于层资源量大小分配方式的Design 4表现更好,这可能是由层内各站点调查目标相对资源量的差异所导致。对于日本鲭、远东拟沙丁鱼相对资源量和总相对资源量来说,A层和B层的相对资源量方差较大(表2),当以相对资源量大小分配站点时(Design 4),抽样误差相对Design 3更低。澳洲鲐相对资源量在A层、B层、C层各站点间的差异很小(表2),满足层内异质性小的分层原则,因而基于层内站点数量的分配方式(Design 3)表现更好。

表2  各层不同调查目标值的方差
Tab.2  Variance of different sampling objectives in the different stratums
调查目标 Survey objectivesA层 Stratum-AB层 Stratum-BC层 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对日本鲭存在高估现象,这可能与日本鲭的空间分布较为集中有

27-28。当以日本鲭为调查对象时,A层的潜在站点数量占比为14.29%,层内潜在站点总资源量占日本鲭总资源量65.64%,这就造成A层抽样站点数量对估计结果的影响较大。Design 4按照各层资源量比例进行站点数量分配,A层分配更多的站点使调查精度更高,由于日本鲭分布相对集中,A层分配更多的站点导致分层随机抽样对日本鲭相对资源量存在高估。因此,分层随机抽样的分层方式和站点分配原则还需进一步优29。总体来看,Design 4的估计效果稳定且精度较高,但对于不同调查目标,Design 4最适的抽样调查站点数量存在不同(表1)。当以澳洲鲐、日本鲭、远东拟沙丁鱼、渔业资源总体为调查对象时,最适调查站点数量分别为28、26、18、20,最大值为28。在调查站点数量为28时,Design 4对4种调查目标的调查精度均相对较好。当以多目标进行采样设计比较时,很少出现某一种采样方案对各调查目标的估计皆最20,虽然Design 4对日本鲭相对资源量存在一定程度的高估,但能够满足大部分调查目标的要求,因此,Design 4可作为西北太平洋渔业资源调查多目标的最优采样设计。

4 不足与展望

在各目标鱼种及总相对资源量真值获取上,本研究使用了Kriging插值法,并未考虑环境及其他因素对物种分布的影响,存在一定的局限性。联合物种分布模型和广义加性模型可以分析不同物种之间的相互作用及对其分布产生的影响,并提供更准确的预测结

30-31,在未来研究中可以考虑使用这些方法获得“真值”并进行不同采样设计比较;在近海海域及江湖流域通常以影响物种分布的环境因素对研究区域进行分2932,但大洋环境因素年际变化33-34,温度、盐度均易受到洋流及热带气旋等影35,如何更合理分层也还需进一步探讨。

参考文献

1

柳昊张文超李文甲. 黄海中部夏季鳀鱼卵、仔稚鱼分布与早期生长特征[J]. 上海海洋大学学报2024331): 172-185. [百度学术] 

LIU HZHANG W CLI W Jet al. Distribution and early growth characteristics of anchovy eggs, larvae and juveniles in the central Yellow Sea in summer[J]. Journal of Shanghai Ocean University2024331): 172-185. [百度学术] 

2

王淼娣王雪辉杜飞雁. 基于LBB方法估算北部湾竹䇲鱼种群参数[J]. 上海海洋大学学报2022311): 212-222. [百度学术] 

WANG M DWANG X HDU F Yet al. Estimation of the population parameters of Trachurus japonicus in the Beibu Gulf based on the LBB method[J]. Journal of Shanghai Ocean University2022311): 212-222. [百度学术] 

3

鲁红月汪金涛张畅. 基于产卵场海表面温度的太平洋褶柔鱼冬生群资源评估[J]. 上海海洋大学学报2022313): 650-657. [百度学术] 

LU H YWANG J TZHANG 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 University2022313): 650-657. [百度学术] 

4

孙菁煜戴小杰朱江峰. 淀山湖鱼类多样性分析[J]. 上海水产大学学报2007165): 454-459. [百度学术] 

SUN J YDAI X JZHU J Fet al. Analysis of the fish species diversity in Dianshan Lake[J]. Journal of Shanghai Fisheries University2007165): 454-459. [百度学术] 

5

ZHAO JYANG K EMA J. Optimization of sampling effort for different fishery groups in the Yangtze river estuary, China[J]. Marine and Coastal Fisheries2022144): e10214. [百度学术] 

6

STARR R MCARR MMALONE Det al. Complementary sampling methods to inform ecosystem-based management of nearshore fisheries[J]. Marine and Coastal Fisheries201021): 159-179. [百度学术] 

7

HOWARD R ACIANNELLI LWAKEFIELD W Wet al. Comparing fishery-independent and fishery-dependent data for analysis of the distributions of Oregon shelf groundfishes[J]. Fisheries Research2023258106553. [百度学术] 

8

SCHEIRER KCHEN YWILSON C. A comparative study of American lobster fishery sea and port sampling programs in Maine: 1998-2000[J]. Fisheries Research2004681/3): 343-350. [百度学术] 

9

ZHAO JCAO JTIAN S Qet al. Evaluating sampling designs for demersal fish communities[J]. Sustainability2018108): 2585. [百度学术] 

10

赵静章守宇林军. 不同采样设计评估鱼类群落效果比较[J]. 应用生态学报2014254): 1181-1187. [百度学术] 

ZHAO JZHANG S YLIN Jet al. A comparative study of different sampling designs in fish community estimation[J]. Chinese Journal of Applied Ecology2014254): 1181-1187. [百度学术] 

11

SHI Y CZHANG X MHE Y Ret al. Stock assessment using length-based Bayesian evaluation method for three small pelagic species in the northwest Pacific ocean[J]. Frontiers in Marine Science20229775180. [百度学术] 

12

刘思源张衡杨超. ENSO事件下西北太平洋远东拟沙丁鱼和日本鲭栖息地协同变化特征[J]. 海洋学报2024461): 39-52. [百度学术] 

LIU S YZHANG HYANG Cet al. Exhibit covariation characteristics in the habitat changes of Sardinops melanostictus and Scomber japonicus in the northwestern Pacific Ocean under ENSO event[J]. Haiyang Xuebao2024461): 39-52. [百度学术] 

13

蔡锴周雨霏麻秋云. 西北太平洋鲐的资源养护和管理研究[J]. 中国水产科学2024314): 439-453. [百度学术] 

CAI KZHOU Y FMA Q Yet al. Research on the conservation and management for chub mackerel (Scomber japonicus) in the Northwest Pacific[J]. Journal of Fishery Sciences of China2024314): 439-453. [百度学术] 

14

马金. 基于计算机模拟的河口渔业资源抽样调查设计比较与优化——以长江口为例[D]. 上海上海海洋大学2020. [百度学术] 

MA J. Simulation-based evaluation of sampling design for different fishery groups in estuary: a case study in the Yangtze estuary[D]. ShanghaiShanghai Ocean University2020. [百度学术] 

15

LIU YCHEN YCHENG J H. A comparative study of optimization methods and conventional methods for sampling design in fishery-independent surveys[J]. ICES Journal of Marine Science2009669): 1873-1882. [百度学术] 

16

CHEN Y. A Monte Carlo study on impacts of the size of subsample catch on estimation of fish stock parameters[J]. Fisheries Research1996263/4): 207-223. [百度学术] 

17

CAO JCHEN YCHANG J Het al. An evaluation of an inshore bottom trawl survey design for American lobster (Homarus americanus) using computer simulations[J]. Journal of Northwest Atlantic Fishery Science20144627-39. [百度学术] 

18

MA JTIAN S QGAO C Xet al. Evaluation of sampling designs for different fishery groups in the Yangtze River estuary, China[J]. Regional Studies in Marine Science202038101373. [百度学术] 

19

MA JLI BZHAO Jet al. Environmental influences on the spatio-temporal distribution of Coilia nasus in the Yangtze River estuary[J]. Journal of Applied Ichthyology2020363): 315-325. [百度学术] 

20

张国晟王晶张崇良. 基于多目标鱼种资源量指数估计的不同采样设计比较[J]. 水产学报2021455): 700-715. [百度学术] 

ZHANG G SWANG JZHANG C Let al. Comparison of sampling designs of fishery-independent survey in estimating abundance indices of multiple target species[J]. Journal of Fisheries of China2021455): 700-715. [百度学术] 

21

金勇进杜子芳蒋妍. 抽样技术[M]. 2版. 北京中国人民大学出版社2008. [百度学术] 

JIN Y JDU Z FJIANG Y. Sampling technique[M]. 2nd ed. BeijingChina Renmin University Press2008. [百度学术] 

22

刘勇. 渔业资源评估抽样调查方法的理论探讨与研究[D]. 上海华东师范大学2012. [百度学术] 

LIU Y. Theoretical study on the sampling methods of survey for fishery stock estimation[D]. ShanghaiEast China Normal University2012. [百度学术] 

23

刘思源张衡杨超. 基于最大熵模型的西北太平洋远东拟沙丁鱼和日本鲭栖息地差异[J]. 上海海洋大学学报2023324): 806-817. [百度学术] 

LIU S YZHANG HYANG Cet 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 University2023324): 806-817. [百度学术] 

24

周茉方星楠余为. 厄尔尼诺和拉尼娜事件下西北太平洋柔鱼栖息地时空分布差异[J]. 上海海洋大学学报2022314): 984-993. [百度学术] 

ZHOU MFANG X NYU Wet 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 University2022314): 984-993. [百度学术] 

25

潘邵媛王学昉田思泉. 不同采样设计对长江口海洋环境监测效果影响的评价[J]. 海洋环境科学2021403): 435-441. [百度学术] 

PAN S YWANG X FTIAN S Qet al. Evaluation of the effect of different sampling designs on marine environmental monitoring in the Yangtze River Estuary[J]. Marine Environmental Science2021403): 435-441. [百度学术] 

26

潘邵媛王学昉田思泉. 长江口中华鲟保护区海洋环境监测浮标站点的优化设计[J]. 海洋学报2021434): 55-64. [百度学术] 

PAN S YWANG X FTIAN S Qet al. The design of the stations of marine environmental monitoring buoys in the Chinese sturgeon nature reserve in the Changjiang River Estuary[J]. Haiyang Xuebao2021434): 55-64. [百度学术] 

27

刘思源张衡杨超. 西北太平洋远东拟沙丁鱼与日本鲭种群动态特征及其与环境因子关系研究进展[J]. 大连海洋大学学报2023382): 357-368. [百度学术] 

LIU S YZHANG HYANG Cet 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 University2023382): 357-368. [百度学术] 

28

刘彦琳. 基于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]. ShanghaiShanghai Ocean University2023. [百度学术] 

29

WANG JXU B DZHANG C Let al. Evaluation of alternative stratifications for a stratified random fishery-independent survey[J]. Fisheries Research2018207150-159. [百度学术] 

30

ZHANG C LXU B DXUE Yet al. Evaluating multispecies survey designs using a joint species distribution model[J]. Aquaculture and Fisheries202053): 156-162. [百度学术] 

31

张涛赵天亚栾静. 基于物种分布模型分析环境因子对海州湾偶见种资源分布的影响[J]. 海洋学报2023457): 69-78. [百度学术] 

ZHANG TZHAO T YLUAN Jet 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 Xuebao2023457): 69-78. [百度学术] 

32

XU B DREN Y PCHEN Yet al. Optimization of stratification scheme for a fishery-independent survey with multiple objectives[J]. Acta Oceanologica Sinica20153412): 154-169. [百度学术] 

33

王雅萌陈芃陈新军. 气候变化下西北太平洋大海洋生态系海表面温度特征分析[J]. 上海海洋大学学报2021305): 874-883. [百度学术] 

WANG Y MCHEN PCHEN X J. Characteristics of sea surface temperature for large marine major ecosystems in the northwest Pacific under climate change[J]. Journal of Shanghai Ocean University2021305): 874-883. [百度学术] 

34

余为温健. 洪堡洋流系统内重要渔业资源对气候和海洋环境变化的响应研究进展[J]. 上海海洋大学学报2022313): 620-630. [百度学术] 

YU WWEN J. Review on the response of important fishery resources to the climatic and environmental variability in the Humboldt Current System[J]. Journal of Shanghai Ocean University2022313): 620-630. [百度学术] 

35

谢晓丽. 基于Argo资料的西北太平洋上层温盐对热带气旋的三维响应特征分析[D]. 湛江广东海洋大学2022. [百度学术] 

XIE X L. Analysis of three-dimensional temperature and salinity response to tropical cyclones in the northwest Pacific observed by Argo floats[D]. ZhanjingGuangdong Ocean University2022. [百度学术]