Asymptotic stability and diffusion-driven pattern formation in a predator-prey system with two chemicals

Titel: Asymptotische Stabiliteit en Diffusie-gedreven Patroonvorming in een Roofdier-Bytvogelsysteem met Twee Chemische Stoffen

Auteurs: Gnanasekaran Shanmugasundaram, Jitraj Saha, Oluwole Daniel Makinde, en Joydev Chattopadhyay.

1. Probleemstelling en Motivatie

Dit onderzoek richt zich op de analyse van een roofdier-bytvogelsysteem dat gekoppeld is aan twee chemische stoffen in een begrensde domein $\Omega \subset \mathbb{R}^n$ ($n \ge 2$). Traditionele ecologische modellen vertrouwen vaak op zelfdiffusie (self-diffusion), wat echter onvoldoende is om complexe ruimtelijke patronen te genereren.

De kern van dit probleem ligt in kruisdiffusie (cross-diffusion) die wordt bemiddeld door chemische signalen:

• Prooi ($u_2$) geeft een chemisch signaal ($v_1$) af dat roofdieren ($u_1$) aantrekt (adversieve beweging).
• Roofdieren ($u_1$) geven een chemisch signaal ($v_2$) af dat prooien afschrikt (repulsieve beweging).

Het doel is om de globale existentie van oplossingen te bewijzen, de asymptotische stabiliteit van evenwichtspunten te analyseren en de mechanismen te begrijpen die leiden tot Turing-patronen (ruimtelijke patronen ontstaan door diffusie-instabiliteit) in dit vier-vergelijkingssysteem.

2. Wiskundig Model

Het systeem wordt beschreven door vier gekoppelde parabolische partiële differentiaalvergelijkingen met homogene Neumann-randvoorwaarden:

$$\begin{cases} u_{1t} = \nabla \cdot (d_{11}\nabla u_1 + d_{12}\nabla v_1) + u_1(\mu_1 - \lambda_1 u_1 + \eta_1 u_2) \\ u_{2t} = \nabla \cdot (d_{21}\nabla u_2 - d_{22}\nabla v_2) + u_2(\mu_2 - \lambda_2 u_2 - \eta_2 u_1) \\ v_{1t} = d_3 \Delta v_1 + a_1 u_2 - b_1 v_1 \\ v_{2t} = d_4 \Delta v_2 + a_2 u_1 - b_2 v_2 \end{cases}$$

Waarbij:

• $u_1, u_2$: Dichtheid van roofdier en prooi.
• $v_1, v_2$: Concentratie van de chemische stoffen (respectievelijk geproduceerd door prooi en roofdier).
• $d_{ij}$: Diffusiecoëfficiënten (zelfdiffusie en kruisdiffusie).
• De voorwaarden $4d_{11}d_3 > d_{12}^2$ en $4d_{21}d_4 > d_{22}^2$ worden gesteld om de paraboliciteit van het systeem te garanderen.

3. Methodologie

De auteurs hanteren een combinatie van rigoureuze wiskundige analyse en numerieke simulatie:

