Anhang A: Beweis von Theorem 4.1
Wenn die Anzahl der Arten streng größer ist als die Anzahl der Ressourcen, \(n>m\), ist das System
Es gibt eine nichttriviale Lösung \((c_1, \dots , c_n)\).
Theorem A.1
Angenommen, dass \(\lim _{\Vert \mathbf {x}\Vert \rightarrow \infty } R_j(\mathbf {x})=-\infty , j=1,\dots ,m\). Nehmen wir weiter an, dass n Spezies gemäß (4.3) interagieren, die Anzahl der Spezies größer ist als die Anzahl der Ressourcen \(n>m\) und die Ressourcen von den Speziesdichten gemäß (2.2) abhängen, so dass sie schließlich erschöpft werden. Nehmen wir weiter an, dass \(g_i(\mathbf {x})=1\) und
für \(j=1,\dots ,m\). Sei \((c_1,\dots ,c_n)\) eine nicht-triviale Lösung von (A.1) und angenommen, dass \(\sum _{i=1}^n c_i\left( \alpha _i+\frac{\sigma _{ii}}{2}\right) \ne 0\). Dann ist für beliebige Anfangsdichten \(\mathbf {x}(0)\in (0,\infty )^n\) mit Wahrscheinlichkeit 1
Beweis von Theorem 4.1
Setzen Sie \(g_i(\mathbf {x})=1\), \(i=1,\dots , n\) und \(\sum _jb_{ij}>b_m>0\) für jedes i und irgendein \(b_m>0\) voraus. Beachten Sie, dass, wenn \(\sum _j b_{ij}=0\), wir \(R_j\) aus der Gleichung entfernen können. Nehmen wir an, dass
Dann, da \(\lim _{\Vert \mathbf {x}\Vert \rightarrow \infty } R_j(\mathbf {x})=-\infty \), haben wir, wenn \(|\mathbf {x}|\) groß ist, dass:
Was zusammen mit der Linearität des Diffusionsteils impliziert, dass Annahme 1.1 aus der Arbeit von HN16 mit \(\mathbf {c}=(1,\dots , 1)\) gilt. Daraus folgt, dass für jeden Startpunkt \(\mathbf {x}(0)\in (0,\infty )^n\) die SDE (4.3) eine eindeutige positive Lösung hat und nach Hening und Nguyen (2018a) (Gleichung (5.22)) mit Wahrscheinlichkeit 1
Indem wir eventuell alle \(c_i\) durch \(-c_i\) ersetzen, können wir annehmen, dass \(\sum _{i=1}^n c_i\left( \alpha _i+\frac{\sigma _{ii}}{2}\right) >0\). Wenn wir dies in Verbindung mit (4.3), (A.1) und dem Itô-Lemma verwenden, sehen wir, dass
Lässt man \(t\rightarrow \infty \) und verwendet, dass \(\lim _{t\rightarrow \infty }\dfrac{\sum _{i=1}^n c_iE_i(t)}{t}=0\) mit Wahrscheinlichkeit 1 ist, erhält man, dass mit Wahrscheinlichkeit 1
In Anbetracht von (A.2) impliziert dies, dass mit der Wahrscheinlichkeit 1
(\quadratisch \)
Anhang B: Beweis von Theorem 4.2
Angenommen, zwei Spezies interagieren gemäß
Die Ressource R hängt linear von den Speziesdichten ab
und \(b_i{\overline{R}}>\alpha _i, i=1,2\). Dann gibt es \(\beta _1, \beta _2>0\), so dass die beiden Arten koexistieren.
Beweis
Betrachten
Wenn die Spezies \(x_2\) abwesend ist, hat die Spezies \(x_1\) die ein-dimensionale Dynamik
Da \(b_1{\overline{R}}>\alpha _1\), können wir mit Hilfe von Hening und Nguyen (2018a) zeigen, dass der Prozess x(t) ein eindeutiges invariantes Maß auf \((0,\infty )\) hat, nämlich \(\mu _1\). Außerdem zeigt (Hening und Nguyen 2018a, Lemma 2.1 ), dass
oder
Die Invasionsrate von \(x_2\) in Bezug auf \(x_1\) kann durch (4.2) berechnet werden als
Ahnlich kann man die Invasionsrate von \(x_1\) bezüglich \(x_2\) berechnen als
Da \(b_i{\overline{R}}-\alpha _i>0\), kann man leicht sehen, dass \( \Lambda _{x_2}>0\), und \(\Lambda _{x_1}>0\) wenn \( \beta _1>\dfrac{b_2a_1(b_1{\overline{R}}-\alpha _1)}{b_2\overline{R}-\alpha _2}-b_1a_1. \) und \( \beta _2>\dfrac{b_1a_2(b_2\overline{R}-\alpha _2)}{b_1{\overline{R}}-\alpha _1}-b_2a_2. \) Wenn beide Invasionsraten positiv sind, erhalten wir nach Hening und Nguyen (2018a), dass die Arten koexistieren. \(\quadratisch \)
Anhang C: Beweis von Theorem 4.3
Wir konstruieren ein SDE-Beispiel für zwei Arten, die um eine abiotische Ressource konkurrieren und koexistieren. Wir bemerken, dass dies allein aufgrund des Terms der zufälligen zeitlichen Umweltvariation geschieht.
Theorem C.1
Angenommen, die Dynamik der beiden Arten ist gegeben durch
wobei f eine stetig differenzierbare Lipschitz-Funktion ist, die \(\lim _{x\rightarrow -\infty }f(x)=-\infty \), \(\dfrac{df(x)}{dx}>0, \dfrac{d^2f(x)}{dx^2}\le 0\) für alle \(x\in \mathbb {R}\) und \(\dfrac{d^2f(x)}{dx^2}<0\) für x in einem Teilintervall von \(\lefty ,\frac{{\overline{R}}}{a_1}\right) \). Sei \(a_1, a_2, \sigma _1, \alpha _1,\sigma _1, {\overline{R}}\) beliebige feste positive Konstanten, die \(f({\overline{R}})>\alpha _1+\dfrac{\sigma _1^2}{2}\) erfüllen. Dann gibt es ein Intervall \((c_0,c_1)\Teilmenge (0,\infty )\), so dass die beiden Arten für alle \(\alpha _2\in (c_0,c_1)\) koexistieren.
Beweis
Die Dynamik der Spezies \(x_1\) in Abwesenheit der Spezies \(x_2\) ist gegeben durch die eindimensionale SDE
Da \(\lim _{x\rightarrow \infty }f({\overline{R}}-a_1x)=-\infty \), und \(f(\overline{R})>\alpha _1+\dfrac{\sigma _1^2}{2}\), hat diese Diffusion ein eindeutiges invariantes Wahrscheinlichkeitsmaß \(\mu \) auf \((0,\infty )\), dessen Dichte streng positiv auf \((0,\infty )\) ist (siehe Borodin und Salminen (2016) oder Mao (1997)). Außerdem stellt man fest, dass \(\lim _{t\rightarrow \infty }\frac{\ln x(t)}{t}=0\) mit der Wahrscheinlichkeit 1 ist (unter Verwendung von Lemma 5.1 von Hening und Nguyen (2018a)) und unter Verwendung der Itô-Formel sieht man, dass
Da f eine konkave Funktion ist und \(\dfrac{d^2f(x)}{dx^2}>0\) für alle x in einem Teilintervall von \(\left( -\infty ,\frac{\overline{R}}{a_1}\right) \) müssen wir durch die Jensensche Ungleichung haben, dass
Die Tatsache, dass die Funktion f streng steigend ist, zwingt zusammen mit (C.2) und (C.3)
wobei \(f^{-1}\) der Kehrwert von f ist – er existiert, weil f streng steigend ist. Infolgedessen kann die Invasionsrate der Spezies \(x_2\) bezüglich \(x_1\) des invarianten Wahrscheinlichkeitsmaßes \(\mu \) mit Hilfe von (4.2) berechnet werden als
Dies impliziert, dass \(\Lambda _{x_2}>0\) wenn und nur wenn
Die Dynamik der Spezies \(x_2\) in Abwesenheit der Spezies \(x_1\) ist
Die positiven Lösungen dieser Gleichung konvergieren zum Punkt \(y^*=\dfrac{\overline{R}}-\alpha _2}{a_2}\) wenn und nur wenn
Die Invasionsrate von \(x_1\) in Bezug auf \(x_2\) ist
Da die Funktion f steigend ist, erhalten wir \(\Lambda _{x_1}>0\) wenn und nur wenn
Beachten Sie, dass \(f^{-1}\left( \alpha _1+\dfrac{\sigma _1^2}{2}\right) <{\overline{R}}\), da unter der Annahme \(f({\overline{R}})>\alpha _1+\dfrac{\sigma _1^2}{2}\). Als Ergebnis erhalten wir unter Verwendung der Ungleichungen (C.5), (C.6) und (C.7), dass \(\Lambda _{x_2}>0, \Lambda _{x_1}>0\) wenn überhaupt nur dann, wenn
Dies impliziert durch Theorem 3.1 oder durch Benaim (2018), dass die beiden Arten koexistieren. \(\quadratisch \)
Anhang D: Beweis von Theorem 5.1
Vermuten Sie, dass
Angenommen, es existiert ein Vektor \((c_1,\dots ,c_n)\), der gleichzeitig eine Lösung des Systems (5.5) für alle \(u\ in \{1,\dots ,N\}\) ist. Dann, mit Wahrscheinlichkeit 1,
außer möglicherweise für den kritischen Fall, wenn
wobei \((\nu _k)_{k\in \mathcal {N}}) das invariante Wahrscheinlichkeitsmaß der Markov-Kette (r(t)) ist.
Beweis
Unter der Bedingung (D.1) gibt es ein \(M>0\), so dass die Menge \(K_M:=\{\mathbf {x}\in \mathbb {R}^n: \Vert \mathbf {x}\Vert \le M\}\) ein globaler Attraktor von (5.4) ist. Daraus folgt, dass die Lösung von (5.4) schließlich in die kompakte Menge \(K_M\) eintritt und diese nie verlässt. Dies zeigt insbesondere, dass der Prozess \(\mathbf {x}(t)\) beschränkt ist. Als Nächstes beachten wir, dass wir annehmen können, dass
Andererseits, wenn \(\sum _{i=1}^n c_i\sum _{k=1}^N\alpha _i(k)\pi _k<0\), können wir \(c_i\) durch \(-c_i\) ersetzen, \(i=1,\dots , n\) und erhalten dann (D.3). Unter Verwendung von (5.1) und der Tatsache, dass \(c_i\) gleichzeitig (5.5) löst, erhalten wir
Aus der Ergodizität der Markov-Kette (r(t)) ergibt sich, dass mit der Wahrscheinlichkeit 1
Da \(\mathbf {x}(t)\) begrenzt ist, impliziert dies, dass mit Wahrscheinlichkeit 1
(\quadratisch \)
Anhang E: Beweis von Theorem 5.2
Nach Benaïm und Lobry (2016), Malrieu und Zitt (2017) und Malrieu und Phu (2016) genügt es, ein Beispiel zu finden, für das die Invasionsraten \(\Lambda _{x_1}, \Lambda _{x_2}\) positiv sind. Wir werden Benaïm und Lobry (2016) folgen, um die Invasionsraten der beiden Arten zu berechnen. Setzen Sie für \(u=1,2\)\(\mu _u=-\alpha _1(u)+b_1(u){\overline{R}}, \nu _u=-\alpha _2(u)+b_2(u){\overline{R}}) \({\overline{a}}_u=\dfrac{b_1(u)a_1(u)}{\mu _u}, \overline{b}_u=\dfrac{b_1(u)a_2(u)}{\mu _u}, \overline{c}_u=\dfrac{b_2(u)a_1(u)}{\nu _u}, \overline{d}(u)=\dfrac{b_2(u)a_2(u)}{\nu _u}\), \(p_u=\dfrac{1}{{\overline{a}}_u}, q_u=\dfrac{1}{{\overline{d}}_u}, \gamma _1=\dfrac{q_{12}}{\mu _u}, \gamma _2=\dfrac{q_{21}}{\nu _u}). Wenn \(p_1\ne p_2\), nehmen wir ohne Verlust der Allgemeinheit an, dass \(p_1<p_2\). Definieren Sie die Funktionen
und
Nach Benaïm und Lobry (2016) haben wir
Der Ausdruck für \(\Lambda _{x_1}\) kann durch Vertauschen von \(\mu _i\) und \(\nu _i\) erhalten werden, \(({\overline{a}}_i,{\overline{c}}_i)\) mit \(({\overline{d}}_i,\overline{b}_i)\), und \(p_i\) mit \(q_i\).
Für das Beispiel aus Abb. 3 und 4 haben wir die Integralgleichung aus (E.1) zusammen mit dem numerischen Integrationspaket von Mathematica, um \(\Lambda _{x_1}>0\) und \(\Lambda _{x_2}>0\) zu finden. Dies impliziert durch Theorem 3.1, Benaïm und Lobry (2016) oder durch Benaim (2018), dass die beiden Arten koexistieren.