水产资源 MSY 评估模拟器 返回
水产资源・MSY

水产资源 MSY 最大持续生产量模拟器 — Schaefer

对于鳕鱼、金枪鱼、鲱鱼、凤尾鱼等主要水产资源,从 Schaefer・Pella-Tomlinson・Fox 的余剩生产模型计算 MSY(最大持续生产量)・B_MSY・F_MSY。输入现有资源量和捕获压,立即显示 Kobe plot 上的象限(健全/重建中/过度捕获/危机)和预防性 TAC(捕获限额)。

参数设置
对象鱼种
显示预设的 r 和 K 参考值
环境容纳力 K
t
海域能维持的最大资源量
内部增长率 r
1/y
资源量低时的增长速度
现存资源量 B
t
当前亲鱼资源量估计值
捕获死亡率 F
1/y
捕获导致的瞬间死亡率
余剩生产模型
盈余曲线的形状假设
战略
TAC(捕获限额)计算基准
计算结果
MSY (t/y)
B_MSY (t)
F_MSY (1/y)
B/B_MSY 比
F/F_MSY 比
资源状态
海洋・鱼群动画

显示海面和鱼群、逻辑斯蒂曲线、渔船作业强度。颜色表示资源状态(绿=健全,黄=重建中,橙=过度捕获,红=危机)。

余剩生产曲线 — Schaefer / Pella-Tomlinson / Fox
Kobe Plot — B/B_MSY × F/F_MSY
理论・主要公式

$$MSY = \frac{rK}{4},\quad B_{MSY} = \frac{K}{2},\quad F_{MSY} = \frac{r}{2}$$

Schaefer 逻辑斯蒂余剩生产模型。r:内部增长率,K:环境容纳力,MSY:最大持续生产量。当 B = K/2 时余剩生产最大,以死亡率 r/2 的捕获能持续漁获 MSY。

$$\frac{dB}{dt} = r\,B\left(1-\frac{B}{K}\right) - F\,B$$

资源量 B 的时间变化。右边第 1 项为余剩生产,第 2 项为捕获。dB/dt = 0 的定常解为平衡资源量,F 超过 F_MSY 时平衡 B 会低于 B_MSY。

水产资源 MSY 评估 — Schaefer・Pella-Tomlinson 模型

🙋
「MSY」经常在渔业新闻中听到,到底是什么数字?是「再捕多了就完蛋」的界线吗?
🎓
问得好。MSY(最大持续生产量,Maximum Sustainable Yield)是「每年捕获这么多资源也不会减少的最大量」。用公式表示是 MSY = rK/4,r 是资源的内部增长率,K 是海域的环境容纳力。用默认的大西洋鳕鱼,r=0.3、K=100 万吨,所以 MSY=7.5 万吨/年。低于这个值资源向回复方向,超过这个值长期会枯竭,这就是分界点。
🙋
那为什么不能每年捕到 MSY 的极限?为什么要「预防性地减到 0.8 MSY」?
🎓
有两个原因。第一个是 r 和 K 的估计值本身误差很大,通常 ±30~50% 很正常,真正的 MSY 可能比想的低。每年都按 MSY 捕获,资源在某个时期就会崩。第二个是环保变动,比如厄尔尼诺现象,某年新生一代鱼的数量会急剧下降。1972 年秘鲁凤尾鱼在一年内从几百万吨跌到几万吨。所以 FAO 规定「MSY 的 70~90% 为上限」,这就是预防性方法的含义。
🙋
右边的 Kobe plot,绿・黄・橙・红 4 个象限出现了。默认显示 B/B_MSY=1.20、F/F_MSY=1.00 为「绿(健全)」,只增加 F 会怎样?
🎓
试试把 F 滑块调到 0.20。F/F_MSY=1.33,象限变成「橙(过度捕获)」。这表示「现在资源量还够,但这样持续捕下去,几年后 B 就会掉到 B_MSY 以下」。再把 F 调到 0.30,同时把 B 滑块降到 30 万吨,就变成「红(危机)」。纽芬兰鳕鱼 1992 年陷入的就是这种状态,导致 30 万人失业,到现在 30 年了还没完全恢复。
🙋
为什么有 Schaefer・Pella-Tomlinson・Fox 3 个余剩生产模型?同样的 MSY 评估,必须分开用吗?
🎓
因为根据鱼种的生活史,最优形状不同。Schaefer 是 B_MSY = K/2 对称的钟形。这对寿命长、繁殖稳定的鳕鱼这样的鱼种很好。Pella-Tomlinson 引入 m 参数可以非对称化,适合秋刀鱼、沙丁鱼这样短命、高繁殖力的鱼种。Fox 是 Gompertz 型,B_MSY ≈ 0.37K,早期达到 MSY,用于保守评估最大持续生产或幼鱼承受高捕获压的资源。工具上面的选项换一下,你能看到余剩生产曲线的形状改变。
🙋
这个余剩生产模型看起来很简单,实际工作中直接这样用吗?没有更高级的方法吗?
🎓
其实 ICCAT 的黑鲔评估、WCPFC 的大眼金枪鱼评估中,主流是 Statistical Catch-at-Age(年龄别资源量分析)和 Stock Synthesis 这样的年龄结构模型。这些内含 Beverton-Holt 或 Ricker 的股票补充关系,能再现加入的亲鱼量依赖性。但数据少的资源,现在国际上还普遍用 Schaefer 系的生物量动态模型。把本工具看作「最初的第一认识」的入门工具。MSC(海洋管理委员会)的可持续海鲜认证内部也在看 B/B_MSY 和 F/F_MSY。

