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

Resumo Técnico: Estabilidade de Alternância de Atraso de Memória e Limiares Ecológicos em um Modelo Predador-Presa de Ordem Fracionária Dimensionalmente Balanceado

1. Definição do Problema e Motivação

O artigo aborda uma falha metodológica crítica em muitos modelos ecológicos de ordem fracionária existentes: a inconsistência dimensional. Em modelos clássicos de ordem inteira, a taxa de variação $dx/dt$ possui dimensões de $[Massa][Tempo]^{-1}$. No entanto, a derivada de Caputo de ordem fracionária $0 < \alpha < 1$, denotada por $D_t^\alpha$, possui inerentemente dimensões de $[Massa][Tempo]^{-\alpha}$. Substituir diretamente a derivada de ordem inteira por uma fracionária sem ajustar os parâmetros do sistema resulta em um descompasso dimensional, tornando o modelo fisicamente sem sentido.

Além disso, os autores argumentam que os sistemas ecológicos são regidos por múltiplos mecanismos interativos frequentemente estudados isoladamente:

• Efeitos de memória: Populações biológicas respondem ao histórico passado, não apenas aos estados atuais.
• Atrasos temporais (delays): O crescimento do predador frequentemente apresenta um atraso em relação ao consumo da presa devido à gestação ou digestão.
• Fatores comportamentais: O medo da predação e o refúgio da presa alteram significamente a dinâmica populacional.
• Gestão: A colheita não linear (ex: pesca) é um mecanismo de controle realista.

O estudo visa formular um modelo predador-presa de ordem fracionária dimensionalmente balanceado que integre esses fatores (medo, refúgio, resposta funcional do tipo Holling III, colheita não linear e atraso temporal) enquanto garante a consistência física.

2. Metodologia

Formulação do Modelo:
Os autores constroem um sistema dimensionalmente balanceado reconstruindo os parâmetros de taxa dependentes do tempo. Em vez de simplesmente substituir $r$ por $r^\alpha$, eles definem coeficientes de taxa fracionários (ex: $r_\alpha = r^\alpha$, $d_\alpha = d^\alpha$, $q_\alpha = q^\alpha$) para garantir que cada termo nas equações diferenciais tenha a dimensão $[Massa][Tempo]^{-\alpha}$. Cuidados especiais são tomados com parâmetros de dimensões mistas (como o coeficiente de autolimitação $\sigma$), redefinindo-os como $\sigma_\alpha = \sigma d^{\alpha-1}$ para manter a homogeneidade.

O sistema resultante (Equação 8) é:
$$\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}$$
onde $H_\alpha$ é a resposta funcional Holling tipo III dimensionalmente balanceada, $\tau$ é o atraso temporal e $\alpha$ é a ordem fracionária.

Estrutura Analítica:

• Bem-postura (Well-posedness): Os autores provam a existência local, unicidade, positividade e limitação última das soluções usando a formulação de integral de Volterra e princípios de comparação fracionária.
• Análise de Equilíbrio: Eles derivam o equilíbrio de apenas presa e o equilíbrio de coexistência. Um número de invasão do predador ($R_0$) é definido para determinar o limiar para a persistência versus extinção do predador.
• Análise de Estabilidade:
  • Sem atraso ($\tau=0$): A estabilidade é analisada usando o critério de Matignon, que exige que os autovalores $\mu$ do Jacobiano satisfaçam $|\arg(\mu)| > \alpha\pi/2$.
  • Com atraso ($\tau>0$): A equação característica é derivada. Os autores resolvem para frequências críticas e valores de atraso críticos ($\tau_j$) onde as raízes cruzam o eixo imaginário. Condições de transversalidade são estabelecidas para determinar a direção da alternância de estabilidade.

Simulação Numérica:
As simulações são realizadas usando um método de preditor-corretor para equações diferenciais de atraso de ordem fracionária de Caputo. Os autores reconstroem explicitamente os coeficientes dimensionalmente balanceados para cada simulação, em vez de usar uma substituição de parâmetros ingênua.

3. Principais Contribuições

1. Consistência Dimensional: O artigo fornece um arcabouço rigoroso para construir modelos ecológicos fracionários que sejam fisicamente consistentes em termos dimensionais, corrigindo a prática comum de substituição direta de derivadas.
2. Classificação Abrangente de Estabilidade: O estudo classifica e verifica numericamente seis cenários distintos de alternância de estabilidade induzida pelo atraso para o equilíbrio de coexistência:
   • Permanece estável
   • Estável-para-instável
   • Estável-instável-estável (Bolha de estabilidade)
   • Permanece instável
   • Instável-para-estável
   • Instável-estável-instável (Bolha de instabilidade)
3. Interação Memória-Atraso: A análise demonstra como a ordem fracionária $\alpha$ (memória) interage com o atraso temporal $\tau$. O atraso crítico $\tau_0(\alpha)$ mostra-se ser uma função de $\alpha$, indicando que os efeitos de memória podem estabilizar ou desestabilizar o sistema dependendo da magnitude do atraso.
4. Limiares de Gestão Ecológica: O artigo vincula a estabilidade matemática à gestão ecológica ao analisar o impacto da colheita não linear ($q$) no número de invasão do predador ($R_0$).

4. Resultados

• Existência e Persistência: O modelo garante soluções positivas e limitadas. Uma condição de limiar $R_0(q, \alpha) > 1$ é estabelecida para a invasão do predador. Os autores provam que o aumento da intensidade de colheita $q$ diminui estritamente $R_0$, podendo levar o predador à extinção.
• Alternância de Estabilidade:
  • Para o sistema sem atraso, a estabilidade depende do traço e do determinante do Jacobiano e da ordem fracionária $\alpha$.
  • Para o sistema com atraso, a introdução de $\tau$ pode induzir bifurcações de Hopf. O estudo confirma que a estabilidade não é monotônica em relação ao atraso; o sistema pode transitar de estável para instável e de volta para estável (ou vice-versa) conforme $\tau$ aumenta.
• Estabilização Induzida pela Memória: Resultados numéricos mostram que, para um atraso fixo que causa oscilações em um modelo de ordem inteira, a redução da ordem fracionária $\alpha$ (aumento da memória) pode suprimir oscilações e estabilizar o equilíbrio de coexistência.
• Sensibilidade: A análise de sensibilidade indica que a dinâmica do sistema é altamente responsiva a variações nas taxas de colheita e nas densidades populacionais iniciais.

5. Significância e Alegações

O artigo alega que sua principal significância reside em preencher a lacuna entre a consistência física e o realismo ecológico na modelagem fracionária. Ao garantir o balanço dimensional, o modelo evita artefatos matemáticos associados à fracionarização ingênua.

Os autores enfatizam que seu trabalho fornece um arcabouço unificado para estudar como memória, atraso, medo, refúgio e colheita influenciam conjuntamente a persistência populacional. Especificamente, eles alegam ter:

• Demonstrado que a memória fracionária pode atuar como um mecanismo estabilizador contra oscilações induzidas pelo atraso.
• Fornecido uma classificação completa dos cenários de alternância de estabilidade, mostrando que o atraso pode ter um papel duplo (desestabilizar estados estáveis e estabilizar estados instáveis).
• Oferecido insights de gestão, mostrando que a intensidade da colheita é um parâmetro crítico que pode ser ajustado para controlar a invasão do predador e a coexistência.

O estudo conclui que trabalhos futuros poderiam estender este arcabouço para incluir difusão espacial, perturbações estocásticas ou estratégias de colheita otimizadas, mas a contribuição atual foca na formulação rigorosa e na análise de estabilidade do sistema com atraso dimensionalmente balanceado.