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

Résumé technique : Sur les effets de la zone de protection et du flux de population dirigé dans la dynamique proie-prédateur

Énoncé du problème
Cet article étudie un modèle spatial proie-prédateur incorporant deux mécanismes écologiques distincts : une zone de protection (refuge) inaccessible aux prédateurs et un flux de population dirigé où les prédateurs se déplacent vers des régions riches en proies. Le modèle est défini sur un domaine borné $\Omega \subset \mathbb{R}^N$, avec un sous-domaine $\Omega_0 \subset \Omega$ désigné comme zone de protection pour la proie. Les proies ($u$) diffusent de manière aléatoire dans tout $\Omega$, tandis que les prédateurs ($v$) sont confinés à $\Omega_1 = \Omega \setminus \Omega_0$ et présentent un mouvement dirigé modélisé par un terme de diffusion non linéaire $\nabla \cdot [u^2 \nabla(v/u)]$. Ce terme représente un mouvement « à longue vue » (far-sighted), où les individus se déplacent en fonction de la densité de proies à destination, contrastant avec les termes de chimiotactisme plus courants (« à courte vue » ou near-sighted) basés sur les gradients locaux ($-\nabla \cdot (v\nabla u)$). Le système est régi par les équations suivantes avec des conditions aux limites 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}$$
où $\lambda$ et $\mu$ sont les taux de croissance intrinsèques de la proie et du prédateur, respectivement, et $\alpha \ge 0$ mesure l'intensité du flux dirigé.

Méthodologie
Les auteurs emploient une combinaison de la théorie des équations d'évolution non linéaires, de l'analyse de bifurcation et de l'analyse asymptotique pour étudier les problèmes dépendants du temps et stationnaires.

1. Existence locale (Well-posedness) : Pour le problème non stationnaire (spécifiquement dans le cas $N=2$), les auteurs reformulent le système comme une équation d'évolution abstraite quaterlinéaire dans l'espace de Banach $L^2(\Omega) \times L^2(\Omega_1)$. Ils utilisent la théorie des semi-groupes analytiques développée par Yagi pour établir l'existence et l'unicité de solutions locales dans le temps. Un défi technique clé abordé est la discontinuité du terme de prédation à travers l'interface $\partial \Omega_0$ et la perte potentielle d'ellipticité de l'équation du prédateur si la densité de la proie devient complexe dans le cadre abstrait. Pour gérer cela, une fonction de coupure (cut-off) est appliquée au coefficient de la partie principale de l'opérateur du prédateur, garantissant une ellipticité uniforme tout en préservant la réalité biologique de solutions non négatives.

2. Analyse stationnaire et bifurcation : Pour le problème de l'état stationnaire, les auteurs analysent la structure globale des solutions positives (états de coexistence). Ils utilisent le théorème de bifurcation locale de Crandall-Rabinowitz pour identifier les points de bifurcation à partir de solutions sémitriviales (où une espèce est éteinte). La direction de la bifurcation (supercritique vs sous-critique) est déterminée en calculant la dérivée de la courbe de bifurcation par rapport au paramètre de bifurcation $\lambda$. Les branches de bifurcation globale sont ensuite étendues en utilisant le théorème de bifurcation globale de Rabinowitz et des estimations a priori dérivées des principes du maximum et de la régularité elliptique.

3. Analyse asymptotique ($\alpha \to \infty$) : Pour comprendre le comportement sous un mouvement dirigé fort, les auteurs effectuent une analyse de limite singulière lorsque le paramètre de flux $\alpha$ tend vers l'infini. Ils dérivent un système elliptique non linéaire limite satisfait par la densité de proie rescalée $w = \lim \alpha u$ et la densité de prédateur $v$. La structure de bifurcation globale de ce système limite est analysée en utilisant la méthode de réduction de Lyapunov-Schmidt près du point d'explosion $\lambda = 0$.

Contributions clés et résultats

• Existence locale (Well-posedness) : L'article établit l'existence de solutions uniques, non négatives et locales dans le temps pour le système non stationnaire, surmontant les difficultés posées par l'interface discontinue du refuge et le terme de flux dirigé non linéaire.
• Bifurcation à partir des états sémitriviaux :
  • Pour $\mu \ge 0$ (taux de croissance du prédateur non négatif), les solutions positives bifurquent de l'état de prédateur seul $(0, \mu)$ à $\lambda = \sigma_1(b\mu \mathbb{1}_{\Omega_1})$. Les auteurs identifient un seuil critique $\alpha^*$ tel que la bifurcation est supercritique (branche stable) si $\alpha < \alpha^*$ et sous-critique (branche initialement instable) si $\alpha > \alpha^*$.
  • Pour $\mu < 0$, la bifurcation se produit à partir de l'état de proie seule $(\lambda, 0)$ à $\lambda = |\mu|/c$, et cette bifurcation est toujours supercritique, quel que soit $\alpha$.
• Multiplicité et points de retournement : Un résultat central est que la combinaison de la zone de protection et du flux dirigé fort ($\alpha > \alpha^*$) induit une structure de point de retournement (bifurcation nœud-selle) dans la branche de bifurcation globale. Cela conduit à la multiplicité des états stationnaires de coexistence pour certaines plages de paramètres. Spécifiquement, pour $\lambda$ légèrement en dessous du seuil de bifurcation, deux états stationnaires positifs distincts peuvent exister.
• Comportement asymptotique quand $\alpha \to \infty$ :
  • Dans les régimes où $\lambda$ est suffisamment grand, un mouvement dirigé fort pousse le prédateur à l'extinction, et le système converge vers l'état de proie seule $(\lambda, 0)$.
  • Cependant, dans les régimes où $\lambda$ est faible (spécifiquement $0 < \lambda \le \sigma_1(b\mu \mathbb{1}_{\Omega_1})$), le système présente des comportements asymptotiques alternatifs. Outre la convergence vers l'état de proie seule, le système peut approcher un état où la densité de proie rescalée $w = \alpha u$ et la densité de prédateur $v$ convergent vers une solution positive du système limite.
  • L'analyse du système limite révèle qu'à mesure que $\lambda \to 0^+$, la densité de proie rescalée $w$ explose à un taux de $1/\lambda$, tandis que la densité de prédateur $v$ tend vers l'extinction.
• Effet synergique : L'article démontre que la zone de protection et le flux dirigé agissent de manière synergique. Alors que des études précédentes sur des modèles sans refuges suggéraient que le flux dirigé avait un impact limité sur la région de coexistence, cette étude montre que la présence d'une zone de protection permet au mouvement dirigé fort d'étendre significativement l'espace des paramètres pour la coexistence et d'induire des structures de bifurcation complexes (points de retournement) qui ne sont pas présentes dans les modèles purement diffusifs ou les modèles sans refuges.

Signification et revendications
Les auteurs affirment que leur travail clarifie l'interaction non triviale entre la protection spatiale et le comportement de mouvement des prédateurs. Ils soulignent qu'un mouvement dirigé fort, intuitivement attendu comme bénéfique pour le prédateur, peut paradoxalement permettre à la proie de persister à des taux de croissance intrinsèques plus faibles en raison de l'émergence de multiples configurations de coexistence stables. L'article fournit une base mathématique rigoureuse pour comprendre comment les règles de mouvement à « longue vue », distinctes de la chimiotaxie, interagissent avec les structures de refuge pour façonner la dynamique globale des systèmes proie-prédateur. Les résultats suggèrent que l'effet combiné de ces mécanismes peut inverser les directions de bifurcation et créer une bistabilité, offrant une vision plus nuancée de la persistance des populations dans des environnements hétérogènes.