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

Technische Samenvatting: Een wiskundig model van dengue-transmissie met incorporatie van ziekenhuiscapaciteit en drempelgebaseerde vernevelingsinterventies

Probleemstelling
Dengue koorts blijft een hardnekkige uitdaging voor de volksgezondheid in tropische regio's, waarbij terugkerende uitbraken erop wijzen dat de huidige interventiestrategieën onvoldoende zijn. Bestaande wiskundige modellen leunen vaak op geïdealiseerde aannames die er niet in slagen kritieke operationele beperkingen te vatten: specifiek gaan ze doorgaans uit van onbeperkte ziekenhuiscapaciteit en de continue toepassing van verneveling (fogging). In de werkelijkheid zijn ziekenhuisbronnen eindig en wordt verneveling zelden continu toegepast vanwege financiële, logistieke en milieutechnische beperkingen. In plaats daarvan wordt verneveling vaak reactief uitgevoerd, getriggerd zodra surveillance-indicatoren (zoals gerapporteerde infecties) een specifieke drempel overschrijden. Dit verschil tussen theoretische modellen en praktische beperkingen motiveert de ontwikkeling van een model dat expliciet eindige ziekenhuiscapaciteit en drempelgestuurde, niet-gladde vernevelingsinterventies incorporeert.

Methodologie
De auteurs ontwikkelen een niet-glad gewone differentiaalvergelijking (ODE) model dat de transmissiedynamiek tussen menselijke en muggenpopulaties beschrijft. De menselijke populatie wordt onderverdeeld in compartimenten voor vatbare ($S$), asymptomatisch geïnfecteerde ($A$), symptomatisch geïnfecteerde ($I$) en herstelde ($R$) individuen, terwijl de muggenpopulatie bestaat uit vatbare ($U$) en geïnfecteerde ($V$) compartimenten.

Het model introduceert drie verschillende werkingsregimes die worden beheerst door twee kritieke drempels:

1. Regime met lage infectie ($I < kC$): Er vindt geen verneveling plaats ($h(I)=0$). Alle symptomatische individuen worden opgenomen in het ziekenhuis en herstel vindt plaats tegen de ziekenhuisrate $\gamma_1$.
2. Regime met matige infectie ($kC \le I < C$): Verneveling wordt geactiveerd ($h(I)=\eta$) zodaal het aantal opgenomen gevallen een fractie $k$ van de totale capaciteit $C$ bereikt. De hospitalisatie blijft volledig ($\gamma_1$), maar de mortaliteit van muggen neemt toe door de verneveling.
3. Overbelastingsregime ($I \ge C$): Verneveling blijft actief, maar de ziekenhuiscapaciteit wordt overschreden. Symptomatische individuen die de capaciteit $C$ overschrijden, ontvangen thuiszorg met een lagere herstelrate $\gamma_2$ ($\gamma_2 < \gamma_1$). Cruciaal is dat deze niet-behandelde individuen infectieus blijven en bijdragen aan de transmissie, wat wordt gemodelleerd door een verhoogde infectierate $\beta_{m2}$ voor het niet-gehospitaliseerde deel.

Het systeem wordt behandeld als een stuksgewijs glad dynamisch systeem. De auteurs voeren een rigoureuze analytische studie uit naar de existentie en lokale stabiliteit van ziektevrije en endemische evenwichten binnen elk regime. Bovendien maken zij gebruik van numerieke continuatietechnieken met de COCO softwarepackage om bifurcatiediagrammen te construeren, waarbij zij het gedrag van het systeem analyseren naarmate sleutelparameters (transmissierate $\beta_h$, vernevelingsrate $\eta$, en activatiedrempel $k$) variëren.

Belangrijkste bijdragen en resultaten

