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

Resumo Técnico: Velocidades de Propagação para Sistemas Presa–Predador em um Ambiente em Deslocamento

Enunciado do Problema
Este artigo investiga a velocidade de propagação do componente predador em uma classe de sistemas de reação-difusão presa-predador sujeitos a um ambiente em deslocamento. O sistema é modelado pelas seguintes equações:
$$\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}$$
onde $u$ e $v$ representam as densidades de presa e predador, respectivamente, $d$ é a taxa de difusão da presa, e $c_1 > 0$ é a velocidade do habitat em deslocamento. As funções de interação $f$ e $g$ dependem da variável móvel $z = x - c_1 t$, representando a heterogeneidade espaço-temporal. O objetivo principal é determinar a velocidade de propagação para a direita do predador ($v$) quando a presa está inicialmente estabelecida de forma uniforme.

Uma dificuldade central abordada é que o sistema completo carece de um princípio de comparação direto devido à natureza não cooperativa da interação (a predação reduz o crescimento da presa, enquanto a densidade da presa impulsiona o crescimento do predador). Além disso, o perfil ambiental $g(z, 1, 0)$ pode ser não monotônico, um cenário comum na modelagem biológica onde as taxas de crescimento dependem unimodalmente da temperatura ou de outros fatores.

Metodologia
O autor emprega uma estratégia de redução para contornar a falta de um princípio de comparação para o sistema acoplado. A principal contribuição metodológica é o estabelecimento de uma estimativa pontual (Lema 4.1), que demonstra que, em regiões onde a densidade do predador é suficientemente pequena, a densidade da presa está próxima de sua capacidade de carga (normalizada para 1).

Especificamente, a prova procede da seguinte forma:

1. Região Invariante e Limites: Utilizando princípios de comparação escalares para cada componente separadamente, o autor estabelece que as soluções permanecem dentro de uma região invariante limitada $[0, 1] \times [0, V]$.
2. Estimativa Pontual: Por contradição e pelo uso de estimativas parabólicas (argumentos de explosão/blow-up), demonstra-se que, para qualquer $\alpha > 0$, existem constantes tais que $1 - u(t, x) \leq \alpha + M_\alpha v(t, x)$ para $t$ grande. Isso implica que, onde $v$ é pequeno, $u \approx 1$.
3. Redução para KPP Escalar: Esta estimativa permite que a equação do predador seja limitada inferiormente por uma equação escalar de Fisher–KPP com uma taxa de crescimento efetiva $g(x - c_1 t, 1, 0)$. O limite superior é derivado de forma semelhante usando a monotonicidade de $g$ em relação a $v$.
4. Análise de Autovalores: A velocidade de propagação é determinada pela análise do autovalor principal generalizado $\Lambda_1$ do operador linearizado associado à taxa de crescimento efetiva $g(z, 1, 0)$ no referencial de deslocamento.

Principais Resultados
O teorema principal (Teorema 2.7) fornece uma fórmula explícita para a velocidade de propagação $c^*_v$ do predador, que depende da velocidade de deslocamento $c_1$, das taxas de crescimento assintóticas $r_\pm = \lim_{z\to\pm\infty} g(z, 1, 0)$, e do autovalor principal generalizado $\Lambda_1$. A velocidade é dada por uma função definida em partes:

• Caso 1 ($0 < c_1 \leq 2\sqrt{r_+}$): $c^*_v = 2\sqrt{r_+}$. A velocidade é determinada linearmente pela taxa de crescimento da borda de frente.
• Caso 2 ($2\sqrt{r_+} \leq c_1 \leq 2\sqrt{\Lambda_1}$): $c^*_v = c_1$. Este é o fenômeno de travamento (locking), onde o predador propaga-se na mesma velocidade que o ambiente em deslocamento.
• Caso 3 ($2\sqrt{\Lambda_1} \leq c_1 \leq 2(\sqrt{r_-} + \sqrt{\Lambda_1 - r_-})$): A velocidade é dada por uma fórmula complexa envolvendo $c_1, \Lambda_1,$ e $r_-$. Isso corresponde ao fenômeno de tração não local (nonlocal pulling), onde habitats favoráveis distantes aumentam a velocidade de propagação além do que a borda de frente local sugeriria.
• Caso 4 ($c_1 \geq 2(\sqrt{r_-} + \sqrt{\Lambda_1 - r_-})$): $c^*_v = 2\sqrt{r_-}$.

O artigo também estabelece comportamentos assintóticos (Corolários 2.10 e 2.13). Mostra que à frente da frente, a presa se aproxima de sua capacidade de carga e o predador desaparece. Atrás da frente, a solução converge para um equilíbrio de coexistência determinado pelo ambiente limitante ($f_\pm, g_\pm$) visto no referencial móvel, desde que os sistemas diferenciais ordinários (ODEs) limitantes admitam funções de Lyapunov estritas.

Significância e Alegações
O artigo alega fornecer uma redução rigorosa de um sistema presa-predador não cooperativo para uma equação de Fisher–KPP escalar associada em um ambiente em deslocamento, mesmo quando o perfil ambiental é não monotônico.

• Extensão da Teoria Escalar: Os resultados confirmam que fenômenos anteriormente descobertos para equações KPP escalares — especificamente o travamento (locking) e a tração não local (nonlocal pulling) — persistem ao nível da invasão do predador em sistemas não cooperativos.
• Tratamento da Não Monotonicidade: A metodologia aborda com sucesso as dificuldades introduzidas por ambientes em deslocamento não monotônicos, que são biologicamente significativos, mas matematicamente desafiadores.
• Confirmação de Conjectura: O resultado principal confirma uma conjectura sobre o caso de travamento estabelecida em trabalho anterior [16, Remark 3], que utilizou uma abordagem de Hamilton–Jacobi para um caso monotônico específico. Este artigo estende a validade do fenômeno de travamento para configurações não monotônicas.
• Limitações: O autor observa que, embora a velocidade de propagação seja caracterizada, a correção de atraso logarítmico (correção de Bramson) para a posição da frente da solução permanece um problema em aberto para sistemas presa-predador, e o argumento atual não aborda isso.

O artigo conclui aplicando estes resultados gerais a um modelo de Lotka–Volterra com uma taxa de conversão em deslocamento, demonstrando o cálculo explícito das velocidades de propagação e a verificação da convergência assintótica para estados de coexistência.