Figures
Abstract
An important question for mosquito population dynamics, mosquito-borne pathogen transmission and vector control is how mosquito populations are regulated. Here we develop simple models with heterogeneity in egg laying patterns and in the responses of larval populations to crowding in aquatic habitats. We use the models to evaluate how such heterogeneity affects mosquito population regulation and the effects of larval source management (LSM). We revisit the notion of a carrying capacity and show how heterogeneity changes our understanding of density dependence and the outcome of LSM. Crowding in and productivity of aquatic habitats is highly uneven unless egg-laying distributions are fine-tuned to match the distribution of habitats’ carrying capacities. LSM reduces mosquito population density linearly with coverage if adult mosquitoes avoid laying eggs in treated habitats, but quadratically if eggs are laid in treated habitats and the effort is therefore wasted (i.e., treating 50% of habitat reduces mosquito density by approximately 75%). Unsurprisingly, targeting (i.e. treating a subset of the most productive pools) gives much larger reductions for similar coverage, but with poor targeting, increasing coverage could increase adult mosquito population densities if eggs are laid in higher capacity habitats. Our analysis suggests that, in some contexts, LSM models that accounts for heterogeneity in production of adult mosquitoes provide theoretical support for pursuing mosquito-borne disease prevention through strategic and repeated application of modern larvicides.
Citation: Smith DL, Perkins TA, Tusting LS, Scott TW, Lindsay SW (2013) Mosquito Population Regulation and Larval Source Management in Heterogeneous Environments. PLoS ONE 8(8): e71247. https://doi.org/10.1371/journal.pone.0071247
Editor: Nikos Vasilakis, University of Texas Medical Branch, United States of America
Received: March 2, 2013; Accepted: June 26, 2013; Published: August 7, 2013
Copyright: © 2013 Smith et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: DLS acknowledges funding from the Bloomberg Family Foundation. DLS and SWL acknowledge funding from National Institutes of Health/National Institute of Allergy and Infectious Diseases (U19AI089674). DLS, TAP, SWL, and TWS acknowledge funding from the Research and Policy for Infectious Disease Dynamics (RAPIDD) program of the Science and Technology Directorate, Department of Homeland Security, and the Fogarty International Center, National Institutes of Health.
Competing interests: Co-author David Smith has served as a PLOS ONE Editorial Board member, but this does not alter the authors’ adherence to all the PLOS ONE policies on sharing data and materials.
Introduction
Dynamic models of malaria transmission have influenced strategic decisions about disease prevention from the time of Ronald Ross in the early 20th century, when larval source management (LSM) was the dominant form of vector control [1]. After early field deployment of DDT demonstrated that indoor residual spraying (IRS) was an extremely effective way to control malaria, George Macdonald’s influential mathematical analysis showed that transmission was highly sensitive to adult mosquito mortality rates [2]. This analysis and emerging theory reinforced the prevailing notion at the time that DDT was a sufficient tool for malaria eradication [3], [4], and IRS was implemented largely to the exclusion of LSM. The legacy of Macdonald’s sixty-year old analysis can be seen in contemporary policy decisions by leading international organizations, including a recent evaluation by the World Health Organization (WHO) that was highly critical of larviciding in sub-Saharan Africa [5]. These recommendations, based largely on the Ross-Macdonald model that lacks dynamic mosquito populations and is ill suited to evaluate LSM, come despite evidence that LSM can achieve similar results and at a similar cost to ITNs [6]. Here, we re-examine the simple models that have motivated such analyses, and we derive some basic lessons for mosquito population dynamic and control to guide policy for LSM.
Several mosquito population dynamic models have been developed that link adult and immature aquatic populations [7]–[11], and a few have explicitly considered LSM [12]. Most models of larval populations, whether simple or complex, make some assumption about density dependence and population regulation. Some have considered the complex structure that arises from having populations of eggs, four larval instars, and pupae [13]. Others have considered the dynamics of systems with predators or resource-based competition [14]. Complicated computer-simulation models have considered the effects of heterogeneity in rainfall and temperature, heterogeneous habitat geometries with variable responses to flushing, and desiccation [9], [10], [14]–[17]. Finally, a few models have considered how the distribution of larval habitat constrains egg laying and affects the adult mosquito population dynamics and pathogen transmission [18]–[20]. It remains unclear how heterogeneity affects the way mosquito populations are regulated and what variation in key processes means for LSM. Here, we present a simple theoretical framework that can be used to understand habitat heterogeneity, the local effects of density dependence, the factors that affect the outcome of LSM in dynamic, heterogeneous environments, and their total effects on pathogen transmission.
Methods
Many factors have been implicated in immature mosquito population dynamics, including egg laying, water temperature, resource limitations, predation, larval development rates, the ephemeral availability of mosquito habitat due to evaporation and desiccation, or filling or flushing habitat from the combination of rainfall and hydrology [21]. Here, we take a simpler approach that focuses narrowly population regulation when aquatic habitats are heterogeneous. The model may not be suitable for some purposes, such as simulating mosquito population dynamics when realistic lags for mosquito development are required (i.e. see [13]), but the models do provide insights into the regulation of mosquito populations, population dynamics in heterogeneous habitats and the effects of LSM. These lessons can, perhaps, serve as a theoretical basis for understanding more complicated models.
The Mosquito Population Dynamic Model
The following model considers the coupled dynamics of aquatic immature and terrestrial adult mosquito populations. We assume the population of larval mosquitoes is subdivided into distinct aquatic habitats. Individual aquatic habitats are hereafter called “pools” to facilitate communication, even though this may not be the best description of many kinds of larval habitats.
Let denote the population density of adult mosquitoes at time
and let
denote the per-capita death rate. The number of larvae in each pool at any given time is denoted
Let
denote the mosquito blood feeding rate,
the number of eggs laid by a mosquito each egg laying cycle, and
the fraction of eggs laid in the
pool. In the
pool, mosquitoes are assumed to mature at rate
and die at the per-capita rate
where
represents a pool-specific increase in per-capita mortality in response to crowding. For
which was assumed for most of our analysis, the relationship gives mean crowding, which is analogous to the classical first-order description of density dependence as described by the logistic growth equation. Under these assumptions, mosquito population dynamics are described by the following equations:
(1)
(2)
Homogeneous environments were defined by letting each pool have identical parameters and by distributing eggs evenly among the pools. Environments were made heterogeneous by varying parameters describing larval dynamics or egg laying (i.e. or
from Eq. 1) strategically to illustrate specific aspects of this system. We constructed completely heterogeneous environments by drawing random numbers for all larval dynamic parameters and for egg laying. Parameter names are summarized in Table S1 along with all the values used for the simulations. The population dynamics in these completely heterogeneous environments depends on some notion of the response to crowding, the distribution of eggs laid, and the potential capacity for adult mosquito production of each pool.
Larval Source Management
LSM was simulated by assuming that control was applied either permanently or repeatedly to a subset of these pools, which were called “treated.” This is done to make the analysis simpler and to illustrate properties of the models, even though there might be real constraints on the ability to completely and permanently nullify mosquito productivity. Coverage was defined as the proportion of aquatic habitats that were treated. In our simulations, LSM was assumed to prevent all larval development and eliminate productivity such that no adults emerged from treated pools. The analysis focused on the relationship between coverage and the “control effect size” on transmission, defined for LSM as the proportional decline in the adult mosquito densities when compared to the same system without control.
The control effect sizes of LSM were simulated under two different assumptions about changes in egg-laying behaviour of adult mosquitoes in response to LSM. First, mosquitoes could continue to lay eggs in the pools that had been treated, such as when modern non-repellent larvicides are applied to aquatic habitats. Second, mosquitoes could avoid treated pools and lay eggs elsewhere, either because larvicides in the water acted as repellents or because the habitat was modified or destroyed.
Control effect sizes of increased LSM coverage were examined for five classes of population dynamic simulations based on different assumptions about crowding and egg laying: (1) a homogeneous environment where all pools have identical parameters and eggs are laid evenly; (2) a simple extension of the homogeneous model in which a fraction of habitats in a homogeneous environment were simply non-productive, so that the fraction of eggs laid in productive habitats summed to less than 1; (3) LSM was applied in random order in a completely heterogeneous environment; (4) LSM was “targeted” by treating subsets of the most productive pools in a completely heterogeneous environment (this was done in a perfectly efficient order, such that as coverage increased, the pools with highest productivity were treated first); and (5) to show a contrast, LSM was then inefficiently targeted in a completely heterogeneous environment by treating subsets of the least productive pools.
Results
Mosquito Population Dynamics
The key dynamic feature of the equations describing mosquito population dynamics is that emerging adult mosquitoes become part of an adult mosquito population and that they distribute eggs among many independent pools (Eq. 2). Because of many factors affecting the distribution of eggs in habitats of differing qualities, including the patterns of blood feeding, different patterns emerge from examining the dynamics of completely heterogeneous systems compared with homogeneous systems.
The “carrying capacity” was defined as the equilibrium density of larvae in a system with only one pool or in a homogeneous system, and it is given by the formula:
In the homogeneous system, carrying capacity determines the maximum productivity of each pool, the emergence rate of adult mosquitoes
When egg-laying patterns are heterogeneous, larval density becomes decoupled from carrying capacity. The number of eggs laid and the mean crowding of each habitat affect larval densities. At the steady state, the number of eggs laid in a pool is and:
The number of adults emerging from any individual pool (i.e., productivity), depends on the number of eggs laid and the functional relationship that determines how larval mosquito mortality increases with crowding (Fig. 1a). (Different rules would likely follow from considering predation or other types of regulation that respond dynamically to larval population density.) These productivity curves show that carrying capacity is but one point along a continuum of adult output rates from a pool in relation to the egg input rate. Larval populations thin in response to crowding, so that the proportion of adults emerging from individual pools decreases with the number of eggs laid, but the number of adults emerging increases. Dynamics of individual pools linked by egg-laying females do not, therefore, generally conform to the rules of logistic growth. In heterogeneous environments, productivity and carrying capacity are therefore not generally given by the same quantity.
a) The functional relationship between the rate of egg-laying and productivity depends on the functional response to crowding. In this model, the relationship is sensitive to the power-law scaling relationship ( blue;
red;
purple). Carrying capacity is given for a single value of egg laying rates, given at the steady state if that pool had existed in isolation. b) In a system with 2 pools linked by egg-laying, where the carrying capacity of pool 1 is approximately 90% of the total (dashed blue line) and pool 2 has the rest (dashed red line), the population totals overall (solid black) are generally below the maximum, unless egg laying is fine-tuned such that the proportion of eggs laid was equal to that pool’s proportion of carrying capacity (vertical grey). c) A comparison of productivity (red) and carrying capacity (black line) for a typical set of heterogeneous aquatic habitats. Productivity equals carrying capacity when the distribution of eggs laid is finely tuned to match the distribution of carrying capacities (i.e.
). d) The ratio of productivity to carrying capacity was computed for 100 sets of heterogeneous aquatic habitat. The green line plots the 1∶1 ratio, when productivity equals carrying capacity. These distributions, plotted here as the median (solid line) and the 10th and 90th quantiles (dashed lines), shows the robust pattern that the habitats with the lowest productivity tend to be under capacity and the few highly productive habitats tend to be over capacity.
Some properties of the general system come from exploring a simple system of two pools with different carrying capacities. By varying the proportion of eggs laid in each pool, productivities of the individual pools and of the whole system were compared. In this system, the total productivity equals the total carrying capacity only when the proportion of eggs laid in each pool is equal to that pool’s carrying capacity as a fraction of the total of all pool’s capacities (i.e., if , Fig. 1b).
Numerical simulations demonstrate that this rule holds in completely heterogeneous systems (Fig. 1c). In that case, total productivity is equal to the total carrying capacity only when the distribution of egg laying is fine-tuned to equal the relative distribution of carrying capacity. Unless the proportion of eggs laid is fine-tuned to match the carrying capacities, larval densities will be different than carrying capacity (Fig. 1c,d), often by a large margin. The net effect of this mismatch is unpredictable, but it will depend strongly upon the proportion of eggs laid in the most productive habitats. For the mathematical assumptions made in this model, productivity was lower than capacity in approximately one-third of the cases, but productivity often exceeded carrying capacity. In at least one case, productivity exceeded capacity by 270%.
Intuitively, the dynamic interplay of mobile adults, distributed aquatic habitats, and the response to crowding means that total productivity is strongly affected by the correlation between the distribution of eggs laid and the distribution of carrying capacities in aquatic habitats. The proportion of eggs that survive to become adults in any one pool is reduced as egg laying increases crowding, but the number of eggs being laid depends on the whole ensemble of aquatic habitats. Pools that receive the most eggs will tend to have population densities that exceed their carrying capacities, while those that have the fewest eggs are usually below capacity (Fig. 1d). Crowding will be uneven and the effects of crowding in just a few pools dominate population regulation.
Productivity and carrying capacity should both be correlated with egg-laying (Fig. 2a–c), but the underlying functional relationship between eggs in and adults out is only revealed by plotting the ratio of egg-laying to carrying capacity against the ratio of productivity to carrying capacity (Fig. 2d). Though enlightening, this relationship may not have any practical use unless it is possible to measure carrying capacity directly, perhaps through surrogate measures such as a pool’s surface area, volume, or key resource levels through bioassays [22]. The notion of a carrying capacity is, therefore, useful both conceptually and theoretically. Capacity is not, however, what is typically observed in individual pools or in populations at the steady state. Instead, productivity is determined by the carrying capacities of the individual pools, heterogeneity in egg laying proportions, and the mismatch between the two patterns.
In completely heterogeneous environments, there may be a poor correlation between a) carrying capacity and productivity; b) egg laying and productivity; and c) egg laying and carrying capacity. d) The crowding law governing density dependence is found by plotting the ratio of eggs laid to carrying capacity against the ratio of productivity to carrying capacity (i.e. ). The constant
was used to scale the x-axis.
Another important principle is that, in the absence of immigration from pools outside the study area, the dynamic feedback between egg-laying and aquatic population dynamics is subject to a threshold phenomenon governing mosquito persistence. A sufficient condition for mosquito persistence is The mosquito population can, in theory, persist if at least one adult male and female mosquito is expected to emerge from a pool from an egg laid by a typical single adult mosquito originating from that pool (Analysis S1).
Larval Source Management
The control effect sizes of LSM depend strongly upon the adult female mosquito’s egg-laying behaviour in response to LSM. The most important difference is whether mosquitoes continue to lay eggs in treated habitat. If mosquitoes avoid laying eggs in the treated habitats, then LSM simply reduces the amount of habitat available. The outcomes tend to be consistent with a common use of Macdonald’s formula for R0 with respect to LSM, which assumes linear responses. The dynamics of LSM with heterogeneous biting and targeting are more complicated, however, and different rules govern systems in which adults continue to lay eggs in treated pools. Some general aspects of LSM are best illustrated in homogeneous systems, but other aspects play out differently in heterogeneous systems.
LSM is more effective when mosquitoes continue to lay eggs in treated habitats. To illustrate why this behaviour changes the control dynamic, consider a simple heterogeneous system in which pools either have all the same carrying capacity, or they produce no adults at all. Holding the number of productive pools fixed, the total productivity of the system declines linearly with the number of unproductive habitats. The existence of unproductive aquatic habitats nearby can thus become “egg sinks” [23], [24] and reduce the proportion of eggs laid in the productive pools (Fig. 3a).
a) Holding the total number of productive pools constant, adult mosquito population density declines as the number of unproductive pools increases and absorb eggs. b) The “egg sink” effect gives a non-linear effect to LSM if adult mosquitoes continue to lay eggs in the treated pools, so that treating 50% of the pools reduces adult density by 75%, and treating 75% of the pools reduces adult density by 95% (red). If adult mosquitoes do not lay eggs in the treated pools, however, then reductions in mosquito density are proportional to the % of habitat treated (blue). c) The change in adult mosquito density due to LSM in highly heterogeneous habitat as a function of the proportion of habitats treated depending on whether the adults lay eggs in treated pools (red) or avoid treated pools (blue), and depending on whether LSM was done in one particular random order (grey spikes), perfectly efficiently targeted (dashed lines), or perfectly inefficiently targeted (dotted lines). The black line represents a linear response with respect to coverage. d) For the same graphs as 3c, the effect sizes are plotted on a semi-log scale to highlight the benefits of LSM at high coverage. The best case for this system, with efficient targeting and egg-sink effects, predicts a hundred-fold (99%) reduction in mosquito density for 60% coverage. These benefits also get larger for higher coverage and show that there is enormous potential for LSM to reduce transmission through targeted repeated application of modern larvicides.
To illustrate the relationship between coverage and control effect sizes, LSM was simulated in a homogeneous system with varying coverage levels and with both types of egg-laying responses. In these simulations, when adult mosquitoes continue to lay eggs in treated habitats, there are two effects of LSM. One effect, the reductions in the amount of productive habitat, is complemented by a second effect, an increase in the amount of habitat that serves as a sink for eggs. The two effects are multiplicative, so control effect sizes scale with LSM coverage in a way that is approximately quadratic: removing 50% of habitat reduces mosquito densities by approximately 75%, and removing 80% of the habitat reduces mosquito densities by approximately 96% (Figure 3b).
Similar results occur when habitat is heterogeneous, but the interpretation of “coverage” must be considered in a more nuanced way. In homogeneous systems, coverage describes reductions in capacity, productivity, and egg laying. In heterogeneous systems, however, the mismatch between productivity, capacity, and egg laying means that varying amounts of these three quantities remain as coverage increases. The control effect sizes of LSM in heterogeneous systems thus depend on both the adult egg-laying behaviour in response to LSM and the order that LSM is applied to the pools (Fig. 3c,d).
Not surprisingly, the control effect sizes of LSM would be substantially larger if LSM were targeted at the most productive pools (Fig. 3c,d), and it would be substantially less efficient if not. The most efficient solution – targeting the most productive pools in rank order of their productivity from most to least – results in sharp increase in control effect sizes for even small coverage levels, regardless of adult mosquito egg-laying behaviour. The analysis here suggests that the control effect sizes are greater than log-linear, such that it is possible to reduce transmission by a hundred-fold with moderate coverage through targeted, repeated application of modern (i.e., non-repellent) larvicides and other modes of LSM that create an egg sink effect among the most productive pools.
Like homogeneous systems, the outcome of LSM in heterogeneous systems is also dependent on egg-laying behaviour of mosquitoes in treated pools. The effects of efficient targeting are similar between both types of egg-laying responses, but control effect sizes are always higher, all else equal, when mosquitoes continue to lay eggs in treated pools. Like the homogeneous systems, the “egg-sink” effect in heterogeneous systems complements the removal effect to further reduce population densities (Fig. 3c,d). The magnitude of the egg sink effect varies, however, because the mismatch between productivity and the fraction of eggs laid means that the egg-sink effect is only approximately linear with respect to LSM coverage.
Control effect sizes are, on the other hand, highly variable as a function of coverage when LSM is applied to pools in a random sequence. For the same random sequence, control effect sizes are always higher when eggs are laid in treated pools (Fig. 3c,d). Even with perfectly inefficient targeting, control effect sizes in relation to coverage are nearly linear when eggs are laid in treated pools.
The outcomes of LSM were surprising, however, for some random sequences and for inefficient targeting in the case when adult mosquitoes avoid treated pools and redistribute eggs elsewhere (Fig. 3c,d). In the case of perfectly inefficient targeting (i.e., when a subset of the least productive pools is treated), eggs are redistributed from less to more productive pools. The effect is counteracted by a reduction in total carrying capacity. The net effects change with coverage and with the particular distribution of pools (Fig. 3c,d). Similarly, for a random sequence of pools, productivity can increase as coverage increases whenever the effect of redistributing eggs to more productive habitats is greater than the loss of capacity. The behavioural responses of mosquitoes thus make it possible for LSM to increase overall mosquito density by forcing adult mosquitoes to redistribute eggs in more productive pools when egg-laying under natural conditions is highly inefficient.
Discussion
We conclude that mosquito egg laying is an important factor for mosquito population dynamics and that it can have strong affects on the outcome of LSM. In our analyses, increasing coverage caused quadratic reductions in mosquito density if mosquitoes continued to lay eggs in treated pools. Therefore LSM has the potential to be a highly effective method of malaria control without extensive coverage. In particular, we predict that moderate coverage targeted at the most productive aquatic habitats can achieve reasonably large reductions in transmission. Some of these conclusions are inconsistent with statements from the recent WHO report, which was based on a Ross-Macdonald model perspective [5]. Moreover, recent evidence demonstrates that, in some circumstances, LSM has been effective in reducing clinical malaria outcomes [25]–[31], and for similar costs to those of IRS and long-lasting insecticidal nets (LLINs) [6], [32]–[34]. Given all these caveats, generalizations about LSM are likely to depend heavily on the local context for pathogen transmission and operational constraints.
Our results highlight the lack of attention paid to heterogeneity in mosquito population dynamics and in considering the outcome of LSM. Habitat heterogeneity and local density dependence change the way that dynamics in mosquito population are regulated and could play a role in creating or controlling transmission hotspots [35], [36]. The concept of carrying capacity in models with homogeneous habitat and logistic growth must be modified in light of the heterogeneous structure of aquatic habitats with local density dependence. Crowding could be highly uneven such that a few habitats would have very high larval densities while others would be scarcely populated. In general, adult mosquito population densities will differ from capacity unless there is fine-tuning in the relationship between the carrying capacities of aquatic habitats and the egg-laying patterns of adult mosquitoes. The rules governing mosquito population densities are more aptly described as a system in which crowding thins the aquatic mosquito populations.
When LSM is integrated into these models, egg-laying behaviour is identified, once again, as an important issue. The effects of LSM are approximately quadratic when mosquitoes continue to lay eggs in treated habitats, and these treated habitats function as egg sinks [23], [24]. This effect is quadratic because LSM has two distinct linear effects that could be created separately by first removing productive habitat, and second replacing that habitat with oviposition traps that absorb just as many eggs. Using larvicides that repel mosquitoes has only the first effect, while using larvicides that do not repel mosquitoes has both effects. The product of these two linear effects creates a non-linear response (i.e. quadratic) much like the one that Macdonald identified in his oft-used analysis [2], [4], [5]. This analysis also raises an operationally relevant question about the repellent effects of modern larvicides at concentrations ordinarily used for field application [37]–[46]. Targeting of these systems [47] can lead to disproportionate efficiencies in the effectiveness of LSM, though practical advice about how to identify productive larval habitats for targeting remains a critical need. In places where the aquatic habitats are in the same places year after year, is possible for control programs to learn and adapt to local systems [47]. Factors that may seem to present technical limitations for LSM – such as the need to target the most productive habitats – can be turned into a long-term operational advantage. It may be possible, for example, to accumulate knowledge about the local mosquito ecology and thereby improve the effectiveness of LSM over time.
The collective results of eleven decades of vector control have been mixed [1], [48]. Overall, our study, along with many others [49], emphasizes the important role of various kinds of heterogeneity in transmission dynamics and control. Heterogeneities can have a strong influence on the ability to measure transmission or predict the outcomes of control programs. A general point to be made is that outcomes probably depend on R0, but they also depend on specific aspects of human and vector behaviours in specific contexts. Interventions that have not been explicitly considered in the Ross-Macdonald model cannot be derived intuitively from the formula for R0. Instead, they must be explicitly modelled and integrated into the underlying theory. Though some simple points can be made about the likely effects of LSM, simple mathematical models can often be misleading unless they identify the appropriate sources of heterogeneity. Application of the theory to LSM in this and other modelling studies [7], [8], [12] is increasingly based on information about local mosquito ecology and its relation to transmission. Given these concerns, even though analysis of mathematical models can help to inform policy, but empirical evidence should perhaps play a stronger role in evaluation of policy and in making policy recommendations. The consideration of LSM, LLINs, IRS, spatial repellents, attractive sugar toxic baits, genetic strategies, oviposition traps, and other new vector control tools designed to reduce transmission of a pathogen by mosquitoes all lead to the realization that there will be some advantages and some disadvantages for each approach, and that intervention success will vary by the transmission context and the efficiency with which programs are implemented. In the context of increasingly widespread insecticide and drugs resistance, limitations in delivery and coverage, and national and international funding constraints, what is urgently needed is programmatic flexibility. Success in ever changing environments will depend on the capacity to select from a suite of options a package of interventions that is best suited for local to national and regional vector-borne diseases prevention goals [50], [51].
Author Contributions
Conceived and designed the experiments: DLS TAP LST TWS SWL. Performed the experiments: DLS. Analyzed the data: DLS TAP. Wrote the paper: DLS TAP LST TWS SWL. Developed the model, created the graphics wrote the first draft: DLS. Analyzed the model: DLS and TAP.
References
- 1.
Ross R (1907) The prevention of malaria in British possessions, Egypt, and parts of America Lancet: 879–887.
- 2. Macdonald G (1952) The analysis of the sporozoite rate. Trop Dis Bull 49: 569–586.
- 3. Macdonald G (1956) Theory of the eradication of malaria. Bull World Health Org 15: 369–387.
- 4. Macdonald G (1956) Epidemiological basis of malaria control. Bull World Health Org 15: 613–626.
- 5.
W.H.O. (2012) The role of larviciding for malaria control in sub-Saharan Africa. Interim Position Statement. Geneva: World Health Organization.
- 6. Worrall E, Fillinger U (2011) Large-scale use of mosquito larval source management for malaria control in Africa: a cost analysis. Malar J 10: 338.
- 7. Dye C (1984) Models for the Population Dynamics of the Yellow Fever Mosquito, Aedes aegypti. J Animal Ecology 53: 247–268.
- 8. Yakob L, Yan G (2010) A network population model of the dynamics and control of African malaria vectors. Trans R Soc Trop Med Hyg 104: 669–675.
- 9. Magori K, Legros M, Puente ME, Focks DA, Scott TW, et al. (2009) Skeeter Buster: a stochastic, spatially explicit modeling tool for studying Aedes aegypti population replacement and population suppression strategies. PLoS Negl Trop Dis 3: e508.
- 10. Focks DA, Daniels E, Haile DG, Keesling JE (1995) A simulation model of the epidemiology of urban dengue fever: literature analysis, model development, preliminary validation, and samples of simulation results. Am J Trop Med Hyg 53: 489–506.
- 11. Focks DA, Haile DG, Daniels E, Mount GA (1993) Dynamic life table model for Aedes aegypti (Diptera: Culicidae): analysis of the literature and model development. J Med Entomol 30: 1003–1017.
- 12. White MT, Griffin JT, Churcher TS, Ferguson NM, Basanez MG, et al. (2011) Modelling the impact of vector control interventions on Anopheles gambiae population dynamics. Parasit Vectors 4: 153.
- 13. Hancock PA, Godfray HCJ (2007) Application of the lumped age-class technique to studying the dynamics of malaria-mosquito-human interactions. Malaria Journal 6: 98.
- 14. Depinay JM, Mbogo CM, Killeen G, Knols B, Beier J, et al. (2004) A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malar J 3: 29.
- 15. Bomblies A, Duchemin J, Eltahir E (2008) Hydrology of malaria: Model development and application to a Sahelian village. Water Resources Research 44: W12445.
- 16. Bomblies A, Duchemin J-B, Eltahir EAB (2009) A mechanistic approach for accurate simulation of village scale malaria transmission. Malaria Journal 8: 223.
- 17. Focks DA, Haile DG, Daniels E, Mount GA (1993) Dynamic life table model for Aedes aegypti (diptera: Culicidae): simulation results and validation. J Med Entomol 30: 1018–1028.
- 18. Le Menach A, McKenzie FE, Flahault A, Smith DL (2005) The unexpected importance of mosquito oviposition behaviour for malaria: Non-productive larval habitats can be sources for malaria transmission. Malaria Journal 4: 23.
- 19. Gu WD, Regens JL, Beier JC, Novak RJ (2006) Source reduction of mosquito larval habitats has unexpected consequences on malaria transmission. Proceedings of the National Academy of Sciences of the United States of America 103: 17560–17563.
- 20. Gu W, Novak RJ (2009) Agent-based modelling of mosquito foraging behaviour for malaria control. Trans R Soc Trop Med Hyg 103: 1105–1112.
- 21. Reiner RC Jr, Perkins TA, Barker CM, Niu T, Chaves LF, et al. (2013) A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970–2010. J R Soc Interface 10: 20120921.
- 22. Arrivillaga J, Barrera R (2004) Food as a limiting factor for Aedes aegypti in water-storage containers. J Vector Ecol 29: 11–20.
- 23. Wong J, Morrison AC, Stoddard ST, Astete H, Chu YY, et al. (2012) Linking oviposition site choice to offspring fitness in Aedes aegypti: consequences for targeted larval control of dengue vectors. PLoS Negl Trop Dis 6: e1632.
- 24. Wong J, Stoddard ST, Astete H, Morrison AC, Scott TW (2011) Oviposition site selection by the dengue vector Aedes aegypti and its implications for dengue control. PLoS Negl Trop Dis 5: e1015.
- 25. Yapabandara AM, Curtis CF, Wickramasinghe MB, Fernando WP (2001) Control of malaria vectors with the insect growth regulator pyriproxyfen in a gem-mining area in Sri Lanka. Acta Trop 80: 265–276.
- 26. Yapabandara AM, Curtis CF (2004) Control of vectors and incidence of malaria in an irrigated settlement scheme in Sri Lanka by using the insect growth regulator pyriproxyfen. J Am Mosq Control Assoc 20: 395–400.
- 27. Sharma SK, Tyagi PK, Upadhyay AK, Haque MA, Adak T, et al. (2008) Building small dams can decrease malaria: a comparative study from Sundargarh District, Orissa, India. Acta Trop 107: 174–178.
- 28. Geissbuhler Y, Kannady K, Chaki PP, Emidi B, Govella NJ, et al. (2009) Microbial larvicide application by a large-scale, community-based program reduces malaria infection prevalence in urban Dar es Salaam, Tanzania. PLoS One 4: e5107.
- 29. Fillinger U, Ndenga B, Githeko A, Lindsay SW (2009) Integrated malaria vector control with microbial larvicides and insecticide-treated nets in western Kenya: a controlled trial. Bull World Health Organ 87: 655–665.
- 30. Fillinger U, Kannady K, William G, Vanek MJ, Dongus S, et al. (2008) A tool box for operational mosquito larval control: preliminary results and early lessons from the Urban Malaria Control Programme in Dar es Salaam, Tanzania. Malar J 7: 20.
- 31. Castro MC, Tsuruta A, Kanamori S, Kannady K, Mkude S (2009) Community-based environmental management for malaria control: evidence from a small-scale intervention in Dar es Salaam, Tanzania. Malar J 8: 57.
- 32. Fillinger U, Lindsay SW (2011) Larval source management for malaria control in Africa: myths and reality. Malar J 10: 353.
- 33.
Lengeler C (2004) Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev: CD000363.
- 34.
Pluess B, Tanser FC, Lengeler C, Sharp BL (2010) Indoor residual spraying for preventing malaria. Cochrane Database Syst Rev: CD006657.
- 35. Bousema T, Griffin JT, Sauerwein RW, Smith DL, Churcher TS, et al. (2012) Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Med 9: e1001165.
- 36.
Chaves LF, Hamer GL, Walker ED, Brown WM, Ruiz MO, et al. (2011) Climatic variability and landscape heterogeneity impact urban mosquito diversity and vector abundance and infection. Ecosphere 2.
- 37. Xue RD, Ali A, Crainich V, Barnard D (2007) Oviposition deterrence and larvicidal activity of three formulations of piperidine repellent (AI3–37220) against field populations of Stegomyia albopicta. J Am Mosq Control Assoc 23: 283–287.
- 38. Xue RD, Barnard DR, Ali A (2003) Laboratory evaluation of 18 repellent compounds as oviposition deterrents of Aedes albopictus and as larvicides of Aedes aegypti, Anopheles quadrimaculatus, and Culex quinquefasciatus. J Am Mosq Control Assoc 19: 397–403.
- 39. Ansari MA, Mittal PK, Razdan RK, Sreehari U (2005) Larvicidal and mosquito repellent activities of Pine (Pinus longifolia, Family: Pinaceae) oil. Journal of Vector Borne Diseases 42: 95–99.
- 40. Karunamoorthi K, Ramanujam S, Rathinasamy R (2008) Evaluation of leaf extracts of Vitex negundo L. (Family: Verbenaceae) against larvae of Culex tritaeniorhynchus and repellent activity on adult vector mosquitoes. Parasitology Research 103: 545–550.
- 41. Pandey SK, Upadhyay S, Tripathi AK (2009) Insecticidal and repellent activities of thymol from the essential oil of Trachyspermum ammi (Linn) Sprague seeds against Anopheles stephensi. Parasitology Research 105: 507–512.
- 42. Rajkumar S, Jebanesan A (2009) Larvicidal and oviposition activity of Cassia obtusifolia Linn (Family: Leguminosae) leaf extract against malarial vector, Anopheles stephensi Liston (Diptera: Culicidae). Parasitol Res 104: 337–340.
- 43. Carroll JF, Tabanca N, Kramer M, Elejalde NM, Wedge DE, et al. (2011) Essential oils of Cupressus funebris, Juniperus communis, and J. chinensis (Cupressaceae) as repellents against ticks (Acari: Ixodidae) and mosquitoes (Diptera: Culicidae) and as toxicants against mosquitoes. Journal of Vector Ecology 36: 258–268.
- 44. Govindarajan M (2011) Larvicidal and repellent properties of some essential oils against Culex tritaeniorhynchus Giles and Anopheles subpictus Grassi (Diptera: Culicidae). Asian Pacific Journal of Tropical Medicine 4: 106–111.
- 45. Govindarajan M, Mathivanan T, Elumalai K, Krishnappa K, Anandan A (2011) Mosquito larvicidal, ovicidal, and repellent properties of botanical extracts against Anopheles stephensi, Aedes aegypti, and Culex quinquefasciatus (Diptera: Culicidae). Parasitology Research 109: 353–367.
- 46. Kumar SV, Mani P, Bastin JTMM, Kumar RA, Ravikumar G (2011) Larvicidal, oviposition deterrent and repellent activity of annona squamosa extracts against hazardous mosquito vectors. International Journal of Pharmacy and Technology 3: 3143–3155.
- 47. Fillinger U, Sombroek H, Majambere S, van Loon E, Takken W, et al. (2009) Identifying the most productive breeding sites for malaria mosquitoes in The Gambia. Malar J 8: 62.
- 48. Najera JA, Gonzalez-Silva M, Alonso PL (2011) Some lessons for the future from the Global Malaria Eradication Programme (1955–1969). PLoS Med 8: e1000412.
- 49. Woolhouse ME, Dye C, Etard JF, Smith T, Charlwood JD, et al. (1997) Heterogeneities in the transmission of infectious agents: implications for the design of control programs. Proc Natl Acad Sci U S A 94: 338–342.
- 50. Eisen L, Beaty BJ, Morrison AC, Scott TW (2009) Proactive vector control strategies and improved monitoring and evaluation practices for dengue prevention. J Med Entomol 46: 1245–1255.
- 51.
Scott TW, Morrison AC (2008) Longitudinal field studies will guide a paradigm shift in dengue prevention. In: Vector-borne Diseases. Understanding the Environmental, Human Health, and Ecological Connections. The National Academies Press: Washington, DC. 132–149.