Memory-Delay Stability Switching and Ecological Thresholds in a Dimensionally Balanced Fractional-Order Predator-Prey Model

Resumen Técnico: Cambio de Estabilidad por Retardo de Memoria y Umbrales Ecológicos en un Modelo de Depredador-Presa de Orden Fraccionario Dimensionalmente Balanceado

1. Planteamiento del Problema y Motivación

El artículo aborda un fallo metodológico crítico en muchos modelos ecológicos de orden fraccionario existentes: la inconsistencia dimensional. En los modelos clásicos de orden entero, la tasa de cambio $dx/dt$ tiene dimensiones de $[Masa][Tiempo]^{-1}$. Sin embargo, la derivada de Caputo de orden fraccionario $0 < \alpha < 1$, denotada como $D_t^\alpha$, posee inherentemente dimensiones de $[Masa][Tiempo]^{-\alpha}$. Reemplazar directamente la derivada de orden entero por una fraccionaria sin ajustar los parámetros del sistema resulta en un desajuste dimensional, lo que hace que el modelo carezca de sentido físico.

Además, los autores argumentan que los sistemas ecológicos están gobernados por múltiples mecanismos interactuantes que a menudo se estudian de forma aislada:

• Efectos de memoria: Las poblaciones biológicas responden a la historia pasada, no solo a los estados actuales.
• Retardos temporales: El crecimiento del depredador suele tener un desfase respecto al consumo de la presa debido a la gestación o la digestión.
• Factores conductuales: El miedo a la depredación y el refugio de la presa alteran significativamente la dinámica poblacional.
• Gestión: La cosecha no lineal (por ejemplo, la pesca) es un mecanismo de control realista.

El estudio pretende formular un modelo de depredador-presa de orden fraccionario dimensionalmente balanceado que integre estos factores (miedo, refugio, respuesta funcional de Holling tipo III, cosecha no lineal y retardo temporal) asegurando la consistencia física.

2. Metodología

Formulación del Modelo:
Los autores construyen un sistema dimensionalmente balanceado reconstruyendo los parámetros de tasa dependientes del tiempo. En lugar de simplemente reemplazar $r$ con $r^\alpha$, definen coeficientes de tasa fraccionarios (por ejemplo, $r_\alpha = r^\alpha$, $d_\alpha = d^\alpha$, $q_\alpha = q^\alpha$) para asegurar que cada término en las ecuaciones diferenciales tenga la dimensión $[Masa][Tiempo]^{-\alpha}$. Se presta especial atención a los parámetros de dimensiones mixtas (como el coeficiente de autolimitación $\sigma$), redefiniéndolos como $\sigma_\alpha = \sigma d^{\alpha-1}$ para mantener la homogeneidad.

El sistema resultante (Ecuación 8) es:
$$\begin{aligned} D_t^\alpha X(t) &= \frac{r_\alpha X(t)}{1 + \kappa Y(t)} \left(1 - \frac{X(t)}{K}\right) - H_\alpha(X(t))Y(t) - \frac{q_\alpha X^2(t)}{c + X(t)} \\ D_t^\alpha Y(t) &= e H_\alpha(X(t-\tau))Y(t) - d_\alpha Y(t) - \sigma_\alpha Y^2(t) \end{aligned}$$
donde $H_\alpha$ es la respuesta funcional de Holling tipo III dimensionalmente balanceada, $\tau$ es el retardo temporal y $\alpha$ es el orden fraccionario.

Marco Analítico:

• Bien posed de la solución: Los autores demuestran la existencia local, unicidad, positividad y acotación última de las soluciones utilizando la formulación integral de Volterra y principios de comparación fraccionaria.
• Análisis de Equilibrio: Derivan el equilibrio de solo presa y el equilibrio de coexistencia. Se define un número de invasión del depredador ($R_0$) para determinar el umbral de persistencia o extinción del depredador.
• Análisis de Estabilidad:
  • Sin retardo ($\tau=0$): La estabilidad se analiza mediante el criterio de Matignon, que requiere que los autovalores $\mu$ del Jacobiano satisfagan $|\arg(\mu)| > \alpha\pi/2$.
  • Con retardo ($\tau>0$): Se deriva la ecuación característica. Los autores resuelven para las frecuencias críticas y los valores de retardo críticos ($\tau_j$) donde las raíces cruzan el eje imaginario. Se establecen condiciones de transversalidad para determinar la dirección del cambio de estabilidad.