1. Existentie en Begrensbaarheid:

   • Gebruikmakend van de theorie van quasi-lineaire parabolische vergelijkingen (Amann's theorie) wordt de lokale existentie van klassieke oplossingen bewezen.
   • Lyapunov-functies en Moser-Alikakos iteratie worden toegepast om de uniformiteit en globale begrensbaarheid ($L^\infty$-norm) van de oplossingen te bewijzen, wat leidt tot globale existentie voor alle $t > 0$.
   • De positiviteit van de oplossingen wordt bewezen via de methode van negatieve delen en Gronwall's ongelijkheid.

2. Asymptotische Stabiliteit:

   • Er worden geschikte Lyapunov-functionalen geconstrueerd om de globale asymptotische stabiliteit van de evenwichtspunten aan te tonen.
   • Twee scenario's worden onderzocht:
     • Co-existentie: Wanneer de predatie-intensiteit $\eta_2$ laag is ($\eta_2 < \frac{\mu_2 \lambda_1}{\mu_1}$), convergeert het systeem exponentieel naar een uniek positief evenwicht.
     • Uitsterven van prooi: Wanneer $\eta_2$ hoog is ($\eta_2 \ge \frac{\mu_2 \lambda_1}{\mu_1}$), convergeert het systeem naar een semi-triviaal evenwicht waarbij de prooi uitsterft.

3. Bifurcatie-analyse (Turing-instabiliteit):

   • Een lineaire stabiliteitsanalyse wordt uitgevoerd rond het homogene evenwichtspunt.
   • Vanwege de complexiteit van de vierde-orde karakteristieke vergelijking (die ontstaat door de koppeling van vier vergelijkingen), worden analytische voorwaarden voor Turing-bifurcatie als te complex beschouwd.
   • In plaats daarvan worden numerieke simulaties (met de Finite Difference Method) gebruikt om de instabiliteitsgebieden te identificeren en de criteria voor Turing-bifurcatie af te leiden.

4. Belangrijkste Resultaten

A. Wiskundige Stabiliteit

• Het systeem bezit een unieke klassieke oplossing die uniform begrensd blijft.
• Co-existentie: Als de kruisdiffusie-coëfficiënten ($d_{12}, d_{22}$) klein genoeg zijn en de predatie-intensiteit onder een bepaalde drempel ligt, is het systeem globaal asymptotisch stabiel rond het co-existentie-evenwicht.
• Extinctie: Als de predatie-intensiteit te hoog is, wordt het systeem stabiel rond een evenwicht waarbij de prooi-populatie uitsterft, ongeacht de initiële condities (binnen een bepaald bereik).

B. Ruimtelijke Patroonvorming

De numerieke simulaties onthullen een rijk scala aan Turing-patronen die worden gedreven door kruisdiffusie:

• Invloed van Predatie-intensiteit ($\eta_2$):
  • Bij lage $\eta_2$ ontstaan uitgebreide labyrint-achtige strepen.
  • Bij toenemende $\eta_2$ gaat het systeem over naar gemengde streep-vlek patronen.
  • Bij hoge $\eta_2$ ontstaan geordende hexagonale vlekken (spots). Dit suggereert dat sterkere predatie leidt tot meer lokale aggregatie van roofdieren en een scherpe ruimtelijke segregatie.
• Invloed van Chemische Signalen:
  • De balans tussen de productie van chemische stoffen door roofdieren en prooi ($\eta_1, \eta_2$) bepaalt de geometrie van de patronen.
  • Een hoge $\eta_1$ (sterke reactie van roofdieren op prooi-chemie) bevordert geïsoleerde co-existentie-vlekken.
• Synchronisatie: Een opmerkelijke bevinding is dat roofdieren en prooi vaak maxima en minima op dezelfde ruimtelijke locaties vertonen. Dit komt doordat beide soorten reageren op dezelfde chemische gradiënten (prooi-chemie trekt roofdieren aan, roofdier-chemie drijft prooi weg), wat leidt tot gesynchroniseerde hotspots in plaats van de typische ruimtelijke verschuiving die men bij klassieke roofdier-prooi-systemen ziet.

5. Betekenis en Conclusie

Dit werk biedt nieuwe inzichten in de rol van chemisch bemiddelde kruisdiffusie in ecologische systemen:

1. Theoretische Bijdrage: Het bewijst de globale stabiliteit en existentie voor een complex vier-vergelijkingssysteem, wat zeldzaam is in de literatuur.
2. Mechanisme van Patroonvorming: Het toont aan dat kruisdiffusie, bemiddeld door zelf geproduceerde chemische stoffen, een krachtig mechanisme is voor het genereren van ruimtelijke heterogeniteit, zelfs wanneer zelfdiffusie alleen stabiel blijft.
3. Ecologische Implicaties: De resultaten suggereren dat de intensiteit van predatie en de gevoeligheid voor chemische signalen cruciale factoren zijn voor de vorming van ecologische niches. Het systeem kan schakelen tussen coherente territoriale structuren, gelokaliseerde hulpbronnenpatches en ongeordende distributies, afhankelijk van de parameterwaarden.

Samenvattend demonstreert het artikel hoe subtiele aanpassingen in biologische interacties en diffusieparameters de dynamiek van ecosysteempatronen fundamenteel kunnen veranderen, met name in systemen waar chemische communicatie een centrale rol speelt.