常见问题

Schaefer 的逻辑斯蒂余剩生产模型 dB/dt = rB(1−B/K) 中,余剩生产在最大点为 MSY = rK/4,此时资源量为 B_MSY = K/2、捕获死亡率为 F_MSY = r/2。例如内部增长率 r=0.30、环境容纳力 K=1,000,000 t 的鳕鱼资源,则 MSY = 75,000 t/年、B_MSY = 500,000 t、F_MSY = 0.15/年。本工具从鱼种别预设值实时计算。
Kobe plot 以 B/B_MSY 为横轴、F/F_MSY 为纵轴,将资源状态分 4 区。绿(B≥B_MSY 且 F≤F_MSY)为健全,黄(B<B_MSY 且 F≤F_MSY)为重建中,橙(B≥B_MSY 且 F>F_MSY)为过度捕获,红(B<B_MSY 且 F>F_MSY)为危机。ICCAT・WCPFC 等国际管理机构的标准可视化手法,在 TAC(捕获限额)谈判中广泛应用。
Schaefer 是最简单的式,假设 B_MSY = K/2 的对称型。Pella-Tomlinson 引入形状参数 m 可实现非对称化,适用于短命、繁殖力强的鱼种(如秋刀鱼、沙丁鱼),MSY 在 K/2 的低侧。Fox 是 Gompertz 型,B_MSY ≈ K/e ≈ 0.37K,早期达 MSY,用于保守评估最大持续生产或幼鱼承受高捕获压的资源。
预防性办法(precautionary approach)是考虑估计误差和气候变动导致的加入量变化,将 TAC 定在 MSY 的基础上加安全裕度的方法。1995 年《联合国公海捕鱼协定》国际通用,FAO 指南通常为 MSY 的 70~90% 为上限。本工具可选 0.8 MSY(预防性)和 0.5 MSY(资源重建中)两种战略。1992 年纽芬兰鳕鱼崩溃是过信 MSY 的典型例子。

现实应用

国际渔业管理机构(ICCAT・WCPFC・NPFC)的 TAC 计算:管理大西洋黑鲔的 ICCAT、太平洋黑鲔的 WCPFC、北太平洋秋刀鱼的 NPFC,每年由各成员国科学委员会估算 B/B_MSY 和 F/F_MSY,制作 Kobe plot,据此谈判次年的 TAC(捕获限额)。2010 年代太平洋黑鲔在「红」象限,通过紧急削减捕获和保护未成熟鱼,现在已回复到接近「绿」象限。

