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

Resumen Técnico: Sobre los efectos de la zona de protección y el flujo poblacional dirigido en la dinámica depredador-presa

Planteamiento del Problema
Este artículo investiga un modelo espacial de depredador-presa que incorpora dos mecanismos ecológicos distintos: una zona de protección (refugio) inaccesible para los depredadores y un flujo poblacional dirigido donde los depredadores se desplazan hacia regiones ricas en presas. El modelo se define en un dominio acotado $\Omega \subset \mathbb{R}^N$, con un subdominio $\Omega_0 \subset \Omega$ designado como zona de protección para la presa. Las presas ($u$) difunden aleatoriamente por todo $\Omega$, mientras que los depredadores ($v$) están confinados a $\Omega_1 = \Omega \setminus \Omega_0$ y exhiben un movimiento dirigido modelado por un término de difusión no lineal $\nabla \cdot [u^2 \nabla(v/u)]$. Este término representa un movimiento de "visión lejana" (far-sighted), donde los individuos se mueven basándose en la densidad de la presa en el destino, contrastando con los términos de quimiotaxis de "visión cercana" (near-sighted) más comunes ($-\nabla \cdot (v\nabla u)$) basados en gradientes locales. El sistema está gobernado por las siguientes ecuaciones con condiciones 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}$$
donde $\lambda$ y $\mu$ son las tasas de crecimiento intrínseco de la presa y los depredadores, respectivamente, y $\alpha \ge 0$ mide la fuerza del flujo dirigido.

Metodología
Los autores emplean una combinación de la teoría de ecuaciones de evolución no lineales, análisis de bifurcación y análisis asintótico para estudiar tanto el problema dependiente del tiempo como el estacionario.

1. Bien-posedidad local en el tiempo: Para el problema no estacionario (específicamente en el caso $N=2$), los autores reformulan el sistema como una ecuación de evolución cuasilineal abstracta en el espacio de Banach $L^2(\Omega) \times L^2(\Omega_1)$. Utilizan la teoría de semigrupos analíticos desarrollada por Yagi para establecer la existencia y unicidad de soluciones locales en el tiempo. Un desafío técnico clave abordado es la discontinuidad del término de depredación a través de la interfaz $\partial \Omega_0$ y la potencial pérdida de elipticidad en la ecuación del depredador si la densidad de la presa se vuelve compleja en el entorno abstracto. Para manejar esto, se aplica una función de corte al coeficiente de la parte principal del operador del depredador, asegurando la elipticidad uniforme mientras se preserva la realidad biológica de las soluciones no negativas.

2. Análisis Estacionario y Bifurcación: Para el problema de estado estacionario, los autores analizan la estructura global de las soluciones positivas (estados de coexistencia). Utilizan el teorema de bifurcación local de Crandall-Rabinowitz para identificar puntos de bifurcación desde soluciones semitriviales (donde una especie está extinta). La dirección de la bifurcación (supercrítica vs. subcrítica) se determina calculando la derivada de la curva de bifurcación con respecto al parámetro de bifurcación $\lambda$. Las ramas de bifurcación global se extienden utilizando el teorema de bifurcación global de Rabinowitz y estimaciones a priori derivadas de principios de máximo y regularidad elíptica.

3. Análisis Asintótico ($\alpha \to \infty$): Para comprender el comportamiento bajo un movimiento dirigido fuerte, los autores realizan un análisis de límite singular a medida que el parámetro de flujo $\alpha$ tiende al infinito. Derivan un sistema elíptico no lineal límite satisfecho por la densidad de la presa reescalada $w = \lim \alpha u$ y la densidad del depredador $v$. La estructura de bifurcación global de este sistema límite se analiza mediante el método de reducción de Lyapunov-Schmidt cerca del punto de explosión (blow-up) $\lambda = 0$.

Contribuciones Clave y Resultados

