- Open Access
Modelling the impact of vector control interventions on Anopheles gambiae population dynamics
Parasites & Vectorsvolume 4, Article number: 153 (2011)
Intensive anti-malaria campaigns targeting the Anopheles population have demonstrated substantial reductions in adult mosquito density. Understanding the population dynamics of Anopheles mosquitoes throughout their whole lifecycle is important to assess the likely impact of vector control interventions alone and in combination as well as to aid the design of novel interventions.
An ecological model of Anopheles gambiae sensu lato populations incorporating a rainfall-dependent carrying capacity and density-dependent regulation of mosquito larvae in breeding sites is developed. The model is fitted to adult mosquito catch and rainfall data from 8 villages in the Garki District of Nigeria (the 'Garki Project') using Bayesian Markov Chain Monte Carlo methods and prior estimates of parameters derived from the literature. The model is used to compare the impact of vector control interventions directed against adult mosquito stages - long-lasting insecticide treated nets (LLIN), indoor residual spraying (IRS) - and directed against aquatic mosquito stages, alone and in combination on adult mosquito density.
A model in which density-dependent regulation occurs in the larval stages via a linear association between larval density and larval death rates provided a good fit to seasonal adult mosquito catches. The effective mosquito reproduction number in the presence of density-dependent regulation is dependent on seasonal rainfall patterns and peaks at the start of the rainy season. In addition to killing adult mosquitoes during the extrinsic incubation period, LLINs and IRS also result in less eggs being oviposited in breeding sites leading to further reductions in adult mosquito density. Combining interventions such as the application of larvicidal or pupacidal agents that target the aquatic stages of the mosquito lifecycle with LLINs or IRS can lead to substantial reductions in adult mosquito density.
Density-dependent regulation of anopheline larvae in breeding sites ensures robust, stable mosquito populations that can persist in the face of intensive vector control interventions. Selecting combinations of interventions that target different stages in the vector's lifecycle will result in maximum reductions in mosquito density.
Intensive anti-malaria campaigns targeting local mosquito populations with either Indoor Residual Spraying (IRS) or Long-lasting Insecticide Treated Nets (LLINs) in a number of endemic areas, particularly in Africa, have resulted in 10-50 fold reductions in the adult mosquito population and have been associated with similar declines in the prevalence of infection and incidence of disease in people living in these areas [1–4]. Understanding the determinants of such declines is crucial in the context of malaria control and elimination. Consideration of adult mosquitoes alone (the tenet of simple malaria models) is insufficient to fully understand the resulting changes in mosquito density. Instead, the entire mosquito lifecycle, including larvae, pupae and adults, and mosquito behaviour while feeding, resting and ovipositing needs to be considered. Indeed, it has recently been argued that a comprehensive understanding of vector ecology is a prerequisite for malaria elimination .
The gonotrophic cycle of a female Anopheles mosquito begins with a blood meal taken from a human or animal host. The mosquito will then rest while the blood meal is digested and the eggs in the ovaries mature. After resting and egg maturation, the mosquito is classified as gravid and begins searching for a suitable breeding site where it will lay 80-100 eggs [6, 7]. The mosquito will then search for another blood meal and repeat the gonotrophic cycle. Larvae emerge from eggs and feed on bacteria, yeasts, protozoa and particulate organic matter in the water. Growth takes place via a series of moults through four morphologically distinct larval instars. Fourth instar larvae moult to become non-feeding pupae which develop into adult mosquitoes. The duration of the larval period depends mainly on temperature and, in tropical areas, lasts 7-15 days [8–10].
For a stable mosquito population to exist in a fixed environment, the population must be regulated. In ecological models this is usually assumed to be via density-dependent regulation or by a limited environmental carrying capacity. Both mechanisms lead to a limitation in the number of mosquitoes that a particular environment can sustain. In aquatic stage mosquito ecology, carrying capacity describes how many mosquito larvae/pupae an environment can support, whereas density dependence describes the intra-specific competition between larvae for food and resources, resulting in increased mortality and extended developmental times for high larval densities. The net result of increasing density-dependent competition is equivalent to decreasing the carrying capacity; however when considering interventions directed at different stages of the mosquito lifecycle, care must be taken to differentiate between these two concepts. Environmental management, e.g. filling in breeding sites, will directly reduce the carrying capacity without involving inter- or intra-specific interactions, whereas IRS and LLINs directed against adult mosquitoes will cause a reduction in oviposition and hence reduce density-dependent competition between larvae in breeding sites resulting in decreased larval mortality.
Two qualitatively different forms of density-dependent competition in the larval stages of insects have been described: contest competition where the number of larvae surviving in a breeding site reaches a limit as initial egg density increases; and scramble competition where the number of surviving larvae actually decreases at higher initial egg densities . Experimental evidence suggests that density-dependent regulation of A. gambiae larvae follows a model of contest competition [12–14], although there are likely to be factors contributing to scramble competition, for example, cannibalism of early instar larvae by late instar larvae  and an increased population of predators as a result of high larval densities. Several mechanisms contribute to density-dependent regulation of mosquito larvae in breeding sites: (i) increased larval mortality at high densities; (ii) increased developmental times at higher densities; and (iii) reduced larval size, leading to smaller emerging adults at higher densities [13, 14]. Although mosquito populations will be regulated by all of these factors [13, 16], in this paper we focus on density-dependent larval mortality.
Vector control interventions can be targeted at several stages of the mosquito's lifecycle. Successful interventions will result in a dramatic drop in mosquito density, although interventions directed against adult mosquitoes will have the added benefit of reducing the probability that an infected mosquito survives sporogony to become infectious. Understanding mosquito population dynamics will aid in the selection of optimal packages of vector control interventions. The impact of vector control interventions on malaria transmission and adult mosquito dynamics has been widely investigated using mathematical models [17–25]. In addition a number of studies include models of the full mosquito lifecycle including eggs, larvae and adults [26–29]. In contrast to these studies we focus on the interaction between vector control interventions and the adult and aquatic stages of the anopheline vector, rather than malaria transmission between humans. Using insights from a number of entomological studies, a simple model of the A. gambiae sensu lato lifecycle, accounting for both the aquatic and adult stages is developed that is suitable for incorporation into existing mathematical models of malaria transmission. The model explores the comparative impact of LLINs, IRS, larvicides and pupacides, and the benefits that can be achieved by combining these interventions. The ultimate goal is to understand the combined impact of a suite of interventions in malaria endemic settings in Sub-Saharan Africa , with particular reference to those in which A. gambiae s.l. is the prevailing species.
Mosquito population dynamics
The aquatic part of the mosquito lifecycle consists of an egg stage, four larval instar stages and a pupal stage (the pre-imaginal stages). To reduce the complexity of the model in a biologically sensible way, the eggs and the first two larval instars are grouped into the 'early larval instar stage' E. The third and fourth larval instars are grouped together as 'late larval instars' L. The pupal stage is denoted by P and the female adult mosquito stage by M. On average each female mosquito will lay β eggs per day, which will hatch into early instar larvae E. These larvae will undergo density-dependent daily mortality at a rate μ E (E,L) which is assumed to depend on the density of both 'early' and 'late' instars. Larvae surviving the developmental period of d E days (where the reciprocal of d E is the rate of progression to the next stage) will develop into late instars L. These larvae will undergo density-dependent daily mortality at a rate μ L (E,L) during the d L days of development. Instars that survive development will become pupae at a rate given by the reciprocal of d L which are subjected to density independent mortality at a constant rate μ P throughout their d P days of development. It is assumed that half of all emerging adult mosquitoes are female . Male adult mosquitoes are ignored, only assuming there are enough males for successful mating with females. Emerging adult female mosquitoes will search for a blood meal and begin their gonotrophic cycle. There is evidence of gonotrophic discordance in A. gambiae whereby more than one blood meal may be required to produce a batch of eggs, but for simplicity we assume gonotrophic concordance in our model. Although there is evidence that adult mosquitoes senesce [32, 33], for the sake of tractability, we make the simplifying assumption that mosquitoes undergo constant daily mortality at rate μ M . The mosquito lifecycle can thus be described by the following set of continuous ordinary differential equations,
Since the model is formulated as a set of ordinary differential equations with constant rates of progression between stages, the duration of each developmental stage will be exponentially distributed. However the duration of larval and pupal development cannot be exponentially distributed as this would allow some infinitesimally short development times. This problem is overcome by stratifying the aquatic stages into the three compartments E, L and P, which results in an approximately gamma-shaped distribution for the developmental times from egg to pupa . An alternative approach is to use the lumped age-class technique described by Hancock and Godfray  where a fixed time is spent in each developmental stage.
A female A. gambiae mosquito will oviposit 80-100 eggs [6, 7] per gonotrophic cycle which lasts approximately 3 days . Although each individual mosquito oviposits periodically, a population of mosquitoes will lay eggs at an approximately continuous rate. The daily oviposition rate is estimated as the expected number of eggs oviposited in a mosquito's lifetime divided by the expected lifetime.
If a female adult mosquito experiences constant mortality μ M throughout her lifetime and oviposits every δ days then the number of eggs laid per day can be estimated as
Where εmax is the maximum number of eggs per oviposition, and βmax is the maximum rate at which mosquitoes oviposit eggs. In practice ε ≤ εmax and β ≤ βmax as some eggs may be oviposited in unsuitable breeding sites and some eggs may fail to hatch.
Mosquito growth rate and reproduction number
If a female mosquito lays ε viable eggs per oviposition cycle, then over her lifetime she can expect to oviposit approximately eggs where is the proportion of mosquitoes surviving one gonotrophic cycle , is the proportion surviving two gonotrophic cycles, etc. At low larval densities when mortality rates of the stages are at their lowest, background μ0 values, only a fraction
of these eggs will survive to become imagos. The mosquito basic reproduction number R0 is defined to be the expected number of female adult mosquitoes produced per mosquito during her reproductive lifespan in the absence of regulatory constraints. In the absence of density-dependent larval death it is given by,
This expression holds only in the absence of density-dependent larval mortality. In reality density-dependent effects will always be present so the actual value of Reff (the effective reproduction number) will be substantially smaller. We anticipate that more realistic estimates for Reff can be obtained using seasonal rainfall patterns and density dependence in mosquito breeding sites. Conditional on the density-regulated larval death rates μ E (E,L) and μ L (E,L) at time t, the mosquito reproduction number Reff(t) at time t is given by,
Notice that in this expression we are assuming that the per-female mosquito fecundity and mortality rates are density independent. The value of Reff(t) will decrease with increasing density of early and late aquatic stages. When Reff(t) becomes equal to 1, each female mosquito replaces herself and the population is at a stable equilibrium.
Density-dependent larval mortality
In a study of larval development under semi-field conditions in Tanzania by Njunwa , the proportion of first instar larvae introduced into artificial breeding sites that survived to emerge as adult mosquitoes was measured. Two versions of the experiment were carried out. In the first, first instar larvae were placed in the breeding site in one batch so that all larvae were the same age throughout their development. In the second, batches of first instar larvae were introduced at staggered times so that the larval population would have a varying age profile. To assess the relationship between density of instars and the daily mortality rate linear, quadratic and logistic functions were fitted to the data (Figure 1). The mortality of larvae inserted into breeding sites in one batch was well described by a logistic function, whereas the mortality of larvae inserted at staggered times was well described by a linear function although a quadratic function gave a slightly better fit (Supplementary Information in Additional file 1). A linear relationship between density of instars and mortality was selected for the model both for the sake of parsimony and because the staggered input of larvae recreates more closely the structure of larval populations observed under natural conditions [36, 37].
Mosquito aquatic stages feed and compete during the larval instar stages, and undergo morphological development but do not feed during the pupal stage. As such it is assumed that density-dependent regulation occurs at the larval instar stages only. For a breeding site containing E(t) early instar larvae and L(t) late instar larvae at time t, the density-dependent death rates μ E (E(t),L(t)) and μ L (E(t),L(t)) of the early and late instar larvae, respectively, at time t are given by,
where and are the death rates at very low densities, K(t) is the environmental carrying capacity at time t, and γ is a factor allowing for the differential effects of density-dependent mortality on late instar larvae compared to early instar larvae. The assumption that the larval death rate is linearly proportional to the number of larvae is equivalent to the assumption made in a more complex model by Depinay et al.  where larval death rate is proportional to the total larval biomass.
Parameters for the model (see Table 1) are estimated by fitting to adult mosquito catch and meteorological data from 8 villages in the Garki District of Nigeria. The Garki Project, sponsored by the World Health Organization, was undertaken in the 1970's to investigate the feasibility of eliminating malaria from an area of intense transmission . Although malaria was not eliminated during the campaign, detailed entomological, meteorological and parasitological data were collected.
Measurements of daily rainfall were available for 8 villages. Three models for calculating the environmental carrying capacity K(t) were tested:
Carrying capacity proportional to the mean rainfall during the past τ days(8)
ii) Carrying capacity proportional to the linearly weighted mean rainfall during the past τ days
iii) Carrying capacity proportional to past rainfall weighted by an exponential distribution with mean 2τ
where rain(t) is the daily rainfall and λ is a fitted scaling factor unique for each village. Similar rainfall models have been used in the past to describe mosquito dynamics . During extended dry periods each of these models predicts a zero carrying capacity with the consequence that no mosquito population can be sustained. This agrees with data from the Garki Project where no mosquitoes were detected during the dry season (Figure 2). However, in order for re-emergence of the mosquito population at the start of the rainy season there must be some small population of mosquitoes during the dry season . When no mosquitoes are caught in pyrethrum spray catches, we assume that a mosquito could have been caught with probability 0.2. The carrying capacity during the dry season is thus set to a level capable of supporting a small mosquito population detectable 20% of the time.
The mean air temperature during the observational period of the Garki project was 24°C, fluctuating between a minimum temperature of 18°C and a maximum temperature of 30°C. Whilst daily and seasonal fluctuations in temperature are important determinants of the transmission cycle, affecting larval development times, the duration of sporogony for malaria and survival of the adult mosquito , temperature data in the Garki project was not consistently recorded over time and therefore we assumed a single fixed value. Larval development takes place in water, the mean temperature of which is usually higher than air temperature. In a series of experiments carried out under field conditions, Paaijmans et al.  investigated the relationship between air temperature and water temperature in typical larval breeding sites. Based on the results of these experiments, a mean air temperature of 24°C corresponds to a mean water temperature of 28°C. Priors for the parameters were based on these temperatures.
Bayesian methods were used to fit the model to field-based adult mosquito collections (conducted to estimate indoor-resting density) and rainfall data from 8 villages in the Garki Project. Prior estimates of the model parameters are taken from the entomological literature (see Table 1). Laboratory-based studies by Bayoh and Lindsay  on the relationship between water temperature and A. gambiae sensu stricto larval and pupal developmental times and survival rates were used to provide prior estimates for the duration of each stage and the death rates at low larval densities. Several quantitative studies on the sampling of populations of Anopheles larvae in natural breeding sites were identified: Service [36, 42], Njunwa  and Munga et al.  on A. gambiae s.l. larvae and Mwangangi et al.  on A. arabiensis larvae. These studies compiled sampled data into larval life tables displaying the size of the larval population stratified by instar and pupae. Data from Service  on larval populations sampled from marshes and borrow pits were used to obtain prior estimates of the daily oviposition rate of adult female mosquitoes and larval death rates at naturally occurring densities.
The data consisted of observations of adult mosquitoes from pyrethrum spray collections in houses in the 8 villages. Spray collections consisted of a mixture of A. gambiae s.s. and A. arabiensis. Data were aggregated to the village level for fitting and a negative binomial likelihood was chosen to capture the overdispersion in the mosquito catch data. Let N(t,i) denote the number of adult mosquitoes collected on day t in village i, and M(t,i) the number of adult mosquitoes predicted by the model. The likelihood for the parameter vector given the data D can be written as
where r is a parameter of the negative binomial distribution inversely measuring the degree of overdispersion in the mosquito catch data. Parameter estimates were obtained as the medians of the posterior distribution sampled using Markov Chain Monte Carlo (MCMC) methods. The chains were run with a burn-in of 1000 iterations (judged to have converged by visual examination of the parameter chains) which were discarded. The subsequent 1,000,000 iterations were used to estimate the posterior distribution for parameters. The best estimate of the mosquito density over time was obtained by simulating the model using the median parameter set from the posterior distribution. The 95% credible intervals around this output were obtained by repeatedly sampling parameter sets from the joint posterior distribution obtained from the final MCMC chain. Model fit was assessed visually and for different models compared using the predicted posterior probability of each model (see Supplementary Information).
We consider the impact on adult mosquito density of two anti-malarial interventions directed against adult mosquitoes (LLINs and IRS) and two interventions targeting the aquatic stages (larvicide and pupacide). The following sections describe how these are incorporated into the mosquito dynamics model.
Insecticide treated nets and indoor residual spraying
An existing model [30, 44] of the effects of LLINs and IRS on the Anopheles vector is incorporated into our model (see Supplementary Information). We assume the vector population has the behavioural traits of the indoor and night-time biting A. gambiae s.s. mosquito. LLINs are assumed to have three effects on the adult mosquito population: (i) directly killing mosquitoes that land on the nets; (ii) repelling and possibly diverting mosquitoes to an animal blood host due to either insecticide irritation or the physical barrier of the net; and (iii) lengthening the duration of the gonotrophic cycle leading to a reduced oviposition rate. IRS was again assumed to have three effects on the adult mosquito population: (i) killing mosquitoes resting on walls after feeding; (ii) repelling mosquitoes from houses; and (iii) increasing the duration of the gonotrophic cycle. It was assumed that houses were sprayed with the pyrethroid lambdacyhalothrin.
When a larvicide such as Bacillus thuringiensis var. israelensis (BTI) is applied to larval breeding sites, mosquito larvae and pupae experience increased mortality [45–47]. It is assumed that the deployment of a larvicide increases the larval and pupal background death rates , and μ P by a factor of υ, so that the number of surviving larvae is decreased. In order to mimic an 88% reduction in the number of Anopheles larvae as observed by Kroeger et al. , a value of υ = 55.2 is chosen. An 88% reduction in the number of observable larvae corresponds to a 99% reduction in the number of pupae emerging as adults, as the observed larvae are young and unlikely to survive until pupation. If a fraction c of breeding sites has been treated with larvicidal agents, then the mosquito population can be described by the following model,
Here EBTI, LBTI and PBTI represent the early instars, late instars and pupae respectively that reside in breeding sites treated with BTI. It is assumed for tractability that larvae are distributed evenly across viable breeding sites as the purpose of this model is to obtain average behaviour for a given location which can subsequently be expanded in more detailed models. This assumption is unlikely to hold in practice as there may be a great deal of heterogeneity in breeding site preference [48–50] although the spatial scales of such heterogeneity remain poorly understood. Gu and Novak  have demonstrated how heterogeneity in breeding sites can be exploited if the most productive sites can be identified and targeted with larvicide, although concerns have been raised about the difficulty in identifying productive breeding sites under operational conditions .
Pupacides preferentially target the pupal stage of the aquatic lifecycle of the Anopheles mosquito. Killing immature stages at the pupal stage has the additional benefit of exploiting the contribution to density-dependent regulation of those larvae that are destined to die at a later stage. The juvenile hormone analogue pyriproxyfen (PPF) kills pupae in larval breeding sites by preventing their successful development into adult mosquitoes [52, 53]. In breeding sites receiving a sufficient quantity of PPF there is a 95% reduction in the number of pupae successfully emerging as adult mosquitoes (Greg Devine, personal communication). The application of PPF with efficacy F PPF to a proportion c of breeding sites can be described using the following model,
Parameters for Intervention Impact
To illustrate the impact of combinations of interventions the coverage parameter, c, is varied between 0 (zero coverage) and 1 (100% coverage). For each intervention, this coverage parameter has a slightly different interpretation. For LLINs coverage is the proportion of people sleeping under nets; for IRS coverage is the proportion of houses sprayed; for larvicides or pupacides coverage is the proportion of breeding sites that are treated. We focus on the impact of interventions after deployment before waning of insecticide efficacy becomes significant.
The model in which the carrying capacity is proportional to exponentially weighted past rainfall gave the best fit to the data (see Supplementary Information). All parameter estimates are for this model. The posterior median and 95% credible intervals (95% CrI) for the parameters are shown in Table 1. Comparisons of the predicted and observed mosquito densities for each of the 8 villages in the Garki District are shown in Figure 2.
Mosquito growth rate and reproduction number
In the absence of density-dependent regulation of larvae in breeding sites our parameters result in an estimate of the basic reproduction number for A. gambiae mosquitoes of R 0 = 67 (95% CrI, 40 - 87). This value should be interpreted as a theoretical upper bound for the mosquito reproduction number as density dependence will always play some role. Estimates for the effective mosquito reproduction number Reff(t) for the villages in the Garki Project are shown in Figure 3. These estimates are significantly less than the upper bound of R0 = 67 as population growth is being restricted by competition during the aquatic stages. The largest peak in R eff will occur at the very start of the rainy season when existing breeding sites expand rapidly and new sites are created. For all 8 villages data collection begins just after the start of the first rainy season so the first peak in R eff is missed, although we still have R eff >1 during the growth stage of the mosquito population. For village 7 we have data for the beginning of two rainy seasons where we estimate R eff = 21.36 (95% CrI, 12.74 - 33.35) and R eff = 35.59 (95% CrI, 25.62 - 46.99). For village 8 we have data for the beginning of one rainy season where we estimate R eff = 14.38 (95% CrI, 9.29 - 24.80). For each village the value of R eff will depend on rainfall patterns and the structure of local breeding habitats. In our model the exact value of the peak in R eff will depend on the number of mosquitoes during the dry season; individual mosquitoes in small populations during the dry season will have a greater chance of reproductive success when the rains begin.
Comparison of constant emergence and density-dependent models
Figure 4 shows a comparison of the predicted changes in adult mosquito density due to vector control interventions for a model of constant mosquito emergence and a model of the full mosquito lifecycle with density-dependent regulation in larval breeding sites. The assumption of many vector control models that mosquitoes emerge from breeding sites at a constant rate leads to an underestimate of the reduction in adult mosquito density. This effect becomes more substantial for efficacious interventions such as IRS that are predicted to kill large numbers of adult mosquitoes and hence cause large reductions in the number of eggs oviposited in breeding sites. However the most important effect of vector control interventions on malaria transmission is still the reduction in the proportion of mosquitoes surviving the extrinsic incubation period as has been routinely highlighted in fundamental studies [54–56].
Impact of combined vector control interventions
Vector control interventions are often deployed in settings where other interventions are already in place, most frequently LLINs [2, 46, 47]. Figure 5 shows the predicted changes in adult mosquito density when additional interventions are combined with LLINs at various levels of coverage. At low levels of LLIN coverage (Figures 5A and 5B) IRS causes the largest reduction in mosquito density. At higher levels of LLIN coverage (Figures 5C and 5D) larvicide and pupacide become comparable with, and in some cases, superior to IRS at reducing mosquito density. However, although larval source management may be effective at reducing mosquito numbers, it does not have the additional benefit that LLINs and IRS have in killing the adult mosquito during the malaria parasite's incubation period.
Combining interventions that target different stages of the mosquito's lifecycle can control vector populations enough to eliminate malaria transmission in some settings [5, 57]. Although LLINs and IRS are highly effective interventions, they both target the adult mosquito when feeding on a human host or resting after a blood meal. Thus, if both interventions are deployed together, the barrier effect of the bed net may divert a mosquito from feeding and thus prevent the mosquito from resting on an insecticide-sprayed wall after taking a blood meal. In contrast, interventions targeting different stages of the mosquito's lifecycle, such as LLINs and larvicide, will not interfere with each other.
The larvicide BTI kills both larvae and pupae in breeding sites, and the pupacide PPF prevents pupae from emerging as adult mosquitoes. BTI and PPF are predicted to cause comparable reductions in mosquito numbers (Figure 5). The effect of increased mortality to all aquatic stages due to BTI is similar to the combined effect of PPF-induced pupal mortality and the exploitation of the contribution to density-dependent mortality of larvae destined to die as they become pupae. It may therefore fall to other factors such as cost, ease of application or impact on the ecosystem to choose between these two approaches.
The primary ecological factor controlling A. gambiae s.l. populations is considered to be density-dependent regulation of larval survival in breeding sites [13, 26]. In contrast, density dependence is unlikely to be a significant factor for the adult stages of A. gambiae s.l. . Our analysis of data from an experimental study in Tanzania suggests an approximately linear relationship between larval density and larval death rates when first instars are introduced in a staggered manner. Using this relationship in a model of mosquito dynamics, and assuming seasonal patterns in the larval carrying capacity are driven primarily by rainfall, it was possible to obtain a good fit to data on adult mosquito dynamics in 8 different villages from the Garki Project in Nigeria. Thus density-dependent regulation of the larval stages via an impact on the larval death rates provides a plausible mechanism for seasonal variation in the adult mosquito population.
The model presented here was developed for a single homogeneous site. In practice, a region's larval breeding habitats comprise multiple discrete sites as opposed to a single well mixing site [48, 49, 59, 60]. This variation in type of breeding site, combined with the heterogeneity in breeding site preference by ovipositing females, leads to variation in the degree of density-dependent regulation, with intense intraspecific larval competition in densely populated sites and little or no competition in sparsely populated sites. This heterogeneity leads to more severe regulation of the mosquito population than a scenario where larvae are evenly distributed across all viable breeding sites. Future work will utilise the framework presented here to develop more explicit spatial models that account for this important source of heterogeneity. A second potential limitation to our model of density-dependent regulation is that eggs are oviposited in batches. This intrinsic clumping may cause competitive processes to be present even at low densities. These limitations hinder the model's ability to make predictions at mosquito densities low enough for temporary elimination of the vector population.
Although increased larval mortality is likely to be the main consequence of high larval densities, extended larval development time and the production of smaller larvae are also important factors [13, 16]. Extended larval development times may delay the peak of adult mosquitoes accompanying the start of the rainy season. Smaller larvae at high densities lead to smaller emerging adult mosquitoes, which can affect the epidemiology of malaria in a number of ways: smaller mosquitoes may need to take two blood meals before their first oviposition ; oviposit less eggs ; exhibit increased mortality ; and may ingest a lower number of malaria parasites as they take a smaller blood meal which is digested more quickly (but see Sinden et al.  for density-dependent processes affecting Plasmodium sporogony itself).
Mosquito populations are sensitively dependent to climatic factors, in particular water and air temperatures and rainfall. A. gambiae larval survival and developmental times are dependent on water temperature [8, 9, 14] with faster development and decreased survival generally being observed at higher temperatures. There is only limited data on the interaction between temperature and larval density . This relationship is likely to be complex and non-linear, but a model of the dynamics of larval development at a fixed temperature as demonstrated here may aid in the understanding of the dynamics at variable temperatures. Seasonal patterns of mosquito abundance closely follow seasonal patterns of rainfall. Here we assumed that the environmental carrying capacity is directly proportional to rainfall, however this may not be the case in general as it has been observed that too much rain can flush mosquito larvae from breeding sites . The exact relationship between rainfall and carrying capacity is likely to be complex and dependent on local hydrological conditions.
By fitting our model to the adult mosquito data in the Garki study, we were able to obtain estimates of the variation in the adult mosquito reproduction number over time. Whilst the theoretical maximum in the absence of density-dependent regulation appears high (R0 = 67), in reality density dependence will mean that this is rarely, if ever, achieved. For the seasonal rainfall patterns in Figure 3 the effective reproduction number occurring at the beginning of the rainy season will depend on rainfall patterns and can peak at up to R eff = 35.59. However, this declines rapidly as density-dependent regulation acts and thus at any point during the malaria transmission season the value is likely to be well below this. These estimates have implications for the design of strategies that aim to reduce the adult mosquito population by reducing population growth, including the range of techniques for genetic modification being developed for Anopheles and other mosquito species [65, 66]. For example the results shown in Figure 3 suggest that the optimal time to introduce genetically modified mosquitoes with the aim of replacement of the wild-type species is at the beginning of the rainy season when introduced mosquitoes can exploit low densities in breeding sites to reproduce efficiently.
A substantial intervention programme based on widespread LLINs and IRS can, in theory, increase adult mosquito mortality so that very few mosquitoes survive the duration of sporogony to become infectious. Eliminating the A. gambiae s.l. vector using interventions directed against the adult stages would require increasing adult mosquito mortality to such an extent that nearly all mosquitoes die before completing their first gonotrophic cycle. Such extreme and early mosquito mortalities are difficult to effect with current interventions and risk exerting strong selection pressures on the reproductive fitness of the mosquito population and the rapid development of insecticide resistance. Furthermore, an anti-malaria campaign based on interventions such as LLINs and IRS that targets adult mosquitoes will cause a large decrease in the mosquito population with an accompanying decrease in densities in larval breeding sites. The few mosquitoes that remain are likely to be bigger as they will have undergone less density-dependent larval competition . These remaining mosquitoes may prove to be efficient transmitters of malaria. Thus the impact on adult mosquito populations, and also on malaria transmission, may be less substantial than predicted here.
The most widespread vector control interventions, LLINs and IRS, both target the adult mosquito when it is seeking to feed/has fed on its human host, thus LLINs and IRS are likely to interact inefficiently when deployed together. Our results show that, to improve upon an anti-malarial programme based on high levels of LLIN coverage, to reduce mosquito density interventions should target the mosquito at non-feeding stages. Thus, a modest level of larviciding or pupaciding may be preferable to extensive IRS campaigns (Figure 5D) if high levels (~80%) of coverage of LLINs have been obtained. A potential drawback of larviciding campaigns is the difficulty in maintaining high levels of coverage of larval breeding sites and the need for re-application of insecticide . A possible exception to this could be auto-dissemination of the pupacide PPF from adult females to larval breeding sites  where PPF coverage on a modest proportion of walls could lead to a high level of coverage of breeding sites.
This model simulates the lifecycle of a single mixing A. gambiae s.l. population. In many locations other species of mosquitoes capable of transmitting malaria such as A. funestus, and other culicines may also be present, and interspecific interactions will be important. A. gambiae s.s. and A. arabiensis larvae compete when inhabiting shared breeding sites [9, 68], but exhibit a preference for different sites, with A. gambiae s.s. preferring small temporary pools such as puddles and hoof prints, and A. arabiensis preferring larger, more permanent pools . If both species share the same breeding sites then an intervention targeting the indoor-biting A. gambiae s.s. may facilitate species replacement with the outdoor-biting and more zoophagic A. arabiensis. However if both species breed in different sites, then a decrease in the A. gambiae s.s. population will not be accompanied by an increase in the A. arabiensis population. Instead, there will be an increase in the proportion of A. arabiensis caught compared to A. gambiae s.s. as observed by Bayoh et al.  in Nyanza Province in Kenya when increased ITN coverage caused a crash in the A. gambiae s.s. population but only minor reductions in the A. arabiensis population.
Many of the key anti-malarial interventions target the Anopheles vector, including LLINs, IRS, larviciding and environmental management of breeding sites. In addition to these existing technologies there are a number of promising vector control interventions in various phases of development, including novel insecticides, late-acting biopesticides, attractive toxic sugar bait and the introduction of genetically modified mosquitoes [53, 65, 69, 70]. Mathematical models of mosquito population dynamics can provide insights into the likely impact of these interventions on the vector population, though consideration of possible interactions between the insecticides and the vector competence for Plasmodium would have to be included if projections about the likely impact on malaria transmission are also to be made. Indeed as multiple interventions are deployed targeting different stages of the mosquito's lifecycle, model predictions may prove useful in identifying combinations of interventions that interact synergistically.
The dynamics of Anopheles gambiae populations are driven by climatic factors, rainfall and temperature, and by density-dependent competition in larval breeding habitats. In addition to providing protection to individuals against the bites of malaria infected mosquitoes, vector control interventions can also have a substantial effect on mosquito population dynamics; large reductions in mosquito numbers are frequently seen following the introduction of insecticide-treated nets or indoor residual spraying. Even greater reductions in mosquito numbers are possible by selecting combinations of interventions that target different stages in the mosquito's lifecycle, for example targeting the aquatic stages with larvicide in areas with high ITN coverage.
Molineaux L: The Garki Project. 1980, Geneva: World Health Organisation
Sharp BL, Ridl FC, Govender D, Kuklinski J, Kleinschmidt I: Malaria vector control by indoor residual insecticide spraying on the tropical island of Bioko, Equatorial Guinea. Malaria J. 2007, 6: 52-10.1186/1475-2875-6-52.
Bayoh MN, Mathias DK, Odiere MR, Mutuku FM, Kamau L, Gimnig JE, Vulule JM, Hawley WA, Hamel MJ, Walker ED: Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province, Kenya. Malaria J. 2010, 9: 62-10.1186/1475-2875-9-62.
Draper CC, Smith A: Malaria in the Pare area of Tanganyika II. Effects of three years' spraying of huts with dieldrin. Trans Roy Soc Trop Med and Hyg. 1960, 54 (4): 342-357. 10.1016/0035-9203(60)90115-2.
Ferguson HM, Dornhaus A, Beeche A, Borgemeister C, Gottlieb M, Mulla MS, Gimnig JE, Fish D, Killeen GF: Ecology: A prerequisite for malaria elimination and eradication. PLoS Med. 7: e1000303-
Hogg JC, Hurd H: The effects of natural Plasmodium falciparum infection on the fecundity and mortality of Anopheles gambiae sl in north east Tanzania. Parasitology. 1997, 114: 325-331. 10.1017/S0031182096008542.
Hogg JC, Thomson MC, Hurd H: Comparative fecundity and associated factors for two sibling species of the Anopheles gambiae complex occurring sympatrically in The Gambia. Med Vet Ent. 1996, 10 (4): 385-391. 10.1111/j.1365-2915.1996.tb00761.x.
Bayoh MN, Lindsay SW: Effect of temperature on the development of the aquatic stages of Anopheles gambiae sensu stricto (Diptera: Culicidae). Bull Entomol Res. 2003, 93 (5): 375-381.
Kirby MJ, Lindsay SW: Effect of temperature and inter-specific competition on the development and survival of Anopheles gambiae sensu stricto and An. arabiensis larvae. Acta Tropica. 2009, 109 (2): 118-123. 10.1016/j.actatropica.2008.09.025.
Gething PW, Van Boeckel TP, Smith DL, Guerra CA, Patil AP, Snow RW, Hay SI: Modelling the global constraints of temperature on transmission of Plasmodium falciparum and P. vivax. Parasit Vectors. 2011, 4: 92-10.1186/1756-3305-4-92.
Bellows TS: The descriptive properties of some models for density dependence. J Animal Ecol. 1981, 50 (1): 139-156. 10.2307/4037.
Njunwa K: Studies on the productivity of Anopheles breeding sites in relation to adult mosquito density. PhD Thesis, London School of Hygiene and Tropical Medicine. 1993
Gimnig JE, Ombok M, Otieno S, Kaufman MG, Vulule JM, Walker ED: Density-dependent development of Anopheles gambiae (Diptera: Culicidae) larvae in artificial habitats. J Med Entomol. 2002, 39 (1): 162-172. 10.1603/0022-2585-39.1.162.
Lyimo EO, Takken W, Koella JC: Effect of rearing temperature and larval density on larval survival, age at pupation and adult size of Anopheles-gambiae. Entomologia Experimentalis Et Applicata. 1992, 63 (3): 265-271. 10.1111/j.1570-7458.1992.tb01583.x.
Koenraadt CJM, Takken W: Cannibalism and predation among larvae of the Anopheles gambiae complex. Med Vet Entomol. 2003, 17 (1): 61-66. 10.1046/j.1365-2915.2003.00409.x.
Russell TL, Lwetoijera DW, Knols BGJ, Takken W, Killeen GF, Ferguson HM: Linking individual phenotype to density-dependent population growth: the influence of body size on the population dynamics of malaria vectors. Proc Roy Soc B. 2011
Smith DL, Hay SI, Noor AM, Snow RW: Predicting changing malaria risk after expanded insecticide-treated net coverage in Africa. Trends Parasitol. 2009, 25 (11): 511-516. 10.1016/j.pt.2009.08.002.
Chitnis N, Schapira A, Smith T, Steketee R: Comparing the effectiveness of malaria vector-control interventions through a mathematical model. Am J Trop Med Hyg. 2010, 83 (2): 230-240. 10.4269/ajtmh.2010.09-0179.
Yakob L, Yan GY: Modeling the effects of integrating larval habitat source reduction and insecticide treated nets for malaria control. PLoS One. 2009, 4: e6921-10.1371/journal.pone.0006921.
Gu WD, Regens JL, Beier JC, Novak RJ: Source reduction of mosquito larval habitats has unexpected consequences on malaria transmission. P Natl Acad Sci USA. 2006, 103 (46): 17560-17563. 10.1073/pnas.0608452103.
Gu WD, Novak RJ: Predicting the impact of insecticide-treated bed nets on malaria transmission: the devil is in the detail. Malaria J. 2009, 8: 256-10.1186/1475-2875-8-256.
Gu WD, Novak RJ: Habitat-based modeling of impacts of mosquito larval interventions on entomological inoculation rates, incidence, and prevalence of malaria. Am J Trop Med Hyg. 2005, 73 (3): 546-552.
Killeen GF, Knols BGJ, Gu WD: Taking malaria transmission out of the bottle: implications of mosquito dispersal for vector-control interventions. Lancet Inf Dis. 2003, 3 (5): 297-303. 10.1016/S1473-3099(03)00611-X.
Smith DL, Dushoff J, McKenzie FE: The risk of a mosquito-borne infection in a heterogeneous environment. PLoS Biol. 2004, 2: 368-10.1371/journal.pbio.0020368.
Worrall E, Connor SJ, Thomson MC: A model to simulate the impact of timing, coverage and transmission intensity on the effectiveness of indoor residual spraying (IRS) for malaria control. Trop Med Int Health. 2007, 12 (1): 75-88.
Depinay JMO, Mbogo CM, Killeen G, Knols B, Beier J, Carlson J, Dushoff J, Billingsley P, Mwambi H, Githure J: A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malaria J. 2004, 3: 29-10.1186/1475-2875-3-29.
Hancock PA, Godfray HCJ: Application of the lumped age-class technique to studying the dynamics of malaria-mosquito-human interactions. Malaria J. 2007, 6: 98-10.1186/1475-2875-6-98.
Ermert V, Fink AH, Jones AE, Morse AP: Development of a new version of the Liverpool Malaria Model. I. Refining the parameter settings and mathematical formulation of basic processes based on a literature review. Malaria J. 2011, 10: 35-10.1186/1475-2875-10-35.
Ermert V, Fink AH, Jones AE, Morse AP: Development of a new version of the Liverpool Malaria Model. II. Calibration and validation for West Africa. Malaria J. 2011, 10: 62-10.1186/1475-2875-10-62.
Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, Bousema T, Drakeley CJ, Ferguson NM, Basanez MG: Reducing Plasmodium falciparum malaria transmission in Africa: A model-based evaluation of intervention strategies. PLoS Med. 2010, 7: e1000324-10.1371/journal.pmed.1000324.
Beier JC: Frequent blood-feeding and restrictive sugar-feeding behavior enhance the malaria vector potential of Anopheles gambiae sl and An-funestus (Diptera: Culicidae) in western Kenya. J Med Entomol. 1996, 33 (4): 613-618.
Clements AN, Paterson GD: The analysis of mortality and survival rates in wild populations of mosquitos. J App Ecol. 1981, 18 (2): 373-399. 10.2307/2402401.
Dawes EJ, Churcher TS, Zhuang S, Sinden RE, Basanez MG: Anopheles mortality is both age- and Plasmodium-density dependent: implications for malaria transmission. Malaria J. 2009, 8 (228):
Wearing HJ, Rohani P, Keeling MJ: Appropriate models for the management of infectious diseases. PLoS Med. 2005, 2: e174-10.1371/journal.pmed.0020174.
Killeen GF, McKenzie FE, Foy BD, Schieffelin C, Billingsley PF, Beier JC: A simplified model for predicting malaria entomologic inoculation rates based on entomologic and parasitologic parameters relevant to control. Am J Trop Med Hyg. 2000, 62 (5): 535-544.
Service MW: Studies on sampling larval populations of Anopheles-gambiae complex. Bulletin of the World Health Organization. 1971, 45 (2): 169-
Munga S, Minakawa N, Zhou G, Githeko AK, Yan G: Survivorship of immature stages of Anopheles gambiae s.l. (Diptera: culicidae) in natural habitats in western Kenya highlands. J Med Entomol. 2007, 44 (5): 758-764. 10.1603/0022-2585(2007)44[758:SOISOA]2.0.CO;2.
Ye Y, Hoshen M, Kyobutungi C, Louis VR, Sauerborn R: Local scale prediction of Plasmodium falciparum malaria transmission in an endemic region using temperature and rainfall. Global Health Action. 2009, 2: 1923-
Mala AO, Irungu LW, Shililu JI, Muturi EJ, Mbogo CC, Njagi JK, Githure JI: Dry season ecology of Anopheles gambiae complex mosquitoes at larval habitats in two traditionally semi-arid villages in Baringo, Kenya. Parasit Vectors. 2011, 4: 25-10.1186/1756-3305-4-25.
Craig MH, Snow RW, le Sueur D: A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitology Today. 1999, 15 (3): 105-111. 10.1016/S0169-4758(99)01396-4.
Paaijmans KP, Imbahale SS, Thomas MB, Takken W: Relevant microclimate for determining the development rate of malaria mosquitoes and possible implications of climate change. Malaria J. 2010, 9: 196-10.1186/1475-2875-9-196.
Service MW: Mortalities of larvae of Anopheles-gambiae giles complex and detection of predators by precipitin test. B Entomol Res. 1973, 62 (3): 359-369. 10.1017/S0007485300003862.
Mwangangi JM, Muturi EJ, Shililu J, Muriu SM, Jacob B, Kabiru EW, Mbogo CM, Githure J, Novak R: Survival of immature Anopheles arabiensis (Diptera: Culicidae) in aquatic habitats in Mwea rice irrigation scheme, central Kenya. Malaria J. 2006, 5: 114-10.1186/1475-2875-5-114.
Le Menach A, Takala S, McKenzie FE, Perisse A, Harris A, Flahault A, Smith DL: An elaborated feeding cycle model for reductions in vectorial capacity of night-biting mosquitoes by insecticide-treated nets. Malaria J. 2007, 6: 10-10.1186/1475-2875-6-10.
Kroeger A, Horstick O, Riedl C, Kaiser A, Becker N: The potential for malaria control with the biological larvicide bacillus-thuringiensis israelensis (bti) in Peru and Ecuador. Acta Tropica. 1995, 60 (1): 47-57. 10.1016/0001-706X(95)00101-J.
Fillinger U, Ndenga B, Githeko A, Lindsay SW: Integrated malaria vector control with microbial larvicides and insecticide-treated nets in western Kenya: a controlled trial. Bulletin of the World Health Organization. 2009, 87 (9): 655-665. 10.2471/BLT.08.055632.
Majambere S, Pinder M, Fillinger U, Ameh D, Conway DJ, Green C, Jeffries D, Jawara M, Milligan PJ, Hutchinson R: Is mosquito larval source management appropriate for reducing malaria in areas of extensive flooding in The Gambia? A cross-over intervention trial. Am J Trop Med Hyg. 2010, 82 (2): 176-184. 10.4269/ajtmh.2010.09-0373.
Sattler MA, Mtasiwa D, Kiama M, Premji Z, Tanner M, Killeen GF, Lengeler C: Habitat characterization and spatial distribution of Anopheles sp. mosquito larvae in Dar es Salaam (Tanzania) during an extended dry period. Malaria J. 2005, 4: 4-10.1186/1475-2875-4-4.
Fillinger U, Sonye G, Killeen GF, Knols BGJ, Becker N: The practical importance of permanent and semipermanent habitats for controlling aquatic stages of Anopheles gambiae sensu lato mosquitoes: operational observations from a rural town in western Kenya. Trop Med Int Health. 2004, 9 (12): 1274-1289. 10.1111/j.1365-3156.2004.01335.x.
Kweka EJ, Zhou G, Lee M-C, Gilbreath TM, Mosha F, Munga S, Githeko AK, Yan G: Evaluation of two methods of estimating larval habitat productivity in western Kenya highlands. Parasit Vectors. 2011, 4: 110-10.1186/1756-3305-4-110.
Killeen GF, Tanner M, Mukabana WR, Kalongolela MS, Kannady K, Lindsay SW, Fillinger U, de Castro MC: Habitat targeting for controlling aquatic stages of malaria vectors in Africa. Am J Trop Med Hyg. 2006, 74 (4): 517-518.
Devine GJ, Killeen GF: The potential of a new larviciding method for the control of malaria vectors. Malaria J. 9: 142-
Devine GJ, Perea EZ, Killeen GF, Stancil JD, Clark SJ, Morrison AC: Using adult mosquitoes to transfer insecticides to Aedes aegypti larval habitats. Proc Natl Acad Sci USA. 2009, 106 (28): 11530-11534. 10.1073/pnas.0901369106.
Macdonald G: The epidemiology and control of malaria. 1957, London: Oxford University Press
Dietz K, Molineaux L: A malaria model tested in the African savannah. Ninth International Congress on Tropical Medicine and Malaria Athens 14-21 October 1973 Volume 1 Abstracts of invited papers. 1973, 297-
Molineaux L, Dietz K, Thomas A: Further epidemiological evaluation of a malaria model. Bulletin of the World Health Organization. 1978, 56 (4): 565-571.
WHO: Global strategic framework for integrated vector management. 2004, Geneva: World Health Organisation
Charlwood JD, Smith T, Kihonda J, Heiz B, Billingsley PF, Takken W: Density-independent feeding success of malaria vectors (Diptera, Culicidae) in Tanzania. B Entomol Res. 1995, 85 (1): 29-35. 10.1017/S0007485300051981.
Majambere S, Fillinger U, Sayer DR, Green C, Lindsay SW: Spatial distribution of mosquito larvae and the potential for targeted larval control in The Gambia. Am J Trop Med Hyg. 2008, 79 (1): 19-27.
Fillinger U, Sombroek H, Majambere S, van Loon E, Takken W, Lindsay SW: Identifying the most productive breeding sites for malaria mosquitoes in The Gambia. Malaria J. 2009, 8: 62-10.1186/1475-2875-8-62.
Lyimo EO, Takken W: Effects of adult body-size on fecundity and the pre-gravid rate of Anopheles-gambiae females in Tanzania. Med Vet Entomol. 1993, 7 (4): 328-332. 10.1111/j.1365-2915.1993.tb00700.x.
Takken W, Klowden MJ, Chambers GM: Effect of body size on host seeking and blood meal utilization in Anopheles gambiae sensu stricto (Diptera: Culicidae): The disadvantage of being small. J Med Entomol. 1998, 35 (5): 639-645.
Sinden RE, Dawes EJ, Alavi Y, Waldock J, Finney O, Mendoza J, Butcher GA, Andrews L, Hill AV, Gilbert SC: Progression of Plasmodium berghei through Anopheles stephensi is density-dependent. PLoS Path. 2007, 3: e195-10.1371/journal.ppat.0030195.
Paaijmans KP, Wandago MO, Githeko AK, Takken W: Unexpected high losses of Anopheles gambiae larvae due to rainfall. PLoS One. 2007, 2: e1146-10.1371/journal.pone.0001146.
Alphey L, Beard CB, Billingsley P, Coetzee M, Crisanti A, Curtis C, Eggleston P, Godfray C, Hemingway J, Jacobs-Lorena M: Malaria control with genetically manipulated insect vectors. Science. 2002, 298 (5591): 119-121. 10.1126/science.1078278.
Phuc HK, Andreasen MH, Burton RS, Vass C, Epton MJ, Pape G, Fu GL, Condon KC, Scaife S, Donnelly CA: Late-acting dominant lethal genetic systems and mosquito control. Bmc Biol. 2007, 5: 11-10.1186/1741-7007-5-11.
Chaki PP, Govella NJ, Shoo B, Hemed A, Tanner M, Fillinger U, Killeen GF: Achieving high coverage of larval-stage mosquito surveillance: challenges for a community-based mosquito control programme in urban Dar es Salaam, Tanzania. Malaria J. 2009, 8: 311-10.1186/1475-2875-8-311.
Schneider P, Takken W, McCall PJ: Interspecific competition between sibling species larvae of Anopheles arabiensis and An. gambiae. Med Vet Entomol. 2000, 14 (2): 165-170. 10.1046/j.1365-2915.2000.00204.x.
Hancock PA, Thomas MB, Godfray HCJ: An age-structured model to evaluate the potential of novel malaria-control interventions: a case study of fungal biopesticide sprays. Proc Roy Soc B. 2009, 276 (1654): 71-80. 10.1098/rspb.2008.0689.
Muller GC, Kravchenko VD, Schlein Y: Decline of Anopheles sergentii and Aedes caspius populations following presentation of attractive toxic (spinosad) sugar bait stations in an oasis. J Am Mosq Cont Ass. 2008, 24 (1): 147-149. 10.2987/8756-971X(2008)24[147:DOASAA]2.0.CO;2.
The authors would like to thank Jo Lines, Steve Lindsay, Bob Sinden and Greg Devine for helpful discussions and three anonymous referees for valuable comments. MTW, NMF and ACG acknowledge support from the MRC Centre grant. JTG, NMF and ACG acknowledge support from the Bill and Melinda Gates Foundation Vaccine Modeling Initiative. MGB and TSC acknowledge support from the European Commission FP7 collaborative project TransMalariaBloc (HEALTH-F3-2008-223736).
The authors declare that they have no competing interests.
MTW, ACG MGB and ACG conceived the analysis. JTG and MTW performed the model fitting and generated the model output. MTW wrote the first draft of the manuscript. All authors contributed to the final version of the manuscript.