Simulación Numérica:
Las simulaciones se realizan utilizando un método de predictor-corrector para ecuaciones diferenciales de retardo de Caputo. Los autores reconstruyen explícitamente los coeficientes dimensionalmente balanceados para cada simulación en lugar de utilizar una sustitución de parámetros ingenua.

3. Contribuciones Clave

1. Consistencia Dimensional: El artículo proporciona un marco rigurooso para construir modelos ecológicos fraccionarios que sean físicamente consistentes desde el punto de vista dimensional, corrigiendo la práctica común de la sustitución directa de la derivada.
2. Clasificación Exhaustiva de la Estabilidad: El estudio clasifica y verifica numéricamente seis escenarios distintos de cambio de estabilidad inducidos por el retardo para el equilibrio de coexistencia:
   • Permanecer estable
   • Estable a inestable
   • Estable-inestable-estable (Burbuja de estabilidad)
   • Permanecer inestable
   • Inestable a estable
   • Inestable-estable-inestable (Burbuja de inestabilidad)
3. Interacción Memoria-Retardo: El análisis demuestra cómo el orden fraccionario $\alpha$ (memoria) interactúa con el retardo temporal $\tau$. El retardo crítico $\tau_0(\alpha)$ se muestra como una función de $\alpha$, indicando que los efectos de la memoria pueden estabilizar o desestabilizar el sistema dependiendo de la magnitud del retardo.
4. Umbrales de Gestión Ecológica: El artículo vincula la estabilidad matemática con la gestión ecológica mediante el análisis del impacto de la cosecha no lineal ($q$) en el número de invasión del depredador ($R_0$).

4. Resultados

• Existencia y Persistencia: El modelo asegura soluciones positivas y acotadas. Se establece una condición de umbral $R_0(q, \alpha) > 1$ para la invasión del depredador. Los autores demuestran que aumentar la intensidad de la cosecha $q$ disminuye estrictamente $R_0$, lo que puede conducir a la extinción del depredador.
• Cambio de Estabilidad:
  • Para el sistema sin retardo, la estabilidad depende de la traza y el determinante del Jacobiano y del orden fraccionario $\alpha$.
  • Para el sistema con retardo, la introducción de $\tau$ puede inducir bifurcaciones de Hopf. El estudio confirma que la estabilidad no es monótona con respecto al retardo; el sistema puede transicionar de estable a inestable y de nuevo a estable (o viceversa) a medida que $\tau$ aumenta.
• Estabilización Inducida por la Memoria: Los resultados numéricos muestran que para un retardo fijo que causa oscilaciones en un modelo de orden entero, reducir el orden fraccionario $\alpha$ (aumentar la memoria) puede suprimir las oscilaciones y estabilizar el equilibrio de coexistencia.
• Sensibilidad: El análisis de sensibilidad indica que la dinámica del sistema es altamente sensible a las variaciones en las tasas de cosecha y en las densidades iniciales de población.

5. Significado y Reivindicaciones

El artículo sostiene que su principal importancia radica en cerrar la brecha entre la consistencia física y el realismo ecológico en la modelización fraccionaria. Al asegurar el balance dimensional, el modelo evita los artefactos matemáticos asociados con la fraccionalización ingenua.

Los autores enfatizan que su trabajo proporciona un marco unificado para estudiar cómo la memoria, el retardo, el miedo, el refugio y la cosecha influyen conjuntamente en la persistencia poblacional. Específicamente, afirman haber:

• Demostrado que la memoria fraccionaria puede actuar como un mecanismo estabilizador contra las oscilaciones inducidas por el retardo.
• Proporcionado una clasificación completa de los escenarios de cambio de estabilidad, mostrando que el retardo puede tener un papel dual (desestabilizando estados estables y estabilizando estados inestables).
• Ofrecido perspectivas de gestión, mostrando que la intensidad de la cosecha es un parámetro crítico que puede ajustarse para controlar la invasión del depredador y la coexistencia.

El estudio concluye que trabajos futuros podrían extender este marco para incluir la difusión espacial, perturbaciones estocásticas o estrategias de cosecha óptimas, pero la contribución actual se centra en la formulación rigurosa y el análisis de estabilidad del sistema con retardo dimensionalmente balanceado.