A Mathematical Model of Dengue Transmission Incorporating Hospital Capacity and Threshold-Based Fogging Interventions

Resumen Técnico: Un Modelo Matemático de Transmisión del Dengue que Incorpora la Capacidad Hospitalaria e Intervenciones de Nebulización Basadas en Umbrales

Planteamiento del Problema
El dengue sigue siendo un desafío persistente para la salud pública en las regiones tropicales, con brotes recurrentes que sugieren que las estrategias de intervención actuales son insuficientes. Los modelos matemáticos existentes suelen basarse en supuestos idealizados que no logran capturar restricciones operativas críticas: específicamente, suelen asumir una capacidad hospitalaria ilimitada y la aplicación continua de la nebulización (fumigación). En la realidad, los recursos hospitalarios son finitos y la nebulización rara vez se aplica de forma continua debido a restricciones financieras, logísticas y ambientales. En su lugar, la nebulización se implementa frecuentemente de forma reactiva, activándose solo cuando los indicadores de vigilancia (como las infecciones reportadas) cruzan un umbral específico. Esta discrepancia entre los modelos teóricos y las restricciones prácticas motiva el desarrollo de un modelo que incorpore explícitamente una capacidad hospitalaria finita y una nebulización no suave activada por umbrales.

Metodología
Los autores desarrollan un modelo de ecuaciones diferenciales ordinarias (EDO) no suaves que describe la dinámica de transmisión entre las poblaciones humana y de mosquitos. La población humana se divide en compartimentos de susceptibles ($S$), infectados asintomáticos ($A$), infectados sintomáticos ($I$) y recuperados ($R$), mientras que la población de mosquitos consta de compartimentos de susceptibles ($U$) e infectados ($V$).

El modelo introduce tres regímenes operativos distintos gobernados por dos umbrales críticos:

1. Régimen de Baja Infección ($I < kC$): No se aplica nebulización ($h(I)=0$). Todos los individuos sintomáticos son hospitalizados y la recuperación ocurre a la tasa hospitalaria $\gamma_1$.
2. Régimen de Infección Moderada ($kC \le I < C$): Se activa la nebulización ($h(I)=\eta$) una vez que el número de casos hospitalizados alcanza una fracción $k$ de la capacidad total $C$. La hospitalización permanece completa ($\gamma_1$), pero la mortalidad de los mosquitos aumenta debido a la nebulización.
3. Régimen de Sobrecarga ($I \ge C$): La nebulización permanece activa, pero la capacidad hospitalaria se excede. Los individuos sintomáticos más allá de la capacidad $C$ reciben tratamiento en el hogar con una tasa de recuperación menor $\gamma_2$ ($\gamma_2 < \gamma_1$). Crucialmente, estos individuos no tratados permanecen infecciosos y contribuyen a la transmisión, modelado por una tasa de infección aumentada $\beta_{m2}$ para la fracción no hospitalizada.

El sistema es tratado como un sistema dinámico suave por partes. Los autores realizan una investigación analítica rigurosa sobre la existencia y la estabilidad local de los equilibrios libres de enfermedad y endémicos dentro de cada régimen. Además, emplean técnicas de continuación numérica utilizando el paquete de software COCO para construir diagramas de bifurcación, analizando el comportamiento del sistema a medida que varían parámetros clave (tasa de transmisión $\beta_h$, tasa de nebulización $\eta$ y umbral de activación $k$).

Contribuciones Clave y Resultados

• Dinámica No Suave y Bifurcaciones: El estudio establece que la interacción entre la capacidad hospitalaria y la nebulización basada en umbrales crea un paisaje dinámico complejo. El análisis identifica bifurcaciones de equilibrio de frontera que ocurren precisamente en los umbrales de conmutación ($I=kC$y$I=C$). Estas bifurcaciones marcan las transiciones entre los tres regímenes operativos.
• Bifurcación de Hopf y Brotes Periódicos: Cuando la tasa de transmisión es lo suficientemente alta como para empujar las infecciones más allá de la capacidad hospitalaria ($I > C$), el sistema experimenta una bifurcación de Hopf. Esto conduce a soluciones periódicas sostenidas, lo que implica que los brotes pueden no establecerse en un estado endémico estable, sino que pueden recurrir en oleadas con picos de infección elevados.
• Bifurcación de Plegado (Fold) y Complejidad del Equilibrio: Cerca del umbral epidémico, se identifica una bifurcación de plegado, que produce equilibrios endémicos inestables adicionales incluso cuando el número básico de reproducción ($R_0$) es menor que 1. Esto indica una estructura de equilibrio más intrincada de lo que sugieren los modelos suaves estándar, donde la enfermedad podría persistir en estados inestables antes de desaparecer.
• Eficacia de Control No Monotónica: La continuación numérica revela que la relación entre los parámetros de control y los resultados de salud pública es no monotónica.
  • Tasa de Nebulización ($\eta$) y Umbral ($k$): Existen valores óptimos locales tanto para la intensidad de la nebulización como para el umbral de activación que minimizan los p la de infección máxima ($I_{Peak}$) y la fracción de casos asintomáticos ($AFrac$).
  • Hallazgos Contraintuitivos: La activación demasiado temprana (bajo $k$) o la nebulización excesivamente intensa no necesariamente producen beneficios epidemiológicos proporcionales. La activación temprana mantiene la nebulización activa durante períodos más largos, pero debido a los retrasos estructurales en la dinámica humana y de los mosquitos, esto puede provocar un desajuste con el ciclo epidémico, causando la acumulación de susceptibles y un mayor rebote una vez que el umbral ya no se cumple. Del mismo modo, la nebulización de alta intensidad suprime fuertemente a los mosquitos pero brevemente, lo que conduce a una rápida recuperación de la población una vez que la intervención cesa, pudiendo causar un resurgimiento más agudo.

Significancia
El artículo argumenta que la capacidad hospitalaria, la intensidad de la nebulización y el umbral de activación no son meros detalles operativos, sino impulsores fundamentales de la dinámica del dengue. Los hallazgos sugieren que las políticas de control efectivas requieren la calibración conjunta de estos factores en lugar de tratar la intervención temprana o intensiva como inherentemente óptima. El modelo proporciona orientación cuantitativa para diseñar políticas de control reactivo del dengue, destacando que una nebulización mal cronometrada o retrasada puede prolongar la transmisión o desencadenar resurgimientos más grandes. Al capturar la naturaleza discontinua de las intervenciones del mundo real, el modelo ofrece un marco más realista para comprender cómo las restricciones de recursos y las políticas reactivas moldean la persistencia y la recurrencia de los brotes de dengue.

Direcciones Futuras
Los autores señalan que trabajos futuros podrían incluir una clasificación formal de la criticidad de la bifurcación de Hopf, el análisis de modo deslizante de las variedades de conmutación (switching manifolds), e la incorporación de forzamiento estacional, efectos estocásticos y heterogeneidad espacial. Adicionalmente, se propone formular la estrategia de intervención como un problema de control óptimo calibrado con datos reales para mejorar la utilidad del modelo como herramienta de apoyo a la decisión.