历史崩溃事件的事后分析:1992 年纽芬兰大西洋鳕鱼崩溃(30 万人失业)、1972 年秘鲁凤尾鱼崩溃(200 万 t → 2 万 t)、1947 年加州沙丁鱼崩溃、纽芬兰狭鳎崩溃等 20 世纪后半的主要资源崩溃事件,用 Schaefer / Fox 系模型事后重现时,都证实「长期处于橙~红象限」。这在现代渔业教科书中被反复引用为 MSY 概念重要性的实证例子。

MSC(海洋管理委员会)认证审计:国际可持续海鲜认证 MSC 的供应链认证中,原料渔业的 B/B_MSY ≥ 0.8 和 F/F_MSY ≤ 1.0 为「可持续渔业」判定标准。超市里看到的蓝色标签鱼后面,正在运行与本工具相同的 Schaefer 系模型进行评估。

气候变动下的资源管理情景分析:海水温上升导致 K(环境容纳力)长期变化时,需要「动态 MSY」而不是固定 MSY。用本工具把 K 降到 70% 后重新计算,能确认 MSY 也同比例下降。挪威海鳕鱼、白令海狭鳎等冷水鱼种的北上位移,目前正进行这种动态 MSY 评估讨论。

常见误解和注意事项

最大陷阱是把「MSY=目标捕获量」理解错了。MSY 是理论上的上限,每年都捕到它的话,一旦估计误差(r 和 K 的不确定性)和环保变动同时朝坏方向摇摆,资源就会崩。FAO・联合国公海捕鱼协定定的预防性办法是,将 TAC 定在 MSY 的 70~90%,更进一步配合 HCR(捕获控制规则),当资源量降到 B_MSY 以下时自动削减捕获。本工具默认选「预防性 0.8 MSY」正是遵循这个国际协议。要把 MSY 理解为「绝对不能超过的天花板」而非「目标」。

其次是过度信任余剩生产模型的「决定论・即时平衡」。现实资源的加入量年变动(recruitment variability,CV 30~100%)才是支配性的,特别是短命的凤尾鱼、玉筋鱼这样的鱼种,「应该的 MSY」的一半都捕不到是常事。需要引入 Beverton-Holt 或 Ricker 的股票补充函数表示亲鱼量与加入量的关系、用 Statistical Catch-at-Age 模型求解年龄结构这样的高度化。本工具的结果应理解为「最初的感觉值」,实际评估要用 Stock Synthesis 或 SS3 这样的专用软件。

最后是只看「F_MSY 以下就安全」,忽视 B(资源量)。长期处于 Kobe plot 的橙象限(B ≥ B_MSY、F > F_MSY),最终 B 也会掉到 B_MSY 以下,转变为黄→红。反之黄象限(B < B_MSY、F ≤ F_MSY)是「重建中」,只要压低捕获,几年内就能回复到绿。不能只看 F 或只看 B,必须在 Kobe plot 上同时确认两者,追踪移动的轨迹(trajectory),这是现代资源管理的基本。

使用指南

  1. 输入环境容纳力 K(吨)。例如日本沙丁鱼资源 K=9,000,000 吨。
  2. 设置内部增长率 r(年^-1)。例如狭鳎约 r=0.35 为标准值。
  3. 输入现存资源量 B(吨)和捕获死亡率 F(年^-1),执行计算后得出 MSY・B_MSY・F_MSY。
  4. 从输出的 B/B_MSY 比和 F/F_MSY 比,在 Kobe plot 上判定资源状态(健全性/过度捕获/回复阶段)。

具体计算例

扇贝资源情况:K=450,000 吨,r=0.42 年^-1,现存资源 B=280,000 吨,捕获死亡率 F=0.18 年^-1,计算得 MSY=47,250 吨/年、B_MSY=225,000 吨、F_MSY=0.21 年^-1。B/B_MSY=1.24、F/F_MSY=0.86,资源量超过 MSY 水准,捕获压力适当,判定为「健全度高」状态。

实际工作注意点