Persistent oscillations and backward bifurcation in a malaria model with varying human and mosquito populations: implications for control

• Published in: Journal of mathematical biology Vol. 70; no. 7; pp. 1581 - 1622
• Main Authors: Ngonghala, Calistus N., Teboh-Ewungkem, Miranda I., Ngwa, Gideon A.
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer Berlin Heidelberg 01.06.2015; Springer Nature B.V
• Subjects: Animals; Applications of Mathematics; Computer Simulation; Culicidae - parasitology; Culicidae - physiology; Female; Host-Parasite Interactions; Humans; Insect Control; Logistic Models; Malaria - prevention & control; Malaria - transmission; Male; Mathematical and Computational Biology; Mathematical Concepts; Mathematics; Mathematics and Statistics; Models, Biological; Plasmodium - pathogenicity; Population Dynamics; Reproduction; 34A; 34C; Sensitivity analysis; 37G; Backward bifurcation; Oscillatory dynamics; Mosquito demography; Malaria transmission and control; 92B; 92D; Reproduction numbers
• ISSN: 0303-6812; 1432-1416
• DOI: 10.1007/s00285-014-0804-9

We derive and study a deterministic compartmental model for malaria transmission with varying human and mosquito populations. Our model considers disease-related deaths, asymptomatic immune humans who are also infectious, as well as mosquito demography, reproduction and feeding habits. Analysis of the model reveals the existence of a backward bifurcation and persistent limit cycles whose period and size is determined by two threshold parameters: the vectorial basic reproduction number R m , and the disease basic reproduction number R 0 , whose size can be reduced by reducing R m . We conclude that malaria dynamics are indeed oscillatory when the methodology of explicitly incorporating the mosquito’s demography, feeding and reproductive patterns is considered in modeling the mosquito population dynamics. A sensitivity analysis reveals important control parameters that can affect the magnitudes of R m and R 0 , threshold quantities to be taken into consideration when designing control strategies. Both R m and the intrinsic period of oscillation are shown to be highly sensitive to the mosquito’s birth constant λ m and the mosquito’s feeding success probability p w . Control of λ m can be achieved by spraying, eliminating breeding sites or moving them away from human habitats, while p w can be controlled via the use of mosquito repellant and insecticide-treated bed-nets. The disease threshold parameter R 0 is shown to be highly sensitive to p w , and the intrinsic period of oscillation is also sensitive to the rate at which reproducing mosquitoes return to breeding sites. A global sensitivity and uncertainty analysis reveals that the ability of the mosquito to reproduce and uncertainties in the estimations of the rates at which exposed humans become infectious and infectious humans recover from malaria are critical in generating uncertainties in the disease classes.