Dynamical stability of mutualistic meta-networks depends on the type of dispersal

Resumen Técnico: La estabilidad dinámica de las metaredes mutualistas depende del tipo de dispersión

Planteamiento del problema
Si bien las interacciones mutualistas son fundamentales para el funcionamiento de los ecosistemas (por ejemplo, la polinización), el impacto de la estructura espacial y la dispersión sobre la estabilidad de estos sistemas aún no se comprende completamente. Estudios previos sobre metaredes mutualistas han asumido mayoritariamente una dispersión difusiva e independiente de la densidad. Sin embargo, las observaciones empíricas indican que la dispersión suele ser dependiente de la densidad e influenciada por las abundancias de otras especies (codifusión). Por ejemplo, la dispersión de semillas de plantas puede disminuir con la densidad debido al sombreado, mientras que los polinizadores pueden dispersarse preferentemente hacia parches con socios mutualistas abundantes. Los mecanismos específicos mediante los cuales estos comportamientos de dispersión no lineales y dependientes de la densidad, particularmente la codifusión, afectan la estabilidad dinámica de las metaredes mutualistas no estaban claros anteriormente.

Metodología
Los autores desarrollaron un Modelo Generalizado para metaredes mutualistas con el fin de investigar estas dinámicas. Este enfoque combina el realismo ecológico de los modelos basados en ecuaciones diferenciales con la eficiencia computacional de los modelos de matrices aleatorias.

• Estructura de la red: El modelo genera metaredes consistentes en $N_P$ parches de hábitat. Dentro de cada parche, se construye una red mutualista de especies de plantas y animales utilizando el algoritmo de Thébault y Fontaine (2010), variando la riqueza de especies ($S$) entre 20 y 100 y la conectancia ($C$) entre 0.1 y 0.7. Las conexiones espaciales entre parches se modelan como grafos geométricos aleatorios con grados de conectividad variables.
• Dinámica intra-parche: La dinámica dentro de un parche sigue un modelo mutualista (basado en Metz et al., 2026) que incorpora competencia por recursos, beneficios mutualistas saturantes y mortalidad superlineal.
• Modelado de la dispersión: El modelo se extiende a las metaredes incorporando términos de dispersión que dependen de las abundancias de las poblaciones. Los autores distinguen entre:
  • Autodifusión: Dispersión dependiente de la propia abundancia de la especie.
  • Codifusión: Dispersión dependiente de la abundancia de otras especies (intra-gremio o inter-gremio).
  • Elasticidades: El modelo utiliza elasticidades (derivadas logarítmicas) para parametrizar estas dependencias, permitiendo tanto retroalimentaciones positivas como negativas (por ejemplo, $\hat{\kappa} > 0$ implica un aumento de la inmigración con una mayor abundancia del objetivo).
• Análisis de estabilidad: Se realiza un análisis de estabilidad lineal sobre el estado estacionario homogéneo. La matriz Jacobiana de la metared de deriva analíticamente en la forma $J = I \otimes Z - L \otimes D$, donde $Z$ es la Jacobiana intra-parche, $L$ es el Laplaciano espacial y $D$ es la matriz de dispersión. La estabilidad se determina mediante los autovalores de las matrices resultantes; un sistema se considera dinámicamente inestable si cualquier autovalor tiene una parte real $\ge 0$. La "proporción de redes estables" (PSW, por sus siglas en inglés) se utiliza como la métrica principal de estabilidad a través de grandes conjuntos de réplicas.

Contribuciones clave y resultados

1. La dispersión reduce la estabilidad en estados homogéneos: Contrario a algunos hallazgos previos donde la dispersión estabiliza los sistemas, la inclusión de comportamientos de dispersión realistas y dependientes de la densidad disminuye significativamente la estabilidad dinámica de los estados estacionarios homogéneos en las metaredes mutualistas. Esta desestabilización aumenta con el número de hábitats y su conectividad.
2. Reversión de las relaciones de complejidad-estabilidad: En modelos de un solo parche, la estabilidad típicamente aumenta con la riqueza de especies y la conectancia. Los autores encuentran que en las metaredes con codifusión, estas relaciones se alteran:
   • La relación positiva entre estabilidad y riqueza de especies se vuelve negativa.
   • La relación entre estabilidad y conectancia se vuelve unimodal (aumenta y luego disminuye).
3. El papel de la codifusión: El estudio identifica la codifusión (dependencia de la dispersión de las abundancias de otras especies) como el principal motor de la inestabilidad.
   • La codifusión inter-gremio generalmente reduce la estabilidad.
   • La codifusión intra-gremio es responsable de la observada relación unimodal con la conectancia.
   • Los sistemas que presentan únicamente autodifusión (dispersión independiente de la densidad o puramente dependiente de la propia densidad) no exhiben estas inestabilidades y mantienen patrones de estabilidad similares a los modelos de un solo parche.
4. Mecanismo de inestabilidad: Las inestabilidades son impulsadas por bucles de retroalimentación positiva. Por ejemplo, si la abundancia de una especie en un parche objetivo aumenta, esto puede desencadenar un aumento en la inmigración de su socio mutualista (debido a una elasticidad de codifusión positiva), lo que a su vez aumenta aún más la abundancia de la primera especie. Esto amplifica las perturbaciones en lugar de amortiguarlas.
5. Dependencia de la tasa de dispersión: El efecto desestabilizador es altamente sensible a la tasa de recambio de la dispersión. Cuando el recambio de la dispersión es bajo (por ejemplo, cada 1,000 individuos se dispersa), la proporción de redes estables se recupera a niveles comparables con los de los sistemas de un solo parche, incluso con alta riqueza de especies y conectividad.
6. Naturaleza de las inestabilidades: Las inestabilidades observadas están predominantemente asociadas con autovalores reales, lo que sugiere el inicio de la formación de patrones espaciales (inestabilidades de tipo Turing) en lugar de oscilaciones temporales puras. Los autores señalan que dicha formación de patrones puede no ser perjudicial y podría potencialmente promover la coexistencia de especies o la diversidad funcional.

Significancia
El artículo proporciona una comprensión detallada de cómo el tipo de comportamiento de dispersión, específicamente la codifusión dependiente de la densidad, altera fundamentalmente el panorama de estabilidad de las metacomunidades mutualistas. Desafía la suposición de que la dispersión estabiliza universalmente las redes ecológicas, demostrando que la dispersión no lineal motivada biológicamente puede desencadenar inestabilidades impulsadas por la difusión. Al utilizar el Modelado Generalizado, el estudio aísla los mecanismos de retroalimentación específicos (codifusión) responsables de estas dinámicas, ofreciendo una visión más matizada de la relación complejidad-estabilidad en sistemas mutualistas extendidos espacialmente. Los hallazgos sugieren que la estabilidad de las redes mutualistas del mundo real puede depender críticamente de la magnitud del recambio de la dispersión y de la naturaleza específica de las interacciones entre especies durante la dispersión.