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

Resumo Técnico: Sobre os efeitos da zona de proteção e do fluxo populacional direcionado na dinâmica predador-presa

Enunciado do Problema
Este artigo investiga um modelo espacial de predador-presa que incorpora dois mecanismos ecológicos distintos: uma zona de proteção (refúgio) inacessível aos predadores e um fluxo populacional direcionado onde os predadores se deslocam para regiões ricas em presas. O modelo é definido em um domínio limitado $\Omega \subset \mathbb{R}^N$, com um subdomínio $\Omega_0 \subset \Omega$ designado como uma zona de proteção para a presa. As presas ($u$) difundem-se aleatoriamente por todo o $\Omega$, enquanto os predadores ($v$) estão confinados a $\Omega_1 = \Omega \setminus \Omega_0$ e exibem um movimento direcionado modelado por um termo de difusão não linear $\nabla \cdot [u^2 \nabla(v/u)]$. Este termo representa um movimento "vidente" (far-sighted), onde os indivíduos se movem com base na densidade de presas no destino, contrastando com os termos de quimiotaxia "videntes de curto alcance" (near-sighted) mais comuns ($-\nabla \cdot (v\nabla u)$). O sistema é governado pelas seguintes equações com condições de contorno de 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}$$
onde $\lambda$ e $\mu$ são as taxas de crescimento intrínsecas da presa e dos predadores, respectivamente, e $\alpha \ge 0$ mede a força do fluxo direcionado.

Metodologia
Os autores empregam uma combinação de teoria de equações de evolução não lineares, análise de bifurcação e análise assintótica para estudar tanto o problema dependente do tempo quanto o estacionário.

1. Bem-postura Local no Tempo: Para o problema não estacionário (especificamente para o caso $N=2$), os autores reformulam o sistema como uma equação de evolução quasilinear abstrata no espaço de Banach $L^2(\Omega) \times L^2(\Omega_1)$. Eles utilizam a teoria de semigrupos analíticos desenvolvida por Yagi para estabelecer a existência e unicidade de soluções locais no tempo. Um desafio técnico fundamental abordado é a descontinuidade do termo de predação através da interface $\partial \Omega_0$ e o potencial de perda de elipticidade na equação do predador caso a densidade da presa se torne complexa no contexto abstrato. Para lidar com isso, uma função de corte (cut-off) é aplicada ao coeficiente da parte principal do operador do predador, garantindo elipticidade uniforme enquanto preserva a realidade biológica de soluções não negativas.

2. Análise Estacionária e Bifurcação: Para o problema do estado estacionário, os autores analisam a estrutura global das soluções positivas (estados de coexistência). Eles utilizam o teorema de bifurcação local de Crandall-Rabinowitz para identificar pontos de bifurcação a partir de soluções semitriviais (onde uma espécie está extinta). A direção da bifurcação (supercrítica vs. subcrítica) é determinada pelo cálculo da derivada da curva de bifurcação em relação ao parâmetro de bifurcação $\lambda$. Os ramos de bifurcação global são então estendidos usando o teorema de bifurcação global de Rabinowitz e estimativas a priori derivadas de princípios de máximo e regularidade elíptica.

3. Análise Assintótica ($\alpha \to \infty$): Para compreender o comportamento sob movimento direcionado forte, os autores realizam uma análise de limite singular conforme o parâmetro de fluxo $\alpha$ tende ao infinito. Eles derivam um sistema elíptico não linear limitante satisfeito pela densidade de presa reescalonada $w = \lim \alpha u$ e pela densidade de predador $v$. A estrutura de bifurcação global deste sistema limitante é analisada usando o método de redução de Lyapunov-Schmidt próximo ao ponto de explosão (blow-up) $\lambda = 0$.

Principais Contribuições e Resultados

• Bem-postura Local no Tempo: O artigo estabelece a existência de soluções únicas e não negativas de tempo local para o sistema não estacionário, superando as dificuldades impostas pela interface descontínua do refúgio e pelo termo de fluxo direcionado não linear.
• Bifurcação de Estados Semitriviais:
  • Para $\mu \ge 0$ (taxa de crescimento do predador não negativa), soluções positivas bifurcam do estado apenas de predadores $(0, \mu)$ em $\lambda = \sigma_1(b\mu \mathbb{1}_{\Omega_1})$. Os autores identificam um limiar crítico $\alpha^*$ tal que a bifurcação é supercrítica (ramo estável) se $\alpha < \alpha^*$ e subcrítica (ramo inicialmente instável) se $\alpha > \alpha^*$.
  • Para $\mu < 0$, a bifurcação ocorre a partir do estado apenas de presa $(\lambda, 0)$ em $\lambda = |\mu|/c$, e esta bifurcação é sempre supercrítica, independentemente de $\alpha$.
• Multiplicidade e Pontos de Retorno (Turning Points): Uma descoberta central é que a combinação da zona de proteção e o fluxo direcionado forte ($\alpha > \alpha^*$) induz uma estrutura de ponto de retorno (bifurcação de nó-sela) no ramo de bifurcação global. Isso leva à multiplicidade de estados estacionários de coexistência para certos intervalos de parâmetros. Especificamente, para $\lambda$ ligeiramente abaixo do limiar de bifurcação, dois estados estacionários positivos distintos podem existir.
• Comportamento Assintótico conforme $\alpha \to \infty$:
  • Em regimes onde $\lambda$ é suficientemente grande, o movimento direcionado forte conduz o predador à extinção, e o sistema converge para o estado apenas de presa $(\lambda, 0)$.
  • No entanto, em regimes onde $\lambda$ é pequeno (especificamente $0 < \lambda \le \sigma_1(b\mu \mathbb{1}_{\Omega_1})$), o sistema exibe comportamentos assintóticos alternativos. Além de convergir para o estado apenas de presa, o sistema pode se aproximar de um estado onde a densidade de presa reescalonada $w = \alpha u$ e a densidade de predador $v$ convergem para uma solução positiva do sistema limitante.
  • A análise do sistema limitante revela que, conforme $\lambda \to 0^+$, a densidade de presa reescalonada $w$ explode a uma taxa de $1/\lambda$, enquanto a densidade de predador $v$ tende à extinção.
• Efeito Sinérgico: O artigo demonstra que a zona de proteção e o fluxo direcionado atuam sinergicamente. Enquanto estudos anteriores sobre modelos sem refúgios sugeriam que o fluxo direcionado tem impacto limitado na região de coexistência, este estudo mostra que a presença de uma zona de proteção permite que o movimento direcionado forte expanda significativamente o espaço de parâmetros para a coexistência e induza estruturas de bifurcação complexas (pontos de retorno) que não estão presentes em modelos puramente difusivos ou modelos sem refúgios.

Significância e Alegações
Os autores alegam que seu trabalho esclarece a interação não trivial entre a proteção espacial e o comportamento de movimento do predador. Eles destacam que o movimento direcionado forte, intuitivamente esperado como benéfico para o predador, pode paradoxalmente permitir que a presa persista em taxas de crescimento intrínseco mais baixas devido à emergência de múltiplas configurações de coexistência estáveis. O artigo fornece uma base matemática rigorosa para entender como regras de movimento "videntes" (far-sighted), distintas da quimiotaxia, interagem com estruturas de refúgio para moldar a dinâmica global de sistemas predador-presa. Os resultados sugerem que o efeito combinado desses mecanismos pode reverter direções de bifurcação e criar bistabilidade, oferecendo uma visão mais matizada da persistência populacional em ambientes heterogêneos.