Spreading speeds for prey--predator systems in a shifting environment: a short proof

Sintesi Tecnica: Velocità di Propagazione per Sistemi Preda–Predatore in un Ambiente Variabile

Enunciato del Problema
Questo articolo investiga la velocità di propagazione della componente predatrice in una classe di sistemi di reazione–diffusione preda–predatore soggetti a un ambiente variabile. Il sistema è modellato dalle seguenti equazioni:
$$\begin{cases} \partial_t u = d\partial_{xx}u + u f(x - c_1 t, u, v), \\ \partial_t v = \partial_{xx}v + v g(x - c_1 t, u, v), \end{cases}$$
dove $u$ e $v$ rappresentano rispettivamente le densità di prede e predatori, $d$ è il tasso di diffusione della preda e $c_1 > 0$ è la velocità dell'habitat variabile. Le funzioni di interazione $f$ e $g$ dipendono dalla variabile mobile $z = x - c_1 t$, che rappresenta l'eterogeneità spazio-temporale. L'obiettivo primario è determinare la velocità di propagazione verso destra del predatore ($v$) quando la preda è inizialmente stabilita in modo uniforme.

Una difficoltà centrale affrontata è che l'intero sistema manca di un diretto principio di confronto a causa della natura non cooperativa dell'interazione (la predazione riduce la crescita della preda, mentre la densità della preda guida la crescita del predatore). Inoltre, il profilo ambientale $g(z, 1, 0)$ può essere non monotono, uno scenario comune nella modellazione biologica in cui i tassi di crescita dipendono in modo unimodale dalla temperatura o da altri fattori.

Metodologia
L'autore utilizza una strategia di riduzione per aggirare la mancanza di un principio di confronto per il sistema accoppiato. Il contributo metodologico centrale è l'instaurazione di una stima puntuale (Lemma 4.1), la quale dimostra che nelle regioni in cui la densità del predatore è sufficientemente piccola, la densità della preda è vicina alla sua capacità portante (normalizzata a 1).

Nello specifico, la dimostrazione procede come segue:

1. Regione Invariante e Limiti: Utilizzando i principi di confronto scalari per ciascuna componente separatamente, l'autore stabilisce che le soluzioni rimangono all'interno di una regione invariante limitata $[0, 1] \times [0, V]$.
2. Stima Puntuale: Per contraddizione e mediante l'uso di stime paraboliche (argomenti di blow-up), si dimostra che per ogni $\alpha > 0$, esistono costanti tali che $1 - u(t, x) \leq \alpha + M_\alpha v(t, x)$ per $t$ grande. Ciò implica che dove $v$ è piccolo, $u \approx 1$.
3. Riduzione al KPP Scalare: Questa stima permette di limitare inferiormente l'equazione del predatore da un'equazione scalare di tipo Fisher–KPP con un tasso di crescita effettivo $g(x - c_1 t, 1, 0)$. Il limite superiore è derivato analogamente utilizzando la monotonia di $g$ rispetto a $v$.
4. Analisi degli Autovettori: La velocità di propagazione viene determinata analizzando l'autovettore principale generalizzato $\Lambda_1$ dell'operatore linearizzato associato al tasso di crescita effettivo $g(z, 1, 0)$ nel sistema di riferimento variabile.

Risultati Principali
Il teorema principale (Teorema 2.7) fornisce una formula esplicita per la velocità di propagazione $c^*_v$ del predatore, che dipende dalla velocità di variazione $c_1$, dai tassi di crescita asintotici $r_\pm = \lim_{z\to\pm\infty} g(z, 1, 0)$ e dall'autovettore principale generalizzato $\Lambda_1$. La velocità è data da una funzione a tratti:

• Caso 1 ($0 < c_1 \leq 2\sqrt{r_+}$): $c^*_v = 2\sqrt{r_+}$. La velocità è determinata linearmente dal tasso di crescita del fronte d'attacco.
• Caso 2 ($2\sqrt{r_+} \leq c_1 \leq 2\sqrt{\Lambda_1}$): $c^*_v = c_1$. Questo è il fenomeno di bloccaggio (locking), in cui il predatore si propaga alla stessa velocità dell'ambiente variabile.
• Caso 3 ($2\sqrt{\Lambda_1} \leq c_1 \leq 2(\sqrt{r_-} + \sqrt{\Lambda_1 - r_-})$): La velocità è data da una formula complessa che coinvolge $c_1, \Lambda_1,$ e $r_-$. Ciò corrisponde al fenomeno di attrazione non locale (nonlocal pulling), dove habitat favorevoli distanti potenziano la velocità di propagazione oltre quanto suggerito dal fronte locale.
• Caso 4 ($c_1 \geq 2(\sqrt{r_-} + \sqrt{\Lambda_1 - r_-})$): $c^*_v = 2\sqrt{r_-}$.

Il lavoro stabilisce anche i comportamenti asintotici (Corollari 2.10 e 2.13). Mostra che davanti al fronte, la preda si avvicina alla sua capacità portante e il predatore svanisce. Dietro il fronte, la soluzione converge a un equilibrio di coesistenza determinato dall'ambiente limite ($f_\pm, g_\pm$) osservato nel sistema di riferimento mobile, a condizione che i sistemi ODE limite ammettano funzioni di Lyapunov strette.

Significatività e Rivendicazioni
L'articolo sostiene di fornire una riduzione rigorosa da un sistema preda–predatore non cooperativo a un'associata equazione scalare Fisher–KPP in un ambiente variabile, anche quando il profilo ambientale è non monotono.

• Estensione della Teoria Scalare: I risultati confermano che i fenomeni precedentemente scoperti per le equazioni KPP scalari — specificamente il bloccaggio (locking) e l'attrazione non locale (nonlocal pulling) — persistono al livello dell'invasione del predatore in sistemi non cooperativi.
• Gestione della Non-Monotonia: La metodologia affronta con successo le difficoltà introdotte da ambienti variabili non monotoni, che sono biologicamente significativi ma matematicamente impegnativi.
• Conferma della Congettura: Il risultato principale conferma una congettura riguardante il caso di bloccaggio citata in un lavoro precedente [16, Remark 3], che aveva utilizzato un approccio Hamilton–Jacobi per un caso specifico monotono. Questo articolo estende la validità del fenomeno di bloccaggio a contesti non monotoni.
• Limitazioni: L'autore osserva che, sebbene la velocità di propagazione sia caratterizzata, la correzione del ritardo logaritmico (correzione di Bramson) per la posizione del fronte della soluzione rimane un problema aperto per i sistemi preda–predatore, e l'attuale argomento non affronta questo aspetto.

L'articolo conclude applicando questi risultati generali a un modello di Lotka–Volterra con un tasso di conversione variabile, dimostrando il calcolo esplicito delle velocità di propagazione e la verifica della convergenza asintotica verso stati di coesistenza.