Jump to content

Consumer-resource model: Difference between revisions

From Wikipedia, the free encyclopedia
Content deleted Content added
Citation bot (talk | contribs)
Alter: template type. Add: eprint, class. Removed parameters. Some additions/deletions were parameter name changes. | Use this bot. Report bugs. | #UCB_CommandLine
correcting according to "Articles for possible copyedit from 2024-06-20 dump" – "equations,where", "as,where", "equation,species", "is,In", "satisfy,for", "satisfy,where", "equations,where", "uninvasibility"; removing {{Copy edit|date=June 2024}}
 
(42 intermediate revisions by 12 users not shown)
Line 1: Line 1:
{{Short description|Class of ecological models}}
{{AfC submission|t||ts=20230912152953|u=Magenta.lily|ns=118|demo=}}<!-- Important, do not remove this line before article has been created. -->

In [[theoretical ecology]] and [[nonlinear dynamics]], '''consumer-resource models''' (CRMs) are a class of [[Ecological model|ecological models]] in which a [[Community (ecology)|community]] of consumer species compete for a common pool of resources. Instead of species interacting directly, all species-species interactions are mediated through resource dynamics. Consumer-resource models have served as fundamental tools in the quantitative development of theories of [[Ecological niche|niche construction]], [[Coexistence theory|coexistence]], and [[Biodiversity|biological diversity]]. These models can be interpreted as a quantitative description of a single [[trophic level]].<ref>{{Cite book |last1=Chase |first1=Jonathan M. |url=https://doi.org/10.7208/chicago/9780226101811.001.0001 |title=Ecological Niches |last2=Leibold |first2=Mathew A. |date=2003 |publisher=University of Chicago Press |isbn=978-0-226-10180-4 |language=en |doi=10.7208/chicago/9780226101811.001.0001}}</ref><ref>{{Cite journal |last=Pimm |first=Stuart L. |date=September 1983 |title=TILMAN, D. 1982. Resource competition and community structure. Monogr. Pop. Biol. 17. Princeton University Press, Princeton, N.J. 296 p. $27.50. |url=http://doi.wiley.com/10.4319/lo.1983.28.5.1043 |journal=Limnology and Oceanography |language=en |volume=28 |issue=5 |pages=1043–1045 |doi=10.4319/lo.1983.28.5.1043|bibcode=1983LimOc..28.1043P }}</ref>



In [[theoretical ecology]] and [[nonlinear dynamics]], '''consumer-resource models''' (CRMs), also known as '''resource-competition models''', are a class of [[Ordinary differential equation|ordinary differential equations]]-based [[Ecological model|ecological models]] in which a [[Community (ecology)|community]] of consumer species compete for a common pool of resources. Instead of species interacting directly, all species-species interactions are mediated through resource dynamics. Consumer-resource models have served as fundamental tools in the quantitative development of theories of [[Ecological niche|niche construction]], [[Coexistence theory|coexistence]], and [[Biodiversity|biological diversity]]. These models can be interpreted as a quantitative description of a single (or multiple) [[trophic level|trophic level(s)]].<ref>{{Cite book |last1=Chase |first1=Jonathan M. |url=https://doi.org/10.7208/chicago/9780226101811.001.0001 |title=Ecological Niches |last2=Leibold |first2=Mathew A. |date=2003 |publisher=University of Chicago Press |isbn=978-0-226-10180-4 |language=en |doi=10.7208/chicago/9780226101811.001.0001}}</ref><ref>{{Cite journal |last=Pimm |first=Stuart L. |date=September 1983 |title=TILMAN, D. 1982. Resource competition and community structure. Monogr. Pop. Biol. 17. Princeton University Press, Princeton, N.J. 296 p. $27.50. |url=http://doi.wiley.com/10.4319/lo.1983.28.5.1043 |journal=Limnology and Oceanography |language=en |volume=28 |issue=5 |pages=1043–1045 |doi=10.4319/lo.1983.28.5.1043|bibcode=1983LimOc..28.1043P }}</ref><ref>{{Cite book |last=Levins |first=Richard |url=https://www.jstor.org/stable/j.ctvx5wbbh |title=Evolution in Changing Environments: Some Theoretical Explorations. (MPB-2) |date=1968 |publisher=Princeton University Press |doi=10.2307/j.ctvx5wbbh |jstor=j.ctvx5wbbh |isbn=978-0-691-07959-2}}</ref><ref name=":4">{{Cite journal |last1=Mancuso |first1=Christopher P |last2=Lee |first2=Hyunseok |last3=Abreu |first3=Clare I |last4=Gore |first4=Jeff |last5=Khalil |first5=Ahmad S |date=2021-09-03 |editor-last=Shou |editor-first=Wenying |editor2-last=Walczak |editor2-first=Aleksandra M |editor3-last=Shou |editor3-first=Wenying |title=Environmental fluctuations reshape an unexpected diversity-disturbance relationship in a microbial community |journal=eLife |volume=10 |pages=e67175 |doi=10.7554/eLife.67175 |issn=2050-084X |pmc=8460265 |pmid=34477107 |doi-access=free }}</ref><ref name=":3">{{Cite journal |last1=Dal Bello |first1=Martina |last2=Lee |first2=Hyunseok |last3=Goyal |first3=Akshit |last4=Gore |first4=Jeff |date=October 2021 |title=Resource–diversity relationships in bacterial communities reflect the network structure of microbial metabolism |url=https://www.nature.com/articles/s41559-021-01535-8 |journal=Nature Ecology & Evolution |language=en |volume=5 |issue=10 |pages=1424–1434 |doi=10.1038/s41559-021-01535-8 |pmid=34413507 |s2cid=256708107 |issn=2397-334X}}</ref><ref>{{Cite journal |last1=Posfai |first1=Anna |last2=Taillefumier |first2=Thibaud |last3=Wingreen |first3=Ned S. |date=2017-01-12 |title=Metabolic Trade-Offs Promote Diversity in a Model Ecosystem |journal=Physical Review Letters |volume=118 |issue=2 |pages=028103 |doi=10.1103/PhysRevLett.118.028103 |pmc=5743855 |pmid=28128613|bibcode=2017PhRvL.118b8103P }}</ref><ref>{{Cite book |last=Stevens |first=Hank |url=https://hankstevens.github.io/Primer-of-Ecology/index.html |title=Preface {{!}} Primer of Ecology using R |language=en}}</ref>


A general consumer-resource model consists of ''{{mvar|M}}'' resources whose abundances are <math>R_1,\dots,R_M</math> and ''{{mvar|S}}'' consumer species whose populations are <math>N_1,\dots,N_S</math>. A general consumer-resource model is described by the system of coupled ordinary differential equations,<math display="block">
A general consumer-resource model consists of ''{{mvar|M}}'' resources whose abundances are <math>R_1,\dots,R_M</math> and ''{{mvar|S}}'' consumer species whose populations are <math>N_1,\dots,N_S</math>. A general consumer-resource model is described by the system of coupled ordinary differential equations,<math display="block">
Line 14: Line 16:
f_\alpha(R_1,\dots,R_M,N_1,\dots,N_S), &&\qquad \alpha = 1,\dots,M
f_\alpha(R_1,\dots,R_M,N_1,\dots,N_S), &&\qquad \alpha = 1,\dots,M
\end{align}
\end{align}
</math>
</math>where <math>g_i</math>, depending only on resource abundances, is the per-capita growth rate of species <math>i</math>, and <math>f_\alpha</math> is the growth rate of resource <math>\alpha</math>. An essential feature of CRMs is that species growth rates and populations are mediated through resources and there are no explicit species-species interactions. Despite this, there remains implicit dependence due to the interacting nature of the system.
where <math>g_i</math>, depending only on resource abundances, is the per-capita growth rate of species <math>i</math>, and <math>f_\alpha</math> is the growth rate of resource <math>\alpha</math>. An essential feature of CRMs is that species growth rates and populations are mediated through resources and there are no explicit species-species interactions. Through resource interactions, there are [[Emergence|emergent]] inter-species interactions.


Originally introduced by [[Robert H. MacArthur]]<ref name=":11">{{Cite journal |last=MacArthur |first=Robert |date=1970-05-01 |title=Species packing and competitive equilibrium for many species |url=https://dx.doi.org/10.1016/0040-5809%2870%2990039-0 |journal=Theoretical Population Biology |volume=1 |issue=1 |pages=1–11 |doi=10.1016/0040-5809(70)90039-0 |issn=0040-5809 |pmid=5527624}}</ref> and [[Richard Levins]],<ref name=":10">{{Cite book |last=Levins |first=Richard |url=https://www.jstor.org/stable/j.ctvx5wbbh |title=Evolution in Changing Environments: Some Theoretical Explorations. (MPB-2) |date=1968 |publisher=Princeton University Press |isbn=978-0-691-07959-2 |doi=10.2307/j.ctvx5wbbh |jstor=j.ctvx5wbbh}}</ref> consumer-resource models have found success in formalizing ecological principles and modeling experiments involving microbial ecosystems.<ref>{{Cite journal |last1=Goldford |first1=Joshua E. |last2=Lu |first2=Nanxi |last3=Bajić |first3=Djordje |last4=Estrela |first4=Sylvie |last5=Tikhonov |first5=Mikhail |last6=Sanchez-Gorostiaga |first6=Alicia |last7=Segrè |first7=Daniel |last8=Mehta |first8=Pankaj |last9=Sanchez |first9=Alvaro |date=2018-08-03 |title=Emergent simplicity in microbial community assembly |journal=Science |language=en |volume=361 |issue=6401 |pages=469–474 |doi=10.1126/science.aat1168 |issn=0036-8075 |pmc=6405290 |pmid=30072533|bibcode=2018Sci...361..469G }}</ref><ref name=":3">{{Cite journal |last1=Dal Bello |first1=Martina |last2=Lee |first2=Hyunseok |last3=Goyal |first3=Akshit |last4=Gore |first4=Jeff |date=October 2021 |title=Resource–diversity relationships in bacterial communities reflect the network structure of microbial metabolism |url=https://www.nature.com/articles/s41559-021-01535-8 |journal=Nature Ecology & Evolution |language=en |volume=5 |issue=10 |pages=1424–1434 |bibcode=2021NatEE...5.1424D |doi=10.1038/s41559-021-01535-8 |issn=2397-334X |pmid=34413507 |s2cid=256708107 |hdl-access=free |hdl=1721.1/141887}}</ref>
== History ==
Add a bunch... Developed and introduced significantly by [[Robert H. MacArthur]] and [[Richard Levins]].
== Models ==
[[File:Consumer Resource Models Resource Dynamics Schematics.png|thumb]]


== Models ==
=== Niche models ===
=== Niche models ===
Niche models are a notable class of CRMs which are described by the system of coupled ordinary differential equations,<ref name=":0">{{Cite journal |last1=Marsland |first1=Robert |last2=Cui |first2=Wenping |last3=Mehta |first3=Pankaj |date=2020-09-01 |title=The Minimum Environmental Perturbation Principle: A New Perspective on Niche Theory |url=http://dx.doi.org/10.1086/710093 |journal=The American Naturalist |volume=196 |issue=3 |pages=291–305 |doi=10.1086/710093 |pmid=32813998 |s2cid=59316948 |issn=0003-0147}}</ref>
Niche models are a notable class of CRMs which are described by the system of coupled ordinary differential equations,<ref name=":0">{{Cite journal |last1=Marsland |first1=Robert |last2=Cui |first2=Wenping |last3=Mehta |first3=Pankaj |date=2020-09-01 |title=The Minimum Environmental Perturbation Principle: A New Perspective on Niche Theory |url=http://dx.doi.org/10.1086/710093 |journal=The American Naturalist |volume=196 |issue=3 |pages=291–305 |doi=10.1086/710093 |pmid=32813998 |s2cid=59316948 |issn=0003-0147|arxiv=1901.09673 }}</ref><ref name=":5">{{cite arXiv |eprint=2403.05497 |class=q-bio.PE |first1=Wenping |last1=Cui |first2=Robert |last2=Marsland III |title=Les Houches Lectures on Community Ecology: From Niche Theory to Statistical Mechanics |date=2024-03-08 |last3=Mehta |first3=Pankaj}}</ref>


: <math>
: <math>
Line 39: Line 40:


==== MacArthur consumer-resource model (MCRM) ====
==== MacArthur consumer-resource model (MCRM) ====
The MacArthur consumer-resource model (MCRM), named after [[Robert H. MacArthur]], is a foundational CRM for the development of niche and coexistence theories.<small>[citation needed]</small> The MCRM is given by the following set of coupled ordinary differential equations:<ref name=":1">{{Cite journal |last1=Cui |first1=Wenping |last2=Marsland |first2=Robert |last3=Mehta |first3=Pankaj |date=2020-07-21 |title=Effect of Resource Dynamics on Species Packing in Diverse Ecosystems |journal=Physical Review Letters |volume=125 |issue=4 |pages=048101 |doi=10.1103/PhysRevLett.125.048101 |pmc=8999492 |pmid=32794828|arxiv=1911.02595 |bibcode=2020PhRvL.125d8101C }}</ref><ref name=":2">{{Cite journal |last1=Cui |first1=Wenping |last2=Marsland |first2=Robert |last3=Mehta |first3=Pankaj |date=2021-09-27 |title=Diverse communities behave like typical random ecosystems |journal=Physical Review E |volume=104 |issue=3 |pages=034416 |doi=10.1103/PhysRevE.104.034416 |pmc=9005152 |pmid=34654170|bibcode=2021PhRvE.104c4416C }}</ref><math display="block">\begin{align}
The MacArthur consumer-resource model (MCRM), named after [[Robert H. MacArthur]], is a foundational CRM for the development of niche and coexistence theories.<ref>{{Cite web |date=1982-08-21 |title=Resource Competition and Community Structure. (MPB-17), Volume 17 {{!}} Princeton University Press |url=https://press.princeton.edu/books/paperback/9780691083025/resource-competition-and-community-structure-mpb-17-volume-17 |access-date=2024-03-18 |website=press.princeton.edu |language=en}}</ref><ref>{{Cite book |last1=Chase |first1=Jonathan M. |url=https://press.uchicago.edu/ucp/books/book/chicago/E/bo3638660.html |title=Ecological Niches: Linking Classical and Contemporary Approaches |last2=Leibold |first2=Mathew A. |publisher=University of Chicago Press |series=Interspecific Interactions |location=Chicago, IL |language=en}}</ref> The MCRM is given by the following set of coupled ordinary differential equations:<ref name=":1">{{Cite journal |last1=Cui |first1=Wenping |last2=Marsland |first2=Robert |last3=Mehta |first3=Pankaj |date=2020-07-21 |title=Effect of Resource Dynamics on Species Packing in Diverse Ecosystems |journal=Physical Review Letters |volume=125 |issue=4 |pages=048101 |doi=10.1103/PhysRevLett.125.048101 |pmc=8999492 |pmid=32794828|arxiv=1911.02595 |bibcode=2020PhRvL.125d8101C }}</ref><ref name=":2">{{Cite journal |last1=Cui |first1=Wenping |last2=Marsland |first2=Robert |last3=Mehta |first3=Pankaj |date=2021-09-27 |title=Diverse communities behave like typical random ecosystems |journal=Physical Review E |volume=104 |issue=3 |pages=034416 |doi=10.1103/PhysRevE.104.034416 |pmc=9005152 |pmid=34654170|bibcode=2021PhRvE.104c4416C }}</ref><ref name=":5" /><math display="block">\begin{align}
\frac{\mathrm dN_i}{\mathrm dt} &=
\frac{\mathrm dN_i}{\mathrm dt} &=
\tau_i^{-1} N_i \left(
\tau_i^{-1} N_i \left(
Line 54: Line 55:
&&
&&
\qquad \alpha = 1,\dots,M,
\qquad \alpha = 1,\dots,M,
\end{align}</math>where <math>c_{i\alpha}</math> is the relative preference of species <math>i</math> for resource <math>\alpha</math> and also the relative amount by which resource <math>\alpha</math> is depleted by the consumption of consumer species <math>i</math>; <math>K_\alpha</math> is the steady-state carrying capacity of resource <math>\alpha</math> in absence of consumption (i.e., when <math>c_{i\alpha}</math> is zero); <math>\tau_i</math> and <math>r_\alpha^{-1}</math> are time-scales for species and resource dynamics, respectively; <math>w_\alpha</math> is the quality of resource <math>\alpha</math>; and <math>m_i</math> is the natural mortality rate of species <math>i</math>. This model is said to have self-replenishing resource dynamics because when <math>c_{i\alpha} = 0</math>, each resource exhibits independent [[logistic growth]]. Given positive parameters and initial conditions, this model approaches a unique [[Invasibility|uninvadable]] [[steady state]] (i.e., a steady state in which the re-introduction of a species which has been driven to extinction or a resource which has been depleted leads to the re-introduced species or resource dying out again).<ref name=":0" /> Steady states of the MCRM satisfy the [[competitive exclusion principle]]: the number of coexisting species is less than or equal to the number of non-depleted resources. In other words, the number of simultaneously-occupiable [[Ecological niche|ecological niches]] is equal to the number of non-depleted resources.
\end{align}</math>where <math>c_{i\alpha}</math> is the relative preference of species <math>i</math> for resource <math>\alpha</math> and also the relative amount by which resource <math>\alpha</math> is depleted by the consumption of consumer species <math>i</math>; <math>K_\alpha</math> is the steady-state carrying capacity of resource <math>\alpha</math> in absence of consumption (i.e., when <math>c_{i\alpha}</math> is zero); <math>\tau_i</math> and <math>r_\alpha^{-1}</math> are time-scales for species and resource dynamics, respectively; <math>w_\alpha</math> is the quality of resource <math>\alpha</math>; and <math>m_i</math> is the natural mortality rate of species <math>i</math>. This model is said to have self-replenishing resource dynamics because when <math>c_{i\alpha} = 0</math>, each resource exhibits independent [[logistic growth]]. Given positive parameters and initial conditions, this model approaches a unique [[Invasibility|uninvadable]] [[steady state]] (i.e., a steady state in which the re-introduction of a species which has been driven to extinction or a resource which has been depleted leads to the re-introduced species or resource dying out again).<ref name=":0" /> Steady states of the MCRM satisfy the [[competitive exclusion principle]]: the number of coexisting species is less than or equal to the number of non-depleted resources. In other words, the number of simultaneously occupiable [[Ecological niche|ecological niches]] is equal to the number of non-depleted resources.


==== Externally supplied resources model ====
Occasionally, the MacArthur consumer resource model is extended to the asymmetric MacArthur consumer resource model (aMCRM) by changing the relative weights by which resource abundances are depleted through consumption:<ref>{{Cite arXiv |last1=Blumenthal |first1=Emmy |last2=Rocks |first2=Jason W. |last3=Mehta |first3=Pankaj |date=2023 |title=Phase transition to chaos in complex ecosystems with non-reciprocal species-resource interactions |class=physics.bio-ph |eprint=2308.15757}}</ref><math display="block">\frac{\mathrm dR_\alpha}{\mathrm dt} =
The externally supplied resource model is similar to the MCRM except the resources are provided at a constant rate from an external source instead of being self-replenished. This model is also sometimes called the linear resource dynamics model. It is described by the following set of coupled ordinary differential equations:<ref name=":1" /><ref name=":5" /><math display="block">\begin{align}
\frac{r_\alpha}{K_\alpha} \left(
K_\alpha - R_\alpha
\right)R_\alpha
-
\sum_{i=1}^S N_i e_{i\alpha}R_\alpha</math>Restricting <math>e_{i\alpha} = c_{i\alpha}</math> recovers the original MCRM.

==== Externally-supplied resources model (eCRM) ====
The externally-supplied resource model (eCRM) is similar to the MCRM except the resources are provided at a constant rate from an external source instead of being self-replenished. This model is also sometimes called the linear resource dynamics model. It is described by the following set of coupled ordinary differential equations:<ref name=":1" /><ref name=":2" /><math display="block">\begin{align}
\frac{\mathrm dN_i}{\mathrm dt} &=
\frac{\mathrm dN_i}{\mathrm dt} &=
\tau_i^{-1} N_i \left(
\tau_i^{-1} N_i \left(
Line 77: Line 71:
&&
&&
\qquad \alpha = 1,\dots,M,
\qquad \alpha = 1,\dots,M,
\end{align}</math>where all the parameters shared with the MCRM are the same, and <math>\kappa_\alpha</math> is the rate at which resource <math>\alpha</math> is supplied to the ecosystem. In the eCRM, in the absence of consumption, <math>R_\alpha</math> decays to <math>\kappa_\alpha</math> exponentially with timescale <math>r_\alpha^{-1}</math>.
\end{align}</math>where all the parameters shared with the MCRM are the same, and <math>\kappa_\alpha</math> is the rate at which resource <math>\alpha</math> is supplied to the ecosystem. In the eCRM, in the absence of consumption, <math>R_\alpha</math> decays to <math>\kappa_\alpha</math> exponentially with timescale <math>r_\alpha^{-1}</math>. This model is also known as a chemostat model.


==== Tilman consumer-resource model (TCRM) ====
==== Tilman consumer-resource model (TCRM) ====
The Tilman consumer-resource model (TCRM), named after [[G. David Tilman]], is similar to the eCRM except the rate at which a species depletes a resource is no longer proportional to the present abundance of the resource. The TCRM is the foundational model for Tilman's [[R* rule (ecology)|R* rule]]. It is described by the following set of coupled ordinary differential equations:<ref name=":1" /><ref name=":2" /><math display="block">\begin{align}
The Tilman consumer-resource model (TCRM), named after [[G. David Tilman]], is similar to the externally supplied resources model except the rate at which a species depletes a resource is no longer proportional to the present abundance of the resource. The TCRM is the foundational model for Tilman's [[R* rule (ecology)|R* rule]]. It is described by the following set of coupled ordinary differential equations:<ref name=":5" /><math display="block">\begin{align}
\frac{\mathrm dN_i}{\mathrm dt} &=
\frac{\mathrm dN_i}{\mathrm dt} &=
\tau_i^{-1} N_i \left(
\tau_i^{-1} N_i \left(
Line 93: Line 87:
&&
&&
\qquad \alpha = 1,\dots,M,
\qquad \alpha = 1,\dots,M,
\end{align}</math>where all parameters are shared with the MCRM. In the TCRM, resource abundances can become nonphysically negative.
\end{align}</math>where all parameters are shared with the MCRM. In the TCRM, resource abundances can become nonphysically negative.<ref name=":0" />

<ref name=":0" /><ref>{{Cite journal |last=MacArthur |first=Robert |date=1970-05-01 |title=Species packing and competitive equilibrium for many species |url=https://dx.doi.org/10.1016/0040-5809%2870%2990039-0 |journal=Theoretical Population Biology |volume=1 |issue=1 |pages=1–11 |doi=10.1016/0040-5809(70)90039-0 |pmid=5527624 |issn=0040-5809}}</ref>


==== Microbial consumer-resource model (MiCRM) ====
==== Microbial consumer-resource model (MiCRM) ====
The microbial consumer resource model describes a microbial ecosystem with externally-supplied resources where consumption can produce metabolic byproducts, leading to potential [[Cross feeding|cross-feeding]]. It is described by the following set of coupled ODEs:<math display="block">\begin{align}
The microbial consumer resource model describes a microbial ecosystem with externally supplied resources where consumption can produce metabolic byproducts, leading to potential [[Cross feeding|cross-feeding]]. It is described by the following set of coupled ODEs:<math display="block">\begin{align}
\frac{\mathrm dN_i}{\mathrm dt} &=
\frac{\mathrm dN_i}{\mathrm dt} &=
\tau_i^{-1} N_i \left(
\tau_i^{-1} N_i \left(
Line 114: Line 106:
&&
&&
\qquad \alpha = 1,\dots,M,
\qquad \alpha = 1,\dots,M,
\end{align}</math>where all parameters shared with the MCRM have similar interpretations; <math>D_{\alpha\beta}</math> is the fraction of the byproducts due to consumption of resource <math>\beta</math> which are converted to resource <math>\alpha</math> and <math>l_\alpha</math> is the "leakage fraction" of resource <math>\alpha</math> governing how much of the resource is released into the environment as metabolic byproducts.<ref>{{Cite journal |last1=Marsland |first1=Robert |last2=Cui |first2=Wenping |last3=Mehta |first3=Pankaj |date=2020-02-24 |title=A minimal model for microbial biodiversity can reproduce experimentally observed ecological patterns |journal=Scientific Reports |language=en |volume=10 |issue=1 |pages=3308 |doi=10.1038/s41598-020-60130-2 |pmid=32094388 |pmc=7039880 |arxiv=1904.12914 |bibcode=2020NatSR..10.3308M |issn=2045-2322}}</ref><ref name=":0" />
\end{align}</math>where all parameters shared with the MCRM have similar interpretations; <math>D_{\alpha\beta}</math> is the fraction of the byproducts due to consumption of resource <math>\beta</math> which are converted to resource <math>\alpha</math> and <math>l_\alpha</math> is the "leakage fraction" of resource <math>\alpha</math> governing how much of the resource is released into the environment as metabolic byproducts.<ref name=":12">{{Cite journal |last1=Marsland |first1=Robert |last2=Cui |first2=Wenping |last3=Mehta |first3=Pankaj |date=2020-02-24 |title=A minimal model for microbial biodiversity can reproduce experimentally observed ecological patterns |journal=Scientific Reports |language=en |volume=10 |issue=1 |pages=3308 |doi=10.1038/s41598-020-60130-2 |pmid=32094388 |pmc=7039880 |arxiv=1904.12914 |bibcode=2020NatSR..10.3308M |issn=2045-2322}}</ref><ref name=":5" />
=== Multi-trophic Model ===
The multi-trophic consumer-resource model extends the typical paradigm of consumer-resource models which describe single [[Trophic level|trophic levels]] to describe multiple trophic levels. A description of an ecosystem with two-level trophic structure consists of <math>M_R</math> primary producers with abundances <math>\{R_\alpha\}_{\alpha = 1}^{M_R}</math>, <math>M_N</math> primary consumers with populations <math>\{N_i\}_{i=1}^{M_N}</math>, and <math>M_X</math> secondary consumers with populations <math>\{X_k\}_{k=1}^{M_X}</math>. In an example ecosystem, the primary producers could be plants, the primary consumers could be herbivores, and the secondary consumers could be carnivores; in this example, plants self-renew or are supplied externally, herbivores only consume plants, and carnivores only consume herbivores. Such a two-level trophic structure can be described by the following set of coupled ODEs:<ref name=":5">{{Cite arXiv |last1=Feng |first1=Zhijie |last2=Marsland |first2=Robert |last3=Rocks |first3=Jason W. |last4=Mehta |first4=Pankaj |date=2023 |title=Emergent competition shapes the ecological properties of multi-trophic ecosystems |class=q-bio.PE |eprint=2303.02983 }}</ref><math display="block">\begin{align}
\frac{\mathrm dX_k}{\mathrm d t}
&=
X_k \left(
\sum_{i=1}^{M_N} d_{ki}N_i - u_k
\right),
&&\qquad
k = 1 ,\dots,M_X,
\\
\frac{\mathrm d N_i}{\mathrm d t}
&=
N_i
\left(
\sum_{\alpha =1 }^{M_R} c_{i\alpha} R_\alpha - m_i - \sum_{k = 1}^{M_X} d_{ki} X_k\right),
&&\qquad
i = 1,\dots,M_N,\\
\frac{\mathrm dR_\alpha}{\mathrm d t}
&=
R_\alpha \left(
K_\alpha - R_\alpha - \sum_{i=1}^{M_N} c_{i\alpha} N_i
\right),
&&\qquad \alpha = 1,\dots,M_R.
\end{align}</math>The form of the multi-trophic CRM has similar structure to the MCRM, specifically in that the primary producers are self-renewing. This model can be generalized further by allowing for any number of intermediate consumer/predator species where species in the <math>n</math>th trophic level consume only those in the <math>(n-1)</math>th trophic level and are consumed only by those in the <math>(n+1)</math>th trophic level. Such intermediate-level trophic interactions could be described by the following set of coupled ODEs:<ref name=":5" /><math display="block">\begin{align}
\frac{\mathrm dN_i^{(n)}}{\mathrm d t} &= N_i^{(n)} \left(
\sum_{j=1}^{M_{(n-1)}} c^{(n)}_{ij} N^{(n-1)}_j
- m^{(n)}_i
- \sum_{k = 1}^{M_{(n+1)}} c_{ki}^{(n+1)} N_k^{(n+1)}
\right),
\qquad i =1,\dots, M_{(n)}.
\end{align}</math>

== Experimental results ==

== Connection to Lotka–Volterra models ==
The (asymmetric) MacArthur consumer resource model can be written in the form of a [[Generalized Lotka–Volterra equation|generalized Lotka–Volterra model]]:<math display="block">\frac{\mathrm dx_i}{\mathrm dt}
=
x_i
\left(
b_i - \sum_{j=1}^{n} A_{ij} x_j
\right),
\qquad
i = 1,\dots,n.</math>For the aMCRM, there are <math>n = M+S</math> species where <math>x_i</math> is identified with <math>N_i</math> for <math>i =1,\dots,S</math> is identified with <math>R_\alpha</math> for <math>i = S+1,\dots,M+S</math>. In this form, the relevant Lotka–Volterra parameters are,<math display="block">\begin{align}
\left[b_i\right]
&=
\begin{bmatrix}
-\left[m_i\right]\\
\left[K_\alpha\right]\\
\end{bmatrix}
=
\begin{bmatrix}
-m_1\\
\vdots\\
-m_S\\
K_1\\
\vdots\\
K_M
\end{bmatrix},
\\
\left[A_{ij}\right]
&=
\begin{bmatrix}
0 & \left[c_{i\alpha}\right]\\
\left[e_{i\alpha}\right]^{\mathrm T} & -I
\end{bmatrix}
=
\begin{bmatrix}
0 & \cdots & 0 & c_{11} & \cdots & c_{1M}\\
\vdots & \ddots & \vdots & \vdots & \ddots & \vdots\\
0 & \cdots & 0 & c_{S1} & \cdots & c_{SM}\\
-e_{11} & \cdots & -e_{S1} & -1 & \cdots & 0\\
\vdots & \ddots & \vdots & \vdots & \ddots & \vdots \\
-e_{1M} & \cdots & -e_{SM} & 0 & \cdots & -1
\end{bmatrix}.
\end{align}</math>Some but not all consumer resource models can be written as generalized Lotka–Volterra models.


== Symmetric interactions and constrained optimization ==
== Symmetric interactions and optimization ==


=== MacArthur's Minimization Principle ===
=== MacArthur's Minimization Principle ===
For the MacArthur consumer resource model (MCRM), [[Robert H. MacArthur|MacArthur]] introduced an optimization principle to identify the uninvadable steady state of the model (i.e., the steady state so that if any species with zero population is re-introduced, it will fail to invade, meaning the ecosystem will return to said steady state). To derive the optimization principle, one assumes resource dynamics become sufficiently fast (i.e., <math>r_\alpha \gg 1</math>) that they become entrained to species dynamics and are constantly at steady state (i.e., <math>{\mathrm d}R_\alpha/{\mathrm d}t = 0</math>) so that <math>R_\alpha</math> is expressed as a function of <math>N_i</math>. With this assumption, one can express species dynamics as,<math display="block">
For the MacArthur consumer resource model (MCRM), [[Robert H. MacArthur|MacArthur]] introduced an optimization principle to identify the uninvadable steady state of the model (i.e., the steady state so that if any species with zero population is re-introduced, it will fail to invade, meaning the ecosystem will return to said steady state). To derive the optimization principle, one assumes resource dynamics become sufficiently fast (i.e., <math>r_\alpha \gg 1</math>) that they become entrained to species dynamics and are constantly at steady state (i.e., <math>{\mathrm d}R_\alpha/{\mathrm d}t = 0</math>) so that <math>R_\alpha</math> is expressed as a function of <math>N_i</math>. With this assumption, one can express species dynamics as,
<math display="block">
\frac{\mathrm dN_i}{\mathrm dt}
\frac{\mathrm dN_i}{\mathrm dt}
=
=
Line 204: Line 122:
-m_i
-m_i
\right],
\right],
</math>
</math>where <math>\sum_{\alpha \in M^\ast}</math> denotes a sum over resource abundances which satisfy <math>R_\alpha = r_\alpha - \sum_{j=1}^S N_j c_{j\alpha} \geq 0</math>. The above expression can be written as <math>\mathrm{d}N_i/\mathrm{d}t=-\tau_i^{-1}N_i \,\partial Q/\partial N_i</math>, where,<math display="block">
where <math>\sum_{\alpha \in M^\ast}</math> denotes a sum over resource abundances which satisfy <math>R_\alpha = r_\alpha - \sum_{j=1}^S N_j c_{j\alpha} \geq 0</math>. The above expression can be written as <math>\mathrm{d}N_i/\mathrm{d}t=-\tau_i^{-1}N_i \,\partial Q/\partial N_i</math>, where,<math display="block">
Q(\{N_i\})
Q(\{N_i\})
=
=
Line 217: Line 136:
</math>
</math>


At un-invadable steady state <math>\partial Q/\partial N_i = 0</math> for all surviving species <math>i</math> and <math>\partial Q/\partial N_i > 0</math> for all extinct species <math>i</math>.<ref name=":0" /><ref>{{Cite journal |last=Arthur |first=Robert Mac |title=Species Packing, and What Competition Minimizes |date=December 1969 |journal=Proceedings of the National Academy of Sciences |language=en |volume=64 |issue=4 |pages=1369–1371 |doi=10.1073/pnas.64.4.1369 |issn=0027-8424 |pmc=223294 |pmid=16591810 |doi-access=free }}</ref><ref>{{Cite journal |last=MacArthur |first=Robert |date=1970-05-01 |title=Species packing and competitive equilibrium for many species |url=https://dx.doi.org/10.1016/0040-5809%2870%2990039-0 |journal=Theoretical Population Biology |volume=1 |issue=1 |pages=1–11 |doi=10.1016/0040-5809(70)90039-0 |pmid=5527624 |issn=0040-5809}}</ref><ref>{{Cite journal |last=Haygood |first=Ralph |date=March 2002 |title=Coexistence in MacArthur-Style Consumer–Resource Models |url=https://linkinghub.elsevier.com/retrieve/pii/S0040580901915667 |journal=Theoretical Population Biology |language=en |volume=61 |issue=2 |pages=215–223 |doi=10.1006/tpbi.2001.1566|pmid=11969391 }}</ref><ref name=":6">{{Cite journal |last1=Case |first1=Ted J. |last2=Casten |first2=Richard G. |date=1979 |title=Global Stability and Multiple Domains of Attraction in Ecological Systems |url=https://www.jstor.org/stable/2459962 |journal=The American Naturalist |volume=113 |issue=5 |pages=705–714 |doi=10.1086/283427 |jstor=2459962 |s2cid=85016468 |issn=0003-0147}}</ref>
At un-invadable steady state <math>\partial Q/\partial N_i = 0</math> for all surviving species <math>i</math> and <math>\partial Q/\partial N_i > 0</math> for all extinct species <math>i</math>.<ref>{{Cite journal |last=Arthur |first=Robert Mac |title=Species Packing, and What Competition Minimizes |date=December 1969 |journal=Proceedings of the National Academy of Sciences |language=en |volume=64 |issue=4 |pages=1369–1371 |doi=10.1073/pnas.64.4.1369 |issn=0027-8424 |pmc=223294 |pmid=16591810 |doi-access=free }}</ref><ref>{{Cite journal |last=MacArthur |first=Robert |date=1970-05-01 |title=Species packing and competitive equilibrium for many species |url=https://dx.doi.org/10.1016/0040-5809%2870%2990039-0 |journal=Theoretical Population Biology |volume=1 |issue=1 |pages=1–11 |doi=10.1016/0040-5809(70)90039-0 |pmid=5527624 |issn=0040-5809}}</ref>


=== Minimum Environmental Perturbation Principle (MEPP) ===
=== Minimum Environmental Perturbation Principle (MEPP) ===
MacArthur's Minimization Principle has been extended to the more general Minimum Environmental Perturbation Principle (MEPP) which maps certain niche CRM models to constrained optimization problems. When the population growth conferred upon a species by consuming a resource is related to the impact the species' consumption has on the resource's abundance through the equation,<math display="block">q_{i\alpha}(\mathbf R) = - a_i(\mathbf R)b_\alpha(\mathbf R) \frac{\partial g_i}{\partial R_\alpha}
MacArthur's Minimization Principle has been extended to the more general Minimum Environmental Perturbation Principle (MEPP) which maps certain niche CRM models to constrained optimization problems. When the population growth conferred upon a species by consuming a resource is related to the impact the species' consumption has on the resource's abundance through the equation,<math display="block">q_{i\alpha}(\mathbf R) = - a_i(\mathbf R)b_\alpha(\mathbf R) \frac{\partial g_i}{\partial R_\alpha}
,</math>species-resource interactions are said to be ''symmetric''. In the above equation <math>a_i</math> and <math>b_\alpha</math> are arbitrary functions of resource abundances. When this symmetry condition is satisfied, it can be shown that there exists a function <math>d(\mathbf R)</math> such that:<ref name=":0" /><math display="block">\frac{\partial d}{\partial R_\alpha}
,</math> species-resource interactions are said to be ''symmetric''. In the above equation <math>a_i</math> and <math>b_\alpha</math> are arbitrary functions of resource abundances. When this symmetry condition is satisfied, it can be shown that there exists a function <math>d(\mathbf R)</math> such that:<ref name=":0" /><math display="block">\frac{\partial d}{\partial R_\alpha}
=
=
-\frac{h_\alpha(\mathbf R)}{b_\alpha (\mathbf R)}.</math>After determining this function <math>d</math>, the steady-state uninvadable resource abundances and species populations are the solution to the [[constrained optimization problem]]:<math display="block">\begin{align}
-\frac{h_\alpha(\mathbf R)}{b_\alpha (\mathbf R)}.</math>After determining this function <math>d</math>, the steady-state uninvadable resource abundances and species populations are the solution to the [[constrained optimization problem]]:<math display="block">\begin{align}
Line 233: Line 152:
0 &\leq N_i ,&& \qquad i =1,\dots,S.
0 &\leq N_i ,&& \qquad i =1,\dots,S.
\end{align}</math>Lines 1, 3, and 4 are the statements of feasibility and uninvadability: if <math>\overline N_i > 0</math>, then <math>g_i(\mathbf R)</math> must be zero otherwise the system would not be at steady state, and if <math>\overline N_i = 0 </math>, then <math>g_i(\mathbf R)</math> must be non-positive otherwise species <math>i </math> would be able to invade. Line 2 is the stationarity condition and the steady-state condition for the resources in nice CRMs. The function <math>d(\mathbf R) </math> can be interpreted as a distance by defining the point in the state space of resource abundances at which it is zero, <math>\mathbf R_0 </math>, to be its minimum. The [[Lagrange multiplier|Lagrangian]] for the [[Duality (optimization)|dual problem]] which leads to the above KKT conditions is,<math display="block">L(\mathbf R,\{N_i\}) =
\end{align}</math>Lines 1, 3, and 4 are the statements of feasibility and uninvadability: if <math>\overline N_i > 0</math>, then <math>g_i(\mathbf R)</math> must be zero otherwise the system would not be at steady state, and if <math>\overline N_i = 0 </math>, then <math>g_i(\mathbf R)</math> must be non-positive otherwise species <math>i </math> would be able to invade. Line 2 is the stationarity condition and the steady-state condition for the resources in nice CRMs. The function <math>d(\mathbf R) </math> can be interpreted as a distance by defining the point in the state space of resource abundances at which it is zero, <math>\mathbf R_0 </math>, to be its minimum. The [[Lagrange multiplier|Lagrangian]] for the [[Duality (optimization)|dual problem]] which leads to the above KKT conditions is,<math display="block">L(\mathbf R,\{N_i\}) =
d(\mathbf R) - \sum_{i = 1}^S N_i g_i(\mathbf R). </math>
d(\mathbf R) - \sum_{i = 1}^S N_i g_i(\mathbf R). </math> In this picture, the unconstrained value of <math>\mathbf R</math> that minimizes <math>d(\mathbf R)</math> (i.e., the steady-state resource abundances in the absence of any consumers) is known as the resource supply vector.

== Geometric perspectives ==
== Geometric perspectives ==
The steady states of consumer resource models can be analyzed using geometric means in the space of resource abundances.<ref name=":4">{{Cite journal |last1=Mancuso |first1=Christopher P |last2=Lee |first2=Hyunseok |last3=Abreu |first3=Clare I |last4=Gore |first4=Jeff |last5=Khalil |first5=Ahmad S |date=2021-09-03 |editor-last=Shou |editor-first=Wenying |editor2-last=Walczak |editor2-first=Aleksandra M |editor3-last=Shou |editor3-first=Wenying |title=Environmental fluctuations reshape an unexpected diversity-disturbance relationship in a microbial community |journal=eLife |volume=10 |pages=e67175 |doi=10.7554/eLife.67175 |issn=2050-084X |pmc=8460265 |pmid=34477107 |doi-access=free}}</ref><ref>{{Cite journal |last1=Tikhonov |first1=Mikhail |last2=Monasson |first2=Remi |date=2017-01-27 |title=Collective Phase in Resource Competition in a Highly Diverse Ecosystem |url=https://link.aps.org/doi/10.1103/PhysRevLett.118.048103 |journal=Physical Review Letters |volume=118 |issue=4 |pages=048103 |doi=10.1103/PhysRevLett.118.048103|pmid=28186794 |arxiv=1609.01270 |bibcode=2017PhRvL.118d8103T }}</ref><ref name=":5" />


=== Zero Net-Growth Isoclines (ZNGIs) ===
=== Zero net-growth isoclines (ZNGIs) ===
For a community to satisfy the uninvisibility and steady-state conditions, the steady-state resource abundances (denoted <math>
[[Nullcline|Nullclines]] [[Isocline]] [[File:Coexistence cones.png|thumb|Demonstration of ZNGIs and coexistence cones in the MacArthur consumer-resource model with <math>M = 2</math> resources and <math>S = 5</math> species. Left: ZNGIs for each species and the infeasible region. Center: Coexistence cones and resource dynamics for various carrying capacities (A, ..., E). Right: Species population dynamics corresponding to each carrying capacity in the center pane.]]
\mathbf R^\star

</math>) must satisfy,
Talk about ZNGIs<ref name=":4" />
<math display="block">
g_i(\mathbf R^\star) \leq 0,
</math>
for all species <math>
i
</math>. The inequality is saturated if and only if species <math>
i
</math> survives. Each of these conditions specifies a region in the space of possible steady-state resource abundances, and the realized steady-state resource abundance is restricted to the intersection of these regions. The boundaries of these regions, specified by <math>
g_i(\mathbf R^\star) = 0
</math>, are known as the zero net-growth [[Isocline|isoclines]] (ZNGIs). If species <math>
i = 1,\dots,S^\star
</math>survive, then the steady-state resource abundances must satisfy, <math>
g_1(\mathbf R^\star),\ldots, g_{S^\star}(\mathbf R^\star) = 0
</math>. The structure and locations of the intersections of the ZNGIs thus determine what species and feasibly coexist; the realized steady-state community is dependent on the supply of resources and can be analyzed by examining coexistence cones.


=== Coexistence cones ===
=== Coexistence cones ===
The structure of ZNGI intersections determines what species can feasibly coexist but does not determine what set of coexisting species will be realized. Coexistence cones determine what species determine what species will survive in an ecosystem given a resource supply vector. A coexistence cone generated by a set of species <math> i = 1,\ldots, S^\star </math> is defined to be the set of possible resource supply vectors which will lead to a community containing precisely the species <math> i =1,\ldots,S^\star </math>.
Talk about [[Coexistence theory|coexistence]] cones

To see the cone structure, consider that in the MacArthur or Tilman models, the steady-state non-depleted resource abundances must satisfy,<math display="block"> \mathbf K = \mathbf R^\star + \sum_{i=1}^S N_i \mathbf C_i,</math> where <math> \mathbf K</math> is a vector containing the carrying capacities/supply rates, and <math> \mathbf C_i = (c_{i1},\ldots,c_{iM^\star})</math> is the <math> i</math>th row of the consumption matrix <math> c_{i\alpha
}</math>, considered as a vector. As the surviving species are exactly those with positive abundances, the sum term becomes a sum only over surviving species, and the right-hand side resembles the expression for a [[convex cone]] with apex <math> \mathbf R^\star</math> and whose generating vectors are the <math> \mathbf C_i</math> for the surviving species <math> i</math>.

== Complex ecosystems ==
== Complex ecosystems ==
In an ecosystem with many species and resources, the behavior of consumer-resource models can be analyzed using tools from statistical physics, particularly [[mean-field theory]] and the [[cavity method]].<ref name=":7">{{Cite journal |last1=Advani |first1=Madhu |last2=Bunin |first2=Guy |last3=Mehta |first3=Pankaj |date=2018-03-20 |title=Statistical physics of community ecology: a cavity solution to MacArthur's consumer resource model |journal=Journal of Statistical Mechanics: Theory and Experiment |volume=2018 |issue=3 |pages=033406 |doi=10.1088/1742-5468/aab04e |issn=1742-5468 |pmc=6329381 |pmid=30636966|bibcode=2018JSMTE..03.3406A }}</ref><ref>{{Cite journal |last1=Marsland |first1=Robert |last2=Cui |first2=Wenping |last3=Mehta |first3=Pankaj |date=2020-02-24 |title=A minimal model for microbial biodiversity can reproduce experimentally observed ecological patterns |journal=Scientific Reports |language=en |volume=10 |issue=1 |pages=3308 |doi=10.1038/s41598-020-60130-2 |pmid=32094388 |pmc=7039880 |arxiv=1904.12914 |bibcode=2020NatSR..10.3308M |issn=2045-2322}}</ref><ref name=":9">{{Citation |last1=Mehta |first1=Pankaj |title=Cross-feeding shapes both competition and cooperation in microbial ecosystems |date=2021-10-10 |arxiv=2110.04965 |last2=Marsland III |first2=Robert}}</ref> In the large ecosystem limit, there is an explosion of the number of parameters. For example, in the MacArthur model, <math>O(SM)</math> parameters are needed. In this limit, parameters may be considered to be drawn from some distribution which leads to a distribution of steady-state abundances. These distributions of steady-state abundances can then be determined by deriving mean-field equations for random variables representing the steady-state abundances of a randomly selected species and resource.


=== MacArthur consumer resource model cavity solution ===
== Steady states and attractors ==
In the MCRM, the model parameters can be taken to be random variables with means and variances:<math>\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.</math>

With this parameterization, in the thermodynamic limit (i.e., <math>M,S \to \infty </math> with <math>S/M = \Theta(1)</math>), the steady-state resource and species abundances are modeled as a random variable, <math>N, R</math>, which satisfy the self-consistent mean-field equations,<ref name=":7" /><math display="block">\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}</math> where <math>\langle N\rangle, \langle N^2\rangle, \langle R\rangle, \rangle R^2\rangle</math> are all moments which are determined self-consistently, <math>Z_R,Z_N</math> are independent [[standard normal]] random variables, and <math>\nu = \langle \partial N/\partial m \rangle</math> and <math>\chi = \langle \partial R/\partial K \rangle</math> are average susceptibilities which are also determined self-consistently.

This mean-field framework can determine the moments and exact form of the abundance distribution, the average susceptibilities, and the fraction of species and resources that survive at a steady state.

Similar mean-field analyses have been performed for the externally supplied resources model,<ref name=":1" /> the Tilman model,<ref name=":2" /> and the microbial consumer-resource model.<ref name=":9" /> These techniques were first developed to analyze the [[random generalized Lotka–Volterra model]].

== See also ==

* [[Theoretical ecology]]
* [[Community (ecology)]]
* [[Competition (biology)]]
* [[Lotka–Volterra equations]]
* [[Competitive Lotka–Volterra equations]]
* [[Generalized Lotka–Volterra equation]]
* [[Random generalized Lotka–Volterra model]]


== References ==
== References ==
<!-- Inline citations added to your article will automatically display here. See en.wikipedia.org/wiki/WP:REFB for instructions on how to add citations. -->
{{reflist}}
{{reflist}}

== Further reading ==

* {{cite arXiv |eprint=2403.05497 |last1=Cui |first1=Wenping |author2=Robert Marsland III |last3=Mehta |first3=Pankaj |title=Les Houches Lectures on Community Ecology: From Niche Theory to Statistical Mechanics |date=2024 |class=q-bio.PE }}
* Stefano Allesina's Community Ecology course lecture notes: https://stefanoallesina.github.io/Theoretical_Community_Ecology/

{{modelling ecosystems|expanded=other}}

[[Category:Ecology]]
[[Category:Ordinary differential equations]]
[[Category:Mathematical modeling]]
[[Category:Biophysics]]
[[Category:Community ecology]]
[[Category:Ecological niche]]
[[Category:Population ecology]]
[[Category:Dynamical systems]]
[[Category:Random dynamical systems]]
[[Category:Theoretical ecology]]

Latest revision as of 04:28, 30 June 2024

In theoretical ecology and nonlinear dynamics, consumer-resource models (CRMs) are a class of ecological models in which a community of consumer species compete for a common pool of resources. Instead of species interacting directly, all species-species interactions are mediated through resource dynamics. Consumer-resource models have served as fundamental tools in the quantitative development of theories of niche construction, coexistence, and biological diversity. These models can be interpreted as a quantitative description of a single trophic level.[1][2]


A general consumer-resource model consists of M resources whose abundances are and S consumer species whose populations are . A general consumer-resource model is described by the system of coupled ordinary differential equations, where , depending only on resource abundances, is the per-capita growth rate of species , and is the growth rate of resource . An essential feature of CRMs is that species growth rates and populations are mediated through resources and there are no explicit species-species interactions. Through resource interactions, there are emergent inter-species interactions.

Originally introduced by Robert H. MacArthur[3] and Richard Levins,[4] consumer-resource models have found success in formalizing ecological principles and modeling experiments involving microbial ecosystems.[5][6]

Models

[edit]

Niche models

[edit]

Niche models are a notable class of CRMs which are described by the system of coupled ordinary differential equations,[7][8]

where is a vector abbreviation for resource abundances, is the per-capita growth rate of species , is the growth rate of species in the absence of consumption, and is the rate per unit species population that species depletes the abundance of resource through consumption. In this class of CRMs, consumer species' impacts on resources are not explicitly coordinated; however, there are implicit interactions.

MacArthur consumer-resource model (MCRM)

[edit]

The MacArthur consumer-resource model (MCRM), named after Robert H. MacArthur, is a foundational CRM for the development of niche and coexistence theories.[9][10] The MCRM is given by the following set of coupled ordinary differential equations:[11][12][8]where is the relative preference of species for resource and also the relative amount by which resource is depleted by the consumption of consumer species ; is the steady-state carrying capacity of resource in absence of consumption (i.e., when is zero); and are time-scales for species and resource dynamics, respectively; is the quality of resource ; and is the natural mortality rate of species . This model is said to have self-replenishing resource dynamics because when , each resource exhibits independent logistic growth. Given positive parameters and initial conditions, this model approaches a unique uninvadable steady state (i.e., a steady state in which the re-introduction of a species which has been driven to extinction or a resource which has been depleted leads to the re-introduced species or resource dying out again).[7] Steady states of the MCRM satisfy the competitive exclusion principle: the number of coexisting species is less than or equal to the number of non-depleted resources. In other words, the number of simultaneously occupiable ecological niches is equal to the number of non-depleted resources.

Externally supplied resources model

[edit]

The externally supplied resource model is similar to the MCRM except the resources are provided at a constant rate from an external source instead of being self-replenished. This model is also sometimes called the linear resource dynamics model. It is described by the following set of coupled ordinary differential equations:[11][8]where all the parameters shared with the MCRM are the same, and is the rate at which resource is supplied to the ecosystem. In the eCRM, in the absence of consumption, decays to exponentially with timescale . This model is also known as a chemostat model.

Tilman consumer-resource model (TCRM)

[edit]

The Tilman consumer-resource model (TCRM), named after G. David Tilman, is similar to the externally supplied resources model except the rate at which a species depletes a resource is no longer proportional to the present abundance of the resource. The TCRM is the foundational model for Tilman's R* rule. It is described by the following set of coupled ordinary differential equations:[8]where all parameters are shared with the MCRM. In the TCRM, resource abundances can become nonphysically negative.[7]

Microbial consumer-resource model (MiCRM)

[edit]

The microbial consumer resource model describes a microbial ecosystem with externally supplied resources where consumption can produce metabolic byproducts, leading to potential cross-feeding. It is described by the following set of coupled ODEs:where all parameters shared with the MCRM have similar interpretations; is the fraction of the byproducts due to consumption of resource which are converted to resource and is the "leakage fraction" of resource governing how much of the resource is released into the environment as metabolic byproducts.[13][8]

Symmetric interactions and optimization

[edit]

MacArthur's Minimization Principle

[edit]

For the MacArthur consumer resource model (MCRM), MacArthur introduced an optimization principle to identify the uninvadable steady state of the model (i.e., the steady state so that if any species with zero population is re-introduced, it will fail to invade, meaning the ecosystem will return to said steady state). To derive the optimization principle, one assumes resource dynamics become sufficiently fast (i.e., ) that they become entrained to species dynamics and are constantly at steady state (i.e., ) so that is expressed as a function of . With this assumption, one can express species dynamics as, where denotes a sum over resource abundances which satisfy . The above expression can be written as , where,

At un-invadable steady state for all surviving species and for all extinct species .[14][15]

Minimum Environmental Perturbation Principle (MEPP)

[edit]

MacArthur's Minimization Principle has been extended to the more general Minimum Environmental Perturbation Principle (MEPP) which maps certain niche CRM models to constrained optimization problems. When the population growth conferred upon a species by consuming a resource is related to the impact the species' consumption has on the resource's abundance through the equation, species-resource interactions are said to be symmetric. In the above equation and are arbitrary functions of resource abundances. When this symmetry condition is satisfied, it can be shown that there exists a function such that:[7]After determining this function , the steady-state uninvadable resource abundances and species populations are the solution to the constrained optimization problem:The species populations are the Lagrange multipliers for the constraints on the second line. This can be seen by looking at the KKT conditions, taking to be the Lagrange multipliers:Lines 1, 3, and 4 are the statements of feasibility and uninvadability: if , then must be zero otherwise the system would not be at steady state, and if , then must be non-positive otherwise species would be able to invade. Line 2 is the stationarity condition and the steady-state condition for the resources in nice CRMs. The function can be interpreted as a distance by defining the point in the state space of resource abundances at which it is zero, , to be its minimum. The Lagrangian for the dual problem which leads to the above KKT conditions is, In this picture, the unconstrained value of that minimizes (i.e., the steady-state resource abundances in the absence of any consumers) is known as the resource supply vector.

Geometric perspectives

[edit]

The steady states of consumer resource models can be analyzed using geometric means in the space of resource abundances.[16][17][8]

Zero net-growth isoclines (ZNGIs)

[edit]

For a community to satisfy the uninvisibility and steady-state conditions, the steady-state resource abundances (denoted ) must satisfy, for all species . The inequality is saturated if and only if species survives. Each of these conditions specifies a region in the space of possible steady-state resource abundances, and the realized steady-state resource abundance is restricted to the intersection of these regions. The boundaries of these regions, specified by , are known as the zero net-growth isoclines (ZNGIs). If species survive, then the steady-state resource abundances must satisfy, . The structure and locations of the intersections of the ZNGIs thus determine what species and feasibly coexist; the realized steady-state community is dependent on the supply of resources and can be analyzed by examining coexistence cones.

Coexistence cones

[edit]

The structure of ZNGI intersections determines what species can feasibly coexist but does not determine what set of coexisting species will be realized. Coexistence cones determine what species determine what species will survive in an ecosystem given a resource supply vector. A coexistence cone generated by a set of species is defined to be the set of possible resource supply vectors which will lead to a community containing precisely the species .

To see the cone structure, consider that in the MacArthur or Tilman models, the steady-state non-depleted resource abundances must satisfy, where is a vector containing the carrying capacities/supply rates, and is the th row of the consumption matrix , considered as a vector. As the surviving species are exactly those with positive abundances, the sum term becomes a sum only over surviving species, and the right-hand side resembles the expression for a convex cone with apex and whose generating vectors are the for the surviving species .

Complex ecosystems

[edit]

In an ecosystem with many species and resources, the behavior of consumer-resource models can be analyzed using tools from statistical physics, particularly mean-field theory and the cavity method.[18][19][20] In the large ecosystem limit, there is an explosion of the number of parameters. For example, in the MacArthur model, parameters are needed. In this limit, parameters may be considered to be drawn from some distribution which leads to a distribution of steady-state abundances. These distributions of steady-state abundances can then be determined by deriving mean-field equations for random variables representing the steady-state abundances of a randomly selected species and resource.

MacArthur consumer resource model cavity solution

[edit]

In the MCRM, the model parameters can be taken to be random variables with means and variances:

With this parameterization, in the thermodynamic limit (i.e., with ), the steady-state resource and species abundances are modeled as a random variable, , which satisfy the self-consistent mean-field equations,[18] where are all moments which are determined self-consistently, are independent standard normal random variables, and and are average susceptibilities which are also determined self-consistently.

This mean-field framework can determine the moments and exact form of the abundance distribution, the average susceptibilities, and the fraction of species and resources that survive at a steady state.

Similar mean-field analyses have been performed for the externally supplied resources model,[11] the Tilman model,[12] and the microbial consumer-resource model.[20] These techniques were first developed to analyze the random generalized Lotka–Volterra model.

See also

[edit]

References

[edit]
  1. ^ Chase, Jonathan M.; Leibold, Mathew A. (2003). Ecological Niches. University of Chicago Press. doi:10.7208/chicago/9780226101811.001.0001. ISBN 978-0-226-10180-4.
  2. ^ Pimm, Stuart L. (September 1983). "TILMAN, D. 1982. Resource competition and community structure. Monogr. Pop. Biol. 17. Princeton University Press, Princeton, N.J. 296 p. $27.50". Limnology and Oceanography. 28 (5): 1043–1045. Bibcode:1983LimOc..28.1043P. doi:10.4319/lo.1983.28.5.1043.
  3. ^ MacArthur, Robert (1970-05-01). "Species packing and competitive equilibrium for many species". Theoretical Population Biology. 1 (1): 1–11. doi:10.1016/0040-5809(70)90039-0. ISSN 0040-5809. PMID 5527624.
  4. ^ Levins, Richard (1968). Evolution in Changing Environments: Some Theoretical Explorations. (MPB-2). Princeton University Press. doi:10.2307/j.ctvx5wbbh. ISBN 978-0-691-07959-2. JSTOR j.ctvx5wbbh.
  5. ^ Goldford, Joshua E.; Lu, Nanxi; Bajić, Djordje; Estrela, Sylvie; Tikhonov, Mikhail; Sanchez-Gorostiaga, Alicia; Segrè, Daniel; Mehta, Pankaj; Sanchez, Alvaro (2018-08-03). "Emergent simplicity in microbial community assembly". Science. 361 (6401): 469–474. Bibcode:2018Sci...361..469G. doi:10.1126/science.aat1168. ISSN 0036-8075. PMC 6405290. PMID 30072533.
  6. ^ Dal Bello, Martina; Lee, Hyunseok; Goyal, Akshit; Gore, Jeff (October 2021). "Resource–diversity relationships in bacterial communities reflect the network structure of microbial metabolism". 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.
  7. ^ a b c d Marsland, Robert; Cui, Wenping; Mehta, Pankaj (2020-09-01). "The Minimum Environmental Perturbation Principle: A New Perspective on Niche Theory". The American Naturalist. 196 (3): 291–305. arXiv:1901.09673. doi:10.1086/710093. ISSN 0003-0147. PMID 32813998. S2CID 59316948.
  8. ^ a b c d e f Cui, Wenping; Marsland III, Robert; Mehta, Pankaj (2024-03-08). "Les Houches Lectures on Community Ecology: From Niche Theory to Statistical Mechanics". arXiv:2403.05497 [q-bio.PE].
  9. ^ "Resource Competition and Community Structure. (MPB-17), Volume 17 | Princeton University Press". press.princeton.edu. 1982-08-21. Retrieved 2024-03-18.
  10. ^ Chase, Jonathan M.; Leibold, Mathew A. Ecological Niches: Linking Classical and Contemporary Approaches. Interspecific Interactions. Chicago, IL: University of Chicago Press.
  11. ^ a b c Cui, Wenping; Marsland, Robert; Mehta, Pankaj (2020-07-21). "Effect of Resource Dynamics on Species Packing in Diverse Ecosystems". Physical Review Letters. 125 (4): 048101. arXiv:1911.02595. Bibcode:2020PhRvL.125d8101C. doi:10.1103/PhysRevLett.125.048101. PMC 8999492. PMID 32794828.
  12. ^ a b Cui, Wenping; Marsland, Robert; Mehta, Pankaj (2021-09-27). "Diverse communities behave like typical random ecosystems". Physical Review E. 104 (3): 034416. Bibcode:2021PhRvE.104c4416C. doi:10.1103/PhysRevE.104.034416. PMC 9005152. PMID 34654170.
  13. ^ Marsland, Robert; Cui, Wenping; Mehta, Pankaj (2020-02-24). "A minimal model for microbial biodiversity can reproduce experimentally observed ecological patterns". 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.
  14. ^ Arthur, Robert Mac (December 1969). "Species Packing, and What Competition Minimizes". Proceedings of the National Academy of Sciences. 64 (4): 1369–1371. doi:10.1073/pnas.64.4.1369. ISSN 0027-8424. PMC 223294. PMID 16591810.
  15. ^ MacArthur, Robert (1970-05-01). "Species packing and competitive equilibrium for many species". Theoretical Population Biology. 1 (1): 1–11. doi:10.1016/0040-5809(70)90039-0. ISSN 0040-5809. PMID 5527624.
  16. ^ Mancuso, Christopher P; Lee, Hyunseok; Abreu, Clare I; Gore, Jeff; Khalil, Ahmad S (2021-09-03). Shou, Wenying; Walczak, Aleksandra M; Shou, Wenying (eds.). "Environmental fluctuations reshape an unexpected diversity-disturbance relationship in a microbial community". eLife. 10: e67175. doi:10.7554/eLife.67175. ISSN 2050-084X. PMC 8460265. PMID 34477107.
  17. ^ Tikhonov, Mikhail; Monasson, Remi (2017-01-27). "Collective Phase in Resource Competition in a Highly Diverse Ecosystem". Physical Review Letters. 118 (4): 048103. arXiv:1609.01270. Bibcode:2017PhRvL.118d8103T. doi:10.1103/PhysRevLett.118.048103. PMID 28186794.
  18. ^ a b Advani, Madhu; Bunin, Guy; Mehta, Pankaj (2018-03-20). "Statistical physics of community ecology: a cavity solution to MacArthur's consumer resource model". 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.
  19. ^ Marsland, Robert; Cui, Wenping; Mehta, Pankaj (2020-02-24). "A minimal model for microbial biodiversity can reproduce experimentally observed ecological patterns". 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.
  20. ^ a b Mehta, Pankaj; Marsland III, Robert (2021-10-10), Cross-feeding shapes both competition and cooperation in microbial ecosystems, arXiv:2110.04965

Further reading

[edit]