On the effects of protection zone and directed population flux in prey-predator dynamics

Sintesi Tecnica: Sull'effetto della zona di protezione e del flusso di popolazione direzionato nella dinamica preda-predatore

Definizione del Problema
Questo articolo investiga un modello spaziale preda-predatore che incorpora due distinti meccanismi ecologici: una zona di protezione (rifugio) inaccessibile ai predatori e un flusso di popolazione direzionato in cui i predatori si muovono verso regioni ricche di prede. Il modello è definito su un dominio limitato $\Omega \subset \mathbb{R}^N$, con un sottodominio $\Omega_0 \subset \Omega$ designato come zona di protezione per la preda. Le prede ($u$) diffondono casualmente in tutto $\Omega$, mentre i predatori ($v$) sono confinati in $\Omega_1 = \Omega \setminus \Omega_0$ ed esibiscono un movimento direzionato modellato da un termine di diffusione non lineare $\nabla \cdot [u^2 \nabla(v/u)]$. Questo termine rappresenta un movimento "a lungo raggio" (far-sighted), dove gli individui si muovono in base alla densità delle prede nella destinazione, contrastando con i più comuni termini di chemiotassi "a corto raggio" (near-sighted) ($-\nabla \cdot (v\nabla u)$) basati su gradienti locali. Il sistema è governato dalle seguenti equazioni con condizioni al contorno di Neumann:
$$\begin{cases} u_t = \Delta u + u(\lambda - u - b \mathbb{1}_{\Omega_1}v), & (x, t) \in \Omega \times (0, T), \\ v_t = \Delta v + \alpha \nabla \cdot \left[ u^2 \nabla \left( \frac{v}{u} \right) \right] + v(\mu + cu - v), & (x, t) \in \Omega_1 \times (0, T), \end{cases}$$
dove $\lambda$ e $\mu$ sono i tassi di crescita intrinseci di preda e predatore, rispettivamente, e $\alpha \ge 0$ misura l'intensità del flusso direzionato.

Metodologia
Gli autori utilizzano una combinazione di teoria delle equazioni di evoluzione non lineari, analisi di biforcazione e analisi asintotica per studiare sia il problema tempo-dipendente che quello stazionario.

1. Ben determinazione locale nel tempo: Per il problema non stazionario (specificamente nel caso $N=2$), gli autori riformulano il sistema come un'equazione di evoluzione astratta quasilineare nello spazio di Banach $L^2(\Omega) \times L^2(\Omega_1)$. Utilizzano la teoria dei semigruppi analitici sviluppata da Yagi per stabilire l'esistenza e l'unicità di soluzioni locali nel tempo. Una sfida tecnica chiave affrontata è la discontinuità del termine di predazione attraverso l'interfaccia $\partial \Omega_0$ e la potenziale perdita di ellitticità nell'equazione del predatore se la densità della preda diventa complessa nell'ambito astratto. Per gestire questo, viene applicata una funzione di taglio (cut-off) al coefficiente della parte principale dell'operatore del predatore, garantendo un'ellitticità uniforme pur preservando la realtà biologica delle soluzioni non negative.

2. Analisi stazionaria e biforcazione: Per il problema dello stato stazionario, gli autori analizzano la struttura globale delle soluzioni positive (stati di coesistenza). Utilizzano il teorema di biforcazione locale di Crandall-Rabinowitz per identificare i punti di biforcazione dalle soluzioni semitriviali (dove una specie è estinta). La direzione della biforcazione (supercritica vs. subcritica) è determinata calcolando la derivata della curva di biforcazione rispetto al parametro di biforcazione $\lambda$. I rami di biforcazione globale vengono poi estesi utilizzando il teorema di biforcazione globale di Rabinowitz e stime a priori derivate dai principi di massimo e dalla regolarità ellittica.

3. Analisi asintotica ($\alpha \to \infty$): Per comprendere il comportamento sotto un forte movimento direzionato, gli autori eseguono un'analisi del limite singolare al tendere del parametro di flusso $\alpha$ all'infinito. Derivano un sistema ellittico non lineare limite soddisfatto dalla densità di preda riscalata $w = \lim \alpha u$ e dalla densità del predatore $v$. La struttura di biforcazione globale di questo sistema limite è analizzata utilizzando il metodo di riduzione di Lyapunov-Schmidt vicino al punto di esplosione (blow-up) $\lambda = 0$.