• Bien-posedidad Local: El artículo establece la existencia de soluciones únicas y no negativas para el sistema local en el tiempo, superando las dificultades planteadas por la interfaz discontinua del refugio y el término de flujo dirigido no lineal.
• Bifurcación desde Estados Semitriviales:
  • Para $\mu \ge 0$ (tasa de crecimiento del depredador no negativa), las soluciones positivas bifurcan desde el estado de solo depredador $(0, \mu)$ en $\lambda = \sigma_1(b\mu \mathbb{1}_{\Omega_1})$. Los autores identifican un umbral crítico $\alpha^*$ tal que la bifurcación es supercrítica (rama estable) si $\alpha < \alpha^*$ y subcrítica (rama inicialmente inestable) si $\alpha > \alpha^*$.
  • Para $\mu < 0$, la bifurcación ocurre desde el estado de solo presa $(\lambda, 0)$ en $\lambda = |\mu|/c$, y esta bifurcación es siempre supercrítica independientemente de $\alpha$.
• Multiplicidad y Puntos de Retorno (Turning Points): Un hallazgo central es que la combinación de la zona de protección y el flujo dirigido fuerte ($\alpha > \alpha^*$) induce una estructura de punto de retorno (bifurcación de tipo saddle-node) en la rama de bifurcación global. Esto conduce a la multiplicidad de estados estacionarios de coexistencia para ciertos rangos de parámetros. Específicamente, para un $\lambda$ ligeramente por debajo del umbral de bifurcación, pueden existir dos estados estacionarios positivos distintos.
• Comportamiento Asintótico cuando $\alpha \to \infty$:
  • En regímenes donde $\lambda$ es suficientemente grande, el movimiento dirigido fuerte impulsa la extinción del depredador, y el sistema converge al estado de solo presa $(\lambda, 0)$.
  • Sin embargo, en regímenes donde $\lambda$ es pequeño (específicamente $0 < \lambda \le \sigma_1(b\mu \mathbb{1}_{\Omega_1})$), el sistema exhibe comportamientos asintóticos alternativos. Además de converger al estado de solo presa, el sistema puede aproximarse a un estado donde la densidad de la presa reescalada $w = \alpha u$ y la densidad del depredador $v$ convergen a una solución positiva del sistema límite.
  • El análisis del sistema límite revela que a medida que $\lambda \to 0^+$, la densidad de la presa reescalada $w$ explota a una tasa de $1/\lambda$, mientras que la densidad del depredador $v$ tiende a la extinción.
• Efecto Sinérgico: El artículo demuestra que la zona de protección y el flujo dirigido actúan sinérgicamente. Mientras que estudios previos sobre modelos sin refugios sugirieron que el flujo dirigido tiene un impacto limitado en la región de coexistencia, este estudio muestra que la presencia de una zona de protección permite que el movimiento dirigido fuerte expanda significamente el espacio de parámetros para la coexistencia e induzca estructuras de bifurcación complejas (puntos de retorno) que no están presentes en modelos puramente difusivos o modelos sin refugios.

Significancia y Reivindicaciones
Los autores afirman que su trabajo esclarece la interacción no trivial entre la protección espacial y el comportamiento de movimiento del depredador. Destacan que el movimiento dirigido fuerte, que intuitivamente se esperaría que beneficie al depredador, puede paradójicamente permitir que la presa persista con tasas de crecimiento intrínseco más bajas debido a la emergencia de múltiples configuraciones de coexistencia estables. El artículo proporciona una base matemática rigurosa para comprender cómo las reglas de movimiento de "visión lejana", distintas de la quimiotaxis, interactúan con las estructuras de refugio para dar forma a la dinámica global de los sistemas depredador-presa. Los resultados sugieren que el efecto combinado de estos mecanismos puede revertir las direcciones de bifurcación y crear bistabilidad, ofreciendo una visión más matizada de la persistencia poblacional en entornos heterogéneos.