個体群の増殖モデル

指数増加とロジスティック増殖、環境収容力の意味を数理モデルで理解する。

この単元の目標
  • 理想条件下の指数増加と資源制限下のロジスティック増殖の違いを説明できる
  • 環境収容力 K が個体数変動にどう効くかを式で読み取れる
  • 増加量が最大になる個体数 K/2 の意味と資源管理への応用を説明できる

2つの増殖モデル

個体群とは、同じ種の個体が同じ場所・同じ時間に生活している集団のことです。食物や空間が無限にある理想条件下では、増加率は個体数に比例するので、個体数は時間とともに J 字型の指数増加を示します。しかし現実の環境では資源が有限で、個体数が大きくなるほど種内競争によって1個体あたりの増加率は下がります。この密度効果を組み込んだのが S 字型のロジスティック増殖モデルで、個体数はやがて環境収容力 K に漸近します。

dNdt=rN(指数増加)\frac{dN}{dt} = rN \quad (\text{指数増加})
dNdt=rN(1−NK)(ロジスティック増殖)\frac{dN}{dt} = rN\left(1 - \frac{N}{K}\right) \quad (\text{ロジスティック増殖})
項目指数増加ロジスティック増殖
前提資源は無限・競争なし資源は有限・種内競争あり
1個体あたりの増加率r で一定r(1 − N/K):N が増えるほど低下する
曲線の形J 字型S 字型
長期的な帰結資源の枯渇で破綻する(現実には続きにくい)一定の r・K の連続時間モデルでは K に単調に近づく

なぜ S 字になるのか

修正項 (1−NK)\left(1 - \frac{N}{K}\right) は、密度による増殖の抑制を単純化したものです。残った資源を直接測定した割合ではありません。r > 0、K > 0 とすると、1個体あたりの増加率は N とともに低下しますが、個体群全体の増加速度 dN/dt は N = K/2 までは大きくなり、それ以降は小さくなります。0 < N < K なら増加し、N > K なら減少して K に近づきます。K は一定環境での平衡個体数で、現実には季節や資源量に応じて変わります。

🔍 増加率はどこで最大になるか——最大持続生産量への道

  1. Step 1: 個体数が少ないとき
    N≪K  ⇒  (1−NK)≈1  ⇒  dNdt≈rNN \ll K \;\Rightarrow\; \left(1 - \frac{N}{K}\right) \approx 1 \;\Rightarrow\; \frac{dN}{dt} \approx rN

    あいた資源が豊富なので、ほぼ指数増加と同じように増えられます。

  2. Step 2: 個体数が K/2 のとき
    dNdt=r⋅K2⋅(1−12)=rK4\frac{dN}{dt} = r \cdot \frac{K}{2} \cdot \left(1 - \frac{1}{2}\right) = \frac{rK}{4}

    個体数が十分あって、なお資源の余裕もある——増加量 dN/dt が最大になるのはここです。

  3. Step 3: 個体数が K に達したとき
    N→K  ⇒  (1−NK)→0  ⇒  dNdt→0N \to K \;\Rightarrow\; \left(1 - \frac{N}{K}\right) \to 0 \;\Rightarrow\; \frac{dN}{dt} \to 0

    資源が使い尽くされ、出生数と死亡数がつり合います。

  • 最大持続生産量: この理想モデルでは N=K2N = \frac{K}{2} で自然増加速度が最大の rK4\frac{rK}{4} になります。現実の漁業では環境変動・推定誤差・年齢構成などがあり、この値の漁獲を続ければ安全とは限りません
  • 密度効果: 密度が上がると餌不足・感染症の流行・縄張り争いなどで出生率が下がり死亡率が上がる
  • 現実の変動: 資源回復の時間遅れや気候・天敵といった外因のために、個体数は K の周りで振動することもあります
例題

r = 0.1(1年あたり)、K = 1000 個体の個体群で、現在の個体数が N = 500 のとき、増加率 dN/dt はいくらか。

解答・解説

dNdt=0.1×500×(1−5001000)=25\frac{dN}{dt} = 0.1 \times 500 \times \left(1 - \frac{500}{1000}\right) = 25 個体/年。ちょうど N = K/2 なので、増加率が最大になっている局面です。

「なぜ K/2 で最も増えるのか」を直感的に言うと、N が小さすぎると増やす側の個体が少なく、N が大きすぎると資源が足りません。両者のバランス点が K/2 です。

🔍 平方完成で最大値を確かめる

  1. rN(1−N/K)=−rK(N−K/2)2+rK4rN(1-N/K)=-\frac{r}{K}(N-K/2)^2+\frac{rK}{4}

    r と K が正なら、二乗の項の係数は負です。

  2. N=K/2⇒(dN/dt)max⁡=rK/4N=K/2\quad\Rightarrow\quad (dN/dt)_{\max}=rK/4

    二乗の項が0のとき最大になります。dN/dt の単位は個体/年など、r の単位は1/年です。

同じ増加速度でも、密度効果は違う

r = 0.2/年、K = 800個体で、N = 200 と N = 600 の増加速度を比較せよ。

解答・解説

N = 200 では0.2×200×0.75=30個体/年、N = 600 では0.2×600×0.25=30個体/年で等しくなります。一方、1個体あたりではそれぞれ0.15/年と0.05/年です。全体の増加速度と1個体あたりの増加率を混同しないことが大切です。

練習問題
  1. 指数増加モデルが自然界で長期間成り立たない理由を、資源との関係から説明せよ。

    ヒント

    個体数が増えたとき、1個体あたりの資源量はどうなるか考える。

    解答・解説

    個体数が増えると1個体あたりの資源(食物・空間)が減り、種内競争で出生率が下がったり死亡率が上がったりするためです。増加率は個体数に比例し続けられず、必ず減速します。

  2. r = 0.2/年、K = 800個体のとき、N = 200 における増加速度 dN/dt を求めよ。

    ヒント

    dNdt=rN(1−NK)\frac{dN}{dt} = rN\left(1 - \frac{N}{K}\right) に代入する。

    解答・解説

    0.2×200×(1−200800)=40×0.75=300.2 \times 200 \times \left(1 - \frac{200}{800}\right) = 40 \times 0.75 = 30 個体/年。これは現在の瞬間的な増加速度であり、1年後に必ず230個体になるという意味ではありません。

  3. 同じモデル(r = 0.2/年、K = 800個体)で N = 1000 のとき、増加速度とその後の変化を説明せよ。

    ヒント

    1 − N/K の符号を見る。

    解答・解説

    0.2×1000×(1−1000/800)=−50個体/年。個体数は減少します。一定の r・K の連続時間モデルでは K を下回って振動するのではなく、上側から K に近づきます。

  4. ロジスティック増殖をする魚群(r = 0.4、K = 10000)を漁獲するとき、持続的に最大の漁獲を得るには個体数をいくつに保てばよいか。そのときの自然増加量も求めよ。

    ヒント

    増加量が最大になる個体数は K/2。

    解答・解説

    個体数を N=K2=5000N = \frac{K}{2} = 5000 に保つ。このときの自然増加量は rK4=0.4×100004=1000\frac{rK}{4} = \frac{0.4 \times 10000}{4} = 1000 個体/年で、これを上回って漁獲すると翌年以降の資源量が減ります。

確認クイズ

問1. ロジスティック増殖で、単位時間あたりの増加量 dN/dt が最大になるのはどのとき?

問2. 環境収容力 K の説明として正しいものはどれ?

問3. 個体数が一時的に K を超えたとき、どうなるか。

📝 メモこの端末のみに保存
メモ欄はページ読み込み後に使えるようになります。