Class of ecological models
理論生態学 および 非線形動態学 において 、 消費者資源モデル (CRM)は、消費者種の群集 が共通の資源プールをめぐって競争する 生態学的モデル の一種である 。種同士が直接相互作用するのではなく、すべての種間相互作用は資源動態を介して媒介される。消費者資源モデルは、 ニッチ構築 、 共存 、 生物多様性 に関する理論の定量的発展において、基本的なツールとして機能してきた。これらのモデルは、単一の 栄養段階 の定量的記述として解釈することができる。 [1] [2]
一般的な消費者資源モデルは、存在量が である M 個の資源と 、 個体数が である S 個 の消費者種から構成されます。一般的な消費者資源モデルは、連立常微分方程式によって記述されます。
ここで 、 は資源の存在量のみに依存し、 は種の 1 人当たりの成長率 、 は資源の成長率です 。CRM の重要な特徴は、種の成長率と個体数が資源を介して媒介され、明示的な種間相互作用がないことです。資源相互作用を通じて、 種間相互作用が
出現します。
R
1
,
…
,
R
M
{\displaystyle R_{1},\dots ,R_{M}}
N
1
,
…
,
N
S
{\displaystyle N_{1},\dots ,N_{S}}
d
N
i
d
t
=
N
i
g
i
(
R
1
,
…
,
R
M
)
,
i
=
1
,
…
,
S
,
d
R
α
d
t
=
f
α
(
R
1
,
…
,
R
M
,
N
1
,
…
,
N
S
)
,
α
=
1
,
…
,
M
{\displaystyle {\begin{aligned}{\frac {\mathrm {d} N_{i}}{\mathrm {d} t}}&=N_{i}g_{i}(R_{1},\dots ,R_{M}),&&\qquad i=1,\dots ,S,\\{\frac {\mathrm {d} R_{\alpha }}{\mathrm {d} t}}&=f_{\alpha }(R_{1},\dots ,R_{M},N_{1},\dots ,N_{S}),&&\qquad \alpha =1,\dots ,M\end{aligned}}}
g
i
{\displaystyle g_{i}}
i
{\displaystyle i}
f
α
{\displaystyle f_{\alpha }}
α
{\displaystyle \alpha }
消費者資源モデルは、ロバート・H・マッカーサー [3] と リチャード・レビンス [ 4] によって最初に導入され 、生態学的原理を定式化し、微生物生態系を含む実験をモデル化することに成功しました。 [5] [6]
モデル
ニッチモデル
ニッチモデルはCRMの注目すべきクラスであり、連立常微分方程式のシステムによって記述される。 [7] [8]
d
N
i
d
t
=
N
i
g
i
(
R
)
,
i
=
1
,
…
,
S
,
d
R
α
d
t
=
h
α
(
R
)
+
∑
i
=
1
S
N
i
q
i
α
(
R
)
,
α
=
1
,
…
,
M
,
{\displaystyle {\begin{aligned}{\frac {\mathrm {d} N_{i}}{\mathrm {d} t}}&=N_{i}g_{i}(\mathbf {R} ),&&\qquad i=1,\dots ,S,\\{\frac {\mathrm {d} R_{\alpha }}{\mathrm {d} t}}&=h_{\alpha }(\mathbf {R} )+\sum _{i=1}^{S}N_{i}q_{i\alpha }(\mathbf {R} ),&&\qquad \alpha =1,\dots ,M,\end{aligned}}}
ここで 、は資源存在量のベクトル略語、 は種の一人当たり成長率 、は 消費がない場合 の種の成長率、は種が消費を通じて 資源の豊富さを枯渇させる 単位種個体群当たりの速度です 。この種類のCRMでは、消費種の資源への影響は明示的に調整されていませんが、暗黙的な相互作用が存在します。
R
≡
(
R
1
,
…
,
R
M
)
{\displaystyle \mathbf {R} \equiv (R_{1},\dots ,R_{M})}
g
i
{\displaystyle g_{i}}
i
{\displaystyle i}
h
α
{\displaystyle h_{\alpha }}
α
{\displaystyle \alpha }
−
q
i
α
{\displaystyle -q_{i\alpha }}
i
{\displaystyle i}
α
{\displaystyle \alpha }
マッカーサー消費者資源モデル(MCRM)
ロバート・H・マッカーサー にちなんで名付けられたマッカーサー消費者資源モデル(MCRM)は、 ニッチ理論と共存理論の発展における基礎的なCRMです。 [9] [10] MCRMは、次の一連の連立常微分方程式で与えられます。 [11] [12] [8] ここで、 は 資源に対する 種の相対的な選好であり、 消費者種の消費によって 資源が枯渇する相対的な量です 。 は消費がない場合(つまり、 がゼロの場合) の資源の定常収容力です。 と は それぞれ種と資源のダイナミクスの時間スケールです。 は資源の質です 。 は種の自然死亡率です。 のとき、各資源が独立した ロジスティック成長 を示す ため、このモデルは自己補充的な資源ダイナミクスを持っていると言われています 。正のパラメータと初期条件が与えられた場合、このモデルは唯一の 侵入不可能な 定常状態 (すなわち、絶滅に追い込まれた種または枯渇した資源の再導入が、再導入された種または資源の絶滅につながる定常状態)に近づく。 [7] MCRMの定常状態は 競争排除原理を満たす。すなわち、共存する種の数は枯渇していない資源の数以下である。言い換えれば、同時に占有可能な 生態学的ニッチ の数は 枯渇していない資源の数に等しい。
d
N
i
d
t
=
τ
i
−
1
N
i
(
∑
α
=
1
M
w
α
c
i
α
R
α
−
m
i
)
,
i
=
1
,
…
,
S
,
d
R
α
d
t
=
r
α
K
α
(
K
α
−
R
α
)
R
α
−
∑
i
=
1
S
N
i
c
i
α
R
α
,
α
=
1
,
…
,
M
,
{\displaystyle {\begin{aligned}{\frac {\mathrm {d} N_{i}}{\mathrm {d} t}}&=\tau _{i}^{-1}N_{i}\left(\sum _{\alpha =1}^{M}w_{\alpha }c_{i\alpha }R_{\alpha }-m_{i}\right),&&\qquad i=1,\dots ,S,\\{\frac {\mathrm {d} R_{\alpha }}{\mathrm {d} t}}&={\frac {r_{\alpha }}{K_{\alpha }}}\left(K_{\alpha }-R_{\alpha }\right)R_{\alpha }-\sum _{i=1}^{S}N_{i}c_{i\alpha }R_{\alpha },&&\qquad \alpha =1,\dots ,M,\end{aligned}}}
c
i
α
{\displaystyle c_{i\alpha }}
i
{\displaystyle i}
α
{\displaystyle \alpha }
α
{\displaystyle \alpha }
i
{\displaystyle i}
K
α
{\displaystyle K_{\alpha }}
α
{\displaystyle \alpha }
c
i
α
{\displaystyle c_{i\alpha }}
τ
i
{\displaystyle \tau _{i}}
r
α
−
1
{\displaystyle r_{\alpha }^{-1}}
w
α
{\displaystyle w_{\alpha }}
α
{\displaystyle \alpha }
m
i
{\displaystyle m_{i}}
i
{\displaystyle i}
c
i
α
=
0
{\displaystyle c_{i\alpha }=0}
外部供給リソースモデル
外部供給型資源モデルはMCRMに類似していますが、資源は自己補充されるのではなく、外部から一定の割合で供給されます。このモデルは線形資源動態モデルとも呼ばれます。これは、以下の連立常微分方程式で記述されます。 [11] [8] ここで、MCRMと共有されるすべてのパラメータは同じであり、は生態系への 資源供給速度です 。eCRMでは、消費がない場合、は 時間スケール で指数関数的 に減少します 。このモデルはケモスタットモデルとも呼ばれます。
d
N
i
d
t
=
τ
i
−
1
N
i
(
∑
α
=
1
M
w
α
c
i
α
R
α
−
m
i
)
,
i
=
1
,
…
,
S
,
d
R
α
d
t
=
r
α
(
κ
α
−
R
α
)
−
∑
i
=
1
S
N
i
c
i
α
R
α
,
α
=
1
,
…
,
M
,
{\displaystyle {\begin{aligned}{\frac {\mathrm {d} N_{i}}{\mathrm {d} t}}&=\tau _{i}^{-1}N_{i}\left(\sum _{\alpha =1}^{M}w_{\alpha }c_{i\alpha }R_{\alpha }-m_{i}\right),&&\qquad i=1,\dots ,S,\\{\frac {\mathrm {d} R_{\alpha }}{\mathrm {d} t}}&=r_{\alpha }(\kappa _{\alpha }-R_{\alpha })-\sum _{i=1}^{S}N_{i}c_{i\alpha }R_{\alpha },&&\qquad \alpha =1,\dots ,M,\end{aligned}}}
κ
α
{\displaystyle \kappa _{\alpha }}
α
{\displaystyle \alpha }
R
α
{\displaystyle R_{\alpha }}
κ
α
{\displaystyle \kappa _{\alpha }}
r
α
−
1
{\displaystyle r_{\alpha }^{-1}}
ティルマン消費者資源モデル(TCRM)
G・デイヴィッド・ティルマン にちなんで名付けられたティルマン消費者資源モデル(TCRM)は、 外部供給資源モデルに類似していますが、種が資源を枯渇させる速度が、もはやその資源の現在の豊富さに比例しないという点が異なります。TCRMはティルマンの R*則 の基礎モデルです。これは、以下の連立常微分方程式で記述されます。 [8] ここで、すべてのパラメータはMCRMと共有されます。TCRMでは、資源の豊富さは非物理的に負になることがあります。 [7]
d
N
i
d
t
=
τ
i
−
1
N
i
(
∑
α
=
1
M
w
α
c
i
α
R
α
−
m
i
)
,
i
=
1
,
…
,
S
,
d
R
α
d
t
=
r
α
(
K
α
−
R
α
)
−
∑
i
=
1
S
N
i
c
i
α
,
α
=
1
,
…
,
M
,
{\displaystyle {\begin{aligned}{\frac {\mathrm {d} N_{i}}{\mathrm {d} t}}&=\tau _{i}^{-1}N_{i}\left(\sum _{\alpha =1}^{M}w_{\alpha }c_{i\alpha }R_{\alpha }-m_{i}\right),&&\qquad i=1,\dots ,S,\\{\frac {\mathrm {d} R_{\alpha }}{\mathrm {d} t}}&=r_{\alpha }(K_{\alpha }-R_{\alpha })-\sum _{i=1}^{S}N_{i}c_{i\alpha },&&\qquad \alpha =1,\dots ,M,\end{aligned}}}
微生物消費者資源モデル(MiCRM)
微生物消費者資源モデルは、外部から供給される資源を持つ微生物生態系を記述する。これらの資源の消費によって代謝副産物が生成され、潜在的な 相互摂食 が生じる可能性がある。これは、以下の連立常微分方程式によって記述される。 ここで、MCRMと共有されるすべてのパラメータは同様の解釈を持つ。 は、資源の消費によって副産物が 資源に変換される割合であり 、は、 代謝副産物として環境に放出される資源の量を決定する 資源の「漏出率」である。 [13] [8]
d
N
i
d
t
=
τ
i
−
1
N
i
(
∑
α
=
1
M
(
1
−
l
α
)
w
α
c
i
α
R
α
−
m
i
)
,
i
=
1
,
…
,
S
,
d
R
α
d
t
=
κ
α
−
r
R
α
−
∑
i
=
1
S
N
i
c
i
α
R
α
+
∑
i
=
1
S
∑
β
=
1
M
N
i
D
α
β
l
β
w
β
w
α
c
i
β
R
β
,
α
=
1
,
…
,
M
,
{\displaystyle {\begin{aligned}{\frac {\mathrm {d} N_{i}}{\mathrm {d} t}}&=\tau _{i}^{-1}N_{i}\left(\sum _{\alpha =1}^{M}(1-l_{\alpha })w_{\alpha }c_{i\alpha }R_{\alpha }-m_{i}\right),&&\qquad i=1,\dots ,S,\\{\frac {\mathrm {d} R_{\alpha }}{\mathrm {d} t}}&=\kappa _{\alpha }-rR_{\alpha }-\sum _{i=1}^{S}N_{i}c_{i\alpha }R_{\alpha }+\sum _{i=1}^{S}\sum _{\beta =1}^{M}N_{i}D_{\alpha \beta }l_{\beta }{\frac {w_{\beta }}{w_{\alpha }}}c_{i\beta }R_{\beta },&&\qquad \alpha =1,\dots ,M,\end{aligned}}}
D
α
β
{\displaystyle D_{\alpha \beta }}
β
{\displaystyle \beta }
α
{\displaystyle \alpha }
l
α
{\displaystyle l_{\alpha }}
α
{\displaystyle \alpha }
対称的な相互作用と最適化
マッカーサーの最小化原則
マッカーサー消費者資源モデル(MCRM)において、 マッカーサーは モデルの侵入不可能な定常状態(個体数がゼロの種が再導入されても侵入に失敗することになり、生態系が前述の定常状態に戻るような定常状態)を特定するための最適化原理を導入した。最適化原理を導出するには、資源動態が十分に速く(つまり、 )、種の動態に同調して常に定常状態(つまり、 )にあると仮定し、 は の関数として表される 。この仮定のもと、種の動態は と表すことができる。
ここで は を満たす資源存在量の合計を表す 。上記の式は と書くことができ 、ここで、
r
α
≫
1
{\displaystyle r_{\alpha }\gg 1}
d
R
α
/
d
t
=
0
{\displaystyle {\mathrm {d} }R_{\alpha }/{\mathrm {d} }t=0}
R
α
{\displaystyle R_{\alpha }}
N
i
{\displaystyle N_{i}}
d
N
i
d
t
=
τ
i
−
1
N
i
[
∑
α
∈
M
∗
r
α
−
1
K
α
w
α
c
i
α
(
r
α
−
∑
j
=
1
S
N
j
c
j
α
)
−
m
i
]
,
{\displaystyle {\frac {\mathrm {d} N_{i}}{\mathrm {d} t}}=\tau _{i}^{-1}N_{i}\left[\sum _{\alpha \in M^{\ast }}r_{\alpha }^{-1}K_{\alpha }w_{\alpha }c_{i\alpha }\left(r_{\alpha }-\sum _{j=1}^{S}N_{j}c_{j\alpha }\right)-m_{i}\right],}
∑
α
∈
M
∗
{\displaystyle \sum _{\alpha \in M^{\ast }}}
R
α
=
r
α
−
∑
j
=
1
S
N
j
c
j
α
≥
0
{\displaystyle R_{\alpha }=r_{\alpha }-\sum _{j=1}^{S}N_{j}c_{j\alpha }\geq 0}
d
N
i
/
d
t
=
−
τ
i
−
1
N
i
∂
Q
/
∂
N
i
{\displaystyle \mathrm {d} N_{i}/\mathrm {d} t=-\tau _{i}^{-1}N_{i}\,\partial Q/\partial N_{i}}
Q
(
{
N
i
}
)
=
1
2
∑
α
∈
M
∗
r
α
−
1
K
α
w
α
(
r
α
−
∑
j
=
1
S
c
j
α
N
j
)
2
+
∑
i
=
1
S
m
i
N
i
.
{\displaystyle Q(\{N_{i}\})={\frac {1}{2}}\sum _{\alpha \in M^{\ast }}r_{\alpha }^{-1}K_{\alpha }w_{\alpha }\left(r_{\alpha }-\sum _{j=1}^{S}c_{j\alpha }N_{j}\right)^{2}+\sum _{i=1}^{S}m_{i}N_{i}.}
侵入不可能な定常状態において、 すべての生存種 と すべての絶滅種が対象となります 。 [14] [15]
∂
Q
/
∂
N
i
=
0
{\displaystyle \partial Q/\partial N_{i}=0}
i
{\displaystyle i}
∂
Q
/
∂
N
i
>
0
{\displaystyle \partial Q/\partial N_{i}>0}
i
{\displaystyle i}
最小環境擾乱原理(MEPP)
マッカーサーの最小化原理は、より一般的な最小環境擾乱原理 (MEPP) に拡張され、特定のニッチ CRM モデルを制約付き最適化問題にマッピングします。 資源を消費することによって種に与えられる個体群増加が、方程式 を通じて種の消費が資源の豊富さに与える影響と関連している場合、 種と資源の相互作用は 対称 であると言われています。 上記の方程式で 、 と は 資源の豊富さの任意の関数です。 この対称条件が満たされている場合、次のような関数が存在することが示されます 。 [7] この関数 を決定した後 、定常状態の侵入不可能な資源の豊富さと種の個体数が、 制約付き最適化問題 の解です。 種の個体数は、 2 行目の制約の ラグランジュ乗数 です。これは、 をラグランジュ乗数とする KKT 条件 を見るとわかります。 行 1、3、4 は実現可能性と侵入不可能性のステートメントです。 の場合 、 は ゼロでなければならず、そうでなければシステムは定常状態にありません。また 、 の場合、 は 非正でなければならず、そうでなければ種が 侵入できるようになります。行 2 は、優れた CRM 内のリソースの定常条件と定常状態条件です。関数 は、 リソースの存在量 の状態空間で がゼロになる点、 を最小値と定義することにより、距離として解釈できます 。 上記の KKT 条件につながる 双対問題 の ラグランジアンは 次のとおりです。この図では、 (つまり、消費者がいない定常状態のリソースの存在量) を最小化するの制約なしの値は、 リソース供給ベクトルとして知られています。
q
i
α
(
R
)
=
−
a
i
(
R
)
b
α
(
R
)
∂
g
i
∂
R
α
,
{\displaystyle q_{i\alpha }(\mathbf {R} )=-a_{i}(\mathbf {R} )b_{\alpha }(\mathbf {R} ){\frac {\partial g_{i}}{\partial R_{\alpha }}},}
a
i
{\displaystyle a_{i}}
b
α
{\displaystyle b_{\alpha }}
d
(
R
)
{\displaystyle d(\mathbf {R} )}
∂
d
∂
R
α
=
−
h
α
(
R
)
b
α
(
R
)
.
{\displaystyle {\frac {\partial d}{\partial R_{\alpha }}}=-{\frac {h_{\alpha }(\mathbf {R} )}{b_{\alpha }(\mathbf {R} )}}.}
d
{\displaystyle d}
min
R
d
(
R
)
s.t.,
g
i
(
R
)
≤
0
,
i
=
1
,
…
,
S
,
R
α
≥
0
,
α
=
1
,
…
,
M
.
{\displaystyle {\begin{aligned}\min _{\mathbf {R} }&\;d(\mathbf {R} )&&\\{\text{s.t.,}}&\;g_{i}(\mathbf {R} )\leq 0,&&\qquad i=1,\dots ,S,\\&\;R_{\alpha }\geq 0,&&\qquad \alpha =1,\dots ,M.\end{aligned}}}
N
i
{\displaystyle N_{i}}
0
=
N
i
g
i
(
R
)
,
i
=
1
,
…
,
S
,
0
=
∂
d
∂
R
α
−
∑
i
=
1
S
N
i
∂
g
i
∂
R
α
,
α
=
1
,
…
,
M
,
0
≥
g
i
(
R
)
,
i
=
1
,
…
,
S
,
0
≤
N
i
,
i
=
1
,
…
,
S
.
{\displaystyle {\begin{aligned}0&=N_{i}g_{i}(\mathbf {R} ),&&\qquad i=1,\dots ,S,\\0&={\frac {\partial d}{\partial R_{\alpha }}}-\sum _{i=1}^{S}N_{i}{\frac {\partial g_{i}}{\partial R_{\alpha }}},&&\qquad \alpha =1,\dots ,M,\\0&\geq g_{i}(\mathbf {R} ),&&\qquad i=1,\dots ,S,\\0&\leq N_{i},&&\qquad i=1,\dots ,S.\end{aligned}}}
N
¯
i
>
0
{\displaystyle {\overline {N}}_{i}>0}
g
i
(
R
)
{\displaystyle g_{i}(\mathbf {R} )}
N
¯
i
=
0
{\displaystyle {\overline {N}}_{i}=0}
g
i
(
R
)
{\displaystyle g_{i}(\mathbf {R} )}
i
{\displaystyle i}
d
(
R
)
{\displaystyle d(\mathbf {R} )}
R
0
{\displaystyle \mathbf {R} _{0}}
L
(
R
,
{
N
i
}
)
=
d
(
R
)
−
∑
i
=
1
S
N
i
g
i
(
R
)
.
{\displaystyle L(\mathbf {R} ,\{N_{i}\})=d(\mathbf {R} )-\sum _{i=1}^{S}N_{i}g_{i}(\mathbf {R} ).}
R
{\displaystyle \mathbf {R} }
d
(
R
)
{\displaystyle d(\mathbf {R} )}
幾何学的視点
消費者資源モデルの定常状態は、資源豊富度の空間における幾何平均を用いて分析することができる。 [16] [17] [8]
ゼロ純成長等傾斜線(ZNGI)
コミュニティが不可視条件と定常状態条件を満たすためには、定常状態の資源存在量( と表記 )が
すべての種 についてを満たす必要があります 。不等式は、種 が 生き残る場合にのみ飽和します。これらの各条件は、可能な定常状態の資源存在量の空間内の領域を指定し、実現される定常状態の資源存在量はこれらの領域の交差に制限されます。 で指定されるこれらの領域の境界は、 ゼロ純成長等 傾斜線 (ZNGI)として知られています。種が 生き残る場合、定常状態の資源存在量は を満たす必要があります 。したがって、ZNGI の交差点の構造と位置によって、どの種とが共存可能かが決まります。実現される定常状態のコミュニティは、資源の供給に依存しており、共存円錐を調べることで分析できます。
R
⋆
{\displaystyle \mathbf {R} ^{\star }}
g
i
(
R
⋆
)
≤
0
,
{\displaystyle g_{i}(\mathbf {R} ^{\star })\leq 0,}
i
{\displaystyle i}
i
{\displaystyle i}
g
i
(
R
⋆
)
=
0
{\displaystyle g_{i}(\mathbf {R} ^{\star })=0}
i
=
1
,
…
,
S
⋆
{\displaystyle i=1,\dots ,S^{\star }}
g
1
(
R
⋆
)
,
…
,
g
S
⋆
(
R
⋆
)
=
0
{\displaystyle g_{1}(\mathbf {R} ^{\star }),\ldots ,g_{S^{\star }}(\mathbf {R} ^{\star })=0}
共存コーン
ZNGI交差の構造は、どの種が共存可能かを決定しますが、どのような共存種集合が実現するかを決定するものではありません。共存円錐は、資源供給ベクトルが与えられた場合、生態系においてどの種がどの種を生き残らせるかを決定します。ある種集合によって生成される共存円錐は、 まさにその種を含むコミュニティにつながる可能性のある資源供給ベクトルの集合として定義されます 。
i
=
1
,
…
,
S
⋆
{\displaystyle i=1,\ldots ,S^{\star }}
i
=
1
,
…
,
S
⋆
{\displaystyle i=1,\ldots ,S^{\star }}
円錐構造を理解するために、マッカーサーモデルまたはティルマンモデルにおいて、定常状態の非枯渇資源存在量は を満たす必要があることを考慮してください。 ここで、 は収容力/供給率を含むベクトルであり、 は 消費行列 の 番目の行 であり 、ベクトルとして扱われます。生き残る種は正の存在量を持つ種だけなので、和の項は生き残る種についてのみの和となり、右辺は 頂点 を持ち 、生成ベクトルが 生き残る種 の である 凸円錐 の式に似ています。
K
=
R
⋆
+
∑
i
=
1
S
N
i
C
i
,
{\displaystyle \mathbf {K} =\mathbf {R} ^{\star }+\sum _{i=1}^{S}N_{i}\mathbf {C} _{i},}
K
{\displaystyle \mathbf {K} }
C
i
=
(
c
i
1
,
…
,
c
i
M
⋆
)
{\displaystyle \mathbf {C} _{i}=(c_{i1},\ldots ,c_{iM^{\star }})}
i
{\displaystyle i}
c
i
α
{\displaystyle c_{i\alpha }}
R
⋆
{\displaystyle \mathbf {R} ^{\star }}
C
i
{\displaystyle \mathbf {C} _{i}}
i
{\displaystyle i}
複雑な生態系
多くの種と資源が存在する生態系では、消費者-資源モデルの挙動は統計物理学、特に 平均場理論 と キャビティ法 のツールを用いて分析できる。 [18] [19] [20] 大規模な生態系の限界では、パラメータの数が爆発的に増加する。例えば、マッカーサーモデルでは パラメータが必要となる。この限界では、パラメータは定常状態の個体数の分布につながる分布から引き出されると考えられる。これらの定常状態の個体数の分布は、ランダムに選択された種と資源の定常状態の個体数を表す確率変数の平均場方程式を導出することで決定できる。
O
(
S
M
)
{\displaystyle O(SM)}
マッカーサー消費者資源モデル空洞解
MCRM では、モデル パラメータは平均と分散を持つランダム変数として扱われます。
⟨
c
i
α
⟩
=
μ
/
M
,
var
(
c
i
α
)
=
σ
2
/
M
,
⟨
m
i
⟩
=
m
,
var
(
m
i
)
=
σ
m
2
,
⟨
K
α
⟩
=
K
,
var
(
K
α
)
=
σ
K
2
.
{\displaystyle \langle c_{i\alpha }\rangle =\mu /M,\quad \operatorname {var} (c_{i\alpha })=\sigma ^{2}/M,\quad \langle m_{i}\rangle =m,\quad \operatorname {var} (m_{i})=\sigma _{m}^{2},\quad \langle K_{\alpha }\rangle =K,\quad \operatorname {var} (K_{\alpha })=\sigma _{K}^{2}.}
このパラメータ化により、熱力学的極限(すなわち、 )において 、 定常状態の資源および種の豊富さはランダム変数 としてモデル化され、 自己無撞着な平均場方程式 [18] を満たします。ここで 、 はすべて自己無撞着に決定されるモーメントであり、 は独立した 標準正規 ランダム変数であり、 および は 平均感受性であり、これらも自己無撞着に決定される。
M
,
S
→
∞
{\displaystyle M,S\to \infty }
S
/
M
=
Θ
(
1
)
{\displaystyle S/M=\Theta (1)}
N
,
R
{\displaystyle N,R}
0
=
R
(
K
−
μ
S
M
⟨
N
⟩
−
R
+
σ
K
2
+
S
M
σ
2
⟨
N
2
⟩
Z
R
+
σ
2
S
M
ν
R
)
,
0
=
N
(
μ
⟨
R
⟩
−
m
−
σ
2
χ
N
+
σ
2
⟨
R
2
⟩
+
σ
m
2
Z
N
)
,
{\displaystyle {\begin{aligned}0&=R(K-\mu {\tfrac {S}{M}}\langle N\rangle -R+{\sqrt {\sigma _{K}^{2}+{\tfrac {S}{M}}\sigma ^{2}\langle N^{2}\rangle }}Z_{R}+\sigma ^{2}{\tfrac {S}{M}}\nu R),\\0&=N(\mu \langle R\rangle -m-\sigma ^{2}\chi N+{\sqrt {\sigma ^{2}\langle R^{2}\rangle +\sigma _{m}^{2}}}Z_{N}),\end{aligned}}}
⟨
N
⟩
,
⟨
N
2
⟩
,
⟨
R
⟩
,
⟩
R
2
⟩
{\displaystyle \langle N\rangle ,\langle N^{2}\rangle ,\langle R\rangle ,\rangle R^{2}\rangle }
Z
R
,
Z
N
{\displaystyle Z_{R},Z_{N}}
ν
=
⟨
∂
N
/
∂
m
⟩
{\displaystyle \nu =\langle \partial N/\partial m\rangle }
χ
=
⟨
∂
R
/
∂
K
⟩
{\displaystyle \chi =\langle \partial R/\partial K\rangle }
この平均場フレームワークは、存在量分布の瞬間と正確な形、平均感受性、および定常状態で生き残る種と資源の割合を決定することができます。
同様の平均場解析は、外部供給資源モデル [11] 、ティルマンモデル [12] 、微生物消費者資源モデル [20]に対しても行われてきました。これらの手法は 、ランダム一般化ロトカ・ヴォルテラモデル を解析するために初めて開発されました 。
参照
参考文献
^ チェイス, ジョナサン・M.; ライボルド, マシュー・A. (2003). 『生態学的ニッチ 』 シカゴ大学出版局. doi :10.7208/chicago/9780226101811.001.0001. ISBN 978-0-226-10180-4 。
^ ピム, スチュアート・L. (1983年9月). 「ティルマン, D. 1982. 資源競争と群集構造. Monogr. Pop. Biol. 17. プリンストン大学出版局, プリンストン, NJ. 296ページ. $27.50. 陸水学・海洋学 . 28 (5): 1043– 1045. 書誌コード :1983LimOc..28.1043P. doi :10.4319/lo.1983.28.5.1043.
^ マッカーサー, ロバート (1970-05-01). 「種のパッキングと多くの種における競争平衡」. 理論個体群生物学 . 1 (1): 1– 11. Bibcode :1970TPBio...1....1M. doi :10.1016/0040-5809(70)90039-0. ISSN 0040-5809. PMID 5527624.
^ レヴィンズ、リチャード (1968). 『変化する環境における進化:理論的探究』(MPB-2) . プリンストン大学出版局. doi :10.2307/j.ctvx5wbbh. ISBN 978-0-691-07959-2 . JSTOR j.ctvx5wbbh.
^ Goldford, Joshua E.; Lu, Nanxi; Bajić, Djordje; Estrela, Sylvie; Tikhonov, Mikhail; Sanchez-Gorostiaga, Alicia; Segrè, Daniel; Mehta, Pankaj; Sanchez, Alvaro (2018-08-03). 「微生物コミュニティの構築における創発的な単純性」. Science . 361 (6401): 469– 474. Bibcode :2018Sci...361..469G. doi :10.1126/science.aat1168. ISSN 0036-8075. PMC 6405290. PMID 30072533 .
^ Dal Bello, Martina; Lee, Hyunseok; Goyal, Akshit; Gore, Jeff (2021年10月). 「細菌群集における資源と多様性の関係は微生物代謝のネットワーク構造を反映する」 . Nature Ecology & Evolution . 5 (10): 1424– 1434. Bibcode :2021NatEE...5.1424D. doi :10.1038/s41559-021-01535-8. hdl : 1721.1/141887 . ISSN 2397-334X. PMID 34413507. S2CID 256708107.
^ abcd Marsland, Robert; Cui, Wenping; Mehta, Pankaj (2020-09-01). 「最小環境擾乱原理:ニッチ理論への新たな視点」. The American Naturalist . 196 (3): 291– 305. arXiv : 1901.09673 . Bibcode :2020ANat..196..291M. doi :10.1086/710093. ISSN 0003-0147. PMID 32813998. S2CID 59316948.
^ abcdef Cui, Wenping; Marsland III, Robert; Mehta, Pankaj (2024-03-08). 「レ・ズッシュ講義:群集生態学:ニッチ理論から統計力学へ」 arXiv : 2403.05497 [q-bio.PE].
^ 「資源競争とコミュニティ構造(MPB-17)、第17巻 | プリンストン大学出版局」 press.princeton.edu 1982年8月21日 2024年3月18日 閲覧 。
^ チェイス、ジョナサン・M.;ライボルド、マシュー・A.「生態学的ニッチ:古典的アプローチと現代的アプローチの連携」種間相互作用、シカゴ、イリノイ州:シカゴ大学出版局。
^ abc Cui, Wenping; Marsland, Robert; Mehta, Pankaj (2020-07-21). 「多様な生態系における種のパッキングに対する資源動態の影響」. Physical Review Letters . 125 (4) 048101. arXiv : 1911.02595 . Bibcode :2020PhRvL.125d8101C. doi :10.1103/PhysRevLett.125.048101. PMC 8999492. PMID 32794828 .
^ ab Cui, Wenping; Marsland, Robert; Mehta, Pankaj (2021-09-27). 「多様なコミュニティは典型的なランダムエコシステムのように振る舞う」. Physical Review E. 104 ( 3) 034416. Bibcode :2021PhRvE.104c4416C. doi :10.1103/PhysRevE.104.034416. PMC 9005152. PMID 34654170 .
^ Marsland, Robert; Cui, Wenping; Mehta, Pankaj (2020-02-24). 「微生物多様性の最小モデルは実験的に観察された生態学的パターンを再現できる」. Scientific Reports . 10 (1): 3308. arXiv : 1904.12914 . Bibcode :2020NatSR..10.3308M. doi :10.1038/s41598-020-60130-2. ISSN 2045-2322. PMC 7039880. PMID 32094388 .
^アーサー、ロバート ・ マック(1969年12月)「種のパッキングと競争が最小化するもの」 米国科学アカデミー紀要 64 (4): 1369–1371 . doi : 10.1073/pnas.64.4.1369 . ISSN 0027-8424. PMC 223294. PMID 16591810 .
^ マッカーサー, ロバート (1970-05-01). 「種のパッキングと多くの種における競争平衡」. 理論個体群生物学 . 1 (1): 1– 11. Bibcode :1970TPBio...1....1M. doi :10.1016/0040-5809(70)90039-0. ISSN 0040-5809. PMID 5527624.
^ Mancuso, Christopher P; Lee, Hyunseok; Abreu, Clare I; Gore, Jeff; Khalil, Ahmad S (2021-09-03). Shou, Wenying; Walczak, Aleksandra M; Shou, Wenying (編). 「環境変動は微生物群集における予期せぬ多様性と撹乱の関係を再構築する」. eLife . 10 e67175. doi : 10.7554/eLife.67175 . ISSN 2050-084X. PMC 8460265. PMID 34477107 .
^ Tikhonov, Mikhail; Monasson, Remi (2017-01-27). 「高度に多様な生態系における資源競争の集団的段階」. Physical Review Letters . 118 (4) 048103. arXiv : 1609.01270 . Bibcode :2017PhRvL.118d8103T. doi :10.1103/PhysRevLett.118.048103. PMID 28186794.
^ ab Advani, Madhu; Bunin, Guy; Mehta, Pankaj (2018-03-20). 「群集生態学の統計物理学:マッカーサーの消費者資源モデルに対する空洞解」. Journal of Statistical Mechanics: Theory and Experiment . 2018 (3): 033406. Bibcode :2018JSMTE..03.3406A. doi :10.1088/1742-5468/aab04e. ISSN 1742-5468. PMC 6329381. PMID 30636966 .
^ Marsland, Robert; Cui, Wenping; Mehta, Pankaj (2020-02-24). 「微生物多様性の最小モデルは実験的に観察された生態学的パターンを再現できる」. Scientific Reports . 10 (1): 3308. arXiv : 1904.12914 . Bibcode :2020NatSR..10.3308M. doi :10.1038/s41598-020-60130-2. ISSN 2045-2322. PMC 7039880. PMID 32094388 .
^ ab Mehta, Pankaj; Marsland III, Robert (2021-10-10), Cross-feeding shapes both competition and collaboration in microbial ecosystems , arXiv : 2110.04965
さらに読む
Cui, Wenping; Robert Marsland III; Mehta, Pankaj (2024). 「レ・ズッシュ講義:群集生態学:ニッチ理論から統計力学へ」 arXiv : 2403.05497 [q-bio.PE].
Stefano Allesina のコミュニティ生態学コースの講義ノート: https://stefanoallesina.github.io/Theoretical_Community_Ecology/