Contributi Chiave e Risultati

• Ben determinazione locale: L'articolo stabilisce l'esistenza di soluzioni uniche, non negative e locali nel tempo per il sistema non stazionario, superando le difficoltà poste dall'interfaccia discontinua del rifugio e dal termine di flusso direzionato non lineare.
• Biforcazione da stati semitriviali:
  • Per $\mu \ge 0$ (tasso di crescita del predatore non negativo), le soluzioni positive biforcano dallo stato solo-predatore $(0, \mu)$ a $\lambda = \sigma_1(b\mu \mathbb{1}_{\Omega_1})$. Gli autori identificano una soglia critica $\alpha^*$ tale che la biforcazione sia supercritica (ramo stabile) se $\alpha < \alpha^*$ e subcritica (ramo inizialmente instabile) se $\alpha > \alpha^*$.
  • Per $\mu < 0$, la biforcazione avviene dallo stato solo-preda $(\lambda, 0)$ a $\lambda = |\mu|/c$, e questa biforcazione è sempre supercritica indipendentemente da $\alpha$.
• Molteplicità e punti di svolta: Un risultato centrale è che la combinazione della zona di protezione e del forte flusso direzionato ($\alpha > \alpha^*$) induce una struttura a punto di svolta (biforcazione saddle-node) nel ramo di biforcazione globale. Ciò porta alla molteplicità degli stati stazionari di coesistenza per determinati intervalli di parametri. Nello specifico, per $\lambda$ leggermente al di sotto della soglia di biforcazione, possono esistere due distinti stati stazionari positivi.
• Comportamento asintotico al tendere di $\alpha \to \infty$:
  • Nei regimi in cui $\lambda$ è sufficientemente grande, il forte movimento direzionato guida il predatore all'estinzione, e il sistema converge allo stato solo-preda $(\lambda, 0)$.
  • Tuttavia, nei regimi in cui $\lambda$ è piccolo (specificamente $0 < \lambda \le \sigma_1(b\mu \mathbb{1}_{\Omega_1})$), il sistema esibisce comportamenti asintotici alternativi. Oltre a convergere allo stato solo-preda, il sistema può avvicinarsi a uno stato in cui la densità di preda riscalata $w = \alpha u$ e la densità del predatore $v$ convergono a una soluzione positiva del sistema limite.
  • L'analisi del sistema limite rivela che quando $\lambda \to 0^+$, la densità di preda riscalata $w$ esplode con un tasso di $1/\lambda$, mentre la densità del predatore $v$ tende all'estinzione.
• Effetto sinergico: L'articolo dimostra che la zona di protezione e il flusso direzionato agiscono sinergicamente. Mentre studi precedenti su modelli senza rifugi suggerivano che il flusso direzionato avesse un impatto limitato sulla regione di coesistenza, questo studio mostra che la presenza di una zona di protezione permette al forte movimento direzionato di espandere significativamente lo spazio dei parametri per la coesistenza e indurre strutture di biforcazione complesse (punti di svolta) che non sono presenti in modelli puramente diffusivi o modelli senza rifugi.

Significatività e Rivendicazioni
Gli autori affermano che il loro lavoro chiarisce l'interazione non banale tra protezione spaziale e comportamento di movimento del predatore. Evidenziano come il forte movimento direzionato, intuitivamente atteso come un beneficio per il predatore, possa paradossalmente permettere alla preda di persistere a tassi di crescita intrinseci inferiori grazie all'emergere di molteplici configurazioni di coesistenza stabili. L'articolo fornisce una base matematica rigorosa per comprendere come le regole di movimento "a lungo raggio", distinte dalla chemiotassi, interagiscano con le strutture di rifugio per plasmare la dinamica globale dei sistemi preda-predatore. I risultati suggeriscono che l'effetto combinato di questi meccanismi può invertire le direzioni di biforcazione e creare bistabilità, offrendo una visione più sfumata della persistenza delle popolazioni in ambienti eterogenei.