Ecological carrying capacity of six species of stock enhancement in Pearl River estuary based on Ecopath model
-
摘要:
增殖放流是渔业资源养护的重要方式。放流前对放流海域进行生态容纳量评估，有计划地实施增殖放流活动，可避免对原有生态系统造成破坏。文章根据2016年珠江口渔业资源数据，构建了由29个功能组组成的基于珠江口的生态系统通道(Ecopath)模型，利用该模型分析了生态系统的总体特征、食物网结构与混合营养效应，估算了适宜于该水域的6种不同增殖放流种类的生态容纳量。结果表明，功能组营养级范围为1~4.2级，6种适宜放流种类营养级介于2.2~3.7，最高营养级功能组为哺乳动物，系统总流量9 092.447 t·(km2·a)–1，系统总能量转化效率12.23%，连接指数0.370，系统杂食指数0.287。食物链通道主要有2类，以碎屑食物链为主。花鲈(Lateolabrax japonicus)、黑鲷(Acanthopagrus schlegelii)、黄鳍鲷(A. latus)、长毛对虾(Penaeus penicillatus)、墨吉对虾(P. monodon)和波纹巴非蛤(Paphia undulata)的最大容纳量分别为0.094 t·km–2、0.500 t·km–2、0.650 t·km–2、1.580 t·km–2、1.610 t·km–2和75.870 t·km–2。
Abstract:Stock enhancement is important for conservation of exploitable resources. Assessment of ecological capacity of the discharged sea area before discharge and planned implementation of proliferation and discharge activities can avoid damage to the original ecosystem. Based on the survey data of the fishery resources of Pearl River estuary in 2016, an Ecopath mass-balance model of the Pearl River estuary ecosystem consisting of 29 functional groups was constructed. Thus, the overall characteristics of the ecosystem, the food-web structure and the mixed trophic impact were analyzed. The ecological carrying capacities of six species of stock enhancement which were suitable for this area were estimated. The results show that the fractional trophic levels ranged from 1 to 4.2, and those of the six species ranged from 2.2 to 3.7, with marine mammals occupying the highest trophic level. The total system throughput of the ecosystem was 9 092.447 t·(km2·a)–1, and the total energy transfer efficiency was 12.23%. The connectance index and system omnivory index were 0.370 and 0.287, respectively. There were two main channels of food chain in the ecsystem, and the energy flow was dominated by grazing food chain. The ecological carrying capacities of biomass for Lateolabrax japonicus, Acanthopagrus schlegelii, A. latus, Penaeus penicillatus, P. monodon and Paphia undulata were 0.094 t·km–2, 0.500 t·km–2, 0.650 t·km–2, 1.580 t·km–2, 1.610 t·km–2 and 75.870 t·km–2, respectively.
-
Keywords:
- Ecopath model /
- ecological carrying capactiy /
- enhancement /
- Pearl River estuary
-
深水金线鱼 (Nemipterus bathybius) 属鲈形目、金线鱼科、金线鱼属，为暖水性底层经济鱼类，在我国主要分布于南海北部和台湾周边海域，是底拖网渔业的主要渔获物之一[1-2]。体型较大的深水金线鱼个体通常以甲壳类动物、小型鱼类和头足类为食，较小的个体以桡足类、介形虫和片脚类为食。深水金线鱼生长速度快且性成熟较早，具有分批产卵、生命周期短等生物学特征，可以承受一定的捕捞压力[3]。随着南海北部近海渔业资源捕捞压力的持续上升，深水金线鱼资源状况也开始呈现过度开发现象，近年来其渔获物多为幼鱼，种群结构以低龄化个体为主[4-6]。南海幅员辽阔，但我国在南海大部分的捕捞活动和产量都集中在北部陆架区。北部陆架区是重要的经济鱼类产卵场，也是重要的渔场[7]。改革开放以来，随着沿海地区经济的快速发展，人民群众对海产品的需求量持续上升，海洋渔船数量和捕捞努力量不断增加[8]，导致了南海北部近海渔业资源的衰退，特别是典型的半封闭海湾[9-10]。Zhang等[11]对北部湾9种主要经济鱼类进行资源评估，结果表明9种鱼类在不同时期均遭受了过度捕捞。《中国渔业统计年鉴》的全国海洋捕捞鱼类产量数据显示，2018—2022年南海金线鱼类产量平均为31.4万t，约占南海海洋鱼类总捕捞产量的16.2%，该类群经济价值高，其中深水金线鱼是主要的优势种类之一。以往对南海北部深水金线鱼种群的资源评估多针对生长、死亡和性成熟等生物学特征，且评估数据多为2010年以前的。对北部湾深水金线鱼种群参数的评估研究[4,11]表明，其生物学参数年际变化明显，种群个体虽生长速度有所增加，但小型化现象突出，捕捞个体多为幼鱼。
渔业资源评估是制定合理管理策略的基础，科学有效的评估结果有助于实现渔业资源的可持续开发。模型是进行渔业资源评估的重要工具，数据缺乏条件下的资源评估模型主要分为单鱼种评估方法、多鱼种评估方法和生态系统模型，目前国际上常用的是单鱼种的数据缺乏方法[12]。相对于其他渔业资源数据，长度数据的收集操作简便、结果读取直接、采集成本低，因此基于长度数据的渔业资源评估方法近些年发展较为迅速。本研究根据2014—2019年南海北部底拖网调查的渔业生物学资料，使用基于长度的贝叶斯生物量评估 (Length-based Bayesian biomass estimation method, LBB) 和基于长度的繁殖潜力比 (Length-based spawning potential ratio, LBSPR) 2种模型，对南海北部深水金线鱼种群资源现状进行评估，以期为该种群的科学管理和可持续利用提供技术支撑。
1. 材料与方法
1.1 实验材料
本研究所用样品数据均来自2014—2019年中国水产科学研究院南海水产研究所在南海北部海域底拖网调查资料，调查站点根据渔业资源按水深分布[13]的特点设计(图1)。所有航次的调查采样与生物学测定均按照GB/T 12763.3—2007《海洋调查规范》进行。采集的深水金线鱼样本冷冻保存，带回实验室进行生物学测量和鉴定，包括体长、体质量、性成熟度、胃饱满度等。每个站位深水金线鱼样本多于50尾时随机选取50尾进行测定，样本不足50尾时则全部测定，体长和体质量的测定分别精确到1 mm和1 g。12个航次共测定深水金线鱼3 059尾，样品信息见表1。
![]()
图 1 南海北部近海研究区域及站点分布Fig. 1 Study area and station distribution in offshore northern South China Sea表 1 南海北部深水金线鱼的样品信息Table 1 Sample information of N. bathybius in northern South China Sea年份
Year月份
Month样本数
Sampling
number体长 Body length/mm 范围
Range均值±标准差
Mean±Standard
deviation2014 7 293 93~193 121.2±16.3 10 225 77~184 119.9±21.6 2015 1 120 73~172 121.8±18.7 4 330 75~189 111.4±19.2 2016 7 304 84~200 121.8±19.4 10 78 74~178 122.1±23.8 2017 1 71 72~176 109.3±18.7 4 197 72~158 115.0±15.3 2018 3 220 75~176 113.8±19.0 8 87 78~182 126.0±24.2 2019 3 524 73~176 106.6±16.8 9 610 65~192 117.3±19.0 1.2 方法
1.2.1 体长-体质量关系
采用幂指数拟合深水金线鱼体长与体质量关系，计算公式为：
$$ W {\mathrm{=}}a{L}^{b} $$ (1) 式中：W为体质量 (g)；L为体长 (mm)；a为生长的肥满度指数；b为幂指数。
1.2.2 生长、死亡参数估计
深水金线鱼的生长关系可通过von Bertalanffy生长方程[4]进行拟合：
$$L_{{t}}{\rm{=}}L_{\rm {inf }}\left\{1{\text{−}}{\mathrm{exp}} \left[{\text{−}}K\left(t{\text{−}}t_0\right)\right]\right\} $$ (2) 式中：Lt为t龄的体长；Linf为渐近体长；K为生长系数；t0为体长为0时的理论年龄。本研究使用R软件包TropFishR[14]中的ELEFAN (Electronic length frequency analysis) 技术估算深水金线鱼的渐近体长 (Linf) 和生长系数 (K)。参数t0采用以下Pauly经验方程[15]进行计算：
$$ \lg\left({\text{−}}t_0\right){\mathrm{=}}{\text{−}}0.392\,2{\text{−}}0.275 \lg L_{\rm {inf }}{\text{−}}1.038 \lg K $$ (3) 根据以下经验公式[16]计算生长性能指数 (Φ)：
$$ \Phi{\mathrm{=}}\lg K{\mathrm{+}}2 \lg L_{\rm {inf }} $$ (4) $$ M_1{\mathrm{=}}\mathrm{e}^{{\text{−}}0.015\,2{\mathrm{+}}0.654\,3 \times \ln ({K}){\text{−}}0.279 \times \ln \left(L_{\text {inf }}\right){\mathrm{+}}0.463\,4 \times \ln (T)}$$ (5) $$M_2{\mathrm{=}}4.118 \times K^{0.73} \times L_{\text {{\mathrm{inf}} }}^{{\text{−}}0.33} $$ (6) 式中：T为环境平均温度，本研究采用27 ℃[19]；取M1、M2的平均值，以获得南海北部深水金线鱼自然死亡率(M)相对准确的值。
1.2.3 50%性成熟体长
统计深水金线鱼繁殖期雌性个体的性腺发育状况，按照I—VI期标准进行划分，性腺发育达到IV期或以上，则认为该个体为性成熟个体[4]。50%性成熟体长(L50)采用最小二乘法 (Least square method, LSM)拟合以下逻辑斯蒂模型得到[20]：
$$ P_i{\mathrm{=}}\frac{1}{1{\mathrm{+}}{\mathrm{exp}} \left[{\text{−}}r\left(L_i{\text{−}}L_{50}\right)\right]} $$ (7) 式中：Pi为体长组Li对应的成熟个体比例；r为模型参数。95%性成熟体长(L95)由1.1×L50估算[21]。
1.3 基于长度数据的评估模型
1.3.1 基于长度的贝叶斯生物量模型
LBB是近年来新开发的一种根据长度-频率数据来评估鱼类种群开发状况的方法，该方法使用贝叶斯蒙特卡罗马尔科夫链 (MCMC)来计算死亡率参数和相对种群数量，适用于可以终生持续生长的物种[22]，该方法可以估算目标物种的渐近体长 (Linf)、50%个体可被渔具捕获的体长值 (Lc)、最适开捕体长 (Lc_opt)、相对自然死亡率 (M/K)、捕捞死亡率相对于自然死亡率 (F/M)、已经开发的生物量相对于未开发的生物量 (B/B0) 和相对种群大小 (B/BMSY，其中MSY为最大可持续产量)。LBB方法只需要长度-频率数据便可进行计算，但也可以通过设置M/K、Linf、L50等先验参数，以获得更准确的结果[23] (表2)。
表 2 LBB和LBSPR方法所需数据目录Table 2 Data for LBB and LBSPR methods参数
Parameter估算值
Estimated value基于长度的繁殖潜力比模型
Length-based spawning potential ratio
(LBSPR) assessment model基于长度的贝叶斯估计模型
Length-based Bayesian biomass
(LBB) assessment model长度-频率数据
Length–frequency data/ √ √ 渐近体长 (Linf)
Asymptotic length/cm23.7 cm √ 〇 生长系数 (K)
Growth coefficient0.326 〇 自然死亡率 (M)
Natural mortality0.762 〇 相对自然死亡率M/K 2.33 √ 〇 50%性成熟体长 (L50)
Body length at 50% sexual maturity/cm11.76 √ 〇 95%性成熟体长 (L95)
Body length at 95% sexual maturity/cm12.94 √ 渐近体长的变异系数
Coefficient of variation of Linf0.1 √ 长度-体质量关系参数 (a)
Length-body mass relationship parameter6.854 0×10−5 〇 长度-体质量关系参数 (b)
Length-body mass relationship parameter2.797 4 〇 注：“√” 代表必需数据；“〇”代表可选数据。 Note: "√" represents required data; "〇" represents optional data. 根据B/BMSY和Lc/Lc_opt的评估结果，确定深水金线鱼的种群开发现状。如果一个种群的B/BMSY 1，可认为该种群处于健康状态；当0.8 ≤ B/BMSY < 1时，认为该种群处于轻度过度开发状态；当0.5 B/BMSY < 0.8时，认为该种群正处于过度捕捞状态；当0.2≤B/BMSY < 0.5时，认为该种群处于严重过度捕捞状态；当B/BMSY < 0.2时，认为该种群资源已经崩溃[24]。如果估计的Lc/Lc_opt< 1，认为种群遭受了“生长过度捕捞”(当鱼的平均捕捞尺寸小于每次捕捞可产生最大产量的尺寸时)，其中Lc_opt为最适开捕长度，表示可保持种群处于健康状态且捕捞量接近MSY的体长值。LBB方法在2018年由Forese等[25]提出，本研究使用的模型R语言代码(LBB-33a)来自基尔亥姆霍兹海洋研究中心 (http://oceanrep.geomar.de/44832/)。
1.3.2 基于长度的繁殖潜力比模型
在基于长度-频率的数据缺乏分析法中，繁殖潜力比 (Spawning potential ratio, SPR) 是最常用的参考点之一，可以为数据缺乏的渔业管理提供具有决策性的建议[26]，SPR定义为在一定捕捞压力下(F>0)，当前种群与未被捕捞种群的繁殖潜力比值[27]。繁殖潜力指在不考虑密度制约的情况下，补充量群体的潜在产卵量。LBSPR是基于SPR这一生物学参考点建立的，广泛应用于数据缺乏种群的资源评估[21,28-29]。LBSPR假设被捕捞鱼类种群的SPR是相对捕捞死亡率 (F/M) 和2个生活史比率 (M/K、L50/Linf) 的函数，与LBB不同的是，LBSPR的运行除长度-频率外还需要其他参数 (即M/K、Linf、L50、L95，如表2所示)，这些参数的估计可通过R软件包TropFishR、M的经验公式和逻辑斯蒂模型计算获得。
在LBSPR模型的运行过程中，使用最大似然估计法来估算50%和95%选择性下的平均体长($S_{L_{50}} $和$S_{L_{95}} $)和相对捕捞死亡率(F/M)，若SPR≤0.2，说明该种群繁殖潜力比过低，种群资源濒临崩溃；当0.2<SPR≤0.35时，说明种群处于充分开发状态；当0.35<SPR<0.4时，则说明种群处于MSY水平[21,26,29]。LBSPR方法在2015年由Hordyk等[30]提出，本研究使用R语言LBSPR软件包完成(https://mirrors.pku.edu.cn/CRAN/web/packages/LBSPR/index.html)。
1.3.3 模型不确定性分析
本文对2个评估模型的Linf、M/K、L50 3种参数进行不确定性分析，以判断不同先验参数在不同区间的敏感性。Linf、M/K和L50根据LBB对3种参数所估算的95%置信区间的数值进行上下浮动，来保证先验信息相对准确，其中Linf分为6组 (数值如表3中第2—第7组所示)，体长数值间隔为2 mm；L50分为6组 (数值如表3中第8—第13组所示)，体长数值间隔为2 mm；M/K分为6组 (数值如表3中第14—第19组所示)，数值间隔为0.02。
表 3 LBB和LBSPR在不同先验参数下的评估结果Table 3 Estimated results of LBB and LBSPR parameters with different priors组号
Group No.先验参数
Prior parameterLBSPR结果
Results of LBSPRLBB结果
Results of LBBLinf/mm L50/mm M/K F/M SPR B/BMSY F/M Lc/Lc_opt 1 237 117.6 2.33 1.71 0.19 0.89 1.00 0.85 2 232 1.73 0.20 0.99 0.90 0.87 3 234 1.76 0.20 0.94 0.97 0.86 4 236 1.71 0.19 0.90 1.00 0.85 5 238 1.84 0.18 0.87 1.10 0.84 6 240 1.73 0.19 0.84 1.10 0.83 7 242 1.94 0.17 0.81 1.10 0.82 8 237 112.6 2.33 1.71 0.21 0.89 1.00 0.85 9 114.6 1.71 0.21 0.89 1.00 0.85 10 116.6 1.71 0.20 0.89 1.00 0.85 11 118.6 1.71 0.19 0.89 1.00 0.85 12 120.6 1.71 0.18 0.89 1.00 0.85 13 122.6 1.71 0.17 0.89 1.00 0.85 14 237 117.6 2.28 1.77 0.19 0.84 1.10 0.83 15 2.30 1.75 0.19 0.87 1.00 0.84 16 2.32 1.73 0.19 0.87 1.00 0.84 17 2.34 1.70 0.20 0.89 1.00 0.84 18 2.36 1.68 0.20 0.92 1.00 0.86 19 2.38 1.66 0.20 0.92 0.99 0.86 2. 结果
2.1 参数估计
2014—2019年采集的深水金线鱼样品体长范围为6.5~20.0 cm，体质量范围为9.3~202.0 g。雄性、雌性和整体的体长-体质量关系 (图2) 分别为：
![]()
图 2 南海北部雄性、雌性、整体深水金线鱼体长-体质量关系Fig. 2 Length-body mass relationship among male, female and all N. bathybius in northern South China Sea$$ \text { 雄性}： W{\rm{=}}6.911\,8 \times 10^{{\text{−}}5} L^{2.794\,8}\left(r^2{\rm{=}}0.920\,1\right) $$ (8) $$ \text { 雌性 } ：W{\mathrm{=}}8.211\,0 \times 10^{{\text{−}}5} L^{2.761\,0}\left(r^2{\mathrm{=}}0.911\,5\right) $$ (9) $$ \text { 整体}： W{\mathrm{=}}6.854\,0 \times 10^{{\text{−}}5} L^{2.797\,4}\left(r^2{\mathrm{=}}0.920\,3\right) $$ (10) 在R软件包TropFishR中使用ELEFAN方法提供的von Bertalanffy生长参数估计，评估的Linf为23.7 cm，K为0.326，t0为 −0.54，估算的生长性能指数(Φ)为2.26。von Bertalanffy生长函数为：
$$ L_t{\mathrm{=}}23.7 \times\{{\text{−}}{\mathrm{exp}} [{\text{−}}0.33 \times(t{\mathrm{+}}0.54)]\} $$ (11) 使用2种经验公式估算的深水金线鱼自然死亡率M分别为0.908和0.616，本研究使用2个结果的平均值0.762。因此，LBB和LBSPR的先验信息参数M/K为2.33 (表2)。
根据非线性回归拟合的LSM 逻辑斯蒂曲线 (图3)，估算的深水金线鱼50%性成熟体长(L50)为11.76 cm，95%性成熟体长(L95)为1.1×L50=12.94 cm。
![]()
图 3 南海北部深水金线鱼性成熟比例与体长拟合的逻辑斯蒂曲线Fig. 3 Logistic model fitted for relationship between body length and percentage of mature N. bathybius in northern South China Sea2.2 LBB模型评估结果
根据LBB评估结果(图4)显示，南海北部深水金线鱼的渐近体长Linf估算值为23.8 cm (95%置信区间为23.4~24.2 cm)，相对自然死亡率M/K为2.35 (95%置信区间为2.22~2.49)，相对生物量B/B0为0.31 (95%置信区间为0.23~0.38)，种群资源量水平B/BMSY为0.89 (95%置信区间为0.66~1.10)，Lc/Lc_opt为0.85，表明深水金线鱼处于轻度过度开发状态和生长性过度捕捞状态，估算的相对捕捞死亡率F/M为1.00 (95%置信区间为0.81~1.20)，表明捕捞死亡率和自然死亡率基本一致，开发率 (E) 可由F/(F+M)计算，结果为0.5，Lc_opt为11.1 cm，捕获的约60%的深水金线鱼个体超过了这一长度。
![]()
图 4 南海北部深水金线鱼基于LBB方法的分析结果注：左图为模型对体长数据的拟合曲线，右图为LBB方法的预测；Lc是50%个体被渔具捕获的体长；Lopt为未开发种群获得最大生物量的捕捞体长。Fig. 4 Estimated values for N. bathybius in northern South China Sea by length-based Bayesian biomass assessment modelNote: The left figure shows the fits of the model to length data, and the right figure shows the predictions of the LBB analysis; Lc. Length of 50% of the individuals captured by the gear; Lopt. Fishing body length for obtaining maximum biomass in undeveloped populations.2.3 LBSPR模型评估结果
根据LBSPR的输出结果 (图5)，南海北部深水金线鱼的相对死亡率F/M为1.71，开发率E为0.63，高于LBB计算出的开发率，SPR值为0.19，表明深水金线鱼处于过度捕捞状态。选择性和成熟度曲线显示性成熟时的长度低于首次捕获时的长度 (图5-b)。在50%和95%选择性 ($S_{L_{50}} $和$S_{L_{95}} $)下，预估的体长分别为9.9 和11.9 cm (图5-c)。假设本次的长度-频率数据处于稳定状态[26]，目前南海深水金线鱼的SPR值为0.19 (图5-c)，与其阈值水平相比 (SPR=0.4)，该种群目前所捕获样品比预期的个体小，尤其是大体长个体数量明显低于预期 (图5-d)。
![]()
图 5 南海北部深水金线鱼基于LBSPR方法的分析结果注：a. 长度-频率分布(柱状图)，预测的种群大小组成(折线图)；b. LBSPR模型拟合的性成熟和选择性曲线，其中L50=11.76 cm和L95=12.94 cm时；c. 选择性参数($S_{L_{50}}$和 $S_{L_{95}} $)、捕捞死亡率与自然死亡率的比值(F/M)和繁殖潜力比的估计值；d. 在目标SPR=0.4时，观测到的长度-频率数据与预期的大小组成。Fig. 5 Length-based spawning potential ratio for N. bathybius in northern South China SeaNote: a. Length-frequency distribution (Pillars), and predicted fished size composition (Line chart); b. Maturity and selectivity curves from the fitted LBSPR model when L50 was 11.76 cm and L95 was 12.94 cm; c. Distribution of mean selectivity parameters ($S_{L_{50}} $ and $S_{L_{95}} $), fishing mortality to natural mortality (F/M), and spawning potential ratio; d. Observed length-frequency data against an expected size composition at a target SPR of 0.4.2.4 模型不确定性分析
LBSPR和LBB 模型的估算结果对参数Linf的设置极为敏感，对参数M/K的设置较为敏感，对参数L50的设置不敏感 (表3)。从第2—第6组Linf的数值变化2 mm，即可引起2种方法估算所得的各个结果的明显变化，第1、第4、第5组参数Linf的数值相差1 mm，各个估算结果同样变化明显，可判断LBSPR和LBB模型的估算结果对参数Linf的设置极为敏感。根据第8—第13组计算结果、第8组和第13组参数L50差值为10 mm，LBB的估算结果并无明显变化，表明LBB对参数L50的设置不敏感；LBSPR的结果则出现了2种情况，根据第8—第13组计算结果，F/M的估算结果未随参数L50的变化发生明显改变，SPR的值在参数L50变化为2 mm条件下变化较小，根据第1、第9、第12组参数L50变化为3 mm条件下的结果对照来看，SPR计算结果变化较为明显，SPR对参数L50设置较为敏感，2种模型的另外4个估算结果无明显变化，LBSPR和LBB 模型的估算结果对参数L50的设置不敏感。根据第14—第19组参数M/K的数值相差0.02的设定，结果表明，2种模型的估算结果变化较为明显，根据第1、第16、第17组，参数M/K变化为0.01的情况下，LBSPR模型估算的F/M结果和LBB模型估算的B/BMSY结果变化明显，另外3个值的估算结果变化小；根据第1、第14、第19组，参数M/K的变化为0.05的情况下，LBSPR和LBB模型的估算结果产生了明显变化，这2种模型的估算结果对参数M/K设置的变化比较敏感。通过对先验参数的不确定性分析，发现LBB和LBSPR模型对参数Linf设置的变化极为敏感，对参数M/K设置的变化比较敏感，对参数L50设置的值不敏感，故在使用这2种基于长度的数据缺乏模型时对参数Linf和M/K应谨慎设置。
3. 讨论
3.1 资源状况
本研究使用LBB和LBSPR 2种基于体长的数据缺乏评估方法，对2014—2019年南海北部深水金线鱼进行评估，LBB模型评估结果显示种群资源量水平(B/BMSY)、Lc/Lc_opt和E分别为0.89、0.85和0.5，LBSPR模型评估结果显示繁殖潜力比 (SPR)为0.19，E为0.63。一般认为SPR=0.4是比较合理的目标参考点，种群在此水平下可长期维持可持续发展，而SPR=0.2常被设定为最低限制参考点，低于此值则表明种群遭受过度捕捞，种群补充量不足，难以维持种群稳定，必须采取严格的管理措施来恢复种群资源量[26,31-32]，LBB模型估算结果表明深水金线鱼正处于过度捕捞状态。Gulland[33]提出鱼类的最适开发率E为0.5，大于0.5为过度开发，2种方法的评估结果略显差异。同一种群，使用不同的数据缺乏评估方法可能产生不同的结果。例如，Xu等[34]使用LBB和LBSPR 2种基于体长的数据缺乏评估方法，对2012—2014年、2019—2021年南海细鳞圆鲹 (Decapterus macarellus) 2个种群进行资源评估，2种方法的结果均表明2019—2021年该种群已经被过度捕捞，而这2种方法在对2012—2014年南海细鳞圆鲹种群状况评估中出现了分歧：LBB模型评估结果显示为健康状态，而LBSPR模型评估结果显示为过度捕捞状态。Liao等[35]结合基于产量 (Monte Carlo Catch Maximum Sustainable Yield Model, CMSY)、长度 (LBB)和丰度指数 (Abundance-based method, AMSY) 3种数据缺乏评估方法对2005年南大西洋长鳍金枪鱼（Thunnus alalunga）进行资源评估，基于产量(CMSY)和基于长度(LBB) 的模型估算结果表明在2005年南大西洋长鳍金枪鱼种群处于健康(从过度捕捞状态中恢复健康)状态，而基于丰度指数 (AMSY) 的模型估算结果表明该种群处于过度捕捞状态。本研究中，LBB模型计算结果表明南海北部深水金线鱼处于轻度过度开发状态和生长性过度捕捞状态，而LBSPR模型计算结果表明其正处于过度捕捞状态。通过结合不同模型，能更完整地展现出种群的实际情况，为渔业管理提供更准确的科学依据。
3.2 南海北部深水金线鱼种群参数的年际变化
大多数鱼类的表型特征在经过较强的捕捞压力胁迫后会产生适应性响应，其中个体规格的变化较为明显[36]。根据陈作志等[4]对北部湾深水金线鱼种群参数的研究 (表4，本研究的种群参数在表4中记为2017年)，在持续多年的捕捞胁迫下，深水金线鱼的群体结构发生了较为明显的变化，渐近体长(Linf)和50%性成熟体长(L50)在减小后均有所回升，拐点分别为1997年和2009年，Linf从1992年242 mm减小至1997年220 mm后，1997—2009年增长到235 mm，2009—2017年略有增长，为237 mm；L50在1992—2009年呈减小趋势，从128.6 mm减小至110.5 mm，2009—2017年增长到117.6 mm。生长系数(K)在1992—2009年呈增大趋势，2009—2017年K减小至0.326，生长加快的现象有所缓解，这一变化趋势与北部湾蓝圆鲹 (D. maruadsi)[37]和南海北部带鱼(Trichiurus japonicus)[38]类似，这也与陈作志等[4]提到的环境胁迫理论相同，即过大的捕捞强度使得资源群体的密度稀疏，剩余群体的饵料丰度和易得性大幅增加，从而改善其生长环境，使得生长加快，这种变化在停止捕捞或捕捞强度下降时会出现逆转。
表 4 不同时期深水金线鱼种群参数Table 4 Population parameters of N. bathybius in different periods参数
Parameter年份 Year 1992 1997 2009 2017 渐近体长 Asymptotic length (Linf)/mm 242 220 235 237 生长系数Growth coefficient (K) 0.39 0.45 0.48 0.326 50%性成熟体长Body length at 50% sexual maturity (L50)/mm 128.6 116.5 110.5 117.6 3.3 模型参数及模型不确定性分析
Liao等[35]研究表明，在使用数据缺乏的方法时，应注意模型参数先验的设置和不同类型数据中隐含的不同趋势。对于LBSPR模型，Linf估计的不确定性反映在估计的F/M和SPR上，其中Linf输入越高，F/M越高，SPR越低，反之亦然(表3)[39]。
LBSPR模型结果的准确性依赖于以下主要假设能否满足[26]：1) 渔具选择性曲线符合逻辑斯蒂型选择性曲线；2) 鱼类在整个生命周期中体型持续增加；3) 体长频率数据具有良好的代表性，能代表所在海区内种群的体长结构特征。
LBB 是一种简单快速的估算鱼类相对种群参数的方法，与其他独立的渔业资源评估方法相比，它不需要年龄、性成熟度、补充量、生长和死亡参数、捕捞努力量等信息，只需要能代表评估对象整个世代的长度频率数据。长度频率数据是渔业生物学测定中最易获取和测量误差最小的数据，不受野外船只晃动和内脏有无的影响[40]。在长度-频率能代表资源各开发阶段的长度组成时，LBB 估算的结果能较好地反映其真实情况[41]。Zhang等[11,42]研究表明LBB模型估计参数(Linf、Lc/Lc_opt、Z/K、B/B0和B/BMSY)对样本量不敏感，电子体长频率分析 (ELEFAN) 方法的先验信息对LBB法有效。与ELEFAN相比，LBB可以估计更多的相关结果 (如B/BMSY和Lc/Lc_opt)。因此，建议使用ELEFAN和LBB方法来拟合数据贫乏的鱼类种群的长频数据，因为它们在估计管理参考点方面是互补的[11]。根据表3，通过对先验参数的不确定性分析，发现LBB和LBSPR模型对参数Linf设置的变化极为敏感，对参数M/K设置的变化较为敏感，对参数L50设置的值不敏感，故在使用这2种基于长度的数据缺乏模型时对Linf和M/K 2种参数应谨慎设置。本研究只针对M/K、L50和Linf 3个参数对LBB和LBSPR的结果影响进行分析，而渔业数据的变化多种多样，在数据采集时，并非只有一种先验参数会发生变化，在以后的研究中将进一步探讨不同误差的情况，为数据有限条件下的渔业资源评估方法优化提供依据。
-
![]()
图 2 珠江口生态系统食物网
Ⅰ~Ⅳ表示营养级等级；数字1~29表示功能组序号，序号所具体代表的功能组见表1；后图同此
Figure 2. Food web of Pearl River estuary ecosystem
Ⅰ−Ⅳ. trophic grades; 1−29. functional groups serial numbers, and the functional group specifically represented by the serial number is shown in Tab.1. The same case in the following figure.
![]()
图 3 珠江口生态系统功能组的混合营养效应
Figure 3. Mixed trophic impact of functional groups in Pearl River Estuary ecosytem
表 1 珠江口生态系统Ecopath模型的功能群及主要种类
Table 1 Functional groups and main species in Pearl River estuary
序号
No.功能组
functional group种类
species1 海洋哺乳动物 中华白海豚 Sousa chinensis 2 鲨鳐类 何氏鳐 Raja hollandi、尖头斜齿鲨 Scoliodon sorrakowah 等 3 其他大型中上层鱼类 长颌宝刀鱼 Chirocentrus nudus、四指马鲅 Eleutheronema tetradactylus、鮻鱼 Liza haematocheila、带鱼 Trichiurus haumela、短带鱼 Trichiurus brevis、二长棘鲷 Paerargyrops edita 等 4 花鲈 花鲈 Lateolabrax maculatus 5 其他小型中上层鱼类 前鳞骨鲻 Mugilophuyseni、白氏银汉鱼 Allanetta bleekeri、六指马鲅 Polynemus sextarius、丽叶鲹 Syngnathus acus Linnaeus 等 6 其他鲱形目鱼类 裘氏小沙丁鱼 Sardinella jussieu、斑鰶 Clupanodon punctatus、鳓 Ilisha elongate、康氏小公鱼 Stolephorus commersoni、赤鼻棱鳀 Thrissa kammalensis、杜氏棱鳀 Thrissa dussumieri、汉氏棱鳀 Thrissa hamiltonii、七丝鲚 Coilia grayii 等 7 凤鲚 凤鲚 Coilia mystus 8 花鰶 花鰶 Clupanodon thrissa 9 其他大型底栖鱼类 长蛇鲻 Saurida elongata、尖吻鳗 Uroconger lepturus、海鳗 Muraenesox cinereus、龙头鱼 Harpodon nehereus、中华海鲶 Arius sinensis、多鳞鱚 Sillago sihama、黄斑篮子鱼 Siganus oramin、斑点鸡笼鲳 Drepane punctata、广东鲂 Megalobrama hoffmanni、大鳞舌鳎 Cynoglossus macrolepidotus、半滑舌鳎 Cynoglossus semilaevis 等 10 石首鱼科 棘头梅童鱼 Collichthys lucidus、皮氏叫姑鱼 Johnius belengeri、截尾白姑鱼 Argyrosomus aneus、勒氏短须石首鱼 Umbrina russelli 等 11 黑鲷 黑鲷 Acanthopagrus schlegelii 12 黄鳍鲷 黄鳍鲷 Acanthopagrus latus 13 小型底层鱼类 眶棘双边鱼 Ambassis gymnocephalus、白肌银鱼 Leucosoma chinensis、鹿斑鲾 Leiognathus ruconius、短吻鲾 Leiognathus brevirostris 等 14 鲀形目鱼类 斑纹东方鲀 Fugu ocellatus、棕腹刺鲀 Gastrophysus spadiceus、黑鳃兔头鲀 Lagocephalus inermis 等 15 鰕虎鱼类 矛尾鰕虎鱼 Chaeturichthys stigmatias、钟馗鰕虎鱼 Triaenopgon barbatus、红狼牙鰕虎鱼 Odontamblyopus rubicundus、孔鰕虎鱼 Trypauchen vagina 等 16 头足类 短蛸 Octopus ocellatus、中国枪乌贼 Loligo chinensis、曼氏无针乌贼 Sepiella japonica 等 17 其他无脊椎动物 毛海参类等 18 蟹类 锯缘青蟹 Scylla serrata、红星梭子蟹 Portunus sanguinolentus、日本蟳 Charybdis japonica、锈斑蟳 Charybdis feriatus、字纹弓蟹 Varuna litterata、直额绒螯蟹 Eriocheir rectus、字纹弓蟹 Varuna litterata 等 19 其他虾类 近缘新对虾 Metapenaeus affinis、周氏新对虾 Metapenaeus joyneri、墨吉对虾 Banana prawn、脊尾白虾 Exopalaemon carinicauda、鲜明鼓虾 Alpheus distinguendus、细螯沼虾 Macrobrachium superbum、广东长臂虾 Palaemon guangdongensis、亨氏仿对虾 Parapenaeopsis hungerford 等 20 长毛对虾 长毛对虾 Penaeus penicillatus 21 斑节对虾 斑节对虾 Penaeus monodon 22 虾蛄类 黑斑口虾蛄 Oratosquilla kempi、猛虾蛄 Harpiosquilla harpax、口虾蛄 Oratosquilla oratoria、眼斑猛虾蛄 Harpiosquilla annandalei 等 23 其他大型底栖动物 体长>1 mm，包括其他棘皮动物、甲壳类、软体动物包括光滑河篮蛤 Potamocorbula laevis、河蚬 Corbicula fluminea、中国圆田螺 Cipangopaludina chinensis、菲律宾偏顶蛤 Modiolus philippinarum、小荚蛏 Siliqua minima 等 24 波纹巴非蛤 波纹巴非蛤 Paphia undulata 25 棘皮动物 海胆类、海参类、海星类、蛇尾类等 26 小型底栖动物 体长<1 mm，包括轮虫类、环节动物、多毛类、软体动物、甲壳类等 27 浮游动物 原生动物、水螅水母类、桡足类、枝角类、端足类、长尾类、短尾类、糠虾类、毛颚类、有尾类、其他浮游幼体等 28 浮游植物 硅藻、甲藻、蓝藻、绿藻类：颗粒直链藻 Melosira granulata、中肋骨条藻 Skeletonema costatum、小球藻 Chlorella sp.、黄丝藻 Tribonema sp.等 29 有机碎屑 颗粒有机碳和溶解有机碳 
下载: 导出CSV
表 2 珠江口生态系统Ecopath模型功能组估算参数
Table 2 Basic input and output parameters of groups estimated by Ecopath for Pearl River estuary ecosystem
序号
No.功能组
functional group营养级
trophic level生物量/t·km−2
B生产量/生物量/a−1
P/B消耗量/生物量/a−1
Q/B转化效率
EE1 海洋哺乳动物 4.159 0.009 0.045 14.77 0.000 2 鲨鳐类 4.152 0.012 0.820 9.500 0.000 3 其他大型中上层鱼类 3.366 0.721 0.940 5.110 0.294 4 花鲈 3.722 0.004 0.270 7.300 0.721 5 其他小型中上层鱼类 3.277 1.370 2.880 11.00 0.391 6 其他鲱形目鱼类 3.112 1.750 4.080 17.35 0.243 7 凤鲚 2.696 0.330 2.825 13.50 0.837 8 花鰶 2.931 0.120 3.650 12.10 0.769 9 其他大型底栖鱼类 3.682 0.012 2.432 27.44 0.616 10 石首鱼科鱼类 3.226 0.160 3.550 11.35 0.608 11 黑鲷 3.680 0.005 0.330 6.600 0.600 12 黄鳍鲷 3.332 0.008 0.450 8.100 0.681 13 小型底层鱼类 3.251 1.770 3.940 15.59 0.896 14 鲀形目鱼类 3.401 0.120 1.780 7.110 0.086 15 鰕虎鱼类 3.098 0.870 4.110 17.20 0.805 16 头足类 3.647 1.24 3.500 11.60 0.674 17 其他无脊椎动物 2.469 1.270 9.600 41.54 0.903 18 蟹类 2.651 0.810 5.650 26.90 0.918 19 其他虾类 2.669 1.089 6.800 41.54 0.820 20 长毛对虾 2.669 0.030 6.500 30.00 0.647 21 斑节对虾 2.669 0.005 6.500 30.00 0.452 22 虾蛄类 3.074 1.700 5.300 18.90 0.591 23 其他大型底栖动物 2.450 5.600 5.100 18.50 0.915 24 波纹巴非蛤 2.200 1.950 3.000 18.500 0.970 25 棘皮动物 2.184 1.570 10.260 41.54 0.561 26 小型底栖动物 2.250 13.76 6.570 36.90 0.789 27 浮游动物 2.000 13.50 36.000 186.0 0.437 28 浮游植物 1.000 15.00 230.000 − 0.672 29 有机碎屑 1.000 200.0 − − 0.374 注：黑体部分为模型估算值 Note: The estimated parameters by model are in bold.
下载: 导出CSV
表 3 珠江口水域生态系统的总体特征参数
Table 3 Summary statistics of net flow in Pearl River estuary ecosystem by Ecopath
特征参数
attribute parameter数值
value数值1
Value 1数值2
Value 2数值3
Value 3数值4
Value 4数值5
Value 5数值6
Value 6总消耗量/t·(km2·a)−1
total consumption3 491.934 3 492.591 3 495.201 3 497.134 3 538.434 3 540.084 4 859.454 总输出量/t·(km2·a)−1
total export1 345.636 1 345.135 1 343.186 1 341.765 1 318.511 1 317.548 473.380 总呼吸量/t·(km2·a)−1
total respiratory flow2 104.364 2 104.865 2 106.814 2 108.235 2 131.489 2 132.452 2 976.620 流向碎屑总量/t·(km2·a)−1
total flows into detritus2 150.513 2 150.012 2 148.063 2 146.746 2 128.968 2 128.204 1 415.009 系统总流量/t·(km2·a)−1
total system throughput9 092.447 9 092.603 9 093.264 9 093.880 9 117.402 9 118.288 9 724.463 总生产量/t·(km2·a)−1
total production4 139.183 4 139.208 4 139.347 4 139.472 4 149.258 4 149.616 4 360.943 总初级生产力/t·(km2·a)−1
calculated total net primary production3 450.000 3 450.000 3 450.000 3 450.000 3 450.000 3 450.000 3 450.000 总初级生产力/总呼吸量
total primary production/total respiration1.639 1.639 1.638 1.636 1.619 1.618 1.159 净系统生产量/t·(km2·a)−1
net system production1 345.636 1 345.135 1 343.186 1 341.765 1 318.511 1 317.548 473.380 总初级生产力/总生物量
total primary production/total biomass53.253 53.179 52.849 52.731 52.009 51.966 24.873 总生物量/t·(km2·a)−1
total biomass (excluding detritus)64.785 64.875 65.280 65.427 66.335 66.390 138.705 联结指数
connectance index0.370 0.370 0.370 0.370 0.370 0.370 0.370 系统杂食指数
system omnivory index0.287 0.284 0.284 0.285 0.288 0.289 0.283 注：“数值”代表当前的系统状态；“数值1”至“数值6”依次表示花鲈、黑鲷、黄鳍鲷、长毛对虾、斑节对虾和波纹巴非蛤达到最大生态容纳量后的系统特征参数 Note: The "value" represents the present status of the system; Value 1−Value 6 represent the status after a large amount of L. japonicus, A. schlegelii, A. latus, P. penicillatus, P. monodon and P. undulata biomass.
下载: 导出CSV
表 4 6种适宜放流种类在珠江口生态系统中的最大生态容纳量
Table 4 Maximum ecological carrying capacity of six species in Pearl River estuary ecosystem
序号
No.适宜放流种类
suitable release species当前生物量/t·km−2
current biomass最大容纳量/t·km−2
maximum capacity转换效率/%
conversion efficiency1 花鲈 0.004 0.094 12.28 2 黑鲷 0.005 0.500 12.46 3 黄鳍鲷 0.008 0.650 12.40 4 长毛对虾 0.030 1.580 12.33 5 斑节对虾 0.005 1.610 12.34 6 波纹巴非蛤 1.950 75.870 11.06 注：当前系统转换效率为12.23%
Note: The current system conversion efficiency is 12.23%.
下载: 导出CSV
-
[1] 袁梦, 汤勇, 徐姗楠, 等. 珠江口南沙海域秋季渔业资源群落结构特征[J]. 南方水产科学, 2017, 13(2): 18-25. doi: 10.3969/j.issn.2095-0780.2017.02.003 [2] 晏磊, 谭永光, 杨吝, 等. 珠江口水域秋季刺网的渔获组成及多样性分析[J]. 南方水产科学, 2016, 12(1): 111-119. doi: 10.3969/j.issn.2095-0780.2016.01.015 [3] 肖瑜璋, 王蓉, 张保学. 珠江口海域海洋渔业资源现状分析与建议[J]. 黑龙江科技信息, 2010(28): 233. doi: 10.3969/j.issn.1673-1328.2010.28.228 [4] 段丽杰. 基于EwE的珠江口渔业和近海生态系统模拟研究[D]. 广州: 中山大学, 2009: 9-15. [5] 程家骅, 姜亚洲. 海洋生物资源增殖放流回顾与展望[J]. 中国水产科学, 2010, 17(3): 610-617. [6] HEYMANS J J, COLL M, LINK J S, et al. Best practice in Ecopath with Ecosim food-web models for ecosystem-based management[J]. Ecol Model, 2016, 331: 173-184. doi: 10.1016/j.ecolmodel.2015.12.007
[7] De MUTSERT K, STEENBEEK J, LEWIS K, et al. Exploring effects of hypoxia on fish and fisheries in the northern Gulf of Mexico using a dynamic spatially explicit ecosystem model[J]. Ecol Model, 2016, 331: 142-150. doi: 10.1016/j.ecolmodel.2015.10.013
[8] 米玮洁, 胡菊香, 赵先富. 生态通道模型及其在水生态系统中的应用探讨[J]. 环境科学与技术, 2012, 35(7): 186-190, 196. doi: 10.3969/j.issn.1003-6504.2012.07.040 [9] 许祯行, 陈勇, 田涛, 等. 基于Ecopath模型的獐子岛人工鱼礁海域生态系统结构和功能变化[J]. 大连海洋大学学报, 2016, 31(1): 85-94. [10] 陈作志, 邱永松, 贾晓平, 等. 捕捞对北部湾海洋生态系统的影响[J]. 应用生态学报, 2008, 19(7): 1604-1610. [11] 杨超杰, 吴忠鑫, 刘鸿雁, 等. 基于Ecopath模型估算莱州湾朱旺人工鱼礁区日本蟳、脉红螺捕捞策略和刺参增殖生态容量[J]. 中国海洋大学学报(自然科学版), 2016, 46(11): 168-177. [12] 吴忠鑫, 张秀梅, 张磊, 等. 基于线性食物网模型估算荣成俚岛人工鱼礁区刺参和皱纹盘鲍的生态容纳量[J]. 中国水产科学, 2013, 20(2): 327-337. [13] 王腾, 张贺, 张虎, 等. 基于营养通道模型的海州湾中国明对虾生态容纳量[J]. 中国水产科学, 2016, 23(4): 965-975. [14] 刘玉, 隋丽杰, 段丽杰, 等. 珠江口EwE模型功能组划分研究[J]. 海洋环境科学, 2008, 27(5): 480-483. doi: 10.3969/j.issn.1007-6336.2008.05.018 [15] 农业部. 农业部关于做好“十三五”水生生物增殖放流工作的指导意见[EB/OL]. [2016-05-25]. http://www.ynagri.gov.cn/dl/news730/20160525/6269010.shtml. [16] CHRISTENSEN V, WALTERS C J. Ecopath with Ecosim: methods, capabilities and limitations[J]. Ecol Model, 2004, 172(2/3/4): 109-139.
[17] CHRISTENSEN V, WALTERS C J, PAULY D. Ecopath with ecosim: a user's guide[M]. Vancouver: University of British Columbia, 2005: 154.
[18] 王增焕, 李纯厚, 贾晓平. 应用初级生产力估算南海北部的渔业资源量[J]. 海洋水产研究, 2005, 26(3): 9-15. [19] PITCHER T J, BUCHARY E, ASUMAILA U R. Spatial simulations of Hong Kong's marine ecosystem: ecological and economic forecasting of marineprotected areas with human-made reefs[J]. Fish Centre Res Rep, 2000, 10(3): 154-158.
[20] PAULY D. On the interrelationships between natural mortality, growth parameters, and mean environmental temperature in 175 fish stocks[J]. ICES J Mar Sci, 1980, 39(2): 175-192. doi: 10.1093/icesjms/39.2.175
[21] PALOMARES M D, PAULY D. Predicting food consumption of fish populations as functions of mortality, food type, morphometrics, temperature and salinity[J]. Mar Freshw Res, 1998, 49(5): 447-453. doi: 10.1071/MF98015
[22] WANG Y, DUAN L, LI S, et al. Modeling the effect of the seasonal fishing moratorium on the Pearl River estuary using ecosystem simulation[J]. Ecol Model, 2015, 312: 406-416. doi: 10.1016/j.ecolmodel.2015.06.011
[23] 顾洪静. 福建九龙江口水域鱼类群落及其资源的研究[D]. 厦门: 集美大学, 2014: 47-62. [24] MORISSETE L, HAMMILL M O, SAVENKOFF C. The trophic role of marine mammals in the northern Gulf of St. Lawrence[J]. Mar Mammal Sci, 2006, 22(1): 74-103. doi: 10.1111/j.1748-7692.2006.00007.x
[25] 黄美珍. 台湾海峡及邻近海域主要无脊椎动物食性特征及其食物关系研究[J]. 海洋科学, 2005, 29(1): 73-80. doi: 10.3969/j.issn.1000-3096.2005.01.015 [26] 林群, 李显森, 李忠义, 等. 基于Ecopath模型的莱州湾中国对虾增殖生态容量[J]. 应用生态学报, 2013, 24(4): 1131-1140. [27] ODUM E P. The strategy of ecosystem development[J]. Science, 1969, 164(3877): 262-270. doi: 10.1126/science.164.3877.262
[28] 莫宝霖, 秦传新, 陈丕茂, 等. 基于Ecopath模型的大亚湾海域生态系统结构与功能初步分析[J]. 南方水产科学, 2017, 13(3): 9-19. doi: 10.3969/j.issn.2095-0780.2017.03.002 [29] LIN H J, SHAO K T, KUO S R, et al. A trophic model of a sandy barrier lagoon at Chiku in southwestern Taiwan[J]. Estuar Coast Shelf Sci, 1999, 48(5): 575-588. doi: 10.1006/ecss.1998.0457
[30] LIN H J, SHAO K T, HWANG J S, et al. A trophic model for Kuosheng Bay in northern Taiwan[J]. J Mar Sci Technol, 2004, 12(5): 424-432.
[31] 杨林林, 姜亚洲, 袁兴伟, 等. 象山港生态系统结构与功能的Ecopath模型评价[J]. 海洋渔业, 2015, 37(5): 399-408. doi: 10.3969/j.issn.1004-2490.2015.05.002 [32] 陈作志, 邱永松. 南海北部生态系统食物网结构、能量流动及系统特征[J]. 生态学报, 2010, 30(18): 4855-4865. [33] 陈丕茂. 广东人工鱼礁区增殖放流种类初探[J]. 南方水产, 2005, 1(1): 11-20. doi: 10.3969/j.issn.2095-0780.2005.01.003 [34] 郭晓奇. 广东海洋生物增殖放流存在的问题与对策[J]. 海洋与渔业, 2014(7): 40-41. doi: 10.3969/j.issn.1672-4046(s).2014.07.038 [35] 佚名. 农业部和广东省联合举办南海生物资源增殖放流[J]. 北京农业, 2009(21): 34. [36] 洪巧巧, 庄平, 杨刚, 等. 长江口中国花鲈食性分析[J]. 生态学报, 2012, 32(13): 4181-4190. -
期刊类型引用(25)
1. 庄平，赵峰，罗刚，张涛，石小涛，冯广朋，王思凯. 水生生物资源增殖放流的发展历程与问题思考. 水生生物学报. 2025(01): 43-54 . 
百度学术
2. 章欣仪，郑春芳，秦松，刘伟成，张川，范青松. 基于Ecopath模型的瓯江口斑鰶、刀鲚和鮻的增殖生态容量评估. 水产学报. 2025(02): 97-108 . 百度学术
3. 肖玉林，施凯，许强. 海南岛热带典型天然牡蛎礁生态系统营养结构与功能评价. 生态学报. 2025(04): 1697-1710 . 百度学术
4. 冯瑞玉，陶峰，郭禹，秦传新，孙金辉，吴一桂，王静. 基于Ecopath模型的增殖放流对南朗水域生态系统的影响. 水产学报. 2025(05): 94-108 . 百度学术
5. 王孟佳，徐开达，王好学，周永东，李鹏飞，朱凯，陈欣怡，陈璐. 浙江近海甲壳类资源增殖放流现状研究. 海洋开发与管理. 2024(01): 136-144 . 百度学术
6. 杨禧越，刘永，李纯厚，唐广隆，张达娟，游奕来，潘淑芳，郑秋实，肖雅元，吴鹏. 珠江口万山群岛黄鳍棘鲷食性随年龄的变化特征. 中国水产科学. 2024(02): 219-231 . 百度学术
7. 韩毓，张杭君. 水产生态容量及在淡水增养殖上的应用研究进展. 水产科学. 2024(04): 675-682 . 百度学术
8. 张聪，孔令宇. 海上油气田开发工程渔业资源补偿措施研究. 环境保护. 2024(15): 26-28 . 百度学术
9. 袁华荣，章守宇，陈丕茂. 海洋牧场建设效益评价研究进展与展望. 南方水产科学. 2024(05): 1-13 . 本站查看
10. 范泽宇，白雪兰，徐聚臣，黄路全，王晓宁，吕亚兵，侯杰，何绪刚. 运用Ecopath模型构建大水面增殖放流方案——以洈水水库为例. 华中农业大学学报. 2023(01): 82-91 . 百度学术
11. 方光杰，周永东，梁君，徐开达，龙翔宇，刘润泽. 岛礁型海洋牧场聚鱼增殖模式综述. 浙江海洋大学学报(自然科学版). 2023(02): 165-172 . 百度学术
12. 冯瑞玉，郭禹，李金明，孙金辉，于刚，吴一桂，秦传新. 基于EnhanceFish模型的鱼类增殖放流策略研究:以中山市南朗水域黄鳍棘鲷增殖放流为例. 渔业科学进展. 2023(05): 1-10 . 百度学术
13. 江满菊，郭禹，秦传新，潘莞倪，于刚，马振华. 黄鳍棘鲷幼鱼对不同开孔形状和直径的人工鱼礁模型的行为响应. 中国水产科学. 2023(12): 1496-1506 . 百度学术
14. 裴精花，陈清华，范金金，刘伟杰，郭照良，隋昊志. 珠江口海域游泳动物群落结构及多样性特征. 南方农业学报. 2023(12): 3727-3738 . 百度学术
15. 洪小帆，陈作志，张俊，江艳娥，龚玉艳，蔡研聪，杨玉滔. 基于Ecopath模型的七连屿礁栖性生物的生态承载力分析. 热带海洋学报. 2022(01): 15-27 . 百度学术
16. 孔业富，尹成杰，王林龙，刘杨，林黎，康斌. 基于Ecopath模型的三门湾生态系统结构与功能. 应用生态学报. 2022(03): 829-836 . 百度学术
17. 袁旸，线薇薇，张辉. 基于生态通道模型的我国渔业资源生态容量研究进展. 海洋科学. 2022(07): 105-119 . 百度学术
18. 陈璐，王好学，徐开达，吕泽砚，黄波，李哲，梁君，周永东，李鹏飞，刘连为. 浙江近岸水域岛礁性鱼类增殖放流现状分析. 浙江海洋大学学报(自然科学版). 2022(05): 459-465+472 . 百度学术
19. 马文刚，尹洪洋，孙春阳，王兆国，魏一凡，冯博轩，奉杰，许强，李秀保，王爱民. 热带典型珊瑚岛礁海洋牧场花刺参底播增殖容量及其生态效应预测. 海洋与湖沼. 2022(06): 1573-1584 . 百度学术
20. 朱克诚，刘宝锁，伞利择，刘波，张楠，郭华阳，郭梁，江世贵，张殿昌. 黃鳍棘鲷放流苗种的遗传质量评估. 广东海洋大学学报. 2021(03): 138-144 . 百度学术
21. 范泽宇，白雪兰，徐聚臣，王晓宁，吕亚兵，侯杰，何绪刚. 基于Ecopath模型的洈水水库生态系统特征及鲢、鳙生态容量分析. 中国水产科学. 2021(06): 773-784 . 百度学术
22. 王书献，张胜茂，戴阳，王永进，隋江华，朱文斌. 利用声呐数据提取磷虾捕捞深度方法研究. 南方水产科学. 2021(04): 91-97 . 本站查看
23. 许龙飞，梁志强，李君轶，郭星辰，姜海波，安苗，邵俭. 野生鲫、鳡形态性状与体质量的通径分析. 山地农业生物学报. 2021(06): 71-75 . 百度学术
24. 林坤，麦广铭，王力飞，王学锋. 2015—2018年珠江口近岸海域鱼类群落结构及其稳定性. 水产学报. 2020(11): 1841-1850 . 百度学术
25. 戴媛媛，吴会民，张韦，王健，缴建华. 基于Ecopath模型的我国海洋渔业生态系统研究概况. 海洋湖沼通报. 2020(06): 150-157 . 百度学术
其他类型引用(15)
计量
- 文章访问数: 5279
- HTML全文浏览量: 2399
- PDF下载量: 137
- 被引次数: 40
粤公网安备 44010502001741号