• Niet-gladde dynamiek en bifurcaties: De studie stelt vast dat de interactie tussen ziekenhuiscapaciteit en drempelgestuurde verneveling een complexe dynamische landschap creëert. De analyse identificeert boundary-equilibrium bifurcaties die precies plaatsvinden bij de schakelingsdrempels ($I=kC$ en $I=C$). Deze bifurcaties markeren de overgangen tussen de drie werkingsregimes.
• Hopf-bifurcatie en periodieke uitbraken: Wanneer de transmissierate voldoende hoog is om infecties voorbij de ziekenhuiscapaciteit te duwen ($I > C$), ondergaat het systeem een Hopf-bifurcatie. Dit leidt tot aanhoudende periodieke oplossingen, wat impliceert dat uitbraken niet noodzakelijkerwijs stabiliseren in een endemische toestand, maar eerder in golven met hoge infectiepieken terugkeren.
• Fold-bifurcatie en evenwichtcomplexiteit: Nabij de epidemische drempel wordt een fold-bifurcatie geïdentificeerd, die extra onstabiele endemische evenwichten produceert, zelfs wanneer het basisreproductiegetal ($R_0$) kleiner is dan 1. Dit duidt op een meer complexe evenwichtsstructuur dan standaard gladde modellen suggereren, waarbij de ziekte potentieel in onstabiele toestanden kan voortbestaan voordat deze uitsterft.
• Niet-monotone effectiviteit van controle: Numerieke continuatie onthult dat de relatie tussen controleparameters en volksgezondheidsuitkomsten niet-monotoon is.
  • Vernevelingsrate ($\eta$) en Drempel ($k$): Er bestaan lokaal optimale waarden voor zowel de intensiteit van de verneveling als de activatiedrempel om de piekinfecties ($I_{Peak}$) en het deel asymptomatische gevallen ($AFrac$) te minimaliseren.
  • Contraintuïtieve bevindingen: Te vroege verneveling (lage $k$) of excessief intense verneveling levert niet noodzakelijkerwijs proportionele epidemiologische voordelen op. Vroege activatie houdt de verneveling langer actief, maar door de structurele vertragingen in de menselijke en muggen-dynamica kan dit leiden tot een mismatch met de epidemische cyclus, wat resulteert in de accumulatie van vatbare personen en een grotere terugslag zodra de drempel niet langer wordt behaald. Eveneens onderdrukt intensieve verneveling de muggenpopulatie sterk maar kortstondig, wat kan leiden tot een snelle populatieherstel zod zodra de interventie stopt, wat mogelijk een scherpere resurgence veroorzaakt.

Significantie
Het artikel betoogt dat ziekenhuiscapaciteit, vernevelingsintensiteit en de activatiedrempel niet louter operationele details zijn, maar fundamentele drijvers van de dengue-dynamiek. De bevindingen suggereren dat effectieve controlebeleid vereist dat deze factoren gezamenlijk worden gekalibreerd, in plaats van te ervan uit te gaan dat vroege of intensieve interventie inherent optimaal is. Het model biedt kwantitatieve begeleiding voor het ontwerpen van reactieve dengue-controlebeleid, waarbij wordt benadrukt dat slecht getimede of vertraagde verneveling de transmissie kan verlengen of grotere resurgenties kan triggeren. Door de discontinue aard van real-world interventies te vangen, biedt het model een realistischer kader voor het begrijpen van hoe middelenrestricties en reactieve beleidsmaatregelen het voortbestaan en de terugkeer van dengue-uitbraken vormgeven.

Toekomstige richtingen
De auteurs merken op dat toekomstig werk een formele classificatie van de Hopf-bifurcatie kritikaliteit, sliding-mode analyse van de schakelmanifolds, en de incorporatie van seizoensgebonden forcering, stochastische effecten en ruimtelijke heterogeniteit zou kunnen omvatten. Daarnaast wordt voorgesteld om de interventiestrategie te formuleren als een optimaal controleprobleem dat wordt gekalibreerd tegen reële data, om de bruikbaarheid van het model als beslissingsondersteunend instrument te vergroten.