Feasibility as a moving target: Fluctuating species interactions lead to universal power law in equilibrium abundances

Título: A Viabilidade como um Alvo Móvel: Interações de Espécies Flutuantes Levam a uma Lei de Potência Universal nas Abundâncias de Equilíbrio

Autores: Cagatay Eskin, Vu Nguyen e Dervis Can Vural
Instituição: Universidade de Notre Dame, EUA

---

1. Problema e Motivação

A ecologia teórica tradicionalmente equipara a persistência de um ecossistema à estabilidade local de um ponto de equilíbrio fixo (baseado no trabalho seminal de Robert May). No entanto, os autores argumentam que a estabilidade matemática é insuficiente se o ponto de equilíbrio estiver localizado em um ortante negativo (ou seja, se as abundâncias calculadas forem negativas), o que implica extinção biológica.

O problema central abordado é a viabilidade (a condição de que todas as abundâncias de equilíbrio sejam positivas) na presença de flutuações ambientais. Diferente de perturbações que apenas deslocam temporariamente as abundâncias em torno de um ponto fixo, flutuações nas interações entre espécies (ex: mudanças sazonais, ontogenéticas ou climáticas) deslocam o próprio ponto de equilíbrio. O objetivo é entender como essas flutuações afetam a probabilidade de um ecossistema "escapar" da viabilidade (entrar em um estado onde uma ou mais espécies teriam abundância negativa) e como isso escala com a complexidade da comunidade.

2. Metodologia

Modelo Teórico

• Framework: Utilizam o modelo de Lotka-Volterra generalizado (gLV) para descrever a dinâmica da biomassa $n_i(t)$.
• Flutuações: As interações na matriz $A_{ij}$ não são constantes, mas flutuam em torno de um valor médio. Essas flutuações são modeladas como um Processo de Ornstein-Uhlenbeck (OU), que incorpora forças restauradoras (tendência a retornar ao valor canônico) e ruído aleatório.
• Equilíbrio Estocástico: O ponto de equilíbrio $\vec{y}$ torna-se uma variável aleatória $\vec{y} = -[A + \delta A]^{-1}\vec{r}$. Os autores derivam analiticamente a distribuição de probabilidade conjunta $h(\vec{y})$ das abundâncias de equilíbrio.
• Análise de Caudas: Investigam o comportamento das caudas da distribuição de probabilidade para grandes desvios, focando na probabilidade de $y_i$ tornar-se negativo.

Validação Empírica e Simulação

• Dados Reais: Analisaram dois grandes conjuntos de dados:
  1. 33 anos de medições de biomassa no Lago Paul (dados de campo).
  2. Um mesocosmo laboratorial com comunidade de plâncton do Mar Báltico (2.300 dias de condições controladas).
• Teste Estatístico: Aplicaram o framework rigoroso de Clauset et al. para testar a aderência a leis de potência, comparando com distribuições exponenciais e log-normais.
• Simulações Numéricas: Realizaram simulações estocásticas (esquema Euler-Maruyama) em 98 redes ecológicas empíricas (82 redes mutualísticas e 16 redes alimentares) para validar as previsões analíticas de taxas de escape e tempo de persistência.

3. Principais Contribuições e Resultados

A. Lei de Potência Universal para Abundâncias

• Descoberta Teórica: Os autores demonstram que, independentemente da estrutura da rede ou do tamanho da comunidade, flutuações de cauda leve (Gaussianas) nas interações levam a uma distribuição de abundância de equilíbrio de cauda pesada descrita por uma lei de potência:
  $$P(y) \sim \frac{1}{y^\alpha}$$
• Expoente Universal: O expoente teórico é $\alpha = 2$. Este valor é universal e não depende da estrutura de interação, do tamanho da comunidade ($N$) ou da força do ruído.
• Implicação: Isso significa que, mesmo com flutuações pequenas, existe uma probabilidade não desprezível de grandes deslocamentos no equilíbrio, empurrando espécies para a extinção (valores negativos).

B. Evidência Empírica

• A análise de 34 espécies nos dados empíricos revelou que a maioria segue uma lei de potência.
• O expoente mediano observado foi $\alpha \approx 2.56$, situando-se no intervalo de 2 a 3, o que está próximo da previsão teórica de 2.
• A discrepância entre o valor teórico (2) e o empírico (~2.56) é explicada por:
  1. Respostas funcionais não-lineares (tipo Holling), que alteram o expoente para $\alpha = m + 1$ (onde $m$ é o expoente de Holling, tipicamente entre 1 e 3).
  2. Tempos de relaxação lentos de algumas espécies, que "podam" as extremidades da distribuição, aumentando o expoente aparente.

C. Fragilidade e Ruído Crítico

• Perda de Viabilidade: Os autores derivam uma fórmula para o ruído crítico ($\sigma_c$) além do qual a perda de viabilidade ocorre com "certeza quase absoluta".
• Escala com o Tamanho: O ruído crítico escala com o tamanho da comunidade ($N$) como:
  $$\sigma_c(N) \sim N^{-1}$$
• Conclusão: Comunidades maiores são significativamente mais frágeis a flutuações de interação. À medida que $N \to \infty$, qualquer ruído infinitesimal pode levar à inviabilidade do sistema, mesmo que a estabilidade local (autovalores) seja mantida.

D. Métricas de Avaliação de Risco

O estudo propõe duas métricas analíticas mensuráveis:

1. Taxa de Escape Específica ($\lambda_i$): Classifica quais espécies são mais propensas a precipitar a perda de viabilidade.
2. Tempo Médio até a Perda de Viabilidade ($\langle T_{loss} \rangle$): Estima a persistência global do sistema.

• Validação: As previsões analíticas dessas métricas concordaram fortemente com as simulações estocásticas nas 98 redes reais (Coeficiente de Correlação de Concordância > 0.90 para redes alimentares).

4. Significado e Impacto

• Mudança de Paradigma: O trabalho desafia a visão tradicional de que a estabilidade local é o principal determinante da persistência. Ele sugere que a viabilidade (estar em um ortante positivo) é o gargalo crítico, especialmente em sistemas complexos sujeitos a flutuações.
• Eventos "Black Swan": Oferece um mecanismo teórico para colapsos súbitos de ecossistemas que parecem estáveis. Um sistema pode manter a estabilidade dinâmica local, mas o ponto de equilíbrio pode "escapar" para a região de extinção devido a flutuações raras, mas grandes, nas interações.
• Gerenciamento de Conservação: As métricas desenvolvidas permitem identificar quais espécies são "pontos de falha" (mais vulneráveis a levar o sistema à inviabilidade) e quantificar a fragilidade de redes complexas em relação ao ruído ambiental.
• Universalidade: A descoberta de uma lei de potência universal com $\alpha=2$ sugere uma propriedade fundamental de sistemas ecológicos complexos sob flutuações, independente de detalhes específicos da rede.

Em resumo, o artigo estabelece que a persistência de ecossistemas complexos é limitada não apenas pela estabilidade dinâmica, mas pela probabilidade de flutuações nas interações deslocarem o equilíbrio para estados biologicamente inviáveis, um risco que aumenta drasticamente com o tamanho da comunidade.